Products: ABAQUS/Standard ABAQUS/Explicit
The underwater shock analysis capability in ABAQUS is provided by the *INCIDENT WAVE loading option. The coupling is accomplished by using *TIE in which ABAQUS calculates both the structural response and the fluid pressure response at the interaction surface.
A one–dimensional continuum model is analyzed as described below.
The ABAQUS model (Figure 1.13.1–1) consists of a single C3D8R continuum element constrained to deform in one dimension. The stiffness and density have been chosen so that the structure alone has a natural frequency of 1 Hz. A single acoustic element is coupled to the structural element. The continuum element has unit cross sectional area and a thickness of 10 units. Two cases are analyzed here. In the first case the acoustic element thickness is set to 0.01 units. In the second case the thickness is increased to 0.1 units. The surface of the acoustic element and the structural element are tied together using *TIE. The fluid density and speed of sound are both set to 1.0. A plane wave pressure pulse is applied to the front face of the continuum element and the back face of the acoustic element using the *INCIDENT WAVE option. The back face of the continuum element is held fixed and a plane wave absorbing boundary condition is specified on the front face of the acoustic element via the *SIMPEDANCE option. The pressure pulse travels along the x-axis and is a step function in time. The sound source is located at (100, 0, 0) for the plane waves, and the stand-off point is located at (10, 0, 0). The analysis is run for 10 seconds. The response of the front face of the continuum element is one-dimensional and oscillatory. The analysis is performed in both ABAQUS/Standard (u1d_std_c3d8r_ac3d8.inp) and ABAQUS/Explicit (u1d_xpl_c3d8r_ac3d8r.inp). There is no damping applied in either the structure or the fluid model except from the radiation boundary condition. However, in ABAQUS/Explicit bulk viscosity is introduced in the model by using *BULK VISCOSITY to introduce damping. The value of linear bulk viscosity is chosen as 0.02, and the value of quadratic bulk viscosity is chosen as 0.5.
The model data are kept the same for the restart analysis. In the initial run the loading is applied for 2 seconds. During the restart run, the initial conditions are read from the restart files for the new *DYNAMIC step. The loading is applied for 2 seconds. A restart analysis is performed for both ABAQUS/Standard and ABAQUS/Explicit.
The response comparison is based on the velocity of the front face of the structure. History plots of velocity versus time for ABAQUS/Standard and ABAQUS/Explicit are shown in Figure 1.13.1–2 and Figure 1.13.1–3. The plots are compared individually to the reference solution. Both ABAQUS/Standard and ABAQUS/Explicit exactly match the direct integration reference solution.
FORTRAN program used to generate the direct integration numerical solution for the one-dimensional continuum problem. The program uses trapezoidal integration.
ABAQUS/Standard analysis: C3D8R/AC3D8 model with acoustic element of thickness 0.01 units.
ABAQUS/Explicit analysis: C3D8R/AC3D8R model with acoustic element of thickness 0.01 units.
Initial ABAQUS/Standard run used for restart analysis.
Restart analysis.
Initial ABAQUS/Explicit run used for restart analysis.
Restart analysis.
ABAQUS/Standard analysis: C3D8R/AC3D8 model with acoustic element of thickness 0.1 units.
ABAQUS/Explicit analysis: C3D8R/AC3D8R model with acoustic element of thickness 0.1 units.
A one-dimensional plate model is analyzed as described below.
The ABAQUS model (Figure 1.13.1–4) consists of a single S4R shell element constrained to translate as a rigid body in one dimension. The planar shell is modeled in the X–Y plane with a length of 1.5 units on all sides. The shell properties are those of steel. The *TIE option is used to couple the structure to the acoustic fluid. The fluid modeled is water. A plane wave underwater pressure pulse is applied to the front face of the structural element and to the back face of the acoustic element. The source is at (0.75, 0.75, 100), and the stand-off point is at (0.75, 0.75, 0). The pulse travels along the z-axis and is a step function in time. The reaction of the plate is rigid body translation in one dimension with constant acceleration. The response comparison is based on the velocity time history of the plate. The analysis is performed in both ABAQUS/Standard (u1d_std_s4r_ac3d8.inp) and ABAQUS/Explicit (u1d_xpl_s4r_ac3d8r.inp).
The ABAQUS S4R/AC3D8 model results are identical to the reference solution. History plots of velocity versus time for ABAQUS/Standard and ABAQUS/Explicit are shown in Figure 1.13.1–5 and Figure 1.13.1–6.
ABAQUS/Standard analysis: S4R/AC3D8 model.
ABAQUS/Explicit analysis: S4R/AC3D8R model.
