ELEMENTOS DE SECCIÓN VARIALE

M.I. Octavio García Domínguez

INTRODUCCIÓN

Se plantea el problema de encontrar la relación entre los desplazamientos y las fuerzas en el extremo libre de una viga en
voladizo de sección variable. Tómese el extremo A.

MA

A
B dyA

EI x Ax
MB
FxA FxB

FyA
FyB
dxA

Esta relación se puede establecer en forma similar a la

ec. fundamental del método de rigideces, es decir: f

FA = KAA  dA
MA MB
La matriz KAA contiene las fuerzas o rigideces debidas a FxB

desplazamientos unitarios aplicados en el extremo A de la viga FyA

FyB

Su obtención es posible hacerla por columnas, calculando las fuerzas asociadas a los desplazamientos unitarios, aplicados de
manera independiente. Sin embargo, no es un problema fácil dado que el elemento estructural tiene sección variable en su
longitud. Para este caso, se asume que la variación de la sección está dada por la variación del peralte, manteniendo
constente el ancho del elemento.

1
Para facilitar el proceso de cálculo de la matriz K AA se plantea la siguiente relación dA = fAA FA donde fAA = KAA
[Link] es la matriz de flexibilidades en el extremo A y contiene, por columnas,
los desplamientos producidos por cargas unitarias en ese extremo.

FxA FyA MA

 f11 f12 f13 

 
 f21 f22 f23 
f f32 f33 
 31

Pág.- 1 25/08/2014 06:42 PM

Este proceso se puede realizar de forma simple, utilizando el método de la carga unidad, cuya ecuación se presenta a
continuación:

   
N n  V v M m T t
Δ = dx   dx   dx   dx
 E A G  Ac  E I  GJ
   


Según el método, las literales mayúsculas (N,V, M y T ) representan las ecs de los diagramas de elementos mecánicos
debidos a la carga real, mientras que las literales minúsculas (n,v,m y t), corresponden a los diagrmas de elementos
mecánicos debidos a una carga unitaria aplicada en el punto y dirección del deplazamiento de interés. Es de interés aclara
que para nuestro caso, la carga real (por definición de la matriz de flexibilidad), es también unitaria, lo cual facilita aún más
los cálculos.
Obtención de la primera columna de [Link] .

Se aplica Fx A =1y se calculará dxA, dy A y φA. De acuerdo con el método, se deberán aplicar fuerzas unitarias en dirección
de cada uno de estos desplazamientos y desarollar las integrales correpondientes:

f , MA = 1
DxA, fxA = 1
FxA

DyA, f yA = 1

N n n n

V v v v

M m m m

FxA =1 N  1 fxA =1 n = 1 fyA =1 n=0 mA =1 n=0

V  0 v =0 v =1 v =0
M  0 m=0 m=x m = 1

Para dxA n   1 v  0 m  0

    
N n  V v M m 1 1
dxA   dx   dx   dx   dx f11 =  dx
 E  Ax G  Ac x  E  I x   E  Ax  E  Ax
    


Para dyA n  0 v  1 m  x

  
N n  V v M m f21 = 0
dyA   dx   dx   dx  0
 E  Ax G  Ac x  E  I x 
  


Para φA n  0 v  0 m  1

  
N n  V v M m f31 = 0
ϕA   dx   dx   dx  0
 E  Ax G  Ac x  E  I x 
  


Pág.- 2 25/08/2014 06:42 PM

Obtención de la segunda columna de [Link] . Se aplica Fy A =1y se calculará dxA, dy A y φA.

f , MA = 1
DxA, fxA = 1

FyA DyA, f yA = 1

n n n
N

v v v
V

m m m
M

FyA =1 N  0 fxA =1 n = 1 fyA =1 n=0 mA =1 n=0

V  1 v =0 v =1 v =0
M  x m=0 m=x m = 1

Para dxA n   1 v  0 m  0

  
N n  V v M m
dxA   dx   dx   dx  0 f12 = 0
 E  Ax G  Ac x  E  I x 
  

Para dyA n  0 v  1 m  x  x

   
N n V v M m  x
2 
  1
dyA   dx   dx   dx   dx   dx f22
 E  Ax G  Ac x  E  I x   E  I x  G  Ac x
    
 

Para φA n  0 v  0 m  1
Sin considerar el
 signo nefativo
 1
   x 
N n  V v M m  E x
ϕA   dx   dx   dx   dx f32 =  dx
 E  Ax G  Ac x  E  I x   I x   E  I x 
    


Obtención de la tercera columna de [Link] . Se aplica MA =1y se calculará dxA, dy A y φA.

f , MA = 1
DxA, fxA = 1
MA

