modeling of space structurs(dome structure)

[mohammadkh]

Hello dear
I have modeled an unit dome in OpenSees. when I want to analyze eigen values,the following error shows. Is there anyone who can help me to solve this problem?
** On entry to DLASCL, parameter number 4 had an illegal value
** On entry to DLASCL, parameter number 4 had an illegal value
ArpackSolver::Error with _saupd info = -9999
Could not build an Arnoldi factorization.IPARAM(5) the size of the current Arnol
di factorization: is 1factorization. The user is advised to check thatenough wor
kspace and array storage has been allocated.
WARNING DirectIntegrationAnalysis::eigen() - EigenSOE failed in solve()
** On entry to DLASCL, parameter number 4 had an illegal value
** On entry to DLASCL, parameter number 4 had an illegal value
ArpackSolver::Error with _saupd info = -9999
Could not build an Arnoldi factorization.IPARAM(5) the size of the current Arnol
di factorization: is 1factorization. The user is advised to check thatenough wor
kspace and array storage has been allocated.
WARNING DirectIntegrationAnalysis::eigen() - EigenSOE failed in solve()
missing function argument at _@_
in expression "pow ( _@_, 0.5)"
(parsing expression "pow ( , 0.5)")
invoked from within
"expr pow ($w12 , 0.5)"
invoked from within
"set w1 [expr pow ($w12 , 0.5)]"


Also, I have attach my model.
wipe all ;
model basic -ndm 3 -ndf 6
file mkdir "verify"

node 1 0. 0. 0.157
node 2 3. 0. 0.
node 3 1.5 2.6 0.
node 4 -1.5 2.6 0.
node 5 -3. 0. 0.
node 6 -1.5 -2.6 0.
node 7 1.5 -2.6 0.

puts "end of node"

fix 2 1 1 1 0 0 0
fix 3 1 1 1 0 0 0
fix 4 1 1 1 0 0 0
fix 5 1 1 1 0 0 0
fix 6 1 1 1 0 0 0

puts "end of fix"

uniaxialMaterial Steel02 1 2400.E4 2.1E10 0.02 10 0.925 0.15

puts "end of material"

section Fiber 1 {;
patch circ 1 8 1 0. 0. 6.096E-3 6.35E-3 0. 0.
};

puts "end of Fiber"

geomTransf Corotational 1 1 0 1
geomTransf Corotational 2 0 1 1
geomTransf Corotational 3 0 0 1
geomTransf Corotational 4 -1 0 1
geomTransf Corotational 5 0 -1 1

puts "end of geom"


element dispBeamColumn 1 1 2 10 1 1
element dispBeamColumn 2 1 3 10 1 2
element dispBeamColumn 3 2 3 10 1 3
element dispBeamColumn 4 1 4 10 1 3
element dispBeamColumn 5 3 4 10 1 3
element dispBeamColumn 6 1 5 10 1 4
element dispBeamColumn 7 4 5 10 1 3
element dispBeamColumn 8 1 6 10 1 5
element dispBeamColumn 9 5 6 10 1 3
element dispBeamColumn 10 1 7 10 1 3
element dispBeamColumn 11 6 7 10 1 3
element dispBeamColumn 12 7 2 10 1 3


puts "end of element"


mass 1 0.00000001 0.000000001 102. 0.000000001 0.00000001 0.000000001

puts "end of mass"


eigen 1


recorder display space 100 20 400 300 -wipe
vup 0 0 1
vpn 0 1 0
vrp 0 4 0
prp -10. -5. -30.
viewWindow -6 5 -6 10.
display 1 2 1
after 1000

set pi 3.1415
set a [eigen 1]
set w12 [lindex $a 0]
#set w22 [lindex $a 1]
#set w32 [lindex $a 2]
set w1 [expr pow ($w12 , 0.5)]
#set w2 [expr pow ($w22 , 0.5)]
#set w3 [expr pow ($w32 , 0.5)]
set T1 [expr (2*$pi)/$w1]
#set T2 [expr (2*$pi)/$w2]
#set T3 [expr (2*$pi)/$w3]

puts "T1 = $T1"
#puts "T2 = $T2"
#puts "T3 = $T3"


# Static Analysis

wipeAnalysis ;
constraints Transformation ;
numberer RCM ;
system SparseGeneral ;
test EnergyIncr 1e-7 25 1 ;
algorithm ModifiedNewton ;
integrator LoadControl 0.01 ;
analysis Static ;
analyze 100 ;
loadConst -time 0.0 ;

[naserhassania]

Hi ,

I am just wondering if you did solve this problem that was causing this error.
If yes what was it because I am getting something similar to it

Regards
Naser

[fmk]

the original post was due to the 2 angles being the same .. resulting in a gibberish section definition.