How can I fix my errors of ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii)?

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

Moderators: silvia, selimgunay, Moderators

Post Reply
kangmur
Posts: 3
Joined: Mon Feb 07, 2022 9:49 pm

How can I fix my errors of ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii)?

Post by kangmur » Tue Feb 08, 2022 11:16 pm

The errors are:
ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (29, 0)
ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (110, 0)
ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (218, 0)
ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (292, 0)

T1 = 1.4199903451566205 s
T2 = 1.2115937173449813 s
Eigen Analysis Done
WARNING: CTestNormDispIncr::setEquiSolnAlgo() - no SOE
Analysis_1 k 1 j 1 lambda 1 Dt 0.001 Tol 1.0e-4 Algorithm NewtonLineSearch -type Bisection
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 30
WARNING NewtonLineSearch::solveCurrentStep() -the LinearSysOfEqn failed in solve()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.001
OpenSees > analyze failed, returned: -3 error flag

I think that it's because of orient of zerolengthsection. Is it right?


kangmur
Posts: 3
Joined: Mon Feb 07, 2022 9:49 pm

Re: How can I fix my errors of ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii)?

Post by kangmur » Wed Feb 09, 2022 6:02 pm

mhscott wrote:
Wed Feb 09, 2022 8:10 am
This post explains the issue: https://portwooddigital.com/2021/09/12/ ... -to-solve/
Thank you for reply.

I read this post and checked my model.
But I cannot find any disconnection, floating, problem of DOF or else.
Can you check my script and give some advice?

=======================================================================
wipe all

model BasicBuilder -ndm 3 -ndf 6

geomTransf PDelta 1 0 0 1
geomTransf PDelta 2 1 0 0

set IDctrlNode 41
set IDctrlDOF 1

file mkdir Results3D

source DisplayModel3D.tcl
source DisplayPlane.tcl
source LibUnits.tcl

#model parameter
set h1 0.200000
set b1 0.200000
set h2 0.250000
set b2 0.150000
set h3 0.250000
set b3 0.150000
set Ec 23000000.000000
set nu 0.200000
set g 9.810000

#node
node 11 0.000000 0.000000 0.000000
node 12 1.500000 0.000000 0.000000
node 13 0.000000 0.000000 0.750000
node 14 1.500000 0.000000 0.750000

node 21 0.000000 1.200000 0.000000
node 22 1.500000 1.200000 0.000000
node 23 0.000000 1.200000 0.750000
node 24 1.500000 1.200000 0.750000

node 31 0.000000 2.400000 0.000000
node 32 1.500000 2.400000 0.000000
node 33 0.000000 2.400000 0.750000
node 34 1.500000 2.400000 0.750000

node 41 0.000000 3.600000 0.000000
node 42 1.500000 3.600000 0.000000
node 43 0.000000 3.600000 0.750000
node 44 1.500000 3.600000 0.750000

#boundary condition
fix 11 1 1 1 1 1 1
fix 12 1 1 1 1 1 1
fix 13 1 1 1 1 1 1
fix 14 1 1 1 1 1 1

#rigiddiaphragm
set RigidDiaphragm ON
set perpDirn 2
node 2 0.750000 1.200000 0.375000
fix 2 0 1 0 1 0 1
node 3 0.750000 2.400000 0.375000
fix 3 0 1 0 1 0 1
node 4 0.750000 3.600000 0.375000
fix 4 0 1 0 1 0 1

rigidDiaphragm $perpDirn 2 21 22 23 24
rigidDiaphragm $perpDirn 3 31 32 33 34
rigidDiaphragm $perpDirn 4 41 42 43 44

