Application of fractional-derivative standard linear solid model to impact response of human frontal bone

The impact of a rigid body on a thin plate with a buffer is investigated in this paper. A buffer is assumed as a linear spring fractional derivative dashpot which exhibits the viscoelastic features. The fractional-derivative standard linear solid model is suggested for describing the shock interaction of the impactor with a circular elastic plate. We assume that a transient wave of transverse shear is generated in the plate and the reflected wave does not have sufficient time to interact with the plate before the impact process is completed. The ray method is used outside the contact spot, but the Laplace transform method is applied within the contact region. The time-dependence of the contact force is determined. A numerical example is carried out by considering crash scenarios in frontal impacts of the human head which could estimate brain injury risks.     

Viscoplastic analysis for direct impact of sports balls

A coefficient of restitution (COR) is used to represent viscoplastic dissipation of energy in the contact region of colliding bodies [6]. A model has been developed that distinguishes between rate-independent plastic and rate-dependent viscous energy dissipation during impact. A nondimensional contact force shape factor for compression α is introduced which can be measured experimentally. The measured α and COR are used to separate the energy dissipated during collision into part due to plastic deformation and another part due to viscous dissipation. In comparison to the viscoelastic compliance model, the viscoplastic compliance model approximates the force profile accurately, especially in terms of the maximum force that occurs during impact.     

飞机撞击建(构)筑物研究进展

飞机撞击重要建(构)筑物会导致灾难性后果。本文从试验研究、理论分析、数值模拟等3个方面对飞机撞击建(构)筑物的国内外研究现状,相关研究难点、需要注意的问题和研究方向及趋势进行总结,包括:缩比试验的系统和验证,飞机撞击力模型,撞击所致局部破坏计算公式,飞机和建(构)筑物的精细化建模,撞击所致振动特性,撞击荷载和火荷载对结构的耦合毁伤效应,一般模型和精细化模型、解耦和耦合方法以及不同数值模拟程序计算结果的对比分析等方面,以期为后续研究提供参考。     

强夯接触力和接触位移的碰撞分析求解

首次采用大变形动力接触有限元法分析了强夯对地基土的冲击碰撞过程 ,锤底接触力和夯沉量过程中直接求得 ,且能模拟夯锤与地基土的多次接触分离过程 ,适应于各种复杂的地基土条件和夯击能条件 ,克服了现有其它方法的不足和缺陷。文章最后给出了系列工程算例 ,验证了本文计算模型的合理型     

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

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

On the accurate determination of contact compliance for impact test modelling

To model the specimen interaction with supports during an impact test, simple formulas for indentation–contact force relation between a beam specimen and a rigid cylindrical indenter have been derived using a mixed analytical/numerical approach. Two types of boundary conditions for the specimen (i) support by a frictionless rigid foundation and (ii) conventional three-point bending have been considered. The first scheme of loading (the compression indentation test, CIT) is sometimes used for quasi-static estimation of the specimen–striker or specimen–support contact compliance instead of the second scheme, which more closely corresponds to the real loading conditions of the specimen during an impact test. It has been found that the indentation (and, therefore, the contact compliance) of the specimen loaded according to the first scheme is up to 19% higher than for the second one. A simple correction of the results of CIT, which allows to estimate the contact compliance accurately has been proposed. Approximate formulas for the linearized contact compliance have been derived for both schemes of loading using three different methods of linearization. The best result has been obtained by the method of the equality of work done by nonlinear and linearized contact forces. An example of modelling of the three-point-bending test using the computed contact stiffness of the anvil is presented.     

柔性梁斜碰撞问题的非线性动力学建模和实验研究

本文研究柔性梁点面斜碰撞问题。用Hertz接触模型处理法向撞击力,分别用Hertz切向接触模型和Coulomb摩擦力模型处理粘滞状态和滑动状态的摩擦力。从精确的应变与位移的关系出发,用绝对节点坐标法建立了柔性梁的动力学方程。为了准确地处理斜碰撞切向运动的复杂状态,提出滑动-粘滞切换的准则,在此基础上,设计了斜碰撞实验,数值对比了法向撞击力和法向速度的时间历程的仿真计算结果与实验结果,验证了Hertz理论在斜碰撞情况下的正确性。另一方面切向速度的实验与理论的结果对照表明滑动-粘滞切换准则的有效性。     

多体系统碰撞动力学中接触力模型的研究进展

