CN109241695B - 基于新颖pso-sade的滚动轴承可靠性多目标优化设计方法 - Google Patents

基于新颖pso-sade的滚动轴承可靠性多目标优化设计方法 Download PDF

Info

Publication number
CN109241695B
CN109241695B CN201811365759.4A CN201811365759A CN109241695B CN 109241695 B CN109241695 B CN 109241695B CN 201811365759 A CN201811365759 A CN 201811365759A CN 109241695 B CN109241695 B CN 109241695B
Authority
CN
China
Prior art keywords
bearing
pso
sade
rolling bearing
optimization design
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.)
Active
Application number
CN201811365759.4A
Other languages
English (en)
Other versions
CN109241695A (zh
Inventor
张义民
高淑芝
徐林涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenyang University of Chemical Technology
Original Assignee
Shenyang University of Chemical Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenyang University of Chemical Technology filed Critical Shenyang University of Chemical Technology
Priority to CN201811365759.4A priority Critical patent/CN109241695B/zh
Publication of CN109241695A publication Critical patent/CN109241695A/zh
Application granted granted Critical
Publication of CN109241695B publication Critical patent/CN109241695B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]

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)
  • Biophysics (AREA)
  • Artificial Intelligence (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (AREA)
  • Biomedical Technology (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Rolling Contact Bearings (AREA)

Abstract

基于新颖PSO‑SADE的滚动轴承可靠性多目标优化设计方法，涉及一种滚动轴承优化设计方法，该方法为一种新颖PSO‑SADE融合算法对轮毂轴承进行多目标可靠性优化设计；以滚动轴承最大疲劳寿命、弹性动力润滑膜厚度、摩擦功率损耗作为目标函数，建立相应的非线性动态约束；其次，提出新颖PSO‑SADE融合算法来解决滚动轴承可靠性多目标优化设计问题；最后对滚动轴承的可靠性多目标优化设计进行了仿真研究；结果表明，滚动轴承多目标优化设计后的寿命相对于标准轴承寿命要更长，更能适用于多种变载荷的运行环境。这也为轴承在更多变载荷条件下进行优化设计提供了一定基础和依据。

Description

基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法
技术领域
本发明涉及一种滚动轴承优化设计方法，特别是涉及一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法。
背景技术
滚动轴承这一机械元件，现已成为国内外大型工业旋转机械设备的关键性零件之一，同时在航空航天、车辆、家用电器、医疗设备等方面也充当着至关重要的角色。所以滚动轴承的优化设计无论是在寿命方面还是可靠性方面都应该去实现其最大化的目的。滚动轴承的内部几何结构较为复杂，同时其内部结构几何参数对于轴承本身的寿命、润滑度、可靠度、摩擦等特性有重要的影响，最终使轴承能够正常、稳定、持续运行，实现长寿命的目的。
发明内容
本发明的目的在于提供一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，本发明方法采用粒子群算法和自适应差分算法进行融合的方法来实现滚动轴承的可靠性多目标优化设计。首先通过粒子群优化算法对滚动轴承进行可靠性多目标优化设计，然后利用SADE对PSO经过迭代后产生的精英集进一步利用使进入下一个种群迭代，并达到最终的目的。通过此次的算法融合，滚动轴承优化设计后的效果非常显著，在针对解决汽车轮毂轴承承受外部载荷能力差的问题上提供了方法依据。
本发明的目的是通过以下技术方案实现的：
一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，所述方法包括以下步骤：
步骤一：首先明确滚动轴承的基本构造、工作环境，以及实际影响轴承疲劳寿命的设计参数和新颖PSO-SADE算法参数；
步骤二：现建立三种深沟球轴承优化设计模型（目标函数）
包括： 2.1 基本额定动载荷
Figure 443778DEST_PATH_IMAGE001
；
2.2 弹性动力润滑最小膜厚度
Figure 308966DEST_PATH_IMAGE002
；
2.3 摩擦功率损耗
Figure 323190DEST_PATH_IMAGE003
；
2.4 约束条件
步骤三：应用提出的新颖PSO-SADE融合算法对深沟球轴承进行优化
3.1融合算法（PSO-SADE）思想:首先初始数据集时，利用其中的PSO算法做一个整体迭代寻优；将最优个体汇聚成一个数据集，然后利用该数据集和SADE结合，求出最优解；
3.2 3.2 算法步骤如下：
（1） 首先初始化种群，以及各个参数；
（2）将其内部分成两个不同的群体，所以为了避免在一个区域内引起种群之间的相互交流，现将粒子群算法位置更新步骤中包含的社会学习部分去除，即c2学习因子去掉，保留c1；
（3）先计算PSO适应度值，然后将迭代完之后的所有最优个体组成一个最优集；
（4）利用该集合和自适应差分算法的变异、交叉、选择进行集中操作；
（5）最后将每一次迭代完之后的最优个体汇聚成一个数据集，求出最优解。
所述的一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，所述设计参数，确定一系列轴承几何、约束参数：
Figure 906618DEST_PATH_IMAGE004
是节圆直径；
Figure 31700DEST_PATH_IMAGE005
是球直径，
Figure 892339DEST_PATH_IMAGE006
是球的数量，
Figure 202098DEST_PATH_IMAGE007
是内滚道槽的曲率半径系数，
Figure 148188DEST_PATH_IMAGE008
是外滚道槽的曲率半径系数；
Figure 201595DEST_PATH_IMAGE009
分别是滚动体的最小、最大可能直径，它们各自决定着
Figure 916741DEST_PATH_IMAGE010
的下界和上界；
Figure 397401DEST_PATH_IMAGE011
是轴承内、外环厚度因子；
Figure 565208DEST_PATH_IMAGE012
是轴承运动约束参数；
Figure 297672DEST_PATH_IMAGE013
是滚动体有效长度的参数；以上皆为设计变量，由
Figure 991959DEST_PATH_IMAGE014
表示：
Figure 518886DEST_PATH_IMAGE015
现确定深沟球轴承的工况条件，即径向力，轴向力，轴承实际转速参数。
所述的一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，所述PSO-SADE算法的关键参数，即种群总个体数N、惯行权重w、学习因子c1、总迭代次数Gmax、变异概率F、交叉概率Cr ；一般情况下，学习因子取1.0；处理多个设计变量问题，种群个体总数取100~200；交叉概率Cr=0.8；变异概率F在此算法中以指数形式逐步递减，初始设为0.5；最大迭代次数Gmax为500；以上参数设置为后续的滚动轴承可靠性多目标优化设计提供数据基础。
所述的一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，所述基本额定动载荷
Figure 564203DEST_PATH_IMAGE016
：首先以滚动轴承基本额定动载荷为目标；基本额定动载荷是最能直接决定轴承寿命大小的量，即额定动载荷越大，轴承承受外部载荷的能力就越强，其寿命也就越长；轴承寿命由等式：
Figure 834778DEST_PATH_IMAGE017
定义，
Figure 258937DEST_PATH_IMAGE018
是可靠度为90%的寿命，也称为基本额定寿命，
Figure 346979DEST_PATH_IMAGE019
（r/min）为轴承转速；
Figure 489379DEST_PATH_IMAGE020
轴承当量动载荷，
Figure 688279DEST_PATH_IMAGE021
是基本额定动载荷；
Figure 966944DEST_PATH_IMAGE022
为球轴承的载荷-寿命指数；
Figure 225887DEST_PATH_IMAGE023
在等式（20）中，
Figure 855583DEST_PATH_IMAGE024
和
Figure 467961DEST_PATH_IMAGE025
分别为径向载荷系数、轴向载荷系数，Fa为轴向载荷，Fr为径向载荷；
为使优化设计后的轴承可靠度能够达到最好，最终的目的就是使基本额定动载荷最大化，由等式（3）可得
Figure 991346DEST_PATH_IMAGE026
Figure 296557DEST_PATH_IMAGE027
——轴承动载荷系数，有
Figure 272603DEST_PATH_IMAGE028
Figure 688672DEST_PATH_IMAGE029
——深沟球轴承的接触角，为0°；
上述（3）式中基本额定载荷相对于滚动体直径的大小，有适用范围。
所述的一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，所述
Figure 207509DEST_PATH_IMAGE030
的存在能够尽量避免轴承实际运行过程中滚道和滚动体接触的面积（点接触），并通过使用优化方法计算找到最大的
Figure 542676DEST_PATH_IMAGE030
，这可以在很大程度上延长轴承的摩损寿命；因此最小膜厚度方程为：
Figure 146963DEST_PATH_IMAGE031
作用在滚动体上最大载荷Q为：
Figure 225778DEST_PATH_IMAGE032
其中
Figure 599122DEST_PATH_IMAGE033
为接触点处x方向上的曲率半径，
Figure 105189DEST_PATH_IMAGE034
为接触点处y方向上的曲率半径，
Figure 196773DEST_PATH_IMAGE035
压力粘度系数，
Figure 954645DEST_PATH_IMAGE036
是等效弹性模量。
所述的一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，所述约束条件包括：找出由轴承摩擦引起的功率损耗；设定一系列非线性约束不等式，找到最合适的轴承数据；得出相应的约束条件。
本发明的优点与效果是：
1.本发明提出一种新颖PSO-SADE融合优化算法并应用于6204深沟球轴承优化设计中，在满足给定的非线性约束条件下，利用其遍历性和非线性，对滚动轴承可靠性优化设计进行研究，提高了算法解决高维度问题的能力，最大限度地提高了滚动轴承的疲劳寿命。
2.本发明采用PSO和SADE两种算法进行结合，一是能够快速迭代寻优，二是计算的结果多样性要比普通的多目标优化算法产生的效果要好。并为后续的领域研究提供了一定依据和基础。
3.本发明算法相对于其他算法的不同点在于：首先初始数据集时，利用其中的PSO算法做一个整体迭代寻优，相当于为下一代的进化操作引入了一个精英集。此集合的产生为SADE算法产生的种群提供了更多的选择，同时在一定程度上减少了运算复杂度。
附图说明
图1 深沟球轴承的二维几何平面图；
图2 深沟球轴承的几何剖视图；
图3 三种算法与标准值的对比曲线；
图4 可靠度对比曲线；
图5 算法步骤如下述流程图。
具体实施方式
下面结合附图所示实施例对本发明进行详细说明。
本发明为一种新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法。
滚动轴承的优化设计变量：
Figure 57861DEST_PATH_IMAGE037
在这些参数中，
Figure 875775DEST_PATH_IMAGE038
是节圆直径，
Figure 313710DEST_PATH_IMAGE039
是球直径，
Figure 633131DEST_PATH_IMAGE040
是球的数量，
Figure 105700DEST_PATH_IMAGE041
是内滚道槽的曲率半径系数，
Figure 828937DEST_PATH_IMAGE042
是外滚道槽的曲率半径系数。
Figure 629534DEST_PATH_IMAGE043
分别是滚动体的最大、最小可能直径，它们各自决定着
Figure 853842DEST_PATH_IMAGE005
的下界和上界。
Figure 56284DEST_PATH_IMAGE044
是轴承内、外环厚度因子。
Figure 340635DEST_PATH_IMAGE045
是轴承运动约束参数。
Figure 628528DEST_PATH_IMAGE046
是滚动体有效长度的参数。其中的每一个设计参数对轴承的疲劳寿命和运行状态都存在一定的影响。如图1和图2分别是对深沟球轴承的几何图形参数做了详细描述。
（1） 基本额定动载荷
Figure 390948DEST_PATH_IMAGE047
：
首先以滚动轴承基本额定动载荷为目标。基本额定动载荷是最能直接决定轴承寿命大小的量，即额定动载荷越大，轴承承受外部载荷的能力就越强，其寿命也就越长。轴承寿命由等式：
Figure 447896DEST_PATH_IMAGE048
定义，
Figure 778515DEST_PATH_IMAGE049
是可靠度为90%的寿命，也称为基本额定寿命，
Figure 678338DEST_PATH_IMAGE050
（r/min）为轴承转速。
Figure 119814DEST_PATH_IMAGE051
此时为当量动载荷，
Figure 31270DEST_PATH_IMAGE052
是基本额定动载荷。
Figure 657423DEST_PATH_IMAGE053
为球轴承的载荷-寿命指数。
Figure 919908DEST_PATH_IMAGE054
在等式（20）中，
Figure 24131DEST_PATH_IMAGE055
和
Figure 100002_DEST_PATH_IMAGE056
分别为径向载荷系数、轴向载荷系数，a为轴向载荷，Fr为径向载荷。
为使优化设计后的轴承可靠度能够达到最好，最终的目的就是使基本额定动载荷最大化，由等式（3）可得
Figure 931038DEST_PATH_IMAGE057
Figure 728093DEST_PATH_IMAGE058
——轴承动载荷系数，有
Figure 212295DEST_PATH_IMAGE059
Figure 385787DEST_PATH_IMAGE060
——深沟球轴承的接触角，为0°。
上述（3）式中基本额定载荷相对于滚动体直径的大小，有适用范围。
（2） 弹性动力润滑最小膜厚度
Figure 271835DEST_PATH_IMAGE061
目前来说，和轴承摩损寿命相关性最密切的属滚动轴承的最小油膜厚度。在轴承规定的标准边界内，
Figure 239791DEST_PATH_IMAGE061
的存在能够尽量避免轴承实际运行过程中滚道和滚动体接触的面积（点接触），并通过使用优化方法计算找到最大的
Figure 211289DEST_PATH_IMAGE061
，这可以在很大程度上延长轴承的摩损寿命。因此最小膜厚度方程为：
Figure 657314DEST_PATH_IMAGE062
作用在滚动体上最大载荷Q为：
Figure 538813DEST_PATH_IMAGE063
其中
Figure 677670DEST_PATH_IMAGE064
为接触点处x方向上的曲率半径，
Figure 136465DEST_PATH_IMAGE065
为接触点处y方向上的曲率半径，
Figure 386180DEST_PATH_IMAGE066
压力粘度系数，
Figure 246820DEST_PATH_IMAGE067
是等效弹性模量，
Figure 431945DEST_PATH_IMAGE068
为当量动载荷。
（3）摩擦功率损耗
Figure 502669DEST_PATH_IMAGE069
在轴承实际工作过程中，其内部的各个零件会直接相互接触并发生相对滑动的现象，所以轴承在实际运行过程中是存在着摩擦阻力，其中任何状态下的摩擦都会造成能量损失，并使轴承的运转性能降低。实际上，运转过程中的滚动轴承因为摩擦会导致温度升高，该损耗可通过摩擦阻力矩来衡量。
总摩擦转矩为
Figure 431442DEST_PATH_IMAGE070
Figure 271222DEST_PATH_IMAGE071
——由外部载荷引起的摩擦转矩经验公式
Figure 627248DEST_PATH_IMAGE072
——润滑剂粘性摩擦产生的力矩经验公式
由轴承摩擦引起的功率损耗方程为
Figure 919689DEST_PATH_IMAGE073
其中P L 的单位是W，M的单位N·mm2。轴承的摩擦功率损耗越低，轴承的使用价值将大大提升。
（4）约束条件
为了满足实际轴承装配可能性和优化算法执行效果，要设定一系列非线性约束不等式，找到最合适的轴承数据。
通过以下公式推导出深沟球轴承实际球的数量和直径范围大小：
Figure 652153DEST_PATH_IMAGE074
最大装配角
Figure 221806DEST_PATH_IMAGE075
大小由等式（19）确定：
Figure 873367DEST_PATH_IMAGE076
而且轴承的装配角和轴承几何关系非常密切，根据上述关系式推导，约束
Figure 794050DEST_PATH_IMAGE077
如下：
Figure 189259DEST_PATH_IMAGE078
轴承滚动体的参数界限也要设定界限，有
Figure 613418DEST_PATH_IMAGE079
Figure 701460DEST_PATH_IMAGE080
和
Figure 843859DEST_PATH_IMAGE081
分别决定着
Figure 777180DEST_PATH_IMAGE082
的最小和最大可能直径，并且
Figure 321425DEST_PATH_IMAGE083
和
Figure 455734DEST_PATH_IMAGE084
分别是轴承外径和内径，从上述不等式中可以得出，相应的约束条件为：
Figure 944485DEST_PATH_IMAGE085
上述约束
Figure 822442DEST_PATH_IMAGE086
和
Figure 345827DEST_PATH_IMAGE087
中的参数
Figure 385458DEST_PATH_IMAGE088
和
Figure 502450DEST_PATH_IMAGE089
，它们都是在一定的区间内随机取值，根据轴承几何参数的确定了范围。
为了能使轴承能够持续稳定运行，轴承节径和平均直径的范围应小于一定的给定值。因此，需要满足以下两个约束条件：
Figure 43153DEST_PATH_IMAGE090
Figure 272973DEST_PATH_IMAGE091
其中
Figure 873719DEST_PATH_IMAGE092
是未知常数，它和轴承的运动状态有关。
在轴承实际运行过程中，内环与滚动体接触产生的弹性形变量较大，所以轴承内环需要施加约束。
Figure 212427DEST_PATH_IMAGE093
如下：
Figure 432187DEST_PATH_IMAGE094
其中
Figure 930165DEST_PATH_IMAGE095
和
Figure 311599DEST_PATH_IMAGE096
分别是内、外滚道凹槽直径。轴承外滚道底部的环厚度应大于等于
Figure 278549DEST_PATH_IMAGE097
，因为此约束对于球直径的影响比较明显，所以有关于球直径和外圈厚度的约束条件
Figure 36420DEST_PATH_IMAGE098
：
Figure 388904DEST_PATH_IMAGE099
为了保证球体在轴承制造时不会被挤出套圈端面，所以球直径不能大于轴承宽度的80%~85%，则有如下约束
Figure 941239DEST_PATH_IMAGE100
：
Figure 379174DEST_PATH_IMAGE101
第9约束是为了使内圈的外径大于内径：
Figure 940737DEST_PATH_IMAGE102
第10约束是为了确保外圈的内径小于外径：
Figure 147727DEST_PATH_IMAGE103
第11个设计约束是指滚动体中产生的最大接触应力，并且将它的大小限制在轴承制造材料所能承受的范围内。对于椭圆接触处的应力由如下等式确定：
Figure 136543DEST_PATH_IMAGE104
Figure 937140DEST_PATH_IMAGE105
是椭形接触处的长半轴，
Figure 161448DEST_PATH_IMAGE106
椭形接触处的短半轴，约束如下：
Figure 363890DEST_PATH_IMAGE107
其中，
Figure 648241DEST_PATH_IMAGE108
是轴承所能允许的最大应力。
在设计轴承时，应注意保持架的实际宽度和钢球直径的关系。约束
Figure 670554DEST_PATH_IMAGE109
如下：
Figure 698553DEST_PATH_IMAGE110
Figure 755502DEST_PATH_IMAGE111
是保持架宽度，
Figure 961487DEST_PATH_IMAGE112
是关于保持架本身的结构材料特性参数。对于200系列球轴承
Figure 861310DEST_PATH_IMAGE113
，且精度为0.1。
内、外滚道的曲率半径和钢球直径的等式关系约束
Figure 302786DEST_PATH_IMAGE114
和
Figure 338875DEST_PATH_IMAGE115
如下：
Figure 840395DEST_PATH_IMAGE116
（5）算法步骤
（1）首先初始化种群，以及各个参数，
（2）将其内部分成两个不同的群体，所以为了避免在一个区域内引起种群之间的相互交流，现将粒子群算法位置更新步骤中包含的社会学习部分去除，即c2学习因子去掉，保留c1。
（3）先计算PSO适应度值，然后将迭代完之后的所有最优个体组成一个最优集，
（4）利用该集合和自适应差分算法的变异、交叉、选择进行集中操作。
（5）最后将每一次迭代完之后的最优个体汇聚成一个数据集，求出最优解。
实施例
本文选取6204型号深沟球轴承，在给定几何约束条件下，将其相关几何参数记录下来，如下：
Figure 837301DEST_PATH_IMAGE117
现在设计在上述提到的三个影响轴承疲劳寿命和摩损寿命的目标函数，即动态承载能力、弹性动力润滑最小膜厚度、摩擦功率损耗。
无论是轴承本身的几何设计参数（节圆直径、球直径等等），还是关于轴承运行的约束常数因子都会被考虑作为轴承优化设计的设计变量。首先，我们先要确定一种类型的轴承，并将它的尺寸参数记录下来，然后在给定设计约束条件下来确定轴承几何参数的范围，最后利用PSO和改进PSO-SADE融合算法进行优化设计，并将得到的结果和标准的轴承动态承载能力分别进行比较，其中两者差值越大的对应算法性能就更优越。
表1中的数据是在优化设计之后对这10个设计变量进行统计，每一个参数集合都是算法在迭代15次之后得到的最优个体，并将其记录在下列表格中。
表1 优化设计参数
Figure 207103DEST_PATH_IMAGE118
从以下表2中看出，以上是三种算法对多种型号轴承进行优化设计之后产生的一系列最优解，然后将最优解和标准值进行比较分析，可以看出PSO-SADE的结果要更优于其它算法迭代后的值。
表2 多种方法下的动态承载能力结果对比
Figure 238644DEST_PATH_IMAGE119
表 3 多目标优化设计结果值
Figure 911065DEST_PATH_IMAGE120
本发明主要是以标准轴承尺寸为轴承内部几何参数进行优化设计，将表2中得到的动态承载能力作为比较对象，进而挑选出最优个体来计算最小膜厚度和摩擦功率损耗值。动态承载能力主要影响着滚动轴承疲劳寿命，弹性动力润滑最小膜厚度的大小影响着轴承的磨损寿命，摩擦功率损耗主要影响着轴承的实际运行状态。最终通过对三者的计算来实现多目标优化设计的目的。
将滚动轴承经PSO-SADE优化后的数据进行统计，然后在MATLAB环境中通过运行相应程序，并将三个目标函数结果做成迭代曲线，进而和已知的GA、PSO结果曲线进行对比分析。由附图2可以看出，随着轴承尺寸（型号）的逐渐递增，对应轴承的动态承载能力逐渐增加，相对应的新颖PSO-SADE优化设计结果要优于其余方法的设计效果。
表 4 设定环境下的轴承可靠度和寿命
Figure 395267DEST_PATH_IMAGE121
随着轴承运行时间的增加，可靠度开始逐渐递减。利用PSO-SADE、PSO和标准值进行仿真得出的可靠度曲线如附图3所示。虽然在本文中该轴承最大寿命设为5000h，但从图3中曲线的走势来看，轴承的失效率也在随着时间的延长而逐渐升高，也就是说轴承的可靠度也在逐渐降低，最终到0。正常情况下，在评估轴承寿命时，可靠度降到50%以下就已经算是非常低的水平了，对于这种轴承只能选择弃用，如果把可靠度
Figure 178546DEST_PATH_IMAGE123
作为寿命评估标准的话，这在很大程度上提高了轴承实际使用价值，对于轴承的选取也很有实际意义。
根据表4可知，通过仿真分析，本发明提出的新颖PSO-SADE计算出的寿命值相对于标准的轴承疲劳寿命值要增加了1.8倍，可靠度提升了11.1倍，而且轴承优化设计后的效果非常显著。在针对解决汽车轮毂轴承承受外部载荷能力差的问题上提供了方法依据。

Claims (6)

1.一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，其特征在于，所述方法包括以下步骤：
步骤一：首先明确滚动轴承的基本构造、工作环境，以及实际影响轴承疲劳寿命的设计参数和新颖PSO-SADE算法参数； 根据实际需要选取新颖PSO-SADE算法的关键参数；
步骤二：现建立三种深沟球轴承优化设计模型（目标函数）
包括： 2.1基本额定动载荷
Figure DEST_PATH_IMAGE002
；
2.2 弹性动力润滑最小膜厚度
Figure DEST_PATH_IMAGE004
；
2.3 摩擦功率损耗
Figure DEST_PATH_IMAGE006
；
2.4 约束条件
步骤三：应用提出的新颖PSO-SADE融合算法对深沟球轴承进行优化
3.1融合算法（PSO-SADE）思想:首先初始数据集时，利用其中的PSO算法做一个整体迭代寻优；将最优个体汇聚成一个数据集，然后利用该数据集和SADE结合，求出最优解；
算法步骤如下：
首先初始化种群，以及各个参数；
（2）将其内部分成两个不同的群体，所以为了避免在一个区域内引起种群之间的相互交流，现将粒子群算法位置更新步骤中包含的社会学习部分去除，即c2学习因子去掉，保留c1；
（3）先计算PSO适应度值，然后将迭代完之后的所有最优个体组成一个最优集；
（4）利用该集合和自适应差分算法的变异、交叉、选择进行集中操作；
（5）最后将每一次迭代完之后的最优个体汇聚成一个数据集，求出最优解。
2.根据权利要求1所述的一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，其特征在于，所述设计参数，确定一系列轴承几何、约束参数：
Figure DEST_PATH_IMAGE008
是节圆直径，
Figure DEST_PATH_IMAGE010
是球直径，
Figure DEST_PATH_IMAGE012
是球的数量，
Figure DEST_PATH_IMAGE014
是内滚道槽的曲率半径系数，
Figure DEST_PATH_IMAGE016
是外滚道槽的曲率半径系数；
Figure DEST_PATH_IMAGE018
分别是滚动体的最大、最小可能直径，它们各自决定着
Figure DEST_PATH_IMAGE020
的下界和上界；
Figure DEST_PATH_IMAGE022
是轴承内、外环厚度因子；
Figure DEST_PATH_IMAGE024
是轴承运动约束参数；
Figure DEST_PATH_IMAGE026
是滚动体有效长度的参数；以上皆为设计变量，由
Figure DEST_PATH_IMAGE028
表示：
Figure DEST_PATH_IMAGE030
现确定深沟球轴承的工况条件，即径向力，轴向力，轴承实际转速参数。
3.根据权利要求1所述的一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，其特征在于，所述PSO-SADE算法的关键参数，即种群总个体数N、惯行权重w、学习因子c1、总迭代次数Gmax、变异概率F、交叉概率Cr ；通常来说，学习因子取1.0；处理多个设计变量问题，种群个体总数取100~200；交叉概率Cr=0.8；变异概率F在此算法中以指数形式逐步递减，初始设为0；最大迭代次数Gmax为500次；以上参数设置为后续的滚动轴承可靠性多目标优化设计提供数据基础。
4.根据权利要求1所述的一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，其特征在于，所述基本额定动载荷
Figure DEST_PATH_IMAGE032
：首先以滚动轴承基本额定动载荷为目标；基本额定动载荷是最能直接决定轴承寿命大小的量，即额定动载荷越大，轴承承受外部载荷的能力就越强，其寿命也就越长；轴承寿命由等式：
Figure DEST_PATH_IMAGE034
定义，
Figure DEST_PATH_IMAGE036
是可靠度为90%的寿命，也称为基本额定寿命，
Figure DEST_PATH_IMAGE038
（r/min）为轴承转速；
Figure DEST_PATH_IMAGE040
此时为当量动载荷，
Figure DEST_PATH_IMAGE042
是基本额定动载荷；
Figure DEST_PATH_IMAGE044
为球轴承的载荷-寿命指数；
Figure DEST_PATH_IMAGE046
在等式（20）中，
Figure DEST_PATH_IMAGE048
分别为径向载荷系数、轴向载荷系数，
Figure DEST_PATH_IMAGE050
为轴向载荷，
Figure DEST_PATH_IMAGE052
为径向载荷；
为使优化设计后的轴承可靠度能够达到最好，最终的目的就是使基本额定动载荷最大化，由等式（3）可得
Figure DEST_PATH_IMAGE054
Figure DEST_PATH_IMAGE056
—轴承动载荷系数，有
Figure DEST_PATH_IMAGE058
Figure DEST_PATH_IMAGE060
—深沟球轴承的接触角，为0°；
上述（3）式中基本额定载荷相对于滚动体直径的大小，有适用范围。
5.根据权利要求1所述的一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，其特征在于，所述
Figure DEST_PATH_IMAGE062
的存在能够尽量避免轴承实际运行过程中滚道和滚动体接触的面积（点接触），并通过使用优化方法计算找到最大的
Figure 45291DEST_PATH_IMAGE062
，这可以在很大程度上延长轴承的摩损寿命；因此最小膜厚度方程为：
Figure DEST_PATH_IMAGE064
作用在滚动体上最大载荷Q为：
Figure DEST_PATH_IMAGE066
其中
Figure DEST_PATH_IMAGE068
为接触点处x方向上的曲率半径，
Figure DEST_PATH_IMAGE070
为接触点处y方向上的曲率半径，
Figure DEST_PATH_IMAGE072
压力粘度系数，
Figure DEST_PATH_IMAGE074
是等效弹性模量。
6.根据权利要求1所述的一种基于新颖PSO-SADE的滚动轴承可靠性多目标优化设计方法，其特征在于，所述约束条件包括：找出由轴承摩擦引起的功率损耗；设定一系列非线性约束不等式，找到最合适的轴承数据；得出相应的约束条件。
CN201811365759.4A 2018-11-16 2018-11-16 基于新颖pso-sade的滚动轴承可靠性多目标优化设计方法 Active CN109241695B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811365759.4A CN109241695B (zh) 2018-11-16 2018-11-16 基于新颖pso-sade的滚动轴承可靠性多目标优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811365759.4A CN109241695B (zh) 2018-11-16 2018-11-16 基于新颖pso-sade的滚动轴承可靠性多目标优化设计方法

Publications (2)

Publication Number Publication Date
CN109241695A CN109241695A (zh) 2019-01-18
CN109241695B true CN109241695B (zh) 2022-10-28

Family

ID=65075590

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811365759.4A Active CN109241695B (zh) 2018-11-16 2018-11-16 基于新颖pso-sade的滚动轴承可靠性多目标优化设计方法

Country Status (1)

Country Link
CN (1) CN109241695B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011069681A (ja) * 2009-09-25 2011-04-07 Ntn Corp 転がり軸受の余寿命推定方法および推定装置
CN105426630A (zh) * 2015-12-13 2016-03-23 北京工业大学 一种基于粒子群算法的角接触球轴承优化设计方法
CN106682688A (zh) * 2016-12-16 2017-05-17 华南理工大学 基于粒子群优化的堆叠降噪自编码网络轴承故障诊断方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011069681A (ja) * 2009-09-25 2011-04-07 Ntn Corp 転がり軸受の余寿命推定方法および推定装置
CN105426630A (zh) * 2015-12-13 2016-03-23 北京工业大学 一种基于粒子群算法的角接触球轴承优化设计方法
CN106682688A (zh) * 2016-12-16 2017-05-17 华南理工大学 基于粒子群优化的堆叠降噪自编码网络轴承故障诊断方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于改进SADE算法的神经网络预测储层物性;李虎等;《测井技术》;20121231;第36卷(第06期);第585-589页 *

Also Published As

Publication number Publication date
CN109241695A (zh) 2019-01-18

Similar Documents

Publication Publication Date Title
CN106649982B (zh) 风力发电机大锥角圆锥主轴承摩擦力矩计算方法
CN110765617B (zh) 基于润滑理论的圆柱滚子轴承滚子对数修形的设计方法
CN103256299B (zh) 一种变桨轴承结构参数的确定方法
Duggirala et al. Design optimization of deep groove ball bearings using crowding distance particle swarm optimization
Parmar et al. Macro geometry multi-objective optimization of planetary gearbox considering scuffing constraint
Dandagwhal et al. Design optimization of rolling element bearings using advanced optimization technique
CN109241695B (zh) 基于新颖pso-sade的滚动轴承可靠性多目标优化设计方法
Ahmad et al. Multi-objective robust optimization of deep groove ball bearings considering manufacturing tolerances based on fatigue and wear considerations
CN110259830A (zh) 一种滚道表面带有微织构的球轴承及其润滑方法
Kalyan et al. Multi-objective optimization in geometric design of tapered roller bearings based on fatigue, wear and thermal considerations through genetic algorithms
CN111222207A (zh) 一种圆锥滚子轴承内圈大挡边最优锥角的设计方法
CN205876965U (zh) 一种耐磨损球轴承
CN109578431B (zh) 一种流体动压组合轴承及轴承中表面微造型的设计方法
Mannan et al. Temperature effect on load distribution, friction, and wear of a grease-lubricated spherical roller bearing (SRB)
Fang et al. Analytical determination of the optimal clearance for the fatigue life of ball bearing under different load conditions
CN105160090A (zh) 一种新型陶瓷球轴承优化设计方法
Grujicić et al. The analysis of impact of intensity of contact load and angular shaft speed on the heat generation within radial ball bearing
Dragoni Optimal design of paired tapered roller bearings under centred radial and axial static loads
CN112347581B (zh) Rv减速机陶瓷球主轴承的选配方法及选配系统和存储介质
CN214499763U (zh) 一种具有抗震功能的大型推力球轴承
Li et al. The design method for surface texture of sliding friction pairs based on machine learning under mixed lubrication
CN114638069B (zh) 一种四点接触球轴承设计方法
Dabbara et al. Multi-Objective Optimization of Deep Groove Ball Bearings Using Fatigue-Wear-Thermal Considerations Through Genetic Algorithms
Horng et al. Investigation of running-in process based on surface roughness parameters, real contact area ratio and tribological properties
CN113536476B (zh) 一种流体动压润滑径向滑动轴承优化制备方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant