CURSO: DINAMICA ESTRUCTURAL APLICADA A LA

INGENIERIA SISMICA CON MATHCAD PRIME

MÓDULO 2: BASES TEORICAS

Numero de barras: nb ≔ 3

Contador: i ≔ 1 ‥ nb

Longitudes de cada barra: ⎡3⎤

⎢ ⎥
l≔⎢4⎥
⎢ ⎥
⎢⎣ 3 ⎥⎦
⎡ 27 ⎤
⎢ ⎥
Rigidez a flexión de cada barra: EI ≔ ⎢ 1 ⎥
⎢ ⎥
⎢⎣ 1 ⎥⎦

La matriz de las barras 1 y 3 en el sistema local será:

⎡ 12 -6 ⎤
⎢ ― ―― ⎥
⎢ l3 l2 ⎥
⎢ i i ⎥ ⎡ 12 -18 ⎤ ⎡ 0.444 -0.667 ⎤
k ≔ EI ⋅ ⎢ ⎥ k =⎢ ⎥ k =⎢ ⎥
⎢ -6 4 ⎥ ⎢⎣ -18 36 ⎥⎦ ⎢⎣ -0.667 1.333 ⎥⎦
i i 1 3

⎢ ――2
―⎥
l
⎢ i l ⎥
i
⎣ ⎦

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 1 de 19
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

La matriz de la barra 2 en el sistema local será:

⎡ 12 -6 12 -6 ⎤
⎢ ― ―― -― ―― ⎥
⎢ l3 l2 l3 l2 ⎥
⎢ 2 2 2 2 ⎥
⎢ ⎥
⎢ -6 4 6 2 ⎥
⎢ 2―― ― ― ―⎥ ⎡ 0.188 -0.375 -0.188 -0.375 ⎤
⎢ l2 l l2 l ⎥ ⎢ ⎥
2 2 2 ⎢ -0.375 1 0.375 0.5 ⎥
k ≔ EI ⋅ ⎢ ⎥=⎢ ⎥
2 2 ⎢ 12 6 12 6 ⎥ ⎢ -0.188 0.375 0.188 0.375 ⎥
⎢ -― ― ― ―⎥ ⎢ ⎥
⎢ l
3
l2 l3 l 2 ⎥ ⎣ -0.375 0.5 0.375 1 ⎦
2 2 2 2
⎢ ⎥
⎢ -6 2 6 4 ⎥
⎢ ――― ― ―⎥
⎢ l2 l l2 l ⎥
⎢⎣ 2 2 2 2 ⎥
⎦
Como convertimos estas matrices al sistema global?

⎡ D1 D2 D3 ⎤
⎣ ⎦
⎡ D1 D2 D3 ⎤ ⎡ D1 D2 D3 ⎤
⎣ ⎦ ⎡ 0 -0.75 0 ⎤ ⎣ ⎦
⎢ ⎥
⎡1 0 0⎤ ⎢0 1 0⎥ ⎡1 0 0⎤
a ≔⎢ ⎥ a ≔⎢ ⎥ a ≔⎢ ⎥
1 ⎢⎣ 0 1 0 ⎥⎦ 2
⎢0 0 0⎥ 3 ⎢⎣ 0 0 1 ⎥⎦
⎢0 0 1 ⎥⎦
⎣

Matriz de rigidez del pórtico:

⎡ 12.444 -18 -0.667 ⎤
nb ⎢ ⎥
K ≔ ∑ ⎛a ⋅ ⎛k ⎞ ⋅ a ⎞ = ⎢ -18
T
37.668 0.781 ⎥
⎜ i ⎜⎝ i⎟⎠ i⎟⎠ ⎢
i=1 ⎝ ⎥
⎢⎣ -0.667 0.781 2.333 ⎥⎦

Desplazamientos en el sistema Q-D:

⎡ 20 ⎤ ⎡ 5.255 ⎤
⎢ ⎥ ⎢ ⎥
Q≔⎢ 0 ⎥ D ≔ K -1 ⋅ Q = ⎢ 2.497 ⎥
⎢ ⎥ ⎢ ⎥
⎢⎣ 0 ⎥⎦ ⎢⎣ 0.665 ⎥⎦
Desplazamientos de cada barra:
⎡ -1.873 ⎤
⎢ ⎥
⎡ 5.255 ⎤ ⎢ 2.497 ⎥ ⎡ 5.255 ⎤
d ≔a ⋅D d =⎢ ⎥ d =⎢ ⎥ d =⎢ ⎥
i i 1 ⎢⎣ 2.497 ⎥⎦ 2
⎢ 0 ⎥
3 ⎢⎣ 0.665 ⎥⎦
⎢ 0.665 ⎥
⎣ ⎦

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 2 de 19
 -1.873
⎢ ⎥
⎢ 2.497
CURSO: DINAMICA ⎥
ESTRUCTURAL APLICADA A LA
d =⎢ ⎥
INGENIERIA2 SISMICA
⎢ 0 CON
⎥ MATHCAD PRIME
⎢ 0.665
MÓDULO 2: BASES ⎥
⎣ TEORICAS
⎦

Fuerzas en los elementos:

⎡ -1.537 ⎤
⎢ ⎥
⎡ 18.108 ⎤ ⎢ 3.532 ⎥ ⎡ 1.892 ⎤
q ≔k ⋅d q =⎢ ⎥ q =⎢ ⎥ q =⎢ ⎥
⎢ 1.537 ⎥
i i i 1 ⎢⎣ -4.685 ⎥⎦ 2 3 ⎢⎣ -2.616 ⎥⎦
⎢ 2.616 ⎥
⎣ ⎦

Desplazamiento horizontal: Δh ≔ ⎛d ⎞ = 5.255

⎜⎝ 1⎟⎠
1

Fuerza aplicada: F ≔ Q = 20
1

F
Rigidez: kl ≔ ―= 3.806 EI
Δh

Usando condensación estática: K0 ≔ K = 12.444

