Analysis error

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

Moderators: silvia, selimgunay, Moderators

orodzarrin
Posts: 16
Joined: Sun Mar 17, 2013 3:26 am
Location: Eastern Mediterranean University

Analysis error

Post by orodzarrin » Fri Oct 25, 2013 9:56 am

these is problem about this commend. i can not find out how to correct it.
pleas help me.

integrator DisplacementControl 5 1 0.1

this is its error :
can't read "nodeTag": no such variable

the number 5 at the first parameter is a joint.

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

Re: Analysis error

Post by fmk » Fri Oct 25, 2013 10:00 am

node 5 needs to be a node in the model. does the command "print -node 5" work?

orodzarrin
Posts: 16
Joined: Sun Mar 17, 2013 3:26 am
Location: Eastern Mediterranean University

Re: Analysis error

Post by orodzarrin » Fri Oct 25, 2013 10:13 am

it is in the model.
no it doesn't work.

brag006
Posts: 173
Joined: Wed Feb 15, 2012 1:26 pm
Location: University of Auckland

Re: Analysis error

Post by brag006 » Fri Oct 25, 2013 12:37 pm

Can you post your node definitions

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

Re: Analysis error

Post by fmk » Fri Oct 25, 2013 2:52 pm

if that command does not work .. you have no node 5 and thus cannot control the analysis using the displacement at node 5.

orodzarrin
Posts: 16
Joined: Sun Mar 17, 2013 3:26 am
Location: Eastern Mediterranean University

Re: Analysis error

Post by orodzarrin » Fri Oct 25, 2013 11:27 pm

######################## Nodes #######################
node 1 0.0 0.0
node 2 5.0 0.0
node 3 9.0 0.0
node 4 0.0 3.0
node 5 5.0 3.0
node 6 9.0 3.0
node 7 0.0 6.0
node 8 5.0 6.0
node 9 9.0 6.0
node 10 0.0 9.0
node 11 5.0 9.0
node 12 9.0 9.0

I used it for other node such as 12 but it doesn't work as well.

brag006
Posts: 173
Joined: Wed Feb 15, 2012 1:26 pm
Location: University of Auckland

Re: Analysis error

Post by brag006 » Sun Oct 27, 2013 6:45 pm

Can you post all your code. Maybe I can help you with the error

orodzarrin
Posts: 16
Joined: Sun Mar 17, 2013 3:26 am
Location: Eastern Mediterranean University

Re: Analysis error

Post by orodzarrin » Sun Oct 27, 2013 10:04 pm

######################### Model Class ########################
wipe

set GMdir "GMfiles"; # ground-motion file directory
source LibUnits.tcl; # define basic and system units




######################### Dimansion ######################
model basic -ndm 2 -ndf 3

puts "ok 1"
######################## Nodes #######################
node 1 0.0 0.0; # First Fix Support #
node 2 5.0 0.0 # second Fix Support #
node 3 9.0 0.0 #-mass 4.0 0.0 0.0
node 4 0.0 3.0 #-mass 4.0 0.0 0.0
node 5 5.0 3.0 #-mass 4.0 0.0 0.0
node 6 9.0 3.0 #-mass 4.0 0.0 0.0
node 7 0.0 6.0 #-mass 0.0 0.0 0.0
node 8 5.0 6.0 #-mass 0.0 0.0 0.0
node 9 9.0 6.0 #-mass 0.0 0.0 0.0
node 10 0.0 9.0 #-mass 0.0 0.0 0.0
node 11 5.0 9.0 #-mass 0.0 0.0 0.0
node 12 9.0 9.0 #-mass 0.0 0.0 0.0

puts "ok 2"
########################## Support ########################
fix 1 1 1 1
fix 2 1 1 1
fix 3 1 1 1

puts "ok 3"
########################## EqualDof #######################1 2 mafsal mikonad#
equalDOF 4 3 1 2
equalDOF 5 14 1 2
equalDOF 5 15 1 2
equalDOF 6 16 1 2
equalDOF 7 17 1 2
equalDOF 8 18 1 2
equalDOF 8 19 1 2
equalDOF 9 20 1 2
equalDOF 10 21 1 2
equalDOF 11 22 1 2
equalDOF 11 23 1 2
equalDOF 12 24 1 2

