Problem with 1.7.5 version of OpenSees in updating material

Lanka · Post by **Lanka** » Tue Feb 05, 2008 10:22 am

Hi Silvia,

I recently installed 1.7.5 version of OpenSees. I am getting following warning for updateMaterialStage commad of PDMY material with 1.7.5 version.

WARNING: updateMaterialStage - MaterialStageParameter:: setDomain() - no effect with material tag $ID

My script worked fine with 1.7.3 version of OpenSees. I have attached the script, Could you check this?

Thanks in advance.

-Lanka.

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

#============================================================
#
# Unit: m, s, kN
#
#============================================================

wipe

file mkdir ./data/rawOut
file mkdir ./data/rawOut_gravity

#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#
# Soil Gravity
#
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

##############################################################################################[1]
model BasicBuilder -ndm 2 -ndf 2

##############################################################################################[2]
# ==============================================
# (2.1) nodes (soil,py,dashpot,boundary,zeronode)
# ==============================================
node 1 0.0000000000 0.0000000000
node 2 1.0000000000 0.0000000000
node 3 0.0000000000 1.1570000000
node 4 1.0000000000 1.1570000000
node 5 0.0000000000 2.3140000000
node 6 1.0000000000 2.3140000000
node 7 0.0000000000 3.4710000000
node 8 1.0000000000 3.4710000000
node 9 0.0000000000 4.6280000000
node 10 1.0000000000 4.6280000000
node 11 0.0000000000 5.7850000000
node 12 1.0000000000 5.7850000000
node 13 0.0000000000 6.9420000000
node 14 1.0000000000 6.9420000000
node 15 0.0000000000 8.0990000000
node 16 1.0000000000 8.0990000000
node 17 0.0000000000 9.2560000000
node 18 1.0000000000 9.2560000000
node 19 0.0000000000 10.4130000000
node 20 1.0000000000 10.4130000000
node 21 0.0000000000 11.5700000000
node 22 1.0000000000 11.5700000000
node 23 0.0000000000 12.7270000000
node 24 1.0000000000 12.7270000000
node 25 0.0000000000 13.8840000000
node 26 1.0000000000 13.8840000000
node 27 0.0000000000 15.0410000000
node 28 1.0000000000 15.0410000000
node 29 0.0000000000 16.1980000000
node 30 1.0000000000 16.1980000000
node 31 0.0000000000 17.3550000000
node 32 1.0000000000 17.3550000000
node 33 0.0000000000 18.5120000000
node 34 1.0000000000 18.5120000000
node 35 0.0000000000 19.6690000000
node 36 1.0000000000 19.6690000000
node 37 0.0000000000 20.8260000000
node 38 1.0000000000 20.8260000000
node 39 0.0000000000 21.9830000000
node 40 1.0000000000 21.9830000000
node 41 0.0000000000 23.1400000000
node 42 1.0000000000 23.1400000000
node 43 0.0000000000 24.2970000000
node 44 1.0000000000 24.2970000000
node 45 0.0000000000 25.4540000000
node 46 1.0000000000 25.4540000000
node 47 0.0000000000 26.6110000000
node 48 1.0000000000 26.6110000000
node 49 0.0000000000 27.7680000000
node 50 1.0000000000 27.7680000000
# ==============================================
# (2.2) mass (soil,py,dashpot,boundary,zeronode)
# ==============================================
##############################################################################################[3]
# ==============================================
# (3.1) define nDmaterial for soil
# ==============================================

nDMaterial PressureDependMultiYield 51 2 2.158 57500 172500 37.000 0.100 80.000 0.500 27.000 0.050 0.600 3.000 5.000 0.003 1.000
updateMaterialStage -material 51 -stage 0
nDMaterial PressureDependMultiYield 52 2 2.158 57500 172500 37.000 0.100 80.000 0.500 27.000 0.050 0.600 3.000 5.000 0.003 1.000
updateMaterialStage -material 52 -stage 0
nDMaterial PressureDependMultiYield 53 2 2.158 57500 172500 37.000 0.100 80.000 0.500 27.000 0.050 0.600 3.000 5.000 0.003 1.000
updateMaterialStage -material 53 -stage 0

##############################################################################################[6]

###############################################################################
# ==============================================
# (3.3) define quadElement for soil (max-10 soils)
# ==============================================

