压电分析问题请教

menghanbo

请问，做压电分析时如何添加边界条件呢？
/clear
/title,Model of Cantilever
/PREP7
/GRAPHICS,FULL
!##############################  parameter  ##############################
tstep=2*1E-6                    ! 加工技术步长
ch=1E-6
*SET,t_Si , 8.0*ch            ! Thickness of Si  
*SET,t_SiO2 , 0.8*ch             ! Thickness of   SiO2
*SET,t_pt1 , 0.2*ch            ! Thickness of up-pt  
*SET,t_PZt , 0.75*ch           ! Thickness of pzt
*SET,t_pt2 , 0.15*ch            ! Thickness of down-pt
BTOL,10E-4
VOL=0
VOH=10

!##############################  Material property   #############################

!********************** 1---Si **********************
ET,1,SOLID45                                 ! 用于硅薄膜层的8节点立方单元类型具有(UX,UY,UZ)的自由度 (solid 45)
MP,EX,1,1.2e11
MP,PRXY,1,0.42
MP,DENS,1,2329
R,1,t_Si
!********************** 2---SiO2 **********************
ET,2,SOLID45
MP,EX,2,0.72e11
MP,PRXY,2,0.17
MP,DENS,2,2200
R,2,t_SiO2
!********************** 3---pt **********************
ET,3,SOLID45
MP,EX,3,1.7e11
MP,PRXY,3,0.39
MP,DENS,3,21450
R,3,t_pt1
!********************** 4---PZT **********************
ET,4,SOLID226,1001
MP,EX,4,0.75e11
MP,PRXY,4,0.31
MP,DENS,4,7500
R,4,t_pzt

