压电材料在transient下拟真的震动量为何会那么小呢?
浏览：100713 回答：1
各位高手
小弟用transient拟真超声波压电马达的transient,
由于轴向(z轴向)所获得震动量仅3.5E-9 mm,与实验上所获得的数值相差很多,
约小了100倍(实际上为1.6~3.5E-7 mm),由于小弟的功力实在太弱,找不出哪里设定错误,所以不知高手可否帮小弟检查一下,看看小弟到底是哪里设定错误或是哪里该设定而未设定,希望有高手可以不吝涉指点小弟,
小弟真是会感激不尽
帖子上仅是部份程序码,附件内有完整程序码
会是因为我在压电片与金属片之间未设定材料连结吗?
或是压电片与金属圆杆之间未设定材料连结吗?
如果是的话,该怎么设定呢?
!此模型为一个超音波压电马达,组成组件为两片圆盘型压电片,一个圆盘型铝金属片(金属片粘在两压电
!片之间),一根铝杆所组合.
!其中压电片与铝材料的材料性质均势根据online help 上的 VM176范例.
!我在第一块压电片的右边加载+20v电压,在第二块压电片的左边加载-20v电压,然后在金属圆片的两侧!!
!加载0 v电压.
!输入电压波型为sin wave,单位为MKS制,Frequency=64840 Hz,希望求得超音波压电马达在TRANSIENT!
!下的震动量,速度,加速度
~PARAIN,tula170m,x_t,,SOLIDS,0,0 !汇入一个叫tula170m,x_t的solidworks档案
/FACET,NORML
/DIST,1,1.08222638492,1
/REP,FAST
FLST,2,4,6,ORDE,2
FITEM,2,1
FITEM,2,-4
VOVLAP,P51X
/UNITS,SI ! 单位
ET,1,SOLID5 !element type 设定，SOLID5是压电材料
ET,2,SOLID45 !element type 设定，SOLID45是两压电片之间的金属片
ET,3,SOLID45 !element type 设定，SOLID45是3D杆件
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7730 !压电材料密度
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,PERX,1,,993.55 !压电材料perx
TB,PIEZ,1,,,0 !压电材料的压电性质，由参考于范例VM176
TBMODIF,1,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,1,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,1,3,-6.1 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,2,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,2,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,2,3,-6.1 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,3,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,3,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,3,3,15.7 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,4,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,4,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,4,3,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,5,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,5,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,5,3,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,6,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,6,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,6,3,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
!*
TB,ANEL,1,1,21,0 !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10 E-10 N/m^2
TBTEMP,0 !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10 E-10 N/m^2
TBDATA,,1.28E011,6.8E010,6.6E010,0,0,0 !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10 E-10 N/m^2
TBDATA,,1.28E011,6.6E010,0,0,0,1.1E011 !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10 E-10 N/m^2
TBDATA,,0,0,0,2.1E010,0,0 !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10 E-10 N/m^2
TBDATA,,2.1E010,0,2.1E010,,, !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,2690 !两压电片之间的铝圆盘的密度，由参考于范例VM176
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,7.03e+10 !两压电片之间的铝圆盘的杨氏系数，由参考于范例VM176
MPDATA,PRXY,2,,0.345 !两压电片之间的铝圆盘的普松比，由参考于范例VM176
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,3,,2690 !铝长杆的密度，由参考于范例VM176
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,3,,7.03e+10 !铝长杆的杨氏系数，由参考于范例VM176
MPDATA,PRXY,3,,0.345 !铝长杆的普松比，由参考于范例VM176
WPSTYLE,,,,,,,,1
/POST26
FINISH
/CLEAR
/COM,ANSYS RELEASE 10.0 UP20050718 17:39:25 09/11/2007
/input,start100,ans,'C:\Program Files\Ansys Inc\v100\ANSYS\apdl\',,,,,,,,,,,,,,,,1
/POST1
INRES,BASIC
FILE,'0911_tula_vm176_transient_64840','rst','.'
SET,LAST
FINISH
/POST26
FILE,'0911_tula_vm176_transient_64840','rst','.'
/UI,COLL,1
NUMVAR,200
SOLU,191,NCMIT
STORE,MERGE
FILLDATA,191,,,,1,1
REALVAR,191,191
/VIEW,1,1
/ANG,1
/REP,FAST
NSLA,S,1
/DIST,1,1.08222638492,1
/REP,FAST
nplot
d,all,ux
d,all,uy
d,all,uz
ALLSEL,ALL
VPLOT
/AUTO,1
/REP,FAST
cmsel,u,al2
/DIST,1,1.08222638492,1
/REP,FAST
cmsel,u,al3
/DIST,1,1.08222638492,1
/REP,FAST
/USER, 1
/VIEW, 1, 0.994024475624 , 0.430303263369E-01, -0.100318158251
/ANG, 1, 9.15554682511
/REPLO
/DIST,1,0.729,1
/REP,FAST
/DIST,1,0.729,1
/REP,FAST
/FOC,1,0.3,,,1
/REP,FAST
/FOC,1,0.3,,,1
/REP,FAST
/FOC,1,0.3,,,1
/REP,FAST
/DIST,1,0.729,1
/REP,FAST
/REPLO
FLST,5,8,5,ORDE,8
FITEM,5,99
FITEM,5,103
FITEM,5,107
FITEM,5,110
FITEM,5,113
FITEM,5,117
FITEM,5,121
FITEM,5,124
ASEL,S, , ,P51X
/DIST,1,1.08222638492,1
/REP,FAST
NSLA,S,1
/DIST,1,1.08222638492,1
/REP,FAST
nplot
FINISH
/SOL
d,all,ux
d,all,uy
d,all,uz
ALLSEL,ALL
VPLOT
/AUTO,1
/REP,FAST
!***************************
!nsel,s,loc,z,0
!cp,63,volt,all
!*get,ntop1,node,,num,min
!nsel,s,loc,z,0.0002
!cp,83,volt,all
!*get, ntop2,node,,num,min
!nsel,s,loc,z,0.0003
!cp,241,volt,all
!*get, ntop3,node,,num,min
!nsel,s,loc,z,0.0005
!cp,553,volt,all
!*get, ntop4,node,,num,min
!/DIST,1,1.08222638492,1
!/REP,FAST
!**************************************
!*
ANTYPE,4 !瞬时分析
!*
TRNOPT,FULL
LUMPM,0
NLGEOM,0
SSTIF,0
NROPT,AUTO, ,
EQSLV, , ,0,
PRECISION,0
MSAVE,0
TOFFST,0,
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'PSIN2'
*SET,_FNCCSYS,0
! /INPUT,PSIN1.func,,,1 !呼叫TPSIN1.func的附檔
*DIM,%_FNCNAME%,TABLE,6,10,1,,,,%_FNCCSYS%
! Begin of equation: 20*SIN(64840*2*3.14*{TIME}) !瞬时分析所所使用的SIN波输入方程式-正电压
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 64840, 0, 0, 0 !64840HZ频率
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 2, 0, 0, -1 !2是方程式里面的
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 3.14, 0, 0, -3 !3.14是pi
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 1, -2, 3, 1
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 9, 1, -1, 0, 0
*SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 20, 0, 0, -1 !20是振幅
*SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,10,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 20*SIN(64840*2*3.14*{TIME}) !瞬时分析所所使用的SIN波输入方程式-正电压
!-->
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'MSIN2'
*SET,_FNCCSYS,0
! /INPUT,MSIN1.func,,,1 !呼叫TMSIN1.func的附檔
*DIM,%_FNCNAME%,TABLE,6,10,1,,,,%_FNCCSYS%
!
! Begin of equation: -20*SIN(64840*2*3.14*{TIME}) !瞬时分析所所使用的SIN波输入方程式-负电压
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0,64840, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 3.14, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 1, -2, 3, 1
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 9, 1, -1, 0, 0
*SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, -20, 0, 0, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,10,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: -20*SIN(64840*2*3.14*{TIME}) !瞬时分析所所使用的SIN波输入方程式-负电压
!-->
/USER, 1
/VIEW, 1, -0.313512525242 , -0.101271202724 , -0.944168438371
/ANG, 1, -142.877426533
/REPLO
cmsel,s,pie2
/DIST,1,1.08222638492,1
/REP,FAST
/PNUM,KP,0
/PNUM,LINE,0
/PNUM,AREA,1
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0
!*
/PNUM,ELEM,0
/REPLOT
!*
/DIST,1,0.729,1
/REP,FAST
/DIST,1,0.729,1
/REP,FAST
/VIEW, 1, -0.332709513458 , -0.116953646678 , -0.935749017730
/ANG, 1, -143.008836487
/REPLO
/VIEW, 1, -0.222645891108 , -0.284052437740E-01, -0.974485479265
/ANG, 1, -142.520128566
/REPLO
FLST,5,8,5,ORDE,8
FITEM,5,63 !选择面积(要输入正电压用的)
FITEM,5,68 !选择面积(要输入正电压用的)
FITEM,5,72 !选择面积(要输入正电压用的)
FITEM,5,76 !选择面积(要输入正电压用的)
FITEM,5,100 !选择面积(要输入正电压用的)
FITEM,5,104 !选择面积(要输入正电压用的)
FITEM,5,108 !选择面积(要输入正电压用的)
FITEM,5,111 !选择面积(要输入正电压用的)
ASEL,S, , ,P51X
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, -0.325053709860 , -0.191269668011 , -0.926151175460
/ANG, 1, -143.228440365
/REPLO
/VIEW, 1, -0.180759442462 , -0.745706314671E-01, -0.980696306143
/ANG, 1, -142.039486003
/REPLO
DA,all,volt,%PSIN2%
ALLSEL,ALL
VPLOT
cmsel,s,pie2
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, 0.580484352477E-01, -0.135661982313E-01, 0.998221587340
/ANG, 1, 155.445074610
/REPLO
FLST,5,8,5,ORDE,8
FITEM,5,64 !选择面积(要输入接地用的)
FITEM,5,69 !选择面积(要输入接地用的)
FITEM,5,73 !选择面积(要输入接地用的)
FITEM,5,77 !选择面积(要输入接地用的)
FITEM,5,101 !选择面积(要输入接地用的)
FITEM,5,105 !选择面积(要输入接地用的)
FITEM,5,109 !选择面积(要输入接地用的)
FITEM,5,112 !选择面积(要输入接地用的)
ASEL,S, , ,P51X
/DIST,1,1.08222638492,1
/REP,FAST
DA,all,volt,0 !所选的面积均给定0电压，就是接地，
ALLSEL,ALL
VPLOT
cmsel,s,pie3
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, 0.283772140130E-01, 0.339233883598E-01, -0.999021489983
/ANG, 1, -144.443788699
/REPLO
!*************************************
ALLSEL,ALL
VPLOT
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/VIEW, 1, -0.125652195359 , 0.118292083422E-01, -0.992003828436
/ANG, 1, -145.492519493
/REPLO
/VIEW, 1, -0.687082056433 , -0.689682939182 , -0.228595037410
/ANG, 1, -168.810091409
/REPLO
cmsel,s,pie3
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, -0.162516996410 , -0.189007171557 , -0.968434052983
/ANG, 1, -135.969666032
/REPLO
/DIST,1,0.729,1
/REP,FAST
/DIST,1,0.729,1
/REP,FAST
/DIST,1,0.729,1
/REP,FAST
!*************************************
FLST,5,8,5,ORDE,8
FITEM,5,79 !选择面积(要输入接地用的)
FITEM,5,83 !选择面积(要输入接地用的)
FITEM,5,86 !选择面积(要输入接地用的)
FITEM,5,90 !选择面积(要输入接地用的)
FITEM,5,114 !选择面积(要输入接地用的)
FITEM,5,118 !选择面积(要输入接地用的)
FITEM,5,122 !选择面积(要输入接地用的)
FITEM,5,125 !选择面积(要输入接地用的)
ASEL,S, , ,P51X
/DIST,1,1.08222638492,1
/REP,FAST
DA,all,volt,0 !所选的面积均给定0电压，就是接地，
ALLSEL,ALL
VPLOT
cmsel,s,pie3
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, 0.259013259927E-01, -0.519267094794E-01, 0.998314949379
/ANG, 1, 120.444088609
/REPLO
FLST,5,8,5,ORDE,8
FITEM,5,16 !选择面积(要输入负电压用的)
FITEM,5,29 !选择面积(要输入负电压用的)
FITEM,5,33 !选择面积(要输入负电压用的)
FITEM,5,61 !选择面积(要输入负电压用的)
FITEM,5,115 !选择面积(要输入负电压用的)
FITEM,5,119 !选择面积(要输入负电压用的)
FITEM,5,123 !选择面积(要输入负电压用的)
FITEM,5,126 !选择面积(要输入负电压用的)
ASEL,S, , ,P51X
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, 0.544737228281E-01, 0.131106423062 , 0.989870556867
/ANG, 1, 120.550708076
/REPLO
/VIEW, 1, 0.413317973760E-01, 0.108713026137 , 0.993213552301
/ANG, 1, 120.458732391
/REPLO
DA,all,volt,%MSIN2%
ALLSEL,ALL
VPLOT
*SET,fre,50000
*SET,bufferlength,20
TIME,20/fre
AUTOTS,on
DELTIM,1/fre/bufferlength
OUTRES,ALL,1
TINTP,,0.25,0.5,0.5
ALPHAD,0.0000001
BETAD,0.0000009
/STATUS,SOLU
SOLVE
SAVE
SAVE
FINISH
小弟用transient拟真超声波压电马达的transient,
由于轴向(z轴向)所获得震动量仅3.5E-9 mm,与实验上所获得的数值相差很多,
约小了100倍(实际上为1.6~3.5E-7 mm),由于小弟的功力实在太弱,找不出哪里设定错误,所以不知高手可否帮小弟检查一下,看看小弟到底是哪里设定错误或是哪里该设定而未设定,希望有高手可以不吝涉指点小弟,
小弟真是会感激不尽
帖子上仅是部份程序码,附件内有完整程序码
会是因为我在压电片与金属片之间未设定材料连结吗?
或是压电片与金属圆杆之间未设定材料连结吗?
如果是的话,该怎么设定呢?
!此模型为一个超音波压电马达,组成组件为两片圆盘型压电片,一个圆盘型铝金属片(金属片粘在两压电
!片之间),一根铝杆所组合.
!其中压电片与铝材料的材料性质均势根据online help 上的 VM176范例.
!我在第一块压电片的右边加载+20v电压,在第二块压电片的左边加载-20v电压,然后在金属圆片的两侧!!
!加载0 v电压.
!输入电压波型为sin wave,单位为MKS制,Frequency=64840 Hz,希望求得超音波压电马达在TRANSIENT!
!下的震动量,速度,加速度
~PARAIN,tula170m,x_t,,SOLIDS,0,0 !汇入一个叫tula170m,x_t的solidworks档案
/FACET,NORML
/DIST,1,1.08222638492,1
/REP,FAST
FLST,2,4,6,ORDE,2
FITEM,2,1
FITEM,2,-4
VOVLAP,P51X
/UNITS,SI ! 单位
ET,1,SOLID5 !element type 设定，SOLID5是压电材料
ET,2,SOLID45 !element type 设定，SOLID45是两压电片之间的金属片
ET,3,SOLID45 !element type 设定，SOLID45是3D杆件
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7730 !压电材料密度
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,PERX,1,,993.55 !压电材料perx
TB,PIEZ,1,,,0 !压电材料的压电性质，由参考于范例VM176
TBMODIF,1,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,1,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,1,3,-6.1 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,2,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,2,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,2,3,-6.1 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,3,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,3,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,3,3,15.7 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,4,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,4,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,4,3,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,5,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,5,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,5,3,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,6,1,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,6,2,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
TBMODIF,6,3,0 !压电材料的压电性质，由参考于范例VM176 piezoelectric MMMMatrix[e] C/m^2
!*
TB,ANEL,1,1,21,0 !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10 E-10 N/m^2
TBTEMP,0 !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10 E-10 N/m^2
TBDATA,,1.28E011,6.8E010,6.6E010,0,0,0 !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10 E-10 N/m^2
TBDATA,,1.28E011,6.6E010,0,0,0,1.1E011 !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10 E-10 N/m^2
TBDATA,,0,0,0,2.1E010,0,0 !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10 E-10 N/m^2
TBDATA,,2.1E010,0,2.1E010,,, !压电材料的压电性质，由参考于范例VM176 stiffness matrix[c]x10
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,2690 !两压电片之间的铝圆盘的密度，由参考于范例VM176
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,7.03e+10 !两压电片之间的铝圆盘的杨氏系数，由参考于范例VM176
MPDATA,PRXY,2,,0.345 !两压电片之间的铝圆盘的普松比，由参考于范例VM176
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,3,,2690 !铝长杆的密度，由参考于范例VM176
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,3,,7.03e+10 !铝长杆的杨氏系数，由参考于范例VM176
MPDATA,PRXY,3,,0.345 !铝长杆的普松比，由参考于范例VM176
WPSTYLE,,,,,,,,1
/POST26
FINISH
/CLEAR
/COM,ANSYS RELEASE 10.0 UP20050718 17:39:25 09/11/2007
/input,start100,ans,'C:\Program Files\Ansys Inc\v100\ANSYS\apdl\',,,,,,,,,,,,,,,,1
/POST1
INRES,BASIC
FILE,'0911_tula_vm176_transient_64840','rst','.'
SET,LAST
FINISH
/POST26
FILE,'0911_tula_vm176_transient_64840','rst','.'
/UI,COLL,1
NUMVAR,200
SOLU,191,NCMIT
STORE,MERGE
FILLDATA,191,,,,1,1
REALVAR,191,191
/VIEW,1,1
/ANG,1
/REP,FAST
NSLA,S,1
/DIST,1,1.08222638492,1
/REP,FAST
nplot
d,all,ux
d,all,uy
d,all,uz
ALLSEL,ALL
VPLOT
/AUTO,1
/REP,FAST
cmsel,u,al2
/DIST,1,1.08222638492,1
/REP,FAST
cmsel,u,al3
/DIST,1,1.08222638492,1
/REP,FAST
/USER, 1
/VIEW, 1, 0.994024475624 , 0.430303263369E-01, -0.100318158251
/ANG, 1, 9.15554682511
/REPLO
/DIST,1,0.729,1
/REP,FAST
/DIST,1,0.729,1
/REP,FAST
/FOC,1,0.3,,,1
/REP,FAST
/FOC,1,0.3,,,1
/REP,FAST
/FOC,1,0.3,,,1
/REP,FAST
/DIST,1,0.729,1
/REP,FAST
/REPLO
FLST,5,8,5,ORDE,8
FITEM,5,99
FITEM,5,103
FITEM,5,107
FITEM,5,110
FITEM,5,113
FITEM,5,117
FITEM,5,121
FITEM,5,124
ASEL,S, , ,P51X
/DIST,1,1.08222638492,1
/REP,FAST
NSLA,S,1
/DIST,1,1.08222638492,1
/REP,FAST
nplot
FINISH
/SOL
d,all,ux
d,all,uy
d,all,uz
ALLSEL,ALL
VPLOT
/AUTO,1
/REP,FAST
!***************************
!nsel,s,loc,z,0
!cp,63,volt,all
!*get,ntop1,node,,num,min
!nsel,s,loc,z,0.0002
!cp,83,volt,all
!*get, ntop2,node,,num,min
!nsel,s,loc,z,0.0003
!cp,241,volt,all
!*get, ntop3,node,,num,min
!nsel,s,loc,z,0.0005
!cp,553,volt,all
!*get, ntop4,node,,num,min
!/DIST,1,1.08222638492,1
!/REP,FAST
!**************************************
!*
ANTYPE,4 !瞬时分析
!*
TRNOPT,FULL
LUMPM,0
NLGEOM,0
SSTIF,0
NROPT,AUTO, ,
EQSLV, , ,0,
PRECISION,0
MSAVE,0
TOFFST,0,
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'PSIN2'
*SET,_FNCCSYS,0
! /INPUT,PSIN1.func,,,1 !呼叫TPSIN1.func的附檔
*DIM,%_FNCNAME%,TABLE,6,10,1,,,,%_FNCCSYS%
! Begin of equation: 20*SIN(64840*2*3.14*{TIME}) !瞬时分析所所使用的SIN波输入方程式-正电压
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 64840, 0, 0, 0 !64840HZ频率
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 2, 0, 0, -1 !2是方程式里面的
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 3.14, 0, 0, -3 !3.14是pi
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 1, -2, 3, 1
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 9, 1, -1, 0, 0
*SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 20, 0, 0, -1 !20是振幅
*SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,10,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 20*SIN(64840*2*3.14*{TIME}) !瞬时分析所所使用的SIN波输入方程式-正电压
!-->
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'MSIN2'
*SET,_FNCCSYS,0
! /INPUT,MSIN1.func,,,1 !呼叫TMSIN1.func的附檔
*DIM,%_FNCNAME%,TABLE,6,10,1,,,,%_FNCCSYS%
!
! Begin of equation: -20*SIN(64840*2*3.14*{TIME}) !瞬时分析所所使用的SIN波输入方程式-负电压
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0,64840, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 3.14, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 1, -2, 3, 1
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 9, 1, -1, 0, 0
*SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, -20, 0, 0, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,10,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: -20*SIN(64840*2*3.14*{TIME}) !瞬时分析所所使用的SIN波输入方程式-负电压
!-->
/USER, 1
/VIEW, 1, -0.313512525242 , -0.101271202724 , -0.944168438371
/ANG, 1, -142.877426533
/REPLO
cmsel,s,pie2
/DIST,1,1.08222638492,1
/REP,FAST
/PNUM,KP,0
/PNUM,LINE,0
/PNUM,AREA,1
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0
!*
/PNUM,ELEM,0
/REPLOT
!*
/DIST,1,0.729,1
/REP,FAST
/DIST,1,0.729,1
/REP,FAST
/VIEW, 1, -0.332709513458 , -0.116953646678 , -0.935749017730
/ANG, 1, -143.008836487
/REPLO
/VIEW, 1, -0.222645891108 , -0.284052437740E-01, -0.974485479265
/ANG, 1, -142.520128566
/REPLO
FLST,5,8,5,ORDE,8
FITEM,5,63 !选择面积(要输入正电压用的)
FITEM,5,68 !选择面积(要输入正电压用的)
FITEM,5,72 !选择面积(要输入正电压用的)
FITEM,5,76 !选择面积(要输入正电压用的)
FITEM,5,100 !选择面积(要输入正电压用的)
FITEM,5,104 !选择面积(要输入正电压用的)
FITEM,5,108 !选择面积(要输入正电压用的)
FITEM,5,111 !选择面积(要输入正电压用的)
ASEL,S, , ,P51X
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, -0.325053709860 , -0.191269668011 , -0.926151175460
/ANG, 1, -143.228440365
/REPLO
/VIEW, 1, -0.180759442462 , -0.745706314671E-01, -0.980696306143
/ANG, 1, -142.039486003
/REPLO
DA,all,volt,%PSIN2%
ALLSEL,ALL
VPLOT
cmsel,s,pie2
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, 0.580484352477E-01, -0.135661982313E-01, 0.998221587340
/ANG, 1, 155.445074610
/REPLO
FLST,5,8,5,ORDE,8
FITEM,5,64 !选择面积(要输入接地用的)
FITEM,5,69 !选择面积(要输入接地用的)
FITEM,5,73 !选择面积(要输入接地用的)
FITEM,5,77 !选择面积(要输入接地用的)
FITEM,5,101 !选择面积(要输入接地用的)
FITEM,5,105 !选择面积(要输入接地用的)
FITEM,5,109 !选择面积(要输入接地用的)
FITEM,5,112 !选择面积(要输入接地用的)
ASEL,S, , ,P51X
/DIST,1,1.08222638492,1
/REP,FAST
DA,all,volt,0 !所选的面积均给定0电压，就是接地，
ALLSEL,ALL
VPLOT
cmsel,s,pie3
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, 0.283772140130E-01, 0.339233883598E-01, -0.999021489983
/ANG, 1, -144.443788699
/REPLO
!*************************************
ALLSEL,ALL
VPLOT
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/VIEW, 1, -0.125652195359 , 0.118292083422E-01, -0.992003828436
/ANG, 1, -145.492519493
/REPLO
/VIEW, 1, -0.687082056433 , -0.689682939182 , -0.228595037410
/ANG, 1, -168.810091409
/REPLO
cmsel,s,pie3
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, -0.162516996410 , -0.189007171557 , -0.968434052983
/ANG, 1, -135.969666032
/REPLO
/DIST,1,0.729,1
/REP,FAST
/DIST,1,0.729,1
/REP,FAST
/DIST,1,0.729,1
/REP,FAST
!*************************************
FLST,5,8,5,ORDE,8
FITEM,5,79 !选择面积(要输入接地用的)
FITEM,5,83 !选择面积(要输入接地用的)
FITEM,5,86 !选择面积(要输入接地用的)
FITEM,5,90 !选择面积(要输入接地用的)
FITEM,5,114 !选择面积(要输入接地用的)
FITEM,5,118 !选择面积(要输入接地用的)
FITEM,5,122 !选择面积(要输入接地用的)
FITEM,5,125 !选择面积(要输入接地用的)
ASEL,S, , ,P51X
/DIST,1,1.08222638492,1
/REP,FAST
DA,all,volt,0 !所选的面积均给定0电压，就是接地，
ALLSEL,ALL
VPLOT
cmsel,s,pie3
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, 0.259013259927E-01, -0.519267094794E-01, 0.998314949379
/ANG, 1, 120.444088609
/REPLO
FLST,5,8,5,ORDE,8
FITEM,5,16 !选择面积(要输入负电压用的)
FITEM,5,29 !选择面积(要输入负电压用的)
FITEM,5,33 !选择面积(要输入负电压用的)
FITEM,5,61 !选择面积(要输入负电压用的)
FITEM,5,115 !选择面积(要输入负电压用的)
FITEM,5,119 !选择面积(要输入负电压用的)
FITEM,5,123 !选择面积(要输入负电压用的)
FITEM,5,126 !选择面积(要输入负电压用的)
ASEL,S, , ,P51X
/DIST,1,1.08222638492,1
/REP,FAST
/VIEW, 1, 0.544737228281E-01, 0.131106423062 , 0.989870556867
/ANG, 1, 120.550708076
/REPLO
/VIEW, 1, 0.413317973760E-01, 0.108713026137 , 0.993213552301
/ANG, 1, 120.458732391
/REPLO
DA,all,volt,%MSIN2%
ALLSEL,ALL
VPLOT
*SET,fre,50000
*SET,bufferlength,20
TIME,20/fre
AUTOTS,on
DELTIM,1/fre/bufferlength
OUTRES,ALL,1
TINTP,,0.25,0.5,0.5
ALPHAD,0.0000001
BETAD,0.0000009
/STATUS,SOLU
SOLVE
SAVE
SAVE
FINISH
shenhaha 
