可破碎颗粒体在动力载荷下的耗能特性

采用离散元的数值方法, 通过连接键将若干小颗粒绑定为一个具有不规则外形的大颗粒体, 设置不同连接键强度模拟了颗粒体在外加动力载荷下破碎过程, 并探讨其中系统能量耗散特性. 计算结果表明, 颗粒体的破碎程度决定了系统能量耗散率, 即内部耗能占外界输入能量的比例. 破碎率越高, 颗粒间相互摩擦和碰撞越剧烈,系统能量耗散率越高. 同时, 在循环载荷下系统内颗粒体破碎绝大部分发生在加载初期, 随着颗粒体的分解破碎速率逐渐减小, 系统耗能能力也随之降低.      

颗粒材料破碎演化路径细观热力学机制

颗粒材料在高应力环境下会发生颗粒破碎现象,颗粒破碎不仅影响颗粒材料的力学特性,同时与大量工程问题密切相关.目前的相关研究主要集中在唯象地描述颗粒破碎的演化以及破碎对力学特性的影响层面,对颗粒破碎演化路径的物理机制研究较少.本文基于热力学框架,采用细观力学中细观-宏观的均匀化方法推导了颗粒体系弹性能和破碎能量耗散,并在最大能量耗散的假设下,在热力学框架内,建立了理想化的无摩擦球体颗粒等向压缩过程的弹性-破碎模型,阐述了颗粒材料破碎演化路径细观热力学机制.由于模型的推导不依赖任何唯象的经验公式,因此模型中包含的参数均有明确的物理意义.模型预测与前人试验结果对比表明,材料的初始级配对弹性压缩模量和破碎应力的影响并不相同:不同分形维数级配对应的弹性体变模量存在极大值,而破碎应力却随着分形维数的增大单调递增;颗粒破碎的演化符合最大能量耗散原理,且颗粒材料的压缩曲线可以分为弹性-破碎-拟弹性3个机制不同的阶段.      

Effect of particle polydispersity on micromechanical properties and energy dissipation in granular mixtures

A series of numerical tests was conducted to study the micromechanical properties and energy dissipation in polydisperse assemblies of spherical particles subjected to uniaxial compression. In general, distributed particle size assemblies with standard deviations ranging from 0% to 80% of the particle mean diameter were examined. The microscale analyses included the trace of the fabric tensor, magnitude and orien- tation of the contact forces, trace of stress, number of contacts and degree of mobilization of friction in contacts between particles. In polydisperse samples, the average coordination numbers were lower than in monodisperse assemblies, and the mobilization of friction was higher than in monodisperse assemblies due to the non-uniform spatial rearrangement of spheres in the samples and the smaller displacements of the particles. The effect of particle size heterogeneity on both the energy density and energy dissipation in systems was also investigated.     

Effect of particle polydispersity on micromechanical properties and energy dissipation in granular mixtures

A series of numerical tests was conducted to study the micromechanical properties and energy dissipation in polydisperse assemblies of spherical particles subjected to uniaxial compression. In general, distributed particle size assemblies with standard deviations ranging from 0% to 80% of the particle mean diameter were examined. The microscale analyses included the trace of the fabric tensor, magnitude and orientation of the contact forces, trace of stress, number of contacts and degree of mobilization of friction in contacts between particles. In polydisperse samples, the average coordination numbers were lower than in monodisperse assemblies, and the mobilization of friction was higher than in monodisperse assemblies due to the non-uniform spatial rearrangement of spheres in the samples and the smaller displacements of the particles. The effect of particle size heterogeneity on both the energy density and energy dissipation in systems was also investigated.     

Numerical study of coarse coal particle breakage in pneumatic conveying

