采用无单元法计算含边沿裂纹功能梯度材料板的应力强度因子.pdfVIP

第17 卷第 5 期 工 程 力 学 Vol.17 No.5 2000 年 10 月 ENGINEERING MECHANICS Oct. 2000 文章编号 1000-4750(2000)05-139-06 采用无单元法计算含边沿裂纹 功能梯度材料板的应力强度因子 陈 建 吴林志 杜善义 (哈尔滨工业大学复合材料研究所 哈尔滨 150001) 摘 要 本文采用无单元法分析了功能梯度材料的断裂力学问题 无单元法采用基于滑动最 小二乘近似的位移插值形式 节点布置变得非常自由 这种插值形式不仅很好地反映了材料 变形 而且使得无单元法在分析功能梯度材料时可以方便地采用各个积分点处的材料特性 数值计算结果表明无单元法在分析功能梯度材料力学行为方面具有较高的效率和精度 关键词 无单元法(EFM) 滑动最小二乘法 功能梯度材料(FGM) 应力强度因子(SIF) 中图分类号 O343.4 文献标识码 A ?? 前言 功能梯度材料(FGM)的可设计材料特性大大改善了它的热力学性能 从而受到工程界 的广泛重视 但由于 FGM 材料的非均匀性使得其力学行为相当复杂,尤其是几何构形或受 力复杂时 试图用解析法分析其力学行为相当困难 常常只能借助数值方法 有限元法通常采用均匀材料单元 在分析非均匀材料时需要不断细化单元 李春雨和 邹振祝[1]认识到均匀材料单元实际上是在单元内部采用零次材料特性插值 即使采用高阶 位移插值对非均匀材料也很难提高计算精度 为此 他们提出了材料特性采用与位移相同 插值形式的多重等参单元的有限元法 实际上 分别取单元内部各积分点处的材料特性形 成单元刚度阵也是多重单元的一种特殊形式 同样会提高有限元法计算精度 近年来兴起的无单元法(EFM)[2,3]非常适合分析类似 FGM 之类的非均匀材料的力学行 为 首先 EFM 的滑动最小二乘近似使得其位移插值摆脱了有限元法那样严格的单元概念 即使同一积分子区域内不同积分点处的位移形状函数及其偏导都不一样 因此要分别形 成各个积分点处的刚度阵 这样 EFM 可以方便地采用各个积分点处的材料特性来形成刚 度矩阵以反映材料的真实特性 此外 EFM 的位移插值形式不仅能够比较真实地反映材料 变形 而且使得节点布置变得非常自由 EFM 可以在相同的积分区域划分形式下根据需要 布置节点 而不需进行单元重分 这非常适合模拟裂纹扩展 2 无单元法(EFM) 考虑一线弹性体 ? 应力边界为 Γ t 位移边界为 Γu 平衡方程和边界条件分别为 收稿日期 1999-12-25 修改日期 2000-03-04 作者简介 陈建(1971) 男 四川省人 博士 140 工 程 力 学 σ ij, j + bi = 0 在?内 (1) σ ij n j = t 在Γt 上 (2) ui = ui 在Γu 上 (3) 其