An improved numerical simulation mode of nonlinear wave with consideration of high order terms

In this paper the 0-1 combined BEM is adopted to subdivide the computational domain boundary,and to discretize the Green’s integral expression based on Laplace equation.The FEM is used to subdivide the wave surface and deduce the surface equation which satisfies the nonlinear boundary conditions on the surface.The equations with potential function and wave surface height as an unknown quantity by application of Taylor expansion approach can be solved by iteration within the time step.In m-time iteration within the computational process of time step(n-1)Δt to nΔt,the results of the previous iteration are taken as the initial value of the two-order unknown terms in the present iteration.Thus,an improved tracking mode of nonlinear wave surface is established,and numerical results of wave tank test indicate that this mode is improved obviously and is more precise than the previous numerical model which ignored the two-order unknown terms of wave surface location and velocity potential function in comparison with the theoretical values.     

Seismic soil–pile–structure kinematic and inertial interaction—A new beam approach

The main purpose of this study is to investigate the accuracy of an advanced beam model for the soil–pile–structure kinematic and inertial interaction and demonstrate its efficiency and advantages compared to other commonly used beam or solid models. Within this context, a Beam on Nonlinear Winkler Foundation model is adopted based on the Boundary Element Method (BEM), accounting for the effects induced by geometrical nonlinearity, rotary inertia and shear deformation, employing the concept of shear deformation coefficients. The soil nonlinearity is taken into consideration by means of a hybrid spring configuration consisting of a nonlinear (p–y) spring connected in series to an elastic spring–damper model. The nonlinear spring captures the near-field plastification of the soil while the spring–damper system (Kelvin–Voigt element) represents the far-field viscoelastic character of the soil. An extensive case study is carried out on a pile-column–deck system of a bridge, found in two cohesive layers of sharply different stiffness and subjected to various earthquake excitations, providing insight to several phenomena. The results of the proposed model are compared with those obtained from a Beam-FE solution as well as from a rigorous fully three-dimensional (3-D) continuum FE scheme.     

二维地基波阻板隔振分析

本文采用了以薄层法层状半空间基本解答作为格林函数的边界元法,对均质弹性半空间和层状半空间中二维波阻板隔振设计进行了详细的参数分析。分析结果表明,均质弹性半空间和层状半空间内采用波阻板隔振均可取得较好的隔振效果;应保证波阻板具有合理的宽度和较小的埋深,才能得到理想的隔振效果;同时增加波阻板的厚度和模量是提高隔振效果的最有效的2种措施;此外,分层性对波阻板隔振效果也具有一定的影响。     

Modeling pressure pipe embedded in two-layer soil by a half-plane BEM

In this paper, an elastostatic half-plane boundary element method (BEM) formulation was applied to analyze the stress behavior of underground pressure pipes, embedded in two-layer soils. In the use of this method, only the boundary of pipe and interfaces were required to be discretized. In this regard, first, a computer code was prepared based on a multi-region substructuring process in the BEM scheme. Then, the efficiency and applicability of the method as well as the prepared algorithm were verified by solving some practical examples and comparing the results with those of the published works. Finally, a parametric study was done to evaluate the effect of pipe depth and determine the soil stress distribution. The studies showed that the half-plane BEM was in good agreement with the existing solutions and its capability was very favorable for elastostatic problems including semi-infinite domain. It is obvious that this method can be practically used to analyze the geotechnical underground buildings in substituting the full-plane BEM formulation.     

Soil-structure interaction using BEM–FEM coupling through ANSYS software package

The main aim of this work is to develop, verify and apply in simulation study an efficient hybrid approach to study seismic response of a soil-structure system taking into account all the important components as: (1) the line time-harmonic source with its specific geophysical properties; (2) the inhomogeneity and heterogeneity of the wave path from the source to the local geological region; (3) the geotechnical properties of the near-field local geological profile and finally (4) the properties of the engineering structure itself. Plane strain state is considered. The hybrid computational tool is based on the boundary element method (BEM¹) for modeling the infinite far-field geological media and finite element method (FEM²) for treating the dynamic behavior of the structure and the near-field finite soil geological region. Each of the two techniques is applied in that part of the whole model where it works more efficiently. The hybrid numerical scheme is realized via the sub-structure approach, direct BEM¹, conventional FEM² and insertion of the BEM¹ model of the seismically active far-field geological media as a macro-finite element (MFE³) in the FEM² commercial program ANSYS. The accuracy and verification study of the proposed method is presented by solution of numerical test examples simulating different seismic scenarios. The obtained results show clearly that the hybrid model is able to demonstrate the sensitivity of the synthetic signals to the source properties, to the heterogeneous character of the wave path, to the relief peculiarities of the local layered geological deposit and to the specific properties of the engineering structure.     

