非线性波消波及波场分布研究

本文基于雷诺平均N-S方程,并结合RNG k-e方程建立了粘性数值波浪水槽,对不同波陡、不同相对水深、不同相对波高的非线性规则波的阻尼消波问题和波场分布进行研究。文中提出了两种描述消波区内部阻尼变化的阻尼函数,分别适用于小波陡情形和高波陡情形。研究结果表明,小波陡组消波区可设为一个波长,阻尼系数取10~4～10~5即可满足消波要求,计算结果与实验结果及造波理论吻合良好;高波陡组消波区可设为两个波长,阻尼系数取10~4～10~5亦可满足消波要求,计算结果与实验结果吻合良好。此外,当波陡较小时,波场内反射情况的小幅改变即可对整个波场造成影响,特别是当水深较浅时这种影响极为明显,需谨慎考虑。当波陡较大时,水波能量较高,整个波场沿水波传播方向可观测到明显的衰减现象,在具体试验中需进行考虑。     

Canonical theory of perturbations using eccentric anomaly as independent variable

It is shown in this paper how to build a canonical transformation of variables, so that the eccentric anomaly becomes the new independent variable. In the case of eccentric elliptical orbits it changes the equations of motion so, that they can be integrated analytically to any order of approximation comparatively easy.     

A study on the reference frame consistency in recent Earth gravitational models

All gravity field functionals obtained from an Earth gravitational model (EGM) depend on the underlying terrestrial reference frame (TRF), with respect to which the EGM’s spherical harmonic coefficients refer to. In order to maintain a coherent framework for the comparison of current and future EGMs, it is thus important to investigate the consistency of their inherent TRFs, especially when their use is intended for high precision studies. Following the methodology described in an earlier paper by Kleusberg (1980), the similarity transformation parameters between the associated reference frames for several EGMs (including the most recent CHAMP/GRACE models at the time of writing this paper) are estimated in the present study. Specifically, the differences between the spherical harmonic coefficients for various pairs of EGMs are parameterized through a 3D-similarity spatial transformation model that relates their underlying TRFs. From the least-squares adjustment of such a parametric model, the origin, orientation and scale stability between the EGMs’ reference frames can be identified by estimating their corresponding translation, rotation and scale factor parameters. Various aspects of the estimation procedure and its results are highlighted in the paper, including data weighting schemes, the sensitivity of the results with respect to the selected harmonic spectral band, the correlation structure and precision level of the estimated transformation parameters, the effect of the estimated differences of the EGMs’ reference frames on their height anomaly signal, and the overall feasibility of Kleusberg’s formulae for the assessment of TRF inconsistencies among global geopotential models.     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical calculations in the orbital determination of an artificial satellite for a long arc

Numerical methods are usually used for the computation of ephemerides with perturbations for the precise orbital determination of an artificial satellite. But their numerical stability will be encountered in a long arc. In this case the use the improved Encke special perturbation methods has been suggested. The results of this paper show that Encke's method does indeed have a certain effectiveness, but cannot yet completely resolve the numerical stability, and the more efficient method is to use the energy integral or its variational relation to control the growth of the along-track error in general numerical calculations so that the aim of stabilization can be achieved.     

FITS、BMP和SCR图象格式及相互转换

本文详细介绍图象的ＦＩＴＳ格式、ＢＭＰ格式和ＳＣＲ格式。ＦＩＴＳ格式已经是天文界的通用格式，几乎所有的天文软件包都支持这一格式，而ＢＭＰ格式在ＰＣ计算机上有广泛的运用，有大量ＰＣ软件支持ＢＭＰ图象的显示，处理和打印。云南天文台的ＳＣＲ格式是ＰＣ机采集ＣＣＤ图象所使用的格式。实现这三种格式的相互转换，就可以自由的将ＣＣＤ图象在ＰＣ机和工作站上进行各种处理。转换的包括图象头和图象数据的转换，其关键在于交换数据的高低字节。     

江苏潮滩沉积物激光粒度仪与移液管-筛析分析结果的对比

激光粒度仪的广泛应用带来了与历史数据的对比问题,因此需要建立激光粒度仪与早期分析结果之间的关系。根据江苏海岸潮滩沉积物样品的激光粒度仪和移液管-筛析法分析,对粒度参数、粒度组分等进行了对比,对两种方法的差异进行了分析。结果表明,对于江苏潮滩沉积物粒度参数中平均粒径的激光粒度仪与移液管-筛析分析结果之间有良好的线性关系;筛分法测得的粗颗粒物质较激光法偏少,而移液管法法测得的细颗粒物质较激光法偏多;将样品分类之后再进行两种方法所获粒度参数的回归分析,相关性得以提高,说明不同粒度组成的沉积物对分析结果的对比有不同的影响。两种方法之间的换算关系不仅与研究区域有关,而且与沉积物本身的粒度组成有关。     

PSO-RBFNN模型及其在岩土工程非线性时间序列预测中的应用

岩土工程受力变形演化是一个典型的非线性问题,其演化的高度非线性和复杂性,很难用简单的力学、数学模型描述,但可用粒子群优化径向基神经网络对岩土工程应力、位移非线性时间序列进行动态实时预测。网络径向基层的单元数通过均值聚类法确定后,所有其它参数：中心位置、形状参数、网络权值,均通过粒子群优化算法在全局空间优化确定。工程实例应用表明,该模型预测结果准确、精度高,有良好的应用前景。     

1.0GPa和常温至1100 ℃条件下角闪石斜长片麻岩的V_p变化:实验测量与理论计算

在1．0GPa、常温至1100℃条件下测量了角闪石斜长片麻岩的纵波速度（Vp）,并统计了不同温度实验产物中各种矿物的体积百分含量。结果显示,1．0GPa下,角闪石斜长片麻岩的Vp首先随温度升高（室温至700℃）缓慢降低约6％,然后（700—850℃）快速降低约6％,再（850～950℃）转而急剧升高15％～25％,最后（T〉950℃）又快速下降。实验产物分析表明在高温高压下α-石英-β石英相变和岩石部分熔融是岩石Vp异常变化的主要因素。由取样产物的矿物含量和弹性参数,计算了各温度条件下岩石的Vp得出与实验测量相同的波速-温度变化趋势,即Vp随温度升高先缓慢降低,接着快速降低后又急剧升高,最后又快速减小。实验测量和理论计算对比研究表明,通过高温高压下岩石中的物相变化观测结果进行岩石波速的计算,是检验岩石弹性波速测量结果和研究地球内部地震波结构的一种有效方法。     

影响低渗透油藏低速非线性渗流的实验研究

以膨胀性的粘土为主要研究对象,通过粘土浆体的流变性实验和填砂管渗流实验,并使用一种新研制的微流量测量工具,研究了低渗岩心在较低流速下的渗流特性;利用微尺度流动理论,Einstein粘度定律等进行了理论分析,对引起低渗透油藏低速非线性渗流的影响因素进行了初步的探索.结论认为,粘土的膨胀性以及其水化后形成浆体的流变性是引起低速非线性渗流的主要因素之一.