අන්තර්ගතය වෙත යන්න

රේඛීය ස්ථිතික විශ්ලේෂණය (ප්‍රත්‍යාස්ථ)

මෙම විශ්ලේෂණය සඳහා 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 හි පෙන්වන සීමා කළ මතුපිටේ විස්ථාපනය සීමා කර, භාර මතුපිටට සංකේන්ද්‍රිත භාරයක් යොදා ආතති විශ්ලේෂණයක් සිදු කරයි. සමස්ත පාලන දත්ත සහ විශ්ලේෂණ පාලන දත්ත පහත පෙන්වයි.

සමස්ත පාලන දත්ත 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

විශ්ලේෂණ වර්ගය, විස්ථාපන මායිම් කොන්දේසි, සංකේන්ද්‍රිත භාර ආදිය නිර්වචනය කර විසඳුම්කාරකය සහ දෘශ්‍යකරණයද පාලනය කරයි.

#  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, 1
 1.0e-08, 1.0, 0.0
## Post Control
!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 ක්‍රියාත්මක කරන්න.

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