Pneumatic conveying of coarse coal particles with various pipeline configurations and swirling intensities was investigated using a coupled computational fluid dynamics and discrete element method. A particle cluster agglomerated by the parallel-bond method was modeled to analyze the breakage of coarse coal particles. The numerical parameters, simulation conditions, and simulation results were experimentally validated. On analyzing total energy variation in the agglomerate during the breakage process, the results showed that downward fluctuation of the total particle energy was correlated with particle and wall collisions, and particle breakage showed a positive correlation with the energy difference. The correlation between the total energy variation of a particle cluster and particle breakage was also analyzed. Particle integrity presented a fluctuating upward trend with pipe bend radius and increased with swirling number for most bend radii. The degree of particle breakage differed with pipeline bending direction and swirling intensity: in a horizontal bend, the bend radius and swirling intensity dominated the total energy variations; these effects were not observed in a vertical bend. The total energy of the particle cluster exiting a bend was generally positively correlated with the bend radius for all conditions and was independent of bending direction.     

基于离散单元法的颗粒阻尼耗能减振特性研究

采用离散单元法并从能量耗散的角度研究颗粒阻尼对系统减振特性的影响。建立了颗粒介质细观下的法向、切向和滚动方向的粘弹性接触模型和能量耗散模型,通过冲击激励和简谐激励下系统振动响应的多参数能量耗散分析来研究颗粒阻尼的耗能机理和减振特性。数值试验表明,颗粒介质可以在一个较宽的振动幅值范围内有效的发挥其阻尼效应,其耗能具有阶梯状周期性的特点。填充率是影响颗粒阻尼耗能减振效果的主要工程可控参数并对系统共振频率产生重大影响,当填充率接近极值时,系统出现无阻尼共振及共振频率超出无颗粒系统固有频率的现象。系统在最优填充率下共振时,颗粒与箱体保持恒定相位差的超振幅稳态运动。较小粒径的颗粒可以提高能量耗散率并使振动系统更快趋向静平衡状态,而恢复系数和摩擦系数则对法向和切向耗能的比值有较大影响。     

Statistics of particle interactions in dense granular material under uniaxial compression

Stress evolution in a dense granular material is closely related to interactions of contacting particles. We investigate statistics related to particle interactions and the relationship between the averaged local relative motion and the macroscopic motion. The validity of the Voigt and Reuss assumptions is examined, and extensions to these assumptions are proposed. Effects of history in the dense granular material are investigated. Statistical samples used in this paper are obtained using three-dimensional numerical simulations of dense granular media under uniaxial cyclical compression. The results show that stresses arise mostly from normal forces between particles, and direct contributions from frictional tangential forces between particles are small. Tangential friction, however, significantly increases the particle contact time, and thus reduces the rate of contact breakage. The contact breakage rate is demonstrated to be a stress relaxation rate. Therefore, stress increases significantly with friction between particles as a result of prolonged relaxation time.     

Breakage mechanics—Part II: Modelling granular materials

The compression of granular materials has been traditionally modelled with the limitations of classical elasto-plasticity. The energy was implicitly assumed to dissipate from the frictional interaction of particles. However, the fact that brittle granular materials crush suggests that energy must also be dissipated from the fracturing of the grains, as in fracture mechanics. The concept of breakage as a thermomechanical internal variable was introduced in Part I [Einav, I., 2006. Breakage mechanics—Part I: theory. J. Mech. Phys. Solids 00,000-000] to describe the fracturing mechanisms. The theory allows to treat ideal theoretical materials that undergo dissipation purely from breakage with no other mechanism allowed for the energy consumption. However, as accounted for in elasto-plasticity, dissipation must also occur from the frictional rearrangement of grains. The combination of the two dissipative mechanisms of breakage and plasticity must therefore be investigated, as we do in this paper. Those two mechanisms are generally coupled, in the sense that one inevitably appears when the other develops. Plastic dissipation emerges as a by-product of breakage dissipation because after grains crush, local rearrangement must occur. This scenario may be termed an ‘active breakage mechanism’, and typifies compression deformations. In shear the plastic dissipation is dominant but breakage appears inevitably from grains abrasion. This scenario may be called a ‘passive breakage mechanism’. Based on the coupling assumption, models are developed for granular materials. In particular, we show that in compression isotropic hardening of sands may appear without involving plastic strains, i.e., independent of frictional dissipation. This interpretation of hardening is different from the one used in classical critical state soil mechanics. However, frictional dissipation leads to plastic straining that are necessary for the models to be predictive in unloading.     

弹塑性微凸体侧向接触相互作用能耗

传统的结合面研究多基于光滑刚性平面与等效粗糙表面接触假设,忽略了结合面上微凸体侧向接触及相邻微凸体之间的相互作用,这导致理论模型与实际结合面存在较大出入.针对承受法向静、动态力的机械结合面,从微观上研究了微凸体侧向接触及相互作用的接触能耗.将法向静、动态力分解为法向分力和切向分力,获取弹性/弹塑性/塑性阶段考虑微凸体侧接触及相互作用的加、卸载法向分力-变形和切向分力-位移的关系.通过力的合成定理,从而获取加、卸载法向合力与总变形之间的关系,由于法向分力产生的塑性变形及切向分力产生的摩擦,导致加载、卸载法向合力-总变形曲线存在迟滞回线.通过对一个加、卸载周期内的法向合力-总变形曲线积分,获得一个周期的微凸体接触能耗,包括应变能耗及摩擦能耗.仿真分析表明:微凸体在3个阶段的能耗均随变形的增大而非线性增大.微凸体侧向接触角度越大,能耗越大,且在弹性阶段最为明显.在弹性阶段,仅存在侧向的摩擦能耗,故结合面在低载荷作用下必须采用双粗糙表面假设.在塑性阶段,由于微凸体接触能耗为应变能耗,且接触角对其能耗影响甚微,故结合面在大载荷作用下可采用单平面假设对其进行研究.相对于KE和Etsion模型,本文提出的模型与Bartier的实验结果更吻合.     

一种滞弹簧耗能的新型离散元滚动阻力模型研究

颗粒间滚动阻力对颗粒体系的稳定性起着重要作用. 在传统的离散元法中, 滚动阻力模型通常由转动弹簧、转动黏壶和摩擦元件表达, 颗粒滚动动能由黏滞力(矩)和摩擦力做功耗散. 由于黏滞力(矩)与滚动速度相关, 临近静止状态的颗粒滚动速度变小, 动能耗散减弱, 传统的离散元模拟得到颗粒由滚动到静止耗费的时间比试验观测的结果要长. 为解决这一问题, 基于摩擦学理论分析了滚动阻力产生的材料滞弹性机理, 将其引入离散元滚动阻力模型, 提出了一种速度无关型动能耗散的滞弹簧, 给出了滞弹簧的弹性恢复力计算公式, 建立了一种新型的离散元滞弹性滚动阻力模型(HDEM). 为验证新型滚动阻力模型的正确性, 通过一个光学物理试验对单个圆形颗粒试件的自由滚动过程进行了测量, 将测量数据与新型的滞弹型离散元模型和传统离散元模型计算结果进行了对比. 结果显示, 基于滞弹性滚动阻力模型HDEM计算结果与试验数据吻合程度更高, 而且模拟得到的颗粒摆动频率更符合试验现象.      

地震作用下摩擦耗能减震结构优化

提出了地震作用下摩擦耗能支撑参数优化的一种新的数学模型。它的特点是在给定的几条地震波作用下,在满足框架的规范层间位移角限值要求下,框架各层安装的耗能支撑刚度之和最小,从而实现安装较少的耗能装置而能达到相同的抗震要求。本文利用遗传算法求解该优化问题,并编制了基于遗传算法的几条地震波分别作用下的摩擦耗能框架结构支撑参数弹塑性动力优化分析程序GAOFF。最后给出了一6层框架在三条地震波分别输入下的优化算例,分析结果表明本文方法是有效的。     

考虑接触界面材料微观结构与势能参数的滑动摩擦计算研究

根据独立振子模型的能量耗散机理,提出了在无磨损界面摩擦中,利用通用界面粘附能量函数计算滑动时接触界面势能的变化从而计算摩擦力和摩擦系数的方法,建立摩擦力和摩擦系数与金属材料表面能与微观结构参数的关系,并利用已有的实验数据进行计算.结果表明,摩擦力与金属材料的表面能基本呈线性关系,与比例参数呈线性负相关关系,与晶格常数基本无关.计算结果与粘附理论公式以及超高真空原子力显微镜试验结果吻合较好.     

