OBTENCIÓN DE LA MATRIZ DE RIGIDEZ DE UN ELEMENTO DE SECCIÓN VARIABLE.

Partiendo de la energía de deformación de un elemento plano se

obtiene la relación entre fuerzas y desplazamientos de un nodo extremo del
elemento, a través de la matriz de flexibilidades. Por otro lado, se deduce la
matriz de rigidez del nodo mencionado, invirtiendo simplemente la matriz de
flexibilidades. Después, aplicando desplazamientos unitarios y por equilibrio de
fuerzas se deduce la matriz de rigidez para ambos extremos del elemento de
sección variable. Finalmente, como un ejemplo de este trabajo se obtiene la
matriz de rigidez de un elemento de sección variable rectangular llena.

Matriz de rigidez de un elemento de sección variable.

Para la obtención de la matriz de rigidez de un elemento de sección

variable con fuerzas o desplazamientos aplicados en los nodos extremos (figura ,
se puede proceder como a continuación se describe.

La energía de deformación para un elemento plano con

comportamiento lineal se puede expresar como:

L ⎡ N 2 k1V y2 M z2 ⎤
W =∫ ⎢ + + ⎥dx
⎣⎢ 2 EA 2GA 2 EI z ⎦⎥
0

F2
F1 A B N x
F3 M

z V
x
Por equilibrio se tiene:
N (x ) = F1
V ( x ) = F2 (6-2)
M z ( x ) = F3 − F2 x

Sustituyendo los valores de N , V y M en la ecuación 6-1 se obtiene:

L ⎡ F12 (F3 − F2 x )2 k1 F22 ⎤

W =∫ ⎢ + + ⎥ dx
⎣⎢ 2 EA( x ) 2 EI z (x ) 2GA( x ) ⎦⎥
0
(6-3)
L ⎡ F12 F32 − 2 F2 F3 x + F22 x 2 k1 F22 ⎤
W =∫ ⎢ + + ⎥ dx
⎣ 2 EA( x ) 2 EI z ( x ) 2GA( x ) ⎦
0

De acuerdo al teorema de Castigliano:

∂W
= δi (6-4)
∂Fi

∂W F1 ⎡ L dx ⎤
= δ1 = ∫
L
∴ dx = ⎢ ∫ ⎥ F1
∂F1 0 EA( x )
⎣ 0 EA( x) ⎦

∂W L ⎡ 2F x − 2F x
2
2k1 F2 ⎤ ⎧ L ⎡ x2 k1 ⎤ ⎫ ⎡ L xdx ⎤
= δ2 = ∫ ⎢ 2 3
+ ⎥ dx = ⎨ ∫ ⎢ + ⎥ dx ⎬ F − ⎢ ∫0 EI ( x ) ⎥ F3
∂F2 0
⎣ 2 EI z ( x ) 2GA ( x ) ⎦ ⎩ ⎣ z
0 EI ( x ) GA( x ) ⎦ ⎭
2
⎣ z ⎦

∂W L ⎡ − 2F x + 2F ⎤ ⎡ L xdx ⎤ ⎡ L dx ⎤
= δ3 = ∫ ⎢ 2 3
⎥ dx = − ⎢ ∫ ⎥ F2 + ⎢ ∫ ⎥ F3 (6-5)
∂F3 0
⎣ 2 EI z ( x ) ⎦ ⎣ 0 EI z ( x )⎦ ⎣ 0 EI z ( x )⎦

Expresando estas relaciones en forma matricial se tiene:

⎡ L dx ⎤
⎢ ∫0 EA( x ) 0 0 ⎥
⎧δ 1 ⎫ ⎢ ⎥ ⎧ F1 ⎫
⎪ ⎪ ⎢ L⎛ x2 k1 ⎞ L xdx ⎥ ⎪ ⎪
⎨δ 2 ⎬ = ⎢ 0 ∫0 ⎜⎜⎝ EI z (x ) + GA(x ) ⎟⎟⎠dx − ∫0 EI z (x )⎥ ⎨F2 ⎬
⎪δ ⎪ ⎢ ⎥⎪F ⎪
⎩ 3⎭ L xdx dx ⎥ ⎩ 3 ⎭
⎢ −∫
L

⎢⎣
0
0 EI ( x )
z
∫0 EI z (x ) ⎥⎦

⎧δ 1 ⎫ ⎧ F1 ⎫ ⎡ f 11 0 0 ⎤ ⎧ F1 ⎫
⎪ ⎪ ⎪ ⎪ ⎢ ⎪ ⎪
⎨δ 2 ⎬ = [ f ]AA ⎨ F2 ⎬ = ⎢ 0 f 22 − f 23 ⎥⎥ ⎨ F2 ⎬ (6-6)
⎪δ ⎪ ⎪F ⎪ ⎢ 0 − f ⎪
f 33 ⎥⎦ ⎩ F3 ⎭⎪
⎩ 3⎭ ⎩ 3⎭ ⎣ 32

De acuerdo a lo anteriormente expuesto (ecuación 6-6), la matriz de

flexibilidades del nodo A es:
 ⎡ L dx ⎤
⎢ ∫0 0 0 ⎥
⎡ f11 0 0 ⎤ ⎢ EA( x ) ⎥
L⎛ k1 ⎞
− f 23 ⎥⎥ = ⎢⎢
x2 L xdx ⎥
[ f ]AA = ⎢⎢ 0 f 22 0 ∫0 ⎜⎜⎝ EI z (x ) + GA(x ) ⎟⎟⎠dx − ∫0 EI z (x )⎥
⎢⎣ 0 − f 32 f 33 ⎥⎦ ⎢ L xdx
⎥
⎢ L dx ⎥
0 −∫ ∫0 EI z (x ) ⎥⎦
⎢⎣ 0 EI ( x )
z

(6-7)

Dada esta matriz de flexibilidades del nodo A , (ecuación 6-7), y la

relación que existe entre las fuerzas y desplazamientos de un mismo nodo
extremo de un elemento (figura 6-2), se puede proceder como sigue, para la
obtención de la matriz de rigidez completa de un elemento de sección
variable.

F2 , δ2 F5 , δ5
F1 , δ1 F6 , δ6
F3 , δ3 x
F4 , δ4
A B
L

Figura 6-2.Fuerzas y desplazamientos en los nodos extremos del elemento.

A partir de la ecuación 6-7, se obtiene para el nodo A la relación

siguiente:

⎡ 1 ⎤
⎢f 0 0 ⎥
⎢ 11 ⎥
[FA ] = [k A ]{δ A } = [ f ]−A1 {δ }A =⎢ 0
⎢
f 33 f 23 ⎥
{δ }A (6-8)
D D⎥
⎢ f 23 f 22 ⎥
⎢ 0 ⎥
⎣ D D⎦

Donde [k AA ]
es la matriz de rigidez del nodo A y el determinante es:
D = Det = f 22 f 33 − f 232 .

a) Dado un desplazamiento unitario en A en la dirección de F1 (figura 6-3),

