HyperWorks Solvers

22.1 - Ditching using SPH (Mono-Domain)

22.1 - Ditching using SPH (Mono-Domain)

Previous topic Next topic No expanding text in this topic  

22.1 - Ditching using SPH (Mono-Domain)

Previous topic Next topic JavaScript is required for expanding text JavaScript is required for the print function  

 

Title

Ditching using SPH (Mono-Domain)

rad_ex_22.1

Number

22.1

Brief Description

Impact of a simple object on water.

Keywords

SPH modeling and hexagonal net
Hydrodynamic viscous fluid law (/MAT/LAW6) and impact on water modeling
Type 7 interface

RADIOSS Options

Rigid body (/RBODY)
Initial velocity (/INIVEL)
Accelerometer (/ACCEL)
Gravity (/GRAV)
Interface (/INTER)
SPH outlet (/SPH/INOUT)

Compared to / Validation Method

Experimental data provided by Politecnico di Milano (Polytechnic University of Milan) [1]
Analytic solution proposed by Von Karman [2]

Input File

Impact_velocity=3.5m/s: <install_directory>/demos/hwsolvers/radioss/22_Ditching/Ditching_Mono_Domain_SPH/v_35/*

Impact_velocity=6.8m/s: <install_directory>/demos/hwsolvers/radioss/22_Ditching/Ditching_Mono_Domain_SPH/v_68/*

Impact_velocity=11m/s: <install_directory>/demos/hwsolvers/radioss/22_Ditching/Ditching_Mono_Domain_SPH/v_110/*

Technical / Theoretical Level

Skilled

Overview


Physical Problem Description

The problem consists of a simple object falling into water simulating the ditching of a helicopter.

Units: mm, ms, KN, GPa, kg.

Impact of a triangular section object on water is performed and the results are compared qualitatively [2], also using the experimental data obtained from the Politecnico di Milano [1].

The computation is performed using several impact velocities: 3.5 m/s, 6.8 m/s and 11 m/s.

The material used for the object follows a linear elastic law (/MAT/LAW1) with the following characteristics:

Initial density: 7.8 x 10-6 kg.mm-3
Young modulus: 206.9 GPa
Poisson ratio: 0.3

The material law for water is a hydrodynamic viscous fluid law (/MAT/LAW6) with the following characteristics:

Initial density: 1 x 10-6 kg.mm-3

Kinematic viscosity: 0Pressure cutoff:  -0.0001 GPa
C0:  0 GPaPressure shift:  0 GPa
C1:  2199 GPaC4:  0 GPa
C2:  5351 GPaC5:  0 GPa
C3:  7324 GPaInitial energy/unit of volume:  0 mJ/mm3

ex22_fig1

Fig 1: Problem data.

Analysis, Assumptions and Modeling Description


Modeling Methodology

The object is modeled using shell elements with an average mesh size of 15 x 15 mm2.

The water is modeled using SPH particles having a hexagonal compact net with a smoothing length "ho" equal to 28.2843 mm. Each particle of the net represents a volume equal to 16 mm3 and weighs 16 g. This part uses 36075 SPH cells.

The size of the water block is adapted to the shape of the object for the purpose of reducing the model’s size and the simulation’s CPU time.

RADIOSS Options Used

Rigid body:

The object is modeled using a rigid body, and a mass of 23.0422 kg is added to the rigid body’s master node (ID: 287002).

Initial velocity:

An initial velocity, in accordance with the Z-axis, is set on the rigid body’s master node and its value is set successively at 3.5 m/s, 6.8 m/s and 11 m/s.

Gravity:

A gravity load  gz = -9.81 m.s-2 is applied to the object.

Accelerometer:

An accelerometer is set on the rigid body’s master node.

SPH outlet:

The parallelepiped water mesh is surrounded on five faces by Outlet SPH boundary frontier absorbing conditions. A control surface is placed at a distance equal to 2 x ho inside the water. This surface shown in green on Fig 2, is oriented so that its normal vector points face the interior of the domain. On this outlet surface, specific non-reflective frontiers (NRF) are applied to the SPH cells.

Interfaces:

One type of contact occurs in the simulation. Contact between the skin structure (shell finite elements) and the water (SPH cells) is modeled using a sliding interface (type 7). The gap between the surface skin and the SPH cells is equal to 3 mm. After optimization, a scale factor on the Penalty stiffness interface equal to 0.1 is used for controlling the interface forces between the rigid object and the water.

rad_ex_fig_22-2

Fig 2: SPH outlet boundary surface non-reflective frontiers (NRF) (green).

Simulation Results and Conclusions


Curves and Animations

Output Pressure

The SPH module does not enable output pressure time history values; only the animation values may be compared to the experimental data. Figures 3 and 4, show a variation in the pressure when the wedge ditches into the SPH water, corresponding to an impact velocity of 6.8 m/s. In order to attain a representative pressure gradient, the pressure range has been fixed at -0.01 GPa for the minimum value and 0.08 GPa for the maximum value. In comparison to the experimental data, the pressure values are basically the same, the wave propagation being well described.

rad_ex_fig_22-3

Fig 3: Ditching simulation (from the beginning to 3 ms).

rad_ex_fig_22-4

Fig 4: Ditching simulation (from 4 ms to 9 ms).

Output Acceleration

For the object, an accelerometer is set on the master node of the corresponding rigid body. The acceleration values expressed in g units are compared to both the experimental values [1] and the analytic solution proposed by Von Karman [2]. The signal is filtered using a CFC 60 (-3db) filter frequency after calculation. The filtering reduces discrepancy between the peaks.

The following diagrams indicate the time history acceleration results at the wedge object’s rigid body’s master node for the three cases of impact speed: 3.5, 6.8 and 11 m/s.

rad_ex_fig_22-5

Fig 5: Deceleration of the wedge for an impact velocity of 3.5 m/s.

rad_ex_fig_22-6

Fig 6: Deceleration of the wedge for an impact velocity of 6.8 m/s.

rad_ex_fig_22-7

Fig 7: Deceleration of the wedge for an impact velocity of 11 m/s.

For these three cases, the SPH approach using the OUTLET SPH boundary conditions indicates a good deceleration. For an impact velocity nearing the 8 m/s of the Helicopter ditching configuration, the deceleration is in correlation with the experimental data [1] and also with the analytic solution proposed by Von Karman [2].

Conclusion

The simulations show that the SPH approach using the OUTLET option, allows the ditching of simple objects to be modeled without any numerical problems.

The SPH and OUTLET results are very close to the experimental test results and also to the analytical solution. In conclusion, to achieve ditching simulations with the correct results, it is necessary to model the water block using the SPH method with the OUTLET boundary conditions.

References

[1] CAST Deliverable 5.5.1 Generic Water Impact Tests performed at Politecnico di Milano (Polytechnic University of Milan)

[2] Olivier Pastore Study and modelization of rigid bodies impact during sea landing phase; Annex 1 Von Karman's Theoretical Models, T. Miloh et al. May.