应力函数与其对偶理论在有限元中的应用.pdfVIP

第36卷第4期 力 学 学 报 Vbl．36．NO．4 2004年7月 ACTAMECHANICASINICA July，2004 应力函数及其对偶理论在有限元中的应用D 黄若煜吴长春 (上海交通大学建工学院，上海200240) 摘要借助于Cosserat连续介质模型，探讨了应力函数和位移对避免有限元C1连续性困难的互补性作用．通 过对应力函数对偶理论的深入分析，为将应力函数列式得到的余能单元转化为具有一般位移自由度的势能单元 提供了严格的理论基础，在此基础上，给出应用应力函数构造有限元的一般方法． 关键词有限元， C1连续性， Cosserat模型，应力函数，对偶理论 引 言 指出了板弯曲与平面弹性的比拟关系；Yu[引、钟万 勰等【4】指出了平面偶应力与R／M板弯曲的比拟关 在一些力学问题中，由于应变是位移的二阶导 系．在这些比拟关系中，互为比拟的双方在数学上 数，变分法要求位移场应当具有至少一阶导数的连 是同构的，而且任_方的应力函数与另一方的位移 续性，这就构成了有限元中众所周知的C1连续性困 互为比拟，比如，板弯曲中的挠度W与平面弹性的 难．薄板弯曲和含偶应力的弹性理论都属于具有这 Airy应力函数比拟，而平面弹性的位移与板弯曲的 种本质性困难的问题． 弯矩(应力)函数比拟【引，等等．因此，应用应力函 总结板弯曲和偶应力理论的有限元现状，可以 数，不仅可以避免C1连续性困难，还可以借用与之 把解决C1连续性困难的方法大致归结为如下3种： 比拟的相应理论体系的现有有限元方法，这对构造 1)放松很难实现的强连续性条件，代之以满足 板弯曲和偶应力理论的有限元无疑是很有意义的． 弱连续条件的非协调元方法．拟协调元、广义协调 余下的困难是：以应力函数为节点自由度的单元不 元属于这类单元． 便于实际应用，应当将其转化为具有一般位移自由 2)基于多变量变分原理的杂交／混合元方法． 度的单元，文献[5—8】对这类问题进行了深入讨论， 3)将当前问题归结为不存在C1连续性困难的 提出了一套系统化的方法，并且以大量数值实践表 一个更广泛问题的特例，在后者的理论框架中进行 明了方法的有效性． 有限元列式，再施加相应的约束获得前者的有限单 在此基础上，一个令人感兴趣的问题是：对一 元．例如：由三维弹性理论的块体单元退化得到的 般具有位移c1连续性要求的力学模型，是否都有可 薄板单元；对R／M板单元列式施加Kirchhoff约束 能应用应力函数来避免有限元列式的C1连续性困 得到的离散Kirchhoff薄板单元；对具有独立旋转自