se pueden obtener las fuerzas en el nodo A y por equilibrio deducir las
fuerzas en el nodo B:
 ⎧ 1 ⎫
⎧ F11 ⎫ ⎡k11 0 0 ⎤ ⎧1⎫ ⎧k11 ⎫ ⎪ f ⎪
⎪ ⎪ ⎢ ⎪ ⎪ ⎪ ⎪ ⎪ 11 ⎪
⎨ F21 ⎬ = ⎢ 0 k 22 k 23 ⎥⎥ ⎨0⎬ = ⎨ 0 ⎬ = ⎨ 0 ⎬ (6-9)
⎪F ⎪ ⎢ 0 k 33 ⎥⎦ ⎪⎩0⎪⎭ ⎪⎩ 0 ⎪⎭ ⎪ 0 ⎪
⎩ 31 ⎭ ⎣ k 32
⎪ ⎪
⎩ ⎭

1
F51 = F61 = 0, F41 = − F11 = − (6-10)
f 11

F21 F51
F11 F61
F31 x
F41
A B
L
δ1= 1

Figura 6-3.Fuerzas provoc adas por un desplazamiento unitario en la direc c ión F1.

b) De la misma manera, se aplica un desplazamiento unitario en A en la

dirección de F2 (figura 6-4), para obtener las fuerzas en el nodo A y por
equilibrio deducir las fuerzas en el nodo B.

