UNIVERSIDAD TÉCNICA DE MANABÍ
FACULTAD DE CIENCIAS MATEMATICAS FISICAS Y
QUIMICAS
DEPARTAMENTO DE CONSTRUCCIONES CIVILES
TAREA#1
ESTUDIANTE:
KAVIANA GENITH ALAVA MACIAS
ASIGNATURA:
INGENIERIA SISMICA
PARALELO:
“D”
SEMESTRE:
NOVENO
DOCENTE:
ING. MARCOS ZEVALLOS
PERIODO:
OCTUBRE 2019 – FEBRERO 2020
EJERCICIO EN CLASE
function [q]=v_libre(zi,w,qo,qpo)
%
% Vibraciones libres en sistemas de un grado de libertad
%
% Por: Roberto Aguiar Falconi
% CEINCI-ESPE
%-------------------------------------------------------------
% v_libre(zi,w,qo,qpo)
%-------------------------------------------------------------
% zi: Vector que contiene dos factores de amortiguamiento para
% los cuales se encuentra la respuesta en el tiempo.
% w : frecuencia natural del sistema de 1 gdl.
% qo: desplazamiento en t=0
% qpo: velocidad en t=0
% tmax: tiempo maximo de la respuesta igual a 0.6 segundos.
t=linspace(0,0.6,1000)';
np=length (zi); %np es el numero de factores de amortiguamiento
for i=1:np
if zi(i)<1
wa=w*sqrt(1-zi(i)*zi(i));
B=qo; A=(qpo+zi(i)*w*B)/wa;
q1=(A*sin(wa*t)+B*cos(wa*t));
q2=exp(-zi(i)*w*t);
q(:,i)=q2.*q1;
elseif zi(i)==1
B=qo; A=qpo+B*w;
q(:,i)=(A*t+B).*exp(-w*t);
else
landa1=-zi(i)*w+w*sqrt(zi(i)*zi(i)-1);
landa2=-zi(i)*w-w*sqrt(zi(i)*zi(i)-1);
C=[1 1; landa1 landa2]; D=[qo; qpo];
X=C\D; A=X(1); B=X(2);
q(:,i)=A*exp(landa1*t)+B*exp(landa2*t);
end
endplot (t,q(:,1), 'r', 'LineWidth',2); hold on; plot (t,q(:,2),
'LineWidth',2)
xlabel ('Tiempo'); ylabel ('Desplazamiento')
title ('Vibracion libre en sistema de 1 gdl')
%---fin---
SOLUCIÓN
>> zi=[0.05;0.25]; w=31.416; qo=2;qpo=10;
>> v_libre(zi,w,qo,qpo)
GRAFICA
o Cuando: zi=[0.00;0.00]
SOLUCIÓN
>> zi=[0.00;0.00]; w=31.416; qo=2;qpo=10;
>> v_libre(zi,w,qo,qpo)
GRAFICA
o Cuando: zi=[0.05;0.05]
SOLUCIÓN
>> zi=[0.05;0.05]; w=31.416; qo=2;qpo=10;
>> v_libre(zi,w,qo,qpo)
GRAFICA
o Cuando: zi=[1.2;1.2]
SOLUCIÓN
>> zi=[1.2;1.2]; w=31.416; qo=2;qpo=10;
>> v_libre(zi,w,qo,qpo)
GRAFICA
DEDUCCIÓN DE LA FORMULA
mg=kδ(1.1)
k ( q +δ )+ c q̇ +m q̈−mg=0
Al sustituir (1.1) en ésta última ecuación, se tiene:
m q̈+c q̇+ kq=0(1.2)
Se conoce que la frecuencia natural W n y el periodo de vibración T, valen:
W n=
√ k
m
T=
2π
Wn
(1.3)
Por otra parte, se define el factor de amortiguamiento ξ como:
c
ξ= (1.4)
2 √ mk
Si la ecuación diferencial (1.2) se divide para m se tiene:
c 2
q+ q+W n q=0
m
Al multiplicar y dividir el término c/m por 2 √ mk y al utilizar la ecuación (1.4) se tiene:
c c 2 √ mk
= =2 ξ W n
m 2 √ mk m
Luego otra forma de presentar la ecuación diferencial del movimiento es:
2
q+ 2 ξ W n q+W n q=0 (1.5)
Solución de la ecuación diferencial
Se plantea la solución de la ecuación diferencial (1.5) de la siguiente forma:
q ( t )=a e λt (1.6)
Donde a es una constante de integración y λ es una variable a determinar. Al derivar la ecuación
(1.6) con respecto al tiempo y reemplazar en (1.5) se tiene:
λt
q̇=a λ e
q̈˙ =a λ e
2 λt
Al reemplazar en la ecuación diferencial, se tiene:
a λ 2 e λt +2 ξ W n a e λt +W 2n a e λ t=0
λt 2 2
a e ( λ +2 ξ W n λ+W n)=0
Para que la última ecuación sea igual a cero es necesario que la cantidad del paréntesis sea cero.
2 2
λ +2 ξ W n λ+W n =0
−2 ξ W n ± √ 4 ξ W n−4 W n
2 2
λ=
2
λ=−ξ W n ±W n √ ξ −1(1.7)
2
Las raíces de λ dependen del valor de ξ ya que el radical puede ser positivo, cero o negativo.
Vibración libre sin amortiguamiento
En este caso ξ = 0, es un caso ideal que significa que la estructura queda vibrando
indefinidamente. Al ser ξ = 0 las raíces que se obtiene de (1.7) son:
λ=± W n √−1
Luego la solución se transforma en:
q ( t )= A cos (W ¿¿ n t )+ B sin(W ¿ ¿ nt )=C sin ( W n t +γ ) ¿ ¿
C=√ A +B (1.8)
2 2
Siendo γ el ángulo de fase y c la amplitud máxima.
