蘇三慶，高 波，王 威，易術(shù)春，馬小平

(西安建筑科技大學(xué)土木工程學(xué)院，陜西 西安 710055)

鋼材在建筑鋼結(jié)構(gòu)、橋梁鋼結(jié)構(gòu)及大跨空間鋼結(jié)構(gòu)等土木工程領(lǐng)域得到廣泛使用，鋼結(jié)構(gòu)的節(jié)點等關(guān)鍵部位在服役過程中承受各種動載或交變荷載作用，容易造成局部應(yīng)力集中和微觀缺陷拓展，這可能導(dǎo)致構(gòu)件在低于材料屈服強度情況下發(fā)生突發(fā)性的斷裂，及時找出結(jié)構(gòu)斷裂前的危險區(qū)域?qū)@些大型工程項目尤為重要．目前，常規(guī)的無損檢測方法存在只能檢測出已成形的宏觀缺陷的弊端，在實際工程應(yīng)用中具有一定的局限性；而金屬磁記憶檢測技術(shù)作為一門新興的無損檢測方法[1-2]，它不需要外磁場的激勵，而以地磁場為背景磁場，通過測量應(yīng)力引起的殘余漏磁場，即可對鐵磁材料的應(yīng)力集中程度和損傷等實現(xiàn)早期檢測診斷，解決了上述問題，所以，金屬磁記憶檢測技術(shù)有望在建筑鋼結(jié)構(gòu)等工程中的早期健康檢測得到廣泛應(yīng)用[3-5]．

現(xiàn)階段金屬磁記憶檢測方法的研究主要集中在探索性實驗和理論模型建立等方面，例如，任吉林等對 20#鋼進(jìn)行了不同程度的拉伸試驗，測量了卸載后的試件表面漏磁場信號，指出磁記憶信號的變化與材料內(nèi)部微觀組織的狀態(tài)密切相關(guān)[6]；基于經(jīng)典Jiles-Artherton理論模型，冷建成等對Q235低碳鋼做了拉伸試驗，提出了考慮塑性變形影響的改進(jìn)模型[7]，通過對45#鋼扭轉(zhuǎn)試驗，徐明秀等驗證了扭轉(zhuǎn)材料的改進(jìn)Jiles-Atherton模型，拓展了該模型的應(yīng)用范圍[8]；楊理踐等采用密度泛函理論的平面波贗勢法，計算了晶體磁性能隨力的定量變化關(guān)系，并對金屬磁記憶效應(yīng)中的力-磁耦合關(guān)系做出定量分析[9]．但是目前缺乏相關(guān)數(shù)值模擬的研究，姚凱[10]等假設(shè)鐵磁材料的應(yīng)力和磁導(dǎo)率之間滿足線性關(guān)系，僅分析了一定應(yīng)力水平下提離值等對漏磁信號幅值的變化趨勢，未做定量化分析．

因此，通過對建筑鋼結(jié)構(gòu)廣泛使用的Q345B鋼的人工淺槽拉伸試件進(jìn)行二維平面有限元數(shù)值模擬，根據(jù)應(yīng)力表征鐵磁材料磁參數(shù)(磁導(dǎo)率)的力磁效應(yīng)關(guān)系式，分析不同應(yīng)力狀態(tài)下空氣中法向漏磁信號B(y)及其梯度的分布規(guī)律，最后通過試驗驗證模擬結(jié)果．

1 基本原理

1.1 麥克斯韋方程組

麥克斯韋方程組是經(jīng)典電磁理論的基礎(chǔ)，它可以解決任何宏觀電磁現(xiàn)象和電磁過程．金屬磁記憶檢測是以地磁場為背景磁場，對給定地理位置座標(biāo)的空間來說，可以認(rèn)為地磁場強度值保持恒定不變．當(dāng)鐵磁構(gòu)件處于無電流源域時，其無源穩(wěn)態(tài)磁場的麥克斯韋方程組為

式中：{B}為磁感應(yīng)強度矢量，{H}磁場強度矢量．

為提高電磁場數(shù)值計算的效率，引入矢量磁位函數(shù)(磁矢勢){A}來減少求解的自由度數(shù)，此磁矢勢滿足

從上述(1)～(3)式可推得:

