劉文光 ，陳國平

(1.南京航空航天大學 航空宇航學院，南京 210016;2.南昌航空大學 航空制造工程學院，南昌 330063)

含裂紋懸臂梁的振動與疲勞耦合分析

劉文光1，2，陳國平1

(1.南京航空航天大學 航空宇航學院，南京 210016;2.南昌航空大學 航空制造工程學院，南昌 330063)

基于Paris方程和同步分析方法考慮振動與疲勞裂紋擴展耦合之影響，提出一種含裂紋梁的振動疲勞壽命分析思路。振動分析過程中，利用線性彈簧等效裂紋段，復彈性模量引入阻尼損耗因子，得到考慮裂紋擴展、激勵頻率和阻尼等因素影響的動應力響應。結(jié)果表明:裂紋擴展、激勵頻率和阻尼等因素對疲勞壽命具有重要的影響。通過振動分析與疲勞裂紋擴展壽命估算同步進行，可進一步提高疲勞壽命估算精度。

振動疲勞;分段分析;裂紋擴展;裂紋懸臂梁

絕大多數(shù)工程結(jié)構(gòu)在使用過程中都不可避免要受到振動載荷的作用，諸如飛機、火箭、汽車、化工容器以及建筑橋梁，等等。振動不僅影響工程結(jié)構(gòu)的正常工作，還會引起振動疲勞破壞，嚴重地影響結(jié)構(gòu)的安全可靠性，直接威脅人類的生命安全，并造成重大的經(jīng)濟損失。因此，一般的工程結(jié)構(gòu)都把抗振動疲勞作為結(jié)構(gòu)的重要設(shè)計準則。

梁是工程中廣泛使用的一種結(jié)構(gòu)，由于各種原因?qū)е铝翰豢杀苊獾卮嬖诟鞣N損傷，大部分以裂紋的形式存在。長期處于振動環(huán)境中，結(jié)構(gòu)裂紋會不斷擴展并導致斷裂。為保證裂紋不發(fā)生疲勞失效，必須對工程結(jié)構(gòu)進行定期檢查。生產(chǎn)實際中，通常根據(jù)經(jīng)驗確定裂紋可檢測長度，而非靠預測方法來確定。

隨著斷裂力學的發(fā)展，提出了很多估算疲勞壽命的方法［1－3］。已知裂紋位置的名義應力或應力場，就能確定從起始裂紋擴展到某一臨界值所要經(jīng)歷的振動次數(shù)。然而，盡管裂紋擴展是因強迫振動引起，疲勞壽命分析中仍然采用靜力學方法計算應力，應力被內(nèi)在地假定與激勵頻率、阻尼無關(guān)，且常忽略裂紋擴展引起的動響應變化。工程實際中，大量結(jié)構(gòu)服役于寬頻激勵的振動環(huán)境中，使結(jié)構(gòu)可能一直處于共振狀態(tài)并導致共振疲勞。由于涉及振動響應，所以研究此類問題時必須考慮阻尼和激勵頻率對疲勞壽命的影響。

為此，有關(guān)激勵頻率對結(jié)構(gòu)振動疲勞壽命的影響已有不少研究［4～9］。例如，1981年 Valanis利用內(nèi)時理論推導了結(jié)構(gòu)疲勞壽命與加載頻率之間的關(guān)系;1993年毛罕平和陳翠英推導了全面反映加載頻率影響的裂紋擴展速率公式。隨后，Dentsoras和Kouvaritakis研究了軸向共振條件下激勵頻率對聚合物材料疲勞裂紋擴展壽命的影響;Shih和Chen分析了激勵頻率對304不銹鋼疲勞裂紋擴展速率的影響，并提出一種考慮激勵頻率的疲勞裂紋擴展模型。

但是，上述研究只考慮了激勵頻率的影響，忽略了阻尼以及裂紋擴展對疲勞壽命的作用。另外，由于裂紋擴展會改變結(jié)構(gòu)動態(tài)特性，所以結(jié)構(gòu)動響應分析時也須考慮裂紋擴展與振動的耦合作用。

因此，本文在結(jié)構(gòu)振動分析時，考慮裂紋擴展對振動的影響，得到考慮激勵頻率、結(jié)構(gòu)阻尼以及裂紋擴展的動態(tài)應力;基于Paris方程和同步分析方法使振動分析與疲勞壽命估算同步進行，進一步提高結(jié)構(gòu)疲勞裂紋擴展壽命估算精度。

1 振動分析

1.1 結(jié)構(gòu)裂紋的處理

