*Heading
** Job name: Transducer_coarse2 Model name: CASE-4
**
** PARTS
**
*Part, name=PIEZO
*End Part
*Part, name=base
*End Part
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*INCLUDE,INPUT=dynamictransducer_mesh.inp
**
** Constraint: TIE-1
*Tie, name=TIE-1, adjust=no, position tolerance=1e-05, no rotation
CORE-TOP, TOP-PIEZO
** Constraint: TIE-2
*Tie, name=TIE-2, adjust=no, position tolerance=1e-05, no rotation
CORE-BOT, PIEZO-BOT
**
*Equation
2
VOLT-TOP, 9, -1.
MVOLT-TOP, 9, 1.
*Equation
2
VOLT-T-CORE, 9, -1.
MVOLT-T-CORE, 9, 1.
*Equation
2
VOLT-B-CORE, 9, -1.
MVOLT-B-CORE, 9, 1.
*Equation
2
VOLT-BOT, 9, -1.
MVOLT-BOT, 9, 1.  
*End Assembly
** 
** MATERIALS
** 
*Material, name=BASE
*Density
1500.,
*Elastic
 6e+09, 0.35
*Material, name=PZT-5H-Y
*Density
7500.,
*Dielectric, type=ORTHO
 1.505e-08, 1.301e-08, 1.505e-08
*Elastic, type=ENGINEERING CONSTANTS
 6.061e+10, 4.831e+10, 6.061e+10, 0.512, 0.289, 0.408, 2.3e+10, 2.35e+10
   2.3e+10,
*Piezoelectric, type=E
0., 0., 0., 7.41e-10, 0., 0., -2.74e-10,  5.93e-10
 -2.74e-10, 0., 0., 0., 0., 0., 0., 0.
        0.,  7.41e-10
