IC0611 Análisis Matricial De Estructuras Reticulares

Especialización en Ingenierı́a Sismo-Resistente

Profesor: Juan Camilo Molina Villegas
Lista de formulas

Transformación de las cargas externas de coordenadas globales a coordenadas locales

La carga externa horizontal (positiva en dirección del eje global x) se transforma a los ejes locales
empleando las siguientes ecuaciones:

pE (x′E ) = QE (x′E ) ∣sin(θE )∣ cos(θE ) 0 < x′E < LE (1a)

qE (x′E ) = −QE (x′E ) ∣sin(θE )∣ sin(θE ) 0 < x′E < LE (1b)

La carga externa vertical (positiva en dirección del eje global y) se transforma a los ejes locales empleado
las siguientes ecuaciones:

pE (x′E ) = RE (x′E ) ∣cos(θE )∣ sin(θE ) 0 < x′E < LE (2a)

qE (x′E ) = RE (x′E ) ∣cos(θE )∣ cos(θE ) 0< x′E < LE (2b)

Transformación de formulación del método de rigidez de coordenadas locales a coordenadas

globales

La formulación del método de rigidez en coordenadas globales es:

T T
{FE } = [TE ] [KE
′
] [TE ] {∆E } + [TE ] {FE′f } = [KE ] {∆E } + {FEf } (3)

Donde:
{FE′f }: Vector de fuerzas de empotramiento del elemento en coordenadas locales.
{FEf } = [TE ]T {FE′f }: Vector de fuerzas de empotramiento del elemento en coordenadas globales.
[KE′
]: Matriz de rigidez del elemento en coordenadas locales.
T
[KE ] = [TE ] [KE ′
] [TE ]: Matriz de rigidez del elemento en coordenadas globales.
[TE ]: Matriz de transformación del elemento.

Resorte

Fuerzas internas
PE (x′E ) = kE (u′E
j − ui )
′E
(4)
Sistema de ecuaciones del método de rigidez con funciones de Green
⎧ ⎫ ⎡ ⎧ ′E ⎫
⎪
⎪F Xi′E ⎪
⎪ ⎢ kE −kE ⎤⎥ ⎪
⎪u ⎪ ⎪
⎨ ⎬=⎢
′E ⎪ ⎢
⎥ ⎨ i′E ⎬ (5)
⎪
⎪F Xj ⎪ ⎣−kE ⎥ ⎪
kE ⎦ ⎪uj ⎪
⎩ ⎭ ⎩ ⎪ ⎭
Barra

Ecuación diferencial gobernante

d2 uE ′
AE EE (xE ) = −pE (x′E ) (6)
dx′2E
Fuerzas internas
duE ′
PE (x′E ) = AE EE (x ) (7)
dx′E E
Funciones de forma

x′E
ψ1E (x′E ) = 1 − (8a)
LE
x′E
ψ4E (x′E ) = (8b)
LE

1
Funciones de Green
⎧
⎪ LE
⎪
⎪
⎪ ψ4E (x′E )ψ1E (ξE
′
) 0 < x′E ≤ ξE
′
⎪
⎪ E EE
A
Gxx (xE , ξE ) = ⎨
E ′ ′
(9)
⎪
⎪
⎪ LE
⎪
⎪
⎪ ψ1E (x′E )ψ4E (ξE
′
) ξE
′
≤ x′E < LE
⎩ AE EE
Sistema de ecuaciones del método de rigidez con funciones de Green

