基于超二次曲面的非球形离散单元模型研究

通过对颗粒体系接触过程的运动学和动力学分析,建立了一种基于超二次曲面的非球形离散单元模型,该模型避免了球形接触模型描述颗粒形状的局限性,使离散单元法更接近物理事实,并在此基础上提出了计算求解模型的数值方法,实现了对复杂形状的颗粒体系的模拟计算。将所建立的数值计算方法进行了编程实现,并对模型和算法进行了算例测试,证实了本文所建立的非球形离散单元模型的可行性和正确性。测试结果表明,本文的模型能够比较准确地模拟复杂颗粒体系的真实运动,可为复杂颗粒体系的模拟研究提供一种新的数值计算方法。     

Calibration of granular material parameters for DEM modelling and numerical verification by blade-granular material interaction

The discrete element method (DEM) is a promising approach to model blade-granular material interactions. The accuracy of DEM models depends on the model parameters. In this study, a calibration process was developed to determine the parameter values. The particle size was the same as the real material and the particle shape was modelled using two spherical particles rigidly clumped together to form a single grain. Laboratory shear tests and compressions tests were used to determine the material internal friction angle and stiffness, respectively. These tests were replicated numerically using DEM models with different sets of particle friction coefficients and particle stiffness values. The shear test results are found to be dependent on both the particle friction coefficient and the particle stiffness. The compression test results show that it is only dependent on the particle stiffness. The combination of shear test and compression test results can be used to determine a unique set of particle friction and particle stiffness values. The calibration process was validated experimentally and numerically by modelling a blade moving through granular material. Results show that the forces acting on the blade can be accurately modelled with DEM and the maximum error is found to be 26%. The relative particle-blade displacements were used to predict the position and shape of the shear lines in front of the blade. A good qualitative correlation was achieved between the experiments and the DEM simulations.     

Characterization of artificial spherical particles for DEM validation studies

This paper describes a study in which advanced particle-scale characterization was carried out on spherical particles that can be used in experimental tests to validate discrete element method (DEM) simulations. Two types of particle, alkaline and borosilicate glass beads, made from two different materials, were considered. The particle shape, stiffness, contact friction properties and surface roughness were measured. The influences of hardness and roughness on the mechanical response of the particles were carefully considered. Compared to the alkaline beads, the borosilicate beads were more spherical and more uniform in size, and they exhibited mechanical characteristics closer to natural quartz sand. While only two material types were studied, the work has the broader implication as a methodology for selecting particles suitable for use in DEM studies and the key parameters that should be considered in the selection process are highlighted.     

Micromechanical study of elastic moduli of loose granular materials

In micromechanics of the elastic behaviour of granular materials, the macro-scale continuum elastic moduli are expressed in terms of micro-scale parameters, such as coordination number (the average number of contacts per particle) and interparticle contact stiffnesses in normal and tangential directions. It is well-known that mean-field theory gives inaccurate micromechanical predictions of the elastic moduli, especially for loose systems with low coordination number. Improved predictions of the moduli are obtained here for loose two-dimensional, isotropic assemblies. This is achieved by determining approximate displacement and rotation fields from the force and moment equilibrium conditions for small sub-assemblies of various sizes. It is assumed that the outer particles of these sub-assemblies move according to the mean field. From the particle displacement and rotation fields thus obtained, approximate elastic moduli are determined. The resulting predictions are compared with the true moduli, as determined from the discrete element method simulations for low coordination numbers and for various values of the tangential stiffness (at fixed value of the normal stiffness). Using this approach, accurate predictions of the moduli are obtained, especially when larger sub-assemblies are considered. As a step towards an analytical formulation of the present approach, it is investigated whether it is possible to replace the local contact stiffness matrices by a suitable average stiffness matrix. It is found that this generally leads to a deterioration of the accuracy of the predictions. Many micromechanical studies predict that the macroscopic bulk modulus is hardly influenced by the value of the tangential stiffness. It is shown here from the discrete element method simulations of hydrostatic compression that for loose systems, the bulk modulus strongly depends on the stiffness ratio for small stiffness ratios.     

The elastic response of a cohesive aggregate—a discrete element model with coupled particle interaction

