感覺自己建的模型沒有什么問題,模型由三部分組成（實體單元），封頭+筒體+環向加強圈。因為是存在曲面，在建模時封頭和筒體連接處存在小空隙。但是能正常出來結果，所以我認為模型基本沒啥問題。加載應該也沒啥問題，就是內壓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 !求解開始