Transform methods for frequency‐dependent complex stiffness to time domain using real or imaginary data only

The author has studied the transform method of complex stiffness, which is strongly dependent on frequency, to the time domain. In this paper, new transform methods that use only the real part or the imaginary part data are proposed. By applying them to example problems, it is confirmed that the accuracy of the proposed methods is almost high in both causal and noncausal cases. These methods can also be thought as a standby for the Hilbert transform, because the undefined part is easily calculated by the recovering complex function from the obtained impulse response. Copyright © 2007 John Wiley & Sons, Ltd.     

A practical method to transform frequency dependent impedance to time domain

In order to perform time history earthquake response analyses with consideration to both the dynamic soil–structure interaction and the non‐linear behaviour of the structure, it is important to transform the soil impedance in the frequency domain to the impulse response in the time domain. In this paper, a new transform method with high practicality is proposed. First, the formulation of the proposed transform method is described. Next, the validity of the method is examined using an example problem whose impulse response is analytically obtained. Then, the impedance of the rigid foundation on 2‐layered soil is transformed to the time domain, and the characteristics of the impulse response are investigated. Finally, time history earthquake response analyses of a structure on the soil using the obtained impulse response are carried out. The validity and the efficiency of the proposed method are confirmed through these investigations. Copyright © 2005 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.     

Practical causal hysteretic damping

A number of experiments indicate that the internal damping corresponding to the energy dissipation of many materials is essentially frequency independent. Accordingly, an analysis model that can express such characteristics (called a hysteretic damping model) in the time domain is needed. Although a great number of investigations into this subject have been carried out, there are a few practical methods. In this paper, a simple hysteretic damping model which satisfies the causality condition is presented using an extension of the complex stiffness transfer method that the author has proposed. Compared with the energy proportional damping model and the Biot model, the applicability and the efficiency of this model to time history response analyses were confirmed well by example problems. Copyright © 2006 John Wiley & Sons, Ltd.     

A time-domain viscous damping model based on frequency-dependent damping ratios

This paper introduces the mathematics and procedures used in developing a time-dependent damping model for integration analyses of structural response. To establish the time-dependent viscous damping model, frequency-dependent damping ratios of the structure under a series of steady-state unit impulses corresponding to actual loads are first calculated. For simplicity, the ratios can be incorporated with the static stiffness of the structure to model approximately the impulse induced damping spectrum. According to the nature of the problem, these ratios can be calculated from the theoretical impedance functions and experimental observations. With the computed damping spectrum, the damping coefficient in the time domain can be obtained with the Fourier transform technique. Adopting the impulse–response method, the damping can be modeled rationally through integration with changing loads. Numerical examples are presented to show the feasibility of this model while the transform criterion is satisfied.     

Time‐domain response of linear hysteretic systems to deterministic and random excitations

The causal and physically realizable Biot hysteretic model proves to be the simplest linear model able to describe the nearly rate‐independent behaviour of engineering materials. In this paper, the performance of the Biot hysteretic model is analysed and compared with those of the ideal and causal hysteretic models. The Laguerre polynomial approximation (LPA) method, recently proposed for the time‐domain analysis of linear viscoelastic systems, is then summarized and applied to the prediction of the dynamic response of linear hysteretic systems to deterministic and random excitations. The parameters of the LPA model generally need to be computed through numerical integrals; however, when this model is used to approximate the Biot hysteretic model, closed‐form expressions can be found. Effective step‐by‐step procedures are also provided in the paper, which prove to be accurate also for high levels of damping. Finally, the method is applied to the dynamic analysis of a highway embankment excited by deterministic and random ground motions. The results show that in some cases the inaccuracy associated with the use of an equivalent viscous damping model is too large. Copyright © 2005 John Wiley & Sons, Ltd.     

Time‐domain viscoelastic analysis of earth structures

