P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.     

Variational iteration solving method for El Nio phenomenon atmospheric physics of nonlinear model

1 IntroductionThe El Ni’o atmospheric physics oscillation is anabnormal phenomenon involved in the tropical Pacificocean- atmosphere interactions. Studies on the El Ni(ophenomenon are very attractive (Lin et al., 2000; Linet al., 2001, 2002; Wang, 2001; Feng et al., 2001;Feng et al., 2002; Liu et al., 2002; Wei and Chen,2003; Xie et al., 2002; Zhu et al., 2002; Pu et al.,2003; Gu et al., 2004; Yu and Liu, 2004). Lin and Mo(2004), Mo and Lin (2004) and Mo et al. (2004) alsostudied a s…     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

Semi-analytical solutions of a contaminant transport equation with nonlinear sorption in 1D

A new method to determine semi-analytical solutions of one-dimensional contaminant transport problem with nonlinear sorption is described. This method is based on operator splitting approach where the convective transport is solved exactly and the diffusive transport by finite volume method. The exact solutions for all sorption isotherms of Freundlich and Langmuir type are presented for the case of piecewise constant initial profile and zero diffusion. Very precise numerical results for transport with small diffusion can be obtained even for larger time steps (e.g., when the Courant-Friedrichs-Lewy (CFL) condition failed).     

有限弹性变形地壳中的地震波

地壳由半无限大的基岩上一层厚度为H^-的表土层组成，入射地震波为垂直的SH波，产生水平地面运动。当浅源大地震发生时，在极震区以外行波传播产生地面运动将使地壳介质有非线性的有限弹性变形。用小参数摄动法使非线性控制方程为线性化的小参数各阶控制方程，得出头两阶线性控制方程的解析解。     

Development of site-specific design ground motions in Western and Eastern North America

This paper presents results recently obtained for generating site-specific ground motions needed for design of critical facilities. The general approach followed in developing these ground motions using either deterministic or probabilistic criteria is specification of motions for rock outcrop or very firm soil conditions followed by adjustments for site-specific conditions. Central issues in this process include development of appropriate attenuation relations and their uncertainties, differences in expected motions between Western and Eastern North America, and incorporation of site-specific adjustments that maintain the same hazard level as the control motions, while incorporating uncertainties in local dynamic material properties. For tectonically active regions, such as the Western United States (WUS), sufficient strong motion data exist to constrain empirical attenuation relations for M up to about 7 and for distances greater than about 10–15 km. Motions for larger magnitudes and closer distances are largely driven by extrapolations of empirical relations and uncertainties need to be substantially increased for these cases.

For the Eastern United States (CEUS), due to the paucity of strong motion data for cratonic regions worldwide, estimation of strong ground motions for engineering design is based entirely on calibrated models. The models are usually calibrated and validated in the WUS where sufficient strong motion data are available and then recalibrated for applications to the CEUS. Recalibration generally entails revising parameters based on available CEUS ground motion data as well as indirect inferences through intensity observations. Known differences in model parameters such as crustal structure between WUS and CEUS are generally accommodated as well. These procedures are examined and discussed.     

河南及邻区地震前中小地震活跃-平静的定量指标研究

以累计频度定量计算方法，分区讨论了河南及邻区1970年以来地震活动非线性度ZL值的时间进程曲线，系统计算了M≥5．0级地震前的ZL值，结果表明：开始出现活跃(平静)异常的时间集中在地震前1—2年，结束异常的时间在地震前一年内，集中于0—4个月。同时提出了异常的定量指标以及发震类型。     

Nonlinear Rate Dependent Model of High Damping Rubber Bearing

The accurate analysis of the response of isolated structures requires the application of appropriate models of isolation devices. The purpose of this paper is to analyse a nonlinear strain rate dependent model of a high damping rubber bearing which simulates the horizontal behaviour of the device under specified vertical load using a nonlinear elastic spring-dashpot element. The effectiveness of the model is checked by fitting the experimental data concerning three different rubber bearings. The results of the study show that the model can simulate the bearing behaviour over a wide shear strain range with small simulation errors. This revised version was published online in July 2006 with corrections to the Cover Date.     

小波网络方法用于时间序列分析中非线性模型的研究

把小波网络方法用于时间序列分析中非线性模型，利用小波网络的非线性逼近能力，构造一个非线性函数，并对温度资料进行处理，得到了较好的拟合效果。