材料力学剪力图弯矩图绘制(有详细的程序)教材.pdfVIP

材料力学剪力图弯矩图绘制（有详细的程序）

说明：

输入变量:

分段数组x

分段点一般在集中力，集中力偶作用出和分布载荷的起末端。

载荷数组MPQ

若梁上的外载荷总数为PN，则用PN行四列的数组MPQ储存载荷，数组MPQ第一列代表载荷

的类型：1为集中力偶，2为集中力，3为分布载荷，第二列代表载荷的大小，第三列代表

集中力，集中力偶或者分布载荷左端与简支梁左端的距离，第四列代表均匀载荷右端与简支

梁左端的距离，当载荷为集中力或者集中力偶时，第四列为0.

符号规定

集中力和均匀载荷向下为正，向上为负，集中力偶顺时针为正，逆时针为负。

输出变量:

内力数组XQM

如果梁被分为NN-1段，则内力数组XQM为NN行，三列的数组，第一列代表梁的横截面的位

置，第二列代表剪力，第三列代表弯矩。

剪力极值及位置QDX

QDX是一个二行二列的数组，第一列代表极值所在的位置，第二列代表极值

弯矩极值及位置MDX

MDX是一个二行二列的数组，第一列代表极值所在的位置，第二列代表极值

1.子程序

1.1集中力偶对弯矩贡献的子函数QMM

1.2集中力对剪力和弯矩贡献的子函数QMP

1.3分布载荷对剪力和弯矩贡献的子函数QMQ

1.4求剪力和弯矩极值的子函数MAX_MIN

1.5绘制剪力图和弯矩图的子函数TU_QM

2.计算分析程序

2.1简支梁QMDJ

2.2左端固定悬臂梁QMDXZ

2.3右端固定悬臂梁QMDXY

2.4左端外伸梁QMDWZ

2.5右端外伸梁QMDWY

2.6两端外伸梁QMDWL

1

1.子程序

1.1集中力偶对弯矩贡献的子函数QMM

functionMM=QMM(n,x1,a,M,MM)

forj=1:n

ifx1(j)==a

n1=j;

end

end

MM(n1:n)=MM(n1:n)+M;

1.2集中力对剪力和弯矩贡献的子函数QMP

function[QQ,MM]=QMP(n,x1,b,P,QQ,MM)

forj=1:n

ifx1(j)==b;

n1=j;

end

end

QQ(n1:n)=QQ(n1:n)-P;

MM(n1:n)=MM(n1:n)-P*(x1(n1:n)-b);

1.3分布载荷对剪力和弯矩贡献的子函数QMQ

function[QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM)

forj=1:n

ifx1(j)c

QQ(j)=QQ(j)-q*(x1(j)-c);

MM(j)=MM(j)-0.5*q*(x1(j)-c)^2;

end

ifx1(j)d

QQ(j)=QQ(j)+q*(x1(j)-d);

MM(j)=MM(j)+0.5*q*(x1(j)-d)^2;

end

end

1.4求剪力和弯矩极值的子函数MAX_MIN

function[QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM)

XQM=[x1,QQ,MM];

[Qmax,i]=max(QQ);

Q1=[Qmax,x1(i)];

[Qmin,i]=min(QQ);

Q2=[Qmin,x1(i)];

[Mmax,i]=max(MM);

M1=[Mmax,x1(i)];

[Mmin,i]=min(MM);

2

M2=[Mmin,x1(i)];

disp(剪力极值及位置)

QDX=[Q1;Q2]

disp(弯矩极值及位置)

MDX=[M1;M2]

t1=findobj(0,Tag,text31);

str=num2str(Q1(1));

set(t1,String,str);

t2=findobj(0,Tag,text39);

str=num2str(Q1(2));

set(t2,String,str);

t3=findobj(0,Tag,text32);

str=num2str(Q2(1));

set(t3,String,str);

t4=findobj(0,Tag,text40);

str=num2str(Q2(2))