** 
** BOUNDARY CONDITIONS
** 
** Name: BC-1 Type: Displacement/Rotation
*Boundary
SPRT-UX, 1, 1
SPRT-UY, 2, 2
SPRT-UZ, 3, 3
** Name: BC-4 Type: Electric potential
*Boundary
MVOLT-T-CORE, 9, 9
MVOLT-B-CORE, 9, 9
** ----------------------------------------------------------------
** 
*Step, name=Step-1, nlgeom, inc=1000
SHAPE FABRICATION STEP 1:  DEFORMATION INDUCED VIA 1000 VOLT LOAD
*Static
0.1, 1., 1e-05, 0.1
** 
*Boundary
MVOLT-TOP, 9, 9, 1000.
MVOLT-BOT, 9, 9, 1000.
** 
*Restart, write, frequency=0
**
*Output, field,freq=999
*Node Output
U, RF, EPOT, RCHG, CECHG
*Node Output, nset=QA_TEST
U, EPOT
*Element Output
S, E, EPG, EFLX
*Element Output, elset=QA_TEST
S,
** 
*Output, history
*Node Output, nset=MVOLT-B-CORE
EPOT, RCHG
*Node Output, nset=MVOLT-BOT
EPOT, RCHG
*Node Output, nset=MVOLT-T-CORE
EPOT, RCHG
*Node Output, nset=MVOLT-TOP
EPOT, RCHG
*Node Output, nset=N-OUT1
U1, U2, U3
*Node Output, nset=N-OUT2
U3, 
*Node Output, nset=QA_TESTH
U1, U2, U3
*El Print, freq=0
*Node Print, freq=0
*End Step
**
** -------------------------------------------------------------------  
** 
*Step, name=Step-2, nlgeom, inc=1000
SHAPE FABRICATION STEP 2: FIX ENDS & REDUCE VOLTAGE TO 0.0
*Static
0.1, 1., 1e-05, 0.1
** 
*Boundary, op=new
MVOLT-T-CORE, 9, 9, 0.0
MVOLT-B-CORE, 9, 9, 0.0
MVOLT-TOP,    9, 9, 0.0
MVOLT-BOT,    9, 9, 0.0
*Boundary, op=new, Fixed
END-1, 1, 3
END-2, 1, 3
*End Step
**
** ------------------------------------------------------------------- 
** 
*Step, name=Step-3, perturbation
Fabricated Shape: Short Circuit Modal Analysis
*Frequency, eigensolver=subspace, normalization=displacement
10, , ,18,100
** 
*Output, field
*Node Output
U,
*End Step
**
** ------------------------------------------------------------------- 
** 
*Step, name=Step-4, nlgeom, inc=1000
STEP 4: FABRICATED SHAPE FOR AN OPEN CIRCUIT (REMOVE VOLTAGE BCs)
*Static
1.0, 1.0, 1e-05, 1.0
** 
*Boundary, op=new
MVOLT-T-CORE, 9, 9, 0.0
MVOLT-B-CORE, 9, 9, 0.0
*Boundary, op=new, Fixed
END-1, 1, 3
END-2, 1, 3
*End Step
**
** ------------------------------------------------------------------- 
** 
*Step, name=Step-5, perturbation
Fabricated Shape: Open Circuit Modal Analysis
*Frequency, eigensolver=subspace, normalization=displacement
10, , ,18,100
** 
*Output, field
*Node Output
U, 
*End Step
** 
** ------------------------------------------------------------------- 
**
*Step, name=Step-6, nlgeom, inc=1000
STATIC TEST OF OPERATIONAL LOAD (SET VOLTAGE TO 200.0)
*Static
0.25, 1., 1e-05, 0.25
** 
*Boundary, op=mod
MVOLT-TOP,    9, 9, 200.0
MVOLT-BOT,    9, 9, 200.0
*End Step
**
** ------------------------------------------------------------------- 
**
*Step, name=Step-7, perturbation
Static Shape (Operational Load): Short Circuit Modal Analysis
*Frequency, eigensolver=subspace, normalization=displacement
10, , ,18,100
** 
*Output, field
*Node Output
U,
*End Step
**
** ------------------------------------------------------------------- 
** 
*Step, name=Step-8, nlgeom, inc=1000
RESET TO ZERO VOLT CONDITION PRIOR TO TRANSIENT
*Static
0.25, 1., 1e-05, 0.25
**  
*Boundary, op=mod
MVOLT-TOP,    9, 9, 0.0
MVOLT-BOT,    9, 9, 0.0
*End Step
** 
**    ############## BEGIN TRANSIENT DYNAMICS ################
**
*Step, name=Step-9, nlgeom, inc=1000, AMPLITUDE=STEP
VOLTAGE PULSE: STEP CHANGE EPOT TO 200V (0.00265 SECOND DURATION) 
*DYNAMIC, HAFTOL=400.0
2.0E-5, 0.0027, 0.5E-5, 2.0E-4 
** 
*Boundary, op=mod
MVOLT-TOP,    9, 9, 200.0
MVOLT-BOT,    9, 9, 200.0
**
*Restart, write, frequency=0 
**
** 
*Output, history, frequency=3
*Node Output, nset=MVOLT-B-CORE
EPOT, RCHG
*Node Output, nset=MVOLT-BOT
EPOT, RCHG
*Node Output, nset=MVOLT-T-CORE
EPOT, RCHG
*Node Output, nset=MVOLT-TOP
EPOT, RCHG
*Node Output, nset=N-OUT1
U1, U2, U3
*Node Output, nset=N-OUT2
U3, 
*Node Output, nset=QA_TESTH
U1, U2, U3
*El Print, freq=0
*Node Print, freq=0
**
*End Step
** 
** ------------------------------------------------------------------- 
**
*Step, name=Step-10, nlgeom, inc=1000, AMPLITUDE=STEP
FREE VIBRATION (0.0165 SECOND DURATION): STEP CHANGE EPOT TO 0V
*DYNAMIC, HAFTOL=400.0
2.0E-5, 0.0165, 0.5E-5, 2.0E-4 
** 
*Boundary, op=mod
MVOLT-TOP,    9, 9, 00.0
MVOLT-BOT,    9, 9, 00.0
*End Step
** 
** ------------------------------------------------------------------- 
**
*Step, name=Step-11, nlgeom, inc=1000, AMPLITUDE=STEP
FREE VIBRATION (0.033 SECOND DURATION): OPEN CIRCUIT CONDITION
*DYNAMIC, HAFTOL=400.0
2.0E-5, 0.033, 0.5E-5, 2.0E-4 
** 
*Boundary, op=new
MVOLT-T-CORE, 9, 9, 0.0
MVOLT-B-CORE, 9, 9, 0.0
*Boundary, op=new, Fixed
END-1, 1, 3
END-2, 1, 3
*End Step