Taller #1 Dinamica Estructural

TALLER # 1
JONATAN MIGUEL VILORIA ROJAS
DINMICA ESTRUCTURAL Y DISEO SSMICO

FABIAN CONSUEGRA Ph. D.
DIRECCIN DE ESPECIALIZACIONES
ESPECIALIZACIN EN ANLISIS Y DISEO DE ESTRUCTURAS
FUNDACIN UNIVERSIDAD DEL NORTE

BARRANQUILLA-COLOMBIA
25 SEPTIEMBRE DE 2015
TABLA DE CONTENIDO
1. SOLUCIN DEL PROBLEMA # 1

1.1Solucin teniendo en cuenta solo la masa de 100 Kg.
impuesta.
1.2Solucin teniendo en cuenta solo la masa de la viga.
1.3Solucin teniendo en cuenta todas las masas.
sobre
2. SOLUCION DEL PROBLEMA # 2

2.1 Modelacin en SAP2000.
3. SOLUCION AL PROBLEMA # 3
3.1 Graficas de la respuesta del sistema con MATLAB.
4.1Evaluacin de los efectos de las condiciones iniciales y el
amortiguamiento en la respuesta.
4.2Graficas combinaciones de u0 y v0 con amortiguacin de 1%
4.3Graficas de combinaciones de u0 y v0 con amortiguamiento de 5%.
5.1Graficas de respuesta del sistema a cada frecuencia de excitacin (F)
cada 0.2 Hz hasta 4 Hz.
5.2Clculo y grafica de la respuesta en frecuencia del sistema.
5.3Calculo del amortiguamiento del sistema con el mtodo de halfpower.
1. SOLUCIN DEL PROBLEMA # 1

DATOS:
fc= 21MPa, L=2.5 m, M=100 Kg
1.1 Solucin de la frecuencia natural teniendo en cuenta solo la masa
de 100 Kg. Sobre impuesta.
( E) .
Hallamos el mdulo de elasticidad
E=3900 f ' c=3900 21 MPa E=17872 MPa
Hallamos la inercia
( I)
de la seccin.
bh3 0.5 m0.5 m3

I=
=
=I =0.005208 m2
12
12
La rigidez
(K ) de una viga empotrada en un extremo y simplemente
apoyada en el otro extremo est dada por:
3
7
768EI
768178720.005208
K=
=K=
3
7L
Hallamos la frecuencia natural en
(rad / s).
K
81700000 N . m
= =
==903.88 rad /s
M
100 Kg.
=
En ciclos por segundo.
f=
903.88 rad /s
=f =
=f =143.85 Hz
2
2
Periodo de natural en segundos.
T=
2
2
=T =
=T =0.00695 s
903.88rad /s
1.2 Solucin de la frecuencia natural teniendo en cuenta solo la masa

de la viga.
=2400
Tenemos que la densidad del concreto por
m ( viga )=2400
Kg
m3
Kg
0.5 m0.5m5 m=m ( viga )=3000 Kg
m3
( E) .
E=17872 MPa
Hallamos la inercia
( I)
de la seccin.
I =0.005208 m2
La rigidez
3
7
768EI
768178720.005208
K=
=K=
3
7L
(rad / s).
K
81700000 N . m
= =
==165.02rad / s
M
3000 Kg.
=
f=
165.02rad / s
=f =
=f =26.26 Hz
2
2
Periodo natural en segundos.
T=
2
2
=T =
=T =0.038 s
165.02 rad /s
1.3 Solucin de la frecuencia natural

masas.
La masa total del sistema seria
teniendo en cuenta todas las
100 Kg. de la carga+ 3000 Kg . viga=3100 Kg .
( E) .
E=17872 MPa
Hallamos la inercia
( I)
de la seccin.
I =0.005208 m2
La rigidez
3
7
768EI
768178720.005208
K=
=K=
3
7L
(rad / s).
K
81700000 N . m
= =
==162.34 rad /s
M
3100 Kg.
=
f=
162.34 rad /s
=f =
=f =25.83 Hz
2
2
T=
Periodo en segundos.
2
2
=T =
=T=0.0038 s
162.34 rad /s
2. SOLUCIN DEL PROBLEMA #2

