Investigaciones y Estudios

HORMIGÓN Y ACERO | 67

Generación automática de esquemas

de bielas y tirantes considerando
criterios constructivos
Automatic generation of strut-and-tie schemes adjusted
for constructability
Jesús Miguel Bairán García(1)

Recibido | Received: 01-10-2011

Aceptado | Accepted: 29-12-2011

Resumen

E l método de bielas y tirantes (BT) es de gran utilidad en el dimensionamiento a rotura de elementos de hor-
migón, particularmente en casos de regiones D donde las hipótesis clásicas de vigas y placas no son aplica-
bles. Aunque el método viene recogido en muchas normativas modernas, existen al menos tres aspectos que li-
mitan un uso más extendido en la práctica: no da información en servicio, la selección de esquemas factibles no
es trivial en casos complejos y la solución no es única. En este artículo se propone un método que genera esque-
mas factibles de BT mediante el uso de análisis lineales y criterios de optimización para modificar la distribu-
ción de propiedades del material en la región de estudio. No es necesario disponer un armado inicial y permi-
te considerar criterios constructivos; en particular, se aborda el de armado ortogonal. Se incluyen casos de estu-
dio y ofrecen conclusiones y perspectivas futuras.

Palabras clave: regiones D, método de bielas y tirantes, criterios constructivos, optimización topológica, campos
de tensiones

Abstract*

S trut-and-Tie (ST) method is very useful in the Ultimate Limit State design of reinforced concrete elements, particular-
ly in the case of D regions where beam and plate theories cannot be applied. In spite of being included in most modern
design codes, there are at least three issues that limit a more extended use of this method in practice: it doesn’t provide in-
formation in service, the selection of suitable schemes is not trivial for complex cases and the solution isn’t unique. In this
paper a novel method to produce suitable ST schemes is proposed by combining linear elastic analyses and optimization cri-
teria. It is based on suitable modifications of the distribution of elastic parameters in the structural model. Without requir-
ing initial distribution of reinforcement, it is capable of considering constructability conditions, particularly the condition
of orthogonal reinforcement layout is considered in this paper. A series of case-studies are presented and analyzed; finally,
conclusions are drawn.

Keywords: D-regions, strut-and-tie method, constructability, topological optimization, stress fields.

* An extensive English language summary of the present article is provided on page 94 of this issue both for the
convenience of non-Spanish-speaking readers and inclusion in databases.

(1) Dr. Ingeniero de Caminos, Canales y Puertos. Prof. Lector. Universitat Politècnica de Catalunya (Barcelona, España).

Persona de contacto / Corresponding author: jesus.miguel.bairan@upc.edu

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689

Investigaciones y Estudios 68 | HORMIGÓN Y ACERO Generación automática de esquemas de bielas y tirantes considerando criterios constructivos

1. INTRODUCCIÓN las situaciones de Estados Límite de Servicio (E.L.S.)

y Estados Límite Últimos (E.L.U.), involucrando cál-
La mayor parte de las estructuras pueden ser asimila- culos en rotura, lo cual sólo puede hacerse con méto-
das a una serie de elementos prismáticos o superficia- dos no-lineales, más complejos y que requieren cono-
les en los que son aplicables las teorías clásicas de vi- cer las armaduras con anterioridad.
gas y placas. En éstas se postulan hipótesis sobre la dis-
tribución de deformaciones en las secciones transver- El método de bielas y tirantes (BT) fue propuesto for-
sales; p.e. las teorías de Navier-Bernoulli y Timoshen- malmente en la década de 1980 por Schlaich et al [1]
ko en el caso de vigas, o su contrapartida de Love- como un método de dimensionamiento de carácter ge-
Kirchoff y Reissner-Mindlin en el caso de placas. Estas neral aplicable a cualquier tipo de estructuras y, como
hipótesis simplifican notablemente el cálculo y se pue- caso particular, a regiones D. En este método se apro-
den desarrollar formulaciones analíticas para analizar xima el comportamiento resistente al de una celosía
o dimensionar en situaciones de servicio o rotura. equivalente que explica el flujo de sus fuerzas y el
equilibrio interno, ya sea de forma estable o indiferen-
Sin embargo, ya sea por las condiciones geométricas o te. La celosía equivalente se compone de elementos a
de aplicación de las cargas, existen regiones de la es- compresión (bielas), elementos a tracción (tirantes) y la
tructura donde no son aplicables las teorías clásicas. intersección entre ellos (nudos). Las bases del método
De esta forma, se identifican dos tipos de regiones en y las condiciones que debe reunir la celosía equivalen-
una estructura. Se conoce como región-B (“Beam re- te se resumen en el apartado 2; en las referencias [2, 3,
gion”), o región viga, las zonas o elementos donde se 4, 5] se puede encontrar una exposición más detallada.
pueden aplicar las teorías clásicas. Las zonas o elemen- En las referencias [3, 4, 5, 6] se incluyen además nume-
tos donde no se pueden aplicar las teorías clásicas se rosos casos prácticos de aplicación.
conocen como zonas de discontinuidad o región-D
(“Disturbed region”). En la Figura 1a se identifican las El método BT puede considerarse como una genera-
diferentes regiones en una estructura porticada típica. lización de otras teorías clásicas, como la analogía de
Las regiones D pueden ser zonas localizadas de estruc- la celosía o los campos de compresiones, donde las
turas lineales, e.g. apoyos a media madera, ménsulas, bielas y tirantes representan las resultantes de los
zonas de anclaje de pretensado o próximas a la aplica- campos de compresión y tracción respectivamente.
ción de cargas, nudos, etc. o bien constituir toda una Esta representación es la considerada en la analogía
estructura global o elemento cuando la geometría no de la celosía de Ritter-Mörsch,tradicionalmente usa-
permite su asimilación a un prisma o placa, como la re- das para explicar los mecanismos resistentes a cor-
presentada en la Figura 1b. En las zonas locales, la ex- tante y torsión en el hormigón fisurado. Actualmente
tensión de la región D está limitada, en virtud del el método BT es aceptado por la gran mayoría de
Principio de Saint-Venant, a dimensiones próximas al normativas de proyecto, entre ellas la EHE-08 [7] en
espesor del elemento. la que se usa como método general de proyecto en
rotura. En la literatura se han propuesto otros méto-
En general, el cálculo de las regiones D debe hacerse dos similares que mantienen el concepto de campo
considerando el problema tridimensional de equili- en la distribución de tensiones de tracción y compre-
brio de un medio continuo. Actualmente existen nu- sión, para aplicaciones similares, ver Ref. [8].
merosas herramientas comerciales de cálculo que
son capaces de resolver dicho problema con cierta fa- Sin embargo, existen al menos tres aspectos que limi-
cilidad, al menos en rango lineal elástico. No obstan- tan un uso más extendido en la práctica: no es apli-
te, las soluciones obtenidas no son directamente apli- cable en servicio, la elección de esquemas factibles
cables para el hormigón armado, ya que éste puede no es trivial en casos complejos y la solución no es
fisurar para cargas relativamente pequeñas. Más única. Los dos últimos puntos son cruciales para la
aún, de acuerdo a las normativas actuales, basadas estimación de los esfuerzos internos a ser resistidos y
en el método de los estados límite, es necesario estimar el posterior dimensionado de las armaduras. Es por

