A computationally efficient non-linear beam model

Aerospace structures with large aspect ratio, such as airplane wings, rotorcraft blades, wind turbine blades, and jet engine fan and compressor blades, are particularly susceptible to aeroelastic phenomena. Finite element analysis provides an effective and generalized method to model these structures; however, it is computationally expensive. Fortunately, the large aspect ratio of these structures is exploitable as these potential aeroelastically unstable structures can be modeled as cantilevered beams, drastically reducing computational time.In this paper, the non-linear equations of motion are derived for an inextensional, non-uniform cantilevered beam with a straight elastic axis. Along the elastic axis, the cross-sectional center of mass can be offset in both dimensions, and the principal bending and centroidal axes can each be rotated uniquely. The Galerkin method is used, permitting arbitrary and abrupt variations along the length that require no knowledge of the spatial derivatives of the beam properties. Additionally, these equations consistently retain all third-order non-linearities that account for flexural-flexural-torsional coupling and extend the validity of the equations for large deformations.Furthermore, linearly independent shape functions are substituted into these equations, providing an efficient method to determine the natural frequencies and mode shapes of the beam and to solve for time-varying deformation.This method is validated using finite element analysis and is extended to swept wings. Finally, the importance of retaining cubic terms, in addition to quadratic terms, for non-linear analysis is demonstrated for several examples.     

带集中质量的旋转柔性曲梁动力学特性分析

本文对带集中质量的平面内旋转柔性曲梁动力学特性进行了研究.基于绝对节点坐标法推导出曲梁单元,其中该曲梁单元采用Green-Lagrangian应变,并根据曲梁变形前后的曲率变化和曲率的精确表达式计算了曲梁单元弹性力所作的虚功.通过虚功原理,利用$\delta$函数和中心刚体与悬臂曲梁之间的固支边界条件,建立了带集中质量的旋转柔性曲梁非线性动力学模型.基于该模型,本文仿真计算了悬臂曲梁的纯弯曲问题和带有刚柔耦合效应的旋转柔性曲梁动力学响应问题,以此分别讨论了所提出曲梁单元的收敛性和动力学模型的正确性.进一步应用D'Alembert原理,将旋转曲梁等效为带离心力的无旋转曲梁,通过线性摄动处理得到系统的特征方程,以此分别研究了旋转角速度、初始曲率和集中质量对曲梁动力学特性的影响.最后重点分析了旋转曲梁的频率转向和振型切换问题,并阐述了两者之间的相互关系.研究结果表明:随着旋转角速度的增大,曲梁的频率特性与直梁的频率特性相近,以及曲梁拉伸变形占主导的模态振型会提前.      

输液曲管平面内振动的波动方法研究

采用Flügge曲梁模拟弯曲管道,推导了管内流体的加速度,在总体轴线不可伸长假定的基础上建立了曲管平面内振动的动力学方程;采用波动方法,获得了曲管内振动波的传播和反射矩阵,提出了计算曲管平面内振动固有频率的数值方法。算例分析中,通过计算两端固定半圆形曲管的临界流速并与已有文献结果对比,验证了本文方法的正确性。最后,计算了两端固定半圆形曲管在四种不同流速下的前四阶固有频率,结果表明,管内流速的增大会降低管道的固有频率,当流速增大到某一特定值时,管道的一阶固有频率消失。     

Generalized variational principle of dynamic analysis on naturally curved and twisted box beams for anisotropic materials

An incomplete generalized variational functional for naturally curved and twisted composite box beams with complete constrained boundaries at two ends is established by means of Lagrange multiplier method. The equations of motion governing the dynamic behavior of the beams and corresponding boundary conditions are derived from the stationary condition of the functional. The non-classical influences relevant to the beams are those due to transverse shear deformations, torsion-related warping and several elastic couplings that can arise in composite beams. In order to demonstrate the correctness of the theory developed the natural frequencies and normal mode shapes of the beams under in-plane free vibration are evaluated and compared with the results using PATRAN’s beam elements.     

Nonlinear analysis of members curved in space with warping and Wagner effects