對二維平面模型，對于磁矢勢{A}有

式中，AZ表示磁矢勢{A}在Z方向的分量，它是電磁方程計算中的唯一變量．

1.2 力磁關(guān)系式

文獻(xiàn)[11]根據(jù)磁介質(zhì)磁化理論分析計算，得到了反映鐵磁材料力磁效應(yīng)的函數(shù)關(guān)系式

從該方程可以推導(dǎo)如下方程式

式中：真空磁導(dǎo)率μ0=4π×10-7H/m，Q345B為低碳鋼，其飽和磁感應(yīng)強度Bm=2.5 T，飽和磁滯伸縮系數(shù)λm=5×10-6，初始相對磁導(dǎo)率μi=260，μσ為應(yīng)力不為零時計算得到的相對磁導(dǎo)率．

2 模擬過程

模擬由靜力學(xué)分析和靜磁學(xué)分析兩階段組成．靜力學(xué)采用線性強化彈塑性模型分析得到鋼板拉伸荷載下的應(yīng)力分布結(jié)果，將此結(jié)果保存在單元表中；在靜磁學(xué)分析中，調(diào)用儲存在單元表中的應(yīng)力分布結(jié)果來修改鐵磁材料的相對磁導(dǎo)率，以此來模擬應(yīng)力對鐵磁材料磁學(xué)屬性的影響[12]．

圖1給出了Q345B鋼試件的尺寸模型，在試件中部的兩側(cè)預(yù)先設(shè)置對稱分布的人工淺槽以模擬應(yīng)力集中，淺槽大小為l×b=1 mm×3 mm，試件厚δ=10 mm，檢測線分布在槽口附近，長度為60 mm，檢測點等間距分布．

圖1 試件尺寸Fig.1 Size of specimen

2.1 建立模型

靜力學(xué)和靜磁學(xué)分析時，都采用相同單元劃分的試件尺寸模型，為考慮應(yīng)力集中的影響，對缺陷所在區(qū)域的網(wǎng)格細(xì)化．靜磁學(xué)分析時還應(yīng)在試件模型外部添加空氣模型，以采集試件外表面空氣中的漏磁信號．建立的試件尺寸模型如圖2所示．

圖2 單元網(wǎng)格劃分Fig.2 Meshing of finite element model

劃分網(wǎng)格的單元大小設(shè)置為 1 mm，在缺陷處包含一層網(wǎng)格加密區(qū)，圖2(a)為外層空氣模型，圖2(b)為內(nèi)層試件模型，單元總計15 154個．

2.2 加載和求解

靜力學(xué)分析中，對模型一端施加全部位移約束，另一端施加均勻拉力的線載荷，建好模型后，端面的約束和載荷會自動傳遞到單元的節(jié)點上．為了使模擬結(jié)果和實驗結(jié)果具有可比性，荷載的加載等級與實驗結(jié)果一致．

在地磁場環(huán)境下進(jìn)行二維平面靜磁學(xué)分析時，除了試件模型外，其外側(cè)還有一層空氣模型，對于空氣模型的外邊界需施加通量垂直邊界條件[13]．根據(jù)力磁關(guān)系式(8)可求得不同應(yīng)力水平下鐵磁材料的磁導(dǎo)率，鐵磁材料單元的磁學(xué)屬性賦值以后，便可得到該模型在地磁場下施加拉力載荷后鐵磁材料的表面漏磁場分布情況，并將計算所得的磁場強度法向分量保存起來．因為計算的單元數(shù)量很多，模擬時通過ANSYS參數(shù)化語言將材料的磁學(xué)屬性賦值給單元，并保存計算得到的漏磁分布結(jié)果．

3 模擬結(jié)果分析

3.1 應(yīng)力的分布規(guī)律

試驗前對Q345B鋼板拉伸測試發(fā)現(xiàn)，當(dāng)試件拉伸到110 kN時，在試件中部出現(xiàn)非常明顯的頸縮現(xiàn)象，此時頸縮對試件內(nèi)部磁疇結(jié)構(gòu)的影響加劇[6]，不適于用應(yīng)力表征磁信號的變化規(guī)律，因此不考慮110 kN以后試件拉伸作用下磁信號的變化情況．

