Mode‐acceleration approach to seismic response of multiply‐supported secondary systems

A mode‐acceleration approach has been proposed for estimating the seismic response of a linear, classically‐damped, multiply‐supported secondary system within the framework of a power spectral density function (PSDF)‐based stochastic approach, while the primary system is linear and classically‐damped. Response transfer functions have been formulated in terms of chosen numbers of fixed‐base modes of the primary and secondary systems. The proposed approach does not involve the determination of combined system properties, and is applicable to the secondary systems with high mass ratios also. Through a few example primary–secondary systems and an example band‐limited white noise excitation, it has been shown that this approach leads to reasonably accurate results when only a few primary and secondary modes are to be considered. The proposed formulation has been used to obtain input data for a decoupled response spectrum analysis of secondary systems. This data accurately accounts for the effects of interaction between the primary and secondary systems. It is shown to lead to substantial reductions in the errors associated with the envelope spectrum method in the case of moderately heavy to heavy secondary systems and when the spatial coupling does not play a major role. Copyright © 2002 John Wiley & Sons, Ltd.     

A pseudo‐force iterative method with separate scale factors for dynamic analysis of structures with non‐proportional damping

Dynamic equilibrium equations of structural systems with non‐proportional damping are coupled through the damping terms. Such coupling invalidates application of the classical modal superposition method. In this paper, a mode‐superposition pseudo‐force method is proposed. The coupled equilibrium equations are solved by an iterative process in which the coupling terms are treated as pseudo‐forces. A scale factor for each mode of the system is obtained by optimizing the iteration convergence. Through these uniquely solved scale factors, the modified modal equations not only converge much faster but also yield results with higher accuracy. A proof of the convergence of the iterative process is also presented. Copyright © 2002 John Wiley & Sons, Ltd.     

2.5D forward modeling and inversion of frequency-domain airborne electromagnetic data

Frequency-domain airborne electromagnetics is a proven geophysical exploration method. Presently, the interpretation is mainly based on resistivity—depth imaging and one-dimensional layered inversion; nevertheless, it is difficult to obtain satisfactory results for two- or three-dimensional complex earth structures using 1D methods. 3D forward modeling and inversion can be used but are hampered by computational limitations because of the large number of data. Thus, we developed a 2.5D frequency-domain airborne electromagnetic forward modeling and inversion algorithm. To eliminate the source singularities in the numerical simulations, we split the fields into primary and secondary fields. The primary fields are calculated using homogeneous or layered models with analytical solutions, and the secondary (scattered) fields are solved by the finite-element method. The linear system of equations is solved by using the large-scale sparse matrix parallel direct solver, which greatly improves the computational efficiency. The inversion algorithm was based on damping least-squares and singular value decomposition and combined the pseudo forward modeling and reciprocity principle to compute the Jacobian matrix. Synthetic and field data were used to test the effectiveness of the proposed method.     

Evaluation of the half‐power bandwidth method to estimate damping in systems without real modes

The objective of this work was to assess the significance of the values of damping obtained applying the half‐power bandwidth method to the frequency response records of the steady‐state response of a system that does not possess real modes either because the damping matrix does not satisfy the orthogonality condition or because its parameters are functions of frequency. A multi‐degree of freedom system with real modes and different types of damping is considered first. A two degree of freedom system with an arbitrary damping matrix, a rigid mass on an elastic foundation subjected to vertical and coupled horizontal/rocking vibrations, and a single degree of freedom model of a building accounting for inertial soil structure interaction effects are considered next in more detail. The results show that the predictions of the method, when applicable compare very well with those provided by approximate formulae and procedures used in practice. Copyright © 2010 John Wiley & Sons, Ltd.     

Response of multiply supported secondary systems to earthquakes in frequency domain

A new formulation of the transfer function has been proposed for the seismic analysis of linear, multiply supported secondary systems. The transfer function for a given response quantity has been formulated by directly using the fixed-base modes of the primary and secondary systems. This approach is exact and does not involve the determination of the combined system properties. Further, it is applicable to the secondary systems with various mass ratios and configurations. A few example primary–secondary systems have been considered to illustrate the proposed formulation in case of different mass ratios. It has also been shown how the proposed formulation can be used to obtain reasonably accurate stochastic estimates of the secondary system responses. © 1998 John Wiley & Sons, Ltd.     

