随机结构非线性动力响应的概率密度演化分析

提出了随机结构非线性动力响应分析的概率密度演化方法．根据结构动力响应的随机状态方程，利用概率守恒原理，建立了随机结构非线性动力响应的概率密度演化方程．结合Newmark-Beta时程积分方法与Lax-Wendroff差分格式，提出了概率密度演化方程的数值分析方法．通过与Monte Carlo分析方法对比，表明所给出的概率密度演化方法具有良好的计算精度和较小的计算工作量．研究表明：随机结构非线性动力响应概率密度具有典型的演化特征，随着时间增长，概率密度曲线分布趋于复杂．     

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

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

载流圆板在磁场中的随机稳定性分析

本文根据大挠度板壳力学基础理论和电磁弹性力学理论,建立了载流圆板的非线性磁弹性随机振动力学模型,采用伽辽金变分法将其变换成非线性常微分动力学方程．通过拟不可积哈密顿系统的平均理论将该方程等价为一个一维伊藤随机微分方程．通过计算该方程的最大Lyapunov 指数判断该系统的局部随机稳定性,并进一步采用基于随机扩散过程的奇异边界理论判断该系统的全局稳定性．最后通过讨论该系统的稳态概率密度函数图的形状变化讨论了该动力系统的随机Hopf分岔的变化规律,并采用数值模拟对理论分析进行了验证．     

Self-adaptive one-dimensional nonlinear finite element method based on element energy projection method

The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.     

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

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

CHAOTIC TRANSIENTS IN A CURVED FLUID CONVEYING TUBE

I. INTRODUCTIONChaotic motion is a kind of reciprocal non-periodic motion caused by a deterministic system. Itis very sensitive to the initial conditions, apparently random and incapable of long-term prediction.Chaotic transient is such a motion that the system will undergo a final steady state, which is one ofseveral or much more possible steady motions of the system. But this final steady state is very sensitiveto the initial conditions of the system. As it is e?ected by many uncertain…     

STOCHASTIC OPTIMAL CONTROL OF HYSTERETIC SYSTEMS UNDER EXTERNALLY AND PARAMETRICALLY RANDOM EXCITATIONS

A stochastic optimal control method for nonlinear hysteretic systems under exter-nally and/or parametrically random excitations is presented and illustrated with an example ofhysteretic column system.A hysteretic system subject to random excitation is first replaced bya nonlinear non-hysteretic stochastic system.An It stochastic differential equation for the to-tal energy of the system as a one-dimensional controlled diffusion process is derived by usingthe stochastic averaging method of energy envelope.A dynamical programming equation is thenestablished based on the stochastic dynamical programming principle and solved to yield the op-timal control force.Finally,the responses of uncontrolled and controlled systems are evaluatedto determine the control efficacy.It is shown by numerical results that the proposed stochasticoptimal control method is more effective and efficient than other optimal control methods.     

Stochastic finite element analysis of local average random fields

This paper discussed the concepts underlying a new approach to representation and analysis of homogeneous random Fields. An orthogonalized isoparametric local average model of random fields is proposed, and the optimum selection of random variables is determined based on the theory of nonlinear programming. The finite element method is used to develop a methodology for random field analysis. Numerical examples showing the accuracy and efficiency of the method are studied.     

Study of differential control method for solving chaotic solutions of nonlinear dynamic system

In this paper, we focus on the need to solve chaotic solutions of high-dimensional nonlinear dynamic systems of which the analytic solution is difficult to obtain. For this purpose, a Differential Control Method (DCM) is proposed based on the Mechanized Mathematics-Wu Elimination Method (WEM). By sampling, the computer time of the differential operator of the nonlinear differential equation can be substituted by the differential quotient of solving the variable time of the sample. The main emphasis of DCM is placed on substituting the differential quotient of a small neighborhood of the sample time of the computer system for the differential operator of the equations studied. The approximate analytical chaotic solutions of the nonlinear differential equations governing the high-dimensional dynamic system can be obtained by the method proposed. In order to increase the computational efficiency of the method proposed, a thermodynamics modeling method is used to decouple the variable and reduce the dimension of the system studied. The validity of the method proposed for obtaining approximate analytical chaotic solutions of the nonlinear differential equations is illustrated by the example of a turbo-generator system. This work is applied to solving a type of nonlinear system of which the dynamic behaviors can be described by nonlinear differential equations.     

