Paris法则 — CAE术语解说

分类：术语集 | 2026-01-15

CAE visualization for paris law - technical simulation diagram

Paris法则

🧑🎓

老师，Paris法则的公式 da/dN = C(ΔK)^m 在教科书中经常出现，这个公式的含义是什么呢？

Paris法则的理论基础

Paris法则的基本概念

🧑🎓

教科书中说「Paris法则是表示裂纹扩展速度的经验法则」，具体是通过什么「经验」获得的公式呢？

🎓

好问题。1960年代，P.C. Paris等人将铝合金和钢材的大量疲劳试验数据绘制成图，发现裂纹扩展速度

$$ da/dN $$

和应力强度因子范围

$$ \Delta K $$

的对数图上，中间领域呈直线。也就是说，

$$ \log(da/dN) $$

和

$$ \log(\Delta K) $$

的关系呈线性。这个「经验」基础导出了Paris法则的基本形式

$$ da/dN = C (\Delta K)^m $$

。

🧑🎓

应力强度因子范围

$$ \Delta K $$

怎样求得？仅从最大应力和最小应力就能确定吗？

🎓

基本上从最大应力和最小应力求出的应力范围

$$ \Delta \sigma = \sigma_{max} - \sigma_{min} $$

、现在的裂纹长度

$$ a $$

和表示几何形状的修正系数

$$ Y $$

来计算。公式是

$$ \Delta K = Y \Delta \sigma \sqrt{\pi a} $$

。这个

$$ Y $$

，在无限板中心裂纹时为1.0，但实际结构物因有限宽度和自由端的影响，会比1.0大。例如，宽度50mm的板端部有10mm的裂纹时，

$$ Y $$

约为1.12。

🧑🎓

材料常数

$$ C $$

和

$$ m $$

，比如普通结构用钢（SM490）大概是什么值呢？

🎓

通过ASTM E647等标准试验方法来测定。作为粗略目标，空气中的钢铁材料的指数

$$ m $$

通常在2到4的范围内。例如，某种SM490钢的数据表明，

$$ C = 6.9 \times 10^{-12} $$

(单位为 m/cycle, MPa√m)、

$$ m = 3.0 $$

。但这仅是一个例子，实际设计中需要从对象材料的试验数据或信任度高的数据库（如日本材料学会疲劳数据表）中选定。

Paris法则的数值计算方法

CAE软件中的处理方法

🧑🎓

在CAE软件中进行疲劳裂纹分析时，Paris法则具体在哪个模块或功能中使用？

🎓

主要在「裂纹扩展分析」或「疲劳裂纹分析」专用模块中使用。例如，Abaqus的「Abaqus/Standard」有「Contour Integral」计算功能可求得

$$ \Delta K $$

，结合「Direct Cyclic」分析步，可执行基于Paris法则的裂纹扩展模拟。Ansys中，「Ansys Mechanical」的「Fracture」工具箱内的「Fatigue Crack Growth」功能相当。

🧑🎓

分析中，随着裂纹扩展需要更新网格吗？还是软件自动完成？

🎓

取决于方法。主流有「网格重新生成法」和「扩展有限元法（XFEM）」。网格重新生成法是裂纹尖端向前推进时，删除周围网格，按新的裂纹形状重新生成网格。许多软件（如FRANC3D等）都实现了自动化，但计算成本较高。另一方面，Ansys和Abaqus中实现的XFEM，裂纹穿过单元而扩展，无需更新网格。但是，裂纹扩展方向的预测算法（如最大主应力方向）的设置很关键。

🧑🎓

Paris法则的公式

$$ da/dN = C (\Delta K)^m $$

在求解器内部怎样积分的，裂纹长度增量是如何计算的？

🎓

通常采用欧拉法等简单的数值积分。例如，某个荷载循环块

$$ \Delta N $$

内的裂纹扩展量

$$ \Delta a $$

从当前裂纹长度

$$ a_i $$

求得的

$$ \Delta K_i $$

计算为，

$$ \Delta a = (da/dN)_i \times \Delta N = C (\Delta K_i)^m \times \Delta N $$

。然后用

$$ a_{i+1} = a_i + \Delta a $$

更新裂纹长度。这个

$$ \Delta N $$

的取值直接关系到分析的精度和效率，软件通常允许用户设定「扩展增量」。

Paris法则的实务应用

分析工作流程和注意事项

🧑🎓

实际要用Paris法则评价飞机零部件的裂纹扩展寿命，最初应该从什么开始？

🎓

首先是「初期裂纹尺寸的定义」。这由规格确定。例如，航空航天领域要求损伤容限设计，初期缺陷尺寸基于制造检测能力来确定。如果可以进行毫米波探伤检查（NDI），根据FAA的公告圆形，初期裂纹长度可能设定为1.27mm（0.05英寸）。这个初始值若设置错误，寿命预测就会与现实相差甚远。

🧑🎓

荷载数据怎样准备？实际飞行数据像随机振动一样，波形不规则呢。

🎓

完全正确。实际运行荷载历史不规则。这时采用「雨流计数法」等方法，将应力范围

$$ \Delta \sigma $$

和平均应力组合统计化。之后，为了考虑平均应力的影响，用有效应力强度因子范围代替

$$ \Delta K $$

，或者在Paris法则中引入Walker公式等平均应力修正。软件（如nCode DesignLife）中，雨流计数和疲劳裂纹扩展计算的联动工作流程已经构建好了。