⎧
⎪ ⎫
⎪ ⎡ ⎤⎧⎪ ⎫
⎪ ⎧
⎪ ⎫
⎪
⎪
⎪ ⎪ ⎢ AE E E AE EE ⎥⎥ ⎪ ⎪ ⎪ ⎪
⎪
⎪ ′E ⎪
⎪
⎪ ⎢ ⎪
⎪
⎪ ′E ⎪
⎪
⎪ ⎪
⎪
⎪ LE E ⎪
⎪
⎪
⎪ ⎪ − ⎥⎪ ⎪ ⎪ (x )p (x )dx⎪
⎪ ⎢⎢ LE ∫
⎪F X ⎪ u ψ ′ ′
⎪ LE ⎥⎥ ⎪ ⎪ ⎪ ⎪ ⎪
⎪ 0 ⎪
⎪
i i 1 E E E
⎨ ⎬=⎢ ⎨ ⎬−⎨ ⎬ (10)
⎪
⎪ ⎪
⎪ ⎢ ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ⎪ ⎢ AE EE AE EE ⎥⎥ ⎪ ⎪ ⎪ ⎪
⎪
⎪ ′E ⎪
⎪
⎪ ⎢− ⎪
⎪
⎪ ⎪
′E ⎪
⎪ ⎪
⎪
⎪ L ⎪
⎪
⎪
⎪
⎪F X ⎪
⎪ ⎢ ⎥⎪ u ⎪ ⎪∫
E
ψ E
(x ′
)p (x ′
)dx⎪
LE ⎦ ⎪ ⎪ ⎪ ⎪
E
⎭ ⎣ LE
j j 4 E E
⎩ ⎩ ⎭ ⎩ ⎭
0

EULER-BERNOULLI
Viga uniforme (Euler-Bernoulli)

Ecuación diferencial gobernante

d4 vE ′
EE IE (xE ) = qE (x′E ) (11)
dx′4E
Cálculo de fuerzas internas

d3 vE ′
VE (x′E ) = −EE IE (xE ) (12a)
dx′3E
d2 vE ′
ME (x′E ) = EE IE (xE ) (12b)
dx′2E

Funciones de forma
2 3
x′E x′
ψ2E (x′E ) = 1 − 3 ( ) + 2( E ) (13a)
LE LE
⎡ ′ 3⎤
⎢x ⎥
2
x′ x′
ψ3E (x′E ) = ⎢⎢ E − 2 ( E ) + ( E ) ⎥⎥ LE (13b)
⎢ LE LE LE ⎥
⎣ ⎦
2 3
x ′
x ′
ψ5E (x′E ) = 3 ( E ) − 2 ( E ) (13c)
LE LE
⎡ 3⎤
⎢ x ⎥
2
′
x′
ψ6E (x′E ) = ⎢⎢− ( E ) + ( E ) ⎥⎥ LE (13d)
⎢ LE LE ⎥
⎣ ⎦
Funciones de Green
⎧
⎪ ⎡ 2 E ′ ⎤
⎢ x′E ψ3 (ξE ) ⎥⎥
3
⎪
⎪
⎪
L3E ⎢− ( x′E
⎪ ⎢ L ) ψ2 (ξE ) + 3 ( L ) 0 < x′E ≤ ξE
E ′
LE ⎥⎥
′
⎪
⎪
⎪ ⎢
⎪ 6EE IE ⎣ E E
⎦
yy (xE , ξE ) = ⎨
GE ′ ′
⎡ 2 E ′ ⎤
(14)
⎪
⎪ ⎢
3
⎥
⎪ L3E ⎢− (1 − E ) ψ E (ξ ′ ) − 3 (1 − E ) 6 (ξE ) ⎥
x x ψ
′ ′
⎪
⎪
⎪ ≤ < LE
⎢ LE ⎥⎥
ξE x′E
′
⎪
⎪ 6EE IE
⎪ ⎢ LE 5 E
LE
⎩ ⎣ ⎦
Sistema de ecuaciones del método de rigidez con funciones de Green

