Static analysis of an infinite beam resting on a tensionless Pasternak foundation

This paper addresses the static response of an infinite beam supported on a unilateral (tensionless) two-parameter Pasternak foundation and subjected to complex transverse loads, including self weight. The transfer displacement function method (TDFM) is employed to determine the initially unknown lengths that remain in contact. In contrast to a Winkler Foundation System (WFS), the lift-off points in a PFS (Pasternak Foundation System) are not necessarily at zero displacement but may be determined sequentially through considering the compatibility conditions at the junctions of contact and non-contact segments. After the response of the whole system including the beam and foundation is expressed through the displacement constants of the initial segment, the contact problem is reduced to two nonlinear algebraic equations with two unknowns. The foundation reactions and the internal actions of the beam may also be determined from the displacement response of the system. Two simple cases are solved to illustrate the influence of the foundation stiffness factors and finally, a third example of a beam with several contact segments is presented to demonstrate the application of the TDFM.     

高速荷载下多孔饱和地基的动力响应

研究高速荷载作用下梁与多孔饱和半空间的动力响应。由Fourier变换求解多孔饱和固体的动力基本方程，根据梁与半空间的接触条件得出多孔饱和半空间上梁的垂直位移的表达式。文中的数值算例考虑了荷载移动速度对梁的动力位移的影响，并与相应的弹性半空间问题作了对比。从算例中可以发现荷载移动速度对动力位移有很大的影响，当移动速度与半空间的表面波速相近时，地面会当产生很大的振动，同时还发现当速度大于介质的剪切波速时，多孔饱和半空间上梁的动力响应与弹性半空间上梁的动力响应有很大的差别。     

Experimental and numerical investigations of the responses of a cantilever beam possibly contacting a deformable and dissipative obstacle under harmonic excitation

In this paper, the dynamics of a cantilever beam subjected to harmonic excitations and to the contact of an obstacle is studied with the help of experimental and numerical investigations. The steel flexible structure is excited close to the free end with a shaker and may come into contact with a deformable and dissipative obstacle. A technique for modeling contact phenomena using piece-wise linear dynamics is applied. A finite-dimensional modal model is developed through a Galerkin projection. Concentrated masses, dampers and forces are considered in the equations of motion in such a way that the boundary conditions are those of a cantilever beam. Numerical studies are conducted by assuming finite-time contact duration to investigate the frequency response of the impacted beam for different driving frequencies. Experimental results have been extrapolated through a displacement laser sensor and a load cell. The comparison between numerical and experimental results show many qualitative and quantitative similarities.The novelty of this paper can be synthetized in (a) the development of experimental results that are in good agreement with the numerical implementation of the introduced model; (b) the development of a comprehensive contact model of the beam with an unilateral, deformable and dissipative obstacle located close to the tip; (c) the possibility of accounting for higher modes for the cantilever beam problem, and hence of analyzing how the response varies when moving the excitation (and/or the obstacle) along the beam, and of investigating the effect of the linearly elastic deformability of the built‐in end of the beam; (d) an easy and intuitive solution to the problem of accounting for spatially singular masses, dampers, springs and forces in the motion equations; (e) the possibility of accounting for finite gap and duration of the contact between beam and obstacle.     

含黏弹性夹芯FG-GRC后屈曲梁的低速冲击响应

