預(yù)腐蝕對(duì)多部位損傷發(fā)生概率的影響

張建宇 袁 軍 時(shí)新紅

(北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院，北京100191)

多部位損傷(MSD，Multiple Side Damage)是指在同一個(gè)結(jié)構(gòu)中的相似細(xì)節(jié)在循環(huán)載荷作用下產(chǎn)生一群長(zhǎng)度相近、互相影響的宏觀相對(duì)小裂紋，或是一條長(zhǎng)裂紋向一群這樣的裂紋擴(kuò)展.這些裂紋的存在會(huì)顯著降低結(jié)構(gòu)的剩余強(qiáng)度，甚至?xí)斐蔀?zāi)難性的后果［1］.2010年11月，美國(guó)聯(lián)邦航空局頒布了最新的適航規(guī)章［2］，明確指出無(wú)論新設(shè)計(jì)的飛機(jī)還是在役飛機(jī)，都應(yīng)給出一個(gè)限制期，應(yīng)通過(guò)分析和試驗(yàn)表明結(jié)構(gòu)在這個(gè)期限內(nèi)不會(huì)發(fā)生MSD.

飛機(jī)結(jié)構(gòu)中存在大量的成排鉚釘孔結(jié)構(gòu)，其使用及停放過(guò)程中的預(yù)腐蝕問(wèn)題是不可避免的.關(guān)于疲勞載荷下的MSD發(fā)生問(wèn)題，文獻(xiàn)［3］提出了一種新的MSD發(fā)生概率計(jì)算方法，可以通過(guò)單一細(xì)節(jié)疲勞壽命的分布特征，估算出多細(xì)節(jié)結(jié)構(gòu)的MSD發(fā)生概率;但是預(yù)腐蝕影響下的MSD發(fā)生概率研究還未見(jiàn)相關(guān)文獻(xiàn)資料.因此研究這些相似細(xì)節(jié)結(jié)構(gòu)在預(yù)腐蝕后發(fā)生多部位損傷的概率具有重要的意義.

本文提出了一種新的預(yù)腐蝕后多細(xì)節(jié)結(jié)構(gòu)MSD發(fā)生概率公式，在單細(xì)節(jié)試驗(yàn)的基礎(chǔ)上得到預(yù)腐蝕影響系數(shù).將不經(jīng)過(guò)預(yù)腐蝕處理進(jìn)行的疲勞試驗(yàn)規(guī)定為常規(guī)試驗(yàn)，通過(guò)常規(guī)試驗(yàn)結(jié)果，結(jié)合預(yù)腐蝕影響系數(shù)，預(yù)測(cè)預(yù)腐蝕后任意細(xì)節(jié)數(shù)、任意壽命下的MSD發(fā)生概率，以及任意細(xì)節(jié)數(shù)、任意可靠度下MSD發(fā)生時(shí)的壽命，進(jìn)而分析預(yù)腐蝕對(duì)MSD發(fā)生概率的影響.

1 預(yù)腐蝕后MSD發(fā)生概率公式推導(dǎo)

在對(duì)預(yù)腐蝕處理后的多細(xì)節(jié)結(jié)構(gòu)MSD發(fā)生概率進(jìn)行計(jì)算時(shí)，根據(jù)文獻(xiàn)［4-5］中的結(jié)論，做出兩個(gè)假設(shè):①假設(shè)在預(yù)腐蝕處理后，單細(xì)節(jié)結(jié)構(gòu)壽命仍然符合對(duì)數(shù)正態(tài)分布.②假設(shè)單細(xì)節(jié)結(jié)構(gòu)對(duì)數(shù)疲勞壽命的分散性只與應(yīng)力水平有關(guān)，預(yù)腐蝕條件對(duì)疲勞壽命的分散性沒(méi)有影響.那么，在計(jì)算MSD發(fā)生概率PMSD時(shí)，無(wú)論是否有預(yù)腐蝕的影響，對(duì)數(shù)疲勞壽命標(biāo)準(zhǔn)差均取常規(guī)條件下相同應(yīng)力水平時(shí)裂紋萌生壽命的對(duì)數(shù)壽命標(biāo)準(zhǔn)差.

預(yù)腐蝕影響系數(shù)C定義［6］為

其中，NC50表示預(yù)腐蝕后，對(duì)應(yīng)應(yīng)力水平下的單細(xì)節(jié)疲勞試驗(yàn)的中值壽命;N50表示常規(guī)單細(xì)節(jié)疲勞試驗(yàn)的中值壽命.

