基于ABAQUS法兰盘感应淬火残余应力场仿真 - 图文

基于ABAQUS法兰盘感应淬火残余应力场仿真 2012年7月

基于ABAQUS法兰盘感应淬火残余应力场仿真

摘要：表面感应淬火是一种工程中常用的热处理工艺，使用感应器来对工件的局部进行加热、迅速冷却，从而使工件表面产生残余应力，抵消工作载荷所产生的部分拉应力。本文就不同淬火温度下产生的残余应力大小的合理性及残余应力对缓解应力集中的作用两方面进行了分析，分析就过对提高产品质量、疲劳寿命及减少安全隐患有很大帮助。

关键字：法兰盘；残余应力；感应淬火；温度场；ABAQUS

法兰盘在制造过程中，将受到来自各种工艺等因素的作用与影响且其影响不能随之而完全消失，仍有部分作用与影响残留在构件内，则这种残留的作用与影响称为残留应力或残余应力。法兰盘是使管子与管子相互连接的零件，法兰连接使用方便，能够承受较大的压力。在工业管道中，法兰连接的使用十分广泛，法兰盘的作用就是使得管件连接处固定并密封，法兰连接是管道施工的重要连接方式。

2.1 思路

设定整个模型的初始温度20C，在分析步中令淬硬层区域的温度升高至某个温度值Thigh，其余区域的温度保持20C。这种温度差异会使高温区域产于压应力，相当于所要模拟的残余压应力。经过几次运算，就可以找到合适的Thigh，使法兰盘内圆角表面压应力与试验结果大致吻合。施加工作载荷时，保持上述温度场不变，就可以模拟在残余应力作用下应力场。 2.2 主要研究内容

①对淬硬层不同温度所产生的残余应力分析；

②对在相同外载荷，残余应力对缓和应力集中的作用进行分析。

o o

1、问题的提出

残余应力是一个不稳定的应力状态，构建受外力时，与残余应力相互作用，使构件某些局部呈现塑性变形。在实际工况中，法兰盘经过表面感应淬火后，淬硬层将会产生残余应力，影响着法兰的强度、寿命。本文就法兰盘不同淬火温度的残余应力场进行仿真分析，并分析此残余应力在缓和应力集中方面所起的作用，最后针对模型得出合理的淬火温度及产生的残余应力。

3、力学模型

本文以内孔直径为24mm，材料的弹性模量为2.1e5MPa，泊松比为0.3，线胀系数为1.35e-5C的法兰盘进行建模，其一端固定，另一端的整个端面受到向下的面载荷

o

2、研究思路及内容

1

基于ABAQUS法兰盘感应淬火残余应力场仿真 2012年7月

p=100MPa。由假定实验测得法兰盘的内圆角表面残余应力约为-430MPa。 3.1三维模型

尽管几何模型和残余应力场都具有轴对称，但端面上的载荷不具有轴对称，因此需要建立三维实体模型（基于对称性，只对1/2模型进行建模），如图1：

根据法兰盘工况，法兰盘的底部使用全约束，法兰盘对称面的各个区域使用ZSYMM。

对称面 淬硬层的分界线

图3 边界条件模型

3.5 设置分析步

本模型属于静态响应，分析步类型为

图1 三维模型

Static，general。

①第一个分析步：淬火硬层区域的温度分别升高至100C 、120C、140C、160C,其余区域的温度保持20C，不施加外载荷，仿真工件的残余应力场。

②第二个分析步：保持上述温度场不变(有残余应力)，施加外载荷p=100MPa。

③第三个分析步：整个工件的温度都变为20 oC（无残余应力），保持外载荷p=100MPa，从而得到没有残余应力时的应

o

o

o

o

o

3.2定义材料、截面属性

设置杨氏模量为210000，泊松比为0.3，线胀系数1.35e-5。 3.3划分网格

由于关心的是应力集中部位的应力状态，所以淬硬层使用C3D20R单元类型。

淬硬层 力场，用来与第二个分析步结果作比较。

图2 网格划分

3.4边界条件

2

基于ABAQUS法兰盘感应淬火残余应力场仿真 2012年7月

600最小主应力σminp5004003002001000100120不同淬火温度140160635号节点

图4 定义温度场载荷模型

图6 第一分析步（无外载荷）

由上图表明最小主应力随着淬火温度的上升而增大，淬火温度越高残余应力就越

4、有限元计算结果分析

由分析处理后可以得到图5法兰盘圆角

大，过大的残余应力（实验得知本模型残余应力为σ

minp

=430MPa）既不能平衡工作载荷

产生的拉应力，还会对法兰盘的刚度、寿命

表面残余应力的曲线图（第一分析步，温度场，无应力场）.

σ

minp

产生影响。当淬火的温度在120C时，

=433MPa与实验相符。因此，将淬火硬

o

o

层区域的温度设为120C是比较合理的。

635号节点

641号节点 图5 取法兰盘内圆角表面路径曲线

从图中并结合模型得出，最小主应力在635号节点处值为σ

minp

图7 淬火层圆角处最小主应力

=433MPa。

