*HEADING
static analysis of an elastic plate
*PARAMETER
** number of nodes in x1-direction
nnx     = 21
** number of nodes in x2-direction
nny    = nnx
** number of nodes in x3-direction
nnz    = nnx
** number of elements in x1-direction
nex    = nnx-1
** number of elements in x2-direction
ney    = nny-1
** number of elements in x3-direction
nez    = nnz-1
** spatial dimensions
lx = 40.0
ly  = 40.0
lz  = 40.0
**
nyp1 = nny+1
nyp2 = nny+2
nn3 = nnx*nny-nny
nnxnny = nnx*nny
nexney = nex*ney
nn3p1 = nn3+1
nn4 = nnx*nny*nnz - nnx*nny
nntot = nnx*nny*nnz
n5 = nnx*nny+1
n6 = nnx*nny+nny+1
n7 = n6+1
n8 = n5+1
nel = nex*ney*nez
midel = nel / 2
midelp1 = midel+1
**
** Nodal coordinates
**
*NODE, NSET=N1
1, 0.0, 0.0, 0.0
*NODE, NSET=N2
<nny>, 0.0, <ly>, 0.0
*NGEN, NSET=LEFT
1,<nny>,1
*NCOPY, CHANGENUMBER=<nn3>, OLDSET=LEFT, NEWSET=RIGHT, SHIFT
<lx>, 0.0, 0.0
0.0,0.0,0.0,0.0,1.0,0.0,0.0
*NFILL, NSET=BOTNODES
LEFT, RIGHT, <nex>, <nny>
*NCOPY, CHANGENUMBER=<nn4>, OLDSET=BOTNODES, NEWSET=TOPNODES, SHIFT
0.0,0.0,<lz>
0.0,0.0,0.0,0.0,1.0,0.0,0.0
*NFILL, NSET=ALLNODES
BOTNODES, TOPNODES, <nez>, <nnxnny>
**
** Elements
**
*ELEMENT, TYPE=C3D8, ELSET=BOT
1,1,<nyp1>,<nyp2>,2,<n5>,<n6>,<n7>,<n8>
*ELGEN, ELSET=EALL
1,<ney>,1,1,<nex>,<nny>,<ney>,<nez>,<nnxnny>,<nexney>
*ELSET, ELSET=HARD, gen
1,<midel>,1
*ELSET, ELSET=SOFT, gen
<midelp1>,<nel>,1
**
** Physical and material properties
**
*SOLID SECTION, ELSET=HARD, MATERIAL=STEEL
1.0
*SOLID SECTION, ELSET=SOFT, MATERIAL=SOFT
1.0
*MATERIAL, NAME=STEEL
*ELASTIC
1.E11, 0.3
*MATERIAL, NAME=SOFT
*ELASTIC
1.E10, 0.3
*STEP, PERT, solver=ddm
Static analysis
*STATIC
*solver controls
1.e-6,,10
*BOUNDARY
BOTNODES, ENCASTRE
*CLOAD
TOPNODES,3,-1.e6
*output, field
*node output
U
*element output
S
*END STEP