基于广义函数空间的不连续梁振动分析

首先运用广义函数建立了轴向力作用下含任意不连续点的弹性基础Euler(欧拉)梁的自由振动的统一微分方程.不连续点的影响由广义函数(Dirac delta函数)引入梁的振动方程.微分方程运用Laplace变换方法求解;与传统方法不同的是,该文方法求得的模态函数为整个不连续梁的一般解.由于模态函数的统一化以及连续条件的退化,特征值的求解得到了极大地简化.最后,以梁-质量块模型和轴向力作用下弹性基础裂纹梁模型为例验证了该文方法的正确性与有效性.     

轴向运动梁的横向振动分析

论述了轴向运动梁横向振动问题以及研究轴向运动梁横向振动问题的方法,指出对轴向运动梁横向振动问题研究中存在的一些错误并进行了更正.针对一端可看作固定边界条件的轴向运动悬臂梁,基于连续体的模态叠加法,推导出含自重效应的轴向运动梁动力响应的计算公式,进行实例计算,并对计算结果进行了详细的讨论,得出影响轴向运动梁振动响应的因素主要有速度和运动方向.     

基于广义弹性理论的微梁固有频率及模态的尺寸效应

经典弹性理论在近代工程技术中得到广泛应用,但其本构关系中不包含任何与尺寸相关的参数,因此不适用于微观结构,不能预测和解释尺寸效应.广义弹性理论增加了偶应力及其对应的曲率张量,完善了对小变形的几何描述,适用于微结构的尺寸效应研究.该文采用广义弹性理论,并结合Hamilton变分原理推导了悬臂微梁的振动微分方程,对微梁的固有频率及其模态进行了分析.结果表明,随着微梁厚度的不断减小,固有频率的尺寸效应与其对应的模态密切相关.扭转和弯曲模态包含了旋转变形,其对应的固有频率显著提高,表现出了显著的尺寸效应;而拉压模态不涉及旋转变形,固有频率未产生明显变化,没有尺寸效应.     

新型空间薄壁梁单元

基于Timoshenko梁理论和Vlasov薄壁杆件约束扭转理论,建立了具有内部结点的新型空间薄壁截面梁单元．通过对弯曲转角和翘曲角采取独立插值的方法,考虑了横向剪切变形,扭转剪切变形及其耦合作用,弯曲变形和扭转变形的耦合以及二次剪应力等因素影响,由Hellinger-Reissner广义变分原理,推得单元刚度矩阵．算例表明所建模型具有良好的精度,可用于空间薄壁杆系结构的有限元分析．     

基于模态自适应的大变形多体系统动力学分析

柔性大变形系统在进行模态降阶时，若模态选取不当，会影响求解精度甚至导致求解结果发散.对此，提出了基于绝对节点坐标法(ANCF)的柔性大变形系统模态自适应选择方法.通过ANCF梁单元建立系统的动力学模型；利用全模态稀疏表示内部区域的坐标；根据Latin超立方抽样构建采样矩阵，作用于动力学方程，以减少方程的数量；以采样后的动力学方程作为约束，构造模态坐标范数优化问题；求解优化问题可以得到具有重大贡献的模态.通过两个实例表明：数值计算结果与常用方法的结果高度吻合并且求解效率显著提升.     

一维纳米准晶层合梁的非局部振动、屈曲与弯曲研究

基于非局部理论,建立了一维纳米准晶层合简支深梁模型,研究了其自由振动、屈曲行为及其弯曲变形问题.采用伪Stroh型公式,导出了纳米梁的控制方程,并通过传递矩阵法获得简支边界条件下纳米准晶层合梁固有频率、临界屈曲载荷及弯曲变形广义位移和广义应力的精确解.通过数值算例,分析了高跨比、层厚比、叠层顺序及非局部效应对一维纳米准晶层合简支梁固有频率、临界屈曲载荷和弯曲变形的影响.结果表明：固有频率和临界屈曲载荷随着非局部参数增大而减小;外层准晶弹性常数更高时,固有频率和临界屈曲载荷更大;叠层顺序对纳米准晶梁的力学行为有较大影响.所得的精确解可为纳米尺度下梁结构的各种数值方法和实验结果提供参考.     

