有限元法剪切锁闭现象的研究论文答辩.ppt

可以看出当梁为细长梁时,希望能得到 。如果 不恒为零,则一定要求矩阵 为奇异阵,使得任何情况下有 。另一种办法在计算 的积分时，不采用精确积分，而用一点积分（缩减积分）来计算， * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 4.1.2 选择性缩减积分 实质上同缩减积分没有区别，只是在处理实际问题时加以改善了，对于可能会发生锁闭的局部采用缩减积分，减少计算量。 4.1.3 非协调元 非协调元是为克服完全积分的，一阶线性单元的剪切闭锁问题的一种尝试。由于剪切闭锁是单元的位移场不能模拟与弯曲相关的变形而引起的，所以在一阶单元中引入一个增强单元变形梯度的附加自由度。这种对变形梯度的增强可以允许变形梯度在一阶单元的单元域上有一个线性变化，如图4.4(a) 所示，而标准的单元数学公式在单元中只能得到一个常数变形梯度，如图4.4（b）所示，这导致与剪力闭锁相关的非零剪切应力。 在弯曲问题中，非协调元可能产生与二次单元相当精度的结果，但计算成本却明显的下降。然而，它们对单元的扭曲很是敏感。 图4.4 变形梯度的变化 图4.5使用位移场增强而非变形梯度增强非协调模式单元间潜在的运动非协调性，ABAQUS/Standard使用后一表达式为其非协调模式单元 图4.6 非协调模式单元的扭曲网格 图4.7 平行和交错扭曲对非协调元的影响 4.2 厚板单元解析试函数法 在有限元的离散法中吸取解析法的优点并加以利用，很早就有文献论述。利用基本解析解作为试函数而构造出的广义协调元 在厚板平衡微分方程中，如果将内力改用位移表示，就得到厚板位移法的基本微分方程。考虑载荷为零时的齐次问题，则有： (4.8) 设F(x,y)是双调和函数，满足双调和方程： (4.9) 其中： (4.10) 则下列为上式微分方程的解： (4.11) 设 是一个特征长度。引入无量纲坐标 ， ： (4.12) (4.13) 代入得到: (4.14) 如果 ,则厚板理论的解析解退化为薄板理论的解析解： (4.15) (4.16) 结论 在有限元分析中，不能忽视切应力的影响，需要对具体问题进行必要的分