圖1 裂紋梁振動系統(tǒng) (梁長l，寬b，高h，裂深a)Fig.1 Cracked Cantilever Beam(Length:l，Width:b;Height:h;Crack Depth:a)

圖1(a)所示為受橫向加載的單邊裂紋懸臂梁。在橫向簡諧振動載荷作用下，裂紋擴展會引起結(jié)構(gòu)單元剛度發(fā)生退化，并改變其動態(tài)特性［10］。為便于結(jié)構(gòu)動力學分析，需要對結(jié)構(gòu)裂紋作近似處理。通常采用線性扭轉(zhuǎn)彈簧與拉壓彈簧等效裂紋段，如圖1(b)所示。等效彈簧的剛度由結(jié)構(gòu)材料、裂紋類型以及裂紋位置等參數(shù)確定。根據(jù)斷裂力學的方法［13］，可推導結(jié)構(gòu)裂紋引起的局部柔度

式中:kT為扭轉(zhuǎn)剛度;kb為拉壓剛度;FⅠ2(x)為I型裂紋修正系數(shù);FⅡ(x)為Ⅱ型裂紋修正系數(shù);ζ=a/h為相對裂紋深度;E'=E/(1－v2)對應于平面應變情況。

根據(jù)裂紋類型以及加載形式［11］，裂紋修正系數(shù)為:

將上述裂紋修正系數(shù)分別代入公式(1)和(2)，可得到等效彈簧剛度的具體表達式。通過變剛度可以把裂紋擴展引進到結(jié)構(gòu)振動分析中。

1.2 結(jié)構(gòu)振動分析

結(jié)構(gòu)共振時，振幅及內(nèi)力值等驟增，此時必須考慮阻尼的影響。若采用復數(shù)阻尼理論，則等截面梁彎曲自由振動微分方程為:

如圖1，當裂紋在結(jié)構(gòu)根部時，結(jié)構(gòu)自由振動邊界條件修改為:

根據(jù)以上邊界條件得到考慮裂紋擴展變化的修正的固有頻率方程:

求解上述方程時，彈性模量只取其中的儲能模量部分。方程(6)表達了結(jié)構(gòu)固有頻率與裂紋深度之間的一一對應關(guān)系，對結(jié)構(gòu)裂紋的監(jiān)測具有十分重要的理論意義。

橫向力激勵下，結(jié)構(gòu)強迫振動邊界條件應修改為:

根據(jù)上述邊界條件可求得結(jié)構(gòu)位移響應:

由于Ⅰ型裂紋是工程中最常見的裂紋形式，且其破壞程度比較嚴重，故以Ⅰ型裂紋為研究重點。根據(jù)HOOKE定律得到結(jié)構(gòu)動應力響應:

結(jié)構(gòu)裂紋位置的動應力表達式為:

按靜力學方法，激勵頻率為0時應力為σ0=F0ly/I，則最大動應力可表示為:

式中:β為動應力放大系數(shù)，它反映了激勵頻率、結(jié)構(gòu)阻尼以及裂紋擴展對動應力的影響。

2 結(jié)構(gòu)振動疲勞裂紋擴展

2.1 動應力強度因子的確定

動應力強度因子KI是振動載荷作用下表征裂紋尖端應力場分布的物理量，通常而言，動應力強度因子與裂紋長度的關(guān)系可按下式確定:

式中:F(ζ)為裂紋修正因子;Δσd為動應力振幅;ΔKI是動應力強度因子的振幅。

假設(shè)加載為脈動循環(huán)，裂紋為開放式裂紋，則動應力強度因子可表示為:

把公式(12)代入公式(15)并依據(jù)裂紋及加載類型，可得到考慮裂紋擴展、激勵頻率和阻尼綜合作用的動應力強度因子

2.2 振動疲勞裂紋擴展速率方程

研究表明，在中等應力狀態(tài)且平均應力很低時，使用Paris方程［12］能很好地描述直裂式裂紋的疲勞擴展速率。Paris方程表示式如下:

式中:C，n是材料常數(shù);其余參數(shù)同上述公式。

把動應力強度因子表達式(16)代入Paris方程得到振動疲勞裂紋擴展速率模型:

2.3 振動疲勞裂紋擴展分析

相同結(jié)構(gòu)內(nèi)含不同大小的裂紋時，其動態(tài)特性必然不同，必定改變結(jié)構(gòu)的動應力響應并導致裂紋尖端應力場強發(fā)生變化。分析認為:結(jié)構(gòu)振動會引起裂紋擴展，而裂紋擴展也會改變結(jié)構(gòu)振動特性，進而影響裂紋擴展特性，兩者之間存在耦合關(guān)系。傳統(tǒng)的疲勞裂紋擴展理論卻忽略了這種耦合作用。為此，本文采用同步分析方法，使結(jié)構(gòu)振動分析與疲勞裂紋擴展壽命估算同步進行，考慮耦合作用對疲勞裂紋擴展壽命的影響，進一步提高疲勞壽命估算精度。

