基于Hamilton变分原理的四边形4结点杂交元.pdfVIP

i I 基于Hamton变分原理的四边形4结点杂交元 张涛1，‘岑松1’2 (1清of；人学航天航空学院1：程力学系，北京100084：2清华人学航天航空学院．AML．北京100084) 摘要： 本义从Hamilton变分原理}|：发，利用泛函驻值条件，给…了其丰}I应的二：维有限厄格式，并钍此基础Ii 选取适当的应力试解，构造了两个4结点杂交单元HAMH4．1和HAMH4．II。数值算例表明：一方面，与传统的4 结点等参元04相比，这些单元精度高；另一方面，与经典的杂交元相比，Hamilton杂交元采用的应力分馈由二三 个变为两个．使计算规模有所降低，显示了Hamilton杂交元的优越性。这为进一步构造高性能有限元模型提供了 新的思路。 关键词： 有限元；Hamilton变分原理；Hamilton杂交元；应力试解 1 引言 一方面，其精度远优于传统的4结点平面单元；另一方面，其应力试解中所需的应力变量个数比传统杂交 元要少，因此有较高的计算效率。本文延续了文献【2】的思路，采用不同的应力试解，构造-『两种新型 Hamilton四边形杂交元，经数值计算，结果体现出其优良的性能。 2 Ham．I ton能量泛函14J 对于二维各向同性的弹性力学平面问题，首先写出基于Hellinge卜Reissner变分原理的能量泛函： 利用平面应力问题本构关系 盯，．=∥盯，+E兰‘ (2) ’ 妙 3 Hamil ton杂交元的一般格式【2l 令单兀的位移场祁肫力场其有如卜的形式： 黔咐 阱M ㈣ 其中q。,fu∥分别为单元的结点位移矢量和应力参数矢量。川和∥伪相应的形函数和应力试解钳i阵。将式 (4)代入式(3)，并不计体力，可得： 17JI，：三．q。TSq。+∥TGq。一i1∥T月∥一g，R (5) ‘苹龠项}I：闻家自然科学坫金资助项t1，教白鄙新世纪优秀人，j‘支持i1划顼II(NUt1’I．07-0477) f’：肯简介：张涛(1984--}，锝，黑龙j．I．人．硕i：，|．}耍从事计算㈣体力学研究： +岑松(1972一}．男．榀建人，剐教授．}I孽导．h委从事结构I：稃干uitt：【|Il^体力学研亢“：一madeducn ccnsong似sIn曲uaJ。 ·I．530· c6a，b， 扎《￡)2删=gJ却。；儿[吒(象+∥霉)+勺(雾+2)]幽=∥。Gg。； JJl壶k(1一∥2)+2(1+胁刍k=∥7琊； (6c’ 由驻值原理加Ⅳ／op=0得 琊一Gg。=0 (7) j∥=H叫Gq8 g。TK8叮。一98T且 (8) 仃Ⅳ=iI K2=S+GTH一1G (9)