非比例阻尼线性体系平稳随机地震响应计算的虚拟激励法

应用复振型分解方法,将非比例阻尼线性体系在地震作用下的动力方程求解问题转化为若干个广义复振子的求解与叠加问题。通过假定地震地面运动为一零均值的平稳随机激励,应用虚拟激励法原理,推导得到了广义复振子动力坐标的解析计算公式,进而得到了以复振型为基础的非比例阻尼线性体系随机地震响应计算的一般实数解析解答。算例证实了这种方法的可靠性及高效率。     

Stochastic Dynamic Response of Multiply Supported Tuned Subsystems to Nonstationary Seismic Excitation

ABSTRACT

Stochastic dynamic response of multiple degree of freedom (MDOF) subsystems that are attached at multiple points to an MDOF supporting structure that is subjected to seismic excitation is obtained by using a modal time domain random vibration approach. An earthquake is modeled as a nonstationary, nonwhite, vector random process that has a realistic spectral shape and finite ground displacement, velocity, and acceleration. Analyses of secondary systems are decoupled from the primary system without loss of accuracy when the cascading assumption is appropriate. This makes the design process of subsystems convenient and efficient. When tuning occurs, cascading assumptions are no longer appropriate. Previously obtained results for the modal properties of combined systems in terms of modal properties of the primary system and subsystems are then utilized. Results account for tuning spatial coupling in terms of stiffness and inertia. Exact analytical expressions are obtained for the elements of the state transition matrix and the evolutionary covariance matrix of various responses. The solution is exact, within the limitations of modal analysis, since all the cross terms between modal coordinates are considered. Finally the capabilities of this approach are illustrated through analysis of an example.     

多自由度柔性结构非线性随机振动响应分析方法

随机等效线性化方法是多自由度结构非线性随机振动分析的有效方法,然而柔性结构隐式的系统方程限制了该方法的应用.本文首先提出了把多自由度柔性结构隐式系统方程显式化的等效非线性系统建立方法,该方法结合模态分析,将系统几何非线性作用等效为模态坐标的高次多项式组合,把多自由度物理系统转化为容易求解的模态系统;在此基础上运用随机等效线性化技术建立柔性结构非线性随机振动响应分析方法.该方法通过引入虚拟激励原理,大幅提高了计算效率,使对多自由柔性结构的非线性随机振动响应分析成为可能.通过算例分析,本文分析结果和精确解与数值模拟解的比较结果表明,本文方法具有较好的计算精度和较高的计算效率,能够应用于实际柔性工程结构的随机振动响应分析.     

Efficient computational procedures for reliability estimates of MDOF-systems

—This paper presents a new numerical code for determining the reliability of multi- degree-of-freedom (MDOF) linear or non-linear models of structure Several examples of application are presented. They are discussed in context with the issue of accuracy and computational efficiency of the new code and of other approximate solution methods such as statistical linearization.     

Drifting response of hysteretic oscillators to stochastic excitation

Under a zero-mean, broad-band, stationary-random load, the symmetric elastic perfectly plastic oscillator and many similar hysteretic systems exhibit a Brownian-like displacement response: asymptotically, the displacement mean vanishes and the displacement variance linearly increases with time. This diverging behavior, often referred to as the drift, is observed even when the excitation power spectrum vanishes at zero frequency, an instance so far lacking a satisfactory modeling within the framework of statistical linearization. The paper presents a linearization-based method which captures the drift in such an instance without requiring any simulation-calibrated parameter. The method combines statistical linearization with stochastic averaging and a generalized van der Pol transformation comprising terms introduced to make allowance for the drift. Model predictions are compared with Monte Carlo estimates for an excitation whose power spectrum vanishes at zero frequency. Good agreement is found for a wide range of excitation levels despite the extremeness of the non-linearity.     

非线性系统非平稳随机响应的时域模态分析法

