abaqus中復合材料的累積損傷及失效

1、-. z.纖維增強材料的累積損傷與失效：Abaqus擁有纖維增強材料的各向異性損傷的建模功能纖維增強材料的損傷與失效概論，19.3.1節(jié)。假設未損傷材料為線彈性材料。因為該材料在損傷的初始階段沒有大量的塑性變形，所以用來預測纖維增強材料的損傷行為。Hashin標準最開場用來預測損傷的產(chǎn)生，而損傷演化規(guī)律基于損傷過程和線性材料軟化過程中的能量耗散理論。另外，Abaqus也提供混凝土損傷模型，動態(tài)失效模型和在粘著單元以及連接單元中進展損傷與失效建模的專業(yè)功能。本章節(jié)給出了累積損傷與失效的概論和損傷產(chǎn)生與演變規(guī)律的概念簡介，并且僅限于塑性金屬材料和纖維增強材料的損傷模型。損傷與失效模型的通用框架Ab

2、aqus提供材料失效模型的通用建模框架，其中允許同一種的材料應用多種失效機制。材料失效就是由材料剛度的逐漸減弱而引起的材料承當載荷的能力完全喪失。剛度逐漸減弱的過程采用損傷力學建模。為了更好的了解Abaqus中失效建模的功能，考慮簡單拉伸測試中的典型金屬樣品的變形。如圖19.1.1-1中所示，應力應變圖顯示出明確的劃分階段。材料變形的初始階段是線彈性變形a-b段，之后隨著應變的加強，材料進入塑性屈服階段b-c段。超過c點后，材料的承載能力顯著下降直到斷裂c-d段。最后階段的變形僅發(fā)生在樣品變窄的區(qū)域。C點說明材料損傷的開場，也被稱為損傷開場的標準。超過這一點之后，應力-應變曲線c-d由局部變形

3、區(qū)域剛度減弱進展決定。根據(jù)損傷力學可知，曲線c-d可以看成曲線c-d的衰減，曲線c-d是在沒有損傷的情況下，材料應該遵循的應力-應變規(guī)律曲線。圖19.1.1-1 金屬樣品典型的軸向應力-應變曲線因此，在Abaqus中失效機制的詳細說明里包括四個明顯的局部：材料無損傷階段的定義如圖19.1.1-1中曲線a-b-c-d損傷開場的標準如圖19.1.1-1中c點損傷開展演變的規(guī)律如圖19.1.1-1中曲線c-d單元的選擇性刪除，因為一旦材料的剛度完全減退就會有單元從計算中移除如圖19.1.1-1中的d點。 關(guān)于這幾局部的內(nèi)容，我們會對金屬塑性材料金屬塑性材料的損傷與失效概論，19.2.1節(jié)和纖維增強材

4、料纖維增強符合材料的損傷與失效概論，19.3.1節(jié)進展分開討論。網(wǎng)格依賴性在連續(xù)介質(zhì)力學中，通常是根據(jù)應力-應變關(guān)系建立材料本構(gòu)模型。當材料表現(xiàn)出導致應變局部化的應變軟化行為時，有限元分析的結(jié)果帶有強烈的網(wǎng)格依賴性，能量的耗散程度取決于網(wǎng)格的精簡程度。在Abaqus中所有可使用損傷演化模型都使用減輕網(wǎng)格依賴性的公式。這是通過在公式中引入特征長度來實現(xiàn)的，特征長度作為一個應力-位移關(guān)系可以表達本構(gòu)關(guān)系中軟化局部，它與單元尺寸有關(guān)系。在此情況下，損傷過程中耗散的能量不是由每個單位體積衡量，而是由每個單位面積衡量。這個能量值作為另外一個材料參數(shù)，用來計算材料發(fā)生完全損傷時的位移。這是與材料斷裂力學中

5、臨界能量釋放率的概念一致的。此公式確保了適宜能量的耗散以及最大程度減輕網(wǎng)格的依賴。19.3 纖維增強復合材料的損傷與失效纖維增強復合材料的損傷與失效：概論 19.3.1節(jié)纖維增強復合材料的損傷初始準則 19.3.2節(jié)纖維增強復合材料的損傷演化規(guī)律與單元移除 19.3.3節(jié)19.3.1纖維增強復合材料的損傷與失效：概論產(chǎn)品：Abaqus/Standard Abaqus/E*plicit Abaqus/CAE參考：*DAMAGE INITIATION*DAMAGE EVOLUTION*DAMAGE STABILIZATION“Hashin damage in “Dening damage, Sec

