Pushover analysis of nonlinear beam column element and zero

[ericven]

I want to set up a single pier model considering the bond slip of longitudinal reinforcement for Pushover analysis.The code is as follows,But I got the following error.Please check it out for me.
##################################
wipe
file mkdir Data3
model Basic Builder -ndm 2 -ndf 3
source units.tcl
set LCol [expr 1.80*$m]
set Weight [expr 173.*$kN]
set Dia [expr 60*$cm]
set PCol [expr $Weight]
set Mass [expr $PCol/$g]
set ACol [expr 3.14159*pow($Dia,2)/4 ]
set IzCol [expr 3.14159*pow($Dia,4)/64]
node 1 0 0
node 2 0 0
node 3 0 $LCol
puts "nodes"
fix 1 1 1 1
fix 2 1 0 0
mass 1 0 0 0
mass 2 0 0 0
mass 3 $Mass $Mass 0
puts "mass"
set ColSecTag1 1
set ColSecTag2 2
set coverCol [expr 25*$mm]
set Db [expr 14*$mm]
set barArea [expr 3.14159*pow($Db/2,2)]
set IDconcU 1
set IDconc 2
set IDreinf 3
set IDBond 4
set fcU [expr -31.9*$MPa];
set fc [expr -39.4*$MPa];
set Ec [expr 26.7*$GPa];
set fc1U [expr $fcU]
set eps1U -0.003
set fc2U [expr 0.2*$fc1U]
set eps2U -0.006
set lambda 0.1
set ftU [expr -0.07*$fc1U]
set EtsU [expr $ftU/0.002]
set fc1c [expr $fc]
set eps1c -0.002
set fc2c [expr 0.2*$fc1c]
set eps2c -0.014
set ftc [expr -0.07*$fc1c]
set Etsc [expr $ftc/0.002]
set Fy2 [expr 396*$MPa];
set Fu2 [expr 532*$MPa];
set Es2 [expr 200*$GPa];
set Esh2 [expr 3*$GPa]
set esh2 0.007
set eult2 0.15
set Cf 0.20
set alpha 0.506
set Cd 0.389
set alpha1 0.4
set Sy 0.356
#set Sy [expr 254*pow($Db*$Fy2*(2*$alpha1+1)/(8437*pow($fc,0.5)),1/$alpha1)+0.34];
set Su [expr 30*$Sy];
set b 0.5;
set R 0.35;

uniaxialMaterial Concrete02 $IDconcU $fc1U $eps1U $fc2U $eps2U $lambda $ftU $EtsU
uniaxialMaterial Concrete02 $IDconc $fc1c $eps1c $fc2c $eps2c $lambda $ftc $Etsc
uniaxialMaterial ReinforcingSteel $IDreinf $Fy2 $Fu2 $Es2 $Esh2 $esh2 $eult2 -CMFatigue $Cf $alpha $Cd
uniaxialMaterial Bond_SP01 $IDBond $Fy2 $Sy $Fu2 $Su $b $R ;
puts "material"
set exR [expr $Dia/2.0]
set inR [expr $Dia/2.0-$coverCol]
set nfTranc 20
set nfRadc 10
set nfTranu 20
set nfRadu 2
section fiberSec $ColSecTag1 {
patch circ $IDconc $nfTranc $nfRadc 0.0 0.0 0.0 $inR 0.0 360.0
patch circ $IDconcU $nfTranu $nfRadu 0.0 0.0 $inR $exR 0.0 360.0
layer circ $IDreinf 12 $barArea 0.0 0.0 $inR 0.0 360.0
}
section fiberSec $ColSecTag2 {
patch circ $IDconc $nfTranc $nfRadc 0.0 0.0 0.0 $inR 0.0 360.0
patch circ $IDconcU $nfTranu $nfRadu 0.0 0.0 $inR $exR 0.0 360.0
layer circ $IDBond 12 $barArea 0.0 0.0 $inR 0.0 360.0
}
puts "section"
set ColTransfTag 1
geomTransf Linear $ColTransfTag
set numIntgrPts 5
element nonlinearBeamColumn 1 2 3 $numIntgrPts $ColSecTag1 $ColTransfTag
element zeroLengthSection 2 1 2 $ColSecTag2
puts "elements"
recorder Node -file Data3/DFree.out -time -node 3 -dof 1 2 3 disp
recorder Node -file Data3/DBase.out -time -node 1 -dof 1 2 3 disp
recorder Node -file Data3/RBase.out -time -node 1 -dof 1 2 3 reaction
#recorder Drift -file Data3/Drift.out -time -iNode 2 -jNode 3 -dof 1 -perpDirn 2
#recorder Element -file Data3/FCol.out -time -ele 2 globalForce
#recorder Element -file Data3/ForceColSec1.out -time -ele 1 section 1 force
#recorder Element -file Data3/DefoColSec1.out -time -ele 1 section 1 deformation
#recorder Element -file Data3/ForceZero.out -time -ele 2 section 1 stressStrain
puts "recorder"
pattern Plain 1 Linear {
load 3 0 -$PCol 0
}
set Tol 1.0e-3
constraints Plain
numberer Plain
system BandSPD
test NormDispIncr $Tol 6
algorithm Broyden
set NstepGravity 10
set DGravity [expr 1./$NstepGravity]
integrator LoadControl $DGravity
analysis Static
analyze $NstepGravity
loadConst -time 0.0
puts "Model Built"
set Push [expr -0.2*$PCol]
pattern Plain 2 Linear {
load 3 $Push 0 0
}
constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr $Tol 10
algorithm Newton
integrator DisplacementControl 3 1 2.000E-001
analysis Static
analyze 100
puts "Pushover"
#######################################

WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 4
WARNING Broyden::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 4
DisplacementControl::newStep(void) - failed in solver
StaticAnalysis::analyze() - the Integrator failed at iteration: 0 with domain at load factor 0
OpenSees > analyze failed, returned: -2 error flag

[fmk]

you need to provide the axis to the zeolength or else change the orientation of the column, i.e.node 3 0 $Lcol 0

[prczhangqian]

Do you have solved it? I try change the fix mode , it was wrong.And I change the zeolength axis to“ element zeroLengthSection 8 9 1 4 -orient 0 1 0 -1 0 0 ” It also was wrong. I don't know how to change the orientation of the column, i.e.node 3 0 $Lcol 0,as @fmk said.

[selimgunay]

Did you try fixing node 2 in the vertical direction?