⎧ ⎫ ⎡ ⎤⎧ ⎫ ⎧ ⎫
⎪
⎪ ⎪
⎪ ⎢ 12E I ⎥⎪ ⎪ ⎪ ⎪
⎪ ⎪ ⎢ ⎥⎪ ⎪ ⎪ ⎪
⎥⎪ ′E ⎪ ⎪ ⎪
6EE IE 12EE IE 6EE IE
⎪
⎪
⎪ ⎪
⎪ ⎢
E E
− ⎪
⎪ ⎪
⎪ ⎪
⎪ LE E
(x )q (x )dx ⎪
⎪
F Yi ⎪ ⎥⎪ v ⎪ ⎪∫ ψ ⎪
′E ′ ′
⎪
⎪ ⎪
⎪ ⎢ L 3 L2E L3E L2E ⎪
⎪ i ⎪
⎪ ⎪
⎪ 0 2 E E E ⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ E ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ ⎥
2EE IE ⎥ ⎪⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪
⎪ ⎢ 6E E IE 4EE IE 6EE IE
⎪ Mi′E ⎪⎪ ⎢ − ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
LE E
⎪
⎪ ⎪
⎪ ⎢ ⎥ ⎪
⎪ θ ′E
⎪
⎪ ⎪
⎪∫ ψ (x ′
)q (x ′
)dx ⎪
⎪
⎪ ⎪ ⎢ L2E LE L2E LE ⎥ ⎪ ⎪ ⎪ i 0 3 E E E ⎪
⎨ ⎬=⎢ ⎨
⎥⎪ ⎪ ⎪ ⎬ − ⎨ ⎬ (15)
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ 12E I ⎥ ⎪ ⎪ ⎪ ⎪
⎪
⎪
⎪ ⎪
′E ⎪ ⎢ 6EE IE 12EE IE 6EE IE ⎥ ⎪
⎪
⎪ ⎪
⎪
′E ⎪
⎪
⎪
⎪ ⎪
⎪
⎪
⎪
⎪F Y ⎪
j ⎪ ⎢−
⎪ ⎢ E E
− − ⎥ ⎪
⎪v ⎪
⎪
j ⎪ ⎪∫0⎪
⎪
LE E
ψ (x ′
)q (x ′
)dx ⎪
⎪
⎪
⎪ ⎪ 3 L2E L3E 2 ⎥ ⎪
⎪ ⎪ ⎪ 5 E E E ⎪
⎪
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ L LE ⎥ ⎪ ⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ E
⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪ ⎢ 4EE IE ⎥ ⎪ ⎪ ⎪ ⎪
Mj′E ⎪ ⎪ ⎪ ⎪ ⎪
⎪ ⎪ 6E E IE 2EE IE 6EE IE ⎪ ⎪ ⎪ ⎪
⎪
⎪ ⎪
⎪ ⎢ − ⎥⎪⎪ θ ′E
⎪
⎪ ⎪
⎪∫
L E
ψ E
(x ′
)q (x ′
)dx ⎪
⎪
⎩ ⎢
⎭ ⎣ L2E LE L2E ⎥
LE ⎦ ⎩ ⎭ ⎩
j 0 6 E E E
⎭

2
Campo de desplazamiento
vE (x′E ) = vE
h
(x′E ) + vE
f
(x′E ) (16)
Donde:
h
vE (x′E ) = ψ2E (x′E )vi′E + ψ3E (x′E )θi′E + ψ5E (x′E )vj′E + ψ6E (x′E )θj′E (17a)
x′E LE
f
vE (x′E ) = ∫ yy (xE , ξE )qE (ξE )dξE + ∫
GII ′ ′ ′ ′
GIyy (x′E , ξE
′
)qE (ξE
′
)dξE
′
(17b)
0 x′E

Viga nouniforme (Euler-Bernoulli)

Ver el artı́culo [1].

Viga uniforme sobre fundación de Winkler (Euler-Bernoulli)

Ecuación diferencial gobernante

La ecuación diferencial gobernante para una viga sobre fundación flexible es:

d4 vE ′
EE IE (xE ) + kE vE (x′E ) = qE (x′E ) (18)
dx′4E

Donde kE es la rigidez por unidad de longitud del suelo.

O dividiendo a ambos lados de esta por EE IE :

d4 vE ′ qE (x′E )
(xE ) + 4λ4E vE (x′E ) = (19)
dxE ′4 EE IE
√
kE
Donde λE = 4
.
4EE IE
Funciones de forma

1
2 (xE ) =
ΨE ′
[(c2 − ch2 ) sin(λE x′E ) sinh(λE x′E ) + (s ⋅ c + sh ⋅ ch) sin(λE x′E ) cosh(λE x′E )
sh − s2
2

−(s ⋅ c + sh ⋅ ch) cos(λE x′E ) sinh(λE x′E ) + (sh2 − s2 ) cos(λE x′E ) cosh(λE x′E )] (20)