動應力作用下，不能直接采用Paris方程估算結(jié)構(gòu)疲勞壽命。為此，必須采用同步分析近似處理。具體做法是，把每振動一周得到的動應力幅值近似為恒定值，代入公式(20)積分得到每振動一次產(chǎn)生的裂紋擴展量。結(jié)構(gòu)受到變振幅載荷的作用時，最終的裂紋深度可由疊加法計算，表示如下:

式中，a0為初始裂紋深度;Δaj是第 j次循環(huán)的裂紋增量;i為總振動循環(huán)次數(shù);a為振動i次后的裂紋總深度。

利用公式(18)積分可計算恒幅值激勵下結(jié)構(gòu)疲勞裂紋擴展增量:

振動循環(huán)的增量可表示為:

考慮 ΔNj很小(本文取 ΔNj=1)，則:

因此，在裂紋中心點的裂紋擴展量可被方程式(23)近似計算出:

2.4 疲勞裂紋擴展失效判據(jù)

為了判斷結(jié)構(gòu)是否失效，文中采用以下判據(jù)

準則1:如果結(jié)構(gòu)裂紋擴展至梁的中面就認為結(jié)構(gòu)已經(jīng)破斷，此時裂紋的擴展不再適合用裂紋尖端的SIF來進行估算。

式中:ac為臨界裂紋長度，本文取ac=h/2。

準則2:如果裂紋尖端應力強度因子大于材料的斷裂韌性時，則認為結(jié)構(gòu)發(fā)生失穩(wěn)斷裂。

式中:Kmax為最大應力強度因子;KC為材料斷裂韌性。

3 結(jié)構(gòu)振動疲勞分析算例

假設(shè)結(jié)構(gòu)幾何尺寸 l=300 mm、h=20 mm、b=2 mm;材料為AISI1050低碳合金鋼［14］，材料參數(shù)為E=2.1E11Pa、ρ=7 860 kg/m3、σb=723.45E6、KC=1.172 2 E9 N/m3/2、ΔKth=9.3241E5 N/m3/2;試驗常數(shù)為 C=3.0 093E －32、m=3.3。初始裂紋深度 a0=2 mm;計算步長取ΔN=1周;激勵幅值為50 N，保持不變。

強迫振動分為共振和非共振兩種形式，影響其破壞的主要因素也有所變化。圖2～圖5分別表明了激勵頻率跟蹤結(jié)構(gòu)固有頻率、共振區(qū)固定激勵頻率、非共振區(qū)固定激勵頻率三種激勵狀態(tài)下結(jié)構(gòu)振動循環(huán)次數(shù)與裂紋擴展之間的內(nèi)在規(guī)律。結(jié)果表明，激勵頻率和阻尼對結(jié)構(gòu)疲勞壽命具有重要的影響，在共振環(huán)境下尤其明顯。相同幅值的外部激勵作用下，第一階模態(tài)共振引起的疲勞破壞更為嚴重，小阻尼對疲勞裂紋擴展速率影響不是很明顯。激勵頻率跟蹤結(jié)構(gòu)固有頻率時的結(jié)構(gòu)疲勞壽命遠遠小于其它強迫振動下的結(jié)構(gòu)疲勞壽命。固定激勵頻率時，激勵頻率越接近結(jié)構(gòu)固有頻率，結(jié)構(gòu)疲勞裂紋擴展壽命越短。圖4、圖5中α表示激勵頻率與無裂紋結(jié)構(gòu)第一階固有頻率之比。

圖2 激勵頻率跟蹤固有頻率的疲勞裂紋擴展曲線(不同模態(tài))Fig.2 Crack Growth vs.Fatigue Life with Resonance Frequency Tracking(Different Mode)

4 結(jié)束語