a) b)

Figura 1. Regiones D en estructuras porticadas (a) y en elementos completos (b)

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García
 Investigaciones y Estudios
Generación automática de esquemas de bielas y tirantes considerando criterios constructivos HORMIGÓN Y ACERO | 69

ello que en la literatura se han propuesto diferentes ner un armado inicial. Esta estrategia tiene la venta-
metodologías para poder obtener esquemas resisten- ja de que en función de los criterios energéticos usa-
tes de forma racional y automática. Entre ellos pue- dos y de distribución de propiedades mecánicas se
den mencionarse trabajos en los que a partir de una pueden obtener esquemas alternativos o introducir
distribución de armaduras conocida (distribución de condicionantes constructivos. En particular, utilizan-
tirantes) se obtiene la distribución del campo de do propiedades de materiales isótropos y ortótropos
compresiones (bielas) necesario para la estabilidad. es posible producir esquemas con armados ortogo-
Este tipo de métodos se basa, en general, en una se- nales o no. Adicionalmente se presentan una serie de
rie de análisis no-lineales, más o menos simplifica- ejemplos de aplicación que permiten evaluar las bon-
dos, en los que al menos se debe considerar la posi- dades del método, analizar su comportamiento y
bilidad de fisuración del hormigón a tracción. Si el prestaciones. Finalmente se ofrecen conclusiones y
armado es adecuado, el análisis no-lineal converge perspectivas.
con una distribución de tensiones en el hormigón
que reproduce el campo de compresiones. Entre los
trabajos que aplican esta estrategia cabe citar las re- 2. BASES DEL MÉTODO DE BIELAS Y TIRANTES
ferencias [9, 10].
El método BT se basa en el Teorema del Límite In-
La estrategia anterior es especialmente útil para veri- ferior de la Teoría de la Plasticidad, el cual establece
ficar armados propuestos. En el caso de que se bus- que si una estructura se carga con una solicitación Q*
que todo el esquema BT, incluyendo la disposición y se puede encontrar una distribución interna de ten-
de tirantes, se ha recurrido a métodos de optimiza- siones que satisfaga las siguientes condiciones: 1) en
ción topológica en los que se busca una forma estruc- ningún punto se supera el límite de plastificación del
tural de mínimo peso o bien de menor energía de de- material, 2) satisface el equilibrio interno y externo y
formación, ver [11]. En la Ref. [12] se propone un mé- 3) satisface las condiciones de contorno, entonces se
todo que, partiendo de una región sólida modeliza- tiene que la carga Q* es un límite inferior de la carga
da mediante elementos finitos, se van eliminando los de rotura de la estructura, Ec. (1).
elementos que contribuyen menos al sistema, uno
por cada paso de iteración, hasta que la forma resul-
tante converge a un esquema que se asemeja a una (1)
celosía; una metodología similar se sigue en [13]. En
la referencia [14] se plantea el problema de obtener la Del teorema anterior puede deducirse que aunque el
forma estructural de menor energía pero, en lugar de método BT no arroja soluciones exactas, éste es siem-
usar una representación en continuo (sólido), consi- pre conservador, lo cual lo hace muy atractivo para
dera directamente la representación en forma de es- el proyecto de estructuras. Sin embargo, la validez
tructura de barras articuladas investigando las dife- del teorema no puede garantizarse matemáticamen-
rentes conectividades posibles entre los nudos. te en el caso de las estructuras de hormigón armado
por las siguientes razones, Ref. [16]:
Las técnicas basadas en la optimización topológica
tienen la ventaja de que no necesitan proponer una • La ductilidad del hormigón es limitada y no pue-
armadura inicial y que el esquema final es una solu- de equipararse a un material plástico perfecto co-
ción óptima en términos de peso o energía. Sin em- mo se asume en la teoría de plasticidad. Además,
bargo, puede arrojar soluciones no convenientes des- algunos esquemas de BT, aunque cumplan los re-
de el punto de vista constructivo o de disposición de quisitos del teorema, no se pueden desarrollar al
armaduras, p.e. con armaduras inclinadas que difi- necesitar mucha ductilidad en el hormigón.
cultan la construcción.
• La resistencia a compresión del hormigón no es
Por otro lado, cabe destacar que se han realizado es- constante, sino que depende de las tensiones y de-
tudios experimentales para evaluar la validez del formaciones laterales.
método BT en el caso del hormigón armado y delimi-
tar su campo de aplicación y puntos débiles. En ese • Existe la posibilidad de rotura por anclaje, la cual
sentido, cabe mencionar la referencia [15], así como es de carácter frágil.
las recopiladas en [5].
La puesta en práctica del método requiere comple-
En este artículo se propone un método de generación mentar los requisitos formales con otros de carácter
de bielas y tirantes, basado en optimización topoló- técnico y de buena práctica que garanticen que la
gica que utiliza análisis lineales elásticos de la región carga que se está estimando es segura [2, 3, 4, 5, 6].
D en estudio. El método se basa en modificaciones Éstos incluyen la definición de la capacidad resisten-
sucesivas de la distribución de propiedades mecáni- te de las bielas, tirantes y nudos para las diferentes
cas en la estructura de acuerdo a criterios energéti- situaciones [7], las comprobaciones de anclajes de ar-
cos. Al final del proceso, el campo de tensiones prin- madura y recomendaciones sobre los ángulos míni-
cipales de compresión y tracción permite definir el mos de inclinación entre bielas y tirantes para garan-
esquema de bielas y tirantes, sin necesidad de dispo- tizar que la compatibilidad de deformaciones del hor-

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García
Investigaciones y Estudios 70 | HORMIGÓN Y ACERO Generación automática de esquemas de bielas y tirantes considerando criterios constructivos

