You do not have permission to edit this page, for the following reason:

You are not allowed to execute the action you have requested.


You can view and copy the source of this page.

x
 
1
==1 Introducción==
2
3
Las estructuras de materiales compuestos de tipo laminado son especialmente sensibles a las tensiones fuera del plano debido a que las propiedades resistentes y de rigidez en la dirección del espesor del laminado son inferiores a las propiedades en el plano. Tensiones en dirección del espesor aparecen incluso aunque las cargas estén contenidas en el plano del laminado, por ejemplo, en puntos cercanos a un borde libre asociadas a las diferencias en el comportamiento entre dos láminas adyacentes con diferentes orientaciones y a la necesidad de verificar la condición de contorno. Estas tensiones pueden ser muy relevantes en la estimación de la carga de rotura de un laminado, por lo que una reducción de las mismas puede incrementar la resistencia del laminado.
4
5
Asociados al efecto de borde libre hay dos modos de fallo controlados por diferentes componentes del tensor de tensiones. Uno de los modos de fallo más críticos en los laminados es la delaminación que es controlado por la componente normal en dirección del espesor [1] o por tensiones tangenciales [2]. Este modo de fallo es responsable de una pérdida significativa de capacidad resistente de la estructura frente a cargas de compresión [2]. Otro de los modos de fallo es el fallo por cortadura que origina fisuras transversales, controlado en este caso por las componentes a cortadura interlaminar del tensor de tensiones. Determinar las componentes de tensión fuera del plano es por tanto de gran relevancia, así como estimar sus valores relativos para determinar el posible modo de fallo que puede aparecer.
6
7
El uso de modelos bidimensionales para abordar el estudio de estructuras de tipo laminado, como puede ser la teoría del laminado, no permite el cálculo de tensiones fuera del plano y, por tanto, puede hacer que se produzca un fallo prematuro en la estructura con resultados catastróficos. Los modelos tridimensionales aplicados de forma general a un problema arbitrario resultan a menudo complejos de abordar. Una aproximación muy extendida es el empleo del método de los elementos finitos [3-4], el cual requiere un gran afinamiento de la malla para recoger el fenómeno de concentración de tensiones. Otras metodologías son el método de los elementos de contorno [5] o las diferencias finitas [6]. Una alternativa es la estimación de las tensiones interlaminares mediante consideraciones de equilibrio que junto a hipótesis simplificativas permite resolver casos de carga simples y laminados simétricos.
8
9
El problema de las tensiones fuera del plano asociado a efectos de borde ha sido muy estudiado en la literatura (ver por ejemplo, los trabajos [7] y [8] que realizan una revisión sobre el estado del arte de este problema) desde los trabajos pioneros de Pipes y Pagano en 1970 [6]. Se han estudiado diferentes secuencias de apilamiento: 0/90, 90/0, +45/-45, cuasi-isótropos, etc [3-5]. Los resultados muestran que las tensiones interlaminares se ven modificadas significativamente con la secuencia de apilamiento. Por tanto, se hace necesario un estudio sistemático de diferentes orientaciones y secuencia de apilamiento que engloble las orientaciones empleadas por diferentes autores.
10
11
En este trabajo se analizan las tensiones interlaminares en laminados simétricos, evaluando la influencia de la orientación de las láminas considerando varias secuencias de apilamiento. Se emplea un modelo basado en diferencias finitas, que se contrasta con resultados de la literatura cientifica obtenidos mediante una metodología diferente [3].
12
13
==2 Descripción del modelo==
14
15
Se ha analizado el problema de un laminado simétrico respecto del plano ‘xy’ compuesto por 4 láminas con secuencia de apilamiento ±ðœ sometido a un estado de carga uniaxial en sus bordes, Figura 1. El laminado tiene un ancho <math display="inline">2b</math>, con <math display="inline">b=</math><math>8{h}_{0}</math>, siendo <math display="inline">{h}_{0}</math> el espesor de una lámina.
16
[[File:Solis_et_al_2018b-image5.png|centre|600x600px]]
17
18
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
19
<span style="text-align: center; font-size: 75%;">'''Figura 1. '''Sistema de ejes geométricos para el laminado y relación con los ejes principales de ortotropía de cada lámina. El plano z=0 corresponde al plano de simetría del laminado.</span></div>
20
21
El laminado está constituido por láminas con comportamiento lineal elástico con simetría ortótropa. La ecuación constitutiva de la lámina en los ejes x- y- z se puede expresar como:
22
{| class="formulaSCP" style="width: 100%; text-align: center;" 
23
|-
24
| <math>\left( \begin{matrix}\begin{matrix}{\sigma }_{x}\\{\sigma }_{y}\\{\sigma }_{z}\end{matrix}\\\begin{matrix}{\tau }_{yz}\\{\tau }_{xz}\\{\tau }_{xy}\end{matrix}\end{matrix}\right) =</math><math>\left( \begin{matrix}\begin{matrix}\overline{{C}_{11}}&\overline{{C}_{12}}&\overline{{C}_{13}}\\\overline{{C}_{12}}&\overline{{C}_{22}}&\overline{{C}_{23}}\\\overline{{C}_{13}}&\overline{{C}_{23}}&\overline{{C}_{33}}\end{matrix}&\begin{matrix}0&0&\overline{{C}_{16}}\\0&0&\overline{{C}_{26}}\\0&0&\overline{{C}_{36}}\end{matrix}\\\begin{matrix}0&0&0\\0&0&0\\\overline{{C}_{16}}&\overline{{C}_{26}}&\overline{{C}_{36}}\end{matrix}&\begin{matrix}\overline{{C}_{44}}&\overline{{C}_{45}}&0\\\overline{{C}_{45}}&\overline{{C}_{55}}&0\\0&0&\overline{{C}_{66}}\end{matrix}\end{matrix}\right) \left( \begin{matrix}\begin{matrix}{\epsilon }_{x}\\{\epsilon }_{y}\\{\epsilon }_{z}\end{matrix}\\\begin{matrix}{\gamma }_{yz}\\{\gamma }_{xz}\\{\gamma }_{xy}\end{matrix}\end{matrix}\right) \quad \quad \quad \quad \, \, (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 1)</math>
25
|}
26
27
28
Las componentes de la matriz de rigidez son función de las constantes elásticas del material y del ángulo de orientación 𝜠de las direcciones de ortotropía.
29
30
En un estado de carga uniaxial, el equilibrio de tensiones predice que las componentes del tensor de tensiones son independientes de la coordenada ‘x’ en puntos suficientemente alejados de la zona de aplicación de la carga (principio de Saint-Venant).
31
32
La independencia de las tensiones con respecto a ‘x’ se traduce en independencia también con respecto a ‘x’ de las componentes del tensor de deformaciones. Asumiendo pequeños desplazamientos y pequeñas deformaciones las relaciones de compatibilidad se escriben:
33
{| class="formulaSCP" style="width: 100%; text-align: center;" 
34
|-
35
| <math>\begin{matrix}{\epsilon }_{x}={u}_{,x}\\{\epsilon }_{y}={v}_{,y}\\\begin{matrix}{\epsilon }_{y}={v}_{,y}\\{\gamma }_{yz}={w}_{,y}+{v}_{,z}\\\begin{matrix}{\gamma }_{xz}={u}_{,z}+{w}_{,x}\\{\gamma }_{xy}={u}_{,y}+{v}_{,x}\end{matrix}\end{matrix}\end{matrix}\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 2)</math>
36
|}
37
38
39
Dado que existe simetría de cargas y geométrica en el laminado bajo estudio, el campo de desplazamiento que resulta de la integración de la ecuación (2) deberá tener la siguiente forma funcional:
40
41
<math display="inline">\begin{matrix}u={\epsilon }_{x}x+U(y,z)\\v=V(y,z)\\w=W(y,z)\end{matrix}\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \, \, (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 3</math>)
42
43
Donde &#x03b5;<sub>x</sub> es la deformación impuesta sobre el laminado bajo estudio. La forma del campo de desplazamientos dada en la ecuación (3) proporciona así un campo de tensiones tal que:
44
{| class="formulaSCP" style="width: 100%; text-align: center;" 
45
|-
46
| <math>{\sigma }_{ij}={\sigma }_{ij}\left( U,V,W\right) ={\sigma }_{ij}\left( y,z\right) \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 4)</math>
47
|}
48
49
50
Lo que muestra que las componentes de tensiones y deformaciones de este problema dependen sólo de las coordenadas ‘y‘ y ‘z‘. Con estas consideraciones las ecuaciones de equilibrio interno se reducen a:
51
{| class="formulaSCP" style="width: 100%; text-align: center;" 
52
|-
53
| <math>\begin{matrix}\frac{\partial {\tau }_{xy}}{\partial y}+\frac{\partial {\tau }_{xz}}{\partial z}=0\\\frac{\partial {\sigma }_{y}}{\partial y}+\frac{\partial {\tau }_{yz}}{\partial z}=0\\\frac{\partial {\tau }_{yz}}{\partial y}+\frac{\partial {\sigma }_{z}}{\partial z}=0\end{matrix}\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \, (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 5)</math>
54
|}
55
56
57
Por tanto, es posible estudiar sólo la parte del laminado contenida en el primer cuadrante (parte izquierda de la Figura 2):
58
[[File:Solis_et_al_2018b-image6.png|centre|600x600px]]
59
''' '''
60
61
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
62
<span style="text-align: center; font-size: 75%;">'''Figura 2. '''Corte transversal del laminado en el plano ‘yz’.</span></div>
63
64
En los bordes del laminado, deberán verificarse las condiciones de contorno de bordes libres de tensiones:
65
{| class="formulaSCP" style="width: 100%; text-align: center;" 
66
|-
67
| <math>\begin{matrix}{\sigma }_{y}={\tau }_{yz}={\tau }_{xy}=0\, \, en\, \, y=b\\{\sigma }_{z}={\tau }_{yz}={\tau }_{xz}=0\, \, en\, \, z=2{h}_{o}\end{matrix}\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 6)</math>
68
|}
69
70
71
Para resolver el problema, se ha considerado útil definir las siguientes variables adimensionales:
72
73
<math display="inline">\overline{y}=\frac{y}{b};\quad \overline{z}=\frac{z}{{h}_{o}};\quad \overline{U}=</math><math>\frac{U}{{h}_{o}};\quad \overline{V}=\frac{V}{b};\quad \overline{W}=\frac{W}{{h}_{o}};\quad R=</math><math>\frac{b}{{h}_{o}}\quad \quad \quad \quad (\, \, SEQ\, Eq\, \ast \, MERGEFORMAT\, 7)</math>
74
75
Combinando las expresiones anteriores es posible formular las ecuaciones que deben verificar las componentes del campo de desplazamientos en cada punto del dominio, ecuacion 8
76
{| class="formulaSCP" style="width: 100%; text-align: center;" 
77
|-
78
| <math>\begin{matrix}R\overline{{C}_{26}}{\overline{V}}_{,\overline{y}\overline{y}}+R\left( \overline{{C}_{36}}+\overline{{C}_{45}}\right) {\overline{W}}_{,\overline{y}\overline{z}}+\overline{{C}_{66}}{\overline{U}}_{,\overline{y}\overline{y}}+{R}^{3}\overline{{C}_{45}}{\overline{V}}_{,\overline{z}\overline{z}}+{R}^{2}\overline{{C}_{55}}{\overline{U}}_{,\overline{z}\overline{z}}=0\\R\overline{{C}_{22}}{\overline{V}}_{,\overline{y}\overline{y}}+R\left( \overline{{C}_{23}}+\overline{{C}_{44}}\right) {\overline{W}}_{,\overline{y}\overline{z}}+\overline{{C}_{26}}{\overline{U}}_{,\overline{y}\overline{y}}+{R}^{3}\overline{{C}_{44}}{\overline{V}}_{,\overline{z}\overline{z}}+{R}^{2}\overline{{C}_{45}}{\overline{U}}_{,\overline{z}\overline{z}}=0\\\overline{{C}_{44}}{\overline{W}}_{,\overline{y}\overline{y}}+R\left( \overline{{C}_{45}}+\overline{{C}_{36}}\right) {\overline{U}}_{,\overline{y}\overline{z}}+{R}^{2}\left( \overline{{C}_{23}}+\overline{{C}_{44}}\right) {\overline{V}}_{,\overline{y}\overline{z}}+{R}^{2}\overline{{C}_{33}}{\overline{W}}_{,\overline{z}\overline{z}}=0\end{matrix}\quad \quad (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 8)\,</math> 
79
|}
80
81
82
Mientras que en los puntos del contorno se deberán verificar las condiciones dadas por las ecuaciones 9 y 10:
83
{| class="formulaSCP" style="width: 100%; text-align: center;" 
84
|-
85
| <math>\begin{matrix}\, \overline{{C}_{22}}{\overline{V}}_{,\overline{y}}+\overline{{C}_{23}}{\overline{W}}_{,\overline{z}}+\frac{1}{R}\overline{{C}_{26}}{\overline{U}}_{,\overline{y}}=-\overline{{C}_{12}}{\epsilon }_{x}\\\frac{1}{R}\overline{{C}_{44}}{\overline{W}}_{,\overline{y}}+R\overline{{C}_{44}}{\overline{V}}_{,\overline{z}}+\overline{{C}_{45}}{\overline{U}}_{,\overline{z}}=0\\\overline{{C}_{26}}{\overline{V}}_{,\overline{y}}+\overline{{C}_{36}}{\overline{W}}_{,\overline{z}}+\frac{1}{R}\overline{{C}_{66}}{\overline{U}}_{,\overline{y}}=-\overline{{C}_{16}}{\epsilon }_{x}\end{matrix}\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 9)\,</math> 
86
|}
87
{|class="formulaSCP" style="width: 100%; text-align: center;" 
88
|-
89
| <math>\, \begin{matrix}\, \overline{{C}_{23}}{\overline{V}}_{,\overline{y}}+\overline{{C}_{33}}{\overline{W}}_{,\overline{z}}+\frac{1}{R}\overline{{C}_{36}}{\overline{U}}_{,\overline{y}}=-\overline{{C}_{13}}{\epsilon }_{x}\\\frac{1}{R}\overline{{C}_{44}}{\overline{W}}_{,\overline{y}}+R\overline{{C}_{44}}{\overline{V}}_{,\overline{z}}+\overline{{C}_{45}}{\overline{U}}_{,\overline{z}}=0\\\frac{1}{R}\overline{{C}_{45}}{\overline{W}}_{,\overline{y}}+R\overline{{C}_{45}}{\overline{V}}_{,\overline{z}}+\overline{{C}_{55}}{\overline{U}}_{,\overline{z}}=0\end{matrix}\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \, (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 10)\,</math> 
90
|}
91
92
93
Para resolver el sistema descrito se emplea un esquema de diferencias finitas que discretizando el dominio mediante un esquema de puntos discretos como el mostrado en la parte derecha de la Figura 2, permita obtener el campo de desplazamientos, deformaciones y tensiones del problema tridimensional completo. Para preservar la continuidad de los desplazamientos y de las tensiones fuera del plano entre dos láminas orientadas con diferentes orientaciones, la discretización se hace de forma que la interfase entre las láminas quede confinada entre dos filas de puntos materiales.
94
95
Para construir el esquema de diferencias finitas se utilizan fórmulas de diferencias centradas para las derivadas segundas y fórmulas de diferencias centradas y regresivas para las derivadas primeras. Asímismo se tienen en cuenta para escribir las ecuaciones de campo, las siguientes relaciones de simetría:
96
{| class="formulaSCP" style="width: 100%; text-align: center;" 
97
|-
98
| <math>\begin{matrix}{U}_{,z}\left( y,0\right) =0\\{V}_{,z}\left( y,0\right) =0\\W\left( y,0\right) =0\end{matrix}\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \, (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 11)\, \,</math> 
99
|}
100
{|class="formulaSCP" style="width: 100%; text-align: center;" 
101
|-
102
| <math>\begin{matrix}{U}_{,y}\left( 0,z\right) =0\\V\left( 0,z\right) =0\\{W}_{,y}\left( 0,z\right) =0\end{matrix}\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \, \, (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 12)\, \,</math> 
103
|}
104
105
106
También para impedir los desplazamientos como sólido rígido se imponen las siguientes condiciones:
107
108
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
109
{| class="formulaSCP" style="width: 100%; text-align: center;" 
110
|-
111
| <math>\begin{matrix}U\left( 0,z\right) =U\left( 0,0\right) =0\\\, V\left( 0,0\right) =0\, \, \\W\left( 0,0\right) =0\end{matrix}\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \, \, (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 13)</math>
112
|}
113
</div>
114
115
Así con todos estos elementos se construye un sistema algebraico del tipo <math display="inline">Ax=</math><math>b</math>. De esta manera, para un dominio discretizado con N puntos, se tiene un sistema con los siguientes elementos:
116
117
Un vector de incógnitas que representa las componentes de los desplazamientos en los puntos del dominio:
118
{| class="formulaSCP" style="width: 100%; text-align: center;" 
119
|-
120
| <math>x=\left( \begin{matrix}\begin{matrix}{U}_{1}\\{V}_{1}\\{W}_{1}\end{matrix}\\\begin{matrix}\begin{matrix}{U}_{2}\\{V}_{2}\\{W}_{2}\end{matrix}\\\begin{matrix}\ldots \\\begin{matrix}{U}_{N}\\{V}_{N}\\{W}_{N}\end{matrix}\end{matrix}\end{matrix}\end{matrix}\right) \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 14)</math>
121
|}
122
123
124
Una matriz que contiene los coeficientes que multiplican a cada una de las incógnitas y que están expuestos en las ecuaciones (8), (9) y (10).
125
126
Y un vector que contiene los términos independientes, descritos nuevamente en las ecuaciones (8), (9) y (10).
127
128
Conviene notar que la estructura del sistema algebraico de ecuaciones a resolver se puede escribir también:
129
{| class="formulaSCP" style="width: 100%; text-align: center;" 
130
|-
131
| <math>A\left( \frac{x}{{\epsilon }_{x}}\right) =\left( \frac{b}{{\epsilon }_{x}}\right) \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \, \, (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 15)</math>
132
|}
133
134
135
Dada la estructura del vector ''b''   se tiene que:
136
{| class="formulaSCP" style="width: 100%; text-align: center;" 
137
|-
138
| <math>\frac{b}{{\epsilon }_{x}}\not =f\left( {\epsilon }_{x}\right) \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad (\, SEQ\, Eq\, \ast \, MERGEFORMAT\, 16)</math>
139
|}
140
141
142
Esto implica que es posible obtener los desplazamientos de los puntos de la malla por unidad de deformación axial y obtener, gracias a la linealidad de los operadores involucrados, y a partir de éstos, las deformaciones y las tensiones también por unidad de deformación. Para obtener el valor real de las variables, para un nivel de deformación aplicado, basta por tanto multiplicar los resultados obtenidos por dicho valor.
143
144
Se ha realizado un análisis de sensibilidad de los resultados del modelo al tamaño de malla. Para ello se ha estudiado un laminado con una relación geométrica b=8h<sub>o</sub>. Se ha seleccionado un laminado con dos secuencias de apilamiento [-45/45]<sub>S</sub> y [0/-45/90/45]<sub>S</sub> de un material cuyas propiedades mecánicas se han tomado de la literatura [3]. Se han considerado cuatro mallas de 184, 363, 616 y 864 puntos.
145
146
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
147
<span style="text-align: center; font-size: 75%;">'''Tabla 1. '''Propiedades mecánicas del material utilizado en el análisis de sensibilidad [3]</span></div>
148
149
{| style="width: 81%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
150
|-
151
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">E<sub>11</sub></span>
152
153
<span style="text-align: center; font-size: 75%;">(GPa)</span>
154
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">E<sub>22</sub></span>
155
156
<span style="text-align: center; font-size: 75%;">E<sub>33</sub></span>
157
158
<span style="text-align: center; font-size: 75%;">(GPa)</span>
159
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">G<sub>12</sub></span>
160
161
<span style="text-align: center; font-size: 75%;">(GPa)</span>
162
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">G<sub>13</sub></span>
163
164
<span style="text-align: center; font-size: 75%;">G<sub>23</sub></span>
165
166
<span style="text-align: center; font-size: 75%;">(GPa)</span>
167
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">&#x03bd;<sub>12</sub></span>
168
169
<span style="text-align: center; font-size: 75%;">&#x03bd;<sub>13</sub></span>
170
171
<span style="text-align: center; font-size: 75%;">&#x03bd;<sub>23</sub></span>
172
|-
173
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">137.9</span>
174
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">14.5</span>
175
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">5.86</span>
176
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">5.86</span>
177
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.21</span>
178
|}
179
180
181
Se han analizado las tensiones fuera del plano del laminado a lo largo del ancho de la interfase entre cada lámina y la distribución de estas tensiones a lo largo del espesor en el borde libre para los dos laminados analizados. En las Figuras 3 y 4 se muestran los resultados correspondientes al laminado <span style="text-align: center; font-size: 75%;">[-45/45]<sub>S</sub></span>, para el otro laminado los resultados son semejantes. Como era de esperar, las tensiones fuera del plano son prácticamente nulas durante alrededor de un 50% del ancho del laminado para el caso representado en la Figura 3. La distancia al borde libre en la que las tensión cae a cero es igual al espesor del laminado [5,6]. La tensión <math display="inline">{\sigma }_{xz}\,</math> crece al acercarse a la interfase entre dos láminas con orientaciones diferentes, figura 4a, mientras que la tensión <math display="inline">{\sigma }_{zz}</math>cambia de signo conforme se acerca a la interfase, figura 4b.
182
183
La variación de las tensiones es independiente del número de puntos de la malla, recogiendose con más precisión los incrementos de las tensiones al aumentar el número de puntos.
184
185
Los resultados muestran que considerar una malla de 616 puntos es suficiente dado que el incremento en la precisión que se consigue con una malla más fina es prácticamente despreciable.
186
[[File:Solis_et_al_2018b-image7-c.png|centre|561x561px]]
187
188
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
189
<span style="text-align: center; font-size: 75%;">a)</span></div>
190
191
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
192
 [[Image:Solis_et_al_2018b-image8-c.png|561x561px]] </div>
