Module

for

Earthquake Induced Vibrations

 

Background.

    In the study of earthquake induced vibrations on multistory buildings, the free transverse oscillations satisfy the equation  

(1)        [Graphics:Images/EarthQuakeModelMod_gr_1.gif],  

where the forces acting on the i-th floor are  
    
        [Graphics:Images/EarthQuakeModelMod_gr_2.gif].  

Consider a building with  n floors each of mass  m  slugs and the horizontal restoring force of   k tons/foot  between floors. Then this system reduces to the form  

(2)        [Graphics:Images/EarthQuakeModelMod_gr_3.gif],
where
        [Graphics:Images/EarthQuakeModelMod_gr_4.gif] .  

Proof  Earthquake Models  Earthquake Models  

Computer Programs  Earthquake Models  Earthquake Models  

 

Example 1.  Consider a building with  n = 6  floors each of mass  m = 1250 slugs (weight of 20 tons)
and the horizontal restoring force of   k = 10,000 lb/ft = 5 tons/foot  between floors.
Then  [Graphics:Images/EarthQuakeModelMod_gr_5.gif],  and this system reduces to the form  

Compute the eigenvalues of matrix  [Graphics:Images/EarthQuakeModelMod_gr_6.gif],
and the natural frequencies  [Graphics:Images/EarthQuakeModelMod_gr_7.gif]  and  periods  P  of oscillation of the building.
Solution 1.

 

More Background.

    A horizontal earthquake oscillation of amplitude  [Graphics:Images/EarthQuakeModelMod_gr_39.gif]  of the form  [Graphics:Images/EarthQuakeModelMod_gr_40.gif]  will produce an acceleration  [Graphics:Images/EarthQuakeModelMod_gr_41.gif], and the opposite internal force on each floor of the building is   [Graphics:Images/EarthQuakeModelMod_gr_42.gif].  The resulting non-homogeneous system is   

    [Graphics:Images/EarthQuakeModelMod_gr_43.gif],   where   [Graphics:Images/EarthQuakeModelMod_gr_44.gif].  

 

Example 2.  Solving the above non-homogeneous system for the coefficient vector  v  for  X[t].  
The vector  v  is the solution to the equation  [Graphics:Images/EarthQuakeModelMod_gr_45.gif].  
Use the earthquake amplitude  e = 0.075 ft = 0.9 in.  for this example.

Solve the linear system using the parameters  [Graphics:Images/EarthQuakeModelMod_gr_46.gif] and  e = 0.075.  

Find the coefficient vector  v  and the vector  X[t].  Plot the vibrations of each floor.
Solution 2.

 

Example 3.  For the above non-homogeneous system the coefficient vector  v  is the solution to the equation  

        [Graphics:Images/EarthQuakeModelMod_gr_73.gif].  

Plot the maximum amplitude of oscillation of the floors vs the parameter  [Graphics:Images/EarthQuakeModelMod_gr_74.gif]  over the interval  [Graphics:Images/EarthQuakeModelMod_gr_75.gif],  
this graph should have six vertical asymptotes corresponding to each value  [Graphics:Images/EarthQuakeModelMod_gr_76.gif]  in the table above.  
Then plot the maximum amplitude as a function of the period  p  in seconds.
Solution 3.

 

Old Lab Project (Earthquake Induced Vibrations  Earthquake Induced Vibrations).  Internet hyperlinks to an old lab project.  

 

Research Experience for Undergraduates

Earthquake Models  Earthquake Models  Internet hyperlinks to web sites and a bibliography of articles.  

 

Download this Mathematica Notebook Earthquake Induced Vibrations

 

Return to Numerical Methods - Numerical Analysis

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004