The centroidal axis of a member that is curved in space is generally a space curve. The curvature of the space curve is not necessarily in the direction of either of the principal axes of the cross-section, but can be resolved into components in the directions of both of these principal axes. Hence, a member curved in space is primarily subjected to combined compressive, biaxial bending and torsional actions under vertical (or gravity) loading. In addition, warping actions in particular may occur in curved members with an open thin-walled cross-section, and as the deformations increase, significant interactions of the compressive, biaxial bending and torsional actions occur and profoundly nonlinear deformations are developed in the nonlinear range of structural response. This makes the nonlinear behaviour of a member curved in space very complicated, making it difficult to obtain a consistent differential equation of equilibrium for the nonlinear analysis of members curved in space. In addition, because torsion is one of the primary actions in these members, when the torsional deformations become large, the Wagner effects including both Wagner moment and the conjugate Wagner strain terms are increasingly significant and need to be included in the nonlinear analysis. This paper takes advantage of the merits of so-called “geometrically exact beam theory” and the weak form formulation of the differential equations of equilibrium in beam theory, and it develops consistent differential equations of equilibrium for the nonlinear elastic analysis of members curved in space with warping and Wagner effects. The application of the nonlinear differential equations of equilibrium to various problems is illustrated.     

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.     

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

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

薄壁曲梁的稳定性研究进展

 曲梁是桥梁、建筑、船舶、航空和航天工程中常见的薄壁构件,根据外载荷与主曲率平面的关系,又被称为拱或水平曲梁．随着工程材料的日益发展,如复合材料、功能梯度材料的引入,曲梁的应用范围更加广泛,进一步推进了薄壁曲梁稳定性问题的研究．本文首先对薄壁梁结构的稳定性行为进行了分类．接着简述了薄壁构件的基本假设,对比了近几十年来薄壁曲梁的基本理论,针对复合材料薄壁曲梁,总结了相应的本构关系,并对各理论间存在的分歧进行了归纳．结合最新的薄壁曲梁研究,根据平衡法、能量法和虚位移（虚功）原理推导出控制微分方程,阐述了相应的求解方法,如解析法、半解析法和数值解法．为验证薄壁曲梁理论的准确性,曲梁承载能力试验验证尤为重要,但目前国内外相关研究还很少,亟待发展．最后讨论了现阶段薄壁曲梁研究的局限性和未来发展的方向．     

Elastic bending analysis of bilayered beams by an alternative two-variable method

An alternative two-variable method is used to reanalyze thermoelastic bending problems of bilayered beams subjected to external moments and internal stresses. The differences among zero-stress axis, zero-strain axis (i.e. neutral axis), bending axis, centroidal axis, and the parameter conditions for null/single/dual zero-stress axes are investigated analytically and numerically. Comparisons of thermoelastic stress predictions by the present model with Stoney's model and Hsueh's model are discussed in a representative case of GaAs top coat/Si substrate wafers. Results showed that the neutral axis does not coincide with the zero-stress axis in the general case, and the numbers and the locations of zero-strain or zero-stress axes depend on not only elastic modulus, thickness and/or thermal expansion coefficient ratios between the film and the substrate but also mechanical/thermal loading ratio.     

New shape functions for non-uniform curved Timoshenko beams with arbitrarily varying curvature using basic displacement functions

Basic Displacement Functions (BDFs) are introduced and derived using the basic principles of structural mechanics. BDFs are then utilized to obtain new shape functions for arbitrarily curved non-uniform beams, including the effects of shear deformation and extensibility of neutral axis. The main virtue of the proposed shape functions is their susceptibility to the variations in cross-sectional area, moment of inertia, and curvature along the axis of the beam element. Competence of the present method in static and free vibration analyses, as well as its applicability to the special case of straight Timoshenko beam has been verified through several numerical examples.     

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.     

Large deflection of curved elastic beams made of Ludwick type material

The large deflection of an axially extensible curved beam with a rectangular cross-section is investigated. The elastic beam is assumed to satisfy the Euler-Bernoulli postulation and be made of the Ludwick type material. Through reasonably simplified integration, the strain and curvature of the axis of the beam are presented in implicit formulations. The governing equations involving both geometric and material nonlinearities of the curved beam are derived and solved by the shooting method. When the initial curvature of the beam is zero, the curved beam is degenerated into a straight beam,and the predicted results obtained by the present model are consistent with those in the open literature. Numerical examples are further given for curved cantilever and simply supported beams, and the couplings between elongation and bending are found for the curved beams.     

功能梯度变曲率曲梁的几何非线性模型及其数值解

基于弹性曲梁平面问题的精确几何非线性理论,建立了功能梯度变曲率曲梁在机械和热载荷共同作用下的无量纲控制方程和边界条件,其中基本未知量均被表示为变形前的轴线坐标的函数。以椭圆弧曲梁为例,采用打靶法求解非线性常微分方程的两点边值问题,获得了两端固定功能梯度椭圆弧曲梁在横向非均匀升温下的热弯曲变形数值解,分析了材料梯度指数、温度参数、结构几何参数等对曲梁受力及变形的影响。     

The elastic theory of thin-walled open cross sections with local deformations

A set of governing equations for nonlinear theory of spatially curved elastic beams of thin-walled open cross section composed of straight rectangular elements is presented explicitly in the Lagrangian form. It is shown that local deformations, i.e. in-plane distortion of the cross section may easily be taken into account by the use of the analytical model proposed by Epstein and Murray. The essential feature which distinguishes the present work from Epstein and Murray's is the use of an auxiliary element when the axial curve of beams is not located on the cross section. This enables us to select arbitrarily the axial curve of rods. For the engineering theory of rods, the simplified governing equations for the nonlinear and linear theories with and without local deformations are derived from the rigorous nonlinear theory by employing the thinness assumption. It is also shown that the reduced linear theory without local deformations agrees with the Vlasov theory.     

考虑翘曲效应的圆柱螺旋弹簧的振动分析

以空间曲梁理论为基础对簧丝截面为矩形的圆柱螺旋弹簧的自由振动特性进行了研究. 在弹簧的运动微分方程中, 所有的位移函数和广义翘曲坐标均定义在横截面的形心主轴上, 同时考虑了翘曲变形对弹簧固有频率的影响. 通过精确地应用符号运算软件MATHEMATICA可以得到振动模态的显式表达式, 固有频率则可用搜索的方法来确定. 在较宽的范围内, 给出了各种参数变化, 如簧丝截面的宽高比(a/ b = 0.6 ~ 1.7)、螺旋角(\bar {\alpha } = 5^\circ ~ 12.5^ \circ)、弹簧工作圈数(n = 6 ~ 12)和圆柱螺旋线半径(R = 4~ 10mm)对固有频率的影响. 为了证明解析法的有效性, 对两端固支和一端固支、一端自由矩形截面圆柱螺旋弹簧的固有频率和振动模态进行了求解, 并同ANSYS三维实体单元(Solid45)的结果和文献的结果进行了比较. 计算表明: 用解析方法得到的解和用数值方法得到的结果吻合得很好.      

Comparison of waveguide properties of curved versus straight planar elastic layers

Dispersion equations are solved for the in-plane and anti-plane wave propagation in planar elastic layer with constant curvature. The classical Lamé formulation of displacements via elastic potentials is applied and appropriate simplifications are employed. The dispersion diagrams in each case are compared with their counterparts for a straight layer, e.g., the classical Rayleigh–Lamb solution. The curvature-induced symmetry-breaking effects are investigated for layers with symmetric boundary conditions. The role of curvature is also investigated in the cases, when the boundary conditions are not symmetrical. The elementary Bernoulli–Euler theory is employed to analyze the wave guide properties of a curved planar elastic beam in its in-plane deformation. The validity range of the Bernoulli–Euler theory is assessed via comparison of dispersion diagrams.     

地震栓波器假频本质的动力稳定性解释

地震栓波器是用来栓测竖向地震波的仪器.本文根据它的实际结构,提出了一个较精确的力学模型——双层三悬丝模型.把悬丝作为曲杆,用动力稳定性理论进行分析,揭示了地震仪栓波器由于竖向地震波以外的因素而产生输出的所谓假频现象产生的原因:位于悬丝所在水平面内的激励在一定条件下将导致悬丝在面外动力失稳而产生强烈振动,从而引起虚假输出.计算结果说明了影响假频的诸因素,为激励速度峰值、激励相对于栓波器的方向和曲杆面内振动的固有频率等的作用,并与实验值作了比较.     

Free vibration analysis of third-order shear deformable composite beams using dynamic stiffness method

The dynamic stiffness method is introduced to investigate the free vibration of laminated composite beams based on a third-order shear deformation theory which accounts for parabolic distribution of the transverse shear strain through the thickness of the beam. The exact dynamic stiffness matrix is found directly from the analytical solutions of the basic governing differential equations of motion. The Poisson effect, shear deformation, rotary inertia, in-plane deformation are considered in the analysis. Application of the derived dynamic stiffness matrix to several particular laminated beams is discussed. The influences of Poisson effect, material anisotropy, slenderness and end condition on the natural frequencies of the beams are investigated. The numerical results are compared with the existing solutions in literature whenever possible to demonstrate and validate the present method.     

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

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

Disturbances to a rotating fluid with a sphere moving along the axis of rotation

I.IntroductionWhenabodyprojectedinthewaterisintranslationalmotionthroughwateritwillcertainlycausesdisturbances.Howtodetectthedisturbancesf'arawayfromthebodyandhowtodetermillebytheintbrlnationobtainedthepositionandvelocityofthebodyandthesiteofprojectionareimportant.Astheprojectedbodyisrotating,thefluidaroundisinauniformlyrotationrelativetothereferenceframerotatinginsynchronismwiththebody,andthisisaquestionofthemotionofabodyinarotatingfluid.SuchproblemswerefirststudiedbyProudman(1916)l']andTay…