基于重力梯度杆和磁铁的小卫星三轴姿态控制

运用卫星低轨道两个主要环境力矩（重力梯度矩和地磁力矩）对圆轨道卫星三轴姿态进行被动控制。利用重力梯度矩实现卫星对地指向：卫星上的永久磁铁获取所需的地磁力矩,稳定偏航姿态。给出卫星的姿态分析,并给出仿真结果。从分析和仿真结果可以看出,此卫星具有结构简单、时刻对地定向、低轨道倾角时卫星姿态稳定精度较高的优点。     

固体导弹发动机药柱长期重力作用下的解析计算

为分析固体导弹发动机在长期重力作用下的位移、应力和应变,用应力函数推导了发动机药柱水平贮存时的弹性和粘弹性解析解,给出了相应的计算公式。根据Burgers模型表征的药柱力学特性,用弹性和粘弹性有限元分析法对发动机药柱进行了数值计算。算例分析表明,解析解与有限元解吻合较好,误差小于1%。     

AIUB-CHAMP02S: The influence of GNSS model changes on gravity field recovery using spaceborne GPS

The gravity field model AIUB-CHAMP02S, which is based on six years of CHAMP GPS data, is presented here. The gravity field parameters were derived using a two step procedure: In a first step a kinematic trajectory of a low Earth orbiting (LEO) satellite is computed using the GPS data from the on-board receiver. In this step the orbits and clock corrections of the GPS satellites as well as the Earth rotation parameters (ERPs) are introduced as known. In the second step this kinematic orbit is represented by a gravitational force model and orbit parameters.     

3D Imaging of the OH mesospheric emissive layer

A new and original stereo imaging method is introduced to measure the altitude of the OH nightglow layer and provide a 3D perspective map of the altitude of the layer centroid. Near-IR photographs of the OH layer are taken at two sites separated by a 645 km distance. Each photograph is processed in order to provide a satellite view of the layer. When superposed, the two views present a common diamond-shaped area. Pairs of matched points that correspond to a physical emissive point in the common area are identified in calculating a normalized cross-correlation coefficient (NCC). This method is suitable for obtaining 3D representations in the case of low-contrast objects. An observational campaign was conducted in July 2006 in Peru. The images were taken simultaneously at Cerro Cosmos (12°09′08.2″ S, 75°33′49.3″ W, altitude 4630 m) close to Huancayo and Cerro Verde Tellolo (16°33′17.6″ S, 71°39′59.4″ W, altitude 2272 m) close to Arequipa. 3D maps of the layer surface were retrieved and compared with pseudo-relief intensity maps of the same region. The mean altitude of the emission barycenter is located at 86.3 km on July 26. Comparable relief wavy features appear in the 3D and intensity maps. It is shown that the vertical amplitude of the wave system varies as exp (Δz/2H) within the altitude range Δz = 83.5–88.0 km, H being the scale height. The oscillatory kinetic energy at the altitude of the OH layer is comprised between 3 × 10⁻⁴ and 5.4 × 10⁻⁴ J/m³, which is 2–3 times smaller than the values derived from partial radio wave at 52°N latitude.     

Shear flow energy redistribution stipulated by the internal-gravity wavy structures in the dissipative ionosphere

The linear mechanism of generation, intensification and further nonlinear dynamics of internal gravity waves (IGW) in stably stratified dissipative ionosphere with non-uniform zonal wind (shear flow) is studied. In case of the shear flows the operators of linear problem are non-selfadjoint, and the corresponding Eigen functions – nonorthogonal. Thus, canonical – modal approach is of less use studying such motions. Non-modal mathematical analysis becomes more adequate for such problems. On the basis of non-modal approach, the equations of dynamics and the energy transfer of IGW disturbances in the ionosphere with a shear flow is obtained. Exact analytical solutions of the linear as well as the nonlinear dynamic equations of the problem are built. The increment of shear instability of IGW is defined. It is revealed that the transient amplification of IGW disturbances due time does not flow exponentially, but in algebraic – power law manner. The effectiveness of the linear amplification mechanism of IGW at interaction with non-uniform zonal wind is analyzed. It is shown that at initial linear stage of evolution IGW effectively temporarily draws energy from the shear flow significantly increasing (by an order of magnitude) own amplitude and energy. With amplitude growth the nonlinear mechanism turns on and the process ends with self-organization of nonlinear solitary, strongly localized IGW vortex structures (the monopole vortex, the transverse vortex chain or the longitudinal vortex street). Accumulation of these vortices in the ionospheric medium can create the strongly turbulent state.     

Lunar geophysical parameters inversion based on gravity/topography admittance and particle swarm optimization

This paper presents a FORTRAN computer program. The program as code will be used for lunar parameter inversions based on gravity/topography admittance. This will be done by assuming that the lunar lithosphere is modeled as a thin elastic spherical shell. The parameters discussed here include; load ratio, crustal thickness, subsurface load depth, crustal density and elastic lithosphere thickness. The admittance of the best-fitting model can be found through automatically adjusting misfits between one theoretical admittance and an observed one. The results in this paper indicate that this research’s theoretical model is reasonable for exploring the best-fitting parameters. In addition, this code is not only able to automatically and simultaneously calculate the global optimum solution of the parameters studied, but also performs well in computational speed. The code can be easily modified to include more parameter inversions; such as the inversion for subsurface density anomaly and the case of considering infilling material in some lunar mare basins.     

