Skip to topic | Skip to bottom

Start of topic | Skip to actions

Fracture of a thin copper plate from water hammer

The simulations on this page try to reproduce experimental results provided by V.S. Deshpande, University Cambridge. A strong pressure wave in water is created by firing a steel projectile on a piston at the end of the tube. The wave propagates through the tube and impinges onto a thin copper plate sealing the other end. Depending on the initial projectile velocity the pressure induces plastic deformation or different rupture patterns.

The experimental facility is described in debth in V.S. Deshpande, A. Heaver, N.A. Fleck, An underwater shock simulator. Royal Society of London Proceedings Series A, 462(2067):1021-1041, 2006.

Basic computational setup


  • Pressure wave generated by solving equation of motion for piston during entire fluid-structure simulation. Piston mass: 238g. The initial velocity v0 is varied to create different loading conditions with initial pressure p0.
  • Full water shock tube of 1.3m length, 64mm diameter considered in fluid. 1d simulations used to provide proper flow field for coupled 3d simulations.
  • Modeling of water with stiffened gas equation of state with \gamma=7.415, p_\infty=296.2\,\rm MPa
  • Multi-dimensional 2nd order upwind finite volume scheme, negative pressures from cavitation avoided by energy correction
  • AMR base level: 350x20x20, 2 additional levels, refinement factor 2,2.
  • Approx. 1.2M cells used in fluid on average instead of 9M (uniform)
  • An unsigned distance function is used to allow the simulation of the rupture event.

Solid with SFC

  • Copper plate of 0.25mm, J2 plasticity model with hardening, rate sensitivity, and thermal softening, cohesive interface model
  • Solid mesh: 4675 nodes, 8896 elements
  • 12 nodes 3.4 GHz Intel Xeon dual processor, Gigabit ethernet network, ca. 800h CPU to t=1ms
  • Source codes: [fluid codes] [solid codes]

p0=34MPa, v0=-22.94m/s

Pressure traces Plate deformation
Pressure traces p0=34MPa Plate deformation p0=34MPa

p0=64MPa, v0=-43.18m/s

Pressure traces Plate deformation
Pressure traces p0=64MPa Plate deformation p0=64MPa

p0=88MPa, v0=-59.38m/s

Pressure traces Plate deformation
Pressure traces p0=88MPa Plate deformation p0=88MPa

p0=130MPa, v0=-87.72m/s

Pressure traces Plate deformation
Pressure traces p0=130MPa Plate deformation p0=130MPa

p0=173MPa, v0=-116.73m/s

Pressure traces Plate deformation
Pressure traces p0=173MPa Plate deformation p0=1730MPa

-- RalfDeiterding - 21 Apr 2008

I Attachment sort Size Date Who Comment
Press34MPa.mpg 5267.8 K 11 Mar 2008 - 20:08 RalfDeiterding Plate deformation, p0=34MPa
Press88MPa.mpg 5437.1 K 11 Mar 2008 - 20:09 RalfDeiterding Pressure loading, p0=88MPa
VFrac88MPa.mpg 6085.2 K 11 Mar 2008 - 20:10 RalfDeiterding Volume fraction, p0=88MPa
Press173MPa.mpg 3063.6 K 11 Mar 2008 - 20:11 RalfDeiterding Pressure loading, p0=173MPa
VFrac173MPa.mpg 3182.0 K 11 Mar 2008 - 20:12 RalfDeiterding Volume Fraction, p0=173MPa
Press64MPa.mpg 6173.6 K 18 Mar 2008 - 18:38 RalfDeiterding Pressure loading, p0=64MPa
VFrac64MPa.mpg 7296.0 K 18 Mar 2008 - 18:39 RalfDeiterding Volume Fraction, p0=64MPa
Press130MPa.mpg 3934.3 K 18 Mar 2008 - 18:39 RalfDeiterding Pressure loading, p0=130MPa
VFrac130MPa.mpg 4110.7 K 18 Mar 2008 - 18:40 RalfDeiterding Volume Fraction, p0=130MPa

You are here: CoupledSimulations > PlateWaterhammerFrac

to top

Copyright © 1997-2018 California Institute of Technology.