含裂紋缺陷結構J積分計算方法研究

含裂紋缺陷結構J積分計算方法研究*

黃裕龍，蔣 瑋，蔣險峰

(大連理工大學機械工程學院，遼寧大連116024)

摘要：裂紋參量J積分一直是斷裂力學領域研究的熱點，它是定量描述裂紋尖端區域應力應變場強度的重要斷裂參量。J積分的獲取有多種方法，歸納總結了幾種常用的含裂紋缺陷結構的J積分計算方法并對其進行分析，最后列舉一實例進行計算，分析了各類方法的特點，為斷裂參量J積分的研究以及計算提供參考。

關鍵詞:斷裂J積分EPRICOD有限元

中圖分類號：O346.1文獻標識碼：A

基金項目：國家自然科學基金資助項目51075050。

作者簡介：黃裕龍(1988-)，男，碩士研究生，主要研究方向為鋼結構裂紋特性。

收稿日期：2015-03-04

Research of cracked defect structure’s J-integral computing methods

HUANG Yulong,JIANG Wei,JIANG Xianfeng

Abstract:Crack parameter J-integral has been a hot issue of fracture mechanics. It is an important fracture parameter which quantitatively describes the stress and strain fields near crack tip. There are varieties of methods to obtain J-integral. This paper summarizes several methods of calculating the J-integral of cracked defect structure and analyzes the result. The paper cites an example at the end, and analyzes the characteristics of various methods, provides a reference for the future study and computing of J-integral.

Keywords:fracture; J-integral; EPRI; COD; finite element method

0引言

J積分是研究裂紋問題的一個重要參量，用以表征裂紋尖端附近的應力應變關系，可用來研究裂紋性質并能判斷裂紋是否失穩擴展，對預防和控制斷裂行為具有重要意義。自從1968年Rice[1]提出J積分以來，它就成為了一個重要的斷裂力學參數，一度成為國內外學者的研究重點。J積分代表作用于裂紋尖端的一個廣義力，也稱為裂紋擴展力或能量釋放率。研究J積分的計算和測量方法，獲取J積分值，將有助于判斷裂紋的擴展和含裂紋結構的可靠性。

目前，對于J積分的計算一般有兩大類：一類是按照簡化模型或者工程估算方法進行計算，另一類是通過有限元分析或者實驗方法獲取J積分值。本文將從這兩類方法入手，對J積分的計算方法進行研究，并在最后給出具體實例進行對比研究，比較各類計算方法的特點。

1J積分的計算方法

設有一均質板，板上有一穿透裂紋，裂紋表面為自由表面，外力使裂紋周圍產生二維的應力、應變場。定義J積分[1]如下：

(1)

圖1　J積分路徑

式中，Γ為始于裂紋下表面任一點，止于裂紋上表面的圍繞裂紋尖端的閉合路徑，如圖1；W為應變能密度：

W=∫σijdεij

(2)

其中σij、εij為Γ路徑上任一點的應力應變值；T是積分路徑Γ邊界上的應力矢量，U是路徑Γ上的位移矢量，dy為y方向的增量，ds為對應于Γ的積分曲線增量。

在線彈性斷裂力學中，J積分與應力強度因子KI有如下關系式：

(3)

1.1EPRI工程估算法

為了達到在工程實際中不借助有限元即可完成J積分計算的目的，美國通用電器(GE)和美國電力研究院(EPRI)在其研究報告NP-1931《彈-塑性斷裂分析的工程方法》中提出了EPRI的J積分工程估算法及其延性斷裂手冊[2-4]。運用該方法，在面對工程實際計算J積分的問題時，無需進行有限元分析就能將J積分計算出來，簡化了J積分的計算，推動了J積分在工程實際中的應用[5]。EPRI的J積分工程估算法的關鍵是將彈-塑性J積分分解成彈性J積分與塑性J積分兩部分相加[6-8]。

對于在單項拉伸的情況下，該方法認為材料為Ramberg - Osgood冪硬化材料，其應力、應變關系符合以下關系式：

(4)

其中，σ、ε材料在某一時刻的應力和應變；σ0、ε0為材料的屈服應力和屈服應變；α為材料的硬化系數；n為材料的硬化指數。

對其進行彈塑性分析，根據EPRI的J積分計算方法，可以將J積分換算為彈性階段的J積分Je和塑性階段的J積分Jp兩者之和，即：

J=Je(ae)+Jp(α，n)

(5)

其中，彈性階段的J積分Je是關于修正以后的有效裂紋長度ae的函數：

(6)

式中，m為形狀系數，它與試件裂紋大小、位置等因素有關；σ為名義應力，即裂紋位置上按照無裂紋計算的應力；ae可以由下式方程求得：

ae=a+φγy

(7)

(8)

對于平面應力，β=2；對于平面應變β=6。

(9)

(10)

