VM288 Input Listing


/COM,ANSYS MEDIA REL. 2024R2 (05/10/2024) REF. VERIF. MANUAL: REL. 2024R2
/VERIFY,VM288
/PREP7
/TITLE,VM288,ELECTRICAL INPUT ADMITTANCE ON A PIEZOELECTRIC TRANDUCER
/COM,        USING ENFORCED MOTION APPROACH
/COM,   
/COM,   REFERENCE:
/COM    KAGAWA AND YAMABUCHI, FINITE ELEMENT SIMULATION OF A COMPOSITE
/COM    PIEZOELECTRIC ULTRASONIC TRANSDUCER, IEEE TRANS. SONICS AND
/COM    ULTRASONICS, VOL. SU-26, NO.2, MARCH 1979
/COM
ET,1,SOLID186,,1                 ! 3-D STRUCTURAL SOLID
MP,DENS,3,7730                   ! NEPEC DENSITY
MP,EX,2,7.03E10                  ! ALUMINUM MODULUS OF ELASTICITY
MP,NUXY,2,.345                   ! ALUMINUM POISSON RATIO
MP,DENS,2,2690                   ! ALUMINUM DENSITY
MP,EX,4,10E9                     ! ADHESIVE MODULUS OF ELASTICITY
MP,DENS,4,1700                   ! ADHESIVE DENSITY
MP,NUXY,4,.38                    ! ADHESIVE POISSON RATIO

ET,2,SOLID226,1001               ! 3-D COUPLED-FIELD SOLID
TB,PIEZ,3                        ! DEFINE PIEZO. TABLE FOR NEPEC
TBDATA,3,-6.10                   ! PIEZO MATRIX CONSTANTS
TBDATA,6,-6.10
TBDATA,9,15.70
MP,PERX,3,993.55                 ! PERMITTIVITY
TB,ANEL,3                        ! DEFINE STRUCTURAL TABLE FOR NEPEC
TBDATA,1,12.80E10,6.8E10,6.6E10  ! INPUT [C] MATRIX FOR NEPEC
TBDATA,7,12.8E10,6.6E10
TBDATA,12,11.0E10
TBDATA,16,2.1E10
TBDATA,19,2.1E10
TBDATA,21,2.1E10

/COM                             ** DEFINE GEOMETRIC PARAMETERS **
R=27.5E-3                        ! DISK RADIUS
HA=15.275E-3                     ! ALUMINUM 1/2 HEIGHT LOCATION
HN=5E-3                          ! NEPEC 1/2 HEIGHT
HB=5.275E-3                      ! ADHESIVE MATERIAL HEIGHT
RDIV=5                           ! NO. ELEMENTS ALONG RADIUS
HADV=3                           ! NO. ELEMENTS ALONG ALUMINUM HEIGHT
HNDV=2                           ! NO. ELEMENTS ALONG NEPEC HEIGHT
HBDV=1                           ! NO. ELEMENTS ALONG ADHESIVE HEIGHT
ZRO=1E-5                         ! DEFINE ZERO FOR KEYPOINT LOCATION
CSYS,1
K,1,ZRO,-5                       ! DEFINE KEYPOINTS FOR MESH, WEDGE ELEMENTS
K,2,R,-5                         ! ARE NOT ALLOWED BY MESH MODULE SO KEYPOINTS
K,3,R,-5,HN                      ! ARE DEFINED NEAR ZERO AND LATER MERGED
K,4,ZRO,-5,HN
K,5,R,-5,HB
K,6,ZRO,-5,HB
K,7,R,-5,HA
K,8,ZRO,-5,HA
KGEN,2,1,8,1,,10
L,2,3                            ! DEFINE LINE SEGMENTS
LESIZE,1,,,HNDV
L,3,5
LESIZE,2,,,HBDV
L,5,7
LESIZE,3,,,HADV
L,2,10
LESIZE,4,,,1
V,11,3,4,12,10,2,1,9             ! CREATE NEPEC VOLUME
VATT,3,,2                        ! ASSIGN MATERIAL ATTRIBUTES
V,13,5,6,14,11,3,4,12            ! CREATE ADHESIVE VOLUME
VSEL,S,VOLU,,2
VATT,4,,1                        ! ASSIGN MATERIAL ATTRIBUTES
V,15,7,8,16,13,5,6,14            ! CREATE ALUMINUM VOLUME
VSEL,S,VOLU,,3
VATT,2,,1                        ! ASSIGN MATERIAL ATTRIBUTES
VSEL,ALL
MSHK,1                           ! MAPPED VOLUME MESH 
MSHA,0,3D                        ! USING HEX  
ESIZE,,RDIV
SHPP,OFF                         ! TURN OFF SHAPE CHECKING TO ALLOW
!                                  FOR WEDGE SHAPE MESH OF SOLIDS
TYPE,1
/OUT,SCRATCH
VMESH,ALL                        ! MESH ALL VOLUMES
NUMMRG,NODE                      ! MERGE NODES TO CREATE WEDGE ELMENTS
!                                  AT AXIS
NSEL,S,LOC,Y,-5
DSYM,SYMM,Y,1                    ! SYMMETRY B.C. AT THETA=-5 DEG.
NSEL,S,LOC,Y,5
DSYM,SYMM,Y,1                    ! SYMMETRY B.C. AT THETA=5 DEG.
NSEL,S,LOC,X,0,.001
DSYM,SYMM,X,1                    ! SYMMETRY B.C. AT X=0
NSEL,S,LOC,Z,0
DSYM,SYMM,Z,1                    ! SYMMETRY B.C. AT Z=0
NSEL,S,LOC,Z,HN                  ! SELECT NODES AT TOP ELECTRODE
CP,1,VOLT,ALL                    ! COUPLE VOLT DOF ON ELECTRODE
*GET,N1,NODE,,NUM,MIN            ! GET NODE ON ELECTRODE
NSEL,S,LOC,Z,0                   ! SELECT NODES AT SYMMETRY PLANE
D,ALL,VOLT,0.0                   ! SET VOLT TO ZERO AT SYMMETRY PLANE
NSEL,ALL
FINISH

