您好,欢迎访问三七文档
当前位置:首页 > 临时分类 > 材料力学上机大作业(matlab编)
一、可实现课题在如图所示的悬臂梁中,杆件为圆杆。杆长为L,直径为D,材料弹性模量为E。输入集中力F大小,作用点a,弯矩M,作用点b,即可求得悬臂梁的挠度曲线图。二、程序代码clearalldisp('请给定材料信息');%输入材料信息L=input('圆杆长度L(/M)=');D=input('圆杆直径D(/M)=');E=input('弹性模量E(/GPa)=');I=double(D^4*3.14/32);disp('请给定受力情况');%输入受力情况F=input('切向集中力大小F(/N)=');a=input('切向集中力作用位置(/M)=');M=input('弯矩大小M(/N*M)=');b=input('弯矩作用位置(/M)=');x1=0:0.01:a;%F引入的挠度vx1=(-F*x1.^2*3*a+F*x1.^3)*(1/(6*E*10^9*I));x2=a:0.01:L;vx2=(-F*a.^2*3*x2+F*a.^3)*(1/(6*E*10^9*I));v11=[vx1,vx2];x11=[x1,x2];x3=0:0.01:b;%M引入的挠度vx3=(-M*x3.^2)*(1/(2*E*10^9*I));x4=b:0.01:L;vx4=(-M*b*x4+M*0.5*b.^2)*(1/(E*10^9*I));x22=[x3,x4];v22=[vx3,vx4];v33=v22+v11;%叠加plot(x11,v33),xlabel('x/M'),ylabel('v(x)/M')abLFMtitle('挠曲线图')gridon;三、使用方法运行代码输入圆杆长度(单位:m)输入圆杆直径(单位:m)输入弹性模量(单位:GPa)输入集中力大小(单位:N)(向下为正,若无请输入0)输入集中力作用位置(单位:m)(若无请输入0)输入弯矩大小(单位:N*m)(逆时针为正,若无请输入0)输入弯矩作用位置(单位:m)(若无请输入0)输出挠曲线图四、运行实例【实例1】圆杆同时受集中力与弯矩作用,输入、输出见下图。【实例2】圆杆只受集中力作用,输入、输出见下图【实例3】圆杆只受弯矩作用,输入、输出见下图
本文标题:材料力学上机大作业(matlab编)
链接地址:https://www.777doc.com/doc-7358900 .html