View source for Micro-mechanics of failure

[[Image:Hierarchy of micromechanics-based analysis procedure for composite structures.png|thumb|right|400px|alt=Hierarchy of micromechanics-based analysis procedure for composite structures|Hierarchy of micromechanics-based analysis procedure for composite structures.]]

The theory of '''micro-mechanics of failure''' aims to explain the [[material failure|failure]] of  [[fibre-reinforced plastic|continuous fiber reinforced composites]] by micro-scale analysis of stresses within each constituent material (such as fiber and matrix), and of the stresses at the interfaces between those constituents, calculated from the macro stresses at the ply level.<ref name=SKHa/>

As a completely mechanics-based failure theory, the theory is expected to provide more accurate analyses than those obtained with phenomenological models such as [[Tsai-Wu failure criterion|Tsai-Wu]]<ref name=SWTsai/> and Hashin<ref name=Hashin1/><ref name=Hashin2/> failure criteria, being able to distinguish the critical constituent in the critical ply in a composite laminate.

[[Image:Comparison between theoretical failure predictions and test data.png|thumb|right|300px|alt=Comparison between theoretical failure predictions and test data|Failure envelopes generated by MMF and the Tsai-Wu failure criterion for a carbon/epoxy UD ply, with test data superimposed. Failed constituent envelopes are predicted by MMF but not by Tsai-Wu.]]

==Basic concepts==
The basic concept of the micro-mechanics of failure (MMF) theory is to perform a hierarchy of micromechanical analyses, starting from mechanical behavior of constituents (the fiber, the matrix, and the interface), then going on to the mechanical behavior of a ply, of a laminate, and eventually of an entire structure.

At the constituent level, three elements are required to fully characterize each constituent:
* The [[Constitutive equation|constitutive relation]], which describes the transient, or time-independent, response of the constituent to external mechanical as well as hygrothermal loadings;
* The [[master curve (material science)|master curve]], which describes the time-dependent behavior of the constituent under creep or fatigue loadings;
* The [[material failure theory|failure criterion]], which describes conditions that cause failure of the constituent.
The constituents and a unidirectional lamina are linked via a proper micromechanical model, so that ply properties can be derived from constituent properties, and on the other hand, micro stresses at the constituent level can be calculated from macro stresses at the ply level.

== Unit cell model ==

[[Image:Schematic illustration of idealized fiber arrays and their corresponding unit cells.png|thumb|right|300px|alt=Schematic illustration of idealized fiber arrays and their corresponding unit cells|Schematic illustration of idealized fiber arrays and their corresponding unit cells.]]

Starting from the constituent level, it is necessary to devise a proper method to organize all three constituents such that the microstructure of a UD lamina is well-described. In reality, all fibers in a UD ply are aligned longitudinally; however, in the cross-sectional view, the distribution of fibers is random, and there is no distinguishable regular pattern in which fibers are arrayed. To avoid such a complication cause by the random arrangement of fibers, an idealization of the fiber arrangement in a UD lamina is performed, and the result is the regular fiber packing pattern. Two regular fiber packing patterns are considered: the square array and the hexagonal array. Either array can be viewed as a repetition of a single element, named unit cell or [[representative volume element]] (RVE), which consists of all three constituents. With periodical boundary conditions applied,<ref name=ZXia/> a unit cell is able to respond to external loadings in the same way that the whole array does. Therefore, a unit cell model is sufficient in representing the microstructure of a UD ply.

== Stress amplification factor (SAF) ==

Stress distribution at the laminate level due to external loadings applied to the structure can be acquired using [[Finite element method|finite element analysis (FEA)]]. Stresses at the ply level can be obtained through transformation of laminate stresses from laminate coordinate system to ply coordinate system. To further calculate micro stresses at the constituent level, the unit cell model is employed. Micro stresses <math>\sigma</math> at any point within fiber/matrix, and micro surface tractions <math>t</math> at any interfacial point, are related to ply stresses <math>\bar{\sigma}</math> as well as temperature increment <math>\Delta T</math> through:<ref name=KKJin/>

:<math>
\begin{array}{lcl}
\sigma_{\mathrm{f}}&=&M_{\mathrm{f}}\bar{\sigma} + A_{\mathrm{f}}\Delta T\\
\sigma_{\mathrm{m}}&=&M_{\mathrm{m}}\bar{\sigma} + A_{\mathrm{m}}\Delta T\\
t_{\mathrm{i}}&=&M_{\mathrm{i}}\bar{\sigma} + A_{\mathrm{i}}\Delta T
\end{array} </math>