1,1

K1 ≔ submatrix ⎛⎝K , 1 , 1 , 2 , 3⎞⎠ = ⎡⎣ -18 -0.667 ⎤⎦

⎡ -18 ⎤
⎛ ⎞
K2 ≔ submatrix ⎝K , 2 , 3 , 1 , 1⎠ = ⎢ ⎥
⎢⎣ -0.667 ⎥⎦

⎡ 37.668 0.781 ⎤
K3 ≔ submatrix ⎛⎝K , 2 , 3 , 2 , 3⎞⎠ = ⎢ ⎥
⎢⎣ 0.781 2.333 ⎥⎦

Rigidez Lateral del Pórtico: kl ≔ ⎛⎝K0 - K1 ⋅ K3 -1 ⋅ K2⎞⎠ = 3.806 EI

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 3 de 19
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

Calcule la matriz de rigidez y masa de la estructura mostrada, considerando un peso de

15 ton en cada piso intermedio:
Datos:
tonnef
Modulo de Elasticidad: E ≔ 2.2 ⋅ 10 6 ―――
m2

Alturas de entrepiso: H1 ≔ 5 m H2 ≔ 3.6 m H3 ≔ 3.6 m

⎡H ⎤ ⎡5 ⎤
⎢ 1⎥ ⎢ ⎥
h ≔ ⎢ H2 ⎥ = ⎢ 3.6 ⎥ m
⎢ ⎥ ⎢ ⎥
⎢⎣ H3 ⎥⎦ ⎢⎣ 3.6 ⎥⎦

Numero de pisos: n ≔ rows ⎛⎝h⎞⎠ = 3 i≔1‥n

Pesos: W1 ≔ 65 tonnef W2 ≔ 75 tonnef W3 ≔ 70 tonnef

⎡ W ⎤ ⎡ 65 ⎤
⎢ 1⎥ ⎢ ⎥ m
W ≔ ⎢ W2 ⎥ = ⎢ 75 ⎥ tonnef g ≔ 9.81 ―
⎢ ⎥ ⎢ ⎥ s2
⎢⎣ W3 ⎥⎦ ⎢⎣ 70 ⎥⎦

⎡ 6.626 ⎤
W ⎢ ⎥ s2
Masas: m' ≔ ―= ⎢ 7.645 ⎥ tonnef ⋅ ―
g ⎢ ⎥ m
⎢⎣ 7.136 ⎥⎦
Matriz de Masas:
⎡ m' 0 0 ⎤
⎢ 1 ⎥ ⎡ 6.626 0 0 ⎤
⎢ ⎥ ⎢ ⎥ s2
M ≔ ⎢ 0 m' 0 ⎥ = ⎢ 0 7.645 0 ⎥ tonnef ⋅ ―
⎢
2
⎥ ⎢0 ⎥ m
0 7.136 ⎥⎦
⎢ 0 0 m' ⎥ ⎢⎣
⎢⎣ 3⎥
⎦

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción totalcm

b1 ≔ 120 o parcial h
sin
1≔ 120 cm
consentimiento del autor. Página 4 de 19
 ⎢ 1 ⎥ 6.626 0 0
⎢ ⎥ ⎢ ⎥ s2
M ≔ CURSO:
⎢ 0 m' DINAMICA
0 ⎥ = ⎢ 0 ESTRUCTURAL
7.645 0 tonnef ⋅ ―A LA
⎥ APLICADA
PRIMEm
2 ⎢ ⎥
⎢INGENIERIA ⎥SISMICA
0 CON
0 MATHCAD
7.136 ⎥⎦
⎢ 0 0 m' ⎥ ⎢⎣
MÓDULO
⎢⎣ 2: BASES
3⎥
⎦ TEORICAS

Sección Rectangular: b1 ≔ 120 cm h1 ≔ 120 cm

b1 ⋅ h1 3
Inercia: I1 ≔ ――― = 17280000 cm 4
12

E ⋅ I1 = 380160 tonnef ⋅ m 2

H1 3 m
Matriz de flexibilidad: f11 ≔ ――― = ⎛⎝1.096 ⋅ 10 -4⎞⎠ ―――
3 ⋅ E ⋅ I1 tonnef
H1 3 H1 2 m
f21 ≔ ――― + ――― ⋅ H2 = ⎛⎝2.28 ⋅ 10 -4⎞⎠ ―――
3 ⋅ E ⋅ I1 2 ⋅ E ⋅ I1 tonnef

H1 3 H1 2 m
f31 ≔ ――― + ――― ⋅ ⎛⎝H2 + H3⎞⎠ = ⎛⎝3.463 ⋅ 10 -4⎞⎠ ―――
3 ⋅ E ⋅ I1 2 ⋅ E ⋅ I1 tonnef
⎛H + H ⎞ 3 m
⎝ 1 2⎠
f22 ≔ ―――― = ⎛⎝5.577 ⋅ 10 -4⎞⎠ ―――
3 ⋅ E ⋅ I1 tonnef
3 2
⎛H + H ⎞ ⎛H + H ⎞ m
⎝ 1 2⎠ ⎝ 1 2⎠
f32 ≔ ―――― + ―――― ⋅ H3 = ⎛⎝9.079 ⋅ 10 -4⎞⎠ ―――
3 ⋅ E ⋅ I1 2 ⋅ E ⋅ I1 tonnef

⎛H + H + H ⎞ 3 m
⎝ 1 2 3⎠
f33 ≔ ―――――= ⎛⎝1.592 ⋅ 10 -3⎞⎠ ―――
3 ⋅ E ⋅ I1 tonnef

⎡ f f f ⎤ ⎡ 1.096 ⋅ 10 -4 2.28 ⋅ 10 -4 3.463 ⋅ 10 -4 ⎤