研究了含黏弹性夹芯的功能梯度石墨烯增强复合材料(functionally graded graphene reinforced composite, FG-GRC)后屈曲梁在低速跌落冲击下的跳跃振荡行为.采用修正Halpin-Tsai细观模型预测FG-GRC的材料宏观属性.使用赫兹点接触模型确定冲击器和梁之间的接触力.提出了考虑轴向预应力的复合材料层本构关系和阻尼层的Kelvin型黏弹性本构.通过一种广义高阶剪切变形锯齿梁模型建立夹芯梁的非线性位移场. 基于Hamilton 能量变分原理, 推导了动力学控制方程组. 通过两步分析,首先获得弹性后屈曲平衡路径作为冲击问题的初始状态. 随后, 结合四阶龙格库塔法,拓展了两步摄动-伽辽金法计算接触力的时程曲线以及后屈曲梁的位移时程曲线.研究了后屈曲梁在单次和两次撞击下双稳态大幅振荡过程的动力学特征.讨论了轴向载荷、冲击速度、黏弹性阻尼特性、冲击器材料等因素对于碰撞接触力以及后屈曲梁动力响应的影响规律.结果表明, 接触力仅对冲击速度较为敏感,一定的结构碰撞参数设计可以在接触力变化不大的情况下,使得后屈曲梁由单势能阱运动转变为双阱大幅振荡.      

考虑局部切向约束接触问题的直接刚度法

针对工程中大量存在的切向滑移受到约束的接触问题,提出了基于Lagrange乘子的点-面及点-点接触直接刚度法,该方法生成的接触协调条件可以直接组装到结构刚度阵中,从而可直接用于考虑存在结构几何非线性及材料非线性的接触问题分析中。采用这一算法进行了考虑拱坝横缝张合效应的地震响应分析,结果表明考虑切向滑移约束对接触中的张开度有明显影响。     

Soft-impact dynamics of deformable bodies

Systems constituted by impacting beams and rods of non-negligible mass are often encountered in many applications of engineering practice. The impact between two rigid bodies is an intrinsically indeterminate problem due to the arbitrariness of the velocities after the instantaneous impact and implicates an infinite value of the contact force. The arbitrariness of after-impact velocities is solved by releasing the impenetrability condition as an internal constraint of the bodies and by allowing for elastic deformations at contact during an impact of finite duration. In this paper, the latter goal is achieved by interposing a concentrate spring between a beam and a rod at their contact point, simulating the deformability of impacting bodies at the interaction zones. A reliable and convenient method for determining impact forces is also presented. An example of engineering interest is carried out: a flexible beam that impacts on an axially deformable strut. The solution of motion under a harmonic excitation of the beam built-in base is found in terms of transverse and axial displacements of the beam and rod, respectively, by superimposition of a finite number of modal contributions. Numerical investigations are performed in order to examine the influence of the rigidity of the contact spring and of the ratio between the first natural frequencies of the beam and the rod, respectively, on the system response, namely impact velocity, maximum displacement, spring stretching and contact force. Impact velocity diagrams, nonlinear resonance curves and phase portraits are presented to determine regions of periodic motion with impacts and the appearance of chaotic solutions, and parameter ranges where the functionality of the non-structural element is at risk.     

考虑胎路多点接触的电动汽车?路面耦合系统振动分析

轮毂电机驱动电动汽车的簧下质量大, 使得轮胎动载荷增加, 且电机激励进一步加剧车轮振动. 同时, 轮胎与路面单点接触的简化模型, 其动力学计算结果与实际存在差别. 鉴于此, 考虑电机的电磁激励、胎路多点接触和非线性地基, 建立了电动汽车?路面系统机电耦合动力学模型, 通过Galerkin法推导了非线性地基梁的垂向振动, 利用积化和公式推导了非线性地基梁中非线性项积分的精确表达式, 提出了路面截断阶数选取的简易方法, 并通过路面位移响应的收敛性进行了验证. 在此基础上, 研究了胎路多点接触、非线性地基、电机激励、车速、路面不平顺幅值等对路面及车辆响应的影响. 结果表明, 非线性地基及多点接触对车辆响应的影响中, 轮胎动载荷的影响最大, 车身加速度和悬架动挠度的影响较小, 且考虑电机激励时, 二者对车辆响应的影响显著增大. 从对路面响应的影响看, 电机激励的影响最大, 非线性地基的影响次之, 多点接触的影响较小. 所建模型及研究方法可为电动汽车的垂向动力学分析提供一种新思路.      

考虑尺度效应的微梁静态弯曲数值分析

