第二章 弹性力学平面问题的有限元法.pptVIP

第二章  弹性力学平面问题的有限元法 空间离散：建立坐标系，选择单元类型，将连续域划分为有限个单元；形成单元总数、结点总数、单元号、结点号、结点坐标、单元内结点次序（单元结点信息）;在结点上加上正确的位移约束。 图示：单元（1，2，3），已知E,h,μ=0，平面应力，坐标（0，1）（0，0）（1，0），写出单刚；若结点次序变化，讨论单刚的变化。 图示：单元（1，2，3），已知E,h,μ=0，平面应力，坐标（0，1）（0，0）（1，0），写出单刚；若结点次序变化，讨论单刚的变化。 1. 由结点平衡法建立有限元方程 总刚度矩阵 3. 利用最小势能原理建立有限元方程 图示：有限元模型，已知E,t,μ=0，平面应力，用刚度集成法求总刚度矩阵K 图示：有限元模型，已知E,t,μ=0，平面应力，用刚度集成法求总刚度矩阵K 图示：有限元模型，已知E,t,μ=0，平面应力，用刚度集成法求总刚度矩阵K K的存贮方式有等带宽存贮，变带宽一维存贮等 图示：有限元模型，已知E,t,μ=0，平面应力，用刚度集成法求总刚度矩阵K 图示：有限元模型，已知E,t,μ=0，平面应力，用刚度集成法求总刚度矩阵K 图示：有限元模型，已知E,t,μ=0，平面应力，用刚度集成法求总刚度矩阵K 图示：已知厚度为t的边长为2 有限元模型，已知E, μ=0，平面应力，求各单元应力及有关结点处的约束反力。 2.2 连续体弹性力学分析有限元法的一般过程 一、连续域的离散化 二、单元特性分析 三、离散单元的装配和集成 四、引入位移约束 五、解方程，求结点位移 六、求单元应力、应变、约束力等 连续体弹性力学分析有限元法的一般过程 连续体弹性力学分析有限元法的一般过程 结构离散化 单元分析 整体分析 位移约束条件的引入 方程求解 有关量计算 结束 2.3 高阶三角形单元 一、 面积坐标 中有任意一点P P点的面积坐标 1.面积坐标定义 x y O i j m P 2. 面积坐标的性质： ①一点的三个坐标值是 线性相关的 ② ③在三个结点上 当Li为常数时，表示 内平行 于jm边的线段 x y O i j m P 一、 面积坐标 3、面积坐标与 整体坐标的关系 x y O i j m P 一、 面积坐标 (i,j,m) 4、面积坐标与微积分运算 一、 面积坐标 二、 高次三角形单元 1. 6结点二次三角形单元 x y O 1 3 4 2 5 6 6结点：三个角结点，三个边中点 12个自由度 位移函数 —形函数矩阵 —形函数 二、 高次三角形单元 1. 6结点二次三角形单元 x y O 1 3 4 2 5 6 P为插值函数的次数 —为通过除结点i外的其余结点的直线方程的左端。 —为 在结点i的取值。 x y O 1 3 4 2 5 6 二、 高次三角形单元 1. 6结点二次三角形单元 —为通过除结点i外的其余结点的直线方程的左端。 —为 在结点i的取值。 x y O 1 3 4 2 5 6 1. 6结点二次三角形单元 1. 6结点二次三角形单元 x y O 1 3 4 2 5 6 位移函数 1. 在单元边界上，位移是按抛物线变化的，每条公共边界上 有3个公共结点，可保证相邻两单元公共边上的位移连续性。 是否收敛? 2. 位移函数里包含了常数项和坐标一次项 二、 高次三角形单元 2. 10结点三次三角形单元 10结点：三个角结点，6个三分点， 一个形心。 20个自由度 x y O 1 3 4 2 5 6 7 8 9 10 位移函数 —形函数矩阵 —形函数 2. 10结点三次三角形单元 x y O 1 3 4 2 5 6 7 8 9 10 角结点 三分点 形心 2.4 四结点矩形单元 单元中心： 引入单元局部坐标： 2.4 四结点矩形单元 矩形单元：点的编号1，2，3，4为单元内点的次序，逆时针（右手）转向。结点坐标分别为（xi,yi）;i=1,…,4 x y x′ y′ 1 2 3 4 则结点坐标： 引入单元局部坐标： x y x′ y′ 1 2 3 4 一、单元位移场 2.4 四结点矩形单元 2.4 四结点矩形单元 是行阵， 是列阵。 例： 解： 单元② ② ① 解： 例： 的性质: 一、K是对称矩阵。 二、K是稀疏矩阵。 三、K的非零元素呈带状分布。 四、K是奇异矩阵。 的存储: 行半带宽 半带宽D D=(所有单元中的最大点号差+1)×2 （平面） D=8 的存储: 原总刚度矩阵K 半带宽刚度矩阵T (2)分析 O x y i j m e i j m （i,j,m） i j m 点号对应