Solution of random structural system subject to nonstationary excitation: transforming the equation with random coefficients to one with deterministic coefficients and random initial conditions

This paper deals with the lower order (first four) nonstationary statistical moments of the response of linear systems with random stiffness and random damping properties subject to random nonstationary excitation modeled as white noise multiplied by an envelope function. The method of analysis is based on a Markov approach using stochastic differential equations (SDE). The linear SDE with random coefficients subject to random excitation with deterministic initial conditions are transformed to an equivalent nonlinear SDE with deterministic coefficients and random initial conditions subject to random excitation. In this procedure, new SDE with random initial conditions, deterministic coefficients and zero forcing functions are introduced to represent the random variables. The joint statistical moments of the response are determined by considering an augmented dynamic system with state variables made up of the displacement and velocity vectors and the random variables of the structural system. The zero time-lag joint statistical moment equations for the augmented state vector are derived from the Itô differential formula. The statistical moment equations are ordinary nonlinear differential equations where hierarchy of moments appear. The hierarchy is closed by the cumulant neglect closure method applied at the fourth order statistical moment level. General formulation is given for multi-degree-of-freedom (MDOF) systems and the performance of the method in problems with nonstationary excitations and large variabilities is illustrated for a single-degree-of-freedom (SDOF) oscillator.     

Lanczos method for dynamic analysis of damped structural systems

In this paper we extend the Lanczos algorithm for the dynamic analysis of structures⁷ to systems with general matrix coefficients. The equations of dynamic equilibrium are first transformed to a system of first order differential equations. Then the unsymmetric Lanczos method is used to generate two sets of vectors. These vectors are used in a method of weighted residuals to reduce the equations of motion to a small unsymmetric tridiagonal system. The algorithm is further simplified for systems of equations with symmetric matrices. By appropriate choice of the starting vectors we obtain an implementation of the Lanczos method that is remarkably close to that in Reference 7, but generalized to the case with indefinite matrix coefficients. This simplification eliminates one of the sets of vectors generated by the unsymmetric Lanczos method and results in a symmetric tridiagonal, but indefinite, system. We identify the difficulties that may arise when this implementation is applied to problems with symmetric indefinite matrices such as vibration of structures with velocity feedback control forces which lead to symmetric damping matrices. This approach is used to evaluate the vibration response of a damped beam problem and a space mast structure with symmetric damping matrix arising from velocity feedback control forces. In both problems, accurate solutions were obtained with as few as 20 Lanczos vectors.     

Non‐iterative equivalent linearization of inelastic SDOF systems for earthquakes in Japan and California

The seismic performance of existing structures can be assessed based on nonlinear static procedures, such as the Capacity Spectrum Method. This method essentially approximates peak responses of an inelastic single‐degree‐of‐freedom (SDOF) system using peak responses of an equivalent linear SDOF model. In this study, the equivalent linear models of inelastic SDOF systems are developed based on the constant strength approach, which does not require iteration for assessing the seismic performance of existing structures. To investigate the effects of earthquake type and seismic region on the equivalent linear models, four ground‐motion data sets—Japanese crustal/interface/inslab records and California crustal records—are compiled and used for nonlinear dynamic analysis. The analysis results indicate that: (1) the optimal equivalent linear model parameters (i.e. equivalent vibration period ratio and damping ratio) decrease with the natural vibration period, whereas they increase with the strength reduction factor; (2) the impacts of earthquake type and seismic region on the equivalent linear model parameters are not significant except for short vibration periods; and (3) the degradation and pinching effects affect the equivalent linear model parameters. We develop prediction equations for the optimal equivalent linear model parameters based on nonlinear least‐squares fitting, which improve and extend the current nonlinear static procedure for existing structures with degradation and pinching behavior. Copyright © 2010 John Wiley & Sons, Ltd.     

Vibration suppression of structures

