PROPAGATION OF IMPACT PULSE IN A CHAIN OF ELASTIC BALLS

采用Hertz 接触理论处理牛顿摆中球与球之间的相互接触,建立球链碰撞模型,观察到碰撞过程中类似孤立波的传播现象。为了便于处理“n 对m”碰撞的计算,编写了带有GUI (graphic user interface) 界面的计算程序,可以研究不同撞击球个数对撞击过程的影响,包括球心位置、球心速度以及相邻小球间的变形量与相互作用力随时间的变化规律。进一步比较了一维球链模型和三维有限元模型的结果,并实验观察了碰撞过程的细节,得到一致性结论。     

Ball spin generation for oblique impacts with a tennis racket

In this paper, we describe an experimental investigation of the oblique impact between a tennis ball and head clamped tennis racket. It was found that the magnitude of the ball rebound spin was not a function of the material, gage or tension of the string used in the tennis racket. Furthermore, it was concluded that all strings exhibit a sufficiently large friction coefficient that the ball begins to roll during impact. There is anecdotal evidence from tennis players that suggests that a high string tension or a rough string surface enable them to impart more spin to the ball. For example, players have been quoted as saying that a high string tension makes the strings “bite” into the ball, giving more spin. The data reported in this study do not support these observations. Analysis of the experimental data has shown that the balls are rebounding from the surface with more spin than would typically be associated with rolling. A second experiment showed that the balls commenced rolling at the mid-point of the impact. This information was used in a theoretical model to show that the spin that acts on the ball during the impact can be higher than the value of the rolling spin at the end of the impact.     

Numerical study of particle motion near a charged collector

The behavior of particles impacting the surface of a charged droplet involves adhesion, rebound, and submersion. In the present study, a numerical model for simulating particle impacts on charged droplets is presented that takes into account the various impact modes. With the droplet considered as a solid boundary, the criterion for rebounding is that the particle’s impact angle is <85°. The simulated trajectories of the particles are verified by comparing with experimental data for low-velocity particles to assess the reliability of the model. For impact angles >85°, particles undergo three distinct modes depending on normal impact velocities. The critical velocity of adhesion/rebound and rebound/submersion is used to identify the mode that the particles are undergoing. The criteria are also verified by comparing with analytical data. The results show that the impact angle of particles increases with increasing Coulomb number and decreases dramatically with increasing Stokes number, both of which lead to a high probability for particle rebound.     

优化形式的刚体系统动力学模拟

基于广义坐标形式的高斯最小拘束原理来研究刚体系统的动力学问题的优化方法. 相比目前动力学建模普遍采用的质点形式的高斯最小拘束原理,广义坐标形式的高斯最小拘束原理因对所选择的广义坐标没有要求,而使得建模过程更简单及具有更强的通用性. 本文分别建立了有约束和无约束条件下问题的优化动力学模型,对问题进行了动力学数值模拟,并与用拉格朗日微分方程处理的模型进行了对比分析,从而验证了所提方法的有效性.     

勒让德伪谱法求解三维刚体摆姿态运动最优控制问题

研究伪谱法求解三维刚体摆姿态运动最优控制问题. 针对三维刚体摆这类含有约束的力学模型,提出了基于勒让德伪谱法的三维刚体摆姿态最优控制方法. 利用插值逼近设计了三维刚体摆姿态运动最优控制算法,得到了三维刚体摆的姿态最优控制轨迹,并结合松弛参数来控制插值点的取值,寻找满足的可行解. 仿真结果表明,基于勒让德伪谱法的最优控制算法使得三维刚体摆能以较小的误差运动到期望的末端姿态,且计算速度快,能够获得精度较高的控制输入量.     

基于修正型罗德里格斯参数的三维刚体摆姿态控制

3D 刚体摆是研究地球静止轨道航天器的一个力学简化模型, 它绕一个固定、无摩擦的支点旋转, 具有3 个转动自由度. 文章给出基于修正型罗德里格斯(Rodrigues) 参数描述的3D 刚体摆的姿态动力学方程, 针对3D 刚体摆姿态和角速度稳定的非线性控制设计问题, 基于无源性控制理论利用能量法设计了3D 刚体摆的系统控制器, 并证明了系统满足无源性. 构造了系统的李雅普诺夫(Lyapunov) 函数, 利用能量法设计出3D 刚体摆的姿态控制律, 并由拉萨尔(LaSalle) 不变集原理证明了该控制律的渐近稳定性. 仿真实验给出了3D 刚体摆在倒立平衡位置的姿态和角速度的渐近稳定性, 仿真实验结果表明基于能量方法的3D 刚体摆姿态控制是有效的.     

