Source: http://docplayer.es/12153090-Calculo-de-placas-a-traves-de-distintas-metodologias.html
Timestamp: 2018-12-17 12:52:15+00:00

Document:
CÁLCULO DE PLACAS A TRAVÉS DE DISTINTAS METODOLOGÍAS - PDF
Download "CÁLCULO DE PLACAS A TRAVÉS DE DISTINTAS METODOLOGÍAS"
Esperanza Lozano Hernández
1 UNIVERSIDAD CARLOS III DE MADRID ESCUELA POLITÉCNICA SUPERIOR Departamento de Mecánica de Medios Continuos y Teoría de Estructuras Ingeniería Técnica Industrial Mecánica Proyecto Fin de Carrera CÁLCULO DE PLACAS A TRAVÉS DE DISTINTAS METODOLOGÍAS Autor: Miguel Ortega Oyarzábal Tutor: Sergio González López
2 Octubre, 2013 Índice Introducción y Objetivos Introducción Objetivos... 1 Introducción a Placas Placas, Estructuras tipo placa y Clasificación Hipótesis básicas del cálculo de placas Proceso de cálculo Integración de tensiones Ecuaciones Cinemáticas (Compatibilidad) Ecuaciones Constitutivas Ecuaciones de equilibrio Métodos de Cálculo Método de Navier Métodos numéricos Método de los Elementos Finitos Método de las Diferencias Finitas Casos y Resultados en Estático Resolución analítica. Método de Navier Carga Uniforme Carga Puntual Carga Triangular Resolución mediante el Método de Elementos Finitos (M.E.F.) Carga Uniforme Carga Triangular Carga Puntual Resolución mediante diferencias finitas Carga Uniforme Carga Triangular Carga Puntual Introducción a la Dinámica de Placas i
3 5.1. Introducción a la dinámica Casos y Resultados en Dinámico Resolución Analítica Carga Uniforme Carga Puntual Carga Triangular Resolución Mediante M.E.F Carga Uniforme Carga Triangular Carga Puntual Resolución mediante Diferencias Finitas Método Dinámico Método Dinámico Estático Comparativa Comparativa de resultados en estático Carga Uniforme Carga Puntual Carga Triangular Comparativa de resultados en dinámico Carga Uniforme Carga Triangular Carga Puntual Conclusiones Conclusiones Trabajos Futuros Bibliografía Anexos Anexo A Guía de uso de Abaqus Anexo B Tabla con resultados numéricos de la flecha de una placa sometida a carga triangular ii
4 Índice de Figuras Fig 1. Ejemplo de placa (izquierda) y de laja (derecha)... 3 Fig 2. Forjado de un edificio... 4 Fig 3. Cubierta plegada... 4 Fig 4. Losa de puente Fig 5. Aproximación a lámina Fig 6. Momentos Flectores creados por x y y... 8 Fig 7. Momentos torsores creados por las tensiones xy y yx... 9 Fig 8. Cortantes creados por las tensiones xz y yz Fig 9. Desplazamiento en el eje X Fig 10. Desplazamiento en el eje Y Fig 11. Esfuerzos en un elemento diferencial de placa Fig 12. Borde empotrado Fig 13. Borde apoyado Fig 14. Borde libre en x=a Fig 15. Métodos de resolución de placas Fig 16. Placa apoyada en sus cuatro lados Fig 17. División de un cuerpo unidimensional mediante nodos Fig 18. Ejemplo de malla bidimensional y tridimensional Fig 19. Grados de libertad en un elemento rectangular Fig 20. Ejemplo de malla mixta Fig 21. Diferentes tipos de malla Fig 22. Aproximación de un contorno irregular mediante una malla Fig 23. Aproximación de una función mediante diferencias finitas Fig 24. Molécula computacional con los coeficientes de cada elemento Fig 25. Borde empotrado Fig 26. Borde apoyado Fig 27. Puntos ficticios en borde libre Fig 28. Borde libre Fig 29. Placa apoyada en sus cuatro lados sometida a carga uniforme Fig 30. Placa apoyada en sus cuatro lados sometida a una carga puntual Fig 31. Placa apoyada en sus cuatro lados sometida a una carga triangular Fig 32. Gráfico de la flecha total en la placa sometida a carga triangular Fig 33. Gráfico con los resultados para m = n = iii
5 Fig 34. Gráfico para el resto de combinaciones de valores de m y n Fig 35. Viga sometida a carga triangular Fig 36. Imagen de Placa Apoyada con Carga Repartida Constante Fig 37. Gráfico Tensión de V. Mises para análisis de sensibilidad de malla Fig 37. Gráfico de la Flecha para análisis de sensibilidad de malla Fig 38. Mallado de 400 elementos Fig 39. Placa deformada Fig 40. Placa con carga triangular Fig 41. Gráfico de la Tensión de V. Mises para análisis de sensibilidad de malla Fig 42. Gráfico de la Flecha para análisis de sensibilidad de malla Fig 43. Placa deformada bajo carga triangular Fig 44. Placa con Carga puntual en dirección negativa al eje Z Fig 45. Gráfico de la Flecha para análisis de sensibilidad de malla Fig 46. Gráfico de la Tensión de V. Mises para análisis de sensibilidad de malla Fig 47. Deformada Carga Puntual Fig 48. Zonas de la malla para el cálculo de la placa Fig 49. Gráfica de iteraciones vs. deformación para carga uniforme Fig 50. Deformación de la placa bajo carga uniforme Fig 51. Comprobación de la solución para carga uniforme Fig 52. Gráfica de iteraciones vs. deformación para carga triangular Fig 53 Placa deformada bajo carga triangular Fig 54. Deformada de la placa bajo carga triangular Fig 55. Comprobación de la solución para carga triangular Fig 56. Gráfico de iteraciones vs deformación para carga puntual Fig 57. Deformada bajo carga puntual Fig 58. Imagen de la matriz de comprobación para espesor de malla 0,2 m Fig 59. Gráfico de iteraciones vs deformación para carga puntual Fig 60. Deformada bajo carga puntual Fig 61. Imagen de la matriz de comprobación para espesor de malla 0,1 m Fig 62. Gráfico de iteraciones vs deformación para carga puntual Fig 63. Deformada bajo carga puntual Fig 64. Imagen de la matriz de comprobación para espesor de malla 0,05 m Fig 65. Placa sometida a una fuerza transversal F(x,y,t) Fig 66. Placa apoyada en sus cuatro lados sometida a carga uniforme Fig 67. Deformación en función del tiempo para el primer término de Galerkin Fig 68. Deformación en función del tiempo para el tercer término de Galerkin iv
6 Fig 69. Deformación en función del tiempo total Fig 70. Placa apoyada en sus cuatro lados sometida a una carga puntual Fig 71. Deformación en función del tiempo para carga puntual Fig 72. Placa apoyada en sus cuatro lados sometida a una carga triangular Fig 73. Deformación en función del tiempo para carga triangular Fig 74. Imagen de Placa Apoyada con Carga Repartida Constante Fig 75. Análisis dinámico carga uniforme 3 segundos Fig 76. Placa con carga triangular Fig 77. Análisis dinámico carga triangular 3 segundos Fig 78. Placa con Carga puntual en dirección negativa al eje Z Fig 79. Análisis dinámico carga puntual 3 segundos Fig 80. Esquema desarrollo diferencias finitas dinámica Fig 81. Deformación punto medio de la placa con Δt=0,001 s Fig 82. Deformación punto medio de la placa con Δt=0,0001 s y Δt=0,00001 s Fig 83. Deformación punto medio de la placa con Δt=0,00005 s y Δt=0,00001 s Fig 84. Deformación en el punto de máxima flecha de la placa con Δt=0,00005 s y Δt=0,00001 s Fig 85. Esquema desarrollo método combinado dinámico estático Fig 86. Deformación del punto medio de la placa Fig 87. Deformación del punto medio de la placa con Δt=0,001 s. Detalle Fig 88. Deformación del punto medio de la placa con Δt=0,0001 s y Δt=0,00001 s. 126 Fig 89. Deformación del punto de máxima flecha de la placa. Y detalle Fig 90. Deformación en el punto de máxima flecha de la placa con Δt=0,0001 s Fig 91. Deformación en el punto de máxima flecha de la placa con Δt=0,00001 s y Δt=0, s Fig 92. Deformación en el punto de máxima flecha de la placa con Δt=0,001 s Fig 93. Deformación del punto medio de la placa con Δt=0,0001 s y Δt=0,00001 s Fig 94. Gráficas obtenidas mediante (de izquierda a derecha y de arriba abajo) método analítico, M.E.F., diferencias finitas (solo dinámico Δt=0,00001 s) Fig 95. Gráficas obtenidas mediante (de izquierda a derecha y de arriba abajo) método analítico, M.E.F., diferencias finitas (solo dinámico Δt=0,00001 s) Fig 96. Gráficas obtenidas mediante (de izquierda a derecha y de arriba abajo) método analítico, M.E.F., diferencias finitas (solo dinámico Δt=0,00001 s) Fig.1. Pantalla Principal de Abaqus Aspecto de la placa después de aplicarle las condiciones de contorno y la carga Cuadro de dialogo para seleccionar la forma de los elementos Cuadro de diálogo para definir el tipo de elemento v
7 Resultado del mallado de la placa Placa deformada con mapa de colores y presentando el máximo y mínimo valor de los movimientos vi
8 Capítulo 1 Introducción y Objetivos 1.1 Introducción A lo largo de la historia las placas se han convertido en un tipo de estructura fundamental en la ingeniería y su empleo es muy amplio en campos como la ingeniería civil. Algunas de sus funciones más usuales son como estructura de cubierta, en muros de contención o en losas de cimentación. También han sido el tipo estructural por antonomasia en los tableros de puentes, en los cuales se ha usado ampliamente. Dada la importancia de las placas como elemento estructural el análisis del comportamiento de las placas es realmente importante. La forma tradicional de solucionar los problemas de placas es mediante el uso de desarrollos en series de Fourier, los cuales están ligados a las condiciones de contorno y al sistema de coordenadas elegido. Actualmente los cálculos de placas se realizan mayormente mediante métodos numéricos como los elementos finitos, o las diferencias finitas. Estos métodos permiten encontrar una solución a problemas que resolver de la manera tradicional sería muy difícil o incluso imposible. Sin embargo, para poder usar estos métodos se hace necesario el uso de ordenadores. En el presente PFC la placa delgada será sometida a diferentes tipos de carga (manteniendo la condiciones de contorno) constante en el tiempo. Y los resultados se analizarán desde el punto de vista estático y dinámico Objetivos En este PFC, puesto que se divide en dos partes bien diferenciadas (estática y dinámica), se puede hablar de dos objetivos principales. Estos objetivos son desarrollar el procedimiento de obtención de la flecha en una placa sometida a cargas perpendiculares de manera estática y dinámica. 1
9 Para cumplir los dos objetivos principales se han de completar los siguientes objetivos secundarios: Obtención de la ecuación diferencial de la flexión de placas. Ecuación de Lagrange. Desarrollo de los métodos empleados en el cálculo de la flecha. Método de Navier, Elementos Finitos y Diferencias Finitas.El método de Navier ofrece una solución analítica. Los otros dos métodos son métodos numéricos. Para la resolución mediante el método de los Elementos finitos se empleará el programa Abaqus y para las Diferencias Finitas se usará como apoyo el programa Matlab. Inclusión de la variable temporal en la ecuación de Lagrange. Aplicación de los métodos analítico (Galerkin), Elementos Finitos (Abaqus) y Diferencias Finitas (Matlab) en dinámica. Solución de los casos planteados (carga uniforme, triangular y puntual), tanto en estática como en dinámica. Finalmente, una vez se han obtenido todas las soluciones, estas se analizan y comparan entre sí. De esta manera se pueden comparar los diferentes métodos empleados en los análisis. 2
10 Capítulo 2 Introducción a Placas 2.1 Placas, Estructuras tipo placa y Clasificación Se denomina placa al elemento estructural plano en el que una dimensión (el espesor) es muy pequeño en comparación con las otras dos dimensiones. Clásicamente se ha venido diferenciando entre placa y laja según la dirección de las fuerzas que actúan sobre ellas. En la placa las fuerzas son normales a su plano medio, mientras que en la laja están contenidas en él. La placa también puede estar sometida a momentos cuyos ejes estén contenidos en el plano medio. Se define plano medio de la placa como aquel que divide el espesor en dos partes iguales. Fig 1. Ejemplo de placa (izquierda) y de laja (derecha). Estructuras reales tipo placa El elemento estructural placa se encuentra muy a menudo en las estructuras reales. Puede aparecer solo como placa, como en los forjados de un edificio, o como una combinación de placa y laja (p.e.: cubiertas plegadas, losas de 3
11 puentes) o en una aproximación a lámina (análogo a la sustitución de un arco por una poligonal). A continuación se muestran unos ejemplos. Fig 2. Forjado de un edificio Fig 3. Cubierta plegada 4
12 Fig 4. Losa de puente. Fig 5. Aproximación a lámina. Clasificación de las placas Las placas se pueden clasificar de diferentes maneras: Según el espesor: o Delgadas: El espesor es menor que la quinta parte de la menor de sus otras dimensiones. o Gruesas: El espesor es mayor que la quinta parte de la menor de sus otras dimensiones. 5
13 En las placas delgadas se supone la deformación por cortante despreciable. Algunos autores denominan placa delgada a aquella en la el espesor es 20 veces menor que la más pequeña de sus otras dimensiones. Según sus características elasto mecánicas: o Isótropas: Las características no dependen de la dirección (igual en cualquier dirección). o Anisótropas: Las características dependen de la dirección. o Ortótropas: Las características dependen de la dirección presentando dos direcciones ortogonales principales. Para el cálculo de placas se tienen en cuenta una serie de hipótesis de partida Hipótesis básicas del cálculo de placas A la hora de calcular placas, se estudiarán placas delgadas e isótropas. Además se tendrán en cuenta unas hipótesis de partida. Estás hipótesis de partida (hipótesis de Kirchoff) facilitan el cálculo de las placas. 1. El material es elástico y lineal. 2. El espesor es mucho más pequeño que las otras dos dimensiones de la placa. 3. La placa inicialmente es plana. 4. Las deformaciones son pequeñas en comparación con el espesor de la placa. 5. Los puntos del plano medio solo tiene movimientos perpendiculares a dicho plano. 6. Los puntos situados en una perpendicular al plano medio antes de la deformación siguen situados en una recta perpendicular al plano medio después de la deformación. 7. Las tensiones normales al plano medio se consideran despreciables. 8. Los esfuerzos en el plano medio producidos por fuerzas en ese plano, se consideran despreciables. 6
14 2.3. Proceso de cálculo Al igual que se hace para el cálculo de vigas, en el cálculo de placas se integran las tensiones en el espesor y se trabaja con esfuerzos lo que permite reducir el problema a un estado bidimensional en vez de uno tridimensional. Datos: Se suponen conocidas las características elasto geométricas de la placa así como las cargas aplicadas. Incógnitas: o Movimientos (sólo se tendrán en cuenta los movimientos del plano medio). o Deformaciones o Esfuerzos (se obtienen integrando las tensiones). Ecuaciones: o Cinemáticas: Relacionan los movimientos con las deformaciones. o Constitutivas: Relacionan los esfuerzos con las deformaciones. o Equilibrio: Relacionan las cargas con los esfuerzos Integración de tensiones Integrando las tensiones se obtienen los esfuerzos por unidad de longitud. Los esfuerzos que pueden aparecer en una placa delgada sometida a una carga transversal son momentos flectores, esfuerzos cortantes y momentos torsores. Así, un elemento diferencial de placa de tamaño dxxdy está sometido a tensiones tal y como se representa en la figura. Donde x y y crean los momentos flectores m x y m y respectivamente. 7
15 Fig 6. Momentos Flectores creados por x y y. Los momentos m x y m y se expresan de la siguiente manera: (2.1) (2.2) Los momentos torsores m xy y m yx son creados por las tensiones xy y yx respectivamente. 8
16 Fig 7. Momentos torsores creados por las tensiones xy y yx. (2.3) (2.4) Debido a las hipótesis de partida xy = yx, de donde m xy = m yx. Y por último, los cortantes q x y q y los crean las tensiones xz y yz respectivamente. 9
17 Fig 8. Cortantes creados por las tensiones xz y yz. (2.5) (2.6) Ecuaciones Cinemáticas (Compatibilidad) Estas ecuaciones relacionan deformaciones y movimientos. Se Estudiará el desplazamiento debido a flexión que tienen en la placa dos puntos A y B separados por una distancia z perpendicular al plano medio de la placa. El punto A pertenece al plano medio, por lo que solo tiene movimiento vertical. La perpendicular al plano medio permanece normal al plano medio deformado. 10
18 Fig 9. Desplazamiento en el eje X. Así se puede ver que ambos puntos se desplazan una distancia w en la dirección del eje Z. También se observa que el punto B sufre un pequeño desplazamiento u en sentido contrario al eje X. Este desplazamiento se expresa de la siguiente manera: (2.7) La deformación en la dirección X resultaría: (2.8) 11
19 Fig 10. Desplazamiento en el eje Y. Igual que se ha hecho en el eje X se hace en el eje Y. En este caso en punto sufrirá un desplazamiento w en el eje Zy un desplazamiento v en el eje Y. (2.9) (2.10) 12
20 La deformación en el plano XY también es posible expresarla en función de las derivadas de la flecha. (2.11) Si se denomina curvaturas a: 2 2 w x 2 y 2 x w y xy 2 w x y (2.12) Se puede expresar las deformaciones de la siguiente manera: x z x y z y xy 2 z xy (2.13) Ecuaciones Constitutivas Estas ecuaciones relacionan las tensiones y las deformaciones. Se asume que el material es elástico por lo que se puede usar la ley de Hooke de manera que se obtienen las tensiones en función de las deformaciones. ( ) (2.14) ( ) (2.15) (2.16) También se puede escribir las tensiones en función de las curvaturas: 13
21 ( ) (2.17) ( ) (2.18) (2.19) Donde E es el módulo de elasticidad y es el módulo de Poisson. Integrando las tensiones para obtener los esfuerzos: ( ) ( ) (2.20) La solución de la integral es: (2.21) Con lo que el momento m x queda expresado: ( ) (2.22) Donde (2.23) 14
22 que representa la rigidez. Y m x quedaría: ( ) (2.24) Procediendo igual con los otros momentos (m y y m xy ) queda: ( ) (2.25) (2.26) Ecuaciones de equilibrio Ahora se plantearan las ecuaciones de equilibrio de los esfuerzos en un diferencial de placa de dimensiones dxxdy y espesor h como el de la figura. Fig 11. Esfuerzos en un elemento diferencial de placa. 15
23 Considerando el equilibrio de fuerzas en el eje Y se obtiene: ( ) ( ) (2.27) ( ) Simplificando: (2.28) Se desprecia el término de tercer orden: (2.29) Y finalmente: (2.30) Operando igual en el eje X: ( ) ( ) (2.31) ( ) Simplificando y despreciando los términos de tercer orden: (2.32) 16
24 Y finalmente el equilibrio en el eje Z: ( ) ( ) (2.33) Simplificando: (2.34) (2.35) Ahora se sustituyen las expresiones de los momentos obtenidas anteriormente en las ecuaciones de equilibrio de los ejes X e Y. [ ( )] [ ] (2.36) [ ] (2.37) Simplificando: [ ] ( ) (2.38) Procediendo igual con la otra expresión de cortante (q y ): 17
25 [ ] [ ( )] (2.39) [ ] (2.40) Simplificando: [ ] ( ) (2.41) Y sustituyendo estas dos expresiones en la del equilibrio en el eje Z nos queda: [ ( )] [ ( )] (2.42) ( ) (2.43) (2.44) Usando el operador Laplaciano : (2.45) Esta es la ecuación diferencial de la flexión de placas. Esta ecuación fue obtenida por Lagrange y es la base de la teoría de placas. 18
26 El problema queda reducido a encontrar una flecha w(x,y) que cumpla la ecuación de Lagrange y las condiciones de contorno. Una vez conocida w(x,y) se pueden averiguar los esfuerzos con las expresiones obtenidas anteriormente en que los esfuerzos estaban en función de las derivadas de la flecha. A continuación se verán algunas de las condiciones de contorno más usuales en el cálculo de placas. 1. Borde empotrado Cuando un lado de una placa está empotrado no hay movimiento alguno a lo largo del mismo (ni desplazamiento, ni giro). Por lo que, si el lado empotrado de la placa es x=a, se cumplirá: (2.46) (2.47) y x=a x Fig 12. Borde empotrado. 19
27 2. Borde apoyado En este caso, la condición implica que el desplazamiento y el momento en el lado apoyado son nulos. (2.48) (2.49) y x=a x Fig 13. Borde apoyado. 3. Borde Libre El borde libre conlleva que los esfuerzos en ese borde son nulos. Por tanto: (2.50) (2.51) 20
28 (2.52) y x Fig 14. Borde libre en x=a. Según Kirchoff las dos últimas condiciones se pueden agrupar en una. Representando el momento torsor m xy como pares de fuerzas sobre la cara lateral de la placa se puede definir un cortante adicional que es equivalente a la variación del momento torsor a lo largo del borde. (2.53) La suma de este cortante adicional más el cortante q x es lo que se conoce como reacción de borde v x. [ ] [ ] (2.54) Con condición de contorno de borde libre esta reacción ha de ser nula. De esta manera: [ ] [ ( )] (2.55) 21
29 [ ] [ ( )] (2.56) 22
30 Capítulo 3 Métodos de Cálculo En este capítulo se verán diferentes métodos de cálculo para las placas. Hay muchos métodos distintos, algunos de los más usuales se reflejan en el siguiente esquema. 23
31 Fig 15. Métodos de resolución de placas 24
32 En este PFC se desarrollaran los cálculos analíticamente con el método de Navier, y numéricamente con los métodos de elementos finitos y de diferencias finitas Método de Navier Este método permite obtener la solución para la ecuación diferencial de la flexión de placas en el caso de tener los cuatro lados apoyados. Se basa en la aplicación de desarrollos en serie de Fourier. y x Fig 16. Placa apoyada en sus cuatro lados. Para una placa cualquiera de dimensiones axb y de espesor h, apoyada en sus cuatro lados y sometida a una carga p(x,y) la flecha ha de cumplir la ecuación diferencial de las placas. (3.1) La solución propuesta por Navier expresa la carga p(x,y) como una doble serie de Fourier con la siguiente expresión: (3.2) 25
33 Igualmente la flecha w(x,y) se expresa como una doble serie de Fourier: (3.3) (3.4) Esta es la solución a la ecuación de Lagrange. Se puede comprobar que esta solución cumple las condiciones de contorno de cuatro lados apoyados, es decir las flechas en los cuatro lados es nula: w 0 w 0 w 0 w 0 x 0 x a Y los momentos en los apoyos son nulos: y 0 y b m 0 m 0 m 0 m 0 x x 0 x x a y y 0 y y b Introduciendo las expresiones de w(x,y) y p(x,y) en (3.1)se pueden despejar los coeficientes W m,n de la flecha en función de los coeficientes P m,n de la carga. ( ) (3.5) Por tanto primero deben calcularse los coeficientes P m,n de la carga. Para ello hay que integrar la expresión de la carga p(x,y), y se obtiene: (3.6) 26
34 En el desarrollo en serie de Fourier hay que coger conjuntos completos de términos, es decir, suma de M+N iguales, por ejemplo: TÉRMINO VALORES DE M VALORES DE N SUMA DE M+N 1 er Término º Término 3 er Término Tabla 1. Ejemplo de valores de m y n. Una vez obtenida la flecha se pueden determinar los esfuerzos aplicando las ecuaciones obtenidas anteriormente en que se relacionaban la flecha y los esfuerzos Métodos numéricos Ingenieros y científicos a menudo se enfrentan a problemas cuya solución mediante métodos analíticos convencionales es de extrema dificultad o incluso imposible. Por ejemplo, un cuerpo cualquiera tridimensional sobre el que actúan una serie de fuerzas externas. Para poder analizar la respuesta exacta de este cuerpo a esas fuerzas se busca una solución aproximada a las ecuaciones que rigen su deformación. Sin embargo, debido a la complejidad geométrica que por regla general tienen los problemas prácticos es extremadamente difícil y a menudo imposible obtener dicha solución. Ante esto el recurso es encontrar una solución numérica al problema. Hay varios métodos 27
35 numéricos válidos. En este proyecto nos centraremos en el Método de Elementos Finitos y el Método de las Diferencias Finitas Método de los Elementos Finitos El de los elementos finitos es un método relativamente reciente. El método fue desarrollado por M.J. Turner en la década de los 50. Sin embargo se considera el verdadero inicio del M.E.F. actual a partir de 1956 a raíz del estudio llevado a cabo por Turner, R.W. Clough, H.C. Martin y L.J. Topp. Este estudio se centraba en la rigidez y deformación de estructuras complejas. En sus inicios el M.E.F. está muy unido a la industria aeroespacial, pues era la única que se podía permitir los equipos necesarios para los cálculos. Originalmente el M.E.F. se utilizaba solamente para el cálculo de estructuras. Sin embargo, con el posterior desarrollo de los ordenadores el M.E.F. se generalizó y se pudo aplicar a más campos diferentes. Unido al desarrollo de los ordenadores se extendió el uso de programas especializados en diferentes campos. El MEF permite obtener una solución numérica aproximada sobre un cuerpo, estructura o dominio sobre el que están definidas ciertas ecuaciones diferenciales que caracterizan el comportamiento físico del problema dividiéndolo en un número elevado de subdominios denominados «elementos finitos». El M.E.F. facilita el cálculo en problemas en que es prácticamente imposible encontrar una solución analítica. Actualmente, gracias a la generalidad del M.E.F. se puede aplicar en una gran cantidad de campos diferentes desde estudios estructurales a problemas termodinámicos o incluso magnéticos. Y se usa en muchos y diversos campos industriales (aviación, automoción, ) o de la construcción. Hoy en día existen varios programas para el uso del M.E.F. algunos de los más conocidos son: SolidWorks ANSYS Abaqus (empleado en este PFC) Nastran En este método el cuerpo se divide en cierto número de elementos o dimensiones finitas, de ahí su nombre. Si el cuerpo tiene n (n= 1,2,3) dimensiones en el espacio, se dividirá en elementos finitos de n dimensiones. Los cuerpos unidimensionales se dividen en elementos finitos mediante nodos. 28
36 Fig 17. División de un cuerpo unidimensional mediante nodos. Para la división de los cuerpos de dos y tres dimensiones se usan líneas y planos. Fig 18. Ejemplo de malla bidimensional y tridimensional. En los cuerpos unidimensionales los elementos finitos resultantes pueden tener longitudes diferentes, así mismo los elementos finitos en los cuerpos de dos y tres dimensiones pueden tener diferentes tamaño y forma. En todos los casos los elementos finitos estarán conectados por nodos. De esta manera el cuerpo objeto del estudio es sustituido por un sistema de elementos finitos conectados por nodos. Una vez se tiene el cuerpo dividido en elementos finitos el siguiente paso es determinar la matriz de rigidez de los elementos individualmente. Luego estas 29
37 se unen para formar la matriz de rigidez del cuerpo de manera que la continuidad de los movimientos y el equilibrio de fuerzas prevalezca en todos los nodos del modelo. Esto nos lleva a la siguiente ecuación matricial: K P (3.7) Donde [K] es la matriz de rigidez del cuerpo, (P) es el vector de las fuerzas externas aplicadas en todos los nodos y (δ) es el vector de los desplazamientos de los nodos. Para unas condiciones de contorno específicas y unas fuerzas aplicadas sobre el cuerpo esta ecuación se puede resolver hallando los desplazamientos de los nodos. Una vez hallados se pueden calcular las tensiones y los esfuerzos. Como se ha dicho anteriormente [K] se forma a partir de las matrices de rigidez de cada uno de los elementos [K i ] en que se divide el cuerpo objeto del estudio. Para averiguar la [K i ] de un elemento se estudian los desplazamientos en los vértices de ese elemento. En este caso se empleará el elemento rectangular, puesto que va a ser el que se utilice posteriormente para resolver los casos propuestos como ejemplos. Fig 19. Grados de libertad en un elemento rectangular. Cada uno de los vértices del elemento tiene tres grados de libertad, un desplazamiento y dos giros. Como son cuatro vértices se puede definir el desplazamiento w(x,y) como un polinomio definido por 12 constantes. De esta manera, 30
38 (3.8) También se puede expresar así { } { } (3.9) Y { } { } (3.10) De esta manera, en forma matricial { }{ } (3.11) Partiendo de la expresión de w(x,y) como un polinomio se pueden obtener los giros. (3.12) (3.13) (3.14) (3.15) 31
39 Particularizando las expresiones de los desplazamientos para uno de los vértices del elemento se obtiene: (3.16) (3.17) (3.18) Donde el subíndice i indica que pertenece al elemento i. Realizando esto mismo en todos los vértices del elemento, se puede expresar de la siguiente manera: { } [ ]{ } (3.19) Donde { }es el vector de todos los desplazamientos de los vértices del elemento i. { } { } (3.20) La matriz [c] es la matriz de coeficientes y solo depende de las coordenadas de los vértices. 32
40 [ ] [ ] Derivando la expresión de w(x,y) se obtiene: (3.21) { } [ ] { } { } Se puede escribir { } [ ]{ } (3.22) Siendo { } el vector de deformaciones. Se puede escribir la expresión general de la matriz de rigidez de un elemento de la siguiente manera: [ ] [ ] [ ] [ ][ ] [ ] (3.23) Donde [d] es la matriz de elasticidad de la placa: 33
41 [ ] [ ] (3.24) Conviene señalar que para estos cálculos se ha utilizado el centroide del elemento como origen de coordenadas. Una vez obtenidas las matrices de rigidez [ ] de todos los elementos en que se ha dividido la placa hay que ensamblarlas para crear la matriz de rigidez de la placa. Debido al tamaño de matrices que se ha de manejar, el trabajo de ensamblaje de la matriz de rigidez global se realiza mediante ordenadores. Resumiendo, los pasos a seguir para resolver un problema mediante el método de elementos finitos serían los siguientes: 1. División del cuerpo objeto de estudio en un sistema de elementos finitos. 2. Obtención de la matriz de rigidez para cada uno de los elementos que representan el cuerpo. 3. Unión de las anteriores para obtener la matriz de rigidez general y la obtención del vector de fuerzas general. 4. Resolver la ecuación anterior para las condiciones de contornos dadas obteniendo el vector de desplazamientos, y finalmente 5. Calcular los esfuerzos y las tensiones a partir de los desplazamientos de los nodos obtenidos. Para que el cálculo mediante este método sea lo más preciso posible es muy importante la elección del tipo de elementos en que se va a dividir el cuerpo objeto del estudio. Cuando queremos decidir cómo dividir el cuerpo nos guiamos en gran medida por su geometría, sobre todo por la forma de sus extremos, tanto externos como internos (p.e.: un agujero en el centro de una placa). Así podemos ver como en cuerpos bidimensionales los elementos triangulares o cuadrangulares son una mejor opción que los elementos rectangulares a la hora de discretizar un cuerpo curvo o con una forma extraña. Igualmente en los cuerpos tridimensionales el tetraedro es mejor que el prisma rectangular. Otras veces dependiendo de la forma del cuerpo lo mejor es hacer una división mixta, es decir, con elementos de diferentes formas. 34
42 Fig 20. Ejemplo de malla mixta Método de las Diferencias Finitas El método de las diferencias finitas es un método numérico aproximado basado en sustituir las derivadas parciales que intervienen en la ecuación de Lagrange por diferencias de valores numéricos (en este caso la flecha) en una serie de puntos determinados. Estos puntos se sitúan en una malla que puede ser triangular, rectangular, hexagonal 35
43 Fig 21. Diferentes tipos de malla. En este caso se optará una malla rectangular y regular (separaciones entre los puntos iguales). Mediante este tipo de malla puede estudiarse cualquier forma de contorno de la placa a estudiar. 36
44 Fig 22. Aproximación de un contorno irregular mediante una malla. Para una comprensión más fácil se aplicará los fundamentos del método a una función genérica f(x) en una dirección. La derivada de f(x) en el punto m se define como el límite de Δy/Δx cuando Δx tiende a cero. ( ) ( ) ( ) (3.25) Fig 23. Aproximación de una función mediante diferencias finitas 37
45 Así se puede expresar ( ) ( ) (3.26) La expresión utilizada para Δy la llamamos diferencia hacia delante. También se pueden usar diferencias centrales ( ) (3.27) O diferencias hacia detrás ( ) (3.28) Esta sería la primera derivada, pero se puede hacer de manera similar para la segunda derivada ( ) [ ( )] [ ( )] [( ) ( ) ] (3.29) ( ) Se procede igual para obtener la tercera y la cuarta derivadas. ( ) (3.30) ( ) (3.31) De la primera a la cuarta derivada mediante diferencias centrales quedan de la siguiente manera: (3.32) ( ) 38
46 ( ) (3.33) ( ) (3.34) ( ) (3.35) En este caso, para calcular placas, lo que se quiere es obtener la flecha w que cumpla la ecuación de Lagrange. (3.36) Donde w no depende solo de x, sino de x y de y. Calculando cada uno de los términos mediante diferencias parciales se obtiene: ( ) (3.37) ( ) (3.38) { ( )} {[ ] [ ] [ ] } (3.39) [ ( ) ] Así la ecuación de Lagrange para un punto cualquiera m, n. 39
47 [ ( ) (3.40) ( ) ] En todas las expresiones anteriores λ=δx=δy. Esto es así porque se ha elegido una malla rectangular y regular. Si representamos esta ecuación gráficamente se obtiene lo que se llama molécula computacional. Fig 24. Molécula computacional con los coeficientes de cada elemento. Esta expresión da la flecha en cada punto de la malla que se ha elegido. El conjunto de todas las expresiones de la flecha en cada punto de la malla crea un sistema de ecuaciones. Este sistema de ecuaciones hay que completarlo con las condiciones de contorno del problema. Una vez obtenido w m,n se pueden hallar los esfuerzos (momento flectores y torsores y los cortantes) en cualquier punto m, n. Las expresiones de los esfuerzos serían: 40
48 [ ] (3.41) [( ) ( )] ( ) [ ] (3.42) [( ) ( )] ( ) ( ) { [ ]} (3.43) ( ) { [ ]} (3.44) {[( ) ( ) ] [( ) ( ) ]} [( ) ( )] ( ) { [ ]} (3.45) {[( ) ( ) ] [( ) ( ) ]} [( ) ( )] 41
49 Las condiciones de contorno más usuales, tal y como se ha visto anteriormente, son borde empotrado, borde apoyado y borde libre. Estas condiciones de contorno, hay que expresarlas en diferencias finitas. Al estar usando diferencias centradas, con lo que para poder expresar las condiciones de contorno en diferencias finitas, será necesario introducir unos puntos ficticios fuera de la placa. 1. Borde empotrado Fig 25. Borde empotrado. Expresando las condiciones de contorno de un borde empotrado (la deformación es nula y su derivada es también nula) en diferencias finitas tenemos: ; (3.46) ( ) ( ) (3.47) De esta manera, (3.48) 42
50 2. Borde apoyado En el caso del apoyo simple la deformación es nula y el momento flector es nulo también. (3.49) ( ) (3.50) Y se obtiene que, (3.51) Fig 26. Borde apoyado. 3. Borde libre En el caso de borde libre todos los esfuerzos son nulos. Si el punto central de las diferencias está en el eje libre, hay que introducir cuatro puntos ficticios fuera de la placa. 43
51 Borde Libre Fig 27. Puntos ficticios en borde libre. Al igual que en los casos anteriores las deformaciones de estos puntos ficticios se pueden expresar en función de las deformaciones de puntos de la malla que están dentro de la placa. Para ello especificamos que la fuerza y el momento en el eje son nulos. ( ) ( ) (3.52) ( ) ( ) ( ) (3.53) Para poder eliminar las deformaciones de los puntos ficticios son necesarias dos ecuaciones más, que son las siguientes: 44
52 ( ) ( ) (3.54) ( ) ( ) (3.55) Eliminando los cuatro puntos ficticios de este sistema de ecuaciones y sustituyendo el resultado en la ecuación de Lagrange se obtiene la representación en diferencias finitas de la condición de borde libre. (3.56) Fig 28. Borde libre. Usualmente en el método de las diferencias finitas para que una solución se considere válida se establece el criterio de convergencia. Un método de diferencias finitas es convergente si la solución de la ecuación de diferencias finitas se aproxima a la solución exacta de la ecuación diferencial parcial cuando los tamaños de los pasos en la malla tienden a cero. 45
53 Cuando Δx 0 Para que un método sea convergente ha de ser consistente y estable. La consistencia de un método con la ecuación diferencial parcial se consigue si la ecuación discreta usada por el método es equivalente a la ecuación diferencial cuando el tamaño de paso tiende a cero. La estabilidad se define como: Cuando una ecuación diferencial parcial tiene una solución acotada, se dice que la ecuación de diferencias asociada es estable si produce una solución acotada y es inestable si produce una solución no acotada. El concepto de la estabilidad está ligado con el crecimiento o decrecimiento de los errores que se introducen en la etapa de cómputo. Así se puede decir que un método particular es estable si el efecto acumulativo de todos los errores de redondeo producidos al aplicar un determinado algoritmo es insignificante. La prueba de que una solución aproximada converge a la solución exacta de una ecuación diferencial parcial es generalmente muy difícil, aún en los casos más simples. En el caso de este PFC habría que aplicar esos criterios a una ecuación diferencial parcial de cuarto grado. 46
54 Capítulo 4 Casos y Resultados en Estático 4.1. Resolución analítica. Método de Navier. Para resolver el problema se va a utilizar el método de Navier, Como se ha visto en la introducción este método se usa cuando la placa tiene sus cuatro lados apoyados (como en el caso del problema) y facilita mucho la labor de cálculo. El método de Navier se basa en expresar la función deformación y la función de la carga como un desarrollo en doble serie de Fourier. De esta manera se obtiene que (4.1) (4.2) Donde ( ) (4.3) Y (4.4) 47
55 Siendo D (4.5) Carga Uniforme Fig 29. Placa apoyada en sus cuatro lados sometida a carga uniforme. Los datos de partida del problema son DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Módulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga N/m 2 En este caso, carga repartida y constante el coeficiente P m,n del método Navier es (4.6) 48
56 Como se vio anteriormente en la introducción hay que coger conjuntos completos de términos. Es decir, sumas de m+n iguales. Para este caso se realizarán los cálculos con tres términos, que ofrece suficiente precisión. Termino Valores de m Valores de n Valor de P m,n Valor de W m,n 1 er Término ,139 6,939e-3 2º Término ,370 9,231e ,370 9,231e ,228 8,193e-6 3 er Término ,126 8,193e ,228 9,497e-6 Tabla 2. Valores de P m,n y W m,n para carga uniforme. Ahora hay que sustituir estos valores de W m,n en la expresión de la flecha para el punto central de la placa (x=2, y=2). En el punto central de la placa se da la flecha máxima. Termino Valores de m Valores de n Valor de Valor de la W m,n flecha w (m) 1 er Término 1 1 6,939e-3 6,939e-3 2º Término 1 3 9,231e-5-9,231e ,231e-5-9,231e ,193e-6 8,193e-6 3 er Término 3 3 9,497e-6 9,497e ,193e-6 8,193e-6 Tabla 3. Valores de la flecha para carga uniforme. La flecha en el punto central será la suma de todos los valores de las flechas de la tabla anterior. w(2,2) 6,780e-3 m 49
57 Carga Puntual En este nuevo caso la carga es una carga puntual en el centro de la placa. Se averiguará la flecha en el centro de la placa, que es el punto en el que se dará la máxima deformación. Fig 30. Placa apoyada en sus cuatro lados sometida a una carga puntual. Para ello también se empleará el método de Navier. Los datos de partida del problema son DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Módulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga 1000 N En este caso la expresión de la que se obtiene el coeficiente P m,n es la siguiente (4.7) 50
58 Donde β y μ son las coordenadas del punto donde es aplicada la carga. En este caso coincide con el punto en el que se quiere hallar la flecha (x= β=2; y= μ=2). Y siendo m,n =1,2,3, Valor de Termino Valores de m Valores de n Valor de P m,n W m,n 1 er Término ,068e-3 2º Término ,271e-5 3 er Término ,271e-5 Tabla 4. Valores de P m,n y W m,n para carga puntual. Sustituyendo estos valores en la expresión en doble serie de Fourier para la flecha se obtiene: Termino Valores de m Valores de n Valor de Valor de la W m,n flecha w (m) 1 er Término 1 1 1,068e-3 1,068e-3 2º Término ,271e-5 4,271e-5 3 er Término ,271e-5 4,271e-5 Tabla 5. Valores de la flecha para carga puntual. La flecha es la suma de los valores obtenidos. w(2,2) 1,153e-3 m Carga Triangular El último caso que se va a analizar es el de una placa apoyada en sus cuatro lados y con una carga triangular. 51
59 Fig 31. Placa apoyada en sus cuatro lados sometida a una carga triangular. Los datos de partida del problema son DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 2 m Espesor (h) 0,02 m Material Acero Características de la Placa Módulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga 250*x N/m 2 Para este caso, aplicando el método de Navier, se obtiene que (4.8) Al ser la carga simétrica respecto a la recta y = 1, la deformación de la placa será simétrica respecto a esa misma recta. Sin embargo la carga no es simétrica en x, por lo que se puede esperar que lo sea la deformación. Este hecho dificulta la obtención del punto de máxima deformación, del que si conocemos su coordenada según el eje y, y = 1. Por ello se obtendrán los valores de la deformación a lo largo de la recta y = 1, desde x = 0 hasta x = 4 con incrementos de x de 0,1 m (Δx = 0,1 m). Para obtener de una manera más cómoda estos valores se utilizará una hoja de cálculo de Excel. Para los cálculos al igual que en los dos casos anteriores se emplearán tres términos 52
60 completos de las series de Fourier. Al hacer los cálculos cada 0,1 m se obtiene la flecha en cuarenta puntos de la recta y = 1. Para facilitar el análisis de los resultados obtenidos se realizarán unos gráficos. La tabla con todos los resultados numéricos se adjuntará al final como un anexo. 0,0006 0,0005 0,0004 0,0003 Flecha Total 0,0002 0, Fig 32. Gráfico de la flecha total en la placa sometida a carga triangular. En contra de lo esperado, la solución obtenida de la deformación de la placa es simétrica. Por lo que el punto de máxima deformación según estos resultados sería el x = 2, y = 1. Sin embargo es obvio que la deformación bajo una carga triangular no puede ser simétrica en la dirección del eje x. w(2,1) 5,275E-04 A continuación se incluyen los resultados para todas las combinaciones de valores de m y n de los tres términos calculados. 53
61 0,0006 0,0005 0,0004 0,0003 m=n=1 0,0002 0, Fig 33. Gráfico con los resultados para m = n = 1. 0, , , , , m=1; n=3 m=3; n=1 m=1; n=5 5 m=n=3 m=5; n=1-0, , ,00004 Fig 34. Gráfico para el resto de combinaciones de valores de m y n. Los resultados obtenidos para el primer término son los que mayor peso tiene posteriormente en el valor final obtenido para flecha, por ello se ha preferido representar los resultados obtenidos en dos gráficos. De esta manera se puede 54
62 apreciar más claramente como los resultados, para cualquiera de los valores de m y n, son simétricos respecto a x = 2. Esta incongruencia puede deberse a la propia elección del sistema de cálculo para este modelo. El método de Navier es una aproximación a la realidad y por eso mismo sus resultados son aproximados. Como analogía a este caso se puede ver la deformación de una viga sometida a una carga triangular. Fig 35. Viga sometida a carga triangular. La carga puntual equivalente (p e = ql/2) estaría aplicada a un tercio de la longitud total de la viga (L), tomado desde el extremo con el valor más alto de carga puesto que es la coordenada del centro de gravedad del triángulo. Planteando las ecuaciones de equilibrio en el extremo B obtenemos las reacciones en ambos extremos (4.9) (4.10) La ley de momentos flectores es (4.11) El momento flector máximo 55
63 (4.12) Una vez obtenida la ley de momentos flectores se puede averiguar la expresión de la deformada sin más que integrar la ley de flectores dos veces respecto de x dos veces. (4.13) Al integrar se obtiene [ ] (4.14) Donde C 1 y C 2 son dos constantes de integración que se obtienen al sustituir en f(x) las condiciones de contorno. En este caso la flecha en los dos extremos de la viga es nula. (4.15) (4.16) Y la deformada queda [ ] (4.17) Para hallar la máxima flecha y el punto en que esta aparece hay que derivar la expresión de la deformada e igualar a cero. [ ] (4.18) Se iguala a cero y nos queda [ ] (4.19) 56
64 Es una ecuación de cuarto grado por lo que tendrá cuatro soluciones. Sin embargo solo una de ellas es válida. La flecha máxima se dará aproximadamente en (4.20) Al contrario de lo que ocurría en el caso de la placa al resolver mediante el método de Navier, en este caso la flecha máxima si aparece desplazada del centro de la viga hacia el extremo donde la carga es mayor Resolución mediante el Método de Elementos Finitos (M.E.F.) Ahora se resolverán los mismos casos que antes pero por el Método de Elementos Finitos. Para ello utilizaremos el programa Abaqus. Con este programa se pueden calcular mediante este método desde cargas estructurales a corrientes eléctricas o térmicas. Puesto que posteriormente se incluye como anexo una guía de uso de Abaqus en este apartado nos limitaremos a visualizar los resultados obtenidos una vez analizados los casos por Abaqus, sin explicar paso por paso todo lo que se ha hecho hasta llegar al resultado final. Sí veremos el análisis de sensibilidad de malla realizado en cada uno de los casos, para encontrar el mallado apropiado para cada caso Carga Uniforme Los datos de partida son los siguientes: DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Modulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga N/m 2 57
65 Fig 36. Imagen de Placa Apoyada con Carga Repartida Constante Tal y como se ha comentado en la introducción al M.E.F. el proceso del mallado es fundamental para la posterior validez de los resultados. Para asegurarse de obtener un mallado lo más adecuado posible se realizará un análisis de sensibilidad de malla. En él se realizan diversos mallados sobre el modelo y posteriormente se analizan los resultados. Para este caso escogeremos una malla regular y de elementos rectangulares, es decir, todos los elementos tienen las mismas dimensiones. Empezando por una malla con pocos elementos y partir de ahí se aumentará el número de los mismos hasta encontrar la malla apropiada. En la Tabla 6 se pueden ver los resultados obtenidos con los diferentes tipos de malla. Nº de Elementos Tamaño (m 2 ) Tensión de V. Mises (Pa) Flecha (m) 4 2x2 6,962E+06 6,401E x1 1,026E+07 6,845E ,8x0,8 1,078E+07 6,195E ,5x0,5 1,207E+07 6,784E ,4x0,4 1,245E+07 6,777E ,25x0,25 1,289E+07 6,771E ,2x0,2 1,295E+07 6,770E ,16x0,16 1,293E+07 6,746E-03 Tabla 6. Análisis de sensibilidad de malla para carga uniforme. 58
66 La diferencia de la flecha obtenida entre un mallado con 25 elementos y uno con 64 elementos varía un 9.5%, mientras que la diferencia entre el mallado de 256 elementos y el de 400 elementos es de 0.014%. O entre el mallado de 400 elementos y el de 625 elementos es de 0.35%. En el gráfico de la tensión de V. Mises la diferencia entre el mallado de 256 elementos y el de 400 elementos es de 0.45% y en el caso del mallado de 400 y el de 625 es de 0.15%. Plasmando estos datos en sendos gráficos es posible hacerse, de manera visual, una idea de lo observado en los datos de la tabla. En ambos gráficos los resultados obtenidos con los mallados de pocos elementos tienen una gran variación. Sin embargo cuando se aumenta significativamente el número de elementos los resultados tienden a estacionarse. 1,400E+07 Tensión de V.Mises 1,200E+07 1,000E+07 8,000E+06 6,000E+06 Tensión de V.Mises 4,000E+06 2,000E+06 0,000E Fig 37. Gráfico Tensión de V. Mises para análisis de sensibilidad de malla. 59
67 6,900E-03 6,800E-03 6,700E-03 6,600E-03 6,500E-03 6,400E-03 6,300E-03 6,200E-03 6,100E-03 Flecha Flecha Fig 37. Gráfico de la Flecha para análisis de sensibilidad de malla. Con estos resultados se puede seleccionar una malla óptima con el menor número de elementos posibles. Si bien al ser el modelo un ejemplo sencillo el tiempo que tarda Abaqus en realizar el cálculo es mínimo. El mallado que elegido para este modelo es el de 400 elementos. Podía haberse elegido el mallado de 256 elementos, pero como el tiempo de cálculo no es un problema, se empleará el mallado de 400. Fig 38. Mallado de 400 elementos 60
68 Una vez definido el mallado se puede pasar a visualizar el resultado que nos ofrece Abaqus. El programa nos enseña la deformada de la placa con un mapa de colores en el que el rojo sería la flecha máxima y el azul la mínima. Abaqus también muestra la escala de ese mapa de colores. La flecha máxima se da en el centro de la placa (lo cual era de esperar). Si bien existe un área en que la flecha tiene valores muy similares. Fig 39. Placa deformada En la Fig. 39 se puede ver el resultado obtenido en Abaqus con la forma deformada de la placa. El valor máximo de la flecha obtenido es el siguiente: w max (m) 6,770E Carga Triangular Los datos de partida para este modelo son los siguientes: 61
69 DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 2 m Espesor (h) 0,02 m Material Acero Características de la Placa Modulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga 250*x N/m 2 Fig 40. Placa con carga triangular. Al igual que en el caso de carga repartida, de manera previa al cálculo se realizará un análisis de sensibilidad de malla en este modelo. Nº de Elementos Tamaño (m 2 ) Tensión de V. Mises máx. (Pa) Flecha máx. (m) 2 2x2 7,87E x1 1,604E+06 4,049E ,5X0,5 2,627E+06 5,449E ,57x0,4 2,960E+06 5,251E ,4x0,4 2,721E+06 5,244E ,25x0,25 3,037E+06 5,616E ,2x0,2 3,118E+06 5,617E ,13x0,2 3,146E+06 5,625E ,13x0,13 3,169E+06 5,624E ,1x0,1 3,168E+06 5,664E-04 Tabla 7. Análisis de sensibilidad de malla para carga triangular. 62
70 Se aprecia claramente como los resultados con una malla de pocos elementos no son muy fiables y son muy diferentes a los obtenidos con un número considerable de elementos. Por ejemplo la variación entre la tensión de V. Mises en las mallas de 32 y 35 elementos es del 12,67%. Cogiendo las mallas de 300 y 450 elementos la diferencia sería de 0,7%. En el caso de la flecha, y para las mismas mallas, en el primer caso la variación es del 3,7% y en el segundo del 0,01%. Gráficamente también se puede ver, como a partir de una malla con 200 elementos los resultados no varían tanto y se estabilizan en torno a un valor. 3,50E+06 Tensión V. Mises 3,00E+06 2,50E+06 2,00E+06 1,50E+06 Tensión V. Mises 1,00E+06 5,00E+05 0,00E Fig 41. Gráfico de la Tensión de V. Mises para análisis de sensibilidad de malla. 63
71 0,0006 Flecha máx. 0,0005 0,0004 0,0003 0,0002 Flecha máx. 0, Fig 42. Gráfico de la Flecha para análisis de sensibilidad de malla. La malla que elegida en esta ocasión es la de 450 elementos. Con lo que la flecha máxima obtenida mediante el M.E.F. es: w max (m) 5,624E-04 Esta deformación máxima se obtiene en el punto de coordenadas x=2.53 m; y=1 m. La deformada que muestra Abaqus con ese mallado se puede ver en la Fig
72 Fig 43. Placa deformada bajo carga triangular Carga Puntual Los datos de partida del problema son los siguientes: DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Módulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga 1000 N 65
73 Fig 44. Placa con Carga puntual en dirección negativa al eje Z. Observando la Fig 44. llaman la atención las dos líneas que cruzan la placa dividiéndola en 4 partes iguales. Esto se debe a que para este problema se ha tenido que realizar un paso que no se ha realizado en ninguno de los otros dos casos. Para poder escoger el punto central de la placa para situar la carga hay que realizar una partición de la placa, la cual se lleva a cabo con esas dos líneas que dividen la placa en cuatro partes. Una vez realizado esto Abaqus permite seleccionar el punto central para situar la placa. Hay que resaltar que si bien se divide la placa en cuatro partes los cálculos posteriores no se ven afectados, puesto que luego se utiliza la placa como un continuo. Procediendo igual que en los dos casos anteriores se realiza un análisis de sensibilidad de la malla. Nº de Elementos Tamaño (m 2 ) Tensión de V. Mises máx. (Pa) Flecha máx. (m) 16 1x1 1,26E-03 2,11E ,66x0,66 1,23E-03 2,65E ,5x0,5 1,22E-03 3,07E ,4x0,4 1,21E-03 3,40E ,25x0,25 1,21E-03 4,12E ,2x0,2 1,21E-03 4,46E ,15x0,15 1,21E-03 4,87E ,13x0,13 1,21E-03 5,09E+06 Tabla 8. Análisis de sensibilidad de malla para carga puntual. 66
74 Fig 45. Gráfico de la Flecha para análisis de sensibilidad de malla. Fig 46. Gráfico de la Tensión de V. Mises para análisis de sensibilidad de malla. En los gráficos (Fig 45 y 46) y en la Tabla 8 se puede ver como para la flecha máxima a partir de 400 elementos el valor de la flecha se estabiliza, mientras que en el caso de la tensión de Von Mises ese valor no se estabiliza, sin embargo si suaviza en gran medida su crecimiento respecto a los primeros valores. Al ser la versión de Abaqus empleada una versión de estudiante, no 67
75 admite más de 1000 elementos. De esta manera parece claro que la malla a elegir en este caso es la de 900 elementos. Fig 47. Deformada Carga Puntual En la figura de la deformada se visualizan también el punto máximo de deformación y el mínimo. El punto de máximo deformación como es obvio está en el centro de la placa. El valor máximo de la placa obtenido por Abaqus es: w max (m) 1,209E Resolución mediante diferencias finitas Para la resolución de los diferentes casos según el método de las Diferencias Finitas, se usará el programa Matlab, el cual permite crear de manera más o menos sencilla pequeños programas ejecutables con los que se pueden visualizar los resultados. 68
76 En este capítulo se mostraran los resultados obtenidos una vez ejecutados los programas, así como se explicará el desarrollo de los mismos. La forma tradicional de resolver el cálculo de placas mediante diferencias finitas implica la resolución de un sistema de ecuaciones. Cuantos más puntos tiene la malla más se complica la resolución del sistema, puesto que se hace más complicada la programación de la matriz de coeficientes. Además este sistema implica la inversión de la matriz de coeficientes, que puede ser un proceso arduo y muy complicado dependiendo de la matriz en sí. Por estas razones la aproximación que se realizará a este problema será desde una perspectiva diferente. Como se ha visto antes la ecuación general de la deformación de una placa desarrollada en diferencias finitas es: [ ( ) (4.21) ( ) ] Despejando el término w m,n se obtiene: [ ( ) (4.22) ( ) ] A esta ecuación hay que añadir las condiciones de contorno de la placa, cuyas ecuaciones se presentaron en el capítulo de la teoría de diferencias finitas. Una vez obtenido todo lo anterior, se realiza un programa que haga un barrido por todos los puntos de la malla aplicando esa ecuación. Para obtener un resultado válido no es suficiente con un solo barrido de los puntos de la malla, hay que hacer una serie de iteraciones. Una vez realizado el barrido el suficiente número de iteraciones para conseguir una solución aceptable se obtendrá una matriz solución, en la que estarán representados los valores de las flechas en todos los nodos de la malla. 69
77 Para averiguar el valor de cada uno de los puntos de la malla, se utiliza la llamada molécula computacional. Debido al uso de esta molécula computacional hay puntos de la malla que están situados fuera de la placa. Son los llamados puntos imaginarios. Pero son fundamentales para poder obtener el valor de la flecha en cada uno de los nodos de la malla. Para realizar el cálculo de la flecha, el programa partirá de una matriz de resultados inicial de valor nulo, todos los elementos serán iguales a cero. Posteriormente se inicia realmente el programa. Para poder calcular de manera más sencilla la flecha en la malla se ha dividido esta en varias zonas diferentes: Zona Central (1): Es la zona central de la placa, donde se aplica la carga. Zonas Izquierda (2) y Derecha (3): Son las zonas a izquierda y derecha de la zona central de la malla en las que hay puntos imaginarios. Zonas Inferior (4) y Superior (5): Son las zonas inferior y superior a la zona central de la malla en las que hay puntos imaginarios. Fig 48. Zonas de la malla para el cálculo de la placa. Se puede observar en la imagen como hay algunos nodos que no están asignados a ninguna zona. En estos nodos a la deformación se le asigna valor nulo. 70
78 En general las zonas serán las mismas en todos los casos, aunque habrá variaciones en alguno de ellos que se expondrán dentro del caso correspondiente. Por comodidad en los cálculos se emplearan unas mallas lo más similar posible a las utilizadas en el M.E.F Carga Uniforme El primer caso estudiado es el de una placa cuadrada apoyada en todo su contorno con carga uniforme distribuida en toda a superficie. Los datos de partida son los siguientes: DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Módulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga N/m 2 Al igual que en el caso de los elementos finitos se realizaba un análisis de sensibilidad de malla para obtener la malla más apropiada para los cálculos, ahora realizaremos un análisis similar para obtener el número de iteraciones necesario para obtener una solución aceptable. De esta manera se puede ver la convergencia de la solución hacia un valor determinado. 71
79 Fig 49. Gráfica de iteraciones vs. deformación para carga uniforme. En la gráfica se puede observar como al principio el valor de deformación máxima en el centro de la placa cambia mucho según se van aumentando las iteraciones. A partir de iteraciones el valor se hace constante y ya no crece más. Ese valor de la deformación será la solución de nuestro estudio. Se escoge el valor correspondiente a iteraciones iteraciones w max (m) 6,758*10-3 Matlab también permite crear una gráfica en la que se ve una superficie deformada con los valores de deformación que se dan en la placa. 72
80 Fig 50. Deformación de la placa bajo carga uniforme. Realmente este análisis solo confirma la convergencia de la solución en un punto de la placa, en este caso el centro. Para ver si la solución es aceptable en todos los puntos hay que realizar otra comprobación adicional. Comprobación de la solución obtenida Realizar una comprobación punto por punto de la malla es una tarea inviable, pues se perdería demasiado tiempo. Por ello esta comprobación ha de comparar todos y cada uno de los datos de la malla de una sola vez. Esta comprobación final consiste en comparar los resultados obtenidos en la última iteración con los resultados de la penúltima iteración de manera que se obtiene una matriz en la que cada uno de sus elementos es el porcentaje de variación de ese elemento en la matriz deformación entre las dos últimas iteraciones. [ ] [ ] [ ] [ ] (4.23) Esta ecuación solo se aplica en la zona de la malla correspondiente a al interior de la placa. Lo que interesa es que los valores de los elementos de la matriz [c] sean lo más pequeños posible. Esta matriz tiene cerca de 400 elementos, lo 73
81 que hace que revisar uno por uno todos los elementos de la matriz sea una tarea ardua. Sin embargo, representando gráficamente la nueva matriz [c] se puede ver rápidamente si la variación es muy grande y en qué punto se da la mayor variación. Fig 51. Comprobación de la solución para carga uniforme. En la gráfica se observa los valores de variación de la última iteración con respecto a la inmediatamente anterior, los cuales son prácticamente nulos (valores de orden ). Y el valor máximo de esta variación en la matriz [c] es 7,6540* La solución para iteraciones también es válida en el resto de puntos de la placa. 74
82 Carga Triangular Los datos de partida para este modelo son los siguientes: DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 2 m Espesor (h) 0,02 m Material Acero Características de la Placa Módulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga 250*x N/m 2 Procediendo igual que en el caso anterior, se realiza un análisis de los resultados en el punto de mayor deformación de la placa respecto al número de iteraciones. El espesor de malla usado en los cálculos es λ=0,1 m. Fig 52. Gráfica de iteraciones vs. deformación para carga triangular. Se observa en la gráfica como a partir de las iteraciones el valor de la flecha comienza a estabilizarse. Se elige de nuevo el valor correspondiente a iteraciones. 75
83 iteraciones w max (m) 5,664*10-4 Si se representa la superficie deformada: Fig 53 Placa deformada bajo carga triangular. El punto de máxima deformación en este caso no será el centro de la placa, si no que estará desplazado hacia el lado donde el valor de la carga sea mayor. El punto donde se da la máxima flecha es x=2.5 m; y =1 m. en la siguiente gráfica se visualiza mejor el desplazamiento del punto de máxima deformación hacia el extremo de la placa donde se da la mayor carga. 76
84 Fig 54. Deformada de la placa bajo carga triangular. Comprobación de la solución obtenida Al igual que en el caso anterior ahora se realizará la comprobación de la solución respecto a la solución para la iteración anterior. Para ello se aplicará la ecuación (4.23). Fig 55. Comprobación de la solución para carga triangular. 77
85 En este caso la máxima diferencia entre los valores de una iteración y la anterior es de 1,0424*10-6. Que es un valor muy bajo. Por lo que la solución obtenida para iteraciones es válida no solo para el punto de máxima deformación, sino también para el resto de la placa Carga Puntual Los datos de partida del problema son los siguientes: DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Módulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga 1000 N Al igual que anteriormente para los cálculos se empleará una malla similar a la utilizada en el cálculo mediante el M.E.F. aunque al ser un caso especial por ser la carga puntual, también se empleará un malla más fina. De este modo se puede establecer una comparación entre las mallas. Además con el programa Matlab no existe una limitación del número de nodos, mientras que con el M.E.F. si hay una limitación del número de elementos a emplear en el estudio. Es importante tener en cuenta que la carga en la ecuación general de las placas expresada en diferencias finitas es solo válida si esta es uniforme. En este caso es una carga aplicada en un punto. Para poder aplicar la ecuación (4.22) se aproximará la carga puntual a una carga uniforme distribuida en una pequeña superficie. Esta pequeña superficie será lo más pequeña que se pueda en función de la malla empleada y ha de ser lo suficientemente grande como para albergar una molécula computacional completa. Matlab, a diferencia de lo que ocurre en Abaqus, no admite un valor de espesor de malla que dé un número no entero de nodos. Por ello, para este caso en particular se realizará más de un estudio de malla. En Abaqus se había escogido un espesor de malla de 0.13 m. Para resolverlo por diferencias finitas se estudiarán primero unas mallas con 0,2 m, 0,1 m y 0,05 m de espesor de malla respectivamente. 78
86 Malla con espesor de malla λ=0,2m Cumpliendo con la condición expresada anteriormente sobre el tamaño de la superficie en la que repartir la carga equivalente, para esta malla se obtiene una superficie de 0,8x0,8 m 2. Realizando el mismo análisis de iteraciones que en los casos anteriores, se obtiene: Fig 56. Gráfico de iteraciones vs deformación para carga puntual. Se observa en el gráfico como a partir de iteraciones el valor máximo de la deformación se estabiliza en 1,73*10-3 m. 79
87 Fig 57. Deformada bajo carga puntual. El valor de la deformación máxima en la placa bajo una carga puntual obtenido mediante diferencias finitas es: iteraciones w max (m) 1,730*10-3 Comprobación de la solución obtenida Al comprobar la solución se obtiene que la variación de la solución entre las últimas iteraciones es del orden de La máxima variación es 8,7398*
88 Fig 58. Imagen de la matriz de comprobación para espesor de malla 0,2 m. Malla con espesor de malla λ=0,1m En el análisis de sensibilidad de malla se obtiene el siguiente resultado: Puntual (λ=0,1) 2,00E-03 1,80E-03 1,60E-03 1,40E-03 1,20E-03 1,00E-03 8,00E-04 6,00E-04 4,00E-04 2,00E-04 0,00E Deformación Fig 59. Gráfico de iteraciones vs deformación para carga puntual. En el gráfico se observa como el valor de la deformación se estabiliza a partir de iteraciones. Si bien se estabiliza en un valor similar al resultado obtenido con la malla anterior, han sido necesarias más iteraciones para lograrlo iteraciones w max (m) 1,836*
89 Fig 60. Deformada bajo carga puntual El caso actual es un caso especial puesto que toda la carga está aplicada en un solo punto. Con una malla muy gruesa, podrían darse discontinuidades, debido a efectos locales producidos por la carga. Por esa razón es normal emplear una malla más fina de manera que se disminuyen esas discontinuidades y puede obtenerse un resultado más fiable. A continuación se empleara una malla de 0,05 m de espesor. Comprobación de la solución obtenida Fig 61. Imagen de la matriz de comprobación para espesor de malla 0,1 m. 82
90 Deformación (m) En el caso de la malla con espesor de 0,1 m, los valores de variación obtenidos son del orden de 10-7 y el máximo es 8,9833*10-7. Malla con espesor de malla λ=0,05m Para esta nueva malla también hay que realizar el análisis de convergencia con el número de iteraciones. Puntual (λ=0,05) 2,00E-03 1,80E-03 1,60E-03 1,40E-03 1,20E-03 1,00E-03 8,00E-04 6,00E-04 4,00E-04 2,00E-04 0,00E Iteraciones Series1 Fig 62. Gráfico de iteraciones vs deformación para carga puntual. Al disminuir el tamaño de la malla aumenta el número de iteraciones a partir de las cuales el resultado se estabiliza. En este caso se estabiliza a partir de 5 millones de iteraciones. Obviamente este aumento de las iteraciones conlleva un aumento del tiempo de cálculo iteraciones w max (m) 1,865*
91 Fig 63. Deformada bajo carga puntual Procediendo igual que en los dos casos anteriores se realiza una comprobación de la solución. Comprobación de la solución obtenida Fig 64. Imagen de la matriz de comprobación para espesor de malla 0,05 m. En este caso los valores de variación de la solución son del orden de Y el valor máximo en toda la malla es 3,7303*10-7. La solución sería válida en todos los puntos de la malla. 84
92 Capítulo 5 Introducción a la Dinámica de Placas 5.1. Introducción a la dinámica En este capítulo se estudia el comportamiento dinámico de una placa bajo cargas transversales a la misma. Fig 65. Placa sometida a una fuerza transversal F(x,y,t). 85
93 Se realizaran las operaciones tomando una placa de longitud a, ancho b y espesor h sometida una carga distribuida en su superficie F(x,y,t). El material de la placa tiene densidad ρ y módulo de elasticidad E. El objetivo de este estudio es averiguar los movimientos de todos los pintos de la placa en función del tiempo, es decir, averiguar la función w(x,y,t). Los lados de la placa estarán sujetos a determinadas condiciones de contorno. Estas normalmente serán alguna de las que se citan a contunuación para el caso de estudio dinámico. 1. Borde apoyado Esta condición nos dice que tanto el desplazamiento en cualquier dirección como el momento flector en el borde apoyado son nulos. Si el borde apoyado es x = 0, (5.1) (5.2) 2. Borde empotrado En este caso está impedido cualquier movimiento, por lo que el desplazamiento y el giro son nulos. Si el borde empotrado es x=0 (5.3) (5.4) 3. Borde libre 86
94 Ahora los desplazamientos no son nulos. Sin embargo todos los esfuerzos (momentos y cortantes) en el borde libre si son nulos. Siendo x = a el borde libre. (5.5) (5.6) (5.7) Si hacemos igual que en el caso estático entonces tenemos que las reacciones de borde han de ser nulas. Recordemos que estas reacciones de borde se obtienen de la suma del cortante y de un cortante adicional obtenido de representar el momento torsor como un par de fuerzas. [ ] [ ( )] (5.8) [ ] [ ( )] (5.9) Además en los estudios dinámicos se parte de unas condiciones iniciales de desplazamiento y velocidad en el instante inicial (t=0) que hay que tener en cuenta a la hora de resolver el problema. Generalmente estas condiciones iniciales son desplazamiento y velocidad nulos en t=0. (5.10) (5.11) Hay que reseñar que en este desarrollo y tal y como se viene haciendo habitualmente en este tipo de cálculos, se despreciará el efecto del peso propio de la estructura ya que es menor que el resto de fuerzas que actúan sobre la estructura. O lo que es lo mismo, los desplazamientos que origina el peso son despreciables frente a los que originan las fuerzas externas. 87
95 Si al igual que en el casó estático se considera un diferencial de placa de tamaño dx x dy se obtiene la ecuación diferencial de la placa. Solo que en este caso también depende del tiempo. ( ) (5.12) La solución a la ecuación () será la suma de la solución de la ecuación homogénea (F(x,y,t) = 0) más una solución particular de la ecuación completa. La ecuación homogénea se puede asemejar a un problema de vibraciones libres, mientras que la particular se asemeja a un problema de vibraciones forzadas. (5.13) Donde w h (x,y,t) es la solución homogénea y w p (x,y,t) la particular. La ecuación homogénea se obtiene cuando las fuerzas externas a la placa son nulas (F(x,y,t)=0). La ecuación diferencial de la placa queda de la siguiente manera ( ) (5.14) Un método muy habitual para resolver la ecuación homogénea es la separación de variables, o método de Galerkin, que consiste en obtener una solución como suma de unas funciones espaciales y temporales. De esta manera (5.15) Donde ϕ i y q i son las funciones espaciales y temporales respectivamente. Introduciendo la ecuación (5.15) en la ecuación (5.14) se obtiene 88
96 ( ) ( ) (5.16) ( ) Esto debe cumplirse para cualquier sumando i-ésimo. Así se obtiene (5.17) Así se obtienen dos ecuaciones dos ecuaciones diferenciales, una para la componente espacial ϕ i (x,y) y otra para la componente temporal q i (x,y). Para poder resolverlas hay que tener en cuenta las condiciones iniciales y las condiciones de contorno. (5.18) (5.19) La componente espacial ϕ i (x,y) es el modo propio de la vibración libre. Y se puede escribir: (5.20) 89
97 Si se asume que el movimiento es armónico: (5.21) Esto se puede expresar de esta otra forma: (5.22) Siendo x 0 e y 0 constantes arbitrarias, ω la frecuencia propia y k n y c m son el número de ondas. Insertando esta expresión el la ecuación general del movimiento [ ] (5.23) Reordenando: (5.24) Los desplazamientos w(x) y w(y) se pueden escribir: (5.25) (5.26) Las expresiones exponenciales pueden expresarse como una sucesión de senos y cosenos: (5.27) (5.28) Donde las constantes A i y B i se obtienen aplicando las condiciones de contorno las expresiones anteriores e imponiendo la condición de que exista al menos 90
98 una solución diferente a la trivial, quedando una de las constantes con un valor libre. Las funciones ϕ n,m (x,y) son los llamados modos propios de vibración y ω la frecuencia propia de vibración. Estos parámetros son fundamentales para saber el movimiento transversal de la placa. Es importante tener en cuenta la ortogonalidad de los modos propios la cual se rige por las siguientes ecuaciones: (5.29) Es habitual normalizar los modos propios fijando un valor de la constante libre, cumpliéndose la siguiente expresión: (5.30) En el caso de placa apoyada en sus cuatro lados las frecuencias y los modos propios serían: (5.31) ( ) (5.32) Una vez obtenidas las frecuencias y los modos propios se analizará la ecuación del movimiento completa. Para ello se usará la transformada de Laplace y el método de Galerkin. 91
99 Se aplica la transformada de Laplace a la ecuación de movimiento de manera que tenemos W(x,y,s) y F(x,y,s) en vez de w(x,y,t) y F(x,y,s) respectivamente. ( ) (5.33) Ahora se aplica el método de Galerkin que como se ha visto antes desarrolla una función como una serie de productos de funciones. En este caso: (5.34) Donde las funciones ϕ j (x,y) son los modos propios obtenidos en la parte de vibraciones libres. Estos modos propios solo deben cumplir las condiciones de contorno y geométricas. De esta manera se puede expresar la ecuación general de la siguiente forma: (5.35) ( ( ) ) Para simplificar la expresión se aplica la ortogonalidad de los modos propios. ( ) ( ) (5.36) ( ) ( ) (5.37) Así se obtiene la siguiente ecuación matricial: (5.38) 92
100 [ ] { } { } Donde: [ ] [ ] [ ] { } { } { } { } (5.39) (5.40) (5.41) (5.42) Donde M y K son matrices cuyos elementos pueden calcularse de la siguiente manera: (5.43) La matriz M resulta ser la matriz identidad. ( ) (5.44) Y la matriz K es también una matriz diagonal. Para obtener las funciones q k (t), hay que aplicar la inversa de Laplace: { } { } [[ ] { }] (5.45) 93
101 Capítulo 6 Casos y Resultados en Dinámico 6.1 Resolución Analítica. Al igual que en el caso estático en este apartado/capítulo se resolverán unos ejemplos de placas pero de manera dinámica. El objetivo es al igual que en el caso estático hallar la flecha, aunque en este caso la deformación dependerá del tiempo. Anteriormente se ha visto que la flecha puede expresarse como un producto de unas funciones espaciales por unas temporales. (6.1) Donde las funciones espaciales son los modos propios de vibración. Estos modos propios solo dependen de las características de la placa y de las condiciones de contorno. Para una placa apoyada en los cuatro lados los modos propios de vibración se expresan de la siguiente manera: (6.2) La frecuencia propia: ( ) (6.3) 94
102 Las funciones temporales se obtienen aplicando la transformada de Laplace. La ecuación matricial mediante la cual obtenemos las funciones temporales es: { } { } [[ ] { }] (6.4) Carga Uniforme Fig 66. Placa apoyada en sus cuatro lados sometida a carga uniforme. Los datos de partida del problema son DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Densidad (ρ) 7850 kg/m 3 Módulo de Young (E) 210*10-9 Coeficiente de Poisson (ν) 0,3 Carga Carga N/m 2 Condiciones iniciales Desplazamiento 0 m Velocidad 0 m/s Del mismo modo que en el caso estático se averiguará la flecha en el centro de la placa(x=2, y=2), que es el punto en que se dará la mayor deformación. 95
103 Generalmente para estos cálculos se suelen coger los cuatro primeros términos con m y n iguales. De manera que los modos propios y las frecuencias propias serían Valor de Termino Valores de m Valores de n Valor de ϕ m,n w m,n 1 er Término º Término er Término o Término Tabla 9. Valor de los modos propios y frecuencias propias para carga uniforme con 4 términos de m y n. Para poder hallar las funciones temporales primero hay que averiguar las matrices M y K que nos dan la matriz P. [ ] [ ] [ ] (6.5) Las matrices M y K son matrices diagonales por lo que la matriz P también será diagonal. [ ] (6.6) [ ] Ahora hay que hallar el vector f. { } { } (6.7) { } ( ) 96
104 El hecho de que las matrices sean diagonales facilita mucho los cálculos, pues se pueden realizar como ecuaciones individuales directamente. (6.8) { } ( ) Una vez obtenido el vector se aplica la inversa de Laplace a dicho vector y obtenemos las funciones q(t). Multiplicando estas funciones por los modos propios de vibración hallados anteriormente se obtiene la flecha en función del tiempo. (6.9) Para una mejor comprensión de los resultados obtenidos se obtendrán mediante Matlab unas gráficas en las que se pueden visualizar los valores de los términos calculados, así como una gráfica a de la deformación obtenida al sumar la deformación de todos los términos. El segundo y cuarto términos son nulos y no aportan nada al valor final de la deformación. 97
105 Fig 67. Deformación en función del tiempo para el primer término de Galerkin. Fig 68. Deformación en función del tiempo para el tercer término de Galerkin. 98
106 Fig 69. Deformación en función del tiempo total. Hasta ahora se han usado cuatro términos de Galerkin (N=4). La deformación total se obtendrá al sumar las deformaciones obtenidas para cada uno de los términos. Deformación por términos Deformación total w 1 =0,01385 m w 2 = 0 m w t =0,01386 m w 3 =1,898 E-05 m w 4 = 0 m Tabla 10. Valores de deformación para cada uno de los términos y deformación total para carga uniforme. JUSTIFICACIÓN DE N=1 PARA LOS CÁLCULOS DE LA DEFORMACIÓN A continuación se va a analizar que aporta cada elemento al resultado global y si se puede pasar de N=4 a N=1 y de esta manera reducir el tiempo de los cálculos. Para ver el peso que tiene cada desplazamiento en el total de la deformación se aplicará la siguiente ecuación: 99
107 (6.10) Para la deformación obtenida con el primer término se obtiene: Teniendo en cuenta este resultado se puede llegar a la conclusión de que no es necesario realizar los cálculos con cuatro términos de Galerkin si no que se pueden realizar usando N=1 e igualmente tener una alta precisión en el resultado. Por esta razón a partir de ahora en los próximos casos solo se empleará N= Carga Puntual En este nuevo caso la carga es una carga puntual en el centro de la placa. Se averiguará la flecha en el centro de la placa, que es el punto en el que se dará la máxima deformación. Fig 70. Placa apoyada en sus cuatro lados sometida a una carga puntual. Los datos de partida del problema son 100
108 DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Densidad (ρ) 7850 kg/m 3 Módulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga 1000 N Condiciones iniciales Desplazamiento 0 m Velocidad 0 m/s En este caso, y al ser la misma placa que en el caso anterior (sometida a diferente carga), los modos propios y las frecuencias son las mismas. Termino Valores de m Valores de n Valor de ϕ m,n Valor de w m,n 1 er Término Tabla 11. Frecuencia y valor propio del primer término para carga puntual. Las matrices M y K también serán las mismas que en el caso anterior. [ ] (6.11) [ ] A la hora de averiguar f(s) en este caso de carga puntual aparece un problema. Este sistema sólo es válido para cargas repartidas, no puntuales. Para poder sortear esta limitación se creará una carga repartida equivalente a la carga puntual del problema, localizada en un área pequeña alrededor del punto central de la placa. El área elegida será cuadrada y con una superficie de 0,04x0,04 m 2. Con estos datos la carga resultante equivalente es de N/m 2. { } (6.12) De esta manera: 101
109 { } ( ) (6.13) Multiplicando una vez más por el resultado obtenido para las funciones espaciales y aplicando la inversa de Laplace se obtiene: Fig 71. Deformación en función del tiempo para carga puntual. De esta gráfica se puede sacar el valor de la máxima deformación. w max (m) 2,135E Carga Triangular El último caso que se va a analizar es el de una placa apoyada en sus cuatro lados y con una carga triangular. 102
110 Fig 72. Placa apoyada en sus cuatro lados sometida a una carga triangular. Los datos de partida del problema son DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 2 m Espesor (h) 0,02 m Material Acero Características de la Placa Densidad (ρ) 7850 kg/m 3 Módulo de Young (E) 210*10 9 Coeficiente de Poisson (ν) 0,3 Carga Carga 250*x N/m 2 Condiciones iniciales Desplazamiento 0 m Velocidad 0 m/s Al igual que en los casos anteriores, hay que averiguar los modos propios y las frecuencias propias, las matrices K y M y el vector f. Se ha averiguado la flecha en el punto de máxima deformación obtenido en el caso estático resuelto mediante el MEF, x=2,533 m, y=1 m. Termino Valores de m Valores de n Valor de ϕ m,n Valor de w m,n 1 er Término 1 1 0, ,5479 Tabla 12. Frecuencia y valor propios para carga triangular. 103
111 [ ] [ ] { } (6.14) Por lo tanto el vector de las funciones temporales q será: { } ( ) (6.15) Una vez multiplicado este vector de funciones temporales por los ϕ m,n y realizada la inversa de Laplace, se puede realizar la gráfica para este caso. Fig 73. Deformación en función del tiempo para carga triangular. En la gráfica se pueden ver representados el máximo valor de la deformación (línea roja) y la media (línea verde). La deformación máxima obtenida es: w max (m) 1,013E
112 COMENTARIOS GENERALES A TODOS LOS CASOS Las gráficas se han realizado mediante la obtención de 300 puntos con el programa Matlab. En estos cálculos no se ha tenido en cuenta el amortiguamiento interno de la placa, por ello en ningún momento los datos tienden a estacionarse en un valor. 6.2 Resolución Mediante M.E.F. ara la resolución mediante elementos finitos de la parte dinámica, al igual que en la parte estática se usará el programa Abaqus. A la hora de seleccionar la malla a utilizar en los diferentes casos se ha decidido mantener la malla empleada en los casos estáticos. Al ser un análisis dinámico, la solución no es un solo valor de la flecha, sino que este valor depende del tiempo. Tal y como se verá en una gráfica más adelante, si se hacen los cálculos con un intervalo de tiempo muy grande entrará en juego el amortiguamiento interno de la propia placa, y los valores obtenidos tenderán a estabilizarse en un valor. Ese valor coincide con el hallado en los casos estáticos. Inicialmente se empleará un valor para el intervalo de tiempo de 3 segundos Carga Uniforme Los datos de partida son los siguientes: DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Densidad (ρ) 7850 kg/m 3 Módulo de Young (E) 210*10-9 Coeficiente de Poisson (ν) 0,3 Carga Carga N/m 2 Condiciones iniciales Desplazamiento 0 m Velocidad 0 m/s 105
113 Fig 74. Imagen de Placa Apoyada con Carga Repartida Constante Al igual que se realizó anteriormente en el problema estático, ahora se averiguará la flecha en el centro de la placa Al realizar el análisis en Abaqus, este nos devuelve una gráfica, en ella se puede ver la evolución de la deformación en función del tiempo del análisis. También se puede obtener y representar el valor máximo de esa deformación. 106
114 Fig 75. Análisis dinámico carga uniforme 3 segundos. En este caso de carga uniforme y con un análisis de 3 segundos el valor máximo de la flecha en el centro de la placa es: w max (m) 1,134*
115 Carga Triangular Los datos de partida para este modelo son los siguientes: DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 2 m Espesor (h) 0,02 m Material Acero Características de la Placa Densidad (ρ) 7850 kg/m 3 Módulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga 250*x N/m 2 Condiciones iniciales Desplazamiento 0 m Velocidad 0 m/s Fig 76. Placa con carga triangular. En este caso no se buscará la deformación en el centro de la placa, pues no es el punto en el que se va a dar la máxima flecha. Se averiguará la deformación en el punto en que se dio el desplazamiento mayor en el caso estático. Es decir, el punto de coordenadas x=2.53 m; y=1 m. El primer análisis ofrece los siguientes resultados: 108
116 Fig 77. Análisis dinámico carga triangular 3 segundos. El valor máximo de la flecha resulta: w max (m) 8,697*
117 Carga Puntual Los datos de partida del problema son los siguientes: DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Densidad (ρ) 7850 kg/m 3 Módulo de Young (E) 210E9 Coeficiente de Poisson (ν) 0,3 Carga Carga 1000 N Condiciones iniciales Desplazamiento 0 m Velocidad 0 m/s Fig 78. Placa con Carga puntual en dirección negativa al eje Z. Al igual que en los casos anteriores inicialmente se realiza un análisis durante 3 segundos, cuyo resultado se plasma en la siguiente gráfica. En ella podemos ver la oscilación de la deformación durante ese tiempo y el máximo valor obtenido. w max (m) 2,332E
118 Fig 79. Análisis dinámico carga puntual 3 segundos. 6.3 Resolución mediante Diferencias Finitas Al igual que en la parte estática, para la solución dinámica mediante este método se usará el programa informático Matlab. Para desarrollar esta solución se partirá del desarrollo realizado en la parte estática. Aunque obviamente habrá que introducirle algunos cambios pues en este caso los movimientos de la placa también depende del tiempo. Recordando la ecuación general de las placas en dinámica: ( ) (6.16) En los casos planteados para su resolución la carga aplicada a la placa es constante en el tiempo: 111
119 ( ) (6.17) En esta ecuación se pueden ver los términos espaciales y temporales. La ecuación estática para la deformación de una placa expresada en diferencias fintas se ha visto anteriormente, así que solo queda expresar el término temporal en diferencias finitas y añadirlo. ( ) (6.18) Así pues la ecuación general de la dinámica de placas quedaría de la siguiente manera: [ ( ) (6.19) ( ) ] ( ) Para realizar los cálculos hay que aplicar las condiciones de contorno, que serán las mismas que en los casos en estático, y las condiciones iniciales. CONDICIONES DE CONTORNO Las condiciones de contorno son las mismas que para el caso estático. Como la placa está poyada en todos sus extremos la flecha y el momento flector en los apoyos son nulos. Expresados en diferencias finitas: (6.20) ( ) (6.21) 112
120 Así: (6.22) Estas condiciones se cumplen para cualquier instante de tiempo. CONDICIONES INICIALES Tal y como se ha comentado anteriormente las condiciones iniciales son que el desplazamiento en la plaza y la velocidad son nulos para t=0. (6.23) Hay que expresar estas condiciones mediante diferencias finitas. La primera de las condiciones es muy sencilla de expresar en diferencias finitas. (6.24) La segunda condición se expresa de la siguiente manera (6.25) Despejando: (6.26) 113
121 Aplicando esta dos condiciones en la ecuación general de las placas en dinámica expresada en diferencias finitas se obtiene que en el instante inicial, o lo que es lo mismo cuando t=0: (6.27) A partir de estas tres ecuaciones ya se puede iniciar la realización del programa en Matlab, pues se conoce el valor de la flecha en tres instantes de tiempo consecutivos, el instante inicial (t=0), el anterior (t=-δt) y el posterior (t= Δt). Matlab no admite instantes negativos y tampoco no enteros. El valor de Δt puede tomar valores negativos (como en el primer instante) y valores no enteros. Por ello se introduce el índice k. A continuación se incluye a modo de aclaración una tabla de relación instantes tiempo. Tiempo Instante t=-δt k=1 t=0 k=2 t= Δt k=3 t=2* Δt k=4 T=3* Δt k=5 Tabla 13. Relación tiempo instante. En la ecuación de la dinámica de placas en diferencias finitas están presentes los valores de la flecha en tres momentos de tiempo. Así pues se puede despejar la flecha de un instante específico en función de la flecha en los otros dos. Para los casos que se van a resolver se ha optado por despejar el instante k+1 en función de los dos instantes anteriores k y k-1. Al despejar el instante k+1 y convertirlo en el instante actual (k) que se quiere averiguar los otros dos instantes se convierten en k-1 y k-2 respectivamente. La ecuación general de las placas quedará de la siguiente manera: 114
122 [ [ ( ) ( ) ]] (6.28) Anteriormente se había calculado el instante inicial (t=0), el anterior y el posterior. Introduciendo el instante inicial (k=2) y el posterior (k=3) en la ecuación anterior, se puede averiguar el instante siguiente (k=4). Una vez obtenido k=4 se introduce en la ecuación de nuevo junto con k=3, de esta manera se averigua k=5. Y así sucesivamente Método Dinámico Fig 80. Esquema desarrollo diferencias finitas dinámica Al igual que en la parte estática, en la parte dinámica también hay separar las diferentes zonas de la placa. Estas zonas serán las mismas que en los casos resueltos en estático. Un factor importante a tener en cuenta a la hora de hacer los cálculos es el valor de Δt. Este valor indica el espesor de la malla temporal. Es como el espesor λ de la malla de la placa, pero aplicado en el tiempo. Cuanto menor sea Δt, más fina será la malla temporal. Pero hay que tener en cuenta que 115
123 cuanto más pequeño se haga Δt habrá que emplear más instantes para poder encontrar una solución relevante. Por ejemplo, si Δt=0.01 s, el índice k debe ser k=300. De esta manera el tiempo de estudio será de tres segundos, el mismo que se ha empleado en los estudios analítico y mediante MEF. Para empezar los cálculos se empleará Δt=0,001 s. Este valor no está escogido al azar. Se ha empleado los resultados obtenidos analíticamente. A partir de ellos se averigua el periodo (T) de la vibración obtenida como resultado y se divide entre 100. Este es el valor que se utilizará para iniciar los cálculos. (6.29) Una vez obtenidos los primeros resultados a partir de este valor de Δt, se variará para intentar obtener una mejora en los resultados. El espesor de la malla a utilizar será el empleado en los cálculos en estática Carga Uniforme Los datos de partida son: DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Densidad (ρ) 7850 kg/m 3 Módulo de Young (E) 210*10 9 Coeficiente de Poisson (ν) 0,3 Carga Carga N/m 2 Condiciones iniciales Desplazamiento 0 m Velocidad 0 m/s Los cálculos se iniciarán con un valor de espesor de malla λ=0,2 m y Δt=0,001 s. Posteriormente se irá reduciendo el valor de Δt. 116
124 Una vez ejecutado el programa hecho en Matlab se obtiene una gráfica en la que se visualiza la deformación del punto central de la placa en función de los instantes. El número de instantes durante los cuales correrá el programa se cogerán de manera que una vez multiplicado por Δt el resultado sea de 3 segundos. Para este caso el número de instantes será También se puede obtener una representación gráfica de la placa deformada en cada uno de los instantes que calcula el programa. Fig 81. Deformación punto medio de la placa con Δt=0,001 s. Se observa como la gráfica no muestra un resultado válido pues se interrumpe a la altura del instante número 200 y no va más allá. Incluso el programa no muestra la gráfica de la deformación de la placa. Esto se puede deber la elección de Δt no es correcta y este valor debe ser menor. Tal y como se ha dicho anteriormente se procede a reducir el valor de Δt. Como segunda aproximación se emplearán los valores Δt=0,0001 s y Δt =0,00001 s. Si posteriormente fuese necesario, se emplearán valores de Δt inferiores a estos. La reducción del Δt conlleva un aumento del número de instantes, de esta manera se mantiene el tiempo del estudio en 3 segundos, a y instantes respectivamente. El objetivo es conseguir que un valor de la deformación en el centro de la placa que sea igual o muy parecido aunque cambie el valor de Δt. 117
125 Fig 82. Deformación punto medio de la placa con Δt=0,0001 s y Δt=0,00001 s. Con la reducción de Δt se consiguen unos resultados más acorde con la realidad. En las gráficas se observa una vibración en la cual la deformación máxima es: w max (m) Δt (s) Nº Instantes 0, ,564*10-3 0, ,565*10-3 0, Tabla 14. Flecha máxima obtenida para diferentes Δt, para carga uniforme. 118
126 El valor de la deformación es prácticamente idéntico para ambos valores de Δt. Por lo que no es necesario reducir más Δt Carga Puntual DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 4 m Espesor (h) 0,02 m Material Acero Características de la Placa Densidad (ρ) 7850 kg/m 3 Módulo de Young (E) 210*10 9 Coeficiente de Poisson (ν) 0,3 Carga Carga N Condiciones iniciales Desplazamiento 0 m Velocidad 0 m/s Tal y como se ha hecho anteriormente se inician los cálculos con Δt=0,001 s y posteriormente se reducirá. El espesor de malla en este caso es λ=0,05 m. Los resultados obtenidos se resumen en la siguiente tabla: w max (m) Δt (s) Nº Instantes 0, , ,727*10-3 0, ,728*10-3 0, Tabla 15. Flecha máxima obtenida para diferentes Δt, para carga puntual. Al igual que pasó en el caso de carga uniforme, con los dos primeros valores de Δt, el programa no da resultados válidos. Sin embargo con los dos siguientes valores de Δt si se obtienen resultados decentes. Los valores obtenidos de flecha máxima son prácticamente idénticos en ambos casos. 119
127 Fig 83. Deformación punto medio de la placa con Δt=0,00005 s y Δt=0,00001 s. En las gráficas obtenidas como resultado se observa una variación del valor de la deformación en función del instante en que es calculado. Una vibración. 120
128 Carga Triangular DATOS DEL PROBLEMA Longitud (a) 4 m Geometría de la Placa Ancho (b) 2 m Espesor (h) 0,02 m Material Acero Características de la Placa Densidad (ρ) 7850 kg/m 3 Módulo de Young (E) 210*10 9 Coeficiente de Poisson (ν) 0,3 Carga Carga 250*x N/m 2 Condiciones iniciales Desplazamiento 0 m Velocidad 0 m/s El valor del espesor de malla que se usará para realizar los cálculos será λ=0,1 m. w max (m) Δt (s) Nº Instantes 0, , ,283*10-3 0, ,283*10-3 0, Tabla 16. Flecha máxima obtenida para diferentes Δt, para carga triangular. Tanto el resultado numérico como la gráfica son idénticos para los dos valores de Δt que ofrecen resultados significativos. 121
129 Fig 84. Deformación en el punto de máxima flecha de la placa con Δt=0,00005 s y Δt=0,00001 s Método Dinámico Iterativo Este otro método está basado en el anterior, con la inclusión de una parte de iteraciones en cada instante. Es decir, en cada instante se realizarán una serie de iteraciones con el objetivo de intentar conseguir un resultado mejor. El valor de la deformación en los tres primeros instantes se sigue calculando aparte, puesto que es constante a lo largo de la placa. En el resto de instantes inicialmente se aplica la ecuación de la dinámica de placas utilizada en el 122
130 método anterior y posteriormente se aplica la ecuación de las placas en estático a la hora de realizar iteraciones. Cuando se hacen las iteraciones en un instante determinado se emplea la ecuación de la estática, puesto que dentro de un instante concreto el tiempo es constante. Una vez que se han hecho todas las iteraciones y se ha obtenido un valor de la deformación se pasa al siguiente instante. Fig 85. Esquema desarrollo método combinado dinámico estático. Por ejemplo, suponiendo que se está averiguando la deformación en el instante k=4; esta deformación está en función de los dos instantes anteriores (k=2 y k=3). Se aplica la ecuación de las placas en dinámica y se averigua un valor de la deformación en k=4. Este valor se introduce en la ecuación de la estática y se itera. Una vez realizadas todas las iteraciones se ha obtenido un nuevo valor de la deformación para el instante k=4. Este nuevo valor es el que se usará para averiguar la deformación en los instantes siguientes. A la hora de programa este método hay que tener cuidado con el número de iteraciones que se van a emplear, porque entran en juego los factores limitantes de tiempo de cálculo y capacidad del ordenador. Si disminuye mucho Δt, aumentará el número de instantes a calcular manteniendo la ventana de tres segundos para el estudio. Si en cada instante se emplean muchas iteraciones entonces el programa puede tardar demasiado tiempo en encontrar 123
Simulación numérica con elementos finitos de un micro-relé
Simulación numérica con elementos finitos de un micro-relé Titulación: Alumno: Directores: Ingeniería Industrial César Alfaro Sánchez Juan Álvaro Fuentes Moreno y Francisco Periago Esparza Cartagena, 7
I.- CONCEPTOS BÁSICOS DEL MÉTODO DE LOS ELEMENTOS FINITOS
I.- CONCEPTOS BÁSICOS DEL MÉTODO DE LOS ELEMENTOS FINITOS 1.1.- Introducción Muchos problemas de importancia práctica que frecuentemente aparecen en ingeniería, resultan de una complejidad matemática tal
A continuación se enumeran las etapas del proyecto y las principales actividades que han tenido lugar dentro de ellas.
DISEÑO Y ANÁLISIS MEDIANTE ELEMENTOS FINITOS DE ESTRUCTURAS DENTALES IMPLANTO-SOPORTADAS. Autor: Frutos Martín, Jesús de. Director: Fraiz Cosano, David. Entidad colaboradora: PIC Dental. RESUMEN DEL PROYECTO
3º Tema.- Síntesis de mecanismos.
Universidad de Huelva ESCUELA POLITECNICA SUPERIOR Departamento de Ingeniería Minera, Mecánica y Energética Asignatura: Ingeniería de Máquinas [570004027] 5º curso de Ingenieros Industriales 3º Tema.-
Máster Universitario en Ingeniería de Caminos, Canales y Puertos Introducción al Análisis Numérico
Máster Universitario en Ingeniería de Caminos, Canales y Puertos Introducción al Análisis Numérico Departamento de Matemática Aplicada Universidad Granada Introducción El Cálculo o Análisis Numérico es
Realidad virtual en arquitectura
T E C N O L O G Í A Realidad virtual en arquitectura ALEJANDRO EZEQUIEL JUÁREZ sánchez Imagen que muestra una de las salas equipadas para la investigación en el diseño arquitectónico. Los elementos estructurales
Grupo de Ingeniería Gráfica y Simulación Escuela Técnica Superior de Ingenieros Industriales Universidad Politécnica de Madrid ALTIMETRÍA DE OBRAS
Grupo de Ingeniería Gráfica y Simulación Escuela Técnica Superior de Ingenieros Industriales Universidad Politécnica de Madrid ALTIMETRÍA DE OBRAS ALTIMETRÍA DE OBRAS 1. INTRODUCCIÓN: PLANTA, TRAZA Y RASANTE...
Software gratuito como herramienta docente para cálculo de estructuras. Enrique Relea Gangas 1, Andrés Martínez Rodríguez 2
Software gratuito como herramienta docente para cálculo de estructuras Enrique Relea Gangas 1, Andrés Martínez Rodríguez 2 1 Campus de La Yutera, Avda. Madrid 44, 34004 Palencia 1 Tel 979 10 83 13 erelea@iaf.uva.es
Hallar gráfica y analíticamente la resultante de los siguientes desplazamientos: hacia el Noroeste), B. (35 m Sur)
VECTORES: OPERACIONES BÁSICAS Hallar gráfica y analíticamente la resultante de los siguientes desplazamientos: hacia el Noroeste), B (0 m Este 30º Norte) y C (35 m Sur) Solución: I.T.I. 94, I.T.T. 05 A
VIGAS Y MARCOS PLANOS ANALIZADOS POR EL MÉTODO DE ELEMENTOS FINITOS JULIAN DAVID TORO DUQUE
VIGAS Y MARCOS PLANOS ANALIZADOS POR EL MÉTODO DE ELEMENTOS FINITOS JULIAN DAVID TORO DUQUE UNIVERSIDAD TECNOLÓGICA DE PEREIRA FACULTAD DE INGENIERÍA MECÁNICA PEREIRA 007 VIGAS Y MARCOS PLANOS ANALIZADOS
8.- MÉTODOS GENERALES: ANÁLISIS MATRICIAL
8.- MÉTODOS GENERALES: ANÁLISIS MATRICIAL 1 8.1 FLEXIBILIDAD Y RIGIDEZ 8.1.1 Concepto de flexibilidad.- a) La ley de Hooke aplicada a una barra de longitud L y sección A que, sometida a un esfuerzo axil
PROBLEMAS DE CINEMÁTICA DE MECANISMOS
TEORÍA DE MÁQUINAS PROBLEMAS DE CINEMÁTICA DE MECANISMOS Antonio Javier Nieto Quijorna Área de Ingeniería Mecánica E.T.S. Ingenieros Industriales Capítulo 1 GRADOS DE LIBERTAD. 1.1. PROBLEMA. En la figura
UNIVERSIDAD DE IBAGUE FACULTAD DE INGENIERIA PROGRAMA DE INGENIERIA MECANICA ASIGNATURA: MECÁNICA ANALÍTICA CÓDIGO: 2120 AREA: MECANICA CICLO: BÁSICA SEMESTRE: III PRERREQUISITOS: ALGEBRA LINEAL, FÍSICA
PATOLOGÍAS DE ORIGEN TÉRMICO EN ESTRUCTURAS
PATOLOGÍAS DE ORIGEN TÉRMICO EN ESTRUCTURAS Ing. Eduardo Pedoja Profesor de Hormigón Armado y Proyecto Facultad de Ingeniería, Universidad de Montevideo Una de las causas más frecuentes de la aparición
INTRODUCCION A LA FISICA NEWTONIANA Manuscrito de cátedra
INTRODUCCION A LA FISICA NEWTONIANA Manuscrito de cátedra ADVERTENCIA: manuscrito en estado de preparación muy preliminar, particularmente en lo que respecta a la secuencia temática, orden y terminación

References: Resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución 
 resolución 
 resolución 
 resolución 
 Resolución 
 Resolución 
 Resolución 
 resolución 
 resolución 
 resolución 
 Resolución 
 Resolución 
 resolución 
 Resolución 
 resolución