A procedure for vibration suppression of elastic structures and systems is introduced in this paper. First, the equations of the response of the system to external forces, based on a discretized linear model, are developed. Next, the vibration suppression is formulated as a minimization problem with control (external) forces as unknowns. Then, the optimality condition is applied to this problem and a set of linear simultaneous equations is obtained for the unknown forces. The scheme is applied to a number of simple examples for illustration. The results show that the procedure is very effective and, at the same time, very flexible compared to methods that are based on continuous models (partial differential equations). The actual (hardware) implementation can be achieved by means of either passive absorbers or active forces (actuators).     

Direct use of peak ground motion parameters for the estimation of inelastic displacement ratio of SDOF systems subjected to repeated far fault ground motions

This study is aimed at developing statistical equations to estimate the inelastic displacement ratio of singledegree- of-freedom systems subjected to far fault repeated earthquakes. In the study, peak ground motion parameters are used to define the scatter of the original data. The ratio of peak ground acceleration to peak ground velocity, and peak ground velocity of the ground motion records and structural parameters such as period of vibration and lateral strength ratio are used in the proposed equations. For the development of the equations, nonlinear time history analyses of single-degree-offreedom systems are conducted. Then, the results are used in a multivariate regression procedure. The equations are verified by comparing the estimated results with the calculated results. The average error and coefficient of variation of the proposed equations are presented. The analyses results revealed that the direct use of peak ground motion parameters for the estimation of inelastic displacement ratio significantly reduced the scatter in the original data and yielded accurate results. From the comparative results it is also observed that results obtained using equations specific to peak ground velocity or peak ground acceleration to peak ground velocity ratio are similar.     

A parallel,fully coupled,fully implicit solution to reactive transport in porous media using the preconditioned Jacobian-Free Newton-Krylov Method

Modeling large multicomponent reactive transport systems in porous media is particularly challenging when the governing partial differential algebraic equations (PDAEs) are highly nonlinear and tightly coupled due to complex nonlinear reactions and strong solution-media interactions. Here we present a preconditioned Jacobian-Free Newton-Krylov (JFNK) solution approach to solve the governing PDAEs in a fully coupled and fully implicit manner. A well-known advantage of the JFNK method is that it does not require explicitly computing and storing the Jacobian matrix during Newton nonlinear iterations. Our approach further enhances the JFNK method by utilizing physics-based, block preconditioning and a multigrid algorithm for efficient inversion of the preconditioner. This preconditioning strategy accounts for self- and optionally, cross-coupling between primary variables using diagonal and off-diagonal blocks of an approximate Jacobian, respectively. Numerical results are presented demonstrating the efficiency and massive scalability of the solution strategy for reactive transport problems involving strong solution-mineral interactions and fast kinetics. We found that the physics-based, block preconditioner significantly decreases the number of linear iterations, directly reducing computational cost; and the strongly scalable algebraic multigrid algorithm for approximate inversion of the preconditioner leads to excellent parallel scaling performance.     

求解薄板频率特性的边界元法解析研究

根据弹性薄板理论和振动理论,运用边界元法(BEM)分析了薄板弯曲自由振动的频率特性,从薄板自由振动的微分方程式出发,采用动态基本解推导出其相应的边界积分方程式,计算中对边界积分方程进行缩减,应用频率扫描的方法研究薄板的动态特性,准确、有效地得到了其相应的固有频率、给出的数值算例证明了本文的方法简单、可行,且具有良好的解析计算精度。     

New approaches for non‐classically damped system eigenanalysis

This paper presents three new approaches for solving eigenvalue problems of non‐classically damped linear dynamics systems with fewer calculations than the conventional state vector approach. In the latter, the second‐order differential equation of motion is converted into a first‐order system by doubling the size of the matrices. The new approaches simplify the approach and reduce the number of calculations. The mathematical formulations for the proposed approaches are presented and the numerical results compared with the existing method by solving a sample problem with different damping properties. Of the three proposed approaches, the expansion approach was found to be the simplest and fastest to compute. Copyright © 2005 John Wiley & Sons, Ltd.     

全埋入单桩基础横向振动模态正交性研究

