eigen error

wangjingze1985 · Post by **wangjingze1985** » Mon Sep 06, 2010 1:03 am

I need to ask for the eigenvalues of a modle with a beam and two Bearings .But opensees give some warnings:
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 1
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 1
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 1
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 1
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 1
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 1
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 1
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 1

My code is:

# SET UP ----------------------------------------------------------------------------
wipe;
model basic -ndm 3 -ndf 6;
logFile errorFile;

node 1 -175 -3 0;
node 2 -85 -3 10;
node 3 0 -3 0;
node 4 85 -3 10;
node 5 175 -3 0;

node 11 -175 -3 0;
node 15 175 -3 0;

# Single point constraints -- Boundary Conditions

fix 11 1 1 1 1 1 1;

fix 15 1 1 1 1 1 1;

fix 2 0 0 0 1 1 1;
fix 3 0 0 0 1 1 1;
fix 4 0 0 0 1 1 1;

mass 1 127375 127375 127375 0 0 0;
mass 2 123325 123325 123325 0 0 0;
mass 3 123700 123700 123700 0 0 0;
mass 4 127375 127375 127375 0 0 0;
mass 5 123325 123325 123325 0 0 0;

# Define ELEMENTS -------------------------------------------------------------
#坐标转换
geomTransf Linear 1 0 -1 0; #zhu liang

#支座定义
#材料选取
uniaxialMaterial Elastic 11 10e+11;
uniaxialMaterial Elastic 12 2.750e+6;
uniaxialMaterial Elastic 13 1.119e+5;

element elasticBeamColumn 1 1 5 5.2 3.45e10 1.0e20 20 48.1 11.56 1;

element zeroLength 5 11 1 -mat 12 12 11 11 11 13 -dir 1 2 3 4 5 6;

element zeroLength 7 15 5 -mat 12 12 11 11 11 13 -dir 1 2 3 4 5 6;

constraints Plain;

numberer RCM;

system BandGeneral;

puts "[eigen 2]";

I think my modle may not be fully constrained,however,I do not know where the error is.Anyone can help me?

thanks!!

vesna · Post by **vesna** » Fri Sep 10, 2010 3:56 pm

Nodes 2, 3, and 4 are hanging nodes. Either delete them or add elements.

wangjingze1985 · Post by **wangjingze1985** » Wed Sep 15, 2010 7:01 pm

thank you very much four your reply. I have sloved the problem.

ksecond · Post by **ksecond** » Thu Mar 22, 2012 5:49 am

wangjingze1985 wrote:
> thank you very much four your reply. I have sloved the problem.
Hello,I came across such problem right now. I'm wondering how can you solve such problem? Thanks a lot ~

kadirtopcu · Post by **kadirtopcu** » Thu Nov 21, 2013 7:24 am

yes, how did you solve that problem?

brag006 · Post by **brag006** » Thu Nov 21, 2013 11:14 am

check you have no nodes hanging. As in nodes with nothing attached to them. Also check your constraints.

khabatfaramarzi · Post by **khabatfaramarzi** » Thu Jan 22, 2015 2:03 pm

please help me for solving this error

OpenSees -- Open System For Earthquake Engineering Simulation
Pacific Earthquake Engineering Research Center -- 2.4.4 (rev 5764)
Modified Ibarra-Medina-Krawinkler Model with Bilinear Hysteretic Response
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 49
WARNING ModifiedNewton::solveCurrentStep() -the LinearSysOfEqn failed in solve()

StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 0.1
OpenSees > analyze failed, returned: -3 error flag
Model Built
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 49
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 49
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 49
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 49
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 49
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 49
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 49
---------------------------------------------------------------------------
End of script <E:\xabat\arshad\final> reached, Press any key to continue
my code
for {set i 1} {$i<10} {incr i 1} {
file mkdir northridge/[expr $i./10]pga
}
set txt [open eigen/eigen.out w+]
model basic -ndm 2 -ndf 3
set EL 2.059E+11
set epsi [expr 0.0]
node 10 0 0
node 11 0 4.2
node 12 0 7.2
node 13 0 10.2
node 14 0 13.2
node 15 0 16.2
node 16 0 19.2
node 17 0 22.2
node 18 0 25.2
node 20 5 0
node 21 5 4.2
node 22 5 7.2
node 23 5 10.2
node 24 5 13.2
node 25 5 16.2
node 26 5 19.2
node 27 5 22.2
node 28 5 25.2
node 30 9 0
node 31 9 4.2
node 32 9 7.2
node 33 9 10.2
node 34 9 13.2
node 35 9 16.2
node 36 9 19.2
node 37 9 22.2
node 38 9 25.2
node 40 14 0
node 41 14 4.2
node 42 14 7.2
node 43 14 10.2
node 44 14 13.2
node 45 14 16.2
node 46 14 19.2
node 47 14 22.2
node 48 14 25.2
fix 10 1 1 1
fix 20 1 1 1
fix 30 1 1 1
fix 40 1 1 1
# mass define-----------------
# mass node st1-----------------
mass 11 8850. 0. 0.
mass 21 15930. 0. 0.
mass 31 15930. 0. 0.
mass 41 8850. 0. 0.
# mass node st2-----------------
mass 12 8850. 0. 0.
mass 22 15930. 0. 0.
mass 32 15930. 0. 0.
mass 42 8850. 0. 0.
# mass node st3-----------------
mass 13 8850. 0. 0.
mass 23 15930. 0. 0.
mass 33 15930. 0. 0.
mass 43 8850. 0. 0.
# mass node st4-----------------
mass 14 8850. 0. 0.
mass 24 15930. 0. 0.
mass 34 15930. 0. 0.
mass 44 8850. 0. 0.
# mass node st5-----------------
mass 15 8850. 0. 0.
mass 25 15930. 0. 0.
mass 35 15930. 0. 0.
mass 45 8850. 0. 0.
# mass node st6-----------------
mass 16 8850. 0. 0.
mass 26 15930. 0. 0.
mass 36 15930. 0. 0.
mass 46 8850. 0. 0.
# mass node st7-----------------
mass 17 8850. 0. 0.
mass 27 15930. 0. 0.
mass 37 15930. 0. 0.
mass 47 8850. 0. 0.
# mass node st8-----------------
mass 18 7950. 0. 0.
mass 28 14310. 0. 0.
mass 38 14310. 0. 0.
mass 48 7950 0. 0.
geomTransf Linear 1
# elsatic section---------------------------------
# SECTION OF COLUMN ST1
element elasticBeamColumn 1 10 11 3.352E-02 $EL 3.206220E-04 1
element elasticBeamColumn 2 20 21 3.675E-02 $EL 4.561E-04 1
element elasticBeamColumn 3 30 31 3.675E-02 $EL 4.561E-04 1
element elasticBeamColumn 4 40 41 3.352E-02 $EL 3.206220E-04 1
# SECTION OF BEAM ST1
element elasticBeamColumn 5 11 21 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 6 21 31 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 7 31 41 2.626E-02 $EL 2.289261E-04 1
# SECTION OF COLUMN ST2
element elasticBeamColumn 8 11 12 3.352E-02 $EL 3.206220E-04 1
element elasticBeamColumn 9 21 22 3.352E-02 $EL 3.206220E-04 1
element elasticBeamColumn 10 31 32 3.352E-02 $EL 3.206220E-04 1
element elasticBeamColumn 11 41 42 3.352E-02 $EL 3.206220E-04 1
# SECTION OF BEAM ST2
element elasticBeamColumn 12 12 22 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 13 22 32 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 14 32 42 2.626E-02 $EL 2.289261E-04 1
# SECTION OF COLUMN ST3
element elasticBeamColumn 15 12 13 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 16 22 23 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 17 32 33 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 18 42 43 2.626E-02 $EL 2.289261E-04 1
# SECTION OF BEAM ST3
element elasticBeamColumn 19 13 23 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 20 23 33 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 21 33 43 2.626E-02 $EL 2.289261E-04 1
# SECTION OF COLUMN ST4
element elasticBeamColumn 22 13 14 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 23 23 24 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 24 33 34 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 25 43 44 2.626E-02 $EL 2.289261E-04 1
# SECTION OF BEAM ST4
element elasticBeamColumn 26 14 24 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 27 24 34 2.626E-02 $EL 2.289261E-04 1
element elasticBeamColumn 28 34 44 2.626E-02 $EL 2.289261E-04 1
# SECTION OF COLUMN ST6
element elasticBeamColumn 29 14 15 2.200E-02 $EL 1.831333E-04 1
element elasticBeamColumn 30 24 25 2.200E-02 $EL 1.831333E-04 1
element elasticBeamColumn 31 34 35 2.200E-02 $EL 1.831333E-04 1
element elasticBeamColumn 32 44 45 2.200E-02 $EL 1.831333E-04 1
# SECTION OF BEAM ST5
element elasticBeamColumn 33 15 25 2.200E-02 $EL 1.831333E-04 1
element elasticBeamColumn 34 25 35 2.200E-02 $EL 1.831333E-04 1
element elasticBeamColumn 35 35 45 2.200E-02 $EL 1.831333E-04 1
# SECTION OF COLUMN ST6
element elasticBeamColumn 36 15 16 1.958E-02 $EL 1.430673E-04 1
element elasticBeamColumn 37 25 26 1.958E-02 $EL 1.430673E-04 1
element elasticBeamColumn 38 35 36 1.958E-02 $EL 1.430673E-04 1
element elasticBeamColumn 39 45 46 1.958E-02 $EL 1.430673E-04 1
# SECTION OF BEAM ST6
element elasticBeamColumn 40 16 26 1.958E-02 $EL 1.430673E-04 1
element elasticBeamColumn 41 26 36 1.958E-02 $EL 1.430673E-04 1
element elasticBeamColumn 42 36 46 1.958E-02 $EL 1.430673E-04 1
# SECTION OF COLUMN ST7
element elasticBeamColumn 43 16 17 1.900E-02 $EL 1.199833E-04 1
element elasticBeamColumn 44 26 27 1.900E-02 $EL 1.199833E-04 1
element elasticBeamColumn 45 36 37 1.900E-02 $EL 1.199833E-04 1
element elasticBeamColumn 46 46 47 1.900E-02 $EL 1.199833E-04 1
# SECTION OF BEAM ST7
element elasticBeamColumn 47 17 27 1.900E-02 $EL 1.199833E-04 1
element elasticBeamColumn 48 27 37 1.900E-02 $EL 1.199833E-04 1
element elasticBeamColumn 49 37 47 1.900E-02 $EL 1.199833E-04 1
# SECTION OF COLUMN ST8
element elasticBeamColumn 50 17 18 1.476E-02 $EL 8.107425E-05 1
element elasticBeamColumn 51 27 28 1.476E-02 $EL 8.107425E-05 1
element elasticBeamColumn 52 37 38 1.476E-02 $EL 8.107425E-05 1
element elasticBeamColumn 53 47 48 1.476E-02 $EL 8.107425E-05 1
# SECTION OF BEAM ST8
element elasticBeamColumn 54 18 28 1.476E-02 $EL 8.107425E-05 1
element elasticBeamColumn 55 28 38 1.476E-02 $EL 8.107425E-05 1
element elasticBeamColumn 56 38 48 1.476E-02 $EL 8.107425E-05 1

