about nDMaterial

[hanjun009]

hi developers:
I am simulating infill wall of frame with nDMaterial in opensees. But the nDMaterial do not work well. when i use nDMaterial Template3Dep (codes followed), one error happen: "elasticisotropicmaterial::get copy--subclass responsibility", why? What Is Wrong?

wipe
model BasicBuilder -ndm 3 -ndf 6

#define material
set mattag 1

#nDMaterial ElasticIsotropic 200 2.0e10 0.25 2500
nDMaterial ElasticIsotropic 200 3e+10 2500
# Material "Material02":
# YieldSurfaceType <ParameterList>
set ys "-DP"
# PotentialSurfaceType <ParameterList>"
set ps "-DP +2.000000E-001"
# EvolutionLawType <ParameterList>"
set el1 "-Leq +1.000000E+000"
# InitialStressTensor Sxx Sxy Sxz Syx Syy Syz Szx Szy Szz
set sts " +1.000000E-002 +0.000000E+000 +0.000000E+000 +0.000000E+000 +1.000000E-002 +0.000000E+000 +0.000000E+000 +0.000000E+000 +1.000000E-002"
# Elasto-PlasticState <nt> ns sc1 sc2... <sts>"
set eps " -NOS 2 +2.000000E-001 +0.000000E+000 -stressp $sts"
# matTag ElmTag ys ps eps <el> <et>
nDMaterial Template3Dep $mattag 200 -YS $ys -PS $ps -EPS $eps -ELS1 $el1

#nDMaterial J2Plasticity $mattag 228e8 128e8 2000000 2000000 0 0
#nDMaterial ElasticIsotropic $mattag 3e+10 0.2
nDMaterial PlateFiber 2 $mattag
section PlateFiber 3 2 0.250
# ELEMENT
set Plate shell

block2D 2 2 100 10000 $Plate 3 {
1 -9 -5 3.9
2 -9 0 3.9
3 -9 0 7.2
4 -9 -5 7.2
}
print node
print ele

fix 100 1 1 1 1 1 1
fix 101 1 1 1 1 1 1
fix 102 1 1 1 1 1 1

fix 103 1 0 0 0 1 1
fix 104 1 0 0 0 1 1
fix 105 1 0 0 0 1 1
fix 106 1 0 0 0 1 1
fix 107 1 0 0 0 1 1
fix 108 1 0 0 0 1 1

pattern Plain 1 "Linear" {
load 108 0 20e2 0 0 0 0
}
#system SparseGeneral -piv
system BandGeneral
numberer RCM
constraints Plain
integrator LoadControl 0.1;
#integrator LoadControl 1.0
test EnergyIncr 1.0e-10 500 1
algorithm Newton
#algorithm KrylovNewton
analysis Static
analyze 10




then When i use J2Plasticity material, I also encountered two problems:1. "nDmaterial J2Plasticity $matTag $K $G $sig0 $sigInf $delta $H", how define the $delta and $H; 2. when i use following codes the opensees do not work,but i change the nDMaterial to ElasticIsotropic,it work well.

wipe
model BasicBuilder -ndm 3 -ndf 6

#define material
set mattag 1

nDMaterial J2Plasticity $mattag 228e8 128e8 2000000 2000000 0 0
#nDMaterial ElasticIsotropic $mattag 3e+10 0.2

nDMaterial PlateFiber 2 $mattag

section PlateFiber 3 2 0.250

# ELEMENT
set Plate shell
block2D 2 2 100 10000 $Plate 3 {
1 -9 -5 3.9
2 -9 0 3.9
3 -9 0 7.2
4 -9 -5 7.2
}
print node
print ele

fix 100 1 1 1 1 1 1
fix 101 1 1 1 1 1 1
fix 102 1 1 1 1 1 1

fix 103 1 0 0 0 1 1
fix 104 1 0 0 0 1 1
fix 105 1 0 0 0 1 1
fix 106 1 0 0 0 1 1
fix 107 1 0 0 0 1 1
fix 108 1 0 0 0 1 1

pattern Plain 1 "Linear" {
load 108 0 20e2 0 0 0 0
}


#system SparseGeneral -piv
system BandGeneral

numberer RCM

constraints Plain

integrator LoadControl 0.1;
#integrator LoadControl 1.0

test EnergyIncr 1.0e-10 500 1

algorithm Newton
#algorithm KrylovNewton
analysis Static
analyze 10

help please
thank you!

[hanjun009]

dear silvia and fmk
i have the matter about nonlinear nDMaterial Template3Dep and J2Plasticity. would you like to spend some time checking up my codes above? i also will be appreciated with you if you can give me one example about it.
thank you very much.
best wishes.

[jack4823119]

Dear Silvia:
I'd like to show great honor to you for your great contribution to earthquake engineering ! Whenever I get a problem and whatever it is, I will always get key help from you!
From a student in China

[fmk]

you cannot add Template3Dep to a shell element .. you are limited to J2 and ElasticIsotropic.

[hanjun009]

dear fmk
thank you for your valuable reply. But the J2 do not work as also.
i use J2Plasticity material, I also encountered two problems:1. "nDmaterial J2Plasticity $matTag $K $G $sig0 $sigInf $delta $H", how define the $delta and $H; 2. when i use following codes the opensees do not work,but i change the nDMaterial to ElasticIsotropic,it work well. what is wrong? please give me some valuable advice again.thank you again!


wipe
model BasicBuilder -ndm 3 -ndf 6

#define material
set mattag 1

nDMaterial J2Plasticity $mattag 228e8 128e8 2000000 2000000 0 0
#nDMaterial ElasticIsotropic $mattag 3e+10 0.2

nDMaterial PlateFiber 2 $mattag

section PlateFiber 3 2 0.250

# ELEMENT
set Plate shell
block2D 2 2 100 10000 $Plate 3 {
1 -9 -5 3.9
2 -9 0 3.9
3 -9 0 7.2
4 -9 -5 7.2
}
print node
print ele

fix 100 1 1 1 1 1 1
fix 101 1 1 1 1 1 1
fix 102 1 1 1 1 1 1

fix 103 1 0 0 0 1 1
fix 104 1 0 0 0 1 1
fix 105 1 0 0 0 1 1
fix 106 1 0 0 0 1 1
fix 107 1 0 0 0 1 1
fix 108 1 0 0 0 1 1

pattern Plain 1 "Linear" {
load 108 0 20e2 0 0 0 0
}


#system SparseGeneral -piv
system BandGeneral

numberer RCM

constraints Plain

integrator LoadControl 0.1;
#integrator LoadControl 1.0

test EnergyIncr 1.0e-10 500 1

algorithm Newton
#algorithm KrylovNewton
analysis Static
analyze 10

[fmk]

the example you posted WORKS with OpenSees2.0.0.

[hanjun009]

dear fmk
thank you very much for your help! Can you also give me some advice about the values of the parameter '$delta' and '$H' in the command 'nDmaterial J2Plasticity $matTag $K $G $sig0 $sigInf $delta $H'?
thank you very much.
best wishes.