/INTER/TYPE21
Block Format Keyword Specific interface between a nondeformable main surface and a secondary surface designed for stamping. All nodes of the main surface must belong to the rigid body.
Description
 A node cannot be a secondary and a main node at the same time.
 The normals to the main segments must be oriented toward the secondary surface.
 For each secondary node, a single impact will be retained, in a way which ensures continuity of the normal force and the tangent force when this impact slides from one segment to a neighboring one.
 Gap may vary according to the variation of shells and 3node shells thickness, on the secondary side.
 Fast search algorithm.
 High speedup with SPMD version.
Format
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

/INTER/TYPE21/inter_ID/unit_ID  
inter_title  
surf_ID_{s}  surf_ID_{m}  I_{stf}  I_{the}  I_{gap}  I_{del}  Invn  I_{adm}  
Fscale_{gap}  Gap_{max}  DEPTH  P_{max}  I_{Tlim}  
St_{min}  St_{max}  P_{skid}  
Stfac  Fric  Gap_{min}  T_{start}  T_{stop}  
I_{BC}  Inacti  VIS_{s}  Bumult  
I_{fric}  I_{filtr}  X_{freq}  sens_ID  fct_ID_{F}  Ascale_{F} 
(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} 
(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

NR_{adm}  P_{adm}  Angl_{adm}  
K_{the}  fct_ID_{K}  Ascale_{K}  T_{int}  I_{the_form}  
F_{rad}  D_{rad}  Fheat  fct_ID_{c}  Dcond  
ID_{rby}  ID_{ref}  Damp  Damp_{r} 
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_ID_{s}  Secondary surface
identifier. (Integer) 

surf_ID_{m}  Main surface
identifier. (Integer) 

I_{stf}  Stiffness definition flag.
(Integer) 

I_{the}  Heat transfer flag.
(Integer) 

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

I_{del}  Secondary node deletion flag.
(Integer) 

Invn  Test flag for inverted
main normals. 21
(Integer) 

I_{adm}  Computing local curvature
flag for adaptive meshing. 3
(Integer) 

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

Gap_{max}  Maximum
gap. (Real) 
$\left[\text{m}\right]$ 
DEPTH  The drawbead depth. 4
(Real) 
$\left[\text{m}\right]$ 
P_{max}  Maximum contact pressure,
due to thickening. 5
Default = 10^{30} (Real) 
$\left[\text{Pa}\right]$ 
I_{Tlim}  Activate tangential force
limitations. 5
(Integer) 

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]$ 
P_{skid}  Maximum contact pressure
used for defining a limit tangential force for skid line output
(/H3D/NODA/SKID_LINE). Default = 10^{30} (Real) 
$\left[\text{Pa}\right]$ 
Stfac  Interface stiffness (if
I_{stf} =
1). Default set to 0.0 (Real) 
$\left[\frac{\text{N}}{\text{m}}\right]$ 
Stiffness scale factor for
the interface (if I_{stf} =
0). Default set to 1.0 (Real) 

Fric  Coulomb friction (if fct_ID_{F} =
0). Default = 0.0 (Real) 

Coulomb friction scale
factor (if fct_ID_{F}
$\ne $
0). Default = 1.0 (Real) 

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

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

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

Bumult  Sorting factor is used to
speed up the sorting algorithm. 10 Default set to 0.20 (Real) 

I_{fric}  Friction formulation flag.
12
13
(Integer) 

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

X_{freq}  Filtering coefficient.
15 (Real) 

sens_ID  Sensor identifier to
activate/deactivate the interface. 20 If an identifier sensor is defined, the activation/deactivation of interface is based on sensor and not on T_{start} or T_{stop}. (Integer) 

fct_ID_{F}  Function identifier for
friction coefficient with temperature. Default = 0 (Integer) 

Ascale_{F}  Abscissa scale factor on
fct_ID_{F}. Default = 1.0 (Real) 
$\left[\text{K}\right]$ 
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) 

NR_{adm}  Number of elements through
a 90° radius (used only if I_{adm} =2). Default = 3 (Integer) 

P_{adm}  Criteria on the percentage
of penetration. Default = 1.0 (Real) 

Angl_{adm}  Angle
criteria. (Real) 
$\left[\mathrm{deg}\right]$ 
K_{the}  Conductive heat exchange
coefficient (if fct_ID_{K} =
0). Default = 0.0 (Real) 
$\left[\frac{\text{W}}{{\text{m}}^{\text{2}}\text{K}}\right]$ 
Heat exchange scale factor
(if fct_ID_{K}≠
0). Default = 0.0 (Real) 

fct_ID_{K}  Function identifier for
heat exchange definition with contact pressure. Default = 0 (Integer) 

Ascale_{K}  Abscissa scale factor on
fct_ID_{K}. Default = 1.0 (Real) 
$\left[\text{Pa}\right]$ 
T_{int}  Interface
temperature. (Real) 
$\left[\text{K}\right]$ 
I_{the_form}  Heat contact formulation flag.
(Integer) 

F_{rad}  Radiation factor. 17 (Real) 
$\left[\frac{\text{W}}{{\text{m}}^{\text{2}}{\text{K}}^{\text{4}}}\right]$ 
D_{rad}  Maximum distance for
radiation computation. (Real) 
$\left[\text{m}\right]$ 
Fheat  Frictional heating factor.
18 (Real) 

fct_ID_{c}  Function identifier for
the conductive heat exchange coefficient definition as a function of
distance. 19 Default = 0 (Integer) 22 

Dcond  Maximum distance for
conductive heat exchange. 23 Default = 0.0 (Real) 
$\left[\text{m}\right]$ 
ID_{rby}  Rigid body
identifier. (Integer) 

ID_{ref}  Reference interface TYPE21
identifier for damping.
(Integer) 

Damp  Translational critical
damping factor. 19
(Real) 

Damp_{r}  Rotational critical
damping factor. (Real) 
Comments
 In case of SPMD, each main segment defined by surf_ID_{m} must be associated to an element (possibly to a void element).
 Contact gap
 If I_{gap} = 0
(constant gap),
Gap is constant over the secondary surface and along the time, equal to Gap_{min}. And a default value for Gap_{min} is computed as $t/2$ , $t$ being the average thickness of the secondary shell elements.
In case of constant gap, Gap_{max} and Fscale_{gap} will not be used.
If contact thickness of the part is not defined in input /PART:
 If I_{gap} = 1,
variable gap over the secondary surface is computed as:
(1) $$\mathrm{max}\left[Ga{p}_{\mathrm{min}},\mathrm{min}\left(\mathit{Fscal}{e}_{gap}\cdot {g}_{s},Ga{p}_{\mathrm{max}}\right)\right]$$  If I_{gap} = 2,
variable gap over the secondary surface and along the time is computed
at each time, as:
(2) $$\mathrm{max}\left[Ga{p}_{\mathrm{min}},\mathrm{min}\left(\mathit{Fscal}{e}_{gap}\cdot {g}_{s},Ga{p}_{\mathrm{max}}\right)\right]$$and will vary along the time according to the variation of shells and 3node shells thickness, on the secondary side.
If contact thickness of the part is defined in input /PART:
 If I_{gap} = 1,
variable gap is computed as:
(3) $$\mathrm{max}\left[Ga{p}_{\mathrm{min}},\mathrm{min}\left(\mathit{Fscal}{e}_{gap}\cdot \frac{{t}_{part}}{2},Ga{p}_{\mathrm{max}}\right)\right]$$  If I_{gap} = 2,
variable gap is computed as:
(4) $$\mathrm{max}\left[Ga{p}_{\mathrm{min}},\mathrm{min}\left(\mathit{Fscal}{e}_{gap}\cdot \frac{t{t}_{0}{t}_{part}}{2},Ga{p}_{\mathrm{max}}\right)\right]$$Where,
${g}_{s}$
: secondary node gap
${g}_{s}=\frac{t}{2}$ , with $t$ is thickness of the secondary element for shell elements.
${g}_{s}=0$ for brick elements.
 ${t}_{part}$ is the contact thickness of the part.
 ${t}_{0}$ is the initial thickness of the secondary node.
If I_{gap} = 1 or 2, the variable gap is always at most equal to Gap_{max} and default value for Gap_{max} will be set to 10^{30} and is always at least equal to Gap_{min} (but there is no default value for Gap_{max}).

${g}_{s}$
: secondary node gap
 If I_{gap} = 0
(constant gap),
 In case of adaptive meshing:
 I_{adm} =
1:If the contact occurs in a zone (main side) whose radius of curvature is lower than the element size (secondary side), the element on the secondary side will be divided (if not yet at maximum level).
 I_{adm} =
2:
If the contact occurs in a zone (main side) whose radius of curvature is lower than NR_{adm} times the element size (secondary side), the element on the secondary side will be divided (if not yet at maximum level).
If the contact occurs in a zone (main side) where the angles between the normals are greater than Angl_{adm}, and the percentage of penetration is greater than P_{adm}, the element on the secondary side will be divided (if not yet at maximum level).
 I_{adm} =
1:
 The interface allows secondary
nodes to cross the main surface; if a secondary node gets into the main surface
from a distance greater than DEPTH, no contact force is computed on the node.
A default value for DEPTH is computed as the maximum of:
 Upper value of the gap (at time 0) among all nodes
 Smallest side length of secondary element
If the input value is not equal to 0, DEPTH will be raised up to the upper value of the gap (at time 0) among all nodes.
Too large of a DEPTH will decrease the performances of search algorithms for contact.
 Maximum contact pressure due to
thickening. P_{max} is used only
if I_{gap} = 2.
 It can be used for limiting the contact force in case of thickening.
 It can be used for limiting the normal contact force in case of
thickening according to the following equation:
(5) $${F}_{n}\le {P}_{\mathrm{max}}\cdot S$$  The tangent contact force is also limited when the flag
I_{Tlim} =
0 by the following equation:
(6) $${F}_{t}\le {P}_{\mathrm{max}}\cdot \frac{S}{\sqrt{3}}$$
 $F=K$
 Stfac can be larger than 1.0.
 Deactivation of the boundary condition is applied to secondary nodes.
 Inacti = 3 may create
initial energy if the node belongs to a spring element.
Inacti = 5 or Inacti = 6, the gap is initially reduced and recovers its computed value as the secondary node depenetrates.
Inacti = 6 is recommended instead of Inacti =5, to avoid high frequency effects into the interface.  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.
 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 $ is Coulomb Friction coefficient.
 If fct_ID_{F} =
0
Fric is Coulomb friction
$\mu =\mathit{Fric}$

If fct_ID_{F} ≠ 0
Fric becomes a scale factor of Coulomb friction coefficient which depends on the temperature.(7) $$\mu =\mathit{Fric}\cdot {\mathrm{f}}_{F}\left(Ascal{e}_{F},{T}_{interface}\right)$$While, ${T}_{interface}$ is the interface temperature which is taken as the mean temperature of secondary and main:(8) $${T}_{interface}=\frac{{T}_{secondary}+{T}_{main}}{2}$$
 If fct_ID_{F} =
0
 For flag I_{fric} > 0,
new friction models are introduced. In this case, the friction
coefficient is set by a function.
(9) $$\mu =\text{\mu}(\rho ,V)$$Where, $p$
 Pressure of the normal force on the main segment
 $V$
 Tangential velocity of the secondary node relative to the main segment
 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 (if I_{fric} = 0) and C_{6} is not used (if I_{fric} = 1) to define
a variable friction coefficient
$\mu $
for new friction formulations.The following formulations are available:
 I_{fric} = 1
(Generalized Viscous Friction Law):
(10) $$\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
(Modified Darmstad Law):
(11) $$\mu =Fric+{C}_{1}\cdot {e}^{\left({C}_{2}V\right)}\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):
$\mu ={C}_{1}+\left({C}_{3}{C}_{1}\right)\cdot \frac{V}{{C}_{5}}\cdot \left(2\frac{V}{{C}_{5}}\right)$ if $V\in \left[0,{C}_{5}\right]$
$\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)$ if $V\in \left[{C}_{5},{C}_{6}\right]$
$\mu ={C}_{2}\frac{1}{\frac{1}{{C}_{2}{C}_{4}}+{\left(V{C}_{6}\right)}^{2}}$ if $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}\ne 0$ ).
 First critical velocity ${V}_{cr1}={C}_{5}$ must be lower than the second critical velocity ${V}_{cr2}={C}_{6}$ ( ${C}_{5}<{C}_{6}$ ).
 The static friction coefficient ${C}_{1}$ and the dynamic friction coefficient ${C}_{2}$ , must be less than the maximum friction ${C}_{3}$ ( ${C}_{1}\le {C}_{3}$ and ${C}_{2}\le {C}_{3}$ ).
 The minimum friction coefficient ${C}_{4}$ must be less than the static friction coefficient ${C}_{1}$ and the dynamic friction coefficient ${C}_{2}$ ( ${C}_{4}\le {C}_{1}$ and ${C}_{4}\le {C}_{2}$ ).
