1、一、概述國內外研究狀況燃氣渦輪發動機的主要發展方向是提高發動機渦輪前燃燒溫度、增加推重比和提升渦輪發動機部件在包括腐蝕和氧化等嚴酷服役環境下的熱效率。隨著發動機燃燒溫度、推重比和熱效率的提高，發動機熱端部件，特別是燃燒室中的燃氣溫度和燃氣壓力不斷提高，而現有的高溫合金和冷卻技術難以滿足需要，熱障涂層（Thermal Barrier Coatings,TBCs）技術1, 2得到了廣泛的重視。熱障涂層是金屬緩沖層或粘結層和耐熱性、隔熱性好的陶瓷熱保護功能涂層組成的“層和型”金屬陶瓷復合涂層系統。表面陶瓷工作層是借助于基體和陶瓷層之間的金屬粘結層與高溫基體結合。此中間過渡層具有優異的耐高溫、抗氧化性

2、能，熱膨脹系數介于基體金屬與表面陶瓷層之間，減少了陶瓷涂層與合金基體之間熱實配問題，可減緩界面應力，提高涂層結合強度、抗熱震性能和工作壽命。隨著燃氣溫度的不斷提高，如何獲得性能更優異、壽命更持久的熱障涂層已經成為研究人員迫切需要解決的重大難題。熱障涂層主要包括雙層系統、多層系統和梯度系統1。這三種結構形式各有特點，針對不同的環境要求，可以采用不同的結構體系。多數實際應用的熱障涂層采用雙層結構如圖1(a)所示，這種結構制備工藝相對簡單、耐熱性強，但由于涂層熱膨脹系數在界面躍變較大，在熱載荷作用下涂層內容易積聚較大的應力，因此抗熱震性能難以得到進一步提高。為了緩解涂層內的熱效應匹配問題，提高涂層整

3、體抗氧化及熱腐蝕能力，發展了多層結構系統（圖1(b)）。其每一層都具有各自的特定功能，外層封閉層和阻擋層主要用于阻擋燃氣腐蝕產物的侵蝕，擴散阻擋層則用于降低氧原子進一步向涂層內擴散。多層體系結構的力學行為更為復雜，制備也相對困難，付諸實際運用很少。日本學者新野正之、平井敏和渡邊龍三等于1987年首次提出了兩種或多種金屬與陶瓷材料復合3，制備一種在結構和組分上呈連續梯度變化的新型梯度化結構材料，即功能梯度材料（Functionally Graded Materials，FGM）的新概念。這種梯度涂層消除了層狀結構的明顯層間界面，使得涂層的力學性能由基體向陶瓷表層連續過渡，從而避免了熱膨脹系數等不

4、匹配造成的陶瓷層過早剝落現象。由于高溫環境下，梯度涂層內彌散分布MCrAlY金屬組元氧化對系統壽命的影響機理尚不清楚，且功能梯度熱障涂層的制備工藝更為復雜，工藝優化參數未統一，在渦輪發動機葉片等高溫部件上得到實際運用還有許多問題需要解決。(a) (b) (c)圖1熱障涂層的基本設計思想就是利用陶瓷的高耐熱性、抗腐蝕性和低導熱性，實現對基體合金材料的保護，普通金屬材料已經無法滿足要求。NASA對幾種可能適用于高溫隔熱涂層使用的陶瓷材料進行的對比研究表明ZrO2的綜合性能最好。ZrO2具有陶瓷材料中最接近金屬材料的熱膨脹系數，導熱系數為Al2O3的1/100左右，具有較高的抗彎強度和斷裂韌性，尤其

5、是部分穩定的ZrO2特有的微裂紋和相變增韌機制，使得抗熱震性能非常好。但ZrO2在通常使用溫度范圍內，會發生馬氏體相變，伴隨有4%左右的體積變化4，會增加涂層內的應力。在熱循環狀態下這種相變將導致ZrO2涂層的碎裂。為了使熱障涂層能夠適應這一相變溫度左右的熱循環工作環境，延長涂層壽命，在ZrO2中添加少量的氧化物穩定劑可以起到控制這種相變發生的作用。目前應用廣泛的穩定劑是Y2O35。尋求更高性能的陶瓷表層材料和更好的氧化物穩定劑，一直是熱障涂層研究中的重要和熱點方向。為了緩解陶瓷涂層和基體的熱不匹配問題，同時也為了提高基體抗氧化性，在基體和陶瓷涂層間加了一層MCrAlY金屬粘結層，這層涂層成分