⎧ ⎫
⎧ F12 ⎫ ⎡k11 0 0 ⎤ ⎧0⎫ ⎧ 0 ⎫ ⎪ 0 ⎪
⎪ ⎪ ⎢ ⎪ ⎪ ⎪ ⎪ ⎪⎪ f ⎪⎪
⎨ F22 ⎬ = ⎢ 0 k 22 k 23 ⎥⎥ ⎨1⎬ = ⎨k 22 ⎬ = ⎨ 33 ⎬ (6-11)
⎪F ⎪ ⎢ 0 D
⎩ 32 ⎭ ⎣ k 32 k 33 ⎥⎦ ⎪⎩0⎪⎭ ⎪⎩k 32 ⎪⎭ ⎪ f 23 ⎪
⎪ ⎪
⎩⎪ D ⎭⎪

f 33 f 23
F12 = 0, F22 = , F32 =
D D
(6-12)
− f 33 f L − f 23
F42 = 0, F52 = − F22 = , F62 = F22 x − F32 = 33
D D
 F52
δ2= 1 F62
F12 x
A F42
F32 B
F22

Figura 6-4.Fuerzas provoc adas por un desplazamiento unitario en la direc c ión F2.

c) Finalmente, se aplica un desplazamiento unitario en A (figura 6-5), en la

dirección F3 , para obtener las fuerzas en el nodo A y por equilibrio
deducir las fuerzas en el nodo B.

⎧ ⎫
⎧ F16 ⎫ ⎡k11 0 0 ⎤ ⎧0⎫ ⎧ 0 ⎫ ⎪ 0 ⎪
⎪ ⎪ ⎢ ⎪ ⎪ ⎪ ⎪ ⎪⎪ f ⎪⎪
⎨ F22 ⎬ = ⎢ 0 k 22 k 23 ⎥⎥ ⎨0⎬ = ⎨k 23 ⎬ = ⎨ 23 ⎬ (6-13)
⎪F ⎪ ⎢ 0 D
⎩ 32 ⎭ ⎣ k 32 k 33 ⎥⎦ ⎪⎩1⎪⎭ ⎪⎩k 33 ⎪⎭ ⎪ f 22 ⎪
⎪ ⎪
⎪⎩ D ⎪⎭

δ3= 1 F53
F63
F13 x
F33
A B F43
F23

Figura 6-5.Fuerzas provoc adas por un desplazam iento unitario en la direc c ión F3.

f 22 f 23
F13 = 0, F33 = , F23 =
D D
(6-14)
f 23 f 23 L − f 22
F43 = 0, F53 = − , F63 =
D D

Por lo tanto la submatriz de rigidez AB es:

 ⎡ 1 ⎤
⎢− f 0 0 ⎥
⎢ 11 ⎥
f 33 f 33 L − f 23 ⎥
k AB =⎢ 0 − (6-15)
⎢ D D ⎥
⎢ f 23 f 23 L − f 22 ⎥
⎢ 0 − ⎥
⎣ D D ⎦

Y por simetría de la matriz de rigidez del elemento se tiene que:

⎡ 1 ⎤
⎢− f 0 0 ⎥
⎢ 11 ⎥
f 33 f 23 ⎥
k BA = k BA =⎢ 0 − −
T
(6-16)
⎢ D D ⎥
⎢ f 33 L − f 23 f 23 L − f 22 ⎥
⎢ 0 ⎥
⎣ D D ⎦

Dado que se conoce la submatriz de rigidez [k BA ] , se conocen también las

fuerzas producidas en el nodo A por los efectos de los desplazamientos unitarios
en B, por tanto se puede por equilibrio deducir las fuerzas y rigideces en el nodo
B, (figura 6-6, 6-7 y 6-8).

