EL-Centro Earthquake 1940 N-S response analysis
by Nasser Abbasi

In[1]:=

Remove["Global`*"] ;

SetDirectory["E:/nabbasi/data/nabbasi_web_Page/my_courses/UCI_COURSES/CREDIT_COURSES/fall_2006/CEE_247/HWs/HW4"] ;

Remove :: rmnsm : There are no symbols matching \"Global`*\". More…

Load El-Centro data and plot the earthquack recorded accelration

In[3]:=

data = Import["el_centro.txt", "Table"] ;

dataLength = 8 Length[data] ;

g = Table[{0, 0}, {dataLength}] ;

t = Table[0, {dataLength}] ;

u = Table[{0, 0}, {dataLength}] ;

v = Table[{0, 0}, {dataLength}] ;

acc = Table[{0, 0}, {dataLength}] ;

delT = 0.02 ; (*sec*)k = 0 ;

currentTime = 0 ;

ListPlot[g, PlotJoined→True, PlotRange→All, AxesLabel→ {"Time sec", "Acc. g"}, PlotLabel->"El centro earthquack 1940 N-S"] ;

[Graphics:HTMLFiles/index_15.gif]

Implementation of Linear Accelration step-by-step algorithm

In[38]:=

Plot results

In[46]:=

ListPlot[u, PlotJoined→True, PlotRange→All, AxesLabel→ {"Time sec", "u(t) inch"}, PlotLabel->"Displacment from El centro earthquack 1940 N-S"] ;

[Graphics:HTMLFiles/index_18.gif]

In[47]:=

ListPlot[v, PlotJoined→True, PlotRange→All, AxesLabel→ {"Time sec", "v(t) in/sec"}, PlotLabel->"speed from El centro earthquack 1940 N-S"] ;

[Graphics:HTMLFiles/index_20.gif]

In[48]:=

[Graphics:HTMLFiles/index_22.gif]

Created by Mathematica  (November 13, 2006) Valid XHTML 1.1!