DATOS:
L=5 m
E=200000 MPa
Seccin cuadrada l=50 cm
Hallamos la inercia de la seccin.
0.5 m
3
0.5 m
bh3
I=
=I =
12
La rigidez
(K ) de una viga empotrada en un extremo y en voladizo en el otro
extremo est dada por:
3 EI
32000000.00508
K= 3 =K =
L
Hallamos la masa de la barra sabiendo que densidad del acero es
=7850 Kg/m3
2
m=( 0.5 m ) 5 m7850
m=
Kg
=mbarra=9812.5 Kg.
m3
9812.5 Kg.
=m =
=m =3270.83 Kg .
( masa
)
3
3
Hallamos
u0
0= 0.00128 m
m.g
32.08 kN
u0=
=u 0=
=u
K
25000 kN .m
f=
K
25000000 N .m
rad
=
==87.42
m .
3270.83 Kg .
s
=
2
T=
rad /s
rad
s
=f =13.913 Hz .
2
87.42
2
2
=T =
=T =0.07187 s
87.42rad / s
2.1
Modelacin en SAP2000.
Teniendo los valores calculados se precede a modelar el sistema en SAP2000
Se calcula el periodo T en el programa
Se calcula el desplazamiento
T =0.07196 s
u0=0.00198 m
Se puede observar que los valores son muy parecidos a los calculados
manualmente.
3. SOLUCIN AL PROBLEMA # 3
Se retoman los valores calculados en el problema # 1 especficamente donde
se analizan todas las masas tanto de la viga como de la carga de 100 Kg.
La masa total del sistema seria
100 Kg. de la carga+ 3000 Kg . viga=3100 Kg .
( E) .
E=17872 MPa
Hallamos la inercia
( I)
de la seccin.
I =0.005208 m
Hallamos la rigidez
(K ) .
3
7
768EI
768178720.005208
K=
=K=
3
7L
(rad / s).
K
81700000 N . m
= =
==162.34 rad /s
M
3100 Kg.
=
f=
162.34 rad /s
=f =
=f =25.83 Hz
2
2
T=
Periodo en segundos.
2
2
=T =
=T=0.0038 s
162.34 rad /s
Hallamos el desplazamiento inicial
u
( 0) teniendo en cuenta la carga de
1000 N justo encima de la masa M.

La carga total seria:
P=( MasaVigag ) + ( masa 100 Kg.g ) +(1000 N ) entonces
P=29.4 kN +0.981 kN +1 kN =P=31.38 kN

u0=
P
31.38 kN
=
=u =0.0003841 m
K 81700 kN . m 0
Calculamos la respuesta del sistema partiendo del reposo quitando

sbitamente la carga M.
Como parte delresposo u 0=0 entonces la respuesta del sistema quedara:

u=u 0 cos ( t ) , se reemplazan los valores de =162.34
u=(0.0003841 m)cos (162.34
rad
y u0 =0.0003841 m
s
rad
t )
s
Se grafica la respuesta del sistema en MATLAB en un intervalo de tiempo de 2

s.
% Frmula vibraciones libres sin amortiguamiento %

clear all, clc, close all
u0
v0
f
w
dt
t
u
v
=
=
=
=
=
=
=
=
0.0003841; % Desplazamiento inicial [m]

