粘弹性结构动力稳定性分析的谐波平衡法

本文分析了线粘弹性结构的长期动力稳定特性。设材料具积分型本构关系，且其松弛模量能用Prony级数描述，将微分-积分型控制方程化成微分型方程，应用谐波平衡法确定动力稳定区域，着重讨论了材料参数及系统振动频率对动力稳定区域的影响发现该类粘弹性结构具有与一般阻尼系统不同的动力稳定特性。文中也将系统平衡法直接应用于微分-积分型控制方程，忽略卷积积分运算后产生的随时间衰减的非谐波项来得到决定动力稳定边界的特征方程，并对两种应用所得结果进行了比较。     

粘弹性桩的混沌运动分析

研究了在轴向载荷和周期性横向载荷共同作用下非线性粘弹性嵌岩桩的混沌运动情况。假定桩和土体分别满足Leaderman非线性粘弹性和线性粘弹性本构关系，得到的运动方程为非线性偏微分．积分方程；利用Galerkin方法将方程简化为非线性常微分一积分方程，同时利用非线性动力系统中的数值方法，进行了数值计算，得到了不同载荷参数、几何参数、材料参数时粘弹性桩发生周期运动、多周期运动及混沌运动的时程曲线、相图、功率谱、Poincare截面图，同时得到了挠度-载荷、挠度-几何参数、挠度-材料参数等分叉图，考察了各种参数的影响。数值结果表明非线性粘弹性桩在一定的条件下可以通过倍周期分叉的方式进入混沌运动状态，且桩的载荷参数、几何参数、材料参数对其运动状态有较大的影响。     

分数阶拟周期Mathieu方程的动力学分析

分数阶微积分有着诸多优异的特点, 目前在动力学领域主要用来提高非线性系统振动特性研究的准确性. 本文在拟周期Mathieu方程的基础上, 引入分数阶微积分理论, 研究了分数阶微分项参数对方程稳定性的影响. 首先, 采用摄动法得到方程稳定区和非稳定区分界线(即过渡曲线)近似表达式, 利用数值方法验证了解析结果的准确性, 图像显示两者吻合较好. 随后, 通过归纳总结不同情况下的过渡曲线近似表达式, 发现在系统中分数阶微分项以等效线性刚度和等效线性阻尼的方式存在. 根据这一特点, 得到了系统等效线性阻尼和等效线性刚度的一般形式, 并且定义了非稳定区域厚度. 最后, 通过数值仿真直观地分析了分数阶微分项参数对方程稳定区域大小和过渡曲线位置的影响. 结果发现, 分数阶微分项不仅具有阻尼特性还具有刚度特性, 并且以等效线性刚度和等效线性阻尼的方式影响着方程稳定区域大小和过渡曲线位置. 合理选择分数阶微分项参数可以使其呈现不同程度的刚度特性或阻尼特性, 方程稳定区域的大小和过渡曲线的位置也因此产生了不同程度的变化.      

具不对称粘弹性支承弹性盘轴转子系统的非线性动力稳定性分析

考虑轴和盘的质量及轴的几何非线性，建立了在周期轴向载荷作用下具不对称粘弹性支承弹性盘轴转子系统的非线性动力学方程．基于Lyapunov运动稳定性理论和Floquet判据，对系统的线性和非线性动力稳定性进行了分析．     

几何缺陷浅拱的动力稳定性分析

研究了几何缺陷对粘弹性铰支浅拱动力稳定性能的影响。从达朗贝尔原理和欧拉-贝努利假定出发推导了粘弹性铰支浅拱在正弦分布突加荷载作用下的动力学控制方程,并采用Galerkin截断法得到了可用龙格-库塔法求解的无量纲化非线性微分方程组。同时引入能有效追踪结构动力后屈曲路径的广义位移控制法,对含几何缺陷浅拱的响应曲线进行几何、材料双重非线性有限元分析。用这两种方法分析了前三阶谐波缺陷对浅拱动力稳定性能的影响,其中动力临界荷载由B-R准则判定。主要结论有:材料粘弹性使浅拱动力临界荷载增大且结构响应曲线与弹性情况差别很大;二阶谐波缺陷影响显著,它使动力临界荷载明显下降且使得浅拱粘弹性动力临界荷载可能低于弹性动力临界荷载。     

