浙大化工机械复杂系统分析与建模课程作业

1、复杂系统分析及建模课程大作业1、计算条件1.1结构示意图容器的结构示意图如图1所示：图1 容器结构示意图其中球型封头与筒体的过渡采用堆焊球封头的方式。1.2计算要求（1）运用有限元法预测爆破压力（2）若容器壁温从200降至室温20，同时内压满足如下条件，求疲劳寿命，并说明热应力利弊。2、容器爆破过程分析图2 厚壁圆筒中压力与变形关系塑性材料制造的压力容器的爆破过程如图2所示。在弹性变形阶段（OA线段），器壁应力较小，产生弹性变形，内压与容积变化量成正比，到A点时容器内表面开始屈服，与A点对应的压力为初始屈服压力；在弹塑性变形阶段（AC线段），随着内压的继续提高，材料从内壁向外壁屈服，此时，一方

2、面因塑性变形而使材料强化导致承压能力提高，另一方面因壁厚不断减薄而使承压能力下降，但材料强化作用大于壁厚减薄作用，到C点时两种作用已接近，C点对应的压力是容器所能承受的最大压力，称为塑性垮塌压力；在爆破阶段（CD线段），容积突然急剧增大，使容器继续膨胀所需要的压力也相应减小，压力降落到D点，容器爆炸，D点所对应的压力为爆破压力。对于内压容器，爆破过程中，内压和容积变化量的关系与材料塑性、容器容积和壁厚有关。对于脆性材料，不会出现弹塑性变形阶段。虽然塑性垮塌压力大于爆破压力，但工程上往往把塑性垮塌压力视为爆破压力。本报告也近似认为塑性垮塌压力就是爆破压力。3、爆破压力的经验计算公式Faupel厚

3、壁圆筒爆破压力计算公式：，,，4、爆破压力数值计算4.1计算简化根据结构的特点，从偏于危险和节约计算成本的角度，建模时只考虑整体的1/4，几何模型如图2所示:图2 几何模型4.2材料属性、单元选择及边界条件材料为Q345R，弹性模量205000MPa，泊松比0.3，热膨胀系数，实际应力应变关系曲线如图3所示：图3 Q345真实应力应变关系曲线采用有限元软件Abaqus6.9进行数值模拟计算，单元采用具有沙漏控制的减缩积分单元C3D8R。有限元模型的边界条件如图4所示：图4边界条件4.3网格收敛性验证：对结构采用两种不同大小的网格进行划分，以爆破压力作为评价标准，考察网格的收敛性。本分析采用的网

4、格模型和较细的网格划分模型的网格图片分别如图5所和图6所示。它们相应的爆破压力计算结果如表1所示。图5 本报告采用网格图6 较细的网格表 1网格模型单元数爆破压力本报告采用的网格1934549.40较细网格2461149.41从表1可知，本报告采用的网格划分符合收敛性条件。4.4数值计算分别采用弧长法和非线性稳定性算法计算容器的爆破压力。4.4.1弧长法采用弧长法得到的结构在爆破时的Mises应力云图如图7所示。容器的爆破压力为49.40MPa。爆破压力与弧长关系曲线如图8所示：图7容器爆破时Mises应力云图图8 爆破压力与弧长关系曲线4.4.2非线性稳定性算法采用ABAQUS的dynami

5、c/implicit模块，运用非线性稳定算法得到的结构的Mises应力云图如图9所示，载荷与弧长关系曲线如图10所示：图9 非线性稳定算法得到的容器Mises应力云图图10 爆破压力曲线图得到容器的爆破压力为49.33MPa。4.5有无开孔及开孔大小对爆破压力的影响为研究有无开孔及开孔大小对爆破压力的影响，考虑如表2所示的几个模型：表 2模型名称模型说明模型1原始模型模型2仅封头开孔模型3无开孔模型4开孔尺寸变大计算结果如表3所示：表 3模型名称爆破压力（MPa）爆破位置模型149.40筒体接管相贯处内壁模型249.41封头接管相贯处内壁模型349.45筒体模型448.93筒体接管相贯处内壁模

