有限元静力分析基本原理.pptVIP

2001年10月1日 Lb-* 3.3.1 单元贡献矩阵 Definition 由图可见，单元1的单元刚度矩阵方程组的分块表示形式为： 扩阶 2001年10月1日 Lb-* 3.3.2 整体刚度矩阵 节点处的变形谐调条件: 节点处的静力平衡条件： 2001年10月1日 Lb-* 3.3.3 整体分析结果 1.单元贡献矩阵 2.整体刚度矩阵 结构刚度矩阵具有奇异性。因为考虑结构的具体支承情况，即结构具有悬空性。 2001年10月1日 Lb-* 3.4 载荷处理 整体分析时的结构刚度方程组是根据外载荷作用在节点上得出得。如果在单元跨间作用由集中力或分布力，则必须用虚功等效原则（即等效前后载荷在任何虚位移方向上的虚功相等）将此跨间载荷移置到节点上。这一移置工作称为载荷处理。 2001年10月1日 Lb-* 3.5 约束处理 1 引入边界条件 2 约束处理方法 2001年10月1日 Lb-* 3.5.1 边界条件 考虑结构的具体支承情况 2001年10月1日 Lb-* 3.5.2 约束处理方法 由于在整体分析时，没有考虑结构的具体支承情况，因此，结构刚度方程组中的整体刚度矩阵[K]在数学上具有奇异性，即其逆矩阵不存在，也就是说，由此方程不能求得唯一的唯一解。这一现象在力学意义上说，是由于在引入支承条件之前，结构还是一个没有支承的悬空结构引起的。 所谓约束处理正是用结构的实际支承情况对结构刚度方程组进行处理，在数学意义上消除整体刚度矩阵[K]的奇异性；在力学意义上，消除结构的悬空性。使结构刚度方程组有唯一解的处理过程。 常用的处理方法有：消行消列法、置大数法和置1法。 置1法原理： [K]中与已知位移对应行主对角元素置1，其他元素置0。 {P}中与已知位移对应元素置已知的位移置。 中与已知位移对应元素仍按未知量处理。 2001年10月1日 Lb-* 3.6 求解 比较有限元解与解析解，作为一级近似，用一维杆受拉力作用的下列公式计算轴向位移: 2001年10月1日 Lb-* 3.7 后处理 1.确定应变 2.确定应力 单元解： 节点解：绕节点平均法等 3.确定主应力 4、确定等效应力 5、进行强度评价 2001年10月1日 Lb-* 4、ANSYS 平面应力求解验证 问题描述：为说明解平面应力问题的有限元方法，我们给出一个问题详细的解。图5-20所示薄板受表面T=1000MPa的拉力作用，确定节点位移和单元应力。板厚为1mm，E=30×106MPa，μ=0.30。 图5- 20 受拉薄板 平面应力问题 练习1-实体求解验证 用一维杆受拉力作用的公式计算轴向位移： 100个单元，121个节点，总刚2*121=242阶方阵 200个单元，363个节点，总刚3*363=1089阶方阵 2001年10月1日 Lb-* 5、ANSYS 静力学问题实例 （1）二维问题对比分析 （2）带孔板分析-应力集中系数 2001年10月1日 Lb-* （1）二维问题对比分析 平面应变问题 问题描述： 一个20×10 ×300mm的正方形条，底端固定，右端受1000Mpa的均布压力，材料的弹性模量为30×106MPa，试用有限元法进行静力分析。 轴对称问题 问题描述：一个截面为20×10mm的正方形，半径为100mm的圆环，底端固定，内壁面受1000MPa的均布压力，材料的弹性模量为30×106MPa，试用有限元法进行静力分析。 薄板弯曲问题 问题描述：一个20×10 ×1mm的正方形板，左端固定，板面受1000MPa的均布拉力，材料的弹性模量为30×106MPa，试用有限元法进行静力分析。 平面应力问题 问题描述：一个20×10 ×1mm的正方形板，左端固定，右端面受1000MPa的均布拉力，材料的弹性模量为30×106MPa，试用有限元法进行静力分析。 结论：同一几何模型可求解不同问题。 2001年10月1日 Lb-* （2）带孔板分析-应力集中系数 ? 如图所示，一个带有直径50mm小圆孔的200×100×5mm板，一端固定，一端承受30MPa的压力，使用系数。已知材料的弹性模量E=2.0×105MP；泊松比ν=0.3。试计算其应力分布并确定孔边的应力集中系数。 命令流 2001年10月1日 Lb-* 有限元静力分析基本原理 Definition 2001年10月1日 Lb-* 1、弹性力学基础 1.1 弹性力学 与材料力学的区别与联系 1.2 弹性力学基本方程 1.2.1 平衡微分方程-应力体力关系 1.2.2 静力边界条件-应力面力关系 1.2.