δ1= 1
F24 F54
F14 F64
F34 x
A F44
B
L

Figura 6-6.Fuerzas provoc adas por un desplazam iento unitario en la direc c ión F4.

1
F44 = F14 = , F24 = F34 = 0, F54 = F64 = 0 (6-17)
f 11
 F65 δ5= 1
F15 x
A B F45
F35
F25 F55

Figura 6-7.Fuerzas provoc adas por un desplazamiento unitario en la direc c ión F5.

f 33 f 23
F15 = 0, F25 = , F35 =
D D
(6-18)
f 33 f 23 − f 33 L
F45 = F15 = 0, F55 = F25 = , F65 = F35 − F25 L =
D D

F56
F66
F16 x
F36
A B F46
F26 δ6= 1

Figura 6-8.Fuerzas provoc adas por un desplazam iento unitario en la direc c ión F6.

f 33 L − f 23 f 23 L − f 22
F16 = 0, F26 = , F36 =
D D

f 23 − f 33 L
F16 = 0, F56 = − F56 = , F33 = F26 L − F36 (6-19)
D

( f 33 L − f 23 )L f 23 L − f 22 f 33 L2 + f 22 − 2 f 23 L
F33 = − =
D D D
La matriz resultante en B es:
⎡ 1 ⎤
⎢f 0 0 ⎥
⎢ 11 ⎥
f 23 − f 33 L
[k BB ] = ⎢⎢ 0 f 33 ⎥
⎥
(6-20)
D D
⎢ f 23 − f 33 L f 22 − 2 f 23 L + f 33 L2 ⎥
⎢ 0 ⎥
⎣ D D ⎦

Ensamblando las submatrices obtenidas, de las ecuaciones 6-8, 6-15, 6-16

y 6-20, se obtiene la matriz de rigideces de un elemento A-B de sección
variable:

⎡ 1 1 ⎤
⎢ f 0 0 − 0 0 ⎥
f11
⎢ 11 ⎥
⎢ 0 f 33 f 23 f 33 f 33 L − f 23 ⎥
0 −
⎢ D D D D ⎥
⎢ f 23 f 22 f 23 f 23 L − f 22 ⎥
[k ] [k ]AB ⎤ ⎢ 0 0 − ⎥
[k ] = ⎡⎢ AA =⎢ D D D D ⎥
⎣[k ]BA [k ]BB ⎥⎦ ⎢− 1 0 0
1
0 0 ⎥
⎢ f11 f11 ⎥
⎢ f 33 f 23 f 33 f 23 − f 33 L ⎥
⎢ 0 − − 0 ⎥
⎢ D D D D ⎥
⎢ f 33 L − f 23 f 23 L − f 22 f 23 − f 33 L f 22 − 2 f 23 L + f 33 L2 ⎥
⎢⎣ 0 D D
0
D D ⎥⎦
(6-21)

Renombrando los términos iguales, esta matriz se puede representar como:

⎡ k11 0 0 − k11 0 0 ⎤
⎢ 0 k 22 k 23 0 − k 22 k 26 ⎥⎥
⎢
⎢ 0 k 23 k 33 0 − k 23 k 36 ⎥
[k ] = ⎢ ⎥ (6-22)
⎢− k11 0 0 k11 0 0 ⎥
⎢ 0 − k 22 − k 23 0 k 22 − k 26 ⎥
⎢ ⎥
⎣⎢ 0 k 26 k 36 0 − k 26 k 66 ⎦⎥

Donde:

1 f 33 f 23 f 22
k11 = k 22 = , k 23 = , k 33 =
f 11 , D D D
(6-23)
f 33 L − f 23 f 23 L − f 22 f 33 L − 2 f 23 L + f 22
2
k 26 = , k 36 = , k 66 =
D D D
6.3 Ejemplo de una viga de sección variable rectangular.

Sea un elemento de sección variable rectangular llena, como el

mostrado en la figura 6-9. En este caso en particular, el peralte varía linealmente
a lo largo de la longitud, y tanto el área como el momento de inercia se
pueden expresar en función de x.