Accelerating the T‐matrix approach to seismic full‐waveform inversion by domain decomposition

A seismic variant of the distorted Born iterative inversion method, which is commonly used in electromagnetic and acoustic (medical) imaging, has been recently developed on the basis of the T‐matrix approach of multiple scattering theory. The distorted Born iterative method is consistent with the Gauss–Newton method, but its implementation is different, and there are potentially significant computational advantages of using the T‐matrix approach in this context. It has been shown that the computational cost associated with the updating of the background medium Green functions after each iteration can be reduced via the use of various linearisation or quasi‐linearisation techniques. However, these techniques for reducing the computational cost may not work well in the presence of strong contrasts. To deal with this, we have now developed a domain decomposition method, which allows one to decompose the seismic velocity model into an arbitrary number of heterogeneous domains that can be treated separately and in parallel. The new domain decomposition method is based on the concept of a scattering‐path matrix, which is well known in solid‐state physics. If the seismic model consists of different domains that are well separated (e.g., different reservoirs within a sedimentary basin), then the scattering‐path matrix formulation can be used to derive approximations that are sufficiently accurate but far more speedy and much less memory demanding because they ignore the interaction between different domains. However, we show here that one can also use the scattering‐path matrix formulation to calculate the overall T‐matrix for a large model exactly without any approximations at a computational cost that is significantly smaller than the cost associated with an exact formal matrix inversion solution. This is because we have derived exact analytical results for the special case of two interacting domains and combined them with Strassen's formulas for fast recursive matrix inversion. To illustrate the fact that we have accelerated the T‐matrix approach to full‐waveform inversion by domain decomposition, we perform a series of numerical experiments based on synthetic data associated with a complex salt model and a simpler two‐dimensional model that can be naturally decomposed into separate upper and lower domains. If the domain decomposition method is combined with an additional layer of multi‐scale regularisation (based on spatial smoothing of the sensitivity matrix and the data residual vector along the receiver line) beyond standard sequential frequency inversion, then one apparently can also obtain stable inversion results in the absence of ultra‐low frequencies and reduced computation times.     

Amplification system for concentrated and distributed energy dissipation devices

Recent analytic, experimental, and practical studies are developing energy dissipation devices combined with amplifying mechanisms (AM) to enhance structural behavior. This research presents the theoretical and experimental development of the eccentric lever‐arm system (ELAS), a new system generically called amplified added damping (AAD), which is a combination of an AM with one or more dampers capable of supporting large deformations. The proposed AM device is a variant of the well‐known lever‐arm system. This work is divided in four parts: (1) kinematics of the ELAS and definition of an equivalent AAD; (2) parametric analysis of a linear single‐story structure with ELAS; (3) numerical analysis of a multi‐degree‐of‐freedom structure with frictional damping with and without AM; and (4) pseudo‐dynamic tests of a full scale asymmetric one story steel structure with and without frictional AAD. Parametric analyses demonstrate that using high‐amplification ratios and low supplemental damping could be a good practice. On the other hand, similar to systems without AMs, dissipation efficiency increases conformably with the stiffness of the secondary structure. As expected, it was observed that deformation was highly concentrated in the flexible edge of the asymmetric test model without damper. Conversely, the structure with frictional AAD clearly showed uniform plane deformation. The implemented AM, which has a large amplifying ratio of α≈11, performed with close accordance with numerical simulations and a high mechanical efficiency (≈95%) using a frictional damper with a very low force capacity. Copyright © 2016 John Wiley & Sons, Ltd.     

Simplified seismic analysis of one‐way asymmetric elastic systems with supplemental damping