In this paper two causal models that approximate the nearly frequency‐independent cyclic behaviour of soils are analysed in detail. The study was motivated by the need to conduct time‐domain viscoelastic analysis on soil structures without adopting the ad hoc assumption of Rayleigh damping. First, the causal hysteretic model is introduced in which its imaginary part is frequency independent the same way that is the imaginary part of the popular non‐causal constant hysteretic model. The adoption of an imaginary part that is frequency independent even at the zero‐frequency limit, in conjunction with the condition that the proposed model should be causal, yields a real part that is frequency dependent and singular at zero frequency. The paper shows that the causal hysteretic model, although pathological at the static limit, is the mathematical connection between the non‐causal constant hysteretic model and the physically realizable Biot model. The mathematical structure of the two causal models is examined and it is shown that the causal hysteretic model is precisely the high‐frequency limit of the Biot model. Although both models have a closed‐form time‐domain representation, only the Biot model is suitable for a time‐domain viscoelastic analysis with commercially available computer software. The paper demonstrates that the simplest, causal and physically realizable linear hysteretic model that can approximate the cyclic behaviour of soil is the Biot model. The proposed study elucidates how the dynamic analysis of soil structures can be conducted rigorously in terms of the viscoelastic properties of the soil material and not with the ad hoc Rayleigh damping approach which occasionally has been criticized that tends to overdamp the higher vibration modes. The study concludes that under pulse‐type motions the Rayleigh damping approximation tends to overestimate displacements because of the inappropriate viscous type of dissipation that is imposed. Under longer motions that induce several cycles, the concept of equivalent viscous damping is more appropriate and the Rayleigh damping approximation results to a response that is comparable to the response computed with a rigorous time‐domain viscoelastic finite element analysis. Copyright © 2000 John Wiley & Sons, Ltd.     

A practical method for estimating dynamic soil stiffness on surface of multi‐layered soil

It is important to estimate the influence of layered soil in soil–structure interaction analyses. Although a great number of investigations have been carried out on this subject, there are very few practical methods that do not require complex calculations. In this paper, a simple and practical method for estimating the horizontal dynamic stiffness of a rigid foundation on the surface of multi‐layered soil is proposed. In this method, waves propagating in the soil are traced using the conception of the cone model, and the impulse response function can be calculated directly and easily in the time domain with a good degree of accuracy. The characteristics of the impedance, that is the transformed value to the frequency domain of the obtained impulse response, are studied using two‐ to four‐layered soil models. The cause of the fluctuation of impedance is expressed clearly from its relation to reflected waves from the lower layer boundary in the model. Copyright © 2005 John Wiley & Sons, Ltd.     

剪力墙的等效周期和等效阻尼比研究

