problem with convergence

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

Moderators: silvia, selimgunay, Moderators

Post Reply
sisa
Posts: 27
Joined: Sun Mar 13, 2011 8:03 am
Location: Bologna

problem with convergence

Post by sisa »

I'm analyzing an infilled frame of 5 storeys and 3 bays. I'm doing a static puhs-over analysis with force proportional to mass. I've problem with the convergence because I've tried to change a lot of parameters (in analysis and in the model) but it doesnt'work. Do you have any suggest? please help me because I have to graduate 3 october and I need to this analysis!!
this is the original model...i've tried to reduce number of fiber in the columns or number of secion in the nonlinearBeamColumn but it doesen't converge again. I'm not very practice with parameter of analysis so perhaps it's possible to change something to reach covergence. I hope you can help me!!
thanks
Annalisa

-------model--------------------
#telaio nudo di 5 piani a 3 campate
wipe
set cartel TelaioTamponato
file mkdir $cartel
model BasicBuilder -ndm 2 -ndf 3

#GEOMETRIA
set Hpil 3300
set Ltr 5000
set Npiani 5

node 1 0 0
node 2 0 $Hpil
node 3 0 [expr 2*$Hpil]
node 4 0 [expr 3*$Hpil]
node 5 0 [expr 4*$Hpil]
node 6 0 [expr 5*$Hpil]
node 7 $Ltr 0
node 8 $Ltr $Hpil
node 9 $Ltr [expr 2*$Hpil]
node 10 $Ltr [expr 3*$Hpil]
node 11 $Ltr [expr 4*$Hpil]
node 12 $Ltr [expr 5*$Hpil]
node 13 [expr 2*$Ltr] 0
node 14 [expr 2*$Ltr] $Hpil
node 15 [expr 2*$Ltr] [expr 2*$Hpil]
node 16 [expr 2*$Ltr] [expr 3*$Hpil]
node 17 [expr 2*$Ltr] [expr 4*$Hpil]
node 18 [expr 2*$Ltr] [expr 5*$Hpil]
node 19 [expr 3*$Ltr] 0
node 20 [expr 3*$Ltr] $Hpil
node 21 [expr 3*$Ltr] [expr 2*$Hpil]
node 22 [expr 3*$Ltr] [expr 3*$Hpil]
node 23 [expr 3*$Ltr] [expr 4*$Hpil]
node 24 [expr 3*$Ltr] [expr 5*$Hpil]

fix 1 1 1 1
fix 7 1 1 1
fix 13 1 1 1
fix 19 1 1 1

set m2 10.483
set m1 16.588

mass 2 $m1 $m1 0
mass 3 $m1 $m1 0
mass 4 $m1 $m1 0
mass 5 $m1 $m1 0
mass 6 $m2 $m2 0

mass 8 $m1 $m1 0
mass 9 $m1 $m1 0
mass 10 $m1 $m1 0
mass 11 $m1 $m1 0
mass 12 $m2 $m2 0

mass 14 $m1 $m1 0
mass 15 $m1 $m1 0
mass 16 $m1 $m1 0
mass 17 $m1 $m1 0
mass 18 $m2 $m2 0

mass 20 $m1 $m1 0
mass 21 $m1 $m1 0
mass 22 $m1 $m1 0
mass 23 $m1 $m1 0
mass 24 $m2 $m2 0

#MATERIALI
#calcestruzzo non confinato
set ClsUnconf 10
set fc -25.
set fcu [expr 0.2*$fc]
set Ec [expr 25000]
set epsc1 -0.002
set epscu -0.004

uniaxialMaterial Concrete01 $ClsUnconf $fc $epsc1 $fcu $epscu

#acciaio
set acciaio 18
set fy 430
set Es [expr 2.15e5]
set b 0.005

uniaxialMaterial Steel02 $acciaio $fy $Es $b 18 0.925 0.15

#CLS CONF1
set ClsConf1 11
set fcc1 -31.10
set epsc1c1 -0.0044
set fcuc1 [expr 0.2*$fcc1]
set epscuc1 -0.01

uniaxialMaterial Concrete01 $ClsConf1 $fcc1 $epsc1c1 $fcuc1 $epscuc1

#CLS CONF2
set ClsConf2 12
set fcc2 -30.12
set epsc1c2 -0.0044
set fcuc2 [expr 0.2*$fcc2]
set epscuc2 -0.0095

uniaxialMaterial Concrete01 $ClsConf2 $fcc2 $epsc1c2 $fcuc2 $epscuc2

#SEZIONI
set cf 30.

#PILASTRO ESTERNO
set SezPilEs1 1
set Hp1 400.
set Bp1 300.

set diamArm 16
set areaArm [expr 3.14*$diamArm**2/4]

set coverZ [expr $Bp1/2]
set coverY [expr $Hp1/2]
set coreY [expr $coverY-$cf]
set coreZ [expr $coverZ-$cf]
set fibreyCore 34;#parallele a z
set fibrezCore 24;#parallele a y
set fibreyCover 34
set fibrezCover 24
set fibraRec1 [expr $coreY-6]

section Fiber $SezPilEs1 {
patch quadr $ClsConf1 $fibrezCore $fibreyCore -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ -$coreY $coreZ;
patch quadr $ClsUnconf 2 $fibreyCore -$coverY -$coverZ -$coreY -$coreZ -$coreY $coreZ -$coverY $coverZ; # copriferro verticale
patch quadr $ClsUnconf 2 $fibreyCore $coreY -$coreZ $coverY -$coverZ $coverY $coverZ $coreY $coreZ
patch quadr $ClsUnconf $fibrezCore 2 -$coreY $coreZ $coreY $coreZ $coverY $coverZ -$coverY $coverZ; #copriferro orizzontale
patch quadr $ClsUnconf $fibrezCore 2 -$coverY -$coverZ $coverY -$coverZ $coreY -$coreZ -$coreY -$coreZ
layer straight $acciaio 2 $areaArm -$coreY $coreZ -$coreY -$coreZ;#armatura superiore
layer straight $acciaio 2 $areaArm $coreY $coreZ $coreY -$coreZ;#armatura inferiore
}

set SezPilEs2 2
set Hp2 350.
set Bp2 300.

set coverZ2 [expr $Bp2/2]
set coverY2 [expr $Hp2/2]
set coreY2 [expr $coverY2-$cf]
set coreZ2 [expr $coverZ2-$cf]
set fibreyCore2 29
set fibrezCore2 24
set fibreyCover2 29
set fibrezCover2 24
set fibraRec2 [expr $coreY2-6]

section Fiber $SezPilEs2 {
patch quadr $ClsConf1 $fibrezCore2 $fibreyCore2 -$coreY2 -$coreZ2 $coreY2 -$coreZ2 $coreY2 $coreZ2 -$coreY2 $coreZ2;
patch quadr $ClsUnconf 2 $fibreyCore2 -$coverY2 -$coverZ2 -$coreY2 -$coreZ2 -$coreY2 $coreZ2 -$coverY2 $coverZ2; # copriferro verticale
patch quadr $ClsUnconf 2 $fibreyCore2 $coreY2 -$coreZ2 $coverY2 -$coverZ2 $coverY2 $coverZ2 $coreY2 $coreZ2
patch quadr $ClsUnconf $fibrezCore2 2 -$coreY2 $coreZ2 $coreY2 $coreZ2 $coverY2 $coverZ2 -$coverY2 $coverZ2; #copriferro orizzontale
patch quadr $ClsUnconf $fibrezCore2 2 -$coverY2 -$coverZ2 $coverY2 -$coverZ2 $coreY2 -$coreZ2 -$coreY2 -$coreZ2
layer straight $acciaio 2 $areaArm -$coreY2 $coreZ2 -$coreY2 -$coreZ2;#armatura superiore
layer straight $acciaio 2 $areaArm $coreY2 $coreZ2 $coreY2 -$coreZ2;#armatura inferiore
}

