No Integrator Warning

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

Moderators: silvia, selimgunay, Moderators

Post Reply
yfeng
Posts: 73
Joined: Thu Aug 28, 2008 11:47 am
Location: North Carolina State University

No Integrator Warning

Post by yfeng »

Hi Vesna,
I wrote a static analysis script for rotational spring. Though I specified integrator to be displacement control, OpenSees kept showing me warning " no Integrator specified". I checked a lot of times and failed to find the problem. Can you take a look at it? Thanks

# unit is in ksi kips

# 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 DataTrial; # set up name of data directory (can remove this)
file mkdir $dataDir; # create data directory


# define GEOMETRY -------------------------------------------------------------
# define structure-geometry paramters

node 1 0 0
node 2 0 0


print node;


# BOUNDARY CONDITIONS
fix 1 1 1 1;
fix 2 1 1 0; # fix axial first



# Definition of parameters for fiber RC section

# MATERIAL parameters

set matTag 1;
set fc 50.; # MPa confined concrete strength
set E [expr 5000*pow($fc,0.5)]; # MPa young's modulous of confined concrete strength

uniaxialMaterial Elastic $matTag $E;

set ColSecTag 1;

# section GEOMETRY -------------------------------------------------------------
set DSec 609.; # Column Diameter
set Across [expr (3.14/4)*pow($DSec,2)];

set IDColTransf 1; # all columns
set ColTransfType Linear ; # options: Linear, PDelta, Corotational
# set ColTransfType PDelta ; # options: Linear, PDelta, Corotational
geomTransf $ColTransfType $IDColTransf; # only columns can have PDelta effects (gravity effects)

# Create the columns using displacement controlled beam-column elements
# tag ndI ndJ secID
element zeroLength 1 1 2 -mat $matTag -dir 6
#print element;

recorder Node -file $dataDir/disp.xls -node 2 -dof 3 disp; # check input dipl
recorder Element -file $dataDir/eleglobal_1.xls -time -ele 1 globalForce;

set IDloadTag 100; # load tag

# -------------define computational parameters for cyclic test---------
source LibAnalysisStaticParameters2.tcl;


source test9US.tcl; # displacement history for test

# create load pattern for monotonic load

set Hload 100; # define the lateral load as a proportion of the weight so that the pseudo time equals the lateral-load coefficient when using linear load pattern
set iPushNode 2; # define nodes where lateral load is applied in static lateral analysis
pattern Plain 200 Linear { ; # define load pattern -- generalized
load $iPushNode 0.0 0.0 $Hload
}

# Set up parameters that are particular to the model for displacement control

set IDctrlNode 2; # node where displacement is read for displacement control
set IDctrlDOF 3; # degree of freedom of displacement read for displacement control

set Dpre 0;
set Dincr 0.1;

foreach Dmax $iDmax {
set Ddiff [expr $Dmax-$Dpre];
set Ncycles [expr int(abs($Ddiff)/$Dincr)];
if {$Ddiff<0} { # avoid the divide by zero
set dx [expr -1*$Dincr];
} else {
set dx $Dincr;
# puts "Spring Simulation not Done";
}
integrator DisplacementControl $IDctrlNode $IDctrlDOF $dx;
analyze $Ncycles;
set Dpre $Dmax;
}; # end of iDmaxCycl
Yuhao Feng
Research Assistant & graduate student
CCEE Dept. North Carolina State University
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: No Integrator Warning

Post by vesna »

What I noticed is that you are missing ; after { at this line:

if {$Ddiff<0} { # avoid the divide by zero
Post Reply