专利名称：：一种预测金属材料大气腐蚀速率的方法
技术领域：
：本发明涉及大气腐蚀的腐蚀预测技术，具体为一种预测金属材料大气腐蚀速率的方法，主要应用于区域性碳钢、低合金钢、铜、铝、锌或其他金属材料的大气腐蚀性预测研
背景技术：
：我国开展大气腐蚀研究三十多年，也做了大量的试验工作，相对于我国辽阔的地域，仍显得不足，特别是区域性腐蚀调查和腐蚀预测方面严重落后。在自然环境下进行暴露腐蚀试验，影响金属材料腐蚀的因素很多。在无污染或环境污染较轻的环境下，主要影响因素有相对温度、相对湿度、降雨量(或降雨时数)、日照时间等。目前，所积累的大量腐蚀数据和环境数据都是分散在全国各个试验点，不但试验点数量少，并且每个试验点的气候环境和污染环境都有显著差别，彼此不能形成整体进行综合分析，就难以进行腐蚀预测研究。
发明内容本发明的目的在于提供一种预测金属材料大气腐蚀速率的方法，解决现有技术中难以进行腐蚀预测研究等问题。本发明的技术方案是—种预测金属材料大气腐蚀速率的方法，步骤如下(1)以金属材料为试验材料，在需要预测地区选择试验点，依据IS0/DIS8565-1987《金属和合金户外大气腐蚀试验的一般规则》进行自然暴露试验，同时监测各试验点的气象因素，试验周期为一年至二年；(2)试验结束后，将试验材料除锈，样品经干燥称重后，计算腐蚀失重，同时统计各个试验点的气象数据；(3)采用多元线性回归分析的统计方法进行回归分析，建立回归方程Y=a+bXi+cX^dXfeXffXs其中，a、b、c、d、e、f为回归系数，Y代表钢的腐蚀失重，X工代表平均温度，X2代表相对湿度，X3代表相对湿度>80%的时数，X4代表日照时数，X5代表降雨量。所述的预测金属材料大气腐蚀速率的方法，在没有环境污染的条件下，影响金属材料腐蚀的主要因素有平均温度、相对湿度、相对湿度>80%的时数、日照时数和表降雨所述的预测金属材料大气腐蚀速率的方法，以该回归方程可以对海南地区或者气候条件与之相近地区的大气腐蚀速率进行预测。所述金属材料为Q235钢、低合金钢、铜、铝、锌或其他金属材料。所述试验材料固定在金属架上，试验材料与金属架通过绝缘处理，防止缝隙腐蚀的发生。所述试验材料面朝正南，与水平面呈45°角，悬挂在距地面3米高处。3所述试验结束后，针对不同的试验材料，采用不同的除锈溶液和试验条件进行清洗腐蚀产物，同时用未腐蚀的钢样校正除锈液对钢基体的腐蚀。本发明的优点及有益效果如下1、金属材料在自然环境下的暴露腐蚀试验，是研究金属大气腐蚀最常用的试验方法，获得的腐蚀数据最准确，能真实体现试验地点的腐蚀状况。由此得出的腐蚀数据是制订相关国家标准的重要依据，也是各项防腐工程设计的理论基础。但是，受到大气自然暴露试验周期很长的制约，不能立即知道某地的腐蚀状况。因此有针对性地选择试验地点进行腐蚀试验，通过对整个区域的腐蚀状况进行综合分析，研究影响金属材料腐蚀的环境因素，并进行回归分析，建立回归方程。通过检测某地的环境因素或者采用历史数据就可以预测某地的腐蚀数据。2、本发明第一次对区域性的大气腐蚀速率进行了预测，并提供了范例。同以往仅能对单一地点的纵向预测进行了扩展，使更大范围的腐蚀预测变得可行。3、本发明通过多元线性回归分析，气象因素中的平均温度、相对湿度、相对湿度大于80X时数、日照时数和降水量五个参数，与Q235钢的腐蚀失重有相关性，在a=0.05的水平上都是显著的。4、一个地区的气象因素通常在一定的范围内波动，该波动范围一般较小。只要知道海南地区某地的气象因素，通过回归方程就可以对该地的Q235钢的腐蚀速率进行预测。多元线性回归分析的原理假设对变量y禾PXl、x2.....xp，收集了n组独立观测值(xal、xa2.....xap，ya)，a=1，2，…，n，(1.1)假定它们之间有如下关系ya=P0+elXal+.+Ppxap+ea，a=1，2，，n，(1.2)其中，|3。，1^，...，~是。+1个待估参数，^，^，...，、是。个可以精确测量或可控制的一般变量的第a组观测值，^，h，...，、是n个相互独立的随机变量，且E(ea)=0，D=(ea)=o2，a=1，2，，n;这就是多元线性回归的数学模型。根据给出的n组观测值(1.l)，采用最小二乘法求出(1.2)中参数P。，！^，...，Pp的估计A，A，...，A;从而，可以获得p元线性回归方程J)=A+Ax,+在x2&xp(1.3)要求在，a，…，》能使00。'A'…，a^=2](j。-A-A、r…-a^)2，<formula>formulaseeoriginaldocumentpage4</formula>其中，_^=在+Axw+…+Ax^经整理，得即《<formula>formulaseeoriginaldocumentpage4</formula>组的解》al/4+lXalA+…+》。1/^=》。'凡'(")从(1.4)中第一式可得将(1.5)代入(1.4)的后p个式子，并经整理后知A，AA是下列方禾'王乂「"i^+/12/2+...+/iPA=/iy'/2lA+/22A+...+/2PA=/2y,(1.6)A后代入(1.5)，可求得在，其中从(1.6)求得A，A，…，》ai,i=1，2，，p，歹,二丄I]少a，X丄i，j=l，2，...，p，(1.7)i=1，2，.，p。(1.8)由上述公式可知，求P元线回归方程中系数的步骤如下1、计算各变量之和;&'，i=4，2，…，P，2、计算各变量之平均值《，i=l，2，...，p，3、计算各变量的平方和及变量间的乘积和》ai》aj，i，j=l，2，…，P;》ai》a，i=l，2，4、按公式计算(1.7)、(1.8)，计算I".,liy，i，j5、解方程组(1.6)，得A，A，…，式；6、按公式(1.5)计算戍。P;1，2，P具体实施方式实施例1:1990年6月1992年5月，在海南地区进行了2年的大气腐蚀性调查。在20个市县设有试验点，基本呈均匀分布。试验材料为Q235钢，样品规格为100mmX50mmX3mm，Q235钢的化学成分见表l。样品经铣床、磨床，加工至光洁度V7或以上，丙酮除油，乙醇脱水，放置干燥器24h后，用感量为0.lmg分析天平称重。表1、A3钢的化学成份(wt%)c<table>tableseeoriginaldocumentpage6</column></row><table><table>tableseeoriginaldocumentpage7</column></row><table>分析各个试验点的气象数据与腐蚀失重的相关关系，采用多元线性回归的统计分析方法，分析影响金属材料腐蚀的主要因素，并建立回归方程。用该回归方程可以对海南地区或者气候条件相近地区Q235钢的腐蚀速率进行腐蚀预测。气象因素对Q235钢腐蚀失重(第一年)影响的回归过程1、计算各变量的平均值(见表3)表3、所有变量的平均值<table>tableseeoriginaldocumentpage7</column></row><table>第一列为回归方程中所有变量的序号，第二列为各个变量的平均值。L(矩阵)表4、回归系数矩阵<table>tableseeoriginaldocumentpage8</column></row><table>INVERSOFL(逆矩阵)表5、回归系数的逆矩阵<table>tableseeoriginaldocumentpage8</column></row><table>回归系数B丄=5.4510;B2=2.8536;B3=0.0402;B4=-0.00701;B5=-0.01839;B0=-293.6729。回归方程:Y=-293.6729+5.4510X^2.8536X2+0.0402X3_0.00701X4_0.01839X5总的平方和SO=3026.871，回归平方和SR=2861.22，剩余平方和SE=165.6506;复相关系数R=0.97225，S估计SGM=3.5696;对方程作检验的F比F=44.9088;FI(l)=13.17865>F001(l，13)=9.07;FI(2)=54.4332>F0.01(1，13)=9.07;FI(3)=9.5393>F0.01(1，13)=9.07;FI(4)=5.1022>F,(l，13)=4.67;FI(5)=45.7070>F0.01(1，13)=9.07;五个变量在a=0.05的水平上都是显著的。表6腐蚀失重的回归值和拟合偏差<table>tableseeoriginaldocumentpage9</column></row><table>本实施例中，腐蚀预测原理及相关数据如下给出Xl、x2.....Xp的一组值xQ1、x。2.....x加，可以求出对应的E(y。)的估计值<formula>formulaseeoriginaldocumentpage10</formula>当n较大，x。j与^(j=1，2，...，P)较接近时，可用下式求得j)。的预测区间(近似的)95%的预测区间为(夕。-2》，j)。+2纟)；99%的预测区间为(》。-3纟，j)。+3纟)；其中，腐蚀预测相关数据如下上述回归方程的标准差》=3.5696;现测得琼海大气腐蚀实验站的气象数据如下年平均温度24.5°C;相对湿度86%;相对湿度>80时数6314小时；年降雨量1881mm;年日照时数2116小时。由回归方程求得琼海Q235钢腐蚀失重的估计值是j)=58.936295%的预测区间为(58.9362-2X3.5696,58.9632+2X3.5696)=(51.7970，66.0754);实测琼海实验站Q235钢的腐蚀失重是65.05，与预测值十分接近。注气象数据的检测点与腐蚀实验站点不是同一地点。实施例结果表明，本发明采用多元线性回归分析方法，分析影响Q235钢腐蚀的主要因素，并建立回归方程，以该回归方程对海南地区Q235钢的腐蚀速率进行预测，具有重现性。采用该方法，可以预测各种钢在各种区域环境中的大气腐蚀。10权利要求一种预测金属材料大气腐蚀速率的方法，其特征在于，步骤如下(1)以金属材料为试验材料，在需要预测地区选择试验点，依据ISO/DIS8565-1987《金属和合金户外大气腐蚀试验的一般规则》进行自然暴露试验，同时监测各试验点的气象因素，试验周期为一年至二年；(2)试验结束后，将试验材料除锈，样品经干燥称重后，计算腐蚀失重，同时统计各个试验点的气象数据；(3)采用多元线性回归分析的统计方法进行回归分析，建立回归方程Y＝a+bX1+cX2+dX3-eX4-fX5其中，a、b、c、d、e、f为回归系数，Y代表钢的腐蚀失重，X1代表平均温度，X2代表相对湿度，X3代表相对湿度＞80％的时数，X4代表日照时数，X5代表降雨量。2.按照权利要求1所述的预测金属材料大气腐蚀速率的方法，其特征在于，在没有环境污染的条件下，影响金属材料腐蚀的主要因素有平均温度、相对湿度、相对湿度>80%的时数、日照时数和表降雨量。3.按照权利要求1所述的预测金属材料大气腐蚀速率的方法，其特征在于，以该回归方程可以对海南地区或者气候条件与之相近地区的大气腐蚀速率进行预测。4.按照权利要求1所述的预测金属材料大气腐蚀速率的方法，其特征在于，所述金属材料为Q235钢、低合金钢、铜、铝、锌或其他金属材料。5.按照权利要求1所述的预测金属材料大气腐蚀速率的方法，其特征在于，所述试验材料固定在金属架上，试验材料与金属架通过绝缘处理，防止缝隙腐蚀的发生。6.按照权利要求1所述的预测金属材料大气腐蚀速率的方法，其特征在于，所述试验材料面朝正南，与水平面呈45°角，悬挂在距地面3米高处。7.按照权利要求1所述的预测金属材料大气腐蚀速率的方法，其特征在于，所述试验结束后，针对不同的试验材料，采用不同的除锈溶液和试验条件进行清洗腐蚀产物，同时用未腐蚀的钢样校正除锈液对钢基体的腐蚀。全文摘要本发明涉及大气腐蚀的腐蚀预测技术，具体为一种预测金属材料大气腐蚀速率的方法，主要应用于区域性碳钢的大气腐蚀性预测研究。以金属材料为试验材料，在需要预测地区选择试验点，依据ISO/DIS8565-1987《金属和合金户外大气腐蚀试验的一般规则》进行自然暴露试验，同时监测各试验点的气象因素，试验周期为一年至二年；试验结束后，将试验材料除锈，样品经干燥称重后，计算腐蚀失重，同时统计各个试验点的气象数据；采用多元线性回归分析的统计方法进行回归分析，建立回归方程，以该回归方程对需要预测地区的大气腐蚀速率进行预测。本发明具有重现性，采用该方法可以预测各种钢在各种区域环境中的大气腐蚀。文档编号G01N5/04GK101782499SQ20091001024公开日2010年7月21日申请日期2009年1月21日优先权日2009年1月21日发明者于国才,王振尧,陈鸿川,韩薇申请人:中国科学院金属研究所