The OpenSees Community

I am also working with Bouc-Wen model. Can you please share the parameters you used?

I have increased the tolerance to a value equal to 1e-5 and increased iteration number to 100 and I got no warning.
I hope the problem have been solved.

Dear all,
I have run the following script for the problem of base isolation but got these warnings.
However I end up get the results for the node displacement and I am not sure if they are correct.

1.What should I edit so that I get no warning?
2. Under these warnings will the results be correct or not ? Why?

I will appreciate for any available help.

WARNING: CTestNormUnbalance::test() - failed to converge
after: 6 iterations current Norm: 7.27596e-012 (max: 1e-012, Norm deltaX: 2.85079e-022)
NewtnRaphson::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.16
OpenSees > analyze failed, returned: -3 error flag
OpenSees > puts "record: $record period: $Tn damping ratio: $dampRatio max disp: $maxD"
record: Elcentro period: 1 damping ratio: 0.02 max disp: 0.006052093951979983
OpenSees > wipe all;


# Code:2DOF system for base isolation structure

OpenSees > wipe all;
OpenSees > set Tn 1
1
OpenSees > set M 980000
980000
OpenSees > set dampRatio 0.02
0.02
OpenSees > set record Elcentro
Elcentro
OpenSees > set alpha1 0.02
0.02
OpenSees > set alpha2 0.008
0.008
OpenSees > set ko1 967221231
967221231
OpenSees > set ko2 9800000000
9800000000
OpenSees > set n 8
8
OpenSees > set gamma 0.5
0.5
OpenSees > set beta 0.5
0.5
OpenSees > set Ao 1
1
OpenSees > set deltaA 0
0
OpenSees > set deltaNu 0
0
OpenSees > set deltaEta 0
0
OpenSees > #number of modes
OpenSees > set numModes 2
2
OpenSees > #set some constants
OpenSees > set g 9.81
9.81
OpenSees > set PI [expr 2.0 * asin(1.0)]
3.141592653589793
OpenSees > #derived quantaties
OpenSees > set Wn [expr 2.0 * $PI / $Tn]
6.283185307179586
OpenSees > set c [expr 2.0*$M*$Wn*$dampRatio]
246300.86404143978
OpenSees > set cb 171918.0031
171918.0031
OpenSees > # create the model
OpenSees > wipe
OpenSees > model basic -ndm 1 -ndf 1
OpenSees > node 1 0.0
OpenSees > node 2 0.0 -mass $M
OpenSees > node 3 1.0 -mass $M
OpenSees > fix 1 1
0
OpenSees > uniaxialMaterial Elastic 1 0.0 $c
OpenSees > uniaxialMaterial Elastic 2 0.0 $cb
OpenSees > uniaxialMaterial BoucWen 3 $alpha1 $ko1 $n $gamma $beta $Ao $deltaA $deltaNu $deltaEta
OpenSees > uniaxialMaterial BoucWen 4 $alpha2 $ko2 $n $gamma $beta $Ao $deltaA $deltaNu $deltaEta
OpenSees > uniaxialMaterial Parallel 5 1 3
OpenSees > uniaxialMaterial Parallel 6 2 4
OpenSees > #element truss eleno -i -j area matt
OpenSees > element zeroLength 1 1 2 -mat 6 -dir 1
OpenSees > element truss 2 2 3 1.0 5
OpenSees > # perform eigen analysis
OpenSees > set lambda [eigen -fullGenLapack $numModes];
890.48215477182214000000 11083.43872482001800000000
OpenSees > # calculate frequencies and periods of the structure
OpenSees > set omega {}
OpenSees > set f {}
OpenSees > set T {}
OpenSees > set pi 3.141593
3.141593
OpenSees > foreach lam $lambda {
lappend omega [expr sqrt($lam)]
lappend f [expr sqrt($lam)/(2*$pi)]
lappend T [expr (2*$pi)/sqrt($lam)]
}
OpenSees > puts "omega are $omega"
omega are 29.840947618529512 105.27791185628645
OpenSees > puts "periods are $T"
periods are 0.21055584696306703 0.05968190182739469
OpenSees > print node 3;