Figura 2. Límites de la región D en el caso de un nudo de un pórtico, diagrama de cuerpo libre y sistema de cargas equivalentes

migón pueda cumplirse sin excesivas demandas de la región de estudio y aislar el cuerpo libre de la mis-
ductilidad en el hormigón [4]. ma, como se representa en las Figuras 2a y 2b. Pos-
teriormente, en virtud del Principio de Saint-Venant,
En la definición del esquema de BT no existen reglas se sustituyen los esfuerzos internos en los planos de
generales, ya que en principio cualquier esquema corte del cuerpo libre por cargas puntuales estática-
compatible con las condiciones geométricas y de mente equivalentes como se aprecia en la Figura 2c.
apoyos que explique el mecanismo de transmisión
de las cargas (a veces denominado flujo de cargas) es
El siguiente paso consiste en realizar un cálculo line-
posible. Algunas recomendaciones de buena práctica
al elástico de la región aislada. Para ello se ha de-
son: intentar no elegir esquemas muy diferentes de
sarrollado un programa de análisis mediante el mé-
la solución elástica, intentar usar celosías isostáticas
todo de elementos finitos. El programa utiliza ele-
para evitar que la solución dependa de la compatibi-
mentos triangulares de 3 nodos de tipo lagrangeano
lidad de deformaciones y elegir los esquemas de me-
con interpolación lineal entre nodos y deformación
nor energía de deformación. Teniendo en cuenta que
constante. En esta etapa la respuesta a nivel material
las bielas (hormigón) suelen ser mucho más rígidas se considera perfectamente lineal, elástica, isótropa y
que los tirantes (armaduras) y que éstos últimos se homogénea. De igual forma, en esta etapa se constru-
dimensionarán para trabajar a una tensión más o me- ye una malla de elementos finitos de la región D que
nos constante (fyd≤ 400 MPa), aproximadamente se será la misma en los pasos siguientes.
obtiene el esquema de menor energía de deforma-
ción con la menor longitud de tirantes.
A partir de la información que proporciona la solu-
A pesar de la relativa libertad que existe en la elec- ción elástica, se definen una serie de factores de efi-
ción del esquema de BT, este paso puede influir nota- cacia de material que se usarán como base para mo-
blemente en el diseño final. En el apartado siguiente dificar las características del material en cada punto.
se propone un método que obtiene esquemas de bie- Los factores de eficacia cuantifican la importancia re-
las y tirantes en regiones de geometría cualquiera. lativa de cada elemento y, opcionalmente, identifica-
rá su modo de trabajo (tracción, compresión). A par-
Se debe reconocer que el término “flujo”se suele usar tir de aquí se inicia un proceso iterativo cuyo esque-
de forma incorrecta ya que el esquema BT no repre- ma general se muestra en la Figura 3. En cada itera-
senta realmente la integración de un campo vecto- ción, tras modificar la distribución de propiedades
rial. Sin embargo, en este artículo se mantendrá este elásticas en la región, se repite el análisis elástico.
término por su amplio uso en este campo. En ese
sentido, se hará referencia al flujo de tracción y com- En cada paso de iteración la nueva distribución de
presión asociado a dos campos vectoriales de tensio- propiedades elásticas tiende a reducir la energía de
nes principales cuyos módulos son respectivamente deformación (W) de la estructura. Se considera que
las magnitudes de las tensiones principales (σ11 y σ22) el proceso converge cuando entre dos iteraciones se-
y su dirección viene definida por las direcciones princi- guidas, el cambio relativo en la energía de deforma-
pales. Las bielas y los tirantes, representan la integra- ción de la estructura, calculado como se indica en la
ción (resultante) de dichos campos vectoriales. Ec. (2), es menor que la tolerancia elegida. En los es-
tudios presentados en este artículo se ha usado una
3. MÉTODO PROPUESTO PARA GENERAR ESQUEMAS tolerancia de 1%.
DE BIELAS Y TIRANTES

