劉兆月，司良英，馬家驥，黎先浩，楊平，滕仁昊，王現(xiàn)輝

基于Stoney公式適用性分析——取向硅鋼絕緣涂層張應(yīng)力的計算

劉兆月1,2，司良英2，馬家驥3，黎先浩2，楊平1，滕仁昊2，王現(xiàn)輝2

（1.北京科技大學(xué) 材料科學(xué)與工程學(xué)院，北京 100083；2.首鋼智新遷安電磁材料有限公司，河北 遷安 064404；3.北京首鋼股份有限公司，河北 遷安 064404）

改善涂層張應(yīng)力狀態(tài)，進而提升取向硅鋼的電磁性能。用掃描電鏡、輝光光譜儀分析了添加磷酸鹽的膠體二氧化硅絕緣涂層的厚度和微觀結(jié)構(gòu)，采用靜態(tài)拉伸試驗得到了取向硅鋼沿不同位向的力學(xué)特性。在此基礎(chǔ)上，根據(jù)取向硅鋼各向異性的特點和絕緣涂層的實際狀態(tài)，用Stoney公式對取向硅鋼涂層應(yīng)力的適用性進行了討論。取向硅鋼的各向異性不滿足Stoney公式的適用條件，實際測試的取向硅鋼的力學(xué)性能結(jié)果表明，取向硅鋼符合正交對稱各向異性材料的一般規(guī)律，據(jù)此給出了涂層對取向硅鋼產(chǎn)生的張應(yīng)力的計算公式，利用該公式對不同涂覆工藝條件下的膠體二氧化硅絕緣涂層的張應(yīng)力進行了計算，涂層厚度在1~2 μm范圍內(nèi)，雙面張應(yīng)力的計算值為4~11 MPa。文中張應(yīng)力的計算公式不僅適用于取向硅鋼，同樣適用于其他基體為各向異性的類似情況。通過張應(yīng)力的計算可知，涂層厚度與張應(yīng)力呈正相關(guān)，張應(yīng)力的增大將會促使取向硅鋼的鐵損降低。因此，開發(fā)大張力的絕緣涂層，是進一步改善取向硅鋼磁性能的有效途徑之一。

取向硅鋼；絕緣涂層；張應(yīng)力；Stoney公式；各向異性

作為電力和電子行業(yè)不可或缺的一種重要的軟磁材料，取向硅鋼廣泛地應(yīng)用在變壓器等電力設(shè)施中。因具備鋒銳的Goss({110}<001>)織構(gòu)，取向硅鋼具有高磁感、低鐵損的優(yōu)勢[1-2]。隨著節(jié)能降耗、淘汰落后產(chǎn)能的需求不斷提升，降低取向硅鋼的鐵損已經(jīng)成為工業(yè)生產(chǎn)的核心目標(biāo)之一。為達到這一目的，不僅需要采用合理的工藝參數(shù)優(yōu)化材料本身的組織特性，同時也可以沿取向硅鋼軋向施加一定的張應(yīng)力來進一步降低鐵損[3-8]。在實際生產(chǎn)過程中，通過各種方式使取向硅鋼表面覆蓋應(yīng)力層（硅酸鎂底層、應(yīng)力絕緣涂層及其他涂層）就可以實現(xiàn)施加張應(yīng)力的目的，取向硅鋼表面的應(yīng)力絕緣涂層不僅能夠有效降低鐵損，還可以改善耐蝕性、提高可加工性及降低變壓器的噪音，從而提升產(chǎn)品的綜合使用性能[9-15]。分析涂層對基體的張應(yīng)力，有助于指導(dǎo)現(xiàn)場工藝優(yōu)化、提升產(chǎn)品性能和開發(fā)新型涂層。因此，需要采用恰當(dāng)?shù)姆椒▉頊y試硅鋼基體受到來自涂層的張應(yīng)力。

