第四章空间问题的有限元.pptVIP

第1页，共34页，编辑于2022年，星期二 4.1 空间问题 （1）单元为块体形状。常用单元：四面体单元、长方体单元、直边六面体单元、曲边六面体单元 、轴对称单元。（2）结点位移3个分量。（3）基本方程比平面问题多。3个平衡方程，6个几何方程，6个物理方程。 第2页，共34页，编辑于2022年，星期二 4结点四面体单元：是空间问题最简单的单元，也是常应变、常应力单元，可以类似平面问题三结点三角形单元进行分析。 8结点长方体单元：可以类似平面四结点矩形单元进行分析。 8结点直边六面体单元：可以类似平面四结点任意四边形等参元分析 。 20结点曲边六面体单元：等参单元,可以类似平面八结点曲边四边形等参元进行分析 。 轴对称单元：一平面单元绕一对称轴旋转形成的空间问题。只需在rz平面划分网格，就像平面问题xy平面中的网格一样，这样这类空间问题可以得到简化。 （环向位移等于零） 第3页，共34页，编辑于2022年，星期二 4.2 四结点四面体单元 1、位移模式 单元结点位移向量 位移函数 第4页，共34页，编辑于2022年，星期二 将上式中第1式应用于4个结点，则 由此可解出代定常数a1~a4再代回到式（4-3）的第1式，可得形函数u: 第5页，共34页，编辑于2022年，星期二 编号约定：当沿i,j,m的方向转动时，n在大拇指所指的方向 第6页，共34页，编辑于2022年，星期二 采用同样的方法，可得 单元位移： 第7页，共34页，编辑于2022年，星期二 2、单元应变 将单元中位移（4-11）代入上式 常量 第8页，共34页，编辑于2022年，星期二 3、单元应力 弹性矩阵[D]: 第9页，共34页，编辑于2022年，星期二 代入单元应变计算公式，整理后： 其中[S]为应力矩阵，且： 常量 第10页，共34页，编辑于2022年，星期二 4、单元刚度矩阵 分块矩阵的形式 第11页，共34页，编辑于2022年，星期二 式中子矩阵[krs]为3×3的矩阵 ： 5、等效结点荷载 体积力与表面力的计算公式与平面三角形单元公式 相似，可以采用静力等效原则简化计算。 第12页，共34页，编辑于2022年，星期二 4.2 八结点六面体单元 第13页，共34页，编辑于2022年，星期二 (Tri-linear functions) 第14页，共34页，编辑于2022年，星期二 4.3 二十结点六面体单元 一、形函数 二、位移模式 第15页，共34页，编辑于2022年，星期二 ANSYS 中 Solid45 第16页，共34页，编辑于2022年，星期二 4.4 空间轴对称问题的有限元法 对空间轴对称问题，常采用圆柱坐标系。r表示径向坐标，z表示轴向坐标，任一对称面为rz面。在有限元分析时，可采用轴对称的环形单元进行。环形单元 可以是任何平面单元，本节以三角形单元为例。 第17页，共34页，编辑于2022年，星期二 1、位移模式 轴对称问题的环向位移恒等于零，径向r位移与轴向z位移不等于零。对于图示情形，依照平面问题的三角形单元分析，取位移模式为 代入结点位移后，可解出a1-a6，再代入上式，得 xr, yz 第18页，共34页，编辑于2022年，星期二 其中形函数： ； 单元中位移 2、单元中应变 根据弹性力学理论，空间轴对称问题的几何方程为 第19页，共34页，编辑于2022年，星期二 将u,w表达式代入上式，整理后 第20页，共34页，编辑于2022年，星期二 式中 其中 [B]矩阵中含有变量r，z，因此它不是常数矩阵，即轴对称问题的三角形环形单元不是常应变单元。 第21页，共34页，编辑于2022年，星期二 3、单元中应力 根据弹性力学理论，空间轴对称问题的应力-应变关系为 弹性矩阵: 第22页，共34页，编辑于2022年，星期二 单元中任意一点的应力： 4、单元刚度矩阵 由于被积函数与θ无关，故在三角形截面的环单元的积分可简化为在三角形截面上的积分。故有： 第23页，共34页，编辑于2022年，星期二 单元刚度矩阵的积分参照图示分区，按下式采用数值积分的方法进行 第24页，共34页，编辑于2022年，星期二