材料力学课程设计--右端外伸梁.doc

材料力学课程设计 MATLAB编程部分说明书 编程题目：端外伸梁 设计者班级：指导教师：成绩： 年 月 日 clear; clc; L2=2 x=[0,2,4,6] MPQ = [2.5 2 0 4 3 -3 2 0 3.5 4 6 0] XQM=QMYW(x,L2,MPQ); 三，次主程序 function XQM=QMYW(x,L2,MPQ) [n,m]=size(x); L=x(m); x1=[]; for i=1:m-1 x1=[x1,linspace(x(i),x(i+1),50)]; end MM=zeros(size(x1)); QQ=zeros(size(x1)); [m,t]=size(MPQ); [t,n]=size(x1); for i=1:m switch MPQ(i,1); case 1 a=MPQ(i,3); M=MPQ(i,2); RA=-M/(L-L2); RB=-RA; QQ=QQ+RA; MM=MM+RA*x1; if a= 0aL MM=QMM(n,x1,a,M,MM); end %if a==0 % MM=MM+M; %end [QQ,MM]=QMP(n,x1,L-L2,-RB,QQ,MM); case 2 P=MPQ(i,2); b=MPQ(i,3); RA=P-P*b/(L-L2); RB=P*b/(L-L2); QQ=QQ+RA; MM=MM+RA*x1; if b0bL [QQ,MM]=QMP(n,x1,b,P,QQ,MM); end [QQ,MM]=QMP(n,x1,L-L2,-RB,QQ,MM); case 3 q=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*(d-c)*(d+c)/[2*(L-L2)]; QQ=QQ+RA; MM=MM+RA*x1; [QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM); [QQ,MM]=QMP(n,x1,L-L2,-RB,QQ,MM); end end [QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM); TU_QM(x1,QQ,MM); 四．运行结果 QDX = 0 4.5000