梁式结构受移动荷载作用非平稳随机振动的DQ-PEM方法

针对梁式结构受移动荷载作用的非平稳随机振动问题,提出了一种综合利用微分求积法和虚拟激励法DQ-PEM的新方法。梁式结构受移动荷载作用的振动控制方程为含Dirac函数的偏微分方程,利用微分求积(DQ)-积分求积法(IQ)法将其振动控制方程转化为不含Dirac函数的常微分方程。同时,将表示荷载位置变化的Dirac函数视为移动荷载的非平稳化函数,再结合虚拟激励法的思想,可得梁式结构在确定性荷载作用下的虚拟响应,进而得到其非平稳随机响应。通过工程算例验证了该方法的准确性与有效性,并进一步讨论了不同速度和不同边界条件下梁式结构受移动荷载作用的随机振动问题。     

非线性滞迟系统的随机振动分析

本文采用非线性滞后函数模型，对于粘弹性系统的随机振动问题，应用等效线性化和方差分析的方法进行了分析研究，给出了白噪声激励下的响应计算解。     

Nonlinear vibration of a two-mass system with nonlinear stiffnesses

An analytical approach is developed for nonlinear free vibration of a conservative, two-degree-of-freedom mass–spring system having linear and nonlinear stiffnesses. The main contribution of the proposed approach is twofold. First, it introduces the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation and, more significantly, the treatment a nonlinear differential system by linearization coupled with Newton’s method and harmonic balance method. New and accurate higher-order analytical approximate solutions for the nonlinear system are established. After solving the nonlinear differential equation, the displacement of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, this higher-order Newton–harmonic balance (NHB) method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation unlike the classical harmonic balance method which results in complicated algebraic equations requiring further numerical analysis. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and the classical harmonic balance methods. Two examples of nonlinear two-degree-of-freedom mass–spring system are analyzed and verified with published result, exact solutions and numerical integration data.     

Axisymmetric spreading of a thin power-law fluid under gravity on a horizontal plane

Separable solutions admitted by a nonlinear partial differential equation modeling the axisymmetric spreading under gravity of a thin power-law fluid on a horizontal plane are investigated. The model equation is reduced to a highly nonlinear second-order ordinary differential equation for the spatial variable. Using the techniques of Lie group analysis, the nonlinear ordinary differential equation is linearized and solved. As a consequence of this linearization, new results are obtained.     

Nonlinear forced vibration analysis of postbuckled beams

A numerical solution methodology is proposed herein to investigate the nonlinear forced vibrations of Euler–Bernoulli beams with different boundary conditions around the buckled configurations. By introducing a set of differential and integral matrix operators, the nonlinear integro-differential equation that governs the buckling of beams is discretized and then solved using the pseudo-arc-length method. The discretized governing equation of free vibration around the buckled configurations is also solved as an eigenvalue problem after imposing the boundary conditions and some complicated matrix manipulations. To study forced and nonlinear vibrations that take place around a buckled configuration, a Galerkin-based numerical method is applied to reduce the partial integro-differential equation into a time-varying ordinary differential equation of Duffing type. The Duffing equation is then discretized using time differential matrix operators, which are defined based on the derivatives of a periodic base function. Finally, for any given magnitude of axial load, the pseudo -arc-length method is used to obtain the nonlinear frequencies of buckled beams. The effects of axial load on the free vibration, nonlinear, and forced vibrations of beams in both prebuckling and postbuckling domains for the lowest three vibration modes are analyzed. This study shows that the nonlinear response of beams subjected to periodic excitation is complex in the postbuckling domain. For example, the type of boundary conditions significantly affects the nonlinear response of the postbuckled beams.     

A PREDICT-CORRECT NUMERICAL INTEGRATION SCHEME FOR SOLVING NONLINEAR DYNAMIC EQUATIONS

A new numerical integration scheme incorporating a predict-correct algorithm forsolving the nonlinear dynamic systems was proposed in this paper. A nonlinear dynamic systemgoverned by the equation v=F(v,t) was transformed into the form as v=Hv f(v,t). Thenonlinear part f(v,t) was then expanded by Taylor series and only the first-order term retained inthe polynomial. Utilizing the theory of linear differential equation and the precise time-integrationmethod, an exact solution for linearizing equation was obtained. In order to find the solution of theoriginal system, a third-order interpolation polynomial of v was used and an equivalent nonlinearordinary differential equation was regenerated. With a predicted solution as an initial value andan iteration scheme, a corrected result was achieved. Since the error caused by linearization couldbe eliminated in the correction process, the accuracy of calculation was improved greatly. Threeengineering scenarios were used to assess the accuracy and reliability of the proposed method andthe results were satisfactory.     

