非保守功能梯度弹性组合曲梁的精确模型及其数值解

基于直法线假设,采用可伸长梁的几何非线性理论,建立了功能梯度材料弹性组合曲梁受切线均布随从力作用下的静态大变形数学模型。该模型不仅计及了轴线伸长,同时也精确地考虑了梁的初始曲率对变形的影响以及轴向变形与弯曲变形之间的耦合效应。用打靶法数值求解了由金属和陶瓷两相材料所构成的一种FGM组合曲梁在沿轴线均布切向随动载荷作用下的非线性平面弯曲问题,给出了不同梯度指标下FGM弹性曲梁随载荷参数大范围变化的平衡路径,并与金属和陶瓷两种单相材料曲梁的相应特性进行了比较。     

Direct method for analysis of flexible cantilever beam subjected to two follower forces

The static analysis of the flexible non-uniform cantilever beams under a tip-concentrated and intermediate follower forces is considered. The angles of inclination of the concentrated forces with respect to the deformed axis of the beam remain unchanged during deformation. The governing non-linear boundary-value problem is reduced to an initial-value problem by change of variables. The resulting problem can be solved without iterations. It is shown that there are no critical loads in the Euler sense (divergence) for any flexural-stiffness distribution and angles of inclination of the follower forces. In particular, if the follower forces are tangential, the rectilinear shape of the non-uniform cantilever beam is the only possible equilibrium configuration. In this paper some equilibrium configurations of the uniform cantilever under normal or tangential follower forces are presented using direct method.     

弹性曲梁在机械和热载荷共同作用下的几何非线性模型及其数值解

基于Euler-Bernoulli梁的几何非线性理论,建立了弹性曲梁在任意分布机械载荷和热载荷共同作用下的几何非线性静平衡控制方程。该模型不仅计及了轴线伸长,同时也精确地考虑了梁的初始曲率对变形的影响以及轴向变形与弯曲变形之间的相互耦合效应。应用打靶法数值求解了半圆形曲梁在横向均匀升温作用下的非线性弯曲问题,数值比较了轴向伸长对曲梁变形的影响。     

Flexural-torsional bifurcations of a cantilever beam under potential and circulatory forces I: Non-linear model and stability analysis

The stability of a cantilever elastic beam with rectangular cross-section under the action of a follower tangential force and a bending conservative couple at the free end is analyzed. The beam is herein modeled as a non-linear Cosserat rod model. Non-linear, partial integro-differential equations of motion are derived expanded up to cubic terms in the transversal displacement and torsional angle of the beam. The linear stability of the trivial equilibrium is studied, revealing the existence of buckling, flutter and double-zero critical points. Interaction between conservative and non-conservative loads with respect to the stability problem is discussed. The critical spectral properties are derived and the corresponding critical eigenspace is evaluated.     

Series solution for large deflections of a cantilever beam with variable flexural rigidity

In this paper large deflection and rotation of a nonlinear Bernoulli-Euler beam with variable flexural rigidity and subjected to a static co-planar follower loading is studied. It is assumed that the angle of inclination of the force with respect to the deformed axis of the beam remains unchanged during deformation. The governing equation of this problem is solved analytically for the first time using a new kind of analytical technique for nonlinear problems, namely the Homotopy Analysis Method (HAM). The present solution can be used for the analysis of a wide range of loads, material/cross section properties and lengths for beams undergoing large deformations. The results obtained from HAM are compared with results reported in previous works. Finally, the load–displacement characteristics of a uniform cantilever beam with different material properties under a follower force applied normal to the deformed beam axis are presented.     

Elasticity solution of clamped-simply supported beams with variable thickness

This paper studies the stress and displacement distributions of continuously varying thickness beams with one end clamped and the other end simply supported under static loads. By introducing the unit pulse functions and Dirac functions, the clamped edge can be made equivalent to the simply supported one by adding the unknown horizontal reactions. According to the governing equations of the plane stress problem, the general expressions of displacements, which satisfy the governing differefitial equations and the boundary conditions attwo ends of the beam, can be deduced. The unknown coefficients in the general expressions are then determined by using Fourier sinusoidal series expansion along the upper and lower boundaries of the beams and using the condition of zero displacements at the clamped edge. The solution obtained has excellent convergence properties. Comparing the numerical results to those obtained from the commercial software ANSYS, excellent accuracy of the present method is demonstrated.     

Dynamic stability and response of fluttered beams subjected to random follower forces