3.1. Descripción general del proceso (2)

Al igual que en el método general BT, el primer paso En la ecuación anterior, q y d son respectivamente
en el estudio de una región D es definir la frontera de los vectores de cargas y desplazamientos nodales del

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García
 Investigaciones y Estudios
Generación automática de esquemas de bielas y tirantes considerando criterios constructivos HORMIGÓN Y ACERO | 71

el volumen de la estructura original. We es la densidad

de energía de cada elemento, como se muestra en la
Ec. (5.1); la cual aparece desarrollada en la Ec. (5.2) pa-
ra problemas 2D, en el sistema de coordenadas glo-
bal y en el de las direcciones principales.

(4)

(5.1)
Figura 3. Diagrama de flujo general del método de generación
de esquemas BT
(5.2)
modelo de la estructura. El esquema BT queda defini-
do al final de la iteración en función de la distribución
de las características mecánicas y el flujo de fuerzas Las propiedades mecánicas de cada elemento (i) se
internas. Variando los criterios con los que se modi- modifican al final de cada iteración como se indica
fican las propiedades mecánicas a lo largo del proce- en la Ec. (6). Se debe notar que la matriz de rigidez
so se pueden obtener esquemas de armados distin- del elemento seguirá siendo de forma isótropa. En
tos. No obstante, resulta evidente que no todos con- estas ecuaciones E y G son respectivamente el módu-
vergerán al mismo valor de energía mínima ya que al lo elástico y el módulo de cortante del material origi-
cambiar los criterios de toma de decisiones se aña- nal y Ee y Ge son las propiedades mecánicas modifi-
den condiciones nuevas al problema. cadas tras la iteración.

En los apartados siguientes se proponen dos criterios

de armado distintos: el primero tiende a minimizar Ee = fe E (6.1)
la energía sin condicionantes en la orientación del
campo de tracciones, que posteriormente será usado Ge = fe G (6.2)
para definir los tirantes. Por lo tanto, puede producir
tirantes inclinados. Este criterio se ha llamado crite- De forma similar, el volumen efectivo de material
rio “isótropo - sin consideraciones constructivas” ya queda definido como se indica en la Ec. (7), donde
que, como se verá, aplica el mismo tratamiento a las V(i) es el volumen de cada elemento.
zonas de compresión y a la de tracción y el tipo de
matriz constitutiva usado es siempre de esta forma.
El segundo criterio, usa matrices constitutivas ortó- (7)
tropas en las zonas de tracción para forzar orienta-
ciones determinadas de los tirantes, por lo que ha si-
do denominado como criterio “ortótropo - con consi- 3.3. Criterio 2: Ortótropo - con consideraciones
deraciones constructivas”. constructivas

El criterio anterior aplica igual tratamiento a los ele-

3.2. Criterio 1: Isótropo - sin consideraciones
mentos dentro de los campos de compresión (bielas)
constructivas
y los de tracción (tirantes). En este apartado se pro-
pone un segundo criterio que intenta discriminar en-
En el criterio isótropo se define el factor de eficiencia
tre los elementos que formarán parte de las zonas
de cada elemento como el cociente entre la densidad
comprimidas de aquellas que están principalmente
de energía (energía por unidad de volumen) del ele-
traccionadas, las cuales terminarán generando tiran-
mento y la densidad de energía de la estructura com-
tes. Una vez identificados los elementos que forma-
pleta, como se indica en la Ec. (3.1). Para garantizar
la estabilidad numérica, el factor de efectividad no rán parte de tirantes, es posible modificar las propie-
superara los límites superior e inferior que se indican dades mecánicas del material en esos puntos para
en la Ec. (3.2). potenciar ciertas direcciones en los campos de ten-
siones. En particular, se usarán matrices de rigidez
ortótropas para dirigir las tracciones en sentido ver-
(3.1) tical u horizontal con la intención de producir esque-
mas de armado que faciliten la construcción.

0,01 ≤ fe ≤ 1,0 (3.2) Para discriminar si un elemento debe ser considera-

do como parte de una biela o un tirante se recurre a
Donde Wm es la densidad de energía de toda la es- la Ec. (5) donde se observa que la energía de defor-
tructura, definida en la Ec. (4). En esta ecuación, V es mación es invariante con el sistema de coordenadas.

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García
Investigaciones y Estudios 72 | HORMIGÓN Y ACERO Generación automática de esquemas de bielas y tirantes considerando criterios constructivos