Application of novel force control strategies to enhance robotic abrasive belt grinding quality of aero-engine blades

Robotic belt grinding has emerged as a finishing process in recent years for machining components with high surface finish and flexibility.The surface machining consistency, however,is difficult to be guaranteed in such a process.To overcome this problem, a method of hybrid force-position control combined with PI/PD control is proposed to be applied in robotic abrasive belt grinding of complex geometries.Voltage signals are firstly obtained and transformed to force information with signal conditioning methods.Secondly, zero drift and gravity compensation algorithms are presented to calibrate the F/T transducer which is installed on the robot end-effector.Next, a force control strategy combining hybrid force-position control with PI/PD control is introduced to be employed in robotic abrasive belt grinding operations where the force control law is applied to the Z direction of the tool frame and the positon control law is used in the X direction of the tool frame.Then, the accuracy of the F/T transducer and the robotic force control system is analyzed to ensure the stability and reliability of force control in the robotic grinding process.Finally, two typical cases on robotic belt grinding of a test workpiece and an aero-engine blade are conducted to validate the practicality and effectiveness of the force control technology proposed.     

空间太阳能电站的准对日定向姿态

针对空间太阳能电站的俯仰姿态运动，提出一种能追踪太阳运动的准对日定向(QSP)姿态方案。此方案的太阳能电池阵列在万有引力梯度力矩的作用下，始终在垂直于太阳光的方向附近作幅值约为18.8°的振动，且几乎不需要姿态控制力矩。准对日定向姿态方案解决了大型太阳能电池阵列对日定向所需的巨大俯仰姿态控制力矩问题。准对日定向姿态的发电效率为对日定向姿态的97.3%，对Abacus空间太阳能电站而言每年可节省燃料约36791 kg。通过数值方法得到了准对日定向姿态的精确初始条件。随后，设计了比例-微分控制器，保证了系统存在初始姿态误差的条件下收敛到准对日定向姿态。最后研究了轨道、姿态和结构振动对准对日定向姿态的影响，并发现准对日定向姿态下的结构振动幅值比对日定向姿态减小约40倍。     

Study on dependence of hydrodynamic characteristics on gravity in a vertical upward channel bubbly flow

Since the bubbly flow has extensive applications in the space field, it is very necessary to comprehend the dependence of hydrodynamic characteristics on gravity. In this paper, the dependence of the microbubble distribution and the liquid turbulence modulation by the microbubbles on gravity was investigated in detail with the Euler–Lagrange method. The liquid Navier–Stokes equation was solved using direct numerical simulations (DNS), and the microbubble motion was tracked with Newtonian motion equation considering drag force, shear lift force, added mass force, pressure gradient force, and wall lift force. The coupling between the gas and the liquid phases regarded the interphase force as a momentum source term in the momentum equation of the liquid. The results showed that the phase profile and the turbulence modulation by the microbubbles strongly depend on the magnitude of gravity. When the influence of gravity is relatively weak, the microbubbles approximately uniformly disperse in the wide central region of the channel, and the average statistics of the liquid turbulence are almost not influenced due to the addition of the microbubbles. However, when the influence of gravity is comparatively important, the majority of the microbubbles accumulate near the wall of the channel, and the injection of the microbubbles modifies the profiles of the liquid average statistics.     

Regional gravity field model in Pakistan area from the combination of CHAMP,GRACE and ground data using least squares collocation: A case study

This study describes a methodology of recovery of the Earth’s gravity field from CHAMP and GRACE satellites data in Pakistan using least squares collocation (LSC) based downward continuation technique. The CHAMP height anomalies and GRACE gravity disturbances derived from the observed satellite data have been used in combination solution using LSC with observed gravity values at the Earth surface. The combined covariance functions of height anomalies and/or gravity disturbances at satellite altitudes and observed gravity anomalies at Earth surface have been used as the basis for combination and downward continuation solution. The variance of predicted gravity anomalies from GRACE gravity disturbances is relatively lower than the corresponding results of gravity anomalies from CHAMP height anomalies. This fact may be attributed partly to the amplification of noise and partly to the unstable inverse transformation process of height anomalies to gravity anomalies. The impact of data error variance has been studied in the context of smoothing and noise reduction in the final solution of downward continuation using least squares collocation. The raising of data error suppresses the noise and as a result a smooth final solution is obtained. The prediction results appear to be dependent on the quality of data and goodness of combined covariance function, which are fairly comparable for the CHAMP and GRACE data. The recovered gravity field from satellite data appears to contribute mainly to medium and long wavelength parts of total gravity field spectrum. Due to flexibility of data handling in least squares collocation, this procedure is applicable to any observable of gravity field being at different altitudes and with different data spacing.