有限元基础讲稿第2讲之4.pptVIP

? Copyright 2000 PTC. Visit Visit ? Copyright 2000 Parametric Technology Corporation Thanks for your listening * ? Copyright 2000 PTC. Visit 第四节 2.1 有限元基本方法 2.2 平面问题 2.3 轴对称问题 2.4 空间问题 2.5 等参元和板壳元 四、空间问题 有 限 元 分 析 空间问题比轴对称问题要复杂得多，本节先讨论其中一种最简单的形式，即用四面体常应变（也是常应力）单元处理的方法。 典型的四面体单元如右图所示。设其结点编号依次为i、j、m、p，空间各结点有三个位移分量，设为u、v、w，写成列阵形式为： 四、空间问题 有 限 元 分 析 整个单元的结点位移列阵为： 设其位移模式为简单线性函数： （xi yi zi）、（xj yj zj）、（xm ym zm）和（xp yp zp） 代入 四、空间问题 有 限 元 分 析 得到关于 的方程组，用矩阵表示为： 与平面三角形常应变单元类似，求得系数 后可得形函数： 四、空间问题 有 限 元 分 析 式中： aj、bj、cj、dj、am bm、cm、dm、ap、bp、cp、dp类似 四、空间问题 有 限 元 分 析 V是四面体ijmp的体积。为了不使体积出现负值，单元的4个顶点的标号i、j、m、p必须按一定的规则排序，在右手坐标系下，要使右手在按照i－j－m的旋向转动时是向p方向前进。 求出形函数后，位移可表示为： 四、空间问题 有 限 元 分 析 矩阵形式 对所有单元都进行如上式的线性插值，不难证明两个单元的交界面上的位移是连续的。 四、空间问题 有 限 元 分 析 单元体内任一点的位移知道后，就可以确定单元体内任一点的应变，由几何方程： 求得应变与位移的关系为： Bj、Bm、Bp 子块依次类推 四、空间问题 有 限 元 分 析 由物理方程 ，在各向同性的线弹性范围内有： 四、空间问题 有 限 元 分 析 代 入 应力矩阵 可见应力矩阵[S]由4个6×3的子块组成，具体形式为： 四、空间问题 有 限 元 分 析 与平面问题中的推导方法类似，空间问题的单元刚度矩阵为： 单元刚度矩阵 仍可以分块表示为： 四、空间问题 有 限 元 分 析 对应子块均为3×3方阵，各子块按下式计算： （s、r=i,j,m,p） 单元刚度矩阵是由单元结点的坐标和单元材料的弹性常数决定的，它是一个常数矩阵。 如果弹性体划分为m个单元和n个结点，经过类似于平面问题的组合，得整个结构平衡方程： 四、空间问题 有 限 元 分 析 总刚度矩阵 可见，[K]中每个元素都表示一个3×3的子块，它是各单元内对应子块之和，即： （r＝1,2,…,n；s＝1,2,…,n） 与平面问题相同，总体刚度矩阵也是对称的、带状稀疏的，在消除刚体位移（加必要的约束）后，它是正定的。 四、空间问题 有 限 元 分 析 3×3子块[Krs]是若干单元中的之和，其中[krs]的个数是m，但这m个[krs]中只有第r号结点和第s号结点落在同一单元内时对应的[krs]不为零，r号结点与s号结点不在同一单元时[krs]都是零这就是总体刚度矩阵是稀疏矩阵的原因。 空间问题中等效载荷列阵的计算比平面问题更为复杂。点线、面、体四种不同形式的载荷，只有作用在结点上的集中载荷直接列入载荷列阵中以外，其余三种均要进行等效移置。 四、空间问题 有 限 元 分 析 线力：即作用在某边上的分布力等效移置到该边所联系的两端结点上，方法与平面问题所介绍的方法相同 面力：当作用在单元面积上的载荷为 ，则结点i上的力为： 式中，Ni是形函数在i点的值，A为载荷作用面积。如果这种载荷是在三角形ijm面上呈线性分布，则作用在i点的力表示为： 四、空间问题 有 限 元 分 析 式中，pi、pj、pm分别为载荷分布在i、j、m点处的强度，Δ为三角形面积。 Qj、Qm 依次类推 体积力：当单位体积上的载荷为 ，则结点i处的力记为： 匀质单元分配到四个结点上，其值为 ，γ为密度或重