algorithm failed error

[rasoulazadi]

hello dear
i write this program and have problem with solution and analyse and has very large periods
please help me to solve this problem
thanks alot

# OpenSees beyrami and sadat
# OpenSees Primer
#
# Units: kg, m, sec

# ------------------------------
# Start of model generation
# ------------------------------
# Create ModelBuilder (with two-dimensions and 3 DOF/node)
model BasicBuilder -ndm 2 -ndf 3
# Create nodes
# ------------
source nodes.txt
# Define materials for nonlinear elements
#uniaxialMaterial Steel02 1 35153481 20390000000 0.03 20 0.925 0.15
#uniaxialMaterial Steel02 2 25310930 20390000000 0.03 20 0.925 0.15
uniaxialMaterial Hysteretic 1 35153481 0.001724 35153481 0.015 44300660 0.17 -35153481 -0.001724 -35153481 -0.015 -44300660 -0.17 1 1 0.0 0.0 0.0




geomTransf Linear 1
geomTransf Linear 2
source Wsection.tcl

#define sections for columns

#tag matID d tw bf tf nfdw nftw nfbf nftf
Wsection 1 1 0.82550 0.050 0.3886 0.0899 10 1 10 1
Wsection 2 1 0.7120 0.003860 0.348 0.0691 10 1 10 1
Wsection 3 1 0.28960 0.02260 0.2896 0.0404 10 1 10 1
Wsection 4 1 0.28960 0.02260 0.2896 0.0404 10 1 10 1
Wsection 5 1 0.4750 0.0150 0.2845 0.0239 10 1 10 1
#section for beams
Wsection 6 1 0.30990 0.009398 0.2052 0.0163 10 1 10 1




source elements2.txt
puts "end definitions"

####################################################################################
# Gravity Analysis #
####################################################################################

set GravSteps 10

constraints Plain;
numberer RCM; # renumber dof's to minimize band-width
puts "end definitions1"


system UmfPack;
test NormDispIncr 1.0e-6 53
puts "end definitions"


algorithm KrylovNewton
puts "end definitions"


integrator LoadControl 0.1
puts "end definitions4"


analysis Static
analyze $GravSteps
puts "end definitions"


puts "Gravity Analyze Perform Successfully..."
loadConst -time 0.0;

####################################################################################
# Eigen Analysis #
####################################################################################
set pi 3.14
set autoval [eigen 3]
puts "
EIGEN VALUES: $autoval"
set T1 [expr 2*$pi/pow([lindex $autoval 0],.5)]
set T2 [expr 2*$pi/pow([lindex $autoval 1],.5)]
set T3 [expr 2*$pi/pow([lindex $autoval 2],.5)]
puts "PERIODS: $T1 , $T2 , $T3
* Finished Eigen Analysis After Gravity Anlalyze... *
"

[fmk]

if the periods are wrong then obviously your model is wrong .. you probably have your units messed up

[rasoulazadi]

hi dear vesna
i am modeling shear wall using shell element
for section define i am using fiber section
i should consider the 0,0 of section in the center or in the i point of shell element

[vesna]

0,0 is usually at the centroid of the section.

[rasoulazadi]

thanks for replay dear vesna
i have an question
how can i consider shear properties of wall if i use fiber sections?

[vesna]

You can use section aggregator command (http://opensees.berkeley.edu/wiki/index ... Aggregator) to add shear force-deformation relationship to a fiber section. However, be aware that although you have P-M interaction you won't have P-M-V interaction.