一种考虑初始缺陷影响的非线性梁单元

在目前广泛应用的梁单元中,尚缺乏全面考虑以下四种因素影响的粱单元：(1)轴力的影响;(2)剪切交形的影响;(3)初始弯曲的影响;(4)弯曲变形对轴向应变的影响,即弓形效应。事实上,以上四种非线性因素都会对钢框架结构的稳定和极限承载力有影响,需同时考虑。本文将致力于推导同时考虑以上四种因素影响的平面梁单元的平衡微分方程,最后得到精确的粱单元刚度矩阵,并研究以上四种因素对钢框架构件及钢框架结构的影响。     

Nonlinear coupling of transverse modes of a fixed–fixed microbeam under direct and parametric excitation

Tuning of linear frequency and nonlinear frequency response of microelectromechanical systems is important in order to obtain high operating bandwidth. Linear frequency tuning can be achieved through various mechanisms such as heating and softening due to DC voltage. Nonlinear frequency response is influenced by nonlinear stiffness, quality factor and forcing. In this paper, we present the influence of nonlinear coupling in tuning the nonlinear frequency response of two transverse modes of a fixed–fixed microbeam under the influence of direct and parametric forces near and below the coupling regions. To do the analysis, we use nonlinear equation governing the motion along in-plane and out-of-plane directions. For a given DC and AC forcing, we obtain static and dynamic equations using the Galerkin’s method based on first-mode approximation under the two different resonant conditions. First, we consider one-to-one internal resonance condition in which the linear frequencies of two transverse modes show coupling. Second, we consider the case in which the linear frequencies of two transverse modes are uncoupled. To obtain the nonlinear frequency response under both the conditions, we solve the dynamic equation with the method of multiple scale (MMS). After validating the results obtained using MMS with the numerical simulation of modal equation, we discuss the influence of linear and nonlinear coupling on the frequency response of the in-plane and out-of-plane motion of fixed–fixed beam. We also analyzed the influence of quality factor on the frequency response of the beams near the coupling region. We found that the nonlinear response shows single curve near the coupling region with wider width for low value of quality factor, and it shows two different curves when the quality factor is high. Consequently, we can effectively tune the quality factor and forcing to obtain different types of coupled response of two modes of a fixed–fixed microbeam.

     

考虑变形耦合的几何非线性空间梁单元

以"精确几何模型梁单元"为代表的很多几何非线性梁单元,在构造过程中分别对描述截面转动的转角和描述截面形心位置的位移进行了独立插值,由此引起了诸如运动学描述冗余和剪切闭锁等困难。其根本原因在于单元形函数没能体现细长梁中的变形耦合关系。本文对这类传统单元进行了改造,通过深入研究单元变形之间的内在联系,提出了一种变形场完全满足Bernoulli梁变形耦合关系的新单元,避免了构造过程中对转动矢量的插值,并通过数值算例检验了单元的有效性。     

Effects of rotary inertia on sub- and super-critical free vibration of an axially moving beam

The most important issue in the vibration study of an engineering system is dynamics modeling. Axially moving continua is often discussed without the inertia produced by the rotation of the continua section. The main goal of this paper is to discover the effects of rotary inertia on the free vibration characteristics of an axially moving beam in the sub-critical and super-critical regime. Specifically, an integro-partial-differential nonlinear equation is modeled for the transverse vibration of the moving beam based on the generalized Hamilton principle. Then the effects of rotary inertia on the natural frequencies, the critical speed, post-buckling vibration frequencies are presented. Two kinds of boundary conditions are also compared. In super-critical speed range, the straight configuration of the axially moving beam loses its stability. The buckling configurations are derived from the corresponding nonlinear static equilibrium equation. Then the natural frequencies of the post-buckling vibration of the super-critical moving beam are calculated by using local linearization theory. By comparing the critical speed and the vibration frequencies in the sub-critical and super-critical regime, the effects of the inertia moment due to beam section rotation are investigated. Several interesting phenomena are disclosed. For examples, without rotary inertia, the study overestimates the stability of the axially moving beam. Moreover, the relative differences between the super-critical fundamental frequencies of the two theories may increase with an increasing beam length.     

