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!
PressureDependMultiYield Material
Moderators: silvia, selimgunay, Moderators