This paper presents a study of non-linear response of a fluttered, cantilevered beam subjected to a random follower force at the free end. The random follower force is characterized as the sum of a post-critical static force and a stationary process with a zero mean. First, the Ritz-Galerkin method is applied to yield a set of discretized system equations. The system equations are then partially uncoupled by a special modal analysis based on normal modes of the corresponding linear, autonomous system at the onset of fluttering. Next, the stochastic averaging method is utilized to get Ito's differential equation governing the amplitude of the fluttered mode. Finally, the probability density function for the amplitude of the fluttered mode is obtained by solving the FPK equation. Numerical results show that the probability density function for the amplitude of the fluttered mode is determined by the sample behavior of the beam near the trivial equilibrium configuration.     

Three-dimensional equilibria of nonlinear pre-curved beams using an intrinsic formulation and shooting

This article describes a shooting method for computing three-dimensional equilibria of pre-curved nonlinear beams with axial and shear flexibility using the intrinsic beam formulation. For distributed and concentrated follower loads acting on a cantilevered beam, the method amounts to a direct solution approach requiring only a single shot (zero iterations) to compute the equilibria. This is possible since the system equations are defined in a local coordinate system that rotates and translates with the beam, akin to the follower loads themselves. A general procedure employing nonconservative follower loads, which invokes the Picard–Lindelöf theorem on uniqueness and existence of solutions, is also introduced for finding all solutions for three-dimensional pre-curved beam problems with conservative loading. This is particularly useful in beam buckling problems where multiple stable and unstable solutions exist. Three-dimensional equilibrium solutions are generated for many loading cases and boundary conditions, including three-dimensional helical beams, and are compared to similar solutions where available in the literature. Excellent agreement is documented in all comparison cases. For buckling examples, the stability of the computed solutions is assessed using a dynamic finite element code based on the same intrinsic beam equations. Due to the ability to avoid iteration, the presented approach may find application in model-based control for practical three-dimensional problems such as the control of manipulators utilized in endoscopic surgeries and the control of spacecraft with robotic arms and long cables.     

Free vibration of functionally graded material beams with surface-bonded piezoelectric layers in thermal environment

Free vibration of statically thermal postbuckled functionally graded material （FGM） beams with surface-bonded piezoelectric layers subject to both temperature rise and voltage is studied. By accurately considering the axial extension and based on the Euler-Bernoulli beam theory, geometrically nonlinear dynamic governing equations for FGM beams with surface-bonded piezoelectric layers subject to thermo-electro- mechanical loadings are formulated. It is assumed that the material properties of the middle FGM layer vary continuously as a power law function of the thickness coordinate, and the piezoelectric layers are isotropic and homogenous. By assuming that the amplitude of the beam vibration is small and its response is harmonic, the above mentioned non-linear partial differential equations are reduced to two sets of coupled ordinary differential equations. One is for the postbuckling, and the other is for the linear vibration of the beam superimposed upon the postbuckled configuration. Using a shooting method to solve the two sets of ordinary differential equations with fixed-fixed boundary conditions numerically, the response of postbuckling and free vibration in the vicinity of the postbuckled configuration of the beam with fixed-fixed ends and subject to transversely nonuniform heating and uniform electric field is obtained. Thermo-electric postbuckling equilibrium paths and characteristic curves of the first three natural frequencies versus the temperature, the electricity, and the material gradient parameters are plotted. It is found that the three lowest frequencies of the prebuckled beam decrease with the increase of the temperature, but those of a buckled beam increase monotonically with the temperature rise. The results also show that the tensional force produced in the piezoelectric layers by the voltage can efficiently increase the critical buckling temperature and the natural frequency.     

Equilibria of axially moving beams in the supercritical regime

Equilibria of axially moving beams are computationally investigated in the supercritical transport speed ranges. In the supercritical regime, the pattern of equilibria consists of the straight configuration and of non-trivial solutions that bifurcate with transport speed. The governing equations of coupled planar is reduced to a partial-differential equation and an integro-partial-differential equation of transverse vibration. The numerical schemes are respectively presented for the governing equations and the corresponding static equilibrium equation of coupled planar and the two governing equations of transverse motion for non-trivial equilibrium solutions via the finite difference method and differential quadrature method under the simple support boundary. A steel beam is treated as example to demonstrate the non-trivial equilibrium solutions of three nonlinear equations. Numerical results indicate that the three models predict qualitatively the same tendencies of the equilibrium with the changing parameters and the integro-partial-differential equation yields results quantitatively closer to those of the coupled equations.     

Nonlinear vibrations of a sandwich piezo-beam system under piezoelectric actuation