Pattern formation in particle systems driven by color field

The structural evolution of systems with two kinds of particles driven in opposite directions, i.e., driven by a color field, is investigated by molecular dynamics simulations. Gaussian thermostat,a common treatment to restrict the thermal velocity of the particles in the systems, has been used so as to account for the dissipation of heat and allow the system to reach a steady state. It has been found that with the increase of the strngth of driving force (F), the system undergoes an obvious structural transition from an initially random mixing state to a state characterized by separate lanes and in each lane only one kind of particles exists. The analysis shows that the reason for the formation of lane structure is not only the increase of F but also the variation of particle friction coefficient. While using Ganssian thermostat the particle friction coefficient becomes a function of F. Increasing F leads to high particle friction coefficient and inevitably results in lane formation for strong enough driving force. When lifting the effect of F on friction coefficient and choosing a constant friction coefficient,our results show that for a given F there always exists a critical value of friction coefficient higher than which the system will develop into lane structure.     

假肢接受腔-人体皮肤界面往复摩擦能量损耗的研究

利用万能材料试验机和多功能数据采集分析系统,在往复摩擦模式下模拟假肢接收腔-残肢皮肤穿戴界面,研究了不同切向力作用下两种常用的假肢接受腔材料与皮肤的摩擦能量损耗.结果表明:接受腔-皮肤界面往复摩擦过程中,随着切向力的增加,皮肤发生的弹性变形增大,皮肤摩擦前后的温差增加,摩擦过程中的能量损耗增大.接受腔材料表面结构直接影响着残疾人穿戴假肢的舒适性.相较于表面粗糙度较大的国内生产接受腔材料,表面粗糙度较小的进口接受腔材料降低了皮肤的摩擦能耗,减轻了皮肤的发热和变形,降低了皮肤损伤的风险,提高了与皮肤界面的适配性.     

Numerical simulation of fluid–particle hydrodynamics in a rectangular spouted vessel

A three-dimensional, Eulerian simulation was developed to describe isothermal, two-phase flow of the continuous (water) and dispersed (solid particles) phases in a rectangular spouted vessel. The mass and momentum conservation equations for each phase were solved using the finite volume technique, which treats each phase separately, while coupling them through drag, turbulence, and energy dissipation due to particle fluctuations. Particle–particle interactions via friction were also included.     

Energy dissipation mechanism and experiment of particle dampers for gear transmission under centrifugal loads

As a passive means of vibration reduction, particle damping is mainly applied to the horizontal or vertical steady field. However, it is seldom applied to centrifugal fields. Under high speed and heavy loading, the vibration of tooth surfaces of gear transmissions becomes more severe shortening gear service life and augmenting noise. Under centrifugal loading, the particle system exhibits different characteristics, for example, particles are extruded at the end farthest from the center. We investigated gears with drilled via holes filled with damping particles. Using the discrete-element method, we developed an energy dissipation model for the particle system accounting for friction and inelastic collisions. Energy dissipation and damping characteristics of this system were analyzed. Experiments were also conducted with the gear system having different particle filling rates. The results show that this filling rate is an important parameter associated with particle damping in a centrifugal field. An unsuitable filling rate would significantly reduce damping effectiveness. With changes in rotation speed and load, the gear transmission system has different optimal filling rates. The results provide guidelines for the application of particle damping in centrifugal fields of gear transmissions.     

Impact comminution of solids due to local kinetic energy of high shear strain rate: I. Continuum theory and turbulence analogy

