研究柔性体撞击问题的子结构离散方法

本文用子结构方法研究刚性小球和均质柔性杆的纵向撞击以及和均质柔性梁的横向撞击问题,导出了用模态坐标表示的动力学方程,通过对刚性小球和柔性杆的纵向撞击的仿真教育处发现,在总单元数相同的情况下,取8个子结构比较合理,在此基础上对刚性小球和均质柔性梁的横向撞击的特性进行研究,发现撞击力在变化过程中会产生上下波动,当梁的弹性模量增加时,撞击力增大,撞击时间缩短。     

复合材料T型接头冰雹高速撞击损伤的数值模拟

利用高速空气炮进行了冰雹撞击复合材料T型接头的实验,研究了不同撞击速度下结构的损伤情况。同时,采用光滑质点流体动力学方法与黏聚区模型相结合的方法,对冰雹撞击复合材料T型接头进行了数值模拟。通过与实验结果的对比,验证了数值模拟模型的有效性。在此基础上,利用数值模型研究了影响复合材料T型接头冰雹撞击损伤的各种因素。计算结果表明:冰雹撞击对复合材料T型接头造成的损伤主要为分层损伤;冰雹的撞击速度、尺寸、入射角等,都对T型接头的损伤程度有很大影响;T型接头的分层面积与冰雹的撞击能量之间呈近似线性关系,分层面积随着撞击能量的增大而增大;冰雹的入射角越大,分层面积与撞击力峰值也越大。     

圆柱杆侵彻过程及应力波传播分析

研究了初始侵彻速度为１０００ｍ／ｓ～２０００ｍ／ｓ的圆锥杆对半无限厚水泥靶的垂直侵彻过程及其应力波传播。针对初始侵彻速度范围内的撞击特点，发展了计算撞击力的刚体－流体撞击模型和计算靶体变形的法向膨胀理论。模拟动量定理和功能守恒定律，建立了圆锥杆垂直侵彻半元限厚水泫靶时的动力学模型。结非线性动力学方程的数值解，确定了撞击力。在些基础上，根据牛顿第二定律建立了圆锥杆的波动方程，研究了应力波在圆锥中的传     

柱形聚碳酸酯弹丸撞击刚性靶的实验研究

用Ф22×64mm的聚碳酸酯(PC)圆柱形弹丸,以172m/s~234m/s的速度撞击刚性靶板,用幅频为4.05×104幅/秒的高速数字相机记录了PC弹丸撞击时的变形过程,用PVDF测力薄膜测定了弹丸头部与刚性靶表面碰靶过程的应力-时间曲线,获得了可供建立分析模型参考的实验数据。采用Syomds-Cowper过应力模型对PC材料应变率敏感性进行了初步分析,根据经典的Taylor理论计算了PC的名义应变率和相应的动态流动应力,其结果与实验数据基本一致。     

商用飞机撞击核电站屏蔽厂房数值模拟

为研究大型商用飞机撞击核电站屏蔽厂房的毁伤特性,建立了Boeing 767商用飞机和双钢板混凝土屏蔽厂房的有限元模型,模拟两者的相互撞击作用过程,得到飞机不同部位(机身、发动机)对屏蔽厂房的撞击力时程曲线,对每一个作用部分给出简化撞击力曲线和作用面积,确定了撞击力分布形式。结果分析表明:飞机轴向网格尺寸对撞击力影响较大;屏蔽厂房被撞击部位变形明显,其他区域变形较小;撞击速度对撞击作用时间影响较小,而对结构响应位移影响很大,撞击力合力随着撞击速度的降低迅速下降。     

撞击荷载作用下单层椭圆抛物面网壳的动力稳定分析

利用ANSYS/LS-DYNA 有限元程序,考虑几何非线性和材料非线性效应,开展单层椭圆抛物面网壳结构在撞击荷载作用下的动力稳定性研究.分析了该类网壳在以一定高度的自由落体撞击作用下的撞击力响应(撞击力、节点位移、杆件轴力和能量)及网壳失稳模态,进一步根据产生的动力响应对其在撞击作用下的弹塑性动力稳定问题进行了判定.结果表明:网壳结构在撞击作用下失稳时,撞击力时程曲线形状基本上呈陡峭的三角形脉冲荷载形式,其最大峰值和脉宽与撞击冲量及网壳所处变形阶段的刚度性能相关;特征节点的位移响应突然增大,失稳部位杆件的轴力大幅减小,结构的变形能和总能量也突然增大,整个结构变形形状发生突变;节点的残余位移、撞击力和持时等特征参数都产生一致性突变.     

Q235钢板对半球形头弹抗侵彻特性

