PressureDependMultiYield Material

[zhoulinlu]

Hi,
Please point out the problem with the following command for me.
wipe;
model Basic -ndm 2 -ndf 3
# Define nodes
node 1 0.0 0.0
node 2 0.0 0.0
# Fix node 1
fix 1 1 1 1
fix 2 0 1 0

set accMul 2 ;# acceleration multiplier
set period 1 ;# Period of applied sinusoidal load

set massProportionalDamping 0.0 ;
set stiffnessProportionalDamping 0.02 ;
set gamma 0.600 ;# Newmark integration parameter
set E 90000.0 ;#shear modulus
set poisson 0.40 ;
set G [expr $E/(2*(1+$poisson))] ;
set B [expr $E/(3*(1-2*$poisson))] ;

nDMaterial PressureDependMultiYield 2 2 2 $G $B 31.4 0.1 80 0.5 26.5 0.17 0.4 10 10 0.015 1.0
nDMaterial FluidSolidPorous 1 2 2 2.2e+6

element zeroLengthND 1 1 2 2

updateMaterialStage -material 2 -stage 0

system ProfileSPD
test NormDispIncr 1.e-12 25 0
constraints Transformation
integrator LoadControl 1 1 1 1
algorithm Newton
numberer RCM

# create the Analysis
analysis Static

#analyze
analyze 2

# switch the material to plastic
updateMaterialStage -material 1 -stage 1
updateMaterialStage -material 2 -stage 1
updateMaterials -material 2 bulkModulus [expr $G*2/3.];

#analyze
analyze 1

############################################################
# NOW APPLY LOADING SEQUENCE AND ANALYZE (plastic)

# rezero time
setTime 0.0
#loadConst -time 0.0D0
wipeAnalysis

pattern UniformExcitation 1 1 -accel "Sine 0 1000 $period -factor $accMul"

# create the Analysis
constraints Transformation;
test NormDispIncr 1.e-6 25 0
algorithm Newton
numberer RCM
system ProfileSPD
rayleigh $massProportionalDamping 0.0 $stiffnessProportionalDamping 0.
integrator Newmark $gamma [expr pow($gamma+0.5, 2)/4]
analysis VariableTransient

recorder Element -ele 1 -time -file stress1.out -dT 0.01 material 2 stress
recorder Element -ele 1 -time -file strain1.out -dT 0.01 material 2 strain

set ok [analyze 2000 0.01]
puts $ok
puts "finished"
wipe;

Here is the error report.

WARNING: MaterialStageParameter::setDomain() - no effect with material tag 2
ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (1, 0)
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at load factor 1
OpenSees > analyze failed, returned: -3 error flag
WARNING: MaterialStageParameter::setDomain() - no effect with material tag 1
WARNING: MaterialStageParameter::setDomain() - no effect with material tag 2
MatParameter::setDomain(Domain *theDomain) - NO RESULT
ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (1, 0)
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at load factor 1
OpenSees > analyze failed, returned: -3 error flag
ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (1, 0)
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.01
OpenSees > analyze failed, returned: -3 error flag
-3
finished

I am in urgent need of your help!