Impulsive loading of ideal fibre-reinforced rigid-plastic beams—III Cantilever beam

The theory outlined in Part I is applied to the problem of a cantilever beam struck transversely at any point by a mass which subsequently adheres to the beam. In the subsequent motion, slope and velocity discontinuities propagate outwards from the point of impact. Solutions for the velocity and deflection of the various segments of the beam are obtained for the case of linear strain-hardening, and simpler approximate solutions are derived for the case of low impact velocity and/or slight strain-hardening. The discontinuity propagating towards the free end of the beam always comes to rest before it reaches this end, but for sufficiently high values of impact mass and velocity, and a strain-hardening parameter, one or more reflections of the discontinuity may occur at the fixed end of the beam and at the point of impact.     

Impulsive loading of ideal fibre-reinforced rigid-plastic beams—II Beam with end supports

The theory outlined in Part I is applied to the problem of a beam, supported at both ends, struck transversely at any point by a mass which subsequently adheres to the beam. For sufficiently long beams, the resulting moving discontinuities in slope and velocity do not propagate to the ends of the beam, and for this case the solution is obtained for the non-linear strain-hardening law used in Part I. A simpler approximate solution is derived for the case of low impact velocity and/or slight strain hardening. For shorter beams, the propagating discontinuities may undergo one or more successive reflections at the end points and at the point of initial impact. The deformation after such reflections is analysed, for linear strain-hardening, in the final two sections. Some particularly simple results are found in the case of central impact.     

Closed form solutions of Euler–Bernoulli beams with singularities

The problem of the integration of the static governing equations of the uniform Euler–Bernoulli beam with discontinuities is studied. In particular, two types of discontinuities have been considered: flexural stiffness and slope discontinuities. Both the above mentioned discontinuities have been modeled as singularities of the flexural stiffness by means of superimposition of suitable distributions (generalized functions) to a uniform one dimensional field. Closed form solutions of governing differential equation, requiring the knowledge of the boundary conditions only, are proposed, and no continuity conditions are enforced at intermediate cross-sections where discontinuities are shown. The continuity conditions are in fact embedded in the flexural stiffness model and are automatically accounted for by the proposed integration procedure. Finally, the proposed closed form solution for the cases of slope discontinuity is compared with the solution of a beam having an internal hinge with rotational spring reproducing the slope discontinuity.     

The growth of plane discontinuities propagating into a homogeneously deformed elastic material

It is shown that, in general, plane acceleration discontinuities propagating into an isotropic elastic material in a state of homogeneous deformation either become infinite in a finite time or decay to zero in an infinite time. Exceptions to this result are transverse discontinuities which propagate along a principal axis of strain without change in strength. Conditions governing the growth of acceleration discontinuities travelling into undeformed material are found to be identical with the thermodynamic conditions derived by Bland [2] for shock propagation. Plane discontinuities of order higher than the second are shown to propagate with constant strength.     

Impulsive loading of ideal fibre-reinforced rigid-plastic beams—I Free beam under central impact

A theory was developed in [1] for the dynamical behaviour under transverse load of ideal fibre-reinforced beams (that is, beams which are inextensible in their longitudinal direction) which exhibit rigid-plastic mechanical response. This theory is here applied to the problem of a beam of finite length, free at both ends, which is struck centrally by a mass which subsequently adheres to the beam. The general solution for the motion of the beam is determined for a fairly wide class of non-linear strain-hardening laws. Simplified approximate solutions are derived for the cases of (a) a heavy striker, (b) a light striker and (c) low impact speed and/or slight strain-hardening.     

On the validity of planar, thick curved beam models derived with respect to centroidal and neutral axes

