A complex variable solution for stress and displacement field around a lined circular tunnel at great depth

By applying the series expansion technique in the complex variable method established by Muskhelishvili, the plane elasticity problem for the stress and displacement field around a lined circular tunnel in conjunction with the consideration of misfit and interaction between the liner and the surrounding geomaterial is dealt with. The tunnel is assumed to be driven in a homogeneous and isotropic geomaterial. The coefficients in the Laurent series expansion of the stress functions are determined. The complex potentials in the liner and the surrounding geomaterial are explicitly derived, respectively. As an example, the case of a lined circular tunnel located in an isotropic initial stress field but subjected to uniform internal pressure is numerically considered. Numerical results indicate that the installation of tunnel liner can reduce the influences of the tunnel excavation on the in situ displacement and stress fields. However, the relative thickness and rigidity of the liner should be in an appropriate range. In addition, the effect of the tunnel excavation upon the displacement field is more significant than that upon the stress field. As far as the stress field in the surrounding geomaterial is concerned, when the ratio between the cover depth of tunnel and the tunnel radius is larger than 5, the results for the stress field in the paper are applicable. When the ratio between the tunnel depth and the tunnel radius is larger than 20, the results are applicable for the displacement field. Copyright © 2008 John Wiley & Sons, Ltd.     

Vertical dynamic response of a rigid foundation embedded in a poroelastic soil layer

A simplified analytical method is presented for the vertical dynamic analysis of a rigid, massive, cylindrical foundation embedded in a poroelastic soil layer. The foundation is subjected to a time‐harmonic vertical loading and is perfectly bonded to the surrounding soil in the vertical direction. The soil underlying the foundation base is represented by a single‐layered poroelastic soil based on rigid bedrock while the soil at the side of the foundation is modeled as an independent poroelastic layer composed of a series of infinitesimally thin layers. The behavior of the soil is governed by Biot's poroelastodynamic theory and its governing equations are solved by the use of Hankel integral transform. The contact surface between the foundation base and the soil is smooth and fully permeable. The dynamic interaction problem is solved following standard numerical procedures. The accuracy of the present solution is verified by comparisons with the well‐known solutions obtained from other approaches for both the elastodynamic interaction problem and poroelastodynamic interaction problem. Numerical results for the vertical dynamic impedance and response factor of the foundation are presented to demonstrate the influence of nondimensional frequency of excitation, soil layer thickness, poroelastic material parameters, depth ratio and mass ratio on the dynamic response of a rigid foundation embedded in a poroelastic soil layer. Copyright © 2008 John Wiley & Sons, Ltd.     

Extensive genetic divergence among Diptychus maculatus populations in northwest China

D. maculates is a kind of specialized Schizothoracinae fish has been locally listed as a protected animal in Xinjiang Province,China. Ili River located in north of Tianshan Mountain and Tarim River located in north of Qinghai-Tibetan Plateau were two main distribution areas of this fish. To investigate the genetic diversity and genetic structure of D. maculates,four populations from Tarim River system and two populations from Ili River system were collected in this study. A 570-bp sequence of the control region was obtained for 105 specimens. Twenty-four haplotypes were detected from six populations,only Kunes River population and Kashi River population shared haplotypes with each other. For all the populations examined,the haplotype diversity(h) was 0.904 8±0.012 6,nucleotide diversity(π) was 0.027 9±0.013 9,and the average number of pairwise nucleotide differences(k) was 15.878 3±7.139 1. The analysis of molecular variance(AMOVA) showed that 86.31% of the total genetic variation was apportioned among populations,and the variation within sampled populations was 13.69%. Genetic differences among sampled populations were highly significant. F st statistical test indicated that all populations were significantly divergent from each other(P 0.01). The largest F st value was between Yurungkash River population and Muzat River population,while the smallest F st value was between Kunes River population and Kashi River population. NJ phylogenetic tree of D-loop haplotypes revealed two main clades. The neutrality test and mismatch distribution analysis suggested that the fish had went through a recent population expansion. The uplift of Tianshan Mountain and movement of Qinghai-Tibetan Plateau might contribute to the wide genetic divergence of D. maculates in northwest China.     

20世纪南极菲尔德斯半岛毛皮海豹食谱比例变化的同位素混合模型比较分析

