Home

Dynamics of Crack Propagation

 V.A.Kalatsky

This page contains selected animations, papers, and presentations of my work on dynamics of crack propagation at
Materials Science Division (MSD),  Argonne National Laboratory.

Goal

        A continuum theory capable of description, in the same unified framework, of the whole phenomenology of the fractures: crack initiation, quasi-stationary propagation, instability, branching.

Articles and Presentations

Preprint of the paper in PDF or PS format (published in PRL 85, 118 [2000]).
A short introduction presented at the review of the University of Chicago, in PDF or PS format.
A short introduction presented at the American Physical Society March meeting in Minneapolis, MN, in PDF or PS format.
A longer colloquium slide presention by V.A.Kalatsky. Presented at the University of Chicago, Jan. 19, 2000, Northwestern University, Feb. 18, 2000, and Northern Illinois University, Feb. 25, 2000, in PDF or PS format.

Animations

        The animations demonstrate results of simulations of the continuum field model of crack propagation in brittle amorphous solids.   The core of the model is the equations of elastic displacements coupled to the order parameter equation which accounts for the dynamics of material defects.   The model captures all important phenomenology of crack propagation: crack initiation, propagation, dynamic fracture instability, sound emission, crack branching and fragmentation.   All data were generated on SGI Origin 2000 (128 R10000 processors) at Mathematics and Computer Science Division of Argonne National Laboratory.   The animations are in three formats: QuickTime movie, SGI movie, and GIF animation. The SGI movies have high quality and VERY LARGE SIZE.

        The following animations are the result of simulations on 4000x800 grid points (2000x400 domain size) with fixed-grips boundary conditions. The one-panel animations show dynamics of the order parameter only. The four-panel animations show dynamics of the order parameter and the stress-tensor components $\sigma_{xx}$ , $\sigma_{yy}$ , and $\sigma_{xy}$ respectively. The last two animations in this series (relative displacement 0.04 and 0.05) are based on data from 2000x400 grid point simulations with extremely large loading.

Relative GIF GIF SGI QuickTime
displacement animation animation movie movie
0.005 delta=0.005 delta=0.005 delta=0.005 delta=0.005
0.01 delta=0.01 delta=0.01 delta=0.01 delta=0.01
0.015 delta=0.015 delta=0.015 delta=0.015 delta=0.015
0.02 delta=0.02 delta=0.02 delta=0.02 delta=0.02
0.025 delta=0.025 delta=0.025 delta=0.025 delta=0.025
0.03 delta=0.03 delta=0.03 delta=0.03 delta=0.03
0.04
delta=0.04 delta=0.04 delta=0.04
0.05
delta=0.05 delta=0.05 delta=0.05

        The next set of animations is the result of simulations on 2400x400 grid points (1200x200 domain size) with some parameters being different from those in the previous set. To prevent retraction of the branches we "froze" the order parameter behind the leading crack tip. The format of the files is QuickTime movie. All movies are accompanied with sound effects. To produced the audio files we used velocity of the crack tip as a frequency modulator. The three panels show the order parameter, hydrostatic pressure ( $\sigma_{xx}$ + $\sigma_{yy}$ ), and shear stress ( $\sigma_{xy}$ ) respectively.

Relative GIF QuickTime
displacement animation movie
0.06 delta=0.06 delta=0.06
0.089 delta=0.089 delta=0.089
0.11 delta=0.11 delta=0.11

This collection was prepared by Valery Kalatsky, January 8, 2000.
 Last updated on October 30, 2003.
Home
  © 2003 Valery A. Kalatsky, Ph.D. All rights reserved.