请教一下大家关于压电陶瓷电压加载的问题（附命令流和模型）

wuming1112

现在我想在三个压电陶瓷中加载两个100sinwt和100coswt的交变电压。(相位相差90度)
压电陶瓷是沿z向极化的。其中第一和第三压电陶瓷极化方向相反，加载100coswt。第二个极化方向与第一个压电陶瓷相同，但要加载100sinwt。
因为刚刚接触压电分析。所以又很多问题不太懂。麻烦大家指点
1、极化方向方面：
压电陶瓷我选用pzt-8，参数按照书上查得。陶瓷伸长方向沿总体笛卡尔坐标系的z轴。这样的话我的电压是不是应该加在两个侧面？
2、压电陶瓷变形引起底板变形。所以我体之间用了vglue对吗？
3、就是电压的加载形式。我直接选择面。这样对吗？
4、就是我现在最不知道的。。
   相位相差90度的电压应该怎么加载呢？
附上和压电陶瓷有关的命令流（自己编的不太规范，大家别见笑）
以及模型

命令流：
et,2,solid5,3
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
MP,DENS,3,7500       !+            
MP,PERX,3,7.124e-09   
MP,PERY,3,7.124e-09   
MP,PERZ,3,5.841e-09   
   
TB,ANEL,3                    
tbdata,1,1.32e11,7.1e10,7.3e10,0,0,0   
tbdata,7,1.32e11,7.3e10,0,0,0        
tbdata,12,11.5E10,0,0,0                 
tbdata,16,3.06E10,0,0                 
tbdata,19,2.56E10,0               
tbdata,21,2.56E10                 
TB,PIEZ,3                 
tbdata,3,-4.1   
tbdata,6,-4.1   
tbdata,9,14.1   
tbdata,14,10.5   
tbdata,16,10.5           
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
MP,DENS,4,7500            !-      
MP,PERX,4,7.124e-09  
MP,PERY,4,7.124e-09   
MP,PERZ,4,5.841e-09
         
TB,ANEL,4        
tbdata,1,1.32e11,7.1e10,7.3e10,0,0,0   
tbdata,7,1.32e11,7.3e10,0,0,0        
tbdata,12,11.5E10,0,0,0                 
tbdata,16,3.06E10,0,0                 
tbdata,19,2.56E10,0               
tbdata,21,2.56E10  
TB,PIEZ,4           
tbdata,3,4.1           
tbdata,6,4.1   
tbdata,9,-14.1         
tbdata,14,-10.5      
tbdata,16,-10.5      
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
BLOCK,0,h1,-b1,b1,-l1,l1,
BLOCK,-h2,0,-b2,b2,-l2,l2,
VGEN,2,2, , , , ,32e-3, ,1
VGEN,2,2, , , , ,-32e-3, ,1
vsel,all
vglue,all
NUMCMP,all
TYPE,   2   
MAT,       2
REAL,   
ESYS,       0   
SECNUM,
ESIZE,0.003,0,  
MSHAPE,0,3D
MSHKEY,1
VSEL, , , ,1,2
VMESH,all

TYPE,   2   
MAT,       3
REAL,   
ESYS,       0   
SECNUM,
ESIZE,0.004,0,  
MSHAPE,0,3D
MSHKEY,1
VSEL, , , , 3
VMESH,all

TYPE,   1   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,
ESIZE,0.015,0,   
SMRT,8  
MSHAPE,1,3D
MSHKEY,0
VSEL, , , ,       4   
VMESH,all   
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ASEL,S,,,18      
NSLA,S,1
cp,18,volt,all
allsel,all
ASEL,S,,,19      
NSLA,S,1
cp,19,volt,all
allsel,all
ASEL,S,,,6      
NSLA,S,1
cp,6,volt,all
allsel,all
ASEL,S,,,7      
NSLA,S,1
cp,7,volt,all
allsel,all
ASEL,S,,,12      
NSLA,S,1
cp,12,volt,all
allsel,all
ASEL,S,,,13     
NSLA,S,1
cp,13,volt,all
allsel,all
da,18,volt,VOL
da,19,volt,VOL
da,6,volt,VOL
da,7,volt,VOL
da,12,volt,VOL
da,13,volt,VOL
fini
save
!~~~~~~~~~~~~~~~~~~~~~~~~模态分析
/SOL
ANTYPE,2
MSAVE,0
MODOPT,LANB,20  
EQSLV,SPAR  
MXPAND,20, , ,0
LUMPM,0
PSTRES,0
MODOPT,LANB,20,0,0, ,OFF
SOLVE
FINISH

xuruikl

我没做过压电
但我觉得你下面那个大长方体的网格有必要改进一些

wuming1112

是指不能自由划分网格吗？

阿政狼