⎢ 11 21 31 ⎥ ⎢ ⎥ m
Ff ≔ ⎢ f21 f22 f32 ⎥ = ⎢ 2.28 ⋅ 10 -4 5.577 ⋅ 10 -4 9.079 ⋅ 10 -4 ⎥ ―――
⎢ ⎥ ⎢ ⎥ tonnef
⎢⎣ f31 f32 f33 ⎥⎦ ⎢⎣ 3.463 ⋅ 10 9.079 ⋅ 10 0.002
-4 -4
⎥⎦

Matriz de Rigidez:
⎡ 9.058 ⋅ 10 4 -6.902 ⋅ 10 4 1.965 ⋅ 10 4 ⎤
⎢ ⎥ tonnef
K ≔ Ff = ⎢ -6.902 ⋅ 10 4 7.758 ⋅ 10 4 -2.923 ⋅ 10 4
-1
⎥ ―――
⎢ 4 4 4
⎥ m
⎢⎣ 1.965 ⋅ 10 -2.923 ⋅ 10 1.302 ⋅ 10 ⎥⎦

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 5 de 19
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

Sea el siguiente campo de desplazamiento admisible en un elemento viga:

Es posible demostrar imponiendo las condiciones de frontera con los parámetros nodales
que las funciones de forma para este elemento son:

Donde 1 y 3 están relacionados a los desplazamientos verticales en los extremos y 2 y 4 a

las rotaciones.
2
3 ⋅ x2 2 ⋅ x3 ⎛ x⎞
N1 ⎛⎝x , i⎞⎠ ≔ 1 - ――+ ―― ⎛ ⎞
N2 ⎝x , i⎠ ≔ x ⋅ ⎜1 - ― ⎟
l2 l3 ⎜ l⎟
i i ⎜⎝ i⎟
⎠
3 ⋅ x2 2 ⋅ x3 x 2 ⎛ x ⎞
⎛ ⎞
N3 ⎝x , i⎠ ≔ ――- ―― ⎛ ⎞
N4 ⎝x , i⎠ ≔ ―⋅ ― ⎜ - 1⎟
l2 l3 l ⎜ l ⎟
i i i ⎜⎝ i ⎟⎠

x ≔sin
Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial 0 , 0.1 ‥l
consentimiento del autor. Página 6 de 19
1
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

Por ejemplo para el elemento 1 del primer ejemplo: x ≔ 0 , 0.1 ‥ l

1

Como vimos en el método de los elementos finitos estas funciones de forma nos permite
obtener la matriz de rigidez de un elemento y las fuerzas nodales equivalentes:

Con estas expresiones podemos deducir la misma matriz de rigidez del elemento viga
deducida según la teoría del análisis estructural aplicando desplazamientos unitarios y
calculando las fuerzas que producen dichos desplazamientos.

Porque obtenemos los mismo valores de la matriz de rigidez?

Podemos usar las mismas funciones de forma para interpolar los

desplazamientos internos para cualquier condición de carga de forma
exacta?

Incluyendo las deformación axiales, se deduce la matriz de masa para un elemento columna
en 2D:

x
N0 ⎛⎝x , i⎞⎠ ≔ 1 - ―
l
i
x
N5 ⎛⎝x , i⎞⎠ ≔ ―
l
i

Por ejemplo para la columna del ejemplo, suponiendo que es de concreto y de sección 30x50cm
Página 7 de 19
Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor.
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

Por ejemplo para la columna del ejemplo, suponiendo que es de concreto y de sección 30x50cm

Masa del elemento: me ≔ ⎛⎝0.3 ⋅ 0.5 ⋅ 2.4⎞⎠ ⋅ l = 1.08 ton

1 1
me
1 ton
Masa distribuida por unidad de longitud: m'e ≔ ―― = 0.36 ――
1
l m2
1
⎡ x ⎤
⎢ 1-― ⎥
⎢ l ⎥
⎢ 1 ⎥
⎢ ⎥
⎢ 3 ⋅ x2 2 ⋅ x3 ⎥
1 - ――+ ――
⎢ ⎥
⎢ l2 l3 ⎥
1 1
⎢ ⎥
2
⎢ ⎛ x ⎞ ⎥
⎢ x ⋅ ⎜1 - ― ⎟ ⎥ l
⎢ ⎜ l ⎟ ⎥ 1

⎛ ⎞ ⎢ ⎜⎝ 1⎟ ⎥ ⌠
N ⎝x⎠ ≔ ⎢ ⎠
⎥ M11 ≔ ⎮ m'e ⋅ ⎛⎝N ⎛⎝x⎞⎠⎞⎠ ⋅ ⎛⎝N ⎛⎝x⎞⎠⎞⎠ d x = 0.36
⎢ x ⎥ ⎮
⌡
1 1 1
―
⎢ ⎥ 0
l
⎢ 1 ⎥ l
1
⎢ ⎥
⎢ 3 ⋅ x2 2 ⋅ x3 ⎥ ⌠
M22 ≔ ⎮ m'e ⋅ ⎛⎝N ⎛⎝x⎞⎠⎞⎠ ⋅ ⎛⎝N ⎛⎝x⎞⎠⎞⎠ d x = 0.401
⎢ ―― - ―― ⎥ 1 2 2
l2 l3 ⎮
⌡
⎢ ⎥
1 1 0
⎢ ⎥ l
1
⎢ x2 ⎛ x ⎞ ⎥
⌠
⎢ ―⋅ ⎜―- 1⎟ ⎥ M33 ≔ ⎮ m'e ⋅ ⎛⎝N ⎛⎝x⎞⎠⎞⎠ ⋅ ⎛⎝N ⎛⎝x⎞⎠⎞⎠ d x = 0.093
⎢ l ⎜l ⎟ ⎥
⎮
⌡
1 3 3
⎢⎣ 1 ⎜⎝ 1 ⎟⎠ ⎥⎦
0

⎡ 140 0 0 70 0 ⎤ 0
⎢ ⎥
⎢ 0 156 22 ⋅ l 0 54 -13 ⋅ l ⎥
⎢ 1 1⎥

