Integración de ANSYS Mechanical APDL y
MATLAB
Enlace de ANSYS y MATLAB para la ejecución de la simulación y obtención de resultados en un archivo formato txt.
Marilyn Coro P., Esteban Guayasamín
Departamento de Ciencias de la Energía y Mecánica, Universidad de las Fuerzas Armadas ESPE
Sangolquí-Ecuador
mncoro@[Link], esguayasamin@[Link]
Abstract [ CITATION ING16 \l 3082 ] . Con el módulo
In this document you will find the code that links básico de MATLAB el desarrollo de algoritmos a
MATLAB with ANSYS, in addition to performing the
través de un lenguaje propio de programación que
analysis of an exercise with the results shown by
MATLAB, performing the base simulation in ANSYS, resulta ser abierto (integrable con C, Excel y
a series of combination of loads is performed for the Bases de Datos), extensible (a través de las
presentation of results and their analysis. funcionalidades que aportan las librerías
especializadas complementarias) y de sintaxis
Keywords: Deflexión, Tensión similar al C (pero sin las dificultades técnicas de
programación que presenta C). [ CITATION
I. INTRODUCCIÓN Cés02 \l 3082 ]
En el campo de análisis de simulación por
elementos finitos ayuda a obtener resultados de ANSYS es un software confiable para realizar
zonas específicas, en algunos de estos casos se simulaciones de ingeniería mecánica, por lo que
generan bases de datos muy extensas que utiliza es importante que los alumnos que cursan esta
gran memoria computacional, la fusión de estos carrera tengan conocimientos básicos de esta, para
dos programas permite el análisis y comparación realizar investigaciones se requiere mucho de un
de múltiples casos, discretizando los resultados software de simulación y de diferentes análisis
para su posterior análisis. que pueden llevar algún tiempo para realizarlos, lo
Se va a usar Ansys Mechanical debido a que que nos lleva a integrar estas dos potentes
permite resolver situaciones muy particulares en herramientas computacionales para reducir el
diseño mecánico como elementos estructurales tiempo en la generación de resultados.
lineares, no lineares y análisis dinámico.
II. METODOLOGÍA Y PROCEDIMIENTOS
[ CITATION 3DC19 \l 3082 ] . ANSYS usa una
combinación de software robusto y pruebas físicas Generación de geometría en ANSYS APDL
en el sitio para ofrecer servicios de simulación y La geometría a analizar consiste en una cercha
modelado líderes en la industria. Con servicios de metálica generada por siete nodos y once
laboratorio de pruebas físicas, validan y afinan las elementos, cada nodo es común para tres
simulaciones para asegurarse de que los modelos elementos (Figura 1)[ CITATION TCh02 \l
están produciendo resultados fiables. Su equipo
también puede proporcionar orientación sobre qué 3082 ].
modelar, cómo modelarlo y qué significan los
resultados. [ CITATION ANS19 \l 3082 ]
Traducción realizada con el traductor
[Link]/Translator MATLAB resulta
ser una plataforma muy confiable para resolver
problemas de ingeniería y científicos. El lenguaje
de MATLAB se basa en matrices, para expresar
las matemáticas computacionales. Los gráficos
facilitan la visualización de los datos y la
obtención de información a partir de ellos.
obtendrán los resultados de los valores
máximos y mínimos de deformación al igual
que los valores de tensión axial para ser
analizados.
Además, implementando el código en
MATLAB, se obtendrán todos los resultados
de ANSYS sin necesidad de abrir el programa,
mediante un archivo txt, generado por
MATLAB.
Código de MATLAB usado
Figura 1 Generación de geometría
clear all
clc
Definición de condiciones de contorno
! "C:\Program Files\ANSYS
Con la geometría establecida, el nodo uno, sus Inc\v180\ansys\bin\winx64\[Link]" -b -i
desplazamientos estarán restringidos en las tres P10-CORO-GUAYASAMIN-360-280-280-
direcciones X, Y, Z, y el nodo 7 se restringe [Link] -o [Link];
los desplazamientos en la dirección Y.
III. RESULTADOS
Propiedades del material En el presente documento se presentarán graficas
Se le asignará un área constante de 3250 mm 2, de los resultados con las cargas mínimas y
además un módulo de elasticidad de 200GPa. máximas, obtenidas de ANSYS, además de un
cuadro de resumen con todos los resultados
Aplicación de cargas obtenidos con las aplicaciones de cargas obtenidos
mediante un archivo txt generado por MATLAB.
Las cargas se aplicarán en los nodos 1, 3, 5, y 7
iniciando con valores de 360 kN en los Figura 3 Código en MATLAB
extremos (Nodo 1 y 7) y 280 kN en los nodos
intermedios (Nodo 3 y 5), se realizarán
variaciones de carga en el rango de [360, 810]
kN y [280, 730] kN, cada uno con variaciones
de 50kN.
Figura 2 Aplicación de cargas Figura 4 Formato de salida generado por MATLAB
Análisis de resultados
Se correrá el programa en ANSYS para
obtener los resultados gráficos de código, y se
Figura 5 Deflexión con cargas de 360 kN y 280 kN Figura 8 Deflexiones con cargas de 810kN y 730kN
Figura 6 Reacciones en los apoyos, Cargas 360kN y 280
kN
Figura 9 Reacciones en los apoyos, Cargas de 810kN y
730kN
Figura 7 Tensiones en cada elemento, con cargas de
360 kN y 280kN
Se registra una deflexión máxima de 8.46 mm.
Tabla 2 Resultados Tensión máxima y mínima
Tensión Tensión
máxima mínima
Carga
Elemento Elemento
3 1
Nodo 1, Nodo
kN/mm2 kN/mm2
7 3,5
360 280 99,48 -99,48
410 330 117,24 -117,24
460 380 135,01 -135,01
510 430 152,77 -152,77
560 480 170,54 -170,54
610 530 188,30 -188,30
660 580 206,06 -206,06
710 630 223,83 -223,83
760 680 241,59 -241,59
Figura 10 Tensiones en cada elemento, cargas de 810 730 259,36 -259,36
810kN y 730kN
Se registra una deflexión de 22.06 mm
Se despreciarán las reacciones del eje x en el nodo
uno ya que son valores muy bajos.
Tabla 3 Resultados de tensiones en cada elemento
(Primera parte)
Combinación Elemento
de cargas
Tabla 1 Resultados de las reacciones
360/280 410/330 460/380 510/430 560/480
Combinació Reacciones
n de cargas Nodo 1 Nodo 7
(kN) Fy(kN) Fy (kN) Tensión Tensión
Tensión Tensión Tensión
(kN/mm2 (kN/mm2
360/280 640,00 640,00 (kN/mm2) (kN/mm2) (kN/mm2)
) )
410/330 740,00 740,00
1 -99,479 -117,24 -135,01 -152,77 -170,54
460/380 840,00 840,00
2 49,736 58,617 67 76,38 85,262
510/430 940,00 940,00
3 99,479 117,24 135,01 152,77 170,54
560/480 1040,00 1040,00
4 -99,472 -117,23 -135 -152,76 -170,52
610/530 1140,00 1140,00
5 6,50E-15 2,75E-14 6,70E-14 -5,39E-14 -2,07E-15
660/580 1240,00 1240,00
6 99,472 117,23 135 152,76 170,52
710/630 1340,00 1340,00
7 -7,06E-14 -1,57E-14 1,76E-14 -1,69E-14 -5,76E-14
760/680 1440,00 1440,00
8 -99,472 -117,23 -135 -152,76 -170,52
810/730 1540,00 1540,00
9 99,479 117,24 135,01 152,77
Ten sió n -Carga
10 49,736 58,617 67,499 76,38 Elem en to 1
11 -99,479 -117,24 -135,01 -152,77
-90.00
300 400 500 600 700 800
-110.00
Tabla 4 Resultados de tensiones en cada elemento -130.00
(Segunda parte) -150.00
-170.00
Combinación
-190.00
de cargas
-210.00
610/530 660/580 710/630 760/680 -230.00
-250.00
-270.00
Elemento
Tensión Tensión
Tensión Tensión
(kN/mm2 (kN/mm2 Figura 12 Tensión mínima-Carga
(kN/mm2) (kN/mm2)
) )
La tensión máxima del elemento 3 incrementa
1 -188,3 -206,06 -223,83 -241,59 hasta llegar a valor de 259,36 kN/mm 2 con una
2 94,143 103,02 111,91 120,79 deflexión de 22.06mm y registra una tención
3 188,3 206,06 223,83 241,59 mínima en el elemento 1de -99,48 kN/mm 2 y una
deflexión de 8.46mm.
4 -188,29 -206,05 -223,81 -241,57
5 2,51E-14 6,46E-14 4,24E-14 -1,96E-13 IV. ANÁLISIS DE RESULTADOS
6 188,29 206,05 223,81 241,57
7 -9,21E-14 -3,41E-14 -9,33E-14 -1,22E-13 Se observa que para el elemento 3 que presenta
8 -188,29 -206,05 -223,81 -241,57 las tensiones máximas, el valor de las mismas
9 188,3 206,06 223,83 241,59 aumenta al aumentar la carga, mientras que en el
elemento 1, mientras se aumenta la carga el valor
10 94,143 103,02 111,91 120,79
de las tensiones decrece.
11 -188,3 -206,06 -223,83 -241,59
En los resultados obtenidos se puede observar que Los valores de las reacciones en los nodos 1 y 7
las tensiones máximas y mínimas tienen el mismo presentan un comportamiento típico de un cuerpo
valor con signos distintos, además se observa un con dos apoyos, dividiendo la carga de forma
aumento en el valor de la tensión que es equitativa en cada uno, presentan los mismos
directamente proporcional al aumento de carga. valores.
En cuanto a las deflexiones se tiene un valor
Te n sión -Carg a mínimo de 8.46 y un máximo de 22.06, utilizando
Ele m en to 3 la relación de L/360 para la longitud de vigas que
270.00 se está analizando es admisible una deflexión
250.00 máxima de 12mm, que cumple en el caso del
230.00 valor mínimo, pero no para el caso de deflexión
210.00 máxima.
190.00
170.00 V. CONCLUSIONES
150.00 Utilizando variaciones de cargas con
130.00 intervalos de 50 kN, desde 360 a 810 kN
110.00 en los extremos y de 280 a 730 kN en
90.00 los nodos medios, los elementos críticos
300 400 500 600 700 800
serán el elemento 1 y elemento 3 debido
a que están cercanos al nodo 1 que tiene
Figura 11 Tensión máxima –Carga
restricción de desplazamiento en todas
las direcciones.
Las cargas son distribuidas
equitativamente para cada nodo 1 y 7 de
apoyos, obteniendo valores crecientes
desde 640 kN a 1540 kN que van
incrementándose en intervalos de 100
kN, resultados que son los esperados
debido a la distribución de cargas del
sistema.
La deflexión máxima obtenida con la
primera combinación de cargas validan el
modelo para su aplicación obteniendo un
valor inferior a 12mm, sin embargo, para
la última combinación de cargas el
modelo ya no es válido ya que supera al
valor permitido con 10mm.
Para lograr la correcta interfaz entre
MATLAB y ANSYS es necesario que
los dos programas a ejecutar se
encuentren en el mismo directorio, de lo
contrario MATLAB no podrá encontrar
el archivo para ejecutar el comando.
VI. REFERENCIAS
3D CAD Portal . (2019). Primer portal
CAD CAM CAE en español . Obtenido
de
[Link]
ANSYS. (2019). Explore Pervasive
Engineering Simulation. Obtenido de
[Link]
INGELSI. (2016). MATLAB El lenguaje
del cálculo técnico . Obtenido de
[Link]
Pérez, C. (2002). MATLAB y sus
aplicaciones en las ciencias y la
ingeniería. Madrid: Pearson Educación,
S.A.
T, C., & Belegundu , A. (2002).
Introducción al elemento finito en
ingenieria . New Jersey : Prentice Hall.
