ansys压电分析求助

mark沙

我的模型建好以后施加了冲击荷载，但是最后电极层的电压都是0，求助各位大神，看应该如何修改
代码如下：FINI
/CLEAR
/PREP7
SMRT,OFF
/TITLE,Analysis of 1-3 cement-based piezoelectric composites

ET,1,SOLID5                      !水泥层材料参数
MP,EX,1,1.38E10
MP,DENS,1,2000
MP,PRXY,1,0.22

ET,2,SOLID5                      !电极层材料参数
MP,EX,2,10e10
MP,DENS,2,8430

ET,3,SOLID5,3                     !压电层材料参数
MP,DENS,3,7500
MP,PERX,1696
MP,PERY,1696
MP,PERZ,1467
TB,PIEZ,3
TBDATA,1,0,0,-6.5
TBDATA,4,0,0,-6.5
TBDATA,7,0,0,23.3
TBDATA,10,0,0,0
TBDATA,13,0,17.44,0
TBDATA,16,17.44,0,0
TB,ANEL,3
TBDATA,1,12.6E10,7.95E10,8.41E10,0,0,0
TBDATA,7,12.6E10,8.41E10,0,0,0
TBDATA,12,11.7E10,0,0,0
TBDATA,16,2.3E10,0,0
TBDATA,19,2.3E10,0
TBDATA,21,2.35E10


L=0.0005
W=0.0005
K,1,0,0
K,2,0,-W
K,3,0,-2*W
k,4,2*L,-2*W
K,5,L,-W
KGEN,2,1,5,1,,,0.004
KGEN,2,6,10,1,,,0.001
KGEN,2,11,15,1,,,0.01
KGEN,2,16,20,1,,,0.001
KGEN,2,21,25,1,,,0.004


L,1,6
L,6,11
L,11,16
L,16,21
L,21,26
LESIZE,1,,,16
LESIZE,2,,,4
LESIZE,3,,,40
LESIZE,4,,,4
LESIZE,5,,,16
ESIZE,,2
MSHK,1
MSHA,0,3D


V,1,2,5,6,7,10               !1号水泥层   1
V,2,3,4,5,7,8,9,10           !2号水泥层   2
V,6,7,10,11,12,15            !3号电极层   3
V,7,8,9,10,12,13,14,15       !4号电极层   4
V,11,12,15,16,17,20          !5号压电层   5
V,12,13,14,15,17,18,19,20    !6号水泥层   6
V,16,17,20,21,22,25          !7号电极层   7
V,17,18,19,20,22,23,24,25    !8号电极层   8
V,21,22,25,26,27,30          !9号水泥层   9
V,22,23,24,25,27,28,29,30    !10号水泥层  10
VGLUE,ALL



VSEL,S,VOLU,,1
VATT,1,,1
VSEL,S,VOLU,,2
VATT,1,,1
VSEL,S,VOLU,,3
VATT,2,,2
VSEL,S,VOLU,,4
VATT,2,,2
VSEL,S,VOLU,,5
VATT,3,,3
VSEL,S,VOLU,,6
VATT,1,,1
VSEL,S,VOLU,,7
VATT,2,,2
VSEL,S,VOLU,,8
VATT,2,,2
VSEL,S,VOLU,,9
VATT,1,,1
VSEL,S,VOLU,,10
VATT,1,,1
VSEL,ALL
VMESH,ALL



*SET,V1,0                        !设置压电耦合
NSEL,S,LOC,Z,0.005
NSEL,R,LOC,X,0,0.0005
NSEL,R,LOC,Y,0,-0.0005
CP,1,VOLT,ALL
*GET,N1,NODE,,NUM,MIN
D,N1,VOLT,V1

NSEL,S,LOC,Z,0.015
NSEL,R,LOC,X,0,0.0005
NSEL,R,LOC,Y,0,-0.0005
CP,2,VOLT,ALL
*GET,N2,NODE,,NUM,MIN


NSEL,S,LOC,Z,0                   !约束位移边界条件
D,ALL,ALL

DSYM,SYMM,X,0                    !设置位移对称边界条件
DSYM,SYMM,Y,0

NSEL,ALL
FINISH


/SOLU
ANTYPE,TRANS
TRNOPT,FULL,,
TINTP,,0.25,0.5,0.5

NSUBST,500,0,0  
OUTRES,ERASE
OUTRES,ALL,1
TIME,1e-4
/VIEW,1,1,1,1
/ANG,1  
/REP,FAST