式中，Lr為外載荷P與結構屈服載荷P0(即塑性極限載荷)之比。

對于平面應力，P0=1.071 ηcσ0；對于平面應變P0=1.455ηcσ0。

(11)

c為試件的未開裂韌帶。

對于塑性階段的J積分Jp，可由下式計算而得：

(12)

h1為J積分全塑性解的函數，與材料、結構、載荷以及裂紋尺寸有關。通過有限元分析，可以得到h1值，編制成全塑性解系數h1手冊。當需要計算J-積分時，只要查表確定h1的值就可完成計算工作，不再需要有限元計算，因此在工程實際應用中利用全塑性解h1是非常簡單的[9]。

1.2COD準則法

J積分是一個回路積分，COD是裂紋尖端張開位移，它們是彈塑性斷裂力學的兩個重要參量，以這兩個參量為基礎建立了彈塑性斷裂力學的重要斷裂準則[10]。在D-M模型的基礎上，J積分和COD有一定關系，如下：

J=σsδ

(13)

式中的σs為屈服強度。

由于D-M模型過于簡單，將塑性區考慮為理想塑性，塑性區受到σs的作用；實際上許多材料都存在硬化現象，不能簡單的認為J積分和COD的關系僅僅如(13)式。人們做了許多的實驗，發現應該在(13)式的基礎上加一個系數k，稱之為COD的降低系數，其值大約在1～3左右。新關系式如下：

J=kσsδ

(14)

Robinson[11]指出，k隨著塑性區的增加而增加，在塑性區較小時，k=1，k最大可以達到2.6。Shih[12]的計算結果表明，用HRR理論和有限元計算結果相差不大，k隨著硬化指數的增加而減小。Xian-KuiZhu[13]等人通過認真總結與歸納，給出了不同標準試件的COD降低系數k的算法，其中就有CT試件的k的算法，如下：

(15)

其中，σts為材料的極限抗拉強度。

對于三點彎曲式樣，有：

對于緊湊拉伸式樣，有

A0=3.62A1=4.21A2=4.33A3=2.00

1.3有限元法

有限元法(Finite Element Method, FEM)是固體力學領域中應用最廣泛的數值模擬方法，它的基本思想是用較簡單的問題代替復雜問題后再求解。應用有限元方法求斷裂力學參數J積分早在80年代Owen就已經開始了[14]，在他的著作中先后就彈性與塑性材料在二維情況下作了較詳盡的論述并編制了計算機程序。目前的有限元軟件如ANSYS、ABAQUS等已經得到了很好的發展，對于各種問題均能夠很好的解決。林國裕[15]利用有限元方法對具有中心裂紋的矩形板作了線彈性和彈塑性分析，得到了積分基本守恒和隨時間的變化規律。李繼紅[16]等人利用彈塑性損傷-應變耦合的有限元方法,研究了在不同強度匹配焊接接頭中，母材區斷裂應變的變化對焊縫區裂紋斷裂力學參量J積分的影響規律。

運用有限元法計算J積分，需要建立對應的模型，下文實例中將用到該方法計算，設置材料屬性為：彈性模量E=195.1 GPa，泊松比μ=0.267，密度ρ=7.93×103kg/m3。劃分網格并施加約束及載荷，提交運算即可算出J積分值。

1.4實驗法

用實驗法求解J積分，多數情況下是用以對理論方法的驗證。實驗法的中心思想是想方設法獲取裂紋尖端附近區域的應力應變值，通常是先獲取應變值，再通過彈塑性斷裂力學應力應變關系求解應力值，進而計算J積分。傳統的實驗方法有應變片法[17]、柔度法[18]、光彈法[19]。傳統的實驗法存在著一定的不足：貼片法繁瑣且該方法測量點有限，柔度法無法用于動態測量，光彈法的實驗設備復雜，而且無法用于非透明材料的測量[20]。因此，迫切需要一種先進的技術來克服傳統實驗方法存在的缺陷，人們對于尋求更有效實驗方法的努力從未停止過。

隨著科技的進步和發展，近些年出現了運用數字圖像相關技術的方法來獲取裂紋J積分值的方法，并一度成為研究的熱點。數字圖像相關由M.A. Sutton等人提出[21-22]，通過跟蹤匹配變形前后所采集圖像的灰度信息來測量物體在各種載荷作用下表面整體的瞬時位移場和應變場。該方法具有非接觸、精度高、光路簡單、受環境影響小、自動化程度高等優點，已經逐漸成為實驗力學領域一種非常重要的光測手段[23]。運用數字圖像相關技術，能夠輕松的獲取裂紋尖端附近區域的實時應變值，再通過應力應變關系求出應力值，即可獲取J積分值。

2實例計算及分析

圖2　CT式樣尺寸