193
194
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
195
<span style="text-align: center; font-size: 75%;">b)</span></div>
196
197
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
198
<span style="text-align: center; font-size: 75%;">'''Figura 3. '''Tensiones fuera del plano a través del ancho en la interfase para distintas mallas en el laminado [-45/45]<sub>S</sub>, a) tensión </span> <math display="inline">{\sigma }_{xz}</math><span style="text-align: center; font-size: 75%;"> , b) tensión </span> <math display="inline">{\sigma }_{zz}</math></div>
199
200
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
201
 [[Image:Solis_et_al_2018b-image9-c.png|600px]] </div>
202
203
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
204
<span style="text-align: center; font-size: 75%;">a)</span></div>
205
206
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
207
 [[Image:Solis_et_al_2018b-image10-c.png|600px]] </div>
208
209
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
210
<span style="text-align: center; font-size: 75%;">b)</span></div>
211
212
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
213
<span style="text-align: center; font-size: 75%;">'''Figura 4. '''Tensiones fuera del plano a través del espesor en el borde libre para distintas mallas en el laminado [-45/45]<sub>S</sub> a) tensión </span> <math display="inline">{\sigma }_{xz}</math><span style="text-align: center; font-size: 75%;"> , b) tensión </span> <math display="inline">{\sigma }_{zz}</math></div>
214
215
Con objeto de verificar los resultados que proporciona el esquema empleado para el análisis de las tensiones fuera del plano, se ha realizado una comparación con los mostrados por Lessard et al. [3] para dos secuencias de apilamiento <math display="inline">{[0/90]}_{s}</math> y <math display="inline">{[90/0]}_{s}</math>. Las propiedades del material analizado son las mismas que las mostradas en la Tabla 1. A modo de ejemplo se muestran sólo los resultados del primero de los laminados.
216
217
Se han estudiado las tensiones <math display="inline">{\sigma }_{zz}</math> y <math display="inline">{\sigma }_{yz}</math> a lo largo del ancho del laminado en la interfase de láminas con distinta orientación (Figuras 5 y 6), y a lo largo del espesor en el borde libre (Figuras 7 y 8), y se han comparado con las curvas análogas del artículo de Lessard et al. [3] para una deformación axial aplicada de un 1%.
218
219
Conviene recordar que la tensión <math display="inline">{\sigma }_{yz}</math> a lo largo del espesor en el borde libre se anula por condición de contorno. Así que las gráficas representadas para esta tensión a lo largo del espesor, figura 8, corresponden a los valores que se han estimado en una zona interior del laminado muy cercana al borde libre (en concreto en los nodos inmediatamente anteriores a los situados en el borde libre).
220
[[File:Solis_et_al_2018b-image11-c.png|centre|561x561px]]
221
222
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
223
<span style="text-align: center; font-size: 75%;">'''Figura 5. '''Tensión  </span> <math display="inline">{\sigma }_{zz}</math><span style="text-align: center; font-size: 75%;">a través del ancho en la interfase. Laminado [0/90]<sub>S</sub>''' '''</span></div>
224
[[File:Solis_et_al_2018b-image12.png|centre|561x561px]]
225
226
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
227
<span style="text-align: center; font-size: 75%;">'''Figura 6. '''Tensión &#x03c3;<sub>yz</sub> a través del ancho en la interfase. laminado [0/90]<sub>S</sub></span></div>
228
229
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
230
 [[Image:Solis_et_al_2018b-image13-c.png|600px]] </div>