#PILASTRI INTERNI
set SezPilInt1 3
set Hp3 600.
set Bp3 350.

set coverZ3 [expr $Bp3/2]
set coverY3 [expr $Hp3/2]
set coreY3 [expr $coverY3-$cf]
set coreZ3 [expr $coverZ3-$cf]
set fibreyCore3 54
set fibrezCore3 29
set fibreyCover3 54
set fibrezCover3 29
set fibraRec3 [expr $coreY3-9]

section Fiber $SezPilInt1 {
patch quadr $ClsConf2 $fibrezCore3 $fibreyCore3 -$coreY3 -$coreZ3 $coreY3 -$coreZ3 $coreY3 $coreZ3 -$coreY3 $coreZ3;
patch quadr $ClsUnconf 2 $fibreyCore3 -$coverY3 -$coverZ3 -$coreY3 -$coreZ3 -$coreY3 $coreZ3 -$coverY3 $coverZ3; # copriferro verticale
patch quadr $ClsUnconf 2 $fibreyCore3 $coreY3 -$coreZ3 $coverY3 -$coverZ3 $coverY3 $coverZ3 $coreY3 $coreZ3
patch quadr $ClsUnconf $fibrezCore3 2 -$coreY3 $coreZ3 $coreY3 $coreZ3 $coverY3 $coverZ3 -$coverY3 $coverZ3; #copriferro orizzontale
patch quadr $ClsUnconf $fibrezCore3 2 -$coverY3 -$coverZ3 $coverY3 -$coverZ3 $coreY3 -$coreZ3 -$coreY3 -$coreZ3
layer straight $acciaio 3 $areaArm -$coreY3 $coreZ3 -$coreY3 -$coreZ3;#armatura superiore
layer straight $acciaio 3 $areaArm $coreY3 $coreZ3 $coreY3 -$coreZ3;#armatura inferiore
}

set SezPilInt2 4
set Hp4 500.
set Bp4 350.

set coverZ4 [expr $Bp4/2]
set coverY4 [expr $Hp4/2]
set coreY4 [expr $coverY4-$cf]
set coreZ4 [expr $coverZ4-$cf]
set fibreyCore4 44
set fibrezCore4 29
set fibreyCover4 44
set fibrezCover4 29
set fibraRec4 [expr $coreY4-9]

section Fiber $SezPilInt2 {
patch quadr $ClsConf2 $fibrezCore4 $fibreyCore4 -$coreY4 -$coreZ4 $coreY4 -$coreZ4 $coreY4 $coreZ4 -$coreY4 $coreZ4;
patch quadr $ClsUnconf 2 $fibreyCore4 -$coverY4 -$coverZ4 -$coreY4 -$coreZ4 -$coreY4 $coreZ4 -$coverY4 $coverZ4; # copriferro verticale
patch quadr $ClsUnconf 2 $fibreyCore4 $coreY4 -$coreZ4 $coverY4 -$coverZ4 $coverY4 $coverZ4 $coreY4 $coreZ4
patch quadr $ClsUnconf $fibrezCore4 2 -$coreY4 $coreZ4 $coreY4 $coreZ4 $coverY4 $coverZ4 -$coverY4 $coverZ4; #copriferro orizzontale
patch quadr $ClsUnconf $fibrezCore4 2 -$coverY4 -$coverZ4 $coverY4 -$coverZ4 $coreY4 -$coreZ4 -$coreY4 -$coreZ4
layer straight $acciaio 3 $areaArm -$coreY4 $coreZ4 -$coreY4 -$coreZ4;#armatura superiore
layer straight $acciaio 3 $areaArm $coreY4 $coreZ4 $coreY4 -$coreZ4;#armatura inferiore
}

set SezPilInt3 5
set Hp5 500.
set Bp5 300.

set coverZ5 [expr $Bp5/2]
set coverY5 [expr $Hp5/2]
set coreY5 [expr $coverY5-$cf]
set coreZ5 [expr $coverZ5-$cf]
set fibreyCore5 44
set fibrezCore5 24
set fibreyCover5 44
set fibrezCover5 24
set fibraRec5 [expr $coreY5-9]

section Fiber $SezPilInt3 {
patch quadr $ClsConf1 $fibrezCore5 $fibreyCore5 -$coreY5 -$coreZ5 $coreY5 -$coreZ5 $coreY5 $coreZ5 -$coreY5 $coreZ5;
patch quadr $ClsUnconf 2 $fibreyCore5 -$coverY5 -$coverZ5 -$coreY5 -$coreZ5 -$coreY5 $coreZ5 -$coverY5 $coverZ5; # copriferro verticale
patch quadr $ClsUnconf 2 $fibreyCore5 $coreY5 -$coreZ5 $coverY5 -$coverZ5 $coverY5 $coverZ5 $coreY5 $coreZ5
patch quadr $ClsUnconf $fibrezCore5 2 -$coreY5 $coreZ5 $coreY5 $coreZ5 $coverY5 $coverZ5 -$coverY5 $coverZ5; #copriferro orizzontale
patch quadr $ClsUnconf $fibrezCore5 2 -$coverY5 -$coverZ5 $coverY5 -$coverZ5 $coreY5 -$coreZ5 -$coreY5 -$coreZ5
layer straight $acciaio 2 $areaArm -$coreY5 $coreZ5 -$coreY5 -$coreZ5;#armatura superiore
layer straight $acciaio 2 $areaArm $coreY5 $coreZ5 $coreY5 -$coreZ5;#armatura inferiore
}

set SezPilInt4 6
set Hp6 450.
set Bp6 300.

set coverZ6 [expr $Bp6/2]
set coverY6 [expr $Hp6/2]
set coreY6 [expr $coverY6-$cf]
set coreZ6 [expr $coverZ6-$cf]
set fibreyCore6 39
set fibrezCore6 24
set fibreyCover6 39
set fibrezCover6 24
set fibraRec6 [expr $coreY6-9]

section Fiber $SezPilInt4 {
patch quadr $ClsConf1 $fibrezCore6 $fibreyCore6 -$coreY6 -$coreZ6 $coreY6 -$coreZ6 $coreY6 $coreZ6 -$coreY6 $coreZ6;
patch quadr $ClsUnconf 2 $fibreyCore6 -$coverY6 -$coverZ6 -$coreY6 -$coreZ6 -$coreY6 $coreZ6 -$coverY6 $coverZ6; # copriferro verticale
patch quadr $ClsUnconf 2 $fibreyCore6 $coreY6 -$coreZ6 $coverY6 -$coverZ6 $coverY6 $coverZ6 $coreY6 $coreZ6
patch quadr $ClsUnconf $fibrezCore6 2 -$coreY6 $coreZ6 $coreY6 $coreZ6 $coverY6 $coverZ6 -$coverY6 $coverZ6; #copriferro orizzontale
patch quadr $ClsUnconf $fibrezCore6 2 -$coverY6 -$coverZ6 $coverY6 -$coverZ6 $coreY6 -$coreZ6 -$coreY6 -$coreZ6
layer straight $acciaio 2 $areaArm -$coreY6 $coreZ6 -$coreY6 -$coreZ6;#armatura superiore
layer straight $acciaio 2 $areaArm $coreY6 $coreZ6 $coreY6 -$coreZ6;#armatura inferiore
}


#TRAVE
set SezTrave 7
set Htr 500.
set Btr 300.
set M1 214.677e6
set M2 231.445e6
set M3 339.953e6
set C1 0.01e-3
set C2 0.04e-3
set C3 0.486e-3

set IDtraveM 16
set IDtraveN 17
set EA 1e10

uniaxialMaterial Hysteretic $IDtraveM $M1 $C1 $M2 $C2 $M3 $C3 -$M1 -$C1 -$M2 -$C2 -$M3 -$C3 1 1 0 0
uniaxialMaterial Elastic $IDtraveN $EA;
section Aggregator $SezTrave $IDtraveN P $IDtraveM Mz;

