Abstract:The crack phase field model has been well established and validated for a variety of complex crack propagation patterns within a homogeneous medium under either tensile or shear loading. However, relatively less attention has been paid to crack propagation under combined tensile and shear loading or crack propagation within composite materials made of two constituents with very different elastic moduli. In this work, we compare crack propagation under such circumstances modelled by two representative formulations, anisotropic and hybrid formulations, which have distinct stiffness degradation schemes upon crack propagation. We demonstrate that the hybrid formulation is more adequate for modeling crack propagation problems under combined loading because the residual stiffness of the damaged zone in the anisotropic formulation may lead to spurious crack growth and altered load-displacement response.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the propagation of cracks under combined shear and tensile loads in composite materials, especially in composite materials where the two constituent materials have very different elastic moduli. Traditionally, the phase - field model has been well - established for simulating complex crack propagation patterns in homogeneous media, but there are relatively few studies on crack propagation under combined shear and tensile loads or in composite materials.
Specifically, the paper compares the performance of two representative formulas (the anisotropic formula and the hybrid formula) in simulating these situations. These two formulas have different stiffness degradation schemes during crack propagation. The research shows that the hybrid formula is more suitable for simulating crack propagation problems under combined loads, because the residual stiffness in the damaged area in the anisotropic formula may lead to false crack propagation and a changed load - displacement response.
In summary, the paper mainly focuses on how to more accurately simulate and predict the crack propagation behavior under combined shear and tensile loads, especially for composite materials with large differences in elastic moduli. This helps to improve the performance prediction and safety assessment of such materials in practical engineering applications.
The key formulas expressed in Chinese are as follows:
1. Definition of the phase - field parameter \(d(\mathbf{x})\):
\[
d(\mathbf{x}) =
\begin{cases}
1, & \text{on } \Gamma \\
0, & \text{otherwise}
\end{cases}
\]
2. The differential equation satisfied by the phase - field parameter \(d(\mathbf{x})\):
\[
\begin{cases}
-l^2 \nabla^2 d = 0, & \text{in } \Omega \\
d(\mathbf{x}) = 1, & \text{on } \Gamma \\
\nabla d(\mathbf{x}) \cdot \mathbf{n} = 0, & \text{on } \partial \Omega
\end{cases}
\]
3. Definition of the strain energy \(E(\mathbf{u}, d)\):
\[
E(\mathbf{u}, d) = \int_{\Omega} \psi(\boldsymbol{\varepsilon}(\mathbf{u}), d) \, dV
\]
4. Strain energy in the isotropic formula:
\[
\psi(\boldsymbol{\varepsilon}(\mathbf{u}), d) = (1 - d)^2 \psi_0(\boldsymbol{\varepsilon})
\]
5. Strain energy in the anisotropic formula:
\[
\psi(\boldsymbol{\varepsilon}(\mathbf{u}), d) = (1 - d)^2 \psi_0^+(\boldsymbol{\varepsilon})+\psi_0^-(\boldsymbol{\varepsilon})
\]
6. Strain energy in the hybrid formula:
\[
\psi(\boldsymbol{\varepsilon}(\mathbf{u}), d) = (1 - d)^2 \psi_0^+(\boldsymbol{\varepsilon})+\psi_0^-(\boldsymbol{\varepsilon})
\]
Through these formulas, the paper shows the performance differences of different formulas in simulating crack propagation and emphasizes the advantages of the hybrid formula.