！！！求教压电与层合板耦合问题 ，盼高手指点！！！

daifengyiqiu

压电与层合板通过节点耦合
这是我的命令流

/com,learning
/title,learning pzt4 and composite plate

/prep7
!#定义参数
th=0.02       !层厚(m)
l=1           !边长(m)
E1=132.4E9    !复合材料特性单位(N/M2)
E2=10.8E9
E3=10.8E9
G12=5.6E9
G13=5.6E9
G23=3.6E9
V12=0.24
V13=0.24
V23=0.49

!选择单元：压电用solid226(8节点--4DOF/node:UX, UY, UZ,VOLT)
!               层复合材料用solid191(8节点--3DOF/node:UX, UY, UZ)
et,1,solid226,1001   !压电单元1
et,2,solid191,,,     !层复合单元2

! solid191的实常数   
r,1,4,,,,,
rmore,
rmore,2,0,th,2,90,th    ! 2为材料号(层复合材料)
rmore,2,90,th,2,0,th

!输入压电应力矩阵(e) 单位C/m2
TB,PIEZ,1               !1为材料号(压电)                        
TBDATA,16,10.5                     ! E61 PIEZOELECTRIC CONSTANT
TBDATA,14,10.5                     ! E52 PIEZOELECTRIC CONSTANT
TBDATA,3,-4.1                      ! E13 PIEZOELECTRIC CONSTANT
TBDATA,6,-4.1                      ! E23 PIEZOELECTRIC CONSTANT
TBDATA,9,14.1                      ! E33 PIEZOELECTRIC CONSTANT

!输入压电刚度矩阵(c) 单位N/m2
TB,ANEL,1                        
TBDATA,1,13.2E-10,7.1E-10,7.3E-10     
TBDATA,7,13.2E-10,7.3E-10
TBDATA,12,11.5E-10
TBDATA,16,3.0E-10
TBDATA,19,2.6E-10
TBDATA,21,2.6E-10

!输入常应变介电矩阵[εS]
MP,PERX,1,804.6                    ! PERMITTIVITY (X AND Y DIRECTION)
MP,PERZ,1,659.7                    ! PERMITTIVITY (Z DIRECTION)

!输入层复合材料特性
mp,ex,2,e1
mp,ey,2,e2
mp,ez,2,e3
mp,prxy,2,v12
mp,pryz,2,v23
mp,prxz,2,v13
mp,gxy,2,g12
mp,gyz,2,g23
mp,gxz,2,g13

!建立几何模型
BLOCK, 0, l, 0, l, 0, 4*th             !层复合板(体)
BLOCK,0,l,0,l,4*th,5*th                   !压电(体)

!分配单元属性
vsel,all
vsel,s,loc,z,0,4*th
vatt,2,1,2,           !层复合材料
vsel,all
vsel,s,loc,z,4*th,5*th
vatt,1,,1,            !压电  

!划分单元
vsel,all
esize,th
vmesh,all

!耦合节点
nsel,all
nsel,s,loc,z,4*th
cp,1,ux,all
cp,1,uy,all
cp,1,uz,all              

!求解
/solu     
antype,static
lsel,all
!位移约束
lsel,s,loc,z,0
nsel,r,loc,z,0
d,all,ux,0
d,all,uy,0
d,all,uz,0

!加电压
asel,all
asel,s,loc,z,5*th
nsel,r,loc,z,5*th
d,all,volt,-100
solve

**********************************
出现如下的错误不知什么原因。 望赐教
*** WARNING ***                         CP=      6.950   TIME= 09:58:19
Node 2 on element 1 is unselected.                                      

*** WARNING ***                         CP=      6.950   TIME= 09:58:19
Node 303 on element 1 is unselected.                                    

*** WARNING ***                         CP=      6.950   TIME= 09:58:19
Node 451 on element 1 is unselected.                                    

*** WARNING ***                         CP=      6.950   TIME= 09:58:19
Node 100 on element 1 is unselected.                                    

*** WARNING ***                         CP=      6.950   TIME= 09:58:19
Node 302 on element 1 is unselected.                                    

*** WARNING ***   SUPPRESSED MESSAGE    CP=      6.950   TIME= 09:58:19
Node 401 on element 1 is unselected.                                    

*** WARNING ***   SUPPRESSED MESSAGE    CP=      6.950   TIME= 09:58:19
Node 450 on element 1 is unselected.                                    

*** WARNING ***   SUPPRESSED MESSAGE    CP=      6.950   TIME= 09:58:19
Node 101 on element 1 is unselected.                                    

*** WARNING ***   SUPPRESSED MESSAGE    CP=      6.950   TIME= 09:58:19
Node 15403 on element 1 is unselected.
。。。很多这样的错误！！！

！！！！望赐教  。我的邮箱 daifeng0301243@126.com