桥梁滞变非线性随机地震响应分析

对于含有滞迟剪切型单元的复杂多自由度系统，采用Bouc—Wen微分方程模型描述非线性抗力构件的滞变特性。以虚拟激励法代替常规的Lyapunov方程法结合等效线性化技术对平稳随机激励下系统的动力响应进行求解，激励谱可为白谱或其他任意形式。不但计算效率高，而且不引入除等效线性化之外的附加误差。对不同周期单自由度体系以数值模拟法验证了算法精度。在此基础上，对一座四跨桥进行了线性／非线性分析对比，得到合理的结果。     

随机激励下滞迟系统的稳态响应闭合解

滞迟系统属于一类典型的强非线性系统,滞迟力不仅取决于系统的瞬时变形,还与变形历程有关.虽然滞迟系统的随机振动问题已被广泛研究,但至今尚未得到滞迟系统随机响应概率密度函数的精确闭合解.本文运用迭代加权残值法获得了高斯白噪声激励下Bouc-Wen滞迟系统稳态响应概率密度函数的近似闭合解.首先,运用等效线性化法求出系统的稳态高斯概率密度函数;然后以此构造权函数,应用加权残值法求得了系统指数多项式形式的非高斯概率密度函数;最后引入迭代的过程,逐步优化权函数,提高计算所得结果的精度.以随机地震激励下钢纤维陶粒混凝土结构的稳态响应作为算例,其中Bouc-Wen模型的参数是基于拟静力学试验数据,并应用最小二乘法辨识获得.与Monte Carlo模拟结果相比,等效线性化法得到的结果精度较差;由加权残值法得到的结果能够表现出非线性特征,但其精度依然无法令人满意;采用迭代加权残值法得到的近似闭合解与Monte Carlo模拟的结果吻合非常好;对于较强随机激励情形,采用渐进迭代加权残值法具有较高的求解效率,所获得的理论解析解具有较高的精度.结果表明,所获得的近似闭合解不仅对于土木工程领域具有重要的实际应用价值,而且还可作为检验其他非线性系统随机响应预测方法的精度的标准.     

非高斯随机激励下滞迟系统响应与首次穿越失效

针对由有界噪声、泊松白噪声和高斯白噪声共同构成的非高斯随机激励,通过Monte Carlo数值模拟方法研究了此激励作用下双线性滞迟系统和Bouc-Wen滞迟系统这两类经典滞迟系统的稳态响应与首次穿越失效时间。一方面,分析了有界噪声和泊松白噪声这两种分别具有连续样本函数和非连续样本函数的非高斯随机激励,在不同激励参数条件下对双线性滞迟系统和Bouc-Wen滞迟系统的稳态响应概率密度、首次穿越失效时间概率密度及其均值的不同影响;另一方面,揭示了在这类非高斯随机激励荷载作用下,双线性滞迟系统的首次穿越失效时间概率密度将出现与Bouc-Wen滞迟系统的单峰首次穿越失效时间概率密度截然不同的双峰形式。     

非高斯随机激励下滞迟系统响应与首次穿越失效

针对由有界噪声、泊松白噪声和高斯白噪声共同构成的非高斯随机激励,通过Monte Carlo数值模拟方法研究了此激励作用下双线性滞迟系统和Bouc-Wen滞迟系统这两类经典滞迟系统的稳态响应与首次穿越失效时间。一方面,分析了有界噪声和泊松白噪声这两种分别具有连续样本函数和非连续样本函数的非高斯随机激励,在不同激励参数条件下对双线性滞迟系统和Bouc-Wen滞迟系统的稳态响应概率密度、首次穿越失效时间概率密度及其均值的不同影响;另一方面,揭示了在这类非高斯随机激励荷载作用下,双线性滞迟系统的首次穿越失效时间概率密度将出现与Bouc-Wen滞迟系统的单峰首次穿越失效时间概率密度截然不同的双峰形式。     

多自由度Duffing系统受多相位平稳随机激励的响应

本文用虚拟激励法结合等效线性法来求解受多相位平稳随机激励的多自由度Ｄｕｆｆｉｎｇ系统，过程简洁，计算高效，给出了在该激励下多自由度Ｄｕｆｆｉｎｇ系统的数值模拟算法，二者进行了比较，得到了比较满意的结果。     

含间隙裂缝的RC结构滞回系统混沌振动分析

研究了含间隙裂缝的钢筋混凝土结构对称滞回非线性问题. 建立了 一种分段线性的对称滞回模型, 利用一次谐波线性方法求解结构系统的等效阻尼和 等效刚度系数,得到了对称滞回非线性系统的等价线性方程. 通过数值分析比较了考虑和不考虑间隙与碰撞 影响的两种情况下系统的混沌动力特性,研究表明: 不考虑间隙与碰撞影响的系统出现周期运 动, 考虑间隙与碰撞影响系统更容易出现混沌运动; 在特定的参数范围内系统一定会出现无 序的混沌运动.     

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

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

