飞机柔性件装配中的偏差分析范文

《机械设计与制造工程杂志》2014年第六期

1柔性件铆接装配过程建模

目前，关于柔性件装配偏差分析的相关研究中，最具代表性的是基于装配过程柔性件线弹性变形假设的影响系数法［4］(MethodofInfluenceCoeffi-cient，MIC)。其基本思想是通过力学模型建立柔性件局部受力与变形的线性关系，然后运用该模型进行装配偏差统计分析，该方法的显著优点是在效率上远远优于直接运用有限元分析工具建立仿真模型进行偏差分析的蒙特卡罗仿真法［1］。关于柔性件装配过程的偏差来源，影响装配准确度的主要因素包括以下几方面:偏差(零件偏差、过程偏差)、装配方案(装配顺序、夹具方案、连接配置)和概念设计(几何体形状、刚度条件、零件数量、接头类型)［5］。

1．1基于线弹性假设的装配偏差分析线弹性假设方法是基于柔性件小变形、线弹性的假设，运用线性有限元的方法研究柔性件在装配过程中装配力与装配变形之间的关系，建立反映装配前后柔性件关键测点偏差之间关系的力学模型。飞机柔性件装配的过程中，由于零件制造偏差、装配夹具偏差以及定位偏差的耦合，相互具有连接关系的零件在装配点会形成一定的位置偏差。为了完成装配，需施加一定的装配力，强制柔性件变形，使装配点的偏差闭合。零件变形可以分为弹性变形和塑性变形，弹性变形是线性的，塑性变形是非线性的。小变形情况下，可认为零件只有弹性变形，这时位移与受力的关系是线性的。柔性件在装配过程中的变形是小变形，可以假设为线弹性的过程，刚度矩阵K不变。在这种假设的前提下，可以运用线性关系式表达力F与变形量δ之间的关系，如式(1)所示。有限元模型的总体刚度矩阵是大规模的稀疏矩阵，不利于计算。本文根据线性分析各阶段的柔性件定位以及约束加载边界条件，利用商业软件ANSYS提取超元刚度矩阵，超元刚度矩阵只包含对装配偏差起作用的关键节点及自由度，大大简化了矩阵，提高了计算效率。

1．2柔性件铆接装配过程分解飞机柔性件铆接，首先采用夹具对工件进行定位夹持，然后用铆枪(铆机)将两柔性件上的铆接点定位到理论连接位置并铆接，随后释放铆枪(铆机)，两板回弹至平衡位置。根据装配中不同的力作用过程可将柔性件的铆接过程分为以下4个步骤。步骤1:零件定位，即参与装配的零件各自定位在夹具上，如图1(a)所示。零件存在制造误差，图中零件在关键点处的制造偏差分别为V1和V2。步骤2:铆枪把薄板零件夹持到名义位置，如图1(b)所示。铆枪对两个薄板分别施加铆接力F1和F2，则有步骤3:铆接两薄板零件，如图1(c)所示。两柔性件通过铆接的方式连接在一起。连接后装配体的总刚度矩阵变为KZ。步骤4:释放铆枪及部分夹具，装配体从名义位置回弹至平衡位置，如图1(d)所示。于是，回弹力FH与回弹位移U之间有以下关系:FH=KZU。回弹力大小应等于两板的定位夹紧力和铆接力，按照以上对装配过程的分析，可建立装配误差与零件输入误差之间的线性关系，通过预测装配体的质量来评价零件制造误差和设计方案是否符合制造准确度要求，提高了仿真分析效率。

2基于线弹性假设的接触计算分析

