Vai al contenuto

Analisi di contatto (Parte 1)

Questa analisi utilizza i dati in tutorial/09_contact_hertz.

Modello di analisi

Si tratta di un problema di contatto di Hertz. La Figura 4.9.1 mostra la geometria da analizzare e la Figura 4.9.2 mostra i dati di mesh.

Voce Descrizione Note Riferimento
Tipo di analisi Analisi statica non lineare (elastoplasticità, contatto) !SOLUTION,TYPE=NLSTATIC !CONTACT
Numero di nodi 408
Numero di elementi 168
Tipo di elemento Elemento esaedrico del primo ordine a 8 nodi !ELEMENT,TYPE=361
Nome del materiale MAT1 !MATERIAL,NAME=MAT1
Proprietà del materiale ELASTIC !ELASTIC
Condizioni al contorno Vincoli, spostamento prescritto
Solutore lineare CG/SSOR !SOLVER,METHOD=CG,PRECOND=1

Geometria del modello di analisi

Fig. 4.9.1 Geometria del modello di analisi

Dati di mesh del modello di analisi

Fig. 4.9.2 Dati di mesh del modello di analisi

Descrizione dell'analisi

Viene eseguita un'analisi di contatto mediante il metodo dei moltiplicatori di Lagrange aumentato, applicando uno spostamento prescritto in direzione di compressione alla superficie superiore di un modello di quarto di disco. I dati di controllo dell'analisi sono riportati di seguito.

Dati di controllo dell'analisi cgs3.cnt

#  Control File for FISTR
## Analysis Control
!VERSION
 3
!SOLUTION, TYPE=NLSTATIC
!WRITE,RESULT
!WRITE,VISUAL
## Solver Control
### Boundary Condition
!BOUNDARY, GRPID=1
 ALL, 3, 3, 0.0
 BOTTOM, 2, 2, 0.0
 CENTER, 1, 1, 0.0
 UPPER,  2, 2, -0.306
!CONTACT_ALGO, TYPE=ALAGRANGE
!CONTACT, GRPID=1
 CP1, 0.0
### STEP
!STEP, SUBSTEPS=5, CONVERG=1.0e-5
 BOUNDARY, 1
 CONTACT,  1
### Material
!MATERIAL, NAME=MAT1
!ELASTIC
 1100.0, 0.0
### Solver Setting
!SOLVER,METHOD=CG,PRECOND=1,ITERLOG=YES,TIMELOG=YES
 1000, 1
 1.0e-10, 1.0, 0.0
## Post Control
!VISUAL,method=PSR
!surface_num=1
!surface 1
!output_type=VTK
!END

Procedura di analisi

Eseguire il comando eseguibile FrontISTR fistr1.

cd FrontISTR/tutorial/09_contact_hertz
fistr1 -t 4

(Esecuzione con 4 thread)

Risultati dell'analisi

Per il quinto sottopasso, viene creato con REVOCAP_PrePost un diagramma della forma deformata con contorni dello spostamento nella direzione y, mostrato nella Figura 4.9.3. Come risultati numerici dell'analisi, di seguito è riportato un estratto del file di log dei risultati.

Risultati dell'analisi: deformazione e spostamento nella direzione y

Fig. 4.9.3 Risultati dell'analisi: deformazione e spostamento nella direzione y

Log dei risultati dell'analisi 0.log

 fstr_setup: OK
#### Result step=     0
 ##### Local Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    0.0000E+00         1  0.0000E+00         1
 //U2    0.0000E+00         1  0.0000E+00         1
 //U3    0.0000E+00         1  0.0000E+00         1
 //E11   0.0000E+00         1  0.0000E+00         1
 //E22   0.0000E+00         1  0.0000E+00         1
 //E33   0.0000E+00         1  0.0000E+00         1
 //E12   0.0000E+00         1  0.0000E+00         1
 //E23   0.0000E+00         1  0.0000E+00         1
 //E31   0.0000E+00         1  0.0000E+00         1
 //S11   0.0000E+00         1  0.0000E+00         1
 //S22   0.0000E+00         1  0.0000E+00         1
 //S33   0.0000E+00         1  0.0000E+00         1
 //S12   0.0000E+00         1  0.0000E+00         1
 //S23   0.0000E+00         1  0.0000E+00         1
 //S31   0.0000E+00         1  0.0000E+00         1
 //SMS   0.0000E+00         1  0.0000E+00         1
 ##### Local Summary @Element :Max/IdMax/Min/IdMin####
 //E11   0.0000E+00         1  0.0000E+00         1
 //E22   0.0000E+00         1  0.0000E+00         1
 //E33   0.0000E+00         1  0.0000E+00         1
 //E12   0.0000E+00         1  0.0000E+00         1
 //E23   0.0000E+00         1  0.0000E+00         1
 //E31   0.0000E+00         1  0.0000E+00         1
 //S11   0.0000E+00         1  0.0000E+00         1
 //S22   0.0000E+00         1  0.0000E+00         1
 //S33   0.0000E+00         1  0.0000E+00         1
 //S12   0.0000E+00         1  0.0000E+00         1
 //S23   0.0000E+00         1  0.0000E+00         1
 //S31   0.0000E+00         1  0.0000E+00         1
 //SMS   0.0000E+00         1  0.0000E+00         1
 ##### Global Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    0.0000E+00         1  0.0000E+00         1
 //U2    0.0000E+00         1  0.0000E+00         1
 //U3    0.0000E+00         1  0.0000E+00         1
 //E11   0.0000E+00         1  0.0000E+00         1
 //E22   0.0000E+00         1  0.0000E+00         1
 //E33   0.0000E+00         1  0.0000E+00         1
 //E12   0.0000E+00         1  0.0000E+00         1
 //E23   0.0000E+00         1  0.0000E+00         1
 //E31   0.0000E+00         1  0.0000E+00         1
 //S11   0.0000E+00         1  0.0000E+00         1
 //S22   0.0000E+00         1  0.0000E+00         1
 //S33   0.0000E+00         1  0.0000E+00         1
 //S12   0.0000E+00         1  0.0000E+00         1
 //S23   0.0000E+00         1  0.0000E+00         1
 //S31   0.0000E+00         1  0.0000E+00         1
 //SMS   0.0000E+00         1  0.0000E+00         1
 ##### Global Summary @Element :Max/IdMax/Min/IdMin####
 //E11   0.0000E+00         1  0.0000E+00         1
 //E22   0.0000E+00         1  0.0000E+00         1
 //E33   0.0000E+00         1  0.0000E+00         1
 //E12   0.0000E+00         1  0.0000E+00         1
 //E23   0.0000E+00         1  0.0000E+00         1
 //E31   0.0000E+00         1  0.0000E+00         1
 //S11   0.0000E+00         1  0.0000E+00         1
 //S22   0.0000E+00         1  0.0000E+00         1
 //S33   0.0000E+00         1  0.0000E+00         1
 //S12   0.0000E+00         1  0.0000E+00         1
 //S23   0.0000E+00         1  0.0000E+00         1
 //S31   0.0000E+00         1  0.0000E+00         1
 //SMS   0.0000E+00         1  0.0000E+00         1

 ...

