/INTER/TYPE25
Block Format Keyword TYPE25 is a general nodes to surface contact interface using the penalty method. The penalty stiffness is constant and therefore the time step is not affected for standard penalty stiffness.
Unlike TYPE7 contact, solid elements have zero contact gap thickness. Three types of contact inputs can be defined: single surface, surface to surface, or nodes to surface. This contact interface can replace interface TYPE3, TYPE5, TYPE7, or TYPE24. This interface is not available with the implicit solution.
Format
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
/INTER/TYPE25/inter_ID/unit_ID | |||||||||
inter_title | |||||||||
surf_ID1 | surf_ID2 | Istf | Igap | Irem_i2 | Idel | ||||
grnd_IDs | Gap_scale | %mesh_size | Gap_max_s | Gap_max_m | |||||
Stmin | Stmax | Igap0 | Ishape |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
Stfac | Fric | Tstart | Tstop | ||||||
IBC | IVIS2 | Inacti | VISs | ||||||
Ifric | Ifiltr | Xfreq | sens_ID | fric_ID |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
C1 | C2 | C3 | C4 | C5 |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
C6 |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
ViscFluid | SigMaxAdh | ViscAdhFact |
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) |
|
surf_ID1 | First surface identifier. 1 (Integer) |
|
surf_ID2 | Second surface
identifier (Integer) |
|
Istf | Interface stiffness
definition flag. 2
(Integer) |
|
Igap | Gap/element option flag.
3
|
|
Irem_i2 | Deactivating flag for the
slave node, if the same contact pair (nodes) has been defined in
interface TYPE2.
|
|
Idel | Node and segment deletion flag.
|
|
grnd_IDs | Nodes group
identifier
1 If defined, node group will be added as slave nodes. (Integer) |
|
Gap_scale | Gap scale factor for all
Igap options. Default = 1.0 (Real) |
|
%mesh_size | Percentage of mesh size
(used only when Igap = 3). Default = 0.4 (Real) |
|
Gap_max_s | Slave maximum gaps. 3 Default = 1030 (Real) |
|
Gap_max_m | Master maximum gaps. 3 Default = 1030 (Real) |
|
Stmin | Minimum stiffness (used
only when Istf >
1 and Istf < 7) 2 (Real) |
|
Stmax | Maximum stiffness (used
only when Istf > 1
and Istf < 7) 2 Default = 1030 (Real) |
|
Igap0 | Gap modification flag for
slave shell nodes on the free edges. 3
(Integer) |
|
Ishape | Flag defining the shape of
the gap along free edges.
(Integer) |
|
Stfac | Interface stiffness scale
factor. 2 Default = 1.0 (Real) |
|
Fric | Coulomb friction
(Real) |
|
Tstart | Start time. 8
(Real) |
|
Tstop | Temporary deactivation
time. 8 Default = 1030 (Real) |
|
IBC | Deactivation flag of
boundary conditions at impact. (Boolean) |
|
Inacti | Initial penetration flag.
(Integer) |
|
VISs | Critical damping
coefficient on interface stiffness. Default = 0.05 (Real) |
|
Ifric | Friction formulation
flag. Only used if fric_ID is not defined.
(Integer) |
|
Ifiltr | Friction filtering flag.
(Integer) |
|
Xfreq | Filtering
coefficient. Default = 1.0 (Real) |
|
sens_ID | Sensor identifier to
activate/deactivate the interface. (Integer) |
|
fric_ID | Friction identifier for
friction definition for selected pairs of parts.
(Integer) |
|
C1 | Friction law coefficient.
5 (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) |
|
IVIS2 | Interface adhesion flag.
10
(Interger) |
|
ViscFluid | Viscosity of the fluid at
the interface. 10 (Real) |
|
SigMaxAdh | Maximum transverse
adhesive stress at interface. 10 (Real) |
|
ViscAdhFact | Tangential viscous
resistant force scaling factor. 10 (Real) |
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 |
Definitions
Field | Contents | SI Unit Example |
---|---|---|
IBCX | Deactivation flag of X
boundary condition at impact.
(Boolean) |
|
IBCY | Deactivation flag of Y
boundary condition at impact.
(Boolean) |
|
IBCZ | Deactivation flag of Z
boundary condition at impact.
(Boolean) |
Comments
- Contact master/slave pairs
can be defined in three ways:
- Single self-impacting surface only: surf_ID1 > 0, and surf_ID2 = 0
- Symmetric surface to surface: surf_ID1 > 0, and surf_ID2 > 0
- Nodes to surface: grnd_IDs > 0, surf_ID1 = 0, and surf_ID2 > 0
grnd_IDs > 0 is used to define node to surface contact type, but it may also be used in other contact types. In that case, the node group will be added simply as supplementary slave nodes, which is useful when users want to add spring element nodes, master node of rigid body, etc. into the contact (as slave nodes).
If the surface is defined with shells, two contact segments (shifted by half thickness (t)) with opposite normal directions will be generated:In case of SPMD, each master segment defined by surf_IDi (i=1, 2) must be associated to an element (possibly to a void element).
In cases where quadratic elements are used, it is recommended to define the surfaces by using /SURF/PART/EXT as in that case, middle nodes of quadratic elements are used in the contact treatment.
The surface definition /SURF/PART/ALL is not available with TYPE25.
- Contact stiffness,
is computed as:
(1) Where, depends on Istf:- Istf = 1000,
- Istf = 2,
- Istf = 3,
- Istf = 4,
- Istf = 5,
: master segment stiffness and computed as:- , when the master segment lies on a shell.
- , when master segment lies on a solid.
- , when master segment is shared by shell and solid.
: Slave node stiffness is an equivalent nodal stiffness considered for interface TYPE25, and computed as:- , when node is connected to a shell element,
- , when node is connected to solid element.
Where,- Segment area
- Volume of the solid
- Bulk modulus
- Thickness of the shell
The Stfac value can be larger than 1.0. There is no limitation value to the stiffness factor (a value larger than 1.0 can reduce the initial time step).
When using /PROP/VOID and /MAT/VOID, material properties and thickness for the VOID material must be entered; otherwise, the contact stiffness of the void elements will be zero. This is especially important if VOID shell elements share elements with solid elements as the stiffness of the shell elements is used in the contact calculation.
- The gap is computed
automatically for each impact as:
- If Igap = 1, variable gap is computed as:
(2) - If Igap=2, variable gap is computed as:
(3) with deactivation of slave nodes when the element size is smaller than gap values:For self-impact contact, when Curvilinear Distance (from a node of the master segment to a slave node) is smaller than (in initial configuration), this slave node will not be taken into account by this master segment, and it will not be deleted from the contact for the other master segments.
- If Igap= 3, variable gap is computed as:
(4) Where,-
: master element gap:
, with is the thickness of the master element for shell elements
, for brick elements
-
: slave node gap:
, if the slave node is not connected to any element or is only connected to brick or spring elements.
, if the slave node is connected to a shell element, with being the largest thickness of the shell elements connected to the slave node.
, if the slave node is connected to truss or beam elements, with being the cross section of the 1D element.
If the gap modification flag for slave shell nodes on the free edges Igap0 is set to 1: is reset to zero if the slave node lies on the free edges of the slave surface. The gap modification flag for slave shell nodes on the free edges has no effect if the slave node is defined through the optional node group (grnod_IDs).
If the slave node is connected to multiple shells and/or beams or trusses, the largest computed slave gap is used.
-
: master element gap:
- : length of the smallest edge of the master segment.
-
: if the slave node belongs to the master
surface,
is the length of the smalledst edge of
master segments connected to the slave node,
=1E+30, otherwise.
In any case, amd are limited separately by Gap_max_m and Gap_max_s before the gap is computed.
If the slave node does not belong to the master surface, the gap remains(5)
- If Igap = 1, variable gap is computed as:
- Ishape
determines the shape of the gap along shell free edges on the master side. The
gap never extends more than the slave node gap out of the free edge. But
Ishape determines if the
shape of this gap is square or round and the contact force (normal)
direction.The gap used for contact at the master free edges and resulting force direction are as:
Ishape =1 is not available with Igap =3 and will then be reset to Ishape =2.
- If fric_ID is defined, the contact friction is defined in /FRICTION
and the friction inputs (Ifric, C1, etc.) in this input card are not used.The friction forces are:
(6) While an adhesion force is computed as:
with
Where, is the Coulomb friction coefficient and is defined as:- For flag Ifric by default:
with (Coulomb friction)
- For flag Ifric > 1, new friction
models are introduced. In this case, the friction coefficient is set
by a function:Where,
- Pressure of the normal force on the master segment
- Tangential velocity of the slave node relative to the master segment
Currently, the coefficients C1 through C6 are used to define a variable friction coefficient for new friction formulations.
The following formulations are available:- Ifric = 1
(Generalized Viscous Friction law):
(7) - Ifric = 2
(Modified Darmstad law):
(8) - Ifric = 3
(Renard law):
if
if
if
Where,- First critical velocity must be different to 0 ( ).
- First critical velocity must be less than the second critical velocity .
- The static friction coefficient and the dynamic friction coefficient , must be less than the maximum friction ( and ).
- The minimum friction coefficient
must be less than the
static friction coefficient
and the dynamic friction
coefficient
(
and
).
Table 1. Units for Friction Formulations Ifric Fric C1 C2 C3 C4 C5 C6 1 2 3
- For flag Ifric by default:
- Friction filteringIf Ifiltr = 1, 2 or 3, the tangential forces are smoothed using a filter:
(9) Where, α coefficient is calculated from:- If Ifiltr = 1: α = Xfreq, simple numerical filter
- If Ifiltr = 2: , standard -3dB filter, with , and T = filtering period
- If Ifiltr = 3: , standard -3dB filter, with Xfreq = cutting frequency
The filtering coefficient Xfreq should have a value between 0 and 1.
- Inacti and Ipen_max,
initial penetration treatment:
- Inacti = 1000: The initial penetrations are ignored: no contact force is applied, but the nodes are not deactivated from the contact; if the node goes out of the contact and later gets back into contact, contact forces are then applied.
- Inacti = -1: Initial forces are applied on all penetrating nodes. High initial penetrations should be avoided, as they might generate high contact forces and lead to high energy error at the beginning of the computation.
- Inacti = 5: The master segment is shifted by the initial penetration value ( ); therefore, at time zero no initial forces are applied.
The master segment position is restored only in case of rebound larger than .
In the opposite case, when slave node continues to penetrate, the penetration is computed as:(10) - Intersections and large initial penetration (Inacti= -1 and
5):
Shells: initial intersections should be avoided, as they will lead to wrong direction of contact force and possible slave nodes anchorage.
- When sens_ID is defined for activation/deactivation of the interface,Tstart and Tstop are not taken into account.
- For output forces:
When the contact type is asymmetric surface to surface, the output normal contact forces in Time History are calculated correctly, if the two surfaces are well separated.
- IVIS2=-1: is used to add
adhesion in the normal direction and viscous resistive forces in the tangential
direction. This can be used to model thermoplastic composite forming.When used, half of the contact gap is considered an adhesive zone and the other half a physical contact zone. Therefore, to maintain the same physical contact gap, the contact thickness should be doubled using Gap_scale.The adhesive force is only applied after slave nodes have entered the physical contact zone and then move back into the adhesion zone. The adhesive force acts to prevent the node from moving out of the adhesion zone and is applied in the normal direction.
(11) Where,- Area of the slave surface
- Penetration into the adhesion zone
- Contact gap as calculated in Comment 3
The adhesive spring ruptures as the node exits the adhesion zone and will be recreated if the node enters the contact zone again.
Viscous resistive forces are applied in the tangential direction when the slave nodes enter into the adhesion zone. A viscous tangential opposing force is applied instead of a friction force and is calculated as:(12) Where,- Area of the slave surface
- Penetration into the adhesion zone
- Contact gap