注册 登录
Simwe仿真论坛(forum.simwe.com)，CAE/CAD/CAM/，FEA/FEM/有限元分析论坛---(手机验证注册) 返回首页

hy1966的个人空间 https://home.simwe.com/?79715 [收藏] [复制] [分享] [RSS]

留言板

facelist doodle涂鸦板

您需要登录后才可以留言 登录 | 注册

321821927 2017-4-18 16:51
老师您好，我在做一个换能器的瞬态分析，一下为我的命令流，分析没有rst文件是怎么回事？


/PREP7  
ET,1,SOLID5
KEYOPT,1,1,3
KEYOPT,1,3,0
KEYOPT,1,5,0
ET,2,SOLID95
!*  
!*  
MP,DENS,1,7600  
MP,PERX,1,600   
MP,PERY,1,900   
MP,PERZ,1,900   
TB,PIEZ,1   
TBDATA,1,14
TBDATA,4,-4.1   
TBDATA,7,-4.1   
TBDATA,11,10.3  
TBDATA,18,10.3  
TB,ANEL,1   
TBDATA,1,13.7E10,8.11E10,8.11E10
TBDATA,7,14.9E10,8.11E10
TBDATA,12,14.9E10   
TBDATA,16,3.13E10   
TBDATA,19,3.4E10
TBDABA,21,3.13E10   
MP,DENS,2,7600  
MP,PERX,2,1000  
MP,PERY,2,1290  
MP,PERZ,2,1290  
TB,PIEZ,2   
TBDATA,1,-14
TBDATA,4,4.1
TBDATA,7,4.1
TBDATA,11,-10.3
TBDATA,18,-10.3
TB,ANEL,2   
TBDATA,1,13.7E10,8.11E10,8.11E10
TBDATA,7,14.9E10,8.11E10
TBDATA,12,14.9E10   
TBDATA,16,3.13E10   
TBDATA,19,3.4E10
TBDABA,21,3.13E10   
MP,EX,3,7.03E10
MP,PRXY,3,0.34  
MP,DENS,3,2700  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,4,,226000000000   
MPDATA,PRXY,4,,0.28
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,4,,7900
K,1,0,0,0,  
K,2,0.0075,0,0,
K,3,0.0075,0,0.023,
K,4,0,0,0.023,  
K,5,0.00325,0,0.023,
K,6,0.00325,0,0.025,
K,7,0.0075,0,0.025,
K,8,0.0075,0,0.027,
K,9,0.00325,0,0.027,
K,10,0.00325,0,0.029,   
K,11,0.0075,0,0.029,
K,12,0.0075,0,0.031,
K,13,0.00325,0,0.031,   
K,14,0,0,0.031,
K,15,0,0,0.039,
K,16,0.0075,0,0.039,
K,17,0.00325,0,0.039,   
K,18,0.00325,0,0.041,   
K,19,0.0075,0,0.041,
K,20,0.0075,0,0.043,
K,21,0.00325,0,0.043,   
K,22,0.00325,0,0.045,   
K,23,0.0075,0,0.045,
K,24,0.0075,0,0.047,
K,25,0.00325,0,0.047,   
K,26,0,0,0.047,
K,27,0,0,0.087,
K,28,0.0075,0,0.087,
LSTR,       1,       2  
LSTR,       2,       3  
LSTR,       3,       5  
LSTR,       5,       4  
LSTR,       4,       1  
LSTR,       1,       2  
LSTR,       3,       7  
LSTR,       7,       6  
LSTR,       6,       5  
LSTR,       6,       9  
LSTR,       9,       8  
LSTR,       8,       7  
LSTR,       8,      11  
LSTR,      11,      10  
LSTR,      10,       9  
LSTR,      11,      12  
LSTR,      12,      13  
LSTR,      13,      14  
LSTR,      10,      13  
LSTR,      14,      15  
LSTR,      15,      17  
LSTR,      17,      16  
LSTR,      16,      12  
LSTR,      16,      19  
LSTR,      19,      18  
LSTR,      18,      17  
LSTR,      18,      21  
LSTR,      21,      20  
LSTR,      20,      19  
LSTR,      20,      23  
LSTR,      23,      22  
LSTR,      22,      21  
LSTR,      22,      25  
LSTR,      25,      24  
LSTR,      24,      23  
LSTR,      25,      26  
LSTR,      26,      27  
LSTR,      27,      28  
LSTR,      28,      24  
LLIST,ALL, , ,  
LPLOT   
FLST,2,5,4  
FITEM,2,5   
FITEM,2,2   
FITEM,2,1   
FITEM,2,4   
FITEM,2,3   
AL,P51X
FLST,2,4,4  
FITEM,2,8   
FITEM,2,3   
FITEM,2,7   
FITEM,2,6   
AL,P51X
FLST,2,4,4  
FITEM,2,9   
FITEM,2,7   
FITEM,2,10  
FITEM,2,11  
AL,P51X
FLST,2,4,4  
FITEM,2,10  
FITEM,2,13  
FITEM,2,14  
FITEM,2,12  
AL,P51X
FLST,2,4,4  
FITEM,2,18  
FITEM,2,16  
FITEM,2,15  
FITEM,2,13  
AL,P51X
FLST,2,6,4  
FITEM,2,17  
FITEM,2,19  
FITEM,2,20  
FITEM,2,21  
FITEM,2,22  
FITEM,2,16  
AL,P51X
FLST,2,4,4  
FITEM,2,24  
FITEM,2,21  
FITEM,2,25  
FITEM,2,23  
AL,P51X
FLST,2,3,4  
FITEM,2,27  
FITEM,2,28  
FITEM,2,26  
AL,P51X
FLST,2,4,4  
FITEM,2,24  
FITEM,2,26  
FITEM,2,27  
FITEM,2,28  
AL,P51X
FLST,2,4,4  
FITEM,2,31  
FITEM,2,27  
FITEM,2,30  
FITEM,2,29  
AL,P51X
FLST,2,4,4  
FITEM,2,32  
FITEM,2,33  
FITEM,2,34  
FITEM,2,30  
AL,P51X
FLST,2,5,4  
FITEM,2,36  
FITEM,2,35  
FITEM,2,33  
FITEM,2,38  
FITEM,2,37  
AL,P51X
FLST,2,1,5,ORDE,1   
FITEM,2,1   
FLST,8,2,3  
FITEM,8,1   
FITEM,8,4   
VROTAT,P51X, , , , , ,P51X, ,90, ,  
LPLOT   
FLST,2,4,5,ORDE,2   
FITEM,2,2   
FITEM,2,-5  
FLST,8,2,3  
FITEM,8,4   
FITEM,8,14  
VROTAT,P51X, , , , , ,P51X, ,90, ,  
FLST,2,2,5,ORDE,2   
FITEM,2,6   
FITEM,2,11  
FLST,8,2,3  
FITEM,8,15  
FITEM,8,27  
VROTAT,P51X, , , , , ,P51X, ,90, ,  
FLST,2,4,5,ORDE,2   
FITEM,2,7   
FITEM,2,-10
FLST,8,2,3  
FITEM,8,15  
FITEM,8,26  
VROTAT,P51X, , , , , ,P51X, ,90, ,  
APLOT   
vglue,all   
FLST,5,4,6,ORDE,4   
FITEM,5,4   
FITEM,5,9   
FITEM,5,12  
FITEM,5,14  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       1, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,4,6,ORDE,4   
FITEM,5,3   
FITEM,5,5   
FITEM,5,10  
FITEM,5,13  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       2, ,   1,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,3,6,ORDE,3   
FITEM,5,1   
FITEM,5,7   
FITEM,5,15  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       3, ,   2,       0   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
EXTOPT,VSWE,AUTO,0  
EXTOPT,ESIZE,10,0   
LPLOT   
FLST,5,1,4,ORDE,1   
FITEM,5,36  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,30, , , , ,1
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,37  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,5, , , , ,1  
!*  
FLST,5,12,4,ORDE,12
FITEM,5,3   
FITEM,5,7   
FITEM,5,10  
FITEM,5,13  
FITEM,5,16  
FITEM,5,-17
FITEM,5,20  
FITEM,5,-21
FITEM,5,24  
FITEM,5,27  
FITEM,5,30  
FITEM,5,33  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,3, , , , ,1  
!*  
FLST,5,8,4,ORDE,8   
FITEM,5,8   
FITEM,5,14  
FITEM,5,18  
FITEM,5,25  
FITEM,5,-26
FITEM,5,31  
FITEM,5,-32
FITEM,5,64  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,2, , , , ,1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,19  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,8, , , , ,1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,5   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,20, , , , ,1
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,1   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,6, , , , ,1  
!*  
!*  
VSWEEP,7,11,44  
!*  
VPLOT   
!*  
VSWEEP,14,8,70  
!*  
VSWEEP,14,10,70
!*  
VPLOT   
!*  
VSWEEP,10,9,57  
!*  
VPLOT   
!*  
VSWEEP,9,8,53   
!*  
VPLOT   
!*  
VSWEEP,13,7,67  
!*  
VPLOT   
!*  
VSWEEP,15,6,73  
!*  
VPLOT   
!*  
VSWEEP,5,5,33   
!*  
VPLOT   
!*  
VSWEEP,4,4,29   
!*  
VPLOT   
!*  
VSWEEP,3,3,25   
!*  
VPLOT   
!*  
VSWEEP,12,2,64  
!*  
APLOT   
VSWEEP,1,1,16   
!*  
APLOT
NSEL,ALL
NSEL,S,LOC,Z,-0.002  
NSEL,a,LOC,Z,-0.006  
CP,1,VOLT,ALL   
*GET,N1,NODE,,NUM,MIN   
NSEL,ALL
ALLSEL,ALL  
NSEL,ALL
NSEL,S,LOC,Z,0  
NSEL,a,LOC,Z,-0.004  
NSEL,a,LOC,Z,-0.008  
CP,2,VOLT,ALL   
*GET,N2,NODE,,NUM,MIN   
NSEL,ALL
FINISH  
/SOLU                  !进入求解器
ALLSEL,ALL
ANTYPE,TRANS           !指定分析类型
TRNOPT,FULL            !指定瞬态分析方法
!NLGEOM,OFF            !-----关闭大变形选项
TINTP,,0.25,0.5,0.5    !压电瞬态分析中，瞬态积分参数
EQSLV,SPAR
OUTPR,ALL,ALL          !控制结果的输出
OUTRES,ALL,ALL         !保存所有的数据到结果文件中
TIMINT,ON,ALL          !使用瞬态效应
MP,DAMP,1,0.0002*2/(2*3.1416*8000) !假想钢的阻尼
MP,DAMP,2,0.015*2/(2*3.1416*8000)  !假想PZT的阻尼
*do,i,1,100,1
t=i/500000
time,t
nsubst,50
v=220*sin(2*3.14159*50046*t)
D,N2,VOLT,0
D,N1,VOLT,1
allsel,all
solve
*enddo
finish
科大声学子 2016-8-3 13:52
老师，您好，有个问题向您请教一下：钛合金支架中间是谐振腔（空腔），支架上下面分别固定一片压电陶瓷片，分析其谐响应，看其导纳曲线，问题是，怎么给上下两片陶瓷片加电压，怎么将由空腔隔开的两个面耦合成等势面，即将两片陶瓷片串联起来，真心求教老师您，打扰您了
1988guo 2014-5-15 21:45
您好，有个问题想请教您一下。一个六面体结构，假设固有频率是1，那么它完全处于空气时自振频率大概也是1,。如果只有两个面接触大气，这种情况求解出的共振频率是不是也是1呢？
nokla119 2014-2-12 11:07
hy1966老师你好，有个生表面波諧振器命令流一直做不出结果，是否能帮我看一下哪里出错了吗?
/prep7  

