ConvergenceTest object failed

[danile]

hi all!

i built a rc frame 3d building (6 storeys) with fiber section. now i apply gravity load, but opensees does error:
"
WARNING: CTestEnergyIncr::test() - failed to converge
after: 10 iterations
NewtnRaphson::solveCurrentStep() -the ConvergenceTest object failed in test()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 0.1
OpenSees > analyze failed, returned: -3 error flag"


it's strange because before building a big rc frame i did some tests on a simple rc 3d frame; and i have no problem.

After i try with elasticbeamcolumn, but i have the same problem...


here down same steps of the script:

1) sections are like that:
"
set ColSecTag91_911 1;
set coverCol 30;
set numBars91_911 3;
set barArea91_911 314;
set numBarsInt91_911 1;
set barAreaInt91_911 314;
# RC section:
set cover40Y [expr $HCol40/2.0];
set cover40Z [expr $BCol40/2.0];
set core40Y [expr $cover40Y-$coverCol]
set core40Z [expr $cover40Z-$coverCol]
set nf40Y 17;
set nf40Z 17;
set nfcover40Y 1;
set nfcover40Z 1;
section fiberSec $ColSecTag91_911 {;
patch rect $IDconcC $nf40Y $nf40Z $core40Y $core40Z -$core40Y -$core40Z;
patch rect $IDconcU $nfcover40Y $nfcover40Z $cover40Y $cover40Z -$cover40Y $core40Z
patch rect $IDconcU $nfcover40Y $nfcover40Z $cover40Y -$core40Z -$cover40Y -$cover40Z
patch rect $IDconcU $nfcover40Y $nfcover40Z $cover40Y $core40Z $core40Y -$core40Z
patch rect $IDconcU $nfcover40Y $nfcover40Z -$core40Y $core40Z -$cover40Y -$core40Z
layer straight $IDreinf $numBars91_911 $barArea91_911 -$core40Y $core40Z -$core40Y -$core40Z;
layer straight $IDreinf $numBars91_911 $barArea91_911 $core40Y $core40Z $core40Y -$core40Z;
layer straight $IDreinf $numBarsInt91_911 $barAreaInt91_911 $core40Y $core40Z -$core40Y $core40Z;
layer straight $IDreinf $numBarsInt91_911 $barAreaInt91_911 $core40Y -$core40Z -$core40Y -$core40Z;
};

set SecTagColTorsion 1001;
set ColSecTag3D91_911 1101;
uniaxialMaterial Elastic $SecTagColTorsion $GJCol40;
section Aggregator $ColSecTag3D91_911 $SecTagColTorsion T -section $ColSecTag91_911;
set ColTransfTag91_911 1;
geomTransf Linear $ColTransfTag91_911 -1 0 0;"


2) beam and column are like that:
"
element nonlinearBeamColumn 91911 91 911 $numIntgrPtsColBase $ColSecTag3D91_911 $ColTransfTag91_911;"


3) gravity is applied on the nodes:

"set Q1 46 ;
set Q2 139 ;
set Q3 93 ;
pattern Plain 1 "Linear" {
load 11 12 21 22 31 32 41 42 51 52 61 62 0 0 -Q1 0 0 0
load 13 14 23 24 33 34 43 44 53 54 63 64 0 0 -Q2 0 0 0
load 15 16 25 26 35 36 45 46 55 56 65 66 0 0 -Q3 0 0 0
}

constraints Transformation;
numberer RCM;
system BandGeneral;
set Tol 1.0e-8;
test EnergyIncr $Tol 10 ;
algorithm Newton;
set NstepGravity 10;
set DGravity [expr 1./$NstepGravity];
integrator LoadControl $DGravity;
analysis Static;
analyze $NstepGravity;
loadConst -time 0.0"



However if i change section and i use elasticbeamcolumn there is the same problem.

could anyone tell me how can i solve this problem ?!?!?!?


thx

[silvia]

try to check the eigenvalues right after building the model.
if you are having problems with gravity, then yes, something's wrong with the model.
if you are having problems during lateral analysis, you may be okey, just use the scripts given in the examples manual -- i think example 5 and greater have them.

[danile]

hi silvia. thank you very much for the fast reply.

ok, eigenvalues are not right - eigenvalues negative!?!? impossible? -
(i found it before appling gravity) and are opensees says to me:
"
-1.777160e-010
-4.939411e-011
1.647267e-011
5.835856e-010
1.024824e-009
1.485272e+000"
domain error: argument not in valid range"

and a similar result is for eigenvalues analysis for elastic frames:

so i have to think that there are some problem, and in a part of the script common from elastic frame and non-elastic frame: node definition, fixing factors, rigidDiaphragm, geometric definition or transftag for 3d elements. not other things. it's true?


thx

[danile]

i found the error: fix

thx

[Jean Domínguez]

[quote="danile"]i found the error: fix

thx[/quote]

How do you fix the "domain not in valid range" error ?