求助压电微悬臂梁的问题

dxdang200

我写的压电微悬臂梁的程序，但是在主梁边缘加载一个FZ方向1E-2N的力求解以后，第一次网格用LESIZE,ALL,e-5, , , ,1, , ,1 ，求解以后后处理里面没有显示出PLOT RESULT
出现，后面几次我加大了网格，求解里面是出现了，但是会提示 the degree of freedom solution is not available....


另外我求解的时候还有一个警告，谁帮我看看是那里出了问题。。。多谢了！偶是新新手哇


/prep7
/title,PIEZ CANTILEVER BEAM
/UNIT,SI    ! Unit setting (International Standard)
!/define PZT material properties
mp,dens,1,7500
tb,ANEL,1                      ! stiffness C, N/m^2
tbdata,1,12.6E10,7.95E10,8.41E10   ! c11,c12,c13
tbdata,7,12.6E10,8.41E10          ! c11,c33
tbdata,12,11.7E10                ! c33
tbdata,16,2.30E10              
tbdata,19,2.30E10              
tbdata,21,2.35E10              
tb,PIEZ,1                      ! PIEZ constant, C/m^2
tbdata,3,-6.5                    ! e31
tbdata,6,-6.5                    ! e31
tbdata,9,23.3                    ! e33
tbdata,11,17.0                   ! e15
tbdata,13,17.0                   ! e15
!/define main beam material properties
mp,dens,2,7800
mp,EX,2,209e9
mp,nuxy,2,0.3
!/define dimension of PZT
L=3e-4
W=3e-4
H =2e-5
!/define dimension of main beam
L1=8E-4
H1=3e-5
block,,L,0,W,0,H                 !bottom PZT
block,0,L1,0,W,H,H+H1              !MAIN BEAM
block,0,L,0,W,H+H1,H+H1+H         !top PZT
vglue,all                           !booleam connection
et,1,solid5,3                       !define the element of PIEZ
et,2,solid45                        !define the element of beam
asel,s,loc,Z,H
cm,CM_1,area
cmplot,CM_1
accat,CM_1
asel,s,loc,Z,H1+H
cm,CM_2,area
accat,CM_2
! meshing
LESIZE,ALL,e-5, , , ,1, , ,1
mat,1 $ type,1         !BOTTOM PZT
Vmesh,1
mat,2 $ type,2           !MAIN BEAM
Vmesh,5
mat,1 $ type,1          !TOP pzt
vmesh,4
Nsel,s,loc,Z,0                  !define top electrode
cp,1,volt,all
*get,n_top,node,0,num,min
Nsel,s,loc,Z,H+H1+H                 !define bottom electrode
cp,2,volt,all
*get,n_bot,node,0,num,min
Nsel,s,loc,Z,H              !couple the area
nsel,r,loc,x,0,L
cm,CM_3,node
Nsel,s,loc,Z,H+H1
Nsel,r,loc,x,0,L
cm,CM_4,node
cmsel,s,cm_3,node   
cmsel,a,cm_4,node
cp,3,volt,all
finish
/solu
ANTYPE,MODAL
MODOPT,LANB,3       ! EXTRACT 3 MODES
MXPAND,3
nsel,s,loc,x,0                              
d,all,ux,0,,,,uy,uz
d,n_top,volt,0              !short circuit the top electrode            
d,n_bot,volt,0              !short circuit the bottom electrode
nsel,all

hy1966

没什么问题，/post1中先读结果，显示！

dxdang200

多谢啦！问个弱智的问题，我是新学的！ 提示the requested epel data is not available是什么意思！

peakz

epel就是Elastic strains,这样你就明白了

dxdang200

像我这样的程序，要读Z方向的形变，与电压输出。。后处理程序应该怎么写呢！我后处理还一窍不通

peakz

path命令后，用map一下，然后选择Z方向的形变，与电压

peakz

path命令后，用map一下，然后选择Z方向的形变与电压