基于state-space-split法的滞回非线性系统随机动力响应分析

采用SSS(state-space-split)法,建立了引入Bouc-Wen滞回模型的杜芬非线性系统在高斯白噪声激励下的概率密度函数(PDF)的近似求解方法,分析了其随机动力响应变化规律.首先,将Bouc-Wen滞回模型引入杜芬非线性系统,分别考虑非线性系统中的几何非线性和材料非线性对动力响应的影响.随后,对该模型进行了等效线性化(EQL)处理,基于等效线性化法结果,介绍了SSS法的简化方法和计算原理,并通过该方法求解了三维FPK方程的近似联合概率密度函数.最后,将该方法应用于三个案例和一个工程实例,通过求解其概率密度函数分布及动力可靠度验证了提出方法的适用性、可行性和优越性.     

旗帜型滞迟系统的首次穿越失效研究

自复位结构作为一种新型的抗震结构模型受到广泛关注,但这类结构的随机振动问题是技术上的难点,而首次穿越失效作为结构可靠性分析的一个重要模型,则属于非线性随机振动研究中最为困难的问题之一。本文针对随机地震激励下旗帜型滞迟系统的首次穿越失效问题。将旗帜型自复位恢复力在广义谐波平衡技术下应用,近似为幅值相关的等效阻尼和等效刚度组合,从而推导出系统运动方程的等效系统;利用随机平均原理导出关于幅值的平均It?微分方程;采用有限差分法并结合适当的边界条件和初始条件,求解与上述方程对应的后向Kolmogorov(BK)方程,得到首次穿越时间的条件概率密度函数和条件可靠度函数。作为算例,分别讨论了两种随机地震激励模型作用时,系统参数对条件可靠性函数与首次穿越时间的条件概率密度函数的影响。另外,通过与对原方程的蒙特卡罗模拟数据的对比,验证了理论解的有效性。     

等效线性化方法的最优性

等效线性化方法的最优性陈立群（鞍山钢铁学院理论力学教研室鞍山１１４００２）对于弱非线性单自由度自由振动系统的第一次近似用能量平衡原理或谐波平衡原理都可以导出等效线性化系统ｌ‘］其中而且可以证明方程（３）是精确到Ｏ（。）的方程．这里将进一步证明以等效线...     

SEISMIC RANDOM RESPONSE ANALYSIS OF HYSTERETIC SYSTEMS WITH PSEUDO EXCITATION METHOD

The stationary random responses of nonlinear shear-type Multi-Degrees-of-Freedom (MDOF) hysteretic system are analyzed by using the Pseudo Excitation Method (PEM) combined with the Equivalent Linearization Method (ELM). The solution of the equivalent linear system is obtained by iteratively solving complex algebraic equations instead of the Lyapunov equations. The efficiency of this method is much higher for practical engineering systems with many degrees-of-freedom. Supported by NNSFC (Project No. 19772009), NKBRSF fund (No. G1999032805) and Doctoral Funding of State Education Ministry (No. 97014120).     

Principal Lyapunov Exponents and Principal Floquet Spaces of Positive Random Dynamical Systems. III. Parabolic Equations and Delay Systems

This is the third part in a series of papers concerned with principal Lyapunov exponents and principal Floquet subspaces of positive random dynamical systems in ordered Banach spaces. The current part focuses on applications of general theory, developed in the authors’ paper Mierczyński and Shen (Trans Am Math Soc 365(10):5329–5365, 2013), to positive continuous-time random dynamical systems on infinite dimensional ordered Banach spaces arising from random parabolic equations and random delay systems. It is shown under some quite general assumptions that measurable linear skew-product semidynamical systems generated by random parabolic equations and by cooperative systems of linear delay differential equations admit measurable families of generalized principal Floquet subspaces, and generalized principal Lyapunov exponents.     

双面约束多点摩擦多体系统的建模和数值方法

提出了一种建立具有固定双面约束多点摩擦的多体系统动力学方程的方法. 用笛卡尔坐标阵描述系统的位形,根据局部方法的递推关系建立系统的约束方程,应用第一类Lagrange方程建立该系统的动力学方程,使得具有摩擦的约束面的法向力与Lagrange乘子一一对应,便于摩擦力的分析与计算,并用矩阵形式给出了摩擦力的广义力的一般表达式. 应用增广法将微分-代数方程组转化为常微分方程组,并用分块矩阵的形式给出,以便于方程的编程与计算.给出了一种改进的试算法,可提高计算效率. 最后给出了一个算例,应用试算法和RK法对算例进行了数值仿真.      

