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

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/500-b1)^2*q/2; else

Fsx2(ff)=0; Mx2(ff)=0; end end end

for pp=1:leng if pp<=xx3; Fsx3(pp)=0; Mx3(pp)=M; else

Fsx3(pp)=0; Mx3(pp)=0; end end flag1=1; flag2=1;

Fsx=Fsx3+Fsx2+Fsx1; Mx=Mx3+Mx2+Mx1; for jj=1:leng

if abs(Fsx(jj))>=abs(Fsx(flag1)) flag1=jj; end

if abs(Mx(jj))>=abs(Mx(flag2)) flag2=jj; end end yy=0*ones(1,leng); yyy2=abs(Fsx); yyyy2=max(yyy2); yyy3=abs(Mx); yyyy3=max(yyy3); axes(handles.axes2); grid on

plot(x,yy,'LineWidth',2); axis([0,L,-2*yyyy2,2*yyyy2]); hold on

plot(x,Fsx,'LineWidth',2); hold on

xxx1=L/500*(flag1-1);

plot(xxx1,Fsx(flag1),'bo','MarkerSize',10); hold on

for kk=1:4:leng

plot([x(kk),x(kk)],[0,Fsx(kk)],'LineWidth',0.4) end

xlabel('x/m'); ylabel('??á|Fs/N'); axes(handles.axes3); grid on

plot(x,yy,'LineWidth',2);

axis([0,L,-1.2*yyyy3,1.2*yyyy3]); hold on

plot(x,Mx,'LineWidth',2); hold on

xxx2=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) end

str1=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); set(handles.edit13,'String',str4);

从程序文本框中输入参数值，输入完毕后点击“计算”按钮，在作图框内便显示出梁的剪力和弯矩图和最大值点，在结果文本框内便得出相应最大值和最大值所在位置。

其他梁的程序设计过程也如此，这里不再叙述

四、程序算例

下面给出各个梁的算例。 1、左固定端悬臂梁计算例题

例1.有一左固定端悬臂梁长l=4m，在梁中间作用有集中力P=4N，梁右端作用有集中力偶m=5N·m，左端作用均布载荷，载荷集度q=2N/m,作用图示如下。请作出梁的剪力弯矩图，并求出剪力弯矩的最大值与其出现位置。

计算结果：

图5 计算结果截图

根据题意可得FA=8N，MB=7N·m，将梁分为AC、BC段进行分析，即可作出梁的剪力、弯矩图，所作图形与程序所作结果相同。从剪力、弯矩图上显然可以看出剪力最大值Fmax=8N，弯矩最大值Mz=7N·m，均与程序相同，故此段程序编写正确。

2、右固定端悬臂梁计算例题

例2.有一右固定端悬臂梁长l=2m，在梁中间作用有集中力F=5N，梁右端作用竖直向下的集中力P，大小为5N，在梁的0.5m到1.5m作用有均布载荷，载荷集度为q=2N/m，梁左端作用有顺时针集中力偶m=10N·m，作用图示如下。请作出梁的剪力弯矩图，并求出剪力弯矩的最大值与其出现位置。

计算结果：

根据题意可得FB=7N，MA=5N·m，将梁分为AD、DC、CE、EB段进行分析即可作出梁的剪力弯矩图，所作图形与程序结果相同。从剪力、弯矩图上显然可以

看出剪力最大值Fmax=7N（出现在EB段，本段程序显示结果在E截面处），弯矩最大值Mz=5N·m（出现在AD段，本段程序显示结果在A截面处），均与程序相同，故此段程序编写正确。

图6 计算结果截图

3、简支梁计算例题

例3. 已知一横梁长l=4m，两端铰支。现在x=2m处作用一大小为4N、方向竖直向下的集中力，并且在整根梁上作用有竖直向下的均布载荷，大小为q=2N/m，作用示意图如下。试做出梁上的剪力和弯矩图，并求出剪力弯矩的最大、最小值及它们出现的位置。

计算结果：

求支反力得FA=6N，FB=6N，作梁的剪力弯矩图与程序所作结果相同，从剪力弯矩图上显然可以得到剪力最大值Fmax=6N，弯矩最大值Mz=8N·m，与程序相同。理论计算绘制图示上可以看出剪力最大值在梁的左右两端点处取到（x=0和x=4m处），程序中只显示了一个（x=4m处）。

图7 计算结果截图

4、左外伸梁计算例题

例4. 已知一左外伸梁长3m，可动铰支座距左端面1m，在x=2m处作用一大小为P=10N，方向竖直向下的集中力，图形如下。试做出梁上的剪力和弯矩图，并求出剪力弯矩的最大、最小值及它们出现的位置。

计算结果：

图8 计算结果截图

求支反力得FA=2.5N FB=7.5N，作梁的剪力弯矩图与程序所作结果基本相

搜索“diyifanwen.net”或“第一范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，第一范文网，提供最新综合文库哈工大材料力学上机报告--绘制梁的剪力弯矩图(2)全文阅读和word下载服务。