材料力学上机大作业(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')

title('挠曲线图')

grid on;

三、使用方法

运行代码

输入圆杆长度（单位：m）

输入圆杆直径（单位：m）

输入弹性模量（单位：GPa）

输入集中力大小（单位：N）（向下为正，若无请输入0）

输入集中力作用位置（单位：m）（若无请输入0）

输入弯矩大小（单位：N*m）（逆时针为正，若无请输入0）输入弯矩作用位置（单位：m）（若无请输入0）

输出挠曲线图

四、运行实例

【实例1】

圆杆同时受集中力与弯矩作用，输入、输出见下图。

【实例2】

圆杆只受集中力作用，输入、输出见下图

【实例3】

圆杆只受弯矩作用，输入、输出见下图