Análisis estático lineal (Elasticidad)¶
Este análisis utiliza los datos de tutorial/01_elastic_hinge.
Objeto del análisis¶
El objeto del análisis es una pieza de bisagra. Su geometría se muestra en la Figura 4.1.1 y los datos de malla se muestran en la Figura 4.1.2.
| Elemento | Descripción | Notas | Referencia |
|---|---|---|---|
| Tipo de análisis | Análisis estático lineal | !SOLUTION,TYPE=STATIC | |
| Número de nodos | 84,056 | ||
| Número de elementos | 49,871 | ||
| Tipo de elemento | Elemento tetraédrico cuadrático de 10 nodos | !ELEMENT,TYPE=342 | Biblioteca de elementos |
| Nombre de la propiedad del material | STEEL | !MATERIAL,NAME=STEEL | Datos de material |
| Comportamiento del material | ELASTIC | !ELASTIC | |
| Condiciones de contorno | Restricciones, carga concentrada | ||
| Solver matricial | CG/SSOR | !SOLVER,METHOD=CG,PRECOND=1 |
Figura 4.1.1 Geometría de la pieza de bisagra
Figura 4.1.2 Datos de malla de la pieza de bisagra
Descripción del análisis¶
Extraiga el código fuente de FrontISTR, desplácese al directorio de este ejemplo y confirme que están presentes los archivos necesarios para el análisis
| Nombre del archivo | Tipo | Función |
|---|---|---|
hecmw_ctrl.dat |
Datos de control general | Especifica los archivos de entrada/salida para los datos de malla y los datos de control del análisis |
hinge.cnt |
Datos de control del análisis | Define el tipo de análisis, las condiciones de contorno de desplazamiento, las cargas concentradas, etc., y especifica el control del solver y del visualizador |
hinge.msh |
Datos de malla | Define la malla de elementos finitos, los datos de material y los datos de sección |
.
Salida:
Se realiza un análisis de tensiones restringiendo el desplazamiento de las superficies restringidas mostradas en la Figura 4.1.1 y aplicando una carga concentrada a la superficie de carga. A continuación se muestran los datos de control general y los datos de control del análisis.
Datos de control general hecmw_ctrl.dat¶
Especifica los archivos de entrada/salida para los datos de malla y los datos de control del análisis.
#
# for solver
#
!MESH, NAME=fstrMSH, TYPE=HECMW-ENTIRE # Especificar los datos de malla única
hinge.msh
!CONTROL, NAME=fstrCNT # Especificar los datos de control del análisis
hinge.cnt
!RESULT, NAME=fstrRES, IO=OUT # Especificar los datos de resultados
hinge.res
!RESULT, NAME=vis_out, IO=OUT # Especificar los datos de visualización
hinge_vis
Datos de control del análisis hinge.cnt¶
Define el tipo de análisis, las condiciones de contorno de desplazamiento, las cargas concentradas, etc., y también especifica el control del solver y del visualizador.
# Control File for FISTR
## Analysis Control
!VERSION # Especificar la versión del formato de archivo
3
!SOLUTION, TYPE=STATIC # Especificar el tipo de análisis
!WRITE,RESULT # Especificar la salida de datos de resultados
!WRITE,VISUAL # Especificar la salida de datos de visualización
## Solver Control
### Boundary Condition
!BOUNDARY
BND0, 1, 3, 0.000000 # Especificar la superficie restringida 1
!BOUNDARY
BND1, 1, 3, 0.000000 # Especificar la superficie restringida 2
!CLOAD
CL0, 1, 0.01000 # Especificar la superficie de carga
### Material
!MATERIAL, NAME=STEEL # Especificar las propiedades del material
!ELASTIC # Definir el material elástico
210000.0, 0.3
!DENSITY # Definir la densidad de masa
7.85e-6
### Solver Setting
!SOLVER,METHOD=CG,PRECOND=1,ITERLOG=YES,TIMELOG=YES # Control del solucionador
10000, 1
1.0e-08, 1.0, 0.0
## Post Control
!VISUAL,method=PSR # Especificar el método de visualización
!surface_num=1 # Número de superficies en una renderización de superficie
!surface 1 # Especificar el contenido de la superficie
!output_type=VTK # Especificar el tipo de archivo de visualización
!END # Indicar el final de los datos de control del análisis
Datos de malla¶
Define la malla de elementos finitos y sus datos de material y sección.
(Extracto)
!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
Procedimiento de análisis¶
Ejecute el comando de FrontISTR fistr1.
(Ejecutar con 4 hilos)
##################################################################
# 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 !!
El análisis ha finalizado cuando se muestra FrontISTR Completed !!.
Resultados del análisis¶
Al finalizar el análisis se crean varios archivos nuevos.
Salida:
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
Los archivos *.res.* son datos de resultados que contienen los resultados del análisis de FrontISTR. Pueden visualizarse con REVOCAP_PrePost y herramientas similares.
Los archivos *_vis_* son datos de visualización y pueden mostrarse con software de visualización de propósito general. En este ejemplo, la salida está en formato VTK, por lo que puede visualizarse con ParaView o software similar.
En la Figura 4.1.3 se muestra un gráfico de contornos de tensión de von Mises creado con REVOCAP_PrePost. Además, a continuación se muestra parte del archivo de registro de resultados del análisis como datos numéricos de los resultados.
Figura 4.1.3 Resultado del análisis de tensión de von Mises
Registro de resultados del análisis 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