求问：pzt-4压电圆管在水中模拟命令流检查

dengnigui

/PREP7  
et,1,plane13
et,2,fluid29
et,3,fluid29
et,4,fluid129

KEYOPT,1,1,7
KEYOPT,1,2,0
KEYOPT,1,3,1
KEYOPT,1,4,0
KEYOPT,1,5,0
KEYOPT,2,2,0
KEYOPT,2,3,1
!*  
KEYOPT,3,2,1
KEYOPT,3,3,1
!*  
KEYOPT,4,3,1

!压电材料参数
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,7500
TB,ANEL,1,1,21,0
TBTEMP,0
TBDATA,,11.5e10,7.43e10,7.43e10,0,0,0   
TBDATA,,13.9e10,7.78e10,0,0,0,13.9e10   
TBDATA,,0,0,0,2.56e10,0,0   
TBDATA,,3.06e10,0,2.56e10,,,
TB,PIEZ,1,,,0   
TBMODIF,1,1,15.1
TBMODIF,1,2,
TBMODIF,1,3,
TBMODIF,2,1,-5.2
TBMODIF,2,2,
TBMODIF,2,3,
TBMODIF,3,1,-5.2
TBMODIF,3,2,
TBMODIF,3,3,
TBMODIF,4,1,
TBMODIF,4,2,12.7
TBMODIF,4,3,
TBMODIF,5,1,
TBMODIF,5,2,
TBMODIF,5,3,
TBMODIF,6,1,
TBMODIF,6,2,
TBMODIF,6,3,12.7
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,PERX,1,,635  
MPDE,PERY,1
MPDATA,PERY,1,,370  
MPDE,PERZ,1
MPDATA,PERZ,1,,370  
SAVE

MP,DENS,2,1000   
MP,SONC,2,1500
MP,MU,2,0  

MP,DENS,3,1000  
MP,SONC,3,1500  
MP,MU,3,1
!实常数
R,1,0.5,0,0,

RECTNG,0.036,0.042,-0.019,0.019
cyl4,,,0,-90,0.5,90
ASBA,2,1,,DELETE,DELETE
RECTNG,0.036,0.042,-0.019,0.019
NUMCMP,ALL
NUMMRG,ALL,,,,LOW
SAVE

csys,1
k,,0.5,30,
k,,0.5,-30,
csys,0
lstr,4,6
lstr,1,5
lstr,8,3
lstr,9,2
lstr,7,4
lstr,7,1
lsel,s,line,,8,13
cm,L1,line
asbl,2,l1
alls

!给面分配属性
ASEL,S,,,1
AATT,1,0,1,0

ASEL,S,,,3
ASEL,A,,,4
ASEL,A,,,5
ASEL,A,,,6
ASEL,A,,,7
ASEL,A,,,8
AATT,2,0,3,0

!划分线
LESIZE,1,,,20,1
LESIZE,2,,,100,1
LESIZE,3,,,20,1
LESIZE,4,,,100,1
LESIZE,6,,,100,-10
LESIZE,7,,,100,-10
LESIZE,8,,,100,10
LESIZE,9,,,100,10
LESIZE,10,,,100,10
LESIZE,11,,,100,10
LESIZE,12,,,100,-10
LESIZE,13,,,100,-10

!划分面不建议用命令
ASEL,S,,,3
ASEL,A,,,4
ASEL,A,,,5
ASEL,A,,,6
ASEL,A,,,7
ASEL,A,,,8
MSHKEY,1
MSHAPE,0,2D
ASEL,S,,,1
AMESH,3,8
alls
!下面划分面不建议用命令
MSHKEY,1
MSHAPE,0,2D
ASEL,S,,,1
AMESH,1
alls


!改变流固接触面单元属性
LSEL,S,,,1
lsel,a,,,2
LSEL,A,,,3
LSEL,A,,,4  
NSLL,S,1  
ESLN,S     
ESEL,u,TYPE,,1
NSLE,S
EMODIF,ALL,TYPE,2
alls

!创建fluid129单元
CSYS,1
LSEL,S,LOC,X,0.5
nsll,s,1
TYPE,4
MAT,2
REAL,1
ESURF
CSYS,0
alls

