基于有限元法的NW型行星架结构优化.pdf

132023,61(11)总第7 11期机械制造专题报基于有限元法的NW型行星架结构优化杜佐飞汤天宝余春祥于海生孟思宇无锡星驱动力科技有限公司江苏无锡214100摘要：介绍了NW型行星架结构及工作原理，基于有限元方法对NW型行星架结构进行了可靠性分析。根据分析结果，针对行星齿轮轴相对变形量与最大螺栓利用率两个不达标项进行结构优化。结构优化后,NW型行星架结构的行星齿轮轴相对变形量和最大螺栓利用率均满足要求。同时对NW型行星架结构进行疲劳损伤校核，确认NW型行星架结构发生疲劳损伤的风险较小，与高周扭转疲劳试验结果吻合。关键词：有限元行星架结构优化中图分类号：TH132.46文献标志码：A文章编号：10 0 0-49 9 8（2 0 2 3）11-0 0 13-0 6Abstract:The NW planetary carrier structure and its working principle were introduced,and thereliability analysis of the NW planetary carrier structure was carried out based on the finite element method.According to the analysis result,the structure was optimized for the two non-conforming items of planetary gearshaft relative deformation and maximum bolt utilization rate.After structural optimization,the planetary gearshaft relative deformation and the maximum bolt utilization rate of the NW planetary carrier structure meet therequirements.At the same time,fatigue damage check was carried out on the NW planetary carrier structure,and it is confirmed that the risk of fatigue damage of the NW planetary carrier structure is small,which isconsistent with the result of high-cycle torsional fatigue test.Keywords:Finite ElementPlanetaryCarrierStructureOptimization1优化背景当前，国内新能源汽车市场已经步人电动化发展时代。电驱系统是新能源整车三大系统中至关重要的机构，直接影响到汽车寿命及驾驶员的安全。电驱系统主要由电机总成、电机控制器和传动总成组成。其中，传动总成中的减速器在传动执行机构中起匹配转速和传递转矩的作用，目的是降低转速，增大转矩，以保证驱动电机的转矩、转速满足整车要求。行星齿轮减速器因行星排轴向空间要求而结构紧凑，具有回程间隙小、精度高、使用寿命长、额定输出扭矩大等优点，被广泛用于电驱系统的传动总成结构。行星架是行星齿轮减速器传动结构中受力较为复杂的零部件，作为输入构件，承受较大的外加扭矩。由此，行星架的结构强度和刚度对行星轮间的载荷分配及传动装置的承载能力、噪声、振动等有很大影响2 。目前，国内学者对减速器进行有限元分析的研究较多，但对NW型行星架进行有限元分析的研究则还比较少。吴迅2 利用有限元法，在差速减速器行星架前后两侧板和销轴不同配合状态下对行星架进行应力分析，结合分析结果对行星架前后两侧板与销轴之间的公差配合进行优化改进。赵丽娟等3 对矿用减速器行星架进行有限元分析，获得行星架组件的变形及应力分布。孙黎等4 对风电齿轮箱低速行星架进行极限强度分析，得到不同加载约束条件下的危险点分布情况。李辉云5 利用ADAMS构建减速器行星齿轮传动系统的虚拟样机，将齿轮接触间的摩擦因数引人仿真计算，得到齿轮间接触力随速度按非线性曲线变化的结论。涂杰等6 研究行星架销轴系统在额定工况下的应力与变形分布情况，为行星架销轴系统的设计和实际应用提供参考。笔者针对某电驱动系统中的NW型行星架结构，在研发初期为确保结构可靠性，采用有限元分析方法进行强度、疲劳分析，并进行结构优化。2NW型行星架结构及工作原理NW型行星架结构如图1所示，主要由电机轴、行星轮、行星齿轮轴、行星架、行星架后端板、外齿圈等组成。行星轮与电机轴齿轮、外齿圈同时接触，成为多接触结构7 O142023,61(11)总第7 11期机械制造专题报导双联齿大行星轮双联齿小行星轮电机输出端电机外齿圈图1NW型行星架结构行星轮为双联齿轮，分布为12 0，分别装配在三个行星齿轮轴上。行星轮同时与电机轴齿轮和外齿圈啮合。行星架是行星轮的支持构件，行星架后端板通过螺栓与行星架本体相连。外齿圈与电驱壳体固定在一起。以上零部件共同组成NW型行星架结构。NW型行星架结构通过电机驱动高速旋转的输入轴，直接带动行星轮转动。由于两者齿数不同，导致转速降低。外齿圈因固定在电驱壳体上不运动，便驱动行星架做输出运动。行星齿轮减速器的扭矩由电机轴直接输入，电机轴齿轮起太阳轮作用，与行星齿轮啮合，将扭矩传递至双联齿大行星轮上。双联齿小行星轮与外齿圈啮合，外齿圈固定在电驱总成壳体上，所以行星齿轮带着行星架绕外齿圈公转。行星架后端板与半轴相连，通过花键将扭矩传递至半轴，进而传递至整车轮毂，为整车提供动力。3有限元建模3.1模型有限元模型基于真实物理模型简化而来8 。本文中用于计算分析的有限元模型和真实物理模型一致，主要包括行星架、行星架后端板、行星齿轮轴、前后滚针轴承、行星轮止推垫片、前后球轴承、连接螺栓等，如图2 所示。滚针轴承对行星轮与行星齿轮轴起支撑作用，对滚针轴承进行简化处理，其余零件均按真实物理模型在三维软件中建模，并按实际工作状态进行装配，再进行有限元网格处理。行星架行星轮止推垫片行星架后端板前球轴承一行星齿轮轴后球轴承前滚针轴承后滚针轴承图2有限元模型整个有限元模型有限元网格数共为9 9 10 6 2，节点数共为17 0 9 13 4。计算采用修正二阶四面体C3D10M单元类型，采用Abaqus软件进行计算及后处理。3.2边界定义减速器行星机构在实际工作状态时，行星架前后球轴承均固定在壳体上，所以壳体对行星架前后球轴承外圈起支撑作用。同时球轴承因自身功能特性，内外圈间可相对滚动。对于计算边界，前后球轴承内外圈之间建立弹簧连接单元，赋予弹簧刚度，并对球轴承外圈进行15 个自由度约束。因行星架后端板花键轴为输出端，所以约束其旋转自由度。3.3载荷定义利用Master软件搭建电驱总成整机模型，如图3所示。根据整机实际工作情况，输出不同扭矩工况下NM型行星架结构有限元分析所需要的计算载荷。输出位置位于行星轮与行星齿轮轴中间的滚针轴承处，因为NW型行星架结构中每个行星齿轮轴上有前后行星轮，所以输出载荷时每个行星齿轮滚针轴承处有前后两组载荷。考虑偏载作用，需在行星齿轮轴止推垫片上加载相应轴向力及扭矩。电驱总成在工作时有正驱和倒拖等工况，正驱扭矩为415 Nm，倒拖扭矩为2 8 4Nm。计算载荷见表1，节点位置如图4所示。在NW型行星架结构设计初期，考虑成本及加工工艺等因素，行星架和行星架后端板分别采用不同材料加工。为了后期维护及拆装方便，行星架和行星架后端板不用焊接连接，而采用螺栓连接形式。螺栓规格为普通螺纹M8、12.9 级。根据计算所得螺栓预紧轴力为12.6 8 kN28.56kN。152023,61(11)总第7 11期机械制造专题报导A图3电驱总成整机模型因为正驱工况载荷大于倒拖工况，所以在后文，仅列出正驱415 Nm工况下的有限元分析结果。3.4材料属性定义对NW型行星架结构进行有限元分析，行星架和行星架后端板两个零部件的材料属性见表2。4有限元分析结果4.1评估内容对NW型行星架结构进行有限元求解计算后，考虑实际功能需求，主要从以下方面进行评估：对行星架及行星架后端板进行强度分析9 ；对行星齿轮表1计算载荷力/kN扭矩/（Nm）工况位置节点X向Y向Z向X向Y向Z向行星齿轮轴1-5.357-7.428-1.707 10 5-31.08924.79040行星齿轮轴行星齿轮轴28.57850.982 4-1.715 105-6.052836.9640前滚针轴承行星齿轮轴3-3.4948.0303-1.712 10-535.518113.47580正驱415 Nm行星齿轮轴12.641 9-11.247-1.708 10-582.124.120.7170行星齿轮轴行星齿轮轴28.4837.7821-1.716 10-5-58.41661.570 70后滚针轴承行星齿轮轴3-11.013.4094-1.713 10 523.65681.1930行星齿轮轴14.979.25.1390.000471922.456 4-24.2750行星齿轮轴行星齿轮轴2-6.9621.79450.000 475110.0731.7150前滚针轴承行星齿轮轴32.0022-6.90290.000472432.356-7.51060倒拖2 8 4Nm行星齿轮轴1-6.8037.40320.000472-56.048-55.0760行星齿轮轴行星齿轮轴2-3.0579.69850.000475176.9246-21.456 60后滚针轴承行星齿轮轴39.78372.20580.000 472 419.79775.801 40表2材料属性弹性模量/屈服强度/抗拉强度/零部件牌号泊松比GPaMPaMPa行星架QT7001740.305420700行星架20MnCrS52100.27810802.550后端板行星齿轮轴2（12 0）行星齿轮轴1（0）Z行星齿轮轴3（-12 0）图4节点位置轴进行相对变形量分析；对螺栓利用率进行分析；根据工程需求对NW型行星架结构进行疲劳校核。相对变形量Drelaive计算式为：Drelative/(x)-x,)+(y1-y2)2(1)式中：x1y1x2y2为行星齿轮轴端面PI、Pz 两点在最大扭矩下的位移结果坐标，如图5 所示。Pi(xiyl)Pi(x2.y2)图5相对变形量计算参数4.2应力云图在正驱415 Nm工况下，行星架最大应力为2 3 3MPa，行星架后端板最大应力为3 8 3 MPa，分别如图6、图7 所示，均小于材料屈服强度。4.3行星齿轮轴相对变形在正驱415 Nm工况下，行星齿轮轴相对变形最大值为0.16 3 mm，变形量较大。行星齿轮轴相对变形见表3，目标为不大于0.0 5 mm。162023,61(11)总第7 11期机械制造专题报导+2.737e+03+3.500e+02+2.900e+02+2.300e+02+1.700e+02+1.100e+02+5.000e+01+1.125e-01(MPa)N:19829应力:2 3 3 MPa4图6行星架应力云图+3.560e+03+6.000e+02+4.900e+02+3.800e+02+2.700e+02+1.600e+02+5.000e+01+4.314e-04(MPa)N:1071221应力：3 5 0 MPaN:1092541应力：3 8 3 MPa图7行星架后端板应力云图表3 行星齿轮轴相对变形节点行星齿轮轴1行星齿轮轴2行星齿轮轴3相对变形量/mm0.1630.1610.1624.4螺栓利用率在正驱415 Nm工况下，行星架拧紧螺栓在实际工作状态下，最大螺栓利用率为9 9.9 7%，螺栓有发生松脱风险。螺栓利用率见表4，目标为低于9 5%。表4螺栓利用率螺栓标号实际切向力/N允许切向力/N最大螺栓利用率115091510.499.93%215191520.299.96%315091510.599.93%415201520.899.96%515091510.399.93%615201520.899.97%通过上述分析，可以看出NW型行星架结构在前期设计状态下，行星架及行星架后端板强度没有问题，但行星齿轮轴变形量较大，这将对行星齿轮在传递扭矩时产生严重的错位误差，进一步导致行星齿轮间产生敲击等问题。螺栓利用率过高，将导致行星架在工作期间螺栓扭紧力矩不够，行星架后端板在传递扭矩时与行星架发生松脱，进而使行星机构无法正确传递扭矩，使整个电驱总成无法正常工作，甚至会产生严重的安全事故。必须针对以上两点进行优化。5优化方案分析5.1优化方案针对上述分析，根据实际工程项目经验，对NW型行星架结构提出如下优化方案：将行星齿轮轴与行星架间由间隙配合调整为小过盈配合；将行星架与行星架窗口区域圆角连成一体，提高此区域的连接刚度，如图8 所示；增大行星齿轮轴的直径，提高弯曲刚度；在螺栓规格保持不变的情况下，螺栓拧紧方式采用扭矩+转角形式，增大螺栓拧紧力，或直接更改螺栓直径及扭矩，增大预紧力。A图8行星架窗口区域改进针对以上优化方案，因为NW型行星架结构设计考量及行星齿轮轴是统一沿用件，所以第条未采纳，其余均采纳。在螺栓规格保持不变的情况下，最小预紧轴力增大至2 6.0 3 kN。行星架与行星齿轮轴由间隙配合调整为-0.0 2 6 mm小过盈配合。对于行星架板与行星架窗口区域，适当提高连接刚度。5.2优化后有限元分析优化后，在正驱415 Nm工况下，行星架最大应力为2 44MPa，行星架后端板最大应力为3 7 1MPa，分别如图9、图10 所示，行星架强度满足要求。+3.530e+03+3.500e+02+2.900e+02+2.300e+02+1.700e+02+1.100e+02+5.000e+011+2.940e-01(MPa)N:2809376应力:2 44MPa八图9优化后行星架应力云图优化后，行星齿轮轴相对变形最大值为0.0 8 1mm，见表5，较优化前改善明显。优化后，最大螺栓利用率为9 3.5%，见表6，尽管安全余量较小，但是满足目标要求。优化后，螺栓扭紧力矩满足要求。172023,61(11)总第7 11期机械制造专题报导+4.141e+03+6.000e+02+4.900e+02+3.800e+02+2.700e+02+1.600e+02+5.000e+01+4.740e-04(MPa)N:2125171应力：3 40 MPaN:2161132应力：3 7 1MPa图10优化后行星架后端板应力云图表5 优化后行星齿轮轴相对变形节点行星齿轮轴1行星齿轮轴2行星齿轮轴3相对变形量/mm0.0800.0810.079表6优化后螺栓利用率螺栓标号实际切向力/N允许切向力/N最大螺栓利用率12.887.13089.493.5%22.658.53108.685.5%32.867.83090.392.8%42.627.53110.184.5%52.8753089.793.1%62.6373.109.284.8%针对行星齿轮轴相对变形量与目标值差距较大这一情况，经与项目专家评审，一致决定可对标同类型竞品机的类似结构。对同类型竞品机进行分析，保持所有外界边界条件一致。经分析，竞品机行星齿轮轴相对变形最大值为0.12 0 mm，见表7。竞品机行星齿轮轴相对变形大于优化后的行星齿轮轴相对变形。表7 竞品机行星齿轮轴相对变形节点行星齿轮轴1行星齿轮轴2行星齿轮轴3相对变形量/mm0.1170.1200.118后期经过对比及评审，一致认为优化后的行星齿轮轴相对变形可以接受。6疲劳分析6.1累积损伤寿命计算根据设计验证需要，对NW型行星架结构进行高周疲劳分析，其中，高周疲劳载荷工况由试验大纲载荷谱而来，见表8。进行高周疲劳分析，先根据各循环工况进行强度计算，再根据每个工况的应力结果，在疲劳软件中叠加循环次数，进行损伤累积计算10-经过高周疲劳计算，NW型行星架结构累积损伤寿命最大为2.9 8 10-，行星架后端板累积损伤寿命表8高周疲劳试验大纲载荷谱序号输出端最大输出端最大循环循环时间/负扭矩/(Nm）正扭矩/（Nm）次数ms01-20815016 26033302-30020816 26033303-1753241626033309-38353240833310-416591407333最大为5.6 110-4，分别如图11、图12 所示。行星架损伤寿命安全裕度较大，无发生疲劳破坏风险。0.2460.0330.00330.000333.3e-0063.3e-006损伤寿命：2.9 8 10-1e-029图11行星架高周疲劳损伤寿命云图损伤寿命：0.0 0 0 12 30.2460.0330.00330.000333.3e-0052.8e-005损伤寿命：5.6 110 1e-029图12行星架后端板高周疲劳损伤寿命云图6.2试验验证对NW型行星架结构进行高周扭转疲劳试验，试验所用工况和试验大纲要求保持一致。试验开始前，对减速器进行正反转各6 0 s润滑试验，再对减速器按试验工况在高周扭转疲劳试验台架上进行扭转疲劳试验。高周扭转疲劳试验台架如图13 所示。试验完成后，对减速器花键、齿、轴及行星架等进行探伤检测。检测结果显示以上零部件均未出现目视可见的磨损和塑性变形，机械结构探伤检测无裂纹，182023,61(11)总第7 11期机械制造罡（编辑n6时问间？专题报导图13高周扭转疲劳试验台架密封结构均未出现损坏，机械件未出现运转损坏。可见，NW型行星架结构经高周扭转疲劳试验后，满足试验检测要求。用有限元法对优化后的NW型行星架结构进行疲劳计算，发生失效风险较小，与试验结果高度吻合，说明用有限元法对前期行星架结构设计具有切实可行的指导意义。7结束语笔者基于有限元法对NW型行星架结构进行仿真校核，根据实际工作需求，从应力、行星齿轮轴相对变形量、最大螺栓利用率三个方面进行评估。通过仿真，可以看出NW型行星架结构可靠性优良，但存在行星齿轮轴相对变形量较大及最大螺栓利用率不足等风险项，这将使NW型行星架结构在工作中有产生齿轮敲击及松脱等风险。针对识别的风险，对NW型行星架结构进行优化，并与竞品机进行对比分析。优化后行星齿轮轴最大变形量优于竞品机，最大螺栓利用率满足目标要求。组织评审，一致认为优化后的NW型行星架结构在运行过程中产生风险较小。最后对优化后的NW型行星架结构进行疲劳损伤寿命校核，得到发生疲劳失效风险较小，与后期试验结果高度吻合。用有限元法对NW型行星架结构在设计初期进行校核，可以有效规避存在的设计风险，为后期优化提供方向，并且可以有效节约设计成本，能够为后期试验提供支撑。参考文献1饶振纲.行星齿轮传动设计M.2版.北京：化学工业出版社,2 0 14.2吴迅.基于有限元应力分析下对差动减速器行星架设计优化J.中国设备工程，2 0 2 0（9）：9 0-9 1.3赵丽娟，陈令国，刘红梅。矿用减速器行星架的有限元分析J.煤矿机械，2 0 0 7,2 8（1）：5 1-5 2.4孙黎，石鹏飞，代海涛，等.风电齿轮箱低速级行星架极限强度分析J.机械制造，2 0 17,5 5（8）：43-46.5李辉云.行星齿轮减速器仿真及有限元研究D.天津：河北工业大学，2 0 13.6涂杰，罗强.大型风电齿轮箱行星架一销轴系统有限元分析J.机床与液压,2 0 18,46(2 1)：16 9-17 2.7徐恺，熊永红，刘威力，等.NCW行星轮系虚拟样机研究J.机械与电子,2 0 0 9,2 7（4）：7 4-7 5.8王成.有限单元法M.北京：清华大学出版社，2 0 0 3.9石亦平，周玉蓉.ABAQUS有限元分析实例详解M.北京：机械工业出版社，2 0 0 6.10许金泉.疲劳力学M.北京：科学出版社，2 0 17.11 拉兰内.疲劳损伤M.张慰，译.北京：国防工业出版社，2021.收仅何时尚：2 0 2 3 一00作者简介：杜佐飞（19 8 7 一），男，工程师，主要研究方向为机械结构可靠性及疲劳耐久性。格劳博中国开放日盛大举行日前，格劳博机床（中国）有限公司在大连生产基地举办了盛大的开放日活动。这是违三年后格劳博中国首次举办的大规模公开展示活动，也是格劳博中国第二工厂和全新应用技术中心落成后的首次公开亮相。以“德国顶尖机床一一源自格劳博中国”为主题，本次开放日在全面展示格劳博中国自主生产能力的同时，通过典型零件现场切削形式直观呈现格劳博设备、自动化、系统业务的高效创新，携手制造生态合作伙伴企业，聚焦新能源汽车、医学工程、航空航天、机械制造、工具模具等多行业解决方案。本次开放日活动吸引了用户企业代表、合作伙伴、研究机构等约2 3 0 位嘉宾参加。格劳博中国表示，将全价值链的制造环节都集中在中国，是格劳博集团加码投资中国的又一个证明。做“最中国的德国公司”，在大连，用户可以获得与德国明德海姆总部一样的体验，甚至在某些层面更加灵活、更加高效。扎根中国、服务好中国用户是格劳博集团的一项长期战略，这一愿景于2 0 12 年随着格劳博机床（中国）有限公司在大连正式成立而得以落地生根。站在第二个十年的起点，格劳博中国已历练成为一家具备完全自主生产实力的本地化公司，为中国乃至亚洲用户提供从设计、研发，生产、装配、工艺调试，再到技术支持、售后等的一系列服务。(本刊)