Error Showing while using nDMaterial ElasticIsotropic comman

[manishjindal4]

I have made a program for brickwall problem and I have used nDMaterial ElasticIsotropic for defining material properties for Brick27N element.
It is showing error that
ElasticIsotropicMaterial::getCopy -- subclass responsibility

I am attaching here the program code with this post. please help me in this problem.
#model dimension and degree of freedom
model basic -ndm 3 -ndf 3

#brick material
nDMaterial ElasticIsotropic 1 3750000000 0.1

#mortar material
nDMaterial ElasticIsotropic 2 545000000 0.19

set nob 12
#nob->number of bricks
set nobl 35
#nobl->number of brick layers
set noer 100
#number of elements in a row
set tnonl [expr (2*$noer+1)*12]
#total number of nodes in a layer(brick+mortar)::(2*$noer+1)*12::2412

for { set k 1 } { $k <= $nobl } { incr k } {
for { set j 1 } { $j <= [expr (2*$noer+1)] } { incr j } {
node [expr ($k-1)*$tnonl+($j-1)*9+1] 0 [expr ($j-1)*15] [expr ($k-1)*0.085]
node [expr ($k-1)*$tnonl+($j-1)*9+2] -0.125 [expr ($j-1)*15] [expr ($k-1)*0.085]
node [expr ($k-1)*$tnonl+($j-1)*9+3] -0.125 [expr ($j-1)*15] [expr ($k-1)*0.085+0.075]
node [expr ($k-1)*$tnonl+($j-1)*9+4] 0 [expr ($j-1)*15] [expr ($k-1)*0.085+0.075]
node [expr ($k-1)*$tnonl+($j-1)*9+5] -0.0625 [expr ($j-1)*15] [expr ($k-1)*0.085]
node [expr ($k-1)*$tnonl+($j-1)*9+6] -0.125 [expr ($j-1)*15] [expr ($k-1)*0.085+0.0375]
node [expr ($k-1)*$tnonl+($j-1)*9+7] -0.0625 [expr ($j-1)*15] [expr ($k-1)*0.085+0.075]
node [expr ($k-1)*$tnonl+($j-1)*9+8] 0 [expr ($j-1)*15] [expr ($k-1)*0.085+0.0375]
node [expr ($k-1)*$tnonl+($j-1)*9+9] -0.0625 [expr ($j-1)*15] [expr ($k-1)*0.085+0.0375]


node [expr ($k-1)*$tnonl+1810+($j-1)*3] -0.125 [expr ($j-1)*15] [expr ($k-1)*0.085+0.08]
node [expr ($k-1)*$tnonl+1811+($j-1)*3] 0 [expr ($j-1)*15] [expr ($k-1)*0.085+0.08]
node [expr ($k-1)*$tnonl+1812+($j-1)*3] -0.0625 [expr ($j-1)*15] [expr ($k-1)*0.085+0.08]
}
}
for { set j 1 } { $j <= [expr (2*$noer+1)] } { incr j } {
node [expr 35*$tnonl+($j-1)*9+3] -0.125 [expr ($j-1)*15] [expr 35*0.085+0.075]
node [expr 35*$tnonl+($j-1)*9+4] 0 [expr ($j-1)*15] [expr 35*0.085+0.075]
node [expr 35*$tnonl+($j-1)*9+7] -0.0625 [expr ($j-1)*15] [expr 35*0.085+0.075]
}