以緊湊拉伸(CT)式樣為例進行說明，材料選用304不銹鋼，厚度為5 mm，對其進行單項拉伸，CT式樣尺寸如圖2所示。為了效驗計算結果的準確性，將各方法的計算結果與(3)式用應力強度因子法計算的理論值做比較。根據上述計算J積分的方法，對該實例進行計算，結果如圖3所示。

由圖3可知，各種J積分計算方法的結果具有相同的增長趨勢且相差不大，通過與理論值的對比，驗證了各種常用J積分計算方法的準確性與可行性。

圖3　J積分計算結果圖

3結語

本文歸納總結了斷裂力學領域中斷裂參量J積分計算的若干常用方法，并通過具體實例計算分析，驗證了各種方法的可行性和準確性。估算法雖是對工程實際問題的簡化計算，但能夠對裂紋構件斷裂的全過程給予定量描述，在工程實際當中運用該方法尤為方便且準確可行；COD準則法簡單有效，但對于COD值的精度要求比較高；有限元方法和實驗法均通過獲取應力應變值計算J積分，不同點是有限元方法通過軟件的內部算法直接獲取J積分值而無需單獨獲取應力應變，實驗法則需要先獲取應力應變才能獲得J積分值。

參考文獻

[1]RICE J.R. Plane strain deformation dear a crack tip in a power-law hardening material [J]. Journal of The Mechanics of Solids.1968,16:1

[2]Kumar V, German M D, Shih C F. Translated by ZHOU Fanhong, et al. An engineering approach for elastic-plastic fracture analysis . Beijing :National Defense Publishing House , 1985 ( In Chinese)

[3]Kumar V, German M D. Wilkening W W. Development on elastic-plastic fracture analysis (EPRI NP - 3607). Palo Alto, California: Research & Development Center of EPRI, 1984

[4]EPRI and Novetech Corporation. Ductile fracture handbook. Vol. 1~3, Palo Alto, California: Research Report Center,1991-1993A

[5]孫亮，李培寧，劉長軍. 任意應力應變關系材料J積分估算的等效原場應力法[J]. 機械強度，2003，25(4)：450-455

[6]ZAHOOR, A. Fracture of circumferentially cracked pipes[J]. ASME J. of Pressure Vessel Technology, 1987, vol. 108: 529-531

[7]ZAHOOR, A. Evaluation of J-Integral estimation scheme for flawed through-wall pipes[J]. Nuclear Engineering and Design, 1987, vol. 100: 1-9

[8]ZAHOOR, A. J-integral analysis for notched round bar in tension[J]. ASNE J. of Pressure Vessel Technology, 1987, vol. 109: 155-158

[9]杜國強. 斜接彎管內側裂紋斷裂參數研究[D]. 大連理工大學碩士學位論文，2009

[10]張莉. 影響淺裂紋試件J-COD關系的兩種因素[J]. 哈爾濱理工大學學報. 2000，5(6): 7-9

[11]Robinson J N. Int. J. of Fract. 1976, 12:5

[12]Shih C F. J. Mech. Phys. Solids, 1981, 29: 305

[13]Xian-Kui Zhu, James A. Joyce. Review of fracture toughness (G, K, J, CTOD, CTOA) testing and standardization [J]. Engineering Fracture Mechanics, 2012, vol85: 1-46

[14]Owen,D.R.J..Engineering fracture mechanics: Numerical methods and applications[M]. Pineridge Press Ltd. Swansea, U.K., 1983

[15]林國裕，董友成. 彈塑性動態J積分的有限元分析[J].華僑大學學報，1989,10(2)：156-161

[16]李繼紅，張建勛. 含損傷焊接接頭裂紋體J積分的有限元數值研究[J].2003,37(9)：937-940

[17]Dallly J W, Sanford R J. Strain gage methods for measuring the opening-mode stress intensity factor [J]. Experimental Mechanics, 1987, 27: 381-388

[18]Bonesteel R M, Pipers D E, Davinroy A T. Compliance and calibration of double cantilever beam (DCB) specimens [J]. Engineering Fracture Mechanics, 1978, 10: 425-428

[19]Wang W C, Hwang C H. Experimental analysis of vibration characteristics of an edge-cracked composite plate by ESPI method [J]. International Journal of Fracture, 1998, 91(4): 311-321

[20]張蕊，賀玲鳳. 數字圖像相關法測量聚碳酸酯板應力強度因子[J]. 工程力學，2012，29(12)：22-27

[21]Peters W H, Ranson W F . Digital imaging techniques in experimental stress analysis [J]. Optics and Precision Engineering. 1982，21(3):427-431

[22]Chu T C, Ranson W F, Sutton M A, et. Applications of digital image correlation techniques to experimental mechanics[J]. Experimental Mechanics , 1985,25(3):232-245

[23]胡浩，梁晉，唐正宗，盧崗. 數字圖像相關法測量金屬薄板焊接的全場變形[J]. 光學精密工程,2012，20(7)，1 636-1 643