The paper describes the use of active structures technology for deformation and nonlinear free vibrations control of a simply supported curved beam with upper and lower surface-bonded piezoelectric layers, when the curvature is a result of the electric field application. Each of the active layers behaves as a single actuator, but simultaneously the whole system may be treated as a piezoelectric composite bender. Controlled application of the voltage across piezoelectric layers leads to elongation of one layer and to shortening of another one, which results in the beam deflection. Both the Euler–Bernoulli and von Karman moderately large deformation theories are the basis for derivation of the nonlinear equations of motion. Approximate analytical solutions are found by using the Lindstedt–Poincaré method which belongs to perturbation techniques. The method makes possible to decompose the governing equations into a pair of differential equations for the static deflection and a set of differential equations for the transversal vibration of the beam. The static response of the system under the electric field is investigated initially. Then, the free vibrations of such deformed sandwich beams are studied to prove that statically pre-stressed beams have higher natural frequencies in regard to the straight ones and that this effect is stronger for the lower eigenfrequencies. The numerical analysis provides also a spectrum of the amplitude-dependent nonlinear frequencies and mode shapes for different geometrical configurations. It is demonstrated that the amplitude–frequency relation, which is of the hardening type for straight beams, may change from hard to soft for deformed beams, as it happens for the symmetric vibration modes. The hardening-type nonlinear behaviour is exhibited for the antisymmetric vibration modes, independently from the system stiffness and dimensions.

     

Linear and non-linear interactions between static and dynamic bifurcations of damped planar beams

The critical and post-critical behavior of a non-conservative non-linear structure, undergoing statical and dynamical bifurcations, is analyzed. The system consists of a purely flexible planar beam, equipped with a lumped visco-elastic device, loaded by a follower force. A unique integro-differential equation of motion in the transversal displacement, with relevant boundary conditions, is derived. Then, the linear stability diagram of the trivial rectilinear configuration is built-up in the parameter space. Particular emphasis is given to the role of the damping on the critical scenario. The occurrence of different mechanisms of instability is highlighted, namely, of divergence, Hopf, double zero, resonant and non-resonant double Hopf, and divergence-Hopf bifurcation. Attention is then focused on the two latter (codimension-two) bifurcations. A multiple scale analysis is carried-out directly on the continuous model, and the relevant non-linear bifurcation equations in the amplitudes of the two interactive modes are derived. The fixed-points of these equations are numerically evaluated as functions of two bifurcation parameters and some equilibrium paths illustrated. Finally, the bifurcation diagrams, illustrating the system behavior around the critical points of the parameter space, are obtained.     

Elastica of a variable-arc-length circular curved beam subjected to an end follower force

This paper presents postbuckling behaviors of a variable-arc-length (VAL) circular curved beam subjected to an end follower force. One end of the VAL circular curved beam is hinged while the other end is supported by a frictionless slot, which is fixed horizontally and vertically but is allowed to rotate corresponding to loading direction. When the VAL circular curved beam is deformed, the total arc-length of the circular curved beam varies. Two approaches have been applied for the solution of this problem. The first approach is an elliptic integrals method based on elastica theory, which yields the exact closed-form solution in terms of the first and second kinds of elliptic integrals. For validation of the results, the shooting method is employed for a numerical solution by developing the set of nonlinear governing differential equations together with boundary conditions, and then integrating them by using the fourth-order Runge–Kutta algorithm. The results from both approaches are in very good agreement. From the results, it is found that the VAL circular curved beam subjected to an end follower force can be deformed in many mode shapes. For the first and third modes, the beam exhibits both stable and unstable configurations, whereas for the second mode only an unstable configuration exists. The influences of initial curvature on the critical load and the deformed configurations are highlighted.     

FREE VIBRATION OF A SIMPLY SUPPORTED BEAM UNDER A NON-CONSERVATIVE DISTRIBUTED LOAD

对受非保守载荷的简支梁在后屈曲附近的自由振动进行了研究. 基于可伸长梁的大变形理论,建立了受沿轴线分布切向非保守力作用的简支梁后屈曲附近自由振动的几何非线性模型. 在小振幅和谐振动假设下,简化得到后屈曲梁线性振动的控制方程. 采用打靶法求解振动问题的控制方程,给出了前三阶固有频率与载荷之间的特征关系曲线. 结果表明：非保守载荷作用下梁的振动响应与保守载荷作用下梁的振动响应有着明显不同.     

Non-linear stability analysis of imperfect thin-walled composite beams

The static non-linear behavior of thin-walled composite beams is analyzed considering the effect of initial imperfections. A simple approach is used for determining the influence of imperfection on the buckling, prebuckling and postbuckling behavior of thin-walled composite beams. The fundamental and secondary equilibrium paths of perfect and imperfect systems corresponding to a major imperfection are analyzed for the case where the perfect system has a stable symmetric bifurcation point. A geometrically non-linear theory is formulated in the context of large displacements and rotations, through the adoption of a shear deformable displacement field. An initial displacement, either in vertical or horizontal plane, is considered in presence of initial geometric imperfection. Ritz's method is applied in order to discretize the non-linear differential system and the resultant algebraic equations are solved by means of an incremental Newton-Rapshon method. The numerical results are presented for a simply supported beam subjected to axial or lateral load. It is shown in the examples that a major imperfection reduces the load-carrying capacity of thin-walled beams. The influence of this effect is analyzed for different fiber orientation angle of a symmetric balanced lamination. In addition, the postbuckling response obtained with the present beam model is compared with the results obtained with a shell finite element model (Abaqus).     

A robust hyper-elastic beam model under bi-axial normal-shear loadings

In this paper, a simple and robust constitutive model is proposed to simulate mechanical behaviors of hyper-elastic materials under bi-axial normal-shear loadings in the finite strain regime. The Mooney–Rivlin strain energy function is adopted to develop a two-dimensional (2D) normal-shear constitutive model within the framework of continuum mechanics. A motion field is first proposed for combined normal and shear deformations. The deformation gradient of the proposed field is calculated and then substituted into right Cauchy–Green deformation tensor. Constitutive equations are then derived for normal and shear deformations. They are two explicit coupled equations with high-level polynomial non-linearity. In order to examine capabilities of the developed hyper-elastic model, uniaxial tensile responses and non-linear stability behaviors of moderately thick straight and curved beams undergoing normal axial and transverse shear deformations are simulated and compared with experiments. Fused deposition modeling technique as a 3D printing technology is implemented to fabricate hyper-elastic beam structures from soft poly-lactic acid filaments. The printed specimens are tested under tensile/compressive in-plane and compressive out-of-plane forces. A finite element formulation along with the Newton–Raphson and Riks techniques is also developed to trace non-linear equilibrium path of beam structures in large defamation regimes. It is shown that the model is capable of predicting non-linear equilibrium characteristics of hyper-elastic straight and curved beams. It is found that the modeling of shear deformation and finite strain is essential toward an accurate prediction of the non-linear equilibrium responses of moderately thick hyper-elastic beams. Due to simplicity and accuracy, the model can serve in the future studies dealing with the analysis of hyper-elastic structures in which two normal and shear stress components are dominant.     

A nonlinear mathematical model for large deflection of incompressible saturated poroelastic beams

Nonlinear governing equations are established for large deflection of incom- pressible fluid saturated poroelastic beams under constraint that diffusion of the pore fluid is only in the axial direction of the deformed beams.Then,the nonlinear bend- ing of a saturated poroelastic cantilever beam with fixed end impermeable and free end permeable,subjected to a suddenly applied constant concentrated transverse load at its free end,is examined with the Gaierkin truncation method.The curves of deflections and bending moments of the beam skeleton and the equivalent couples of the pore fluid pressure are shown in figures.The results of the large deflection and the small deflection theories of the cantilever poroelastic beam are compared,and the differences between them are revealed.It is shown that the results of the large deflection theory are less than those of the corresponding small deflection theory,and the times needed to approach its stationary states for the large deflection theory are much less than those of the small deflection theory.     

Nonlinear dynamic stability analysis of Euler–Bernoulli beam–columns with damping effects under thermal environment

In this study, a unified nonlinear dynamic buckling analysis for Euler–Bernoulli beam–columns subjected to constant loading rates is proposed with the incorporation of mercurial damping effects under thermal environment. Two generalized methods are developed which are competent to incorporate various beam geometries, material properties, boundary conditions, compression rates, and especially, the damping and thermal effects. The Galerkin–Force method is developed by implementing Galerkin method into force equilibrium equations. Then for solving differential equations, different buckled shape functions were introduced into force equilibrium equations in nonlinear dynamic buckling analysis. On the other hand, regarding the developed energy method, the governing partial differential equation for dynamic buckling of beams is also derived by meticulously implementing Hamilton’s principles into Lagrange’s equations. Consequently, the dynamic buckling analysis with damping effects under thermal environment can be adequately formulated as ordinary differential equations. The validity and accuracy of the results obtained by the two proposed methods are rigorously verified by the finite element method. Furthermore, comprehensive investigations on the structural dynamic buckling behavior in the presence of damping effects under thermal environment are conducted.