一种具有正交应力模式的杂交膜元.pdfVIP

7l 数学技术应用科学 Mathematics Science Techndogy Applied 一种具有正交应力模式的杂交膜元· 王宇，曹睿 (四川大学数学学院，四川成都610064) 摘要：本文提出了一种具有正交应力模式的四节点四边形杂交元膜元．由于假设应力的正交性，该元与通 常四节点四边形杂交元相比，提高了单元形成过程中的运算速度．而数值试验表明，该元与一些高精度杂交元 如EcQ4具有相当的精度，同时运算量有所降低，是一种简单快速实用的杂交元． 关键词：有限元；正交应力模式；四边形杂交膜元 中图分类号：O241．82文献标识码：A文章编号：7—5416-2336—9(2006)一0071—04 有限元方法自提出以来。得到了长足的发展及 广泛的应用．随着有限元法研究的深入，特别是混 阵，F(孑)=妻(v+v’)孑为应变张量，(盯6，≯) 合杂交元方法的发展，一大批具有高精度的低阶膜 二 ∈矿×扩．变分原理对应的鞍点问题是：寻找(盯‘， 元相继被提出，如P—S元…，ECQ4元K1等．这些 元通过优选应力模式，使低阶位移元达到了较高的 ≯)∈矿×扩使得 精度，由于其列式简单，易于实现的特点得到了广 口(矿“，丁)一6(7-，≯)=0，V丁∈矿(2) 泛的应用．特别的，由【3’中提出的一类等价的宏元 V孑∈扩 b(o-h,孑)=八孑)， (3) HQM4／HQ4，其应力模式满足自平衡条件，因而得 到了较好的精度．本文在其基础上，通过在高阶位 其中u=(u，)7为位移向量，盯分片独立假设的应 力张量． 移部分加上这一因子得到了HQ4m．该元与HQ4 相比，达到了更高的精度，对网格畸变更加不敏感． 对于一般的具有顶点e(x。，Yi)(i=1，…，4) 由于其低阶应力与高阶应力部分正交，使得在单元 的四边形单元K我们有如下的等参变换：FK：定= 形成过程中计算量大大减小于通常得四边形杂交 [一1，1]2一÷K 元膜元并便于并行处理． 其中亭，，7是等参坐标，Ni(i=1，…，4)形式如下． 1 变分格式 7l=U{K}为对区域力的正则剖分，K是相应 单元．杂交膜元基于的变分原理可写成如下形式： a。60 厂1 1 名1 ，，l inf 口l bl 1 I一1 1 茗2 儿 sul!{6(r，孑)一了1口(r，r)一灭孑)}，(1) - 了E驴7‘r’ a2 b2 2了ll —l 石3 儿 L．．1．．1 a3 63 x4 ’4 口(or，r)=厶丁。D～o-dO，6(r，孑))=Lr‘ 让研为通常的双线性协调位移空间，则单元 8(石)m以i)=驴·孑m+虫；·孑山．