Out-of-plane responses of a circular curved Timoshenko beam due to a moving load

For the cases of using the finite curved beam elements and taking the effects of both the shear deformation and rotary inertias into consideration, the literature regarding either free or forced vibration analysis of the curved beams is rare. Thus, this paper tries to determine the dynamic responses of a circular curved Timoshenko beam due to a moving load using the curved beam elements. By taking account of the effect of shear deformation and that of rotary inertias due to bending and torsional vibrations, the stiffness matrix and the mass matrix of the curved beam element were obtained from the force–displacement relations and the kinetic energy equations, respectively. Since all the element property matrices for the curved beam element are derived based on the local polar coordinate system (rather than the local Cartesian one), their coefficients are invariant for any curved beam element with constant radius of curvature and subtended angle and one does not need to transform the property matrices of each curved beam element from the local coordinate system to the global one to achieve the overall property matrices for the entire curved beam structure before they are assembled. The availability of the presented approach has been verified by both the existing analytical solutions for the entire continuum curved beam and the numerical solutions for the entire discretized curved beam composed of the conventional straight beam elements based on either the consistent-mass model or the lumped-mass model. In addition to the typical circular curved beams, a hybrid curved beam composed of one curved-beam segment and two identical straight-beam segments subjected to a moving load was also studied. Influence on the dynamic responses of the curved beams of the slenderness ratio, moving-load speed, shear deformation and rotary inertias was investigated.     

计及剪切影响的RC梁的非线性分析

为提高钢筋混凝土RC梁的计算效率和精度,提出了一种基于梁截面弯矩-曲率关系的宏观有限元方法,可用于各种跨高比RC梁的材料非线性分析。首先假定了混凝土和钢筋的非线性应力-应变关系,然后引入经过修正的Rodriguez截面模型,根据边界顶点把截面划分成若干梯形单元,利用quasi-Newton法求解由两个变量耦合而成的截面非线性平衡方程,由此建立RC截面的弯矩-曲率关系。在此基础上利用Timoshenko梁弯曲理论建立考虑横向剪切变形影响的RC梁的有限元分析模型。通过对试验梁的分析对比验证了所提出的分析方法的适用性。     

Refined Series Methodology for the Fully Coupled Thin-Walled Laminated Beams Considering Foundation Effects

The refined power series solutions are presented for the coupled static analysis of thin-walled laminated beams resting on elastic foundation. For this purpose, the elastic strain energy considering the material and structural coupling effects and the energy including the foundation effects are constructed. The equilibrium equations and the force-displacement relationships are derived from the extended Hamilton's principle, and the explicit expressions for displacement parameters are presented based on power series expansions of displacement components. Finally, the member stiffness matrix is determined by using the force-displacement relationships. For comparison, the finite element model based on the Hermite cubic interpolation polynomial is presented. In order to verify the accuracy and the superiority of the laminated beam element developed by this study, the numerical solutions are presented and compared with results obtained from the regular finite beam elements and the ABAQUS's shell elements. The influences of the fiber angle change and the boundary conditions on the coupled behavior of laminated beams with mono-symmetric I-sections are investigated.     

压电层合梁静动力学分析的高精度梁单元

给出了一个对复合材料压电层合梁进行数值分析的高精度压电层合梁单元。基于Shi三阶剪切变形板理论的位移场和Layer-wise理论的电势场,用力-电耦合的变分原理及Hamilton原理推导了压电层合梁单元列式。采用拟协调元方法推导了一个可显式给出单元刚度矩阵的两节点压电层合梁单元,并应用于压电层合梁的力-电耦合弯曲和自由振动分析。计算结果表明,该梁单元给出的梁挠度和固有频率与解析解吻合良好,并优于其它梁单元的计算结果,说明了本文所给压电层合梁单元的可靠性和准确性。研究结果可为力-电耦合作用下压电层合梁的力学分析提供一个简单、精确且高效的压电层合梁单元。     

Free vibration characteristics of sectioned unidirectional/bidirectional functionally graded material cantilever beams based on finite element analysis