h2 − h1
h(x ) = h1 + x
L

(6-24)

3
⎛ h −h ⎞ bh( x)3 b ⎛ h −h ⎞
A( x) = bh( x) = b⎜ h1 + 2 1 x ⎟ e I ( x) = = ⎜ h1 + 2 1 x ⎟
⎝ L ⎠ 12 12 ⎝ L ⎠

h2

h1
b

b
x

Figura 6-9.Elem ento de sec c ión variable rec tangular llena.

De acuerdo con el capítulo anterior se puede obtener la matriz de

rigideces a partir de la matriz de flexibilidades, o de los términos de la matriz de
flexibilidades, f11 , f 22 , f 23 y f 33 :

L L L
dx dx 1 dx
f11 = ∫
EA( x ) ∫0 ∫
= = (6-25)
⎛ h − h1 ⎞ Eb 0 h − h1
0
Eb⎜ h1 + 2 x⎟ h1 + 2 x
⎝ L ⎠ L

Haciendo un cambio de variable:

 h2 − h1 h2 − h1
u = h1 + x, du = dx (6-26)
L L

L
dx = du , límites : x = 0 → u = h1 ; x = L → u = h2
h2 − h1

La integral se convierte en:

⎛ L ⎞
⎜ ⎟du
1 ⎜⎝ h2 − h1 ⎟⎠
h2
⎤
h2 h2

∫ =
L
∫
du
=
L
In(u )⎥ =
L
[In(h2 ) − In(h1 )]
Eb h1 u Eb(h2 − h1 ) h1 u Eb(h2 − h1 ) ⎦h Eb (h2 − h1 )
1

(6-27)

Por consiguiente se tiene que el primer término de la matriz de flexibilidades es:

L ⎡ ⎛ h2 ⎞⎤
f11 = ⎢ In⎜ ⎟⎟⎥ (6-28)
Eb(h2 − h1 ) ⎣ ⎜⎝ h1 ⎠⎦
El siguiente término es:

L L
x 2 dx k dx
f 22 = ∫ +∫ 1 (6-29)
0
EI ( x ) 0 GA( x )

Tomando el primer término de la integral:

L L
x 2 dx 12 x 2 dx
∫0 EI (x ) Eb ∫0 ⎛ h − h ⎞ 3
= (6-30)
⎜ h1 + 2 1
x⎟
⎝ L ⎠
Y haciendo un cambio de variables:

h2 − h1 u − h1 u 2 − 2uh1 + h12
u = h1 + x, x= L, x2 = L2 ,
L h2 − h1 (h2 − h1 ) 2

h2 − h1 L
du = dx, dx = du (6-31)
L h2 − h1

Con límites: sí x = 0 → u = h1 y si x = L → u = h2

La integral se convierte en:

 ⎛ u 2 − 2uh1 + h12 ⎞⎛ L ⎞
⎜ ⎟⎜ ⎟
12 L2 2 ⎜⎝ (h2 − h1 )
h 2 ⎟⎜ h − h ⎟du h2
(u )
− 2uh1 + h12 du
⎠⎝ 2 1 ⎠
2
12 L3
Eb h∫1 ∫
=
u3 Eb(h2 − h1 )
3
h1 u3

12 L3 ⎡ h2 u 2 du h2 h2
du ⎤ 12 L3 ⎛ 2 1 ⎞⎤
h2

In(u ) + 2h1 − h1
udu 1
= ⎢ ∫ 3 − 2h1 ∫ 3 + h1 ∫ 3 ⎥ = 3 ⎜ 2 ⎟⎥
2

Eb(h2 − h1 ) ( − ) ⎝ ⎠⎦ h1
3
⎢⎣ h1 u h1 u h1 u ⎥⎦ Eb h2 h1
u 2u

12 L3 ⎡ ⎛ h2 ⎞ ⎛ 1 1⎞ ⎛ 1 1 ⎞⎤
= ⎢ In⎜⎜ ⎟⎟ + 2h1 ⎜⎜ − ⎟⎟ − h12 ⎜⎜ 2 − 2 ⎟⎟⎥
Eb(h2 − h1 )
3
⎢⎣ ⎝ h1 ⎠ ⎝ h2 h1 ⎠ ⎝ 2h2 2h1 ⎠⎥⎦
(6-32)

