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

MATLAB程序------材料力学简支梁剪力、弯矩及强度校核、截面设计 绝对原创！！！ MATLAB实现简支梁剪力弯矩的计算、绘图、强度校核、截面设计。 想怎么搞怎么搞！ 我的E-mail:??? fnhpeia2006@126.com 不知道你能否看得懂没有哦。 你自己能看得懂就我服你。 如果本文对你有帮助，我感到非常高兴，证明写得还不错我的这个程序。 OVER。 如需转载，请链接来源于本文。 为了各位更好的理解我的程序， 我已经在相册中上传了图片加以说明，详见《MATLAB材料力学----辅助图片》相册。 我运算的例子用的是：刘鸿文主编的第四版的《材料力学》上册第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\n\n);N=load(F:\MATLAB\work\aaa1.txt);?????????? %读取****save N;??????????????????????????????????????YP=load(F:\MATLAB\work\aaa2.txt);????????? %读取****save YP;OQ=load(F:\MATLAB\work\aaa3.txt);????????? %读取****save OQ;SU=load(F:\MATLAB\work\aaa4.txt);????????? %读取****save SU;JS=load(F:\MATLAB\work\aaa5.txt);save JS; 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的矩阵用以存储所有载荷作用时各分段点的剪力 % 计算所有集中力在各分段点的剪力及弯矩 if N(1)==0?? F1=F1;?? M1=M1;??? elseif N(1)==1??????????????? %判断所受集中力为一个?????? F1(1)=YP(2)*(N(4)-YP(1))/N(4);????? %第一个分段点的剪力?????????????????????????? for n=1:1000??????????? if N(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));?? %集中力作用位置前面分段点的弯矩??????????? else?????????????? F1(n+1)=-YP(2)*YP(1)/N(4);????????? %集中力作用位置后面分段点的剪力?????????????? M1(n+1)=YP(2)*YP(1)*(1-N(4)*n/(1000*N(4)));??