puts "ok 4"
########################## rigid link beam ################
rigidLink beam 4 5 # connect node 3 to node 4 via a rigid link-beam #
rigidLink beam 4 6
rigidLink beam 7 8
rigidLink beam 7 9
rigidLink beam 10 11
rigidLink beam 10 12

puts "ok 5"
########################## Material ##fy=235359.605##E=200100000##b=0.02##
uniaxialMaterial Steel01 1 24 2 0.02

puts "ok 6"
########################## section #########################aval (y) bad (z) #
section Fiber 1 {
patch quad 1 5 4 0.05 -0.05 0.05 0.05 0.06 0.05 0.06 -0.05
patch quad 1 3 10 -0.05 -0.005 -0.05 0.005 0.05 0.005 0.05 -0.005
patch quad 1 5 4 -0.06 -0.05 -0.06 0.05 -0.05 0.05 -0.05 -0.05
}

puts "ok 7"
#

################################## Column 1 ## Beam 2 #
geomTransf PDelta 1
geomTransf Linear 2

puts "ok 8"

puts "ok 9"
############################################## Nonlinear Beam - column Element ######

element nonlinearBeamColumn 1 1 4 5 1 1
element nonlinearBeamColumn 2 2 5 5 1 1
element nonlinearBeamColumn 3 3 6 5 1 1
element nonlinearBeamColumn 4 4 5 5 1 2
element nonlinearBeamColumn 5 5 6 5 1 2
element nonlinearBeamColumn 6 4 7 5 1 1
element nonlinearBeamColumn 7 5 8 5 1 1
element nonlinearBeamColumn 8 6 9 5 1 1
element nonlinearBeamColumn 9 7 8 5 1 2
element nonlinearBeamColumn 10 8 9 5 1 2
element nonlinearBeamColumn 11 7 10 5 1 1
element nonlinearBeamColumn 12 8 11 5 1 1
element nonlinearBeamColumn 13 9 12 5 1 1
element nonlinearBeamColumn 14 10 11 5 1 2
element nonlinearBeamColumn 15 11 12 5 1 2

puts "ok 10"
######################## Path time series ##### baray moshakhas kardan bare zelzele#######292#


pattern Plain 1 Linear { ; # define LoadPattern 1. impose load in a linear manner
load 5 100.0 0.0 20.0; # apply force and moment at node 5
}

puts "ok 11"
################################################ Record #############303#
recorder Node -file dataDir/NodeDis5.out -time -node 5 -dof 1 disp
recorder Node -file dataDir/Node6.eig -time -node 6 -dof 1 "eigen 1"
recorder Node -file dataDir/RNode.out -time -node 1 2 -dof 1 reaction
recorder Drift -file dataDir/drift.out -time -inode 3 4 -jnode 5 6 -dof 1 perpDrin 2
recorder Element -file dataDir/ele2local.out -time -ele 2 localForce
recorder Element -file dataDir/ele2local.out -time -ele 2 globalForce
recorder Element -file ele1trusslocal.out -time -ele 1 localForce
recorder display push 10 10 500 500 #? it has got problem #
prp 4.0 6.0
vup 0 1 0
display 1 5 200
recorder plot node6.out push 525 10 400 400 -column 2 1;


puts "ok 12"
##################################################### analysis ##########318# pushover

constraints Transformation
numberer RCM
system BandGeneral
test NormDispIncr 1.0e-6 100 3
algorithm Newton
integrator DisplacementControl 3 1 0.1 ### in baray push hast#integrator DisplacementControl $nodeTag $dofTag $dU1(meghdar taghir makani ke bayad bedahim 20cm ast- 200*0.0001=20cm) <$Jd $minDu $maxDu>
analysis Static
analyze 100

puts "ok 13"
################################################

pinash
Posts: 5
Joined: Sun Oct 27, 2013 10:08 pm
Location: tehran

Re: Analysis error

