求助压电仿真，怎么加载正弦电压？

renmi1985

我在做声表面波传感器的压电分析，需要在耦合电极上加载5v的正弦电压，还请各位大侠帮帮忙啊。。。谢谢！

damao5277

可以用场路耦合吗？ 1# renmi1985

tolly8210

你应该是做稳态分析吧? 如果是，只要加一个RMS值就可以了

renmi1985

你好，我是菜鸟，能具体讲一下吗？

renmi1985

这是我的命令流，可是叉指电极的建模都出问题了，还请指点一下啊。。。
/COM First test based on vm176 from Ansys Help
/com last changed 060612
! data from <3D FE simulation of H2 SAW gas sensor>
/FILNAME, SAW_sensor
/PREP7
SMRT,OFF
/TITLE, SAW sensor

/COM  ***
ET,1,SOLID98,3                   ! 3-D COUPLED-FIELD SOLID
ET,2,solid92
/COM  ***  MATERIAL PROPERTIES ***
MP,EX,1,1.21e11                  ! palladium modulus of elastic
MP,DENS,1,12000                  ! palladium density
MP,PRXY,1,0.39    ! poisson's ratio of palladium
MP,EX,2,7.03E10                  ! ALUMINUM MODULUS OF ELASTICITY
MP,NUXY,2,.345                   ! ALUMINUM POISSON RATIO
MP,DENS,2,2690                   ! ALUMINUM DENSITY
!***set the piezoelectric data***
density=7500    !density
reps11=730
reps33=635                      !permitivity
c11=13.9e10    !stiffness matrix
c12=0.778e11
c13=0.743e11
c33=1.15e11
c44=0.256e11
c66=0.306e11
e15=12.7
e31=-5.2
e33=15.1


MP,DENS,3,7500                   ! PZT4 DENSITY
TB,PIEZ,3                        ! DEFINE PIEZO. TABLE FOR LiNbO3 (FROM THE INTERNET)
TBDATA,3,e31
TBDATA,6,e31
TBDATA,9,e33
TBDATA,14,e15
TBDATA,16,e15
MP,PERX,3,reps11               ! Permittivity of LiNbO3  
MP,PERY,3,reps11  
MP,PERZ,3,reps33  
TB,ANEL,3,,,0                       ! DEFINE STRUCTURAL TABLE FOR LiNbO3
TBDATA,1, c11,c12,c13   
TBDATA,7,c11,c13
TBDATA,12,c33
TBDATA,16,c66  
TBDATA,19,c44   
TBDATA,21,c44        !

/COM  *** DEFINE GEOMETRIC PARAMETERS ***
ADIV=15                           ! element size
L=523.2e-6                        ! Length and width of the substrate
H=138e-6                       ! depth of the substrate
LP=138e-6                      ! length of the palladium film
WP=65e-6                          ! width of the palladium film
HP=1e-6                           ! height of the palladium film
LQ=5.75e-6                          ! quarter of lambda
LPth=138e-6                    ! length of the delay path


/COM  *** DEFINE THE MODEL ***
BLOCK,0,H,-L/2,L/2,-L/2,L/2  ! Define the piezoelectric body

!BLOCK,-HP,0,-LP/2,LP/2,-WP/2,WP/2  ! Define the palladium film


!BLC4,(-LPth/2-LQ),(-LP/2),LQ,LP  ! Define the IDT
!BLC4,(-LPth/2-3*LQ),(-LP/2),LQ,LP
!BLC4,(-LPth/2-5*LQ),(-LP/2),LQ,LP
!BLC4,(-LPth/2-7*LQ),(-LP/2),LQ,LP
!BLC4,(-LPth/2-9*LQ),(-LP/2),LQ,LP
!BLC4,(-LPth/2-11*LQ),(-LP/2),LQ,LP
!BLC4,(LPth/2),(-LP/2),LQ,LP
!BLC4,(LPth/2+2*LQ),(-LP/2),LQ,LP
!BLC4,(LPth/2+4*LQ),(-LP/2),LQ,LP
!BLC4,(LPth/2+6*LQ),(-LP/2),LQ,LP
!BLC4,(LPth/2+8*LQ),(-LP/2),LQ,LP
!BLC4,(LPth/2+10*LQ),(-LP/2),LQ,LP
VSEL,ALL
VGLUE,1,2

