后置旅游车底盘车架的静力学分析及强度、刚度计算

1、材料力学课程设计材料力学课程设计hz140tr2后置旅游车底盘车架的静力学分析及强度、刚度计算 专业：工程机械专业 姓名： 学号： 目录材料力学课程设计1一、材料力学课程设计的目的1二、材料力学课程设计的任务和要求1、设计计算说明书的要求12、分析讨论及说明部分的要求23、程序计算部分的要求2三、设计题目2、求支座反力3、画出车架的内力图5、画出各截面上弯曲正应力最大值沿轴线方向的变化曲线74、求出最大绕度的值，画出车架扰曲线的大致形状85、按等截面梁重新设定截面尺寸12四、程序计算131、语言程序132、matlab程序14五、设计体会161、体会162，建议17参考文献17材料力学课程设计

2、一、 材料力学课程设计的目的本课程设计的目的是在于系统学习完材料力学之后，能结合工程中的实际问题，运用材料力学的基本理论和计算方法，独立的计算工程中的典型零部件，以达到综合运用材料力学的知识解决实际问题的目的。同时，可以使学生将材料力学的理论和现代计算方法及手段融为一体。既从整体上掌握了基本理论和现代的计算方法，又提高了分析问题、解决问题的能力，既是为以前所学知识（高等数学、工程图学、理论力学、算法语言、计算机和材料力学等）的综合应用，又为后继课程（零件、专业课等）打下基础，并初步掌握工程中的设计思想和设计方法，对实际工作能力有所提高。具体的有以下六项：1.使学生的材料力学知识系统化完整化；2

3、.在全面复习的基础上，运用材料力学知识解决工程中的实际问题；3.由于选题力求结合专业实际，因而课程设计可以把材料力学知识与专业需要结合起来；4.综合运用以前所学习的各门课程的知识（高等数学、工程图学、理论力学、算法语言、计算机和材料力学等），使相关学科的只是有机的联系起来；5.初步了解和掌握工程实践中的设计思想和设计方法；6.为后续课程的教学打下基础。二、 材料力学课程设计的任务和要求参加设计者要系统复习材料力学课程的全部基本理论和方法，独立分析、判断设计题目的已知条件和所求问题，画出受力分析计算简图和内力图，列出理论依据并导出计算公式，独立编制计算机程序，通过计算机给出计算结果并完成了设计说

4、明书。、设计计算说明书的要求设计计算说明书是该题目设计思想、设计方法和设计结果的说明，要求书写工整，语言简洁，条理清晰、明确，表达完整。具体内容应包括：1设计题目的已知条件、所求及零件图。2画出构件的受力分析计算简图，按比例标明尺寸、载荷及支座等。3. 静不定结构要画出所选择的基本静定系统及与之相应的全部求解过程。4. 画出全部内力图，并表面可能的各危险截面。5. 各危险截面上各种应力的分布规律图及由此判定各危险点处的应力状态图。6. 各危险点的主应力大小及主平面的位置。7. 选择强度理论并建立强度条件。8. 列出全部计算过程的理论依据、公式推导过程及必要的说明。9. 对变形及刚度分析要写明所

5、用的能量法计算过程及必要的内力图和单位力图。10. 疲劳强度计算部分要说明循环特征，max，min，r，m，a的计算，所查各系数的，,，依据，疲劳强度校核过程及结果，并绘出构件的持久曲线。2、分析讨论及说明部分的要求1. 分析计算结果是否合理，并讨论其原因、改进措施。2. 提出改进设计的初步方案及设想。3. 提高强度、刚度及稳定性的措施及建议。3、程序计算部分的要求1. 程序框图。2. 计算机程序（含必要的说明语言及标识符说明）。3. 打印结果（数据结果要填写到设计计算说明书上）。三、 设计题目hz140tr2后置旅游车底盘车架简化后如下图所示。满载时受重力fa作用，后部受重力fb作用，乘客区

