哈工大材料力学上机报告--绘制梁的剪力弯矩图

1、Harbin Institute of Technology材料力学上机作业课程名称： 材料力学 设计题目： 绘制梁的剪力弯矩图 院 系： 机电工程学院 班 级： 分 析 者： 学 号： 指导教师： 张桂莲 设计时间： 2014.6.3 哈尔滨工业大学一、设计题目题目4 绘制梁的剪力弯矩图输入：1.梁的总长度l2.支撑条件及量的各区段长度输入（左、右固定端悬臂梁；简支梁；左、右、双外伸梁）3.各载荷大小、作用位置及方向（qi、ai、bi；pj、cj、mk、dk）输出：1. 结构构型图（图示）2. 剪力、弯矩（图示）3. 输出剪力、弯矩的最大值及截面位置。二、程序介绍1、简介本程序编译平台为MA

2、TLAB，能够对六种梁（左固定端悬臂梁、右固定端悬臂梁、简支梁、左外伸梁、右外伸梁、双外伸梁）在集中载荷、均布载荷及弯矩作用下的剪力、弯矩进行求解并绘制剪力、弯矩图。本程序基于MATLAB的图形用户界面（GUI）功能，操作简单、运算稳定、界面简约美观，有良好的交互性。2、程序使用方法 此程序运行环境为安装有MATLAB R2011a或以上版本的计算机。要运行此程序，首先打开MATLAB，在其中打开“KaiShiJieMian.m”文件，在函数编辑窗口里点击运行按钮（或按F5键），程序便开始启动。启动画面如下：图1 启动界面 点击“进入”即可进入选择梁的类型界面，在此界面中，使用者根据实际情况选

3、择要研究的梁的类型，点击即可进入相应梁的计算界面。例如想要计算双外伸梁，则点击“6、双外伸梁”即可进入计算界面，若想返回开始界面，则点击“返回开始界面”即可。图2 梁的选择界面图3 梁的计算界面 此程序可以计算梁在集中力、均布力和集中力偶的作用下梁的剪力和弯矩，并求出剪力的最大值和弯矩的最大值。以双外伸梁为例，对照左侧参数输入提示，在相应文本框里输入梁的基本参数，如弹性模量E、长度L，和载荷大小及位置，如均布力q、均布力起始位置b1、结束位置b2。输入完毕后点击“计算”按钮，在作图框内便显示出梁的剪力和弯矩图和最大值点，在结果文本框内便得出相应最大值和最大值所在位置。若要改变参数，先点击“清除

4、”按钮，使做图框和结果文本框初始化，再更改相应参数，再次点击“计算”。三、程序设计过程 以左固定端悬臂梁为例，首先选取坐标系，确定载荷大小、类型和位置坐标，然后分别计算单个载荷所对应的支反力和力矩，用公式计算位置x处对应的剪力和弯矩值，最后把各个载荷产生的在位置x处的剪力和弯矩值相加，便得出总的剪力和弯矩值。要注意的是，集中力和均布力取向下为正，外力偶取逆时针为正。所用公式如下：图4 左固定端悬臂梁计算示意图 设F产生的支反力为FOy1,力偶为MO1，对应的剪力为Fsx1,弯矩为Mx1；q产生的支反力为FOy2,力偶为MO2，对应的剪力为Fsx2,弯矩为Mx2；M产生的支反力为FOy3,力偶为

5、MO3，对应的剪力为Fsx3,弯矩为Mx3。则：FOy1=F;MO1=FOy1*a;FOy2=q*(b2-b1);MO2=q*(b2-b1)*(b2+b1)/2;FOy3=0;MO3=M;对于F： 当x=a时， Fsx1=-FOy1; Mx1=FOy1*x-MO1; 当ax=L时 Fsx1=0; Mx1=0;对于q： 当x=b1时 Fsx2=-FOy2; Mx2=FOy2*x-MO2; 当b1x=b2时 Fsx2=-FOy2 + q*(x-b1); Mx2(ff)=FOy2*x-MO2-(x-b1)2*q/2; 当b2x=L时 Fsx2=0; Mx2=0;对于M： 当x=c时 Fsx3=0;