分析残余应力的缓解应力集中方面的作用，在第一个分析步中，内圆角641号节点应力为σ

minp

不同淬火温度下残余应力的大小各异，从而对法兰盘带来不同的影响。下图表明100C、 120C、 140C 、160C温度下，淬火温度的产生的残余应力的大小。（法兰盘圆角处635号节点残余应力曲线图。

o

o

o

o

= -238MPa。

3

基于ABAQUS法兰盘感应淬火残余应力场仿真 2012年7月

图8 温度场、应力场下最大主应力

图10 635、641节点最大、最小主应力曲线

在第二个分析步中（存在残余应力），635号节点的最大主应力σ641号节点最大应力为σ

minp

由上图，再结合分析模型可见残余应力显著降低了应力集中处的最大主应力，在此节点处，最大主应力降低的量676-469=209MPa，大致等于残余应力的值238MPa。

= -1MPa，

minp

= -469MPa。

5结论

本文通过对法兰盘感应淬火残余应力场仿真分析了不同淬火温度下的残余应力及残余应力缓解应力集中两方面内容，得出

图9 应力场下最大主应力（无温度场）

以下结论：

1）残余应力主要发生在法兰盘淬硬层内圆角表面处。

2）最小主应力随着淬火温度的上升而增大，淬火温度越高残余应力就越大，针对本模型淬火温度应该约为120C比较合理。

3）在合理的温度下，残余应力显著降低了应力集中处的最大主应力，适当的残余应力能缓解法兰盘的拉应力，对提高法兰盘的寿命起到了很大作用。

o

第三个分析步（没有残余应力），此641号节点的最大主应力σ节点最大应力σ

minp

minp

= 676MPa，635号

= -178MPa。

通过有限元分析不同淬火温度下的残余应力及残余应力缓解应力集中两方面等处最小、最大主应力图（635、641号节点），如图10.

4

基于ABAQUS法兰盘感应淬火残余应力场仿真 2012年7月

上述方法是比较简单，与实际工况还存

在一定差距，精确度还不够高，对结果有一

定的影响。在后续的工作中，结合工况，需

要建立复杂的法兰盘传热分析和热弹塑性

分析。

参考文献

[1] 付学义，刘莉.残余应力的危害及对策

[J].包钢科技，2006.11.

[2] 石亦平，周玉蓉.ABAQUS有限元分析实

例详解.机械工业出版社[M]。2012.4.

[3] 郭秋艳，庞浩，邓磊.基于ANSYS的轿

车轮毂轴承内法兰盘的力矩刚性分析

[J]。机械设计及制作.2010.7.

[4] 刘展.ABAQUS6.6基础教程与实例详解.

中国水利水电出版社[M].2008.2.

命令流：

*Heading 34768, -28.4204025, 62.9329185, ** Job name: Job-1 Model name: Model-1 0. ** Generated by: Abaqus/CAE 6.10-1 34769, -29.1488552, 61.7691116, *Preprint, echo=NO, model=NO, 0. history=NO, contact=NO 34770, -29.7671928, 60.5814781, ** 0. ** PARTS *Element, type=C3D15 ** 1, 675, 665, 733, 3304, 3294, 3362, *Part, name=flange 9024, 9023, 9022, 9025, 9026, 9027, 9029, *Node 9028, 9030 1, -26.7999992, 66.2205734, 2, 725, 724, 743, 3354, 3353, 3372, 0. 9033, 9032, 9031, 9034, 9035, 9036, 9038, 2, -1.74696732, 53.5455704, 9037, 9039 0. 3, 3304, 3294, 3362, 3397, 3387, 3455, 3, 4.9000001, 53.5455704, 9025, 9026, 9027, 9040, 9041, 9042, 9044, 0. 9043, 9045 ------------------------------------ 4, 3354, 3353, 3372, 3447, 3446, 3465, 34767, -26.7999992, 64.9721298, 9034, 9035, 9036, 9046, 9047, 9048, 9050, 0. 9049, 9051 ------------------------------------ 58, 5865, 5864, 5883, 5958, 5957, 5976, 9358, 9359, 9360, 9370, 9371, 9372, 9374, 9373, 9375

6 展望

5

基于ABAQUS法兰盘感应淬火残余应力场仿真 2012年7月

59, 5908, 5898, 5966, 768, 758, 826, 9364, 9365, 9366, 9376, 9377, 9378, 9380, 9379, 9381 60, 5958, 5957, 5976, 818, 817, 836, 9370, 9371, 9372, 9382, 9383, 9384, 9386, *Nset, nset=Set-fix, instance=flange-1 9, 10, 21, 22, 78, 79, 80, 81, 82, 83, 84, 85, 236, 237, 238, 239

240, 241, 242, 243, 244, 9385, 9387

*Element, type=C3D20R

61, 701, 737, 84, 83, 3330, 3366, 896, 897, 9391, 9390, 9389, 9388, 9392, 9393, 9394,

9395, 9397, 9396, 9398, 9399

62, 673, 738, 703, 736, 3302, 3367, 3332, 3365, 9403, 9402, 9401, 9400, 9404, 9405, 9406,

9407, 9409, 9408, 9410, 9411

