Hamilton系统数值计算的新方法   总被引：7，自引：0，他引：7  

廖新浩  刘林 《天文学进展》1996,14(1):3-11

系统地介绍了近年来对Ｈａｍｉｌｔｏｎ系统数值计算新建立的辛算法和线性对称多步法，并对它们在动力天文中的应用作了一简要回顾。  相似文献   

Consolidation of double‐layered ground with vertical drains     

Xiao‐Wu Tang  Katsutada Onitsuka 《国际地质力学数值与分析法杂志》2001,25(14):1449-1465

Vertical drains are usually installed in subsoil consisting of several layers. Due to the complex nature of the problem, over the past decades, the consolidation properties of multi‐layered ground with vertical drains have been analysed mainly by numerical methods. An analytical solution for consolidation of double‐layered ground with vertical drains under quasi‐equal strain condition is presented in this paper. The main steps for the computation procedure are listed. The convergence of the series solution is discussed. The comparisons between the results obtained by the present analytical method and the existing numerical solutions are described by figures. The orthogonal relation for the system of double‐layered ground with vertical drains is proven. Finally, some consolidation properties of double‐layered ground with vertical drains are analysed. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

On the stability of particular solutions of the singly averaged Hill problem     

M. A. Vashkov’yak  N. M. Teslenko 《Astronomy Letters》2005,31(12):844-852

We consider the particular solutions of the evolutionary system of equations in elements that correspond to planar and spatial circular orbits of the singly averaged Hill problem. We analyze the stability of planar and spatial circular orbits to inclination and eccentricity, respectively. We construct the instability regions of both particular solutions in the plane of parameters of the problem.  相似文献   

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction     

M. Zeghal  U. El Shamy 《国际地质力学数值与分析法杂志》2004,28(14):1361-1383

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

The formation and stability of the configuration of the planetary system HD 126611     

Zhang Xu-dong  Zhou Ji-lin 《Chinese Astronomy and Astrophysics》2006,30(4):420-430

Recent Doppler velocity measurements have revealed the existence of two planets orbiting the star HD 12661 on medium-eccentricity orbits. The inner planet has a period of 263.6 d and a mass of 2.3 Jupiter masses, and the outer planet has a period of 1444.5d and a mass of 1.57 Jupiter masses. The stability of this system requires the two planets to be in a state of mean-motion orbit resonances. By numerical method we have studied the orbit migration and stability of the system in its early ages under the action of the proto-stellar disk, and calculated the probabilities of the planets being captured into the mean -motion resonances during their migrations. It is found that at present the two planets are probably situated at the edge of the 11:2 mean-motion resonance and are in chaotic motions. This result may be helpful to clarify the arguments on the present configuration. Besides, it is indicated that very probably, after the formation of the system, the gaseous disk has almost disappeared before the planets migrated to the present configuration.  相似文献   

On the orbits and masses of the satellites of the Pluto-Charon system     

Man Hoi Lee  S.J. Peale 《Icarus》2006,184(2):573-583