TEC结构的三维非线性瞬态温度场分析

热电制冷器(TEC)以其体积小、作用速度快及无噪音等机械制冷无法替代的优点在航空航天和电子工业等领域得到了越来越广泛的应用。本文根据TEC的导热特点,推导了TEC结构稳态温度场的解析解,建立了其瞬态非线性温度场分析的微分方程。利用伽辽金法导出TEC结构热分析的有限元方程,对非线性热分析的有限元方程进行了求解,得到了TEC的稳态温度场和瞬态响应温度场。算例结果表明,本文提出的TEC结构热分析有限元模型具有较高的精度,能够有效地分析TEC的非线性瞬态温度场。     

泊松白噪声激励下强非线性系统的半解析瞬态解

自然界与工程中都普遍存在着随机扰动, 且大多数呈现出固有的非高斯性质, 若采用高斯激励建模可能会导致巨大的误差. 泊松白噪声作为一种典型且重要的非高斯激励模型, 已引起了广泛的关注. 目前, 泊松白噪声激励下系统的动态特性分析主要集中于稳态响应的研究, 而针对瞬态响应的求解难度仍较大, 需进一步发展. 本文引入径向基神经网络, 提出了一种泊松白噪声激励下单自由度强非线性系统瞬态响应预测的高效半解析方法. 首先将广义Fokker-Plank-Kolmogorov (FPK) 方程的瞬态解表示为一组含时变待定权值系数的高斯径向基神经网络; 然后采用有限差分法离散时间导数项, 并结合随机取样技术构造含时间递推式的损失函数; 最后通过拉格朗日乘子法使得损失函数最小化获得时变最优权值系数. 作为算例, 探究了两个经典强非线性系统, 并采用蒙特卡罗模拟方法对解析结果加以验证. 结果表明: 本文方法所获得的瞬时概率密度函数与蒙特卡罗模拟数据吻合地较好, 并且算法具备较高的计算效率. 在系统响应的整个演化过程中, 本文所提方法能够非常有效地捕捉到系统响应在各个时刻下的复杂非线性特征. 此外, 本文方法所获得的高精度半解析瞬态解, 不仅可作为基准解检验其他非线性随机振动分析方法的精度, 对于结构的优化设计也存在巨大的潜在应用价值.      

Stochastic finite element analysis of non-linear plane trusses

—This study considers the responses of geometrically and materially non-linear plane trusses under random excitations. The stress-strain law in the inelastic range is based on an explicit differential equation model. After a total Lagrangian finite element discretization, the nodal displacements satisfy a system of stochastic non-linear ordinary differential equations with right-hand-sides given by random functions of time. The exact solution of the above stochastic differential equation is generally difficult to obtain. To seek an approximate solution with good accuracy and reasonable computational effort, the stochastic linearization method is used to find the first and second statistical moments (i.e. the mean vector and the one-time covariance matrix) of the nodal displacements. Results of simple structures under Gaussian white-noise excitation indicate that the proposed method has good accuracy (generally underestimates the r.m.s. stationary response by 5–14%) and requires only a small fraction of the computation time of the time-history Monte-Carlo method.     

随机杆系结构几何非线性分析的递推求解方法

建立了随机静力作用下考虑几何非线性的随机杆系结构的随机非线性平衡方程. 将和 位移耦合的随机割线弹性模量以及随机响应量表示为非正交多项式展开式，运用传统的摄动方法获 得了关于非正交多项式展式的待定系数的确定性的递推方程. 在求解了待定系数后，利用非 正交多项式展开式和正交多项式展开式的关系矩阵，可以很方便地得到未知响应量的二阶统计矩. 两杆结构和平面桁架拱的算例结果表明，当随机量涨落较大时，递推随机有限元方法比基于 二阶泰勒展开的摄动随机有限元方法更逼近蒙特卡洛模拟结果，显示了该方法对几何非线性 随机问题求解的有效性.