材料力学上机大作业梁的强度校核Word格式文档下载.docx
《材料力学上机大作业梁的强度校核Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《材料力学上机大作业梁的强度校核Word格式文档下载.docx（38页珍藏版）》请在冰点文库上搜索。
2.点击开始并选择支撑条件为左固定端悬臂梁，点击确定。
如图
3.选择梁的截面形状并输入尺寸参数，点击确定。
4.输入受力条件，点击确定进行计算。
5.程序运行结果如图
计算结果与实际运算结果相同
2.求解《新编材料力学》P246例12-8。
1.运行程序计入主界面，点击开始。
2.选择支撑条件为右外伸梁，点击确定。
3.选择梁的截面形状并输入尺寸，点击确定。
如图，（不妨先设b=10，h=30。
注意，应避免b或h为零以防止出现除数为零的情况使程序出错）
4.输入受力情况及右支架位置并点击确定，如图（如果未能输入右支架位置，可能使程序出错）
5.得到计算结果
计算结果与答案b=30mm，h=60mm一致，不过以截面抗弯系数代替截面具体尺寸。
其中小数点后的微量误差受穷举法所采用数量级的影响。
值得注意的是，由于本程序涉及的截面形状及载荷种类较多，而且载荷位置不固定，因此截面设计和许用载荷计算难以得到最终的具体结果。
时间仓促，为减小工作量，本程序将许用载荷计算简化为梁截面所能承受的最大弯矩的计算。
将截面设计略微简化为截面最小抗弯系数的设计。
验证表明，计算结果同样精确。
[主要算法代码和部分窗体截图]：
部分窗体截图：
Form9窗体内源代码：
PrivateSubCommand1_Click（）‘进入实际工作界面’
Form0.Show
Me.Hide
EndSub
PrivateSubCommand2_Click（）‘若点击“结束”按钮则结束程序’
End
Form0窗体内源代码：
PubliclAsDouble‘声明变量’
PublicqAsDouble
Publicf1AsDouble
Publicf2AsDouble
PublicfAsDouble
Publicm1AsDouble
Publicm2AsDouble
PublicmAsDouble
PublicxAsDouble
PublicxaAsDouble
PublicxbAsDouble
Publicxq1AsDouble
Publicxq2AsDouble
Publicxf1AsDouble
Publicxf2AsDouble
Publicxm1AsDouble
Publicxm2AsDouble
Publicx1AsDouble
Publicx2AsDouble
Publicx3AsDouble
Publicx4AsDouble
Publicx5AsDouble
Publicx6AsDouble
PublicfaAsDouble
PublicfbAsDouble
PublicwAsDouble
PublicbAsDouble
PublichAsDouble
Publicd1AsDouble
Publicd2AsDouble
Publicd3AsDouble
PublicwzAsDouble
PublicsylAsDouble
PublicxylAsDouble
PublicpiAsDouble
PubliczdzhAsDouble
PrivateSubCommand1_Click（）‘进入截面形状选择及尺寸参数输入界面’
Form7.Show
PrivateSubCommand2_Click（）‘若点击“结束”按钮则终止程序’
PrivateSubForm_Load（）
PrivateSubOption1_Click（）‘选择支撑条件’
Option1.Value=True
PrivateSubOption2_Click（）
Option2.Value=True
PrivateSubOption3_Click（）
Option3.Value=True
PrivateSubOption4_Click（）
Option4.Value=True
PrivateSubOption5_Click（）
Option5.Value=True
PrivateSubOption6_Click（）
Option6.Value=True
Form7窗体内源代码：
‘本窗体为截面形状选择及尺寸输入界面’
PrivateSubCommand1_Click（）
IfForm0.Option1.Value=TrueThen‘根据用户选择的不同支撑条件调用不
jz.Show同窗体’
ElseIfForm0.Option2.Value=TrueThen
zxb.Show
ElseIfForm0.Option3.Value=TrueThen
yxb.Show
ElseIfForm0.Option4.Value=TrueThen
zws.Show
ElseIfForm0.Option5.Value=TrueThen
yws.Show
Else
sws.Show
EndIf
PrivateSubCommand2_Click（）
PrivateSubForm_Load（）
PrivateSubOption1_Click（）‘选择不同的截面形状’
Form1（简支梁）窗体内源代码：
‘本窗体完成支撑条件—简支梁部分计算’
DimiAsDouble
DimkAsDouble
pi=3.1415926‘给变量赋值’
l=Val（Form7.l.Text）
IfForm7.Option1.Value=TrueThen‘计算截面抗弯系数’
b=Val（Form7.b.Text）
h=Val（Form7.h.Text）
wz=b*h^2/6
ElseIfForm7.Option2.Value=TrueThen
d1=Val（Form7.d1.Text）
wz=pi*d1^3/32
ElseIfForm7.Option3.Value=TrueThen
d2=Val（Form7.d2.Text）
d3=Val（Form7.D.Text）
wz=pi*d3^3*（1-d2^4/d3^4）/32
ElseIfForm7.Option4.Value=TrueThen
wz=Val（Form7.wz.Text）
f1=Val（f11.Text）‘通过用户输入，给变量赋值’
f2=Val（f21.Text）
xf1=Val（xf11.Text）
xf2=Val（xf21.Text）
m1=Val（m11.Text）
m2=Val（m21.Text）
xm1=Val（xm11.Text）
xm2=Val（xm21.Text）
q=Val（q1.Text）
xq1=Val（xq11.Text）
xq2=Val（xq21.Text）
Ifxf1>
xf2Then‘确定f1，f2的位置前后，便于后续计算’
k=f1
f1=f2
f2=k
k=xf1
xf1=xf2
xf2=k
Ifxm1>
xm2Then‘确定m1，m2的位置前后，便于后续计算’
k=m1
m1=m2
m2=k
k=xm1
xm1=xm2
xm2=k
l=l/1000
xf1=xf1/1000
xf2=xf2/1000
xm1=xm1/1000
xm2=xm2/1000
xq1=xq1/1000
xq2=xq2/1000
fb=（m1+m2+f1*xf1+f2*xf2+q*（xq2-xq1）*（xq1+xq2）/2）/l‘求支反力’
fa=f1+f2+q*（xq2-xq1）-fb
m=0
k=0
Fori=0TolStepl/1000000‘该循环用于求截面最大弯矩’
m=fa*i
Ifxf1<
iThen
m=m-f1*（i-xf1）
Ifxf2<
m=m-f2*（i-xf2）
Ifxm1<
m=m+m1
Ifxm2<
m=m+m2
Ifxq1<
iAndxq2>
m=m-q*（i-xq1）^2/2
Ifxq2<
m=m-q*（xq2-xq1）*（i-（xq2+xq1）/2）
Ifm^2>
k^2Then
k=m
Nexti
m=k‘求得最大截面弯矩’
syl=m/wz‘计算最大正应力’
syl=（syl^2）^0.5*10^3
Form8.Label3.Caption="
实际最大正应力为"
&
syl&
"
Mpb"
‘显示结果’
xyl=Val（Form7.xyl.Text）
zdzh=wz*xyl/10^3‘许用载荷计算’
w=m/xyl‘截面设计计算’
w=（w^2）^0.5
Ifsyl>
xylThen
Form8.Label1.Caption="
由"
>
"
xyl&
vbCrLf&
vbclf&
实际应力大于许用应力，不安全。
可改变梁的载荷使梁截面上最大弯矩小于等于"
zdzh&
N.m。
或者改变梁截面形状和尺寸，使其截面抗弯系数大于等于"
w&
cm^3"
<
实际应力小于许用应力，安全"
Form8.Show
EndSub‘显示结果，结束程序’
PrivateSubCommand2_Click（）‘完成“返回上一级”功能’
‘此后的代码与简支梁类似，仅有些许差别，因此不再注释’
Form2（左悬臂梁）窗体内源代码
PrivateSubLabel23_Click（）
pi=3.1415926
IfForm7.Option1.Value=TrueThen
xa=0
xb=l
f1=Val（f11.Text）
xf2Then
xm2Then
Fori=0TolStepl/1000000
If（l-xf1）<
m=m-f1*（i+xf1-l）
If（l-xf2）<
m=m-f2*（i+xf2-l）
If（l-xm1）<
If（l-xm2）<
If（l-xq2）<
iAnd（l-xq1）>
m=m-q*（i+xq2-l）^2/2
If（l-xq1）<
m=m-q*（xq2-xq1）/2*（i-（2*l-xq1-xq2）/2）
m=k
syl=m/wz
xyl=Val（Form7.xyl.Text）
zdzh=wz*xyl/10^3
w=m/xyl
Form3（右悬臂梁）窗体内源代码
Form4（左外伸梁）窗体内源代码
b=
升级会员 