ANSYS谐响应分析求助

north

使用ANSYS进行谐响应分析，得不出结果，不知道是什么原因。
PS:已经做过模态分析，并得出固有振动模态的频率。下面是我的命令流，希望大家都够帮我找到原因。
/PREP7


*SET,R,0.001                               !radius of the circle
*SET,A,0.0025                             !rectangle half length of side
*SET,B,0.0005                             !piezoceramics width
*SET,L,0.015                               !the stator's length


!*  
ET,1,PLANE13
!*  
ET,2,SOLID5  
!*  
!*   
MPTEMP,1,0  
MPDATA,EX,1,,1.10E11               
MPDATA,PRXY,1,,0.37  
MPTEMP,1,0  
MPDATA,DENS,1,,8900
!*  
!*   
  
MPTEMP,1,0
MPDATA,DENS,2,,7600
TB,PIEZ,2                        ! DEFINE PIEZO. TABLE FOR NEPEC
TBDATA,3,-6.10                   ! PIEZO MATRIX CONSTANTS   
TBDATA,6,-6.10  
TBDATA,9,15.70  
MP,PERX,2,8.7969E-9              ! PERMITTIVITY
TB,ANEL,2                        ! 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
!*  



K,3,A,A,,   
K,9,A,-A,,  
CSYS,1  
K,8,(R+A)/2,-45,,   
K,7,R,-45,,
K,5,(R+A)/2,0,,
K,4,R,0,,   
K,1,R,45,,  
K,2,(R+A)/2,45,,
LSTR,       1,       2  
LESIZE,1,,,1,,,,,1
LSTR,       2,       3  
LESIZE,2,,,4,,,,,1  
LSTR,       3,       9  
LESIZE,3,,,10,,,,,1
LSTR,       9,       8  
LESIZE,4,,,4,,,,,1  
LSTR,       8,       7
LESIZE,5,,,1,,,,,1  
LARC,       7,       1,       4
LESIZE,6,,,10,,,,,1
LARC,       8,       2,       5
LESIZE,7,,,10,,,,,1
A,1,2,8,7
ASEL,S,AREA,,1
AATT,       1, ,   1,       0,   
A,2,3,9,8  
ASEL,S,AREA,,2
AATT,       1, ,   1,       0,
ASEL,ALL

ASEL,ALL
AGEN,4,1,2, , ,90, , ,0


NUMMRG,ALL, , , ,LOW   
NUMCMP,KP   
NUMCMP,LINE


CSYS,0

K,15,-A-B,A,,   
K,16,-A-B,-A,,  
K,17,-A,-A-B,,  
K,18,A,-A-B,,  





LSTR,       15,      11
LESIZE,21,,,3,,,,,1
LSTR,       15,      16
LESIZE,22,,,10,,,,,1
LSTR,       16,      14
LESIZE,23,,,3,,,,,1
LSTR,       14,      17
LESIZE,24,,,3,,,,,1
LSTR,       17,      18
LESIZE,25,,,10,,,,,1
LSTR,       18,       9
LESIZE,26,,,3,,,,,1
A,15,11,14,16  
ASEL,S,AREA,,9
AATT,       2, ,   1,       0,
A,14,17,18,9
ASEL,S,AREA,,10
AATT,       2, ,   1,       0,
ASEL,ALL





ESIZE,0,0.5,
MSHAPE,0,2D
MSHKEY,1
AMESH,ALL

CSYS,0              !ADD
TYPE,   2   
EXTOPT,ESIZE,10,0,  
EXTOPT,ACLEAR,0
  
EXTOPT,ATTR,0,0,0   
MAT,1   
REAL,_Z4
ESYS,0  
  

VEXT,ALL, , ,0,0,L,,,,  
  
FINISH  

/SOL
!*  
ANTYPE,2
!*  
!*  

!*  
MODOPT,LANB,5
EQSLV,SPAR  
MXPAND,5, , ,0  
LUMPM,0
PSTRES,0
!*  
MODOPT,LANB,5,100,100000, ,OFF,,2
/STATUS,SOLU
SOLVE   


*GET,F1,MODE,1,FREQ
*GET,F2,MODE,2,FREQ
*GET,F3,MODE,3,FREQ

FINISH  

!/ANG,1,-45,ZM,1
!/REPLOT




/SOLU   
CSYS,0              !ADD

ANTYPE,3
NSEL,S,LOC,X,-A
NSEL,R,LOC,Y,-A,A
D,ALL,VOLT,0,  
NSEL,ALL
NSEL,S,LOC,X,-A-B
D,ALL,VOLT,100,
NSEL,ALL   




!NSEL,S,LOC,Y,-A
!D,ALL,VOLT,100,
!NSEL,ALL   
!NSEL,S,LOC,Y,-A-B
!D,ALL,VOLT,-100,      
!NSEL,ALL

HROUT,OFF              
HARFRQ,F2-400,F2+400   
NSUBST,10,            
KBC,0                  
ALPHAD,0,              
BETAD,0,               
DMPRAT,0.002,         
SOLVE   
FINISH  

north

查看了体积的属性，发现了问题，难道通过拉伸后，材料的属性并没有赋予体积，难怪没结果

north

问题解决了。