6、型2、模型3和模型4在爆破时相应的Mises应力云图分别如图11、12和13所示。图11模型2爆破时Mises应力云图图12模型3爆破时Mises应力云图图13模型4爆破时Mises应力云图在前面的计算结果可以看出，在本文采用的计算模型的前提下，有无开孔，及开孔大小对爆破压力影响很小，但是对爆破位置却影响很大。4.6倒角大小对爆破压力的影响在原始模型的基础上，研究筒体与接管相贯处内倒角大小对爆破压力的影响，考虑如下的几个模型：表 4模型名称模型说明模型1原始模型倒角R=10mm模型5倒角R=15mm模型6倒角R=5mm计算结果如表5所示：表 5爆破压力（MPa）爆破位置49.40筒体接管相贯处

7、内壁49.399筒体接管相贯处内壁49.41筒体接管相贯处内壁模型5和模型6在爆破时相应的Mises应力云图分别如图14和15所示。图14模型5爆破时Mises应力云图图15模型6爆破时Mises应力云图对于该容器来说，倒角的大小变化对爆破压力影响很小，但对爆破位置影响很大。5、极限载荷数值计算5.1极限分析方法极限分析的基本概念：在加载过程中，结构中的高应力区首先进入塑性，当载荷继续增加时塑性区便不断扩大，同时还出现应力重分布现象。当载荷增大到某一极限值时，由理想塑性材料制成的结构将变成不稳定的几何可变机构，从而丧失承载能力，出现不可限制的塑性流动，此时载荷不变但应变能无限增加，这种状态称为

8、塑性极限状态，相应的载荷称为极限载荷（全屈服载荷）。假如只计算极限载荷而不计及极限状态到达以前的变形过程，不考虑加载的历史，这种分析方法就称为极限分析法。极限分析基本假设：材料为理想弹塑性，变形足够小，满足比例加载条件。极限载荷确定的方法：采用清华大学陆明万教授等提出的零曲率法，该方法具有分散性小的特点。5.2材料属性材料基本参数：材料为Q345R，弹性模量205000MPa，泊松比0.3，初始屈服应力。5.3极限载荷的理论计算值筒体整体屈服时，按Mises屈服失效判据，极限载荷理论计算值为：5.4数值计算分别采用Newton法和弧长法计算容器的极限载荷。5.4.1Newton法在容器内壁所有

9、与介质接触的表面施加压力载荷p=60MPa，采用Newton法计算得到的应力云图如图16所示：图16 极限分析应力云图载荷-最大点等效塑性应变曲线，如图17所示：图17 LPF曲线图采用零斜率准则确定极限载荷： 60×0.610404=36.625.4.2弧长法在容器内壁所有与介质接触的表面施加压力载荷p=60MPa，采用risk方法计算得到的应力云图如图18所示：图18 极限分析容器的Mises应力云图载荷-最大点等效塑性应变曲线，如图19所示：图19 LPF-曲线图采用零斜率准则确定极限载荷： 60×0.592061=35.526、应力分类法6.1应力分类的概念一次应力

10、超过材料屈服极限时，将会引起过量的总体塑性变形而造成结构破坏。分析设计一次加二次应力强度的控制值采用了，即以结构是否安定为判断依据。6.2材料属性材料基本参数：材料为Q345R，弹性模量205000MPa，泊松比0.3。6.3名义弹性应力计算方法材料为线弹性，取单位压力1，计算得到的应力云图如图20所示：图20 名义弹性应力分析应力云图及分析路径选取关键路径并对其进行应力评定，此处仅列出最危险路径的应力线性化数据：表 6应力分类计算值许用值11.92238.520.0119.1447724.92采用应力分类法计算得到的允许压力为20.01。7、疲劳寿命计算7.1基本概念工作循环：由初始状态进入

