Anales de Edificación

Vol. 1, Nº2, 27-41 (2015)

Received: 20-07-2015 ISSN: 2444-1309
Accepted:10-08-2015 Doi: xxxxxxxx/xxxxxx

Diagramas momento-curvatura para secciones de

hormigón armado. Determinacion del valor de ductilidad
local de una sección de hormigón armado. Empleo del
programa matlab en aplicaciones estructurales a nivel
seccional (parte II).
Moment-curvature diagrams for reinforced concrete
sections. Sectional ductility ratio calculation. Applications
related to section study using matlab (part II).
Rodríguez Morales, Sergio

INITEC Industrial Técnicas Reunidas CO. (Spain, [email protected])

Resumen—En el entorno de la ingeniería de estructuras, resulta de interés el estudio de comportamiento a nivel seccional. Entre
los diversos procedimientos que se disponen a tal efecto, encontramos el “Método Multicapa”. El método introducido en un
artículo previo, será de nuevo tratado en este segundo texto para estudiar el concepto de ductilidad a nivel de sección, en elementos
de hormigón armado, y con la idea fundamental de mostrar las ventajas del método en relación al comportamiento No lineal de los
materiales.

Palabras clave— Ductilidad, curvatura; momento-curvatura; hormigón; matlab; acero; método multicapa; rotula plástica.

Abstract- Sectional study it is an important matter from structural engineering point of view. Among the different available
procedures in the present technical paper Multi-layer Method will be treated. The Method which was presented in the previous
article, once again will be used in this document to deal in this time with ductility concept in reinforced concrete sections. Also the
advantages of the method shall be exposed regarding to Nonlinear material behavior.

Index Terms— Ductility; curvature; momento-curvature; concrete; matlab; steel; multilayer method; plastic hinge.

 interés el estudio de comportamiento a nivel seccional. Entre

I. INTRODUCCIÓN los diversos procedimientos que se disponen a tal efecto,
encontramos el “Método Multicapa”. El método introducido
E N EL entorno de la ingeniería de estructuras, resulta de
en un artículo previo, será de nuevo tratado en este segundo
texto para estudiar el concepto de ductilidad a nivel de
Sergio Rodríguez Morales, At present finishing Master of Science in
Dynamics and Seismic Engineering at ETSIIM-UPM Master in Building
sección, en elementos de hormigón armado, y con la idea
Structures UPM-EE, Arquitecto Técnico UPM, and Associate Professor at fundamental de mostrar las ventajas del método en relación al
ETSEM-UPM during 2010 to 2012 years. Civil Engineering Leader at comportamiento No lineal de los materiales.
INITEC Industrial Técnicas Reunidas CO. (e-mail: [email protected],
[email protected]). Durante el proyecto sísmico de una estructura los métodos

Anales de Edificación, Vol. 1, Nº2, 27-41 (2015). ISSN: 2444-1309

28 Diagramas momento-curvatura para secciones de hormigón armado... (parte II).
Moment-curvature diagrams for reinforced concrete sections... (part II).

de análisis (como por ejemplo el método del empuje

incremental o “pushover” en sus diferentes posibilidades)
requieren la definición de diagramas momento-curvatura de
las secciones de los elementos que la componen. Entre la
información que aportan estos diagramas, encontramos el
concepto de ductilidad a nivel seccional. La ductilidad se
define como el cociente entre la curvatura última y la
curvatura elástica o asociada al comienzo de la plastificación
de la armadura de refuerzo.
En el presente artículo se pretende determinar los valores de
curvatura asociados al momento flector resultado de la
primera plastificación de la armadura y al momento último de
la sección, y por lo tanto el valor de la ductilidad a nivel
seccional. La obtención de estos parámetros se llevará a cabo Fig. 1. Figura 3.7 de Seismic Design Criteria Caltrans. Tanto si se
mediante un código programado en Matlab. quiere usar la curva real como la idealizada es necesario
También es objeto de este artículo describir someramente determinar en la curva real el punto donde se produce la
las relaciones entre la ductilidad a nivel de sección o local con plastificación de la armadura.
la ductilidad a nivel global, las cuales resultan de gran interés ni de su rigidez.
a la hora de enfrentarse al diseño de una estructura sometida a La ductilidad local o ductilidad a nivel de sección se
la acción sísmica. determina mediante un diagrama momento curvatura, siendo
esta el coeficiente entre el valor de curvatura última y el de
II. CONCEPTO DE DUCTILIDAD
curvatura en el momento de la primera plastificación de la
En general, la ductilidad se puede definir como la propiedad armadura:
de los materiales de experimentar deformaciones sin llegar a
(1)
romperse debida a la aplicación de una carga mantenida en el
tiempo. Desde un punto vista estructural y de acuerdo a
bibliografía especializada (Paulay & Priestly, Caltrans), un Un diagrama momento-curvatura se debe basar en el
elemento estructural dúctil es aquel que ha sido diseñado comportamiento lo más real posible de los materiales que
intencionadamente para admitir deformaciones inelásticas tras constituyen la sección. La Agencia para el Transporte del
muchos ciclos de carga debidos a la acción sísmica, sin que Estado de California (Caltrans) propone idealizar la curva con
esto implica una degradación significativa ni de su resistencia un diagrama elástico-plástico perfecto, de tal forma que la

Fig. 2. Trabajo en el rango inelástico de la estructura frente al tratamiento elástico.

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

S. Rodriguez 29

Fig. 3. Los dos primeros esquemas representan posibles mecanismos de colapso, mientras que el tercer mecanismo representa la mejor
situación posible para la formación de rótulas plásticas, en extremos de dinteles y en el arranque de los pilares en cimentación.