for { set k 1 } { $k <= $nobl } { incr k } {
for { set j 1 } { $j <= $noer } { incr j } {

set a1 [expr ($j-1)*18+($k-1)*$tnonl+19]
set a2 [expr ($j-1)*18+($k-1)*$tnonl+20]
set a3 [expr ($j-1)*18+($k-1)*$tnonl+2]
set a4 [expr ($j-1)*18+($k-1)*$tnonl+1]
set a5 [expr ($j-1)*18+($k-1)*$tnonl+22]
set a6 [expr ($j-1)*18+($k-1)*$tnonl+21]
set a7 [expr ($j-1)*18+($k-1)*$tnonl+3]
set a8 [expr ($j-1)*18+($k-1)*$tnonl+4]
set a9 [expr ($j-1)*18+($k-1)*$tnonl+23]
set a10 [expr ($j-1)*18+($k-1)*$tnonl+11]
set a11 [expr ($j-1)*18+($k-1)*$tnonl+5]
set a12 [expr ($j-1)*18+($k-1)*$tnonl+10]
set a13 [expr ($j-1)*18+($k-1)*$tnonl+25]
set a14 [expr ($j-1)*18+($k-1)*$tnonl+12]
set a15 [expr ($j-1)*18+($k-1)*$tnonl+7]
set a16 [expr ($j-1)*18+($k-1)*$tnonl+13]
set a17 [expr ($j-1)*18+($k-1)*$tnonl+26]
set a18 [expr ($j-1)*18+($k-1)*$tnonl+24]
set a19 [expr ($j-1)*18+($k-1)*$tnonl+6]
set a20 [expr ($j-1)*18+($k-1)*$tnonl+8]
set a21 [expr ($j-1)*18+($k-1)*$tnonl+27]
set a22 [expr ($j-1)*18+($k-1)*$tnonl+15]
set a23 [expr ($j-1)*18+($k-1)*$tnonl+9]
set a24 [expr ($j-1)*18+($k-1)*$tnonl+17]
set a25 [expr ($j-1)*18+($k-1)*$tnonl+16]
set a26 [expr ($j-1)*18+($k-1)*$tnonl+14]
set a27 [expr ($j-1)*18+($k-1)*$tnonl+18]

set b1 [expr ($j-1)*18+($k-1)*$tnonl+22]
set b2 [expr ($j-1)*18+($k-1)*$tnonl+21]
set b3 [expr ($j-1)*18+($k-1)*$tnonl+3]
set b4 [expr ($j-1)*18+($k-1)*$tnonl+4]
set b5 [expr ($j-1)*18+$k*$tnonl+19]
set b6 [expr ($j-1)*18+$k*$tnonl+20]
set b7 [expr ($j-1)*18+$k*$tnonl+2]
set b8 [expr ($j-1)*18+$k*$tnonl+1]
set b9 [expr ($j-1)*18+($k-1)*$tnonl+25]
set b10 [expr ($j-1)*18+($k-1)*$tnonl+12]
set b11 [expr ($j-1)*18+($k-1)*$tnonl+7]
set b12 [expr ($j-1)*18+($k-1)*$tnonl+13]
set b13 [expr ($j-1)*18+$k*$tnonl+23]
set b14 [expr ($j-1)*18+$k*$tnonl+11]
set b15 [expr ($j-1)*18+$k*$tnonl+5]
set b16 [expr ($j-1)*18+$k*$tnonl+10]
set b17 [expr ($j-1)*9+($k-1)*$tnonl+1809+8]
set b18 [expr ($j-1)*9+($k-1)*$tnonl+1809+7]
set b19 [expr ($j-1)*9+($k-1)*$tnonl+1809+1]
set b20 [expr ($j-1)*9+($k-1)*$tnonl+1809+2]
set b21 [expr ($j-1)*9+($k-1)*$tnonl+1809+9]
set b22 [expr ($j-1)*9+($k-1)*$tnonl+1809+4]
set b23 [expr ($j-1)*9+($k-1)*$tnonl+1809+3]
set b24 [expr ($j-1)*9+($k-1)*$tnonl+1809+5]
set b25 [expr ($j-1)*18+$k*$tnonl+14]
set b26 [expr ($j-1)*18+($k-1)*$tnonl+16]
set b27 [expr ($j-1)*9+($k-1)*$tnonl+1809+6]

element Brick27N [expr $j+2*$noer*$k] $a1 $a2 $a3 $a4 $a5 $a6 $a7 $a8 $a9 $a10 $a11 $a12 $a13 $a14 $a15 $a16 $a17 $a18 $a19 $a20 $a21 $a22 $a23 $a24 $a25 $a26 $a27 1 0 0 0 2000
element Brick27N [expr $j+(2*$k-1)*$noer] $b1 $b2 $b3 $b4 $b5 $b6 $b7 $b8 $b9 $b10 $b11 $b12 $b13 $b14 $b15 $b16 $b17 $b18 $b19 $b20 $b21 $b22 $b23 $b24 $b25 $b26 $b27 2 0 0 0 2200
}
}

for { set j 1 } { $j <= [expr (2*$noer+1)] } { incr j } {
fix [expr ($j-1)*9+1] 1 1 1
fix [expr ($j-1)*9+2] 1 1 1
fix [expr ($j-1)*9+5] 1 1 1
}

for { set k 1 } { $k <= [expr ($tnonl*$nobl } { incr k } {
recorder Node -file [format "output/node%i.out" $k] -node $k -dof 1 2 3 disp
}
pattern Plain 1 Linear {
load 100 2000. 0. 0. 0. 0. 0.
load 101 2000. 0. 0. 0. 0. 0.
load 102 2000. 0. 0. 0. 0. 0.
}
recorder display model 120 120 800 800 -wipe
prp 0. 3000 3000
vup 0 0 1
display 1 3 1000
constraints Plain
numberer RCM
system UmfPack
test NormDispIncr 1.0e-8 10
integrator LoadControl .002
algorithm Linear
analysis Static
analyze 500


Thanking you for giving your valuable time to this


With Regards

Manish Jindal
Final Year Undergraduate Student
School of Infrastructure
Indian Institute of Technology Bhubaneswar
India
+91-9040506605
manish.iitbbs@gmail.com
manishjindal4@gmail.com

[jdgome]

Hi Manish. Have you solved this trouble? I have similar issue

[fmk]

the 27 noded brick was removed years ago.

[jdgome]

Mr McKenna, Do you have a 2.2 Opensees version available to share? I need u-p-U brick element developed by prof. Boris Jeremic. And current version does not have it. My contact email is jgomezz@eafit.edu.co