An improved numerical method for constructing Halo/Lissajous orbits in a full solar system model

An improved numerical method that can construct Halo/Lissajous orbits in the vicinity of collinear libration points in a full solar system model is investigated. A full solar system gravitational model in the geocentric rotating coordinate system with a clear presentation of the angular velocity relative to the inertial coordinate system is proposed. An alternative way to determine patch points in the multiple shooting method is provided based on a dynamical analysis with Poincare′sections. By employing the new patch points and sequential quadratic programming, Halo orbits for L_1, L_2, and L_3 points as well as Lissajous orbits for L_1 and L_2 points in the EarthMoon system are generated with the proposed full solar system gravitational model to verify the effectiveness of the proposed method.     

基于正交试验的航行体弹道仿真优化设计

为确保水下发射航行体再入水后弹道安全,分析了航行体弹道影响因素,构建了带有阻尼板航行体的弹道数学模型,将发射后航行体与发射载体的最近距离作为安全性判定标准,采用正交试验方法,对水下航行体的阻尼板数量、尺寸及张开角度3个参数进行了仿真优化分析。仿真结果表明:当阻尼板数量为6块、长度为800 mm、张开角度为80°时,水下发射航行体质心z方向侧移量最大,航行体尾端与发射载体假想壁距离最大,安全性最高。     

全模颤振风洞试验三索悬挂系统多体动力学分析

全模颤振风洞试验需要通过软支撑系统模拟飞行器的自由飞行状态并调整模型姿态达到配平状态。参考NASA双索悬挂方案,提出了一种两电机驱动的三索悬挂系统,利用后方两索的同向/反向联动实现模型俯仰和滚转姿态的调整,利用弹簧刚度以及钢绳张力设计实现支撑频率要求。基于柔性多体动力学方法,建立了包括飞行器刚体模型、柔性索、滑轮、弹簧、气动力模型、伺服电机控制在内的复杂系统动力学模型,其中,利用任意拉格朗日-欧拉（ALE）变长度索单元描述钢绳,利用不约束物质坐标的索结点约束描述钢绳与滑轮相互作用,利用索结点物质输运速度约束描述伺服电机绞盘,利用飞行力学的气动力模型描述吹风下的气动力。基于该模型,通过小扰动响应辨识研究了弹簧刚度、钢绳张力、连接点位置等因素对支撑频率的影响规律,并分析了系统姿态调整能力,俯仰调整范围达到-12.5°~12.5°,滚转调整范围达到-45°~45°。采用滑轮处电位计测量的钢绳相对位移作为反馈信号,基于设计的控制律利用多体动力学求解器与Simulink对风洞吹风下的姿态调整过程进行仿真,模型达到配平状态,获得了吹风下的索拉力和伺服电机功率,为系统设计提供基础。     

一种针对局部结构的有限元模型修正方法

针对实际结构有限元模型（FEM）的建模误差通常仅存在于局部区域,提出了一种对局部结构单独进行模型修正的方法。首先,根据频响函数（FRF）解耦理论得到由残余结构频响函数与包含待修正参数的局部结构动刚度所重构的整体结构频响函数的拟合值,然后通过迭代优化使其与测量值的残差最小化,从而得到参数的极大似然估计。在此基础上,将残差关于参数的灵敏度以局部结构动力学矩阵表示,建立了模型修正的基本方程,利用整体结构的测试数据即可直接对分离出来的局部结构进行模型修正。最后,对喷气式飞机和三角机翼飞机分别进行了数值模拟和实验研究,验证了所提方法的可行性和有效性。结果表明,所提方法可以成功地用于仅局部区域含有建模误差的实际结构有限元模型的修正,修正后的有限元模型的动态特性与实际结构有较好的一致性。     

高超声速滑翔目标自适应跟踪方法

针对机动模式复杂多变的高超声速滑翔目标跟踪问题,提出了一种机动频率自适应跟踪方法。采用介于常速度和常加速度模型之间的Singer模型来表征目标气动力加速度的变化,从而建立跟踪系统的状态方程。根据地基雷达量测量获得系统的量测方程,鉴于距离和角度信息的量级相差较大将其由球形量测量转换为位置量测量。为了适应高超声速滑翔目标灵活多样的机动模式,基于正交性原理和无迹卡尔曼滤波算法实现了Singer模型中机动频率参数的自适应。利用滤波信息计算得到能够反映状态模型误差大小的调整因子,用于放大Singer模型中的机动频率,进而调整状态方程的过程噪声以降低模型误差。通过对2种典型机动轨迹的跟踪仿真,并与交互式多模型等方法进行比较,结果表明所提方法的跟踪精度高、计算量小,能够较好地适应阶跃机动和连续幅值变化的机动。     

基于自适应Kriging代理模型的交叉熵重要抽样法

