Particle shape consideration in numerical simulation of assemblies of irregularly shaped particles

The mechanical behavior of granular materials depends much on the shape of the constituent particles. Therefore appropriate modeling of particle, or grain, shape is quite important. This study employed the method of direct modeling of grain shape (Matsushima & Saomto, 2002), in which, the real shape of a grain is modeled by combining arbitrary number of overlapping circular elements which are connected to each other in a rigid way. Then, accordingly, a discrete-element program is used to simulate the assembly of grains. In order to measure the effects of grain shape on mechanical properties of assembly of grains, three types of grains—high angular grains, medium angular grains and round grains are considered where several biaxial tests are conducted on assemblies with different grain types. The results show that the angularity of grains greatly affects the behavior of granular soil.     

NUMERICAL SIMULATION STUDY ON ROCK BREAKING MECHANISM AND PROCESS UNDER HIGH PRESSURE WATER JET

The numerical simulation method to study rock breaking process and mechanism under high pressure water jet was developed with the continuous mechanics and the FEM theory. The rock damage model and the damage-coupling model suited to analyze the whole process of water jet breaking rock were established with continuum damage mechanics and micro damage mechanics. The numerical results show the dynamic response of rock under water jet and the evolvement of hydrodynamic characteristic of jet during rock breaking is close to reality, and indicates that the body of rock damage and breakage under the general continual jet occurs within several milliseconds, the main damage form is tensile damage caused by rock unload and jet impact, and the evolvement of rock damage shows a step-change trend. On the whole, the numerical results can agree with experimental conclusions, which manifest that the analytical method is feasible and can be applied to guide the research and application of jet breaking rock theory.     

动力荷载对围压卸载下岩石动态变形模量的影响

利用SHPB岩石动静组合加载实验系统,研究在不同轴压水平下围压以1 MPa/s速率卸载至预加值50%时矽卡岩受频繁冲击作用的动态变形模量变化规律。实验结果表明：高轴压促使岩石内部微裂纹萌发与扩展,降低了岩石抵抗外部冲击的能力。围压的侧向约束阻碍岩石内部裂纹的横向扩展,但在围压卸载时会加剧岩石内部的损伤,这是由于高轴压下,围压卸载导致岩石内部应力重新分布。轴压与围压卸载共同影响着冲击作用下的岩石动态变形模量,通过岩样在冲击荷载下的能量耗散分析岩石动态变形模量的变化规律,有助于了解深部岩体开挖的破坏机制。

     

频繁动力扰动对围压卸载中高储能岩体的动力学影响

基于冬瓜山铜矿深部巷道围岩开挖过程中面临的高应力和频繁爆破扰动问题,利用改进的SHPB动静组合加载系统,开展了频繁动力扰动对围压卸载中高储能岩体动力学影响的研究。研究结果表明,围压卸载中的矽卡岩受到动力扰动时,其动态峰值应力和弹性模量随动力扰动次数非线性变化。围压卸载中的高储能矽卡岩受到动力扰动时会释放能量。轴压促使岩样内微裂隙轴向发育,造成岩样抵抗动力扰动能力减弱;围压减缓岩样内微裂隙轴向发育,造成岩样抵抗动力扰动能力增强。动力扰动对微裂隙扩展有促进作用,使围压卸载中的岩样由拉伸破坏向剪切破坏转变。     

Numerical simulation of polygonal particles moving in incompressible viscous fluids

A two-dimensional coupled lattice Boltzmann immersed boundary discrete element method is introduced for the simulation of polygonal particles moving in incompressible viscous fluids. A collision model of polygonal particles is used in the discrete element method. Instead of a collision model of circular particles, the collision model used in our method can deal with particles of more complex shape and efficiently simulate the effects of shape on particle–particle and particle–wall interactions. For two particles falling under gravity, because of the edges and corners, different collision patterns for circular and polygonal particles are found in our simulations. The complex vortexes generated near the corners of polygonal particles affect the flow field and lead to a difference in particle motions between circular and polygonal particles. For multiple particles falling under gravity, the polygonal particles easily become stuck owing to their corners and edges, while circular particles slip along contact areas. The present method provides an efficient approach for understanding the effects of particle shape on the dynamics of non-circular particles in fluids.     

深埋地下洞室断裂型岩爆机理的数值模拟

针对深部岩体中由断层、节理等不连续性结构面引发的岩爆地质灾害, 根据深埋地下隧洞中潜在发震断裂的分布特征和几何形态建立数值分析模型, 采用离散元单元法模拟存在刚性平直断裂的深部围岩的开挖响应, 并分别考察开挖接近并通过断裂附近时围岩应力状态的变化特征。通过探讨断裂的存在对围岩应力状态改变的作用机理, 揭示出断裂型岩爆是开挖面附近一定范围内存在的断裂构造在高应力作用下发生错动, 导致能量突然释放, 对围岩造成强烈冲击作用的结果, 基本与地震的断层“粘滑”机制相类似。     

Numerical simulation of tetrahedral particle mixing and motion in rotating drums

