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
2
== Resumen ==
3
4
El presente trabajo presenta los conceptos básicos y de aplicación de la estrategia que combina los métodos de elementos finitos (FEM) y de elementos discretos (DEM) para el estudio de la propagación de fracturas en estructuras de hormigón. El cálculo de la estructura, modelada como un continuo, se inicia con el FEM y se hace uso del DEM para iniciar y hacer crecer las grietas que puedan aparecer en la estructura. Esta metodología ha sido propuesta por los autores en dos y tres dimensiones. Recientemente se ha agregado el uso de elementos unidimensionales de acero embebidos en el continuo para modelar el efecto resistente del armado en estructuras de hormigón. En el trabajo se presentan diferentes ejemplos de aplicación al estudio de la rotura de piezas y estructuras de hormigón en masa y armado, así como la rotura múltiple de una estructura histórica de mampostería debida a un terremoto.
5
6
'''Keywords''': Método de elementos discretos, Método de elementos finitos, Mecánica de fractura, Estrategia FEM-DEM
7
8
==1.- Introducción==
9
10
Recientemente muchos autores han realizado interesantes desarrollos con la finalidad de poder definir el inicio y crecimiento de una fractura dentro de una estructura, modelada como un medio continuo, ya sea frágil o dúctil<div id="citeF-1"></div><div id="citeF-2"></div><div id="citef-3"></div><div id="citeF-4"></div>[[#cite-1|[1-]][[#cite-4|4]]].
11
12
Uno de los trabajos más recientes discretiza el medio continuo por medio de elementos discretos de forma circular (en 2 dimensiones) o esférica (en 3 dimensiones)<div id="citeF-5"></div>[[#cite-5|[5]]]. Sin embargo, la dificultad inherente para calibrar los parámetros del material en el método de los elementos discretos (DEM, pos sus siglas en inglés), así como la necesidad de contar con un gran número de elementos discretos<div id="citeF-2"></div><div id="citeF-5"></div>[[#cite-2|[2]],[[#cite-5|5]]], pone en duda su efectividad, a pesar de que cualitativamente el número y dirección de las grietas y los resultados numéricos obtenidos de la carga última de la estructura son bastante aceptables.
13
14
Este trabajo hace uso de una estrategia que utiliza el método de elementos finitos (FEM, por sus siglas en inglés) para discretizar la estructura inicial y seguir su comportamiento bajo cargas crecientes hasta el inicio de la primera fisura. Tras ello se utiliza una técnica de eliminación de los elementos finitos dañados y se introducen elementos discretos en los labios de las fisuras <div id="citeF-6"></div><div id="citeF-7"></div><div id="citeF-8"></div>[[#cite-6|[6-]][[#cite-8|8]]]. La denominada técnica FEM-DEM se ha utilizado con éxito en 2 y 3 dimensiones <div id="citeF-9"></div><div id="citeF-10"></div>[[#cite-9|[9]],[[#cite-10|10]]] para estudiar el comportamiento no lineal y el fallo de estructuras de hormigón y mampostería bajo diferentes solicitaciones cuasi-estáticas y dinámicas.
15
16
Otros aspectos importantes inherentes a la formulación FEM-DEM son el uso de un campo de esfuerzos suavizado, la conservación de la masa y el uso de un algoritmo simple para asegurar el contacto post-fractura entre las paredes de la grieta. Adicionalmente, la armadura de acero se considera embebida en la malla de elementos finitos, lo cual simplifica los cálculos.
17
18
==2.- Formulación FEM-DEM==
19
20
En general, la estrategia que se sigue para el cálculo no lineal de estructuras con modelos de elementos finitos de sólido bidimensionales (2D) o tridimensionales (3D)<div id="citeF-11"></div>[[#cite-11|[11]]] pasa para evaluar a lo largo del tiempo la respuesta de la estructura discretizada con una malla de elementos finitos. En dicha malla las conectividades inter-elementales emulan los enlaces entre los elementos que se degradan de forma progresiva, generalmente mediante un sencillo modelo de daño isótropo. Dichos enlaces interelementales pueden interpretarse como las conexiones entre una colección de elementos discretos que reemplazarían los nodos de la malla de elementos finitos. Esta analogía es el punto de partida de la técnica FEM-DEM<div id="citeF-9"></div><div id="citeF-10"></div>[[#cite-9|[9]],[[#cite-10|10]]].
21
22
La estrategia FEM-DEM se puede resumir en los 5 pasos siguientes:
23
24
:*  Discretización del continuo que modela la estructura mediante el FEM.
25
26
:*  Obtención del campo de tensiones sobre la estructura.
27
28
:*  Obtención del daño en el interior y los lados de los elementos.
29
30
:*  Discretización mediante el DEM.
31
32
:*  Integración temporal en sub-pasos.
33
34
Debido a la naturaleza de los ejemplos, se describirán los pasos anteriores enfocados al caso 3D utilizando tetraedros lineales de cuatro nodos para la discretización con el FEM.
35
36
Una de las claves de la técnica FEM-DEM es el procedimiento para pasar el daño producido en un elemento finito, con la consiguiente degradación de la rigidez elemental, a una discretización de los labios de la fisura por elementos discretos de forma circular (en 2D) o esférica (en 3D).
37
38
Cuando el daño que se induce en un elemento es mayor que un cierto valor, se considera que la rigidez del elemento se encuentra tan disminuida que es posible eliminarlo, para ello se utiliza una técnica de eliminación de elementos<div id="citeF-6"></div><div id="citeF-7"></div><div id="citeF-8"></div>[[#cite-6|[6-]][[#cite-8|8]]]. En ese momento, se crean nuevos elementos discretos en los vértices del tetraedro eliminado (en 3D), lo que permite la apertura de la grieta en el continuo discretizado por los elementos finitos, mientras que los labios de la fisura quedandefinidos por los elementos discretos. A medida que la fisura crece, e incluso se ramifica, algunos elementos discretos pueden separarse de la malla de elementos finitos creando una disgregación del continuo. El hecho de utilizar elementos discretos para definir las grietas permite, de manera natural, considerar la apertura y cierre de estas sin añadir procedimientos adicionales.
39
40
==3.- Discretización mediante el FEM==
41
42
El FEM es probablemente la técnica numérica más popular para modelar el inicio y la propagación de fisuras en materiales friccionales (hormigón, rocas, mampostería, cerámicas, etc.) <div id="citeF-12"></div><div id="citeF-13"></div><div id="citeF-14"></div><div id="citeF-15"></div><div id="citeF-16"></div>[[#cite-12|[12-]][[#cite-16|16]]]. Sin embargo, la mayoría de los procedimientos basados en el FEM para la predicción de la aparición y evolución de grietas utilizan formulaciones de elementos muy sofisticadas, que en ocasiones requieren un remallado en la vecindad de las posibles grietas <div id="citeF-14"></div><div id="citeF-15"></div>[[#cite-15|15]],[[#cite-16|16]]]. El enfoque seguido en este trabajo utiliza el FEM para modelar la estructura mediante elementos de sólido (2D o 3D), cuya eventual fractura se describe mediante el DEM. En trabajos previos de la técnica FEM-DEM hemos utilizado el sencillo triángulo de tres nodos en 2D y el tetraedro de cuatro nodos en 3D<div id="citeF-9"></div><div id="citeF-10"></div>[[#cite-9|[9]],[[#cite-10|10]]].
43
44
Considerando, por ejemplo, el caso 3D, inicialmente todo el dominio de la estructura se discretiza con una malla de tetraedros, como en cualquier análisis lineal 3D por el FEM. La no linealidad se modela introduciendo una degradación de la rigidez elemental mediante la clásica teoría de daño isótropo <div id="citeF-13"></div><div id="citeF-17">[[#cite-13|13]],[[#cite-17|17]]].
45
46
El efecto resistente de la armadura de acero se modela mediante elementos unidimensionales (1D) embebidos en la malla de tetraedros. Ello implica localizar los puntos de intersección entre los elementos 1D que modelan el acero y todos los tetraedros que discretizan el dominio. Para cada tetraedro intersectado se determinan los puntos de intersección, los cuales corresponden a los nodos del segmento de la armadura embebida en el tetraedro. Como es posible describir el desplazamiento de dichosnodos en función de los nodos del tetraedro, la matriz de rigidez del elemento 1D de acero es función de dichos desplazamientos.
47
48
Esta técnica no incrementa el número de ecuaciones a resolver y permite de una forma muy sencilla tomar en cuenta la rigidez asociada a la armadura.
49
50
===3.1 Definición del campo de tensiones===
51
52
El considerar un campo de tensiones adecuado es fundamental para poder definir correctamente el nivel del daño en la estructura y la trayectoria de las fisuras. Dado que en el FEM las tensiones entre elementos son discontinuas <div id="citeF-11"></div>[[#cite-11|[11]]], es esencial realizar un alisado de las mismas para obtener un campo continuo de las tensiones entre elementos. En problemas 2D, se obtiene un campo de tensiones continuo alisando los valores de las tensiones, calculadas inicialmente en el baricentro de cada elemento, para obtener su valor el centro de los lados de cada elemento. Esta estrategia corresponde al método de la parcela súper-convergente propuesta por Zienkiewicz, and Zhu <div id="citeF-18"></div>[[#cite-18|[18]]] y evita agregar términos de estabilización al campo de tensiones como es necesario en otros procedimientos alternativos <div id="citeF-12"></div><div id="citeF-13"></div><div id="citeF-14"></div>[[#cite-12|[12-]][[#cite-14|14]]].
53
54
En análisis 3D se procede de manera similar, ya que el campo de tensiones entre elementos es discontinuo, por lo que se utiliza el mismo mecanismo de suavizado. Es decir, las tensiones en cada una de las aristas de un elemento tetraédrico es la media aritmética de las tensiones evaluadas dentro de todos los elementos que comparten dicha arista.
55
56
Al igual que en el caso 2D <div id="citeF-9">[[#cite-9|[9]]], se ha utilizado el criterio de fallo de Rankin para definir el inicio de una fractura en materiales frágiles tomando en cuenta la tensión principal mayor <math display="inline">{\sigma }_{I}</math> <div id="citeF-13"></div><div id="citeF-17"></div>[[#cite-13|[13]],[[#cite-17|[17]]]. El crecimiento del daño se realiza mediante la siguiente ley exponencial que degrada progresivamente la rigidez del material <div id="citeF-19">[[#cite-19|[19]].
57
58
{| class="formulaSCP" style="width: 100%; text-align: center;" 
59
|-
60
| <math>d=1-\frac{{\sigma }_{f}}{{\sigma }_{I}}exp\left[ A\left( 1-\frac{{\sigma }_{I}}{{\sigma }_{f}}\right) \right] \, \, con\quad A=</math><math>{\left[ \frac{G{E}_{0}}{l{\sigma }_{f}^{2}}-\frac{1}{2}\right] }^{-1}(1)</math>
61
|}
62
63
donde <math display="inline">\, l</math><es la longitud de la arista del elemento finito dañado, <math display="inline">{\sigma }_{f}</math> es el límite elástico a tracción del material, <math display="inline">G</math> la energía de fractura y <math display="inline">{E}_{0}</math> el módulo de Young del material no dañado. Resulta obvio que entre dos pasos de tiempo el daño de un lado o una arista no puede disminuir, ya que ello implicaría la reparación espontánea del material, lo cual es termodinámicamente inadmisible.
64
65
===3.2 Definición del daño elemental===
66
67
Una vez que el daño ha sido evaluado en cada lado o arista de un elemento (es decir, un enlace virtual entre elementos discretos) el daño sobre el conjunto del elemento se calcula como el máximo daño existente en todos los planos que cortan un triángulo (2D) o un tetraedro (3D) (Figura 1). En cada tetraedro existen cuatro planos de corte que aíslan un vértice (Figura 1b) y otros tres que aíslan dos vértices (Figura 1c). De esta manera, el daño sobre el plano de corte se define como el valor medio del daño de cada arista del elemento que corta dicho plano.
68
69
Si el daño elemental sobrepasa un cierto umbral entonces se elimina el triángulo o tetraedro de la malla [6-8] y se crean nuevos elementos discretos en los vértices. El tamaño de estos elementos discretos se define de manera que se conserve la masa del elemento eliminado.
70
71
[[Image:Draft_Samper_634227043-image1.jpeg|600px]]
72
73
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
74
<span style="text-align: center; font-size: 75%;">Figura 1. Planos de corte. a) Elemento triangular b) Tetraedro aislando un vértice c) Tetraedro aislando dos vértices.</span></div>
75
76
=5.- Discretización mediante DEM=
77
78
Cuando un tetraedro es totalmente eliminado (''i.e.'' su rigidez es despreciable) se crean cuatro nuevos elementos discretos en los vértices de dicho elemento finito. La creación de dichos elementos discretos queda condicionada a que no hayan sido creados con anterioridad debido a la eliminación de algún elemento finito vecino. En este trabajo se han utilizado círculos y esferas para representar los elementos discretos en 2D y 3D, respectivamente. La masa de cada elemento discreto corresponde a la masa nodal que es compatible con la obtenida mediante elementos finitos. El radio de un nuevo elemento discreto esférico se obtiene de tal manera que sea el máximo que garantice el contacto con las esferas vecinas, pero sin crear solapamientos. Existen otros algoritmos para generar elementos discretos [2, 20], sin embargo, el propuesto, a pesar de su simplicidad, ha dado buenos resultados.
79
80
Una vez que se crea un elemento discreto, las fuerzas en las interfaces de contacto se usan para definir la interacción de dicho elemento con los adyacentes. Estas fuerzas se deben únicamente a un mecanismo de contacto en las direcciones normal y tangencial al plano de contacto entre esferas (en 3D), considerando el radio mínimo entre las dos partículas en contacto [5].
81
82
En problemas 3D, la fuerza normal de contacto se genera en el punto de contacto entre dos esferas y viene dada por.
83
{| class="formulaSCP" style="width: 100%; text-align: center;" 
84
|-
85
| <math>{F}_{n}^{ij}=\frac{{A}^{ij}}{{d}_{ij}}{E}_{0}{u}_{n}\, \, con\quad {A}^{ij}=</math><math>\pi {r}_{c}^{2}\quad \, \, (2)</math>
86
|}
87
88
89
<span style="text-align: center; font-size: 75%;">donde </span> <math display="inline">{A}^{ij}</math><span style="text-align: center; font-size: 75%;"> es el área de la superficie de contacto entre las dos esferas, </span> <math display="inline">{r}_{c}</math><span style="text-align: center; font-size: 75%;">es el radio de la menor de las dos esferas que interactúan en la interfaz ''i j''. </span> <math display="inline">{u}_{n}</math><span style="text-align: center; font-size: 75%;"> es el solapamiento en la dirección de los centros de las dos esferas y </span> <math display="inline">{d}_{ij}</math><span style="text-align: center; font-size: 75%;"> es la distancia entre dichos centros [5].</span>
90
91
<span style="text-align: center; font-size: 75%;">La fuerza tangencial </span> <math display="inline">{F}_{s}</math><span style="text-align: center; font-size: 75%;"> en el punto de contacto se descompone en dos direcciones ortogonales </span> <math display="inline">{s}_{1}</math><span style="text-align: center; font-size: 75%;">y </span> <math display="inline">{s}_{2}</math><span style="text-align: center; font-size: 75%;">contenidas en el plano normal a la dirección </span> <math display="inline">{d}_{ij}</math><span style="text-align: center; font-size: 75%;">. Para cada dirección </span> <math display="inline">{s}_{i}</math><span style="text-align: center; font-size: 75%;"> la fuerza tangencial se describe por</span>
92
{| class="formulaSCP" style="width: 100%; text-align: center;" 
93
|-
94
| <math>{F}_{s}=min\left\{ \begin{matrix}{u}_{{s}_{i}}\frac{{E}_{0}}{2(1+\upsilon )}\\\mu {F}_{n}\left( \frac{{u}_{{s}_{i}}}{{u}_{s}}\right) \end{matrix}\right. \quad \quad (3)</math>
95
|}
96
97
98
<span style="text-align: center; font-size: 75%;">donde </span> <math display="inline">\mu</math> <span style="text-align: center; font-size: 75%;"> es el coeficiente de rozamiento entre esferas y </span> <math display="inline">\upsilon</math> <span style="text-align: center; font-size: 75%;"> es el coeficiente de Poisson. Más detalles de las ecuaciones (2) y (3) se pueden encontrar en [5].</span>
99
100
En general, el número total de elementos discretos generados en un análisis es solo una pequeña fracción del número total de nodos en la malla de elementos finitos [10]. Por lo tanto, los algoritmos de búsqueda para localizar las fuerzas de contacto entre elementos discretos solo representan un pequeño porcentaje del tiempo total de cálculo.
101
102
Es importante mencionar que la estrategia del uso de elementos discretos para definir el contacto entre los labios de una fisura permite resolver problemas estructurales en donde existe apertura y cierre de múltiples grietas, como se muestra en el análisis sísmico de un edificio histórico de mampostería, presentado al final de este artículo.
103
104
=6.- Integración temporal en sub-pasos=
105
106
Uno de los principales problemas en el análisis 3D en el tiempo es el gran número de ecuaciones que surgen de la discretización. Sumando a estos la dificultad para considerar un gran período de tiempo, se concluye que en problemas dinámicos la integración implícita en el tiempo es la mejor estrategia de solución. Sin embargo, este tipo de integraciones son complejas de utilizar en el DEM debido a que es imposible identificar y cuantificar correctamente los contactos y las fuerzas de contacto entreelementos discretos.
107
108
La implementación para el estudio de problemas transitorios seguida en esta investigación corresponde a un esquema de integración temporal de sub-etapas en el cual los elementos finitos se calculan mediante un esquema implícito de Newmark y los elementos discretos mediante un esquema explícito de diferencias centradas. Normalmente el incremento de tiempo del esquema implícito suele ser 100 veces mayor que el explícito.
109
110
<span style="text-align: center; font-size: 75%;">La ventaja de esta estrategia es en que el número de elementos discretos suele ser mucho menor que el de elementos finitos que discretizan la estructura, por lo que la integración explícita realizada dentro de un paso de tiempo implícito es bastante rápida. Antes de comenzar un nuevo paso de tiempo </span> <math display="inline">{\Delta }_{I}t</math><span style="text-align: center; font-size: 75%;"> en el esquema implícito, se calcula el mismo período de tiempo con un esquema explícito sobre los elementos discretos usando un incremento de tiempo  </span> <math display="inline">{\Delta }_{E}t</math><span style="text-align: center; font-size: 75%;">.</span>
111
112
El contacto entre elementos discretos se cuantifica por la suma de los impulsos que se producen a lo largo del análisis explícito y se expresa como una fuerza sobre los elementos discretos, aplicada en el tiempo ''t'' + 1 correspondiente al esquema implícito tal como se muestra en la Figura 2.
113
114
[[Image:Draft_Samper_634227043-image2.jpeg|600px]]
115
116
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
117
<span style="text-align: center; font-size: 75%;">Figura2 Integración temporal en sub-etapas utilizado para el análisis conjunto de elementos finitos y discretos.</span></div>
118
119
<span style="text-align: center; font-size: 75%;">Se ha observado que no es aconsejable que los intervalos de tiempo </span> <math display="inline">{\Delta }_{I}t</math><span style="text-align: center; font-size: 75%;"> y </span> <math display="inline">{\Delta }_{E}t</math><span style="text-align: center; font-size: 75%;"> tengan una relación mayor que 1: 500 ya que pueden haber discrepancias entre la solución explícita y la implícita. Teniendo esto en cuenta, la estrategia de solución en sub-pasos implementada permite obtener excelentes resultados, como se muestra en los ejemplos del siguiente apartado.</span>
120
121
=7.- Ejemplos=
122
123
En este apartado se muestran varios ejemplos a fin de mostrar el buen comportamiento de la estrategia FEM-DEM descrita. El primer ejemplo corresponde al estudio 3D de una probeta normalizada en un ensayo a tracción. El segundo ejemplo es una viga de hormigón bi-entallada donde predomina la fractura en modo mixto. El tercer ejemplo corresponde al ensayo de tracción indirecta, ampliamente usado en mecánica de rocas. El cuarto ejemplo consiste en un ensayo de cortante en hormigón propuesto por Luong [21]. El quinto ejemplo es un ensayo de una probeta a compresión simple. El sexto ejemplo es el estudio de un forjado reticular de hormigón armado, afectado por un asentamiento diferencial. El acero se considera embebido y solidario con los desplazamientos de los elementos finitos utilizados.
124
125
<span style="text-align: center; font-size: 75%;">Finalmente, se presenta el análisis sísmico de la nave central de la iglesia del monasterio de Poblet , formada por elementos de mampostería, aplicándole un sismo de 6 </span> <math display="inline">{M}_{w}</math><span style="text-align: center; font-size: 75%;">. En este ejemplo se puede observar como la estrategia FEM-DEM modela el efecto de la apertura y cierre de múltiples grietas.</span>
126
127
==7.1 Ensayo normalizado de tracción==
128
129
El primer ejemplo corresponde al análisis de fractura de una probeta de hormigón sujeta a fuerzas de tracción. El objetivo principal es mostrar la independencia del tamaño de la malla en la generación de la grieta, de manera que la energía utilizada en la fractura sea independiente del tamaño elemental. La geometría se define de acuerdo con la norma D638 de la Sección Norteamericana de la Asociación Internacional de Ensayo de Materiales (ASTM, American Society for Testing and Materials) [22]. En la Figura 3 se muestran las tres mallas de elementos tetraédricos de 4 nodos utilizadas, así como las condiciones de contorno. La probeta se carga imponiendo un campo de velocidad constante de tracción en sus extremos, representados por la zona oscura en la figura.
130
131
El estudio se ha realizado utilizando la metodología FEM-DEM en 3D antes descrita. Con el fin de localizar la fractura, sólo se permite que una banda de elementos se rompa al nivel de la tensión de fallo, usando el modelo de daño mencionado. Los resultados obtenidos se analizan dibujando los desplazamientos de los puntos PA y PB mostrados en la Figura 3.
132
133
<span style="text-align: center; font-size: 75%;">El módulo de Young, el coeficiente de Poisson y la densidad son respectivamente </span> <math display="inline">{E}_{0}</math><span style="text-align: center; font-size: 75%;"> = 30 &#x00d7; 10<sup>9</sup>Pa, </span> <math display="inline">\nu</math> <span style="text-align: center; font-size: 75%;"> = 0.2, </span> <math display="inline">\gamma</math> <span style="text-align: center; font-size: 75%;">= 1.0 &#x00d7; 10<sup>3</sup> N/m<sup>3</sup>, la tensión máxima de tracción </span> <math display="inline">{\sigma }_{f}</math><span style="text-align: center; font-size: 75%;"> = 10 &#x00d7; 10<sup>3</sup>Pa y la energía de fractura ''G'' = 7.5 &#x00d7; 10<sup>-</sup><sup>3</sup> J/m<sup>2</sup>. </span>
134
135
La probeta se deforma aplicando una velocidad constante de tracción de 0.5&#x00d7;10<sup>-</sup><sup>7</sup> m/s en ambos extremos . La Figura 4 muestra la geometría dañada para las tres mallas del FEM consideradas.
136
137
Obsérvese que cuando se produce una fractura, se crean elementos discretos en los labios de la única fisura que aparece en este ejemplo, como se explica en los apartados anteriores. El tamaño de las esferas no es uniforme debido a que no todos los elementos llegan al daño máximo de manera simultánea.
138
139
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
140
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image3.jpeg|396px]] </span></div>
141
142
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
143
<span style="text-align: center; font-size: 75%;">Figura 3. Ensayo normalizado 3D de tracción. Malla de elementos finitos y dimensiones según la norma ASTMD638.</span></div>
144
145
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
146
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image4.jpeg|318px]] </span></div>
147
148
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
149
<span style="text-align: center; font-size: 75%;">Figura 4. Ensayo normalizado de tracción. Zona agrietada incluyendo los elementos discretos generados para las tres mallas de elementos finitos consideradas.</span></div>
150
151
[[Image:Draft_Samper_634227043-image5.jpeg|600px]]
152
153
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
154
<span style="text-align: center; font-size: 75%;">Figura 5. Ensayo normalizado 3D de tracción. Curva carga-desplazamiento de los puntos PA y PB de la muestra.</span></div>
155
156
Con objeto de evaluar la apertura de la grieta se analiza el desplazamiento de los puntos PA y PB situados a la derecha y en el centro de la probeta, respectivamente (Figura 3). La Figura 5 muestra la relación carga-desplazamiento en estos puntos. Para las tres mallas consideradas la evolución del desplazamiento es muy similar y de acuerdo con los resultados esperados [9]. Debido a que los elementos por donde se abre la grieta tienen un tamaño diferente para cada malla, el desplazamiento del punto PB en la región elástica se hace más pequeño a medida que se reduce el tamaño del elemento.
157
158
==7.2 Viga bi-entallada a flexión.==
159
160
El siguiente ejemplo corresponde al ensayo de una viga de hormigón en masa bi-entallada. El análisis se realiza mediante las hipótesis de tensión plana y es un buen ejemplo de fractura en modo mixto. La viga se sostiene en dos puntos y se somete a flexión aplicando un desplazamiento impuesto mediante control de velocidad equivalente a 1 mm/s en los dos puntos representados en la Figura 6. Igualmente, en dicha figura también se muestran la geometría y las dimensiones de la probeta.
161
162
<span style="text-align: center; font-size: 75%;">La viga presenta dos puntos singulares en la punta de las dos entallas en donde las tensiones de tracción son altas y el daño comienza en esta zona. Las propiedades del material son </span> <math display="inline">{E}_{0}</math><span style="text-align: center; font-size: 75%;"> = 30 &#x00d7; 10<sup>9</sup>Pa, </span> <math display="inline">\nu</math> <span style="text-align: center; font-size: 75%;"> = 0.2, </span> <math display="inline">{\sigma }_{f}</math><span style="text-align: center; font-size: 75%;"> = 2 MPa y ''G''= 1 &#x00d7; 10<sup>2</sup> J/m<sup>2</sup>. El problema ha sido resuelto con la técnica FEM-DEM en 2D.</span>
163
164
En la Figura 7 se muestra un detalle de las tres diferentes mallas utilizadas, formadas por 1165 nodos y 2202 elementos triangulares lineales para la malla gruesa, 1847 nodos y 3480 elementos para la malla intermedia y 5747 nodos y 11206 elementos para la malla fina. El análisis se ha realizado tanto forma cuasiestática como de forma dinámica, respetando la velocidad de aplicación de la carga. En ambos casos los resultados han sido muy similares como se observa en la Figura 9.
165
166
La Figura 8 muestra la dirección de las fisuras para las tres mallas analizadas las cuales coinciden con los experimentos numéricos [14]. La Figura 9 muestra la relación entre la reacción y el desplazamiento impuesto en cualquiera de los dos puntos representados en la figura 6 (los resultados son idénticos para los dos puntos). Los gráficos son concordantes con los resultados obtenidos por Cervera et. al. [14].
167
168
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
169
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image6.jpeg|522px]] </span></div>
170
171
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
172
<span style="text-align: center; font-size: 75%;">Figura 6. Viga bi-entallada a flexión. Geometría y condiciones de contorno.</span></div>
173
174
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
175
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image7.jpeg|516px]] </span></div>
176
177
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
178
<span style="text-align: center; font-size: 75%;">Figura 7. Viga bi-entallada a flexión. Mallas utilizadas.</span></div>
179
180
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
181
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image8.jpeg|504px]] </span></div>
182
183
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
184
<span style="text-align: center; font-size: 75%;">Figura 8. Viga bi-entallada a flexión. Superficies de igual desplazamiento y apertura de las grietas. </span></div>
185
186
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
187
<span style="text-align: center; font-size: 75%;">a) Malla gruesa, b) malla intermedia, c) Malla fina.</span></div>
188
189
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
190
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image9.jpeg|600px]] </span></div>
191
192
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
193
<span style="text-align: center; font-size: 75%;">Figura 9. Viga bi-entallada a flexión. Relación entre la fuerza y el desplazamiento en cualquiera de los puntos referenciados en la Figura 6. Los resultados obtenidos son comparados con los aquellos dados en [14].</span></div>
194
195
==7.3 Ensayo de tracción indirecta==
196
197
El ensayo brasileño de tracción indirecta (BTS, de Brazilian Tensile Strength) es un procedimiento sencillo para evaluar la resistencia a la tracción de hormigón y geomateriales. La probeta de hormigón analizada es un cilindro de 0.2m de diámetro (''D'') y 0.1m de espesor (''t''), sujeto a una carga diametralmente opuesta (Figura 10). El valor de la resistencia a la tracción se calcula mediante la siguiente expresión [23, 24].
198
199
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
200
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image10.jpeg|372px]] </span></div>
201
202
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
203
<span style="text-align: center; font-size: 75%;">Figura 10. Ensayo de tracción indirecta. Dimensiones de la muestra, condiciones de contorno y mallas usadas de elementos finitos tetraédricos de cuatro nodos</span></div>
204
{| class="formulaSCP" style="width: 100%; text-align: center;" 
205
|-
206
| <math>{\sigma }_{f}=\frac{2P}{\pi tD}\quad \, \, (4)</math>
207
|}
208
209
210
<span style="text-align: center; font-size: 75%;">donde ''P'' es el valor de la carga aplicada. Las propiedades del material son </span> <math display="inline">{E}_{0}</math><span style="text-align: center; font-size: 75%;"> = 21 &#x00d7; 10<sup>9</sup>Pa, </span> <math display="inline">\nu</math> <span style="text-align: center; font-size: 75%;">= 0.2,  </span> <math display="inline">\gamma</math> <span style="text-align: center; font-size: 75%;">= 7.8 &#x00d7; 10<sup>3</sup>N/m<sup>3</sup>, </span> <math display="inline">{\sigma }_{f}</math><span style="text-align: center; font-size: 75%;"> = 10 KPa y ''G''= 1 &#x00d7; 10<sup>2</sup> J/m<sup>2</sup>. Lo que proporciona una carga máxima de fallo de ''P'' = 314.16 N.</span>
211
212
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
213
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image11.jpeg|384px]] </span></div>
214
215
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
216
<span style="text-align: center; font-size: 75%;">Figura 11. Ensayo de tracción indirecta. Grieta y elementos discretos generados. a) Malla gruesa, b) malla intermedia, c) malla fina.</span></div>
217
218
Para realizar los análisis se han usado tres mallas de 9338, 31455 y 61623 elementos tetraédricos lineales como se muestra en la Figura 10. El ensayo se realiza imponiendo una velocidad constante vertical en la parte superior de la probeta.
219
220
<span style="text-align: center; font-size: 75%;">La Figura 11 muestra la grieta y los elementos discretos generados. Se puede observar que el patrón de fisuración es similar para las tres mallas y de acuerdo con el resultado esperado. La Figura 12 muestra la curva de carga-desplazamiento. Los valores obtenidos para la resistencia a la tracción máxima para las mallas gruesa, media y fina son respectivamente: 10.693 KPa, 10.351 KPa y 10.235 KPa, correspondientes a un rango de entre 6% y 2% de error frente al valor esperado de </span> <math display="inline">{\sigma }_{f}</math><span style="text-align: center; font-size: 75%;"> = 10 KPa.</span>
221
222
Es remarcable la insensibilidad de la curva carga-desplazamiento al tipo de malla de elementos finitos utilizada. Esta “objetividad” de los resultados numéricos frente al tamaño de la malla del FEM es una de las características esenciales de la técnica FEM-DEM.
223
224
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
225
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image12.jpeg|402px]] </span></div>
226
227
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
228
<span style="text-align: center; font-size: 75%;">Figura 12. Ensayo de tracción indirecta. Relación fuerza-desplazamiento para las tres mallas utilizadas.</span></div>
229
230
==7.4 Ensayo de cortante==
231
232
<span id='_GoBack'>El ensayo de cortante se diseña para aplicar un esfuerzo de cizalladura sobre una probeta de modo que experimente una falla por deslizamiento a lo largo de un plano paralelo a las fuerzas aplicadas. Por lo general las fuerzas de cizallamiento provocan que una de las superficies de fallo de un material se mueva en una dirección y la otra superficie en dirección opuesta, de manera que el material se encuentra sometido a un estado de corte. Este ejemplo se realizó en 80 minutos usando un procesador a 2.5 MHz</span>
233
234
El ensayo que aquí se considera tiene como objetivo determinar la resistencia a cortante del hormigón y ha sido propuesto por Luong [21]. La probeta tiene forma tubular y su eje coincide con el eje z = 0. Tiene varias entallas y se somete a una carga central en una de sus caras y otra excéntrica en la cara opuesta, de manera que se generen tensiones cortantes paralelas al eje z = 0, como se describe en la Figura 13 en la que también se muestran las condiciones de contorno impuestas. La profundidad de la entalla es de 10mm y su ancho es de 4mm. La carga se aplica imponiendo a la placa superior una velocidad constante de 1mm/s hasta llegar a la fractura.
235
236
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
237
<span style="text-align: center; font-size: 75%;">.</span></div>
238
239
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
240
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image13.jpeg|600px]] </span></div>
241
242
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
243
<span style="text-align: center; font-size: 75%;">Figura 13. Ensayo de cortante. Geometría, cargas y condiciones de contorno.</span></div>
244
245
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
246
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image14.jpeg|450px]] </span></div>
247
248
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
249
<span style="text-align: center; font-size: 75%;">Figura 14. Ensayo de cortante. a) Definición geométrica de la probeta b) Malla de elementos finitos.</span></div>
250
251
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
252
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image15.jpeg|474px]] </span></div>
253
254
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
255
<span style="text-align: center; font-size: 75%;">Figura 15. Ensayo de cortante. a) Vista superior b) Perspectiva de la grieta obtenida por la técnica FEM-DEM. </span></div>
256
257
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
258
<span style="text-align: center; font-size: 75%;">c) Rotura en probetas de laboratorio.</span></div>
259
260
La definición geométrica y la malla de elementos finitos utilizada se muestran en la Figura 14. El plano de cizallamiento ha sido discretizado con 4 elementos finitos tetraédricos de 4 nodos con el fin de captar adecuadamente el gradiente de tensiones en esta zona.
261
262
<span style="text-align: center; font-size: 75%;">Las propiedades del hormigón son </span> <math display="inline">{E}_{0}</math><span style="text-align: center; font-size: 75%;"> = 35 &#x00d7; 10<sup>9</sup>Pa, </span> <math display="inline">\nu</math> <span style="text-align: center; font-size: 75%;"> = 0.22, </span> <math display="inline">\gamma</math> <span style="text-align: center; font-size: 75%;"> = 7.8 &#x00d7; 10<sup>3</sup>  N/m<sup>3</sup>, </span> <math display="inline">{\sigma }_{c}</math><span style="text-align: center; font-size: 75%;">  = 30 MPa, </span> <math display="inline">{\sigma }_{f}</math><span style="text-align: center; font-size: 75%;"> = 3 MPa y G = 75 &#x00d7; 10<sup>-3</sup> J/m<sup>2</sup>. En general el esfuerzo máximo de corte en el hormigón es aproximadamente 1/5 o 1/6 del esfuerzo a compresión [21].</span>
263
264
La Figura 15 muestra la grieta obtenida por el experimento numérico con la técnica FEM-DEM y se compara con los resultados de los ensayos realizados en el laboratorio. La Figura 16 muestra la relación fuerza-desplazamiento obtenida numéricamente y en la que se muestra claramente la rama elástica hasta que los elementos comienzan a dañar. Considerando el tamaño de la probeta, de esta gráfica se desprende que el esfuerzo cortante último alcanza un valor de 5.10 MPa. Este resultado es coherente y próximo al valor esperado [21].
265
266
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
267
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image16.jpeg|462px]] </span></div>
268
269
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
270
<span style="text-align: center; font-size: 75%;">Figura 16. Ensayo de cortante. Relación fuerza-desplazamiento obtenida por el ensayo numérico.</span></div>
271
272
==7.5 Ensayo de compresión simple==
273
274
Uno de los ensayos más frecuentes realizados en probetas de hormigón es el de compresión simple [25]. El ensayo se realiza de acuerdo con la instrucción de hormigón estructural (EHE-08) [26] la cual hace referencia a la norma UNE-EN 12390-3.2009 [27] en donde se especifica las dimensiones de las probetas y las condiciones de ensayo.
275
276
El experimento numérico que se presenta se realiza en 2D utilizando probetas hexaédricas en tensión plana y para el caso 3D probetas cilíndricas normalizadas.
277
278
El ensayo 2D consiste en comprimir la probeta cuyas dimensiones y condiciones de carga se muestra en la Figura 17 y en la cual también se presentan las mallas utilizadas. El objetivo del ensayo no solo se trata de encontrar la tensión máxima a compresión que es capaz de resistir la probeta sino también la forma de las fisuras en función de las distintas condiciones de contorno que se pueden dar en el ensayo.
279
280
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
281
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image17.jpeg|474px]] </span></div>
282
283
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
284
<span style="text-align: center; font-size: 75%;">Figura 17. Ensayo de compresión simple. Probeta y mallas utilizadas en el ensayo 2D.</span></div>
285
286
<span style="text-align: center; font-size: 75%;">Las propiedades del hormigón utilizadas son. </span> <math display="inline">{E}_{0}</math><span style="text-align: center; font-size: 75%;"> = 30.0 &#x00d7; 10<sup>9</sup>Pa, </span> <math display="inline">\nu</math> <span style="text-align: center; font-size: 75%;"> = 0.20, </span> <math display="inline">\gamma</math> <span style="text-align: center; font-size: 75%;"> = 24 &#x00d7; 10<sup>3</sup>N/m<sup>3</sup>,</span> <math display="inline">\, {\sigma }_{c}</math><span style="text-align: center; font-size: 75%;"> =  2000 MPa, </span> <math display="inline">{\sigma }_{t}</math><span style="text-align: center; font-size: 75%;">= 20.0 MPa y G = 105 &#x00d7; 10<sup>-3</sup> J/m<sup>2</sup> En este caso la tensión máxima de compresión, independientemente de las condiciones de apoyo es de 20 &#x00d7; 10<sup>6</sup>Pa. En la Figura 18 se muestra que los resultados obtenidos con ambas probetas se acercan mucho al resultado esperado.</span>
287
288
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
289
<span style="text-align: center; font-size: 75%;">''' [[Image:Draft_Samper_634227043-image18.jpeg|390px]] '''</span></div>
290
291
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
292
<span style="text-align: center; font-size: 75%;">Figura 18. Ensayo de compresión simple. Relación tensión-deformación para el análisis 2D.</span></div>
293
294
Los resultados más interesantes en este caso se muestran en las Figuras 19 a 21. En todas ellas la figura a) corresponde al patrón de grietas generado, la b) son las superficies de igual desplazamiento y la figura c) el daño elemental. Estos resultados se pueden comparar con la figura d) que corresponde a los resultados teóricos [28] o bien con la figura e) que corresponde a los resultados obtenidos con LS-DYNA [29] para el caso bi-empotrado (Figuras 19 y 20) o con un ensayo real (Figura 21)
295
296
La Figura 19 corresponde al caso en que las superficies de carga tienen un desplazamiento horizontal nulo utilizando la malla de elementos finitos estructurada. Los resultados muestras claramente dos grietas a 45º claramente definidas, y corroboradas por el patrón de daño mostrado. El resultado numérico mostrado en la figura e) muestra el mismo patrón. Sin embargo, observando los resultados teóricos [28] el daño se genera sobre dos bandas a 45º formadas por grietas verticales.
297
298
La Figura 20 muestra el mismo caso, con el desplazamiento horizontal de las superficies de carga impedidas, pero utilizando la malla no estructurada. Como se puede observar en la Figura 20 a) el patrón de grietas se asemeja más a los resultados teóricos esperados, con grietas verticales localizadas en dos bandas a 45º. Es importante mencionar que ningún código numérico al que los autores han tenido acceso es capaz de presentar este tipo de resultados, pues todos muestran correctamente la banda de daño (Figura e) pero no así las fisuras verticales sobre dicha banda. El por qué la malla no estructurada presenta un resultado más próximo al teórico está relacionado con la no uniformidad en la discretización de dominio. Las pequeñas variaciones numéricas inducidas por los elementos de distinto tamaño son suficientes para alterar la estabilidad de los resultados que se obtienen con una malla estructurada, encontrado así un resultado diferente.
299
300
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
301
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image19.jpeg|504px]] </span></div>
302
303
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
304
<span style="text-align: center; font-size: 75%;">Figura 19. Ensayo de compresión simple. Probeta biempotrada, malla estructurada.</span></div>
305
306
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
307
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image20.jpeg|516px]] </span></div>
308
309
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
310
<span style="text-align: center; font-size: 75%;">Figura 20. Ensayo de compresión simple. Probeta biempotrada, malla no estructurada.</span></div>
311
312
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
313
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image21.jpeg|504px]] </span></div>
314
315
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
316
<span style="text-align: center; font-size: 75%;">Figura 21. Ensayo de compresión simple. Probeta simplemente apoyada, malla no estructurada.</span></div>
317
318
La Figura 21 presenta los resultados obtenidos con la malla no estructurada en el caso en qué las superficies de carga tengan libertad de movimiento horizontal. Teóricamente [28], como se muestra en la Figura 21 d), la fisura debe ser vertical, sin embargo, como se observa en la Figura 21 e), los ensayos de laboratorio tienden a presentar varias fisuras verticales [27]. Los resultados mostrados en la Figura 21 a) concuerdan con los resultados esperados.
319
320
El mismo ejemplo ha sido resuelto en 3D utilizando una probeta circular (Figura 17) y las mismas propiedades mecánicas. La malla está formada por 926 nodos y 4030 tetraedros. Los resultados obtenidos se muestran en la Figura 22 a) para el caso en que una de las superficies de carga se encuentre impedida en su desplazamiento horizontal, mientras que la otra no. La Figura 22 b) muestra el caso en que ambas superficies tengan el desplazamiento horizontal impedido. En ambos ejemplos la tensión máxima se aproxima bastante a la esperada. No solo eso, a pesar que la malla sea tan gruesa, las grietas generadas coinciden con los resultados teóricos. Cabe comentar que este ejemplo no requiere más de 5 minutos de cálculo en un procesador a 2.5 GHz.
321
322
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
323
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image22.jpeg|444px]] </span></div>
324
325
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
326
<span style="text-align: center; font-size: 75%;">Figura 22. Ensayo de compresión simple. Probeta 3D para distintos casos de apoyo. a) Bases empotrada-libre, b) Bases bi-empotradas.</span></div>
327
328
==7.6 Losa reticular==
329
330
El uso de forjados reticulares como elementos estructurales está ampliamente consensuado y permite construir geometrías arquitectónicas más libres. El forjado analizado corresponde al techo de una vivienda. La planta es un cuadrado de 22m de lado con una superficie plana ovalada en el centro y los lados en parteaguas con un desnivel de 2m. En la Figura 23 se muestra la planta estructural en donde se puede observar claramente la armadura de acero. Básicamente se trata de un forjado reticular de 0.25m de espesor con una capa de compresión de 0.05m y casetones de 0.40 &#x00d7; 0.40 &#x00d7; 0.20m. A lo largo de las nervaduras se colocan varillas de acero del No. 3 (ϕ = 0.0095m).
331
332
Adicionalmente existe una serie de cerramientos (D1) y vigas (T1 y T2) así como cuatro capiteles de columna. Los cerramientos tienen una sección de 0.25 &#x00d7; 0.20m armados con 4 varillas de acero del No. 3. Las vigas T1 y T2 tienen un espesor de 0.25m y un ancho de 0.20m para la viga T1 y de 0.25m para la viga T2. Ambas están armadas con 4 varillas de acero del No. 4 (ϕ = 0.0127m). Finalmente, el armado de los capiteles se realiza con varillas del No. 3 [mailto:@ @] 0.15m en ambos lechos y en ambas direcciones. La Figura 24 a) muestra el armado de toda la losa, así como un detalle del armado de los capiteles en la Figura 24 b).
333
334
<span style="text-align: center; font-size: 75%;">Las propiedades del hormigón son. </span> <math display="inline">{E}_{0}</math><span style="text-align: center; font-size: 75%;"> = 21 &#x00d7; 10<sup>9</sup>Pa, </span> <math display="inline">\nu</math> <span style="text-align: center; font-size: 75%;"> = 0.20, </span> <math display="inline">\gamma</math> <span style="text-align: center; font-size: 75%;"> = 24 &#x00d7; 10<sup>3</sup> N/m<sup>3</sup>, </span> <math display="inline">{\sigma }_{c}</math><span style="text-align: center; font-size: 75%;">= 20 MPa, </span> <math display="inline">{\sigma }_{f}</math><span style="text-align: center; font-size: 75%;">= 2 MPa y G = 100 &#x00d7; 10<sup>-3</sup> J/m<sup>2</sup>.</span>
335
336
La placa se encuentra apoyada sobre los muros mostrados en la Figura 25. A efectos del análisis dichos muros se consideran lo suficientemente rígidos y empotrados en su base.
337
338
Las cargas consideradas son el peso propio de la estructura, las cargas muertas 1620 N/m<sup>2</sup> y las cargas vivas 400 N/m<sup>2</sup>. Dichos valores consideran el peso del plafón, las instalaciones y acabados finales, sin considerar cargas de nieve o viento. Adicionalmente se impone un asentamiento diferencial de 0.0035 m de los puntos de la base del muro como se indica en la Figura 25.
339
340
[[Image:Draft_Samper_634227043-image23.jpeg|600px]]
341
342
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
343
<span style="text-align: center; font-size: 75%;">Figura 23. Losa reticular. Plano estructural.</span></div>
344
345
[[Image:Draft_Samper_634227043-image24.jpeg|600px]]
346
347
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
348
<span style="text-align: center; font-size: 75%;">Figura 24. Losa reticular. a) Parrilla del armado de la losa. b) Detalle del armado en los capiteles.</span></div>
349
350
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
351
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image25.jpeg|288px]] </span></div>
352
353
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
354
<span style="text-align: center; font-size: 75%;">Figura 25. Losa reticular. Localización de los muros de carga. El muro con un asentamiento diferencial se encuentra marcado como DIFF.</span></div>
355
356
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
357
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image26.jpeg|438px]] </span></div>
358
359
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
360
<span style="text-align: center; font-size: 75%;">Figura 26. Losa reticular. a) Vista del conjunto. b) Detalle de las nervaduras cercanas al lucernario.</span></div>
361
362
El problema se ha resuelto utilizando una malla de 1,369,192 tetraedros y 368,969 nodos. Por otra parte, se han definido 5,084 elementos lineales para representar las varillas de acero. La Figura 26 a) representa el conjunto del forjado y muros a analizar mientras que la Figura 26 b) corresponde a un detalle cercano al lucernario donde se aprecian las nervaduras de la placa, así como algunas de las vigas T2.
363
364
la Figura 27 a) muestra los resultados del desplazamiento de la estructura después de que el asentamiento diferencial de la base sea de 0.0035 m. En la Figura 27 b) se aprecia un corte muy definido en la cubierta superior. Se observa que el desgarro en la cubierta es perpendicular a la dirección del muro, poniendo de manifiesto la gran rigidez que tiene la cubierta superior.
365
366
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
367
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image27.jpeg|600px]] </span></div>
368
369
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
370
<span style="text-align: center; font-size: 75%;">Figura 27. Losa reticular. a) Vista superior de la estructura con iso-superficies de igual desplazamiento. b) Vista superior de las grietas generadas.</span></div>
371
372
En la Figura 28 a) se puede apreciar la deformada de la estructura vista desde abajo. Es interesante observar que el lucernario sufre los desplazamientos máximos, incluso mayores a los que sufre el muro sujeto al asentamiento diferencial. La Figura 28 b) muestra una vista inferior de la losa donde se aprecia el daño producido en la unión del muro con ésta.
373
374
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
375
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image28.jpeg|600px]] </span></div>
376
377
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
378
<span style="text-align: center; font-size: 75%;">Figura 28. Losa reticular. a) Vista inferior de la estructura deformada b) Vista inferior de las grietas generadas.</span></div>
379
380
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
381
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image29.jpeg|336px]] </span></div>
382
383
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
384
<span style="text-align: center; font-size: 75%;">Figura 29. Losa reticular. Elementos del refuerzo de acero en la zona cercana al lucernario.</span></div>
385
386
En la Figura 29 se observan los esfuerzos axiales sobre los elementos de acero pertenecientes a la viga T2 en la zona cercana al lucernario. En esta figura se puede observar los elementos en tracción y compresión debido a gran flexión a la que está sujeta esta parte de la estructura.
387
388
La Figura 30 muestra las zonas dañadas de la estructura utilizando una visualización de la estructura deformada con un factor de escala de 5000. El recuadro a) muestra una vista inferior de la losa en donde se puede observar la separación entre el muro y la losa, mientras que el recuadro b) muestra una vista superior mostrando el desgarro de la cubierta. En vista de los resultados obtenidos es claramente notorio la gran rigidez que posee la losa
389
390
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
391
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image30.jpeg|600px]] </span></div>
392
393
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
394
<span style="text-align: center; font-size: 75%;">Figura 30. Losa reticular. Daño en la estructura, a) vista inferior. b) vista superior</span></div>
395
396
==7.7 Iglesia cisterciense del monasterio de Poblet==
397
398
La iglesia, perteneciente al Real Monasterio de Santa Maria de Poblet, ubicado en Vimbod (Tarragona), es un monumento histórico artístico declarado por la UNESCO como patrimonio de la humanidad. Construido a partir de la segunda mitad del siglo XII, el templo adopta una planta basilical orientado su ábside al Este. La integran tres naves de siete tramos, con crucero, ábside central, girola y capillas absidiales. La nave central tiene unas dimensiones de 85m de longitud, 21m de anchura y 28m de altura, mientras que las laterales alcanzan los 18m de altura. La diferencia de altura respecto a la nave central no se resuelve mediante arbotantes del estilo gótico, sino con contrafuertes. La nave central plenamente románica, está cubierta con bóveda de cañón apuntada, con arcos fajones en cada tramo como se muestra en la Figura 31.
399
400
<span style="text-align: center; font-size: 75%;">Para simplificar el estudio se analiza únicamente una sección de la nave central que se considera simétrica respecto al eje normal a la dirección de dicha nave. La malla que discretiza la iglesia tiene 9294 nodos y 32487 elementos tetraédricos de cuatro nodos como se muestra en la Figura 32 a). Las propiedades del material utilizado se corresponden con una piedra caliza propia de la región, con  </span> <math display="inline">{E}_{0}</math><span style="text-align: center; font-size: 75%;"> = 35 &#x00d7; 10<sup>9</sup>Pa, </span> <math display="inline">\nu</math> <span style="text-align: center; font-size: 75%;">= 0.22, </span> <math display="inline">\gamma</math> <span style="text-align: center; font-size: 75%;"> = 25 &#x00d7; 103 N/m<sup>3</sup>, </span> <math display="inline">{\sigma }_{c}</math><span style="text-align: center; font-size: 75%;"> = 780 MPa, </span> <math display="inline">{\sigma }_{f}</math><span style="text-align: center; font-size: 75%;"> = 5 MPa y G = 100 &#x00d7; 10<sup>-3</sup> J/m<sup>2</sup>. En cuanto a las condiciones de apoyo, se considera empotrada toda la superficie inferior del modelo.</span>
401
402
Las cargas aplicadas corresponden a un movimiento oscilatorio impuesto en dirección normal al eje de la nave central, como se muestra en la Figura 32 b). Los valores mostrados corresponden al movimiento oscilatorio en dirección S 00E de los registros en desplazamientos obtenidos del acelerograma del sismo “El Centro” del 18 de mayo de 1940 [30]. Estos valores se aplican al modelo con una reducción del 50%, de forma que simule un sismo de magnitud6.0 M<sub>W</sub>.
403
404
Aunque el registro de movimientos dura más de 50 segundos la estructura solo es capaz de soportar en pie 4 segundos, colapsando de forma total a partir de ese instante. En la Figura 33 se observa a intervalos de un segundo las grietas formadas en la Iglesia. Inicialmente ocurre un cizallamiento de las columnas y de la pared más rígida, para posteriormente, cortar las bases de los arcos. A partir del segundo cinco se observa como la estructura colapsa en su totalidad reduciéndose a escombros.
405
406
Este ejemplo muestra claramente la capacidad de la técnica FEM-DEM para simular la aparición y propagación de múltiples fracturas en estructuras de mampostería, así como el eventual colapso de la estructura bajo cargas dinámicas.
407
408
=8.- Conclusiones=
409
410
Se ha descrito en el presente artículo las líneas generales de la metodología FEM-DEM propuesta por los autores en [9, 10] para la predicción de la aparición y evolución de fisuras en estructuras de hormigón. Los ejemplos que se han presentado muestran las posibilidades de la técnica FEM-DEM para el cálculo no lineal de estructuras de hormigón en masa y armado, así como en estructuras de mampostería. La técnica FEM-DEM es también aplicable al estudio de la fractura en macizos rocosos [31, 32].
411
412
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
413
<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Samper_634227043-image31.jpeg|534px]] </span></div>
414
415
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
416
<span style="text-align: center; font-size: 75%;">Figura 31. Iglesia cisterciense del monasterio de Poblet. a) Vista interior de la nave central. b) Planta y sección analizada</span></div>
417
418
[[Image:Draft_Samper_634227043-image32.jpeg|600px]]
419
420
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
421
<span style="text-align: center; font-size: 75%;">Figura 32. Iglesia cisterciense del monasterio de Poblet. a) Malla de elementos finitos de una sección perteneciente a la nave principal. b) Excitación temporal en desplazamientos</span></div>
422
423
[[Image:Draft_Samper_634227043-image33.jpeg|600px]]
424
425
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
426
<span style="text-align: center; font-size: 75%;">Figura 33. Iglesia cisterciense del monasterio de Poblet sometida a un terremoto. Evolución de las fracturas en la nave central en distintos tiempos</span></div>
427
428
=9.- Agradecimientos=
429
430
Los autores agradecen la colaboración del Ing. Juan José Cuellar Ornelas al facilitar de forma desinteresada la información usada en el análisis del forjado. Los resultados aquí presentados han sido obtenidos utilizando los programas FEM2DEM y DEMPACK ([http http].[http:////www.cimne.com/dempack //www.cimne.com/dempack]) en los que se ha implementado la metodología FEM-DEM descrita.
431
432
=REFERENCES=
433
434
<div id="cite-1"></div>
435
[[#citeF-1|[1]]] P. A. Cundall and O. D. L. Strack, A discrete numerical model for granular assemblies, ''Geotechnique '' '''29''' (1979), no. 1, 47 –65.
436
437
<div id="cite-2"></div>
438
[[#citeF-2|[2]]] C. Labra, ''Advances in the development of the discrete element method for excavation processes'', Ph.D. Thesis, Barcelona, 2012.
439
440
<div id="cite-3"></div>
441
[[#citeF-3|[3]]] J. Rojek, E.. Oñate, F. Zarate, and J. Miquel, Modelling of rock, soil and granular materials using spherical elements, ''2nd European Conference on Computational Mechanics ECCM- 2001'', Cracow (2629 June 2001).
442
443
<div id="cite-4"></div>
444
[[#citeF-4|[4]]] J. Williams and R. O’Connor, Discrete element simulation and contact problem, ''Archives of Computer Methods in Engineering'' '''6''' (1999), no. 4, 279 –304.
445
446
<div id="cite-5"></div>
447
[[#citeF-5|[5]]] E. Oñate, F. Zárate, J. Miquel, M. Santasusana, M.A. Celigueta, F. Arrufat, R. Gandijota, K. Valiullin, and L. Ring, A local constitutive model for the discrete element method. Application to geomaterials and concrete, ''Computational Particle Mechanics'' '''2''' (2015), 139–160.
448
449
<div id="cite-6"></div>
450
[[#citeF-6|[6]]] S. Katagiri and S. Takada, Development of fem-dem combined method for fracture analysis of a continuos media, ''Memoirs of the Graduate School of Science and Technology, Kobe University Japan''. '''20A''' (2002 –03), 65 –79.
451
452
<div id="cite-7"></div>
453
[[#citeF-7|[7]]] A. Munjiza, ''The combined finite-discrete element method'', isbn 0-470-84199-0, Wiley, 2004.
454
455
<div id="cite-8"></div>
456
[[#citeF-8|[8]]] Shmauder S. Wulf J. and Fischmeister H.F., Finite element modelling of crack propagation in ductile fracture, ''Computation Materials Science'' '''1''' (1993), 297 –301.
457
458
<div id="cite-9"></div>
459
[[#citeF-9|[9]]] F. Zárate and E. Oñate, A simple FEM-DEM technique for fracture prediction in materials and structures, ''Computational particle mechanics'' '''2''' (2015), no. 3, 301–314.
460
461
<div id="cite-10"></div>
462
[[#citeF-10|[10]]] F. Zárate, A. Cornejo, and E. Oñate, A three-dimensional FEM-DEM technique for predicting the evolution of fracture in geomaterials and concrete, ''Computational particle mechanics ''(2017), DOI 10.1007/s40571-017-0178-z
463
464
<div id="cite-11"></div>
465
[[#citeF-11|[11]]] E. Oñate. Cálculo de estructuras por el método de elementos finitos. Análisis estático lineal. Vol. 1 Sólidos (CIMNE, 2017, en castellano); Vol 2. Placas y láminas (Springer 2009, en inglés)
466
467
<div id="cite-12"></div>
468
[[#citeF-12|[12]]] M. Cervera, M. Chiumenti, and C. Agelet de Saracibar, Shear band localization via local j2 continuum damage mechanics, ''Comput. Methods Appl. Mech. Engrg''. '''193''' (2004), 849 –880.
469
470
<div id="cite-13"></div>
471
[[#citeF-13|[13]]] M. Cervera, M. Chiumenti, and R. Codina, Mixed stabilized finite element methods in nonlinear solid mechanic’s part I. Formulation., ''Computer Methods in Applied Mechanics and Engineering'' '''199''' (2010), 2559 –2570.
472
473
<div id="cite-14"></div>
474
[[#citeF-14|[14]]] M. Cervera, M. Chiumenti, and R. Codina, Mesh objective modelling of cracks using continuous linear strain and displacements interpolations, ''Int. J. Numer. Meth. Engng''. '''87''' (2011), 962 –987.
475
476
<div id="cite-15"></div>
477
[[#citeF-15|[15]]] P.R. Johnson, N. Petrinic, and E. Sli, Element –splitting for simulation of fracture in 3d solid continua, ''VIII International Conference on Computational Plasticity'', Barcelona. (2005).
478
479
<div id="cite-16"></div>
480
[[#citeF-16|[16]]] L. Mishnaevsky Jr, N. Lippmann, and S.; Schmauder, Computational modelling of crack propagation in real microstructures of steels and virtual testing of artificially designed materials, ''International Journal of Fracture'' '''120''' (2003), 581 –600.
481
482
<div id="cite-17"></div>
483
[[#citeF-17|[17]]] J. Lopez, S. Oller, E. Oñate, and J. Lubliner, A homogeneous constitutive model for masonry, ''Int. J. Numer. Meth. Engng''. '''46''' (1999), 16511671.
484
485
<div id="cite-18"></div>
486
[[#citeF-18|[18]]] O.C Zienkiewicz and J.Z. Zhu, The superconvergent patch recovery (spr) and adaptive finite element refinement., ''Comput. Methods App. Mech. Engng''. '''101''' (1992), 207 –224.
487
488
<div id="cite-19"></div>
489
[[#citeF-19|[19]]] E. Oñate “Desarrollos y aplicaciones de modelos de fractura en la Escuela de Ingenieros de Caminos de Barcelona'''', ''Centro Internacional de Métodos Numéricos en Ingeniería, CIMNE ''(2000).
490
491
<div id="cite-20"></div>
492
[[#citeF-20|[20]]] C. Labra and E. Oñate, High-density sphere packing for discrete element method simulations, ''Commun. Numer. Meth. Engng''. '''25''' (2009), no. 7, 837 –849.
493
494
<div id="cite-21"></div>
495
[[#citeF-21|[21]]] M. Luong, Tensile and shear strengths of concrete and rock, ''Engineering Fracture Mechanics'' '''1-3''' (1990), no. 35, 127–135.
496
497
<div id="cite-22"></div>
498
[[#citeF-22|[22]]] Astm standard d638 - 10, 2003, “standard test method for tensile properties of plastics,” astm international, west Conshohocken, PA, 2003, doi. 10.1520/d0638-10, [http://www.astm.org www.astm.org].
499
500
<div id="cite-23"></div>
501
[[#citeF-23|[23]]] F. L. L. B. Carneiro, A new method to determine the tensile strength of concrete, ''Proceedings of the 5th meeting of the brazilian association for technical rules.'', 1943, pp. 126 –129. (In Portuguese).
502
503
<div id="cite-24"></div>
504
[[#citeF-24|[24]]] C.G. Rocco, G.V. Guinea Tortuero, J. Planas Roselló y M. Elices Calafat. Efecto del tamaño de probeta sobre la resistencia a la tracción medida con el ensayo brasileño. ''Hormigón y acero'' '''204''' (1997), 47-63
505
506
<div id="cite-25"></div>
507
[[#citeF-25|[25]]] J. Rodríguez del Viso, J.R. Carmona y G. Ruiz López. Efecto de la forma y el tamaño de la probeta en la resistencia a compresión en hormigón de alta resistencia. Hormigón y Acero 248 (2008), 77-86
508
509
<div id="cite-26"></div>
510
[[#citeF-26|[26]]] Norma a Instrucción de Hormigón Estructural (EHE-08)
511
[https://www.fomento.gob.es/MFOM/LANG_CASTELLANO/ORGANOS_COLEGIADOS/MASORGANOS/CPH/instrucciones/EHE_es/ https.//www.fomento.gob.es/MFOM/LANG_CASTELLANO/ORGANOS_COLEGIADOS/MASORGANOS/CPH/instrucciones/EHE_es/]
512
513
<div id="cite-27"></div>
514
[[#citeF-27|[27]]] ''UNE-EN 12390-3.2009 Ensayos de hormigón endurecido. Parte 3. Determinación de la resistencia a compresión de probetas''. (2009) Asociación Española de Normalización y Certificación AENOR
515
516
<div id="cite-28"></div>
517
[[#citeF-28|[28]]] Z. P. Bazant, J. Planas, ''Fracture and Size Effect in Concrete and Other Quasibrittle Materials ''CRC Press LLC (1998). ISBN 9780849382840
518
519
<div id="cite-29"></div>
520
[[#citeF-29|[29]]] Y.D. Murray, A Abu-Odeh and R. Bligh., ''Evaluation of LS-DYNA Concrete Material Model 159 Publication No. FHWA-HRT-05-063'', (2007) US Department of Transportation.
521
522
<div id="cite-30"></div>
523
[[#citeF-30|[30]]] Ground accelerogram from el-centro, imperial valley irrigation district (comp s00e). [http://www.eng.ucy.ac.cy/petros/earthquakes/eq1.txt http.//www.eng.ucy.ac.cy/petros/earthquakes/eq1.txt].
524
525
<div id="cite-31"></div>
526
[[#citeF-31|[31]]] J. M. González, F. Zárate, E. Oñate  [https://link.springer.com/article/10.1007/s40571-017-0174-3 Pulse fracture simulation in shale rock reservoirs: DEM and FEM–DEM approaches]. ''Computational particle mechanics'' '''Online''' (2017), 1–17.
527
528
<div id="cite-32"></div>
529
[[#citeF-32|[32]]] E. Oñate, F. Zárate, M. A. Celigueta, J. M. González, J. Miquel, J. M. Carbonell, F. Arrufat, S. Latorre, M. Santasusana, [https://link.springer.com/chapter/10.1007/978-3-319-60885-3_15 Advances in the DEM and Coupled DEM and FEM Techniques in Non Linear Solid Mechanics] ''Advances in Computational Plasticity,'' Springer (2018), 309-335
530

Return to Onate Zarate 2019a.

Back to Top

Document information

Published on 01/01/2018

DOI: 10.1016/j.hya.2018.05.002
Licence: CC BY-NC-SA license

Document Score

0

Times cited: 1
Views 72
Recommendations 0

Share this document