መስመራዊ ስታቲክ ትንተና (ኤላስቲክ፣ ፓራሌል)¶
የስታቲክ ትንተና (ኤላስቲክ)ን በ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¶
#
# ለክፍፍል ፕሮግራም
#
!MESH, NAME=part_in,TYPE=HECMW-ENTIRE # በhecmw_part1 የሚከፈለው የመነሻ ሜሽ ውሂብ
hinge.msh
!MESH, NAME=part_out,TYPE=HECMW-DIST # በhecmw_part1 ከተከፈለ በኋላ የፋይሉ ስም
hinge_4
#
# ለሶልቨር
#
!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 # በ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-መንገድ ፓራሌል፣ በ1 OpenMP ትሬድ ማስኬድ)
##################################################################
# 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.* የውጤት ውሂብ ሲሆን፣ ለእያንዳንዱ የFrontISTR MPI ኖድ የትንተና ውጤቶችን ይይዛል። በREVOCAP_PrePost ወዘተ ሊታይ ይችላል።
*_vis_* የምስላዊ ማሳያ ውሂብ ሲሆን፣ በአጠቃላይ የምስላዊ ማሳያ ሶፍትዌር ሊታይ ይችላል። በዚህ ምሳሌ በVTK ቅርጸት ስለሚወጣ፣ ParaView ወዘተ በመጠቀም ይመልከቱ።
ስዕል 4.2.2 የሚዘስ ስትሬስ ትንተና ውጤት (በ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