Figure 1 : Load‑displacement response (phase‑field vs. LEFM). Figure 2 : Phase‑field contour at (F = 0.9F_c) (crack tip radius ≈ 3(\ell)). A DCB specimen (length 0.2 m, thickness 0.01 m) is subjected to a symmetric opening displacement. The energy release rate calculated from the phase‑field solution
Given uⁿ, φⁿ: 1. Update history field Hⁿ⁺¹ ← max(Hⁿ, ψ⁺(ε(uⁿ))) 2. Solve displacement problem → uⁿ⁺¹ (with φⁿ fixed) 3. Solve phase‑field problem → φⁿ⁺¹ (with uⁿ⁺¹ fixed) 4. Check convergence: ‖uⁿ⁺¹‑uⁿ‖ + ‖φⁿ⁺¹‑φⁿ‖ < ε_tol 5. If not converged → repeat steps 2‑4 The linearised systems are assembled using (e.g., via the Sacado package) to obtain consistent tangent operators. 3.4. Load Control & Arc‑Length For softening problems, displacement control can cause snap‑back. We implement an arc‑length (Riks) method that controls the total work increment: Working Model 2d Crack-
All source files are provided in the supplementary material (GitHub repository github.com/YourGroup/2DPhaseFieldCrack ). 4.1. Benchmark 1 – Single‑Edge Notched Tension (SENT) Geometry : rectangular plate (L=1.0) m, (H=0.5) m, notch length (a_0=0.2) m. Material : (E=30) GPa, (\nu=0.2), (G_c=2.7) kJ/m(^2). Parameters : (\ell = 2.5,h_\min) (where (h_\min) is the smallest element size after refinement). Figure 1 : Load‑displacement response (phase‑field vs
[ \psi^+(\boldsymbol\varepsilon) ;\rightarrow; H(\mathbfx) . \tag4 ] 3.1. Finite‑Element Discretisation Both fields are approximated using quadratic Lagrange shape functions on an unstructured triangular mesh: A DCB specimen (length 0
[ \mathbfu^h(\mathbfx) = \sum_i=1^N_n \mathbfN_i(\mathbfx) , \mathbfu i, \qquad \phi^h(\mathbfx) = \sum i=1^N_n N_i(\mathbfx) , \phi_i, \tag5 ]
[ \Delta W = \int_\Gamma_N \mathbft\cdot \Delta\mathbfu,\mathrmdS . \tag7 ]