LC1 MC516
LC1 MC516
LC1 MC516
2022-II
ÍNDICE
1. Introducción
2. Enunciado del problema
3. Solución analítica
4. Resultados analíticos
5. Solución manual
6. Solución con Matlab
7. Solución con ANSYS
8. Conclusiones
9. Bibliografía
1. Introducción
El presente informe tiene como finalidad dar a conocer de manera detallada tres diferentes formas de solución
de un problema basado en la teoría y los principios del cálculo por elementos finitos, en donde se analizarán los
esfuerzos y deformaciones que surgen debido a la aplicación de una carga, así como la fuerza reactiva generada por
dicha carga.
El primer método consistirá en aplicar el método débil paso a paso, valiéndose de las leyes y métodos en las
que se basa dicha formulación, tanto de forma escalar como de forma matricial. De esta manera se obtendrán
resultados que luego serán comparados con los obtenidos por los softwares de Matlab y ANSYS.
En el segundo método se obtendrán resultados más rápidos y precisos que serán calculados por medio del
software de Matlab, corroborando así los resultados obtenidos en el primer método.
Finalmente, en el tercer método se mostrarán los pasos detallados para realizar la simulación del problema por
medio del software ANSYS 2022-R2.
2. Enunciado del problema
Se tiene una barra con una sección transversal
variable que soporta una carga , tal como se muestra en
la figura 1.1. La barra está fija en un extremo y lleva la
carga P en el otro extremo. Se designa el ancho de la
barra en la parte superior por , en la parte inferior por ,
un espesor de y longitud . El módulo de elasticidad de
la barra se denotará por . Se desea aproximar la
deformación de la barra en varios puntos a lo largo de
su longitud cuando esté sujeta a la carga P. Asumiendo
que la carga aplicada es considerablemente mayor que
el peso de la barra, en el análisis se despreciará el peso
de la barra.
Figura 1.1
3. Solución analítica
3.1. Discretización en elementos finitos
El primer paso consiste en subdividir la placa en
nodos y elementos, para este modelo se usarán cinco
nodos y cuatro elementos, tal como se muestra en la
figura 1.2. La placa dada se modela usando cuatro
segmentos individuales (elementos), cada uno de los
cuales tiene una sección transversal uniforme. El área
de la sección transversal de cada elemento está
representada por un área promedio de las secciones
transversales en los nodos que componen (definen) el
elemento.
Figura 1.2
3.2. Aproximación del comportamiento de cada
elemento
Para estudiar el comportamiento de un elemento
típico, como el de la figura 1.3 (a), se considera la
deflexión de un miembro sólido con una sección
transversal uniforme que tiene una longitud cuando se
somete a una fuerza .
Para este elemento, se conoce que su esfuerzo
promedio “” y deformación unitaria normal promedio
“” se definen como:
Reordenando las ecuaciones de equilibrio dadas en forma matricial y separando las fuerzas de reacción y
externa aplicada de las fuerzas internas, se tiene que la matricial de rigidez global es:
3.5. Condición de contorno
Debido a que la barra está fija en la parte superior, el desplazamiento del nodo es cero. Por lo tanto, solo hay
cuatro valores de desplazamiento nodal desconocidos, , , y . La fuerza de reacción en el nodo , , también se
desconoce; en total, hay cinco incógnitas. Debido a que hay cinco ecuaciones de equilibrio, se debería poder
resolver todas las incógnitas. Sin embargo, para eliminar la necesidad de considerar simultáneamente la fuerza de
reacción desconocida y poder centrarse primero en los desplazamientos desconocidos, se utiliza la condición de
frontera conocida y se reemplaza la primera fila de la matriz de rigidez global con una fila que dice .
3.6. Fase de solución algebraica
Para la solución numérica de los desplazamientos nodales, en el problema se van a considerar los siguientes
datos: (aluminio), , , , y .
De estos datos, se desprende la siguiente tabla aplicada para cada elemento y nodo:
La variación del área de la sección transversal de la barra en la dirección se puede expresar como:
Usando esta ecuación, se puede calcular que el área transversal en cada nodo es:
A continuación, el coeficiente de rigidez equivalente para cada elemento se calcula a partir de las
ecuaciones:
A partir de esto, las matrices elementales son:
Ensamblando las matrices elementales, se tiene que la matriz de rigidez global es:
Aplicando la condición de frontera y la carga , se obtiene:
Otra forma es a partir del equilibrio estático de toda la barra, la cual requiere que:
Finalmente, para la obtención de los esfuerzos normales promedio en cada elemento, se pueden determinar a
partir de la ecuación:
Empleando la ecuación, la tensión normal promedio para cada elemento es:
Por otro lado, para el problema dado, independientemente de dónde se corte una sección a través de la barra, la
fuerza interna en la sección es igual a , por ende:
4. Resultados analíticos
Los valores obtenidos para la reacción y los esfuerzos debido a la carga aplicada son:
5. Solución manual:
6. Solución con Matlab
Compilando el Script:
7. Solución con ANSYS
8. Conclusiones
• Un modelo con una mayor cantidad de nodos y elementos adicionales a los considerados en el problema,
hubiese aumentado la precisión de los resultados obtenidos.
• Existe más de una forma posible de hallar las deformaciones de los nodos, los esfuerzos promedios normales de
los elementos y la fuerza reactiva de la barra.
• Matlab ayuda a resolver el problema con facilidad; debido a que, se le puede indicar la cantidad de elementos
que se va a utilizar, y modificar los datos del problema.
9. Bibliografía
• Saeed Moavenni, “Teoría y aplicación de análisis por elementos finitos”, cuarta edición.
• Galeas R. A. & Galeas J. A. MATLAB Y SIMULINK Aplicado a la ingeniería (2017). Primera edición.