DERIVATION OF NONLINEAR WAVE EQUATION FOR FLEXURAL MOTIONS OF AN ELASTIC BEAM TRAVELLING IN AN AIR-FILLED TUBE

Asymptotic analysis is carried out to derive a nonlinear wave equation for flexural motions of an elastic beam of circular cross-section travelling along the centre-axis of an air-filled, circular tube placed coaxially. Both the beam and tube are assumed to be long enough for end-effects to be ignored and the aerodynamic loading on the lateral surface of the beam is considered. Assuming a compressible inviscid fluid, the velocity potential of the air is sought systematically in the form of power series in terms of the ratios of the tube radius to a wavelength and of a typical deflection to the radius. Evaluating the pressure force acting on the lateral surface of the beam, the aerodynamic loading including the effects of finite deflection as well as of air's compressibility and axial curvature of the beam are obtained. Although the nonlinearity arises from the kinematical condition on the beam surface, it may be attributed to the presence of the tube wall. With the aerodynamic loading thus obtained, a nonlinear wave equation is derived, whereas linear theory is assumed for the flexural motions of the beam. Some discussions are given on the results.     

Analytical formulation of stresses in curved composite beams

Summary This paper outlines an analytical method for computing normal and shear stresses generated in a curved laminated beam under bending loads. Each cross section is assumed to be symmetrical and loads are applied in the plane of symmetry. We build a statically admissible stress field in order to plot normal and shear stress distributions. Received 5 March 1997; accepted for publication 18 September 1997     

Calculation of the motion of a gas behind the front of a luminous detonation wave taking account of the lateral expansion of the plasma column

Nonsteady-state gasdynamic processes are considered in a plasma column (hot channel), formed behind the front of a shock wave moving toward a laser beam. A quasi-one-dimensional approximation is used-the parameters in the channel are assumed to be compensated with respect to cross section, but depend on the time and distance along the axis. Motion in the cold dense casing surrounding the channel is assumed to be one-dimensional and cylindrically symmetrical. The solutions of the corresponding systems of equations in partial derivatives permit the parameters to be determined approximately both in the case when the mean free path of the radiation is small in comparison with the radius of the beam (luminous detonation) and also in the case when the mean free path is comparable with the radius. Examples are given of the corresponding numerical calculations. It is shown that in the latter case, a cycle of incomplete absorption can be achieved when behind the shock wave front, moving with constant velocity up to the Jouguet plane, only a part of the radiation energy incident on the front is released.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 18–28, May–June, 1976.     

夹层梁低速冲击下的接触定律与结构响应

在考虑面板面内拉伸刚度等非线性因素的基础上,将夹层梁上面板视为弹塑性地基上的梁,由其平衡方程推导准静态压入接触定律. 进而研究了底面固定于刚性平面的夹层梁的低速冲击响应,验证了所推导的准静态压入接触定律在低速冲击下依然足够精确. 最后,用能量守恒模型,离散模型和连续模型分析了两端简支夹层梁的低速冲击响应以及接触持续时间.同时讨论了冲头与夹层梁质量比和初速度对结构响应以及分析模型适用性的影响. 理论分析结果与数值模拟结果比较表明:推导的接触定律是准确有效的. 给定初速度,质量比很小时,只有连续模型有效,而当质量比较大时3种模型均有效;给定质量,初速度较大时只有连续模型比较有效. 连续模型预测的接触持续时间与初速度无关.      

局部裂纹损伤简支梁的曲率模态特性

将裂缝损伤简化成矩形凹槽,采用delta函数表示简支梁的裂纹损伤位置,得到了全梁范围内截面转动惯量和单位长度质量的表达式,建立了局部裂缝损伤简支梁的横向自由振动方程.利用摄动方法给出了裂纹摄动项的一般表达式,根据摄动项和完整梁都同时满足边界条件的特点,将一阶和二阶摄动项都表示成完整梁模态的线性组合,结合delta函数的性质,最终获得了受损简支梁的特征值和模态振型的解析表达式.最后,通过数值计算得到结构模态参数,对比了一阶摄动和二阶摄动对计算结果的影响,分析了不同阶固有频率和模态曲率的变动量,为简支梁的损伤监控和检测提供了理论依据.     

