实体单元建模出现的问题? 20
感觉自己建的模型没有什么问题,模型由三部分组成(实体单元),封头+筒体+环向加强圈。因为是存在曲面,在建模时封头和筒体连接处存在小空隙。但是能正常出来结果,所以我认为模型基本没啥问题。加载应该也没啥问题,就是内压0.4MPa。最后结果出来最大应力在环向加强圈上,比较奇怪,觉得应该实在筒体上。所以又质疑自己的模型了。求大佬帮忙查看一下,看看问题出在哪里。
下面是两张结果图,第一张是全部实体应力图,第二张是单看加强圈上的应力图。发现最大应力在加强圈上
这是代码
!***************环境设置*********************
FINI
/CLEAR
/FILN,Test
/PREP7
!********参数设置********
!********筒体尺寸********
T=10 !全局厚度
P=0.4 !内压
Tank_Width=2396 !筒体宽度
Tank_Height=2396 !筒体高度
R1=Tank_Width/2
R2=Tank_Height/2
R3=R1+T
R4=R2+T
L=6020 !筒体总长度(含封头)
R0=3000 !建模辅助圆周半径
YJ_1=500 !筒体内圆角
YJ_2=YJ_1+T !筒体外圆角
!********封头尺寸********
H=500 !封头总深度
H1=H+T
!********工字钢********
Z1=-1080 !第一条加强圈的位置
W1=80 !底边缘长度
W2=80 !顶边缘长度
W3=140 !总深度
T1=5.5 !底边缘厚度
T2=5.5 !顶边缘厚度
T3=5.5 !腰厚度
!********定义材料和单元属性********
MP,EX,1,2E5
MP,NUXY,1,0.3
ET,1,SOLID185
!********建模********
!********建立出筒体轮廓********
K,1,0,0,0
K,2,R2-R0,0,0
K,3,0,R1-R0,0
K,4,R2,0,0
K,5,0,R1,0
K,6,R4,0,0
K,7,0,R3,0
CIRCLE,2,R0
CIRCLE,3,R0
CIRCLE,2,R0+T
CIRCLE,3,R0+T
LDELE,2,4,1,1
LDELE,6,8,1,1
LDELE,10,12,1,1
LDELE,14,16,1,1 !删除多余的线
L,4,6
L,5,7
LSBL,1,5,,,KEEP2
LSBL,5,6
LSBL,9,13,,,KEEP2
LSBL,13,6
LDELE,1,8,7,1
LFILLT,5,9,YJ_1
LFILLT,4,7,YJ_2 !分别给内外筒体构造圆角
LCOMB,1,9
LCOMB,1,5 !合并外筒体的线
LCOMB,4,6
LCOMB,4,7 !合并内筒体的线
ALLS,ALL
LGLUE,ALL
KDELE,2,3,1 !删除多余关键点
!********椭圆轮廓线********
K,,0,0,H
K,,0,0,H+T
CSYS,4 !激活工作平面
WPRO,0,90,0 !将工作平面绕X轴旋转90度
CSWPLA,11,1,H/R2,1 !在工作平面上建立一个局部坐标系-柱坐标系,Y/X=H/R2,Z/X=1
L,4,2 !构造出内椭圆线
CSWPLA,12,1,H1/R3,1 !在工作平面上建立一个局部坐标系-柱坐标系,Y/X=H1/R4,Z/X=1
L,16,3 !构造出外椭圆线
WPRO,0,-90,0 !将工作平面恢复,绕X轴旋转-90度
WPRO,0,0,90 !将工作平面绕Y轴旋转90度
WPRO,90,0,0 !将工作平面绕Z轴旋转90度
CSWPLA,13,1,H/R1,1 !在工作平面上建立一个局部坐标系-柱坐标系,Y/X=H/R2,Z/X=1
L,5,2 !构造出内椭圆线
CSWPLA,13,1,H1/R3,1 !在工作平面上建立一个局部坐标系-柱坐标系,Y/X=H1/R3,Z/X=1
L,21,3 !构造外椭圆线
WPRO,-90,0,0
WPRO,0,0,-90 !恢复工作平面
CSYS,0 !激活笛卡尔坐标系,即XYZ平面
L,2,3
ALLS,ALL
LGLUE,ALL
KGEN,2,4,,,0,0,Z1+W1
L,4,6 !构造辅助线10
AL,1,7,5,6
AL,6,2,9,3
AL,9,4,7,8
AL,5,2,4
AL,1,3,8
ASEL,S,,,1,5,1
VA,ALL
ALLS,ALL
VDRAG,1,,,,,,10 !第一段筒体
LDELE,10,,,1 !删除辅助线
ASEL,S,,,1
ASEL,A,,,6,10
AGLUE,ALL
ALLS,ALL
AGLUE,2,7
ALLS,ALL
AGLUE,3,9 !将封头与筒体侧面粘接
!内表面不能粘接,A4,A8;A5,A6
ALLS,ALL
KGEN,2,11,,,0,0,-W1
L,11,6
VDRAG,10,,,,,,10 !第二段筒体
LDELE,10,,,1 !删除辅助线
ASEL,S,,,10,15,1
AGLUE,ALL
ALLS,ALL
KGEN,2,15,,,0,0,Z1+W1
L,15,6
VDRAG,15,,,,,,10 !第三段筒体
LDELE,10,,,1 !删除辅助线
ASEL,S,,,15,20,1
AGLUE,ALL
ALLS,ALL
KGEN,2,20,,,0,0,-W1
L,20,6
VDRAG,20,,,,,,10 !第四段筒体
LDELE,10,,,1 !删除辅助线
ASEL,S,,,20,25,1
AGLUE,ALL
ALLS,ALL
KGEN,2,25,,,0,0,2*W1-2*Z1-L/2
L,25,6
VDRAG,25,,,,,,10 !第五段筒体
LDELE,10,,,1 !删除辅助线
ASEL,S,,,25,30,1
AGLUE,ALL
ALLS,ALL
VSEL,S,,,2,6,1
VGLUE,ALL
ALLS,ALL
!********第一条工字钢********
KGEN,2,11,,,-T1
KGEN,2,6,,,,,T3/2-W1/2
KGEN,2,31,,,T1+T2-W3
KGEN,2,32,,,,,W2/2-T3/2
KGEN,2,33,,,-T2
KGEN,2,34,,,,,-W2
KGEN,2,35,,,T2
KGEN,2,36,,,,,W2/2-T3/2
KGEN,2,37,,,W3-T1-T2
KGEN,2,38,,,,,T3/2-W1/2
L,11,6
L,6,31
L,31,32
L,32,33
L,33,34
L,34,35
L,35,36
L,36,37
L,37,38
L,38,39
L,39,15
LSEL,S,,,51,60,1
LSEL,A,,,10,23,13
LGLUE,ALL
AL,ALL !形成工字钢面
ALLS,ALL
VDRAG,31,,,,,,16 !拉伸成工字筋加强圈
VDELE,7
ADELE,43,,,1
ADELE,42
ADELE,44
ADELE,32
LDELE,62,,,1
LDELE,83,84,1,1
LGLUE,60,24,81,82
LGLUE,10,16,63,61
LSEL,S,,,62,82,2
LSEL,A,,,25
LGLUE,ALL
ALLS,ALL
AL,60,24,81,62
AL,10,16,63,82
LSEL,S,,,62,82,2
LSEL,A,,,25
AL,ALL
ALLS,ALL
ASEL,S,,,31,43
ASEL,A,,,13
VA,ALL
ALLS,ALL
!********第二条工字钢********
KGEN,2,20,,,-T1
KGEN,2,40,,,,,T3/2-W1/2
KGEN,2,51,,,T1+T2-W3
KGEN,2,52,,,,,W2/2-T3/2
KGEN,2,53,,,-T2
KGEN,2,54,,,,,-W2
KGEN,2,55,,,T2
KGEN,2,56,,,,,W2/2-T3/2
KGEN,2,57,,,W3-T1-T2
KGEN,2,58,,,,,T3/2-W1/2
L,20,40
L,40,51
L,51,52
L,52,53
L,53,54
L,54,55
L,55,56
L,56,57
L,57,58
L,58,59
L,59,25
LSEL,S,,,83,92,1
LSEL,A,,,39,61,22
LGLUE,ALL
AL,ALL !形成工字钢面
ALLS,ALL
VDRAG,44,,,,,,32 !拉伸成工字筋加强圈
VDELE,8
ADELE,45,,,1
ADELE,46
ADELE,56
ADELE,57
LDELE,93,95,1,1
LGLUE,96,97,40,92
LGLUE,61,32,115,116
LSEL,S,,,98,114,2
LSEL,A,,,93,94,1
LSEL,A,,,41
LGLUE,ALL
ALLS,ALL
AL,40,93,97,92
AL,61,32,94,115
LSEL,S,,,98,114,2
LSEL,A,,,93,94,1
LSEL,A,,,41
AL,ALL
ALLS,ALL
ASEL,S,,,44,56
ASEL,A,,,23
VA,ALL
ALLS,ALL
!********划分网格********
ALLS,ALL
VATT,1,,1 !定义单元及材料
!********筒体划分网格********
LSEL,S,,,1,5,4
LSEL,A,,,11,16,5
LSEL,A,,,19,24,5
LSEL,A,,,27,32,5
LSEL,A,,,35,40,5
LSEL,A,,,43,48,5
LESIZE,ALL,,,70,,,,,1
LSEL,S,,,12,17,5
LSEL,A,,,13,15,2
LSEL,A,,,28,33,5
LSEL,A,,,29,31,2
LESIZE,ALL,,,60,,,,,1
LSEL,S,,,44,49,5
LSEL,A,,,45,47,2
LESIZE,ALL,,,25,,,,,1
LSEL,S,,,20,25,5
LSEL,A,,,21,23,2
LSEL,A,,,36,41,5
LSEL,A,,,37,39,2
LESIZE,ALL,,,10,,,,,1
LSEL,S,,,46,50,4
LSEL,A,,,38,42,4
LSEL,A,,,30,34,4
LSEL,A,,,22,26,4
LSEL,A,,,14,18,4
LSEL,A,,,6,7,1
LESIZE,ALL,,,2,,,,,1
MSHAPE,0,3D
MSHKEY,1
VSEL,S,,,2,6,1
VMESH,ALL
ALLS,ALL
!********封头划分网格********
LSEL,S,,,2,3,1
LSEL,A,,,4,8,4
LESIZE,ALL,,,60,,,,,1
LSEL,S,,,6,7,1
LSEL,A,,,9
LESIZE,ALL,,,2,,,,,1
MSHAPE,0,3D
MSHKEY,1
ALLS,ALL
VMESH,1
ALLS,ALL
!********工字筋划分网格********
LSEL,S,,,63,81,2
LSEL,A,,,97,115,2
LESIZE,ALL,,,70,,,,,1
LSEL,S,,,52,58,6
LSEL,A,,,66,78,12
LSEL,A,,,84,90,6
LSEL,A,,,100,112,12
LESIZE,ALL,,,15,,,,,1
LSEL,S,,,87,106,19
LSEL,A,,,55,72,17
LESIZE,ALL,,,10,,,,,1
LSEL,S,,,98,114,16
LSEL,A,,,102,110,8
LSEL,A,,,83,91,8
LSEL,A,,,85,89,4
LSEL,A,,,64,80,16
LSEL,A,,,68,76,8
LSEL,A,,,51,59,8
LSEL,A,,,53,57,4
LESIZE,ALL,,,4,,,,,1
LSEL,S,,,62,82,20
LSEL,A,,,70,74,4
LSEL,A,,,54,56,2
LSEL,A,,,10,60,50
LSEL,A,,,93,94,1
LSEL,A,,,104,108,4
LSEL,A,,,86,88,2
LSEL,A,,,61,92,31
LESIZE,ALL,,,2,,,,,1
ALLS,ALL
VSWEEP,7
ALLS,ALL
VSWEEP,8
ALLS,ALL
!********加载********
/SOLU !进入求解器
ASEL,S,LOC,X,0
ASEL,A,LOC,Y,0
ASEL,A,,,30
DA,ALL,SYMM !对选中的面进行对称约束
ALLS,ALL
ASEL,S,,,4,8,4
ASEL,A,,,13,23,10
ASEL,A,,,18,28,10
SFA,ALL,1,PRES,P !对内表面加内压0.4MPa
ALLS,ALL
SOLVE !求解开始