#zerolength node
node 111 0.000000 0.000000 0.000000
node 121 1.500000 0.000000 0.000000
node 131 0.000000 0.000000 0.750000
node 141 1.500000 0.000000 0.750000
node 211 0.000000 1.200000 0.000000
node 221 1.500000 1.200000 0.000000
node 231 0.000000 1.200000 0.750000
node 241 1.500000 1.200000 0.750000
node 311 0.000000 2.400000 0.000000
node 321 1.500000 2.400000 0.000000
node 331 0.000000 2.400000 0.750000
node 341 1.500000 2.400000 0.750000
node 213 0.000000 1.200000 0.000000
node 223 1.500000 1.200000 0.000000
node 233 0.000000 1.200000 0.750000
node 243 1.500000 1.200000 0.750000
node 313 0.000000 2.400000 0.000000
node 323 1.500000 2.400000 0.000000
node 333 0.000000 2.400000 0.750000
node 343 1.500000 2.400000 0.750000
node 413 0.000000 3.600000 0.000000
node 423 1.500000 3.600000 0.000000
node 433 0.000000 3.600000 0.750000
node 443 1.500000 3.600000 0.750000

node 212 0.000000 1.200000 0.000000
node 224 1.500000 1.200000 0.000000
node 312 0.000000 2.400000 0.000000
node 324 1.500000 2.400000 0.000000
node 412 0.000000 3.600000 0.000000
node 424 1.500000 3.600000 0.000000
node 232 0.000000 1.200000 0.750000
node 244 1.500000 1.200000 0.750000
node 332 0.000000 2.400000 0.750000
node 344 1.500000 2.400000 0.750000
node 432 0.000000 3.600000 0.750000
node 444 1.500000 3.600000 0.750000
node 215 0.000000 1.200000 0.000000
node 225 1.500000 1.200000 0.000000
node 315 0.000000 2.400000 0.000000
node 325 1.500000 2.400000 0.000000
node 415 0.000000 3.600000 0.000000
node 425 1.500000 3.600000 0.000000
node 236 0.000000 1.200000 0.750000
node 246 1.500000 1.200000 0.750000
node 336 0.000000 2.400000 0.750000
node 346 1.500000 2.400000 0.750000
node 436 0.000000 3.600000 0.750000
node 446 1.500000 3.600000 0.750000

#bilin for column
uniaxialMaterial Bilin 9111 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1
uniaxialMaterial Bilin 9112 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1
uniaxialMaterial Bilin 9113 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1
uniaxialMaterial Bilin 9114 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1
uniaxialMaterial Bilin 9121 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1
uniaxialMaterial Bilin 9122 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1
uniaxialMaterial Bilin 9123 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1
uniaxialMaterial Bilin 9124 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1
uniaxialMaterial Bilin 9131 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1
uniaxialMaterial Bilin 9132 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1
uniaxialMaterial Bilin 9133 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1
uniaxialMaterial Bilin 9134 1533.333333 0.050000 0.050000 7288.863000 -7288.863000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.042000 0.042000 0.079292 0.079292 0.200000 0.200000 0.147645 0.147645 1 1

#bilin for beam
uniaxialMaterial Bilin 9211 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1
uniaxialMaterial Bilin 9212 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1
uniaxialMaterial Bilin 9213 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1
uniaxialMaterial Bilin 9214 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1
uniaxialMaterial Bilin 9221 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1
uniaxialMaterial Bilin 9222 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1
uniaxialMaterial Bilin 9223 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1
uniaxialMaterial Bilin 9224 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1
uniaxialMaterial Bilin 9231 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1
uniaxialMaterial Bilin 9232 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1
uniaxialMaterial Bilin 9233 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1
uniaxialMaterial Bilin 9234 2246.093750 0.050000 0.050000 4062.290000 -4062.290000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074495 0.074495 0.121727 0.121727 0.200000 0.200000 0.202511 0.202511 1 1

#section aggregator
uniaxialMaterial Elastic 9911 5266670.000000
uniaxialMaterial Elastic 9811 2289860.000000
uniaxialMaterial Elastic 9711 1000000000.000000
section Aggregator 1111 9911 P 9111 Mz 9811 Vy 9711 T

uniaxialMaterial Elastic 9912 5266670.000000
uniaxialMaterial Elastic 9812 2289860.000000
uniaxialMaterial Elastic 9712 1000000000.000000
section Aggregator 1112 9912 P 9112 Mz 9812 Vy 9712 T