#ELEMENTI
geomTransf Linear 200;#colonne
geomTransf Linear 201;#trave
set CtrlSezP 5
set CtrlSezT 5

#pilastri esterni
element nonlinearBeamColumn 50 1 2 $CtrlSezP $SezPilEs1 200
element nonlinearBeamColumn 51 2 3 $CtrlSezP $SezPilEs1 200
element nonlinearBeamColumn 52 3 4 $CtrlSezP $SezPilEs2 200
element nonlinearBeamColumn 53 4 5 $CtrlSezP $SezPilEs2 200
element nonlinearBeamColumn 54 5 6 $CtrlSezP $SezPilEs2 200
element nonlinearBeamColumn 55 19 20 $CtrlSezP $SezPilEs1 200
element nonlinearBeamColumn 56 20 21 $CtrlSezP $SezPilEs1 200
element nonlinearBeamColumn 57 21 22 $CtrlSezP $SezPilEs2 200
element nonlinearBeamColumn 58 22 23 $CtrlSezP $SezPilEs2 200
element nonlinearBeamColumn 59 23 24 $CtrlSezP $SezPilEs2 200
#pilastri interni
element nonlinearBeamColumn 60 7 8 $CtrlSezP $SezPilInt1 200
element nonlinearBeamColumn 61 8 9 $CtrlSezP $SezPilInt2 200
element nonlinearBeamColumn 62 9 10 $CtrlSezP $SezPilInt3 200
element nonlinearBeamColumn 63 10 11 $CtrlSezP $SezPilInt4 200
element nonlinearBeamColumn 64 11 12 $CtrlSezP $SezPilInt4 200
element nonlinearBeamColumn 65 13 14 $CtrlSezP $SezPilInt1 200
element nonlinearBeamColumn 66 14 15 $CtrlSezP $SezPilInt2 200
element nonlinearBeamColumn 67 15 16 $CtrlSezP $SezPilInt3 200
element nonlinearBeamColumn 68 16 17 $CtrlSezP $SezPilInt4 200
element nonlinearBeamColumn 69 17 18 $CtrlSezP $SezPilInt4 200
#travi
element nonlinearBeamColumn 30 2 8 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 31 3 9 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 32 4 10 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 33 5 11 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 34 6 12 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 35 8 14 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 36 9 15 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 37 10 16 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 38 11 17 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 39 12 18 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 40 14 20 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 41 15 21 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 42 16 22 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 43 17 23 $CtrlSezT $SezTrave 201
element nonlinearBeamColumn 44 18 24 $CtrlSezT $SezTrave 201

#INFILL
#MATERIALE
set r 0.2
set pincheps 0.8;#regolano il ramo di ricarico: pincheps di quanto si riducono le deformazioni rispetto al caso senza pinching
set pinchsig 0.1;#pincheps di quanto si riducono le tensioni rispetto al caso senza pinching
set beta 0.8;

#TAMPONAMENTO A-E
set IDmuraturaA 21
set epsyA 2.16e-4
set sigmauA 2.65
set epsuA 8.9e-3
set sigmarA [expr $r*$sigmauA]

uniaxialMaterial Hysteretic $IDmuraturaA 0.000001 6e-6 0.00000001 6.5e-4 -$sigmauA -$epsyA -$sigmarA -$epsuA $pincheps $pinchsig 0 0 $beta

#TAMPONAMENTO B-C-D
set IDmuraturaC 22
set epsyC 1.89e-4
set sigmauC 2.32
set epsuC 8.09e-3
set sigmarC [expr $r*$sigmauC]

uniaxialMaterial Hysteretic $IDmuraturaC 0.000001 6e-6 0.00000001 6.5e-4 -$sigmauC -$epsyC -$sigmarC -$epsuC $pincheps $pinchsig 0 0 $beta

#SEZIONI
set fibrea 5
set fibreb 5
#TAMPONAMENTO A
set tA 240
set wA 623
set SezBiellaA 23

set aA [expr $tA/2]
set bA [expr $wA/2]
section Fiber $SezBiellaA {patch rect $IDmuraturaA $fibrea $fibreb -$bA -$aA $bA $aA}

#TAMPONAMENTO B
set tB 200
set wB 690
set SezBiellaB 24

set aB [expr $tB/2]
set bB [expr $wB/2]
section Fiber $SezBiellaB {patch rect $IDmuraturaC $fibrea $fibreb -$bB -$aB $bB $aB}

#TAMPONAMENTO C
set tC 150
set wC 697
set SezBiellaC 25

set aC [expr $tC/2]
set bC [expr $wC/2]
section Fiber $SezBiellaC {patch rect $IDmuraturaC $fibrea $fibreb -$bC -$aC $bC $aC}

#TAMPONAMENTO D
set tD 150
set wD 676
set SezBiellaD 26

set aD [expr $tD/2]
set bD [expr $wD/2]
section Fiber $SezBiellaD {patch rect $IDmuraturaC $fibrea $fibreb -$bD -$aD $bD $aD}

#TAMPONAMENTO E
set tE 150
set wE 676
set SezBiellaE 27

set aE [expr $tE/2]
set bE [expr $wE/2]
section Fiber $SezBiellaE {patch rect $IDmuraturaA $fibrea $fibreb -$bE -$aE $bE $aE}

element trussSection 120 1 8 $SezBiellaA
element trussSection 121 2 9 $SezBiellaA
element trussSection 122 3 10 $SezBiellaE
element trussSection 123 4 11 $SezBiellaE
element trussSection 124 5 12 $SezBiellaE
element trussSection 125 2 7 $SezBiellaA
element trussSection 126 3 8 $SezBiellaA
element trussSection 127 4 9 $SezBiellaE
element trussSection 128 5 10 $SezBiellaE
element trussSection 129 6 11 $SezBiellaE

element trussSection 130 7 14 $SezBiellaB
element trussSection 131 8 15 $SezBiellaB
element trussSection 132 9 16 $SezBiellaC
element trussSection 133 10 17 $SezBiellaD
element trussSection 134 11 18 $SezBiellaD
element trussSection 135 8 13 $SezBiellaB
element trussSection 136 9 14 $SezBiellaB
element trussSection 137 10 15 $SezBiellaC
element trussSection 138 11 16 $SezBiellaD
element trussSection 139 12 17 $SezBiellaD

element trussSection 140 13 20 $SezBiellaA
element trussSection 141 14 21 $SezBiellaA
element trussSection 142 15 22 $SezBiellaE
element trussSection 143 16 23 $SezBiellaE
element trussSection 144 17 24 $SezBiellaE
element trussSection 145 14 19 $SezBiellaA
element trussSection 146 15 20 $SezBiellaA
element trussSection 147 16 21 $SezBiellaE
element trussSection 148 17 22 $SezBiellaE
element trussSection 149 18 23 $SezBiellaE

puts "modello costruito"
recorder Node -file $cartel/D.out -time -node 6 -dof 1 disp;#spostamento per diagramma
recorder Node -file $cartel/D'.out -node 12 -dof 1 disp
recorder Node -file $cartel/taglioBase1.out -node 1 -dof 1 reaction
recorder Node -file $cartel/taglioBase7.out -node 7 -dof 1 reaction
recorder Node -file $cartel/taglioBase13.out -node 13 -dof 1 reaction
recorder Node -file $cartel/taglioBase19.out -node 19 -dof 1 reaction
recorder Drift -file $cartel/Drift.out -iNode 1 -jNode 6 -dof 1 -perpDirn 2
recorder Node -file $cartel/Mbase.out -node 1 7 13 19 -dof 3 reaction

#------------------GRAVITY LOAD--------------------------
#carico distribuito sulle travi
set qtraveCop 29.85
set qtravePiano 40.59

