平面问题有限元解法(公式推导讲解)精编版.pptVIP

* 平面有限元解法——整体劲度矩阵 每个单元，节点的局部编码和整体编码对应关系已经确定，每个单元劲度矩阵中任一子矩阵在整体劲度矩阵中的位置及其力学意义也就明确了。如单元Ⅰ的kii，即k33，它的四个元素表示当结构的节点3沿x或y方向有单位位移时，在节点3的x方向或y方向引起的节点力。 暂时不考虑位移边界条件，把所分析结构的整体节点平衡方程组列出： 整体劲度矩阵写成6×6的矩阵，它的每个子块是2×2的矩阵，实际它是一个12×12的矩阵。如K23，它的四个元素表示当结构的节点3沿x或y方向有单位位移时，在节点2的x方向或y方向引起的节点力。 * 平面有限元解法——整体劲度矩阵续 由于节点3和节点2在结构中是通过Ⅰ和Ⅲ这两个单元相联系，因而K23应是单元Ⅰ 的k23和单元 Ⅲ的k23之和。同理，可以找到各单元劲度矩阵中所有子矩阵在整体劲度矩阵K中的位置，得到整体劲度矩阵。 式中k的上标1,2,3,4表示是哪一个单元的劲度矩阵中的子矩阵，空白处是2×2的零矩阵。 * 平面有限元解法——整体劲度矩阵续 对于单元Ⅰ、Ⅱ、Ⅳ，根据公式，可求得A=0.5m2， 将上式中各子块的具体数值代入整体刚度矩阵K表达式中，得出整体刚度矩阵。 对于单元Ⅲ，根据公式，可求得A=0.5m2， 把μ＝0，t＝1m，代入单元的劲度矩阵，得两种单元的劲度矩阵k都是： (37’) * 平面有限元解法——整体劲度矩阵续 整体刚度矩阵K (38’) * 平面有限元解法——位移边界条件 位移边界条件为： (39’) 因此，整体结点的位移列阵就简化为： 与这6个零位移分量相应的6个平衡方程不必建立，因此，将整体刚度矩阵中，第1、3、7、8、10、12各行以及同序号的各行划去，因而整体劲度矩阵K简化为： * 平面有限元解法——整体节点载荷列阵 确定了每个单元的节点载荷列阵： (38) 根据各单元的节点局部编码与整体编码的关系，确定三个子块FLi,FLj,FLm在FL中的位置。 由于该结构只是在节点1受有向下1N/m的载荷，因而，非零元素子块，只有 在考虑了边界条件后，整体载荷列阵为： * 平面有限元解法——求解整体节点载荷列阵 求解化简后的整体刚度矩阵： (39) 求解以后，得节点位移： * 平面有限元解法——求解应力转换矩阵 应用单元的应力转换矩阵S，求出各单元中的应力： 根据μ＝0，以及已求出的A、b和c的值，再由式(21)和(22)得出应力转换矩阵如下，对于单元Ⅰ、Ⅱ、Ⅳ ： 对于单元Ⅲ * 平面有限元解法——求解各单元中的应力（续） 应用单元的应力转换矩阵S，求出各单元中的应力： Pa 单元Ⅰ 单元Ⅱ Pa * 平面有限元解法——求解各单元中的应力 应用单元的应力转换矩阵S，求出各单元Ⅲ、Ⅳ中的应力： Pa 单元Ⅲ 单元Ⅳ Pa * 平面有限元解法——计算机编程解题步骤 划分单元格，并按照一定的规律将所有的结点和单元格分别编上号码，需要注意单元编码和整体编码的对应关系。 选定一个直角坐标系。按照计算程序的要求，填写各种输入信 息有：每个结点的坐标值，即x1,y1,x2,y2 …等；材料的弹性常数值；各种载荷信息，即载荷点的点号及载荷的大小等；约束信息，即哪些结点哪个方向上的位移为零或为某个已知值。将这些信息按照计算机程序规定的格式输入。 计算程序中对输入的各种信息进行加工、运算，一般均有如下几步：输入初始数据，形成整体刚度矩阵K；形成整体载荷列阵FL；求解线性代数方程组，解得结构的整体结点位移阵列δ；计算各单元的应力分量及主应力、主向；打印计算结果。 * 平面有限元解法——计算机程序界面 * 平面有限元解法——计算机程序计算结果 * 通用有限元计算程序ANSYS计算结果 * 通用有限元计算程序ANSYS计算结果 * THE END * * 圣维南原理的应用 例，厚度δ=1的梁中，左右两端x=±l，的边界面是正、负x面，其上作用有一般分布的面力± 。按照严格的应力边界条件，应力分量在边界上满足： 上式要求在边界上y值不同的各点，应力分量与对应的面力分量必须处处相等，这种严格的条件是较难满足的。 当lh时， x=±l 是梁的边界的一小部分，可以应用圣维南原理，利用静力等效条件来代替，即，使应力的主矢量和主矩分别等于对应的面力的主矢量和主矩。 * 圣维南原理的应用 应力的主矢量和主矩的绝对值分别等于面力的主矢量和主矩的绝对值； 面力的主矢量和主矩的方向就是应力的主矢量和主矩的方向。 * 按位移法求解平面问题 以上几节已经建立了弹性力学平面问题的 基本方程和边界条件，即：平衡微分方程、几何方程和物理方程，以及位移的边界条件和应力的边界条件。 求解弹性力学平面问题即求解3个应力分量、3个形变分量及2个位移分量的未知函数。通