HyperWorks Solvers

Time Step

Time Step

Previous topic Next topic Expand/collapse all hidden text  

Time Step

Previous topic Next topic JavaScript is required for expanding text JavaScript is required for the print function  
hmtoggle_plus1Why is the cycle time step very small?

First, check which object is responsible for the low time step in the listing file (Runname_nnnn.out) of RADIOSS Engine. It can be an element, a node or an interface.

For an element, check the related material (especially its Young modulus and density in case of an elastic-plastic material; and its viscosity in case of a visco-elastic material). There must not be an error in the units system that this data is given in.

Check the size of the element, since

where, lc is the characteristic length of the element, c is the material sound speed, and dtsca is the scale factor.

For a node, check the characteristics of connected elements. If the node is on the master side or the slave side of an interface, this interface must be verified.

If a master node of a rigid body gives the time step, check the rigid body inertia reported in the listing file (Runname_0000.out), at:

Rigid Body Initialization (for a rigid body made of only a few nodes, it may be necessary to set a spherical inertia).

For an interface, the gap of the interface must be verified if some failure happens on the master or the slave side of the interface.
hmtoggle_plus1Using beam elements (/PROP/BEAM) the time step given by the Starter in listing file (Runname_0000.out) and the time step given by the Engine on this element are different; is this normal?

For this time step computation, RADIOSS Starter does not take into account the damping coefficients; while RADIOSS Engine does.

Particularly for thick beams, due to flexural damping, the time step computed by Engine may be quite different from the one computed by Starter.

Note:The beam formulation in property type /PROP/BEAM is not well-suited for thick beams.
hmtoggle_plus1How is the time step computed on master nodes of rigid bodies?

For a rigid body and nodal time step, the time step for stability is computed at the master node (the nodal time step at slave nodes does not need to be considered).

The time step computed at the master node of the rigid body is:

Where,

M is the mass of the rigid body

l is the lower principal inertia of the rigid body inertia matrix around the master node

K is the stiffness reported to the master node

S is the slave node

Kr  is the rotation stiffness reported to the master node

S is the slave node

D being the distance from slave to master node.

The time step at the master node of the rigid body is computed by reporting stiffness at slave nodes to the master node.

The mass and inertia of the rigid body depends on the flag ICoG in /RBODY.

If ICoG =4, the mass that is input in the rigid body option is considered; whereas for ICoG =2, the mass that is input in the rigid body option is added to the slave nodes mass.

hmtoggle_plus1The cycle time step reported in listing file (Runname_nnnn.out), the time step is found on a master node of a rigid body; is this normal and how to increase them?

Yes, it is possible. First check that the inertia which has been given to the rigid body is correct.

Check the rigid body inertia as it is reported in the listing file (Runname_0000.out), at:

RIGID BODY INITIALIZATION

If the conditioning of the inertia matrix is bad, either because the slave nodes are aligned or because there are only a few slave nodes, it is recommended to set Ispher =1 for the rigid body in order to use a spherical inertia.

hmtoggle_plus1How is it possible that the time step is controlled by a node even while using elementary time step(i.e did not specify /DT/NODA)?

CYCLE

TIME

TIME-STEP

ELEMENT

ERROR

I-ENERGY

K-ENERGY T

K-ENERGY R

EXT-WORK

MAS.ERR

100

0.1846E0-02

0.1846E-04

NODA   123

0.0%

0.4261E-04

0.6069E-02

0.3660E-25

0.6112E-02

0.000

The input file probably includes an Interface Type 7, Type 10 or Type 11.

faq_timestep5

For an impact between a slave node and a master segment, the stiffness of the penalty spring for this impact is applied to both slave and master nodes, in order to ensure stability.

More precisely, Interface Type 7 associates a penalty spring and damping, so that an equivalent stiffness to the whole system is considered.

The nodal stability time step is computed as follows:

where, M is the nodal mass and K is the nodal stiffness

is the sum of the stiffness applied to the node by all elements and contact interfaces.

If the deck includes Interface Type 7, Type 10 or Type 11, this time step is always computed at all nodes of the deck, in order to ensure stability even if an elementary time step is used.

When the time step is taken onto a node, this node is probably impacting (otherwise nodal time step computed with faq_timestep5C is greater than elementary time step).

hmtoggle_plus1Using /DT/INTER/DEL in RADIOSS Engine file (Runname_0001.rad).  How is it possible that the time step is given in the listing file (Runname_0001.out) by an interface, and is lower than the minimum value given in /DT/INTER/DEL?

faq_timestep6

The interface kinematic time step avoids the slave nodes to cross the master surface.

1.A kinematic time step is computed for the whole interface, as follows:

where vr is an estimation of the maximum relative speed between master nodes and slave nodes.

This time step computation prevents a slave node from crossing the master surface while ensuring the performances of searching the contacts through the interface.

2.For the nodes which lie within the gap of a master segment, the kinematic time step is computed more precisely for this slave/master couple, as follows:

If this time step is lower than the value given in option /DT/INTER/DEL, the slave node is suppressed from the interface and the following message is written to the listing file of RADIOSS Engine:

**WARNING MINIMUM TIME STEP 0.3001E-03 IN INTERFACE 1

REMOVE SLAVE NODE 526210 FROM INTERFACE

3.For the nodes which lie within the gap of a master segment, the stability time step is computed for this slave/master couple. If it is lower than the minimum time step given in option /DT/INTER/DEL, the slave node is also suppressed from the interface.

In certain situations, the kinematic time step of the interface may be lower than the minimum time step which was specified in option /DT/INTER/DEL. If such a situation is encountered from the beginning of the computation, the value of the Gap may be checked.

hmtoggle_plus1What is the difference between the messages "Delete Slave Node" and "Remove Slave Node" concerning interfaces Type 7 in the listing file (Runname_0001.out)?

These messages are related to the option /DT/INTER/DEL. They both consist in the slave node suppression from the interface. This node will no longer impact any master segment for this interface.

faq_timestep7

The message "Remove Slave Node" is written when the interface kinematic time step is lower than the minimum time step; which was specified in option /DT/INTER/DEL. For the slave/master couple which is considered, the kinematic time step is:

If it is such that , then the slave node is suppressed from the interface.

The message "Delete Slave Node" is written when the stability time step of the interface for the slave/master couple; which is considered is lower than the minimum time step which was specified in option /DT/INTER/DEL.

The stability time step for the slave/master couple is computed as follows:

where, M is the minimum nodal mass among the slave node and the nodes of the master segment, and Ki is the interface stiffness for the slave/master couple which is considered (Ki is an equivalent stiffness to the penalty spring and the damping which are applied between the slave node and the master segment).

When the stability time step for the slave/master couple is such that , then the slave node is suppressed from the interface.

It is possible that those messages would be written several times within the same cycle for the same node and the same interface (if the node impacts several master segments within the same cycle and several contacts get a lower time step than the minimum time step which was specified in option /DT/INTER/DEL).