采用时域模态分析和统计线性化法，得到了一个计算非线性多自由度系统非平稳随机响应的方法。该方法是基于统计线性化参数在一系列微小时间间隔内保持不变，而在这些微小时间间隔的分界点突然改变的假定。考虑了等效线性化系统的时变性；获得了响应协方差矩阵的递推关系。给出了两个算例，并将计算结果与相应的数字模拟结果进行了比较。结果证明该方法是具有满意精度和有效，而且可用于时变系统     

Explicit computational method of dynamic response for non-viscously damped structure systems

Non-viscous damping models in which the damping forces depend on the past history of velocities via convolution integrals over some kernel functions have been raised in many engineering fields. This paper describes an explicit computational method of dynamic response for the non-viscously damped structure systems. The explicit formula is derived using the differential property of convolution and the central difference formula of acceleration. The explicit computational procedure of dynamic response is given in detail. Finally, the dynamic responses of MDOF structure system with double exponential model dampers and SDOF structure system with Gaussian model damper are computed using the proposed explicit method. The accuracy and efficiency are discussed by comparison with other two developed methods.     

Response and bifurcation analysis of a MDOF rotor system with a strong nonlinearity

A new HB (Harmonic Balance)/AFT (Alternating Frequency Time) method is further developed to obtain synchronous and subsynchronous whirling response of nonlinear MDOF rotor systems. Using the HBM, the nonlinear differential equations of a rotor system can be transformed to algebraic equations with unknown harmonic coefficients. A technique is applied to reduce the algebraic equations to only those of the nonlinear coordinates. Stability analysis of the periodic solutions is performed via perturbation of the solutions. To further reduce the computational time for the stability analysis, the reduced system parameters (mass, damping, and stiffness) are calculated in terms of the already known harmonic coefficients. For illustration, a simple MDOF rotor system with a piecewise-linear bearing clearance is used to demonstrate the accuracy of the calculated steady-state solutions and their bifurcation boundaries. Employing ideas from modern dynamics theory, the example MDOF nonlinear rotor system is shown to exhibit subsynchronous, quasi-periodic and chaotic whirling motions.     

An Approximate Analysis of Dry-Friction-Induced Stick-Slip Vibrations by a Smoothing Procedure

This paper deals with a systematic procedure to find both stable and unstable periodic stick-slip vibrations of autonomous dynamic systems with dry friction. In this procedure, the discontinuous friction forces are approximated by smooth functions. Using the simple shooting method with a stiff-ODE solver, in combination with a path following algorithm, branches of periodic solutions are computed for a changing design variable. For testing purposes, both 1 and 2-DOF autonomous block-on-belt models and a 1-DOF autonomous drill string model from literature are investigated. Comparison of the results shows that the smoothing procedure accurately describes the behavior of the discontinuous systems. The proposed procedure can also easily be applied to more complex MDOF models, as well as to nonautonomous dynamic systems.     

考虑参数不确定性的非线性梁随机振动分析

非线性系统的随机振动分析一直是结构动力学领域中的难点,已有一些研究表明基于矩等效的线性化方法在功率谱预测上会得到不恰当的分析结果;另一方面,由于不确定性在实际工程中普遍存在,如果同时考虑非线性和不确定性,更是显著增加了问题难度。本文以具有非线性非理想边界梁为研究对象,基于梁模型的动力学微分方程推导了对应的广义频响函数,并应用Volterra级数理论建立了非线性系统随机振动的谱分析方法,最后,结合蒙特卡洛抽样方法计算了具有参数不确定性非线性梁响应功率谱的均值和方差,讨论了不确定性对结构随机振动响应统计特征的影响。     

An approach to the problem of closure in non-linear stochastic mechanics

An approach combining the method of moment equations and the statistical linearization technique is proposed for analysis of the response of non-linear mechanical systems to random excitation. The adaptive statistical linearization procedure is developed for obtaining a more accurate mean square of responses. For these, a Duffing oscillator and an oscillator with cubic non-linear damping subject to white noise excitation are considered. It is shown that the adaptive statistical linearization proposed yields good accurate results for both weak and strong non-linear stochastic systems.Presented at the First European Solid Mechanics Conference, September 9–13, 1991. Munich, Germany     