相對于其他測試方法，曲率法測量范圍廣、方法簡單，對基體及涂層的種類和狀態(tài)的要求較低，避免了因涂層局部缺陷引起的測量偏差大等問題[16]。1909年，Stoney[17]推導(dǎo)并提出了用于曲率法測量涂覆層殘余應(yīng)力的計算公式，并用此公式計算了電沉積鎳薄膜的殘余應(yīng)力。該公式自提出以來得到了廣泛應(yīng)用，已有許多學(xué)者以此為基礎(chǔ)，計算了包括半導(dǎo)體、陶瓷、電解沉積金屬及氧化膜等在內(nèi)的涂（覆）層殘余應(yīng)力數(shù)據(jù)[18-23]。曾有學(xué)者采用Stoney公式對陶瓷膜對基體的張應(yīng)力進行分析[15]，得到了涂層的張應(yīng)力數(shù)據(jù)。然而，Stoney公式是否適用于取向硅鋼涂層尚需細致分析。另外，對絕緣涂層產(chǎn)生的應(yīng)力進行計算和分析的相關(guān)研究報道較少。

本文以Stoney公式[17-18]為基礎(chǔ)，結(jié)合取向硅鋼及涂層特性，對公式進行了適用性的討論，并給出了適用于取向硅鋼的計算公式。以此為基礎(chǔ)，對取向硅鋼絕緣涂層的應(yīng)力進行了分析，簡要討論了張應(yīng)力與涂層厚度之間的關(guān)系，并分析了張應(yīng)力對鐵損的影響。

1 試驗

1.1 樣品制備

本文的研究對象為0.27 mm厚度的3%Si-Fe高磁感取向硅鋼，經(jīng)過熱軋、冷軋、脫碳滲氮、涂覆氧化鎂及高溫退火后，發(fā)生二次再結(jié)晶，在表面生成一層硅酸鎂底層，然后涂覆應(yīng)力絕緣涂層，最終獲得表面狀態(tài)良好的高磁感取向硅鋼。絕緣涂層以膠體SiO2為主，并加入了磷酸鹽等無機鹽類。分別采用5組不同涂布量參數(shù)，在同樣的燒結(jié)條件下，在樣品表面涂覆了應(yīng)力絕緣涂層。

1.2 曲率法測試原理

采用曲率法進行張應(yīng)力測試，樣品僅存在一側(cè)涂層的情況下獲得樣品彎曲的曲率。取尺寸為30 mm× 300 mm的樣品（長邊方向即為應(yīng)力測試方向，本文為軋向），單側(cè)存在絕緣涂層時發(fā)生彎曲，此時樣品軋面可以近似看作圓柱側(cè)面的一部分，長邊則可以認為是圓柱底面圓的部分弧長。通過簡單的幾何學(xué)就能計算出這個底面圓的半徑，可以得到所需的曲率半徑

1.3 性能測試及組織分析

為保證準(zhǔn)確性，以下測量樣品均取自待測樣品周圍部位。

1）利用靜態(tài)拉伸試驗實測無絕緣涂層的樣品（表面存在硅酸鎂底層）的彈性模量和泊松比數(shù)據(jù)。

2）利用掃描電鏡對取向硅鋼表面涂層的形貌進行分析。

3）利用輝光光譜儀測量絕緣涂層的厚度，其分析區(qū)域大小是直徑為3 mm的圓，相對于微觀方法更具代表性。由于取向硅鋼還存在硅酸鎂底層，涂層狀態(tài)較為復(fù)雜，以成品取向硅鋼的分析結(jié)果為例（見圖1）對分析方法進行說明：在距離表面1.3 μm處，Mg元素和Si元素含量發(fā)生明顯變化，絕緣涂層中的其他元素逐漸消失，則可認為此處是硅酸鎂底層與絕緣涂層的交界位置，即本例中絕緣涂層的厚度大約為1.3 μm。

4）利用愛潑斯坦方圈對硅鋼的鐵芯損耗（1.7/50）進行測試。

2 結(jié)果與討論

2.1 Stoney公式的原理及演變

一般認為，涂層的殘余應(yīng)力起源于基體和涂層的熱膨脹失配作用，作用原理如圖2所示。一般情況下，樣品上下同時存在相近厚度的涂層，由于受力平衡因而在宏觀上表現(xiàn)為平整的外觀。但當(dāng)樣品僅有一側(cè)存在應(yīng)力涂層時，樣品就會發(fā)生彎曲，如圖3所示。