Node: 3
Coordinates : 1
Disps: 0
Velocities : 0
commitAccels: 0
Mass :
980000

Rayleigh Factor: alphaM: 0
Rayleigh Forces: 0
Eigenvectors:
1 -0.0977529

ID : 0

OpenSees > print node 2;

Node: 2
Coordinates : 0
Disps: 0
Velocities : 0
commitAccels: 0
Mass :
980000

Rayleigh Factor: alphaM: 0
Rayleigh Forces: 0
Eigenvectors:
0.0977529 1

ID : 1

OpenSees > set dT 0.02;
0.02
OpenSees > set nPts 1500;
1500
OpenSees > set factor $g
9.81
OpenSees > timeSeries Path 1 -filePath Elcentro.dat -dt $dT -factor $g
OpenSees > pattern UniformExcitation 1 1 -accel 1
OpenSees > # create the analysis
OpenSees > constraints Plain
OpenSees > integrator Newmark 0.5 [expr 1.0/6.0]
OpenSees > system ProfileSPD
OpenSees > test NormUnbalance 1.0e-12 6 0
OpenSees > algorithm Newton
OpenSees > numberer RCM
OpenSees > analysis Transient
OpenSees > rayleigh 0.7334 0.0026 0. 0.
OpenSees > set t 0.0
0.0
OpenSees > set maxT [expr (1+$nPts)*$dT];
30.02
OpenSees > set ok 0.0
0.0
OpenSees > set maxD 0.0
0.0
OpenSees > recorder Node -file disp3.out -time -node 3 -dof 1 disp
0
OpenSees > recorder Node -file disp2.out -time -node 2 -dof 1 disp
1
OpenSees > recorder Node -file disp1.out -time -node 1 -dof 1 disp
2
OpenSees > recorder Node -file vel3.out -time -node 3 -dof 1 vel
3
OpenSees > recorder Node -file accel3.out -time -node 3 -dof 1 accel
4
OpenSees > recorder Drift -file drift4.out -time -iNode 2 -jNode 3 -dof 1 -perpDirn 1
5
OpenSees > while {$ok == 0 && $t < $maxT} {
set ok [analyze 1 $dT]
if {$ok != 0} {
test NormDispIncr 1.0e-12 100 0
algorithm ModifiedNewton -initial
set ok [analyze 1 0.001]
test NormDispIncr 1.0e-12 10
algorithm Newton
}
set time [getTime]
set d [nodeDisp 3 1]
if {$d > $maxD} {
set maxD $d
} elseif {$d < [expr -$maxD]} {
set maxD [expr -$d]
}
set t [expr $t + $dT]
}
WARNING: CTestNormUnbalance::test() - failed to converge
after: 6 iterations current Norm: 7.27596e-012 (max: 1e-012, Norm deltaX: 2.85079e-022)
NewtnRaphson::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.16
OpenSees > analyze failed, returned: -3 error flag
OpenSees > puts "record: $record period: $Tn damping ratio: $dampRatio max disp: $maxD"
record: Elcentro period: 1 damping ratio: 0.02 max disp: 0.006052093951979983
OpenSees > wipe all;

Thank you very much for your kind response, let me work on it.
What about the zero length between node 2 and 4,where i want to include the impact material?
Can you please comment something!
Thank you.

Hi everyone,
I have tried to run the following script for the problem of 2DOF system considering impact between the base and the moat wall under earthquake loading.
The program(opensees) stopped running,i don't know why.
Can someone help me to check where is the problem and suggest what do.
Thank you.

