File Exchange

image thumbnail

OpenSeismoMatlab

version 3.2 (1.13 MB) by George Papazafeiropoulos
Open source software for strong ground motion data processing

48 Downloads

Updated 13 Aug 2019

View Version History

View License

OpenSeismoMatlab can calculate the following strong motion data processing output: (1) Time histories of acceleration, velocity and displacement, with or without baseline correction, (2) Resampled time histories of acceleration, velocity and displacement with modified time step size, (3) Peak ground acceleration, velocity and displacement (PGA, PGV, PGD), (4) Total cumulative energy and its normalized time history, (5) Significant duration according to Trifunac & Brady (1975), (6) Arias intensity, (7) Linear elastic response spectra, (8) Constant ductility inelastic response spectra, (9) Fourier amplitude spectrum, (10) Mean period (Tm)

Cite As

George Papazafeiropoulos (2021). OpenSeismoMatlab (https://www.mathworks.com/matlabcentral/fileexchange/67069-openseismomatlab), MATLAB Central File Exchange. Retrieved .

Comments and Ratings (24)

@Ziad Bennour:
Thanks for the rating. Please ensure that the file "elcentro.dat" is inside your working directory.

Ziad Bennour

Hey, I'm getting this error with Matlab 2021a

Error using textscan
Invalid file identifier. Use fopen to generate a valid
file identifier.

Error in Verification (line 9)
text=textscan(fid,'%f %f');

Any reason why?

毅斌 张

Thanks for your consideration. We are open to such proposals.

Glenn Thompson

This could be a good set of tools to incorporate within GISMO, which is object-oriented and has functions for reading and writing different seismic data formats.

Caglar Kizilpinar

Caglar Kizilpinar

long li

@Ming Yin Samson LAU:
Thank you for your interest in this submission. Let me know if you have any other issues!!!

Ming Yin Samson LAU

Hi sorry seems i figured out what's the problem
i did not download the el dat file and put it under the same path so the file could not be opened
thanks you for your work, you help a lot for my paper

Ming Yin Samson LAU

Hi i am a new user to matlab since i was doing my jornal paper on ground motion selection and scling
After downloading your add-on i tried to use it inside matlab,
while i copy your first few script on load earthquake data, it shows
Error using textscan
Invalid file identifier. Use fopen to generate a valid file identifier.
Am i doing anything wrong? Please help me

younes farez

Fabrice Hollender

@ michele paoletti: The code has been updated so that this issue is resolved. Thank you for your suggestions and apologies for the delayed response.

Cong-Thuat DANG

Hamdullah Livaoglu

Ruiyang Zhang

Ruiyang Zhang

michele paoletti

Thanks for answer George. td is considered as time duration of the seismic wave. For example if we use definition of Trifunac duration td is the time where the energy go from 5% to 95% of the total energy. Reducing time duration respect your formula we take in account only energy of the wave and not the noise energy (although very low) that is present first and then the seismic wave. The result obtained is a slightly lower value respect your formula. I'd like to know what you think about.

OpenSeismoMatlab calculates the Arias intensity in the whole time duration of the input earthquake acceleration time history. This is the reason for which there are not any limits in the code segment. What do you mean by td?

michele paoletti

Hi, Arias Intensity must be calculated in the interval from 0 to td, while this line code that you post not have this limit:
aint2 = cumsum(A.^2)*pi*dt/(2*g);
Why?

Thank you very much Ricardo A. The code has been fixed and enhanced based on your comments. You can now download the latest version. Regards!

Ricardo A.

2 errors:
line 188 should be: seismic.vel = cumtrapz(time,xgtt);
line 190 should be: seismic.disp = cumtrapz(time,seismic.vel);
previously, the velocity was multiplied by an extra 'dt' and the displacement by an extra 'dt^2'. Thanks for your work!

Ricardo A.

Also, maybe you can add these 2 parameters calculations:
% MEAN FREQUENCY
fi = f(f>0.25 & f<20);
Ci = U(f>0.25 & f<20);
Fm = ((Ci(:)'.^2)*(fi(:)))/(Ci(:)'*Ci(:));
seismic.Fm = Fm;

% PREDOMINANT PERIOD & FREQUENCY
maxi=0;
indice=0;
for i=1:length(seismic.PSa);
if seismic.PSa(i)>maxi
maxi=seismic.PSa(i);
indice=seismic.Period(i);
else
maxi=maxi;
indice=indice;
endif
end
seismic.PredSpec=maxi;
seismic.PredPeriod=indice;
seismic.PredFreq=1/maxi;
maxi=0;
indice=0;

MATLAB Release Compatibility
Created with R2017b
Compatible with any release
Platform Compatibility
Windows macOS Linux
Acknowledgements

Inspired: ArtifSeismoMatlab

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

OpenSeismoMatlab