situación de la rama horizontal plástica quede definida tras obstante existen diversas ecuaciones que relacionan ambos
representar de forma conjunta la curva ideal y la curva real de parámetros. Se presentan a continuación algunas de estas
tal forma que se iguale el área tanto por encima como por relaciones de empleo en estructuras aporticadas de hormigón
debajo de la recta horizontal ideal. La rama elástica ideal será armado:
aquella que pase por el punto que marca la plastificación de la Paulay & Priestly
armadura de refuerzo en la curva real. En la figura 1 se
muestra el concepto descrito. (2)
Cabe también hablar de ductilidad desde un punto de vista
global, entendiendo esta como la habilidad de la estructura
donde de forma simplificada se puede adoptar la longitud de la
para trabajar en el rango inelástico permitiendo la disipación
rótula plástica como la mitad del canto de la sección:
de energía, mediante la formación de rótulas plásticas. La
ductilidad global se define normalmente en términos de
desplazamiento. La ductilidad global está estrechamente (3)
relacionado con el parámetro definido en las distintas
normativas (en la Normativa Americana se denomina Factor
de modificación de respuesta R, o el Factor de
Comportamiento q en el Eurocódigo) cuya función es reducir
el espectro de Respuesta Elástico. El concepto de trabajo en el
rango inelástico se describe en la figura 2.
En estructuras aporticadas la formación de rótulas plásticas
se debe localizar en los extremos de las vigas o en los
arranques de los pilares a nivel de cimentación
exclusivamente, ya que de lo contrario la estabilidad de la
estructura puede estar comprometida durante un evento
sísmico (ver figura 3 y figura 4). Este mecanismo de
disipación se garantiza cumpliendo ciertas relaciones
resistentes a nivel de flexión entre los elementos que
concurren a los nudos de la estructura (concepto de columna
fuerte-viga débil). Además, para permitir el trabajo de la
Fig. 4. Rotura de soporte de hormigón armado en edificio de
estructura en el rango inelástico es necesario garantizar la viviendas, durante el terremoto de Lorca en 2011. Situación
ductilidad a nivel global y local. Los pilares son elementos insegura a evitar y que pretende evitarse aplicando el concepto
caracterizados por bajos valores de ductilidad, ya que en estos de columna fuerte-viga débil.
los valores de curvatura son inferiores a los de las vigas,
debido a la presencia predominante de fuerzas axiles. La
relación entre ductilidad global y local es difícil de obtener, no Eurocódigo 8, Parte 1: Reglas Generales, Acciones
sísmicas y reglas para edificación. ENV 1998-1-1:1994: en su

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

30 Diagramas momento-curvatura para secciones de hormigón armado... (parte II).
Moment-curvature diagrams for reinforced concrete sections... (part II).

apartado 5.2.3.4 Condición de ductilidad Local encontramos reproducen a continuación, siendo necesario consultar el texto
en los comentarios la siguiente aproximación considerara original citado en las referencias para su estudio en
como conservadora: profundidad:

(4)

Si se emplea está última expresión se puede deducir que

para disponer de un valor de ductilidad global de la estructura
igual a 3, el menor valor de ductilidad que se debe aportar a
nivel seccional no puede ser inferior a 5.

III. ECUACIONES CONSTITUTIVAS DE LOS

MATERIALES
La ductilidad está estrechamente relacionado con la
ductilidad de los materiales que componen la estructura. En
bibliografía de referencia y en códigos internacionales de En el caso de una sección rectangular de hormigón armado
reconocido prestigio podemos encontrar una gran variedad de sin armadura de confinamiento, el diagrama tensión-
modelos de comportamiento de los materiales hormigón y deformación es el que aparece reflejado en la figura 5.
acero, descritos por medio de diagramas tensión-deformación. Mientras que en la figura 6 se muestran superpuestos los dos
A continuación se muestran algunos de ellos. modelos de hormigón confinado y sin confinar. Se puede
apreciar entonces el incremento en resistencia y en
deformación última que aporta la armadura de confinamiento
A. Diagrama tensión-deformación de acuerdo a EC-2 a la sección.
En el apartado 3.1.5. del mencionado código se encuentra la
ecuación constitutiva del hormigón para cálculo no-lineal:

En el modelo anterior no se tiene en cuenta la posibilidad de

confinar el hormigón mediante una armadura transversal
adecuado. Como es por todos sabido el confinamiento por
Fig. 5. Curva tensión-deformación según el modelo de Mander
medio de armadura transversal o camisas de acero no solo para un hormigón de resistencia 25 MPa sin confinamiento.
aporta resistencia si no que incrementa sustancialmente los
valores de deformación en rotura, lo que conlleva un aporte
adicional de ductilidad en la sección. El confinamiento se C. Modelo King para acero
puede implementar mediante el modelo de comportamiento De igual forma, se estable una curva tensión-deformación
que se describe a continuación. para el acero pasivo empleado en la sección. A tales efectos se
emplea el modelo de King (figura 7), que queda definida por
las siguientes expresiones:
B. Modelo Mander para Hormigón
En este modelo se describe el comportamiento tensión-
deformación para el hormigón armado bajo la hipótesis de
suponer este no confinado y confinado. Algunas de las
ecuaciones más representativas de dicho modelo se

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

S. Rodriguez 31

un plano de deformación en la sección se verificará el

equilibrio a nivel de fuerzas horizontales. El plano irá girando
con respecto a un punto fijo o charnela, generando fuerzas
resultantes en el acero y en el hormigón coherentes con el
plano de deformación impuesto. Una vez obtenido el
donde: equilibrio horizontal, se obtendrá el momento flector con
respecto al eje de la sección. El proceso se describe en detalle
en los puntos que se enumeran a continuación:

1) Es necesario definir el plano de deformaciones que actúa

en toda la altura de la sección, para ello es necesario
suponer que la sección permanece plana después de la
deformación y que por lo tanto la relación de
deformaciones es lineal. El sistema de referencia escogido
se sitúa en la fibra superior de la sección y la ecuación de
la recta que describe la deformación es la siguiente recta:

2) En el caso de la curvatura elástica es necesario fijar como

charnela del plano de deformaciones la posición de la
armadura traccionada. El valor de la deformación es este
punto es el correspondiente al de la plastificación de
acero. En lo referente a la curvatura asociada al momento
flector último, se escogerá como charnela el extremo
Fig. 6. Curva tensión-deformación del modelo confinado y el no comprimido de la sección, con el valor último
confinado superpuestos.
deformación el correspondiente a la rotura del hormigón.

3) El estado tensional de los materiales hormigón y acero

