is it wrong? help me please

[ocean2000]

i analyze a column just want to make sure the opensees's result. it's very simple. the script as follow.the disp of node 2 is:commitDisps: -2.78508e-016 -0.114927 1.71292e-019
if you try the model in the sap2000V9,concrete and steel material default,delete the per weight, and set the P hige in the column base.then add load for nonlinear anlaysis(consider p-delta too) ,the disp of node 2 in global z is 0.15108(very differ form 0.114927).
--------------------------------------------------
because i try to analyze a rc frame,but find the result is wrong,the opensees get very bigger force than experiment.
so please give me some advice. and what's wrongs in my script follow.
thanks all
--------------------------------------------------
#unit N,mm,c
model basic -ndm 2 -ndf 3
#define nodes
#nodetag x y
node 1 0 0
node 2 0 3000

#define restraint
fix 1 1 1 1

#define material
uniaxialMaterial Concrete01 2 -27.579 -0.00222 -5.5158 -0.0178
uniaxialMaterial Steel01 3 413.6855 199947.98 0.00

#define section
source RCsection.tcl
RCsection 1 200 200 25 2 2 3 3 200 4 4 1 6

#define element
geomTransf PDelta 1
set nP 5
# Columns
# tag ndI ndJ nPts secID transf
element nonlinearBeamColumn 1 1 2 $nP 1 1

#define grcity loads
set p1 -50000
pattern Plain 1 Linear {
load 2 0 $p1 0
}
#------------------
#modle builded
#------------------
#anlayze gracity
system BandGeneral
constraints Transformation
numberer RCM
test NormDispIncr 1.0e-6 10 3
algorithm Newton
integrator LoadControl 0.1 1 0.1 0.1
analysis Static
analyze 10
puts "Gravity load analysis completed"
print node 2
--------------------------------------------------