DyA, fyA = 1

n n n
N

v v v
V

m m m
M

MA =1 N  0 fxA =1 n = 1 fyA =1 n=0 mA =1 n=0

V  0 v =0 v =1 v =0
M   1 m=0 m=x m = 1

Pág.- 3 25/08/2014 06:42 PM

 Para dxA n   1 v  0 m  0

  
N n  V v M m
dxA   dx   dx   dx  0 f13 = 0
 E  Ax G  Ac x  E  I x 
  


Para dyA n  0 v  1 m  x  x


 1
   x 
N n  V v M m  E x Sin considerar el
dyA   dx   dx   dx    dx f23 =  dx
 E  Ax G  Ac x  E  I x   I x   E  I x  signo nefativo
    


Para φA n  0 v  0 m  1

    
N n  V v M m 1 1
ϕA   dx   dx   dx   dx f33 =  dx
 E  Ax G  Ac x  E  I x   E  I x   E  I x 
    


1
Determinadas la matriz de flexibilidades fAA, la matriz de rigidez K AA resulta KAA = fAA

dXA dYA ϕA

 1 0  FXA
0
f Donde :
 f11 0 0   11 
  1
 f33 f23  F Det  f22  f33  f23  f32
fAA =  0 f22 f23  KAA = fAA = 0  YA
Det Det
 0 f   
 32 f33   f32 f22 
MA
 0 Det 
 Det

Ahora la matriz K BA se obtiene por equilibrio

MA
f

dyA
MB MA
1/f11 1/f11

FyA FyA
dxA
FyB FyB

1 1
Para dxA tenemos:
 Fx = 0 f11
 FBA = 0 FxB =
f11

Para dy A tenemos:

f33  f33
 Fy = 0
Det
 FyB = 0 FyB =
Det

Pág.- 4 25/08/2014 06:42 PM

 f33 f32 f33  L  f32
 MB = 0 MB 
Det
L 
Det
=0 MB =
Det

Para A tenemos:

f23  f23
 Fy = 0
Det
 FYB = 0 FyB =
Det

f23 f22 f23  L  f22

 MB = 0 MB 
Det
L 
Det
=0 MB =
Det

dXA dYA ϕA
Entonces K BA queda cmo sigue:
1
FXB  0 0
f11
 
  f33  f23 
FYB  0  = KBA
Det Det
 
 f33  L  f32 f23  L  f22 
 0 
MB  Det Det 

Si se tiene una sección constante la matrices fAA y k AA, adquieren la siguiente forma:

L L 2 L
 1 L  x L  1 L
f11   dx  f32   dx  f23  f32 f33   dx 
 E A A E  E I 2 E I  E I E I
0 0 0

L
 L
 x
2
 1 L
3
L L
3

6  1  ν I 
f22   dx   dx   f22   1  C  Donde : c 
E I  G  Ac 3  E  I Ac G 3 E  I Ac  L
2
0 0

 L 0

0
 1 0 0 
f
 f11 0 0   A E   k11 0 0   11 
   L
3
L
2     f33 f23 
fAA =  0 f22 f23  = 0  1  c  kAA =  0 k22 k23  = 0 
3 E  I 2 E I Det Det
 0 f     0 k   
 32 f33    32 k33 
L   f22 
2 f32
L
 0 E  I   0 Det 
 2 E I  Det

 L 3   L   L2  
2
L  4 C  1
4
L
4
L  1  4C 
4

Det    1  C         factor  Det =  4  1  C   3 =
 3 EI   EI   2 EI   12  EI
2
12 EI 2 12 EI 2

Pág.- 5 25/08/2014 06:42 PM

  Det   A E
f 0 0  0 0 
 11   L 
L  1  4C 
4  L L
2  1  12  E  I 6 E I 
kAA   0  0 
Det    Det L  4  C  1
3
L  4 C  1
2 
12 E  I 2  E I 2 E I
  
 2 3   6 E I 4  E  I  C  1 
L L
 1  C   
 0  0 L  4  C  1
2 L  4 C  1 
 2  E  I 3 E  I  

Entonces K BA queda cmo sigue:

 A E 
 1 0 0   L
 0 0 
  
f11
 12  E  I 6 E I
  f33  f23   0   
 L  4  C  1 L  4 C  1 
3 2
KBA   0
 Det Det
  
L  C  1  
2 2  L
3 3
 f33  L  f32 f23  L  f22   12  E  I    
 0   0 6 E I 2 E I 3 E  I  
 Det Det   
L  4 C  1 L  4 C  1
2 4
 

  A E 0 0

 L 
  12  E  I 6  E  I 