1
3 (xE ) =
ΨE ′
[(s ⋅ c − sh ⋅ ch) sin(λE x′E ) sinh(λE x′E ) + sh2 sin(λE x′E ) cosh(λE x′E )
λE (sh − s )
2 2

−s2 cos(λE x′E ) sinh(λE x′E )] (21)

1
5 (xE ) =
ΨE [2s ⋅ sh ⋅ sin(λE x′E ) sinh(λE x′E ) − (s ⋅ ch + c ⋅ sh) sin(λE x′E ) cosh(λE x′E )
′
sh − s
2 2

+(s ⋅ ch + c ⋅ sh) cos(λE x′E ) sinh(λE x′E )] (22)

1
6 (xE ) =
ΨE [(c ⋅ sh − s ⋅ ch) sin(λE x′E ) sinh(λE x′E ) + s ⋅ sh ⋅ sin(λE x′E ) cosh(λE x′E )
′
λE (sh2 − s2 )
−s ⋅ sh ⋅ cos(λE x′E ) sinh(λE x′E )] (23)

Donde:
s = sin(λE LE ).
c = cos(λE LE ).
sh = sinh(λE LE ).
ch = cosh(λE LE ).

Funciones de Green
⎧
⎪
⎪GI (x′E , ξE′
) = Z2E (x′E )ΨE
2 (ξE ) + Z3 (xE )Ψ3 (ξE )
′ E ′ E ′
0 < x′E ≤ ξE
′
Gyy (x′E , ξE
′
) = ⎨ yy (24)
⎪
⎪GII (x′ , ξ ′ ) = Z5E (x′E )ΨE
5 (ξE ) + Z6 (x)Ψ6 (ξE )
′ E E ′
ξE
′
≤ x′E < LE
⎩ yy E E

3
Donde:
sinh(λE x′E ) cos(λE x′E ) − sin(λE x′E ) cosh(λE x′E )
Z2E (x′E ) = (25a)
4EE IE λ3E
sin(λE x′E ) sinh(λE x′E )
Z3E (x′E ) = (25b)
2EE IE λ2E
sinh[λE (LE − x′E )] cos[λE (LE − x′E )] − sin[λE (LE − x′E )] cosh[λE (LE − x′E )]
Z5E (x′E ) = (25c)
4EE IE λ3E
sin[λE (LE − x′E )] sinh[λE (LE − x′E )]
Z6E (x′E ) = − (25d)
2EE IE λ2E

Para facilitar el cálculo de los campos de desplazamiento empotrados, se recomienda emplear la siguiente
forma expandida de Z5E (x′E ) y Z6E (x′E ) en lugar de la presentada en las formulas anteriores:

(s ⋅ sh + c ⋅ ch) sin(λE x′E ) cosh(λE x′E ) − (s ⋅ ch + c ⋅ sh) sin(λE x′E ) sinh(λE x′E )
Z5E (x′E ) =
4EE IE λ3E
(s ⋅ sh − c ⋅ ch) cos(λE x′E ) sinh(λE x′E ) − (s ⋅ ch − c ⋅ sh) cos(λE x′E ) cosh(λE x′E )
+ (26)
4EE IE λ3E

c ⋅ sh ⋅ sin(λE x′E ) cosh(λE x′E ) − c ⋅ ch ⋅ sin(λE x′E ) sinh(λE x′E )

Z6E (x′E ) =
2EE IE λ2E
s ⋅ ch ⋅ cos(λE x′E ) sinh(λE x′E ) − s ⋅ sh ⋅ cos(λE x′E ) cosh(λE x′E )
+ (27)
2EE IE λ2E

Sistema de ecuaciones del método de rigidez con funciones de Green

