『求助』压电材料分析，有问题问高手！！！1

做了一个悬臂梁的压电分析，具体命令流如下，但却出现了个怪问题，上面的压电材料无反映，高手指教了，是否是cp定义得 有问题？
finish
/clear
/prep7
! difinition de la geometrie

l1=20
l2=40
e=4
ep=1
l=180

k,1,0,-e/2
k,2,l1,-(e/2+ep)
k,3,l1+l2,-(e/2+ep)
k,4,l1,-e/2
k,5,l1+l2,-e/2
k,6,l,-e/2
k,7,l,e/2
k,8,l1,e/2
k,9,l1+l2,e/2
k,10,l1,e/2+ep
k,11,l1+l2,e/2+ep
k,12,0,e/2

l,1,4
l,4,5
l,2,3
l,4,2
l,3,5
l,5,6
l,6,7
l,7,9
l,8,9
l,8,12
l,8,10
l,10,11
l,11,9
l,1,12

et,1,plane183
keyopt,1,3,2
et,2,plane223,1001
keyopt,2,3,2

!formation de l'aire

al,1,2,6,7,8,9,10,14
al,2,3,4,5
al,9,11,12,13

!difinir materiaux

mp,ex,1,4.5e6
mp,PRXY,1,0.35

![K]

mp,perx,2,804.6
mp,pery,2,804.6
mp,perz,2,659.7

![E]

tb,piez,2
tbdata,3,-4.1
tbdata,6,-4.1
tbdata,9,-14.1
tbdata,11,10.5
tbdata,14,10.5

![c]

tb,anel,2
tbdata,1,13.2e10,7.1e10,7.3e10
tbdata,7,13.2e10,7.3e10
tbdata,12, 11.5e10
tbdata,16,3.0e10
tbdata,19,2.6e10
tbdata,21,2.6e10

!maillage

lsel,s,,,11
lsel,a,,,13
lsel,a,,,4
lsel,a,,,5
lesize,all,,,2

lsel,s,,,7
lsel,a,,,14
lesize,all,,,4

lsel,s,,,1
lsel,a,,,10
lesize,all,,,20

lsel,s,,,2
lsel,a,,,3
lsel,a,,,9
lsel,a,,,12
lesize,all,,,40

lsel,s,,,6
lsel,a,,,8
lesize,all,,,80

type,1
mat,1
amesh,1

type,2
mat,2
amesh,2
type,2
mat,2
amesh,3

! Define top electrode

lsel,s,,,12
nsll,s,1
cp,1,volt,all
allsel

lsel,s,,,9
nsll,s,1
cp,2,volt,all
allsel

! Define bottom electrode

lsel,s,,,2
nsll,s,1
cp,3,volt,all
allsel

lsel,s,,,3
nsll,s,1
cp,4,volt,all
allsel

!dinition de la condition limite

lsel,s,,,14
dl,14,,ux,,0
dl,14,,uy,,0

/solu

antype,static

!definiton de la force

FK,7,FY,-100

solve