*Heading ** Job name: blocks-surface Model name: Blocks ** ** PARTS ** *Part, name=Blocks *End Part *Part, name="RIGID SURFACE" *End Part ** ** ASSEMBLY ** *Assembly, name=Assembly ** *Instance, name=Blocks-1, part=Blocks *Node 1, 0., 1., 0. 2, 0., 0., 0. 3, 0., 1., 1. 4, 0., 0., 1. 5, 1., 1., 0. 6, 1., 0., 0. 7, 1., 1., 1. 8, 1., 0., 1. *Element, type=C3D8E 1, 5, 6, 8, 7, 1, 2, 4, 3 *Orientation, name=Ori-4 1., 0., 0., 0., 1., 0. 3, 0. ** Region: (block:Picked), (Material Orientation:Picked) *Elset, elset=_I1, internal 1, ** Section: block *Solid Section, elset=_I1, orientation=Ori-4, material=PZT-5H 1., *End Instance ** ***Instance, name=Blocks-2, part=Blocks ***Node ** 1, 0., 1., 0. ** 2, 0., 0., 0. ** 3, 0., 0., 1. ** 4, 0., 1., 1. ** 5, 1., 0., 0. ** 6, 1., 0., 1. ** 7, 1., 1., 1. ** 8, 1., 1., 0. ** 9, 0.499939, 0.5, 0.5 ** 10, 0.24997, 0.25, 0.25 ** 11, 0.5, 0., 0. ** 12, 0.74997, 0.25, 0.25 ** 13, 1., 0.5, 0. ** 14, 0.74997, 0.75, 0.25 ** 15, 0.5, 0.5, 0. ** 16, 0.74997, 0.25, 0.75 ** 17, 1., 0., 0.5 ** 18, 0.5, 0., 0.5 ** 19, 0.74997, 0.75, 0.75 ** 20, 1., 1., 0.5 ** 21, 1., 0.5, 0.5 ** 22, 1., 0.5, 1. ** 23, 0.24997, 0.75, 0.75 ** 24, 0.5, 1., 0.5 ** 25, 0.5, 1., 1. ** 26, 0.5, 0.5, 1. ** 27, 0.5, 0., 1. ** 28, 0.24997, 0.25, 0.75 ** 29, 0., 0.5, 1. ** 30, 0., 0.5, 0.5 ** 31, 0.24997, 0.75, 0.25 ** 32, 0., 1., 0.5 ** 33, 0.5, 1., 0. ** 34, 0., 0.5, 0. ** 35, 0., 0., 0.5 ***Element, type=C3D10E ** 1, 9, 5, 2, 8, 12, 11, 10, 14, 13, 15 ** 2, 9, 5, 6, 2, 12, 17, 16, 10, 11, 18 ** 3, 9, 8, 7, 5, 14, 20, 19, 12, 13, 21 ** 4, 9, 5, 7, 6, 12, 21, 19, 16, 17, 22 ** 5, 9, 8, 4, 7, 14, 24, 23, 19, 20, 25 ** 6, 9, 4, 6, 7, 23, 26, 16, 19, 25, 22 ** 7, 9, 3, 6, 4, 28, 27, 16, 23, 29, 26 ** 8, 9, 1, 3, 4, 31, 30, 28, 23, 32, 29 ** 9, 9, 1, 4, 8, 31, 32, 23, 14, 33, 24 **10, 9, 8, 2, 1, 14, 15, 10, 31, 33, 34 **11, 9, 2, 3, 1, 10, 35, 28, 31, 34, 30 **12, 9, 2, 6, 3, 10, 18, 16, 28, 35, 27 ***Orientation, name=Ori-4 ** 1., 0., 0., 0., 1., 0. **3, 0. **** Region: (block:Picked), (Material Orientation:Picked) ***Elset, elset=_I1, internal, generate ** 1, 12, 1 **** Section: block ***Solid Section, elset=_I1, orientation=Ori-4, material=PZT-5H **1., ***End Instance ** *Instance, name=Blocks-3, part=Blocks *Node 1, 0., 1., 0. 2, 0., 0., 0. 3, 0., 1., 1. 4, 0., 0., 1. 5, 1., 1., 0. 6, 1., 0., 0. 7, 1., 1., 1. 8, 1., 0., 1. 9, 1., 1., 0.5 10, 1., 0.5, 1. 11, 1., 0., 0.5 12, 1., 0.5, 0. 13, 0., 0.5, 0. 14, 0., 0., 0.5 15, 0., 0.5, 1. 16, 0., 1., 0.5 17, 0.5, 0., 0. 18, 0.5, 1., 0. 19, 0.5, 0., 1. 20, 0.5, 1., 1. *Element, type=C3D20RE 1, 5, 6, 8, 7, 1, 2, 4, 3, 12, 11, 10, 9, 13, 14, 15, 16, 18, 17, 19, 20 *Orientation, name=Ori-4 1., 0., 0., 0., 1., 0. 3, 0. ** Region: (block:Picked), (Material Orientation:Picked) *Elset, elset=_I1, internal 1, ** Section: block *Solid Section, elset=_I1, orientation=Ori-4, material=PZT-5H 1., *End Instance ** *Instance, name=Blocks-4, part=Blocks *Node 1, 0., 1., 0. 2, 0., 0., 0. 3, 0., 1., 1. 4, 0., 0., 1. 5, 1., 1., 0. 6, 1., 0., 0. 7, 1., 1., 1. 8, 1., 0., 1. *Element, type=C3D6E 1, 1, 4, 2, 5, 8, 6 2, 1, 3, 4, 5, 7, 8 *Orientation, name=Ori-4 1., 0., 0., 0., 1., 0. 3, 0. ** Region: (block:Picked), (Material Orientation:Picked) *Elset, elset=_I1, internal 1,2 ** Section: block *Solid Section, elset=_I1, orientation=Ori-4, material=PZT-5H 1., *End Instance ** *Instance, name=Blocks-5, part=Blocks *Node 1, 0., 1., 0. 2, 0., 0., 0. 3, 0., 0., 1. 4, 0., 1., 1. 5, 1., 0., 0. 6, 1., 0., 1. 7, 1., 1., 1. 8, 1., 1., 0. 9, 0.499939, 0.5, 0.5 *Element, type=C3D4E 1, 9, 5, 2, 8 2, 9, 5, 6, 2 3, 9, 8, 7, 5 4, 9, 5, 7, 6 5, 9, 8, 4, 7 6, 9, 4, 6, 7 7, 9, 3, 6, 4 8, 9, 1, 3, 4 9, 9, 1, 4, 8 10, 9, 8, 2, 1 11, 9, 2, 3, 1 12, 9, 2, 6, 3 *Orientation, name=Ori-4 1., 0., 0., 0., 1., 0. 3, 0. ** Region: (block:Picked), (Material Orientation:Picked) *Elset, elset=_I1, internal, generate 1, 12, 1 ** Section: block *Solid Section, elset=_I1, orientation=Ori-4, material=PZT-5H 1., *End Instance ** *Instance, name="RIGID SURFACE-1", part="RIGID SURFACE" 1., 0., 0. 1., 0., 0., 1., 1., 0., 90. *Node 1, 0., 0., -1. *Nset, nset="RIGID SURFACE-1-REFPT_", internal 1, *Surface, type=CYLINDER, name=RigidSurface_, internal START, 0., 0. LINE, 0., 2. *Rigid Body, ref node="RIGID SURFACE-1-REFPT_", analytical surface=RigidSurface_ *End Instance ** *Instance, name="RIGID SURFACE-2", part="RIGID SURFACE" 0., 0., 1. *Node 1, 0., 0., -1. *Nset, nset="RIGID SURFACE-2-REFPT_", internal 1, *Surface, type=CYLINDER, name=RigidSurface_, internal START, 0., 2. LINE, 0., 0. *Rigid Body, ref node="RIGID SURFACE-2-REFPT_", analytical surface=RigidSurface_ *End Instance ** *Instance, name="RIGID SURFACE-3", part="RIGID SURFACE" 1., -1.11022302462516e-16, 0. 1., -1.11022302462516e-16, 0., 1.57735026918963, 0.577350269189626, 0.577350269189626, 120. *Node 1, 0., 0., -1. *Nset, nset="RIGID SURFACE-3-REFPT_", internal 1, *Surface, type=CYLINDER, name=RigidSurface_, internal START, 0., 2. LINE, 0., 0. *Rigid Body, ref node="RIGID SURFACE-3-REFPT_", analytical surface=RigidSurface_ *End Instance *Nset, nset=SUPPORT, instance="RIGID SURFACE-1" 1, *Nset, nset=SUPPORT, instance="RIGID SURFACE-2" 1, *Nset, nset=SUPPORT, instance="RIGID SURFACE-3" 1, *Nset, nset=BOTTOM, instance=Blocks-1 1, 2, 5, 6 ***Nset, nset=BOTTOM, instance=Blocks-2 ** 1, 2, 5, 8, 11, 13, 15, 33, 34 *Nset, nset=BOTTOM, instance=Blocks-3 1, 2, 5, 6, 12, 13, 17, 18 *Nset, nset=BOTTOM, instance=Blocks-5 1, 2, 5, 8 *Nset, nset=BOTTOM, instance=Blocks-4 1, 2, 5, 6 *Nset, nset=TOP, instance=Blocks-1 3, 4, 7, 8 ***Nset, nset=TOP, instance=Blocks-2 ** 3, 4, 6, 7, 22, 25, 26, 27, 29 *Nset, nset=TOP, instance=Blocks-3 3, 4, 7, 8, 10, 15, 19, 20 *Nset, nset=TOP, instance=Blocks-5 3, 4, 6, 7 *Nset, nset=TOP, instance=Blocks-4 3, 4, 7, 8 *Elset, elset=__G99_S3, internal, instance=Blocks-1 1, ***Elset, elset=__G99_S3, internal, instance=Blocks-2 ** 1,10 *Elset, elset=__G99_S3, internal, instance=Blocks-3 1, *Elset, elset=__G99_S3, internal, instance=Blocks-5 1,10 *Elset, elset=__G99_S5, internal, instance=Blocks-4 1, *Surface, type=ELEMENT, name=_G99, internal __G99_S3, S3 __G99_S5, S5 *Elset, elset=__G106_S2, internal, instance=Blocks-1 1, *Elset, elset=__G106_S2, internal, instance=Blocks-3 1, *Elset, elset=__G106_S1, internal, instance=Blocks-4 1,2 ***Elset, elset=__G106_S3, internal, instance=Blocks-2 ** 8, 11 *Elset, elset=__G106_S3, internal, instance=Blocks-5 8, 11 *Surface, type=ELEMENT, name=_G106, internal __G106_S2, S2 __G106_S3, S3 __G106_S1, S1 *Elset, elset=__G112_S4, internal, instance=Blocks-1 1, *Elset, elset=__G112_S4, internal, instance=Blocks-3 1, *Elset, elset=__G112_S4, internal, instance=Blocks-4 1, ***Elset, elset=__G112_S3, internal, instance=Blocks-2 ** 2, 12 *Elset, elset=__G112_S3, internal, instance=Blocks-5 2, 12 *Surface, type=ELEMENT, name=_G112, internal __G112_S4, S4 __G112_S3, S3 *End Assembly ** ** MATERIALS ** *Material, name=PZT-5H *Density 7500., *Dielectric, type=ORTHO 1.505e-08, 1.505e-08, 1.301e-08 *Elastic, type=ENGINEERING CONSTANTS 6.061e+10, 6.061e+10, 4.831e+10, 0.289, 0.512, 0.512, 2.35e+10, 2.3e+10 2.3e+10, *Piezoelectric, type=E 0., 0., 0., 0., 7.41e-10, 0., 0., 0. 0., 0., 0., 7.41e-10, -2.74e-10, -2.74e-10, 5.93e-10, 0. 0., 0. ** ** INTERACTION PROPERTIES ** *Surface Interaction, name=IntProp-1 1., *Surface Behavior, no separation, pressure-overclosure=HARD ** ** BOUNDARY CONDITIONS ** ** Name: BC-3 Type: Electric potential *Boundary BOTTOM, 9, 9 ** Name: BC-4 Type: Displacement/Rotation *Boundary SUPPORT, 1, 1 SUPPORT, 2, 2 SUPPORT, 3, 3 SUPPORT, 4, 4 SUPPORT, 5, 5 SUPPORT, 6, 6 ** ** INTERACTIONS ** ** Interaction: Int-1 *Contact Pair, interaction=IntProp-1, small sliding, adjust=0.01 _G99, "RIGID SURFACE-1".RigidSurface_ ** Interaction: Int-2 *Contact Pair, interaction=IntProp-1, small sliding, adjust=0.01 _G106, "RIGID SURFACE-2".RigidSurface_ ** Interaction: Int-3 *Contact Pair, interaction=IntProp-1, small sliding, adjust=0.01 _G112, "RIGID SURFACE-3".RigidSurface_ ** ---------------------------------------------------------------- ** ** STEP: Step-1 ** *Step, name=Step-1, nlgeom Apply charge on TOP *Static 1., 1., 0.01, 1. ** ** BOUNDARY CONDITIONS ** *CONTROLS,PARAMETERS=FIELD,FIELD=DISPLACEMENT , , ,100. ** Name: BC-8 Type: Electric potential *Boundary TOP, 9, 9, 1. ** ** OUTPUT REQUESTS ** *Restart, write, frequency=0 ** ** FIELD OUTPUT: F-Output-1 ** *Output, field, frequency=9999 *ELEMENT OUTPUT S, EFLX *NODE OUTPUT EPOT, *El Print, freq=0 *Node Print, freq=0 *End Step ** ---------------------------------------------------------------- ** ** STEP: Step-2 ** *Step, name=Step-2, nlgeom Rotate 90 degrees about x-axis *Static 0.1, 1., 1e-05, 1. ** ** BOUNDARY CONDITIONS ** ** Name: BC-3 Type: Electric potential *Boundary, op=NEW BOTTOM, 9, 9 ** Name: BC-5 Type: Displacement/Rotation *Boundary, op=NEW SUPPORT, 1, 1 SUPPORT, 2, 2 SUPPORT, 3, 3 SUPPORT, 4, 4, 1.5708 SUPPORT, 5, 5 SUPPORT, 6, 6 ** Name: BC-8 Type: Electric potential *Boundary, op=NEW TOP, 9, 9, 1. ** ** *End Step ** ---------------------------------------------------------------- ** ** STEP: Step-3 ** *Step, name=Step-3, nlgeom Rotate 90 degrees about y-axis *Static 0.05, 1., 1e-05, 1. ** ** BOUNDARY CONDITIONS ** ** Name: BC-6 Type: Displacement/Rotation *Boundary SUPPORT, 1, 1 SUPPORT, 2, 2 SUPPORT, 3, 3 SUPPORT, 4, 4, 1.5708 SUPPORT, 5, 5, 1.5708 SUPPORT, 6, 6 ** *End Step ** ---------------------------------------------------------------- ** ** STEP: Step-4 ** *Step, name=Step-4, nlgeom Rotate 90 degrees about z-axis *Static 0.05, 1., 1e-05, 1. ** ** BOUNDARY CONDITIONS ** ** Name: BC-7 Type: Displacement/Rotation *Boundary SUPPORT, 1, 1 SUPPORT, 2, 2 SUPPORT, 3, 3 SUPPORT, 4, 4, 1.5708 SUPPORT, 5, 5, 1.5708 SUPPORT, 6, 6, 1.5708 ** *End Step