⎧
⎪ ⎫ ⎡ ⎤ ⎧ ′E ⎫ ⎧ ′ ⎫
F Yi′E ⎪ k26 ⎥⎥ ⎪vi ⎪ ⎪∫0 Ψ2 (xE )qE (xE )dxE ⎪
LE
⎪
⎪ ⎪
⎪ ⎢k22 k23 k25 ⎪
⎪ ⎪
⎪ ⎪
⎪
E ′ ′
⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎪ ⎪ ⎪ ⎪
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ ⎥⎪⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪
⎪ i ⎪ ⎢⎢k32
⎪ ⎪ k36 ⎥⎥ ⎪ ⎪ ⎪ ⎪
L
M ′E
k33 k35 ⎪ i ⎪ ⎪∫0
θ ′E E
Ψ E
3 (x ′
)q
E E (x ′
E )dx E⎪
′
⎨ ⎬=⎢ ⎥⎨ ⎬ − ⎨ L ⎬ (28)
⎪
⎪ ′E ⎪
⎪ ⎢k ⎥
k56 ⎥ ⎪⎪ ′E ⎪
⎪ ⎪
⎪ (x )q (x )dx ′ ⎪
⎪
⎪
⎪F Y ⎪
⎪ ⎢ 52 k53 k55 ⎪v
⎪ ⎪ ⎪
⎪ ⎪
⎪∫
E
Ψ E ′ ′
⎪
⎪
⎪
⎪
j
⎪
⎪ ⎢ ⎥⎪ j
⎪ 0 5 E E E E ⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥⎪⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ k66 ⎥⎦ ⎪
⎪ ⎪
⎪ ⎪
⎪
LE E
(x )q (x )dx E⎪
⎪
M
⎩ j ⎭ ⎣k62 k63 k65 ⎩ j ⎭ ⎩∫0
θ Ψ
′E ′E ′ ′ ′
6 E E E ⎭
Donde:
s ⋅ c + sh ⋅ ch
k22 = k55 = 4EE IE λ3E (29a)
sh2 − s2
s + sh2
2
k23 = k32 = −k56 = −k65 = 2EE IE λ2E 2 (29b)
sh − s2
s ⋅ ch + c ⋅ sh
k25 = k52 = −4EE IE λ3E (29c)
sh2 − s2
s ⋅ sh
k26 = k62 = −k35 = −k53 = 4EE IE λ2E 2 (29d)
sh − s2
sh ⋅ ch − s ⋅ c
k33 = k66 = 2EE IE λE (29e)
sh2 − s2
s ⋅ ch − c ⋅ sh
k36 = k63 = 2EE IE λE (29f)
sh2 − s2
Campo de desplazamiento
vE (x′E ) = vE
h
(x′E ) + vE
f
(x′E ) (30)
Donde:
h
vE (x′E ) = ψ2E (x′E )vi′E + ψ3E (x′E )θi′E + ψ5E (x′E )vj′E + ψ6E (x′E )θj′E (31a)
x′E LE
f
vE (x′E ) = ∫ yy (xE , ξE )qE (ξE )dξE + ∫
GII ′ ′ ′ ′
GIyy (x′E , ξE
′
)qE (ξE
′
)dξE
′
(31b)
0 x′E

Pórtico plano uniforme (Euler-Bernoulli)

Ecuaciones diferenciales gobernantes

Las ecuaciones diferenciales gobernantes de este problema son (6) y (11).

4
Funciones de forma

Las funciones de forma de este elementos son ψ1E (x′E ) a ψ6E (x′E ) presentadas en las ecuaciones (8) y (13).

Funciones de Green

Las funciones de forma a emplear en este tipo de elementos son Gxx (x′E , ξE
′
) y Gyy (x′E , ξE
′
) presentadas
en las ecuaciones (9) y (14) respectivamente.

Sistema de ecuaciones del método de rigidez con funciones de Green

