PCB热循环测试中循环次数与通孔铜壁疲劳裂纹扩展长度的Paris拟合

一、Paris公式的理论基础与物理适用性

Paris公式（帕里斯公式）是断裂力学中描述疲劳裂纹扩展速率的经典经验方程，其基本形式为：da/dN = C·(ΔK)^m。其中a为裂纹长度，N为应力循环次数，ΔK为应力强度因子范围，C和m为与材料相关的常数。在双对数坐标系中，lg(da/dN)与lg(ΔK)呈线性关系，该线性区间对应疲劳裂纹扩展的第II阶段（稳态扩展阶段）。

通孔铜壁在热循环疲劳中的裂纹扩展符合Paris公式的物理假设。热循环过程中，PCB基板与铜镀层之间因热膨胀系数（CTE）失配产生循环热应力，孔壁铜层经历周期性拉伸-压缩应变。当裂纹尖端应力强度因子范围ΔK超过门槛值ΔKth后，裂纹萌生并进入稳定扩展期，适用Paris公式描述。

研究表明，Paris公式已被广泛应用于电子封装领域的疲劳寿命预测。De?planque等在功率晶体管焊点疲劳寿命预测中采用基于Paris定律的裂纹扩展有限元方法。Samet等通过双悬臂梁试样的循环加载实验，测定环氧模塑料/铜界面的Paris常数，验证了该方法对界面疲劳裂纹的适用性。倒装芯片底充胶/芯片界面的分层裂纹扩展速率同样可用Paris半经验方程描述。

二、热循环测试中裂纹扩展数据的获取方法

**测试样品与条件设定**

热循环测试样品应为多层PCB通孔结构，推荐孔直径0.2-0.3mm，板厚1.6-2.4mm，厚径比8:1至12:1。热循环条件可参考IPC-9701标准：温度范围-40℃至125℃（或-55℃至125℃），升降温速率10-15℃/min，驻留时间10-15分钟。

**裂纹长度测量方法**

采用“单组试样法”获取裂纹扩展数据，即在多个特定循环次数节点（如0、200、500、1000、1500、2000次）分别取出样品进行金相切片。每组节点至少3个样品，测量孔壁铜层裂纹长度，包括最大裂纹长度（沿孔壁纵向延伸）、裂纹深度（从孔口或孔中向镀层内部扩展）以及裂纹周向占比（裂纹环绕孔周的弧度比例）。

为确保测量精度，取样间隔不大于试验总循环次数的10%。裂纹长度测量应在200-500倍金相显微镜下进行，记录裂纹从起始点到最远延伸端的距离。当裂纹扩展至孔壁周向的50%或贯穿孔壁铜层厚度时，认为该通孔已失效。

三、ΔK的计算方法

对于通孔铜壁，裂纹尖端的应力强度因子范围ΔK需通过有限元或解析方法计算。

**简化解析模型**

将通孔简化为含裂纹的圆筒结构，在热循环温差ΔT作用下，孔壁承受的环向应力范围Δσ为：Δσ = (α_sub - α_Cu)·ΔT·E_Cu，其中α_sub为基板Z轴CTE（典型值50-70ppm/°C），α_Cu为铜的CTE（约17ppm/°C），ΔT为热循环温差，E_Cu为铜的弹性模量（约110GPa）。

考虑孔壁几何和裂纹几何的修正因子F(a,W)，ΔK可表示为：ΔK = Δσ·√(πa)·F(a,W)。其中a为裂纹长度，W为孔壁铜层厚度（典型值20-25μm）或孔壁剩余韧带宽度。

**有限元计算方法**

更精确的ΔK计算可采用有限元方法。建立轴对称或三维通孔模型，施加热循环载荷，提取裂纹尖端的J积分或能量释放率G。ΔK与ΔG的关系为：ΔK = √(ΔG·E/(1-ν²))（平面应变条件）。Paris公式也可用能量释放率范围ΔG表达：da/dN = C_G·(ΔG)^m_G，其中m_G ≈ m/2。

四、Paris参数的拟合方法

**数据转换与线性回归**

在稳定扩展阶段（裂纹长度a约5-20μm），每个循环次数节点对应一组测量数据（N_i, a_i）。计算各节点的裂纹扩展速率：取相邻两节点的裂纹长度差Δa_i = a_{i+1} - a_i，对应循环次数差ΔN_i = N_{i+1} - N_i，则扩展速率(da/dN)_i = Δa_i / ΔN_i。

计算各节点的ΔK_i，对Paris公式两边取对数：lg(da/dN) = lg(C) + m·lg(ΔK)。以lg(ΔK)为横坐标，lg(da/dN)为纵坐标绘制散点图，进行线性回归分析。

**回归结果解读**

回归直线的斜率即为指数m，截距即为lg(C)。焊料合金的典型Paris参数范围为：m值2-6，C值1e-8至1e-6（单位取决于ΔK单位）。环氧/铜界面的Paris参数也存在类似关系。

拟合优度R²应大于0.95，否则可能需增加数据点或检查数据一致性。统计显著性P值应小于0.05，确保参数估计可靠性。置信区间评估可给出C和m的95%置信区间，用于寿命预测中的不确定性分析。

五、典型参数范围与工程验证

基于焊点及互连结构的裂纹扩展研究，Da/W为裂纹长度与特征尺寸之比。相关系数r通常>0.98，表明Paris公式对通孔铜壁疲劳裂纹扩展具有良好的拟合效果。剩余标准差SD < 0.05，说明模型预测精度较高。

在工程应用中，利用拟合得到的Paris参数可预测给定热循环次数N下的预期裂纹长度a_N，为可靠性评估提供定量依据。也可根据允许的裂纹长度反推安全循环次数，指导产品寿命设计和筛选测试规格制定。

六、结论

Paris公式可有效描述PCB通孔铜壁在热循环测试中的疲劳裂纹扩展行为。通过金相切片获取不同循环次数的裂纹长度数据，计算裂纹扩展速率da/dN，结合有限元或解析方法计算裂纹尖端应力强度因子范围ΔK，进行双对数线性回归即可拟合得到材料常数C和m。该拟合结果为PCB通孔的热疲劳寿命预测和可靠性设计提供了量化依据。