me ⎢ 0 22 ⋅ l 4⋅l 2
0 13 ⋅ l -3 ⋅ l 2 ⎥
1 ⎢ 1 1 1 1 ⎥
Mele ≔ ―― ⋅⎢ ⎥
420 ⎢ 70 0 0 140 0 0 ⎥
⎢ ⎥
⎢ 0 54 13 ⋅ l 0 156 -22 ⋅ l ⎥
1 1
⎢ ⎥
2
⎢ 0 -13 ⋅ l -3 ⋅ l 0 -22 ⋅ l 4 ⋅ l 2 ⎥
⎢⎣ 1 1 1 1 ⎥⎦

⎡ 0.36 0 0 0.18 0 ⎤ 0
⎢ ⎥
⎢0 0.401 0.17 0 0.139 -0.1 ⎥
⎢ ⎥
0 0.17 0.093 0 0.1 -0.069 ⎥
Mele = ⎢
⎢ 0.18 0 0 0.36 0 0 ⎥
⎢ ⎥
⎢0 0.139 0.1 0 0.401 -0.17 ⎥
⎢ ⎥
⎢⎣ 0 -0.1 -0.069 0 -0.17 0.093 ⎥⎦

La matriz de masa en el sistema global se puede ensamblar similar a la matriz de rigidez, ahora
que conocemos la masa y rigidez del sistema es posible estudiar el comportamiento dinámico.

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 8 de 19
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

Datos:

kgf
Resistencia del Concreto: f'c ≔ 320 ――
cm 2

Dimensiones de las columnas:

Cuadrada: b ≔ 60 cm h ≔ 60 cm

Circular: D ≔ 25 cm

Alturas de entrepiso: h1 ≔ 4 m h2 ≔ 3 m

Dimensiones totales en planta: B ≔ 12 m L ≔ 12 m

tonnef tonnef
Cargas en primer piso: QD1 ≔ 0.5 ――― QL1 ≔ 0.4 ―――
m2 m2
tonnef tonnef
Cargas en segundo piso: QD2 ≔ 0.4 ――― QL2 ≔ 0.15 ―――
m2 m2

Porcentaje de carga viva: λ ≔ 0.5

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 9 de 19
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

Desarrollo:

Modulo de Elasticidad del Concreto según E-060:

‾‾‾‾‾‾‾kgf kgf
Ec ≔ 15000 ⋅ f'c ⋅ ――= ⎛⎝2.683 ⋅ 10 5 ⎞⎠ ――
cm 2 cm 2

Inercia de las columnas:

Cuadrada: h3 ⎛
Icuad ≔ b ⋅ ― = ⎝1.08 ⋅ 10 6 ⎞⎠ cm 4
12
⎛0.5 ⋅ D⎞ 4
Circular: ⎝ ⎠
Icir ≔ π ⋅ ―――= ⎛⎝1.917 ⋅ 10 4 ⎞⎠ cm 4
4

Rigidez lateral de las columnas:

Biemprotradas:

Cuadradas:
Ec ⋅ Icuad tonnef
1 nivel: k1cu ≔ 12 ⋅ ――― = ⎛⎝5.434 ⋅ 10 3 ⎞⎠ ―――
h1 3 m
Ec ⋅ Icuad tonnef
2 nivel: k2cu ≔ 12 ⋅ ――― = ⎛⎝1.288 ⋅ 10 4 ⎞⎠ ―――
h2 3 m
Circulares:
Ec ⋅ Icir tonnef
1 nivel: k1ci ≔ 12 ⋅ ――― = 96.471 ―――
h1 3 m
Ec ⋅ Icir tonnef
2 nivel: k2ci ≔ 12 ⋅ ――― = 228.672 ―――
h2 3 m

Articuladas en primer nivel:

Cuadradas: Circulares:
Ec ⋅ Icuad tonnef Ec ⋅ Icir tonnef
k'1cu ≔ 3 ⋅ ――― = ⎛⎝1.358 ⋅ 10 3 ⎞⎠ ――― k'1ci ≔ 3 ⋅ ――― = 24.118 ―――
h1 3 m h1 3 m

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 10 de 19
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
Ec ⋅ Icuad tonnef Ec ⋅ Icir tonnef
k'1cu ≔ 3 ⋅ ――― = ⎛⎝1.358 ⋅ 10INGENIERIA
3⎞
⎠ ――― SISMICA k'1ciCON
≔ 3 ⋅MATHCAD
――― PRIME
= 24.118 ―――
h1 3 MÓDULOm 2: BASES TEORICAS h1 3 m

Numero de columnas:

Cuadradas Empotradas: N1 ≔ 0

Cuadradas Apoyadas: N2 ≔ 4

Circulares Empotradas: N3 ≔ 5

Circulares Apoyadas: N4 ≔ 0

Rigidez de entrepiso:

tonnef
Primer Nivel: K1 ≔ N1 ⋅ k1cu + N2 ⋅ k'1cu + N3 ⋅ k1ci + N4 ⋅ k'1ci = 5916.001 ―――
m
tonnef
Segundo Nivel: K2 ≔ N1 ⋅ k2cu + N2 ⋅ k2cu + N3 ⋅ k2ci + N4 ⋅ k2ci = 52662.368 ―――
m

Rigidez Lateral del Edificio:

⎡ K + K -K ⎤ ⎡ 58578.369 -52662.368 ⎤ tonnef

1 2 2
K≔⎢ ⎥=⎢ ⎥ ―――
⎢ -K2
⎣ K ⎥
2 ⎦ ⎣⎢ -52662.368 52662.368 ⎥⎦ m

Peso Sísmico de entrepiso:

Área en Planta: A ≔ B ⋅ L = 144 m 2

Primer Nivel: P1 ≔ A ⋅ ⎛⎝QD1 + λ ⋅ QL1⎞⎠ = 100.8 tonnef