⎧ ⎫ ⎡ ⎤⎧ ⎫ ⎧ ⎫
⎪
⎪ ⎪
⎪ ⎢ AE EE ⎥⎪ ⎪ ⎪ ⎪
⎪
⎪ ′E ⎪
⎪ ⎢ AE EE ⎥⎪⎪ ⎪
′E ⎪
⎪
⎪ ⎪
⎪
⎪
⎪ i ⎪⎪ ⎢ − ⎥⎪⎪ ⎪
⎪ ⎪
⎪ LE E
(x )p (x )dx⎪
⎪
⎥⎪ i ⎪ ⎪ ⎪
F X 0 0 0 0 u ∫ ψ ′ ′
⎪
⎪ ⎪
⎪ ⎢ LE ⎪ ⎪ ⎪ E ⎪
⎥⎪ ⎪ ⎪ ⎪
0 1 E E
⎪ ⎪ ⎢ LE ⎪ ⎪ ⎪ ⎪
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ ⎥⎪⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥ ⎪
⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ′E ⎪
⎪
⎪ ⎢ 12EE IE 6EE IE 12EE IE 6EE IE ⎥ ⎪⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪
⎪
⎪ F Y ⎪
⎪ ⎢ 0 0 − ⎥ ⎪
⎪ v i ⎪
′E
⎪ ⎪
⎪ ∫
L E
ψ E
(x ′
)q (x′
)dx ⎪
⎪
⎪
⎪
i ⎪
⎪ ⎢ L3E L2E L3E LE ⎥ ⎪
2 ⎥ ⎪ ⎪
⎪ ⎪
⎪ 0 2 E E E ⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎪
⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥ ⎪
⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ′E ⎪ ⎪
⎪ ⎢ 6EE IE 4EE IE 6EE IE 2EE IE ⎥ ⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ M ⎪
⎪ ⎢ 0 0 − ⎥⎪⎪
⎪ θ i ⎪
′E ⎪
⎪ ⎪
⎪
⎪ ∫
LE E
ψ (x ′
)q (x′
)dx ⎪
⎪
⎪
⎪ i
⎪ ⎢ L2E L2E LE ⎥ ⎪⎥ ⎪ ⎪ 0 3 E E E
⎪
⎬=⎢
LE
⎨ ⎢ ⎥⎪⎨ ⎬− ⎨ ⎬
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ ⎥⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ′E ⎪
⎪ ⎢ AE EE AE EE ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪ ⎪ ⎢− ⎥⎪ ⎪ ⎪ ψ4E (x′E )pE (x′E )dx⎪
⎥⎪ j ⎪ ⎪ ⎪
L
⎪F Xj ⎪ ⎪ ⎪ ⎪∫ ⎪
E
⎪ ⎪ ⎢ 0 0 0 0 ⎪ u′E ⎪ ⎪ ⎪
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ LE LE ⎥⎪⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪
0 ⎪
⎪
⎪
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ ⎥⎪⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ⎪ ⎢ 12EE IE 6EE IE 12EE IE 6EE IE ⎥ ⎪ ⎪ ⎪ ⎪
⎪ F Yj′E ⎪⎪ ⎢ − − − ⎥ ⎪
⎪ ′E ⎪ ⎪
⎪ ⎪ (x )q E (xE )dx ⎪
⎪
L E
⎪ ⎪ ⎪ ⎪ ⎪ ∫ ⎪
E ′ ′
⎪ ⎪ ⎢ 0 0 ⎥ v ψ
LE ⎥ ⎪ ⎪ ⎪ ⎪
j 5
⎪ ⎪ ⎪ ⎪ ⎪ ⎪
E
⎢
0
⎪ ⎪ L3E L2E L3E 2
⎪ ⎪ ⎪ ⎪
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ ⎥⎪⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ 6EE IE 2EE IE 6EE IE 4EE IE ⎥ ⎪
⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ Mj′E ⎪ ⎪ ⎢ − ⎥⎪ ′E ⎪
⎪ θj ⎪ ⎪ ⎪⎪∫0 ψ6 (xE )qE (xE )dx⎪
⎪ LE E ′ ′
⎪
⎪ ⎭ ⎢
⎪ ⎥ ⎪ ⎪
0 0
⎩ ⎣ L2E LE L2E LE ⎦ ⎩ ⎭ ⎩ ⎭
(32)
Matriz de transformación

