电压命令流是怎么样的，

1060295719

我的电压怎么一直是0呢

senja

建议楼主具体详细点描述问题。如果把你自己的命令以文件的形式传一份上来就更好了

1060295719

这是我的命令流，我现在存在几个疑惑，请亲给我解答，谢谢。我本以此图建立模型。在我的命令流想以串联形式连接上下压电晶体。但是，我现在考虑上下压电晶体极化是否正确，电路连接是否正确（就是定义上下极和中间连接）。能否告知如何是串联，如何是并联？请帮我检查一下命令流并给予回答。谢谢.
FINISH
/CLEAR,START
! 定义压电复合梁几何模型
L=50.8e-3
W=31.8e-3
H =0.26e-3    !压电片几何尺寸
H1=0.14e-3    !主结构几何尺寸
H2=3.5e-3      !物块
L1=6e-3
*SET,R1,1E-10 ! 电阻值取足够时，相当于导线
*SET,R2,1000

/prep7
!定义材料类型
et,1,solid5,3               !定义压电单元
et,2,solid45                !定义主结构单元
ET,3,CIRCU94             ! 电路单元(电阻)
R,1,R1,
R,2,R2,
! PZT-5H 材料特性参数
mp,DENS,1,7800                 ! 密度, kg/m^3
EMUNIT,EPZRO,8.854E-12       !介电常数
MP,PERX,1,1.730E3  
MP,PERY,1,1.730E3  
MP,PERZ,1,1.500E3
tb,ANEL,1                      ! 弹性劲度系数, N/m^2
tbdata,1,12.1E10,7.54E10,7.52E10  ! c11,c12,c13
tbdata,7,12.1E10,7.52E10      ! c22,c23
tbdata,12,11.1E10              ! c33
tbdata,16,2.11E10              ! c44
tbdata,19,2.11E10               ! c55
tbdata,21,2.26E10               ! c66
tb,PIEZ,1                      ! 压电（应力）常数, C/m^2
TBDATA,3,-9.1912
TBDATA,6,-9.1912  
TBDATA,9,12.938  
TBDATA,11,12.32  
TBDATA,13,12.32           
!定义主结构的材料参数
mp,dens,2,9000
mp,EX,2,105e9
mp,nuxy,2,0.33
!定义质量块的材料参数
mp,dens,3,8984.73
mp,EX,3,105E9
mp,nuxy,3,0.33
local,11                ! 建立下层局部坐标 +Z 方向
local,12, , , , , , ,180! 建立上层局部坐标 -Z 方向
local,13,,L,0,-H2,,,,
csys,11                 ! 激活局部坐标系11 +y 方向
block,0,L,0,W,O,H1
block,0,L,0,W,0,-H
block,0,L,0,W,H1,H1+H
block,L,L+L1,0,W,-H2,H2
vglue,all                  
NUMCMP,all
!进行网格划分
LESIZE,ALL,1e-3, , , ,1, , ,1,
mat,2 $ type,2 $ esys,11  !对中间结构网格划分
Vsweep,3
allsel
mat,3 $ type,2 $ esys,11
vsweep,4
allsel
mat,1 $ type,1 $ esys,11  !对上层压电片网格划分设置元素的材料属性指针
Vmesh,1
allsel
mat,1 $ type,1 $ esys,12  !对下层压电片网格划分
vmesh,2
allsel
NUMMRG,NODE
NSEL,ALL
NUMCMP,NODE
nsel,s,loc,z,-H           !定义下层电极
nsel,r,loc,x,0,L
*get,n_bot,node,0,num,min
cp,n_bot,volt,all
allsel
nsel,s,loc,z,H+H1         !定义上层电极
nsel,r,loc,x,0,L
*get,n_top,node,0,num,min
cp,n_top,volt,all
allsel
ASEL,S,LOC,Z,0    !中间导线连接，接地     
nsel,R,loc,z,0           
nsel,r,loc,x,0,L
*GET,CM_3,NODE,,NUM,min
allsel
ASEL,S,LOC,Z,0        
nsel,R,loc,z,0           
nsel,r,loc,x,0,L
cp,CM_3,volt,all
d,CM_3,volt,0
allsel
ASEL,S,LOC,Z,H1
nsel,R,loc,z,H1
nsel,r,loc,x,0,L
*GET,Cm_4,NODE,,NUM,min
cp,CM_4,volt,all
allsel
TYPE,3
REAL,1
E,CM_3,Cm_4
allsel
TYPE,3
REAL,2
E,n_top,n_bot
allsel
!定义约束
ASEL,S,LOC,X,0
NSEL,R,LOC,X,0
NSEL,R,LOC,Z,0,H1
D,ALL,UX,,,,,UY,UZ
ALLSEL,ALL
!谐响应
/SOLU
ANTYPE,HARMIC,NEW !谐响应分析
HROPT,FULL ! 完全法
HROUT,OFF ! 输出结果为幅值和相角
OUTPR,BASIC,1
NSUBST,50 ! 指定载荷步中所需要的子步数
HARFRQ,35,60 ! 频率范围
KBC,1 ! 选择加载方式
DMPRAT,0.027 !结构阻尼
ACEL,0,0,46000,  !加速度
SOLVE
/POST26
NSOL,2,3215,U,Z,displace     !位移
/AXLAB,Y,displace
PLVAR,2,
DERIV,3,2,1, ,vel/g, , ,1/10,   !速度
/AXLAB,Y,vel/g
/XRANGE,35,60   
/YRANGE,0,1,1
PLVAR,3
nsol,4,n_top,volt
PLVAR,4
FINI

1060295719

附件，理论短路与开路频率分别为45.54和48.5Hz，我也不知道哪里错了，请帮忙指正