材力电算大作业梁的变形matlab.docx

材料力学电算大作业题目名称：梁的变形计算作者班号：作者学号：作者姓名：指导教师：完成时间：评语：成绩（满分10分）：签名：时间：1、程序框图：说明：当a≠0,b=0时，结构为简支梁；当a=0,b≠0时，结构为悬臂梁；当a≠0,b≠0时，结构为外伸梁。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(请输入惯性矩I(/cm2): );if (F0|a0|b0|c0|E=0|I=0|c(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)*(a^2-x(i).^2-(a-c)^2)/(6*E*I*10^5*a);u=-F*c*(a-c)*(2*a-c)/(6*E*I*a);w=-F*(a-c)*(28*a*c-c^2)^(3/2)/(15*E*I*a);q(i)=F*(a-c)*(a^2*x(i)-x(i)^3-(a-c)^2*x(i))/(6*E*I*a);z(i)=F*(a-c)*(a^2-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+(a^2-(a-c)^2)*x(i)-x(i).^3)/(6*E*I*10^5*a);u=-F*(a-c)*(a+c)/(6*E*I*a);w=-F*(a-c)*(28*a*c-c^2)^(3/2)/(15*E*I*a);q(i)=F*(a-c)*(a^2*x(i)-x(i)^3-(a-c)^2*x(i))/(6*E*I*a);z(i)=F*(a-c)*(a^2-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)endend endif a==0 b~=0 n=length(x); v=zeros(1,n);for i=1:(b/0.01+1)if x(i)=c v(i)=-F*x(i).^2*(3*c-x(i))/(6*E*I*10^5);else v(i)=-F*c^2*(3*x(i)-c)/(6*E*I*10^5);endendq(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)endif a~=0 b~=0 n=length(x); v=zeros(1,n);if cafori=1:((a+b)/0.01+1)if x(i)=c v(i)=-F*(a-c)*x(i)*(a^2-x(i).^2-(a-c)^2)/(6*E*I*a*10^5);elseif x(i)c x(i)=a v(i)=-F*(a-c)*(a/(a-c)*(x(i)-c).^3+(a^2-(a-c)^2)*x(i)-x(i).^3)/(6*E*I*a*10^5);else r=tan(F*c*(a-c)*(a+c)/(6*E*I*a*10^5));