Contact Interfaces

WARNING ID 94 and WARNING ID 477.

This Starter output warning occurs in contact interfaces that use the constant gap formulation, Igap=1000 and Gapmin is defined such that (1)
G a p min > l min 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4raiaadg gacaWGWbWaaSbaaSqaaiGac2gacaGGPbGaaiOBaaqabaGccqGH+aGp daWcaaqaaiaadYgadaWgaaWcbaGaciyBaiaacMgacaGGUbaabeaaaO qaaiaaikdaaaaaaa@4173@

Where, l min MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaaciGGTbGaaiyAaiaac6gaaeqaaaaa@39E6@ is the minimum main element edge length.

WARNING ID :           94
** WARNING IN INTERFACE GAP
DESCRIPTION :  
   -- INTERFACE ID : x
   -- INTERFACE TITLE : xxx
   INPUT GAPmin xxx > (MINIMUM MAIN ELEMENT EDGE LENGTH) / 2
   FOR SELF CONTACT, ARTIFICIAL STIFFENING CAN OCCUR IF AN ELEMENT IS COMPRESSED
   TO PREVENT THIS WARNING, DEFINE GAPmin < xxx
  POSSIBLE IMPACT NUMBER:        xxx  (<=           2 *NSN)

The Gapmin: will be calculated, if it is not defined in the model.

This warning is important for self-impacting interfaces and can be ignored if the interface is not self-imp-impacting. Figure 1 shows an element with edge length, L MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamitaaaa@36C7@ , and a self-impacting interface where the G a p = L 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4raiaadg gacaWGWbGaeyypa0ZaaSaaaeaacaWGmbaabaGaaGOmaaaaaaa@3B41@ . If this element is compressed more than 50%, the red node enters the gap of the neighboring element and a self-contact is then detected which is not physical. This leads to over-stiffening of the structure and non-physical behavior.


Figure 1.

Therefore, it is recommended to define Ga p min < l min 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4raiaadg gacaWGWbWaaSbaaSqaaiGac2gacaGGPbGaaiOBaaqabaGccqGH8aap daWcaaqaaiaadYgadaWgaaWcbaGaciyBaiaacMgacaGGUbaabeaaaO qaaiaaikdaaaaaaa@416F@ .

If G a p min > L MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4raiaadg gacaWGWbWaaSbaaSqaaiGac2gacaGGPbGaaiOBaaqabaGccqGH+aGp caWGmbaaaa@3D7F@ , a self-contact will be detected which generates contact forces from the beginning of the computation. Such a situation can be accepted if it is local enough, but not if it is a frequent situation over the self-impacting interface. It is possible to find the source of this message by using a pre-processor and selecting the elements through size criteria.

If variable gap Igap=1 or 2 is used, then the following warning can occur:
WARNING ID :          477
** WARNING IN INTERFACE GAP
DESCRIPTION :  
   -- INTERFACE ID : x
   -- INTERFACE TITLE : xxx
   MAXIMUM VARIABLE GAP xxx > (MINIMUM MAIN ELEMENT EDGE LENGTH)/2
   FOR SELF CONTACTS WHERE THE VARIABLE GAP IS LARGE AND THE ELEMENT EDGE LENGTH IS SMALL
   ARTIFICIAL STIFFENING CAN OCCUR DURING ELEMENT COMPRESSION
SOLUTION :     
TO PREVENT THIS ISSUE, USE Igap=3 TO AUTOMATICALLY REDUCE THE GAP BASED ON ELEMENT SIZE
  POSSIBLE IMPACT NUMBER:        xxx  (<=           1 *NSN)

This warning message is like the constant contact gap warning except in this case the maximum variable gap for the entire contact is used in the check, max ( G a p v a r i a b l e ) > l min 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaciyBaiaacg gacaGG4bGaaiikaiaadEeacaWGHbGaamiCamaaBaaaleaacaWG2bGa amyyaiaadkhacaWGPbGaamyyaiaadkgacaWGSbGaamyzaaqabaGcca GGPaGaeyOpa4ZaaSaaaeaacaWGSbWaaSbaaSqaaiGac2gacaGGPbGa aiOBaaqabaaakeaacaaIYaaaaaaa@4A3B@ . This warning is not as precise because it is very possible that the maximum variable contact Gap is in a different location then the minimum element length and therefore no problem exists in the model.

However as suggested in the warning message, using the variable contact gap Igap=3 will automatically reduce the contact gap in any area where the variable contact gap, calculated as F s c a l e g a p ( g s + g m ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGgbGaam 4CaiaadogacaWGHbGaamiBaiaadwgadaWgaaWcbaGaam4zaiaadgga caWGWbaabeaakiabgwSixpaabmaabaGaam4zamaaBaaaleaacaWGZb aabeaakiabgUcaRiaadEgadaWgaaWcbaGaamyBaaqabaaakiaawIca caGLPaaaaaa@47AA@ , is large as compared to the element edge length by adding the following term to the variable gap calculation.(2)
% m e s h _ s i z e ( g s _ l + g m _ l ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaGGLaGaam yBaiaadwgacaWGZbGaamiAaiaac+facaWGZbGaamyAaiaadQhacaWG LbGaeyyXIC9aaeWaaeaacaWGNbWaaSbaaSqaaiaadohacaGGFbGaam iBaaqabaGccqGHRaWkcaWGNbWaaSbaaSqaaiaad2gacaGGFbGaamiB aaqabaaakiaawIcacaGLPaaaaaa@4C05@
Where,
g m_l MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGNbWaaS baaSqaaiaad2gacaGGFbGaamiBaaqabaaaaa@3A3C@
Length of the smaller edge of element connected to the main
g s_l MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGNbWaaS baaSqaaiaadohacaGGFbGaamiBaaqabaaaaa@3A42@
Length of the smaller edge of element connected to the secondary node
%mesh_size MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaGGLaGaam yBaiaadwgacaWGZbGaamiAaiaac+facaWGZbGaamyAaiaadQhacaWG Lbaaaa@3F7A@
Percentage of mesh size, default=0.4.
So the variable GAP calculation for Igap=3 becomes:(3)
max{ Ga p min ,min[ Fscal e gap ( g s + g m ),%mesh_size( g s_l + g m_l ),Ga p max ] } MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaciGGTbGaai yyaiaacIhadaGadaqaaiaadEeacaWGHbGaamiCamaaBaaaleaaciGG TbGaaiyAaiaac6gaaeqaaOGaaiilaiGac2gacaGGPbGaaiOBamaadm aabaGaamOraiaadohacaWGJbGaamyyaiaadYgacaWGLbWaaSbaaSqa aiaadEgacaWGHbGaamiCaaqabaGccqGHflY1daqadaqaaiaadEgada WgaaWcbaGaam4CaaqabaGccqGHRaWkcaWGNbWaaSbaaSqaaiaad2ga aeqaaaGccaGLOaGaayzkaaGaaiilaiaacwcacaWGTbGaamyzaiaado hacaWGObGaai4xaiaadohacaWGPbGaamOEaiaadwgacqGHflY1daqa daqaaiaadEgadaWgaaWcbaGaam4Caiaac+facaWGSbaabeaakiabgU caRiaadEgadaWgaaWcbaGaamyBaiaac+facaWGSbaabeaaaOGaayjk aiaawMcaaiaacYcacaWGhbGaamyyaiaadchadaWgaaWcbaGaciyBai aacggacaGG4baabeaaaOGaay5waiaaw2faaaGaay5Eaiaaw2haaaaa @748A@

and the warning message will no longer be output.

Infinite domain detected.

This message can be written during the simulation in the Engine output file.

Such an error message may appear during sorting of contact interfaces when the distance between some nodes of the model becomes infinite.

One has to look at the model behavior to understand why some node coordinates are infinite. One common explanation is "flying" nodes after failure of elements belonging to interface surfaces. In this case using Idel =2 in contact interfaces should fix this problem.

Parameter Bumult.

To detect the secondary nodes and main segments which are in contact, a spatial sort is necessary. Even if using high performance sort algorithms, a full sort at each cycle would not be a fast enough solution. So during sorting, impacts within some security distance up to the gap are retained.

The frequency of interface sorting, that is to say the security distance up to the GAP for storing the candidates to the impact, relies on the parameter Bumult.

It is recommended to not change its value (which is by default 0.20), since it has been determined in order to optimize the performances, based upon our experience.

Recommended tied contact option Spotflag.

For tied contact /INTER/TYPE2, Spotflag=27 should be used for tied contact involving solid elements. Spotflag=28 should be used with tied contact involving shells. Both of these tied contacts use the kinematic method, but will automatically switch to penalty formulation if kinematic conditions are already defined on a secondary node of the tied contact.