equalDOF 11 21 2 3
equalDOF 11 31 2 3
equalDOF 11 41 2 3
equalDOF 12 22 2 3
equalDOF 12 32 2 3
equalDOF 12 42 2 3
equalDOF 13 23 2 3
equalDOF 13 33 2 3
equalDOF 13 43 2 3
equalDOF 14 24 2 3
equalDOF 14 34 2 3
equalDOF 14 44 2 3
equalDOF 15 25 2 3
equalDOF 15 35 2 3
equalDOF 15 45 2 3
equalDOF 16 26 2 3
equalDOF 16 36 2 3
equalDOF 16 46 2 3
equalDOF 17 27 2 3
equalDOF 17 37 2 3
equalDOF 17 47 2 3
equalDOF 18 28 2 3
equalDOF 18 38 2 3
equalDOF 18 48 2 3
set Pi [expr 4*atan(1)]
set eig [eigen 8]
puts $txt "Eigen Values are: $eig"
puts $txt "1 Mode Period is:[expr 2*3.1415/pow([lindex $eig 0],0.5)]"
puts $txt "2 Mode Period is:[expr 2*3.1415/pow([lindex $eig 1],0.5)]"
puts $txt "3 Mode Period is:[expr 2*3.1415/pow([lindex $eig 2],0.5)]"
puts $txt "4 Mode Period is:[expr 2*3.1415/pow([lindex $eig 3],0.5)]"
puts $txt "5 Mode Period is:[expr 2*3.1415/pow([lindex $eig 4],0.5)]"
puts $txt "6 Mode Period is:[expr 2*3.1415/pow([lindex $eig 5],0.5)]"
puts $txt "7 Mode Period is:[expr 2*3.1415/pow([lindex $eig 6],0.5)]"
puts $txt "8 Mode Period is:[expr 2*3.1415/pow([lindex $eig 7],0.5)]"
# column hinges at BASE
node 101 0.0 0.0
node 201 5.0 0.0
node 301 9.0 0.0
node 401 14.0 0.0
# column hinges at bot ST1
node 113 0.0 4.2
node 213 5.0 4.2
node 313 9.0 4.2
node 413 14.0 4.2
# column hinges at top ST1
node 111 0.0 4.2
node 211 5.0 4.2
node 311 9.0 4.2
node 411 14.0 4.2
# column hinges at bot ST2
node 123 0.0 7.2
node 223 5.0 7.2
node 323 9.0 7.2
node 423 14.0 7.2
# column hinges at top ST2
node 121 0.0 7.2
node 221 5.0 7.2
node 321 9.0 7.2
node 421 14.0 7.2
# column hinges at bot ST3
node 133 0.0 10.2
node 233 5.0 10.2
node 333 9.0 10.2
node 433 14.0 10.2
# column hinges at top ST3
node 131 0.0 10.2
node 231 5.0 10.2
node 331 9.0 10.2
node 431 14.0 10.2
# column hinges at bot ST4
node 143 0.0 13.2
node 243 5.0 13.2
node 343 9.0 13.2
node 443 14.0 13.2
# column hinges at top ST4
node 141 0.0 13.2
node 241 5.0 13.2
node 341 9.0 13.2
node 441 14.0 13.2
# column hinges at bot ST5
node 153 0.0 16.2
node 253 5.0 16.2
node 353 9.0 16.2
node 453 14.0 16.2
# column hinges at top ST5
node 151 0.0 16.2
node 251 5.0 16.2
node 351 9.0 16.2
node 451 14.0 16.2
# column hinges at bot ST6
node 163 0.0 19.2
node 263 5.0 19.2
node 363 9.0 19.2
node 463 14.0 19.2
# column hinges at top ST6
node 161 0.0 19.2
node 261 5.0 19.2
node 361 9.0 19.2
node 461 14.0 19.2
# column hinges at bot ST7
node 173 0.0 22.2
node 273 5.0 22.2
node 373 9.0 22.2
node 473 14.0 22.2
# column hinges at top ST7
node 171 0.0 22.2
node 271 5.0 22.2
node 371 9.0 22.2
node 471 14.0 22.2
# column hinges at bot ST8
node 183 0.0 25.2
node 283 5.0 25.2
node 383 9.0 25.2
node 483 14.0 25.2
# beam hinges at ST1
node 112 [expr 0 + $epsi] 4.2
node 214 [expr 5 - $epsi] 4.2
node 212 [expr 5 + $epsi] 4.2
node 314 [expr 9 - $epsi] 4.2
node 312 [expr 9 + $epsi] 4.2
node 414 [expr 14 - $epsi] 4.2
# beam hinges at ST2
node 122 [expr 0 + $epsi] 7.2
node 224 [expr 5 - $epsi] 7.2
node 222 [expr 5 + $epsi] 7.2
node 324 [expr 9 - $epsi] 7.2
node 322 [expr 9 + $epsi] 7.2
node 424 [expr 14 - $epsi] 7.2
# beam hinges at ST3
node 132 [expr 0 + $epsi] 10.2
node 234 [expr 5 - $epsi] 10.2
node 232 [expr 5 + $epsi] 10.2
node 334 [expr 9 - $epsi] 10.2
node 332 [expr 9 + $epsi] 10.2
node 434 [expr 14 - $epsi] 10.2
# beam hinges at ST4
node 142 [expr 0 + $epsi] 13.2
node 244 [expr 5 - $epsi] 13.2
node 242 [expr 5 + $epsi] 13.2
node 344 [expr 9 - $epsi] 13.2
node 342 [expr 9 + $epsi] 13.2
node 444 [expr 14 - $epsi] 13.2
# beam hinges at ST5
node 152 [expr 0 + $epsi] 16.2
node 254 [expr 5 - $epsi] 16.2
node 252 [expr 5 + $epsi] 16.2
node 354 [expr 9 - $epsi] 16.2
node 352 [expr 9 + $epsi] 16.2
node 454 [expr 14 - $epsi] 16.2
# beam hinges at ST6
node 162 [expr 0 + $epsi] 19.2
node 264 [expr 5 - $epsi] 19.2
node 262 [expr 5 + $epsi] 19.2
node 364 [expr 9 + $epsi] 19.2
node 362 [expr 9 + $epsi] 19.2
node 464 [expr 14 - $epsi] 19.2
# beam hinges at ST7
node 172 [expr 0 + $epsi] 22.2
node 274 [expr 5 - $epsi] 22.2
node 272 [expr 5 + $epsi] 22.2
node 374 [expr 9 - $epsi] 22.2
node 372 [expr 9 + $epsi] 22.2
node 474 [expr 14 - $epsi] 22.2
# beam hinges at ST8
node 182 [expr 0 + $epsi] 25.2
node 284 [expr 5 - $epsi] 25.2
node 282 [expr 5 + $epsi] 25.2
node 384 [expr 9 - $epsi] 25.2
node 382 [expr 9 + $epsi] 25.2
node 484 [expr 14 - $epsi] 25.2
#--------------------------------------------------------------------
# uniaxialMaterial Bilin $matTag $K0 $as_Plus $as_Neg $My_Plus $My_Neg $Lamda_S $Lamda_C $Lamda_A $Lamda_K $c_S $c_C $c_A $c_K $theta_p_Plus $theta_p_Neg $theta_pc_Plus $theta_pc_Neg $Res_Pos $Res_Neg $theta_u_Plus $theta_u_Neg $D_Plus $D_Neg <$nFactor>
uniaxialMaterial Bilin 1 10692593.10 2.65E-02 2.65E-02 1829241 -1829241 5.877 5.877 5.877 5.877 1 1 1 1 0.059 0.059 0.343 0.343 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 2 15211814.63 2.65E-02 2.65E-02 2307155 -2307155 3.917 3.917 3.917 3.917 1 1 1 1 0.047 0.047 0.262 0.262 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 3 10692593.10 2.65E-02 2.65E-02 1829241 -1829241 5.877 5.877 5.877 5.877 1 1 1 1 0.059 0.059 0.343 0.343 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 4 15211814.63 2.65E-02 2.65E-02 2307155 -2307155 3.917 3.917 3.917 3.917 1 1 1 1 0.047 0.047 0.262 0.262 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 5 9427176.80 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.06 0.06 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 6 11783971.00 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.056 0.056 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 7 9427176.80 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.06 0.06 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 8 10692593.1 2.65E-02 2.65E-02 1829241 -1829241 5.877 5.877 5.877 5.877 1 1 1 1 0.053 0.053 0.343 0.343 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 9 10692593.1 2.65E-02 2.65E-02 1829241 -1829241 5.877 5.877 5.877 5.877 1 1 1 1 0.053 0.053 0.343 0.343 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 10 10692593.1 2.65E-02 2.65E-02 1829241 -1829241 5.877 5.877 5.877 5.877 1 1 1 1 0.053 0.053 0.343 0.343 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 11 10692593.1 2.65E-02 2.65E-02 1829241 -1829241 5.877 5.877 5.877 5.877 1 1 1 1 0.053 0.053 0.343 0.343 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 12 9427176.80 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.06 0.06 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 13 11783971.00 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.056 0.056 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 14 9427176.80 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.06 0.06 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 15 20949281.77 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.051 0.051 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 16 20949281.77 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.051 0.051 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 17 20949281.77 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.051 0.051 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 18 20949281.77 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.051 0.051 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 19 9427176.80 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.06 0.06 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 20 11783971.00 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.056 0.056 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 21 9427176.80 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.06 0.06 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 22 20949281.77 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.051 0.051 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 23 20949281.77 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.051 0.051 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 24 20949281.77 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.051 0.051 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 25 20949281.77 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.051 0.051 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 26 9427176.80 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.06 0.06 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 27 11783971.00 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.056 0.056 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 28 9427176.80 2.65E-02 2.65E-02 1383971 -1383971 4.096 4.096 4.096 4.096 1 1 1 1 0.06 0.06 0.278 0.278 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 29 16758731.76 2.65E-02 2.65E-02 1154863 -1154863 3.033 3.033 3.033 3.033 1 1 1 1 0.047 0.047 0.217 0.217 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 30 16758731.76 2.65E-02 2.65E-02 1154863 -1154863 3.033 3.033 3.033 3.033 1 1 1 1 0.047 0.047 0.217 0.217 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 31 16758731.76 2.65E-02 2.65E-02 1154863 -1154863 3.033 3.033 3.033 3.033 1 1 1 1 0.047 0.047 0.217 0.217 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 32 16758731.76 2.65E-02 2.65E-02 1154863 -1154863 3.033 3.033 3.033 3.033 1 1 1 1 0.047 0.047 0.217 0.217 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 33 7541429.29 2.65E-02 2.65E-02 1154863 -1154863 3.033 3.033 3.033 3.033 1 1 1 1 0.056 0.056 0.217 0.217 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 34 9426786.62 2.65E-02 2.65E-02 1154863 -1154863 3.033 3.033 3.033 3.033 1 1 1 1 0.056 0.056 0.217 0.217 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 35 7541429.29 2.65E-02 2.65E-02 1154863 -1154863 3.033 3.033 3.033 3.033 1 1 1 1 0.052 0.052 0.217 0.217 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 36 13092247.59 2.65E-02 2.65E-02 955689 -955689 3.130 3.130 3.130 3.130 1 1 1 1 0.052 0.052 0.235 0.235 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 37 13092247.59 2.65E-02 2.65E-02 955689 -955689 3.130 3.130 3.130 3.130 1 1 1 1 0.052 0.052 0.235 0.235 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 38 13092247.59 2.65E-02 2.65E-02 955689 -955689 3.130 3.130 3.130 3.130 1 1 1 1 0.052 0.052 0.235 0.235 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 39 13092247.59 2.65E-02 2.65E-02 955689 -955689 3.130 3.130 3.130 3.130 1 1 1 1 0.052 0.052 0.235 0.235 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 40 5891511.41 2.65E-02 2.65E-02 955689 -955689 3.130 3.130 3.130 3.130 1 1 1 1 0.061 0.061 0.235 0.235 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 41 7364389.27 2.65E-02 2.65E-02 955689 -955689 3.130 3.130 3.130 3.130 1 1 1 1 0.057 0.057 0.235 0.235 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 42 5891511.41 2.65E-02 2.65E-02 955689 -955689 3.130 3.130 3.130 3.130 1 1 1 1 0.061 0.061 0.235 0.235 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 43 10979805.10 2.65E-02 2.65E-02 859526 -859526 4.082 4.082 4.082 4.082 1 1 1 1 0.059 0.059 0.276 0.276 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 44 10979805.10 2.65E-02 2.65E-02 859526 -859526 4.082 4.082 4.082 4.082 1 1 1 1 0.059 0.059 0.276 0.276 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 45 10979805.10 2.65E-02 2.65E-02 859526 -859526 4.082 4.082 4.082 4.082 1 1 1 1 0.059 0.059 0.276 0.276 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 46 10979805.10 2.65E-02 2.65E-02 859526 -859526 4.082 4.082 4.082 4.082 1 1 1 1 0.059 0.059 0.276 0.276 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 47 4940912.29 2.65E-02 2.65E-02 859526 -859526 4.082 4.082 4.082 4.082 1 1 1 1 0.07 0.07 0.276 0.276 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 48 6176140.37 2.65E-02 2.65E-02 859526 -859526 4.082 4.082 4.082 4.082 1 1 1 1 0.065 0.065 0.276 0.276 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 49 4940912.29 2.65E-02 2.65E-02 859526 -859526 4.082 4.082 4.082 4.082 1 1 1 1 0.07 0.07 0.276 0.276 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 50 7419194.70 2.65E-02 2.65E-02 621744 -621744 3.13 3.13 3.13 3.13 1 1 1 1 0.061 0.061 0.252 0.252 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 51 7419194.70 2.65E-02 2.65E-02 621744 -621744 3.13 3.13 3.13 3.13 1 1 1 1 0.061 0.061 0.252 0.252 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 52 7419194.70 2.65E-02 2.65E-02 621744 -621744 3.13 3.13 3.13 3.13 1 1 1 1 0.061 0.061 0.252 0.252 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 53 7419194.70 2.65E-02 2.65E-02 621744 -621744 3.13 3.13 3.13 3.13 1 1 1 1 0.061 0.061 0.252 0.252 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 54 3338637.62 2.65E-02 2.65E-02 621744 -621744 3.13 3.13 3.13 3.13 1 1 1 1 0.072 0.072 0.252 0.252 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 55 4173297.02 2.65E-02 2.65E-02 621744 -621744 3.13 3.13 3.13 3.13 1 1 1 1 0.067 0.067 0.252 0.252 0 0 0.07 0.07 1 1
uniaxialMaterial Bilin 56 3338637.62 2.65E-02 2.65E-02 621744 -621744 3.13 3.13 3.13 3.13 1 1 1 1 0.072 0.072 0.252 0.252 0 0 0.07 0.07 1 1
# Column hinges ------------------------------
element zeroLength 101 10 101 -mat 1 -dir 6
element zeroLength 102 20 201 -mat 2 -dir 6
element zeroLength 103 30 301 -mat 3 -dir 6
element zeroLength 104 40 401 -mat 4 -dir 6
element zeroLength 105 11 113 -mat 1 -dir 6
element zeroLength 106 11 111 -mat 8 -dir 6
element zeroLength 107 21 213 -mat 2 -dir 6
element zeroLength 108 21 211 -mat 9 -dir 6
element zeroLength 109 31 313 -mat 3 -dir 6
element zeroLength 110 31 311 -mat 10 -dir 6
element zeroLength 111 41 413 -mat 4 -dir 6
element zeroLength 112 41 411 -mat 11 -dir 6
element zeroLength 113 12 123 -mat 8 -dir 6
element zeroLength 114 12 121 -mat 15 -dir 6
element zeroLength 115 22 223 -mat 9 -dir 6
element zeroLength 116 22 221 -mat 16 -dir 6
element zeroLength 117 32 323 -mat 10 -dir 6
element zeroLength 118 32 321 -mat 17 -dir 6
element zeroLength 119 42 423 -mat 11 -dir 6
element zeroLength 120 42 421 -mat 18 -dir 6
element zeroLength 121 13 133 -mat 15 -dir 6
element zeroLength 122 13 131 -mat 22 -dir 6
element zeroLength 123 23 233 -mat 16 -dir 6
element zeroLength 124 23 231 -mat 23 -dir 6
element zeroLength 125 33 333 -mat 17 -dir 6
element zeroLength 126 33 331 -mat 24 -dir 6
element zeroLength 127 43 433 -mat 18 -dir 6
element zeroLength 128 43 431 -mat 25 -dir 6
element zeroLength 129 14 143 -mat 22 -dir 6
element zeroLength 130 14 141 -mat 29 -dir 6
element zeroLength 131 24 243 -mat 23 -dir 6
element zeroLength 132 24 241 -mat 30 -dir 6
element zeroLength 133 34 343 -mat 24 -dir 6
element zeroLength 134 34 341 -mat 31 -dir 6
element zeroLength 135 44 443 -mat 25 -dir 6
element zeroLength 136 44 441 -mat 32 -dir 6
element zeroLength 137 15 153 -mat 29 -dir 6
element zeroLength 138 15 151 -mat 36 -dir 6
element zeroLength 139 25 253 -mat 30 -dir 6
element zeroLength 140 25 251 -mat 37 -dir 6
element zeroLength 141 35 353 -mat 31 -dir 6
element zeroLength 142 35 351 -mat 38 -dir 6
element zeroLength 143 45 453 -mat 32 -dir 6
element zeroLength 144 45 451 -mat 39 -dir 6
element zeroLength 145 16 163 -mat 36 -dir 6
element zeroLength 146 16 161 -mat 43 -dir 6
element zeroLength 147 26 263 -mat 37 -dir 6
element zeroLength 148 26 261 -mat 44 -dir 6
element zeroLength 149 36 363 -mat 38 -dir 6
element zeroLength 150 36 361 -mat 45 -dir 6
element zeroLength 151 46 463 -mat 39 -dir 6
element zeroLength 152 46 461 -mat 46 -dir 6
element zeroLength 153 17 173 -mat 43 -dir 6
element zeroLength 154 17 171 -mat 50 -dir 6
element zeroLength 155 27 273 -mat 44 -dir 6
element zeroLength 156 27 271 -mat 51 -dir 6
element zeroLength 157 37 373 -mat 45 -dir 6
element zeroLength 158 37 371 -mat 52 -dir 6
element zeroLength 159 47 473 -mat 46 -dir 6
element zeroLength 160 47 471 -mat 53 -dir 6
element zeroLength 161 18 183 -mat 50 -dir 6
element zeroLength 162 28 283 -mat 51 -dir 6
element zeroLength 163 38 383 -mat 52 -dir 6
element zeroLength 164 48 483 -mat 53 -dir 6
# Beam hinges ------------------------------
element zeroLength 201 11 112 -mat 5 -dir 6
element zeroLength 202 21 214 -mat 5 -dir 6
element zeroLength 203 21 212 -mat 6 -dir 6
element zeroLength 204 31 314 -mat 6 -dir 6
element zeroLength 205 31 312 -mat 7 -dir 6
element zeroLength 206 41 414 -mat 7 -dir 6
element zeroLength 207 12 122 -mat 12 -dir 6
element zeroLength 208 22 224 -mat 12 -dir 6
element zeroLength 209 22 222 -mat 13 -dir 6
element zeroLength 210 32 324 -mat 13 -dir 6
element zeroLength 211 32 322 -mat 14 -dir 6
element zeroLength 212 42 424 -mat 14 -dir 6
element zeroLength 213 13 132 -mat 19 -dir 6
element zeroLength 214 23 234 -mat 19 -dir 6
element zeroLength 215 23 232 -mat 20 -dir 6
element zeroLength 216 33 334 -mat 20 -dir 6
element zeroLength 217 33 332 -mat 21 -dir 6
element zeroLength 218 43 434 -mat 21 -dir 6
element zeroLength 219 14 142 -mat 26 -dir 6
element zeroLength 220 24 244 -mat 26 -dir 6
element zeroLength 221 24 242 -mat 27 -dir 6
element zeroLength 222 34 344 -mat 27 -dir 6
element zeroLength 223 34 342 -mat 28 -dir 6
element zeroLength 224 44 444 -mat 28 -dir 6
element zeroLength 225 15 152 -mat 33 -dir 6
element zeroLength 226 25 254 -mat 33 -dir 6
element zeroLength 227 25 252 -mat 34 -dir 6
element zeroLength 228 35 354 -mat 34 -dir 6
element zeroLength 229 35 352 -mat 35 -dir 6
element zeroLength 230 45 454 -mat 35 -dir 6
element zeroLength 231 16 162 -mat 40 -dir 6
element zeroLength 232 26 264 -mat 40 -dir 6
element zeroLength 233 26 262 -mat 41 -dir 6
element zeroLength 234 36 364 -mat 41 -dir 6
element zeroLength 235 36 362 -mat 42 -dir 6
element zeroLength 236 46 464 -mat 42 -dir 6
element zeroLength 237 17 172 -mat 47 -dir 6
element zeroLength 238 27 274 -mat 47 -dir 6
element zeroLength 239 27 272 -mat 48 -dir 6
element zeroLength 240 37 374 -mat 48 -dir 6
element zeroLength 241 37 372 -mat 49 -dir 6
element zeroLength 242 47 474 -mat 49 -dir 6
element zeroLength 243 18 182 -mat 54 -dir 6
element zeroLength 244 28 284 -mat 54 -dir 6
element zeroLength 245 28 282 -mat 55 -dir 6
element zeroLength 246 38 384 -mat 55 -dir 6
element zeroLength 247 38 382 -mat 56 -dir 6
element zeroLength 248 48 484 -mat 56 -dir 6
# Eigenvalue analysis----------------------------
# Uniform loading ST1 to ST7 dead and live
pattern Plain 1 Constant {
eleLoad -ele 5 6 7 12 13 14 19 20 21 26 27 28 33 34 35 40 41 42 47 48 49 -type -beamUniform -3540 0.0;
}
# Uniform loading ST8 dead and live
pattern Plain 2 Constant {
eleLoad -ele 54 55 56 -type -beamUniform -3180 0.0;
}
#define gravity loads and perform gravity analysis

constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.e-8 60
algorithm ModifiedNewton
integrator LoadControl 0.1
analysis Static
analyze 10
loadConst -time 0.0
# maintain constant gravity loads and reset time to zero
puts "Model Built"
recorder Node -file nortridge/[expr $i./10]pga/node11dispX[expr $i./10]disp11.out -time -node 11 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node21dispX[expr $i./10]disp21.out -time -node 21 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node31dispX[expr $i./10]disp31.out -time -node 31 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node41dispX[expr $i./10]disp41.out -time -node 41 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node12dispX[expr $i./10]disp12.out -time -node 12 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node22dispX[expr $i./10]disp22.out -time -node 22 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node32dispX[expr $i./10]disp32.out -time -node 32 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node42dispX[expr $i./10]disp42.out -time -node 42 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node13dispX[expr $i./10]disp13.out -time -node 13 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node23dispX[expr $i./10]disp23.out -time -node 23 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node33dispX[expr $i./10]disp33.out -time -node 33 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node43dispX[expr $i./10]disp43.out -time -node 43 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node14dispX[expr $i./10]disp14.out -time -node 14 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node24dispX[expr $i./10]disp24.out -time -node 24 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node34dispX[expr $i./10]disp34.out -time -node 34 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node44dispX[expr $i./10]disp44.out -time -node 44 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node15dispX[expr $i./10]disp15.out -time -node 15 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node25dispX[expr $i./10]disp25.out -time -node 25 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node35dispX[expr $i./10]disp35.out -time -node 35 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node45dispX[expr $i./10]disp45.out -time -node 45 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node16dispX[expr $i./10]disp16.out -time -node 16 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node26dispX[expr $i./10]disp26.out -time -node 26 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node36dispX[expr $i./10]disp36.out -time -node 36 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node46dispX[expr $i./10]disp46.out -time -node 46 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node17dispX[expr $i./10]disp17.out -time -node 17 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node27dispX[expr $i./10]disp27.out -time -node 27 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node37dispX[expr $i./10]disp37.out -time -node 37 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node47dispX[expr $i./10]disp47.out -time -node 47 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node18dispX[expr $i./10]disp18.out -time -node 18 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node28dispX[expr $i./10]disp28.out -time -node 28 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node38dispX[expr $i./10]disp38.out -time -node 38 -dof 1 disp
recorder Node -file nortridge/[expr $i./10]pga/node48dispX[expr $i./10]disp48.out -time -node 48 -dof 1 disp
for {set j 1} {$j<=15} {incr j} {

# define damping
set xDamp 0.05; # damping ratio
set MpropSwitch 1.0;
set KcurrSwitch 0.0;
set KcommSwitch 1.0;
set KinitSwitch 0.0;
set nEigenI 1; # mode 1
set nEigenJ 3; # mode 3
set lambdaN [eigen [expr $nEigenJ]]; # eigenvalue analysis for nEigenJ modes

set lambdaI [lindex $lambdaN [expr $nEigenI-1]]; # eigenvalue mode i
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]]; # eigenvalue mode j
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM [expr $MpropSwitch*$xDamp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)]; # M-prop. damping; D = alphaM*M
set betaKcurr [expr $KcurrSwitch*2.*$xDamp/($omegaI+$omegaJ)]; # current-K; +beatKcurr*KCurrent
set betaKcomm [expr $KcommSwitch*2.*$xDamp/($omegaI+$omegaJ)]; # last-committed K; +betaKcomm*KlastCommitt
set betaKinit [expr $KinitSwitch*2.*$xDamp/($omegaI+$omegaJ)]; # initial-K; +beatKinit*Kini
rayleigh $alphaM $betaKcurr $betaKinit $betaKcomm; # RAYLEIGH damping

