/INTER/LAGDT/TYPE7
Block Format Keyword Describes the interface TYPE7 with constant minimum time step. Which means same behavior as interface TYPE7 with possible switch to Lagrange multiplier formulation, if minimum time step defined with /DT/INTER/CST is reached.
Description
 Same limitation as interface TYPE7 with Lagrange Multiplier formulation.
 Friction is not working after switching into Lagrange Multiplier formulation.
 Not yet compatible with SPMD.
Format
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

/INTER/LAGDT/TYPE7/inter_ID/unit_ID  
inter_title  
grnd_ID_{s}  surf_ID_{m}  I_{stf}  I_{gap}  I_{bag}  I_{del}  
Fscale_{gap}  Gap_{max}  
St_{min}  St_{max}  
Stfac  Fric  Gap_{min}  T_{start}  T_{stop}  
I_{BC}  Inacti  VIS_{s}  VIS_{F}  Bumult  
I_{fric}  I_{filtr}  X_{freq}  I_{form} 
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

C_{1}  C_{2}  C_{3}  C_{4}  C_{5} 
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

C_{6} 
Definitions
Field  Contents  SI Unit Example 

inter_ID  Interface
identifier (Integer, maximum 10 digits) 

unit_ID  Unit Identifier (Integer, maximum 10 digits) 

inter_title  Interface
title (Character, maximum 100 characters) 

grnd_ID_{s}  Secondary nodes group
identifier (Integer) 

surf_ID_{m}  Main surface
identifier (Integer) 

I_{stf}  Stiffness definition flag.
4
(Integer) 

I_{gap}  Gap/element option flag.
(Integer) 

I_{bag}  Airbag vent holes closure
flag in case of contact.
(Integer) 

I_{del}  Node and segment deletion
flag. 2
(Integer) 

Fscale_{gap}  Gap scale
factor. Default = 1.0 (Real) 

Gap_{max}  Maximum gap.
(Real) 
$\left[\text{m}\right]$ 
St_{min}  Minimum
stiffness. (Real) 
$\left[\frac{\text{N}}{\text{m}}\right]$ 
St_{max}  Maximum
stiffness. Default = 10^{30} (Real) 
$\left[\frac{\text{N}}{\text{m}}\right]$ 
Stfac  Interface stiffness, if
I_{stf} =
1 Default = 0 (Real) 
$\left[\frac{\text{N}}{\text{m}}\right]$ 
Stiffness scale factor for
the interface, if I_{stf} =
0 Default = 1.0 (Real) 

Fric  Coulomb
friction (Real) 

Gap_{min}  Minimum gap for impact
activation (Real) 
$\left[\text{m}\right]$ 
T_{start}  Start
time (Real) 
$\left[\text{t}\right]$ 
T_{stop}  Time for temporary
deactivation (Real) 
$\left[\text{t}\right]$ 
I_{BC}  Deactivation flag of
boundary conditions at impact (Boolean) 

Inacti  Deactivation flag of
stiffness in case of initial penetrations. 8
(Integer) 

VIS_{s}  Critical damping
coefficient on interface stiffness. Default set to 0.05 (Real) 

VIS_{F}  Critical damping
coefficient on interface friction. 16
Default set to 1.0 (Real) 

Bumult  Sorting factor. 11
12
Default set to 0.20 (Real) 

I_{fric}  Friction formulation flag.
15
(Integer) 

I_{filtr}  Friction filtering flag.
16
(Integer) 

X_{freq}  Filtering
coefficient. Should have a value between 0 and 1. (Real) 

I_{form}  Friction penalty
formulation type.
(Integer) 

C_{1}  Friction law
coefficient (Real) 

C_{2}  Friction law
coefficient (Real) 

C_{3}  Friction law
coefficient (Real) 

C_{4}  Friction law
coefficient (Real) 

C_{5}  Friction law
coefficient (Real) 

C_{6}  Friction law
coefficient (Real) 
Flags for Deactivation of Boundary Conditions: IBC
(1)1  (1)2  (1)3  (1)4  (1)5  (1)6  (1)7  (1)8 

I_{BCX}  I_{BCY}  I_{BCZ} 
Definitions
Field  Contents  SI Unit Example 

I_{BCX}  Deactivation flag of X
boundary condition at impact
(Boolean) 

I_{BCY}  Deactivation flag of Y
boundary condition at impact
(Boolean) 

I_{BCZ}  Deactivation flag of Z
boundary condition at impact
(Boolean) 
Comments
 For the flag I_{bag}, refer to the monitored volume option (Monitored Volumes (Airbags)).
 Flag I_{del} = 1 has a CPU cost higher than I_{del} = 2.
 If I_{gap} = 2, the
variable gap is computed as:
(1) $$\text{max}\left\{{\mathit{Gap}}_{\text{min}},\text{min}\left[{\mathit{Fscale}}_{\mathit{gap}}\cdot \left({g}_{s}+{g}_{m}\right),{\mathit{Gap}}_{\text{max}}\right]\right\}$$The values given in Line 4 are ignored if I_{gap} ≠ 2.
 Contact stiffness computed
