Simulation of the full two rigid body problem using polyhedral mutual potential and potential derivatives approach

Herein we investigate the coupled orbital and rotational dynamics of two rigid bodies modelled as polyhedra, under the influence of their mutual gravitational potential. The bodies may possess any arbitrary shape and mass distribution. A method of calculating the mutual potential’s derivatives with respect to relative position and attitude is derived. Relative equations of motion for the two body system are presented and an implementation of the equations of motion with the potential gradients approach is described. Results obtained with this dynamic simulation software package are presented for multiple cases to validate the approach and illustrate its utility. This simulation capability is useful both for addressing questions in dynamical astronomy and for enabling spacecraft missions to binary asteroid systems.     

搅拌桩复合地基平面模拟的简化方法探讨

水泥土搅拌桩(DMM)在高速公路工程中得到广泛应用,人们也希望能对该复合地基在路堤荷载下的地基的变形深入了解,其数值分析是一种重要的方法。但是,采用三维有限元进行数值模拟在高速公路工程中容易造成数组过大,计算时间并且收敛性比较差,因此把三维问题简化为二维问题是众望所归。在把三维问题简化二维平面问题简化计算时,最重要的是参数的折减。基于水泥土搅拌桩复合地基的特性,结合理论分析提出了一种参数折减的方法。并依据连盐高速公路实测的数据验证了该参数简化方法的可行性。     

The Data Depth-weight-kernel Estimation of the Model Error of Satellite Orbit Determination

It is an objective fact that there exists error in the satellite dynamic model and it will be transferred to satellite orbit determination algorithm, forming a part of the connotative model error. Mixed with the systematic error and random error of the measurements, they form the unitive model error and badly restrict the precision of the orbit determination. We deduce in detail the equations of orbit improvement for a system with dynamic model error, construct the parametric model for the explicit part of the model and nonparametric model for the error that can not be explicitly described. We also construct the partially linear orbit determination model, estimate and fit the model error using a two-stage estimation and a kernel function estimation, and finally make the corresponding compensation in the orbit determination. Beginning from the data depth theory, a data depth weight kernel estimator for model error is proposed for the sake of promoting the steadiness of model error estimation. Simulation experiments of SBSS are performed. The results show clearly that the model error is one of the most important effects that will influence the precision of the orbit determination. The kernel function method can effectively estimate the model error, with the window width as a major restrict parameter. A data depth-weight-kernel estimation, however, can improve largely the robustness of the kernel function and therefore improve the precision of orbit determination.     

Effect of a thermokarst lake on soil physical properties and infiltration processes in the permafrost region of the Qinghai-Tibet Plateau,China

Changes in the hydrological processes in alpine soil constitute one of the several key problems encountered with studying watershed hydrology and ecosystem stability against the background of global warming. A typically developing thermokarst lake was chosen as a subject for a study using model simulation based on observations of soil physical properties, infiltration processes, and soil moisture. The results showed that the selected thermokarst lake imposed certain changes on the soil infiltration processes and, with the degree of impact intensifying, the initial infiltration rate decreased. The greatest reduction was achieved in the area of moderate impact. However, the stable infiltration rate and cumulative infiltration gradually increased in the surface layer at a depth of 10 and 20 cm, both decreasing initially and then increasing, which is correlated significantly with soil textures. Moreover, the cumulative infiltration changed in line with steady infiltration rate. Based on a comparative analysis, the Horton model helps better understand the effect on the soil infiltration processes of the cold alpine meadow close to the chosen thermokarst lake. In conclusion, the formation of the thermokarst lake reduced the water holding capacity of the alpine meadow soil and caused the hydraulic conductivity to increase, resulting in the reduction of runoff capacity in the area of the thermokarst lake.     

Prediction of Water Temperature as Affected by a Pre‐Constructed Reservoir Project Based on MIKE11

Reservoir construction can lead to much more water stored in front of the dam and significantly increase heat storage capacity of the reservoir waters, thus resulting in different distribution pattern of water temperature in reservoir area compared to river. Especially for large reservoir, the obvious stratification of water temperature will appear in the reservoir with deeper water levels. Meanwhile, the low water temperature will be observed in the downstream river due to the operation of the reservoir. The vertical numerical simulation model for reservoir from MIKE 11 was used to predict the changes of water temperature of Wuxikou Reservoir to check the effects of the reservoir construction on water temperature. The water temperature prediction model was developed to simulate the water temperature of the reservoir and the discharged outflow water. The predicted results can contribute to assessing the feasibility of the pre‐constructed project based on the environmental influence of water temperature.     

复波场的实部和虚部对同时源波形反演的影响

