[理学]弹性力学 第6章_平面问题的有限单元法.ppt

§6.7 结构的整体分析 结点平衡方程组 是单元结点的局部编号； 是整体结点的整体编号。 代入式 ，可表示为 将式 按整体结点编号排列，得整个结构的平衡方程组。 §6.7 结构的整体分析 结点平衡方程组 考虑结构的约束条件后，从式 求出 ，就可以求出各单元的位移和应力。 －－整体结点位移列阵， －－整体结点荷载列阵， －－整体劲度矩阵。 §6.7 结构的整体分析 结点平衡方程组 例1 列出图示结构i 结点的平衡条件。 §6.7 结构的整体分析 结点平衡方程组 结点i联系四个单元①～④，对i点列方程： §6.7 结构的整体分析 结点平衡方程组 整体劲度矩阵是由单元劲度矩阵的元素集合而成，称总刚的组装。 组装原则：对号入座 整体劲度矩阵的特点： （1）对称性；（2）奇异性；（3）稀疏性；（4）非零元素成带状分布 §6.7 结构的整体分析 结点平衡方程组 例2 （见书中P.121） 设有对角受压的正方形薄板，荷载沿厚度均匀分布，为2N/m，试用有限单元法求解薄板内的应力。 §6.7 结构的整体分析 结点平衡方程组 解： （1）建立有限单元法的计算模型： 考虑对称性： 4个单元，6个结点 平面应力问题 §6.7 结构的整体分析 结点平衡方程组 （2）列出整体结构的结点平衡方程组K?=FL §6.7 结构的整体分析 结点平衡方程组 （3）建立结构的整体结点荷载列阵FL 整体结点位移列阵： §6.7 结构的整体分析 结点平衡方程组 （4）求出各单元的劲度矩阵k §6.7 结构的整体分析 结点平衡方程组 应用式（6-38），可得两种单元的劲度矩阵都是： §6.7 结构的整体分析 结点平衡方程组 （5）求出整体劲度矩阵K 建立结点的整体编码与单元内的局部编码之间的关系： §6.7 结构的整体分析 结点平衡方程组 §6.7 结构的整体分析 结点平衡方程组 整体劲度矩阵的组装： §6.7 结构的整体分析 结点平衡方程组 §6.7 结构的整体分析 结点平衡方程组 （6）结点平衡方程组（a）的简化与求解 由于有位移边界条件u1=u2=u4=v4=v5=v6=0，未知的整体结点位移列阵简化为： 整体荷载列阵简化为： §6.7 结构的整体分析 结点平衡方程组 整体平衡方程： §6.7 结构的整体分析 结点平衡方程组 求解后，得结点位移： §6.7 结构的整体分析 结点平衡方程组 （7）由位移求应力： 对于单元Ⅰ，Ⅱ，Ⅳ： 对于单元Ⅲ： §6.7 结构的整体分析 结点平衡方程组 各单元的应力如下： §6.7 结构的整体分析 结点平衡方程组 §6.3 单元的位移模式与解答的收敛性 因为当单元 时，单元中的位移和应变都趋近于基本量－－刚体位移和常量位移。 （1）位移模式必须能反映单元的刚体位移。 （2）位移模式必须能反映单元的常量应变。 §6.3 单元的位移模式与解答的收敛性 可见刚体位移项在式（a）中均已反映。 与刚体位移相比， 将式（a）写成 §6.3 单元的位移模式与解答的收敛性 （3）位移模式应尽可能反映位移的连续性。 即应尽可能反映原连续体的位移连续 性。在三角形单元内部，位移为连续；在两单元边界ij 上， 之间均为线性变化，也为连续。 对式（a）求应变，得： 可见常量应变也已反映。 §6.3 单元的位移模式与解答的收敛性 （1）和（2）是必要条件，而加上（3）就为充分条件。 为了保证FEM的收敛性： §6.3 单元的位移模式与解答的收敛性 思考题 1.应用泰勒级数公式来选取位移模式，为什么必须从低次项开始选取？ 2.试考虑：将结构力学解法引入到求解连续体的问题时，位移模式的建立是一个关键性工作，它使得单元(连续体)内部的分析工作都有可能进行了。 §6.4 单元的应变列阵和应力列阵 其中， 单元中的位移函数用位移模式表示为 §6.4 单元的应变列阵和应力列阵 应用几何方程，求出单元的应变列阵： §6.4 单元的应变列阵和应力列阵 再应用物理方程，求出单元的应力列阵： B ——称为应变矩阵，用分块矩阵表示， S——称为应力转换矩阵，写成分块形式为 §6.4 单元的应变列阵和应力列阵 对于线性位移模式，求导后得到的应变和应力，均成为常量，因此，称为常应变（应力）单元。应变和应力的误差量级是 其精度比位移低一阶，且相邻单元的应力是跳跃式的。 §6.4 单元的应变列阵和应力列阵 思考题 1.如果在位移模式中取到泰勒级数中的二次幂项，略去 高阶小量，试考虑位移、应变和应