Elastic Response Spectra Script

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

silvia
Posts: 3912
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia » Sat Apr 26, 2008 8:00 pm

that's a dynamics question.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104

ebe79442114
Posts: 73
Joined: Tue Feb 20, 2007 6:48 am
Location: IIEES
Contact:

Post by ebe79442114 » Sun May 18, 2008 11:54 am

I wrote a program with Matlab language, that correct the row acceleration data and draw the response spectra.I released it under Gnu general public licence.
It draws response spectra with three methods.
address:http://rapidshare.com/files/62694518/Se ... ar.gz.html

mkumar07
Posts: 14
Joined: Mon Mar 10, 2008 5:37 am
Location: LONDON

Results not matching with SeismoSignal

Post by mkumar07 » Mon Jun 22, 2009 8:40 am

Hi there!

I have performed the elastic time history of 124 ground motions using this code; the results are seriously contrasting to Seismosignal (another software). I am very sure results of Seismosignal are right.

But, I cannot find any mistake whatsoever in this code (except for the damping) but cannot see anyother mistake.

Please help.

Kumar, M.
Mukesh Kumar
PhD student,
Imperial College of London, UK

silvia
Posts: 3912
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia » Mon Jun 22, 2009 10:51 pm

look for a systematic difference/error.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104

mkumar07
Posts: 14
Joined: Mon Mar 10, 2008 5:37 am
Location: LONDON

No consistent error

Post by mkumar07 » Tue Jun 23, 2009 1:37 am

Dear Silvia,

I couldn't find a systematic error at all; the difference in both software is random. It is kind of strange that such a simple code which seems logically correct is causing this difference in results.

Thanx for your prompt response.

Kumar, M.
Mukesh Kumar
PhD student,
Imperial College of London, UK

foadEng
Posts: 28
Joined: Sun Jan 24, 2010 12:51 pm
Location: Iran

Post by foadEng » Mon May 17, 2010 12:14 am

Hi
The topic here is about elastic response spectra. I wanted to know about inelastic response spectra. Can we generate inelastic response spectra with this algorithm?
Than you very much!
Foad Mohajeri Nav
Iran University of Science and Technology

felipe_cordero
Posts: 18
Joined: Wed Jun 02, 2010 2:05 pm
Location: Universidad de Chile

observation

Post by felipe_cordero » Tue Oct 26, 2010 12:13 pm

if you have a too many ground motion files, the script fails because it opens every time the tmpEnvelop files.

It happens when we use;

set Dcontent [read [open tmpEnvNodalDisp.out ]];
set Di [lindex $Dcontent 2];
puts Sd$Di
set iDi "$iDi $Di";

I recommend to use the following lines:

set in [open tmpEnvNodalDisp.out];
set Dcontent [read $in];
set Di [lindex $Dcontent 2];
#puts Sd$Di
set iDi "$iDi $Di";
close $in;

This way closes the tmp file once it is used.

Thanks any way for the script, It has been very useful for me.

naba077
Posts: 34
Joined: Mon May 28, 2012 9:38 pm
Location: Saitama University

Re: Elastic Response Spectra Script

Post by naba077 » Wed Jan 16, 2013 1:29 am

I'm sorry, but when I run the program, it says that
"can't read dt: no such variable"
What does that mean?
Where should I define dt and what is dt?

naba077
Posts: 34
Joined: Mon May 28, 2012 9:38 pm
Location: Saitama University

Re: Elastic Response Spectra Script

Post by naba077 » Wed Jan 16, 2013 1:44 am