Post by pinash » Sun Oct 27, 2013 10:57 pm

may be
do not use rigid Diaphragm in 2d model
it works in 3d model
use equalldof for Rigid Diaphragm
: equalDOF 10 20 1:
why you use ridid link?
Last edited by pinash on Mon Oct 28, 2013 10:51 am, edited 1 time in total.

orodzarrin
Posts: 16
Joined: Sun Mar 17, 2013 3:26 am
Location: Eastern Mediterranean University

Re: Analysis error

Post by orodzarrin » Mon Oct 28, 2013 1:45 am

I modeled 2D. I tried to use Ridig Diaphram and equalDOF 10 20 1
but it doesn't work.

orodzarrin
Posts: 16
Joined: Sun Mar 17, 2013 3:26 am
Location: Eastern Mediterranean University

Re: Analysis error

Post by orodzarrin » Mon Oct 28, 2013 7:28 am

could you run all of the above codes by yourself and check the problem please?
I need ur help.

pinash
Posts: 5
Joined: Sun Oct 27, 2013 10:08 pm
Location: tehran

Re: Analysis error

Post by pinash » Mon Oct 28, 2013 10:33 am

first
add ; in your Description text
reform

rigidLink beam 4 5; # connect node 3 to node 4 via a rigid link-beam #

use ; # befor your Description text

if didnt correct then:
delete text:
puts "ok 4"
########################## rigid link beam ################
rigidLink beam 4 5 # connect node 3 to node 4 via a rigid link-beam #
rigidLink beam 4 6
rigidLink beam 7 8
rigidLink beam 7 9
rigidLink beam 10 11
rigidLink beam 10 12

then reform
########################## rigid Diaphragm ################
equalDOF 4 6 1; # Rigid Diaphragm; Forces the horizontal displacement of column nodes to be equal
equalDOF 7 9 1
equalDOF 10 12 1

brag006
Posts: 173
Joined: Wed Feb 15, 2012 1:26 pm
Location: University of Auckland

Re: Analysis error

Post by brag006 » Mon Oct 28, 2013 12:01 pm

There are many mistakes in your code. You need to make the following changes:
1. Commenting: In many places your commenting is being treated as code by the interperter. Please read on how to comment in tcl. Also use notepad++ it will help you with writing your code.
2. EqualDof: in this command you are calling many nodes that you haven't even created. I don't quite understand what you are doing there. I assume you are trying to constraint all the floor nodes? if so here is the code you should use:

equalDOF 4 5 1 2
equalDOF 4 6 1 2
equalDOF 7 8 1 2
equalDOF 7 9 1 2
equalDOF 10 11 1 2
equalDOF 10 12 1 2

3. If you are using equalDOF why are you also using rigidLink? Use one or the other not both
4. the first line is usually: model BasicBuilder -ndm 2 -ndf 3
5. Recorders. You need to review how to define recorders a few of them are incorrect. I have commented out the ones causing problems

So use the code below and it should run but I am not sure about the quality of the results:

---------------------------------------------------------------------------------------------------------------------

# ######################## Model Class ########################
wipe

set GMdir "GMfiles"; # ground-motion file directory
source LibUnits.tcl; # define basic and system units




# ######################## Dimansion ######################
# model basic -ndm 2 -ndf 3
model BasicBuilder -ndm 2 -ndf 3;

puts "ok 1"
# ####################### Nodes #######################
node 1 0.0 0.0
node 2 5.0 0.0
node 3 9.0 0.0
node 4 0.0 3.0
node 5 5.0 3.0
node 6 9.0 3.0
node 7 0.0 6.0
node 8 5.0 6.0
node 9 9.0 6.0
node 10 0.0 9.0
node 11 5.0 9.0
node 12 9.0 9.0

puts "ok 2"
# ######################### Support ########################
fix 1 1 1 1
fix 2 1 1 1
fix 3 1 1 1

puts "ok 3"
# ######################### EqualDof #######################1 2 mafsal mikonad#
equalDOF 4 5 1 2
equalDOF 4 6 1 2
equalDOF 7 8 1 2
equalDOF 7 9 1 2
equalDOF 10 11 1 2
equalDOF 10 12 1 2