结构弹塑性反应可用等效线性化方法中的等效刚度和等效阻尼比进行计算。为研究钢筋混凝土剪力墙的等效刚度和等效阻尼比,通过分析剪力墙的力-位移简化四折线滞回模型,推导了其屈服点和峰值点的等效刚度、等效周期和等效阻尼比的计算公式。为验证其准确性和适用性,采用提出的峰值点等效刚度、等效阻尼比公式和FEMA 440建议的公式,分别计算了42片钢筋混凝土剪力墙试件的等效刚度和等效阻尼比,并进行了对比分析。结果表明:提出的计算公式和FEMA 440建议的公式所得到的等效周期和等效阻尼比值与试验值比较吻合,因此,所提出的计算剪力墙的等效周期和等效阻尼比的公式较为合理。     

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 (…     

Stiffness‐damping simultaneous identification using limited earthquake records

A new method of stiffness‐damping simultaneous identification of building structures is proposed using limited earthquake records. It is shown that when horizontal accelerations are recorded at the floors just above and below a specific storey in a shear building model, the storey stiffness and the damping ratio can be identified uniquely. The viscous damping coefficient and the linear hysteretic damping ratio can also be identified simultaneously in a numerical model structure. The accuracy of the present identification method is investigated through the actual limited earthquake records in a base‐isolated building. It is further shown that an advanced identification technique for mechanical properties of a Maxwell‐type model can be developed by combining the present method with a perturbation technique. Copyright © 2000 John Wiley & Sons, Ltd.     

Simplified seismic design approach for steel portal frame piers with hysteretic dampers

A simplified seismic design procedure for steel portal frame piers installed with hysteretic dampers is proposed, which falls into the scope of performance‐based design philosophy. The fundamental goal of this approach is to design a suite of hysteretic damping devices for existing and new bridge piers, which will assure a pre‐defined target performance against future severe earthquakes. The proposed procedure is applicable to multi‐degree‐of‐freedom systems, utilizing an equivalent single‐degree‐of‐freedom methodology with nonlinear response spectra (referred to as strength‐demanded spectra) and a set of formulae of close‐form expressions for the distribution of strength and stiffness produced in the structure by the designed hysteretic damping devices. As an illustrative example, the proposed procedure is applied to a design of a simple steel bridge pier of portal frame type with buckling‐restrained braces (one of several types of hysteretic dampers). For the steel portal frame piers, an attempt is made to utilize not only the displacement‐based index but also the strain‐based index as pre‐determined target performance at the beginning of design. To validate this procedure, dynamic inelastic time‐history analyses are performed using the general‐purpose finite element program ABAQUS. The results confirm that the proposed simplified design procedure attains the expected performance level as specified by both displacement‐based and strain‐based indices with sufficient accuracy. Copyright © 2006 John Wiley & Sons, Ltd.     

A rate‐dependent constitutive model of high damping rubber bearings: modeling and experimental verification

In this study, a constitutive model of high damping rubber bearings (HDRBs) is developed that allows the accurate representation of the force–displacement relationship including rate‐dependence for shear deformation. The proposed constitutive model consists of two hyperelastic springs and a nonlinear dashpot element and expresses the finite deformation viscoelasticity laws based on the classical Zener model. The Fletcher–Gent effect, manifested as high horizontal stiffness at small strains and caused by the carbon fillers in HDRBs, is accurately expressed through an additional stiffness correction factor α in the novel strain energy function. Several material parameters are used to simulate the responses of high damping rubber at various strain levels, and a nonlinear viscosity coefficient η is introduced to characterize the rate‐dependent property. A parameter identification scheme is applied to the results of the multi‐step relaxation tests and the cyclic shear tests, and a three‐dimensional function of the nonlinear viscosity coefficient η with respect to the strain, and strain rate is thus obtained. Finally, to investigate the accuracy and feasibility of the proposed model for application to the seismic response assessment of bridges equipped with HDRBs, an improved real‐time hybrid simulation (RTHS) test system based on the velocity loading method is developed. A single‐column bridge was used as a test bed and HDRBs was physically tested. Comparing the numerical and RTHS results, advantage of the proposed model in the accuracy of the predicted seismic response over comparable hysteretic models is demonstrated. Copyright © 2016 John Wiley & Sons, Ltd.     

Using Bézier curves to model gradual stiffness transitions in nonlinear elements: Application to self‐centering systems

A number of techniques are available for modelling nonlinear elements, but most available hysteretic rules do not capture the gradual stiffness changes that are typical of physical systems. In particular, there has not previously been a hysteretic rule with rounded hysteretic corners that could be used to model self‐centering elements, where multiple stiffness changes occur within one loading cycle. This paper presents a new hysteretic rule that allows the gradual stiffness transitions that occur in real systems to be modelled. In this paper, the rule is formulated for flag‐shaped hystereses, but it is shown that the same model also produces hystereses that can be used to model systems that are not self‐centering. The same technique could also be applied to round the corners of different backbone hystereses. A previous study has shown how abrupt stiffness changes can cause very large acceleration spikes, particularly in self‐centering systems. This paper shows that acceleration spikes due to stiffness changes may be reduced by designing systems to change stiffness more gradually, and that this typically has little effect on other aspects of the seismic response. When modelling structural systems, especially if they are self‐centering, sharp‐cornered hysteretic models may be used for initial analysis, but round‐cornered hysteretic models should be considered when using nonlinear rotational springs or when accelerations are of particular importance. Copyright © 2011 John Wiley & Sons, Ltd.     

An improved capacity spectrum method for ATC‐40

In order to account for the non‐linear behavior of structures via non‐linear static procedure, the capacity spectrum method has been adopted by ATC‐40 for evaluation and retrofit of reinforced concrete buildings. For elastic‐perfectly‐plastic SDOF systems, the accuracy of the capacity spectrum method depends only on the acceleration response spectrum chosen to form the demand spectrum and the adopted model for calculating the equivalent viscous damping ratios. According to this method, the pseudo‐acceleration response spectrum (PS_a) is used to create the demand diagram. It is found that the ATC‐40 procedure, using its Type A hysteretic model, may be inaccurate especially for systems with damping ratios greater than 10% and periods longer than 0.15sec. In order to improve the accuracy of the capacity spectrum method, this study proposes to use the real absolute acceleration response spectrum (S0._a) instead of the PS_a to establish the demand diagram. The step‐by‐step procedure of the improved method and examples are implemented in this paper to illustrate the calculations of earthquake‐induced deformations. In addition, three selected models of equivalent viscous damping are also compared in this paper to assess the accuracy of the model used in the ATC‐40 procedure. Results show that the WJE damping model may be used by the capacity spectrum method to reasonably predict the inelastic displacements when the ductility demand (μ) of the structures is less than 4, whereas the damping model proposed by Kowalsky can be implemented when μ>4.0. Alternatively, the damping model proposed by Kowalsky may be used to calculate the equivalent viscous damping for the entire range of ductility. Copyright © 2003 John Wiley & Sons, Ltd.     

线性土-结构动力相互作用时域-频域联合解法

提出一种线性土-结构动力相互作用时域-频域联合解法．首先，用近场波动数值模拟解耦技术求得在短时脉冲作用下采用Rayleigh阻尼系统的时域解；再对时域解进行富立叶变换得到相应频域解；然后根据阻尼与系统动力反应结果的关系，利用泰勒级数展开技术得到具有复阻尼系统的频域解．这一方法充分利用了时域解耦显式算法的优点，提高了线性土 结构动力相互作用分析的计算效率．      

Free‐vibration analysis of a three‐dimensional soil–structure system

A procedure which involves a non‐linear eigenvalue problem and is based on the substructure method is proposed for the free‐vibration analysis of a soil–structure system. In this procedure, the structure is modelled by the standard finite element method, while the unbounded soil is modelled by the scaled boundary finite element method. The fundamental frequency, and the corresponding radiation damping ratio as well as the modal shape are obtained by using inverse iteration. The free vibration of a dam–foundation system, a hemispherical cavity and a hemispherical deposit are analysed in detail. The numerical results are compared with available results and are also verified by the Fourier transform of the impulsive response calculated in the time domain by the three‐dimensional soil–structure–wave interaction analysis procedure proposed in our previous paper. The fundamental frequency obtained by the present procedure is very close to that obtained by Touhei and Ohmachi, but the damping ratio and the imaginary part of modal shape are significantly different due to the different definition of damping ratio. This study shows that although the classical mode‐superposition method is not applicable to a soil–structure system due to the frequency dependence of the radiation damping, it is still of interest in earthquake engineering to evaluate the fundamental frequency and the corresponding radiation damping ratio of the soil–structure system. Copyright © 2001 John Wiley & Sons, Ltd.     

Extraction of hysteretic properties of seismically isolated bridges from quick‐release field tests

A simple procedure for identifying hysteretic properties of seismically isolated bridges from full‐scale quick‐release tests is presented in this paper. An analytical solution for the quick‐release response of a SDOF system with a bilinear spring is derived. Based on the solution, some characteristics of such systems are obtained. A time domain optimization method is employed to identify the hysteretic properties of the lead–rubber bearings installed in seismically isolated bridges. The total damping effects of the isolation system are expressed as a combination of the rate‐independent (hysteretic) damping and the linear viscous damping. The Menegotto–Pinto (MP) model and bilinear model are used to represent the force–displacement relation of the lead–rubber bearings. In both the longitudinal and transverse directions the bridges have been idealized as single degree of freedom (SDOF) systems. Time histories recorded from the field quick‐release tests on two bridges are used for the examples presented herein. The hysteretic loops of the isolators obtained from laboratory tests are compared with those obtained using the optimization method, and they agree well. In conclusion, the procedure shown in this paper can be used to identify the essential in situ hysteretic characteristics of isolation bearings from quick‐release field testing. Copyright © 2001 John Wiley & Sons, Ltd.     

A linearization scheme for hysteretic systems subjected to random excitation

The problem of predicting the response of hysteretic yielding systems under random excitation is considered. The paper presents a linearization scheme which is motivated by extensive simulation studies of the bilinear hysteretic system. The response process is decomposed into a low frequency drift component and a component with a frequency near the linearized natural frequency of the system. Effective stiffness and damping parameters are presented. The accuracy of the various steps of the scheme is examined. The approach gives results which are in close agreement with simulation estimates even for large-system non-linearity.     

Evaluating foundation mass,damping and stiffness by the least‐squares method

This paper discusses how to use the three‐dimensional (3D) time‐domain finite‐element method incorporating the least‐squares method to calculate the equivalent foundation mass, damping and stiffness matrices. Numerical simulations indicate that the accuracy of these equivalent matrices is acceptable when the applied harmonic force of 1+sine is used. Moreover, the accuracy of the least‐squares method using the 1+sine force is not sensitive to the first time step for inclusion of data. Since the finite‐element method can model problems flexibly, the equivalent mass, damping and stiffness matrices of very complicated soil profiles and foundations can be established without difficulty using this least‐squares method. Copyright © 2003 John Wiley & Sons, Ltd.