CN106294957B - 一种提高回转支承高温环境下承载能力的方法 - Google Patents
一种提高回转支承高温环境下承载能力的方法 Download PDFInfo
- Publication number
- CN106294957B CN106294957B CN201610625931.XA CN201610625931A CN106294957B CN 106294957 B CN106294957 B CN 106294957B CN 201610625931 A CN201610625931 A CN 201610625931A CN 106294957 B CN106294957 B CN 106294957B
- Authority
- CN
- China
- Prior art keywords
- pivoting support
- bearing capacity
- raceway
- contact
- radius
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Rolling Contact Bearings (AREA)
Abstract
本发明公开了一种提高回转支承高温环境下承载能力的方法,属于回转支承设计领域。本发明利用有限元建模模拟高温环境下回转支承的工况,求解出滚动体和滚道接触区域的应力分布情况,并分析单一结构因素对回转支承承载能力的影响,再利用正交试验分析多个结构因素共同改变时对回转支承承载能力的影响,得出回转支承承载能力最强的参数组合,从而提高回转支承的承载能力。本发明可以模拟更接近实际的工况,计算结果更加精确,且大大地减少了工作量,更加高效、快速、经济。
Description
技术领域
本发明属于回转支承设计领域,更具体地说,涉及一种提高回转支承高温环境下承载能力的方法。
背景技术
回转支承即转盘轴承,其结构组成主要包括滚动体和内、外滚圈等。因其可以承载较大倾覆力矩和轴、径向载荷的共同作用,同时兼具运转方便灵活、结构简单紧凑、装拆方便快捷等应用优势,所以逐渐在各种机械领域得到广泛运用。在一些高温工况(如工频炉、铁水钢包、熔炼炉转包等)中使用的回转支承,不仅要承受高温产生的巨大热应力,同时还要承载回转平台等物体自重的共同作用,一旦超过其许可的承载应力值,极容易造成回转支承滚道损坏,从而影响整个机械以及整条作业线的正常运转。
回转支承的承载能力与滚动体和滚道接触区域的应力场密切相关。目前回转支承接触区域的应力分析主要有两种方法:基于Herts接触理论的解析法和有限元法;基于Herts接触理论的解析法需要对模型进行大量的简化,并且不考虑材料的塑性变形、滚动体与滚道之间的摩擦及滚圈的表面处理等等,因此这种解析法的误差较大。
有限元法是近些年来随着计算机技术的快速发展而形成的一种利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟的数值分析方法,这种方法可以模拟更接近实际的工况,计算结果更加精确。当前利用有限元软件已能够比较方便地对常温下的回转支承承载性能进行分析,然而要想准确地反映高温环境下回转支承的承载性能则比较困难。主要是由于回转支承模型一般较大,需要将温度场与应力场进行耦合来模拟更接近实际的工况,而加上热力耦合算法求解速度相对较慢,如果想得到较为精确的结果,需要花费大量的计算时间。且对高温回转支承的优化改进需考虑更多因素,这样就需要对很多组不同模型进行热力耦合分析,将会浪费大量的精力和时间在有限元求解计算中,是相当不值得的。
经检索,在常温下利用有限元法设计回转支承的方案已有公开,如中国专利号ZL201210102961.4,授权公告日为2014年5月7日,发明创造名称为:一种风力发电机主轴调心滚子轴承及其设计方法;该申请案的轴承设计过程包括以下步骤:a.建立坐标系;b.计算主轴受力;c.计算调心滚子轴承受力;d-1.计算接触应力;d-2.取下限值A1;d-3.取上限值A2;d-4.利用差值法进行计算;d-5.计算满足|Ui-Vi|<3Ui/100且|Ui-Vi|<3Vi/100的比例系数Ai,设计出符合要求的轴承。该申请案有效减小了滚子与滚道的接触应力,采用有限元和差值算法,将两排滚子接触应力的差值控制在较小的3%以内,解决了滚子受力不均衡问题,提高了轴承使用寿命。但正如上文分析,该申请案并不适用于高温环境下回转支承的承载性能计算,仍需有针对性的重新设计一种提高回转支承高温环境下承载能力的方案。
发明内容
1.发明要解决的技术问题
为了准确地计算高温环境下回转支承接触区域的应力分布并对其承载能力进行优化,本发明提供了一种高温环境下承载能力强的回转支承及提高其承载能力的方法;本发明利用有限元建模模拟高温环境下回转支承的工况,求解出接触区域的应力分布情况,并综合考虑回转支承接触角、滚道半径与滚球半径比值、滚球半径、淬硬层深度和滚道硬度等因素的影响,利用正交试验分析得出其承载能力最强的参数组合,从而提高回转支承的承载能力。
2.技术方案
为达到上述目的,本发明提供的技术方案为:
本发明的一种提高回转支承高温环境下承载能力的方法,利用有限元建模模拟高温环境下回转支承的工况,求解出滚动体和滚道接触区域的应力分布情况,并分析单一结构因素对回转支承承载能力的影响,再利用正交试验分析多个结构因素共同改变时对回转支承承载能力的影响,得出回转支承承载能力最强的参数组合,从而提高回转支承的承载能力。
更进一步地,所述的结构因素包括回转支承接触角、滚道半径与滚球半径比值、滚球半径、淬硬层深度和滚道硬度。
更进一步地,所述的利用有限元建模模拟高温环境下回转支承的工况,求解出滚动体和滚道接触区域的应力分布情况的具体过程为:
1)根据回转支承实际尺寸参数在SolidWorks软件中建立回转支承模型;
2)将SolidWorks软件创建好的回转支承模型导入到ABAQUS有限元软件中;
3)在ABAQUS软件中建立两个静力分析步,在第一个分析步中,对回转支承施加轴向载荷,使滚动体和滚圈之间建立稳定的接触关系;在第二个分析步中施加温度场;
4)在ABAQUS软件的接触设置里选用表面与表面接触类型,且设置滚动体表面为主接触面,内、外圈滚道面为从接触面,同时设置摩擦属性切向行为为“罚”函数,法向行为为“硬接触”;
5)对回转支承内圈施加全约束,对外圈两侧面施加循环对称约束边界,滚动体释放其沿回转支承轴向的自由度,在外圈上端面施加轴向载荷;
6)对滚动体和内外滚圈接触区域进行局部网格细化,网格选用六面体网格,单元类型选择C3D8R;
7)在ABAQUS软件的作业模块创建回转支承分析作业,提交计算,得到滚动体和滚道接触区域的应力分布分析结果。
更进一步地,步骤1)建立回转支承模型过程中,只创建单个滚动体和滚道的接触模型,并在分析时对相关截面添加循环对称边界约束。
更进一步地,步骤1)建立回转支承模型过程中,仅创建回转支承整体模型的1/84进行分析,且在建模时忽略倒角棱边、油孔,螺栓安装孔类对分析影响不大的特征。
更进一步地,分析单一结构因素对回转支承承载能力的影响时,控制其他结构因素的数值固定,对该结构因素设置多组数值,并分别建立常温及高温条件下,各数值对应的回转支承模型,进行仿真分析滚动体和滚道接触区域的应力分布。
更进一步地,在轴向载荷为3168kN,工作环境温度为300℃的工况下,对单排四点接触球式回转支承进行单一结构因素分析;
当分析滚球和滚道接触角对回转支承承载能力的影响时,控制其他结构因素的数值固定,令滚球和滚道接触角分别为30°、35°、40°、45°、50°、55°和60°;
分析滚道半径和滚球半径比对回转支承承载能力的影响时,控制其他结构因素的数值固定,令滚道半径和滚球半径比分别为1.01、1.02、1.03、1.04、1.05、1.06、1.07和1.08;
分析滚球半径对回转支承承载能力的影响时,控制其他结构因素的数值固定,令滚球半径分别为22.5mm、23mm、23.5mm和24mm;
分析淬硬层深度对回转支承承载能力的影响时,控制其他结构因素的数值固定,令淬硬层深度分别为5mm、5.5mm、6mm;
分析滚道硬度对回转支承承载能力的影响时,控制其他结构因素的数值固定,令滚道硬度HRC为55、57、59、61。
更进一步地,利用正交试验分析多个结构因素共同改变对回转支承承载能力的影响时,对每个结构因素选择4个水平,分别为:滚球和滚道接触角30°、40°、50°和60°;滚道半径和滚球半径比1.02、1.04、1.06和1.08;滚球半径22.5mm、23mm、23.5mm和24mm;淬硬层深度5mm、5.5mm、6mm、6.5mm;滚道硬度55、57、59、61。
本发明的一种高温环境下承载能力强的回转支承,在300℃温度环境下,回转支承接触角为60°、滚道半径与滚球半径比值为1.04、滚球半径为22.5mm、淬硬层深度为5.5mm、滚道硬度为61HRC。
3.有益效果
采用本发明提供的技术方案,与已有的公知技术相比,具有如下显著效果:
(1)本发明的一种提高回转支承高温环境下承载能力的方法,针对高温环境下回转支承接触应力大的情况,从回转支承的接触角、滚道半径与滚球半径比值、滚球半径、淬硬层深度和滚道硬度等因素考虑,分析各单一因素对其承载能力的影响,最后通过采用正交试验分析得出多个因素共同改变时回转支承承载能力最强的参数组合,从而提高了回转支承的承载能力,该方案大大地减少了工作量,更加高效、快速、经济;
(2)本发明的一种提高回转支承高温环境下承载能力的方法,利用ABAQUS有限元软件模拟高温回转支承的工况,求解出接触区域的应力分布情况,可以模拟更接近实际的工况,计算结果更加精确;
(3)本发明的一种提高回转支承高温环境下承载能力的方法,在建立回转支承模型时,只创建单个滚动体和滚道的接触模型,并在分析时对相关截面添加循环对称边界约束,同时只对回转支承整体模型的1/84进行分析,大大简化了计算,节约了时间。
附图说明
图1为回转支承的结构示意图;
图2为某型回转支承的结构参数示意图;
图3为回转支承的SolidWorks建模图;
图4为回转支承1/84简化模型图;
图5为载荷和边界条件设置示意图;
图6为回转支承网格划分图;
图7为300℃环境下某型回转支承内圈应力分析结果图;
图8为300℃环境下某型回转支承外圈应力分析结果图;
图9中的(a)和(b)分别为常温、高温下外圈接触应力随接触角的变化图;
图10中的(a)和(b)为常温、高温下外圈接触应力随半径比的变化图;
图11中的(a)和(b)为常温、高温下外圈接触应力随滚球半径的变化图;
图12为各因素不同水平的试验指标值结果图;
图13为结构改进后回转支承内外圈上的接触应力图。
示意图中的标号说明:
1、外圈;2、密封带;3、滚动体;4、加油嘴;5、内圈。
具体实施方式
为进一步了解本发明的内容,结合附图和实施例对本发明作详细描述。
本发明的一种提高回转支承高温环境下承载能力的方法,其核心思想为:在SolidWorks软件中建立回转支承模型并导入ABAQUS软件中,利用ABAQUS有限元软件仿真分析高温回转支承在一定轴向载荷下的承载能力,分析求解回转支承在高温下的接触应力。由于高温下的回转支承接触应力较大,且经发明人长期实践总结得出,回转支承的失效形式主要是滚道受到损坏和断齿,其中滚道损坏的比例达到98%以上,滚道质量的好坏直接影响回转支承的性能。而回转支承承载能力主要与其滚道结构的刚性和滚道区域的应力分布有关,不同的滚道结构参数所对应的应力分布不一样。因此,发明人从回转支承的接触角、滚道半径与滚球半径比值、滚球半径、淬硬层深度和滚道硬度等因素考虑,又由于多因素热力耦合分析计算量太大,为了简化分析过程和节省计算时间,分析各单一因素对其承载能力的影响,最后通过采用正交试验分析法综合考虑多个因素共同改变时对回转支承承载能力的影响,得出其承载能力最强的参数组合和应满足的热处理工艺条件组合,从而提高回转支承的承载能力。本发明可以模拟更接近实际的工况,计算结果更加精确。
实施例1
参看图1,本实施例针对某型单排四点接触球式回转支承进行研究,其结构主要由外圈1、密封带2、滚动体3、加油嘴4和内圈5等组成。其主要结构参数如图2所示,研究工况为:轴向载荷为3168kN,工作环境温度为300℃。
本实施例的具体实施步骤如下:
1)根据图2提供的回转支承尺寸参数在SolidWorks软件中建立如图3所示回转支承的模型。
由于回转支承是对称形状,当其承受轴向载荷时各个滚动体与滚道接触区域所受的载荷均相同,所以从简化模型和方便分析计算的角度,只创建单个滚动体和滚道的接触模型,并在分析时对相关截面添加循环对称边界约束。
利用SolidWorks三维软件创建回转支承整体模型的1/84进行分析,为了节省计算时间,在建模时忽略倒角棱边、油孔,螺栓安装孔等对分析影响不大的特征。创建的1/84回转支承模型如图4所示。
2)将SolidWorks软件创建好的回转支承模型导入到ABAQUS中,通过属性模块创建并定义回转支承的材料属性,并指派到相应的部件实例。
3)选用静力分析步,由于滚动体和滚圈间的复杂接触关系,同时施加轴向载荷和温度场的情况,需建立两个分析步,在第一个分析步中,对回转支承施加轴向载荷,使滚动体和滚圈之间建立稳定的接触关系;在第二个分析步中施加温度场。
4)如图1所示,滚动体与内、外滚圈通过四点接触形成两对接触,当其受到轴向载荷等作用时会由点接触转化为面面接触,因此在ABAQUS接触设置里选用表面与表面接触类型,且设置滚动体表面为主接触面,内、外圈滚道面为从接触面。同时设置摩擦属性切向行为为“罚”函数,设定摩擦系数为0.05,法向行为为“硬接触”。
5)在轴向载荷作用在回转支承的过程中,其内圈是固定不动的,载荷施加在外圈上端面通过滚动体与滚道的接触传递至内圈,因此对回转支承内圈施加全约束,对外圈两侧面施加循环对称约束边界,滚动体释放其沿回转支承轴向的自由度。在外圈上端面施加轴向载荷,大小为Fa=3168kN,换算到1/84模型上的压强为8.8Mpa,具体如图5所示。
6)对滚动体和内外滚圈接触区域进行局部网格细化,以保证计算的准备可靠性,网格选用六面体网格,单元类型选择C3D8R,并对网格进行检查,确保没有错误和警告,划分好的网格如图6所示。
7)在作业模块创建回转支承分析作业,提交计算,得到分析结果如图7、图8所示,300℃环境下其最大S,Mises应力和接触应力分别为2376MPa、4109MPa。
8)由上面的分析结果可以看出,回转支承在高温环境下的接触应力很大。而影响回转支承承载能力的主要因素为其接触角、滚道半径与滚球半径比值、滚球半径、淬硬层深度和滚道硬度等,所以为了对现有回转支承进行优化改进,分别从这5个因素考虑,分析各单一因素对回转支承承载能力的影响,具体如①~⑤所述。
①滚球和滚道接触角对其承载能力的影响
当滚道半径和滚球半径的数值固定时,选择不同大小的滚球和滚道接触角也可能对回转支承承载性能产生影响。以该型单排四点接触球式回转支承为例,滚道半径为10.4mm,滚球半径为22.5mm保持不变,令滚球和滚道的接触角分别为30°、35°、40°、45°、50°、55°、60°建立7组模型,分别进行仿真分析。得出回转支承内外圈上的接触应力如表1所示:
表1不同接触角下内外圈上的接触应力
由表1及图9可以看出:在常温下,接触角在45°~60°范围时,外圈接触应力比较小,因此通常在常温下接触角一般选用45°,但在高温下接触角越大,接触应力越小,而且变化比较明显,因此仅从接触角考虑,如果采用60°的接触角会明显减小接触应力的值,内外圈上的接触应力降低19%左右,降低幅度比较明显。
②滚道半径和滚球半径比值大小对其承载能力的影响
保持回转支承的其他参数(滚球半径、接触角为45°)不变,只改变滚道半径,让滚道半径和滚球半径比值分别为1.01、1.02……1.08,相对应的滚道半径为22.725mm、22.95mm、……24.3mm,建立8组模型进行仿真分析,分析结果如下表:
表2不同半径比下内外圈上的接触应力
由表2和图10可以看出常温下,当滚道半径和滚球半径的比值从1.01逐渐增大到1.08时,回转支承与滚道的距离逐渐变大,能够承受应力的滚道面积越来越小,从而导致滚道上的接触应力越来越大,比值为1.01时滚圈上的接触应力最小,但通常从加工方面、额定静容量和寿命方面考虑通常滚道半径和滚球半径比值取1.04。
300℃高温下,当滚道半径和滚球半径的比值从1.01逐渐增大时,回转支承与滚道的距离虽然逐渐变大,但是变化程度比较小,导致接触应力先减小后基本保持不变。因此,在高温下,回转支承适合采用与常温下相同的滚道半径和滚球半径比值比较合理。
③当滚道半径和滚球半径取1.04时,滚球半径对承载能力的影响
当滚道半径和滚球半径的比值取1.04时,选用不同大小的滚球半径对回转支承承载性能产生影响。
表3不同滚球半径下内外圈上的承接触应力
由表3和图11可以看出,当滚道半径和滚球半径固定时,随着滚球半径的增大,无论是在常温条件还是在高温环境下,内外圈上的接触应力都会增大。增大滚球半径,对降低滚球和滚道之间的接触应力有着重要作用,但是如果为了满足更大的滚球直径而增大回转支承滚圈径向深度和轴向高度,会使回转支承整体变大。而在此处只是增大滚球与滚道的直径,并没有增加滚圈径向深度和轴向高度,因此导致接触应力反而增大,因此在该结构下采用直径为45mm滚球最为适宜。
④淬硬层深度对其承载能力的影响
在回转支承结构中,一定的淬硬层深度是保证回转支承滚道不产生剥落的重要因素,当外载荷作用于回转支承时,滚球与滚道的点对点接触就变成了面对面接触,接触面为椭圆形,滚道除了受到压应力外,还会受到剪切应力的作用。因此淬硬层深度对回转支承的承载能力有着重要影响,需探究滚道不同淬硬层深度对回转支承承载能力的影响。
分别建立3组模型,让内外滚圈淬硬层深度分别为5、5.5、6mm,进行仿真分析,结果如表4所示:
表4不同淬硬层深度下内外圈上的接触应力
通过表4可以看出,随着内外圈滚道淬硬层深度的增加,滚道上的接触应力逐渐减小,但是对于承受一般负荷的回转支承,其滚道淬硬层深度可按支承滚道直径大小确定,大致规定如下:直径小于1500mm,滚道淬硬层深度为3~4mm;直径在1500~5000mm之间时,淬硬层深度为3~5mm;直径大于5000mm,淬硬层深度5~6mm。
⑤滚道硬度的影响
相关标准规定回转支承的最低硬度为HRC55,实际加工中淬火硬度能够达到HRC57左右,而当硬度低于HRC53时,即使留有1.2的安全系数,在使用时也已经不安全了。尤其是当硬度只有HRC50时,即使1.7倍的安全系数也非常危险。因此,硬度不够极易造成回转支承失效。下面针对HRC为55、57、59、61四个硬度值,进行仿真分析,结果见表5:
表5不同滚道硬度下内外圈上的接触应力
9)以上只考虑了单个因素的变化对回转支承承载能力的影响,但不确定多个因素同时改变时,回转支承的承载能力会发生什么变化。而对高温回转支承的优化改进会考虑多个因素,需要对很多组不同模型进行热力耦合分析,将花费大量的时间精力去求解计算。正交试验是一种高效,快速的多因素试验方法,能够同时对多因素、多水平进行比较,通过试验常能找出最佳的结构参数,因此采用正交试验分析法综合考虑多个因素共同改变时对回转支承承载能力的影响,以找到一组最优的结构参数值,使得在高温环境下回转支承滚道上的接触应力最小,承载能力最强。
影响回转支承结构的因素:(i)滚球和滚道接触角;(ii)滚道半径和滚球半径比值;(iii)滚球半径;(iv)淬硬层深度;(v)滚道硬度。在这里考虑这五个因素。每个因素选取4个水平,如表6所示:
表6影响因素及各因素水平
试验指标:回转支承滚道上接触应力的大小,要求接触应力越小越好。
本实施例为四水平五因素试验,从尽量减少试验工作量角度出发,选择正交表L16(45)。
表7正交表L16(45)
A | B | C | D | E | |
实验1 | 1 | 1 | 1 | 1 | 1 |
实验2 | 1 | 2 | 2 | 2 | 2 |
实验3 | 1 | 3 | 3 | 3 | 3 |
实验4 | 1 | 4 | 4 | 4 | 4 |
实验5 | 2 | 2 | 1 | 3 | 4 |
实验6 | 2 | 1 | 2 | 4 | 3 |
实验7 | 2 | 4 | 3 | 1 | 2 |
实验8 | 2 | 3 | 4 | 2 | 1 |
实验9 | 3 | 3 | 1 | 4 | 2 |
实验10 | 3 | 4 | 2 | 3 | 1 |
实验11 | 3 | 1 | 3 | 2 | 4 |
实验12 | 3 | 2 | 4 | 1 | 3 |
实验13 | 4 | 4 | 1 | 2 | 3 |
实验14 | 4 | 3 | 2 | 1 | 4 |
实验15 | 4 | 2 | 3 | 4 | 1 |
实验16 | 4 | 1 | 4 | 3 | 2 |
①排表头
分别把滚球和滚道接触角、滚道半径和滚球半径比值、滚球半径,淬硬层深度和滚道硬度放在正交表L16(45)上方的五个纵列顶端上,如表8所示,每列上放一种,从左到右依次放入接触角A、滚道半径和滚球半径比值B及滚球半径C,淬硬层深度D和滚道硬度E,最后这五个因素相应的水平填入表格内,表中各列的数字表示这些因素的水平,完成后如表8所示:
表8试验方案表
所在列 | A | B | C | D | E | |
因素 | 接触角 | R/比值 | 滚球半径 | 淬硬层深度 | 硬度 | 实验结果 |
实验1 | 30 | 1.02 | 22.5 | 5 | 55 | |
实验2 | 30 | 1.04 | 23 | 5.5 | 57 | |
实验3 | 30 | 1.06 | 23.5 | 6 | 59 | |
实验4 | 30 | 1.08 | 24 | 6.5 | 61 | |
实验5 | 40 | 1.02 | 23 | 6 | 61 | |
实验6 | 40 | 1.04 | 22.5 | 6.5 | 59 | |
实验7 | 40 | 1.06 | 24 | 5 | 57 | |
实验8 | 40 | 1.08 | 23.5 | 5.5 | 55 | |
实验9 | 50 | 1.02 | 23.5 | 6.5 | 57 | |
实验10 | 50 | 1.04 | 24 | 6 | 55 | |
实验11 | 50 | 1.06 | 22.5 | 5.5 | 61 | |
实验12 | 50 | 1.08 | 23 | 5 | 59 | |
实验13 | 60 | 1.02 | 24 | 5.5 | 59 | |
实验14 | 60 | 1.04 | 23.5 | 5 | 61 | |
实验15 | 60 | 1.06 | 23 | 6.5 | 55 | |
实验16 | 60 | 1.08 | 22.5 | 6 | 57 |
②采用直观分析法来分析试验:
在ABAQUS中建立模型进行分析,得到16个外圈滚道接触应力的值填入L16(45)正交表的最右侧一列内,如表9所示。现从有限元分析得到的这16个接触应力值的结果出发,来分析试验结果。
表9承载性能试验结果表
所在列 | A | B | C | D | E | |
因素 | 接触角 | R/比值 | 滚球半径 | 淬硬层深度 | 硬度 | 实验结果 |
实验1 | 1 | 1 | 1 | 1 | 1 | 5029 |
实验2 | 1 | 2 | 2 | 2 | 2 | 5231 |
实验3 | 1 | 3 | 3 | 3 | 3 | 5387 |
实验4 | 1 | 4 | 4 | 4 | 4 | 5423 |
实验5 | 2 | 1 | 2 | 3 | 4 | 4846 |
实验6 | 2 | 2 | 1 | 4 | 3 | 4532 |
实验7 | 2 | 3 | 4 | 1 | 2 | 4787 |
实验8 | 2 | 4 | 3 | 2 | 1 | 4838 |
实验9 | 3 | 1 | 3 | 4 | 2 | 4123 |
实验10 | 3 | 2 | 4 | 3 | 1 | 4098 |
实验11 | 3 | 3 | 1 | 2 | 4 | 3602 |
实验12 | 3 | 4 | 2 | 1 | 3 | 4167 |
实验13 | 4 | 1 | 4 | 2 | 3 | 3426 |
实验14 | 4 | 2 | 3 | 1 | 4 | 3212 |
实验15 | 4 | 3 | 2 | 4 | 1 | 3341 |
实验16 | 4 | 4 | 1 | 3 | 2 | 3579 |
均值1 | 5267.500 | 4356.000 | 4185.500 | 4298.750 | 4326.500 | |
均值2 | 4750.750 | 4268.250 | 4396.250 | 4274.250 | 4430.000 | |
均值3 | 3997.500 | 4279.250 | 4390.000 | 4477.500 | 4378.000 | |
均值4 | 3389.500 | 4501.750 | 4433.500 | 4354.750 | 4270.750 | |
极差 | 1878.000 | 233.500 | 248.000 | 203.250 | 159.250 |
③试验结果的分析
(i)各因素对指标的影响
极差表征的是A、B、C、D、E五个因素对试验指标外圈接触应力值的影响程度。其越大,表征该因素对外圈接触应力值的影响程度也越大。
第1列的极差值最大,说明因素A接触角对回转支承外圈滚道上接触应力的影响最大;因此可以由极差值大小排列顺序得出五个因素的主次顺序:主→次:A→C→B→D→E。
(ii)各因素水平的选取
在该试验中,要求试验指标回转支承外圈滚道上的接触应力越小越好,因此取各因素A、B、C、D、E中四个水平中最小的那个即A4,B2,C1,D2,E4。
(iii)随各因素改变指标的变化规律
将试验中各因素不同水平的试验指标值画在图12中。
10)通过上述研究可以看出,在300℃的高温环境下,为了提高回转支承的承载能力,可以把回转支承的接触角改为60°,同时滚道与滚球半径比仍采用1.04保持不变,淬硬层深度为5.5mm,滚道硬度61HRC,此时滚动体直径为45mm。再以优化后的参数重新建模,分析回转支承在300℃高温环境下的承载情况。优化后回转支承内、外圈滚圈上的接触应力如图13所示,可以看出回转支承结构和热处理工艺条件改进后,在300℃的高温环境下,内、外圈上的接触应力明显减小,内圈上的接触应力相比改进前降低了约900MPa,大约降幅27%左右,外圈上的接触应力降低了约1000MPa,降幅24%左右。因此回转支承在工作温度300℃左右时,可以采用此结构来降低内外圈滚道上的接触应力,提高承载能力,结果见图13。
以上示意性的对本发明及其实施方式进行了描述,该描述没有限制性,附图中所示的也只是本发明的实施方式之一,实际的结构并不局限于此。所以,如果本领域的普通技术人员受其启示,在不脱离本发明创造宗旨的情况下,不经创造性的设计出与该技术方案相似的结构方式及实施例,均应属于本发明的保护范围。
Claims (7)
1.一种提高回转支承高温环境下承载能力的方法,其特征在于:利用有限元建模模拟高温环境下回转支承的工况,求解出滚动体和滚道接触区域的应力分布情况,并分析单一结构因素对回转支承承载能力的影响,再利用正交试验分析多个结构因素共同改变时对回转支承承载能力的影响,得出回转支承承载能力最强的参数组合,从而提高回转支承的承载能力;
所述的利用有限元建模模拟高温环境下回转支承的工况,求解出滚动体和滚道接触区域的应力分布情况的具体过程为:
1)根据回转支承实际尺寸参数在SolidWorks软件中建立回转支承模型;
2)将SolidWorks软件创建好的回转支承模型导入到ABAQUS有限元软件中;
3)在ABAQUS软件中建立两个静力分析步,在第一个分析步中,对回转支承施加轴向载荷,使滚动体和滚圈之间建立稳定的接触关系;在第二个分析步中施加温度场;
4)在ABAQUS软件的接触设置里选用表面与表面接触类型,且设置滚动体表面为主接触面,内、外圈滚道面为从接触面,同时设置摩擦属性切向行为为“罚”函数,法向行为为“硬接触”;
5)对回转支承内圈施加全约束,对外圈两侧面施加循环对称约束边界,滚动体释放其沿回转支承轴向的自由度,在外圈上端面施加轴向载荷;
6)对滚动体和内外滚圈接触区域进行局部网格细化,网格选用六面体网格,单元类型选择C3D8R;
7)在ABAQUS软件的作业模块创建回转支承分析作业,提交计算,得到滚动体和滚道接触区域的应力分布分析结果。
2.根据权利要求1所述的一种提高回转支承高温环境下承载能力的方法,其特征在于:所述的结构因素包括回转支承接触角、滚道半径与滚球半径比值、滚球半径、淬硬层深度和滚道硬度。
3.根据权利要求2所述的一种提高回转支承高温环境下承载能力的方法,其特征在于:步骤1)建立回转支承模型过程中,只创建单个滚动体和滚道的接触模型,并在分析时对相关截面添加循环对称边界约束。
4.根据权利要求3所述的一种提高回转支承高温环境下承载能力的方法,其特征在于:步骤1)建立回转支承模型过程中,仅创建回转支承整体模型的1/84进行分析,且在建模时忽略倒角棱边、油孔,螺栓安装孔类对分析影响不大的特征。
5.根据权利要求4所述的一种提高回转支承高温环境下承载能力的方法,其特征在于:分析单一结构因素对回转支承承载能力的影响时,控制其他结构因素的数值固定,对该结构因素设置多组数值,并分别建立常温及高温条件下,各数值对应的回转支承模型,进行仿真分析滚动体和滚道接触区域的应力分布。
6.根据权利要求5所述的一种提高回转支承高温环境下承载能力的方法,其特征在于:在轴向载荷为3168kN,工作环境温度为300℃的工况下,对单排四点接触球式回转支承进行单一结构因素分析;
当分析滚球和滚道接触角对回转支承承载能力的影响时,控制其他结构因素的数值固定,令滚球和滚道接触角分别为30°、35°、40°、45°、50°、55°和60°;
分析滚道半径和滚球半径比对回转支承承载能力的影响时,控制其他结构因素的数值固定,令滚道半径和滚球半径比分别为1.01、1.02、1.03、1.04、1.05、1.06、1.07和1.08;
分析滚球半径对回转支承承载能力的影响时,控制其他结构因素的数值固定,令滚球半径分别为22.5mm、23mm、23.5mm和24mm;
分析淬硬层深度对回转支承承载能力的影响时,控制其他结构因素的数值固定,令淬硬层深度分别为5mm、5.5mm、6mm;
分析滚道硬度对回转支承承载能力的影响时,控制其他结构因素的数值固定,令滚道硬度为55、57、59、61。
7.根据权利要求6所述的一种提高回转支承高温环境下承载能力的方法,其特征在于:利用正交试验分析多个结构因素共同改变对回转支承承载能力的影响时,对每个结构因素选择4个水平,分别为:滚球和滚道接触角30°、40°、50°和60°;滚道半径和滚球半径比1.02、1.04、1.06和1.08;滚球半径22.5mm、23mm、23.5mm和24mm;淬硬层深度5mm、5.5mm、6mm、6.5mm;滚道硬度55、57、59、61。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610625931.XA CN106294957B (zh) | 2016-08-02 | 2016-08-02 | 一种提高回转支承高温环境下承载能力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610625931.XA CN106294957B (zh) | 2016-08-02 | 2016-08-02 | 一种提高回转支承高温环境下承载能力的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106294957A CN106294957A (zh) | 2017-01-04 |
CN106294957B true CN106294957B (zh) | 2019-04-23 |
Family
ID=57664035
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610625931.XA Expired - Fee Related CN106294957B (zh) | 2016-08-02 | 2016-08-02 | 一种提高回转支承高温环境下承载能力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106294957B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110516290B (zh) * | 2019-07-12 | 2022-10-21 | 徐州美驰车桥有限公司 | 一种球回转支承的有限元分析方法 |
CN113326563A (zh) * | 2021-05-24 | 2021-08-31 | 扬州大学 | 一种基于反馈筛选的智能拨穗机最佳工作参数确定方法 |
CN114091306B (zh) * | 2021-11-18 | 2024-04-26 | 大连华锐重工集团股份有限公司 | 基于有限元模拟的回转支撑轴承安装平面挠曲量分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104235177A (zh) * | 2012-09-12 | 2014-12-24 | 马鞍山统力回转支承有限公司 | 渐变硬度齿轮式回转支承齿轮淬火加工工艺 |
FR3025805A1 (fr) * | 2014-09-12 | 2016-03-18 | Emag Holding Gmbh | Procede de trempe de pieces de grande taille |
-
2016
- 2016-08-02 CN CN201610625931.XA patent/CN106294957B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104235177A (zh) * | 2012-09-12 | 2014-12-24 | 马鞍山统力回转支承有限公司 | 渐变硬度齿轮式回转支承齿轮淬火加工工艺 |
FR3025805A1 (fr) * | 2014-09-12 | 2016-03-18 | Emag Holding Gmbh | Procede de trempe de pieces de grande taille |
Non-Patent Citations (3)
Title |
---|
冷拉棒材锻造性能工艺研究;王永丰;《解决方案》;20111231(第8期);164-165 |
单排四点接触球式回转支承承载能力分析;李庚袁;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20150715(第07期);C029-90 |
回转支承承载性能分析与优化;郑岩;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20150615(第06期);C029-35 |
Also Published As
Publication number | Publication date |
---|---|
CN106294957A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106294957B (zh) | 一种提高回转支承高温环境下承载能力的方法 | |
CN105117536B (zh) | 一种rpv含裂纹类缺陷的简化弹塑性断裂力学分析方法 | |
CN103970965A (zh) | 燃气涡轮发动机加速寿命试验试车方法 | |
He et al. | Hardened raceway calculation analysis of a three-row roller slewing bearing | |
Ju et al. | Fatigue design of offshore wind turbine jacket-type structures using a parallel scheme | |
He et al. | Optimization design of structural parameters of single-row four-point contact ball slewing bearing | |
CN111985141A (zh) | 一种转盘轴承表面硬化滚道许用接触应力的确定方法 | |
QIAN et al. | FE analysis for radial spread behavior in three-roll cross rolling with small-hole and deep-groove ring | |
Wang et al. | A three-dimensional constitutive model of high strength structural steels | |
CN109948215B (zh) | 一种热冲压工艺制定方法 | |
CN104732020A (zh) | 千吨级全地面起重机用六排滚子转盘轴承的校核方法 | |
CN103148095B (zh) | 全陶瓷轴承及其制造方法、基本结构参数的确定方法 | |
Wang et al. | Alternative approach for imperfection-tolerant design optimization of stiffened cylindrical shells via energy barrier method | |
Liu et al. | Multiaxial fatigue life prediction model with relative stress gradient and size effect | |
CN106768547A (zh) | 一种计算增生装置滚轮接触力的方法 | |
Fang et al. | Effect of weld on axial buckling of cylindrical shells | |
CN107092726B (zh) | 一种提高飞轮转子储能密度的结构设计方法及装置 | |
Zhu et al. | Development of vacuum vessel design and analysis module for CFETR integration design platform | |
CN109915327B (zh) | 确定钢筋混凝土-钢组合风机塔架固有频率的方法 | |
Hu et al. | Analysis on effect of injection residual stress on impact resistance of composite wheel made of long glass fiber reinforced thermoplastic | |
Huang et al. | A simple method for estimating the roundness of minimum zone circle | |
Xu et al. | Reliability analysis of tractor multi‐way valves based on the improved weighted grey relational method | |
CN109959555B (zh) | 考虑厚度和裂纹深度影响的中心裂纹拉伸试样的孕育期预测方法 | |
CN205898548U (zh) | 具有温度变形自补偿功能的混凝土温度应力试验机 | |
Song | Research on Mechanical Design and Manufacturing of Rationalization of Machining Manufacturing Process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190423 |