உள்ளடக்கத்திற்கு செல்க

நேரியல் நிலையியல் பகுப்பாய்வு (மீள்தன்மை)

இந்தப் பகுப்பாய்வைச் செய்ய tutorial/01_elastic_hinge தரவைப் பயன்படுத்தவும்.

பகுப்பாய்வுப் பொருள்

பகுப்பாய்வுப் பொருள் ஒரு கீல் பகுதி. அதன் வடிவம் படம் 4.1.1-லும் மெஷ் தரவு படம் 4.1.2-லும் காட்டப்பட்டுள்ளது.

உருப்படி விளக்கம் குறிப்பு மேற்கோள்
பகுப்பாய்வு வகை நேரியல் நிலையியல் பகுப்பாய்வு !SOLUTION,TYPE=STATIC
கணுக்களின் எண்ணிக்கை 84,056
உறுப்புகளின் எண்ணிக்கை 49,871
உறுப்பு வகை 10-கணு இரண்டாம் வரிசை நான்முகக் உறுப்பு !ELEMENT,TYPE=342 உறுப்பு நூலகம்
பொருள் பெயர் STEEL !MATERIAL,NAME=STEEL பொருள் தரவு
பொருள் பண்பு ELASTIC !ELASTIC
எல்லை நிபந்தனை கட்டுப்பாடு, செறிவூட்டப்பட்ட சுமை
அணி தீர்வு முறை CG/SSOR !SOLVER,METHOD=CG,PRECOND=1

கீல் பகுதியின் வடிவம்

படம் 4.1.1 கீல் பகுதியின் வடிவம்

கீல் பகுதியின் மெஷ் தரவு

படம் 4.1.2 கீல் பகுதியின் மெஷ் தரவு

பகுப்பாய்வு அமைப்பு

FrontISTR மூலக் குறியீட்டை விரித்து, இந்த எடுத்துக்காட்டின் அடைவிற்கு நகர்ந்து, பகுப்பாய்விற்கு தேவையான கோப்புகள்

கோப்பு பெயர் வகை பங்கு
hecmw_ctrl.dat மொத்த கட்டுப்பாட்டு தரவு மெஷ் தரவும் பகுப்பாய்வு கட்டுப்பாட்டு தரவும் பயன்படுத்தும் உள்ளீட்டு/வெளியீட்டு கோப்புகளை குறிப்பிடுகிறது
hinge.cnt பகுப்பாய்வு கட்டுப்பாட்டு தரவு பகுப்பாய்வு வகை, இடப்பெயர்ச்சி எல்லை நிபந்தனை, செறிவூட்டப்பட்ட சுமை போன்றவற்றை வரையறுத்து, தீர்வான் மற்றும் காட்சிப்படுத்தி கட்டுப்பாடுகளையும் குறிப்பிடுகிறது
hinge.msh மெஷ் தரவு வரையறுக்கப்பட்ட உறுப்பு மெஷையும் அதன் பொருள் தரவு மற்றும் பிரிவு தரவையும் வரையறுக்கிறது

இருப்பதை உறுதிப்படுத்தவும்.

tar xvf FrontISTR.tar.gz
cd FrontISTR/tutorial/01_elastic_hinge
ls

இயக்க முடிவு:

hecmw_ctrl.dat  hinge.cnt  hinge.msh

படம் 4.1.1-ல் காட்டப்பட்ட கட்டுப்பாட்டு பரப்பின் இடப்பெயர்ச்சியை கட்டுப்படுத்தி, கட்டாயப் பரப்பில் செறிவூட்டப்பட்ட சுமை வழங்கும் தகைவு பகுப்பாய்வு செய்யப்படுகிறது.

பகுப்பாய்வு நடைமுறை

FrontISTR இயக்கக் கட்டளை fistr1-ஐ இயக்கவும்.

fistr1 -t 4

(4 த்ரெட்களுடன் இயக்குதல்)

