An approximate method for evaluating the shear band thickness in saturated soils

A formula for the thickness of a shear band formed in saturated soils under a simple shear or a combined stress state has been proposed. It is shown that the shear band thickness is dependent on the pore pressure properties of the material and the dilatancy rate, but is independent of the details of the combined stress state. This is in accordance with some separate experimental observations. Copyright © 2004 John Wiley & Sons, Ltd.     

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).     

Testing earthquake prediction methods: «The West Pacific short-term forecast of earthquakes with magnitude MwHRV ≥ 5.8»

Analyzing the tables and probability maps posted by Yan Y. Kagan and David D. Jackson in April 2002–September 2004 at http://scec.ess.ucla.edu/~ykagan/predictions_index.html and the catalog of earthquakes for the same period, the conclusion is drawn that the underlying method could be used for prediction of aftershocks, while it does not outscore random guessing when main shocks are considered.     

用逐步回归预测棉铃虫发生期和发生量

收集了1973年至2000年的气象，农作物、棉铃虫虫害等196年因子，以运城、汾阳和临汾为山西省代表站，利用最优二分割法把发生量分为10级，5级和不分级3种情况；使用逐步回归计算了219个模型，从中选出27个棉铃虫二、三代发生期，发生量最优模型，在农气服务中应用，预报准确率为95％，效果良好。     

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

地壳由半无限大的基岩上一层厚度为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.