Random vibration of viscoelastic system under broad-band excitations

An analytical method is proposed to study the response of a viscoelastic system with strongly non-linear stiffness force and under broad-band random excitations. The random excitations can be additive, or multiplicative, or both, and they can be stationary or non-stationary with evolutionary spectra. With the proposed method, contributions of the viscoelastic force to both damping and stiffness are taken into account separately, and then the extended version of the stochastic averaging, called the quasi-conservative averaging, is applied to the system to derive the averaged equation of energy envelope. Probability density functions of system responses, such as the total energy, the amplitude, and the state variables, can then be obtained analytically. The accuracy of the method is substantiated by comparing the analytical results with those from Monte Carlo simulations. Effects of parameters in the viscoelastic force and in the non-linear stiffness force on the system responses are also investigated.     

粘弹性随机有限元

以近似不可压缩粘弹增量有限元和摄动法为基础，利用增量法处理遗传积分，应用参数摄动考虑随机性，采用局部平均方法对随机场进行离散，通过相关结构分解减少计算量，发展了一种粘弹性随机有限元方法。研究表明，尽管粘弹性材料本构关系具有时间相依性，其随机摄动格式并不存在“长期项”的影响。应用该方法进行粘弹性结构的随机模拟，程序实施简单，计算效率较高，精度较高。     

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.     

Stability and chaotic motion in columns of nonlinear viscoelastic material

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｐｒｏｂｌｅｍｏｆｔｈｅｄｙｎａｍｉｃａｌｓｔａｂｉｌｉｔｙｏｆｖｉｓｃｏｅｌａｓｔｉｃｓｔｒｕｃｔｕｒｅｓｉｓｍｏｒｅｃｏｍｐｌｉｃａｔｅｄｔｈａｎｔｈａｔｏｆｅｌａｓｔｉｃｓｔｒｕｃｔｕｒｅｓ,ｓｉｎｃｅｔｈｅｓｔｒｅｓｓ＿ｓｔｒａｉｎｒｅｌａｔｉｏｎｏｆｖｉｓｃｏｅｌａｓｔｉｃｍａｔｅｒｉａｌｓｏｆｔｅｎｌｅａｄｓｔｏａｎｉｎｔｅｇｒｏ＿ｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｏｆｍｏｔｉｏｎ .Ｉｎｔｈｅｅａｒｌｙｒｅｓｅａｒｃｈｅｓｏｎｌｉｎｅａｒｖｉｓｃｏ…     

Stationary response of Duffing oscillator with hardening stiffness and fractional derivative

The stationary response of Duffing oscillator with hardening stiffness and fractional derivative under Gaussian white noise excitation is studied. First, the term associated with fractional derivative is separated into the equivalent quasi-linear dissipative force and quasi-linear restoring force by using the generalized harmonic balance technique, and the original system is replaced by an equivalent nonlinear stochastic system without fractional derivative. Then, the stochastic averaging method of energy envelope is applied to the equivalent nonlinear stochastic system to yield the averaged Itô equation of energy envelope, from which the corresponding Fokker–Planck–Kolmogorov (FPK) equation is established and solved to obtain the stationary probability densities of the energy envelope and the amplitude envelope. The accuracy of the analytical results is validated by those from the Monte Carlo simulation of original system.     

变截面悬链线无铰拱应变影响线的解析解

针对变截面悬链线无铰拱应变影响线尚无解析解的现状,通过弹性中心法对其力法方程进行简化,利用Ritter截面变化规律简化变截面悬链线无铰拱的曲线积分,从而推导出变截面悬链线无铰拱应变影响线的闭合解表达式,再将解析结果与有限元分析结果进行对比研究,并对轴力参数展开对比分析。研究结果表明,本文推导得到变截面悬链线无铰拱应变影响线的解析解,数值解析解同有限元结果间最大相对误差小于2%,其轴力影响随矢跨比和测点位置变化而变化,本文公式具有较高的工程精度和计算分析参考价值。     

Dynamical Stability of Viscoelastic Column with Fractional Derivative Constitutive Relation