231
232
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
233
<span style="text-align: center; font-size: 75%;">'''Figura 7. '''Tensión  </span> <math display="inline">{\sigma }_{zz}</math><span style="text-align: center; font-size: 75%;">a través del espesor del laminado. a) laminado [0/90]<sub>S</sub></span></div>
234
[[File:Solis_et_al_2018b-image14-c.png|centre|561x561px]]
235
236
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
237
<span style="text-align: center; font-size: 75%;">'''Figura 8. '''Tensión &#x03c3;<sub>yz</sub> a través del espesor del laminado. Laminado [0/90]<sub>S </sub></span></div>
238
239
Se observa que los resultados obtenidos mediante el modelo desarrollado reproducen las tendencias mostradas en [3] y aproximan razonablemente los valores extremos de las curvas. Las variaciones de las tensiones interlaminares observadas en la figura 5 y 6 también son semejantes a las observadas por otros autores [4 y 9].
240
241
==3 Resultados==
242
243
Con el modelo desarrollado se ha estudiado la variación del valor de las tensiones interlaminares en laminados sometidos a una carga uniaxial en su plano, considerando diferentes secuencias de apilamiento en las que se ha variado la orientación de las láminas. Se ha analizado la respuesta de cuatro secuencias de apilamiento, dos con 4 láminas <span style="text-align: center; font-size: 75%;">[𝛳/−𝛳]𝑠 </span>y <span style="text-align: center; font-size: 75%;">[2𝛳/−𝛳]𝑠</span>, y otros dos con 6 láminas <span style="text-align: center; font-size: 75%;">[𝛳<sub>2</sub>/−𝛳]𝑠 y [𝛳/2𝛳/−𝛳]𝑠.</span> En todos los casos se ha analizado la variación de la orientación entre 0º y 90º. Estas secuencias de apilamiento se han definido de la siguiente forma, por ejemplo para una orientación de &#x03b8;=30º, los laminados considerados serán <span style="text-align: center; font-size: 75%;">[30/−30]𝑠 </span>y <span style="text-align: center; font-size: 75%;">[60/−30]𝑠 </span>los de 4 láminas, y <span style="text-align: center; font-size: 75%;">[30<sub>2</sub>/−30]𝑠 y [30/60/−30]𝑠 </span>los de 6 láminas.
244
245
En las Figuras 9 y 10 se muestran los valores de las diferentes componentes de las tensiones interlaminares en puntos cercanos al borde libre a la altura de la interfase, al variar el ángulo &#x03b8; para los laminados <span style="text-align: center; font-size: 75%;">[𝛳/−𝛳]𝑠 </span>y <span style="text-align: center; font-size: 75%;">[2𝛳/−𝛳]𝑠, </span>respectivamente. Se observa que en el laminado <span style="text-align: center; font-size: 75%;">[𝛳/−𝛳]𝑠 </span>la componente predominante para todas las orientaciones &#x03b8; es &#x03c3;<sub>xz</sub>, tomando un valor máximo para &#x03b8;=15.5º. La componente &#x03c3;<sub>zz</sub> tiene valores un orden de magnitud inferiores, siendo máximo para &#x03b8;=32º.Dado que esta componente es de compresión en todo el rango de <span style="text-align: center; font-size: 75%;">𝛳, </span>no contribuye a la aparición de delaminación. Sin embargo, en el laminado <span style="text-align: center; font-size: 75%;">[2𝛳/−𝛳]𝑠 </span>los valores de las componentes&#x03c3;<sub>zz</sub> y &#x03c3;<sub>xz </sub>tienen valores similares entre 40º y 90º, siendo los valores de la componente &#x03c3;<sub>xz</sub> mucho mayores para &#x03b8; entre 0 y 40º. Entre 15º y 60º, la componente &#x03c3;<sub>zz</sub> tiene valores negativos, no contribuyendo por tanto al modo de fallo de delaminación. La componente &#x03c3;<sub>xz </sub>se hace máxima para un ángulo en torno a&#x03b8;=10.1º, con valores similares a los del laminado <span style="text-align: center; font-size: 75%;">[𝛳/−𝛳]𝑠.</span>
246
[[File:Solis_et_al_2018b-image15-c.png|centre|561x561px]]
247
248
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
249
<span style="text-align: center; font-size: 75%;">'''Figura 9. '''Tensiones interlaminares en el borde libre a la altura de la interfase para distintas orientaciones de las fibras para un laminado tipo [𝛳/−𝛳]𝑠.</span></div>
250
[[File:Solis_et_al_2018b-image16-c.png|centre|561x561px]]
251
252
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
253
<span style="text-align: center; font-size: 75%;">'''Figura 10. '''Tensiones interlaminares en el borde libre a la altura de la interfase para distintas orientaciones de las fibras para un laminado tipo </span> <math display="inline">{[2\theta ,-\theta ]}_{s}</math></div>
254
255
En las Figuras 11 y 12 se muestran los valores de las diferentes componentes de las tensiones interlaminares en el borde libre a la mitad del espesor desde el plano medio, al variar el ángulo &#x03b8; para los laminados <span style="text-align: center; font-size: 75%;">[𝛳<sub>2</sub>/−𝛳]𝑠 y [𝛳/2𝛳/−𝛳]𝑠, </span>respectivamente. Las tensiones se han obtenido a esa altura para comparar el mismo lugar geométrico en todos los laminados analizados.
256
257
En la Figura 11 se observa que en el laminado <span style="text-align: center; font-size: 75%;">[𝛳<sub>2</sub>/−𝛳]𝑠 </span>la componente&#x03c3;<sub>zz </sub>es prácticamentenula para cualquier orientación. La componente &#x03c3;<sub>xz </sub>es máxima para una orientación &#x03b8;=16º, pero con valores mucho menores de los observados en los láminados de 4 láminas.
258
259
En la Figura 12 se observa la variación de las componentes de las tensiones interlaminares para el laminado <span style="text-align: center; font-size: 75%;">[𝛳/2𝛳/−𝛳]𝑠. </span>En este laminado, la componente&#x03c3;<sub>zz </sub>solo se anula para orientacionesde 20º y 60º, teniendo valores positivos para ángulos entre 20º y 60º, contribuyendo para estas orientaciones a ladelaminación del laminado. Entre 60º y 90º esta componente es mayor que la componente&#x03c3;<sub>xz</sub>. La componente &#x03c3;<sub>xz</sub> tiene valores similares a los obtenidos en el laminado <span style="text-align: center; font-size: 75%;">[𝛳<sub>2</sub>/−𝛳]𝑠 </span>para los angulos entre 0 y 60, siendo máximo en este caso para 19º.
260
[[File:Solis_et_al_2018b-image17-c.png|centre|561x561px]]
261
262
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
263
<span style="text-align: center; font-size: 75%;">'''Figura 11. '''Tensiones interlaminares en el borde libre a la mitad del espesor desde el plano medio del laminado para distintas orientaciones de las fibras para un laminado tipo </span> <math display="inline">{[{\theta }_{2},-\theta ]}_{s}</math></div>
264
[[File:Solis_et_al_2018b-image18-c.png|centre|561x561px]]
265
266
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
267
<span style="text-align: center; font-size: 75%;">'''Figura 12. '''Tensiones intralaminares en el borde libre a la mitad del espesor desde el plano medio para distintas orientaciones de las fibras para un laminado tipo </span> <math display="inline">{[\theta ,2\theta ,-\theta ]}_{s}</math></div>
268
269
Observamos como la tensión <math display="inline">{\sigma }_{yz}</math> es nula por condición de contorno en todos los laminados analizados.
270
271
==4 Conclusiones==
272
273
Se han analizado las tensiones interlaminares que aparecen debido al efecto de borde libre en laminados simétricos sometidos a cargas uniaxiales en su plano. Se ha puesto de manifiesto que la secuencia de apilamiento no sólo afecta al valor de las tensiones interlaminares sino también a su valor relativo.
274
275
En laminados con el mismo número de láminas orientadas a +&#x03b8; y -&#x03b8; la tensión predominante es &#x03c3;<sub>xz</sub> para todas las orientaciones. Este comportamiento también se produce cuando se aumenta el número de láminas orientadas a +&#x03b8;. Sin embargo si la secuencia de apilamiento es del tipo <math display="inline">{[2\theta ,-\theta ]}_{s}</math><span style="text-align: center; font-size: 75%;"> o </span> <math display="inline">{[\theta ,2\theta ,-\theta ]}_{s}</math>existe un rango de orientaciones en la que la componente&#x03c3;<sub>zz</sub> es mayor que &#x03c3;<sub>xz</sub>. Aunque el valor máximo de la componente normal sea inferior a la de cortadura.
276
277
==Agradecimientos==
278
279
Los autores agradecen la financiación recibida para el desarrollo de este trabajo del Ministerio de Economía y Finanzas de España en el marco del proyecto DPI2013-42240-R.
280
281
==Referencias==
282
283
[1] C.T. Sun, G.D. Chu, ''Journal of Composite Materials,'' '''25''', pág. 142-161 (1991).
284
285
[2] S. Sihn, R.Y. Kim, K. Kawabe, S.W. Tsai, ''Composites Science and Technology, '' '''67''', pág. 996–1008 (2007).
286
287
[3] L.B. Lessard, A.S. Schmidt, ''International Journal Solids Structures, '' '''33''', 15, pág. 2243-2259 (1996).
288
289
[4] V-T. Nguyen, J.F. Caron, ''Computers and Structures, '' '''84''', pág. 1538–1546 (2006).
290
291
[5] J. Lindemann, W. Becker, ''Mechanics of Composite Materials,'' '''36''', 3, pág. 208-214 (2000).
292
293
[6] R.B. Pipes, N.J. Pagano, ''Journal of Composite Materials,'' '''4''', pág. 538-548 (1970).
294
295
[7] T. Kant, K. Swaminathan, ''Composite Structures,'' '''49''', pág. 65–75 (2000).
296
297
[8] C. Mittelstedt, W. Becker, ''Journal of Composite Materials''; '''38''', 12, pág. 1037-1062 (2004).
298
299
[9] A.S.D. Wang, F.W. Crossman, ''Journal of Composite Materials'', '''11''', pág. 92-106 (1977).
300

Return to Solis et al 2018b.

Back to Top
GET PDF

Document information

Published on 14/10/18
Accepted on 14/10/18
Submitted on 14/10/18

Volume 02 - Comunicaciones Matcomp17 (2018), Issue Núm. 4 - Comportamiento en servicio de los materiales compuestos (2), 2018
DOI: 10.23967/r.matcomp.2018.10.003
Licence: Other

Document Score

0

Views 2
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?