吉林大学材料力学课程设计

1、材料力学课程设计吉林大学材料力学课程设计HZ140TR2后置旅游车底盘车架的静力学分析及强度、刚度计算2011-10-21班级： 机械学号： 4109姓名： 目录材料力学课程设计1、题目：HZ140TR2后置旅游车底盘车架的静力学分析及强度、刚度计算1一、材料力学课程设计的目的1二、材料力学课程设计的任务和要求1、设计计算说明书的要求12、分析讨论及说明部分的要求23、程序计算部分的要求2三、设计题目要求3四、问题求解41.计算各端支座反力42.画出车架的内力图63.画出各截面上弯曲正应力最大值沿轴线方向的变化曲线74.用能量法求解车架的最大挠度fmax及画出车架挠曲线的形状85.按等截面梁重

2、新设计车架截面尺寸9五、程序设计（MATLAB）101.流程图与程序102.程序结果16六、课程设计总结18七、参考资料18材料力学课程设计、 题目：HZ140TR2后置旅游车底盘车架的静力学分析及强度、刚度计算一、 材料力学课程设计的目的本课程设计的目的是在于系统学习完材料力学之后，能结合工程中的实际问题，运用材料力学的基本理论和计算方法，独立的计算工程中的典型零部件，以达到综合运用材料力学的知识解决实际问题的目的。同时，可以使学生将材料力学的理论和现代计算方法及手段融为一体。既从整体上掌握了基本理论和现代的计算方法，又提高了分析问题、解决问题的能力，既是为以前所学知识（高等数学、工程图学、

3、理论力学、算法语言、计算机和材料力学等）的综合应用，又为后继课程（零件、专业课等）打下基础，并初步掌握工程中的设计思想和设计方法，对实际工作能力有所提高。具体的有以下六项：1.使学生的材料力学知识系统化完整化；2.在全面复习的基础上，运用材料力学知识解决工程中的实际问题；3.由于选题力求结合专业实际，因而课程设计可以把材料力学知识与专业需要结合起来；4.综合运用以前所学习的各门课程的知识（高等数学、工程图学、理论力学、算法语言、计算机和材料力学等），使相关学科的只是有机的联系起来；5.初步了解和掌握工程实践中的设计思想和设计方法；6.为后续课程的教学打下基础。二、 材料力学课程设计的任务和要求

4、参加设计者要系统复习材料力学课程的全部基本理论和方法，独立分析、判断设计题目的已知条件和所求问题，画出受力分析计算简图和内力图，列出理论依据并导出计算公式，独立编制计算机程序，通过计算机给出计算结果并完成了设计说明书。、设计计算说明书的要求设计计算说明书是该题目设计思想、设计方法和设计结果的说明，要求书写工整，语言简洁，条理清晰、明确，表达完整。具体内容应包括：1设计题目的已知条件、所求及零件图。2画出构件的受力分析计算简图，按比例标明尺寸、载荷及支座等。3. 静不定结构要画出所选择的基本静定系统及与之相应的全部求解过程。4. 画出全部内力图，并表面可能的各危险截面。5. 各危险截面上各种应力

5、的分布规律图及由此判定各危险点处的应力状态图。6. 各危险点的主应力大小及主平面的位置。7. 选择强度理论并建立强度条件。8. 列出全部计算过程的理论依据、公式推导过程及必要的说明。9. 对变形及刚度分析要写明所用的能量法计算过程及必要的内力图和单位力图。10. 疲劳强度计算部分要说明循环特征，max，min，r，m，a的计算，所查各系数的，,，依据，疲劳强度校核过程及结果，并绘出构件的持久曲线。2、分析讨论及说明部分的要求1. 分析计算结果是否合理，并讨论其原因、改进措施。2. 提出改进设计的初步方案及设想。3. 提高强度、刚度及稳定性的措施及建议。3、程序计算部分的要求1. 程序框图。2.

6、 计算机程序（含必要的说明语言及标识符说明）。3. 打印结果（数据结果要填写到设计计算说明书上）。三、 设计题目要求HZ140TR2后置旅游车底盘车架简化后如下图所示。满载时受重力FA作用，后部受重力FB作用，乘客区均布载荷为q（含部分车身重），梁为变截面梁。计算过程中忽略圆角影响，并把梁抽象为等厚闭口薄壁矩形截面的阶梯梁。材料的弹性模量E，许用应力及有关数据有表7-1给出。本次课程设计选取计算数据FB=4400 N，q=15300 N/m。表7-1 固定数据 1.1 1.6 3.1 1.6 2.1 0.1 0.06 0.120.08 0.11 0.07 0.005210 160 26801.

