利用Schamberg模型计算凸体的气动力矩

Schamberg的“超热力学自由分子流”粒子——表面作用模型已用于计算凸体的阻力系数。本文进一步给出按此模型计算凸体气动力矩的公式,并推荐反射束半顶角φ_0与入射角θ_1之间的函数表达式。     

《中国空间科学技术》1987年总目录

期号页数1洲防291盯4500月性 空间飞行器系统工程空间站发展的探讨任意偏心率卫星的大气阻力的短周期摄动空间站概述航天飞行器入轨道误差分析载人航天器轨道运行段的救生问题整星稳态温度的热网络分析方法地球引力场带谐调和项引起的卫星轨道摄动Q曰1上,.︸一刁自0 QU进任61火︸O咬1 LJ3 4 4 56 飞行动力学航天飞行器的两大数学问题三层壁喷管的不稳定导热用Nocilla模型计算卫星气动力系数和热流密度分布求解Navier一Stokes方程的指数型差分法用Schamberg方法计算卫星阻力系数四层壁的不稳定导热3 41匕6 测控技术和数据处理模式识别技…     

卫星姿控发动机高空羽流场工程分析

发展卫星姿控小推力发动机高空羽流场的工程估算方法,计算分析其主要燃气组分的分布特征.采用网格技术将发动机喷管出口截面剖分形成多独立点源,建立了非均匀出口流动条件下的分析模型.与自由分子单一点源模型相比,该模型的计算结果更接近实验值.对MBB 10 N双组元发动机多组分羽流场的研究表明,压强场和内能场的变化较剧烈;不同于轻分子组分,重分子组分主要集中于轴线附近,其密度场呈较强方向性.      

光压摄动对卫星姿态轨道耦合的影响分析

随着卫星对地测量精度要求的不断提高, 对卫星轨道的精度要求也随之提高. 目前Topex, Jason-1, Jason-2等一系列海洋测地卫星的轨道计算精度已经达到厘米量级, 相应对卫星动力学模型的要求也越来越精细. 以Topex海洋测地卫星为背景, 考虑卫星帆板有规律的运动, 将其几何形状简化为高精度轨道计算中比较通用的Boxing-Wing模型, 计算了Topex卫星的Boxing-Wing模型在轨运行中受到的太阳光压力及光压力矩. 考虑卫星姿态和轨道耦合的情况下, 计算了太阳光压力及光压力矩对Topex卫星轨道半长轴和卫星姿态的影响. 通过一个轨道周期的计算可知, 光压对卫星轨道半长轴的影响大约为9cm, 对卫星滚动角和俯仰角的影响在6°左右, 因此, 在高精度的轨道计算和姿态控制中这个影响是应该考虑的.      

卫星全景象片的几何校正

本文给出了卫星全景相片的几何校正模型。该模型由三部分组成:1.推导了卫星全景相片的几何校正模型;2.给出了计算卫星姿态运动和卫星轨道运动的计算公式;3.用系统辨识方法估算卫星姿态运动公式。     

几何和流动参数对上翘后体阻力系数的影响

利用N-S方程对9个上翘后体模型进行了气动力计算,主要研究后体几何参数和流动参数对上翘后体阻力系数的影响.研究结果表明, 后体的压差阻力系数分别随上翘角、收缩比的增加及迎角的减小而明显增加;后体摩擦阻力系数分别随后体的长细比的增加和雷诺数的减小而增加;后体越扁平,其压差阻力系数越大;在跨音速时,波阻系数也与上翘角有关,上翘角增加会导致波阻系数进一步加大.     

基于双频GPS观测的简化动力学最小二乘批处理精密定轨

星载双频GPS载波相位和伪距观测量已成为低轨卫星获取精确三维位置和速度信息的主要方式. 本文以非差消电离载波相位和伪距组合作为观测量,应用简化动力学最小二乘批处理方法进行地球低轨卫星的精密定轨,并给出完整定轨流程. 采用逐段常量的经验加速度对动力学模型误差进行补偿,描述了经验加速度敏感矩阵及稀疏带状矩阵求逆的有效计算方法. 利用GRACE-A卫星GPS观测数据对定轨位置精度进行分析,结果显示,三维位置定轨精度优于5cm,经验加速度在径向、切向和法向上的补偿水平不超过40nm·s-2,大气阻力系数和辐射光压系数的估计值符合物理实际,星载接收机钟差大致呈线性并具有短周期小波动.      

仿壁虎刚毛阵列对卫星表面吸附能力模型与计算

壁虎脚的粘着机理对航天机器人脚掌的研制和开发具有重要启发意义.通过单根刚毛与光滑表面黏着力模型与计算、卫星粗糙表面分布模型与接触概率与计算、刚毛阵列与卫星粗糙表面黏着力模型与计算等，分析仿壁虎刚毛阵列对卫星表面的吸附能力，发现在一定条件下仿壁虎刚毛阵能够在太空失重环境中黏着在卫星表面，对航天舱内爬行机器人、管道机器人的研制和开发具有重要的启示意义.     

一种敏捷卫星积分时间计算方法