11、新状态，随后又回到初始状态开始点的过程。包括启动停止循环、正常工作循环和设计中必须考虑的任何紧急状态或异常情况由起始到恢复的循环。本次计算时仅仅考虑正常工作循环。应力循环：指应力由初始值开始，经过代数最大值和代数最小值，然后又返回初始值的循环。一个工作循环可以引起一个或多个应力循环。根据容器的工作条件，一个工作循环包括两个应力循环。疲劳评定是以结构应力循环引起的应力差波值为基础的。根据工作循环和应力循环的定义，结合本题给定的条件。此工作循环包括2个应力循环。疲劳评定方法：当所考虑点的主应力方向在循环中变化时，应按如下步骤确定交变应力强度幅：（1）确定所考虑点在整个应力循环中与时间相对应的包括总

12、体和局部结构不连续以及热效应所引起的六个应力分量；（2）选取循环条件的极端点（代数最大值或最小）所对应的时刻，记下这一时刻的各应力分量；（3）循环中的每一时刻，从每个与该时刻对应的应力分量中减去i时刻的相应的应力分量记下波动应力分量；（4）在循环的每一时刻，计算由六个波动应力分量所导出的主应力波动范围。这些波动主应力的方向虽然在循环中变化，但其编号应保持不变；（5）计算在整个应力循环中，相对于时间的波动主应力差；（6）确定各波动主应力差的最大波动范围，计算交变应力强度幅。7.2内加热对容器疲劳的影响分析对于厚壁圆筒在内压与内加热同时作用下的应力特点来说，如果由内压引起的应力与内加热温差所引起的

13、热应力同时存在，在弹性变形前提下筒壁的总应力则为两种应力的叠加。内加热情况下内壁应力叠加后得到改善，而外壁应力有所恶化。外加热时则相反，内壁应力恶化，而外壁应力得到很大改善。通常内压作用下应力最大点都出现在接管和筒体相贯处的内壁，疲劳评定时也通常选取此点。因此从疲劳的角度分析，内加热时对冲内压的容器的疲劳是有益的。相对应的是对与容器受外压的情况，内加热对容器的疲劳是不利的。7.3采用ABAQUS软件对容器进行的疲劳计算与评定分析方法：考虑热应力影响时，采用间接耦合方法，即先计算温度场再计算应力场的方法。温度场计算时仅仅考虑热传导。考虑几何非线性效应。采用固定时间步长以0.2小时为单位增量单位。

14、对容器内壁施加随时间变化的温度载荷，外壁保持温度恒定为20，且仅考虑热传导效应。不考虑温度时的容器疲劳评定不考虑温度时容器在第1小时和第3小时的应力云图如图21所示：图21 第1小时和第3小时对应的容器名义弹性应力云图以节点1389为疲劳评定点，相关应力数据如表7所示：表 7应力分量任意时刻应力(MPa)最大或最小应力(MPa)应力差(MPa)S1119.1852-19.185238.3704S22381.043-381.043762.086S3324.1375-24.137548.275S1230.7763-30.776361.5526S137.20228-7.2022814.40456S2

15、3-5.992245.99224-11.9845相应三个主应力值为767.44，57.72，25.56(MPa)主应力波动差值为767.44-25.56=741.88MPa对应的交变应力幅为370.1MPa采用JB4732疲劳载荷计算公式：对应的应力循环次数约为3472次。采用ansys的疲劳模块，得到：评定点在事件中各载荷点的应力值：LIST FATIGUE STRESSES FOR LOCATIONS 1 TO 5 IN STEPS OF 1 EVENTS 1 TO 3 - LOADS 1 TO 10 * POST1 FATIGUE STRESSES * LOCATION 1 NODE 1

16、389 EVENT 1 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 1 SX,SY,SZ 19.180 381.04 24.140 SXY,SYZ,SXZ,TEMP 30.780 -5.9900 7.2000 0.0000 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 2 SX,SY,SZ -19.180 -381.04 -24.140 SXY,SYZ,SXZ,TEMP -30.780 5.9900 -7.2000 0.0000 疲劳评定 PERFORM FATIGUE CALCU