A STUDY ON THE COUNTER-INTUITIVE BEHAVIORS OF PIN-ENDED BEAMS UNDER PROJECTILE IMPACT

The counter-intuitive behaviors of pin-ended beams under the projectile impact are investigated with ANSYS/LS-DYNA in this paper. It studies in detail their displacement-time history curves, final deformed shapes, energy relationships and projectile impact velocity ranges related to their counter-intuitive behaviors. The influences of the impact positions on their counterintuitive behaviors are also discussed. The results show that no matter where the impact position on the beam is, the counter-intuitive behaviors of pinned beams will occur as long as the impacting velocity lies within a proper range. Corresponding to the occurring of the counter-intuitive behaviors, the rebounding number in the displacement history curves of the beams decreases from a few times to zero with an increase of the impact velocity. The final deformation modes of the beam corresponding to the counter-intuitive behaviors will appear in symmetrical and unsym-metrical forms no matter where the impact position is; the impact velocity of the first-occurring of the counter-intuitive behaviors of the beam increases slowly with the deviation of the impact position away from the mid-span.     

Focusing of an ultrasonic beam by a curved interface

We calculate the wavefields near the caustics, and their cusps, formed when a well collimated, ultrasonic beam scatters from a concave fluid-solid interface. The radius of curvature of the interface is assumed to be sufficiently large and the angle of incidence sufficiently small that only reflection and transmission of the beam need be considered. The incident beam is modeled as a two-dimensional wavefield whose initial profile is rectangular. The aperture is assumed to be sufficiently large (in wavelengths) that a well collimated beam is radiated, and the interface is assumed to lie close enough to the aperture that it is struck by a wavefield that has not yet evolved into a cylindrical wave. The scattered wavefields are represented as multiple integrals and are evaluated using a combination of asymptotic and numerical analysis. Special attention is given to the sometimes competing effects of the shadow boundary of each scattered beam with its corresponding caustic, and cusp if one is formed.     

径向基点插值法在旋转柔性梁动力学中的应用

将无网格径向基点插值法用于旋转柔性梁的动力学分析. 利用无网格方法对柔性梁的变形场进行离散,考虑梁的纵向拉伸变形和横向弯曲变形,并计入横向弯曲变形引起的纵向缩短,即非线性耦合项,运用第二类拉格朗日方程推导得到系统刚柔耦合动力学方程. 将无网格径向基点插值法的仿真结果有限元法和假设模态法进行比较分析,说明假设模态法的局限性,并表明其作为一种柔性体离散方法在刚柔耦合多体系统动力学的研究中具有可推广性,并讨论了径向基形状参数的影响. 同时运用3 种求解系统动力学方程的方法:纽马克方法、4阶龙格库塔法、亚当姆斯预报校正法,并比较各方法的计算效率, 结果表明纽马克方法最快.      

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

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

基于车桥耦合理论的梁式桥阻尼比识别研究

提出一种基于车桥耦合动力学理论的梁式桥阻尼比识别方法. 首先按照动力学理论将测试车设计为单自由度体系, 然后利用安装在测试车上的传感器采集信号, 从测试车与桥梁接触点响应信号中得到梁式桥响应的信号, 基于车桥耦合动力学原理滤波处理得到包含梁式桥第一阶频率的信号, 最后假定梁式桥阻尼比值, 通过假定的梁式桥阻尼比值获取假定的梁式桥第一阶振型, 不断循环直至假定的阻尼比值下计算的第一阶振型最大值点居中, 即为识别的梁式桥真实阻尼比. 本文首先从车桥耦合动力学理论推导上说明了该方法的可行性, 然后考虑在不同车速与非恒定车速、路面粗糙度、环境噪音等影响因素下进行车桥耦合动力学模型分析, 最后通过实桥试验进行了初步验证. 研究结果表明: 该方法能一定程度上克服外界不利因素的影响, 达到识别梁式桥阻尼比的目的, 为识别梁式桥阻尼比提供一种更优方法, 其具有参数设置较少、操作简单方便以及更高测试精度等优点, 同时有助于推动基于车桥耦合的车桥耦合动力学理论技术在梁式桥模态参数识别工作中的实际工程应用.      

