基于ISIGHT的双列圆锥滚子轴承优化设计.docx

1、基于ISIGHT的双列圆锥滚子轴承优化设计基于ISIGHT的双列圆锥滚子轴承优化设计 基于ISIGHT的双列圆锥滚子轴承优化设计 李 震，郑林征，闫振勇，李 铎 （大连理工大学 机械工程学院，辽宁 大连 116024） 摘 要：在考虑轴承受到径向载荷、轴向载荷和弯矩载荷的复杂受载情况下，建立了双列圆锥滚子轴承的优化设计模型，以轴承整体寿命为优化目标，以每列轴承滚子个数、圆锥滚子平均直径、圆锥滚子有效接触长度和轴承的节圆直径为设计变量，基于ISIGHT设计优化软件对设计模型进行了试验设计（DOE）分析，利用多岛遗传算法对优化设计模型进行了求解。将该模型应用于3MW风电主轴承，得到了试验设计分析下

2、四个设计变量对轴承整体寿命影响的Pareto图和主效应图，参数优化后的风电主轴承整体寿命提高了83.8%。 关键词：双列圆锥滚子轴承；优化设计；ISIGHT；多岛遗传算法；试验设计 1 引言 圆锥滚子轴承能够承受较大的径向力、轴向力和弯矩，广泛应用于机床主轴、汽车前后轮、火车轮轴、建筑机械、各种减速机和航空航天领域，双列圆锥滚子轴承由于其能够承受大弯矩的结构特点在风电设备领域得到了广泛的应用。 目前轴承的优化设计主要集中于深沟球轴承和圆柱滚子轴承，关于圆锥滚子轴承的研究还比较少。文献1首先研究了双列圆锥滚子轴承的载荷分布情况，文献2提出了高转速下圆锥滚子轴承的分析方法，文献3基于向量运算综合研

3、究了双列圆锥滚子轴承的受力及润滑状况，文献4利用遗传算法对圆锥滚子轴承进行了优化设计，文献5-6基于改进的遗传算法和自适应算法对圆锥滚子轴承进行了优化设计。 ISIGHT作为一款功能强大的多学科优化设计软件，能够多快好省地改进产品的性能，有效降低产品的设计成本，在机械设计领域得到了广泛的应用7-9。兆瓦级风电主轴承是典型的双列圆锥滚子轴承，由于其庞大的体积和复杂的加工制造性，其设计、制造成本一直居高不下。将ISIGHT优化软件应用于3MW风电主轴承，通过ISIGHT中的试验设计（DOE）模块和优化设计模块对其主要参数进行了优化设计，优化后的主轴承的寿命提高了83.8%，为企业的设计制造提供了数

4、据支持。 2 双列圆锥滚子轴承优化模型 2.1 双列圆锥滚子轴承几何模型 双列圆锥滚子轴承的结构图及主要外形参数，如图1所示。图中：d轴承内圈公称直径；D轴承外圈公称直径；Dr轴承节圆直径；T轴承总宽度；C轴承外圈宽度；Dm两列滚子在轴向方向的中心间距；公称接触角。 图1 双列圆锥滚子轴承结构图Fig.1 Schematic View of a Double-Row Tapered Roller Bearing 一般情况下，圆锥滚子轴承的破损主要是滚子与内、外圈滚道的疲劳破损，文献4-6中均选用轴承的疲劳寿命作为优化设计目标。对于具有线接触特征的圆锥滚子轴承10，其疲劳寿命计算公式为： 式中：

5、L轴承的基本额定寿命；Qc滚道的额定动载荷；Qe与滚子实际载荷分布相关的当量载荷。 在轴承承载相同时，Qc越大，其额定寿命L越长，因此通常将Qc作为轴承优化的目标函数。但上述优化方法仅考虑了轴承仅受轴向载荷和径向载荷时的受载情况，并未考虑弯矩载荷的影响，而在风电主轴承的实际受载中，弯矩载荷占了很大的比重，其对轴承内部载荷分布有很大的影响，对于额定寿命的影响也是不能忽略的。因此针对其受载特点，对双列圆锥滚子轴承的基本额定寿命提出了一种新的优化模型。 对于双列圆锥滚子轴承，滚道额定动载荷为： 式中：Lm1、Lm2每列轴承的修正额定寿命。 2.2 目标函数 上述模型的计算过程考虑了轴承内、外圈相对于

6、外载荷静止或者转动的影响，并且在计算的过程中考虑了每个圆锥滚子的受力情况，因此可以将轴承整体受力中的径向力、轴向力和弯矩综合考虑，相比于单纯的通过额定动载荷的判断方法更合理，因此本次优化设计中选取双列圆锥滚子轴承修正后的疲劳寿命Lm作为优化的目标函数。因此，优化的目标函数可以表示为： 通过前面的讨论得，目标函数Lm是Dm，Le，Z，Dr和的函数。式中：Dm圆锥滚子平均直径，mm；Le圆锥滚子与内外圈滚道 的有效接触长度，mm；Z每列轴承的滚子数；Dr圆锥滚 子轴承的节圆直径，mm；公称接触角。 2.3 设计变量 设计变量是模型优化过程中需要优化的参数。由上述计算过程可知，目标函数Lm是Dm，L

7、e，Z，Dr和函数，因此优化可选取五个设计变量，即：圆锥滚子的平均直径Dm，圆锥滚子与内外圈滚道的有效接触长度Le，每列轴承的滚子数Z，圆锥滚子轴承的节圆直径Dr，公称接触角。根据企业提供的技术规范，公称接触角现在已成规范系列，因此公称接触角在优化设计中设置为定值，设计变量可以表示为： 2.4 约束条件 根据文献的研究4可知，圆锥滚子的平均直径Dm应该满足以下约束条件： 式中：KDmin、KDmax滚子轴径系数的极小值和极大值，其大小为： 0.3 KDmin0.4，0.5KDmax0.6。 轴承的节圆直径和平均直径应满足以下约束条件才能使轴承平稳转动，即： 式中：e约束常数，其取值范围为：0.

8、01e0.07。 圆锥滚子的有效接触长度Le在轴承的优化中是一个非常重要的参数，为了保证滚子大端不超出外圈端部，应满足： 在滚动轴承中，圆锥滚子的有效接触长度Le通常要大于等于圆锥滚子的平均直径，因此其满足如下约束条件： 轴承外圈底部的厚度应Dm，式中：常数，其取值范围为：0.40.5，因此其满足如下约束条件： 在圆锥滚子轴承中，滚子在内、外圈滚道间滚动，为了避免滚子间的碰撞，两个滚子之间必须留有一定的空隙，因此每列滚子的个数应该满足以下约束条件： 同时，为了满足滚子与内外圈间的几何关系，滚子的最小个数还应满足以下约束条件： 2.5 优化原理 ISIGHT参数优化软件由于其良好的集成性、兼容性

9、以及强大的数据后处理，在优化设计领域得到了广泛的应用。以ISIGHT软件为优化平台，集成双列圆锥滚子轴承疲劳寿命计算的MATLAB程序，通过ISIGHT和MATLAB联合完成对双列圆锥滚子轴承的优化。ISIGHT中的试验设计方法（DOE）能够帮助设计人员快速有效的得到优化参数的有效信息，通过DOE分析结果中的Pareto图，可以清晰直观的得到各优化参数对于目标函数的贡献百分比，通过分析结果中的主效应图，可以得到各优化参数对于优化目标函数的主效应，能够清晰地表达各参数对于目标函数的影响趋势。 ISIGHT与MATLAB联合优化过程，如图2所示。在联合优化过程中，首先利用ISIGHT的试验设计（D

10、OE）模块来分析各优化参数对于目标函数影响的重要程度，然后利用优化模块中多岛遗传算法来对模型进行优化求解，得到最终的优化结果。 图2 ISIGHT与MATLAB联合优化图Fig.2 Optimal Flow Diagram with ISIGHT and MATLAB 3 风电主轴承优化 某型号3MW风电主轴承，其结构参数，如表1所示。载荷工况为企业提供的均载载荷工况。 表1 3MW风电主轴承结构参数Tab.1 Structural Parameters of 3MW Bearing结构参数 值 单位滚子与外圈接触角 45.4 deg滚子与内圈接触角 42.6 deg滚子与挡边接触角 42.3