##################################################################
#                         FrontISTR
#
##################################################################
---
version:    5.1.0
git_hash:   acab000c8c633b7b9d596424769e14363f720841
build:
  date:     2020-10-05T07:39:55Z
  MPI:      enabled
  OpenMP:   enabled
  option:   "-p --with-tools --with-refiner --with-metis --with-mumps --with-lapack --with-ml --with-mkl "
  HECMW_METIS_VER: 5
execute:
  date:       2020-10-07T10:01:16+0900
  processes:  1
  threads:    4
  cores:      4
  host:
    0: flow-p06
---
...
 Step control not defined! Using default step=1
 fstr_setup: OK
 Start visualize PSF 1 at timestep 0

 loading step=    1
 sub_step= 1,   current_time=  0.0000E+00, time_inc=  0.1000E+01
 loading_factor=    0.0000000   1.0000000
### 3x3 BLOCK CG, SSOR, 1
      1    1.903375E+00
      2    1.974378E+00
      3    2.534627E+00
...
...
   2967    1.080216E-08
   2968    1.004317E-08
   2969    9.375729E-09
### Relative residual = 9.39429E-09

### summary of linear solver
      2969 iterations      9.394286E-09
    set-up time      :     1.953022E-01
    solver time      :     5.704201E+01
    solver/comm time :     5.145826E-01
    solver/matvec    :     2.306329E+01
    solver/precond   :     2.632665E+01
    solver/1 iter    :     1.921253E-02
    work ratio (%)   :     9.909789E+01

 Start visualize PSF 1 at timestep 1
### FSTR_SOLVE_NLGEOM FINISHED!

 ====================================
    TOTAL TIME (sec) :     59.99
           pre (sec) :      0.71
         solve (sec) :     59.29
 ====================================
 FrontISTR Completed !!

FrontISTR Completed !! காட்டப்பட்டால் பகுப்பாய்வு முடிந்தது.

பகுப்பாய்வு முடிவுகள்

பகுப்பாய்வு முடிந்ததும் சில புதிய கோப்புகள் உருவாக்கப்படுகின்றன.

ls

இயக்க முடிவு:

0.log       hecmw_ctrl.dat  hinge.res.0.0            hinge_vis_psf.0001
FSTR.dbg.0  hecmw_vis.ini   hinge.res.0.1            hinge_vis_psf.0001.pvtu
FSTR.msg    hinge.cnt       hinge_vis_psf.0000
FSTR.sta    hinge.msh       hinge_vis_psf.0000.pvtu

*.res.* என்பது முடிவு தரவு; அதில் FrontISTR பகுப்பாய்வு முடிவுகள் சேமிக்கப்படுகின்றன. REVOCAP_PrePost போன்றவற்றால் அவற்றைக் காணலாம்.

*_vis_* என்பது காட்சிப்படுத்தல் தரவு; பொதுப் பயன்பாட்டு காட்சிப்படுத்தல் மென்பொருளில் இதைக் காணலாம். இந்த எடுத்துக்காட்டில் VTK வடிவம் பயன்படுத்தப்படுவதால் ParaView போன்றவற்றைப் பயன்படுத்தலாம்.

von Mises தகைவின் கண்டூர் வரைபடத்தை REVOCAP_PrePost மூலம் உருவாக்கி படம் 4.1.3-ல் காட்டியுள்ளோம். மேலும் எண்ணியல் முடிவாக பகுப்பாய்வு முடிவு பதிவு கோப்பின் ஒரு பகுதி கீழே தரப்பட்டுள்ளது.

von Mises தகைவு பகுப்பாய்வு முடிவு