En el caso particular de las direcciones principales, se dez del elemento toma una forma ortótropa con di-
puede descomponer en dos términos: uno que englo- recciones principales del material predefinidas se-
ba la energía absorbida por las tensiones principales gún la dirección de armado deseado, p.e. para forzar
menores, que tienden a ser de compresión (Ec. 8), y las direcciones de las armaduras en sentido vertical y
otro que representa la energía absorbida por las ten- horizontal, la matriz de rigidez será como se indica
siones principales mayores, que tienden a ser de trac- en la Ec. (11).
ción (Ec. 9). La suma de ambos términos es igual a la
densidad de energía completa del elemento (Ec. 5.2).
(11)
(8)

(9) Donde fx y fy modulan la rigidez del elemento en las

dos direcciones de armado potenciando la existencia
de un solo tirante, horizontal o vertical, cuando el
Los casos particulares en que las tensiones principa-
trabajo tiene una dirección preferente. Por otro lado,
les son del mismo signo: compresión biaxial (σ11 es
el coeficiente feT afecta la matriz de rigidez en todas
de compresión) y tracción biaxial (σ22 es de tracción) las direcciones por igual en función de la importan-
no requieren un tratamiento especial ya que se con- cia del elemento. Estos factores se obtienen en fun-
siderarán como elementos del campo de compresio- ción de la energía de deformación, calculada en el
nes o de tracciones respectivamente. paso anterior como se indica en las las ecuaciones si-
La cantidad de trabajo en tracción del elemento se guientes:
cuantifica mediante el factor fT definido en la Ec. (10).
Así pues, es posible definir un umbral a* para deci- (12)
dir si un elemento formará parte de un campo de
tracciones o de compresiones. A partir de estudios de
sensibilidad iniciales se recomienda que este umbral (13)
se tome entre 0.7 y 0.9. En este artículo se ha tomado
un umbral de 0.8 para todos los ejemplos estudiados.
(14)
(10)

(15)
Elementos en el campo de compresiones: fT≤a*

Si la fracción de energía de tracción no supera el um- De las expresiones (14) y (15) se observa que los fac-
bral, (fT≤a*),se considera que el elemento trabaja fun- tores fx y fy son cantidades positivas cuya suma es
damentalmente a compresión y que forma parte de siempre la unidad. Por último, el coeficiente multi-
un campo de compresiones. En este caso se incluyen plicador a toda la matriz constitutiva (feT) se define
las zonas solicitadas a compresión biaxial y aquellos como se indica en la Ec. (16).
en los que la biela está sometida a tracción ortogonal
a la dirección de compresión principal. Este último
representa el caso de bielas en abanico o en forma de (16)
botella. Esta información podrá usarse en una posi-
ble etapa de comprobación de la resistencia de la bie- Se observa que la formulación propuesta incluye co-
la, aunque esto no es motivo del presente artículo. mo caso particular las zonas armadas con dos tiran-
tes oblicuos en las zonas de tracción biaxial (fx>0 y
Las propiedades mecánicas de los elementos dentro
fy>0). En las zonas de tracción en dos direcciones en
de los campos de compresiones serán modificadas,
la que prevalezca una dirección sobre la otra se ten-
de acuerdo al factor de eficiencia (Ec. 3), de la misma
drá que uno de los factores ( fxó fy ) tiende a cero y se
forma que en el criterio 1 (Ec. 6). Por lo que las zonas
tendrá un solo tirante.
de compresión tendrán comportamiento isótropo y
la dirección principal del campo podrá ser inclinada.
4. EJEMPLOS DE APLICACIÓN
Elementos en el campo de tracciones: fT>a* ó σ22>0
En este apartado se desarrollan y analizan tres ejem-
Si fT supera el umbral a*, o bien si la solicitación es de plos de aplicación con el primer objetivo de mostrar
tracción en las dos direcciones principales (biaxial), las capacidades del método propuesto y también pa-
se considera que el elemento trabaja fundamental- ra estudiar el funcionamiento y evolución del proce-
mente a tracción y que forma parte de un campo de so de generación de BT. En la Figura 4 se muestran
tracciones (tirante). En estos casos, la matriz de rigi- los esquemas estáticos de los tres ejemplos conside-

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García
 Investigaciones y Estudios
Generación automática de esquemas de bielas y tirantes considerando criterios constructivos HORMIGÓN Y ACERO | 73

Figura 4. Casos de estudio considerados (cotas en mm)

rados. Los dos primeros consisten en vigas de gran El método parte de la solución elástica del problema,
canto con una abertura, teniendo el primero una car- considerando material lineal e isótropo. La distribu-
ga superior y el segundo la carga colgada. El último ción de tensiones principales de tracción (σ11) y com-
ejemplo es una ménsula corta con una relación de as- presión (σ22) se muestra en la Figura 5 para los tres
pecto (L/h) de 0,43. En todos los casos se analizarán casos anteriores. En todos los casos se han usado ele-
los esquemas de armado obtenidos con los dos crite- mentos triangulares de tres nodos con interpolación
rios expuestos en el apartado 3. lineal. El número de elementos usado es de 1.624 para

σ11 σ22

σ11 σ22

σ11 σ22

Figura 5. Solución lineal elástica inicial del estado tensional para los tres casos de estudio (MPa)

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García
Investigaciones y Estudios 74 | HORMIGÓN Y ACERO Generación automática de esquemas de bielas y tirantes considerando criterios constructivos

Figura 6. Campos de tensiones principales (MPa) y distribución de factor de eficiencia (adimensional)

para el Ejemplo 1 - Criterio 1

