一种基于敏感因子加权的机构动力学仿真方法.pdf

第 44 卷 第 3 期航天返回与遥感 2023 年 6 月 SPACECRAFT RECOVERY&REMOTE SENSING 51 收稿日期：2022-10-20 基金项目：国家自然科学基金（62227812）引用格式：于婷婷,连华东,杨柳,等.基于敏感因子加权的机构动力学仿真方法J.航天返回与遥感,2023,44(3):51-61.YU Tingting,LIAN Huadong,YANG Liu,et al.A Simulation Method of the Dynamic Characteristics for Space Optical Mechanisms Based on Weighted FactorsJ.Spacecraft Recovery&Remote Sensing,2023,44(3):51-61.(in Chinese)一种基于敏感因子加权的机构动力学仿真方法 于婷婷 连华东 杨柳 赵鑫 梁凤超 范建凯 林喆（北京空间机电研究所，北京 100094）摘 要 为了解决批产高精度调焦机构有限元分析与扫频试验结果存在偏差的问题，提高空间复杂机构设计能力，文章提出了一种基于敏感因子加权系数的机构动力学有限元分析法。该方法基于 Hertzian接触理论建立了运动副连接的等效刚度模型，确定了敏感因子及其刚度贡献的关键方向，根据刚度与基频关系以及相关试验数据，为模型中敏感因子分配了加权系数，实现了机构的高精度仿真。利用该仿真方法对某空间相机高精度调焦机构进行仿真分析，并与扫频试验进行对比，结果表明前三阶基频振动特性一致且误差均小于 2.8%，验证了该方法的有效性和计算效率，可为高精度空间机构的仿真分析提供一定参考。关键词 空间机构 刚度 加权 有限元分析 振动特性 空间相机 中图分类号:V423.9 文献标志码:A 文章编号:1009-8518(2023)03-0051-11 DOI:10.3969/j.issn.1009-8518.2023.03.006 A Simulation Method of the Dynamic Characteristics for Space Optical Mechanisms Based on Weighted Factors YU Tingting LIAN Huadong YANG Liu ZHAO Xin LIANG Fengchao FAN Jiankai LIN Zhe（Beijing Institute of Space Mechanics&Electricity,Beijing 100094,China）Abstract To reduce of the deviation between analysis and test results and improve the design performance for space complex mechanisms,a dynamic finite element analysis method for space focus adjustment mechanisms is proposed based on sensitive factor weighted coefficients.This method establishes an equivalent stiffness model of the motion pair joints based on the Hertzian contact theory,and clarifies the key direction of the sensitive factor and its stiffness contribution.According to the relationship between stiffness,the fundamental frequency and relevant test data,the weighted coefficients of the sensitive factors in the model are determined,thus achieving high-precision simulation.Using this simulation method,the performance of a high-precision focus adjuststment mechanism in a space camera is analyzed with and without weighted coefficients,and the results are compared with the sweep test ones.The comparisons show that for the first three fundamental frequencies,the vibration characteristics are consistent and the fundamental frequency error is less than 2.8%,verifying the effectiveness of the method and providing a reference for the simulation analysis of high-precision space mechanisms.52 航 天 返 回 与 遥 感 2023 年第 44 卷 Keywords space mechanism;stiffness;weighted;finite element analysis;vibration characteristics;space camera 0 引言 空间光学相机在轨工作寿命期内，由于发射阶段过载冲击、在轨环境温度变化等因素的影响，会产生离焦从而影响成像品质。调焦机构作为当前空间相机中重要组成部分之一1，能够校正偏离的焦面，使相机的光学性能满足成像要求。目前调焦机构通常采用固体润滑角接触球轴承实现结构支撑，利用步进电机进行驱动，并通过滚珠丝杠副及导轨滑块组件实现定位导向。由于机构连接部组件复杂，运动副连接中存在间隙、摩擦等现象，且预紧力直接影响刚度的变化，导致活动机构在有限元分析时，不能像对简单结构件仿真一样逐一处理接触应力或等效刚度等，如此造成活动机构的有限元分析工作量巨大、程序运算慢，甚至计算出的结果与实际力学特性差距很大，不能为机构的设计迭代提供指导。尤其空间相机中的机构负载均是光学元件，机构有限元分析的不准确则不能为光学元件提供有力保障。因此，如何高效准确的模拟预测机构刚度特别是运动副连接刚度，并快速准确的进行有限元处理，对于高精度空间运动机构的设计与力学特性分析具有重要意义。针对复杂机构刚度的分析方法已开展了许多相关研究，如 Gosselin C2建立了并联机构工作空间内的刚度映射和条件数映；赵铁石等3基于影响系数法，并借助虚功原理，建立了包含弹性变形及刚度连续变化过程的并联机构连续刚度非线性映射通用模型；李嘉等4利用有限元技术，针对两类移动副和一类球副的不同结构组合，分析了结构参数对其刚度性能的影响。上述研究从理论方面扩展了刚度模型的分析方法，其中许多研究针对刚度模型采用了大量的分析计算及复杂的有限元分析手段。但在工程实践中，特别是针对高精度、部组件连接关系复杂的空间精密机构，高效准确完成优化模型的设计，以简单的方式建立连接副刚度的近似模型，并给出可靠准确的仿真结果，具有重要的工程意义。本文针对光学遥感相机中一维直线调焦机构的动力学特性开展了有限元分析研究，通过判断刚度模型敏感因子，分析敏感因子的主要贡献刚度方向，并采用加权系数分配法5-6将不敏感因子的刚度贡献等效至敏感因子刚度，将复杂机构的刚度模型简单化，基于刚度参数创建仿真连接关系，提高了仿真模型的准确性；最后通过试验对仿真模型进行了验证，结果表明增加权系数后的仿真结果精度大大提高。该方法提高了复杂活动机构设计迭代优化的效率，高效准确的仿真结果为后续试验验证奠定了有力基础，可应用于同类传动方式的运动机构的设计与动力学特性研究。1 机构敏感因子刚度分析 刚度是影响机构动态性能的重要评价指标，本文以某空间遥感相机中的调焦机构为例，基于对机构运动副的刚度分析，开展了对机构的有限元分析研究。如图 1 所示，调焦机构由步进电机通过减速齿轮副驱动精密微型滚珠丝杠实现旋转变直线运动，其中滚珠丝杠副主要承受机构运动方向的轴向力；丝杠两端使用面对面安装角接触球轴承支撑并轴向定位预紧，增加机构轴向承载力；为实现高精度一维导向运动，机构中采用了一对可承受法向载荷的 HSR 导轨滑块组件。根据对相同构型机构的试验数据整理发现，前 3 阶扫频试验的运动方向均相同，分别表现在机构的轴向、径向和法向（即图 1 中的 X、Y、Z 方向），且关键因子刚度的变化直接线性影响机构的前 3 阶基频。调焦机构中主要通过滚珠丝杠副和轴承承受 X 向载荷，导轨滑块组件和轴承同时承受 Y 向载荷，但导轨滑块主要承受 Z 向载荷，因此可以确定在 X、Y、Z 三个方向主要贡献刚度的因子为滚珠丝杠副、轴承和导轨滑块副，并分析出了 3 个敏感因子的主要贡献方向。由于活动机构在不同状态下的刚度变现不 第 3 期 于婷婷 等:基于敏感因子加权的机构动力学仿真方法 53 同，研究方法为了保证机构状态一致，机构力学试验模拟发射过程进行了位置锁定，此时可认为机构刚度保持不变。在机构的有限元分析中，调焦机构在 X、Y、Z 向的刚度 KX、KY、KZ是关键指标，直接影响机构力学试验结果。其中，KX主要取决于滚珠丝杠副及角接触球轴承系统的刚度，KY取决于导轨滑块及角接触球轴承系统的刚度，KZ取决于导轨滑块系统的刚度。因此基于 Hertzian 接触理论5-6，在假定材料均匀、各向同性且完全弹性，以及接触表面的摩擦力忽略不计的条件下，结合机构中敏感因子的各方向刚度的串并联关系，并忽略结构主体结构刚度影响，可近似得出调焦机构在 X、Y、Z 向的刚度 KX、KY、KZ为 111121211111212122223232333322222XYZKKKKKKKKKK+=+=（1）式中 K11为轴系滚珠丝杠副轴向刚度；K21为轴承轴向刚度；K22为轴承径向刚度；K23为轴承法向刚度；K32为导轨滑块径向刚度；K33为导轨滑块法向刚度；11为滚珠丝杠副刚度权系数；21为轴承轴向刚度权系数；22为轴承径向刚度权系数；32为导轨滑块径向刚度权系数；33为导轨滑块法向刚度权系数。系统刚度模拟如图 2 所示。图 1 调焦机构模型示意 图 2 系统刚度模拟示意 Fig.1 Model of the focus adjustment mechanism Fig.2 System stiffness matrix 对系统各敏感因子进行刚度分析计算时，虽然忽略了其他结构件的刚度，但为提高系统刚度模型的计算精度、保证仿真的准确性，需将不敏感因素对系统力学特性的影响，通过刚度参数拟合到敏感因子权系数 中；确定敏感因子后，对 3 个敏感因子关键方向的刚度展开求解及权系数计算。由于敏感因子结构内均包含滚珠，刚度计算复杂，因此在依据 Hertzian 接触理论的基础上还参考了近年较为准确而简便的算法6-7。1.1 滚珠丝杠副刚度 调焦机构中承载轴向载荷的活动副包含滚珠丝杠副与轴承，其中滚珠丝杠副的预紧方式一般有 3 种：双螺母垫片预紧（通过改变两螺母间垫片厚度调整预紧力）、单螺母变位导程预紧（通过在导珠管之间改变某一导程施加预紧）、单螺母增大钢球预紧（通过大直径钢球调整预紧）。图 3 为滚珠丝杠副综合性能退化曲线，可以看出：预紧力越大，接触刚度增加，定位精度高，伴随寿命降低。因此在调焦机构关键运动副设计中，综合考虑产品力学性能和服役寿命要求，保证钢球固体润滑 MoS2膜层的寿命，设计滚珠与滚道存在 5 m 间隙，理论是无预紧状态，但实际装配后，由于滚道与钢球存在尺寸偏差，会出现微预紧力，存在图 4 所示的 2 点式接触。54 航 天 返 回 与 遥 感 2023 年第 44 卷 图 3 滚珠丝杠副综合性能退化曲线示意 图 4 钢球预紧方向示意 Fig.3 Degeneration curves of the integrated performances of the bass screw Fig.4 Directions of the pretighten forces on the balls 根据滚珠丝杠副在机构中的安装方式及受力方式，轴向刚度 K11主要由滚珠丝杠轴向刚度 KS、滚珠螺母轴向刚度 KN、螺母及轴承安装基座轴向刚度 KH三方面组成7。1）滚珠丝杠轴向刚度 KS与安装方式有关，本文调焦机构采用固定-固定的安装方式，因此，在刚度计算过程中必须考虑丝杠螺母到固定支撑点产生的最大轴向位移，KS为7 3S10/()KAELL=（2）式中 A为以丝杠螺纹滚道底径为直径的横截面积（mm2）；E为弹性模量（MPa）；L为两端支撑轴承间距（mm）；为丝杠螺母至固定点间距（mm），当=L/2时产生最大轴向位移。2）钢球与螺纹滚道在2点式接触过程中会产生弹性和塑性变形，在钢球与滚道面承受最大接触应力处，只产生1/10 000钢球直径的塑性变形，所以仅分析接触面发生弹性变形条件下的刚度7，即 3Npra/(0.1)KFC=（3）式中 为行程范围内根据精度等级确定的行程变动量系数，取值范围为0.5,0.6；Fpr为丝杠副的预紧力（Nm）；Ca为额定动载荷（N）。3）考虑滚珠丝杠副及支撑轴承安装基座的轴向刚度，可按照结构件单轴方向受载计算KH，即 HK=F/L（4）式中 F为支撑基座轴向方向承受的载荷。综合上述因素，可以得到敏感因子滚珠丝杠副的轴向刚度K11为7 ()11SNSHHNSHN+/()KK KK KK KK K K=（5）根据式（5），并结合式（2）（4）及滚珠丝杠副手册中相关参数，可求得轴向接触刚度K11=851.32 N/m。由于机构设计过程中，滚珠丝杠副两端通过一组面对面角接触球轴承支撑，并通过一组卸载径向力装置与活动板连接，已经最大限度避免了滚珠螺母承受径向载荷或扭矩载荷的情况，故这里忽略径向、法向刚度K12和K13。1.2 轴承刚度 调焦机构中滚珠丝杠副两端采用了面对面方式的角接触球轴承，如图5所示通过预紧端盖压紧轴承外圈，由于实际装配过程需在安装一端轴承端盖后，才可以对另一端进行预紧，施加轴向力Fa使内圈压紧，轴承内外圈产生不等的相对位移la，两端预紧力转变为Fa1、Fa2，此时轴向预紧力Fa可表示为8 251/23/2aa1a2ba1a2a2a13.462(sin)()max(),0)FFFzDllll=+（6）式中 Db为角接触球轴承滚珠直径（mm）；为角接触球轴承给接触角（）；la1、la2分别为未预紧与预紧端轴承内外圈位移量（mm）；z为钢球数量。第 3 期 于婷婷 等:基于敏感因子加权的机构动力学仿真方法 55 通过轴向位移与轴向力的关系，可计算出轴向刚度K21，即 21aa/KFl=（7）在施加预紧力初始阶段，轴承刚度与预紧力呈正比，与寿命呈反比关系（如图6所示）；之后伴随预紧力的增大，刚度逐渐增大，而寿命逐渐减小，当预紧力到达临界点时，轴承刚度趋于稳定值，寿命加剧减小。图 5 轴承轴向预紧时力与位移关系 图 6 轴承预紧后参数变化关系 Fig.5 Relationship between the bearing force and displacement under axial preloading Fig.6 Parameter change relationship after bearing preload 对于轴向预紧作用下轴承的刚度，已有许多简便算法，针对角接触球轴承的轴向刚度K21及径向刚度K22，本文参考文献8-11进行计算，将轴向预紧力看作施加轴向载荷Fa，从而在轴承内外圈产生相对位移la，得到K21及K22的近似值：251/31/321ba2/321/31/322ba6.866(sin)1.7164cos/sinKzDFKzDF=（8）调焦机构中使用的角接触球轴承为708C轴承，为满足刚度与寿命的要求，采用定压预紧的方式施加小预紧力，对708C轴承单个轴承施加轴向载荷进行计算，通过预紧增大摩擦力矩（23）103 Nm，保证轴承刚度约为 988.7 N/m。对于径向刚度K22，在后续计算中减少变量，因此假定K22保持不变。由于机构中轴承支撑于滚珠丝杠副两端，且滚珠丝杠副与负载间存在间隙无法承载法向力，因此分析中可以忽略轴承法向刚度K23。1.3 导轨滑块刚度 导轨在调焦机构中主要承载径向载荷，可以忽略轴向刚度K31。当导轨副承受径向和法向载荷时，4列滚道中任一钢球都与导轨、滑块形成接触副，具体受力情况如图7所示。调焦机构中支撑负载的径向刚度与法向刚度主要取决于导轨的结构参数和滚珠的弹性变形。在负载水平匀速运动过程中，滑块与导轨在容许载荷范围内产生弹性变形，此时产生的径向载荷与变形量之比即为径向刚度。进行导轨滑块刚度分析时，将导轨外滚道、滑块内滚道与钢球的局部接触区域均视为弹性体，载荷作用在一个小的椭圆形区域，从而产生局部弹性变形。依据THK导轨技术手册提供的关键参数12-16，计算导轨滑块组件的径向刚度K32与法向刚度K3317，具体公式为 图 7 受径向载荷时导轨滑块副受力示意 Fig.7 Load schematic of the rail under radical loads 56 航 天 返 回 与 遥 感 2023 年第 44 卷 11.51/21/21/23233423 2F()/()()KKkDR=（9）式中 为滚珠变形量（mm）；为滚珠与沟槽接触的主曲率之和；22122112(1)(1)/()REEE E=+，其中E1、E2分别为滚珠和导轨的弹性模量（GPa），1、2分别为滚珠和导轨的泊松比；D为椭圆长半轴系数；F(k)为第一类完全椭圆积分；k为椭圆离心率。调焦机构应用的HSR12导轨滑块副包含4列滚动面，各钢球列接触角均为45，拟定在径向与法向具有相同的额定载荷。根据导轨参数，通过式（9）可计算出导轨单个滚珠在径向与法向的刚度，K32=K33=899104 N/m。2 敏感因子加权系数分配 加权分配法考虑到各分系统对整个系统的影响程度，以分系统刚度引起整个系统刚度特性发生变化的概率为依据，来计算各影响因子权系数。章节1中已获得3个敏感因子关键方向的刚度，基于大量相同构型机构的历史试验数据，可计算各因子权系数。但考虑到机构是由多个结构件和高精度运动副装配组成，权系数不仅评价因子的敏感程度，还需要考虑其他不敏感因素的影响刚度，因此本文在采用加权分配法的同时，也将不敏感因子刚度的影响综合到了权系数中，保证整体刚度接近实际值。由于式（1）中未知权系数大于求解方程数，因此先根据机构3个方向的刚度特性确定3个权系数。其中滚珠丝杠副只承受轴向载荷，其轴向刚度值在实际中装配易于控制，误差范围较小，因此定义权系数11=1；轴承通过预紧来增加轴向刚度，对于刚度KY、KZ，轴承对整体结构的径向刚度影响小，因此可以假定权系数22=1。因此，对于X、Y、Z三个方向的不敏感因素刚度误差，将其中的轴向刚度误差等效至系数21，径向刚度误差等效至系数32，法向（Y向）误差等效至系数33。相同构型的调焦机构具有相同的敏感因子。根据基频与刚度影响呈线性关系的条件，对相同构型机构的扫频试验结果及仿真结果进行对比计算，可以获得KX、KY、KZ，将这3个数值和之前计算出的敏感因子刚度值及系数11、22、23等一并代入式（1），求解出21、32、33。经计算，21=1.05，32=1.2，33=1.33。利用上述结果，可得到分配权系数后的敏感因子刚度模型。3 仿真与试验验证 通过分析敏感因子对机构刚度的影响，获得各组件刚度信息，在此基础上建立调焦机构有限元仿真模型，并对机构敏感因子刚度的准确性进行验证，为后续加权系数的确定提供可靠数据。考虑到机构中细节特征复杂、异形结构件较多，几何建模过程中需要花费较多时间进行几何清理和网格划分，为了快速准确发现模型问题及优化模型，本文通过Ansys Workbench建立机构仿真模型，实现针对复杂机械系统的快速仿真验证。在仿真验证过程中，首先利用Workbench自动创建的刚性连接模型得到仿真结果，由于该结果并不完全符合直线机构力学特性规律，因此初始模型不够准确，证明机构内活动部件设置模拟连接副的重要性；通过对敏感因子设置模拟连接副，建立准确的仿真模型，并对有无加权系数的情况进行仿真；最后与试验结果进行对比分析，验证增加权系数的意义。仿真分析时涉及的结构材料属性如表1所示。为了将仿真误差降到最低，仿真模型中将网格的精细程度以及单元阶数的设置，选择10节点四面体单元，保证分析结果符合收敛性要求。第 3 期 于婷婷 等:基于敏感因子加权的机构动力学仿真方法 57 表 1 结构件材料 Tab.1 Material properties 主要结构件 材料名称 弹性模型 E/GPa 泊松比 密度/（gcm-3）机构基座 TC4 113 0.34 2.85 THK-滑块 马氏体/奥氏体不锈钢 184 0.30 7.90 THK-导轨 滚珠丝杠副丝杠 不锈钢 184 0.30 7.90 滚珠丝杠副螺母 角接触球轴承 GCr15 212 0.29 7.81 3.1 仿真模型的建立与优化 根据一维直线运动机构的力学特性可知，机构的第一阶基频表现为X向的位移运动，原因是由于机构沿X向一维运动，X向刚度主要依靠轴承与滚珠丝杠副的轴向刚度；第二阶基频表现为法向方向的运动；第三阶基频表现为径向方向的运动。通过该规律初步验证仿真模型的准确性。首先，通过Workbench中自动创建的刚性连接关系进行简单仿真，验证自动生成连接关系的仿真模型是否准确，初始仿真结果如图8所示。结果表明：第一阶频率发生在机构与负载连接的柔性卸载节处；第二阶频率发生在柔性卸载节中连接杆上，表现为负载的Y向平移；由于整体结构中悬臂结构的特性，在第三阶频率出现了机构整体绕Z向的旋转运动。上述仿真结果中，最明显的一阶基频并不是沿X向的平动，与力学特性规律存在明显差别，因此该仿真模型通过刚性连接的方法误差太大，需要根据运动副刚度及其连接关系，对模型中各部件连接方式进行修改完善。（a）第一阶频率 210.85 Hz（b）第二阶频率 248.37 Hz（c）第三阶频率 253.09 Hz（a）First order frequency 210.85 Hz（b）Second order frequency 248.37 Hz（c）Third order frequency 253.09 Hz 图 8 初始仿真结果 Fig.8 Initial simulation results 综合上述分析，采用模拟运动副方法，用弹簧单元来模拟连接刚度，创建更为准确的连接关系。3个敏感因子的模拟连接关系如下。1）滚珠丝杠副：滚珠丝杠副螺杆的转动变为螺母的平动，在仿真模型中，放开螺母与螺杆轴向转动和轴向平移自由度，并施加轴向刚度约束K11。2）轴承：由于电机驱动齿轮副带动丝杠轴旋转，轴承内圈相对外圈旋转，在模型中释放轴承与基座之间轴向转动和平动约束，分别在X向施加刚度K21，在径向施加刚度K22，模拟轴承的连接刚度。3）导轨滑块组件：作为机构中支撑和导向部件，导轨滑块组件刚度主要影响径向和法向基频。在模型中连接关系设置为：释放X、Y、Z方向平动自由度，并在Y向和Z向分别施加刚度K32和K33。58 航 天 返 回 与 遥 感 2023 年第 44 卷 在建立了模拟连接副并设置敏感因子的刚度K11、K21、K22、K32和K33后求解，仿真结果如图9所示，为更直观的判断机构模态仿真结果的运动方向，图9中将负载与活动板隐藏。（a）第一阶频率 104.5 Hz（b）第二阶频率 126.32 Hz（c）第三阶频率 138.6 Hz（a）First order frequency 104.5 Hz（b）Second order frequency 126.32 Hz（c）Third order frequency 138.6 Hz 图 9 建立连接副后的仿真结果（无加权系数）Fig.9 Simulation result after establishing the connection pair（without weighted coefficients)可以看出，第一阶频率为104.5 Hz，运动表现为活动部件（活动板、滑块及丝杠螺母）沿X向平动；第二阶频率为126.32 Hz，运动表现为机构整体沿Y向的摆动；第三阶频率为138.6 Hz，运动表现为机构整体在Z向摆动。前3阶机构频率特性与一维直线机构的力学特性规律一致，因此可确认仿真模型建立正确，同时也获得了无加权系数情况下的仿真结果。在该仿真模型基础上进行权系数分配，保证后续提高基频分析精度。3.2 分配加权系数后的仿真 基于准确的机构仿真模型，对模拟连接副中的刚度值依次分配加权系数后，再次进行仿真分析验证，验证敏感因子刚度值分配权系数的必要性。由图10可以看到，第一阶频率为115.88 Hz，运动表现为运动部件（活动板、滑块及丝杠螺母）的X向平动；第二阶频率为139.6 Hz，机构整体Y向摆动；第三阶频率为143.8 Hz，机构整体表现为Z向摆动。各阶运动方向依旧与力学特性规律一致。（a）第一阶频率 115.88 Hz（b）第二阶频率 139.6 Hz（c）第三阶频率 143.8 Hz（a）First order frequency 115.88 Hz（b）Second order frequency 139.6 Hz（c）Third order frequency 143.8 Hz 图 10 分配加权系数后的仿真结果 Fig.10 Simulation results with the weighted coefficients 3.3 扫频试验验证 为了验证模拟连接副配合权系数仿真方法的精度，对机构进行了动力学试验，图11中给出了X、Y、Z三个方向的扫频试验结果，以及关键测点的频率响应曲线，可以看到3个方向的基频分别为116.2、142.5、147.9 Hz，试验结果与同构型一维直线机构的力学特性规律保持一致，且曲线中频率点未出现断崖式下跌，说明机构中各组件不存在异常情况，试验结果可靠。表2给出了试验结果与有无加权系数模态仿真结果的对比。可以看出，在考虑了连接刚度的加权系数后，得到的仿真结果与试验结果误差最大不超过 第 3 期 于婷婷 等:基于敏感因子加权的机构动力学仿真方法 59 2.8%，而无加权系数仅考虑连接刚度时的误差最小为6%，最大可至10%。因此，通过分配准确的连接刚度加权系数，在此基础上进行机构仿真计算，可以大大提高机构的仿真精度，为后续机构设计中提高基频、保护负载起到了重要作用。（a）X 向频率响应曲线（a）X-Direction frequency response curve （b）Y 向频率响应曲线（b）Y-Direction frequency response curve （c）Z 向频率响应曲线（b）Z-Direction frequency response curves 图 11 频率响应曲线 Fig.11 Frequency response curve 60 航 天 返 回 与 遥 感 2023 年第 44 卷 表 2 仿真与试验结果对比验证结果 Tab.2 Comparison between simulation and test results 方向 试验结果/Hz 无加权系数 分配加权系数 仿真结果/Hz 误差/%仿真结果/Hz 误差/%X 116.2 104.5 10.0 115.88 0.3 Y 142.5 126.3 11.4 139.6 2.0 Z 147.9 138.6 6.2 143.8 2.8 4 结束语 针对批产调焦机构仿真结果偏差较大问题，为了提高动力学特别是频率仿真精度，以便为设计改进和力学试验提供有力参考，本文利用Hertzian接触理论在对接触区域弹性分析的基础上，分析了滚珠丝杠副、角接触球轴承及导轨滑块等运动副在3个不同方向的刚度。根据同类构型调焦产品试验结果和仿真数据，确定了仿真模型中不同连接刚度所需的加权系数。考虑了连接刚度加权系数后，仿真精度大大提高，与试验结果误差最大不超过2.8%。该仿真方法深度剖析了一维运动机构各部组件的连接关系及对整体的刚度贡献值，提高了复杂活动机构的设计与仿真验证的效率与准确性，可应用于同类传动方式的运动机构的设计与动力学特性研究。参考文献(References)1 谭伟,王殿中,何红艳,等.“高分七号”卫星相机高精度快速在轨调焦方法J.航天返回与遥感,2020,41(2):78-86.TAN Wei,WANG Dianzhong,HE Hongyan,et al.Reserch on High Accurate Refocusing Method of GF-7 CameraJ.Spacecraft Recovery&Remote Sensing,2020,41(2):78-86.(in Chinese)2 GOSSELIN C.Stiffness Mapping for Parallel ManipulatorsJ.IEEE Transactions on Robotics and Automation,1990,6(3):337-382.3 赵铁石,赵延治,边辉,等.空间并联机构连续刚度非线性映射J.机械工程学报,2008,44(8):20-25,32.ZHAO Tieshi,ZHAO Yanzhi,BIAN Hui,et al.Continuous Stiffness Nonlinear Mapping of Spatial Parallel MechanismJ.Chinese Journal of Mechanical Engineering,2008,44(8):20-25,32.(in Chinese)4 李嘉,陈恳,董怡,等.并联柔性铰机器人的静刚度研究J.清华大学学报(自然科学报),1999,39(8):16-20.LI Jia,CHEN Ken,DONG Yi,et al.Study on the Static Stiffness of the Parallel Flexure Joint RobotJ.Journal of Tsinghua University(Science&Technology),1999,39(8):16-20.(in Chinese)5 王乐,李冬,余慕春,等.基于动态加权系数和多目标进化的模型修正方法J.振动与冲击,2020,39(4):284-290.WANG Le,LI Dong,YU Muchun,et al.A Method of Model Updating Based on Dynamic Weighting Coefficients and Multi-Objective EvolutionJ.Journal of Vibration and Shock,2020,39(4):284-290.(in Chinese)6 程鹏辉,吴萌源.空间相机调焦平台的动力学特性研究J.光学精密工程,2019,27(3):602-609.CHENG Penghui,WU Mengyuan.Research on Dynamic Characteristics of Space Camera Focusing PlatformJ.Optics and Precision Engineering,2019,27(3):602-609.(in Chinese)7 JCSGY 精密滚珠丝杠副手册M.北京:北京工研精机股份有限公司,2008.JCSGY BALL SCREWSM.Beijing:Beijing Precision Mechinery&Research CO.LTD.(in Chinese)8 李为民,王海涛.轴向定位预紧轴承刚度计算J.河北工业大学学报,2001,30(2):15-19.LI Weimin,WANG Haitao.Rigidity Calculation of Axial Position Preload BearingsJ.Journal of Hebei University of Technology,2001,30(2):15-19.(in Chinese)9 GUNDUZ A,DREYER J T,SINGH R.Effect of Bearing Preloads on the Modal Characteristics of Ashaft-Bearing Assembly:Experiments on Double Row Angular Contact Ball BearingsJ.Mechanical Systems and Signal Processing,2012,31:176-195.10 OZTURK E,KUMAR U,TURNER S,et al.Investigation of Spindle Bearing Preload on Dynamics and Stability Limit in MillingJ.CIRP Annals-Manufacturing Technology,2012,61(1):343-346.11 杨静文.端盖螺栓拧紧工艺对轴承预紧载荷影响研究D.大连:大连理工大学,2018.YANG Jingwen.Research on Influence of Tightening Technology of Bearing Caps Bolts on Bearing PreloadD.Dalian:Dalian University of Technology,2018.(in Chinese)第 3 期 于婷婷 等:基于敏感因子加权的机构动力学仿真方法 61 12 THK 公司.THK 综合产品目录-B 产品尺寸规格M.出版地不详:THK 公司,2017.THK CO,LTD.T HK Catalog No 511C-B Product SpecificationsM.S.l.:THK CO,LTD.,2017.(in Chinese)13 YI Y S,KIM Y Y,CHOI J S,et al.Dynamic Analysis of a Linear Motion Guide Having Rolling Elements for Precision Positioning DevicesJ.Mechanical Science and Technology,2008,22(1):50-60.14 KATO S,IGARASHI Y.Effects of Ball Groupings on Ball Passage Vibrations of a Linear Guideway Type Ball Bearing,Pitching and Yawing Ball Passage VibrationsJ.Journal of Tribology,2007,129(1):188-193.15 宁星超.滚动直线导轨动态特性研究D.大连:大连理工大学,2016.NING Xingchao.Study on the Dynami Characteristics of Rolling Linear GuideD.Dalian:Dalian University of Technology,2016.(in Chinese)16 安源,杜