Прескокнете до содржината

Контактна анализа (дел 2)

За оваа анализа се користат податоците од tutorial/10_contact_2tubes.

Предмет на анализа

Анализата е проблем на втиснување на цилиндри; геометријата на предметот на анализа е прикажана на слика 4.10.1, а податоците за мрежата на слика 4.10.2.

Ставка Содржина Забелешки Референца
Тип на анализа Нелинеарна статичка анализа (еластичност, контакт) !SOLUTION,TYPE=NLSTATIC !CONTACT
Број на контактни точки 4,000
Број на елементи 2,888
Тип на елемент 8-јазолен линеарен хексаедарски елемент !ELEMENT,TYPE=361
Име на материјал M1 !MATERIAL,NAME=M1
Материјални својства ELASTIC !ELASTIC
Гранични услови Ограничување, зададено поместување
Матричен решавач Директен метод !SOLVER,METHOD=MUMPS

Облик на предметот на анализа

Слика 4.10.1: Облик на предметот на анализа

Податоци за мрежата на предметот на анализа

Слика 4.10.2: Податоци за мрежата на предметот на анализа

Содржина на анализата

Со методот на Lagrange-ови множители се изведува контактна анализа при која на принудната површина прикажана на слика 4.10.1 ѝ се задава зададено поместување во насока на втиснување. Контролните податоци за анализата се прикажани подолу.

Контролни податоци за анализата 2tubes.cnt

#  Control File for FISTR
## Analysis Control
!VERSION
 3
!SOLUTION, TYPE=NLSTATIC
!WRITE,RESULT
!WRITE,VISUAL
## Solver Control
### Boundary Condition
!BOUNDARY, GRPID=1
  X0, 1, 3, 0.0
  Y0, 2, 2, 0.0
  Z0, 3, 3, 0.0
!BOUNDARY, GRPID=2
  X1, 1, 1, 0.0
!BOUNDARY, GRPID=3
  X1, 1, 1, -1.0
!CONTACT_ALGO, TYPE=SLAGRANGE
!CONTACT, GRPID=1, INTERACTION=FSLID, NPENALTY=1.0e+2
  CP1, 0.0, 1.0e+5
### STEP
!STEP, SUBSTEPS=4, CONVERG=1.0e-5
 BOUNDARY, 1
 BOUNDARY, 3
 CONTACT, 1
### Material
!MATERIAL, NAME=M1
!ELASTIC
 2.1e+5, 0.3
### Solver Setting
!SOLVER,METHOD=MUMPS
## Post Control
!VISUAL,method=PSR
!surface_num=1
!surface 1
!output_type=VTK
!END

Постапка за анализа

Извршете ја командата fistr1 на FrontISTR.

cd FrontISTR/tutorial/10_contact_2tubes
fistr1 -t 4

(Извршување со 4 нишки.)

Резултати од анализата

За резултатите од четвртиот подчекор, со REVOCAP_PrePost е создаден приказ на деформацијата со контури на Mises-овото напрегање и е прикажан на слика 4.10.3. Дел од дневникот на резултатите од анализата е прикажан подолу како нумерички податоци од резултатите.

Резултати од анализата: деформација и Mises-ово напрегање

Слика 4.10.3: Резултати од анализата: деформација и Mises-ово напрегање

Дневник на резултатите од анализата 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=     4
 ##### Local Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    8.7136E-04        32 -1.0021E+00      2006
 //U2    8.7633E-03       104 -7.0546E-03      2006
 //U3    8.7634E-03         4 -7.0546E-03      1901
 //E11   7.5985E-04      1901 -4.1250E-04       105
 //E22   9.8385E-04         2 -9.3128E-04      2058
 //E33   9.8386E-04       102 -9.3120E-04      3843
 //E12   5.3544E-04       133 -2.8361E-04       278
 //E23   1.2509E-03      1901 -1.4195E-03         4
 //E31   5.3554E-04        33 -2.8366E-04      1678
 //S11   7.7172E+01       103 -8.9880E+01       101
 //S22   2.0158E+02         2 -2.2918E+02      1905
 //S33   2.0159E+02       102 -2.2920E+02      2010
 //S12   4.3210E+01       133 -2.2913E+01       278
 //S23   1.0057E+02      1901 -1.1446E+02         4
 //S31   4.3219E+01        33 -2.2918E+01      1678
 //SMS   3.0051E+02      1901  3.1670E+00      2454
 ##### Local Summary @Element :Max/IdMax/Min/IdMin####
 //E11   4.3715E-04      2809 -3.1965E-04      1293
 //E22   8.8407E-04        73 -7.9864E-04      1369
 //E33   8.8408E-04      1365 -7.9867E-04      2809
 //E12   5.3226E-04      1324 -1.9554E-04       188
 //E23   1.0836E-03      2149 -1.0974E-03       684
 //E31   5.3223E-04        32 -1.9559E-04      1176
 //S11   6.6688E+01      1296 -7.9313E+01      1293
 //S22   1.8533E+02        73 -2.1547E+02      1369
 //S33   1.8533E+02      1365 -2.1548E+02      2809
 //S12   4.2951E+01      1324 -1.5800E+01       188
 //S23   8.7338E+01      2149 -8.8751E+01       684
 //S31   4.2949E+01        32 -1.5804E+01      1176
 //SMS   1.9064E+02      1369  3.2782E+00      2126
 ##### Global Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    8.7136E-04        32 -1.0021E+00      2006
 //U2    8.7633E-03       104 -7.0546E-03      2006
 //U3    8.7634E-03         4 -7.0546E-03      1901
 //E11   7.5985E-04      1901 -4.1250E-04       105
 //E22   9.8385E-04         2 -9.3128E-04      2058
 //E33   9.8386E-04       102 -9.3120E-04      3843
 //E12   5.3544E-04       133 -2.8361E-04       278
 //E23   1.2509E-03      1901 -1.4195E-03         4
 //E31   5.3554E-04        33 -2.8366E-04      1678
 //S11   7.7172E+01       103 -8.9880E+01       101
 //S22   2.0158E+02         2 -2.2918E+02      1905
 //S33   2.0159E+02       102 -2.2920E+02      2010
 //S12   4.3210E+01       133 -2.2913E+01       278
 //S23   1.0057E+02      1901 -1.1446E+02         4
 //S31   4.3219E+01        33 -2.2918E+01      1678
 //SMS   3.0051E+02      1901  3.1670E+00      2454
 ##### Global Summary @Element :Max/IdMax/Min/IdMin####
 //E11   4.3715E-04      2809 -3.1965E-04      1293
 //E22   8.8407E-04        73 -7.9864E-04      1369
 //E33   8.8408E-04      1365 -7.9867E-04      2809
 //E12   5.3226E-04      1324 -1.9554E-04       188
 //E23   1.0836E-03      2149 -1.0974E-03       684
 //E31   5.3223E-04        32 -1.9559E-04      1176
 //S11   6.6688E+01      1296 -7.9313E+01      1293
 //S22   1.8533E+02        73 -2.1547E+02      1369
 //S33   1.8533E+02      1365 -2.1548E+02      2809
 //S12   4.2951E+01      1324 -1.5800E+01       188
 //S23   8.7338E+01      2149 -8.8751E+01       684
 //S31   4.2949E+01        32 -1.5804E+01      1176
 //SMS   1.9064E+02      1369  3.2782E+00      2126