[发明专利]双覆盖层结构涡流检测方法

权利要求书

1.一种双覆盖层结构涡流检测方法，其特征在于，包括以下步骤：

S1、建立双覆盖层加基体的涡流检测的理论模型，获得试样的阻抗的理论值Z_C；

S2、建立与理论模型的参数一致的实体模型，所述试样的阻抗的实测值Z_Cmeasured；

S3、建立评价函数

其中，Z_C为步骤S1中所述试样的阻抗的理论值，Z_Cmeasured为步骤S2中所述试样的阻抗的实测值；m为测量频率的个数，即扫频涡流检测测量数据点数，β为变量，β＝(z₁，d₁，σ₁，μ₁，d₂，σ₂，μ₂，d₃，σ₃，μ₃)，z₁为提离距离，表示线圈的下端面至试样的上表面沿线圈的轴线方向的距离，d₁为第一层的厚度，μ₁为第一层的相对磁导率，σ₁为第一层的电导率，d₂为第二层的厚度，μ₂为第二层的相对磁导率，σ₂为第二层的电导率，d₃为第三层的厚度，μ₃为第三层的相对磁导率，σ₃为第三层的电导率，f为测量频率；

S4、计算评价函数F(β，f)的雅克比矩阵J(β)和海赛矩阵H(β)，得到评价函数F(β，f)的梯度函数

以及海赛矩阵

H(β)＝J^TJ (3)

其中，F′(β)为F(β，f)的关于β的梯度函数，ε(β)＝Z_C-Z_Cmeasured，ω为电磁激励的角频率；

S5、设置参数β_k，采用LM算法的迭代步进h_lm进行迭代优化，获得第k次迭代的局部最小值β_local-k；其中，

h_lm＝-(H(β)+μI)^-1g＝-(J^TJ+μI)^-1g (4)

μ为阻尼因子，μ＞0；

S6、将局部最小值β_local-k带入海赛矩阵H(β)，计算对应的海赛矩阵H_k；

S7、对H_k进行Cholesky分解，判断H_k是否为正定矩阵；若是，则执行步骤S8；若否，则执行步骤S9；

S8、将局部最小值β_local-k存储为一个局部最小确认值，执行步骤S9；

S9、根据迭代次数k≤K_max判断是否需要进行下一次迭代，其中K_max是设定的最大迭代次数；若是，则执行步骤S10；若否，则执行步骤S11；

S10、重新调整参数β_k为β_k+1，β_k+1＝β_k+Δβ，Δβ根据各变量的参数范围进行选择，然后执行步骤S5；

S11、结束迭代，获得多个局部最小确认值；

S12、判断多个局部最小确认值是否为唯一解；若是，则执行步骤S13；若否，则执行步骤S14；

S13、该局部最小确认值为最优解β_opt；

S14、计算F_k(β)，获得F_k(β)最小时对应的局部最小确认值，该局部最小确认值为最优解β_opt。