找回密码
 注册
Simdroid-非首页
查看: 268|回复: 1

[压电] 换能器瞬态分析 求大神讲解

[复制链接]
发表于 2017-4-18 16:44:58 | 显示全部楼层 |阅读模式 来自 天津
换能器分析
一直出错,
下面是我的命令流

/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  
 楼主| 发表于 2017-4-18 16:48:17 | 显示全部楼层 来自 天津
Simdroid开发平台
求解那弄错了,是这样的
/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  
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

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

GMT+8, 2026-1-7 04:49 , Processed in 0.024193 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

快速回复 返回顶部 返回列表