Nonlinear finite element analysis of inflatable beams made from orthotropic woven fabric

This paper was devoted to the three-dimensional nonlinear finite element analysis of inflatable beams. The beams under consideration are made of modern textile materials and can be used as a load-bearing beams or arches when inflated. A 3D Timoshenko beam with a homogeneous orthotropic woven fabric (OWF) was proposed. The model took into account the geometric nonlinearities and the follower force resulting from the inflation pressure. The use was made of the usual total Lagrangian form of the virtual work principle to perform the nonlinear equilibrium equations which were discretized by the finite element method. Two kinds of solutions were then investigated: finite elements solutions for linearized problems which were obtained by the means of the linearization around the prestressed reference configuration of the nonlinear equations and nonlinear finite element solutions which were performed by the use of an optimization algorithm based on the Quasi-Newton method. As an example, the bending problem of a cantilever inflated beam under concentrated load was considered and the deflection results improve the existing theoretical models. As these beams are made from fabric, the beam models were validated through their comparison with a 3D thin-shell finite element model. The influence of the material effective properties and the inflation pressure on the beam response was also investigated through a parametric study. The finite elements solutions for linearized problems were found to be close to the theoretical results existing in the literature. On the other hand, the results for the nonlinear finite element model were shown to be close to the results for the linearized finite elements model in the case of high mechanical properties and the nonlinear finite element model was used to improve the linearized model when the mechanical properties of the fabric are low.     

横向载荷作用下轴向运动梁的屈曲失稳以及非线性振动特性

梁的轴向运动会诱发其产生横向振动并可能导致屈曲失稳,对结构的安全性和可靠性产生重大的影响。本文重点研究了横向载荷作用下轴向运动梁的屈曲失稳及横向非线性振动特性。基于Hamilton变分原理,建立了横向载荷作用下轴向运动梁的动力学方程,获得了梁的后屈曲构型。使用截断Galerkin法,将控制方程改写成Duffing方程的形式。用同伦分析方法确定载荷作用下轴向运动梁的非线性受迫振动的封闭形式的表达式。结果表明,后屈曲构型对轴向速度和初始轴向应力有明显的依赖性。通过同伦分析法得出非线性基频的显式表达式,获得了初始轴向力会影响非线性频率随初始振幅和轴向速度的线性关系。另外,轴向外激励的方向也会改变系统固有频率。     

A nonlinear planar beam formulation with stretch and shear deformations under end forces and moments

A new nonlinear planar beam formulation with stretch and shear deformations is developed in this work to study equilibria of a beam under arbitrary end forces and moments. The slope angle and stretch strain of the centroid line, and shear strain of cross-sections, are chosen as dependent variables in this formulation, and end forces and moments can be either prescribed or resultant forces and moments due to constraints. Static equations of equilibria are derived from the principle of virtual work, which consist of one second-order ordinary differential equation and two algebraic equations. These equations are discretized using the finite difference method, and equilibria of the beam can be accurately calculated. For practical, geometrically nonlinear beam problems, stretch and shear strains are usually small, and a good approximate solution of the equations can be derived from the solution of the corresponding Euler–Bernoulli beam problem. The bending deformation of the beam is the only important one in a slender beam, and stretch and shear strains can be derived from it, which give a theoretical validation of the accuracy and applicability of the nonlinear Euler–Bernoulli beam formulation. Relations between end forces and moments and relative displacements of two ends of the beam can be easily calculated. This formulation is powerful in the study of buckling of beams with various boundary conditions under compression, and can be used to calculate post-buckling equilibria of beams. Higher-order buckling modes of a long slender beam that have complex configurations are also studied using this formulation.     

含电学边界的压电层合梁的非线性弯曲波

