夏世林，磨季云

(武漢科技大學(xué) 理學(xué)院，湖北 武漢 430065)

0 引 言

加筋板結(jié)構(gòu)是船舶、航空航天領(lǐng)域最為常用的結(jié)構(gòu)類型之一。對于此類結(jié)構(gòu)，當(dāng)結(jié)構(gòu)受力達(dá)到材料屈服強(qiáng)度后，其仍有一定的承載能力，直至載荷水平達(dá)到其極限強(qiáng)度。在加筋板結(jié)構(gòu)設(shè)計中，極限狀態(tài)（ULS）設(shè)計已經(jīng)被公認(rèn)為優(yōu)于傳統(tǒng)的許用應(yīng)力方法，因為前者能更真實(shí)反映結(jié)構(gòu)的安全極限[1]。加筋板結(jié)構(gòu)在服役年限中不可避免遭受不同種類的損傷，裂紋損傷就是影響結(jié)構(gòu)失效行為的重要因素之一。在加筋板結(jié)構(gòu)中，裂紋經(jīng)常在焊接和應(yīng)力集中區(qū)域萌生，例如在加強(qiáng)筋和底板的交接處。在加筋板壽命周期中，這些裂紋可能在不同載荷狀態(tài)下擴(kuò)展，結(jié)果形成不同的尺寸和位置。

按載荷形式的不同，含裂紋加筋板的剩余強(qiáng)度問題可分為壓縮載荷下加筋板的屈曲強(qiáng)度、拉伸載荷下加筋板的極限強(qiáng)度和剪切載荷下加筋板的剩余強(qiáng)度。對于此類問題，傳統(tǒng)的有限元方法是將裂紋建立成類似于狹縫的幾何模型，并在裂紋尖端進(jìn)行網(wǎng)格加密，從而將斷裂力學(xué)問題簡化為一般靜力學(xué)問題。Sadek等[2]用數(shù)值模擬的方法分析了含裂紋加筋板在壓縮應(yīng)力下的屈曲特性，分析中考慮了裂紋的演化發(fā)展；Cao等[3]根據(jù)歐洲規(guī)范EN1990，EN1991和EN1993，對環(huán)向和經(jīng)向制造裂紋的圓形鋼筒倉的屈曲性能進(jìn)行了評估；Duc等[4]通過數(shù)值模擬計算，建立了基于裂紋有限元分析的靜穩(wěn)定板模型。Brighenti等[5]研究了裂紋損傷對拉伸、壓縮、剪切載荷下的矩形板屈曲強(qiáng)度的影響，其中裂紋長度、裂紋方向和材料泊松比是變化的，結(jié)果表明裂紋方向垂直拉伸和壓縮載荷方向時是最危險的工況。Margaritis Y等[9]研究了含裂紋加筋板的極限和失效響應(yīng)，詳細(xì)檢查了裂紋閉合對結(jié)構(gòu)響應(yīng)的作用。Paik等[10]對完整和含裂紋加筋板的極限剪切強(qiáng)度進(jìn)行了數(shù)值研究，并基于完整加筋板方程發(fā)展了預(yù)報含裂紋加筋板的極限剪切強(qiáng)度方程。Saad-Eldeen等[11]對含中心橢圓開口和封閉裂紋的薄壁鋼板剩余強(qiáng)度進(jìn)行了試驗研究，試驗中的試件是從服役期中的實(shí)際結(jié)構(gòu)構(gòu)件中取得的，而不是新建試件。

與傳統(tǒng)的一般靜力學(xué)分析方法相比，J積分方法具有以下優(yōu)勢：一是不需要建立裂紋模型，而是從斷裂力學(xué)機(jī)理出發(fā)，直接定義裂尖和裂紋方向，這更加符合裂紋模型的實(shí)際情況；二是J積分方法避開了直接求解裂端塑性應(yīng)力場的困難。Cherepanov[6]和Rice[7]各自獨(dú)立提出了一種與積分路徑無關(guān)的J積分，用于綜合度量裂端應(yīng)力應(yīng)變場強(qiáng)度。它對彈塑性斷裂力學(xué)發(fā)展起到了重要的作用。

