Run Time History and Pushover analysis
Moderators: silvia, selimgunay, Moderators
-
- Posts: 2
- Joined: Fri Jan 21, 2011 7:07 am
- Location: Università degli studi di Padova
Run Time History and Pushover analysis
I have a problem: I can't write a program wich is able to run a Time History analysis and, after, a Pushover. The file filePath A10000.acc is the same I've found in the Examples Manual. The model built is the same for the analysis.
Here the tcl I've written:
#"Time history analysis"
wipe;
model BasicBuilder -ndm 2 -ndf 3;
#Definisco i due nodi d'estremità, con altezza della pila 10 m
node 1 0 0;
node 2 0 10000;
#Condizioni divincolo: incastro al nodo 1
fix 1 1 1 1;
# nodal masses:
mass 2 5.18 1.e-9 0.; # node#, Mx My Mz, Mass=Weight/g.
# definizione dei materiali
uniaxialMaterial ElasticPP 1 32e3 0.00175;
#Sezione circolare, d=1000 mm
# Circular Patch command
section fiberSec 10 {
patch circ 1 12 4 0.0 0.0 0.0 500 0 360
}
geomTransf PDelta 1
set NGauss 4;
element nonlinearBeamColumn 1 1 2 $NGauss 10 1 ;
# define GRAVITY -------------------------------------------------------------
timeSeries Linear 1
pattern Plain 1 1 {
load 2 0. -1e6 0.; # node#, FX FY MZ -- superstructure-weight
}
constraints Plain; # how it handles boundary conditions
numberer Plain; # renumber dof's to minimize band-width (optimization), if you want to
system BandGeneral; # how to store and solve the system of equations in the analysis
algorithm Linear; # use Linear algorithm for linear analysis
integrator LoadControl 0.1; # determine the next time step for an analysis, # apply gravity in 10 steps
analysis Static # define type of analysis static or transient
analyze 10; # perform gravity analysis
loadConst -time 0.0; # hold gravity constant and restart time
# DYNAMIC ground-motion analysis -------------------------------------------------------------
# create load pattern
set G 9810
timeSeries Path 2 -dt 0.005 -filePath A10000.acc $G;
pattern UniformExcitation 2 1 -accel 2;
set freq [expr [eigen 1]**0.5];
set dampRatio 0.02;
rayleigh 0. 0. 0. [expr 2*$dampRatio/$freq];
wipeAnalysis;
constraints Plain;
numberer Plain;
system BandGeneral;
algorithm Linear;
integrator Newmark 0.5 0.25;
analysis Transient;
analyze 3995 0.01;
puts "Done TH!"
This have no problem, then I try to run the following (without "wipe" command):
# RECORDER--------------------------------------------------------------------------------------------------
recorder Node -file RBase.txt -time -node 1 -dof 1 2 3 reaction; # support reaction
recorder Node -file Disp2.txt -time -node 2 -dof 1 2 3 disp; # dx,dy,rz nodo2
# Carico laterale push over F= 1e3 N -------------------------------------------------------------
# Lateral load pattern
pattern Plain 3 Linear {
load 2 1000 0 0 ;
}
# pushover: diplacement controlled static analysis
# Spostamento di 1 mm per 1000 step; spost totale D=500 mm
integrator DisplacementControl 2 1 1;
analyze 500;
puts "Done PO!"
But the program says me that there are problems with "Transient analysis".
If anyone could help me. Thank you.
Here the tcl I've written:
#"Time history analysis"
wipe;
model BasicBuilder -ndm 2 -ndf 3;
#Definisco i due nodi d'estremità, con altezza della pila 10 m
node 1 0 0;
node 2 0 10000;
#Condizioni divincolo: incastro al nodo 1
fix 1 1 1 1;
# nodal masses:
mass 2 5.18 1.e-9 0.; # node#, Mx My Mz, Mass=Weight/g.
# definizione dei materiali
uniaxialMaterial ElasticPP 1 32e3 0.00175;
#Sezione circolare, d=1000 mm
# Circular Patch command
section fiberSec 10 {
patch circ 1 12 4 0.0 0.0 0.0 500 0 360
}
geomTransf PDelta 1
set NGauss 4;
element nonlinearBeamColumn 1 1 2 $NGauss 10 1 ;
# define GRAVITY -------------------------------------------------------------
timeSeries Linear 1
pattern Plain 1 1 {
load 2 0. -1e6 0.; # node#, FX FY MZ -- superstructure-weight
}
constraints Plain; # how it handles boundary conditions
numberer Plain; # renumber dof's to minimize band-width (optimization), if you want to
system BandGeneral; # how to store and solve the system of equations in the analysis
algorithm Linear; # use Linear algorithm for linear analysis
integrator LoadControl 0.1; # determine the next time step for an analysis, # apply gravity in 10 steps
analysis Static # define type of analysis static or transient
analyze 10; # perform gravity analysis
loadConst -time 0.0; # hold gravity constant and restart time
# DYNAMIC ground-motion analysis -------------------------------------------------------------
# create load pattern
set G 9810
timeSeries Path 2 -dt 0.005 -filePath A10000.acc $G;
pattern UniformExcitation 2 1 -accel 2;
set freq [expr [eigen 1]**0.5];
set dampRatio 0.02;
rayleigh 0. 0. 0. [expr 2*$dampRatio/$freq];
wipeAnalysis;
constraints Plain;
numberer Plain;
system BandGeneral;
algorithm Linear;
integrator Newmark 0.5 0.25;
analysis Transient;
analyze 3995 0.01;
puts "Done TH!"
This have no problem, then I try to run the following (without "wipe" command):
# RECORDER--------------------------------------------------------------------------------------------------
recorder Node -file RBase.txt -time -node 1 -dof 1 2 3 reaction; # support reaction
recorder Node -file Disp2.txt -time -node 2 -dof 1 2 3 disp; # dx,dy,rz nodo2
# Carico laterale push over F= 1e3 N -------------------------------------------------------------
# Lateral load pattern
pattern Plain 3 Linear {
load 2 1000 0 0 ;
}
# pushover: diplacement controlled static analysis
# Spostamento di 1 mm per 1000 step; spost totale D=500 mm
integrator DisplacementControl 2 1 1;
analyze 500;
puts "Done PO!"
But the program says me that there are problems with "Transient analysis".
If anyone could help me. Thank you.
Re: Run Time History and Pushover analysis
you need to issue a "wipeAnalysis" command after the dynamic analysis and then build a static analysis for the pushover.
there is no space between wipe and Analysis, it's all one word.
there is no space between wipe and Analysis, it's all one word.
-
- Posts: 2
- Joined: Fri Jan 21, 2011 7:07 am
- Location: Università degli studi di Padova
Re: Run Time History and Pushover analysis
But I want to compare the Pushover analysis (Pushover Curve) before and after the Time History analysis, I want to conserve the results of the Time history and run the Pushover to obtain new and lower results. It is possible?
Re: Run Time History and Pushover analysis
i believe the wipeAnalysis only removes the analysis parameters, not the results and status
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Re: Run Time History and Pushover analysis
Hi Elena,
did you solve the problem?. I am also interested in comparing the performance of a structure before an earthquake (with a pushover analysis) and after that (time history and a subsequent pushover). Now I'll try to solve the problem. However, let me know; maybe we can compare the methods we implemented (I hope).
See you,
Andrea
did you solve the problem?. I am also interested in comparing the performance of a structure before an earthquake (with a pushover analysis) and after that (time history and a subsequent pushover). Now I'll try to solve the problem. However, let me know; maybe we can compare the methods we implemented (I hope).
See you,
Andrea
Scientists study the world as it is; engineers create the world that never has been.
Re: Run Time History and Pushover analysis
Hi Elena,
I tried an example (on a simple cantilever beam) and it works.
At the end of the tcl file with the specifications of time history analysis just put wipeAnalysis, and after that you can run a pushover analysis. In your output you will see first the results of TH and then those related to the static analysis.
Andrea
I tried an example (on a simple cantilever beam) and it works.
At the end of the tcl file with the specifications of time history analysis just put wipeAnalysis, and after that you can run a pushover analysis. In your output you will see first the results of TH and then those related to the static analysis.
Andrea
Scientists study the world as it is; engineers create the world that never has been.