A model is presented for the deformation of a cohesive aggregate of elastic particles that incorporates two important effects of large-sized inter-particle junctions. A finite element model is used to derive a particle response rule, for both normal and tangential relative deformations between pairs of particles. This model agrees with the Hertzian contact theory for small junctions, and is valid for junctions as large as half the nominal particle size. Further, the aggregate model uses elastic superposition to account for the coupled force–displacement response due to the simultaneous displacement of all of the neighbors of each particle in the aggregate. A particle stiffness matrix is developed, relating the forces at each junction to the three displacement degrees of freedom at all of the neighboring-particle junctions. The particle response satisfies force and moment equilibrium, so that the model is properly posed to allow for rigid rotation of the particle without introducing rotational degrees of freedom. A computer-simulated sintering algorithm is used to generate a random particle packing, and the stiffness matrix is derived for each particle. The effective elastic response is then estimated using a mean field or affine displacement calculation, and is also found exactly by a discrete element model, solving for the equilibrium response of the aggregate to uniform-strain boundary conditions. Both the estimate and the exact solution compare favorably with experimental data for the bulk modulus of sintered alumina, whereas Hertzian contact-based models underestimate the modulus significantly. Poisson's ratio is, however, accurately determined only by the full equilibrium discrete element solution, and shown to depend significantly on whether or not rigid particle rotation is permitted in the model. Moreover, this discrete element model is sufficiently robust, so it can be applied to problems involving non-homogeneous deformations in such cohesive aggregates.     

自适应无网格热弹塑性接触模型研究

提出一种自适应无网格热弹塑性接触求解模型,求解接触问题的线性规划-增量初应力法与基于应变能梯度的自适应无网格法相结合,给出了模型计算理论和算法实现.通过圆柱体与弹塑性平面热弹塑性接触算例对模型进行验证.对是否考虑材料应变硬化,是否考虑摩擦力和热输入,是否考虑材料屈服强度温度相关等情况的两种算例进行了讨论.结果表明,该模型能有效地求解考虑不同情况下的热弹塑性接触问题,在较真实地模拟接触状况的同时,具有较高的计算精度和计算效率.     

单峰接触研究及其在分形表面接触中的应用

基于有限元方法,建立了弹塑性单峰的接触模型.粗糙峰为理想的弹塑性材料,为了考虑不同的材料特性对微凸体变形的影响,分别对9种不同的材料进行了分析.根据有限元计算结果,分析了接触面积,平均接触压力和接触力与变形干涉量之间的关系,并进行了经验公式的拟合.单峰接触所经历的4个不同的阶段,以及不同阶段之间的转化点均作了明确的表达.然后,根据分形理论,将单峰接触模型扩展到了三维的粗糙表面的接触,并提出了一个计算接触表面法向刚度的模型.通过与实验数据和以往模型的结果对比,证明本文中所提出的模型具有较高的精度.     

具有随机参数的含裂纹板弯曲应力强度因子的统计分析

本文首次应用随机有限元法研究了具有随机参数的含裂纹板裂纹尖端弯曲应力强度因子的统计性质。文中首先给出了杂交模式的裂纹尖端奇异单元的刚度矩阵,然后基于随机场的局部平均理论和一阶泰勒展开得到了应力强度因子均值和方差的计算公式。作为数例,详细讨论了杨氏模量、泊松比及板厚度的不确定性对应力强度因子的影响。     

Dynamic coupling of fluid and structural mechanics for simulating particle motion and interaction in high speed compressible gas particle laden flow

