基于蒙特卡洛模拟的某机型前...作动筒运动行程公差控制研究_张丽娟.pdf

教练机TRAINER教练机 2022.NO.40 引 言前护板作动筒是控制飞机护板收放的重要液压系统件，一般由外筒、活塞杆、活塞、端盖等多个零件组成，在飞机起飞、降落时通过液压系统驱动活塞杆伸出、缩进从而起到开闭护板的作用。作动筒的运动行程指的就是活塞杆从极限缩进位置运动到极限伸出位置时的总长度，它直接影响了舱门的开闭运动的可靠性，是作动筒重要的性能指标。但是，前护板作动筒装配结构复杂，其运动行程受多个零件结构特征影响，每个零件在制造过程中存在一定的制造公差，所有制造公差在装配后累积在一起则极易导致行程超差。解决这种由零件制造公差累积形成的装配超差问题，传统的方法为预装配法，在不安装橡胶件的状态下对作动筒进行装配，随后测量运动行程，根据测量结果返工零件，再进行正式装配。某机型前护板作动筒是典型的滚珠锁作动筒，其整体结构相对复杂，在某一批次56件作动筒装配工作中发现有11件在预装配时行程超差，合格率约80.4%，返工率较高，导致批生产的装配周期难以控制。因此，本文选择其为研究对象，收集每个零件中影响行程的相关尺寸的质量数据，随后研究这些尺寸的分布规律并运用蒙特卡洛模拟法对最终装配后行程的尺寸波动情况进行模拟，找出超差原因并实施改进方案，从而大幅度提升该作动筒的一次装配合格率。1 蒙特卡洛模拟法蒙特卡洛模拟法也称为随机模拟法，它是一种以概率论和数理统计为基础，通过对随机变量的统计模拟试验来获得问题近似解的方法。它的主要步骤是：基于蒙特卡洛模拟的某机型前护板作动筒运动行程公差控制研究张丽娟，黄振孜，李 尧（航空工业洪都，江西 南昌，330024）摘要：本文选择某型教练机前护板作动筒作为研究对象，测量分析相关零件特征尺寸分布情况后，运用蒙特卡洛模拟法对最终运动行程的尺寸分布情况进行模拟，根据模拟结果实施工艺改进和工程验证，从而实现作动筒运动行程公差控制。关键词：蒙特卡洛模拟；作动筒；装配公差；运动行程Research on Tolerance Control of Nose Fairing Door ActuatorTravel of a Certain Aircraft Based on Monte Carlo SimulationZhang Lijuan,Huang Zhenzi,Li Yao(AVIC HONGDU,Nanchang,Jiangxi,330024)Abstract:This paper chooses the nose fairing door actuator of a certain type of trainer aircraft as the research object.After measuring and analyzing the characteristic size distribution of the relevant parts,the sizedistribution of the final travel is simulated by using Monte Carlo simulation method.According to the simulation result,process improvement and engineering verification are carried out to achieve the tolerance controlon the actuator travel.Key words:Monte Carlo simulation;Actuator;Assembly tolerance;Travel66教练机 2022.NO.4首先建立简单且便于实现的概率分布模型，分布模型的参数正好对应所求问题的解，再对分布模型进行足够数量的随机模拟，最后建立各种统计量的估计，获得所求解的统计估计值。2 某机型前护板作动筒的运动行程公差分析2.1影响因素分析对作动筒装配件的结构形式及装配过程分析表明，最终的运动行程受到上端盖、活塞杆、外筒、轴套等零件的结构特征影响，装配关系简图见图1，其中，A1为外筒有效长度，A2为两侧调节套的总长，A3为杆接头长度，A4为活塞杆止动面到内孔端面的长度，A5为轴套长度。图1装配关系示意图2.2函数关系建立根据装配公差尺寸链计算的方式可得出最终装配行程与各相关特征的函数关系：行程L=A1+2*A2+A6-A3-A4-A5，最终行程L的尺寸及公差要求为：L=101.10.5。2.3质量数据收集对每个零件抽样40件进行测量，根据质量数据形成如表1所示的各特征的尺寸分布情况，其中尺寸及制造要求是指工艺人员对相应特征提出工艺制造尺寸及公差要求，而尺寸分布情况是根据质量数据绘制的尺寸分布直方图来体现。从表1可以看出，所有的尺寸分布均接近正态分布，其中A1、A4的分布并不是预想的“以极限尺寸平均数为均值的正态分布”，其尺寸误差分布中心相对偏离了公差带中心。2.4尺寸分布模型分析从表1的尺寸分布图形可以大致看出，各特征的尺寸接近正态分布，为了进一步验证尺寸分布的正态性，选用了W检验法进行验证。W检验法全称为Shapiro-Wilk检验，是一种常用的数据正态性检测方法，最终统计量W越接近1则说明数据与正态分布的拟合程度越高，W超过规定的分位数W则认为数据符合正态分布。以A1为例，对样本数据进行分析计算，其中X（k）代表将数据由大到小排列好后第k个数值，a（k）则是计算统计量W需要的系数，其中a（k）可以通过查询W检验法相关的系数表得出。表1各特征的尺寸分布情况尺寸及制造要求尺寸分布情况尺寸及制造要求尺寸分布情况A1（152.50.2）A4（410.2）A2（15.3+0.150）A5（23.50-0.2）A3（320-0.2）A6（14.50.1）L=20k=1a（k）（x（n+1-k)-x（k）0.46640i=1（xi-x?）20.23W=L2i=140（xi-x?）20.942查表可得，在n=40，=0.05的情况下，W0.05=0.940，0.9400.9421，故可以认为A1的尺寸分布近似服从正态分布，相应运用W检测法可以计算出其余几个尺寸的分布情况，同时可以根据样本参数计算出其分布的均值与标准差，计算结果见表2，可以看出，表中的尺寸分布均近似服从正态分布。表2 W检验法计算结果及关键参数计算结果尺寸A1A2A3A4A5A6W0.9420.9600.9680.9660.9570.942均值152.615.431.9241.0323.4214.53标准差0.0770.0070.0370.0840.0380.032.5利用蒙特卡洛模拟进行仿真根据验证分析的结果并结合表2中的数据，将相关信息输入Oracle Crystal Ball软件中即可进行蒙特专题研究67教练机TRAINER教练机 2022.NO.4卡洛模拟。模拟次数设定为50万次，最终模拟情况见图2，可以看出，模拟的结果是其行程保持在101.10.5以内的概率约76.76%，而根据56件产品中11件行程超差的情况计算出行程满足要求的概率约为80.4%，与仿真结果相近。图2蒙特卡洛模拟情况3 基于仿真情况的工艺改进及验证3.1工艺改进根据仿真结果可以看出，公差累积会导致行程超差，且超差均为行程偏长，因此，根据尺寸链计算方式，可以通过减小增环尺寸或增加减环尺寸实现缩短行程。同时考虑到可制造性，应当尽量选择可调整空间比较大的特征进行公差带调整，以满足装配要求。从表1及表2中可以看出，A2的标准差明显小于其余5个尺寸的标准差，而制造公差带要求却与其余几个尺寸相差不大，进一步计算后发现，其余几个尺寸的公差带基本控制在4至6范围内，近似满足公差带设计的原则，而A2尺寸的制造要求的公差带接近20，可改进空间较大。深入分析后发现，由于A2特征对应的零件壁厚较薄，只能通过磨削加工保证厚度，而平面磨床属于高精度加工设备，其理论加工精度为6级，加工能力远超过工艺给定的公差要求，需要进一步确定调整A2加工公差的合理性。最终，根据公差带大小应不小于6值的原则，调整后的A2尺寸要求为15.330.03，即改进后的A2尺寸应当近似服从均值为15.33，标准差为0.007的正态分布。将改进后的方案通过蒙特卡洛进行验证，其验证结果见图3，可以看出，改进后的行程合格率接近97%，相对改进前有明显改善。3.2工程验证实施改进后，共装配了48套作动筒，一次装配行程超差2件，合格率为95.8%，相比改进前的80.4%的合格率有了明显改进。图3改进后的蒙特卡洛模拟4 结 语蒙特卡洛模拟法可以模拟出多个零件尺寸特征相互影响下的最终作动筒运动行程公差的分布情况，借助这种方法能够有效地对装配公差进行预测并实施控制。相对传统方法，它提高了产品一次交检合格率，减少了返工率。该方法可对解决类似液压系统件的装配公差控制问题提供参考。参考文献1 王锋.基于蒙特卡洛模拟法的动车组裙板装配优化研究J.中国设备工程,2020（6）：2.2 熊美，胡莹.蒙特卡洛仿真法在飞轮装配公差设计中的应用J.技术新工艺，2017（6）：4.3 赵东平，张辉，郭津呈.基于关键特性的飞机装配公差分析方法J.西安航空学院学报，2020，38（1）：6.4 赵嫚，茅健.面向制造过程的装配公差模型研究J.工程设计学报，2014，21（6）：5.（收稿日期：2022-10-26）作者简介张丽娟，女，1988年5月出生，2011年毕业于沈阳航空航天大学，工程师。68