為便于分析，用線性彈簧近似取代裂紋單元、復彈性模量引入阻尼;通過結(jié)構(gòu)振動分析，得到考慮裂紋擴展、阻尼以及激勵頻率等影響的動應力，與常規(guī)疲勞裂紋擴展分析相比，應力的計算更符合實際;基于疲勞裂紋擴展方程和同步分析方法，考慮振動與裂紋擴展的耦合關(guān)系對疲勞壽命的影響，使振動分析與疲勞裂紋擴展壽命估算同步進行，提出了一種含裂紋結(jié)構(gòu)的振動疲勞分析思路，進一步提高結(jié)構(gòu)疲勞壽命估算精度。分析結(jié)果表明，裂紋擴展、結(jié)構(gòu)阻尼以及激勵頻率對疲勞壽命具有重要的影響，在共振區(qū)尤其明顯;跟蹤第一階共振頻率的結(jié)構(gòu)疲勞裂紋擴展過程十分迅速，小阻尼對疲勞擴展影響不顯著;跟蹤固有頻率的結(jié)構(gòu)疲勞裂紋擴展壽命遠小于其它強迫振動作用下的結(jié)構(gòu)疲勞裂紋擴展壽命;遠離共振環(huán)境后，阻尼對結(jié)構(gòu)疲勞裂紋擴展壽命的影響可以忽略不計，但疲勞裂紋擴展速率隨激勵頻率離固有頻率越近疲勞壽命越短。

［1］Ponomarev P V，Lopatin A D.Calculation of the fatigue fracture underthe influence ofdynamic loads［J］.International Applied Mechanics，1972，8(6):613 －617.

［2］Shih Y S，Wu G Y.Effect of vibration on fatigue crack growth of an edge crack for a rectangular plate ［J］.International Journal of Fatigue，2002，24:557 566.

［3］Schlums D H.Fatigue testing and crack analysis of resonating structures［D］.Swiss Federal Institute of Technology Zurich，2001.

［4］Yeh J T，Hong H S.Effect of branch frequency on dynamic fatigue behavior of slowly notched polyethylene polymers［J］.Journal of Polymer Research，1994，4(1):375－383.

［5］ Dentsoras A J，Kouvaritakis E P. Effects of vibration frequency on fatigue crack propagation of a polymer at resonance ［J］.Engineering Fracture Mechanics，1995，50(4):467－473.

［6］ Valanis K C.On the effect of frequency on fatigue life［C］.Presented at the winter annual meeting of the American Society ofMechanicalEngineering， Washington D.C，November，1981:15 －20.

［7］ Shih Y S，Chen J J.The frequency effect on the fatigue crack growth rate of 304 stainless steel［J］.Nuclear Engineering and Design，1999，191:225 230.

［8］毛罕平，陳翠英.考慮動應力時加載頻率對裂紋擴展速率的影響［J］.江蘇工學院學報，1993，14(6):26－33.

［9］沈 為，高大興，蔣松山，等.加載頻率對疲勞裂紋擴展速率的影響［M］.武漢:湖北人民出版社，1980，2:186－195.

［10］ Papaeconomou N，Dimarogonas A.Vibration of cracked beams［J］.Computational Mechanics，1989，5:88 －94.

［11］Tada H，Paris P C，Irwin G.R.The stress analysis of cracks handbook ［M］.DelResearch Corporation Hellertown，Pennsylvania.1973.

［12］ Paris P C， Erdogan F. A criticalanalysisofcrack Propagation laws［J］.Journal of Basic Engineering，1963，85:528－534.

［13］ Ewalds H L，Wanhill R J.Fracture mechanics［M］.Edward Arnold Ltd，1984.

［14］ Dentsoras A J，Dimarogonas A D.Resonance controlled fatigue crack propagation in a beam under longitudinal vibration ［J］.International Journal of Fracture，1983，23:15－22.

Coupling analysis for vibration and fatigue of a cracked cantilever beam

LIU Wen-guang1，2，CHEN Guo-ping1

(1.College of Aerospace Engineering，Nanjing University of Aeronautics and Astronautics，Nanjing 20016，China;2.School of Aeronautical Manufacturing Engineering，Nanchang Hangkong University，Nanchang 330063，China)

An analytical method about vibration and fatigue of a cracked cantilever beam was proposed based on timing analysis and Paris equation.Coupling relationship between vibration and fatigue crack growth was considered by using timing analysis methods.Crack part was replaced with a linear spring and damping loss factor was introduced with complex elastic modulus during vibration analysis.So crack growth，excitation frequency and damping were considered in vibration analysis.Results indicated that crack length，excitation frequency and damping have important effects on the fatigue life of the beam;in the resonance region，the rate of crack growth is obvious and the speed decreases with increase in damping.The prediction precision of fatigue life could be further improved using this new analytical method.

vibration and fatigue;piecewise analysis;crack propagation;cracked cantilever beam

V216.3;V224

A

南昌航空大學?；鹳Y助(EC200906236)

2009－12－02 修改稿收到日期:2010－03－22

劉文光 男，博士生，講師，1978年5月生