Segundo Nivel: P2 ≔ A ⋅ ⎛⎝QD2 + λ ⋅ QL2⎞⎠ = 68.4 tonnef

Masas de entrepiso:

P1 s2
Primer Nivel: M1 ≔ ―= 10.275 tonnef ⋅ ―
g m
P2 s2
Segundo Nivel: M2 ≔ ―= 6.972 tonnef ⋅ ―
g m

Matriz de Masa del Edificio:

⎡ M 0 ⎤ ⎡ 10.275 0 ⎤ s2
1
M≔⎢ ⎥=⎢ ⎥ tonnef ⋅ ―
⎢ 0 M2 ⎥ ⎢⎣ 0 6.972 ⎥⎦ m
⎣ ⎦

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 11 de 19
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

Datos:
tonnef
Modulo de Elasticidad: E ≔ 2.2 ⋅ 10 6 ―――
m2
Alturas de entrepiso: h1 ≔ 4 m h2 ≔ 3.5 m h3 ≔ 3 m
⎡h ⎤ ⎡4 ⎤
⎢ 1⎥ ⎢ ⎥
h ≔ ⎢ h2 ⎥ = ⎢ 3.5 ⎥ m
⎢ ⎥ ⎢ ⎥
⎢⎣ h3 ⎥⎦ ⎢⎣ 3 ⎥⎦

Numero de pisos: n ≔ rows ⎛⎝h⎞⎠ = 3

Pesos: W1 ≔ 300 tonnef W2 ≔ 300 tonnef W3 ≔ 180 tonnef

⎡ W ⎤ ⎡ 300 ⎤
⎢ 1⎥ ⎢ ⎥ m
W ≔ ⎢ W2 ⎥ = ⎢ 300 ⎥ tonnef g ≔ 9.81 ―
⎢ ⎥ ⎢ ⎥ s2
⎢⎣ W3 ⎥⎦ ⎢⎣ 180 ⎥⎦

⎡ 30.581 ⎤
W ⎢ ⎥ s2
Masas: m' ≔ ―= ⎢ 30.581 ⎥ tonnef ⋅ ―
g ⎢ ⎥ m
⎢⎣ 18.349 ⎥⎦
Sección Rectangular: b1 ≔ 50 cm h1 ≔ 50 cm

b1 ⋅ h1 3
Inercia: I1 ≔ ――― = 520833.3333 cm 4
12

E ⋅ I1 = 11458.333 tonnef ⋅ m 2

Numero de columnas: N1 ≔ 12 N2 ≔ 12 N3 ≔ 8
⎡ N ⎤ ⎡ 12 ⎤
⎢ 1⎥ ⎢ ⎥
N ≔ ⎢ N2 ⎥ = ⎢ 12 ⎥
⎢ ⎥ ⎢ ⎥
⎢⎣ N3 ⎥⎦ ⎢⎣ 8 ⎥⎦

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 12 de 19
 CURSO: DINAMICA ⎡ESTRUCTURAL
N1 ⎤ ⎡ 12 ⎤ APLICADA A LA
⎢ ⎥ ⎢ ⎥
INGENIERIA SISMICA
N ≔ ⎢ N2CON MATHCAD
⎥ = ⎢ 12 ⎥ PRIME
⎢ ⎥ ⎢ ⎥
MÓDULO 2: BASES TEORICAS
8 ⎢⎣ N3 ⎥⎦ ⎢⎣ ⎥⎦

⎡ 2148.4375 ⎤
12 ⋅ E ⋅ I1 ⎢ ⎥ tonnef
Rigidez de columnas: (Bi-empotrado) k ≔ ―――= ⎢ 3206.9971 ⎥ ―――
h3 ⎢ ⎥ m
⎢⎣ 5092.5926 ⎥⎦

i≔1‥n

⎛ tonnef ⎞
f ⎛⎝x , y⎞⎠ ≔ 0 ⎜――― ⋅ s2 ⎟ M ≔ matrix ⎛⎝n , n , f⎞⎠
⎜ m ⎟
⎝ ⎠

M ≔ m'
i,i i

Rigidez de entrepisos: K' ≔ k ⋅ N

i i i
⎡ 25781.25 ⎤
⎢ ⎥ tonnef
K' = ⎢ 38483.965 ⎥ ―――
⎢ ⎥ m
⎢⎣ 40740.741 ⎥⎦
a) Matriz de rigidez y masa de la estructura:

tonnef
j≔1‥n-1 g ⎛⎝x , y⎞⎠ ≔ 0 ―――
m
KE ≔ matrix ⎝n , n , g⎠ KE ≔ if ⎛i < n , K' + K' , K' ⎞
⎛ ⎞
i,i ⎜⎝ i i+1 i⎟
⎠
KE ≔ -K' KE ≔ -K'
j,j+1 j+1 j+1,j j+1

⎡ 64265.215 -38483.965 0 ⎤
⎢ ⎥ tonnef
KE = ⎢ -38483.965 79224.706 -40740.741 ⎥ ―――
⎢ ⎥ m
⎢⎣ 0 -40740.741 40740.741 ⎥⎦
⎡ 30.581 0 0 ⎤
⎢ ⎥ tonnef
M=⎢ 0 30.581 0 ⎥ ――― ⋅ s2
⎢ ⎥ m
⎢⎣ 0 0 18.349 ⎦⎥

n
tonnef 2
Mtotal ≔ ∑ m' = 79.511 ――― ⋅s
i=1
i
m

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 13 de 19
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

Datos:
tonnef
Modulo de Elasticidad: E ≔ 2.2 ⋅ 10 6 ―――
m2

Alturas de entrepiso: h1 ≔ 3 m h2 ≔ 3 m h3 ≔ 3 m
⎡h ⎤ ⎡3⎤
⎢ 1⎥ ⎢ ⎥
h ≔ ⎢ h2 ⎥ = ⎢ 3 ⎥ m
⎢ ⎥ ⎢ ⎥
⎢⎣ h3 ⎥⎦ ⎢⎣ 3 ⎥⎦