根據(jù)中值壽命與對(duì)數(shù)壽命均值的關(guān)系，預(yù)腐蝕影響系數(shù)C可以通過(guò)對(duì)數(shù)壽命均值表示為

其中，μC表示預(yù)腐蝕后的單細(xì)節(jié)疲勞試驗(yàn)的對(duì)數(shù)壽命均值;μ表示單細(xì)節(jié)常規(guī)疲勞試驗(yàn)對(duì)數(shù)壽命均值.

預(yù)腐蝕后的單細(xì)節(jié)疲勞壽命N服從對(duì)數(shù)正態(tài)分布，則對(duì)數(shù)壽命lgN服從正態(tài)分布，即lg N～N(μC，σC).預(yù)腐蝕后單細(xì)節(jié)對(duì)數(shù)壽命的概率密度函數(shù)fC(x)就可以用正態(tài)分布的形式表示為

由式(2)可推出，預(yù)腐蝕后對(duì)數(shù)壽命均值為

根據(jù)假設(shè)②，fC(x)可以用常規(guī)試驗(yàn)單細(xì)節(jié)對(duì)數(shù)壽命均值μ、標(biāo)準(zhǔn)差σ及C表示:

假設(shè)飛機(jī)相似多細(xì)節(jié)結(jié)構(gòu)各細(xì)節(jié)相互獨(dú)立，相互之間不會(huì)發(fā)生影響.飛機(jī)相似多細(xì)節(jié)結(jié)構(gòu)各細(xì)節(jié)之間應(yīng)力差異很小，所處的工作環(huán)境相同，可以認(rèn)為這些細(xì)節(jié)承受相同大小的應(yīng)力水平，并且預(yù)腐蝕條件相同.從初始到循環(huán)數(shù)N時(shí)，各個(gè)細(xì)節(jié)裂紋萌生概率為

結(jié)合文獻(xiàn)［3］中常規(guī)試驗(yàn)MSD發(fā)生概率公式，提出預(yù)腐蝕后的PMSD計(jì)算公式為其中，N為對(duì)應(yīng)的疲勞壽命;n為細(xì)節(jié)數(shù).通過(guò)上述公式可以看出，只要得到了C，就可以通過(guò)單細(xì)節(jié)常規(guī)試驗(yàn)的結(jié)果，預(yù)測(cè)預(yù)腐蝕后任意細(xì)節(jié)數(shù)孔板壽命為N時(shí)的MSD發(fā)生概率，或任意細(xì)節(jié)數(shù)孔板某可靠度下發(fā)生MSD時(shí)的疲勞壽命.

2 預(yù)腐蝕影響系數(shù)C的測(cè)定

通過(guò)對(duì)7B04-T74δ2.0鋁合金板材進(jìn)行單細(xì)節(jié)的常規(guī)疲勞試驗(yàn)與預(yù)腐蝕后的疲勞試驗(yàn)，得到預(yù)腐蝕影響系數(shù)C.試驗(yàn)機(jī)為INSTRON-8803電液伺服材料試驗(yàn)機(jī).常規(guī)疲勞試驗(yàn)采用載荷控制的方法，疲勞載荷為等幅譜，波形為正弦波，加載頻率為 15 Hz，應(yīng)力比 R=0.1，應(yīng)力水平 S=87.5 MPa.

單細(xì)節(jié)試件是指只含有一個(gè)危險(xiǎn)部位，當(dāng)此部位產(chǎn)生疲勞裂紋或者斷裂時(shí)，認(rèn)為該試件發(fā)生破壞.但是為了節(jié)省資金及試驗(yàn)的時(shí)間可以采用串聯(lián)雙細(xì)節(jié)試件進(jìn)行單細(xì)節(jié)試驗(yàn).

所謂串聯(lián)雙細(xì)節(jié)試件是指在同一個(gè)試件中，串聯(lián)布置兩個(gè)細(xì)節(jié)，兩孔間的距離足夠遠(yuǎn)，使得加載時(shí)兩個(gè)細(xì)節(jié)受載荷情況完全獨(dú)立，在其中一個(gè)孔斷裂后，可對(duì)另一孔繼續(xù)加載至斷裂，從而得到兩個(gè)有效數(shù)據(jù)，節(jié)省了試驗(yàn)時(shí)間和試驗(yàn)成本.圖1為串聯(lián)雙細(xì)節(jié)試件形式及加載方向.

圖1 串聯(lián)雙細(xì)節(jié)孔板(單位:mm)