KBA   0 
L  4 C  1 L  4 C  1
3 2
 
 6 E I 2  E  I  1  2  C  
 0 L  4  C  1
2 L  4 C  1 


De las matrices de rigdez KBA y KAA se puede obtener la matriz de rigidez angular para un elemento cualquiera

 4 E  I  C  1 
2 E I 1  2 C 

L  4  C  1 L  4 C  1
a) Para sección constante k   
 2  E  I  1  2  C  4  E  I  C  1 
 L  4 C  1 L  4  C  1 


b) Para sección variable faltaría obtener rBB

θA = 1 θB = 1

 f22 f23  L  f22 

 
 rAA rAB   Det Det 
k= =
r r
 BA BB   f23  L  f22 
 rBB 
 Det 

Pág.- 6 25/08/2014 06:42 PM

Por trabajos recíprocos: El momento
asociado a un desplazamiento vertical
unitario es igual al cortante asociado a
un giro unitario.

De manera análoga para el extremo B se tiene:

Resumiendo, la matriz de rigidez angular resulta :

 f22 f23  L  f22 

 
 AA AB  
r r Det Det

    2 
r r
 BA BB   23 f  L  f f
22 33  L  2  f23  L  f22

 Det Det 

Pág.- 7 25/08/2014 06:42 PM

 Obtención de la 2a columna de [K AA] en función de la matriz de rigidez angular:

El cortante en la deformada de la rigidez angular

se obtiene por estática como lasuma de momentos
entre la longitud

Por trabajos recíprocos, el momento debido a un

desplazamiento vertical unitario es igual al cortante
debido a una rotación unitaria.

El cortante en la deformada de rigidez angular y por

tanto el momento en la deformada de rigidez lineal
dependen de la rigidez angular en el extremo
correspondiente. Por lo anterior, le 2a columna de [KAA]
se puede calcular como se indica:

Comprobando :

2
rAA  rAB f22  f23  L  f22 f23 rBB  rAB f33  L  2f23  L  f22  f32  L  f22 f33  L  f23
= = = =
L L  Det Der L L  Det Der

EJEMPLOS : Determine las matrices [fAA], [K AA] y la matriz de rigidez angular [k] de los siguiientes elementos, c=0.

EJEMPLO 1

L
2 L 2 2
 x  x 5 L 5 L
f32   dx   dx  f23  f32 
E I  2 E I 16  E  I 16  E  I
0
L
2 L
L
L 2 2  2 3
L  x  x 3 L
2  f22   dx   dx 
 1 1 3 L E I 2 E I 16  E  I
f33   dx   dx  0 L
E I  2 E I 4 E  I 
0 2
L
2

 L 0 0 
 f11 0 0   A E 
   3 L
3
5 L
2 
11 L
4
fAA   0 f22 f23    0   Det  f22  f33  f23  f32 
16  E  I 16  E  I 2 2
 0 f    256  E I
 32 f33   2 
5 L 3 L
 0  
 16  E  I 4 E  I 

Pág.- 8 25/08/2014 06:42 PM

  1 0   A E
0 0 0 
f
 11   L 
 f33 f23   192  E  I 80  E  I 
kAA   0  0 11 L
3
11 L
2 
Det Det
   
 f32 f22  
0
80  E  I 48 E  I 
 0 Det  
2 11 L 
 Det 11 L 

 f22 f23  L  f22   48 E  I

  32  E  I 
 rAA rAB   Det Det
   11 L 11 L 
    2   
 rAB rBB   f23  L  f22 f33  L  2  f23  L  f22   32  E  I 80  E  I

 
11 L 11 L 
 Det Det

Para la 2a col de kAA y KBA

rAA  rAB 80  E  I 80
MA    7.27
L 2 11
11 L

rBB  rAB 112  E  I 112

MB    10.18
L 2 11
11 L

MA  MB 192  E  I 192
VA    17.45
L 3 11
11 L

EJEMPLO 2.

L L
4 2 L 2 2
 x  x  x 13 L 13 L
f32   dx   dx   dx  f23  f32 
E I 2 E I  3 E  I 64 E  I 64 E  I
0 L
 L
4 2
L L
4 2 L
 1  1  1 13 L
f33   dx   dx   dx 
E I 2 E I  3 E  I 24 E  I
0 L
 L
4 2

L L
L
4 2 2 2  2 3
 x  x  x 139  L
f22   dx   dx   dx 
E I 2 E I 3 E  I 1152  E  I
0 L L
 
4 2

Pág.- 9 25/08/2014 06:42 PM

  L 0 0

 f11 0 0   A E 
   139  L