接触碰撞行为作为大自然与多体系统中的常见现象,其接触力模型对于多体系统的碰撞行为机理研究与性能预测至关重要.静态弹塑性接触模型与考虑能量耗散的连续接触力模型是研究接触碰撞行为的两类不同方法,在多体系统碰撞动力学中存在诸多共性与差异.本文分别从上述两类接触模型的发展历程入手,详细介绍了两类模型的区别与联系.首先,根据阻尼项分母中是否含有初始碰撞速度将连续接触力模型分为黏性接触力模型与迟滞接触力模型,讨论了能量指数与Hertz接触刚度之间的关系,阐述了现有连续接触力模型在计算弹塑性材料接触碰撞行为时存在的问题.其次,着重介绍了分段连续的准静态弹塑性接触力模型(可连续从完全弹性转换到完全塑性接触阶段),分析了利用此类弹塑性接触力模型计算碰撞行为的技术特点.同时,以恢复系数为桥梁和借助线性化的弹塑性接触刚度,避免了Hertz刚度对弹塑性接触刚度的计算误差,根据碰撞前后多体系统的能量与动能守恒推导了弹塑性接触模型等效的迟滞阻尼因子.探索了连续接触力模型与准静态弹塑性接触力模型之间的内在联系,数值计算结果定量说明了人为阻尼项代表的能量耗散与弹塑性接触力模型中加卸载路径代表的能量耗散具有等效性.另外...     

Design and modeling of an endoscopic piezoelectric tactile sensor

The paper presents a proof-of-concept design of a tactile sensor capable of measuring compliance of a contact tissue/sensed object. The main objective of this study is to design and model a piezoelectric sensor capable of measuring the total applied force on the sensed object as well as its compliance. The sensor consists of a rigid and compliant cylindrical element. Determination of the compliance of sensed objects is based on the ratio of force experienced by the rigid cylinder to the total force applied to the sensor. To obtain this force ratio, a circular PVDF film is sandwiched between rigid cylinder and plate to measure the force applied on the rigid element and a rectangular PVDF film is sandwiched between the two base plates to measure the total force applied on the sensor. The detailed design of the senor was performed using finite element analysis. A prototype was fabricated and tested and it has been shown that good agreement exists between the finite element results and experimental values. The proposed sensor exhibits high force sensitivity and good linearity and offers the potential for future miniaturization in order to be integrated with the commercial endoscope graspers used in minimally invasive surgery.     

An analytical elastic-perfectly plastic contact model

A new formulation for elastic-perfectly plastic contact in the normal direction between two round surfaces that is solely based on material properties and contact geometries is developed. The problem is formulated as three separate domains: the elastic regime, mixed elastic–plastic behavior, and unconstrained (fully plastic) flow. Solutions for the force–displacement relationship in the elastic regime follow from Hertz’s classical solution. In the fully plastic regime, two well supported assumptions are made: that there is a uniform pressure distribution and there is a linear force–deflection relationship. The force–displacement relationship in the intermediate, mixed elastic–plastic regime is approximated by enforcing continuity between the elastic and fully plastic regimes. Transitions between the three regimes are determined based on empirical quantities: the von Mises yield criterion is used to determine the initiation of mixed elastic–plastic deformation, and Brinell’s hardness for the onset of unconstrained flow. Unloading from each of these three regimes is modeled as an elastic process with different radii of curvature based on the regime in which the maximum force occurred. Simulation results explore the relationship between the impact velocity and coefficient of restitution. Further comparisons are made between the model, experimental results found in the literature, and other existing elastic–plastic models. The new model is well supported by the experimental measurements of compliance curves for elastic–plastic materials and of coefficients of restitution from impact studies, and in elastic-perfectly plastic regimes is demonstrated to be more accurate than existing models found in the literature.     

Simulation of impacts in geartrains using different approaches

Gear hammering in diesel engines is a well-known phenomenon in geared drives, exhibiting not only noise but also influencing the performance and durability of diesel engines. Gear hammering is characterised by flanks in contact that lift off and cause impacts when the contact reestablishes, which induces high, sharp dynamic loads. The knowledge of these contact forces is very important for the design of gears. Since contact forces in meshing gears are extremely difficult and expensive to measure, the simulation of these forces plays an important role. Nowadays, these contact simulations are usually carried out within overall models of entire engines using commercial multibody programs that provide submodels for gear contacts, usually based on rigid-body models. However, to reduce inertia effects, gears in geartrains are often designed with very thin bodies, whose elastic compliance influences the contact behaviour to a large extent. For a closer insight into the dynamic behaviour, and especially the influence of thin gear bodies during impact, a typical gear pairing is selected and impacts between one tooth pair are investigated for different boundary and initial conditions with three different models. Besides a multibody model, similar to those used in commercial multibody programs, a fully nonlinear finite-element model and a modally reduced model in combination with a local force law is used. The results of the different approaches are benchmarked in terms of accuracy and numerical effort.     

A general method for impact dynamic analysis of a planar multi-body system with a rolling ball bearing joint

