charlie_x 发表于 2007-5-16 21:30:36

原创:Cymbal压电振子瞬态分析

!进入前处理
finish
/PREP7
!建立cymbal模型
*SET,t1,0.0010 !      !定义模型几何参数
*SET,t2,0.0003 !
*SET,t3,0.0003 !
*SET,r,0.015!
*SET,r0,0.008!
*SET,b,0.001!
*SET,H,0.002
*SET,r10,0.0050!内径
K,1,,r10,t1
K,2,,r10,,
K,3,,r,0
K,4,,r,0+t1
K,5,,r,t3+t1
K,6,,r-b,t3+t1
K,7,,r-b,t1
K,8,,r0,H+t2+t1
K,9,,r0,H+t1
K,10,,0,H+t2+t1
K,11,,0,H+t1
LSTR,       1,       2
LSTR,       2,       3
LSTR,       3,       4
LSTR,       4,       1
LSTR,      4,      7
LSTR,      7,      9
LSTR,      9,       11
LSTR,       11,       10
LSTR,       10,      8
LSTR,      8,      6
LSTR,      6,      5
LSTR,      5,      4
FLST,2,4,4
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,4   
AL,P51X
FLST,2,8,4
FITEM,2,5   
FITEM,2,6   
FITEM,2,7   
FITEM,2,8   
FITEM,2,9
FITEM,2,10
FITEM,2,11   
FITEM,2,12   
AL,P51X
FLST,2,2,5,ORDE,2
FITEM,2,1   
FITEM,2,2
VROTAT,P51X, , , , , ,10,11 ,90, ,
VSEL,ALL
VGLUE,ALL
NUMMRG,ALL
NUMCMP,ALL
ET,1,SOLID98
!ET,2,SOLID98
!定义黄铜材料参数
mat,1
ET,1,SOLID98
MP,EX,1,10.8E10                  ! bronze MODULUS OF ELASTICITY
MP,NUXY,1,.34                  ! bronze POISSON RATIO   
MP,DENS,1,8900                   ! bronze DENSITY
!定义PZT－4（来自ANSYS VM175）压电振子特性参数
mat,2
ET,2,SOLID98
MP,DENS,2,7600                     ! DENSITY
MP,PERX,2,7.124E-9               ! PERMITTIVITY (X AND Y DIRECTION)MP,PERZ,3,5.841E-9
!PERMITTIVITY (Z DIRECTION)
TB,PIEZ,2                        ! DEFINE PIEZ. TABLE
TBDATA,16,12.7                     ! E61 PIEZOELECTRIC CONSTANT
TBDATA,14,12.7                     ! E52 PIEZOELECTRIC CONSTANT
TBDATA,3,-5.2                      ! E13 PIEZOELECTRIC CONSTANT
TBDATA,6,-5.2                      ! E23 PIEZOELECTRIC CONSTANT
TBDATA,9,15.1                      ! E33 PIEZOELECTRIC CONSTANT
TB,ANEL,2                        ! DEFINE STRUCTURAL TABLE
TBDATA,1,13.9E10,7.43E10,7.78E10   ! INPUT MATRIX   
TBDATA,7,13.9E10,7.78E10
TBDATA,12,11.5E10   
TBDATA,16,3.06E10
TBDATA,19,2.56E10
TBDATA,21,2.56E10
!vglue,all
!VSEL, , , ,      1
vplot
!FLST,5,4,6,ORDE,5   
!FITEM,5,1
!FITEM,5,3
!FITEM,5,5
!FITEM,5,7
VSEL,s , , , 1
VATT,1, ,1,       0
VSEL,s , , , 2
VATT,2, ,   2,       0
vsel,all
smrt,on
esize,3
vmesh,all
smrt,off
!nsel,s,loc,z,0
!cp,11,volt,all
!*get,ntop1,node,,num,min
asel,s,loc,x,0
da,all,SYMM
asel,s,loc,y,0
da,all,SYMM
asel,s,loc,z,0
da,all,SYMM
alls
asel,s,,,3
aplot
!nsel,s,loc,z,0
nsla,s,1
cp,11,volt,all
*get,ntop1,node,,num,min
nsel,s,,,ntop1
nplot
asel,s,,,13
aplot
!nsel,s,loc,z,0
nsla,s,1
cp,21,volt,all            !在指定节点处生成并修改耦合自由度集
*get,nbot1,node1,,num,min !提取变量0的取值
nsel,s,,,nbot1
nplot
!!!!!!!!瞬态分析开始

/solu
alls
ANTYPE,TRANS       !指定分析类型为瞬态动力分析
trnopt,full      !分析方法为完全法      
timint,on,all      !打开时间积分效果
betad,1.1e-5       !设置阻尼系数beta
alphad,18.672
NSUBST,1, , ,1   !子步数设置
!载荷步数
!循环圈数
NQUAN=75
NDIAN=4
NT=NQUAN*NDIAN
f=6216!频率
zhouqi=1/f!周期
!DT=zhouqi/NDIAN
!con=360/NDIAN
!Ni=0
!NDO=NT+10
*do,Ni,1,NT
i101=abs(cos(3.14159265*Ni/2))
!i102=Ni-1
*IF,i101,LT,0.1,then
i102=Ni
t101=(Ni-1)*zhouqi/2+Ni*zhouqi/2*0.001
c=100*cos(3.14159265*(i102-1)/2)
TIME,t101
*else
i102=Ni-1
t101=Ni*zhouqi/2
TIME,t101
c=150*cos(3.14159265*(i102-1)/2)
*endif
!c=100*cos(3.14159265*i102)
d,ntop1,volt,c          !加载电压
d,nbot1,volt,0
TINTP,,0.25,0.5,0.5,,,,,
! DELTA=0.5, and THETA=0.5
solve
*enddo

[ 本帖最后由 kadeli 于 2007-5-25 10:08 编辑 ]

ykczj 发表于 2008-11-11 12:45:03

好贴怎么没人顶呢

judishi 发表于 2008-11-15 15:21:14

应该是很不错的，注释更多一点就好了，
谢谢提供分享

lizhanhui 发表于 2008-12-4 16:32:34

谢谢提供分享!!!!

everyeffort 发表于 2009-12-24 10:07:32

谢谢分享！可是我的是超声振动！

Josephysis 发表于 2009-12-25 10:55:19

非常感谢，目前我也是做压电分析，难度不小啊，资料只有一点一点找

pinychen2006 发表于 2012-5-29 15:32:44

感谢楼主的分享，顶！

ashuiamstrong 发表于 2012-6-1 17:20:19

楼主的注释比较少
页: [1]
查看完整版本: 原创:Cymbal压电振子瞬态分析