This study investigates the effectiveness of the modal analysis using two‐degree‐of‐freedom (2DOF) modal stick to deal with the seismic analysis of one‐way asymmetric elastic systems with supplemental damping. The 2DOF modal stick possessing the non‐proportional damping property enables the modal translation and rotation to not be proportional even at elastic state. The analytical results of one‐storey and three‐storey buildings obtained by the proposed method are compared with those obtained by direct integration of the equation of motion and conventional approximate method, which neglects the off‐diagonal elements in the transformed damping matrix. It is found that the proposed simplified method, compared to conventional approximate methods, can significantly improve the accuracy of the analytical results and, at the same time, without obviously increasing computational efforts. Copyright © 2006 John Wiley & Sons, Ltd.     

Damper placement optimization in a shear building model with discrete design variables: a mixed‐integer second‐order cone programming approach

Supplemental damping is known as an efficient and practical means to improve seismic response of building structures. Presented in this paper is a mixed‐integer programming approach to find the optimal placement of supplemental dampers in a given shear building model. The damping coefficients of dampers are treated as discrete design variables. It is shown that a minimization problem of the sum of the transfer function amplitudes of the interstory drifts can be formulated as a mixed‐integer second‐order cone programming problem. The global optimal solution of the optimization problem is then found by using a solver based on a branch‐and‐cut algorithm. Two numerical examples in literature are solved with discrete design variables. In one of these examples, the proposed method finds a better solution than an existing method in literature developed for the continuous optimal damper placement problem. Copyright © 2013 John Wiley & Sons, Ltd.     

Spectral analysis of systems with non-classical damping using classical mode superposition technique

A spectral method for random vibration analysis of a structural system with non-proportional damping is presented using classical (undamped) mode superposition technique. The method obtains the frequency response function of the system by solving the dynamic equilibrium equations in generalized co-ordinates through an iterative process. The iterative solution is written in closed form and the proof for convergence of the iterative process is given. Numerical examples show the convergence characteristics of the process and an excellent accuracy of the obtained results. The method turns out to be computationally more efficient than the conventional methods of spectral analysis using damped mode shapes and frequencies.     

结构地震反应的频域精细传递矩阵法

传统的传递矩阵法需要对控制微分方程进行求解,获得相应的传递矩阵。公式繁琐、复杂。文中提出将传递矩阵法与精细积分法中的指数矩阵运算技巧结合起来,在频域内对结构进行动力分析。与传统的传递矩阵法相比,无需对微分方程进行求解,只需按照迭代公式进行计算,就可以得到所需要的传递矩阵。这种方法公式简单,理论上可实现任意精度,而且计算效率较高,能够快速、高精度的进行结构的地震反应分析。算例显示了精细传递矩阵法的有效性。     

Three‐dimensional modelling of controlled‐source electromagnetic surveys using an edge finite‐element method with a direct solver

A fully three‐dimensional finite‐element algorithm has been developed for simulating controlled‐source electromagnetic surveys. To exploit the advantages of geometric flexibility, frequency‐domain Maxwell's equations of the secondary electric field were discretised using edge‐based finite elements while the primary field was calculated analytically for a horizontally layered‐earth model. The resulting system of equations for the secondary field was solved using a parallel version of direct solvers. The accuracy of the algorithm was successfully verified by comparisons with integral‐equations and iterative solutions, and the applicability to models containing large conductivity contrasts was verified against published data. The advantages of geometry‐conforming meshes have been demonstrated by comparing different mesh systems to simulate an inclined sheet model. A comparison of the performance between direct and iterative solvers demonstrated the superior efficiency of direct solvers, particularly for multisource problems.     

Non-classical damping in dynamic analysis of base-isolated structures with internal equipment

It has been shown that the use of base isolation not only attenuates the response of a primary structural system but also reduces the response of a secondary system mounted on or within the main structure. The isolation system, superstructure and equipment may be made of different materials with significantly different energy dissipation characteristics such that the damping matrix for the combined system is non-classical and can only be approximately expressed by modal damping ratios if the classical mode method is used for analysis. The object of this paper is to evaluate the accuracy of this procedure in approximating the responses of base-isolated structures and internal equipment. The complex mode method can provide exact solutions to problems with non-classical damping and is used here to find the exact response of the isolation-superstructure-equipment system. The entire system is assumed to be linear elastic with viscous damping and the superstructure is assumed to be proportionally damped so that the deformation of the superstructure can be expressed in terms of its classical modes. Recognizing that the ratio of the equipment mass to the structural mass and the ratio of the stiffness of the isolation system to the superstructural stiffness are both small, perturbation methods are used to find the response. This study shows that the response of base-isolated structures can be determined by the classical mode method to some degree of accuracy, but the higher frequency content is distorted. The equipment response derived by the classical mode method is much smaller than the exact solution so that the complex mode method should be applied to find equipment response.     

