steel02, $sigInit, not seem to working

[sriram]

Hi silvia and others,
i am trying to put some initial stress in the rebar using $siginit option in steel02. I tried using a truss element to test the material model and i am getting follwing error.. I am posting the my code below.. please let me know if $siginit is not working in steel 02. or if there is anything wrong in my code. I have tried same script with out $siginit, it works.

thanks.
sriram.

My Script:
model basic -ndm 2 -ndf 3
set IDPTBAR 9;
#defining the material parameters for PTbbar
set Fy 245.;
set Es 28800.;
set bratio 0.01;
#definig the material properties
uniaxialMaterial Steel02 $IDPTBAR $Fy $Es $bratio 15 0.925 0. 0. 0. 0. 50.;
node 1 0 0
node 4 1 0
fix 1 1 1 1
fix 4 0 1 1
#geomTransf PDelta 1;
element truss 1 1 4 1 $IDPTBAR
# Define reference load in X direction
pattern Plain 1 "Linear" {
# Fx Fy Mz
load 4 1.0 0 0
}
recorder Node -file testing.out -time -node 4 -dof 1 2 3 disp
set du 0.0001
system SparseGeneral -piv
numberer Plain
constraints Plain
integrator DisplacementControl 4 1 $du
test NormDispIncr 1.0e-5 1000 1
algorithm KrylovNewton
analysis Static
analyze 400


OUTPUT:

Starting OpenSees..


OpenSees -- Open System For Earthquake Engineering Simulation
Pacific Earthquake Engineering Research Center -- 2.0.0

(c) Copyright 1999,2000 The Regents of the University of California
All Rights Reserved
(Copyright and Disclaimer @ http://www.berkeley.edu/OpenSees/copyright.html)



WARNING SuperLU::solve(void)- Error 1 returned in factorization dgstrf()
WARNING KrylovNewton::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 2.88
OpenSees > analyze failed, returned: -3 error flag
---------------------------------------------------------------------------
End of script <C:\tcleditor\bin\material> reached, Press any key to continue