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

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

2、ing Zhou Li（College of Civil Engineering & Mechanics, XiangTanUniversity, 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,

3、obtained 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 ag

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

5、44+222+4x 4x y y =0显然满足。 （2）体力不计，求得应力分量表达式：212F x =2=-3xyy hy =0xy3F=-2h4y 2 1-h 2（3）由应力分量求解应变分量112F x =(x -y =-3xyE Eh 1y =(y -x =0E 2(1+ 6(1+ F 4y 2xy =1-2xy = E 2Eh h （4）边界条件：y =±h2的主要边界上，应该满足应力边界条件如下：a. 在(y y =±h2=0 xy3F =-2h4y 2 1-h 2=0b.在x =0, x =l 应用圣维南原理，可列出三个积分的应力边界条件如下：(h 2-h 2h

6、2-h 2x x =0, ldy =0（a ）(h 2-h 2z z =0ydy =0，(h 2-h 2x x =lydy =-Fl(bxy x =0, ldy =-F(c对于如图所示矩形板和坐标系，当板内发生上述应力时，由应力边界条件式可知上边、下边无面力；而左边界上受有铅直力；右边界上有按线性变化的水平面合力为一力偶，和铅直面力。所以，能解决悬臂梁在自由端受集中力作用的问题。3 塑性解析解由弹性阶段的应力、应变分量关系可知，矩形截面偏压柱中纵向截面中任意一点的应力状态和应变状态都是相同的。材料为理想弹塑性材料，=0.3，屈服应力为10000MPa 。根据Mises 屈服条件，有：当构件变形

7、进入塑性阶段后，屈服条件：(1-2 2+(2-3 2+(3-1 2-2s 2=0在平面应力状态下，有一个主应力为零，假定3=0 则Mises 屈服条件变为：212-12+2=s 2在直角坐标系中： 1，2=弹性解析解中：y =0x +y2综合以上各式，可得：2x =s2=E ，s ；由理想弹塑性模型（见图3-1），可知，当s 时，当=s时，s 。 图2 理想弹塑性模型弹性解析解中：212F如该处的应力达到材料的x =2=-3xy 选取y=-5这一路径上的点，y h 屈服应力（y =s ），F =x h 360x则有当x=100时，得出临界力F=1.67KN，从而可以得知，在临界条件下，x =0

8、. 04286。 理论解与有限元解的比较假定一组数据：h=10mm, L=100mm,b=1mm F=1KN，E=210000Mpa，=0.3那么可以得到计算结果如下：表1 路径一的各节点理论值(固定端 表2 路径一的各节点有限元值 图3 路径一各点正应力理论解与有限元拟合 图4 路径一各点剪应力理论解与有限元拟合 图5 路径一各点正应变理论解与有限元拟合表3 路径二的各节点理论值(上边界 表4 路径二的各节点有限元值(上边界 图6路径二各点正应力理论解与有限元拟合 图7 路径二各点剪应力理论解与有限元拟合 图8 路径二各点正应变理论解与有限元拟合从上图可以看出，abaqus 模拟的结果与理论

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

10、.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 总结 通过对上述问题进行弹塑性分析和有限元的模拟， 我们可以知道受集中荷载 的悬臂梁的最大应力在固