双周期平行四边形排列裂纹反平面问题的有限元分析.pdfVIP

双周期平行四边形排列裂纹反平面问题的有限元分析· 皮赞严鹏蒋持平 (北京航空航天大学固体力学研究所， 100083) 摘要发展了双周期平行四边形排列裂纹反平面问题的有限元方法。通过对单位胞元引入 周期边界条件，在裂纹尖端采用奇异单元，解决了有限元分析这类问题的效率和精度问题． 将现有限于对矩形阵列的双周期裂纹的研究推广到了非对称的平行四边形阵列的研究，揭 示了这种非对称性对断裂参数与材料的有效性质的影响规律． 关键词：双周期有限元裂纹反平面问题应力强度因子 一、 己I喜 J1日 双周期裂纹是研究裂纹相互作用的重要力学模型。大量工程材料和结构，特别是复 合材料和岩石结构，提供了双周期裂纹模型的实际应用背景。 J运用复变函数的共形映照技术与椭圆函数理论，获得了双周期对称 最近，Tong等Il 排列裂纹反平面问题的含椭圆函数的封闭解，为评估数值方法的精度提供了参照解，但 是解析公式的应用范围有局限性。 本文对一般双周期平行四边形排列裂纹的胞元引入周期边界条件¨J，解决了由基本 胞元精确模拟无穷双周期裂纹体的问题。同时在裂纹尖端使用奇异单元，进一步提高了 计算精度。各种有解析解对照的情形，断裂参数计算的相对误差均在0．2％以下，实现 了计算便捷、效率与精度的统一。 模型及周期性边界条件 2．1双周期平行四边形排列裂纹的反平面剪切问题的描述 f @@@@@@ ■=昌E=曹e=jEj皇 一_，，；．，-一 三．．-．，。 0ppQQO (口) 图1 双周期裂纹：(a)双周期平行四边形排列：(b)代表性体积单元(单胞) 图Ka)是双周期平行四边形排列裂纹横截面示意图，基体为各向同性线弹性材料， 远场受反平面切应力f作用。图l(b)为一个单位胞元平行四边形横截面示意图。 ’国家自然科学基助项目 —-137—— 2．2单位胞元的周期性边界条件 根据Xia等【2J的研究，我们可以从问题的周期性来解决这个问题。考虑更…般的空间 三周期分布问题，为了保证相邻胞元不会相互嵌入和相互分离，各重复的胞元边界位移 是双准周期的，即胞元～对平行边界上各对应点的位移差是～个常数，能够由平均应变 表示。 三、算例及分析 3．1裂纹的分布参数对应力强度因子的影响规律 取单位胞元平行四边形横截面(如图1所示)底边上的高形口=l，w,／a=3，4，5，6，10， 无量纲应力强度因子墨。I。=K。I／(f√磊)随曰角度的变化画在图2(a)中；取山h／a=10， Wl／a=3，4，5，6，10，无量纲应力强度因子KIl+随口角度的变化画在图2(b)中； 曲 ，蓝 ∞ 基 t越 {搴 ，∞ ?18 ，∞ §辩 n蜘 O-S 4幻 O讳 07口 0辫 臀 鄹 W 饲 蜀 孵 ∞ w 西 研 ∞ ∞ ∞ Ⅺ 萄 0 0 (a)h／a=l时，Kin*随口的变化曲线 ‘(b)h／a=10时，Kill*随秒的变化曲线 图2双周期平行四边形阵列裂纹倾斜角对应力强度因子的影响 本例中角度的变化可以看作双周期裂纹的各行裂纹问的相互错动。如果叠间距在其 相互影响范罔内时，双周期裂纹的错动效应对应力强度因子影响很大；反之，如果叠间 距在其相互影响范围外时，双周期裂纹的错动效应对