HyperWorks Solvers

/INTER/TYPE19

/INTER/TYPE19

Previous topic Next topic Expand/collapse all hidden text  

/INTER/TYPE19

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

Block Format Keyword

/INTER/TYPE19 - Interface Type 19

Description

This is a combination of interface TYPE7 and TYPE11, with common input based on the same slave/master surfaces. Slave node group for interface TYPE7, as well as slave and master line segments used by equivalent TYPE11 interface are virtually generated from these input surfaces.

Format

(1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

(9)

(10)

/INTER/TYPE19/inter_ID/unit_ID

inter_title

surf_IDs

surf_IDm

Istf

 

Igap

 

Ibag

Idel

Icurv

 

Fscalegap

Gap_max

 

 

 

 

 

 

Stmin

Stmax

 

 

dtmin

 

 

 

Insert if Icurv = 1 or 2

(1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

(9)

(10)

node_ID1

node_ID2

 

 

 

 

 

 

 

 

 

Required Fields

(1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

(9)

(10)

Stfac

Fric

Gapmin

Tstart

Tstop

IBC

 

 

Inacti

VISS

VISF

Bumult

Ifric

Ifiltr

Xfreq

Iform

sens_ID

 

 

 

 

 

Insert if Ifric > 0 (Optional)

(1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

(9)

(10)

C1

C2

C3

C4

C5

 

Read this input only if Ifric > 1 (Optional)

(1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

(9)

(10)

C6

 

 

 

 

 

 

 

 

hmtoggle_plus1Flag Definition

Field

Contents

SI Unit Example

inter_ID

Interface identifier

(Integer, maximum 10 digits)

 

unit_ID

Optional unit identifier

(Integer, maximum 10 digits)

 

inter_title

Interface title

(Character, maximum 100 characters)

 

surf_IDs

Slave surface identifier

(Integer)

 

surf_IDm

Master surface identifier

(Integer)

 

Istf

Stiffness definition flag (Comment 6)

(Integer)

= 0: default, set to the value defined in /DEFAULT/INTER/TYPE19

= 1: Stfac is a stiffness value

= 2, 3, 4 and 5: Stfac is a stiffness scale factor and the stiffness is computed from both master and slave characteristics

= 6: Stfac is a stiffness scale factor (default in case /DEFAULT/INTER/TYPE19 is not defined)

For node to 3 and 4 node segments contacts, the stiffness is computed according to the master side characteristics;

For 2 node segments to 2 node segments contacts, Istf = 5 is used (Stfac is a stiffness scale factor and the stiffness is computed from both master and slave characteristics).

 

Igap

Gap/element option flag

(Integer)

= 0: default, set to the value defined in /DEFAULT/INTER/TYPE19

= 1: gap varies accordingly to the characteristics of the impacted master surface and the impacting slave node

= 2: variable gap + gap scale correction of the computed gap

= 1000: gap is constant and equal to the minimum gap (default in case /DEFAULT/INTER/TYPE19 is not defined)

 

Ibag

Airbag vent holes closure flag in case of contact

Default = 0  (Integer)

= 0: default, set to the value defined in /DEFAULT/INTER/TYPE19

= 1: closure

= 2: no closure (default in case /DEFAULT/INTER/TYPE19 is not defined)

 

Idel

Node and segment deletion flag (Comment 4)

Default = 0 (Integer)

= 0: default, set to the value defined in /DEFAULT/INTER/TYPE19

= 1: when all the elements (4-node shells, 3-node shells, or solids) associated to one segment are deleted, the segment is removed from the master side of the interface. It is also removed in case of explicit deletion using RADIOSS Engine keyword /DEL in the Engine file.

Additionally, non-connected nodes are removed from the slave side of the interface.

= 2: when a 4-node shell, a 3-node shell or a solid element is deleted, the corresponding segment is removed from the master side of the interface. It is also removed in case of explicit deletion using RADIOSS Engine keyword /DEL in the Engine file.

Additionally, non-connected nodes are removed from the slave side of the interface.

= 3: no deletion (default)

= -1: same as = 1, except non-connected nodes are not removed from the slave side of the interface.

= -2: same as = 2, except non-connected nodes are not removed from the slave side of the interface.

=1000: no deletion (default in case /DEFAULT/INTER/TYPE19 is not defined)

 

Icurv

Slave gap with curvature (Comments 9, 10 and 11 )

(Integer)

= 0: no curvature

= 1: spherical curvature

= 2: cylindrical curvature

= 3: automatic bicubic surface

 

Fscalegap

Gap scale factor

Default = 1.0  (Real)

 

Gap_max

Maximum gap

= 0: there is no maximum value for the gap.

(Real)

Stmin

Minimum stiffness

(Real)

Stmax

Maximum stiffness

Default = 1030  (Real)

dtmin

Limiting nodal time step (Comment 22)

(Real)

node_ID1

First node identifier

(Integer)

 

node_ID2

Second node identifier

(Integer)

 

Stfac

Interface stiffness (if Istf = 1)

Default = 1.0

(Real)

Stiffness scale factor for the interface (if Istf 1)

Default = 0.0

(Real)

 

Fric

Coulomb friction

(Real)

 

Gapmin

Minimum gap for impact activation

(Real)

Tstart

Start time

(Real)

Tstop

Time for temporary deactivation

(Real)

IBC

Deactivation flag of boundary conditions at impact

(Boolean)

 

Inacti

Deactivation flag of stiffness in case of initial penetrations (Comment 14)

(Integer)

= 0: default, set to the value defined in /DEFAULT/INTER/TYPE19

= 1: deactivation of stiffness on nodes

= 2: deactivation of stiffness on elements

= 3: change node coordinates to avoid initial penetrations

= 5: gap is variable with time and initial gap is computed as follows:

, with the initial penetration

= 6: gap is variable with time but initial penetration is computed as follows (the node is slightly depenetrated):

= 1000: no action (default in case /DEFAULT/INTER/TYPE19 is not defined)

 

VISS

Critical damping coefficient on interface stiffness

Default set to 0.05  (Real)

 

VISF

Critical damping coefficient on interface friction

Default set to 1.0  (Real)

 

Bumult

Sorting factor

Default set to 0.20  (Real)

 

Ifric

Friction formulation flag (Comments 18 and 19)

(Integer)

= 0: static Coulomb friction law

= 1: generalized viscous friction law

= 2: Darmstad friction law

= 3: Renard friction law

 

Ifiltr

Friction filtering flag (Comment 20)

(Integer)

= 0: no filter is used

= 1: simple numerical filter

= 2: standard -3dB filter with filtering period

= 3: standard -3dB filter with cutting frequency

 

Xfreq

Filtering coefficient (Comment 20)

(Real)

 

Iform

Friction penalty formulation type

(Integer)

= 0: default, set to the value defined in /DEFAULT/INTER/TYPE19

= 1: viscous (total) formulation (default in case /DEFAULT/INTER/TYPE19 is not defined)

= 2: stiffness (incremental) formulation

 

sens_ID

Sensor identifier to Activate/Deactivate the interface (Comment 1)

(Integer)

If an identifier sensor is defined, the activation/deactivation of interface is based on sensor and not with Tstart or Tstop.

 

C1

Friction law coefficient

(Real)

 

C2

Friction law coefficient

(Real)

 

C3

Friction law coefficient

(Real)

 

C4

Friction law coefficient

(Real)

 

C5

Friction law coefficient

(Real)

 

C6

Friction law coefficient

(Real)

 

hmtoggle_plus1 Flags for Deactivation of Boundary Conditions: IBC

(1)-1

(1)-2

(1)-3

(1)-4

(1)-5

(1)-6

(1)-7

(1)-8

 

 

 

 

 

IBCX

IBCY

IBCZ

 

Field

Contents

SI Unit Example

IBCX

Deactivation flag of X boundary condition at impact

=0: free DOF

=1: fixed DOF

(Boolean)

 

IBCY

Deactivation flag of Y boundary condition at impact

=0: free DOF

=1: fixed DOF

(Boolean)

 

IBCZ

Deactivation flag of Z boundary condition at impact

=0: free DOF

=1: fixed DOF

(Boolean)

 

hmtoggle_plus1Comments
1.When sens_ID is defined for activation/deactivation of the interface, Tstart and Tstop are not taken into account.
2.In case of SPMD, each master segment defined by surf_IDm must be associated to an element (possibly to a void element).
3.For flag Ibag, refer to the monitored volume option (/MONVOL keyword).
4.Flag Idel = 1 has a CPU cost higher than Idel = 2.
5.If Igap = 2, the variable gap is computed as:

while,

gm: master element gap

, with t: thickness of the master element for shell elements

gm = 0 for brick elements

gs: slave node gap:

gs = 0 if the slave node is not connected to any element or is only connected to brick or spring elements.

, with t: largest thickness of the shell elements connected to the slave node.

for truss and beam elements, with S being the cross section of the element.

Igap = 2 is not available in interface TYPE11.

6.Contact stiffness

For node to 3-node and 4-node segments or 2-node segments to 2-node segments contacts computation as follows:

Where,

Kn is computed from both master segment stiffness Km and slave node stiffness Ks:

Istf = 2,

Istf = 3,

Istf = 4,

Istf = 5,

Km is master segment stiffness and computed as follows:

When master segment lies on a shell or is shared by shell and solid:

When master segment lies on a solid:

Where,

S is the segment area,

V is the volume of the solid,

B is the Bulk Modulus

Ks is an equivalent nodal stiffness considered for interface TYPE7, and computed as:

when node is connected to a shell element:

when node is connected to solid element:

There is no limitation value to the stiffness factor Stfac (but a value can be larger than 1.0 to reduce the initial time step)

7.The values given in Line 4 are ignored, if Igap 2.
8.The values given in Line 5 are ignored, if Istf < 1.
9.Spherical curvature (Icurv = 1) is defined with node_ID1 (center of the sphere).
10.The node_ID2 given in Line 6 is ignored, if Icurv = 1.
11.Cylindrical curvature (Icurv = 2) is defined with node_ID1 and node_ID2 (on the axis of the cylinder).
12.A default value for Gapmin is computed as the minimum of:

While,

: master surface gap

t: average thickness of the master elements for shell elements.

l: average side length of the master brick elements.

lmin: the smallest side length of all master segments (shell or brick).

: slave surface gap: computation identical to ; except that it is applied on slave side elements.

13.The gap is computed for each impact as:

with:

gm: master element gap:

, with t is thickness of the master element for shell elements

for brick elements

gs: slave node gap:

if the slave node is not connected to any element or is only connected to brick or spring elements.

, with t is the largest thickness of the shell elements connected to the slave node.

, with S being the cross section of the truss and beam element.

If the slave node is connected to multiple shells and/or beams or trusses, the largest computed slave gap is used.

The variable gap is always at least equal to Gapmin.

14.Inacti = 3 may create initial energy if the node belongs to a spring element.

Inacti = 5 is recommended for airbag simulation deployment.

Inacti = 6 is recommended instead of Inacti =5, in order to avoid high frequency effects into the interface.

Inacti_flag_6

15.The sorting factor, Bumult is used to speed up the sorting algorithm.
16.The sorting factor, Bumult is machine dependent.
17.One node can belong to the two surfaces at the same time.
18.For friction formulation
If the friction flag Ifric = 0 (default), the old static friction formulation is used:

While,

with is Coulomb Friction coefficient.

For flag Ifric > 0, new friction models are introduced. In this case, the friction coefficient is set by a function

Where,

p is the pressure of the normal force on the master segment,

V is the tangential velocity of the slave node.

19.Currently, the coefficients C1 ~ C6 are used to define a variable friction coefficient for new friction formulations.
Ifric = 1 (generalized viscous friction law):

Ifric = 2  (Darmstad law):

Ifric = 3  (Renard law):

 if

 if

if

Where,

       

   

     

First critical velocity must be different to 0 ().
First critical velocity must be lower than the second critical velocity ().
The static friction coefficient and the dynamic friction coefficient , must be lower than the maximum friction ( and ).
The minimum friction coefficient , must be lower than the static friction coefficient and the dynamic friction coefficient ( and )

 

20.Friction filtering

If Ifiltr ≠ 0, the tangential forces are smoothed using a filter:

where coefficient is calculated from:

If Ifiltr = 1 , simple numerical filter
If Ifiltr = 2 , standard -3dB filter, with , and T is filtering period
If Ifiltr = 3 , standard -3dB filter, with Xfreq is cutting frequency

The filtering coefficient Xfreq should have a value between 0 and 1.

21.Friction penalty formulation Iform
If Iform = 1, (default) viscous formulation, the friction forces are:

While an adhesion force is computed as follows:

with

If Iform = 2, stiffness formulation, the friction forces are:

While an adhesion force is computed as follows:

with

Where,

Vt is contact tangential velocity

22.Slave segment is deactivated from the contact when the segment kinematic time step calculated for this contact becomes smaller than dtmin.

See Also:

Interface TYPE7 in User’s Guide

Interface TYPE11 in User’s Guide

Penalty Method in User's Guide

Interface TYPE7 in Theory Manual

Interface TYPE7 in FAQ

Interface TYPE11 in FAQ