地震动下考虑各向异性土体-盾构隧道数值模拟

针对层状土体-盾构隧道地震反应分析,引入了横观各向同性弹塑性模型,建立了横观各向同性介质的双渐近-多向透射边界条件。针对盾构隧道抗震设计的特点,基于横观各向同性弹塑性模型,研制了考虑层状土体各向异性和施工开挖效应,适合于地下结构动力计算的弹塑性动力有限元程序。在程序中对不同的材料采用了不同的本构关系和单元形式,并采用了关联流动法则和多种屈服准则,可同时进行横观各向同性土体与地下结构的平面应力、平面应变和轴对称问题的静动力有限元分析。最后,利用所研制的程序对上海地铁二号线石门一路站附近区间隧道在不同超越概率地震波输入下的地震反应进行了计算。结果表明,在层状土体-地铁盾构隧道的抗震分析中考虑土体各向异性的影响是必要的,所提出的计算模型是合理易行的。     

Homogenization of anisotropic elastoplastic behaviors of a porous polycrystal with interface effects

This paper is devoted to develop a theoretical framework to predict the macroscopic transversely isotropic elastoplastic behavior of clay‐like material, which is viewed as a porous polycrystal. We consider evolutions of two local plastic mechanisms of grains and interface simultaneously, for which a Schmid criterion is used for the strength of sheet‐like grains and a Tresca criterion for the strength of interfaces between particles. By adapting the standard incremental method, we propose firstly a classic self‐consistent model, which does not consider the effect of interface, then a generalized self‐consistent model in which the solid phase is represented by laminated (or isotropic) spherical grains surrounded by interfaces. Comparisons of numerical predictions between these two methods are performed and have demonstrated the validity of the generalized self‐consistent model taking account of interface effects. Numerical simulations of uniaxial compression tests have shown that the macroscopic elastoplastic behavior of polycrystalline (clay‐like) material can be successfully predicted by the way of considering the two local plastic mechanisms at microscopic scale. Copyright © 2013 John Wiley & Sons, Ltd.     

Combined action of uniform flow and oscillating flow around marine riser at low Keulegan-Carpenter number

With the increase of petroleum and gas production in deep ocean, marine risers of circular cylinder shape are widely used in the offshore oil and gas platform. In order to research the hydrodynamic performance of marine risers, the dynamic mesh technique and User-Defined Function （UDF） are used to simulate the circular cylinder motion. The motion of a transversely oscillat-ing circular cylinder in combination of uniform flow and oscillating flow is simulated. The uniform flow and oscillating flow both are in x direction. SIMPLE algorithm is used to solve the Navier-Stokes equations. The User-Defined Function is used to control the cylinder transverse vibration and the inlet flow. The lift and drag coefficient changing with time and the map of vorticity isolines at different phase angle are obtained. Force time histories are shown for uniform flow at Reynolds number （Re） of 200 and for the com-bination of uniform and oscillating flows. With the increase of amplitude of oscillating flow in combined flow, the change of lift am-plitude is not sensitive to the the change of cylinder oscillating frequency. Lift amplitude increases with the increase of oscillating flow amplitude in the combined flow, but there is no definite periodicity of the lift coefficient. The drag and inertia force coefficients change when the maximum velocity of the oscillating flow increases in the combined flow. The vortex shedding near the circular cylinder shows different characteristics.     

Elastic solutions of displacements for a transversely isotropic full space with inclined planes of symmetry subjected to a point load

This work presents analytical solutions for displacements caused by three‐dimensional point loads in a transversely isotropic full space, in which transversely isotropic planes are inclined with respect to the horizontal loading surface. In the derivation, the triple Fourier transforms are employed toyield integral expressions of Green's displacement; then, the triple inverse Fourier transforms and residue calculus are performed to integrate the contours. The solutions herein indicate that the displacements are governed by (1) the rotation of the transversely isotropic planes (?), (2) the type and degree of material anisotropy (E/E′, ν/ν′, G/G′), (3) the geometric position (r, φ, ξ) and (4) the types of loading (P_x, P_y, P_z). The solutions are identical to those of Liao and Wang (Int. J. Numer. Anal. Methods Geomechanics 1998; 22 (6):425–447) if the full space is homogeneous and linearly elastic and the transversely isotropic planes are parallel to the horizontal surface. Additionally, a series of parametric study is conducted to demonstrate the presented solutions, and to elucidate the effect of the aforementioned factors on the displacements. The results demonstrate that the displacements in the infinite isotropic/transversely isotropic rocks, subjected to three‐dimensional point loads could be easily determined using the proposed solutions. Also, these solutions could realistically imitate the actual stratum of loading situations in numerous areas of engineering. Copyright © 2007 John Wiley & Sons, Ltd.     

横观各向同性饱和地基竖向振动的衰减特性

