north 发表于 2012-5-9 16:57

ANSYS谐响应分析求助

使用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 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 发表于 2012-5-10 18:02

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

north 发表于 2012-5-11 10:42

问题解决了。
页: [1]
查看完整版本: ANSYS谐响应分析求助