通過以上文獻(xiàn)可以看出，以前的研究主要關(guān)注含裂紋結(jié)構(gòu)在典型拉伸載荷下的失效行為或塑性屈曲強(qiáng)度評估。含裂紋結(jié)構(gòu)在軸向拉伸下的剩余極限強(qiáng)度在近年來得到了廣泛關(guān)注，但是該領(lǐng)域的研究仍然是有限、不充分的。因此，需要進(jìn)行更多的研究來更深理解含裂紋缺陷結(jié)構(gòu)構(gòu)件的失效行為。本文的主旨在于發(fā)現(xiàn)含裂紋損傷的加筋板在單調(diào)軸向拉伸下的極限強(qiáng)度行為。假定板橫向裂紋位于加筋板內(nèi)截面，或只位于板單元上，或位于板單元和加強(qiáng)筋上。需要指出的是，本文不考慮裂紋擴(kuò)展和斷裂相關(guān)的臨界裂紋長度。

加筋板材料為Q345B高強(qiáng)度鋼，為了達(dá)到目的，進(jìn)行了含一系列不同尺寸、位置裂紋加筋板在單軸拉伸載荷下的非線性彈塑性大變形分析。

1 力學(xué)分析方法

1.1 J積分及其有限元列式

J積分理論避開了直接計算裂紋尖端附近的彈塑性應(yīng)力、應(yīng)變場，而用遠(yuǎn)場J積分作為表示裂紋尖端應(yīng)力應(yīng)變集中特征的平均參量。對于二維問題，J積分可定義為[3]：

其中： Γ 為由裂紋下表面某點(diǎn)到裂紋上表面某點(diǎn)的簡單的積分路徑；W為彈塑性應(yīng)變能密度；Ti為作用于積分回路單位周長上的主應(yīng)力；ui為積分回路邊界上的位移；ds為積分回路線的弧長。

在線彈性斷裂力學(xué)中，J積分與應(yīng)力強(qiáng)度因子KI有如下關(guān)系式：

式中：對于平面應(yīng)力，彈性模量E′=E，對于平面應(yīng)變E′=E/1-μ2。

可以解析證明，J積分與積分路線 Γ 的選取無關(guān)，即J積分滿足路徑守恒。因此，可選取應(yīng)力應(yīng)變場較易求解的線路得到J積分值，而此值與積分路徑非?？拷讯说慕Y(jié)果是相同的。J積分有著明確的物理意義，它代表著裂紋沿x方向擴(kuò)展單位長度時的能量釋放率。

由于J積分與選擇的路徑無關(guān)，并能描述裂紋尖端的應(yīng)變狀態(tài)。這一特點(diǎn)使J積分計算可以在離開高應(yīng)力、高應(yīng)變梯度的裂紋尖端區(qū)外的任意回線進(jìn)行。對于二維體的J積分列式可以寫為[11]：

其中：

將平面上定義的二維J積分沿著裂紋縫線逐點(diǎn)積分可得裂紋線的三維J積分為：

式中：L為裂紋縫線長度；dA為包繞裂紋縫線的圓柱曲面上的面積微元， dA=dsdΓ；n為微元dA的外法線方向； λ (s)為裂紋縫線s點(diǎn)處的長度。

J積分的計算必須采用數(shù)值積分，如果回線經(jīng)過單元積分點(diǎn)，可以采用有限元中的高斯積分策略。從遠(yuǎn)離裂紋尖端的應(yīng)力場、應(yīng)變場計算J積分，只能在無體力、初始應(yīng)變和裂紋表面拉力情況下能用。當(dāng)有體力或初始應(yīng)變時，J積分要加入封閉表面的積分。此時路徑無關(guān)性將喪失，必須從完整的應(yīng)力和應(yīng)變結(jié)果中計算，包括裂紋尖端的高應(yīng)力區(qū)。

J積分的主要優(yōu)點(diǎn)是不需要用很細(xì)的網(wǎng)格對裂紋尖端區(qū)進(jìn)行分析。路徑無關(guān)性的非線彈性材料也存在，真實(shí)的彈塑性材料與非線彈性材料類似。J積分也可以用于描述彈塑性材料的裂紋尖端的應(yīng)力場和應(yīng)變場，材料用增量塑性理論描述并不允許有卸載，以免W的加載歷程相關(guān)性，J值的增量可以從每個載荷步的應(yīng)力和應(yīng)變增量計算得到。

1.2 精度驗證算例

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

圖1 CT試樣尺寸Fig.1 Dimension of the CT specimen

圖2 J積分計算結(jié)果圖Fig.2 The scheme of J-integral calculation result

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

1.3 有限元模型

1.3.1 幾何和材料