7、计算前簧固定端C处，前簧滑板D处、后簧固定端F处、后簧滑板G处的支反力。2.画出车架的内力图。3.画出各截面上弯曲正应力最大值沿轴线方向的变化曲线。4.用能量法求出车架最大挠度的值及所发生的截面，画出车架挠曲线的大致形状。5.若壁厚t不变，取h/b=1.5，按等截面梁重新设计车架截面尺寸。四、 问题求解1.计算各端支座反力这个设计题目的力学模型是多支座的超静定梁，从所给的模型可以判断为三次超静定。由于该模型仅在竖直方向上受到外力作用，水平方向上不受外力作用，因此可以把第三个固定端铰支座看作是活动端铰支座，经过适当处理可以把底盘车架的力学模型再度简化为多跨连续梁的力学模型，形成了二次超静定的连续

8、梁力学模型。计算这种梁时最好选用这样的静定基：在中间支座处将梁切开，安上铰链并加上成对的力矩X1、X2、Xm代替被去掉的在相邻两跨之间的约束。本题设计的静定基如下图。简化后的多跨连续梁力学计算模型根据材料力学的知识，对每一个中间支座可以列出一个三弯矩方程，Mn-1ln+2Mnln+ln+1+Mn+1ln+1=-6(wnanln+wn+1an+1ln+1)可列出的方程数目恰好等于中间支座的数目，也就是等于超静定的次数。解方程组求出力矩后，可以把力矩作为载荷，每一跨梁可以看作是静定的简支梁，进而便可求出梁的内力、应力、变形。由于MATLAB具有强大的数值计算功能和图形绘制功能，编程易于实现，本次设

9、计采用MATLAB编程求解问题，具体程序见五、程序设计。将力FA、FB向C、G点简化，得到Mc=-FAl1=-2948Nm，MG=-FBl5=-9240Nm。对跨度l2和l3写出三弯矩方程，Mcl2+2MD(l2+l3)+MFl3=-6(w2a2l2+w3b3l3)对跨度l3和l4写出三弯矩方程，MDl3+2MF(l3+l4)+MGl4=-6(w3a3l3+w4b4l4)a2=b2=l2 2=0.8m, a3=b3=l3 2=1.55m, a4=b4=l4 2=0.8m2=112ql24=5222.4Nm2, 3=112ql34Nm2=37984Nm, 4=112ql44Nm2=5222.4N

10、m将以上数据代入两个三弯矩方程，解方程组可求得弯矩MD=-10388.4Nm, MF=-8790.4Nm下面用简化的模型对每一段梁求解支反力对梁l2求解支反力：由计算简图可列出平衡方程MC=MD+FD1l2-12ql22-MC=0MD=MC+FC-FAl2-12ql22-MD=0代入数据，可得此方程组的解为FC=10269.8NFD1=16890.2N对梁l3求解支反力，同样可以求出：FD2=24230.5NFF1=23199.5N对梁l4求解支反力，同样可以求出：FF2=11959.0NFG=16921.0N综上，各支座支反力如下FC=10270NFD=FD1+FD2=41121NFF=FF

11、1+FF2=35159NFG=16921N2.画出车架的内力图以A点为原点，由1问得到的结果，对车架的每一段求剪力方程，得到第一段到第五段梁的剪力方程（单位：N）分别为：FQ1=-2680FQ2=24419.77-15300xFQ3=65540.47-15300xFQ4=100699-15300xFQ5=4400由此可以画出剪力图对剪力方程积分，利用边界条件，求得第一段到第五段梁的弯矩方程(单位：Nm)分别为：M1=-2680xM2=-7650x2+24420x-20553.M3=-7650x2+65540.5x-131579M4=-7650x2+100699x-335499M5=4400x-

