基于定性分析与定量预测的金属材料大气腐蚀预测方法与流程

本发明涉及一种定性分析与定量预测有机统一的金属材料大气腐蚀预测方法，属于环境适应性评价
技术领域：
。
背景技术：
：金属材料大气腐蚀是一个普遍而严重的问题，容易引起各型装备失效，带来巨大经济损失和社会危害。反映金属材料真实大气环境适应性的自然环境试验某些情况下耗时很长，可能数年，甚至数十年，虽然所获数据十分宝贵却很难跟上材料研究和装备研制的高速发展；而实验室加速腐蚀试验与装备实际服役工况存在不可避免的差异。由于金属材料自然环境试验的时间跨度长，空间分布广，腐蚀影响因素复杂多变，其大气腐蚀数据通常具有小样本、数据分布蕴含层次结构、高噪声等特征，这一特点制约了很多传统数据分析技术的应用，特别是大数据分析技术的应用，增加了金属材料大气腐蚀准确预测的困难性。目前，国内外是通过部分或全部原始腐蚀数据训练建立幂函数、指数函数、神经网络、多元线性回归等模型后进行金属材料大气腐蚀预测。然而，金属材料大气腐蚀通常初期发展迅速，随着暴露时间的延长，金属材料大气腐蚀有不断减小的趋势。金属材料大气腐蚀原始数据只能反映腐蚀阶段性的历史状态，这个状态与金属未来的腐蚀发展状态可能存在偏差，容易导致定量预测结果与直观定性分析结果大相径庭。此时我们还需要根据经验，将金属材料大气腐蚀未来的发展趋势施加到现有模型中，从而实现未来更准确的预测。此外，基于既有腐蚀原始数据建立数学模型虽然常常获得较好的模拟精度，但是这并不表示我们已经掌握了金属材料大气腐蚀演变的客观规律，在实际腐蚀预测过程中经常出现拟合效果很好但是对腐蚀发展趋势预测精度不高的情况。因此，提出一种定性分析与定量预测有机统一的金属材料大气腐蚀预测方法，将金属长期腐蚀的发展趋势有机融入腐蚀模型构建过程中，提高金属材料大气腐蚀的预测准确度和可靠性是非常必要且具有重要意义。技术实现要素：针对上述现有技术中的不足之处，本发明旨在提供一种基于定性分析与定量预测的金属材料大气腐蚀预测方法，通过这样一种将定性分析结果定量化的“跨界”方法，能够克服现有金属材料大气腐蚀预测方法建模时完全依赖原始数据的不足，避免“就模型而模型”的预测误区，提高金属材料大气腐蚀预测结果的准确度和可靠性。为实现上述目的，本发明的技术方案：一种基于定性分析与定量预测的金属材料大气腐蚀预测方法，其包括如下步骤：(1)通过开展金属材料大气环境试验，获取金属材料大气腐蚀失重原始数据序列w0：w0＝{w0(t1),w0(t2),...,w0(tm)}，式中，w0(t1),w0(t2),…w0(tm)分别代表金属材料在大气环境试验时间为t1、t2、…tm时采集的大气腐蚀失重原始数据，t1、t2、…tm为互不相同的任意正实数，m表示在整个大气环境试验周期内采集的金属材料大气腐蚀失重原始数据个数，m为不小于3的正整数；(2)对金属材料大气腐蚀失重原始数据序列w0进行弱化变速转换，构建大气腐蚀失重弱化变速数据序列ws：ws＝w0d＝{w0(t1)d,w0(t2)d,...,w0(tm)d}，式中，w0(t1)d,w0(t2)d,…w0(tm)d分别代表金属材料在大气环境试验时间为t1、t2、…tm时的大气腐蚀失重弱化变速数据，其中，式中，k＝1,2,……,m，α为权重因子；(3)基于金属材料大气腐蚀失重弱化变速数据序列ws进行幂函数建模，得到金属材料大气腐蚀预测值w(t)＝atn，其中，t代表大气环境试验时间，w(t)代表金属材料在大气环境试验时间为t时的大气腐蚀失重预测数据，a、n是基于金属材料大气腐蚀失重弱化变速数据序列ws进行幂函数拟合建模后得到的常数。进一步还包括步骤：(4)利用自然环境试验得大气腐蚀失重实测数据，并将该实测数据与金属材料大气腐蚀预测结果进行对比验证，计算其相对预测误差，计算公式如下：其中，w(t)表示金属材料在大气环境试验时间为t时的大气腐蚀失重预测数据，wr(t)表示金属材料在大气环境试验时间为t时的大气腐蚀失重实测数据。进一步的，按gjb8893-2017《军用装备自然环境试验方法》开展金属材料大气环境试验，获得金属材料大气腐蚀失重数据。进一步的，所述α取值范围为0～1。进一步的，所述α＝0.5。本发明的技术效果：采用弱化变速转换与幂函数拟合建模相结合的方法，将金属长期腐蚀的发展趋势施加到幂函数模型中，能够实现金属材料大气腐蚀发展定性分析与定量预测的有机统一，克服现有金属材料大气腐蚀预测方法建模时完全依赖原始数据的不足，避免“就模型而模型”的预测误区，提高金属材料大气腐蚀预测结果的准确度和可靠性。此方面简便易行，对于利用金属材料短周期大气腐蚀数据预测中长期大气腐蚀发展趋势，提高金属材料大气腐蚀预测准确度和可靠性具有重要支撑作用。附图说明图1是本发明的流程框图。具体实施方式下面结合具体及附图来进一步详细说明本发明。一种基于定性分析与定量预测的金属材料大气腐蚀预测方法，包括如下步骤：(1)按照gjb8893-2017《军用装备自然环境试验方法》开展金属材料大气环境试验，获取金属材料大气腐蚀失重原始数据序列w0：wo＝{wo(t1),wo(t2),...,wo(tm)}，式中，wo(t1),wo(t2),…wo(tm)分别代表金属材料在大气环境试验时间为t1、t2、…tm时采集的大气腐蚀失重原始数据；t1、t2……tm为互不相同的任意正实数；m表示在整个大气环境试验周期内采集的金属材料大气腐蚀失重原始数据个数，m为不小于3的正整数；这里整个大气环境试验周期为任意正实数，在试验开始时就已设定。(2)对金属材料大气腐蚀失重原始数据序列wo进行弱化变速转换，构建大气腐蚀失重弱化变速数据序列ws：ws＝wod＝{wo(t1)d,wo(t2)d,...,wo(tm)d}，式中，wo(t1)d,wo(t2)d,…wo(tm)d分别代表金属材料在大气环境试验时间为t1、t2、…tm时的大气腐蚀失重弱化变速数据，其中，式中，k＝1,2,……,m，α为权重因子，α优选但不限于为0.5。(3)基于金属材料大气腐蚀失重弱化变速数据序列ws进行幂函数建模，得到金属材料大气腐蚀预测值w(t)＝atn，其中，t代表大气环境试验时间，w(t)代表金属材料在大气环境试验时间为时的大气腐蚀失重预测数据；a、n是基于金属材料大气腐蚀失重弱化变速数据序列ws进行幂函数拟合建模后得到的常数。实际工作中，可用excel或者origin等软件进行幂函数曲线拟合建模。(4)利用自然环境试验的大气腐蚀失重实测数据，并将该实测数据与金属材料大气腐蚀预测结果进行对比验证，计算其相对预测误差，计算公式如下：其中，w(t)表示金属材料在大气环境试验时间为t时的大气腐蚀失重预测数据，wr(t)代表金属材料在大气环境试验时间为t时的大气腐蚀失重实测数据。需要特别说明的是，本方法中，在获取金属材料大气腐蚀失重原始数据序列w0时，时间t1、t2……tm为互不相同的任意时间值，且t1<t2<……<tm，比如：可以是第1年、第2年、第3年，也可以是第2年、第3年、第5年、第9年，也可以是第1个月、第3个月、第4个月、第6个月、第9个月等，但t1、t2……tm所选时间单位需与金属材料大气腐蚀预测值w(t)＝atn中的代表大气环境试验时间t的时间单位一致。为了更好地理解本发明中的金属材料大气腐蚀预测方法，下面采用实例来说明本案中一种定性分析与定量预测有机统一的金属材料大气腐蚀预测方法。实例：一种定性分析与定量预测有机统一的高强钢大气腐蚀预测方法(1)按照gjb8893-2017《军用装备自然环境试验方法》在海南万宁试验站开展试验周期为5年的高强钢大气环境试验，获取高强钢在大气环境试验时间为1年、2年、3年、5年的腐蚀失重原始数据。采用大气环境试验时间为1年、2年、3年的腐蚀失重原始数据为样本预测大气环境试验为5年时的腐蚀失重。第5年腐蚀失重实测数据作为验证数据。高强钢大气环境试验时间为1年、2年、3年时的腐蚀失重原始数据如表1所示。表1高强钢大气环境试验时间为1年、2年、3年时的腐蚀失重原始数据(g/m2)1年2年3年万宁2940.00005628.31807319.3100高强钢大气环境腐蚀失重原始数据序列w0为：wo＝{wo(t1),wo(t2),wo(t3)}＝{2940.0000，5628.3180，7319.3100}(2)对高强钢原始腐蚀失重数据进行弱化变速转换(α＝0.5，m＝3)。wo(3)d＝0.5×7319.3100+0.5×7319.3100＝7319.3100高强钢大气腐蚀失重弱化变速数据序列ws为：ws＝wod＝{wo(t1)d,wo(t2)d,wo(t3)d}＝{5321.964,6782.463，7319.310}(3)基于高强钢大气腐蚀失重弱化变速数据序列ws，利用excel建立幂函数模型为，w(t)＝5375t0.296(r2＝0.979)利用幂函数模型预测高强钢大气环境户外暴露5年的大气腐蚀失重，t＝5时，w(5)＝8655.13(g/m2)(4)高强钢大气环境试验5年的大气腐蚀失重实测数据为wr(5)＝7993.3570(g/m2)根据高强钢大气环境试验5年腐蚀失重实测数据，计算预测相对误差为为对比本方法与常用的单一幂函数方法预测效果差异性，采用表1中的高强钢大气腐蚀失重原始数据序列wo，利用excel建立幂函数模型为，w(t)＝2992t0.841(r2＝0.991)t＝5时，w(5)＝11582.31(g/m2)根据高强钢大气环境试验5年腐蚀失重实测数据，计算其相对误差为通过应用实例可以看出，采用单一幂函数方法进行高强钢大气腐蚀预测时出现拟合效果很好但是预测误差非常大的情况。采用弱化变速转换与幂函数建模相结合的方法，高强钢大气腐蚀预测误差从44.90％降为8.28％，高强钢大气腐蚀预测准确度和可靠性大为提升。以上对本发明实施例所提供的技术方案进行了详细介绍，本文中应用了具体个例对本发明实施例的原理以及实施方式进行了阐述，以上实施例的说明只适用于帮助理解本发明实施例的原理；同时，对于本领域的一般技术人员，依据本发明实施例，在具体实施方式以及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。当前第1页12