ＩｎｔｒｏｄｕｃｔｉｏｎＭｏｒｅａｎｄｍｏｒｅｐｒｏｂｌｅｍｓｉｎｓｃｉｅｎｃｅａｎｄｅｎｇｉｎｅｅｒｉｎｇｈａｖｅｃｏｎｃｅｒｎｅｄｔｈｅｓｔａｂｉｌｉｔｙａｎａｌｙｓｉｓｏｆｔｈｅｖｉｓｃｏｅｌａｓｔｉｃｓｔｒｕｃｔｕｒｅｓｄｕｅｔｏｔｈｅｉｒｂｒｏａｄａｐｐｌｉｃａｔｉｏｎｓ .Ｉｎ [1]ＣｅｄｅｒｂａｕｍａｎｄＭｏｎｄａｐｐｌｉｅｄｔｈｅｍｕｌｔｉｐｌｅ＿ｓｃａｌｅｓｍｅｔｈｏｄｔｏｔｒｅａｔｓｔａｂｉｌｉｔｙｏｆａｖｉｓｃｏｅｌａｓｔｉｃｃｏｌｕｍｎｕｎｄｅｒａｐｅｒｉｏｄｉｃａｘｉａｌｌ…     

Interaction between a screw dislocation and a viscoelastic piezoelectric bimaterial interface

The complex variable method is employed to derive analytical solutions for the interaction between a piezoelectric screw dislocation and a Kelvin-type viscoelastic piezoelectric bimaterial interface. Through analytical continuation, the original boundary value problem can be reduced to an inhomogeneous first-order partial differential equation for a single function of location z = x + iy and time t defined in the lower half-plane, which is free of the screw dislocation. Once the initial, steady-state and far-field conditions are known, the solution to the first order differential equation can be obtained. From the solved function, explicit expressions are then derived for the stresses, strains, electric fields and electric displacements induced by the piezoelectric screw dislocation. Also presented is the image force acting on the screw dislocation due to its interaction with the Kelvin-type viscoelastic interface. The derived solutions are verified by comparing with existing solutions for the simplified cases, and various interesting features are observed, particularly for those associated with the image force.     

Nonlinear vibrations and stability of elastic and viscoelastic systems under random stationary loads

The paper deals with numerical analysis of nonlinear vibrations of viscoelastic systems under a stochastic action in the form of a Gaussian stationary process with rational spectral density. The analysis is based on numerical simulation of the original stationary process, numerical solution of the differential equations describing the motion of the system, and computation of the maximum Lyapunov exponent if the stability of this motion is studied. An example of a plate subjected to a random stationary load applied in its plane is used to consider specific issues concerning the application of the proposed method and the peculiarities of the behavior of geometrically nonlinear elastic and viscoelastic stochastic systems. Special attention is paid to the interaction of a deterministic periodic action and a stochastic action from the viewpoint of stability of the system motion. It is shown that in some cases imposing a “colored” noise may stabilize an unstable system subjected to a periodic load.     

NON-AXISYMMETRIC MIXED BOUNDARY VALUE PROBLEM FOR DYNAMIC INTERACTION BETWEEN SATURATED POROUS FOUNDATION AND RECTANGULAR PLATE$^{\bf 1)}$

在同一界面的不同区域具有多种边界条件, 称之为混合边界, 这是一个熟知的力学问题. 对这类问题进行精确分析时, 必须要进行混合边值问题的求解. 而对于一般的三维非轴对称情形, 混合边值问题的求解往往存在数学困难. 本文利用Hilbert定理和双重Fourier变换, 给出了一种求解三维非轴对称混合边值问题的解析方法, 利用该方法对具有混合透水边界的饱和多孔地基上矩形板的振动弯曲进行了解析研究(板与地基接触面为不透水边界, 其余为透水边界). 首先, 基于Kirchhoff理论和Biot多孔介质理论建立矩形板与饱和多孔地基的动力控制方程, 进行耦合求解. 针对板土接触面和非接触面的混合边值问题, 采用双重Fourier变换构造出两对二维对偶积分方程, 以接触应力和接触面孔隙压力为基本未知量, 用Jacobi正交多项式将未知量展开, 再利用Schmidt法对二维对偶积分方程完成求解, 最终推导出板土系统在动力作用下的位移和应力解析式. 通过将本文计算模型退化为单一弹性地基, 与已有研究结果进行对比, 验证了本文方法的正确性和有效性. 最后, 通过数值算例, 对饱和多孔地基上矩形板的动力响应及参数影响做出分析和讨论. 此外, 本文提出的解析法具有一般性, 可广泛应用于复杂接触问题和多场耦合问题的求解.     

饱和多孔地基与矩形板动力相互作用的非轴对称混合边值问题

