how can i get shear force in the base???

[babibibak]

Hi silvia
i have modeled a one story building which has one bay in each direction.
i want to do pushover analysis.unfortunately i get different reaction in each base. and the sum of the base reactions are much more than base shear according to Etabs model.

wipe
model basic -ndm 3 -ndf 6
#######Units : Kg , cm
set h 300
set b 500
node 1 0 0 0
node 2 $b 0 0
node 3 0 $b 0
node 4 $b $b 0
node 101 0 0 $h
node 102 $b 0 $h
node 103 0 $b $h
node 104 $b $b $h
node 105 [expr $b/2] [expr $b/2] $h
puts "All nodes is created"
####################################
fix 1 1 1 1 1 1 1
fix 2 1 1 1 1 1 1
fix 3 1 1 1 1 1 1
fix 4 1 1 1 1 1 1
fix 105 0 0 1 1 1 0
puts "All supports are fixed"
####################################
rigidDiaphragm 3 105 101 102 103 104
#equalDOF 105 101 1 2 3 4 5 6
#equalDOF 105 102 1 2 3 4 5 6
#equalDOF 105 103 1 2 3 4 5 6
#equalDOF 105 104 1 2 3 4 5 6
puts "Rigid Diaphragm is created"
####################################
###Cover Concrete###################
uniaxialMaterial Concrete01 1 -290 -0.002 0 -0.006
###Core Concrete####################
uniaxialMaterial Concrete01 2 -380 -3.88e-3 -210 -0.202
###Reinforcement Steel##############
uniaxialMaterial Steel01 3 4200 2.1e6 0.02
puts "All material is existed"
####################################
set cover 5
set height 400
set width 400
set y [expr 400/2]
set z [expr 400/2]
set Ac 2.01
section Fiber 3 {
patch rect 2 5 5 [expr $y-$cover] [expr $z-$cover] [expr $cover-$y] [expr $cover-$z]
patch rect 1 5 2 $y $z [expr -$y] [expr $z-$cover]
patch rect 1 5 2 $y [expr $cover-$z] [expr -$y] [expr -$z]
patch rect 1 2 5 $y [expr $z-$cover] [expr $y-$cover] [expr $cover-$z]
patch rect 1 2 5 [expr $cover-$y] [expr $z-$cover] [expr -$y] [expr $cover-$z]
layer straight 3 3 $Ac [expr $y-$cover] [expr $z-$cover] [expr $cover-$y] [expr $z-$cover]
layer straight 3 2 $Ac [expr $y-$cover] 0 [expr $cover-$y] 0
layer straight 3 3 $Ac [expr $y-$cover] [expr $cover-$z] [expr $cover-$y] [expr $cover-$z]
}
puts "Column is defined"
#####################################
set Ab 2.01
section Fiber 4 {
patch rect 2 5 5 [expr $y-$cover] [expr $z-$cover] [expr $cover-$y] [expr $cover-$z]
patch rect 1 5 2 $y $z [expr -$y] [expr $z-$cover]
patch rect 1 5 2 $y [expr $cover-$z] [expr -$y] [expr -$z]
patch rect 1 2 5 $y [expr $z-$cover] [expr $y-$cover] [expr $cover-$z]
patch rect 1 2 5 [expr $cover-$y] [expr $z-$cover] [expr -$y] [expr $cover-$z]
layer straight 3 4 $Ab [expr $y-$cover] [expr $z-$cover] [expr $cover-$y] [expr $z-$cover]
layer straight 3 4 $Ab [expr $y-$cover] [expr $cover-$z] [expr $cover-$y] [expr $cover-$z]
}
puts "Beam is defined"
#####################################
set GA [expr 840000*1333.3]
set GJ [expr 840000*360960]
uniaxialMaterial Elastic 4 $GA
uniaxialMaterial Elastic 5 $GJ
section Aggregator 1 5 T 4 Vy 4 Vz -section 3
section Aggregator 2 5 T 4 Vy 4 Vz -section 4
#####################################
geomTransf Linear 1 0 1 0
set n 5
element nonlinearBeamColumn 1 1 101 $n 1 1
element nonlinearBeamColumn 2 2 102 $n 1 1
element nonlinearBeamColumn 3 3 103 $n 1 1
element nonlinearBeamColumn 4 4 104 $n 1 1
puts "All columns are existed"
######################################
geomTransf Linear 2 0 0 1
element nonlinearBeamColumn 5 101 102 $n 2 2
element nonlinearBeamColumn 6 101 103 $n 2 2
element nonlinearBeamColumn 7 102 104 $n 2 2
element nonlinearBeamColumn 8 103 104 $n 2 2
puts "All beams are existed"
######################################
set gama 2400e-6
set g 9806
set Tslab 20
set Wcol [expr $gama*$width*$height*$h]
set Wbeam [expr $gama*$width*$height*$b]
set Wslab [expr $gama*$Tslab*$b*$b]
set M [expr ($Wcol/2 + $Wbeam/2 + $Wbeam/2 + $Wslab/4)/$g]
#mass 101 $M $M 0 0 0 0
#mass 102 $M $M 0 0 0 0
#mass 103 $M $M 0 0 0 0
#mass 104 $M $M 0 0 0 0
mass 105 [expr 4*$M] [expr 4*$M] 0 0 0 [expr $M*(pow($b,20)+pow($b,2))*(1/3)] ;
set f [expr -3125]
pattern Plain 1 Constant {
load 101 0 0 $f 0 0 0
load 102 0 0 $f 0 0 0
load 103 0 0 $f 0 0 0
load 104 0 0 $f 0 0 0
}
puts "Gravity loads are applied"
######################################
recorder display animation 10 10 500 500 -file 1
prp 4 900 60
vup 0 0 2
fill 2
display 100 5 3
#after 10000
######################################
system BandGeneral
constraints Transformation
numberer RCM
test EnergyIncr 1e-6 20
algorithm Newton
integrator LoadControl 0.1
analysis Static
analyze 10
puts "Gravity analysis is done"
######################################
#set pi 3.141592654 ;
#set N 3 ;
#set Landa [eigen $N] ;
#for {set i 1} {$i <= $N} {incr i 1} { ;
#set T [expr 2*$pi/pow([lindex $Landa [expr $i-1]],0.5)] ;
#puts "********************************************" ;
#puts "T(sec)=" ;
#puts "$T" ;
#puts " " ;
#} ;
######################################
loadConst -time 0.0
pattern Plain 2 Linear {
load 105 42230.0 0.0 0.0 0.0 0.0 0.0
}
recorder Node -file RBASE.txt -node 1 2 3 4 -dof 1 reaction
recorder Node -file DFREE.txt -node 101 102 103 104 -dof 1 disp
set dU 0.01
set maxU 1.0
set numSteps [expr int($maxU/$dU)]
constraints Transformation
numberer RCM
system BandGeneral
test EnergyIncr 0.1 1000 0
algorithm Newton
integrator DisplacementControl 105 1 $dU
analyze $numSteps
puts "Pushover analysis is done"


where is the problem?

[silvia]

i don't know, the models may be different.