Monitoring structural damage of components using an effective modulus approach

This paper describes a novel nondestructive damage detection method that was developed to study the influence of a crack on the dynamic properties of a cantilever beam subjected to bending. Experimental measurements of transfer functions for the cracked cantilever beam revealed a change in the natural frequency with increasing crack length. A finite element model of a cracked element was created to compute the influence of severity and location of damage on the structural stiffness. The proposed model is based on the response of the cracked beam element under a static load. The change in beam deflection as a result of the crack is used to calculate the reduction in the global component stiffness. The reduction of the beam stiffness is then used to determine its dynamic response employing a modal analysis computational model. Euler–Bernoulli and Timoshenko beam theories are used to quantify the elastic stiffness matrix of a finite element. The transfer functions from both theories compare well with the experimental results. The experimental and computational natural frequencies decreased with increasing crack length. Furthermore the Euler–Bernoulli and Timoshenko beam theories resulted in approximately the same decrease in the natural frequency with increasing crack length as experimentally measured.     

考虑约束扭转的薄壁梁单元刚度矩阵

推导了薄壁空间梁单元刚度矩阵 ,考虑了双向弯曲及截面约束扭转对杆件轴向变形的影响 ;计算了截面的翘曲变形 ,以及二次剪应力对翘曲变形的影响 ,可适用于任意截面 (包括开口、闭口和混合剖面 )的薄壁杆件。计算结果表明 ,考虑约束扭转的薄壁梁单元刚度矩阵有相当好的精确度 ,可以用于薄壁杆件的静动力分析。     

薄壁杆系结构的梁元分析模型

本文导出了用于薄壁杆系结构弹性分析的薄壁梁元分析模型，在空间梁元分析模型^[3]的基础上，采用了一种改进的位移模式，考察了薄壁杆件可能发生的拉压，剪切，弯曲，扭转和翘曲等各变形形式以及它们的耦合效应，得出了相应的单元形函数，同时从工程应变的定义出发，采用Taylor级数展开的方法，建立了单元的五阶近似正交变表达式，并建立了相应的薄壁单元刚度方程，从而得出了局部坐标系下单元刚度矩阵的显式，根据本文所导出的薄壁梁元分析模型，编制了相应的结构计算程序，通过算例验证了本文所推导的单元刚度矩阵，同时通过与传统空间梁元计算模型计算结果的比较，阐述了截面翘曲对薄壁杆系结构的影响。     

Calculation of higher natural frequencies of simply supported elastic composite beams using Riccati matrix method

A composite beam is composed of an upper slab and a lower beam connected at the interface by shear transmitting studs. In this paper, an improved and efficient numerical model for the calculation of higher natural frequencies of an elastic composite beam is presented. The numerical model uses the Riccati transfer matrix method. First, the exact field transfer matrix for an element of the beam is represented using the combination of the Analog Beam Method and the transfer matrix method. Second, applying Riccati method to the beam system, the natural frequencies can be easily calculated. The advantage of the present model is to overcome the numerical instabilities of the ordinary transfer matrix method, especially when calculating the higher natural frequencies of structures. A numerical example is given to illustrate and compare the results with those available from other methods. Finally, a parametric study is given to examine the effect of various parameters of the elastic composite beam on its free vibration behavior.     

Free vibration and dynamic stability of rotating thin-walled composite beams

The dynamic stability behavior of thin-walled rotating composite beams is studied by means of the finite element method. The analysis is based on Bolotin’s work on parametric instability for an axial periodic load. The influence of fiber orientation and rotating speeds on the natural frequencies and the unstable regions is studied for symmetrically balanced laminates. The regions of instability are obtained and expressed in non-dimensional terms. The “modal interchange” phenomenon arising in rotating beams is described. The dynamic stability problem is formulated by means of linearizing a geometrically nonlinear total Lagrangian finite element with seven degrees of freedom per node. This finite element formulation is based on a thin-walled beam theory that takes into account several non-classical effects such as anisotropy, shear flexibility and warping inhibition.     

Warping effect in free vibration analysis of unidirectional composite non-cylindrical helical springs

The differential equations of motion for unidirectional composite non-cylindrical helical springs including warping, which consist of 14 first-order partial differential equations with variable coefficients, are first derived based on arbitrary spatially curved anisotropic beam theory. An analytical formula for the warping function of Saint-Venant’s torsion of unidirectional composite beams with rectangular cross-section is also obtained. The natural frequencies are determined using improved Riccati transfer matrix method. The element transfer matrix is calculated by the use of the Scaling and Squaring method and Pad’e approximations. Comparisons are made with the EF-results on the natural frequencies of the springs, made from rectangular wire, with inclusion of the warping effect. Information is given on the effect on the natural frequencies of the ratio of radii of the minimum cylinder to the maximum cylinder, the helix pitch angle and the number of active turns. Numerical results reveal that the warping deformation has a significant influence on the natural frequencies, which should be considered in the free vibration analysis of such springs.     

