压电分析 问题

wangbaowei

下面是我分析的log 文件  在最后的top12 电极上加电压的时候 老是提示说什么0到0 无效

还有 压电片  加上电压之后没有能使结构变形  ， 而我在结构上加一个力的时候  整个结构都会变形   不知道 什么原因  请懂这方面的人帮忙解决  谢谢。
附件里也有我用的log 文件



WPSTYLE,,,,,,,,0
*SET,tstep , 10.0E-6            ! Technology step   
*SET,l_PZt , 50E-6              ! Length of PZt on plate
*SET,w_PZt , 50E-6              ! Width of PZt layer
*SET,t_PZt , 2.0E-6             ! Thickness of piezoelectric actuator ceramic  

*SET,lplate , 4*l_pzt             ! Length of plate 200E-6
*SET,wplate , 4*w_pzt            ! Width of plate  200E-6
*SET,tplate , 20.0E-6            ! Thickness of plate

*SET,l_Pt , l_pzt               ! Length of bottom Pt layer on plate
*SET,w_Pt , w_pzt               ! Width of bottom Pt layer  
*SET,t_Pt , 0.1E-6            ! Thickness of lower platinum electrode

/prep7  
ET,1,5,0
ET,2,45
ET,3,41
MP,DENS,1,7400  
MP,PERX,1,4.8E-9           ! Permittivity of piezoceramic   
MP,PERY,1,4.8E-9
MP,PERZ,1,4.8E-9
TB,PIEZ,1   
TBDATA,3,-1.86  
TBDATA,6,-1.86  
TBDATA,9,9.0
TBDATA,11,9.8   
TBDATA,13,9.8   
TB,ANEL,1   
TBDATA,1,13.5E10,6.74E10,6.81E10
TBDATA,7,13.5E10,6.81E10
TBDE,ANEL,1,,,  
TB,ANEL,1,1,21,0
TBTEMP,0
TBDATA,,1.35E+011,6.74E+010,6.81E+010,0,0,0
TBDATA,,1.35E+011,6.81E+010,0,0,0,11.3E10   
TBDATA,,0,0,0,2.22E10,0,0   
TBDATA,,2.22E10,0,3.355E10,,,   
TBDE,ANEL,1,,,  
TB,ANEL,1,1,21,0
TBTEMP,0
TBDATA,,1.35E+011,6.74E+010,6.81E+010,0,0,0
TBDATA,,1.35E+011,6.81E+010,0,0,0,1.13E+011
TBDATA,,0,0,0,2.22E+010,0,0
TBDATA,,2.22E+010,0,3.355E+010,,,   
MP,MURX,1,0
MP,KXX,1,0  
MP,EX,2,120e9   
MP,NUXY,2,0.42  
MP,DENS,2,2329  
MP,EX,3,95e9               ! Young's modulus of Pt  
MP,NUXY,3,0.33             ! Poisson's ration of of Pt  
MP,DENS,3,21450
R,1,t_Pt   
R,2,t_Pt  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,RSVX,3,,1.06e-7  
FINISH  
/PREP7

BLOCK,0,l_pzt,0,w_pzt,0,tplate,  ! Model structure volume 1-16
FLST,3,1,6,ORDE,1   
FITEM,3,1   
VGEN,4,P51X, , ,w_pzt, , , ,0   
FLST,3,4,6,ORDE,2   
FITEM,3,1   
FITEM,3,-4  
VGEN,4,P51X, , , ,l_pzt, , ,0   

BLOCK,0,l_pzt,l_pzt,2*l_pzt,tplate,tplate+t_pzt,   ! Modeling piezoelectric layers volume 17-28
BLOCK,w_pzt,2*w_pzt,0,l_pzt,tplate,tplate+t_pzt,
/PNUM,KP,0  
/PNUM,LINE,0
/PNUM,AREA,0
/PNUM,VOLU,1
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0   
/PNUM,ELEM,0
/REPLOT

FLST,3,1,6,ORDE,1           
FITEM,3,17  
VGEN,2,P51X, , , ,l_pzt, , ,0   
FLST,3,2,6,ORDE,2   
FITEM,3,17  
FITEM,3,19  
VGEN,4,P51X, , ,w_pzt, , , ,0   
FLST,3,1,6,ORDE,1   
FITEM,3,18  
VGEN,2,P51X, , ,w_pzt, , , ,0   
FLST,3,2,6,ORDE,2   
FITEM,3,18  
FITEM,3,26  
VGEN,2,P51X, , , ,3*l_pzt, , ,0

