Damage Mechanics Challenge#

In 2022, we took part in the Damage Mechanics Challenge, where the assignment was to most accurately predict a complex, mixed-mode fracture in an anisotropic, 3D-printed rock material. To this end, experimental data from various tests were provided [1,2], based on which a computational approach was to be calibrated. We chose the well-established phase-field approach for brittle fracture and focused on developing a calibration procedure to properly set the material parameters. The main idea was to recast the parameter calibration into an optimization problem, with the objective of minimizing the deviation between the numerical and experimental force-displacement response. Our detailed approach and findings can be found in the paper [3]. Below we provide AR models to visualize the central 3D aspect of the obtained crack surfaces for various three point bending tests. These include the HC (centered notch), HB (off-centered notch), HA (far off-centered notch) and H45 (angled initial notch) test, which we used to calibrate the parameters of our model, as well as the DMC test, which is the actual challenge prediction. All the experimental data shown is postprocessed data taken from [1,2].

Note that only the elements with \(d \geq 0.25\) are shown in the outline of the respective three point bending test geometry.

DMC Test#

colormap_coolwarm
\(1\)
\(d\)
\(0\)

colormap_turbo
\(1.90\)mm
\(\Delta s\)
\(0\)

HC Test#

colormap_coolwarm
\(1\)
\(d\)
\(0\)

colormap_turbo
\(1.11\)mm
\(\Delta s\)
\(0\)

HB Test#

colormap_coolwarm
\(1\)
\(d\)
\(0\)

colormap_turbo
\(1.80\)mm
\(\Delta s\)
\(0\)

HA Test#

colormap_coolwarm
\(1\)
\(d\)
\(0\)

colormap_turbo
\(4.68\)mm
\(\Delta s\)
\(0\)

H45 Test#

colormap_coolwarm
\(1\)
\(d\)
\(0\)

colormap_turbo
\(1.18\)mm
\(\Delta s\)
\(0\)

Material Properties#

To visualize the material constants we calibrated for the assumption of orthotropic material, we employ the strategy and implementation presented in [4]. Below, we show the Young’s modulus \(E\) as well as the average and maximum of the Poisson’s ratio \(\nu\) along a director \(\boldsymbol{d}\).

colormap_turbo
\(2220\) MPa
\(E\)
\(280\) MPa
colormap_turbo
\(0.53\)
\(\text{mean}(\nu)\)
\(0.24\)
colormap_turbo
\(1.46\)
\(\max(\nu)\)
\(0.33\)

References#

[1] Jiang, L., Pyrak-Nolte, L., Yoon, H., Bobet, A. and Morris, J. (2021): Calibration data set for damage mechanics challengeon brittle-ductile material, doi.org/10.4231/QF39-Q924
[2] Jiang, L., Pyrak-Nolte, L., Yoon, H., Bobet, A. and Morris, J. (2021): Digital image, x-ray ct , xrd and ultrasonic data sets for damage mechanics challenge on brittle-ductile material, doi.org/10.4231/2E8M-W085
[3] J. Heinzmann, P. Carrara, C. Luo, M. Manav, A. Mishra, S. Nagaraja, H. Oudich, F. Vicentini, D. De Lorenzis (2024): Calibration of a Phase-field Model of Brittle Fracture with Experimental Data on Anisotropic Rock (submitted)
[4] J. Nordmann, M. Aßmus and H. Altenbach (2018): Visualising elastic anisotropy: theoretical background and computational implementation, Continuum Mechanics and Thermodynamics 30: 689–708, doi.org/10.1007/s00161-018-0635-9