MATLAB程序------材料力学简支梁剪力、弯矩及强度校核、截面设计绝对原创！！！MATLAB实现简支梁剪力弯矩的计算、绘图、强度校核、截面设计。我运算的例子用的是：刘鸿文主编的第四版的《材料力学》上册的例5.2为例。计算机所得结果与例子一致，程序编制成功！程序如下：%程序最后完成于2009.06.07韦**编%计算简支梁在多个集中力，均布载荷，力偶同时作用时的剪力，弯矩并绘制其图形%N存放aaa1.txt的数据%YP存放aaa2.txt的数据%OQ存放aaa3.txt的数据%SU存放aaa4.txt的数据%JS存放aaa5.txt的数据%fff存放各分段点的剪力%mmm存放各分段点的弯矩%Fsmax存放最大剪力%Mmax存放最大弯矩clc;fprintf('\n\n\n');N=load('F:\MATLAB\work\aaa1.txt');%读取****saveN;YP=load('F:\MATLAB\work\aaa2.txt');%读取****saveYP;OQ=load('F:\MATLAB\work\aaa3.txt');%读取****saveOQ;SU=load('F:\MATLAB\work\aaa4.txt');%读取****saveSU;JS=load('F:\MATLAB\work\aaa5.txt');saveJS;F1=zeros(1,1001);%建立一个1×1001的矩阵用以存储集中力作用时各分段点的剪力F2=zeros(1,1001);%建立一个1×1001的矩阵用以存储均布载荷作用时各分段点的剪力F3=zeros(1,1001);%建立一个1×1001的矩阵用以存储力偶作用时各分段点的剪力F=zeros(1,1001);%建立一个1×1001的矩阵用以存储所有载荷作用时各分段点的剪力M1=zeros(1,1001);%建立一个1×1001的矩阵用以存储集中力作用时各分段点的剪力M2=zeros(1,1001);%建立一个1×1001的矩阵用以存储均布载荷作用时各分段点的剪力M3=zeros(1,1001);%建立一个1×1001的矩阵用以存储力偶作用时各分段点的剪力M=zeros(1,1001);%建立一个1×1001的矩阵用以存储所有载荷作用时各分段点的剪力%计算所有集中力在各分段点的剪力及弯矩ifN(1)==0F1=F1;M1=M1;elseifN(1)==1%判断所受集中力为一个F1(1)=YP(2)*(N(4)-YP(1))/N(4);%第一个分段点的剪力forn=1:1000ifN(4)*n/1000<=YP(1)%判断分段点位于集中力作用位置前面F1(n+1)=YP(2)*(N(4)-YP(1))/N(4);%集中力作用位置前面分段点的剪力M1(n+1)=YP(2)*(N(4)-YP(1))*N(4)*n/(1000*N(4));%集中力作用位置前面分段点的弯矩elseF1(n+1)=-YP(2)*YP(1)/N(4);%集中力作用位置后面分段点的剪力M1(n+1)=YP(2)*YP(1)*(1-N(4)*n/(1000*N(4)));%集中力作用位置后面分段点的弯矩endendelsefori=1:N(1)F1(1)=F1(1)+YP(i,2)*(N(4)-YP(i,1))/N(4);%各集中力作用时第一个分段点的剪力累和forn=1:1000ifN(4)*n/1000<=YP(i,1)%判断分段点位于各集中力作用位置前面F1(n+1)=F1(n+1)+YP(i,2)*(N(4)-YP(i,1))/N(4);%各集中力各自作用时在其作用位置前面分段点的剪力累和M1(n+1)=M1(n+1)+YP(i,2)*(N(4)-YP(i,1))...*N(4)*n/(1000*N(4));%各集中力各自作用时在其作用位置前面分段点的弯矩累和elseF1(n+1)=F1(n+1)-YP(i,2)*YP(i,1)/N(4);%各集中力各自作用时在其作用位置后面分段点剪力累和M1(n+1)=M1(n+1)+YP(i,2)*YP(i,1)...*(1-N(4)*n/(1000*N(4)));%各集中力各自作用时在其作用位置后面分段点的弯矩累和endendendend%计算所有均布载荷在各分段点的剪力及弯矩ifN(2)==0F2=F2;M2=M2;elseifN(2)==1%判断所受均布载荷为一个F2(1)=OQ(3)*OQ(2)*(N(4)-OQ(1)-OQ(2)/2)/N(4);%第一个分段点的剪力forn=1:1000ifN(4)*n/1000<=OQ(1)%判断分段点位于均布载荷作用位置前面F2(n+1)=OQ(3)*OQ(2)*(N(4)-OQ(1)-OQ(2