月球车地球敏感器图像地心位置提取算法

地球是月球上可视半径最大的天体,且在天球上运动范围较小,不存在升降现象.当月球车在月球对地面区域活动时,利用地球敏感器对地球成像可实现月球车长期自主天文导航.地球图像地心位置提取是利用地球敏感器进行天文导航的关键技术之一,直接决定了地球敏感器的观测精度.本文通过研究地球敏感器镜头投影模型分析地球成像规律,提出一种不受地相变化约束的地球敏感器图像地心位置提取算法,采用取半搜索法和循环搜索法两步实现地球真实边缘线的充分筛选,并基于此拟合地心位置.半物理仿真实验校验结果表明,本算法能针对不同地相图像有效提取地球中心,外符合平均精度约为9.78"~16.68",在实验条件随机改变的情况下,地心位置外符合精度标准差互差最大不超过0.98".      

红外地球敏感器测量值修正算法及其应用研究

 圆锥扫描式红外地球敏感器常用于测算地心矢量,而地球扁率影响地心矢量精度。研究了基于地球扁率的地心矢量的修正算法,并将该方法应用于双圆锥红外地球敏感器,给出了考虑扁率的地心矢量计算方法。仿真结果表明,该方法能有效地提高地心矢量的计算精度。     

空间交会V—bar接近冲量机动运动分析

从相对运动与绝对运动两方面全面阐述空间交会V—bar接近冲量机动变轨的运动规律。交会接近段的起点与终点为绝对运动转移轨道通径的两个端点，转移轨道半长轴略大于机动变轨前的轨道半径，转移轨道偏心率以及追踪星对目标星缩减的地心角取决于径向速度增量，并给出视界角设计方法。     

地球卫星运动中坐标系附加摄动与参考系选择问题

针对地球赤道面的摆动,曾提出过两种空间地心赤道坐标系,即历元地心平赤道坐标系和混合形式的轨道坐标系.在后一种坐标系中虽然可避免计算因地球引力位变化所引起的坐标系附加摄动的麻烦,但这又带来了两种坐标系之间的不统一和转换的麻烦.针对实际现状,直接在历元地心平赤道坐标系中给出相应的坐标系附加摄动解,所有与此有关的转换问题,均可在历元地心平赤道坐标系中进行,而无需为了避免坐标系附加摄动的计算去通过轨道坐标系来实现.从而建议,不必再引进轨道坐标系,在卫星定轨和预报中,所有的转换问题均可采用同一坐标系,即历元(目前是J2000.0)地心平赤道坐标系.      

精确动力学模型中的行星引力辅助轨道设计

轨道设计分为初步设计和精确模型迭代两步,初步设计基于等效脉冲模型,用圆锥曲线拼接法确定时间窗口和引力辅助产生的速度脉冲。精确模型中引力辅助看作是一个连续的过程,将简化模型得到的引力辅助双曲线轨道化为行星心目标B平面参数,以地心逃逸速度作为设计变量,通过微分修正的方法进行求解。通过算例对比分析了简化模型和精确模型设计结果之间的关系,结果表明,引力辅助脉冲等效模型精度较好。     

地磁场模型简化计算方法

为了实现地磁导航在近地轨道卫星上的应用,需要解决传统的高斯球谐函数递推法在计算地磁场强度时存在方法复杂、计算量大、实时性难以保证等问题.针对上述问题,提出一种用地球磁场估计函数代替球谐函数的地磁场强度计算方法.该方法利用伪地心下的磁偶极子模型代替主磁场模型,通过多项式拟合离线得到高度范围内固定经纬度网格点的伪中心距系数,进而利用插值法和偶极子模型得到任意点的地磁场强度.仿真结果表明,轨道高度为300～500km,经纬度网格间隔为0.5°时,地球磁场估计函数法的导航精度与地球主磁场模型WMM精度一致的同时,降低了计算负担,提高了算法计算效率.      

蒙皮加筋圆柱壳扭转频率的三种计算模型

为了给火箭等加筋圆柱壳结构的简化计算提供理论基础,研究了加筋圆柱壳面积等效方法中的扭转刚度计算方法,说明了不能直接根据等效厚度计算扭转刚度的物理机理,得到了桁条对圆柱壳扭转刚度的贡献可以忽略的结论.根据等效厚度以及桁条与蒙皮连接关系的不同处理方法,给出了计算扭转固有振动频率的等效梁模型、变厚度壳模型和壳-梁模型,分析了桁条在3种模型中对刚度和质量(包括转动惯量)的不同作用.为了消除根据等效厚度计算扭转刚度和扭转固有振动频率的偏差,分别给出了扭转刚度和扭转固有振动频率的修正方法.数值结果验证了修正方法和所得结论的正确性.      