末端质量作用下变长度轴向运动梁的振动特性

对带集中质量,变长度(或速度)轴向运动梁的振动特性采用两种精确方法求解.首先,对变长度轴向运动Euler(欧拉)梁横向自由振动方程进行化简,通过复模态分析得到本征方程,并在有集中质量的边界条件下得到频率方程,用数值方法求解固有频率和模态函数.然后,采用有限元方法建立运动梁自由振动的方程,求解矩阵方程得到复特征值和复特征向量,结合形函数得到复模态位移.最后,将两种方法的计算结果进行了分析和对比.数值算例的结果表明:不同的轴向运动速度和集中质量对变长度轴向运动梁的振动特性有显著影响,两种计算方法的结果接近且均有效.     

匀速转动内悬臂梁系统的建模与分岔研究

应用广义Hamilton原理建立旋转刚环内接悬臂梁系统在材料本构关系为线性时Rayleigh梁大挠度非线性动力学模型,并讨论当假定梁为Euler-Bernoulli型时的蜕化结果.研究表明,上述两种模型对系统的临界分岔值及分岔解静态没有影响.据此,采用了假设模态法解析地研究了匀速转动的Euler-Bernoulli梁模型的分岔行为,得到两种在物理上存在的分岔类型.为校验解析分析的结果,使用了有限元及打靶法两种数值研究工具.计算结果肯定了研究结论.     

两自由度非对称三次系统非线性模态的奇异性质

利用非线性模态子空间的不变性和摄动技术,研究两自由度非对称三次系统在奇异条件下系统的性质.重点考虑子系统之间线性耦合退化时的奇异性质.对于非共振情形,所得到的解析结果表明,系统出现单模态运动以及振动局部化现象,这种现象的强弱不但与非线性耦合刚度有关,而且与非对称参数有关.并解析地得到了参数的门槛值;对于1:1共振情形,模态随非线性耦合刚度和非对称参数的变化会出现分岔,得到了参数分岔集以及模态的分岔曲线.     

纵向数据广义部分线性模型的惩罚GMM估计

对纵向数据的部分线性模型,通常的做法是用样条方法或者核方法逼近非参数部分,然后再用广义估计方程的估计方法去估计参数部分.本文使用P-样条拟合非参数函数,对不同的矩条件用不同的广义矩方法对模型的参数和非参数进行估计,并且给出了估计量的大样本性质;并用计算机模拟和实例证明了当模型中存在不同的矩条件时,采用不同的惩罚广义矩方法可以显著地提高估计精度.     

Dynamic modeling and analysis of rotating FG beams for capturing steady bending deformation

A dynamic analysis of rotating functionally gradient (FG) beams is presented for capturing the steady bending deformation by using a novel floating frame reference (FFR) formulation. Usually, the cross section of bending beams should rotate round the point at the neutral axis while centrifugal inertial forces are supposed to act on centroid axis. Due to material inhomogeneity of FG beams, centroid and neutral axes may be in different positions, which leads to the eccentricity of centrifugal forces. Thus, centrifugal forces can be divided into three componets: transverse component, axial component and force moment acting on the points of the neutral axis, in which transverse component and force moment can make the beam produce the steady bending deformation. However, this speculation has not been presented and discussed in existing literatures. To this end, a novel FFR formulation of rotating FG beams is especially developed considering centroid and neutral axes. The FFR and its nodal coordinates are used to determine the displacement field, in which kinetic and elastic energies can be accurately formulated according to centroid and neutral axes, respectively. By using the Lagrange's equations of the second kind, the nonlinear dynamic equations are derived for transient dynamics problems of rotating FG beams. Simplifying the nonlinear dynamic equations obtains the equilibrium equations about inertial and elastic forces. The equilibrium equations can be solved to capture the steady bending deformation. Based on the steady bending state, the nonlinear dynamic equations are linearized to obtain eigen-frequency equations. Transient responses obtained from the nonlinear dynamic equations and frequencies obtained from the eigen-frequency equations are compared with available results in existing literatures. Finally, effects of material gradient index and angular speed on the steady bending deformation and vibration characteristics are investigated in detail.     

