Skip to content

രേഖീയ സ്റ്റാറ്റിക് വിശകലനം (ഇലാസ്റ്റിക്)

ഈ വിശകലനം നടത്താൻ 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-ൽ കാണിച്ചിരിക്കുന്ന നിയന്ത്രിത ഉപരിതലത്തിലെ സ്ഥാനചലനം നിയന്ത്രിച്ച്, ബലം പ്രയോഗിക്കുന്ന ഉപരിതലത്തിൽ കേന്ദ്രീകൃത ലോഡ് നൽകി സ്ട്രെസ് വിശകലനം നടത്തുന്നു.

വിശകലന നടപടിക്രമം

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 പോലുള്ള സോഫ്റ്റ്‌വെയർ ഉപയോഗിച്ച് പ്രദർശിപ്പിക്കുക.

REVOCAP_PrePost ഉപയോഗിച്ച് Mises സ്ട്രെസിന്റെ കോൺടൂർ ചിത്രം സൃഷ്ടിച്ച് ചിത്രം 4.1.3-ൽ കാണിക്കുന്നു. കൂടാതെ വിശകലന ഫലങ്ങളുടെ സംഖ്യാത്മക ഡാറ്റയായി, വിശകലന ഫല ലോഗ് ഫയലിന്റെ ഒരു ഭാഗം താഴെ കാണിക്കുന്നു.

Mises സ്ട്രെസ് വിശകലന ഫലം

ചിത്രം 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