🧑🎓

分析结果输出「还能持续多少周期才破坏」的寿命，怎样验证这个结果的可信度？

🎓

绝对的验证只有实物试验。但有验证分析合理性的检查清单。1) 使用的材料常数

$$ C, m $$

是否与分析对象的环境（湿度、温度）相符。2) 裂纹扩展过程中，

$$ \Delta K $$

是否低于材料的阈值

$$ \Delta K_{th} $$

，或者接近断裂韧性

$$ K_{IC} $$

（Paris法则不适用的区域）。3) 裂纹扩展路径在物理上是否合理（例如主要是Mode I开口型）。不进行这些检查而仅凭寿命数值是危险的。

Paris法则的软件比较

主要软件的功能和特点

🧑🎓

Ansys、Abaqus、NASGRO（前NASA代码）这三个软件中，使用Paris法则进行裂纹扩展分析的方法有什么不同？

🎓

Paris法则的核心公式相同，但集成度和数据库差别很大。Ansys Mechanical将「Fatigue Crack Growth」工具集成到通用结构分析中，用XFEM实现与网格无关的分析。Abaqus也类似，Python脚本的定制性很强。另一方面，NASGRO（现在是商用软件）是为航空航天领域专门开发的，最大优势是拥有庞大的材料裂纹扩展数据库和包含Paris法则在内的更高级的扩展公式（NASGRO方程）。可以使用符合规格（ESDU、MMPDS）的数据。

🧑🎓

COMSOL Multiphysics也能做疲劳裂纹分析吗？那么Paris法则的实现是怎样的呢？

🎓

COMSOL的「Structural Mechanics Module」中的「Fracture Mechanics」功能可以计算J积分和应力强度因子。但据我所知，内置的「疲劳裂纹扩展求解器」不像Ansys或Abaqus那样直接。相反，「用Paris法则定义变量

$$ da/dN $$

，随时间（或循环）更新裂纹长度参数」的过程需要用户利用COMSOL的「Events」界面或「ODE（常微分方程）界面」自己构建。灵活性高，但需要专业知识。

🧑🎓

免费或开源软件（CalculiX、Code_Aster等）能进行Paris法则分析吗？

🎓

可以，但大部分操作依赖于手动或脚本。例如，Code_Aster（Salome-MECA平台）具有非常强大的裂纹力学功能，

$$ K $$

的计算是可以的。但是，要执行基于Paris法则自动扩展裂纹、更新网格的一系列过程，需要正确组合专用命令（在DEFI_MATERIAU中定义扩展公式，使用PROPAGATION_MAILLAGE等）。离商用软件的图形界面一键操作还很远，需要对理论和软件都深入理解。

Paris法则的故障排除

常见错误和对应方法

🧑🎓

分析执行时，裂纹突然向意外方向弯曲扩展了。可能的原因是什么？

🎓

主要有两个原因。第一，裂纹尖端应力场的计算精度不足。特别是裂纹尖端的网格太粗时，应力强度因子

$$ K_{II} $$

（面内剪切模式）的计算值会变得不稳定，扩展方向预测出现偏差。第二，扩展方向判定基准的问题。许多软件采用最大主应力方向（

$$ \theta = 2 \arctan\left(\frac{1}{4}\frac{K_I}{K_{II}} \pm \sqrt{(\frac{K_I}{K_{II}})^2 + 8}\right) $$

之类的公式），但在复杂多轴应力状态下，这个基准本身可能不适用。首先应细化裂纹尖端的网格，确认是否使用了奇异单元（singularity element）。

🧑🎓

寿命预测结果与过去的实验数据相比，极端地短（或长）。材料常数是从实验数据中取得的，为什么会这么大的偏差？

🎓

Paris法则公式

$$ da/dN = C (\Delta K)^m $$

仅表示对数图上的直线部分。偏差最常见的原因是分析中计算的

$$ \Delta K $$

和实验中使用的标准试验片（CT试验片等）的

$$ \Delta K $$

「是否用同一计算方法求得」的忽视。特别是几何形状修正系数

$$ Y $$

是否正确。而且，随着分析对象结构物中裂纹的扩展，

$$ Y $$

会大幅变化。正确反映这个

$$ Y $$

变化的「裂纹扩展曲线」与常数

$$ Y $$

的简单计算，寿命相差可能达数倍。

🧑🎓

用XFEM进行分析，裂纹扩展计算在某个循环数时突然停止，出现「不收敛」的错误。

🎓

XFEM中收敛失败很常见。直接原因是裂纹尖端到达单元内特定位置（如高斯积分点附近）时，水平集函数更新中出现数值不稳定性。应对措施有几个。1) 减小裂纹扩展增量

$$ \Delta a $$

。若默认值为0.001m，试试改成0.0005m。2) 改变求解器设置。Abaqus/Standard中，在「步」设置中启用「允许更大的增量尺寸」，或提高迭代矩阵的更新频率。3) 根本方法是把裂纹扩展路径附近的网格改为均匀细密的网格。XFEM号称与网格无关，但粗网格容易导致扩展行为不稳定。

Paris法则 — CAE术语解说

Paris法则

Paris法则的理论基础

Paris法则的基本概念

Paris法则的数值计算方法

CAE软件中的处理方法

Paris法则的实务应用

分析工作流程和注意事项

Paris法则的软件比较

主要软件的功能和特点

Paris法则的故障排除

常见错误和对应方法

相关主题