Impact affects the dynamic characteristics of mechanical multi-body systems and damages those rotating parts, such as the joint rolling element bearings, which are high-precision, defect intolerant components. Based on multi-body dynamic theory, Hertzian contact theory, and a continuous contact model, this study proposed a modelling method that can describe the dynamic behaviour of planar mechanical multi-body systems containing a rolling ball bearing joint under impact. In this method, the rigid bodies and bearing joint were connected according to their joint force constraints; the impact constraint between the multi-body system and the target rigid body was constructed using a continuous contact force model. Based on this method, the reflection relationship between the external impacts of the mechanical multi-body system and the variation law governing the dynamic load on the rolling bearing joint were revealed. Subsequently, an impact multi-body system, which was composed of a sliding–crank mechanism containing a rolling ball bearing joint and the target rigid body with an elastic support, was analysed to explore the dynamic response of such a complex discontinuous dynamic system andthe relevant relationship governing the dynamic load on the rolling bearing joint. In addition, a multi-body dynamic simulation software was used to build a virtual prototype of the impact slider–crank system. Compared with the theoretical model, the prototype had an additional deep groove ball bearing. That is to say, the prototype model took account of the specific geometric structural characteristics and the complex contact relationship of the inner and outer races, rolling balls, and bearing cage. Finally, the effectiveness of the theoretical method proposed in this study was verified by comparative analysis of the results. The results suggested that the external impact of a mechanical multi-body system was prone to induce sudden changes in the equivalent reaction force on its bearing joint and the dynamic load carried on its rolling balls. This study provided an effective method for exploring the distribution characteristics of dynamic loads on rolling ball bearing joints under working impact load conditions. Moreover, it offered support for the parameter optimisation of geometric structure, performance evaluation, and dynamic design of the rolling ball bearings.     

基于LuGre摩擦模型的接触约束法旋转柔性梁斜碰撞研究

基于接触约束法和LuGre摩擦模型对在重力场作用下作大范围旋转运动的柔性梁系统和斜坡发生含摩擦斜碰撞的动力学问题进行研究. 首先运用刚柔耦合的多体系统动力学理论对大范围运动的柔性梁进行离散化和动力学建模, 在碰撞时采用冲量动量法求出跳跃速度, 其次在法向上引入接触约束求解出碰撞力, 在切向上采用LuGre摩擦模型分两种方式求解摩擦力, 第一种是在滑动时摩擦力由摩擦系数和碰撞力计算得出, 黏滞状态下引入切向约束计算拉格朗日乘子反应实际摩擦力, 根据黏滞/滑动切换判断计算出碰撞过程摩擦力(与Coulomb摩擦模型计算摩擦力一致); 第二种根据LuGre摩擦模型摩擦系数和法向碰撞力计算其摩擦力, 从而在碰撞时无需黏滞/滑动切换, 采用相同的摩擦力计算公式. 通过与Coulomb摩擦模型对比发现, LuGre摩擦模型描述碰撞切向摩擦过程更精确, LuGre摩擦模型黏滞时建立约束方程和碰撞采用统一的摩擦力公式这两种建模方式描述的斜碰撞动力学特性没有区别, 进而说明采用法向接触约束和LuGre摩擦模型具有满足碰撞非嵌入情况、避免黏滞/滑动切换、描述摩擦力相对准确的优势.      

Characteristics of local compliance of an elastic body under a small punch indented into the plane part of its boundary

An asymptotic solution of the contact problem of an elastic body indented (without friction) by a circular punch with a flat base is obtained under the assumption of a small relative size of the contact zone. The resulting formulas involve integral characteristics of the elastic body, which depend on its shape, dimensions, fixing conditions, Poisson's ratio, and location of the punch center. These quantities have the mechanical meaning of the coefficients of local compliance of the elastic body. Relations that, generally, reduce the number of independent coefficients in the asymptotic expansion are obtained on the basis of the reciprocal theorem. Some coefficients of local compliance at the center of an elastic hemisphere are calculated numerically. The asymptotic model of an elastic body loaded by a point force is discussed.     

低速冲击激励下嵌入黏弹性阻尼芯层的纤维金属混杂层合板动态响应预测模型

本文首次从解析角度建立了低速冲击激励下嵌入黏弹性阻尼芯层的纤维金属混杂层合板动态响应预测模型. 首先,结合经典层合板理论和冯$\cdot$卡门假设,建立了嵌入黏弹性芯层的纤维金属混杂层合板弹性损伤本构关系. 然后,将层合板受冲击时的变形分成接触和拉伸两个区域,在接触区域内,对金属层采用 Von Mises 失效准则,纤维层采用 Tsai-Hill 失效准则和对黏弹性层采用指数 Drucker-Prager 失效准则判断层合板损伤情况. 考虑不同材料层对冲击动态响应的贡献来修正两个变形区域的位移公式,进而计算结构因弹性变形产生的应变能,以及接触区域因塑性变形消耗的能量,实现每次失效事件发生后各层材料的能量、位移和冲击接触力的理论求解,并给出了结构动态响应分析的具体流程图. 最后,以嵌入 Zn33 黏弹性芯层的 TA2 钛合金混杂 T300 碳纤维/树脂层合板为研究对象,开展落锤冲击实验. 验证结果表明,理论预测与测试获得的冲击接触力、位移响应以及冲击载荷-位移曲线吻合较好,且关注的峰值点计算误差最大不超过 9%,进而验证了所提出的理论模型的有效性.      

