Hello, I have a problem.
I'm doing a nonlinear dynamic analysis with accelerogram.
During the analysis, I get:
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.16
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.16
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.17
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.17
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.2
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.2
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.21
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.21
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.94
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.94
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.95
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 20.95
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 21.01
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 21.37
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 21.68
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 21.68
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 21.91
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 21.91
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 21.92
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 21.92
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 22.17
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 22.17
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 23.19
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 23.31
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 24.09
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 24.1
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 24.1
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 25
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 25.04
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 25.06
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 25.09
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 25.09
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 25.11
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 25.13
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 25.15
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 26.5
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 26.5
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 27.08
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 27.16
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 29.54
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 30.53
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 30.53
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 30.55
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 30.62
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 30.62
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 31.15
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 31.15
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 31.2
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 31.2
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 31.31
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 31.31
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 31.32
OpenSees > analyze failed, returned: -3 error flag
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 31.32
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
Ground Motion Done. End Time: 34.860000
---------------------------------------------------------------------------
End of script <D:\Gallo\Tesi\OpenSees2.2.2.exe> reached, Press any key to contin
ue
Carry only the last part of the analysis.
The length of time history has a duration of 34.85s with a step of 0.01s.
The results I get are affected by error??
Thanks.
CTestEnergyIncr::test() - failed to converge
Moderators: silvia, selimgunay, Moderators
-
GALLO ANTONIO
- Posts: 23
- Joined: Tue Nov 16, 2010 6:14 am
- Location: Sapienza
Re: CTestEnergyIncr::test() - failed to converge
Yes, at times that are shown (20.16, 20.17,20.2,....) you did not achieve the convergence. Thus, from time 20.16 and on, your results are not reliable.
-
GALLO ANTONIO
- Posts: 23
- Joined: Tue Nov 16, 2010 6:14 am
- Location: Sapienza
Re: CTestEnergyIncr::test() - failed to converge
Thanks for your reply.
So how can I solve the problem of convergence?
Thanks.
So how can I solve the problem of convergence?
Thanks.
Re: CTestEnergyIncr::test() - failed to converge
For this you can make some changes in your model (if you are using materials with a sharp transition from elastic to plastic use materials that have smooth transition, and so on...), or if you want to keep your model as it is you will have to play using different algorithms, changing number of iteration for a step, changing analyzing time-step increment (dt), maybe tolerance, and so on...
-
GALLO ANTONIO
- Posts: 23
- Joined: Tue Nov 16, 2010 6:14 am
- Location: Sapienza
Re: CTestEnergyIncr::test() - failed to converge
Thanks for your reply.
Gentilmemte if I send you my model
give me a hand to solve the problem?
wipe; # clear memory of all past model definitions
model BasicBuilder -ndm 2 -ndf 3; # Define the model builder, ndm=#dimension, ndf=#dofs
set dataDir TelaioFiberRCSectionEQ; # set up name of data directory (you can remove this)
file mkdir $dataDir; # create data directory
set GMdir "../GMfiles/"; # ground-motion file directory
# define UNITS ----------------------------------------------------------------------------
set m 1.; # define basic units -- output units
set kN 1.; # define basic units -- output units
set sec 1.; # define basic units -- output units
set LunitTXT "m"; # define basic-unit text for output
set FunitTXT "kN"; # define basic-unit text for output
set TunitTXT "sec"; # define basic-unit text for output
set kPa [expr $kN/pow($m,2)];
set ps [expr $kN/pow($m,3)]; # peso specifico
set m2 [expr $m*$m]; # m^2
set m4 [expr $m*$m*$m*$m]; # m^4
set g [expr 9.81*$m/pow($sec,2)]; # gravitational acceleration
proc BuildRCrectSection {id HSec BSec coverH coverB coreID coverID steelID numBarsTop barAreaTop numBarsBot barAreaBot numBarsIntTot barAreaInt nfCoreY nfCoreZ nfCoverY nfCoverZ} {
################################################
# BuildRCrectSection $id $HSec $BSec $coverH $coverB $coreID $coverID $steelID $numBarsTop $barAreaTop $numBarsBot $barAreaBot $numBarsIntTot $barAreaInt $nfCoreY $nfCoreZ $nfCoverY $nfCoverZ
################################################
set coverY [expr $HSec/2.0]; # The distance from the section z-axis to the edge of the cover concrete -- outer edge of cover concrete
set coverZ [expr $BSec/2.0]; # The distance from the section y-axis to the edge of the cover concrete -- outer edge of cover concrete
set coreY [expr $coverY-$coverH]; # The distance from the section z-axis to the edge of the core concrete -- edge of the core concrete/inner edge of cover concrete
set coreZ [expr $coverZ-$coverB]; # The distance from the section y-axis to the edge of the core concrete -- edge of the core concrete/inner edge of cover concrete
set numBarsInt [expr $numBarsIntTot/2]; # number of intermediate bars per side
# Define the fiber section
section fiberSec $id {
# Define the core patch
patch quadr $coreID $nfCoreZ $nfCoreY -$coreY $coreZ -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ
# Define the four cover patches
patch quadr $coverID 2 $nfCoverY -$coverY $coverZ -$coreY $coreZ $coreY $coreZ $coverY $coverZ
patch quadr $coverID 2 $nfCoverY -$coreY -$coreZ -$coverY -$coverZ $coverY -$coverZ $coreY -$coreZ
patch quadr $coverID $nfCoverZ 2 -$coverY $coverZ -$coverY -$coverZ -$coreY -$coreZ -$coreY $coreZ
patch quadr $coverID $nfCoverZ 2 $coreY $coreZ $coreY -$coreZ $coverY -$coverZ $coverY $coverZ
# define reinforcing layers
layer straight $steelID $numBarsInt $barAreaInt -$coreY $coreZ $coreY $coreZ; # intermediate skin reinf. +z
layer straight $steelID $numBarsInt $barAreaInt -$coreY -$coreZ $coreY -$coreZ; # intermediate skin reinf. -z
layer straight $steelID $numBarsTop $barAreaTop $coreY $coreZ $coreY -$coreZ; # top layer reinfocement
layer straight $steelID $numBarsBot $barAreaBot -$coreY $coreZ -$coreY -$coreZ; # bottom layer reinforcement
}; # end of fibersection definition
}; # end of procedure
# define GEOMETRY -------------------------------------------------------------
# define structure-geometry paramters
set LCol [expr 3*$m]; # column height
set LBeam [expr 3*$m]; # beam length
# calculate locations of beam/column intersections:
set X1 0.;
set X2 [expr $X1 + $LBeam];
set X3 [expr $X2 + $LBeam];
set X4 [expr $X3 + $LBeam];
set Y1 0.;
set Y2 [expr $Y1 + $LCol];
set Y3 [expr $Y2 + $LCol];
set Y4 [expr $Y3 + $LCol];
# define nodal coordinates
node 11 $X1 $Y1
node 12 $X2 $Y1
node 13 $X3 $Y1
node 14 $X4 $Y1
node 21 $X1 $Y2
node 22 $X2 $Y2
node 23 $X3 $Y2
node 24 $X4 $Y2
node 31 $X1 $Y3
node 32 $X2 $Y3
node 33 $X3 $Y3
node 34 $X4 $Y3
node 41 $X1 $Y4
node 42 $X2 $Y4
node 43 $X3 $Y4
node 44 $X4 $Y4
# Set up parameters that are particular to the model for displacement control
set IDctrlNode 41; # node where displacement is read for displacement control
set IDctrlDOF 1; # degree of freedom of displacement read for displacement control
set NStory 3; # number of stories above ground level
set NBay 3; # number of bays
set LBuilding $Y4; # total building height
# BOUNDARY CONDITIONS
fix 11 1 1 1
fix 12 1 1 1
fix 13 1 1 1
fix 14 1 1 1
# Define SECTIONS -------------------------------------------------------------
set SectionType FiberSection; # options: Elastic FiberSection
# define section tags:
set ColSecTag 1
set BeamSecTag 2
# Section Properties:
set HCol [expr 0.5*$m]; # square-Column width
set BCol $HCol
set HBeam [expr 0.5*$m]; # Beam depth -- perpendicular to bending axis
set BBeam [expr 0.5*$m]; # Beam width -- parallel to bending axis
if {$SectionType == "Elastic"} {
# material properties:
set fc 4000*$psi; # concrete nominal compressive strength
set Ec [expr 57*$ksi*pow($fc/$psi,0.5)]; # concrete Young's Modulus
# column section properties:
set AgCol [expr $HCol*$BCol]; # rectuangular-Column cross-sectional area
set IzCol [expr 0.5*1./12*$BCol*pow($HCol,3)]; # about-local-z Rect-Column gross moment of inertial
# beam sections:
set AgBeam [expr $HBeam*$BBeam]; # rectuangular-Beam cross-sectional area
set IzBeam [expr 0.5*1./12*$BBeam*pow($HBeam,3)]; # about-local-z Rect-Beam cracked moment of inertial
section Elastic $ColSecTag $Ec $AgCol $IzCol
section Elastic $BeamSecTag $Ec $AgBeam $IzBeam
} elseif {$SectionType == "FiberSection"} {
# MaterialeRC.tcl: define a library of Reinforced-Concrete materials
##########################################################
#
# General Material parameters
set G $Ubig; # make stiff shear modulus
set J 1.0; # torsional section stiffness (G makes GJ large)
set GJ [expr $G*$J];
# -----------------------------------------------------------------------------------------------------# confined and unconfined CONCRETE
# nominal concrete compressive strength
set fck 25;
set fcm [expr $fck+8];
set fc [expr -$fck*1000*$kPa]; # CONCRETE Compressive Strength, kPa (+Tension, -Compression)
set Ec [expr 30500000*$kPa]; # Concrete Elastic Modulus
set nu 0.2;
set Gc [expr $Ec/2./[expr 1+$nu]]; # Torsional stiffness Modulus
# confined concrete
set Kfc 1.3; # ratio of confined to unconfined concrete strength
set Kres 0.2; # ratio of residual/ultimate to maximum stress
set fc1C [expr $Kfc*$fc]; # CONFINED concrete (mander model), maximum stress
set eps1C [expr 2.*$fc1C/$Ec]; # strain at maximum stress
set fc2C [expr $Kres*$fc1C]; # ultimate stress
set eps2C [expr 20*$eps1C]; # strain at ultimate stress
set lambda 0.1; # ratio between unloading slope at $eps2 and initial slope $Ec
# unconfined concrete
set fc1U $fc; # UNCONFINED concrete (todeschini parabolic model), maximum stress
set eps1U -0.003; # strain at maximum strength of unconfined concrete
set fc2U [expr $Kres*$fc1U]; # ultimate stress
set eps2U -0.01; # strain at ultimate stress
# tensile-strength properties
set ftC [expr -0.14*$fc1C]; # tensile strength +tension
set ftU [expr -0.14*$fc1U]; # tensile strength +tension
set Ets [expr $ftU/0.002]; # tension softening stiffness
# set up library of materials
if { [info exists imat ] != 1} {set imat 0}; # set value only if it has not been defined previously.
set IDconcCore 1
set IDconcCover 2
uniaxialMaterial Concrete02 $IDconcCore $fc1C $eps1C $fc2C $eps2C $lambda $ftC $Ets; # Core concrete (confined)
uniaxialMaterial Concrete02 $IDconcCover $fc1U $eps1U $fc2U $eps2U $lambda $ftU $Ets; # Cover concrete (unconfined)
# -----------------------------------------------------------------------------------------------------# REINFORCING STEEL parameters
#
set Fy [expr 450000*$kPa]; # STEEL yield stress
set Es [expr 200000000.*$kPa]; # modulus of steel
set Bs 0.01; # strain-hardening ratio
set R0 18; # control the transition from elastic to plastic branches
set cR1 0.925; # control the transition from elastic to plastic branches
set cR2 0.15; # control the transition from elastic to plastic branches
set IDSteel 3
uniaxialMaterial Steel02 $IDSteel $Fy $Es $Bs $R0 $cR1 $cR2
# FIBER SECTION properties
# Column section geometry:
set cover [expr 0.05*$m]; # rectangular-RC-Column cover
set numBarsTopCol 6; # number of longitudinal-reinforcement bars on top layer
set numBarsBotCol 6; # number of longitudinal-reinforcement bars on bottom layer
set numBarsIntCol 3; # TOTAL number of reinforcing bars on the intermediate layers
set barAreaTopCol [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set barAreaBotCol [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set barAreaIntCol [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set numBarsTopBeam 6; # number of longitudinal-reinforcement bars on top layer
set numBarsBotBeam 6; # number of longitudinal-reinforcement bars on bottom layer
set numBarsIntBeam 3; # TOTAL number of reinforcing bars on the intermediate layers
set barAreaTopBeam [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set barAreaBotBeam [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set barAreaIntBeam [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set nfCoreY 20; # number of fibers in the core patch in the y direction
set nfCoreZ 20; # number of fibers in the core patch in the z direction
set nfCoverY 20; # number of fibers in the cover patches with long sides in the y direction
set nfCoverZ 20; # number of fibers in the cover patches with long sides in the z direction
# rectangular section with one layer of steel evenly distributed around the perimeter and a confined core.
BuildRCrectSection $ColSecTag $HCol $BCol $cover $cover $IDconcCore $IDconcCover $IDSteel $numBarsTopCol $barAreaTopCol $numBarsBotCol $barAreaBotCol $numBarsIntCol $barAreaIntCol $nfCoreY $nfCoreZ $nfCoverY $nfCoverZ
BuildRCrectSection $BeamSecTag $HBeam $BBeam $cover $cover $IDconcCore $IDconcCover $IDSteel $numBarsTopBeam $barAreaTopBeam $numBarsBotBeam $barAreaBotBeam $numBarsIntBeam $barAreaIntBeam $nfCoreY $nfCoreZ $nfCoverY $nfCoverZ
} else {
puts "No section has been defined"
return -1
}
# define ELEMENTS
# set up geometric transformations of element
# separate columns and beams, in case of P-Delta analysis for columns
set IDColTransf 1; # all columns
set IDBeamTransf 2; # all beams
set ColTransfType Linear ; # options, Linear PDelta Corotational
geomTransf $ColTransfType $IDColTransf ; # only columns can have PDelta effects (gravity effects)
geomTransf Linear $IDBeamTransf
# Define Beam-Column Elements
set np 5; # number of Gauss integration points for nonlinear curvature distribution-- np=2 for linear distribution ok
# columns
element nonlinearBeamColumn 111 11 21 $np $ColSecTag $IDColTransf; # level 1-2
element nonlinearBeamColumn 112 12 22 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 113 13 23 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 114 14 24 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 121 21 31 $np $ColSecTag $IDColTransf; # level 2-3
element nonlinearBeamColumn 122 22 32 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 123 23 33 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 124 24 34 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 131 31 41 $np $ColSecTag $IDColTransf; # level 3-4
element nonlinearBeamColumn 132 32 42 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 133 33 43 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 134 34 44 $np $ColSecTag $IDColTransf
# beams
element nonlinearBeamColumn 221 21 22 $np $BeamSecTag $IDBeamTransf; # level 2
element nonlinearBeamColumn 222 22 23 $np $BeamSecTag $IDBeamTransf;
element nonlinearBeamColumn 223 23 24 $np $BeamSecTag $IDBeamTransf;
element nonlinearBeamColumn 231 31 32 $np $BeamSecTag $IDBeamTransf; # level 3
element nonlinearBeamColumn 232 32 33 $np $BeamSecTag $IDBeamTransf;
element nonlinearBeamColumn 233 33 34 $np $BeamSecTag $IDBeamTransf;
element nonlinearBeamColumn 241 41 42 $np $BeamSecTag $IDBeamTransf; # level 4
element nonlinearBeamColumn 242 42 43 $np $BeamSecTag $IDBeamTransf;
element nonlinearBeamColumn 243 43 44 $np $BeamSecTag $IDBeamTransf;
# Define GRAVITY LOADS, weight and masses
# calculate dead load of frame, assume this to be an internal frame (do LL in a similar manner)
# calculate distributed weight along the beam length
set GammaConcrete [expr 25*$ps]; # Reinforced-Concrete floor slabs
set Tslab [expr 0.2*$m]; # 6-inch slab
set Lslab [expr 2*$LBeam/2]; # assume slab extends a distance of $LBeam1/2 in/out of plane
set Qslab [expr $GammaConcrete*$Tslab*$Lslab];
set QBeam [expr $GammaConcrete*$BBeam*$HBeam]; # W-section weight per length
set QdlBeam [expr $Qslab + $QBeam]; # dead load distributed along beam.
set QdlCol [expr $GammaConcrete*$BCol*$HCol]; # W-section weight per length
set WeightCol [expr $QdlCol*$LCol]; # total Column weight
set WeightBeam [expr $QdlBeam*$LBeam]; # total Beam weight
# assign masses to the nodes that the columns are connected to
# each connection takes the mass of 1/2 of each element framing into it (mass=weight/$g)
mass 21 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2)/$g] 0. 0.; # level 2
mass 22 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 23 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 24 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2)/$g] 0. 0.;
mass 31 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2)/$g] 0. 0.; # level 3
mass 32 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 33 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 34 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2)/$g] 0. 0.;
mass 41 [expr ($WeightCol/2 +$WeightBeam/2)/$g] 0. 0.; # level 4
mass 42 [expr ($WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 43 [expr ($WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 44 [expr ($WeightCol/2 +$WeightBeam/2)/$g] 0. 0.;
# calculate total Floor Mass
set WeightFloor2 [expr $WeightCol*4/2+$WeightCol*4/2+3*$WeightBeam]; # level 2 weight
set WeightFloor3 [expr $WeightCol*4/2+$WeightCol*4/2+3*$WeightBeam];
set WeightFloor4 [expr $WeightCol*4/2+3*$WeightBeam];
set WeightTotal [expr $WeightFloor2 + $WeightFloor3 + $WeightFloor4]; # total frame weight
set MassFloor2 [expr $WeightFloor2/$g];
set MassFloor3 [expr $WeightFloor3/$g];
set MassFloor4 [expr $WeightFloor4/$g];
set MassTotal [expr $MassFloor2+$MassFloor3+$MassFloor4]; # total frame mass
# LATERAL-LOAD distribution for static pushover analysis
# calculate distribution of lateral load based on mass/weight distributions along building height
# Fj = WjHj/sum(WiHi) * Weight at each floor j
set sumWiHi [expr $WeightFloor2*$Y2 + $WeightFloor3*$Y3 + $WeightFloor4*$Y4]; # denominator
set Fj2 [expr $WeightFloor2*$Y2/$sumWiHi*$WeightTotal]; # total for floor 2
set Fj3 [expr $WeightFloor3*$Y3/$sumWiHi*$WeightTotal]; # total for floor 3
set Fj4 [expr $WeightFloor4*$Y4/$sumWiHi*$WeightTotal]; # total for floor 4
set Fi2 [expr $Fj2/4]; # per node on floor 2
set Fi3 [expr $Fj3/4]; # per node on floor 3
set Fi4 [expr $Fj4/4]; # per node on floor 4
set iFi "$Fi2 $Fi3 $Fi4"; # vectorize
# Define RECORDERS -------------------------------------------------------------
recorder Node -file $dataDir/DFreelevel1.out -time -node 21 22 23 24 -dof 1 disp; # displacements of free node
recorder Node -file $dataDir/DFreelevel2.out -time -node 31 32 33 34 -dof 1 disp; # displacements of free node
recorder Node -file $dataDir/DFreelevel3.out -time -node 41 42 43 44 -dof 1 disp; # displacements of free node
recorder Node -file $dataDir/DBase.out -time -node 11 12 13 14 -dof 1 2 3 disp; # displacements of support nodes
recorder Node -file $dataDir/RBase.out -time -node 11 12 13 14 -dof 1 2 3 reaction; # support reaction
# define GRAVITY -------------------------------------------------------------
# GRAVITY LOADS # define gravity load applied to beams and columns -- eleLoad applies loads in local coordinate axis
pattern Plain 101 Linear {
eleLoad -ele 221 222 223 -type -beamUniform -$QdlBeam; ; # beams level 2 (in -ydirection)
eleLoad -ele 231 232 233 -type -beamUniform -$QdlBeam;
eleLoad -ele 241 242 243 -type -beamUniform -$QdlBeam
eleLoad -ele 111 112 113 114 -type -beamUniform 0 -$QdlCol; # columns level 1-2 (in -xdirection)
eleLoad -ele 121 122 123 124 -type -beamUniform 0 -$QdlCol;
eleLoad -ele 131 132 133 134 -type -beamUniform 0 -$QdlCol;
}
# Gravity-analysis parameters -- load-controlled static analysis
set Tol 1.0e-8; # 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 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
#eigen generalized symmSparseArpack 3 ;
analyze $NstepGravity; # apply gravity
# ------------------------------------------------- maintain constant gravity loads and reset time to zero
loadConst -time 0.0
puts "Model Built"
# DYNAMIC EQ ANALYSIS --------------------------------------------------------
# Uniform Earthquake ground motion (uniform acceleration input at all support nodes)
set GMdirection 1; # ground-motion direction
set GMfile "EQ7.acc" ; # ground-motion filenames
set GMfact 3.0;
# ground-motion scaling factor
# set up ground-motion-analysis parameters
set DtAnalysis [expr 0.01]; # time-step Dt for lateral analysis
set TmaxAnalysis [expr 20.00]; # maximum duration of ground-motion analysis -- should be 50*$sec
# DYNAMIC ANALYSIS PARAMETERS
constraints Transformation ;
numberer Plain
system SparseGeneral -piv
set Tol 1.e-8; # 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 EnergyIncr; # Convergence-test type
test $TestType $Tol $maxNumIter $printFlag;
set algorithmType ModifiedNewton
algorithm $algorithmType;
set NewmarkGamma 0.5; # Newmark-integrator gamma parameter (also HHT)
set NewmarkBeta 0.25; # Newmark-integrator beta parameter
integrator Newmark $NewmarkGamma $NewmarkBeta
analysis Transient
# Uniform EXCITATION: acceleration input
set IDloadTag 400; # load tag
set dt 0.01; # time step for input ground motion
set GMfatt [expr $GMfact*9.81]; # data in input file is in g Unifts -- ACCELERATION TH
set AccelSeries "Series -dt $dt -filePath $GMfile -factor $GMfatt"; # time series information
pattern UniformExcitation $IDloadTag $GMdirection -accel $AccelSeries ; # create Unifform excitation
set Nsteps [expr int($TmaxAnalysis/$DtAnalysis)];
set ok [analyze $Nsteps $DtAnalysis]; # actually perform analysis; returns ok=0 if analysis was successful
if {$ok != 0} { ; # if analysis was not successful.
# change some analysis parameters to achieve convergence
# performance is slower inside this loop
# Time-controlled analysis
set ok 0;
set controlTime [getTime];
while {$controlTime < $TmaxAnalysis && $ok == 0} {
set ok [analyze 1 $DtAnalysis]
set controlTime [getTime]
set ok [analyze 1 $DtAnalysis]
if {$ok != 0} {
puts "Trying Newton with Initial Tangent .."
test NormDispIncr $Tol 1000 0
algorithm Newton -initial
set ok [analyze 1 $DtAnalysis]
test $TestType $Tol $maxNumIter 0
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying Broyden .."
algorithm Broyden 8
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch .8
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmType
}
}
}; # end if ok !0
puts "Ground Motion Done. End Time: [getTime]"
Gentilmemte if I send you my model
give me a hand to solve the problem?
wipe; # clear memory of all past model definitions
model BasicBuilder -ndm 2 -ndf 3; # Define the model builder, ndm=#dimension, ndf=#dofs
set dataDir TelaioFiberRCSectionEQ; # set up name of data directory (you can remove this)
file mkdir $dataDir; # create data directory
set GMdir "../GMfiles/"; # ground-motion file directory
# define UNITS ----------------------------------------------------------------------------
set m 1.; # define basic units -- output units
set kN 1.; # define basic units -- output units
set sec 1.; # define basic units -- output units
set LunitTXT "m"; # define basic-unit text for output
set FunitTXT "kN"; # define basic-unit text for output
set TunitTXT "sec"; # define basic-unit text for output
set kPa [expr $kN/pow($m,2)];
set ps [expr $kN/pow($m,3)]; # peso specifico
set m2 [expr $m*$m]; # m^2
set m4 [expr $m*$m*$m*$m]; # m^4
set g [expr 9.81*$m/pow($sec,2)]; # gravitational acceleration
proc BuildRCrectSection {id HSec BSec coverH coverB coreID coverID steelID numBarsTop barAreaTop numBarsBot barAreaBot numBarsIntTot barAreaInt nfCoreY nfCoreZ nfCoverY nfCoverZ} {
################################################
# BuildRCrectSection $id $HSec $BSec $coverH $coverB $coreID $coverID $steelID $numBarsTop $barAreaTop $numBarsBot $barAreaBot $numBarsIntTot $barAreaInt $nfCoreY $nfCoreZ $nfCoverY $nfCoverZ
################################################
set coverY [expr $HSec/2.0]; # The distance from the section z-axis to the edge of the cover concrete -- outer edge of cover concrete
set coverZ [expr $BSec/2.0]; # The distance from the section y-axis to the edge of the cover concrete -- outer edge of cover concrete
set coreY [expr $coverY-$coverH]; # The distance from the section z-axis to the edge of the core concrete -- edge of the core concrete/inner edge of cover concrete
set coreZ [expr $coverZ-$coverB]; # The distance from the section y-axis to the edge of the core concrete -- edge of the core concrete/inner edge of cover concrete
set numBarsInt [expr $numBarsIntTot/2]; # number of intermediate bars per side
# Define the fiber section
section fiberSec $id {
# Define the core patch
patch quadr $coreID $nfCoreZ $nfCoreY -$coreY $coreZ -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ
# Define the four cover patches
patch quadr $coverID 2 $nfCoverY -$coverY $coverZ -$coreY $coreZ $coreY $coreZ $coverY $coverZ
patch quadr $coverID 2 $nfCoverY -$coreY -$coreZ -$coverY -$coverZ $coverY -$coverZ $coreY -$coreZ
patch quadr $coverID $nfCoverZ 2 -$coverY $coverZ -$coverY -$coverZ -$coreY -$coreZ -$coreY $coreZ
patch quadr $coverID $nfCoverZ 2 $coreY $coreZ $coreY -$coreZ $coverY -$coverZ $coverY $coverZ
# define reinforcing layers
layer straight $steelID $numBarsInt $barAreaInt -$coreY $coreZ $coreY $coreZ; # intermediate skin reinf. +z
layer straight $steelID $numBarsInt $barAreaInt -$coreY -$coreZ $coreY -$coreZ; # intermediate skin reinf. -z
layer straight $steelID $numBarsTop $barAreaTop $coreY $coreZ $coreY -$coreZ; # top layer reinfocement
layer straight $steelID $numBarsBot $barAreaBot -$coreY $coreZ -$coreY -$coreZ; # bottom layer reinforcement
}; # end of fibersection definition
}; # end of procedure
# define GEOMETRY -------------------------------------------------------------
# define structure-geometry paramters
set LCol [expr 3*$m]; # column height
set LBeam [expr 3*$m]; # beam length
# calculate locations of beam/column intersections:
set X1 0.;
set X2 [expr $X1 + $LBeam];
set X3 [expr $X2 + $LBeam];
set X4 [expr $X3 + $LBeam];
set Y1 0.;
set Y2 [expr $Y1 + $LCol];
set Y3 [expr $Y2 + $LCol];
set Y4 [expr $Y3 + $LCol];
# define nodal coordinates
node 11 $X1 $Y1
node 12 $X2 $Y1
node 13 $X3 $Y1
node 14 $X4 $Y1
node 21 $X1 $Y2
node 22 $X2 $Y2
node 23 $X3 $Y2
node 24 $X4 $Y2
node 31 $X1 $Y3
node 32 $X2 $Y3
node 33 $X3 $Y3
node 34 $X4 $Y3
node 41 $X1 $Y4
node 42 $X2 $Y4
node 43 $X3 $Y4
node 44 $X4 $Y4
# Set up parameters that are particular to the model for displacement control
set IDctrlNode 41; # node where displacement is read for displacement control
set IDctrlDOF 1; # degree of freedom of displacement read for displacement control
set NStory 3; # number of stories above ground level
set NBay 3; # number of bays
set LBuilding $Y4; # total building height
# BOUNDARY CONDITIONS
fix 11 1 1 1
fix 12 1 1 1
fix 13 1 1 1
fix 14 1 1 1
# Define SECTIONS -------------------------------------------------------------
set SectionType FiberSection; # options: Elastic FiberSection
# define section tags:
set ColSecTag 1
set BeamSecTag 2
# Section Properties:
set HCol [expr 0.5*$m]; # square-Column width
set BCol $HCol
set HBeam [expr 0.5*$m]; # Beam depth -- perpendicular to bending axis
set BBeam [expr 0.5*$m]; # Beam width -- parallel to bending axis
if {$SectionType == "Elastic"} {
# material properties:
set fc 4000*$psi; # concrete nominal compressive strength
set Ec [expr 57*$ksi*pow($fc/$psi,0.5)]; # concrete Young's Modulus
# column section properties:
set AgCol [expr $HCol*$BCol]; # rectuangular-Column cross-sectional area
set IzCol [expr 0.5*1./12*$BCol*pow($HCol,3)]; # about-local-z Rect-Column gross moment of inertial
# beam sections:
set AgBeam [expr $HBeam*$BBeam]; # rectuangular-Beam cross-sectional area
set IzBeam [expr 0.5*1./12*$BBeam*pow($HBeam,3)]; # about-local-z Rect-Beam cracked moment of inertial
section Elastic $ColSecTag $Ec $AgCol $IzCol
section Elastic $BeamSecTag $Ec $AgBeam $IzBeam
} elseif {$SectionType == "FiberSection"} {
# MaterialeRC.tcl: define a library of Reinforced-Concrete materials
##########################################################
#
# General Material parameters
set G $Ubig; # make stiff shear modulus
set J 1.0; # torsional section stiffness (G makes GJ large)
set GJ [expr $G*$J];
# -----------------------------------------------------------------------------------------------------# confined and unconfined CONCRETE
# nominal concrete compressive strength
set fck 25;
set fcm [expr $fck+8];
set fc [expr -$fck*1000*$kPa]; # CONCRETE Compressive Strength, kPa (+Tension, -Compression)
set Ec [expr 30500000*$kPa]; # Concrete Elastic Modulus
set nu 0.2;
set Gc [expr $Ec/2./[expr 1+$nu]]; # Torsional stiffness Modulus
# confined concrete
set Kfc 1.3; # ratio of confined to unconfined concrete strength
set Kres 0.2; # ratio of residual/ultimate to maximum stress
set fc1C [expr $Kfc*$fc]; # CONFINED concrete (mander model), maximum stress
set eps1C [expr 2.*$fc1C/$Ec]; # strain at maximum stress
set fc2C [expr $Kres*$fc1C]; # ultimate stress
set eps2C [expr 20*$eps1C]; # strain at ultimate stress
set lambda 0.1; # ratio between unloading slope at $eps2 and initial slope $Ec
# unconfined concrete
set fc1U $fc; # UNCONFINED concrete (todeschini parabolic model), maximum stress
set eps1U -0.003; # strain at maximum strength of unconfined concrete
set fc2U [expr $Kres*$fc1U]; # ultimate stress
set eps2U -0.01; # strain at ultimate stress
# tensile-strength properties
set ftC [expr -0.14*$fc1C]; # tensile strength +tension
set ftU [expr -0.14*$fc1U]; # tensile strength +tension
set Ets [expr $ftU/0.002]; # tension softening stiffness
# set up library of materials
if { [info exists imat ] != 1} {set imat 0}; # set value only if it has not been defined previously.
set IDconcCore 1
set IDconcCover 2
uniaxialMaterial Concrete02 $IDconcCore $fc1C $eps1C $fc2C $eps2C $lambda $ftC $Ets; # Core concrete (confined)
uniaxialMaterial Concrete02 $IDconcCover $fc1U $eps1U $fc2U $eps2U $lambda $ftU $Ets; # Cover concrete (unconfined)
# -----------------------------------------------------------------------------------------------------# REINFORCING STEEL parameters
#
set Fy [expr 450000*$kPa]; # STEEL yield stress
set Es [expr 200000000.*$kPa]; # modulus of steel
set Bs 0.01; # strain-hardening ratio
set R0 18; # control the transition from elastic to plastic branches
set cR1 0.925; # control the transition from elastic to plastic branches
set cR2 0.15; # control the transition from elastic to plastic branches
set IDSteel 3
uniaxialMaterial Steel02 $IDSteel $Fy $Es $Bs $R0 $cR1 $cR2
# FIBER SECTION properties
# Column section geometry:
set cover [expr 0.05*$m]; # rectangular-RC-Column cover
set numBarsTopCol 6; # number of longitudinal-reinforcement bars on top layer
set numBarsBotCol 6; # number of longitudinal-reinforcement bars on bottom layer
set numBarsIntCol 3; # TOTAL number of reinforcing bars on the intermediate layers
set barAreaTopCol [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set barAreaBotCol [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set barAreaIntCol [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set numBarsTopBeam 6; # number of longitudinal-reinforcement bars on top layer
set numBarsBotBeam 6; # number of longitudinal-reinforcement bars on bottom layer
set numBarsIntBeam 3; # TOTAL number of reinforcing bars on the intermediate layers
set barAreaTopBeam [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set barAreaBotBeam [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set barAreaIntBeam [expr 0.00020106*$m*$m]; # longitudinal-reinforcement bar area
set nfCoreY 20; # number of fibers in the core patch in the y direction
set nfCoreZ 20; # number of fibers in the core patch in the z direction
set nfCoverY 20; # number of fibers in the cover patches with long sides in the y direction
set nfCoverZ 20; # number of fibers in the cover patches with long sides in the z direction
# rectangular section with one layer of steel evenly distributed around the perimeter and a confined core.
BuildRCrectSection $ColSecTag $HCol $BCol $cover $cover $IDconcCore $IDconcCover $IDSteel $numBarsTopCol $barAreaTopCol $numBarsBotCol $barAreaBotCol $numBarsIntCol $barAreaIntCol $nfCoreY $nfCoreZ $nfCoverY $nfCoverZ
BuildRCrectSection $BeamSecTag $HBeam $BBeam $cover $cover $IDconcCore $IDconcCover $IDSteel $numBarsTopBeam $barAreaTopBeam $numBarsBotBeam $barAreaBotBeam $numBarsIntBeam $barAreaIntBeam $nfCoreY $nfCoreZ $nfCoverY $nfCoverZ
} else {
puts "No section has been defined"
return -1
}
# define ELEMENTS
# set up geometric transformations of element
# separate columns and beams, in case of P-Delta analysis for columns
set IDColTransf 1; # all columns
set IDBeamTransf 2; # all beams
set ColTransfType Linear ; # options, Linear PDelta Corotational
geomTransf $ColTransfType $IDColTransf ; # only columns can have PDelta effects (gravity effects)
geomTransf Linear $IDBeamTransf
# Define Beam-Column Elements
set np 5; # number of Gauss integration points for nonlinear curvature distribution-- np=2 for linear distribution ok
# columns
element nonlinearBeamColumn 111 11 21 $np $ColSecTag $IDColTransf; # level 1-2
element nonlinearBeamColumn 112 12 22 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 113 13 23 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 114 14 24 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 121 21 31 $np $ColSecTag $IDColTransf; # level 2-3
element nonlinearBeamColumn 122 22 32 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 123 23 33 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 124 24 34 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 131 31 41 $np $ColSecTag $IDColTransf; # level 3-4
element nonlinearBeamColumn 132 32 42 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 133 33 43 $np $ColSecTag $IDColTransf
element nonlinearBeamColumn 134 34 44 $np $ColSecTag $IDColTransf
# beams
element nonlinearBeamColumn 221 21 22 $np $BeamSecTag $IDBeamTransf; # level 2
element nonlinearBeamColumn 222 22 23 $np $BeamSecTag $IDBeamTransf;
element nonlinearBeamColumn 223 23 24 $np $BeamSecTag $IDBeamTransf;
element nonlinearBeamColumn 231 31 32 $np $BeamSecTag $IDBeamTransf; # level 3
element nonlinearBeamColumn 232 32 33 $np $BeamSecTag $IDBeamTransf;
element nonlinearBeamColumn 233 33 34 $np $BeamSecTag $IDBeamTransf;
element nonlinearBeamColumn 241 41 42 $np $BeamSecTag $IDBeamTransf; # level 4
element nonlinearBeamColumn 242 42 43 $np $BeamSecTag $IDBeamTransf;
element nonlinearBeamColumn 243 43 44 $np $BeamSecTag $IDBeamTransf;
# Define GRAVITY LOADS, weight and masses
# calculate dead load of frame, assume this to be an internal frame (do LL in a similar manner)
# calculate distributed weight along the beam length
set GammaConcrete [expr 25*$ps]; # Reinforced-Concrete floor slabs
set Tslab [expr 0.2*$m]; # 6-inch slab
set Lslab [expr 2*$LBeam/2]; # assume slab extends a distance of $LBeam1/2 in/out of plane
set Qslab [expr $GammaConcrete*$Tslab*$Lslab];
set QBeam [expr $GammaConcrete*$BBeam*$HBeam]; # W-section weight per length
set QdlBeam [expr $Qslab + $QBeam]; # dead load distributed along beam.
set QdlCol [expr $GammaConcrete*$BCol*$HCol]; # W-section weight per length
set WeightCol [expr $QdlCol*$LCol]; # total Column weight
set WeightBeam [expr $QdlBeam*$LBeam]; # total Beam weight
# assign masses to the nodes that the columns are connected to
# each connection takes the mass of 1/2 of each element framing into it (mass=weight/$g)
mass 21 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2)/$g] 0. 0.; # level 2
mass 22 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 23 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 24 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2)/$g] 0. 0.;
mass 31 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2)/$g] 0. 0.; # level 3
mass 32 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 33 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 34 [expr ($WeightCol/2 + $WeightCol/2 +$WeightBeam/2)/$g] 0. 0.;
mass 41 [expr ($WeightCol/2 +$WeightBeam/2)/$g] 0. 0.; # level 4
mass 42 [expr ($WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 43 [expr ($WeightCol/2 +$WeightBeam/2 +$WeightBeam/2)/$g] 0. 0.;
mass 44 [expr ($WeightCol/2 +$WeightBeam/2)/$g] 0. 0.;
# calculate total Floor Mass
set WeightFloor2 [expr $WeightCol*4/2+$WeightCol*4/2+3*$WeightBeam]; # level 2 weight
set WeightFloor3 [expr $WeightCol*4/2+$WeightCol*4/2+3*$WeightBeam];
set WeightFloor4 [expr $WeightCol*4/2+3*$WeightBeam];
set WeightTotal [expr $WeightFloor2 + $WeightFloor3 + $WeightFloor4]; # total frame weight
set MassFloor2 [expr $WeightFloor2/$g];
set MassFloor3 [expr $WeightFloor3/$g];
set MassFloor4 [expr $WeightFloor4/$g];
set MassTotal [expr $MassFloor2+$MassFloor3+$MassFloor4]; # total frame mass
# LATERAL-LOAD distribution for static pushover analysis
# calculate distribution of lateral load based on mass/weight distributions along building height
# Fj = WjHj/sum(WiHi) * Weight at each floor j
set sumWiHi [expr $WeightFloor2*$Y2 + $WeightFloor3*$Y3 + $WeightFloor4*$Y4]; # denominator
set Fj2 [expr $WeightFloor2*$Y2/$sumWiHi*$WeightTotal]; # total for floor 2
set Fj3 [expr $WeightFloor3*$Y3/$sumWiHi*$WeightTotal]; # total for floor 3
set Fj4 [expr $WeightFloor4*$Y4/$sumWiHi*$WeightTotal]; # total for floor 4
set Fi2 [expr $Fj2/4]; # per node on floor 2
set Fi3 [expr $Fj3/4]; # per node on floor 3
set Fi4 [expr $Fj4/4]; # per node on floor 4
set iFi "$Fi2 $Fi3 $Fi4"; # vectorize
# Define RECORDERS -------------------------------------------------------------
recorder Node -file $dataDir/DFreelevel1.out -time -node 21 22 23 24 -dof 1 disp; # displacements of free node
recorder Node -file $dataDir/DFreelevel2.out -time -node 31 32 33 34 -dof 1 disp; # displacements of free node
recorder Node -file $dataDir/DFreelevel3.out -time -node 41 42 43 44 -dof 1 disp; # displacements of free node
recorder Node -file $dataDir/DBase.out -time -node 11 12 13 14 -dof 1 2 3 disp; # displacements of support nodes
recorder Node -file $dataDir/RBase.out -time -node 11 12 13 14 -dof 1 2 3 reaction; # support reaction
# define GRAVITY -------------------------------------------------------------
# GRAVITY LOADS # define gravity load applied to beams and columns -- eleLoad applies loads in local coordinate axis
pattern Plain 101 Linear {
eleLoad -ele 221 222 223 -type -beamUniform -$QdlBeam; ; # beams level 2 (in -ydirection)
eleLoad -ele 231 232 233 -type -beamUniform -$QdlBeam;
eleLoad -ele 241 242 243 -type -beamUniform -$QdlBeam
eleLoad -ele 111 112 113 114 -type -beamUniform 0 -$QdlCol; # columns level 1-2 (in -xdirection)
eleLoad -ele 121 122 123 124 -type -beamUniform 0 -$QdlCol;
eleLoad -ele 131 132 133 134 -type -beamUniform 0 -$QdlCol;
}
# Gravity-analysis parameters -- load-controlled static analysis
set Tol 1.0e-8; # 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 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
#eigen generalized symmSparseArpack 3 ;
analyze $NstepGravity; # apply gravity
# ------------------------------------------------- maintain constant gravity loads and reset time to zero
loadConst -time 0.0
puts "Model Built"
# DYNAMIC EQ ANALYSIS --------------------------------------------------------
# Uniform Earthquake ground motion (uniform acceleration input at all support nodes)
set GMdirection 1; # ground-motion direction
set GMfile "EQ7.acc" ; # ground-motion filenames
set GMfact 3.0;
# ground-motion scaling factor
# set up ground-motion-analysis parameters
set DtAnalysis [expr 0.01]; # time-step Dt for lateral analysis
set TmaxAnalysis [expr 20.00]; # maximum duration of ground-motion analysis -- should be 50*$sec
# DYNAMIC ANALYSIS PARAMETERS
constraints Transformation ;
numberer Plain
system SparseGeneral -piv
set Tol 1.e-8; # 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 EnergyIncr; # Convergence-test type
test $TestType $Tol $maxNumIter $printFlag;
set algorithmType ModifiedNewton
algorithm $algorithmType;
set NewmarkGamma 0.5; # Newmark-integrator gamma parameter (also HHT)
set NewmarkBeta 0.25; # Newmark-integrator beta parameter
integrator Newmark $NewmarkGamma $NewmarkBeta
analysis Transient
# Uniform EXCITATION: acceleration input
set IDloadTag 400; # load tag
set dt 0.01; # time step for input ground motion
set GMfatt [expr $GMfact*9.81]; # data in input file is in g Unifts -- ACCELERATION TH
set AccelSeries "Series -dt $dt -filePath $GMfile -factor $GMfatt"; # time series information
pattern UniformExcitation $IDloadTag $GMdirection -accel $AccelSeries ; # create Unifform excitation
set Nsteps [expr int($TmaxAnalysis/$DtAnalysis)];
set ok [analyze $Nsteps $DtAnalysis]; # actually perform analysis; returns ok=0 if analysis was successful
if {$ok != 0} { ; # if analysis was not successful.
# change some analysis parameters to achieve convergence
# performance is slower inside this loop
# Time-controlled analysis
set ok 0;
set controlTime [getTime];
while {$controlTime < $TmaxAnalysis && $ok == 0} {
set ok [analyze 1 $DtAnalysis]
set controlTime [getTime]
set ok [analyze 1 $DtAnalysis]
if {$ok != 0} {
puts "Trying Newton with Initial Tangent .."
test NormDispIncr $Tol 1000 0
algorithm Newton -initial
set ok [analyze 1 $DtAnalysis]
test $TestType $Tol $maxNumIter 0
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying Broyden .."
algorithm Broyden 8
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch .8
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmType
}
}
}; # end if ok !0
puts "Ground Motion Done. End Time: [getTime]"
-
GALLO ANTONIO
- Posts: 23
- Joined: Tue Nov 16, 2010 6:14 am
- Location: Sapienza
Re: CTestEnergyIncr::test() - failed to converge
Where the time history is as follows:
0.004389876
0.004138238
0.003538426
0.002956079
0.002300461
0.000958434
0.001340644
0.004399356
0.006974449
0.008374172
0.00788987
0.010585962
0.011566041
0.009968489
0.010941475
0.007110497
0.006148309
0.006578479
0.008128808
0.013749401
0.019004095
0.014562379
0.0069579
0.005822071
0.008795188
0.005871725
0.002216147
0.006665813
0.009231095
0.010979309
0.011716109
0.011303527
0.010068906
0.007797371
0.003007371
0.002053698
0.010031558
0.026374172
0.046171248
0.052403793
0.052932024
0.058454975
0.058768585
0.067543961
0.066608369
0.053094506
0.037442166
0.032747816
0.029241055
0.022598119
0.022125896
0.014702437
-0.002984791
-0.004900933
-0.013001241
-0.023455277
-0.013545915
-0.002640109
0.004623368
0.003935909
-0.009988242
-0.026583966
-0.044086028
-0.043933187
-0.025314456
-0.023492301
-0.040920228
-0.043290861
-0.040735353
-0.047085773
-0.031774867
-0.003795079
0.004417128
0.017984247
0.036742911
0.037315238
0.022797396
0.016780613
0.012980396
0.018816974
0.031083621
0.021892844
0.030406939
0.003875493
-0.049668845
-0.096326455
-0.099464335
-0.075044557
-0.068058118
-0.070133775
-0.060190443
-0.06692908
-0.050251424
-0.022884341
0.030718228
0.047963366
0.028263319
0.019636098
-0.009945349
-0.006615837
0.042356066
0.021618096
-0.025133919
-0.039636467
-0.036101315
-0.063346729
-0.109972924
-0.100563735
-0.052185543
-0.025679411
-0.01020145
0.019292211
0.055217765
0.075294882
0.044998959
0.042072207
0.034597144
0.034453262
0.043323968
0.060132377
0.035380051
0.046251744
0.038187187
0.011047371
0.009233863
0.034305055
0.021512978
0.024317242
0.025100769
0.055347566
0.082517222
0.099676646
0.11055956
0.095017403
0.091345496
0.08531516
0.091366097
0.058829732
-0.02147487
-0.062186196
-0.07657595
-0.07898993
-0.047280844
-0.023482857
-0.018443087
-0.005183474
0.00067654
-0.002617953
-0.044077925
-0.114724182
-0.141991869
-0.159872487
-0.217377454
-0.224519059
-0.178097233
-0.175022557
-0.176696211
-0.136746496
-0.123677075
-0.085693039
-0.013966632
0.035090212
0.050612431
0.045808882
0.072029606
0.064596757
0.035395753
0.10086754
0.111724965
0.039512467
-0.004864808
0.034604825
0.057628743
0.085553154
0.099744268
0.100784816
0.086710095
0.031881507
-0.020293508
-0.055310227
-0.023913372
-0.033108454
-0.054750212
-0.038600005
-0.019289847
0.025195153
0.031383514
0.031031655
0.00732173
-0.036314026
-0.015704393
0.04125287
0.006782964
-0.004442574
0.032867163
0.086281635
0.206558928
0.238221645
0.206363082
0.054481227
-0.034699071
-0.077240109
-0.068054847
-0.011507249
0.027141836
0.068018384
0.042309061
0.018113339
0.035914503
0.042923145
-0.005328159
0.003197319
-0.010312952
-0.03201396
0.039041474
0.085358322
0.069810122
0.019421455
-0.008638076
-0.015821826
-0.026125092
-0.09495388
-0.134317532
-0.104278386
-0.043131333
-0.041233797
-0.025282316
-0.024528721
-0.015251717
0.005746465
0.041982014
0.036497783
0.012848216
-0.011986255
-0.053300429
-0.05255511
-0.072286196
-0.067120001
-0.08280351
-0.095052712
-0.074766047
-0.068195023
-0.022107238
0.051868174
0.042095818
0.045275349
0.087857679
0.053255603
0.007478693
-0.006987691
-0.032102179
-0.054716956
-0.049552917
-0.014691803
0.039705373
0.105895624
0.116678119
0.049160235
0.063530371
0.131679595
0.164959475
0.185520545
0.175114229
0.18942976
0.187566891
0.183584258
0.129206046
0.086754017
0.094620831
0.067868114
0.091821954
0.154070199
0.169293538
0.17330642
0.160751805
0.116935693
0.074370682
0.068655297
0.099135317
0.062715508
0.025828622
-0.001267581
-0.031657115
-0.082366355
-0.075164802
-0.069023564
-0.077361226
-0.18009603
-0.238762692
-0.220637485
-0.179838777
-0.228515774
-0.228194639
-0.138895318
-0.140812069
-0.12339206
-0.064908393
-0.078718267
-0.062946007
-0.00103504
0.007666875
0.006007117
-0.009354099
-0.049894448
-0.016770985
-0.012860426
-0.079803064
-0.141280085
-0.110163137
-0.108224563
-0.111272588
-0.040706016
0.002670593
0.035435319
0.062670484
0.060045082
0.078587346
0.047201227
-0.012014098
-0.064838246
-0.089429051
-0.105027899
-0.113477737
-0.127376944
-0.121771365
-0.118962914
-0.160575449
-0.204382554
-0.253239125
-0.268444806
-0.249026641
-0.195515603
-0.161946267
-0.143055901
-0.087561198
-0.128968492
-0.154228136
-0.089996606
-0.08030273
-0.07514099
-0.069628678
-0.088012442
-0.1297618
-0.103878252
-0.019945446
0.038701799
0.029131671
0.038021017
0.102432355
0.128920987
0.069437489
0.010602042
-0.000483863
0.024478998
0.038328841
0.040727396
0.032824844
-0.015539099
-0.0503456
-0.014558623
-0.00748615
-0.034387209
-0.010439106
-0.013288285
-0.030585457
0.036311362
0.030015795
0.03252903
0.025640277
-0.008299512
0.011420015
0.036547035
0.007053625
-0.012425179
-0.042136468
-0.021312159
0.007412418
-0.050826088
-0.100183763
-0.115574025
-0.105981104
-0.08848273
-0.10301844
-0.125949264
-0.131981537
-0.087121747
-0.062966332
-0.044291142
-0.065409757
-0.016360315
0.055094779
0.099504277
0.169816658
0.148446813
0.053961568
-0.023430729
-0.071043953
-0.027465085
0.013452497
-0.010729436
-0.005107792
-0.001314264
-0.019754402
0.010598665
-0.001760871
0.053848825
0.047590174
0.024093928
0.023419751
0.020481698
0.066653065
0.114673659
0.071195804
0.104845189
0.154766858
0.146247268
0.198285833
0.268434793
0.276665241
0.220919564
0.150192663
0.070541821
0.022955595
0.037346646
0.052325886
0.02868638
-0.004585291
-0.040150806
-0.035184436
0.015591319
0.059459414
0.017683357
-0.059171733
-0.056212962
-0.058371
-0.09174601
-0.049149513
-0.077831045
-0.113683842
-0.075383417
-0.076165192
-0.132980496
-0.227275655
-0.184773013
-0.094276026
-0.070428468
-0.110784523
-0.137689441
-0.109232582
-0.079811096
-0.00835235
-0.025476513
-0.060736451
-0.027161319
-0.04330321
-0.019583208
0.006330533
0.013584677
0.008081361
0.017939458
0.054331191
0.040236071
0.02795046
0.059853215
0.084023401
0.160413429
0.190545768
0.203166395
0.185712293
0.175010711
0.133382872
0.11645788
0.113888092
0.065955505
0.115141883
0.194195032
0.207252294
0.162922472
0.162107989
0.129412577
0.037308112
0.001289477
-0.040498033
-0.055742372
-0.075794503
-0.094154865
-0.141408473
-0.164765954
-0.160227999
-0.186974078
-0.169996709
-0.135958284
-0.038093958
0.056442861
0.101021044
0.1292101
0.105391487
0.029447995
0.025793614
0.040120065
0.071422867
0.052875195
0.008043284
-0.008641883
-0.009572518
0.063596122
0.123738088
0.091491133
0.053688508
0.044477422
0.032545213
0.053915787
0.024620835
0.012838631
0.033009417
0.044526979
0.003437634
0.037034113
0.06787841
0.127609536
0.182717964
0.177482486
0.19061695
0.190846875
0.145910621
0.128924981
0.132716596
0.118908569
0.101858824
0.119268984
0.138784766
0.15850693
0.098713771
0.003999926
-0.055424936
-0.102571249
-0.148509979
-0.092618912
-0.01123844
0.043049876
0.050097577
0.039061014
0.051775746
0.039264865
-0.013775383
-0.038203597
-0.009675848
-0.046133373
0.009924787
0.053420857
-0.010197366
-0.058407374
-0.050644141
-0.04937207
-0.09228991
-0.121362038
-0.188260898
-0.17861031
-0.151900783
-0.106594056
-0.082105756
-0.068192683
-0.056719352
-0.072887599
-0.089877538
-0.139100552
-0.164130613
-0.135225073
-0.125551999
-0.147427708
-0.153776288
-0.082781613
-0.006815014
0.027425965
0.028930675
0.036038376
0.094926327
0.083093755
0.022690196
-0.02233386
-0.041778479
-0.008570619
0.05895713
0.114532627
0.087312065
0.090907335
0.10739664
0.039083377
0.040633582
0.047920797
0.059957188
0.064590298
0.051486667
0.062788874
0.083931699
0.093663625
0.100586951
0.048265848
-0.003041288
-0.0725343
-0.149264693
-0.164445162
-0.180927724
-0.165441051
-0.12300463
-0.072564319
0.006662579
-0.033231743
-0.042434316
0.018175816
0.054569922
0.047878623
-0.018910987
-0.01950888
-0.049443353
-0.142873123
-0.220553383
-0.272845298
-0.259322971
-0.235190585
-0.225628868
-0.216409728
-0.195022687
-0.13942261
-0.116300173
-0.142197013
-0.168836325
-0.161538079
-0.156031296
-0.115752056
-0.080743782
0.004815917
0.084683761
0.044227317
-0.030619025
-0.012500649
0.075133547
0.178499013
0.222744703
0.209623083
0.28842479
0.350096583
0.299970031
0.250269502
0.196404546
0.171851456
0.167426467
0.246651828
0.268478811
0.218463749
0.199080229
0.198228478
0.13845484
0.097431853
0.068056002
0.049522586
0.053581707
0.028624281
0.013652835
-0.025569383
-0.123640411
-0.166638419
-0.148583293
-0.066730812
-0.00221677
0.05905636
0.04027288
0.055814914
0.07306575
0.082259335
0.116849631
0.106711403
0.057184469
0.008970197
-0.014485115
-0.00316432
0.000828466
-0.017712196
-0.076118156
-0.08145947
-0.02739669
0.017521003
0.050027005
0.07927978
0.135163993
0.165154636
0.184301138
0.161105707
0.175331637
0.19975163
0.175987929
0.166578814
0.158831343
0.171066731
0.148071662
0.063943021
-0.002897062
0.004464262
0.029098753
0.062108021
0.086413361
0.108353168
0.139394611
0.179077923
0.198447466
0.173068792
0.101427764
0.002875132
-0.122601017
-0.155247912
-0.201861262
-0.347599924
-0.350593895
-0.333272249
-0.288043499
-0.24214229
-0.220248953
-0.151501939
-0.134753287
-0.168050721
-0.186690286
-0.140635937
-0.082204916
-0.050225515
-0.006857778
-0.001305159
-0.022105826
-0.049406536
-0.107557513
-0.132905975
-0.070751436
-0.047414422
-0.08698988
-0.054657914
-0.01220632
0.053120159
0.116007313
0.117257677
0.108856715
0.112051934
0.124573335
0.127473161
0.214074612
0.262830645
0.238227755
0.200920805
0.245996028
0.24187395
0.119590208
0.010628896
-0.038126141
-0.078206673
-0.150070593
-0.197204709
-0.219218016
-0.291837543
-0.323337704
-0.315931797
-0.249574289
-0.199295342
-0.189535275
-0.169362649
-0.104707181
-0.104726233
-0.116100237
-0.102004662
-0.052015252
0.004629215
0.003209735
0.048084844
0.101253636
0.084979415
0.054007139
0.01546963
0.024102589
0.033831432
0.042306863
0.084896609
0.08665695
0.048973523
-0.011994101
-0.013486257
0.016142447
-0.025257366
-0.077199802
-0.095984653
-0.113123663
-0.139941931
-0.14556393
-0.155254737
-0.090928674
-0.041798696
-0.04139502
-0.012623093
0.06074708
0.095009409
0.164344758
0.241050079
0.235544905
0.253403604
0.266828448
0.226403415
0.175562873
0.205246001
0.261426687
0.247781545
0.188658938
0.205553874
0.240348041
0.207255423
0.168956161
0.141112164
0.091098554
0.05015352
0.020443449
-0.027955413
-0.059405155
-0.025786044
0.018179359
0.022815011
-0.00149671
-0.06096765
-0.07277213
-0.023560327
0.002005536
0.045602906
0.056839135
0.042840108
0.044762373
0.048839685
0.047889937
-0.020861952
0.024696078
0.099641278
0.055086851
-0.020806495
-0.003074558
0.009929988
-0.019656239
-0.078323901
-0.075673029
-0.085612819
-0.087740131
-0.032917436
-0.009666321
-0.050624903
-0.106706128
-0.138963774
-0.178197011
-0.185912162
-0.1700259
-0.210941702
-0.237469852
-0.249348566
-0.226435363
-0.194474086
-0.135895357
-0.012829147
0.085702844
0.079072699
0.051635355
0.056965314
0.006308983
-0.032660849
-0.026627831
-0.099196374
-0.104629628
-0.011109727
0.021691293
0.004987638
-0.026379542
-0.056380115
-0.027281251
0.04713764
0.093570001
0.084913194
0.069676168
0.064123839
0.036445986
0.021342762
-0.0249002
-0.012444232
-0.009573135
0.043087974
0.04736663
-0.013752058
0.021025168
0.024697382
0.001431413
0.063563272
0.075505018
0.048893183
0.035958856
-0.035261717
-0.113736808
-0.15111433
-0.122489884
-0.152192876
-0.183135793
-0.161948785
-0.159887642
-0.077011727
-0.044141281
-0.047557507
0.036005914
0.145909786
0.116012923
0.088287607
0.103065729
0.02945335
-0.020842232
-0.086875081
-0.140062749
-0.135303915
-0.088159755
-0.088365674
-0.058485456
-0.08154016
-0.091838963
-0.031374894
-0.014583507
-0.006194479
-0.017967548
-0.070953757
-0.063599363
-0.021371901
0.007267625
0.029860955
0.003405939
0.058771208
0.153367743
0.126692235
0.211100966
0.186623365
0.143556833
0.086655468
0.043308221
0.039305814
-0.001117805
-0.046166774
-0.043482538
-0.032891743
-0.084246591
-0.163034618
-0.180114925
-0.132728443
-0.020666506
0.023993358
0.057101935
0.097326137
0.067969896
0.037049096
-0.013445113
-0.036026098
-0.006059064
-0.009637177
-0.041904949
-0.074324578
-0.040994078
0.003820812
0.046555184
0.094806388
0.130894855
0.123802692
0.093020894
0.07041081
0.051300399
0.104150183
0.139077574
0.109171256
0.087856971
0.090874746
0.043669015
-0.026082417
-0.070718139
-0.090271153
-0.054686405
-0.011102407
-0.008730405
-0.030453941
-0.040173266
-0.081435867
-0.105581425
-0.100748003
0.002363411
0.085470378
0.065989852
0.00703148
-0.062507607
-0.118311398
-0.185841426
-0.261439562
-0.204148144
-0.121601216
-0.129147977
-0.120874442
-0.092257306
-0.083272561
-0.101775713
-0.087224834
-0.084016688
-0.03715989
-0.037122864
-0.021511286
0.019027874
0.030642711
0.094267048
0.09768638
0.043278068
0.005073365
-0.034308966
-0.086377203
-0.120629065
-0.096118823
-0.08076106
-0.066955946
-0.037291087
-0.012826661
-0.008682664
-0.021895736
-0.014300513
-0.022467403
-0.089986406
-0.139530033
-0.13871792
-0.084680237
-0.11034327
-0.201704487
-0.30289501
-0.324801236
-0.325032264
-0.26242283
-0.202887952
-0.201770455
-0.17418848
-0.088007741
-0.003059228
0.032020666
0.068806455
0.057686087
0.046247415
0.043956716
0.024245201
0.063935123
0.135090232
0.102608033
0.094928138
0.107743539
0.150220573
0.100394621
0.062285397
0.121374965
0.147366598
0.126741931
0.109674454
0.115365371
0.09204036
0.086267754
0.072023772
0.079529129
0.120979413
0.113247909
0.143327937
0.154793218
0.090492867
0.131235674
0.180338159
0.151288286
0.079263762
0.048090935
0.048317637
0.004398799
-0.027263688
-0.02934574
0.004469914
0.019339245
0.048417483
0.029959729
-0.034633715
-0.095585838
-0.1089551
-0.076651089
-0.048219785
0.012914742
0.019559389
0.026984775
0.019400971
0.012449955
-0.048830461
-0.062843598
-0.052444413
-0.086460374
-0.063545302
-0.050699491
-0.050006673
-0.06862852
-0.022757636
0.015960988
-0.02019505
0.011616307
-0.026016768
-0.099318758
-0.077964425
-0.108473547
-0.187139869
-0.184651673
-0.210732415
-0.196854904
-0.173544139
-0.131256714
-0.061608061
-0.049864002
-0.064820744
-0.080722339
-0.112692863
-0.133500084
-0.140668735
-0.075191788
-0.039797798
-0.011747966
-0.040329039
-0.040258296
0.041721933
0.05634002
0.064489082
0.143608868
0.194600448
0.185345992
0.140834123
0.102522954
0.101124063
0.1007782
0.075491637
0.030221513
-0.05723872
-0.14431864
-0.147751987
-0.085598432
-0.014448578
0.050872028
0.071015991
0.039098445
0.017263288
0.027767744
0.09608493
0.138993308
0.178548291
0.194834143
0.235413656
0.28559199
0.278368473
0.241123348
0.195185706
0.192364916
0.218187928
0.15810214
0.115447119
0.10493084
0.116331413
0.058630478
-0.008568316
-5.12826E-05
-0.009841596
-0.046769582
-0.027373793
-0.036681421
-0.042938203
0.040492132
0.080075823
0.051116094
0.021781823
-0.002876388
0.002587873
0.044006161
0.086347304
0.14829576
0.155643985
0.129751131
0.138760313
0.142003611
0.11821457
0.053108834
0.057142496
0.046112221
0.008114048
-0.015750512
-0.054158125
-0.037632015
0.056993362
0.118817538
0.125830159
0.110355258
0.086432345
0.094824009
0.159679249
0.183181509
0.160581425
0.107249901
0.106441706
0.081070699
0.090700857
0.10587208
0.113140173
0.069613628
-0.021600178
-0.086729333
-0.074569538
-0.03670546
-0.025952043
-0.025341792
-0.017940046
-0.041510284
-0.078825884
-0.09567564
-0.077603213
-0.087469675
-0.119701006
-0.125730723
-0.142107084
-0.164294183
-0.179353222
-0.114884704
-0.066498727
-0.030455578
0.005234221
0.026741901
0.054554891
0.046620645
0.065481253
0.072705068
0.013877688
-0.033211131
-0.064842984
-0.092232481
-0.041157357
-0.019848563
-0.042716473
-0.048254486
-0.04876975
-0.084143624
-0.137380674
-0.096445225
-0.065584227
-0.06128988
-0.046271015
-0.065170251
-0.025356736
0.047375001
0.120975353
0.169803083
0.184528857
0.176969528
0.176275581
0.144725472
0.055322904
0.044384498
0.038965102
-0.022405334
-0.090040676
-0.155185744
-0.179570258
-0.1929162
-0.176576242
-0.197918341
-0.1887687
-0.173831031
-0.17545785
-0.199336097
-0.187993795
-0.155047566
-0.1429822
-0.113755479
-0.06258145
-0.004203648
0.019878723
0.012724577
0.023923703
0.054913577
0.117840946
0.146716267
0.1926357
0.242313862
0.230213448
0.195051193
0.100078009
0.025579626
0.011733852
0.023595551
0.0457788
0.032576013
-0.020310951
-0.007586265
0.050038945
0.093871713
0.099238396
0.079350919
0.052090898
0.036560528
-0.006635576
-0.029518725
-0.043425117
-0.086670779
-0.13623324
-0.100596316
-0.047373924
-0.044101406
0.008511082
0.033392049
-0.019582191
-0.075836107
-0.097073376
-0.092587486
-0.060327388
-0.047091443
-0.070770912
-0.057603657
-0.04551338
-0.064614408
-0.085414372
-0.112792738
-0.113121934
-0.085036613
-0.060319513
-0.041246455
-0.031668466
-0.107986391
-0.112640284
-0.126404926
-0.144257054
-0.143416584
-0.099573471
-0.064424969
-0.072886191
-0.055253938
-0.031040922
-0.026168277
0.029154763
0.07232888
0.067512482
0.032398552
-0.02976661
-0.026691768
0.011687158
0.063227437
0.110953569
0.095487751
0.116631307
0.148710206
0.174526289
0.196413264
0.164830521
0.116937436
0.044913549
-0.012157983
-0.014610712
0.022349847
0.033670813
-0.003500326
-0.051819861
-0.062017012
-0.048514999
-0.032721531
-0.058883019
-0.045337498
-0.011548349
0.015238414
0.021446735
0.029786469
0.029372938
0.005049209
-8.92665E-05
-0.008858737
0.07297837
0.148917764
0.134035587
0.109646492
0.079387926
0.033775467
0.027442873
0.0365283
0.019201942
0.000399996
0.020126875
0.052452151
0.065390579
0.071265906
0.062632121
0.045868997
0.032044679
0.035117988
0.084628746
0.072853141
0.034492847
0.004638158
-0.063108526
-0.079899095
-0.105983891
-0.104891621
-0.094010435
-0.076798163
-0.091009684
-0.086414464
-0.025358805
0.019589148
-0.002057199
-0.022010269
0.004198944
0.004878161
-0.005804646
-0.012077836
-0.068757474
-0.115264922
-0.104688063
-0.060936268
-0.040783819
-0.058612399
-0.055561446
-0.044478625
-0.006305431
0.004055913
0.027495543
0.054763116
0.067134947
0.113152228
0.103698328
0.082746521
0.070382051
0.055444241
0.062819354
0.066205807
0.066834673
0.083738498
0.067149684
0.052876808
0.046564113
0.039922059
0.009053633
-0.024505723
-0.05568723
-0.10449715
-0.093162492
-0.119510747
-0.13942644
-0.123777717
-0.080447532
-0.024892837
0.052430492
0.054330289
0.041896131
0.006368777
-0.008169903
0.022541611
0.071777977
0.098785311
0.091523729
0.07057748
0.109296747
0.127644554
0.080217056
0.048270755
0.059440546
0.071018711
0.056448117
0.025894012
0.041673891
0.045030132
0.070245594
0.12100216
0.104996741
0.096089333
0.09670987
0.082546346
0.106361158
0.083328053
-0.005539254
-0.099174887
-0.170026198
-0.205078214
-0.168906301
-0.172321305
-0.187197581
-0.167367041
-0.122949839
-0.06306003
-0.09026286
-0.15338622
-0.163922742
-0.194680452
-0.233320549
-0.248506129
-0.272894621
-0.293170184
-0.279560387
-0.244370505
-0.206411213
-0.169187814
-0.14188166
-0.079930775
-0.008377491
0.008218165
0.01769235
0.04569504
0.071883738
0.067376502
0.058656618
0.023486832
0.018636774
0.063103743
0.100482665
0.089651354
0.091695234
0.079580471
0.081637211
0.122598253
0.132956788
0.128254235
0.156787097
0.135799587
0.083332852
0.062552385
0.075667046
0.093988888
0.128990531
0.139926955
0.0831232
0.081486769
0.09054058
0.048574034
0.006128532
-0.001259237
0.01326308
-0.002466368
0.005632223
-0.018772101
-0.016583471
0.04889857
0.069049336
0.020649565
-0.020003486
0.000437578
-0.022624671
-0.094176002
-0.139690563
-0.105684429
-0.072725736
-0.105897695
-0.109052286
-0.05650625
-0.0549553
-0.071305379
-0.048579913
-0.021856338
-0.006317633
-0.013667645
-0.000936741
0.007055105
-0.030498981
-0.0185654
0.004094329
0.019820202
0.028838472
0.026697433
0.016591512
0.027620619
0.037026867
0.017242439
0.005761677
0.023047147
0.021282975
0.007326302
-0.03437832
-0.043663051
-0.047580905
-0.061081339
-0.058429554
-0.037997548
-0.043777779
-0.0734144
-0.081928037
-0.039530244
-0.00230236
0.029496288
0.037657645
0.048246033
0.057037286
0.048817176
0.046927944
0.025978422
0.008005762
0.019946394
0.004353557
0.009419938
-0.002877715
-0.041263059
-0.047108956
-0.028951081
-0.008516276
0.006929583
0.025201714
0.001747303
0.011650991
0.04008254
0.069369845
0.101948604
0.120137997
0.126352102
0.098372251
0.065066561
0.04918826
0.024085037
0.023863215
0.034472294
0.005103803
-0.019013701
-0.026102366
-0.037205752
-0.025471214
-0.023486705
-0.032141589
-0.037703287
-0.073522769
-0.078104511
-0.080326885
-0.057987127
-0.039427996
-0.025249777
-0.019407732
-0.040926628
-0.040955961
-0.043384865
-0.071322098
-0.104807861
-0.111549892
-0.062767997
-0.020546854
-0.005817012
0.008082216
0.027768714
0.023588169
0.016892197
0.008117477
0.003389583
-0.002834801
0.002523888
0.009303325
0.008895005
-0.000706777
0.007101542
-0.010399385
-0.02764667
-0.017987559
-0.004245691
0.011980807
0.031838935
0.023661477
-0.012129479
-0.043302063
-0.055932086
-0.045428749
-0.030915681
-0.052222736
-0.043274909
-0.013687124
0.00367392
0.007124259
0.022140902
0.052083537
0.05316804
0.04882703
0.046905242
0.054063022
0.049249094
0.058346495
0.053572547
0.05175728
0.069569513
0.077176526
0.074752957
0.089129068
0.089525983
0.069997057
0.064898834
0.04371262
0.019497918
0.020732908
-0.004969482
-0.026236393
-0.041958012
-0.055189811
-0.041575506
-0.018029019
-0.005238523
-0.002384826
-0.002878704
-0.009190834
-0.001673016
0.014990809
0.019261038
0.011834977
0.014429924
-0.006194397
-0.013690857
-0.01815436
-0.027352409
-0.011725085
-0.007274496
-0.004866365
-0.000498055
-0.024686567
-0.038879536
-0.008724798
0.020850286
0.030796662
0.042091474
0.029409859
0.019482784
0.015027554
0.025668813
0.021417856
-0.004021549
-0.008313496
0.003285972
0.015137304
0.020884132
0.020320788
0.018038591
0.0111319
0.00482287
-3.32985E-06
0.013206283
0.021511154
0.008566726
-0.014607481
-0.035937171
-0.041226365
-0.032938849
-0.036250491
-0.031613734
-0.027737988
-0.015449633
-0.025691589
-0.037606847
-0.023819452
-0.008332955
-0.02272051
-0.024265362
-0.023423638
-0.023754835
-0.008664452
0.010672833
0.009229788
-0.011631425
-0.024956243
-0.023508839
-0.030030534
-0.014050701
-0.023040876
-0.03041411
-0.007875336
0.015406552
0.017798344
0.028288897
0.04858486
0.048115905
0.034862392
0.02364371
0.013462433
-0.000195919
-0.014769358
-0.021335134
-0.029456133
-0.02629525
-0.003338757
0.007424269
0.013872341
0.025195401
0.04352821
0.040133093
0.029292908
0.026523091
0.011537333
-0.009784389
-0.001986326
-0.0032865
-0.004223608
-0.002127193
-0.009363564
-0.008282364
-0.003192856
0.000500254
0.000902361
-0.000334126
-0.014894703
-0.013578523
0.002373168
0.016081851
0.013832772
0.027076416
0.033990778
0.036758874
0.03900085
0.033391785
0.025377925
0.006117424
-0.007982285
0.001321586
0.012323154
0.000274451
-0.003341291
0.004245928
0.02437862
0.020422032
0.018299138
0.013248277
-0.004205177
-0.00940933
0.007648413
0.037467923
0.042979371
0.031528763
0.029529661
0.033087172
0.036130879
0.025318846
0.013500556
0.016358124
0.019501613
0.006909644
-0.010932402
-0.023442391
-0.023995096
-0.020144735
-0.022617867
-0.031533118
-0.031856138
-0.018209685
-0.017285462
-0.010917584
-0.005178685
0.001658815
0.012611457
0.00584742
-0.010092684
-0.010470636
-0.000809887
0.005100795
0.009637169
0.01450677
0.011180992
0.003989922
0.002567968
-0.001744421
-0.012866669
-0.026271375
-0.034311958
-0.03018233
-0.025640557
-0.023231883
-0.023810105
-0.025842901
-0.018591259
-0.010650679
-0.008763921
-0.009135894
-0.001373824
0.006697438
-0.001906455
-0.007491163
-0.006014751
0.001089316
0.006058067
0.017044736
0.015292175
0.007066178
0.00374424
-0.002431477
-0.004618892
-0.006236679
-0.001577474
0.007022256
0.010737851
0.012828908
0.01042438
0.009834042
0.006932662
0.004736242
0.008002718
0.006544638
0.006877357
0.004591276
0.001104655
0.003414839
0.008876175
0.016261095
0.024845826
0.032212958
0.030769438
0.029917592
0.026902232
0.021839717
0.017675372
0.007775919
0.001524836
0.002055481
0.007216346
0.012307151
0.017424561
0.020510783
0.013047465
0.005065563
0.002401456
-0.001621276
-0.001340179
-0.008335399
-0.003247947
0.004523735
0.005820044
0.008328085
0.008836355
0.002405991
-0.005828646
-0.011412126
-0.010274104
-0.01045389
-0.008889354
-0.013874169
-0.01940581
-0.014897407
-0.013080622
-0.012278072
-0.013192892
-0.011359205
-0.008470882
-0.007730011
-0.004104892
-0.000214774
0.00229679
0.008195203
0.014533549
0.016879939
0.016947389
0.018964881
0.020868294
0.018920347
0.012795542
0.008139875
0.007960005
0.006814465
0.004858091
0.002352709
0.000415702
-0.000690862
-0.000260775
0.000460456
-0.000123262
0.001940525
-0.000397992
-0.001015657
-0.001861165
-0.002120086
-5.65154E-05
-0.000420958
-0.000272386
0.002476442
0.004417398
0.002179623
-0.000833663
-0.002700218
-0.003614433
-0.002216537
-0.001820936
-0.003399915
-0.002240567
-0.000178083
0.001156171
0.002909538
0.004093854
0.0034255
0.002472007
0.002117974
0.000695984
4.32592E-05
-0.00034121
0.000462581
0.000567747
0.000642022
0.000763812
0.001109343
0.001665021
0.001976526
0.001763458
0.00160097
0.001316621
0.000807511
0.000856272
0.000980047
0.001368015
0.001536564
0.001563827
0.001883912
0.0020244
0.004389876
0.004138238
0.003538426
0.002956079
0.002300461
0.000958434
0.001340644
0.004399356
0.006974449
0.008374172
0.00788987
0.010585962
0.011566041
0.009968489
0.010941475
0.007110497
0.006148309
0.006578479
0.008128808
0.013749401
0.019004095
0.014562379
0.0069579
0.005822071
0.008795188
0.005871725
0.002216147
0.006665813
0.009231095
0.010979309
0.011716109
0.011303527
0.010068906
0.007797371
0.003007371
0.002053698
0.010031558
0.026374172
0.046171248
0.052403793
0.052932024
0.058454975
0.058768585
0.067543961
0.066608369
0.053094506
0.037442166
0.032747816
0.029241055
0.022598119
0.022125896
0.014702437
-0.002984791
-0.004900933
-0.013001241
-0.023455277
-0.013545915
-0.002640109
0.004623368
0.003935909
-0.009988242
-0.026583966
-0.044086028
-0.043933187
-0.025314456
-0.023492301
-0.040920228
-0.043290861
-0.040735353
-0.047085773
-0.031774867
-0.003795079
0.004417128
0.017984247
0.036742911
0.037315238
0.022797396
0.016780613
0.012980396
0.018816974
0.031083621
0.021892844
0.030406939
0.003875493
-0.049668845
-0.096326455
-0.099464335
-0.075044557
-0.068058118
-0.070133775
-0.060190443
-0.06692908
-0.050251424
-0.022884341
0.030718228
0.047963366
0.028263319
0.019636098
-0.009945349
-0.006615837
0.042356066
0.021618096
-0.025133919
-0.039636467
-0.036101315
-0.063346729
-0.109972924
-0.100563735
-0.052185543
-0.025679411
-0.01020145
0.019292211
0.055217765
0.075294882
0.044998959
0.042072207
0.034597144
0.034453262
0.043323968
0.060132377
0.035380051
0.046251744
0.038187187
0.011047371
0.009233863
0.034305055
0.021512978
0.024317242
0.025100769
0.055347566
0.082517222
0.099676646
0.11055956
0.095017403
0.091345496
0.08531516
0.091366097
0.058829732
-0.02147487
-0.062186196
-0.07657595
-0.07898993
-0.047280844
-0.023482857
-0.018443087
-0.005183474
0.00067654
-0.002617953
-0.044077925
-0.114724182
-0.141991869
-0.159872487
-0.217377454
-0.224519059
-0.178097233
-0.175022557
-0.176696211
-0.136746496
-0.123677075
-0.085693039
-0.013966632
0.035090212
0.050612431
0.045808882
0.072029606
0.064596757
0.035395753
0.10086754
0.111724965
0.039512467
-0.004864808
0.034604825
0.057628743
0.085553154
0.099744268
0.100784816
0.086710095
0.031881507
-0.020293508
-0.055310227
-0.023913372
-0.033108454
-0.054750212
-0.038600005
-0.019289847
0.025195153
0.031383514
0.031031655
0.00732173
-0.036314026
-0.015704393
0.04125287
0.006782964
-0.004442574
0.032867163
0.086281635
0.206558928
0.238221645
0.206363082
0.054481227
-0.034699071
-0.077240109
-0.068054847
-0.011507249
0.027141836
0.068018384
0.042309061
0.018113339
0.035914503
0.042923145
-0.005328159
0.003197319
-0.010312952
-0.03201396
0.039041474
0.085358322
0.069810122
0.019421455
-0.008638076
-0.015821826
-0.026125092
-0.09495388
-0.134317532
-0.104278386
-0.043131333
-0.041233797
-0.025282316
-0.024528721
-0.015251717
0.005746465
0.041982014
0.036497783
0.012848216
-0.011986255
-0.053300429
-0.05255511
-0.072286196
-0.067120001
-0.08280351
-0.095052712
-0.074766047
-0.068195023
-0.022107238
0.051868174
0.042095818
0.045275349
0.087857679
0.053255603
0.007478693
-0.006987691
-0.032102179
-0.054716956
-0.049552917
-0.014691803
0.039705373
0.105895624
0.116678119
0.049160235
0.063530371
0.131679595
0.164959475
0.185520545
0.175114229
0.18942976
0.187566891
0.183584258
0.129206046
0.086754017
0.094620831
0.067868114
0.091821954
0.154070199
0.169293538
0.17330642
0.160751805
0.116935693
0.074370682
0.068655297
0.099135317
0.062715508
0.025828622
-0.001267581
-0.031657115
-0.082366355
-0.075164802
-0.069023564
-0.077361226
-0.18009603
-0.238762692
-0.220637485
-0.179838777
-0.228515774
-0.228194639
-0.138895318
-0.140812069
-0.12339206
-0.064908393
-0.078718267
-0.062946007
-0.00103504
0.007666875
0.006007117
-0.009354099
-0.049894448
-0.016770985
-0.012860426
-0.079803064
-0.141280085
-0.110163137
-0.108224563
-0.111272588
-0.040706016
0.002670593
0.035435319
0.062670484
0.060045082
0.078587346
0.047201227
-0.012014098
-0.064838246
-0.089429051
-0.105027899
-0.113477737
-0.127376944
-0.121771365
-0.118962914
-0.160575449
-0.204382554
-0.253239125
-0.268444806
-0.249026641
-0.195515603
-0.161946267
-0.143055901
-0.087561198
-0.128968492
-0.154228136
-0.089996606
-0.08030273
-0.07514099
-0.069628678
-0.088012442
-0.1297618
-0.103878252
-0.019945446
0.038701799
0.029131671
0.038021017
0.102432355
0.128920987
0.069437489
0.010602042
-0.000483863
0.024478998
0.038328841
0.040727396
0.032824844
-0.015539099
-0.0503456
-0.014558623
-0.00748615
-0.034387209
-0.010439106
-0.013288285
-0.030585457
0.036311362
0.030015795
0.03252903
0.025640277
-0.008299512
0.011420015
0.036547035
0.007053625
-0.012425179
-0.042136468
-0.021312159
0.007412418
-0.050826088
-0.100183763
-0.115574025
-0.105981104
-0.08848273
-0.10301844
-0.125949264
-0.131981537
-0.087121747
-0.062966332
-0.044291142
-0.065409757
-0.016360315
0.055094779
0.099504277
0.169816658
0.148446813
0.053961568
-0.023430729
-0.071043953
-0.027465085
0.013452497
-0.010729436
-0.005107792
-0.001314264
-0.019754402
0.010598665
-0.001760871
0.053848825
0.047590174
0.024093928
0.023419751
0.020481698
0.066653065
0.114673659
0.071195804
0.104845189
0.154766858
0.146247268
0.198285833
0.268434793
0.276665241
0.220919564
0.150192663
0.070541821
0.022955595
0.037346646
0.052325886
0.02868638
-0.004585291
-0.040150806
-0.035184436
0.015591319
0.059459414
0.017683357
-0.059171733
-0.056212962
-0.058371
-0.09174601
-0.049149513
-0.077831045
-0.113683842
-0.075383417
-0.076165192
-0.132980496
-0.227275655
-0.184773013
-0.094276026
-0.070428468
-0.110784523
-0.137689441
-0.109232582
-0.079811096
-0.00835235
-0.025476513
-0.060736451
-0.027161319
-0.04330321
-0.019583208
0.006330533
0.013584677
0.008081361
0.017939458
0.054331191
0.040236071
0.02795046
0.059853215
0.084023401
0.160413429
0.190545768
0.203166395
0.185712293
0.175010711
0.133382872
0.11645788
0.113888092
0.065955505
0.115141883
0.194195032
0.207252294
0.162922472
0.162107989
0.129412577
0.037308112
0.001289477
-0.040498033
-0.055742372
-0.075794503
-0.094154865
-0.141408473
-0.164765954
-0.160227999
-0.186974078
-0.169996709
-0.135958284
-0.038093958
0.056442861
0.101021044
0.1292101
0.105391487
0.029447995
0.025793614
0.040120065
0.071422867
0.052875195
0.008043284
-0.008641883
-0.009572518
0.063596122
0.123738088
0.091491133
0.053688508
0.044477422
0.032545213
0.053915787
0.024620835
0.012838631
0.033009417
0.044526979
0.003437634
0.037034113
0.06787841
0.127609536
0.182717964
0.177482486
0.19061695
0.190846875
0.145910621
0.128924981
0.132716596
0.118908569
0.101858824
0.119268984
0.138784766
0.15850693
0.098713771
0.003999926
-0.055424936
-0.102571249
-0.148509979
-0.092618912
-0.01123844
0.043049876
0.050097577
0.039061014
0.051775746
0.039264865
-0.013775383
-0.038203597
-0.009675848
-0.046133373
0.009924787
0.053420857
-0.010197366
-0.058407374
-0.050644141
-0.04937207
-0.09228991
-0.121362038
-0.188260898
-0.17861031
-0.151900783
-0.106594056
-0.082105756
-0.068192683
-0.056719352
-0.072887599
-0.089877538
-0.139100552
-0.164130613
-0.135225073
-0.125551999
-0.147427708
-0.153776288
-0.082781613
-0.006815014
0.027425965
0.028930675
0.036038376
0.094926327
0.083093755
0.022690196
-0.02233386
-0.041778479
-0.008570619
0.05895713
0.114532627
0.087312065
0.090907335
0.10739664
0.039083377
0.040633582
0.047920797
0.059957188
0.064590298
0.051486667
0.062788874
0.083931699
0.093663625
0.100586951
0.048265848
-0.003041288
-0.0725343
-0.149264693
-0.164445162
-0.180927724
-0.165441051
-0.12300463
-0.072564319
0.006662579
-0.033231743
-0.042434316
0.018175816
0.054569922
0.047878623
-0.018910987
-0.01950888
-0.049443353
-0.142873123
-0.220553383
-0.272845298
-0.259322971
-0.235190585
-0.225628868
-0.216409728
-0.195022687
-0.13942261
-0.116300173
-0.142197013
-0.168836325
-0.161538079
-0.156031296
-0.115752056
-0.080743782
0.004815917
0.084683761
0.044227317
-0.030619025
-0.012500649
0.075133547
0.178499013
0.222744703
0.209623083
0.28842479
0.350096583
0.299970031
0.250269502
0.196404546
0.171851456
0.167426467
0.246651828
0.268478811
0.218463749
0.199080229
0.198228478
0.13845484
0.097431853
0.068056002
0.049522586
0.053581707
0.028624281
0.013652835
-0.025569383
-0.123640411
-0.166638419
-0.148583293
-0.066730812
-0.00221677
0.05905636
0.04027288
0.055814914
0.07306575
0.082259335
0.116849631
0.106711403
0.057184469
0.008970197
-0.014485115
-0.00316432
0.000828466
-0.017712196
-0.076118156
-0.08145947
-0.02739669
0.017521003
0.050027005
0.07927978
0.135163993
0.165154636
0.184301138
0.161105707
0.175331637
0.19975163
0.175987929
0.166578814
0.158831343
0.171066731
0.148071662
0.063943021
-0.002897062
0.004464262
0.029098753
0.062108021
0.086413361
0.108353168
0.139394611
0.179077923
0.198447466
0.173068792
0.101427764
0.002875132
-0.122601017
-0.155247912
-0.201861262
-0.347599924
-0.350593895
-0.333272249
-0.288043499
-0.24214229
-0.220248953
-0.151501939
-0.134753287
-0.168050721
-0.186690286
-0.140635937
-0.082204916
-0.050225515
-0.006857778
-0.001305159
-0.022105826
-0.049406536
-0.107557513
-0.132905975
-0.070751436
-0.047414422
-0.08698988
-0.054657914
-0.01220632
0.053120159
0.116007313
0.117257677
0.108856715
0.112051934
0.124573335
0.127473161
0.214074612
0.262830645
0.238227755
0.200920805
0.245996028
0.24187395
0.119590208
0.010628896
-0.038126141
-0.078206673
-0.150070593
-0.197204709
-0.219218016
-0.291837543
-0.323337704
-0.315931797
-0.249574289
-0.199295342
-0.189535275
-0.169362649
-0.104707181
-0.104726233
-0.116100237
-0.102004662
-0.052015252
0.004629215
0.003209735
0.048084844
0.101253636
0.084979415
0.054007139
0.01546963
0.024102589
0.033831432
0.042306863
0.084896609
0.08665695
0.048973523
-0.011994101
-0.013486257
0.016142447
-0.025257366
-0.077199802
-0.095984653
-0.113123663
-0.139941931
-0.14556393
-0.155254737
-0.090928674
-0.041798696
-0.04139502
-0.012623093
0.06074708
0.095009409
0.164344758
0.241050079
0.235544905
0.253403604
0.266828448
0.226403415
0.175562873
0.205246001
0.261426687
0.247781545
0.188658938
0.205553874
0.240348041
0.207255423
0.168956161
0.141112164
0.091098554
0.05015352
0.020443449
-0.027955413
-0.059405155
-0.025786044
0.018179359
0.022815011
-0.00149671
-0.06096765
-0.07277213
-0.023560327
0.002005536
0.045602906
0.056839135
0.042840108
0.044762373
0.048839685
0.047889937
-0.020861952
0.024696078
0.099641278
0.055086851
-0.020806495
-0.003074558
0.009929988
-0.019656239
-0.078323901
-0.075673029
-0.085612819
-0.087740131
-0.032917436
-0.009666321
-0.050624903
-0.106706128
-0.138963774
-0.178197011
-0.185912162
-0.1700259
-0.210941702
-0.237469852
-0.249348566
-0.226435363
-0.194474086
-0.135895357
-0.012829147
0.085702844
0.079072699
0.051635355
0.056965314
0.006308983
-0.032660849
-0.026627831
-0.099196374
-0.104629628
-0.011109727
0.021691293
0.004987638
-0.026379542
-0.056380115
-0.027281251
0.04713764
0.093570001
0.084913194
0.069676168
0.064123839
0.036445986
0.021342762
-0.0249002
-0.012444232
-0.009573135
0.043087974
0.04736663
-0.013752058
0.021025168
0.024697382
0.001431413
0.063563272
0.075505018
0.048893183
0.035958856
-0.035261717
-0.113736808
-0.15111433
-0.122489884
-0.152192876
-0.183135793
-0.161948785
-0.159887642
-0.077011727
-0.044141281
-0.047557507
0.036005914
0.145909786
0.116012923
0.088287607
0.103065729
0.02945335
-0.020842232
-0.086875081
-0.140062749
-0.135303915
-0.088159755
-0.088365674
-0.058485456
-0.08154016
-0.091838963
-0.031374894
-0.014583507
-0.006194479
-0.017967548
-0.070953757
-0.063599363
-0.021371901
0.007267625
0.029860955
0.003405939
0.058771208
0.153367743
0.126692235
0.211100966
0.186623365
0.143556833
0.086655468
0.043308221
0.039305814
-0.001117805
-0.046166774
-0.043482538
-0.032891743
-0.084246591
-0.163034618
-0.180114925
-0.132728443
-0.020666506
0.023993358
0.057101935
0.097326137
0.067969896
0.037049096
-0.013445113
-0.036026098
-0.006059064
-0.009637177
-0.041904949
-0.074324578
-0.040994078
0.003820812
0.046555184
0.094806388
0.130894855
0.123802692
0.093020894
0.07041081
0.051300399
0.104150183
0.139077574
0.109171256
0.087856971
0.090874746
0.043669015
-0.026082417
-0.070718139
-0.090271153
-0.054686405
-0.011102407
-0.008730405
-0.030453941
-0.040173266
-0.081435867
-0.105581425
-0.100748003
0.002363411
0.085470378
0.065989852
0.00703148
-0.062507607
-0.118311398
-0.185841426
-0.261439562
-0.204148144
-0.121601216
-0.129147977
-0.120874442
-0.092257306
-0.083272561
-0.101775713
-0.087224834
-0.084016688
-0.03715989
-0.037122864
-0.021511286
0.019027874
0.030642711
0.094267048
0.09768638
0.043278068
0.005073365
-0.034308966
-0.086377203
-0.120629065
-0.096118823
-0.08076106
-0.066955946
-0.037291087
-0.012826661
-0.008682664
-0.021895736
-0.014300513
-0.022467403
-0.089986406
-0.139530033
-0.13871792
-0.084680237
-0.11034327
-0.201704487
-0.30289501
-0.324801236
-0.325032264
-0.26242283
-0.202887952
-0.201770455
-0.17418848
-0.088007741
-0.003059228
0.032020666
0.068806455
0.057686087
0.046247415
0.043956716
0.024245201
0.063935123
0.135090232
0.102608033
0.094928138
0.107743539
0.150220573
0.100394621
0.062285397
0.121374965
0.147366598
0.126741931
0.109674454
0.115365371
0.09204036
0.086267754
0.072023772
0.079529129
0.120979413
0.113247909
0.143327937
0.154793218
0.090492867
0.131235674
0.180338159
0.151288286
0.079263762
0.048090935
0.048317637
0.004398799
-0.027263688
-0.02934574
0.004469914
0.019339245
0.048417483
0.029959729
-0.034633715
-0.095585838
-0.1089551
-0.076651089
-0.048219785
0.012914742
0.019559389
0.026984775
0.019400971
0.012449955
-0.048830461
-0.062843598
-0.052444413
-0.086460374
-0.063545302
-0.050699491
-0.050006673
-0.06862852
-0.022757636
0.015960988
-0.02019505
0.011616307
-0.026016768
-0.099318758
-0.077964425
-0.108473547
-0.187139869
-0.184651673
-0.210732415
-0.196854904
-0.173544139
-0.131256714
-0.061608061
-0.049864002
-0.064820744
-0.080722339
-0.112692863
-0.133500084
-0.140668735
-0.075191788
-0.039797798
-0.011747966
-0.040329039
-0.040258296
0.041721933
0.05634002
0.064489082
0.143608868
0.194600448
0.185345992
0.140834123
0.102522954
0.101124063
0.1007782
0.075491637
0.030221513
-0.05723872
-0.14431864
-0.147751987
-0.085598432
-0.014448578
0.050872028
0.071015991
0.039098445
0.017263288
0.027767744
0.09608493
0.138993308
0.178548291
0.194834143
0.235413656
0.28559199
0.278368473
0.241123348
0.195185706
0.192364916
0.218187928
0.15810214
0.115447119
0.10493084
0.116331413
0.058630478
-0.008568316
-5.12826E-05
-0.009841596
-0.046769582
-0.027373793
-0.036681421
-0.042938203
0.040492132
0.080075823
0.051116094
0.021781823
-0.002876388
0.002587873
0.044006161
0.086347304
0.14829576
0.155643985
0.129751131
0.138760313
0.142003611
0.11821457
0.053108834
0.057142496
0.046112221
0.008114048
-0.015750512
-0.054158125
-0.037632015
0.056993362
0.118817538
0.125830159
0.110355258
0.086432345
0.094824009
0.159679249
0.183181509
0.160581425
0.107249901
0.106441706
0.081070699
0.090700857
0.10587208
0.113140173
0.069613628
-0.021600178
-0.086729333
-0.074569538
-0.03670546
-0.025952043
-0.025341792
-0.017940046
-0.041510284
-0.078825884
-0.09567564
-0.077603213
-0.087469675
-0.119701006
-0.125730723
-0.142107084
-0.164294183
-0.179353222
-0.114884704
-0.066498727
-0.030455578
0.005234221
0.026741901
0.054554891
0.046620645
0.065481253
0.072705068
0.013877688
-0.033211131
-0.064842984
-0.092232481
-0.041157357
-0.019848563
-0.042716473
-0.048254486
-0.04876975
-0.084143624
-0.137380674
-0.096445225
-0.065584227
-0.06128988
-0.046271015
-0.065170251
-0.025356736
0.047375001
0.120975353
0.169803083
0.184528857
0.176969528
0.176275581
0.144725472
0.055322904
0.044384498
0.038965102
-0.022405334
-0.090040676
-0.155185744
-0.179570258
-0.1929162
-0.176576242
-0.197918341
-0.1887687
-0.173831031
-0.17545785
-0.199336097
-0.187993795
-0.155047566
-0.1429822
-0.113755479
-0.06258145
-0.004203648
0.019878723
0.012724577
0.023923703
0.054913577
0.117840946
0.146716267
0.1926357
0.242313862
0.230213448
0.195051193
0.100078009
0.025579626
0.011733852
0.023595551
0.0457788
0.032576013
-0.020310951
-0.007586265
0.050038945
0.093871713
0.099238396
0.079350919
0.052090898
0.036560528
-0.006635576
-0.029518725
-0.043425117
-0.086670779
-0.13623324
-0.100596316
-0.047373924
-0.044101406
0.008511082
0.033392049
-0.019582191
-0.075836107
-0.097073376
-0.092587486
-0.060327388
-0.047091443
-0.070770912
-0.057603657
-0.04551338
-0.064614408
-0.085414372
-0.112792738
-0.113121934
-0.085036613
-0.060319513
-0.041246455
-0.031668466
-0.107986391
-0.112640284
-0.126404926
-0.144257054
-0.143416584
-0.099573471
-0.064424969
-0.072886191
-0.055253938
-0.031040922
-0.026168277
0.029154763
0.07232888
0.067512482
0.032398552
-0.02976661
-0.026691768
0.011687158
0.063227437
0.110953569
0.095487751
0.116631307
0.148710206
0.174526289
0.196413264
0.164830521
0.116937436
0.044913549
-0.012157983
-0.014610712
0.022349847
0.033670813
-0.003500326
-0.051819861
-0.062017012
-0.048514999
-0.032721531
-0.058883019
-0.045337498
-0.011548349
0.015238414
0.021446735
0.029786469
0.029372938
0.005049209
-8.92665E-05
-0.008858737
0.07297837
0.148917764
0.134035587
0.109646492
0.079387926
0.033775467
0.027442873
0.0365283
0.019201942
0.000399996
0.020126875
0.052452151
0.065390579
0.071265906
0.062632121
0.045868997
0.032044679
0.035117988
0.084628746
0.072853141
0.034492847
0.004638158
-0.063108526
-0.079899095
-0.105983891
-0.104891621
-0.094010435
-0.076798163
-0.091009684
-0.086414464
-0.025358805
0.019589148
-0.002057199
-0.022010269
0.004198944
0.004878161
-0.005804646
-0.012077836
-0.068757474
-0.115264922
-0.104688063
-0.060936268
-0.040783819
-0.058612399
-0.055561446
-0.044478625
-0.006305431
0.004055913
0.027495543
0.054763116
0.067134947
0.113152228
0.103698328
0.082746521
0.070382051
0.055444241
0.062819354
0.066205807
0.066834673
0.083738498
0.067149684
0.052876808
0.046564113
0.039922059
0.009053633
-0.024505723
-0.05568723
-0.10449715
-0.093162492
-0.119510747
-0.13942644
-0.123777717
-0.080447532
-0.024892837
0.052430492
0.054330289
0.041896131
0.006368777
-0.008169903
0.022541611
0.071777977
0.098785311
0.091523729
0.07057748
0.109296747
0.127644554
0.080217056
0.048270755
0.059440546
0.071018711
0.056448117
0.025894012
0.041673891
0.045030132
0.070245594
0.12100216
0.104996741
0.096089333
0.09670987
0.082546346
0.106361158
0.083328053
-0.005539254
-0.099174887
-0.170026198
-0.205078214
-0.168906301
-0.172321305
-0.187197581
-0.167367041
-0.122949839
-0.06306003
-0.09026286
-0.15338622
-0.163922742
-0.194680452
-0.233320549
-0.248506129
-0.272894621
-0.293170184
-0.279560387
-0.244370505
-0.206411213
-0.169187814
-0.14188166
-0.079930775
-0.008377491
0.008218165
0.01769235
0.04569504
0.071883738
0.067376502
0.058656618
0.023486832
0.018636774
0.063103743
0.100482665
0.089651354
0.091695234
0.079580471
0.081637211
0.122598253
0.132956788
0.128254235
0.156787097
0.135799587
0.083332852
0.062552385
0.075667046
0.093988888
0.128990531
0.139926955
0.0831232
0.081486769
0.09054058
0.048574034
0.006128532
-0.001259237
0.01326308
-0.002466368
0.005632223
-0.018772101
-0.016583471
0.04889857
0.069049336
0.020649565
-0.020003486
0.000437578
-0.022624671
-0.094176002
-0.139690563
-0.105684429
-0.072725736
-0.105897695
-0.109052286
-0.05650625
-0.0549553
-0.071305379
-0.048579913
-0.021856338
-0.006317633
-0.013667645
-0.000936741
0.007055105
-0.030498981
-0.0185654
0.004094329
0.019820202
0.028838472
0.026697433
0.016591512
0.027620619
0.037026867
0.017242439
0.005761677
0.023047147
0.021282975
0.007326302
-0.03437832
-0.043663051
-0.047580905
-0.061081339
-0.058429554
-0.037997548
-0.043777779
-0.0734144
-0.081928037
-0.039530244
-0.00230236
0.029496288
0.037657645
0.048246033
0.057037286
0.048817176
0.046927944
0.025978422
0.008005762
0.019946394
0.004353557
0.009419938
-0.002877715
-0.041263059
-0.047108956
-0.028951081
-0.008516276
0.006929583
0.025201714
0.001747303
0.011650991
0.04008254
0.069369845
0.101948604
0.120137997
0.126352102
0.098372251
0.065066561
0.04918826
0.024085037
0.023863215
0.034472294
0.005103803
-0.019013701
-0.026102366
-0.037205752
-0.025471214
-0.023486705
-0.032141589
-0.037703287
-0.073522769
-0.078104511
-0.080326885
-0.057987127
-0.039427996
-0.025249777
-0.019407732
-0.040926628
-0.040955961
-0.043384865
-0.071322098
-0.104807861
-0.111549892
-0.062767997
-0.020546854
-0.005817012
0.008082216
0.027768714
0.023588169
0.016892197
0.008117477
0.003389583
-0.002834801
0.002523888
0.009303325
0.008895005
-0.000706777
0.007101542
-0.010399385
-0.02764667
-0.017987559
-0.004245691
0.011980807
0.031838935
0.023661477
-0.012129479
-0.043302063
-0.055932086
-0.045428749
-0.030915681
-0.052222736
-0.043274909
-0.013687124
0.00367392
0.007124259
0.022140902
0.052083537
0.05316804
0.04882703
0.046905242
0.054063022
0.049249094
0.058346495
0.053572547
0.05175728
0.069569513
0.077176526
0.074752957
0.089129068
0.089525983
0.069997057
0.064898834
0.04371262
0.019497918
0.020732908
-0.004969482
-0.026236393
-0.041958012
-0.055189811
-0.041575506
-0.018029019
-0.005238523
-0.002384826
-0.002878704
-0.009190834
-0.001673016
0.014990809
0.019261038
0.011834977
0.014429924
-0.006194397
-0.013690857
-0.01815436
-0.027352409
-0.011725085
-0.007274496
-0.004866365
-0.000498055
-0.024686567
-0.038879536
-0.008724798
0.020850286
0.030796662
0.042091474
0.029409859
0.019482784
0.015027554
0.025668813
0.021417856
-0.004021549
-0.008313496
0.003285972
0.015137304
0.020884132
0.020320788
0.018038591
0.0111319
0.00482287
-3.32985E-06
0.013206283
0.021511154
0.008566726
-0.014607481
-0.035937171
-0.041226365
-0.032938849
-0.036250491
-0.031613734
-0.027737988
-0.015449633
-0.025691589
-0.037606847
-0.023819452
-0.008332955
-0.02272051
-0.024265362
-0.023423638
-0.023754835
-0.008664452
0.010672833
0.009229788
-0.011631425
-0.024956243
-0.023508839
-0.030030534
-0.014050701
-0.023040876
-0.03041411
-0.007875336
0.015406552
0.017798344
0.028288897
0.04858486
0.048115905
0.034862392
0.02364371
0.013462433
-0.000195919
-0.014769358
-0.021335134
-0.029456133
-0.02629525
-0.003338757
0.007424269
0.013872341
0.025195401
0.04352821
0.040133093
0.029292908
0.026523091
0.011537333
-0.009784389
-0.001986326
-0.0032865
-0.004223608
-0.002127193
-0.009363564
-0.008282364
-0.003192856
0.000500254
0.000902361
-0.000334126
-0.014894703
-0.013578523
0.002373168
0.016081851
0.013832772
0.027076416
0.033990778
0.036758874
0.03900085
0.033391785
0.025377925
0.006117424
-0.007982285
0.001321586
0.012323154
0.000274451
-0.003341291
0.004245928
0.02437862
0.020422032
0.018299138
0.013248277
-0.004205177
-0.00940933
0.007648413
0.037467923
0.042979371
0.031528763
0.029529661
0.033087172
0.036130879
0.025318846
0.013500556
0.016358124
0.019501613
0.006909644
-0.010932402
-0.023442391
-0.023995096
-0.020144735
-0.022617867
-0.031533118
-0.031856138
-0.018209685
-0.017285462
-0.010917584
-0.005178685
0.001658815
0.012611457
0.00584742
-0.010092684
-0.010470636
-0.000809887
0.005100795
0.009637169
0.01450677
0.011180992
0.003989922
0.002567968
-0.001744421
-0.012866669
-0.026271375
-0.034311958
-0.03018233
-0.025640557
-0.023231883
-0.023810105
-0.025842901
-0.018591259
-0.010650679
-0.008763921
-0.009135894
-0.001373824
0.006697438
-0.001906455
-0.007491163
-0.006014751
0.001089316
0.006058067
0.017044736
0.015292175
0.007066178
0.00374424
-0.002431477
-0.004618892
-0.006236679
-0.001577474
0.007022256
0.010737851
0.012828908
0.01042438
0.009834042
0.006932662
0.004736242
0.008002718
0.006544638
0.006877357
0.004591276
0.001104655
0.003414839
0.008876175
0.016261095
0.024845826
0.032212958
0.030769438
0.029917592
0.026902232
0.021839717
0.017675372
0.007775919
0.001524836
0.002055481
0.007216346
0.012307151
0.017424561
0.020510783
0.013047465
0.005065563
0.002401456
-0.001621276
-0.001340179
-0.008335399
-0.003247947
0.004523735
0.005820044
0.008328085
0.008836355
0.002405991
-0.005828646
-0.011412126
-0.010274104
-0.01045389
-0.008889354
-0.013874169
-0.01940581
-0.014897407
-0.013080622
-0.012278072
-0.013192892
-0.011359205
-0.008470882
-0.007730011
-0.004104892
-0.000214774
0.00229679
0.008195203
0.014533549
0.016879939
0.016947389
0.018964881
0.020868294
0.018920347
0.012795542
0.008139875
0.007960005
0.006814465
0.004858091
0.002352709
0.000415702
-0.000690862
-0.000260775
0.000460456
-0.000123262
0.001940525
-0.000397992
-0.001015657
-0.001861165
-0.002120086
-5.65154E-05
-0.000420958
-0.000272386
0.002476442
0.004417398
0.002179623
-0.000833663
-0.002700218
-0.003614433
-0.002216537
-0.001820936
-0.003399915
-0.002240567
-0.000178083
0.001156171
0.002909538
0.004093854
0.0034255
0.002472007
0.002117974
0.000695984
4.32592E-05
-0.00034121
0.000462581
0.000567747
0.000642022
0.000763812
0.001109343
0.001665021
0.001976526
0.001763458
0.00160097
0.001316621
0.000807511
0.000856272
0.000980047
0.001368015
0.001536564
0.001563827
0.001883912
0.0020244
-
adamdressller80
- Posts: 2
- Joined: Tue Jan 04, 2011 8:17 am
- Contact:
Re: CTestEnergyIncr::test() - failed to converge
vesna wrote:
> Yes, at times that are shown (20.16, 20.17,20.2,....) you did not achieve
> the convergence. Thus, from time 20.16 and on, your results are not
> reliable.
thats bad it failed
> Yes, at times that are shown (20.16, 20.17,20.2,....) you did not achieve
> the convergence. Thus, from time 20.16 and on, your results are not
> reliable.
thats bad it failed