采用Timoshenko梁模型作为挠曲波在全埋入单桩基础中的传播模型,根据频域中全埋入单桩基础横向振动时的运动微分方程,结合各节点的力平衡、位移协调方程及局部对偶坐标系中内力及位移的对偶变换关系,严格推导全埋入单桩基础横向振动的模态正交性条件。以回传射线矩阵法为基础,求解全埋入单桩基础横向振动的自振频率、衰减系数及不同自振频率所对应的模态,并通过具体算例验证全埋入单桩基础横向振动模态的正交性条件公式推导的正确性。     

A parcel formulation of Hamiltonian layer models

Starting from the three-dimensional hydrostatic primitive equations, we derive Hamiltonian N-layer models with isentropic tropospheric and isentropic or isothermal stratospheric layers. Our construction employs a new parcel Hamiltonian formulation which describes the fluid as a continuum of Hamiltonian ordinary differential equations bound together by integral transport laws. In particular, we show that this parcel Hamiltonian structure is compatible with the stacking of layers under isentropic or isothermal constraints. The appeal of the parcel formulation is the simplification of various calculations, in particular the derivation of the continuum Poisson bracket and the proof of the Jacobi identity. A comparison and connection is made between the Hamiltonian dynamics of fluid parcels and the Hamiltonian system of partial differential equations. The parcel formulation can be seen as a precursor and tool for the study of Hamiltonian numerical schemes.     

土-桩-框架结构非线性相互作用的精细数值模型及其验证

利用有限元软件ABAQUS,建立了土-桩-框架结构非线性相互作用（SSI）的二维精细有限元模型,分别采用记忆型粘塑性嵌套面模型和损伤塑性模型模拟土体和混凝土材料,采用梁单元和rebar单元模拟RC桩基及其内部纵筋,采用接触面对法模拟桩土接触效应,取得了良好的计算效果。将自由场、框架、土-桩-框架结构模型的分析结果和其它成熟的计算软件进行对比,验证了数值模型的有效性。分析发现：桩基外侧靠近承台处的土体的非线性反应很强烈,而桩基内部土体的非线性反应较小,很大程度上只是跟随群桩一起运动。由于桩土动力接触,桩顶的加速度反应可能超出上部结构,并且桩顶的加速度时程曲线上有非常明显的＂针＂状突变。随着地震动强度的增加,上部框架逐渐表现出单自由度体系的动力特征,加速度反应谱有从多个波峰退化为单一波峰的趋势。     

A nearly analytic symplectic partitioned Runge–Kutta method based on a locally one-dimensional technique for solving two-dimensional acoustic wave equations

In this paper, we develop a new nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique for numerically solving two-dimensional acoustic wave equations. We first split two-dimensional acoustic wave equation into the local one-dimensional equations and transform each of the split equations into a Hamiltonian system. Then, we use both a nearly analytic discrete operator and a central difference operator to approximate the high-order spatial differential operators, which implies the symmetry of the discretized spatial differential operators, and we employ the partitioned second-order symplectic Runge–Kutta method to numerically solve the resulted semi-discrete Hamiltonian ordinary differential equations, which results in fully discretized scheme is symplectic unlike conventional nearly analytic symplectic partitioned Runge–Kutta methods. Theoretical analyses show that the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique exhibits great higher stability limits and less numerical dispersion than the nearly analytic symplectic partitioned Runge–Kutta method. Numerical experiments are conducted to verify advantages of the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique, such as their computational efficiency, stability, numerical dispersion and long-term calculation capability.     

Dynamic nonlinear response of pile foundations under vertical vibration—Theory versus experiment