⎡ cos(θE ) sin(θE ) 0 0 0 0⎤
⎢ ⎥
⎢− sin(θE ) cos(θE ) 0⎥
⎢ 0 0 0 ⎥
⎢ 0⎥
⎢ 0 0 1 0 0 ⎥
[TE ] = ⎢ ⎥ (33)
⎢ 0 0 0 cos(θE ) sin(θE ) 0⎥
⎢ ⎥
⎢ − sin(θE ) cos(θE ) 0⎥
⎢ 0 0 0 ⎥
⎢ 1⎥
⎣ 0 0 0 0 0 ⎦

Campos de desplazamiento
uE (x′E ) = uhE (x′E ) + ufE (x′E ) (34)
Donde:

uhE (x′E ) = ψ1E (x′E )u′E

i + ψ4 (xE )uj
E ′ ′E
(35a)
x′E LE
ufE (x′E ) = ∫ xx (xE , ξE )pE (ξE )dξE + ∫
GII ′ ′ ′ ′
GIxx (x′E , ξE
′
)pE (ξE
′
)dξE
′
(35b)
0 x′E

vE (x′E ) = vE
h
(x′E ) + vE
f
(x′E ) (36)

h
vE (x′E ) = ψ2E (x′E )vi′E + ψ3E (x′E )θi′E + ψ5E (x′E )vj′E + ψ6E (x′E )θj′E (37a)
x′E LE
f
vE (x′E ) = ∫ yy (xE , ξE )qE (ξE )dξE + ∫
GII ′ ′ ′ ′
GIyy (x′E , ξE
′
)qE (ξE
′
)dξE
′
(37b)
0 x′E

Pila uniforme sin fricción en el fuste (Euler-Bernoulli)

Ecuaciones diferencial gobernantes

Las ecuaciones diferenciales de este elemento son (6) y (18).

Funciones de forma

Las funciones de forma de este problema son ψ1E (x′E ), ΨE

2 (xE ), Ψ3 (xE ), ψ4 (xE ), Ψ5 (xE ) y Ψ6 (xE ) presentadas
′ E ′ E ′ E ′ E ′

en las ecuaciones (8) y entre la ecuación (20) y (23).

Funciones de Green

Las funciones de Green a emplear con este elemento son Gxx (x′E , ξE
′
) presentada en la ecuación (9) y Gyy (x′E , ξE
′
)
presentada en (24).

5
Sistema de ecuaciones del método de rigidez con funciones de Green