圖2 涂層對取向硅鋼應(yīng)力作用的局部示意圖

圖3 僅有單面涂層情況下取向硅鋼的彎曲情況示意圖

Stoney[17]于1909年提出了計算薄膜殘余應(yīng)力的公式，見式(1)。該公式利用彎曲形變的曲率半徑“”來表征薄膜的殘余應(yīng)力大小，通過基體和薄膜之間力的關(guān)系建立等式。考慮到薄膜應(yīng)力的不均勻性（見圖2），實際上Stoney公式計算得到的是應(yīng)力的平均值。

式中：f表示涂層的殘余應(yīng)力，f為涂層厚度，s為基體的彈性模量，s為基板厚度，為基體彎曲的曲率。

Stoney公式自提出以來便得到了廣泛應(yīng)用。目前，有很多關(guān)于對Stoney公式在不同應(yīng)用條件下適用性及修正的討論[19-20,24-25]：這些公式考慮了包括基體厚度、基體的彎曲程度等在內(nèi)的因素。其中值得一提的是，關(guān)于雙軸應(yīng)力的修正[17]。Stoney公式即式(1)僅僅考慮了單軸應(yīng)力，然而，實際樣品中，計算樣品受到的應(yīng)力均應(yīng)考慮雙軸狀態(tài)：除樣品厚度外的長（）和寬（）方向。

在基體受到雙軸應(yīng)力的情況下，當(dāng)基體和薄膜均為各向同性時，涂層對基體產(chǎn)生的張應(yīng)力與基體應(yīng)變之間的關(guān)系可以表示為：

式中：s為基體因涂層受到的應(yīng)力，s為基體的泊松比，s則表示基體發(fā)生的應(yīng)變。由于薄膜的各向同性，則可認為各方向的張應(yīng)力值相等，根據(jù)式(2)，用s/(1-s)代替(1)式中的s，Stoney公式就變?yōu)楦毡榈男问?3)[18]：

2.2 Stoney公式適用性討論

Stoney在提出公式時就認為該公式在涂層非常薄的情況下才能夠適用[17]。事實上，應(yīng)用Stoney公式必須滿足以下假設(shè)[24]：（i）與橫向尺寸相比，薄膜和襯底厚度都很??；（ii）薄膜厚度遠小于襯底厚度；（iii）基材均勻、各向同性、線性彈性好，薄膜材料各向同性；（iv）襯底邊緣附近的邊緣效應(yīng)不重要，所有物理量在與界面平行的位置發(fā)生變化時都不變；（v）厚度方向上的所有應(yīng)力成分在整個材料中消失，且應(yīng)變和旋轉(zhuǎn)無窮小。因此，若要采用Stoney公式計算涂層的殘余應(yīng)力或基體受到的應(yīng)力，必須結(jié)合取向硅鋼及其涂層的實際情況對公式的適用性進行討論，具體如下所述。

采用掃描電鏡分析取向硅鋼涂層的典型結(jié)構(gòu)。實際上，取向硅鋼的絕緣涂層下方還存在一層硅酸鎂底層，因此對涂覆前后樣品的狀態(tài)分別進行分析。圖4為取向硅鋼為涂覆絕緣涂層前截面的微觀狀態(tài)，可以看到硅酸鎂底層的顆粒緊密地附著在基體上，厚度大約為1 μm，平整度較差。圖5展示取向硅鋼成品絕緣涂層的微觀形貌和結(jié)構(gòu)，可以看到絕緣涂層緊密地附著在硅酸鎂底層上，二者沒有明顯的界限，涂層狀態(tài)均勻致密，涂覆涂層后樣品表面變得平整。

圖4 硅酸鎂底層的截面形貌