求解非线性方程组的混合遗传算法

非线性方程组的求解是数值计算领域中最困难的问题。大多数的数值求解算法例如牛顿法的收敛性和性能特征在很大程度上依赖于初始点。但是对于很多非线性方程组,选择好的初始点是一件非常困难的事情。本文结合遗传算法和经典算法的优点,提出了一种用于求解非线性方程组的混合遗传算法。该混合算法充分发挥了遗传算法的群体搜索和全局收敛性,有效地克服了经典算法的初始点敏感问题;同时在遗传算法中引入经典算法(Powell法、拟牛顿迭代法)作局部搜索,克服了遗传算法收敛速度慢和精度差的缺点。选择了几个典型非线性方程组,从收敛可靠性、计算成本和适用性等指标分析对比了不同算法。计算结果表明所设计的混合遗传算法有着可靠的收敛性和较高的收敛速度和精度,是求解非线性方程组的一种成功算法。     

Integration of linear and some model non-linear equations of motion of incompressible fluids

We suggest a new exact method that allows one to construct solutions to a wide class of linear and some model non-linear hydrodynamic-type systems. The method is based on splitting a system into a few simpler equations; two different representations of solutions (non-symmetric and symmetric) are given. We derive formulas that connect solutions to linear three-dimensional stationary and non-stationary systems (corresponding to different models of incompressible fluids in the absence of mass forces) with solutions to two independent equations, one of which being the Laplace equation and the other following from the equation of motion for any velocity component at zero pressure. To illustrate the potentials of the method, we consider the Stokes equations, describing slow flows of viscous incompressible fluids, as well as linearized equations corresponding to Maxwell's and some other viscoelastic models. We also suggest and analyze a differential-difference fluid model with a constant relaxation time. We give examples of integrable non-linear hydrodynamic-type systems. The results obtained can be suitable for the integration of linear hydrodynamic equations and for testing numerical methods designed to solve non-linear equations of continuum mechanics.     

Numerical solution of compressible flow in a channel and blade cascade

This paper deals with a numerical solution of compressible flows. In the case of Euler equations, a numerical solver is presented on a structured quadrilateral grid. The Advection Upstream Splitting Method (AUSM) scheme is used and the spatial accuracy is improved by linear reconstruction with slope limiters. The influence of those limiters are then tested in cases of transonic flow through a channel and a blade cascade.     

Finite Dimensional State Representation of Linear and Nonlinear Delay Systems

We consider the question of when delay systems, which are intrinsically infinite dimensional, can be represented by finite dimensional systems. Specifically, we give conditions for when all the information about the solutions of the delay system can be obtained from the solutions of a finite system of ordinary differential equations. For linear autonomous systems and linear systems with time-dependent input we give necessary and sufficient conditions and in the nonlinear case we give sufficient conditions. Most of our results for linear renewal and delay differential equations are known in different guises. The novelty lies in the approach which is tailored for applications to models of physiologically structured populations. Our results on linear systems with input and nonlinear systems are new.     

Nonholonomic systems with non-linear constraint equations

The development of a form of Lagrange's equations applicable with nonholonomic systems with non-linear constraint equations is presented and discussed. The analysis is based upon, and is an extension of. a method developed by the authors for nonholonomic systems with linear constraint equations in the generalized coordinate derivatives. The method is illustrated with the problem of the “balancing pole”.     

基于改进降维法的可靠度分析

单变量维数缩减法可以高效、准确地进行结构响应矩的分析.与传统的一阶可靠度算法FORM(First Order Reliability Method),二阶可靠度算法SORM(Second Order Reliability Method)相比,该方法不需要响应的导数,也不需要迭代搜索最可能失效点.然而近期的研究发现,该...     

On the use of linear graph theory in multibody system dynamics

Multibody dynamics involves the generation and solution of the equations of motion for a system of connected material bodies. The subject of this paper is the use of graph-theoretical methods to represent multibody system topologies and to formulate the desired set of motion equations; a discussion of the methods available for solving these differential-algebraic equations is beyond the scope of this work. After a brief introduction to the topic, a review of linear graphs and their associated topological arrays is presented, followed in turn by the use of these matrices in generating various graph-theoretic equations. The appearance of linear graph theory in a number of existing multibody formulations is then discussed, distinguishing between approaches that use absolute (Cartesian) coordinates and those that employ relative (joint) coordinates. These formulations are then contrasted with formal graph-theoretic approaches, in which both the kinematic and dynamic equations are automatically generated from a single linear graph representation of the system. The paper concludes with a summary of results and suggestions for further research on the graph-theoretical modelling of mechanical systems.