Übersicht:
Das Schädigungsmodell simuliert den Schädigungsbeginn und das fortschreitende Wachstum für elastisch-sprödes, orthotropes Material, z. B. einen Verbundwerkstoff. Das Modell soll in erster Linie der Simulation von faserverstärkten Verbundwerkstoffen dienen. Das Schädigungsmodell beschreibt drei Phasen des Materialverhaltens:
- das Verhalten des unbeschädigten Materials, das linear-elastisch sein muss
- der Schädigungsbeginn
- das fortschreitende Schädigungswachstum
Im Modell ist die Schädigung in Materialien durch die Degradierung der Steifheit gekennzeichnet. Es wird davon ausgegangen, dass die Fasern der Verbundwerkstoffe parallel sind und in der Richtung der Materialachse 1 des Elements verlaufen. Es werden vier verschiedene Ausfallarten berücksichtigt:
- Faserreißen bei Spannung
- Faserbeulen und -knicken bei Komprimierung
- Matrixriss bei Querschubspannung und Scherung
- Matrixbruch bei Querschubkomprimierung und Scherung
Das Materialverhalten wird definiert als σ = C d ε, wobei C d die elastische Matrix ist, die die Schädigung widerspiegelt und folgendermaßen berechnet wird:
Dabei gilt:
- D=1-(1-d f )(1-dm)ν 21 ν 12 ν
- d f stellt die Faserschädigung dar.
- d m stellt die Matrixschädigung dar.
- d s stellt die Scherschädigung dar.
- E 1 ist das Elastizitätsmodul in der Faserrichtung.
- E 2 ist das Elastizitätsmodul in der Richtung senkrecht zu den Fasern.
- G ist das Schermodul.
- ν 21 und ν 12 sind Poisson-Zahlen.
Schädigungsbeginn
Bei den Kriterien für den Schädigungsbeginn werden vier unterschiedliche Schädigungsmodi berücksichtigt, die sich folgendermaßen darstellen:
|
Modus I: Faserspannung ( |
|
|
Modus II: Faserkomprimierung ( |
|
|
Modus III: Matrixspannung ( |
|
|
Modus IV: Matrixkomprimierung ( |
|
Dabei gilt:
- f I , f II , f III , f IV sind die Variablen, die mit den Kriterien für den Schädigungsbeginn verbunden sind. Ein Wert von 1.0 oder höher gibt an, dass das Anfangskriterium erreicht wurde.
- X T stellt die Zugfestigkeit in Längsrichtung dar.
- X C stellt die Druckfestigkeit in Längsrichtung dar.
- Y T stellt die Zugfestigkeit in Querrichtung dar.
- Y C stellt die Druckfestigkeit in Querrichtung dar.
- S L stellt die Scherfestigkeit in Längsrichtung dar.
- S T stellt die Scherfestigkeit in Querrichtung dar.
- α stellt den Koeffizienten dar, der bestimmt, inwiefern die Scherspannung zum Kriterium für die Faserfestigkeit beiträgt, und 0<=α <=1 festlegt.
, 
und 
sind die Komponenten des effektiven Spannungszugs,
, der sich aus =Mσ berechnet. σ ist die Nominalspannung und M ist der Schädigungsoperator, der wie folgt definiert ist:
Die internen Schädigungsvariablen d f , d m und d s werden von den Schädigungsvariablen d ft , d fc , d mt und d mc entsprechend der vier Modi abgeleitet:
|
d f = d
ft
, wenn d
fc
, wenn |
|
d m = d
mt
, wenn d
mc
, wenn |
d s = 1-(1-d ft )(1-d fc )(1-d mt )(1-d mc )
Die Gleichungen für die Schädigungsvariablen d ft , d fc , d mt , und d mc sind unten angegeben.
Schädigungsentwicklung:
Basierend auf dem Schädigungsbeginn im vorherigen Abschnitt wird hier die Entwicklung der Schädigungsvariablen für jeden Modus definiert. Um die Netzabhängigkeit bei der Materialerweichung zu verringern, werden Äquivalentspannung und Verschiebung eingeführt.
Äquivalentspannungen und Verschiebungen sind für jeden der vier Schädigungsmodi definiert:
|
Faserzugmodus:
|
|
Faserdruckmodus:
|
|
Matrixzugmodus:
|
|
Matrixdruckmodus:
|
wobei δ ij = L c ε ij und L c , die charakteristische Länge, als Quadratwurzel des Integrationspunktbereichs berechnet wird. Der Operator <x> ist der Macauley-Operator, der als <x> = (x +|x|)/2 definiert ist.
Basierend auf den Äquivalentspannungen und Verschiebungen wird die Schädigungsvariable wie folgt berechnet:
wobei i den vier Modi I bis IV entspricht (oder d ft , d fc , d mt und d mc ). Die Ergebnisse liegen im Bereich von 0 bis 1, wobei 0 angibt, dass keine Schädigung vorliegt und das Element eine volle Last tragen kann, und 1 gibt an, dass eine vollständige Schädigung vorliegt und das Element keine Last in der entsprechenden Richtung tragen kann.
δ 0 eq und σ 0 eq sind die äquivalente Anfangsverschiebung und Spannung, bei der das Anfangskriterium erreicht wurde. Sie sind wie folgt definiert:
.
δ f eq = 2G c i /σ 0 eq ist die äquivalente Verschiebung, bei der das Material in diesem Ausfallmodus vollständig beschädigt ist, wobei G c i , die durch den Ausfallmodus i verbrauchte Energie (Spannung oder Komprimierung) die Schädigungsentwicklung bestimmt.
Die maximale Größe des finiten Elements für jedes Schädigungsgesetz ist L c <=2*E i G c i /(X i )2, wobei E i , G c i und X i jeweils das Elastizitätsmodul, die Bruchenergie und die Bruchfestigkeit darstellen. i entspricht wiederum den vier Schädigungsmodi von I bis IV (oder entsprechende Bezeichnung abhängig von der Variablen).
Stabilisierung bei Schädigung und viskose Regularisierung:
Ein verallgemeinertes Duvaut-Lions-Schema zur viskosen Regularisierung wird eingeführt, um Konvergenzschwierigkeiten durch Materialerweichung und Steifigkeitsdegradation in einer impliziten Analyse zu überwinden. Im Regularisierungsschema wird eine viskose Schädigungsvariable definiert als
η i ist der Viskositätskoeffizient, der die Relaxationszeit des viskosen Systems darstellt, und d i ist die Schädigungsvariable, die im nichtviskosen Modell evaluiert wird, wobei i einen Schädigungsmodus angibt (I bis IV, oder ft, fc, mt und mc entsprechend der Variablen). Die regularisierte Schädigungsvariable wird aktualisiert als:
Das beschädigte Verhalten des viskosen Materials wird als σ = C d (dν)ε angegeben. Ein kleiner Viskositätsparameter verbessert in der Regel die Konvergenzrate des Modells im Erweichungsschema, ohne die Ergebnisse zu beeinflussen.
Andere:
Weitere Analysen von Interesse sind:
Ergänzende freie Energiedichte:
Schädigungsenergie:
Viskoser Energieverbrauch:
Referenzen:
- Hashin, Z., A fatigue failure criterion for fiber reinforced materials, Journal of Composite Materials, 7 (1973), 448-464.
- Lapczyk, I. and Hurtado, J.A., Progressive damage modeling in fiber-reinforced materials, Composites: Part A, 38 (2007), 2333-2341.
- Iannucci, L. and Ankersen, J., An energy based damage model for thin laminated composites, Composites Science and Technology, 66 (2006), 934-951.
- Maimi, P., Camanho, P.P., Mayugo, J.A., and Davila, C.G., A continuum damage model for composite laminates: Part I - constitutive model, Mechanics of Materials, 39 (2007), 897-908.
- Maimi, P., Camanho, P.P., Mayugo, J.A., and Davila, C.G., A continuum damage model for composite laminates: Part II - computational implementation and validation, Mechanics of Materials, 39 (2007), 909-919.