有限元计算中疏密网格过渡方法研究_钟红.pdfVIP

第24 卷第6期 Vol. 2 4, No . 6 2007 年12 月 ChineseJournal of Computational Mechanics De ember 2007 : 1007-4708 2007 06-0887-05 * 钟 红 , 林 皋, 胡志强 , 116023 : 工程计算中出于节省计算量的目的, 往往需要在 一个有限元模型中布置粗 不同的网格为保证计算 结果的准确性, 必须保证网格突变情况下的位移协调问题本文工作之 一是在强天驰界面过渡单元的基础上, 引入虚拟节点和子单元, 在子单元中应用节理元思想, 提出了基于最小势能原理的弹簧节理单元法简化了积 分运算, 避免了精度要求极高的坐标转换, 从而提高了方法的精度和实用性; 二是提出了基于位移约束的主从自 由度法, 简便实用, 只需简单的矩阵运算即可实现两种方法均实现了不同尺寸网格间位移的协调性和刚度的 匹配, 从而使之满足有限元收敛准则, 且生成的刚度阵具有对称性及带状性算例证明两种方法精度良好, 并可 方便地应用于求解大规模工程问题 : 疏密网格; 界面协调; 弹簧节理单元; 主从自由度法 : O241. 82 : A [ 1] : 1 , , , , , , [ 2] - , , , , , ns- 1 s 2, 3, , , , , n 3 [ 3] , , , , - , , , , , [ 4] , , , , : 2005- 12-22; : 2006-07-031 , 1% : 座机电话号码 1 100% , * : 1981- , , E-mail: zh onghongm@ eyau. om ; 1929- , ,1 , 888 计 算 力 学 学 报 第24 卷 , 1 vp l2vm + l 1vn 2b [ 5] 2l , , 2l l1 + l2 , 1i,j , m p mp ij , , , x x , , u - um + 1+ up l 1 l 1 x x u - ui + 1+ uj , l 1 l 1 , , x $u u - u - um - ui + l 1 [ 6] 1+ x x , up - uj - um - ui + l 1 l 1 , x l2 l 1 1+ um + un - uj , l 1 2l 2l x 1 1 2 l 1 ui - l1 l 1 + x uj + 2l l 1 + x un + 1 , 1 l1 l2 - 2lx + l2 x um 3a 2ll 1 , , x 1 1 , $v v - v l 1 vi - l1 l 1 + x vj + 2 1 1 l 1 + x vn + l 1l2 - 2lx + l2x vm j 2l 2ll 1 p , , 3b 3a 3b : $D [ M] D e 4a : $u D e [ ui v i uj vj uk v k u n v n um v m ] T $D $v 4b 1 x 1 1 0 - l1+ x 0 0 0 l 1+ x l 1 l 1 2l u v x y p [ M] x 1 m n , 0 0 - l1+ x 0 0 0 l 1 l1 1 up l2 um + l 1un 2a 1 2l 0 l 1l2- 2lx+ l2x 0 2ll 1 1 1 1 1 2 2 l + x 0 l l - 2lx + l x 2l 2ll 1 4 1 2 Fig. 1 Coupling of finite elements and spring j oint elem ents Fig. 2 T he spring joint element 第6 期 钟 红, 等: 有限元计算中疏密网格过渡方法研究 889 , K sj K sj , , , $D y 0 $u $v , kx ky kx ky x y , [ 4] : kx 0 L c 1- K 5 kx ky k L L E 12 0 ky b 1+ 1- 2 L , E , b , c , 10 ~ 100 0 1 2 1 2 5 , , c 40 F - l 2kx $u + 2 ky $v dx Q ~ 200 1 0 kx 0 $u - l $u $v dx 3 2 Q 0 ky $v 0 1 eT T e 3