IDA analysis

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

Moderators: silvia, selimgunay, Moderators

Post Reply
hongqi
Posts: 16
Joined: Tue Jul 22, 2008 12:23 am

IDA analysis

Post by hongqi » Thu Jun 27, 2013 6:59 pm

hi ,i wrote the code of IDA,i meet two problem ,1) in the loop , when $i=2 the ground start time is 28.01,but it should be 0.01; 2) the root displacement is different the $gmf=2.998 which wipe the loop, can you tell me what's wrong with the code?
#“IDA CODE”
#Units: N,m .Kg
#Define model builder
file mkdir ida
puts "modle built]"
model basic -ndm 2 -ndf 3
#---------------------------------------
node 1 0.0 0.0
node 2 0.0 1.0
node 3 0.0 2.0
node 4 0.0 3.0
node 5 0.0 4.0
node 6 0.0 5.0
node 7 0.0 6.0
node 8 0.0 7.0
node 9 0.0 8.0
node 10 0.0 9.0
node 11 0.0 10.0
fix 1 1 1 1 ;
set p 5000000.0; #--DEAD LOAD--------
puts "mass"
mass 1 4044.7 1e-9 0.0
mass 2 8009.4 1e-9 0.0
mass 3 8009.4 1e-9 0.0
mass 4 8009.4 1e-9 0.0
mass 5 8009.4 1e-9 0.0
mass 6 8009.4 1e-9 0.0
mass 7 8009.4 1e-9 0.0
mass 8 8009.4 1e-9 0.0
mass 9 8009.4 1e-9 0.0
mass 10 8009.4 1e-9 0.0
mass 11 513896.59 1e-9 0.0

#Cover concrete
uniaxialMaterial Concrete02 1 -25000000.0 -0.002 -5000000.0 -0.004 0.05 2500000.0 3125000.0
#Core concrete
uniaxialMaterial Concrete02 2 -26819372.69 -0.0027 -56891687.56 -0.016 0.05 2500000.0 3125000.0
#Steel model
uniaxialMaterial Steel02 3 388000000.0 2.0e11 0.01 18.5 0.925 0.15

#Define cross-section for nonlinear column
#----------------------------------------
set cover 0.05
set D 2.00
set d2 0.036
set r [expr $D/2]
set r1 [expr $r-$cover]
set r2 [expr $r1-$d2/2]
#Column section
section Fiber 1 {
patch circ 2 40 40 0.0 0.0 0.0 $r1 0.0 360.0
patch circ 1 40 40 0.0 0.0 $r1 $r 0.0 360.0
layer circ 3 75 101.7e-5 0.0 0.0 $r2 0.0 360.0
}
geomTransf Linear 1
set np 5
#element nonlinearBeamColumn 1 1 2 $np 1 1
element dispBeamColumn 1 1 2 $np 1 1
element dispBeamColumn 2 2 3 $np 1 1
element dispBeamColumn 3 3 4 $np 1 1
element dispBeamColumn 4 4 5 $np 1 1
element dispBeamColumn 5 5 6 $np 1 1
element dispBeamColumn 6 6 7 $np 1 1
element dispBeamColumn 7 7 8 $np 1 1
element dispBeamColumn 8 8 9 $np 1 1
element dispBeamColumn 9 9 10 $np 1 1
element dispBeamColumn 10 10 11 $np 1 1