In this work, structural finite element analyses of particles moving and interacting within high speed compressible flow are directly coupled to computational fluid dynamics and heat transfer analyses to provide more detailed and improved simulations of particle laden flow under these operating conditions. For a given solid material model, stresses and displacements throughout the solid body are determined with the particle–particle contact following an element to element local spring force model and local fluid induced forces directly calculated from the finite volume flow solution. Plasticity and particle deformation common in such a flow regime can be incorporated in a more rigorous manner than typical discrete element models where structural conditions are not directly modeled. Using the developed techniques, simulations of normal collisions between two 1 mm radius particles with initial particle velocities of 50–150 m/s are conducted with different levels of pressure driven gas flow moving normal to the initial particle motion for elastic and elastic–plastic with strain hardening based solid material models. In this manner, the relationships between the collision velocity, the material behavior models, and the fluid flow and the particle motion and deformation can be investigated. The elastic–plastic material behavior results in post collision velocities 16–50% of their pre-collision values while the elastic-based particle collisions nearly regained their initial velocity upon rebound. The elastic–plastic material models produce contact forces less than half of those for elastic collisions, longer contact times, and greater particle deformation. Fluid flow forces affect the particle motion even at high collision speeds regardless of the solid material behavior model. With the elastic models, the collision force varied little with the strength of the gas flow driver. For the elastic–plastic models, the larger particle deformation and the resulting increasingly asymmetric loading lead to growing differences in the collision force magnitudes and directions as the gas flow strength increased. The coupled finite volume flow and finite element structural analyses provide a capability to capture the interdependencies between the interaction of the particles, the particle deformation, the fluid flow and the particle motion.     

复合板旋转约束刚度推导及箱型结构屈曲分析

基于结构稳定性理论,推导出正交各向异性纤维树脂增强复合材料箱型结构各板连接处旋转约束刚度的近似表达式．与以往的经验公式相比,考虑了材料的正交各向异性和截面尺寸的影响．通过对箱型结构的正交各向异性比、截面属性以及纵横比对旋转约束刚度的影响进行参数研究,验证了旋转约束刚度新公式的精确性和适用范围．通过比较采用不同近似表达式得到的结果与有限元计算结果的对比,表明本文的旋转约束刚度新公式可以更精确地用来计算箱型结构的临界屈曲载荷．     

包含几何误差的机械结合面法向刚度研究

对包含几何误差的机械结合面进行离散化,离散后的微表面的基准平面高度满足结合面几何误差分布.每个微表面内,微凸体的高度只受粗糙度的影响.基于接触理论建立了微表面的法向刚度模型,通过对微表面模型集成获得了结合面的法向刚度模型.通过对所建模型的数值仿真,揭示了结合面法向刚度与间隙的非线性关系,几何误差的幅值和波长对法向刚度的影响以及非线性刚度对结合面振动特性的影响.计算结果表明:法向刚度随着间隙的减少而迅速增加,几何误差会导致结合面宏观上的局部接触和应力集中;在相同干涉量下,法向刚度随着几何误差幅值的增加而增加,但与结合面的波长没有关系;非线性刚度会导致结合面固有频率的下降和振动位移的不对称.     

A size-dependent Kirchhoff micro-plate model based on strain gradient elasticity theory

A size-dependent Kirchhoff micro-plate model is developed based on the strain gradient elasticity theory. The model contains three material length scale parameters, which may effectively capture the size effect. The model can also degenerate into the modified couple stress plate model or the classical plate model, if two or all of the material length scale parameters are taken to be zero. The static bending, instability and free vibration problems of a rectangular micro-plate with all edges simple supported are carried out to illustrate the applicability of the present size-dependent model. The results are compared with the reduced models. The present model can predict prominent size-dependent normalized stiffness, buckling load, and natural frequency with the reduction of structural size, especially when the plate thickness is on the same order of the material length scale parameter.     

Discrete element method for high-temperature spread in compacted powder systems

The discrete element method is applied to investigate high-temperature spread in compacted metallic particle systems formed by high-velocity compaction. Assuming that heat transfer only occurs at contact zone between particles, a discrete equation based on continuum mechanics is proposed to investigate the heat flux. Heat generated internally by friction between moving particles is determined by kinetic equations. For the proposed model, numerical results are obtained by a particle-flow-code-based program. Temperature profiles are determined at different locations and times. At a fixed location, the increase in temperature shows a logarithmic relationship with time. Investigation of three different systems indicates that the geometric distribution of the particulate material is one of the main influencing factors for the heat conduction process. Higher temperature is generated for denser packing, and vice versa. For smaller uniform particles, heat transfers more rapidly.     

Micromorphic modeling of granular dynamics

