删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

采煤机摇臂系统行星轮系疲劳可靠性灵敏度设计

本站小编 Free考研考试/2020-03-23

张义民, 王婷, 黄婧
东北大学 机械工程与自动化学院, 辽宁 沈阳 110819
收稿日期: 2015-09-01
基金项目: 国家重点基础研究发展计划项目(2014CB046303);国家自然科学基金资助项目(51135003).
作者简介: 张义民(1958-), 男, 吉林长春人, 东北大学教授, 博士生导师。

摘要: 由于采煤机摇臂系统的行星轮系结构比较复杂,采用集中质量参数法对其进行有限元建模,分析行星轮与太阳轮的动态接触应力以确定其失效模式.主要考虑行星轮与太阳轮结构尺寸的随机性,利用BP神经网络的非线性映射功能,模拟得到疲劳寿命与随机参数的关系表达式.采用最大可能点摄动法进行可靠性设计.最后通过对行星轮与太阳轮的可靠性灵敏度设计,得到了各参数均值和方差对结构可靠性的影响情况.Monte-Carlo仿真试验验证了所提方法的正确性,为行星轮与太阳轮可靠性灵敏度设计提供了理论参考.
关键词:行星轮与太阳轮有限元建模BP神经网络最大可能点摄动法可靠性灵敏度设计
Fatigue Reliability-Based Sensitivity Design of Planet Gear for Shearer Rocker Arm System
ZHANG Yi-min, WANG Ting, HUANG Jing
School of Mechanical Engineering & Automation, Northeastern University, Shenyang 110819, China
Corresponding author: ZHANG Yi-min, E-mail: neuwangting@126.com
Abstract: As the planet gears structure of shearer rocker arm system is very complicated, a finite element model was built using the lumped-mass parameter method, and the dynamic contact stress of the planet gear and sun gear was analyzed to determine its failure mode. Considering the randomness of structural dimension for the planet gear and sun gear, a relational expression between the fatigue life and random parameters was obtained by using the nonlinear mapping function of BP neural network. The MPPPM was applied to the reliability design. The effect of the parameters' mean and variance on the structure reliability was achieved by the planet gear and sun gear reliability-based sensitivity design.The simulation experiment using the Monte-Carlo showed the correctness of the as-proposed method, and provided a theoretical reference for the reliability-based sensitivity design of the planet gear and sun gear.
Key Words: planet gear and sun gearfinite element modelBP neural networkMPPPM (most probable point perturbation method)reliability-based sensitivity design
摇臂系统是采煤机重要部件之一,其主要作用就是把动力传递给截割部的滚筒.现代采煤机的摇臂系统大部分都由直齿轮系和行星轮系构成.其中,行星轮系主要承担减速和传递大扭矩的作用,所承受的载荷复杂,且处于恶劣的地下工况,一旦发生故障将无法继续工作,而齿轮的接触疲劳是影响其使用寿命的主要因素之一[1],齿面处易发生点蚀、磨损和胶合等疲劳失效模式,从而造成齿轮的严重失效.因此对其进行结构可靠性设计是十分必要的.文献[2-4]以行星轮系的体积和最小为目标函数,进行可靠性优化设计,但未以齿轮的接触疲劳为目标进行结构可靠性设计;文献[5]采用罚函数法对行星轮系进行可靠性优化设计,但没有进行可靠性灵敏度设计;文献[6]分析了行星轮系的系统参数对其固有频率和振动模态的灵敏度,但没有对其进行疲劳可靠性灵敏度设计.
考虑到齿轮传动系统的结构特点及工程实际下的疲劳寿命,很难确定极限状态函数,而数值模拟工作量太大,为此采用传统的响应面法.与多项式拟合相比,BP神经网络具有更强的非线性逼近能力,并逐渐应用于各种机构的分析当中[7].本文采用BP神经网络模拟随机设计变量与失效寿命之间的表达式,然后应用随机摄动法对已训练的神经网络进行计算,完成对其结构可靠性灵敏度设计.与蒙特卡洛方法的结果进行对比,验证了此方法的精确性和高效性.
1 行星轮与太阳轮的有限元分析1.1 行星轮系建模研究对象为某型采煤机摇臂系统的行星轮系,行星轮系的基本参数如表 1所示.由于太阳轮所在轴为齿轮轴,且内齿圈与外壳固定在一起,故其轴承刚度ks=kr=0.齿轮和齿轮轴的材料参数:密度ρ为7 850 kg/m3,弹性模量E为2.1×1011Pa,泊松比μ为0.3.
表 1(Table 1)
表 1 行星轮系的基本参数Table 1 Basic parameters of planet gears
基本参数 太阳轮s 内齿圈r 行星架c 行星轮p
质量/kg 10.770 83.189 702.813 23.495
8.754 96.212 612.136 19.063
135.315 575.01 387 211.430
3.418×109 4.665×109 - -
- - 3.529×109 1.215 79×109
- 2×1010 - -
啮合角/(°) 26.364 67 19.999 98 - -