6、Mx3=M; 当cx=L时 Fsx3=0; Mx3=0;最后，设合剪力为Fsx,合弯矩为Mx,则：Fsx=Fsx1+Fsx2+Fsx3;Mx=Mx1+Mx2+Mx3;对应程序如下所示：E=109*str2double(get(handles.edit1,String);L=str2double(get(handles.edit2,String);F=str2double(get(handles.edit3,String);a=str2double(get(handles.edit4,String);q=str2double(get(handles.edit5,String);b1=str2do

7、uble(get(handles.edit6,String);b2=str2double(get(handles.edit7,String);M=str2double(get(handles.edit8,String);c=str2double(get(handles.edit9,String);x=0:L/500:L;FOy1=F;MO1=FOy1*a;FOy2=q*(b2-b1);MO2=q*(b2-b1)*(b2+b1)/2;FOy3=0;MO3=M;leng=length(x);xx1=round(a/L*leng);xx21=round(b1/L*leng);xx22=round(b

8、2/L*leng);xx3=round(c/L*leng);for ii=1:leng if ii=xx1 Fsx1(ii)=-FOy1; Mx1(ii)=FOy1*(ii-1)*L/500-MO1; else Fsx1(ii)=0; Mx1(ii)=0; endend for ff=1:leng if ff=xx21 Fsx2(ff)=-FOy2; Mx2(ff)=FOy2*(ff-1)*L/500-MO2; else if ff=xx22 Fsx2(ff)=-FOy2 + q*(ff-1)*L/500-b1); Mx2(ff)=FOy2*(ff-1)*L/500-MO2-(ff-1)*L/

9、500-b1)2*q/2; else Fsx2(ff)=0; Mx2(ff)=0; end endendfor pp=1:leng if pp=abs(Fsx(flag1) flag1=jj; end if abs(Mx(jj)=abs(Mx(flag2) flag2=jj; endend yy=0*ones(1,leng);yyy2=abs(Fsx);yyyy2=max(yyy2);yyy3=abs(Mx);yyyy3=max(yyy3);axes(handles.axes2);grid onplot(x,yy,LineWidth,2);axis(0,L,-2*yyyy2,2*yyyy2);

10、hold onplot(x,Fsx,LineWidth,2);hold onxxx1=L/500*(flag1-1);plot(xxx1,Fsx(flag1),bo,MarkerSize,10);hold onfor kk=1:4:leng plot(x(kk),x(kk),0,Fsx(kk),LineWidth,0.4)endxlabel(x/m);ylabel(Fs/N);axes(handles.axes3);grid onplot(x,yy,LineWidth,2);axis(0,L,-1.2*yyyy3,1.2*yyyy3);hold onplot(x,Mx,LineWidth,2)

11、;hold onxxx2=L/500*(flag2-1);plot(xxx2,Mx(flag2),bo,MarkerSize,10);for ll=1:4:leng plot(x(ll),x(ll),0,Mx(ll),LineWidth,0.4)endstr1=num2str(yyyy2);str2=num2str(xxx1);set(handles.edit10,String,str1);set(handles.edit11,String,str2);str3=num2str(yyyy3);str4=num2str(xxx2);set(handles.edit12,String,str3);

12、set(handles.edit13,String,str4);从程序文本框中输入参数值，输入完毕后点击“计算”按钮，在作图框内便显示出梁的剪力和弯矩图和最大值点，在结果文本框内便得出相应最大值和最大值所在位置。 其他梁的程序设计过程也如此，这里不再叙述四、程序算例下面给出各个梁的算例。1、左固定端悬臂梁计算例题例1.有一左固定端悬臂梁长l=4m，在梁中间作用有集中力P=4N，梁右端作用有集中力偶m=5Nm，左端作用均布载荷，载荷集度q=2N/m,作用图示如下。请作出梁的剪力弯矩图，并求出剪力弯矩的最大值与其出现位置。计算结果：图5 计算结果截图 根据题意可得FA=8N，MB=7Nm，将梁分为

13、AC、BC段进行分析，即可作出梁的剪力、弯矩图，所作图形与程序所作结果相同。从剪力、弯矩图上显然可以看出剪力最大值Fmax=8N，弯矩最大值Mz=7Nm，均与程序相同，故此段程序编写正确。2、右固定端悬臂梁计算例题例2.有一右固定端悬臂梁长l=2m，在梁中间作用有集中力F=5N，梁右端作用竖直向下的集中力P，大小为5N，在梁的0.5m到1.5m作用有均布载荷，载荷集度为q=2N/m，梁左端作用有顺时针集中力偶m=10Nm，作用图示如下。请作出梁的剪力弯矩图，并求出剪力弯矩的最大值与其出现位置。计算结果： 根据题意可得FB=7N，MA=5Nm，将梁分为AD、DC、CE、EB段进行分析即可作出梁的

14、剪力弯矩图，所作图形与程序结果相同。从剪力、弯矩图上显然可以看出剪力最大值Fmax=7N（出现在EB段，本段程序显示结果在E截面处），弯矩最大值Mz=5Nm（出现在AD段，本段程序显示结果在A截面处），均与程序相同，故此段程序编写正确。图6 计算结果截图3、简支梁计算例题例3. 已知一横梁长l=4m，两端铰支。现在x=2m处作用一大小为4N、方向竖直向下的集中力，并且在整根梁上作用有竖直向下的均布载荷，大小为q=2N/m，作用示意图如下。试做出梁上的剪力和弯矩图，并求出剪力弯矩的最大、最小值及它们出现的位置。计算结果： 求支反力得FA=6N，FB=6N，作梁的剪力弯矩图与程序所作结果相同，从剪

15、力弯矩图上显然可以得到剪力最大值Fmax=6N，弯矩最大值Mz=8Nm，与程序相同。理论计算绘制图示上可以看出剪力最大值在梁的左右两端点处取到（x=0和x=4m处），程序中只显示了一个（x=4m处）。图7 计算结果截图4、左外伸梁计算例题例4. 已知一左外伸梁长3m，可动铰支座距左端面1m，在x=2m处作用一大小为P=10N，方向竖直向下的集中力，图形如下。试做出梁上的剪力和弯矩图，并求出剪力弯矩的最大、最小值及它们出现的位置。计算结果：图8 计算结果截图 求支反力得FA=2.5N FB=7.5N，作梁的剪力弯矩图与程序所作结果基本相同（程序所作剪力弯矩图中有断点存在，这是由于本程序绘制剪力弯

16、矩图采用的是散点法，当剪力弯矩值有突变时程序计算的因数值不连续而产生“跳跃”现象），从剪力弯矩图上显然可以得到剪力最大值Fmax=7.5N，弯矩最大值Mz=7.5Nm，与程序相同。5、右外伸梁计算例题例5. 已知一右外伸梁长1.5米，右支点距右端面0.5米，在两支架中间作用一大小为P=5N，方向竖直向上的集中力，图形如下，试做出梁上的剪力和弯矩图，并求出剪力弯矩的最大、最小值及它们出现的位置。计算结果：图9 计算结果截图 求支反力得FA=2.5N FB=2.5N，作梁的剪力弯矩图与程序所作结果相同（程序所作剪力弯矩图中有断点存在，这是由于本程序绘制剪力弯矩图采用的是散点法，当剪力弯矩值有突变时

17、程序计算的因数值不连续而产生“跳跃”现象），从剪力弯矩图上显然可以得到剪力最大值Fmax=2.5N，弯矩最大值Mz=1.25Nm，与程序计算结果相同。理论计算绘制图示上可以看出剪力和弯矩最大值应该在x=0.5m处取到，而程序计算的最大值分别在x=0.999m和x=0.501m取到，这是由于在程序编写过程中把梁在集中作用出分成了两段，集中力作用以左部分所有的x都小于0.5m，集中力作用以右部分所有x都大于0.5m，因而显示结果会出现与理论计算的微小差别。6、双外伸梁计算例题例6. 已知一双外伸梁长4米，固定铰支座距左端面1m,可动铰支座距右端面1m，在距离梁的左端0.5m处作用一大小为F=4N、

18、方向竖直向下的集中力，两铰支座之间作用有向下的均布载荷，大小为2N/m，图形如下。试做出梁上的剪力和弯矩图，并求出剪力弯矩的最大、最小值及它们出现的位置。计算结果： 求支反力得FA=1N FB=7N，作梁的剪力弯矩图与程序所作结果相同（程序所作剪力图中有断点存在，这是由于本程序绘制剪力弯矩图采用的是散点法，当剪力弯矩值有突变时程序计算的因数值不连续而产生“跳跃”现象），从剪力弯矩图上显然可以得到剪力最大值Fmax=4N，弯矩最大值Mz=-2Nm，与程序计算结果相同。理论计算绘制图示上可以看出剪力和弯矩最大值分别在x=0.5mx=1m处取得，与程序计算结果也相同，因而此段程序运行正确。图10 计

19、算结果截图五、程序缺陷 由于本程序绘制剪力弯矩图采用的是散点法，因而在剪力或弯矩有突变时图示上会产生断点，未连接上。可以采用“line”指令，在数据有突变时直接直线连接。六、体会与收获在这次上机编程作业中，我虽然花费了很长时间，但是收获很多。在此过程中，我对MATLAB语言有了更加熟悉的运用，这将成为我一生的财富，MATLAB也必将成为我日后学习工作的强大工具。在看到自己做出的最终结果以及它运行时的高效、简洁，我不禁有了一丝成就感。我为能有一次这样锻炼的机会而感到幸运，希望以后还会更多这样的机会充实自己。对于材料力学这门课程，我想也会因为这次程序的编写而记忆的更加深刻，我一定要争取多掌握一些这

20、方面的知识，然后在工作中充分的应用进去！当然，作为初学者，制作出的软件不免瑕疵很多，还烦请老师批评指正。但我相信，在今后的学习中，我对MATLAB的应用能力会愈加娴熟。七、程序源代码1.左悬臂梁的源程序E=getappdata(BianXing_ZuoGuDingDuan,E);L=getappdata(BianXing_ZuoGuDingDuan,L);F=getappdata(BianXing_ZuoGuDingDuan,F);a=getappdata(BianXing_ZuoGuDingDuan,a);q=getappdata(BianXing_ZuoGuDingDuan,q);b1=g

21、etappdata(BianXing_ZuoGuDingDuan,b1);b2=getappdata(BianXing_ZuoGuDingDuan,b2);M=getappdata(BianXing_ZuoGuDingDuan,M);c=getappdata(BianXing_ZuoGuDingDuan,c);dyuan=str2double(get(handles.edit1,String);Dyuanhuan=str2double(get(handles.edit2,String);dyuanhuan=str2double(get(handles.edit3,String);hjuxing

22、=str2double(get(handles.edit4,String);bjuxing=str2double(get(handles.edit5,String);Izgongzixing=str2double(get(handles.edit6,String);Iz1=pi*dyuan4/64;Iz2=pi*(Dyuanhuan4-dyuanhuan4)/64;Iz3=bjuxing*hjuxing3/12;Iz4=Izgongzixing;Iz=Iz1+Iz2+Iz3+Iz4;x=0:L/500:L;leng=length(x);xx1=round(a/L*leng);xx21=roun

23、d(b1/L*leng);xx22=round(b2/L*leng);xx3=round(c/L*leng);for ii=1:leng if ii=abs(vx(flag1) flag1=jj; endend yy=0*ones(1,leng);yyy2=abs(vx);yyyy2=max(yyy2);axes(handles.axes5);grid onplot(x,yy,LineWidth,2);axis(0,L,-2*yyyy2,2*yyyy2);hold onplot(x,100000*vx,LineWidth,2);hold onxxx1=L/500*(flag1-1);plot(

24、xxx1,vx(flag1),bo,MarkerSize,10);hold onfor kk=1:4:leng plot(x(kk),x(kk),0,vx(kk),LineWidth,0.4)endxlabel(x/m);ylabel(Fs/N);2.右悬臂梁的源程序E=109*str2double(get(handles.edit1,String);L=str2double(get(handles.edit2,String);F=str2double(get(handles.edit3,String);a=str2double(get(handles.edit4,String);q=str2

25、double(get(handles.edit5,String);b1=str2double(get(handles.edit6,String);b2=str2double(get(handles.edit7,String);M=str2double(get(handles.edit8,String);c=str2double(get(handles.edit9,String);x=0:L/500:L;FOy1=F;MO1=FOy1*a;FOy2=q*(b2-b1);MO2=q*(b2-b1)*(b2+b1)/2;FOy3=0;MO3=M;leng=length(x);xx1=round(a/

26、L*leng);xx21=round(b1/L*leng);xx22=round(b2/L*leng);xx3=round(c/L*leng);for ii=1:leng if ii=xx1 Fsx1(ii)=-FOy1; Mx1(ii)=FOy1*(ii-1)*L/500-MO1; else Fsx1(ii)=0; Mx1(ii)=0; endend for ff=1:leng if ff=xx21 Fsx2(ff)=-FOy2; Mx2(ff)=FOy2*(ff-1)*L/500-MO2; else if ff=xx22 Fsx2(ff)=-FOy2 + q*(ff-1)*L/500-b1

27、); Mx2(ff)=FOy2*(ff-1)*L/500-MO2-(ff-1)*L/500-b1)2*q/2; else Fsx2(ff)=0; Mx2(ff)=0; end endendfor pp=1:leng if pp=abs(Fsx(flag1) flag1=jj; end if abs(Mx(jj)=abs(Mx(flag2) flag2=jj; endend yy=0*ones(1,leng);yyy2=abs(Fsx);yyyy2=max(yyy2);yyy3=abs(Mx);yyyy3=max(yyy3);axes(handles.axes2);grid onplot(x,y

28、y,LineWidth,2);axis(0,L,-2*yyyy2,2*yyyy2);hold onplot(x,Fsx,LineWidth,2);hold onxxx1=L/500*(flag1-1);plot(xxx1,Fsx(flag1),bo,MarkerSize,10);hold onfor kk=1:4:leng plot(x(kk),x(kk),0,Fsx(kk),LineWidth,0.4)endxlabel(x/m);ylabel(Fs/N);axes(handles.axes3);grid onplot(x,yy,LineWidth,2);axis(0,L,-1.2*yyyy

29、3,1.2*yyyy3);hold onplot(x,Mx,LineWidth,2);hold onxxx2=L/500*(flag2-1);plot(xxx2,Mx(flag2),bo,MarkerSize,10);for ll=1:4:leng plot(x(ll),x(ll),0,Mx(ll),LineWidth,0.4)endxlabel(x/m);ylabel(M/N.m);str1=num2str(yyyy2);str2=num2str(xxx1);set(handles.edit10,String,str1);set(handles.edit11,String,str2);str

30、3=num2str(yyyy3);str4=num2str(xxx2);set(handles.edit12,String,str3);set(handles.edit13,String,str4);% - Executes on button press in togglebutton2.function togglebutton2_Callback(hObject, eventdata, handles)% hObject handle to togglebutton2 (see GCBO)% eventdata reserved - to be defined in a future v

31、ersion of MATLAB% handles structure with handles and user data (see GUIDATA)close(gcf)a=gcf;XuanZeLiangDeLeiXing% Hint: get(hObject,Value) returns toggle state of togglebutton2% - Executes on button press in togglebutton3.function togglebutton3_Callback(hObject, eventdata, handles)% hObject handle t

32、o togglebutton3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)a=gcf;BianXing_ZuoGuDingDuanE=109*str2double(get(handles.edit1,String);L=str2double(get(handles.edit2,String);F=str2double(get(handles.edit3,String)

33、;a=str2double(get(handles.edit4,String);q=str2double(get(handles.edit5,String);b1=str2double(get(handles.edit6,String);b2=str2double(get(handles.edit7,String);M=str2double(get(handles.edit8,String);c=str2double(get(handles.edit9,String);setappdata(BianXing_ZuoGuDingDuan,E,E);setappdata(BianXing_ZuoG

34、uDingDuan,L,L);setappdata(BianXing_ZuoGuDingDuan,F,F);setappdata(BianXing_ZuoGuDingDuan,a,a);setappdata(BianXing_ZuoGuDingDuan,q,q);setappdata(BianXing_ZuoGuDingDuan,b1,b1);setappdata(BianXing_ZuoGuDingDuan,b2,b2);setappdata(BianXing_ZuoGuDingDuan,M,M);setappdata(BianXing_ZuoGuDingDuan,c,c); % Hint:

35、 get(hObject,Value) returns toggle state of togglebutton3% - Executes on button press in togglebutton5.function togglebutton5_Callback(hObject, eventdata, handles)% hObject handle to togglebutton5 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with ha

36、ndles and user data (see GUIDATA)axes(handles.axes2);cla;axes(handles.axes3);cla;str1=num2str(0);str2=num2str(0);set(handles.edit10,String,str1);set(handles.edit11,String,str2);str3=num2str(0);str4=num2str(0);set(handles.edit12,String,str3);set(handles.edit13,String,str4);3.简支梁的源程序E=109*str2double(g

37、et(handles.edit1,String);L=str2double(get(handles.edit2,String);F=str2double(get(handles.edit3,String);a=str2double(get(handles.edit4,String);q=str2double(get(handles.edit5,String);b1=str2double(get(handles.edit6,String);b2=str2double(get(handles.edit7,String);M=str2double(get(handles.edit8,String);

38、c=str2double(get(handles.edit9,String);x=0:L/500:L;FAy1=F-F*a/L;FBy1=F*a/L;FBy2=q*(b2-b1)*(b2+b1)/(2*L);FAy2=q*(b2-b1)-FBy2;FAy3=M/L;FBy3=M/L;leng=length(x);xx1=round(a/L*leng);xx21=round(b1/L*leng);xx22=round(b2/L*leng);xx3=round(c/L*leng);for ii=1:leng if ii=xx1 Fsx1(ii)=-FAy1; Mx1(ii)=FAy1*(ii-1)

39、*L/500; else Fsx1(ii)=FBy1; Mx1(ii)=FBy1*(L-(ii-1)*L/500); endend for ff=1:leng if ff=xx21 Fsx2(ff)=-FAy2; Mx2(ff)=FAy2*(ff-1)*L/500; else if ff=xx22 Fsx2(ff)=q*(ff-1)*L/500-b1)-FAy2; Mx2(ff)=FAy2*(ff-1)*L/500-(ff-1)*L/500-b1)2*q/2; else Fsx2(ff)=FBy2; Mx2(ff)=FBy2*(L-(ff-1)*L/500); end endendfor pp

40、=1:leng if pp=abs(Fsx(flag1) flag1=jj; end if abs(Mx(jj)=abs(Mx(flag2) flag2=jj; endend yy=0*ones(1,leng);yyy2=abs(Fsx);yyyy2=max(yyy2);yyy3=abs(Mx);yyyy3=max(yyy3);axes(handles.axes2);grid onplot(x,yy,LineWidth,2);axis(0,L,-1.2*yyyy2,1.2*yyyy2);hold onplot(x,Fsx,LineWidth,2);hold onxxx1=L/500*(flag

41、1-1);plot(xxx1,Fsx(flag1),bo,MarkerSize,10);hold onfor kk=1:4:leng plot(x(kk),x(kk),0,Fsx(kk),LineWidth,0.4)endxlabel(x/m);ylabel(Fs/N);axes(handles.axes4);grid onplot(x,yy,LineWidth,2);axis(0,L,-1.2*yyyy3,1.2*yyyy3);hold onplot(x,Mx,LineWidth,2);hold onxxx2=L/500*(flag2-1);plot(xxx2,Mx(flag2),bo,Ma

42、rkerSize,10);for ll=1:4:leng plot(x(ll),x(ll),0,Mx(ll),LineWidth,0.4)endxlabel(x/m);ylabel(M/N.m);str1=num2str(yyyy2);str2=num2str(xxx1);set(handles.edit10,String,str1);set(handles.edit11,String,str2);str3=num2str(yyyy3);str4=num2str(xxx2);set(handles.edit14,String,str3);set(handles.edit15,String,st

43、r4);4.左外伸梁的源程序E=109*str2double(get(handles.edit1,String);L=str2double(get(handles.edit2,String);F=str2double(get(handles.edit4,String);a=-str2double(get(handles.edit5,String);q=str2double(get(handles.edit6,String);b1=-str2double(get(handles.edit7,String);b2=-str2double(get(handles.edit8,String);M=st

44、r2double(get(handles.edit9,String);c=-str2double(get(handles.edit10,String);d=-str2double(get(handles.edit3,String);x=0:-L/500:-L;leng=length(x);xx1=round(-a/L*leng);xx21=round(-b1/L*leng);xx22=round(-b2/L*leng);xx3=round(-c/L*leng);xx4=round(-d/L*leng);if -a=-d FBy1=F*a/d; FAy1=F-FBy1; for ii=1:leng if ii=xx1 Fsx1(ii)=-FAy1; Mx1(ii)=FAy1*(ii-1)*L/500; else if ii=xx4 Fsx1(ii)=FBy1; Mx1(ii)=FBy1*(-d-(ii-1)*L/500); else Fsx1(ii)=0; Mx1(ii)=0; end end endelse FBy1=F*a/d; FAy1=FBy1-F; for ii=1:leng if ii=xx4 Fsx1(ii)=FAy1; Mx1(ii)=-FAy1*(ii-1)*L