基于提出的横观各向同性饱和多孔介质Biot波动方程的一般解,研究了饱和半空间地基在竖向点源简谐激振荷载作用下地表振动的衰减特征,分析了激振频率以及横观各向同性饱和土介质的各向异性参数和孔隙渗透系数对地表振动特征的影响。计算结果表明,低频和高频激振时,地表位移衰减特性存在明显差异;在饱和土的各向异性参数中,纵向和水平方向动态渗透系数比值和刚度系数比值对地表位移衰减影响最大,这也说明采用各向同性饱和介质的动力学模型不能准确地描述具有明显各向异性特性的饱和土地基的动力特性。     

2-D与3-D横向各向同性介质中波慢度面研究

首先简述了2-D横向各向同性介质中的波慢度面方程及其似P波、似SV波与似SH波速度关系式;接着推导了3-D横向各向同性介质中波慢度面方程,导出了3-D条件下地震波速度随波动传播方向变化的关系式;计算并讨论了2-D与3-D横向各向同性介质中地震波波慢度面及地震波波速,并发现,似P波、似SV波及似SH波波慢度面一般为不规则椭圆或椭球;似P波、似SV波与似SH波速度不仅与θ、(?)有关,而且,似SV波与似SH波的传播速度不同。     

Determining elastic constants of transversely isotropic rocks using Brazilian test and iterative procedure

The elastic constants of rocks are the basic parameters for rock mechanics, and play a very important role in engineering design. There are many laboratory methods to determine the elastic constants of transversely isotropic rocks, and the Brazilian test is a popular method. This paper presented a method combination of the Brazilian test, back calculation, and iterative procedure to evaluate the five independent elastic constants of transversely isotropic rocks in laboratory. The strain data at the centre of discs were obtained using Brazilian test. The stresses at the centre of discs were computed using numerical programs. By using back calculation, the temporary elastic constants were computed after the stresses and stains were substituted into elastic mechanics equations. After iterative procedure, the convergent values of the elastic constants can be obtained. One numerical example and three experimental cases were proposed to show the applicability of this method. The convergent values of the five independent elastic constants can be obtained in no more than 10 iterative cycles. The results coming from numerical analysis method exhibited satisfactory outcome in accordance with those of generalized reduced gradient method. The merits of this method include convenient specimen preparation of the Brazilian test, simple iterative procedure, and readily available commercially numerical programs, so that this method can be easily popularized in research and engineering analysis. Copyright © 2007 John Wiley & Sons, Ltd.     

横观各向同性地基中管桩扭转振动响应解析解

考虑地基沉积过程中产生的竖向和水平向力学性质的差异,对横观各向同性地基中管桩扭转振动频域响应进行了理论研究。基于横观各向同性材料的本构关系以及桩-土耦合扭转振动,建立了桩土系统定解问题,通过Laplace变换和分离变量法求得了桩周土和桩芯土扭转振动位移形式解。通过桩-土接触面的连续条件,求得了管桩扭转频域响应解析解,并得到了桩顶复动刚度和速度导纳的表达式。将所得解退化到横观各向同性地基中实心桩解以及均匀地基中管桩解,并与已有文献进行了对比,验证了解的合理性。通过数值算例,分析了桩周土和桩芯土的横观各向同性力学参数对桩顶扭转复刚度及速度导纳的影响。     

多层横观各向同性地基轴对称固结的传递矩阵解

从横观各向同性地基轴对称Biot固结的基本方程出发,通过关于t的Laplace变换和关于r的Hankel变换,得到关于z的一阶常微分方程组。然后,对变换域内的基本未知量进行线性化处理,建立了变换域内的基本状态变量在z = 0处和任意深度处z的显式关系。利用传递矩阵法,结合层间连续性条件和边界条件,得到了多层横观各向同性地基的Biot固结轴对称问题的解答。该解答能避免随着层数增加而需要求解大型方程组的困难,明显地提高了计算效率。     

墙下条形基础与层状横观各向同性地基共同作用

运用对偶积分方程来求解层状横观各向同性地基与墙下条形基础的共同作用问题。从直角坐标平面应变问题控制方程出发,通过傅里叶（Fourier）变换和层间连续性条件,可以得到层状横观各向同性地基的传递矩阵解。基于该传递矩阵解,并利用条形基础与地基接触的混合边值条件,推导出一组关于基础挠度和地基反力的对偶积分方程。考虑墙下条形基础受到竖向集中荷载的情况,利用弹性薄板理论先求解出条形基础挠度;随后应用雅可比（Jacobi）正交多项式和级数展开的方法,将对偶积分方程转化为线性代数方程组进行求解。编制了相应的计算程序,其计算结果与有限元软件ABAQUS的结果基本吻合,从而验证了所提理论的正确性。算例分析表明,板土相对刚度与地基成层性对地基反力、地表沉降和沿z轴竖向正应力有很大的影响。