表 1 行星轮系的基本参数 Table 1 Basic parameters of planet gears

本文参照文献[8],采用mass21单元将太阳轮、齿圈、行星架及4个行星轮分别简化成集中质量单元;采用matrix27单元将模型中所有轴承及由行星轮分别与太阳轮、齿圈构成的齿轮副简化成弹簧单元.图 1为行星轮系的集中参数有限元模型,由固定坐标系x-y和旋转坐标系ζ-η组成.
图 1(Fig. 1)
图 1 行星轮系的有限元模型Fig.1 Finite element model of planetary gears

1.2 有限元分析结果通过行星轮系的集中参数系统有限元模型及实验测得的电机动态轴输入载荷,对系统进行动力学特性分析, 如图 2所示.分析得到行星轮与太阳轮的动态齿面接触应力较大,4个行星轮中最大的动态齿面接触应力如图 3所示.
图 2(Fig. 2)
图 2 电机轴动态输入载荷Fig.2 Dynamic input load of motor shaft

图 3(Fig. 3)
图 3 行星轮与太阳轮的动态接触应力Fig.3 Dynamic contact stress of planet gear and sun gear

图 3可知,行星轮与太阳轮啮合时的齿面接触应力最大峰值为1 782 MPa,大于材料的接触疲劳极限1 519 MPa.基于上述分析可知,在截割简单煤层工况下,采煤机截割部传动系统中,齿轮的失效形式主要是由于接触应力过大产生的齿面点蚀,应对其进行结构可靠性设计优化.
2 接触疲劳寿命拟合2.1 随机变量的抽样齿轮在实际设计、加工制造与装配过程中产生的误差,造成其几何尺寸具有随机性,并服从正态分布.均值根据设计要求确定,方差根据准则确定,具体取值如表 2所示.对这5个随机变量采取拉丁超立方抽样技术,设置样本容量为300,即可得到所需的样本.
表 2(Table 2)
表 2 随机变量Table 2 Random variables
随机变量 均值 标准差 备注
B2/mm 130 0.216 7 行星轮系的工作齿宽
a2/mm 193.5 0.007 7 行星轮与太阳轮的中心距
a3/mm 193.5 0.012 0 行星轮与内齿圈的中心距
α2/(°) 20 0.033 3 行星轮系啮合的压力角
x2/mm 0.659 2 6.776 7e-4 行星轮的变位系数


表 2 随机变量 Table 2 Random variables