#### Result step=     5
 ##### Local Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    1.3079E-02        31 -4.1442E-02        54
 //U2    4.9597E-03      1008 -3.0604E-01        69
 //U3    0.0000E+00         1  0.0000E+00         1
 //E11   7.1171E-03      1049 -4.7627E-02         1
 //E22   2.0644E-03      1012 -5.3428E-02        50
 //E33   9.3301E-03        50 -4.0103E-03        54
 //E12   4.3231E-02      1000 -2.7666E-02        54
 //E23   1.3028E-12      1046 -1.2515E-12         2
 //E31   3.5660E-13        48 -2.1028E-12        51
 //S11   7.9455E+00      1049 -4.9987E+01         1
 //S22   2.2819E+00      1012 -5.3282E+01        50
 //S33   1.1359E+01        50 -4.4937E+00        54
 //S12   2.0802E+01      1000 -1.4668E+01        54
 //S23   6.7682E-10      1046 -6.3666E-10         2
 //S31   2.0235E-10        48 -1.1497E-09        51
 //SMS   6.0944E+01      1000  1.0198E-01        69
 ##### Local Summary @Element :Max/IdMax/Min/IdMin####
 //E11   4.6311E-03      1024 -3.5443E-02         8
 //E22   1.7091E-03      1031 -5.4262E-02         8
 //E33   1.5903E-04         7 -1.6631E-05       101
 //E12   1.8590E-02      1009 -2.6073E-02        14
 //E23   1.5204E-13      1000 -4.1521E-13        16
 //E31   1.3902E-13         5 -8.2391E-13         7
 //S11   5.2516E+00      1024 -3.8099E+01         8
 //S22   1.8872E+00      1031 -5.5262E+01         8
 //S33   2.2824E-01         7 -1.8263E-02       101
 //S12   9.8516E+00      1009 -1.3352E+01        14
 //S23   7.8767E-11      1000 -2.1075E-10        16
 //S31   7.7553E-11         5 -4.4780E-10         7
 //SMS   4.9470E+01         8  2.4786E-01      1007
 ##### Global Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    1.3079E-02        31 -4.1442E-02        54
 //U2    4.9597E-03      1008 -3.0604E-01        69
 //U3    0.0000E+00         1  0.0000E+00         1
 //E11   7.1171E-03      1049 -4.7627E-02         1
 //E22   2.0644E-03      1012 -5.3428E-02        50
 //E33   9.3301E-03        50 -4.0103E-03        54
 //E12   4.3231E-02      1000 -2.7666E-02        54
 //E23   1.3028E-12      1046 -1.2515E-12         2
 //E31   3.5660E-13        48 -2.1028E-12        51
 //S11   7.9455E+00      1049 -4.9987E+01         1
 //S22   2.2819E+00      1012 -5.3282E+01        50
 //S33   1.1359E+01        50 -4.4937E+00        54
 //S12   2.0802E+01      1000 -1.4668E+01        54
 //S23   6.7682E-10      1046 -6.3666E-10         2
 //S31   2.0235E-10        48 -1.1497E-09        51
 //SMS   6.0944E+01      1000  1.0198E-01        69
 ##### Global Summary @Element :Max/IdMax/Min/IdMin####
 //E11   4.6311E-03      1024 -3.5443E-02         8
 //E22   1.7091E-03      1031 -5.4262E-02         8
 //E33   1.5903E-04         7 -1.6631E-05       101
 //E12   1.8590E-02      1009 -2.6073E-02        14
 //E23   1.5204E-13      1000 -4.1521E-13        16
 //E31   1.3902E-13         5 -8.2391E-13         7
 //S11   5.2516E+00      1024 -3.8099E+01         8
 //S22   1.8872E+00      1031 -5.5262E+01         8
 //S33   2.2824E-01         7 -1.8263E-02       101
 //S12   9.8516E+00      1009 -1.3352E+01        14
 //S23   7.8767E-11      1000 -2.1075E-10        16
 //S31   7.7553E-11         5 -4.4780E-10         7
 //SMS   4.9470E+01         8  2.4786E-01      1007