二自由度结构受迫非线性振动响应的数值计算

本文研究了由飞机尾部简化而成的二自由度结构受迫非线性振动,由于对复杂的运动方程不作进一步简化,方程中含有强非线性耦合项,首先采用数字滤波技术对系统受确定性激励的响应进行了计算;然后在数字模拟计算的基础上用嫁接的ITO积分方法对系统受随机白噪音激励的响应进行了计算。结果表明,ITO积分的综合嫁接方法有更高的计算效率,结果并表明系统的非线性耦合部分对响应有重大影响,是不能加以忽略的,本文的研究为多自由度非线性结构的响应计算提供了有用的思想和方法。     

Optimal bounded control of stochastically excited MDOF nonlinear viscoelastic systems

A new procedure for designing optimal bounded control of stochastically excited multi-degree-of-freedom (MDOF) nonlinear viscoelastic systems is proposed based on the stochastic averaging method and the stochastic maximum principle. First, the system is formulated as a quasi-integrable Hamiltonian system with viscoelastic terms and each viscoelastic term is replaced approximately by an elastically restoring force and a visco-damping force based on the randomly periodic behavior of the motion of quasi-integrable Hamiltonian system. Thus, a stochastically excited MDOF nonlinear viscoelastic system is converted to an equivalent quasi-integrable Hamiltonian system without viscoelastic terms. Then, by applying stochastic averaging, the system is further reduced to a partially averaged system of less dimension. The adjoint equation and maximum condition for the optimal control problem of the partially averaged system are derived by using the stochastic maximum principle, and the optimal bounded control force is determined from the maximum condition. Finally, the probability and statistics of the stationary response of optimally controlled system are obtained by solving the Fokker–Plank–Kolmogorov equation (FPK) associated with the fully averaged Itô equation of the controlled system. An example is worked out to illustrate the proposed procedure and its effectiveness.     

用虚拟激励法求解非比例阻尼线性体系的非平稳随机地震响应

应用复振型分解方法,将非比例阻尼线性体系在地震作用下的动力方程求解问题转化为若干个广义复振子的求解与叠加问题。通过假定地震地面运动为一零均值的非平稳随机激励,应用虚拟激励法原理,推导得到了广义复振子动力坐标计算的一般公式,进而得到了非比例阻尼线性体系非平稳随机地震响应计算的一般解答。由于可以选择少量共轭复振型的影响进行计算,对于大型复杂非比例阻尼结构,其随机地震响应计算工作量可以大幅度减小。算例证实了这种方法的可靠性及可行性。     

Pseudo-fault induction in engineering structures: Multi-input multi-output control

In a previous paper, it was demonstrated that any linear system could be made to respond to harmonic excitation as if a static nonlinearity of specified type and position were present, this response being obtained at a single predetermined point on the structure. The method requires the excitation of the linear structure by an additional or auxiliary input. In the present paper, the theory is extended to allow the possibility of producing a specified nonlinear response at more than one point on a linear structure. It is shown that N responses can be obtained by specifying N or less auxiliary inputs. The theory is also extended to provide for polynomial damping in addition to stiffness nonlinearity. The theory is validated using numerical simulation of MDOF lumped-parameter systems.     

Stochastic stability for nonlinear systems driven by Lévy noise

This paper is to investigate the stochastic stability for nonlinear systems with Lévy process based on Lyapunov exponents. A method of equivalent linearization is proposed to reduce and simplify the original systems. And the mean square responses are carried out to verify the effectiveness of the proposed approach. Then the Lyapunov exponents will be defined and derived to explore the stochastic stability, and two examples are presented to demonstrate the procedure of equivalent linearization and stochastic stability is considered for these two special examples. The results show that the technique of equivalent linearization can be used to study nonlinear systems excited by Lévy noise.     

