Geometric Crack Modeling Module

• Geometric crack modeling modules are computational frameworks that simulate crack propagation using explicit geometric representations and lattice-based algorithms.
• They integrate elastic stress solvers, threshold criteria, and wavelet-based roughness analyses to quantify crack front morphology and regime transitions.
• The module supports both displacement-control and gradient loading to calibrate material properties and match experimentally observed scaling exponents.

A geometric crack modeling module is a computational framework for simulating, analyzing, and predicting the growth and morphology of cracks in solids through explicit geometric or lattice-based representations. These modules encode crack topology, discretization of material interactions, and algorithms for crack advance, branching, and statistical characterization. The goal is to capture both the physical mechanisms (elastic interactions, threshold disorder) and the geometric/statistical features (front roughness, scaling exponents) of interfacial crack propagation.

1. Model Geometry and Physical Setup

The essential geometric framework, as presented in "A model for stable interfacial crack growth" (Gjerden et al., 2012), consists of a constrained configuration where two blocks—one perfectly rigid, one with finite Young's modulus $E$ and Poisson ratio $\nu$—are separated by a square lattice of $N=L^2$ linear elastic fibers (discrete springs), each spanning a patch $a \times a$. The crack propagates along the interfacial plane between the blocks, with the lattice discretization indexed in two dimensions: $i \equiv (x_i, y_i)$.

Boundary conditions are periodic in the direction parallel to the crack front ($x$), simulating an infinite system. Depending on the system's stiffness, bi-periodic boundary conditions (for $x$ and $y$) are employed for numerical stability in conjugate-gradient solvers for soft regimes.

Loading is implemented either by imposing a global displacement $D$ (displacement-control), where intact fiber $i$ carries force $f_i = -k_i(u_i-D)$, or by assigning fiber breaking thresholds $t_i = t_0 + g\,y_i$ (threshold-gradient loading), simulating an external gradient of driving force.

2. Governing Equations and Failure Criteria

The crack modeling module is constructed from the following field equations and key physical assumptions:

• Stress-strain relations: Each fiber is linear elastic, with force law

$f_i = -k_i(u_i-D)$

where $k_i$ is the local stiffness ($k_i=1$ for intact, $k_i=0$ for broken fibers). Displacement $u_i$ is not local but coupled: the elastic block transmits force via a discretized Green's function,