pattern Plain 100 Linear {
eleLoad -ele 34 39 44 -type -beamUniform -$qtraveCop
eleLoad -ele 30 31 32 33 35 36 37 38 40 41 42 43 -type -beamUniform -$qtravePiano
}


set Ppiles 8662.
set Ppilint 12375.

set qpiles [expr $Ppiles/$Hpil]
pattern Plain 101 Linear {
eleLoad -ele 50 51 52 53 54 55 56 57 58 59 -type -beamUniform 0 -$qpiles
}

set qpilint [expr $Ppilint/$Hpil]
pattern Plain 103 Linear {
eleLoad -ele 60 61 62 63 64 65 66 67 68 69 -type -beamUniform 0 -$qpilint
}

# Gravity-analysis parameters -- load-controlled static analysis
set Tol 1.0e-5; # convergence tolerance for test
variable constraintsTypeGravity Plain; # default;
if { [info exists RigidDiaphragm] == 1} {
if {$RigidDiaphragm=="ON"} {
variable constraintsTypeGravity Lagrange; # large model: try Transformation
}; # if rigid diaphragm is on
}; # if rigid diaphragm exists
constraints $constraintsTypeGravity ; # how it handles boundary conditions
numberer RCM; # renumber dof's to minimize band-width (optimization), if you want to
system UmfPack;#BandGeneral ; # how to store and solve the system of equations in the analysis (large model: try UmfPack)
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 "analisi carichi verticali fatta"

#---------------------PUSH OVER ANALYSIS------------------------------
set IDctrlNode 6
set IDctrlDOF 1

#analisi statica non lineare
set F2 $m1
set F3 $m1
set F4$m1
set F5 $m1
set F6 $m2
puts "le forze iniziali sono $F2 $F3 $F4 $F5 $F6"

#caratteristiche dell'analisi pushover
set driftmax 0.025
set Dmax 400;#[expr $driftmax*$Hpil*$Npiani]
set Nincr 400
set Dincr [expr $Dmax/$Nincr]

#carico
pattern Plain 102 Linear {
load 2 $F2 0.0 0.0
load 3 $F3 0.0 0.0
load 4 $F4 0.0 0.0
load 5 $F5 0.0 0.0
load 6 $F6 0.0 0.0
}

constraints Plain
numberer Plain
system BandGeneral
set Tol 1.e-6; # Convergence Test: tolerance
set maxNumIter 10; # Convergence Test: maximum number of iterations that will be performed before "failure to converge" is returned
set printFlag 0; # Convergence Test: flag used to print information on convergence (optional) # 1: print information on each step;
set TestType NormDispIncr; # Convergence-test type
test $TestType $Tol $maxNumIter $printFlag;
set algorithm Newton
algorithm $algorithm
integrator DisplacementControl $IDctrlNode $IDctrlDOF $Dincr
analysis Static

# --------------------------------- perform Static Pushover Analysis

set ok [analyze $Nincr]; # this will return zero if no convergence problems were encountered
if {$ok != 0} {puts "provo a rifare"}
if {$ok != 0} {
# if analysis fails, we try some other stuff, performance is slower inside this loop
set ok 0;
set controlDisp 0.0;
set D0 0.0; # analysis starts from zero
set Dstep [expr ($controlDisp-$D0)/($Dmax-$D0)]
while {$Dstep < 1.0 && $ok == 0} {
set controlDisp [nodeDisp $IDctrlNode $IDctrlDOF ]
set Dstep [expr ($controlDisp-$D0)/($Dmax-$D0)]
set ok [analyze 1 ]
if {$ok != 0} {
puts "Trying Newton with Initial Tangent .."
test NormDispIncr $Tol 2000 0
algorithm Newton -initial
set ok [analyze 1 ]
test $TestType $Tol $maxNumIter 0
algorithm $algorithm
}
if {$ok != 0} {
puts "Trying Broyden .."
algorithm Broyden 8
set ok [analyze 1 ]
algorithm $algorithm
}
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch .8
set ok [analyze 1 ]
algorithm $algorithm
}
}; # end while loop
}; # end if ok !0

if {$ok != 0} {puts "PROBLEMA Control Disp=[nodeDisp $IDctrlNode $IDctrlDOF]"}

if {$ok == 0} {puts "Pushover Done. Control Disp=[nodeDisp $IDctrlNode $IDctrlDOF]"}
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: problem with convergence

Post by vesna »

if you look at the norms of displacement at the point you encounter convergence problem you will notice that they do not change no matter what algorithm you use. inside of that analysis step the element (or elements) oscillate back and forth and can not reach convergence. this problem can not be fixed by changing the analysis objects. you will have to go to your model and fix elements that are causing the problem. You can identify the elements that are causing the problem by looking at delatX vector and finding components that are oscillating. To get deltaX vector you need to use flag 4 inside of "test" command.

I would also use more number of increments (Nincr) for pushover analysis (at least 10 times more).
sisa
Posts: 27
Joined: Sun Mar 13, 2011 8:03 am
Location: Bologna

Re: problem with convergence

Post by sisa »