含摩擦与碰撞平面多刚体系统动力学线性互补算法

基于接触力学理论和线性互补问题的算法, 给出了一种含接触、碰撞以及库伦干摩擦, 同时具有理想定常约束(铰链约束) 和非定常约束(驱动约束) 的平面多刚体系统动力学的建模与数值计算方法. 将系统中的每个物体视为刚体, 但考虑物体接触点的局部变形, 将物体间的法向接触力表示成嵌入量与嵌入速度的非线性函数,其切向摩擦力采用库伦干摩擦模型. 利用摩擦余量和接触点的切向加速度等概念, 给出了摩擦定律的互补关系式; 并利用事件驱动法, 将接触点的黏滞-滑移状态切换的判断及黏滞状态下摩擦力的计算问题转化成线性互补问题的求解. 利用第一类拉格朗日方程和鲍姆加藤约束稳定化方法建立了系统的动力学方程, 由此可降低约束的漂移, 并可求解该系统的运动、法向接触力和切向摩擦力, 还可以求解理想铰链约束力和驱动约束力. 最后以一个类似夯机的平面多刚体系统为例, 分析了其动力学特性, 并说明了相关算法的有效性.      

蜻蜓滑翔时柔性褶皱前翅气动特性分析

蜻蜓是自然界优秀的飞行家,滑翔是其常见且有效的飞行模式.蜻蜓优异的飞行能力来源于其翅膀的巧妙结构,褶皱是蜻蜓翅膀上最为显著的结构之一,不仅提高了翅膀的刚度,还改变了其气动特性,而飞行过程中柔性翅膀会产生变形是蜻蜓翅膀的另一特性.为揭示蜻蜓在滑翔时,柔性褶皱前翅的变形,探究褶皱和柔性的共同作用对其气动特性的影响,基于逆向工程,依据前人的测量数据和研究成果,通过三维建模软件建立了蜻蜓三维褶皱前翅的计算流体力学(computational fluiddynamics,CFD)模型和计算结构力学(computational structuralmechanics,CSD)模型,并通过模态分析验证了此模型有足够的精度.基于CFD方法和CFD/CSD双向流固耦合计算方法分别对蜻蜓滑翔飞行时刚性和柔性褶皱前翅的气动特性进行了数值模拟,结果表明,柔性褶皱前翅受气动载荷后,翅脉和翅膜产生形变,柔性前翅上下表面压力差相较于刚性前翅减小了,从而其升力和阻力也减小了,而在大攻角时,变形后的前缘脉诱导出比刚性前翅更强的前缘涡.因此在攻角小于10$^\circ$时刚性前翅的气动特性优于柔性前翅,继续增大攻角,柔性前翅的气动特性则优于刚性前翅.前翅受载后气动响应时间短,翅尖的变形最大,仅仅产生了垂直于翅膀所在平面方向上的变形,而没有发生扭转,翼根处受到应力最大,褶皱上凸部分承受蜻蜓滑翔时前翅的主要载荷.      

链轮驱动系统的多刚体建模及分析

采用多刚体方法对链条链轮系统的驱动过程进行了建模和动力学分析．链节看作刚体，链轮看作等边刚性多边形，辊子看作两个链节之间的旋转副．啮合约束当作可变约束处理，当啮合力为正时，链节和链轮之间是固定约束，当力变负时，链节脱离链轮．利用龙格库塔法求解系统的拉格朗日方程．最后计算并分析了不同工况下的驱动过程，解释了系统的锁死现象，并进一步分析了两种避免锁死的处理方法．     

STUDY ON ANALYSIS METHOD OF “THREE BOX” IN HEAT AND MASS TRANSFER PROCESS OF RESERVOIR MEDIA

油藏多孔介质孔隙组成及结构变化多样,一些特性参数很难全部获得,精确描述和分析困难;另外,多孔介质内渗流过程水力条件和作用机理复杂,存在热流固耦合作用,目前的一些分析方法和研究模型具有一定的局限性.提出了油藏多孔介质的表征单元体（representative elementary volume,REV）描述表征方法;基于表征单元体建立了多孔介质的黑箱模型、灰箱模型和白箱模型,据此提出了多孔介质的“黑箱→灰箱→白箱”分析过程.基于黑箱模型和灰箱模型推导了REV导热系数计算公式、给出了REV热质传递过程的热平衡方程.结合中国油藏热采情况,对多孔介质导热系数变化规律和蒸汽驱热质传递特性进行了分析,得到了一些有意义的结果.该工作为多孔介质热质传递过程分析提供了新思路和新方法.     