A new beam element for analyzing geometrical and physical nonlinearity

Based on Timoshenko＇s beam theory and Vlasov＇s thin-walled member theory, a new model of spatial thin-walled beam element is developed for analyzing geometrical and physical nonlinearity, which incorporates an interior node and independent interpolations of bending angles and warp and takes diversified factors into consideration, such as traverse shear deformation, torsional shear deformation and their coupling, coupling of flexure and torsion, and the second shear stress. The geometrical nonlinear strain is formulated in updated Lagarange （UL） and the corresponding stiffness matrix is derived. The perfectly plastic model is used to account for physical nonlinearity, and the yield rule of von Mises and incremental relationship of Prandtle-Reuss are adopted. Elastoplastic stiffness matrix is obtained by numerical integration based on the finite segment method, and a finite element program is compiled. Numerical examples manifest that the proposed model is accurate and feasible in the analysis of thin-walled structures.     

轴向受载的Bernoulli—Euler薄壁梁固有振动和稳定性的理论研究

通过直接求解轴向受载的单对称均匀Bernoulli-Euler薄壁梁单元弯扭耦合振动的运动微分方程,推导了其动态传递矩阵,讨论了轴向载荷的变化对薄壁梁弯扭耦合振动固有频率的影响,并由此得到零频率振动（弹性屈出）发生时相应的轴向载荷,数值结果表明本文方法在其适应范围内是精确有效的。     

非对称混杂边界轴向运动Timoshenko梁橫向振动分析

研究两端带有扭转弹簧且弹簧系数均可任意变化的非对称混杂边界下的轴向运动Timoshenko梁的横向振动.利用非对称混杂边界条件推导对应任意弹簧系数的系统超越方程以及特征函数.运用数值方法计算系统的固有频率及其相应的模态函数,并研究确定梁的刚度、轴向速度以及边界处扭转弹簧的刚度的影响.通过数值算例,比较7imoshenko梁、瑞利梁、剪切梁和欧拉梁的固有频率随轴向速度的变化,分析转动惯量和剪切变形的影响.     

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.     

简化有限元方法的波纹管模态分析

采用简化方法对单层U型波纹管进行有限元模态分析。在长度,质量,对轴线转动惯量,体积,轴向弹簧比率以及周向弹簧比率不变的前提下,把波纹管简化成直壁薄管,用直壁薄管模型代替波纹管模型进行有限元模态计算。简化了计算模型,减少了计算量,提高了计算效率。本文给出了明确的等效直壁薄管的几何参数,物性参数的求解公式,改变了Broman人为设定直壁薄管厚度的方法,取得了更加准确的计算结果,使利用直壁薄管模型求解波纹管固有频率具有通用意义。进行了直壁薄管模型得到模态振型与波纹管模型得到模态振型之间的比较,认为直壁薄管模型可以求解波纹管的振型。     

Aerothermoelastic behavior of supersonic rotating thin-walled beams made of functionally graded materials

In this study, a thin-walled beam made of functionally graded material (FGM) which is used as rotating blades in turbomachinery under aerothermoelastic loading is investigated. The governing equations, which are based on first-order shear deformation theory, include the effects of the presetting angle, the secondary warping, temperature gradient through the wall thickness of the beam and also the rotational speed. Moreover, quasi-steady aerodynamic pressure loadings are determined using first-order piston theory, and steady beam surface temperature is obtained from gas dynamics theory. Then, the blade partial differential equations are transformed into a set of ordinary differential equations using the extended Galerkin method. Finally, having solved the resulting structural–fluid–thermal eigenvalue system of equations, the effects of Mach number and geometric parameters on natural frequencies are presented. The results demonstrate that the natural frequencies decrease under aerothermoelastic loading at high Mach numbers.     

拱桥悬臂施工过程中面内特征值的传递矩阵法

采用传递矩阵法对拱桥悬臂施工过程中的面内特征值问题进行求解,建立了该类桥型施工过程中面内竖弯刚度的评估方法.首先,将索和拱分别视为无垂度的张紧弦和欧拉伯努利梁,基于传递矩阵法基本理论推导了系统的总传递矩阵,考虑拱和索的边界条件以及索拱节点的位移连续性条件得到系统的特征值方程,进而计算出系统的频率和模态.同时,采用有限元...     

Free vibration analyses of axially loaded laminated composite beams based on higher-order shear deformation theory

The dynamic stiffness matrix method is introduced to solve exactly the free vibration and buckling problems of axially loaded laminated composite beams with arbitrary lay-ups. The Poisson effect, axial force, extensional deformation, shear deformation and rotary inertia are included in the mathematical formulation. The exact dynamic stiffness matrix is derived from the analytical solutions of the governing differential equations of the composite beams based on third-order shear deformation beam theory. The application of the present method is illustrated by two numerical examples, in which the effects of axial force and boundary condition on the natural frequencies, mode shapes and buckling loads are examined. Comparison of the current results to the existing solutions in the literature demonstrates the accuracy and effectiveness of the present method.     