Here <math>\sigma</math>, <math>\bar{\sigma}</math>, and <math>t</math> are column vectors with 6, 6, and 3 components, respectively. Subscripts serve as indications of constituents, i.e. <math>{\mathrm{f}}</math> for fiber, <math>{\mathrm{m}}</math> for matrix, and <math>{\mathrm{i}}</math> for interface. <math>M</math> and <math>A</math> are respectively called stress amplification factors (SAF) for macro stresses and for temperature increment. The SAF serves as a conversion factor between macro stresses at the ply level and micro stresses at the constituent level. For a micro point in fiber or matrix, <math>M</math> is a 6×6 matrix while <math>A</math> has the dimension of 6×1; for an interfacial point, respective dimensions of <math>M</math> and <math>A</math> are 3×6 and 3×1. The value of each single term in the SAF for a micro material point is determined through [[Finite element method|FEA]] of the unit cell model under given macroscopic loading conditions. The definition of SAF is valid not only for constituents having [[Linear elasticity|linear elastic]] behavior and constant [[Coefficient of thermal expansion|coefficients of thermal expansion (CTE)]], but also for those possessing complex [[Constitutive equation|constitutive relations]] and variable [[Coefficient of thermal expansion|CTEs]].

== Constituent failure criteria ==

=== Fiber failure criterion ===

Fiber is taken as transversely isotropic, and there are two alternative failure criteria for it:<ref name=SKHa/> a simple maximum stress criterion and a quadratic failure criterion extended from [[Tsai-Wu failure criterion]]:

:<math>
\begin{array}{lcl}
\text{Maximum stress failure criterion:}-X^\prime_{\mathrm{f}} < \sigma_1 < X_{\mathrm{f}}\\
\text{Quadratic failure criterion: }\displaystyle\sum_{j=1}^6\displaystyle\sum_{i=1}^6 F_{ij}\sigma_i\sigma_j + \displaystyle\sum_{i=1}^6 F_i\sigma_i = 1
\end{array}</math>

The Coefficients involved in the quadratic failure criterion are defined as follows:
:<math>
F_{11} = \cfrac{1}{X_{\mathrm{f}}X^\prime_{\mathrm{f}}}\ ,\ F_{22} = F_{33} = \cfrac{1}{Y_{\mathrm{f}}Y^\prime_{\mathrm{f}}}
</math>
:<math>
F_{44} = \cfrac{1}{S_{\mathrm{f}4}^2}\ ,\ F_{55} = F_{66} = \cfrac{1}{S_{\mathrm{f}6}^2}
</math>
:<math>
F_{1} = \cfrac{1}{X_{\mathrm{f}}} - \cfrac{1}{X_{\mathrm{f}}^\prime}\ ,\ F_{2} = F_{3} = \cfrac{1}{Y_{\mathrm{f}}} - \cfrac{1}{Y_{\mathrm{f}}^\prime}
</math>
:<math>
F_{12} = F_{21} = F_{13} = F_{31} = -\cfrac{1}{2\sqrt{X_{\mathrm{f}} {X}_{\mathrm{f}}^\prime Y_{\mathrm{f}}Y_{\mathrm{f}}^\prime}}\ ,\ F_{23} = F_{32} = -\cfrac{1}{2Y_{\mathrm{f}}Y_{\mathrm{f}}^\prime}
</math>
where <math>X_{\mathrm{f}}</math>, <math>X_{\mathrm{f}}^\prime</math>, <math>Y_{\mathrm{f}}</math>, <math>Y_{\mathrm{f}}^\prime</math>, <math>S_{\mathrm{f}4}</math>, and <math>S_{\mathrm{f}6}</math> denote longitudinal tensile, longitudinal compressive, transverse tensile, transverse compressive, transverse (or through-thickness) shear, and in-plane shear strength of the fiber, respectively.

Stresses used in two preceding criteria should be micro stresses in the fiber, expressed in such a coordinate system that 1-direction signifies the longitudinal direction of fiber.

=== Matrix failure criterion ===

The polymeric matrix is assumed to be isotropic and exhibits a higher strength under uniaxial compression than under uniaxial tension. A modified version of [[Von Mises yield criterion|von Mises failure criterion]] suggested by Christensen<ref name=Christen/> is adopted for the matrix:

:<math>
\begin{array}{lcl}
\cfrac{\sigma_{Mises}^2}{C_{\mathrm{m}}T_{\mathrm{m}}} + \left(\cfrac{1}{T_{\mathrm{m}}} - \cfrac{1}{C_{\mathrm{m}}}\right)I_1 = 1 
\end{array}</math>