稳定同位素比值分析已成功用于示踪动物的食物来源研究,近年来已发展了多个同位素混合模型用于计算混合物中多来源物质的组成比例。本文以过去百年来南极菲尔德斯半岛南极毛皮海豹(Arctocephalus gazella)毛发的C、N同位素数据为基础,运用同位素混合模型中的欧几里得距离法、线性模型和贝叶斯模型进行计算、分析和比较,得到毛皮海豹食谱中南极磷虾(Euphausia superb)、南极电灯鱼(Electrona antarctica)和尼氏裸灯鱼(Gymnoscopelus nicholsi)所占比例随时间的变化趋势。各模型计算得出的结果可分为两组,分别是欧几里得距离法公式(1)和(2)的结果,表现为近百年来毛皮海豹食谱中Euphausia superb比例与其毛发稳定N同位素比值呈同步上升的趋势;第二组为欧几里得距离法公式(3)、线性模型和贝叶斯模型的结果,Euphausia superb比例与毛皮海豹毛稳定N同位素比值呈反相关关系。综合对比研究区域气候、海冰环境变化背景和已有的南极磷虾种群密度调查资料分析,第二组模型结果与实际情况较为相符,可用于讨论食物比例的时间变化趋势。近百年来,菲尔德斯半岛毛皮海豹食谱中Euphausia superb的比例显著下降,很可能是由于区域气候快速变暖和海冰的消退造成的。本文的计算结果指示西南极半岛气候变暖已显著影响到区域海洋食物链变化。     

Seismic data reconstruction based on iterative linear expansion of thresholds

Based on the compressive sensing,a novel algorithm is proposed to solve reconstruction problem under sparsity assumptions.Instead of estimating the reconstructed data through minimizing the objective function,the authors parameterize the problem as a linear combination of few elementary thresholding functions,which can be solved by calculating the linear weighting coefficients.It is to update the thresholding functions during the process of iteration.The advantage of this method is that the optimization problem only needs to be solved by calculating linear coefficients for each time.With the elementary thresholding functions satisfying certain constraints,a global convergence of the iterative algorithm is guaranteed.The synthetic and the field data results prove the effectiveness of the proposed algorithm.     

大棚蔬菜最低气温预报及冻害防御技术

在大棚蔬菜低温冻害防御中,通过搭建拱棚加盖不同层塑料内膜的大棚温度调控模式应用广泛。通过对不同层塑料内膜覆盖的大棚小气候环境温度观测,分析了棚内外日最低气温的相关关系,建立了棚内小气候环境日最低气温的一元线性回归预测模型,应用气象台天气预报中的日最低气温预报开展了大棚小气候环境的日最低气温预报,根据预报结果和蔬菜冻害指标,确定大棚增温防冻需要加盖塑料内膜的层数,提出大棚蔬菜内膜覆盖方式的冻害防御对策,为设施农业提供安全生产保障依据。     

The evolution of the z distribution of normal neutron stars in the Galaxy

Under the two initial 1‐D one parameter velocity distribution forms (one is normal, the other is exponential), the z direction scale height evolution of normal neutron stars in the Galaxy is studied by numerical simulation. We do statistics for the cases at different time segments, also do statistics for the cumulative cases made of each time segment. The results show in the cumulative cases the evolution curves of the scale heights are smoother than in the each time segment, i.e., the cumulation improve the signal‐to‐noise ratio. Certainly the evolution cases are different at different Galactic disk locations, which also have very large difference from the average cases in the whole disk. In the initial stages of z evolution of normal neutron stars, after the beginning transient states, the cumulative scale heights increase linearly with time, and the cumulative scale height increasing rates have linear relationship with the initial velocity distribution parameters, which have larger fluctuation in the vicinity of the Sun than in the whole disk. We utilize the linear relationship of the cumulative scale height increasing rates vs. the initial velocity distribution parameters in the vicinity of the Sun to make comparison with the observation near the Sun. The results show if there is no magnetic decay, then the deserved initial velocity parameters are obvious lower than the present well known results from some authors; whereas if introducing magnetic decay, for the 1‐D normal case we can make consistence among concerning results using magnetic decay time values which are supported by some authors, while for the 1‐D exponential case the results show the lackness of young pulsar samples in the larger z in the vicinity of the Sun (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Secular period decreasing of detached chromospherically active binaries

The long-term orbital period changes of a large sample of detached chromospherically active binaries (CABs) were studied. Eleven such systems were found to be undergoing secular period decreases with the rates of ?6.3×10^?9 to ?1.1×10^?6 days per year. The period decreasing rates are found to vary depending on the orbital period. The longer the orbital period is, the more rapidly the period decreases. Following Stepien (Mon. Not. R. Astron. Soc. 274:1019, 1995), the period decreasing rate predicted by angular momentum loss (AML) caused by magnetic wind is computed for each system. A comparison between the observed and calculated period decreasing rates shows that the former values are obviously larger than the latter by 1–3 orders of magnitude. It suggests that the magnetic wind is not likely the determinant mechanism driving the AML in these systems. Finally, the orbital angular momentum (AM) and the rate of AML, $\dot{J}$ , are computed for each system. It shows that the AM have a similar change with the orbital period like dP/dt does, but $\log\dot{J}/J$ presents no strict changing with the kinematical ages.