படம் 4.1.3 von 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=     1
 ##### Local Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    3.9115E-02     82452 -7.1083E-04     65233
 //U2    7.4504E-05       354 -5.8813E-04       696
 //U3    5.9493E-04        84 -5.8751E-03     61080
 //E11   1.3777E-03       130 -1.3653E-03     77625
 //E22   4.9199E-04        61 -5.4370E-04       102
 //E33   6.8634E-04     51036 -6.1176E-04     30070
 //E12   7.1556E-04     27808 -6.8093E-04     27863
 //E23   5.3666E-04        56 -5.4347E-04        82
 //E31   7.2396E-04     36168 -9.6621E-04       130
 //S11   3.8626E+02       130 -3.6387E+02     28580
 //S22   1.6628E+02       130 -1.5743E+02     28580
 //S33   1.6502E+02     30033 -1.5643E+02     28580
 //S12   5.7795E+01     27808 -5.4998E+01     27863
 //S23   4.3345E+01        56 -4.3896E+01        82
 //S31   5.8474E+01     36168 -7.8040E+01       130
 //SMS   2.8195E+02     77625  1.2755E-02     75112
 ##### Local Summary @Element :Max/IdMax/Min/IdMin####
 //E11   1.0731E-03     10485 -1.2123E-03     41779
 //E22   3.9143E-04     33536 -4.1389E-04     22892
 //E33   5.9415E-04     44563 -5.0497E-04     47965
 //E12   5.3264E-04      9163 -5.0405E-04      9161
 //E23   3.9226E-04     33024 -4.1464E-04     23465
 //E31   5.7633E-04     43142 -4.8019E-04      9571
 //S11   2.7231E+02      9180 -2.9763E+02     41779
 //S22   1.0792E+02      9180 -1.0656E+02     41779
 //S33   1.3921E+02     44569 -1.1431E+02     47974
 //S12   4.3021E+01      9163 -4.0712E+01      9161
 //S23   3.1683E+01     33024 -3.3490E+01     23465
 //S31   4.6550E+01     43142 -3.8785E+01      9571
 //SMS   2.4057E+02     41779  3.1383E-02     38687
 ##### Global Summary @Node    :Max/IdMax/Min/IdMin####
 //U1    3.9115E-02     82452 -7.1083E-04     65233
 //U2    7.4504E-05       354 -5.8813E-04       696
 //U3    5.9493E-04        84 -5.8751E-03     61080
 //E11   1.3777E-03       130 -1.3653E-03     77625
 //E22   4.9199E-04        61 -5.4370E-04       102
 //E33   6.8634E-04     51036 -6.1176E-04     30070
 //E12   7.1556E-04     27808 -6.8093E-04     27863
 //E23   5.3666E-04        56 -5.4347E-04        82
 //E31   7.2396E-04     36168 -9.6621E-04       130
 //S11   3.8626E+02       130 -3.6387E+02     28580
 //S22   1.6628E+02       130 -1.5743E+02     28580
 //S33   1.6502E+02     30033 -1.5643E+02     28580
 //S12   5.7795E+01     27808 -5.4998E+01     27863
 //S23   4.3345E+01        56 -4.3896E+01        82
 //S31   5.8474E+01     36168 -7.8040E+01       130
 //SMS   2.8195E+02     77625  1.2755E-02     75112
 ##### Global Summary @Element :Max/IdMax/Min/IdMin####
 //E11   1.0731E-03     10485 -1.2123E-03     41779
 //E22   3.9143E-04     33536 -4.1389E-04     22892
 //E33   5.9415E-04     44563 -5.0497E-04     47965
 //E12   5.3264E-04      9163 -5.0405E-04      9161
 //E23   3.9226E-04     33024 -4.1464E-04     23465
 //E31   5.7633E-04     43142 -4.8019E-04      9571
 //S11   2.7231E+02      9180 -2.9763E+02     41779
 //S22   1.0792E+02      9180 -1.0656E+02     41779
 //S33   1.3921E+02     44569 -1.1431E+02     47974
 //S12   4.3021E+01      9163 -4.0712E+01      9161
 //S23   3.1683E+01     33024 -3.3490E+01     23465
 //S31   4.6550E+01     43142 -3.8785E+01      9571
 //SMS   2.4057E+02     41779  3.1383E-02     38687