#start of the code
wipe all;
set Tn 1.0
set M 100000
set dampRatio 0.02
set alpha1 0.02
set alpha2 0.008
set ko1 966240.8
set ko2 9800000
set n 8
set gamma 0.5
set beta 0.5
set Ao 1
set deltaA 0
set deltaNu 0
set deltaEta 0
set K1 1600000
set K2 326000
set delta -0.00254
set gap -0.7
#number of modes
set numModes 2
#set some constants
set g 9.810
set PI [expr 2.0 * asin(1.0)]
#derived quantaties
set Wn [expr 2.0 * $PI / $Tn]
set c [expr 2.0*$M*$Wn*$dampRatio]
# create the model
wipe
model basic -ndm 1 -ndf 1
node 1 0.0
node 2 0.0 -mass $M
node 3 1.0 -mass $M
node 4 0.0
fix 1 1
uniaxialMaterial Elastic 1 0.0 $c
uniaxialMaterial BoucWen 2 $alpha1 $ko1 $n $gamma $beta $Ao $deltaA $deltaNu $deltaEta
uniaxialMaterial BoucWen 3 $alpha2 $ko2 $n $gamma $beta $Ao $deltaA $deltaNu $deltaEta
uniaxialMaterial ImpactMaterial 4 $K1 $K2 $delta $gap
uniaxialMaterial Parallel 5 1 3
uniaxialMaterial Parallel 6 1 2
uniaxialMaterial Parallel 7 1 4
#element truss eleno -i -j area matt
element truss 1 1 2 1.0 5
element truss 2 2 3 1.0 6
element truss 3 2 4 1.0 7
# perform eigen analysis
set lambda [eigen -fullGenLapack $numModes];
# calculate frequencies and periods of the structure
set omega {}
set f {}
set T {}
set pi 3.141593
foreach lam $lambda {
lappend omega [expr sqrt($lam)]
lappend f [expr sqrt($lam)/(2*$pi)]
lappend T [expr (2*$pi)/sqrt($lam)]
}
puts "omega are $omega"
puts "periods are $T"
print node 3;
print node 2;
set dT 0.02;
set nPts 1500;
set factor $g
timeSeries Path 1 -filePath Elcentro.dat -dt $dT -factor $g
pattern UniformExcitation 1 1 -accel 1
# create the analysis
constraints Plain
integrator Newmark 0.5 [expr 1.0/6.0]
system ProfileSPD
test NormUnbalance 1.0e-12 6 0
algorithm Newton
numberer RCM
analysis Transient
rayleigh 0.7334 0.0026 0. 0.
set t 0.0
set maxT [expr (1+$nPts)*$dT];
set ok 0.0
set maxD 0.0
recorder Node -file disp3.out -time -node 3 -dof 1 disp
recorder Node -file vel3.out -time -node 3 -dof 1 vel
recorder Node -file accel3.out -time -node 3 -dof 1 accel
recorder Drift -file drift4.out -time -iNode 2 -jNode 3 -dof 1 -perpDirn 1
while {$ok == 0 && $t < $maxT} {
set ok [analyze 1 $dT]
if {$ok != 0} {
test NormDispIncr 1.0e-12 100 0
algorithm ModifiedNewton -initial
set ok [analyze 1 .01]
test NormDispIncr 1.0e-12 10
algorithm Newton
}
set time [getTime]
set d [nodeDisp 3 1]
if {$d > $maxD} {
set maxD $d
} elseif {$d < [expr -$maxD]} {
set maxD [expr -$d]
}
set t [expr $t + $dT]
}
puts "record: $record period: $Tn damping ratio: $dampRatio max disp: $maxD"
wipe all;

@Prafullamalla I sent you the files for opensees navigator hope you received them.

Ok, I will try to send you set up files.

I was able to plot the graphs for responses and the results seem to be OK,thank you once again.Have you tried to do this example using opensees navigator?

Well explained,i understand now.Thank you.I will plot graph.

Yes i checked the output files with results(did not plot the graph),does it means that convergence error should be ignored and there is no effect to the results? How do you plot graph?

No I tried it first on opensees navigator, but when I tried it in opensees it had a convergence error.

