On engineering equivalent linear analysis of soils in the frequency domain

This study investigated the performance, in the frequency domain, of an engineering equivalent linear method that is widely used in soil dynamics. A stochastic linear method using an extended differential equation by Wen and Baber was employed as the base line for evaluation. Stationary responses of single-degree-of-freedom non-linear systems under filtered white noise from both approaches were compared. For both a smooth nearly elasto-plastic material and soils, the study found that the engineering equivalent linear analysis underestimated rms displacements, but provided good estimates of rms acceleration responses. Furthermore, the frequency response functions of acceleration from both approaches were close at all levels of response.     

A hybrid analysis method for linear dynamic soil-structure interaction in time and frequency domain

IntroductionThe analysis of dynamic soil-structure interaction for important engineering project is still based on linear model (including equivalent linear model) with complex damping, and traditional frequency domain method (Lysmer, et al, 1975, 1981; DING, et al, 1999). Namely, first calculating frequency domain solution by Fourier transform, and then calculating time domain solution by Fourier inverse transform. The motion equation of a system in frequency domain is usually written as (…     

Dynamic response analysis in the frequency domain

Analysis in the frequency domain using discrete Fourier transforms is an efficient means of calculating the dynamic response of linear systems. In fact, for systems with frequency dependent parameters and also in those cases where the complex frequency response functions are more easily determined, frequency domain analysis may be the only effective means of determining the dynamic response. The use of discrete transforms along with finite summation requires that the forcing function and the impulse function be converted into periodic forms. This modification may introduce unacceptable errors in the results of analysis, unless appropriate steps are taken to avoid or minimize the effect of aliasing or overlapping. For single-degree-of-freedom systems, procedures that will eliminate the effect of aliasing have already been developed. However, problems related to frequency domain analysis still exist for multi-degree-of-freedom systems with non-proportional damping, in analysis through substructuring and in those cases where a continuum solution is involved. A new procedure which addresses these problems and is applicable to both single- and multi-degree-of-freedom system as well as to analysis through substructuring is presented here.     

An incremental mode-superposition for non-linear dynamic analysis

This paper uses an incremental mode-superposition procedure to compute the inelastic dynamic response of multi-degree-of-freedom systems. A damping matrix proportional to the instantaneous properties is used throughout the analysis. The non-linear response of several shear type plane and space frames with elastic-plastic and bilinear column properties subjected to ground excitation was computed by both the incremental mode-superposition and the direct integration of the coupled equations of motion. When all modes are considered, the responses computed by the incremental mode-superposition are identical to those from the direct integration. Fewer modes can also be used to compute the response with reasonable accuracy by performing the modal truncation for each time increment. The study shows that incorporating instantaneous damping in non-linear dynamic analysis is relatively simple and requires less computational time than the direct integration.     

A hybrid time frequency domain formulation for non-linear soil–structure interaction

Methods that combine frequency and time domain techniques offer an attractive alternative for solving Soil–Structure-interaction problems where the structure exhibits non-linear behaviour. In the hybrid-frequency-time-domain procedure a reference linear system is solved in the frequency domain and the difference between the actual restoring forces and those in the linear model are treated as pseudo-forces. In the solution scheme explored in this paper, designated as the hybrid-time-frequency-domain (HTFD) procedure, the equations of motion are solved in the time domain with due consideration for non-linearities and with the unbounded medium represented by frequency-independent springs and dampers. The frequency dependency of the impedance coefficients is introduced by means of pseudo-forces evaluated in the frequency domain at the end of each iteration. A criterion of stability for the HTFD approach is derived analytically and its validity is sustained numerically. As is often the case, the criterion takes the form of a limit of unity on the spectral radius of an appropriately defined matrix. Inspection of the terms in this matrix shows that convergence can be guaranteed by suitable selection of the reference impedance. The CPU times required to obtain converged solutions with the HTFD are found, in a number of numerical simulations, to be up to one order of magnitude less than those required by the alternative hybrid-frequency-time-domain approach. © 1998 John Wiley & Sons, Ltd.     

频率域波形反演中与频率相关的影响因素分析

波动方程深度偏移是解决复杂地质体成像的关键技术,基于波动方程的速度建模为其提供更为精确的速度模型.频率域波形反演是目前研究最为广泛的波动方程速度建模方法之一,它推动了波形反演在勘探尺度下的应用.本文通过对频率域波形反演的实现,分析对比了其有效执行过程中与频率相关的影响因素.介绍了时间域的多尺度反演方法在频率域的一种实现方式,对比分析了输入数据的频点带宽和应用的子波频带范围不同时对反演结果的影响.本文通过设计的山地地质模型对频率域波形反演进行了测试和对比,得到的结论为频率域波形反演的有效计算提供了依据和参考.     

An imposed force summation method for non-linear dynamic analysis

A method of obtaining the response of a non-linear system is proposed, which involves imposing the solution upon the equations of motion for a linear system and solving using selected modes of vibration. Contribution from the remaining modes of the system is included in terms of their quasi-static response. The method is developed and illustrated with reference to a system with material non-linearity and special reference is given to the dynamic analysis of embankment dams.     

An operational modal analysis method in frequency and spatial domain

A frequency and spatial domain decomposition method （FSDD） for operational modal analysis （OMA） is presented in this paper, which is an extension of the complex mode indicator function （CMIF） method for experimental modal analysis （EMA）. The theoretical background of the FSDD method is clarified, Singular value decomposition is adopted to separate the signal space from the noise space. Finally, an enhanced power spectrum density （PSD） is proposed to obtain more accurate modal parameters by curve fitting in the frequency domain. Moreover, a simulation case and an application case are used to validate this method.     

Viscous damping formulation and high frequency motion propagation in non-linear site response analysis

Non-linear time domain site response analysis is widely used in evaluating local soil effects on propagated ground motion. This approach has generally provided good estimates of field behavior at longer periods but has shortcomings at relatively shorter periods. Viscous damping is commonly employed in the equation of motion to capture damping at very small strains and employs an approximation of Rayleigh damping using the first natural mode only. This paper introduces a new formulation for the viscous damping using the full Rayleigh damping. The new formulation represents more accurately wave propagation for soil columns greater than 50 m thick and improves non-linear site response analysis at shorter periods. The proposed formulation allows the use of frequency dependent viscous damping. Several examples, including a field case history at Treasure Island, California, demonstrate the significant improvement in computed surface response using the new formulation.     

Three-dimensional dynamic soil–structure interaction analysis in the time domain

A new numerical procedure is proposed for the analysis of three-dimensional dynamic soil–structure interaction in the time domain. In this study, the soil is modelled as a linear elastic solid, however, the methods developed can be adapted to include the effects of soil non-linearities and hysteretic damping in the soil. A substructure method, in which the unbounded soil is modelled by the scaled boundary finite-element method, is used and the structure is modelled by 8–21 variable-number-node three-dimensional isoparametric or subparametric hexahedral curvilinear elements. Approximations in both time and space, which lead to efficient schemes for calculation of the acceleration unit-impulse response matrix, are proposed for the scaled boundary finite-element method resulting in significant reduction in computational effort with little loss of accuracy. The approximations also lead to a very efficient scheme for evaluation of convolution integrals in the calculation of soil–structure interaction forces. The approximations proposed in this paper are also applicable to the boundary element method. These approximations result in an improvement over current methods. A three-dimensional Dynamic Soil–Structure Interaction Analysis program (DSSIA-3D) is developed, and seismic excitations (S-waves, P-waves, and surface waves) and externally applied transient loadings can be considered in analysis. The computer program developed can be used in the analysis of three-dimensional dynamic soil–structure interaction as well as in the analysis of wave scattering and diffraction by three-dimensional surface irregularities. The scattering and diffraction of seismic waves (P-, S-, and Rayleigh waves) by various three-dimensional surface irregularities are studied in detail, and the numerical results obtained are in good agreement with those given by other authors. Numerical studies show that the new procedure is suitable and very efficient for problems which involve low frequencies of interest for earthquake engineering. Copyright © 1999 John Wiley & Sons Ltd     

地震计阻尼和自振频率的频域测定

利用地震观测系统中的阶跃标定响应记录 ,在频域基于阶跃信号的幅频特性确定地震计的阻尼常数和自振周期。介绍了该方法的原理 ,采用非线性最小二乘拟合方法和具体的实际应用     

Time domain dynamic analysis of piles under impact loading

In this paper, time domain dynamic analysis of piles under impact loading is presented. For this purpose a hybrid boundary element technique is implemented. Linear beam column finite elements are used to model the piles and resulting governing equations are solved using an implicit integration scheme. The continuum is assumed to be elastic and an efficient step-by-step time integration scheme is implemented by using an approximate half space integral formulation. By enforcing displacement equilibrium conditions at each time step, a system of equations is generated which yields the solution. Results of this time domain formulation under linear material behavior are compared with Laplace domain results to validate the methods.     

Seismic displacement demand prediction in non-linear domain: Optimization of the N2 method

In Europe, computation of displacement demand for seismic assessment of existing buildings is essentially based on a simplified formulation of the N2 method as prescribed by Eurocode 8(EC8). However, a lack of accuracy of the N2 method in certain conditions has been pointed out by several studies. This paper addresses the assessment of effectiveness of the N2 method in seismic displacement demand determination in non-linear domain. The objective of this work is to investigate the accuracy of the N2 method through comparison with displacement demands computed using non-linear timehistory analysis(NLTHA). Results show that the original N2 method may lead to overestimation or underestimation of displacement demand predictions. This may affect results of mechanical model-based assessment of seismic vulnerability at an urban scale. Hence, the second part of this paper addresses an improvement of the N2 method formula by empirical evaluation of NLTHA results based on EC8 ground-classes. This task is formulated as a mathematical programming problem in which coefficients are obtained by minimizing the overall discrepancy between NLTHA and modified formula results. Various settings of the mathematical programming problem have been solved using a global optimization metaheuristic. An extensive comparison between the original N2 method formulation and optimized formulae highlights benefits of the strategy.     

Sonification of geophysical data through time–frequency analysis: theory and applications

In this paper, we introduce a new method of geophysical data interpretation based on simultaneous analysis of images and sounds. The final objective is to expand the interpretation workflow through multimodal (visual–audio) perception of the same information. We show how seismic data can be effectively converted into standard formats commonly used in digital music. This conversion of geophysical data into the musical domain can be done by applying appropriate time–frequency transforms. Using real data, we demonstrate that the Stockwell transform provides a very accurate and reliable conversion. Once converted into musical files, geophysical datasets can be played and interpreted by using modern computer music tools, such as sequencers. This approach is complementary and not substitutive of interpretation methods based on imaging. It can be applied not only to seismic data but also to well logs and any type of geophysical time/depth series. To show the practical implications of our integrated visual–audio method of interpretation, we discuss an application to a real seismic dataset in correspondence of an important hydrocarbon discovery.     

Criterion of stability and implementation issues of hybrid frequency-time-domain procedure for non-linear dynamic analysis

A criterion of stability pertaining to the hybrid frequency-time-domain procedure is derived. It stems from the initial-value theorem on one side and from the condition of convergence which applies to the harmonic counterpart of the procedure on the other side. The criterion depends on the dynamic stiffnesses of both the original non-linear system and the pseudo-linear system evaluated at the Nyquist frequency. Since at this frequency both stiffnesses are usually controlle by the mass, the criterion of stability is satisfied for most structural systems. The procedure must, however, be applie sequentially to individual time segments. The validity of this criterion is confirmed by a numerical investigation performe for a SDOF system. The solution of non-linear dynamic problems requires consideration of implementation issues which are also discusse in the paper. The analysis of a non-linear soil-structure-interaction system in which the soil's stiffness coefficients an directly defined in the frequency domain also demonstrates the accuracy of the hybrid frequency-time-domain procedure.     

System parameters of soil foundations for time domain dynamic analysis

Soil-structure interaction analysis is usually carried out in the frequency domain, because the compliance functions of the half-space are known only in the frequency domain. Since non-linear analysis cannot be carried out in the frequency domain, a system with frequency independent parameters is used to represent the half-space soil medium so that a nonlinear analysis in the time domain becomes possible. The objective of this paper is to propose a system with lumped parameters, which are independent of frequency, to represent the half-space soil medium. The proposed frequency independent system consists of a number of real discrete structure elements; thus the existing dynamic analysis programs may be adoptable with little modification. In this paper, the parameters are found by minimizing the sum of the squares of deviations between the steady-state responses of the theoretical half-space model and those of the lumped parameter system over a specified frequency range. Once the parameters have been found, the lumped parameter system can be used in practical applications for time domain dynamic analysis of either linear or non-linear structures. In comparison with the dynamic response of the theoretical half-space model, the lumped parameter system yields satisfactory results.     

Storm‐coverage effect on dynamic flood‐frequency analysis: empirical data analysis

In this study, a dynamic flood‐frequency analysis model considering the storm coverage effect is proposed and applied to six sub‐basins in the Pyungchang River basin, Korea. The model proposed is composed of the rectangular pulse Poisson process model for rainfall, the Soil Conservation Service curve number method for infiltration and the geomorphoclimatic instantaneous unit hydrograph for runoff estimation. Also, the model developed by Marco and Valdes is adopted for quantifying the storm‐coverage characteristics. By comparing the results from the same model with and without the storm‐coverage effect consideration, we could quantify the storm‐coverage effect on the flood‐frequency analysis. As a result of that, we found the storm‐coverage effect was so significant that overestimation of the design flood was unavoidable without its consideration. This also becomes more serious for larger basins where the probability of complete storm coverage is quite low. However, for smaller basins, the limited number of rain gauges is found to hamper the proper quantification of the storm‐coverage characteristics. Provided with a relationship curve between the basin size and the storm coverage (as in this study), this problem could be overcome with an acceptable accuracy level. Copyright © 2003 John Wiley & Sons, Ltd.     

Non-linear soil-structure-interaction analysis using dynamic stiffness or flexibility of soil in the time domain

The basic equation of motion to analyse the interaction of a non-linear structure and an irregular soil with the linear unbounded soil is formulated in the time domain. The contribution of the unbounded soil involves convolution integrals of the dynamic-stiffness coefficients in the time domain and the corresponding motions. Alternatively, a flexibility formulation for the contribution of the unbounded soil using the dynamic-flexibility coefficients in the time domain, together with the direct-stiffness method for the structure and the irregular soil can be applied. The dynamic-stiffness or flexibility coefficient in the time domain is calculated as the inverse Fourier transform of the corresponding value in the frequency domain. The dynamic-stiffness coefficient's asymptotic behaviour for high frequencies determines the singular part whose transformation exists only in the sense of a distribution. As the dynamic-flexibility coefficient converges to zero for the frequency approaching infinity, the corresponding coefficient in the time domain is simpler to calculate, as no singular part exists. The salient features of the dynamic-stiffness and flexibility coefficients in the time domain are illustrated using a semi-infinite rod with exponentially increasing area. The dynamic-flexibility coefficients in the time domain are calculated for a rigid circular disc resting on the surface of an elastic halfspace and of a layer built-in at its base. Material damping is also introduced using the three-parameter Kelvin and the Voigt models.     

Comparison between non-linear dynamic and static seismic analysis of structures according to European and US provisions

Several procedures for non-linear static and dynamic analysis of structures have been developed in recent years. This paper discusses those procedures that have been implemented into the latest European and US seismic provisions: non-linear dynamic time-history analysis; N2 non-linear static method (Eurocode 8); non-linear static procedure NSP (FEMA 356) and improved capacity spectrum method CSM (FEMA 440). The presented methods differ in respect to accuracy, simplicity, transparency and clarity of theoretical background. Non-linear static procedures were developed with the aim of overcoming the insufficiency and limitations of linear methods, whilst at the same time maintaining a relatively simple application. All procedures incorporate performance-based concepts paying more attention to damage control. Application of the presented procedures is illustrated by means of an example of an eight-storey reinforced concrete frame building. The results obtained by non-linear dynamic time-history analysis and non-linear static procedures are compared. It is concluded that these non-linear static procedures are sustainable for application. Additionally, this paper discusses a recommendation in the Eurocode 8/1 that the capacity curve should be determined by pushover analysis for values of the control displacement ranging between zero and 150% of the target displacement. Maximum top displacement of the analyzed structure obtained by using dynamic method with real time-history records corresponds to 145% of the target displacement obtained using the non-linear static N2 procedure.