磁电热弹耦合材料三维多裂纹超奇异积分法.docVIP

第４０卷第１期２００８年１月 Ｃｈｉｎｅｓｅ 力学学报 Ｖｂｌ．４０．Ｎｏ．１ Ｊａｎ．，２００８ ＪｏｕｒｎａｌｏｆＴｈｅｏｒｅｔｉｃａｌａｎｄＡｐｐｌｉｅｄＭｅｃｈａｎｉｃｓ 磁电热弹耦合材料三维多裂纹超奇异积分法 朱伯靖＋，ｔ秦太验ｕ’ ＋（中国农业大学理学院，北京１０００８３） 十（中国科学院研究生院地球科学学院，北京１０００４９） 摘要用超奇异积分方程法将多场耦合载荷作用下磁电热弹耦合材料内含任意形状和位置三维多裂纹问题转化为求解一以广义位移间断为未知函数的超奇异积分方程组问题，退化得到内含任意形状平行三维多裂纹问题的超奇异积分方程组；推导出平行三维多裂纹问题的裂纹前沿广义奇异应力场解析表达式、定义了广义（应力、 应变能）强度因子和广义能量释放率；应用有限部积分概念及体积力法，为超奇异积分方程组建立了数值求解 方法，编制了ＦＯＲＴＲＡＮ程序，以平行双裂纹为例，通过典型算例，研究了广义（应力、应变能）强度因子随 裂纹位置、裂纹形状及材料参数变化规律，得到裂纹断裂评定准则．最后，分析了裂纹ｆ，Ｔｑ：扰、屏蔽作用及其 在工程实际中的应用． 关键词三维多裂纹，磁电热弹耦合材料，超奇异积分，应力强度因子中图分类号：０３４６．１，０３０２ 文献标识码： Ａ 文章编号：０４５９－１８７９（２００８）０１—００４６—１３ 引言 从Ｂｒｅｗｓｔｅｒ，Ｃｕｒｉｅ兄弟，Ｖａｓｌｅｋ到Ｄａｖｉｓ，Ｖａｎ— 合材料内含任意形状和位置三维多裂纹问题转化为求解一以广义位移间断为未知函数的超奇异积分方程组问题，退化得到内含任意形状平行三维多裂纹问题的超奇异积分方程组；应用主部分析法，通过分析平行多裂纹问题的裂纹前沿广义应力奇性指数，得到了广义奇性应力场解析表达式，定义了广义（应力、应变能）强度因子和能量释放率；应用有限部积分概念和体积力法，对超奇异积分方程组建立了数值方法，编写了广义（应力、应变能）强度因子计算专 ｄｅｎｂｏｏｍａｇａａｒｄ，Ｔｉｎｋｈａｍ，Ａｂｏｕｄｉ，磁电热弹耦合问题一直受到广泛关注［１￣６】；Ａｂｏｕｄｉ，Ｇａｏ等［７￣１１】对磁电热弹耦合材料进行了基础性理论研究； Ｚｈｅｎ 和Ｚｈｏｕ建立了电磁材料结构非线性力学行为的完整理论框架，预测了铁磁材料典型实验现象磁力表达式［１２，１３】，给出了多场耦合非线性本构关系［１４，１５】，对电磁材料结构行为给出定量预测［１６，１７】；Ｚｈｏｕ等对磁电热弹耦合作用下磁通跳跃现象给出成功预测［１８，１９】；Ｓｉｈ［２０￣２３】最早提出应变能强度因子理论， 用ＦＯＲＴＲＡＮ程序；以典型平行三维双裂纹为例， 研究了广义（应力、应变能）强度因子与裂纹形状、裂纹间位置关系，得到裂纹断裂评定准则，分析了多裂纹间干扰和屏蔽作用及其在工程实际中的应用． 并把它应用到复合材料断裂力学研究中，做了许多 开创性工作…；杨卫［２５，２６】，Ｌｉ［２７】，刘又文［２８】等学 者也从不同方面对磁电热弹耦合材料研究做出了突出贡献．以上这些成果为磁电热弹耦合材料断裂问题研究奠定了良好的理论基础．但鉴于数值方法、理论模型、数学力学和物理学上的困难，特别是数值 １数学模型及超奇异积分方程组 如图１所示，设磁电热弹耦合材料空间区域口中含有任意形状和位置的三维多裂纹．选择整体坐标系为空间直角坐标系耽（ｉ＝１，２，３），在裂纹面＆（ｉ＝ 方法（计算量与精度）、理论模型（三维裂纹扩展模型）、数学（公式庞杂）及力学和物理学（问题复杂） 上的困难，两相或多相耦合场作用方面研究虽然取得一些成果［２０，２９“３７］，但三维多裂纹干扰及相互作用问题的研究较少． 本文根据广义位移基本解【３８】，将磁电热弹耦 ２００７－０６—１３收到第１稿，２００７－０９—２４收到修改稿１）Ｆ—ｍａｉｌ：ｔｙｑｉｎ＠ｃａｕ．ｅｄｕ．ｃａ １，２，…，ｎ）表面上分别作用力载荷Ｐ；（Ｊ＝１，２，３），电 载荷ｑ‘，磁载荷ｂ‘和热载荷Ｐｉ．为了研究问题方便，建立局部坐标系ｏｉ，裂纹ｉ位于局部系坐标轴茁ｊ ｚ； 构成的平面内，同时与局部系坐标轴。§垂直，裂纹 中心点在局部系中坐标为０。（鬈），局部系坐标轴茁： 万方数据　 第１期朱伯靖等：磁电热弹耦合材料三维多裂纹超奇异积分法４７ ２（∈２） 图１磁电热弹耦合材料内含任意形状和位置三维多裂纹示意图 Ｆｉｇ．１Ｍｕｌｔｉｐｌｅｉｎｔｅｒａｃｔｉｎｇｄｉｓｌｏｃａｔｉｏｎｓｉｎｅｌｅｃｔｒｏｍａｇｎｅｔｏｔｈｅｒｍｏｅｌａｓｔｉｃｍｕｌｔｉｐｈａｓｅｃｏｍｐｏｓｉｔｅｓ 与整体系坐标轴ｚｉ夹角０ｉ（ｚｔ，ｚｉ）．利用Ｓｏｍｉｇｌｉａｎａ广义位移巩（ｐ）可表示为 公式，磁电热弹耦合材料中任意一点ｐ（ｘｌ，Ｘ２，Ｘ３）处 ＵＩ＝｜ＵｏｊｆｊｄＶ一｜∞ＪＵ３＋Ｕ｛ＪＴｏ、）ｄ８＋ＪｎＪｒ （Ｚ０：嗡巧一嗡吩ｄｓｔ扣！％巧一嗡吩ｄ、