铆接装配中的接触问题属于边界非线性问题。由于基于线弹性假设的柔性件装配过程建模难以将接触力影响直接纳入力与位移的简单线性关系中，相关研究也很少考虑装配中柔性件的相互接触问题，然而在装配存在较大零件干涉的情况下，薄板件可能出现相互穿透这种与实际装配情况相悖的情况，严重影响分析结果的准确性和可靠性。铆接装配工艺通常要求连接件在指定连接位置下形状闭合，即需求解给定铆接位置条件下的装配力。铆接完成后，进行接触计算，将参与装配的两个柔性件分为上侧面和下侧面，并规定上侧面的节点不可穿透下侧面的单元，以约束两柔性件的干涉行为。接触计算主要解决装配中接触状态分析、接触力与节点位移关系、接触平衡计算等问题。基于接触的装配偏差分析流程如图2所示。

2．1装配中动态接触状态分析基于三角网格模型进行几何干涉分析，通过计算网格模型中上侧板节点与下侧板单元的距离矢量，判断上侧板节点干涉情况及干涉量。如图3所示，已知上侧面节点P在笛卡尔直角坐标系中的坐标，下侧面相应的最近的单元ΔABC的法矢n，节点A，B，C的坐标，那么可知上侧面节点与相应下侧面单元的距离d=AP•n，其中AP为以A为起点、P为终点的向量。当d＞0时，节点P处于ΔABC的法矢正方向一侧，未发生干涉;当d=0时，节点P处于ΔABC面内，未发生干涉;当d＜0时，节点P处于ΔABC的法矢负方向一侧，发生干涉。

