A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

有关边界层问题中一类奇异积分方程新的存在性结果

对边界层理论新结果中出现的一类奇异积分方程w（t）=∫1 t（1-s）（λ＋λs＋s）/w（s）ds＋（1-t）∫t 0s/w（s）ds,t∈（0,1）进行讨论，并得出了上述方程在λ∈（-1/2,0）上正解存在性的新结果。     

地震前后垂直形变场动态演化的量化指标

提出一种表达垂直形变场动态演化过程的量化指标——区域应变率、应变集中度。在此基础上，对南北地震带各水准监测区近30年的垂直形变资料进行了实际计算．并结合具体震例进行了对比研究。结果表明：量化指标在一定程度上能够反映地震前后形变场的动态演化过程和地壳运动状态，对地震的中长期预报有一定的积极意义。     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

Boundary curves for families of planar orbits

For monoparametric familiesf(x,y)=c of planar orbits, created by a planar potentialV(x,y), we introduce the notion of the family boundary curves (FBC). All members of the familyf(x,y)=c are traced in an allowable region of thexy plane, defined by the corresponding FBC, with total energyE=E(c) varying along the family. Family boundary curves are also found for two-parametric familiesf(x,y,b)=c. The relation of equilibrium points and asymptotic orbits, possibly possessed by the potentialV(x,y), to be FBC is studied.     

差分观测与伪距观测的定位精度

严格论证了在对空间或地面目标的定位中,利用伪距测量观测模型的几何因子大于用台站间一次差分观测模型的几何因子.此外,台站间完全差分测量观测的定位精度优于不完全测量情形的定位精度.这些结果给人们在实际工作中采取何种定位方式提供了可靠的理论依据.     

Ultra-Deep Optical Spectroscopy with PMAS. Using the Nod-and-Shuffle Technique

PMAS, the Potsdam Multi-Aperture Spectrophotometer, is a new integral field spectrograph in the optical, which is optimized for good transmission and high image quality from 350 nm to 1 μm. We present our plan to implement a CCD charge-shuffle mode to allow for beam switching with a very high degree of sky subtraction accuracy for faint object 3-D spectroscopy. This revised version was published online in July 2006 with corrections to the Cover Date.     

The origin of the correlation between the spin parameter and the baryon fraction of galactic disks

The puzzling correlation between the spin parameter λ of galactic disks and the disk-to-halo mass fraction f _disk is investigated. It is shown that such a correlation arises naturally from large uncertainties in determining the virial masses of dark matter halos. This result indicates that halo properties derived from fits to observed rotation curves are still very uncertain. An analysis of λ versus f _disk as function of the adopted halo virial mass shows that for reasonable halo concentrations f _disk ≈ 0.01 - 0.07 which is significantly smaller than the universal baryon fraction. Most of the available gas either never settled into the galactic disks or was ejected subsequently. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Homoclinic Tangle of A 1:2 Resonance in a 2-D Hamiltonian System

We study in great detail the geometry of the homoclinic tangle, with respect to the energy, corresponding to an unstable periodic orbit of type 1:2, on a surface of section representing a 2-D Hamiltonian system. The tangle consists of two resonance areas, in contrast with the tangles of type-l or -{l, m, k, x = 0} considered in previous studies, that consist of only one resonance area. We study the intersections of the inner and outer lobes of the same resonance area and of the two resonance areas. The intersections of the lobes follow certain rules. The detailed study of these rules allows us to derive quantitative relations about the number of intersections and to understand the complex behavior of the higher order lobes by studying the lower order lobes. We find 1st, 2nd, 3rd, etc. order intersections formed by lobes making 1, 2, 3, etc. turns around an island. After a sufficiently high order of iterations a lobe may intersect its image and thus produce a Poincaré recurrence. Numerical results for a wide interval of energies are presented. The number of intersections changes through tangencies. In any finite interval of the energy between two tangencies of 1st order, an infinite number of higher order tangencies occur and thus, according to the Newhouse theorem, there exist nearby islands of stability.