A regular tetrahedron is the simplest three-dimensional structure and has the largest non-sphericity. Mixing of tetrahedral particles in a thin drum mixer was studied by the soft-sphere-imbedded pseudo-hard particle model and compared with that of spherical particles. The two particle types were simulated with different rotation speeds and drum filling levels. The Lacey mixing index and Shannon information entropy were used to explore the effects of sphericity on the mixing and motion of particles. Moreover, the probability density functions and mean values and variances of motion velocities, including translational and rotational, were computed to quantify the differences between the motion features of tetrahedra and spheres. We found that the flow regime depended on the particle shape in addition to the rotation speed and filling level of the drum. The mixing of tetrahedral particles was better than that of spherical particles in the rolling and cascading regimes at a high filling level, whereas it may be poorer when the filling level was low. The Shannon information entropy is better than the Lacey mixing index to evaluate mixing because it can reflect the real change of flow regime from the cataracting to the centrifugal regime, whereas the mixing index cannot.     

多孔介质在压力梯度作用下的热质耦合数值模拟

多孔介质干燥导致热质耦合传输过程。本文基于连续介质力学的宏观尺度,对多孔介质的热、湿和气三者耦合迁移进行数值模拟,研究压力梯度对热质传输的影响。多孔介质传质机理主要为水汽和空气的对流和扩散传输、吸附水在含湿量梯度作用下的自由扩散和其在温度梯度即Soret效应驱动下的流动。采用Galerkin加权余量的有限元方法,提出了...     

Numerical and experimental investigations of micromechanical processes during wire sawing

In the following paper the mechanisms of the slurry-based wire sawing process at micro level are studied by numerical simulations. For the simulation a new approach in simulating the wire sawing process is presented. A discrete element method (DEM) with polyhedral particles is used to simulate the movement of the abrasive particles as well as the wire during the process. The results are verified by direct observations of the sawing channel of an experimental setup. The presented numerical model is applied for a detailed study of the influence of different process parameters.     

Numerical simulation study on particle breakage behavior of granular materials in confined compression tests

In order to study the fragmentation law, the confined compression experiment of granular assemblies has been conducted to explore the particle breakage characteristic by DEM approach in this work. It is shown that contact and contact force during the loading process gradually transform from anisotropy to isotropy. Meanwhile, two particle failure modes caused by different contact force states are analyzed, which are single-through-crack failure and multi-short-crack failure. Considering the vertical distribution of the number of cracks and the four characteristic stress distributions (the stress related to the maximum contact force, the major principal stress, the deviatoric stress and the mean stress), it is pointed out that the stress based on the maximum contact force and the major principal stress can reflect the distribution of cracks accurately. In addition, the size effect of particle crushing indicates that small size particles are prone to break. The lateral pressure coefficient of four size particles during the loading process is analyzed to explain the reason for the size effect of particle breakage.     

Study of particle rearrangement during powder compaction by the Discrete Element Method

This paper presents simulations of cold isostatic and closed die compaction of powders based on the Discrete Element Method. Due to the particulate nature of powders, densification of the compact proceeds both through the plastic deformation at the particle contact and the mutual rearrangement of particles. The relative weight of each mechanism on the macroscopic deformation process depends on the contact law, the relative density, and the type of stress exerted on the particles (shear or pressure). 3D computer simulations have been carried out to investigate the role of these parameters on the deformation mechanisms of powder compacts. The effect of rearrangement is studied by comparing simulations that use a homogeneous strain field solution for which local rearrangement is omitted and simulations that include local rearrangement. It is shown that local rearrangement has some effect on average quantities such as the average coordination number, the average contact area and the macroscopic stress. The effect on averaged quantities is much stronger for closed die compaction than for isostatic compaction. However the main effect of local rearrangement is to widen the distribution of the parameters that define the contact (contact area in particular). The results of these simulations are compared to available experimental data and to statistical models that use a homogeneous strain field assumption.     

贾木那大桥水压定向爆破及其数值模拟

为了分析及研究薄壁空腔桥墩的水压爆破机理,以孟加拉国贾木那大桥水压定向爆破拆除工程为研究对象,采用显示动力学分析软件模拟了桥墩的水压爆破破碎过程,根据数值模拟结果,得出了水压爆破主要是利用爆炸冲击波、爆轰气体、反射拉伸波及高速水流对桥墩壁面进行破坏,在此基础上,探讨了桥墩水压定向爆破的药包布置方式、起爆顺序等技术问题。     

Equivalent boundary model of lunar soil drilling simulation by DEM

Aiming to solve the computational cost problem in the discrete element simulation for lunar soil drilling sampling, an equivalent boundary method was proposed. A high-accuracy DEM model of lunar soil was established firstly. As the novel alterable constitutive law, the accuracy of the model was verified to meet the performance of real lunar soil very much both in shear strength indices and elastic–plastic behavior. A common drill bit in the geological exploration field for sampling soil was chosen as the simulation object. In preanalysis, it was known that with the increase of drilling depth, the stress concentration area was always near the drill bit, while the affected area of the lunar soil was a cylindrical area around the drill pipe, which extended towards the drilling direction instead of extending around it. Then a big boundary drilling simulation scene was established to investigate the flow direction of lunar soil particles. The motion law of particles and the velocity field information were obtained, and a U-shape chain was described around the drill bit. Finally an equivalent boundary was set near the U-shaped chain, and the size was determined by comparing the soil stress in the fierce collision zone and around the reference boundary. This method could be a reference for other lunar soil drilling researches with other drills of different sizes.     