預(yù)腐蝕的方法為3.5%NaCl溶液連續(xù)浸泡240 h.預(yù)腐蝕按照《JB7901—1999金屬材料實(shí)驗(yàn)室均勻腐蝕全浸試驗(yàn)方法》進(jìn)行.預(yù)腐蝕后的疲勞試驗(yàn)同樣為載荷控制，等幅譜加載，波形為正弦波，應(yīng)力比R=0.1，加載頻率為15 Hz，應(yīng)力水平S=87.5 MPa.

單細(xì)節(jié)常規(guī)疲勞試驗(yàn)對(duì)數(shù)壽命均值與對(duì)數(shù)壽命標(biāo)準(zhǔn)差如表1所示.

預(yù)腐蝕處理后進(jìn)行疲勞試驗(yàn)，得到預(yù)腐蝕后對(duì)數(shù)壽命均值，結(jié)果如表2所示.

表1 單細(xì)節(jié)常規(guī)試驗(yàn)結(jié)果

表2 預(yù)腐蝕后單細(xì)節(jié)疲勞試驗(yàn)結(jié)果

通過(guò)上述單細(xì)節(jié)常規(guī)疲勞試驗(yàn)及單細(xì)節(jié)預(yù)腐蝕后疲勞試驗(yàn)結(jié)果，利用式(2)可得到經(jīng)過(guò)3.5%NaCl浸泡240 h后，應(yīng)力水平S=87.5 MPa所對(duì)應(yīng)的預(yù)腐蝕影響系數(shù):

3 MSD發(fā)生概率計(jì)算

3.1 多細(xì)節(jié)試驗(yàn)

多細(xì)節(jié)試驗(yàn)同樣采用7B04-T74δ2.0鋁合金板材，分別進(jìn)行三細(xì)節(jié)孔板常規(guī)環(huán)境疲勞試驗(yàn)及五細(xì)節(jié)孔板預(yù)腐蝕-疲勞試驗(yàn).試驗(yàn)中所用的試驗(yàn)機(jī)及加載方式均與單細(xì)節(jié)試驗(yàn)相同，應(yīng)力水平仍取S=87.5 MPa.三細(xì)節(jié)、五細(xì)節(jié)試驗(yàn)件尺寸分別見(jiàn)圖2、圖3.三細(xì)節(jié)孔板不做預(yù)腐蝕處理，直接進(jìn)行疲勞試驗(yàn)，共9件.五細(xì)節(jié)孔板先用3.5%NaCl溶液浸泡240 h，再進(jìn)行疲勞試驗(yàn)，共3件.

圖2 三細(xì)節(jié)孔板試件(單位:mm)

圖3 五細(xì)節(jié)孔板試件(單位:mm)

在疲勞試驗(yàn)過(guò)程中，分別記錄每一試件、每一孔的裂紋萌生壽命，得到第二條裂紋萌生時(shí)的疲勞壽命，即為多細(xì)節(jié)結(jié)構(gòu)發(fā)生MSD時(shí)對(duì)應(yīng)的疲勞壽命NMSD.

各個(gè)試件的存活率［4］為

破壞率(即MSD發(fā)生概率)為

其中，i表示各試件按出現(xiàn)多裂紋時(shí)對(duì)數(shù)壽命由小到大排列的順序號(hào);n表示該組試驗(yàn)的有效數(shù)據(jù)個(gè)數(shù).

三細(xì)節(jié)常規(guī)疲勞試驗(yàn)結(jié)果如表3所示，五細(xì)節(jié)預(yù)腐蝕后疲勞試驗(yàn)結(jié)果如表4所示.

表3 三細(xì)節(jié)常規(guī)試驗(yàn)結(jié)果

表4 預(yù)腐蝕后五細(xì)節(jié)試驗(yàn)結(jié)果

3.2 MSD發(fā)生概率

利用單細(xì)節(jié)常規(guī)疲勞試驗(yàn)結(jié)果μ=5.260 4，σ=0.1845，預(yù)腐蝕影響系數(shù)C=0.7506，及細(xì)節(jié)數(shù)n=3，n=5，通過(guò)文獻(xiàn)［3］中的方法及式(7)分別計(jì)算得到常規(guī)疲勞試驗(yàn)及預(yù)腐蝕影響后，三細(xì)節(jié)和五細(xì)節(jié)孔板的PMSD-lgN曲線.

