FeynKrack: A continuum model for quasi-brittle damage through Feynman-Kac killed diffusion

Ved Prakash,Upadhyayula M. M. A. Sai Gopal,Sanhita Das,Ananth Ramaswamy,Debasish Roy

2024-12-01

Abstract:Continuum damage mechanics (CDM) is a popular framework for modelling crack propagation in solids. The CDM uses a damage parameter to quantitatively assess what one loosely calls `material degradation'. While this parameter is sometimes given a physical meaning, the mathematical equations for its evolution are generally not consistent with such physical interpretations. Curiously, degradation in the CDM may be viewed as a change of measures, wherein the damage variable appears as the Radon-Nikodym derivative. We adopt this point of view and use a probabilistic measure-valued description for the random microcracks underlying quasi-brittle damage. We show that the evolution of the underlying density may be described via killed diffusion as in the Feynman-Kac theory. Damage growth is then interpreted as the reduction in this measure over a region, which in turn quantifies the disruption of bonds through a loss of force-transmitting mechanisms between nearby material points. Remarkably, the evolution of damage admits an approximate closed-form solution. This brings forth substantive computational ease, facilitating fast yet accurate simulations of large dimensional problems. By selecting an appropriate killing rate, one accounts for the irreversibility of damage and thus eliminates the need for ad-hoc history-dependent routes typically employed, say, in phase field modelling of damage. Our proposal FeynKrack (a short form for Feynman-Kac crack propagator) is validated and demonstrated for its efficacy through several simulations on quasi-brittle damage. It also offers a promising stochastic route for future explorations of non-equilibrium thermodynamic aspects of damage.

Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the prediction of the complex micro - crack network in quasi - brittle materials (such as concrete, rock, ceramics, etc.) before fracture. Specifically, the paper proposes a new continuous damage mechanics model - FeynKrack. This model describes the development of random micro - cracks inside the material through the "killing diffusion" process in the Feynman - Kac theory, thereby achieving efficient simulation of quasi - brittle damage. The following are the key problems that the paper attempts to solve:

1. **Enhancing the physical meaning and mathematical basis of the model**: Existing quasi - brittle damage models often lack sufficient physical motivation or mathematical basis, and it is difficult to replace repetitive and costly experiments with high precision. The FeynKrack model enhances the physical meaning and mathematical rigor of the model by introducing probability measure values to describe the random distribution of micro - cracks.

2. **Handling the non - local effects of materials**: Traditional continuous damage mechanics models usually assume that damage is local, but in reality, the damage of a material point may be affected by its neighborhood. The FeynKrack model considers the interaction between micro - cracks by introducing non - local effects, and more accurately describes the damage process of the material.

3. **Eliminating artificial history - dependent paths**: Many existing damage models (such as phase - field models) need to introduce artificial history - dependent paths to handle the irreversibility of damage, which increases the complexity and computational cost of the model. The FeynKrack model naturally handles the irreversibility of damage by choosing an appropriate "killing rate", simplifying the implementation of the model.

4. **Providing an efficient numerical simulation method**: The damage evolution equation of the FeynKrack model has a closed - form solution, which greatly reduces the amount of calculation and makes it possible to quickly and accurately simulate large - scale problems.

### Main contributions of the paper

1. **Probability measure value description**: The paper interprets the damage parameter as the Radon - Nikodym derivative of the probability measure, and describes the distribution of random micro - cracks through the probability measure value, providing a new perspective for damage evolution.

2. **Feynman - Kac killing diffusion**: The paper uses the "killing diffusion" process in the Feynman - Kac theory to describe damage evolution. This process not only conforms to the physical mechanism but also has a closed - form solution, improving computational efficiency.

3. **Handling of non - local effects**: The paper introduces non - local effects through an integral operator, considers the interaction between micro - cracks, and enables the model to more accurately describe the damage process of the material.

4. **Numerical verification and application**: The paper verifies the effectiveness of the FeynKrack model through multiple numerical simulations, including crack growth prediction in three - point bending tests, L - shaped plates and square plates with inclined cracks, three - dimensional solids with through - thickness cracks, and Brazilian test specimens with inclined notches.

### Formula summary

- **Radon - Nikodym derivative of probability measure**:
\[
g=\frac{P(d\Lambda_0)}{d\Lambda_0}=\frac{dP(\Lambda_0)}{d\Lambda_0}
\]

- **Feynman - Kac formula**:
\[
q(x,\tau)=\mathbb{E}_P\left[q_t(\xi)\exp\left(-\int_\tau^t G(z,s)\,ds\right)\mid\xi(\tau)=x\right]
\]

- **Transition density function**:
\[
g(|x' - x|)=\frac{1}{(2\pi\beta^2(t)\Delta t)^{n/2}}\exp\left(-\frac{|x' - x|^2}{2\beta^2(t)\Delta t}-G(x,t)\Delta t\right)
\]

- **Non - local strain rate**:
\[
\dot{\epsilon}^* = R(\dot{\epsilon})+\int_{\Lambda_0'}\epsilon\dot{g}(|x' - x|)\,d\Lambda_0'