怎么更改这个命令流使接收到信号随弹性模量的增加而增加
FINISH/UNITS,SI
/CLEAR
/CONFIG,NRES,200000
/PREP7 !进入前处理器
/TITLE,piezoelectric
ET,1,SOLID185 !定义单元类型
ET,2,SOLID226,1001 !PZT-5A
!!!!设置水泥材料参数
MP,EX,1,30E9!!!!!!！！！！！！！！！！弹性模量！！！！！！！！！
MP,DENS,1,2400 !定义1号单元材料属性
MP,PRXY,1,0.3
MP,BETD,1,4E-8 !定义水泥材料阻尼
!!!!设置压电片材料参数
MP,DENS,2,7800 !定义压电片相关参数PZT-5A
MP,PRXY,2,0.35!定义压电陶瓷材料泊松比
EMUNIT,EPZRO,8.85E-12!定义真空介电常数
MP,PERX,2,916 !定义介电常数 k11=916,k33=830
MP,PERY,2,916 !相对介电常数k11(Z轴极化)
MP,PERZ,2,830 !k33
!定义压电应力系数
e33=15.8 !单位为c/m2
e31=-5.4
e15=12.3
!定义各向异性弹性刚度
c11=12.1E10 !单位为N/m2
c12=7.54E10
c13=7.52E10
c33=11.1E10
c44=2.11E10
c66=2.26E10
TB,PIEZ,2 !创建压电应力系数
TBDATA,3,e31 !e31
TBDATA,6,e31 !e31
TBDATA,9,e33 !e33
TBDATA,14,e15 !e15
TBDATA,16,e15 !e15
TB,ANEL,2,,,0 !创建各向异性弹性刚度
TBDATA,1,c11,c12,c13 !C11,C12,C13
TBDATA,7,c11,c13 !C11,C13
TBDATA,12,c33 !C33
TBDATA,16,c66 !C66
TBDATA,19,c44 !C44
TBDATA,21,c44 !C44
BLOCK,0,0.1,0,0.1,0,0.1!!!V1
BLOCK,0.044,0.056,0.044,0.056,0.03,0.0308!!!V2
BLOCK,0.044,0.056,0.044,0.056,0.0692,0.07!!!V3
VSBV,1,2,,DELETE,KEEP
VSBV,4,3,,DELETE,KEEP
VGLUE,ALL
!切分
wpoff,0.044,0,0 !!移动工作平面
wpro,,,90.000000 !!!旋转工作平面，绕z轴,x,y...对着看逆时针为正
VSBW, all!!!!按工作平面切分
wpoff,0,0,0.012
VSBW, all
wpoff,0,0.044,0
wpro,,-90.000000,
VSBW, all
wpoff,0,0,0.012
VSBW, all
WPCSYS,-1,0 !!!工作平面复位
wpoff,0,0,0.03
VSBW, all
wpoff,0,0,0.0008
VSBW, all
wpoff,0,0,0.0384
VSBW, all
wpoff,0,0,0.0008
VSBW, all
WPCSYS,-1,0
VSEL,S,VOLU,,2,3
VATT,2,,2
VSEL,s,VOLU,,4,53
VSEL,r,VOLU,,1!选择
VATT,1,,1!!!定义单元属性
VSEL,ALL
!压电片划分
lsel,s,loc,x,0.1/2-0.006+0.00001,0.01/2+0.006-0.00001
lsel,a,loc,y,0.1/2-0.006+0.00001,0.01/2+0.006-0.00001
lesize,all,,,10, ,1, , ,1,
allsel,all
lsel,s,loc,z,0.03+0.00001,0.0308-0.00001
lsel,a,loc,z,0.0692+0.00001,0.07-0.00001
lesize,all,,,1,,1, , ,1,
allsel,all!!!显示几何元素
!水泥划分
lsel,s,loc,x,0.00001,0.1/2-0.006-0.00001
lsel,a,loc,x,0.1/2+0.006+0.00001,0.1-0.00001
lsel,a,loc,y,0.00001,0.1/2-0.006-0.00001
lsel,a,loc,y,0.1/2+0.006+0.00001,0.1-0.00001
lesize,all,,,5, ,1, , ,1,
allsel,all
lsel,s,loc,z,0.00001,0.03-0.00001
lsel,a,loc,z,0.0308+0.00001,0.0692-0.00001
lsel,a,loc,z,0.07+0.00001, 0.1-0.00001
LESIZE,ALL, , ,10, ,1, , ,1,
allsel,all
VSEL,ALL !网格划分
MSHKEY,1
MSHAPE,0,3D
VMESH,ALL
!施加约束
!前后约束
asel,s,loc,z,0
asel,a,loc,z,0.1
nsla,r!!!选择与所选面相关的节点,R---从当前子集中再选择一组面，形成新的当前子集
cm,yueshu,node
d,yueshu,all,
allsel,all
!左右约束
asel,s,loc,x,0
asel,a,loc,x,0.1
nsla,r
cm,yueshu1,node
d,yueshu1,all,
allsel,all
*SET,V1,0
NSEL,S,LOC,Z,0.0308 !耦合第一个压电片的电压
NSEL,R,LOC,X,0.044,0.056
NSEL,R,LOC,Y,0.044,0.056
CP,1,VOLT,ALL !耦合电压
*GET,N1,NODE,,NUM,MIN
D,N1,VOLT,V1 !施加电压荷载
NSEL,S,LOC,Z,0.03
NSEL,R,LOC,X,0.044,0.0506
NSEL,R,LOC,Y,0.044,0.0506
CP,2,VOLT,ALL
*GET,N2,NODE,,NUM,MIN
NSEL,S,LOC,Z,0.0692 !耦合第二个压电片的电压
NSEL,R,LOC,X,0.044,0.056
NSEL,R,LOC,Y,0.044,0.056
CP,3,VOLT,ALL
*GET,N3,NODE,,NUM,MIN
D,N3,VOLT,V1 !施加电压荷载
NSEL,S,LOC,Z,0.07
NSEL,R,LOC,X,0.044,0.056
NSEL,R,LOC,Y,0.044,0.056
CP,4,VOLT,ALL
*GET,N4,NODE,,NUM,MIN
ALLSEL,ALL
!施加表荷载
*dim,V2,table,5001,1,1 !定义表数组,为5001行,1列,1平面
!*dim,V2,table,5001,1,1,time,voltage1,voltage2
!横坐标数据标签为时间,纵坐标为电压(数据标签绝不能和表名相同，即voltage1不能定义为dianya)
!若定义了参数前两个参数，第三个参数必须定义(若定义time,voltage1，必定义voltage2 )
*tread,V2,chen1,txt,, !读入表数据！！chen1直接放入工作目录即可
*vplot,V2(1,0),V2(1,1)
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 !使用瞬态效应
time,0.2 ! 第一载荷步
nsubst,5001,5001,5001
ALLSEL,ALL
D,N2,volt,%V2%
KBC,1 !定义荷载不增加形式
SOLVE
FINISH
一般贴关键代码就好
页:
[1]