Numero de pisos: n ≔ rows ⎛⎝h⎞⎠ = 3

Pesos: W1 ≔ 80 tonnef W2 ≔ 80 tonnef W3 ≔ 30 tonnef

⎡ W ⎤ ⎡ 80 ⎤
⎢ 1⎥ ⎢ ⎥ m
W ≔ ⎢ W2 ⎥ = ⎢ 80 ⎥ tonnef g ≔ 9.81 ―
⎢ ⎥ ⎢ ⎥ s2
⎢⎣ W3 ⎥⎦ ⎢⎣ 30 ⎥⎦
⎡ 8.155 ⎤
W ⎢ ⎥ s2
Masas: m' ≔ ―= ⎢ 8.155 ⎥ tonnef ⋅ ―
g ⎢ ⎥ m
⎢⎣ 3.058 ⎥⎦
Sección Rectangular: b1 ≔ 40 cm h1 ≔ 40 cm

b1 ⋅ h1 3
Inercia: I1 ≔ ――― = 213333.3333 cm 4
12
E ⋅ I1 = 4693.333 tonnef ⋅ m 2

Numero de columnas: N1 ≔ 12 N2 ≔ 12 N3 ≔ 6
⎡ N ⎤ ⎡ 12 ⎤
⎢ 1⎥ ⎢ ⎥
N ≔ ⎢ N2 ⎥ = ⎢ 12 ⎥
⎢ ⎥ ⎢ ⎥
⎢⎣ N3 ⎥⎦ ⎢⎣ 6 ⎥⎦
⎡ 2085.9259 ⎤
12 ⋅ E ⋅ I1 ⎢ ⎥ tonnef
Rigidez de columnas: (Bi-empotrado) k ≔ ―――= ⎢ 2085.9259 ⎥ ―――
h3 ⎢ ⎥ m
⎢⎣ 2085.9259 ⎥⎦
⎡ 521.4815 ⎤
Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total
3 ⋅ Eo⋅ Iparcial
⎢ sin consentimiento del autor. Página 14 de 19
1
⎥ tonnef
 ⎡ 2085.9259
CURSO: DINAMICA ESTRUCTURAL ⎤
APLICADA A LA
INGENIERIA SISMICA
12 ⋅ E ⋅ I
CON1 MATHCAD tonnef
⎢ ⎥
k ≔ ――― = ⎢ 2085.9259PRIME
⎥ ―――
h3
MÓDULO 2: BASES TEORICAS ⎢ ⎥ m
2085.9259 ⎢⎣ ⎥⎦
⎡ 521.4815 ⎤
3 ⋅ E ⋅ I1 ⎢ ⎥ tonnef
Rigidez de columnas: (Articulado) k' ≔ ――― = ⎢ 521.4815 ⎥ ―――
h3 ⎢ ⎥ m
⎢⎣ 521.4815 ⎥⎦

i≔1‥n ii ≔ 2 ‥ n

⎛ tonnef ⎞
f ⎛⎝x , y⎞⎠ ≔ 0 ⎜――― ⋅ s2 ⎟ M ≔ matrix ⎛⎝n , n , f⎞⎠
⎜ m ⎟
⎝ ⎠

M ≔ m'
i,i i

Rigidez de entrepisos: K' ≔ k ⋅ N K' ≔ k' ⋅ N

ii ii ii 1 1 1

⎡ 58578.369 -52662.368 ⎤ tonnef

K=⎢ ⎥ ―――
⎢⎣ -52662.368 52662.368 ⎥⎦ m

a) Matriz de rigidez y masa de la estructura:

tonnef
j≔1‥n-1 g ⎛⎝x , y⎞⎠ ≔ 0 ―――
m

KE ≔ matrix ⎛⎝n , n , g⎞⎠ KE ≔ if ⎛i < n , K' + K' , K' ⎞

i,i ⎜⎝ i i+1 i⎟
⎠
KE ≔ -K' KE ≔ -K'
j,j+1 j+1 j+1,j j+1

⎡ 31288.889 -25031.111 0 ⎤
⎢ ⎥ tonnef
KE = ⎢ -25031.111 37546.667 -12515.556 ⎥ ―――
⎢ ⎥ m
⎢⎣ 0 -12515.556 12515.556 ⎥⎦

⎡ 8.155 0 0 ⎤
⎢ ⎥ tonnef
M=⎢0 8.155 0 ⎥ ――― ⋅ s2
⎢ ⎥ m
⎢⎣ 0 0 3.058 ⎥⎦

n
tonnef 2
Mtotal ≔ ∑ m' = 19.368 ――― ⋅s
i=1
i
m

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 15 de 19
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

COORDENADAS DE LOS VERTICES EN LAS DIRECCION DE LOS PORTICOS PRINCIPALES:

PLANTAS:

PRIMERA SEGUNDA PORTICOS PORTICOS

PLANTA TIPICA PLANTA TIPICA PRIMERA PLANTA SEGUNDA PLANTA

Vertice1 X1 Y1 Vertice2 X2 Y2 Portico inicio1 fin1 Portico inicio2 fin2

⎛ m⎞ ⎛ m⎞ ⎛ m⎞ ⎛ m⎞
⎝ ⎠ ⎝ ⎠ ⎝ ⎠ ⎝ ⎠

1 0 0 1 0 0 1 1 3 1 1 2

2 7 0 2 7 0 2 4 6 2 3 4

3 16 0 3 0 6 3 7 8 3 5 6

4 0 6 4 7 6 4 1 7 4 1 5

5 7 6 5 0 12 5 2 8 5 2 6

6 11.5 6 6 7 12 6 3 8

7 0 12

8 7 12

Alturas de entrepiso y asignación de plantas:

Piso Planta hpiso

⎛ m⎞
⎝ ⎠
1 1 4 Numero de pisos: Npisos ≔ rows ⎛⎝Piso⎞⎠ = 2 i ≔ 1 ‥ Npisos

