I have encountered a very strange problem while running a batch process.
I have a 2D model (3Story SAC Building-nonlinear time-history problem) consisting nonlinearBeamColumn elements, subjected to a series of earthquake time-histories (accel.). The script is executed in a single opensees session.
The run terminates itself with the following error message:
WARNING - ForceBeamColumn2d::update - failed to get compatible element forces & deformations for element: 1104(dW: << 2.63762e-005)
The strange point comes as:
After I checked my model and coding, I have tried to execute the same "single" earthquake acceleration time-hist. several times (exactly 5 times) following one another in the same batch process, under the same opensees session:
It solved the very same problem at first 2 runs, but failed at the rest; giving the same warning message at the very begining of the time-hist. analysis!!!!
To recheck myself I have executed the batch process from the command prompt (i.e. c:\opensees run.tcl) and it solved without any problem or warning!!!
If there is problem with the model, opensees should not have solved it at 1st two runs. Why does it fail at the rest ???? This problem becomes very important when I try to run incremental dynamic analysis!!
The code is a little bit long so I could not post it on this message. In case anyone (especially Silvia) is interested I may send it to her/his email.
My email is dutkutug@gmail.com
Urgent help needed!! - Batch Process in Opensees -
Moderators: silvia, selimgunay, Moderators
Urgent help needed!! - Batch Process in Opensees -
Deniz Utkutug
METU, Turkey
METU, Turkey
can you zip up the directory and send me the file(s) so i can have a look at them.
fmckenna@ce.berkeley.edu
fmckenna@ce.berkeley.edu
the similar problem on batch analysis
I have this problem too. I use opensees batch analysis. I use foreach command. On each loop i analysis one earthquake and then change Pga and run the analysis again. but always after 3 or 4 analysis i got this message
" BandArpakSolver::Noshift could be applied during inplicit Arnoldi update,try increasing NCV
domain error:argument not in valid range"
Then after this error occurred, when I change my foreach, that opensees continue the job, exactly on the same opensees session, it do the same analysis without any problem. Again after for or five analysis i get the same message again.
I use wipe command at the end of each analysis on the foreach loop. I use fiber section with nonlinear beamcolumn element.
I use three foreach loops, one for changing earthquake and the other for earthquake intencity and the third one for changing model.
I will be grateful if someone has a way to overcome this problem.
" BandArpakSolver::Noshift could be applied during inplicit Arnoldi update,try increasing NCV
domain error:argument not in valid range"
Then after this error occurred, when I change my foreach, that opensees continue the job, exactly on the same opensees session, it do the same analysis without any problem. Again after for or five analysis i get the same message again.
I use wipe command at the end of each analysis on the foreach loop. I use fiber section with nonlinear beamcolumn element.
I use three foreach loops, one for changing earthquake and the other for earthquake intencity and the third one for changing model.
I will be grateful if someone has a way to overcome this problem.
I shortened the script. this script call model.tcl that for each analysis generate nodes, elements, masses and constraints.
###Run.tcl Start
wipe
# Define earthquake excitation
# ----------------------------
foreach pga {0.05 0.1 0.15 0.25 0.35 0.45 0.55 0.7} {
set pg [expr int($pga*100)]
set fac8 [expr $pga/0.1690589*980]
set fac9 [expr $pga/0.2277347*980]
set fac10 [expr $pga/0.2285088*980]
set fbimp165x "Path -filePath fbimp165x.txt -dt 0.01 -factor $fac8 "
set fbpalm535x "Path -filePath fbpalm535x.txt -dt 0.02 -factor $fac9 "
set fbcape811x "Path -filePath fbcape811x.txt -dt 0.02 -factor $fac10 "
foreach earthquake {8 9 10} {
if {$earthquake==8} {
set earthx $fbimp165x
set earthtime 30.0
}
puts "0"
if {$earthquake==9} {
set earthx $fbpalm535x
set earthtime 20
}
if {$earthquake==10} {
set earthx $fbcape811x
set earthtime 30
}
puts "PGA=$pga $pg"
puts "1"
foreach model {1 2 3} {
#changing model properties and reading geometry configuration from file
switch $model {
1 { set wallWidth1 126.4753
set E1 294903.2
set wallWidth2 126.4753
set E2 294903.2
set wallWidth3 126.4753
set E3 294903.2
puts " $wallWidth1 , $wallWidth2 , $$wallWidth3"
}
2 { set wallWidth1 105.8304
set E1 458137.1
set wallWidth2 126.4753
set E2 294903.2
set wallWidth3 144.219
set E3 181041.8
puts " $wallWidth1 , $wallWidth2 , $$wallWidth3"
}
3 { set wallWidth1 111.3551
set E1 435407.5
set wallWidth2 126.4753
set E2 294903.2
set wallWidth3 139.9973
set E3 186501.2
puts " $wallWidth1 , $wallWidth2 , $$wallWidth3"
}
}
puts "2"
source model.tcl
puts "3"
pattern UniformExcitation 2 2 -accel $earthx
# -----------------------
# End of model generation
# -----------------------
# ----------------------------
# Start of analysis generation
# ----------------------------
# Create the convergence test
# tol maxIter printFlag
test EnergyIncr 1.0e-4 100 3
# Create the solution algorithm
algorithm Newton
# Create the system of equation storage and solver
system SparseGeneral -piv
# Create the constraint handler
constraints Transformation
#damping 5%
# ------------PeriodFreq&Damping.tcl-----------------------------------------------------------------
# determine Natural Period, Frequency & damping parameters for SDOF
set PI 3.141592654
set xDamp 0.05; # damping ratio (0.02-0.05-typical)
set lambda [eigen 1]
set omega [expr pow($lambda,0.5)]
set Tperiod [expr 2*$PI/$omega]; # period (sec.)
puts $Tperiod
#set alphaM 0; # stiffness-prop. RAYLEIGH damping parameter; D = alphaM*M
set alphaM [expr 2*$xDamp*$omega]
set betaK 0; # stiffness proportional damping; +beatK*KCurrent
#set betaKcomm [expr 2*$xDamp/$omega]; # mass-prop. RAYLEIGH damping parameter; +betaKcomm*KlastCommitt
set betaKcomm 0.0
set betaKinit 0; # initial-stiffness proportional damping +beatKinit*Kini
# Create the time integration scheme
# gamma beta
integrator Newmark 0.5 0.25 $alphaM $betaK $betaKinit $betaKcomm
# Create the DOF numberer
numberer RCM
# Create the transient analysis
analysis Transient
# --------------------------
# End of analysis generation
# --------------------------
# ----------------------------
# Start of recorder generation
# ----------------------------
# Record DOF 1 and 2 displacements at nodes 9, 14, and 19
#recorder Node -file armin1.out -time -node 1 2 3 4 5 -dof 1 2 3 4 5 6 disp
recorder Node -file disp_$pg$model$earthquake.out -time -node 1 3 2 -dof 2 disp
recorder Node -file rot_$pg$model$earthquake.out -time -node 3 -dof 6 disp
# --------------------------
# End of recorder generation
# --------------------------
# --------------------
# Perform the analysis
# --------------------
# Analysis duration of 20 seconds
# numSteps dt
set numstep 2
set numstep [expr int($earthtime/0.001)]
analyze $numstep 0.001
puts "analysis model $model done"
wipe
}
puts "analysis $earthquake = $earthx done"
}
}
####run.tcl end
###Run.tcl Start
wipe
# Define earthquake excitation
# ----------------------------
foreach pga {0.05 0.1 0.15 0.25 0.35 0.45 0.55 0.7} {
set pg [expr int($pga*100)]
set fac8 [expr $pga/0.1690589*980]
set fac9 [expr $pga/0.2277347*980]
set fac10 [expr $pga/0.2285088*980]
set fbimp165x "Path -filePath fbimp165x.txt -dt 0.01 -factor $fac8 "
set fbpalm535x "Path -filePath fbpalm535x.txt -dt 0.02 -factor $fac9 "
set fbcape811x "Path -filePath fbcape811x.txt -dt 0.02 -factor $fac10 "
foreach earthquake {8 9 10} {
if {$earthquake==8} {
set earthx $fbimp165x
set earthtime 30.0
}
puts "0"
if {$earthquake==9} {
set earthx $fbpalm535x
set earthtime 20
}
if {$earthquake==10} {
set earthx $fbcape811x
set earthtime 30
}
puts "PGA=$pga $pg"
puts "1"
foreach model {1 2 3} {
#changing model properties and reading geometry configuration from file
switch $model {
1 { set wallWidth1 126.4753
set E1 294903.2
set wallWidth2 126.4753
set E2 294903.2
set wallWidth3 126.4753
set E3 294903.2
puts " $wallWidth1 , $wallWidth2 , $$wallWidth3"
}
2 { set wallWidth1 105.8304
set E1 458137.1
set wallWidth2 126.4753
set E2 294903.2
set wallWidth3 144.219
set E3 181041.8
puts " $wallWidth1 , $wallWidth2 , $$wallWidth3"
}
3 { set wallWidth1 111.3551
set E1 435407.5
set wallWidth2 126.4753
set E2 294903.2
set wallWidth3 139.9973
set E3 186501.2
puts " $wallWidth1 , $wallWidth2 , $$wallWidth3"
}
}
puts "2"
source model.tcl
puts "3"
pattern UniformExcitation 2 2 -accel $earthx
# -----------------------
# End of model generation
# -----------------------
# ----------------------------
# Start of analysis generation
# ----------------------------
# Create the convergence test
# tol maxIter printFlag
test EnergyIncr 1.0e-4 100 3
# Create the solution algorithm
algorithm Newton
# Create the system of equation storage and solver
system SparseGeneral -piv
# Create the constraint handler
constraints Transformation
#damping 5%
# ------------PeriodFreq&Damping.tcl-----------------------------------------------------------------
# determine Natural Period, Frequency & damping parameters for SDOF
set PI 3.141592654
set xDamp 0.05; # damping ratio (0.02-0.05-typical)
set lambda [eigen 1]
set omega [expr pow($lambda,0.5)]
set Tperiod [expr 2*$PI/$omega]; # period (sec.)
puts $Tperiod
#set alphaM 0; # stiffness-prop. RAYLEIGH damping parameter; D = alphaM*M
set alphaM [expr 2*$xDamp*$omega]
set betaK 0; # stiffness proportional damping; +beatK*KCurrent
#set betaKcomm [expr 2*$xDamp/$omega]; # mass-prop. RAYLEIGH damping parameter; +betaKcomm*KlastCommitt
set betaKcomm 0.0
set betaKinit 0; # initial-stiffness proportional damping +beatKinit*Kini
# Create the time integration scheme
# gamma beta
integrator Newmark 0.5 0.25 $alphaM $betaK $betaKinit $betaKcomm
# Create the DOF numberer
numberer RCM
# Create the transient analysis
analysis Transient
# --------------------------
# End of analysis generation
# --------------------------
# ----------------------------
# Start of recorder generation
# ----------------------------
# Record DOF 1 and 2 displacements at nodes 9, 14, and 19
#recorder Node -file armin1.out -time -node 1 2 3 4 5 -dof 1 2 3 4 5 6 disp
recorder Node -file disp_$pg$model$earthquake.out -time -node 1 3 2 -dof 2 disp
recorder Node -file rot_$pg$model$earthquake.out -time -node 3 -dof 6 disp
# --------------------------
# End of recorder generation
# --------------------------
# --------------------
# Perform the analysis
# --------------------
# Analysis duration of 20 seconds
# numSteps dt
set numstep 2
set numstep [expr int($earthtime/0.001)]
analyze $numstep 0.001
puts "analysis model $model done"
wipe
}
puts "analysis $earthquake = $earthx done"
}
}
####run.tcl end
Last edited by Armin on Wed Dec 07, 2005 1:47 am, edited 1 time in total.
Continue : problem on bach process
I think batch process will stop when it arive to eigen command
dear Sylvia
model.tcl file:
********************************
model BasicBuilder -ndm 3 -ndf 6
set h 300.0;
set by 1800.0;
set bx 2000.0;
node 1 [expr -$bx/2] 0 [expr $h]
node 2 [expr $bx/2] 0 [expr $h]
node 4 [expr -$bx/2] 0 0
node 5 [expr $bx/2] 0 0
node 6 0 0 0
node 7 0 [expr -$by/2] [expr $h]
node 8 0 [expr $by/2] [expr $h]
node 9 0 [expr -$by/2] 0
node 10 0 [expr $by/2] 0
node 3 0 0 [expr $h]
fix 4 1 1 1 1 1 1
fix 5 1 1 1 1 1 1
fix 6 1 1 1 1 1 1
fix 9 1 1 1 1 1 1
fix 10 1 1 1 1 1 1
rigidDiaphragm 3 3 1 2 7 8
fix 3 0 0 0 0 0 0
uniaxialMaterial Steel01 3 200 $E1 1e-5
uniaxialMaterial Steel01 4 200 $E2 1e-5
uniaxialMaterial Steel01 5 200 $E3 1e-5
set wallThickness 5.0
puts "wall width1= $wallWidth1"
set y1 [expr $wallWidth1/2.0]
set z1 [expr $wallThickness/2.0]
section Fiber 1 {
patch rect 3 50 2 [expr -$y1] [expr -$z1] [expr $y1] [expr $z1]
}
set wallThickness 5.0
puts "wall width2= $wallWidth2"
set y1 [expr $wallWidth2/2.0]
set z1 [expr $wallThickness/2.0]
section Fiber 2 {
patch rect 4 50 2 [expr -$y1] [expr -$z1] [expr $y1] [expr $z1]
}
set wallThickness 5.0
set y1 [expr $wallWidth3/2.0]
set z1 [expr $wallThickness/2.0]
section Fiber 3 {
patch rect 5 50 2 [expr -$y1] [expr -$z1] [expr $y1] [expr $z1]
}
set PDelta "OFF"
if {$PDelta == "ON"} {
geomTransf LinearWithPDelta 1 0 1 0
geomTransf LinearWithPDelta 2 1 0 0
} else {
geomTransf Linear 1 0 1 0
geomTransf Linear 2 1 0 0
geomTransf Linear 3 0 0 1
}
set GJ 1.0e-5;
uniaxialMaterial Elastic 10 $GJ
section Aggregator 4 10 T -section 1
section Aggregator 5 10 T -section 2
section Aggregator 6 10 T -section 3
set np 8
element nonlinearBeamColumn 1 4 1 $np 4 2
element nonlinearBeamColumn 2 5 2 $np 6 2
element nonlinearBeamColumn 3 6 3 $np 5 2
element nonlinearBeamColumn 4 9 7 $np 5 1
element nonlinearBeamColumn 5 6 3 $np 5 1
element nonlinearBeamColumn 6 10 8 $np 5 1
set g 980; # Gravitational constant
set m [expr 306000/$g]
set i [expr $m*($bx*$bx+$by*$by)/12.0]
mass 3 $m $m 0 0 0 $i
***************************
best regards
model.tcl file:
********************************
model BasicBuilder -ndm 3 -ndf 6
set h 300.0;
set by 1800.0;
set bx 2000.0;
node 1 [expr -$bx/2] 0 [expr $h]
node 2 [expr $bx/2] 0 [expr $h]
node 4 [expr -$bx/2] 0 0
node 5 [expr $bx/2] 0 0
node 6 0 0 0
node 7 0 [expr -$by/2] [expr $h]
node 8 0 [expr $by/2] [expr $h]
node 9 0 [expr -$by/2] 0
node 10 0 [expr $by/2] 0
node 3 0 0 [expr $h]
fix 4 1 1 1 1 1 1
fix 5 1 1 1 1 1 1
fix 6 1 1 1 1 1 1
fix 9 1 1 1 1 1 1
fix 10 1 1 1 1 1 1
rigidDiaphragm 3 3 1 2 7 8
fix 3 0 0 0 0 0 0
uniaxialMaterial Steel01 3 200 $E1 1e-5
uniaxialMaterial Steel01 4 200 $E2 1e-5
uniaxialMaterial Steel01 5 200 $E3 1e-5
set wallThickness 5.0
puts "wall width1= $wallWidth1"
set y1 [expr $wallWidth1/2.0]
set z1 [expr $wallThickness/2.0]
section Fiber 1 {
patch rect 3 50 2 [expr -$y1] [expr -$z1] [expr $y1] [expr $z1]
}
set wallThickness 5.0
puts "wall width2= $wallWidth2"
set y1 [expr $wallWidth2/2.0]
set z1 [expr $wallThickness/2.0]
section Fiber 2 {
patch rect 4 50 2 [expr -$y1] [expr -$z1] [expr $y1] [expr $z1]
}
set wallThickness 5.0
set y1 [expr $wallWidth3/2.0]
set z1 [expr $wallThickness/2.0]
section Fiber 3 {
patch rect 5 50 2 [expr -$y1] [expr -$z1] [expr $y1] [expr $z1]
}
set PDelta "OFF"
if {$PDelta == "ON"} {
geomTransf LinearWithPDelta 1 0 1 0
geomTransf LinearWithPDelta 2 1 0 0
} else {
geomTransf Linear 1 0 1 0
geomTransf Linear 2 1 0 0
geomTransf Linear 3 0 0 1
}
set GJ 1.0e-5;
uniaxialMaterial Elastic 10 $GJ
section Aggregator 4 10 T -section 1
section Aggregator 5 10 T -section 2
section Aggregator 6 10 T -section 3
set np 8
element nonlinearBeamColumn 1 4 1 $np 4 2
element nonlinearBeamColumn 2 5 2 $np 6 2
element nonlinearBeamColumn 3 6 3 $np 5 2
element nonlinearBeamColumn 4 9 7 $np 5 1
element nonlinearBeamColumn 5 6 3 $np 5 1
element nonlinearBeamColumn 6 10 8 $np 5 1
set g 980; # Gravitational constant
set m [expr 306000/$g]
set i [expr $m*($bx*$bx+$by*$by)/12.0]
mass 3 $m $m 0 0 0 $i
***************************
best regards
Stop analysis during Batch process
I have a similar problem too. Using foreach command, I keep running the analyses by changing earthquake records. Between each analysis, I use wipe command and recreate the model for the next earthquake record.
However, Opensees stops always after 3 earthquake records and I got a small windows with error message in the middle of 4th analysis :
"opensees.exe has encountered a problem and needs to close..."
This seems to be nothing to do with modeling since, when I rerun it by changing foreach, opensees can finish the job with the earthquake record where it stopped previously. However, again after 3 analyses, the same message windows pops up during 4th analysis.
I use v1.7.3 on Windows plaform. I found similar postings but I have not found how to fix it. Does anyone know why? Besides I cannot post my scripts here because it is so long.
However, Opensees stops always after 3 earthquake records and I got a small windows with error message in the middle of 4th analysis :
"opensees.exe has encountered a problem and needs to close..."
This seems to be nothing to do with modeling since, when I rerun it by changing foreach, opensees can finish the job with the earthquake record where it stopped previously. However, again after 3 analyses, the same message windows pops up during 4th analysis.
I use v1.7.3 on Windows plaform. I found similar postings but I have not found how to fix it. Does anyone know why? Besides I cannot post my scripts here because it is so long.
still same problem
thanks Silvia,
But I still have same problem even if I run it from DOS prompt.
But I still have same problem even if I run it from DOS prompt.
I have the same problem while running sampling analysis. After certain number of runs (~50 for a fairly large model, ~10000 for a small model) I get the same message
"opensees.exe has encountered a problem and needs to close..."
I would like to know why this occurs. The problem is independent of the seed value and convergence issues.
I ran it from DOS prompt.
Thanks.
"opensees.exe has encountered a problem and needs to close..."
I would like to know why this occurs. The problem is independent of the seed value and convergence issues.
I ran it from DOS prompt.
Thanks.