6、tion 12.8.3 of the Abaqus/CAE Users Manual, in the online HTML version of this manual概論Abaqus具有為纖維增強復合材料的漸進損傷和破壞建模的能力。它能預測各向異性彈-脆性材料的損傷產(chǎn)生與演化規(guī)律。材料模型要求以下定義：未損傷時的材料屬性必須是線彈性的損傷初始產(chǎn)生準則損傷演化規(guī)律單元選擇性移除單向板的損傷根本概念損傷的特點是材料剛度的逐漸減小。這在纖維增強復合材料的分析中有很重要的作用。很多這樣的材料表現(xiàn)出彈-脆性行為，也就是材料在小變形的情況下就開場發(fā)生損傷。所以在建立此種材料的模型時，材料塑性被忽略。假

7、設纖維增強復合材料中的纖維是平行的，如圖19.3.1-1所示。我們必須在用戶定義的局部坐標系中定義材料屬性。單向?qū)游挥?-2平面內(nèi)，1方向表示纖維方向。我們要用定義正交線彈性材料的方法來定義材料未損傷時的行為。最簡單的方法是定義平面應力的正交材料。然而，材料行為也可以采用定義工程常數(shù)或直接定義彈性剛度矩陣的方法來定義。圖19.3.1-1 單向?qū)覣baqus支持的各向異性損傷模型基于Matzenmiller et. al(1995),Hashin and Rotem(1973),Hashin(1980),and Camanho and Davila(2002)的工作。四種不同的失效模型：拉伸載荷

8、作用下的纖維斷裂壓縮載荷下的纖維屈曲和扭結(jié)橫向拉伸和剪切載荷下的基體斷裂橫向壓縮和剪切載荷下的基體破碎在Abaqus中，損傷萌生是由Hashin1980和Rotem1973提出的損傷初始準則來決定的，準則中的失效面是由有效應力空間來表示的可以有效承受力載荷的面上的應力。這些準則的細節(jié)將在19.3.2節(jié)“纖維增強復合材料的損傷產(chǎn)生中討論。材料的應力根據(jù)下式計算，式中表示應變，表示彈性矩陣并反映任何損傷，有以下形式：式中，反映當前纖維損傷狀態(tài)，反映當前基體損傷狀態(tài)，反映當前剪切損傷狀態(tài)，為纖維方向的楊氏模量，為垂直于纖維方向上的楊氏模量，為剪切模量，為泊松比。決定損傷彈性矩陣的演化將在19.3.3

9、節(jié)“纖維增強復合材料的損傷演化與單元移除中詳細介紹，19.3.3節(jié)還將介紹：處理嚴重損傷的選擇最大變形與單元的選擇性移除在19.3.3節(jié)“纖維增強復合材料的損傷演化與單元移除粘滯阻力在19.3.3節(jié)“纖維增強復合材料的損傷演化與單元移除中的“粘滯阻力單元纖維增強復合材料的損傷模型必須采用平面應力單元，包括平面應力單元、殼單元、連續(xù)殼單元和薄膜單元。其他參考Hashin, Z., and A. Rotem, “A Fatigue Criterion for Fiber-Reinforced Materials, Journal of posite Materials, vol. 7, pp. 4

10、48464, 1973.Hashin, Z., “Failure Criteria for Unidirectional Fiber posites, Journal of Applied Mechanics,vol. 47, pp. 329334, 1980.Matzenmiller, A., J. Lubliner, and R. L. Taylor, “A Constitutive Model for Anisotropic Damage in Fiber-posites, Mechanics of Materials, vol. 20, pp. 125152, 1995.Camanho

11、, P. P., and C. G. Davila, “Mi*ed-Mode Decohesion Finite Elements for the Simulation of Delamination in posite Materials, NASA/TM-2002211737, pp. 137, 2002.19.3.2 纖維增強復合材料的損傷萌生產(chǎn)品：Abaqus/Standard Abaqus/E*plicit Abaqus/CAE參考：*DAMAGE INITIATION“Hashin damage in “Dening damage, Section 12.8.3 of the Ab

