悬臂梁在均布荷载作用下有限元分析

1、悬臂梁承受集中荷载作用问题的弹塑性分析何方平 邹里 （湘潭大学 土木工程与力学学院，湖南 湘潭 411105）摘要 本文针对曲杆在水平力作用下的受力性能，结合弹性力学基本方程和塑性力学中Mises屈服条件，得到了弹性阶段应力、位移之间的关系，以及材料发生塑性变形时，处于临界状态点的应力、应变值。同时，利用有限元分析软件ABAQUS，进行了数值模拟，分析结果与理论值吻合较好，证明所建立的有限元模型是合理的。关键词：悬臂梁；集中荷载THE ELASTIC-PLASTIC ANALYSIS OF THE CANTILEVER BEAM UNDER concentrated loadHe Fang-P

2、ing Zhou Li（College of Civil Engineering & Mechanics, XiangTan University, Xiangtan 411105, China）【Abstract】This article in view of the force performance of CANTILEVER BEAM UNDER concentrated load, combined with elastic mechanics basic equations and the plastic mechanics Mises yield conditions, obta

3、ined the elastic stage between stress and displacement, and the relationship between material happen plastic deformation, a critical state points of stress and strain value. At the same time, the finite element analysis software ABAQUS, the numerical simulation and analysis results and a good agreem

4、ent with the theoretical value, show that the established finite element model is reasonable.Keywords: CANTILEVER BEAM concentrated load 题目：试考察应力函数能满足相容方程，并求出应力分量（不记体力），画出例题3-2图所示矩形体边界上的面力分布（在次要边界上表示出面力的主矢量和主矩），指出该应力函数所能解决的问题。 图1 1 弹性力学解（1）考察相容条件，将应力函数代入相容方程=0显然满足。（2）体力不计，求得应力分量表达式： （3）由应力分量求解应变分量（4

5、）边界条件：a. 在的主要边界上，应该满足应力边界条件如下： b.在应用圣维南原理，可列出三个积分的应力边界条件如下： （a） ， (b) (c) 对于如图所示矩形板和坐标系，当板内发生上述应力时，由应力边界条件式可知上边、下边无面力；而左边界上受有铅直力；右边界上有按线性变化的水平面合力为一力偶，和铅直面力。所以，能解决悬臂梁在自由端受集中力作用的问题。3 塑性解析解由弹性阶段的应力、应变分量关系可知，矩形截面偏压柱中纵向截面中任意一点的应力状态和应变状态都是相同的。材料为理想弹塑性材料，屈服应力为10000MPa。根据Mises屈服条件，有：当构件变形进入塑性阶段后，屈服条件：在平面应力状

6、态下，有一个主应力为零，假定则Mises屈服条件变为： 在直角坐标系中：弹性解析解中：=0综合以上各式，可得：由理想弹塑性模型（见图3-1），可知，当时，；当时，。图2 理想弹塑性模型弹性解析解中： 选取y=-5这一路径上的点，如该处的应力达到材料的屈服应力（），则有当x=100时，得出临界力F=1.67KN，从而可以得知，在临界条件下，。理论解与有限元解的比较假定一组数据：h=10mm, L=100mm,b=1mm F=1KN，E=210000Mpa，=0.3那么可以得到计算结果如下：表1 路径一的各节点理论值(固定端)编号xy/MPa/MPA1100-5600000.028572100-4

7、480054.000.022863100-3360096.000.017144100-22400126.000.011435100-11200144.000.00571610000150.00071001-1200144.00-0.00571481002-2400126.00-0.01142991003-360096.00-0.017143101004-480054.00-0.022857111005-60000-0.028571表2 路径一的各节点有限元值编号xy/MPa/MPA1100-558920.000190.0282292100-44713.652.5960.0221263100-3

8、3535.293.5040.0166974100-22356.8122.7240.0115665100-11178.4140.2560.005834610000145.80.000012571001-1185.6139.968-0.0057881002-2371.2122.472-0.0116891003-3556.893.312-0.01687101004-4742.452.488-0.02229111005-59280.00023-0.02803 图3 路径一各点正应力理论解与有限元拟合 图4 路径一各点剪应力理论解与有限元拟合 图5 路径一各点正应变理论解与有限元拟合表3 路径二的各节点

9、理论值(上边界)编号xy/MPa/MPA15-530000.00143210-560000.00286320-5120000.00571430-5180000.00857540-5240000.01143650-5300000.01429760-5360000.01714870-5420000.02 980-5480000.022861090-5540000.025711195-5570000.02714 表4 路径二的各节点有限元值(上边界)编号xy/MPa/MPA15-5296.1-4.52435E-160.00144210-5592.2-3.25815E-160.002896320-51

10、211.76-2.14842E-160.005788430-51817.64-1.05703E-160.008681540-52349.62.90567E-180.011573650-529371.11613E-160.014466760-53524.42.21109E-160.017359870-54176.483.32787E-160.020251 980-54773.124.52772E-160.0231441090-55513.45.59642E-160.0260351195-55819.76.71073E-160.027321图6 路径二各点正应力理论解与有限元拟合图7 路径二各点剪

11、应力理论解与有限元拟合图8 路径二各点正应变理论解与有限元拟合从上图可以看出，abaqus模拟的结果与理论结果比较吻合，通过表1与表2、表3与表4的数据，可以得到路径一上的的理论值绝对值平均值为36000，模拟平均值为35460，两者误差为0.13%。的理论平均值为90，模拟平均值为87.55636，两者误差为2.7%，正应变的理论平均值绝对值为0.015584，有限元模拟值为0.015351，两者误差为1.5%。路径二上的的理论平均值为3000，模拟平均值为2969.722，两者误差为1.0%。的理论平均值为0，模拟平均值为1.13919E-16，两者误差很小，应变的理论平均值为0.0142

12、85，有限元模拟值为0.01445，两者误差为1.1%。分析误差原因可能是因为有限元分析本来就是近似分析带来的系统误差等，但是误差均在允许范围之内，所以abaqus模拟值是合理的。4 弹塑性有限元分析(1) 定义单元类型：通过查阅abaqus的单元库对单元的性质和应用范围进行了解,选择八结点双向二次平面应力四边形单元，缩减积分。(2) 定义材料参数：主要输入的是弹性模量，泊松比，假定弹性模量为210GPa，泊松比为0.3。(3) 建立模型：该模型为平面应力模型，通过abaqus建立二维几何平面应力模型，在abaqus中生成有限元模型。 (4) 划分网格利用abaqus的网格划分工具，设置网格，采用结构方式划分网格，划分网格后的模型如图2所示。 有限元网格划分(5) 加载数据在悬臂梁自由端施加一个大小为1kN的集中力，设定加载时间为1秒，分1个子步，每个子步1秒。(6) 求解(7) 后处理求解后，abaqus得到了结构在当前边界条件下受力的详细情况，如各节点和单元的应力应变值，受力后结构的变形情况，等等。现列举几个重要结果如下。 应力变形变形云图 应变云图5总结通过对上述问题