nDMaterial Template3Dep

[arbarbosa]

I am trying to run example 8.1 using the nDmaterial Template3Dep. However, when reading the nDmaterial Template3Dep the run stops as if waiting for more parameters... (I am using OpenSees v1.7.3)

I am using the example for material model ConcreteEPmodelDP.tcl . I had to make a change since it had missing the a old material tag, but it seems I am still missing something that I cannot find.

Please find attached the files I am using. Help is appreciated.
Thanks,
Andre

# DPtest.tcl
#
# Adapted from OpenSees Example 8.1
# OpenSees Primer
#
# Units: kN, m, sec
# ----------------------------
# Start of model generation
# ----------------------------
set ConcreteMaterial 1
# Create ModelBuilder with 3 dimensions and 3 DOF/node
model basic -ndm 3 -ndf 3
# create the material
nDMaterial ElasticIsotropic3D 1 30000 0.20 2.4
if {$ConcreteMaterial==1} {
puts "Trying to read Template3Dep..."
source ConcreteEPmodelDP.tcl
puts "Finished reading Template3Dep..."
}
#
# Define geometry
# ---------------
# define some parameters
set eleArgs "1"
set element stdBrick
set nz 10
set nx 5
set ny 2
set nn [expr ($nz+1)*($nx+1)*($ny+1) ]
# mesh generation
block3D $nx $ny $nz 1 1 $element $eleArgs {
1 -0.35 -0.1 0
2 0.35 -0.1 0
3 0.35 0.1 0
4 -0.35 0.1 0
5 -0.35 -0.1 1.70
6 0.35 -0.1 1.70
7 0.35 0.1 1.70
8 -0.35 0.1 1.70
}
set load 0.1
# Constant point load
pattern Plain 1 Linear {
load $nn 0.0 $load 0.0
}
# boundary conditions
fixZ 0.0 1 1 1
# --------------------------------------------------------------------
# Start of static analysis (creation of the analysis & analysis itself)
# --------------------------------------------------------------------
# Load control with variable load steps
# init Jd min max
integrator LoadControl 1.0 1 10 100.0
# Convergence test
# tolerance maxIter displayCode
test NormUnbalance 1.0e-10 20 1
# Solution algorithm
algorithm Newton
# DOF numberer
numberer RCM
# Cosntraint handler
constraints Plain
# System of equations solver
system ProfileSPD
# Analysis for gravity load
analysis Static
# Perform the analysis
analyze 5
# --------------------------
# End of static analysis
# --------------------------
# ----------------------------
# Start of recorder generation
# ----------------------------
recorder Node -file Node.out -time -node $nn -dof 2 disp
recorder plot Node.out CenterNodeDisp 625 10 625 450 -columns 1 2
recorder display ShakingBeam 0 0 300 300 -wipe
prp -100 100 120.5
vup 0 0 1
display 1 0 1
# --------------------------
# End of recorder generation
# --------------------------
# ---------------------------------------
# Create and Perform the dynamic analysis
# ---------------------------------------
# Remove the static analysis & reset the time to 0.0
wipeAnalysis
setTime 0.0
# Now remove the loads and let the beam vibrate
remove loadPattern 1
# add some mass proportional damping
# rayleigh 0.01 0.0 0.0 0.0
# Create the transient analysis
test EnergyIncr 1.0e-10 20 1
algorithm Newton
numberer RCM
constraints Plain
integrator Newmark 0.5 0.25
analysis Transient
# Perform the transient analysis (20 sec)
# numSteps dt
analyze 100 0.01




#
# ConcreteEPmodelDP.tcl - Boris Jeremic Material Examples
#
# Template elasto-plastic material model for concrete, Drucker-Prager yield surface
# Young's modulus E=30000 kpa,
# Poisson's ratio v=0.18,
# Friction angle = 40deg
# Cohesion k = 0.0

# Aug. 12, 2002 Feng Xiong
# Boris Jeremic (@ucdavis.edu)

# Yield surface
set DPys "-DP"

# Potential surface
set DPps "-DP 0.1"

# Scalar evolution law: linear hardening coef=1.0
set ES1 "-Leq 1.10"

# Initial stress
set stressp "0.1 0 0 0 0.1 0 0 0 0.1"

# EPState
#------------E--------Eo--------v----rho-----------------alpha----k
set EPS "30000.0 30000.0 0.18 2.4 -NOD 0 -Nos 2 0.315 0.0 -stressp $stressp"

# Create nDMaterial using Template Elastic-Plastic Model
nDMaterial Template3Dep 2 1 -YS $DPys -PS $DPps -EPS $EPS -ELS1 $ES1

[Boris]

I am not sure what 1.7.3 has inside, you need to talk to Frank, I follow version by date... To be honest we have not used Template3Dep for some time now as we have a NewTemplate3Dep... (see my lecture notes for commands).

Also the main developer of that new system is finishing his thesis so there will be a lot of examples availablel in report that is comming up soon (couple of months...).

Boris

[arbarbosa]

Thank you for the reply.

I just ran (10/26/2006) cvs and compiled the opensees to have the updated files. The Template3D still does not work for me.

I also tried the NewTemplate3Dep and I get the following message before the run stops

"Warning: TclNewTemplate3Dep - unknown keyword/command :"

Any ideas so I can try to solve this?

Thank you again.
Andre
--------------------------
Note: I used the following example as shown in 12.3.8.3 of http://sokocalo.engr.ucdavis.edu/~jerem ... ode14.html

## Example 03, for New template3Dep material model
## Drucker-Prager, Perfectly-Plastic

set rho 2.4
set E 3.0e4
set v 0.2
set a 0.3
set k 0.0
set d 0.1
set MC "$rho $E $v $a $k $d"
set MP "MaterialConstant 6 $MC"
set ES "Isotropic 2 3"
set YF "Drucker-Prager 0 4 0 5"
set PF "Drucker-Prager 0 6"

nDMaterial NewTemplate3Dep 1 -MaterialParameter $MP \
-ElasticState $ES \
-YieldFunction $YF \
-PlasticFlow $PF

[WAELL]

Hello Andre:

It seams that we have the same problems with "template3Dep" and "Newtemplate3Dep" commands. So….if you found any solution, please e-mail me and I will do the same for you.
My e-mail address is waell_ab@yahoo.com

Thanks
Wael