2.2 计算接触疲劳寿命为了精确估算疲劳寿命,采用雨流计数法对已求得的齿轮接触线各单元的应力载荷谱进行统计分析,得到各单元的各级应力均值、应力幅值及相应的应力循环数.由于应力均值对疲劳累积损伤的影响,本文利用Goodman直线[9]对雨流计数的结果修正为对称循环载荷谱.本文中行星轮与太阳轮的材料为18Cr2Ni4WA,抗拉强度Su=1 180 MPa.
Goodman直线表达式为
(1)
式中:Sij为单元的对称循环应力;Saij为单元j的第i级应力幅值;Smij为单元j的第i级应力均值.
查阅文献[10],得到材料为18Cr2Ni4WA的齿轮接触疲劳S-N曲线,其方程为
(2)
式中:N为应力循环次数;S为对应N的条件疲劳极限;m为由试验确定的指数;C为由试验确定的常数.曲线方程参数值如下:m为11.001 1,lgC为42.700 8,接触疲劳应力σb为1 519.4 MPa,疲劳曲线拐点对应循环次数为5×107.
当疲劳载荷具有多个载荷水平时,选用Miner线性疲劳累积损伤理论,对各级应力循环造成的损伤进行累加,由此完成对行星轮与太阳轮接触疲劳寿命的估算:
(3)
式中:k为载荷水平的个数;ni为第i个载荷水平的循环数;Ni为第i个载荷水平下的疲劳寿命;Nf为安全运行的次数.
2.3 BP神经网络拟合疲劳寿命本文采用BP神经网络, 以随机变量B2a2a3α2x2为网络的输入参数,疲劳寿命Nf为输出参数,建立5-12-1三层神经网络.其中,隐含层的传递函数选用Logsig函数,输出层的传递函数选用Purelin函数,训练函数选用Trainlm函数,则疲劳寿命与随机变量[B2a2a3α2x2]=[X1X2X3X4X5]之间的函数关系可以近似为
(4)
式中:b1为输入层与隐含层间的阈值;b2为隐含层与输出层间的阈值;w1为输入层与隐含层间的权重;w2为隐含层与输出层间的权重.通过计算,求出该式中4个参数的取值如下所示.
为提高训练、收敛速度及BP神经网络的稳定性,每次仿真都对训练样本的输入和输出参数进行归一化处理:
(5)
(6)
经过BP神经网络训练后,得到行星轮与太阳轮接触寿命的训练样本为300个,训练后的均方误差为1.852 92×10-2(目标误差为2×10-2),神经网络计算结果和有限元分析结果的最大误差为7.937 25×10-2,可知采用BP神经网络模拟的功能函数较为精确.
3 结构可靠性灵敏度设计3.1 行星轮与太阳轮的结构可靠性分析利用神经网络拟合出的疲劳寿命Nf与随机矢量X=[X1X2X3X4X5]的功能函数为Nf(X),则行星轮与太阳轮的极限状态函数为
(7)
式中:N0为给定的疲劳寿命.当极限状态函数g(X)≤0时,行星轮与太阳轮为失效状态;当g(X) > 0时,行星轮与太阳轮为安全状态.g(X)=0表示一个五维极限状态曲面,通常称为失败面.
Xg(X)表示为
(8)
(9)
式中:0 < |ε| < < 1为一小参数;下标d为确定成分; 下标p为随机成分,且均值为零.要求确定部分要远大于随机部分.
g(X)取一阶矩、二阶矩、三阶矩和四阶矩,根据Kronecker[11]代数、相应的随机分析理论及矩阵值函数的Taylor展开式,有
(10a)
(10b)
(10c)
(10d)
式中,Var (X),C3(X),C4(X)为随机参数的二阶矩、三阶矩和四阶矩矢量.
可靠性指标定义为
(11)
求得可靠性指标后,通过四阶矩技术和Edgeworth级数将未知的状态分布函数展开成标准的正态分布函数,从而确定行星轮与太阳轮的结构可靠度.Edgeworth级数可表示为
(12)
式中:Hj(y)为j阶Hermite多项式,递推关系为
(13)
此系统的可靠度为
(14)
R > 1时,采用下述经验公式进行修正:
(15)
由于行星轮与太阳轮接触寿命的分散性较大,若N0=100 d,可求得行星轮与太阳轮系的可靠度指标及可靠度: β=4.484 866 896 233 15,R=0.997 103 650 217 132.
为了验证摄动法的正确性,本文使用MATLAB软件进行了Monte-Carlo仿真试验,随机抽样次数为1×106,得到相应的可靠度RMCS=1,可以看出采用摄动法求得的行星轮与太阳轮的可靠性是比较准确的.
3.2 行星轮与太阳轮的结构可靠性灵敏度设计行星轮与太阳轮的可靠度对基本随机参数矢量X的均值和方差的灵敏度分别为
(16)
(17)
式中,
(18)
(19)
(20)
(21)
(22)
(23)
(24)
R > 1时,可靠性指标β的灵敏度为
(25)
式中,
将式(25)替换式(16)和式(17)中的?R(β)/,即可求得均值和方差的灵敏度.
把已知条件和计算结果代入式(16)和式(17),得到各随机参数的均值和方差灵敏度:
此结果与采用一次二阶矩法分析的结果在数值上虽有差异,但总体趋势是一致的.对于均值灵敏度,行星轮的压力角α2和行星轮的齿宽B2的灵敏度值较高,且灵敏度符号为正值,可知在一定范围内,行星轮系的可靠度随α2B2的增加而增大.而其余尺寸随机参数对其影响可忽略, 因此,在满足行星轮系体积、质量和成本要求的前提下,适当增大行星轮的压力角与齿宽,行星轮系的可靠性会有所提高.
4 结论本文对采煤机摇臂传动系统的行星轮系进行了有限元建模和行星轮与太阳轮的动态接触疲劳应力分析,由分析结果知,齿面动态接触应力峰值确定了行星轮与太阳轮的失效模式.采用BP神经网络方法模拟出的疲劳寿命表达式适用于复杂结构、精度要求较高的可靠性灵敏度设计.Monte-Carlo仿真试验验证了所提方法的正确性.本文方法适用于滚筒式电牵引采煤机,截割功率≤2×300 kW,牵引功率≤2×40 kW.结果表明,采用本文方法得到了采煤机摇臂系统各个结构尺寸参数的均值和方差对行星轮系可靠性的影响情况,具备了较高的应用价值.
参考文献
[1]Wang J, Wang Y, Huo Z. Finite element residual stress analysis of planetary gear tooth[J].Advances in Mechanical Engineering, 2013(5) : 1–12.
[2]Wang S J, Lyu J. Reliability optimum design of NWG-type oil-submerged planetary gear reducer[J].Key Engineering Materials, 2005, 297/298/299/300 : 1895–1900.
[3]Li X, Jiang S, Li S, et al. Nonlinear transient engagement characteristics of planetary gear train[J].Journal of Vibroengineering, 2013, 15(2) : 933–941.
[4]Yang J, Zhang C, Lin Z, et al. Dynamic load analysis of bearings in a three-ring planetary gear transmission[J].Journal of Mechanical Engineering, 2004, 17(2) : 827–829.
[5] Ericson T, Parker R G.Planetary gear modal properties and dynamic response:experiments and analytical simulation[C]//Proceedings of the ASME 2011 International.Washington D C:American Society of Mechanical Engineers, 2011:331-343.
[6]Lin J, Parker R G. Sensitivity of planetary gear natural frequencies and vibration modes to model parameters[J].Journal of Sound and Vibration, 1999, 228(1) : 109–128.DOI:10.1006/jsvi.1999.2398
[7]朱丽莎, 张义民, 卢昊, 等. 基于神经网络的转子振动可靠性灵敏度分析[J].计算机集成制造系统, 2012, 18(1) : 149–155.
( Zhu Li-sha, Zhang Yi-min, Lu Hao, et al. Reliability sensitivity analysis of rotor vibration based on neural network[J].Computer Integrated Manufacturing Systems, 2012, 18(1) : 149–155.)
[8]Lin J, Parker R G. Analytical characterization of the unique properties of planetary gear free vibration[J].Journal of Vibration and Acoustics, 1999, 121(3) : 316–321.DOI:10.1115/1.2893982
[9]陈传尧. 疲劳与断裂[M].武汉: 华中科技大学出版社, 2002: 3-33.
( Chen Chuan-yao. Fatigue and fracture[M].Wuhan: Huazhong University of Science and Technology, 2002: 3-33.)
[10]Wu Z F, Wang T, Zhang R. Experimental study of contact fatigue property for 18Cr2Ni4WA gears[J].Acta Armamentarii, 2010, 31(5) : 598–602.
[11]Vetter W J. Matrix calculus operations and Taylor expansions[J].SIAM Review, 1973, 15(2) : 352–356.DOI:10.1137/1015034