高精度的成像需要准确的速度场信息,波形反演被认为是目前具有最高分辨率的速度反演方法之一.计算效率是目前全波形反演需要考虑的一个主要问题.为了很好地解决计算效率的问题,本文引入了一种高效的无串扰同时源反演方法,并详细介绍了其原理与计算流程.同时,基于此多震源同时反演方法,本文拓展出实波场反演、虚波场反演及复波场反演的三种反演策略,进而分析复波场的实部和虚部对全波形反演的影响.实验表明,相比于复波场全波形反演,无论是实波场反演还是虚波场反演,反演分辨率有所降低;相比于实波场反演,虚波场反演对初始模型的依赖性较小,目标函数的非线性较弱;最后,通过使用组合反演策略,即初始阶段采用虚波场反演,中后期阶段采用复波场反演,不仅可以降低反演的非线性,而且能够保证高精度建模.

     

A simplified method for simulation of strong ground motion using finite rupture model of the earthquake source

We present a simplified method to simulate strong ground motion for a realistic representation of a finite earthquake source burried in a layered earth. This method is based on the stochastic simulation method of Boore (Boore, D. M., 1983, Bull. Seism. Soc. Am. 73, 1865–1894) and the Empirical Greens Function (EFG) method of Irikura (Irikura, K., 1986, Proceedings of the 7th Japan Earthquake symposium, pp. 151–156). The rupture responsible for an earthquake is represented by several subfaults. The geometry of subfaults and their number is decided by the similarity relationships. For simulation of ground motion using the stochastic simulation technique we used the shapping window based on the kinetic source model of the rupture plane. The shaping window deepens on the geometry of the earthquake source and the propagation characteristics of the energy released by various subfaults. The division of large fault into small subfaults and the method for accounting their contribution at the surface is identical to the EGF. The shapping window has been modified to take into account the effect of the transmission of energy released form the finite fault at various boundaries of the layered earth model above the source. In the present method we have applied the correction factor to adjust slip time function of small and large earthquakes. The correction factor is used to simulate strong motion records having basic spectral shape of ² source model in broad frequency range. To test this method we have used the strong motion data of the Geiyo earthquake of 24th March 2001, Japan recorded by KiK network. The source of this earthquake is modelled by a simple rectangular rupture of size 24 × 15 km, burried at a depth of 31 km in a multilayered earth model. This rupture plane is divided into 16 rectangular subfaults of size 6.0 × 3.75 km each. Strong motion records at eight selected near-field stations were simulated and compared with the observed records in terms of the acceleration and velocity records and their response spectrum. The comparison confirms the suitability of proposed rupture model responsible for this earthquake and the efficacy of the approach in predicting the strong motion scenario of earthquakes in the subduction zone. Using the same rupture model of the Geiyo earthquake, we compared the simulated records from our and the EGF techniques at one near-field station. The comparison shows that this technique gives records which matches in a wide frequency range and that too from simple and easily accessible parameters of burried rupture.     

重庆地区夏季一次降水过程及增雨潜力的数值模拟分析

利用国家气象中心GRAPES人工增雨云系模式,选取2008年7月4日重庆地区一次降水过程进行数值模拟,分析了重庆地区降雨天气的水汽分布、云系宏微观分布、云中微物理转化和增雨潜力等特征。结果表明：本次降水大气过程中,重庆地区水汽含量极为丰富,水汽分布与地形分布呈明显的对应关系,低层水汽输送较大,整层水汽通量较高,有明显水汽辐合,云中液态水对地面降水影响很大。西南气流和地形共同作用为重庆地区液态水的形成提供了有利条件,在东北部山区迎风坡处大量水汽累积抬升,易形成丰富的液态水。重庆东北部地区水汽向云水转化较强,过冷液态水含量丰富,冰晶含量少,0℃层附近水汽垂直通量较大,降水效率较低,有较大的增雨潜力。     

利用ARC/INFO进行三维地形模型的建立、处理和动态模拟

介绍了利用 ARC/ INFO建立三维地形模型的方法 ,以及利用其新版本的部分新功能处理三维地形模型、进行动态模拟的方法 ,与以往处理三维模型必须借助外部三维图形开发软件包的方法相比 ,具有高速、灵活、可视性强、支持强大分析功能等优点。形成的三维模型可随时与 ARC/ INFO的其它数据模型交换信息 ,与外部进程实时通讯 ,为使用者提供了二维分析工具所不能提供的信息 ,为开发视景仿真系统等与表面分析、立体显示有关的应用系统提供了强有力的支持。     

锦屏电站坝厂区卸荷岩体力学参数研究

选择具有代表性的区块,根据结构面的性状和规模拟定模型层次和尺寸,建立岩体数值仿真模型,通过计算机数值仿真,采用卸荷非线性力学理论对锦屏电站坝厂区岩体进行卸荷岩体宏观力学特性及参数的数值分析;同时,用RMR法和BP神经网络法对岩体宏观力学特性及参数进行分析,综合确定岩体力学参数,为设计部门提供设计参考依据.