Matlab : Plot Pseudo Response Spectrum for El Centro Data

Matlab : Plot Pseudo Response Spectrum for El Centro Data
#matlab #response #elcentro
The purpose of this video is to show you how to plot the pseudo-acceleration and pseudo-velocity response spectrum of EL CENTRO earthquake data by using Matlab code. The method to solve the equation of motion is based on interpolation of excitation which is popular for solving this type of equation. At the end of the video, you can learn how to write a function to solve the equation of motion and implement it in the main Matlab program to extract the data for earthquake record and use it to plot the graphs.
This tutorial is explained using imperial unit.
Matlab code file: bit.ly/2YMq2zi
Plot Pseudo Response Spectrum for El Centro Data using Python:
• Python : Plot Pseudo R...
Optimization playlist: • Optimization
- PSO algorithm: • Python Code for Partic...
- Size optimization using PSO: • Size Optimization of T...
- Size and shape optimization using PSO: • PSO and Python for siz...
Truss structural analysis using Python:
- 2D truss: • Python: Direct Stiffne...
- 3D truss: • Python: Direct Stiffne...
Direct stiffness method: • Direct Stiffness Metho...
Truss structural analysis using Matlab:
- Part 1: • Matlab : Direct Stiffn...
- Part 2: • Matlab : Direct Stiffn...
-------------------------------------------------------------------
Donate me: / engineerhunter
Join membership: / @engineerhunter
TH-cam: / engineerhunter
Facebook: / engineerhunterse
Instagram: / engineerhunter
TikTok: www.tiktok.com...
Contact me: [e-mail]

[@EngineerHunter]

You can download the Matlab code file on my Patreon page. bit.ly/2YMq2zi

[@Abdullahtaha-lh4hg]

great video , however, what should I do if I want to plot the 3 figures in 4-Way logarithmic figure ??

[@user-hd6ii6qs2v]

What is your El Centro acceleration unit? gal or g? Thanks for your video teaching!

[@user-hd6ii6qs2v]

oh i know it's g

[@EngineerHunter]

Yes, it’s g. Thanks for watching.

[@IEmmaReader]

If i wanted my code to return the entire response for one loop instead of the max value, what would i change?

[@EngineerHunter]

You can remove max() function from the loop to return all “u” or “v” value for one period. And if you want to calculate only at one period T, you don’t need to use main loop for i=1:length(T).

[@faisalaziz6331]

thanks for the video but I have a question if we wanna find system acceleration how do we do it and how do we code it in matlab thanks

[@georgiclarke2769]

Hi there, when I do this and plot, my units for displacement, displacement and acceleration are 10^-4 out? Any idea why this could be? When importing the raw acceleration data I have the acceleration in g units and time in seconds. Any help would be greatly appreciated

[@EngineerHunter]

Be sure to check all the units for each parameters. I have no idea why without seeing your code.

[@Nikesh977]

Hi, thanks for the video, Can you modify the code to show response for multiple damping ratios. and show them in same graph? eg. use xi = [xi1, xi2, xi3, xi4]. Thanks.

[@anubhab_IITG]

Thankyou for the video. It was very helpful. I am an absolute beginner in terms of learning a matlab, can u tell me what is the meaning when u used i = length (e) -1 ? What does it mean? And why u used for loop two times? And how can I strengthen up writing a code and to get acquainted with all related stuffs. I am a student of earthquake engineering.

[@EngineerHunter]

i = 1:length(E)-1 means that loop will run through the acceleration list/vector. (This is for solving differential equation).
Loop two times become each system has different Natural Period (Tn from 0.01 to 3 sec). (This is for finding pseudo response spectra).

[@MuhammetBilalAyan]

How can I do that in excel? can you share me any resource about that?

[@furkanylmaz8511]

sir, could you share the article seemed in tutorial ?

[@EngineerHunter]

Thanks for your question. Both Numerical method and Definition of Response Spectra that I mentioned in is this video were based on book “Dyanamics of Structures” by Anil K. Chopra. You can check it out.

[@carlosalexandreseruti1921]

Is anyone has this code in Mathcad, please?

[@abaich90]

Sir, did you find this code in MathCAD?

[@ali.al-askary]

i have use the same code it says errors.

[@EngineerHunter]

Could you point out which line causes the program to error?

[@jesusvalenzuela9484]

For me it says error on the u = zeroes(length(E),1)

[@georgen9755]

pirated versions with same code it says errors