CFD结合降阶模型预测阵风响应

传统的阵风响应主要在频域内进行分析,气动载荷基于线性方法计算,不能考虑黏性和跨声速流动影响. 飞机设计需考虑不同频率和不同形状阵风的响应,基于CFD的阵风响应预测由于计算工况太多,工作量巨大. 本文发展了一种CFD结合非定常气动力ARMA(autoregressive-moving-averagemodel)降阶模型的阵风响应分析方法,CFD只要针对给定频率和形状的一种阵风响应进行计算,对获得的气动力时间历程运用线性最小二乘法参数辨识ARMA降阶模型的系数,则对任意频率和形状的阵风,代入降阶模型即可确定该阵风的响应,大大提高了计算效率. 为效验发展的方法,先计算NACA0012翼型在低马赫数0.11的阵风响应,通过对比CFD、ARMA降阶模型及早期发展的不可压阵风响应预测方法的结果,验证了方法的有效性. 再对比CFD、ARMA在跨声速马赫数0.8的阵风响应预测结果,证实所发展的方法对跨声速阵风响应预测亦是有效的.      

Lyapunov function construction for nonlinear stochastic dynamical systems

Though the Lyapunov function method is more efficient than the largest Lyapunov exponent method in evaluating the stochastic stability of multi-degree-of-freedom (MDOF) systems, the construction of Lyapunov function is a challenging task. In this paper, a specific linear combination of subsystems’ energies is proposed as Lyapunov function for MDOF nonlinear stochastic dynamical systems, and the corresponding sufficient condition for the asymptotic Lyapunov stability with probability one is then determined. The proposed procedure to construct Lyapunov function is illustrated and validated with several representative examples, where the influence of coupled/uncoupled dampings and excitation intensities on stochastic stability is also investigated.     

A highly-efficient method for stationary response of multi-degree-of-freedom nonlinear stochastic systems

Analytical and numerical studies of multi-degree-of-freedom(MDOF) nonlinear stochastic or deterministic dynamic systems have long been a technical challenge.This paper presents a highly-efficient method for determining the stationary probability density functions(PDFs) of MDOF nonlinear systems subjected to both additive and multiplicative Gaussian white noises. The proposed method takes advantages of the sufficient conditions of the reduced Fokker-Planck-Kolmogorov(FPK) equation when constructing the trial solution. The assumed solution consists of the analytically constructed trial solutions satisfying the sufficient conditions and an exponential polynomial of the state variables, and delivers a high accuracy of the solution because the analytically constructed trial solutions capture the main characteristics of the nonlinear system. We also make use of the concept from the data-science and propose a symbolic integration over a hypercube to replace the numerical integrations in a higher-dimensional space, which has been regarded as the insurmountable difficulty in the classical method of weighted residuals or stochastic averaging for high-dimensional dynamic systems. Three illustrative examples of MDOF nonlinear systems are analyzed in detail. The accuracy of the numerical results is validated by comparison with the Monte Carlo simulation(MCS) or the available exact solution. Furthermore, we also show the substantial gain in the computational efficiency of the proposed method compared with the MCS.     

A nonlinear finite element formulation for axisymmetric torsion of biphasic materials

In this paper, we present a finite element formulation for describing the large deformation torsional response of biphasic materials, with specific application to prediction of nonlinear coupling between torsional deformation and fluid pressurization in articular cartilage. Due to the use of a cylindrical coordinate system, a particular challenge arises in the linearization of the weak form. The torsional axisymmetric case considered gives rise to additional geometric terms, which are important for the robustness of the numerical implementation and that would not be present in a Cartesian formulation. A detailed derivation of this linearization process is given, couched in the context of a variational formulation suitable for finite element implementation. A series of numerical parametric studies are presented and compared to experimental measurements of the time dependent response of cartilage.