圖3(a)展示了試件在85 kN拉伸荷載作用下的應(yīng)力云圖，從中可以看出，應(yīng)力幅值的最大值出現(xiàn)在缺口的內(nèi)側(cè)，缺口外邊緣處應(yīng)力幅值最小，并且缺口處應(yīng)力云圖呈現(xiàn)“魚尾”狀分布．圖3(b)給出了試件名義屈服前后，沿著試件長度方向缺口部位的應(yīng)力分布情況，越靠近缺陷位置，應(yīng)力幅值越大，應(yīng)力集中程度越高，圖中出現(xiàn)的雙峰現(xiàn)象與應(yīng)力云圖分布規(guī)律吻合，X坐標(biāo)軸中間位置的應(yīng)力幅值與計算得到的缺陷處名義應(yīng)力值一致．

圖3 槽口處應(yīng)力集中分布Fig.3 Stress concentration of pits in specimen

3.2 磁信號的分布規(guī)律

由金屬磁記憶檢測法的原理[14]可知，在軸向拉伸應(yīng)力作用下鐵磁體內(nèi)部的磁疇組織發(fā)生了定向和不可逆的重新取向，從而導(dǎo)致應(yīng)力或應(yīng)變集中區(qū)域表面的漏磁場發(fā)生變化．當(dāng)拉應(yīng)力增加時，鐵磁性材料中磁疇的磁矩方向趨向一致，最終在應(yīng)力或應(yīng)變集中處積聚大量磁化能形成較強的漏磁場，表現(xiàn)為法向分量B(y)過零點，梯度值K出現(xiàn)極值．圖4分別表示彈、塑性階段應(yīng)力集中區(qū)磁感應(yīng)強度B(y)和其梯度值K隨應(yīng)力的變化規(guī)律．

圖4 有限元模擬的彈、塑性階段法向漏磁信號及其梯度值的分布Fig.4 Normal component with its gradient of MFL value at elastic & plastic stage of the modeling

分析圖 4(a)和圖 4(c)可以看出，宏觀上，彈性階段試件表面的磁信號曲線形貌相近，其大小分布與載荷強度存在正相關(guān)關(guān)系，與磁信號法向分量的分布圖過坐標(biāo)軸原點呈中心對稱不同，梯度值左右兩側(cè)的分布關(guān)于x軸的中線呈軸對稱。細(xì)觀上，進(jìn)一步觀察可見，磁信號法向分量的絕對值隨載荷增加而增大，其大小在應(yīng)力σ=333.3 MPa時比σ=166.7 MPa時增大近一倍，達(dá)到0.4 mT。同時，梯度的幅值也隨載荷增加而增大，在缺陷處梯度值的變化幅度明顯，在x軸中點處出現(xiàn)極值，為?0.04 mT/mm。此外，在試件中部即x=30 mm處，磁信號法向分量始終過零點，不受載荷大小的變化而發(fā)生位置改變，這些與理論分析的結(jié)果符合．進(jìn)入塑性階段后，試件缺陷處名義應(yīng)力從σ=354.2 MPa增大到σ=416.7 MPa，從圖4(b)和圖4(d)中可以看出，應(yīng)力集中區(qū)域的磁信號的變化比彈性階段更劇烈，磁信號法向分量及其梯度值的曲線發(fā)展趨勢保持一致，但幅值變化顯著．需要指出的是，塑性階段缺陷附近處磁信號法向分量及其梯度值的分布范圍沿試件長度方向伸長，這與彈性階段變化稍有區(qū)別，因為塑變對試件的磁疇的影響加劇，造成試件表面的漏磁場強度在缺陷處出現(xiàn)波動，非線性變化趨勢增強．

基于上述分析可知，磁感應(yīng)強度變化較好地反映了彈性階段磁記憶信號隨應(yīng)力集中程度的分布規(guī)律，在試件屈服后，塑變引起的磁疇結(jié)構(gòu)變化顯著，雖然對用應(yīng)力表征磁記憶信號造成干擾，但可利用彈性階段和塑性階段磁信號變化迥異的特點，實現(xiàn)應(yīng)力集中程度及損傷狀態(tài)的早期識別．

4 試驗結(jié)果驗證