uniaxialMaterial Elastic 9913 5266670.000000
uniaxialMaterial Elastic 9813 2289860.000000
uniaxialMaterial Elastic 9713 1000000000.000000
section Aggregator 1113 9913 P 9113 Mz 9813 Vy 9713 T

uniaxialMaterial Elastic 9914 5266670.000000
uniaxialMaterial Elastic 9814 2289860.000000
uniaxialMaterial Elastic 9714 1000000000.000000
section Aggregator 1114 9914 P 9114 Mz 9814 Vy 9714 T

uniaxialMaterial Elastic 9921 5266670.000000
uniaxialMaterial Elastic 9821 2289860.000000
uniaxialMaterial Elastic 9721 1000000000.000000
section Aggregator 1121 9921 P 9121 Mz 9821 Vy 9721 T

uniaxialMaterial Elastic 9922 5266670.000000
uniaxialMaterial Elastic 9822 2289860.000000
uniaxialMaterial Elastic 9722 1000000000.000000
section Aggregator 1122 9922 P 9122 Mz 9822 Vy 9722 T

uniaxialMaterial Elastic 9923 5266670.000000
uniaxialMaterial Elastic 9823 2289860.000000
uniaxialMaterial Elastic 9723 1000000000.000000
section Aggregator 1123 9923 P 9123 Mz 9823 Vy 9723 T

uniaxialMaterial Elastic 9924 5266670.000000
uniaxialMaterial Elastic 9824 2289860.000000
uniaxialMaterial Elastic 9724 1000000000.000000
section Aggregator 1124 9924 P 9124 Mz 9824 Vy 9724 T

uniaxialMaterial Elastic 9931 5266670.000000
uniaxialMaterial Elastic 9831 2289860.000000
uniaxialMaterial Elastic 9731 1000000000.000000
section Aggregator 1131 9931 P 9131 Mz 9831 Vy 9731 T

uniaxialMaterial Elastic 9932 5266670.000000
uniaxialMaterial Elastic 9832 2289860.000000
uniaxialMaterial Elastic 9732 1000000000.000000
section Aggregator 1132 9932 P 9132 Mz 9832 Vy 9732 T

uniaxialMaterial Elastic 9933 5266670.000000
uniaxialMaterial Elastic 9833 2289860.000000
uniaxialMaterial Elastic 9733 1000000000.000000
section Aggregator 1133 9933 P 9133 Mz 9833 Vy 9733 T

uniaxialMaterial Elastic 9934 5266670.000000
uniaxialMaterial Elastic 9834 2289860.000000
uniaxialMaterial Elastic 9734 1000000000.000000
section Aggregator 1134 9934 P 9134 Mz 9834 Vy 9734 T

uniaxialMaterial Elastic 9611 6846680.000000
uniaxialMaterial Elastic 9511 2976820.000000
uniaxialMaterial Elastic 9411 1000000000.000000
section Aggregator 1221 9611 P 9211 Mz 9511 Vy 9411 T

uniaxialMaterial Elastic 9612 6846680.000000
uniaxialMaterial Elastic 9512 2976820.000000
uniaxialMaterial Elastic 9412 1000000000.000000
section Aggregator 1222 9612 P 9212 Mz 9512 Vy 9412 T

uniaxialMaterial Elastic 9613 6846680.000000
uniaxialMaterial Elastic 9513 2976820.000000
uniaxialMaterial Elastic 9413 1000000000.000000
section Aggregator 1223 9613 P 9213 Mz 9513 Vy 9413 T

uniaxialMaterial Elastic 9614 6846680.000000
uniaxialMaterial Elastic 9514 2976820.000000
uniaxialMaterial Elastic 9414 1000000000.000000
section Aggregator 1224 9614 P 9214 Mz 9514 Vy 9414 T

uniaxialMaterial Elastic 9621 6846680.000000
uniaxialMaterial Elastic 9521 2976820.000000
uniaxialMaterial Elastic 9421 1000000000.000000
section Aggregator 1231 9621 P 9221 Mz 9521 Vy 9421 T

