ansys压电能量收集器

ssgg1026

ansys悬臂梁压电能量收集器命令流，然后我想在质量块右端加一个轴向力，进行模态分析，但是不管加多大的力，结果不发生变化，请问哪里操作不对吗？还会弹出这个对话框。

/prep7
/RGB,INDEX,100,100,100, 0
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15
/REPLOT
et,1,solid227,1001      
et,2,solid187                  
mp,DENS,1,7500        
emunit,epzro,8.854e-12        
mp,perx,1,1704.4                  
mp,pery,1,1704.4   
mp,perz,1,1433.6
tb,ANEL,1                     
tbdata,1,12.72E10,8.021E10,8.467E10   
tbdata,7,12.72E10,8.467E10      
tbdata,12,11.744E10            
tbdata,16,2.299E10      
tbdata,19,2.299E10         
tbdata,21,2.347E10           
tb,PIEZ,1                  
tbdata,3,-6.6228               
tbdata,6,-6.6228            
tbdata,9,23.2403              
tbdata,14,17.0345               
tbdata,16,17.0345               
mp,dens,2,2700                 
mp,EX,2,70e9           
MP,DENS,3,8900  
MP,EX,3,207E9
csys,0
L=85e-3
b=10e-3
hp=0.3e-3
hs=0.3e-3
Lt=4e-3   
block,0,L,0,b,hp,hs+hp
block,0,L,0,b,0,hp      
block,0,L,0,b,hp+hs,2*hp+hs  
block,L,L+Lt,0,b,hp+1/2*hs-1/2*Lt,hp+1/2*hs+1/2*Lt
k,50,L+Lt,b/2,(2*hp+hs)/2
vglue,all
VSEL,S, , ,       1   
ESIZE,0.005,0,   
VATT,       2, ,   2,       0  
Vmesh,all            
VSEL,S, , ,       5
vplot
VATT,       1, ,   1,       0
ESIZE,0.0025,0,            
MSHAPE,1,3D
MSHKEY,0
VMESH,all  
VSEL,S, , ,       6
local,11,,,,,,180           
VATT,       1, ,   1,      11
ESIZE,5e-3,0,              
MSHAPE,1,3D
MSHKEY,0  
VMESH,all  
VSEL,S, , ,      7
VATT,       3, ,   2,       0
VMESH,all  
csys,0
nsel,s,loc,z,0            
nsel,u,loc, x,L,L+Lt
cp,1,volt,all         
*get,n_bot,node,0,num,min         
nsel,all   
nsel,s,loc,z,2*hp+hs            
nsel,u,loc, x,L,L+Lt
cp,2,volt,all
*get,n_top,node,0,num,min
nsel,all
nsel,s,loc,z,hp                  
nsel,u,loc, x,L,L+Lt
cm,cm_1,node
nsel,s,loc,z,hp+hs
nsel,u,loc, x,L,L+Lt
cm,cm_2,node
cmsel,s,cm_1,node
cmsel,a,cm_2,node
cp,3,volt,all
nsel,all
/solu
antype,modal            
modopt,lanb,3
mxpand,3
nsel,s,loc,x,0
d,all,ux,0,,,,uy,uz