Operator‐splitting method for real‐time substructure testing

It has been shown that the operator‐splitting method (OSM) provides explicit and unconditionally stable solutions for quasi‐static pseudo‐dynamic substructure testing. However, the OSM provides only an explicit target displacement but not an explicit target velocity, so that it is essentially an implicit method for real‐time substructure testing (RST) when the velocity‐dependent restoring force is considered. This paper proposes a target velocity formulation based on the forward difference of the predicted displacements so as to render the OSM explicit for RST. The stability and accuracy of the resulting OSM‐RST algorithm are investigated. It is shown that the OSM‐RST is unconditionally stable so long as the non‐linear stiffness and damping are of the softening type (i.e. the tangent stiffness and damping never exceed the initial values). The stability of the OSM‐RST for structures with infinite tangent damping coefficient or stiffness is also proved, and the stability of the method for MDOF structures with a non‐classical damping matrix is demonstrated by an energy criterion. The effects of actuator delay and compensation are analysed based on the bilinear approximation of the actuator step response. Experiments on damped SDOF and MDOF structures verify that the stability of the OSM‐RST is preserved when the experimental substructure generates velocity‐dependent reaction forces, whereas the stability of real‐time substructure tests based on the central difference method is worsened by the damping of the specimen. Copyright © 2005 John Wiley & Sons, Ltd.     

Forced vibration test of a building with semi‐active damper system

The authors developed a semi‐active hydraulic damper (SHD) and installed it in an actual building in 1998. This was the first application of a semi‐active structural control system that can control a building's response in a large earthquake by continuously changing the device's damping coefficient. A forced vibration test was carried out by an exciter with a maximum force of 100 kN to investigate the building's vibration characteristics and to determine the system's performance. As a result, the primary resonance frequency and the damping ratio of a building that the SHDs were not jointed to, decreased as the exciting force increased due to the influence of non‐linear members such as PC curtain walls. The equivalent damping ratio was estimated by approximating the resonance curves using the steady‐state response of the SDOF bilinear hysteretic system. After the eight SHDs were jointed to the building, the system's performance was identified by a response control test for steady‐state vibration. The elements that composed the semi‐active damper system demonstrated the specified performance and the whole system operated well. Copyright © 2000 John Wiley & Sons, Ltd.     

A response‐based decoupling criterion for multiply‐supported secondary systems

It is often infeasible to carry out coupled analyses of multiply‐supported secondary systems for earthquake excitations. ‘Approximate’ decoupled analyses are then resorted to, unless the response errors due to those are significantly high. This study proposes a decoupling criterion to identify such cases where these errors are likely to be larger than an acceptable level. The proposed criterion is based on the errors in the primary system response due to decoupling and has been obtained by assuming (i) the input excitation to be an ideal white noise process, (ii) cross‐modal correlation to be negligible, and (iii) the combined system to be classically damped. It uses the modal properties of the undamped combined system, and therefore, a perturbation approach has been formulated to determine the combined system properties in case of light to moderately heavy secondary systems. A numerical study has been carried out to illustrate the accuracy achieved with the proposed perturbation formulation. The proposed decoupling criterion has been validated with the help of two example primary‐secondary systems and four example excitation processes. Copyright © 2002 John Wiley & Sons, Ltd.     

基于复模态的有限元模型修正算法