# for normal gravity analysis
set gY -9.81
set gX 0.0
set p 0.0
set eleThick 1.0
set eleDensity 0.0

#---------------
set matID1 51
set matID2 52
set matID3 53
set matID4 54
set matID5 55
set matID6 56
set matID7 57
set matID8 58
set matID9 59
set matID10 60
#---------------

element quad 1 34 36 35 33 $eleThick PlaneStrain $matID1 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 2 36 38 37 35 $eleThick PlaneStrain $matID1 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 3 38 40 39 37 $eleThick PlaneStrain $matID1 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 4 40 42 41 39 $eleThick PlaneStrain $matID1 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 5 42 44 43 41 $eleThick PlaneStrain $matID1 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 6 44 46 45 43 $eleThick PlaneStrain $matID1 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 7 46 48 47 45 $eleThick PlaneStrain $matID1 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 8 48 50 49 47 $eleThick PlaneStrain $matID1 $p $eleDensity $gX [expr $gY/(1/2.158)]

element quad 9 18 20 19 17 $eleThick PlaneStrain $matID2 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 10 20 22 21 19 $eleThick PlaneStrain $matID2 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 11 22 24 23 21 $eleThick PlaneStrain $matID2 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 12 24 26 25 23 $eleThick PlaneStrain $matID2 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 13 26 28 27 25 $eleThick PlaneStrain $matID2 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 14 28 30 29 27 $eleThick PlaneStrain $matID2 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 15 30 32 31 29 $eleThick PlaneStrain $matID2 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 16 32 34 33 31 $eleThick PlaneStrain $matID2 $p $eleDensity $gX [expr $gY/(1/2.158)]

element quad 17 15 16 18 17 $eleThick PlaneStrain $matID3 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 18 13 14 16 15 $eleThick PlaneStrain $matID3 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 19 11 12 14 13 $eleThick PlaneStrain $matID3 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 20 9 10 12 11 $eleThick PlaneStrain $matID3 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 21 7 8 10 9 $eleThick PlaneStrain $matID3 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 22 5 6 8 7 $eleThick PlaneStrain $matID3 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 23 3 4 6 5 $eleThick PlaneStrain $matID3 $p $eleDensity $gX [expr $gY/(1/2.158)]
element quad 24 1 2 4 3 $eleThick PlaneStrain $matID3 $p $eleDensity $gX [expr $gY/(1/2.158)]
#-----------------------------------------------

##############################################################################################[4]

##############################################################################################[5]

##############################################################################################[7]

# ===============================================
# (7.1) define equal DOF (soil to soil)
# ==============================================
equalDOF 1 2 1 2
equalDOF 3 4 1 2
equalDOF 5 6 1 2
equalDOF 7 8 1 2
equalDOF 9 10 1 2
equalDOF 11 12 1 2
equalDOF 13 14 1 2
equalDOF 15 16 1 2
equalDOF 17 18 1 2
equalDOF 19 20 1 2
equalDOF 21 22 1 2
equalDOF 23 24 1 2
equalDOF 25 26 1 2
equalDOF 27 28 1 2
equalDOF 29 30 1 2
equalDOF 31 32 1 2
equalDOF 33 34 1 2
equalDOF 35 36 1 2
equalDOF 37 38 1 2
equalDOF 39 40 1 2
equalDOF 41 42 1 2
equalDOF 43 44 1 2
equalDOF 45 46 1 2
equalDOF 47 48 1 2
equalDOF 49 50 1 2

# ==============================================
# (8.1) define the fixity
# ==============================================
fix 1 1 1
fix 2 1 1

##############################################################################################[7]
# gravity application (elastic)
##############################################################################################

# ===============================
# Recorder (gravity application)
# ===============================

set f1 rawOut_Gravity

recorder Node -file ./data/$f1/displacement_gStep.out -time -node 1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47 49 -dof 2 disp
recorder Node -file ./data/$f1/acceleration_gStep.out -time -node all -dof 1 2 accel

##############################################################################################[7]
# gravity application (elastic)
##############################################################################################
# ==============================================
# (9.1) create the analysis method
# ==============================================

system ProfileSPD
test NormDispIncr 1e-4 5 0
algorithm Newton
constraints Transformation
integrator LoadControl 1 1 1 1
numberer Plain
analysis Static
# ==============================================
# (9.2) analyize
# ==============================================
analyze 2