uniaxialMaterial Elastic 9622 6846680.000000
uniaxialMaterial Elastic 9522 2976820.000000
uniaxialMaterial Elastic 9422 1000000000.000000
section Aggregator 1232 9622 P 9222 Mz 9522 Vy 9422 T

uniaxialMaterial Elastic 9623 6846680.000000
uniaxialMaterial Elastic 9523 2976820.000000
uniaxialMaterial Elastic 9423 1000000000.000000
section Aggregator 1233 9623 P 9223 Mz 9523 Vy 9423 T

uniaxialMaterial Elastic 9624 6846680.000000
uniaxialMaterial Elastic 9524 2976820.000000
uniaxialMaterial Elastic 9424 1000000000.000000
section Aggregator 1234 9624 P 9224 Mz 9524 Vy 9424 T

uniaxialMaterial Elastic 9631 6846680.000000
uniaxialMaterial Elastic 9531 2976820.000000
uniaxialMaterial Elastic 9431 1000000000.000000
section Aggregator 1241 9631 P 9231 Mz 9531 Vy 9431 T

uniaxialMaterial Elastic 9632 6846680.000000
uniaxialMaterial Elastic 9532 2976820.000000
uniaxialMaterial Elastic 9432 1000000000.000000
section Aggregator 1242 9632 P 9232 Mz 9532 Vy 9432 T

uniaxialMaterial Elastic 9633 6846680.000000
uniaxialMaterial Elastic 9533 2976820.000000
uniaxialMaterial Elastic 9433 1000000000.000000
section Aggregator 1243 9633 P 9233 Mz 9533 Vy 9433 T

uniaxialMaterial Elastic 9634 6846680.000000
uniaxialMaterial Elastic 9534 2976820.000000
uniaxialMaterial Elastic 9434 1000000000.000000
section Aggregator 1244 9634 P 9234 Mz 9534 Vy 9434 T

#column element
element elasticBeamColumn 3011 111 213 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1
element elasticBeamColumn 3012 121 223 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1
element elasticBeamColumn 3013 131 233 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1
element elasticBeamColumn 3014 141 243 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1

element elasticBeamColumn 3021 211 313 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1
element elasticBeamColumn 3022 221 323 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1
element elasticBeamColumn 3023 231 333 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1
element elasticBeamColumn 3024 241 343 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1

element elasticBeamColumn 3031 311 413 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1
element elasticBeamColumn 3032 321 423 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1
element elasticBeamColumn 3033 331 433 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1
element elasticBeamColumn 3034 341 443 0.040000 23000000.000000 9583333.333333 0.000267 0.000133 0.000133 1

#beam element
element elasticBeamColumn 3121 212 224 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 1
element elasticBeamColumn 3221 215 236 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 2
element elasticBeamColumn 3122 232 244 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 1
element elasticBeamColumn 3222 225 246 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 2

element elasticBeamColumn 3131 312 324 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 1
element elasticBeamColumn 3231 315 336 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 2
element elasticBeamColumn 3132 332 344 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 1
element elasticBeamColumn 3232 325 346 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 2

element elasticBeamColumn 3141 412 424 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 1
element elasticBeamColumn 3241 415 436 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 2
element elasticBeamColumn 3142 432 444 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 1
element elasticBeamColumn 3242 425 446 0.037500 23000000.000000 9583333.333333 0.000266 0.000070 0.000195 2

#element for node and zerolength node(column)
element zeroLengthSection 4111 11 111 1111 -orient 0 1 0 1 0 0
element zeroLengthSection 5111 213 21 1111 -orient 0 1 0 0 0 1
element zeroLengthSection 4112 12 121 1112 -orient 0 1 0 1 0 0
element zeroLengthSection 5112 223 22 1112 -orient 0 1 0 0 0 1
element zeroLengthSection 4113 13 131 1113 -orient 0 1 0 1 0 0
element zeroLengthSection 5113 233 23 1113 -orient 0 1 0 0 0 1
element zeroLengthSection 4114 14 141 1114 -orient 0 1 0 1 0 0
element zeroLengthSection 5114 243 24 1114 -orient 0 1 0 0 0 1