The influence of nonlinearity on the dynamic response of cast-in-situ reinforced concrete piles subjected to strong vertical excitation was studied. Forced vibration test of single piles (L/d=10, 15, 20) and 2×2 pile groups (s/d=2, 3, 4 for each L/d) were conducted in the field for two different embedded conditions of pile cap. From the measured nonlinear response curves, the effective pile–soil system mass, stiffness and damping were determined and the nonlinear response curves were back-calculated using the theory of nonlinear vibration. The test results were compared with the continuum approach of Novak with dynamic interaction factor approach using both linear and linear-equivalent numerical methods. Reasonable match between the measured and predicted response was found for linear-equivalent methods by introducing a weak boundary-zone around the pile to approximately account for the nonlinear behaviour of pile–soil system. The test data were used to establish the empirical relationship in order to estimate the extent of soil separation around the pile with soil under vertical vibration.     

无矩阵迭代法在膜结构风振耦合分析中的应用

提出在采用浸入物体法（IOM）对膜结构和空气流体建模时,可以采用带有预定条件的无矩阵Newton-Krylov迭代算法求解浸入物体法,并引入了预定条件矩阵。将提出的无矩阵迭代方法应用于一双坡型膜结构的风振耦合分析中,得出了结构的风压和风速分布,并对带有预定条件和不带预定条件的无矩阵迭代算法进行了对比。结果表明,将带有预定条件的无矩阵迭代算法应用于膜结构风振的耦合分析中,可以得到准确结果,并使计算效率大大提高。     

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

Based on the Hilbert–Huang spectral analysis, a method is proposed to identify multi‐degree‐of‐freedom (MDOF) linear systems using measured free vibration time histories. For MDOF systems, the normal modes have been assumed to exist. In this method, the measured response data, which are polluted by noises, are first decomposed into modal responses using the empirical mode decomposition (EMD) approach with intermittency criteria. Then, the Hilbert transform is applied to each modal response to obtain the instantaneous amplitude and phase angle time histories. A linear least‐square fit procedure is proposed to identify the natural frequency and damping ratio from the instantaneous amplitude and phase angle for each modal response. Based on a single measurement of the free vibration time history at one appropriate location, natural frequencies and damping ratios can be identified. When the responses at all degrees of freedom are measured, the mode shapes and the physical mass, damping and stiffness matrices of the structure can be determined. The applications of the proposed method are illustrated using three linear systems with different dynamic characteristics. Numerical simulation results demonstrate that the proposed system identification method yields quite accurate results, and it offers a new and effective tool for the system identification of linear structures in which normal modes exist. Copyright © 2003 John Wiley & Sons, Ltd.     

三维井间电磁场的正反演计算

The method of modified local nonlinear iteration （MLNI） for volume integral equations is introduced and used to calculate forward and inverse problems of 3 D cross hole electromagnetic fields. The large scale scatters are divided into two kinds, the near region scatters and the far region scatters, whose locations and sizes are determined by the position of the field point. The influence of near region scatters is computed by the method of local nonlinear approximation, while that of the far region scatters is treated as external exciting source and is calculated reiteratively. The method has the advantages, such as faster calculating speed, less computer memory requirements and more efficient convergency. The nonlinear inversion algorithm based on MLNI is applied to inversion problems which can concentrate the imaging region to a particular range instead of whole 3 D domain. Since the algorithm considers the effect of the nonlinear term in each element of the Jacobi matrix, it is more accurate and stable in comparison with the convenitional iterative Born inversion method. The double cell decomposition method is applied during the first rough imaging process considering the limitation of the information and computer memory, then the imaging range is reduced and the second more accurate imaging process can be carried out. Numerical results indicate that MLNI is an efficient method to compute the scattering fields of large scale scatters and that the application of the method to inversion can give high resolution 3 D cross hole conductivity images.     

Random vibration of hysteretic systems under bi-directional ground motions

A random vibration method is developed for the response analysis of hysteretic structural systems under stochastic two-dimensional earthquake excitations. The biaxial hysteretic restoring force is modelled by coupled non-linear differential equations; the response statistics are obtained using the equivalent linearization technique. The validity of the proposed model is appraised using available biaxial loading tests of reinforced concrete columns. A parametric study was performed to examine the significance of the effects of biaxial interaction under earthquake excitations. A practical method to evaluate the extreme response statistics is also presented.     

