《工程力学数值方法》作业解说.doc

《工程力学数值方法》 课程作业 姓名：李艳 学号： 专业: 道路与铁道工程 指导教师： 院系： 2014年6月29日 目录1 问题概述 1 2 圆孔孔口应力集中问题的理论计算 2 3 圆孔孔口应力集中的ANSYS有限元分析 4 3.1 建立有限元模型 4 3.2 施加荷载并求解 4 3.3 查看分析结果并验证结果的正确性 4 3.4 用理论解对数值解进行验证并得出结论 7 4 圆孔孔口应力集中的边界元分析 8 4.1边界单元间接法的基本原理 8 4.2边界单元间接法的数值方法 9 4.3 孔口应力集中问题 11 5 总结 13 1 问题概述 在工程结构设计中，常常会根据需要在构件上进行开孔，由于孔口的存在势必会破坏金属材料的连续性，引起主应力弯曲绕行，应力作用线在孔口附近密集弯曲，造成局部应力增大的现象。由圣维南原理可知，在远离小孔的地方，孔口局部应力集中的影响将消失。孔口附近的应力远大于无孔时的应力，也远大于距孔口较远处的应力，这种现象称为孔口应力集中。应力集中是由于开口后的应力扰动引起的, 最大与最小应力一般发生在孔边上，且应力扰动主要发生在距孔边1.5倍孔口尺寸的范围内。孔口应力集中与孔口形状有关，圆形孔口的应力集中程度比较低，故在工程上被广泛使用。 在此，研究小孔口问题，即孔口的尺寸远小于弹性体的尺寸，并且孔边距弹性体的边界比较远(约大于1.5倍的孔口尺寸)，以排除弹性体边界条件对孔口应力分布的影响。孔口应力集中的存在会使材料的脆性加大，构件的强度降低，尤其会严重影响交变应力下构件的持久疲劳极限。因此，有必要对孔口应力集中问题进行深入的研究和探讨，为工程施工提供理论指导和依据。 本文首先从理论上对孔口应力集中问题进行分析，然后利用有限元软件ANSYS14.5该问题进行模拟分析,并用理论解对ANSYS计算的数值解进行验证,并将理论解与数值解进行分析比较,总结出孔口应力集中问题的特征。 2 圆孔孔口应力集中问题的理论计算 图2-1所示一矩形薄板只在左右两边受有均布的拉力，板中有一圆孔，半径为r，矩形薄板的尺寸远大于孔口尺寸，且孔边距矩形板边界较远。 图2-1 圆孔的矩形薄板 可按“小孔口问题”求解，采用极坐标坐标原点取在圆孔中心，利用弹性力学的方法可得到该矩形薄板应力状态的基尔斯解答如下: 为了研究孔口应力集中情况，在此分情况进行讨论： （1）沿孔边()的环向应力为： 而径向应力，其分布如图-2（）所示。 （2）沿y轴() 的横截面面上的环向应力为： 其分布规律如图-2（）所示。 （3）沿x轴() 纵向截面上的环向应力为： 其分布规律如图-2（）所示。 （4）孔边最大应力出现在处，即： 最小应力出现在处，即： 小圆孔的应力集中系数为： (a) (b) 图2-2 孔口附近应力分布 3 圆孔孔口应力集中的ANSYS有限元分析 ANSYS是一种大型的通用有限元分析软件软件主要包括三个部分前处理模块、分析计算模块后处理模块。 假设矩形板长为0mm，宽为0mm，圆孔半径为20mm，材质为钢板，采用线性静力分析方法。利用大型通用有限元软件ANSYS14.进行分析计算，其具体步骤如下: 3.1 建立有限元模型 (1) 定义单元类型。采用的单元类型为Plane82。 (2) 生成几何模型。由于结构和荷载的对称性，取1/4结构进行分析计算，建立一个长为0mm，宽为0mm的正方形，以左小角为圆心建立半径为20mm的圆，进行减运算得到有限元分析模型。 (3) 定义材料属性。确定薄板的弹性模量E=2.06×10，泊松比为0.3，密度7850kg/m3。 (4) 划分网格。每条直线边平均划分为10段，圆弧段平均划分为20段。分别采用三角形单元与四边形单元进行网格划分，如图3-1所示。 3.2 施加荷载并求解 (1) 定义约束并施加荷载。固定薄板左侧边线，约束左侧边线上所有节点的x方向自由度；在薄板右侧边施加均布荷载，荷载大小为100N/m。 (2) 设置分析选项并求解。选择分析类型为静力分析，并进行求解。 a 分析模型 b 三角形单元 c 四边形单元 3-1 1/4薄板分析模型及有限元网格划分模型 3.3 查看分析结果并验证结果的正确性 划分单元时采用了三角形和四边形单元两种划分方法，相应产生两种分析结果，括号前为三角形