element zeroLengthSection 4121 21 211 1121 -orient 0 1 0 1 0 0
element zeroLengthSection 5121 313 31 1121 -orient 0 1 0 0 0 1
element zeroLengthSection 4122 22 221 1122 -orient 0 1 0 1 0 0
element zeroLengthSection 5122 323 32 1122 -orient 0 1 0 0 0 1
element zeroLengthSection 4123 23 231 1123 -orient 0 1 0 1 0 0
element zeroLengthSection 5123 333 33 1123 -orient 0 1 0 0 0 1
element zeroLengthSection 4124 24 241 1124 -orient 0 1 0 1 0 0
element zeroLengthSection 5124 343 34 1124 -orient 0 1 0 0 0 1

element zeroLengthSection 4131 31 311 1131 -orient 0 1 0 1 0 0
element zeroLengthSection 5131 413 41 1131 -orient 0 1 0 0 0 1
element zeroLengthSection 4132 32 321 1132 -orient 0 1 0 1 0 0
element zeroLengthSection 5132 423 42 1132 -orient 0 1 0 0 0 1
element zeroLengthSection 4133 33 331 1133 -orient 0 1 0 1 0 0
element zeroLengthSection 5133 433 43 1133 -orient 0 1 0 0 0 1
element zeroLengthSection 4134 34 341 1134 -orient 0 1 0 1 0 0
element zeroLengthSection 5134 443 44 1134 -orient 0 1 0 0 0 1

#element for node and zerolength node(beam)
element zeroLengthSection 6121 21 212 1221 -orient 1 0 0 0 1 0
element zeroLengthSection 7121 21 215 1221 -orient 0 0 1 1 0 0
element zeroLengthSection 6122 224 22 1222 -orient 1 0 0 0 1 0
element zeroLengthSection 7122 22 225 1222 -orient 0 0 1 1 0 0
element zeroLengthSection 6123 23 232 1223 -orient 1 0 0 0 0 1
element zeroLengthSection 7123 236 23 1223 -orient 0 0 1 1 0 0
element zeroLengthSection 6124 244 24 1224 -orient 1 0 0 0 1 0
element zeroLengthSection 7124 246 24 1224 -orient 0 0 1 1 0 0

element zeroLengthSection 6131 31 312 1231 -orient 1 0 0 0 1 0
element zeroLengthSection 7131 31 315 1231 -orient 0 0 1 1 0 0
element zeroLengthSection 6132 324 32 1232 -orient 1 0 0 0 1 0
element zeroLengthSection 7132 32 325 1232 -orient 0 0 1 1 0 0
element zeroLengthSection 6133 33 332 1233 -orient 1 0 0 0 0 1
element zeroLengthSection 7133 336 33 1233 -orient 0 0 1 1 0 0
element zeroLengthSection 6134 344 34 1234 -orient 1 0 0 0 1 0
element zeroLengthSection 7134 346 34 1234 -orient 0 0 1 1 0 0

element zeroLengthSection 6141 41 412 1241 -orient 1 0 0 0 1 0
element zeroLengthSection 7141 41 415 1241 -orient 0 0 1 1 0 0
element zeroLengthSection 6142 424 42 1242 -orient 1 0 0 0 1 0
element zeroLengthSection 7142 42 425 1242 -orient 0 0 1 1 0 0
element zeroLengthSection 6143 43 432 1243 -orient 1 0 0 0 0 1
element zeroLengthSection 7143 436 43 1243 -orient 0 0 1 1 0 0
element zeroLengthSection 6144 444 44 1244 -orient 1 0 0 0 1 0
element zeroLengthSection 7144 446 44 1244 -orient 0 0 1 1 0 0