mp,DENS,1,7700               
mp,perx,1,1704.4                              
mp,pery,1,1704.4
mp,perz,1,1433.6   
tb,ANEL,1                     
tbdata,1,1.27205E11,8.02122E10,8.46702E010  
tbdata,7,1.27205E11,8.46702E10      
tbdata,12,1.17436E11               
tbdata,16,2.29885E10               
tbdata,19,2.29885E010              
tbdata,21,2.34742E10               
tb,PIEZ,1                        
tbdata,3,-6.62281                  
tbdata,6,-6.62281                  
tbdata,9,23.2403                  
tbdata,11,17.0345                  
tbdata,13,17.0345                  

mp,dens,2,2700
EX,2,70e9
nuxy,2,0.35  


L=40e-6
W=4e-6      
LI=1e-7
WI=1e-6     
blc4,0,-40e-6,W,L,
blc4,0.5e-6,0,WI,LI,
blc4,2.5e-6,0,WI,LI,
aglue,all

et,1,plane223,1001,,2,               
et,2,plane13,7

mat,1 $ type,1
mshkey,1,   
amap,4,1,2,3,4
mat,2 $ type,2  
mshkey,0,
amesh,3  
mat,2 $ type,2
mshkey,0,
amesh,2

Cp,1,UX,40,98,
Cp,2,UX,63,102
Cp,3,UX,62,103
Cp,4,UX,61,104
Cp,5,UX,60,105
Cp,6,UX,59,106
Cp,7,UX,58,107
Cp,8,UX,57,108
Cp,9,UX,56,109
Cp,10,UX,55,110
Cp,11,UX,54,111
Cp,12,UX,53,112
Cp,13,UX,52,113
Cp,14,UX,51,114
Cp,15,UX,50,115
Cp,16,UX,49,116
Cp,17,UX,48,117
Cp,18,UX,47,118
Cp,19,UX,46,119
Cp,20,UX,45,120
Cp,21,UX,44,121
Cp,22,UX,43,122
Cp,23,UX,42,123
Cp,24,UX,41,124
Cp,25,UX,1,2