对于复杂失效域和小失效概率耦合的可靠性分析问题,本文提出了一种交叉熵重要抽样（CE-IS）方法结合自适应Kriging （AK）代理模型的求解方法（CE-IS-AK）。所提方法基于交叉熵原理,用混合高斯模型逐步逼近最优重要抽样密度函数,并采用AK模型协助逼近过程中混合高斯模型的参数的更新,从而提高了CE-IS方法的计算效率。另外,本文还改进了CE-IS方法的收敛准则,避免了方法的冗余迭代,扩大了方法的适用范围。由于在CE-IS方法中引入了AK模型,因此,本文方法所构建的重要抽样函数在保证精度的基础上提高了效率。相较于AK-MCS方法,本文方法中引入了重要抽样的思想,因此在Kriging训练点数目基本相同的情况下,大幅缩减小失效概率计算时样本池规模,并且由于利用了混合高斯模型,因而对多失效域具有较好的适用性。算例分析也证明了本文所提方法的优越性。     

基于随机粒子追踪方法的环形燃烧室点火动态模拟

为了深入研究低排放燃烧室点火联焰规律,在全新的环形模型燃烧室中开展了点火模拟和试验研究。点火模拟采用随机粒子追踪方法,能够基于时均冷态流场的仿真结果快速模拟火焰传播过程。环形燃烧室包含16个中心分级旋流器,仅向预燃级通入丙烷,用于模拟航空发动机低排放燃烧室点火状态下的空气燃油分级。试验采用PIV技术测量3个头部区域流场,利用高速相机拍摄火焰CH*/C2*基团化学发光信号。对多个流量和当量比条件下的联焰过程、联焰时间和传焰速率进行了分析,试验和模拟的结果均表明：环形燃烧室内火焰双向传播,燃烧室内外环流速度差异导致了双向火焰传播速度差,传焰速率随燃烧室湍流速度和当量比的增加而增加。点火模型很好地捕捉了环形燃烧室点火动态,所得传焰速率也符合湍流火焰传播规律,表明该模型具有较强预测能力。     

多重应力下滚动轴承剩余寿命预测

为了实现多重应力下滚动轴承的剩余寿命预测,有效利用不用应力下的退化数据,提出了一种基于加速模型和贝叶斯（Bayesian）理论的滚动轴承剩余寿命预测方法。通过拟合优度检验和威布尔（Weibull）概率图检验法对滚动轴承试验中的数据进行有效性分析。利用switching Kalman filters(SKF)判断滚动轴承各时刻的退化状态。当滚动轴承进入加速退化时,用指数模型拟合轴承退化过程,利用广义线性对数模型表示退化模型参数与应力的关系,根据修正后的轴承实时退化数据利用贝叶斯算法更新模型参数,得到滚动轴承剩余寿命的概率密度函数,从而实现滚动轴承剩余寿命预测。采用XJTU-SY轴承数据集进行验证,预测结果的均方根误差在20 min以内,证明该方法能够有效预测滚动轴承的剩余寿命。     

民机襟翼交联机构吸能元件的吸能特性

采用有限元分析与冲击试验相结合的方式对膨胀管与波纹管的吸能特性进行探究：通过数值仿真分析对吸能元件的试验过程进行初步掌握,同时根据冲击试验结果对试验设计的不合理之处进行修正;然后开展两种吸能元件的动态冲击试验,探究其吸能特性,最终比较分析有限元仿真分析与落锤冲击试验的试验结果,选择效果理想的吸能元件应用到襟翼交联机构上。结果表明：无论是从吸能行程、吸能稳定性还是吸能结构的设计考虑,波纹管都是最适合应用到交联机构中的吸能元件,该吸能元件比吸能大,吸能过程平缓稳定,吸收能量之后只有塑性形变但结构不发生破坏,同时吸能行程短,吸能结构简单,便于安装拆卸。     

Stochastic modelling considering ionospheric scintillation effects on GNSS relative and point positioning

Global Navigation Satellite Systems (GNSS), in particular the Global Positioning System (GPS), have been widely used for high accuracy geodetic positioning. The Least Squares functional models related to the GNSS observables have been more extensively studied than the corresponding stochastic models, given that the development of the latter is significantly more complex. As a result, a simplified stochastic model is often used in GNSS positioning, which assumes that all the GNSS observables are statistically independent and of the same quality, i.e. a similar variance is assigned indiscriminately to all of the measurements. However, the definition of the stochastic model may be approached from a more detailed perspective, considering specific effects affecting each observable individually, as for example the effects of ionospheric scintillation. These effects relate to phase and amplitude fluctuations in the satellites signals that occur due to diffraction on electron density irregularities in the ionosphere and are particularly relevant at equatorial and high latitude regions, especially during periods of high solar activity. As a consequence, degraded measurement quality and poorer positioning accuracy may result.