基于修正偶应力理论及表面弹性理论,本文提出了一种新的双曲线剪切变形梁模型,用于均匀微尺度梁的静态弯曲分析。该理论可以直接利用本构关系获得横向剪切应力,满足梁顶部和底部的无应力边界条件,避免了引入剪切修正因子。根据广义Young-Laplace方程建立了梁的内部与表面层的应力连续性条件,单一的变量场可以描述梁的位移模式。通过在位移场中考虑表面层厚度以及表面层的应力连续条件,可以使新模型能够更准确地预测微尺寸和表面能相关的尺度效应。通过Hamilton原理推导出了梁的控制方程和边界条件。应变能除了考虑经典弹性理论,还要考虑微结构效应和表面能。Navier-type的解析解适用于简支边界条件,而基于拉格朗日插值的微分求积法(DQEM)可以研究在不同边界条件下的力学响应。把该数值解与Navier方法得出的解析解作了对比,得出:微尺度梁在考虑表面能或微尺寸效应、不同载荷和梁高变化下的响应一致;当不考虑微结构相关性和表面能效应时,该模型退化为经典的欧拉梁模型。     

Transient wave analysis of a cantilever Timoshenko beam subjected to impact loading by Laplace transform and normal mode methods

This study applies two analytical approaches, Laplace transform and normal mode methods, to investigate the dynamic transient response of a cantilever Timoshenko beam subjected to impact forces. Explicit solutions for the normal mode method and the Laplace transform method are presented. The Durbin method is used to perform the Laplace inverse transformation, and numerical results based on these two approaches are compared. The comparison indicates that the normal mode method is more efficient than the Laplace transform method in the transient response analysis of a cantilever Timoshenko beam, whereas the Laplace transform method is more appropriate than the normal mode method when analyzing the complicated multi-span Timoshenko beam. Furthermore, a three-dimensional finite element cantilever beam model is implemented. The results are compared with the transient responses for displacement, normal stress, shear stress, and the resonant frequencies of a Timoshenko beam and Bernoulli–Euler beam theories. The transient displacement response for a cantilever beam can be appropriately evaluated using the Timoshenko beam theory if the slender ratio is greater than 10 or using the Bernoulli–Euler beam theory if the slender ratio is greater than 100. Moreover, the resonant frequency of a cantilever beam can be accurately determined by the Timoshenko beam theory if the slender ratio is greater than 100 or by the Bernoulli–Euler beam theory if the slender ratio is greater than 400.     

On the asymptotic boundary conditions of an anisotropic beam via virtual work principle

In this research, an efficient and effective method is proposed to derive the boundary conditions of an anisotropic beam in the asymptotic sense. We first set up the constrained virtual work by introducing the Lagrange multiplier on the displacement prescribed boundary. The macroscopic beam and microscopic cross-section equations with the boundary conditions are simultaneously obtained by taking the asymptotic expansion on the displacement vector. In this way, the three-dimensional characteristics of the beam are asymptotically smeared into the macroscopic beam equations and the beam boundary conditions. The boundary conditions obtained are then compared to those from the decay analysis method. The beam bending slope boundary condition obtained in the frame work of variational principle is different from the well-known average condition. This new boundary condition is more accurate than the average one for a sandwich beam. This is further demonstrated and discussed via the examples of a cantilever beam loaded at the end.     

Design of nonlinear springs for attaining a linear response in gap-closing electrostatic actuators

Gap-closing electrostatic actuators are inherently nonlinear and their dynamic range is often limited by the pull-in instability. To overcome this, we propose a nonlinear spring that counteracts the nonlinear effects of electrostatic attraction. The nonlinear spring is designed to extend the stable range of the actuator and to enforce a linear electromechanical response. We present a method for designing elastic springs with monotonically increasing stiffness. The mechanism we propose is effective shortening of a straight clamed-guided beam flexure, by wrapping it over a cam. We consider two specific cases. The first case assumes the wrapped section of the beam flexure fully conforms to the cam shape. The second case assumes that there is a single contact point between the beam flexure and the cam. To validate the concept we have designed and measured the response of a nonlinear spring with a prescribed force–displacement law. Experimental measurements of a macro-scale spring are in good agreement with the model predictions.     

