This example is intended as a verification test for the Eulerian-Lagrangian coupling in the VTF. An elastic steel panel is hit by a planar Mach-1.21 shock wave in air. The computation is two-dimensional and the beam simulated by a straightforward one-dimensional implicit finite-difference solver that solves the Euler-Bernoulli beam equation
The load q(x,t) is evaluated on-the-fly as the difference of the hydrodynamic pressures on both sides of the beam. The BeamELCCoupledSolver-object is an excellent example for the required interface to utilize a given solid mechanics solvers for coupled fluid-structure interaction simulations.The computations correspond to the simulations in J. Giordano, G. Jourdan, Y. Burtschell, M. Medale, D.E. Zeitoun, L. Houas, Shock wave impacts on deforming panel, an application of fluid-structure interaction. Shock Waves 14 (1-2): 103-110, 2005.
Fluid initial conditions:
Verification | Dynamic Impact | |
rhol [kg/m^3] | 1.2 | 1.6458 |
ul [m/s] | 0.0 | 112.61 |
pl [kPa] | 200.0 | 156.18 |
rhor [kg/m^3] | 1.2 | 1.2 |
ur [m/s] | 0.0 | 0.0 |
pr [kPa] | 100.0 | 100.0 |
Location relative to beam | 0.0 | -0.02m |
-- RalfDeiterding - 12 Jul 2007
I | Attachment | Size | Date | Who | Comment |
---|---|---|---|---|---|
Schl50.mpg | 7521.7 K | 12 Feb 2007 - 16:59 | RalfDeiterding | Beam length 50mm, schlieren plot of fluid density in whole domain | |
SchlL50.mpg | 10617.0 K | 12 Feb 2007 - 17:00 | RalfDeiterding | Beam length 50mm, schlieren plot of fluid density in whole domain, zoom into beam, shows beam motion clearly | |
Dist50.mpg | 3254.0 K | 12 Feb 2007 - 17:02 | RalfDeiterding | Beam length 50mm, change of distribution to fluid processors | |
Schl40.mpg | 9925.2 K | 12 Feb 2007 - 17:03 | RalfDeiterding | Beam length 40mm, schlieren plot of fluid density in whole domain | |
SchlL40.mpg | 10176.0 K | 12 Feb 2007 - 17:03 | RalfDeiterding | Beam length 40mm, schlieren plot of fluid density in whole domain, zoom into beam |