/COM  *** MESHING ***
VSEL,all
SHPP,LSTET,ON
VATT,3,,1   ! check
lsel,s,line,,1
lsel,a,line,,3
lsel,a,line,,6
lsel,a,line,,8,12                 ! select the lines to be divided
lesize,all,,,25                  ! divide the selected lines to 20 division
lsel,s,line,,2
lsel,a,line,,7
lsel,a,line,,4,5
lesize,all,,,6
MSHK,0                           ! free VOLUME MESH
vmesh,1  
               ! mesh
arefine,5                       ! refine the mesh on area 7 and 13
VSEL,ALL
VSEL,S,VOLU,,2
VATT,1,,2
MSHK,0
vMESH,2                          ! number of node not enough.
!asel,s,area,,7,18
!aesize,all,LQ/2
!amesh,all
!nummrg,node,LQ/4


/COM  *** Boundary condition; bearing ***
NSEL,ALL
/COM fix the area
DA,6,UX,0
DA,6,UY,0
DA,6,UZ,0
!select the nodes on the electrode to be zero
NSEL,S,LOC,x,0
NSEL,R,LOC,Y,-LP/2,LP/2
NSEL,R,LOC,z,-LPTH/2-LQ,-LPTH/2
*GET,NULLIDT,NODE,,NUM,MIN
CP,1,VOLT,ALL  
NSEL,ALL
NSEL,S,LOC,x,0
NSEL,R,LOC,z,-LPth/2-5*LQ,-LPth/2-4*LQ
NSEL,R,LOC,Y,-LP/2,LP/2
CP,high,VOLT,ALL
NSEL,ALL
NSEL,S,LOC,x,0
NSEL,R,LOC,z,-LPth/2-9*LQ,-LPth/2-8*LQ
NSEL,R,LOC,Y,-LP/2,LP/2
CP,high,VOLT,ALL
NSEL,ALL

NSEL,S,LOC,x,0
NSEL,R,LOC,z,-LPth/2-3*LQ,-LPth/2-2*LQ
NSEL,R,LOC,Y,-LP/2,LP/2
*GET,INPUTIDT,NODE,,NUM,MIN
CP,2,VOLT,ALL    !select and couple the input IDT which is not xero
NSEL,ALL
NSEL,S,LOC,x,0
NSEL,R,LOC,z,-LPth/2-7*LQ,-LPth/2-6*LQ
NSEL,R,LOC,Y,-LP/2,LP/2
CP,HIGH,VOLT,ALL
NSEL,ALL
NSEL,S,LOC,x,0
NSEL,R,LOC,z,-LPth/2-11*LQ,-LPth/2-10*LQ
NSEL,R,LOC,Y,-LP/2,LP/2
CP,HIGH,VOLT,ALL
NSEL,ALL

NSEL,S,LOC,x,0
NSEL,R,LOC,z,LPth/2,LPth/2+LQ
NSEL,R,LOC,Y,-LP/2,LP/2
*GET,OUTPUT1,NODE,,NUM,MIN   ! get the min. node mumber which is on the output IDT
CP,3,VOLT,ALL
NSEL,ALL
NSEL,S,LOC,x,0
NSEL,R,LOC,z,LPth/2+4*LQ,LPth/2+5*LQ
NSEL,R,LOC,Y,-LP/2,LP/2
CP,HIGH,VOLT,ALL
NSEL,ALL
NSEL,S,LOC,x,0
NSEL,R,LOC,z,LPth/2+8*LQ,LPth/2+9*LQ
NSEL,R,LOC,Y,-LP/2,LP/2
CP,HIGH,VOLT,ALL
NSEL,ALL
NSEL,S,LOC,x,0
NSEL,R,LOC,z,LPth/2+2*LQ,LPth/2+3*LQ
NSEL,R,LOC,Y,-LP/2,LP/2
*GET,OUTPUT2,NODE,,NUM,MIN
CP,4,VOLT,ALL
NSEL,ALL
NSEL,S,LOC,x,0
NSEL,R,LOC,z,LPth/2+6*LQ,LPth/2+7*LQ
NSEL,R,LOC,Y,-LP/2,LP/2
CP,HIGH,VOLT,ALL
NSEL,ALL
NSEL,S,LOC,x,0
NSEL,R,LOC,z,LPth/2+10*LQ,LPth/2+11*LQ
NSEL,R,LOC,Y,-LP/2,LP/2
CP,HIGH,VOLT,ALL
save
FINISH

tolly8210

叉指电极 可以不用真的建叉指，你可以等效成d31模式，只要把压电，介电等参数改一下，具体看
http://research-tlh.appspot.com/docs/journal1.pdf

renmi1985

可是到后面的划分网格还是出错，可以加QQ聊吗？524360734，真心想向您请教！谢谢了！