哈工大材力上机 matlab,材料力学上机大作业(matlab编)

一、可实现课题

在如图所示的悬臂梁中，杆件为圆杆。杆长为L，直径为D，材料弹性模量为E。输入集中力F大小，作用点a，弯矩M，作用点b，即可求得悬臂梁的挠度曲线图。

二、程序代码

clear all

disp('请给定材料信息'); %输入材料信息

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')