Shear deformation effect in non-linear analysis of composite beams of variable cross section

In this paper the non-linear analysis of a composite Timoshenko beam with arbitrary variable cross section undergoing moderate large deflections under general boundary conditions is presented employing the analog equation method (AEM), a BEM-based method. The composite beam consists of materials in contact, each of which can surround a finite number of inclusions. The materials have different elasticity and shear moduli with same Poisson's ratio and are firmly bonded together. The beam is subjected in an arbitrarily concentrated or distributed variable axial loading, while the shear loading is applied at the shear center of the cross section, avoiding in this way the induction of a twisting moment. To account for shear deformations, the concept of shear deformation coefficients is used. Five boundary value problems are formulated with respect to the transverse displacements, the axial displacement and to two stress functions and solved using the AEM. Application of the boundary element technique yields a system of non-linear equations from which the transverse and axial displacements are computed by an iterative process. The evaluation of the shear deformation coefficients is accomplished from the aforementioned stress functions using only boundary integration. Numerical examples are worked out to illustrate the efficiency, the accuracy, the range of applications of the developed method and the influence of the shear deformation effect.     

不确定边界条件下悬臂梁裂纹参数识别方法研究

基于Bernoulli-Euler梁振动理论,以等效弹簧来模拟裂纹引起的局部软化效应和由非完全固支边界条件引起的转角效应。推导了悬臂梁在不确定边界条件下确定其振动频率的特征方程,直接利用该特征方程,提出一种有效估计裂纹参数的优化方法,通过计算测量频率和理论频率之间的误差目标函数最小化即可识别裂纹参数-裂纹位置和深度。最后,应用两个实例-理想固支边界条件下和非完全固支边界条件下的悬臂梁实验来说明本文方法的有效性。实验结果表明:只需梁结构前三阶频率即可识别裂纹位置和深度。对于理想边界条件下的裂纹参数识别,在测量频率存在小误差情况下,该方法仍能给出比较满意的结果,对于非完全固支边界条件下的裂纹参数识别,利用本文方法能得到比Narkis的方法更精确的裂纹位置识别结果。同时本文方法还能给出比较满意的裂纹深度识别结果。     

A SEMI-ANALYTICAL METHOD FOR THE VIBRATION OF AND SOUND RADIATION FROM A TWO-DIMENTIONAL BEAM-STIFFENED PLATE

A semi-analytical method based on space harmonics to investigate the vibration of and sound radiation from an infinite,fluid-loaded plate is presented.The plate is reinforced with two sets of orthogonally and equally spaced beam stiffeners,which are assumed to be line forces.The response of the stiffened plate to a convected harmonic pressure in the wave-number space is obtained by adopting the Green’s function and Fourier transform methods.Using the boundary conditions and space harmonic method,we establish the relationship between the stiffener forces and the vibration displacement of the plate.In this paper,the stiffener forces are expressed in terms of harmonic amplitudes of the plate displacement,which are calculated by using a numerical reduction technique.Finally,the Fourier inverse transform is employed to find expressions of the vibration and sound radiation in physical space.Agreements with existing results prove the validity of this approach and more numerical results are presented to show that this method converges rapidly.     

表面裂纹问题瞬态响应的边界元分析

本文将拉氏变换-边界元法用于表面裂纹问题的瞬态响应分析。文中讨论了拉氏反演参数的选择和动态应力强度因子的计算方法。作为程序的考核和离散方案的选择,分别地计算了水平柱体一端固定、另一端受p(t)=poH(t)载荷时的位移响应和具有贯穿裂纹的厚板在两种离散方案时的动态应力强度因子响应。最后,还计算了若干载荷工况的半圆表面裂纹板应力强度因子的瞬态响应,获得了有效的数值结果。     

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.     