12、41800由此可以画出弯矩图 3.画出各截面上弯曲正应力最大值沿轴线方向的变化曲线最大弯曲正应力max=MymaxIz=MWz 其中Iz=bh3-(b-2t)(h-2t)312，ymax=h2WZ= Izymax=bh3-(b-2t)(h-2t)36h代入数据求得：WZ1=3.925×10-5WZ2=6.2597×10-5WZ3=5.0258×10-5用公式max=MWz分别求得各段梁的最大弯曲正应力函数，画出图形如下4.用能量法求解车架的最大挠度fmax及画出车架挠曲线的形状已经求得各段的弯矩方程，确定各段最大挠度发生截面后在该截面作用单位力1，可以采用图乘法，

13、为了准确地画出挠曲线图，可以采用莫尔积分求解。第一段和第五段最大挠度截面发生在自由端顶点所在截面，中间三段梁的最大挠度截面近似可以认为发生在各段中间，此处给出准确结果的求解方法。采用公式EIv''(x)=M(x)两次积分可以求出各段挠度。v(x)=M(x)EI+Cx+D，转角=v'(x)，代入边界条件，求出了各段挠曲线方程，进而可以求出最大值和所对应的截面，结果如下v1x=0m, v2x=2.11m, v3x=4.27m, v4x=6.61m, v5x=9.5m下面对各截面施加单位力，用莫尔积分公式=M(x)M(x)EIdx分别求出各段对应截面处的最大挠度为f1max=

14、-0.0053m, f2max=0.0011m, f3max=-0.0087m, f4max=0.002m, f5max=-0.0358m综合可知，最大挠度fmax=-0.0358m画出的挠曲线的图形为5.按等截面梁重新设计车架截面尺寸根据前面已经求出的弯矩方程和弯矩图比较各截面最大弯矩，求出最大的弯矩Mmax=10388 NmMmaxWZ< ， 其中WZ=bh3-(b-2t)(h-2t)36h ，h=1.5b取MmaxWZ=，代入解方程，得到b=0.0814m，h=0.1221m 五、 程序设计（MATLAB）1.流程图与程序首先创建一个M函数F_Liang.m:%求解单个梁模型两端支

15、反力function Fr1,Fr2=F_Liang(F1,F2,M1,M2,l,q)A=l,0;0,l;C=M2+q*l2/2-M1+F1*l;M1+q*l2/2-M2+F2*l;X=AC;Fr1=X(1);Fr2=X(2);主程序部分main.m:%材料力学课程设计第一题所有问题的求解%1.以下求解未知力矩%参考数据q=15300,FB=4400;clc;clear;l1=1.1;l2=1.6;l3=3.1;l4=1.6;l5=2.1;FA=2680;q=input('输入q的值（单位N/m）：');FB=input('输入FB的值（单位N）：');clos

16、e all;MC=-FA*l1;MG=-FB*l5;w2=q*l23/12;w3=q*l33/12;w4=q*l43/12;A=2*(l2+l3),l3;l3,2*(l3+l4);C=-6*(w2/2+w3/2)-MC*l2;-6*(w3/2+w4/2)-MG*l4;RE=AC;fprintf('n1.求解未知力矩的结果：n')MD=RE(1);MF=RE(2);fprintf('MD=%.1fNm , MF=%.1fNmnn',MD,MF);%以下求解支反力，需要调用M函数F_Liang.mdisplay('第一段梁求解结果:');FC,FD1

17、=F_Liang(FA,0,MC,MD,l2,q);fprintf('FC=%.1fN , FD1=%.1fNnn',FC,FD1);display('第二段梁求解结果:');FD2,FF1=F_Liang(0,0,MD,MF,l3,q);fprintf('FD2=%.1fN , FF1=%.1fNnn',FD2,FF1);display('第三段梁求解结果:');FF2,FG=F_Liang(0,FB,MF,MG,l4,q);fprintf('FF2=%.1fN , FG=%.1fNnn',FF2,FG);FD=

18、FD1+FD2;FF=FF1+FF2;display('综合得到多跨梁各支座的支反力如下：')fprintf('FC=%.0fNnFD=%.0fNnFF=%.0fNnFG=%.0fNnnn',FC,FD,FF,FG);%2.以下画出内力图%求解剪力方程，画出剪力图FQ1max=-FA;FQ1min=-FA;FQ2max=FC-FA;FQ2min=FC-FA-l2*q;FQ3max=FC-FA-l2*q+FD;FQ3min=FC-FA-l2*q+FD-l3*q;FQ4max=FC-FA-l2*q+FD-l3*q+FF;FQ4min=FC-FA-l2*q+FD-l3