在同一界面的不同区域具有多种边界条件, 称之为混合边界, 这是一个熟知的力学问题. 对这类问题进行精确分析时, 必须要进行混合边值问题的求解. 而对于一般的三维非轴对称情形, 混合边值问题的求解往往存在数学困难. 本文利用Hilbert定理和双重Fourier变换, 给出了一种求解三维非轴对称混合边值问题的解析方法, 利用该方法对具有混合透水边界的饱和多孔地基上矩形板的振动弯曲进行了解析研究(板与地基接触面为不透水边界, 其余为透水边界). 首先, 基于Kirchhoff理论和Biot多孔介质理论建立矩形板与饱和多孔地基的动力控制方程, 进行耦合求解. 针对板土接触面和非接触面的混合边值问题, 采用双重Fourier变换构造出两对二维对偶积分方程, 以接触应力和接触面孔隙压力为基本未知量, 用Jacobi正交多项式将未知量展开, 再利用Schmidt法对二维对偶积分方程完成求解, 最终推导出板土系统在动力作用下的位移和应力解析式. 通过将本文计算模型退化为单一弹性地基, 与已有研究结果进行对比, 验证了本文方法的正确性和有效性. 最后, 通过数值算例, 对饱和多孔地基上矩形板的动力响应及参数影响做出分析和讨论. 此外, 本文提出的解析法具有一般性, 可广泛应用于复杂接触问题和多场耦合问题的求解.      

Dynamic stability of a viscoelastic rotating shaft under parametric random excitation

This paper deals with dynamic stability of a viscoelastic rotating shaft subjected to a parametric random axial compressive thrust, by using moment Lyapunov exponents and the largest Lyapunov exponents as indicators. The equation of motion for the shaft is derived, which is a system of gyroscopic stochastic differential equations. The method of stochastic averaging is used to decouple the governing equations into Itô equations, from which the moment Lyapunov exponent is obtained by using mathematical transformations only. The largest Lyapunov exponent is obtained through its relation with moment Lyapunov exponents. The effects of various parameters on the stochastic dynamic stability are discussed. The approximate analytical results are confirmed by Monte Carlo simulation.     

First passage of stochastic fractional derivative systems with power-form restoring force

In this paper, the first-passage failure of stochastic dynamical systems with fractional derivative and power-form restoring force subjected to Gaussian white-noise excitation is investigated. With application of the stochastic averaging method of quasi-Hamiltonian system, the system energy process will converge weakly to an Itô differential equation. After that, Backward Kolmogorov (BK) equation associated with conditional reliability function and Generalized Pontryagin (GP) equation associated with statistical moments of first-passage time are constructed and solved. Particularly, the influence on reliability caused by the order of fractional derivative and the power of restoring force are also examined, respectively. Numerical results show that reliability function is decreased with respect to the time. Lower power of restoring force may lead the system to more unstable evolution and lead first passage easy to happen. Besides, more viscous material described by fractional derivative may have higher reliability. Moreover, the analytical results are all in good agreement with Monte-Carlo data.     

粘弹性问题的插值型无单元伽辽金比例边界法

作为一种最近发展起来的半解析数值方法,插值型无单元伽辽金比例边界法不仅无需基本解,且在处理应力奇异性问题和无限域问题时十分有效．为了更有效地求解粘弹性问题,对插值型无单元伽辽金比例边界法应用于此类问题进行了研究,并发展了相应的算法. 通过时域分段展开,将时空耦合的初边值问题转化为一系列递推形式的边值问题,然后采用插值型无单元伽辽金比例边界法进行自适应计算．在径向保持解析特性的基础上,环向采用无单元伽辽金法离散可简化前处理和后处理工作量．此外,改进的插值型移动最小二乘法形函数具有插值性,有效地解决了本质边界条件不能直接施加的困难．最后给出了数值算例,并验证了所提方法的有效性和正确性.     

宽带噪声作用下黏弹性板的矩Lyapunov指数

主要研究了在超音速流中受宽带噪声作用的黏弹性板随机振动系统的矩Lyapunov指数.首先, 采用vonKarman板弯理论, 活塞理论以及Galerkin近似法建立了两个自由度耦合的系统运动的随机微分方程. 其次, 应用随机平均法将四维系统降为二维系统. 接着, 对系统采用极坐标变换,通过Girsanov定理和Feynmann-Kac公式得到后向微分算子. 通过对特征函数进行正交Fourier余弦级数展开得到系统矩Lyapunov指数的近似解析式. 并通过MonteCarle仿真得到系统矩Lyapunov指数的数值解验证了近似解析式的可信性. 最后研究了系统参数、气动力参数以及随机噪声谱密度对黏弹性板稳定性的影响.      

Reliability evaluation and control for wideband noise-excited viscoelastic systems

