Անցնել պարունակությանը

Translation update pending

The Armenian translation may be out of date because the Japanese source has changed. The existing translation remains available while it is being updated.

Գծային ստատիկ վերլուծություն (առաձգականություն, զուգահեռ)

Ստատիկ (առաձգական) վերլուծությունը 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-ով ցանցը տիրույթների բաժանելու համար

։

tar xvf FrontISTR.tar.gz
cd FrontISTR/tutorial/02_elastic_hinge_parallel
ls

Կատարման արդյունք՝

hecmw_ctrl.dat  hecmw_part_ctrl.dat  hinge.cnt  hinge.msh

Ամրացված մակերեսի տեղաշարժը սահմանափակվում է, իսկ բեռնավորման մակերեսին տրվում է կենտրոնացված բեռ և կատարվում է լարման վերլուծություն։

Ընդհանուր կառավարման, վերլուծության կառավարման և տիրույթի բաժանման կառավարման տվյալները բերված են ստորև։

Ընդհանուր կառավարման տվյալներ 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

#  Control File for FISTR
## Analysis Control
!VERSION                   # նշել ֆայլի ձևաչափի տարբերակը
 3
!SOLUTION, TYPE=STATIC     # նշել վերլուծության տեսակը
!WRITE,RESULT              # նշել արդյունքների տվյալների արտածումը
!WRITE,VISUAL              # նշել տեսապատկերման տվյալների արտածումը
## Solver Control
### Boundary Condition
!BOUNDARY
 BND0, 1, 3, 0.000000      # նշել ամրացված մակերես 1-ը
!BOUNDARY
 BND1, 1, 3, 0.000000      # նշել ամրացված մակերես 2-ը
!CLOAD
 CL0, 1, 0.01000           # նշել բեռնված մակերեսը
### Material
!MATERIAL, NAME=STEEL      # նշել նյութի հատկությունները
!ELASTIC                   # սահմանել առաձգական նյութը
 210000.0, 0.3
!DENSITY                   # սահմանել զանգվածային խտությունը
 7.85e-6
### Solver Setting
!SOLVER,METHOD=CG,PRECOND=1,ITERLOG=YES,TIMELOG=YES  # լուծիչի կառավարում
 10000, 2
 1.0e-08, 1.0, 0.0
## Post Control
!VISUAL,method=PSR         # նշել տեսապատկերման եղանակը
!surface_num=1             # մակերեսների քանակը մեկ մակերեսային արտապատկերման մեջ
!surface 1                 # նշել մակերեսի բովանդակությունը
!output_type=VTK           # նշել տեսապատկերման ֆայլի տեսակը
!END                       # նշում է վերլուծության կառավարման տվյալների ավարտը

Տիրույթի բաժանման կառավարման տվյալներ hecmw_part_ctrl.dat

!PARTITION,TYPE=NODE-BASED,METHOD=PMETIS,DOMAIN=4,UCD=part.inp

Վերլուծության ընթացակարգ

FrontISTR-ը MPI-ով գործարկելու համար նախ ցանցային hinge.msh տվյալները բաժանեք 4 տիրույթի։

hecmw_part1

Կատարման արդյունք՝

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 ֆայլերը։

ls

Կատարման արդյունք՝

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-ով։

mpirun -np 4 fistr1 -t 1

(կատարում 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 !!, վերլուծությունն ավարտված է։

Վերլուծության արդյունքներ

Վերլուծության ավարտից հետո կստեղծվեն մի քանի նոր ֆայլեր։

ls

Կատարման արդյունք՝

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-ով և նման ծրագրերով։

Mises-ի լարման վերլուծության արդյունքը (VTK դիտում)

Նկար 4.2.2 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