thanks vesna, if I increase Nincr it runs!
I have another question...I've see that moments of beams are very different from right to left section..is it normal? I tried Exemple 5 of OpenSees manual and they aren't differen, so I think there is a problem in my model...do you have idea of what could it be? for example i put localForce recorder of beam 35 (first storey) where you can see that two moments are very different...
--------LOCAL FORCE BEAM 35----------------------------------------------
-1333.16 10147.5 8.46602e+006 1333.16 10147.5 -8.46602e+006
-2602.58 20295 1.69407e+007 2602.58 20295 -1.69407e+007
-3809 30442.5 2.5421e+007 3809 30442.5 -2.5421e+007
-4961.43 40590 3.39048e+007 4961.43 40590 -3.39048e+007
-6068.08 50737.5 4.23907e+007 6068.08 50737.5 -4.23907e+007
-7136.12 60885 5.08778e+007 7136.12 60885 -5.08778e+007
-8168.07 71032.5 5.93653e+007 8168.07 71032.5 -5.93653e+007
-9167.41 81180 6.78528e+007 9167.41 81180 -6.78528e+007
-10135.4 91327.5 7.63399e+007 10135.4 91327.5 -7.63399e+007
-11077.4 101475 8.48265e+007 11077.4 101475 -8.48265e+007
-10614.4 100863 8.33263e+007 10614.4 102087 -8.63848e+007
-10135.9 100250 8.18238e+007 10135.9 102700 -8.79488e+007
-9683.83 99641.7 8.03369e+007 9683.83 103308 -8.95032e+007
-9232.68 99031.6 7.88476e+007 9232.68 103918 -9.10646e+007
-8804.48 98421 7.73595e+007 8804.48 104529 -9.26297e+007
-8363.14 97804.8 7.58587e+007 8363.14 105145 -9.42096e+007
-7873.47 97188.3 7.43608e+007 7873.47 105762 -9.57941e+007
-7375.82 96569.6 7.28593e+007 7375.82 106380 -9.73862e+007
-6870.51 95949 7.13549e+007 6870.51 107001 -9.8985e+007
-6322.31 95333.2 6.98761e+007 6322.31 107617 -1.00585e+008
-5816.83 94690.2 6.82901e+007 5816.83 108260 -1.02214e+008
-5277.22 94003.4 6.6529e+007 5277.22 108947 -1.03887e+008
-4735.47 93322.9 6.4786e+007 4735.47 109627 -1.05547e+008
-4288.32 92712.8 6.3319e+007 4288.32 110237 -1.0713e+008
-3806.23 92095.2 6.18329e+007 3806.23 110855 -1.08732e+008
-3290.73 91469.8 6.0322e+007 3290.73 111480 -1.10348e+008
-2772.6 90842.4 5.88063e+007 2772.6 112108 -1.11969e+008
-2310.29 90203 5.72665e+007 2310.29 112747 -1.13626e+008
-1838.93 89561.6 5.57216e+007 1838.93 113388 -1.15288e+008
-1357.83 88917.1 5.41687e+007 1357.83 114033 -1.16958e+008
-862.139 88271.2 5.26123e+007 862.139 114679 -1.18631e+008
-347.652 87626.6 5.10599e+007 347.652 115323 -1.20302e+008
166.856 86979.1 4.95e+007 -166.856 115971 -1.2198e+008
673.646 86328.6 4.79328e+007 -673.646 116621 -1.23665e+008
1190.26 85675.9 4.63602e+007 -1190.26 117274 -1.25356e+008
1719.39 85017.8 4.4775e+007 -1719.39 117932 -1.27061e+008
2287.7 84354.5 4.31789e+007 -2287.7 118595 -1.28781e+008
2864.72 83686 4.15703e+007 -2864.72 119264 -1.30515e+008
3441 83012.1 3.99484e+007 -3441 119938 -1.32263e+008
4023.82 82330.8 3.83083e+007 -4023.82 120619 -1.34029e+008
4621.72 81646 3.66602e+007 -4621.72 121304 -1.35805e+008
5235.21 80955.8 3.49995e+007 -5235.21 121994 -1.37595e+008
5864.37 80261.2 3.33282e+007 -5864.37 122689 -1.39397e+008
6493.17 79562.9 3.16484e+007 -6493.17 123387 -1.41209e+008
7140.27 78860.3 2.99586e+007 -7140.27 124090 -1.43032e+008
7799.52 78154.6 2.82607e+007 -7799.52 124795 -1.44863e+008
8434.88 77414.9 2.64503e+007 -8434.88 125535 -1.46751e+008
9057.26 76662.1 2.46008e+007 -9057.26 126288 -1.48665e+008
9690.79 75908.1 2.27468e+007 -9690.79 127042 -1.50582e+008
10326.4 75152.3 2.0887e+007 -10326.4 127798 -1.52501e+008
10961 74393.2 1.90183e+007 -10961 128557 -1.54427e+008
11596 73631.8 1.71431e+007 -11596 129318 -1.56359e+008
12230.8 72865.3 1.52543e+007 -12230.8 130085 -1.58303e+008
12860.5 72088.8 1.3339e+007 -12860.5 130861 -1.6027e+008
13497.8 71308.6 1.14148e+007 -13497.8 131641 -1.62247e+008
14133.2 70524.7 9.48049e+006 -14133.2 132425 -1.64232e+008
14768.1 69737.6 7.53705e+006 -14768.1 133212 -1.66224e+008
15406.6 68951.1 5.59725e+006 -15406.6 133999 -1.68217e+008
16055.1 68163.1 3.6554e+006 -16055.1 134787 -1.70215e+008
16708.1 67379.6 1.72983e+006 -16708.1 135570 -1.72207e+008
17352 66628.5 -82256.5 -17352 136321 -1.7415e+008
18000.4 65875.2 -1.90078e+006 -18000.4 137075 -1.76098e+008
18639.4 65119.4 -3.72434e+006 -18639.4 137831 -1.78054e+008
19283.9 64365.3 -5.54387e+006 -19283.9 138585 -1.80005e+008
19931 63608.9 -7.37088e+006 -19931 139341 -1.8196e+008
20587.5 62851.6 -9.20229e+006 -20587.5 140098 -1.83915e+008
21246.9 62093 -1.10393e+007 -21246.9 140857 -1.85871e+008
21889.4 61334.9 -1.28763e+007 -21889.4 141615 -1.87824e+008
22511.9 60569.8 -1.47298e+007 -22511.9 142380 -1.89796e+008
23127.1 59803.9 -1.65876e+007 -23127.1 143146 -1.91768e+008
23759.5 59039.8 -1.8445e+007 -23759.5 143910 -1.93731e+008
24395.8 58275.4 -2.0306e+007 -24395.8 144675 -1.95692e+008
25029 57513.1 -2.21656e+007 -25029 145437 -1.97644e+008
25652.5 56753.1 -2.40233e+007 -25652.5 146197 -1.99586e+008
26268 55993.4 -2.5884e+007 -26268 146957 -2.01524e+008
26883.6 55236.1 -2.77442e+007 -26883.6 147714 -2.0345e+008
27489.4 54484.3 -2.9598e+007 -27489.4 148466 -2.05356e+008
28091.4 53738.7 -3.14448e+007 -28091.4 149211 -2.07237e+008
28680.6 53000.5 -3.32847e+007 -28680.6 149949 -2.09088e+008
29265.1 52273.8 -3.5111e+007 -29265.1 150676 -2.10895e+008
29847.3 51563 -3.69169e+007 -29847.3 151387 -2.12643e+008
30421.9 50874.4 -3.86932e+007 -30421.9 152076 -2.1431e+008
31000.5 50212.9 -4.04333e+007 -31000.5 152737 -2.15877e+008
31595.9 49615.9 -4.20801e+007 -31595.9 153334 -2.17215e+008
32189.8 49076 -4.36438e+007 -32189.8 153874 -2.18351e+008
32780.4 48577.4 -4.51499e+007 -32780.4 154373 -2.19338e+008
33384.4 48095.1 -4.66507e+007 -33384.4 154855 -2.20249e+008
33970 47622.9 -4.81424e+007 -33970 155327 -2.21118e+008
34550.7 47153.3 -4.96355e+007 -34550.7 155797 -2.21973e+008
35133.3 46686.6 -5.1123e+007 -35133.3 156263 -2.22819e+008
35720.9 46220.4 -5.26124e+007 -35720.9 156730 -2.23661e+008
36318.8 45753.7 -5.4106e+007 -36318.8 157196 -2.24501e+008
36949.7 45340.1 -5.5523e+007 -36949.7 157610 -2.25151e+008
37578.6 44988.4 -5.68516e+007 -37578.6 157962 -2.25581e+008
38132.9 44634.2 -5.81858e+007 -38132.9 158316 -2.26018e+008
38646.2 44278.1 -5.95344e+007 -38646.2 158672 -2.2645e+008
39165 43909.3 -6.09369e+007 -39165 159041 -2.26892e+008
39643.7 43537.5 -6.23572e+007 -39643.7 159413 -2.2733e+008
40113.6 43164.2 -6.37841e+007 -40113.6 159786 -2.2777e+008
40535.1 42790.3 -6.52224e+007 -40535.1 160160 -2.28201e+008
40969.3 42411.4 -6.66701e+007 -40969.3 160539 -2.28648e+008
41420.3 42110.4 -6.80125e+007 -41420.3 160840 -2.2881e+008
41918.5 41916.8 -6.92044e+007 -41918.5 161033 -2.28586e+008
42381.6 41711.6 -7.0429e+007 -42381.6 161238 -2.28388e+008
42791.5 41447.8 -7.1757e+007 -42791.5 161502 -2.28379e+008
43114.8 41121.3 -7.31805e+007 -43114.8 161829 -2.28588e+008
43463.7 40733.1 -7.47643e+007 -43463.7 162217 -2.28945e+008
43787.8 40310.8 -7.64224e+007 -43787.8 162639 -2.29399e+008
44076.8 39859.4 -7.81512e+007 -44076.8 163091 -2.29927e+008
44332.4 39405.3 -7.98937e+007 -44332.4 163545 -2.30455e+008
44494.8 38982.8 -8.1558e+007 -44494.8 163967 -2.30903e+008
44694.7 38555.4 -8.3239e+007 -44694.7 164395 -2.31359e+008
44866.3 38123.5 -8.49417e+007 -44866.3 164826 -2.31816e+008
45009.2 37685.1 -8.66737e+007 -45009.2 165265 -2.32276e+008
45153.7 37236.1 -8.84449e+007 -45153.7 165714 -2.32749e+008
45552.3 36823.2 -9.01392e+007 -45552.3 166127 -2.3312e+008
45782.8 36436.3 -9.17851e+007 -45782.8 166514 -2.33408e+008
46015.5 36047.8 -9.34387e+007 -46015.5 166902 -2.33697e+008
46354.8 35655.5 -9.50903e+007 -46354.8 167294 -2.34007e+008
46588 35267.2 -9.67421e+007 -46588 167683 -2.34297e+008
46780.8 34872 -9.84125e+007 -46780.8 168078 -2.34602e+008
46961.2 34478.9 -1.00075e+008 -46961.2 168471 -2.34906e+008
47139.8 34078.7 -1.01769e+008 -47139.8 168871 -2.35213e+008
47328.5 33674.9 -1.03478e+008 -47328.5 169275 -2.35522e+008
47697.1 33276.2 -1.05141e+008 -47697.1 169674 -2.35853e+008
48067.9 32876.9 -1.06807e+008 -48067.9 170073 -2.36184e+008
48422.9 32476.8 -1.08479e+008 -48422.9 170473 -2.36512e+008
48778 32078.3 -1.10144e+008 -48778 170872 -2.36839e+008
49191.1 31681.4 -1.11794e+008 -49191.1 171269 -2.37173e+008
49598.4 31284.4 -1.13446e+008 -49598.4 171666 -2.37507e+008
50009.5 30886.6 -1.15102e+008 -50009.5 172063 -2.3784e+008
50421.5 30488.3 -1.16759e+008 -50421.5 172462 -2.38174e+008
50839.5 30088.4 -1.18423e+008 -50839.5 172862 -2.3851e+008
51248.5 29688.1 -1.20091e+008 -51248.5 173262 -2.38843e+008
51672.4 29297.4 -1.21716e+008 -51672.4 173653 -2.39172e+008
52080.8 28907.1 -1.23342e+008 -52080.8 174043 -2.39498e+008
52480.9 28517.1 -1.24968e+008 -52480.9 174433 -2.39822e+008
52883.8 28122.9 -1.26611e+008 -52883.8 174827 -2.40149e+008
53260.4 27726.3 -1.2827e+008 -53260.4 175224 -2.40473e+008
53629.6 27328.3 -1.29937e+008 -53629.6 175622 -2.40796e+008
54003.2 26927.7 -1.31615e+008 -54003.2 176022 -2.41122e+008
54370.7 26526.3 -1.33297e+008 -54370.7 176424 -2.41446e+008
54728.7 26121.7 -1.34996e+008 -54728.7 176828 -2.41771e+008
55079.3 25715.9 -1.36701e+008 -55079.3 177234 -2.42094e+008
55432.7 25310.1 -1.38407e+008 -55432.7 177640 -2.42418e+008
55781.1 24903.4 -1.40117e+008 -55781.1 178047 -2.42741e+008
56079.1 24494.3 -1.41847e+008 -56079.1 178456 -2.43057e+008
56380.7 24080.9 -1.43595e+008 -56380.7 178869 -2.43375e+008
56678.3 23666.7 -1.45348e+008 -56678.3 179283 -2.43693e+008
56889.8 23255.4 -1.47108e+008 -56889.8 179695 -2.4399e+008
56968 22841.8 -1.48902e+008 -56968 180108 -2.44264e+008
57090.2 22420.6 -1.50726e+008 -57090.2 180529 -2.44547e+008
56658.9 22001.4 -1.52634e+008 -56658.9 180949 -2.44733e+008
53529.6 21599.2 -1.54929e+008 -53529.6 181351 -2.4445e+008
52775.9 21123.1 -1.57151e+008 -52775.9 181827 -2.44609e+008
54038.5 20383.6 -1.60204e+008 -54038.5 182566 -2.45254e+008
55194.8 19759.2 -1.62774e+008 -55194.8 183191 -2.45805e+008
55652.1 19255.3 -1.6493e+008 -55652.1 183695 -2.46169e+008
56147.1 18746.7 -1.67102e+008 -56147.1 184203 -2.4654e+008
56377.2 18270.3 -1.69236e+008 -56377.2 184680 -2.46787e+008
56844.4 17791.5 -1.71353e+008 -56844.4 185159 -2.47065e+008
55782.3 17378.2 -1.73447e+008 -55782.3 185572 -2.47037e+008
55182.8 16923 -1.75641e+008 -55182.8 186027 -2.47119e+008
55364.5 16506.8 -1.77515e+008 -55364.5 186443 -2.47326e+008
55300.1 16043.1 -1.79646e+008 -55300.1 186907 -2.47514e+008
54976.9 15692.3 -1.81315e+008 -54976.9 187258 -2.47599e+008
55276.2 15281.9 -1.83137e+008 -55276.2 187668 -2.47828e+008
55209.1 14877.4 -1.84999e+008 -55209.1 188073 -2.47989e+008
55407.9 14508.4 -1.86656e+008 -55407.9 188442 -2.48177e+008
55379.1 14153.2 -1.88291e+008 -55379.1 188797 -2.48318e+008
55658.4 13747.3 -1.90097e+008 -55658.4 189203 -2.48541e+008
55737.6 13366.2 -1.91825e+008 -55737.6 189584 -2.48719e+008
55744.3 13022.3 -1.93402e+008 -55744.3 189928 -2.48862e+008
55932.9 12666 -1.95001e+008 -55932.9 190284 -2.49044e+008
56328.8 12307.2 -1.96573e+008 -56328.8 190643 -2.49266e+008
56739.4 11920.5 -1.98267e+008 -56739.4 191030 -2.49505e+008
56401.5 11575.1 -1.99903e+008 -56401.5 191375 -2.49597e+008
56411.4 11253.1 -2.01376e+008 -56411.4 191697 -2.49734e+008
56654.4 10922.3 -2.02845e+008 -56654.4 192028 -2.49918e+008
56896.2 10583.8 -2.04348e+008 -56896.2 192366 -2.50108e+008
57322.1 10218.2 -2.05938e+008 -57322.1 192732 -2.50346e+008
57436.3 9876.14 -2.07473e+008 -57436.3 193074 -2.50521e+008
57334.9 9575.26 -2.0886e+008 -57334.9 193375 -2.50639e+008
57178 9273.5 -2.10261e+008 -57178 193676 -2.50747e+008
57474.3 8953 -2.11665e+008 -57474.3 193997 -2.50945e+008
57679.5 8640.52 -2.13047e+008 -57679.5 194309 -2.51125e+008
57883.8 8328.07 -2.14425e+008 -57883.8 194622 -2.5131e+008
58208.3 8011.54 -2.15799e+008 -58208.3 194938 -2.51518e+008
58476 7710.3 -2.1711e+008 -58476 195240 -2.51713e+008
58669.6 7431.4 -2.18329e+008 -58669.6 195519 -2.51889e+008
58833.2 7160.41 -2.19512e+008 -58833.2 195790 -2.5206e+008
58715.3 6916.48 -2.20619e+008 -58715.3 196034 -2.52174e+008
59011.7 6673.23 -2.21643e+008 -59011.7 196277 -2.52366e+008
59433.2 6431.9 -2.22629e+008 -59433.2 196518 -2.52587e+008
59692.9 6221.62 -2.23498e+008 -59692.9 196728 -2.52769e+008
59960 6024.03 -2.24303e+008 -59960 196926 -2.52952e+008
60392.4 5817.94 -2.25111e+008 -60392.4 197132 -2.53175e+008
60554 5645.63 -2.2581e+008 -60554 197304 -2.53337e+008
60669.2 5484.91 -2.26463e+008 -60669.2 197465 -2.53488e+008
60888.7 5325.54 -2.27091e+008 -60888.7 197624 -2.53656e+008
61147.5 5163.5 -2.27722e+008 -61147.5 197786 -2.53835e+008
61328.7 5004.6 -2.28346e+008 -61328.7 197945 -2.54006e+008
61425.9 4850.8 -2.28957e+008 -61425.9 198099 -2.54164e+008
61529.2 4696.49 -2.29566e+008 -61529.2 198254 -2.54327e+008
61707.9 4541.42 -2.30161e+008 -61707.9 198409 -2.54507e+008
62095.4 4427.77 -2.30514e+008 -62095.4 198522 -2.54722e+008
62246.9 4324.63 -2.30861e+008 -62246.9 198625 -2.54891e+008
62271.7 4232.72 -2.31179e+008 -62271.7 198717 -2.55032e+008
62278.8 4143.19 -2.3149e+008 -62278.8 198807 -2.55169e+008
62393 4044.14 -2.31824e+008 -62393 198906 -2.5533e+008
62606.1 3935.32 -2.32181e+008 -62606.1 199015 -2.55517e+008
62750.6 3828.5 -2.32533e+008 -62750.6 199122 -2.55699e+008
62909.7 3722.35 -2.32883e+008 -62909.7 199228 -2.5588e+008
63198.7 3608.19 -2.33255e+008 -63198.7 199342 -2.56079e+008
63561.9 3484.01 -2.33667e+008 -63561.9 199466 -2.56288e+008
63940.2 3492.22 -2.3344e+008 -63940.2 199458 -2.56474e+008
64399.9 3490.99 -2.33247e+008 -64399.9 199459 -2.56673e+008
64881.3 3477.25 -2.33111e+008 -64881.3 199473 -2.56877e+008
65528.9 3394.51 -2.3328e+008 -65528.9 199555 -2.57122e+008
65956.2 3295.55 -2.33569e+008 -65956.2 199654 -2.57328e+008
66145.9 3196.22 -2.33906e+008 -66145.9 199754 -2.57488e+008
66327.4 3082.77 -2.34312e+008 -66327.4 199867 -2.57649e+008
66458 2963.67 -2.34758e+008 -66458 199986 -2.57799e+008
66573.1 2843.3 -2.35213e+008 -66573.1 200107 -2.57945e+008
66695.4 2735.32 -2.35607e+008 -66695.4 200215 -2.58092e+008
66741.8 2634.23 -2.35981e+008 -66741.8 200316 -2.58223e+008
66770.9 2533.74 -2.36355e+008 -66770.9 200416 -2.58352e+008
66814.2 2431.25 -2.36734e+008 -66814.2 200519 -2.58485e+008
66829.5 2337.43 -2.37073e+008 -66829.5 200613 -2.58615e+008
66825.3 2259.08 -2.37336e+008 -66825.3 200691 -2.58744e+008
66814.8 2179.72 -2.37603e+008 -66814.8 200770 -2.58874e+008
66717.3 2125.21 -2.37772e+008 -66717.3 200825 -2.58977e+008
66564.4 2079.91 -2.37908e+008 -66564.4 200870 -2.59068e+008
66360.3 2037.76 -2.38037e+008 -66360.3 200912 -2.59149e+008
65848.2 2023.87 -2.38089e+008 -65848.2 200926 -2.59167e+008
65780.6 1978.72 -2.38207e+008 -65780.6 200971 -2.59275e+008
65616.8 1939.79 -2.38312e+008 -65616.8 201010 -2.59364e+008
65279.9 1914.01 -2.38384e+008 -65279.9 201036 -2.59421e+008
64701.1 1908.27 -2.38405e+008 -64701.1 201042 -2.59429e+008
64719.8 1850.8 -2.38564e+008 -64719.8 201099 -2.59557e+008
64696.5 1797.11 -2.38713e+008 -64696.5 201153 -2.59676e+008
64667.8 1744.04 -2.38861e+008 -64667.8 201206 -2.59794e+008
64628.3 1691.9 -2.39007e+008 -64628.3 201258 -2.59909e+008
64565.8 1641.63 -2.39148e+008 -64565.8 201308 -2.60019e+008
64467.7 1594.67 -2.39281e+008 -64467.7 201355 -2.60121e+008
64338.7 1550.83 -2.39407e+008 -64338.7 201399 -2.60214e+008
64207.7 1507.63 -2.39531e+008 -64207.7 201442 -2.60306e+008
64040.6 1467.93 -2.39646e+008 -64040.6 201482 -2.60389e+008
63858.8 1429.61 -2.39759e+008 -63858.8 201520 -2.60468e+008
63641.4 1394.12 -2.39864e+008 -63641.4 201556 -2.6054e+008
63484 1353.69 -2.39982e+008 -63484 201596 -2.60624e+008
63313.6 1314.6 -2.40098e+008 -63313.6 201635 -2.60704e+008
63323.7 1259.86 -2.40253e+008 -63323.7 201690 -2.60823e+008
63179.6 1216.98 -2.40381e+008 -63179.6 201733 -2.60909e+008
63026.3 1174.03 -2.40507e+008 -63026.3 201776 -2.60998e+008
62885.7 1128.96 -2.40637e+008 -62885.7 201821 -2.61093e+008
62769 1081.4 -2.40771e+008 -62769 201869 -2.61197e+008
62629.6 1036.78 -2.40899e+008 -62629.6 201913 -2.61292e+008
62474.1 991.543 -2.41027e+008 -62474.1 201958 -2.6139e+008
62322.3 946.404 -2.41155e+008 -62322.3 202004 -2.61488e+008
62231.4 898.128 -2.41293e+008 -62231.4 202052 -2.61591e+008
62145.3 849.589 -2.41432e+008 -62145.3 202100 -2.61695e+008
62062.4 800.289 -2.41571e+008 -62062.4 202150 -2.61802e+008
62008.4 747.718 -2.41718e+008 -62008.4 202202 -2.61918e+008
61991.2 691.126 -2.41874e+008 -61991.2 202259 -2.62046e+008
61994.5 631.9 -2.42034e+008 -61994.5 202318 -2.62181e+008
62008.9 571.574 -2.42196e+008 -62008.9 202378 -2.62322e+008
61826.7 526.295 -2.42314e+008 -61826.7 202424 -2.6243e+008
61630 482.215 -2.42429e+008 -61630 202468 -2.62535e+008
60367.1 519.546 -2.4231e+008 -60367.1 202430 -2.62467e+008
59252.9 546.002 -2.42222e+008 -59252.9 202404 -2.62423e+008
58671.8 531.457 -2.42252e+008 -58671.8 202419 -2.62465e+008
58567.6 477.338 -2.42395e+008 -58567.6 202473 -2.62593e+008
59434.8 321.648 -2.42821e+008 -59434.8 202628 -2.62945e+008
60127 192.24 -2.43173e+008 -60127 202758 -2.63241e+008
60639.4 88.5157 -2.43454e+008 -60639.4 202861 -2.63478e+008
61031.8 -6.88383 -2.4371e+008 -61031.8 202957 -2.63699e+008
60993.9 -63.8124 -2.43858e+008 -60993.9 203014 -2.63836e+008
61406.1 -155.777 -2.44107e+008 -61406.1 203106 -2.64047e+008
61843.3 -250.03 -2.44363e+008 -61843.3 203200 -2.64263e+008
61904.2 -312.226 -2.44527e+008 -61904.2 203262 -2.64409e+008
61914.2 -369.998 -2.44678e+008 -61914.2 203320 -2.64547e+008
62477.1 -475.021 -2.44962e+008 -62477.1 203425 -2.64789e+008
62682.6 -547.533 -2.45155e+008 -62682.6 203498 -2.64958e+008
62881.7 -620.864 -2.45348e+008 -62881.7 203571 -2.65131e+008
63259.6 -707.857 -2.45577e+008 -63259.6 203658 -2.65337e+008
63755 -804.429 -2.45829e+008 -63755 203754 -2.65568e+008
63397.3 -832.101 -2.45891e+008 -63397.3 203782 -2.65645e+008
63813.1 -919.365 -2.46118e+008 -63813.1 203869 -2.65854e+008
64124.5 -999.157 -2.46325e+008 -64124.5 203949 -2.66046e+008
64526.7 -1086.43 -2.46551e+008 -64526.7 204036 -2.66256e+008
65426.3 -1214.46 -2.46885e+008 -65426.3 204164 -2.66563e+008
65647.8 -1286.4 -2.4707e+008 -65647.8 204236 -2.66737e+008
65923 -1363.42 -2.4727e+008 -65923 204313 -2.66922e+008
66148.2 -1436.85 -2.47454e+008 -66148.2 204387 -2.67105e+008
66375.2 -1504.46 -2.47618e+008 -66375.2 204454 -2.67279e+008
66695.7 -1574.97 -2.47776e+008 -66695.7 204525 -2.67474e+008
67316.5 -1673.93 -2.48016e+008 -67316.5 204624 -2.67729e+008
67569 -1738.53 -2.48155e+008 -67569 204689 -2.67912e+008
67757 -1798.98 -2.48288e+008 -67757 204749 -2.68082e+008
68025 -1865.72 -2.48438e+008 -68025 204816 -2.68266e+008
68462 -1946.46 -2.48626e+008 -68462 204896 -2.68481e+008
68632.4 -2004.89 -2.48753e+008 -68632.4 204955 -2.68647e+008
68816.2 -2065.26 -2.48883e+008 -68816.2 205015 -2.68818e+008
68946.5 -2121.94 -2.49001e+008 -68946.5 205072 -2.68983e+008
69149.8 -2184.85 -2.49136e+008 -69149.8 205135 -2.69163e+008
69538.7 -2260.76 -2.49314e+008 -69538.7 205211 -2.69365e+008
69715.7 -2319.18 -2.4944e+008 -69715.7 205269 -2.69531e+008
69987.3 -2385.79 -2.4959e+008 -69987.3 205336 -2.69714e+008
70196.3 -2447.2 -2.49725e+008 -70196.3 205397 -2.69886e+008
70410.3 -2509.2 -2.49861e+008 -70410.3 205459 -2.7006e+008
70712.1 -2578.43 -2.50018e+008 -70712.1 205528 -2.7025e+008
70865.6 -2634.86 -2.50139e+008 -70865.6 205585 -2.70411e+008
71059.6 -2693.86 -2.50265e+008 -71059.6 205644 -2.70579e+008
71312.4 -2756.07 -2.50398e+008 -71312.4 205706 -2.70758e+008
71773.3 -2836.25 -2.50582e+008 -71773.3 205786 -2.70974e+008
71962.7 -2892.67 -2.50699e+008 -71962.7 205843 -2.71139e+008
72052.6 -2940.33 -2.50792e+008 -72052.6 205890 -2.71285e+008
72219.7 -2994.72 -2.50904e+008 -72219.7 205945 -2.71445e+008
72409.7 -3051.08 -2.51021e+008 -72409.7 206001 -2.71609e+008
72599.5 -3107.51 -2.51138e+008 -72599.5 206058 -2.71774e+008
72787.9 -3163.84 -2.51255e+008 -72787.9 206114 -2.71939e+008
72748.3 -3201.11 -2.51316e+008 -72748.3 206151 -2.72065e+008
73126.8 -3273.46 -2.51476e+008 -73126.8 206223 -2.72266e+008
73465.8 -3342.35 -2.51632e+008 -73465.8 206292 -2.72455e+008
73761.7 -3409.23 -2.51779e+008 -73761.7 206359 -2.72642e+008
73976.3 -3469.14 -2.51907e+008 -73976.3 206419 -2.72814e+008
74154.9 -3524.99 -2.52024e+008 -74154.9 206475 -2.72976e+008
74324.8 -3579.88 -2.52139e+008 -74324.8 206530 -2.73136e+008
74488.3 -3634.02 -2.52251e+008 -74488.3 206584 -2.73294e+008
74639.6 -3688.54 -2.52365e+008 -74639.6 206639 -2.73453e+008
74688.3 -3732.42 -2.5245e+008 -74688.3 206682 -2.73587e+008
74871 -3787.72 -2.52566e+008 -74871 206738 -2.73748e+008
75148.6 -3851.15 -2.52705e+008 -75148.6 206801 -2.73926e+008
75425.9 -3914.55 -2.52843e+008 -75425.9 206865 -2.74105e+008
75742.5 -3981.39 -2.52991e+008 -75742.5 206931 -2.74291e+008
75941.2 -4039.04 -2.53112e+008 -75941.2 206989 -2.74458e+008
76164.1 -4099.12 -2.5324e+008 -76164.1 207049 -2.74631e+008
76380.8 -4158.46 -2.53365e+008 -76380.8 207108 -2.74802e+008
76551.2 -4213.61 -2.5348e+008 -76551.2 207164 -2.74963e+008
76686.4 -4265.64 -2.53586e+008 -76686.4 207216 -2.75117e+008
76797.1 -4315.13 -2.53685e+008 -76797.1 207265 -2.75266e+008
76865.7 -4361.04 -2.53773e+008 -76865.7 207311 -2.75407e+008
76860.8 -4401.26 -2.53844e+008 -76860.8 207351 -2.75537e+008
77066.1 -4459.4 -2.53966e+008 -77066.1 207409 -2.75706e+008
77328.6 -4521.47 -2.54101e+008 -77328.6 207471 -2.75881e+008
77389.8 -4566.86 -2.54189e+008 -77389.8 207517 -2.76021e+008
77438.1 -4611.17 -2.54273e+008 -77438.1 207561 -2.76158e+008
77522.1 -4658.67 -2.54365e+008 -77522.1 207609 -2.76303e+008
77610.1 -4706.47 -2.54459e+008 -77610.1 207656 -2.76448e+008
77711.3 -4755.25 -2.54555e+008 -77711.3 207705 -2.76596e+008
77824.8 -4804.99 -2.54654e+008 -77824.8 207755 -2.76746e+008
77948.9 -4855.63 -2.54756e+008 -77948.9 207806 -2.76897e+008
78075.2 -4906.45 -2.54858e+008 -78075.2 207856 -2.77049e+008
78207.8 -4957.86 -2.54962e+008 -78207.8 207908 -2.77202e+008
78366.1 -5012.1 -2.55073e+008 -78366.1 207962 -2.77362e+008
78532.6 -5066.54 -2.55185e+008 -78532.6 208017 -2.77523e+008
79287.5 -5178.25 -2.55455e+008 -79287.5 208128 -2.77811e+008
79365.7 -5223.85 -2.55542e+008 -79365.7 208174 -2.77952e+008
79347.1 -5261.42 -2.55607e+008 -79347.1 208211 -2.78075e+008
79302.8 -5296.85 -2.55666e+008 -79302.8 208247 -2.78193e+008
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: problem with convergence