根據(jù)涂層的結(jié)構(gòu)、成分和特點可知，涂層內(nèi)部的殘余應(yīng)力是各向同性的。盡管取向硅鋼具有鋒銳的Goss({110}<001>)織構(gòu)，存在明顯的各向異性，但因晶體結(jié)構(gòu)的高對稱性（體心立方α-Fe），晶粒尺寸較大（與單晶類似），因此可認為熱膨脹系數(shù)不存在各向異性[26-27]，故可認為最終作用在硅鋼基體上各方向的張應(yīng)力為各向同性。取向硅鋼表面的涂層厚度在1~4 μm左右（參看圖1和圖5），而鋼板的板厚則為0.20~0.30 mm，二者特點滿足Stoney公式對于涂層性質(zhì)和厚度的假設(shè)條件（ii）。而關(guān)于基體的假設(shè)條件中，當(dāng)被測試樣品在狀態(tài)良好、尺寸恰當(dāng)?shù)臈l件下，條件（i）、（iv）和（v）容易得到滿足，但是取向硅鋼存在明顯的各向異性。圖6為采用靜態(tài)拉伸方法獲得的（部分）力學(xué)性能數(shù)據(jù)，其中圖6a是沿TD方向拉伸的縱向和橫向的應(yīng)力-應(yīng)變曲線，圖6b則為沿RD方向拉伸的縱向和橫向的應(yīng)力-應(yīng)變曲線，結(jié)果顯示，橫向和縱向的彈性模量和泊松比有明顯區(qū)別。因此，取向硅鋼不滿足條件（iii），無法直接使用式(3)來進行計算，需要針對取向硅鋼的特點，結(jié)合Stoney公式的原理，進行更深入的分析和推導(dǎo)。

圖6 沿兩個方向測試的取向硅鋼靜態(tài)拉伸曲線

根據(jù)廣義的胡克定律，具有彈性正交對稱性的各向異性材料各應(yīng)變分量與應(yīng)力分量之間的關(guān)系如下：

已有學(xué)者[28]對Si單晶材料表面薄膜殘余應(yīng)力的計算公式進行推導(dǎo)，公式中采用了材料某些特定方向的S（彈性柔度系數(shù)，第一個下標(biāo)為應(yīng)變方向，第二個下標(biāo)為樣品受到的應(yīng)力方向，其他力學(xué)性能參數(shù)的下標(biāo)含義相同）來分別表達Si(001)單晶和Si(111)單晶表面膜層的殘余應(yīng)力，公式的表達形式與基體的晶體取向及相應(yīng)的力學(xué)參數(shù)密切相關(guān)，形式多樣且較為復(fù)雜。

由于取向硅鋼基體所受應(yīng)力的計算過程中不涉及剪切應(yīng)力，根據(jù)式(4)，方向或者方向張應(yīng)力的計算公式為：

涂層的張應(yīng)力主要作用在基體的、方向，因此計算過程中可不考慮厚度方向的σ和ε。根據(jù)彈性柔順系數(shù)的定義，以方向為例，式(5)可變?yōu)槭?7)的形式，方向可以類推。其中，基體受到的張應(yīng)力為各向同性，故σ與σ相等，故式(7)統(tǒng)一表達為s。

式(7)可做如下理解，方向的應(yīng)變由兩部分組成：僅由方向受到的應(yīng)力引起的應(yīng)變與方向的應(yīng)力導(dǎo)致方向的應(yīng)變。在式(7)中，ν為方向受力引起的方向發(fā)生應(yīng)變，即沿方向進行拉伸試驗所得到的泊松比數(shù)據(jù)。

根據(jù)式(4)，正交對稱性材料的彈性柔順系數(shù)S=S，結(jié)合彈性柔順系數(shù)的定義，可得：

通過測量得到的未涂覆絕緣涂層樣品的彈性模量和泊松比數(shù)據(jù)結(jié)果如表1所示。

表1 取向硅鋼力學(xué)性能數(shù)據(jù)

Tab.1 Mechanical properties of grain-oriented silicon steel

在彈性各向同性的材料中，鋼鐵材料的泊松比一般在0.20~0.30范圍內(nèi)。但因取向硅鋼具有各向異性，泊松比數(shù)據(jù)的“范圍”發(fā)生了變化，已有許多學(xué)者證明[29-31]：各向異性材料的泊松比的限制條件與各向同性材料有很大的不同。LEMPRIERE[32]推導(dǎo)出了具有各向異性的立方晶體材料的泊松比受彈性模量的比例限制，如式(9)所示。

