0% encontró este documento útil (0 votos)
63 vistas108 páginas

Tutorial: Análisis Computacional de Problemas Termo-Mecánicos en El Hormigón

Descargar como pdf o txt
Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1/ 108

TUTORIAL

Análisis computacional de problemas


termo-mecánicos en el hormigón

Universidad Nacional de La Plata


Facultad de Ingeniería

Ing. Rafael Lamattina


Año 2014
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

1 OBJETIVOS
Este tutorial se desarrolla con el fin de transmitir los conocimientos
básicos para el análisis computacional de problemas termo-mecánicos
relacionados al hormigón. De esta forma se pretende realizar un aporte al
diseño seguro de las estructuras de hormigón sometidas a este tipo de
fenómenos.

En particular se desarrolla un ejemplo de un problema típico para poder


utilizar las herramientas proporcionadas por el código comercial ABAQUS de
análisis por el Método de los Elementos Finitos (FEM en inglés).

Se considera que el lector posee un manejo básico del programa y


conocimientos acerca de las propiedades y características más importantes del
hormigón, previamente a la realización de este tutorial.

2 CONCEPTOS PREVIOS
2.1 Fisuración de origen térmico
El hormigón es un material compuesto fabricado a partir del cemento, el
agua, agregados gruesos y finos, y en general se le incorporan aditivos para
cumplir ciertos requisitos en estado fresco o endurecido.

Al entrar en contacto los componentes del cemento y el agua durante el


mezclado, se producen las denominadas reacciones de hidratación, que tienen
la característica de ser exotérmicas. En consecuencia, durante el proceso de
fraguado y endurecimiento del hormigón, se genera calor debido a las
reacciones químicas y aumenta la temperatura en el material.

Las estructuras se encuentran en contacto con los fluidos del ambiente


como el aire, agua, etc., que en general se encuentran a menor temperatura en
esta etapa temprana del material. Debido a esto, las zonas externas de la
estructura alcanzarán una menor temperatura que las zonas internas,
generándose así gradientes térmicos entre ambos puntos de la pieza.

2
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Gráfico 1: Esquema del “salto térmico” que se produce en las estructuras de hormigón masivo.

Los gradientes térmicos se traducen en deformaciones, y estas últimas


en tensiones. Según sea el grado de restricción de la estructura, se pueden
alcanzar valores tensionales que superen la resistencia a tracción del material,
generándose de esta forma las fisuras que tienen origen en el comportamiento
térmico del hormigón.

Las fisuras en general se originan en la zona del “núcleo” de la


estructura de hormigón masivo (MSC por sus siglas en inglés: Mass Concrete
Structure), o en la superficie de la misma. Estos dos fenómenos de fisuración
se denominan gradiente de masa y gradiente de superficie respectivamente.

Además se debe mencionar el hecho de que este material evoluciona de


un estado fresco donde la mezcla tiene cierto grado de fluidez, a un estado
endurecido una vez finalizado el proceso de fraguado. Esto implica que las
propiedades mecánicas varían fuertemente durante dicho proceso lo que tiene
fundamental importancia en los fenómenos de fisuración en edad temprana. En
particular, es interesante analizar la evolución de la resistencia a tracción (que
depende directamente de la ganancia de resistencia a compresión), el módulo
de elasticidad, etc.

Gráfico 2: Ganancia de resistencia a compresión con la edad.

3
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

2.2 Hormigón masivo


El American Concrete Institute (ACI) define al hormigón masivo como
todo volumen de hormigón cuyas dimensiones son lo suficientemente grandes
como para que sea necesario tomar medidas en relación a la generación del
calor de hidratación del cemento y a los consiguientes cambios de volumen de
forma tal de minimizar la fisuración.

2.2.1 Estructuras de hormigón masivo

De la definición anterior, se desprende que este tipo de estructuras son


de gran envergadura y gran parte de su resistencia es aportada por su masa.
De acuerdo a esto, no se trata de hormigón armado convencional ya que no
existen armaduras de acero en toda la estructura y las cuantías pueden ser
bajas.

Entre las estructuras de hormigón masivo se pueden nombrar: muros de


gravedad (contención de suelos o líquidos), fundaciones de turbinas
termoeléctricas o maquinaria que transmite cargas dinámicas a la estructura,
presas de hormigón convencional, presas de Hormigón Compactado a Rodillo
(HCR), losas de espesor y extensión importantes, pilas de gran diámetro en
puentes, muertos de anclaje de cables en puentes colgantes, etc.

Es importante entonces, realizar un estudio térmico de este tipo de


estructuras, ya que las fisuras pueden tener consecuencias graves como
generar filtraciones, afectar la durabilidad, comprometer el comportamiento
dinámico de las estructuras en servicio, etc.

Figura 1: Galería de inspección de una presa con filtraciones.

4
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 2: Muro de contención con filtraciones.

2.2.1.1 Presas de Hormigón Compactado a Rodillo

Las presas HCR fueron ampliamente utilizadas en las últimas décadas.


Para su construcción se utiliza un hormigón de bajo asentamiento, de
resistencia media/baja, y en general el cemento posee cantidades importantes
de adiciones minerales para disminuir el calor de hidratación.

La metodología constructiva consiste en realizar tongadas (capas de


hormigón) de pequeño espesor, aproximadamente 0,30 m y luego se compacta
el material utilizando maquinaria de uso vial.

Figura 3: Construcción de una presa HCR.

Los parámetros de diseño con los que puede trabajar el proyectista para
disminuir los problemas térmicos de este tipo de presas, son entre otros:

• El espesor de las tongadas de hormigón.


• El intervalo de tiempo que transcurre entre dos tongadas
consecutivas.
• La temperatura de colocación del hormigón.

5
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

• El período de hormigonado (en relación con el clima en la zona de


emplazamiento y la temperatura ambiente).
• La composición del cemento con el cual se fabrica el hormigón.
• Enfriamiento mediante colocación de conducciones con líquidos
refrigerantes.
• Ejecución de juntas de expansión/contracción.

3 HERRAMIENTA DE CÁLCULO
Como ya se explicó en el apartado OBJETIVOS, para resolver los
problemas planteados se utiliza el software ABAQUS v6.11.

A continuación se desarrollarán ciertos aspectos que es importante tener


en cuenta antes de abordar el problema particular relacionado con este tutorial.

3.1 Ecuaciones de Transferencia de Calor


En este apartado se presenta el desarrollo matemático que involucra
tanto las ecuaciones que describen el fenómeno de transferencia de calor,
como la formulación necesaria de las mismas para aplicar el FEM.

3.1.1 Relación constitutiva

Las observaciones experimentales indican que el flujo de calor en un


material es proporcional al gradiente de temperaturas, siendo el factor de
proporcionalidad la conductividad térmica del material.

 ∂θ
q = −K ⋅ Ec.(1)
∂x

Siendo q el vector de flujo de calor, K la matriz de conductividad y θ las


temperaturas.

3.1.2 Ecuación de balance térmico -Forma Diferencial-

Se parte de una formulación Lagrangiana del problema, analizando un


volumen arbitrario de un sólido como el que se muestra a continuación:

6
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 4: Volumen de control del material en estudio.

En la figura anterior, se indica el volumen de control con la letra V y la


superficie de control con la letra S. Se transcribe a continuación el balance de
energía térmica en el cuerpo:

∫ ρ ⋅ U ⋅ dV = ∫ q ⋅ dS + ∫ r ⋅ dV
V S V
Ec. (2)

En esta ecuación q es el flujo de calor que ingresa/egresa a través de la


superficie de control, r es el calor generado en el volumen de control, Ů es la
variación temporal de la energía interna, y ρ es la densidad del material.

Para obtener la forma diferencial de la ecuación de balance térmico,


primero se reemplaza el flujo de calor superficial por su forma de producto
escalar de dos vectores. El signo negativo indica un flujo de calor saliente:

q = −q • nˆ Ec. (3)

Reemplazando la ec. (3) en la ec. (2) de balance térmico:


• 
∫ ⋅ dV = − ∫ q • n̂ ⋅ dS + ∫ r ⋅ dV
ρ ⋅U Ec. (4)
V S V

Para simplificar el desarrollo matemático, sin quitarle generalidad al


planteo, se considerará de aquí en más un problema unidireccional de
transferencia de calor.

Aplicando el Teorema de la Divergencia en el primer término a la


izquierda de la igualdad en la ec. (4):

 ∂q
∫ q • nˆ ⋅ dS = ∫ ∂x ⋅ dV
S V
Ec. (5)

Reemplazando la ec. (5) en la ec. (4):

7
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

• ∂q
∫ ρ ⋅ U ⋅ dV = −∫
V V
∂x
⋅ dV + ∫ r ⋅ dV
V
Ec. (6)

Dado que el volumen es arbitrario, se puede derivar la expresión


anterior, obteniéndose:
• ∂q
ρ ⋅U = − +r Ec. (7)
∂x

Esta última es la forma diferencial de la ecuación de balance térmico.

3.1.3 Ecuación de balance térmico -Forma Variacional-

La expresión diferencial vista anteriormente (denominada “fuerte”), no


puede ser integrada numéricamente. Para ello se debe “relajar” la ecuación de
la siguiente forma:

Se multiplica la ec. (6) por una temperatura variacional arbitraria δθ


(válida para toda variación admisible compatible) y se integra en el volumen:
• ∂q
∫ δθ ⋅ ρ ⋅U ⋅ dV = −∫ δθ ⋅
V V
∂x
⋅ dV + ∫ δθ ⋅ r ⋅ dV
V
Ec. (8)

Operando con el primer término a la derecha de la igualdad:

∂q ∂q ∂δθ ∂δθ
− δθ ⋅ = −δθ ⋅ − q ⋅ + q⋅ Ec. (9)
∂x ∂x ∂x ∂x

Aplicando la regla de la cadena en los dos primeros términos a la


derecha de la igualdad en la ecuación anterior:

∂q ∂(q ⋅ δθ ) ∂δθ
− δθ ⋅ =− +q⋅ Ec. (10)
∂x ∂x ∂x

Reemplazando la ec. (10) en la ec. (8):

• % ∂(q ⋅ δθ ) ∂δθ (
∫ δθ ⋅ ρ ⋅U ⋅ dV = ∫ '&− ∂x
+q⋅ * ⋅ dV + ∫ δθ ⋅ r ⋅ dV
∂x )
Ec. (11)
V v V

Aplicando el Teorema de la Divergencia:

∂(q ⋅ δθ ) 
∫− ⋅ dV = ∫ −nˆ • q ⋅ δθ ⋅ dS Ec. (12)
V
∂x S

Convirtiendo el flujo de calor superficial a la forma escalar dada por la


ec. (3):

8
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

∫ −n̂ • q ⋅ δθ ⋅ dS = ∫ q ⋅ δθ ⋅ dS Ec. (13)


S S

Reemplazando la expresión (13) en la ec. (11):


• ∂δθ
∫ δθ ⋅ ρ ⋅ U ⋅ dV = ∫ q ⋅ δθ ⋅ dS + ∫ q ⋅
V S V
∂x
⋅ dV + ∫ δθ ⋅ r ⋅ dV
V
Ec. (14)

Reacomodando los términos de la ecuación anterior:


• ∂δθ
∫ δθ ⋅ ρ ⋅ U ⋅ dV − ∫ q ⋅
V V
∂x
⋅ dV = ∫ q ⋅ δθ ⋅ dS + ∫ δθ ⋅ r ⋅ dV
S V
Ec. (15)

Esta última es la forma variacional de la ecuación de balance térmico.

3.1.4 Aproximación a través del Método de los Elementos


Finitos

Este método se basa en la discretización especial y temporal del


problema y se procede integrando numéricamente las ecuaciones diferenciales
que describen matemáticamente el fenómeno físico.

Para ello, se parte de la expresión débil hallada anteriormente y se


discretiza el volumen en elementos finitos, de forma tal que la temperatura en
cualquier punto del material se calcula utilizando funciones de interpolación:

θ = N N ⋅θ N
Ec. (16)

N
Donde N N (x) son las funciones de interpolación y θ las temperaturas
nodales.

Utilizando la aproximación de Galerkin, las temperaturas variacionales


también son interpoladas a partir de las mismas funciones:

δθ = N N ⋅ δθ N
Ec. (17)

Introduciendo las interpolaciones (16) y (17) en la forma variacional de la


ecuación de balance térmico (15), resulta:


∂(N N ⋅ δθ N )
∫ N N ⋅ δθ N ⋅ ρ ⋅U ⋅ dV − ∫ q ⋅ = ∫ q⋅ N N
⋅ δθ N ⋅ dS + ∫N N
⋅ δθ N ⋅ r ⋅ dV
V V ∂x S V

Ec. (18)

9
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Reemplazando en el segundo término a la izquierda de la igualdad, la


expresión del flujo de calor “q” dado por la Ley de Fourier ec. (1):

 ∂θ ∂(N M ⋅ θ M )
q = −K ⋅ = −K ⋅ Ec. (19)
∂x ∂x

Introduciendo la ec. (19) en la ec. (18) y reacomodando los términos:

⎧ • ∂N M ∂N N ⎫ ⎧ ⎫
δθ N ⋅ ⎨∫ N N ⋅ ρ ⋅ U ⋅ dV + θ M ⋅ ∫ K ⋅ ⋅ ⋅ dV ⎬ = δθ N ⋅ ⎨∫ N N ⋅ r ⋅ dV + ∫ N N ⋅ q ⋅ dS ⎬
⎩V V
∂x ∂x ⎭ ⎩V S ⎭

Ec. (20)

Como las temperaturas variacionales son arbitrarias, se pueden


simplificar obteniéndose como resultado:

N

M ∂N M ∂N N
∫ N ⋅ ρ ⋅U ⋅ dV + θ ⋅∫K ⋅ ⋅ ⋅ dV = ∫ N N ⋅ r ⋅ dV + ∫ N N ⋅ q ⋅ dS
V V
∂x ∂x V S
Ec. (21)

Esta ecuación es la base para la resolución de problemas termo-