粘弹性-弹性层合微悬臂梁的自由振动

吸附膜的粘弹性性质对微梁生物传感器的固有频率有显著影响.首先,在欧拉梁假设下,采用线性粘弹性积分型本构关系和拉普拉斯变换方法,建立了动态识别技术中粘弹性-弹性层合微悬臂梁自由振动的基本方程;其次,采用空域分离变量法和时域微分求导法,获得了积分-偏微分系统的固有频率,并采用求解代数方程的卡尔丹公式和不等式的性质,在材料参数和几何参数张成的高维空间获得了齐次通解的结构;最后,研究了微梁的几何尺寸、吸附膜的粘弹性参数、膜基厚度比和模量比对微梁自由振动的影响.结果表明:吸附膜的粘弹性阻尼效应使得微梁的稳态固有频率低于瞬态固有频率;随着吸附膜松弛时间的减小,微梁瞬态固有频率漂移与稳态固有频率漂移之间的差别逐渐增大;通过控制膜基厚度比或模量比等参数可以使微梁振动进入弱阻尼振动区域.     

黏弹性板的大挠度蠕变屈曲

研究了具有初始小挠度受轴向压载黏弹性板的蠕变屈曲问题，在建立控制方程时，利用了von Karman非线性应变-位移关系，并考虑了初始挠度，用标准线性固体模型描述材料的黏弹性特性，在求解非线性积分方程时，利用梯形公式计算记忆积分式，将非线性积分方程化为非线性代数方程进行数值求解，得到了结构的蠕变变形过程，又将问题退化到小挠度情况进行研究，得到了挠度随时间扩展的解析解，分析了瞬时失稳临界载荷、持久临界载荷的物理意义，讨论了考虑几何非线性对黏弹性板蠕变屈曲的影响。     

不可压饱和粘弹性多孔介质的变分原理及其无网格模拟

基于饱和多孔介质理论,在固相和液相微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,建立了流体饱和粘弹性多孔介质动力响应的若干Gurtin型变分原理,包括Hu-Washizu变分原理.利用所建立的变分原理,导出了流体饱和粘弹性多孔介质动力响应无网格数值模拟的离散控制方程,此方程是一个关于时间的对称微分方程组,便于分析计算.作为数值例子,研究了流体饱和粘弹性多孔柱体的一维动力响应,数值结果揭示了流体饱和粘弹性多孔柱体中波的传播特性以及固相粘性的影响.     

在有限变形条件下损伤粘弹性梁的动力学行为

本文在有限变形条件下，根据损伤粘弹性材料的一种卷积型本构关系和温克列假设，建立了粘弹性基础上损伤粘弹性Timoshenko梁的控制方程。这是一组非线性积分——偏微分方程。为了便于分析，首先利用Galerkin方法对该方程组进行简化，得到一组非线性积分一常微分方程。然后应用非线性动力学中的数值方法，分析了粘弹性地基上损伤粘弹性Timoshenko梁的非线性动力学行为，得到了简化系统的相平面图、Poincare截面和分叉图等。考察了材料参数和载荷参数等对梁的动力学行为的影响。特别，考察了基础和损伤对粘弹性梁的动力学行为的影响。     

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

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

Nonlinear wave propagation in damaged hysteretic materials using a frequency domain-based PM space formulation

