求助：为什么输出电压是零呢？

lchylinchy

请问一下我这个程序为什么没有电压产生呢？
某个石英英叉放在空气中，在石英英叉的两个叉指上加力，并且移动该力（主要是为了找到该力在什么位置产生的电压最大），但是我做出了的结果电压为零，
建立的模型如图（红线围着的部分是英叉，周围是空气，橙色线部分是加力的位置）
请问是为什么啊？
谢谢
程序如下
fini
/clear
/prep7
/title,quartz volt anslysis
!et,1,PLANE13,7  !结构单元类型
et,1,PLANE223,1001
et,2,29         !流体单元类型与结构接触
et,3,129   !吸收边界
et,4,29,,1,0      !流体单元类型与结构没有接触 在xyz方向没有位移
r,3,6e-3,0,0
!press=1e-8     !力的大小
C11=86.74e3                       
C12=6.99e3  
C13=11.91e3
C14=17.91e3
C33=107.2e3
C44=57.94e3                    
                     
tb,ANEL,1,,,1              ! 弹性系数矩阵（tb定义非线性材料单元表，同tbdata一起连用）
!最后一个1表示输入的弹性矩阵为弹性形式，如果是0表示输入的弹性矩阵为刚度形式
tbda,1,C11,C12,C13,0,C14
tbda,7,C11,C13,0,-C14
tbda,12,C33
tbda,16,(C11-C12)/2,0,C14
tbdata,19,C44
tbdata,21,C44
tb,PIEZ,1,,,0                  ! Piezoelectric strees matrix
tbda,1,0.171
tbda,4,-0.171
tbda,11,-0.171
tbda,13,-0.0406
tbda,17,0.0406
tb,DPER,1,,,0                   ! Permittivity（介电常数） at constant stress
tbdata,1,4.43,4.43
!tblist,all
mp,dens,1,2.649e-15       !石英密度
mp,dens,2,2.35E-3    !空气密度
MP,SONC,2,1100       !空气声速
mp,mu,2,0.75           !石英英叉的吸收系数
mp,dens,3,2.35E-3    !空气密度
MP,SONC,3,1100       !空气声速
!结构模型
blc4,-0.75e-3,-2.7e-3,1.5e-3,1.8e-3
blc4,-0.75e-3,-0.9e-3,0.6e-3,3.6e-3
blc4,0.15e-3,-0.9e-3,0.6e-3,3.6e-3
aadd,all
cyl4,0,0,6e-3,0,,360
aovlap,1,4
!esize,0.1e-3,0,  !指定单元划分大小
lesize,1,,,15
lesize,15,,,3
lesize,7,,,6
lesize,11,,,6
lesize,4,,,18
lesize,2,,,18
lesize,14,,,36
lesize,13,,,36
lesize,12,,,36
lesize,6,,,36
mshape,0,2d
mshkey,0
!smrt,1
asel,s,,,4
aatt,1,,1
amesh,all
asel,s,,,2
aatt,2,,2
smrt,1
mshkey,0
mshape,0,2d
amesh,all
save,pingm122_mesh,db
!将不与结构相连的单元类型改为4
esel,s,type,,1
nsle,s
esln,s,0
nsle,s
esel,inve
nsle,s
emodif,all,type,4
esel,all
nsel,all
!定义吸收边界
csys,1
nsel,s,loc,x,6e-3
type,3
real,3
mat,2
esurf
esel,all
nsel,all
csys,0
!压电耦合
!定义上层电极
lsel,s,,,6,12,6
NSLL,S,1
cp,1,volt,all
*get,n_top,node,0,num,min
!d,n_top,volt,0
esel,all
nsel,all
!定义下层电极
nsel,s,loc,x,-0.75e-3
nsel,r,loc,y,-0.9e-3,2.7e-3
cp,2,volt,all
*get,n_bot1,node,0,num,min
d,n_bot1,volt,0
esel,all
nsel,all
nsel,s,loc,x,0.75e-3
nsel,r,loc,y,-0.9e-3,2.7e-3
cp,3,volt,all
*get,n_bot2,node,0,num,min
d,n_bot2,volt,0
esel,all
nsel,all
!定义流固界面
csys,0
nsel,s,loc,x,-0.75e-3,0.75e-3
nsel,r,loc,y,-2.7e-3
esel,s,type,,2
sf,all,fsi,1i
nsel,s,loc,x,-0.75e-3
nsel,r,loc,y,-2.7e-3,2.7e-3
esel,s,type,,2
sf,all,fsi,1
nsel,s,loc,x,0.75e-3
nsel,r,loc,y,-2.7e-3,2.7e-3
esel,s,type,,2
sf,all,fsi,1
Lsel,s,,,6,7,1
Lsel,a,,,11,12,1
Lsel,a,,,15
NSLL,S,1
esel,s,type,,2
sf,all,fsi,1
esel,all
nsel,all

!边界条件
nsel,s,loc,x,-0.75e-3,0.75e-3
nsel,r,loc,y,-2.7e-3
d,all,ux,0,,,uy
alls
fini
/solu
/psf,pres,norm,2,0,1
*do,t,1,36   !（移动36次）
outres,all,all
nsel,s,loc,x,-0.15e-3
nsel,r,loc,y,-0.9e-3,(t-9)*(1e-4)（对其中某一个位置加力）
sf,all,pres,10
alls
solve
sfdele,all,pres  !删除上述施加的面载荷
*enddo
alls
fini

biomechanics

sf,all,fsi,1i
nsel,r,loc,y,-0.9e-3,(t-9)*(1e-4)（对其中某一个位置加力）
这两处命令流有误，其次没定义fluid29的实常数reference pressure
再检查声固耦合

lchylinchy

谢谢啊，那应该怎样改啊?

biomechanics

sf,all,fsi,1i(改为1）
nsel,r,loc,y,-0.9e-3,(t-9)*(1e-4) !（对其中某一个位置加力） 少了注释符号
cp2和cp3感觉没耦合完，你的单位也有问题，如果说长度是m的话，密度等的单位就不对
d,n_bot1,volt,0为什么这点的电压为零，根据是什么

lchylinchy

d,n_bot1,volt,0 是为了将这个节点设置为低电压点

tzhdacy

一般是流固耦合问题

phyee88

好复杂的问题啊

eager0626

真没搞懂啊