MP,PERX,4,804.6                                                 ! Permittivity of piezoceramic  (压电陶瓷的介电常数）
MP,PERY,4,804.6
MP,PERZ,4,659.7
TB,PIEZ,4   
TBDATA,3,-5.2                                                   !（给当前数据表定义数据）
TBDATA,6,-5.2
TBDATA,9,15.1
TBDATA,14,12.7   
TBDATA,16,12.7
!********************** 5---pt **********************
ET,5,SOLID45
MP,EX,5,1.7e11
MP,PRXY,5,0.39
MP,DENS,5,21450
R,5,t_pt2

!********************** end **********************
!下面定义几何模型参数和建立几何模型

!********************** 1--Si **********************
NUMSTR,KP,1
NUMSTR,LINE,1
NUMSTR,AREA,1
k,1,100*ch,0,0
k,2,100*ch,-25*ch,0
k,3,200*ch,-25*ch,0
k,4,200*ch,-75*ch,0
k,5,100*ch,-75*ch,0
k,6,100*ch,-100*ch,0
k,7,50*ch,-100*ch,0
k,8,50*ch,-600*ch,0
ksymm,1,1,8
a,1,2,3,4,5,6,7,8,16,15,14,13,12,11,10,9
VOFFST,1,-t_Si
!********************** 2--SiO2 **********************
WPOFFS,0,0,t_Si               
AGEN,2,2,,,,,,,,
VOFFST,19,-t_SiO2
!********************** 3--pt1 **********************
WPOFFS,0,0,t_SiO2
AGEN,2,20,,,,,,,,
VOFFST,37,-t_pt1
!********************** 4--PZT **********************
WPOFFS,0,0,t_pt1
AGEN,2,38,,,,,,,,
VOFFST,55,-t_PZt
!********************** 5--pt2 **********************
WPOFFS,0,0,t_PZt
AGEN,2,56,,,,,,,,
VOFFST,73,-t_pt2
WPAVE,0,0,0                                                    !工作平面回到原始位置
WPSTYLE,,,,,,,,0                                               !不显示工作平面   
vglue,all
NUMCMP,all
/VIEW,1,1,1,1                                                  ! 变换视角
/AUTO,1                                                        ! 自动适应窗口大小  
/TRIAD,OFF                                                     ! 关闭坐标系显示
VPLOT   
!**************************************** 划分网格 ********************************************************
/COM,                                                          ! 对于四种材料网格剖分  
ESIZE,5*tstep                                                  ! 设置单元尺寸  

/COM,                                                          ! 压电层与硅膜层耦合网格剖分  
VSEL,S,LOC,Z,t_Si+t_SiO2+t_pt1 ,t_Si+t_SiO2+t_pt1+t_pzt        ! 选择压电层体
CM,v_pzt,VOLU                                                  ! 定义压电层为一个组，名字为v_pzt  
VATT,4,4,4                                                      ! 定义压电层材料类型及单元类型  

VSEL,S,LOC,Z,0,t_Si                                            ! 选择硅膜层体  
CM,v_si,VOLU                                                   ! 定义硅膜层为一个组，名字为v_si
VATT,1,1,1                                                      ! （定义体的属性）定义压电层材料类型及单元类型 （VATTT, MAT, REAL, TYPE, ESYS）

VSEL,S,LOC,Z,t_Si,t_Si+t_SiO2                                ! 选择二氧化硅膜层体
CM,v_SiO2,VOLU                                                 ! 定义二氧化硅膜层为一个组，名字为v_SiO2
VATT,2,2,2                                                      ! （定义体的属性）定义二氧化硅材料类型及单元类型

VSEL,S,LOC,Z,t_Si+t_SiO2,t_Si+t_SiO2+t_pt1                                ! 选择下电极层体
CM,v_pt1,VOLU                                                             ! 定义下电极层为一个组，名字为v_SiO2
VATT,3,3,3                                                                ! （定义体的属性）定义pt1材料类型及单元类型

VSEL,S,LOC,Z,t_Si+t_SiO2+t_pt1+t_pzt,t_Si+t_SiO2+t_pt1+t_pzt+t_pt2        ! 选择上电极层体
CM,v_pt2,VOLU                                                             ! 定义上电极层为一个组，名字为v_pt2  
VATT,5,5,5                                                                ! 定义pt2层材料类型及单元类型                                               

CMSEL,all,volum                                           ! 添加选中压电层 （CMSEL，Type, Name, Entity 选择元素和集合的子集，a表示添加）

VSWEEP,ALL                                                     ! （体扫掠生成单元）对选中体划分网格（VSWEEP, VNUM, SRCA, TRGA, LSMO,含义：从邻近的面穿过体通过扫掠的方法填充没划分网格的体）

/COM,                                                          ! 下电极网格划分及节点组件定义  
ASEL,S,,,19                                                    ! 选择下电极层  
nsla,s,1
CM,bot_elec,NODE                                               ! 定义下电极节点组件  
CP,1,VOLT,ALL                                                  ! 耦合电压自由度在下电极面  
*GET,nel1,NODE,0,NUM,MIN                                       ! 获得所选节点的最低节点序号赋值给nel1变量

ASEL,S,,,22                                                    ! 选择上电极层  
nsla,s,1  
CM,top_elec,NODE                                               ! 定义下电极节点组件
CP,2,VOLT,ALL                                                  ! 耦合电压自由度在上电极面 （2为耦合组编号）
*GET,nel2,NODE,0,NUM,MIN                                        ! 获得所选节点的最低节点序号赋值给nel2变量
ALLSEL                                                         ! 选择所有实体  
FINISH                                                         ! 退出当前处理器  
                                                   
!**************************************** 定义边界条件  *************************************************
!**************************************** boundary condition  *************************************************
ASEL,S,LOC,Y,O                                                            !面 0

ASEL,S,LOC,X,200E-6                                                       !面1                              
ASEL,R,LOC,Y,-25E-6,O

ASEL,S,LOC,X,-200E-6                                                      !对称面1
ASEL,R,LOC,Y,-25E-6,O

ASEL,S,LOC,Y,-25E-6                                                       !面2
ASEL,R,LOC,X,100E-6,200E-6

ASEL,S,LOC,Y,-25E-6                                                       !对称面2
ASEL,R,LOC,X,-100E-6,-200E-6

ASEL,S,LOC,X,200E-6                                                       !面3
ASEL,R,LOC,Y,-75E-6,-25E-6

ASEL,S,LOC,X,-200E-6                                                      !对称面3
ASEL,R,LOC,Y,-75E-6,-25E-6

ASEL,S,LOC,Y,-75E-6                                                       !面4
ASEL,R,LOC,X,-200E-6,-100E-6

ASEL,S,LOC,Y,-75E-6                                                       !对称面4
ASEL,R,LOC,X,100E-6,200E-6

ASEL,S,LOC,X,100E-6                                                       !面5
ASEL,R,LOC,Y,-100E-6,-75E-6

ASEL,S,LOC,X,-100E-6                                                      !对称面5
ASEL,R,LOC,Y,-100E-6,-75E-6

ASEL,S,LOC,Y,-100E-6                                                      !面6
ASEL,R,LOC,X,50E-6,100E-6

ASEL,S,LOC,Y,-100E-6                                                      !对称面6
ASEL,R,LOC,X,-100E-6,-50E-6

DA,AREA,ALL                                                                 
DALIST,Area                                                                !列表显示面上的DOF约束                                                        

ALLSEL,ALL                                                              
FINISH                                                                  
!*******************************************  求解  ***********************************************
/solu
antype,static
ALLS

d,bot_elec,volt,VOL
d,top_elec,volt,VOH

time,1
nsubst,5                                                                       !子步数为5
outres,all,all                                                                 !输出数据库选项
kbc,0                                                                          !设置载荷变化方式
nsel,all                                                                        !选择节点
esel,all                                                                        !选择元素
save
solve
finish

/POST1                                                                          ! 进入通用后处理器  
/TITLE,Electric potential distribution in model for static voltage of 10V       !定义标题（Z方向位移）
PLNSOL,VOLT                                                                     ! 显示电压分布云图  
/WAIT,10                                                                        ! 延时10S  
/TITLE,Vertical displacements in the piezo micropump.                           !定义标题（Z方向位移）
PLNSOL,U,Z                                                                      ! 显示Z方向位移云图  
FINISH                                                                          ! 退出当前处理器


这是我的命令流，哪位能帮我看看哪个地方出现问题了呢？小弟感激不尽．期待您的解答．