2．2节点间接触力的等效转化上侧面节点与相应下侧面单元的相互作用力应满足力与力矩平衡条件。如图4所示，P''''为上侧面节点P在下侧面单元ΔABC的投影点，分别以ΔABC三条边为转矩，按力与力矩平衡条件分别求出ΔABC的各节点所受的接触反力。设P点受单元ΔABC作用力为R，ΔABC顶点A，B，C受到。

2．3节点干涉量与接触力的线性作用关系发生接触的两个物体必须满足3个条件［6－8］，即:(1)单边接触条件，两薄板之间间隙必须大于等于零，即d≥0;(2)压力条件，两薄板之间相互作用力必须为负，即为正压力，p≤0;(3)补充条件，受压力的节点，与对应单元处于恰好接触状态，即p×d=0。为了消除已经产生的干涉，在干涉节点上施加接触力，使之回到正常位置，那么选择节点的次序、施加接触力的大小以及判定接触状态正常与否的条件将是接触算法的重点。在小变形的装配过程中，零件的几何形状、夹具布局、零件的节点和单元间的相互位置关系近似不变，为简化接触问题建模和提高效率，不妨认为接触位置迭代过程中节点受力与节点间隙的作用关系保持不变。若上侧面可能发生接触的区域共n个节点，设任意节点i受到接触力Ri作用，该区域受到接触力RUp(RUp={R1，R2，…，Rn})作用，相应的下侧面接触区域单元中各节点受到等效接触力RDown作用，则RDown与RUp满足力与力矩平衡条件。在两柔性件夹紧、铆枪压紧以后，两薄板在接触区域发生干涉，设第j个节点与相应下侧面单元的初始间隙为d0j，在该节点上施加单位力，相应下侧面单元也将产生力的响应，从而引起上侧面n个节点与相应。

3接触平衡位置的迭代算法

接触平衡算法以接触区域节点干涉量与接触力的线性关系为依据，以单边接触条件、压力条件和补充条件为约束，通过动态地调整各节点的接触力和节点位置，多次迭代求解后得到零件平衡位置，平衡后的零件间接触力也是后续回弹力计算的基础。应用接触算法解决干涉问题的基本思路是通过对接触区域存在最大干涉量的节点施加接触力来消除干涉，同时更新网格节点位置，再按照接触的条件进行检测，若仍存在干涉则重复上述过程，通过多次迭代求解接触力和接触平衡位置。接触算法流程如图5所示。初始阶段，接触力向量R清零(即两板不存在接触力)。首先，寻找上侧面接触区域中发生干涉最严重的节点i，即最大穿透点，dmin=di。在该节点施加绝对值尽量小的接触力Ri，其大小刚好使接触区域某一节点的干涉消除。Ri被记录和累加于接触力向量R中。根据式(6)更新两板之间的间隙向量d。检验上侧面所有存在接触力的节点，是否间隙为非正(条件3)。如果条件不成立，则记录下不满足条件的节点号m，并在i点和m点分别施加一定的接触力改变量ΔR(与R的正负相反)。为了避免接触力向量中出现拉力，必须检验所有接触力是否均为非正值(条件2)。如果不成立，则需要将不满足条件的节点间接触力设置为零。然后重新计算间隙向量，直至满足条件2和3。最后，检验间隙向量d是否全部为非负(条件1)，假如满足条件，则算法结束，否则回到第一步寻找发生干涉最严重的节点，重新循环。以上算法仅适用于铆接力和穿透较小的情况，这是因为只有在小变形情况下，线性响应假设才能成立，式(6)才可用。存在较大装配变形和干涉的情况下，上述算法依据的线性响应方法可能导致算法的不可靠或难以收敛。解决方法是将铆接装配过程分为若干小变形阶段，然后采用部分析因设计方法来建立不同阶段的节点作用力和间隙的线性关系，这种分段线性逼近的方法可以更好地体现接触的动态非线性特点。

4实验验证

为了验证本文提出的接触力建模方法，采用VisualC++6．0和OpenGL技术实现了上述接触力建模和装配偏差分析算法，并在柔性件装配实验台上进行了装配偏差分析实验。一方面，通过软件分别在考虑接触和不考虑接触两种情况下进行柔性件装配仿真，记录关键节点的位置信息。另一方面，在实验平台上进行装配实验，测量并记录零件上关键点位置信息，如图6所示。该装配实验台固定两薄板一端，铆接区域两侧采用过约束的夹头压紧，在两板重叠区域进行铆接。零件制造误差、夹具定位误差引起装配体在装配耦合位置的偏差，该偏差在夹具和铆枪的作用下被消除。夹具和铆枪释放后，装配件在回弹力的作用下到达最终平衡位置。图7所示为根据该装配过程建立的有限元网格模型及其约束，两个柔性件采用相同的定位方式，分别有两个固定点(用▲标示)，限制x，y，z向自由度;3个定位夹紧点(用△标记)，限制z向自由度;两个铆接点(用●标记)，限制6个方向自由度。定位夹紧点在回弹阶段被释放。两柔性件的材料属性为:弹性模量E=72GPa，泊松比v=0．3。采用shell63三角形单元(ANSYS单元类型)。两薄板尺寸均为340mm×300mm×0．85mm，各包含158个单元和87个节点。采用vmax/L=3．7/340的实验情况。在考虑接触和不考虑接触情况下分别进行薄板装配仿真，两柔性件装配回弹后的状态如图8所示。可见，在不考虑接触情况下，零件间发生了明显的穿透现象。对比装配仿真与装配实验回弹后的节点位置信息，结果见表1。由表中可知，采用仿真分析时，考虑接触情况下节点位置相对于实验结果的平均偏差为－0．176%，而不考虑接触情况下的节点偏差相对于实验结果的平均偏差为－6．711%。由此可见，考虑接触时装配偏差分析结果比不考虑接触的装配偏差分析更接近实际的装配结果，接触现象对于装配仿真过程的精确度有很大影响。

5结束语

本文提出的飞机柔性件铆接中的接触力建模和装配偏差分析方法，着重解决基于几何干涉分析的接触状态检测、接触平衡求解算法，使得仿真过程更加接近真实装配过程。但是目前的研究仅通过较为简单的柔性平板件验证了飞机柔性件装配接触力模型，要提高该模型的适用性，还需要采用更加复杂的飞机曲面零件为研究对象，拓展接触力模型的应用范围。

作者：张雪谭昌柏王志国单位：南京航空航天大学机电学院