as:
 For I_{stf} = 0,
stiffness
$K={K}_{m}$
 For I_{stf} > 1,
stiffness
(2) $$K=\text{max}\left[S{t}_{\text{min}},\text{min}\left(S{t}_{\text{max}},{K}_{n}\right)\right]$$Where,
${K}_{n}$
is computed from both main
segment stiffness
${K}_{m}$
and secondary node stiffness
${K}_{s}$
I_{stf} = 2, ${K}_{n}=\frac{{K}_{m}+{K}_{s}}{2}$
I_{stf} = 3, ${K}_{n}=\text{max}\left({K}_{m},{K}_{s}\right)$
I_{stf} = 4, ${K}_{n}=\text{min}\left({K}_{m},{K}_{s}\right)$
I_{stf} = 5, ${K}_{n}=\frac{{K}_{m}\cdot {K}_{s}}{{K}_{m}+{K}_{s}}$

${K}_{m}$ is main segment stiffness and computed as:
When main segment lies on a shell or is shared by shell and solid(3) $${K}_{m}=\mathit{Stfac}\cdot 0.5\cdot E\cdot t$$When main segment lies on a solid:(4) $${K}_{m}=B\cdot \frac{{S}^{2}}{V}$$Where, $S$
 Segment area
 $V$
 Volume of the solid
 $B$
 Bulk modulus
 ${K}_{s}$
 An equivalent nodal stiffness considered for
interface TYPE7, and computed as:
 When node is connected to a shell
element:
(5) $${K}_{s}=\mathit{Stfac}\cdot \frac{1}{2}\cdot E\cdot t$$  When node is connected to a solid
element:
(6) $${K}_{s}=\mathit{Stfac}\cdot B\cdot \sqrt[3]{V}$$
 When node is connected to a shell
element:

${K}_{n}$
is computed from both main
segment stiffness
${K}_{m}$
and secondary node stiffness
${K}_{s}$
There is no limitation to the value of stiffness factor (but a value larger than 1.0 can reduce the initial time step).
 For I_{stf} = 0,
stiffness
 The values given in Line 5 are ignored if I_{stf} < 1.
 A default value for Gap_{min} is computed as the minimum
of:
(7) $$Ga{p}_{\text{min}}=\text{min}\left(t,\frac{l}{10},\frac{{l}_{\text{min}}}{2}\right)$$Where, $t$
 Average thickness of the main shell elements
 $l$
 Average side length of the main brick elements
 ${l}_{\mathrm{min}}$
 The smallest side length of all main segments (shell or brick)
 The gap is computed for each
impact as:
(8) $${\mathit{Fscale}}_{gap}\cdot \left({g}_{s}+{g}_{m}\right)$$Where,
${g}_{m}$
: main element gap:
(9) $${g}_{m}=\frac{t}{2}$$with $t$ : thickness of the main element for shell elements
${g}_{m}$ = 0 for brick elements

${g}_{s}$
: secondary node gap:${g}_{s}$ = 0 if the secondary node is not connected to any element or is only connected to brick or spring elements.
(10) $${g}_{s}=\frac{t}{2}$$With $t$ being the largest thickness of the shell elements connected to the secondary node.
${g}_{s}=\frac{1}{2}\sqrt{S}$ for truss and beam elements, with $S$ being the cross section of the element.
If the secondary node is connected to multiple shells and/or beams or trusses, the largest computed secondary gap is used.
The variable gap is always at least equal to Gap_{min}.

${g}_{m}$
: main element gap:
 Deactivation of the boundary condition is applied to secondary nodes group (grnd_ID_{s})
 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.  The sorting factor, Bumult is used to speed up the sorting algorithm.
 The default value for Bumult is automatically increased to 0.30 for models which have more than 1.5 million nodes and to 0.40 for models with more than 2.5 million of nodes.
 One node can belong to the two surfaces at the same time.
 There is no limitation value to the stiffness factor (but a value larger than 1.0 can reduce the initial time step).
 For Friction Formulation
 If the friction flag I_{fric} = 0
(default), the old static friction formulation is used:
${F}_{t}\le \mu \cdot {F}_{n}$ with $\mu =\mathit{Fric}$ ( $\mu $ is Coulomb Friction coefficient)
 For flag I_{fric} > 0,
new friction models are introduced. In this case, the friction coefficient
is set by a function
$\mu =\text{\mu}(\rho ,V)$
Where,
 $\rho $
 Pressure of the normal force on the main segment
 $V$
 Tangential velocity of the secondary node
 If the friction flag I_{fric} = 0
(default), the old static friction formulation is used:
 Currently, the coefficients
C_{1} through C_{6} are used to define a variable
friction coefficient
$\mu $
for new friction formulations.The following formulations are available:
 I_{fric} =
