在复合材料中间挖去一尖锐孔后填入别的材料，整体进行...
我在一块复合材料中间挖去一个孔（挖去的孔有尖锐部分），然后在空中填入别的材料，最后对整体进行拉伸，对拉伸后的云图查看发现在空的尖锐部分没有出现应力集中现象，麻烦哪位大神能帮我看一下吗？（该复合材料有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!绘制等效应力等值线图
页:
[1]