表1中測試的數(shù)據(jù)符合式(8)和(9)的條件，這說明上述關(guān)于力學(xué)性能的分析是適用于取向硅鋼的。在式(8)中，分別取、為、方向，代入式(7)即可得到基體沿方向或方向受到張應(yīng)力的計算公式：

式中：s表示基體受到來自涂層的殘余應(yīng)力，s為基體厚度，E為基體沿方向的彈性模量，則表示基體沿方向彎曲的曲率半徑。雖然基體受到的應(yīng)力為各向同性的，但是由于力學(xué)性能的各向異性，沿不同方向得到的曲率值卻各不相同。就本文來說，軋向的曲率半徑明顯小于橫向，且更易測量。

本文在Stoney公式和公式(3)的基礎(chǔ)上，引入了各向異性材料沿不同位向的力學(xué)性能，修正了公式(3)，解決了取向硅鋼這種各向異性材料受到來自基體的張應(yīng)力計算的問題。對于大多數(shù)材料來說，能夠通過靜態(tài)拉伸試驗或其他測量方法獲得力學(xué)性能數(shù)據(jù)，測量較為簡單，只要將力學(xué)性能數(shù)據(jù)代入式(10)中就能夠得到基體的張應(yīng)力數(shù)據(jù)。因此，式(10)不僅適用于計算取向硅鋼表面膜層的殘余應(yīng)力，還可以用于其他條件類似的材料。當(dāng)然，使用公式的前提是必須滿足Stoney公式的適用條件。

綜上，計算取向硅鋼絕緣涂層沿某一方向?qū)w產(chǎn)生的張應(yīng)力s的計算公式為(10)，需要注意的是，本公式考慮的情況是單側(cè)涂層的張應(yīng)力，當(dāng)取向硅鋼兩側(cè)均存在相同的涂層，因力的方向一致，可認為基體所受的應(yīng)力為雙側(cè)應(yīng)力之和。

2.3 張應(yīng)力與涂層厚度及鐵損之間的關(guān)系

張應(yīng)力可以明顯改變?nèi)∠蚬桎摰蔫F損。研究證明[3-5]，張應(yīng)力能夠細化磁疇結(jié)構(gòu)，消除匕首疇和90°閉合疇等結(jié)構(gòu)，在一定范圍內(nèi)，對取向硅鋼施加張應(yīng)力將降低鐵損，鐵損的下降幅度和張應(yīng)力數(shù)值呈線性相關(guān)。在實際生產(chǎn)過程中，人們利用涂層對基體產(chǎn)生張應(yīng)力來實現(xiàn)鐵損降低的目的，如硅酸鎂底層、TiN類陶瓷涂層、硅溶膠-磷酸鹽絕緣涂層等[8-13]，涂層種類、涂布方式和涂層涂布量都將對鐵損的降低幅度產(chǎn)生影響，因此可以通過改變涂層狀態(tài)來進一步優(yōu)化取向硅鋼的鐵損，這對高磁感取向硅鋼來說更重要。

本文計算得到了曲率半徑，測試了樣品的力學(xué)性能，代入式(10)中計算得到樣品受到來自絕緣涂層的張應(yīng)力數(shù)據(jù)，每組試驗樣品分別測試了五組數(shù)據(jù)，具體見表2，可以看到同組樣品的張應(yīng)力計算數(shù)據(jù)較為穩(wěn)定，一致性較好。

表2 張應(yīng)力測試結(jié)果

Tab.2 Test result of tensile stress

圖7為涂層厚度和基體受到的張應(yīng)力數(shù)據(jù)（雙面涂層的作用），可知在同樣的燒結(jié)工藝條件下，隨著涂層厚度的增大，基體受到的張應(yīng)力增加。圖8為基體受到來自涂層的張應(yīng)力對樣品鐵損的影響，其中，縱坐標(biāo)的數(shù)據(jù)為同一樣品涂覆應(yīng)力絕緣涂層后鐵損降低的數(shù)值。從圖8中可以看到，隨著張應(yīng)力的增大，樣品的1.7/50的降低幅度越來越大，本文得到的張應(yīng)力與鐵損的數(shù)值關(guān)系與文獻數(shù)據(jù)[6,8]基本相當(dāng)。雖然隨著張應(yīng)力的升高鐵損不斷得到優(yōu)化，但是對于本文中所采用的涂層類別和涂覆工藝下，絕緣涂層厚度超過1.8 μm時，涂層易于脫落，附著性變差，厚度達到2.0 μm時，附著性極差，已經(jīng)不滿足實際使用需求（隨著涂層成分類別的變化，涂層厚度的極限值可能發(fā)生變化）。因此，該種類涂層的厚度應(yīng)控制在合理的范圍內(nèi)，同時滿足磁性能、使用性能及生產(chǎn)成本的要求。另外，在實際生產(chǎn)中，應(yīng)該在滿足使用需求的情況下，開發(fā)張力更大的涂層，進一步提升產(chǎn)品的磁性能。

圖7 絕緣涂層厚度與張應(yīng)力的關(guān)系

圖8 涂層的張應(yīng)力對鐵損的影響

3 結(jié)論

1）本文討論了Stoney公式對取向硅鋼的適用性，實際測試了取向硅鋼的力學(xué)性能，認為取向硅鋼符合正交對稱各向異性材料的一般規(guī)律，在此基礎(chǔ)上提出了取向硅鋼涂層應(yīng)力的計算公式，解決了各向異性Stoney公式應(yīng)用中的問題，此公式還可用于其他滿足Stoney公式適用條件且基體為各向異性的情況。

2）本文的絕緣涂層為添加了磷酸鹽的膠體二氧化硅涂層，計算了其對取向硅鋼產(chǎn)生的張應(yīng)力，涂層厚度在1~2 μm范圍內(nèi)，雙面張應(yīng)力的計算值為4~11 MPa。

3）隨著絕緣涂層厚度的增大，涂層給予基體的張應(yīng)力增加，而隨著張應(yīng)力的增加，鐵損逐漸降低。在滿足實際使用需求的前提下，開發(fā)大張力的涂層，是優(yōu)化取向硅鋼性能的途徑之一。

[1] MATSUO M. Texture control in the production of grain oriented silicon steels[J]. Transactions of the iron & steel institute of Japan, 1989, 29(10): 809-827.

[2] 黎先浩, 孟小濤, 趙鵬飛, 等. 高磁感取向硅鋼研發(fā)現(xiàn)狀與展望[J]. 中國冶金, 2019, 29(1): 1-7. LI X H, MENG X T, ZHAO P F, et al. Present status and future prospect of high permeability grain-oriented silicon steel[J].China metallurgy,2019, 29(1): 1-7.

[3] GNIEWEK J J. Effects of tensile stress on the domain structure in grain-oriented 3.25% silicon steel[J]. Journal of applied physics, 1963, 34(12): 3618-3622.

[4] HOUZE G L. Effect of longitudinal tensile stress on the static and 60 Hz domain structures of grain-oriented siliconsteel[J]. Journal of applied physics, 1969, 40(3): 1090-1091.

[5] PERRYMAN R. Effects of tensile and compressive stress on 3% grain-oriented silicon iron[J]. Journal of physics D: Applied physics, 1975, 8(15): 1901-1909.

[6] OVERSHOTT K J, Foot G. The effect of tensile stress on the power loss of 3% grain-oriented silicon-iron[J]. IEEE transactions on magnetics, 1982, 18(6): 1496-1498.

[7] HEGAB S, Al-Mutawa S. Strain-magnetization properties and domain structures of silicon steel sheets[J]. Electrical engineering in Japan, 2010, 149(2): 10-17.

[8] MOSES A J, Pike E C. Influence of surface coating on transverse magnetostriction in silicon iron[J]. Electrical engi-neers proceedings of the institution, 1977, 124(6): 583-585.

[9] INOKUT Y. Grain oriented silicon steel sheet with a ceramicfilm characterized by ultra-low iron loss[J]. Vacuum, 1996, 47(6-8): 857-862.