FLST,2,16,6,ORDE,2   ! Glue volume 1-16  
FITEM,2,1   
FITEM,2,-16
VGLUE,P51X  
FLST,2,2,6,ORDE,2   ! Glue volume
FITEM,2,18  
FITEM,2,30  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,26  
FITEM,2,31  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,17  
FITEM,2,29  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,19  
FITEM,2,33  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,20  
FITEM,2,35  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,21  
FITEM,2,38  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,22  
FITEM,2,36  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,23  
FITEM,2,39  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,24  
FITEM,2,37  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,25  
FITEM,2,40  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,27  
FITEM,2,41  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,28  
FITEM,2,42  
VGLUE,P51X  

WPSTYLE,,,,,,,,1           ! Modeling electrode
wpoff,0,w_pt,tplate
CSYS,4  
RECTNG,0,l_pt,0,w_pt,     ! Bottom electrode plane 58-69
RECTNG,l_pt,2*l_pt,-w_pt,0, ! Bottom electrode plane 59


FLST,3,1,5,ORDE,1   
FITEM,3,58  
AGEN,2,P51X, , , ,l_pt, , ,0   
FLST,3,2,5,ORDE,2   
FITEM,3,58  
FITEM,3,60  
AGEN,4,P51X, , ,w_pt, , , ,0
FLST,3,1,5,ORDE,1   
FITEM,3,59  
AGEN,2,P51X, , ,w_pt, , , ,0
FLST,3,2,5,ORDE,2   
FITEM,3,59  
FITEM,3,67  
AGEN,2,P51X, , , ,3*w_pt, , ,0  

FLST,3,12,5,ORDE,2     ! Top electrode plane 70-82
FITEM,3,58  
FITEM,3,-69
AGEN,2,P51X, , , , ,t_pzt, ,0   

CSYS,0  
WPAVE,0,0,0
CSYS,0  

BLOCK,0,l_pzt,0,-tplate,tplate,-2*l_pzt,
BLOCK,lplate,lplate+tplate,0,l_pzt,tplate,-2*l_pzt,
FLST,3,1,6,ORDE,1   
FITEM,3,14  
VGEN,4,P51X, , ,l_pzt, , , ,0   
FLST,3,1,6,ORDE,1   
FITEM,3,15  
VGEN,4,P51X, , , ,l_pzt, , ,0   
BLOCK,lplate,lplate+tplate,0,-tplate,tplate,-2*l_pzt,   
FLST,3,5,6,ORDE,5   
FITEM,3,14  
FITEM,3,16  
FITEM,3,29  
FITEM,3,-30
FITEM,3,36  
VGEN,2,P51X, , , ,lplate+tplate, , ,0   
FLST,3,6,6,ORDE,6   
FITEM,3,15  
FITEM,3,31  
FITEM,3,33  
FITEM,3,35  
FITEM,3,-36
FITEM,3,41  
VGEN,2,P51X, , ,-lplate-tplate, , , ,0  

FLST,2,20,6,ORDE,9  
FITEM,2,14  
FITEM,2,-16
FITEM,2,29  
FITEM,2,-31
FITEM,2,33  
FITEM,2,35  
FITEM,2,-42
FITEM,2,44  
FITEM,2,-48
VGLUE,P51X  

FLST,2,2,6,ORDE,2   
FITEM,2,1   
FITEM,2,14  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,2   
FITEM,2,50  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,3   
FITEM,2,53  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,32  
FITEM,2,54  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,15  
FITEM,2,32  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,10  
FITEM,2,52  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,11  
FITEM,2,55  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,43  
FITEM,2,56  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,43  
FITEM,2,61  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,13  
FITEM,2,60  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,12  
FITEM,2,59  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,34  
FITEM,2,37  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,34  
FITEM,2,65  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,5   
FITEM,2,64  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,4   
FITEM,2,63  
VGLUE,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,1   
FITEM,2,62  
VGLUE,P51X  

FLST,5,20,6,ORDE,13    ! Meshing structure
FITEM,5,14  
FITEM,5,-16
FITEM,5,29  
FITEM,5,-30
FITEM,5,32  
FITEM,5,-33
FITEM,5,35  
FITEM,5,-42
FITEM,5,44  
FITEM,5,49  
FITEM,5,51  
FITEM,5,57  
FITEM,5,-58
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y   
CMSEL,S,_Y1
VATT,       2,       1,   2,       0
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  

ESIZE,tstep,0,  
FLST,5,20,6,ORDE,13
FITEM,5,14  
FITEM,5,-16
FITEM,5,29  
FITEM,5,-30
FITEM,5,32  
FITEM,5,-33
FITEM,5,35  
FITEM,5,-42
FITEM,5,44  
FITEM,5,49  
FITEM,5,51  
FITEM,5,57  
FITEM,5,-58
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VSWEEP,_Y1  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  


