Single quadrilateral element with a block2D command

[Carlota Gama]

Hello,

I'm trying to do the example 15 of the OpenSees Geotechnical Simulation Capabilities but with an uniform excitation at the base and create a mesh of quadrilateral elements. I've tryed every thing but still can run the problem with the block2D command, it gives this error: FourNodeQuad::FourNodeQuad -- failed to get a copy of material model. Does anyone know what could be wrong in the scrip?

set accMul 1 ;# acceleration multiplier
set massDen 1.50 ;# solid mass density
set fluidDen 1.0 ;# fluid mass density
set massProportionalDamping 0.0 ;
set stiffnessProportionalDamping 0.001 ;
set fangle 31.40 ;#friction angle
set ptangle 26.50 ;#phase transformation angle
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))] ;
set press 0.0 ;# isotropic consolidation pressure on quad element(s)
set deltaT 0.010 ;# time step for analysis
set numSteps 3600 ;# Number of analysis steps
set gamma 0.600 ;# Newmark integration parameter
set pi 3.1415926535 ;
set inclination 0 ;
set unitWeightX [expr ($massDen-$fluidDen)*9.81*sin($inclination/180.0*$pi)] ;# unit weight in X direction
set unitWeightY [expr -($massDen-$fluidDen)*9.81*cos($inclination/180.0*$pi)] ;# unit weight in Y direction

#############################################################

#create the ModelBuilder
model basic -ndm 2 -ndf 2

# define material and properties
nDMaterial PressureDependMultiYield 2 2 $massDen $G $B $fangle .1 80 0.5 \ $ptangle 0.17 0.4 10 10 0.015 1.0
nDMaterial FluidSolidPorous 1 2 2 2.2e+6

set Quad quad

set eleArgs "1 PlaneStrain2D 2"

set nx 2
set ny 2
set bn [expr $nx + 1 ]

block2D $nx $ny 1 1 $Quad $eleArgs {
1 0 0
2 1 0
3 1 1
4 0 1
}


# fix the base
fix 1 1 1
fix 2 1 1

#tie nodes 3 and 4
equalDOF 3 4 1 2

puts "Modelo construído"

#############################################################

# GRAVITY APPLICATION (elastic behavior)
# create the SOE, ConstraintHandler, Integrator, Algorithm and Numberer
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

puts "Gravidade aplicada."

# 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

wipeAnalysis

set GM "Series -dt 0.02 -filePath acelerograma1.txt -factor $accMul"

pattern UniformExcitation 1 1 -accel $GM

# 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 Node -file Data/disp_011.txt -time -node 1 2 3 4 -dof 1 -dT 0.01 disp
recorder Node -file Data/acce_011.txt -time -node 1 2 3 4 -dof 1 2 -dT 0.01 accel
recorder Element -ele 1 -time -file Data/stress1_011.txt -dT 0.01 material 1 stress
recorder Element -ele 1 -time -file Data/strain1_011.txt -dT 0.01 material 1 strain
recorder Element -ele 1 -time -file Data/stress3_011.txt -dT 0.01 material 3 stress
recorder Element -ele 1 -time -file Data/strain3_011.txt -dT 0.01 material 3 strain
recorder Element -ele 1 -time -file Data/press1_011.txt -dT 0.01 material 1 pressure
recorder Element -ele 1 -time -file Data/press3_011.txt -dT 0.01 material 3 pressure

#analyze
set startT [clock seconds]
analyze $numSteps $deltaT
set endT [clock seconds]
puts "Execution time: [expr $endT-$startT] seconds."

wipe #flush ouput stream

[Carlota Gama]

Hi again,

Maybe someone can tell me if the block2D command works with the PressureDependMultiYield material and with a PlainStrain2D problem, is this possible?