1 (Generalized viscous friction
law):
(11) $$\mu =\mathit{Fric}+{C}_{1}\cdot p+{C}_{2}\cdot V+{C}_{3}\cdot p\cdot V+{C}_{4}\cdot {p}^{2}+{C}_{5}\cdot {V}^{2}$$  I_{fric} =
2 (Darmstad law):
(12) $$\mu ={C}_{1}\cdot {e}^{({C}_{2}V)}\cdot {p}^{2}+{C}_{3}\cdot {e}^{\left({C}_{4}V\right)}\cdot p+{C}_{5}\cdot {e}^{\left({C}_{6}V\right)}$$  I_{fric} =
3 (Renard law):
(13) $$\mu ={C}_{1}+\left({C}_{3}{C}_{1}\right)\cdot \frac{V}{{C}_{5}}\cdot \left(2\frac{V}{{C}_{5}}\right)\hspace{0.17em}\text{if}\hspace{0.17em}V\in \left[0,{C}_{5}\right]$$(14) $$\mu ={C}_{3}\left(\left({C}_{3}{C}_{4}\right)\cdot {\left(\frac{V{C}_{5}}{{C}_{6}{C}_{5}}\right)}^{2}\cdot \left(32\cdot \frac{V{C}_{5}}{{C}_{6}{C}_{5}}\right)\right)\hspace{0.17em}\text{if}\hspace{0.17em}V\in \left[{C}_{5},{C}_{6}\right]$$(15) $$\mu ={C}_{2}\frac{1}{\frac{1}{{C}_{2}{C}_{4}}+{\left(V{C}_{6}\right)}^{2}}\hspace{0.17em}\text{if}\hspace{0.17em}V\ge {C}_{6}$$Where,
${C}_{1}={\mu}_{s}$
${C}_{2}={\mu}_{d}$
${C}_{3}={\mu}_{\mathrm{max}}$
${C}_{4}={\mu}_{\mathrm{min}}$
${C}_{5}={V}_{cr1}$
${C}_{6}={V}_{cr2}$
 First critical velocity ${V}_{cr1}={C}_{5}$ = must be different to 0 (C_{5} ≠ 0).
 First critical velocity ${V}_{cr1}={C}_{5}$ must be lower than the second critical velocity ${V}_{\mathit{cr}2}={C}_{6}({C}_{5}<{C}_{6})$ .
 The static friction coefficient C_{1} and the dynamic friction coefficient C_{2}, must be lower than the maximum friction coefficient C_{2} (C_{4} ≤ C_{1} and C_{4} ≤ C_{2}).
 I_{fric} =
1 (Generalized viscous friction
law):
 Friction FilteringIf I_{filtr} ≠ 0 , the tangential forces are smoothed using a filter:
(16) $${F}_{t}=\alpha \cdot {{F}^{\prime}}_{t}+\left(1\alpha \right)\cdot {{F}^{\prime}}_{t}{}^{1}$$Where α coefficient is calculated from: If I_{filtr} = 1 ➤ $\alpha ={X}_{freq}$ , simple numerical filter
 If I_{filtr} = 2 ➤ $\alpha =\frac{2\cdot \pi}{{X}_{\mathit{freq}}}$ standard 3dB filter, with ${X}_{\mathit{freq}}=\frac{dt}{T}$ and $T$ is filtering period
 If I_{filtr} = 3 ➤ $\alpha =2\cdot \pi \cdot {X}_{\mathit{freq}}\cdot \mathit{dt}$ standard 3dB filter, with X_{freq} is cutting frequency
The filtering coefficient X_{freq} should have a value between 0 and 1.
 Friction penalty formulation
I_{form}
 If I_{form} =
1, (default) viscous formulation, the friction
forces are:
(17) $${F}_{t}=\mathrm{min}\left(\mu {F}_{n},{F}_{\text{adh}}\right)$$  While an adhesion force is computed as:
(18) $${F}_{\mathit{adh}}=C\cdot {V}_{t}\text{\hspace{0.17em}}\text{with}\text{\hspace{0.17em}}\text{\hspace{0.17em}}C={\mathit{VIS}}_{F}\cdot \sqrt{2\mathit{Km}}$$  If I_{form} =
2, stiffness formulation, the friction forces
are:
(19) $${F}_{t}^{\mathit{new}}=\mathrm{min}\left(\mu {F}_{n},{F}_{\mathit{adh}}\right)$$  While an adhesion force is computed as:
(20) $${F}_{\mathit{adh}}={F}_{t}^{\mathit{old}}+\mathrm{\text{\Delta}}{F}_{t}\text{\hspace{0.17em}}\text{with}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{\text{\Delta}}{F}_{t}=K\cdot {V}_{t}\cdot {\delta}_{t}$$Where, ${V}_{t}$ is contact tangential velocity.
 If I_{form} =
1, (default) viscous formulation, the friction
forces are: