应用三维有限单元法计算应力强度因子.docxVIP

应用三维有限单元法计算应力强度因子 摘要 描述了两种基于有限单元计算面形裂纹应力强度因子的方法，建议了一?种创造三维有限单 元网格的途径。计算方法的精度通过和其它解析解或数值解的比较得到了说明。 关键词 应力强度因子 有限元 损伤容限设计 断裂评定 无论在损伤容限设计还是在缺陷评处阶段，工程师们需要知道止在分析的构件中裂纹的应力 强度因子，因为判断含裂纹构件的断裂，或者计算剩余疲劳寿命人多依赖丁?这一参量。因此，在 断裂力学发展中，如何求収应力强度因子一直是一个重要的课题。当前已有许多方法可用來计算 应力强度因子，较为典型的有解析法、边界配位法、有限单元法、边界元素法、体力法、权函数 法和线弹簧模型。利川这些方法，人量的应力强度因子解已经获得，已出版的应力强度手册⑷ 中收编了许多典型的解。尽管如此，工程师们仍然会感到H己所需要的应力强度因子解很难找到, 这是因为要解决的工程问题往往是一些受复杂载荷的构件，包含的裂纹也往往是一些不规则裂 纹。 木文简单介绍了两种基于三维有限单元法计算而形裂纹应力强度因子的方法。冇限单元法已 经成为工程设计分析领域中一个强有力的计算?工具，它能模拟非常复杂的构件。基于有限元的应 力强度因子计算方法，白然也将具有卓越的工程能力。除了计算方法的介绍以外，还将简单描述 一种简化网格的生成方法。最后提供了一些所得到的典型应力强度因子解，并和大家熟知的解进 行了比较，以说明木文所描述的方法的可靠性。 1计算方法简介 图1表示了笔者 图1裂纹前?沿单元网格 式。网榕由3个半环共12个20节点三维等参单元组成，每个半环有4个单?元，其中第一个半环 单元的中节点被移金1/4点位置，以模拟裂纹尖端应力和应变场的奇异性皿5] o根据线弹性断裂 理论S3,裂纹尖端的位移场可表示为 壮OS伶(1~加+血誇) 壮OS伶(1~加+血誇) (1) 吩0 (3) 式中，n和t分别为裂纹前沿的法线方向和切线方向；z为SfiT裂纹平面方向(见图2)； u为 材料剪切模量。 图2曲线裂纹前沿坐标系 1. 1 1/4点位移法 根据式(2),若裂纹农血上3二180° )某一点的垂玄于裂纹平血的位移已甸，那么应力强度 因子 用从有限单元法求出的1/4点位移6(^)(见图1)求应力强度因子的方法即为1/4点位移法，即 K = — 4 V r(l/4) 2(1 — P) ⑸ 显然这一方法很简单，但要求裂纹尖端附近的应力应变场能较好地被模拟。 1.2三维J积分法 J积分法实际上是一种能量方法，近來被广泛用來计算应力强度因子，因为两者可按下式转 换 K = yEJ /(1 — V2)(平面应变)(6) 用冇限元计算J积分的方法通常冇2种：回路积分法与虚拟裂纹扩展法。后者最早由Parks [9] 和Hellen[,°-独立提出，主要是通过移动有限元模型的节点位置來模拟裂纹扩展。最近deLorenzi tIb 121已经根据连续介质力学成功地推导了面形裂纹的能量释放计算公式 式中符号请参见文献［11］。商用冇限元分析软件ABAQUS C经收编了上述基于虚拟裂纹扩展 原理”的三维J积分方法，可以直接获得面形裂纹前沿各角节点和中节点的J积分值。 2网格生成 为了方便三维网格的牛成，笔者建议了图3所示的分块方法。整个网格由裂纹块和非裂纹块 图3结构分块示意图 组成，裂纹块包含了裂纹，英网格需要进行精心划分，因为它对应力强度因子有相当大的影响。 相对于裂纹块，非裂纹块显得不是很重要，一般网格可以划分得粗一些，但必须使载荷和约束彳j 一个比较合理的分配。裂纹块和非裂纹块网格的组合可采用“多点位移约束法”，使其位移在其 接触而上基本保持连续，从而使组合后的网格保持变形协调。对丁?复杂含裂纹构件，只要合理地 设计裂纹块和非裂纹块即可大大简化网格的生成。这种网格形式也特别适合于疲劳裂纹扩展的模 拟计算，疲劳裂纹扩展后的新网格重新牛成可以只对裂纹进行，而保持非裂纹块网格不变? 3】。 图4表示了一个笔者建议的三维裂纹块网格的生成过程。首先生成一个8节点等参单元二维 网格，在这个网榕里定义开裂农而和裂纹前沿，然后将其扩展成一个三维网格，单元的扩展有2 种，见图4。一个非裂纹前沿单元扩展成一个20节点单元，而一个裂纹前沿单元扩展成6个20 节点单元。围绕裂纹前沿，网格星放射形。最后组合成的裂纹尖端网格具有图1所示的结构。显 然，这种网格扩展方法进一步简化了三维网格的建立。 图4三维裂纹块网格的形成 3计算程序 笔者最近发展了如图5所示的计算应力强度因子的计算机程序。对于一些常见的工程裂纹, 如平板屮的半椭圆衣而裂纹、压力容器中的半椭圆表血裂纹、平板孔边的角裂纹、圆棒中的表面 裂纹等，只耍输入一些关键的几何参数和载荷参数，程丿芋将口动输出应力强度因子结果。该程序 口前主要和