基于ANSYS/LS-DYNA的钢板混凝土墙冲击实验的有限元分析

运用经典的显式非线性动力分析软件ANSYS/LS-DYNA,分析了1/7.5缩尺飞机模型撞击钢板混凝土墙的冲击实验。选用两种不同的混凝土材料本构模型(Winfrith模型、CSCM模型)模拟混凝土的非线性破坏过程,将钢板混凝土墙的破坏模式以及飞机模型的残余速度等与实验结果进行了对比。结果表明,有限元分析结果与实验吻合较好,且Winfrith材料模型能够更好地模拟混凝土的大应变、高应变率的非线性性能,验证了钢板混凝土墙和飞机材料本构模型的选取以及整个分析方法的适用性和有效性。     

A DYNAMIC RESPONSE PREDICTION MODEL OF FIBER-METAL HYBRID LAMINATED PLATES EMBEDDED WITH VISCOELASTIC DAMPING CORE UNDER LOW-VELOCITY IMPACT EXCITATION 1)

本文首次从解析角度建立了低速冲击激励下嵌入黏弹性阻尼芯层的纤维金属混杂层合板动态响应预测模型. 首先,结合经典层合板理论和冯$\cdot$卡门假设,建立了嵌入黏弹性芯层的纤维金属混杂层合板弹性损伤本构关系. 然后,将层合板受冲击时的变形分成接触和拉伸两个区域,在接触区域内,对金属层采用 Von Mises 失效准则,纤维层采用 Tsai-Hill 失效准则和对黏弹性层采用指数 Drucker-Prager 失效准则判断层合板损伤情况. 考虑不同材料层对冲击动态响应的贡献来修正两个变形区域的位移公式,进而计算结构因弹性变形产生的应变能,以及接触区域因塑性变形消耗的能量,实现每次失效事件发生后各层材料的能量、位移和冲击接触力的理论求解,并给出了结构动态响应分析的具体流程图. 最后,以嵌入 Zn33 黏弹性芯层的 TA2 钛合金混杂 T300 碳纤维/树脂层合板为研究对象,开展落锤冲击实验. 验证结果表明,理论预测与测试获得的冲击接触力、位移响应以及冲击载荷-位移曲线吻合较好,且关注的峰值点计算误差最大不超过 9%,进而验证了所提出的理论模型的有效性.     

锻压操作机夹持长棒料随动过程动态受力分析

研究了锻压操作机夹持长棒料锻件锻压随动过程中夹钳的受力状态。利用拉格朗日方程及有限元方法建立了包含操作机及其锻件的总体动力学方程。操作机系统建立为端部悬挂质量弹簧的Euler-Bernoulli梁模型,随动过程中因为梁上节点位移激励引起操作机系统的瞬态振动。计算结果与LS-DYNA显式有限元模型结果进行了比对,吻合良好。结果表明,对于长棒料锻件,液压机锻压的位置对操作机的动态受力有着十分显著的影响。     

基于Hopkinson杆和人工神经网络的三轴冲击力传感器同步解耦标定方法

三分量冲击力载荷的同步激励与输入输出间的精准建模是三轴冲击力传感器标定所面临的主要挑战。为了实现对三轴冲击力传感器的有效标定,使其能够准确测量空间中的三维冲击力载荷。首先,基于Hopkinson杆与矢量分解原理建立了一种高幅值（104 N量级）、窄脉宽（10?4 s量级）可计量三分量冲击力载荷的同步激励方法,实现了对三轴冲击力传感器的同步加载。然后,基于最小二乘原理与矩阵微分构建了三轴冲击力传感器的线性标定模型,并通过改变子弹结构与冲击气压揭示了线性解耦标定模型中传感器主灵敏度系数与轴间耦合灵敏度系数并非固定常数而均与冲击力载荷脉冲构型（幅值、脉宽）相关的冲击特性。最后,将能够反映载荷构型信息的传感器各轴输出电压脉冲的幅值与脉宽作为影响因素,并以神经元的形式添加到人工神经网络（artificial neural network, ANN）的输入层,建立了基于ANN的三轴冲击力传感器输出电压与输入载荷间的代理模型,实现了数据驱动的三轴冲击力传感器非线性解耦标定。结果表明,相对最小二乘模型,ANN标定精度更高,采用ANN进行三轴冲击力传感器标定具有可行性和有效性。