丁玉波

（上海飛機設(shè)計研究院，上海 201210）

基于應(yīng)變不變量失效理論的自動校核程序設(shè)計與實現(xiàn)

丁玉波

（上海飛機設(shè)計研究院，上海 201210）

復(fù)合材料的微觀力學(xué)評價方法是近年來的研究熱點。在對Gosse J.H等提出的應(yīng)變不變量失效準(zhǔn)則（SIFT）進(jìn)行研究的基礎(chǔ)上，通過PCL實現(xiàn)了基于該準(zhǔn)則的復(fù)合材料評價自動化程序。該程序同時考慮了機械載荷和熱載荷情況，并建立了代表體元（RVE）的p型單元微觀力學(xué)模型。利用該模型在周期性邊界條件下的力學(xué)響應(yīng)因子對宏觀力學(xué)模型進(jìn)行修正。同時討論了微觀模型修正系數(shù)的收斂性。最后，通過開孔拉伸實例對該程序進(jìn)行了實算，計算結(jié)果合理。本自動校核程序為基于微觀力學(xué)的復(fù)合材料失效評價提供了一種便捷工具。

復(fù)合材料；PCL；代表體元；應(yīng)變不變量；p型單元; 周期邊界條件

0 引言

幾年前，波音公司的研究員Gosse J.H等建立了一種基于物理意義的破壞準(zhǔn)則——應(yīng)變不變量失效準(zhǔn)則[1,2]。該準(zhǔn)則利用材料的臨界固有屬性，如同判定金屬屈服與否的Von-Mises準(zhǔn)則一樣，當(dāng)?shù)刃?yīng)力等于金屬的屈服應(yīng)力時，屈服就會發(fā)生。如果已知屈服應(yīng)力，那么不同幾何形狀，載荷和邊界條件的金屬結(jié)構(gòu)的初始屈服都可以分析。應(yīng)變不變量失效準(zhǔn)則與此類似，如果知道復(fù)合材料的某些臨界屬性，那么就可以判定不同幾何形狀，載荷和邊界條件的復(fù)合材料結(jié)構(gòu)的損傷起始和擴展。

應(yīng)變不變量準(zhǔn)則由于其廣泛的適用性和明確的物理內(nèi)涵而成為近幾年來國內(nèi)外學(xué)者研究的熱點[3～6]。本文在對應(yīng)變不變量失效準(zhǔn)則進(jìn)行研究的基礎(chǔ)上，利用MSC.patran的編程語言PCL編制了可用于評估復(fù)合材料強度的自動化程序，為基于微觀力學(xué)的復(fù)合材料評價提供了一種便捷工具。同時，對于后續(xù)建立在微觀力學(xué)基礎(chǔ)上的復(fù)合材料損傷擴展分析具有實際意義。

1 應(yīng)變不變量失效理論（SIFT）簡介

多年前，科學(xué)家根據(jù)雙軸實驗得出了樹脂的破壞包線，如圖1所示。

由包線可知，在第一象限內(nèi)，這一包線與大多數(shù)延展性金屬所具有的由形狀改變引起的橢圓形破壞包線有一個明顯的差異。因此，聚合物的臨界變形不僅僅只有一種破壞機制即形狀改變。實際上存在兩種破壞機制：一種是體積改變，一種是形狀改變。而與兩種破壞機制相對應(yīng)的兩種臨界材料屬性是體應(yīng)變和等效應(yīng)變。

圖1 聚合物的破壞包線[1,2]

體積應(yīng)變和等效應(yīng)變都是三個主應(yīng)變的函數(shù)，同時都可以表示成應(yīng)變不變量的形式。應(yīng)變不變量的表達(dá)式為：

在式(1)～式(3)中，

ε1、ε2和ε3分別為材料的第一、第二和第三主應(yīng)變；

J1、J2和J3分別為材料的第一、第二和第三應(yīng)變不變量。

體應(yīng)變是三個應(yīng)變不變量之和，表達(dá)式為：

等效應(yīng)變的表達(dá)式為：

因此，應(yīng)變不變量失效準(zhǔn)則表達(dá)式為：

2 應(yīng)變不變量失效理論（SIFT）的程序設(shè)計

應(yīng)變不變量失效準(zhǔn)則是一種基于多尺度的復(fù)合材料評價理論，需要建立宏觀和微觀兩級模型來實現(xiàn)整個校核過程。因此，程序設(shè)計的具體流程如圖2所示。

3 應(yīng)變不變量失效理論（SIFT）的程序?qū)崿F(xiàn)

這里需要著重介紹一下程序?qū)崿F(xiàn)的基本原理和基本方法，具體見下文。

圖2 SIFT程序計算流程

3.1 微觀力學(xué)模型的建立

根據(jù)顯微觀測，單向復(fù)合材料剖面顯微觀察到的纖維排列如圖3(a)所示。通常利用代表體元（RVE）技術(shù)，把實際的纖維排列按方形（如圖3(b)所示）或六邊形排列（如圖3(c)所示）進(jìn)行簡化。本文選取方形排列的邊長為1的立方單胞模型進(jìn)行程序設(shè)計。

圖3 顯微觀測結(jié)果及其簡化

在使用應(yīng)變不變量破壞準(zhǔn)則即式(6)、式(7)之前，首先要對各項同性的宏觀應(yīng)變場進(jìn)行微觀力學(xué)修正。這里的宏觀應(yīng)變場指通過層壓板理論得到的應(yīng)變場。本文通過采用代表體元法提取微觀修正因子，來對宏觀應(yīng)變場進(jìn)行微觀力學(xué)修正。常用的代表體元法模型分單胞和多胞兩種。這里采用具有一定h細(xì)化的p型單元的單胞模型進(jìn)行分析，模型如圖4所示。