本文的分析對象為含2根加強(qiáng)筋（截面對稱）和1個板單元的加筋板，其長細(xì)比范圍為4～6，如圖3所示。在圖3中，板的尺寸用a×b×t表述，加強(qiáng)筋的尺寸用hw×tw表述。底板和加強(qiáng)筋的詳細(xì)尺寸見表1。

圖3 加筋板的模型尺寸Fig.3 Model dimension of the stiffened-panel

表1 加筋板的幾何特征Tab.1 Geometrical features of the stiffened-panel

其中長細(xì)比β是控制加筋板屈曲/失效行為的重要參數(shù)，定義如下[13]：

其中：b為板寬；t為板厚； σs為材料屈服強(qiáng)度；E為楊氏模量。

本文使用的材料為Q345B鋼。一般來說，鋼結(jié)構(gòu)具有應(yīng)變硬化正切模量Et，其典型取值范圍為楊氏模量的5%～15%。顯然應(yīng)變硬化對加筋板的的非線性行為有一定影響。影響的程度是包括板長細(xì)比在內(nèi)的眾多因素的函數(shù)。在本文中，板的材料行為建模成應(yīng)變硬化率為E/65的雙線性彈塑性行為（見圖4），該應(yīng)變硬化率的值是Khedmati[14]通過大量彈塑性大變形分析得到的。詳細(xì)材料屬性如表2所示，材料的應(yīng)力-應(yīng)變行為如圖4所示。

1.3.2 裂紋損傷研究工況

為了考慮裂紋的不同位置，本文應(yīng)用了3種不同的模型，如圖5所示。這些模型按裂紋位置分類，如裂紋只在底板上，裂紋在底板和加強(qiáng)筋上，裂紋只在加強(qiáng)筋上。

裂紋在加筋板中的深度取決于多種因素，其中裂紋深度等于板厚時最為危險。在工程中，在相對較薄的加筋板中，裂紋一般被認(rèn)為貫穿板厚。因此，在本文研究中，底板和加強(qiáng)筋上的裂紋都假設(shè)貫穿厚度。在所有的3種模型中，不論裂紋在底板上還是加強(qiáng)筋上，都沿板的橫向。為了考慮不同裂紋長度，假設(shè)加筋板上的裂紋長度c是對應(yīng)板寬的某一特定比例，其比值c/b的變化范圍為0.05～0.15，而圖5（b）和圖5（c）中在加強(qiáng)筋上的裂紋長度等于加強(qiáng)筋高度。表3描述了本文考慮的不同幾何參數(shù)。

表2 Q345B鋼材料屬性匯總Tab.2 A summary of Q345B steel material properties

圖4 材料的理想彈塑性應(yīng)力-應(yīng)變行為Fig.4 Ideal elastic-plastic stress-strain behavior of the material

圖5 3種裂紋模型Fig.5 3 types of crack models

表3 FE模型幾何參數(shù)值匯總Tab.3 A summary of FE models’ geometrical parameters

1.3.3 載荷和邊界條件

在進(jìn)行FE分析之前，在FE模型上施加合適的邊界條件盡可能合理模擬結(jié)構(gòu)實(shí)際響應(yīng)是十分重要的。一個常規(guī)的方法是在極限狀態(tài)評估中在橫向邊緣施加簡支或夾持邊界條件。本文主要關(guān)注的是加筋板在單軸拉伸載荷下的極限強(qiáng)度。通常認(rèn)為加筋板橫剖面會受到均勻的拉伸外載荷力，可采用位移控制的方式來實(shí)現(xiàn)加筋板均勻軸向載荷的逐步加載。圖6描述了含裂紋加筋板模型的載荷和邊界條件，其中t為時間。

圖6 含裂紋加筋板載荷和邊界條件Fig.6 Loads and BCs of the stiffened-panel with cracks

1.3.4 單元和網(wǎng)格尺寸

結(jié)構(gòu)極限強(qiáng)度計算是非線性問題。在本文的研究中，含裂紋加筋板的失效行為和極限強(qiáng)度評估使用有限元軟件Abaqus，評估中同時考慮幾何和材料非線性。在Abaqus包含可使用的單元庫中，選擇8節(jié)點(diǎn)S8R單元來離散加筋板模型。S8R適合屬于加筋板結(jié)構(gòu)范疇的適度薄壁殼結(jié)構(gòu)。其適合線性、大應(yīng)變和大轉(zhuǎn)動非線性應(yīng)用。

