MATLAB程序------材料力学简支梁剪力、弯矩和强度校核、截

《MATLAB程序------材料力学简支梁剪力、弯矩和强度校核、截》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、MATLAB程序------材料力学简支梁剪力、弯矩及强度校核、截面设计绝对原创！！！MATLAB实现简支梁剪力弯矩的计算、绘图、强度校核、截面设计。想怎么搞怎么搞！我的E-mail:   fnhpeia2006@126.com不知道你能否看得懂没有哦。你自己能看得懂就我服你。如果本文对你有帮助，我感到非常高兴，证明写得还不错我的这个程序。OVER。如需转载，请链接来源于本文。为了各位更好的理解我的程序，我已经在相册中上传了图片加以说明，详见《MATLAB材料力学----辅助图片》相册。我运算的例子用的是：刘鸿文主编的第四版的《材料

2、力学》上册第145页的例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=load('F:MATLAB

3、workaaa1.txt');          %读取****saveN;                                      YP=load('F:MATLABworkaaa2.txt');         %读取****saveYP;OQ=load('F:MATLABworkaaa3.txt');         %读取****saveOQ;SU=load('F:MATLABworkaaa4.txt');         %读取****saveSU;JS=load('F:MATLA

4、Bworkaaa5.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的矩阵用以存储所

5、有载荷作用时各分段点的剪力M1=zeros(1,1001);                %建立一个1×1001的矩阵用以存储集中力作用时各分段点的剪力   M2=zeros(1,1001);                %建立一个1×1001的矩阵用以存储均布载荷作用时各分段点的剪力M3=zeros(1,1001);                %建立一个1×1001的矩阵用以存储力偶作用时各分段点的剪力M=zeros(1,1001);                 %建立一个1×1001的矩阵用以存储所有载荷作用时各分段

6、点的剪力%计算所有集中力在各分段点的剪力及弯矩ifN(1)==0  F1=F1;  M1=M1;   elseifN(1)==1               %判断所受集中力为一个      F1(1)=YP(2)*(N(4)-YP(1))/N(4);     %第一个分段点的剪力                          forn=1:1000           ifN(4)*n/1000<=YP(1)         %判断分段点位于集中力作用位置前面              F1(n+1)=YP(2)*(N(4)-

7、YP(1))/N(4);  %集中力作用位置前面分段点的剪力                                M1(n+1)=YP(2)*(N(4)-YP(1))*N(4)*n/(1000*N(4));  %集中力作用位置前面分段点的弯矩           else              F1(n+1)=-YP(2)*YP(1)/N(4);         %集中力作用位置后面分段点的剪力              M1(n+1)=YP(2)*YP(1)*(1-N(4)*n/(1000*N(4)));    %集

8、中力作用位置后面分段点的弯矩            end        end           else         fori=1:N(1)           F1(1)=F1(1)+YP(i,2)*(N(4)-YP(i,