##############################################################################################[10]

# ==============================================
# (10.1) updateMaterialStage
# ==============================================
updateMaterialStage -material 51 -stage 1
updateMaterialStage -material 52 -stage 1
updateMaterialStage -material 53 -stage 1

puts "==============================================="
puts "stage I. (soil gravity) - done!"

##############################################################################################[18]
# ===============================
# (18.1) define the display
# ===============================
recorder display "OpenSees Real Time" 10 10 700 700 -wipe
prp 0 0 100
#vrp 1 1 1
vup 0 1 0
vpn 0 0 1
display 1 1 200

#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#
# Earthquake
#
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

puts " ................................... "
puts " ......... earthquake ! .......... "
puts " ................................... "

wipeAnalysis
loadConst -time 0.0

set accelSeries1 "Path -filePath ./baseMotions/MIL03_09_base_ax.txt -dt 0.0127 -factor 9.81"
pattern UniformExcitation 2 1 -accel $accelSeries1

#set accelSeries2 "Path -filePath ./baseMotions/MIL03_02_base_az.txt -dt 0.0127 -factor 9.81"
#pattern UniformExcitation 3 2 -accel $accelSeries2

# ==============================================
# (9.1) recorder (earthquake step)
# ==============================================
remove recorders

# ===============================
# (18.1) define the display
# ===============================
recorder display "OpenSees Real Time" 10 10 700 700 -wipe
prp 0 0 100
#vrp 1 1 1
vup 0 1 0
vpn 0 0 1
display 1 1 200

set f1 rawOut

recorder Node -file ./data/$f1/displacement.txt -time -node 11 19 29 39 45 49 -dof 1 2 disp
recorder Node -file ./data/$f1/acceleration.txt -time -node 11 19 29 39 45 49 -dof 1 2 accel
recorder Element -ele 8 -file ./data/$f1/modulus1.out -time material 1 backbone 17
recorder Element -ele 8 -file ./data/$f1/modulus2.out -time material 2 backbone 17
recorder Element -ele 8 -file ./data/$f1/modulus3.out -time material 3 backbone 17
recorder Element -ele 8 -file ./data/$f1/modulus4.out -time material 4 backbone 17
recorder Element -file ./data/$f1/stress2.out -time -eleRange 1 0 material 2 stress
recorder Element -file ./data/$f1/stress3.out -time -eleRange 1 0 material 3 stress
recorder Element -file ./data/$f1/stress4.out -time -eleRange 1 0 material 4 stress
recorder Element -file ./data/$f1/strain1.out -time -eleRange 1 0 material 1 strain
recorder Element -file ./data/$f1/strain2.out -time -eleRange 1 0 material 2 strain
recorder Element -file ./data/$f1/strain3.out -time -eleRange 1 0 material 3 strain
recorder Element -file ./data/$f1/strain4.out -time -eleRange 1 0 material 4 strain
recorder Element -file ./data/$f1/pressure1.out -time -eleRange 1 0 material 1 pressure
recorder Element -file ./data/$f1/pressure2.out -time -eleRange 1 0 material 2 pressure
recorder Element -file ./data/$f1/pressure3.out -time -eleRange 1 0 material 3 pressure
recorder Element -file ./data/$f1/pressure4.out -time -eleRange 1 0 material 4 pressure

# ===============================
# (19.1) analysis method
# ===============================
constraints Penalty 1.0e12 1.0e12
test NormDispIncr 1e-3 50 0
algorithm Newton
numberer Plain
system ProfileSPD
integrator Newmark 0.6 0.3025 0.0 0.00 0.00 0.002
analysis Transient

set startT [clock seconds]
analyze 4096 0.0127 ;# analyze (plasitc)
set endT [clock seconds]
puts "stage III. (earthquake) - done!"
puts "==========================================="
puts "Execution time: [expr $endT-$startT] seconds."

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

#################################################################################
# ===============================
# (20.1) wipe all
# ===============================
wipe ;#flush ouput stream

Post by **fmk** » Tue Feb 05, 2008 11:20 am

place the first updates after the element definitions. b.t.w. they are not needed, as by default the stage is set to 0.

Lanka · Post by **Lanka** » Tue Feb 05, 2008 11:44 am

thanks frank.