在FE分析中，計算精度和速度在很大程度上取決于單元大小或網(wǎng)格密度。網(wǎng)格尺寸必須足夠精細(xì)來捕捉含裂紋加筋板的屈曲/塑性失效行為。另一方面，更為精確的網(wǎng)格引起更大的計算代價。從這一點(diǎn)來說，必須進(jìn)行精度和計算時間的平衡。對于含裂紋加筋板，精確網(wǎng)格必須分配在局部裂紋區(qū)域來模擬更早的屈服和由于應(yīng)力集中引起的高應(yīng)力梯度。

進(jìn)行不同網(wǎng)格尺寸的收斂性研究決定相對合適的非線性有限元分析網(wǎng)格模型。最終，沿板寬方向b和板長方向a的殼單元數(shù)目指定為20。對于加強(qiáng)筋腹板在高度方向hw單元數(shù)量為1。同時在裂紋附近進(jìn)行網(wǎng)格細(xì)化。裂紋區(qū)域網(wǎng)格細(xì)化的含裂紋加筋板網(wǎng)格劃分實(shí)例如圖7所示。

圖7 含裂紋加筋板網(wǎng)格實(shí)例Fig.7 Mesh example of the stiffened-panel with cracks

2 分析和討論

使用Abaqus有限元軟件進(jìn)行一系列彈塑性大變形分析研究含裂紋加筋板在單軸拉伸載荷下的失效行為。在分析中，使用牛頓-拉普生方法來解決非線性控制方程，考慮上文描述的3種裂紋模型。而且改變裂紋長度和板厚來研究其對含裂紋加筋板極限強(qiáng)度特性的影響。

2.1 長細(xì)比的影響

考慮表3中所列的不同板厚來研究其對含裂紋加筋板極限強(qiáng)度特性的影響。選擇裂紋長度c/b=0.1和裂紋縱向位置固定于加筋板的中橫剖面。將加載試驗中采集的軸向荷載和軸向拉伸變形進(jìn)行換算，即將軸向荷載除以實(shí)測的試件截面積，得到截面的平均應(yīng)力值；將軸向拉伸變形除以試件的實(shí)測初始長度，得到試件的平均應(yīng)變值[15]，圖8給出了不同長細(xì)比的模型（a）應(yīng)力-應(yīng)變曲線（裂紋長度c/b=0.1）。圖中也包含了相關(guān)無裂紋損傷完整加筋板作為參考模型。應(yīng)力-應(yīng)變曲線的峰值即為加筋板的極限拉伸強(qiáng)度，當(dāng)載荷水平達(dá)到該值時，含裂紋加筋板的應(yīng)力云圖（以β=5為例）如圖9所示?？梢钥闯?，當(dāng)載荷水平達(dá)到含裂紋加筋板極限強(qiáng)度時，加筋板中橫剖面幾乎所有單元的等效米塞斯應(yīng)力都達(dá)到材料的抗拉強(qiáng)度，加筋板失去抗拉能力。

為了測量裂紋損傷對加筋板極限拉伸強(qiáng)度的影響，將結(jié)果與相關(guān)無裂紋損傷完整加筋板進(jìn)行對比是必要的。從這一點(diǎn)出發(fā)，引入極限強(qiáng)度縮減系數(shù)，如下式：

圖8 不同長細(xì)比的模型（a）應(yīng)力-應(yīng)變曲線（裂紋長度c/b=0.1）Fig.8 The stress-strain curve of Model (a) with different slenderness (crack length c/b=0.1)

圖9 含裂紋加筋板的應(yīng)力云圖（β=5）Fig.9 The stress distribution of the stiffened-panel with cracks（β=5）

其中： σu為含裂紋加筋板極限強(qiáng)度； σu0為相關(guān)完整加筋板極限強(qiáng)度。

3種長細(xì)比下極限強(qiáng)度縮減系數(shù)分別為0.914，0.902，0.895，含裂紋加筋板極限拉伸強(qiáng)度隨加筋板長細(xì)比的增大略有減小，說明當(dāng)加筋板厚度越大，裂紋損傷對加筋板極限強(qiáng)度的影響越明顯?？傮w來看，加筋板長細(xì)比對極限強(qiáng)度的影響不大，這是因為在本文的算例中，當(dāng)加筋板板厚增大時，裂紋深度也隨之增大，表明含貫穿型裂紋的加筋板在單軸拉伸載荷下的剩余強(qiáng)度對加筋板長細(xì)比不敏感。

2.2 裂紋長度的影響

