TimeSeries UniformExcitation

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

Moderators: silvia, selimgunay, Moderators

Post Reply
talledo.diego
Posts: 20
Joined: Sat Jun 27, 2009 2:59 am
Location: IUAV - Venice

TimeSeries UniformExcitation

Post by talledo.diego »

I'm looking example "Ex1a.Canti2D.EQ.tcl" and also Getting Started.
There is an error in:

[code]
# create load pattern
source ReadSMDFile.tcl
ReadSMDFile BM68elc.th BM68elc.acc dt
set accelSeries "Series -dt $dt -filePath BM68elc.acc -factor 1";
pattern UniformExcitation 2 1 -accel $accelSeries[/code]

I tried also to put each line into OpenSees interpreter:

[code]
OpenSees > wipe
OpenSees > source example.tcl
3
OpenSees > source GravityLoads.tcl
OpenSees > source ReadSMDFile.tcl
OpenSees > ReadSMDFile BM68elc.th BM68elc.acc dt
OpenSees > put $dt
.01000
OpenSees > set accelSeries "Series -dt 0.01 -filePath BM68elc.acc -factor 1";
Series -dt 0.01 -filePath BM68elc.acc -factor 1
OpenSees > pattern UniformExcitation 2 1 -accel $accelSeries;
expected integer but got "-dt"
OpenSees >
[\code]

I think the error is that interpreter expects the Tag of TimeSeries after pattern command. But it is strange because with pattern command I can use "Linear" time serie without errors.

What do you think about it?

Regards
Diego Alejandro Talledo
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Post by vesna »

J just ran this file and it works fine. Save the file Ex1a.Canti2D.EQ.tcl and the ground motion file in the same folder. Let me know if you still have problems with running this file.
talledo.diego
Posts: 20
Joined: Sat Jun 27, 2009 2:59 am
Location: IUAV - Venice

Post by talledo.diego »

File Ex1a.Canti2D.EQ.tcl works fine. And this is very strange for me: this file gives as output

"OpenSees > source GS_Dyn.tcl
expected integer but got "500.0"
OpenSees >
"

My file (GS_Dyn.tcl):
[code]
model basic -ndm 2 -ndf 3
# nodal coordinates:
node 1 0 0
node 2 504 0
node 3 0 432
node 4 504 432

# bondary conditions:
fix 1 1 1 1
fix 2 1 1 1
fix 3 0 0 0
fix 4 0 0 0

# nodal masses:
mass 3 5.18 0. 0.
mass 4 5.18 0. 0.

# transformation:
geomTransf Linear 1

# element connectivity:
element elasticBeamColumn 1 1 3 3600 4227 1080000 1
element elasticBeamColumn 2 2 4 3600 4227 1080000 1
element elasticBeamColumn 3 3 4 5760 4227 4423680 1

# recorders
recorder Node -file Node3.out -time -node 3 -dof 1 2 disp
recorder Element -file Element1.out -time -ele 1 force

#Gravity loads
pattern Plain 1 Linear {
load 3 0.0 -2000 -168074
load 4 0.0 -2000 168074
}

constraints Transformation
numberer RCM
system BandGeneral
test NormDispIncr 1.0e-6 6
algorithm Newton
integrator LoadControl 0.1
analysis Static

analyze 10
loadConst -time 0.0

# Dynamic analysis
# create load pattern
set accelSeries "Series -dt 0.01 -filePath BM68elc.acc -factor 1";
pattern UniformExcitation 2 1 -accel $accelSeries

# set damping factors
rayleigh 0. 0. 0. [expr 2*0.02/pow([eigen 1],0.5)]

# create the analysis
wipeAnalysis
constraints Plain
numberer RCM
system UmfPack
test NormDispIncr 1.0e-8 10
algorithm Newton
integrator Newmark 0.5 0.25
analysis Transient

analyze [expr 10/0.02] 0.02
[/code]

I really don't understand what is happening.
Another thing I just discovered. If I try to give by hand all commands of file Ex1a.Canti2D.EQ.tcl it gives same error I posted in my previous message (expected integer but got "-dt" )

I'm using right now version 2.2.1
I will try tomorrow with older versions.

thanks!
Diego Alejandro Talledo
fmk
Site Admin
Posts: 5884
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley
Contact:

Post by fmk »

10/0.02 gives 500.0 .. the analyze command requires an integer .. you need to convert the number to an integer .. use
[expr int(10/0.02)] instead.
talledo.diego
Posts: 20
Joined: Sat Jun 27, 2009 2:59 am
Location: IUAV - Venice

Post by talledo.diego »

Thank you Frank. The error was this.
Sorry, but when I tried inputing all lines manually, when I saw this, I thought it was an error. But it's not, becuase if I go ahead the analysis is ok. I think it is a warning, is it right?

[code]
OpenSees > model basic -ndm 2 -ndf 3
OpenSees > node 1 0 0
OpenSees > node 2 504 0
OpenSees > node 3 0 432
OpenSees > node 4 504 432
OpenSees > fix 1 1 1 1
OpenSees > fix 2 1 1 1
OpenSees > fix 3 0 0 0
OpenSees > fix 4 0 0 0
OpenSees > mass 3 5.18 0. 0.
OpenSees > mass 4 5.18 0. 0.
OpenSees > geomTransf Linear 1
OpenSees > element elasticBeamColumn 1 1 3 3600 4227 1080000 1
OpenSees > element elasticBeamColumn 2 2 4 3600 4227 1080000 1
OpenSees > element elasticBeamColumn 3 3 4 5760 4227 4423680 1
OpenSees > pattern Plain 1 Linear {
load 3 0.0 -2000 -168074
load 4 0.0 -2000 168074
}
OpenSees > constraints Transformation
OpenSees > numberer RCM
OpenSees > system BandGeneral
OpenSees > test NormDispIncr 1.0e-6 6
OpenSees > algorithm Newton
OpenSees > integrator LoadControl 0.1
OpenSees > analysis Static
OpenSees > analyze 10
0
OpenSees > loadConst -time 0.0
OpenSees > set accelSeries "Series -dt 0.01 -filePath BM68elc.acc -factor 1"
Series -dt 0.01 -filePath BM68elc.acc -factor 1
OpenSees > pattern UniformExcitation 2 1 -accel $accelSeries
expected integer but got "-dt"
OpenSees >
[/code]
Diego Alejandro Talledo
Post Reply