材料力学课程设计龙门刨床.doc

材料力学课程设计说明书 -7.5-23龙门刨床班 级：机械5班姓 名： 学 号：指导老师：完成时间：2016.09.19目录一设计的目的、任务及要求二设计题目三题目求解1. 约束力求解2. 校核门架的的强度3. 求门架上加力点的水平、垂直位移四C语言程序部分五CATIA有限元分析验证计算六设计感想七参考文件一材料力学课程设计的目的、任务及要求。1.设计目的本课程设计是在系统学完材料力学课程之后，结合工程实际中的问题，运用材料力学的基本理论和计算方法，独立的计算工程中的典型零部件，以达到综合运用材料力学知识解决工程实际问题的目的。同时可以使学生间材料力学的理论和现代计算方法及手段融为一体，既从整体上掌握了基本理论和现代计算方法，又提高了分析问题、解决问题的能力；既是对以前所学知识（高等数学、工程图学、理论力学、算法语言、计算机和材料力学等）的综合运用，又为后续课程（机械设计、专业课等）的学习打下基础，并初步掌握工程设计思想和设计方法，使实际工作能力提高。具体有以下六项：（1）使所学的材料力学知识系统化、完整化。（2）在系统全面复习的基础上，运用材料力学知识解决工程实际中的问题。（3）由于选题力求结合专业实际，因而课程设计可以把材料力学知识与专业需要结合起来。（4）综合运用以前所学的各门课程的知识，使相关学科的知识有机地联系起来。（5）初步了解和掌握工程实践中的设计思想和设计方法。为后续课程的教学打下基础。2.设计的任务和要求参加设计者要系统股息材料力学课程的全部基本理论和方法，独立分析、判断设计题目的已知条件和所求问题，画出手里分析计算简图和内力图，列出理论依据并到处计算公式，独立编制计算程序，通过计算机给出计算结果，并完成设计计算说明书。设计计算说明书是该题目设计思路、设计方法和设计结果的说明，要求书写工整，语言简练，条理清晰、明确，表达完整。二设计题目1.已知某龙门刨床门架示意图如图1所示，可简化为图2所示的钢架，尺寸如图所示。危险工况有最大切削力F，另外有两种力矩Me1、Me2，作用位置及作用方式如图所示。门架材料为灰铸铁（HT250）。（1）校核门架的强度。（取安全系数n=3）（2）求门架上加力点的水平、垂直位移。已知数据：F=230kN Me1=10kNm Me2=45kNm。查手册知灰铸铁HT250的E=100GPa，G=40GPa。 图1 图2 图3 图42.题目分析：（1）结构中存在两个空间的固定端，和一个空间的封闭结构，所以是12次超静定。（2）求解时应运用书11.5节中利用结构对称、载荷对称、平面-空间系统的知识进行求解。（3）利用力法正则方程，分别求解F、Me1、Me2作用下的多余约束力。（4）对横梁及立柱进行校核。（5）利用图形互乘法求解位移。三题目求解1.约束力求解（1）当F单独作用时该结构是关于y轴是对称结构，F单独作用时为对称载荷，截面反对称内力素都为零。其次，此结构为平面-空间系统，所以在任意截面上作用于结构平面内的内力也为零。所以，只存在一个作用于x-z平面的弯矩。将结构沿着z-y平面截开，取左半部分为静定基，分别施加多余约束X1、X2于横梁上。令X1=X2=1，画出力矩图如下： 图1 图2 图3由二次超静定力法正则方程可知： 11X1+12X2+1F2=0 21X1+22X2+2F2=0对于非圆截面立柱，利用插值法求得1=0.1793，2=0.0.1845，所以t柱=1h1b13-2h2b23=7.6910-4m4y梁=b1h1312-b2h2312=3.1410-5m411=1McEy梁+McxGt柱=2.4310-7m22=1McEy梁+McxGt柱=2.1710-7m由位移互等定理可得：12=21=McxGt柱=2.6010-8m1F2=-McxGt柱=-1.8010-3m2F2=-1McEy梁-McxGt柱=-8.3910-3m带入方程组，解得 X1=3.29kNm（方向与单位力方向相同） X2=38.24kNm(方向与单位力方向相同)（2）当Me1单独作用时该结构是关于y轴是对称结构，Me1单独作用时为反对称载荷，截面对称内力素都为零。进行受力分析，固定端不存在垂直于结构平面的作用力，所以将结构沿着z-y平面截开，只存在一个作用于y方向的剪力。取左半部分为静定基，分别施加多余约束X3、X4于横梁上。令X3=X4=1，画出力矩图如下： 由二次超静定力法正则方程可知： 33X3+34X4+3Me12=0 43X3+44X4+4Me12=0Iz梁=h1b1312-h2b2312=1.3810-5m4Iz柱=h1b1312-h2b2312=3.7310-4m433=1Mc1EIz梁+McEIz柱=6.7610-8m44=1Mc1EIz梁+McEIz柱=5.9910-8m由位移互等定理可得：34=43=McEIz柱=7.72*10-9m3F=McEIz柱=-6.4310-5m4F=1 Mc1EIz梁+McEIz柱=-7.1710-4m带入方程组，解得 X3=-421.51N(方向与单位力方向相反) X4=12.02kN(方向与单位力方向相同)（3）当Me2单独作用时该结构是关于y轴是对称结构，Me2单独作用时为对称载荷，截面反对称内力素都为零。进行受力分析，固定端不存在平行于结构平面的作用力，所以将结构沿着z-y平面截开，只存在一个作用于水平面的弯矩。取左半部分为静定基，分别施加多余约束X5、X6于横梁上。令X5=X6=1，画出力矩图如下：由二次超静定力法正则方程可知： 55X5+56X6+5Me22=0 65X5+66X6+6Me22=0 由图可知，显然5F=6F= 0，而55，56，65，66都不为零，欲使方程组成立，则有 X5=X6=0 Nm2. 校核门架的强度。（1）顶部梁架的校核由前面计算可知，顶部梁架受到X1，X3的作用，为斜弯曲的情况，危险截面在左侧端面，力矩图如图所示：Mzmax=-X30.6=252.91Nm Mymax=X1=3.29kNm 经过分析，灰铸铁的抗拉性能最差，所以截面p点受到两个方向的最大拉应力，为危险点。 p=MzmaxyIz梁+MymaxzIy梁=11.58MPa 查资料得HT250的b=250Mpa，n=3,所以=bn=83.33MPa 由第一强度理论， p 所以，下端横梁是不安全的。（3）立柱的校核原理同下端横梁的校核，立柱受到 F2，Me12,Me22, X1， X2，X3以及X4的作用,为两弯、一扭和拉伸的复杂组合。弯矩，扭矩以及轴力图如下：弯矩图：扭矩图： 轴力图： 由以上图可知，危险截面在端面。 校核A点： A点只存在正应力。 FN=X3+X4=11.60kNA=0.30.4-0.260.36=0.0264m2Mxmax=F20.8-Me22=69.50kNmMzmax=-Me12+X40.6+X30.6=1.96kNmIx柱=h1b1312-h2b2312=5.8910-4m4max=MxmaxzIx柱+MzmaxxIz柱+FNA=38.20106MPa由第一强度理论，max,所以，A点安全。 校核B点： B点同时存在正应力和切应力。 By=MzmaxxIz柱+FNA=0.94MPa由教科书P74页薄壁杆件扭转切应力计算公式可知：My=F20.6-X1-X2=27.47kNm=My2min=27.4710320.40.30.02=5.72MPa max=By2+By22+2=6.21MPa由第一强度理论，max,所以，B点安全。校核C点：cy=MxmaxzIx柱+FNA=37.70MPa=My2min=5.72MPa max=Cy2+Cy22+2=38.55MPa由第一强度理论，max,所以，C点安全。综上，立柱安全。3. 求门架上加力点的水平、垂直位移。（1）求水平位移求Z方向位移：在门架上加力点沿着Z方向施加一个单位力，力矩图如图： 由之前计算得到的内力，分别画出力矩图： x柱=b1h1312-b2h2312=5.8910-4m4运用图形互乘法求解：z=iMciEIi+iMxciGIti =12F20.6323EIy梁 -X20.6223EIy梁 + 12F20.8323EIx柱-Me220.8223EIx柱+F20.620.8GIt柱-X10.60.8GIt柱-X20.60.8GIt柱=0.39m (方向与单位力方向相同) 求X方向位移：在门架上加力点沿着X方向施加一个单位力，力矩图如图：运用图形互乘法求解： x=iMciEIi=Me120.8212EIz柱+X30.820.612EIZ柱+X40.820.612EIz柱=-1.6810-2m(方向与单位力方向相反)（2）求水平位移求Y方向位移：在门架上加力点沿着Y方向施加一个单位力，力矩图如图： 运用图形互乘法求解： y=iMciEIi=-Me120.6212EIz梁+23X40.6312EIZ梁-Me120.80.6EIz柱-0.60.80.6X3EIz柱+0.60.80.6X4EIz柱=1.6710-4m(方向与单位力方向相同)四C语言程序#include#define A E*Ix /定义变量#define B E*Iy#define C E*Iz1#define D E*Iz2#define H G*It#define J F/2#define K Me1/2#define N Me2/2void main()int i;float h1,h2,b1,b2,p1,p2,pA,pB,pC,p=250/3;double E,G,a,l,Me1,Me2,F,N1,N2,I5,X1, X2,X3,X4,x,y,z,Ix,Iy,Iz1,Iz2,It;E=100e9; /为变量赋值G=40e9; a=0.6; l=0.8; Me1=1.0e4; Me2=4.5e4; F=2.3e5;N1=0.1793; N2=0.1845; X1=3.29e3;X2=3.824e4; X3=-0.42151e3; X4=1.202e4;for(i=0;i3)I+i=N1*h1*b1*b1*b1-N2*h2*b2*b2*b2;elseI+i=(h1*b1*b1*b1-h2*b2*b2*b2)/12;printf(I%d=%.2en,i,Ii);printf(nInput the p1n); /为p1赋值scanf(%f,&p1);if(p1p) /与许用应力比较，判断是否安全printf(n顶部梁架安全。n);else printf(n顶部梁架不安全。n);printf(nInput the p2n); /为p2赋值scanf(%f,&p2);if(p2p) /与许用应力比较，判断是否安全printf(n下端横梁在M安全。n);else printf(n下端横梁在M不安全。n);printf(nInput the pA,pB,pCn); /为pA,pB,pC赋值scanf(%f,%f,%f,&pA,&pB,&pC); if(pAp&pBp&pCp) /与许用应力比较，判断是否安全printf(n立柱安全。n);else printf(n立柱不安全。n);printf(n门架上加力点在X、Y、Z轴上的位移依次为x、y、zn);Ix=I1; Iy=I2; Iz1=I4; /为变量赋值Iz2=I3; It=I5;x=K*l*l/(2*C)-X4*a*l*l/(2*C)+X3*a*l*l/(2*C); /代入公式求解x方向位移printf(nx=%.2en,x);y=-K*l*a/C-K*a*a/(2*D)-X3*a*l*a/C+X4*a*l*a/C+X4*a*a*a/(3*D); /代入公式求解y方向位移printf(y=%.2en,y);z=J*l*2*l*l/(2*3*A)+J*a*2*a*a/(2*3*B)+J*l*a*a/H-X2*l*a/H-X1*l*a/H-N*l*l*2/(3*A)-X2*a*a*2/(3*B); /代入公式求解z方向位移printf(z=%.2en,z);五CATIA有限元分析验证计算通过CATIA建模，进行有限元分析，可以验证上述计算结果的真实性。首先建模，设置材料，修改材料属性。 施加外力，观察变形。 图一 变形图生成应力图像： 图二 应力张量图 图三 范米塞斯应力图由以上验证可知，底部横梁左端面应力最大，也最危险，这与前面的计算结果吻合。六设计感想原以为不太费力的一次课程设计却花费了这么长的时间，真的没想到。没真正入手之前，永远不知道自己有多少个细节没有掌握，永远不知道知识点有多少个漏洞。我在完成课程设计的过程中了解到了一般零件的设计校核过程，使自己对实际工程问题有了一定的认识。在设计过程中