微曲输流管道振动固有频率分析与仿真北大核心CSCD

首次建立了基于Timoshenko梁理论的微曲输流管道横向振动的动力学模型,并分析了流体流动影响下微曲管道横向自由振动的固有特征.采用广义Hamilton原理,导出了考虑流体影响的微曲管道横向振动的控制方程,通过Galerkin截断对控制方程离散化,再由广义本征值问题得到管道横向振动的固有频率,并研究了液体流速和弯曲幅度对管道横向固有振动特征的影响.发展了基于等效刚度和等效阻尼方法的考虑流体影响的微曲管道振动分析的有限元仿真计算方法,并通过有限元软件实现数值仿真,验证了Galerkin截断的分析结果以及所建立的Timoshenko微曲管道动力学模型的有效性.研究表明,流体的流速以及管道的弯曲幅度对管道横向振动固有频率均有显著影响.     

An exact transfer matrix expression for bending vibration analysis of a rotating tapered beam

This paper presents a transfer matrix expression that can be used to determine the eigenpairs of a rotating beam with a cross section height that linearly decreases along the length of the beam element. The proposed method considers the effect of centrifugal force, including the effects of the axial force, hub radius, and taper ratio. Differential equations are solved for the in-plane bending vibration using the Frobenius method for a power series. The effect of the rotational speed on the eigenpairs of a rotating tapered beam is first investigated, followed by an examination of the contribution rates of the bending strain and additional strain energies generated by centrifugal forces for each mode by analyzing the variation of the energies computed from the strain and kinetic energies. To compute these contribution rates, we used a shape function that was defined by the displacement at both ends of the beam elements. The effect of tapering on the eigenfrequencies of the transverse vibration of rotating beams is analyzed by using various examples, and the contribution rates are examined by using taper ratios of 0 and 0.5.     

Mixed finite element models for free vibrations of thin-walled beams

Simple mixed finite element models are developed for the free vibration analysis of curved thin-walled beams with arbitrary open cross section. The analytical formulation is based on a Vlasov's type thin-walled beam theory which includes the effects of flexural-torsional coupling, and the additional effects of transverse shear deformation and rotary inertia. The fundamental unknowns consist of seven internal forces and seven generalized displacements of the beam. The element characteristic arrays are obtained by using a perturbed Lagrangian-mixed variational principle. Only C⁰ continuity is required for the generalized displacements. The internal forces and the Lagrange multiplier are allowed to be discontinuous at interelement boundaries.

Numerical results are presented to demonstrate the high accuracy and effectiveness of the elements developed. The standard of comparison is taken to be the solutions obtained by using two-dimensional plate/shell models for the beams.     

Dynamic response of clamped axially moving beams: Integral transform solution

The generalized integral transform technique (GITT) is employed to obtain a hybrid analytical-numerical solution for dynamic response of clamped axially moving beams. The use of the GITT approach in the analysis of the transverse vibration equation leads to a coupled system of second order differential equations in the dimensionless temporal variable. The resulting transformed ODE system is then solved numerically with automatic global accuracy control by using the subroutine DIVPAG from IMSL Library. Excellent convergence behavior is shown by comparing the vibration displacement of different points along the beam length. Numerical results are presented for different values of axial translation velocity and flexural stiffness. A set of reference results for the transverse vibration displacement of axially moving beam is provided for future co-validation purposes.     

The Vibrational Behavior of Coupled Bladed Disks with Variable Rotational Speed

