工程力学创新实践报告.docVIP

石家庄铁道大学创新实践 2009 级 工程力学系 专 业 工程力学 学 号 学生姓名 任 飞 指导教师 郭树起 完成日期 2012年7月2日 变截面梁的理论解与有限元解 摘要 本次创新实践意在通过对选题的分析，相关资料的查阅，问题求解的过程，将力学灵活运用增强创新能力。通过对变截面梁问题进行理论求解和有限元求解，并对两种解法的解进行比较。重温了材料力学，弹性力学的相关知识，进一步加深了运用MATLAB求解有限元问题的能力，最终得材料力学与有限元解的相对差值为1.69%，基本符合要求。 关键词：变截面梁 有限元 MATLAB 力学创新 第一部分：选题及要求 （一）题目：选题33 如图1示意图所示，变截面梁上端固定，下端自由，圆截面直径线性减少到上端部的1/3,下端受集中力，并考虑重力，求其位移。上端直径为30cm,集中力为30KN.选择梁长为2.1m。 （二）题目要求 分别计算题目中结构的理论解，并通过MATLAB编程计算结构的有限元解。对两种情况的解进行分析比较。 （三）选题意义 通过对结构进行理论分析和有限元分析，扩展学习了材料力学，弹性力学的相关知识，锻炼了使用MATLAB求解有限元问题的流程。相关资料的查阅有助于对课程更深刻的认识，培养了探究性学习和创新的能力。 图1 示意图 第二部分：理论求解 （一）按照材料力学理论进行求解 如图2微元分析图所示，图形为截面为变截面梁。截面半径随着L的增大而线性减小，故可采用积分的形式进行位移求解。面积线性变化系数设为。上端截面面积为A1,下端截面为A2。分别考虑P和自重作用下的位移，之后采用叠加法进行组合计算即可。 P单独作用 微元段对应的截面的面积=，则全梁的伸长量为，积分后得到 =1.0491e-005 自重单独作用 微元段对应的截面面积=，对应的圆台的体积为 。与之对应的，则全梁的伸长量为。由于被积函数的形式复杂，采用数值积分的方式用MATLAB进行编程计算结果为 4.2752e-007 图2微元分析图 采用简单叠加法叠加 采用简单叠加法进行叠加，将P作用与自重作用相叠加得出最终的杆端伸长位移为1.0918e-005 附材料力学求解位移采用的MATLAB代码如下： clc clear syms F L E A1 A2 x p g F=30000; L=2.1; E=2.1*10^11; A1=pi*0.15^2; A2=pi*0.05^2; g=9.8; p=7800; l1=F*L/(E*(A1-A2))*(log(A1)-log(A2)) lx=p*g/(3*E)*x*(A2*L/(A2*L+(A1-A2)*x)+(A2*L/(A2*L+(A1-A2)*x))^0.5+1); l2x=@(x)eval(lx); l2=quad(l2x,0,2.1,0.00001) l=l1+l2 % answer l=1.0918e-005 （二）按照弹性力学理论进行求解 材料力学假设的不足 通过一个简单的应力分析来说明材料力学中的假设有不足之处，有相当的近似性。 如图3局部楔形体受力图所示的应力情况，材料力学假设截面正应力均匀分布，且不考虑切应力。对于图示边缘的楔形体，显然只考虑正应力是无法满足平衡条件的，因而材料力学的假设有不足之处，需要采用弹性力学理论进行更符合实际情况的分析。但弹性力学的分析是过程很繁杂，只做初步的探讨和分析。 图3 局部楔形体受力分析图 弹性力学模型的建立和分析 选取图4所示的楔形体模型，作为应力分析模型。 本题中锥角为5.45°误差为-0.16%,则根据材料力学结论可推算出弹性力学近似的位移L3=1.0918e-005*（1+0.0016）= 1.0935e-005 第三部分：MATALB有限元编程求解 有限元前期分析 本题为一端固定约束，一端自由的变截面杆状模型。物体是轴对称物体，约束也是轴对称约束，但杆端力为集中力不满足轴对称荷载的要求。因而采用圣维南原理将下端集中力转化为分布荷载，即可满足轴对称载荷的要求，形成轴对称系统。可以采用轴对称问题的方法进行有限元求解。 轴对称