Post by vesna »

Based on your results it looks like you apply gravity on beams as element load (in the first step of the output the moment is the same with the opposite sign). Once you start applying lateral load moments will start changing since at one end you are adding and at the other end you are subtracting.
sisa
Posts: 27
Joined: Sun Mar 13, 2011 8:03 am
Location: Bologna

Re: problem with convergence

Post by sisa »

you're right vesna...the problem is that my professor isn't sure that results are correct so I'm tring to control that there aren't error. also in columns moments at the bottom and at the top are very different and during analysis (fist apply gravity load and than lateral loads) axial loads in columns change a lot. is it normal? is it possible to make the axial load costant during analysis in order to control that the results are the same obtained with an other software? really I dont' understand if there is a problem in the model, I think no, but with STRAUSS (an other software) the pushover curve results differet and the difference is in the moment of columns that with STRAUSS are very similar at the bottom and at the top...
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: problem with convergence

Post by vesna »

With this command you are keeping gravity load constant (copied from your code):

# ------------------------------------------------- maintain constant gravity loads and reset time to zero
loadConst -time 0.0

The moments in the column at the top and bottom are usually not that different. I would suggest you to make very simple model: one-storey one-bay frame in order to test your model. I would also do the analysis without gravity to test the model for lateral load only.
Post Reply