los ejemplos 1 y 2 y de 1.032 para el ejemplo 3. Las rior del elemento que suben hacia los apoyos, con
propiedades elásticas iniciales consideradas en todos mayor inclinación en la zona del apoyo izquierdo pa-
los casos son: módulo de elasticidad E=30 GPa y co- ra salvar la abertura de la puerta. En la zona del apo-
eficiente de Poisson ν = 0,2. yo derecho el tirante sube con una inclinación más
suave para colaborar en el equilibrio vertical del
4.1. Ejemplo 1: Viga de gran canto con abertura y campo de compresiones que allí se desarrolla entre el
carga superior apoyo y la carga superior. Este esquema puede no
ser conveniente desde el punto de vista de la dispo-
El criterio de armado 1 aplicado en el primer caso de sición de armadura en obra. El objetivo del criterio 2
estudio converge tras 15 iteraciones. En la Figura 6 se es precisamente tener en cuenta estos aspectos para
muestra la distribución de los campos de tensiones potenciar esquemas de armados ortogonales, lo cual
de tracción y compresión junto a la distribución del se aprecia en la misma figura. El proceso de iteración
factor de eficacia que determina la forma estructural de la topología de BT según este criterio ha convergi-
optimizada. Puede aprecierse que la distribución de do en 9 iteraciones.
eficiencias propuestas tiende a concentrar las tensio-
nes en bandas de tracción y compresión; las resultan- La Figura 7 muestra la capacidad del criterio 2 de ge-
tes de dichas bandas representan los tirantes y las nerar esquemas con tirantes orientados de forma ver-
bielas respectivamente. tical y horizontal mientras que las bielas se desarrollan
en el continuo hormigón adquieren la orientación ne-
En la Figura 7 se indica el flujo de fuerzas internas cesaria para obtener un esquema BT factible. Se resalta
obtenido de acuerdo a los criterios de armado 1 (isó- que las posiciones y orientaciones de los tirantes son
tropo) y 2 (armadura ortogonal). El flujo de fuerzas un resultado del modelo propuesto, es decir, no es ne-
define con bastante precisión el esquema de bielas y cesario que el usuario disponga un esquema de arma-
tirantes mediante zonas de compresión y tracción y do, lo cual le da al método un valor adicional de cara
las respectivas inclinaciones de los campos en dichas al estudio de las posiciones óptimas de las armaduras
zonas. En la misma figura se ha representado una in- en regiones D singulares.
terpretación del esquema BT obtenido directamente
a partir del flujo de fuerzas tras la optimización. Las Figuras 8 y 9 muestran, para la aplicación de los
criterios 1 y 2 respectivamente, la evolución de la
En el esquema obtenido según el criterio 1 se desta- energía elástica total y la densidad de energía duran-
ca la presencia de tirantes inclinados en la cara infe- te los procesos de iteración, así como del volumen

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García
 Investigaciones y Estudios
Generación automática de esquemas de bielas y tirantes considerando criterios constructivos HORMIGÓN Y ACERO | 75

Figura 7. Flujo de fuerzas internas e interpretación como esquema BT para el Ejemplo 1 (adimensional)

Figura 8. Evolución del problema durante el proceso iterativo para el Ejemplo 1 – Criterio 1: a) doble de la energía total (mJ), b) doble
de la densidad de energía (mJ/mm3), c) volumen efectivo (mm3), d) doble de la densidad de energía de tracción (mJ/mm3)

Figura 9. Evolución del problema durante el proceso iterativo para el Ejemplo 1 – Criterio 2: a) doble de la energía total (mJ), b) doble
de la densidad de energía (mJ/mm3), c) volumen efectivo (mm3), d) doble de la densidad de energía de tracción (mJ/mm3)

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García
Investigaciones y Estudios 76 | HORMIGÓN Y ACERO Generación automática de esquemas de bielas y tirantes considerando criterios constructivos

Figura 10. Campos de tensiones principales (MPa) y distribución de factor de eficiencia (adimensional) para el Ejemplo 2

efectivo del material (Ec. 7) y de la energía absorbida A medida que se minimiza la energía de deforma-
por la tracción. Se destaca que en ambos criterios la ción, se observa que el volumen efectivo de material
energía total del sistema tiende a reducirse a medida también se reduce. Si bien los mínimos energéticos y
que avanza la iteración. Al final del proceso se obtie- de volumen efectivo no coinciden exactamente en la
ne, por lo tanto, una estructura optimizada en la que misma iteración, el volumen efectivo final es en ge-
se ha minimizado la energía y el volumen del mate- neral muy próximo al valor mínimo registrado du-
rial. No obstante, es evidente que ambos criterios rante toda la iteración. Como el volumen efectivo es
convergen a valores distintos, dadas las condiciones menor, la densidad de energía (energía por unidad
adicionales del criterio 2. de volumen) de toda la estructura crece, indicando
En general, se tendrá que la solución del criterio 2, si un mejor aprovechamiento del material efectivo.
bien produce el esquema de armado deseado, absor-
be más energía que la solución generada con el crite-
rio 1 al tener que satisfacer las condiciones construc- 4.2. Ejemplo 2: Viga de gran canto con abertura y
tivas adicionales. Por lo tanto, la elección del criterio carga colgada
2 de generación de esquemas BT implica asumir el
compromiso de obtener soluciones que no son ener- Las soluciones de los criterios 1 y 2 para el problema
géticamente las más eficientes pero que, además de de una viga de gran canto con carga colgada se re-
ser factibles, son constructivamente más adecuadas. presentan en las Figuras 10 para ambos criterios de