No it is not convergence error.I can send you an image for the error message ,and please if you don't mind can you text me your email to my email address nyangip@yahoo.com, I would like to discuss in detail and I may have some question on your model that I still need some explanation, I am in Wuhan.

I have tried to model this 3DOF system using opensees navigator and I get error message.can you help me please to figure out the problem

# Coded by Prafulla Malla, Sichuan University
# Structural Response 3DOF
# Report Vibration Control
# X--*************--X --************--X----************--X
# 1Fix（m,k,c) 2 m,k,c) 3 （m,k,c) 4 Top Floor
wipe all;
set Tn 1.0
set K 2.0e8
set M 4.0e5
set dampRatio 0.02
#number of modes
set numModes 3

#set some constants
set g 9.810
set PI [expr 2.0 * asin(1.0)]

#derived quantaties
set Wn [expr 2.0 * $PI / $Tn]
set c [expr 2.0*$M*$Wn*$dampRatio]

# create the model
wipe
model basic -ndm 1 -ndf 1

node 1 0.0
node 2 1.0 -mass $M
node 3 2.0 -mass $M
node 4 3.0 -mass $M

fix 1 1

uniaxialMaterial Elastic 1 $K 0.0
uniaxialMaterial Elastic 2 0.0 $c
uniaxialMaterial Parallel 3 1 2

#element truss eleno -i -j area matt
element truss 1 1 2 1.0 3
element truss 2 2 3 1.0 3
element truss 3 3 4 1.0 3


# perform eigen analysis
#-----------------------------
set lambda [eigen -fullGenLapack $numModes];

# calculate frequencies and periods of the structure
#---------------------------------------------------
set omega {}
set f {}
set T {}
set pi 3.141593

foreach lam $lambda {
lappend omega [expr sqrt($lam)]
lappend f [expr sqrt($lam)/(2*$pi)]
lappend T [expr (2*$pi)/sqrt($lam)]
}
puts "omega are $omega"
puts "periods are $T"
print node 4;
print node 3;
print node 2;

set dT 0.02;
set nPts 1500;
set factor $g



timeSeries Path 1 -filePath Elcentro.dat -dt $dT -factor $g
pattern UniformExcitation 1 1 -accel 1

# create the analysis
constraints Plain
integrator Newmark 0.5 [expr 1.0/6.0]
system ProfileSPD
test NormUnbalance 1.0e-12 6 0
algorithm Newton
numberer RCM
analysis Transient
rayleigh 0.7334 0.0026 0. 0.

set t 0.0
set maxT [expr (1+$nPts)*$dT];
set ok 0.0
set maxD 0.0
recorder Node -file disp4.out -time -node 4 -dof 1 disp
recorder Node -file vel4.out -time -node 4 -dof 1 vel
recorder Node -file accel4.out -time -node 4 -dof 1 accel
recorder Drift -file drift4.out -time -iNode 3 -jNode 4 -dof 1 -perpDirn 1
recorder Node -file mode.out -nodeRange 2 4 -dof 1 eigen ;

while {$ok == 0 && $t < $maxT} {
set ok [analyze 1 $dT]
if {$ok != 0} {
test NormDispIncr 1.0e-12 100 0
algorithm ModifiedNewton -initial
set ok [analyze 1 .01]
test NormDispIncr 1.0e-12 10
algorithm Newton
}

set time [getTime]
set d [nodeDisp 4 1]

if {$d > $maxD} {
set maxD $d
} elseif {$d < [expr -$maxD]} {
set maxD [expr -$d]
}
set t [expr $t + $dT]
}

puts "record: $record period: $Tn damping ratio: $dampRatio max disp: $maxD"
wipe all;
Prafullamalla

Posts: 34
Joined: Mon Feb 02, 2015 7:32 pm

selimgunay wrote:
> I will send you an email.

Still need your help,i could not be able to figure out this.
I will appreciate.

Thank you for sharing such a good example.Will you kindly share a typical example for base isolated 2DOF system subjected to ground motion(earthquake),i was trying to do the same and blocked somewhere.Will appreciate for your help.