serán coherentes con el plano de deformación impuesto.
El caso del hormigón se genera una curva de resistencia
cuya extensión define la profundidad del bloque
comprimido. Las siguientes ecuaciones describen el
comportamiento en donde la deformación en función de la
profundidad y con respecto a la fibra superior.

Fig. 7. Curva tensión-deformación del acero según el modelo de

King.

IV. CÁLCULO DE LA DUCTILIDAD DE UNA

SECCIÓN DE HORMIGÓN ARMADO
En este caso el número de fibras consideradas en el cálculo
La metodología para obtener los valores de curvaturas es el número de puntos en los que se quiere evaluar la
requeridos consistirá en un proceso iterativo en la que fijado deformación y por lo tanto la resistencia del hormigón.

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

32 Diagramas momento-curvatura para secciones de hormigón armado... (parte II).
Moment-curvature diagrams for reinforced concrete sections... (part II).

Fig. 8. Posibles planos de deformación hasta conseguir el equilibrio horizontal durante la primera plastificación de la armadura de
refuerzo.

La resultante de este volumen de resistencia se obtiene

mediante la integración de esta curva entre 0 y la profundidad
del bloque comprimido en todo lo ancho de la sección.
En el caso del acero operamos de forma similar al
hormigón, con la excepción que la resultante de fuerzas de la
armadura se obtiene tan solo multiplicando la tensión del
acero por el área de refuerzo. Con respecto a la posición de las
fuerzas en el acero se suponen aplicadas en el CdG de las
armaduras.

El momento asociado a cada valor de curvatura se puede

calcular mediante la siguiente ecuación:

Fig. 9. Integración del volumen de resistencia a compresión

4) Las curvaturas o el giro que experimenta la sección plana
aportado por el bloque de hormigón.
queda definida por la siguiente relación:

De cara a determinar el momento de la sección asociado al

plano de deformaciones es necesario determinar el punto de
aplicación del bloque comprimido. Para ello será necesario
obtener el CdG de la ley de resistencia del hormigón, o lo que V. CASOS PRÁCTICOS DE APLICACIÓN
es lo mismo calcular el cociente entre el momento estático y el
En este apartado se obtendrán los valores de ductilidad de
área encerrada bajo la mencionada ley.
dos secciones de hormigón armado, una de ellas una viga

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

S. Rodriguez 33

plana y la otra una viga de canto, mediante los códigos

diseñados a tal efecto en el entorno Matlab. Ambas secciones
tendrán una resistencia a flexión similar. Los datos necesarios
para el cálculo de la ductilidad de ambas secciones se
presentan a continuación:

1) El hormigón tendrá una resistencia media a compresión

de 30 MPa, donde la máxima deformación tendrá un valor
de 0.2% y la deformación última un valor de 0.4%. La
armadura transversal no será capaz de proporcionar un
confinamiento adecuado al hormigón, por lo que se puede
emplear el modelo de comportamiento definido en el EC-
2.
2) Armadura pasiva tendrá un límite elástico de 420 MPa y
un valor de rotura 520 MPa. Se empleará como diagrama
de tensión-deformación el modelo de King, donde el valor
de deformación del final del escalón de cedencia del acero Fig. 11. Diagrama Tensión-deformación del hormigón.
tiene un valor 1.25% y la deformación en la rotura es de
10%.
3) La geometría de las secciones de hormigón armado tendrá
las siguientes dimensiones 300x600 (HxB) y 600x30
(HxB). La posición de la armadura longitudinal y su
cantidad están definidos en la figura 10.

Fig. 12. Modelo de comportamiento del material acero de

acuerdo a King.
Fig. 10. Definición geométrica y armaduras de las secciones
estudiadas.
C. Curvatura Última
VI. SOLUCIONES Con respecto a las curvaturas asociadas a los momentos
A. Ecuaciones constitutivas últimos se puede apreciar que la viga plana no es capaz de
alcanzar la máxima deformación en el hormigón, por lo que el
Los materiales acero y hormigón son iguales en ambas
valor de curvatura es inferior al de la viga de canto (Tabla II).
secciones, siendo las relaciones tensión-deformación de ambos
En la figura 13 se ha representado el plano de deformación
las representadas en las figuras 11 y 12.
de la viga plana, y en la figura 14 la curva que describe la
resistencia del hormigón del bloque comprimido.
B. Curvatura Elástica
De acuerdo al programa Matlab los resultados obtenidos
D. Ductilidad
para ambas secciones cuando se produce la primera
A la vista de los resultados los valores de ductilidad para
plastificación de la armadura se recogen en la tabla I.
ambas secciones son los siguientes:
Como era de esperar en la viga de mayor canto se produce
Viga Plana 300x600 (HXB)
antes la plastificación de la armadura que en la viga plana, lo
Viga de Canto 600x300 (HxB)
que implica un valor de curvatura menor.

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

34 Diagramas momento-curvatura para secciones de hormigón armado... (parte II).
Moment-curvature diagrams for reinforced concrete sections... (part II).

TABLA I
RESULTADOS DE AMBAS SECCIONES CUANDO SE PRODUCE LA PRIMERA PLASTIFICACIÓN DE LA ARMADURA

TABLA II
RESULTADOS DE AMBAS SECCIONES CUANDO SE PRODUCEN LOS MOMENTOS ÚÑTIMOS

VII. CONCLUSIONES 3) En el caso de un evento sísmico y siempre que se diseñe

En el presente artículo se ha vuelto demostrar la gran la estructura para trabajar en el rango inelástico, parece
utilidad del método multicapa en todo aquello relacionado con conveniente el empleo de secciones del mayor canto
el estudio de las secciones. Entornos como Matlab facilitan la posible y incorporar en el diseño de la estructura el
generación de cálculos iterativos característicos de este tipo de concepto de columna fuerte-viga débil.
problemas, además de facilitar la generación de salidas
gráficas que facilitan enormemente la compresión de los
resultados.
Desde el punto de vista estricto de la aplicación en
hormigón tratado en este segundo artículo se pueden extraer
las siguientes conclusiones:

1) El tratamiento no-lineal del comportamiento de los