FLST,5,16,6,ORDE,5  
FITEM,5,1   
FITEM,5,-13
FITEM,5,31  
FITEM,5,34  
FITEM,5,43  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
CMSEL,S,_Y1
VATT,       2,       1,   2,       0
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  

ESIZE,tstep,0,  
FLST,5,16,6,ORDE,5  
FITEM,5,1   
FITEM,5,-13
FITEM,5,31  
FITEM,5,34  
FITEM,5,43  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y   
VSWEEP,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  

FLST,5,12,6,ORDE,2    ! Mesh piezoelectric layers
FITEM,5,17  
FITEM,5,-28
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y   
CMSEL,S,_Y1
VATT,       1,       1,   1,       0
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  

ESIZE,tstep,0,
FLST,5,12,6,ORDE,2  
FITEM,5,17  
FITEM,5,-28
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y   
VSWEEP,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  


ASEL,S,LOC,Z,tplate  ! Mesh bottom electrode A58-69
APLOT   
FLST,5,12,5,ORDE,2  
FITEM,5,58  
FITEM,5,-69
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
CMSEL,S,_Y1
AATT,       3,       1,   3,       0,   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  

FLST,2,12,5,ORDE,2  
FITEM,2,58  
FITEM,2,-69
AESIZE,P51X,tstep,  
MSHAPE,1,3D
MSHKEY,0
MSHAPE,0,2D

FLST,5,12,5,ORDE,2  
FITEM,5,58  
FITEM,5,-69
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y   
AMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  

ALLSEL,ALL  
ASEL,S,LOC,Z,tplate+t_pzt   
APLOT   
FLST,5,12,5,ORDE,4  
FITEM,5,70  
FITEM,5,-75
FITEM,5,77  
FITEM,5,-82
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
CMSEL,S,_Y1
AATT,       3,       2,   3,       0,   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  

FLST,2,12,5,ORDE,4  
FITEM,2,70  
FITEM,2,-75
FITEM,2,77  
FITEM,2,-82
AESIZE,P51X,tstep,  
FLST,5,12,5,ORDE,4  
FITEM,5,70  
FITEM,5,-75
FITEM,5,77  
FITEM,5,-82
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y   
AMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  


ALLSEL,ALL  
ASEL,S,LOC,Z,tplate   
APLOT   
ASEL,S, , ,      58    ! Building bottom electrode 1
ALLSEL,BELOW,AREA   
cm,bot1,node
cp,1,volt,all   
*get,nel1,node,,num,min

ASEL,S, , ,      61   ! Building bottom electrode 2
ALLSEL,BELOW,AREA   
cm,bot2,node
cp,2,volt,all   
*get,nel2,node,,num,min

ASEL,S, , ,      63   ! Building bottom electrode 3
ALLSEL,BELOW,AREA   
cm,bot3,node
cp,3,volt,all   
*get,nel3,node,,num,min

ASEL,S, , ,      65   ! Building bottom electrode 4
ALLSEL,BELOW,AREA   
cm,bot4,node
cp,4,volt,all   
*get,nel4,node,,num,min

ASEL,S, , ,      60  ! Building bottom electrode 5
ALLSEL,BELOW,AREA   
cm,bot5,node
cp,5,volt,all   
*get,nel5,node,,num,min
ASEL,S, , ,      62   ! Building bottom electrode 6
ALLSEL,BELOW,AREA   
cm,bot6,node
cp,6,volt,all   
*get,nel6,node,,num,min
ASEL,S, , ,      64    ! Building bottom electrode 7
ALLSEL,BELOW,AREA   
cm,bot7,node
cp,7,volt,all   
*get,nel7,node,,num,min
ASEL,S, , ,      66   ! Building bottom electrode 8
ALLSEL,BELOW,AREA   
cm,bot8,node
cp,8,volt,all   
*get,nel8,node,,num,min
ASEL,S, , ,      59   ! Building bottom electrode 9
ALLSEL,BELOW,AREA   
cm,bot9,node
cp,9,volt,all   
*get,nel9,node,,num,min
ASEL,S, , ,      67   ! Building bottom electrode 10
ALLSEL,BELOW,AREA
cm,bot10,node   
cp,10,volt,all  
*get,nel10,node,,num,min
ASEL,S, , ,      68   ! Building bottom electrode 11
ALLSEL,BELOW,AREA   
cm,bot11,node   
cp,11,volt,all  
*get,nel11,node,,num,min
ASEL,S, , ,      69   ! Building bottom electrode 12
ALLSEL,BELOW,AREA   
cm,bot12,node   
cp,12,volt,all  
*get,nel12,node,,num,min