puts "damping ok!"
#-------------------------------
# change these parameters for each record:
# TmaxAnalysis , DtAnalysis , "Series -dt .. - filepath .."
set TmaxAnalysis 59.98 ; # maximum duration of ground-motion analysis
set DtAnalysis 0.02; # time-step Dt for lateral analysis
#increment gmf
set accel "Series -dt 0.02 -filePath northridge.txt -factor [expr $i.*9.81*12.52/10]"

pattern UniformExcitation [expr $i+100] 1 -accel $accel

#------Dynamic Analysis Parameters
variable constraintsTypeDynamic Transformation;
constraints $constraintsTypeDynamic ;

variable numbererTypeDynamic RCM
numberer $numbererTypeDynamic

variable systemTypeDynamic UmfPack; # try UmfPack for large problems
system $systemTypeDynamic;
variable testTypeDynamic EnergyIncr; # Convergence-test type
# Define the ranges of tolerances to try for convergence. Note that the maximum tolerance
# used in the analsyis will be output in the output files.
test $testTypeDynamic $TolDynamic $maxNumIterDynamic $printFlagDynamic;
set testTolerance 1.0e-6
set testMinTolerance1 1.0e-5;
set testMinTolerance2 1.0e-4;
set testMinTolerance3 1.0e-3;
set testMinTolerance4 1.0e-2;
set testMinTolerance5 1.0e-1;