利用轻气炮进行了半球形头杆弹正撞击单层板和等厚接触式三层板的实验, 得到了这两种结构靶体的初始-剩余速度曲线以及弹道极限。采用ABAQUS/EXPLICIT数值模拟软件对杆弹撞击金属板的过程进行了数值模拟研究, 通过对比数值模拟和实验结果, 验证了数值模拟材料模型和参数的有效性。研究了靶体结构对抗侵彻特性的影响, 并分析了弹体对靶体的撞击过程。研究结果表明:多层板的弹道极限高于等厚单层板。单层板主要失效模式为剪切, 而多层板的主要失效模式为整体的蝶形变形和局部的盘式隆起。对于多层板, 靶板具体的失效模式与其在靶中位置相关。     

弹体高速侵彻钢筋混凝土靶试验研究

为研究结构弹体对钢筋混凝土靶的高速侵彻破坏效应,利用口径35 mm弹道炮开展了1 030～1 520 m/s速度范围内的高速侵彻试验,获得了弹体的撞击速度、破坏形态、剩余长度、剩余质量和靶体中的侵彻深度及成坑尺寸等试验数据,分析了侵彻深度和侵彻机理随速度的变化关系。结果表明:在1 030～1 390 m/s的速度范围内,弹体头部磨蚀,磨蚀程度随侵彻速度增加而加剧,侵彻深度随撞击速度近似线性增大;撞击速度在1 390～1 480 m/s范围内,弹体头部严重磨蚀,侵彻深度随撞击速度增加而减小;撞击速度大于1 480 m/s后,弹体严重破碎,侵彻深度急剧下降。针对结构弹体高速侵彻过程中的破坏特点,将侵彻速度划分为刚体侵彻区、准刚体侵彻区、侵蚀体侵彻区和破碎体侵彻区,可为钻地弹结构设计提供参考。     

高温作用下钢管混凝土构件侧向撞击性能

通过耦合ABAQUS有限元软件中的隐式静态分析和显示动态分析，提出钢管混凝土构件在火灾与撞击联合作用下的数值计算方法，分别对已有钢管混凝土构件的温度场试验、火灾下轴向撞击试验和常温下侧向撞击试验进行数值模拟，以验证本文方法的合理性。在此基础上建立了钢管混凝土构件在不同温度下的侧向撞击有限元模型，分别对不同温度下的挠度和撞击力时程曲线进行对比，采用极值后平均撞击力和吸能系数对高温作用下构件的抗侧向撞击性能进行量化分析，并分析了600 ℃下构件撞击全过程。结果表明:温度对钢管混凝土构件的侧向撞击性能影响明显，随着温度升高，构件跨中挠度大幅增加，撞击时程变长；高温下构件的撞击力时程曲线与常温下差异明显，高温下曲线可分为震荡阶段、下降阶段和卸载阶段；构件主要通过整体弯曲变形吸收落锤的动能，随着温度升高，极值后平均撞击力和吸能系数逐渐降低，表明构件的抗撞击性能逐渐降低，当温度超过400 ℃后，构件抗撞击性能损失严重。     

带模拟装药弹体高速冲击岩石靶时的断裂特性

针对高速侵彻过程中的弹体破碎断裂问题，本文中设计2种不同壁厚的试验弹，进行约1 000 m/s着速的高强度岩体侵彻试验，试验表明:在该高着速条件下，两种结构的试验弹体均发生完全破碎且未能有效侵入岩石靶，而岩石靶体仅在表层产生粉碎性破坏；另外，高速侵彻岩石靶的弹体头部破碎情况与侵彻金属薄靶有所区别。在试验基础上，利用Autodyn-3D建立了弹体侵彻岩石靶的物理模型，结合SPH算法与Mott失效模型对弹体破坏过程进行了数值模拟，可有效地揭示弹体破碎机理，并进一步讨论模拟装药和小范围内不同高速对弹体破坏的影响。试验结果和建立的数值模型可为研究高速侵彻中弹体结构安全提供参考。     

The dynamic response of clamped rectangular Y-frame and corrugated core sandwich plates

The dynamic response of fully clamped, monolithic and sandwich plates of equal areal mass has been measured by loading rectangular plates over a central patch with metal foam projectiles. All plates are made from AISI 304 stainless steel, and the sandwich topologies comprise two identical face-sheets and either Y-frame or corrugated cores. The resistance to shock loading is quantified by the permanent transverse deflection at mid-span of the plates as a function of projectile momentum. At low levels of projectile momentum both types of sandwich plate deflect less than monolithic plates of equal areal mass. However, at higher levels of projectile momentum, the sandwich plates tear while the monolithic plates remain intact. Three-dimensional finite element (FE) calculations adequately predict the measured responses, prior to the onset of tearing. These calculations also reveal that the accumulated plastic strains in the front face of the sandwich plates exceed those in the monolithic plates. These high plastic strains lead to failure of the front face sheets of the sandwich plates at lower values of projectile momentum than for the equivalent monolithic plates.     

