无网格法与有限元法的耦合在土木工程中的应用研究.pdfVIP

134 应用科学 蒜霸 无网格法与有限元法的耦合在土木工程中的应用研究 宋伟峰，陈 雷 (安徽淮南矿业集团选煤厂，安徽淮南 232047) 摘 要 采用无网格Galerkin方法与有限元 (EFG—FE)耦合的方法来计算裂纹问题。这种耦合的方法不仅解决了无网格Galerkin法力学边界 条件施加的难点，而且还克服r无网格Galerkin法耗时较多的缺点。数值计算结果表明了方法的有效性。 关键词 无网格Ga】erkm法；EF’G—FE 中圈分类号 TU348 文献标识码 A 文章编号 1673—9671一2(010)081—0134—01 本文采用有限元与无网格耦合的方法，只在裂尖附近区域布置无网 P()=[1，，Y，√。其次，由于在线弹性问题中，J是裂尖处的能 格元，而在其他区域采用一般的有限元。这样，不仅可以充分利用无网 格元精度高，处理裂纹扩展灵活的特点，而且可以容易地处理力学边界 量释放率，因此在求出J后，可由式 (7)求得在平面应力状态下的应力 条件，并提高了求解的效率。 强度因子KI：Ki：J =K ／E (7) 1无网格Galerkin法及其和有限元法的耦合 MLS方法构造出来的形函数并不具有插值性质。本文采用与有限元 计算结果列于表2，并于精确解比较 ，Err为应力强度因子的相对误 耦合的方法来处理这一问题。 差，E ： ％，其中精确解 +为2．358。可以看出，积分回 过渡单元的近似函数为： “()=“ ()+尺()[“ “(X)一“ ()] (1) 路的选择对线性基和二次基计算结果的影响较大，对部分扩展基计算结 式中，R(x)称为坡度函数，它保证了过渡单元实现从无阿格边界 到有限元边界的光滑过渡。它可采用与有限元形函数相同的定义，即： 果的影响较小，部分扩展基由于引入了√项，能较准确地捕捉裂纹尖 端的奇异性。 R()= N，() … J ＼ ￡ 1 xEFgFt； ：上!iL!：L 由式 (1)可以得到过渡单元的形函数为 N1()=[1一R(x)IN，()4-R(x)中， )，x在过渡单元内。 (3) 式中，N，(x)是有限元形函数； ，( 黾EFG方法的形函数。 而过渡单元形函数的导数为 裂纹 K Nr，=一R()Ⅳ，()+1【一 ()]Ⅳf_f()+R() ，() ●——一 +R() ¨() (4) a 式中R ) (川 。 W ， ， 在过渡单元内 XErFF 所以，对于E