Reliability of first-passage type for wideband noise-excited viscoelastic systems and the quasi-optimal bounded control strategy for maximizing system reliability are investigated. The viscoelastic term is approximately replaced by equivalent damping and stiffness separately. By using the stochastic averaging method based on the generalized harmonic functions, the averaged Itô stochastic differential equation is obtained for the system amplitude. The associated backward Kolmogorov equation is derived and solved to obtain the system reliability. By applying the dynamic programming principle to the averaged system, the quasi-optimal bounded control is devised by maximizing system reliability. The application of the proposed analytical procedures and the effectiveness of the control strategy are illustrated through one example.     

Stochastic response of an axially moving viscoelastic beam with fractional order constitutive relation and random excitations

A convenient and universal residue calculus method is proposed to study the stochastic response behaviors of an axially moving viscoelastic beam with random noise excitations and fractional order constitutive relationship, where the random excitation can be decomposed as a nonstationary stochastic process, Mittag-Leffler internal noise, and external stationary noise excitation. Then, based on the Laplace transform approach, we derived the mean value function, variance function and covariance function through the Green's function technique and the residue calculus method, and obtained theoretical results. In some special case of fractional order derivative α , the Monte Carlo approach and error function results were applied to check the effectiveness of the analytical results, and good agreement was found. Finally in a general-purpose case, we also confirmed the analytical conclusion via the direct Monte Carlo simulation.     

含随机参数的非线性黏弹性有限元计算

进行了粗粒土与结构接触面单调和循环加载试验，基于宏细观测量结果, 扩展了 损伤概念以 描述该类接触面在受载过程中的物态演化, 及由于物态演化导致的力学特性从初始状态到最终 稳定状态的连续变化过程. 揭示了接触面损伤的细观物理基础主要是接触面内土的颗粒破碎 和剪切压密这两种物态演化；指出接触面的剪胀体应变可以划分为可逆性和不可逆性剪胀体 应变两部分，其中不可逆性剪胀体应变可作为接触面损伤发展的宏观量度，因此其归一化 形式可作为一种损伤因子的定义；提出了建立粗粒土与结构接触面一种损伤本构关系的基本思路.     

The ultratough peeling of elastic tapes from viscoelastic substrates

The peeling of an elastic thin tape from a flat smooth viscoelastic substrate is investigated. Based on a Green function approach and on the translational invariance, a closed form analytical solution is proposed, which takes into account the viscoelastic dissipation in the substrate material.We find that peeling is prevented from taking place, only when the external force is smaller than the one predicted by Kendall's formula for elastic tapes on rigid substrates. However, we also find that, regardless of the value of the applied force, steady state detachment may occur when the elastic tape is sufficiently stiff. In this case, the constant peeling velocity can be modulated by properly defining the geometrical parameters and the material properties of tape and viscoelastic foundation. On the other hand, for relatively high peeling angles or compliant tapes a threshold value of the peeling force is found, above which the steady-state equilibrium is no longer possible and unstable detachment occurs.The present study contributes to shed light on the behavior of pressure sensitive adhesives in contact with viscoelastic substrates like the human skin. At the same time, it can be considered a first step towards a better understanding of the effect of viscoelastic dissipation on the fracture behavior of solids.     

黏弹性接触界面端附近的奇异应力场

研究蠕变加载条件下线黏弹性材料接触界面端附近的奇异应力场问题.考虑接触界面的摩擦,假设界面端的滑移方向不改变,相对滑移量微小,且其与位移同量级,由此线性化局部边界条件,根据对应原理得到Laplace变换域中的界面端应力场,导出时域中奇异应力场的卷积积分表达式.对卷积积分核函数进行数值反演,考虑接触材料的两类组合,一是持久模量具有量级上的差异,另一是持久模量接近相同.算例结果证实核函数可以用准弹性法求得的解析式较准确地近似.在此基础上,利用积分中值定理,并引入各应力分量的修正系数,得到黏弹性奇异应力场的简化式.结合核函数的数值反演结果分析修正系数表达式的取值范围,得到如下结论,若两相接触材料的持久模量相差很大,可以采用准弹性解的解析式较准确地描述界面端的奇异应力场;一般情况下,应力场不存在统一的奇异值和应力强度系数,当采用类似于准弹性解的表达式近似给出黏弹性应力场时,可以估计此近似描述的误差限.文中最后采用有限元分析黏弹性板端部嵌入部位的应力场,算例包括了黏弹性板与弹性金属支承、黏弹性板与黏弹性垫层所形成的滑移接触界面端,利用黏弹性有限元的数值结果验证理论分析所得结论的有效性.