A blended approach for evolving phase fields using peridynamics: Cyclic loading in quasi-brittle fracture

This paper presents a novel thermodynamically consistent, mesh-free phase field theory that integrates peridynamics with plasticity and damage mechanics to accurately predict cyclic loading behavior and fracture patterns in quasi-brittle materials.

The Gist

Imagine you are watching a piece of concrete crack under pressure. In the real world, this isn't just a clean snap like breaking a dry twig. It's messy. The material squishes, stretches, heats up, and leaves behind a "bruise" (damage) that doesn't heal when you let go. It remembers the pain.

For decades, scientists have struggled to write a computer program that can predict exactly how and where this happens, especially when the material is being pushed and pulled repeatedly (like a bridge swaying in the wind or a car hitting a pothole over and over).

This paper introduces a new, clever way to do this. Think of it as a hybrid recipe that mixes two different cooking styles to bake the perfect "fracture cake."

The Two Ingredients: The "Rubber Sheet" and the "Memory Glass"

The authors combine two powerful ideas:

1. Peridynamics (The Rubber Sheet): Imagine the material isn't a solid block, but a giant, invisible rubber sheet made of billions of tiny dots connected by springs. In traditional physics, a dot only talks to its immediate neighbors. In this "Peridynamic" view, a dot can talk to neighbors a little further away. If a spring breaks, the dot stops talking to that neighbor, and a crack forms naturally. No complex rules needed; the crack just appears where the springs snap.
2. Phase Fields (The Memory Glass): Now, imagine the material has a "memory." If you stretch a piece of clay, it stays stretched a little bit even after you let go. This is called plasticity. The "Phase Field" is like a variable that tracks how "tired" or "damaged" a specific spot is. It's a number between 0 (completely broken) and 1 (brand new).

The Big Innovation: One Law to Rule Them All

Usually, scientists have to write two separate sets of rules: one for how the material moves (Newton's laws) and another for how the damage grows. It's like telling a driver how to steer the car and separately telling the engine how to break down.

This paper says: "Let's just use one rule."

They created a "Blended Approach" where the damage (the phase field) is part of the material's memory, not a separate instruction.

• The Analogy: Imagine a crowd of people holding hands (the bonds). If someone gets pushed too hard, their grip weakens (damage). If they get pushed back and forth (cyclic loading), their hands get sweaty and slippery (plasticity), and they might never hold on as tight again.
• In this new model, the "grip strength" of the hand automatically adjusts based on how hard and how often they were pushed. You don't need a separate rule to say "hands get slippery." The math of the push naturally makes the hands slippery.

How It Handles the "Messy" Stuff (Quasi-Brittle Materials)

Materials like concrete, brick, and stone are "quasi-brittle." They are tough but not unbreakable.

• Tension (Pulling): When you pull concrete, it cracks and loses strength.
• Compression (Squeezing): When you squeeze it, it actually gets stronger or stays stiff.

The authors' model uses two different "memory glasses" (phase fields):

1. One for Pulling: Tracks when the material is stretched too far and starts to crack.
2. One for Squeezing: Tracks when the material is crushed.

Crucially, the model understands that if you pull a piece of concrete, it gets damaged. If you then squeeze it, the crack might close up, but the material is still "bruised" and won't be as strong as before. The model captures this hysteresis (the lag between loading and unloading) perfectly, just like a real rubber band that gets loose after being stretched too many times.

The "Magic" of the Simulation

When the authors ran their computer simulations, they didn't just guess where cracks would go. They let the physics decide.

• No Mesh: Traditional simulations use a grid (like graph paper) to calculate forces. If the crack doesn't follow the grid lines, the math gets messy. This method is mesh-free. It's like having a cloud of dust where the particles move freely. The crack can go in any direction, curving and twisting naturally.
• Size Matters: They tested this on small beams and huge beams. Real concrete behaves differently depending on its size (a small beam is stronger per inch than a giant one). This model predicted that "size effect" automatically, without the scientists having to manually tell the computer, "Hey, make the big one weaker."

The Bottom Line

This paper presents a new "universal translator" for fracture mechanics. It takes the messy, real-world behavior of materials like concrete—where they stretch, squish, remember their damage, and crack in unpredictable ways—and translates it into a single, elegant set of mathematical rules based on Newton's laws.

In simple terms: They figured out how to make a computer simulation where the material "learns" from its own damage. If you pull it, it gets weak. If you push it, it gets stiff. If you do it over and over, it gets tired. And when it finally breaks, the computer predicts the crack path exactly as it would happen in real life, without needing a human to draw the lines.

This is a huge step forward for designing safer bridges, buildings, and aircraft that can withstand the repeated stresses of the real world.

Technical Summary

Here is a detailed technical summary of the paper "A blended approach for evolving phase fields using peridynamics: Cyclic loading in quasi-brittle fracture" by Hayden Bromley and Robert Lipton.

1. Problem Statement

The identification of a suitable field theory to solve the fracture initial boundary value problem (IBVP) for quasi-brittle materials (e.g., concrete, masonry) under cyclic loading remains a significant challenge in continuum mechanics. Traditional approaches face specific limitations:

• Classical Phase Field Methods: Typically couple a time-evolved phase field (governed by a Ginzburg-Landau or Ambrosio-Tortorelli equation) with the displacement field (governed by Newton's second law). This often requires separate evolution equations and variational constraints, which can complicate the thermodynamic consistency and dynamic coupling.
• Peridynamics: While non-local and capable of handling discontinuities, standard bond-based peridynamics often struggles to independently calibrate elastic moduli, strength, and fracture toughness without introducing ad-hoc parameters.
• Cyclic Loading: Existing models often fail to accurately capture the complex hysteresis, plastic dissipation, and damage accumulation associated with cyclic loading in quasi-brittle materials, particularly the distinction between tensile and compressive damage states.

The authors aim to develop a unified, mathematically well-posed framework that predicts damage and fracture in quasi-brittle materials incorporating irreversible (plastic) deformation and damage-dependent elastic softening, driven exclusively by Newton's second law.

2. Methodology

The paper proposes a "Blended Approach" that integrates a non-local constitutive law (Peridynamics) with a history-dependent two-point phase field.

A. Theoretical Framework

• Governing Equation: The displacement field $u(t, x)$ evolves solely according to Newton's second law:
  $$\rho \ddot{u}(t, x) + L_\epsilon[u](t, x) = b(t, x)$$
  where $L_\epsilon$ is a non-local integral operator representing internal forces. There is no separate evolution equation for the phase field; the phase field is a constitutive variable determined by the displacement history.
• Two-Point Phase Fields: The model introduces two independent phase fields, $\gamma_+$ (tension) and $\gamma_-$ (compression), defined on pairs of points $(x, y)$.
  • These fields represent the degradation of elastic stiffness ($0 \le \gamma \le 1$).
  • They are history-dependent, tracking the maximum strain experienced by a bond up to time $t$.
• Constitutive Law with Memory:
  • Strain Decomposition: Total strain is additively decomposed into elastic and irreversible (plastic) parts.
  • Strength Envelope: The stress-strain behavior is defined by a strength envelope (yielding, hardening, and failure) and a family of unloading laws.
  • Unloading Ratio ($\beta$): A key parameter introduced to calibrate the ratio between plastic strain associated with elastic softening and total irreversible strain. This allows the model to distinguish between brittle fracture ($\beta \approx 0$) and quasi-brittle fracture ($\beta > 0$).
  • Tension vs. Compression: Bonds can fail in tension (zero stiffness) but continue to resist compression elastically, mimicking the behavior of cracks closing under compression.

B. Mathematical Properties

• Well-Posedness: The authors prove the existence and uniqueness of the solution to the IBVP using the Banach fixed-point theorem in a suitable Banach space.
• Thermodynamic Consistency: The model satisfies the energy balance law. The rate of energy dissipation ($\dot{D}$) is proven to be non-negative, consistent with the laws of thermodynamics.
• Fracture Energy Recovery: Using geometric measure theory, the authors demonstrate that the energy required to form a flat crack converges to the Griffith fracture energy ($G_c \times \text{Area}$) as the non-locality length scale vanishes, independent of the specific length scale used in the simulation.

C. Numerical Implementation

• Mesh-Free Method: The continuum equations are discretized using a mesh-free approach based on quadrature (nodes with no geometric connections).
• Calibration: Material parameters (Young's modulus, tensile strength, critical energy release rate) are calibrated directly from experimental data. The non-locality length scale ($\epsilon$) is treated as a numerical resolution parameter, not a material property.
• Dynamic Relaxation: Quasi-static results are obtained using a dynamic relaxation method with damping.

3. Key Contributions

1. Unified Evolution: A formulation where fracture evolution is driven entirely by Newton's second law, eliminating the need for a separate phase field evolution equation (e.g., Ginzburg-Landau).
2. Two-Point History Dependence: The introduction of two independent phase fields ($\gamma_+, \gamma_-$) that capture distinct damage states for tension and compression, enabling the modeling of cyclic loading and hysteresis.
3. Thermodynamic Rigor: A rigorous proof that the model satisfies energy balance and positive energy dissipation, derived directly from the equations of motion.
4. Parameter Independence: The non-locality length scale ($\epsilon$) is decoupled from material properties (strength, fracture energy), allowing it to be used purely for numerical resolution.
5. Size Effect Prediction: The model naturally captures the deterministic size effect in quasi-brittle materials (where nominal strength decreases with increasing structural size) without ad-hoc size-dependent rules.

4. Results

The authors validated the model against experimental data for high-strength concrete and other quasi-brittle materials across several benchmark problems:

• Mesh Convergence: Demonstrated convergence of the numerical solution as the grid spacing decreases and the horizon-to-grid ratio ($m$) increases.
• Cyclic Loading (Mode I): Successfully reproduced the hysteresis loops (Force vs. Crack Mouth Opening Displacement) observed in three-point bending tests, matching experimental data from Chen and Liu (2023). The model captured plastic strain localization and the process zone ahead of the crack.
• Mixed-Mode Fracture: Simulated mixed-mode fracture in notched beams (Jenq and Shah, 1988) and L-shaped panels (Winkler, 2001). The predicted crack paths and force-displacement curves showed excellent agreement with experimental observations, including the transition from Mode I to mixed-mode sliding.
• Size Effect: The model accurately predicted the reduction in nominal strength for concrete beams of varying sizes (García-Álvarez et al., 2012), aligning with the Bažant size-effect law. The results matched experimental data and other peridynamic models (Hobbs et al., 2022).

5. Significance

This work represents a significant advancement in computational fracture mechanics for quasi-brittle materials:

• Physical Fidelity: By deriving fracture and damage directly from Newton's laws and thermodynamic principles, the model avoids the phenomenological assumptions often found in traditional phase field methods.
• Cyclic Loading Capability: It provides a robust framework for simulating fatigue and cyclic loading, which is critical for infrastructure safety but difficult to model with standard brittle fracture theories.
• Computational Efficiency: The mesh-free nature and the elimination of a separate phase field PDE simplify the numerical implementation while maintaining high accuracy.
• Predictive Power: The ability to capture size effects and complex crack paths (including mixed-mode) using a single set of calibrated material constants suggests the model is ready for predictive engineering applications in concrete and geomaterials.

In summary, the paper presents a mathematically rigorous, thermodynamically consistent, and computationally efficient "blended" approach that successfully bridges the gap between peridynamics and phase field methods to solve complex fracture problems in quasi-brittle materials under cyclic loading.