#mass
mass 21 0.012769 0.000000 0.012769 0.000000 0.000000 0.000000
mass 31 0.012769 0.000000 0.012769 0.000000 0.000000 0.000000
mass 41 0.009371 0.000000 0.009371 0.000000 0.000000 0.000000
mass 22 0.012769 0.000000 0.012769 0.000000 0.000000 0.000000
mass 32 0.012769 0.000000 0.012769 0.000000 0.000000 0.000000
mass 42 0.009371 0.000000 0.009371 0.000000 0.000000 0.000000
mass 23 0.012769 0.000000 0.012769 0.000000 0.000000 0.000000
mass 33 0.012769 0.000000 0.012769 0.000000 0.000000 0.000000
mass 43 0.009371 0.000000 0.009371 0.000000 0.000000 0.000000
mass 24 0.012769 0.000000 0.012769 0.000000 0.000000 0.000000
mass 34 0.012769 0.000000 0.012769 0.000000 0.000000 0.000000
mass 44 0.009371 0.000000 0.009371 0.000000 0.000000 0.000000

puts "Model Built"

recorder Node -file Results3D/Disp1.out -time -node 41 -dof 1 disp

#period
set pi [expr 2.0*asin(1.0)]
set nEigen 3
set lambdaN [eigen [expr $nEigen]]
set lambdaI [lindex $lambdaN 0]
set lambdaJ [lindex $lambdaN 1]
set w1 [expr pow($lambdaI,0.5)]
set w2 [expr pow($lambdaJ,0.5)]
set T1 [expr 2.0*$pi/$w1]
set T2 [expr 2.0*$pi/$w2]
puts "T1 = $T1 s"
puts "T2 = $T2 s"

set period "Results3D/Periods.txt"
set Periods [open $period "w"]
foreach t $T1 {puts $Periods "$t"}
close $Periods

puts "Eigen Analysis Done"

pattern Plain 100 Linear {
eleLoad -ele 3121 -type -beamUniform 0 -11.74
eleLoad -ele 3122 -type -beamUniform 0 -11.74
eleLoad -ele 3131 -type -beamUniform 0 -11.74
eleLoad -ele 3132 -type -beamUniform 0 -11.74
eleLoad -ele 3141 -type -beamUniform 0 -11.74
eleLoad -ele 3142 -type -beamUniform 0 -11.74
eleLoad -ele 3221 -type -beamUniform 0 -11.74
eleLoad -ele 3222 -type -beamUniform 0 -11.74
eleLoad -ele 3231 -type -beamUniform 0 -11.74
eleLoad -ele 3232 -type -beamUniform 0 -11.74
eleLoad -ele 3241 -type -beamUniform 0 -11.74
eleLoad -ele 3242 -type -beamUniform 0 -11.74

}

#time history analysis
set k 1
#scaling ID
set j 1

#rayleigh damping coefficients
set xi1 0.05
set xi2 0.05
set pi [expr acos(-1)]
set T1 $T1
set T2 $T2
set omega1 [expr 2*$pi/$T1]
set omega2 [expr 2*$pi/$T2]
#set aR [expr 2*($omega1*$omega2*($omega2*$xi1-$omega1*$xi2))/($omega2**2-$omega1**2)]
#set bR [expr 2*($omega2*$xi2-$omega1*$xi1)/($omega2**2-$omega1**2)]
set aR [expr 0.02*(2*$omega1*$omega2)/($omega1+$omega2)]
set bR 0

#ground motion input
set lambda 1
set dt 0.01
set numstep 4000
set GMfile "BM68elc.acc"
set g 9.81
#TimeSeries "TimeSeries01": dt filePath cFactor
set accelSeries "Series -dt $dt -filePath $GMfile -factor [expr $lambda*$g]"