相关话题/可靠性 系统

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于多变量阶跃响应的风洞系统辨识
    张廷丰,毛志忠,袁平,张弼东北大学信息科学与工程学院,辽宁沈阳110819收稿日期:2015-07-08基金项目:国家自然科学基金资助项目(61333006)。作者简介:张廷丰(1978-),男,辽宁朝阳人,东北大学博士研究生,辽宁工业大学讲师。通信作者:毛志忠(1961-),男,山东莱州人,东北大 ...
    本站小编 Free考研考试 2020-03-23
  • 具有离散和分布时滞系统的鲁棒有限时间能量-峰值控制
    郑连伟东北大学理学院,辽宁沈阳110819收稿日期:2015-12-28基金项目:辽宁省自然科学基金资助项目(2014020022)。通信作者:郑连伟(1963-),男,辽宁新民人,东北大学副教授,博士.E-mail:zhenglianwei@sina.com。摘要:研究了不确定线性时滞系统在有限时 ...
    本站小编 Free考研考试 2020-03-23
  • 自然循环汽化冷却系统水-水引射的行为
    冯明杰,王恩刚,冯慧聪东北大学材料电磁过程研究教育部重点实验室,辽宁沈阳110819收稿日期:2015-04-29基金项目:国家自然科学基金资助项目(51171041)。作者简介:王恩刚(1962-),男,辽宁沈阳人,东北大学教授,博士生导师。通信作者:冯明杰(1971-),男,河南禹州人,东北大学 ...
    本站小编 Free考研考试 2020-03-23
  • 基于Aspen Plus的化学链制氧系统模拟分析
    于庆波,任慧来,彭家燕,王坤东北大学冶金学院,辽宁沈阳110819收稿日期:2015-04-29基金项目:国家自然科学基金资助项目(51576035);国家科技支撑计划项目(2013BAA03B03)。通信作者:于庆波(1966-),男,山东莱阳人,东北大学教授,博士生导师.。摘要:在热量分析和压力 ...
    本站小编 Free考研考试 2020-03-23
  • 基于有限元法的人体座椅系统随机响应分析
    董瑞春,郭立新东北大学机械工程与自动化学院,辽宁沈阳110819收稿日期:2015-07-13基金项目:国家自然科学基金资助项目(51275082,51135003);中央高校基本科研业务费专项资金资助项目(NI30403009)。作者简介:董瑞春(1985-),男,山东日照人,东北大学博士研究生。 ...
    本站小编 Free考研考试 2020-03-23
  • 错过法学院招收推免生的预报名系统,贵院法硕(非法学)还会从全国推免系统中选择推免生参加复试吗?
    提问问题:错过法学院招收推免生的预报名系统,贵院法硕(非法学)还会从全国推免系统中选择推免生参加复试吗?学院:凯原法学院提问人:18***38时间:2018-09-2114:14提问内容:老师您好!我一直以来十分心仪上交,一直等着贵校法学院法硕的推免报名。但由于没能及时得到相关信息,以至于错过了法学 ...
    本站小编 上海交通大学 2019-11-25
  • 系统生物医学研究院联系电话
    提问问题:系统生物医学研究院联系电话学院:生命科学技术学院(含系统生物医学研究院)提问人:18***86时间:2018-09-2109:27提问内容:请问下系统生物医学研究院联系电话?回复内容:见简章http://yzb.sjtu.edu.cn/info/1022/2263.htm ...
    本站小编 上海交通大学 2019-11-25
  • 上交819信号系统与信号处理在考试时能带计算器吗?
    提问问题:上交819信号系统与信号处理在考试时能带计算器吗?学院:电子信息与电气工程学院提问人:13***96时间:2018-09-1913:40提问内容:上交819信号系统与信号处理在考试时能带计算器吗?去年的学长说好像可以带,那么如果能带的话有什么限制吗?回复内容:见准考证上的说明 ...
    本站小编 上海交通大学 2019-11-25
  • 交大推免系统不用预报名吗?
    提问问题:交大推免系统不用预报名吗?学院:安泰经济与管理学院提问人:13***85时间:2018-09-1909:30提问内容:老师好,请问推免生只在研招网上报名就可以了,不用在交大招生系统上预报名了是吗?回复内容:学校没有预报名,具体咨询报考学院。http://yzb.sjtu.edu.cn/zs ...
    本站小编 上海交通大学 2019-11-25
  • 计算机系统考优惠
    提问问题:计算机系统考优惠学院:电子信息与电气工程学院提问人:15***82时间:2015-09-2314:03提问内容:请问统考优惠,进复试的结果什么时候出来?回复内容:到时候直接咨询学院。 ...
    本站小编 上海交通大学 2019-11-25