[10] PUZHEVICH R V, KORZUNIN G S, KARENINA L S, et al. Mutual influence of electrical insulating coatings and structural barriers on the level of magnetic losses in aniso-tropic electrical sheet steel[J]. Physics of metals & metallo-graphy, 2009, 108(6): 543-547.

[11] PUZHEVICH R B, DUKHNOV A G, SAMSIKOV E A, et al. Effect of the ground layer on the magnetic properties of highly permeable silicon steel[J]. Metal science & heat treatment, 1986, 28(6): 435-439.

[12] SHIMANAK H, ICHIDA T, KOBAYASHI S, et al. Effect of several stress inducing coatings on magnetostriction pro-perty of a high induction 3% silicon steel[J]. IEEE transac-tions on magnetics, 2003, 15(6): 1595-1597.

[13] SHINGAKI Y, OKABE S. Influence of crystal orientation on magnetic properties in 3% silicon steel with ultra-high tension TiN coating[J]. IEEE transactions on magnetics, 2012, 48(4): 1469-1472.

[14] INOKUTI Y, SUZUKI K. Relation between iron loss and tensile stress induced in silicon steel sheets by TiN, TiCN and TiC coatings[J]. Journal of the vacuum society of Japan, 2009, 46(3): 329-332.

[15] 汪伊麗, 梁天權(quán), 沈曉明, 等. 熱障涂層界面應(yīng)力測試研究進展[J]. 材料導(dǎo)報, 2016, 30(S2): 495-498.WANG Y L, LIANG T Q, SHEN X M, et al. Research pro-gress of interfacial stress measurement in thermal barrier coatings[J]. Materials review, 2016, 30(S2): 495-498.

[16] STONEY G G. The tension of metallic films deposited by electrolysis[J]. Proceedings of the Royal Society of London, 1909, 82(309): 40-43.

[17] JACCODINE R J, SCHLEGEL W A. Measurement of stra-ins at Si-SiO2interface[J]. Journal of applied physics, 1966, 37(6): 2429-2434.

[18] FLINN P A, GARDNER D S, NIX W D. Measurement and interpretation of stress in aluminum-based metalliza-tion as a function of thermal history[J]. IEEE transactions on electron devices, 1987, 34(3): 689-699.

[19] CHIU C C. Residual stresses in ceramic coatings as deter-mined from the curvature of a coated strip[J]. Materials science & engineering A, 1992, 150(1): 139-148.

[20] ARDIGO M R, AHMED M, BESNARD A L. Stoney for-mula: Investigation of curvature measurements by optical profilometer[J]. Advanced materials research, 2014, 996(8): 361-366.

[21] SHIRI S, ASHTIJOO P, ODESHI A, et al. Evaluation of Sto-ney equation for determining the internal stress of DLC thin films using an optical profiler[J]. Surface & coatings technology, 2016, 308(11): 98-100.

[22] WEN J, WEI Y, CHENG Y T. Examining the validity of Stoney-equation for, in-situ, stress measurements in thin film electrodes using a large-deformation finite-element proce-dure[J]. Journal of power sources, 2018, 387(5): 126-134..

[23] KIM J S, PAIK K W, Oh S H. The multilayer-modified Stoney′s formula for laminated polymer composites on a silicon substrate[J]. Journal of applied physics, 1999, 86(10): 5474-5479.

[24] FREUND L B, FLORO J A, CHASON E. Extensions of the Stoney formula for substrate curvature to configura-tions with thin substrates or large deformations[J]. Applied physics letters, 1999, 74(14): 1987-1989.

[25] KLEIN C A. How accurate are Stoney′s equation and recent modifications[J]. Journal of applied physics, 2000, 88(9): 5487-5489.

[26] SIROTIN Y I, SHASKOL'SKAYA M P. Fundamentals of crystal physics[M]. Moscow: Science, 1979.

[27] ZAMKOVSKAYA A, MAKSIMOVA E. Anisotropic ther-mal expansion in crystals of different categories[j]. Jour-nal of physics conference series, 2015, 58: 661.

[28] JANSSEN G C A M, ABDALLA M M, KEULEN F V, et al. Celebrating the 100th anniversary of the Stoney equation for film stress: Developments from polycrystalline steel strips to single crystal silicon wafers[J]. Thin solid films, 2009, 517(6): 1858-1867.