三細(xì)節(jié)孔板預(yù)腐蝕及常規(guī)條件下的PMSD計(jì)算所得PMSD曲線及三細(xì)節(jié)孔板常規(guī)環(huán)境下疲勞試驗(yàn)結(jié)果如圖4所示.

圖4 三細(xì)節(jié)孔板計(jì)算及試驗(yàn)結(jié)果

從圖4可以看出，三細(xì)節(jié)孔板常規(guī)環(huán)境PMSD預(yù)測(cè)結(jié)果與試驗(yàn)結(jié)果相比，預(yù)測(cè)趨勢(shì)與試驗(yàn)結(jié)果一致.另外預(yù)腐蝕對(duì)MSD發(fā)生概率有促進(jìn)作用，預(yù)腐蝕后相同壽命對(duì)應(yīng)的PMSD明顯增大.這也證實(shí)了所建立的多部位損傷發(fā)生概率估算方法的合理性與可行性.

五細(xì)節(jié)孔板預(yù)腐蝕及常規(guī)條件下的PMSD計(jì)算所得PMSD曲線及五細(xì)節(jié)孔板預(yù)腐蝕后疲勞試驗(yàn)結(jié)果如圖5所示.

圖5 五細(xì)節(jié)孔板計(jì)算及試驗(yàn)結(jié)果

由圖5可以看出，五細(xì)節(jié)孔板預(yù)腐蝕后PMSD預(yù)測(cè)結(jié)果與試驗(yàn)結(jié)果相比，預(yù)腐蝕后同壽命對(duì)應(yīng)的PMSD值同樣明顯增大.另外在PMSD較小時(shí)，公式計(jì)算結(jié)果略偏于危險(xiǎn);在PMSD較大時(shí)，公式計(jì)算結(jié)果略偏于保守，并且計(jì)算結(jié)果與預(yù)腐蝕后疲勞試驗(yàn)結(jié)果一致.

4 預(yù)腐蝕對(duì)MSD發(fā)生概率的影響

結(jié)合預(yù)腐蝕處理后多細(xì)節(jié)孔板疲勞試驗(yàn)結(jié)果，采用固定相關(guān)變量，例如循環(huán)數(shù)N、相似細(xì)節(jié)數(shù)n、可靠度R'的方法，比較預(yù)腐蝕前后MSD發(fā)生概率或MSD發(fā)生時(shí)疲勞壽命的計(jì)算所得值，分析研究預(yù)腐蝕的影響.

4.1 預(yù)腐蝕后細(xì)節(jié)數(shù)n對(duì)PMSD的影響

由式(6)可知，式(7)中單一細(xì)節(jié)裂紋萌生概率pC(N)與細(xì)節(jié)數(shù)n無(wú)關(guān)，只與試驗(yàn)中的加載方式、預(yù)腐蝕方法以及對(duì)應(yīng)的疲勞壽命相關(guān).因此在應(yīng)力水平、預(yù)腐蝕條件確定的條件下，多細(xì)節(jié)孔板中任意一個(gè)細(xì)節(jié)達(dá)到壽命數(shù)N時(shí)萌生裂紋的概率pC(N)均為定值.因此可以把pC(N)看作常數(shù)，將PMSD對(duì)細(xì)節(jié)數(shù)n求導(dǎo)，可以得到

則在n＞1且0＜ pC(N)＜1時(shí)，PMSD對(duì)細(xì)節(jié)數(shù)n的導(dǎo)數(shù)大于0，即當(dāng)p(N)確定的情形下，多部位損傷發(fā)生概率PMSD是關(guān)于細(xì)節(jié)數(shù)n的增函數(shù).因此，在相同的試驗(yàn)條件、相同的疲勞壽命時(shí)，PMSD隨細(xì)節(jié)數(shù)n的增大而增大.

圖6 N=105時(shí)PMSD隨n的變化趨勢(shì)

在應(yīng)力水平取87.5 MPa時(shí)，給定載荷循環(huán)數(shù)N=105，取細(xì)節(jié)數(shù) n=3，5，7，9，…，計(jì)算 N=105時(shí)不同細(xì)節(jié)孔板的MSD發(fā)生概率，分別得到預(yù)腐蝕前后MSD發(fā)生概率隨細(xì)節(jié)數(shù)n的變化趨勢(shì)，如圖6所示.由圖6可看出，N=105且細(xì)節(jié)數(shù)較小時(shí)，相同細(xì)節(jié)數(shù)下，預(yù)腐蝕影響下的MSD發(fā)生概率明顯偏大，預(yù)腐蝕后的PMSD-n曲線斜率較大，PMSD增大較快.在細(xì)節(jié)數(shù)n較大的階段，預(yù)腐蝕與常規(guī)試驗(yàn)的PMSD均趨近于1，MSD必然發(fā)生.

