ansys压电例子VM176别解

shmon

在看ansys verifiation manual176时，我做了自己的方法，但是总是出现elements with negative and positive charges 的 warning， 然后就invalid！
我的命令如下：尽请高手指教


/COM,ANSYS MEDIA REL. 120 (02/19/2009) REF. VERIF. MANUAL: REL. 120
/VERIFY,VM176
/PREP7
SMRT,OFF
/TITLE, VM176, FREQUENCY RESPONSE OF ELECTRICAL INPUT ADMITTANCE FOR A
/COM           PIEZOELECTRIC TRANSDUCER
/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
ET,1,SOLID226,1001                    ! 3-D COUPLED-FIELD 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
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,R,-5                       ! DEFINE KEYPOINTS FOR MESH, WEDGE ELEMENTS
K,2,ZRO,-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
K,9,
K,10,0,0,HA
   
A,1,2,4,3   
A,3,4,6,5   
A,5,6,8,7
VROTAT,1,2 ,3 , , , ,9,10 ,10, ,                                         !VROTAT, NA1, NA2, NA3, NA4, NA5, NA6, PAX1, PAX2, ARC, NSEG
LESIZE,1,,,10
LESIZE,3,,,10
LESIZE,6,,,10
LESIZE,9,,,10
LESIZE,15,,,5
LESIZE,18,,,5
LESIZE,23,,,5
LESIZE,28,,,5
LESIZE,4,,,6
LESIZE,7,,,1
LESIZE,10,,,10
ALLSEL
vsel,,,,1                                                   !VSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP
vatt,3
MSHK,1                           ! 1 FOR MAPPED VOLUME MESH
MSHA,0,3D
VMESH,ALL
ALLSEL
vsel,,,,2
vatt,4
MSHK,1                           ! 1 FOR MAPPED VOLUME MESH
MSHA,0,3D
VMESH,ALL
ALLSEL
vsel,,,,3
vatt,2
MSHK,1                           ! 1 FOR MAPPED VOLUME MESH
MSHA,0,3D
VMESH,ALL
ALLSEL
!SHPP,OFF                         ! TURN OFF SHAPE CHECKING TO ALLOW                        ! 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
D,N1,VOLT,-0.5                   ! APPLY -0.5 VOLT TO TOP 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

/SOLU
EQSLV,SPARSE    ! USING SPARSE MATRIX SOLVER   
ANTYPE,HARMIC                    ! PERFORM HARMONIC ANALYSIS
OUTRES,ALL,ALL                   ! STORE EVERY SUBSTEP
HARFRQ,5000,35000                ! SOLVE FOR FREQ=20KHZ AND 35KHZ
NSUBST,2
KBC,1                            ! STEP BOUNDARY CONDITIONS  
EQSLV,ICCG                       ! ICCG SOLVER
SOLVE
HARFRQ,39000,45000               ! SOLVE FOR FREQ=42KHZ AND 45KHZ
SOLVE
HARFRQ,46000,54000               ! SOLVE FOR FREQ=50KHZ AND 54KHZ
SOLVE
FINISH