And also how to define these variables?
"
set iGroundFileSMSR {"IV79cal" "IV79chi" "IV79e01" "IV79e12" "IV79e13" "IV79qkp" "IV79wsm" "LV80kod" "LV80srm" "MH84agw" "MH84g02" "MH84g03" "MH84gmr" "PM73phn" "PS86psa" "WN87cas" "WN87cat" "WN87flo" "WN87w70" "WN87wat" }
#small magnitude, large distance
set iGroundFileSMLR {"BO42elc" "CO83c05" "CO83c08" "IV79cc4" "IV79cmp" "IV79dlt" "IV79nil" "IV79pls" "IV79vct" "LV80stp" "MH84cap" "MH84hch" "MH84sjb" "PS86h06" "PS86ino" "WN87bir" "WN87cts" "WN87har" "WN87sse" "WN87stc" }

#large magnitude, small distance
set iGroundFileLMSR {"IV40elc" "LD92yer" "LP89agw" "LP89cap" "LP89g03" "LP89g04" "LP89gmr" "LP89hch" "LP89hda" "LP89svl" "NR94cnp" "NR94far" "NR94fle" "NR94glp" "NR94hol" "NR94stc" "SF71pel" "SH87bra" "SH87icc" "SH87wsm" }
#large magnitude, large distance
set iGroundFileLMLR {"BM68elc" "LD92ind" "LD92psa" "LP89a2e" "LP89fms" "LP89hvr" "LP89sjw" "LP89slc" "NR94bad" "NR94cas" "NR94cen" "NR94jab" "NR94lh1" "NR94loa" "NR94lv2" "NR94php" "NR94pic" "NR94sor" "NR94sse" "NR94ver" }
set iGroundFileLMLR {"SH87wsm" "BM68elc" "LD92ind" "LD92psa" "LP89a2e" "LP89fms" "LP89hvr" "LP89sjw" "LP89slc" "NR94bad" "NR94cas" "NR94cen" "NR94jab" "NR94lh1" "NR94loa" "NR94lv2" "NR94php" "NR94pic" "NR94sor" "NR94sse" "NR94ver" }
"
When I assume dt and run the program it displays that
"Warning- Pathseries :: Pathseries ()- could not open file GMFiles/......g3"

Would you also please upload a example ground motion file so that I'll confirm that my ground motion file is OK.

fawadnajam
Posts: 1
Joined: Sat Apr 07, 2012 5:56 am
Location: AIT

Re: Elastic Response Spectra Script

Post by fawadnajam » Wed Mar 26, 2014 12:24 am

How can this script be modified to get inelastic spectrum? I want to model elastoplastic and flag-shaped behaviors and want to get peak displacements of SDOFs.

seraphpan
Posts: 3
Joined: Wed Sep 11, 2013 8:18 pm
Location: UBC

Re: Elastic Response Spectra Script

Post by seraphpan » Mon Sep 28, 2015 2:25 pm

fawadnajam wrote:
> How can this script be modified to get inelastic spectrum? I want to model elastoplastic
> and flag-shaped behaviors and want to get peak displacements of SDOFs.

Hi Fawadnajam,

How's your work now? Have you successfully plot the inelastic spectra? I am trying to do similar stuff and kinda confused for a long time.

Thanks

Seraph

vantanhy
Posts: 1
Joined: Thu Oct 08, 2015 7:02 pm
Contact:

Re: Elastic Response Spectra Script

Post by vantanhy » Thu Oct 08, 2015 7:10 pm

Dear Silvia
About the damps I found large difference according to K-prop factor and M-prop factor. Which one should I use?
When T>1s acceleration values uing M-prop factor are larger when uing the other softwares, also larger than the elastic response spectrum. How to deal with this problem?

Thank you!

fmk
Site Admin
Posts: 5866
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley
Contact:

Re: Elastic Response Spectra Script

Post by fmk » Fri Oct 09, 2015 10:22 am

there will be large differences if you are doing inelastic response spectra and you specify mass proportional or stiffness proportional (current stiffness) damping .. using mass proportional the C is constant, using stifness proportional the C is changing .. hence the damping forces are different .. most inelastic spectrum would be determined using a constant C matrix, i.e. mass proportional (or initial stiffness)

Post Reply