针对地下结构地震响应分析中无限地基辐射阻尼问题，引入复模态情况下的具有非简化的堆积阻尼矩阵的阻尼模型，并针对具有集中质量阵的阻尼模型提出了合并与质量有关的阻尼和堆积阻尼的思想，并据此提出了一种修正此类有限元模型的两步法，首先从复模态参数中提取实模态参数，采用基于模态残余力的识别算法修正刚度矩阵，然后根据复模态参数和已得的刚度矩阵来识别阻尼模型中的刚度参与系数和质量阻尼堆积阻尼联合矩阵。     

Performance spectra based method for the seismic design of structures equipped with passive supplemental damping systems

A new direct performance‐based design method utilizing design tools called performance‐spectra (P‐Spectra) for low‐rise to medium‐rise frame structures incorporating supplemental damping devices is presented. P‐Spectra are graphic tools that relate the responses of nonlinear SDOF systems with supplemental dampers to various damping parameters and dynamic system properties that structural designers can control. These tools integrate multiple response quantities that are important to the performance of a structure into a single compact graphical format to facilitate direct comparison of different potential solutions that satisfy a set of predetermined performance objectives under various levels of seismic hazard. An SDOF to MDOF transformation procedure that defines the required supplemental damping properties for the MDOF structure to achieve the response defined by the target SDOF system is also presented for hysteretic, linear viscous and viscoelastic damping devices. Using nonlinear time‐history analyses of idealized shear structures, the accuracy of the transformation procedure is verified. A seismic performance upgrade design example is presented to demonstrate the usefulness of the proposed method for achieving design performance goals using supplemental damping devices. Copyright © 2012 John Wiley & Sons, Ltd.     

Elastic properties of double-porosity rocks using the differential effective medium model

An approach to determining the effective elastic moduli of rocks with double porosity is presented. The double‐porosity medium is considered to be a heterogeneous material composed of a homogeneous matrix with primary pores and inclusions that represent secondary pores. Fluid flows in the primary‐pore system and between primary and secondary pores are neglected because of the low permeability of the primary porosity. The prediction of the effective elastic moduli consists of two steps. Firstly, we calculate the effective elastic properties of the matrix with the primary small‐scale pores (matrix homogenization). The porous matrix is then treated as a homogeneous isotropic host in which the large‐scale secondary pores are embedded. To calculate the effective elastic moduli at each step, we use the differential effective medium (DEM) approach. The constituents of this composite medium – primary pores and secondary pores – are approximated by ellipsoidal or spheroidal inclusions with corresponding aspect ratios. We have applied this technique in order to compute the effective elastic properties for a model with randomly orientated inclusions (an isotropic medium) and aligned inclusions (a transversely isotropic medium). Using the special tensor basis, the solution of the one‐particle problem with transversely isotropic host was obtained in explicit form. The direct application of the DEM method for fluid‐saturated pores does not account for fluid displacement in pore systems, and corresponds to a model with isolated pores or the high‐frequency range of acoustic waves. For the interconnected secondary pores, we have calculated the elastic moduli for the dry inclusions and then applied Gassmann's tensor relationships. The simulation of the effective elastic characteristic demonstrated that the fluid flow between the connected secondary pores has a significant influence only in porous rocks containing cracks (flattened ellipsoids). For pore shapes that are close to spherical, the relative difference between the elastic velocities determined by the DEM method and by the DEM method with Gassmann's corrections does not exceed 2%. Examples of the calculation of elastic moduli for water‐saturated dolomite with both isolated and interconnected secondary pores are presented. The simulations were verified by comparison with published experimental data.     

Problems encountered from the use (or misuse) of Rayleigh damping

Rayleigh damping is commonly used to provide a source of energy dissipation in analyses of structures responding to dynamic loads such as earthquake ground motions. In a finite element model, the Rayleigh damping matrix consists of a mass‐proportional part and a stiffness‐proportional part; the latter typically uses the initial linear stiffness matrix of the structure. Under certain conditions, for example, a non‐linear analysis with softening non‐linearity, the damping forces generated by such a matrix can become unrealistically large compared to the restoring forces, resulting in an analysis being unconservative. Potential problems are demonstrated in this paper through a series of examples. A remedy to these problems is proposed in which bounds are imposed on the damping forces. Copyright © 2005 John Wiley & Sons, Ltd.