11、 deg滚子大端直径 79.7 mm两列滚子质心轴向间距 160.3 mm轴承节圆直径 2062 mm滚子有效接触长度 96.03 mm每列滚子个数 85 根据技术要求，此次优化中等效接触角的取值为定值=44。因此实际的优化参数为：Dm，Le，Z，Dr。将风电主轴承的零件图简化得到的外圈有效直径D=2156mm，内圈有效直径为d=2156mm，根据约束条件、约束常数的取值以及企业提供的技术参考，得到各优化变量的取值：圆锥滚子的平均直径：50Dm75；轴承的节圆直径：1790Dr2372；滚子有效接触长度：45Le100；每列轴承的滚子个数：43Z90 将四个参数的大小代入ISIGHT优化模型的

12、分析框架中，通过数据处理得到优化变量对寿命影响的Pareto图和各优化变量对轴承寿命影响的主效应，图如图3、图4所示。 图3 优化变量对寿命L影响的Pareto图Fig.3 Pareto Plot for Response L 图4 优化变量对轴承寿命影响的主效应图Fig.4 Main Effects Plot for L of the Design Parameters 由图3可知，在Dm，Le，Z，Dr四个优化变量中，圆锥滚子与内外圈滚道的有效接触长度Le是影响轴承整体寿命L的最关键因素，圆锥滚子的平均直径Dm次之，轴承的节圆直径Dr的影响最小，同时还可以得知：四个优化变量对于目标函数L均

13、为正相关。从图4中四个优化变量对于轴承整体寿命影响的主效应图中四条曲线的变化趋势也可以得到与图3同样的结论，因此在双列圆锥滚子轴承的设计过程中，可以通过增加圆锥滚子与内外圈滚道的有效接触长度或者增大圆锥滚子的平均直径来快速提高轴承的整体寿命。 在多岛遗传算法中，虽然种群个体的增大会提高遗传算法的搜索质量和防止种群在成熟前收敛，但是过多的种群个体会使整体的收敛速度降低。根据经验得知种群数的适当增大和岛数的适当减小能在保证优化结果的前提下减少迭代次数，从而节省优化时间。优化过程中轴承整体寿命随迭代次数的响应，如图5所示。 图 5目标函数优化过程Fig.5 Optimization Process

14、of the Objective Function 图5中的黑点表示迭代过程中的可行域，红点表示不可行域，最终定格的绿点表示最终优化的最优结果。将图5中得到的轴承整体寿命的优化数据以及四个优化变量的数据整理，所得结果如表2所示。从表2中可以看出，优化后每列轴承滚子数Z，圆锥滚子有效接触长度Le，轴承的节圆直径Dr以及圆锥滚子平均直径Dm均增大，轴承整体疲劳寿命由135725h变成了249462h，风电主轴承整体寿命相比优化前提高了83.8%，优化效果明显。 表2 优化前后轴承参数对比Tab.2 Comparison of Optimized Bearings Parameters参数 3MW风

15、电主轴承原参数 ISIGHT优化后的轴承参数每列轴承滚子个数：Z 85 88圆锥滚子接触长度：Le 96.03 107.47圆锥滚子平均直径：Dm 68 74.7轴承节圆直径：Dr 2062 2081轴承整体疲劳寿命：L 135725 249462 4 结论 建立了双列圆锥滚子轴承的优化设计模型，以轴承整体寿命的最大值为优化目标，同时考虑了轴向力、径向力和弯矩的综合影响；基于ISIGHT优化设计软件对双列圆锥滚子轴承优化模型的四个优化变量进行了试验设计（DOE）分析，同时应用多岛遗传算法完成了对其优化变量的优化设计；将该模型应用于某型号3MW风电主轴承，得到了四个优化变量Dm，Le，Z，Dr对

16、于轴承整体寿命L影响的Pareto图和主效应图，得出了圆锥滚子有效接触长度Le和圆锥滚子平均直径Dm是影响3MW风电主轴承整体寿命的关键因素，经过多岛遗传算法优化后风电主轴承的整体寿命相比于优化前提高了83.8%，说明了该优化设计模型的可行性以及正确性，为兆瓦级风电主轴承的设计制造提供了指导。 参考文献 1 AndrasonS.Load distribution in a taper roller bearing arrangement considering misalignmentJ.Tribology，1973，6（3）：84-92. 2Liu J Y.Analysis of taper

