求助：圆柱体与长方体相交的网格划分！

chq8486

小弟是新手，最近接触到一个模型（长方体上面有一个圆柱体），见下图，在划分网格的时候一直不成功请各位高手指教，谢谢！我的命令流我也附上了，请高手帮我修改下，谢谢！

/FILNAM,cylinder
/TITLE,cylinder&block
/UNITS,si,            
tstep=1e-6
l=14e-6
w=14e-6
t=3e-6
r=5e-6
tc=2e-6
tt=1e-6
tb=1e-6
/PREP7
ET,1,5,3
ET,2,45
ET,3,41
MP,DENS,1,7500
MP,PERX,1,804.6
MP,PERY,1,804.6
MP,PERZ,1,659.7
TB,PIEZ,1
TBDATA,3,-4.1   
TBDATA,6,-4.1
TBDATA,9,14.1
TBDATA,14,10.5
TBDATA,16,10.5
TB,ANEL,1
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
MP,EX,2,120e9            
MP,NUXY,2,0.42         
MP,DENS,2,2329
MP,EX,3,95e9               
MP,NUXY,3,0.33            
MP,DENS,3,21450           
R,1,tt1                  
R,2,tb2                 
FINISH   
/prep7
BLC5,0,0,l,w,t
CYL4,0,0,0,,r,,t+tC
BTOL,2e-007
VOVLAP,ALL
/VIEW,1,1,1,1     
/AUTO,1                        
/TRIAD,OFF                       
VPLOT                          
!------------以下分网就有问题了，
ESIZE,tstep
/COM,    !Mesh block  
VSEL,S,LOC,Z,0,t
TYPE,2
MAT,2
CM,blcs,VOLU
VMESH,ALL
/COM,    !Mesh cylinder
VSEL,S,LOC,Z,t,t+tc
TYPE,1
MAT,1
CM,cylp,VOLU
*GET,Piel,VOLU,1,VSEL
VMESH,ALL
/COM,    !Mesh shell-1
TYPE,3
MAT,3
ASEL,S,LOC,Z,t
ASEL,R,LOC,X,0,l
ASEL,R,LOC,Y,0,w  
REAL,1
AMESH,ALL
/COM,    !Mesh shell-2
TYPE,3
MAT,3
ASEL,S,LOC,Z,t+tc
ASEL,R,LOC,X,0,l
ASEL,R,LOC,Y,0,w  
REAL,2
AMESH,ALL
ALLSEL   
FINISH

kuangulugu

切割，扫略画

jiandan32

或者用接触来做