2 2 3

COORDENADAS DE PORTICOS:
Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 16 de 19
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

COORDENADAS DE PORTICOS:
PRIMERA PLANTA SEGUNDA PLANTA
Numero de Pórticos: Nporticos1 ≔ rows ⎛⎝inicio1⎞⎠ = 6 Nporticos2 ≔ rows ⎛⎝inicio2⎞⎠ = 5
‖ | ‖ |
Coordenadas de los pórticos: X1c ≔ ‖ for j ∊ 1 ‥ Nporticos1 | | X2c ≔ ‖ for j ∊ 1 ‥ Nporticos2 | |
| |
‖ ‖ || ‖ ‖ ||
‖ ‖ X1c ← X1 | ‖ ‖ X2c ← X2 |
1,j inicio1 | 1,j inicio2 |
‖ ‖ j |
| ‖ ‖ j |
|
‖ ‖ || ‖ ‖ ||
‖ ‖ X1c2 , j ← X1fin | | ‖ ‖ X2c2 , j ← X2fin | |
1 2
‖ ‖‖ j
|| ‖ ‖‖ j
||
‖ | ‖ |
‖‖ X1c || ‖‖ X2c ||

‖ | ‖ |
Y1c ≔ ‖ for j ∊ 1 ‥ Nporticos1 | | Y2c ≔ ‖ for j ∊ 1 ‥ Nporticos2 | |
| |
‖ ‖ || ‖ ‖ ||
‖ ‖ Y1c ← Y1 | ‖ ‖ Y2c ← Y2 |
1,j inicio1 | 1,j inicio2 |
‖ ‖ j |
| ‖ ‖ j |
|
‖ ‖ || ‖ ‖ ||
‖ ‖ Y 1c ← Y 1 ‖ ‖ Y 2c ← Y 2
2,j fin1 || 2,j fin2 ||
‖ ‖‖ j
|| ‖ ‖‖ j
||
‖ | ‖ |
‖‖ Y1c || ‖‖ Y1c ||

Calculo de las matrices de compatibilidad de los pórticos con respecto al origen:

‖ |
SC ≔ ‖ SC ← matrix ⎛⎝Nporticos1 , 3 , f⎞⎠ |
‖ | |
‖ for j ∊ 1 ‥ Nporticos1 ||
‖ ‖ ||
‖ ‖ xij ← X1inicio ||
1
‖ ‖ j
||
‖ ‖x ←X ||
‖ ‖ fj 1
fin1 ||
‖ ‖ j
||
‖ ‖y ←Y ||
i 1
‖ ‖ j inicio1 ||
‖ ‖ j
||
‖ ‖ yf ← Y 1 ||
‖ ‖ j fin1
j ||
‖ ‖ ||
‖ ‖ ‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾
⎛x - x ⎞ 2 + ⎛y - y ⎞ 2 ||
d ←
‖ ‖ ⎜ fj ij⎟
⎝ ⎠
⎜ fj
⎝
i
j⎟
⎠ || i' ≔ 1 ‥ Nporticos1
‖ ‖ ||
‖ ‖ yf - yi ||
‖ ‖ j j
|| COS ≔ identity ⎛⎝Npisos⎞⎠ ⋅ SC
‖ ‖ SC ← ――― || i' i' , 2
‖ ‖ j,1 d ||
‖ ‖ || SEN ≔ identity ⎛⎝Npisos⎞⎠ ⋅ SC
‖ ‖ xf - xi ||
i' i' , 1
j j
‖ ‖ SC ← ||
‖ ‖ j,2 ―――
|| Rpor ≔ SC
‖ ‖ d i' i' , 3
||
‖ ‖ SC ← identity ⎛N ⎞
‖ ‖ j,3 ⎝ pisos⎠ || aport ≔ augment ⎛COS , SEN , Rpor ⎞
|| i' ⎜ i' i' i'⎟
‖ ‖ ⎝ ⎠
|||
‖ ‖ for k ∊ 1 ‥ Npisos |||
‖ ‖ ‖
‖ ‖ ‖ ⎛⎛0 m - y ⎞ ⋅ SC - ⎛0 m - x ⎞ ⋅ SC ⎞ | | | NGLglob ≔ Npisos ⋅ 3 = 6
⎜⎜ i
j⎟ j,2 ⎜ i
j⎟ j , 1⎟ |
||
‖ ‖ ‖⎛ ⎞ ⎝⎝ ⎠ ⎝ ⎠ ⎠ |||
‖ ‖ ‖ ⎜SCj , 3⎟ ← ―――――――――――― |||
‖ ‖ ‖⎝ ⎠k , k m |
‖ ‖ ‖ || || |
‖ |
‖ SC |

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial

⎡
1 0 0sin0consentimiento
0 0 del autor.
⎤ Página 17 de 19
aport =⎢ ⎥
 ‖ ‖ ‖ 0 m - yi ⋅ SC - 0 m - xi ⋅ SC ||
⎜⎜ j⎟ j,2 ⎜ j⎟ j , 1⎟ |
‖ ‖ ⎝⎝
‖ ⎛SC ⎞ ← ―――――――――――― ⎠ ⎝ ⎠ ⎠ |||
‖ ‖ ‖ ⎜⎝ j , 3⎟⎠ CURSO: DINAMICA ESTRUCTURAL ||| APLICADA A LA
‖ ‖ ‖ k,k
m |
‖ INGENIERIA SISMICA CON |
| | |MATHCAD PRIME
|
‖ ‖
‖ |
‖ SC MÓDULO 2: BASES TEORICAS |

Rigidez Lateral de los pórticos: ⎡1 0 0 0 0 0⎤

aport = ⎢ ⎥
1 ⎢⎣ 0 1 0 0 0 0 ⎥⎦
⎡ 3.673 ⋅ 10 3 -⎛⎝1.755 ⋅ 10 3 ⎞⎠ ⎤
Kp ≔ ⎢ ⎥
1 ⎢⎣ -⎛⎝1.755 ⋅ 10 3 ⎞⎠ 1.356 ⋅ 10 3 ⎥⎦