3
13 L 
2
2665  L
4
fAA   0 f22 f23    0   Det  f22  f33  f23  f32 
1152  E  I 64 E  I 2 2
 0 f    110592  E I
 32 f33  
13 L 
2
13 L
 0  
 64 E  I 24 E  I 

 1 0 0   A E 0 0 
f
 11   L 
 f33 f23   4608  E  I 1728  E  I 
kAA   0  0 205  L
3
205  L
2 
Det Det
   
 f32 f22  
0
1728  E  I 13344 E  I 
 0 Det  
2 2665  L 
 Det 205  L 

 f22 f23  L  f22 

   13344 E  I 1824  E  I 
 rAA rAB   Det Det
 simplify  2665  L 533 L 
    2 
 
r r
 AB BB   f23  L  f22 f33  L  2  f23  L  f22   1824  E  I 5664  E  I

 Det Det   533 L 533 L 

 5.007  E  I 3.422  E  I 
 rAA rAB   L L  1824 13344 5664
     533
 3.422
2665
 5.007
533
 10.627
 rAB rBB   3.422  E  I 10.627  E  I

 L L 

Para la 2a col de kAA y KBA

rAA  rAB 8.429  E  I

MA  simplify 
L 2
L

rBB  rAB 14.049  E  I

MB  simplify 
L 2
L

MA  MB 22.478 E  I
VA  simplify 
L 3
L

Pág.- 10 25/08/2014 06:42 PM

 EJEMPLO 3.

A
Para el miembro de sección
h0
B 3h 0
variable determinar la matriz
de rigidez angular

3
b b  ho
Ao = b  ho Io =
12
b
L

ho ho  x
La ecuación lineal que gobierna el aumento del peralte es: y  ho =  x  L solve y 
L L

 ho  b  h o  L  2  x  2 x 
entonces la sección está dada por: Ax  b  ho  2   x  simplify  Ax  Ao  1  
 L  L  L 

b  h o  L  2  x 
3 3 3 3
 ho  2 x 
I x  Io  1 
1
El momento de inercia es: I x   b  ho  2   x  simplify  
12  L  12  L
3  L 

Resolviendo para las

flexibilidades:
L
 1 L  ln3
f11   dx simplify 
 E  Ax 2  Ao  E
0

L

L   9  ln3  8 
2 3
x
f22   dx simplify 
 E  I x  72  E  Io
0

L 2 2
 x L L
f32   dx  f23  f32 
 E  I x  18 E  Io 18 E  Io
0

L
 1 2 L 2 L
f33   dx  if L  0  L  0   2  L  L f33 
 E  I x  9  E  Io 9  E  Io
0
undefined otherwise

 L  ln3 0 0

 2  Ao  E 
 11
f 0 0   
  L   9  ln3  8
3
L
2
fAA   0 f22 f23  simplify   0  
 72  E  Io 18 E  Io 
 0 f f   
 32 33  2
 L 2 L 

 0 18 E  Io 9  E  Io 
 

L   ln3  1
4
Det  f22  f33  f23  f32 simplify 
2 2
36  E  Io

Pág.- 11 25/08/2014 06:42 PM

  1   2  Ao  E 
0 0 0 0
f
 11   L  ln3 
  8  E  Io 2  E  Io 
f33 f23 
kAA   0  simplify  0 
L   ln3  1 L   ln3  1 
3 2
 Det Det
 
 
 f32 f22  2  E  Io E  Io   9  ln3  8 
 0  0 
 Det Det   2  L   ln3  1  
L   ln3  1 
2
 

 f22 f23  L  f22   E  Io   9  ln3  8 3  E  Io   ln27  4 

    
 rAA rAB   Det Det
 simplify 2  L   ln3  1  2  L   ln3  1 
     
r r 2   3 E  Io   ln27  4 9  E  Io  ln3 
 AB BB   f23  L  f22 f33  L  2  f23  L  f22 
 
 Det Det   2  L   ln3  1 2  L   ln3  1  

 E  Io E  Io   9  ln3  8
 9.57  10.71   9.57
 rAA rAB  L L  2   ln3  1
   
 rAB rBB   10.71 E  Io E  Io 
 50.13 
 L L 

Para la 2a col de kAA y kBA

rAA  rAB 20.28 E  Io rBB  rAB 60.84 E  Io

MA  simplify  MB  simplify 
L 2 L 2
L L

60.84 E  Io 20.28 E  Io

MA  MB 2 2 81.12  E  Io
L L
VA   simplify 
L L 3
L

Pág.- 12 25/08/2014 06:42 PM