非线性科学己成为近代科学发展的一个重要标志,特别是非线性动力学和非线性波的研究对于解决自然科学各领域中遇到的复杂现象和问题有着极其重要的意义.本文研究了含电学边界条件的压电层合梁的非线性弯曲波传播特性.首先,考虑几何非线性效应和压电耦合效应,利用哈密顿原理建立了一维无限长矩形压电层合梁弯曲波的非线性方程.其次,采用Ja...     

广义边界梁的动力学建模与动态特征

对于广义边界条件Euler-Bernoulli梁,采用相对描述方式建立了可描述梁整体运动和相对变形的几何非线性及其线性化动力学模型,应用线性变换得到了该类梁的线性经典动力学方程,得到了广义边界条件下梁的横向振动代数特征方程、特征函数及特征值的退化表达式.算例分析了边界小扰动对固支-固支梁横向振动特征的影响规律.     

梁中非线性弯曲波传播特性的研究

研究了梁中的非线性弯曲波的传播特性,同时考虑了梁的大挠度引起的几何非线性效应和 梁的转动惯性导致的弥散效应,利用Hamilton变分法建立了梁中非线性弯曲波的波动方程. 对该方程进行了定性分析,在不同的条件下,该方程在相平面上存在同宿轨道或异宿轨道, 分别对应于方程的孤波解或冲击波解. 利用Jacobi椭圆函数展开法,对该非线性方程进行 求解,得到了非线性波动方程的准确周期解及相对应的孤波解和冲击波解,讨论了这些解存 在的必要条件,这与定性分析的结果完全相同. 利用约化摄动法从非线性弯曲波动方程中导 出了非线性Schr\"{o}dinger方程,从理论上证明了考虑梁的大挠度和转动惯性时梁中存在 包络孤立波.     

Geometrically exact planar beams with initial pre-stress and large curvature: Static configurations,natural frequencies,and mode shapes

Within this paper, an analytical formulation is provided and used to determine the natural frequencies and mode shapes of a planar beam with initial pre-stress and large variable curvature. The static configuration, mode shapes, and natural frequencies of the pre-stressed beam are obtained by using geometrically exact, Euler–Bernoulli beam theory. The beam is assumed to be not shear deformable and inextensible because of its slenderness and uniform, closed cross-section, as well as the boundary conditions under consideration. The static configuration and the modal information are validated with experimental data and compared to results obtained from nonlinear finite-element analysis software. In addition to the modal analysis about general static configurations, special consideration is given to an initially straight beam that is deformed into semi-circular and circular static configurations. For these special circular cases, the partial differential equation of motion is reduced to a sixth-order differential equation with constant coefficients, and solutions of this system are examined. This work can serve as a basis for studying slender structures with large curvatures.     

General analytic solution of dynamic response of beams with nonhomogeneity and variable cross-section

In this paper, a new method, the step-reduction method, is proposed to investigate the dynamic response of the Bernoulli-Euler beams with arbitrary nonhomogeneity and arbitrary variable cross-section under arbitrary loads. Both free vibration and forced vibration of such beams are studied. The new method requires to discretize the space domain into a number of elements. Each element can be treated as a homogeneous one with uniform thickness. Therefore, the general analytical solution of homogeneous beams with uniform cross-section can be used in each element. Then, the general analytic solution of the whole beam in terms of initial parameters can be obtained by satisfying the physical and geometric continuity conditions at the adjacent elements. In the case of free vibration, the frequency equation in analytic form can be obtained, and in the case of forced vibration, a final solution in analytical form can also be obtained which is involved in solving a set of simultaneous algebraic equations with only     

Geometric nonlinear formulation for curved beams with varying curvature

Based on exact Green strain of spatial curved beam, the nonlinear strain-displacement relation for plane curved beam with varying curvature is derived. Instead of using the previous straight beam elements, curved beam elements are used to approximate the curved beam with varying curvature. Based on virtual work principle, rigid-flexible coupling dynamic equations are obtained. Physical experiments were carried out to capture the large overall motion and the strain of curved beam to verify the present rigid-flexible coupling formulation for curved beam based on curved beam element. Numerical results obtained from simulations were compared with those results from the physical experiments. In order to illustrate the effectiveness of the curved beam element methodology, the simulation results of present curved beam elements are compared with those obtained by previous straight beam elements. The dynamic behavior of a slider-crank mechanism with an initially curved elastic connecting rod is investigated. The advantage of employing generalized-α method is pointed out and the special nonlinear dynamic characteristics of the curved beam are concluded.     