Advancements in manufacturing technology, including the rapid development of additive manufacturing (AM), allow the fabrication of complex functionally graded material (FGM) sectioned beams. Portions of these beams may be made from different materials with possibly different gradients of material properties. The present work proposes models to investigate the free vibration of FGM sectioned beams based on onedimensional (1D) finite element analysis. For this purpose, a sample beam is divided into discrete elements, and the total energy stored in each element during vibration is computed by considering either the Timoshenko or Euler-Bernoulli beam theory. Then, Hamilton's principle is used to derive the equations of motion for the beam. The effects of material properties and dimensions of FGM sections on the beam's natural frequencies and their corresponding mode shapes are then investigated based on a dynamic Timoshenko model (TM). The presented model is validated by comparison with three-dimensional (3D) finite element simulations of the first three mode shapes of the beam.     

Bending of beams on three-parameter elastic foundation

The bending of a Timoshenko beam resting on a Kerr-type three-parameter elastic foundation is introduced, its governing differential equations are formulated and analytically solved, and the solutions are discussed and applied to particular problems. Parametric analyses of elastically supported beams of infinite and finite length are carried out and comparisons are made between one, two or three-parameter foundation models and more accurate 2D finite element models. In order to estimate the necessary soil parameters, an analytical procedure based on the modified Vlasov model is proposed. The presented solutions and applications show the superiority of the Kerr-type foundation model compared to one or two-parameter models.     

Axial-flexural coupled vibration and buckling of composite beams using sinusoidal shear deformation theory

A finite element model based on sinusoidal shear deformation theory is developed to study vibration and buckling analysis of composite beams with arbitrary lay-ups. This theory satisfies the zero traction boundary conditions on the top and bottom surfaces of beam without using shear correction factors. Besides, it has strong similarity with Euler–Bernoulli beam theory in some aspects such as governing equations, boundary conditions, and stress resultant expressions. By using Hamilton’s principle, governing equations of motion are derived. A displacement-based one-dimensional finite element model is developed to solve the problem. Numerical results for cross-ply and angle-ply composite beams are obtained as special cases and are compared with other solutions available in the literature. A variety of parametric studies are conducted to demonstrate the effect of fiber orientation and modulus ratio on the natural frequencies, critical buckling loads, and load-frequency curves as well as corresponding mode shapes of composite beams.     

Nonlinear Vibration of Plane Structures by Finite Element and Incremental Harmonic Balance Method

A nonlinear steady state vibration analysis of a wide class of planestructures is analyzed. Both the finite element method and incrementalharmonic balance method are used. The usual beam element is adopted inwhich the nonlinear effect arising from longitudinal stretching has beentaken into account. Based on the geometric nonlinear finite elementanalysis, the nonlinear dynamic equations including quadratic and cubicnonlinearities are derived. These equations are solved by theincremental harmonic balance (IHB) method. To show the effectiveness andversatility of this method, some typical examples for a wide variety ofvibration problems including fundamental resonance, super- andsub-harmonic resonance, and combination resonance of plane structuressuch as beams, shallow arches and frames are computed. Most of theseexamples have not been studied by other researchers before. Comparisonwith previous results are also made.     

Non-local finite element analysis of damped beams

Non-local viscoelastic beam models are used to analyse the dynamics of beams with different boundary conditions using the finite element method. Unlike local damping models the internal force of the non-local model is obtained as weighted average of state variables over a spatial domain via convolution integrals with spatial kernel functions that depend on a distance measure. In the finite element analysis, the interpolating shape functions of the element displacement field are identical to those of standard two-node beam elements. However, for non-local damping, nodes remote from the element do have an effect on the energy expressions, and hence on the damping matrix. The expressions of these direct and cross damping matrices may be obtained explicitly for some common spatial kernel functions and Euler–Bernoulli beam theory. Alternatively numerical integration may be applied to obtain solutions. Examples are given where the eigenvalues are compared to the exact solution for a pinned–pinned beam to demonstrate the convergence of the finite element method. The results for beams with other boundary conditions are used to demonstrate the versatility of the finite element technique.     

基于初应力法的作大范围运动柔性梁动力学建模理论研究

研究了初应力法的作大范围运动柔性梁的建模理论.根据连续介质理论,考虑应变-位移中的非线性项,用一致质量有限元法对柔性梁进行离散,基于Jourdain速度变分原理导出定轴转动下大范围运动为自由的柔性梁刚-柔耦合动力学方程.从其刚柔耦合动力学方程出发,考虑在大范围运动已知情况下的结构动力学方程.通过引入准静态概念,把其结构动力学方程转化为准静态方程.对纵向和横向变形节点坐标进行坐标分离,解出与纵向变形相关的准静态方程,得到准静态时的纵向应力表达式,从而获得附加刚度项.并对此非惯性系下作大范围运动柔性梁的结构动力学方程进行数值仿真,对零次近似模型、一次近似模型、初应力法动力学模型的仿真结果进行分析,揭示三种模型的动力学性质的差异.     

基于能量平衡原理的有限单元耦合方法

结构非线性数值计算分析应真实反映局部损伤破坏细节,以作为损伤演化全过程分析的依据。对同类构件,有限单元耦合方法可以解决破坏细节与整体模拟的空间尺度差异问题。基于能量平衡原理,建立了梁与实体单元、梁与壳单元以及壳与实体单元的耦合方程,适用于结构的损伤数值计算。对某RC框架结构原位推覆试验的损伤数值分析表明,有限单元耦合模型能正确反映整体结构的承载力和变形性能,并且能准确反映局部损伤破坏细节。     

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

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

Multi-step transversal and tangential linearization methods applied to a class of nonlinear beam equations

A novel and continuously parameterized form of multi-step transversal linearization (MTrL) method is developed and numerically explored for solving nonlinear ordinary differential equations governing a class of boundary value problems (BVPs) of relevance in structural mechanics. A similar family of multi-step tangential linearization (MTnL) methods is also developed and applied to such BVP-s. Within the framework of MTrL and MTnL, a BVP is treated as a constrained dynamical system, i.e. a constrained initial value problem (IVP). While the MTrL requires the linearized solution manifold to transversally intersect the nonlinear solution manifold at a chosen set of points across the axis of the independent variable, the essential difference of the present MTrL method from its previous version [Roy, D., Kumar, R., 2005. A multistep transversal linearization (MTL) method in nonlinear structural dynamics. J. Sound Vib. 17, 829–852.] is that it has the flexibility of treating nonlinear damping and stiffness terms as time-variant damping and stiffness terms in the linearized system. The resulting time-variant linearized system is then solved using Magnus’ characterization [Magnus, W., 1954. On the exponential solution of differential equations for a linear operator. Commun. Pure Appl. Math., 7, 649–673.]. Towards numerical illustrations, response of a tip loaded cantilever beam (Elastica) is first obtained. Next, the response of a simply supported nonlinear Timoshenko beam is obtained using a variationally correct (VC) model for the beam [Marur, S., Prathap, G., 2005. Nonlinear beam vibration problems and simplification in finite element model. Comput. Mech. 35(5), 352–360.]. The new model does not involve any simplifications commonly employed in the finite element formulations in order to ease the computation of nonlinear stiffness terms from nonlinear strain energy terms. A comparison of results through MTrL and MTnL techniques consistently indicate a superior quality of approximations via the transversal linearization technique. While the usage of tangential system matrices is common in nonlinear finite element practices, it is demonstrated that the transversal version of linearization offers an easier and more general implementation, requires no computations of directional derivatives and leads to a consistently higher level of numerical accuracy. It is also observed that higher order versions of MTrL/MTnL with Lagrangian interpolations may not work satisfactorily and hence spline interpolations are suggested to overcome this problem.     

Analytical models for the inflation of a polymeric tube

In the blow moulding processes, a tube-shaped parison (also called preform) made of polymer is inflated inside a mould in order to obtain the desired bottle shape. The free inflation of a cylindrical tube (a simplified kinematic of inflation for the preform before contacting the mould) has been investigated in order to develop analytical models and then provide reliable validation material for finite element software. The solutions calculated for a tube made of a Newtonian fluid and a Lodge's rubberlike liquid material using both exact volumic approach and thin shell approximation have been analysed and compared.     

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

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

Euler–Bernoulli beams with multiple singularities in the flexural stiffness

Euler–Bernoulli beams under static loads in presence of discontinuities in the curvature and in the slope functions are the object of this study. Both types of discontinuities are modelled as singularities, superimposed to a uniform flexural stiffness, by making use of distributions such as unit step and Dirac's delta functions. A non-trivial generalisation to multiple different singularities of an integration procedure recently proposed by the authors for a single singularity is presented in this paper. The proposed integration procedure leads to closed form solutions, dependent on boundary conditions only, which do not require enforcement of continuity conditions along the beam span. It is however shown how, from the solution of the clamped-clamped beam, by considering suitable singularities at boundaries in the flexural stiffness model, responses concerning several boundary conditions can be recovered. Furthermore, solutions in terms of deflection of the beam are obtained for imposed displacements at boundaries providing the so called shape functions. The above mentioned shape functions can be adopted to insert beams with singularities as frame elements in a finite element discretisation of a frame structure. Explicit expressions of the element stiffness matrix are provided for beam elements with multiple singularities and the reduction of degrees of freedom with respect to classical finite element procedures is shown. Extension of the proposed procedure to beams with axial displacement and vertical deflection discontinuities is also presented.     

Nonlinear stability and vibration of pre/post-buckled microstructure-dependent FGPM actuators

The present study deals with the study of the nonlinear stability and small free vibration of microstructure-dependent functionally graded piezoelectric material (FGPM) beams in pre/post-buckling regimes. The Timoshenko beam theory with various inplane and out-of-plane boundary conditions are considered under different types of mechanical and thermal loads. The beam is assumed to be under inplane mechanical, thermal, and electrical excitations. Each thermo-electro-mechanical property of the beam is graded across the thickness (i.e., height) of the beam, based on a power law model. The von Kármán type geometric nonlinearity is included to account for the large deflection behavior of the beam under inplane loads. The modified couple stress theory is included to account for the size effects. A weak-form, displacement-based, finite element formulation is developed to discretize the equations of motion. The resulting system of nonlinear algebraic equations is solved using Newton’s iterative method. The numerical results of frequencies and lateral deflections as a function of load parameters reveal the existence of bifurcation or critical states in some cases. The effects of load type, microstructural dependency, boundary conditions, beam geometry, composition rule of the constituents, and actuator voltage are discussed through the various parametric studies.     

Higher-order impact modeling of sandwich structures with flexible core

In this study, a higher-order impact model is presented to simulate the response of a soft-core sandwich beam subjected to a foreign object impact. A free vibration problem of sandwich beams is first solved, and the results are validated by comparing with numerical finite element modeling results of ABAQUS and the solution by Frostig and Baruch [Frostig, Y., Baruch, M., 1994. Free vibration of sandwich beams with a transversely flexible core: a high order approach. Journal of Sound and Vibration 176(2), 195–208]. Then a foreign object impact process is incorporated in the higher-order model, and the contact force and deflection history as well as the propagation of transverse normal, shear, and axial stresses during the impact are analyzed and discussed. The validity of the model in the impact response predictions is demonstrated by comparing with finite element solutions of LS-DYNA. The calculated stresses caused by a foreign object impact are then used to assess failure locations, failure time, and failure modes in sandwich beams, which are shown to compare well with the available experimental results. The effects of impact mass, initial velocity, core stiffness, and core height on the impact stresses generated in the beams are discussed. The influences of impact mass and initial velocity on the contact force history are close to those by the linearized impact solution, but the proposed higher-order impact model captures the non-linear impact process and different generated stresses. Compared to the fully backed sandwich case, the core height shows a great influence over the impact process of a simply supported sandwich system, in which the global behavior of the sandwich is dominant; while the core stiffness shows minor effect over the impact process. The higher-order impact model of sandwich beams developed in the study provides accurate predictions of the generated stresses and impact process and can be used effectively in design analysis of anti-impact structures made of sandwich materials.     

Differential constitutive equations of incompressible media with finite deformations

A method is proposed for constructing a system of constitutive equations of an incompressible medium with nonlinear dissipative properties with finite deformations. A scheme of the mechanical behavior of a material is used, in which the points are connected by horizontally aligned elastic, viscous, plastic, and transmission elements. The properties of each element of the scheme are described with the use of known equations of the nonlinear elasticity theory, the theory of nonlinear viscous fluids, and the theory of plastic flow of the material under conditions of finite deformations of the medium. The system of constitutive equations is closed by equations that express the relation between the deformation rate tensor of the material and the deformation rate tensor of the plastic element. Transmission elements are used to take into account a significant difference between macroscopic deformations of the material and deformations of elements of the medium at the structural level. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 3, pp. 158–170, May–June, 2009.     

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.