MPS damage model

This model extends the model based on the Microprestress-Solidification theory, described in previous section and summarized in Tab. 37 for ``MPS'', by tensile cracking.

This extension uses the isotropic damage model with Rankine definition of the equivalent deformation defined as the biggest principal effective stress divided by the elastic modulus. The softening law cannot deal directly with the strain because due to creep strain increases and this would lead to further softening without additional loading.

Two different approaches are implemented. The first one, default, (#1) reduces the stiffness only in the directions of tension (in case the tensile strength is exceeded). A full stiffness is restored in compression and after unloading from tension.

The other approach (#2) is the standard isotropic damage model which reduces the stiffness equally in all directions independently of loading. This approach leads to faster convergence because the secant stiffness can be used instead of the incremental viscoelastic stiffness which must be used in the first approach. The second approach becomes useful when the loading is monotonic or when the benefit of the accelerated computation prevails over the consequences of the reduced/underestimated stiffness in compression.

Proper energy dissipation is guaranteed by the crack-band approach.

The following algorithm is used to compute the stress vector (in each time step and until the iteration criteria are met):

  1. Compute effective stress

    $\displaystyle \sigma_{\mathrm{eff},k+1} = \sigma_{\mathrm{eff},k} + \bar{E}_k \...
... \Delta \varepsilon_k'' - \Delta \varepsilon_{sh,k} - \Delta \varepsilon_{T,k})$ (130)

    where $ \sigma_{\mathrm{eff},k}$ is the effective stress vector in the preceding step, $ \bar{E}_k$ is the incremental stiffness, $ \boldsymbol{D}_\nu$ is a unit stiffness matrix and $ \Delta \varepsilon_k$, $ \Delta \varepsilon_k''$, $ \Delta \varepsilon_{sh,k}$, $ \Delta \varepsilon_{T,k}$ are the increments of the total strain, strain due to creep, shrinkage strain (sum of drying and autogenous shrinkage) and thermal strain, respectively.

  2. Compute principal effective stresses $ \sigma_{\mathrm{eff,1}}$, $ \sigma_{\mathrm{eff,2}}$, $ \sigma_{\mathrm{eff,3}}$

  3. Evaluate equivalent deformation

    $\displaystyle \tilde \varepsilon = \mathrm{max}(\sigma_{\mathrm{eff,1}}, \sigma_{\mathrm{eff,2}}, \sigma_{\mathrm{eff,3}})/E$ (131)

  4. If the stress exceeds the material strength, initialize and fix the fracture parameters $ f_t$ and $ G_F$ (exponential softening assumed).

    $\displaystyle \varepsilon_{c,f} = \frac{G_f}{f_t \: h}$ (132)

    where $ h$ is the characteristic length of the finite element in the direction of the biggest principal stress.

  5. Evaluate corresponding damageif $ \tilde \varepsilon > \varepsilon_0$

    $\displaystyle \omega = 1 - \frac{\varepsilon_0}{\tilde \varepsilon} \exp \left(...
...{\tilde \varepsilon - \varepsilon_0}{\varepsilon_{c,f} - \varepsilon_0} \right)$ (133)

  6. Compute principal nominal stresses

  7. Construct the stress vector in the original configuration

    $\displaystyle \boldsymbol {\sigma} = \boldsymbol{T} \boldsymbol {\sigma}_{princ}$ (134)

    where $ \boldsymbol{T}$ is the stress transformation matrix.

The tensile strength and fracture energy can be defined either as fixed, time-independent values (parameters $ ft$, $ gf$) or as variable, depending on the equivalent hydration time $ t_e$. The latter case is activated by keyword $ timeDepFracturing$; tensile strength and fracture energy are then linked to the current value of the mean compressive strength using (slightly modified) formulae from the fib Model Code 2010.

For $ 0 \: \mathrm{MPa} \leq f_{cm}(t) \leq 20$ MPa the tensile strength is computed from a linear function

$\displaystyle f_{tm}(t) = 0.015 \cdot 12^{2/3} \cdot f_{cm}(t) = 0.07862 f_{cm}(t)$ (135)

which starts at the origin and at $ f_{cm}(t) = 20$ MPa connects to expression

$\displaystyle f_{tm}(t) = 0.3 (f_{cm}(t) - 8 \: \mathrm{MPa})^{2/3}$ (136)

which is valid for $ 20 \: \mathrm{MPa} \leq f_{cm}(t) \leq 58$ MPa. Finally, for $ f_{cm}(t) > 58$ MPa

$\displaystyle f_{tm}(t) = 2.12 \: \ln \left( 1 + 0.1 f_{cm}(t) \right )$ (137)

Following the guidelines from the Model Code 2010 (section 5.1.5.2), the fracture energy $ G_f$ can be estimated from the mean compressive strength at 28 days using an empirical formula

$\displaystyle G_{f,28} = 73 \times \left( f_{cm,28} \right)^{0.18}$ (138)

In this function, the mean compressive strength must be provided in MPa and the resulting fracture energy is in N/m.

The growth in fracture energy is approximately proportional to the tensile strength. The current value of fracture energy is very simply computed as

$\displaystyle G_f(t) = G_{f,28} \times f_{tm}(t) / f_{tm,28}(28)$ (139)

The development of the mean compressive strength in time is adopted entirely from the Model Code 2010 (see section 5.1.9.1)

$\displaystyle f_{cm}(t) = f_{cm,28} \exp{\left( fib\_s \: \left( 1 - \sqrt{28/t_e} \right) \right) }$ (140)

where $ fib\_s$ is a cement-type dependent coefficient (0.38 for cement 32.5 N; 0.25 for 32.5 R and 42.5 N; 0.2 for 42.5 R and 52.5 N and R), $ t_e$ is the equivalent age and $ f_{cm,28}$ is the mean compressive strength at the age of 28 days.

For some concretes the prediction formulae do not provide accurate predictions, therefore it is possible to scale the evolution of tensile strength and fracture energy by providing their 28-day values ft28 and gf28.

The model description and parameters are summarized in Tab. 38.


Table 38: MPS damage-summary.
Description MPS damage model for concrete creep with cracking
  (additionaly all parameters from Table 37)
Record Format MPSDamMat [ ft(rn) #] [ gf(rn) #] [ timeDepFracturing ] [ fib_s(rn) #] [ ft28(rn) #] [ gf28(rn) #] [ damlaw(in) #] [ isotropic ] [ checksnapback(in) #]
Parameters - ft tensile strength (constant in time)
  - gf fracture energy (constant in time)
  - timeDepFracturing string activating fib MC 2010 prediction for tensile strength and fracture energy + their time evolution
  - ft28 manual override for the fib MC 2010 prediction of hydration-dependent tensile strength
  - gf28 manual override for the fib MC 2010 prediction of hydration-dependent fracture energy
  - fib_s cement-type dependent coefficient (compulsory with timeDepFracturing)
  - damageLaw traction-separation law: 0 = exponential softening (default), 1 = linear softening, 6 = damage disabled
  - isotropic string activating same reduction of stiffness after the onset of cracking both in in compression and tension (enables to use secant stiffness matrix for faster convergence)
  - checkSnapBack switch for snap-back checking; 1 = activated (default), 0 = deactiavated
Supported modes 3dMat, PlaneStress, PlaneStrain


Borek Patzak
2018-01-02