Numerical simulation and testing verification of the interaction between track and sandy ground based on discrete element method

Tractive effort of tracked vehicles plays an important role in military and agricultural fields. In order to solve the problem of low precision in numerical simulation of the interaction between track and sandy ground, a systematic and accurate discrete element modeling method for sandy road was proposed. The sandy ground was modeled according to the mechanical parameters measured by soil mechanics tests. The interaction coefficients of sandy soil were measured by the repose angle test and triaxial compression test combined with the corresponding simulation. On this basis, a discrete element interaction model of track-sandy ground was established, which can be used to test the tractive effort of track. Numerical simulation calculation of track model at different speeds was carried out, and the simulation results were compared with the results of indoor soil bin test for verification. The verification results show that the interaction between track and sandy ground based on DEM simulation is consistent with the actual soil bin test. The discrete element modeling method in this paper can be used to model the track and sandy ground accurately, and the simulation model can be used to test the tractive effort of tracked vehicle.     

Numerical simulation of a fluid-saturated soil under a dynamically loaded pile

The behaviour of the soil under a dynamically loaded pile toe is studied. The soil is modelled as a fluid saturated porous continuum. The constitutive behaviour of the solid skeleton is described by the elasto-plastic model of Drücker-Prager. The wave propagation is simulated with a dynamical finite-element program.A two-phase model of soil gives useful information about effective stress and pore pressure in the soil. In saturated soil the main wave under the pile toe propagates more downards than in dry soil, due to the higher compressional stiffness in saturated soil. The plastic zone under the pile toe propagates with the velocity of the fast compressional wave. The pore fluid influences the plasticity strongly and can be expected to affect pile driving too.The distribution of effective stress and pore pressure under the pile toe depends on the permeability of the soil and cannot be calculated uniquely from a single-phase calculation. Therefore, a nonlinear soil cannot be modelled correctly as a conventional single-phase material.     

内压对电场下晶内微裂纹演化的影响

基于表面扩散的经典理论及其弱解描述,建立曲率、内压和电场共同作用下的有限单元法,分析了金属材料内部晶内微裂纹的演化过程。对于形态比为β的微裂纹,存在一临界电场值χc,当χ<χc时,微裂纹逐渐圆柱化;当χ>χc时,微裂纹分节为上下或左右两个较小的微裂纹。内压加速有助于微裂纹分节。对于β>1的微裂纹,当β和χ一定时存在一临界内压值qc,当q>qc时,随着电场和内压的增大会率先在裂纹的右端发生裂纹扩展。     

用边界元法分析轧制压力和摩擦力

本文给出了一种新的求解轧制力和摩擦力分布的方法，利用测量手段对轧制过程中的轧件进行测量，然后再利用边界元法对轧辊进行分析，从而求解出轧制中轧制力和摩擦力。这种方法避开了求解轧制力对轧件进行弹塑性分析，也解决了用测力销测量带来的困难，使轧制问题的分析求解得到了简化并提高了求解精度。     

Numerical simulation of natural and mixed convection flows by Galerkin‐characteristic method

A numerical investigation is performed to study the solution of natural and mixed convection flows by Galerkin‐characteristic method. The method is based on combining the modified method of characteristics with a Galerkin finite element discretization in primitive variables. It can be interpreted as a fractional step technique where convective part and Stokes/Boussinesq part are treated separately. The main feature of the proposed method is that, due to the Lagrangian treatment of convection, the Courant–Friedrichs–Lewy (CFL) restriction is relaxed and the time truncation errors are reduced in the Stokes/Boussinesq part. Numerical simulations are carried out for a natural convection in squared cavity and for a mixed convection flow past a circular cylinder. The computed results are compared with those obtained using other Eulerian‐based Galerkin finite element solvers, which are used for solving many convective flow models. The Galerkin‐characteristic method has been found to be feasible and satisfactory. Copyright © 2006 John Wiley & Sons, Ltd.     

基于离散元的杆件断裂全过程数值模拟

杆件的断裂会涉及到大变形、非线性以及不连续等问题,通常的数值计算方法模拟这种复杂力学行为具有局限性。本文基于颗粒离散元法DEM,将接触粘结处的分布式弹簧用梁纤维进行等效,提出了一种适于结构弹塑性分析的DEM纤维梁模型,然后在此基础上构建了构件断裂模拟算法以及纤维破环准则。将该模型应用于悬臂梁结构,模拟了悬臂梁从弹性到弹塑性阶段,再到断裂破坏的全过程,数值模拟得到的结构响应和截面开裂破坏形态均较合理。最后将该方法应用于单层网壳倒塌破坏模拟,并与网壳振动台倒塌试验进行对比,结果表明,数值模拟得到的杆件断裂过程及结构倒塌模式与试验现象一致,验证了该模型的正确性和适用性。