0;
% Velocidad inicial [m/s]
25.83;
% frecuencia natural [Hz]
2*pi*f;
% frecuencia angular [rad/s]
0.01;
0:dt:2;
% Tiempo s
u0*cos(w*t) + v0/w*sin(w*t);
-w*u0*sin(w*t) + v0*cos(w*t);
figure,subplot(211),plot(t,u),title('Vibraciones libres sin

amortiguamiento')
xlabel('t [seg]'),ylabel('Despl. [m]')
subplot(212),plot(t,v),title('Vibraciones libres sin
amortiguamiento')
xlabel('t [seg]'),ylabel('Velocidad [m/s]')
vcal = diff(u)/dt;
figure,plot(t,v,t(1:end-1),vcal),
title('Comparacin velocidades')
xlabel('t [seg]'),ylabel('Velocidad [m/s]'),

legend('analit.','calc.')
3.1 Graficas de las respuestas del sistema.
4. SOLUCIN AL PROBLEMA #4
A continuacion se muestran todas las graficas de las respuestas con todas las
u0
combinaciones posibles de los valores para
v 0,
con amortiguacion de
1% y 5% respectivamente, se tom un tiempo de 30 s para generar los

graficos.
4.1 Evaluacin de los efectos de las condiciones iniciales y el

amortiguamiento en la respuesta.
Despus de analizar las grficas se pudo observar primeramente que
los porcentajes de amortiguamiento marcan una gran diferencia en
las respuestas de los sistemas, se puede notar fcilmente que el
sistema con menos amortiguacin en el tiempo que se estudi (30 s)
no llego al reposo, mientras que la respuestas con un 5% de
amortiguacin con un tiempo de (20 s) ya se encontraban
aproximadamente en estado de reposo.
Las condiciones iniciales de
u0
v0
en las respuestas muestran
que, cuando el desplazamiento inicial es menor que la velocidad la

curva del decremento logaritmico de la amortiguacion se separa mas
de los picos del desplazamiento, esta separacin va disminuyendo
acorde como los valores del desplazamiento se igualan a los de la
velocidad.
En las respuestas donde el desplazamiento no es tan diferente de la
velocidad se puede observar que la envolvente de la amortiguacion
coincide total mente con el drecremento del desplazamiento.
Hallamos los valores correspondientes a las constantes A y B.
t=0 y u 0 es igual a :
u0=e ( 0) ( Asin d ( 0 ) +Bcos d ( 0 ) )+
u0=1( B)+
B=u0 +
P0
1
( ( 1 2 ) sin ( 0 ) ( 2 ) cos ( 0 ) )
2
2
2
k ( 1 ) + ( 2 )
P0
1
(( 2 ))
2
k ( 1 ) 2+ (2 )2
P0
2
2
k ( 1 )2 + ( 2 )2
Se deriva la funcin para hallar el valor de la constante A

t
u=e
( Asin d t +Bcos d t )+
P0
1
2
( ( 1 ) sin t ( 2 ) cos t )
2
2
2
k ( 1 ) + ( 2 )
u =et [ ( d Acos d td B sin d t ) ] + ( e t ) [ Asin d tBcos d t ] +
u0=d AB+
A=
P0
1
[
( 1 2) ]
2
2
2
k ( 1 ) + ( 2 )
P0 [ ( 1 2 ) ]
1
[u0 +
]
d
k ( 1 2 )2 + ( 2 )2
P0
1
[ ( 1 2 ) cos
2
2
2
k ( 1 ) + ( 2 )
Con los valores de las constantes obtenidos se modelan en Matlab, teniendo en

cuenta todos los valores de entrada que se piden en el enunciado.
5.1 Graficas de respuesta del sistema a cada frecuencia de

excitacin (F) cada 0.2 Hz hasta 4 Hz.
Frecuencia de 0.2 Hz
Frecuencia
de 4 Hz
5.2
Clculo y grafica de la respuesta en frecuencia del sistema.
Se tomaron los desplazamientos mximos en cada grafica dependiente de

cada una de las frecuencias de excitacin.
5.3 Calculo del amortiguamiento del sistema con el mtodo de
half-power.
Segn el mtodo y la grfica se tiene lo siguiente:
Dmax= 0.4156
Dmax
=0.293 8
2
Segn la grafica
1=11.23, d=12.86, 2=14.5
2 1
1
=2 =
n
Q
14.511.23
=0.25 4
12.86
2 =0.254 despejamos y tenemos que =0.127

Taller #1 Dinamica Estructural

Cargado por

Copyright:

Formatos disponibles

Taller #1 Dinamica Estructural

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Taller #1 Dinamica Estructural

Cargado por

Copyright:

Formatos disponibles

TALLER # 1

JONATAN MIGUEL VILORIA ROJAS

DINMICA ESTRUCTURAL Y DISEO SSMICO

FUNDACIN UNIVERSIDAD DEL NORTE

1. SOLUCIN DEL PROBLEMA # 1

2. SOLUCION DEL PROBLEMA # 2

1. SOLUCIN DEL PROBLEMA # 1

Hallamos el mdulo de elasticidad

E=3900 f ' c=3900 21 MPa E=17872 MPa

bh3 0.5 m0.5 m3

(K ) de una viga empotrada en un extremo y simplemente

apoyada en el otro extremo est dada por:

Hallamos la frecuencia natural en

En ciclos por segundo.

Periodo de natural en segundos.

1.2 Solucin de la frecuencia natural teniendo en cuenta solo la masa

Tenemos que la densidad del concreto por

Hallamos el mdulo de elasticidad

(K ) de una viga empotrada en un extremo y simplemente

apoyada en el otro extremo est dada por:

Hallamos la frecuencia natural en

En ciclos por segundo.

Periodo natural en segundos.

1.3 Solucin de la frecuencia natural

teniendo en cuenta todas las

100 Kg. de la carga+ 3000 Kg . viga=3100 Kg .

Hallamos el mdulo de elasticidad

(K ) de una viga empotrada en un extremo y simplemente

apoyada en el otro extremo est dada por:

Hallamos la frecuencia natural en

En ciclos por segundo.

2. SOLUCIN DEL PROBLEMA #2

(K ) de una viga empotrada en un extremo y en voladizo en el otro

extremo est dada por:

m=( 0.5 m ) 5 m7850

Teniendo los valores calculados se precede a modelar el sistema en SAP2000

Se calcula el periodo T en el programa

100 Kg. de la carga+ 3000 Kg . viga=3100 Kg .

Hallamos el mdulo de elasticidad

Hallamos la frecuencia natural en

En ciclos por segundo.

Hallamos el desplazamiento inicial

1000 N justo encima de la masa M.

P=( MasaVigag ) + ( masa 100 Kg.g ) +(1000 N ) entonces

P=29.4 kN +0.981 kN +1 kN =P=31.38 kN

Calculamos la respuesta del sistema partiendo del reposo quitando

Como parte delresposo u 0=0 entonces la respuesta del sistema quedara:

u=(0.0003841 m)cos (162.34

Se grafica la respuesta del sistema en MATLAB en un intervalo de tiempo de 2

% Frmula vibraciones libres sin amortiguamiento %

0.0003841; % Desplazamiento inicial [m]

figure,subplot(211),plot(t,u),title('Vibraciones libres sin

xlabel('t [seg]'),ylabel('Velocidad [m/s]'),

combinaciones posibles de los valores para

1% y 5% respectivamente, se tom un tiempo de 30 s para generar los

4.1 Evaluacin de los efectos de las condiciones iniciales y el

en las respuestas muestran

que, cuando el desplazamiento inicial es menor que la velocidad la

Hallamos los valores correspondientes a las constantes A y B.