El segundo término es una integral similar a f11 :

L L
k1 dx k1 dx
∫0 GA(x ) = G ∫ A(x )
0
(6-33)

Por lo tanto, simplificando términos se tiene:

12 L3 ⎡ ⎛ h2 ⎞ h1 h12 3⎤ k1 L ⎡ ⎛ h2 ⎞⎤
f 22 = ⎢ In ⎜ ⎟ + 2 − − ⎥ + ⎢ In⎜ ⎟⎟⎥ (6-34)
Eb(h2 − h1 ) ⎣ ⎜⎝ h1 ⎟⎠ h2 2h22 2 ⎦ Gb(h2 − h1 ) ⎣ ⎜⎝ h1
3
⎠⎦

El siguiente término de la matriz de flexibilidades f 23 es:

L L
xdx 12 xdx
f 23 = ∫ = ∫ (6-35)
EI ( x ) Eb 0 ⎛ h2 − h1 ⎞
3

⎜ h1 + x⎟
0

⎝ L ⎠

Haciendo un cambio de variables:

h2 − h1 u − h1
u = h1 + x x= L
L h2 − h1
h −h L
du = 2 1 dx dx = du (6-36)
L h2 − h1

límites : x = 0 → u = h1 ; x = L → u = h2
La integral se convierte en:

⎛ u − h1 ⎞⎛ L ⎞
⎜⎜ L ⎟⎟⎜⎜ ⎟⎟du
12
h2
⎝ h2 − h1 ⎠⎝ h2 − h1 ⎠ = 12 L2
h2
(u − h1 )du = 12 L2 ⎡ h2 udu h2
du ⎤
Eb h∫1 ∫ ⎢∫ 3 − h1∫ 3 ⎥
u3 Eb(h2 − h1 ) Eb(h2 − h1 )
2
h1 u3 2
⎢⎣ h1 u h1 u ⎥⎦
(6-37)
⎡ 1 1 h1 ⎛ 1 1 ⎞⎤
h2
12 L2 ⎛ 1 1 ⎞⎤ 12 L3
= ⎜ − + h1 2 ⎟⎥ = ⎢ − + ⎜⎜ 2 − 2 ⎟⎟⎥
Eb(h2 − h1 ) ⎝ u 2u ⎠⎦ h1 Eb(h2 − h1 )3
2
⎢⎣ h1 h2 2 ⎝ h2 h1 ⎠⎦⎥

Simplificando términos:

12 L2 ⎡ 1 1 h ⎤
f 23 = ⎢ − + 12 ⎥ (6-38)
Eb(h2 − h1 )
2
⎣ 2h1 h2 2h2 ⎦

Por último, el término f 33 se obtiene como sigue:

L L
dx 12 dx
f 33 = ∫ = ∫ (6-39)
EI ( x ) Eb 0 ⎛ h2 − h1 ⎞
3

⎜ h1 + x⎟
0

⎝ L ⎠
Haciendo el cambio de variable:

h2 − h1
u = h1 + x, límites : x = 0 → u = h1 ; x = L → u = h2
L

h2 − h1 L
du =
dx, dx = du (6-40)
L h2 − h1
La integral se convierte en:

⎛ L ⎞
⎜ ⎟du
12 ⎜⎝ h2 − h1 ⎟⎠ ⎛ 1 1 ⎞
h2
12 L ⎛ 1 ⎞⎤
h2 h2
12 L du 6L
Eb h∫1 Eb(h2 − h1 ) h∫1 u 3 Eb(h2 .h1 ) ⎝ 2u 2 ⎠⎦ h1 Eb(h2 − h1 ) ⎜⎝ h12 h22 ⎟⎠
= = ⎜ − ⎟ ⎥ = ⎜ − ⎟ (6-41)
u3
Simplificando términos se tiene:

