*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 , 0.0, , 0.0 *NGEN, NSET=LEFT 1,,1 *NCOPY, CHANGENUMBER=, OLDSET=LEFT, NEWSET=RIGHT, SHIFT , 0.0, 0.0 0.0,0.0,0.0,0.0,1.0,0.0,0.0 *NFILL, NSET=BOTNODES LEFT, RIGHT, , *NCOPY, CHANGENUMBER=, OLDSET=BOTNODES, NEWSET=TOPNODES, SHIFT 0.0,0.0, 0.0,0.0,0.0,0.0,1.0,0.0,0.0 *NFILL, NSET=ALLNODES BOTNODES, TOPNODES, , ** ** Elements ** *ELEMENT, TYPE=C3D8, ELSET=BOT 1,1,,,2,,,, *ELGEN, ELSET=EALL 1,,1,1,,,,,, *ELSET, ELSET=HARD, gen 1,,1 *ELSET, ELSET=SOFT, gen ,,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