梁-桩-土竖向耦合振动特性分析

针对桩与连续梁组合这一基础形式,研究了梁—桩—土竖向耦合振动特性. 首先,将桩假设为一维黏弹性杆件,采用平面应变模型模拟桩侧成层土对桩的动力作用. 同时,桩顶对梁的支撑简化为竖向点支撑. 然后,分别根据一维杆件纵向振动理论和Timoshenko梁理论给出桩和梁的竖向振动控制方程并求解,进而借助离散傅里叶逆变换获得在梁上瞬时激振下梁和土层以上桩段的时域响应半解析解. 通过与有限元模拟结果对比,验证了解的合理性. 在此基础上,讨论了梁几何参数和桩身缺陷对梁—桩—土竖向耦合振动特性的影响.      

中心刚体-变截面梁系统的动力学特性研究

对在平面内做大范围转动的中心刚体-变截面梁系统的动力学进行了研究.考虑柔性梁横向弯曲变形和纵向伸长变形, 且在纵向位移中计及由于横向变形而引起的纵向缩短项, 即非线性耦合变形项. 采用假设模态法描述变形, 运用第二类Lagrange方程推导得到系统刚柔耦合动力学方程. 在此基础上对做大范围旋转运动的中心刚体-楔形梁以及中心刚体-梯形梁模型的动力学进行了详细研究. 研究表明: 梁宽比、梁高比以及梯形梁变截面位置都对系统的动力学特性有很大影响.      

Chaotic and asymmetrical beam response to impulsive load

Both symmetrical and asymmetrical final displacements are observed for elastic–plastic beams under symmetrical impulsive loading. A three-degree-of-freedom Shanley-type model is developed in this study, which is capable of revealing chaotic and asymmetrical responses of an elastic–plastic beam by introducing initial imperfections. To identify the asymmetrical displacement, the beam response is decomposed into three vibration modes. Corresponding modal participation factors are derived based on the displacement of the three-degree-of-freedom beam model. Phase plane trajectories, Poincaré maps and power spectral density diagrams are derived to illustrate both the symmetrical and asymmetrical chaotic vibrations. Numerical simulations using a general-purpose FE code LS-DYNA are carried out for an elastic–plastic beam subjected to impulsive load. The simulation results indicate that the elastic–plastic beam demonstrates chaotic and asymmetrical vibration when the applied impulsive load exceeds a critical value, which agrees with experimental observations.     

圆柱壳在两种非轴对称同步移动载荷作用下的动力响应

本文根据工程实例计算的需要,研究了有限长弹性圆柱薄壳在两种非轴对称同步移动载荷作用下的动力响应问题。两种非轴对称同步移动载荷作用是指非轴对称移动的集中载荷,以及同步移动且作用范围随移动位置增加的均布载荷的共同作用。建立了在上述两种不同类型载荷作用下的具有对称形式的动力学微分方程组;分别采用Dirac函数与Heaviside函数表示移动的均布载荷与集中载荷,设定位移函数的基础上,应用Galerkin法及Laplace变换,求得了圆柱薄壳应力与位移动态响应的解析解;通过具体算例,将所得到解析解的计算结果与ANSYS数值解进行了对比分析,验证了解析解的可靠性。     

Chaos in a single equilibrium point system: Finite deformations