Cp,1,UY,40,98,
Cp,2,UY,63,102
Cp,3,UY,62,103
Cp,4,UY,61,104
Cp,5,UY,60,105
Cp,6,UY,59,106
Cp,7,UY,58,107
Cp,8,UY,57,108
Cp,9,UY,56,109
Cp,10,UY,55,110
Cp,11,UY,54,111
Cp,12,UY,53,112
Cp,13,UY,52,113
Cp,14,UY,51,114
Cp,15,UY,50,115
Cp,16,UY,49,116
Cp,17,UY,48,117
Cp,18,UY,47,118
Cp,19,UY,46,119
Cp,20,UY,45,120
Cp,21,UY,44,121
Cp,22,UY,43,122
Cp,23,UY,42,123
Cp,24,UY,41,124
CP,25,UY,1,2

Cp,1,UZ,40,98
Cp,2,UZ,63,102
Cp,3,UZ,62,103
Cp,4,UZ,61,104
Cp,5,UZ,60,105
Cp,6,UZ,59,106
Cp,7,UZ,58,107
Cp,8,UZ,57,108
Cp,9,UZ,56,109
Cp,10,UZ,55,110
Cp,11,UZ,54,111
Cp,12,UZ,53,112
Cp,13,UZ,52,113
Cp,14,UZ,51,114
Cp,15,UZ,50,115
Cp,16,UZ,49,116
Cp,17,UZ,48,117
Cp,18,UZ,47,118
Cp,19,UZ,46,119
Cp,20,UZ,45,120
Cp,21,UZ,44,121
Cp,22,UZ,43,122
Cp,23,UZ,42,123
Cp,24,UZ,41,124
CP,25,UZ,1,2