19、*q+FF-l4*q;FQ5min=FC-FA-l2*q+FD-l3*q+FF-l4*q+FG;FQ5max=FQ5min;FQ1=-FA;FQ2=subs(x)FQ2max-q*(x-l1);FQ3=subs(x)FQ3max-q*(x-l1-l2);FQ4=subs(x)FQ4max-q*(x-l1-l2-l3);FQ5=FQ5max;disp('2.第一段到第五段梁的剪力方程（单位：N）分别为：')fprintf('FQ1=%sn',char(vpa(FQ1,7)fprintf('FQ2=%sn',char(vpa(FQ2,7)fprint

20、f('FQ3=%sn',char(vpa(FQ3,7)fprintf('FQ4=%sn',char(vpa(FQ4,7)fprintf('FQ5=%snn',char(vpa(FQ5,7)figure(1);hold on;ezplot(num2str(FQ1),0,l1);ezplot(FQ2,l1,l1+l2);ezplot(FQ3,l1+l2,l1+l2+l3);ezplot(FQ4,l1+l2+l3,l1+l2+l3+l4);ezplot(num2str(FQ5),l1+l2+l3+l4,l1+l2+l3+l4+l5);xlim(0,l1

21、+l2+l3+l4+l5);ylim(-5e4,5e4);plot(l1,l1,FQ1min,FQ2max,'-');plot(l1+l2,l1+l2,FQ2min,FQ3max,'-');plot(l1+l2+l3,l1+l2+l3,FQ3min,FQ4max,'-');plot(l1+l2+l3+l4,l1+l2+l3+l4,FQ4min,FQ5max,'-');plot(0,l1+l2+l3+l4+l5,0,0,'k');text(0.25,FQ1max-0.35e4,num2str(FQ1max);text

22、(l1-0.45,FQ2max+0.35e4,num2str(FQ2max);text(l1+l2-0.45,FQ2min-0.35e4,num2str(FQ2min);text(l1+l2-0.45,FQ3max+0.35e4,num2str(FQ3max);text(l1+l2+l3-0.45,FQ3min-0.35e4,num2str(FQ3min);text(l1+l2+l3-0.45,FQ4max+0.35e4,num2str(FQ4max);text(l1+l2+l3+l4-0.45,FQ4min-0.35e4,num2str(FQ4min);text(8,FQ5max+0.35e

23、4,num2str(FQ5max);xlabel('x /m');ylabel('FQ /N');title('剪力图');%求解弯矩方程，画出弯矩图syms x;M1=int(num2str(-FA),0,x);M1r=subs(M1,x,l1);M2=int(FQ2,l1,x)+M1r;M2r=subs(M2,x,l1+l2);M3=int(FQ3,l1+l2,x)+M2r;M3r=subs(M3,x,l1+l2+l3);M4=int(FQ4,l1+l2+l3,x)+M3r;M4r=subs(M4,x,l1+l2+l3+l4);M5=int(

24、num2str(FQ5max),l1+l2+l3+l4,x)+M4r;M5r=subs(M5,x,l1+l2+l3+l4+l5);x2,M2max=fminbnd(eval('(x)',char(-M2),l1,l1+l2);M2max=-M2max;x3,M3max=fminbnd(eval('(x)',char(-M3),l1+l2,l1+l2+l3);M3max=-M3max;x4,M4max=fminbnd(eval('(x)',char(-M4),l1+l2+l3,l1+l2+l3+l4);M4max=-M4max;disp('

25、第一段到第五段梁的弯矩方程(单位：Nm)分别为：');fprintf('M1=%sn',char(vpa(expand(M1),7);fprintf('M2=%sn',char(vpa(expand(M2),5);fprintf('M3=%sn',char(vpa(expand(M3),6);fprintf('M4=%sn',char(vpa(expand(M4),6);fprintf('M5=%snn',char(vpa(expand(M5),7);fprintf('剪力图和弯矩图如图所示nnn&

26、#39;);figure(2);hold on;box off;ezplot(M1,0,l1);ezplot(M2,l1,l1+l2);ezplot(M3,l1+l2,l1+l2+l3);ezplot(M4,l1+l2+l3,l1+l2+l3+l4);ezplot(M5,l1+l2+l3+l4,l1+l2+l3+l4+l5);plot(0,l1+l2+l3+l4+l5,0,0,'k');xlim(0,l1+l2+l3+l4+l5);ylim(-1.2e4,1e4);plot(l1,l1,0,M1r,'-');plot(l1+l2,l1+l2,0,M2r,'

27、;-');plot(l1+l2+l3,l1+l2+l3,0,M3r,'-');plot(l1+l2+l3+l4,l1+l2+l3+l4,0,M4r,'-');text(l1-0.45,M1r-0.6e3,num2str(M1r);text(x2-0.3,0.6e3,num2str(M2max);text(x2,-0.7e3,'bfdownarrow');text(l1+l2-0.3,M2r-0.6e3,num2str(M2r);text(x3-0.3,M3max+0.6e3,num2str(M3max);text(l1+l2+l3-0.3

28、,M3r-0.6e3,num2str(M3r);text(x4-0.3,M4max+0.6e3,num2str(M4max);text(l1+l2+l3+l4-0.15,M4r-0.6e3,num2str(M4r);title('弯矩图');xlabel('x /m');ylabel('Mz /Nm');%3.以下画出最大弯曲正应力的变化曲线fprintf('3.最大弯曲正应力的变化曲线如图所示nnn');h1=0.1;h2=0.12;h3=0.11;b1=0.06;b2=0.08;b3=0.07;t=0.005;E=210e9;

29、Wz=(b,h)(b*h3-(b-2*t)*(h-2*t)3)/(6*h);Wz1=Wz(b1,h1);Wz2=Wz(b2,h2);Wz3=Wz(b3,h3);sigma1=abs(M1/Wz1);sigma2=abs(M2/Wz2);sigma3=abs(M3/Wz2);sigma4=abs(M4/Wz2);sigma5=abs(M5/Wz3);figure(3);hold on;ezplot(sigma1,0,l1);ezplot(sigma2,l1,l1+l2);ezplot(sigma3,l1+l2,l1+l2+l3);ezplot(sigma4,l1+l2+l3,l1+l2+l3+l

30、4);ezplot(sigma5,l1+l2+l3+l4,l1+l2+l3+l4+l5);plot(0,l1+l2+l3+l4+l5,0,0,'k');xlim(0,l1+l2+l3+l4+l5);ylim(0,2.5e8);Y12=max(subs(sigma1,x,l1),subs(sigma2,x,l1);Y23=max(subs(sigma4,x,l1+l2+l3+l4),subs(sigma5,x,l1+l2+l3+l4);plot(l1,l1,0,Y12,'-');plot(l1+l2+l3+l4,l1+l2+l3+l4,0,Y23,'-&#

31、39;);xlabel('x /m');ylabel('sigma_max /Pa');title('最大弯曲正应力变化曲线');%4.以下求最大挠度%求各段挠度syms C1 D1 C2 D2 C3 D3 C4 D4 C5 D5;v2=int(int(M2*2/(Wz2*h2*E)+C2*x+D2;C2,D2=solve(subs(v2,x,l1),subs(v2,x,l1+l2);v3=int(int(M3*2/(Wz2*h2*E)+C3*x+D3;C3,D3=solve(subs(v3,x,l1+l2),subs(v3,x,l1+l2+l3

32、);v4=int(int(M4*2/(Wz2*h2*E)+C4*x+D4;C4,D4=solve(subs(v4,x,l1+l2+l3),subs(v4,x,l1+l2+l3+l4);theta2=subs(diff(v2,x);theta2C=subs(theta2,x,l1);theta4=subs(diff(v4,x);theta4G=subs(theta4,x,l1+l2+l3+l4);v1=int(int(M1*2/(Wz1*h1*E)+C1*x+D1;theta1=subs(diff(v1,x);theta1C=subs(theta1,x,l1);C1,D1=solve(theta

33、1C-theta2C,subs(v1,x,l1);v5=int(int(M5*2/(Wz3*h3*E)+C5*x+D5;theta5=subs(diff(v5,x);theta5G=subs(theta5,x,l1+l2+l3+l4);C5,D5=solve(theta5G-theta4G,subs(v5,x,l1+l2+l3+l4);v1=subs(v1);v2=subs(v2);v3=subs(v3);v4=subs(v4);v5=subs(v5);%画出挠曲线figure(4);hold on;box off;ezplot(v1,0,l1);ezplot(v2,l1,l1+l2);ezp

34、lot(v3,l1+l2,l1+l2+l3);ezplot(v4,l1+l2+l3,l1+l2+l3+l4);ezplot(v5,l1+l2+l3+l4,l1+l2+l3+l4+l5);plot(0,l1+l2+l3+l4+l5,0,0,'k');xlim(0,l1+l2+l3+l4+l5);ylim(-0.05,0.01);xlabel('x /m');ylabel('v /m');title('挠曲线图');v1x,v1min=fminbnd(eval('(x)',char(v1),0,l1);v2x,v2ma

35、x=fminbnd(eval('(x)',char(-v2),l1,l1+l2);v2max=-v2max;v3x,v3min=fminbnd(eval('(x)',char(v3),l1+l2,l1+l2+l3);v4x,v4max=fminbnd(eval('(x)',char(-v4),l1+l2+l3,l1+l2+l3+l4);v4max=-v4max;v5x,v5min=fminbnd(eval('(x)',char(v5),l1+l2+l3+l4,l1+l2+l3+l4+l5);text(0.2,v1min-0.003

36、,num2str(v1min);text(v2x,v2max+0.003,num2str(v2max);text(v3x,v3min-0.003,num2str(v3min);text(v4x,v4max+0.003,num2str(v4max);text(v5x-1,v5min-0.003,num2str(v5min);fprintf('4.挠曲线图如图所示nn');%用能量法求各段最大挠度M_1=x;M_21=-(l1+l2-v2x)/l2*(x-l1);M_22=(v2x-l1)/l2*(x-l1-l2);M_31=-(l1+l2+l3-v3x)/l3*(x-l1-l2)

37、;M_32=(v3x-l1-l2)/l3*(x-l1-l2-l3);M_41=-(l1+l2+l3+l4-v4x)/l4*(x-l1-l2-l3);M_42=(v4x-l1-l2-l3)/l4*(x-l1-l2-l3-l4);M_5=-x+l1+l2+l3+l4+l5;f11=int(2*M1*M_1/(E*Wz1*h1),x,0,l1);f10=-theta1C*l1;f1=f11+f10;f21=int(2*M2*M_21/(E*Wz2*h2),x,l1,v2x);f22=int(2*M2*M_22/(E*Wz2*h2),x,v2x,l1+l2);f2=f21+f22;f31=int(2*

38、M3*M_31/(E*Wz2*h2),x,l1+l2,v3x);f32=int(2*M3*M_32/(E*Wz2*h2),x,v3x,l1+l2+l3);f3=f31+f32;f41=int(2*M4*M_41/(E*Wz2*h2),x,l1+l2+l3,v4x);f42=int(2*M4*M_42/(E*Wz2*h2),x,v4x,l1+l2+l3+l4);f4=f41+f42;f51=int(2*M5*M_5/(E*Wz3*h3),x,l1+l2+l3+l4,l1+l2+l3+l4+l5);f50=theta4G*l5;f5=f50+f51;f1=double(subs(f1);f2=do

39、uble(subs(f2);f3=double(subs(f3);f4=double(subs(f4);f5=double(subs(f5);A=abs(f1),abs(f2),abs(f3),abs(f4),abs(f5);fmax=max(A);num=find(fmax=A);fmax=eval('f',num2str(num);fmax_x=eval('v',num2str(num),'x');fprintf('各段的最大挠度为：nf1max= %.5fmnf2max= %.6fmnf3max= %.5fmnf4max= %.5f

40、mnf5max= %.4fmnn',f1,f2,f3,f4,f5);fprintf('综合可知车架的最大挠度为：nfmax= -%.5fmn发生的截面位置为：nx= %.5fmnnn',fmax,fmax_x);%5.按等截面梁设计车架截面尺寸syms b;Wzl=Wz(b,1.5*b);h=1.5*b;Yingli=160e6;Mmax=max(abs(M1r),abs(M2r),abs(M3r),abs(M4r),abs(M5r),abs(M2max),abs(M3max),abs(M4max);b=solve(Mmax-Wzl*Yingli,b);b=double(b);b=max(real(b);h=subs(h);fprintf('5.按等截面梁重新设计车架截面尺寸为：nb= %.4fmnh= %.4fmnn',b,h);2.程序结果运行后输入数据15300、4400，回车后输出结果（斜体表示，图形结果见四）：输入q的值（单位N