Working Model 2d Crack- ❲Verified Source❳

The load‑displacement curve obtained with the phase‑field model matches the analytical LEFM prediction for the critical stress intensity factor (K_IC= \sqrtE G_c). The computed (F_c= 4.58) kN is within 2 % of the analytical value. The crack path follows the straight line of the notch, confirming the absence of mesh bias.

[ \eta_e = \int_\Omega_e \ell |\nabla\phi^h|^2 ,\mathrmdV . \tag6 ]

where (N_n) is the number of nodes. Quadratic interpolation is essential to resolve the steep gradients of (\phi) within the diffusive crack zone. A goal‑oriented error estimator based on the phase‑field gradient is used: Working Model 2d Crack-

[ \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 ] [ \eta_e = \int_\Omega_e \ell |\nabla\phi^h|^2 ,\mathrmdV

[ \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:

Elements with (\eta_e > \eta_\texttol) are refined (bisected) and coarsening is applied where (\eta_e < 0.1,\eta_\texttol). This strategy concentrates degrees of freedom only where the crack evolves, keeping the global problem size modest. A monolithic coupling (solving (\mathbfu) and (\phi) simultaneously) is possible but computationally expensive. Instead, we adopt the staggered scheme (Miehe et al., 2010) that is unconditionally stable for quasi‑static loading: A goal‑oriented error estimator based on the phase‑field

The regularisation length (\ell) controls the width of the diffusive crack zone ((\approx 3\ell)). When (\ell\to0), (\Pi) (\Gamma)-converges to the classical Griffith functional. Stationarity of (\Pi) with respect to admissible variations (\delta\mathbfu) and (\delta\phi) yields the coupled Euler‑Lagrange equations :