结构力学课程设计1.docVIP

一、课程设计题目: 外伸梁在多个集中力偶载荷作用下的剪力图、弯矩图的实现。 二、课程设计要求： ① 给出相应题目的计算机编译源程序以及相应的程序说明； ② 给出相应题目的算例； ③ 结果以截图方式呈现。 三、源程序 ：（注：程序由MATLAB语言编写） clear L=input(输入外伸梁两支座间长度,L=); %输入格式： L1=input(输入外伸梁左端长度，L1=); L2=input(输入外伸梁右端长度，L2=); x=input(输入集中力偶大小,f=); [n m]=size(x); xf=input(集中力偶作用点,xf=);%力的坐标向量 j=input(集中力偶方向,j=);%力的夹角向量 function XQM=QMDWL(x,L1,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)); [PN,t]=size(MPQ); [t,n]=size(x1); for i=1:PN switch MPQ(i,1) case 1 M=MPQ(i,2); a=MPQ(i,3); RA=-M/(L2-L1); RB=-RA; if a0 aL MM=QMM(n,x1,a,M,MM); end if a==0 MM=MM+M; end case 2 P=MPQ(i,2); b=MPQ(i,3); LL=L2-L1; bb=b-L1; RA=P*(LL-bb)/LL; RB=P*bb/LL; [QQ,MM]=QMP(n,x1,L1,-RA,QQ,MM); [QQ,MM]=QMP(n,x1,L2,-RB,QQ,MM); if b0 bL [QQ,MM]=QMP(n,x1,b,P,QQ,MM); end if b==0 QQ=QQ-P; MM=MM-P*x1; end case 3 q=MPQ(i,2); c=MPQ(i,3); d=MPQ(i,4); b=(c+d)*0.5; P=(d-c)*q; LL=L2-L1; bb=b-L1; RA=P*(LL-bb)/LL; RB=P*bb/LL; [QQ,MM]=QMP(n,x1,L1,-RA,QQ,MM); [QQ,MM]=QMP(n,x1,L2,-RB,QQ,MM); [QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM); end end [QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM); TU_QM(x1,QQ,MM); disp(梁的有限元分析结果) disp(位置-----------剪力----------弯矩) function MM=QMM(n,x1,a,M,MM) %集中力偶对弯矩贡献的子函数 for j=1:n if x1(j)==a n1=j; end end MM(n1:n)=MM(n1:n)+M; function [QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM) %求剪力和弯矩极值的子函数 XQM=[x1,QQ,MM]; [Qmax,i]=max(QQ); Q1=[Qmax,x1(i)]; [Qmin,i]=min(QQ); Q2=[Qmin,x1(i)]; [Mmax,i]=max(MM); M1=[Mmax,x1(i)]; [Mmin,i]=min(MM); M2=[Mmin,x1(i)]; disp(剪力极值及位置) QDX=[Q1;Q2] disp(弯矩极值及位置) MDX=[M1;M2] t1=findobj(0,Tag,text31); str=num2str(Q1(1)); set(t1,String,str); t2=findobj(0,Tag,text39); str=num2str(Q1(2)); set(t2,String,str); t3=findobj(0,Tag,text32); str=num2str(Q2(1)); set(t3,String,str); t4=findobj(0,Tag,text40); str=num2str(