လိုင်းနီးယား တည်ငြိမ်ခွဲခြမ်းစိတ်ဖြာမှု (အီလတ်စတစ်စတီ, အပြိုင်)¶
တည်ငြိမ် ခွဲခြမ်းစိတ်ဖြာမှု (အီလတ်စတစ်စတီ) ကို 4 ခွဲ အပြိုင် ဖြင့် လုပ်ဆောင်ရန် tutorial/02_elastic_hinge_parallel ရှိ ဒေတာ ကို အသုံးပြုသည်။
ခွဲခြမ်းစိတ်ဖြာမည့် အရာ¶
| အချက် | အကြောင်းအရာ | မှတ်ချက် | ကိုးကားချက် |
|---|---|---|---|
| ခွဲခြမ်းစိတ်ဖြာမှုအမျိုးအစား | လိုင်းနီးယား တည်ငြိမ်ခွဲခြမ်းစိတ်ဖြာမှု | !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.2.1 နိုဒ် တစ်ခုချင်းစီ၏ တွက်ချက်မှု ဒိုမိန်း
ခွဲခြမ်းစိတ်ဖြာမှု အကြောင်းအရာ¶
FrontISTR ကုဒ် ကို ဖြည်ထုတ် လုပ်ပြီး ဤ ဥပမာ ဒိုင်ရက်ထရီ သို့ ရွှေ့ကာ ခွဲခြမ်းစိတ်ဖြာမှု အတွက် လိုအပ်သော အောက်ပါ ဖိုင် များ
| ဖိုင်အမည် | အမျိုးအစား | အခန်းကဏ္ဍ |
|---|---|---|
| hecmw_ctrl.dat | အလုံးစုံထိန်းချုပ်ဒေတာ | |
| hinge.cnt | ခွဲခြမ်းစိတ်ဖြာမှု ထိန်းချုပ်ဒေတာ | |
| hinge.msh | မက်ရှ်ဒေတာ | |
| hecmw_part_ctrl.dat | ဒိုမိန်းခွဲခြမ်းခြင်း ထိန်းချုပ်ဒေတာ | hecmw_part1 ဖြင့် မက်ရှ်ဒေတာ ကို ဒိုမိန်းခွဲခြမ်းခြင်း လုပ်ရာတွင် သုံးသော ထိန်းချုပ်ဒေတာ |
ရှိကြောင်း စစ်ဆေးပါ။
အထွက်:
ကန့်သတ်ထားသော မျက်နှာပြင် ၏ ရွှေ့ပြောင်းမှု ကို ကန့်သတ် လုပ်ပြီး ဝန် မျက်နှာပြင် သို့ စုစည်းဝန် ပေးသော စထရက်စ် ခွဲခြမ်းစိတ်ဖြာမှု ကို လုပ်ဆောင်သည်။
အလုံးစုံထိန်းချုပ်ဒေတာ၊ ခွဲခြမ်းစိတ်ဖြာမှု ထိန်းချုပ်ဒေတာ နှင့် ဒိုမိန်းခွဲခြမ်းခြင်း ထိန်းချုပ်ဒေတာ ကို အောက်တွင် ပြထားသည်။
အလုံးစုံထိန်းချုပ်ဒေတာ hecmw_ctrl.dat¶
#
# for partitioner
#
!MESH, NAME=part_in,TYPE=HECMW-ENTIRE # hecmw_part1 ဖြင့် အပိုင်းခွဲမည့် မူလ မက်ရှ်ဒေတာ
hinge.msh
!MESH, NAME=part_out,TYPE=HECMW-DIST # hecmw_part1 ဖြင့် အပိုင်းခွဲပြီးနောက် ဖိုင်အမည်
hinge_4
#
# for solver
#
!MESH, NAME=fstrMSH, TYPE=HECMW-DIST # အပိုင်းခွဲပြီးသော မက်ရှ်ဒေတာကို သတ်မှတ်ရန်
hinge_4
!CONTROL, NAME=fstrCNT # ခွဲခြမ်းစိတ်ဖြာမှု ထိန်းချုပ်ဒေတာကို သတ်မှတ်ရန်
hinge.cnt
!RESULT, NAME=fstrRES, IO=OUT # ရလဒ်ဒေတာကို သတ်မှတ်ရန်
hinge.res
!RESULT, NAME=vis_out, IO=OUT # ပုံဖော်ဒေတာကို သတ်မှတ်ရန်
hinge_vis
ခွဲခြမ်းစိတ်ဖြာမှု ထိန်းချုပ်ဒေတာ hinge.cnt¶
# FISTR အတွက် ထိန်းချုပ်မှုဖိုင်
## ခွဲခြမ်းစိတ်ဖြာမှု ထိန်းချုပ်ခြင်း
!VERSION # ဖိုင်ဖော်မတ်ဗားရှင်းကို သတ်မှတ်ရန်
3
!SOLUTION, TYPE=STATIC # ခွဲခြမ်းစိတ်ဖြာမှုအမျိုးအစားကို သတ်မှတ်ရန်
!WRITE,RESULT # ရလဒ်ဒေတာထုတ်ပေးမှုကို သတ်မှတ်ရန်
!WRITE,VISUAL # ပုံဖော်ဒေတာထုတ်ပေးမှုကို သတ်မှတ်ရန်
## ဖြေရှင်းကိရိယာ ထိန်းချုပ်ခြင်း
### နယ်နိမိတ်အခြေအနေ
!BOUNDARY
BND0, 1, 3, 0.000000 # ကန့်သတ်မျက်နှာပြင် 1 ကို သတ်မှတ်ရန်
!BOUNDARY
BND1, 1, 3, 0.000000 # ကန့်သတ်မျက်နှာပြင် 2 ကို သတ်မှတ်ရန်
!CLOAD
CL0, 1, 0.01000 # ဝန်သက်ရောက်မျက်နှာပြင်ကို သတ်မှတ်ရန်
### ပစ္စည်း
!MATERIAL, NAME=STEEL # ပစ္စည်းဂုဏ်သတ္တိကို သတ်မှတ်ရန်
!ELASTIC # အီလတ်စတစ်ပစ္စည်းကို သတ်မှတ်ရန်
210000.0, 0.3
!DENSITY # ဒြပ်ထုသိပ်သည်းဆကို သတ်မှတ်ရန်
7.85e-6
### ဖြေရှင်းကိရိယာ ဆက်တင်
!SOLVER,METHOD=CG,PRECOND=1,ITERLOG=YES,TIMELOG=YES # ဖြေရှင်းကိရိယာကို ထိန်းချုပ်ရန်
10000, 2
1.0e-08, 1.0, 0.0
## နောက်ဆက်တွဲလုပ်ငန်း ထိန်းချုပ်ခြင်း
!VISUAL,method=PSR # ပုံဖော်နည်းလမ်းကို သတ်မှတ်ရန်
!surface_num=1 # မျက်နှာပြင်ပုံဖော်မှုတစ်ခုအတွင်း မျက်နှာပြင်အရေအတွက်
!surface 1 # မျက်နှာပြင်အကြောင်းအရာကို သတ်မှတ်ရန်
!output_type=VTK # ပုံဖော်ဖိုင်အမျိုးအစားကို သတ်မှတ်ရန်
!END # ခွဲခြမ်းစိတ်ဖြာမှု ထိန်းချုပ်ဒေတာ၏ အဆုံးကို ဖော်ပြရန်
ဒိုမိန်းခွဲခြမ်းခြင်း ထိန်းချုပ်ဒေတာ hecmw_part_ctrl.dat¶
ခွဲခြမ်းစိတ်ဖြာမှု လုပ်ငန်းစဉ်¶
FrontISTR ကို MPI ဖြင့် လုပ်ဆောင်ရန် ပထမဦးစွာ မက်ရှ်ဒေတာ hinge.msh ကို ဒိုမိန်း 4 ခုသို့ ခွဲသည်။
အထွက်:
Oct 07 11:04:52 Info: Reading mesh file...
Oct 07 11:04:52 Info: Starting domain decomposition...
Oct 07 11:04:52 TH(0/8) Info: Creating local mesh for domain #0 ...
Oct 07 11:04:52 TH(2/8) Info: Creating local mesh for domain #1 ...
Oct 07 11:04:52 TH(6/8) Info: Creating local mesh for domain #2 ...
Oct 07 11:04:52 TH(7/8) Info: Creating local mesh for domain #3 ...
Oct 07 11:04:52 Info: Domain decomposition done
hinge_4.x နှင့် part.inp ဟု အမည်ရှိသော ဖိုင် အသစ်များ ဖန်တီးသည်။
အထွက်:
hecmw_ctrl.dat hecmw_part_ctrl.dat hinge.msh hinge_4.1 hinge_4.3
hecmw_part.log hinge.cnt hinge_4.0 hinge_4.2 part.inp
ထို့နောက် FrontISTR လုပ်ဆောင်မှု အမိန့် fistr1 ကို MPI ဖြင့် လုပ်ဆောင်သည်။
(MPI 4 ခွဲ အပြိုင်, OpenMP သရက်ဒ် 1 ခုဖြင့် လုပ်ဆောင်ခြင်း)
##################################################################
# 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-07T11:07:21+0900
processes: 4
threads: 1
cores: 4
host:
0: flow-p06
1: flow-p06
2: flow-p06
3: 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, 2
1 2.183567E+00
2 2.423900E+00
3 2.939117E+00
...
2084 1.158654E-08
2085 1.032414E-08
2086 9.436273E-09
### Relative residual = 9.43589E-09
### summary of linear solver
2086 iterations 9.435886E-09
set-up time : 4.695220E-02
solver time : 7.103976E+01
solver/comm time : 1.929294E+01
solver/matvec : 1.544405E+01
solver/precond : 3.243278E+01
solver/1 iter : 3.405549E-02
work ratio (%) : 7.284205E+01
Start visualize PSF 1 at timestep 1
### FSTR_SOLVE_NLGEOM FINISHED!
====================================
TOTAL TIME (sec) : 72.42
pre (sec) : 0.29
solve (sec) : 72.13
====================================
FrontISTR Completed !!
FrontISTR Completed !! ပြသလျှင် ခွဲခြမ်းစိတ်ဖြာမှု ပြီးဆုံးသည်။
ခွဲခြမ်းစိတ်ဖြာမှု ရလဒ်¶
ခွဲခြမ်းစိတ်ဖြာမှု ပြီးဆုံးသည့်အခါ ဖိုင် အသစ်အချို့ ဖန်တီးသည်။
အထွက်:
0.log FSTR.dbg.3 hinge.cnt hinge.res.2.1 hinge_vis_psf.0000
1.log FSTR.msg hinge.msh hinge.res.3.0 hinge_vis_psf.0000.pvtu
2.log FSTR.sta hinge.res.0.0 hinge.res.3.1 hinge_vis_psf.0001
3.log hecmw_ctrl.dat hinge.res.0.1 hinge_4.0 hinge_vis_psf.0001.pvtu
FSTR.dbg.0 hecmw_part.log hinge.res.1.0 hinge_4.1 part.inp
FSTR.dbg.1 hecmw_part_ctrl.dat hinge.res.1.1 hinge_4.2
FSTR.dbg.2 hecmw_vis.ini hinge.res.2.0 hinge_4.3
*.res.* များသည် ရလဒ်ဒေတာ ဖြစ်ပြီး MPI နိုဒ် တစ်ခုချင်းစီ၏ FrontISTR ခွဲခြမ်းစိတ်ဖြာမှုရလဒ် ပါဝင်သည်။ REVOCAP_PrePost စသည့် ကိရိယာ များဖြင့် ပြသနိုင်သည်။
*_vis_* များသည် မြင်သာပြဒေတာ ဖြစ်ပြီး အထွေထွေသုံး မြင်သာပြခြင်း ဆော့ဖ်ဝဲ ဖြင့် ပြသနိုင်သည်။ ဤ ဥပမာ တွင် VTK ဖော်မတ် ဖြင့် အထွက် လုပ်ထားသဖြင့် ParaView စသည့် ဆော့ဖ်ဝဲ ကို အသုံးပြု၍ ပြသသည်။
ပုံ 4.2.2 von Mises စထရက်စ် ခွဲခြမ်းစိတ်ဖြာမှုရလဒ် (VTK ဖြင့် ပြသ)
ခွဲခြမ်းစိတ်ဖြာမှု ရလဒ် မှတ်တမ်း 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 14 0.0000E+00 14
//U2 0.0000E+00 14 0.0000E+00 14
//U3 0.0000E+00 14 0.0000E+00 14
//E11 0.0000E+00 14 0.0000E+00 14
//E22 0.0000E+00 14 0.0000E+00 14
//E33 0.0000E+00 14 0.0000E+00 14
//E12 0.0000E+00 14 0.0000E+00 14
//E23 0.0000E+00 14 0.0000E+00 14
//E31 0.0000E+00 14 0.0000E+00 14
//S11 0.0000E+00 14 0.0000E+00 14
//S22 0.0000E+00 14 0.0000E+00 14
//S33 0.0000E+00 14 0.0000E+00 14
//S12 0.0000E+00 14 0.0000E+00 14
//S23 0.0000E+00 14 0.0000E+00 14
//S31 0.0000E+00 14 0.0000E+00 14
//SMS 0.0000E+00 14 0.0000E+00 14
##### Global Summary @Element :Max/IdMax/Min/IdMin####
//E11 0.0000E+00 14 0.0000E+00 14
//E22 0.0000E+00 14 0.0000E+00 14
//E33 0.0000E+00 14 0.0000E+00 14
//E12 0.0000E+00 14 0.0000E+00 14
//E23 0.0000E+00 14 0.0000E+00 14
//E31 0.0000E+00 14 0.0000E+00 14
//S11 0.0000E+00 14 0.0000E+00 14
//S22 0.0000E+00 14 0.0000E+00 14
//S33 0.0000E+00 14 0.0000E+00 14
//S12 0.0000E+00 14 0.0000E+00 14
//S23 0.0000E+00 14 0.0000E+00 14
//S31 0.0000E+00 14 0.0000E+00 14
//SMS 0.0000E+00 14 0.0000E+00 14
#### Result step= 1
##### Local Summary @Node :Max/IdMax/Min/IdMin####
//U1 7.4906E-03 50978 -6.9368E-04 51989
//U2 6.1150E-05 2426 -1.2781E-04 40202
//U3 1.8720E-04 2712 -5.5024E-03 67858
//E11 1.3024E-03 28591 -1.2844E-03 28580
//E22 2.1102E-04 24593 -2.0594E-04 205
//E33 6.8634E-04 51036 -5.9603E-04 67794
//E12 5.2186E-04 209 -5.4559E-04 28641
//E23 2.5417E-04 28521 -3.2956E-04 187
//E31 7.2396E-04 36168 -9.5847E-04 28591
//S11 3.6844E+02 28591 -3.6387E+02 28580
//S22 1.5841E+02 28591 -1.5743E+02 28580
//S33 1.5959E+02 51036 -1.5643E+02 28580
//S12 4.2150E+01 209 -4.4067E+01 28641
//S23 2.0529E+01 28521 -2.6618E+01 187
//S31 5.8474E+01 36168 -7.7415E+01 28591
//SMS 2.5134E+02 28591 3.5157E-02 12726
##### Local Summary @Element :Max/IdMax/Min/IdMin####
//E11 1.0115E-03 9959 -9.7870E-04 9942
//E22 1.5457E-04 13694 -1.4639E-04 7189
//E33 5.9415E-04 44563 -5.0497E-04 47965
//E12 3.5513E-04 9579 -3.9083E-04 9600
//E23 2.0608E-04 9563 -2.3480E-04 9589
//E31 5.7633E-04 43142 -4.8019E-04 9571
//S11 2.6038E+02 9573 -2.5605E+02 9571
//S22 1.0144E+02 9573 -1.0100E+02 9571
//S33 1.3921E+02 44569 -1.1431E+02 47974
//S12 2.8683E+01 9579 -3.1567E+01 9600
//S23 1.6645E+01 9563 -1.8964E+01 9589
//S31 4.6550E+01 43142 -3.8785E+01 9571
//SMS 2.0350E+02 9959 3.4602E-02 28
##### 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