6、可以依據使用條件的不同而調整，不受基體成分的限值，而且厚度也可以調控。MCrAlY涂層的抗氧化機理一般是通過高溫氧化環境中，在表面首先形成Al2O3保護性氧化層以阻止涂層的進一步氧化，達到保護基體的目的。MCrAlY粘結層的成分對粘結層在熱循環過程中熱氧化層的生長速度、成分、完整性、與基體的結合力和剝落行為有決定作用，因此其選擇對熱障涂層的使用壽命非常重要。目前飛機發動機葉片使用較多的是抗氧化和抗熱腐蝕綜合性能較好的NiCoCrAlY粘結層6。先進TBCs體系的首要設計目標是提高耐久性，尤其是防止陶瓷頂層碎裂。改進涂層性能、預測涂層壽命，首先必須要了解熱障涂層在熱循環過程中的失效行為和失效機制

7、。導致涂層碎裂的因素很多，包括外來粒子入侵、陶瓷頂層開裂、粘結層與陶瓷層界面開裂等。前兩類因素通常導致TBCs連續或部分退化，而界面開裂導致TBC的主要部分甚至整體碎裂，更應引起重點關注。雖然陶瓷層中的相變、燒結等是促進開裂的重要因素，但與粘結層相關的界面開裂問題對TBCs的壽命影響更大。不管用何種工藝制備的TBCs系統，在高溫應用環境下，在陶瓷涂層和粘結層界面上都將形成另外一層新的氧化層（Thermally Grown Oxide，TGO），其主要成分是Al2O3。隨著循環次數的增加，氧化層厚度會逐漸增厚，而TGO和陶瓷層的界面會限制這種體積變化，其內容易形成微裂縫和應力高度集中，從而其附近

8、很容易導致陶瓷涂層剝落。因此現今使用的多種TBCs體系雖然沒有一個統一確定的失效部位，但熱障涂層的失效通常發生在TBC界面及其附近，TGO是整個體系中最薄弱的環節。熱障涂層的高溫失效行為非常復雜，影響涂層使用壽命的因素也非常多樣，包括涂層中孔隙與微裂紋的數量、類型及其分布情況，熱循環過程中應力的大小和分布情況，粘結層的抗氧化性能，陶瓷層相結構的熱穩定性能，各層界面的化學相容性，熱循環導致粘結層發生的蠕變行為，金屬基體的粗糙度以及涂層材料的本征性能等等。但是，以上種種因素對涂層使用壽命的影響歸根結底都是通過影響熱應力大小、方向和分布而發揮作用的。在眾多材料和加工等影響變量中，TGO附近的局部應力

9、和界面強度是導致失效的兩個最重要的因素7-9。應力場分析是熱障涂層破壞機理研究的基礎，主要包括制備過程中的殘余應力及熱循環過程的應力生長，如TGO形成過程中和之后系統中產生的熱生長應力，了解其特性對TBCs系統的各種壽命預測模型的建立具有重要意義。能精確地定量分析結合材料界面區域的應力場分布，是準確分析和評價涂層材料界面的力學行為以達到對此類涂層材料或結構形式的結合強度和壽命評價及其優化設計目的的關鍵。近幾十年來諸多研究者從不同途徑如理論分析、數值計算、實驗測試等入手開展了很多研究。對于結合材料如薄膜涂層結構形式的研究，由于界面的存在，理論分析比均質材料復雜很多。比較早的研究可以追溯到1955

10、年Rongved10分析了二維和三維的兩個半無限體結合材料內部受集中力作用的問題，通過引入一個和力作用點關于界面對稱的鏡像點，在滿足界面連續性條件的基礎上求得了理論解，可以退化到半無限體的Mindlin解11。Dundurs和Hetenyi12求解了兩個半無限體以無摩擦接觸方式結合內部受點力作用的理論解。由于必須滿足界面上面力和位移連續性的要求，因此尋求該問題理論解時常常采用積分變換法。Bogy13在分析不同材料參數和幾何參數條件下的界面端奇異應力場特性時采用了梅林變換；Burmister14-16利用Fourier變換求解了雙層彈性體問題，丁皓江17進一步擴展至多層彈性體問題。另外還有諸多學

11、者例如Fares和Li18、Yu和Sanday19、Walpole20等求出了不同情況下集中力作用于各向同性結合異材問題的格林函數。盡管格林函數可以通過數學方法如快速Fourier變換等計算出，但由于涉及到無窮積分，解本身已很復雜，并且準確性也有待驗證，而且積分變換常遇到求解逆變換的問題，通常只在一些特殊情況下才能求得顯示解。對于單層或多層薄膜涂層這種層狀系統，許多學者通過多種途徑尋求不同問題的格林函數以進一步分析更復雜的問題，或者作為在無限域或半無限域具有特殊優勢的邊界元法的基本解進行數值分析。例如狀態空間法21, 22，但由于涉及的數學算子通常需要特定的邊界條件才能運算，受限較多；借助傅立

12、葉或Hankel變換的矩陣轉換法應用較多，由于分析中僅保持了一次項，避免了傅立葉或Hankel的逆變換，但精確度受層厚影響較大；另外還有轉制矩陣法、剛度矩陣法、柔度矩陣法、薄層理論法等等23-26，不再一一描述。還有一些學者采用和電磁場理論中相似的一種途徑鏡像點方法來研究異材結合材料問題，例如上面提到的Rongved10針對集中力作用于兩半無限體結合材料內部時的研究。從上面可以看出，目前針對含有界面的涂層系統理論分析還存在一些問題，例如解的完備性不足，或者針對性太強應用價值低等等。涂層測試方法可分為定性檢測和定量檢測兩大類：定性方法有柵格試驗、杯突試驗、熱震試驗、銼磨法和超聲法等；定量測試有粘

13、結拉伸法、壓痕法、剪切試驗和劃痕試驗等。定性法以經驗判斷和相對比較為主，結果一般很難給出力學參量。定量測量的主要困難在于尋求試驗的試件形狀和加載方式，使得界面上能夠產生不同的應力狀態，即在不同的應力狀態比下發生破壞。目前精確測定應力有效的手段是光激發熒光譜技術（Photo-stimulated Luminescence Spectroscopy, PSLS）。由于數值計算更容易修改涂層厚度和材料物性參數，國內外在數值模擬分析涂層結構應力場方面做了許多研究，主要采用有限元法、邊界元法和界面元法等數值方法對涂層構件應力場進行了計算分析。比較早的如1974年Kennedy和Ling27模擬分析了球形

14、壓頭作用于層狀體表面的壓痕測試實驗，隨后一些相似的分析有Kral和Komvopoulos28-30的硬質表層效應研究，Djabellaa和Arnella31對平面和空間赫茲分布壓力作用下的膜基系統進行了有限元分析。很多研究者還借助有限元途徑從不同角度分析了影響薄膜涂層材料界面結合強度的一些因素，如Peng和Bhushan32通過建立三維數值模型進行了粗糙表面的彈塑性層狀半無限體在不同環境下的接觸分析，研究了層基材料剛度、層厚、外力對接觸分析結果的影響，并提出了減小摩擦、粘滯和磨損的優化參數；Tian和Saka33通過有限元途徑分析了不同摩擦系數時雙層薄膜受滑動接觸的應力和應變分布情形，認為表面

15、變形、初始屈服點以及界面應力應變都和摩擦系數緊密相關。另外一些數值計算途徑如邊界元法，對某些無限域以及半無限域問題，由于邊界元自身特點在計算效率和準確性相對有限元法提高很多，同樣在薄膜涂層材料方面有許多的研究。美國國家航空航天局（NASA）在超級計算機上利用 Marc 軟件仿真熱障涂層的溫度場及應力場，進而探究熱障涂層失效機理。但在實際中，數值計算尚存在計算精度和效率較低的問題，而且計算費時費力效率低。另外由于界面上可能存在的界面端奇異性及應力不連續性等，往往還需對界面上的應力分布即應力場有一個全面的把握，而這項工作卻是難以僅靠數值計算來完成的，必要的理論分析結合是不可避免的。涂覆熱障涂層的部

16、件能否安全可靠服役，很大程度上取決于TBC的壽命。建立針對某一具體服役條件下熱障涂層的壽命預測模型，對渦輪發動機的設計和發展有著重要的指導意義。但是影響涂層壽命的因素非常復雜，除了基體材料和粘結層的成分、制備工藝、微觀組織結構外，還有使用環境溫度和腐蝕程度、碰撞損傷等等。建立壽命評估及預測模型需要了解熱障涂層系統的力學行為，并能估測涂層失效的時間。壽命預測系統的建立一般包括判斷關鍵的失效機理、建立應力或應變彈性場模型以及根據應力狀態和相關的失效標準建立最終的壽命評估準則。目前提出的各種壽命預測模型主要是根據實驗結果的經驗或半經驗公式34。熱障涂層服役過程中，不斷的經歷加熱冷卻加熱的工作環境。每

17、次熱循環中，界面氧化、高溫燒結及熱不匹配等因素引起的殘余應力不斷對涂層系統造成損傷。當損傷變量累積達到某一程度后，定義為熱障涂層失效。研究人員從損傷變量的變化和循環累積的角度進行壽命評估，如Busso等34利用有限元建立熱障涂層“波浪”模型，假定大氣等離子噴涂TBCs系統內的殘余應力主要由高溫氧化、燒結、材料參數不匹配引起，疲勞損傷主要收系統內的離面應力的影響，而提出的基于連續損傷力學的壽命預測模型。Liu等35結合熱震試驗及有限元分析，提出隨著循環次數的增加，熱循環應力不斷增加，所帶來的損傷不斷累積，最終也建立了一種熱障涂層復雜載荷可用單一參數來表達的壽命預測模型。另外還有基于特定參量監測和

18、界面氧化等不同角度提出的壽命預測模型36, 37。但熱障涂層壽命預測模型的建立是一個長期而又復雜的過程，要完全定量分析熱障涂層的失效過程并對其壽命進行預測和評估需要諸多重視并長期堅持。參考文獻：1.Padture, N. P.; Gell, M.; Jordan, E. H., Thermal Barrier Coatings for Gas-Turbine Engine Applications. Science 2002, 296 (5566), 280-284.2.周益春; 劉奇星; 楊麗, et al., 熱障涂層的破壞機理與壽命預測. 固體力學學報 2010, (5), 504-53

19、1.3.Viswanathan, R.; Scheirer, S. T. In Materials technology for advanced land based gas turbines, The Proceeding of 7th International Conference on Creep and Fatigue at Elevated Temperatures, Tsukuba, Japan, Shinanomachi Rengakan Bldg, 2001; Tsukuba, Japan, Shinanomachi Rengakan Bldg, 2001; pp 7-21

20、.4.張平, 熱噴涂材料. 國防工業出版社: 北京, 2006.5.Ouyang, J. H.; Sasaki, S., Microstructure and tribological characteristics of ZrO2Y2O3 ceramic coatings deposited by laser-assisted plasma hybrid spraying. Tribology International 2002, 35 (4), 255-264.6.劉鳳嶺, 面向21世紀的熱障涂層結構設計. 材料保護 2000, (1), 86-88.7.Martena, M.; Bot

21、to, D.; Fino, P., et al., Modelling of TBC system failure: Stress distribution as a function of TGO thickness and thermal expansion mismatch. Engineering Failure Analysis 2006, 13 (3), 409-426.8.Rabiei, A.; Evans, A. G., Failure mechanisms associated with the thermally grown oxide in plasma-sprayed

22、thermal barrier coatings. Acta Materialia 2000, 48 (15), 3963-3976.9.Busso, E. P.; Evans, H. E.; Qian, Z. Q., et al., Effects of breakaway oxidation on local stresses in thermal barrier coatings. Acta Materialia 2010, 58 (4), 12421251.10.Rongved, L. In Force interior to one of two joined semi-infini

23、te solids, 2nd Midwestern Conference on Solid Mechanics, 1955; 1955; pp 1-13.11.Mindlin, R. D., Force at a point in the interior of a semi-innite solid. Journal of Physics 1936, 7, 195-202.12.Dundurs, J.; Hetenyi, M., Transmission of force between two semi-infinite solids. ASME, Journal of Applied M

24、echanics 1965, 32, 671-674.13.Bogy, D. B., Edge bonded dissimilar orthogonal elastic wedges under normal and shear loadings. Journal of Applied Mechanics 1968, 35, 146-154.14.Burmister, D. M., The General Theory of Stresses and Displacements in Layered Soil Systems. II. Journal of Applied Physics 19

25、45, 16 (3), 126-127.15.Burmister, D. M., The General Theory of Stresses and Displacements in Layered Soil Systems. III. Journal of Applied Physics 1945, 16 (5), 296-302.16.Burmister, D. M., The General Theory of Stresses and Displacements in Layered Systems. I. Journal of Applied Physics 1945, 16 (2

26、), 89-94.17.丁皓江, 橫觀各向同性彈性力學. 浙江大學出版社: 杭州, 1997; p 240.18.Fares, N.; Li, V. C., General image method in a plane-layered elastostatic medium. Journal of Applied Mechanics 1988, 55, 781-785.19.Yu, H. Y.; Sanday, S. C. In Elastic fields in joined half-spaces due to nuclei of strain, Proceedings of The R

27、oyal Society of London, 1991; 1991; pp 503-519.20.Walpole, L. J., An elastic singularity in joined half-spaces. International Journal of Engineering Science 1996, 34 (6), 629-638.21.Bahar, L. Y., A state space approach to elasticity. Journal of the Franklin Institute 1975, 299 (1), 33-41.22.Rao, N.

28、S. V. K.; Das, Y. C., A Mixed Method in Elasticity. ASME Journal of Applied Mechanics 1977, 44, 51-56.23.Bahar, L. Y., Transfer Matrix Approach to Layered Systems. Journal of the Engineering Mechanics Division 1972, 98, 1159-1172.24.Chen, W. T., Computation of stresses and displacements in a layered

29、 elastic medium. International Journal of Engineering Science 1971, 9, 775-800.25.Seale, S. H.; Kausel, E., Point Loads in Cross-Anisotropic Layered Halfspaces. Journal of Engineering Mechanics 1989, 115 (3), 509-542.26.Small, J. C.; Booker, J. R., Finite layer analysis of layered elastic materials

30、using a flexibility approach. Part 2 - Circular and rectangular loadings. International Journal for Numerical Methods in Engineering 2005, 23, 959-978.27.Kennedy, F. E.; Ling, F. F., Elasto-plastic indentation of a layered medium. Journal of Engineering Materials and Technology 1974, 96, 97-103.28.K

31、ral, E. R.; Komvopoulos, K., Three-Dimensional Finite Element Analysis of Subsurface Stresses and Shakedown Due to Repeated Sliding on a Layered Medium. Journal of Applied Mechanics 1996, 63 (4), 967-973.29.Kral, E. R.; Komvopoulos, K., Three-Dimensional Finite Element Analysis of Surface Deformatio

32、n and Stresses in an Elastic-Plastic Layered Medium Subjected to Indentation and Sliding Contact Loading. Journal of Applied Mechanics 1996, 63 (2), 365-375.30.Kral, E. R.; Komvopoulos, K., Three-Dimensional Finite Element Analysis of Subsurface Stress and Strain Fields Due to Sliding Contact on an

33、Elastic-Plastic Layered Medium. Journal of Tribology 1997, 119 (2), 332-341.31.Djabella, H.; Arnell, R. D., Finite Element Analysis of the Contact Stresses in Elastic Coating /Substrate Under Normal and Tangential Load. Thin Solid Films 1993, 223 (1), 87-97.32.Peng, W.; Bhushan, B., A Numerical Thre

34、e-Dimensional Model for the Contact of Layered Elastic/Plastic Solids With Rough Surfaces by a Variational Principle. Journal of Tribology 2001, 123 (2), 330-342.33.Tian, H.; Saka, N., Finite Element Analysis of an Elastic-Plastic Two-Layer Half-Space Sliding Contact. Wear 1991, 148 (2), 261-285.34.

35、Busso, E. P.; Lin, J.; Sakurai, S., A mechanistic study of oxidation-induced degradation in a plasma-sprayed thermal barrier coating system.: Part II: Life prediction model. Acta Materialia 2001, 49 (9), 15291536.35.Liu, Y.; Persson, C.; Wigren, J., Experimental and numerical life prediction of ther

36、mally cycled thermal barrier coatings. Journal of Thermal Spray Technology 2004, 13 (3), 415-424.36.Renusch, D.; Schütze, M., Measuring and modeling the TBC damage kinetics by using acoustic emission analysis. Surface and Coatings Technology 2007, 202 (4-7), 740744.37.Meier, S. M.; Sheffler, K. D.; Nissley, D. M. Thermal barrier coating life prediction model development, phase 2. Final report 1991.二、研