Errors with constant strain triangular Tri31 Element

[aryanrezaeirad]

Hi,
I am trying to use Tri31 Element in a linear static analysis, however, I get error even in a very simple case.
Three nodes should be defined for element in a counter-clockwise direction. (For example nodes: 1 2 3). After building the element, it gives the following error:
"setDomain() element 1 does not have 2 or 3 DOF at node 1" (node 1 which is the first defined node of element)

The following is a simple example. Thanks for the help in advance.






wipe all

# unit N, mm, second

model BasicBuilder -ndm 3 -ndf 6;

recorder Node -file Disp.out -node 3 -dof 1 2 3 disp

node 1 2710.18 375 112.935
node 2 2498.26 0 2.64421
node 3 2103.94 0 760.3

fix 1 1 1 1 0 0 0
fix 2 1 1 1 0 0 0

set Ex 10000.0; # modulus of elasticity
set Ey 2400.0
set Ez 130.0
set vxy 600.0
set vyz 22.0; # Poisson
set vzx 60.0
set Gxy 0.09
set Gyz 0.68; # modulus of shear
set Gzx 0.85
set rho [expr 51.0/100000000.0]

nDMaterial ElasticOrthotropic 1 $Ex $Ey $Ez $vxy $vyz $vzx $Gxy $Gyz $Gzx $rho
nDMaterial PlaneStress 2 1

element tri31 [expr 1000+10+1] 1 2 3 21.0 PlaneStress 2 1.09211606409936E-09 510.0 1.0 1.0

set gravityload 0.000829509813904852
set fi [expr (1/3.)*$gravityload]


pattern Plain 200 Linear {

load 3 0.0 0.0 [expr -($gravityload)] 0.0 0.0 0.0

}

set Dmax 40.0
set Dincr 0.01
set Nsteps [expr int($Dmax/$Dincr)];



constraints Plain;
numberer RCM;
system BandGeneral;
test NormUnbalance 1.0e-6 400;
algorithm Newton;
integrator DisplacementControl 3 3 $Dincr;
analysis Static;
set Nsteps [expr int($Dmax/$Dincr)];
set ok [analyze $Nsteps];
puts "Pushover complete";

[fmk]

can you not read the error message! .. your model has 6 dof per node!!!!!!!