為了研究裂紋長度對加筋板失效行為的影響，選擇板厚為10 mm，裂紋位于加筋板中橫剖面的工況，改變底板上裂紋的長度。加強(qiáng)筋上沒有裂紋。圖10演示了裂紋長度對應(yīng)力-應(yīng)變關(guān)系的影響?？芍鸭y長度的增加影響加筋板結(jié)構(gòu)響應(yīng)，極限強(qiáng)度和剩余強(qiáng)度都隨裂紋長度的增大而減小。而且，從圖中還可以看出極限應(yīng)變（失效點(diǎn)應(yīng)變）隨裂紋長度的增大而減?。?種工況下極限應(yīng)變分別為0.158，0.152，0.142）。3種裂紋長度下極限強(qiáng)度縮減系數(shù)分別為0.939，0.902，0.793，含裂紋加筋板極限拉伸強(qiáng)度隨裂紋長度的增大而減小，且減小的幅度逐漸增大，說明含裂紋加筋板在單軸拉伸載荷下的極限強(qiáng)度對裂紋長度較為敏感。

圖10 裂紋長度對應(yīng)力-應(yīng)變關(guān)系的影響Fig.10 The influence of crack length on the stress-strain behavior

2.3 裂紋位置的影響

為了檢查裂紋位置對含裂紋加筋板極限強(qiáng)度行為的影響，設(shè)計了裂紋只出現(xiàn)在底板上、裂紋只出現(xiàn)在加強(qiáng)筋上、裂紋同時出現(xiàn)在底板和加強(qiáng)筋上3組仿真試驗。圖11和圖12分別為裂紋在加強(qiáng)筋上和裂紋同時在底板和加強(qiáng)筋上時的應(yīng)力云圖，圖13為3種工況下含裂紋加筋板應(yīng)力-應(yīng)變曲線。3種工況下極限強(qiáng)度縮減系數(shù)分別為0.939（裂紋只出現(xiàn)在底板上）、0.975（裂紋只出現(xiàn)在加強(qiáng)筋上）、0.924（裂紋同時出現(xiàn)在底板和加強(qiáng)筋上），說明加強(qiáng)筋上的裂紋對含裂紋加筋板極限強(qiáng)度的影響小于底板上的裂紋，而裂紋同時出現(xiàn)在底板和加強(qiáng)筋上時對含裂紋加筋板極限拉伸強(qiáng)度的影響最大。這是因為裂紋只位于底板和加強(qiáng)筋上時，裂紋區(qū)域在單軸拉伸載荷下的變形都受到加強(qiáng)筋或底板的限制，而裂紋同時位于底板和加強(qiáng)筋上時，裂紋在底板和加強(qiáng)筋相交的位置能夠自由變形。

圖11 裂紋在加強(qiáng)筋上時的應(yīng)力云圖Fig.11 The stress distribution when the crack locates on the stiffener

圖12 裂紋同時在底板和加強(qiáng)筋上時的應(yīng)力云圖Fig.12 The stress distribution when the crack locates on both the stiffener and the plate

圖13 三種工況下含裂紋加筋板應(yīng)力-應(yīng)變曲線Fig.13 The stress-strain curve of the stiffened-panel with cracks under 3 different cases

3 結(jié) 語

本文建立了一系列不同長細(xì)比、不同裂紋長度、不同裂紋位置的含裂紋加筋板有限元模型，基于J積分理論對含裂紋加筋板在單軸拉伸載荷下的極限強(qiáng)度進(jìn)行了計算分析，得到以下結(jié)論：

1）長細(xì)比對含貫穿型裂紋的加筋板在單軸拉伸載荷下的剩余強(qiáng)度影響不大；

2）含裂紋加筋板極限拉伸強(qiáng)度隨裂紋長度的增大而減小，且減小的幅度逐漸增大，極限拉伸強(qiáng)度對裂紋長度較為敏感；

3）加強(qiáng)筋上的裂紋對含裂紋加筋板極限強(qiáng)度的影響小于底板上的裂紋，而裂紋同時出現(xiàn)在底板和加強(qiáng)筋上時對含裂紋加筋板極限拉伸強(qiáng)度的影響最大。

本文在進(jìn)行含裂紋加筋板極限強(qiáng)度計算時沒有考慮裂紋擴(kuò)展，這或多或少與實(shí)際情況存在一定的差異。裂紋擴(kuò)展對含裂紋加筋板極限強(qiáng)度的影響，是下一步研究的方向。