Figura 11. Flujo de fuerzas internas e interpretación como esquema BT para el Ejemplo 2 (adimensional)

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García
 Investigaciones y Estudios
Generación automática de esquemas de bielas y tirantes considerando criterios constructivos HORMIGÓN Y ACERO | 77

Figura 12. Campos de tensiones principales (MPa) y distribución de factor de eficiencia (adimensional) para el Ejemplo 3

armado. Asimismo, los flujos de fuerzas e interpreta- La evolución de las características energéticas y de
ción del esquema BT se representan en la Figura 11. volumen efectivo sigue tendencias similares a las co-
En el caso del criterio 1, el esquema BT resultante po- mentadas anteriormente, por lo que, en general, los
see tirantes inclinados en la cara inferior que parten esquemas BT propuestos tienden a minimizar la
del punto de aplicación de la carga. En el extremo iz- energía de deformación total a la vez que satisfacen
quierdo, los tirantes se concentran para pasar por las los condicionantes de armado.
esquinas inferior izquierda y superior derecha de la
abertura y terminar anclados en las bielas superiores.
En el extremo derecho, se forman dos tirantes que se 4.3. Ejemplo 3: Ménsula corta
prolongan hasta la biela de compresión para conse-
guir que la misma se incline hasta el apoyo. En la Figura 12 se comparan las distribuciones de ten-
siones principales y factores de eficacia obtenidos para
La aplicación del criterio de armado 2 produce un ti- el problema de la ménsula corta con relación de aspec-
rante vertical sobre la carga que se desarrolla hasta el to (L/h) de 0,43 según los criterios de armado 1 y 2. El
extremo superior de la viga, punto en el que se abren flujo de fuerzas internas y la interpretación de las mis-
las bielas compresión. En el lado derecho, la biela se mas como esquema BT se representa en la Figura 13
desarrolla recta hasta el apoyo, mientras que en el iz- para ambos criterios de armado.
quierdo, para salvar la abertura, la biela tiene dos
cambios de dirección en los puntos en los que se ge- La solución del criterio 1 consiste en disponer un ti-
neran tirantes verticales y horizontales bordeando la rante inclinado en el interior del pilar para desviar
abertura. las fuerzas horizontales producidas en la ménsula

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García
Investigaciones y Estudios 78 | HORMIGÓN Y ACERO Generación automática de esquemas de bielas y tirantes considerando criterios constructivos

Figura 13. Flujo de fuerzas internas e interpretación como esquema BT para el Ejemplo 3 (adimensional)

hasta las verticales dentro del pilar. Esta armadura Al resolver siempre problemas lineales, el método ha
inclinada produce dos bielas en el interior del nudo demostrado ser estable en cuanto a la convergencia;
que trabajan a un nivel de tensión relativamente ba- más aún, el número de iteraciones requerido es en
jo. Por otro lado, el criterio 2 genera dos bielas de general muy reducido, si se compara con otros méto-
compresión claramente delimitadas y rectas; las cua- dos de optimización topológica debido, por un lado,
les permiten que se generen solo dos tirantes en el a la propia eficacia de los criterios de toma de deci-
nudo (uno vertical y otro horizontal), ver Figura 13. sión y, por otro, a que en cada iteración se actúa so-
En este caso, la energía de ambas soluciones conver- bre todo el dominio mediante la distribución del fac-
gen a valores más parecidos que en los ejemplos an- tor de eficiencia.
teriores, siendo la energía de la solución 2 un 18%
mayor que la de la solución 1. El método no requiere una disposición inicial de los
tirantes o la solución de problemas no-lineales, por
lo que resulta de utilidad también para encontrar la
5. CONCLUSIONES mejor ubicación de las armaduras en regiones D sin-
gulares en las que no se tenga experiencia previa.
Se ha desarrollado un método numérico para la ge-
neración de esquemas de bielas y tirantes de regio- Las modificaciones de las propiedades elásticas se
nes D mediante un proceso iterativo de una serie de llevan a cabo mediante criterios racionales que utili-
problemas lineales-elásticos, que pueden resolverse zan unos indicadores basándose en la distribución
mediante métodos numéricos estándares. En el caso de la densidad de energía en el problema resuelto en
concreto de este trabajo, se desarrolló un programa el paso anterior. Así pues, el parámetro fundamental
de elementos finitos planos (tipo membrana) con ele- es el factor de eficiencia. Otros parámetros que se uti-
mentos triangulares de deformación constante. lizan para guiar la toma de decisiones son la densi-
dad de energía absorbida por las tensiones principa-
La estrategia fundamental del método es la modifi- les máxima y mínima por separado así como las ab-
cación de las propiedades elásticas del material de sorbidas por las tensiones en las direcciones vertical
acuerdo a criterios objetivos basados en la distribu- y horizontal. El esquema BT resultante tiende a pro-
ción de energía de deformación. Los criterios pro- ducir soluciones que reducen tanto la energía de de-
puestos tienden a mínima energía de deformación. formación como el volumen del material efectivo

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García
 Investigaciones y Estudios
Generación automática de esquemas de bielas y tirantes considerando criterios constructivos HORMIGÓN Y ACERO | 79