Here <math>{T}_{\mathrm{m}}</math> and <math>{C}_{\mathrm{m}}</math> represent matrix tensile and compressive strength, respectively; whereas <math>\sigma_{Mises}</math> and <math>{\mathrm{I}}_1</math> are [[Stress (mechanics)#Invariants of the stress deviator tensor|von Mises equivalent stress]] and [[Stress (mechanics)#Principal stresses and stress invariants|the first stress invariant]] of micro stresses at a point within matrix, respectively.

=== Interface failure criterion ===

The fiber-matrix interface features traction-separation behavior, and the failure criterion dedicated to it takes the following form:<ref name=Camanho/>

<math>
\begin{array}{lcl}
\left(\cfrac{\left\langle{t}_{n}\right\rangle}{{Y}_{n}}\right)^2 + \left(\cfrac{{t}_{s}}{{Y}_{s}}\right)^2 = 1
\end{array}</math>

where <math>{t}_{n}</math> and <math>{t}_{s}</math> are normal (perpendicular to the interface) and shear (tangential to the interface) interfacial tractions, with <math>{Y}_{n}</math> and <math>{Y}_{s}</math> being their corresponding strengths. The angle brackets ([[Macaulay brackets]]) imply that a pure compressive normal traction does not contribute to interface failure.

== Further extension of MMF ==
=== Hashin’s Failure Criteria ===
	These are interacting failure criteria where more than one stress components have been used to evaluate the different failure modes. These criteria were originally developed for unidirectional polymeric composites, and hence, applications to other type of laminates and non-polymeric composites have significant approximations. Usually Hashin criteria are implemented within two-dimensional classical lamination approach for point stress calculations with ply discounting as the material degradation model. Failure indices for Hashin criteria are related to fibre and matrix failures and involve four failure modes. The criteria are extended to three-dimensional problems where the maximum stress criteria are used for transverse normal stress component.
	The failure modes included in Hashin’s criteria are as follows.

# Tensile fibre failure for σ11 ≥ 0
# Compressive fibre failure for σ11 < 0
# Tensile matrix failure for σ22 + σ33 > 0
# Compressive matrix failure for σ22 + σ33 < 0
# Interlaminar tensile failure for σ33 > 0
# Interlaminar compression failure for σ33 < 0

where, σij denote the stress components and the tensile and compressive allowable strengths for lamina are denoted by subscripts T and C, respectively. XT, YT, ZT denotes the allowable tensile strengths in three respective material directions. Similarly, XC, YC, ZC denotes the allowable compressive strengths in three respective material directions. Further, S12, S13 and S23 denote allowable shear strengths in the respective principal material directions.

Endeavors have been made to incorporate MMF with multiple progressive damage models and fatigue models for strength and life prediction of composite structures subjected to static or dynamic loadings.

== See also ==
* [[Composite material]]
* [[Strength of materials]]
* [[Material failure theory]]
* [[Tsai-Wu failure criterion]]
* [[Christensen failure criterion]]

== References ==
<references>

<ref name=SKHa>
  Ha, S.K., Jin, K.K. and Huang, Y. (2008). Micro-Mechanics of Failure (MMF) for Continuous Fiber Reinforced Composites, ''Journal of Composite Materials'', '''42'''(18): 1873–1895.
</ref>

<ref name=SWTsai>
  Tsai, S.W. and Wu, E.M. (1971). A General Theory of Strength for Anisotropic Materials, ''Journal of Composite Materials'', '''5'''(1): 58–80.
</ref>

<ref name=Hashin1>
  Hashin, Z. and Rotem, A. (1973). A Fatigue Failure Criterion for Fiber Reinforced Materials, ''Journal of Composite Materials'', '''7'''(4): 448–464.
</ref>

<ref name=Hashin2>
  Hashin, Z. (1980). Failure Criteria for Unidirectional Fiber Composites, ''Journal of Applied Mechanics'', '''47'''(2): 329–334.
</ref>

<ref name=ZXia>
  Xia, Z., Zhang, Y. and Ellyin, F. (2003). A Unified Periodical Boundary Conditions for Representative Volume Elements of Composites and Applications, ''International Journal of Solids and Structures'', '''40'''(8): 1907–1921.
</ref>

<ref name=KKJin>
  Jin, K.K., Huang, Y., Lee, Y.H. and Ha, S.K. (2008). Distribution of Micro Stresses and Interfacial Tractions in Unidirectional Composites, ''Journal of Composite Materials'', '''42'''(18): 1825–1849.
</ref>

<ref name=Christen>
  Christensen, R.M. (2007). A Comprehensive Theory of Yielding and Failure for Isotropic Materials, ''Journal of Engineering Materials and Technology'', '''129'''(2): 173–181.
</ref>

<ref name=Camanho>
  Camanho, P.P. and Dávila, C.G. (2002). Mixed-Mode Decohesion Finite Elements for the Simulation of Delamination in Composite Materials, NASA/TM-2002-211737: 1–37.
</ref>
</references>

{{DEFAULTSORT:Micro-Mechanics Of Failure}}
[[Category:Mechanics]]
[[Category:Solid mechanics]]
[[Category:Mechanical failure]]