Gå til innhold

Translation update pending

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

Lineær statisk analyse (elastisitet)

Denne analysen bruker dataene i tutorial/01_elastic_hinge.

Analyseobjekt

Analyseobjektet er en hengselkomponent. Geometrien vises i figur 4.1.1 og nettdataene i figur 4.1.2.

Element Beskrivelse Merknader Referanse
Analysetype Lineær statisk analyse !SOLUTION,TYPE=STATIC
Antall noder 84,056
Antall elementer 49,871
Elementtype 10-noders kvadratisk tetraederelement !ELEMENT,TYPE=342 Elementbibliotek
Materialnavn STEEL !MATERIAL,NAME=STEEL Materialdata
Materialegenskaper ELASTIC !ELASTIC
Randbetingelser Fastholding, konsentrert last
Matriseløser CG/SSOR !SOLVER,METHOD=CG,PRECOND=1

Geometri for hengselkomponenten

Figur 4.1.1 Geometri for hengselkomponenten

Nettdata for hengselkomponenten

Figur 4.1.2 Nettdata for hengselkomponenten

Analyseinnhold

Pakk ut FrontISTR-kildekoden, gå til katalogen for dette eksemplet, og kontroller at filene som kreves for analysen

Filnavn Type Rolle
hecmw_ctrl.dat Overordnede styringsdata Angir inn- og utdatafilene for nettdataene og analysestyringsdataene
hinge.cnt Analysestyringsdata Definerer analysetype, forskyvningsrandbetingelser, konsentrerte laster osv., og angir også styring av løseren og visualiseringsverktøyet
hinge.msh Nettdata Definerer endelig-element-nettet samt material- og seksjonsdataene

finnes.

tar xvf FrontISTR.tar.gz
cd FrontISTR/tutorial/01_elastic_hinge
ls

Utdata:

hecmw_ctrl.dat  hinge.cnt  hinge.msh

Det utføres en spenningsanalyse der forskyvningen på fastholdingsflatene vist i figur 4.1.1 begrenses og en konsentrert last påføres lastflaten. De overordnede styringsdataene og analysestyringsdataene vises nedenfor.

Overordnede styringsdata hecmw_ctrl.dat

Angir inn- og utdatafilene for nettdataene og analysestyringsdataene.

#
# for solver
#
!MESH, NAME=fstrMSH, TYPE=HECMW-ENTIRE # Angi ett enkelt nettdata
 hinge.msh
!CONTROL, NAME=fstrCNT                 # Angi analysestyringsdata
 hinge.cnt
!RESULT, NAME=fstrRES, IO=OUT          # Angi resultatdata
 hinge.res
!RESULT, NAME=vis_out, IO=OUT          # Angi visualiseringsdata
 hinge_vis

Analysestyringsdata hinge.cnt

Definerer analysetype, forskyvningsrandbetingelser, konsentrerte laster osv., og angir også styring av løseren og visualiseringsverktøyet.

#  Kontrollfil for FISTR
## Analysestyring
!VERSION                   # Angi filformatversjonen
 3
!SOLUTION, TYPE=STATIC     # Angi analysetypen
!WRITE,RESULT              # Angi utdata for resultatdata
!WRITE,VISUAL              # Angi utdata for visualiseringsdata
## Løserstyring
### Randbetingelse
!BOUNDARY
 BND0, 1, 3, 0.000000      # Angi fastholdingsflate 1
!BOUNDARY
 BND1, 1, 3, 0.000000      # Angi fastholdingsflate 2
!CLOAD
 CL0, 1, 0.01000           # Angi lastflaten
### Materiale
!MATERIAL, NAME=STEEL      # Angi materialegenskaper
!ELASTIC                   # Definer elastisk materiale
 210000.0, 0.3
!DENSITY                   # Definer massetetthet
 7.85e-6
### Løserinnstillinger
!SOLVER,METHOD=CG,PRECOND=1,ITERLOG=YES,TIMELOG=YES  # Styr løseren
 10000, 1
 1.0e-08, 1.0, 0.0
## Etterbehandlingsstyring
!VISUAL,method=PSR          # Angi visualiseringsmetoden
!surface_num=1              # Antall flater i 1 flaterendering
!surface 1                  # Angi innholdet i flaten
!output_type=VTK            # Angi typen visualiseringsfil
!END                        # Angi slutten på analysestyringsdataene

Nettdata

Definerer endelig-element-nettet samt material- og seksjonsdataene.

(Kun utdrag)

!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

Analyseprosedyre

Kjør FrontISTR-kommandoen fistr1.

fistr1 -t 4

(Kjøres med 4 tråder)

##################################################################
#                         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 !!

Analysen er fullført når FrontISTR Completed !! vises.

Analyseresultater

Når analysen er fullført, opprettes flere nye filer.

ls

Utdata:

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.* er resultatdatafiler som inneholder analyseresultatene fra FrontISTR. De kan vises med REVOCAP_PrePost og lignende verktøy.

*_vis_* er visualiseringsdata som kan vises med generell visualiseringsprogramvare. I dette eksemplet skrives dataene ut i VTK-format, så de kan vises med ParaView eller lignende.

Et konturplott av von Mises-spenningen opprettet med REVOCAP_PrePost vises i figur 4.1.3. Som numeriske analyseresultater vises også en del av loggfilen nedenfor.

Analyseresultat for von Mises-spenning

Figur 4.1.3 Analyseresultat for von Mises-spenning

Logg for analyseresultater 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