The purpose of this paper is to examine a highly nonlinear model of a slender beam which yields chaotic solutions for some forcing amplitudes. The study is unique in that the governing partial differential equations are solved directly, and that the model lends itself to a more physical analysis of the beam than traditional chaotic models. In addition, the analysis will provide proof that a beam experiencing moderate deformations without stops or an initial axial force can exhibit chaotic motion. The model represents a simply-supported. Euler-Bernoulli beam subjected to a transverse load. The forcing function is sinusoidally distributed in space with an amplitude which also varies sinusoidally in time and is assumed to reach a maximum sufficient to allow nonlinearities associated with finite deformations to become important. During motion, even though displacements are large, the beam is assumed to attain only small strain levels and thus is assumed to be linearly elastic. The results indicate that for most levels of the forcing function the response of the beam is periodic. However, the steady state motion is not sinusoidal in time and in fact exhibits some bifurcated motions. At a certain level of the forcing amplitude, an asymmetry is observed and the periodicity of the motion breaks down as the beam experiences a period doubling cascade which culminates in a chaotic motion. The progression from periodic to chaotic motion is presented through a series of phase plane and Poincané plots, and physical variables such as bending moment are examined.     

开口薄壁梁的扭转理论与应用

以Vlasov薄壁构件理论为基础, 推导了开口薄壁构件一阶扭转理论. 该理论考虑了翘曲剪应力对截面转角的影响, 截面的转角分为自由翘曲转角和约束剪切转角, 在约束扭转中, St.Venant扭矩仅仅与自由翘曲转角有关, 而翘曲扭矩仅与约束剪切转角有关. 利用半逆解方法求出了约束扭转中薄壁构件的St.Venant扭矩表达公式; 依据能量方法, 建立了约束剪切转角和翘曲扭矩之间的关系, 并提出了翘曲剪切系数概念, 给出了一阶扭转理论的微分方程. 为了有效求解微分方程, 给出了求解微分方程的初参数法方程和相应的影响函数矩阵; 当St.Venant扭矩可以忽略时, 得到与一阶弯曲理论(Timoshenko梁理论)相似的一阶扭转理论简化形式. 最后利用算例证明了一阶扭转理论和简化理论的有效性.      

CLASSICAL AND HOMOGENIZED EXPRESSIONS FOR BUCKLING SOLUTIONS OF FUNCTIONALLY GRADED MATERIAL LEVINSON BEAMS

The relationship between the critical buckling loads of functionally graded material(FGM) Levinson beams(LBs) and those of the corresponding homogeneous Euler-Bernoulli beams(HEBBs) is investigated. Properties of the beam are assumed to vary continuously in the depth direction. The governing equations of the FGM beam are derived based on the Levinson beam theory, in which a quadratic variation of the transverse shear strain through the depth is included.By eliminating the axial displacement as well as the rotational angle in the governing equations,an ordinary differential equation in terms of the deflection of the FGM LBs is derived, the form of which is the same as that of HEBBs except for the definition of the load parameter. By solving the eigenvalue problem of ordinary differential equations under different boundary conditions clamped(C), simply-supported(S), roller(R) and free(F) edges combined, a uniform analytical formulation of buckling loads of FGM LBs with S-S, C-C, C-F, C-R and S-R edges is presented for those of HEBBs with the same boundary conditions. For the C-S beam the above-mentioned equation does not hold. Instead, a transcendental equation is derived to find the critical buckling load for the FGM LB which is similar to that for HEBB with the same ends. The significance of this work lies in that the solution of the critical buckling load of a FGM LB can be reduced to that of the HEBB and calculation of three constants whose values only depend upon the throughthe-depth gradient of the material properties and the geometry of the beam. So, a homogeneous and classical expression for the buckling solution of FGM LBs is accomplished.     

光滑节点插值法:计算固有频率下界值的新方法