The dynamic response of composite sandwich beams to transverse impact

The dynamic response of glass fibre–vinylester composite beams is measured by impacting the beams at mid-span with metal foam projectiles. The beams exist in composite monolithic form, and in sandwich configuration with composite face-sheets and a core made from PVC foam or end-grain balsa wood. High-speed photography is used to measure the transient transverse deflection of the beams and to record the dynamic modes of deformation and failure. For both monolithic and sandwich configurations, a flexural wave travels from the impact site towards the supports. Ultimate failure of the monolithic and sandwich beams is by tensile tearing of the faces. The sandwich beams also exhibit cracking of the core, and face-sheet delamination. The dynamic strength of the beams is quantified by the maximum transient transverse deflection at mid-span of the beams as a function of projectile momentum. It is demonstrated that sandwich beams can outperform monolithic beams of equal mass. The trade-off between core strength and core thickness is such that a low density PVC foam core outperforms a higher density PVC foam core. End-grain balsa wood has a superior stiffness and strength to that of PVC foam in compression and in shear. Consequently, sandwich beams with a balsa core outperform beams with a PVC foam core for projectiles of low momentum. The order reverses at high values of projectile momentum: the sandwich beams with a balsa wood core fail prematurely in longitudinal shear by splitting along the grain.     

An analytical model for high velocity impacts on thin CFRPs woven laminated plates

An analytical model to study the impact process of a spherical projectile penetrating at high velocity into a carbon/epoxy plain woven laminate is developed in this work. The model is based on an energy balance, where the kinetic energy of the projectile is absorbed by the laminate by three different mechanisms: laminate crushing, linear momentum transfer and tensile fiber failure. A non-homogeneous differential equation is obtained. A subsequent simplification using regular perturbation analysis gives a closed-form solution that allows the approximative calculation of the residual velocity and hence the ballistic limit. The model is validated with the results of experimental tests in which the residual velocity is measured by means of high speed cameras.     

球形弹丸正撞击薄板防护屏碎片云特性研究

以质量、动量和能量守恒为基础,结合平面激波理论和热力学理论,对球形铝弹丸正撞击薄板铝防护屏形成的碎片云特性进行了建模,模型计算结果与实验结果吻合较好。利用该模型对多种工况下碎片云特性进行了计算,结果表明:(1)碎片云质心速度和膨胀速度随撞击速度和弹丸直径的增加而增大,随防护屏厚度增大而减小;膨胀半角随撞击速度和防护屏厚度的增大而增大,随弹丸直径的增大而减小;(2)速度和膨胀半角变化曲线具有相似性;(3)对于给定的弹丸和防护屏材料,碎片云中受到激波加载部分的材料各相态质量分数只与弹丸撞击速度有关。这些规律与实验规律吻合。     

钨合金弹体对混凝土靶的超高速侵彻机理

为研究钨合金弹体超高速侵彻混凝土靶的相关机理,构建了适用于超高速撞击的金属强度模型、失效模型和混凝土的本构模型,对93钨合金弹体超高速撞击混凝土靶问题进行了数值模拟。开展了钨合金弹体超高速撞击混凝土靶实验,分析了靶板成坑特性,研究了侵彻总深度和残余弹体长度随撞击速度的变化规律,理论分析了长杆钨弹超高速撞击混凝土的侵彻模型和混凝土靶内的应力波传播。得到以下主要结论:（1）利用金属及混凝土的新本构模型获得的超高速撞击混凝土靶的破坏形貌数值模拟结果与实验结果一致;（2）超高速撞击条件下混凝土靶成坑为“弹坑＋弹洞”形,成坑体积与弹体动能近似成正比;（3）超高速撞击条件下,侵彻深度随弹速提高呈现先增大后减小的现象,高速段侵深降低是弹体经历销蚀侵彻后“刚体侵彻阶段”减少造成的;（4）建立的钨合金超高速撞击混凝土侵彻分析模型,可用来预估侵彻深度、残余弹长、蘑菇头直径等参数;（5）采用建立的超高速撞击混凝土靶内应力波传播理论模型得到的计算结果与实验结果吻合较好。     

Evolution of the micromechanical properties of impacted granular materials