This paper provides micromorphic modeling of a granular material. Micromorphic modeling treats an individual particle as a microelement and the particle composition in a representative volume element as a macroelement. By specifying the volume of a macroelement, continuum volume-type quantities such as mass density, body force, body couple, kinetic energy density, internal energy density, specific heat supply, etc., are determined by taking the averages of their discrete counterparts in a macroelement. The discrete expressions for the divergence of surface-type quantities (fluxes) are obtained with the help of discrete–continuum analogy for the discrete balance equations. We demonstrate that the discrete formulation of stress tensor in the dynamic condition, which involves both contributions from body forces and relative particle accelerations in a macroelement, can be simply expressed in terms of contact forces and branch vectors. This study constructs complete discrete-type and continuum-type balance equations for a granular material in a macroelement and at a macroscopic point, using the discrete–continuum correspondence for these field quantities.     

Contact stresses: a short survey of models and methods of computations

Contact stresses are identified as normal and tangential forces between contacting solids. The normal stresses are modeled using unilateral and complementary conditions, elastic response and normal compliance. Friction laws describe the tangential traction. Friction of materials depends on pressure, sliding velocity, surface temperature, time of contact, surface roughness and presence of wear debris. Phenomenological, micro-mechanical and atomic-scale models as well as non-classical models of anisotropic and heterogeneous friction are important steps in the development of friction modeling. Sophisticated friction models are desirable in vibrating systems, materials processing, rolling contacts, rubber and polymers, geomechanics, bioengineering and living systems. Main numerical methods in contact mechanics are: finite element method, boundary element method and discrete element method. To include specific contact constraints, the following computing techniques are applied: Lagrange multipliers, penalty function, perturbated and augmented Lagrangian methods, mathematical programming methods. The advances of adhesion and impact modeling are outlined in this paper.     

On the compression of a stack of truncated elastomeric cones as a nonlinearly responsive spring

In an effort to construct a design tool for a mechanical spring featuring highly nonlinear spring stiffness, compression of truncated elastomeric cones has been studied using nonlinear finite element analyses involving neo-Hookean material law and contact elements. Series of finite element models of various geometric aspect ratios of truncated cones were calculated to form a fundamental database of the design tool. It was found that the compressive stiffness of the rubber cone can be non-dimensionalized with respect to the elastic modulus and a characteristic length of the cone. While the stiffness of the truncated rubber cone appears more linear between 0 and 5% of the compression ratio, the stiffness increases exponentially with progressing compression at higher compression ratios. Regression equations of the non-dimensional axial force and spring stiffness were obtained with reasonable accuracy, compared with the original finite element data.     

基于块体离散元法的多面体接触重叠算法

近年来随着非连续介质力学方法的发展,离散元法成为了颗粒材料的物理力学特性研究的重要工具。而对于任意形状复杂块体,精确的接触检测算法一直是离散元法的难题之一。本文基于几何对偶理论,在耦合模拟器CoSim (Coupling Simulator)的块体离散元框架下开发了多面体接触重叠算法。该算法融合了GJK (Gilbert-Johnson-Keerthi)和快速凸包等算法,能够准确计算重叠多面体,进而从接触重叠体中提取接触点、法线方向、接触面积和嵌入深度等接触特性。通过颗粒碰撞测试和砌体结构破坏试验,验证了该算法的准确性和适用性,能够解决多面体的接触问题。     

粗粒土离散元非球趋真颗粒模型

颗粒形状对粗粒土的物理力学特性有着显著的影响。离散元法广泛应用于粗粒土宏观物理力学特性的细观机理研究。为了考虑颗粒形状的影响,亟待发展计算高效的离散元非球趋真颗粒模型。本文基于X射线CT扫描技术并结合数字图像处理技术,对光滑和棱角性两类典型粗粒土（鹅卵石与碎石）进行三维重构,并提出了两类趋真颗粒模型,分别采用扩展超椭球模型和球多面体模型进行趋真逼近;开展了两类颗粒试样的3D打印和单轴压缩试验,分析了配位数和局部孔隙率分布等细观特性;基于离散元开源程序SudoDEM开展了两类试样的离散元模拟,并将模拟细观分析结果与物理试验进行了对比。结果表明,提出的两类趋真颗粒模型能够较好地对粗粒土颗粒进行离散元建模。