将光滑节点插值法用于悬臂梁的静力学,并首次用于旋转柔性梁的频率分析. 采用梯度光滑技术,用线性插值形函数描述梁的位移场,求解4 阶微分方程. 在静力学分析中,将该方法所得梁中各点位移与假设模态法、有限元法及解析解的结果对比,可知该方法虽用简单的线性插值形函数描述梁的位移场,但精度却很高. 进一步研究表明,采用模态高于9 阶的假设模态法会使刚度阵条件数变差,导致结果发散. 在频率分析中,与有限元法、假设模态法和解析解对比,表明该方法一个重要特性:能提供固有频率的下界值,而有限元法和假设模态法只能提供固有频率的上界值,说明该方法结合有限元法在处理无解析解的问题时可以从上下界最大程度的逼近真实解,提高精度. 光滑节点插值法具有形函数结构简单、独立变量少且能提供固有频率下界值的特性,因此,具有较高的推广及应用价值.      

A NEW APPROACH TO FREE VIBRATION ANALYSIS OF A BEAM WITH A BREATHING CRACK BASED ON MECHANICAL ENERGY BALANCE METHOD

In this paper, a new approach to free vibration analysis of a cracked cantilever beam is proposed. By considering the effect of opening and closing the crack during the beam vibration, it is modeled as a fatigue crack. Also, local stiffness changes at the crack location are considered to be a nonlinear amplitude-dependent function and it is assumed that during one half a cycle, the frequencies and mode shapes of the beam vary continuously with time. In addition, by using the experimental tests, it is shown that the local stiffness at the crack location varies continuously between the two extreme values corresponding to the fully closed and the fully open cases of the crack. Then, by using the mechanical energy balance the dynamic response of the cracked beam is obtained at every time instant. The results show that for a specific crack depth, by approaching the crack location to the fixed end of the beam, more reduction in the fundamental frequency occurs. Furthermore, for a specific crack location, the fundamental frequency diminishes and the nonlinearity of the system increases by increasing the crack depth. In order to validate the results, the variations of the fundamental frequency ratio against the crack location are compared with experimental results.     

STEADY-STATE RESPONSE OF A TIMOSHENKO BEAM ON AN ELASTIC HALF-SPACE UNDER A MOVING LOAD

By introducing the equivalent stiffness of an elastic half-space interacting with a Timoshenko beam, the displacement solution of the beam resting on an elastic half-space subjected to a moving load is presented. Based on the relative relation of wave velocities of the half-space and the beam, four cases with the combination of different parameters of the half-space and the beam, the system of soft beam and hard half-space, the system of sub-soft beam and hard half-space, the system of sub-hard beam and soft half-space, and the system of hard beam and soft half-space are considered. The critical velocities of the moving load are studied using dispersion curves. It is found that critical velocities of the moving load on the Timoshenko beam depend on the relative relation of wave velocities of the half-space and the beam. The Rayleigh wave velocity in the half-space is always a critical velocity and the response of the system will be infinite when the load velocity reaches it. For the system of soft beam and hard half-space, wave velocities of the beam are also critical velocities. Besides the shear wave velocity of the beam, there is an additional minimum critical velocity for the system of sub-soft beam and hard half-space. While for systems of (sub-) hard beams and soft half-space, wave velocities of the beam are no longer critical ones. Comparison with the Euler-Bernoulli beam shows that the critical velocities and response of the two types of beams are much different for the system of (sub-) soft beam and hard half-space but are similar to each other for the system of (sub-) hard beam and soft half space. The largest displacement of the beam is almost at the location of the load and the displacement along the beam is almost symmetrical if the load velocity is smaller than the minimum critical velocity (the shear wave velocity of the beam for the system of soft beam and hard half-space). The largest displacement of the beam shifts behind the load and the asymmetry of the displacement along the beam increases with the increase of the load velocity due to the damping and wave radiation. The displacement of the beam at the front of the load is very small if the load velocity is larger than the largest wave velocity of the beam and the half space. The results of the present study provide attractive theoretical and practical references for the analysis of ground vibration induced by the high-speed train.