!创建流固耦合标识符
LSEL,S,,,1
lsel,a,,,2
lsel,a,,,3
lsel,a,,,4
NSLL,S,1
ESEL,S,TYPE,,2
SF,ALL,FSI,1
alls

!施加载荷以及耦合
lsel,s,,,2
nsll,s
CP,1,VOLT,ALL   
*GET,N1,NODE,0,NUM,MIN  
D,N1,VOLT,10

lsel,s,,,4
nsll,s
CP,2,VOLT,ALL   
*GET,N2,NODE,0,NUM,MIN  
D,N2,VOLT,0

!求解
/solu
antype,harmic
hropt,full
harfrq,5000,30000
nsubst,100
kbc,1
DMPRAT,0.042


!后处理
!导纳曲线
!显示实部G和虚部B
/POST26
RFORCE,2,N1,AMPS   
CFACT,,,,1  
PROD,3,2,1,,MHOS,,,-2*3.1416/10   ! 计算导纳Y=-jwQ/V ,10指的是电压
realvar,4,3,,,YR
imagin,5,3,,,YI
PLCPLX,2       !显示实部
PRCPLX,0       !实部+虚部输出
PLVAR,4,5

!画导纳圆
/POST26
RFORCE,2,N1,AMPS   
CFACT,,,,1  
PROD,3,2,1,,MHOS,,,-2*3.1416   ! 计算导纳Y=-jwQ/V
realvar,4,3,,,YR
imagin,5,3,,,YI
xvar，4
PLVAR,5

dengnigui

我用这个命令进行计算，可以运行，但是后处理里面我不知道该提取那些参数，怎么提取，莫老师的书我看过，但不会操作，大家帮我看看我的命令流正确不！！谢谢！！

dengnigui

问题描述：径向极化的圆管在水中，取其截面来进行模拟，从三维变到二维，减少运算次数。

dengnigui

求版主回复！！谢谢

dengnigui

求版主回复！！谢谢

dengnigui

请求hy1966 老师回答？谢谢了

hy1966

dengnigui 发表于 2011-11-4 17:29
我用这个命令进行计算，可以运行，但是后处理里面我不知道该提取那些参数，怎么提取，莫老师的书我看过，但 ...

这个命令好像没法运行，划分水面的时候有错误，当然你注明了“划分面不建议用命令”，所以没法运行下去
步骤看起来是那么回事，在水的最外层设置129单元，把仅靠辐射面的水单元变成有结构属性的单元，并设置流固耦合标志，然后施加负载求解。
后处理部分提取导纳，作导纳圆图也是正确的，当然既然在水中模拟，一般总想知道声学特性吧，莫老师的书已经介绍的比较清楚了

dengnigui

hy1966 发表于 2011-11-13 09:10
这个命令好像没法运行，划分水面的时候有错误，当然你注明了“划分面不建议用命令”，所以没法运行下去
...

非常感谢您的回复！对于径向极化的圆管换能器，在分析的时候分析它的一个截面，这时候需不需要建立单元坐标系为柱坐标系啊？

hy1966

dengnigui 发表于 2011-11-16 10:45
非常感谢您的回复！对于径向极化的圆管换能器，在分析的时候分析它的一个截面，这时候需不需要建立单元坐 ...

径向极化圆管换能器，使用2-D轴对称模型时，不需要设置特别的坐标系，有时为了用半径、角度选择实体，也设置系统坐标系为柱坐标系，只是Z=0.

dengnigui

hy1966 发表于 2011-11-19 08:55
径向极化圆管换能器，使用2-D轴对称模型时，不需要设置特别的坐标系，有时为了用半径、角度选择实体，也 ...

恩  谢谢您 ！！这几天一直没上网

tzhdacy

转2D模型记住需要打开轴对称指标

419321723

dengnigui 发表于 2011-11-4 17:32
问题描述：径向极化的圆管在水中，取其截面来进行模拟，从三维变到二维，减少运算次数。 ...

从三维编导二维，为什么谐振频率会变呢？

419321723

hy1966 发表于 2011-11-19 08:55
径向极化圆管换能器，使用2-D轴对称模型时，不需要设置特别的坐标系，有时为了用半径、角度选择实体，也 ...

hy1966老师，建模时从从三维变到二维，属性什么都不变，为什么谐振频率会变呢？而且差别很大？

hy1966

因为设置了轴对称，模态少了
同一模态变化不大