资源描述:
《梁计算机分析程序》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、简支梁计算机分析程functionXQM=QMDJ(x’MPQ)[n,m]=size(x);L=x(m);xl=[];fori=l:m-lxl=[xl,linspace(x(i),x(i+l),50)];end%离散简支梁MM=zeros(size(xl));QQ=zeros(size(xl));%将离散点的剪力和弯矩数组QQ、MM赋零[m,t]=size(MPQ);[t,n]=size(xl);fori=l:m;switchMPQ(i,l)case1a=MPQ(i,3);M=MPQ.(i,2);RA=-M/L;QQ二QQ+RA;MM=MM+RA*xl;ifa>
2、0&a0&b3、MIN(xl,QQ,MM);TU_QM(xbQQ,MM);分析实例clear;clc;X二[024]MPQ=[2,30,2,0]XQM=QMDJ(xzMPQ);20100-10'2000.511.522.533.54弯矩图右端外伸梁计算机程序functionXQM=QMYW(x丄乙MPQ)[n,m]=size(x);L=x(m);xl=[];fori=l:m-lxl=[xl,linspace(x(i),x(i+l)/50)];endMM=zeros(size(xl));QQ=zeros(size(xl));[m,t]=size(MPQ);[t,n]=size(x
4、l);fori=l:mswitchMPQ(i,l);case1a=MPQ(i,3);M=MPQ(i,2);RA=-M/(L-L2);RB=RA;QQ二QQ+RA;MM=MM+RA*xl;MM=QMM(n,xl,a,M,MM);[QQ/MM]=QMP(n,xl/L-L2,RB,QQ/MM);endcase2P=MPQ.(i,2);b=MPQ(i,3);RA=P-Pb/(L-L2);RB=-PL/(L-L2);QQ二QQ+RA;MM=MM+RA*xl;[QQ,MM]=QMP(n,xl,b,P,MM);[QQ/MM]=QMP(n,xl,L-L2/RB,QQ/MM);c
5、ase3q=MPQ(i,2);c=MPQ(i,3);d=MPQ(i,4);RA=q*(d-c)-q*(d+c)*(d-c)/[2*(L-L2)];RB=-q*(c+d)*(d-c)/[2*(L-L2)];QQ二QQ+RA;MM=MM+RA*xl;[QQMM]二QMQ(n,x「cdq,QQ,MM);[QQ/MM]=QMP(n,xl,L-L2,RB/QQ/MM);endend[QDX,MDX,XQM]=MAX_MIN(xl,QQ,MM);TU_QM(xl,Q.QzMM);验证程序实例1clear;clc;L2=2x=[0,2,4,6]MPQ=[1,-3,2,0330
6、,42,560]XQM二QMYW(x丄2,MPQ);剪力图弯矩图验证程序实例2clear;clc;L2=2x=[0,l,5,6,8]MPQ=[140502108031015]XQM二QMYW(x丄2,MPQ);20剪力團弯矩團左端固定计算机分析程序functionXQM=QMZG(x,MPQ)[n,m]=size(x);L=x(m);X1=[];fori=l:m-lxl=[xl,linspace(x(i),x(i+l),50)];end%离散简支梁MM=zeros(size(xl));QQ=zeros(size(xl));%将离散点的剪力和弯矩数组QQ、MM赋零
7、[m,t]=size(MPQ);[t,n]=size(xl);fori=l:m;switchcase1a=MPQ(i,3);M=MPQ(iz2);RA=0;MA=-M;QQ二QQ+RA;MM=MM+RA*xl+MA;ifa>0&a0&b<=LQQ二QQ+RA;MM=MA+MM+RA*xl;[QQ,MM]=QMP(n,xl,b,P,QQ,MM);endcase3q=MPQ(i,2);c=MPQ(
8、i,3);d=MPQ(i