求助:压电瞬态分析

jackyfoo

做压电材料的高人帮帮忙啊

我搞不懂压电的瞬态分析了, 毕业论文要用啊, 十分着急, 请高人指教!!!
压电材料加载电压发生形变的瞬态分析是怎么设置的啊???  
我的程序是这样的:不论时间time语句怎么改,得出的曲线都是一个样子.

FINISH
/CLEAR
/PREP7
CYLIND,0.004, ,0.0002,,0,360,   
CYLIND,0.005, ,-0.0002,,0,360,   
vglue,all
ET,1,SOLID98
ET,2,SOLID187
keyopt,1,1,3
MP,DENS,1,7500
MP,PERX,1,1695
MP,PERY,1,1695
MP,PERZ,1,1469
!orthotropic dielectric matrix
TB,PIEZ,1,,,0   
TBMODIF,1,1,0   
TBMODIF,1,2,0   
TBMODIF,1,3,-6.5
TBMODIF,2,1,0   
TBMODIF,2,2,0   
TBMODIF,2,3,-6.5
TBMODIF,3,1,0   
TBMODIF,3,2,0   
TBMODIF,3,3,23  
TBMODIF,4,1,0   
TBMODIF,4,2,0   
TBMODIF,4,3,0   
TBMODIF,5,1,
TBMODIF,5,2,17  
TBMODIF,5,3,0   
TBMODIF,6,1,17  
TBMODIF,6,2,0   
TBMODIF,6,3,0   
!piezoelectric stress matrix
TB,ANEL,1,1,21,0
TBTEMP,0
TBDATA,,12.6e10,7.95e10,8.41e10,0,0,0   
TBDATA,,12.6e10,8.41e10,0,0,0,11.7e10   
TBDATA,,0,0,0,2.32e10,0,0   
TBDATA,,2.3e10,0,2.3e10,,,
MP,DENS,2,8800
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,2,,11.7e10
MPDATA,PRXY,2,,0.374
esize,0.0004
MSHKEY,0
MSHAPE,1,3d
TYPE,1
MAT,1
VMESH,1
TYPE,2
MAT,2
VMESH,3

asel,s,,,3,4
asel,a,,,7,8
da,all,ux,0
da,all,uy,0
da,all,uz,0
alls
DA,1,VOLT,50
DA,2,VOLT,0
FINISH
/SOL
ANTYPE,4
SOLCONTROL,0
TRNOPT, FULL, , , , ,HHT
NLGEOM,ON
OUTRES,,ALL
OUTPR,NSOL,1
TIME,0.1
DELTIM,0.01
KBC,0                          ! RAMP LOAD OVER LOAD STEP
TIMINT,OFF,STRUC            ! SUPPRESS STRUCTURAL DYNAMICS
nsubst,10
TINTP,0, , ,0.5, , ,
/OUTPUT,pzt13-3,txt
SOLVE

/POST26
NSOL,2,3024,UZ !mid
PLVAR,2