The modeling of high velocity impact into brittle or quasibrittle solids is hampered by the unavailability of a constitutive model capturing the effects of material comminution into very fine particles. The present objective is to develop such a model, usable in finite element programs. The comminution at very high strain rates can dissipate a large portion of the kinetic energy of an impacting missile. The spatial derivative of the energy dissipated by comminution gives a force resisting the penetration, which is superposed on the nodal forces obtained from the static constitutive model in a finite element program. The present theory is inspired partly by Grady's model for expansive comminution due to explosion inside a hollow sphere, and partly by analogy with turbulence. In high velocity turbulent flow, the energy dissipation rate gets enhanced by the formation of micro-vortices (eddies) which dissipate energy by viscous shear stress. Similarly, here it is assumed that the energy dissipation at fast deformation of a confined solid gets enhanced by the release of kinetic energy of the motion associated with a high-rate shear strain of forming particles. For simplicity, the shape of these particles in the plane of maximum shear rate is considered to be regular hexagons. The particle sizes are assumed to be distributed according to the Schuhmann power law. The condition that the rate of release of the local kinetic energy must be equal to the interface fracture energy yields a relation between the particle size, the shear strain rate, the fracture energy and the mass density. As one experimental justification, the present theory agrees with Grady's empirical observation that, in impact events, the average particle size is proportional to the (−2/3) power of the shear strain rate. The main characteristic of the comminution process is a dimensionless number B_a (Eq. (37)) representing the ratio of the local kinetic energy of shear strain rate to the maximum possible strain energy that can be stored in the same volume of material. It is shown that the kinetic energy release is proportional to the (2/3)-power of the shear strain rate, and that the dynamic comminution creates an apparent material viscosity inversely proportional to the (1/3)-power of that rate. After comminution, the interface fracture energy takes the role of interface friction, and it is pointed out that if the friction depends on the slip rate the aforementioned exponents would change. The effect of dynamic comminution can simply be taken into account by introducing the apparent viscosity into the material constitutive model, which is what is implemented in the paper that follows.     

基于颗粒流的混凝土材料数值实验研究

混凝土材料是一种典型的非均质材料,其力学行为和破坏过程很复杂.本文利用PFC颗粒流软件的内嵌fish语言编写了骨料生成程序、边界条件控制程序,程序模拟了MTS伺服功能,实现了虚拟的混凝土单轴压缩试验.相对于一般有限元程序投放的骨料,文中骨料可以破裂.在生成的数值试件基础上,设计了单级配和全级配混凝土单轴压缩试验,并对其破裂形态、裂纹扩展过程、应力应变曲线和破裂过程能量变化规律进行了研究.利用该虚拟试验平台,可对混凝土的动、静力学特性进行系统研究.     

A slip-line plasticity analysis of sliding friction of rough surfaces exhibiting self-affine (fractal) behavior

A plasticity analysis of sliding friction of rough (fractal) surfaces sliding against smooth surfaces was developed based on a slip-line model of a rigid spherical asperity (wear particle) plowing and cutting through a soft semi-infinite medium. Solutions of the fraction of fully plastic asperity microcontacts responsible for the evolution of friction and energy dissipation were obtained in terms of the total normal load (global interference), interfacial adhesion characteristics, topography (fractal) parameters of the hard surface, and elastic–plastic material properties of the soft surface. This was accomplished by incorporating the slip-line model of a single microcontact into a friction analysis of sliding surfaces demonstrating multi-scale roughness. Numerical results provide insight into the effects of global interference (normal load), fractal parameters (surface roughness) of the hard surface, interfacial shear strength (adhesion), and material properties of the soft surface on plastic deformation at the microcontact level, global coefficient of friction, and frictional energy dissipated during sliding.     

皮肤往复摩擦条件下的损伤研究

选用成年健康活体兔皮肤替代残肢皮肤,在UMT-Ⅱ多功能摩擦磨损试验机上通过模拟残肢皮肤/接收腔界面间的黏着和中间过渡摩擦行为,对兔皮肤的摩擦系数、Ft-D曲线、能量损耗以及这两种摩擦行为所引起的兔皮肤摩擦损伤机理进行了研究.结果表明:在0.5或3.5 N法向力作用下,随着往复摩擦位移从±2.5 mm增加到±5 mm,兔皮肤的摩擦行为从黏着状态转变到中间过渡状态,往复摩擦过程中的能量损耗增加;由于中间过渡状态下的兔皮肤摩擦系数较高,能量损耗较大,因而造成的皮肤损伤程度比黏着状态严重;与人体手掌或脚掌皮肤相比,兔皮肤的角质层较薄,很难形成摩擦水泡.