《材力电算大作业梁的变形matlab.docx》由会员分享，可在线阅读，更多相关《材力电算大作业梁的变形matlab.docx（10页珍藏版）》请在三一文库上搜索。
1、材料力学电算大作业题目名称：梁的变形计算作者班号：作者学号：作者姓名：指导教师：完成时间：评语：成绩（满分10分）： 签名： 时间：1、程序框图：说明： 当a0,b=0时，结构为简支梁； 当a=0,b0时，结构为悬臂梁； 当a0,b0时，结构为外伸梁。2、源程序：clear%清屏clc;%清理变量disp(请输入所要求梁的信息)F=input(请输入集中力F（/N）: );a=input(请输入梁长度a(/m): );b=input(请输入梁长度b(/m): );c=input(请输入集中力作用位置c(/m): );E=input(请输入弹性模量E(/GPa): );I=input(请输入惯性
2、矩I(/cm2): );if (F0|a0|b0|c0|E=0|I(a+b)|(a=0 & b=0) A=你的输入有误; disp(A)elseif F=0|c=a|c=0 B=你的输入有误该梁将不会弯曲; disp(B)else x=0:0.01:(a+b);%构造0到a+b的a+b个点（每隔0.01画一个点）的x坐标轴 if a=0 & b=0 n=length(x); v=zeros(1,n);%v是一个n列的零矩阵 for i=1:(a/0.01+1) if x(i)=c v(i)=-F*(a-c)*x(i)*(a2-x(i).2-(a-c)2)/(6*E*I*105*a);u=-F*
3、c*(a-c)*(2*a-c)/(6*E*I*a);w=-F*(a-c)*(28*a*c-c2)(3/2)/(15*E*I*a);q(i)=F*(a-c)*(a2*x(i)-x(i)3-(a-c)2*x(i)/(6*E*I*a);z(i)= F*(a-c)*(a2-3*x(i)2-(a-c)2)/(6*E*I*a);disp(x处挠度：)disp(q(i)disp(x处转角：)disp(z(i)disp(最大转角)disp(u)disp(最大挠度)disp(w) else v(i)=-F*(a-c)*(a/(a-c)*(x(i)-c).3+(a2-(a-c)2)*x(i)-x(i).3)/(6
4、*E*I*105*a);u=-F*(a-c)*(a+c)/(6*E*I*a);w=- F*(a-c)*(28*a*c-c2)(3/2)/(15*E*I*a);q(i)=F*(a-c)*(a2*x(i)-x(i)3-(a-c)2*x(i)/(6*E*I*a);z(i)= F*(a-c)*(a2-3*x(i)2-(a-c)2)/(6*E*I*a);disp(x处挠度：)disp(q(i)disp(x处转角：)disp(z(i)disp(最大转角)disp(u)disp(最大挠度)disp(w) end end end if a=0 & b=0 n=length(x); v=zeros(1,n);
5、for i=1:(b/0.01+1) if x(i)=c v(i)=-F*x(i).2*(3*c-x(i)/(6*E*I*105); else v(i)=-F*c2*(3*x(i)-c)/(6*E*I*105); end endq(i)=F*(6*a*x(i)-3*x(i)2)/(6*E*I);z(i)=6*F*(a-x(i)/(6*E*I);u=-F*c*c/(2*E*I);w=-F*c*c*(3*b-c)/(6*E*I);disp(x处挠度：)disp(q(i)disp(x处转角：)disp(z(i)disp(最大转角)disp(u)disp(最大挠度)disp(w) end if a=0
6、 & b=0 n=length(x); v=zeros(1,n); if ca for i=1:(a+b)/0.01+1) if x(i)c & x(i)a syms t d=-F*(t-a)*(c-a)*(3*t-a)-(t-a).2)/(6*E*I*105); e=subs(diff(d),t,c); for i=1:(a+b)/0.01+1) if x(i)a & x(i)c g=-F*(c-a)*(c-a)*(3*c-a)-(c-a).2)/(6*E*I*105); v(i)=g+(x(i)-c)*e; end end endv(i)=-F*(a-c)*x(i)*(a2-x(i).2-
7、(a-c)2)/(6*E*I*105*a);u=-F*c*(a-c)*(2*a-c)/(6*E*I*a);w=-F*(a-c)*(28*a*c-c2)(3/2);q(i)=F*(6*a*x(i)-3*x(i)2)/(6*E*I);z(i)=6*F*(a-x(i)/(6*E*I);disp(x处挠度：)disp(q(i)disp(x处转角：)disp(z(i)disp(最大转角)disp(u)disp(最大挠度)disp(w) end figure plot(x,v) title(挠曲线图); ylabel(v /m),xlabel(x /m); end3、模拟1、 算例： 取F=20，a=2,b=0,c=1,E=2,I=2运行结果为：取F=20，a=0,b=3,c=2,E=2,I=2输出结果为：取F=20，a=2,b=2,c=1,E=2,I=2运行结果为：总结：本算法运用迭代公式，先以梁的有关输入长度变量判断梁的类型和梁是否弯曲，然后再根据算法求出梁的某截面挠度和转角，最大挠度和最大转角和挠曲线的形状，本程序是计算梁受到恒定力时的挠度和转角，对于受到分布力的梁，程序类似，故没有写上。