印刷线路板设计参数对镀通孔可靠性的影响

基于镀通孔(PTH, Plated Through-Hole)应力分布模型,定量地研究了PTH线路板厚度与PTH直径之比(高径比)、PTH半径与镀层厚度之比、线路板有效作用半径与PTH半径之比等几何设计参数以及线路板玻璃化温度对PTH镀层应变及寿命的影响.通过引入等效温度载荷,给出了不同玻璃化温度下的应力分布模型.研究表明,在屈服温度载荷附近时,PTH应变和寿命对几何尺寸非常敏感.在固定线路板有效作用半径与PTH半径之比后,PTH高径比影响较为明显.对于一定厚度的线路板,PTH半径越小,镀层最大应变越大,寿命越低.在玻璃化温度以上,应变急剧增加.研究还表明应变集中系数并非固定不变,而是在弹性和塑性范围内有所变化.      

蒙皮加筋圆柱壳弯曲频率的三种计算模型

研究了加筋圆柱壳弯曲刚度的面积等效和惯性矩等效方法,分别分析了面积等效时截面惯性矩的精度和惯性矩等效时面积的精度.给出了计算蒙皮加筋圆柱壳弯曲固有振动频率和模态的等效梁模型、壳-梁模型和变厚度模型,这三种模型虽然皆可以用于火箭动态特性分析,但不同模型中的桁条作用机理不同.给出了与三维壳-梁模型和变厚度模型对应的一维解析解,据此解释了模型的物理机理.给出了薄壁壳弯曲频率计算中的临界厚度概念,指出了等效厚度方法的不足.     

三维四向编织复合材料弹性性能的有限元预报

在四步法方型三维编织复合材料细观结构单胞模型的基础上,假设纱线具有六边形横截面,应用ANSYS软件,建立了实体有限元模型,该模型表面与预制件表面平行,更加符合三维编织复合材料的实际结构,且有利于力学性能的分析,同时,该模型中还考虑了织物的编织角与内部编织角的几何关系.基于该模型,计算了材料的弹性常数,分析了编织角和纤维体积含量对弹性常数的影响规律.结果表明,数值计算与实验结果吻合较好.此外,还确定了材料内部的应力场分布,为材料的强度预测奠定了基础.      

仅地心矢量测量的卫星转对地定向姿态控制方法

摘要： 针对角速度过大且无其陀螺测量的姿态异常情况,开展仅依赖地心矢量的对地姿态快速恢复控制方法研究.基于地心矢量测量信息给出了天平角及天平角矢量定义,提出一种综合天平角抑制与能量耗散的对地姿态恢复控制方法;针对具有星体角动量偏置的闭环控制非线性系统,采用几何方法求解获取各平衡点,并由其特性分析表明通过参数合理选取可使得期望稳定对地姿态为系统唯一稳定平衡点;利用数学仿真方法验证了所提出方法的有效性.     

基于代理模型的制导火箭炮发射诸元计算方法

针对制导火箭炮发射诸元的快速计算问题,提出了一种结合大样本数据和代理模型计算发射诸元的新方法。运用代理模型建立射角、无控弹道侧偏与炮位纬度、炮位高程、射向、射程、目标点高程及药温之间的函数关系,并根据射程和无控弹道侧偏的预测值对射向进行修正。仿真结果表明,高阶多项式响应面、相关函数为高斯函数的Kriging、高阶单项式径向基函数、核函数为高斯函数的最小二乘支持向量机、激活函数为正弦函数的超限学习机以及由上述单一代理模型构建的组合代理模型均具有较高的预测精度,各种单一代理模型对射角和无控弹道侧偏的预测时间均小于1 ms,证明了基于代理模型的射角和无控弹道侧偏预测方法切实可行,且通过对射向进行修正有效减小了由于地球自转引起的无控弹道侧偏。      

基于紫外敏感器和星敏感器的卫星自主导航