敏捷卫星成像过程中姿态变化灵活,积分时间变化剧烈,传统卫星积分时间计算模型已无法满足精度要求。因此,提出一种改进的积分时间计算模型,在星地斜距计算和摄影点地速计算两方面作了改进。星地斜距计算时利用地表数字高程模型修正理论的斜距值,摄影点地速计算中增加姿态角速度对卫星本体坐标系角速度矢量的影响,并设计了整个算法在星上数管分系统中的实现方案。最后在星载嵌入式环境下进行了仿真计算,并与地面仿真计算结果进行了比对分析。两者比对偏差在0.01‰以内,可以满足相机成像质量要求。     

摆镜摆动对星载CCD相机成像的影响

利用Ｌａｇｒａｎｇｅ方程,建立了摆镜摆动与卫星姿态振动相互影响的仿真模型,并用数值计算方法计算了卫星的姿态振动,仿真分析了对星载线阵ＣＣＤ成像质量的影响。     

Modeling satellite drag coefficients with response surfaces

Satellite drag coefficients are a major source of uncertainty in predicting the drag force on satellites in low Earth orbit. Among other things, accurately predicting the orbit requires detailed knowledge of the satellite drag coefficient. Computational methods are an important tool in computing the drag coefficient but are too intensive for real-time and predictive applications. Therefore, analytic or empirical models that can accurately predict drag coefficients are desired. This work uses response surfaces to model drag coefficients. The response surface methodology is validated by developing a response surface model for the drag coefficient of a sphere where the closed-form solution is known. The response surface model performs well in predicting the drag coefficient of a sphere with a root mean square percentage error less than 0.3% over the entire parameter space. For more complex geometries, such as the GRACE satellite, the Hubble Space Telescope, and the International Space Station, the model errors are only slightly larger at about 0.9%, 0.6%, and 1.0%, respectively.     

Drag and energy accommodation coefficients during sunspot maximum

Conditions appropriate to gas-surface interactions on satellite surfaces in orbit have not been successfully duplicated in the laboratory. However, measurements by pressure gauges and mass spectrometers in orbit have revealed enough of the basic physical chemistry that realistic theoretical models of the gas-surface interaction can now be used to calculate physical drag coefficients. The dependence of these drag coefficients on conditions in space can be inferred by comparing the physical drag coefficient of a satellite with a drag coefficient fitted to its observed orbital decay. This study takes advantage of recent data on spheres and attitude stabilized satellites to compare physical drag coefficients with the histories of the orbital decay of several satellites during the recent sunspot maximum. The orbital decay was obtained by fitting, in a least squares sense, the semi-major axis decay inferred from the historical two-line elements acquired by the US Space Surveillance Network. All the principal orbital perturbations were included, namely geopotential harmonics up to the 16th degree and order, third body attraction of the Moon and the Sun, direct solar radiation pressure (with eclipses), and aerodynamic drag, using the Jacchia-Bowman 2006 (JB2006) model to describe the atmospheric density. After adjusting for density model bias, a comparison of the fitted drag coefficient with the physical drag coefficient has yielded values for the energy accommodation coefficient as well as for the physical drag coefficient as a function of altitude during solar maximum conditions. The results are consistent with the altitude and solar cycle variation of atomic oxygen, which is known to be adsorbed on satellite surfaces, affecting both the energy accommodation and angular distribution of the reemitted molecules.     

Drag coefficient modeling for grace using Direct Simulation Monte Carlo

Drag coefficient is a major source of uncertainty in predicting the orbit of a satellite in low Earth orbit (LEO). Computational methods like the Test Particle Monte Carlo (TPMC) and Direct Simulation Monte Carlo (DSMC) are important tools in accurately computing physical drag coefficients. However, the methods are computationally expensive and cannot be employed real time. Therefore, modeling of the physical drag coefficient is required. This work presents a technique of developing parameterized drag coefficients models using the DSMC method. The technique is validated by developing a model for the Gravity Recovery and Climate Experiment (GRACE) satellite. Results show that drag coefficients computed using the developed model for GRACE agree to within 1% with those computed using DSMC.     

操纵面加装Gurney襟翼对无人机纵向气动特性的影响

通过风洞测力实验研究了矩形Gurney襟翼对多操纵面飞翼布局无人机纵向气动特性的影响.升降副翼和襟副翼加装Gurney襟翼都会增加全机升力系数以及低头力矩系数.通过比较发现,升降副翼加装Gurney襟翼的增升效果好于襟副翼,外侧操纵面控制效果好于内侧.升降副翼和襟副翼分别组合加装Gurney襟翼会进一步提高增升效果,并且升力系数随着Gurney襟翼高度的增加而增加.鸭翼上加装Gurney襟翼可以减小全机低头力矩系数,因此鸭翼和升降副翼、襟副翼组合加装Gurney襟翼在线性段提高全机升力系数的同时,可以保持低头力矩系数增量基本为0,该特性对改善飞行器性能尤为有利.     

回转体倾斜入水空泡及弹道特性实验

