*HEADING EC3FEFS2: PRESSURE AND HYDROSTATIC PRESSURE FOR 3-D SOLIDS [C3D15E] *PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES *RESTART,WRITE *NODE,NSET=NALL 1, 0.0, 0.0, 0.0 2, 2.0, 0.0, 0.0 3, 2.0, 2.0, 0.0 4, 0.0, 2.0, 0.0 5, 0.0, 0.0, 1.0 6, 2.0, 0.0, 1.0 7, 2.0, 2.0, 1.0 8, 0.0, 2.0, 1.0 9, 1.0, 0.0, 0.0 10, 2.0, 1.0, 0.0 11, 1.0, 2.0, 0.0 12, 0.0, 1.0, 0.0 13, 1.0, 0.0, 1.0 14, 2.0, 1.0, 1.0 15, 1.0, 2.0, 1.0 16, 0.0, 1.0, 1.0 17, 0.0, 0.0, 0.5 18, 2.0, 0.0, 0.5 19, 2.0, 2.0, 0.5 20, 0.0, 2.0, 0.5 21, 1.0, 1.0, 0.0 22, 1.0, 1.0, 1.0 *ELEMENT,TYPE=C3D15E, ELSET=EALL 1, 1,2,4,5,6,8,9,21,12,13,22,16,17,18,20 2, 2,3,4,6,7,8,10,11,21,14,15,22,18,19,20 *SOLID SECTION,MATERIAL=A1, ELSET=EALL *surface,NAME=SIDE 1,S2 2,S2 *MATERIAL,NAME=A1 *ELASTIC,TYPE=ISOTROPIC 30.0E6,0.3 *density 200 *PIEZOELECTRIC,TYPE=S 0.,0.,0.,0.,0.,0.,0.,0. 0.,0.,0.,0.,0.,0.,0.,0. 0.,0. *DIELECTRIC,TYPE=ISO 1.0E-3, *BOUNDARY 1,1,3 2,2 4,3 5,1 1,9 2,9 3,9 4,9 9,9 10,9 11,9 12,9 21,9 *STEP,PERTURBATION *STATIC *DLOAD 1, P1, 1000.0 1, P2, 1000.0 1, P3, 1000.0 1, P5, 1000.0 2, P1, 1000.0 2, P2, 1000.0 2, P3, 1000.0 2, P4, 1000.0 *DECHARGE 1, ES1, 1000.0 1, ES2, 1000.0 1, ES3, 1000.0 1, ES5, 1000.0 2, ES1, 1000.0 2, ES2, 1000.0 2, ES3, 1000.0 2, ES4, 1000.0 *CLOAD 2,1, -166.6667 3,1, 166.6667 3,2, 166.6667 3,3, 333.3333 4,1, 166.6667 4,2, -166.6667 5,2, -166.6667 5,3, -333.3333 6,1, -166.6667 6,2, 166.6667 7,1, 166.6667 7,2, 166.6667 7,3, 333.3333 8,1, 166.6667 8,2, -166.6667 9,1, 1333.3333 9,2, 666.6667 9,3, 666.6667 10,1, 666.6667 10,3, -666.6667 11,2, 666.6667 11,3, -666.6667 12,1, 666.6667 12,2, 1333.3333 12,3, 666.6667 13,2, -666.6667 13,3, 666.6667 14,1, -666.6667 14,2, -1333.3333 14,3, -666.6667 15,1, -1333.3333 15,2, -666.6667 15,3, -666.6667 16,1, -666.6667 16,3, 666.6667 17,1, 666.6667 17,2, 666.6667 17,3, 1333.3333 18,1, 666.6667 18,2, -666.6667 19,1, -666.6667 19,2, -666.6667 19,3, -1333.3333 20,1, -666.6667 20,2, 666.6667 21,1, 1333.3333 21,2, 1333.3333 22,1, -1333.3333 22,2, -1333.3333 *CECHARGE 1,,333.333 2,,333.333 3,,333.333 4,,333.333 5,,333.333 6,,333.333 7,,333.333 8,,333.333 9,,-1333.333 10,,-1333.333 11,,-1333.333 12,,-1333.333 13,,-666.667 14,,-666.667 15,,-666.667 16,,-666.667 17,,-1333.333 18,,-1333.333 19,,-1333.333 20,,-1333.333 21,,-1333.333 *EL PRINT NFORC, *SECTION PRINT,NAME=SID5678,SURFACE=SIDE,AXES=LOCAL ,1.,1,1. ,0.,1.,1., ,1.,0.,1 SOAREA,SOF SOCF, SOE, *SECTION FILE,NAME=SID5678,SURFACE=SIDE,AXES=LOCAL ,1.,1,1. ,0.,1.,1., ,1.,0.,1. SOAREA,SOF SOCF, SOE, *ENERGY PRINT *NODE PRINT U,COORD,EPOT RCHG,RF *EL PRINT S,COORD, E, LOADS, CHRGS, EFLX,EPG ENER, ELEN, *OUTPUT,FIELD *ELEMENT OUTPUT S,COORD,E, EFLX,EPG ENER, ELEN, *NODE OUTPUT U,COORD,EPOT,RCHG,RF *OUTPUT,HISTORY *ENERGY OUTPUT,VAR=ALL *END STEP *step,pert *steady state dynamics, direct 50., *DLOAD 1, P1, 2000.0 1, P2, 2000.0 1, P3, 2000.0 1, P5, 2000.0 2, P1, 2000.0 2, P2, 2000.0 2, P3, 2000.0 2, P4, 2000.0 *DECHARGE 1, ES1, 2000.0 1, ES2, 2000.0 1, ES3, 2000.0 1, ES5, 2000.0 2, ES1, 2000.0 2, ES2, 2000.0 2, ES3, 2000.0 2, ES4, 2000.0 *node print u, rf, *el print s, e, *output, field, op=new *node output u, rf, *element output s, e, *end step