#---------gravity analysis--------------------------
pattern Plain 1 "Linear" {
load 11 0.0 [expr -$p] 0.0
}
#Define analysis parameters
integrator LoadControl 0.1
system BandGeneral
test NormDispIncr 1.0e-6 20 ;
numberer Plain
constraints Plain
algorithm Newton
# algorithm NewtonLineSearch 0.75;
analysis Static
#Do one analysis for constant axial load
recorder Node -file ida/reaction1.out -time -node 1 -dof 1 2 3 reaction
analyze 10
puts "Gravity load analysis completed"
loadConst -time 0.0
#--------time-history -----------------------------
set xDamp 0.05;
set nEigenI 1;
set nEigenJ 2;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM [expr $xDamp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaKcurr [expr 2.*$xDamp/($omegaI+$omegaJ)];
rayleigh $alphaM $betaKcurr 0 0
#-----------factor----------------------
for {set i 1} {$i<=2} {incr i} {
set sa 0.208;
set gm [expr 0.02+0.6036*($i-1)] ;
set gmf [expr $gm/$sa];
puts "gmf=[expr $gmf]"
file mkdir x[expr $gmf] ;
puts *****************************
puts "groundmotion start time:[getTime]"

puts "recorder"
recorder Node -file x[expr $gmf]/X_11.out -time -node 11 -dof 1 disp
recorder Element -file x[expr $gmf]/jbele_1.out -time -ele 1 localForce
recorder Element -file x[expr $gmf]/ztele_1.out -time -ele 1 globalForce
set IDloadTag [expr $i+100];
set iGMfile "GM1X.txt";
set iGMdirection "1";
set iGMfactor [expr ($gmf)*9.806]
set TmaxAnalysis 28.0;
set DtAnalysis 0.01;
set AccelSeries "Series -dt $DtAnalysis -filePath $iGMfile -factor $iGMfactor";
puts "$AccelSeries"
pattern UniformExcitation $IDloadTag $iGMdirection -accel $AccelSeries;

variable constraintsTypeDynamic Transformation;
constraints $constraintsTypeDynamic ;
variable numbererTypeDynamic RCM
numberer $numbererTypeDynamic;
variable systemTypeDynamic UmfPack;
system $systemTypeDynamic;
variable maxNumIterDynamic 400;
variable TolDynamic 1.e-3
variable testTypeDynamic EnergyIncr
test $testTypeDynamic $TolDynamic $maxNumIterDynamic
variable algorithmTypeDynamic Newton;
algorithm $algorithmTypeDynamic
variable NewmarkGamma 0.5
variable NewmarkBeta 0.25
variable integratorTypeDynamic Newmark;
integrator $integratorTypeDynamic $NewmarkGamma $NewmarkBeta
variable analysisTypeDynamic VariableTransient
analysis $analysisTypeDynamic
set steps [expr int($TmaxAnalysis/$DtAnalysis)];
set ok [analyze $steps $DtAnalysis];
if {$ok != 0} { ; # analysis was not successful.
# ----------------------------------------------------
set ok 0;
set controlTime [getTime];
while {$controlTime < $TmaxAnalysis && $ok == 0} {
set controlTime [getTime]
set ok [analyze 1 $DtAnalysis]
if {$ok != 0} {
puts "Trying Newton with Initial Tangent .."
test NormDispIncr $TolDynamic 1000 0
algorithm Newton -initial
set ok [analyze 1 $DtAnalysis]
test $testTypeDynamic $TolDynamic $maxNumIterDynamic 0
algorithm $algorithmTypeDynamic
}
if {$ok != 0} {
puts "Trying Broyden .."
algorithm Broyden 8
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmTypeDynamic
}
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch .6
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmTypeDynamic
}
if {$ok != 0} {
puts "Trying ModifiedNewton .."
algorithm ModifiedNewton
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmTypeDynamic
}
if {$ok != 0} {
puts "Trying BFGS.."
algorithm BFGS
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmTypeDynamic
}
}
};

puts "Ground Motion DoneX[expr $gmf]. End Time: [getTime]"
puts ****************************************************************
}

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

Re: IDA analysis

Post by fmk » Sun Jun 30, 2013 12:52 pm

the second analysis starts at the end of the first. if you want each to start from 0 with the model at the initial state you need to move the for loop to the start of the file and put in a wipe.

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

Re: IDA analysis

Post by fmk » Sun Jun 30, 2013 12:52 pm

the second analysis starts at the end of the first. if you want each to start from 0 with the model at the initial state you need to move the for loop to the start of the file and put in a wipe.

hongqi
Posts: 16
Joined: Tue Jul 22, 2008 12:23 am

Re: IDA analysis

Post by hongqi » Sun Jun 30, 2013 5:49 pm

to fmk:
thanks,i move the loop to the start of the code ,the result is equal.




fmk wrote:
> the second analysis starts at the end of the first. if you want each to
> start from 0 with the model at the initial state you need to move the for
> loop to the start of the file and put in a wipe.

hongqi
Posts: 16
Joined: Tue Jul 22, 2008 12:23 am

Re: IDA analysis

