发明内容

(1)本发明的目的：

本发明的目的是提供一种考虑腐蚀动力学稳定性的大气腐蚀预测方法，即一种考虑腐蚀动力学稳定性的大气腐蚀模型预测方法，它是针对大多数大气腐蚀预测时没有充分考虑环境因素、缺少环境因素分析的讨论、缺少模型预测的稳健性考虑的一系列问题，提供一种考虑腐蚀动力学稳定性的大气腐蚀预测方法，它是一种包含大气环境因素敏感性分析、多因素大气腐蚀模型建立、模型稳定性准则确定、模型参数估计和腐蚀损失预测的大气腐蚀预测方法。通过考虑到不同区域的不同环境因素差异性，构建多因素大气腐蚀预测模型，基于极大似然估计理论开展模型参数估计，考虑腐蚀动力学稳定性对模型的稳健性进行讨论，最终给出大气腐蚀预测方法。

(2)技术方案：基于上述理论和思路，本发明提供一种考虑腐蚀动力学稳定性的大气腐蚀预测方法，即一种考虑腐蚀动力学稳定性的大气腐蚀模型预测方法，其具体的实施步骤如下：

步骤一：大气环境因素敏感性分析

首先，根据收集的大气腐蚀外场试验数据和环境数据，采用灰色关联度对环境因素敏感性进行分析；例如，大气环境中某金属实测的平均腐蚀速率和最大蚀坑深度为参考数列，则有参考数列为：

X₀(k)＝{x₀(1),x₀(2),…,x₀(n)}, (1)

这里腐蚀数据量为n，X₀(k)代表大气环境中的平均腐蚀速率或腐蚀损失量，k＝1,2,…,n；则有比较数列为：

X_i(k)＝{x_i(1),x_i(2),…,x_i(n)}. (2)

这里环境因素种类为m，X_i(k)(i＝1,2,…,m)代表大气环境因素(温度、相对湿度、降雨量)；因此，参考数列X₀(k)与比较数列X_i(k)的关联系数ξ_0i(k)可表示为：

这里记ρ为分辨系数，分辨系数的值越小，灰色相关系数的范围就越大，但这不会影响比较序列的最终优先级，因此一般根据信息最小的原则取值为0.5；

通过灰色关联系数ξ_0i(k)，可得到：

这里γ_0i(X₀,X_i)为参考数列X₀(k)与比较数列X_i(k)的灰色关联度；对于参考数列X₀(k)和比较数列X_i(k)，关联度可直接表示为γ_i，代表参考数列和比较数列之间的相关程度，则从大到小进行排序则可得到环境因素的影响排序；通常来说当γ_i＞0.6，即认为该序列具有较好的关联性；

步骤二：多因素大气腐蚀模型建立

在腐蚀预测模型中，腐蚀损失是一般是时间的函数；最广泛使用的腐蚀预测模型为幂律模型：

C＝utⁿ， (5)

其中C表示单位面积腐蚀造成的质量损失(μm)，t是腐蚀时间(年)；参数u和质量损失时间指数n通过实测数据进行估计，并且第一年的质量损失用参数u来表示；

通过步骤一的大气环境因素敏感性分析，可以得到对腐蚀损失影响最显著的环境因素排序，在幂律模型的基础上，根据环境因素的影响大小选择最相关的多个环境因素进行建模：

y＝α·t^β·f(δ), (6)

其中y表示腐蚀损失，α,β为参数；在该模型中，f(δ)视为环境因素的函数，代替公式(5)中的质量损失时间指数n表征的环境因素对腐蚀造成的影响，即此时参数β不再视为表征该模型中腐蚀产物防护性能的系数；

步骤三：模型稳定性准则确定

传统的腐蚀预测模型的确定方法，一般考虑模型对于数据的拟合优度情况，然而拟合得好并不能保证预测得好；因此，本发明提出一种新的模型确定准则，过程如下：

①模型稳定性分析

考虑到数据分散性的影响，模型(6)可以表示为:

Y(t)＝α·tβ·f(δ)+ε,ε～N(0,σ²), (7)

