限元(2b).pptVIP

第二章 平面问题有限元法续 ? 结构离散 ? 位移函数 ? 单元分析 ? 总刚集成 ? 载荷移置 ? 约束处理 ? 求解线性方程组 ? 计算其它物理量 ? 计算结果处理 ?结果显示、打印与分析 结构离散??单元分析??整体分析?? 载荷处理??约束处理 2.6 约束处理（1） 必要性:已知[K]和{R}, {q}仍然不可解，因为: （1）数学上：IKI=0, [K] 为奇异矩阵，[K] -1不存在 , （2）物理上：无约束结构可以有刚体运动，位移不确定； 处理方法： （1）置一置零法； （2）乘大数法。 得到总刚方程 问题；仍无法求解。必须进行约束处理。 2.6 约束处理（续1） 置一置零法（处理边界位移为0的情况，这是一种精确解） 1.将[K]矩阵中零位移分量所对应行和列的主对角元素置为1，而其它元素皆为0； 2.在{R}列阵中，将0位移分量所对应的节点载荷也变为0； 2.6 约束处理（续2 ） （1）刚性约束边界条件： 2.6 约束处理（续3 ） （2）强迫位移约束边界条件 强迫位移定义：边界位移不为零。 处理方法:乘大数法（近似解法） 若节点位移是已知值，则将该位移分量所对应的[K]阵中的主对角元素置为大数，再将载荷列阵{R}中对应的分量置为大数乘以已知的节点位移值，其余不变。 2.6 约束处理（续4） 2.7、求解线性方程组 线性方程组的解法分为两大类：直接解法，迭代解法。 在方程组的阶数不是特别高时，通常采用直接解法。当方程组的阶数过高时，为避免舍入误差和消元时有效数损失等对计算精度的影响，可以选择迭代方法。 直接解法以高斯消去法为基础，包括高斯消去法、等带宽高斯消去法、三角分解法，以及适用于大型方程组求解的分块算法和波前法等。 迭代算法有高斯-赛德尔迭代、超松弛迭代和共轭梯度法等。 2.8、计算其它物理量 应力矩阵 应变矩阵 弹性矩阵 2.9、应力应变计算结果处理 2.10、后处理：结果显示、打印、分析 应力： （1）米塞斯应力； （2）主应力； 后处理： 应力排序，云图，等值线； 变形； 变形图，变形动画。 简例 1.结构离散（6节点，4单元） 2.逐元计算单刚[k]e；单元1，2，4同。 3.组装形成[K] 4.形成{P}T={Px1 P1 || Px2 0 ||0 0|| Px4 Py4|| 0 Py5 || 0 Py6} 简例（续1） 5.约束处理：采用置一法 约束条件：U1＝U2＝U4＝V4＝V5＝V6＝0 约束行列：1，3，7，8，10，12 A、对[K]处理:第1行 K11=1，Ki1=K1i=0（i1) 第3行K33=1,其余的Ki3和K3i＝0 ……… ……… 简例（续2） B、对{P}处理：第1，3，7，8，10，12元素归0。 6.解方程，得V1，V2，U3，V3，U5，U6 最终求得各节点位移 7.逐元求应力 A、逐元求 ，从 对号取出 B、逐元求[D]e和[B]e C、逐元求应力,