Post by hongqi » Sun Jun 30, 2013 6:10 pm

to fmk :
i have an question about reyleigh damping factors,most exeamples show that the value of $betaKcurr and $betaKinit is 0,so i want to know how to define $betaKcurr and $betaKinit ,and What is the difference in using the $betaKcurr ,$betaKinit and $betaKcomm.


fmk wrote:
> the second analysis starts at the end of the first. if you want each to
> start from 0 with the model at the initial state you need to move the for
> loop to the start of the file and put in a wipe.

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

Re: IDA analysis

Post by fmk » Mon Jul 01, 2013 9:45 pm

these are the factor applied to the current tangent, the last committed tangent and the initial tangent. most programs use the factor on the initial tangent.

hongqi
Posts: 16
Joined: Tue Jul 22, 2008 12:23 am

Re: IDA analysis

Post by hongqi » Wed Jul 03, 2013 4:29 am

to fmk
thank you very much .when i check out the results and encountered a problem .
there are two recorder of support reaction in my code ,one is " recorder Node -file ida/reaction1.out -time -node 1 -dof 1 2 3 reaction ",the other is " recorder Element -file x[expr $gmf]/ztele_1.out -time -ele 1 globalForce" , at the same dt increment ,i find the shear of node 1 is different in recorders ,but they should be equal ,maybe the mass at roof is larger than orther nodes?


fmk wrote:
> these are the factor applied to the current tangent, the last committed
> tangent and the initial tangent. most programs use the factor on the
> initial tangent.

punitbarwani
Posts: 19
Joined: Tue Aug 06, 2013 10:48 pm
Location: Indore, India

Comparision Of Opensees And IDARC

Post by punitbarwani » Wed Aug 07, 2013 12:15 am

# --------------------------------------------------------------------------------------------------
# Example4. 2D Portal Frame-- Dynamic EQ input analysis
#
# execute this file after you have built the model, and after you apply gravity
#

# Uniform Earthquake ground motion (uniform acceleration input at all support nodes)
set GMdirection 1; # ground-motion direction
set GMfile "H-E12140" ; # ground-motion filenames
set GMfact 0.005; # ground-motion scaling factor

# set up ground-motion-analysis parameters
set DtAnalysis [expr 1.0*$sec]; # time-step Dt for lateral analysis
set TmaxAnalysis [expr 40.0*$sec]; # maximum duration of ground-motion analysis -- should be 50*$sec

# ----------- set up analysis parameters
source LibAnalysisDynamicParameters.tcl; # constraintsHandler,DOFnumberer,system-ofequations,convergenceTest,solutionAlgorithm,integrator

# define DAMPING--------------------------------------------------------------------------------------
# apply Rayleigh DAMPING from $xDamp
# D=$alphaM*M + $betaKcurr*Kcurrent + $betaKcomm*KlastCommit + $beatKinit*$Kinitial
set xDamp 0.05; # 2% damping ratio
set lambda [eigen 1]; # eigenvalue mode 1
set omega [expr pow($lambda,0.5)];
set alphaM 0.; # M-prop. damping; D = alphaM*M
set betaKcurr 0.; # K-proportional damping; +beatKcurr*KCurrent
set betaKcomm [expr 2.*$xDamp/($omega)]; # K-prop. damping parameter; +betaKcomm*KlastCommitt
set betaKinit 0.; # initial-stiffness proportional damping +beatKinit*Kini
rayleigh $alphaM $betaKcurr $betaKinit $betaKcomm; # RAYLEIGH damping

# --------------------------------- perform Dynamic Ground-Motion Analysis
# the following commands are unique to the Uniform Earthquake excitation
set IDloadTag 400; # for uniformSupport excitation
# read a PEER strong motion database file, extracts dt from the header and converts the file
# to the format OpenSees expects for Uniform/multiple-support ground motions
source ReadSMDFile.tcl; # read in procedure Multinition
# Uniform EXCITATION: acceleration input
set inFile $GMdir/$GMfile.at2
set outFile $GMdir/$GMfile.g3; # set variable holding new filename (PEER files have .at2/dt2 extension)
ReadSMDFile $inFile $outFile dt; # call procedure to convert the ground-motion file
set GMfatt [expr $g*$GMfact]; # data in input file is in g Unifts -- ACCELERATION TH
set AccelSeries "Series -dt $dt -filePath $outFile -factor $GMfatt"; # time series information
pattern UniformExcitation $IDloadTag $GMdirection -accel $AccelSeries ; # create Unifform excitation

set Nsteps [expr int($TmaxAnalysis/$DtAnalysis)];
set ok [analyze $Nsteps $DtAnalysis]; # actually perform analysis; returns ok=0 if analysis was successful

if {$ok != 0} { ; # analysis was not successful.
# --------------------------------------------------------------------------------------------------
# change some analysis parameters to achieve convergence
# performance is slower inside this loop
# Time-controlled analysis
set ok 0;
set controlTime [getTime];
while {$controlTime < $TmaxAnalysis && $ok == 0} {
set controlTime [getTime]
set ok [analyze 1 $DtAnalysis]
if {$ok != 0} {
puts "Trying Newton with Initial Tangent .."
test NormDispIncr $Tol 1000 0
algorithm Newton -initial
set ok [analyze 1 $DtAnalysis]
test $testTypeDynamic $TolDynamic $maxNumIterDynamic 0
algorithm $algorithmTypeDynamic
}
if {$ok != 0} {
puts "Trying Broyden .."
algorithm Broyden 8
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmTypeDynamic
}
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch .8
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmTypeDynamic
}
}
}; # end if ok !0


puts "Ground Motion Done. End Time: [getTime]"
##############################################################################



file 2:-
--------------------------------------------------------------------------------------------------
# Example4. 2D Portal Frame-- Build Model
# nonlinearBeamColumn element, inelastic fiber section
#
#
# ^Y
# |
# 3____________(3)___________4 __
# | | |
# | | |
# | | |
# (1) (2) LCol
# | | |
# | | |
# | | |
# =1= =2= _|_ -------->X
# |----------LBeam-----------|
#

# SET UP ----------------------------------------------------------------------------
wipe; # clear memory of all past model definitions
model BasicBuilder -ndm 2 -ndf 3; # Define the model builder, ndm=#dimension, ndf=#dofs
set dataDir Data; # set up name of data directory
file mkdir $dataDir; # create data directory
set GMdir "GMfiles"; # ground-motion file directory
source LibUnits.tcl; # define basic and system units

# define GEOMETRY -------------------------------------------------------------
set LCol [expr 6000*$mm]; # column length
set LBeam [expr 6000*$mm]; # beam length
set Weight [expr 5*$N/$mm]; # superstructure weight
# define section geometry
set HCol [expr 400.*$mm]; # Column Depth
set BCol [expr 400.*$mm]; # Column Width
set HBeam [expr 500.*$mm]; # Beam Depth
set BBeam [expr 300.*$mm]; # Beam Width

# calculated parameters
set PCol [expr $Weight]; # nodal dead-load weight per column
set Mass [expr $PCol/$g]; # nodal mass
set MCol [expr 1./12.*($Weight/$LBeam)*pow($LBeam,2)]; # beam-end moment due to distributed load.
# calculated geometry parameters
set ACol [expr $BCol*$HCol]; # cross-sectional area
set ABeam [expr $BBeam*$HBeam];
set IzCol [expr 1./12.*$BCol*pow($HCol,3)]; # Column moment of inertia
set IzBeam [expr 1./12.*$BBeam*pow($HBeam,3)]; # Beam moment of inertia

# nodal coordinates:
node 1 0 0; # node#, X, Y
node 2 $LBeam 0
node 3 0 $LCol
node 4 $LBeam $LCol

# Single point constraints -- Boundary Conditions
fix 1 1 1 0; # node DX DY RZ
fix 2 1 1 0; # node DX DY RZ
#fix 3 0 0 1
#fix 4 0 0 1

# nodal masses:
mass 3 $Mass 0. 0.; # node#, Mx My Mz, Mass=Weight/g, neglect rotational inertia at nodes
mass 4 $Mass 0. 0.

# Define ELEMENTS & SECTIONS -------------------------------------------------------------
set ColSecTag 1; # assign a tag number to the column section
set BeamSecTag 2; # assign a tag number to the beam section
# define section geometry
set coverCol [expr 40.*$mm]; # Column cover to reinforcing steel NA.
set numBarsCol 4; # number of longitudinal-reinforcement bars in each side of column section. (symmetric top & bot)
set barAreaCol [expr 1256*$mm*$mm]; # area of longitudinal-reinforcement bars
#set coverBeam [expr 40.*$mm]; # Column cover to reinforcing steel NA.
#set numBarsBeam 8; # number of longitudinal-reinforcement bars in each side of column section. (symmetric top & bot)
#set barAreaBeam [expr 2512*$mm*$mm]; # area of longitudinal-reinforcement bars

# MATERIAL parameters -------------------------------------------------------------------
set IDconcU 1; # material ID tag -- unconfined cover concrete
set IDreinf 2; # material ID tag -- reinforcement
# nominal concrete compressive strength
set fc [expr -25*$MPa]; # CONCRETE Compressive Strength, ksi (+Tension, -Compression)
set Ec [expr 5700*$MPa*sqrt(-$fc)]; # Concrete Elastic Modulus
# unconfined concrete
set fc1U $fc; # UNCONFINED concrete (todeschini parabolic model), maximum stress
set eps1U -0.003; # strain at maximum strength of unconfined concrete
set fc2U [expr 0.2*$fc1U]; # ultimate stress
set eps2U -0.05; # strain at ultimate stress
set lambda 0.1; # ratio between unloading slope at $eps2 and initial slope $Ec
# tensile-strength properties
set ftU [expr -0.14*$fc1U]; # tensile strength +tension
set Ets [expr $ftU/0.002]; # tension softening stiffness
# -----------
set Fy [expr 415*$MPa]; # STEEL yield stress
set Es [expr 200000.*$MPa]; # modulus of steel
set Bs 0.01; # strain-hardening ratio
set R0 18; # control the transition from elastic to plastic branches
set cR1 0.925; # control the transition from elastic to plastic branches
set cR2 0.15; # control the transition from elastic to plastic branches
uniaxialMaterial Concrete02 $IDconcU $fc1U $eps1U $fc2U $eps2U $lambda $ftU $Ets; # build cover concrete (unconfined)
uniaxialMaterial Steel02 $IDreinf $Fy $Es $Bs $R0 $cR1 $cR2; # build reinforcement material

# FIBER SECTION properties -------------------------------------------------------------
# symmetric section
# y
# ^
# |
# --------------------- -- --
# | o o o | | -- cover
# | | |
# | | |
# z <--- | + | H
# | | |
# | | |
# | o o o | | -- cover
# --------------------- -- --
# |-------- B --------|
#
# RC section:
set coverY [expr $HCol/2.0]; # The distance from the section z-axis to the edge of the cover concrete -- outer edge of cover concrete
set coverZ [expr $BCol/2.0]; # The distance from the section y-axis to the edge of the cover concrete -- outer edge of cover concrete
set coreY [expr $coverY-$coverCol]
set coreZ [expr $coverZ-$coverCol]
set nfY 16; # number of fibers for concrete in y-direction
set nfZ 16; # number of fibers for concrete in z-direction
section fiberSec $ColSecTag {; # Define the fiber section
patch quadr $IDconcU $nfZ $nfY -$coverY $coverZ -$coverY -$coverZ $coverY -$coverZ $coverY $coverZ; # Define the concrete patch
layer straight $IDreinf $numBarsCol $barAreaCol -$coreY $coreZ -$coreY -$coreZ; # top layer reinfocement
layer straight $IDreinf $numBarsCol $barAreaCol $coreY $coreZ $coreY -$coreZ; # bottom layer reinforcement
}; # end of fibersection definition

# BEAM section:
section Elastic $BeamSecTag $Ec $ABeam $IzBeam; # elastic beam section

# define geometric transformation: performs a linear geometric transformation of beam stiffness and resisting force from the basic system to the global-coordinate system
set ColTransfTag 1; # associate a tag to column transformation
set BeamTransfTag 2; # associate a tag to beam transformation (good practice to keep col and beam separate)
set ColTransfType Linear ; # options, Linear PDelta Corotational
geomTransf $ColTransfType $ColTransfTag ; # only columns can have PDelta effects (gravity effects)
geomTransf Linear $BeamTransfTag ;

# element connectivity:
set numIntgrPts 5; # number of integration points for force-based element
element nonlinearBeamColumn 1 1 3 $numIntgrPts $ColSecTag $ColTransfTag; # self-explanatory when using variables
element nonlinearBeamColumn 2 2 4 $numIntgrPts $ColSecTag $ColTransfTag;
element nonlinearBeamColumn 3 3 4 $numIntgrPts $BeamSecTag $BeamTransfTag;

# Define RECORDERS -------------------------------------------------------------
recorder Node -file $dataDir/DFree.out -time -node 3 4 -dof 1 2 3 disp; # displacements of free nodes
recorder Node -file $dataDir/DBase.out -time -node 1 2 -dof 1 2 3 disp; # displacements of support nodes
recorder Node -file $dataDir/RBase.out -time -node 1 2 -dof 1 2 3 reaction; # support reaction
recorder Drift -file $dataDir/Drift.out -time -iNode 1 2 -jNode 3 4 -dof 1 -perpDirn 2 ; # lateral drift
recorder Element -file $dataDir/FCol.out -time -ele 1 2 globalForce; # element forces -- column
recorder Element -file $dataDir/FBeam.out -time -ele 3 globalForce; # element forces -- beam
recorder Element -file $dataDir/ForceColSec1.out -time -ele 1 2 section 1 force; # Column section forces, axial and moment, node i
recorder Element -file $dataDir/DefoColSec1.out -time -ele 1 2 section 1 deformation; # section deformations, axial and curvature, node i
recorder Element -file $dataDir/ForceColSec$numIntgrPts.out -time -ele 1 2 section $numIntgrPts force; # section forces, axial and moment, node j
recorder Element -file $dataDir/DefoColSec$numIntgrPts.out -time -ele 1 2 section $numIntgrPts deformation; # section deformations, axial and curvature, node j
recorder Element -file $dataDir/ForceBeamSec1.out -time -ele 3 section 1 force; # Beam section forces, axial and moment, node i
recorder Element -file $dataDir/DefoBeamSec1.out -time -ele 3 section 1 deformation; # section deformations, axial and curvature, node i
recorder Element -file $dataDir/ForceBeamSec$numIntgrPts.out -time -ele 3 section $numIntgrPts force; # section forces, axial and moment, node j
recorder Element -file $dataDir/DefoBeamSec$numIntgrPts.out -time -ele 3 section $numIntgrPts deformation; # section deformations, axial and curvature, node j

# define GRAVITY -------------------------------------------------------------
set WzBeam [expr $Weight/$LBeam];
timeSeries Linear 1
pattern Plain 1 Linear {
eleLoad -ele 3 -type -beamUniform -$WzBeam ; # distributed superstructure-weight on beam
}
# Gravity-analysis parameters -- load-controlled static analysis
set Tol 1.0e-8; # convergence tolerance for test
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
test NormDispIncr $Tol 6 ; # determine if convergence has been achieved at the end of an iteration step
algorithm Newton; # use Newton's solution algorithm: updates tangent stiffness at every iteration
set NstepGravity 10; # apply gravity in 10 steps
set DGravity [expr 1./$NstepGravity]; # first load increment;
integrator LoadControl $DGravity; # determine the next time step for an analysis
analysis Static; # define type of analysis static or transient
analyze $NstepGravity; # apply gravity
# ------------------------------------------------- maintain constant gravity loads and reset time to zero
loadConst -time 0.0

puts "Model Built"
#######################################################################################




file 3:-
# --------------------------------------------------------------------------------------------------
# dynamic-analysis parameters
# I am setting all these variables as global variables (using variable rather than set command)
# so that these variables can be uploaded by a procedure
#


# Set up Analysis Parameters ---------------------------------------------
# CONSTRAINTS handler -- Determines how the constraint equations are enforced in the analysis (http://opensees.berkeley.edu/OpenSees/m ... al/617.htm)
# Plain Constraints -- Removes constrained degrees of freedom from the system of equations
# Lagrange Multipliers -- Uses the method of Lagrange multipliers to enforce constraints
# Penalty Method -- Uses penalty numbers to enforce constraints
# Transformation Method -- Performs a condensation of constrained degrees of freedom
variable constraintsTypeDynamic Transformation;
constraints $constraintsTypeDynamic ;

# DOF NUMBERER (number the degrees of freedom in the domain): (http://opensees.berkeley.edu/OpenSees/m ... al/366.htm)
# determines the mapping between equation numbers and degrees-of-freedom
# Plain -- Uses the numbering provided by the user
# RCM -- Renumbers the DOF to minimize the matrix band-width using the Reverse Cuthill-McKee algorithm
variable numbererTypeDynamic RCM
numberer $numbererTypeDynamic

# SYSTEM (http://opensees.berkeley.edu/OpenSees/m ... al/371.htm)
# Linear Equation Solvers (how to store and solve the system of equations in the analysis)
# -- provide the solution of the linear system of equations Ku = P. Each solver is tailored to a specific matrix topology.
# ProfileSPD -- Direct profile solver for symmetric positive definite matrices
# BandGeneral -- Direct solver for banded unsymmetric matrices
# BandSPD -- Direct solver for banded symmetric positive definite matrices
# SparseGeneral -- Direct solver for unsymmetric sparse matrices (-piv option)
# SparseSPD -- Direct solver for symmetric sparse matrices
# UmfPack -- Direct UmfPack solver for unsymmetric matrices
variable systemTypeDynamic BandGeneral; # try UmfPack for large problems
system $systemTypeDynamic

# TEST: # convergence test to
# Convergence TEST (http://opensees.berkeley.edu/OpenSees/m ... al/360.htm)
# -- Accept the current state of the domain as being on the converged solution path
# -- determine if convergence has been achieved at the end of an iteration step
# NormUnbalance -- Specifies a tolerance on the norm of the unbalanced load at the current iteration
# NormDispIncr -- Specifies a tolerance on the norm of the displacement increments at the current iteration
# EnergyIncr-- Specifies a tolerance on the inner product of the unbalanced load and displacement increments at the current iteration
# RelativeNormUnbalance --
# RelativeNormDispIncr --
# RelativeEnergyIncr --
variable TolDynamic 1.e-8; # Convergence Test: tolerance
variable maxNumIterDynamic 10; # Convergence Test: maximum number of iterations that will be performed before "failure to converge" is returned
variable printFlagDynamic 0; # Convergence Test: flag used to print information on convergence (optional) # 1: print information on each step;
variable testTypeDynamic EnergyIncr; # Convergence-test type
test $testTypeDynamic $TolDynamic $maxNumIterDynamic $printFlagDynamic;
# for improved-convergence procedure:
variable maxNumIterConvergeDynamic 2000;
variable printFlagConvergeDynamic 0;

# Solution ALGORITHM: -- Iterate from the last time step to the current (http://opensees.berkeley.edu/OpenSees/m ... al/682.htm)
# Linear -- Uses the solution at the first iteration and continues
# Newton -- Uses the tangent at the current iteration to iterate to convergence
# ModifiedNewton -- Uses the tangent at the first iteration to iterate to convergence
# NewtonLineSearch --
# KrylovNewton --
# BFGS --
# Broyden --
variable algorithmTypeDynamic ModifiedNewton
algorithm $algorithmTypeDynamic;

# Static INTEGRATOR: -- determine the next time step for an analysis (http://opensees.berkeley.edu/OpenSees/m ... al/689.htm)
# LoadControl -- Specifies the incremental load factor to be applied to the loads in the domain
# DisplacementControl -- Specifies the incremental displacement at a specified DOF in the domain
# Minimum Unbalanced Displacement Norm -- Specifies the incremental load factor such that the residual displacement norm in minimized
# Arc Length -- Specifies the incremental arc-length of the load-displacement path
# Transient INTEGRATOR: -- determine the next time step for an analysis including inertial effects
# Newmark -- The two parameter time-stepping method developed by Newmark
# HHT -- The three parameter Hilbert-Hughes-Taylor time-stepping method
# Central Difference -- Approximates velocity and acceleration by centered finite differences of displacement
variable NewmarkGamma 0.5; # Newmark-integrator gamma parameter (also HHT)
variable NewmarkBeta 0.25; # Newmark-integrator beta parameter
variable integratorTypeDynamic Newmark;
integrator $integratorTypeDynamic $NewmarkGamma $NewmarkBeta

# ANALYSIS -- defines what type of analysis is to be performed (http://opensees.berkeley.edu/OpenSees/m ... al/324.htm)
# Static Analysis -- solves the KU=R problem, without the mass or damping matrices.
# Transient Analysis -- solves the time-dependent analysis. The time step in this type of analysis is constant. The time step in the output is also constant.
# variableTransient Analysis -- performs the same analysis type as the Transient Analysis object. The time step, however, is variable. This method is used when
# there are convergence problems with the Transient Analysis object at a peak or when the time step is too small. The time step in the output is also variable.
variable analysisTypeDynamic Transient
analysis $analysisTypeDynamic
###########################################################################################################



file 4:-
# --------------------------------------------------------------------------------------------------
# LibUnits.tcl -- define system of units
#
#

# define UNITS ----------------------------------------------------------------------------
set mm 1.; # define basic units -- output units
set N 1.; # define basic units -- output units
set sec 1.; # define basic units -- output units
set LunitTXT "mm"; # define basic-unit text for output
set FunitTXT "N"; # define basic-unit text for output
set TunitTXT "sec"; # define basic-unit text for output
set MPa [expr $N/pow($mm,2)];
set mm2 [expr $mm*$mm]; # mm^2
set mm3 [expr $mm*$mm*$mm]; # mm^3
set cm [expr $mm*10]; # centimeter, needed for displacement input in MultipleSupport excitation
set m [expr $mm*1000];
set KN [expr $N*1000];
set PI 3.14; # define constants
set g 9810*$mm/pow($sec,2); # gravitational acceleration
#################################################################################################



file 5:-
###########################################################################
# ReadSMDFile $inFilename $outFilename $dt #
###########################################################################
# read gm input format
#
#
#
# A procedure which parses a ground motion record from the PEER
# strong motion database by finding dt in the record header, then
# echoing data values to the output file.
#
# Formal arguments
# inFilename -- file which contains PEER strong motion record
# outFilename -- file to be written in format G3 can read
# dt -- time step determined from file header
#
# Assumptions
# The header in the PEER record is, e.g., formatted as follows:
# PACIFIC ENGINEERING AND ANALYSIS STRONG-MOTION DATA
# IMPERIAL VALLEY 10/15/79 2319, EL CENTRO ARRAY 6, 230
# ACCELERATION TIME HISTORY IN UNITS OF G
# NPTS= 3930, DT= .00500 SEC

proc ReadSMDFile {inFilename outFilename dt} {
# read gm input format

# Pass dt by reference
upvar $dt DT

# Open the input file and catch the error if it can't be read
if [catch {open $inFilename r} inFileID] {
puts stderr "Cannot open $inFilename for reading"
} else {
# Open output file for writing
set outFileID [open $outFilename w]

# Flag indicating dt is found and that ground motion
# values should be read -- ASSUMES dt is on last line
# of header!!!
set flag 0

# Look at each line in the file
foreach line [split [read $inFileID] \n] {

if {[llength $line] == 0} {
# Blank line --> do nothing
continue
} elseif {$flag == 1} {
# Echo ground motion values to output file
puts $outFileID $line
} else {
# Search header lines for dt
foreach word [split $line] {
# Read in the time step
if {$flag == 1} {
set DT $word
break
}
# Find the desired token and set the flag
if {[string match $word "DT="] == 1} {set flag 1}
}
}
}

# Close the output file
close $outFileID

# Close the input file
close $inFileID
}
};
###########################################################################



file 6:-
puts " -- Uniaxial Inelastic Material, Fiber RC-Section, Nonlinear Model --"
puts " -- Multiple-support Earthquake Excitation --"
source Ex4.Portal2D.build.InelasticFiberSection.tcl
source Ex4.Portal2D.analyze.Dynamic.EQ.Uniform.tcl

############################################################################



this is modelling, is there any problem with modelling, plz anyone can help me to compare it with idarc resutl of Story drift or Drift.....................

can anyone help me to compare drift result of OPENSEES and IDARC of a portal frame, analysis type is dynamic with time history everything in modelling is done but some due to some parameters the results are different
My email id is " punits95@hotmail.com " & " coolpunitrocks@gmail.com " plz plz plz contact me here or at mail i will send u all result and modelling of both OPENSEES and IDARC.....
thank you!!

Post Reply