materiales es factible, y aporta resultados más precisos
que las simplificaciones habituales tales como el
diagrama parábola rectángulo simplificado del hormigón
o el comportamiento elástico-plástico perfecto del acero.
Los valores de ductilidad a nivel de sección que están
estrechamente ligados a la ductilidad de los materiales
que la componen, son representados de una forma más
rigurosa.
Fig. 14. Curva de resistencia del hormigón en el bloque
2) La capacidad de rotación de la viga de canto es mayor que
comprimido, para la viga de canto estudiada.
la de la viga plana por lo que el valor de la ductilidad de
la primera es claramente mayor. Este aspecto queda
reflejado en la figura 15 donde se comparan los dos
diagramas momento-curvatura ideales de las secciones de
hormigón armado estudiadas.

Fig. 15. Diagramas Momento-Curvatura idealizados de las

secciones de hormigón armado estudiadas.
Fig. 13. Plano de deformaciones correspondiente al momento
último de la viga plana.

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

S. Rodriguez 35

APÉNDICE A %
figure;area(ecun,fcu(ecun),'FaceColor','c')
%%%%%%%%%%%%%%%%%%%%%%%%%%555 grid on
% legend(': Unconfined Concrete');grid on;set(gca,'Layer','top');
%Momento último cálculo de momento curvatura ylabel('Stress [MPa]','FontSize',16); xlabel('Strain','FontSize',16);
% title('Ecuación constitutiva Hormigón EC-2 sin confinar','FontSize',14)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% hold off
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%2.2 Material Acero
%%%%%%%%%%%%% %--------------------
% CALCULO DE VALOR DEFORMACIÓN DE
% FIBRA DE HORMIGÓN r = esu - esh;
% CUANDO LA ARMADURA HA COMENZADO m = ((fsu/fyk)*((30*r+1)^2)-60*r-1)/(15*(r^2));
% A PLASTIFICAR. ey = fyk/Es;
%
%Objeto: El programa pretende hallar el valor de deformación de la fibra fs=@(es)((m*(abs(es)-esh)+2)./(60*(abs(es)-esh)+2)+(abs(es)-esh)*(60-
%extrema de hormigón cuando la armadura ha alcanzado el límite elástico. m)./(2*((30*r+1)^2)))*(-1)*fyk.*(-esu<=es&es<-esh)+(-1)*fyk.*(-
%Se realizará un cálculo No-Lineal, donde apoyándose en las ecuaciones esh<=es&es<-ey)+Es*es.*(-
%constitutivas del hormigón (Mander) y del acero (King), se realizarán ey<=es&es<0)+Es*es.*(0<=es&es<ey)+fyk.*(ey<=es&es<esh)+((m*(es-
%varios tanteos hasta que la solución converja, es decir que consiga esh)+2)./(60*(es-esh)+2)+(es-esh)*(60-
%el equilibrio m)./(2*((30*r+1)^2)))*fyk.*(esh<=es&es<=esu)
% es1=-esu:0.0001:-esh; es2=-esh:0.0001:-ey; es3=-
clc ey:0.0001:0;es4=0:0.0001:ey; es5=ey:0.0001:esh; es6=esh:0.00001:esu;
clear figure;
close all hold on
format short g area(es1,fs(es1),'FaceColor','g')
%1. INTRODUCCIÓN DE DATOS hold on
area(es2,fs(es2),'FaceColor','g')
%1.1. Definición geométrica: hold on
B= 300 %Ancho de la sección [mm] area(es3,fs(es3),'FaceColor','g')
H=600 %Canto total de la sección [mm] hold on
rec=50 %Recubrimiento de la sección [mm] area(es4,fs(es4),'FaceColor','g')
%% hold on
area(es5,fs(es5),'FaceColor','g')
As=3*201 %Area de armadura pasiva TRACCIONADA [mm^2] hold on
Asc=3*201 %Area de amrmadura pasiva COMPRIMIDA [mm^2] area(es6,fs(es6),'FaceColor','g')
hold on
%1.2. Definición de los materiales title('Steel Stress Vs. Strain "King Model"','FontSize',16);grid
%ACERO on;set(gca,'Layer','top');
fyk=420 %Límite elástico del acero (MPa) ylabel('Stress [MPa]','FontSize',16); xlabel('Strain','FontSize',16);
fsu=520 %Límite de rotura del acero (MPa) hold off
Es=210000 %Modulo de elasticidad del acero (MPa) %%%%%%%%%%%%%%%%%%
esu=0.10 %Deformación última unitaria del acero
esh=0.0125 %Fin de la meseta de cedencia del acero %%%%%%%%%%%%%%%%%
dels=0.003 %Cálcula el número de puntos para evaluar el acero %%%%%%%%%%%%%%%%%%
%HORMIGÓN %3. Determinación de MOMENTO ÚLTIMO DE LA SECCIÓN
%Deformación del hormigón ec2=2*ec1 y variando el acero
fcm=30 %Resistencia del hormigón empleado en MPa ec1b=ec2 %Es aquí donde hay que meter el valor a tantear HORMIGÓN
%Sin confinar %esb=linspace(ey,esu,10) %Es aquí donde meter el valor deformación
ec1=0.002 %Deformación hormigón ACERO a tantear
ec2=2*ec1 %Deformación asociada a tensión hormigón 0 %Primero usar esta aproximaxión. Una vez localizado
%Deformación asociada a rotura hormigón sin confinar %el punto acotar con mayor precisión.
Ec=39000 %Modulo de deformación del hormigón
%hormigón confinado %esb=linspace(0.01963,0.02,5)
ecc1=0.0167 esb=0.043
%1.3. Numero de iteracciones d=H-rec
n=5 %Nº de puntos para buscar ec asociada a 1ª plastificación dim=length(esb)
nc=3 %Nº de puntos par el diagrama Momento-Curvatura
%2. GENERACIÓN DE LAS ECUACIONES CONSTITUTIVAS DE LOS Resultados_myk=zeros(dim,9);
MATERIALES Resultados_acero_comp_myk=zeros(dim,3);
%%2.1. Hormigón sin confinar for i=1:dim
%--------------------------- %
alpha=1. %
fcc=fcm*alpha %Plano de deformaciones
exilonb = @(y)ec1b-((ec1b+esb(i))/(H-rec)).*y;
ecun=0:0.0001:ec2; %Deformación hormigón sin confirnar figure
ecc=ec1*(1+5*(fcc/fcm-1)); fplot(exilonb,[0 H]);
k=1.1*Ec*ec1/fcm; title('Plano de deformaciones')
axis on
a=@(ecun)k*ecun/ec1-(ecun/ec1).^2; grid on
b=@(ecun)(1+(k-2)*ecun/ec1); %
fcu=@(ecun)fcm*a(ecun)./b(ecun); %Profundidad bloque comprimido en metros m

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

