有限元分析 (FEA)课件.pptVIP

Instructors Guide Introduction to ANSYS 5.3 Chapter II: Finite Element Analysis (FEA), Lesson 1: The Finite Element Analysis (FEA) Method 第四章 轴对称与空间问题 4.1 轴对称问题 三角形截面环状单元刚度矩阵 四、单元的刚度矩阵 应用实例 4.2 空间问题 4.2.1 四面体单元的刚度矩阵 设定位移函数 4.3 空间等参元与空间轴对称等参元 *Introduction to ANSYS - Release 5.5 (001128) 一、柱面坐标 A（r, z, ? ) 二、应变分量和位移分量的关系 对于轴对称问题，构件上r,z值相同的一个圆周上各点的位移、应变和应力都与?无关，只是r,z的函数。 对于应变，在r,z子午面内，有 对于轴对称问题，径向位移u还要引起环向应变，有 在轴对称问题中，应变有四个分量 三、应力和应变的关系 （4-1） （4-2） [D]为轴对称问题中的弹性矩阵。 （4-3） 一、设定位移函数 依照平面问题的三角形板单元，采用线性位移函数，即设。 依照平面问题的作法，将各节点的坐标值和节点位移代入上式，可得 （4-5） （4-4） 其中 A是环形单元在子午面rz上的截面积。 二、由位移函数求单元的应变 将位移函数代入位移与应变关系方程（式4-1 ）中，可得 三、单元的应力矩阵 （4-6） 由刚度矩阵的普遍公式 对三角形截面环状单元 代入普遍公式，可得 可以把单元刚度矩阵写成下列分块形式 （4-10） 其中的子矩阵为 在实用上，为了简化计算和消除对称轴上由于r=0所引起的麻烦,常把各个单元中的r及z近似地当作常量，并且分别等于各单元形心的坐标，即 此时 经过这样简化后，就可把各个单元近似地当作常应变单元 （4-11） （4.11）式成为 写成显式的形式为 只要网格划分不是很稀，这样的近似计算引起的误差不大，而计算工作却大大简化 图4.6所示为一无限长的受内压的轴对称圆筒，该圆筒置于内径为120mm的刚性圆孔中，试求圆筒内径处的位移。结构的材料参数为：E=200GPa 结构的结点位移列阵 结构的结点外载列阵 建立整体刚度方程 边界条件的处理及刚度方程求解 与三节点三角形板单元相类似，采用线性位移函数，即设。 写成矩阵形式，则为 （4-25） 将4个节点坐标分别代入式6.1-1中的第一式，得 （6.1-2） 由此解出 把所求得的{a}代入u（x,y,z）的表达式，得 式中 用同样方法，可以得到其余两个位移分量 令 称为四面体单元的形状函数 写成矩阵形式为 二、应变矩阵 （4-35） 将（4.35）式代入几何方程 ,得到单元体内的应变 式中 （4.37） 由物理方程 ,得到应力列阵 其中 上式表明[B]中的元素都是常量，因此单元中的应变也必是常量。因此采用线性位移模式的四面体单元称为常应变四面体单元。 （4.38） 4.2.2 单元刚度矩阵 [S]称为应力矩阵。显然，在每个单元中的应力也是常量 写成分块形式 （4.43） 其中子矩阵 其中子矩阵 单元刚度矩阵是由单元结点的坐标和单元材料的弹性常数所决定，它是一个常数矩阵 六面体单元 典型的空间等参元为空间8～20节点六面体等参元 母单元与子单元的坐标变换式和位移模式可统一写成如下形式 （4.54） （4.55） 式中n为单元结点数。当 n=8时，指的是8结点等参单元，它的形函数 （i=1，2，…，8） 其中 当 n=20时，指的是20结点等参单元，对于单元的20个结点，仿照上述方法，可以分别写出20个结点的形函数 按几何关系式 * * * *