???????????????????????????о?

针对连续Bernoulli-Euler和Timoshenko梁单元的动态刚度矩阵,分析了在使用连续梁单元 进行结构动态特性分析中的数值问题. 基于连续梁单元的运动方程,导出了连续 Bernoulli-Euler和Timoshenko梁单元的动态刚度矩阵. 分析了影响动态刚度矩阵中双曲函 数自变量的各个独立变量及其产生的影响,并给出了初估连续梁单元合理长度的方法. 使用 单一连续Bernoulli-Euler和Timoshenko梁单元的动态刚度矩阵分别进行了悬臂梁频响曲线 的数值求解. 研究表明,在合理选择连续梁单元的长度时,大多数工程结构的动态特性分析 中都不会产生数值问题.     

框架结构屈曲的精确有限元求解

基于屈曲微分控制方程的一般解,构造了Euler梁在轴力作用下的精确形函数,建立了用于框架结构屈曲分析的精确有限单元,得到了单元刚度矩阵和几何刚度矩阵的显式表达,并提出了基于常规特征值计算的迭代算法以确定屈曲载荷及相应失稳模态的精确解. 研究表明, 对于线性稳定性分析而言,常规框架有限单元可视为精确有限单元的一种近似. 若采用精确单元,无需进行网格细分就可以获得精确的屈曲载荷和失稳模态. 数值算例证明了新单元和算法的效率和精度.      

A new finite element of spatial thin-walled beams

Based on the theories of Timoshenko＇s beams and Vlasov＇s thin-walled members, a new spatial thin-walled beam element with an interior node is developed. By independently interpolating bending angles and warp, factors such as transverse shear deformation, torsional shear deformation and their Coupling, coupling of flexure and torsion, and second shear stress are considered. According to the generalized variational theory of Hellinger-Reissner, the element stiffness matrix is derived. Examples show that the developed model is accurate and can be applied in the finite element analysis of thinwalled structures.     

A quaternion-based formulation of Euler–Bernoulli beam without singularity

This paper proposes a singularity-free beam element with Euler–Bernoulli assumption, i.e., the cross section remains rigid and perpendicular to the tangent of the centerline during deformation. Each node of this two-nodal beam element has eight nodal coordinates, including three global positions and one normal strain to describe the rigid translation and flexible deformation of the centerline, respectively, four Euler parameters or quaternion to represent the attitude of cross section. Adopting quaternion instead of Eulerian angles as nodal variables avoids the traditionally encountered singularity problem. The rigid cross section assumption is automatically satisfied. To guarantee the perpendicularity of cross section to the deformed neutral axes, the position and orientation coordinates are coupled interpolated by a special method developed here. The proposed beam element allows arbitrary spatial rigid motion, and large bending, extension, and torsion deformation. The resulting governing equations include normalization constraint equations for each quaternion of the beam nodes, and can be directly solved by the available differential algebraic equation (DAE) solvers. Finally, several numerical examples are presented to verify the large deformation, natural frequencies and dynamic behavior of the proposed beam element.     

On inner shearing constraints for a direct beam model coarsely describing warping

We investigate inner shearing constraints for a direct one-dimensional beam model coarsely describing warping. In particular, we study how they affect the field equations for the elastic buckling of open thin-walled beams. We show that the distinction between the axes of the shear centres and of the centroids is crucial for the kinematics of the beam, and to derive expression of reactive actions, hence of buckling critical conditions. We also show how a suitable definition of the warping inertial action and properly introduced inner shearing constraints yield the field equations of the literature for the beam in a general dynamic setting.     

剪力滞效应对箱形梁挠度影响的研究

从剪力滞翘曲应力的轴向平衡条件出发,选取双室箱梁的合理翘曲位移函数,引入相应于剪力滞翘曲变形的惯性矩和惯性积等几何特性,用能量变分法建立薄壁箱梁剪力滞效应分析的控制微分方程。通过求解控制微分方程,导出集中荷载和均布荷载作用下简支箱梁和悬臂箱梁的挠度公式及有限梁段单元刚度矩阵,模型试验和ANSYS壳单元计算结果证实了其正确性。结合简支、悬臂和连续箱梁数值算例,具体分析剪力滞效应对箱梁挠度的提高程度。结果表明,无论在集中荷载还是均布荷载作用下,剪力滞效应对简支箱梁的挠度均有显著提高。在集中荷载作用下,剪力滞效应对连续箱梁挠度的提高可达14%;对于跨宽比约为4.0～6.0的简支箱梁,可将按初等梁计算的跨中挠度乘以提高系数1.05～1.11;计算悬臂箱梁的挠度时,一般可以忽略剪力滞效应的影响。