2-有限元静力基本原理详解.ppt

2.6 结构有限元法分析的目的及适用条件 5、静力分析分析实例 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、进行强度评价 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阶方阵 工程应用三步曲 前处理：定类型，画模型，设属性，分网格。 求 解：添约束，加载荷，查错误，求结果。 后处理：列结果，绘图形，显动画，下结论。 2.4 ANSYS结构分析实例 2.4.1 自重作用下的悬臂梁分析 2.4.2 角速度作用下的轮子分析 2.4.3 板梁混合结构分析 2.4.4 扭矩作用下的空心轴分析 鸟巢命令流 2001年10月1日 Lb-* 二、结构静力有限元分析 1、材料力学的局限 2、 弹性力学基础 3、平面应力问题的有限元分析原理入门 4、ANSYS静力分析实践入门 2001年10月1日 Lb-* 1、材料力学的局限 带孔板单向拉伸问题的应力确定。 应力集中系数 名义应力 2001年10月1日 Lb-* 2、 弹性力学基础 2.1 弹性力学基本方程 2.2 弹性力学求解 2.3 主应力 2.4 强度理论及相当应力 2.5 弹性力学问题简化 2001年10月1日 Lb-* 2.1 弹性力学基本方程 （1）平衡微分方程—应力体力关系 （2）静力边界条件—应力面力关系 （3）几何方程—应变位移关系 （4）物理方程—应力应变关系 2001年10月1日 Lb-* Z 总和后整理便得到X方向的平衡方程: 同理得到x、y方向的平衡方程: X Y Z 其中X、Y、Z 为物体的体力分量。 （1）平衡微分方程——应力体力关系 x y z o 2001年10月1日 Lb-* 设斜面ACD为边界面，其外法线n的方向为(l1,l2,l3)，面积为ΔS，边界外力分量为(px，py，,pz)，则三角形ABC、 ABD 、 BCD的面积分别为ΔS在各相应方向上的投影。 l1ΔS, l2ΔS, l3ΔS, N x y z o （2）静力边界条件——应力面力关系 由x方向的平衡得到： XNΔS = l1ΔSσx+l2ΔSτyx+l3ΔSτzx 即 XN = l1σx