# Define the iteration information used for difference situations and tests
set testIterations 50; # Changed on 1-7-05 for Corotational transformation
set testInitialIterations 1000
set testLowIter 10; # Used to try each test in the loop
set ratioForInitialAlgo 200; # This is the ratio of testIterations that is allowed for -initial test
set testHighIter 1000; # Used to try to make it converge at the very end

# Define the initial test, this will be changed later during convergance work.
test $testType $testTolerance $testIterations 0
# Set initialize the maximum tolerance used - for output to know about what tolerance was used to obtain convergence.
set maxTolUsed $testTolerance;

variable algorithmTypeDynamic ModifiedNewton;
algorithm $algorithmTypeDynamic;

variable NewmarkGamma 0.5; # Newmark-integrator gamma parameter (also HHT)
variable NewmarkBeta 0.25; # Newmark-integrator beta parameter
variable integratorTypeDynamic Newmark;
integrator $integratorTypeDynamic $NewmarkGamma $NewmarkBeta
variable analysisTypeDynamic Transient
analysis $analysisTypeDynamic

set numIncrForsAnalysis 1; # This is the number of increments that are used before I enter my convergence loop
# BE SURE TO DO 1 so it will pop out and I can check the displ. for possible
# collapse and for singularity AT EVERY STEP (changed on 1-21-05).
set dt [expr $DtAnalysis];
set dtMin [expr $DtAnalysis / 100.0];
set dtMax [expr $DtAnalysis];