MOMENTOS DE EMPOTRAMIENTO
Los momentos de empotramiento de una viga de sección variable, no solo son función de la carga sino también de las
rigideces angulares en los extremos de la barra. El problema se puede atacar por superposición, descomponiendo la viga
empotrada en 2 vigas simplemente apoyadas, la primera con la carga aplicada y la segunda sólo con los momentos de
empotramiento que se desean conocer.

Por efecto de la carga aplicada, la primer viga sufrirá giros (θAo y θBo) en sus extremos. Si se aplican momentos en los
extremos de la segunda viga, capaces de restituir los giros a su condición inicial (giros nulos), los momentos aplicados será
precisamente los momentos de empotramiento buscados.

El problema se centra ahora en determinar los momentos que produzcan los giros (θAo y θBo) en sentido contrario, hasta
hacerlos nulos. Con base en lo anterior, la segunda viga simplemente apoyada es equivalente a una viga empotrada, cuyos
momentos pueden obtenerse mediante el principio de superposición. Para ello se considera el estudio de 2 vigas en
voladizo con rigideces angulares conocidas, a las que se les aplicará un momento para generar el giro requerido en cada
extremo. Sumando los momentos en cada viga resulta

MA  rAA θAo  rAB  θBo

MB  rBB  θBo  rAB  θAo

El problema querdará resuelto si se conocen los giros (θAo y θBo) producidos por la carga externa sobre la viga
simplemente apoyada.
Haciendo uso del principio de la carga unidad, puede obenerse una expresión gerenarl para el cálculo de los giros.

Cálculo de θB

Pág.- 13 25/08/2014 06:42 PM

 x
m 
L
L
 Mo
1 
θBo   x dx
L  E  I x 
0

Para θA

L

 Mo  1  
x

θAo  
 L  dx

0
E I x 

L
 M
 o
θAo  dx  θBo
 E  I x 
0

EJEMPLO 1. Obtener los momentos de empotramiento

 48 EI 32  EI 
 rAA rAB   11 L 11 L 
    
 rAB rBB   32  EI 80  EI

 11 L 11 L 
2
w Lx w x
Mo  
2 2

L L
L L
2 Mo  Mo 3 2 M  M 3
1  1  7 L w  o  o 9 L w
θBo    x  dx    x  dx  θAo   dx   dx  θBo 
L EI L 2  EI 256  EI EI 2EI 256  EI
0  0 
L L
2 2

2 2
13 L  w 17  L  w
MA  rAA θAo  rAB  θBo  MB  rBB  θBo  rAB  θAo 
176 176

2 2
w L w L
MA = MB =
13.54 10.35

Pág.- 14 25/08/2014 06:42 PM

 EJEMPLO 2. Obtener los momentos de empotramiento

x P
Para una viga con carga puntual en el centro Mp1  P  y Mp2   L  x 
2 2
2
w Lx w x
Mo  
2 2

Sumando los momento debidos a las cargas

P x P  L  x 
Mo1  Mp1  y Mo2  Mp2 
2 2

 E  Io E  Io 
 9.57  10.71 
3  rAA rAB  L L 
I x  Io  1 
2 x 
    
  rAB rBB   10.71 E  Io E  Io 
L 
 50.13 
 L L 
L

Mo 16  L  w  15  L  w  ln3
3 3
1 
θBo   x dx simplify  
L  E  I x  96  E  Io
0

L
 M
L  w   7  ln3  8
3
 o
θAo  dx  θBo simplify  
 E  I x  32  E  Io
0

simplify 2 1
MA  rAA θAo  rAB  θBo  0.03916  L  w  25.536
float 4 0.03916

Pág.- 15 25/08/2014 06:42 PM

 simplify 2 1
MB  rBB  θBo  rAB  θAo  0.1466  L  w  6.821
float 4 0.1466

EJEMPLO 3. Obtener los momentos de empotramiento

 5.007  EI 3.422  EI 
 rAA rAB   L L 
    
 rAB rBB   3.422  EI 10.627  EI

 L L 
2
w Lx w x
Mo  
2 2

L L
L
4 Mo 2 Mo  Mo 3
1  1  1  631 L  w
θBo    x  dx    x  dx    x  dx 
L EI L 2EI L 3EI 36864 EI
0 L 
 L
4 2

L L
L
4 M 2 M  M 3
 o  o  o 889  L  w
θAo   dx   dx   dx  θBo 
EI 2EI 3EI 36864 EI
0 L 
 L
4 2

simplify 2 1
MA  rAA θAo  rAB  θBo  0.06217  L  w  16.085
float 4 0.06217

simplify 2 1
MB  rBB  θBo  rAB  θAo  0.09938 L  w  10.062
float 4 0.09938

Pág.- 16 25/08/2014 06:42 PM