ALLSEL,ALL  
ASEL,S,LOC,Z,tplate+t_pzt   
APLOT                  
ASEL,S, , ,      70   ! Building top electrode 1
ALLSEL,BELOW,AREA   
cm,top1,node
cp,13,volt,all  
*get,nel13,node,,num,min
ASEL,S, , ,      73   ! Building top electrode 2
ALLSEL,BELOW,AREA   
cm,top2,node
cp,14,volt,all  
*get,nel14,node,,num,min
ASEL,S, , ,      75   ! Building top electrode 3
ALLSEL,BELOW,AREA
cm,top3,node
cp,15,volt,all  
*get,nel15,node,,num,min
ASEL,S, , ,      78   ! Building top electrode 4
ALLSEL,BELOW,AREA   
cm,top4,node
cp,16,volt,all  
*get,nel16,node,,num,min
ASEL,S, , ,      72   ! Building top electrode 5
ALLSEL,BELOW,AREA   
cm,top5,node
cp,17,volt,all  
*get,nel17,node,,num,min
ASEL,S, , ,      74   ! Building top electrode 6
ALLSEL,BELOW,AREA   
cm,top6,node
cp,18,volt,all  
*get,nel18,node,,num,min
ASEL,S, , ,      77   ! Building top electrode 7
ALLSEL,BELOW,AREA   
cm,top7,node
cp,19,volt,all  
*get,nel19,node,,num,min
ASEL,S, , ,      79   ! Building top electrode 8
ALLSEL,BELOW,AREA   
cm,top8,node
cp,20,volt,all  
*get,nel20,node,,num,min
ASEL,S, , ,      71   ! Building top electrode 9
ALLSEL,BELOW,AREA   
cm,top9,node
cp,21,volt,all  
*get,nel21,node,,num,min
ASEL,S, , ,      80   ! Building top electrode 10
ALLSEL,BELOW,AREA   
cm,top10,node   
cp,22,volt,all  
*get,nel22,node,,num,min
ASEL,S, , ,      81   ! Building top electrode 11
ALLSEL,BELOW,AREA   
cm,top11,node   
cp,23,volt,all  
*get,nel23,node,,num,min
ASEL,S, , ,      82   ! Building top electrode 12
ALLSEL,BELOW,AREA   
cm,top112,node  
cp,24,volt,all  
*get,nel24,node,,num,min

ALLSEL,ALL  
ASEL,S,LOC,Z,-2*l_pzt   
APLOT   
FINISH  
/SOL
FLST,2,20,5,ORDE,20
FITEM,2,83  
FITEM,2,89  
FITEM,2,235
FITEM,2,301
FITEM,2,306
FITEM,2,310
FITEM,2,315
FITEM,2,319
FITEM,2,324
FITEM,2,327
FITEM,2,331
FITEM,2,335
FITEM,2,340
FITEM,2,345
FITEM,2,349
FITEM,2,354
FITEM,2,358
FITEM,2,360
FITEM,2,364
FITEM,2,369  
/GO
DA,P51X,ALL,   ! Constraint all DOFs

ALLSEL,ALL  
/REPLOT

d,bot1,volt,0   
d,bot2,volt,0   
d,bot3,volt,0   
d,bot4,volt,0   
d,bot5,volt,0   
d,bot6,volt,0   
d,bot7,volt,0   
d,bot8,volt,0   
d,bot9,volt,0   
d,bot10,volt,0  
d,bot11,volt,0  
d,bot12,volt,0  

d,top1,volt,50  
d,top2,volt,100
d,top3,volt,100
d,top4,volt,50  
d,top5,volt,50  
d,top6,volt,100
d,top7,volt,100
d,top8,volt,50  
d,top9,volt,50  
d,top10,volt,50
d,top11,volt,50
d,top12,volt,50

!ANTYPE,0
!/STATUS,SOLU
!SOLVE   
!FINISH

wangbaowei

有没有人帮帮忙  谢谢  我这个很着急的

charlie_x

联系我吧.
jmgao@ujs.edu.cn
QQ:188874724

feasky

注意极化方向

wangbaowei

我这是沿Z向极化 的， 哪位明白的 帮忙解决一下。谢谢。

hy1966

cm,top112,node  
cp,24,volt,all  
*get,nel24,node,,num,min

没什么大错，好象这里定义的是top112，后面负载加的是top12，是空域，不允许

lizhanhui

谢谢！！
非常不错！！

tg3222138

谢谢！！
非常不错！！