set convLogFilename [open convLogFileOUT.out w]
set convPlotFilename [open convPlotFileOUT.out w]

set Nsteps [expr int($TmaxAnalysis/$DtAnalysis)];
set ok [analyze $Nsteps $DtAnalysis]; # actually perform analysis; returns ok=0 if analysis was successful

if {$ok != 0} { ; # analysis was not successful.
# --------------------------------------------------------------------------------------------------
# change some analysis parameters to achieve convergence
# performance is slower inside this loop
# Time-controlled analysis
set ok 0;
set currentTime [getTime];
while {$currentTime < $TmaxAnalysis && $ok == 0} {
test $testType $testTolerance $testIterations 0

# Altered here for Cordova Frank
set ok [analyze 1 $DtAnalysis]
set ok [analyze $numIncrForsAnalysis $dt $dtMin $dtMax]

# Keep track of the maximum tolarance used in this step - for the convPlotFile. This is later increased it necessary.
set maxTolUsedInCurrentStep [expr $testTolerance]

#### Change things for convergence ###############
# If it's not ok, try to decrease dT, but keep the toerance the samecall another file for this (just to keep this file clean)
set currentTolerance [expr $testTolerance]
set currentNumIterations [expr $testLowIter]
set currentDt [expr $DtAnalysis/10];
source SolutionAlgorithmSubFile.tcl

# If it's not ok, try to decrease dT a bit more, but keep the toerance the samecall another file for this (just to keep this file clean)
set currentTolerance [expr $testTolerance]
set currentNumIterations [expr $testLowIter]
set currentDt [expr $DtAnalysis/20];
source SolutionAlgorithmSubFile.tcl

# If it's not ok, try to decrease dT a bit more, but keep the toerance the samecall another file for this (just to keep this file clean)
set currentTolerance [expr $testTolerance]
set currentNumIterations [expr $testLowIter]
set currentDt [expr $DtAnalysis/40];
source SolutionAlgorithmSubFile.tcl

# If it's not ok, try to decrease dT a bit more, but keep the toerance the samecall another file for this (just to keep this file clean)
set currentTolerance [expr $testTolerance]
set currentNumIterations [expr $testLowIter]
set currentDt [expr $DtAnalysis/80];
source SolutionAlgorithmSubFile.tcl

# If it's not ok, go to a more relaxed tolerance1...call another file for this (just to keep this file clean)
set currentTolerance [expr $testMinTolerance1]
set currentNumIterations [expr $testLowIter]
set currentDt [expr $DtAnalysis/10]; # This was /20, so maybe change back?
source SolutionAlgorithmSubFile.tcl

# If it's not ok, go to a more relaxed tolerance1...call another file for this (just to keep this file clean)
set currentTolerance [expr $testMinTolerance1]
set currentNumIterations [expr $testLowIter]
set currentDt [expr $DtAnalysis/20]; # This was /20, so maybe change back?
source SolutionAlgorithmSubFile.tcl

# If it's not ok, go to a more relaxed tolerance1...call another file for this (just to keep this file clean)
set currentTolerance [expr $testMinTolerance1]
set currentNumIterations [expr $testLowIter]
set currentDt [expr $DtAnalysis/40]; # This was /20, so maybe change back?
source SolutionAlgorithmSubFile.tcl

# If it's not ok, go to a more relaxed tolerance1...call another file for this (just to keep this file clean)
set currentTolerance [expr $testMinTolerance1]
set currentNumIterations [expr $testLowIter]
set currentDt [expr $DtAnalysis/80]; # This was /20, so maybe change back?
source SolutionAlgorithmSubFile.tcl

# If it's not ok, go to a more relaxed tolerance2...call another file for this (just to keep this file clean)
set currentTolerance [expr $testMinTolerance2]
set currentNumIterations [expr $testLowIter]
set currentDt [expr $DtAnalysis/10];
source SolutionAlgorithmSubFile.tcl

# If it's not ok, go to a more relaxed tolerance3...call another file for this (just to keep this file clean)
# Decrease dT more
set currentTolerance [expr $testMinTolerance3]
set currentNumIterations [expr $testHighIter]
set currentDt [expr $DtAnalysis/20];
source SolutionAlgorithmSubFile.tcl

# If it's not ok, go to a more relaxed tolerance4...call another file for this (just to keep this file clean)
# Increase the number of iterations
set currentTolerance [expr $testMinTolerance4]
set currentNumIterations [expr $testHighIter]
set currentDt [expr $DtAnalysis/20];
source SolutionAlgorithmSubFile.tcl

# If it's not ok, go to a more relaxed tolerance5...call another file for this (just to keep this file clean)
set currentTolerance [expr $testMinTolerance5]
set currentNumIterations [expr $testHighIter]
set currentDt [expr $DtAnalysis/20];
source SolutionAlgorithmSubFile.tcl

# Write a line to the convPlotFile to give the maximumTolerance used at this step (to plot convIndex over time)
set convPlotFileLine "$currentTime $maxTolUsedInCurrentStep"
puts $convPlotFilename $convPlotFileLine

# Update the time
set currentTime [getTime]

}
}
puts "Ground Motion[expr $i./10]pga done!.End Time: [getTime]"
}