/SOLUTION
ANTYPE,MODAL                     ! PERFORM MODAL SOLVE
MODCONT,,ON                      ! TURN ON ENFORCED MOTION CALCULATION
MODOPT,LANB,3,,                  ! BLOCK-LANCZOS EIGENSOLVER
MXPAND,,,,YES,                   ! EXPAND MODES AND COMPUTE ELEMENT RESULTS
D,N1,VOLT,1                      ! ENFORCED MOTION ID TO APPLY VOLT IN MSUP HARMONIC
SOLVE
FINISH

/POST1
/OUT,
SET,LIST                         ! LIST MODAL FREQUENCIES AND ENFORCED STATIC MODE
FINISH

/OUT,SCRATCH
/SOLUTION
ANTYPE,HARMIC                    ! PERFORM HARMONIC ANALYSIS
HROPT,MSUP                       ! USE MODE-SUPERPOSITION METHOD
OUTRES,ALL,ALL                   ! STORE EVERY SUBSTEP
HARFRQ,5000,35000                ! SOLVE FOR FREQ=5KHZ AND 35KHZ
NSUBST,20                        ! SPECIFY NUMBER OF SUBSTEPS
KBC,1                            ! STEP BOUNDARY CONDITIONS  
DVAL,1,U,-0.5                    ! APPLY ENFORCED MOTION TO TOP ELECTRODE
SOLVE
HARFRQ,39000,45000               ! SOLVE FOR FREQ=39KHZ AND 45KHZ
SOLVE
HARFRQ,46000,54000               ! SOLVE FOR FREQ=46KHZ AND 54KHZ
SOLVE
FINISH

/SOLUTION
EXPASS,ON                        ! EXPANSION PASS
NUMEXP,ALL,,,                    ! EXPAND SOLUTION FOR ALL SUBSTEPS,COMPUTE ELEMENT RESULTS
SOLVE
FINISH

/POST26
RFORCE,2,N1,CHRG                 ! STORE CHARGE ON ELECTRODE
PI2=(3.14159*2.)
PROD,3,2,1,,MHOS,,,PI2           ! CALCULATE ADMITTANCE (10 DEG. SLICE)
PROD,4,3,,,MMHO,,,36000          ! CALCULATE TOTAL ADMITTANCE (MMHOS)
/OUT,
PRVAR,4                          ! PRINT ELECTRICAL ADMITTANCE VS. FREQUENCY
/OUT,SCRATCH
*GET,F1,VARI,4,RTIME,20000
*GET,F2,VARI,4,RTIME,35000
*GET,F3,VARI,4,RTIME,40000
*GET,F4,VARI,4,RTIME,45000
*GET,F5,VARI,4,RTIME,50000
*GET,F6,VARI,4,RTIME,54000
*STATUS,PARM
*DIM,LABEL,CHAR,6,2
*DIM,VALUE,,6,3
LABEL(1,1) = 'Y MMHOS ','Y MMHOS','Y MMHOS ','Y MMHOS ','Y MMHOS ','Y MMHOS '
LABEL(1,2) = '@20 kHz','@35 kHz','@40 kHz','@45 kHz','@50 kHz','@54 kHz'
*VFILL,VALUE(1,1),DATA,.41,.9,1.80,-1.1,.39,.65
*VFILL,VALUE(1,2),DATA,F1,F2,F3,F4,F5,F6
*VFILL,VALUE(1,3),DATA,ABS(F1/.41),ABS(F2/0.9),ABS(F3/1.80),ABS(F4/(-1.1)),ABS(F5/.39),ABS(F6/0.65)
/COM
/OUT,vm288,vrt
/COM, ----------------------------VM288 RESULT COMPARISON-------------------------------------
/COM,
/COM,                   |   TARGET  |   MECHANICAL APDL |   RATIO
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,'   ',F10.2,'  ',F14.2,'   ',1F15.3)
/COM,
/COM, ------------------------------------------------------------------------------------------
/OUT
*LIST,vm288,vrt
/out,scratch
FINISH