针对回转体低速倾斜入水过程空泡的生成机理和演化特性开展研究。通过不同入水速度和入水角度回转体倾斜入水对比实验,研究了入水速度和入水角度对入水空泡、回转体速度、俯仰角及阻力系数的影响规律。研究结果表明:空泡深闭合时刻的空泡长度、最大空泡直径和空泡闭合点深度随着回转体入水速度的增加而增大;空泡深闭合发生时间随着入水速度的增加而小幅增大;同一入水时刻回转体俯仰角和俯仰角速度均随入水速度的增加而减小;随着回转体入水角度的增加,同一时刻入水空泡直径和回转体速度逐渐减小,阻力系数逐渐增大。      

上游双扰动对圆柱气动特性影响的数值模拟

通过用商用计算流体力学软件FLUENT 6.0中的SIMPLE方法求解二维非定常层流模型,对雷诺数Re=200时上游双扰动体对下游圆柱气动特性影响进行了数值模拟研究.计算结果表明当双扰动体相隔较近(H/d=1.1,1.33)时,其后形成单一的Karman涡街,它们对下游圆柱的干扰作用类似于单扰动体时的情况,但是流动结构发生变化的临界间距Lcri较之有效宽度相同的单扰动体的要大,因此可以推断出在该种布置下双扰动体和圆柱绕流会出现2种流动模式.在空穴流动模式下面,下游圆柱的阻力和升力脉动最多减小94%和79%,系统阻力(包括圆柱和扰动体)也能减小33%.而当双扰动体尾流为双Karman涡街时(双扰动体相隔较远H/d=3.33,5.0),圆柱迎风面部分区域仍受到来流的冲击作用,气动力减小有限.     

尾翼对民机后体流动特性的影响

在1.5 m低速风洞中采用表面油流、烟线/激光片光显示和部件测力等实验方法对2种民机布局后体的绕流流型和气动特性进行了研究.实验的风速为40 m/s,雷诺数为1.191×105,迎角范围为-10°~20°.实验结果表明:随着迎角由负到正的变化,2种布局模型后体绕流都经历了上分离—无分离—下分离流型的变化,分离流型属开式分离并由尾部起始向机身发展;小迎角下2种布局模型后体阻力的差别主要由尾翼的摩擦阻力引起,而较大迎角下尾翼绕流分离引起的压差阻力起主要作用;垂尾分离随迎角增加而不断减弱,从而导致正常式布局后体阻力随迎角变化表现为非对称.      

底凹结构减阻效应数值分析

为研究弹丸底凹结构的减阻机理,使用三维定常CFD方法对M910弹丸的流场特性进行了数值模拟。给出了零升阻力系数随马赫数的变化规律,所得结果与实验数据符合很好。在此基础上,为M910弹丸引入底凹结构并进行数值模拟。对比了不同弹底结构的底部流场特性,对底凹结构减阻效应的产生机理进行了分析。结果表明:亚声速下,底凹结构在底凹腔体内引入了高压“死水区”,并以“屈从”的流体边界代替了原固体底面,从而改变了尾部涡街的形成位置、形状和强度,最终增大底部压力,减小弹丸阻力;跨声速下,由于尾部涡街远离弹丸底面,固体底面与流体边界面的作用相同,使得底凹不再具有减阻效果;超声速下,底凹结构的减阻机理与底排弹丸减阻机理类似,即底凹结构中的流体为弹丸底部回流区添加质量从而达到减阻作用。      

锯齿形格尼襟翼气动性能的实验研究

用低速风洞测力试验和襟翼处绕流的PIV测量试验研究锯齿形格尼襟翼在不同偏角下的增升效益.结果表明:锯齿形格尼襟翼能明显提高翼型的升力系数和大升力系数下的翼型升阻比,对于给定的襟翼弦长,存在一个最佳的襟翼偏角,在此偏角下,翼型升阻比不仅在大升力系数下有明显提高,而且在中小升力系数时升阻比也有一定的提高.PIV测量表明从锯齿形格尼襟翼的齿边向上卷起的流向涡使上翼面后部气流向翼型表面吸附,推迟了上翼面气流的分离.     

大气阻力参数修正对低轨空间目标轨道预报精度的改进

对于低轨空间目标, 大气阻力是影响轨道预报精度的主要摄动力. 本文提出了一种 基于空间环境数据和神经网络模型的空间目标大气阻力参数修正方法, 基于目 标的历史两行元根数, 通过模拟得到外推一天轨道预报中预报结果与观测数据 符合最好的阻力调制系数, 分析表明其与太阳F_10.7指数和地磁Ap指数具有很好的相关性. 根据已有数据, 构建神经网络模型, 实现对阻力调制系数 的补偿计算, 从而改进低轨目标外推一天的轨道预报. 结果表明, 神经网络模 型相比两行元根数能够更及时地对空间环境变化进行响应. 将该方案应用于天 宫一号和国际空间站的外推一天轨道预报, 验证了方案的正确性和普适性, 对 地磁扰动引起的较大预报误差改进效果更好, 误差能够降低50%~60%; 平均而言, 预报精度可以提高约30%, 改进成功率达到80%左右.