mecánicos desacoplados unidimensionales mediante elementos finitos.
Cuando la energía interna del sólido no depende del campo de deformaciones
o de tensiones, se puede escribir:
• dU dU dθ dU • •
U= = ⋅ = ⋅θ = c ⋅θ Ec. (22)
dt dθ dt dθ

Siendo c el calor específico del material.

Considerando la ec.(21), la ec. (22) puede escribirse en forma matricial:



C ⋅θ + K ⋅θ = Q Ec. (23)

Donde C, es la matriz de capacidad calórica:

C NM = ∫ ρ ⋅ c ⋅ N N ⋅ N M ⋅ dV Ec. (24)
V

K es la matriz conductividad, definida de la siguiente forma:

∂N M ∂N N
K NM = ∫ K⋅ ⋅ ⋅ dV Ec. (25)
V ∂x ∂x

Q es el vector de flujo de calor:

10
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Q N = ∫ N N ⋅ q ⋅ dS + ∫ N N ⋅ r ⋅ dV Ec. (26)
S V

y θ es el vector de temperaturas nodales.

Para condiciones estacionarias se tiene que Ů = 0, con lo cual el primer


término de la izquierda en la ec. (21) se anula.

Para problemas transitorios, además de espacialmente, se debe integrar


la ecuación en forma temporal. Esto se consigue utilizando el siguiente
algoritmo de diferencias regresivas:

• U t + Δt − U t
U t + Δt = Ec. (27)
Δt

Reemplazando este algoritmo (27) en la ecuación (21) resulta:

