HyperWorks Solvers

Example 4 - Airbag

Example 4 - Airbag

Previous topic Next topic No expanding text in this topic  

Example 4 - Airbag

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

rad_ex_4_airbag

Summary


This example deals with the deployment of a chambered airbag modeled by monitored volumes using communications. The airbag is initially folded along four fold lines. The fabric is meshed with shell elements which undergo an elastic orthotropic behavioral test. Perfect gas is injected into a central chamber via an inflator with the air flow through the connected chambers being simulated. The chambers inflate while the airbag is deploying.

In the self-impacting interface definition, the action of the Inacti flag to deactivate stiffness in the case of initial penetration is studied in order to significantly increase the time step. An adequate gap enables to pass from a kinematic interface time step to a higher element time step.

 

Title

Airbag

rad_ex_4_airbag

Number

4.1

Brief Description

A chambered airbag folded along four fold lines is deployed.

Keywords

Orthotropic shell
Monitored volumes and communicating airbags
Hierarchy organization

RADIOSS Options

Monitored volume with communications (/MONVOL/COMMU1)
Interface (/INTER/ with Inacti flag)

Input File

Inactiv_0_Gap0.1: <install_directory>/demos/hwsolvers/radioss/04_Airbag/Inacti0_Gap01/AIRFIX*

Inactiv_5_Gap0.3: <install_directory>/demos/hwsolvers/radioss/04_Airbag/Inacti5_Gap03/AIRBAG*

Inactiv_5_Gap1.5: <install_directory>/demos/hwsolvers/radioss/04_Airbag/Inacti5_Gap15/AIRBAG2*

Technical / Theoretical Level

Beginner

Overview


Aim of the Problem

The purpose of this example is to deal with monitored volumes using communications on a simple airbag model. Methods for increasing the time step are considered.

Physical Problem Description

A 30-liter airbag is folded along the four fold lines. The following examples illustrate the airbag folded and deployed.

rad_ex_fig_4.1

Fig 1: Folded airbag Fig 2: Deployed airbag

The fabric thickness is 0.33 mm and is modeled using an elastic orthotropic material law (/MAT/LAW19) with the following properties:

Density: 0.85x10-3 g/mm3
Young’s Modulus: 500 MPa in both directions
Shear Modulus: 10 MPa
Reduction factor: 0.001

The property set is /PROP/SH_ORTH (shell orthotropic, type 9), using one integration point.

Analysis, Assumptions and Modeling Description


Modeling Methodology

rad_ex_fig_4-3

Fig 3: Overall mesh of the structure (folded and deployed).

The model is divided into two subsets: the fabric layers and the communication surfaces.

The fabric surface is then divided into 9 subsets, one for each monitored volume. Each "monitored volume" is further divided into two parts. All the parts of the layer of fabric have the same Type and MID.

The same properties apply for the communication surfaces.

rad_ex_fig_4-4

Fig 4: Folder airbag with communications.

The airbag is modeled using 9 communicating volumes in order to simulate the air flow through the folds and the behavioral differences within the airbag when unfolding. The communicating surfaces between the volumes are simulated using dummy membranes. The dummy membranes are modeled using shells with fictitious material (/MAT/LAW0).

RADIOSS Options Used

A monitored volume is defined as a surface area having one or more shell property sets and where the surface must be closed. The monitored volume used is a COMMU1 type for airbags using communications (chambered, with communications, of the folder airbag type). For further details about monitored volumes, see the RADIOSS Theory Manual.

The main properties for this type are:

Volumetric damping factor: 0.001 g.mm-1ms-1 .
External pressure: 0.1 MPa
Constant perfect gas: 1.4
Specific heat at constant pressure: 926 mJ/g (This is the specific heat coefficient related to mass)
Temperature:  780 K.
Communication area: total (Acom =1 and Scom =0)
Time to deflate vent hole: 1030 ms

The gas molecular weight and specific heat coefficients are defined in /MAT/GAS (using MASS type):

Molecular weight of gas MW: 30.09204 g/kmol (AIR)
Specific heat at constant pressure: 926 mJ/g (specific heat coefficient related to mass)

Specific input for the central chamber one (inflator):

Vent hole membrane surface area is 1000 mm2 (Avent =0) and is immediately activated.
Relative vent deflation pressure: 0.0002
Number of injectors: 1 (Njet =1; Ijet =0)

Using /PROP/INJECT1 to described mass inject:

Final injected mass is 46 g injected into the central chamber (FscaleM and FscaleT =1). Two functions define the mass and temperature of the injected gas compared with time (function identifiers: fct_IDM and fct_IDT).

Time (ms)

0

2

4

5

6

8

11

12

15

19

28

30

106

Mass (g)

0

6

11

14

17

22

29

31

36

41

45

46

46

Injected mass function.

Time (ms)

0

106

Temperature (K)

780

780

Temperature of injected gas function.

Interface

Taking into account the fabric is self-impacting with itself, an self-impacting interface must be used. The interface’s Block Format definition is made: defining the master surface (/SURF/PART), then defining the slave nodes for all nodes on this surface (/GRNOD/SURF).

The distance between the fabric layers before unfolding is very small. In order to avoid initial penetration, the gap required is approximately 0.1 mm, thus enabling the time step to considerably decrease when such a gap is chosen.

By using Inacti =5, a 0.3 mm gap is chosen. Any initial penetration below 0.2 mm (two-thirds of the input gap) is ignored (it is strongly recommended to verify that no initial penetration is above this value).

Simulation Results and Conclusions


Curves and Animations

In order to demonstrate the interest of the Inacti flag, the same model was run with Inacti at a value of 0, with a gap of 0.1 mm (no initial penetration).

rad_ex_fig_4-5

Fig 5: Comparison between option Inacti = 5 and Inacti = 0

Using Inacti = 5, the minimum time step is around 10-3 ms. When not using this option, the minimum time step is around 2x10-4 ms. For the full model, the number of cycles may be divided up into 10 or more. Furthermore, the model is numerically less sensitive.

The time step is monitored by the interface time step (kinematic) for up to 40 ms despite the unfolding and the fact that there is no energy contact from 7.8 ms. In order to transfer into the element time step and to reduce computation time, it is advisable to increase the gap so the kinematic step becomes higher than the element step.

Time-stepkinetic < 0.9 x GAP / Nodal_velocityrelative  (using scale factor = 0.9)

The time step is only low during the unfolding phase (before 10 ms) with a gap equal to 1.5 mm.

 

Inacti flag = 5

Inacti flag = 0

GAP = 0.3 mm

GAP = 1.5 mm

GAP = 0.1 mm

Error on energy

-16.3%

-19.6%

-15.5%

Elapsed Time [s]

57.5

47.8

78.58

Airbag deploy

completely deploy

completely deploy

incompletely deploy

rad_ex_fig_4-6

Fig 6: Time step obtained with GAP = 0.3 mm and GAP = 1.5 mm (Inacti = 5).

rad_ex_fig_4-7

Fig 7: Contact energy with GAP = 0.3 mm and GAP = 1.5 mm (Inacti = 5).

It is obvious that a gap of 1.5 mm generates an increase in the contact force. However, the additional error on energy remains quite low and is acceptable.

rad_ex_fig_4-8

Fig 8: Time history of pressure.

rad_ex_fig_4-9

Fig 9: Time history of volume.

Animations

rad_ex_fig_4-10

Fig 10: Central chamber is inflating.

rad_ex_fig_4-11

Fig 11: All chambers are inflating.

rad_ex_fig_4-13

Fig 12: Airbag is deployed