6、均布载荷为q（含部分车身重），梁为变截面梁。计算过程中忽略圆角影响，并把梁抽象为等厚闭口薄壁矩形截面的阶梯梁。材料的弹性模量e，许用应力及有关数据有表7-1给出。其他数据选取设计计算数据第一组中的fb=4610n, q=12000n/m 。表7-1 固定数据 1.1 1.6 3.1 1.6 2.1 0.1 0.06 0.12t /me/gpa/mpa/n0.08 0.11 0.07 0.005210 160 26801.计算前簧固定端c处，前簧滑板d处、后簧固定端f处、后簧滑板g处的支反力。2.画出车架的内力图。3.画出各截面上弯曲正应力最大值沿轴线方向的变化曲线。4.用能量法求出车架最大挠度

7、的值及所发生的截面，画出车架挠曲线的大致形状。5.若壁厚t不变，取h/b=1.5，按等截面梁重新设计车架截面尺寸。、求支座反力为了求解该题，需要用到连续梁问题中的三弯矩方程。为减小跨度很大直梁的弯曲变形和应力，常在其中间安置若干中间支座的这类结构称为连续梁。连续梁中是在两端的固定端铰支座和活动铰支座之间加上若干个活动铰支座组成的力学模型。在hz140tr2后置旅游车底盘车架的力学模型中，该模型是三次超静定结构。该模型仅在竖直方向上受到外力作用，水平方向上不受外力作用，因此可以把第三个固定端铰支座看作是活动端铰支座，这样就把底盘车架的力学模型再度简化为连续梁的力学模型。连续梁力学模型中的静定次数

8、等于连续梁力学模型中间活动铰支座的个数，所以该模型是二次超静定的连续梁力学模型。解决静不定结构的方法是把超静定结构化作静定基本机构，然后再在静定基本结构上添加外荷载及多余约束力得到静不定结构的相当系统。然后在相当系统中列出平衡方程，列出位移或者受力的补充方程，联立求解支座反力。fbm0m1m3m2m2m1设想将每个中间支座上的梁切开，并装上铰链，将连续梁变成若干个简支梁，每个简支梁都是一个静定基本系，这相当于把每个支座上梁的内约束解除，即将其内力弯矩m、m作为多余约束力，则每个支座上方的铰链两侧截面上需加上大小相等、方向相反的一对力偶矩，与其相应的位移是两侧截面的相对转角。于是多余约束处的变形

9、协调条件是梁中间支座处两侧截面的相对转角为零。fa图简化后的hz140tr2后置旅游车底盘车架的力学模型图经过简化的连续梁模型中，标出由于梁的内力解除需要维持与原来系统相当的变形而产生的力矩m、m，其中中间支座之间的力矩m1,m2为多余约束。把f1和f2分别移到支座c和支座g处，则: m0=-fal1=-2680n1.1m=-2948nm; m3=-fal5=-4610n2.1m=-9681nm取静定基的每个跨度皆为简支梁。这些简支梁在原来的外载荷作用下的弯矩图如下图所示。为了计算时便于利用三弯矩方程进行计算，设cd段为l1=l2=1.6m；df段为l2=l3=3.1m；fg段为l3=l4=1

10、.6m 。图3 弯矩图由此计算得：w1=20l212ql2dl=226q(l2)3=112ql13=4096nm2同理可得：w2=112ql23=29791nm2 w3=112ql33=4096nm2 有上图可知，各个力矩图形的形心位置a1=l12=0.8m ; a2=l22=1.55m ; b2=l22=1.55m ; b3=l32=0.8m由三弯矩方程：mn-1ln+2mnln+ln+1+mn+1ln+1=-6(nanln+n+1bn+1ln+1)对跨度l1和l2写出三弯矩方程为：m0l1+2m1l1+l2+m2l2=-6(1a1l1+2b2l2) 对跨度l2和l3写出三弯矩方程为：m1l

