单位分解有限元方法求解应力强度因子.pdfVIP

·218· 数学-力学·物理学·高新技术研究进展——2006(11)卷 单位分解有限元方法求解应力强度因子 、 射欣 (浙江大学机能学院力学系固体力学研究所，浙江杭州，310027) 【摘要】本文采用基于单位分解概念的有限元方法分析了像断裂力学一类的奇异性问题。该 ’‘ 方法能很自然方便地将奇点附近的奇性函数加入增强函数空间中，并保证所有自由 度都定义在节点上，保持方程良好的带状特性。锝到高耩度的数值解。并可以不经 过后处理过程(如有限元中的J积分)得到断裂力学问题中的重要物理量——应力 强度园子。 【关键词】单位分解；有限元；P型；应力强度园子 0引言 传统的有限元方法在高精度求解断裂力学问题时，主要有两种方式：局部加密网格方法 和奇异函数法。前者通过在奇点附近不断加密网格，甚至采用自适应技术，以取得所需精度， 这种方法有效但费事，且纯h型加密收敛性慢，P型、h-p型可以提高收敛性【l】，但对奇异性 很强的问题，得不到可以接受的结果【l】。后者在有限元空间中引入了奇异函数，能得到较高的 精度，但实际上很少采用，主要是其定义了非节点自由度，不便实施。除上面两种方法外， 的精度和高收敛性。但这种方法需要进行局部区域映射变换。这一处理并不容易。本文研究 基于单位分解概念的有限元方法口】，能很自然方便地将奇点附近的奇性函数加入场量近似函数 空间中，并保证所有自由度都定义在节点上，保持方程直好的带状特性，得到高精度的数值 解。并可以不经过后处理过程(如有限元中的．，积分)得到像断裂力学问题中的重要物理量 ——应力强度因子或其他问题(流体)的涡流强度因子。 增强型函数的实现 在奇性问题中，往往在奇点附近的局部区域内，解析解的 r-，，—lI』 形式是已知的。如图I所示的带凹角的奇性问题．在奇点附近， ’ 解析形式具有典型的形式例 图1平面裂纹问题 作者简丹t剃段(196}一)，博士．浙江太掌力学暮固体力学研究所副教授，研究方向为计算力学，计算随机动 力学及其应用． 单位分解有限元方法求解应力强度固子 ·219· Ⅳ=f110+∑4一z(印f ① L 仁I．。 J 。(r’目)是以奇点为原点的极坐标，岛是特征值，且啦≤嘞≤∞≤…，抓目)是目的解析函数， 是特征函数。蜘，丑是展开系数，图1是常见的奇性阀题，及奇性指数伪。当O≤岛≤l时， 奇性出现，嘶在(0，1)内越小，奇性越强，熟知而重要的是裂纹问题，Ql=O．5。 纯粹的解析方法往往只在奇点邻域具有好的逼近精度。但在整个求解区域精度可能很差。 并且解析方法常常只对规则域的偏微分方程有效，对复杂域往往无能为力，需要数值方法。 不考虑场量理性成分(物理含义)的纯多项式逼近，虽有其优点，但也不可避免存在一些不 尽如人意的地方Hjl。在人们对问题本质有了深入的了解，将这一了解加入数值分析中，将会 取得满意的结果。这便是半解析形式的数值方弦各受关注的原因【4j】。下面以角点奇异性为例 来说明如何进行非多项式或增强型多项式增强函数的实现。 。 ，平面弹性问题含角点(凹角)奇异性问题。假设物体体力忽略不计，应变能有限及角点 两边是零面力，且不存在刚体位移(约束了)。。可以得到其解析形式解口】，分对称形式和反对 称形式两部分。位移解 _÷》薅象鸳塑壁氅簪塑墼坠 ～一善象垫塑坐桨警掣；}+ ② 垆一+善等掣霉擎等等垡坐翌 砉象叠塑盟糕雾璺鲤塑 上式中G是剪切模量；k是依赖泊松比v的常数。 平面应力：