Effect of surface stress and surface-induced stress on behavior of piezoelectric nanobeam

A new continuum model is developed to study the influence of surface stress on the behaviors of piezoelectric nanobeams. Different from existing piezoelectric surface models which only consider the surface properties, the proposed model takes surfaceinduced initial fields into consideration. Due to the fact that the surface-induced initial fields are totally different under various boundary conditions, two kinds of beams, the doubly-clamped beam and the cantilever beam, are analyzed. Furthermore, boundary conditions can affect not only the initial state of the piezoelectric nanobeam but also the forms of the governing equations. Based on the Euler-Bernoulli beam theory, the nonlinear Green-Lagrangian strain-displacement relationship is applied. In addition, the surface area change is also considered in the proposed model. The governing equations of the doubly-clamped and cantilever beams are derived by the energy variation principle. Compared with existing Young-Laplace models, the proposed model for the doubly-clamped beam is similar to the Young-Laplace models. However, the governing equation of the cantilever beam derived by the proposed model is very different from that derived by the Young-Laplace models. The behaviors of piezoelectric nanobeams predicted by these two models also have significant discrepancies, which is owing to the surface-induced initial fields in the bulk beam.     

附磁阶梯变厚度悬臂梁压电俘能器的理论建模及分析

论文建立了一种附磁阶梯变厚度压电悬臂梁的动力学模型并分析了系统的俘能特性。基于Euler-Bernoulli梁理论分段建立系统能量函数并引入非线性磁势能,利用Lagrange方程建立了系统机电耦合动力学方程;利用数值方法分析了磁间距对系统振动特性的影响,此外还研究了系统单稳态和双稳态响应,探讨了厚度比、长度比、磁间距和外激励幅值对系统动力学响应和俘能特性的影响。结果表明,磁间距是影响系统势能的主要因素,调节磁间距可使系统产生单稳态和双稳态响应,从而有效提高俘能器俘能特性;与传统等截面悬臂梁压电俘能器相比,通过优化结构参数,附磁阶梯变厚度悬臂梁压电俘能器能够发生明显的非线性振动现象,实现宽频带振动能量采集。     

柔性梁大变形计算的一种新方法

研究了柔性梁大变形问题。常规Lagrangian有限元格式在处理大变形问题时,由于其单元插值函数不满足位移场的协调性要求,从而需要划分较多的单元,才能得到较好的结果。本文首先推导了Lagrangian坐标描述下的位移场变量满足的协调关系式,利用此关系式给出了位移场协调的非线性单元插值函数。基于虚功原理导出了梁大变形问题的非线性控制方程,数值计算结果证明了本文方法的正确性和有效性。     

Random excitation of coupled oscillators with single and simultaneous internal resonances

The primary objective of this paper is to examine the random response characteristics of coupled nonlinear oscillators in the presence of single and simultaneous internal resonances. A model of two coupled beams with nonlinear inertia interaction is considered. The primary beam is directly excited by a random support motion, while the coupled beam is indirectly excited through autoparametric coupling and parametric excitation. For a single one-to-two internal resonance, we used Gaussian and non-Gaussian closures, Monte Carlo simulation, and experimental testing to predict and measure response statistics and stochastic bifurcation in the mean square. The mean square stability boundaries of the coupled beam equilibrium position are obtained by a Gaussian closure scheme. The stochastic bifurcation of the coupled beam is predicted theoretically and experimentally. The stochastic bifurcation predicted by non-Gaussian closure is found to take place at a lower excitation level than the one predicted by Gaussian closure and Monte Carlo simulation. It is also found that above a certain excitation level, the solution obtained by non-Gaussian closure reveals numerical instability at much lower excitation levels than those obtained by Gaussian and Monte Carlo approaches. The experimental observations reveal that the coupled beam does not reach a stationary state, as reflected by the time evolution of the mean square response. For the case of simultaneous internal resonances, both Gaussian and non-Gaussian closures fail to predict useful results, and attention is focused on Monte Carlo simulation and experimental testing. The effects of nonlinear coupling parameters, internal detuning ratios, and excitation spectral density level are considered in both investigations. It is found that both studies reveal common nonlinear features such as bifurcations in the mean square responses of the coupled beam and modal interaction in the neighborhood of internal resonances. Furthermore, there is an upper limit for the excitation level above which the system experiences unbounded response in the neighborhood of simultaneous internal resonances.     

