လိုင်းနီးယား တည်ငြိမ် ခွဲခြမ်းစိတ်ဖြာမှု (အီလတ်စတစ်စတီ)¶
ဤ ခွဲခြမ်းစိတ်ဖြာမှု ကို လုပ်ဆောင်ရန် 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 | မက်ရှ် ဒေတာ | ဖိုင်နိုက်-အီလီမင့် မက်ရှ် ကို သတ်မှတ်ပြီး ပစ္စည်း ဒေတာ နှင့် အပိုင်း ဒေတာ ကို သတ်မှတ်သည်။ |
ရှိကြောင်း စစ်ဆေးပါ။
လုပ်ဆောင်မှုရလဒ်:
ပုံ 4.1.1 တွင် ပြထားသော ကန့်သတ်ချက် မျက်နှာပြင် ၏ ရွှေ့ပြောင်းမှု ကို ကန့်သတ်ပြီး ဝန် မျက်နှာပြင် ပေါ်သို့ စုစည်းဝန် တင်ကာ စထရက်စ် ခွဲခြမ်းစိတ်ဖြာမှု ကို လုပ်ဆောင်သည်။ အလုံးစုံထိန်းချုပ်ဒေတာ နှင့် ခွဲခြမ်းစိတ်ဖြာမှု ထိန်းချုပ်ဒေတာ ကို အောက်တွင် ပြထားသည်။
အလုံးစုံထိန်းချုပ်ဒေတာ hecmw_ctrl.dat¶
မက်ရှ် ဒေတာ နှင့် ခွဲခြမ်းစိတ်ဖြာမှု ထိန်းချုပ်ဒေတာ ၏ ထည့်သွင်း/အထွက် ဖိုင်များကို သတ်မှတ်သည်။
#
# for solver
#
!MESH, NAME=fstrMSH, TYPE=HECMW-ENTIRE # တစ်ခုတည်း မက်ရှ် ဒေတာ ကို သတ်မှတ်သည်
hinge.msh
!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, 1
1.0e-08, 1.0, 0.0
## နောက်ဆက်တွဲလုပ်ငန်း ထိန်းချုပ်ခြင်း
!VISUAL,method=PSR # မြင်သာဖော်ပြမှု နည်းလမ်း ကို သတ်မှတ်သည်
!surface_num=1 # မျက်နှာပြင် ပုံဖော်ပြသမှု တစ်ခုအတွင်း မျက်နှာပြင် အရေအတွက်
!surface 1 # မျက်နှာပြင် အကြောင်းအရာကို သတ်မှတ်သည်
!output_type=VTK # မြင်သာဖော်ပြမှု ဖိုင်အမျိုးအစား ကို သတ်မှတ်သည်
!END # ခွဲခြမ်းစိတ်ဖြာမှု ထိန်းချုပ်ဒေတာ ၏ အဆုံးကို ပြသည်
မက်ရှ် ဒေတာ¶
ဖိုင်နိုက်-အီလီမင့် မက်ရှ် ကို သတ်မှတ်ပြီး ၎င်း၏ ပစ္စည်း ဒေတာ နှင့် အပိုင်း ဒေတာ ကို သတ်မှတ်သည်။
(တစ်စိတ်တစ်ပိုင်းသာ)
!HEADER
HECMW_Msh File generated by REVOCAP
!NODE
1, -1.22042, 2.23355, 1.65220
2, -1.27050, -3.10529, 1.59209
...
!ELEMENT, TYPE=342
1, 1157, 3549, 3321, 3739, 12629, 12627, 12626, 12628, 12631, 12630
2, 8207, 3321, 3549, 3739, 12629, 12633, 12632, 12634, 12630, 12631
...
!MATERIAL, NAME=STEEL, ITEM=2
!ITEM=1, SUBITEM=2
210000.0, 0.3
!ITEM=2, SUBITEM=1
7.85e-6
!SECTION, TYPE=SOLID, EGRP=Solid0, MATERIAL=STEEL
!EGROUP, EGRP=Solid0
1
2
...
!END
ခွဲခြမ်းစိတ်ဖြာမှု လုပ်ငန်းစဉ်¶
FrontISTR လုပ်ဆောင်မှု ကွန်မန်း fistr1 ကို လုပ်ဆောင်ပါ။
(သရက်ဒ် 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 !! ဟု ပြသပါက ခွဲခြမ်းစိတ်ဖြာမှု ပြီးဆုံးပါပြီ။
ခွဲခြမ်းစိတ်ဖြာမှု ရလဒ်¶
ခွဲခြမ်းစိတ်ဖြာမှု ပြီးဆုံးသောအခါ ဖိုင်အချို့ကို အသစ်ဖန်တီးမည်။
လုပ်ဆောင်မှုရလဒ်:
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 စသည်တို့ကို အသုံးပြု၍ ကြည့်ရှုသည်။
Mises စထရက်စ် ကွန်တို ကို REVOCAP_PrePost ဖြင့် ဖန်တီးပြီး ပုံ 4.1.3 တွင် ပြထားသည်။ ထို့အပြင် ခွဲခြမ်းစိတ်ဖြာမှုရလဒ် ၏ ကိန်းဂဏန်း ဒေတာ အဖြစ် ခွဲခြမ်းစိတ်ဖြာမှု-ရလဒ် လော့ဂ် ဖိုင် ၏ တစ်စိတ်တစ်ပိုင်းကို အောက်တွင် ပြထားသည်။
ပုံ 4.1.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= 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