Cp,1,VOLT,40,98,
Cp,2,VOLT,63,102
Cp,3,VOLT,62,103
Cp,4,VOLT,61,104
Cp,5,VOLT,60,105
Cp,6,VOLT,59,106
Cp,7,VOLT,58,107
Cp,8,VOLT,57,108
Cp,9,VOLT,56,109
Cp,10,VOLT,55,110
Cp,11,VOLT,54,111
Cp,12,VOLT,53,112
Cp,13,VOLT,52,113
Cp,14,VOLT,51,114
Cp,15,VOLT,50,115
Cp,16,VOLT,49,116
Cp,17,VOLT,48,117
Cp,18,VOLT,47,118
Cp,19,VOLT,46,119
Cp,20,VOLT,45,120
Cp,21,VOLT,44,121
Cp,22,VOLT,43,122
Cp,23,VOLT,42,123
Cp,24,VOLT,41,124
CP,25,VOLT,1,2
   
ALLSEL,ALL
FINISH

/solu
antype,modal
modopt,lanb,10
mxpand,10
DL,1,,all,0
Asel,all
solve
finish
mk181 2013-7-25 16:27
hy1966老师你好，我有个命令流一直做不出结果，你能帮我看一下哪里出错了吗
/clear
/prep7  
*SET,tcem,0.25E-3  
*SET,tela,0.05E-3
*SET,tpzt,0.40E-3
*SET,r,16E-3  
ET,1,PLANE223,1001   
ET,2,PLANE82,,0