弦支穹顶结构动力反应分析

以天津开发区商务中心大堂的弦支穹顶为研究对象,分析了弦支穹顶结构体系的自振特性,分别运用随机模拟风振分析方法和时间历程分析方法,对其进行了风振和地震反应分析,得到了结构在动力荷载作用下的响应时程,并对分析结果进行了频谱分析和统计分析。研究发现弦支穹顶结构的自振频率呈密集型分布,且振型复杂;结构的风振响应基本以受迫振动为主,没有出现明显的峰值共振现象;结构的地震响应在前几阶基频处出现了较为明显的峰值共振现象;从振动的幅值角度看,风荷载的动力作用效应相对于地震荷载要显著。     

Bivariate Gaussian mixture–based equivalent linearization method for stochastic seismic analysis of nonlinear structures

To address challenges in stochastic seismic analysis of nonlinear structures, this paper further develops a recently proposed Gaussian mixture–based equivalent linearization method (GM‐ELM). The GM‐ELM uses a Gaussian mixture distribution model to approximate the probabilistic distribution of a nonlinear system response. Using properties of the Gaussian mixture model, GM‐ELM can decompose the non‐Gaussian response of a nonlinear system into multiple Gaussian responses of linear single–degree of freedom oscillators. With the set of the equivalent linear systems identified by GM‐ELM, response statistics as crossing rate and first‐passage probability can be computed conveniently using theories of linear random vibration analysis. However, the original version of GM‐ELM may lead to an inaccurate estimate because of the heuristic parameters of the linear system introduced to supplement insufficient information. To overcome this limitation and define unique equivalent linear systems, this paper proposes a further developed version of GM‐ELM, which uses a mixture of bivariate Gaussian densities instead of univariate models. Moreover, to facilitate the use of elastic response spectra for estimating the mean peak responses of a nonlinear structure, a new response spectrum combination rule is proposed for GM‐ELM. Two numerical examples of hysteretic structural systems are presented in this paper to illustrate the application of the bivariate GM‐ELM to nonlinear stochastic seismic analysis. The analysis results obtained by the bivariate GM‐ELM are compared with those obtained by the univariate GM‐ELM, the conventional equivalent linearization method, the tail equivalent linearization method, and Monte Carlo simulation. The supporting source code and data are available for download at https://github.com/yisangri/GitHub‐bGM‐ELM‐code.git     

A nonlinear SDOF model for the simplified evaluation of the displacement demand of low-rise URM buildings

The determination of displacement demands for masonry buildings subjected to seismic action is a key issue in the performance-based assessment and design of such structures. A technique for the definition of single-degree-of-freedom (SDOF) nonlinear systems that approximates the global behaviour of multi-degree-of-freedom (MDOF) 3D structural models has been developed in order to provide useful information on the dependency of displacement demand on different seismic intensity measures. The definition of SDOF system properties is based on the dynamic equivalence of the elastic properties (vibration period and viscous damping) and on the comparability with nonlinear hysteretic behaviour obtained by cyclic pushover analysis on MDOF models. The MDOF systems are based on a nonlinear macroelement model that is able to reproduce the in-plane shear and flexural cyclic behaviour of pier and spandrel elements. For the complete MDOF models an equivalent frame modelling technique was used. The equivalent SDOF system was modelled using a suitable nonlinear spring comprised of two macroelements in parallel. This allows for a simple calibration of the hysteretic response of the SDOF by suitably proportioning the contributions of flexure-dominated and shear-dominated responses. The comparison of results in terms of maximum displacements obtained for the SDOF and MDOF systems demonstrates the feasibility and reliability of the proposed approach. The comparisons between MDOF and equivalent SDOF systems, carried out for several building prototypes, were based on the results of time-history analyses performed with a large database of natural records covering a wide range of magnitude, distance and local soil conditions. The use of unscaled natural accelerograms allowed the displacement demand to be expressed as a function of different ground motion parameters allowing for the study of their relative influence on the displacement demand for masonry structures.