共查询到20条相似文献，搜索用时 343 毫秒
1.
Erfan G. Nezami Youssef M. A. Hashash Dawei Zhao Jamshid Ghaboussi 《国际地质力学数值与分析法杂志》2006,30(8):783-801
With the increasing demand for discrete element simulations with larger number of particles and more realistic particle geometries, the need for efficient contact detection algorithms is more evident. To date, the class of common plane (CP) methods is among the most effective and widely used contact detection algorithms in discrete element simulations of polygonal and polyhedral particles. This paper introduces a new approach to obtain the CP by employing a newly introduced concept of ‘shortest link’. Among all the possible line segments that connect any point on the surface of particle A to any point on the surface of particle B, the one with the shortest length defines the shortest link between the two particles. The perpendicular bisector plane of the shortest link fulfils all the conditions of a CP, suggesting that CP can be obtained by seeking the shortest link. A new algorithm, called shortest link method (SLM), is proposed to obtain the shortest link and subsequently the CP between any two polyhedral particles. Comparison of the analysis time between SLM and previously introduced algorithms demonstrate that SLM results in a substantial speed up for polyhedral particles contact detection. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
2.
Discrete element modeling of dynamic behaviors of railway ballast under cyclic loading with dilated polyhedra 
下载免费PDF全文
下载免费PDF全文 The dynamic behaviors of railway ballast under cyclic loading are simulated with discrete element method (DEM). Dilated polyhedra are constructed based on the Minkowski sum operator in order to resemble the irregular shapes of ballast particles. The polyhedral particle generation, contact detection between particles and contact laws are presented. Ballast box tests with periodic lateral boundaries are conducted to simulate the dynamics of the sleeper and ballast particles. The settlement and effective stiffness of ballast bed are investigated under cyclic loadings with five distinct frequencies. The settlement of ballast bed is significant in the first several cycles and increases with the number of cycles gradually. The higher frequency loading generates larger displacement in the same simulation time. The effective stiffness of ballast bed increases gradually. To study the effect of particles' sharpness, dilated polyhedra with different dilating radii and spherical particles are also developed. Simulation results show the sharper the ballast particles are, the smaller the produced settlement. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
3.
Irregularly shaped (IRS) particles widely exist in many engineering and industrial fields. The macro physical and mechanical properties of the particle system are governed by the interaction between the particles in the system. The interaction between IRS particles is more complicated because of their complex geometric shape with extremely irregular and co‐existed concave and convex surfaces. These particles may interlock each other, making the sliding and friction of IRS particles more complex than that of particles with regular shape. In order to study the interaction of IRS particles more efficiently, a refined method of constructing discrete element model based on computed tomography scanning of IRS particles is proposed. Three parameters were introduced to control the accuracy and the number of packing spheres. Subsequently, the inertia tensor of the IRS particle model was optimized. Finally, laboratory and numerical open bottom cylinder tests were carried out to verify the refined modeling method. The influence of particle shape, particle position, and mesoscopic friction coefficient on the interaction of particles was also simulated. It is noteworthy that with the increase of mesoscopic friction coefficient, the fluidity of IRS particle assembly decreases, and intermittent limit equilibrium state may appear. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
4.
Interaction between solid particles and fluid is of fundamental interest to scientists and engineers in many different applications—cardiopulmonary flows, aircraft and automobile aerodynamics, and wind loading on buildings to name a few. In geomechanics, particle shape significantly affects both particle-particle and particle-fluid interaction. Herein, we present a generalized method for modeling the interaction of arbitrarily shaped polyhedral particles and particle assemblages with fluid using a coupled discrete element method (DEM) and lattice Boltzmann method (LBM) formulation. The coupling between DEM and LBM is achieved through a new algorithm based on a volume-fraction approach to consider three-dimensional convex polyhedral particles moving through fluid. The algorithm establishes the interaction using linear programming and simplex integration and is validated against experimental data. This approach to modeling the interaction between complex polyhedral particles and fluid is shown to be accurate for directly simulating hydrodynamic forces on the particles. 相似文献
5.
基于颗粒流的堆石料三轴剪切试验研究 总被引：2，自引：0，他引：2
近年来频发的地质灾害使得深入研究大坝堆石体受力变形机制显得愈发迫切和重要。针对堆石颗粒具有形状极其不规则、咬合作用突出的特点,编写了一个简单易用的程序来随机生成形状不规则的数值颗粒。以正四面体为核,按照晶胞繁衍的方式,生成堆石的仿真颗粒集合体。模拟排水剪切条件下的大三轴试验,通过分析数值试验过程中细观参量的变化规律,探讨三轴剪切条件下堆石体变形的细观机制。分析表明,按照晶胞繁衍法生成的数值颗粒与圆颗粒相比能更好地模拟堆石,由颗粒簇形成的咬合力能形成更真实的力-变形关系;数值试验得到的应力-应变曲线和体变曲线与室内三轴试验结果基本一致;颗粒的破碎速率是联系堆石体宏、细观力学性质的重要纽带,通过分析堆石体颗粒破碎4个阶段细观参量的变化规律,深化了对堆石体变形机制的认识。 相似文献
6.
结合边界元法和离散元法,提出一种可以进行计算颗粒内部应力和破碎路径的方法。该方法利用离散元法求解颗粒的相互作用和每个颗粒上的荷载。然后利用边界元法计算颗粒的应力分布,为了实现动态平衡,将颗粒的加速度视为恒定大小的体力。但体力导致边界积分方程中出现域积分,故采用直线积分法将域积分转化为边界积分,以保证边界元法降维的优势。为了提高边界元的计算效率,对于几何形状相似的颗粒,以其中一个颗粒作为模板颗粒,只需要计算模板颗粒在局部坐标系中的系数矩阵,其他相似颗粒可以通过局部和全局坐标系之间的映射获得。在得到应力后,基于Hoek-Brown准则来判断颗粒是否破碎。此外,将破坏路径简化为直线,并采用最小二乘法拟合得到破坏路径。 相似文献
7.
8.
Gowri Srinivasan Elizabeth Keating John David Moulton Zora V. Dash Bruce A. Robinson 《Computational Geosciences》2012,16(3):551-563
A convolution-based particle tracking (CBPT) method was recently developed for calculating solute concentrations (Robinson
et al., Comput Geosci 14(4): 779–792, 2010). This method is highly efficient but limited to steady-state flow conditions. Here, we present an extension of this method
to transient flow conditions. This extension requires a single-particle tracking process model run, with a pulse of particles
introduced at a sequence of times for each source location. The number and interval of particle releases depends upon the
transients in the flow. Numerical convolution of particle paths obtained at each release time and location with a time-varying
source term is performed to yield the shape of the plume. Many factors controlling transport such as variation in source terms,
radioactive decay, and in some cases linear processes such as sorption and diffusion into dead-end pores can be simulated
in the convolution step for Monte Carlo-based analysis of transport uncertainty. We demonstrate the efficiency of the transient
CBPT method, by showing that it requires fewer particles than traditional random walk particle tracking methods to achieve
the same levels of accuracy, especially as the source term increases in duration or is uncertain. Since flow calculations
under transient conditions are often very expensive, this is a computationally efficient yet accurate method. 相似文献
9.
黏土颗粒形态不仅反映黏土的矿物组分,更是影响其物理力学性质的重要因素之一。为了研究物质组成对软黏土微宏观性质的影响,采用离散元方法对不同颗粒形态的软黏土试样进行三轴压缩模拟试验。首先,基于扫描电镜图像量化颗粒形态,对天然状态下黏土颗粒的方向角和凹凸度进行统计,引入球度和凹凸度作为颗粒形态的特征参数;然后,基于原生矿物的单粒结构和黏土矿物的片状结构特征,构造球体单粒及圆柱体、正方体、长方体的片状簇体;最后,基于三轴试验离散元模拟方法,分析软黏土颗粒形态对其宏观力学及微观特性的影响。结果表明:片状颗粒试样比球体颗粒试样的初始模量高,抗剪强度大,随加载其排列趋于水平向分布;加载初期,颗粒球度对初始弹性模量影响较明显,初始弹性模量随着球度增大而逐渐减小;加载后期,颗粒凹凸度对抗剪强度指标影响作用逐渐凸显,试样内摩擦角和黏聚力随着凹凸度增大而逐渐减小;微观结构上,颗粒形状对颗粒位移和旋转也有较大影响。 相似文献
10.
常见的地质灾害如滑坡、泥石流、岩崩等通常都涉及不同形状的颗粒物质运动,这些形状不同的颗粒又多具有不同的尺寸和含量。基于典型的颗粒柱坍塌试验,首先根据试验方法确定了离散元模拟所需的各项参数,然后采用随机多面体方法生成了可控制长细比的大颗粒,利用离散元法就不同大颗粒含量下形态变化对二元颗粒柱坍塌特性的影响开展研究,研究结果表明:(1)利用离散元法可以较好地重现室内试验中小球和多面体组成的二元颗粒系统的颗粒柱坍塌过程;(2)在不同长细比的不规则大颗粒和小球组成的二元颗粒柱系统中,当大颗粒含量高于临界含量值20%时,二元颗粒柱坍塌持续的时间随非球形大颗粒长细比的增加而增加;(3)在不同长细比的不规则大颗粒和小球组成的二元颗粒柱中,当大颗粒含量高于临界含量值20%时,在相同百分比的大颗粒含量下,大颗粒长细比的增加会提高大颗粒平均配位数以及降低颗粒的运动能力,大颗粒间形成更强的互锁作用,降低了颗粒柱的整体流动性,使其最终堆积高度更高、最大跑出距离更短以及更小的归一化动能峰值。(4)在不同长细比的不规则大颗粒和小球组成的二元颗粒柱中,小颗粒可以较为明显降低大颗粒间摩擦及互锁作用,增加流动性,降低大骨料形态对坍塌过程的影响。 相似文献
11.
在改进均匀多面体重力场正演公式基础上, 利用二阶张量的坐标变换实现对多面体重力场梯度的求解, 推导了新的多面体重力梯度和磁场的正演公式, 给出了新的统一的均匀多面体重力场、梯度及磁场正演表达式形式, 并用理论模型进行了检验.同时, 应用新的多面体重力场梯度正演公式推导出新的长方体重力场垂直梯度理论表达式.本文给出的均匀多面体重力场、梯度及磁场正演表达式形式统一, 重磁场联合正演中可相互利用其计算过程中的结果, 避免重复计算以提高正演计算效率. 相似文献
12.
根据颗粒离散元Kelvin 接触力计算模型,分析了圆形颗粒体模拟材料力学特性应具备的条件,在此基础上提出了一种新颗粒模型构建方法。该方法首先在复杂模型域内随机生成种子,然后利用相切条件逐步扩展填充整个区域。填充过程中借助局部Delaunay三角化网格控制新颗粒的生成,采用复杂几何体距离控制颗粒与模型边界的相对位置,对靠近模型边界的颗粒进行容忍性优化填充,从而增加模型颗粒与边界的耦合性。同时对模型孔隙进行再填充,保证每个填充颗粒至少与3个颗粒相切,提高了模型内颗粒间的耦合性和模型的密度。最后采用任意多边形控制材料边界,将模型材料的设置简化为判断点是否在多边形内,简化了复杂模型材料属性的设置过程。结果表明：与膨胀颗粒生成法相比,该方法生成模型重叠量小、颗粒间及颗粒-边界相互耦合、填充率高。因此,颗粒黏结力破坏后不会造成飞溢现象,可适用于任意连通域模型的生成,能更好地实现复杂岩土细观介质变形破坏机制的模拟与研究。 相似文献
13.
考虑到颗粒形状对粗粒料的力学特性有重大影响,提出了一种新的表征颗粒形状的方法,即在椭圆上随机选取一系列点连接成多边形颗粒,表征狭长扁平的颗粒。新方法较圆上取点的方法能代表更多类型的颗粒形状,适用范围更广。提出了一种新的粗粒料投放算法,即先缩小颗粒,然采用随机算法将缩小的颗粒投放至给定区域,对颗粒划分好网格后,将颗粒放大到原来的大小,然后采用有限元-离散元(FEM/DEM)方法计算稳定后即生成了相应的试样。通过将上述颗粒生成及投放算法与FEM/DEM结合,应用于粗粒料的数值模拟。分析表明,FEM/DEM是研究粗粒料力学性质的较好方法,对复杂的颗粒形状也可简单建模,且因在颗粒内部划分了有限元网格,复杂的接触判断及接触力计算转化为标准统一的三角形和三角形之间的接触判断及接触力计算,所有的计算均可标准化、统一化。同时因为颗粒是可以变形的,依然保留了连续介质力学中应力和应变的概念,无须像PFC那样需通过测量圆来间接表示某点的应力、应变。最后,通过粗粒料的侧限压缩试验的数值模拟,展现了文中提出的一整套解决方案在模拟粗粒料方面的巨大潜力。 相似文献
14.
Mechanical trapping (or straining) of fine particles is a key mechanism in many filtration systems. For example, the performance of rapid sand filters depends in part on mechanical trapping of larger fine particles, while relying on adsorptive processes to trap very small fine particles and microbes. The ability to trap these particles is directly related to the construction of the packed bed used for filtration in this system. Thus, the ability to model the effect of the inner structure of the packed bed can lead to more efficient design for improved filtration. Because of its significant efficiency, gravitational sphere packing is employed in this work to simulate a bed of mono‐sized randomly packed spheres. The simulated bed provides a way to visualize the pore network and estimate the pore size distribution associated with the void space between particles. Furthermore, by subsequently introducing fine particles into the bed, we evaluate the mass‐rate of fine particles passing through and possibly saturating the packed bed. Results show that fine particles between 15% and 25% of the coarse particle size can be physically strained within the randomly packed bed. These results differ significantly from the results obtained assuming a periodically spaced bed. The technique therefore provides an efficient yet accurate alternative for understanding how fine particles pass through a coarse particulate medium. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
15.
Pb on groundwater particles,Door County,Wisconsin 总被引：4，自引：0，他引：4
N. C. Woo 《Environmental Geology》1994,24(2):150-156
Pb was detected on groundwater particles larger than 0.45m. Pb concentration (micrograms per liter) showed positive correlations with particle concentration (milligrams per liter) at five sampled wells, although the slopes of the best-fit lines were different. X-ray fluorescence scans suggested that the particles primarily consist of K-Al silicate minerals, and that Pb is absorbed onto Fe oxyhydroxide coatings on Mg-rich silicates. Particle compositions varied at wells, but remained consistent at each well through the sampling period of a month. The particulate nature of Pb in groundwater could cause the spatial and temporal variations due to: (1) the unique composition of particles at wells, (2) the irregular pathways for groundwater particles, (3) the unpredictable amount of Pb-carrying particles entering the groundwater system during recharge events, and (4) the irregular sampling of groundwater particles subsequently dispersed along the advective flow path. 相似文献
16.
In order to simulate undrained conditions using the discrete element method, a constant sample volume is often assumed. There are well-recognised problems with these constant-volume triaxial simulations, particularly of dense samples, which inhibit quantitative comparison with laboratory experiments. In this paper, four possible explanations for these problems with conventional constant-volume simulations of ideal spherical particles are explored, each of which has a physical basis: particle crushing, the presence of highly compressible air within the sample, or the reduction in stiffness due to particle surface asperities or non-spherical particle shapes. These options are explored independently and in combination through implementation in the open-source LAMMPS code. In situations where a significant amount of particle crushing occurs, it is important to incorporate this in the simulations so that stresses are not over-estimated. There is experimental evidence that irregular particles have lower Young’s moduli than the Hertzian spheres often used in DEM. In the absence of particle crushing, the most effective method to achieve more realistic stress–strain responses is to reduce the particle shear modulus substantially. This approach has the added computational benefit of enabling an increase in the simulation time-step.
相似文献17.
Separation of small nonferrous particles using a two successive steps eddy-current separator with permanent magnets 总被引：1，自引：0，他引：1
The paper presents a method for separating the small metallic nonferrous particles from two component metallic nonferrous mixtures using a new type of dynamic eddy-current separator with permanent magnets. The so-called two successive steps eddy-current separator (TSECS) consists of a horizontal rotary drum covered with permanent magnets, alternately N–S and S–N oriented. The separation process takes place in two stages, first the strongly conducting particles are separated on the upper part of the drum, and then the remaining undecided and poorly conducting particles are separated at the lower part of the magnetic drum. The experimental results and comments regarding the values obtained for separation efficiency, grade and recovery for wastes consisting in Cu–Pb and Cu–Al mixtures are given. The obtained results are in good agreement with the theoretical analysis. 相似文献
18.
Smoothed particle hydrodynamics (SPH) is a meshfree, Lagrangian particle method which has advantages in handling solids with extremely large deformation. Like any other numerical methods, cares must be taken to ensure its desirable accuracy and stability through considering several correction techniques in calculation. The selection of values for parameters in those correction approaches is a key step in SPH simulation, which is always difficult for new beginners to deal well with effectively. This paper examines the common inconsistency and instability problems in SPH method and studies its computational efficiency when applied to hydrodynamics problems with material strength like soil column collapse. We analyzed in detail how the correction techniques mitigate these inconsistency and instability problems. Also, the numerical testing results associate with different values for the parameters used in the correction techniques are provided for better understanding the influence of these parameters and for finding out the desirable values. It is found that (1) the SPH method is easily subjected to an inconsistency problem in the boundary area due to the boundary deficiency, and it can be treated well by adopting “virtual particles” contributing to the particle summations. (2) The numerical oscillation in SPH simulation can be mitigated effectively by artificial viscosity with the suggested parameter values. (3) The tension cracking treatment, artificial viscosity and artificial stress work well in removing the tensile instability problem in SPH method. In addition, the nearest neighboring particle searching (NNPS) algorithm, spacing ratio, smoothing length and time step influence the efficiency and accuracy of SPH method significantly. It is shown that SPH method with suggested parameters values can produce a very good result compared with the experimental result. 相似文献
19.
自然条件下,颗粒介质大多以非规则单元形态存在。非规则几何形态对颗粒介质的宏观力学性能有很大影响。针对颗粒单元的不同几何形态,采用团颗粒单元对离散介质的直剪试验过程进行了离散元数值计算,详细地讨论了颗粒形态对离散介质剪切强度的影响。该非规则颗粒由不同形态、不同数目、镶嵌尺寸、组合方位和颗粒大小的球形颗粒进行随机构造,其在局部与整体坐标之间的转动、力矩和方位关系通过4元素方法进行确定,基本球体颗粒之间的作用力采用具有Mohr-Coulomb摩擦定侓的Hertz-Mindlin非线性接触模型,并考虑了非线性法向粘滞力的影响。通过构造7种具有相同的质量概率分布的不同形态的团颗粒,在不同法向应力下,对团颗粒的直剪试验进行了离散元模拟,分析了不同形态颗粒的剪切强度。通过对不同形态颗粒介质剪切强度的数值分析,进一步揭示了非规则颗粒间的咬合互锁效应,为分析非规则颗粒的宏观动力特性提供了依据。 相似文献
20.
通过定义距离势函数,提出一种适用于空间任意凸多面体单元的三维距离势函数离散单元法。该方法采用归一化的计算方式,将势函数表征为接触体间的距离函数,并基于此建立接触力计算方程,明确了势函数的物理意义,使接触力计算更加合理,无需对各种可能的接触形式进行差异化处理。新方法克服了原有势函数物理意义不明、接触力计算受单元形式影响等重要缺陷,并突破了四面体单元的限制,可采用空间任意凸多面体单元。通过若干算例说明新方法的正确性和有效性。数值模拟结果表明新方法能够很好处理空间任意多面体单元复杂接触变换过程,准确处理复杂非连续介质的运动过程。 相似文献