Table 1. Units for Friction Formulations I_{fric} Fric C_{1} C_{2} C_{3} C_{4} C_{5} C_{6} 1 $\left[\frac{1}{\text{P}\text{a}}\right]$ $\left[\frac{\text{s}}{\text{m}}\right]$ $\left[\frac{\text{s}}{\text{Pa}\cdot \text{m}}\right]$ $\left[\frac{1}{{\text{Pa}}^{2}}\right]$ $\left[\frac{{\text{s}}^{2}}{{\text{m}}^{2}}\right]$ 2 $\left[\frac{1}{{\text{Pa}}^{2}}\right]$ $\left[\frac{\text{s}}{\text{m}}\right]$ $\left[\frac{1}{\text{P}\text{a}}\right]$ $\left[\frac{\text{s}}{\text{m}}\right]$ $\left[\frac{\text{s}}{\text{m}}\right]$ 3 $\left[\frac{\text{m}}{\text{s}}\right]$ $\left[\frac{\text{m}}{\text{s}}\right]$  I_{fric} = 1
(Generalized Viscous Friction Law):
 The formulation for friction is a
stiffness (incremental) formulation, and the friction forces
are:
(12) $${F}_{t}^{new}=\mathrm{min}\left(\mu {F}_{n},{F}_{adh}\right)$$While and adhesion force is computed as:
${F}_{adh}={F}_{t}^{old}+\mathrm{\text{\Delta}}{F}_{t}$ with $\mathrm{\text{\Delta}}{F}_{t}=K\cdot {V}_{t}\cdot dt$
Where, ${V}_{t}$ is the tangential relative velocity of the secondary node with the main segment.
 Friction filtering:If I_{filtr} ≠ 0 , the tangential forces are smoothed using a filter:
(13) $${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}_{freq}}$ , standard 3dB filter, with ${X}_{freq}=\frac{dt}{T}$ , and $T$ is the filtering period
 If I_{filtr} = 3: $\alpha =2\cdot \pi \cdot {X}_{freq}\cdot dt$ , standard 3dB filter, with X_{freq} is cutting frequency
The filtering coefficient X_{freq} should have a value between 0 and 1.
 Heat exchangeBy I_{the}= 1 (heat transfer activated) to consider heat exchange and heat friction in contact.
 If I_{the_form} =0, then heat exchange is between shell and constant temperature contact T_{int}.
 If I_{the_form} =1, then heat exchange is between all contact pieces.
T_{int} is used only when I_{the_form}= 0. In this case. The temperature of main side assumed to be constant (equal to T_{int}). If I_{the_form}= 1, then T_{int} is not taken into account. So the nodal temperature of main side will be considered.
If I_{the} = 2, Heat transfer is computed using thermal conductance K_{the} only for the secondary side. The temperature of the main side is not assumed to be a constant but is calculated from the temperature field defined on each main node. These nodal temperatures can vary over time and space which are defined using /IMPTEMP.
Thermal conduction
I_{the} = 1 needs the material of the secondary side to be a thermal material using finite element formulation for heat transfer (/HEAT/MAT).
Thermal conduction is computed when the secondary node falls into gap:(14) $$gap=\mathrm{max}\left[Ga{p}_{\mathrm{min}},\mathrm{min}\left(\mathit{Fscal}{e}_{gap}\cdot {g}_{s},Ga{p}_{\mathrm{max}}\right)\right]$$Heat exchange coefficient If fct_ID_{K} = 0, then K_{the} is heat exchange coefficient and heat exchange depends only on heat exchange surface.
 If fct_ID_{K} ≠
0, K_{the} is a scale factor
and heat exchange depends on contact pressure:
(15) $$\mathrm{K}={K}_{the}\cdot {\mathrm{f}}_{K}\left(Ascal{e}_{K},P\right)$$While ${\mathrm{f}}_{K}$ is the function of fct_ID_{K}.
 Radiation:Radiation is considered in contact if ${F}_{rad}\ne 0$ and the distance, ${d}_{}$ , of the secondary node to the main segment is:
