ANSYS压电片建模

zltmsn

那位用ANSYS对压电传感器建模过？？大家讨论下。
以下是我自己写的命令流：
FINISH
/CLEAR
/FILNAME,XX
FINISH
/CLEAR,NOSTART
/PREP7
ET,1,SOLID45
MP,EX,1,1.5E11 !完好的弹性模型
MP,EY,1,8.8E9
MP,EZ,1,8.8E9
MP,GXY,1,4.5E9
MP,GYZ,1,4.5E9
MP,GXZ,1,4.5E9
MP,PRXY,1,0.3
MP,PRYZ,1,0.025
MP,PRXZ,1,0.025
MP,DENS,1,7800
BLC4,,,1.2,2.4,4
VSWEEP,1
!!====================
!!加约束
!!====================
ASEL,,,,2
NSLA
D,ALL,ALL,0
ALLSEL
！定义压电作动器
ET,2,SOLID5,0 !------定义5号单元
keyopt,2,1,3 !开启压电自由度选项
MP,DENS,2,7.64E3 !------定义密度
MP,NUXY,2,0.27 !-----定义泊松比
MP,PERX,2,893 !-----定义相对介电常数
MP,PERY,2,893
MP,PERZ,2,654
TB,ANEL,2 !-----定义PZT刚度矩阵
TBTEMP,0
TBDATA,,1.26E+012,7.95E+011,8.41E+011,0,0,0
TBDATA,,1.26E+012,8.41E+011,0,0,0,1.17E+012
TBDATA,,0,0,0,2.33E+011,0,0
TBDATA,,2.3E+011,0,2.3E+011,,,  
TB,PIEZ,2,,,0   !-----定义压电陶瓷压电矩阵
TBDATA,1,0,,-12.3 !-----定义压电陶瓷压电应力矩阵（在Z轴极化）
TBDATA,4,0,,-12.3
TBDATA,7,0,,20.2
TBDATA,10,,17.5,0
TBDATA,13,17.5,0,0
TBDATA,16,0,0,0
EGEN,2,500,178,,,1,1,,,,,0.4
ESEL,,,,178,181,3
NSLE
CPINTF,UX,0.03
CPINTF,UY,0.03
CPINTF,UZ,0.03
ALLSEL
EGEN,2,500,54,,,1,1,,,,,0.4
ESEL,,,,54,182,128
NSLE
CPINTF,UX,0.03
CPINTF,UY,0.03
CPINTF,UZ,0.03
ALLSEL
SAVE
save
FINISH
/SOLU
ANTYPE,4
TRNOPT,FULL
NLGEOM,On !-----关闭大变形选项
OUTPR,ALL,ALL
OUTRES,ALL,ALL
TIMINT,ON
AUTOTS,ON
T_END=10!------定义载荷向量长度
TINTP,,0.25,0.5,0.5 !-----
!当3=0.005时，计算出a= 0.0496,b= 4.76e-4
!当3=0.008时，计算出a= 0.0828,b= 7.1e-4
!当3=0.01 时，计算出a= 0.0605,b= 1.56E-3
!当3=0.02 时，计算出a= 0.1694,b= 2.36e-3
ALPHAD,0, !-----定义质量阻尼
BETAD,0,  !-----定义刚度阻尼
*DIM,V_H,ARRAY,T_END !-----定义载荷向量1
V_H(1)=1,-1,1,-1,1,-1,1,-1,1,-1
*DIM,V_L,ARRAY,T_END !-----定义载荷向量2
V_L(1)=0,0,0,0,0,0,0,0,0,0
T_INT=0.1 !-----定义时间步长
*DO,I,1,T_END,1
NSUBST,10!-----定义载荷步数目
KBC,1 !-----设置阶跃加载方式
TIME,I*T_INT !-----时间步长
NSEL,,,,599,700,101!压电片上表面
NSEL,A,,,622,639,17
CM,N_H,NODE
D,N_H,VOLT,0.05*V_H(I,1) !-----加方波形式的电压信号
ALLSEL
NSEL,,,,517,807,290 !压电片下表面
NSEL,A,,,590,624,34
CM,N_L,NODE
D,N_L,VOLT,0.05*V_L(I,1) !-----接地
ALLSEL
LSWR,I !-----把第I个载荷数据写入文件
ALLSEL
*ENDDO
ALLSEL
LSSOLVE,1,T_END


这个命令流是模拟压电片作为作动器和传感器的，下半部分是瞬时动力学分析。但是作动器的功能实现了，传感器的电信号我一直得不到。那位高人给我指点一下

ansyslong

不知道你解决了没有，帮你顶一下，同关注。

vipking

NSLE
CPINTF,UX,0.03
CPINTF,UY,0.03
CPINTF,UZ,0.03


不好意思啊，刚刚接触ANSYS，求教楼主
这几行是怎么耦合的？我看效果好像是不是把重合的点给耦合了？那个0.03是哪里的参数啊

zltmsn

我这里是选择容差在0.03范围内的节点进行耦合的，重合的节点之前我已经融合为一个节点了，因此耦合的节点应该都不是同一个节点

享受生活

zltmsn 发表于 2010-7-29 15:51

                               
登录/注册后可看大图

我这里是选择容差在0.03范围内的节点进行耦合的，重合的节点之前我已经融合为一个节点了，因此耦合的节点应 ...

怎么搞压电的人这么少啊，我也想知道啊，初学者啊