理学院毕业论文答辩范例.ppt

有限元法剪切锁闭现象的研究 专 业：工程力学 答辩人： XXX 导 师： XXX 武汉科技大学理学院 2019 年 6 月 16 日 论文题目： 2 答辩提纲 ? 课题背景及任务 ? 有限元分析的误差 ? 有限单元特征及剪切闭锁产生 ? 剪切闭锁的对策 ? 结论及下一步工作 ? 致谢 3 课题背景及任务 当梁不是细长梁时，梁变形后的横截面垂直于中面的假设不再成 （ Kirchhoff 假设）立，需考虑剪切变形的影响，于是 Timoshenko 提 出了剪切变形的梁理论，如果对扰度函数与转角函数进行独立插值， 且考虑剪切变形的影响，构造出的单元为 Timoshenko 梁单元 （ Timoshenko beam element ）。此单元在退化为细长梁单元时，会 导致剪切闭锁现象。经典薄板理论以 Kirchhof 假设为基础，忽略了横 向剪切的影响，故计算结果与实验值相比，总是低估了挠度，高估了 自然频率。由于横向剪切变形对厚板及复合材料的板单元影响较大， 故忽略剪切变形的薄板理论已不再适用，于是 Mindlin 等人提出了板的 剪切变形理论 , 以此构造出的为 Mindlin 板单元（ Mindlin plate element ）。与 Timoshenko 梁单元类似，当板很薄时，会发生剪切闭 锁与零能模式。对于块体，由于形函数的阶数过低，而单元受到复杂 应力作用，如产生弯曲扭转等，导致形函数不能真实表达。本文讨论 剪切闭锁产生的原因，与对策。 4 有限元分析的误差 ? 结构体由于本身存在有自然的连接关系即自然节点 , 所以它们的离散化 叫 自然离散 ? 人为的在连续体内部与边界上划分节点 , 以单元连续的形式来逼近原来 复杂的几何形状 , 这种过程叫 逼近性离散 (approximated discretization) ? 控制误差的 h 方法（ h-version ,h-method ）和 p 方法 (p-version ,p- method) 。 5 2.1 求解精度的估计 ? 在这里我们考察平面单元的求解精度与收敛速度，单元的位移场 可 以展开为以下级数： (2.1) 对于满足完备性（ completeness ）和协调性（ compatibility ）要求的 协调元，当 时，有限元分析的结果是单调收敛的如图 (2.1) 图 2.1 有限元分析结果的收敛情况（纵坐标是总势能） ( ) ( ) i i i u u u u x y x y ? ? ? ? ? ? ? ? ? ? K 6 ? 还可以就两次网格划分所计算的结果进行外推以估计结果的准确值如 第一次网格划分的结果是 , 然后进一步将各单元尺寸减半进行网格划 分 , 得到结果为 . 假设该单元的收敛速度是 , 则其准确解可以按如下方 法估计 : ? (2.2) ? 具体对平面 3 节点三角形单元有 s=2 ，上式可化为： （ 2.3 ） ? 可由此估计出准确解： (2.4) ? 由有限个单元的试函数来逼近整体域的场函数所引起的误差，即离散 误差。 ? (1) (2) ( ) (( / 2) ) s i i s i i u u O h u u O h ? ? ? (1) 2 (2) 2 ( ) 4 (( / 2) ) i i i i u u O h u u O h ? ? ? ? (2) 1 (4 ) 3 i i i u u u ? ? 7 2.2 提高精度的ｈ方法和ｐ方法 ? （１） ｈ方法 （ h-version ） : 不改变各单元 上的形状函数，只通过渐渐加密有限元的 网格使结果趋于准确解 ? （２） ｐ方法 （ p-version ） : 同ｈ方法相反 ，它是不改变单元网格，而采用较高阶的 多项式来进行插值。 8 3. 有限元的单元特征及剪切闭锁的产生 ? Timoshenko 梁单元的剪切闭锁现象 ? Mindlin 板单元的剪切闭锁现象 ? 块体单元的闭锁现象 9 3.1 Timoshenko 梁单元 当梁不是细长梁时，梁变形后的横截面垂直于中性层的假设不再成立（ Kirchhoff 假设）这时需要考虑梁的剪切变形，考虑剪切变形的几何描 述如图 3.1 所示。 图 3.1 具有剪切变形影响的梁变形 变形前的截面 中性层 垂直于中性 层的截面 变形后的截 面 变形前的截面 中性层 β γ d ω /dx (a) 变形前 (b) 变形前 10 ? 设剪切变形为 ，则 （ 3.1 ） ? 设单元的总挠度函数 和截面转角函数 的单元插值模式为（线性插值 ）： （ 3.2 ） （ 3.3 ） 其中 ， 分别为单元节点 1 和 2 的挠则 （ 3.4 ） 由（ 3.1 ）可得到剪应变为： （ 3.5 ） d dx ? ? ? ? ? 1 1 2 2