Inverse Estimation of Cohesive Fracture Properties of Asphalt Mixtures Using an Optimization Approach

Tensile cracking in asphalt pavements due to vehicular and thermal loads has become an experimental and numerical research focus in the asphalt materials community. Previous studies have used the discrete element method (DEM) to study asphalt concrete fracture. These studies used trial-and-error to obtain local fracture properties such that the DEM models approximate the experimental load-crack mouth opening displacement response. In the current study, we identify the cohesive fracture properties of asphalt mixtures via a nonlinear optimization method. The method encompasses a comparative investigation of displacement fields obtained using both digital image correlation (DIC) and heterogeneous DEM fracture simulations. The proposed method is applied to two standard fracture test geometries: the single-edge notched beam test, SE(B), under three-point bending, and the disk-shaped compact tension test, DC(T). For each test, the Subset Splitting DIC algorithm is used to determine the displacement field in a predefined region near the notch tip. Then, a given number of DEM simulations are performed on the same specimen. The DEM is used to simulate the fracture of asphalt concrete with a linear softening cohesive contact model, where fracture-related properties (e.g., maximum tensile force and maximum crack opening) are varied within a predefined range. The difference between DIC and DEM displacement fields for each set of fracture parameters is then computed and converted to a continuous function via multivariate Lagrange interpolation. Finally, we use a Newton-like optimization technique to minimize Lagrange multinomials, yielding a set of fracture parameters that minimizes the difference between the DEM and DIC displacement fields. The optimized set of fracture parameters from this nonlinear optimization procedure led to DEM results which are consistent with the experimental results for both SE(B) and DC(T) geometries.     

On the displacement potential solution of plane problems of structural mechanics with mixed boundary conditions

The present paper describes the advancement of displacement potential approach in relation to solution of plane problems of structural mechanics with mixed mode of boundary conditions. Both the conditions of the plane stress and the plane strain are considered for analyzing the displacement and stress fields of the structural problem. Using the finite difference technique based on the present displacement potential approach for the case of the plane stress and the plane strain conditions, firstly an elastic cantilever beam subjected to a pure shear at its tip is solved and these two solutions (plane stress and plane strain) are compared with Timoshenko and Goodier cantilever beam bending solutions (Theory of elasticity, 2nd edn. McGraw-Hill, New York, 1951); secondly the above-mentioned displacement potential approach for the case of the plane stress and the plane strain conditions are applied to solve a one-end fixed square plate subjected to a combined loading at its tip. Effects of plane stress and plane strain on the elastic field of the plate are discussed in a comparative fashion. Limitations of Timoshenko and Goodier cantilever beam bending solutions (Theory of elasticity, 2nd edn. McGraw-Hill, New York, 1951) over the displacement potential approach for the case of the plane stress and the plane strain conditions are not only discussed but also the superiority of the present displacement potential approach for the case of the plane stress and the plane strain conditions are reflected in the present research work.     

Energy method solution for the postbuckling response of an axially loaded bilaterally constrained beam

Bistable and multistable structures have shown great usefulness in many applications such as MEMS actuation and energy harvesting. Bistability of structures can be achieved through buckling. Confining a buckled beam between two lateral constraints allows it to buckle into higher modes as the axial load increases. This paper presents a theoretical study of the postbuckling response of a bilaterally constrained elastica subjected to gradually increased axial load. Equilibrium states are determined using an energy method. Under small deformation assumptions, the total potential energy is minimized under the defined constraints. The presented model allows for an accurate representation of the flatting behavior and the increase in the length of contact areas with the lateral constraints before the sudden snapping between equilibrium states. Mode transitions are manifested by jumps in the response curves. Previously developed models based on geometry and symmetries overestimate the required forces for higher equilibrium modes and do not match experimental observations. Results are validated with experimental force–displacement measurements under both force- and displacement-control. The kinetic energy released during buckling mode transitions is determined by a dynamic analysis.