Some computational aspects for analysis of low- and high-velocity impact of deformable bodies

In this paper, we present some numerical algorithms devoted to analysis of low- and high-velocity impact of deformable bodies. For modeling of contact phenomena involved in impact problems, an implicit method is used rather than the most popular penalty explicit method. For the appropriate time integration of the discretized equation of motion, a first order algorithm is suggested rather than the most common second order algorithms such as Newmark, Wilson, etc. Numerical results are presented for dynamic analysis such as low-velocity impact of an elastic plate on a rigid surface with rebounding and high-velocity impact of an elastic-plastic copper bar on a rigid wall.     

Asymptotic modeling of the impact of a spherical indenter on an elastic half-space

The impact of a rigid sphere on a homogeneous, isotropic elastic half-space in the absence of friction and adhesion is considered. The influence of the superseismic stage immediately following the moment of first contact upon the impact process is investigated in the frame of the Hertzian impact theory. The first order asymptotic approximation for the contact force in a three-dimensional dynamic contact problem with the slowly moving contact zone boundary is obtained and the corresponding asymptotic model of impact is developed. The motion of the indenter as it indents and rebounds from the elastic medium is analytically described. Explicit formulas are derived for the peak indentation depth, contact time, and rebound velocity as functions of the initial impact velocity, indenter mass, and characteristics of the elastic half-space.     

磁驱动准等熵压缩下LY12铝的强度测量

高压高应变率加载下材料的强度研究一直是冲击动力学的一个难题,目前动态载荷下材料的高压强度测量主要是基于平板撞击技术,冲击温升和应变率效应对材料强度的影响难以分离. 基于小型磁驱动加载装置CQ-4,开展了磁驱动准等熵压缩下LY12 铝的声速和强度测量的实验研究,讨论了考虑加载-卸载过程时磁驱动压缩实验的负载电极设计、实验样品设计、数据处理与分析等内容,并获得了12 GPa 压力范围沿加载-卸载路径的声速变化和峰值压力点的强度数据.      

多孔泡沫牺牲层的动态压溃及缓冲吸能机理研究

本文对强动载荷下多孔泡沫牺牲层的动态压溃行为及缓冲吸能机理进行了研究. 基于刚性-理想塑性-锁定(R-PP-L)及刚性-塑性硬化(R-PH)两类多孔泡沫材料本构, 建立了强动载荷下多孔泡沫牺牲层动态响应的理论分析模型, 分析了一维冲击波在多孔泡沫牺牲层中的传播规律; 利用Voronoi方法建立了多孔泡沫牺牲层的二维细观有限元模型, 获得了冲击载荷下多孔泡沫牺牲层的变形模式和动态响应曲线, 讨论了多孔泡沫材料的层间界面效应对多孔泡沫牺牲层缓冲吸能的影响. 研究结果表明, 考虑多孔泡沫材料塑性硬化影响的理论分析模型能够预测入射波在远端的反射及对多孔泡沫牺牲层的二次压缩过程和端部应力增强现象; 相比较存在界面的多孔泡沫牺牲层, 连续设计的多孔泡沫牺牲层可增强其缓冲吸能能力, 但在界面处增加设计刚性面板则能够降低界面胞元不完整对缓冲吸能的影响; 相同冲量载荷下, 端部应力峰值随冲击能量增大而增大, 而端部冲击波的反射可能是端部应力增强的主要诱因.      

Analysis of a rigid body obliquely impacting granular matter

The model of the oblique rigid body impact with a granular matter is studied. The force acting on the body is a linear superposition of a static (velocity-independent) friction force and a dynamic (velocity-dependent) resistance force. The impact of a sphere, a mathematical and a compound pendulum are modeled and simulated using different initial impact velocity conditions and different impact angles. We analyze how rapidly the rigid body impacting a granular media slows upon collision. For most of the analyzed cases the rigid body under high-force impact (higher initial velocity) comes to rest faster in a granular matter than the same body under low-force impacts (lower initial velocity). Researchers were able to explain this interesting phenomena, not shared by solids or liquids, for the vertical impact of spheres. The simulations for some configurations with small initial impact angles show that as the speed at which the rigid body impacts the media increases, the later it will come to rest.

     