#analysis parameters
#definition of the tolerance and of the maximum number of iteration
set Tol 1.0e-4
set maxIter 5000
#convergence test
#test NormDispIncr $tol $iter
test NormDispIncr $Tol $maxIter 0 0
#define load pattern
#pattern UniformExcitation $patternTag $dir -accel $tsTag <vel0 $vel0> <-fact $cFactor>
pattern UniformExcitation 3 1 -accel $accelSeries
#Constraint Handlerconstraints Transformation
constraints Transformation
#Integrator
#integrator Newmark $gamma $beta
integrator Newmark 0.5 0.25
#Rayleigh Damping
#rayleigh $alphaM $betaK $betaKinit $betaKcomm
rayleigh $aR $bR +0.000000E+000 +0.000000E+000
#Solution Algorithm (default minEta=0.1 & maxEta=10, default tol=0.8 & maxIter=10)
#algorithm NewtonLineSearch <-type $typeSearch> <-tol $tol> -maxIter $maxIter <-minEta $minEta> <-maxEta $maxEta>
algorithm NewtonLineSearch
numberer RCM
system BandGeneral
analysis Transient
#print node

set begin [clock clicks -milliseconds]

#analysis
#Initial variables before starting dynamic analysis
set ok 0
set t 0
set startTime [clock clicks -milliseconds]
set finalt [expr [getTime]+[expr $numstep*$dt]]
set DtAnalysis [expr 0.1*$dt]
puts "Analysis_1 k $k j $j lambda $lambda Dt $DtAnalysis Tol $Tol Algorithm NewtonLineSearch -type Bisection"
while {$ok == 0 && $t <= $finalt} {
set ok [analyze 1 $DtAnalysis]
set t [getTime]
set currentTime [getTime]
}
set finishTime [clock clicks -milliseconds]
set timeSeconds [expr ($finishTime-$startTime)/1000]
set timeMinutes [expr ($timeSeconds/60)]
set timeHours [expr ($timeSeconds/3600)]
set timeMinutes [expr ($timeMinutes - $timeHours*60)]
set timeSeconds [expr ($timeSeconds - $timeMinutes*60 - $timeHours*3600)]

if {$ok == 0} {
puts "################################################"
puts " Analysis k $k j $j lambda $lambda completed SUCCESSFULLY";
puts "Time in hours-min-sec: $timeHours:$timeMinutes:$timeSeconds"
puts "################################################"
} else {
puts "################################################"
puts "Analysis k $k j $j lambda $lambda FAILED";
puts "Time in hours-min-sec: $timeHours:$timeMinutes:$timeSeconds"
puts "################################################"
}

wipe all

mhscott
Posts: 868
Joined: Tue Jul 06, 2004 3:38 pm
Location: Corvallis, Oregon USA
Contact:

Re: How can I fix my errors of ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii)?

Post by mhscott » Thu Feb 10, 2022 5:47 am

No, sorry, I don't have any more advice.

kangmur
Posts: 3
Joined: Mon Feb 07, 2022 9:49 pm

Re: How can I fix my errors of ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii)?

Post by kangmur » Thu Feb 10, 2022 8:12 am

mhscott wrote:
Thu Feb 10, 2022 5:47 am
No, sorry, I don't have any more advice.
Thanks a lot!

mhelsaye
Posts: 1
Joined: Sun May 23, 2021 8:11 am

Re: How can I fix my errors of ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii)?

Post by mhelsaye » Fri Mar 04, 2022 2:36 pm

Hello,

Have you found your error?
I am facing the same issue. The model run when I apply some constraints to the zerosectionelement node such as equalDoF in the X and Y direction but I am not sure about the results

jiaojiao
Posts: 1
Joined: Thu Apr 13, 2023 1:40 am

Re: How can I fix my errors of ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii)?

Post by jiaojiao » Sun Apr 23, 2023 1:16 am

Hello,

Have you found your error?
I am facing the same issue.

destanylynch
Posts: 1
Joined: Wed Apr 26, 2023 12:32 am

Re: How can I fix my errors of ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii)?

Post by destanylynch » Wed Apr 26, 2023 12:34 am

When I put constraints to the zerosectionelement node, such as equalDoF in the X and Y directions, the model runs, but I'm not sure how the results will turn out. How should I do next.

getting over it

SamanTarokh
Posts: 1
Joined: Wed Jun 07, 2023 1:33 pm

Re: How can I fix my errors of ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii)?

Post by SamanTarokh » Wed Jun 07, 2023 1:37 pm

check the Direct of your material with direct of your zerolength elements

Post Reply