CSYS,1
MP,DENS,3,7500                     ! DENSITY
MP,PERX,3,804.6
MP,PERy,3,804.6                   ! PERMITTIVITY (X AND Y DIRECTION)
MP,PERZ,3,659.7                    ! PERMITTIVITY (Z DIRECTION)
TB,PIEZ,3                          ! DEFINE PIEZ. TABLE
TBDATA,16,10.5                     ! E61 PIEZOELECTRIC CONSTANT
TBDATA,14,10.5                     ! E52 PIEZOELECTRIC CONSTANT
TBDATA,3,-4.1                      ! E13 PIEZOELECTRIC CONSTANT
TBDATA,6,-4.1                      ! E23 PIEZOELECTRIC CONSTANT
TBDATA,9,14.1                      ! E33 PIEZOELECTRIC CONSTANT
TB,ANEL,3                          ! DEFINE STRUCTURAL TABLE
TBDATA,1,13.2E10,7.1E10,7.3E10     ! INPUT [C] MATRIX
TBDATA,7,13.2E10,7.3E10
TBDATA,12,11.5E10
TBDATA,16,3.0E10
TBDATA,19,2.6E10
TBDATA,21,2.6E10

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,4,,1.29e10
MPDATA,PRXY,4,,0.23

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,5,,3.02E10
MPDATA,PRXY,5,,0.35



PCIRC,r,r+tcem,30,90
PCIRC,r+tcem,r+tcem+tela,30,90
PCIRC,r+tcem+tela,r+tcem+tela+tpzt,30,90   
PCIRC,r+tcem+tela+tpzt,r+tcem+2*tela+tpzt,30,90
PCIRC,r+tcem+2*tela+tpzt,r+2*tcem+2*tela+tpzt,30,90

ASEL,S,AREA,,1
AATT,4,,2
ASEL,S,AREA,,5
AATT,4,,2
ASEL,S,AREA,,2
AATT,5,,2
ASEL,S,AREA,,4,
AATT,5,,2
ASEL,S,AREA,,3
AATT,3,,1

aglue,1,2,3,4,5

ASEL,ALL
lsla,s
lesize,all,,,60
lsel,r,loc,y,30
lsel,a,loc,y,90
lesize,all,,,4,,1
mshape,0,2d  
mshkey,1
SHPP,OFF
amesh,all            



CSYS,1  
NSEL,S,LOC, Y,90
D,ALL,UX,0  
D,ALL,UY,0  

LSEL,S,LINE,,11  
NSLL,S,1
CP,1,VOLT,ALL                       
*GET,BOTPZT,NODE,0,NUM,MIN
D,BOTPZT,VOLT,0

LSEL,S,LINE,,9
NSLL,S,1
CP,2,VOLT,ALL                       
*GET,TOPPZT,NODE,0,NUM,MIN       
D,TOPPZT,VOLT,100