In this work, a new and simple numerical approach to simulate nonlinear wave propagation in purely hysteretic elastic solids is presented. Conversely to classical time discretization method, which fully integrates the nonlinear equation of motion, this method utilizes a first-order approximation of the nonlinear strain in order to separate linear and nonlinear contributions. The problem for the nonlinear displacements is then posed as a linear one in which the solid is enforced with nonlinear forces derived from the linear strain. In this manner, a frequency analysis can be easily conducted, leading directly to a well-known frequency spectrum for the nonlinear strain. A mesoscale approach known as Preisach–Mayergoyz space (PM space) is used for the chacterization of the nonlinear elastic region of the solid. A meshless element free Galerkin method is implemented for the discretized equations of motion. Nevertheless, a mesh-based method can be still used as well without loss of generality. Results are presented for bidimensional isotropic plates both in plane stress and in plane strain subjected to harmonic monotone excitation.     

Free and forced nonlinear oscillations of a two-layer composite beam with interface slip

Free and forced flexural nonlinear vibrations of a two-layer beam are investigated. Each beam is assumed to have Euler?CBernoulli kinematics and free-free boundary conditions. The interface allows only nonlinear elastic slip between adjacent sides of the beams, so that the transversal displacement is unique. Free vibrations are considered first by the multiple time scale method, which allows to determine the amplitude dependent nonlinear natural frequencies of the system. It is shown that the nonlinear coefficient of the backbone curve is positive, so that hardening/softening behavior of the interface generates hardening/softening behavior of the whole structure. The modifications of the linear normal modes for moderate excitation amplitudes have been computed. Forced and damped nonlinear oscillations are then considered by the same mathematical method, and the nonlinear frequency response curves are obtained.     

A dynamic method for measuring the critical loads of elastic flat plates

A dynamic method is described for determining the linear buckling loads of elastic, perfectly flat, rectangular plates. The proposed method does not require the application of in-plane loads; it requires only vibrational excitation of the plate. The buckling load is determined from the measured normal modes of vibration. The method is applicable to isotropic as well as anisotropic plates with any type of edge support. The accuracy of the dynamic method was evaluated by tests in which buckling loads of aluminum and graphite fiber-reinforced-epoxy composite plates were determined both by the dynamic method and by imposing static in-plane loads on the plates. The results of the dynamic and static tests agree closely. A. Segall (on leave from RAFAEL, Israel)     

Elastic stability of inhomogeneous thin plates on an elastic foundation

We study the elastic stability of infinite inhomogeneous thin plates on an elastic foundation under in-plane compression. The elastic stiffness constants depend on the coordinate variable in the thickness direction of the plate. The elastic foundation is represented as a Winkler-type model characterized by linear and nonlinear spring constants. First we derive the Föppl–von Kármán equations by taking variations of the elastic strain energy. Next we develop the linear stability analysis of the plate under uniform in-plane compression and explicitly derive the critical loads and wave numbers for particular three cases. The effects of the material inhomogeneity, material orthotropy and loading orthotropy on the critical states are examined independently. Finally, we perform a weakly nonlinear analysis of the plate at the onset of the buckling instability. With the multiple scales method, the amplitude equations for the unstable modes that provide insight into the mode type and its amplitude are derived and then the effect of the material inhomogeneity on buckling modes are evaluated qualitatively.     

粘弹性矩形板的混沌和超混沌行为

从薄板Karman理论的基本假设出发；利用线性粘弹性理论中的Boltzman叠加原理，建立了粘弹性薄板非线性动力学分析的初边值问题，其运动方程是一组非线性积分──微分方程．在空间域上利用Galerkin平均化法之后，得到了变型的非线性积分──微分型的Duffing方程．综合利用动力系统中的多种方法，揭示了粘弹性矩形板在横向周期激励下的丰富的动力学行为，如不动点、极限环、混沌、奇怪吸引子、超混沌等，其中，混沌和超混沌是交替出现的．     

Nonlinear vibration of viscoelastic laminated composite plates

