考虑区间叠加不确定间隙的空间可展机构动力学研究.pdf

第 40 卷 2023 年第 4 期上海航天（中英文）AEROSPACE SHANGHAI（CHINESE&ENGLISH）考虑区间叠加不确定间隙的空间可展机构动力学研究孟凡斌1，张渤洋2，侯月婷3，王震南4（1.天津大学 精密仪器与光电子工程学院，天津 300070；2.河北工业大学 机械工程学院，天津 300401；3.天津航天机电设备研究所 天津市微低重力环境模拟技术重点实验室，天津 300301；4.哈尔滨工程大学 机电工程学院，黑龙江 哈尔滨 150001）摘要:本文建立了考虑不确定间隙空间可展机构动力学模型，通过对不确定间隙进行描述，结合传统接触模型提出含不确定区间间隙的接触力模型，并将其嵌入展开单元动力学模型，通过切比雪夫区间求解对含区间系统动力学模型进行求解，针对区间求解中启动阶段估计不准、区间叠加误差过大、区间整合的问题与展开单元叠加问题，提出了相适应的区间扩张、局部均值分解与瞬时叠加方法降低区间求解误差，并给出区间间隙累加下动力学响应边界。结果表明：区间间隙大小主要对系统运动精度产生较大影响，而区间间隙数量对系统运动稳定度影响较为严重。区间叠加方法相比于传统拟合方法可进一步提高空间可展机构预示精度，相较于传统拟合结果可提高 7.74%。关键词:可展机构；不确定间隙；多体系统；动力学；信号分解中图分类号:V 414 文献标志码:A DOI:10.19328/ki.20968655.2023.04.012Study on Spatially Deployable Mechanism Dynamics Considering Interval Superposition Uncertain GapsMENG Fanbin1，ZHANG Boyang2，HOU Yueting3，WANG Zhenan4（1.College of Precision Instrument and Optoelectronics Engineering，Tianjin University，Tianjin 300072，China；2.School of Mechanical Engineering，Hebei University of Technology，Tianjin 300401，China；3.Tianjin Institute of Aerospace Mechanical and Electrical Equipment，Tianjin Key Laboratory of Microgravity and Hypogravity Environment Simulation Technology，Tianjin 300301，China；4.School of Mechanical and Electrical Engineering，Harbin Engineering University，Harbin 150001，Heilongjiang，China）Abstract:A spatially deployable mechanism dynamics model considering uncertain gaps is established，and a contact force model with uncertain interval gaps is proposed by describing the uncertain gaps and combining the traditional contact model.The proposed contact force mode is embedded into the unfolding unit dynamics model，and the dynamics model of the system with intervals is solved by Chebyshev interval solving.In view of the problems of inaccurate estimation in the starting stage，excessive interval superposition error，interval integration，and unfolding unit superposition in interval solving，the appropriate interval expansion，local mean decomposition，and instantaneous superposition methods are proposed，and the dynamic response boundaries under interval gap accumulation are given.The results show that the size of the interval gaps has a significant effect on the accuracy of the system motion，and the number of interval gaps has a serious effect on the stability of the system motion.Compared with the traditional fitting method，the interval superposition method can further improve the prediction accuracy of the spatially extensible mechanism，and has a 7.74%increase.Key words:deployable mechanism；uncertain gap；multi-body system；dynamics；signal decomposition收稿日期：20230509；修回日期：20230518基金项目：国家自然科学基金（U21B2075）；中央高校基本科研业务费专项（3072022CFJ0702）；天津市微低重力环境模拟技术重点实验室开放课题（WDZL202102）作者简介：孟凡斌（1983），男，博士，高级工程师，主要研究方向为光电对抗。88第 40 卷 2023 年第 4 期张渤洋，等：考虑区间叠加不确定间隙的空间可展机构动力学研究0引言 空间展开机构1-2高精度展开是卫星天线3、太阳能列阵4等航天部组件稳定运行的前提。但是受设计公差、加工和装配误差，以及部组件运行过程中间隙接触、挤压、摩擦造成的关节变形等5-6因素影响，间隙大小呈现区间不确定现象，造成系统动力学特性难以精确预示的问题。在含区间不确定的多体系统动力学建模过程中，将不确定参数7-9认定为随机变量与区间变量，是研究含不确定参数多体系统动力学常用方法。其中，将不确定参数认定为随机变量常采用蒙特卡洛方法或者多项式混沌展开逼近方法10。由于航天机构场地与成本的限制，无法进行大量试验来补充试验数据库，采用传统的不确定性分析模型对空间展开机构整体进行可靠性分析难以实现。为解决小样本航天机构中认知不确定分析问题，将航天机构中不确定参数进行处理。针对区间不确定量化11-12方法，对代理模型13（包括多项式模型）、Kringing模型以及神经网络14-17等进行量化分析，其中通常选取多项式模型为代理模型。WU 等18利用泰勒多项式展开与子区间技术相结合的方法与切比雪夫区间方法研究了含区间参数多体系统的响应分析，从而有效避免区间误差叠加导致长时间计算精度不准的问题，降低含区间参数求解中的“过度估计”。邓汉卿19针对空间展开索网天线对机构自由度进行分析，通过运动螺旋系提出广义机构的逆运动综合方法，并给出展开天线的测地线索网生成方法。肖剑20基于极限状态函数、有限元方法与反向传播（Back Propagation，BP）神经网络，提出针对空间展开绳系的高精度全局代理模型。XIANG 等21基于切比雪夫区间求解方法对考虑不确定间隙参数与旋转关节摩擦系数的多体系统进行了灵敏度分析，弥补了预测与测量之间的偏差，提高了机构响应参数变化的灵敏度。MA 等22建立了基于实验拟合参数的改进薄层单元（TLEs）方法来模拟界面上接头的动态特性，并基于切比雪夫多项式分析了含区间参数转子系统的时域响应。然而，上述分析方法在进行不确定性分析时，由于区间参数取不同样本点时，得到的动力学响应出现不同的频率与相位差，且随着时间的增大而逐步扩大，导致瞬时相位累积效应。由此，CUI 等24结合不确定分析方法与信号分解，将系统的时域响应分解为若干振动项与趋势项的叠加，并在每一时刻下建立相应代理模型，解决了传统不确定分析长时域信号瞬时相位累积效应。本文通过对不确定间隙进行描述，结合传统接触模型提出含不确定区间间隙接触模型，并将其嵌入系统动力学模型，通过切比雪夫扩张函数对含不确定区间系统动力学模型进行求解，提出一种相适应的局部均值分解方法、瞬时叠加方法与区间扩张方法进行区间叠加处理，揭示系统区间累加下动力学响应边界分布规律。1不确定间隙描述 在传统含间隙多体系统动力学分析过程中，普遍将间隙认定为旋转关节孔半径与轴半径的相减，推导过程基于初始最大间隙所得。在展开机构运转过程中，由于多间隙的耦合作用，针对旋转关节间隙处发生接触碰撞-碰撞分离-再次碰撞的运动过程中的碰撞分离过程，旋转关节产生碰撞分离过程并不是每次都可以到达初始间隙最大值再次发生下次碰撞，如图 1 所示。每一次发生接触碰撞的初始间隙是一个不确定量，其不确定长度可以表示为R=RxRyT（1）式中：Rx为旋转关节初始间隙不确定量横向投影距离；Ry为旋转关节初始间隙不确定量纵向投影距离。图 1旋转关节轴心任意时刻位置Fig.1Diagram of the position of the axis of the rotary joint at any moment89第 40 卷 2023 年第 4 期上海航天（中英文）AEROSPACE SHANGHAI（CHINESE&ENGLISH）将旋转关节间隙与相连接的连杆长度简化为开环机械臂两段，利用欧氏运动群上的卷积探究开环机械臂的运动传播问题，给出两端机械臂运动工作空间密度为1(g)与2(g)，g=(R，x)SE(n)，对这 2 个工作空间密度进行卷积，得到整体机械臂运动范围密度为 1，2(R，x)=A SO()ny Rn1(A，y)2(ATR，AT(x-y)dAdy（2）式 中：SO(n)为 机 械 臂 横 向 运 动 区 域 量，R 为 全域量。在旋转关节运转空间环境中，假设间隙角在运转角度0，2内均匀分布：=12，0，2（3）间隙角在 X方向累计分布函数可以表示为 P(cos X)=P(arccos X 2-arccos X)=P(2-arccos X)-P(arccos X)=1-arccos X，X-1，1 (4)同理，间隙角在 Y 方向累计分布函数可以表示为 P(sin Y)=1-P(arcsin Y 2-arcsin Y)=arcsin Y，Y-1，1 (5)对间隙角在 X 与 Y 方向上的累计分布函数分别求导，可得间隙角在 X 与 Y 方向上的概率密度函数x与y：x=P(cos X)=11-X2（6）y=P(sin Y)=11-Y2（7）利用欧氏运动群上的卷积公式对式（6）与式（7）间隙角概率密度函数进行卷积积分，可得间隙角末端点在间隙圆内分布概率密度函数为 P(X，Y)=-X-Yxyd cos d sin =12-11-cos sin 1-X21-Y2d（8）根据卷积在Rn的性质，可以得到1，2，n=i=1ni（9）1，2，n=i=1ni（10）结合式（9）与式（10），间隙角末端点分布概率密度函数均值与可表示为=-+xi(x)dx（11）i=-+(x-i)(x-i)Ti(x)dx=2iQn （12）式中：Qn为n n的单位矩阵。对于运动范围不受限制的旋转关节，其间隙角末端点分布可能为间隙圆内任意一点，根据中心极限定理与迭代关系，利用高斯概率函数构建间隙角末端点在 X方向与 Y方向分布概率密度函数：(x，x，x)=1()2i12exp (-12(x-)Ti-1(x-)（13）(y，y，y)=1()2i12exp (-12(y-)Ti-1(y-)（14）结合式（13）与式（14）间隙角末端点在 X方向与Y 方向分布概率密度函数，可得间隙角末端点全域联合分布概率密度函数为 f(x，y)=12xyexp(-12()x-x22x+()x-y22y-2()x-x()x-yxy)（15）在不考虑嵌入量的情况下，间隙角末端点分布最大值为间隙圆边界R，因此在考虑机构运动行程与分布区间大小，给出间隙角末端点在间隙圆概率分布密度函数：f(xi，yi)=12xyexp，-12()x-x22x+()x-y22y-2()x-x()x-yxy，xi2+yi2 R20，otherwise（16）90第 40 卷 2023 年第 4 期张渤洋，等：考虑区间叠加不确定间隙的空间可展机构动力学研究2含不确定机构动力学建模 在考虑实际接触过程中，随着空间展开机构多单元的叠加特点，导致不同位置处的间隙为一未知量。通过对不同位置处间隙采用区间变量的形式进行统一描述，对空间展开机构进行高精度建模分析。根据 Winkler 弹性基础模型25，接触体在接触位置产生一定变形，如图 2所示，接触区域径向位移与最大法向穿透深度的关系满足刚性接触体的几何边界条件，其变形量可表示为un()=R+sin(R2cos +R12-sin2R22)（17）式中：R1、R2分别为接触体半径；R为初始间隙；为半接触角度。轴孔接触模型半接触角度处的法向接触应力为p()=Eun()R1+c（18）式中：E为接触体弹性模量；c为接触区域弹性层厚度。在接触区域上任取一弹性微元，对式（18）在半接触区域内进行积分，r由于轴孔初始间隙R和轴孔嵌入量均远小于轴半径R1，接触点与固定坐标系原点之间的距离近似于接触体半径R1。则接触区域承受的接触力为FN=ER1R1+c0un()cos d（19）将式（17）代入式（19），忽略高阶项，接触域单位接触长度接触力为 FN=2R1ER1+c R+sin22R2+cos 2R21-sin2R22（20）结合图 1 旋转关节在发生三接触状态的过程中，其每一次碰撞点初始位置都不尽相同，导致接触体瞬时接触力产生差异。其考虑不确定区间间隙的接触力模型可进一步表示为 FN=2R1ER1+cmin R+sin22R2+cos 2R21-sin2R22（21）式中：R为不确定初始间隙值，结合式（16）可表示为R=R2-R1-|f(x)|。采用拉格朗日方程对空间展开机构进行动力学建模，空间展开机构展开单元动力学方程可表示为ddt(L i)-Li=Fi+Fc，i=1，2，n（22）式中：L为拉格朗日函数，定义为空间展开单元动能与势能的差值；i为广义坐标；Fi为广义坐标对应的广义力；Fc为接触体的接触力。3考虑不确定参数切比雪夫求解方法 切比雪夫区间分析方法作为应用最为广泛的非嵌入式多体系统动力学解算分析方法，采用切比雪夫函数进行区间求解，不仅在计算含区间不确定的多体系统动力学方程时有较高的计算精度，而且其扩张函数对考虑不确定间隙导致系统动力学方程维度增加时，区间计算产生的累加效应能够起到一定的抑制作用，可以用作较长时间分析含区间不确定的多体系统动力学响应特性。3.1切比雪夫扩张函数基于 Weierstrass 理论，在定义区间a，b上连续函数采用多项式的形式逼近，对于一维问题，在定义区间a，b上的 p阶切比雪夫级数为Cn(x)=cos(n)，=arccos 2x-()a+bb-a （23）式中：变量x a，b 且 0，。图 2接触体二维接触模型Fig.2Two-dimensional contact model for contact bodies91第 40 卷 2023 年第 4 期上海航天（中英文）AEROSPACE SHANGHAI（CHINESE&ENGLISH）对于一维单变量连续函数f(x)a，b 可用切比雪夫扩张函数近似表示为f(x)12f0+i=1pfiCi()x（24）式中：fi为多项式的系数。将每一个一维多项式表示为张量积，从而给出n维切比雪夫多项式定义为Ci1，i2，in(x1，x2，xn)=cos(i11)cos(i22)cos(inn)，in=0，1，n （25）式中：in为多项式的阶数；n=arccos 2xn-()an+bnbn-an 多元函数f(x1，x2，xn)用 n次切比雪夫多项式近似可表示为f(x1，x2，xn)=j1=0njk=0n()12pfj1，j2，jkcos(j11)cos(j11)cos(jkk)（26）式中：fj1，j2，jk为 n次多元函数的系数，基于 Mehler积分可知 j1，j2，jk=(2)k00f()x1，x2，xkcos(j11)cos(j22)cos(jkk)d1d2dk(2)kj1=1mjk=nmf()xj11，xj22，xjkkcos(j1j11)cos(j2j22)cos(jkjkk)(27)式中：m 为多项式插值点数，一般取m=p+1可保证近似精度。jkk=()2jk-12m，jk=1，2，m（28）xjkk=ak+bk2+bk-ak2cos jkk，jk=1，2，m （29）综上所述，含区间不确定性的空间展开机构动力学方程，可通过切比雪夫多项式的形式，构建响应函数的代理模型，并采用区间计算的方法计算系统动力学方程的区间解。3.2切比雪夫扩张函数求解方法本节针对多体系统动力学模型，提出切比雪夫扩张函数对动力学模型进行求解，通过切比雪夫正交多项式构建代理模型，完成对动力学模型求解，切比雪夫求解流程如图 3所示。在图 3中，n为含区间间隙空间展开单元动力学方程维数，y0为动力学响应初始值，t0为初始时间，tend为结束时间；h为切比雪夫区间求解的迭代步长；k为切比雪夫多项式阶数。采用切比雪夫多项式构建含区间间隙空间展开单元动力学方程代理模型时，将具有含不确定区间参数的空间展开单元动力学方程，转换为仅具有特定区间参数的动力学方程并进行求解，根据多个所求数值解构造切比雪夫代理模型，基于区间算法来计算含不确定区间参数的空间，展开机构动力学方程的区间解。其区间解可用如下形式表示：yj+1(cos j1，cos j2，cos jn)=y=f()t，y，i，j1，j2，jny()tj=yj()cos j1，cos j2，cos jn（30）式中：yj+1为tj+1时刻含不确定区间参数的空间展开单元动力学方程数值解；y(tj)为tj时刻插值点处含不确定区间参数的空间展开单元动力学方程数值解。4小区间叠加方法 针对区间求解方法中区间叠加误差过大问题、分解的小区间整合成大区间的问题与启动阶段估计不准问题，提出局部均值分解方法、瞬时叠加方法与区间扩张方法，为提高切比雪夫求解精度大区间分解的小区间整合成大区间下的区间边界，并再此基础上降低启动阶段估计精度。通过小区间累加方法在切比雪夫求解方法上进行进一步处理，获得较高精度区间边界值。图 3切比雪夫求解动力学方程流程Fig.3Flow chart of Chebyshev dynamic equation solving92第 40 卷 2023 年第 4 期张渤洋，等：考虑区间叠加不确定间隙的空间可展机构动力学研究4.1局部均值分解方法针对切比雪夫区间求解方法区间叠加导致误差累积效应，在长时域内求解仅保持短时间有效的缺点，对切比雪夫求解获得动力学响应边界x(，t)采用局部均值分解的方法，将响应边界分解为若干个响应幅值与响应相位乘积之和的形式，对动力学响应边界x(，t)中响应幅值a(，t)与响应相位(，t)建立代理模型，解决传切比雪夫区间求解方法相位累计效应问题。其中响应幅值通常是缓变甚至单调的信号，响应相位指幅值为 1的函数，代理模型指计算结果与原模型非常接近，并便于求解的模型。将所求动力学响应边界x(，t)在每一时刻下建立瞬时响应幅值模型与响应相位模型，完成对初始动力学响应边界的分解。具体步骤如下：取初始残差函数ri(t)等于初始动力学响应边界x(，t)，提取初始动力学响应边界中的极大值点ai与极小值点bi，对提取的所有极大值点ai进行局部均值拟合，得到上包络线vi(t)；同理，得到下包络线ui(t)，将上包络线vi(t)与下包络线ui(t)进行瞬时均值处理，得到动力学响应边界瞬时均值ci(t)为ci(t)=ui()t+vi()t2（31）将响应边界x(，t)除以瞬时均值ci(t)的绝对值并赋值给瞬时响应相位gi(t)：gi(t)=x()w，t|ci()t|，i=1，2，n-1（32）动力学响应瞬时幅值函数在时域内累加可表示为Ci(t)=ci1(t)ci2(t)cin-1(t)=n=1n-1cin()t（33）动力学响应瞬时相位函数在时域内累加可表示为i(t)=arccosgi()t|gi()tmax|（34）i(t)=i1(t)i2(t)in-1(t)=n=1n-1in()t （35）对初始信号全部采样点处理完成后，i(t)在时域内累加构建为动力学区间响应相位i(t)，ci(t)时域内累加构建为动力学区间响应幅值Ci(t)，将动力学区间响应幅值函数与动力学区间响应相位函数相乘，可得第 i 阶乘积函数动力学区间响应代理模型PFi(t)：PFi(t)=Gi(t)|Ci(t)|（36）新的残差函数ri(t)可表示为ri(t)=x(，t)-PFi(t)（37）综上所述，通过局部均值分解方法，将动力学区间响应边界分解为动力学响应幅值函数与动力学响应相位函数乘积之和的形式，具体表达形式为x(，t)=i=1nCi(t)cos i(t)+ri(t)（38）4.2瞬时叠加方法基于上节提出的局部均值分解方法，结合传统响应面分析方法，提出多体系统动力学瞬时叠加方法，对切比雪夫区间求解过程中大区间分解后的小区间进行叠加，获得初始区间下系统动力学响应边界，具体可分为以下步骤。1）输入动力学响应边界、不确定间隙r，在不确定间隙r中选取 n 个样本，记作：r1，r2，rn，n为样本序号，n 个不确定参数对各样本对应的时域内响应分别为x(r1，t)，x(r2，t)，x(rn，t)。2）采用 4.1 节提到的局部均值分解方法，将时域内的不确定参数样本响应分解为 n个动力学响应幅值函数与动力学响应相位函数乘积之和与相对应的残差函数。其中分解结果如下：x(ri，t)=j=1nCi(ri，t)cos i()ri，t+rj(ri，t)（39）3）在每一时刻下，基于多项式函数拟合建立切比雪夫区间求解响应的瞬时响应幅值函数与瞬时响应相位函数，以及残差函数的代理模型，得到相应代理模型i(t)、ci(t)、ri(t)，分别逼近真实解。4）计算每一时刻下，动力学响应边界分解的瞬时响应幅值函数代理模型与瞬时响应相位函数代理模型的耦合形式ci(t)cos i(t)。5）在耦合形式代理模型中加入残差函数代理模型，计算完成形式的代理模型，得到完全的谐波叠加形式的代理模型为givent：x(a，t)=i=1nci(ai，t)cos i()ai，t+ri(ai，t)（40）6）计 算 代 理 模 型 上 边 界 与 下 边 界。如 图 4所示。93第 40 卷 2023 年第 4 期上海航天（中英文）AEROSPACE SHANGHAI（CHINESE&ENGLISH）4.3区间延拓方法在切比雪夫区间求解过程中，在启动阶段出现区间求解误差较大问题，启动阶段前没有相对应的区间与启动阶段进行关联，单一区间求解状态下容易出现误差较大的情况。启动阶段以前，区间通过余弦拓展方法进行扩张，在求解过程中启动阶段有2 个甚至多个区间联合求解，解决启动阶段区间误差较大问题。通过切比雪夫区间求解获得系统动力学响应边界以后，可以根据零点附近动力学响应变化趋势，在动力学响应零点以前区间增加相应的相位、幅值与频率适当的余弦延拓函数，使动力学响应区间代理模型可以在零点处有更多区间与相应响应边界函数，从而提高启动阶段所得上下包络线精度，余弦延拓函数及参数定义如下：xi(t)=ci(t)cos(2tTi+i(t)+ri（41）式中：ci(t)为代理模型瞬时幅值函数；i(t)为代理模型瞬时相位函数；ri为代理模型误差；Ti为延拓周期数，在数值上等于延拓周期与相邻极大值与极小值之前时间差的 n倍。Ti=n|x(ji)-x(ki)|（42）通过上述余弦延拓方法，可保证代理模型在展开过程任意有效时域内都有相邻区间关联求解，从而降低启动过程中单一区间求解带来的动力学响应区间边界误差。5系统区间累加动力学响应边界分析 以 平 板 式 合 成 孔 径 雷 达（Synthetic Aperture Radar，SAR）天线展开机构为研究对象，如图 5 所示，通过区间求解来分析区间间隙对空间展开天线展开过程的动力学特性影响，选取区间间隙大小为0.01，0.1 mm，连杆刚度系数为 5105 N/m，动摩擦系数为 0.1。多体系统动力学特性边界计算流程如图6所示。对多体系统旋转关节间隙进行不确定参数量化取值，将所得区间值代入区间接触模型进行计算并嵌入到含不确定间隙多体系统动力学模型，通过切比雪夫扩张函数进行求解，并对所得结果通过小区间叠加的区间处理方法进行处理，得到动力学特性边界。0.06，0.1 mm 区间间隙标记点小区间叠加区间结果与切比雪夫区间计算对比如图 7所示。对区间计算结果进行区间划分，并统计每个区间下概率密度，同时由区间计算给出该间隙下标记点误差分布区间。采用小区间叠加方法可明显降低区间边界处由代理模型长时域区间迭代产生的误差累计效应。将区间叠加计算平均结果与传统拟合结果进行对比如图 8 所示，由图 8 可知在 0.06，0.1 mm间隙下区间叠加计算平均结果与传统拟合结果相图 4瞬时叠加方法流程Fig.4Flow chart of the instantaneous superposition method图 5空间展开机构仿真Fig.5Simulation schematic of the spatial unfolding mechanism图 6多体系统动力学特性边界计算流程Fig.6Flow chart of the dynamics boundary calculation for multi-body system94第 40 卷 2023 年第 4 期张渤洋，等：考虑区间叠加不确定间隙的空间可展机构动力学研究比提高 7.74%，0.03，0.06 mm 区间叠加计算平均结果与传统拟合结果相比提高 3.73%，可进一步提高空间可展机构展开预示精度。采用局部瞬时分解方法、瞬时叠加方法与区间扩张方法对不同区间间隙大小与数量下区间边界进行处理，得到大区间下位移偏差区间分布状况，如图 9和图 10所示。由图 9可知，随着间隙区间的增大，系统位移偏差基准值出现逐渐增大的趋势，但是区间范围扩张并不明显。由图 10可知，小区间叠加方法可以有效解决小区间与大区间之间的转换、区间累加效应与边界估计不准的问题。随着生效区间间隙数量的增多，系统位移偏差基准值趋于稳定，但是区间范围产生明显增大。说明间隙大小对系统运动精度产生较大影响，而间隙数量对系统运动稳定度影响较为严重，在对空间展开机构展开过程分析中间隙对系统的影响时，不仅要合理控制各旋转关节间隙大小，同时还要关注由于多级机构叠加从而导致系统间隙数量倍增的状况。6结束语针对空间展开机构动力学响应预示不准问题，基于不确定区间间隙改进传统共形接触模型，并构建含区间不确定空间可展机构动力学模型，提出一种相适应的切比雪夫区间求解方法。得到主要结论如下：1）针对区间求解过程中区间叠加误差过大、区间整合困难与启动阶段估计不准等不足，分别提出局 部 均 值 分 解 方 法、瞬 时 叠 加 方 法 与 区 间 扩 张方法。2）区间间隙大小主要对系统运动精度产生较大影响，而区间间隙数量对系统运动稳定度影响较为严重。3）区间叠加方法相比于传统拟合方法可进一步提高空间可展机构预示精度，相较于传统拟合结果可提高 7.74%。图 7 0.06，0.1 mm 区间间隙标记点区间结果与区间计算结果对比Fig.7Comparison of the interval results of the marked gap points with interval calculation results in the 0.06，0.1 mm interval图 8不同区间间隙大小下标记点区间叠加与传统拟合对比Fig.8Comparison of the interval superposition of marked points with conventional fitting under different interval gap sizes图 9不同间隙大小下叠加区间边界Fig.9Boundaries of the stacking intervals with different gap sizes图 10不同区间间隙数量下叠加区间边界Fig.10Boundaries of the stacking intervals with different gap numbers95第 40 卷 2023 年第 4 期上海航天（中英文）AEROSPACE SHANGHAI（CHINESE&ENGLISH）参考文献1 LI Y Y，WEI J F，DAI L.Structural design and dynamic analysis of new ultra-large planar deployable antennas in space with locking systemsJ.Aerospace Science and Technology，2020，106：106082.2 GUO J W，ZHAO Y S，XU Y D，et al.A novel modular deployable mechanism for the truss antenna：assembly principle and performance analysis J.Aerospace Science and Technology，2020，105：105976.3 秦远田，张子昊，朴丽花.构架式天线柔性动力学建模与展开精度分析 J.上海航天（中文），2023，40（2）：13-22，31.4 庞梦非，朱春艳，张美艳，等.具有不确定性连接刚度的太 阳 电 池 阵 基 频 分 析J.上 海 航 天，2017，34（6）：103-108.5 BAI Z F，ZHAO J J，CHEN J，et al.Design optimization of dual-axis driving mechanism for satellite antenna with two planar revolute gap jointsJ.Acta Astronautica，2018，144：80-89.6 YANG L X，ZHANG X M，HUANG Y J.Dynamic analysis of open-loop mechanisms with multiple spatial revolute gap joints J.Proceedings of the Institution of Mechanical Engineers Part C-Journal of Mechanical Engineering Science，2019，233（2）：593-610.7 冯伟，刘保国，丁浩，等.非参数不确定动力学建模研究综述 J.振动与冲击，2020，39（5）：1-9.8 余德文，赵强强，陈飞飞，等.利用结构等效转换的空间展开机构装配误差建模与灵敏度分析 J.西安交通大学学报，2019，53（4）：15-23.9 张海瑞，王浩，王尧，等.飞行器总体不确定性建模与优化设计方法 J.宇航学报，2023，44（4）：486-495.10 杨旭锋.含认知不确定性参数的高效高精度可靠性分析方法研究 D.西安：西北工业大学，2016.11 SCHWER L E.An overview of the PTC 60/V&V 10：guide for verification and validation in computational solid mechanics J.Engineering with computers，2007，23（4）：33-40.12 HYUNG B L，URMILA G，SAMI B，et al.Development and use of engineering standards for computational fluid dynamics for complex aerospace systemsC/46th AIAA Fluid Dynamics Conference，WDC，Reston，Virginia：AIAA，2016：56-63.13 王智明，郑雷，张露熹，等.基于混沌多项式展开的机械结构可靠性分析 J.机械强度，2022，44（4）：852-858.14 LI J，CHEN J，WEI J，et al.Developing an instantaneous response surface method tIRS for time dependent reliability analysis J.Acta Mechanica Solida Sinica，2019，32（4）：446462.15 SU W H，CHOU C S，XIU D.Deep learning of biological models from data：applications to ode modelsJ.Bulletin of Mathematical Biology，2021，83（3）：119.16 HAN S，CHOI H S，CHOI J，et al.A DNNbased datadriven modeling employing coarse sample data for realtime flexible multibody dynamics simulationsJ.Computer Methods in Applied Mechanics and Engineering，2021，373：113480.17 NI B Y，JIANG C，HUANG Z L.Discussions on non-probabilistic convex modelling for uncertain problemsJ.Applied Mathematical Modelling，2018，59：307904.18 YASEEN Z M，ALDLEMY M S，SADEGH M O.Non-gradient probabilistic Gaussian global-best harmony search optimization for first-order reliability methodJ.Engineering with Computers，2019（4）：1-12.19 WU J，LUO Z，ZHANG Y，et al.Interval uncertain method for multibody mechanical systems using Chebyshev inclusion functions J.International Journal for Numerical Methods in Engineeri