allsel  
finish  
/SOLU   
ANTYPE,STATIC   
NSEL,ALL
SOLVE   
FINISH
乔氏佳人 2013-7-4 11:36
求助，请高人指点！由于要求取不同R（压电陶瓷半径）下的多组结果，需要循环建模。在单一R下的命令如下，能帮忙加几条循环语句吗？谢谢
!!!PZT-4、HARD AL、*PZT-4
FINISH
/CLEAR
!----------MATERIAL AND ELEMENT----------
/PREP7
ET,1,PLANE13,7,,1  !压电单元
ET,2,PLANE42,,,1   !结构单元
ET,3,PLANE13,7,,1  !压电单元
MP,DENS,1,7500     !PZT DENSITY
MP,DENS,2,2790     !HARD AL DENSITY
MP,DENS,3,7500     !*PZT DENSITY
MP,PERX,1,730      !PZT PERMITTIVITY(X DIRECTION)
MP,PERY,1,635      !PZT PERMITTIVITY(Y DIRECTION)
MP,PERZ,1,730      !PZT PERMITTIVITY(Y DIRECTION)
MP,PERX,3,730      !PZT PERMITTIVITY(X DIRECTION)
MP,PERY,3,635      !PZT PERMITTIVITY(Y DIRECTION)
MP,PERZ,3,730      !PZT PERMITTIVITY(Y DIRECTION)
MP,EX,2,7.15E10    !HARD AL YOUNG'S MODULUS
MP,PRXY,2,0.34     !HARD AL POISSON RATION
TB,PIEZ,1  !DEFINE PIZ. TABLE(压电常数)
TBDATA,2,-5.2
TBDATA,5,15.1
TBDATA,8,-5.2
TBDATA,10,12.7
TBDATA,15,12.7
TB,PIEZ,3  !DEFINE PIZ. TABLE(*压电常数)
TBDATA,2,5.2
TBDATA,5,-15.1
TBDATA,8,5.2
TBDATA,10,-12.7
TBDATA,15,-12.7
TB,ANEL,1  !DEFINE STRUCTURE TABLE(弹性常数)
TBDATA,1,13.9E10,7.43E10,7.78E10
TBDATA,7,11.5E10,7.43E10
TBDATA,12,13.9E10
TBDATA,16,2.56E10
TBDATA,19,2.56E10
TBDATA,21,3.06E10
TB,ANEL,3  !DEFINE STRUCTURE TABLE(*弹性常数)
TBDATA,1,13.9E10,7.43E10,7.78E10
TBDATA,7,11.5E10,7.43E10
TBDATA,12,13.9E10
TBDATA,16,2.56E10
TBDATA,19,2.56E10
TBDATA,21,3.06E10
!----------PARAMETER AND Geometrical model----------
R=0.005   !压电陶瓷圆片的半径
A=R+0.0006   !金属片的半径
H=0.0012  !压电振子的厚度
M=0.0004  !金属片的厚度
FREQ_MIN=20000
FREQ_MAX=40000
SUBSTEP=20
*DIM,DIS,,SUBSTEP
*DIM,FREQ,,SUBSTEP
INCR=(FREQ_MAX-FREQ_MIN)/SUBSTEP
/PREP7
K,1,0,0,0
K,2,R,0,0
K,3,R,-(H-M)/2,0
K,4,0,-(H-M)/2,0
K,5,A,-(H-M)/2,0
K,6,A,-(H+M)/2,0
K,7,0,-(H+M)/2,0
K,8,R,-(H+M)/2,0
K,9,R,-H,0
K,10,0,-H,0
A,1,2,3,4
A,4,5,6,7
A,7,8,9,10
ALLSEL,ALL
AGLUE,ALL
!----------MESH----------
ALLSEL,ALL
LESIZE,ALL,0.2E-3,,,,,,,0
ALLSEL,ALL
TYPE,1
MAT,1
MSHAPE,0,3D
MSHKEY,1
AMESH,1
ALLSEL,ALL
TYPE,2
MAT,2
MSHAPE,0,3D
MSHKEY,2
AMESH,4
ALLSEL,ALL
TYPE,3
MAT,3
MSHAPE,0,3D
MSHKEY,1
AMESH,3
!----------LOAD AND CONSTRAINT----------
!!APPLY STRUCTURAL BC
LSEL,S,LINE,,6
NSLL,S,1
D,ALL,UX,0
D,ALL,UY,0
!!(短路)APPLY ELECTRIC BC
ALLSEL,ALL
NSEL,S,LOC,Z,0
NSEL,A,LOC,Z,-(H-M)/2
NSEL,A,LOC,Z,-(H+M)/2
NSEL,A,LOC,Z,-H
CM,Z1,NODE
CP,1,VOLT,ALL
*GET,NBOT,NODE,0,NUM,MIN
!----------MODAL ANYSYS----------
/SOLU
ANTYPE,MODAL
OUTERS,ALL,ALL
MSAVE,0
MODOPT,LANB,1,FREQ_MIN,FREQ_MAX,ON,OFF
SOLVE
乔氏佳人 2013-7-4 11:34
hy1966: 在划分网格前先粘结实体，这样界面上的单元会共节点，不要指望先单个实体划分好，用NUMMRG命令
谢谢，再问一个问题啊！
乔氏佳人 2013-7-1 22:33
命令是这样的：!----------Material and element----------
/PREP7
ET,1,SOLID5,3
ET,2,SOLID45
ET,3,SOLID5,3
MP,DMPR,1,0.005  !阻尼系数
MP,DMPR,3,0.005  !阻尼系数
MP,DENS,1,7500 !PZT DENSITY
MP,DENS,2,2790 !HARD AL DENSITY
MP,DENS,3,7500 !*PZT DENSITY
MP,PERX,1,730 !PZT PERMITTIVITY(X AND Y DIRECTION)
MP,PERZ,1,635 !PZT PERMITTIVITY(Z DIRECTION)
MP,PERX,3,730 !*PZT PERMITTIVITY(X AND Y DIRECTION)
MP,PERZ,3,635 !*PZT PERMITTIVITY(Z DIRECTION)
MP,EX,2,7.15E10 !HARD AL YOUNG'S MODULUS
MP,PRXY,2,0.34  !HARD AL POISSON RATION
TB,PIEZ,1  !DEFINE PIZ. TABLE(压电常数)
TBDATA,3,-5.2  !E13 PIEZOELECTRIC CONSTANT
TBDATA,6,-5.2  !E23 PIEZOELECTRIC CONSTANT
TBDATA,9,15.1  !E33 PIEAOELECTRIC CONSTANT
TBDATA,14,12.7 ! E52 PIEZOELECTRIC CONSTANT
TBDATA,16,12.7 ! E61 PIEZOELECTRIC CONSTANT
TB,PIEZ,3  !DEFINE *PIZ. TABLE(*压电常数)
TBDATA,3,5.2  !E13 *PIEZOELECTRIC CONSTANT
TBDATA,6,5.2  !E23 *PIEZOELECTRIC CONSTANT
TBDATA,9,-15.1  !E33 *PIEAOELECTRIC CONSTANT
TBDATA,14,-12.7 !E52 *PIEZOELECTRIC CONSTANT
TBDATA,16,-12.7 !E61 *PIEZOELECTRIC CONSTANT
TB,ANEL,1  !DEFINE STRUCTURE TABLE(弹性常数)
TBDATA,1,13.9E10,7.78E10,7.43E10
TBDATA,7,13.9E10,7.43E10
TBDATA,12,11.5E10
TBDATA,16,2.56E10
TBDATA,19,3.06E10
TBDATA,21,3.06E10
TB,ANEL,3  !DEFINE STRUCTURE TABLE(*弹性常数)
TBDATA,1,13.9E10,7.78E10,7.43E10
TBDATA,7,13.9E10,7.43E10
TBDATA,12,11.5E10
TBDATA,16,2.56E10
TBDATA,19,3.06E10
TBDATA,21,3.06E10
!----------PARAMETER AND Geometrical model----------
R=0.005 !压电陶瓷圆片的半径
A=0.006 !金属片的半径
H=0.0012!压电振子的厚度
M=0.0004!金属片的厚度
FREQ_MIN=10000
FREQ_MAX=200000
SUBSTEP=190
*DIM,DIS,,SUBSTEP
*DIM,FREQ,,SUBSTEP
INCR=(FREQ_MAX-FREQ_MIN)/SUBSTEP
/PREP7
K,1,0,0,0
K,2,0,0,(H-M)/2
K,3,0,0,(H+M)/2
CIRCLE,1,R
AL,1,2,3,4
VOFFST,1,(H-M)/2
CIRCLE,2,A
AL,13,14,15,16
VOFFST,7,M
CIRCLE,3,R
AL,25,26,27,28
VOFFST,13,(H-M)/2
!----------MESH----------
LESIZE,ALL,0.2E-3,,,,,,,0
ALLSEL,ALL
TYPE,1
MAT,1
MSHAPE,0,3D
MSHKEY,1
VMESH,1
ALLSEL,ALL
TYPE,2
MAT,2
MSHAPE,0,3D
MSHKEY,1
VMESH,2
ALLSEL,ALL
TYPE,3
MAT,3
MSHAPE,0,3D
MSHKEY,1
VMESH,3
ALLSEL,ALL
NUMMRG,ALL,1E-10,,,LOW
NUMCMP,ALL
!----------LOAD AND CONSTRAINT----------
!!APPLY STRUCTURAL BC
ASEL,S,AREA,,9,12,1,1
NSLA,S,1
D,ALL,UX,0
D,ALL,UY,0
ALLSEL,ALL
ASEL,S,AREA,,7,8,1,1
NSLA,S,1
CSYS,1  !坐标系转换
NSEL,R,LOC,X,5.4E-3,5.6E-3
D,ALL,UZ,0
ALLSEL,ALL
VGLUE,ALL
!!APPLY ELECTRIC BC
ALLSEL,ALL
NSEL,S,LOC,Z,0
CM,Z1,NODE
CP,1,VOLT,ALL
*GET,NBOT,NODE,0,NUM,MIN
ALLSEL,ALL
NSEL,S,LOC,Z,H
CM,Z2,NODE
CP,2,VOLT,ALL
*GET,NTOP,NODE,O,NUM,MIN
NSEL,ALL
D,NBOT,VOLT,0
D,NTOP,VOLT,200
ALLSEL,ALL
!----------SOLUTION----------
/SOLU
ANTYPE,HARMIC
HROPT,FULL
EQSLV,SPARCE
OUTERS,ALL,ALL
HARFRQ,FREQ_MIN,FREQ_MAX
NSUBST,SUBSTEP
KBC,1
SOLVE
FINISH
!----------Z CURVE----------
/POST26
RFORCE,2,NTOP,AMPS   
CFACT,,,,1  
PROD,3,2,1,,MHOS,,,-2*3.1416/200   !计算导纳Y=-jwQ/V
REALVAR,4,3,,,ZR                   !取实部
PLVAR,4
谢谢啊！不胜感激
乔氏佳人 2013-7-1 22:26
版主好，小的在此给你请安了！真情求教，我做的实物是一个三叠片的压电振子，中间是金属圆片（硬铝），上下为极化方向相反的压电陶瓷圆片（PZT-4），激励时在压电振子的上下两面加电压，想用ANSYS仿真一下，求其导纳曲线？学了一段时间ANSYS命令，也做了一下，没做出来，求指教。命令流给你发过来，帮我看看啊
wsxadnh 2013-3-9 20:01
hy1966老师您好！
       我是在上个月曾给您去信的求教关于换能器的chigw，谢谢你给我的回复，我按照您的指点做了改正，但收效不大
我的采用的是自由网格划分(图2),无论是网格划分的疏密,结果均相差不大,我试着将水域部分改为映射网格（六面体）划分，但却划分失败,总是显示Volume 6 has invalid topology for mapped brick meshing，半径我取了0.5m但收效也不大。
     问题还可能边界条件的设置和对称约束的命令流(我瞎猜的)
不胜感激！盼望您的指正！
whhitsyz 2012-4-2 00:39
版主您好，请教一个简单的问题，对于一个平板施加一个超声振动，想模拟该振动在平板里的声场分布，是应该利用FLUID30单元进行，还是仅仅靠谐响应分析来进行？ qq543190924
zhanfei 2011-11-13 12:40
你好，能不能帮个忙，学习这个软件很久了，自己仿真出来的结果总是跟理论值相差太远。我的模型是：一个U型磁芯，上面绕上线圈。U型磁铁上部一个悬臂梁，我主要是观察悬臂梁在磁力的吸引下能偏转多少。用3D模型，不知你能否帮我写一下命令流？现在毕业设计上要用，但是一直没找到磁场跟结构厂耦合的例子，希望你能够帮帮忙。我邮箱qjjggg@163.com
.

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-26 21:45 , Processed in 0.016824 second(s), 7 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

返回顶部