1 N $• • '
M ∂N M ∂N N
∫N ⋅ ρ ⋅ &U t+Δt −Ut )dV + θ t+Δt ⋅ ∫ kt+Δt ⋅ ⋅ ⋅ dV − ∫ N N ⋅ rt+Δt ⋅ dV + ∫ N N ⋅ qt+Δt ⋅ dS = 0
Δt V
% ( V ∂x ∂x V S

Ec. (28)

Este es un sistema de ecuaciones no lineal que se resuelve a través del


Método de Newton modificado. Se denomina modificado porque la matriz del
Jacobiano no se halla en forma exacta.

Una vez obtenido el campo de temperaturas nodales, se puede obtener


el campo de deformaciones del sólido a partir del coeficiente de expansión
térmica del material α:

ε th = α ⋅ (θ − θ I ) Ec. (29)

Donde θ es la temperatura en el tiempo t + ∆t, y θI es la temperatura


inicial en el tiempo t.

Las deformaciones elásticas del material se obtienen por diferencia entre


las deformaciones totales y las deformaciones de origen térmico:

ε el = ε − ε th Ec. (30)

Conocido el campo de deformaciones, las tensiones pueden calcularse


utilizando la relación constitutiva de Hooke:

σ = E ⋅ ε el Ec. (31)

11
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

3.2 Tipos de análisis


El software cuenta con dos procedimientos de cálculo que tienen
características diferentes y que serán explicadas a continuación.

3.2.1 Fully Coupled Thermal-Stress Analysis (FC)

Este tipo de análisis se lleva a cabo cuando el campo de tensiones


depende del campo de temperaturas y a la vez, el campo de temperaturas
depende del campo de tensiones:

σ = f (T )∧T = f (σ ) Ec. (32)

Estos problemas se presentan por ejemplo en el caso de deformaciones


plásticas en piezas metálicas frente a cargas dinámicas.

La ventaja que presenta es que sólo es necesario realizar un único


modelo en el cual se incorporan las propiedades térmicas y mecánicas del
hormigón, las condiciones de borde también térmicas y mecánicas, y el
programa resuelve en cada paso de cálculo el campo de temperaturas y
tensiones de la estructura.

Como desventaja, se debe mencionar que el costo computacional es


mucho mayor, dado que el tiempo de cálculo y el tamaño de los archivos crece
exponencialmente con la cantidad de nodos y elementos y en mayor medida
que el segundo tipo de análisis que se mencionará a continuación.

Con lo cual este tipo de procedimientos es recomendable para la


resolución de problemas no demasiado complejos, bidimensionales, con
geometrías simples, mallas “gruesas”, etc.

3.2.2 Sequentially Coupled Thermal-Stress Analysis (SC)

Este tipo de análisis corresponde a los problemas en donde el campo de


tensiones depende del campo de temperaturas, pero este último no depende
del primero:

σ = f (T )∧T ≠ f (σ ) Ec. (33)

La ventaja de los procedimientos SC respecto a los FC, radica


principalmente en el costo computacional, ya que los archivos que genera el

12
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

programa ocupan mucho menor espacio y el cálculo se resuelve más


rápidamente.

La desventaja frente a los análisis FC es que en este caso, es necesario


realizar dos modelos por separado. Uno de ellos es el modelo térmico en
donde se introducen las propiedades térmicas del material y las condiciones de
borde térmicas. En este primer análisis el programa calcula las temperaturas
nodales y éstas luego deben ser incorporadas como input en un modelo
mecánico. En este último se deben incorporar las propiedades mecánicas del
hormigón y las condiciones de borde mecánicas de la estructura y el programa
calcula como output las tensiones y deformaciones en los elementos.

En conclusión, este tipo de análisis presenta un preproceso más


complicado pero es aconsejable en aquellos problemas que sean más
complejos, con una gran cantidad de nodos y elementos.

3.3 Análisis dimensional


Para reducir el costo computacional de un modelo numérico, siempre es
conveniente realizar un análisis bidimensional del problema físico. Pero al
adoptar una geometría 2D frente a una 3D (que es más exacta respecto la
realidad del fenómeno), se debe asegurar el cumplimiento de las hipótesis de
la elasticidad plana.

Por ejemplo en el caso en que una dimensión de la estructura sea


mucho mayor que las otras dos y que las cargas sean perpendiculares a dicha
dirección, se puede realizar un análisis bidimensional a través de la teoría de
deformaciones planas:

! σ xx σ xy 0 $
# &
[T ] XYZ = # σ yx σ yy 0 &
# &
#" 0 0 υ (σ xx + σ yy ) &
%
! ε xx (x, y) ε xy (x, y) 0 $
# &
[ D] XYZ = # ε yx (x, y) ε yy (x, y) 0 &
# &
#" 0 0 0 &%

Figura 5: Hipótesis de la deformación plana. Ec. (34)

13
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4 PRIMERA PARTE: ANÁLISIS LINEAL


En esta primera parte del tutorial, se resolverá paso a paso un problema
termo-mecánico de una estructura de hormigón masivo utilizando un análisis
de primer orden. Esto implica que la relación entre tensiones y deformaciones
es lineal elástica y que las tensiones obtenidas como output, pueden superar la
resistencia a tracción/compresión del material.

En este tipo de análisis no se pueden obtener resultados en cuanto a las


posibilidades y características de propagación de una fisura de origen térmico.
Pero por otro lado, se puede observar en qué puntos de la estructura y en qué
instantes de tiempo es de esperarse que se supere la resistencia a tracción del
material y aparezca una fisura. Es decir, que estos modelos aportan una idea
bastante precisa del campo de temperaturas en la estructura, y una idea
aproximada del campo de tensiones.

4.1 Planteo del problema


Se analizará el comportamiento termo-mecánico de un muro de
hormigón perteneciente a un sistema de esclusas. Las dimensiones de un
módulo de muro son de 2,00 m de ancho por 15,00 m de alto por 30,00 m de
largo.

El procedimiento constructivo de esta esclusa se basa en 3 etapas que


dividen la altura total del muro en “capas” de 5.00 m. Se considera que el
tiempo entre cada tongada es de 72 hs. A los 60 días de iniciada la
construcción, se considera que la esclusa está operativa, con lo cual se
procede al llenado de la misma. El tiempo total del análisis será de 100 días
(tiempo en el cual se considera que será posible alcanzar el equilibrio térmico
de la estructura). Se considera que el hormigonado comienza en los meses
más fríos.

El material con el que se construye el muro es un hormigón de Cemento


Portland Normal y resistencia característica 35 MPa (H-35). El mismo está
fundado sobre un macizo de roca cuyas propiedades serán vistas en otro
apartado. Se considera empotrado en la zona inferior y en la zona lateral
existen anclajes dispuestos cada 1 m.

14
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 6: Imágenes del funcionamiento de las esclusas del canal de Panamá.

Se toma como temperatura ambiente los valores medios mensuales


correspondientes a la ciudad de Panamá. Para la temperatura del agua se
tienen en cuenta las medias mensuales de la Temperatura de la Superficie del
Mar (TSM) correspondientes a Isla Grande. La temperatura de colocación es la
misma para las tres tongadas de hormigón.

El modelo a desarrollar es 2D y el procedimiento de cálculo SC.

Figura 7: Perspectiva y vista del módulo de muro considerado.

4.2 Administración de archivos


Es conveniente tener un criterio unificado para nombrar los archivos que
hacen referencia a un mismo trabajo, con el objetivo de tener un control preciso
de la información y de las modificaciones que se realicen. De la misma forma,
se recomienda que todos estos archivos estén contenidos en una misma
carpeta.

En este caso se propone una forma de nomenclatura, pudiendo el lector


adoptar otra, siempre que sea consistente y completa.

15
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Lo primero que se debe hacer al abrir el programa es guardar el archivo


en un directorio cuya ruta sea fácil de recordar. El nombre del archivo debe
contener una referencia al trabajo en sí, y se puede agregar la fecha o la
versión o el tipo de análisis realizado. En este caso se nombrará el archivo de
la siguiente forma:

MOD ESC LIN.cae

(MOD: Modelo, ESC: Esclusa, LIN: Lineal)

En la parte superior de la barra de herramientas del programa puede


leerse la ruta del directorio utilizado.

En el caso de contar con numerosos archivos CAE se puede generar un


archivo en Excel en donde se indiquen las características más importantes de
los modelos y qué modificaciones se van realizando. Por ejemplo:

REF MOD ESC.xlsx

(REF: Referencia, MOD: Modelos, ESC: Esclusa)

4.3 Unidades
Como es sabido, el programa no cuenta con unidades de medida
programadas, con lo cual es de vital importancia que el usuario trabaje durante
el preproceso con unidades consistentes. Esto implica que se utilicen siempre
las mismas unidades al introducir el input en el modelo: por ejemplo si se
ingresa la densidad en kg/m3 no se puede ingresar el módulo de elasticidad en
MPa ya que ambas unidades no son consistentes.

En este caso se trabajará con las siguientes unidades:

MODELO+TÉRMICO MODELO+MECÁNICO
ENTIDAD UNIDAD ENTIDAD UNIDAD
Masa kg Peso N
Longitud mm Longitud mm
Tiempo s Tiempo s
Calor kcal Temperatura ºC
Temperatura ºC

Tabla 1: Sistema de unidades utilizado en el tutorial.

16
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4.4 Modelo Térmico


En primer lugar se asigna el nombre del modelo, en este caso se
propone:

MOD TER ESC 001

(MOD:Modelo, TER: Térmico, ESC: Esclusa, 001: Versión).

Al igual que otros dos archivos en Excel para manejar la información de


entrada y salida del programa:

INPUT MOD ESC.xlsx

(INPUT: Datos de entrada, MOD: Modelo, ESC: Esclusa).

OUTPUT MOD ESC.xlsx

(OUTPUT: Datos de salida, MOD: Modelo, ESC: Esclusa).

4.4.1 Input

4.4.1.1 Parts

A continuación se presenta un esquema de la geometría 2D del análisis


térmico. Se pueden observar las tres partes que corresponden al procedimiento
constructivo del muro (M1,M2,M3), y una cuarta parte que corresponde al
terreno de fundación (T). Las dimensiones dadas al terreno que rodea a la
estructura suelen ser 0.5L a 2L, siendo L una dimensión de la estructura.

Figura 8: Partes que conforman el modelo térmico.

17
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Primero, se generan las 3 partes correspondientes al muro:

Parts à Name: M1à Modeling space: 2D Planarà Type: Deformableà Base Feature:

Shellà Continue…à Create Lines (Rectangle) à Pick a starting corner for the
rectangle -- or enter X,Y: 0,0à Enterà Pick the opposite corner for the rectangle -- or

enter X,Y: 2000,5000à Cancel the procedure à Done

Lo mismo se debe realizar para las partes M2 y M3, teniendo en cuenta


que las coordenadas a introducir son las que figuran en la siguiente tabla:

STARTING"CORNER OPPOSITE"CORNER
PART
X"[mm] Y"[mm] X"[mm] Y"[mm]
M1 0,00 0,00 2000,00 5000,00
M2 0,00 5000,00 2000,00 10000,00
M3 0,00 10000,00 2000,00 15000,00

Tabla 2: Coordenadas de las partes correspondientes al muro.

Figura 9

18
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 10

Luego se debe crear la parte correspondiente al terreno que rodea al


muro:

Parts à Name: Tà Modeling space: 2D Planarà Type: Deformableà Base Feature:

Shellà Continue…à Create Line (Connected) à Pick a starting point for the line -- or
enter X,Y: 0,0à Pick an end point for the line -- or enter X,Y: 4500,0à …à Cancel the

procedure à Done

Donde se indica …, se deben seguir ingresando las coordenadas de los


puntos que conforman la parte:

PUNTO X'[mm] Y[mm]

1"(Inicio) 0,00 0,00


2,00 4500,00 0,00
3,00 4500,00 02500,00
4,00 02500,00 02500,00
5,00 02500,00 15000,00
6,00 0,00 15000,00
7"(Fin) 0,00 0,00

Tabla 3: Coordenadas de la parte correspondiente al terreno de fundación.

19
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 11

4.4.1.2 Materials

Conductividad térmica (k):

Es una propiedad física de los materiales que mide la capacidad de


conducción de calor. Cuando se calienta la materia, la energía cinética
promedio de sus moléculas aumenta, incrementándose su nivel de agitación.
La conducción de calor, que a nivel macroscópico puede modelizarse mediante
la ley de Fourier, a nivel molecular se debe a la interacción entre las moléculas
que intercambian energía cinética sin producir movimientos globales de
materia. Por tanto la conducción térmica difiere de la convección térmica en el
hecho de que en la primera no existen movimientos macroscópicos de materia.
En el hormigón esta propiedad depende del tamaño y volumen de los poros, de
la conductividad térmica del agregado, etc.

Para estos casos, se considera un valor constante para todas las


edades y para todas las temperaturas alcanzadas en el material debidas al
calor de hidratación. El ACI proporciona un método para el cálculo de esta
propiedad: CRD-C 44-63 “Method for calculation of termal conductivity of
concrete”, (basado en el CRD-C 37-73 y el CRD-C 124-73).

20
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Calor específico (c):

Es una magnitud física que se define como la cantidad de calor que hay
que suministrar a la unidad de masa de una sustancia o sistema termodinámico
para elevar su temperatura en una unidad.

Para el hormigón, esta propiedad también es considerada constante. El


ACI proporciona el siguiente artículo para su determinación: CRD-C 124-73
“Method of test for specific heat of aggregates, concrete, and other materials
(Method of Mixtures)”.

Densidad (δ):

Es una magnitud escalar referida a la cantidad de masa en un


determinado volumen de una sustancia. Depende también del tipo y proporción
de agregados utilizados. El método indicado para su determinación es el ASTM
C138/C138M-09 “Standard Test Method for Density (Unit Weight), Yield, and
Air Content (Gravimetric) of Concrete”.

A continuación se presenta la forma de ingresar las propiedades


térmicas de los materiales en el programa:

VALOR
PROPIEDAD SÍMBOLO MACIZO: UNIDAD
HORMIGÓN
ROCOSO
Densidad ρ 2,51E'06 2,60E'06 kg/mm3
Conductividad:térmica k 5,79E'07 5,16E'07 kcal/mm.s.ºC
Calor:específico c 0,23 0,22 kcal/kg.ºC

Tabla 4: Propiedades térmicas correspondientes al input.

Materials à Name: HORMIGON

Generalà Densityà Distribution: Uniformà Number of field variables: 0à Mass Density:
2,51E-06

Thermalà Conductivityà Type: Isotropicà Number of field variables: 0à Conductivity:


5,79E-07

Thermalà Specific Heatà Type: Constant Volumeà Number of field variables:


0à Specific Heat: 0,23

Lo mismo se debe hacer para el material correspondiente al macizo


rocoso del terreno de fundación: MACIZO.

21
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 12

4.4.1.3 Sections

A continuación es necesario generar las secciones correspondientes a


cada material denominadas:

SEC HOR

(SEC: Sección, HOR/MAC: Hormigón/Macizo).

Section à Name: SEC HORà Category: Solidà Type: Homogeneousà Continue…à

Material: HORMIGONà Ok

Figura 13

22
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 14

Lo siguiente es asignar a cada parte la sección correspondiente:

Parts à M1à Section Assigments à Seleccionar parteà Doneà Section: SEC


HORà Thickness, Assigment: From Sectionà Ok

Se debe repetir esta última operación para las 4 partes que conforman el
modelo térmico.

Figura 15

23
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 16

4.4.1.4 Assembly

Se deben ensamblar las partes creadas para poder generar más


adelante las interacciones entre las mismas:

Assembly à Instances à M1à Instance Type: Independent (mesh on


instance)à Apply

Figura 17

24
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4.4.1.5 Steps

Para dividir el modelo en steps o etapas se deben considerar cambios


en las condiciones de borde, cambios en las cargas actuantes, cambios en la
geometría, etc. En este caso se tienen en principio 3 steps correspondientes a
cada tongada de hormigón (cambia la carga térmica y la geometría). Luego se
tiene un cuarto step debido a que cambia una condición de borde, (la superficie
lateral derecha del muro pasa de estar en contacto con el aire, a estar en
contacto con el agua). Como aproximación se considera que el llenado de la
esclusa es instantáneo.

A fin de organizar la duración de cada step, se confecciona la siguiente


tabla:

Duración.
Duración.........Duración. Suceso...... Suceso.......
STEP ti.[s] tf.[s] Acumulada.
[s] [días] Inicio Fin
[días]
Step%1 0.00 259200.00 259200.00 3.00 3.00 M1 M2
Step%2 259200.00 518400.00 259200.00 3.00 6.00 M2 M3
Step%3 518400.00 5184000.00 4665600.00 54.00 60.00 M3 Cond.0Borde
Step%4 5184000.00 8640000.00 3456000.00 40.00 100.00 Cond.0Borde Eq.0Térmico

Tabla 5: Organización del modelo en steps.

Un concepto que se debe tener en cuenta es el rango de incremento de


tiempo que se debe introducir en el programa para que el mismo realice las
iteraciones. Este incremento no debe ser muy alto como para generar
problemas de convergencia, y no debe ser muy pequeño de modo de que la
cantidad de pasos de cálculo para completar el tiempo total sea excesiva.

Si se utilizan elementos de segundo orden (no es este caso) existe un


límite inferior dado por la siguiente ecuación:

ρ ⋅c
Δt ≥ ⋅ Δl 2 Ec. (35)
6⋅k

- ∆t es el mínimo incremento de tiempo que se puede aplicar en el paso


de cálculo.

- ρ, c, k son las propiedades del material ya vistas

- ∆l es la distancia entre nodos de un elemento en donde exista el mayor


gradiente de temperatura.

25
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

El límite superior se debe introducir en función de los requerimientos


particulares del problema. Es decir, si en el OUTPUT es necesario contar con
las temperaturas hora por hora, el incremento máximo será de 3600 segundos.
En nuestro caso se tomará 1 día como máximo incremento.

También se debe introducir el parámetro ∆Tmáx que es el máximo


cambio de temperatura que se admite para un paso de cálculo.

El usuario deberá ir modificando estos tres parámetros hasta lograr la


convergencia del análisis numérico. Es conveniente comenzar con los valores
otorgados por defecto en el programa.

Steps à Name: Step-1à Insert new step after: Initialà Procedure type: Generalà Heat
Transferà Continue…

Basicà Description: …à Response: Transientà Time period: 259200

Incrementationà Type: Automaticà Maximum number of increments: 1000à Increment


size. Initial: 0,001. Minimum: 0,001. Maximum: 86400à Max. allowable temperature
change per increment:12

à Ok.

Figura 18

26
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 19

Figura 20

El mismo procedimiento se debe realizar para los demás steps, teniendo


en cuenta la duración de cada uno.

27
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4.4.1.6 Interactions

Convección Hormigón-Aire:

En la superficie de la estructura ocurre la transferencia de calor entre el


material y los fluidos del ambiente. Esta interacción está gobernada por el
denominado coeficiente de convección entre el hormigón y el aire (h), que
depende de múltiples factores como la rugosidad de la superficie, la velocidad
del aire, la temperatura del aire, etc. En este caso se considera un valor
constante:

PARÁMETRO SÍMBOLO VALOR UNIDAD

Coeficiente6de6convección6Hº<Aire h 3,33E%09 kcal/mm2.s.ºC


Tabla 6

q = −h ⋅ (θ − θ∞ ) Ec. (36)

Figura 21

A continuación se presenta un esquema de la interacción Hº-Aire en


cada step. Se debe recordar que en el último step, la esclusa se encuentra
llena de agua:

Figura 22

En la siguiente tabla se resume la información que figura en el esquema.

28
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

INTERACCIÓN Step01 Step02 Step03 Step04

CONV%HA%M11 ✓ ✓ ✓ ✗
CONV%HA%M12 ✓ ✗ ✗ ✗
CONV%HA%M21 ✗ ✓ ✓ ✗
CONV%HA%M22 ✗ ✓ ✗ ✗
CONV%HA%M31 ✗ ✗ ✓ ✗
CONV%HA%M32 ✗ ✗ ✓ ✓

Tabla 7: Participación de las interacciones en cada step.

La nomenclatura elegida es la que sigue:

CONV HA M11

(CONV: Convección, HA: Hormigón-Aire, M1: Parte correspondiente,


1/2: Superficie lateral / Superficie superior).

En primer lugar se deben definir las superficies correspondientes a cada


interacción, que serán llamadas con el mismo nombre:

Assembly à Surfaces à Name: CONV HA M11à Continue…à Select regions for the
surface (Individually)à Done

Figura 23

29
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 24

Se debe repetir esta operación para todas las demás superficies.

NOTA: para poder eliminar del viewport algunas de las instances y poder
seleccionar las superficies deseadas existen los siguientes elementos:

Oculta todas las instances menos la seleccionada.

Oculta la instance seleccionada.

Muestra todas las instances.

Lo siguiente es generar la curva de amplitud correspondiente a la


temperatura del aire. Para ello se debe confeccionar una tabla en Excel en el
archivo INPUT MOD ESC.xlsx:

30
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Tmedia======
Mes
[ºC] Temperaturas"Medias"Mensuales"
Ene 28,00 (Panamá)"
Feb 29,00 30,00%
Mar 29,00 29,00%
Abr 28,50 28,00%
May 28,50

T"[ºC]"
27,00%
Jun 27,00 26,00%
Jul 26,00
25,00%
Ago 26,50
24,00%
Sep 26,00
23,00%
Oct 26,50 Ene% Feb% Mar% Abr% May% Jun% Jul% Ago% Sep% Oct% Nov% Dic%
Nov 26,00
Mes"
Dic 25,50

Tabla 8 Gráfico 3

AMP"TA
t""""""""""""""""""" Tamb""""""""""""
Tamb............ t...................
t.acumulado...................
t...................
t.acumulado................. [s] [ºC]
Mes
[ºC] [días] [días] [s] [s]
0.00 26.00
Sep 26.00 30.00 30.00 2592000.00 2592000.00 2592000.00 26.50
Oct 26.50 31.00 61.00 2678400.00 5270400.00 5270400.00 26.00
Nov 26.00 30.00 91.00 2592000.00 7862400.00 7862400.00 25.50
Dic 25.50 31.00 122.00 2678400.00 10540800.00 8640000.00 26.31

Tabla 9 Tabla 10: Curva de amplitud

Temperatura"Ambiente"durante"el"
hormigonado"
26.75$

26.50$

26.25$
T"[ºC]"

26.00$

25.75$

25.50$

25.25$
0.E+00$

1.E+06$

2.E+06$

3.E+06$

4.E+06$

5.E+06$

6.E+06$

7.E+06$

8.E+06$

9.E+06$

Tiempo"[s]"

Gráfico 4

La Tabla 10 es la que se debe ingresar en el programa. Para generar la


curva de amplitud se debe hacer lo siguiente:

Amplitudes à Name: AMP TAà Type: Tabularà Continue…à Time span: Total
timeà Smoothing: Use solver defaultà Amplitude Data: Copiar y pegar la tabla de
Excelà Ok

La opción Total time hace referencia a que el tiempo que se ingresa al


copiar la tabla de Excel, será interpretado en el programa como tiempo total del
modelo. Esto permite que la curva de amplitud se ingrese en un step y el
programa siga “leyendo” sus valores en los steps siguientes.

31
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

En cambio, si se ingresa la opción Step time, se estará indicando que los


tiempos de la curva de amplitud corresponden a tiempos del step. Con lo cual,
una vez que finaliza el step, el valor de la tabla de amplitud se mantiene
constante en los siguientes steps lo que sería incorrecto en este caso.

Figura 25

Figura 26

32
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 27

Lo siguiente será generar las interacciones correspondientes, con lo cual


se deberá hacer lo siguiente:

Interactions à Name: CONV HA M11à Step: Step-1 (step en el que “aparece” la


interacción)à Types for selected step: Surface film condictionà Continue…

Surfaces…à Name: CONV HA M11à Continue…

Definition: Embedded coefficientà Film coefficient: 3.33E-09à Film coefficient amplitude:


(Instantaneous)à Sink definition: Uniformà Sink temperature: 1à Sink amplitude: AMP
TA (Seleccionar la amplitud correspondiente a la temperatura ambiente generada en el
paso anterior)à Ok

Figura 28

33
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 29

Figura 30

El valor “Sink temperature” está relacionado con la forma de


confeccionar la curva de amplitud (“Sink amplitude”). Esta última puede estar
dada en función de valores relativos: se dividen todas las temperaturas de la
tabla por el máximo valor (por ejemplo). En dicho caso en la opción “Sink
temperature” en vez del valor 1, se debería ingresar el valor de la temperatura
máxima. Es decir, que es un valor que multiplica a la curva de amplitud
ingresada.

El mismo proceso se debe realizar para generar todas las interacciones


entre el hormigón y el aire.

Una vez finalizado esto, resta desactivar las interacciones en aquellos


“steps” donde las mismas no deben ser tenidas en cuenta. Para ello es de
utilidad observar la Tabla 7 y luego proceder de la siguiente forma:

34
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Expandir el “árbol” en la interacción que se desea editarà Click derecho en el step


donde no debe existir la interacciónà Deactivate

Figura 31

Conducción Hormigón-Hormigón:

Esta interacción hace referencia a la transferencia de calor que existe


entre capas de hormigón. La formulación es a partir de dos superficies a
distintas temperaturas y entre ellas una interfaz denominada “gap”, a través de
la cual se conduce el calor:

Figura 32

La expresión matemática para este intercambio de calor es la que


sigue:

q = kb ⋅ (θ A − θ B ) Ec. (37)

En donde kb es la propiedad denominada “gap conductance”


(conductividad térmica de la interfaz y θA/B las temperaturas correspondientes a
cada superficie).

Esta interfaz tiene un espesor que tiende a cero, por lo tanto puede
considerarse que la masa contenida en la misma es nula y en consecuencia la
energía interna U puede ser despreciada.

35
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

En la bibliografía consultada y entre los diferentes papers científicos, no


pudo encontrarse un ejemplo en donde se utilizara un valor de kb para modelar
un problema térmico de una estructura de hormigón masivo. Debido a esto, se
utilizó una expresión que está propuesta en el manual del programa, tomando
en cuenta algunas hipótesis que se describirán a continuación:

kliner
kb = Ec. (38)
tliner

En donde kliner es la conductividad del material contenido en la interfaz, y


tliner es el espesor de la misma. Las hipótesis consideradas para calcular kb son:

• El material que forma la interfaz es hormigón, con lo cual resulta


kliner = k.
• El espesor de la interfaz es unitario: tliner = 1 mm.

PARÁMETRO SÍMBOLO VALOR UNIDAD

Conductividad<de<la<interfaz<HºDHº kb 5.79E(07 kcal/mm2.s.ºC

Tabla 11

Al igual que en el caso anterior, lo primero que se debe hacer es generar


en el ensamble las superficies correspondientes a este tipo de interacciones.
En este caso se deben crear dos superficies por cada contacto, una
denominada “master” y otra “slave”. La superficie “master” se considera que
pertenece a la tongada preexistente, mientras que la “slave” corresponde a la
nueva tongada.

La nomenclatura se explica de la siguiente forma:

CONT M12 MASTER

(CONT: Contacto, M12: Partes involucradas en la interfaz,


MASTER/SLAVE: Superficie correspondiente).

Para la generalidad de los problemas, la bibliografía recomienda que las


superficies “slave” sean las que posean mallas más refinadas, o el material
menos rígido si las mallas de ambas superficies son iguales, y deben poseer
elementos de primer orden.

36
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Luego de definir las superficies, se procede a generar las propiedades


que tendrá este tipo de interacciones y que se denominarán:

CONT HH

(CONT: Contacto, HH: Hormigón-Hormigón).

Interaction Properties à Name: CONT HHà Type: Contactà Continue…

Contact Property Options

Mechanicalà Tangential behaviorà Friction formulation: Rough

Mechanicalà Normal behaviorà Pressure-overclosure: “Hard” Contactà Constraint


enforcement method: Defaultà Allow separation after contact (destildar).

Thermalà Thermal conductanceà Definition: Tabularà Use only clearance-dependency


dataà Clearance dependency: Copiar y pegar la Tabla 12à Ok

Conductance******** Clearance******************************
[kcal/mm2.s.ºC] [mm]
5.79E&07 0
0 10

Tabla 12

7.E'07%

6.E'07%

5.E'07%
Gap$Conductance$$
[kcal/mm2.s.ºC]$

4.E'07%

3.E'07%

2.E'07%

1.E'07%

0.E+00%
0% 1% 2% 3% 4% 5% 6% 7% 8% 9% 10% 11%
Clearance$[mm]$

Gráfico 5

Como se observa, las propiedades de la interacción se dividen en


mecánicas y térmicas. Las mecánicas a su vez se dividen en el
comportamiento normal y tangencial. Es preciso aclarar que aunque el modelo
térmico no necesitaría propiedades mecánicas en la definición de la interacción
para calcular las temperaturas y el flujo de calor, es necesario introducirlas para
que el programa interprete la posición de las distintas partes del modelo y las
mantenga durante todo el análisis.

37
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

La opción “Rough” implica que no existe deslizamiento relativo entre


ambas superficies una vez que las mismas entran en contacto. La opción
“Hard” está relacionada con el problema de la penetración de la superficie
“master” en la “slave”. Se debe destildar la opción “allow separation after
contact” para que el programa interprete que una vez que las superficies están
en contacto no se deben separar.

En cuanto a la conductividad de la interfaz, la forma de introducirla es en


función de la separación entre ambas superficies. Cuando éstas están en
contacto, la conductividad alcanza el máximo valor (Kb) y a medida que la
separación entre ambas aumenta, la conductividad tiende a cero.

Por último se deben generar las interacciones cuya nomenclatura es:

CONT M12

(CONT: Contacto, M12: Partes involucradas en el contacto).

Interactions à CONT M12à Step: Step-2 (step en el que “aparece” la


interacción)à Types for selected step: surface-to-surface (Standard)à Continue…

Select the master surface from the dialog: CONT M12 MASTERà Continue…

Choose the slave type: surfaceà Select the slave surface from the dialog: CONT M12
SLAVEà Continue…

Sliding formulation: Small slidingà Discretization methodà Surface to surfaceà Use


supplementary contact points: Selectively

Slave Adjustment: Adjust only to remove overclosure

Surface smoothing: Do not smooth

Clearanceà Initial clearance: Uniform value across slave surface: 0

Contact interaction property: CONT HHà Contact controls: (Default)

Ok

38
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 33

Figura 34

39
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 35

La opción “node-to-surface” en “Discretization method” implica que las


condiciones de contacto se aplican entre los nodos de la superficie “slave” y la
superficie “master”. Por otro lado, la opción “surface-to-surface” considera la
forma de ambas superficies al definir las “constraints” del contacto.

Las opciones “small sliding”/“finite sliding” se refieren a la capacidad de


deslizamiento, rotación relativos entre los nodos de las superficies “master” y
“slave”. Se selecciona la primera de ellas por tener un menor costo
computacional.

La opción “Adjust only to remove overclosure” está relacionada con la


interpenetración de ambas superficies.

40
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Conducción Hormigón-Macizo:

Este tipo de interacción es idéntica a la anterior, por lo tanto se debe


proceder de la misma forma.

Para definir las superficies del contacto pertenecientes al macizo, que


serán de tipo “slave”, primero se deben particionar los “edges” (lados)
correspondientes al mismo. Para ello se debe hacer lo siguiente:

Primero se deben crear los puntos en donde se desea “partir” los lados.
El siguiente esquema es útil para comprender la operación a realizar:

Figura 36

Create Datum Point: Enter coordinates à Coordinates for datum point (x,y,z):
Ingresar las coordenadas correspondientes a cada punto (Tabla 13)

Coordinate) Coordinate) Coordinate))))


Datum)))))))
X))))))))))))))) Y))))))))))))))) Z)))))))))))))))
Point
[mm] [mm] [mm]
1 2000.00 0.00 0.00
2 0.00 5000.00 0.00
3 0.00 10000.00 0.00

Tabla 13

41
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 37

Luego, para realizar las particiones se procede de la siguiente forma:

Toolsà Partition…à Type: Edgeà Method: Select mid point / datum pointà Select an
edge to partition (Seleccionar el lado que se desea particionar)à Select a point on the
edge (Seleccionar el datum point correspondiente).

Figura 38

42
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Luego de realizar las particiones, se procede a generar las superficies


como en la interacción hormigón-hormigón. En este caso la nomenclatura será:

CONT M1T MASTER

(CONT: Contacto, M1T: Superficies de partes involucradas en la interfaz,


MASTER/SLAVE: Si se trata de la superficie de hormigón o del macizo
respectivamente).

Las superficies correspondientes a una interacción de tipo “contact pair”,


no pueden intersectarse, es decir no pueden compartir nodos. Con lo cual, para
la superficie correspondiente a la “instance” M1, se deben seleccionar ambos
“edges”.

Las propiedades de la interacción se denominarán:

CONT HM

(CONT: Contacto, HM: Hormigón-Macizo)

La hipótesis para establecer el valor de la conductividad térmica en la


interfaz es la siguiente:

1
⋅ (k + k )
kliner 2 H M
kb = = Ec. (39)
tliner 1

Es decir que se toma una conductividad promedio entre el hormigón y el


macizo rocoso las cuales fueron presentadas en la Tabla 4 (apartado 4.4.1.2
Material).

PARÁMETRO SÍMBOLO VALOR UNIDAD

Conductividad0de0la0interfaz0Hº8M kb 5.48E(07 kcal/mm2.s.ºC

Tabla 14

Por último, se generan las interacciones como en el caso anterior, cada


una en el step que corresponda con la siguiente nomenclatura:

CONT M1T1

(CONT: Contacto, M1T1: Superficies de partes involucradas en la


interfaz).

43
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

NOTA: No se tiene en cuenta en este caso la radiación solar que puede


incidir sobre la estructura porque son datos a los que es complicado acceder y
no modifican sustancialmente los resultados. En ninguno de los trabajos
publicados de problemas térmicos relacionados al calor de hidratación que se
han consultado se encontraron referencias acerca de la consideración de la
radiación solar.

4.4.1.7 Loads

La carga térmica que se debe introducir en el programa, es el flujo de


calor volumétrico correspondiente al calor de hidratación del cemento.

En general, se parte de la Curva de Elevación Adiabática de


temperatura, que se obtiene de un ensayo en el cual se mide cada
determinado intervalo de tiempo la temperatura que va ganando el material en
condiciones adiabáticas (sin intercambio de calor con el medio). Estas
condiciones son representativas de las zonas “nucleares” de las estructuras,
pero se acepta su utilización para todo el conjunto por su simplicidad y porque
no implica la generación de errores importantes en los resultados. El ACI
propone un ensayo en el artículo CRD-C 38-73 “Method of test for temperatura
rise in concrete”.

Figura 39: Esquema de un calorímetro adiabático

44
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Gráfico 6: Curva de elevación adiabática de temperatura resultado de un ensayo.

Si en esta curva, las temperaturas son por kg de cemento, para obtener


la curva correspondiente al hormigón, se deben multiplicar las mismas por el
C.U.C. (Contenido Unitario de Cemento).

Luego, a partir de las propiedades térmicas del material, se transforma la


Curva de Elevación Adiabática de temperatura, en Curva de Flujo de Calor
Volumétrico, utilizando la siguiente expresión:

dT dH " 1 %
= ⋅$ ' Ec. (40)
dt dt #ρ ⋅c&

- dT/dt es la variación de la temperatura en función del tiempo.

- dH/dt es la variación del calor de hidratación por unidad de volumen en


función del tiempo, es decir el flujo de calor volumétrico.

Como el dominio en el que está definida la curva de temperaturas es


discreto, la expresión anterior resulta:

ΔH ΔT
= ⋅ρ ⋅c Ec. (41)
Δt Δt

45
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Gráfico 7: Curvas de flujo de calor volumétrico para hormigones con distinto tipo de cemento.

Cuanto menor sea el intervalo de tiempo entre mediciones, y mayor sea


el número de puntos de los que se disponga en la curva, más precisa será la
simulación numérica.

En la siguiente imagen se muestra la tabla que se debe confeccionar


para definir la curva de calor que se introduce en el programa:

Figura 40

Otra forma de trabajar es a partir de una curva de ajuste del calor de


hidratación por kg de cemento dada por la siguiente expresión:

# τ &α
−% (
$ t+12 '
H (t) = H 0 ⋅ e Ec. (42)

Los parámetros H0, τ, α son otorgados por el fabricante del cemento a


partir de ensayos experimentales. En este caso se tiene:

46
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

PARÁMETRO SÍMBOLO VALOR UNIDAD

CAPACIDAD CALÓRICA DEL CEMENTO Q0 143.09 Kcal/kg


COEFICIENTE DE AJUSTE ζ 116.00 hs
COEFICIENTE DE AJUSTE α 0.40 +

Tabla 15

El procedimiento para determinar las curvas es inverso al anterior. Los


intervalos de tiempo utilizados pueden ser de 1/8 día para la edad temprana del
hormigón (7-14-28 días) y luego pasar a ser de 1 día cuando el flujo de calor se
estabiliza y tiende a un valor constante.

Figura 41

Figura 42

47
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Curva%de%Flujo%de%Calor%Volumétrico%
6.E'10%

5.E'10%

DH/Dt%[kcal/mm3.s]%
4.E'10%

3.E'10%

2.E'10%

1.E'10%

0.E+00%

1.E+07%
0.E+00%

1.E+06%

2.E+06%

3.E+06%

4.E+06%

5.E+06%

6.E+06%

7.E+06%

8.E+06%

9.E+06%
t%[s]%

Gráfico 8

Curva"de"Elevación"Adiabá8ca"de"
Temperatura"
80"

70"

60"
T"[ºC]"

50"

40"

30"

20"

10"

0"
1.E+07"
0.E+00"

1.E+06"

2.E+06"

3.E+06"

4.E+06"

5.E+06"

6.E+06"

7.E+06"

8.E+06"

9.E+06"

t"[s]"

Gráfico 9

El siguiente paso, consiste en confeccionar tres curvas de flujo de calor,


una para cada tongada de hormigón que posee el problema, ya que el
desarrollo de la misma se da para distintos tiempos. Estas tres curvas entonces
deben estar en función del tiempo total del análisis y las tres deben finalizar en
el día 100 (8640000 segundos).

48
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Curvas%de%Calor%para%las%3%tongadas%
6.E'10%

DH/Dt%[kcal/mm3.s]%
5.E'10%

4.E'10%
H1%
3.E'10% H2%

H3%
2.E'10%

1.E'10%

0.E+00%

2.E+05%

4.E+05%

6.E+05%

8.E+05%
0.E+00%

1.E+06%

1.E+06%

1.E+06%
t%[s]%

Gráfico 10

Luego de generar en una tabla las curvas de calor correspondientes a


cada tongada se deben crear las curvas de amplitud en el programa, de la
misma forma que se hizo para la temperatura ambiente. En este caso se
propone la siguiente nomenclatura:

AMP H1, AMP H2, AMP H3

Lo siguiente es generar la carga, para ello se deben seguir los siguientes


pasos:

Loads à Name: H1à Step: Step-1à Category: Thermalà Types for Selected Step:
Body heat fluxà Continue…

Select bodies for the load: M1 (Seleccionar en el ensamble la instance


correspondiente)à Done

Distribution: Uniformà Magnitude: 1à Amplitude: AMP H1 (Seleccionar la curva de


amplitud correspondiente)à Ok

49
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 43

Figura 44

Figura 45

50
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Lo mismo se debe repetir para las cargas H2 y H3.

4.4.1.8 BCs -Boundary Conditions-

Temperatura “estable” del macizo de fundación:

Al definir la geometría, se indicó que el macizo de fundación se


modelaba a cada lado de la estructura una distancia 0.5 a 2L. Esta es una
dimensión razonable a partir de la cual no existe transferencia de calor, es
decir que no existe influencia del calor cedido por el hormigón.

A partir de esta distancia, algunos autores toman la hipótesis de que la


temperatura del terreno es constante e igual a la temperatura media anual del
ambiente que en este caso resulta:
Tmedia======
Mes
[ºC]
Ene 28.00
Feb 29.00
Mar 29.00
Abr 28.50
May 28.50
Jun 27.00
Jul 26.00
Ago 26.50
Sep 26.00
Oct 26.50
Nov 26.00
Dic 25.50
Tmedia=
27.21
anual=[ºC]

Tabla 16

Para ingresar esta condición de borde se debe proceder de la siguiente


manera:

BCs à Name: TEMà Step: Step-1à Category: Otherà Types for Selected Step:
Temperatureà Continue…

Selec regions for the boundary condition (Seleccionar todas las superficies límite del
macizo manteniendo apretado alt)

Distribution: Uniformà Magnitude: 27.21à Amplitude: (Instantaneous)

51
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 46

Figura 47

52
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 48

Temperatura del agua luego del llenado de la esclusa:

Como se describió en el planteo del problema, se considera que a partir


del día 60 la esclusa se encuentra sumergida. Con lo cual las superficies
laterales derechas del muro tendrán la misma temperatura que el agua.

TSM;;;;;;;;;
Mes
[ºC]
Ene 27.50
Feb 27.00
Mar 27.50
Abr 28.00
May 28.50
Jun 28.50
Jul 28.00
Ago 28.00
Sep 28.00
Oct 28.00
Nov 28.00
Dic 27.50
TSM;media;
27.88
[ºC]

Figura 49

Como se puede observar, la temperatura del agua durante el período de


tiempo del análisis no varía significativamente, con lo cual se tomará el valor
TSM = 28 ºC.

La manera de proceder es la misma que en el caso anterior, sólo que en


este punto se deben seleccionar las superficies del muro en contacto con el
agua y en el step correspondiente (Step-4).

53
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4.4.1.9 Predefined Fields

Temperatura de colocación del hormigón:

Como se indicó en un principio, las tres tongadas de hormigón tienen la


misma temperatura de colocación: Tc=20ºC.

Predefined Fields à Name: TCà Step: Initialà Category: Otherà Types for Selected
Step: Temperatureà Continue…

Select regions for the field or press Done to use calculated temperaturas (Seleccionar las
instances M1, M2, M3)à Done

Distribution: Direct specificationà Section variation: Constant through


regionà Magnitude: 20ºC

Figura 50

54
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 51

Temperatura inicial del macizo de fundación:

Lo mismo se debe hacer para definir la temperatura inicial en el terreno


que rodea la esclusa, que en este caso se considera igual a la temperatura
media anual: Ti=27.21ºC.

4.4.1.10 Mesh

Para definir la malla se deben tener en cuenta tres aspectos:

La forma de los elementos:

Pueden ser cuadriláteros, triangulares, de mallado regular, libre, etc.

Module: Mesh

Assign Mesh Controls à Select the regions to be assigned mesh controls


(Individually)à Done

Element Shape: Quadà Technique: Structuredà Ok

55
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 52

Figura 53

Lo mismo se debe hacer con la “instance” correspondiente al macizo,


pero en este caso será: Element Shape=Tri, Technique=Free.

El tamaño de los elementos:

Para ello debemos establecer en los “edges” de cada “instance” los


puntos donde deben existir nodos, luego el programa malla automáticamente la
“instance” a partir de la información ingresada.

56
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Module: Mesh

Seed Edges à Select the regions to be assigned local sedes (Individually).


(Seleccionar los 4 lados de las instances M1, M2, M3)à Done

Basicà Method: By sizeà Bias: Noneà Sizing Controls. Approximate element size:
250à Curvature control (Destildar)

Constraintsà Seed constraints: Do not allow the number of elements to change

Ok

Done

Figura 54

57
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 55

Para el caso del macizo se debe proceder de una manera distinta para
densificar la malla en la zona próxima al muro. Para ello, seleccionar los
“edges” que están en contacto con el muro y el inferior a la derecha del mismo,
y asignarles el mismo tamaño de elementos como se hizo en el caso anterior
(250 mm).

Luego seleccionar los “edges” opuestos, es decir donde la temperatura


es estable y asignarles un tamaño de elementos de 750 mm.

Por último seleccionar los “edges” restantes y asignarles:

Method: By sizeà Bias: Singleà Sizing Controls. Minimum size: 250à Maximum size:
750.

Flip bias: Invierte el sentido de densificación.

Constraintà Seed Constraintsà Allow the number of elements to increase or decrease.

58
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 56

Lo que resta es generar la malla:

Mesh Part Instance à Select the part instances to be Meshed (Seleccionar


todas)à Done

Figura 57

59
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

El tipo de elementos:

Hace referencia al tipo de formulación de los mismos: si se trata de un


análisis mecánico, de transferencia de calor, térmico-eléctrico, etc. ,si la
interpolación es lineal, de segundo orden, etc.

Module: Mesh

Assign Element Type à Select the regions to be assigned element types


(Seleccionar todas las instances)à Done

Element Library: Standardà Family: Heat Transferà Geometric order: Linear

DC2D4: a 4-node linear heat transfer quadrilateral

DC2D3: a 3-node linear heat transfer triangle

Figura 58

60
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 59

4.4.1.11 Element Birth and Death Technique

Para modelar las distintas etapas constructivas de este problema se


debe editar el “keywords”. En esta sección es donde figura el “input” del
modelo, es decir todos los comandos necesarios para que el programa
resuelva el problema.

El primer paso es crear 3 sets de elementos (ESM1, ESM2, ESM3),


cada uno correspondiente a cada tongada de hormigón (M1, M2, M3):

Assembly à Sets à Name: ESM1à Type: Elementà Continue…à Select the


elements for the set (Seleccionar todos los elementos de la instance
correspondiente)à Done

Figura 60

61
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 61

Luego, se procede a editar el “keywords”:

Toolbarà Modelà Edit keywordsà MOD TER ESC 001

En la ventana, “scrollear” hacia abajo hasta hallar el input correspondiente al primer


step. Luego de la definición del mismo (nombre, duración, intervalos máximos y
mínimos), escribir lo siguiente:

*MODEL CHANGE, REMOVE, TYPE=ELEMENT

ESM2

*MODEL CHANGE, REMOVE, TYPE=ELEMENT

ESM3

Esto implica que en el primer step, se ingresa un comando para que el


programa no tenga en cuenta en la simulación todos los elementos de las
tongadas M2 y M3.

Lo siguiente es modificar el input para que estos elementos sí sean


tenidos en cuenta en los steps correspondientes:

En la ventana, “scrollear” hacia abajo hasta hallar el input correspondiente al segundo


step. Luego de la definición del mismo (nombre, duración, intervalos máximos y
mínimos), escribir lo siguiente:

62
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

*MODEL CHANGE, ADD, TYPE=ELEMENT

ESM2

Repetir la operación en el Step-3 para el Element Set ESM3

Figura 62

Figura 63

63
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4.4.1.12 Ask for nodal temperature field (.fil file)

Como se explicó en el apartado 3.2.2, este problema se resuelve a partir


de un procedimiento SC. Con lo cual es necesario generar un segundo modelo
mecánico en donde se introduce como input las temperaturas obtenidas como
output en este modelo térmico.

En este apartado se describe el comando necesario para que el


programa escriba un archivo especial con todas las temperaturas nodales para
cada paso de cálculo. Este archivo (de extensión .fil), es el que luego se
introduce en el modelo mecánico.

Es importante destacar, que el archivo de extensión .prt de ambos


modelos (térmico y mecánico), debe ser el mismo. Esto implica que las partes
en ambos modelos sean las mismas y que la posición y nomenclatura de los
nodos sea coincidente.

Lo primero es crear 3 sets de nodos para cada parte (NSM1, NSM2,


NSM3). Para ello se procede de la misma forma que se explicó para los sets de
elementos, sólo que en este caso se debe tildar la opción “Node” en vez de
“Element” y se deben seleccionar los nodos de la “instance” correspondiente en
cada caso.

Luego se debe editar el “keywords” de la siguiente forma:

Toolbarà Modelà Edit keywordsà MOD TER ESC 001

En la ventana, “scrollear” hacia abajo hasta hallar el input correspondiente al primer


step. Antes de finalizar el mismo (*End Step), escribir lo siguiente:

*OUTPUT, FIELD

*NODE FILE, NSET=NSM1

NT

Para el caso del Step-2, se debe ingresar el comando para escribir el archivo de las
temperaturas correspondientes a M1 y M2, con lo cual se debe escribir:

*OUTPUT, FIELD

*NODE FILE, NSET=NSM1

NT

*OUTPUT, FIELD

64
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

*NODE FILE, NSET=NSM2

NT

Lo mismo se debe repetir para los restantes steps.

Figura 64

La temperatura del macizo no es necesaria para el modelo mecánico,


con lo cual no se tiene en cuenta para la escritura del archivo “.fil”.

4.4.1.13 Job

Para finalizar este modelo térmico se debe crear el “Job”, en este punto
el programa permite establecer las características del procesamiento de datos,
el porcentaje de memoria física que se desea utilizar, la cantidad de
procesadores destinados a la resolución del problema, etc.

Jobs à Name: JOBTERESC001à Source: Modelà MODTERESC001à Continue…à Ok

65
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 65

Figura 66

Para correr el análisis, se debe hacer click derecho sobre el “Job”


correspondiente y seleccionar la opción “Submit”.

66
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4.4.2 Output

4.4.2.1 Generalidades

Para acelerar el postproceso, es conveniente generar sets de nodos en


las zonas críticas de la estructura donde se desea obtener información, antes
de comenzar el “job”. Estas zonas críticas se corresponden con los fenómenos
de gradiente de masa y de superficie como se mencionaba en un principio.

En general se toman ejes ortogonales para las distintas partes de la


estructura, de forma tal de analizar los gradientes térmicos para diferentes
instantes de tiempo.

A continuación se presentan los contornos de temperaturas para


distintos tiempos en el modelo:

Figura 67

67
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 68

Figura 69

Para guardar una imagen del “viewport” se procede de la siguiente


forma:

Fileà Printà Print: Current viewportà Print viewport decorations (if visible)à Destination:
fileà Format: PNGà Ok

68
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 70

Para crear una animación de la evolución de las temperaturas, hacer lo


siguiente:

Animate: Time History

Animateà Save As…à Format: .AVIà Capture: Current viewportà Capture viewport
decorations (if visible)à Frame Rate: 24à Ok

Figura 71

4.4.2.2 Historia de temperaturas nodales

Puede resultar de interés evaluar el pico de temperaturas para distintas


partes de la estructura. Para ello se debe hacer lo siguiente:

Create X,Y data

69
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

ODB field outputà Continue…

Variablesà Position: Unique Nodalà NT11: Nodal temperature

Elements/Nodesà Method: pick from viewportà Edit selection: Seleccionar los nodos en
donde se quieran obtener datos manteniendo apretada la tecla altà Save

Figura 72

Figura 73

Luego se debe crear un archivo .txt con los datos creados en el paso
anterior. A partir de este archivo se pueden copiar los datos a Excel.

Reportà XY…

X,Y Data: Seleccionar los datos correspondientes

Setupà Name: Cambiar el nombre por defectoà Ok

70
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 74

Luego, en el archivo OUTPUT MOD ESC.xlsx se puede trabajar con los


datos para crear mejores gráficos que los que permite ABAQUS.

Historia"de"temperaturas"3er"tongada"
50"

45"

40"

35"

30"
T"[ºC]"

25"

20"

15" Nodo"int."
10" Nodo"ext."

5"

0"
0" 10" 20" 30" 40" 50" 60" 70" 80" 90" 100" 110"
t"[días]"

Gráfico 11

Se puede observar que para la tongada nro. 3, la temperatura se


mantiene constante con un valor de 20ºC (temperatura de colocación) hasta el
día 6 que se corresponde con el tercer “step”. Para el nodo en correspondencia
con la superficie, en el día 60 se produce un salto en la temperatura debido al
llenado de la esclusa con agua (cambio en la condición de borde).

Ambas zonas de la estructura presentan temperaturas que tienden al


equilibrio térmico a unos 28ºC. Mientras que en la zona central, la temperatura

71
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

alcanza los 47ºC, en la zona de la superficie la temperatura máxima es de


32ºC.

Para la generalidad de los problemas relacionados con el calor de


hidratación del hormigón en condiciones normales, los resultados de la
simulación no debieran superar la temperatura máxima de la Curva de
Elevación Adiabática. Esto implica un primer punto de control para establecer si
se ha cometido algún error al cargar el input en el modelo, o en las unidades
utilizadas, etc.

4.4.2.3 Gradientes térmicos

También puede ser de utilidad, extraer las temperaturas a lo largo de


todo el espesor del muro para distintos instantes de tiempo. Para ello se debe
generar un “path”:

Paths

Name: PATHM1X1à Type: Node Listà Continue…

Part instance: M1à Node Labels: 91:99:1à Ok

En la opción “Node Labels”, el 100 indica el “start node”, el 108 el “end


node”, y el 1 el intervalo de nodos del “path”. Cuando sea difícil establecer el
intervalo o no se puedan leer bien las etiquetas de los nodos, puede utilizarse
la opción “Add After” en la cual se deben seleccionar los nodos del camino en
el “viewport” de manera individual.

Figura 75

72
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 76

Lo siguiente es crear los datos correspondientes. Para ello se procede


de la siguiente forma:

X,Y Data

Pathà Continue…

Path: PATHM1X1à Model Shape: Undeformedà X Values: True distanceà Y Values.


Step/Frame: Step-1à Frame: Seleccionar el tstep deseado (1 día, 2 días, etc)à Ok

Save As…à Name: PATHM1X1DIA1

Figura 77

73
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 78

Figura 79

74
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 80

Cuando se debe trabajar con varios días correspondientes a distintos


“steps”, es conveniente trabajar con una tabla como la que sigue:

t%%%%%%%%%%%%%%%%%% t*%%%%%%%%%%%%%%%%%%
Step%Time%%%%%%%%%%%%%%%%%%
Día STEP FRAME%%%%%%%%%%%%%%%%%
[s] [s] [s]

1 8.64E+04 Step-1 8.64E+04 9.67E+04 30


3 2.59E+05 Step-2 0.00E+00 0.00E+00 0
7 6.05E+05 Step-3 8.64E+04 9.22E+04 31
14 1.21E+06 Step-3 6.91E+05 6.77E+05 38
28 2.42E+06 Step-3 1.90E+06 1.88E+06 52
100 8.64E+06 Step-4 3.46E+06 3.46E+06 68

Tabla 17

En la primer columna figuran los días en los que se desea extraer datos
de temperaturas. En la columna a la derecha figura el tiempo total
correspondiente en segundos. Luego el “step” al que pertenecen dichos
instantes de tiempo. En la cuarta columna se presenta el tiempo que
correspondería al “step” (restando al tiempo total la duración del step anterior).
Luego se debe buscar en el programa los “step time” que más se aproximen a
los valores anteriores.

Por último se debe crear el archivo .txt de la misma forma que para el
apartado anterior.

75
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Gradientes"Horizontales"1er"Tongada"
55"

50"

45"

40"

35"

T"[ºC]"
30"

25"

20" Día"1" Día"3"

15"
Día"7" Día"14"
10"

5" Día"28" Día"100"

0" 0.1"

0.5"

0.7"

0.8"

0.9"

1.1"

1.5"

1.7"

1.8"

1.9"
0.0"

0.2"

0.3"

0.4"

0.6"

1.0"

1.2"

1.3"

1.4"

1.6"

2.0"
X"[m]"

Gráfico 12

En este gráfico se muestran los gradientes para una recta horizontal de


nodos ubicada en el centro de la primer tongada de hormigón. Se puede
observar que los gradientes aumentan siendo el máximo mostrado en el gráfico
el día 7, y luego comienzan a disminuir hasta alcanzar el equilibrio térmico (día
100).

Gradientes"Ver3cales"1er"Tongada"
5.0#

4.5#

4.0#

3.5# Día#1#

Día#3#
3.0#
Y"[m]"

Día#7#
2.5#
Día#14#
2.0#
Día#28#
1.5#
Día#100#
1.0#

0.5#

0.0#
0#

5#

10#

15#

20#

25#

30#

35#

40#

45#

50#

55#

T"[ºC]"

Gráfico 13

En el gráfico anterior se muestran las temperaturas en distintos instantes


de tiempo para una recta vertical de nodos que pasa por el centro de la primer
tongada de hormigón. Se observa que los gradientes aumentan (día 1 – día 3)
y luego dejan de ser simétricos respecto del centro de la tongada de hormigón
debido a que aparece la segunda tongada y comienza a producirse la
transferencia de calor entre ambas. Por último los gradientes van disminuyendo
hasta alcanzar el equilibrio térmico (día 100).

76
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4.5 Modelo Mecánico


En este apartado se discutirá la manera de generar el segundo modelo
secuencialmente acoplado, con el objetivo de obtener las tensiones y
deformaciones térmicas.

Lo primero que se debe hacer es copiar el modelo térmico y asignarle


otro nombre. En este caso se propone:

MOD MEC ESC 001

(MOD: Modelo, MEC: Mecánico, ESC: Esclusa, 001: Versión).

Es importante hacer una aclaración antes de comenzar con el desarrollo


de esta parte del tutorial. Para simular la secuencia constructiva en el modelo
mecánico se puede proceder de dos formas:

Una de ellas es trabajando con “Contact-Pairs” de la misma forma que


en el modelo térmico anteriormente explicado, incorporando la denominada
“Elements Birth and Death Technique”. Se ha observado que el campo de
tensiones y deformaciones no resulta continuo en las interfaz, con lo cual se
descartó su utilización.

La otra opción es utilizar “Tie-Constraints” para establecer la continuidad


de tensiones y deformaciones en la interfaz. En consecuencia se debe ingresar
una curva de amplitud para el módulo de elasticidad. Esta curva debe tener un
valor muy bajo (cercano al valor nulo sin que genere problemas de
convergencia) en aquellos instantes de tiempo en los cuales la tongada no
debe “existir” en la simulación. Esto permite que las tensiones en esa misma
tongada sean nulas en dichos instantes de tiempo.

La utilización de la primera de las técnicas mencionadas es similar a lo


visto en el modelo térmico. En este tutorial se explicará la segunda opción,
teniendo el lector la total libertad de continuar investigando la forma más
conveniente y precisa de simular los contactos entre las partes y la simulación
de las etapas constructivas de la estructura.

Continuando con el tutorial, se deben eliminar todos los aspectos


relacionados al modelo térmico en un orden tal que no impliquen problemas de
escritura en el input:

77
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

1) Predefined Fieldsà Borrar la temperatura inicial del macizo. (No se debe eliminar
la temperatura inicial de las tongadas de hormigón porque el archivo “.fil”
escribe las temperaturas a partir del Step-1, y no tiene en cuenta las del step
Initial (en el cual está definida la temperatura de colocación).

2) BCsà Borrar la temperatura estable del macizo y la temperatura del agua.

3) Loadsà Eliminar todas las cargas térmicas (H1, H2, H3).

4) Amplitudesà Eliminar todas las curvas de amplitud (las tres correspondientes al


calor de hidratación y la temperatura ambiente).

5) Interactionsà Debido a que no se trabajará con “Contact-Pairs”, eliminar todas


las interacciones definidas.

6) Interaction Propertiesà Eliminar todas las propiedades de contacto definidas en


el modelo térmico.

7) Stepsà Eliminar todos los steps, ya que no deben ser definidos como de tipo
“Heat Transfer”.

8) Surfacesà Eliminar las correspondientes a las interacciones de convección con


el aire. No borrar las correspondientes a los “Contact-Pairs” ya que serán de
utilidad para definir los “Tie-Constraints”.

9) Instancesà Eliminar la “Part Instance” correspondiente al macizo.

10) Eliminar la designación de secciones en las partes: Partsà M1à Section


Assignmentsà Click derechoà Delete…

11) Sectionsà Eliminar todas las secciones.

12) Materialsà Borrar todos los materiales definidos.

13) Partsà Eliminar la parte correspondiente al macizo.

14) Borrar todo el texto editado en el “Keywords” (*MODEL CHANGE y *OUTPUT,


FIELD).

4.5.1 Input

4.5.1.1 Parts

Como se explicó anteriormente, la geometría en este caso consiste


únicamente en las partes correspondientes al muro de hormigón, ya que se
considera que no existe interacción mecánica (fricción por ejemplo) entre éste y
el macizo.

78
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4.5.1.2 Materials

En la siguiente tabla se presentan las propiedades mecánicas del


hormigón a la edad de 28 días utilizadas para resolver este tutorial. Más abajo
se describe brevemente cada una en particular.

PROPIEDAD SÍMBOLO VALOR UNIDADES

RESISTENCIA MEDIA A LA COMPRESIÓN A 28 DÍAS f'cm28 37.00 MPa


RESISTENCIA MEDIA A LA TRACCIÓN A 28 DÍAS f'tm28 3.60 MPa
MÓDULO DE ELASTICIDAD ESTÁTICO A 28 DÍAS Ec28 33100.00 MPa
COEFICIENTE DE EXPANSIÓN TÉRMICA α 1.05E105 1/ºC
COEFICIENTE DE POISSON u 0.25 1
PESO DE LA UNIDAD DE VOLUMEN P.U.V. 2.51E105 N/mm3

Tabla 18

Resistencia a Compresión (f’c):

Esta es una propiedad fundamental del hormigón que se utiliza para el


proyecto de las estructuras (durabilidad, vida útil, etc.), el dimensionamiento de
las mismas frente a los esfuerzos característicos (resistencia característica a 28
días), etc.

Depende de la relación agua/cemento de la mezcla, del tipo de cemento


utilizado, etc.

Se realizan ensayos destructivos para su medición en laboratorios


(ASTM C39 / C39M - 14 “Standard Test Method for Compressive Strength of
Cylindrical Concrete Specimens”).

Figura 81

Durante el proceso de fraguado, a medida que se forman los distintos


productos de hidratación, la pasta de cemento y el hormigón en forma global

79
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

ganan resistencia. La siguiente es una curva de ajuste extraída del Código


Modelo CEB-90, para describir el comportamiento de la resistencia media a la
compresión con el tiempo para condiciones standard:
) # &,
+ % 28 (.
+α ⋅%1− t (.
+* %$ t1 ('.-
f 'cm = f 'cm28 ⋅ e Ec. (43)

- α es un coeficiente que depende del tipo de cemento utilizado en la


mezcla.

- t1 es un parámetro de ajuste, (en este caso se toma igual a 1 día).

Tipo%de%cemento α
Endurecimiento%rápido%y%alta%resistencia 0.20
Endurecimiento%normal/rápido 0.25
Endurecimiento%lento 0.38

Tabla 19

Evolución%de%la%Resistencia%a%Compresión%
45"

40"

35"

30"
f'cm%[MPa]%

25"

20"

15"

10"

5"

0"
0" 10" 20" 30" 40" 50" 60" 70" 80" 90" 100"
t%[días]%

Gráfico 14

Resistencia a Tracción (f’t):

La resistencia a tracción se relaciona directamente con el origen y


propagación de las fisuras de origen térmico como se explicó al principio del
tutorial.

Como es sabido, alcanza valores mucho menores que la resistencia a


compresión (aproximadamente 10%). Tiene una evolución similar a esta última
y existen distintos ensayos para su determinación entre los que se pueden
nombrar:

80
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

• Ensayo por tracción directa: No se utiliza por la dificultad práctica


que existe en imprimirle una carga axial de tracción a una
probeta.
• Ensayo de tracción por flexión en 3 puntos (ASTM C 78): Se
ensaya una viga de hormigón simplemente apoyada y cargada en
2 puntos ubicados a los tercios. Se mide la resistencia a tracción
utilizando la expresión de Navier:

M P⋅L
fr = = Ec. (44)
W b ⋅ d2

Figura 82

• Ensayo de tracción por compresión diametral (ASTM C 496): Se


ensaya una probeta cilíndrica cargada en forma
aproximadamente lineal sobre una generatriz . La resistencia a
tracción se calcula luego con la siguiente expresión:

2⋅P
f 't = Ec. (45)
π ⋅l ⋅d

Figura 83

Distintos reglamentos ofrecen expresiones para estimar la resistencia a


tracción del hormigón:

81
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

1
ACI 209R fr = gr ⋅ ( wc ⋅ f 'c ) 2
Ec. (46)

fr = Módulo de rotura [MPa].

wc = Peso unitario [kg/m3].

f’c = Resistencia a compresión [MPa].

gr = 0.012 a 0.021
1
ACI 209R f 't = gt ⋅ ( wc ⋅ f 'c ) 2
Ec. (47)

fr = Resistencia a tracción directa [MPa].

wc = Peso unitario [kg/m3].

f’c = Resistencia a compresión [MPa].

gt = 0.0069
2
CEB 90 f 'tm = 0.30 ⋅ ( f 'ck ) 3 Ec. (48)

f’tm = Resistencia media a tracción [MPa].

f’ck = Resistencia característica a compresión [MPa]

La siguiente expresión (Mehta y Monteiro) permite calcular la resistencia


a tracción para cualquier edad del hormigón, en función de la resistencia a la
compresión para dicha edad:
2
" f ' (t) %
f 'tm (t) = 1, 7 ⋅ $ 3 cm ' Ec. (49)
# 12 &

Evolución%de%la%Resistencia%a%Tracción%
4.5#

4.0#

3.5#

3.0#
f'cm%[MPa]%

2.5#

2.0#

1.5#

1.0#

0.5#

0.0#
0# 10# 20# 30# 40# 50# 60# 70# 80# 90# 100#
t%[días]%

Gráfico 15

82
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Módulo de Elasticidad (E):

El módulo de elasticidad relaciona las tensiones con las deformaciones,


compone junto con la inercia de las secciones resistentes la rigidez de las
estructuras. Depende de la edad del hormigón, del contenido unitario y tipo de
agregados, de la relación agua/cemento, etc.

Como es sabido, la relación entre tensiones y deformaciones en el


hormigón es lineal únicamente para bajos valores de la carga. Con lo cual es
complicado definir un único valor para el módulo de elasticidad, y establecer la
forma de medirlo.

!
Gráfico 16

Según el tipo de cargas (instantáneas, dinámicas, permanentes, etc.)


con el cual se esté trabajando, se deberían utilizar algunos de los siguientes
módulos de elasticidad definidos para el hormigón:

• Módulo de Elasticidad Estático


• Módulo de Elasticidad Secante
• Módulo de Elasticidad Cuerda
• Módulo de Elasticidad Tangente en un Punto
• Módulo de Elasticidad Tangente Inicial

83
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

!
Gráfico 17

El artículo “Test for Static Modulus of Elasticity and Poisson’s Ratio of


Concrete in Compression (ASTM C 469)”, proporciona un método para el
cálculo del módulo de elasticidad cuerda:

Se utiliza una probeta cilíndrica de 15 cm x 30 cm y se la carga en


compresión midiendo las deformaciones longitudinales y transversales. El
módulo de elasticidad cuerda se calcula de la siguiente forma:

S2 − S1
Ec = Ec. (50)
ε 2 − 0.00005

Donde S2 = tensión correspondiente al 40% de la resistencia última a


compresión.

S1 = tensión correspondiente a una deformación unitaria de 50E-06.

ε2 = deformación unitaria correspondiente a S2.

Distintos reglamentos proponen fórmulas para la estimación de Ec sin


necesidad de realizar ensayos:

ACI 318 Ec = w1.5 ⋅ 33⋅ f 'c Ec. (51)

E = Módulo de elasticidad estático [psi].

w = Peso de la unidad de volumen [lb/ft3].

f’c = Resistencia característica a compresión [psi].

84
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

f 'cm
CEB 90 Ec = 2.15 ×10 4 ⋅ 3 Ec. (52)
10

E = Módulo de elasticidad estático [MPa].

f’cm = Resistencia media a compresión [MPa].

CIRSOC 201/05 Ec = 4700 ⋅ f 'c Ec. (53)

E = Módulo de elasticidad estático [MPa].

f’c = Resistencia característica a compresión [MPa].

Para establecer la variación del módulo de elasticidad con la edad, se


puede recurrir a la siguiente expresión dada por el Código Modelo CEB 90:
0,5
/ ) #
+ % 28 ( .
&0,5 , 3
1 0,26+1−% ( .1
+ %$ t t (' . 1
1 * 1 -
Ec (t) = Ec28 ⋅ 0e 4 Ec. (54)
1 1
12 15

También para estructuras de gran envergadura como presas, o pilas de


puente, etc. puede ser necesario considerar el efecto del creep. En breves
palabras, el creep es el aumento de la deformación en el hormigón bajo una
carga constante. Depende de del volumen de agregados, de la humedad
ambiente, de la carga impuesta, del momento en que se produce dicha carga,
entre otros factores y está relacionado con la migración del agua químicamente
adsorbida en la pasta de cemento.

Gráfico 18

85
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Para tener en cuenta el efecto del creep, se introduce el denominado


Módulo de Elasticidad Sostenido, cuyo valor es menor que el Módulo de
Elasticidad Estático:

1
Ecsost (t) = ⋅ E (t) Ec. (55)
! log(t +1) $ c
#"0.94 + 4.3 &%

En el siguiente gráfico figura la variación con el tiempo de los módulos


de elasticidad anteriormente definidos:

Evolución"del"Módulo"de"Elas5cidad"
40000"

35000"

30000"
E"[MPa]"

25000"

20000"

15000" Ec"

10000" Ec"sost"

5000"

0"
0" 10" 20" 30" 40" 50" 60" 70" 80" 90" 100"
t"[días]"

Gráfico 19

Coeficiente de Poisson (µ):

Es el parámetro elástico que relaciona las deformaciones longitudinales


de una pieza, con las transversales.

Según el ASTM C 469, el coeficiente de Poisson se calcula de la


siguiente forma:

εt2 − εt1
µ= Ec. (56)
εt2 − 0.00005

εt2 = Deformación transversal producida por S2.

εt1 = Deformación transversal producida por S1.

86
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

!
Figura 84

!
Gráfico 20

Coeficiente de Expansión Térmica (α):

Este coeficiente relaciona las temperaturas en el material con las


deformaciones. Depende en gran medida del coeficiente de expansión térmica
del agregado grueso como figura en el siguiente esquema:

87
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Gráfico 21

A continuación se explicará la forma de ingresar las propiedades


mecánicas en el modelo computacional. Si se consideran constantes para
cualquier edad del hormigón, la forma de introducirlas en el input es similar a lo
visto en el modelo térmico.

En este tutorial se considerará la variación temporal del módulo de


elasticidad, con lo cual se deben realizar previamente el siguiente preproceso:

Para evitar problemas de convergencia, los tiempos en que se define el


módulo de elasticidad variable deben coincidir con los tiempos en los que se
definen las temperaturas nodales (archivo “.fil”). Con lo cual la primer tarea será
extraer alguna variable del modelo térmico en función del tiempo. Lo que nos
importará de este archivo serán los “time step”.

Lo siguiente será eliminar aquellos tiempos (en segundos) que se


encuentran repetidos, ya que al ingresar una curva de amplitud pueden existir
problemas. Se debe observar en el comienzo de cada “step”, ya que en estos
tiempos es donde suelen aparecer valores repetidos como se aprecia en la
siguiente tabla:

88
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

t""""""""""""""""""" t""""""""""""""""""""
[días] [s]
0.00 0.00
0.00 0.00
0.00 0.00
0.00 0.00
0.00 0.01
0.00 0.01
0.00 0.02
0.00 0.03
0.00 0.07
0.00 0.13
0.00 0.26
0.00 0.51
0.00 1.03
0.00 2.05
0.00 4.10
0.00 8.19
0.00 16.39
0.00 32.77
0.00 65.54
0.00 131.07
0.00 262.15
0.01 524.29
0.01 1048.58
0.02 2097.15
0.05 4194.30
0.10 8388.61
0.19 16777.20
0.38 32773.00
0.56 48768.80
0.75 64764.50
1.12 96756.10
1.86 160739.00
2.86 247139.00
3.00 259200.00
3.00 259200.00
3.00 259200.00
3.00 259200.00
3.00 259200.00
3.00 259200.00
3.00 259200.00
3.00 259200.00
3.00 259200.00
3.00 259200.00
3.00 259200.00
3.00 259200.00
3.00 259200.00
3.00 259201.00
3.00 259201.00
3.00 259202.00
3.00 259204.00
3.00 259208.00
3.00 259216.00
3.00 259233.00
3.00 259266.00
3.00 259331.00

Tabla 20

Luego de eliminar los tiempos repetidos, se debe dividir a la tabla por


cada step. A continuación se debe calcular el tiempo correspondiente a la edad
del hormigón de cada tongada (en segundos) con el objetivo de calcular luego,
el módulo de elasticidad para dichos tiempos. El tiempo correspondiente a la
edad de la primer tongada es coincidente con el tiempo total del análisis. En
cambio si se analiza la segunda tongada, por ejemplo, al tiempo total se deberá
restar los 3 primeros días (259200 segundos), ya que en dichos instantes la
esta tongada no existe. Algo similar se debe realizar con la tercer tongada,
restando al tiempo total los 6 primeros días (518400 segundos).

89
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

t*%%%%%%%%%%%%%%% t*%%%%%%%%%%%%%%% t*%%%%%%%%%%%%%%% t*%%%%%%%%%%%%%%% t*%%%%%%%%%%%%%%% t*%%%%%%%%%%%%%%%


t%%%%%%%%%%%%%%%%%%%% 1er% 1er% 2da% 2da% 3er% 3er%
[s] Tongada%%%%%%%%%%%%%%%%%%%%
Tongada%%%%%%%%%%%%%%%%%%%%
Tongada%%%%%%%%%%%%%%%%%%%%
Tongada%%%%%%%%%%%%%%%%%%%%
Tongada%%%%%%%%%%%%%%%%%%%%
Tongada%%%%%%%%%%%%%%%%%%%%
[s] [días] [s] [días] [s] [días]
0.00 0.00 0.00 # # # #
0.01 0.01 0.00 # # # #
0.02 0.02 0.00 # # # #
0.03 0.03 0.00 # # # #
0.07 0.07 0.00 # # # #
0.13 0.13 0.00 # # # #
0.26 0.26 0.00 # # # #
0.51 0.51 0.00 # # # #
1.03 1.03 0.00 # # # #
2.05 2.05 0.00 # # # #
4.10 4.10 0.00 # # # #
8.19 8.19 0.00 # # # #
16.39 16.39 0.00 # # # #
32.77 32.77 0.00 # # # #
STEP%1

65.54 65.54 0.00 # # # #


131.07 131.07 0.00 # # # #
262.15 262.15 0.00 # # # #
524.29 524.29 0.01 # # # #
1048.58 1048.58 0.01 # # # #
2097.15 2097.15 0.02 # # # #
4194.30 4194.30 0.05 # # # #
8388.61 8388.61 0.10 # # # #
16777.20 16777.20 0.19 # # # #
32773.00 32773.00 0.38 # # # #
48768.80 48768.80 0.56 # # # #
64764.50 64764.50 0.75 # # # #
96756.10 96756.10 1.12 # # # #
160739.00 160739.00 1.86 # # # #
247139.00 247139.00 2.86 # # # #
259200.00 259200.00 3.00 # # # #
259201.00 259201.00 3.00 1.00 0.00 # #
259202.00 259202.00 3.00 2.00 0.00 # #
259204.00 259204.00 3.00 4.00 0.00 # #
259208.00 259208.00 3.00 8.00 0.00 # #
259216.00 259216.00 3.00 16.00 0.00 # #
259233.00 259233.00 3.00 33.00 0.00 # #
259266.00 259266.00 3.00 66.00 0.00 # #
259331.00 259331.00 3.00 131.00 0.00 # #
259462.00 259462.00 3.00 262.00 0.00 # #
259724.00 259724.00 3.01 524.00 0.01 # #
STEP%2

260249.00 260249.00 3.01 1049.00 0.01 # #


261297.00 261297.00 3.02 2097.00 0.02 # #
263394.00 263394.00 3.05 4194.00 0.05 # #
267589.00 267589.00 3.10 8389.00 0.10 # #
275977.00 275977.00 3.19 16777.00 0.19 # #
284366.00 284366.00 3.29 25166.00 0.29 # #
301143.00 301143.00 3.49 41943.00 0.49 # #
317920.00 317920.00 3.68 58720.00 0.68 # #
351475.00 351475.00 4.07 92275.00 1.07 # #
418584.00 418584.00 4.84 159384.00 1.84 # #
485692.00 485692.00 5.62 226492.00 2.62 # #
518400.00 518400.00 6.00 259200.00 3.00 # #
518401.00 518401.00 6.00 259201.00 3.00 1.00 0.00
518402.00 518402.00 6.00 259202.00 3.00 2.00 0.00
518404.00 518404.00 6.00 259204.00 3.00 4.00 0.00
518408.00 518408.00 6.00 259208.00 3.00 8.00 0.00
518416.00 518416.00 6.00 259216.00 3.00 16.00 0.00
518433.00 518433.00 6.00 259233.00 3.00 33.00 0.00
518466.00 518466.00 6.00 259266.00 3.00 66.00 0.00
518531.00 518531.00 6.00 259331.00 3.00 131.00 0.00
518662.00 518662.00 6.00 259462.00 3.00 262.00 0.00
518924.00 518924.00 6.01 Tabla 21
259724.00 3.01 524.00 0.01
519449.00 519449.00 6.01 260249.00 3.01 1049.00 0.01
520497.00 520497.00 6.02 261297.00 3.02 2097.00 0.02
Lo siguiente, es calcular el módulo de elasticidad para la edad de cada
522594.00 522594.00 6.05 263394.00 3.05 4194.00 0.05
526789.00 526789.00 6.10 267589.00 3.10 8389.00 0.10
hormigón (columna t* en días) utilizando la expresión dada en la Ec. (54).
535177.00 535177.00 6.19 275977.00 3.19 16777.00 0.19
543566.00 543566.00 6.29 284366.00 3.29 25166.00 0.29
Todos aquellos valores que sean menores a 0.1 MPa pueden generar
560343.00 560343.00
577120.00 577120.00
6.49
6.68
301143.00
317920.00
3.49
3.68
41943.00
58720.00
0.49
0.68
610675.00 610675.00 7.07 351475.00 4.07 92275.00 1.07
problemas de convergencia en el modelo. Con lo cual se deben identificar
677784.00 677784.00 7.84 418584.00 4.84 159384.00 1.84
764184.00 764184.00 8.84 504984.00 5.84 245784.00 2.84
850584.00 850584.00 9.84 591384.00 6.84 332184.00 3.84
936984.00 936984.00 10.84 677784.00 7.84 418584.00 4.84
1023380.00 1023380.00 11.84 764180.00 8.84 504980.00 5.84 90
1109780.00 1109780.00 12.84 850580.00 9.84 591380.00 6.84
1196180.00 1196180.00 13.84 936980.00 10.84 677780.00 7.84
1282580.00 1282580.00 14.84 1023380.00 11.84 764180.00 8.84
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

estos valores y cambiarlos por el menor valor que sigue en la tabla como se
muestra a continuación:

Ec##########################
Ec##########################
Ec##########################
1er# 2da# 3er#
Tongada#####################
Tongada#####################
Tongada#####################
[Mpa] [Mpa] [Mpa]
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
0.14 0.14 0.14
5.51 0.14 0.14
73.19 0.14 0.14
455.75 0.14 0.14
1660.99 0.14 0.14
4144.83 0.14 0.14
7912.74 0.14 0.14
12337.16 0.14 0.14
15088.43 0.14 0.14
17030.44 0.14 0.14
19677.84 0.14 0.14
22764.42 0.14 0.14
25098.32 0.14 0.14
25339.81 0.14 0.14
25339.83 0.14 0.14
25339.85 0.14 0.14
25339.89 0.14 0.14
25339.97 0.14 0.14
25340.12 0.14 0.14
25340.45 0.14 0.14
25341.09 0.14 0.14
25342.35 0.14 0.14
25344.89 0.14 0.14
25349.97 5.50 0.14
25360.12 73.28 0.14
25380.31 455.68 0.14
25420.38 1660.80 0.14
25499.32 4145.05 0.14
25652.42 7912.66 0.14
25799.57 10537.48 0.14
26077.46 14044.52 0.14
26335.69 16364.43 0.14
26801.89 19373.32 0.14
27577.66 22715.77 0.14
28202.18 24647.28 0.14
28465.74 25339.81 0.14
28465.75 25339.83 0.14
28465.75 25339.85 0.14
28465.77 25339.89 0.14
28465.80 25339.97 0.14
28465.86 25340.12 0.14
28465.99 25340.45 0.14
28466.25 25341.09 0.14
28466.75 25342.35 0.14
28467.76 25344.89 0.14
28469.78 25349.97 5.50
28473.81 25360.12 73.28
28481.86 25380.31 455.68
28497.90 Tabla 25420.38 22 1660.80
28529.72 25499.32 4145.05
28592.32 25652.42 7912.66
28653.60 25799.57 10537.48
28772.38 26077.46 14044.52
28886.41 26335.69 16364.43 91
29101.42 26801.89 19373.32
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Lo siguiente es confeccionar una tabla en donde figuran E, µ y en la que


se introduce la denominada “Field variable”. Ésta es una variable que
relacionará la propiedad mecánica (módulo de elasticidad) con la curva de
amplitud del tiempo que se mencionará más adelante. Esta tabla es la que se
debe ingresar dentro de las propiedades mecánicas como se describirá más
adelante.

1er#Tongada

Ec###################
μ FV1
[Mpa]

0.14 0.25 1.00


0.14 0.25 2.00
0.14 0.25 3.00
0.14 0.25 4.00
0.14 0.25 5.00
0.14 0.25 6.00
0.14 0.25 7.00
0.14 0.25 8.00
0.14 0.25 9.00
0.14 0.25 10.00
0.14 0.25 11.00
0.14 0.25 12.00
0.14 0.25 13.00
0.14 0.25 14.00
0.14 0.25 15.00
0.14 0.25 16.00
0.14 0.25 17.00
5.51 0.25 18.00
73.19 0.25 19.00
455.75 0.25 20.00
1660.99 0.25 21.00

Tabla 23

Por último se debe realizar una última tabla en donde están vinculados
el tiempo y la “Field variable”. Esta es la curva de amplitud que se debe
ingresar en el programa.

AMP$FV1

t$$$$$$$$$$$$$$$$$$$$$$
FV1
[s]

0.00 1.00
0.01 2.00
0.02 3.00
0.03 4.00
0.07 5.00
0.13 6.00
0.26 7.00
0.51 8.00
1.03 9.00
2.05 10.00

Tabla 24

92
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Evolución#del#Módulo#de#Elas5cidad#para#
cada#hormigón#
35000"

30000"

25000"

Ec#[MPa]#
20000" 1er#Tongada#
2da#Tongada#
15000"
3er#Tongada#
10000"

5000"

0"
0" 5" 10" 15" 20" 25" 30"
t#[días]#

Gráfico 22

Para ingresar las propiedades se debe proceder de la siguiente forma:

Materials à Name: H1

Generalà Mass Density: 2.51E-05

Mechanicalà Elasticityà Elasticà Type: Isotropicà Number of field variables: 1à Copiar
y pegar la tabla de propiedades mecánicas correspondiente a la primer tongada.

Mechanicalà Expansionà Type: Isotropicà Expansion coeff. alpha: 1.05E-05.

Figura 85

93
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Se debe repetir la operación para crear los materiales correspondientes


a las demás tongadas (H2 y H3), con sus correspondientes evoluciones
temporales del módulo de elasticidad.

Luego se debe generar una curva de amplitud denominada AMP FV1,


de la misma forma que se describió en el modelo térmico para generar la
temperatura ambiente por ejemplo.

4.5.1.3 Sections

En este apartado se deben crear las tres secciones correspondientes a


cada material (SECH1, SECH2, SECH3) y asignarlas a las partes que les
correspondan de la misma forma que se explicó para el modelo térmico.

4.5.1.4 Assembly

No es necesario realizar cambios en el ensamble de las partes en el


modelo mecánico para este caso.

4.5.1.5 Steps

Se deben generar los 4 “steps” al igual que en el modelo térmico con la


misma duración. En este caso, el tipo de análisis será “Static, General”, se
pude tildar la opción “Nlgeom” (nonlinear geometry) para tener en cuenta
grandes desplazamientos en subsecuentes “steps”.

En cuanto a los incrementos, tanto el inicial como el mínimo deben ser


menores que aquellos ingresados en el modelo térmico por cuestiones de
convergencia. Por ejemplo en este caso se considera: 1E-05. Para el máximo
incremento se debe considerar en qué intervalos de tiempo se desea obtener
valores tensionales, (en este caso se sigue considerando 1 día, 86400
segundos).

4.5.1.6 Interactions

Como se explicó en un comienzo, no se utilizarán interacciones para


conectar las distintas partes del modelo, sino que se utilizarán las denominadas
“Tie-Constraints”.

94
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4.5.1.7 Constraints

Para dar continuidad al campo de tensiones y deformaciones en los


contactos entre las distintas partes, se debe proceder de la siguiente forma:

Constraints

Name: COM12à Continue…

Choose the master type: Surfaceà CONT M12 MASTER (Seleccionar la superficie
“master”)à Choose the slave type: Surfaceà CONT M12 SLAVE (Seleccionar la
superficie “slave”)

Discretization method: Analysis defaultà Position tolerance: Use computed defaultà Ok.

Figura 86

Figura 87

95
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4.5.1.8 Loads

En este tipo de análisis termo-mecánicos no se tienen en cuenta cargas


mecánicas como pueden ser el peso propio, sobrecargas útiles, sismos etc., ya
que se desean obtener valores tensionales únicamente relacionados al
desarrollo del calor de hidratación.

4.5.1.9 BCs -Boundary Conditions-

Como se describió en un principio, el muro posee anclajes sobre el


macizo rocoso ubicados cada 1 m. Lo primero será generar 3 “sets” de nodos
en correspondencia con los anclajes de cada una de las partes: ANC M1, ANC
M2, ANC M3. Se debe prestar atención al definir los nodos, ya que pueden
existir problemas denominados “overconstraining” si coinciden nodos
asignados al mismo tiempo a una “tie-constraint” y a una “BC”.

Figura 88

Para asignar las correspondientes condiciones de borde (cuya


nomenclatura es la misma que los “sets” de nodos) se debe hacer lo siguiente:

96
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

BC’s

Name: ANC M1à Step: Step-1 (Seleccionar el step correspondiente en el que


“aparece”ccada parte)à Category: Mechanicalà Types for Selected Step:
Displacement/Rotationà Continue…

Sets…à ANC M1à Continue…

Distribution: Uniformà U1=0, U2=0à Amplitude: (Ramp)à Ok

Figura 89

Figura 90

97
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 91

La mismo se debe hacer para las condiciones de borde restantes en los


demás “steps”.

La parte inferior del muro se considerará empotrada. Para ello se debe


generar otra “BC” que se llamará EMP M1. En este caso en vez de seleccionar
un set de nodos, se elegirá desde el “viewport” el “edge” inferior de la
“instance” M1. El “step” correspondiente es el Step-1.

4.5.1.10 Predefined Fields

En este apartado se explicará cómo introducir las temperaturas nodales


(archivo .fil), que se produjo como output en el modelo térmico.

Predefined Fields

Name: NT STEP1à Step: Step-1 (Seleccionar el step correspondiente)à Category:


Otherà Types for Selected Step: Temperatureà Continue…

Select regions for the field: Seleccionar la/s instance/s correspondientes (en este caso
será la instance M1, en el step 2 serán las instances M1 y M2, y en los steps 3 y 4 serán
las instances M1, M2 y M3à Done

Distribution: From results or output database fileà File name:


JOBTERESC001.filà Begin Step: 1 (Si se trata de las temperaturas nodales para el step-
2, corresponde ingresar en este punto el número 2)à End Step: 1 (vale la misma
aclaración hecha para la opción anterior)à Mesh compatibility: Compatibleà Ok

98
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 92

Figura 93

99
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 94

Las opciones “Begin step” y “End step” se refieren al “step” a partir del
cual comienza y finaliza respectivamente la lectura de temperaturas nodales
por parte del programa. Las opciones “Begin increment” y “End increment” se
refieren a lo mismo pero para el caso de los incrementos, al no ingresar esta
información el programa siempre leerá desde el primer incremento hasta el
último en cada step.

Lo mismo se debe repetir para las temperaturas nodales restantes:

NT STEP2, NT STEP3, NT STEP4

4.5.1.11 Mesh

Ya se explicó que las mallas de ambos modelos térmico y mecánico


deben ser las mismas (tamaño y forma de los elementos, número de nodos y
de elementos, etc.). Sin embargo, se debe cambiar el tipo de elementos de la
librería para poder realizar el análisis mecánico.

Module: Mesh

Assign Element Type à Select the regions to be assigned element types


(Seleccionar todas las instances)à Done

Element Library: Standardà Family: Plane strainà Geometric order: Linear

100
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

CPE4R: A 4-node bilinear plane strain quadrilateral, reduced integration, hourglass


control.

CPE3: A 3-node linear plane strain triangle.

4.5.1.12 Field Variable

Para concluir el input del modelo mecánico, se deben relacionar dos


ítems definidos anteriormente (Propiedades Mecánicas-Field Variable y
Tiempo-Field Variable) a partir de la edición del “keywords”:

Modelà Edit Keywordsà MOD MEC ESC 001à Scrollear hasta hallar la definición de
cada step (duración, incremenetos máximos y mínimos, etc.)à Escribir el siguiente
comando para cada instance:

*FIELD, VARIABLE=(Number of Field Variable), Amplitude=(Amplitude Name)

(Instance Name).(Element Set Name), 1

Para este caso resulta:

*FIELD, VARIABLE=1, AMPLITUDE= AMP FV1

M1-1._PickedSet3,1

*FIELD, VARIABLE=1, AMPLITUDE= AMP FV1

M2-1._PickedSet3,1

*FIELD, VARIABLE=1, AMPLITUDE= AMP FV1

M3-1._PickedSet3,1

Los datos acerca del nombre de la “instance” y del set de elementos


correspondiente a la sección asignada a cada parte se encuentran en el mismo
“keywords” cerca del comienzo:

Figura 95

101
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 96

4.5.1.13 Job

De la misma forma que se procedió para el modelo térmico, se debe


crear el “job” para este modelo cuya nomenclatura será: JOBMECESC001. Por
último se debe correr la simulación.

4.5.2 Output

En las siguientes imágenes se presentan las tensiones S11 para los


distintos “steps”. En primer lugar se debe verificar que los resultados
tensionales sean concordantes con los resultados térmicos. Es decir, que no
haya tensiones en el primer “step” para partes que “aparecen” en el segundo
“step".

Las zonas en las cuales se extraerán los datos tensionales para su


análisis son aquellas que se corresponden con los fenómenos de gradiente de
masa y de superficie. Además, es conveniente analizar en el “viewport” para
los distintos “time step”, si existe alguna singularidad (tensiones de tracción),
generada por alguna condición de borde u otra causa.

102
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Figura 97

Figura 98

Figura 99

103
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

4.5.2.1 Historia de tensiones

A continuación se presentan las tensiones obtenidas como resultado del


modelo mecánico para las zonas superiores de cada parte:

Historia"de"tensiones"para"la"1er"tongada"
4.5%

4.0%

3.5%

3.0%
S"[MPa]"

2.5%

2.0% S11"
1.5% f'tm"
1.0%

0.5%

0.0%

!0.5%
0%

10%

20%

30%

40%

50%

60%

70%

80%

90%

100%
t"[días]"

Gráfico 23

En el siguiente gráfico se muestra una “ampliación” de los primeros días


de la evolución tensional para ver el instante en el que se estaría superando la
resistencia a tracción (aproximadamente a partir del día 1):

4.0#

3.5#

3.0#

2.5#
S"[MPa]"

2.0#

1.5#
S11"
1.0# f'tm"

0.5#

0.0#
0.5#

1.5#

2.5#
0.0#

1.0#

2.0#

3.0#

t"[días]"

Gráfico 24

Para la segunda y tercer tongada se produce algo similar como se


puede observar en los siguientes gráficos:

104
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Historia"de"tensiones"para"la"2da"tongada"
4.5%

4.0%

3.5%

3.0%

S"[MPa]"
2.5%

2.0% S11"
1.5% f'tm"
1.0%

0.5%

0.0%

!0.5%
0%

10%

20%

30%

40%

50%

60%

70%

80%

90%

100%
t"[días]"

Gráfico 25

Historia"de"tensiones"para"la"3er"tongada"
4.5#

4.0#

3.5#

3.0#
S"[MPa]"

S11"
2.5#
f'tm"
2.0#

1.5#

1.0#

0.5#

0.0#
0#

10#

20#

30#

40#

50#

60#

70#

80#

90#

100#

t"[días]"

Gráfico 26

Para extraer las tensiones en función del tiempo se procede de la


siguiente forma:

Create X,Y Data

Sourceà ODB field output

Variablesà Position: Centroidà S:Stress Componentsà S11, SMÁX, SMISES

Elements/Nodesà Edit Selectionà Method: Pick from viewport

Los valores tensionales más precisos serán aquellos que se obtienen en


los puntos de integración, luego siguen los valores en el centroide del elemento

105
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

y por último, los menos “confiables” son aquellos valores que figuran en los
nodos.

4.5.2.2 Perfil de tensiones

Puede ser conveniente contar con valores tensionales en los mismos


perfiles y en los mismos instantes de tiempo que el modelo térmico. Se debe
proceder de la misma forma que se explicó anteriormente: generar los “path”,
establecer los “time step” correspondientes, etc. Se debe recordar que en el
“viewport”, la variable que se está mostrando es la que se extraerá como dato
(S11, por ejemplo).

t%%%%%%%%%%%%%%%%%% t*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Step%Time%%%%%%%%%%%%%%%%%%
Día STEP FRAME%%%%%%%%%%%%%%%%%
[s] [s] [s]

1 8.64E+04 Step-1 8.64E+04 9.68E+04 56


3 2.59E+05 Step-2 0.00E+00 0.00E+00 0
7 6.05E+05 Step-3 8.64E+04 9.68E+04 56
14 1.21E+06 Step-3 6.91E+05 6.49E+05 63
28 2.42E+06 Step-3 1.90E+06 1.94E+06 78
100 8.64E+06 Step-4 3.46E+06 3.46E+06 96

Tabla 25

No se pueden definir “path” con elementos, con lo cual los perfiles


tensionales corresponderán a valores nodales.

Peril#Horizontal#de#Tensiones#####################
1er#Tongada#
0.25%
0.20%
0.15%
0.10%
S11#[MPa]#

0.05%
0.00%
!0.05%
!0.10% Día#1#
Día#3#
!0.15% Día#7#
Día#14#
!0.20%
Día#28#
!0.25% Día#100#
!0.30%
0.25%

0.75%

1.25%

1.75%

2.25%
0.00%

0.50%

1.00%

1.50%

2.00%

X#[m]#

Gráfico 27

106
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

Peril"Ver0cal"de"Tensiones""""""""""""""""""""""""
1er"Tongada"
5.00#
4.50#
4.00#
Día"1"
3.50#
Día"3"

Y"[m]"
3.00# Día"7"
Día"14"
2.50# Día"28"
Día"100"
2.00#
1.50#
1.00#
0.50#
0.00#
)10#
)9#
)8#
)7#
)6#
)5#
)4#
)3#
)2#
)1#
0#
1#
2#
3#
4#
5#
S11"[MPa]"

Gráfico 28

107
TUTORIAL
Análisis computacional de problemas termo-mecánicos en el hormigón

5 SEGUNDA PARTE: ANÁLISIS NO LINEAL

108

También podría gustarte

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy