ANEXOS ASESOR: Dr. Rosendo Franco Rodríguez CO-ASESOR: Dr. Quino Valverde Guzmán Lima, Junio del 2016 PONTIFICIA UNIVERSIDAD CATÓLICA DEL PERÚ FACULTAD DE CIENCIAS E INGENIERÍA Tesis para optar el Título de Ingeniero Mecánico, que presenta el bachiller: ANÁLISIS ESTRUCTURAL Y MODAL DE LA HÉLICE DE UN AEROGENERADOR DE 3kW DE POTENCIA MEDIANTE SIMULACIÓN NUMÉRICA COMPUTACIONAL RENATTO MARCELLO YUPA VILLANUEVA 1 ANEXO 1 GUÍA GENERAL DE SIMULACIÓN NUMÉRICA En este capítulo se explicará detalladamente el procedimiento general a seguir para las simulaciones numéricas a realizar. Los pasos que se detallan a continuación, son comunes para la simulación de flujo externo en ANSYS CFX, para la simulación estructural estática en ANSYS MECHANICAL y para la simulación estructural dinámica en ANSYS MECHANICAL. 1. Definición del problema Para realizar el análisis de sistema físico usando simulación numérica, es un requerimiento definir el problema, es decir fijarnos los objetivos. Es de vital importancia este paso, ya que condiciona totalmente a los tres pasos siguientes. Para poder definir el problema debemos enfocarnos en aspectos fundamentales, que los mencionaré a continuación.  Delimitar correctamente el sistema físico a estudiar.  Definir qué variables del problema son las que se quiere hallar.  Incorporar las condiciones conocidas del sistema físico, es decir, toda la información que se conozca sobre el sistema físico delimitado (datos de entrada). Una vez que se definió correctamente el problema, se puede proceder a los siguientes pasos. Figura A1.1. Procedimiento general de simulación numérica. 2. Pre-procesamiento La etapa del pre-procesamiento es la primera etapa de la simulación numérica que se realiza en el computador, específicamente en el software utilizado, el cual para fines de la presente tesis es ANSYS. En esta etapa, se configuran todos los parámetros necesarios y modelos físicos a utilizar. Es de vital importancia la correcta realización del pre-procesamiento en concordancia con el objetivo de la simulación, ya que, de no ser así, los resultados obtenidos en la etapa de solución serán incorrectos (a pesar que la configuración del pre-procesamiento sea incorrecta, la simulación puede llegar a converger totalmente; sin embargo, los resultados no serán correctos). Por tal motivo, es necesario que el usuario adquiera conocimientos previos del software y de la teoría del método de los elementos finitos, para poder discernir con criterio entre lo correcto e incorrecto durante la 2 simulación. La etapa del pre-procesamiento se subdivide en cinco pasos consecutivos, los cuales se explican a continuación. 2.1. Definición de datos de entrada En esta etapa se recopila toda la información que se tiene sobre el problema a estudiar, ya que serán los datos de entrada para la simulación numérica. Es importante resaltar que se debe discernir de entre todos los datos, los datos relevantes. Asimismo, en esta etapa se deben hacer todas las simplificaciones posibles del sistema físico a estudiar, siempre y cuando los resultados finales de la simulación se mantengan de acorde con la realidad. 2.2. Diseño CAD del modelo: Definición del dominio de interés El diseño CAD del modelo se realiza tomando en consideración el dominio de interés (parte del sistema físico delimitado donde el software resuelve las ecuaciones gobernantes). Es importante la correcta definición del dominio de interés, ya que únicamente sobre este dominio se realiza la simulación numérica. La definición del dominio depende principalmente del sistema físico delimitado y del objetivo del problema. El diseño CAD del modelo se realiza con ayuda de un software CAD. ANSYS trae algunos modeladores CAD 3D tales como SpaceClaim o Design Modeler, los cuales son muy sofisticados. La elección del modelador CAD depende del dominio del usuario. Sin embargo, para fines de la presente tesis, se eligió usar Autodesk Inventor como software CAD, ya que presenta grandes facilidades al dibujante al modelar en 3D. No existe ninguna restricción al crear el modelo 3D en Autodesk Inventor y exportarlo a ANSYS. Al realizar el modelo 3D, se recomienda simplificar al máximo la geometría del dominio de interés. De esta manera, se aligera los posteriores procesos durante la simulación como es la etapa de solución, por ejemplo. Asimismo, una vez que el modelo 3D es exportado a ANSYS, es recomendable “limpiar” la geometría del mismo, es decir, corregir errores de dibujo tales como cruces de aristas, etc. Esto se realiza con el fin de reducir la complejidad de las mallas y, al mismo tiempo, hacer de la geometría lo más cercana a la realidad. Asimismo, la facilidad en el manejo del software CAD es imprescindible en esta etapa. 2.3. Parámetros de mallado Luego de tener definido y creado el dominio de interés, se procede al mallado. Se define mallado como el proceso de dividir el dominio de interés en subdominios mucho más pequeños (elementos finitos), donde a su vez se resolverán las ecuaciones gobernantes del modelo. Sin embargo, antes de realizar el mismo, se debe definir los parámetros del mallado, tales como tipo de mallado (geometría 3 del elemento a usar), método de mallado a utilizar, etc. En general, durante el mallado se debe cumplir con las siguientes indicaciones:  La selección de la forma del elemento (tetraédrico, hexagonal, etc.) debe ser completamente entendido por el solver del ANSYS.  Los elementos finitos deben cubrir al dominio de interés en su totalidad.  Todos los elementos finitos deben estar correctamente interconectados.  Por ningún motivo debe existir elementos superpuestos.  La distribución de los elementos finitos debe ser la adecuada para capturar las características geométricas. Los principales parámetros de mallado (método de mallado y tipo de elemento) se explican a continuación: Tipos de elementos finitos Existen varios tipos de elementos finitos, entre los más importantes se pueden mencionar:  Elementos hexaédricos: Elementos de 6 caras y 3 secciones de cuadriláteros. Pueden ser de 8 nodos para elementos lineales, y 20 nodos para elementos cuadráticos. En la Figura A1.2, se puede apreciar lo mencionado. Figura A1.2. Elementos hexaédricos: lineal de 8 nodos (izquierda) y cuadrático de 20 nodos (derecha).  Elementos tetraédricos: Elementos de 4 caras triangulares y 6 secciones triangulares. Cada sección del tetraedro se construye haciendo pasar un plano desde un borde hasta el punto en el borde opuesto más cercano. Pueden ser de 4 nodos para elementos lineales y 10 nodos para elementos cuadráticos. En la Figura A1.3, se puede apreciar lo antes mencionado. Figura A1.3. Elementos tetraédricos: lineal de 4 nodos (izquierda) y cuadrático de 10 nodos (derecha).  Elementos prismáticos: Elementos de 5 caras, 2 de ellas triangulares y 3 de ellas cuadriláteros. Pueden tener 6 nodos para elementos lineales y 15 nodos para elementos cuadráticos. En la Figura A1.4, se puede apreciar lo antes mencionado. 4 Figura A1.4. Elementos prismáticos: lineal de 6 nodos (izquierda) y cuadrático de 15 nodos (derecha). Finalmente, es importante precisar que la principal diferencia entre los elementos lineales y los elementos cuadráticos es la cantidad de nodos que poseen los elementos. A mayor cantidad de nodos por elemento, se consigue una mejor interpolación de resultados de las variables evaluadas en el elemento finito. Métodos de mallado en ANSYS El software ANSYS tiene diversos métodos de mallado 3D, los principales se detallan a continuación:  Mallado automático: Es un método de mallado que combina el mallado tetraédrico patch conforming con el mallado de barrido, los cuales se explican posteriormente. Este método identifica todas las partes que pueden ser malladas con el método de barrido y procede a mallar de tal manera. Luego, todas las partes restantes son malladas con el método de tetraedros patch conforming.  Mallado tetraédrico: Es un método de mallado que genera elementos tetraédricos, para los cuales existen dos algoritmos de mallado:  Patch Conforming: Este algoritmo es preferentemente usado para geometrías limpias, es decir, geometrías reparadas y simplificadas. Al usar este algoritmo, el mallado se realiza a partir de las caras y bordes, finalizando en los volúmenes grandes. En términos coloquiales, al usar este algoritmo se consigue que los elementos finitos se adapten a la geometría, y no al revés.  Patch Independent: En este algoritmo, en contraste con el anterior descrito, primero se genera una malla desde el volumen central y hacia los bordes y caras, de tal manera que no necesariamente se respeta la geometría, en los detalles. En términos coloquiales, al usar este algoritmo se consigue que la geometría se adapte a la forma de los elementos finitos, y no al revés.  Mallado hexagonal: Es un método de mallado que genera elementos hexaédricos. La malla hexaédrica se caracteriza por presentar un número reducido de elementos, con lo cual se mejora el tiempo de convergencia del modelo. Asimismo, se consigue una reducción del error numérico con lo que se mejora la precisión del cálculo. Estas comparaciones son válidas 5 para mallas tetraédrica y hexaédrica de primer orden (lineales). Existen dos métodos de mallados hexaédricos, los cuales se describen a continuación.  Hexaedros dominantes: El método hexaedro dominante se caracteriza por crear una malla de elementos hexaédricos dominante a los elementos tetraédricos. En términos coloquiales, este método de mallado crea elementos hexaédricos en su mayoría (donde la geometría del modelo permite la creación de elementos hexaédricos), y crea elementos tetraédricos (donde la geometría no permite que se creen elementos hexaédricos).  Mallado de barrido: El mallado de barrido genera únicamente elementos hexaédricos los cuales se forman a partir de elementos cuadriláteros que se forman en la superficie a partir de la cual se barre. Posteriormente, se genera una malla 3D de elementos hexaédricos que se barren desde la superficie inicial. Los elementos de la superficie de llegada deben tener la misma topología que los elementos de la superficie inicial. Cabe mencionar que circunstancialmente se pueden crear elementos prismáticos en el barrido. La siguiente figura muestra lo antes mencionado. Estos métodos combinan los diferentes tipos de elementos finitos, de acuerdo a las necesidades del problema. Cada método presenta ventajas y desventajas frente a los otros métodos. La elección depende del usuario y del tipo de problema a resolver. Asimismo, además de configurar el método de mallado a utilizar y los tipos de elementos, se debe definir las configuraciones globales y locales del mallado. Configuraciones globales Las configuraciones globales de malla se utilizan para el ajuste global en la estrategia de mallado, que incluye todas las funciones de sizing, inflation, smoothing, parámetros de entrada, etc. Son muy útiles para resolver las principales características del modelo con un mínimo de insumos. A continuación, una breve descripción de los principales parámetros en la configuración global:  Physics based settings: Se establece la física de preferencia, es decir, si se trata de un problema CFD o un problema de mecánica de sólidos.  Global sizing controls: Se establecen configuraciones globales de tamaño de elemento.  Inflation: Inflation es una herramienta que sirve para crear un mallado por capas cercano a una superficie o superficies de interés. Esta herramienta es usada principalmente para problemas de interacción fluido-estructura donde se requiere información en la capa límite del fluido. En este punto, se establece la configuración global de la herramienta en el mallado. En la siguiente figura se aprecia con detalle lo antes mencionado. 6  Patch confirming options: Se establecen configuraciones globales sobre el algoritmo patch conforming para mallados tetraédricos.  Statistics: Es una herramienta que sirve para visualizar las estadísticas de mallado obtenidas luego de generar la malla. Posteriormente, se explica a detalle este punto. Configuraciones locales En las configuraciones locales se establece configuraciones específicas de malla para ciertas zonas o superficies de interés que las requiera. A continuación, una breve descripción de los principales parámetros en la configuración local:  Sizing: Es una herramienta que sirve para definir tamaños de elementos finitos para zonas determinadas de la malla global.  Contact sizing: Es una herramienta que sirve para definir tamaño de elementos para zonas de contacto.  Mapped face meshing: Es una herramienta que sirve para estructurar elementos en una cara determinada.  Inflation: Sirve para crear elementos por capas, como se explicó anteriormente. Se configura esta herramienta en la configuración local para superficies específicas. Una vez definido todos estos parámetros, se puede pasar a generar la malla. 2.4. Mallado El siguiente paso es generar el mallado. Con toda la configuración realizada en el paso anterior, solo queda por generar rápidamente el mallado. Es importante mencionar en este punto, que se puede visualizar y evaluar el tipo, la calidad y las propiedades del mallado, luego de la realización del mismo. En caso la malla no sea la adecuada, se puede refinar y dar ajustes para llegar a la calidad de malla adecuada. A continuación, se describe las estadísticas de mallado más importantes.  Element quality: La calidad de elemento proporciona una calidad métrica la cual varía entre 0 y 1. Esta métrica se basa en la relación entre el volumen y la longitud del borde para un elemento dado. Un valor de 1 indica un cubo o un cuadrado perfecto, mientras que un valor de 0 indica que el elemento tiene volumen cero o negativo.  Aspect ratio: La relación de aspecto es una proporción entre el ancho y la altura para un elemento. La relación de aspecto ideal es 1 (para elementos 2D es un triángulo equilátero o un cuadrado).  Jacobian ratio: La relación de jacobiano se calcula y se prueba para todos los elementos excepto triángulos y tetraedros que: 7  Son lineales (no tienen nodos intermedios).  Tienen nodos intermedios, los cuales están perfectamente centrados. Una relación alta indica que el espacio de elementos intermediarios de mapas y el espacio real está convirtiendo, computacionalmente, al modelo en poco fiable.  Parallel Deviation: Es una medida de qué tan inclinados están 2 lados paralelos y opuestos de elemento con respecto a otros 2 lados paralelos y opuestos en el mismo elemento. La siguiente figura ayuda a concebir este concepto. Figura A1.5. Desviación paralela para un elemento 2D cuadrilátero. La desviación paralela ideal es 0º, sin embargo en la imagen se muestran desviaciones paralelas de 70º, 100º, 150º y 170º.  Skewness: La asimetría es una de las principales medidas de la calidad de una malla. Esta mide qué tan cerca a la ideal (triángulo equilátero o cuadrado) está una cara del elemento.  Orthogonal Quality: La calidad ortogonal es una medida de la calidad de malla que se calcula de la siguiente manera:  Se traza un vector desde el centroide del elemento hacia el centroide de cara una de las caras del elemento.  Se traza un vector normal a cada cara del elemento.  La relación entre qué tan parecidos son dichos vectores es la calidad ortogonal. 2.5. Parámetros de control de la simulación (condiciones de borde) Los parámetros de control son los que determinan el tipo de modelamiento, los alcances que tiene la simulación y el grado de precisión del cálculo. En otras palabras, los parámetros de control de la simulación son las condiciones de frontera o borde que se aplican tanto en los límites del dominio de interés, como dentro del mismo. Condiciones de borde más utilizadas en problemas CFD  Material del dominio de interés  Tipo de turbulencia  Condiciones de ingreso, salida, etc.  Velocidad inicial  Temperatura inicial  Tipo de condición de borde, etc. 8 Condiciones de borde más utilizadas en problemas de mecánica de sólidos  Presión  Momentos  Fuerzas puntuales o distribuidas  Velocidad  Contactos  Apoyos  Condiciones de temperatura La correcta configuración de los parámetros de control de la simulación van de acuerdo al objetivo del problema e influyen totalmente en la simulación. Por ello, esta etapa debe ser de sumo cuidado. Cada parámetro elegido debe ser estudiado previamente. Como se mencionó anteriormente, que la simulación converja no asegura la correcta solución del problema. 3. Solución Transcurrida la etapa del pre-procesamiento, donde se configuró todo el modelo, se procede a la solución. En esta etapa, el proceso seguido es el siguiente:  Las ecuaciones diferenciales que rigen el modelo se integran sobre todos los elementos finitos o volúmenes de control (en caso de usar el método de los elementos finitos o el método de los volúmenes finitos) en toda la región de interés.  Luego, el sistema de ecuaciones integrales se convierten en un sistema de ecuaciones algebraicas mediante la generación de un conjunto de aproximaciones para los términos de las ecuaciones integrales.  Finalmente, el sistema de ecuaciones algebraicas se resuelven de manera iterativa. Para problemas de CFD: Durante la última sub-etapa, se menciona que los sistemas de ecuaciones algebraicas se resuelven mediante método iterativos. Al usar dicho método, el solucionador calcula de forma inherente los errores relativos de la iteración. Conforme se acerca más a la realidad, el error relativo es muy pequeño. Entonces, para CFD, se dice que la simulación converge cuando los errores relativos de las variables monitoreadas, llegan al valor de 1 x 10-4, por defecto para ANSYS CFX. El usuario también puede monitorear las variables creadas durante la etapa de pre-procesamiento, puede modificar los criterios de para de la simulación (máximo número de iteraciones, error relativo permisible, etc.), en el caso de ANSYS CFX. 9 4. Post-procesamiento En esta última etapa, el objetivo es obtener la información (cualitativa y cuantitativa) del dominio de interés simulado, para poder mostrarla de manera atractiva, posteriormente. Para ello, el solucionador (solver) generó un archivo de solución al converger el modelo simulado y lo transfiere al post- procesamiento. En esta etapa, se puede realizar procesamiento de gráficos interactivos, generación de reportes, gráficos y tablas, etc. En general, la etapa del post-procesamiento en ANSYS presenta las siguientes características:  Una interfaz gráfica de usuario que incluye un panel visor en el que se representa toda la producción gráfica del post-procesamiento de ANSYS.  Soporte para una variedad de objetos gráficos y geométricos para la creación de parcelas de post-procesamiento, visualización de la malla, y para definir las ubicaciones para el cálculo cuantitativo.  Variables escalares y vectoriales definidas por el usuario.  Congelación de variables.  Generación de un reporte de resultados.  Generación de una variedad de objetos gráficos donde la visibilidad, la transparencia, el color y la visualización de línea / caras puede ser controlado. Asimismo, es importante mencionar en esta etapa que el usuario debe analizar con precaución y criterio ingenieril los resultados obtenidos en las simulaciones numéricas. Para tal caso, normalmente, se estila realizar alguna comprobación de los resultados, mediante cálculos teóricos simples o experimentación. 10 ANEXO 2 TEORÍA BÁSICA DE MATERIALES COMPUESTOS 1. Fundamento teórico 1.1. Definición Es llamado material compuesto que está compuesto de dos o más fases. Este compuesto posee propiedades mecánicas superiores a la de los materiales constituyentes. Generalmente una de las fases es la más rígida y fuerte, la cual es llamada refuerzo, mientras que la fase menos rígida y débil es llamada matriz. Figura A2.1. Material compuesto 1.2. Anisotropía de los materiales compuestos Los materiales compuestos poseen ventajas únicas sobre los materiales monolíticos, tales como alta resistencia, alta rigidez, baja densidad. La base del desempeño superior estructural de los materiales compuestos radica en la alta resistencia específica, la alta rigidez específica y en el carácter anisotrópico y heterogéneo del material. Esto último proporciona al material compuesto varios grados de anisotropía para la configuración óptima del material. Algunas propiedades del material, tales como la densidad, calor especifico, no poseen ninguna direccionalidad asociada con estas y son descritas por una cantidad escalar para materiales isotrópicos y anisotrópicos. Por otro lado, propiedades tales como la rigidez, coeficiente de Poisson, resistencia, conductividad térmica y conductividad eléctrica están asociadas con la dirección y están en función de la orientación en materiales anisotrópicos. Las mayores diferencias ocurren entre las propiedades en la dirección longitudinal (fibras) y transversal (normal a las fibras). En la siguiente figura A2.2 los subíndices 1, 2, y 3 indican respectivamente las propiedades que presenta el material compuesto en los ejes longitudinales (dirección de las fibras), transversales (normal a las fibras), y fuera del plano. 11 Figura A2.2. Orientación de los ejes en función de las fibras. 1.3. Características y configuración de la lámina y laminado Una lámina es una capa plana (o curvada) de una matriz con fibras unidireccionales o tejido de fibras. La lámina es un material ortotrópico, el cual se define como aquel material que presenta tres ejes perpendiculares entre sí, y presenta propiedades constantes definidas en cada una de estas tres direcciones (no necesariamente iguales), con ejes principales en la dirección de las fibras (longitudinal), transversal a las fibras (plano transversal) y normal al plano de la lámina. Estos ejes principales se les designa como 1,2 y 3 respectivamente y se ilustra a continuación. Figura A2.3. Ejes principales de un laminado. Un laminado consiste de dos o más laminas unidireccionales u orientadas en distintas direcciones unidas unas a otras. 12 Figura A2.4. Secuencia de laminado. Los laminados compuestos se designan indicando el número, tipo, orientación, y secuencia de apilamiento de las láminas. A continuación se muestran algunos ejemplos de designaciones de laminados.  Unidireccionales: [0 / 0 / 0 / 0 ] = [04]  Cruzado simétrico [0 / 90 / 90 / 0 ] = [0 / 90]𝑠  Simétrico orientado [+45 / −45 / −45 / +45 ] = [±45]𝑠 [60/−60/60/−60/−60/60/−60/60] = [±60]2𝑠  Multidireccional [0/0/45/−45/0/0/0/0/−45/45/0/0] = [02/±45/02]𝑠  Híbrido [0𝐺/45𝐾/−45𝐾/90𝐶/−45𝐾/45𝐾/0𝐺] = [0𝐺/±45/90̅̅̅̅ 𝐶]𝑠 Donde los subíndices y símbolos significan losiguiente: Subíndice = Múltiplo de láminas o número de laminas s= Secuencia simétrica Barra superior denota que el laminado es simétrico respecto al plano medio de dicha lámina. En el caso de laminados híbridos el índice superior G, K, y C denotan fibras de vidrio, Kevlar (aramida), y carbono (grafito), respectivamente. 13 1.4. Propiedades básicas de la lámina La lámina unidireccional es considerada la parte de construcción básica de cualquier laminado o estructura compuesta. Las propiedades necesarias para su análisis y diseño de estas estructuras son las propiedades básicas de la lámina. Para el algoritmo que se desarrollará más adelante se desprecia el esfuerzo o deformación que se pueda generar en el eje z, normal al plano del laminado, las cuales son algunas de las consideraciones en las cuales está basado el modelo matemático. Es por ello que las propiedades son las siguientes:  𝐸1, 𝐸2 = 𝑀ó𝑑𝑢𝑙𝑜 𝑑𝑒 𝑦𝑜𝑢𝑛𝑔 𝑎 𝑙𝑜 𝑙𝑎𝑟𝑔𝑜 𝑑𝑒 𝑙𝑎𝑠 𝑑𝑖𝑟𝑒𝑐𝑐𝑖𝑜𝑛𝑒𝑠 𝑝𝑟𝑖𝑛𝑐𝑖𝑝𝑎𝑙𝑒𝑠 𝑑𝑒𝑙 𝑝𝑙𝑎𝑛𝑜.  𝐺12 = 𝑀ó𝑑𝑢𝑙𝑜 𝑑𝑒 𝑐𝑜𝑟𝑡𝑒 𝑒𝑛 𝑒𝑙 𝑝𝑙𝑎𝑛𝑜 1 − 2 (𝑖𝑔𝑢𝑎𝑙 𝑎𝑙 𝐺21).  𝜐 12 , 𝜐21 = 𝐶𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑃𝑜𝑖𝑠𝑜𝑛 (𝑒𝑙 𝑝𝑟𝑖𝑚𝑒𝑟 𝑠𝑢𝑏𝑖𝑛𝑑𝑖𝑐𝑒 𝑖𝑛𝑑𝑖𝑐𝑎 𝑙𝑎 𝑑𝑖𝑟𝑒𝑐𝑐𝑖𝑜𝑛 𝑑𝑒 𝑙𝑎 𝑐𝑎𝑟𝑔𝑎 𝑦 𝑒𝑙 𝑠𝑒𝑔𝑢𝑛𝑑𝑜 𝑑𝑒𝑛𝑜𝑡𝑎 𝑙𝑎 𝑑𝑖𝑟𝑒𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑙𝑎 𝑑𝑒𝑓𝑜𝑟𝑚𝑎𝑐𝑖ó𝑛, 𝑒𝑠𝑡𝑜𝑠 𝑑𝑜𝑠 𝑣𝑎𝑙𝑜𝑟𝑒𝑠 𝑛𝑜 𝑠𝑜𝑛 𝑖𝑔𝑢𝑎𝑙𝑒𝑠).  𝐹1𝑡, 𝐹2𝑡 = 𝑅𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎 𝑎 𝑙𝑎 𝑡𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑎 𝑙𝑜 𝑙𝑎𝑟𝑔𝑜 𝑑𝑒 𝑙𝑎𝑠 𝑑𝑖𝑟𝑒𝑐𝑐𝑖𝑜𝑛𝑒𝑠𝑝𝑟𝑖𝑛𝑐𝑖𝑝𝑎𝑙𝑒𝑠 𝑒𝑛 𝑒𝑙 𝑝𝑙𝑎𝑛𝑜.  𝐹12 = 𝑅𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎 𝑎𝑙 𝑐𝑜𝑟𝑡𝑒 𝑒𝑛 𝑒𝑙 𝑝𝑙𝑎𝑛𝑜 1 − 2(𝑒𝑠𝑡𝑎 𝑣𝑎𝑙𝑜𝑟 𝑒𝑠 𝑖𝑔𝑢𝑎𝑙 𝑎 𝐹21). 1.5. Propiedades de los materiales constituyentes 1.5.1. Fibras Existe una gran variedad de fibras como refuerzo para materiales compuestos. Las características que se desean en estos refuerzos de fibras son alta resistencia, alta rigidez y baja densidad relativa. Las fibras de vidrio son las que más comúnmente se usan cuando se requiere un buen desempeño mecánico en los materiales compuestos, debido a su alta resistencia a la tracción y bajo costo. A continuación se muestra la figura 9 que muestra la curva de esfuerzo – deformación de fibras típicas usadas como refuerzos. 14 Figura A2.5. Curva esfuerzo – deformación de fibras típicas 1.5.2. Matrices Los tipos de matrices que se usan en materiales compuestos: polimérica, metálica, cerámica, carbono. Las matrices que comúnmente se usan son las poliméricas, las cuales pueden ser termoestables (epoxi, poliamida, poliéster) o termoplásticas (polisulfona, poliéster-éter cetona). Las matrices de poliéster poseen una buena estabilidad en sus propiedades mecánicas, químicas y eléctricas, ofrecen además una buena resistencia a la corrosión. Las matrices de poliéster más usadas son del tipo ortoftálico y la isoftálica. Las matrices de poliéster ortoftálicas son las matrices de uso general, es decir los elementos que están incluidos en su polimerización son ácidos saturados de bajo costo, es por eso que presentan bajas propiedades mecánicas y químicas comparadas con la matriz de poliéster isoftálica. La matriz de poliéster isoftálica está basada en el ácido isoftálico y en el compuesto orgánico anhídrido maleico. La incorporación del ácido isoftálico crea una resina de alto peso molecular, dando buenas propiedades químicas y mecánicas. En la siguiente figura A2.6 se muestran las propiedades mecánicas que presentan estos dos tipos de poliéster. Figura A2.6. Propiedades mecánicas de los dos tipos de poliéster. 1.6. Propiedades de materiales compuestos típicos 15 Los materiales compuestos que incorporan los constituyentes discutidos anteriormente muestran una amplia gama de características. El comportamiento de los materiales compuestos unidireccionales en la dirección de la fibra, en especial la rigidez, es usualmente dominada por las propiedades de las fibras. Las curvas de esfuerzo-deformación típicas de materiales compuestos unidireccionales sometidos a tracción en la dirección de la fibra se muestran a continuación, comparadas con la de aluminio. En estas curvas se puede observar que mientras la rigidez aumenta, la deformación a la rotura disminuye. Además que todos los materiales compuestos mostrados presentan un comportamiento casi lineal además de tener un esfuerzo a la rotura mayor al aluminio. La siguiente figura A2.7 muestra los materiales sometidos a ensayos de tensión. Figura A2.7. Resultados de materiales sometidos a ensayos de tensión. El comportamiento de materiales compuestos unidireccionales en la dirección transversal a las fibras, en especial la resistencia, está dominado por la matriz. Las curvas de esfuerzo-deformación de materiales compuestos unidireccionales típicos sometidos a tracción en la dirección transversal a las fibras se muestran a continuación. Todos de esos materiales exhiben un comportamiento casi lineal, con bajas deformaciones a la rotura y bajas resistencias, también se observa que los materiales que poseen una matriz epoxi, tienen resistencia transversal casi iguales. 16 Figura A2.8. Curva de esfuerzo deformación de materiales compuestos unidireccionales sometidos a tensión. Se muestra en la siguiente figura A2.9 una lista de propiedades que presentan algunos de los materiales compuestos típicos. Las propiedades mostradas en la figura A2.9 mencionada anteriormente fueron obtenidas a una temperatura ambiental (24 °C) y condición de humedad cero. Figura A2.9. Propiedades de algunos materiales compuestos. 17 ANEXO 3 DESARROLLO DEL ALGORITMO MATEMÁTICO DE LA TEORÍA CLÁSICA DE LAMINADO 1. Relación esfuerzo-deformación en una material ortotrópico. Como se mostró en la sección anterior un material compuesto reforzado con fibra de vidrio unidireccional posee 3 ejes principales (1, 2 y 3), el eje 1 es el eje donde las fibras están orientadas, el eje 2 es el eje transversal a las fibras, donde generalmente las propiedades del compuesto están determinadas por la matriz, y el eje 3 que es el eje normal a los ejes 1 y 2 (ejes en el plano). Esta definición hace que el material se denomine como ortotrópico, como se definió anteriormente un material ortotrópico es aquel que posee tres ejes perpendiculares entre sí, y presenta propiedades constantes definidas a lo largo de estos, además estos ejes son los ejes principales del material. El estado de esfuerzos principales para un volumen definido se muestra en la Figura A3.1. Figura A3.1. Estado de esfuerzos principales. Para el estado de esfuerzos anterior se cumple lo siguiente: 𝜎𝑖𝑗 = 𝜎𝑗𝑖 𝜀𝑖𝑗 = 𝜀𝑗𝑖 (Ec.3.1) 𝜎 11 = 𝜎1 ; 𝜎22 = 𝜎2 ; 𝜎33 = 𝜎3 ; 𝜎23 = 𝜏4 ; 𝜎31 = 𝜏5 ; 𝜎12 = 𝜏6 (Ec.3.2) Por lo tanto las matrices de elasticidad que caracterizan el comportamiento de este material ortotrópico son las siguientes: 18 [ 𝜎1 𝜎2 𝜎3 𝜏4 𝜏5 𝜏6] = [ 𝐶11 𝐶12 𝐶13 0 0 0 𝐶12 𝐶22 𝐶23 0 0 0 𝐶13 𝐶23 𝐶33 0 0 0 0 0 0 𝐶44 0 0 0 0 0 0 𝐶55 0 0 0 0 0 0 𝐶66] . [ 𝜀1 𝜀2 𝜀3 𝛾4 𝛾5 𝛾6] (Ec.3.3.) [ 𝜀1 𝜀2 𝜀3 𝛾4 𝛾5 𝛾6] = [ 𝑆11 𝑆12 𝑆13 0 0 0 𝑆12 𝑆22 𝑆23 0 0 0 𝑆13 𝑆23 𝑆33 0 0 0 0 0 0 𝑆44 0 0 0 0 0 0 𝑆55 0 0 0 0 0 0 𝑆66] . [ 𝜎1 𝜎2 𝜎3 𝜏4 𝜏5 𝜏6] (Ec.3.4.) En las ecuaciones (2.3) y (2.4), las constantes 𝐶𝑖𝑗 y 𝑆𝑖𝑗 , son parámetros elásticos que dependen de las constantes de ingeniería (módulo de elasticidad, coeficiente de poisson, módulo cortante). Para cada una de las matrices anteriores se cumple lo siguiente: 𝐶𝑖𝑗 = 𝐶𝑗𝑖 𝑆𝑖𝑗 = 𝑆𝑗𝑖 (Ec.3.5) Se muestra en las ecuaciones anteriores que un material ortotrópico puede ser caracterizado por nueve constantes elásticas independientes. La matriz [𝑆𝑖𝑗] es inversa de [𝐶𝑖𝑗]. Se pueden hacer las siguientes observaciones de las relaciones de esfuerzo-deformación que se muestran en 2.3 y 2.4. 1. No existe interacción alguna entre los esfuerzos normales 𝜎1, 𝜎2, 𝜎3 y las deformaciones cortantes 𝛾1, 𝛾2 , 𝛾3 , es decir esfuerzos normales actuando en las direcciones principales del material solo generan deformaciones normales. 2. No existe interacción entre los esfuerzos cortantes 𝜏1, 𝜏2, 𝜏3 , y las deformaciones normales 𝜀1, 𝜖2, 𝜀3 es decir esfuerzos cortantes actuando en los planos principales del material solo generan deformaciones cortantes. 3. No existe interacción entre esfuerzos cortantes y deformaciones cortantes en diferentes planos, es decir un esfuerzo cortante actuando en un plano principal produce una deformación cortante en ese plano. 2.1.1 Material ortotrópico bajo esfuerzos en el plano En la mayoría de aplicaciones estructurales los materiales compuestos son usados en la forma de laminados delgados, sometidas a cargas en el plano del laminado. Entonces el material compuesto (o laminado) puede ser considerado a estar bajo una condición de tensión plana, con todos los componentes fuera del plano del laminado (dirección en el eje 3) nulos. 19 𝜎3 = 0 𝜏23 = 𝜏4 = 0 𝜏13 = 𝜏5 = 0 (Ec.3.6) La matriz de esfuerzos-deformaciones (Ecuación 2.3) se reduce a la siguiente forma: [ 𝜎1 𝜎2 0 0 0 𝜏6] = [ 𝐶11 𝐶12 𝐶13 0 0 0 𝐶12 𝐶22 𝐶23 0 0 0 𝐶13 𝐶23 𝐶33 0 0 0 0 0 0 𝐶44 0 0 0 0 0 0 𝐶55 0 0 0 0 0 0 𝐶66] . [ 𝜀1 𝜀2 𝜀3 𝛾4 𝛾5 𝛾6] (Ec.3.7) Resolviendo la matriz se obtiene lo siguiente: 𝜎1 = 𝐶11𝜀1 + 𝐶12𝜀2 + 𝐶13𝜀3 𝜎2 = 𝐶12𝜀1 + 𝐶22𝜀2 + 𝐶23𝜀3 0 = 𝐶13𝜀1 + 𝐶23𝜀2 + 𝐶33𝜀3 𝛾4 = 𝛾5 = 0 𝜏6 = 𝐶66𝛾6 (Ec.3.8) Eliminando la deformación 𝜀3 de la Ecuación (2.9), se tiene lo siguiente: 𝜎1 = (𝐶11 − 𝐶13𝐶13 𝐶33 ) 𝜀1 + (𝐶12 − 𝐶13𝐶23 𝐶33 )𝜀2 𝜎1 = 𝑄11𝜀1 + 𝑄12𝜀2 𝜎2 = (𝐶12 − 𝐶23𝐶13 𝐶33 ) 𝜀1 + (𝐶22 − 𝐶23𝐶23 𝐶33 )𝜀2 𝜎2 = 𝑄12𝜀1 + 𝑄22𝜀2 𝜏6 = 𝐶66𝛾6 = 𝑄66𝛾6 (Ec.3.9) Las ecuaciones (2.10) se pueden mostrar resumidas en la siguiente matriz: [ 𝜎1 𝜎2 𝜏6 ] = [ 𝑄11 𝑄12 0 𝑄12 𝑄22 0 0 0 𝑄66 ] . [ 𝜀1 𝜀2 𝛾6 ] (Ec.3.10) La matriz inversa de la ecuación (2.11) es la siguiente: [ 𝜀1 𝜀2 𝛾6 ] = [ 𝑆11 𝑆12 0 𝑆12 𝑆22 0 0 0 𝑆66 ] . [ 𝜎1 𝜎2 𝜎3 ] (Ec.3.11) 20 Como se puede observar en las matrices anteriores, las relaciones de esfuerzo-deformación bajo cargas en el plano pueden expresarse en cuatro parámetros elásticos independientes 𝑄11, 𝑄12, 𝑄22 , 𝑄66 ó 𝑆11, 𝑆12, 𝑆22 , 𝑆66 2. Algoritmo matemático del material Para poder usar el algoritmo matemático se siguen los siguientes pasos. 1. Introducir las propiedades básicas de la lámina ( 𝑬𝟏 ; 𝑬𝟐 ; 𝑮𝟏𝟐 ; 𝝂𝟏𝟐) Todo laminado analizado en el presente trabajo está compuesto por un número determinado de láminas, cada una teniendo una determinada orientación, sin embargo cada una de estas láminas presenta los mismo valores de estos parámetros en sus direcciones principales, se halla estas propiedades básicas de la lámina referente a sus ejes principales mediante ensayos, denominados por los subíndices 1 y 2, que significan respectivamente, el eje longitudinal de las fibras y el eje transversal a estas. 2. Hallar la matriz de rigidez [𝑸]𝟏𝟐 de cada lámina referido a sus ejes principales (1,2) usando la ecuación (A.13). Para este trabajo todas láminas son del mismo material (matriz de poliéster isoftálica reforzada con fibra de vidrio), y debido a esto todas las láminas presentan la misma matriz de rigidez referida a sus ejes principales, sin embargo para posteriores cálculos se necesita tener la matriz referida a los ejes del sistema de referencia del laminado (x , y), y es acá donde se verá el efecto de la orientación de las láminas. 3. Habiendo seleccionado ya los ejes de referencia del laminado (x , y) , e indicado las cotas superiores (𝒉𝒌) e inferiores (𝒉𝒌−𝟏 ). de cada lámina que conforman el laminado, se procede a calcular las matrices de rigidez del laminado ([𝑨] ; [𝑩] ; [𝑫] , usando la ecuación (A.33). 4. Se procede a calcular las matrices “inversas” del laminado ([𝒂] ; [𝒃] ; [𝒄] ; [𝒅]) según la ecuación (A.38). 5. Se calcula las deformaciones del plano de referencia [𝜺𝟎] 𝒙,𝒚 y curvaturas [𝒌]𝒙,𝒚 , previamente habiendo ingresado la carga mecánica [𝑵]𝒙,𝒚 y momentos [𝑴]𝒙,𝒚 a los que está sometido el laminado, según las ecuaciones (A.36) y (A.37). 21 6. Se calcula las deformaciones de cada lámina “k”( [𝜺𝒌] 𝒙,𝒚 ) referido a los ejes de referencia del laminado (x, y), según la ecuación (A.21). Para analizar la deformación real de cada lámina, se introduce la coordenada “z” de interés para su análisis. Para laminados simétricos bajo carga en el plano (los cuales se analizan en el presente trabajo), se usa la coordenada denominada . Para los casos de laminados sometidos a flexión se analizará los esfuerzos máximos en cada lámina, por lo que la coordenada , será medida desde el plano medio del laminado hasta la cota superior de cada lámina. 29 7. Si los ejes principales de una lámina “k” no coinciden con los ejes de referencia del laminado (x, y), las deformaciones halladas en el paso anterior no serán las principales ([𝜺𝒌] 𝒙,𝒚 ), por lo cual se deberá usar la siguiente ecuación (A.40) para poder hallarlas. 8. Calcular los esfuerzos principales (1, 2) en cada lámina [𝝈𝒌] 𝟏,𝟐 , según la ecuación (2.11). 9. Aplicar una teoría de falla (s). Existen muchas teorías de falla para materiales compuestos (laminados), pero por practicidad se analiza la teoría del esfuerzo máximo y deformación máxima, los cuales se entienden que si algún valor de los esfuerzos o deformaciones principales halladas de cada lámina supera el valor máximo permisible, el material falla en dicha lámina. 3. Consideraciones básicas del modelo matemático del material Relaciones entre los parámetros elásticos y las constantes de ingeniería Los parámetros elásticos 𝐶𝑖𝑗 y 𝐶𝑗𝑖 , que se explicaron líneas atras, se deben expresar en términos de las constantes de ingeniería tales como el módulo de elasticidad, coeficiente de poisson, etc. Estas relaciones son obtenidas realizando experimentos elementales imaginarios que se describen a continuación. Si se sometiera al material a un esfuerzo 𝜎1 , se cumple: 𝜀1 = 𝑆11𝜎1 𝜀2 = 𝑆12𝜎1 𝛾6 = 0 (A.1) De consideraciones de ingeniería la ecuación (A.1) queda de la siguiente forma: 𝜀1 = 𝜎1 𝐸1 𝜀2 = − 𝜐12 𝐸1 𝜎1 22 𝛾6 = 0 (A.2) Cabe señalar que los subíndices en el coeficiente de poisson denotan la dirección del esfuerzo y deformación, respectivamente. De las ecuaciones (A.1) y (A.2) se obtienen las siguientes relaciones: 𝑆11 = 1 𝐸1 ; 𝑆12 = − 𝜐12 𝐸1 (A.3) De manera similar, si el material estuviese sometido a un esfuerzo 𝜎2 , tenemos: 𝜀1 = 𝑆12𝜎2 = − 𝜐21 𝐸2 𝜎2 𝜀2 = 𝑆22𝜎2 = 𝜎2 𝐸2 𝛾6 = 0 (A.4) De donde obtenemos las siguientes relaciones: 𝑆12 = − 𝜐21 𝐸1 ; 𝑆22 = 1 𝐸2 (A.5) Si sometiéramos al material a esfuerzo cortante puro 𝜏6 , tenemos: 𝜀1 = 𝜀2 = 0 𝛾6 = 𝑆66𝜏6 = 𝜏6 𝐺12 (A.6) De donde se obtiene: 𝑆66 = 1 𝐺12 (A.7) La matriz de los componentes hallados en (A.3), (A.5) y (A.7) queda de la siguiente forma: [𝑆𝑖𝑗] = [ 1 𝐸1 − 𝜐21 𝐸2 0 − 𝜐21 𝐸2 1 𝐸2 0 0 0 1 𝐺2] (A.8) La matriz (2.12) quedaría de la siguiente manera: 23 [ 𝜀1 𝜀2 𝛾6 ] = [ 1 𝐸1 − 𝜐21 𝐸2 0 − 𝜐21 𝐸2 1 𝐸2 0 0 0 1 𝐺2] . [ 𝜎1 𝜎2 𝜏6 ] (A.9) Se puede observar que se cumple lo siguiente: 𝜐12 𝐸1 = 𝜐21 𝐸2 (A.10) Para expresar la matriz (2.11) en términos de las constantes de ingeniería, se debe hallar la inversa de la matriz [𝑆𝑖𝑗], y el resultado sería el siguiente: [𝑆𝑖𝑗] −1 = [ 𝐸1 1 − 𝜐12. 𝜐21 𝜐21. 𝐸1 1 − 𝜐12. 𝜐21 0 𝜐21. 𝐸1 1 − 𝜐12. 𝜐21 𝐸2 1 − 𝜐12. 𝜐21 0 0 0 𝐺12] (A.11) [ 𝜎1 𝜎2 𝜏6 ] = [ 𝐸1 1 − 𝜐12. 𝜐21 𝜐21. 𝐸1 1 − 𝜐12. 𝜐21 0 𝜐21. 𝐸1 1 − 𝜐12. 𝜐21 𝐸2 1 − 𝜐12. 𝜐21 0 0 0 𝐺12] . [ 𝜀1 𝜀2 𝛾6 ] (A.12) Dónde: 𝑸𝟏𝟏 = 𝐸1 1 − 𝜐12. 𝜐21 𝑸𝟏𝟐 = 𝜐21. 𝐸1 1 − 𝜐12. 𝜐21 = 𝜐12. 𝐸2 1 − 𝜐12. 𝜐21 𝑸𝟐𝟐 = 𝐸2 1 − 𝜐12. 𝜐21 𝑸𝟔𝟔 = 𝐺12 (A.13) 3.1. Consideraciones básicas El comportamiento global de un material compuesto laminado está en función de las propiedades de cada lámina que lo conforman, así como la distribución de estas láminas dentro del laminado. El modelo matemático planteado a continuación se basa en la teoría clásica de laminación, la cual predice el comportamiento del laminado bajo las siguientes consideraciones: 24  Cada lámina del laminado es quasi-homogénea y ortotrópica.  El laminado es delgado con sus dimensiones laterales más grandes que su espesor, y está cargado solamente en su plano. El laminado y sus láminas están en un estado de tensión plana (𝜎𝑧 = 𝜏𝑥𝑧 = 𝜏𝑦𝑧 = 0).  Todos los desplazamientos son pequeños comparados con el espesor del laminado ( |𝑢| , |𝑣| , |𝑤| ≪ ℎ) .  Los desplazamientos son continuos a través del laminado.  Los desplazamientos en el plano varían linealmente a través del espesor del laminado. Como se verá en este capítulo los desplazamientos 𝑢 , 𝑣 en las direcciones x- e y- respectivamente son lineales y están en función de z.  Las deformaciones cortantes transversales 𝜸𝑥𝑧 y 𝜸𝑦𝑧 son despreciables.  Relaciones de esfuerzo-deformación son lineales.  Las distancias normales del medio de la superficie permanecen constantes, es decir la deformación transversal normal 𝜺𝒛 es despreciable (comparado con las otras deformaciones en el plano 𝜺𝒙 y 𝜺𝒚 . 3.2. Relaciones de desplazamiento-deformación Figura A3.2. Sección del laminado antes (ABCD) y después (A´B´C´D´) de la deformación [1] La Figura 15 muestra una sección del laminado normal al eje y- antes y después de sufrir una deformación. El plano xy es equidistante de la superficie superior e inferior del laminado, y se le denomina plano de referencia. Los desplazamientos del plano de referencia 𝑢0 y 𝑣0 en los ejes x- e y- respectivamente y el desplazamiento 𝑤 fuera del plano en la dirección del eje z, son solamente funciones de x e y. 𝑢0 = 𝑢0 (𝑥, 𝑦) 𝑣0 = 𝑣0 (𝑥, 𝑦) 25 𝑤 = 𝑓 (𝑥, 𝑦) (A.14) Las rotaciones de los ejes x- e y- son: 𝛼𝑥 = 𝜕𝑤 𝜕𝑥 𝛼𝑦 = 𝜕𝑤 𝜕𝑦 (A.15) Los componentes de desplazamiento en el plano del punto B de coordenada 𝑧𝑏 como se muestra en la Figura A.1 son: 𝑢𝑏 = 𝑢0 − 𝛼𝑥 𝑧𝑏 𝑣𝑏 = 𝑣0 − 𝛼𝑦 𝑧𝑏 (A.16) Reemplazando las ecuaciones (A.15) en (A.16), y poniendo la ecuación de forma general, se tiene: 𝑢 = 𝑢0 − 𝑧 𝜕𝑤 𝜕𝑥 𝑣 = 𝑣0 − 𝑧 𝜕𝑤 𝜕𝑦 (A.17) Donde 𝑧 es la coordenada variable de un punto general de la sección transversal. Para pequeños desplazamientos, se obtienen las siguientes relaciones de desplazamiento-deformación de la clásica teoría de elasticidad y derivando las ecuaciones (A.17). 𝜀𝑥 = 𝜕𝑢 𝜕𝑥 = 𝜕𝑢0 𝜕𝑥 − 𝑧 𝜕2𝑤 𝜕𝑥2 𝜀𝑦 = 𝜕𝑣 𝜕𝑦 = 𝜕𝑣0 𝜕𝑦 − 𝑧 𝜕2𝑤 𝜕𝑦2 𝛾𝑥𝑦 = 𝛾𝑠 = 𝜕𝑢 𝜕𝑦 + 𝜕𝑣 𝜕𝑥 = 𝜕𝑢0 𝜕𝑦 + 𝜕𝑣0 𝜕𝑥 − 2𝑧 𝜕2𝑤 𝜕𝑥𝜕𝑦 𝜀𝑧 = 𝛾𝑥𝑧 = 𝛾𝑦𝑧 = 0 (A.18) Las componentes de deformación en el plano de referencia se expresan de la siguiente manera: 𝜀𝑥 0 = 𝜕𝑢0 𝜕𝑥 𝜀𝑦 0 = 𝜕𝑣0 𝜕𝑦 𝛾𝑥𝑦 0 = 𝛾𝑠 0 = 𝜕𝑢0 𝜕𝑦 + 𝜕𝑣0 𝜕𝑥 26 (A.19) Las curvaturas que sufre el laminado se puede expresar también de la siguiente manera: ĸ𝑥 = − 𝜕2𝑤 𝜕𝑥2 ĸ𝑦 = − 𝜕2𝑤 𝜕𝑦2 ĸ𝑥𝑦 = ĸ𝑠 = − 2𝜕2𝑤 𝜕𝑥𝜕𝑦 (A.20) Entonces, reemplazando las ecuaciones (A.19) y (A.20) en la ecuación (A.18), se obtiene lo siguiente: [ 𝜀𝑥 𝜀𝑦 𝛾𝑠 ] = [ 𝜀𝑥 0 𝜀𝑦 0 𝛾𝑠 0 ] + 𝑧. [ ĸ𝑥 ĸ𝑦 ĸ𝑠 ] (A.21) De esta manera podemos relacionar las deformaciones para cualquier punto en el laminado en función de las deformaciones del plano de referencia y las curvaturas del laminado. 3.3. Relaciones de esfuerzo-deformación de la lámina dentro del laminado Se considera una lámina individual “k” dentro de un laminado multidireccional cuyo plano medio está a una distancia 𝑧𝑘 del plano de referencia del laminado (Figura A3.3). Figura A3.3. Lámina “k” dentro del laminado [1] Las relaciones de esfuerzo-deformación para esta lámina referida a sus ejes principales se muestran en las siguientes matrices: [ 𝜎1 𝜎2 𝜏6 ] 𝑘 = [ 𝑄11 𝑄12 0 𝑄21 𝑄22 0 0 0 𝑄66 ] 𝑘 . [ 𝜀1 𝜀2 𝛾6 ] 𝑘 (A.22) La matriz mostrada en (A.22) se transforma al sistema de ejes coordenados del laminado: 27 [ 𝜎𝑥 𝜎𝑦 𝜏𝑠 ] 𝑘 = [ 𝑄𝑥𝑥 𝑄𝑥𝑦 𝑄𝑥𝑠 𝑄𝑦𝑥 𝑄𝑦𝑦 𝑄𝑦𝑠 𝑄𝑠𝑥 𝑄𝑠𝑦 𝑄𝑠𝑠 ] 𝑘 . [ 𝜀𝑥 𝜀𝑦 𝛾𝑠 ] 𝑘 (A.23) Sustituyendo las expresiones de deformaciones de (A.21) en (A.23), obtenemos: [ 𝜎𝑥 𝜎𝑦 𝜏𝑠 ] 𝑘 = [ 𝑄𝑥𝑥 𝑄𝑥𝑦 𝑄𝑥𝑠 𝑄𝑦𝑥 𝑄𝑦𝑦 𝑄𝑦𝑠 𝑄𝑠𝑥 𝑄𝑠𝑦 𝑄𝑠𝑠 ] 𝑘 . [ 𝜀𝑥 0 𝜀𝑦 0 𝛾𝑠 0 ] + 𝑧. [ 𝑄𝑥𝑥 𝑄𝑥𝑦 𝑄𝑥𝑠 𝑄𝑦𝑥 𝑄𝑦𝑦 𝑄𝑦𝑠 𝑄𝑠𝑥 𝑄𝑠𝑦 𝑄𝑠𝑠 ] 𝑘 . [ ĸ𝑥 ĸ𝑦 ĸ𝑠 ] (A.24) De las ecuaciones (A.21) y (A.24), se puede observar que mientras las deformaciones varían linealmente a través del espesor del laminado, los esfuerzos no son lineales. Esto es debido a la variación discontinua de la matriz de rigidez transformada [𝑄]𝑥,𝑦 de lámina a lámina dentro del laminado, lo que haría también que los esfuerzos varían discontinuamente de lámina a lámina. Esto se ilustra tomando como ejemplo un laminado compuesto por 4 láminas como se muestra en la Figura A3.4, sometido bajo esfuerzo uniaxial en la dirección –x. Se observa que para una determinada variación lineal de la deformación a través del espesor, la cual puede ser debido a una carga axial o carga de flexión, la variación del módulo 𝐸𝑥 de lámina a lámina puede causar la variación discontinua del esfuerzo. Figura A3.4. Ilustración de la variación lineal de las deformaciones y la variación discontinua del esfuerzo en un laminado multidireccional [1] Los esfuerzos promedio en cada lámina quedan en función de las deformaciones del plano de referencia del laminado [𝜀0]𝒙,𝒚 , las curvaturas [ĸ]𝒙,𝒚 del laminado, la locación del plano medio de la lámina 𝑧𝑘 respecto al plano de referencia del laminado y la matriz de rigidez transformada [𝑄]𝑥,𝑦 de la lámina como se muestra en la ecuación (A.24). 3.4. fuerza y momentos resultantes de un laminado Como se dijo anteriormente debido a la variación discontinua de los esfuerzos de lámina a lámina, es más conveniente buscar expresiones relacionando fuerza y momentos con la deformación que 28 sufre el laminado. Los esfuerzos que actúan en una lámina “k” del laminado dados por la ecuación (A.24) pueden ser reemplazados por fuerzas y momentos resultantes como se muestra en la Figura A3.5 y según las siguientes ecuaciones: Figura A3.5. Lámina “k” sometida a fuerzas y momentos resultantes [1] 𝑁𝑥 𝑘 = ∫ 𝜎𝑥 𝑡/2 −𝑡/2 𝜕𝑧 𝑁𝑦 𝑘 = ∫ 𝜎𝑦 𝑡/2 −𝑡/2 𝜕𝑧 𝑁𝑥𝑦 𝑘 = 𝑁𝑠 𝑘 = ∫ 𝜏𝑠 𝑡/2 −𝑡/2 𝜕𝑧 (A.25) 𝑀𝑥 𝑘 = ∫ 𝜎𝑥 𝑡/2 −𝑡/2 𝑧 𝜕𝑧 𝑀𝑦 𝑘 = ∫ 𝜎𝑦 𝑡/2 −𝑡/2 𝑧 𝜕𝑧 𝑀𝑥𝑦 𝑘 = 𝑀𝑠 𝑘 = ∫ 𝜏𝑠 𝑡/2 −𝑡/2 𝑧 𝜕𝑧 (A.26) Dónde: 𝑧 = la coordenada variable de un punto en la sección transversal. 𝑡 = Espesor de la lámina. 𝑁𝑥 𝑘 , 𝑁𝑦 𝑘 = fuerzas normales por unidad de longitud. 𝑁𝑠 𝑘 = fuerza cortante por unidad de longitud. 𝑀𝑥 𝑘 , 𝑀𝑦 𝑘 = momentos de flexión por unidad de longitud. 𝑀𝑠 𝑘 = momento de torsión por unidad de longitud. 29 Ahora, para un laminado de varias láminas, la fuerza total y los momentos resultantes son obtenidos sumando los efectos de todas las láminas. Entonces, para un laminado compuesto por n-láminas como se muestra en la Figura 19, la fuerza y los momentos resultantes son obtenidos según: [ 𝑁𝑥 𝑁𝑦 𝑁𝑠 ] = ∑ ∫ [ 𝜎𝑥 𝜎𝑦 𝜏𝑠 ] 𝑘 ℎ𝑘 ℎ𝑘−1 𝜕𝑧 𝑛 𝑘=1 (A.27) [ 𝑀𝑥 𝑀𝑦 𝑀𝑠 ] = ∑ ∫ [ 𝜎𝑥 𝜎𝑦 𝜏𝑠 ] 𝑘 ℎ𝑘 ℎ𝑘−1 𝑧𝜕𝑧 𝑛 𝑘=1 (A.28) Donde ℎ𝑘 y ℎ𝑘−1 son las coordenadas 𝑧 de la superficie superior e inferior respectivamente de la lámina 𝑘 . Figura A3.6. Laminado multidireccional, indicando las coordenadas de cada lámina [1] 3.5. Relaciones generales de carga-deformación Sustituyendo la ecuación (A.24) de esfuerzos para una lámina en las ecuaciones (A.27) y (A.28), se obtiene lo siguiente: [ 𝑁𝑥 𝑁𝑦 𝑁𝑠 ] = ∑ 𝑛 𝑘=1 {[ 𝑄𝑥𝑥 𝑄𝑥𝑦 𝑄𝑥𝑠 𝑄𝑦𝑥 𝑄𝑦𝑦 𝑄𝑦𝑠 𝑄𝑠𝑥 𝑄𝑠𝑦 𝑄𝑠𝑠 ] 𝑘 . [ 𝜀𝑥 0 𝜀𝑦 0 𝛾𝑠 0 ]∫ 𝜕𝑧 ℎ𝑘 ℎ𝑘−1 + [ 𝑄𝑥𝑥 𝑄𝑥𝑦 𝑄𝑥𝑠 𝑄𝑦𝑥 𝑄𝑦𝑦 𝑄𝑦𝑠 𝑄𝑠𝑥 𝑄𝑠𝑦 𝑄𝑠𝑠 ] 𝑘 . [ ĸ𝑥 ĸ𝑦 ĸ𝑠 ]∫ 𝑧 ℎ𝑘 ℎ𝑘−1 𝜕𝑧 } (A.29) [ 𝑀𝑥 𝑀𝑦 𝑀𝑠 ] = ∑ 𝑛 𝑘=1 {[ 𝑄𝑥𝑥 𝑄𝑥𝑦 𝑄𝑥𝑠 𝑄𝑦𝑥 𝑄𝑦𝑦 𝑄𝑦𝑠 𝑄𝑠𝑥 𝑄𝑠𝑦 𝑄𝑠𝑠 ] 𝑘 . [ 𝜀𝑥 0 𝜀𝑦 0 𝛾𝑠 0 ]∫ 𝑧𝜕𝑧 ℎ𝑘 ℎ𝑘−1 + [ 𝑄𝑥𝑥 𝑄𝑥𝑦 𝑄𝑥𝑠 𝑄𝑦𝑥 𝑄𝑦𝑦 𝑄𝑦𝑠 𝑄𝑠𝑥 𝑄𝑠𝑦 𝑄𝑠𝑠 ] 𝑘 . [ ĸ𝑥 ĸ𝑦 ĸ𝑠 ]∫ 𝑧2 ℎ𝑘 ℎ𝑘−1 𝜕𝑧 } (A.30) En las expresiones anteriores la matriz de rigidez [𝑄]𝑘 𝑥,𝑦 , las deformaciones del plano de referencia [𝜀0]𝑥,𝑦 , y las curvaturas [ĸ]𝑥,𝑦 son sacadas de la integral, debido a que estas no dependen de 𝑧. De 30 estas tres cantidades solo la matriz de rigidez [𝑄]𝑘 𝑥,𝑦 es única para cada lámina, mientras que las deformaciones del plano de referencia [𝜀0]𝑥,𝑦 , y las curvaturas [ĸ]𝑥,𝑦 dependen del laminado entero y son los mismos para todas las láminas. Es por eso que [𝜀0]𝑥,𝑦 , y [ĸ]𝑥,𝑦 pueden ser factorizados afuera de la sumatoria como se muestra a continuación: [𝑁]𝑥,𝑦 = [∑[𝑄] 𝑘 𝑥,𝑦 ∫ 𝜕𝑧 ℎ𝑘 ℎ𝑘−1 𝑛 𝑘=1 ] [𝜀0]𝑥,𝑦 + [∑[𝑄] 𝑘 𝑥,𝑦 ∫ 𝑧𝜕𝑧 ℎ𝑘 ℎ𝑘−1 𝑛 𝑘=1 ] [ĸ]𝑥,𝑦 [𝑁]𝑥,𝑦 = [∑[𝑄] 𝑘 𝑥,𝑦 (ℎ𝑘 − ℎ𝑘−1) 𝑛 𝑘=1 ] [𝜀0]𝑥,𝑦 + [ 1 2 ∑[𝑄]𝑘 𝑥,𝑦 (ℎ𝑘 2 − ℎ𝑘−1 2) 𝑛 𝑘=1 ] [ĸ]𝑥,𝑦 [𝑁]𝑥,𝑦 = [𝐴]𝑥,𝑦 [𝜀 0]𝑥,𝑦 + [𝐵]𝑥,𝑦 [ĸ]𝑥,𝑦 (A.31) [𝑀]𝑥,𝑦 = [∑[𝑄] 𝑘 𝑥,𝑦 ∫ 𝑧𝜕𝑧 ℎ𝑘 ℎ𝑘−1 𝑛 𝑘=1 ] [𝜀0]𝑥,𝑦 + [∑[𝑄] 𝑘 𝑥,𝑦 ∫ 𝑧2𝜕𝑧 ℎ𝑘 ℎ𝑘−1 𝑛 𝑘=1 ] [ĸ]𝑥,𝑦 [𝑀]𝑥,𝑦 = [ 1 2 ∑[𝑄]𝑘 𝑥,𝑦 (ℎ𝑘 2 − ℎ𝑘−1 2) 𝑛 𝑘=1 ] [𝜀0]𝑥,𝑦 + [ 1 3 ∑[𝑄]𝑘 𝑥,𝑦 (ℎ𝑘 3 − ℎ𝑘−1 3) 𝑛 𝑘=1 ] [ĸ]𝑥,𝑦 [𝑀]𝑥,𝑦 = [𝐵]𝑥,𝑦 [𝜀 0]𝑥,𝑦 + [𝐷]𝑥,𝑦 [ĸ]𝑥,𝑦 (A.32) Dónde: 𝐴𝑖𝑗 = ∑ 𝑄𝑖𝑗 𝑘 (ℎ𝑘 − ℎ𝑘−1) 𝑛 𝑘=1 𝐵𝑖𝑗 = 1 2 ∑ 𝑄𝑖𝑗 𝑘 (ℎ𝑘 2 − ℎ𝑘−1 2) 𝑛 𝑘=1 𝐷𝑖𝑗 = 1 3 ∑ 𝑄𝑖𝑗 𝑘 (ℎ𝑘 3 − ℎ𝑘−1 3) 𝑛 𝑘=1 𝑖, 𝑗 = 𝑥, 𝑦, 𝑠. (A.33) Entonces, en forma completa las ecuaciones (A.31) y (A.32) quedan de la siguiente manera: [ 𝑁𝑥 𝑁𝑦 𝑁𝑠 ] = [ 𝐴𝑥𝑥 𝐴𝑥𝑦 𝐴𝑥𝑠 𝐴𝑦𝑥 𝐴𝑦𝑦 𝐴𝑦𝑠 𝐴𝑠𝑥 𝐴𝑠𝑦 𝐴𝑠𝑠 ] . [ 𝜀𝑥 0 𝜀𝑦 0 𝛾𝑠 0 ] + [ 𝐵𝑥𝑥 𝐵𝑥𝑦 𝐵𝑥𝑠 𝐵𝑦𝑥 𝐵𝑦𝑦 𝐵𝑦𝑠 𝐵𝑠𝑥 𝐵𝑠𝑦 𝐵𝑠𝑠 ] . [ ĸ𝑥 ĸ𝑦 ĸ𝑠 ] (A.34) [ 𝑀𝑥 𝑀𝑦 𝑀𝑠 ] = [ 𝐵𝑥𝑥 𝐵𝑥𝑦 𝐵𝑥𝑠 𝐵𝑦𝑥 𝐵𝑦𝑦 𝐵𝑦𝑠 𝐵𝑠𝑥 𝐵𝑠𝑦 𝐵𝑠𝑠 ] . [ 𝜀𝑥 0 𝜀𝑦 0 𝛾𝑠 0 ] + [ 𝐷𝑥𝑥 𝐷𝑥𝑦 𝐷𝑥𝑠 𝐷𝑦𝑥 𝐷𝑦𝑦 𝐷𝑦𝑠 𝐷𝑠𝑥 𝐷𝑠𝑦 𝐷𝑠𝑠 ] . [ ĸ𝑥 ĸ𝑦 ĸ𝑠 ] (A.35) 31 Todas las matrices anteriores son simétricas, por lo tanto: 𝐴𝑖𝑗 = 𝐴𝑗𝑖 𝐵𝑖𝑗 = 𝐵𝑗𝑖 𝐷𝑖𝑗 = 𝐷𝑗𝑖 𝑖, 𝑗 = 𝑥, 𝑦, 𝑠. Las relaciones anteriormente descritas son expresadas en términos de 3 matrices de rigidez del laminado [𝐴] , [𝐵] , y [𝐷], las cuales están en función de la geometría, propiedades del material y la secuencia de las láminas como se definió en la ecuación (A.33). Estos son los parámetros elásticos del laminado multidireccional con el siguiente significado: 𝐴𝑖𝑗: Es la matriz de rigidez la cual relaciona cargas en el plano con deformaciones en el plano. 𝐵𝑖𝑗: Es la matriz de rigidez la cual relaciona las cargas en el plano con las curvaturas que sufre el laminado, así como los momentos con las deformaciones en el plano. Si 𝐵𝑖𝑗 ≠ 0, fuerzas en el plano producen deformaciones de flexión y torsionales. 𝐷𝑖𝑗: Es la matriz de rigidez la cual relaciona momentos con las curvaturas. 3.6. Relaciones inversas de carga-deformación Como se dijo anteriormente los laminados multidireccionales se caracterizan de presentar discontinuidades de esfuerzo de lámina a lámina, es por eso que se prefiere trabajar con las deformaciones, las cuales son continuas a través del espesor del laminado. Por esta razón se requiere invertir las relaciones de Carga-Deformación de las ecuaciones (A.34) y (A.35), y expresar las deformaciones y curvaturas como funciones de las fuerzas o momentos aplicados como a continuación: [ 𝜀𝑥 0 𝜀𝑦 0 𝛾𝑠 0 ] = [ 𝑎𝑥𝑥 𝑎𝑥𝑦 𝑎𝑥𝑠 𝑎𝑦𝑥 𝑎𝑦𝑦 𝑎𝑦𝑠 𝑎𝑠𝑥 𝑎𝑠𝑦 𝑎𝑠𝑠 ] . [ 𝑁𝑥 𝑁𝑦 𝑁𝑠 ] + [ 𝑏𝑥𝑥 𝑏𝑥𝑦 𝑏𝑥𝑠 𝑏𝑦𝑥 𝑏𝑦𝑦 𝑏𝑦𝑠 𝑏𝑠𝑥 𝑏𝑠𝑦 𝑏𝑠𝑠 ] . [ 𝑀𝑥 𝑀𝑦 𝑀𝑠 ] (A.36) [ ĸ𝑥 ĸ𝑦 ĸ𝑠 ] = [ 𝑐𝑥𝑥 𝑐𝑥𝑦 𝑐𝑥𝑠 𝑐𝑦𝑥 𝑐𝑦𝑦 𝑐𝑦𝑠 𝑐𝑠𝑥 𝑐𝑠𝑦 𝑐𝑠𝑠 ] . [ 𝑁𝑥 𝑁𝑦 𝑁𝑠 ] + [ 𝑑𝑥𝑥 𝑑𝑥𝑦 𝑑𝑥𝑠 𝑑𝑦𝑥 𝑑𝑦𝑦 𝑑𝑦𝑠 𝑑𝑠𝑥 𝑑𝑠𝑦 𝑑𝑠𝑠 ] . [ 𝑀𝑥 𝑀𝑦 𝑀𝑠 ] (A.37) Donde las matrices “inversas” del laminado ( [𝒂]; [𝒃]; [𝒄] 𝑦 [𝒅] ) se obtienen según lo siguiente:  [𝑎] = [𝐴−1] − {[𝐵∗]. [𝐷∗−1]}. [𝐶∗]  [𝑏] = [𝐵∗]. [𝐷∗−1]  [𝑐] = −[𝐶∗]. [𝐷∗−1]  [𝑑] = [𝐷∗−1] 32 (A.38) Donde:  [𝐴−1] = 𝑖𝑛𝑣𝑒𝑟𝑠𝑎 𝑑𝑒 𝑙𝑎 𝑚𝑎𝑡𝑟𝑖𝑧 [𝐴]  [𝐵∗] = −[𝐴−1]. [𝐵]  [𝐶∗] = [𝐵]. [𝐴−1]  [𝐷∗] = [𝐷] − {[𝐵]. [𝐴−1]}. [𝐵] 3.7.propiedades de láminas a diferentes orientaciones Matriz de rigidez de la lámina Cuando se presentan láminas cuyos ejes principales no coinciden con los ejes de referencia del laminado, se deberá rotar la matriz de rigidez principal ([𝑄12]) reflejada a los ejes de referencia del laminado mediante las ecuaciones que se muestran en (A.39).  𝑄𝑥𝑥 = 𝑚 4. 𝑄11 + 𝑛 4. 𝑄22 + 2𝑚 2𝑛2. 𝑄12 + 4𝑚 2𝑛2. 𝑄66  𝑄𝑦𝑦 = 𝑛 4. 𝑄11 + 𝑚 4. 𝑄22 + 2𝑚 2𝑛2. 𝑄12 + 4𝑚 2𝑛2. 𝑄66  𝑄𝑥𝑦 = 𝑚 2𝑛2. 𝑄11 + 𝑚 2𝑛2. 𝑄22 + (𝑚 4 + 𝑛4). 𝑄12 − 4𝑚 2𝑛2. 𝑄66  𝑄𝑥𝑠 = 𝑚 3𝑛. 𝑄11 − 𝑚𝑛 3. 𝑄22 + (𝑚𝑛 3 − 𝑚3𝑛). 𝑄12 + 2(𝑚𝑛 3 − 𝑚3𝑛). 𝑄66  𝑄𝑦𝑠 = 𝑚𝑛 3. 𝑄11 − 𝑚 3𝑛. 𝑄22 + (𝑚 3𝑛 − 𝑚𝑛3). 𝑄12 + 2(𝑚 3𝑛 − 𝑚𝑛3). 𝑄66  𝑄𝑠𝑠 = 𝑚 2𝑛2. 𝑄11 + 𝑚 2𝑛2. 𝑄22 − 2𝑚 2𝑛2. 𝑄12 + (𝑚 2 − 𝑛2)2. 𝑄66 (A.39) Dónde: [𝑚 = 𝑐𝑜𝑠 θ𝑘 ; n = sen θ𝑘], (θ𝑘) es el ángulo orientado que se encuentra la lámina “k” medido positivo en sentido antihorario desde el eje x. Deformaciones de la lámina Las láminas orientadas un ángulo (θ𝑘) respecto del eje –x, deberá usarse las siguientes ecuaciones (A.40), para poder hallar las deformaciones principales en la lámina a partir de las deformaciones halladas en los ejes (x, y). [ 𝜀1 𝜀2 1 2 𝛾6 ] = [𝑇]. [ 𝜀𝑥 𝜀𝑦 1 2 𝛾𝑠 ] (A.40) Dónde: [𝑇] = [ 𝑚2 𝑛2 2𝑚𝑛 𝑛2 𝑚2 −2𝑚𝑛 −𝑚𝑛 𝑚𝑛 𝑚2 − 𝑛2 ] [𝑚 = 𝑐𝑜𝑠 θ𝑘 ; n = sen θ𝑘] 33 ANEXO 4 SIMULACIÓN NUMÉRICA DE VERIFICACIÓN USANDO MATERIAL ORTOTRÓPICO MEDIANTE LA TEORÍA CLÁSICA DE LAMINADO 1. Análisis usando la Teoría Clásica de Laminado Se realiza simulaciones a laminados simétrico y asimétrico a fin de validar la simulación objetivo (hélice del aerogenerador). Estas simulaciones de validación con sus respectivas validaciones analíticas y así verificar la teoría clásica de laminado. Finalmente realizar un análisis y discusión de resultados de la simulación de la hélice y la simulación de validación. 2. Validación con simulación numérica usando ANSYS APDL Estas simulaciones, como ya se mencionó líneas atrás se basa en un laminado simétrico (CASO 1) y un laminado antisimétrico (CASO 2), que se describe a continuación.  CASO 1: Para este caso CASO 1, el laminado es simétrico el cual es denotado de la siguiente manera: [𝟎/𝟗𝟎/𝟗𝟎/𝟎] o [𝟎/𝟗𝟎]𝒔.  CASO 2: Para este CASO 2. El laminado es antisimétrico el cual es denotado de la siguiente manera: [𝟎/𝟗𝟎]. 2.1. Simulación de laminado simétrico (CASO 1)  Etapa de Pre – procesamiento En esta etapa se define los parámetros de entrada que servirá para realizar la simulación numérica como propiedades del material y la geometría a utilizar, por otro lado también se realiza la discretización o mallado del modelo y finalmente las condiciones de borde, las cuales representan las restricciones del modelo. Definición de datos de entrada Como se mencionó líneas atrás este laminado es simétrico como se aprecia en la figura. Además este laminado de fibra de carbono/epoxy posee las siguientes propiedades elásticas. Las dimensiones de la lámina son de 100 mm de largo con 50 mm de ancho y una altura o espesor por lámina de 0.125 mm como se muestra en la figura A4.1. Tabla A4.1. Propiedades del laminado simétrico [𝟎/𝟗𝟎]𝒔. Propiedades Magnitud 𝐸1 140 kN mm 2⁄ 𝐸2 10 kN mm 2⁄ 𝐺12 5 kN mm 2⁄ 𝜈12 0.3 34 Figura A4.1. Valores de las dimensiones del laminado [𝟎/𝟗𝟎]𝒔 Del ANSYS APDL se define las láminas que conforma el laminado. La figura A4.2 muestra detalladamente la orientación de las fibras mediante el ángulo Theta, también muestra el número de láminas, que en este caso son 4, que están superpuestas unas a otras y finalmente el número de materiales “2” (en realidad es uno); para efectos del software se definen 2 materiales con sus propiedades ortotrópicas cruzadas. 𝐸1 = 140 kN mm 2⁄ y para el material 2, 𝐸1 = 10 kN mm 2⁄ . Figura A4.2. Valores de las dimensiones del laminado [𝟎/𝟗𝟎]𝒔 Parámetros de mallado El mallado se hizo con un elemento plano llamado Shell 181 con 4 nodos, usado para placas. El refinamiento se hizo usando Refine at nodes de grado 2. El mallado posee 162 elementos. Con un tamaño promedio de 5.5 mm por cada elemento. En la figura A4.4 se puede apreciar el mallado del laminado, en donde se ve las 4 láminas que conforman todo el laminado. Figura A4.3. Elemento Shell 181 con 4 nodos. 35 Figura A4.4. Mallado del laminado. Parámetros de control de simulación (Condiciones de borde) En esta etapa de generación del modelo por elementos finitos se procede a poner las restricciones y solicitaciones para el modelo, es decir los grados de libertad y el estado de cargas con lo cual estará sometido el laminado.  Condición de soporte fijo: Esta condición se basa en fijar la lámina mediante un comando en ANSYS APDL restringiendo el desplazamiento de la lámina en el eje Y. Al hacer esto (restringir el modelo) se disminuye los grados de libertad del laminado.  Condición de carga: Para este caso el laminado es sometido a tracción pura con una fuerza de 5000N a lo largo de las caras paralelas al plano Y-Z. Para las validaciones analíticas de esta simulación se puede denotar esta fuerza como 100 𝐍 𝐦𝐦⁄ . Finalmente se aprecia en la figura A5.5 las dos condiciones que se necesitan para la simulación estructural usando ANSYS APDL. Figura A4.5. Condición de soporte fijo y carga.  Etapa de solución El método de elementos finitos es un método computacional, que permite estudiar el comportamiento de un material sometido a un estado de cargas por medio de una 36 aproximación. Luego de una serie de manipulaciones matemáticas, para el caso de elasticidad lineal y estático, el sistema se resuelve según la siguiente ecuación. 𝑲𝒒 = 𝑸 Donde K es la matriz de rigidez generalizada, Q es el vector de fuerzas globales generalizado y q es el vector de fuerzas de desplazamientos generalizados.  Etapa de Post - Procesamiento En esta última etapa de la simulación se obtienen, visualizan y analizan los resultados de la simulación. Específicamente se determinan las deformaciones y esfuerzos que se presentan en el laminado. Los resultados se pueden evaluar en forma cuantitativa y cualitativa para así facilitar el entendimiento del lector. Esfuerzos longitudal (eje X) y transversal (eje Y) En las figuras A4.6 y A4.7 se pueden apreciar tanto los esfuerzos longitudinales como los transversales, correspondientes a las direcciones principal 1 y principal 2, respectivamente. Los esfuerzos máximos son 351 MPa (dirección longitudinal) y 25.3 MPa (dirección transversal). Figura A4.6. Esfuerzo en la dirección principal 1 (Eje X) Figura A4.7. Esfuerzo en la dirección principal 2 (Eje Y) 37 2.2. Simulación de laminado asimétrico (CASO 2)  Etapa de Pre – procesamiento Al igual que para el CASO 1 la etapa de pre-pocesamiento es prácticamente la misma para el CASO 2, ya que la propiedades de las láminas unidireccionales son las mismas, la geometría del laminado y de cada lámina son las mismas e incluso las condiciones de contorno. Lo único que cambia son las disposición de las láminas, ya que este laminado al ser [𝟎 /𝟗𝟎 ] es antisimétrico. Definición de datos de entrada Las propiedades de cada lámina son las mismas, mas no la configuración de las mismas para el laminado. La configuración se aprecia en la figura A4.8 y las propiedades en la tabla A4.1. Figura A4.8. Valores de las dimensiones del laminado [𝟎𝟐/𝟗𝟎𝟐]. Debido a que los parámetros de malla y los parámetros de control de simulación (Condiciones de borde) son iguales al CASO 1, no vale la pena volver a nombrarlos. Por ello, se obvia estas etapas del pre – procesamiento y se muestra los resultados en el post – procesamiento.  Etapa de Post – Procesamiento Esfuerzos longitudal (eje X) y transversal (eje Y) En las figuras A5.9 y A5.11 se puede apreciar tanto los esfuerzos longitudinales como los transversales, correspondientes a las direcciones principal 1 y principal 2, respectivamente. Los esfuerzos máximos son 543 MPa (dirección longitudinal) y -61.6 MPa (dirección transversal - compresión). 38 Figura A4.9. Esfuerzo en dirección principal 1. Figura A4.10. Esfuerzo en dirección principal 1 (vista en plano Y-Z). Figura A4.11. Esfuerzo en dirección principal 2. Figura A4.12. Esfuerzo en dirección principal 2 (vista en plano Y-Z). 39 ANEXO 5 CÁLCULO ANALÍTICO USANDO LA TCL DE LOS LAMINADOS [𝟎/𝟗𝟎]Y [𝟎/𝟗𝟎]𝒔 Se parte de esta igualdad para desarrollar los cálculos analíticos. 𝛎𝟏𝟐 𝐄𝟏⁄ = 𝛎𝟐𝟏 𝐄𝟐⁄ Entonces, 𝛎𝟐𝟏 = 𝟎. 𝟑 𝐱 𝟏𝟎 𝟏𝟒𝟎⁄ = 𝟎. 𝟎𝟐𝟏 Los términos de la matriz de rigidez reducida son: 𝐐𝟏𝟏 = 𝐄𝟏 (𝟏 − 𝛎𝟏𝟐𝛎𝟐𝟏)⁄ = 𝟏𝟒𝟎/(𝟏 − 𝟎. 𝟑 𝐱 𝟎. 𝟎𝟐𝟏) = 𝟏𝟒𝟎. 𝟗 𝐤𝐍/𝐦𝐦 𝟐 𝐐𝟏𝟏 = 𝐄𝟐 (𝟏 − 𝛎𝟏𝟐𝛎𝟐𝟏)⁄ = 𝟏𝟎/(𝟏 − 𝟎. 𝟑 𝐱 𝟎. 𝟎𝟐𝟏) = 𝟏𝟎. 𝟏 𝐤𝐍/𝐦𝐦 𝟐 𝐐𝟏𝟏 = 𝐆𝟏𝟐 = 𝟓. 𝟎 𝐤𝐍/𝐦𝐦 𝟐 𝐐𝟏𝟏 = 𝛎𝟐𝟏𝐄𝟏 (𝟏 − 𝛎𝟏𝟐𝛎𝟐𝟏)⁄ = (𝟎. 𝟎𝟐𝟏 𝐱 𝟏𝟒𝟎)/(𝟏 − 𝟎. 𝟑 𝐱 𝟎. 𝟎𝟐𝟏) = 𝟑. 𝟎 𝐤𝐍/𝐦𝐦 𝟐 Por lo tanto, la matriz queda definida como: 𝐐 = [ 𝟏𝟒𝟎. 𝟗 𝟑. 𝟎 𝟎 𝟑 𝟏𝟎. 𝟏 𝟎 𝟎 𝟎 𝟓. 𝟎 ] Siguiente, los términos de la matriz de rigidez para la lámina de ángulo 𝟎∘ que se obtienen es: 𝐦 = 𝐜𝐨𝐬(𝟎∘) = 𝟏; 𝐦𝟐 = 𝟏; 𝐦𝟒 = 𝟏 𝐧 = 𝐬𝐞𝐧(𝟎∘) = 𝟎; 𝐧𝟐 = 𝟎; 𝐧𝟒 = 𝟎 𝐦𝟐𝐧𝟐 = 𝟎; 𝟐𝐦𝟐𝐧𝟐 = 𝟎; 𝟒𝐦𝟐𝐧𝟐 = 𝟎 𝐦𝟐 − 𝐧𝟐 = 𝟏 𝐦𝟒 + 𝐧𝟒 = 𝟏 𝐦𝟑𝐧 = 𝟎; 𝐦𝐧𝟑 = 𝟎 𝐦𝟑𝐧 − 𝐦𝐧𝟑 = 𝟎; 𝟐(𝐦𝟑𝐧 − 𝐦𝐧𝟑) = 𝟎 𝐦𝐧𝟑 − 𝐦𝟑𝐧 = 𝟎; 2(𝐦𝐧𝟑 − 𝐦𝟑𝐧) = 𝟎 Sustituyendo los valores trigonométricos se obtiene: ?̅?𝟏𝟏 = 𝐦 𝟒𝐐𝟏𝟏 + 𝐧 𝟒𝐐𝟐𝟐 + 𝟐𝐦 𝟐𝐧𝟐𝐐𝟏𝟐 + 𝟒𝐦 𝟐𝐧𝟐𝐐𝟑𝟑 = (𝟏 𝐱 𝟏𝟒𝟎. 𝟗) + (𝟏 𝐱 𝟏𝟎. 𝟏) + (𝟎 𝐱 𝟑. 𝟎) + (𝟎 𝐱 𝟓. 𝟎) = 𝟏𝟒𝟎. 𝟗 𝐤𝐍/𝐦𝐦𝟐 = 𝐐𝟏𝟏 ?̅?𝟐𝟐 = 𝐧 𝟒𝐐𝟏𝟏 + 𝐦 𝟒𝐐𝟐𝟐 + 𝟐𝐦 𝟐𝐧𝟐𝐐𝟏𝟐 + 𝟒𝐦 𝟐𝐧𝟐𝐐𝟑𝟑 = (𝟎 𝐱 𝟏𝟒𝟎. 𝟗) + (𝟏 𝐱 𝟏𝟎. 𝟏) + (𝟎 𝐱 𝟑. 𝟎) + (𝟎 𝐱 𝟓. 𝟎) = 𝟏𝟎. 𝟏 𝐤𝐍/𝐦𝐦𝟐 = 𝐐𝟐𝟐 ?̅?𝟑𝟑 = 𝐦 𝟐𝐧𝟐𝐐𝟏𝟏 + 𝐦 𝟐𝐧𝟐𝐐𝟐𝟐 − 𝟐𝐦 𝟐𝐧𝟐𝐐𝟏𝟐 + (𝐦 𝟐 − 𝐧𝟐)𝐐𝟑𝟑 40 = (𝟎 𝐱 𝟏𝟒𝟎. 𝟗) + (𝟎 𝐱 𝟏𝟎. 𝟏) + (𝟎 𝐱 𝟑. 𝟎) + (𝟏𝟐 𝐱 𝟓. 𝟎) = 𝟓. 𝟎 𝐤𝐍/𝐦𝐦𝟐 = 𝐐𝟐𝟐 ?̅?𝟏𝟐 = 𝐦 𝟐𝐧𝟐𝐐𝟏𝟏 + 𝐦 𝟐𝐧𝟐𝐐𝟐𝟐 + (𝐦 𝟒 + 𝐧𝟒)𝐐𝟏𝟐 − 𝟒𝐦 𝟐𝐧𝟐𝐐𝟑𝟑 = (𝟎 𝐱 𝟏𝟒𝟎. 𝟗) + (𝟎 𝐱 𝟏𝟎. 𝟏) + (𝟏 𝐱 𝟑. 𝟎) + (𝟎 𝐱 𝟓. 𝟎) = 𝟑. 𝟎 𝐤𝐍/𝐦𝐦𝟐 = 𝐐𝟐𝟐 Entonces, la matriz de rigidez para la lámina de 𝟎∘ queda con los siguientes términos. (?̅?)𝟎∘ = [ 𝟏𝟒𝟎. 𝟗 𝟑. 𝟎 𝟎 𝟑 𝟏𝟎. 𝟏 𝟎 𝟎 𝟎 𝟓. 𝟎 ] Como lo mencionamos si fuera el caso de un laminado [𝟎𝟐]𝒔. Se obtendría para cada lámina lo siguiente: Tabla A5.1. Componentes de la matriz de rigidez para cada lámina[02]𝑠. Lámina 𝜽∘ ?̅?𝟏𝟏 ?̅?𝟐𝟐 ?̅?𝟑𝟑 ?̅?𝟏𝟐 ?̅?𝟏𝟑 ?̅?𝟐𝟑 1 0 140.9 10.1 5.0 3.0 0 0 2 0 140.9 10.1 5.0 3.0 0 0 3 0 140.9 10.1 5.0 3.0 0 0 4 0 140.9 10.1 5.0 3.0 0 0 Para validar analíticamente el CASO 1 (Simulación de validación) debemos hacer el estudio para el laminado [𝟎/𝟗𝟎]𝒔. Entonces, se plantea los términos de la matriz de rigidez para la lámina de 𝟗𝟎∘ 𝐦 = 𝐜𝐨𝐬𝟗(𝟎∘) = 𝟎; 𝐦𝟐 = 𝟎; 𝐦𝟒 = 𝟎 𝐧 = 𝐬𝐞𝐧(𝟗𝟎∘) = 𝟏; 𝐧𝟐 = 𝟏; 𝐧𝟒 = 𝟏 𝐦𝟐𝐧𝟐 = 𝟎; 𝟐𝐦𝟐𝐧𝟐 = 𝟎; 𝟒𝐦𝟐𝐧𝟐 = 𝟎 𝐦𝟐 − 𝐧𝟐 = −𝟏 𝐦𝟒 + 𝐧𝟒 = 𝟏 𝐦𝟑𝐧 = 𝟎; 𝐦𝐧𝟑 = 𝟎 𝐦𝟑𝐧 − 𝐦𝐧𝟑 = 𝟎; 𝟐(𝐦𝟑𝐧 − 𝐦𝐧𝟑) = 𝟎 𝐦𝐧𝟑 − 𝐦𝟑𝐧 = 𝟎; 2(𝐦𝐧𝟑 − 𝐦𝟑𝐧) = 𝟎 Sustituyendo los valores trigonométricos se obtiene: ?̅?𝟏𝟏 = 𝐦 𝟒𝐐𝟏𝟏 + 𝐧 𝟒𝐐𝟐𝟐 + 𝟐𝐦 𝟐𝐧𝟐𝐐𝟏𝟐 + 𝟒𝐦 𝟐𝐧𝟐𝐐𝟑𝟑 = (𝟎 𝐱 𝟏𝟒𝟎. 𝟗) + (𝟏 𝐱 𝟏𝟎. 𝟏) + (𝟎 𝐱 𝟑. 𝟎) + (𝟎 𝐱 𝟓. 𝟎) = 𝟏𝟎. 𝟏 𝐤𝐍/𝐦𝐦𝟐 = 𝐐𝟏𝟏 ?̅?𝟐𝟐 = 𝐧 𝟒𝐐𝟏𝟏 + 𝐦 𝟒𝐐𝟐𝟐 + 𝟐𝐦 𝟐𝐧𝟐𝐐𝟏𝟐 + 𝟒𝐦 𝟐𝐧𝟐𝐐𝟑𝟑 = (𝟏 𝐱 𝟏𝟒𝟎. 𝟗) + (𝟎 𝐱 𝟏𝟎. 𝟏) + (𝟎 𝐱 𝟑. 𝟎) + (𝟎 𝐱 𝟓. 𝟎) = 𝟏𝟒𝟎. 𝟗 𝐤𝐍/𝐦𝐦𝟐 = 𝐐𝟐𝟐 ?̅?𝟑𝟑 = 𝐦 𝟐𝐧𝟐𝐐𝟏𝟏 + 𝐦 𝟐𝐧𝟐𝐐𝟐𝟐 − 𝟐𝐦 𝟐𝐧𝟐𝐐𝟏𝟐 + (𝐦 𝟐 − 𝐧𝟐)𝐐𝟑𝟑 = (𝟎 𝐱 𝟏𝟒𝟎. 𝟗) + (𝟎 𝐱 𝟏𝟎. 𝟏) + (𝟎 𝐱 𝟑. 𝟎) + (𝟏𝟐 𝐱 𝟓. 𝟎) = 𝟓. 𝟎 𝐤𝐍/𝐦𝐦𝟐 = 𝐐𝟐𝟐 ?̅?𝟏𝟐 = 𝐦 𝟐𝐧𝟐𝐐𝟏𝟏 + 𝐦 𝟐𝐧𝟐𝐐𝟐𝟐 + (𝐦 𝟒 + 𝐧𝟒)𝐐𝟏𝟐 − 𝟒𝐦 𝟐𝐧𝟐𝐐𝟑𝟑 = (𝟎 𝐱 𝟏𝟒𝟎. 𝟗) + (𝟎 𝐱 𝟏𝟎. 𝟏) + (𝟎 𝐱 𝟑. 𝟎) + ((−𝟏) 𝟐𝐱 𝟓. 𝟎) = 𝟑. 𝟎 𝐤𝐍/𝐦𝐦𝟐 = 𝐐𝟐𝟐 41 ?̅?𝟏𝟑 = 𝐦 𝟑𝐧 𝐐𝟏𝟏 − 𝐦 𝐧𝟑𝐐𝟐𝟐 + (𝐦𝐧 𝟑 − 𝐦𝟑𝐧)𝐐𝟏𝟐 + 𝟐(𝐦𝐧 𝟑 − 𝐦𝟑𝐧)𝐐𝟑𝟑 = (𝟎 𝐱 𝟏𝟒𝟎. 𝟗) − (𝟎 𝐱 𝟏𝟎. 𝟏) + (𝟎 𝐱 𝟑. 𝟎) + (𝟎 𝐱 𝟓. 𝟎) = 𝟎 𝐤𝐍/𝐦𝐦𝟐 = 𝐐𝟏𝟑 ?̅?𝟐𝟑 = 𝐦𝐧 𝟑 𝐐𝟏𝟏 − 𝐦 𝟑𝐧𝐐𝟐𝟐 + (𝐦 𝟑𝐧 − 𝐦𝐧𝟑)𝐐𝟏𝟐 + 𝟐(𝐦 𝟑𝐧 − 𝐦𝐧𝟑)𝐐𝟑𝟑 = (𝟎 𝐱 𝟏𝟒𝟎. 𝟗) − (𝟎 𝐱 𝟏𝟎. 𝟏) + (𝟎 𝐱 𝟑. 𝟎) + (𝟎 𝐱 𝟓. 𝟎) = 𝟎 𝐤𝐍/𝐦𝐦𝟐 = 𝐐𝟐𝟑 Entonces, la matriz de rigidez para la lámina de 9𝟎∘ queda con los siguientes términos. (?̅?)𝟗𝟎∘ = [ 𝟏𝟎. 𝟏 𝟑. 𝟎 𝟎 𝟑. 𝟎 𝟏𝟒𝟎. 𝟗 𝟎 𝟎 𝟎 𝟓. 𝟎 ] Entonces para el laminado laminado [𝟎/𝟗𝟎]𝒔 se obtiene lo siguiente: Tabla 5.2. Componentes de la matriz de rigidez para cada lámina [0/90]𝑠 Lámina 𝛉∘ ?̅?𝟏𝟏 ?̅?𝟐𝟐 ?̅?𝟑𝟑 ?̅?𝟏𝟐 ?̅?𝟏𝟑 ?̅?𝟐𝟑 1 0 140.9 10.1 5.0 3.0 0 0 2 90 10.1 140.9 5.0 3.0 0 0 3 90 10.1 140.9 5.0 3.0 0 0 4 0 140.9 10.1 5.0 3.0 0 0 Como el laminado [𝟎/𝟗𝟎]𝒔 es simétrico entonces la matriz de acoplamiento tensión-flexión [𝑩] = 0, por lo tanto solo se necesita calcular las matrices [𝑨] y [𝑫]; denominadas matriz de rigidez de tensión plana y matriz de rigidez de flexión. Para realizar los cálculos de la matriz de rigidez de tensión plana [𝑨] y la matriz de rigidez de flexión [𝑫], se necesita de 𝒕𝒑 𝒚 (𝒕𝒑?̅?𝒑 𝟐 + 𝒕𝒑 𝟑/𝟏𝟐), donde la figura 97 sirve de apoyo. Los valores para cada lámina están dados en la siguiente tabla. Tabla A5.3. Datos para cálculos de matrices A y D. Lámina 𝛉∘ 𝐭𝐩 ?̅?𝐩 𝐭𝐩?̅?𝐩 𝟐 + 𝐭𝐩 𝟑/𝟏𝟐 1 0 0.125 -0.01875 0.00456 2 90 0.125 -0.0625 0.00065 3 90 0.125 0.0625 0.00065 4 0 0.125 0.1875 0.00456 Para determinar los componentes de la matriz [𝑨] y [𝑫] se usan las siguientes ecuaciones: 𝐀𝐢𝐣 = ∑ 𝐭𝐩(?̅?𝐢𝐣)𝐩 𝐍 𝐩=𝟏 𝐃𝐢𝐣 = ∑(𝐭𝐩?̅?𝐩 𝟐 + 𝐭𝐩 𝟑 𝟏𝟐 )(?̅?𝐢𝐣)𝐩 𝐍 𝐩=𝟏 Con los valores de la tabla A6.3 y las ecuaciones anteriores se procede a calcular los componentes de cada matriz ([𝑨] 𝒚 [𝑫] ). 𝐀𝟏𝟏 = 𝟐[(𝟎. 𝟏𝟐𝟓 𝐱 𝟏𝟒𝟎. 𝟗)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟏𝟐𝟓 𝐱 𝟏𝟒𝟎. 𝟗)𝐋á𝐦𝐢𝐧𝐚 𝟐] = 𝟑𝟕. 𝟖 𝐤𝐍/𝐦𝐦 42 𝐀𝟐𝟐 = 𝟐[(𝟎. 𝟏𝟐𝟓 𝐱 𝟏𝟎. 𝟏)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟏𝟐𝟓 𝐱 𝟏𝟒𝟎. 𝟗)𝐋á𝐦𝐢𝐧𝐚 𝟐] = 𝟑𝟕. 𝟖 𝐤𝐍/𝐦𝐦 𝐀𝟑𝟑 = 𝟐[(𝟎. 𝟏𝟐𝟓 𝐱 𝟓. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟏𝟐𝟓 𝐱 𝟓. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟐] = 𝟐. 𝟓 𝐤𝐍/𝐦𝐦 𝐀𝟏𝟐 = 𝟐[(𝟎. 𝟏𝟐𝟓 𝐱 𝟑. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟏𝟐𝟓 𝐱 𝟑. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟐] = 𝟏. 𝟓 𝐤𝐍/𝐦𝐦 𝐀𝟏𝟑 = 𝐀𝟐𝟑 = 𝟎; como ?̅?𝟏𝟑 = ?̅?𝟐𝟑 = 𝟎 para lámina de 𝟎 ∘ y 𝟗𝟎∘ 𝐃𝟏𝟏 = 𝟐[(𝟎. 𝟎𝟎𝟒𝟓𝟔 𝐱 𝟏𝟒𝟎. 𝟗)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟏𝟎. 𝟏)𝐋á𝐦𝐢𝐧𝐚 𝟐] = 𝟏. 𝟐𝟗𝟖𝟏 𝐤𝐍/𝐦𝐦 𝐃𝟏𝟏 = 𝟐[(𝟎. 𝟎𝟎𝟒𝟓𝟔 𝐱 𝟏𝟎. 𝟏)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟏𝟒𝟎. 𝟗)𝐋á𝐦𝐢𝐧𝐚 𝟐] = 𝟎. 𝟐𝟕𝟓𝟑 𝐤𝐍/𝐦𝐦 𝐃𝟏𝟏 = 𝟐[(𝟎. 𝟎𝟎𝟒𝟓𝟔 𝐱 𝟓. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟓. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟐] = 𝟎. 𝟎𝟓𝟐𝟏 𝐤𝐍/𝐦𝐦 𝐃𝟏𝟏 = 𝟐[(𝟎. 𝟎𝟎𝟒𝟓𝟔 𝐱 𝟑. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟑. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟐] = 𝟎. 𝟎𝟑𝟏𝟑 𝐤𝐍/𝐦𝐦 𝐃𝟏𝟑 = 𝐃𝟐𝟑 = 𝟎; como ?̅?𝟏𝟑 = ?̅?𝟐𝟑 = 𝟎 para lámina de 𝟎 ∘ y 𝟗𝟎∘ Por lo tanto las matrices de rigidez de tensión plana y de flexión quedan de la siguiente manera: 𝐀 = [ 𝟑𝟕. 𝟖 𝟏. 𝟓 𝟎 𝟏. 𝟓 𝟑𝟕. 𝟖 𝟎 𝟎 𝟎 𝟐. 𝟓 ] kN/mm ; 𝑫 = [ 𝟏. 𝟐𝟗𝟖𝟏 𝟎. 𝟎𝟑𝟏𝟑 𝟎 𝟎. 𝟎𝟑𝟏𝟑 𝟎. 𝟐𝟕𝟓𝟑 𝟎 𝟎 𝟎 𝟎. 𝟎𝟓𝟐𝟏 ] kN/mm Teniendo las matrices, se procede a calcular las matrices inversas de [𝑨] 𝒚 [𝑫], denotadas como [𝒂] 𝒚 [𝒅], respectivamente. A continuación se calculas las componentes estas matrices inversas. Matriz inversa de [𝑨] y [𝑫] 𝐀𝟏𝟏 = 𝟑𝟕. 𝟖 𝐀𝟐𝟐 = 𝟑𝟕. 𝟖 𝐀𝟑𝟑 = 𝟐. 𝟓 𝐀𝟏𝟐 = 𝟏. 𝟓 Además, 𝐀𝐀 = 𝐀𝟏𝟏𝐀𝟐𝟐 − 𝐀𝟏𝟐 𝟐 = (𝟑𝟕. 𝟖)(𝟑𝟕. 𝟖) − 𝟏. 𝟓𝟐 = 𝟏𝟒𝟐𝟔. 𝟔 𝐚𝟏𝟏 = 𝐀𝟐𝟐 𝐀𝐀⁄ = 𝟑𝟕. 𝟖 𝟏𝟒𝟐𝟔. 𝟔⁄ = 𝟎. 𝟎𝟐𝟔𝟓 𝐚𝟐𝟐 = 𝐀𝟏𝟏 𝐀𝐀⁄ = 𝟑𝟕. 𝟖 𝟏𝟒𝟐𝟔. 𝟔⁄ = 𝟎. 𝟎𝟐𝟔𝟓 𝐚𝟑𝟑 = 𝟏 𝐀𝟑𝟑⁄ = 𝟏 𝟐. 𝟓⁄ = 𝟎. 𝟒𝟎𝟎𝟎 𝐚𝟏𝟐 = −𝐀𝟏𝟐 𝐀𝐀⁄ = −𝟏. 𝟓 𝟏𝟒𝟐𝟔. 𝟔⁄ = −𝟎.𝟎𝟎𝟏𝟏 𝐃𝐃 = 𝐃𝟏𝟏𝐃𝟐𝟐 − 𝐃𝟏𝟐 𝟐 = (𝟏. 𝟐𝟗𝟖𝟏)(𝟎. 𝟐𝟕𝟓𝟑) − (𝟎. 𝟎𝟑𝟏𝟑)𝟐 = 𝟎. 𝟑𝟓𝟔𝟒 𝐝𝟏𝟏 = 𝐃𝟐𝟐 𝐃𝐃⁄ = 𝟎. 𝟐𝟕𝟓𝟑 𝟎. 𝟑𝟓𝟔𝟒⁄ = 𝟎. 𝟕𝟕 𝐝𝟐𝟐 = 𝐃𝟏𝟏 𝐃𝐃⁄ = 𝟏. 𝟐𝟗𝟖𝟏 𝟎. 𝟑𝟓𝟔𝟒⁄ = 𝟎. 𝟎𝟐𝟔𝟓 𝐝𝟑𝟑 = 𝟏 𝐃𝟑𝟑⁄ = 𝟏 𝟎. 𝟎𝟓𝟐𝟏⁄ = 𝟏𝟗. 𝟏𝟗 𝐝𝟏𝟐 = −𝐃𝟏𝟐 𝐃𝐃⁄ = −𝟎. 𝟎𝟑𝟏𝟑 𝟎. 𝟑𝟓𝟔𝟒⁄ = −𝟎. 𝟎𝟗 𝐀𝟏𝟏 = 𝟑𝟕. 𝟖 𝑫𝟏𝟏 = 𝟏. 𝟐𝟗𝟖𝟏 43 𝐀𝟐𝟐 = 𝟑𝟕. 𝟖 𝐃𝟐𝟐 = 𝟎. 𝟐𝟕𝟓𝟑 𝐀𝟑𝟑 = 𝟐. 𝟓 𝐃𝟑𝟑 = 𝟎. 𝟎𝟓𝟐𝟏 𝐀𝟏𝟐 = 𝟏. 𝟓 𝐃𝟏𝟐 = 𝟎. 𝟎𝟑𝟏𝟑 [𝒂] = [ 𝟎. 𝟎𝟐𝟔𝟓 −𝟎. 𝟎𝟎𝟏𝟏 𝟎 −𝟎. 𝟎𝟎𝟏𝟏 𝟎. 𝟎𝟐𝟔𝟓 𝟎 𝟎 𝟎 𝟎. 𝟒𝟎𝟎𝟎 ] kN/mm [𝐝] = [ 𝟎. 𝟕𝟕 −𝟎. 𝟎𝟗 𝟎 −𝟎. 𝟎𝟗 𝟑. 𝟔𝟒 𝟎 𝟎 𝟎 𝟏𝟗. 𝟏𝟗 ] kN/mm Usando las siguientes ecuaciones: [𝑴] = 𝟎 [ 𝛆𝐱 𝟎 𝛆𝐲 𝟎 𝛄𝐬 𝟎 ] = [ 𝐚𝐱𝐱 𝐚𝐱𝐲 𝐚𝐱𝐬 𝐚𝐲𝐱 𝐚𝐲𝐲 𝐚𝐲𝐬 𝐚𝐬𝐱 𝐚𝐬𝐲 𝐚𝐬𝐬 ] . [ 𝐍𝐱 𝐍𝐲 𝐍𝐬 ] + [ 𝐛𝐱𝐱 𝐛𝐱𝐲 𝐛𝐱𝐬 𝐛𝐲𝐱 𝐛𝐲𝐲 𝐛𝐲𝐬 𝐛𝐬𝐱 𝐛𝐬𝐲 𝐛𝐬𝐬 ] . [ 𝐌𝐱 𝐌𝐲 𝐌𝐬 ] Por lo tanto, [ 𝛆𝐱 𝟎 𝛆𝐲 𝟎 𝛄𝐬 𝟎 ] = [ 𝟎. 𝟎𝟐𝟔𝟓 −𝟎. 𝟎𝟏𝟏 𝟎 −𝟎. 𝟎𝟏𝟏 𝟎. 𝟎𝟐𝟔𝟓 𝟎 𝟎 𝟎 𝟎. 𝟒𝟎𝟎𝟔 ] . [ 𝐍𝐱 𝟎 𝟎 ] 𝐱 𝟏𝟎−𝟑 Que da lo siguiente: 𝜺𝒙 𝟎 = 𝟐𝟔𝟓𝟎 𝒙 𝟏𝟎−𝟔 𝜺𝒚 𝟎 = −𝟏𝟏𝟎 𝒙 𝟏𝟎−𝟔 𝜺𝒛 𝟎 = 𝟎 Estos son los valores de deformación a través de todo el espesor del laminado en las cuatro láminas como se ve en la figura A5.1. Para hacer el estudio en cada lámina se usa las ecuaciones de transformación, las cuales son dependientes de la orientación de cada lámina. [ 𝛆𝟏 𝛆𝟐 𝛄𝟏𝟐 ] = [ 𝐦𝟐 𝐧𝟐 𝐦𝐧 𝐧𝟐 𝐦𝟐 −𝐦𝐧 −𝟐𝐦𝐧 𝟐𝐦𝐧 𝐦𝟐 − 𝐧𝟐 ] . [ 𝛆𝐱 𝟎 𝛆𝐲 𝟎 𝛄𝐬 𝟎 ] 𝐱 𝟏𝟎−𝟔 Figura A5.1. Diagrama de deformaciones según ejes X, Y y Z. 44 Lámina 1 a 𝟎∘ 𝐦 = 𝐜𝐨𝐬(𝟎∘) = 𝟏 y 𝐧 = 𝐬𝐞𝐧(𝟎∘) = 𝟎 [ 𝛆𝟏 𝛆𝟐 𝛄𝟏𝟐 ] = [ 𝟏 𝟎 𝟎 𝟎 𝟏 𝟎 𝟎 𝟎 𝟏 ] . [ 𝛆𝐱 𝟎 𝛆𝐲 𝟎 𝛄𝐬 𝟎 ] 𝐱 𝟏𝟎−𝟔 [ 𝛆𝟏 𝛆𝟐 𝛄𝟏𝟐 ] = [ 𝟐𝟔𝟓𝟎 𝒙 𝟏𝟎−𝟔 −𝟏𝟏𝟎 𝒙 𝟏𝟎−𝟔 𝟎 ] [ σ1 σ2 τ12 ] Lámina 1 = [ 140.9 3.0 0 3.0 10.1 0 0 0 5.0 ] Lámina 1 . [ 2650 −110 0 ] Lámina 1 x10−6 x 103 N/mm2 De donde se obtiene lo siguientes esfuerzos: [ 𝛔𝟏 𝛔𝟐 𝛕𝟏𝟐 ] 𝐋á𝐦𝐢𝐧𝐚 𝟏 = [ 𝟑𝟕𝟑 𝟕 𝟎 ] 𝐋á𝐦𝐢𝐧𝐚 𝟏 𝐍/𝐦𝐦𝟐 Lámina 2 a 𝟗𝟎∘ 𝐦 = 𝐜𝐨𝐬(𝟗𝟎∘) = 𝟏 y 𝐧 = 𝐬𝐞𝐧(𝟗𝟎∘) = 𝟏 [ 𝛆𝟏 𝛆𝟐 𝛄𝟏𝟐 ] = [ 𝟎 𝟏 𝟎 𝟏 𝟎 𝟎 𝟎 𝟎 −𝟏 ] . [ 𝛆𝐱 𝟎 𝛆𝐲 𝟎 𝛄𝐬 𝟎 ] 𝐱 𝟏𝟎−𝟔 [ 𝛆𝟏 𝛆𝟐 𝛄𝟏𝟐 ] = [ −𝟏𝟏𝟎 𝒙 𝟏𝟎−𝟔 𝟐𝟔𝟓𝟎 𝒙 𝟏𝟎−𝟔 𝟎 ] [ σ1 σ2 τ12 ] Lámina 2 = [ 140.9 3.0 0 3.0 10.1 0 0 0 5.0 ] Lámina 2 . [ −110 2650 0 ] Lámina 2 x10−6 x 103 N/mm2 De donde se obtiene lo siguientes esfuerzos: [ 𝛔𝟏 𝛔𝟐 𝛕𝟏𝟐 ] 𝐋á𝐦𝐢𝐧𝐚 𝟐 = [ −𝟕 𝟐𝟔 𝟎 ] 𝐋á𝐦𝐢𝐧𝐚 𝟐 𝐍/𝐦𝐦𝟐 Figura A5.2. Diagrama de deformaciones totales en direcciones principales. 3.7.1.1.1. Cálculos analíticos para CASO 2 El CASO 2 se basa en un laminado compuesto por 2 láminas unidireccionales, la primera orientada a 𝟎∘ y la segunda a 𝟗𝟎∘. A continuación se plantea la metodología de cálculo a fin 45 de entender el comportamiento de estos tipos de materiales (ortotrópicos) para este caso en particular. Para el CASO 1 se calculó la matriz de rigidez correspondiente a cada lámina. Para este caso al ser las mismas láminas, sus matrices de rigidez también son iguales para cada lámina. Ver la tabla A5.4. Tabla A5.4. Componentes de la matriz de rigidez para cada lámina. Lámina 𝛉∘ ?̅?𝟏𝟏 ?̅?𝟐𝟐 ?̅?𝟑𝟑 ?̅?𝟏𝟐 ?̅?𝟏𝟑 ?̅?𝟐𝟑 1 0 140.9 10.1 5.0 3.0 0 0 2 90 10.1 140.9 5.0 3.0 0 0 En forma matricial: (?̅?)𝟎∘ = [ 𝟏𝟒𝟎. 𝟗 𝟑. 𝟎 𝟎 𝟑 𝟏𝟎. 𝟏 𝟎 𝟎 𝟎 𝟓. 𝟎 ] 𝐤𝐍/𝐦𝐦𝟐 ; (?̅?)𝟗𝟎∘ = [ 𝟏𝟎. 𝟏 𝟑. 𝟎 𝟎 𝟑. 𝟎 𝟏𝟒𝟎. 𝟗 𝟎 𝟎 𝟎 𝟓. 𝟎 ] 𝐤𝐍/𝐦𝐦𝟐 Al ser el laminado antisimétrico aparece la matriz de acoplamiento [𝐁]. Por tanto procedemos a calcular [𝐀] , [𝐁] 𝐲 [𝐃]. 𝑨𝒊𝒋 = ∑ 𝒕𝒑(?̅?𝒊𝒋)𝒑 𝑵 𝒑=𝟏 𝑩𝒊𝒋 = −∑ 𝒕𝒑?̅?𝒑(?̅?𝒊𝒋)𝒑 𝑵 𝒑=𝟏 𝑫𝒊𝒋 = ∑ (𝒕𝒑?̅?𝒑 𝟐 + 𝒕𝒑 𝟑 𝟏𝟐 )(?̅?𝒊𝒋)𝒑 𝑵 𝒑=𝟏 Para la realización de los cálculos los datos necesarios para las ecuaciones mostradas se muestran en la tabla A5.5. Tabla A5.5 Datos para los cálculos de matriz A, B y C. Lámina 𝛉∘ 𝐭𝐩 ?̅?𝐩 𝐭𝐩?̅?𝐩 𝐭𝐩?̅?𝐩 𝟐 + 𝐭𝐩 𝟑/𝟏𝟐 1 0 0.125 -0.0625 -0.00781 0.00065 2 90 0.125 0.0625 0.00781 0.00065 𝐀𝟏𝟏 = (𝟎. 𝟏𝟐𝟓 𝐱 𝟏𝟒𝟎. 𝟗)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟏𝟐𝟓 𝐱 𝟏𝟎. 𝟏)𝐋á𝐦𝐢𝐧𝐚 𝟐 = 𝟏𝟖. 𝟗 𝐤𝐍/𝐦𝐦 𝐀𝟐𝟐 = (𝟎. 𝟏𝟐𝟓 𝐱 𝟏𝟎. 𝟏)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟏𝟐𝟓 𝐱 𝟏𝟒𝟎. 𝟗)𝐋á𝐦𝐢𝐧𝐚 𝟐 = 𝟏𝟖. 𝟗 𝐤𝐍/𝐦𝐦 𝐀𝟑𝟑 = (𝟎. 𝟏𝟐𝟓 𝐱 𝟓. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟏𝟐𝟓 𝐱 𝟓. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟐 = 𝟏. 𝟑 𝐤𝐍/𝐦𝐦 𝐀𝟏𝟐 = (𝟎. 𝟏𝟐𝟓 𝐱 𝟑. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟏𝟐𝟓 𝐱 𝟑. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟐 = 𝟎. 𝟖 𝐤𝐍/𝐦𝐦 𝐀𝟏𝟑 = 𝐀𝟐𝟑 = 𝟎; como ?̅?𝟏𝟑 = ?̅?𝟐𝟑 = 𝟎 para lámina de 𝟎 ∘ y 𝟗𝟎∘ 𝐁𝟏𝟏 = −[(−𝟎. 𝟎𝟎𝟕𝟖𝟏 𝐱 𝟏𝟒𝟎. 𝟗)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟕𝟖𝟏 𝐱 𝟏𝟎. 𝟏)𝐋á𝐦𝐢𝐧𝐚 𝟐] = 𝟏. 𝟎 𝐤𝐍/𝐦𝐦 𝐁𝟐𝟐 = −[(−𝟎. 𝟎𝟎𝟕𝟖𝟏 𝐱 𝟏𝟎. 𝟏)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟕𝟖𝟏 𝐱 𝟏𝟒𝟎. 𝟗)𝐋á𝐦𝐢𝐧𝐚 𝟐] = −𝟏. 𝟎 𝐤𝐍/𝐦𝐦 𝐁𝟑𝟑 = −[(−𝟎. 𝟎𝟎𝟕𝟖𝟏 𝐱 𝟓. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟕𝟖𝟏 𝐱 𝟓. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟐] = 𝟎 𝐤𝐍/𝐦𝐦 𝐁𝟏𝟐 = −[−𝟎. 𝟎𝟎𝟕𝟖𝟏 𝐱 𝟑. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟕𝟖𝟏 𝐱 𝟑. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟐] = 𝟎 𝐤𝐍/𝐦𝐦 𝐁𝟏𝟑 = 𝐃𝟐𝟑 = 𝟎; como ?̅?𝟏𝟑 = ?̅?𝟐𝟑 = 𝟎 para lámina de 𝟎 ∘ y 𝟗𝟎∘ 46 𝐃𝟏𝟏 = (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟏𝟒𝟎. 𝟗)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟏𝟎. 𝟏)𝐋á𝐦𝐢𝐧𝐚 𝟐 = 𝟎. 𝟎𝟗𝟖𝟐 𝐤𝐍/𝐦𝐦 𝐃𝟏𝟏 = (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟏𝟎. 𝟏)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟏𝟒𝟎. 𝟗)𝐋á𝐦𝐢𝐧𝐚 𝟐 = 𝟎. 𝟎𝟗𝟖𝟐 𝐤𝐍/𝐦𝐦 𝐃𝟏𝟏 = (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟓. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟑. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟐 = 𝟎. 𝟎𝟎𝟔𝟓 𝐤𝐍/𝐦𝐦 𝐃𝟏𝟏 = (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟑. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟏 + (𝟎. 𝟎𝟎𝟎𝟔𝟓 𝐱 𝟓. 𝟎)𝐋á𝐦𝐢𝐧𝐚 𝟐 = 𝟎. 𝟎𝟎𝟑𝟗 𝐤𝐍/𝐦𝐦 𝐃𝟏𝟑 = 𝐃𝟐𝟑 = 𝟎; como ?̅?𝟏𝟑 = ?̅?𝟐𝟑 = 𝟎 para lámina de 𝟎 ∘ y 𝟗𝟎∘ 𝐀 = [ 𝟏𝟖. 𝟗 𝟎. 𝟖 𝟎 𝟎. 𝟖 𝟏𝟖. 𝟗 𝟎 𝟎 𝟎 𝟏. 𝟑 ] kN/mm 𝐁 = [ 𝟏 𝟎 𝟎 𝟎 −𝟏 𝟎 𝟎 𝟎 𝟎 ] kN/mm 𝑫 = [ 𝟎. 𝟎𝟗𝟖𝟐 𝟎. 𝟎𝟎𝟑𝟗 𝟎 𝟎. 𝟎𝟎𝟑𝟗 𝟎. 𝟎𝟗𝟖𝟐 𝟎 𝟎 𝟎 𝟎. 𝟎𝟎𝟔𝟓 ] kN/mm Usando las ecuaciones se tiene lo siguiente: [ 100 0 0 ] = [ 18900 800 0 800 18900 0 0 0 1300 ] . [ 𝜀𝑥 0 𝜀𝑦 0 𝛾𝑠 0 ] + [ 1000 0 0 0 −1000 0 0 0 0 ] . [ ĸ𝑥 ĸ𝑦 ĸ𝑠 ] [ 𝑀𝑥 𝑀𝑦 𝑀𝑠 ] = [ 1000 0 0 0 −1000 0 0 0 0 ] . [ 𝜀𝑥 0 𝜀𝑦 0 𝛾𝑠 0 ] + [ 98.2 3.9 0 3.9 98.2 0 0 0 6.5 ] . [ ĸ𝑥 ĸ𝑦 ĸ𝑠 ] De esas dos ecuaciones se plantean las siguientes: 100 = 18900 𝛆𝐱 𝟎 + 800 𝛆𝐲 𝟎 + 1000 ĸ𝐱 (a) 0 = 800 𝛆𝐱 𝟎+ 18900 𝛆𝐲 𝟎- 1000 ĸ𝐱 (b) 0 = 1300 𝛆𝐱𝐲 𝟎 (c) 0 = 1000 𝛆𝐱 𝟎+ 98.2 ĸ𝐱 + 3.9 ĸ𝐲 (d) 0 = - 1000 𝛆𝐲 𝟎+ 3.9 ĸ𝐱 + 98.2 ĸ𝐲 (e) 0 = 6.5 ĸ𝐱𝐲 (f) De (f): ĸ𝐱𝐲 = 0 De (c): 𝛆𝐱𝐲 𝟎 = 0 Entonces quedan las ecuaciones (a), (b), (c) y (d). Las cuales poseen 4 incognitas. Por tanto se obtiene lo siguiente: 𝛆𝐱 𝟎 = 11.517 x 𝟏𝟎−𝟑 mm 𝛆𝐲 𝟎 = -0.523 x 𝟏𝟎−𝟑 mm ĸ𝐱 = -117.256 x 𝟏𝟎 −𝟑 1/mm ĸ𝐲 = -0.667 x 𝟏𝟎 −𝟑 1/mm 47 Lámina 1 a 𝟎∘ 𝛆𝐱 = 𝛆 𝟎 − 𝐳𝐤 𝛆𝐱 = [𝟏𝟏𝟓𝟏𝟕 − (−𝟎. 𝟎𝟔𝟐𝟓)(−𝟏𝟏𝟕. 𝟐𝟓𝟔)] 𝒙 𝟏𝟎 −𝟑 = 𝟒. 𝟏𝟖𝟖𝟓 𝒙 𝟏𝟎−𝟑 = 𝛆𝟏 𝛆𝐲 = [−𝟎. 𝟓𝟐𝟑 − (−𝟎. 𝟎𝟔𝟐𝟓)(−𝟎. 𝟔𝟔𝟕)] 𝒙 𝟏𝟎 −𝟑 = −𝟎. 𝟓𝟔𝟒𝟕 𝒙 𝟏𝟎−𝟑 = 𝛆𝟐 𝛆𝐱𝐲 = 𝟎 = 𝛆𝟏𝟐 [ σ1 σ2 τ12 ] Lámina 1 = [ 140.9 3.0 0 3.0 10.1 0 0 0 5.0 ] Lámina 1 . [ 4.1885 −0.5647 0 ] Lámina 1 x10−6 x 103 N/mm2 De donde se obtiene lo siguientes esfuerzos: [ 𝛔𝟏 𝛔𝟐 𝛕𝟏𝟐 ] 𝐋á𝐦𝐢𝐧𝐚 𝟏 = [ 𝟓𝟖𝟖 𝟕 𝟎 ] 𝐋á𝐦𝐢𝐧𝐚 𝟏 𝐍/𝐦𝐦𝟐 Lámina 2 a 𝟗𝟎∘ 𝛆𝐱 = 𝛆 𝟎 − 𝐳𝐤 𝛆𝐱 = [𝟏𝟏𝟓𝟏𝟕 − (𝟎. 𝟎𝟔𝟐𝟓)(−𝟏𝟏𝟕. 𝟐𝟓𝟔)] 𝒙 𝟏𝟎 −𝟑 = 𝟏𝟖. 𝟖𝟖𝟒 𝒙 𝟏𝟎−𝟑 = 𝛆𝟐 𝛆𝒚 = [−𝟎. 𝟓𝟐𝟑 − (𝟎. 𝟎𝟔𝟐𝟓)(−𝟎. 𝟔𝟔𝟕)] 𝒙 𝟏𝟎 −𝟑 = −𝟎.𝟒𝟖𝟏𝟑 𝒙 𝟏𝟎−𝟑 = 𝛆𝟏 𝛆𝐱𝐲 = 𝟎 = 𝛆𝟏𝟐 [ σ1 σ2 τ12 ] Lámina 2 = [ 140.9 3.0 0 3.0 10.1 0 0 0 5.0 ] Lámina 2 . [ −0.4813 18.8455 0 ] Lámina 2 x10−6 x 103 N/mm2 De donde se obtiene lo siguientes esfuerzos: [ 𝛔𝟏 𝛔𝟐 𝛕𝟏𝟐 ] 𝐋á𝐦𝐢𝐧𝐚 𝟐 = [ −𝟏𝟏 𝟏𝟖𝟗 𝟎 ] 𝐋á𝐦𝐢𝐧𝐚 𝟐 𝐍/𝐦𝐦𝟐