Decoupled seismic analysis of an earth dam

The seismic stability of an earth dam is evaluated via the decoupled displacement analysis using the accelerograms obtained by ground response analysis to compute the earthquake-induced displacements. The response analysis of the dam is carried out under both 1D and 2D conditions, incorporating the non-linear soil behaviour through the equivalent linear method. Ten artificial and five real accelerograms were used as input motions and four different depths were assumed for the bedrock.1D and 2D response analyses were in a fair agreement with the exception of the top third of the dam where only a 2D modelling of the problem could ensure that the acceleration field is properly described. The acceleration amplification ratio obtained in the 2D analyses was equal to about 2 in all the cases considered, consistently with data from real case histories.The maximum permanent displacements computed by the sliding block analysis were small, being less than 10% of the service freeboard; a satisfactory performance of the dam can then be envisaged for any of the seismic scenarios considered in the analyses.     

System identification of instrumented bridge systems

Several recorded motions for seven bridge systems in California during recent earthquakes were analysed using parametric and non‐parametric system identification (SI) methods. The bridges were selected considering the availability of an adequate array of accelerometers and accounting for different structural systems, materials, geometry and soil types. The results of the application of SI methods included identification of modal frequencies and damping ratios. Excellent fits of the recorded motion in the time domain were obtained using parametric methods. The multi‐input/single‐output SI method was a suitable approach considering the instrumentation layout for these bridges. Use of the constructed linear filters for prediction purposes was also demonstrated for three bridge systems. Reasonable prediction results were obtained considering the various limitations of the procedure. Finally, the study was concluded by identifying the change of the modal frequencies and damping of a particular bridge system in time using recursive filters. Copyright © 2003 John Wiley & Sons, Ltd.     

Application of extended distinct element method with lattice model to collapse analysis of RC bridges

This paper proposes a simple lattice model for collapse analysis of RC bridges subjected to earthquakes by using the extended distinct element method (EDEM). In the model, a concrete element consists of lumped masses connected to one another by springs, and a reinforcement bar is represented by a discrete model or an integrated model. The proposed lattice model is simple but its parameters are reasonably defined. It has fewer element nodes and connecting springs, which will be of benefit by shortening the CPU time. The processes to determine the initial stiffness of concrete and steel springs, the parameters of the constitutive model and the fracture criteria for springs are described. A re‐contact spring model is also proposed to simulate the re‐contact of the concrete after fracture of springs; and a general grid searching method is used to decrease the CPU time for judging re‐contact after fracture. The lattice model is assessed by numerical simulations and experiments. As an application, a damaged single‐column pier subjected to the Kobe Earthquake in 1995 is analysed by EDEM with the proposed model. The simulation results indicate that the proposed model predicts well qualitatively the collapse process of RC bridges. Copyright © 2003 John Wiley & Sons, Ltd.     

System identification of linear structures based on Hilbert–Huang spectral analysis. Part 2: Complex modes

A method, based on the Hilbert–Huang spectral analysis, has been proposed by the authors to identify linear structures in which normal modes exist (i.e., real eigenvalues and eigenvectors). Frequently, all the eigenvalues and eigenvectors of linear structures are complex. In this paper, the method is extended further to identify general linear structures with complex modes using the free vibration response data polluted by noise. Measured response signals are first decomposed into modal responses using the method of Empirical Mode Decomposition with intermittency criteria. Each modal response contains the contribution of a complex conjugate pair of modes with a unique frequency and a damping ratio. Then, each modal response is decomposed in the frequency–time domain to yield instantaneous phase angle and amplitude using the Hilbert transform. Based on a single measurement of the impulse response time history at one appropriate location, the complex eigenvalues of the linear structure can be identified using a simple analysis procedure. When the response time histories are measured at all locations, the proposed methodology is capable of identifying the complex mode shapes as well as the mass, damping and stiffness matrices of the structure. The effectiveness and accuracy of the method presented are illustrated through numerical simulations. It is demonstrated that dynamic characteristics of linear structures with complex modes can be identified effectively using the proposed method. Copyright © 2003 John Wiley & Sons, Ltd.     

压裂井经济分析的优化方法

水力压裂是油气井增产的一项重要技术措施。压裂效果的经济评价，实际上是一个压裂设计的优化过程。研究经济问题即为研究压裂优化问题。提出了压裂井的施工费用计算方法和贴现率、净收入、净现值的计算公式，并对压裂井进行了经济分析，给出了压裂设计优化的步骤。建立了八屋—后五家户气田压裂经济评价模型，并应用到实际施工现场。使油气单井产量比过去增加2～3倍。     

电法在濒海及岛屿地区寻找地下淡水中的应用

采用电法的基本原理，研究各种不同岩性和成淡水的电性差畀，划分地下成淡水分界面。在严重缺水的濒海及岛屿地区找到了较好的地下淡水，解决了守岛部队吃水难的问题。     

A hybrid indirect boundary element—discrete wave number method applied to simulate the seismic response of stratified alluvial valleys

A hybrid indirect boundary element – discrete wavenumber method is presented and applied to model the ground motion on stratified alluvial valleys under incident plane SH waves from an elastic half-space. The method is based on the single-layer integral representation for diffracted waves. Refracted waves in the horizontally stratified region can be expressed as a linear superposition of solutions for a set of discrete wavenumbers. These solutions are obtained in terms of the Thomson–Haskell propagators formalism. Boundary conditions of continuity of displacements and tractions along the common boundary between the half-space and the stratified region lead to a system of equations for the sources strengths and the coefficients of the plane wave expansion. Although the regions share the boundary, the discretization schemes are different for both sides: for the exterior region, it is based on the numerical and analytical integration of exact Green's functions for displacements and tractions whereas for the layered part, a collocation approach is used. In order to validate this approach results are compared for well-known cases studied in the literature. A homogeneous trapezoidal valley and a parabolic stratified valley were studied and excellent agreement with previous computations was found. An example is given for a stratified inclusion model of an alluvial deposit with an irregular interface with the half-space. Results are displayed in both frequency and time domains. These results show the significant influence of lateral heterogeneity and the emergence of locally generated surface waves in the seismic response of alluvial valleys.     

镍铜硫化物矿石中磁黄铁矿固溶体的退火及其选矿意义

磁黄铁矿固治体从硫化物熔体结晶后，在缓慢冷却过程中经历了显著的退火。出治和出治体的租化是固治体退火的两种方式。叶片状的单斜磁黄铁矿和“火焰状”的镍黄铁矿原始出治相在降温过程中均可发生退火和租化。分布于磁黄铁矿等矿物粒间或包于磁黄铁矿粒内的粒状镍黄铁矿，不只是高温出治的直接产物，有一部分可能是由火焰状出治体租化而成的。磁黄铁矿中单斜变体的出治和租化可使矿石的磁性发生改变，镍黄铁矿出治体的租化使含镍矿物的粒度加大。因而，退火作用对矿石的选矿工艺性能有着显著影响。     

Continuous wavelet-like filter for a spherical surface and its application to localized admittance function on Mars

We have developed a 2D isotropic continuous wavelet-like transform for a spherical surface. The transform is simply defined as the surface convolution between the original field and a kernel, based on the zeroth-order Bessel function with a spherical correction. This spherical correction violates the geometric similarity for the various scales of the kernels, which becomes more apparent at longer wavelengths. We found numerically that this transform is practically equivalent to a Gaussian bandpass filter in the spherical harmonic domain. We have applied this wavelet-like transform on the recently acquired Martian gravity and topography fields. Using a ratio constructed locally from these two fields, we have constructed a map describing the lateral variations of the localized admittance function on Mars.