Most dynamic analyses of planar curved beams found in the literature are carried out based on a curved beam model which assumes that the neutral axis coincides with the centroidal axis of the curved beam. This assumption leads to governing equations of motion which are relatively simple with analysis results that have acceptable accuracy for shallow curved beams. However, when a curved beam is not shallow and/or its cross section is not doubly symmetric, the offset distance between the neutral and centroidal axes may be large enough to influence the in-plane dynamics of the curved beam even for small motion. In this paper, the validity of this underlying assumption for modeling a linear curved beam is examined. To this end, two sets of equations of motion governing the in-plane dynamics of a planar curved beam are derived, in a consistent manner for comparison, based on the linear strain-displacement relations and Hamilton’s principle. The first set of equations is derived from the displacement components measured with reference to the neutral axis of the curved beam while the second set is derived with respect to the centroidal axis of the cross section. The curved beam is considered extensional and the effects of rotary inertia and radial shear deformation are included. In addition to the curvature parameter that characterizes the wave motion for both curved beam models, an eccentricity parameter is introduced in the first model to account for the offset between the neutral and centroidal axes. The dynamic behavior predicted by each curved beam model is compared in terms of the dispersion relations, frequency spectra, cutoff frequencies, natural frequencies and modeshapes, and frequency responses. In order to ensure that the comparison is accurate, the wave propagation technique is applied to obtain exact wave solutions. It is shown that, when the curvature parameter is not small, the underlying assumption has a substantial impact on the accuracy of the linear dynamic analysis of a curved beam.     

Thermal effect on the deformation of a flexible beam with large kinematical driven overall motions

The research investigates the transient longitudinal and transverse deformation of a planar flexible beam with large overall motions in a temperature field. With the increase of temperature, longitudinal deformation is caused by the thermal expansion in axial direction. Due to the coupling of longitudinal and transverse deformation, the transverse deformation is induced, which is significant in cases where temperature increases rapidly in a very short period. Furthermore, the transverse temperature gradient, which is caused by the temperature variation in the transverse direction, may lead to transverse deformation. Considering the thermal strain, equations of motion of a flexible beam with arbitrary large overall motion are derived based on virtual work principle. The high order terms of the strain tensor are taken into account, such that the geometric nonlinear deformation terms are included in the dynamic equations. Simulation results of a rotating beam are shown to reveal the thermal effect and nonlinear effect on the dynamic performance of the beam.     

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.     

Numerical investigation of the slope discontinuities in large deformation finite element formulations

In many multibody system applications, the system components are made of structural elements that can have different orientations, leading to slope discontinuities. In this paper, a numerical investigation of a new procedure that can be used to model structures with slope discontinuities in the finite element absolute nodal coordinate formulation (ANCF) is presented. This procedure can be applied to model slope discontinuities in the case of commutative rotations of gradient deficient elements that are used for modeling thin beam and plate structures. An important special case to which the proposed procedure can be applied is the case of all planar gradient deficient ANCF finite elements. The use of the proposed method leads to a constant orthogonal element transformation that describes an arbitrary initial configuration. As a consequence, one obtains, in the case of large commutative rotations and large deformations, a constant mass matrix for structures which have complex geometry. The procedure used in this investigation to model slope discontinuities requires the use of the concept of the intermediate finite element coordinate system. For each finite element, a new set of gradient coordinates that define, at the discontinuity node, the element deformation with respect to the intermediate element coordinate system is introduced. These new gradient coordinates are assumed to be equal for the two finite elements at the point of intersection. That is, the change of the gradients of two elements at the intersection point from their respective intermediate initial reference configuration is assumed to be the same. This procedure leads to a set of linear algebraic equations that define the orthogonal transformation matrix for the finite element. Numerical examples are presented in order to demonstrate the use of the proposed procedure for modeling slope discontinuities.     

Using the Dirac approach to simulate the rolling of a wheeled vehicle

The rolling of a wheeled vehicle is considered in the case when the turning angles of the front wheels about the vertical axis are small. The small relative slip is taken into account in the adopted model of contact between the wheels and the road surface. It is shown that, when the stiffness of the wheels tends to infinity, the system of equations of motion may become nonclassical and its form is specified by the no-slip conditions along the longitudinal direction of motion and by the primary Dirac constraints arising because of the degeneracy of the Lagrangian.     

旋转中心刚体-FGM梁刚柔热耦合动力学特性研究

对旋转中心刚体-功能梯度材料(functionally graded material,FGM)梁刚柔热耦合动力学特性进行研究.FGM梁为物理性能参数沿厚度方向呈幂律分布的欧拉伯努利梁.考虑柔性梁的横向弯曲变形和轴向拉伸变形, 并计入横向弯曲变形引起的纵向缩短,即非线性耦合变形量.考虑变截面空心梁在外部高温、内冷通道冷却情况下的热力耦合对系统动力学特性的影响,求解得到FGM梁沿厚度方向分布的温度场, 进而在本构关系中计入热应变.采用假设模态法对柔性梁变形场进行离散,运用第二类拉格朗日方程推导得到系统的刚柔热耦合动力学方程,并编制动力学仿真软件, 然后通过仿真算例对系统的动力学问题进行研究.结果表明:不同截面梁动力学响应差异较大, 因此需对实际系统合理建模;大范围运动已知时, 考虑热冲击载荷的FGM梁将有效抑制横向弯曲变形,而大范围运动恒定时随热冲击的叠加会出现高频振荡; 大范围运动未知时,外力矩和热冲击载荷相互作用产生热力耦合效应, 导致系统呈现高频振荡,同时与中心刚体大范围旋转运动产生刚柔热耦合效应.      

