有限元与有限差分法基础.ppt

单元平衡方程组装过程 3 ② U6 1’ 1 U2 U1 2’ u1 u2 U5 u3 u4 单元平衡方程组装过程 再组装单元 ② 总体刚度方程 K 称为总体刚度矩阵 U 称为位移向量 F 称为载荷向量 总体刚度矩阵K的特性 对称性 奇异性 稀疏性 非零元素带状分布 约束处理过程 为什么要约束处理 ? 总体平衡方程组是奇异的 消除无限制的刚体运动 使总体平衡方程组存在唯一一组解 约束处理过程—边界条件 边界条件分类 力(载荷)边界条件 位移边界条件 集中载荷力 表面分布力 自重力 热交换引起的温度载荷 固定位移约束 强制位移约束 关联位移约束 约束处理过程—模型简化 x y 约束处理过程—模型简化 y x x y 约束处理过程—约束方程 1 2 3 4 5 6 7 8 9 10 11 12 y x 约束处理过程—约束处理方法 位移约束处理方法 赋0赋1法 乘大数法 约束处理过程—赋0赋1法 强制位移约束条件处理 U4=C 约束处理过程—赋0赋1法 强制位移约束条件处理 U4=C 约束处理过程—赋0赋1法 有6个方程，5个未知数，如果约束方程可以消除有限元平衡方程组的奇异性，则取任意5个方程联立求解，都会得到方程组的唯一一组解。 系数矩阵由原来的对称的变成了非对称的，这对于大规模有限元方程组求解是十分不利的，采用相同的求解方法，在求解时间和矩阵存贮容量方面都增加了一倍。 约束处理过程—赋0赋1法 为了保证系数矩阵的对称性，去掉方程组第4行 约束处理过程—赋0赋1法 引入强制位移约束方程 U4=C，使方程组求解时直接将自由度U4求出 约束处理过程—赋0赋1法 固定位移约束条件处理 U4=0 约束处理过程—赋0赋1法 基本原理 利用初等变换对求解方程组进行相同的行列变换，既保证方程组解不会改变，又可以保持方程组系数矩阵的对称性。 在进行初等变换时，只要保证对方程组系数矩阵做相同的行列变换，就可以保持方程组系数矩阵的对称性。 约束处理过程—乘大数法 乘大数法 基本原理 利用矩阵的初等变换不改变方程组解的思想。 约束处理过程—乘大数法 强制位移边界条件 约束处理过程—乘大数法 强制约束方程 A是一个大数，是系数矩阵中对角线元素K44的1010倍量级以上 为什么要乘以大数A ？ 放大位移约束方程的优势 约束处理过程—乘大数法 强制位移边界条件 约束处理过程—乘大数法 固定位移边界条件 C = 0 约束后的方程组简化为 单元模型构造—整体坐标系法 二维3节点三角形单元 建立整体坐标系oxy 单元模型构造—整体坐标系法 1. 假设插值多项式 2. 首先，利用节点值求 ?0 、 ? 1 和 ? 2 二维3节点三角形单元 单元模型构造—整体坐标系法 A为单元面积 单元模型构造—整体坐标系法 3. 将 ?0 、 ? 1 和 ? 2 代入插值多项式，按u1、u2、u3合并同类项 单元模型构造—整体坐标系法 4. 同理可得 单元模型构造—整体坐标系法 5. 单元插值多项式为 单元模型构造—整体坐标系法 6. 单元插值多项式写成矩阵形式（常用） 单元模型构造—整体坐标系法 7. 单元插值多项式的另一种矩阵形式（不常用） 单元模型构造—整体坐标系法 4节点四面体单元 单元模型构造—整体坐标系法 1. 假设插值多项式 2. 插值多项式为 单元模型构造—整体坐标系法 （i=1,2,3,4） 循环轮换脚标1、2、3、4，相应可以得到a2，b2 ， c2 ， d2 、 a3 ， b3 ， c3 ， d3 、 a4 ， b4 ， c4 ， d4 单元模型构造—整体坐标系法 3. 单元插值多项式写成矩阵形式（常用） 单元模型构造—整体坐标系法 4. 单元插值多项式另一种矩阵形式（不常用） 单元模型构造—整体坐标系法 从理论上讲，整体坐标系法可以求任意单元的形函数，但计算过程太复杂 只能求一维2节点线单元、二维3节点三角形单元和三维4节点四面体单元3种简单单元的形函数 复杂的或二次以上的单元必须采用局部坐标系法求 位移场 u 是形函数 Ni 的线性组合，因此形函数Ni同样具有插值多项式的特性 单元刚度矩阵—2节点线单元 一维2节点线单元 单元插值关系 单元几何关系 单元本构关系 N=[N1 N2] De=E 单元刚度矩阵—2节点线单元 单元刚度矩阵 A为单元截面积；l为单元长度 矩阵B 单元刚度矩阵—三角形单元 二维3角形单元 单元插值关系 单元刚度矩阵—三角形单元 单元几何关系 单元刚度矩阵—三角形单元 单元本构关系 平面应力问题 单元刚度矩阵—三角形单元 矩阵B 单元刚度矩阵—三角形单元 单元刚度矩阵 h为