Why the groundmotion script and rayleigh script are not work

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

Moderators: silvia, selimgunay, Moderators

Post Reply
YoungXu
Posts: 14
Joined: Wed May 21, 2014 2:12 am
Location: University of XinJiang

Why the groundmotion script and rayleigh script are not work

Post by YoungXu »

Today, I found that the groundmotion and rayleigh damping are not work, for example, I enter the Sine wave and the rayleigh for the cantilever, but the cantilever's displacement plot always sine wave. So, I want to know why? Thanks everybody!

#units: mm, N, MPa, ton
model BasicBuilder -ndm 2 -ndf 3
file mkdir Cantilever
set h 180
set b 170
set d 4
set t 5
set E 200000 ;#200GPa=200,000MPa=2*10^11N/m^2
set Fy 235 ;#235MPa=2.35*10^8N/m^2
set L 1800
set A [expr $h*$b-($b-$d)*($h-2*$t)]
set I [expr ($b*pow($h,3)-($b-$d)*pow($h-2*$t,3))/12]
set density [expr 7.8*pow(10,-6)]
node 1 0 0
node 2 0 $L
fix 1 1 1 1
mass 2 [expr $L*$md*$A/1000] [expr $L*$density*$A/1000] 0
geomTransf Linear 1
uniaxialMaterial Steel01 1 235 200000 0.01 0 1 0 1
element elasticBeamColumn 1 1 2 $A $E $I 1

recorder Node -file Cantilever/Node2_Dsp.out -time -node 2 -dof 1 disp
recorder Node -file Cantilever/Node2_Vel.out -time -node 2 -dof 1 vel
recorder Node -file Cantilever/Node2_Acc.out -time -node 2 -dof 1 accel
recorder Element -file Cantilever/Elmt1_Frc.out -time -ele 1 forces

rayleigh 7.9 [expr 1.3*pow(10,-5)] 0 0

set dt 0.02
set scale [expr 1]
timeSeries Path 1 -filePath Sine.txt -dt $dt -factor [expr $scale]
pattern UniformExcitation 1 1 -accel 1

system BandGeneral
numberer RCM
constraints Transformation
test EnergyIncr 1.0e-20 10
integrator AlphaOS 0.67
algorithm Linear
analysis Transient
analyze 1000 0.02
wipe
exit
Post Reply