⎡ 3.629472 ⋅ 10 3 -⎛⎝1.736735 ⋅ 10 3 ⎞⎠ ⎤
Kp ≔ ⎢ ⎥
1 ⎢⎣ -⎛⎝1.736735 ⋅ 10 3 ⎞⎠ 1.343481 ⋅ 10 3 ⎥⎦
⎡ 1 0 0 0 -6 0 ⎤
aport = ⎢ ⎥
2 ⎢⎣ 0 1 0 0 0 -6 ⎥⎦

⎡ 3.753 ⋅ 10 3 -⎛⎝1.749 ⋅ 10 3 ⎞⎠ ⎤
Kp ≔ ⎢ ⎥
2 ⎢⎣ -⎛⎝1.749 ⋅ 10 3 ⎞⎠ 1.38 ⋅ 10 3 ⎥⎦

⎡ 3.710107 ⋅ 10 3 -⎛⎝1.730349 ⋅ 10 3 ⎞⎠ ⎤
Kp ≔ ⎢ ⎥
2 ⎢⎣ -⎛⎝1.730349 ⋅ 10 3 ⎞⎠ 1.366683 ⋅ 10 3 ⎥⎦

⎡ 1 0 0 0 -12 0⎤
aport = ⎢ ⎥
3 ⎢⎣ 0 1 0 0 0 -12 ⎥⎦

⎡ 3.335 ⋅ 10 3 -⎛⎝1.779 ⋅ 10 3 ⎞⎠ ⎤
Kp ≔ ⎢ ⎥
3 ⎢⎣ -⎛⎝1.779 ⋅ 10 3 ⎞⎠ 1.322 ⋅ 10 3 ⎥⎦

⎡ 3.295599 ⋅ 10 3 -⎛⎝1.760345 ⋅ 10 3 ⎞⎠ ⎤
Kp ≔ ⎢ ⎥
3 ⎢⎣ -⎛⎝1.760345 ⋅ 10 3 ⎞⎠ 1.310438 ⋅ 10 3 ⎥⎦

⎡0 0 1 0 0 0⎤ ⎡0 0 1 0 7 0⎤
aport = ⎢ ⎥ aport = ⎢ ⎥
4 ⎢⎣ 0 0 0 1 0 0 ⎥⎦ 5 ⎢⎣ 0 0 0 1 0 7 ⎥⎦

⎡ 5.331 ⋅ 10 3 -⎛⎝2.986 ⋅ 10 3 ⎞⎠ ⎤
Kp ≔ ⎢ ⎥
4 ⎢⎣ -⎛⎝2.986 ⋅ 10 3 ⎞⎠ 2.357 ⋅ 10 3 ⎥⎦

⎡ 5.261566 ⋅ 10 3 -⎛⎝2.946903 ⋅ 10 3 ⎞⎠ ⎤
Kp ≔ ⎢ ⎥ Kp ≔ Kp
4 ⎢⎣ -⎛⎝2.946903 ⋅ 10 3 ⎞⎠ 2.326165 ⋅ 10 3 ⎥⎦ 5 4

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 18 de 19
 CURSO: DINAMICA ESTRUCTURAL APLICADA A LA
INGENIERIA SISMICA CON MATHCAD PRIME
MÓDULO 2: BASES TEORICAS

⎡ -0.6 0 0.8 0 12.8 0 ⎤

aport = ⎢ ⎥
6 ⎢⎣ 0 -0.6 0 0.8 0 12.8 ⎥⎦

Kp ≔ ⎡⎣ 1.347798 ⋅ 10 3 ⎤⎦
6

Matriz de Rigidez del Edificio:

‖ ||
KEDIFICIO ≔ ‖ for i ∊ 1 ‥ Nporticos1
||
‖ ‖ ⎛ ⎛ ⎞ ⎞||
‖ ‖ aglobi ← submatrix ⎜aporti , 1 , cols ⎜Kpi⎟ , 1 , NGLglob⎟ | |
‖ ‖ ⎝ ⎝ ⎠ ⎠||
‖ ‖ Kpg ← a T ⋅ K ⋅ a ||
glob p glob
‖ ‖ i i i i || |
‖ ‖ |
‖ Nporticos1 |
‖ ∑ Kpg |
‖ i=1 i |
‖ |

⎡ 1.112 ⋅ 10 4 -5.227 ⋅ 10 3 -646.943 0 -7.216 ⋅ 10 4 3.151 ⋅ 10 4 ⎤

⎢ ⎥
⎢ -5.227 ⋅ 10 3 4.021 ⋅ 10 3 0 0 3.151 ⋅ 10 4 -2.393 ⋅ 10 4 ⎥
⎢ ⎥
⎢ -646.943 0 1.139 ⋅ 10 4 -5.894 ⋅ 10 3 5.063 ⋅ 10 4 -2.063 ⋅ 10 4 ⎥
KEDIFICIO = ⎢ ⎥
0 0 -5.894 ⋅ 10 3 4.652 ⋅ 10 3 -2.063 ⋅ 10 4 1.628 ⋅ 10 4
⎢ ⎥
⎢ -7.216 ⋅ 10 4 3.151 ⋅ 10 4 5.063 ⋅ 10 4 -2.063 ⋅ 10 4 1.087 ⋅ 10 6 -4.602 ⋅ 10 5 ⎥
⎢ ⎥
⎢⎣ 3.151 ⋅ 10 4 -2.393 ⋅ 10 4 -2.063 ⋅ 10 4 1.628 ⋅ 10 4 -4.602 ⋅ 10 5 3.519 ⋅ 10 5 ⎥⎦

Realizado por: Alexis Pompilla Yábar, no se permite la reproducción total o parcial sin consentimiento del autor. Página 19 de 19