17、ed roller bearings considering high speed and combined loadingJ.Journal of Tribology，1976，98（4）：564-572. 3Bercea I，Cretu S，Nelias D.Analysis of double-row tapered roller bearings，Part I-ModelJ.Tribology Transactions，2003，46（2）：228-239. 4Tiwari R，Sunil K K，Reddy R S.An optimal design methodology of t

18、apered roller bearings using genetic algorithmsJ.International Journal for Computational Methods in Engineering Science and Mechanics，2012，13（2）：108-127. 5何绍武，邬义杰，周刚.基于改进遗传算法的圆锥滚子轴承优化设计方法的研究J.组合机床与自动化加工技术，2006（9）：1-3.（He Shao-wu，Wu Yi-jie，Zhou Gang.Research on the methods of optimal design of taper

19、roller bearing based on improved genetic algorithmJ.Molular Machine Toolamp;Automatic Manufacturing Technique，2006（9）：1-3.） 6赵章荣，隋晓梅，何绍武.基于改进自适应遗传算法的圆锥滚子轴承优化设计J.机床与液压，2008，36（8）：21-23.（Zhao Zhang-rong，Sui Xiao-mei，He Shao-wu.Optimal design of taper roller bearing based on improved adaptive genetic a

20、lgorithmJ.Machine Toolamp;Hydraulics，2008，36（8）：21-23.） 7王峰，吴永海.基于ISIGHT的陶瓷轴承优化设计J.机械设计与制造，2007（12）：96-97.（Wang Feng，Wu Yong-hai.Optimization design of the ceramic bearing based on ISIGHTJ.Machinery Designamp;Manufacture，2007（12）：96-97.） 8李刚，韩海兰，解瑞春.基于CarSim与ISIGHT的悬架K特性优化设计方法研究J.机械设计与制造，2015（1）：79-

21、82.（Li Gang，Han Hai-lan，Xie Rui-chun.Study on optimization design method of suspension k characteristics based on CarSim and ISIGHT softwareJ.Machinery Designamp;Manufacture，2015（1）：79-82.） 9王南，刘庆阳，周莎莎.基于ISIGHT平台NSGA-方法的3-PRS并联机构多目标优化J.机械设计与制造，2015（12）：198-201.（Wang Nan，Liu Qing-yang，Zhou Sha-sha.Mu

22、lti-objective optimization of 3-PRS parallel mechanism based on NSGA-II method of ISIGHT platformJ.Machinery Designamp;Manufacture，2015（12）：198-201.） 10罗继伟，罗天宇.滚动轴承分析计算与应用M.北京：机械工业出版社，2009：79-83.（Luo Ji-wei，Luo Tian-yu.Analysis Calculation and Application of Roller BearingsM.Beijing：China Machine Pr

23、ess，2009：79-83.） Optimal Design of Double-Row Tapered Roller Bearings Based on ISIGHT Software LI Zhen，ZHENG Lin-zheng，YAN Zhen-yong，LI Duo（School of Mechanical Engineering，Dalian University of Technology，Liaoning Dalian 116024，China） Abstract：An optimal design model of double-row tapered roller bea

24、rings was proposed with the consideration of radial load，axial load and thrust load.The fatigue life of the bearing was set as the objective function，the number of tapered rollers in one row，the average diameter of tapered roller，the effective contact length of the roller and the bearing pitch diame

25、ter were set as the design parameters.The design of experiments was conducted and the design model was solved based on the multi-island GA method.We get the Pareto plot and the main effects plot of the four design parameters for the objective function when the model is applied to 3MWwind power main

26、bearing.The fatigue life the bearing has improved 83.8 percent after the optimization. Key Words：Double-Row Tapered Roller Bearings；Optimization Design；ISIGHT；Multi-Island GA；Design of Experiments 中图分类号：TH16；TH122 文献标识码：A 文章编号：1001-3997（2017）11-0030-04 来稿日期：2017-05-20 基金项目：辽宁省科技创新重大专项项目大轴重铁路货车滚动轴承（2015106009）；辽宁省科技创新重大专项项目5MW及以上大功率风力发电机组配套轴承（201303006） 作者简介：李 震，（1975-），男，辽宁辽阳人，博士研究生，讲师，主要研究方向：轴承设计；郑林征，（1988-），男，河北石家庄人，硕士研究生，主要研究方向：轴承设计