实体单元建模出现的问题? 20

浏览:1062

感觉自己建的模型没有什么问题,模型由三部分组成(实体单元),封头+筒体+环向加强圈。因为是存在曲面,在建模时封头和筒体连接处存在小空隙。但是能正常出来结果,所以我认为模型基本没啥问题。加载应该也没啥问题,就是内压0.4MPa。最后结果出来最大应力在环向加强圈上,比较奇怪,觉得应该实在筒体上。所以又质疑自己的模型了。求大佬帮忙查看一下,看看问题出在哪里。 

下面是两张结果图,第一张是全部实体应力图,第二张是单看加强圈上的应力图。发现最大应力在加强圈上

方形封头(壁厚10,构造圆3000).jpg方形封头(壁厚10,构造圆3000)单看工字筋加强圈.jpg

这是代码

!***************环境设置*********************

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 !求解开始



邀请回答 我来回答

当前暂无回答

回答可获赠 200金币

没解决?试试专家一对一服务

换一批