(16) $$\mathit{Gap}<d<{D}_{\mathit{rad}}$$While ${D}_{rad}$ is the maximum distance for radiation computation. The default value for ${D}_{rad}$ is computed as the maximum of: Upper value of the Gap (at time 0) among all nodes
 Smallest side length of secondary element
It is recommended not to set the value too high for ${D}_{rad}$ , which may reduce the performance of Radioss Engine.
A radiant heat transfer conductance is computed as:(17) $${h}_{\mathit{rad}}={F}_{\mathit{rad}}\left({{T}_{m}}^{2}+{{T}_{s}}^{2}\right)\cdot \left({T}_{m}+{T}_{s}\right)$$with(18) $${F}_{\mathit{rad}}=\frac{\sigma}{\frac{1}{{\epsilon}_{1}}+\frac{1}{{\epsilon}_{2}}1}$$Where, $\sigma =5.669\times {10}^{8}\left[\frac{\text{W}}{{\text{m}}^{2}{\text{K}}^{4}}\right]$
 Stefan Boltzman constant
 ${\epsilon}_{1}$
 Emissivity of secondary surface
 ${\epsilon}_{2}$
 Emissivity of main surface
 Heat Friction
 Frictional energy is converted into heat when I_{the} > 0 for interface
 Fheat is defined as the fraction of this energy which is converted into heat and transferred to the secondary side.
 The frictional heat Q_{Fric}
is so defined for a stiffness formulation:
(19) $${Q}_{\mathit{Fric}}=Fheat\cdot \frac{\left({F}_{adh}{F}_{t}\right)}{K}\cdot {F}_{t}$$
 Critical damping factors allow
for the reduction of dynamic effects, especially for those tools where a loading
is applied. This can be used to model the hydraulic press system:
A damping force (resp. torque) is applied to the tool:
${F}_{d}=Cv$ (resp.) ${m}_{d}={C}_{r}{v}_{r}$
With $C=Damp\cdot \sqrt{2\cdot mass\cdot K}$ resp. ${C}_{r}=Dam{p}_{r}\cdot \sqrt{2\cdot I\cdot {K}_{r}}$
Where,
$C$ is a percentage of the critical damping with the tool mass mass (resp. inertia $I$ ).
$K$ is the total interface stiffness (resp. rotational stiffness ${K}_{r}$ ).
$v$ (resp. ${v}_{r}$ ) is the translational (resp. rotational) velocity of the tool, if ID_{ref} is equal to 0; otherwise, it is the relative velocity with respect to tool of interface ID_{ref}.
 When sens_ID is defined for activation/deactivation of the interface, T_{start} and T_{stop} are not taken into account.
 If Invn =
1, Radioss will check if main
normal are oriented towards the blank or not. Not well oriented normal may cause
big penetrations and wrong results. If this happens, it is important to stop
computation when this problem is detected. This option used to check if
secondary node is contacting the main segment from the correct side (Figure 6) at the time of first impact.
If not, the computation is stopped with an error message.
 Main normal is well oriented: First impact is from correct side (penetration is smaller than Gap).
 Main normal is not well oriented: First impact is from wrong side (penetration is almost equal to DEPTH).
This option is only available for shells and must be used with a special care, because the computation can be stopped even if the normals are well oriented when there are large initial penetrations.
 When
fct_ID_{c} ≠0, the heat
transfer coefficient can change as function of distance d when
Gap < d ≤ Dcond.
Abscises and ordinates of this function must be between 0 and 1.
The heat transfer coefficient is computed as:(20) $$K={K}_{the}\left(P=0\right)*fct\text{}\_I{D}_{c}\left(\frac{dGap}{DcondGap}\right)$$The maximum value of $K$ is equal to the value of K_{the} when K_{the} is constant. Otherwise, in case of K_{the} depending on pressure, the maximum is equal to value of K_{the} for contact pressure $P$ =0.
$K$ drops to zero when distance is equal to Dcond.
 When
fct_ID_{c}≠0, the heat
transfer is computed as:
 Conductive heat transfer when d < Gap
 Conductive and radiative heat transfer when Gap < d ≤ Dcond.
 Radiative heat transfer when Dcond < d ≤ Drad
 When
fct_ID_{c} ≠ 0 and
Dcond = 0, then
Dcond=Drad.
When Frad ≠ 0, fct_ID_{c} ≠ 0, and Drad = 0, then Drad = Dcond.