puts "ok 5"
# ######################### Material ##fy=235359.605##E=200100000##b=0.02##
uniaxialMaterial Steel01 1 24 2 0.02

puts "ok 6"
# ######################### section #########################aval (y) bad (z) #
section Fiber 1 {
patch quad 1 5 4 0.05 -0.05 0.05 0.05 0.06 0.05 0.06 -0.05
patch quad 1 3 10 -0.05 -0.005 -0.05 0.005 0.05 0.005 0.05 -0.005
patch quad 1 5 4 -0.06 -0.05 -0.06 0.05 -0.05 0.05 -0.05 -0.05
}

puts "ok 7"
#

# ################################# Column 1 ## Beam 2 #
geomTransf PDelta 1
geomTransf Linear 2

puts "ok 8"

puts "ok 9"
# ############################################# Nonlinear Beam - column Element ######

element nonlinearBeamColumn 1 1 4 5 1 1
element nonlinearBeamColumn 2 2 5 5 1 1
element nonlinearBeamColumn 3 3 6 5 1 1
element nonlinearBeamColumn 4 4 5 5 1 2
element nonlinearBeamColumn 5 5 6 5 1 2
element nonlinearBeamColumn 6 4 7 5 1 1
element nonlinearBeamColumn 7 5 8 5 1 1
element nonlinearBeamColumn 8 6 9 5 1 1
element nonlinearBeamColumn 9 7 8 5 1 2
element nonlinearBeamColumn 10 8 9 5 1 2
element nonlinearBeamColumn 11 7 10 5 1 1
element nonlinearBeamColumn 12 8 11 5 1 1
element nonlinearBeamColumn 13 9 12 5 1 1
element nonlinearBeamColumn 14 10 11 5 1 2
element nonlinearBeamColumn 15 11 12 5 1 2

puts "ok 10"
# ####################### Path time series ##### baray moshakhas kardan bare zelzele#######292#


pattern Plain 1 Linear { ; # define LoadPattern 1. impose load in a linear manner
load 5 100.0 0.0 20.0; # apply force and moment at node 5
}

puts "ok 11"
# ############################################### Record #############303#
recorder Node -file NodeDis5.out -time -node 5 -dof 1 disp
# recorder Node -file Node6.out -time -node 6 -dof 1 "eigen 1"
recorder Node -file RNode.out -time -node 1 2 -dof 1 reaction
# recorder Drift -file drift.out -time -inode 3 4 -jnode 5 6 -dof 1 perpDrin 2
recorder Element -file ele2local.out -time -ele 2 localForce
recorder Element -file ele2local.out -time -ele 2 globalForce
recorder Element -file ele1trusslocal.out -time -ele 1 localForce


puts "ok 12"
# #################################################### analysis ##########318# pushover

constraints Transformation
numberer RCM
system BandGeneral
test NormDispIncr 1.0e-6 100 3
algorithm Newton
integrator DisplacementControl 5 1 0.1
analysis Static
analyze 100

puts "ok 13"
wipe

----------------------------------------------------------------------------------------------------------------------------------

orodzarrin
Posts: 16
Joined: Sun Mar 17, 2013 3:26 am
Location: Eastern Mediterranean University

Re: Analysis error

Post by orodzarrin » Wed Oct 30, 2013 9:46 pm

thank you very much.it works

shariat
Posts: 2
Joined: Mon Apr 29, 2013 7:44 am
Location: qiau

Re: Analysis error

Post by shariat » Fri Nov 01, 2013 7:33 am

I have a problem about analyze my model. i can not find out how to correct it.
pleas help me.
WARNING SuperLU::solve<void>- Error 1 returned in factorization dgstrf<>
WRANING NewtonRaphson::solve CurrentStep<> -the LinearSysOfEqn failed in solve<>
StaticAnalysis::analyze<> -the Algorithm failed at iteration: 0 with domain at
load factor 1
OpenSees > analyze failed, returned: -3 flag

Post Reply