$$u_i = \sum_{j\,\text{intact}} G_{ij} f_j,\qquad G_{ij} = \frac{1-\nu^2}{\pi E a^2} \int_{-a/2}^{a/2}\int_{-a/2}^{a/2} \frac{1}{|\vec{r}_i - \vec{r}_j'|}\,dx'\,dy'$$

which, in practice, is implemented via Love's formula for a square patch.

• Fracture criterion: For each fiber, local tensile strain is

$s_i = f_i/k_i + D - u_i$

which is compared to the breaking threshold $t_i$. In computation, the fiber with largest $R_i = f_i / t_i$ is chosen for irreversible failure (extremal-dynamics).

• Crack advance: There is no explicit calculation of a stress-intensity factor; crack evolution occurs by local failures, triggered when $(u_i-D)$ meets the threshold $t_i$ at a site.

3. Geometric and Statistical Front Analysis

To extract and quantify crack morphology, the module tracks a front height function $h(x)$ by computing the upper envelope of the largest connected cluster of broken fibers. This mapping removes overhangs using a solid-on-solid (SOS) criterion ($h(x+1)-h(x) \geq -1$), yielding a well-defined single-valued front.

The module then measures self-affine scaling properties: $$\langle [h(x+\ell)-h(x)]^2 \rangle \sim \ell^{2\zeta}$$ Alternatively, wavelet-based roughness analysis computes the mean modulus of wavelet coefficients $W(a) \sim a^{\zeta+1/2}$ over scale $a$.

Typical roughness exponents—$\zeta_1 \sim 0.45$ (stiff regime), $\zeta_2 \sim 0.30$ (soft regime)—are observed, matching but below the corresponding experimental values ($\zeta_1 \sim 0.60$, $\zeta_2 \sim 0.35$).

4. Scaling Mechanism and Physical Regimes

The geometric crack modeling module reveals two distinct self-affine regimes, each arising from different physical mechanisms:

• Small-scale regime ($\ell \ll \ell_c$): Dominated by fracture coalescence. Microcrack islands nucleate ahead of the front and merge (coalescence model). Roughness exponent $\zeta \approx 0.6$ experimentally, $\zeta_1 \approx 0.45$ in the model.
• Large-scale regime ($\ell \gg \ell_c$): The crack front behaves as a fluctuating elastic line driven through a random medium (line model). $\zeta_2 \approx 0.30$ (model); $\zeta_2 \approx 0.35$ (experiment).

The crossover length $\ell_c$ between regimes increases with block stiffness $E$ and decreases with gradient $g$, approximately following

$\ell_c \sim E^\alpha / g^\beta$

where $\alpha,\,\beta > 0$ are non-universal parameters dependent on system details.

5. Computational Architecture and Algorithms

Implementation is structured around efficient data handling and algorithmics:

Core data structures:

• $k[i]$: status array (intact = 1, broken = 0)
• $t[i]$: threshold array
• $f[i]$: current fiber forces
• $u[i]$: local displacements
• $h[x]$: front height extraction

Algorithm pseudocode:

initialize: for each i: k[i] ← 1 t[i] ← t₀ + g·y_i + random_noise set global D ← 0 choose loading scheme loop until all fibers broken or end condition: 1. solve (I + K·G)·f = K·D [CG + FFT acceleration] 2. compute R[i] = f[i] / t[i] for all intact i 3. find i* = argmax R[i] 4. k[i*] ← 0 // break 5. if treadmill used: shift fields downward, reinitialize top row 6. extract h(x) from envelope of largest broken cluster 7. accumulate statistics on h(x)

Stress solver: Fast conjugate-gradient (CG) with FFT-based acceleration. Stress transfer is diagonal in Fourier space due to the Green's function structure. Preconditioning is required at high stiffness; double precision is recommended to control round-off.

Crack-advance logic: The extremal-dynamics rule—breaking the fiber of largest $R[i]$ at each step—naturally advances the crack front without explicit computation of $K_I$.

6. Regime Calibration and Practical Integration

For robust quantitative modeling, system parameters are selected according to physical and computational considerations:

• System size: $L \geq 256$ recommended to reduce finite-size effects, especially in extracting roughness exponents.
• Elastic modulus $E$ and Poisson ratio $\nu$: Wide sampling across $E/E_\text{ref}=10^{-3}$ to $10^{3}$, typical $\nu=0.25$–$0.35$.
• Threshold gradient $g$: Chosen so that $W \approx L/10$–$L/5$.
• CG tolerance: $10^{-6}$–$10^{-8}$ per solve.
• Computational complexity: Each CG iteration is $O(N\log N)$; total cost $\sim N_\text{broken}\,N\log N$. Preconditioning required for numerical stability in high-$E$ regimes.

Module extensions include:

• Calibration of $g$ and threshold noise against experimental data on fracture toughness.
• Substitution of non-uniform $k_i$ for modeling material heterogeneity.
• Introduction of anisotropic $G_{ij}$ for layered media.
• Extraction and fitting of $\ell_c(E,g)$ via wavelet or height function analysis.

7. Scientific Impact and Comparison to Experiments

The geometric crack modeling module systematically unifies mechanisms for stable interfacial crack growth, matching experimentally observed dual-regime self-affine scaling. Its algorithms are suitable for direct implementation, numerical parameter scans, and quantitative comparison to physical measurements. Notably, the module resolves the regime crossover and provides a computational laboratory for testing statistical properties, scaling laws, and the influence of elasticity and disorder on crack front morphology (Gjerden et al., 2012).

This framework's modularity and explicit geometrization make it adaptable to broader classes of interfacial fracture problems, including layered composites, anisotropic systems, and real materials with calibrated disorder. Its statistical characterization tools (SOS mapping, wavelet analysis) are directly transferable to experimental imaging data, facilitating rigorous theory–experiment comparison.