12、aqus/CAE Users Manual, in the online HTML version of this manual概論纖維增強材料的損傷建模功能：要求材料未損傷時為線彈性參考“線彈性行為第17.2.1節(jié)基于Hashin的理論Hashin和Rotem，1973，和Hashin,1980考慮四種不同的失效模型：纖維拉伸、纖維壓縮、基體斷裂和基體破碎可以與19.3.3節(jié)“纖維增強復合材料的損傷演化與單元移除中提到的損傷演化模型一起使用。損傷萌生 損傷萌生是在材料硬點退化開場。在Abaqus中纖維增強復合材料的損傷萌生準則基于Hashin的理論。這些準則考慮了四種不同的損傷萌生機制：纖維

13、拉伸，纖維壓縮，基體斷裂和基體破碎。損傷萌生準則有下面的一般形式：在上面的方程中*T表示縱向拉伸強度; *C表示縱向抗壓強度; YT表示橫向拉伸強度; YC表示橫向抗壓強度; SL表示縱向剪切強度; ST表示橫向剪切強度; 是一個系數(shù)用于決定剪應力對纖維拉伸損傷準則的影響；是有效應力*量的分量，是用來評估萌生標準并按下式計算：上式中是名義應力，M是損傷矩陣：df，dm和ds是內(nèi)部損傷變量分別代表纖維，基體和剪切損傷，這是由損傷變量dft，dfc，dmt，dmc推導出的，用于對應先前所討論的四個模式，如下：在任何損傷萌生和演化之前，損傷控制矩陣M為單位矩陣，因此。一旦至少有一個模型已經(jīng)發(fā)生損傷萌

14、生和演化，損傷控制矩陣在損傷萌生準則中就具有重要意義見“纖維增強復合材料的損傷演化和元素去除，第19.3.3損傷演化的討論。有效應力用于表示有效承載力載荷的損傷面上的應力。上面介紹的損傷萌生準則，可以通過設置=0.0和ST=YC/2獲得Hashin和ROTEM1973提出的模型，設置=1.0獲得Hashin1980提出的模型。與每個損傷萌生準則有關(guān)的輸出變量纖維拉伸，纖維壓縮，基體拉伸，基體壓縮用來表示是否已經(jīng)到達標準。值為1.0或更高則說明萌生準則已滿足。如果只定義損傷萌生模型而沒有定義相關(guān)的演化規(guī)律，萌生準則將只影響輸出。因此，可以在不建立損傷演化模型的情況下使用損傷萌生準則來評估材料的特

15、性。輸入文件使用方法：使用以下選項來定義Hashin損傷萌生標準：*DAMAGEINITIATION,CRITERION=HASHIN,ALPHA=*T，*C，YT，YC，SL, STAbaqus / CAE用法：Property module: material editor: MechanicalDamage for Fiber-Reinforced positesHashin Damage單元：損傷萌生準則必須與平面應力單元一起使用，其中包括平面應力單元，殼單元，連續(xù)殼單元，薄膜單元。輸出：在Abaqus中除了標準輸出標識符“Abaqus /標準輸出變量標識符，“第4.2.1節(jié)，纖維增強

16、復合材料損傷模型中材料*點萌生損傷時涉及的輸出變量還包括：DMICRT 所有損傷引發(fā)的標準組件。HSNFTCRT 在分析過程中纖維拉伸萌生準則所經(jīng)歷的最大值。HSNFCCRT 在分析過程中纖維壓縮萌生準則所經(jīng)歷的最大值。HSNMTCRT 在分析過程基體斷裂萌生準則所經(jīng)歷的最大值。HSNMCCRT 在分析過程基體破碎萌生準則所經(jīng)歷的最大值。對于上述變量說明在損傷模型中損傷萌生準則是否已滿足，如果這個值小于1.0，則說明損傷萌生準則還沒有得到滿足，假設這個值大于等于1.0表示該準則已滿足。如果定義了一個損傷演化模型，則這個變量的最大值不會超過1.0。如果你沒有定義損傷演化模型，這個變量可以大于1.

17、0，這個值說明已超出準則的數(shù)量。其他參考Hashin, Z., and A. Rotem, “A Fatigue Criterion for Fiber-Reinforced Materials, Journal of posite Materials, vol. 7, pp. 448464, 1973.Hashin, Z., “Failure Criteria for Unidirectional Fiber posites, Journal of Applied Mechanics,vol. 47, pp. 329334, 198019.3.3 纖維增強復合材料的損傷演化與單元去除產(chǎn)品：A

18、baqus/Standard Abaqus/ E*plicit Abaqus/ CAE參考文獻 *DAMAGE EVOLUTION “Damage evolution in “Defining damage, Section 12.8.3 of the Abaqus/CAE Users Manual, in the online HTML version of this manual概論 Abaqus中的纖維增強復合材料的損傷演化能力：假定損傷的特點是材料的剛度逐步退化，導致材料失效;要求未損傷材料為線彈性行為;考慮到四個不同的失效模式：纖維拉伸失效，纖維壓縮失效，基體斷裂失效，基體破碎;使用

19、四個損傷變量來描述每個失效模式的損傷;必須結(jié)合Hashin的損傷萌生準則一起使用;基于損傷過程中的能量消耗;提供幾種失效時行為的選擇，包括從網(wǎng)格中去除單元；;可以結(jié)合本構(gòu)方程的粘滯阻力一起使用來提高軟化機制的收斂速度。損傷演化 上一節(jié)討論了纖維增強復合材料在平面應力作用下的損傷萌生。本節(jié)將討論材料在定義損傷演化模型情況下的損傷萌生前的行為。損傷萌生前，材料是線性彈性，有一個平面應力各向異性材料的剛度矩陣。此后，材料行為依據(jù)下式計算：這里是應變和Cd是損傷彈性矩陣，其形式是：其中，反映了纖維損傷的現(xiàn)狀，反映當前基體的損傷狀態(tài)，反映剪切損傷的現(xiàn)狀，是在纖維方向的彈性模量，是垂直方向的彈性模量，是剪

20、切模量，和是泊松比。損傷變量df，dm和ds是由損傷變量dft，dfc，dmt，dmc推導出的，對應先前所討論的四個失效模式，如下：和是有效應力*量的分量。有效應力*量主要用來評估損傷萌生準則見“纖維增強復合材料的損傷開場第19.3.2節(jié)關(guān)于有效應力*量計算的說明。每種模型中損傷變量的演化為減輕材料軟化的過程中的網(wǎng)格依賴性，Abaqus在計算中引入特征長度，使本構(gòu)關(guān)系表示為應力-位移關(guān)系。損傷變量發(fā)生變化得到圖19.3.3-1所示的應力-位移曲線，損傷開場之前正斜率的應力-位移曲線對應線性彈性材料的狀態(tài);損傷萌生后的負斜率，是根據(jù)相應損傷變量的變化得到的，損傷變量根據(jù)下面所示的方程計算。四種失

21、效模式下的等效位移和應力，定義如下：纖維拉伸：纖維壓縮：基體斷裂：基體破碎：對于平面應力單元，特征長度Lc作為一個集成點面積的平方根計算。上述方程中的符號代表麥考括號運算符，定義為： 當時，損傷萌生即后的行為，如圖19.3.3-1所示，一個特定的模型的損傷變量由下式給出其中是到達萌生準則時的初始等效位移，是材料完全損傷時的位移。其關(guān)系顯示在圖19.3.3-2。 圖19.3.3-2 損傷變量與等效位移函數(shù)關(guān)系圖 各種模型的值取決于材料的彈性剛度和強度參數(shù)，二者是損傷萌生定義的一局部見“纖維增強復合材料的損傷開場，第19.3.2節(jié)。對于每個失效模式，必須定義導致失效的能量耗散Gc，其對應于圖19.

22、3.3-3中三角形OAC的面積。各種模式下的值取決于各自的Gc值。當材料從局部損壞狀態(tài)下卸載，如在圖19.3.3-3B點，沿直線路徑走回等效應力與等效位移的原點，在重新加載情況下，以一樣的路徑回到B點，如以下圖所示。 圖19.3.3-3 線性損傷演化輸入文件的使用：使用下面的選項定義損傷演化規(guī)律：*DAMAGE EVOLUTION,TYPE=ENERGY,SOFTENING=LINEAR其中分別代表纖維拉伸、纖維壓縮、基體斷裂和基體破碎四種損傷中的能量耗散。Abaqus/CAE的使用 ：Property module: material editor: MechanicalDamage for

23、 Fiber-Reinforced positesHashin Damage:SuboptionsDamage Evolution: Type: Energy: Softening: Linear最大衰減和單元的選擇性移除我們已經(jīng)掌握了在嚴重損傷時Abaqus處理單元的行為。默認情況下，材料*一點的損傷變量的上限為。我們可以按照21.1.4節(jié)“部件控制中的“單元移除和材料損傷演化過程中的最大衰減的控制中所討論的，來減小上限值。在Abaqus/Standard中，默認情況下，一旦所有失效模型的所有材料硬點到達，單元就會被移除參考21.1.4節(jié)“部件控制中的“單元移除和材料損傷演化過程中的最大衰減

24、的控制。在Abaqus/E*plicit中，當纖維失效模型拉伸或者壓縮相關(guān)的損傷變量到達時，假定材料硬點失效，并且當一個單元的所有集成位置的剖分硬點都滿足此條件時，單元將從網(wǎng)格中移除。例如，在任何一個集成位置處，通過厚度剖分硬點的殼單元在從網(wǎng)格中移除之前，必須是失效狀態(tài)。一個單元移除之后，單元的輸出變量STATUS就置為0，他不再支持隨后產(chǎn)生的變形。然而，在Abaqus中單元依然存在，在Abaqus/CAE的視圖模式是可見的。可以通過創(chuàng)立一個包括所有STATUS=1.0的單元的顯示組來關(guān)閉已經(jīng)移除的單元的顯示參考Abaqus/CAE用戶手冊中的第57.2.2節(jié)“創(chuàng)立和編輯顯示組的選擇方法。你也

25、可選擇性的設置，即使所有的損傷變量都到達，單元依然保存在模型中。在此情況下，一旦所有損傷變量到達最大值，剛度變成一個常數(shù)參照此節(jié)之前的表述。在Abaqus/Standard中單元移除的相關(guān)困難單元從模型中移除之后，其節(jié)點即使沒有與任何激活的單元連接，仍然存在于模型之中。當求解進展時，這些節(jié)點可能會產(chǎn)生非物理位移而導致外推算法，Abaqus/Standard用此來加快求解速度參考“非線性問題的收斂準則節(jié)。這些非物理性位移可以通過關(guān)閉外推法來防止。另外，如果將點載荷施加在不與激活單元連接的節(jié)點上，將會因為沒有剛度承載此載荷造成收斂困難的問題。粘性阻力在隱式分析中，比方Abaqus/Standard

26、，材料表現(xiàn)的軟化行為和剛度衰減經(jīng)常導致嚴重的收斂困難，?？梢允褂谜硿枘釁?shù)來克制一些收斂困難的問題，粘滯阻尼將軟化材料的切向剛度矩陣對于足夠小的時間增量是正的。在此粘性方法中，通過演化方程定義了粘性損傷變量：式中是表示粘性系統(tǒng)弛豫時間的粘性系數(shù)，d是在非粘性骨干模型中變化的損傷變量。粘性材料的損傷行為按下式給出：式中損傷塑性矩陣是根據(jù)每種失效模型的損傷變量的粘性值來計算。粘性系數(shù)很小的粘性阻尼可以提高處于軟化過程中模型的收斂速度。根本的思想是粘性系統(tǒng)的求解與非粘性系統(tǒng)比擬顯得緩和，因為，式中t表示時間。在Abaqus/E*plicit中，粘滯阻尼也是可用的。粘滯阻尼降低了損傷增長率，導致斷裂能量增加和形變率增加，可以作為一種建立依賴變化率的材料行為模型的有效方法。在Abaqus/Standard中，整個模型或者一個單元上與粘滯阻尼相關(guān)的能量集使用輸出變量ALLCD是可見的。粘滯阻尼系數(shù)的定義對于不同的失效模型可以定