The dynamic behavior of laminated composite plates undergoing moderately large deflection is investigated by considering the viscoelastic properties of the material. Based on von Karman's nonlinear deformation theory and Boltzmann's superposition principle, nonlinear and hereditary type governing equations are derived through Hamilton's principle. Finite element analysis and the method of multiple scales are applied to examine the effect of large amplitude on the dissipative nature as well as on the natural frequency of viscoelastic laminated plates. Numerical experiments are performed for the nonlinear elastic case and linear viscoelastic case to check the validity of the procedure presented in this paper. Limitations of the method are discussed also. It is shown that the geometric nonlinearity does not affect the dissipative characteristics in the cases that have nonlinearity of perturbed order.     

Nonlinear Plane Waves in Finite Deformable Infinite Mooney Elastic Materials

For infinite perfectly elastic Mooney materials, nonlinear plane waves are examined in both two and three dimensions. In two dimensions, longitudinal and shear plane waves are examined, while in three dimensions, longitudinal and torsional plane waves are considered. These exact dynamic deformations, applying to the incompressible perfectly elastic Mooney material, can be viewed as extensions of the corresponding static deformations first derived by Adkins [1] and Klingbeil and Shield [2]. Furthermore, the Mooney strain-energy function is the most general material admitting nontrivial dynamic deformations of this type. For two dimensions the determination of plane wave solutions reduces to elementary mathematical analysis, while in three dimensions an integral of the governing system of highly nonlinear ordinary differential equations is determined. In the latter case, solutions corresponding to particular parameter values are shown graphically. This revised version was published online in June 2006 with corrections to the Cover Date.     

Bifurcation and chaos of the circular plates on the nonlinear elastic foundation

According to the large amplitude equation of the circular plate on nonlinear elastic foundation , elastic resisting force has linear item , cubic nonlinear item and resisting bend elastic item. A nonlinear vibration equation is obtained with the method of Galerkin under the condition of fixed boundary. Floquet exponent at equilibrium point is obtained without external excitation. Its stability and condition of possible bifurcation is analysed. Possible chaotic vibration is analysed and studied with the method of Melnikov with external excitation . The critical curves of the chaotic region and phase figure under some foundation parameters are obtained with the method of digital artificial.     

Axisymmetric shells and plates on tensionless elastic foundations

This paper first describes a finite element method for the large deflection analysis of axisymmetric shells and plates on a nonlinear tensionless elastic foundation. Through the use of discrete data points, any form of nonlinear elastic foundation behaviour can be easily modelled. The analysis is then validated by comparison with existing results for circular plates and beams as the only existing results for shells on tensionless foundations are found to be in error. Following this verification, the analysis is applied to investigate the behaviour of shallow spherical shells subject to a central concentrated load on tensionless linear elastic foundations. A number of insightful conclusions regarding the behaviour of such structure-foundation systems are drawn. The numerical results for shells are believed to be the first correct results, which may be useful in benchmarking results from other sources in the future.     

Large free vibration of thin plates: Hierarchic finite Element Method and asymptotic linearization

The aim of this work is to study the free dynamic response of thin plates characterized by geometrical nonlinearities. To achieve this task, the equation of motion of the plate is first carried out through modeling by hierarchical finite element method whose interpolating shape functions are sinusoidal. Then, the study of the nonlinear vibrations was carried out by the development of asymptotic linearization and equivalent linearization methods in modal space. The nonlinear angular frequencies are successively deduced by exciting the corresponding vibrating mode of the structure. The confrontation of these results to those obtained by the iterative method in the physical space and to those found in the literature, showed a very good agreement between the various methods. From the elementary nonlinear frequencies we showed that there exists an equivalent linear dynamical system characterized by only one equivalent linear stiffness matrix. Numerical experiments were carried out on beams and thin plates of various dimensions ratios and boundary conditions. These numerical test simulations, whether in time space or frequency space, have showed that the nonlinear elastic energy is restored by the equivalent linear dynamical system. Nevertheless, we have to say that the dynamic effects of modes above the excited one are neglected.