考虑剪切效应的旋转FGM楔形梁刚柔耦合动力学建模与仿真

本文对一类中心刚体-柔性梁系统在大范围转动下的刚柔耦合动力学问题进行了研究. 柔性梁为功能梯度材料(functionally graded materials, FGM)楔形变截面梁,材料体积分数在梁轴向呈幂律分布变化. 以弧长坐标来描述柔性FGM梁的几何位移关系,分别使用倾角和拉伸应变变量描述柔性梁的横向弯曲和纵向拉伸变形,并计及剪切效应. 采用假设模态法离散变形场,运用第二类拉格朗日方程进行方程推导,得到系统考虑剪切效应的刚柔耦合动力学模型. 基于全新的刚柔耦合动力学建模理论,研究不同轴向材料梯度分布的FGM楔形梁,通过数值仿真计算,分析讨论不同的转速、梯度分布规律以及变截面参数对系统动力学特性的影响. 结果表明,剪切效应对大高跨比的FGM楔形梁的变形影响较为明显,不容忽略;材料梯度分布规律和截面参数的选取均会对旋转FGM楔形梁的动力学响应和频率产生较大影响. 本文提出的考虑剪切效应的倾角刚柔耦合动力学模型是对以往非剪切模型的进一步完善,可应用于工程中的 Timoshenko梁结构的动力学问题求解.      

Towards linear modal analysis for an L-shaped beam: Equations of motion

We consider an L-shaped beam structure and derive all the equations of motion considering also the rotary inertia terms. We show that the equations are decoupled in two motions, namely the in-plane bending and out-of-plane bending with torsion. In neglecting the rotary inertia terms the torsional equation for the secondary beam is fully decoupled from the other equations for out-of-plane motion. A numerical modal analysis was undertaken for two models of the L-shaped beam, considering two different orientations of the secondary beam, and it was shown that the mode shapes can be grouped into these two motions: in-plane bending and out-of-plane motion. We compared the theoretical natural frequencies of the secondary beam in torsion with finite element results which showed some disagreement, and also it was shown that the torsional mode shapes of the secondary beam are coupled with the other out-of-plane motions. These findings confirm that it is necessary to take rotary inertia terms into account for out-of-plane bending. This work is essential in order to perform accurate linear modal analysis on the L-shaped beam structure.     

CRACK PROPAGATION IN STRUCTURES SUBJECTED TO PERIODIC EXCITATION

In the present paper,a simple mechanical model is developed to predict the dynamic response of a cracked structure subjected to periodic excitation,which has been used to identify the physical mechanisms in leading the growth or arrest of cracking.The structure under consideration consists of a beam with a crack along the axis,and thus,the crack may open in Mode I and in the axial direction propagate when the beam vibrates.In this paper,the system is modeled as a cantilever beam lying on a partial elastic foundation,where the portion of the beam on the foundation represents the intact portion of the beam.Modal analysis is employed to obtain a closed form solution for the structural response.Crack propagation is studied by allowing the elastic foundation to shorten(mimicking crack growth)if a displacement criterion,based on the material toughness,is met.As the crack propagates,the structural model is updated using the new foundation length and the response continues.From this work,two mechanisms for crack arrest are identified.It is also shown that the crack propagation is strongly influenced by the transient response of the structure.     

正交各向异性板中非主轴方向的Lamb波

基于线性三维弹性理论,采用勒让德正交多项式展开法,推导了波沿正交各向异性材料非主轴方向传播时的Lamb波耦合波动方程,并对耦合波动方程进行了数值求解。为验证该方法的适用性和正确性,首先将此方法应用于各向同性材料,并与已知的数据结果进行了比较;然后以单向纤维增强复合材料为例,计算了耦合Lamb波沿不同的非主轴方向传播时的相速度频散曲线,并分别研究了传播方向改变时低阶模态Lamb波和高阶模态Lamb波频散特性的变化。最后,针对潜在用于各向异性复合材料结构健康监测的耦合Lamb波低阶模态,给出了其在不同传播方向时的相速度分布和群速度分布。同时,结合低阶模态Lamb波的位移分布特性和材料的各向异性特点,阐释了S0模态对波的传播方向变化最为敏感的原因。     