試驗材料選用廣泛應(yīng)用的Q345B鋼，考慮退火處理后試件內(nèi)部組織結(jié)構(gòu)的變化會造成試件磁導(dǎo)率相應(yīng)的改變，以及考慮到實際工程情況，因此未對試件進(jìn)行高溫退火處理．試驗中采用在線測量方法，實驗前在試件的檢測范圍內(nèi)畫好檢測線，檢測點間距為5 mm，采集數(shù)據(jù)時，沿著畫好的檢測線依次逐點在線測量試件表面漏磁場法向分量B(y)值，保持探頭垂直緊貼于試件表面，測量時應(yīng)遠(yuǎn)離強磁場的干擾．拉伸試驗在WAW-2000D電液伺服萬能試驗機上進(jìn)行，磁記憶信號檢測設(shè)備選用高精度WT10B 型數(shù)字高斯計，最大量程為200 mT，分辨率為0.01 mT．

對淺槽試件進(jìn)行軸向拉伸測試時，發(fā)現(xiàn)淺槽試件分別加載到85 kN 和110 kN時發(fā)生屈服和出現(xiàn)明顯的頸縮現(xiàn)象．由于頸縮會干擾應(yīng)力和磁記憶信號間的變化規(guī)律，因此對頸縮不予考慮，細(xì)化后淺槽試件的加載等級為：0、20、40、60、80、85、90、100 kN．

圖5(a)和(c)、圖5(b)和(d)分別表示拉伸荷載下試件彈性、塑性階段應(yīng)力集中程度對法向磁感應(yīng)強度B(y)和其梯度值K的影響規(guī)律[5]．

圖5 試驗的彈、塑性階段法向漏磁信號及其梯度值的分布Fig.5 Normal component with its gradient of MFL value at elastic & plastic stage of the experiment

由圖5可知，隨著拉應(yīng)力的增加，磁場強度的絕對值隨之增大，磁場強度的梯度值也隨之增大，在試件中間位置磁場強度值B(y)過零點，磁場強度的梯度值K出現(xiàn)最大值。模擬得到的結(jié)果在試件缺陷處有明顯的上下波動現(xiàn)象，而在試驗中未觀察到此現(xiàn)象，一方面與試驗時信號采集點密度以及儀器分辨率等因素有關(guān)，另一方面試驗時試件兩側(cè)端部夾持產(chǎn)生的漏磁場也會影響檢測準(zhǔn)確度，這兩方面是造成試驗中捕捉到的磁記憶信號與模擬的結(jié)果在缺陷處略有區(qū)別的主要因素．總體上，模擬的計算結(jié)果和試驗數(shù)據(jù)的變化趨勢相符合，漏磁場分布區(qū)域及幅值范圍相差不大，證實了模擬結(jié)果的可靠性，本文有限元模型可用于建筑中鋼板件的拉伸磁記憶有限元分析．

5 結(jié)論

(1) 根據(jù)已知的力磁效應(yīng)表征函數(shù)，對 Q345B建筑鋼板件進(jìn)行了有限元建模分析，得出不同受力階段應(yīng)力集中對鐵磁構(gòu)件磁記憶信號的變化規(guī)律，可為工程早期損傷的識別提供參考．

(2) 有限元模擬結(jié)果較好反映了彈性階段磁信號與應(yīng)力間的變化規(guī)律，進(jìn)入塑性階段以后，應(yīng)變的較快增長對用應(yīng)力來表征磁信號變化造成干擾，但可根據(jù)彈塑性階段磁信號變化迥異的特點，實現(xiàn)損傷早期發(fā)現(xiàn)．

(3) 通過模擬得到的漏磁場分布區(qū)域、幅值分布范圍與實驗結(jié)果相吻合，該數(shù)值模擬方法可用于建筑中鋼板件的拉伸磁記憶有限元分析．

References

[1] DUBOV A A. Development of a metal magnetic memory method[J]. Chemical & Petroleum Engineering, 2012,47:837-839.

[2] WU Dabo, XU Minqiang, LI Jianwei. Study on physical mechanism of metal magnetic memory technique[J].Applied Mechanics & Materials, 2010, 34-35:841-844.

[3] WANG Wei, YI Shuchun, SU Sanqing. Experimental investigation of stress and damage characterization of steel beam buckling using magnetic memory signals[J].Structural Design of Tall & Special Buildings, 2016,25(11):508-518.

[4] 王威, 曾發(fā)榮, 蘇三慶, 等. 基于磁記憶的受彎鋼梁力-磁效應(yīng)試驗[J].材料科學(xué)與工程學(xué)報,2016,34(1):109-114.WANG Wei, ZENG Farong, SU Sanqing, et al.Experimental research on magnetic mechanical effectiveness of bending steel beam based on metal magnetic memory[J]. Materials Science & Engineering,2016,34 (1):109-114.

[5] 王威, 易術(shù)春, 蘇三慶, 等. 建筑鋼試件拉伸應(yīng)力與其磁記憶效應(yīng)漏磁場梯度的對應(yīng)關(guān)系研究[J].西安建筑科技大學(xué)學(xué)報(自然科學(xué)版),2015,47(3):341-346.WANG Wei, YI Shuchun, SU Sanqing, et al. Research on the relationship between magnetic flux leakage gradient by metal magnetic memory effect and tensile stress of building steel specimen[J]. J. Xi?an Univ. of Arch. & Tech.(Natural Science Edition), 2015,47(3):341-346.

[6] 任吉林, 王東升, 宋凱, 等. 應(yīng)力狀態(tài)對磁記憶信號的影響[J]. 航空學(xué)報,2007, 28(3): 724-728.REN Jilin, WANG Dongsheng, SONG Kai, et al.Influence of stress state on magnetic memory signal[J].Acta Aeron Nautical Et Astronautica Sinica, 2007,28(3):724-728.

[7] LENG Jiancheng, LIU Yang, ZHOU Guoqiang, et al.Metal magnetic memory signal response to plastic deformation of low carbon steel[J]. Ndt & E International,2013, 55(3):42-46.

[8] XU Mingxiu, XU Minqiang, LI Jianwei et al. Metal magnetic memory field characterization at early fatigue damage based on modified Jiles-Atherton model[J].Journal of Central South University, 2012,19(6):1488-1496.

[9] 楊理踐, 劉斌, 高松巍,等. 金屬磁記憶效應(yīng)的第一性原理計算與實驗研究[J].物理學(xué)報,2013,62(8):086201-086207.YANG Lijian, LIU Bin, GAO Songwei, et al.First-principles calculation and experimental study of metal magnetic memory effects[J]. Acta Physica Sinica,2013, 62(8): 086201-086207.

[10] 姚凱,王正道,鄧博,等. 金屬磁記憶技術(shù)的數(shù)值研究[J].工程力學(xué),2011(9):218-222.YAO Kai, WANG Zhengdao, DENG Bo, et al. Numerical study on metal magnetic memory technique[J]. Engineering Mechanics,2011,28(9):218-222.

[11] 王社良, 王威, 蘇三慶, 等. 鐵磁材料相對磁導(dǎo)率變化與應(yīng)力關(guān)系的磁力學(xué)模型[J].西安科技大學(xué)學(xué)報,2005,25(3):288-291.WANG Sheliang, WANG Wei, SU Sanqing, et al.Magneto-mechanical model on differential permeability and stress of ferromagnetic material[J].Journal of Xi?an University of Science and Technology, 2005,25(3):288-291.

[12] 閻照文.ANSYS10.0工程電磁分析技術(shù)與實例詳解[M].北京:中國水利水電出版社,2006.YAN Zhaowen.ANSYS10.0Engineering electromagnetic analysis technique and detailed examples[M]. Beijing:China Water Resources and Hydropower Publishing House, 2006.

[13] 胡仁喜. ANSYS13.0多物理耦合場有限元分析從入門到精通[M]. 北京:機械工業(yè)出版社,2012.HU Renxi. ANSYS 13.0 multi-physics coupling finite element analysis (From entry to the master)[M]. Beijing:Mechanical Industry Press, 2012.

[14] 任吉林, 林俊明, 池永濱, 等. 金屬磁記憶檢測技術(shù)[M]. 北京: 中國電力出版社,2000.REN Jilin, LIN Junming, CHI Yongbin, etal. Metal magnetic memory testing technology[M]. Beijing: China Electric Power Press,2000.