6L ⎛ 1 1 ⎞
f 33 = ⎜⎜ 2 − 2 ⎟⎟ (6-42)
Eb(h2 − h1 ) ⎝ h1 h2 ⎠

En resumen, los valores de los términos de la matriz de flexión son:

L ⎛ ⎛ h2 ⎞⎞
f11 = ⎜ In⎜ ⎟⎟ ⎟
Eb(h2 − h1 ) ⎜⎝ ⎜⎝ h1 ⎠⎠
⎟

12 L3 ⎡ ⎛ h2 ⎞ h h2 3⎤ k1 L ⎛ ⎛ h2 ⎞⎞
f 22 = ⎢ In⎜⎜ ⎟⎟ + 2 1 − 1 2 − ⎥ + ⎜ In⎜ ⎟⎟ ⎟ (6-43)
Eb(h2 − h1 )
3
⎣ ⎝ h1 ⎠ h2 2h2 2 ⎦ Gb(h2 − h1 ) ⎜⎝ ⎜⎝ h1 ⎠⎠
⎟

12 L2 ⎡ 1 1 h ⎤
f 23 = ⎢ − + 12 ⎥
Eb(h2 − h1 )
2
⎣ 2h1 h2 2h2 ⎦

6L ⎛ 1 1 ⎞
f 33 = ⎜⎜ 2 − 2 ⎟⎟
Eb(h2 − h1 ) ⎝ h1 h2 ⎠

Está claro que teniendo los valores de las flexibilidades se puede obtener
la matriz de rigidez del elemento mediante las ecuaciones 6-23.

Se hace la observación que para otro tipo de secciones, como la sección

hueca o la sección tipo I, las integrales para obtener los términos de la matriz de
flexibilidades se complican, ya que tanto las funciones de área como de los
momentos de inercia son polinomios, que están en el denominador de la
integral; por lo tanto se puede obtener la matriz de rigideces de este tipo de
secciones mediante la composición de secciones macizas y la adición de las
matrices de rigideces de este tipo de secciones macizas, como se muestra en
las figuras siguientes:
 SECCIÓN HUECA S SECCIÓN LLENA A SECCIÓN LLENA B

h(x)
t1
= h(x) - h(x) - 2t2

t2

b b b - 2t1

[KS] = [KA] - [KB]

Figura 6-10.Elemento de sec c ión variable huec a o tubular c uadrada.

SECCIÓN HUECA S SECCIÓN LLENA A SECCIÓN LLENA B

h(x)
t1
= h(x) - h(x) - 2t2

t2

b b b - t1

[KS] = [KA] - [KB]

Figura 6-11.Elem ento de sec c ión variable tipo I.

Es necesario aclarar que en esta composición de secciones debe

hacerse la adición o diferencia con las matrices de rigidez y no con las matrices
de flexibilidad, ya que la flexibilidad es inversamente proporcional al área y al
momento de inercia, por lo cual no pueden sumarse o restarse áreas para
determinar la flexibilidad de un elemento, pero sí es válido sumarlas o restarlas
cuando se trata de rigideces, ya que la rigidez es directamente proporcional al
área y al momento de inercia.

Como comprobación de esta composición de áreas, se hicieron

integrales numéricas para casos específicos y se verificaron los resultados
obtenidos mediante la resta de matrices de rigideces de secciones macizas. Los
resultados obtenidos por ambos métodos son similares.
6.4 Conclusiones.

La matriz de rigidez de un elemento de sección variable ha sido deducida

a través de la matriz de flexibilidades. En particular se ejemplifica la deducción
para un elemento de sección rectangular maciza o llena. Esta solución puede
fácilmente ser extendida a otro tipo de secciones como las secciones huecas o
las tipo I, de uso muy frecuente sobre todo en elementos de acero, sin
necesidad de deducir matrices específicas de cada sección, ya que utilizando
la matriz de rigidez de una sección maciza, se pueden obtener otro tipo de
sección con combinaciones de la misma. Finalmente, cabe señalar que el
método se puede fácilmente sistematizar e implementar en programas de
análisis de uso común.