Low pressure steam turbine blades are subjected to high static and dynamic loads during operation. These loads strongly depend on the turbine's rotational speed, leading to entirely new load conditions. To avoid high dynamic stresses due to the forced vibrations, a coupling of the blades, such as shrouds or snubber coupling, is applied to reinforce the structure. In this work the influence of the rotational speed on the vibration behavior of shrouded blades is investigated. Two fundamental phenomena are considered: the stress stiffening and the spin softening effect. Both effects are caused by centrifugal forces and affect the structural mechanical properties, i.e. the stiffness matrix K , of the rotating system. Since the rotational speed Ω appears quadratically, it is possible to derive the stiffness matrix as a second order matrix polynomial in Ω² [3]. In the case of shrouded blades, contact forces between neighboring blades must be taken into account. The contact status and the pressure distribution in particular is strongly influenced by the rotational speed, respectively, centrifugal forces, caused by the untwisting and radial deformation of the blades. For the calculation, a three dimensional structural mechanical model including a spatial contact model is considered. The solution of the nonlinear equations of motion is based on the well known Multiharmonic Balance Method [2]. Here, the nonlinear forces are computed in the time domain and transferred in the frequency domain by the use of the Fast Fourier Transformation (FFT), also known as the Alternating Frequency Time method (AFT) [1]. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Spatial Geometric Nonlinearity Spline Beam Element With Nodal Parameters Containing Strains北大核心CSCD

Many slender rods in engineering can be modeled as Euler-Bernoulli beams. For the analysis of their dynamic behaviors, it is necessary to establish the dynamic models for the flexible multi-body systems. Geometric nonlinear elements with absolute nodal coordinates help solve a large number of dynamic problems of flexible beams, but they still face such problems as shear locking, nodal stress discontinuity and low computation efficiency. Based on the theory of large deformation beams’ virtual power equations, the functional formulas between displacements and rotation angles at the nodes were established, which can satisfy the deformation coupling relationships. The generalized strains to describe geometric nonlinear effects in this case were derived. Some parameters of boundary nodes were replaced by axial strains and sectional curvatures to obtain a more accurate and concise constraint method for applying external forces. To improve the numerical efficiency and stability of the system’s motion equations, a model-smoothing method was used to filter high frequencies out of the model. The numerical examples verify the rationality and effectiveness of the proposed element. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Multi-scale modeling of beam-like structures: A new boundary condition concept for the RVE

In this work a coupled two-scale beam model using Timoshenko beam elements [1] with finite displacements on the macro scale and fully non-linear 3D brick elements on the micro scale is proposed. The calculation is carried out with the so-called FE² concept. To achieve the coupling between the beam and the brick elements, the algorithm from [2] is adapted. Within the degenerated concept of the Timoshenko beam, the introduction of a pure shear deformation leads to significant problems concerning the equilibrium condition on the micro scale. Applying this deformation mode on the RVE with periodic boundary conditions results in a rigid body rotation. Using linear displacement boundary conditions instead, the wrapping deformation is suppressed on the boundary, leading to a length dependency in the torsional deformation mode. In addition, the shear forces introduce a bending moment, which depends on the length of the RVE and adds spurious normal stresses and a length dependency of the shear stiffness. To overcome these problems, periodic boundary conditions are applied and the displacement assumptions are modified such that the shear deformation is achieved with force pairs on both ends of the RVE. The resulting model leads to length independent results in tension, bending and torsion and a domain which is able to produce a pure shear stress state. Consequently, only this domain of the model should be homogenized which can be accomplished by modifying the variations in the algorithm [2]. The concept is validated by simple linear and non-linear test problems. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

一类指定应力问题的变分原理与应用

针对有限元分析中对应力或内力有指定条件的问题,引入非弹性应变作为实现指定应力条件的附加未知量,在小变形条件下描述了指定应力条件应当满足的弹性力学控制方程;以位移和未知非弹性应变作为独立变量,建立了具有指定应力条件问题的势能变分原理和虚功方程;以位移、弹性应变、未知非弹性应变和应力为独立变量,建立了一个含四类变量的广义变...