这里对模型(6)增加误差项ε，这个误差项服从均值为0和标准差为σ的正态分布；因此，材料腐蚀寿命的累积分布函数F(t)可以表示为：

F(t)＝P(Y(t)≤γ)＝Φ((γ-α·t^β·f(δ))/σ²), (8)

这里γ为给定的腐蚀损失失效阈值，Φ(·)表示标准正态分布的累积分布函数(CDF)；

事实上，腐蚀过程会导致金属材料性能的退化，而对于退化过程中的幂律模型可以等价于腐蚀动力学中的基础的幂律模型，因此，该模型看作是一个腐蚀等效模型；考虑退化过程的机理等同性，则有:

其中当环境应力水平分别为δ_i和δ_j时，对于给定的可靠度R，这里材料的可靠寿命分别为t_R,i和t_R,j；因此，可靠寿命之比为：

K_ij＝t_R,j/t_R,i, (10)

K_ij定义为在第i次环境应力水平和j次环境应力水平下的腐蚀当量系数，其应当只与环境应力水平有关；

根据公式(9)和公式(10)可得：

因为K_ij只与环境应力水平有关，与工作时间和可靠性无关，则可以得到：

σ_i＝σ_j,β_i＝β_j, (12)

该条件为模型的腐蚀动力学稳定性；

②模型稳定性的确定准则

如①所述，通过公式(12)可知，在腐蚀动力学稳定性的条件下，形如公式(6)的腐蚀模型应该有一个稳定的参数β；当腐蚀动力学保持稳定时，腐蚀当量系数应当是由环境应力水平决定的常数，与其他参数无关；然而，在样本量小的情况下，在对试验参数β进行估计时会有误差，为了模型的预测更加稳定，这个误差的差异性应当尽可能地小，即好的模型β的方差应更加稳定；因此，考虑通过变异系数来表示模型的稳定性大小，即：

CV_i＝σ_βi/μ_βi, (13)

这里σ_βi为参数β的标准差，μ_βi为参数β的均值；

步骤四：模型参数估计

在给出了模型的确定准则后，需要对模型进行参数估计；为了简化参数估计的过程，考虑到数据的分散性的同时对模型(6)对数化处理，则有：

lny_i＝g(lnδ)_i+ε₁, (14)

这里ε₁～N(0,σ₁ ²)，然后通过极大似然估计得出参数β的估计值，为了进行比较，通过极大似然函数与Fisher信息矩阵可以求出参数β的均值与方差，从而得出参数β的变异系数大小；其中，Fisher信息矩阵是指通过极大似然估计后得到的参数β的Fisher信息量的矢量化定义；

步骤五：腐蚀损失预测

根据上述步骤确定腐蚀损失模型，并通过4年的腐蚀损失数据以及环境数据得到腐蚀模型参数，从而进行8年腐蚀损失的预测；

通过以上步骤，本发明提出了考虑多种环境因素的大气腐蚀模型，与大多数模型相比，该模型不仅参考了腐蚀预测时间，还重点关注了与腐蚀损失相关的环境因素变化；在考虑腐蚀动力学稳定性的基础上，提出了一种新的腐蚀模型选择准则，有助于寻找更合理的腐蚀预测模型，本方法计算简便，容易实现，且更加符合工程实际。

(3)优点和功效：

本发明是一种考虑腐蚀动力学稳定性的大气腐蚀预测方法，其优点为：

①通过分析影响大气环境的主要因素，提出了考虑多种环境因素的大气腐蚀模型，与大多数模型相比，该模型的优点在于它不仅参考了腐蚀预测时间，还重点关注了与腐蚀损失相关的环境因素变化，这些因素会影响材料和部件的寿命。在考虑腐蚀动力学稳定性的基础上，提出了一种新的腐蚀模型选择准则。利用极大似然估计方法可以计算出参数的变异系数，给出模型的稳健性比较，有助于寻找更合理的腐蚀预测模型，因为对于一般形式的腐蚀模型，在腐蚀动力学稳定性的条件下，参数应该有更稳健的结果，这意味着其估计值之间的差异应该较小。

②本发明提出的方法计算简便，容易实现，且更加符合工程实际，方便工程技术人员掌握使用，方法科学，工艺性好，便于应用推广。