有限元教程4.pptVIP

四、平面问题有限元法 根据虚功相等： 根据三角形单元的位移模式，单元内任意点的位移可以通过结点位移表示： 练习1： 作用于三角形重心上的集中力为： 求等效结点力向量? 体积力的等效结点力 练习2： x y g 1 三角形单元受重力作用，求i结点等效结点力 表面力的等效结点力 x y ds 在ij边上，取ds其上分布力的合力当作集中力处理 {ps} 将此力向结点等效 将此力向结点等效 所有面力等效结点力叠加 练习3： x y ij边上作用着均布力 ，求三个结点等效荷载 c ij边上任意一点 s l 同理： 根据形函数性质，在ij边上： 练习4：如果某边上分布力呈三角形分布时，等效结点力如何分配？ 4.1.7 总刚及整体荷载向量的形成 单元平衡方程： ① ② 1 2 3 4 ⑧ ① ② ③ ④ ⑤ ⑥ ⑦ 1 2 3 4 5 6 7 8 9 将单刚划分成2×2的子块，按照下标位置送入总刚。 用单刚元素表示出总刚某一元素！！！ 整体荷载向量 每个单元有6个结点荷载，按照结点编码送入总体荷载向量 4.1.8 刚度矩阵性质 1、对称性 2、奇异性 3、正定性 4、稀疏、带状 解释? i j i j 1、对称性 功的互等 2、奇异性 令： 根据互等定理： 行、列线性相关→奇异 物理意义：单元（整体结构）在保持平衡的条件下可以有任意刚体位移 3、正定性 意义：使结点产生位移所需要的力在该位移上所作的功恒正 4、稀疏、带状 在总刚中，不相关结点在总刚中对应位置元素为0 半带宽一维存贮技术（只保存包括主对角线元素在内的半个带状区域的数值及地点）： 节省内存、提高计算效率 4.1.9 边界条件处理 一、删行删列法 二、乘大数法 三、置1法 4.1.10 方程组求解 直接解法： 迭代法： 高斯消元法； Choleshy 分解法； 三角分解法； 分块法； 波前法 Jacobi迭代法； Gauss-Sidel迭代； 超松弛迭代 * §4.1 平面三角形常应变单元 §4.2 四结点矩形单元 §4.3 平面高次单元 主要内容： §4.1 平面三角形常应变单元 背景： 杆、梁几何形状规则，受力后变形规则，位移场可用简单函数表示； 形状复杂、边界条件复杂时，很难得到应变及应力场的解析解。 三角形常应变单元特点：对于边界条件适应性强；精度低、网格密 4.1.1位移模式选择 4.1.2 形函数及其性质 4.1.3几何矩阵 4.1.4 应力矩阵 4.1.5 单刚 4.1.6 等效结点荷载 4.1.7 总刚及整体荷载向量的形成 4.1.8 刚度矩阵性质 4.1.9 边界条件处理 4.1.10 方程组求解 4.1.1位移模式选择 位移模式－单元内位移的分布函数。 选择位移模式遵循的原则： 4.1.2 形函数及其性质 选择位移模式： x y 边界条件：－结点位移 代入位移模式 待定系数 其中： A的伴随矩阵 ？ 右手准则 注意：单元编号按照逆时针方向进行，否则△为负值？ 记 其中： 全部为节点坐标函数，可计算出每个代数余子式具体数值 根据位移模式： 同理： 形函数矩阵 形函数性质？ 1、本结点为1，其他结点为0 2、单元内任意点（x,y）形函数之和为1 x y ij的直线方程： －A32 A33 3、在对边上，形函数值为0 在ij边上，有： 代入m结点形函数： 在ij边上任意点： ? 在ij边上，只要结点位移确定了，则边界上位移具有唯一性。 不同单元之间即不会开裂，也不会重叠－－－协调单元 4、形函数与面积坐标 单元内任意点（x,y）形函数之和为1 4.1.3几何矩阵 4.1.4 应力矩阵 平面应力问题本构方程： 写为矩阵形式 弹性矩阵 平面应变问题本构方程： 弹性矩阵 应力矩阵 4.1.5 单刚 最小势能原理：在外力作用下，弹性体产生变形。变形不是任意产生的，应使整个系统的总势能最低。 三角形常应变单元中，[B]、[D]中元素全部为常数 4.1.6 等效结点荷载 有限元基本思路，求解方程组得到基本未知数δ 结点力等效原则→静力等效（虚功相等） 1、结点集中力； 2、单元内集中力的等效； 3、体积力的等效结点力； 4、表面力的等效节点力。 单元内集中力的移置 设在单元内集中力作用下，M点产生虚位移 ；相应单元结点产生虚位移为： *