Post by **fmk** » Thu Jan 22, 2015 3:16 pm

all the nodes after the first eigen command need some restraint in translational x and y!

homashanehsaz · Post by **homashanehsaz** » Tue Feb 03, 2015 11:57 am

hi. sorry if I post here
actually i do not know where should i post my problem but my problem has the same error:
warning bandgenlinlapacksolver
and this is my file:
# units: N, m

# Remove existing model
wipe

# Create ModelBuilder (with two-dimensions and 2 DOF/node)
model BasicBuilder -ndm 2 -ndf3

# Create nodes
# ------------
# Create nodes & add to Domain - command: node nodeId xCrd yCrd
node 1 0.0 0.0
node 2 1.0 0.0

# Set the boundary conditions - command: fix nodeID xResrnt? yRestrnt?
fix 1 0 1 1
fix 2 1 1 1

# Define materials for truss elements
# -----------------------------------
# Create Elastic material prototype - command: uniaxialMaterial Elastic matID E
uniaxialMaterial ElasticPP 1 200000.0 0.0012

# create section for StressStrain diagram

section Uniaxial 1 1 P

# Define elements

# Create truss elements - command: element truss trussID node1 node2 A matID
element truss 1 1 2 1.0 1

# Define loads
# ------------
#

# create a Linear TimeSeries with a tag of 1
timeSeries Linear 1
# Create a Plain load pattern associated with the TimeSeries,
# command: pattern Plain $patternTag $timeSeriesTag { load commands }

pattern Plain 2 1 {
# Create the nodal load - command: load nodeID xForce yForce
load 1 1.0 0.0 0.0
}

# Create a recorder to monitor nodal displacements

# create a Recorder object for the nodal displacements at node 1
recorder Node -file disp1.out -time -node 1 -dof 1 disp

#recorder display force-didp 20 20 500 500

#recorder plot disp1.out time-didp 200 200 500 500 -columns 1 2

# create a Recorder for element forces, one for global system and the other for local system
recorder Element -file eleGlobal.out -time -ele 1 localforces
recorder Element -file eleLocal.out -time -ele 1 basicForces

# get StressStrain diagram
recorder Element -file ele1sec1StressStrain.out -time -ele 1 section 1 fiber .1 .1 stressStrain
# recorder Element -file ele1sec1StressStrain.out -time -ele 1 section 1 fiber x stressStrain
# recorder Element -file StressStrain.out -time ele 1 sec 1

## dasturate analys
# Create the constraint handler, a Plain handler is used as homo constraints
#constraints Plain
constraints Transformation

# Create the DOF numberer, the reverse Cuthill-McKee algorithm
numberer RCM
# dasture azmune hamgerayee safhe 85,
test NormDispIncr 1.0e-8 150
# Create the system of equation# systeme hale moadelat ra miguyad
#system BandSPD
system BandGeneral
#system SparseGeneral
#system SparseGeneral -piv
#system UmfPack
#system ProfileSPD

# Create the solution algorithm, a Linear algorithm is created
#algorithm Linear
algorithm Newton

set j 0.000002
for {set i 1} {$i<=40} {incr i 1} {
set j [expr -1*$j]
integrator DisplacementControl 1 1 $j
# create the analysis object
analysis Static
set ok [ analyze 1500]

puts $i
}

puts "node 1 displacement: [nodeDisp 1]"
print node 1
print element

recorder display DisShape 800 10 400 400 -wipe
vup 0 1 0
vpn 0 0 1
prp 2 1.5 10
display 1 5 100

I could not find where the problem is. I would be appreciate it if anybody can help me.

The OpenSees Community

eigen error

eigen error

Re:

Re: eigen error

Re: eigen error

Re: eigen error

Re: eigen error

Re: eigen error