ANTI-BLAST ANALYSIS OF CELLULAR SACRIFICIAL CLADDING

运用一维冲击波模型和三维细观有限元模型分析了多胞牺牲层的抗爆炸行为. 基于刚性-塑性硬化(R-PH)的多胞材料模型,建立了一维冲击波模型,得到了多胞牺牲层中冲击波传播的控制方程. 揭示了冲击波在多胞牺牲层中的传播特性,并阐述了附加质量和爆炸载荷强度两个参数对牺牲层设计的重要影响. 比较了基于刚性-理想塑性-锁定(R-PP-L) 模型和基于刚性-塑性硬化(R-PH) 模型的多胞牺牲层的结构设计,指出了两种模型的适用范围. 通过基于三维Voronoi 技术的细观有限元方法验证了基于R-PH 模型的多胞牺牲层结构的设计准则.     

页岩气储层暂堵转向压裂直井段暂堵球运移特性研究

考虑附加质量力、相间阻力、Basset力、重力、浮力等力学条件,建立了暂堵球运移过程动量守恒模型,通过四阶龙格-库塔方法对其求解,并以长宁区块气井直井段为例进行研究。结果表明:在直井段,暂堵球呈现先加速后减速,最后匀速运动的规律,暂堵球加速度ab≤0阶段,暂堵球运移速度达到最大峰值;随着暂堵球直径、密度及排量增大,暂堵球运移速度最大峰值均呈现增大趋势,暂堵球密度从1.13g/cm3增大到1.38g/cm3,暂堵球运移速度最大峰值从15.057m/s增加到17.19m/s;暂堵球雷诺数与暂堵球运移速度具有较好伴随性,暂堵球加速度ab>0的阶段,随着暂堵球直径增大,暂堵球雷诺数增大,暂堵球直径7mm同直径13mm相比,雷诺数峰值从214.18增大到397.76;随着排量的增大,暂堵球运移速度呈现增大的趋势,当排量为10m3/min,暂堵球到达直井段末端2536m时,运移速度为14.42m/s,运移总时间为2.835min。     

Simultaneous Impact of a Two-Link Chain

The study of the planar collisions of a rigid two-link chain with an external surface is presented. There are two points of contact during the impact process with friction, one for each link. For the normal impact, a nonlinear collision force is employed for the elasto-plastic deformation, and the Coulomb’s law models the frictional effects. The impact of the chain is studied numerically, and the results are compared with experiment data obtained from a high-speed camera. Special attention is given to simultaneous impact, slippage, and rebounds at the contact points.     

基于深度学习和近场动力学的层合板冲击工况识别

在冲击荷载作用下复合材料会产生断裂和分层等损伤。基于损伤数据对冲击工况进行识别，对改善复合材料的设计和确保其安全使用具有重要意义。基于此，本文提出一种基于深度学习和近场动力学（PD）理论的层合板冲击工况识别方法。首先使用改进的表面修正系数PD理论建立复合材料层合板刚体冲击损伤演化分析PD模型，PD模型数值模拟结果结合噪声数据增强技术构建层合板的冲击工况数据库；基于深度学习-卷积神经网络（CNN），对不同工况下的冲击损伤演化数据进行训练，实现对未知冲击工况的识别。结果显示，对于钢球冲击速度和角度的识别准确率均高于90%。     

Physical mechanics investigations of moving contact lines

移动接触线,指两种互不相溶的流体在固体表面形成移动的三相接触区域.移动接触区域跨越多个尺度,其中三相物质之间的相互作用影响着整个流场的动力学特征.由于在能源、航天、生物等领域中的重要应用和迅速发展,移动接触线在新的应用背景下发展了新的难题.标度分析是度量接触线自相似扩展的重要手段.本文以移动接触线的标度关系为主线,介绍了"力-电-热-化学"多场耦合环境下,亲水内角、微柱阵列、可溶解固体、水力压裂滞后区等复杂几何结构的刚性/柔性固体表面,采用物理力学方法对于移动接触线动力学属性研究的进展.通过跨尺度实验研究、大规模分子动力学模拟和分子动理论/水动力学理论相结合的方法,发现了类固体前驱膜、单分子前驱水链、锯齿形接触线等新现象.从原子尺度的界面结构到连续尺度的流动特性,讨论了移动接触线自相似扩展的标度关系,以及其驱动来源、能量耗散、边界条件等物理机制和规律,为多物理场中的"Huh-Scriven佯谬"探索了解答,为移动接触线的前景和应用提出了展望.