Two small satellites of Pluto, S/2005 P1 (hereafter P1) and S/2005 P2 (hereafter P2), have recently been discovered outside the orbit of Charon, and their orbits are nearly circular and nearly coplanar with that of Charon. Because the mass ratio of Charon-Pluto is ∼0.1, the orbits of P2 and P1 are significantly non-Keplerian even if P2 and P1 have negligible masses. We present an analytic theory, with P2 and P1 treated as test particles, which shows that the motion can be represented by the superposition of the circular motion of a guiding center, the forced oscillations due to the non-axisymmetric components of the potential rotating at the mean motion of Pluto-Charon, the epicyclic motion, and the vertical motion. The analytic theory shows that the azimuthal periods of P2 and P1 are shorter than the Keplerian orbital periods, and this deviation from Kepler's third law is already detected in the unperturbed Keplerian fit of Buie and coworkers. In this analytic theory, the periapse and ascending node of each of the small satellites precess at nearly equal rates in opposite directions. From direct numerical orbit integrations, we show the increasing influence of the proximity of P2 and P1 to the 3:2 mean-motion commensurability on their orbital motion as their masses increase within the ranges allowed by the albedo uncertainties. If the geometric albedos of P2 and P1 are high and of order of that of Charon, the masses of P2 and P1 are sufficiently low that their orbits are well described by the analytic theory. The variation in the orbital radius of P2 due to the forced oscillations is comparable in magnitude to that due to the best-fit Keplerian eccentricity, and there is at present no evidence that P2 has any significant epicyclic eccentricity. However, the orbit of P1 has a significant epicyclic eccentricity, and the prograde precession of its longitude of periapse with a period of 5300 days should be easily detectable. If the albedos of P2 and P1 are as low as that of comets, the large inferred masses induce significant short-term variations in the epicyclic eccentricities and/or periapse longitudes on the 400-500-day timescales due to the proximity to the 3:2 commensurability. In fact, for the maximum inferred masses, P2 and P1 may be in the 3:2 mean-motion resonance, with the resonance variable involving the periapse longitude of P1 librating. Observations that sample the orbits of P2 and P1 well on the 400-500-day timescales should provide strong constraints on the masses of P2 and P1 in the near future.  相似文献   

Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application     

Daichao Sheng  David W. Smith  Scott W. Sloan  Antonio Gens 《国际地质力学数值与分析法杂志》2003,27(9):767-790

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

Nonlinear impact monitoring: line of variation searches for impactors     

Andrea Milani  Steven R. Chesley  Giacomo Tommei 《Icarus》2005,173(2):362-384

  相似文献   

CARBON DYNAMICS OF WETLAND IN THE SANJIANG PLAIN   总被引：2，自引：0，他引：2  

SONG Chang-chun    WANG Yi-yong  WANG Yue-si  YAN Bai-xing WANG De-xuan  ZHAO Zhi-chun  LOU Yan-jing 《中国地理科学(英文版)》2003,13(3):228-231

1INTRODUCTIONWetlandsplayanimportant roleintheprocessofcar-bonstorage.Thetotalcarbonstoredindifferentkindsofwetlandsisabout15%-35%ofthetotalcarboninthegloballandsoils(POSTetal.,1982;GORHAM,1991).Inaddition,wetlandsaresignificantnaturalsources fortheatmospheric CH4 (MOORE,1994).It isestimatedthatabout110×1012gCH4 originates fromanaerobicdecompositioninthenaturalwetlands,CH4 emission fromthenaturalwetlandsis15%-30%oftheglobalCH4 emission andtheCH4 emission from thepeat land at hi…  相似文献   

Use of multi‐platform,multi‐temporal remote‐sensing data for calibration of a distributed hydrological model: an application in the Arno basin,Italy     

Lorenzo Campo  Francesca Caparrini  Fabio Castelli 《水文研究》2006,20(13):2693-2712

Images from satellite platforms are a valid aid in order to obtain distributed information about hydrological surface states and parameters needed in calibration and validation of the water balance and flood forecasting. Remotely sensed data are easily available on large areas and with a frequency compatible with land cover changes. In this paper, remotely sensed images from different types of sensor have been utilized as a support to the calibration of the distributed hydrological model MOBIDIC, currently used in the experimental system of flood forecasting of the Arno River Basin Authority. Six radar images from ERS‐2 synthetic aperture radar (SAR) sensors (three for summer 2002 and three for spring–summer 2003) have been utilized and a relationship between soil saturation indexes and backscatter coefficient from SAR images has been investigated. Analysis has been performed only on pixels with meagre or no vegetation cover, in order to legitimize the assumption that water content of the soil is the main variable that influences the backscatter coefficient. Such pixels have been obtained by considering vegetation indexes (NDVI) and land cover maps produced by optical sensors (Landsat‐ETM). In order to calibrate the soil moisture model based on information provided by SAR images, an optimization algorithm has been utilized to minimize the regression error between saturation indexes from model and SAR data and error between measured and modelled discharge flows. Utilizing this procedure, model parameters that rule soil moisture fluxes have been calibrated, obtaining not only a good match with remotely sensed data, but also an enhancement of model performance in flow prediction with respect to a previous calibration with river discharge data only. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献