在复合材料中间挖去一尖锐孔后填入别的材料，整体进行...

神诸葛佳

我在一块复合材料中间挖去一个孔（挖去的孔有尖锐部分），然后在空中填入别的材料，最后对整体进行拉伸，对拉伸后的云图查看发现在空的尖锐部分没有出现应力集中现象，麻烦哪位大神能帮我看一下吗？（该复合材料有8层，每层厚度为0.125mm，铺层方向都是90度）
并且拉伸后，复合材料整体弯曲了，按理说不应该啊，因为铺层都是90度，不应该弯曲才对啊，有哪位大神能够够为我解答吗？是我设置的材料类型有问题？还是网格划分有问题？还是有别的原因？
下面是我的命令流以及应力云图
/PREP7  
/VIEW,1,1,1,1
*SET,H,50   !设置复合材料长度，
*SET,P,20   !设置复合材料长度方向网格划分段数为P

!!!!!!!!!!!!设置材料属性!!!!!!!!!!!!
ET,1,SOLID186   
ET,2,SOLID45   
ET,3,SOLID185   
ET,4,SOLSH190   !复材基体

KEYOPT,4,2,0
KEYOPT,4,6,0
KEYOPT,4,8,1    !设置输出参数，输出各层数据

MP,EX,1,72E3    !嵌入的第1种材料弹性模量
MP,PRXY,1,0.17  !泊松比  
MP,EX,2,1.4E3   !嵌入的第2种材料弹性模量
MP,PRXY,2,0.48  
MP,EX,3,3.5E3   !嵌入的第3种材料弹性模量
MP,PRXY,3,0.35  

MP,EX,4,181E3   !复合材料基体弹性模量
MP,EY,4,10.3E3  
MP,EZ,4,10.3E3  
MP,PRXY,4,0.28   
MP,PRYZ,4,0.3  
MP,PRXZ,4,0.28   
MP,GXY,4,7.17E3  
MP,GYZ,4,3.78E3  
MP,GXZ,4,7.17E3   

/RGB,INDEX,100,100,100, 0   !选择白色背景
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15  
/REPLOT


!!!设置复合材料每层信息
SECT,1,SHELL,,  
SECDATA, 0.025,4,90.0,3  !复合材料共八层，每层厚度都是0.125mm，都为90度，
SECDATA, 0.025,4,90.0,3  
SECDATA, 0.025,4,90.0,3  
SECDATA, 0.025,4,90.0,3  
SECDATA, 0.025,4,90.0,3  
SECDATA, 0.025,4,90.0,3  
SECDATA, 0.025,4,90.0,3  
SECDATA, 0.025,4,90.0,3  

BLOCK,0,H,-6.25,6.25,0,1  !复合材料尺寸

!!!!在复合材料中挖孔，并填入别的材料!!!!!!
/VIEW,1,1
WPROT,0,0,90              
CYL4,0,0,0.625, , , ,H  
WPOFF,-1,0,0               
CYL4,0,0,0.625, , , ,H   
VINV,2,3                  
WPOFF,0.5,0,0           
CYL4,0,0,0.125, , , ,H      
CYL4,0,0,0.0625, , , ,H     
VPTN,ALL            

!为了划分网格方便，切分嵌入内部的材料
VSEL,S, , ,3
VSEL,A, , ,7
wprot,0,90,0
VSBW,ALL
wprot,0,0,90
VSBW,ALL

ALLSEL,ALL       !选择所有实体
VPLOT            !显示所有实体
WPCSYS,-1
VGLUE,ALL              !对所有实体进行粘合   
NUMMRG,ALL             !合并相同位置上重复定义的关键点、线、面等
NUMCMP,ALL             !压缩所有实体编号

!!!!划分网格
!划分嵌入的第1种材料
/VIEW,1,1
VSEL,S, , ,1     !选择体
VSEL,A, , ,4
VSEL,A, , ,6
VSEL,A, , ,7
MAT,1            !选择材料属性
TYPE,1           !选择单元类型
MSHAPE,0,3D      
MSHKEY,1         

ESIZE,0.02
LSEL,S,LENGTH,,H   !每条侧线网格划分数为P
LESIZE,ALL,,,P
VSWEEP,ALL


!划分嵌入的第2种材料
VSEL,S, , ,8
VSEL,A, , ,9
VSEL,A, , ,10
VSEL,A, , ,11
MAT,2   
TYPE,2  

ESIZE,0.02
VSWEEP,ALL


!划分嵌入的第3种材料
VSEL,S, , ,2
VSEL,A, , ,3
MAT,3   
TYPE,3

ESIZE,0.02
VSWEEP,ALL

!划分最外面的材料
VSEL,S, , ,5
MAT,4   
TYPE,4

ESIZE,0.1
VSWEEP,ALL

ALLSEL,ALL       !选择所有实体
WPCSYS,-1
NUMMRG,ALL  !合并相同位置上重复定义的关键点、线、面等
NUMCMP,ALL             !压缩所有实体编号
FINISH  

!求解
/SOL                !进入求解器
/VIEW,1,1,1,1   
EPLOT               !显示所有网格划分单元
ANTYPE,0            !制定求解类型（静态）

NSEL,S,LOC,X,0
D,ALL, ,0, , , ,ALL,

NSEL,S,LOC,X,H
SF,ALL,PRES,-10          !施加载荷

ALLSEL,ALL
SOLVE
/POST1  
PLDISP,2      !绘制变形图
PLNSOL,U,SUM  !绘制合位移等值线图
PLNSOL,S,EQV  !绘制等效应力等值线图