⎧
⎪ ⎫
⎪ ⎡ A E ⎤⎧ ⎫ ⎧ ⎫
⎪ ⎢ E E ⎥⎪ ⎪ ⎪ ′ ⎪
F Xi′E ⎪ ⎪ ⎪ ⎪ ⎪
⎪ ⎪ AE EE ⎪ ⎪ ⎪
⎪ ⎪ ⎢ − 0 ⎥⎥⎪ i ⎪ ⎪ ∫0 ψ1 (xE )pE (xE )dxE ⎪ ⎪
LE E
⎪ ⎪ ⎢ LE 0 0 0 ⎪u′E ⎪ ⎪
′ ′
⎪
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ LE ⎥⎪⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ ⎥⎪⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪ ′E ⎪ ⎢ ⎥ ⎪ ′E ⎪ ⎪ ′ ⎪
⎪ ⎪ ⎪ ⎪ ⎪ Ψ2 (xE )qE (xE )dxE ⎪ ⎪
LE
⎪
⎪ F Yi ⎪ ⎪ ⎢ 0 k22 k23 0 k25 k26 ⎥ ⎪
⎪ vi ⎪ ⎪ ⎪
⎪∫
E ′ ′
⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥⎪⎪ ⎪
⎪ ⎪
⎪ 0 ⎪
⎪
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ′E ⎪⎪ ⎢ ⎥⎪⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪
⎪
⎪ ⎪ ⎢ ⎥ ⎪ ⎪ ⎪ LE E
(x )q (x )dx ⎪
⎪ Mi ⎪ ⎪ ⎢ k36 ⎥ ⎪ ⎪ ⎪∫0 E⎪
′E ′ ′ ′
0 k32 k33 0 k35 ⎪ θ i ⎪ ⎪ Ψ 3 E ⎪
⎬=⎢ ⎥⎨
E E
⎨ ⎢ ⎥⎪ ⎬ − ⎨ ⎬ (38)
⎪
⎪
⎪ ⎪
⎪ ⎢ AE EE AE EE ⎥⎪ ⎪
⎪ ⎪
⎪ ⎪
′ ⎪
⎪F Xj′E ⎪⎪ ⎢− ⎥
0 ⎥⎪ ⎪ ′E ⎪
uj ⎪ ⎪
⎪ ∫0 ψ4 (xE )pE (xE )dxE ⎪
⎪
LE E
⎪ ⎪ ⎪ ⎪ ⎪ ⎪
′ ′
⎪ ⎪ ⎢ 0 0 0 ⎪ ⎪ ⎪ ⎪
⎪
⎪
⎪ ⎪
⎪
⎪ ⎢ LE LE ⎥⎪⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪ ⎪
⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ⎪ ⎢ ⎥⎪ ⎪ ⎪ ⎪
⎪ F Yj ⎪ ⎪ ⎢ k56 ⎥ ⎪
⎪ ⎪
⎪ ⎪
⎪ (x )q (x )dxE⎪
⎪
L E
⎪ ⎪ ⎥⎪ j ⎪ ⎪∫0 ⎪
′E ′E E ′ ′ ′
⎪ ⎪ 0 k52 k53 0 k55 ⎪ v ⎪ ⎪ Ψ ⎪
⎪
⎪ ⎪
⎪ ⎢ ⎪
⎪ ⎪
⎪ ⎪
⎪
5 E E E
⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥⎪⎪ ⎪
⎪ ⎪
⎪ ⎪
⎪
⎪
⎪ ⎪
⎪ ⎢ ⎥⎪⎪ ⎪
′E ⎪ ⎪
⎪ ′ ⎪
⎪ M ′E ⎪ ⎪ ⎢ ⎥ ⎪ ⎪ ⎪ Ψ6 (xE )qE (xE )dxE ⎪ ⎪
L
⎪ k66 ⎥ ⎪ θj ⎪ ⎪ E
⎪
⎭ ⎢ ∫
E ′ ′
⎩ j ⎪
⎪ ⎣
0 k62 k63 0 k65
⎦⎩⎪ ⎪
⎭ ⎩ ⎪ 0 ⎪
⎭
Donde los valores de los términos kij se presentan en las ecuaciones (29).

Matriz de transformación

La matriz de transformación de este elemento es la presentada en la ecuación (33).

Campos de desplazamiento
uE (x′E ) = uhE (x′E ) + ufE (x′E ) (39)
Donde:

uhE (x′E ) = ψ1E (x′E )u′E

i + ψ4 (xE )uj
E ′ ′E
(40a)
x′E LE
ufE (x′E ) = ∫ xx (xE , ξE )pE (ξE )dξE + ∫
GII ′ ′ ′ ′
GIxx (x′E , ξE
′
)pE (ξE
′
)dξE
′
(40b)
0 x′E

vE (x′E ) = vE
h
(x′E ) + vE
f
(x′E ) (41)

h
vE (x′E ) = ΨE
2 (xE )vi + Ψ3 (xE )θi + Ψ5 (xE )vj + Ψ6 (xE )θj
′ ′E E ′ ′E E ′ ′E E ′ ′E
(42a)
x′E LE
f
vE (x′E ) = ∫ yy (xE , ξE )qE (ξE )dξE + ∫
GII ′ ′ ′ ′
GIyy (x′E , ξE
′
)qE (ξE
′
)dξE
′
(42b)
0 x′E

TIMOSHENKO
Viga uniforme (Timoshenko)

Ver el artı́culo [2].

Referencias
[1] Juan Camilo Molina-Villegas, Jorge Eliecer Ballesteros Ortega, and Giovanni Martı́nez Martı́nez.
Closed-form solution for non-uniform Euler–Bernoulli beams and frames. Engineering Structures,
292:116381, 2023.
[2] Juan Camilo Molina-Villegas and Jorge Eliecer Ballesteros Ortega. Closed-form solution of Timoshen-
ko frames using the green’s function stiffness method. International Journal of Solids and Structures,
269:112180, 2023.