饱和土半解析边界元法及在双排桩被动隔振中的应用

针对均质饱和地基中双排桩远场被动隔振问题,首先,基于饱和多孔介质的边界元法,建立了以薄层法（TLM）基本解作为动力Green函数的饱和土半解析边界元法,该方法可有效地分析饱和半空间的土-结构动力相互作用问题。在此基础上,根据双排桩与周围土体的邻接条件,推导了双排桩对Rayleigh波散射的三维边界元方程;运用上述边界元方程,对双排桩的远场被动隔振问题进行了研究。结果表明：双排桩能够有效地降低屏障后的位移振幅,其隔振效果要优于单排桩;排间净距对隔振效果影响不大,而桩间净距则对双排桩隔振效果起控制作用。     

A fast multipole boundary element method for three-dimensional potential flow problems

A fast multipole methodology (FMM) is developed as a numerical approach to reduce the computational cost and memory requirements in solving large-scale problems. It is applied to the boundary element method (BEM) for three-dimensional potential flow problems. The algorithm based on mixed multipole expansion and numeric, al integration is implemented in combination with an iterative solver. Numerical examinations, on Dirichlet and Neumann problems, are carried out to demonstrate the capability and accuracy of the present method. It has been shown that the method has evident advantages in saving memory and computing time when used to solve huge-scale problems which may be prohibitive for the traditional BEM implementation.     

饱和地基中空沟近场隔振的现场试验与二维半解析边界元分析

动力机器运行和车辆行驶等会产生振动污染,危及邻近建筑物安全和干扰精密仪器设备正常运行等。这些振动污染可通过在地基中设置空沟的方式来降低或消除。针对饱和地基上明置动力机器基础的环境振动影响及空沟近场隔振问题,进行了饱和地基上空沟近场隔振的现场试验,并对试验结果进行了无量纲化分析;基于饱和土半解析边界元法,分别推导了动力机器基础环境振动影响和空沟近场隔振的边界元方程;在此基础上,详细研究了空沟对动力机器基础振动影响的隔振效果,分析了空沟深度、宽度和距振源距离对其隔振效果的影响。结果表明:空沟能够有效的降低动力机器基础的环境振动影响;空沟宽度对其隔振效果影响相对较小,而空沟深度对其隔振效果影响较大,为获得较好的隔振效果,空沟深度建议取1倍Rayleigh波波长;空沟距振源距离对其隔振效果也有较大影响,距离越远则隔振效果也越好,当被保护建筑距振源较远时,建议空沟在被保护建筑附近设置。此外,在某些特殊情况下,空沟隔振系统会由于共振现象而出现隔振效果劣化的现象,在工程设计中应予以注意。     

Coupling of BEM with a large displacement and rotation algorithm

In stability analysis of rock blocks, the deformability of the blocks can conveniently be simulated using the boundary element method (BEM). However, all boundary conditions are given as stresses. Thus, the displacement solution is not unique. In this paper, an algorithm is proposed to remove rigid body motions in the solution of the boundary form of Somigliana identity discretized by the direct BEM formulation. The algorithm is applied to the calculation of the normal stiffness of rock blocks and coupled with BS3D, large displacement and rotation algorithm for the general stability of rock blocks. Copyright © 2010 John Wiley & Sons, Ltd.     

Deformations resulting from the movements of a shear or tensile fault in an anisotropic half space

Earlier solutions (Bull. Seismol. Soc. Amer. 1985; 75 :1135–1154; Bull. Seismol. Soc. Amer. 1992; 82 :1018–1040) of deformations caused by the movements of a shear or tensile fault in an isotropic half-space for finite rectangular sources of strain nucleus have been extended for a transversely isotropic half-space. Results of integrating previous solutions (Int. J. Numer. Anal. Meth. Geomech. 2001; 25 (10): 1175–1193) of deformations due to a shear or tensile fault in a transversely isotropic half-space for point sources of strain nucleus over the fault plane are presented. In addition, a boundary element (BEM) model (POLY3D:A three-dimensional, polygonal element, displacement discontinuity boundary element computer program with applications to fractures, faults, and cavities in the Earth's crust. M.S. Thesis, Stanford University, Department of Geology, 1993; 62) is given. Different from similar researches (e.g. Thomas), the Akaike's view on Bayesian statistics (Akaike Information Criterion Statistics. D. Reidel Publication: Dordrecht, 1986) is applied for inverting deformations due to a fault to obtain displacement discontinuities on the fault plane. An example is given for checking displacements predicted by proposed analytical expressions. Another example is generated for the use of proposed BEM model. It demonstrates the effectiveness of this model in exploring displacement behaviours of a fault. Copyright © 2004 John Wiley & Sons, Ltd.