换能器瞬态分析 求大神讲解
换能器分析
一直出错，
下面是我的命令流
/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)
allsel,all
solve
*enddo
finish
http://forum.simwe.com/static/image/common/fav.gif
求解那弄错了，是这样的
/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
页:
[1]