INFORME FINAL “ANÁLISIS ESTRUCTURAL POR EL METODO DE ELEMENTOS FINITOS DE UNA EMBARCACIÓN FLUVIAL PETROLERA DE 20 TONELADAS” JEFE DE PROYECTO : Ing. JORGE VERA ERMITAÑO Alumno : ANTHONY BRITO TICLAVILCA TENORIO Alumno : FELIX ANTONIO MANRIQUE LOPEZ 2014 INDICE 1 ANALISIS ESTRUCTURAL POR EL METODO DE ELEMENTOS FINITOS DE UNA EMBARCACION PETROLERA DE 20 TONELADAS Pág. INTRODUCCION ……………………………………………… 3 1. Marco teórico ……………………………………………… 4 2. Modelación Estructural 3. Simulación Numérica ……………………………………………… 16 …………………………………………… 17 4. Generación y validación del Software Ansys ………………………. 16 5. Carga Rodante ………………………………………………. 26 6. Modelación Estructural ……………………………………………….. 27 7. Simulación Numérica ……………………………………………….. 29 8. Generacion de malla ………………………………………………… 31 9. Condiciones de Contorno ……………………………………………….. 33 10. Resultados ………………………………………………………………… 35 11. Conclusiones y Bibliografía ……………………………………………… 36 2 ANALISIS ESTRUCTURAL POR EL METODO DE ELEMENTOS FINITOS DE UNA EMBARCACION PETROLERA DE 20 TONELADAS INTRODUCCION El proyecto que a continuación se muestra consiste en el estudio, mediante el método de elementos finitos, de la resistencia estructural y el comportamiento del casco de embarcaciones fluviales, sin esta información técnica es probable que produzcan desperdicios económicos, tanto en la inversión inicial, al usar miembros estructurales sobredimensionados, debidos a que los costos de construcción se basan en el peso del acero, así también en los gastos de operación de la embarcación. Calcular la estructura de un barco, presenta grandes dificultades, debido a la gran complejidad tanto de la estructura como del tipo de cargas que se presentan, es por esto que podemos aplicar criterios probabilísticos basados en experiencia, como lo son los métodos de las sociedades clasificadoras, permitiéndonos evaluar la integridad estructural de una embarcación con un alto margen de seguridad. Para llevar a cabo el estudio se utilizará el software de elementos finitos ANSYS, haciendo previamente una revisión completa del método de elementos finitos y de las teorías pertinentes de resistencia de materiales, tanto para el caso de vigas como planchas. El objetivo general de este trabajo es analizar los niveles de esfuerzos y deformaciones producidos en una embarcación fluvial en condiciones características de carga. Con este trabajo se lograra un adecuado adiestramiento y una buena experiencia en el desarrollo de modelos estructurales de embarcaciones 3 fluviales utilizando el método de elementos finitos y de esta manera, en un futuro se podrán modelar estructuras de embarcaciones más complejas. 1) Marco teórico El Método de los Elementos Finitos (MEF), es un método numérico para la resolución de sistemas de ecuaciones en derivadas parciales. Su desarrollo desde los años cincuenta hasta la actualidad ha sido constante y actualmente puede considerarse como el método de análisis numérico más extendido en la mayoría de los ámbitos de la ingeniería. Un programa de cálculo por Elementos Finitos consta de tres módulos bien diferenciados: el Pre-procesador, el Módulo de análisis y el Post-procesador. Fig.1. Discretizacion del Modelo En el Pre-procesador se hace el modelado, define propiedades del material, se define tipo de elemento, características del elemento y la discretización (malla de elementos finitos) realizada para resolver el problema. En el Módulo del análisis (Solver), en este primero se elige el tipo de análisis a realizar, luego se 4 aplican las cargas y condiciones de contorno del problema a analizar para finalmente resolverlo. En el Post-procesador se pueden visualizar los resultados del problema, se pueden presentar los desplazamientos de los nodos, los mapas de tensiones y deformaciones, la configuración deformada del sólido, gráficas de evolución de cualquiera de las variables analizadas, etc. 1.1) Análisis Estructural del Buque A) Tensiones y Deflexiones Primarias La estructura global del buque es denominada Viga – Buque (ESTRUCTURA PRIMARIA). Tensiones y deflexiones primarias son aquellas tensiones y deflexiones obtenidas por el momento resultante de la viga – buque. La estructura del fondo es normalmente más robusta, consiguiendo que el eje neutro se encuentre normalmente en la parte inferior de la sección transversal. Por lo tanto la condición de máxima tensión de Arrufo y Quebranto es normalmente determinado en la cubierta principal. Los miembros que determinan la resistencia estructural longitudinal o primaria son: Cubiertas, doble fondo, fondo, longitudinales de fondo, longitudinales de cubierta y quillas. B) Tensiones y Deflexiones Secundarias Ciertas partes del buque pueden ser consideradas apoyadas en la estructura restante, debido a que la parte adyacente presenta mayor rigidez que la parte considerada. Un ejemplo será el hecho de que los mamparos y costados en la dirección longitudinal presenta una rigidez mucho mayor que la del panel reforzado, el cual está delimitado por el mamparo y el costado. A este panel reforzado se le llama ESTRUCTURA SECUNDARIA. Las tensiones y deflexiones de este panel reforzado, se les llama tensiones y deflexiones Secundarias. C) Tensiones y Deflexiones Terciarias La tensión y deflexión en una unidad del enchapado es sometida a un cargamento lateral son denominados TENSIONES Y DEFLEXIONES TERCIARIAS. 5 Las cargas laterales son soportadas a través de una distribución lineal de tensiones de tracción y compresión en su espesor, pero esta se da para pequeñas deflexiones. Para deflexiones mayores aparecen las llamadas TENSIONES DE MEMBRANA. Para la condición de Empotrado Se considera deflexión pequeña cuando: 𝑊 ≤ 0.5 𝑡 Para la condición de libre apoyo Se considera deflexión pequeña cuando: 𝑊 ≤1 𝑡 Dónde: W: Deflexión local , t : Espesor 1.2) Panel Reforzado Son estructuras que están formadas por una placa y vigas que actúan como refuerzos, generalmente en direcciones ortogonales. El objetivo de las vigas es soportar el cargamento transversal que actúa sobre el enchapado, manteniendo en su posición y geometría de forma que soporte las cargas existentes en el plano de estructura y transmitir estas cargas hacia la estructura del buque. Los paneles reforzados, más allá de su estanqueidad, confiere alta rigidez en el sentido longitudinal y transversal de la viga – buque. Estos paneles deben de ser diseñados para soportar las tensiones primarias, secundarias y terciarias que son el producto de la flexión de la viga – buque, cargas locales y cargas laterales en la placa respectivamente. Placa Colaborante: Es el nombre que adopta la placa cuando está trabaja como brida de los refuerzos longitudinales y transversales aumentando de esta manera la rigidez de módulo de sección. Entonces como consecuencia de la aplicación de esta placa: Ancho efectivo < Distancia entre refuerzos 6 Concepto de Shear Lag y Ancho Efectivo: La teoría de Vigas establece que las tensiones de deflexión son DP a la distancia del punto considerado hasta el eje neutro, y está se debe a que sus secciones permanecen planas. En muchos modelos físicos la flexión no es causada por flexión pura y si por cargas transversales que son absorbidas por el alma de la viga. Bajo efecto de estas cargas el alma es curvada induciendo deformaciones máximas en las bridas. Además estas deformaciones se originan en el alma y solo llegan a la brida por causa del cizallamiento. Fig.2. Efecto “Shear Lag” La flexión y el corte causan contracciones y distorsiones en el plano de la platabanda (brida), esta distorsión del corte provoca que los extremos de la platabanda se deformen más que la parte central, debido a esto resulta que en la parte central hay menores tensiones a comparación a sus extremos, en consecuencia se hace menos EFECTIVA. Fig.3. Efecto “Shear Lag” en vigas abiertas y tipo casco Las platabandas al deformarse en su plano longitudinal, provocan la distorsión de la sección transversal (Alabeo), este efecto no es contemplado por la TEORIA DE VIGA SIMPLE. 7 Debido a la influencia del corte, las tensiones de flexión lejos del alma avanzan sobre esta última (SHEAR LAG). Mientras que el ancho del enchapado, que cuando se calcula el Momento de Inercia de la sección, dará la tensión correcta en la intersección alma – platabanda, para la teoría de la viga simple (ANCHO EFECTIVO). Llegamos a la conclusión de que la mejor manera de tener en cuenta el efecto de corte es a través del ANCHO EFECTIVO. Fig.4. Longitud Efectiva de perfiles en vigas abierta y tipo casco 1.3) Método De La Teoría De Vigas 1. Estimar los niveles de esfuerzo de un modo “rápido” y “aceptable”. Es necesario: Calcular el cargamento: Cargas concentradas y/o cargas distribuidas. Escoger condición de apoyo que mejor representa el comportamiento de la viga en función de las cargas aplicadas. Seleccionar escantillones. Seleccionar el punto de análisis. Calcular el momento flector máximo producido por la carga aplicada. Determinar la inercia de la sección transversal, incluyendo el ancho efectivo de la plancha colaborante que sea solidaria a la viga. Estimación del esfuerzo producido por la flexión pura de la viga utilizando la conocida expresión: σFLEXION = MxY t 8 Cuando es aplicada la teoría de vigas para estimar los esfuerzos en paneles reforzados los siguientes efectos son despreciados: Efectos cortantes Rigidez torsional del panel Efecto de Poisson Unión con los otros elementos 1.4) Teoría De La Plancha Ortotropica En este método el panel reforzando es representado por una Plancha Ortotropica equivalente cuya rigidez a la flexión y torsión varia en dos direcciones ortogonales. El grado de aproximación depende fuertemente del número de refuerzos en cada dirección y de la uniformidad del mismo. Si el espesor de la plancha ortotrópica rectangular es relativamente pequeño en comparación a las otras dimensiones, (a, b/t>40), el esfuerzo normal y el efecto de deformación por corte pueden ser despreciados Las relaciones constitutivas para una plancha ortotrópica son: 𝜎𝑋 𝐸𝑋 1 { 𝜎𝑌 } = [𝐸𝑌 𝜈𝑋 1 − 𝜈𝑋 𝜈𝑌 𝜏𝑋𝑌 0 𝐸𝑋 𝜈𝑌 𝐸𝑌 0 𝜀𝑋 0 0 ] { 𝜀𝑌 } 𝐺(1 − 𝜈𝑋 𝜈𝑌 ) 𝛾𝑋𝑌 Schade (1941) introdujo curvas de diseño basado en planchas ortotropicas. Asumió que la deflexión es pequeña en comparación con la espesura de la plancha W(max) ( t ≤ 0.5 ) y utilizó cuatro tipos diferentes de paneles reforzados: Tipo A: Fuertemente reforzado en ambas direcciones Tipo B: Refuerzos repetitivos en una sola dirección y solo uno en la otra dirección. Tipo C: Refuerzos en una sola dirección. Tipo D: Plancha sin refuerzos. 9 Para cada uno de los tipos anteriores se pueden admitir cuatro casos de condiciones de frontera: Todos los lados simplemente apoyados. Los lados menores empotrados y los mayores simplemente apoyados. Los lados mayores empotrados y los menores simplemente apoyados. Todos los lados empotrados. Parámetros iniciales a ser conocidos: Presión lateral, geometría del panel (largo, ancho, separación entre refuerzos largos y cortos, Inercias, etc.) Para cada uno de los tipos de paneles se debe calcular varios parámetros con el fin de utilizar las gráficas propuestas por el Dr. Shade. Por ejemplo para un panel “Tipo A”, las expresiones para cada parámetro son: Rigidez unitaria dado que representan el momento de inercia por largo unitario 𝑖𝑎 = 𝑖𝑏 = 𝐼𝑛𝑎 (𝐼𝑎 − 𝐼𝑛𝑎 ) +2 𝑆𝑎 𝑏 𝐼𝑛𝑏 (𝐼𝑏 − 𝐼𝑛𝑏 ) +2 𝑆𝑏 𝑎 “razón de aspecto virtual”, puesto que representa la razón de aspecto del panel alterado por la razón de rigidez unitaria en cada una de las direcciones. 𝜌= 𝑎 4 𝑖𝑏 √ 𝑏 𝑖𝑎 “coeficiente de torsión” pues representa la relación entre la inercia del material sujeto a corte y la inercia del material sujeto a flexión. 𝐼𝑝𝑎 𝐼𝑝𝑏 𝜂= √ 𝐼𝑛𝑎 𝐼𝑛𝑏 Dónde: a = Longitud del panel b = Ancho del panel (a> b) 10 𝐼𝑛𝑎 y 𝐼𝑛𝑏 = Momento de inercia incluyendo ancho efectivo de la placa de los refuerzos repetitivos longitudinales y transversales respectivamente. Sa(Sb) = respectivamente. Separación entre refuerzos longitudinales y transversales ia y ib = Rigidez unitaria en dirección longitudinal (transversal) 𝛾 = Coeficiente de Poisson t = Espesor de placa 𝐼𝑝𝑎 y 𝐼𝑝𝑏 = Momento de Inercia de la placa 𝐼𝑎 𝑦 𝐼𝑏 = Momento de inercia incluyendo ancho efectivo de la placa de refuerzos longitudinales repetitivos (transversales) Por otro lado Schade (1941) indica posibles variaciones del método de la plancha ortotrópica para lugares con concentración de rigidez. Fig.5. Curvas de Shade en paneles reforzados En este siglo Galindo (2007) realizó varios experimentos para identificar el grado de exactitud del método de la plancha ortotrópica para posibles concentraciones de rigidez. 11 Realizó modelos con elementos finitos de diversos paneles reforzados, variando la razón de aspecto y la intersección o no de reforzadores en el centro del panel. Fig.6. Factores de deflexión 1.5) Teoría de la Viga Buque 12 A) Obtención de la Cargas 13 14 B) Obtención de los Esfuerzos primarios 15 2) Modelación Estructural 3D Modelos estructurales son considerados como modelos de la plantilla de cualquier estructura del edificio. Se ofrece una visualización exacta de la estructura del edificio propuesto y modelado en 3D es únicamente el proceso de creación de una representación matemática, armazón de alambre de cualquier objeto en tres dimensiones utilizando un software especializado. Se utiliza para crear el diseño en 3D precisa y rápida los modelos. 3D de modelado estructural es una parte de los servicios de ingeniería estructural y puede ser identificado como el proceso de elaboración de una representación matemática, armazón de alambre de cualquier objeto 3D utilizando técnicas de computación especializados. Se utiliza en el desarrollo de la construcción de modelos excelentes, personalizados modelos estructurales de viviendas y el modelado de estructuras para distintos tipos de edificios, tales como edificios residenciales, comerciales e industriales. Fig.7. Modelo Estructural 3D La creación de un equilibrado modelo estructural 3D es una colección de los pasos detallados, a partir de una Sketch-Up de objeto. Modelación estructural tiene que ver con la creación de modelos de construcción 3D precisos y fiables según los criterios requeridos. Modelado en 3D de estructuras se utiliza para crear modelos exactos de construcción estructural de elementos constructivos y estructurales diferentes modelos en 3D de diferentes estructuras de los edificios juegan un papel decisivo para los constructores, los contratistas y arquitectos para construir sus estructuras de los edificios propuestos con mayor eficacia. 16 De modelado estructural ayuda a encontrar volúmenes específicos del componente y que permite a los constructores o contratistas que entender acerca de la cantidad de materiales y componentes serán necesarios. También es útil para la toma de vistas en 3D exteriores, vistas en 3D de interiores con las especificaciones de dimensiones, los modelos de paisaje de diseño estructural y los modelos de ecuaciones estructurales. El modelo estructural para la barcaza fluvial será el modelo de tipo superficie, ya que los espesores calculados según la norma se dan al momento de la simulación. Fig.8. Modelo Isométrico 3) Simulación Numérica El modelo estructural final es exportado al software ANSYS con formato IGES, para poder comenzar con la simulación numérica. Los métodos numéricos en Ingeniería proporcionan una herramienta que permite traducir modelos matemáticos en procedimientos computacionales, cuyos resultados pueden ser contrastados con las soluciones analíticas, en los casos en los que éstas existan. En este trabajo se muestra una forma de introducir al estudiante de ingeniería en el estudio de las ecuaciones diferenciales en derivadas parciales y su modelado numérico. En particular la ecuación de difusión será el objeto de estudio. La elección de la misma se efectúa basándose en sus múltiples aplicaciones en problemas de la Mecánica. Usando el método de diferencias finitas para determinar la solución numérica de la ecuación de difusión, se abordarán básicamente dos 17 planteos: método explícito y método implícito. Para ambos procedimientos es necesario diseñar algoritmos computacionales. Este trabajo pretende mostrar cómo contribuir a la formación matemática del ingeniero mediante el uso de métodos numéricos. Cómo éstos introducen a los alumnos en la programación para la solución de Ecuaciones Diferenciales Parciales. Finalmente muestra la comparación de los resultados obtenidos con el uso de distintos métodos numéricos y la solución analítica hallada en los casos posibles. Estos son los 7 pasos para tener la solución numérica completa, de acuerdo a las condiciones dadas en el problema. Aquí están incluidas las 3 etapas para toda simulación numérica. - Pre-procesamiento: El primer paso en FEA, pre-procesamiento, es construir un modelo de elementos finitos de la estructura a ser analizada. En muchos paquetes de FEA se requiere de la entrada de una descripción topológica de las características geométricas de la estructura. Ésta puede ser 1D, 2D, o 3D. El objetivo principal del modelo es replicar de manera realista los parámetros importantes y características del modelo real. La manera más sencilla para conseguir similaridad en el análisis es utilizar planos pre existente, modelos CAD, o datos importados de un ambiente FEA. Una vez se ha creado la geometría, se utiliza un procedimiento para definir y dividir el modelo en "pequeños" elementos. En general, un modelo de elementos finitos está definido por una malla, la cual está conformada por elementos y nodos. Los nodos representan puntos en los cuales se calcula el desplazamiento (análisis estructural). Los paquetes de FEA enumeran los nodos como una herramienta de identificación. Los elementos están determinados por conjuntos de nodos, y definen propiedades localizadas de masa y rigidez. Los elementos también están definidos por la numeración de la malla, la cual permite referenciar la correspondiente 18 deflexión o esfuerzo (en análisis estructural) para una localización específica. - Módulo de análisis: En la siguiente etapa en el proceso de análisis de elementos finitos se lleva a cabo una serie de procesos computacionales que involucran fuerzas aplicadas, y las propiedades de los elementos de donde producir un modelo de solución. Tal análisis estructural permite la determinación de efectos como lo son las deformaciones, estiramiento o estrés que son causados por fuerzas estructurales aplicadas como lo son la fuerza, la presión y la gravedad. - Post-procesamiento: Estos resultados entonces pueden ser estudiados utilizando herramientas visuales dentro del ambiente de FEA para ver y para identificar completamente las implicaciones del análisis. Herramientas numéricas y gráficas permiten la localización precisa de información como esfuerzos y deformaciones a ser identificadas. 4) Generación de Malla Una malla es un conjunto de celdas contiguas que permite representar en forma discreta el dominio de un problema a resolver numéricamente. La confiabilidad de la solución numérica obtenida depende de la calidad de la malla generada. Los criterios de calidad a usar dependen del método numérico escogido y del tipo de problema a resolver. Fig.9. Enmallado de cuadernas En esta presentación se discutirá sobre los métodos para modelar y describir la geometría del dominio, los tipos de mallas más usados, algoritmos para generar mallas, mostrando con más detalle uno para 19 generar mallas mixtas, criterios de refinamiento y de calidad usados, el estado de la investigación actual en el área. Las mallas geométricas se pueden clasificar en estructuradas y no estructuradas. Las mallas estructuradas se caracterizan estar compuestas de celdas de un tamaño similar y del mismo tipo, por ejemplo, triángulos o rectángulos, en dos dimensiones, tetraedros o hexaedros en tres dimensiones. Estas mallas son fáciles de generar pero no permiten modelar eficientemente problemas geométricamente complejos. Por otro lado, las mallas no estructuradas permiten el uso de celdas de distinto tipo y/o de diferente tamaño. Estas mallas requieren de algoritmos y estructuras de datos más complejos que las anteriores, pero permiten modelar geometrías complejas y optimizar el número de celdas usadas de acuerdo a la necesidad de la aplicación. Actualmente, la generación de una malla adecuada en dos dimensiones se considera un problema resuelto, pero no así en tres dimensiones. En tres dimensiones se trabaja en el diseño algoritmos eficientes, tanto secuenciales como paralelos, en todas las etapas del proceso de generación de una malla de calidad, para poder resolver cada vez problemas más complejos. Fig.10. Algoritmos de malla 5) Manual Ansys 20 Modelar una viga bi-apoyada de longitu l = 3,8 + 0,2*N (m), sujeta a una carga distribuida de q = 100*(0,9 + 0,1*N) (kN/m), utilizando 20 elementos en la viga. A sección recta da viga es formato I con alas iguales de espesor de 2 cm e ancho 0,9 + 0,1*N (m) y alma de espesor de 2 cm y alto 0,45 + 0,05*N (m). El material es acero común. Seleccione y presente en una plantilla MS-Excel (tabla y gráfico) los siguientes parámetros de respuesta: 1. Momentos flectores y esfuerzos cortantes verticales en todos los nodos del modelo (a partir de los valores seleccionados para los elementos). 2. Esfuerzo en todos los nodos del modelo en los puntos de la sección. Solucionando: I. Abierto el worckbench ,Ir a static de Analys systems. II. Seleccionar unidades ir a milímetros y ok. 21 III. Aquí vas a graficar la sección de la viga I con una longitud L. Click em ZX 22 Ir a modeling plano Zx Sketc1 (ver figura) 23 IV. Aquí vas a graficar la sección de la viga I. 24 Dimensionar la viga Click em line body 25 Hasta aquí se acabó la geometría de la viga, guardar en una carpeta y salir del DM, y regresamos a la primera ventana. V. Hacer 2 click em Model , aqui vamos a general el mallado de la viga. 26 27 Tamaño del elemento L/20 Des pues generar la malla de la siguiente forma_ 28 Malla de la viga con 20 elementos, guardar y salir de esta ventana. VI. Aquí vamos a colocar las condiciones de contorno y las cargas 2 Click en setup 29 1 click aquí,es para seleccionar los puntos de los extremos de la viga y colocar la condicion de contorno 30 Colocar em x=y=z=0 Rotation X=free Rotation Y e Z = 0 APLICAR CARGA LINEAL 31 DATOS PARA LA RESPUESTA, DESLOCAMENTO, FUERZAS Y MOMENTO 32 DEFORMATIO/TOTAL TOOLS/BEAM SOLVE (RESUELVE EL PROBLEMA) 33 34 Adicionar bending moment e shear force la cual me daran las fuerzas cortantes y momentos em cada nodo,luego mandar SOLVE de nuevo. Para conocer el momento flector en cada nodo hacer lo siguiente. AL EXPORTAR EL VA GENERAR UM ARCHIVO TEXTO Y UM EXCEL DE LA FORMA SIGUIENTE(DE FORMA AUTOMATICA EL ABRE EXCEL SI NO ES ASI EL ESTA EN LA CARPETA DONDE GUARDAS EL ARCHIVO) 35 LO MISMO HACEMOS PARA LAS FUERZAS Y ESFUERSOS. 6) Calculo Estructural Estructuras fundamentales de la barcaza Los esfuerzos más importantes del grupo esfuerzos estructurales, son los esfuerzos longitudinales. Estos esfuerzos son derivados del comportamiento del buque como viga, esto es, de lo que se llama viga buque. Los elementos de la viga buque deben ser dimensionados de modo que resistan los esfuerzos a que están sometidos con un adecuado margen de seguridad. Este dimensionamiento se realizara en cada caso de acuerdo con la teoría de elasticidad y la resistencia de materiales, pero con distintas hipótesis simplificativas, ya que la viga buque se presenta mucho mas complicada que las vigas de la teoría. La viga buque está constituida esencialmente por el forro externo cerrado superiormente por una cubierta. Esta viga hueca, si bien apta para resistir esfuerzos longitudinales, está altamente deformable en sentido transversal, 36 dado que su espesor es muy pequeño con respecto a sus dimensiones. Así que tal viga hueca debe ser soportada o reforzada por elementos transversales y longitudinales que forman un verdadero esqueleto sobre el que se fija el forro. Este esqueleto está formado por una serie de marcas transversales unidas entre sí por elementos longitudinales. Estos marcos transversales están constituidos por cuadernas al costado, baos en la cubierta y varengas en el fondo. Los elementos longitudinales y las son a su vez la quilla vertical, los longitudinales y las carlingas en el fondo, las esloras en la cubierta. Al cruzarse un elemento transversal con una longitudinal uno de ellos se interrumpe para dar paso al otro, por lo que este último se llama continuo y el primero intercostal. El forro y la cubierta están formados por chapas dispuestas en hiladas longitudinales llamadas tracas. Entre las tracas, algunas son de importancia especial: la quilla plana; las tracas adyacentes a la misma, llamadas de aparadura; la traca o tracas de pantoque; la traca más alta de forro o traca de cinta y la traca externas de cada cubierta, llamada trancanil. Cálculos estructurales por norma ABS Se realizó el cálculo estructural según la norma ABS, y así poder obtener todo el dimensionamiento de la embarcación. Todo esto se resume en los siguientes planos: plano estructural, plano de detalles y plano de cuadernas y mamparos. Para el cálculo de plancha de fondo, se usa el capítulo 1.8 de la norma ABS, como se muestra en lo siguiente: 37 Para el cálculo de puntales, se usa el capítulo 1.6 de la norma ABS, como se muestra en lo siguiente. 38 El diseño de los escantillones de las consolas según el capítulo 1.7.1, deberán ser, como mínimo según la siguiente figura: 39 De acuerdo a las normas ABS. Se obtuvo los siguientes perfiles con espesores: Casco PL. 6 mm Vagra PL. 110 mm Eslora PL. 100 mm Mamparos PL. 4 mm Ángulos L 50x50x6 mm Cuaderna L 50x50x6 mm Varengas L 50x50x6 mm Bao L 50x50x6 mm 40 Fig.5. PLANO DE MAMPAROS En estos planos de cuadernas y mamparos, nos indica la forma de los mamparos de colisión y a que distancia van estar separados respecto a popa y a proa. También nos indica cómo van estar unidas las cuadernas, varengas y baos. La distancia entre cuadernas es de 500 mm. El plano de mamparos nos muestra un total de 4 mamparos. En la parte de cubierta muestra 3 manhole. En el plano detalle según fig.2, nos muestra como está unido la varengas y la cuadernas mediante cartelas o consolas, según la norma ABS. 41 Fig.5. PLANO DE DETALLES 7) Generación y Validación del Software ANSYS Existen 3 grandes teoría para el análisis estructural, que son las siguientes: - CPT (Classical Plate Theory): Teoría clásica de placa. - FSDT (First order Shear Deformation Theory): Teoría de primer orden. - HSDT (Higher order Shear Deformation Theory): Teoria de orden superior. 42 Este Software está programado bajo la teoría de primer orden. x Aquí se muestra los siguientes pasos para validar la teoría: Displacement field (Campo de desplazamientos) Elasticity theory to get strains (Teoría de elasticidad para obtener def. Unit.) Law of Hook to get stresses (Ley de Hooke para obtener esfuerzos) The principle of virtual work, then use the discretization (Trabajo y Energia) Apply boundary conditions (Condiciones de contorno) Solve it! Fig.6. Teoria de Primer Orden 43 - Comparando resultados FSDT vs HSDT Se obtuvo los siguientes resultados: 4.1) Generación y Validación de una Viga Se modelara una viga bi-apoyada de longitud 5600mm, sujeta a una carga distribuida de 180 KN/m, utilizando 20 elementos en la viga. A sección recta de la viga en formato I con alas iguales de espesor de 2cm y ancho de 180cm, con un alma de espesor de 2cm y una altura de 900mm. El material es acero común. 44 DATOS DE LA VIGA: L= q= t= alto = ancho = densidad = 5600 180 20 900 1800 7850 mm KN/m mm mm mm Kg/m3 Primero: Definiremos el tipo de elemento con el comando “et, 1, beam3”. Segundo: Se define la sección con características indicadas arriba, con el menú preprocesor -> section -> beam -> common sections, donde se abre la siguiente tabla y la rellenamos con los datos obtenidos arriba. Fig.7. Sección de la Viga Tercero: Tomamos nota de los datos de área, momento de inercia con respecto al eje “y” y centro de área en el eje Z, que son: A 89200 mm2 Iyy 0.150*1011 mm4 XZ 450 mm Cuarto: Ahora se define las propiedades del material con el menú preprocesor -> material props -> material models -> structural -> linear elastic -> isotropic, donde definimos el módulo de elasticidad y el módulo de poisson, luego definimos la densidad como se muestra: 45 Quinto: Se agrega las constantes reales propias de la geometría de la sección, con el comando “r, 1, 89200, 0.15e11, 900”. Sexto: Se define una línea gráficamente con el menú, preprocesor -> modeling -> Keypoints -> in active CS, donde definimos las coordenadas (0, 0, 0) y (5600, 0, 0) como se muestran en la figura: Septimo: Luego dibujamos la línea que representara a la viga en el modelamiento, a partir de los keypoints establecidos con el menú preprocesor -> modeling -> créate -> lines -> in active coord. Octavo: Una vez definida la geometría, el tipo de elemento (que define la malla), el tipo de sección y sus constantes reales y por ultimo las propiedades del material procedemos al mallado para lo cual entramos al menú preprocessor -> meshing -> meshtool, una vez abierto ingresamos al element atributes y en la pantalla ingresamos lo siguiente: 46 Donde se puede ver que el tipo de material es el BEAM3, el número de material es 1 (donde se almacena los datos de elasticidad y densidad), el número de constantes reales es 1 (donde se almacenan los datos de la geometría de la sección) y el número de sección es 1(donde se almacenan los datos de las dimensiones de la sección). Noveno: Se define el tamaño y el número de elementos. Decimo: Se comienza con el enmallado, luego se coloca las cargas y restricciones entrando al menú preprocesor -> Loads -> new analysis como se ve en la imagen: Fig.8. Distribución de Fuerzas en la Viga Onceavo: Se define la dirección de la gravedad y por último se va al menú solution para poder correr el modelo, obteniendo los siguientes resultados: En esta tabla se pueden ver los valores de los esfuerzos máximos, las fuerzas cortantes y el momento flector en cada nodo. Comparándolos en una tabla y una gráfica de Excel con los resultados obtenidos por la teoría clásica de Resistencia de Materiales se obtienen los siguientes resultados: 47 - Para las fuerzas cortantes en la viga: NODO 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 V (N) TEO RESIS -654049.2 -601725.88 -549402.56 -497079.24 -444755.92 -392432.6 -340109.28 -287785.96 -235462.64 -183139.32 -130816 -78492.68 -26169.36 26153.96 78477.28 130800.6 183123.92 235447.24 287770.56 340093.88 392417.2 V (N) FEM -654040 -601720 -549400 -497070 -444750 -392430 -340100 -287780 -235460 -183130 -130810 -78485. -26162. 26162. 78485. 130810 183130 235460 287780 340100 392430 ERROR % 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.01 0.00 0.01 0.03 -0.03 -0.01 -0.01 0.00 -0.01 0.00 0.00 0.00 Se pueden observar en la gráfica de abajo ambos resultados ploteados V (N) vs nodos: 600000 400000 200000 0 -200000 V teo resis 0 5 10 15 20 25 V fem -400000 -600000 -800000 Viendo la gráfica directa del programa ANSYS obtenemos lo siguiente: 48 Fig.9. Diagrama de fuerzas Cortantes Por lo tanto, se logró validar el análisis numérico y teórico para una viga mediante el programa Ansys. 4.2) Generación y Validación de una Placa Una placa reforzada en dos direcciones con vigas puede ser analizada con un modelamiento de una malla para las vigas y elemento de placa. Si las vigas son uniforme y múltiples, tienen un alto grado de interacción, o si hay una pared doble, se puede utilizar el concepto de placa equivalente donde las vigas son repartida, por toda la superficie de la placa. Porque las vigas normalmente son diferentes en las dos direcciones. …….. (1) Esta ecuación es la que describe el comportamiento elástico de los paneles estructurales. Se hará el estudio de numérico – analítico de un panel estructural, de acuerdo a las curvas de Shade. Modelaremos una placa de longitud de 1900 mm y un ancho de 1500 mm con un espesor de 8mm. Se distribuye en toda la placa un carga de 0.1 N/mm 2 y se utilizara refuerzos ½ perfil IPE 180. 49 Fig.10. Modelación del panel Luego será simulada con la condición con todos los bordes apoyados. Fig.11. Diagrama de Esfuerzo en el panel Ahora analizaremos este panel según la teoría de cálculo de Shade. Nos encontramos en el caso Tipo A, para planchas con refuerzos cruzados todos los que se requieran. Fig.12. Caso Tipo A Una vez entendido esto pasare a identificar la simbología expuesta al costado de las gráficas. a (b) = Largo (ancho) de la plancha. P = Presión uniforme. 50 Sa (Sb) = Espaciamiento de los refuerzos largo (corto). Ina (Inb) = Momento de Inercia, incluyendo la longitud efectiva de la plancha asociada de refuerzos largos repetitivos. Para la condición de todos los bordes apoyados se usara la siguiente tabla. Fig.13. Curvas de Shade Esta tabla se utiliza para calcular los desplazamientos de planchas solo en el punto medio, de manera de acceder a la tabla es la siguiente, tenemos que entender que existen 4 tipos de condiciones, en la cuales no encontramos en el Tipo A. En el estudio de paneles estructurales según los cálculos de Shade se obtuvieron los siguientes resultados. METODO ANALITICO Esfuerzo Trnsv. Placa 157.6 MPa Esfuerzo Trnsv. refuerzo 142.9 MPa METODO NUMERICO Esfuerzo Trnsv. Placa 158 MPa. Esfuerzo Trnsv. refuerzo 143.7 MPa. Por lo tanto, se logró validar el análisis numérico y teórico para una placa mediante el programa Ansys. 51 8) Carga Rodante La ecuación diferencial (1) ha sido escrita con referencia a la presión normal distribuida como una función continua en la placa. Fig.14. Vista Longitudinal Esta barcaza fluvial va a transportar como máximo 4 camionetas con carga máxima. El peso bruto de cada camioneta es de 4 Toneladas, de la cuales será distribuida en sus 4 llantas, cada llanta soportara 1 tonelada el cual será transmitido a la cubierta de la embarcación. Fig.15. Distribución de Carga El área proyectada de cada llanta a la cubierta en la condición de carga máxima, es un rectángulo de 350mm x 250mm aproximadamente, en esta área será distribuida la carga de cada llanta. 9) Modelación Estructural Modelos estructurales son considerados como modelos de la plantilla de cualquier estructura del edificio. Se ofrece una visualización exacta de la estructura del 52 edificio propuesto y modelado en 3D es únicamente el proceso de creación de una representación matemática, armazón de alambre de cualquier objeto en tres dimensiones utilizando un software especializado. Se utiliza para crear el diseño en 3D precisa y rápida los modelos. 3D de modelado estructural es una parte de los servicios de ingeniería estructural y puede ser identificado como el proceso de elaboración de una representación matemática, armazón de alambre de cualquier objeto 3D utilizando técnicas de computación especializados. Se utiliza en el desarrollo de la construcción de modelos excelentes, personalizados modelos estructurales de viviendas y el modelado de estructuras para distintos tipos de edificios, tales como edificios residenciales, comerciales e industriales. Fig.16. Modelado Estructural 3D La creación de un equilibrado modelo estructural 3D es una colección de los pasos detallados, a partir de una Sketch-Up de objeto. Modelación estructural tiene que ver con la creación de modelos de construcción 3D precisos y fiables según los criterios requeridos. Modelado en 3D de estructuras se utiliza para crear modelos exactos de construcción estructural de elementos constructivos y estructurales diferentes modelos en 3D de diferentes estructuras de los edificios juegan un papel decisivo para los constructores, los contratistas y arquitectos para construir sus estructuras de los edificios propuestos con mayor eficacia. De modelado estructural ayuda a encontrar volúmenes específicos del componente y que permite a los constructores o contratistas que entender acerca de la cantidad de materiales y componentes serán necesarios. También es útil para la toma de vistas en 3D exteriores, vistas en 3D de interiores con las especificaciones de dimensiones, los modelos de paisaje de diseño estructural y los modelos de ecuaciones estructurales. 53 El modelo estructural para la barcaza fluvial será el modelo de tipo superficie, ya que los espesores calculados según la norma se dan al momento de la simulación. Fig.17. Modelado Isométrico 10) Simulación Numérica de la Barcaza El modelo estructural final es exportado al software ANSYS con formato IGES, para poder comenzar con la simulación numérica. Los métodos numéricos en Ingeniería proporcionan una herramienta que permite traducir modelos matemáticos en procedimientos computacionales, cuyos resultados pueden ser contrastados con las soluciones analíticas, en los casos en los que éstas existan. En este trabajo se muestra una forma de introducir al estudiante de ingeniería en el estudio de las ecuaciones diferenciales en derivadas parciales y su modelado numérico. En particular la ecuación de difusión será el objeto de estudio. La elección de la misma se efectúa basándose en sus múltiples aplicaciones en problemas de la Mecánica. Usando el método de diferencias finitas para determinar la solución numérica de la ecuación de difusión, se abordarán básicamente dos planteos: método explícito y método implícito. Para ambos procedimientos es necesario diseñar algoritmos computacionales. Este trabajo pretende mostrar cómo contribuir a la formación matemática del ingeniero mediante el uso de métodos numéricos. Cómo éstos introducen a los alumnos en la programación para la solución de Ecuaciones Diferenciales Parciales. Finalmente muestra la comparación de los resultados obtenidos con el 54 uso de distintos métodos numéricos y la solución analítica hallada en los casos posibles. Estos son los 7 pasos para tener la solución numérica completa, de acuerdo a las condiciones dadas en el problema. Aquí están incluidas las 3 etapas para toda simulación numérica. - Pre-procesamiento: El primer paso en FEA, pre-procesamiento, es construir un modelo de elementos finitos de la estructura a ser analizada. En muchos paquetes de FEA se requiere de la entrada de una descripción topológica de las características geométricas de la estructura. Ésta puede ser 1D, 2D, o 3D. El objetivo principal del modelo es replicar de manera realista los parámetros importantes y características del modelo real. La manera más sencilla para conseguir similaridad en el análisis es utilizar planos pre existente, modelos CAD, o datos importados de un ambiente FEA. Una vez se ha creado la geometría, se utiliza un procedimiento para definir y dividir el modelo en "pequeños" elementos. En general, un modelo de elementos finitos está definido por una malla, la cual está conformada por elementos y nodos. Los nodos representan puntos en los cuales se calcula el desplazamiento (análisis estructural). Los paquetes de FEA enumeran los nodos como una herramienta de identificación. Los elementos están determinados por conjuntos de nodos, y definen propiedades localizadas de masa y rigidez. Los elementos también están definidos por la numeración de la malla, la cual permite referenciar la correspondiente deflexión o esfuerzo (en análisis estructural) para una localización específica. 55 - Módulo de análisis: En la siguiente etapa en el proceso de análisis de elementos finitos se lleva a cabo una serie de procesos computacionales que involucran fuerzas aplicadas, y las propiedades de los elementos de donde producir un modelo de solución. Tal análisis estructural permite la determinación de efectos como lo son las deformaciones, estiramiento o estrés que son causados por fuerzas estructurales aplicadas como lo son la fuerza, la presión y la gravedad. - Post-procesamiento: Estos resultados entonces pueden ser estudiados utilizando herramientas visuales dentro del ambiente de FEA para ver y para identificar completamente las implicaciones del análisis. Herramientas numéricas y gráficas permiten la localización precisa de información como esfuerzos y deformaciones a ser identificadas. 11) Generación de Malla Una malla es un conjunto de celdas contiguas que permite representar en forma discreta el dominio de un problema a resolver numéricamente. La confiabilidad de la solución numérica obtenida depende de la calidad de la malla generada. Los criterios de calidad a usar dependen del método numérico escogido y del tipo de problema a resolver. Fig.18 Enmallado de Cuadernas 56 En esta presentación se discutirá sobre los métodos para modelar y describir la geometría del dominio, los tipos de mallas más usados, algoritmos para generar mallas, mostrando con más detalle uno para generar mallas mixtas, criterios de refinamiento y de calidad usados, el estado de la investigación actual en el área. Las mallas geométricas se pueden clasificar en estructuradas y no estructuradas. Las mallas estructuradas se caracterizan estar compuestas de celdas de un tamaño similar y del mismo tipo, por ejemplo, triángulos o rectángulos, en dos dimensiones, tetraedros o hexaedros en tres dimensiones. Estas mallas son fáciles de generar pero no permiten modelar eficientemente problemas geométricamente complejos. Por otro lado, las mallas no estructuradas permiten el uso de celdas de distinto tipo y/o de diferente tamaño. Estas mallas requieren de algoritmos y estructuras de datos más complejos que las anteriores, pero permiten modelar geometrías complejas y optimizar el número de celdas usadas de acuerdo a la necesidad de la aplicación. Actualmente, la generación de una malla adecuada en dos dimensiones se considera un problema resuelto, pero no así en tres dimensiones. En tres dimensiones se trabaja en el diseño algoritmos eficientes, tanto secuenciales como paralelos, en todas las etapas del proceso de generación de una malla de calidad, para poder resolver cada vez problemas más complejos. Fig.19 Algoritmos de Malla 57 12) Condiciones de Contorno Se aplicó las siguientes condiciones de contorno a la embarcación. - Presiones hidrostáticas, teniendo en cuenta el calado de diseño para su máxima carga. Fig.20 Modelo Teórico - Distribución de presiones Fig.21 Modelo Numérico – Distribución de presiones - Modelando los resortes en los nodos donde se intersecan los refuerzos con los mamparos. Estos resortes nos permiten restringir el movimiento en el eje vertical y también rigidiza la embarcación para que pueda rodar. En el programa ANSYS, se le conoce como el tipo de elemento “Combin 14”. 58 Fig.22 Modelamiento de Resortes - Modelamiento de cargas de las 4 camionetas con carga máxima. El peso bruto de cada camioneta era de 4 Toneladas. Esta carga es distribuida en la parte apoyada de las llantas sobre la cubierta. En cada camioneta hay 4 apoyos, en total hay 16 apoyos, cada uno soporta aproximadamente 10 KN. Fig.23 Distribución de Carga 13) Resultados Se muestra los siguientes diagramas de esfuerzos y deformaciones en toda la embarcación (comportamiento estructural de la barcaza). 59 Fig.24 Diagrama de Esfuerzos Haciendo un corte en la parte de la sección media, vemos como se distribuye los esfuerzos en el casco del fondo por medio de los puntales. Fig.25 Distribución de Esfuerzos En el estudio de la embarcación fluvial se obtuvo como esfuerzo máximo de Von Misses de 203.18 MPa. en las esquinas de algunos perfiles por ser puntos duros (concentración de esfuerzos) y en un análisis particular de los perfiles se soluciona el problema dando cierto radio de curvatura en las esquinas, pero el esfuerzo máximo que debemos resaltar es de 70 MPa. como máximo en un análisis global. Por lo tanto, nuestros resultados están por debajo del esfuerzo de fluencia del material (Acero estructural A36). 60 14) Conclusiones Se determinó la teoría en la cual trabaja el software ANSYS. Se validó el análisis estructural de la placa ortotropica y se determinó el número de nodos mínimo en cada refuerzo, para los futuros enmallados. Se localizó y calculo el esfuerzo de trabajo máximo de la barcaza. Se logró saber cómo es el comportamiento estructural de la embarcación fluvial. 15) Bibliografía Libro de estructura del Ing. Williams (Ex catedrático de la UNI). Zienkiewicks, O.C. y Taylor, R.L. (1994), “El Metodo de los Elementos Finitos”(FEM). American Bureau of Shipping (2008), “Rules for Building and Classing, Steel Vessels, Part 3, Hull Construction and equipment”, American Bureau of Shipping, ABS. “Geometry of Ships” Autor: John Letcher. “Strength of Ships and Ocean Structures” Autor: Alaa Mansour & Don Liu. ABS, “Safehull Finite Element Analysis of Hulls Structures” 61 https://sites.google.com/site/danieljimenez94/mywork/tutorialanalisisestatic odedeformacionesensolidworks Tutoriales de ansys. (http://www.kxcad.net/ansys/) http://www.padtinc.com/blog/the-focus/ansys-mesh-connections-anothertool-for-meshing-surface-assemblies. http://www.simutechgroup.com/FEA/fea-tips-tricks-macros-buoyance floating.html. 62
© Copyright 2024