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

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
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
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
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
• 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
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
• ∂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
∫ −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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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.
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
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
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
Figura 33
Figura 34
39
TUTORIAL
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
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
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
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).
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
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
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
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
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
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
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
Figura 43
Figura 44
Figura 45
50
TUTORIAL
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
Figura 46
Figura 47
52
TUTORIAL
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
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
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
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
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
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
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
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
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
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
*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
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
NT
*OUTPUT, FIELD
64
TUTORIAL
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
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
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
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
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
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
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
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
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
Figura 78
Figura 79
74
TUTORIAL
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
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
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
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
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
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
• 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
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
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
!
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
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
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
!
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
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
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
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
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
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
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
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
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
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
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
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
Figura 92
Figura 93
99
TUTORIAL
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
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
M2-1._PickedSet3,1
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
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
Figura 97
Figura 98
Figura 99
103
TUTORIAL
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
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
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
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
5 SEGUNDA PARTE: ANÁLISIS NO LINEAL
108

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

Cargado por

Información del documentohacer clic para expandir la información del documento

Copyright:

Formatos disponibles

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

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

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

Cargado por

Copyright:

Formatos disponibles

TUTORIAL

Análisis computacional de problemas

Universidad Nacional de La Plata

Ing. Rafael Lamattina

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

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

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

Las estructuras se encuentran en contacto con los fluidos del ambiente

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

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

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

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

2.2 Hormigón masivo

2.2.1 Estructuras de hormigón masivo

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

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

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

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

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.

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

Figura 3: Construcción de una presa HCR.

• El espesor de las tongadas de hormigón.

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

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

3.1 Ecuaciones de Transferencia de Calor

3.1.1 Relación constitutiva

Las observaciones experimentales indican que el flujo de calor en un

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Aplicando el Teorema de la Divergencia:

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

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

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

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

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

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

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

Utilizando la aproximación de Galerkin, las temperaturas variacionales