[29] HAYES M, SHUVALOV A. On the extreme values of Young′s Modulus, the Shear Modulus, and Poisson′s ratio for cubic materials[J]. Journal of applied mechanics, 1998, 65(3): 786-787.

[30] TING T C T, CHEN T Y. Poisson′s ratio for anisotropic elastic materials can have no bounds[J]. Quarterly journal of mechanics & applied mathematics, 2005, 58(1): 73-82.

[31] WALTON R I, SMITH C W, LETHBRIDGE Z A D, et al. Elastic anisotropy and extreme Poisson′s ratios in single crystals[J]. Acta materialia, 2010, 58(19): 6444-6451.

[32] LEMPRIERE B M. Poisson′s ratio in orthotropic mater-ials[J]. AIAA journal, 1968, 6(11): 2226-2277.

Based on Applicability Analysis of Stoney Equation: Calculation of Tensile Stress Induced in Grain-oriented Silicon Steel by Insulation Coating

1,2,2,3,2,1,2,2

(1.School of Materials Science and Engineering, University of Science and Technology Beijing, Beijing 100083, China; 2.Shougang Zhixin Qian'an Electromagnetic Material Co. Ltd, Qian'an 064404, China; 3.Beijing Shougang Co., Ltd, Qian'an 064404, China)

The work aims to optimize the tensile stress of the coating, so as to improve the magnetic properties of grain-oriented silicon steel. The thickness and microstructure of colloidal silica insulation coating with phosphate were analyzed by SEM and GDS and the mechanical properties of substrate along different direction were acquired by tensile test. Based on the anisotropy of grain-oriented silicon steel and the real state of coating, the applicability of Stoney equation to the grain-oriented silicon steel coating was discussed. The anisotropy of grain-oriented silicon steel did not meet the applicable conditions of Stoney equation and the mechanical properties of grain-oriented silicon steel tested in practice showed that the grain-oriented silicon steel conformed to the general law of orthotropic material. Based on that, the equation was given to calculate the tensile stress induced in the substrate of oriented silicon steel by the colloidal silica insulation coating under different technology parameters. The tensile stress was 4~11 MPa when the coating thickness was within 1~2 microns. The equation for calculating the tensile stress is not only applicable to grain-oriented silicon steel, but also applicable to similar cases where other substrates are anisotropic. Through calculation of tensile stress, the stress and thickness of the insulation coating have a positive correlation and the iron loss reduce with increasing of the stress induced in substrate by the coating. Therefore, the research and development of insulation coating with higher tensile stress is the effective method to improve the magnetic properties of grain-oriented silicon steel.

grain-oriented silicon steel; insulation coating; tensile stress; Stoney equation; anisotropy

2019-08-10；

2019-11-03

LIU Zhao-yue (1987—), Female, Doctoral candidate, Research focus: research and development of electrical steel.

楊平（1959—），男，博士，教授，博士研究生導(dǎo)師，主要研究方向為晶體材料的形變、再結(jié)晶、相變及織構(gòu)。郵箱：yangp@mater.ustb.edu.cn

Corresponding author：YANG Ping (1959—), Male, Doctor, Professor, Doctoral supervisors, Research focus: deformation, recrystallization, phase transition and texture of crystalline materials. e-mail: yangp@mater.ustb.edu.cn

劉兆月, 司良英, 馬家驥, 等. 基于Stoney公式適用性分析——取向硅鋼絕緣涂層張應(yīng)力的計算[J]. 表面技術(shù), 2020, 49(6): 146-151.

TG174.4;TB34

A

1001-3660(2020)06-0146-06

10.16490/j.cnki.issn.1001-3660.2020.06.017

2019-08-10；

2019-11-03

劉兆月（1987—），女，博士研究生，主要研究方向為電工鋼材料的研發(fā)。

LIU Zhao-yue, SI Liang-ying, MA Jia-ji, et al. Based on applicability analysis of stoney equation: calculation of tensile stress induced in grain-oriented silicon steel by insulation coating[J]. Surface technology, 2020, 49(6): 146-151.