The impact of a projectile on a granular material induces important changes in the micromechanical properties of the impacted material. These changes are studied using a Discrete Element Method model of the impact. The numerical results show that the impact first entails an energy propagation from the impact point to the limit of the sample through the existing force chains. A significant kinetic energy and a total breakage of the existing contact force chains are the main consequences of the energy propagation. During the long recovery balance phase observed after the energy propagation phase, frictional processes cause the kinetic energy dissipation. The motions of the particles and the numerous contact openings first prevent the formation of stable force chains. However, for long durations after the beginning of the impact, contact openings stop. The balance recovery phase therefore finally results in the creation of new stable contact force networks.     

修正金属本构模型在超高速撞击模拟中的应用

为更加准确地计算93钨合金弹超高速撞击Q345钢板问题,构建了修正的金属本构模型。引入GRAY三相物态方程描述材料相态变化,采用Johnson-Cook强度模型描述撞击后期材料的力学行为。结合封加波损伤演化模型以及Johnson-Cook失效模型描述不同应力三轴度下材料的拉伸、剪切失效行为;引入曹祥提出的断裂演化模型,描述材料失效后应力归零的过程。通过对比超高速撞击数值模拟结果与实验结果,验证了本构模型的适用性,并进一步分析了典型弹靶撞击条件下破片群的空间分布特征。研究结果表明:基于修正金属本构模型获得的超高速撞击靶板穿孔直径、弹体侵蚀长度、破片群扩展速度结果与实验结果一致;GRAY三相物态方程能够相对准确地给出弹体撞击首层靶板以及剩余弹体、破片群撞击第2层靶板时弹靶材料的熔化情况;封加波损伤演化模型能够准确判断超高速撞击过程中靶板是否产生层裂破坏;综合封加波损伤演化模型、Johnson-Cook失效模型以及曹祥提出的断裂演化模型后,数值模拟获得的破片群撞击后效靶板的穿孔面积与累积数量的统计曲线结果与实验结果一致;获得了典型条件下的柱形93钨弹体超高速撞击Q345靶板破片群空间分布结果,破片群的前端具有较高的质量、轴向动量以及横向动量（绝对值）。     

Penetration force measurement of thin plates by laser doppler anemometry

An experimental method is developed by applying laser-Doppler anemometry (LDA) to measure the velocity history of a projectile impact on a target. The force history was obtained by assuming the projectile to be rigid and by differentiating a polynomial which best fitted the obtained velocity data. This method was then applied to the impact of a hemispherically tipped striker on 1-mm thick aluminum plates. The initial velocity of the projectile ranged from several m/s to 93.5 m/s. When these velocities were low enough not to produce cracks on the plate, double-force peaks were observed in each test, and the value corresponding to the higher peak was proportional to the initial velocity of the projectile. Further, this peak force occurred when the projectile just started to rebound, and spring-back behavior of the plate always occurred. Once the cracks on the target started to form, the peak impact force dropped, and the double-peak phenomenon disappeared. It was also found that the peak force was at approximately the same level for impact speeds above the ballistic limit of the plate. On the other hand, a power law was found to perfectly fit the relationship between the absorbed kinetic energy and the initial velocity of the projectile when it was lower than the ballistic limit.     

弹体侵彻厚混凝土靶迎弹面成坑效应

为研究弹体侵彻厚混凝土靶的迎弹面成坑效应，总结了侵彻实验中的成坑现象，分析了经验公式对成坑深度、成坑直径和成坑角等成坑效应的预测效果；考虑了撞击速度、靶板强度、配筋以及弹体直径和质量等因素的影响，采用量纲分析方法建立了新型成坑效应计算公式及成坑阶段耗能计算公式；基于新型成坑效应计算公式，对成坑效应的影响因素和成坑耗能进行了参数化分析。结果表明:无量纲成坑深度受靶板强度、配筋率和弹体质量的影响较大；对于钢筋混凝土，成坑深度随撞击速度提升呈先增大后减小再增大的变化规律；在常见的侵彻速度和质量范围内，成坑角为15°～24°，质量对成坑角影响较小；迎弹面成坑耗能占弹体总动能的10%～25%，且配筋率和靶板强度对成坑耗能比例的影响较小；弹体质量越小，成坑阶段耗能占比越大。新型成坑效应计算公式对成坑深度、直径和角度的计算结果与实验数据吻合较好，可为侵彻弹体设计和工程防护提供参考。     

冰弹丸高速发射方法研究

冰弹丸的高速发射技术是认识冰弹丸高速撞击靶板损伤效应的瓶颈技术。利用二级轻气炮作为发射装置,采用液氮杜瓦制冷系统实现温度控制,选择铝质膜片对气压进行控制,利用高速摄像机对冰弹丸的完整性进行观测,并利用丝网测速仪对弹丸速度进行测量,实现了对冰弹丸的高速发射。通过实验验证了该方法的可行性,得到了完整冰弹丸在1516m/s下撞击铝合金靶板的实验结果。