17、LATION AT LOCATION 1 NODE 0 * POST1 FATIGUE CALCULATION * LOCATION 1 NODE 1389 EVENT/LOADS 1 1 AND 1 2 PRODUCE ALTERNATING SI (SALT) = 371.94 WITH TEMP = 0.0000 CYCLES USED/ALLOWED = 0.1000E+05/ 3445. = PARTIAL USAGE = 2.90242 CUMULATIVE FATIGUE USAGE = 2.90242即采用ansys疲劳模块计算得到的疲劳次数为3445次7.3.2考虑温度时的疲

18、劳计算与评定考虑温度时容器在第0.2小时的温度分布云图如图22所示。第1小时、第3小时、第5小时、第7小时对应的应力云图分别如图2326所示。图22 第0.2小时对应的容器温度云图图23 第1小时对应的容器应力云图图24 第3小时对应的容器应力云图图25 第5小时对应的容器应力云图图26 第7小时对应的容器应力云图以节点1389为疲劳评定点，相关应力数据如表8所示：表 8应力分量任意时刻应力(MPa)最大或最小应力(MPa)应力差(MPa)S1118.72-21.5340.25S22376.18-405.36781.54S3323.24-28.6151.85S1230.33-32.9963.3

19、2S136.96E+00-8.4115.37S235.956.2-0.25相应的三个主应力值为786.91，60.89，25.84(MPa)最大主应力波动差值为786.91-25.84=761.07MPa采用JB4732疲劳载荷计算公式：对应的应力循环次数约为3190次。采用ansys的疲劳模块，得到：评定点在事件中各载荷点的应力值：LIST FATIGUE STRESSES FOR LOCATIONS 1 TO 5 IN STEPS OF 1 EVENTS 1 TO 3 - LOADS 1 TO 10 * POST1 FATIGUE STRESSES * LOCATION 1 NODE 13

20、89 EVENT 1 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 1 SX,SY,SZ 18.720 376.18 23.240 SXY,SYZ,SXZ,TEMP 30.330 5.9500 6.9600 0.0000 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 2 SX,SY,SZ -21.530 -405.36 -28.610 SXY,SYZ,SXZ,TEMP -32.990 6.2000 -8.4100 0.0000 疲劳评定PERFORM FATIGUE CALCULAT

21、ION AT LOCATION 1 NODE 0 * POST1 FATIGUE CALCULATION * LOCATION 1 NODE 1389 EVENT/LOADS 1 1 AND 1 2 PRODUCE ALTERNATING SI (SALT) = 380.54 WITH TEMP = 0.0000 CYCLES USED/ALLOWED = 0.1000E+05/ 3203. = PARTIAL USAGE = 3.12215 CUMULATIVE FATIGUE USAGE = 3.12215即采用ansys疲劳模块计算得到的疲劳次数为3203次7.4采用ansys软件对容器

22、进行的疲劳计算与评定分析方法：采用ansys直接软件12.0进行分析。考虑温度时采用直接耦合法。考虑材料非线性、几何非线性效应。采用固定时间步长以0.2小时为单位增量单位。对容器内壁施加随时间变化的温度载荷，外壁保持温度恒定为20，且仅考虑热传导效应。有限元模型：单元为solid5。分析单元为6993个，通过收敛性验证。网格模型如图27所示：图27 容器有限元网格图网格质量检查：SUMMARIZE SHAPE TESTING FOR ALL SELECTED ELEMENTS - <<<<<< SHAPE TESTING SUMMARY >>&g

23、t;>>> <<<<<< FOR ALL SELECTED ELEMENTS >>>>>> - - | Element count 6993 SOLID5 | - Test Number tested Warning count Error count Warn+Err % - - - - - Jacobian Ratio 6993 0 0 0.00 % Warping Factor 6993 0 0 0.00 % Any 6993 0 0 0.00 % -边界条件：边界条件如图28所示：图28 容器边界