Nonlinear dynamics of an inclined beam subjected to a moving load

In this paper, the nonlinear dynamic response of an inclined pinned-pinned beam with a constant cross section, finite length subjected to a concentrated vertical force traveling with a constant velocity is investigated. The study is focused on the mode summation method and also on frequency analysis of the governing PDEs equations of motion. Furthermore, the steady-state response is studied by applying the multiple scales method. The nonlinear response of the beam is obtained by solving two coupled nonlinear PDEs governing equations of planar motion for both longitudinal and transverse oscillations of the beam. The dynamic magnification factor and normalized time histories of mid-pint of the beam are obtained for various load velocity ratios and the outcome results have been illustrated and compared to the results with those obtained from traditional linear solution. The appropriate parametric study considering the effects of the linear viscous damping, the velocity of the traveling load, beam inclination angle under zero or nonzero axial load are carried out to capture the influence of the effect of large deflections caused by stretching effects due to the beam’s immovable ends. It was seen that quadratic nonlinearity renders the softening effect on the dynamic response of the beam under the act of traveling load. Also in the case where the object leaves the inclined beam, its planar motion path is derived and the targeting accuracy is investigated and compared with those from the rigid solution assumption. Moreover, the stability analysis of steady-state response for the modes equations having quadratic nonlinearity was carried out and it was observed from the frequency response curves that for the considered parameters in the case of internal-external primary resonance, both saturation phenomenon and jump phenomenon can be predicted for the longitudinal excitation.     

双层纤维压电智能薄板几何非线性建模与分析

纤维压电MFC(Micro-Fiber Composite)的强致动力和高柔性等特点具有广泛的应用前景,但是材料组成结构复杂给建模带来了难处。基于Reissner-Mindlin假设,采用冯卡门非线性、中等转角及大转角几何非线性等理论,建立了MFC压电智能结构的多种几何非线性有限元模型。同时该模型考虑了压电纤维角度变化对结构形变的影响。分别对两种不同结构的MFC进行了建模与仿真,分别是MFC-d31和MFC-d33,前一种主要利用压电d31效应,而后一种主要利用压电d33效应。随后,通过一压电悬臂梁结构实验数据验证了模型的准确性。最后,利用所建模型对一种双层纤维压电智能薄板结构进行了几何非线性的计算与仿真。     

On the limitations of linear beams for the problems of moving mass-beam interaction using a meshfree method

This paper deals with the capabilities of linear and nonlinear beam theories in predicting the dynamic response of an elastically supported thin beam traversed by a moving mass. To this end, the discrete equations of motion are developed based on Lagrange’s equations via reproducing kernel particle method (RKPM). For a particular case of a simply supported beam, Galerkin method is also employed to verify the results obtained by RKPM, and a reasonably good agreement is achieved. Variations of the maximum dynamic deflection and bending moment associated with the linear and nonlinear beam theories are investigated in terms of moving mass weight and velocity for various beam boundary conditions. It is demonstrated that for majority of the moving mass velocities, the differences between the results of linear and nonlinear analyses become remarkable as the moving mass weight increases, particularly for high levels of moving mass velocity. Except for the cantilever beam, the nonlinear beam theory predicts higher possibility of moving mass separation from the base beam compared to the linear one. Furthermore, the accuracy levels of the linear beam theory are determined for thin beams under large deflections and small rotations as a function of moving mass weight and velocity in various boundary conditions.