63, 12, 88, 704, 13, 351, 1156, 3333, 380, 9415, 9414, 9413, 9412, 9416, 9417, 9418,

9419, 9421, 9420, 9422, 9423

---------------------------------- 7769, 8981, 1297, 1298, 8980, 2688, 36, 35, 2687, 34388, 9650, 34387, 34372, 34770, 9649, 34769,

34754, 34753, 9644, 9653, 34756

7770, 8967, 1296, 1297, 8981, 2674, 37, 36, 2688, 34188, 9642, 34388, 34369, 34570, 9635, 34770,

34751, 34572, 9646, 9644, 34753

*Nset, nset=_PickedSet2, internal, generate

1, 34770, 1

*Elset, elset=_PickedSet2, internal, generate

1, 7770, 1 ** Section: Section-1

*Solid Section, elset=_PickedSet2, material=Material-1 ,

*End Part ** **

** ASSEMBLY **

*Assembly, name=Assembly **

*Instance, name=flange-1, part=flange *End Instance **

245, 246, 247, 248, 249, 250, 251, 252, 253, 254, 255

256, 257, 258, 259, 260, 261, 262, 263, 264, 265, 266, 267, 268, 269, 270, 271

------------------------------------ 23027, 23028, 23029, 23030, 23031, 23032, 23033, 23034, 23035, 23036, 23037, 23038, 23039, 23040, 23041, 23042

23043, 23044, 23045, 23046, 23047, 23048, 23049, 23050, 23051, 23052, 23053, 23054

*Elset, elset=Set-LOWt, instance=flange-1, generate 1, 3870, 1

*Elset, elset=_Surf-1oad_S4, internal, instance=flange-1

79, 102, 206, 229, 333, 356, 460, 483, 587, 610, 714, 737, 841, 864, 968, 991

1095, 1118, 1222, 1245, 1349, 1372, 1476, 1499, 1603, 1626, 1730, 1753, 1857, 1880, 1984, 2007

2111, 2134, 2238, 2261, 2365, 2388, 2492, 2515, 2619, 2642, 2746, 2769, 2873, 2896, 3000, 3023

3127, 3150, 3254, 3277, 3381, 3404, 3508, 3531, 3635, 3658, 3762, 3785 *Elset, elset=_Surf-1oad_S3, internal, instance=flange-1, generate 89, 3772, 127

*Surface, type=ELEMENT, name=Surf-1oad _Surf-1oad_S4, S4 _Surf-1oad_S3, S3 *End Assembly **

** MATERIALS **

*Material, name=Material-1 *Elastic 210000., 0.3 *Expansion 1.35e-05, **

** BOUNDARY CONDITIONS **

** Name: BC-1 Type:

Symmetry/Antisymmetry/Encastre *Boundary

Set-fix, ENCASTRE

6

基于ABAQUS法兰盘感应淬火残余应力场仿真 2012年7月

** Name: BC-2 Type:

Symmetry/Antisymmetry/Encastre *Boundary

Set-symm, ZSYMM **

** PREDEFINED FIELDS **

** Name: Predefined Field-2 Type: Temperature

*Initial Conditions, type=TEMPERATURE Set-LOWt, 20. ** Name: Predefined Field-HIGHt Type: **

** OUTPUT REQUESTS **

*Restart, write, frequency=0 **

** FIELD OUTPUT: F-Output-1 **

*Output, field, variable=PRESELECT **

** HISTORY OUTPUT: H-Output-1 **

*Output, history, variable=PRESELECT Temperature

*Initial Conditions, type=TEMPERATURE Set-highT, 20. **

------------------------------------ **

** STEP: Step-NOLOAD **

*Step, name=Step-NOLOAD *Static

1., 1., 1e-05, 1. **

** PREDEFINED FIELDS ** ** Name: Predefined Field-HIGHt Type: Temperature *Temperature Set-highT, 120. **

** OUTPUT REQUESTS **

*Restart, write, frequency=0 **

** FIELD OUTPUT: F-Output-1 **

*Output, field, variable=PRESELECT **

** HISTORY OUTPUT: H-Output-1 **

*Output, history, variable=PRESELECT *End Step **

------------------------------------ **

** STEP: Step-HIGHWITHLOAD **

*Step, name=Step-HIGHWITHLOAD *Static

1., 1., 1e-05, 1. **

** LOADS **

** Name: Load-surface Type: Surface traction *Dsload

Surf-1oad, TRSHR, 100., 0., -1., 0.

*End Step **

------------------------------------**

** STEP: Step-LOWWITHLOAD **

*Step, name=Step-LOWWITHLOAD *Static

1., 1., 1e-05, 1. **

** PREDEFINED FIELDS ** ** Name: Predefined Field-HIGHt Type: Temperature *Temperature Set-highT, 20. **

** OUTPUT REQUESTS **

*Restart, write, frequency=0 **

** FIELD OUTPUT: F-Output-1

**

*Output, field, variable=PRESELECT **

** HISTORY OUTPUT: H-Output-1 **

*Output, history, variable=PRESELECT *End Step

7

本文来源：https://www.bwwdw.com/article/y62g.html