11、2+2m2l2+l3+m3l3=-6(2a2l2+3b3l3） 且l1=l3由方程和联立解得：m1=-8179.6nm m2=-6469.6nm求的m1和m2后，相当系统的各部分受力如图: 图4 各部分的受力图根据简单的解简单的简支梁的方法，列出竖直方向的受力平衡方程和一个支点处的力矩平衡方程结果可以求出fc,fd1,fd2,ff2,ff1,fg的值，从而得到固定端c处，前簧滑板d处，后簧固定端f处、后簧滑板g处的支反力。fc=9010.25n ;fd=32021.36n ;ff=25641.3n ;fg=16217.1n、画出车架的内力图相当系统中的各部分都可以看成静定梁，而且载荷和端截面上

12、的弯矩（多余约束力）都是已知的。对每一跨度都可以求出支反力和弯矩图，这些剪力图和弯矩图与原超静定系统的对应部分的剪力图和弯矩图相同，所以这些剪力图和弯矩图连起来就是连续梁的剪力图和弯矩图。先写出各段的的剪力方程组：fs1=-fa fs2=-fa+fc-qx-l0 fs3=fd2-qx-l0-l1 fs4=ff2-qx-l0-l1-l2fs5=fb 带入数据：得到剪力方程：fs1=-2680 0x1.1fs2=-12000x+19530 1.1x2.7fs3=-12000x+51551 2.7x5.8fs4=-12000x+77193 5.8x7.4fs5=4610 7.4x9.5 根据建立方程

13、书写matlab程序，使用matlab画出剪力图，结果如下图： 图5 剪力图根据各段的剪力方程利用matlab积分写出弯矩方程，再利用matlab画出弯矩图，结果如下图：y0=-2680x 0x1.1 y1=-6000x2+19530x-17171 1.1x2.7y2=-6000x2+51551x-103630 2.7x5.8 y3=-6000x2+77193x-252350 5.8x7.4 y4=4610x-43795 7.4x9.5使用matlab精确画出弯矩图的形状，主要是可以观察弯矩图的形状，以此对车架模型的受力情况做出正确的判断。图6 弯矩图为了使读图更方便，使用列表方式把关键的支点

14、剪力发生突变的位置和弯矩的最值点和弯矩的拐点列表写出。关键点上的剪力和弯矩大小列表点/m01.11.62752.74.2959剪力/n-2680-26800-12870,191510弯矩/ nm0-2948-1278.5-81807102点/m5.86.43287.49.5剪力/n-18048,75930-11607,46104610弯矩/ nm-6470-4067.6-96810、画出各截面上弯曲正应力最大值沿轴线方向的变化曲线根据公式： iz=bh3-b-2th-2t312可以求出各段的惯性矩的值： iz1=1.962510-6m3 iz2=3.77583310-6 m3 iz3=2.76

15、416710-6 m3 根据上一题的弯矩方程和公式max=myiz求出各段的最大弯曲正应力方程，有matlab画出弯曲正应力最大值沿轴线方向的变化曲线。图7 最大正应力图特殊点最大正应力值列表：特殊点位置m01.1-1.1+1.62752.74.2959应力值107n/m207.514.692.0313.0011.29特殊点位置m5.86.4287.4-7.4+9.5应力值107n/m210.286.4615.3819.2604、求出最大绕度的值，画出车架扰曲线的大致形状求出车架上特殊点的挠度，其中最大的就是车架最大挠度所在截面。为了便于计算，作出每一个载荷作用下的弯矩图，然后利用图乘法和叠加

16、原理求其总和。根据上图，做出每个载荷单独作用时的弯矩图：fa单独作用时mfamax=fa(l1+l2)fb单独作用时mfbmax=fb(l4+l5)nc单独作用时ng单独作用时cd部分均布载荷单独作用时mqcdmax=ql222df段均布载荷单独作用时mqdfmax=ql328fg段单独作用时mqfgmax=ql422（1）求a点挠度mfamax=l1+l2在a端加单位力，弯矩图如上图所示。由图乘法可知：=i=1nwimcieii fa单独作用下a点挠度：fa1=1eiz112l1l1l1+l2mfamax23l1l1+l2mfamax+1eiz2l2l1l1+l2mfamaxl1+l22l1

17、+l2mfamax+1eiz212l22l1+l2mfamaxl1+23l2l1+l2mfamax+12l1mfamax23mfamax fb单独作用下a点挠度：fa2=1eiz212l3mfamax13mfbmaxnc单独作用下a点挠度：fa3=1eiz212l2ncl2l1+2l23l1+l2mfamax+12l3ncl223mfamaxng单独作用下a点挠度fa4=1eiz212l3ngl413mfamaxcd部分均布载荷单独作用时a点挠度：fa5=1eiz213mqcdmaxl2l1+3l24l1+l2mfamax12l3mfamax23mqcdmaxdf段均布载荷单独作用时a点挠度：

18、fa6=1eiz212mfamax23mqcdmaxfg段均布载荷单独作用时a点挠度：综上得： fa=fa1+fa2+fa3+fa4+fa5+fa6+fa7=5.3mm同理可以计算得：cd中点e的绕度fe=-1.1mmdf中点0的绕度fo=9.2mmfg的中点k的绕度fk=-2.1mmb点的绕度fb=42.3mm由以上计算，可以得到车架在b端得挠度最大 fmax=42.3mm 车架挠曲线如下图所示，单位mm.扰曲线的大致形状5、按等截面梁重新设定截面尺寸重新设计截面梁的截面尺寸时的主要判断依据是判断截面上的最大弯曲正应力是否小于许用应力。根据弯矩图可以判断最大弯矩发生的截面，最大弯矩已知，想要

19、改变最大正应力的大小，就需要根据hb=1.5重新设计截面的尺寸。根据弯曲正应力的强度条件max=|m|maxwz由弯矩图可知，最大弯矩发生在后簧滑板g处的截面：mzmax=9681nm |m|maxwz wz=(bh3-(b-2t)(h-2t)3)/12h/2=bh26-(b-2t)(h-2t)36h h=0.1238 , b=0.0825四、 程序计算1、语言程序#include#define fa 2680main()double l0=1.1,l1=1.6,l2=3.1,l3=1.6,l4=2.1;double a1=l1/2,a2=l2/2,b2=l2/2,b3=l3/2;double

20、 m0,m1,m2,m3,w1,w2,w3,h,g;double fc,fd,ff,fd1,fd2,ff1,ff2,fg,fb,q,fz;double lz1,lz2,lz3;double b1=0.06,h1=0.1,b2=0.08,h2=0.12,b3=0.07,h3=0.11,t=0.005; /*定义变量的值*/int n;lz1=(b1*h1*h1*h1-(b1-2*t)*(h1-2*t)*(h1-2*t)*(h1-2*t)/12;lz2=(b2*h2*h2*h2-(b2-2*t)*(h2-2*t)*(h2-2*t)*(h2-2*t)/12;lz3=(b3*h3*h3*h3-(b3-

21、2*t)*(h3-2*t)*(h3-2*t)*(h3-2*t)/12; /*计算个截面的惯性矩*/printf(lz1=%enlz2=%enlz3=%en,lz1,lz2,lz3);for(n=1;n=10;n+)printf(enter q,fb:n); /*输入参变量fb和q*/scanf(%lf,%lf,&q,&fb);printf(q=%e,fb=%en,q,fb);w1=q*l1*l1*l1/12;w2=q*l2*l2*l2/12;w3=q*l3*l3*l3/12; /*计算各段力矩所围的面积*/m0=-fa*l0;m3=-fb*l4; /*计算m0和m3的值*/h=-6*(w2*a

22、2/l2+w3*b3/l3)-m3*l3;g=-6*(w1*a1/l1+w2*b2/l2)-m0*l1; /*简化方程的计算和表达*/m1=(h*l2-g*2*(l1+l2)/(l2*l2-4*(l1+l2)*(l1+l2);m2=(h-m1*l2)/(2*(l1+l2); /*l1=l3简化计算得到m1和m2*/fd1=-(-m0-q*l1*l1/2+m1)/l1; /*-m0=fa*l0*/fc=fa+q*l1-fd1;ff1=-(-m3-q*l3*l3/2+m2)/l3; fg=fb+q*l3-ff1;fd2=-(-q*l2*l2/2+m1-m2)/l2;ff2=q*l2-fd2;fd=

23、fd1+fd2;ff=ff1+ff2;fz=fc+fd+ff+fg;printf(the result is:n);printf(m1=%enm2=%en,m1,m2);printf(fc=%enfd=%enff=%enfg=%enfz=%fn,fc,fd,ff,fg,fz);printf(fd1=%enfd2=%enff2=%enff1=%enfz=%fn,fd1,fd2,ff2,ff1,fz);c语言运行结果：2、matlab程序、图5 剪力图的matlab程序：fa=2680;fb=4610;fc=9010.25;fd=32021.36;ff=25641.3;fg=16217.1;l0=

24、1.1;l1=1.6;l2=3.1;l3=1.6;l4=2.1;q=12000;syms xy0=-2680;y1=-2680+9010.25-12000*(x-1.1)y2=19151-12000*(x-2.7)y3=7592.9-12000*(x-5.8)y4=4610; !有参数求出各段的弯矩方程，由方程求关键点plot(0,1.1,1.1,2.7,2.7,5.8,5.8,7.4,7.4,9.5,9.5,-2680,-2680,6330,-12870,19151,-18048,7592,-11607,4610,4610,0)title(剪力图)xlabel(x/m)ylabel(fs/k

25、n)、图六 弯矩图的matlab程序：x0=0:0.1:1.1;y0=-2680*x0;x1=1.1:0.05:2.7;y1=-6000*x1.2+19530*x1-17171;x2=2.7:0.05:5.8;x3=5.8:0.05:7.4;x4=7.4:0.05:9.5;y2=-6000*x2.2+51551*x2-103630;y3=-6000*x3.2+77193*x3-252350;y4=4610*x4-43795;plot(x0,y0,x1,y1,x2,y2,x3,y3,x4,y4)title(弯矩图)xlabel(x/m)ylabel(m/kn.m)、利用弯矩方程求最大正应力的方程

26、的matlab程序：iz1=1.9625*10-6;iz2=3.775833*10-6;iz3=2.764167*10-6;syms x;y0=-2680*x;y1=-6000*x2+19530*x-17171;y2=-6000*x2+51551*x-103630;y3=-6000*x2+77193*x-252350;y4=4610*x-43795;s0=abs(y0)*0.05/iz1;s1=abs(y1)*0.06/iz2;s2=abs(y2)*0.06/iz2;s3=abs(y3)*0.06/iz2;s4=abs(y4)*0.55/iz3;s1=expand(s1)s0=expand(s

27、0)s2=expand(s2)s3=expand(s3)s4=expand(s4)、图7 最大正应力图clear alliz1=1.9625*10-6;iz2=3.775833*10-6;iz3=2.764167*10-6;syms x;y0=-2680*x;y1=-6000*x2+19530*x-17171;y2=-6000*x2+51551*x-103630;y3=-6000*x2+77193*x-252350;y4=4610*x-43795;s0=abs(y0)*0.05/iz1;s1=abs(y1)*0.06/iz2;s2=abs(y2)*0.06/iz2;s3=abs(y3)*0.0

28、6/iz2;s4=abs(y4)*0.055/iz3;s0s1s2s3s4fplot(316398554352266229317632*abs(x)/4633822111315839,0,1.1)hold onfplot(1770887431076116955136*abs(6000*x2 - 19530*x + 17171)/111442920025707125,1.1,2.7)hold onfplot(1770887431076116955136*abs(6000*x2 - 51551*x + 103630)/111442920025707125,2.7,5.8)hold onfplot(

29、1770887431076116955136*abs(6000*x2 - 77193*x + 252350)/111442920025707125,5.8,7.4)hold onfplot(3246626956972881084416*abs(4610*x - 43795)/163167619923179225,7.4,9.5)title(最大正应力图)xlabel(x/m)ylabel(n/m)五、 设计体会1、体会通过材料力学课程设计首先我体验到了自己知识的缺乏，仅是学习课本会做课本上的习题是解决不了工程实际问题的。至少目前是这样，目前所学的知识仅仅是可以去做做课后习题。在做课程设计的时候，我使用了很多工具，我查了很多资料，利用网络和图书管的网上资源可以查到自己想要找的参考书籍。另外在材料力学课程设计的过程中，解决工程实际问题的一般方法是把实际问题抽象成力学模型，并通过对力学模型的分析建立平衡方程，从而得到计算结果。 总体感觉在材料力学课程设计的时候涉及的知识方面很广，从cad画出一般的受力图，到matlab根据方程画出内力图，并使用matlab绘图工具进行编辑。从使用力学知识进行力学分析，到使用word编辑文档。每一项想要做好都不是意见简单