卫星利用紫外敏感器和星敏感器进行自主导航的方法中地心矢量的测量精度是影响导航精度的重要因素之一,而地心矢量的测量又受到地球扁率的影响。在考虑地球扁率的前提下,研究了地球扁率对地心矢量测量的影响,给出了基于卫星姿态的地心矢量的补偿方法。仿真结果表明该补偿方法具有较高的补偿精度,并且能有效地提高卫星自主导航精度。     

电离层薄层高度对电离层模型化的影响

利用IRI2012模型分析了电离层薄层高度的时空变化规律,提出了基于应用中STEC的电离层改正误差分析理论,分析了电离层薄层高度变化的相关影响.结果表明,电离层薄层高度变化对电离层穿刺点位置、投影映射函数值、电离层建模结果、电离层模型精化和电离层模型精度评估结果的影响较大.高度截止角为10°时,电离层薄层高度变化导致电离层穿刺点的经纬度差异最大可达3.2°,投影映射函数最高可引入约15.46%的误差,电离层建模结果差异和建模实用误差最高分别达9.71%,3.64%,采用不同薄层高度数据的电离层模型参数拟合和模型精化结果最大可引入约9.26%的误差,采用不同电离层薄层高度数据进行模型精度评定时最大可引入约9.62%的误差.根据这些研究结果可知:在实际应用中应采用电离层薄层高度模型,并选取较大的卫星高度截止角来减小薄层高度变化引入的误差;采用固定高度时,区域电离层建模采用与实际电离层薄层一致的固定高度;进行精度评估时,参考数据的电离层薄层高度与需要精度评估的电离层模型薄层高度相等.      

月地转移轨道快速设计与特性分析

对采用直接大气再入方式的月地转移轨道,考虑大气再入界面参数和地面落点位置约束,提出了一种基于双二体模型的快速设计方法。该方法分为内外两层迭代循环,内层循环使月心段轨道和地心段轨道在月球影响球边界处连续,并采用Lambert问题与Newton-Raphson法相集合的方法求解满足再入角约束的地心段轨道参数;外层循环通过调整地心段轨道倾角和轨道置入时间使月地转移轨道满足地面落点位置约束。分析表明,存在四种类型的月地转移轨道满足大气再入界面约束,分别为降 降型、降 升型、升 降型和升 升型。在此基础上,对四种类型月地转移轨道的近地点地心距、置入分布点、再入点分布等特性进行了分析。仿真结果验证了所提出方法的有效性。     

光电二极管的地球反照光校正及卫星姿态估计

光电二极管作为体积小、成本低的太阳敏感器,可以结合地球敏感器进行卫星三轴姿态估计,但是地球反照光对其具有严重的影响,从而限制了其应用,然而地球反照光数学模型应用复杂。针对此问题,首先建立了一种简化的光电二极管量测模型,将地球反照光设成动态偏置项补偿在光电二极管模型中,并将误差建模为混合高斯噪声。然后应用滑窗估计和随机权重策略动态估计和更新模型中的参数。同时,采用多比例因子分别估计各光电二极管的地球反照光影响,并引入Huber影响函数处理异常值,从而提高了算法鲁棒性和参数估计精度。实验结果表明,采用地球反照光校正的光电二极管量测模型和无迹卡尔曼滤波（UKF）算法可实现卫星的高精度姿态估计,三轴姿态精度能达到0.2°~0.3°。      

The atmospheric radiation response to solar-particle-events.

High-energy solar particles, produced in association with solar flares and coronal mass ejections, occasionally bombard the earth's atmosphere. resulting in radiation intensities additional to the background cosmic radiation. Access of these particles to the earth's vicinity during times of geomagnetic disturbances are not adequately described by using static geomagnetic field models. These solar fluxes are also often distributed non uniformly in space, so that fluxes measured by satellites obtained at great distances from the earth and which sample large volumes of space around the earth cannot be used to predict fluxes locally at the earth's surface. We present here a method which uses the ground-level neutron monitor counting rates as adjoint sources of the flux in the atmosphere immediately above them to obtain solar-particle effective dose rates as a function of position over the earth's surface. We have applied this approach to the large September 29-30, 1989 ground-level event (designated GLE 42) to obtain the magnitude and distribution of the solar-particle effective dose rate from an atypically large event. The results of these calculations clearly show the effect of the softer particle spectra associated with solar particle events, as compared with galactic cosmic rays, results in a greater sensitivity to the geomagnetic field, and, unlike cosmic rays, the near-absence of a "knee" near 60 degrees geomagnetic latitude.