4.2 n對(duì)預(yù)腐蝕后MSD發(fā)生時(shí)對(duì)數(shù)壽命的影響

給定應(yīng)力水平 S=87.5 MPa，可靠度R'=0.95，計(jì)算發(fā)生MSD時(shí)對(duì)數(shù)壽命lgNMSD隨細(xì)節(jié)數(shù)n的變化，如圖7所示.

圖7 R'=0.95時(shí)lgNMSD隨n的變化趨勢(shì)

由圖7可看出，對(duì)某一給定的可靠度，無(wú)論是否經(jīng)過(guò)預(yù)腐蝕處理，發(fā)生 MSD時(shí)的對(duì)數(shù)壽命lgNMSD均隨細(xì)節(jié)數(shù)n的增大而減小，且當(dāng)細(xì)節(jié)數(shù)較大時(shí)，變化趨于平緩，并接近于一條平行于x軸的漸近線.

5 結(jié)論

利用預(yù)腐蝕后疲勞試驗(yàn)測(cè)得的預(yù)腐蝕影響系數(shù)，計(jì)算了預(yù)腐蝕后的含相似多細(xì)節(jié)結(jié)構(gòu)中MSD發(fā)生概率，同時(shí)討論了預(yù)腐蝕對(duì)MSD發(fā)生概率的影響，可得到如下結(jié)論:

1)采用預(yù)腐蝕后PMSD公式計(jì)算得到的預(yù)腐蝕后的PMSD-lgN曲線，計(jì)算所得曲線可以較好地反映實(shí)際情況.預(yù)腐蝕與常規(guī)試驗(yàn)的PMSD曲線，當(dāng)PMSD較小時(shí)，計(jì)算結(jié)果均略偏于危險(xiǎn)，PMSD較大時(shí)，計(jì)算結(jié)果均略偏于保守.

2)隨細(xì)節(jié)數(shù)的增加，相同壽命時(shí)，預(yù)腐蝕與常規(guī)試驗(yàn)的PMSD均不斷增大，預(yù)腐蝕后的PMSD曲線斜率更大.

3)預(yù)腐蝕與常規(guī)試驗(yàn)相同可靠度下的NMSD均隨細(xì)節(jié)增大而不斷減小，預(yù)腐蝕與常規(guī)試驗(yàn)的lgNMSD-n曲線變化趨勢(shì)相同.預(yù)腐蝕后相同細(xì)節(jié)數(shù)的lgNMSD變小，預(yù)腐蝕明顯加速M(fèi)SD的發(fā)生.

References)

［1］ Draft FAA Advisory Circular No 120-YY Widespread fatigue damage on metallic structure as originally released in 2006 for public comment［S］

［2］ Federal Aviation Administration.Ageing airplane program:widespread fatigue damage，final rule［R］.FAA，14CFR Part 25，26，121，2010

［3］ Zhang Jianyu，Bao Rui，Zhang Xiang，et al.A probabilistic estimation method of multiple site damage occurrences for aircraft structures［J］.Procedia Engineering，2010，2(1):1115-1124

［4］高鎮(zhèn)同.疲勞應(yīng)用統(tǒng)計(jì)學(xué)［M］.北京:國(guó)防工業(yè)出版社，1986 Gao Zhentong.Fatigue applied statistics［M］.Beijing:National Defense Industry Press，1986(in Chinese)

［5］賀小帆，劉文珽，王忠波，等.疲勞壽命預(yù)腐蝕影響系數(shù)曲線通用性研究［J］.機(jī)械強(qiáng)度，2004，26(1):110-104 He Xiaofan，Liu Wenting，Wang Zhongbo，et al.Research on the generality of the pre-corrosion influence coefficient curve of fatigue life［J］.Journal of Mechanical Strength，2004，26(1):110-104(in Chinese)

［6］趙學(xué)峰，王富永，趙海軍.預(yù)腐蝕疲勞壽命影響系數(shù)及S-N曲線研究［J］.機(jī)械強(qiáng)度，2008，30(6):977-981 Zhao Xuefeng，Wang Fuyong，Zhao Haijun.Research on the influence coefficient and S-N curve for the pre-corrosion fatigue life［J］.Journal of Mechanical Strength，2008，30(6):977-981(in Chinese)