圖4 單胞模型

3.2 應(yīng)變修正系數(shù)的確定

為建立單向復(fù)合材料宏觀-微觀應(yīng)變轉(zhuǎn)換關(guān)系，David[8]等曾提出如下計算公式：

由式(8)可知，建立宏觀-微觀應(yīng)變關(guān)系的關(guān)鍵是要確定Mij(i=1～6，j=1～6)和Ai(i=1～6)的值。為實現(xiàn)這一目的，分兩種情況進(jìn)行分析：

1)為確定Mij(i=1～6，j=1～6)的值，須令△T=0。然后依次令且于是，確定了單胞模型的力邊界條件[8～10]，表1給出了1時的邊界條件，其余類比可得。

(2)為確定Ai(i=1～6)的值，則△T≠0，須使于是，確定了單胞模型的熱邊界條件[8～10]，表2給出了溫度差△T下的邊界條件。

3.3 應(yīng)變修正系數(shù)收斂性的驗證

表1 力邊界條件

表2 熱邊界條件

圖5 各邊的種子數(shù)

圖6 M21的收斂曲線

圖7 M22的收斂曲線

圖8 M23的收斂曲線

圖9 M31的收斂曲線

圖10 M32的收斂曲線

圖11 M33的收斂曲線

圖12 M44的收斂曲線

圖13 M55的收斂曲線

圖14 M66的收斂曲線

圖15 A1的收斂曲線

圖16 A2的收斂曲線

圖17 A3的收斂曲線

4 程序應(yīng)用實例

將自動校核程序在開孔拉伸模型中進(jìn)行應(yīng)用，微觀力學(xué)模型材料屬性如表3所示，宏觀模型材料屬性可根據(jù)纖維體積分?jǐn)?shù)和微觀材料屬性計算得到。開孔拉伸模型修正前后得到的J1和εequivalent分別如圖18～圖21所示，經(jīng)分析比對可知修正結(jié)果無誤。最后，再根據(jù)試驗確定的和εcr，按照式(6)、式(7)來判定結(jié)構(gòu)是否vonmises失效。

表3 纖維和基體材料屬性

圖18 修正前J1的云圖

【】【】

圖19 修正前εequivalent的云圖

圖20 修正后J1云圖

圖21 修正前εequivalent的云圖

5 結(jié)論

通過本文的研究得到以下結(jié)論：

1）利用PCL語言實現(xiàn)整個校核過程的自動化，為基于微觀力學(xué)的復(fù)合材料失效評價提供了一種便捷工具；

2）基于p單元的微觀力學(xué)模型具有良好的收斂性，在PCL程序中使用這種單元取得了較好的效果；

3）通過PCL編制自動校核程序?qū)崿F(xiàn)了初步的復(fù)合材料失效評估，但是其中失效判據(jù)中的臨界值還需要通過大量試驗來確定；

4）基于SIFT的自動校核程序的實現(xiàn)，為后續(xù)進(jìn)行復(fù)合材料的損傷擴展研究具有實際意義。

[1] Gosse J.H.,Christensen S.Strain Invariant Failure Criteria for Polymers in Composites[Z].AIAA-2001-1184,2001.

[2] Jon H. Gosse, Strain Invariant Failure Criteria for Fiber Reinforced Polymeric Composite Materials[C].Proceedings of 13th ICCM conference,2001.

[3] Tay T.E., Yudhanto A, Tan VBC. A Micro Macro Approach to Modeling Progressive Damage in Composite Structures[J]. International Journal of Damage Mechanics,2008,17(5):5-2.

[4] Cai Hongneng, Miyano Yasushi, Nakada Masayuki. Long-term open-hole compression strength of CFRP laminates based on strain invariant failure theory[J].Journal of Thermoplastic Composite Materials,2009,22:63-81.

[5] L. J. Hart-Smith.Application of the strain invariant failure theory (SIFT) to metals and fiber-polymer composites[J],Philosophical Magazine,2010,90,4263-4331.

[6] 黎增山,關(guān)志東,等.國產(chǎn)碳纖維增強樹脂基復(fù)合材料應(yīng)變不變量性能[J].復(fù)合材料學(xué)報,2011,28(5):192-1968.

[7] Carl T.Herakovich.Mechanics of Fibrous Composites[M],Wiley, 1998,10-11.

[8] David L.B.,Gosse J.H.,Jeffrey A.W.David L.Buchanan, Micromechanical Enhancement of the Macroscopic Strain State for Advanced Composite Materials[J].Composites Science and Technology,2009,69,1974-1978.

[9] Xia Z., Zhang Y.,Ellyin F.A Unified Periodical Boundary Conditions for Representative Volume Elements of Composites and Applications[J].International Journal of Solids and Structures,2003,40(8):1907-1921.

[10] Jin K., Huang Y.,Lee Y.,Ha S. Distribution of Micro Stresses and Interfacial Tractions in Unidirectional Composites[J].Journal of Composite Materials,2008,42(18):1825-1849.

Design and achievement of automatic check program based on SIFT

DING Yu-bo

O343；TB332

A

1009-0134(2017)04-0011-06

2017-02-23

丁玉波（1983 -），男，內(nèi)蒙古通遼人，助理工程師，碩士研究生，研究方向為飛機結(jié)構(gòu)強度分析和適航驗證技術(shù)。