36 Diagramas momento-curvatura para secciones de hormigón armado... (parte II).
Moment-curvature diagrams for reinforced concrete sections... (part II).

bloqueb=ec1b*(H-rec)/(ec1b+esb(i)); Resultados_myk(i,9)=(Mcbc1+Ms)/1e6;
Resultados_myk(i,2)=bloqueb; %Curvatura
Resultados_myk(i,10)=ec1b/(Resultados_myk(i,2)/1e6)
%Curva de resistencia del material hormigón %

ecc=ec1*(1+5*(fcc/fcm-1)); % brazo=@(y)0.5*H-y;
k=1.1*Ec*ec1/fcm; %
y=0:0.001:bloqueb; Mcbc1b=quad(@(y)sigma_no_cob(y).*B.*brazo(y),0,Resultados_myk(i,2))+
a=@(y)k*exilonb(y)./ec1-(exilonb(y)./ec1).^2; C_acerob*(H/2-rec);
b=@(y)(1+(k-2)*exilonb(y)./ec1); % Mcbc1b=Mcbc1b/1e6;
fcu=@(y)fcm*a(y)./b(y); % Msb=Resultados_myk(i,6)*(0.5*H-rec)/1e6;
sigma_no_cob=@(y)fcu(y); %la resistencia en N/mm^2 % Resultados_myk(i,9)=Mcbc1b+Msb;
figure % Resultados_myk(i,10)=ec1b/(Resultados_myk(i,2)/1e6);
fplot(sigma_no_cob,[0 bloqueb]);
area(y,sigma_no_cob(y),'FaceColor','c'); end
grid on end
title('Stress Vs. C Sol. B"','FontSize',16);grid on;set(gca,'Layer','top');
ylabel('Stress Unc [Mpa]' ,'FontSize',16); xlabel('Compress. c disp('*****************OUTPUT-
[mm]','FontSize',16); RESULTADOS********************')
disp('******************************************************')
end disp('B)')
%figure disp('Se han pivotando los planos de deformación desde ec:_')
%fplot(sigma_no_co,[0 bloque]); ec1b
%axis on disp('Resultados Momento último max deformación hormigón y Def Acero
%grid on variando')
% disp(' Ycdg Horm - Prof.Bloque(mm)- Res.Comp(N)- Def. Acero
sigma,s(MPa) Trac. Arm(N) Check C/T -M(mKN) Curvatura')
for i=1:dim Resultados_myk
%Curva de resistencia del material acero disp('Si valor Check es "0" fuera intervalo 0.9<C/T<1.1 y Momento no
escb=ec1b-(ec1b+esb(i))/(H-rec).*rec; válido')
Resultados_acero_comp_myk(i,1)=escb; disp('Si valor Check es "1" Diferencia fuerzas horizontales pequeña dentro
0.9<C/T<1.1')
sigma_acero_compb=fs(escb); disp('******************************************************')
Resultados_acero_comp_myk(i,2)=sigma_acero_compb; disp('Información relativa a armadura comprimida')
C_acerob=sigma_acero_compb*Asc; disp(' Def As` - Stress Acero(MPa) - Res.Comp(N)')
Resultados_acero_comp_myk(i,3)=C_acerob; Resultados_acero_comp_myk

%C_acerob, compresion armadura

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Resultante del bloque comprimido en N y en mm
%%%%%%%%%%%%%
horm_compb=(quad(@(y)sigma_no_cob(y).*B,0,bloqueb));
% CALCULO DE VALOR DEFORMACIÓN DE
Resultados_myk(i,3)=horm_compb+C_acerob;
% FIBRA DE HORMIGÓN
%Chequeo con armadura
% CUANDO LA ARMADURA HA COMENZADO
%Capacidad de armadura en KN
% A PLASTIFICAR.
Resultados_myk(i,4)=esb(i);
%
Resultados_myk(i,5)=fs(esb(i));
%Objeto: El programa pretende hallar el valor de deformación de la fibra
Tb=As*fs(esb(i));
%extrema de hormigón cuando la armadura ha alcanzado el límite elástico.
Resultados_myk(i,6)=Tb;
%Se realizará un cálculo No-Lineal, donde apoyándose en las ecuaciones
%%%%%%
%constitutivas del hormigón (Mander) y del acero (King), se realizarán
%varios tanteos hasta que la solución converja, es decir que consiga
if (0.90<Resultados_myk(i,3)/Tb&Resultados_myk(i,3)/Tb<1.10)
%el equilibrio
Resultados_myk(i,7)=1;
%
Resultados_myk(i,8)=Resultados_myk(i,3)/Tb;
clc
else
clear
Resultados_myk(i,7)=0;
close all
Resultados_myk(i,8)=Resultados_myk(i,3)/Tb;
format shortg
end
%1. INTRODUCCIÓN DE DATOS
%%%%%
%1.1. Definición geométrica:
if Resultados_myk(i,7)==0;
B= 300 %Ancho de la sección [mm]
Resultados_myk(i,9)=0;
H=600 %Canto total de la sección [mm]
else
rec=50 %Recubrimiento de la sección [mm]
%
%%
%Conocida la resultante se obtiene el CDG de la distribución de
%presiones del hormigón
As=3*201 %Area de armadura pasiva TRACCIONADA [mm^2]
A=quad(@(y)sigma_no_cob(y).*y,0,bloqueb)
Asc=3*201 %Area de amrmadura pasiva COMPRIMIDA [mm^2]
B=quad(@(y)sigma_no_cob(y),0,bloqueb)
Ycdg=A/B
%1.2. Definición de los materiales
Resultados_myk(i,1)=Ycdg
%ACERO
%Momento por encima de la mitad de la sección
fyk=420 %Límite elástico del acero (MPa)
Mcbc1=((H/2)-Ycdg)*(Resultados_myk(i,3)-C_acerob)+(0.5*H-
fsu=520 %Límite de rotura del acero (MPa)
rec)*C_acerob
Es=210000 %Modulo de elasticidad del acero (MPa)
%Momento por debajo de la mitad de la sección
esu=0.10 %Deformación última unitaria del acero
Ms=Resultados_myk(i,6)*(0.5*H-rec);
esh=0.0125 %Fin de la meseta de cedencia del acero
%Momento último

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

S. Rodriguez 37

dels=0.003 %Cálcula el número de puntos para evaluar el acero

%%3.1 Planos de deformación
%HORMIGÓN %--------------------

fcm=30 %Resistencia del hormigón empleado en MPa %ec=0.00001:0.0001:ec1

%Sin confinar %Se recomienda ir tanteando en este intervalo las posibles deformaciones
ec1=0.002 %Deformación hormigón %ec=linspace(0.00077,0.00078,n)
ec2=2*ec1 %Deformación asociada a tensión hormigón 0 %Una vez cazada la deformación tantear con un solo valor
%Deformación asociada a rotura hormigón sin confinar ec=0.00045
Ec=39000 %Modulo de deformación del hormigón dim=length(ec)
%hormigón confinado
ecc1=0.0167 eyk=fyk/Es
d=H-rec
%1.3. Numero de iteracciones
n=10 %Nº de puntos para buscar ec asociada a 1ª plastificación Resultados=zeros(dim,9)
nc=3 %Nº de puntos par el diagrama Momento-Curvatura Resultados_acero_comp=zeros(dim,3)
%2. GENERACIÓN DE LAS ECUACIONES CONSTITUTIVAS DE LOS
MATERIALES for i=1:dim

%%2.1. Hormigón sin confinar %

%--------------------------- ec(i);
alpha=1. Resultados(i,1)=ec(i);
fcc=fcm*alpha %
%Plano de deformaciones
ecun=0:0.0001:ec2; %Deformación hormigón sin confirnar exilon = @(y)ec(i)-(ec(i)+eyk)/(H-rec).*y;
ecc=ec1*(1+5*(fcc/fcm-1)); %fplot(exilon,[0 H]);
k=1.1*Ec*ec1/fcm; %axis on
%grid on
a=@(ecun)k*ecun/ec1-(ecun/ec1).^2; %
b=@(ecun)(1+(k-2)*ecun/ec1); %Profundidad bloque comprimido en metros m
fcu=@(ecun)fcm*a(ecun)./b(ecun); bloque=ec(i)/(ec(i)+eyk)*d;
% Resultados(i,2)=bloque;
figure;area(ecun,fcu(ecun),'FaceColor','c')
grid on %Curva de resistencia del material hormigón
legend(': Unconfined Concrete');grid on;set(gca,'Layer','top'); %Resistencia a N/mm^2
ylabel('Stress [MPa]','FontSize',16); xlabel('Strain','FontSize',16);
title('Ecuación constitutiva de Mander Hormigón sin confinar','FontSize',14) ecc=ec1*(1+5*(fcc/fcm-1));
hold off k=1.1*Ec*ec1/fcm;
y=0:0.001:bloque;
%%2.2 Material Acero a=@(y)k*exilon(y)./ec1-(exilon(y)./ec1).^2;
%-------------------- b=@(y)(1+(k-2)*exilon(y)./ec1);
fcu=@(y)fcm*a(y)./b(y);
r = esu - esh; sigma_no_co=@(y)fcu(y); %la resistencia en N/mm^2
m = ((fsu/fyk)*((30*r+1)^2)-60*r-1)/(15*(r^2)); figure
ey = fyk/Es; fplot(sigma_no_co,[0 bloque]);
area(y,sigma_no_co(y),'FaceColor','c');
fs=@(es)((m*(abs(es)-esh)+2)./(60*(abs(es)-esh)+2)+(abs(es)-esh)*(60- grid on
m)./(2*((30*r+1)^2)))*(-1)*fyk.*(-esu<=es&es<-esh)+(-1)*fyk.*(- title('Stress Vs. C Sol. B"','FontSize',16);grid on;set(gca,'Layer','top');
esh<=es&es<-ey)+Es*es.*(- ylabel('Stress Unc [Mpa]' ,'FontSize',16); xlabel('Compress. c
ey<=es&es<0)+Es*es.*(0<=es&es<ey)+fyk.*(ey<=es&es<esh)+((m*(es- [mm]','FontSize',16);
esh)+2)./(60*(es-esh)+2)+(es-esh)*(60-
m)./(2*((30*r+1)^2)))*fyk.*(esh<=es&es<=esu)
es1=-esu:0.0001:-esh; es2=-esh:0.0001:-ey; es3=-
ey:0.0001:0;es4=0:0.0001:ey; es5=ey:0.0001:esh; es6=esh:0.00001:esu; %figure
figure; %fplot(sigma_no_co,[0 bloque]);
hold on %axis on
area(es1,fs(es1),'FaceColor','g') %grid on
hold on %
area(es2,fs(es2),'FaceColor','g') %Curva de resistencia del material acero
hold on es=ec(i)-(ec(i)+eyk)/(H-rec).*rec;
area(es3,fs(es3),'FaceColor','g') Resultados_acero_comp(i,1)=es;
hold on
area(es4,fs(es4),'FaceColor','g')
hold on sigma_acero_comp=fs(es);
area(es5,fs(es5),'FaceColor','g') Resultados_acero_comp(i,2)=sigma_acero_comp;
hold on C_acero=sigma_acero_comp*Asc;
area(es6,fs(es6),'FaceColor','g') Resultados_acero_comp(i,3)=C_acero;
hold on
title('Steel Stress Vs. Strain "King Model"','FontSize',16);grid %Resultante del bloque comprimido en N y en m
on;set(gca,'Layer','top'); Res_horm(i)=(quad(@(y)sigma_no_co(y).*B,0,bloque))
ylabel('Stress [MPa]','FontSize',16); xlabel('Strain','FontSize',16); Resultados(i,3)=Res_horm(i)+C_acero;
hold off %Chequeo con armadura
%%%%%%%%%%%%%%%%%% %Capacidad de armadura en KN
T=As*fyk
%%3. PLANOS DE DEFORMACIÓN Y EQUILIBRIO DE LA SECCIÓN Resultados(i,4)=T;

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

38 Diagramas momento-curvatura para secciones de hormigón armado... (parte II).
Moment-curvature diagrams for reinforced concrete sections... (part II).

if (0.90<Resultados(i,3)/T&Resultados(i,3)/T<1.10) structure. One of the possible information which can be

Resultados(i,5)=1;
Resultados(i,6)=Resultados(i,3)/T; obtained by a moment-curvature diagram is the ductility ratio,
else which is defined as the relationship between the ultimate
Resultados(i,5)=0;
curvature and elastic curvature. The elastic curvature is
Resultados(i,6)=Resultados(i,3)/T; associated with the first yielding of tensile steel reinforcement.
end In this technical paper curvature related to ultimate and
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
elastic bending moment shall be calculated using a code
programed in Matlab. Once curvatures shall be calculated,
if Resultados(i,5)==0;
Resultados(i,7)=0; ductility ratio for two sections shall de estimated.
else The aim of this technical paper is also described briefly the
%Conocida la resultante se obtiene el CDG de la distribución de
%presiones del hormigón relationship between local and global ductility, which are of
A=quad(@(y)sigma_no_co(y).*y,0,bloque) great significance when a structure must be designed subjected
B=quad(@(y)sigma_no_co(y),0,bloque)
to seismic loading.
Ycdg=A/B
Resultados(i,9)=Ycdg
%Momento por encima de la mitad de la sección
Mcbc1=(0.5*H-Ycdg)*Res_horm
Mcbc2=(0.5*H-rec)*C_acero II. BRIEF DUCTILITY DESCRIPTION
%Momento por debajo de la mitad de la sección Ductility shall be defined as material property to suffer
Ms=Resultados(i,4)*(0.5*H-rec);
%Momento elastico deformations without reaching the rupture, due to a sustaining
Resultados(i,7)=(Mcbc1+Mcbc2+Ms)/1e6; load. From a structural point of view and in accordance with
%Curvatura
Resultados(i,8)=ec(i)/(Resultados(i,2)/1e6) technical reference and international codes (Paulay & Priestly,
end Caltrans), ductility shall be defined as the structure ability to
end
undergo inelastic deformations without losing neither strength
%%%%%%%%%%%%%%% capacity nor stiffness. Ductility is also related to cyclic
clc loading due to earthquake events.
disp('*****************OUTPUT-
RESULTADOS********************') Local ductility ratio is provided by a moment-curvature
disp('******************************************************') analysis, and is the result of the relationship between curvature
disp('A)')
disp('Resultados finales al alcanzar el acero a tracción fyk') associated to ultimate bending moment and curvature related
disp(' Def ec - Prof.Bloque(mm)- Res.Comp(N)- Trac. Arm(N)- Check - to elastic bending moment. The following equation described
C/T -M(mKN) Curvatura (1/Km) Ycgd(mm)')
Resultados the concept:
disp('Si valor Check es "0" fuera intervalo 0.9<C/T<1.1 y Momento no
válido')
disp('Si valor Check es "1" Diferencia fuerzas horizontales pequeña dentro
0.9<C/T<1.1') Moment-Curvature shall be based on the most real behavior
disp('******************************************************')
of the structural materials. As example Caltrans proposed to
disp('Información relativa a armadura comprimida')
disp(' Def As` - Stress Acero(MPa) - Res.Comp(N)') create an idealized diagram instead a real moment-curvature
Resultados_acero_comp diagram. In the ideal diagram horizontal branch shall be
located in the proper place where areas in the diagrams are
balanced. Also elastic branch curve should pass through the
APÉNDICE B (VERSIÓN EN INGLÉS)
point marking the first reinforcing bar yield. See figure 1 to
clarify the described concept.
I. INTRODUCTION
Ductility can be also defined including the whole structure,
Sectional study it is an important matter from structural or the possibility to work in inelastic range in which is
engineering point of view. Among the different available possible to dissipate energy through plastic hinge formation.
procedures in the present technical paper Multi-layer Method Global ductility usually is expressed in displacement terms.
will be treated. The Method which was presented in the Global displacement is directly related to code parameters as
previous article, once again will be used in this document to Modification Factor Response o behavior factor, cited in
deal in this time with ductility concept in reinforced concrete American and European codes respectively. Inelastic and
sections. Also the advantages of the method shall be exposed elastic range structural working is represented graphically in
regarding to No-linear material behavior. figure 2.
When a structure subjected to earthquake loading several In RC frames hinge formation shall be located exclusively
analysis methods for calculation are available, and some of in support point beams and in the lowest portion of the
them as time-history or pushover analysis require a section columns right above foundation, because otherwise stability
moment-curvature diagram of the elements forming the structure could not be guaranteed (Figure 3). Energy

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

S. Rodriguez 39

dissipation mechanism is able fulfilling bending capacity relationship. Down below the most important equations of the
strength conditions between columns and beams reaching the model are specified, and reader can read further explanations
same joint (strong column-weak beam concept). Besides, about the model in the article pointed in references.
structural working in inelastic range is only possible providing
certain ductility capacity to the structural elements conforming
the structure. Regarding to columns these are defined as low
ductility elements, since curvature values are rather little than
curvature values in beams due mainly to axial forces.
Relationship between local and global ductility are difficult to
precise, nevertheless the following equations are usually
common is RC Frame structures:
Paulay & Priestly

For a regular RC section without special transversal

In a simple way Plastic hinge length shall be adopted as reinforcement stress stain relationships is represented in figure
follows: 5. In figure 6 confined and unconfined concrete models are
reflected together. Compressive strength capacity as strain
Eurocode 8, Part 1: General Rules, Seismic Actions and increasing are achieved due to a convenient confined stress
rules for buildings. ENV 1998-1-1:1994, iIn section 5.2.3.4, in provided by appropriate transversal reinforcement.
the commentaries is possible to read the following King Model for Steel reinforcement: Nonlinear stress strain
conservatively equation: diagram is defined by the following equations:

According to the former equation to reach a ductility value

of 3, it is necessary to provide structural elements which local
ductility ratio equal to 5.

III. NONLINEAR MATERIAL MODELS FOR

DUCTILE MEMBERS where:
Ductility is directly related to material behavior. In
technical references and international design codes are
specified a great number of different concrete and steel
reinforcing material models, all of them defined by their Stress
Strain diagrams.
IV. CALCULATION PROCEDURE FOR OBTAINING
Nonlinear concrete model in accordance with Eurocode 2:
DUCTILITY RATIO
In section 3.1.5., in European concrete code is defined the
concrete material model to be employed in nonlinear analysis:
Iterative process shall be carried out in order to obtain
curvature and associated bending moment. For this purpose is
necessary to fix a strain point located in tensile reinforcement
or in top concrete face, and then rotate deformation plane up
to achieve horizontal forces balance. Once horizontal balance
is achieved is able to get bending moment due to established
stress strain plane. The procedure is explained in detail below:
1) Section will stay plane after deformation and this
In this model is not possible the concrete to be confined due assumption allows to associate strain and stress in a linear
to an adequate transversal reinforcement or a steel casing. Is way. Reference point is concrete compressed surface and
well known that a confined concrete increase its compressive the equation which described section strain in each point
strength capacity as well as ultimate strain which provides is the following one:
extra local ductility ratio to the section. Confined concrete can
be modelled according to the following criteria.
Mander´s Model for concrete material: Unconfined and
confined concrete material can be treated with this strain stress

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

40 Diagramas momento-curvatura para secciones de hormigón armado... (parte II).
Moment-curvature diagrams for reinforced concrete sections... (part II).

2) In order to obtain elastic curvature is necessary to fix V. EXAMPLES

stress strain diagram at reinforcement tension level, Ductility ratio shall be estimated for two kinds of sections,
assuming in that point a strain corresponding to steel yield one a flat section and the second a depth section. Curvature
limit. Regarding to ultimate curvature fixed point shall be values shall be estimated using a tool programed in Matlab.
top of concrete under compression and deformation shall The information necessary to solve the problem is presented
be breaking strain of concrete. below:
3) Unconfined concrete strain is function of the depth of the 1) Compression strength for concrete is 30 MPa, maximum
section. The following equations describe this aspect: deformation in concrete shall be 0.2% and ultimate
deformation shall be 0.4%. Transversal reinforcement
shall not provide enough confine stress, so Eurocode
concrete material mode can be adopted in calculations.
2) Steel reinforcement at yield limit shall be 420 MPa, and
520 MPa shall be adopted as ultimate value. King´s model
shall be used as material behavior model, assuming an
ultimate deformations equal to 10% and latest horizontal
branch deformation equal to 1.25%.
3) Geometrical and reinforcement definition for studied
Number of layers is in accordance with the number of sections are defined in figure 10.
points where it is required the deformation in the concrete
material to be evaluated. VI. SOLUTIONS
By means of integrating concrete force resultant can be
A. Stress Strain relationships
achieved. Integration shall be performed taking into account
the area enclosed below concrete resistance profile between 0 RC and Steel reinforcement material stress strain
and the compressed depth. relationships are shown in figures 11 and 12.
B. Elastic Curvature
According to matlab the result for the elastic stage can be
consulted in the table I. As it was expected steel reinforcing
yield limit is reached by depth beam before flat beam, so
curvature value it is lower in depth beam.
It is also necessary to establish where concrete resultant force
C. Ultimate Curvature
is located, for that matter usual center of gravity procedure
Regarding the flat beam is not possible to reach the
shall be carried out:
maximum deformation (0.4%) for concrete, therefore
curvature value is clearly lower in the flat beam (Table II).
As example of output results, strength compression profile
in concrete and section deformation along the depth of the
beam are represented in figures 13 y 14.
For Steel reinforcement once strain is known, resultant
D. Ductility Ratio
forces provided by top and bottom rebar are easy to obtain,
According to previous results finally ductility ratios are
just multiplying steel area per stress.
presented below:
Flat beam 300x600 (HXB)
Depth beam 600x300 (HxB)

Bending moment at each stage shall be calculated according

to this equation: VII. CONCLUSIONS
Once again it was proved the utility of the multilayer
Finally curvature could be calculated dividing concrete top method regarding to structural section studies. Iterative
deformations by the length of concrete subjected under numerical methods typical in this kind of problems are
compression called “bloque” in the following equation: available to be performed by programs as matlab, providing
also graphical output to ease results understanding.

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309

S. Rodriguez 41

Regarding to specific concrete problem performed in this REFERENCES

technical article the following conclusions can be drawn: Bentz, E.; Collins, M.P. Response 2000, Sectional analysis
1) Nonlinear material behavior provides a more accurate programe developed at University of Toronto.
results than typical simplifications adopted in regular
Bisch, P.; Carvalho, E.; Degee, H. Eurocode 8: Seismic
calculations as perfect elastic plastic behavior in steel or Design of Buildings Worked Examples.
rectangular compression block in concrete. Ductility ratio
Caltrans Seismic Design Criteria, Version 1.7. Abril de 2013.
are clearly affected by these simplifications.
2) Rotation capacity in depth beams are clearly bigger than EN 1998-1: Eurocode 8: Design of structures for Earthquake
curvature in flat beam. That aspect is reflected in figure Resistance.
number 15, where the moment-curvature diagrams for the EN 1992-1-1: Design of Concrete Structures. General Rules
two sections under study are reflected together. In depth and Rules for Buildings.
beams ductility ratios is clearly bigger than flat beams. Mander, Priestly and Park. Theorical Stress-Strain Model for
3) In case a structure is needed to work under inelastic range Confined Concrete. Journal of Structural Engineering,
during an earthquake, is recommend to design the Vol. 114, No. 8, August 1988.
structure using beams with the biggest depth as possible
Paulay, T. Seismic Design of Reinforced Concrete and
due to their important ductility ratios and apply Strong Masonry Buildings. M.J.N. Priestley.
column-weak beam criteria.

Anales de Edificación, Vol. 1, Nº 2, 27-41 (2015). ISSN: 2444-1309