金屬管道瞬變電磁檢測有限元分析

XX大學(xué)畢業(yè)設(shè)計（論文）題 目： 金屬管道瞬變電磁檢測有限元分析 學(xué) 院： 測試與光電工程學(xué)院專業(yè)名稱： 測控技術(shù)與儀器班級學(xué)號： 學(xué)生姓名： 指導(dǎo)教師： 二Oxx 年 六月 金屬管道瞬變電磁檢測有限元分析摘要：眾所周知，金屬管道已被廣泛應(yīng)用于石油、天然氣、水等資源運輸，但是隨著管道使用年限的增金屬管道的侵蝕已不可避免。侵蝕將會極大的縮短管道的使用年限。瞬變電磁法利用瞬變電磁手段評估金屬管道的剩余壁厚，是一種建立在電磁感應(yīng)原理基礎(chǔ)上的時間域人工源電磁探測方法。因為瞬變電磁檢測法具有簡單易行、信息量豐富、耦合噪聲小等優(yōu)點，并可實現(xiàn)在役、非開挖檢測，瞬變電磁檢測法已經(jīng)被越來越多的應(yīng)用于埋地金屬管道檢測。 ANSYS是專業(yè)的有限元分析軟件，其功能強大在多領(lǐng)域多變工程問題的求解有著廣泛的應(yīng)用。它可以非常方便的對各種實際中的工程問題進行建模和求解，并可以直觀地體現(xiàn)目前還不易觀測、任何試驗都無法看到的發(fā)生在構(gòu)造內(nèi)部的一些物理現(xiàn)象，且ANSYS的仿真結(jié)果與實際十分接近，于是可以利用ANSYS仿真瞬變電磁法對金屬管道的檢測，為瞬變電磁檢測金屬管道的機理提供依據(jù)。關(guān)鍵詞：瞬變電磁檢測，金屬管道，有限元分析，ANSYS Metal pipes Finite Element Analysis of Transient Electromagnetic DetectionAbstract:As we all know, the metal pipe has been widely used in petroleum, natural gas, water, transportation and other resources, but with the age increase pipeline corrosion of metal pipes have been inevitable. Erosion will greatly shorten the life of the pipeline. TEM method using transient electromagnetic means to assess the remaining wall thickness of metal pipes, is based on the principle of electromagnetic induction time domain electromagnetic method for detecting an artificial source establishment. Because transient electromagnetic detection method is simple, informative, coupling noise, etc., and can be implemented in-service, non-excavation detection, transient electromagnetic detection method has been used more and more buried metal piping testing. ANSYS is a professional finite element analysis software, its powerful and varied in many fields of engineering problems to solve with a wide range of applications. It can be very convenient for a variety of practical engineering modeling and solving problems, and can be directly reflected yet easy observation, we can not see any trial takes place in the internal structure of a number of physical phenomena, and ANSYS simulation results of the actual very close, so you can use ANSYS simulation of transient electromagnetic method to detect metal pipes, provide the basis for the mechanism of transient electromagnetic detection of metal pipes.Keywords: Transient Electromagnetic method, Metal pipes, Finite Element Analysis，ANSYS目 錄1 緒論1.1 瞬變電磁檢測發(fā)展及研究現(xiàn)狀（1）1.2 埋地金屬管道瞬變電磁檢測法的特點和技術(shù)優(yōu)勢（2）2 瞬變電磁理論分析2.1 瞬變電磁法檢測基本原理（4）2.2 時域電磁場的理論基礎(chǔ)（4）2.2.1 微分形式的麥克斯韋方程（4）2.2.2 積分形式的麥克斯韋方程和電磁場的邊界條件（5）2.3 瞬變電磁檢測法信號的分析（7）2.3.1 瞬變電磁信號的動態(tài)范圍（7）2.3.2 瞬變電磁信號的衰減特性和頻帶（7）2.3.3 瞬變電磁信號的取樣（8）2.4 瞬變電磁法常用的激發(fā)場波形（8）3 有限元分析方法及ANSYS簡介3.1 有限元分析方法簡介（10）3.2 有限元分析法的特點（10）3.3 ANSYS簡介（10）4 埋地金屬管道的ANSYS三維仿真分析4.1 創(chuàng)建埋地金屬管道三維模型的物理環(huán)境（12）4.2 三維仿真模型的建立（13）4.3 三維仿真模型的網(wǎng)格劃分（16）4.3 施加載荷及邊界條件（17）4.4 求解及后處理過程（18）4.5 三維管道模型仿真結(jié)果分析（20）4.6 不同壁厚的金屬管道仿真結(jié)果分析（21）5 總結(jié)參考文獻（25）致 謝（26）金屬管道瞬變電磁檢測有限元分析1 緒論1.1 瞬變電磁檢測發(fā)展及研究現(xiàn)狀瞬變電磁檢測技術(shù)的應(yīng)用始于20世紀30年代，當時瞬變電磁檢測法是由前蘇聯(lián)科學(xué)家1提出被用于完善解決地質(zhì)結(jié)構(gòu)問題。從這以后前蘇聯(lián)科學(xué)家經(jīng)過近三十年堅持不懈的努力，TEM（Transient Electromagnetic Method)的一維正演和反演被成功推到出。從這以后TEM發(fā)展十分迅速，隨著瞬變電磁檢測法的詮釋理論和實踐理論的建立，瞬變電磁檢測法已經(jīng)進入了實用階段2。 國外，自從二十世紀70年代以來，通過利用瞬變電磁法，前蘇聯(lián)科學(xué)家開展了大量的實驗工作。除此之外，大洋洲和北美等國的科學(xué)家也在瞬變電磁的理論和實踐兩方面做了研究，這使得瞬變電磁法的應(yīng)用得到了迅速的發(fā)展。此時，TEM已經(jīng)在軍事檢測、考古檢測、油氣勘探、工程勘察和環(huán)境調(diào)查等多方面有著廣泛的應(yīng)用34；與此同時，在TEM儀器設(shè)備方面也取得了很大的成功，一些著名的全球儀器儀表公司也先后推出了不同種類的瞬變電磁檢測儀器。例如美國ZONGE公司生產(chǎn)的GDP-32系統(tǒng)；澳大利亞SIROTEM儀器等。這些儀器的最大的特性就是瞬變電磁檢測系統(tǒng)的接收機能夠?qū)崿F(xiàn)全智能化，這使得TEM系統(tǒng)能夠形成多功能的工作平臺。國內(nèi)，對于瞬變電磁檢測法的探索興起于上個世紀七十年代，中南大學(xué)和地礦部物化探研究所等機構(gòu)率先對瞬變電磁法進行了研究。在一維和二維正演5和反演、瞬變電磁理論的研究、野外實驗和儀器的研制等方面它們做出了很大的貢獻。與此同時，有關(guān)TEM方面的書籍也隨之涌現(xiàn)出來，例如牛之璉編著的時間域電磁法原理6和王長清、祝西里編著的瞬變電磁場的理論和計算等。近年來由于國內(nèi)市場的需要，國內(nèi)TEM儀器的發(fā)展也十分迅速。廊坊物化探研究所研制的IGGETEM系列瞬變電磁儀，中南大學(xué)的SD系列瞬變電磁儀器，吉林大學(xué)研制的ATEM系列瞬變電磁儀等。經(jīng)歷了二十多年的發(fā)展，國內(nèi)的瞬變電磁檢測法已經(jīng)在災(zāi)害、環(huán)境監(jiān)測和工程檢測等方面有了廣泛的應(yīng)用。在工程勘探范疇的應(yīng)用研究，國內(nèi)對此方面的研究變得越來越活躍，利用瞬變電磁法檢測地下的礦物質(zhì)一直以來都是很重要的途徑。另外在檢測石油運輸管道腐蝕78狀況方面，國內(nèi)的王淑英和黃桂柏，通過對輸油管道的實地檢測，成功的驗證了瞬變電磁法能夠檢測埋地輸油管道的腐蝕狀況，但是瞬變電磁檢1測法在對埋地金屬管道小面積的比較嚴重的點狀腐蝕的檢測靈敏度還比較低，檢測的準確性較低。盡管雜散電流和外界磁場等外界條件對檢測的結(jié)果具有較大的影響，對于埋地金屬管道的瞬變電磁檢測這已經(jīng)有了長足的進步。瞬變電磁檢測的方法在對埋地金屬管道腐蝕狀況的檢測過程中可以試驗不停輸和不開挖的在役檢測，對金屬管道以及金屬管道內(nèi)傳輸?shù)奈镔|(zhì)對沒有要求，檢測的結(jié)果都能夠直觀的顯示埋地金屬管道的腐蝕狀況。盡管在廣大科技工作者的共同努力下，瞬變電磁法有了很大的發(fā)展，但我們必須認識到在瞬變電磁的理論和計算，TEM儀器研制等方面與國外還是有著不小的差距，對此我們?nèi)孕枰冻龈嗟暮顾团?，進一步縮小與國外的差距。1.2埋地金屬管道瞬變電磁檢測法的特點和技術(shù)優(yōu)勢周圍土壤的腐蝕和雜散電流的腐蝕是造成埋地管道腐蝕的主要外界條件，長期的腐蝕會造成埋地金屬管道穿孔現(xiàn)象的發(fā)生，這將會極大的縮短埋地金屬管道的使用壽命。因此，埋地金屬管道腐蝕狀況的檢測技術(shù)水平的高低對延長管道的使用壽命以及保證工業(yè)順利生產(chǎn)有著非常重要的影響。在實施瞬變電磁法時，具有空間的可分性和時間的可分性9?？臻g的可分性具體指的是，無論是什么性質(zhì)的脈沖都能夠被分解為不同的頻率的余玄信號之和，不同的延時時刻所檢測到的頻率是不同的，所以在不同的時刻場在埋地金屬管道內(nèi)的傳播速度不同，檢測的深度也會有所不同。時間的可分性具體指的是，由于檢測是在電源信號關(guān)斷以后的時間段進行的，并且也只能解析脈沖間隙時間段內(nèi)的數(shù)據(jù)，因此也就去除了一次磁場的干擾。因為上述兩種可分性，瞬變電磁法檢測埋地金屬管道與其它的方法相比較具有的特點有：（1） 在評估埋地金屬管道的腐蝕程度時，瞬變電磁檢測法是在埋地金屬管道所覆蓋的土壤上方直接檢測埋地金屬管道的缺陷，因此瞬變電磁檢測法能夠?qū)崿F(xiàn)不開挖、在役檢測，這樣便可極大的縮短更換和修復(fù)的時間，提高檢測的效率。（2） 瞬變電磁檢測法檢測埋地金屬管道的噪聲信號主要來自于自然電磁場和人文電磁場，提高發(fā)射功率不僅可以增大信噪比10，而且也能夠提高檢測的靈敏度。從而，可以達到提高檢測深度的目的。（3） 瞬變電磁檢測法對激勵線圈的方位和形狀并沒有嚴格的要求，可依據(jù)實地的檢測條件調(diào)整檢測裝置的結(jié)構(gòu)，因此TEM檢測埋地金屬管道具有靈活、工作簡單，工作效率高等優(yōu)點。 總而言之，瞬變電磁法檢測技術(shù)是一種非接觸式信號加載信號的新興檢測技術(shù)，1將其應(yīng)用于檢測埋地金屬管道，能夠?qū)崿F(xiàn)埋地金屬管道檢測的不開挖、在役檢測，具有簡單易行、檢測精度高等特點，該技術(shù)具有非常廣的應(yīng)用前景。2 瞬變電磁理論分析2.1 瞬變電磁法檢測基本原理 瞬變電磁法基于電磁感應(yīng)原理，它以電磁的差異作為依據(jù)，通過不接入地的回線裝置向地下發(fā)射一次磁場，在一次場的間歇時間段內(nèi)，通過接受裝置測量出二次磁場隨時間的變化，以此來達到探測地下介質(zhì)的分布特征和性質(zhì)。如圖2-1所示。圖2.1 瞬變電磁法檢測原理2.2 時域電磁場的理論基礎(chǔ)2.2.1微分形式的麥克斯韋方程 自上個世紀以來，麥克斯韋（Maxwell)歸納了電磁學(xué)的各項研究成果，創(chuàng)建了論述宏觀電磁場運動規(guī)律的方程組，由此奠定了宏觀電磁場的理論基礎(chǔ)，一個多世紀來，成千上萬的電磁場工程和科學(xué)實驗實踐并沒有發(fā)現(xiàn)與麥克斯韋方程組相違背的例子，這便使得人們堅定不移地相信，通過麥克斯韋電磁場理論，用來解決不同類型的宏觀電磁場問題是可行的。在麥克斯韋方程組中，電磁場特性主要由四個物理參量來表示，因為電磁場本質(zhì)上為矢量場，所以上述的4個物理參量均是矢量，分別用E、D、H和B表示,其中：E稱為電場強度（單位為V/m），D稱為電通量密度或電位移矢量（單位為C/m2),H稱為磁場強度（單位為A/m),B稱為磁通量密度或磁感應(yīng)強度（單位為Wb/m2)。從本質(zhì)上講，帶電粒子運動產(chǎn)生的電流以及其空間的分布是宏觀電磁場產(chǎn)生的根本原因。一個電子帶有的電荷量是電量的最基本的單位，所以電荷量的變化是不連續(xù)的，當然其空間分布也是不連續(xù)的。電磁場也被定義為分立的光子，它也是不連續(xù)的。宏觀尺度是宏觀電磁場成立的必要條件，因此場量和電量的分立特性在宏觀尺度下都是可以忽略的，在空間和時間上他們也都被認為是連續(xù)的。在絕大多數(shù)情況下，表示電磁場的各物理參量一般都是時間t和位矢r的函數(shù)。在連續(xù)媒質(zhì)的空間中，當宏觀尺度條件能夠被滿足時，則可以假設(shè)場是r和t的連續(xù)函數(shù)并且此時具有連續(xù)的導(dǎo)數(shù)。作為電磁場的源或被電磁場所誘導(dǎo)電荷密度（單位為C/m3)和電流密度J（單位是A/m2)也假定為r和t的連續(xù)函數(shù)。在上述條件下，麥克斯韋方程組有如下形式：（2.2.1）（2.2.2）（2.2.3）（2.2.4）其中和J還應(yīng)滿足表征電荷守恒定律的連續(xù)性方程，（2.2.5）對于時變電磁場，以上5個方程中只有3個是獨立的。例如，方程（2.2.3）和（2.2.4）可由方程（2.2.1）、（2.2.2）和（2.2.5）導(dǎo)出。當電磁場的運動規(guī)律能夠被數(shù)學(xué)形式表征出來時，這其實就已經(jīng)對電磁場中所包含的物理參量作了某種性質(zhì)上的假定，這樣便能夠保證完成所必需的數(shù)學(xué)計算。一般來說，在電磁場中被假定的物理參量并不是無窮的，并且它們有足夠的可微性和連續(xù)性，因此它們可以自由地交換積、微分的順序，特殊的例子只有在媒質(zhì)不連續(xù)的情況下才會出現(xiàn)。在遇到特殊的例子時應(yīng)進行特殊處理。2.2.2 積分形式的麥克斯韋方程和電磁場的邊界條件 麥克斯韋方程組的微分形式表征了電磁場的區(qū)域性質(zhì)，僅僅在介質(zhì)連續(xù)的區(qū)域才適用，因為微分運算對被積函數(shù)的要求比積分對函數(shù)的要求更高?？紤]空間中的一個有限區(qū)域，其體積為V，表面為S，A為S的一部分，A的邊界為為C。對于V內(nèi)有定義的矢量函數(shù)F，有如下的高斯（Gauss)定理和斯托克斯（Stokes)定理成立， (2.2.6) (2.2.7) 對方程（2.2.1）和（2.2.2）在A上進行積分，再利用方程（2.2.7），對方程（2.2.3）和（2.2.4）在V上進行積分，在利用方程（2.2.6），就可以得到 (2.2.8)(2.2.9)(2.2.10)(2.2.11) 這就是麥斯威爾方程的積分形式。 麥克斯韋方程組的積分形式是對電磁場的全區(qū)域的描述，它表征的是在某一個區(qū)域中電磁場的總體性質(zhì)。確切的來說，只有在分布論的意義上積分形式和微分形式才可以互相轉(zhuǎn)換。由于他們所表達的電磁場特性的角度不同，因此各有不同的作用。因為麥克斯韋方程組的積分形式適用于介質(zhì)的不連續(xù)處，所以將其應(yīng)用于導(dǎo)出介質(zhì)不連續(xù)處場量應(yīng)滿足的關(guān)系是可行的，這就是電磁場的邊界條件。求解微分形式的麥克斯韋方程組必須滿足磁場的邊界條件。假定有介質(zhì)1和介質(zhì)2構(gòu)成的突變交界面，分別用下角標1和2表征相對應(yīng)介質(zhì)中的場量，n表示交界面上由媒質(zhì)2指向媒質(zhì)1的法向單位矢量，則由方程（2.2.8)（2.2.11）可導(dǎo)出交界面上場量應(yīng)滿足的邊界條件11，(2.2.12) (2.2.13) (2.2.14) (2.2.15) 其中Js為交界面上的面電流密度，s為交界面上的面電荷密度。一般來說，我們會把一些條件理想化，這樣就會使得一些問題得到簡化。例如，當導(dǎo)體的導(dǎo)電性能良好時，我們可以把導(dǎo)體的導(dǎo)電率近似的看作為無限大，并成為理想導(dǎo)體。如果用E、H、D和B表示導(dǎo)體外面的場量，則對理想導(dǎo)體可得如下的邊界條件：(2.2.16)(2.2.17)(2.2.18)(2.2.19)如果求解區(qū)域是開放的，則需要給出電磁場在無窮遠處需滿足的條件?；趫瞿艿挠邢扌砸?，場必須滿足輻射條件。2.3 瞬變電磁檢測法信號的分析 依據(jù)法拉第電磁感應(yīng)定律，當發(fā)射線圈中的激勵電流被瞬間關(guān)斷后，地下的導(dǎo)體介質(zhì)將會產(chǎn)生一個和一次磁場方向相同的感應(yīng)磁場12，這個感應(yīng)磁場也稱為二次磁場。二次磁場的信號強度隨著時間推移逐漸擴散衰減，這個擴散衰減的過程可以分為早、中、晚三個階段。埋地金屬管道某實驗點的瞬變電磁響應(yīng)如圖2-2所示。圖 2-2 埋地金屬管道某實驗點的瞬變電磁響應(yīng)曲線2.3.1 瞬變電磁信號的動態(tài)范圍 瞬變電磁信號的動態(tài)范圍13很大，從一開始的1105uv變到后期的0.1uv。信號的動態(tài)范圍與地下導(dǎo)體介質(zhì)電阻率有關(guān)，不同的電阻率介質(zhì)具有不同的瞬變電磁響應(yīng)。低電阻率的導(dǎo)體介質(zhì)初始瞬變電磁響應(yīng)幅度不大，但衰減速度比較慢；高電阻率的導(dǎo)體介質(zhì)初始瞬變電磁響應(yīng)幅度很大，但衰減速度比較快。信號的動態(tài)范圍不僅與介質(zhì)的電阻率有關(guān)，而且與瞬變電磁法所檢測的試件埋深、大小和形狀等因素有關(guān)。因為瞬變電磁信號的動態(tài)范圍大，所以這就對儀器的分辨率具有較高的要求。特別是當所檢測的信號處于后期的信號時。由于外部的噪聲信號很容易覆蓋較弱的瞬變電磁信號，因此對于后期的較弱的電磁信號需要采用有效的措施進行放大處理或者采用有效的信號處理技術(shù)。2.3.2 瞬變電磁信號的衰減特性和頻帶 瞬變電磁檢測法的信號衰減速度14在不同的時期差異比較大，在早、中期這 兩個階段，信號的衰減速度比較快。然而相對于早、中兩個階段的衰減速度，晚期的信號衰減速度則顯得比較慢。依據(jù)瞬變電磁的這種衰減性質(zhì)，要想精確的測定信號的衰減性質(zhì)，在對試驗點采樣時，必須要符合在相應(yīng)的時間段內(nèi)要取出足夠多的采樣點，采樣點的類型決定了采樣點的時間間隔。在早期和中期兩個階段，由于信號衰減的速度比較快，因此取樣點的時間間隔應(yīng)取得足夠小才能夠辨別出信號的衰減特性。對于晚期的信號特點，因為晚期這一階段的信號衰減速度比較慢，所以在這一階段取樣點的時間間隔應(yīng)取得寬一些。瞬變電磁檢測法的信號頻帶較寬，其頻率范圍最高可以達到104赫茲。由于信號的頻帶范圍很寬，這其中也包括各種雜質(zhì)電磁噪聲信號，因此必須對檢測的信號數(shù)據(jù)進行濾波處理，并采用有效的方法用于增強儀器的抗干擾能力，除此之外，還應(yīng)采用有效的手段對晚期較弱的電磁信號進行放大用于提高信噪比。2.3.3 瞬變電磁信號的取樣 模擬積分取樣、模擬積分-數(shù)字化疊加取樣和數(shù)字化疊加取樣是目前各種類型的TEM檢測系統(tǒng)所采用的三種取樣方式15。模擬積分取樣指的是分時間段通過多個積分器在電路上采集電信號，這種方法的優(yōu)點是采樣的準確度比較高，但缺點是對電路的構(gòu)造要求比較嚴格，電路設(shè)計的難度比較大。采用該方法所取出的信號的特點是采樣的次數(shù)、采樣的窗口寬度和采樣的時間由時標信號確定，通常取樣的次數(shù)和取樣的時間是固定的；數(shù)字化疊加取樣指的是經(jīng)過取樣保持電路后，首先對每個取樣點的信號進行數(shù)字化處理，經(jīng)處理的信號直接送入計算機內(nèi)，在計算機內(nèi)進行內(nèi)處理后，最后再進行軟件積分。這種方法的特點是電路簡單、使用比較靈活，在計算機的內(nèi)部可以通過各種方法進行處理，用于提高儀器的性能和檢測信號的能力；模擬集成數(shù)字覆蓋采樣是上述兩種方法相結(jié)合的抽樣方法，結(jié)合這兩種抽樣方法的優(yōu)點，其特點是早期延遲數(shù)據(jù)分辨率較低，數(shù)據(jù)精度晚延遲較高。2.4 瞬變電磁法常用的激發(fā)場波形 三角形波、方波、半正弦波、梯形波和偽隨機波等激勵信號是瞬變電磁法激勵信號波形常用的多種周期性脈沖信號。依據(jù)傅里葉頻譜解析理論我們知道，無論是哪一種脈沖信號都可以被分解為一系列具有不同頻率的諧波信號。因此，無論是哪一種脈沖激勵信號在導(dǎo)體介質(zhì)中產(chǎn)生的信號都可以被看作是一系列具有不同頻率的正弦信號所產(chǎn)生的響應(yīng)之和。圖2-3為目前常用的激勵信號波形。 圖 2.3 常用的激勵信號理論波形和實際波形 (a)雙極性矩形 （b)雙極性梯形 (c)雙極性半正弦 (d)實際發(fā)射波形 通過ANSYS有限元分析14軟件對埋地金屬管道建立三維模型15時，脈沖波形采用單個階躍波。令激勵信號序列的周期趨于無窮大，于是雙極性矩形激勵信號就能被轉(zhuǎn)化為非周期的單階躍波形。然而在實際的工作過程當中，因為激勵線圈的電感、電容特性，輸出的波形不可能是理想的，通常波形應(yīng)呈指數(shù)上升，下降也應(yīng)呈斜階躍下降，如圖2-3(d)所示。3有限元分析方法及ANSYS簡介3.1 有限元分析方法簡介 在確定分析對象的基本特性，并建立起其模型以后，有限元分析方法作為一種分析及計算方法，它可以被概括為如下三點： （1）將所要分析的對象分解為若干個單元，并通過單元邊界上的節(jié)點把各個單元聯(lián)接為一個整體。 （2）將全求解域內(nèi)待求的未知場變量用各個單元內(nèi)假定的近似函數(shù)來分片的表征。而用未知場函數(shù)在單元各個節(jié)點上的數(shù)值及與其對應(yīng)的插值函數(shù)來表征每個單元內(nèi)的近似函數(shù)。因為場函數(shù)在鏈接單元的接地啊上具有相同的數(shù)值，所以把它們作為基本未知量用于數(shù)值求解。這樣便使得原待求函數(shù)的無窮多自由度的求解問題就被轉(zhuǎn)化為場函數(shù)節(jié)點值的有限自由度求解問題。 （3）利用和原問題數(shù)學(xué)模型等效的加權(quán)余量法或變分原理，建立常微分方程組或代數(shù)方程組用于求解基本未知量。此方程組由規(guī)范化的矩陣形式所表征，此后可以用相應(yīng)的數(shù)值分析方法求解該方程并得到原問題的解答。3.2 有限元分析法的特點 （1）復(fù)雜幾何構(gòu)形的適應(yīng)性。因為單元在空間上可以被定義為一維、二維或者三維的，且每一種單元可以采用不同的連續(xù)方式也可以具有不同的形狀。 （2）各種物理問題的適應(yīng)性。因為全求解域的未知場函數(shù)被單元內(nèi)的近似函數(shù)所表征，其并未設(shè)定限制場內(nèi)方程的形式應(yīng)滿足相對應(yīng)的條件，也未設(shè)置各個單元所對應(yīng)的方程必須有相同的形式，所以它適用于各種物理問題，例如動力問題、屈曲問題、粘彈性問題等，而且其還可以被應(yīng)用于不同類型復(fù)雜物理現(xiàn)象相互耦合的問題。 （3）嚴格理論基礎(chǔ)上的可靠性。由于邊界條件和微分方程的等效微分形式是加權(quán)余量法或有限元方程的變分原理，因此當原問題的數(shù)學(xué)模型是無誤的，同時用于解析有限元方程的數(shù)值運算方法是可靠穩(wěn)定的，則隨著單元尺寸的縮小或是隨著單元自由度數(shù)的增加，有限元分析方法求解的近似程度將不斷地被提高。 （4）計算機實現(xiàn)的高效性。因為有限元分析方法的每一步都能夠被表示為規(guī)范化的矩陣形式，這樣便使得求解方程可以被看作是標準的矩陣代數(shù)問題，這種情況下特別符合計算機的執(zhí)行和編程。計算機硬件技術(shù)的不斷發(fā)展以及數(shù)值運算方法的更新?lián)Q代，大規(guī)模的復(fù)雜工程的有限元分析已經(jīng)被看作是工程技術(shù)領(lǐng)域的常規(guī)工作。3.3 ANSYS簡介ANSYS是以有限元分析法為依據(jù)的CAE(Computer-aided Engineering)軟件，它是由美國ANSYS公司所研發(fā)。ANSYS無論是在PC機、工作臺亦或是巨型計算機都具有較好的兼容性，ANSYS文件在其所有的系列產(chǎn)品和工作平臺上都可以運行。這就使得 ANSYS用戶可求解多領(lǐng)域多變工程問題。它可以把電場、磁場、聲場的分析融為一體，ANSYS現(xiàn)已被廣泛應(yīng)用于航天、土木工程、水利水電工程等領(lǐng)域。 ANSYS 主要由以下三個部分組成： （1）前處理模塊(PREP7)前處理模塊完成有限元分析法的前處理過程。其實質(zhì)是為麥克斯韋方程組及相關(guān)方程組設(shè)置參數(shù)和初始條件。模塊提供了強大的實體建模及網(wǎng)格劃分工具，主要包括三部分內(nèi)容：材料屬性的定義及配置、實體建模和單元網(wǎng)格劃分。 （2）加載和求解模塊(SOLUTION)加載和求解模塊完成有限元分析法的計算過程。用于對有限元模型選擇分析類型、施加外界條件和載荷、確定載荷步、選擇求解器并進行求解。 （3）后處理模塊(POST1 和 POST26)。后處理模塊完成有限元分析法的后處理過程?？梢杂脠D形方式顯示模型結(jié)構(gòu)內(nèi)部的計算結(jié)果，比如：等值線顯示、梯度顯示、矢量顯示、透明及半透明顯示等，也可以圖片、曲線形式顯示或輸出計算結(jié)果。4 埋地金屬管道的ANSYS三維仿真分析4.1 創(chuàng)建埋地金屬管道三維模型的物理環(huán)境（1） 過濾圖形界面：從主菜單中選擇Main MenuPreferences,彈出“Preferences for GUI Filtering”對話框，選中“Magnetic-Nodal”來對后面的分析進行菜單及相應(yīng)的圖形界面過濾。如圖4.1所示： 圖 4.1 選擇圖形界面（2）定義工作標題：執(zhí)行菜單欄中的 Utility MenuFileChange Title,在彈出的對話框中輸入“3D LYBS”，單擊“OK”.（3）指定工作名：執(zhí)行菜單欄中的 Utility MenuFileChange Jobname,彈出一個對話框，在“Enter new name”后面輸入“LYBS_3D”，單擊“OK”按鈕。如圖4.2所示： 圖 4.2 定義工作名 （4）選擇單元類型：從主菜單中選擇 Main MenuPreprocessorElement TypeAdd/Edit/Delete,彈出“Element Types”單元類型對話框，單擊“Add”按鈕，彈出“Library of Element Types”單元類型對話框。在該對話框中左面滾動欄中選擇“Magnetic-vector”,在右邊的滾動欄中選擇“Brick 8 node 97”,單擊“OK”按鈕，定義了一個“SOLID97”單元。如圖4.3所示： 圖 4.3 選擇單元類型 （5）定義材料屬性：從主菜單中選擇Main MenuPreprocessorMaterial PropsMaterial Models,彈出“Define Material Model Behavior”對話框。在右邊的欄中連續(xù)單擊“ElectromagneticsRelative PermeabilityConstant”后，又彈出“Permeability for Material Number 1”對話框。在該對話框中“MURX”后面的輸入欄輸入“1”，單擊“OK”按鈕。最后單擊“MaterialExit”結(jié)束。這樣便定義了空氣的相對磁導(dǎo)率，如圖4.4所示： 圖 4.4 定義空氣的相對磁導(dǎo)率 重復(fù)上述的步驟，分別定義管道、線圈、土壤的相對磁導(dǎo)率為350、1、5。同理可在主菜單中的“Material Models”中分別定義線圈的電阻率為3E-8歐姆/米，管道的電阻率為1E-7歐姆/米，土壤的電阻率為50歐姆/米。4.2 三維仿真模型的建立 自頂向下和自底向上是ANSYS中兩種實體建模的方法。 （1）自頂向下自頂向下的建立模型指的是按照從體到面、從面到線、從線到點的順序進行建模，因為線是由點構(gòu)成，面是由線構(gòu)成，體是由面構(gòu)成，所以稱這個順序為自頂向下建模。 （2）自底向上自底向上建模與自頂向上建模正好相反，是按照從點到線，從線到面，從面到體的順序建立模型，因為線是由點構(gòu)成，面是由線構(gòu)成，而體是由面構(gòu)成，所以稱這個順序為自底向上建模。本次三維仿真模型的建立采用自頂向下的建模方法。依據(jù)自頂向下的建模方法，第一步先定義模型的最高級的圖元，在三維仿真模型中需要定義線圈、賣地金屬管道、空氣以及土壤。當上述模型中的內(nèi)容被定義以后，ANSYS會自動匹配相關(guān)的關(guān)鍵點、線和面。用戶圖形界面（Graphic User Interface 簡稱GUI）操作方式和參數(shù)化語言設(shè)計（ANSYS Parametric Design Language 簡稱APDL）方式是ANSYS的兩種操作方式。GUI方式不要求掌握命令的使用格式和編程語言的使用規(guī)則等，用戶可以只通過鼠標便可在圖形界面上進行操作。GUI操作方式對于比較簡單的有限元分析模型，操作速度或許會更快些。而對于比較復(fù)雜的、規(guī)模比較大的有限元分析模型，GUI操作方式的缺點就會有所顯現(xiàn)。由于對于一個模型的分析經(jīng)常需要多次的反復(fù)分析，當需要對修改后的模型進行分析時，若采用GUI操作方式就會出現(xiàn)大量的重復(fù)步驟。大量的計算時間被這些重復(fù)的工作占據(jù)，因此在這種情況下GUI操作方式就顯得比較繁瑣，效率也比較低。APDL參數(shù)化語言設(shè)計操作方式，可以自動地完成部分常規(guī)分析操作。它是一種依據(jù)參數(shù)化編程方式和操作命令來建立分析模型的腳本語言。智能分析，自動地完成多數(shù)GUI操作等都是APDL編寫的腳本程序進行參數(shù)化建模的特點。特別的是APDL操作方式甚至能夠符合GUI操作方式無法滿足的要求。APDL操作方式可以分析屬性和控制相關(guān)的設(shè)計，并允許輸入復(fù)雜的數(shù)據(jù)。比如定義材料的屬性、輸入模型的尺寸、定義邊界條件的位置、施加載荷和調(diào)整網(wǎng)格劃分的密度等。通過APDL操作方式不僅能夠為問題的計算與求解提供更加簡單高效的手段，而且極大的擴展的有限元分析的范圍。由于本次三維仿真模型的建立相對來說比較復(fù)雜，因此采用APDL操作方式。首先建立線圈的模型，定義線圈的匝數(shù)N=10，線圈的外徑為0.102m，線圈的內(nèi)徑為0.1m，線圈的厚度為0.0035m，輸入APDL命令流后得到的模型如圖4.5所示： 圖 4.5 定義線圈的模型 其次再建立管道的模型，定義管道的長度為1.5m，外徑為0.65m，內(nèi)徑為0.52m，管道的埋深為0.5m，輸入APDL命令流后得到的模型如圖4.6所示： 圖 4.6 定義管道的模型 最后，再建立土壤和空氣的模型，定義空氣層的高度為0.2m，定義土壤的模型為長1.5m，寬0.5m，高為0.75m的長方體，土壤將包裹在管道的周圍。輸入APDL命令流后，添加土壤后及最終建立的模型分別如圖4.7、4.8所示： 圖 4.7 添加入土壤后的模型 圖 4.8 埋地金屬管道模型4.3 三維仿真模型的網(wǎng)格劃分在ANSYS有限元分析法中，計算時間和計算精度等因素決定所選取網(wǎng)格的大小。當所要求的單元越小，則網(wǎng)格的劃分就越細，同時計算結(jié)果的精度也越高。盡管網(wǎng)格劃分的越細計算得到的結(jié)果精度也越高，但由于網(wǎng)格單元數(shù)的急劇增加，由此所增加的計算量將會十分龐大，同時計算的時間也會延長。對于此次埋地金屬管道的三維仿真模型采用的是自由劃分和映射劃分相互結(jié)合的方式。采用這種劃分方式，一方面仿真計算時間減少了；另一方面，仿真的計算精度也得到了提高。在三維模型中，線圈和管道采用的是六面體網(wǎng)格，通過掃略方式生產(chǎn)體網(wǎng)格，40個節(jié)點被等分在線圈周長上，同時線圈的厚度方向也等分3個節(jié)點。管道的軸向選中四條直線，管道長為1.5m,在管道軸向方向每0.015m設(shè)置一個節(jié)點。為了減少計算時間，同時由于實際空氣和土壤的結(jié)構(gòu)比較復(fù)雜，空氣和土壤均采用四面體網(wǎng)格進行自由網(wǎng)格劃分，自由網(wǎng)格劃分的劃分等級為1。如圖4.9、4.10、4.11分別是線圈的網(wǎng)格劃分、管道的網(wǎng)格劃分及整個埋地金屬管道模型的網(wǎng)格劃分方式： 圖 4.9 線圈網(wǎng)格劃分圖 圖 4.10 管道的網(wǎng)格劃分圖 4.11 埋地金屬管道模型的網(wǎng)格劃分4.3 施加載荷及邊界條件諾依曼邊界條件和狄利克萊邊界條件是ANSYS分析軟件中的兩類邊界條件。一般來說狄利克萊是作為約束條件被提出，它指的是磁力線與模型邊界相互平行的條件，也叫做第一類邊界條件。諾依曼邊界條件指的是磁力線垂直于模型邊界條件，也叫做第二類邊界條件。對于此次埋地金屬管道三維仿真模型中，僅僅在模型的外表面添加磁力線平行邊界條件，即狄利克萊邊界條件。（各條邊界線上設(shè)置AZ=0）同時通過輸入APDL命令流設(shè)置加載在線圈上的電壓為12伏特。則施加載荷及邊界條件以后的三維仿真模型如圖4.12所示： 圖 4.12 施加載荷及邊界條件后模型圖4.4 求解及后處理過程當線圈被加載上電壓為12V的激勵信號后，輸入SOLVE命令后，線圈材料的每個單元節(jié)點上將被ANSYS自動加載12V的激勵電壓并開始計算。當所有求解完成后，再通過POST1通用后處理器查看激勵電壓關(guān)斷前、關(guān)斷后（t0.03001s為關(guān)斷后）壁厚為13mm且埋深為0.5m的管道的磁通密度分布云圖。不同時刻管道的磁通密度分布云圖如下圖所示。（圖中不同的顏色代表不同磁通密度的大小，其中紅色的表示數(shù)值最大，藍色的表示數(shù)值最小。） t=0.015s t=0.03001s t=0.03825s t=0.06325s t=0.10001s 圖4.13同一管道不同時刻磁通密度分布云圖4.5 三維管道模型仿真結(jié)果分析 由圖4.13金屬管道管壁的磁通密度分布云圖可以得知： (1) 在激勵電壓信號未關(guān)斷之前，線圈正下方的管道管壁磁場強度較強，而管道兩端管壁的磁場強度較弱，且在激勵電壓信號未關(guān)斷之前這種狀態(tài)都比較穩(wěn)定。 (2) 激勵電壓信號關(guān)斷初期，仍然是線圈正下方的管道管壁磁場強度較強，而管道兩端管壁的磁場強度較弱。從激勵電壓信號關(guān)斷中、后期的磁通密度分布云圖可知，管道中間部分的磁場強度逐漸減小，磁場強度有向管道兩端擴散的趨勢，因此管道兩端的磁場強度逐漸增大。最后，管道中間部分的磁場強度首先減少到最小，磁場大部分集中在管道的兩端。4.6 不同壁厚的金屬管道仿真結(jié)果分析管道腐蝕程度的不同主要表現(xiàn)在管道壁厚的不同，分別建立壁厚為10mm、13mm、16mm和18mm的三維仿真模型，分別計算在關(guān)斷激勵電壓信號后的同一時刻（t=0.06325s)的磁通密度分布云圖，如圖4.14所示。由圖中可以得知，在關(guān)斷電壓激勵信號后，隨著管道壁厚的增加，磁場強度由管道中間部分向管道兩端擴散的速度也越來越快。 壁厚d=10mm壁厚d=13mm 壁厚d=16mm 壁厚d=18m圖 4.14不同壁厚磁通密度分布云圖 在完成上述個模型的ANSYS三維仿真后，通過POST26時間歷程處理器計算出各個模型的線圈上的感應(yīng)電流并導(dǎo)出所有計算出的電流值，再利用MATLAB將導(dǎo)出的線圈上的感應(yīng)電流值轉(zhuǎn)化為感應(yīng)電壓值，并畫出圖像，如圖所示： 圖4.15不同腐蝕程度的MATLAB圖形 圖4.16不同腐蝕程度放大后的圖形 從上述的圖形中可以看出，當管道的壁厚不相同但金屬管道的外徑相同時，線圈的瞬變響應(yīng)曲線是可分的，管壁16mm的感應(yīng)電壓值衰減最慢，管壁為10mm的感應(yīng)電壓值衰減最快，而管壁為13mm的感應(yīng)電壓值的衰減速度處于上述兩者之間。即腐蝕程度越大的管道，感應(yīng)電壓值的衰減速度越快。5 總結(jié) 本文研究了金屬管道瞬變電磁檢測的有限元分析，所做的工作總結(jié)如下：（1） 理論準備。針對所研究的課題，了解熟悉了瞬變電磁檢測的發(fā)展及研究現(xiàn)狀和埋地金屬管道瞬變電磁檢測法的特點和技術(shù)優(yōu)勢。通過學(xué)習(xí)研究瞬變電磁法檢測埋地金屬管道技術(shù)以及ANSYS仿真原理，為金屬管道瞬變電磁檢測的有限元分析做好準備。（2） 仿真實驗。首先通過GUI操作方式與APDL命令流相互結(jié)合的方式，建立了包含有空氣、線圈、土壤和管道的埋地金屬管道的三維仿真模型；其次，將12V的電壓激勵信號加載在線圈上，利用ANSYS通用后處理器計算出激勵電壓信號關(guān)斷前后不同時刻金屬管道的磁通密度分布云圖。最后，比較了激勵電壓信號關(guān)斷后壁厚分別為10mm、13mm、16mm和18mm在同一時刻的磁通密度分布云圖以及不同腐蝕程度的管道的瞬變感應(yīng)電壓值變化曲線。（3） 依據(jù)仿真實驗可以得出：一方面，在激勵電壓信號關(guān)斷的初期，線圈正下方的管道管壁的磁場強度較強，并有向管道兩端擴散的趨勢。管道中間部分的磁場強度下降較快，最后，管道中間部分的磁場強度最弱，磁場大部分集中在管道的兩端；另一方面，當埋深、所加激勵信號等條件一致而只有壁厚不同時，隨著壁厚的增加，磁場向管道兩端擴散的速度也越來越快；金屬管道埋深等條件相同、且只有壁厚不相同同時，感應(yīng)電壓曲線的尾支顯然是可分的，且當腐蝕程度越大（管壁越薄）感應(yīng)電壓的衰減速度越快。參考文獻1 冷元寶等.我國堤壩隱患及漏磁探測技術(shù)現(xiàn)狀及展望J.水利水電科技進展，2002，22（2）：