24、条件示意图不考虑温度时的容器疲劳评定不考虑温度时容器在第1小时和第3小时的应力云图如图29所示：图29 第1或第3小时对应的容器应力云图评定点在事件中各载荷点的应力值：LIST FATIGUE STRESSES FOR LOCATIONS 1 TO 5 IN STEPS OF 1 EVENTS 1 TO 5 - LOADS 1 TO 10 * POST1 FATIGUE STRESSES * LOCATION 1 NODE 4792 EVENT 1 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 1 SX,SY,SZ 0.98587E+07

25、 -0.58100E+07 0.35948E+09 SXY,SYZ,SXZ,TEMP 0.21767E+07 0.22379E+06 -0.40174E+08 0.78886E-30 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 2 SX,SY,SZ 0.20798E-04 -0.95407E-06 0.35215E-05 SXY,SYZ,SXZ,TEMP 0.31515E-05 -0.42409E-05 0.63062E-05 0.78886E-30 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT

26、 1 LOAD 3 SX,SY,SZ -0.98585E+07 0.58101E+07 -0.35948E+09 SXY,SYZ,SXZ,TEMP -0.21767E+07 -0.22376E+06 0.40174E+08 0.78886E-30 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 4 SX,SY,SZ -0.20616E-04 0.26279E-04 0.22150E-05 SXY,SYZ,SXZ,TEMP -0.31923E-05 -0.43002E-05 -0.68344E-05 0.78886E-30 (TOTAL ST

27、RESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 5 SX,SY,SZ -0.98585E+07 0.58101E+07 -0.35948E+09 SXY,SYZ,SXZ,TEMP -0.21767E+07 -0.22376E+06 0.40174E+08 0.78886E-30 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 6 SX,SY,SZ -0.84352E-05 0.51755E-05 0.52036E-05 SXY,SYZ,SXZ,TEMP -0.35773E-05 0.49984E-05 -

28、0.82783E-05 0.78886E-30 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 7 SX,SY,SZ 0.98586E+07 -0.58100E+07 0.35948E+09 SXY,SYZ,SXZ,TEMP 0.21766E+07 0.22379E+06 -0.40174E+08 0.78886E-30 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 8 SX,SY,SZ 0.19627E-04 -0.12463E-04 0.13471E-05 SXY,SYZ,SXZ,

29、TEMP 0.76218E-05 -0.51002E-07 0.58006E-05 0.78886E-30疲劳评定：PERFORM FATIGUE CALCULATION AT LOCATION 1 NODE 0 * POST1 FATIGUE CALCULATION * LOCATION 1 NODE 4792 EVENT/LOADS 1 3 AND 1 7 PRODUCE ALTERNATING SI (SALT) = 0.37026E+09 WITH TEMP = 0.78886E-30 CYCLES USED/ALLOWED = 0.1000E+05/ 3495. = PARTIAL

30、USAGE = 2.86082 CUMULATIVE FATIGUE USAGE = 2.86082考虑温度时的容器疲劳评定考虑温度时容器在第0.2小时的温度分布云图如图30所示。第1小时、第3小时、第5小时、第7小时对应的应力云图分别如图3134所示。图30 第0.2小时对应的容器温度云图图31 第1小时对应的容器应力云图图32 第3小时对应的容器应力云图图33 第5小时对应的容器应力云图图34 第7小时对应的容器应力云图评定点在事件中各载荷点的应力值：LOCATION 1 NODE 4792 EVENT 1 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVE

31、NT 1 LOAD 1 SX,SY,SZ 0.12759E+08 -0.95293E+07 0.33127E+09 SXY,SYZ,SXZ,TEMP -0.23766E+07 0.10794E+06 -0.35301E+08 177.50 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 2 SX,SY,SZ 0.24866E+07 -0.31887E+07 -0.24179E+08 SXY,SYZ,SXZ,TEMP -0.39028E+07 -99587. 0.41768E+07 155.00 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 3 SX,SY,SZ -0.77863E+07 0.31533E+07 -0.37963E+09 SXY,SYZ,SXZ,TEMP -0.54289E+07 -0.30631E+06 0.43655E+08 132.50 (TOTAL STRESS - ITEMS 1-7) LOCATION 1 EVENT 1 LOAD 4 SX,SY,SZ 0.1657