mientras que la densidad de energía de toda la es- [5] REINECK, K-H. (Ed.) Examples for the design of
tructura aumenta, lo cual indica que el volumen de structural concrete with strut-and-tie models. ACI SP-
material usado tiende a ser aprovechado cada vez 208. Michigan: American Concrete Institute, 2002.
más. Se observó que el mínimo de energía y de volu- 242 p. ISBN: 9780870310867.
men de material no se dan de forma simultánea, si
bien la diferencia entre el volumen para la mínima [6] REINECK, K-H.; NOVAK, L. (Eds.) Further exam-
energía y el volumen mínimo encontrado durante el ples for the design of structural concrete with strut-and-
proceso iterativo es muy pequeña. tie models. ACI SP-273. Michigan: American Concrete
Institute, 2010.
La principal ventaja de la estrategia propuesta está
en que es posible conseguir esquemas de armado [7] COMISIÓN PERMANENTE DEL HORMIGÓN.
que sean ortogonales para satisfacer así posibles con- Instrucción de Hormigón Estructural – EHE-08.
dicionantes constructivos si las propiedades mecáni- Madrid: Ministerio de Fomento, 2008.
cas del medio se modifican de forma adecuada. Para
ello se han desarrollado dos criterios de armado, uno [8] MUTTONI, A.; FERNÁNDEZ RUIZ, M. “Di-
que permite generar tirantes inclinados y otro que mensionamiento y verificación del hormigón estruc-
obliga a que estos se formen de forma ortogonal. El tural mediante el método de los campos de tensio-
criterio propuesto para producir armados ortogona- nes”. Hormigón y Acero. 2007, nº 243, p. 93-102.
les utiliza un material elástico ortótropo en los cam-
pos de tracciones o tirantes, de forma que los campos [9] MIGUEL SOSA, P.F.; FERNÁNDEZ PRADA,
de tracción tienen una dirección preestablecida. Para M.A.; BONET SENACH, J.L.; PASTOR VILLANUE-
discriminar entre las zonas que serán tratadas como VA, J.M. “Modelo automático para determinar el es-
tirantes frente a las que se considerarán bielas, el al- quema de bielas y tirantes de una región D de forma
goritmo utiliza el trabajo de tracción del elemento y cualquiera”. II Congreso ACHE de Puentes y Es-
compara cuando este ha sobrepasado un umbral de tructuras. Madrid: Asociación Científico-Técnica de
la densidad de energía absorbida en el punto. Esta Hormigón Estructural, 2002. p. 857-865. ISBN 84-
estrategia ha resultado ser satisfactoria como se evi- 89670-37-4.
dencia en los casos de estudio presentados.
[10] FERNÁNDEZ RUIZ, M.; MUTTONI, A. “On de-
Esta metodología abre la posibilidad de estudio de velopment of suitable stress fields for structural con-
soluciones de BT más eficientes así como la inclusión crete”. ACI Structural Journal, 2007, vol. 104, nº 4, p.
de otros tipos de condicionantes, tanto de tipo cons- 495-502.
tructivo como referentes al comportamiento en servi-
cio, funcionales, entre otros. [11] LIANG, Q.Q. “Performance-based optimization of
structures”. New York: Taylor & Francis Group, 2005.
(Spon Research) ISBN: 0-415-33594-9.
Agradecimientos
[12] KWAK, H.G.; NOH, S.H. “Determination of
Este trabajo ha sido realizado con el apoyo del strut-and-tie models using evolutionary structural
Ministerio de Ciencia e Innovación mediante del optimization”. Engineering Structures. 2006, vol. 28,
proyecto de investigación BIA- 2009-11764. nº 10, p. 1440-1449.

[13] BRUGGI, M. “Generating strut-and-tie patterns

REFERENCIAS for reinforced concrete structures using topology op-
timization”. Computers and Structures. 2009, vol. 87,
[1] SCHLAICH, J.; SCHÄFER, K.; JENNEWEIN, M. nº 23-24, p. 1483-1495.
“Towards a consistent design of structural concrete”.
[14] BIONDINI, F.; BONTEMPI, F.; MALERBA, P.G.
PCI Journal. 1987, vol. 32, nº 3, p. 74-150.
“Stress path adapting strut-and-tie models in
[2] ROMO, J. “Cálculo de secciones y elementos es- cracked and uncracked RC elements”. Structural
tructurales. El método de bielas y tirantes”. La EHE Engineering and Mechanics. 2001, vol. 12, nº 6, p. 685-
explicada por sus autores. 338 p. Madrid: Leynfor, 2000. 698.
ISBN: 8495560003.
[15] MAXWELL, B.; BREEN, J. “Experimental evalu-
[3] MIGUEL, P.; FERNÁNDEZ, M.; BONET, J.; ation of strut-and-tie model applied to deep beam
MARTÍ, K.; NAVARRO, J.; CASTRO, M.C. Proyecto de with opening”. ACI Structural Journal. 2000, vol. 97,
estructuras de hormigón mediante el método de bielas y ti- nº 1, p. 142-148.
rantes. Valencia: Ediciones VJ, 2009.
[16] BAIRÁN, J.M. “Nonlinear modelling of rein-
[4] ACHE. Comisión 1, GT1/3. Método de las bielas y forced concrete structures”. Earthquake Engineering by
tirantes. Madrid: ACHE, 2003. 164 p. (Monografía M- the Beach Workshop, July 2-4 2009, Capri, Italy.
6). ISBN: 84-89670-38-2. [Disponible en: www.reluis.it].

Volumen 63, nº 264, 67-79 | abril-junio 2012 | ISSN: 0439-5689 J.M. Bairán García