A non-linear model for the dynamics of open cross-section thin-walled beams—Part II: forced motion

The discrete equations developed in Part I are here used to analyze the non-linear dynamics of an inextensional shear indeformable beam with given end constraints. The model takes into account the non-linear effects of warping and of torsional elongation. Non-linear 3D oscillations of a beam with a cross-section having one symmetry axis is examined. Only terms of higher magnitude are retained in the equations, which exhibit quadratic, cubic and combination resonances. A harmonic load acting in the direction of the symmetry axis and in resonance with the corresponding natural frequency, is considered. Steady-state solutions and their stability are studied; in particular the effects of non-linear warping and of torsional elongation on the response are highlighted.     

A finite element approximation of steady flow through a rotating non-aligned straight tube

A finite element solution is developed for a penalty function formulation of the equations which govern the steady motion of a Newtonian fluid through a pipe that rotates about an axis not parallel to its own. The motion in this system is driven by the Coriolis acceleration, which has components in the axial direction as well as in the transverse plane of the pipe. The relative magnitudes of these components significantly affect the qualitative and quantitative nature of the primary and secondary flow field. The present results compare favourably with those of previously reported experimental and theoretical studies over a wide range of flow regimes.     

Dynamic steady crack growth in elastic-plastic solids—propagation of strong discontinuities

The near-tip field of a mode I crack growing steadily under plane strain conditions is studied. A key issue is whether strong discontinuities can propagate under dynamic conditions. Theories which impose rather restrictive assumptions on the structure of an admissible deformation path through a dynamically propagating discontinuity have been proposed recently. Asymptotic solutions for dynamic crack growth, based on such theories, do not contain any discontinuities. In the present work a broader family of deformation paths is considered and we show that a discontinuity can propagate dynamically without violating any of the mechanical constitutive relations of the material. The proposed theory for the propagation of strong discontinuities is corroborated by very detailed finite element calculations. The latter shows a plane of strong discontinuity emanating from the crack tip (with its normal pointing in the direction of crack advance) and moving with the tip. Elastic unloading ahead of and/or behind the plane of discontinuity and behind the crack tip have also been observed.The numerical investigation is performed within the framework of a boundary layer formulation whereby the remote loading is fully specified by the first two terms in the asymptotic solution of the elasto-dynamic crack tip field, characterized by K₁, and T. It is shown that the family of near-tip fields, associated with a given crack speed, can be arranged into a one-parameter field based on a characteristic length, L_g, which scales with the smallest dimension of the plastic zone. This extends a previous result for quasi-static crack growth.     

Stabilization of beam parametric vibrations with shear deformations and rotary inertia effects

The purpose of this theoretical work is to present a stabilization problem of beam with shear deformations and rotary inertia effects. A velocity feedback and particular polarization profiles of piezoelectric sensors and actuators are introduced. The structure is described by partial differential equations with time-dependent coefficient including transverse and rotary inertia terms, general deformation state with interlaminar shear strains. The first order deformation theory is utilized to investigate beam vibrations. The beam motion is described by the transverse displacement and the slope. The almost sure stochastic stability criteria of the beam equilibrium are derived using the Liapunov direct method. If the axial force is described by the stationary and continuous with probability one process the classic differentiation rule can be applied to calculate the time-derivative of functional. The particular problem of beam stabilization due to the Gaussian and harmonic forces is analyzed in details. The influence of the shear deformations, rotary inertia effects and the gain factors on dynamic stability regions is shown.     

Dynamics of a semi-infinite beam on unilateral springs: Touch-down points motion and detached bubbles propagation

The dynamics of a semi-infinite Bernoulli-Euler beam laid on a bed of unilateral elastic springs is governed by a moving-boundary problem, since the positions of the touch-down points, those points which separate the detached beam parts from the laid ones, are unknown. This problem is solved numerically by means of a self-made finite element code and some numerical results are shown and discussed. The nonlinear and non-smooth effects of the touch-down points motion on the beams dynamics are analyzed. The presence of detached bubbles, which appear, propagate and disappear in the beam, is investigated, and new complex motions are highlighted.