考虑颗粒形状和破碎的胶结钙质砂力学行为离散元模拟研究

质砂作为一种建筑材料,近年来广泛应用于我国南海岛礁工程建设中。本文通过建立考虑钙质砂真实颗粒形状和颗粒破碎的胶结钙质砂离散元模型,研究了二维剪切条件下试样的宏微观力学行为,包括应力-应变行为、颗粒破碎、胶结破坏、位移场和裂纹随剪应变的演化规律,讨论了颗粒形状、颗粒粒径范围、颗粒强度和水泥胶结强度对胶结钙质砂力学行为的影响规律。结果表明,钙质砂颗粒粒径区间越宽,胶结钙质砂的强度越高。同一级配条件下,考虑真实颗粒形状的胶结钙质砂试样比圆颗粒试样的强度更高,试样总体颗粒破碎率也更高。钙质砂颗粒的强度越高,胶结钙质砂的性能越好。但是提高水泥的强度对胶结钙质砂力学性能的影响并不显著。本文的研究结果可为实际工程中钙质砂的加固提供理论依据。     

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

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

Predicting the Dynamic Compressive Strength of Carbonate Rocks from Quasi-Static Properties

In this study, the split Hopkinson pressure bar testing method was used to quantify the dynamic strength characteristics of rocks with short cylindrical specimens. Seventy dynamic compression tests were conducted on 10 different carbonate rocks with the split Hopkinson pressure bar apparatus. Experimental procedure for testing dynamic compressive strength and elastic behaviour of rock material at high strain rate loading is presented in the paper. Pulse-shaper technique was adopted to obtain dynamic stress equilibrium at the ends of the sample and to provide nearly a constant strain rate during the dynamic loading. In addition to dynamic tests, the physical properties covering bulk density, effective porosity, P-wave velocity and Schmidt hardness of rocks, and mechanical properties such as quasi-static compressive strength and tensile strength were determined through standard testing methods. Multiple linear regression analyses were carried out to investigate the variation of dynamic compressive strength depending on physical and mechanical properties of rocks and loading rate. A three parameter model was found to be simple and provided the best surface fit to data. It was found that dynamic compressive strength of rocks increases with increase in loading rate and/or increase in rock property values except porosity. Statistical tests of regression results showed that quasi-static compressive strength and Schmidt hardness are most significant rock properties to adequately predict the dynamic compressive strength value among the other properties. However, P-wave velocity, quasi-static tensile strength of rocks could also be used to estimate the dynamic compressive strength value of rocks, as well, except the bulk density and effective porosity.     

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

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

Zr基非晶合金JH-2模型的构建及应用

为构建Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金在高压、大应变、高应变率状态下的材料模型,采用根据实验数据理论推导和数值模拟对比反馈的方法,对材料的Johnson-Holmquist本构模型（JH-2模型）参数进行了研究:材料的静水压力-体应变关系通过平板冲击实验数据和理论推导得到;无损材料强度与应变、应变率的关系通过轴向压缩实验数据确定;材料损伤参数与破碎材料强度参数的关系通过平板冲击实验数据确定;破碎材料强度参数通过数值模拟与实验结果对比的反馈法得到。将材料模型应用于平板冲击和破片侵彻的数值模拟,通过数值模拟与实验结果对比的方式,验证材料模型的准确性。结果表明,平板冲击实验中,材料的自由面粒子速度曲线与数值模拟结果吻合度较高;破片侵彻实验中,破片对钢靶的侵彻深度、开坑孔径与数值模拟结果的一致性较好,构建的材料模型较准确反映了材料的动态力学特性。     

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

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

CuAlNi单晶形状记忆合金压缩特性的实验研究

具有相变伪弹性特性的CuAlNi单晶是目前应用最广泛的形状记忆合金之一。这种材料被广泛应用在工程、生物和医学科学等领域。由于单晶是各向异性，没有多晶中晶粒之间的相互作用，因而在特定晶向上的力学性能稳定。但是这种材料的一些基本性质，如压缩状态下马氏体的发生、生长和传播等还没有人详细研究。本文主要研究CuAlNi单晶在特定晶向上的变形过程，它的应力—应变特性，并利用特殊的显微成像系统首次获得了沿[110]方向压缩时二维马氏体的发生、生长和传播过程。     

基于水平集-离散元方法的球谐函数颗粒材料缓冲性能分析

在自然环境与工业领域中,颗粒材料是一种常见的缓冲材料,其中大量形态各异的非球形颗粒表现出复杂的力学特性并应用于不同工程领域。本文采用球谐函数构造不同球面度和表面凹凸特性的非规则颗粒,通过水平集方法计算球谐函数颗粒间的接触点和碰撞力,并对冲击过程中球形和凹形颗粒的缓冲性能进行离散元分析。数值结果表明,颗粒床厚度、冲击速度和颗粒形状显著影响球谐函数颗粒材料的缓冲性能。颗粒床底部的冲击力峰值随着颗粒床厚度和表面凹凸性的增加而降低,同时冲击力峰值随着冲击速度和颗粒球面度的增加而增加。与球形颗粒相比,球谐函数颗粒具有凹凸表面和多接触点特性,这有利于冲击荷载向四周扩展并提高凹形颗粒的缓冲效果。     

基于临界状态的砂土本构模型研究

砂土孔隙比及所受压力是其力学特性的重要影响因素. 本文基于砂土临界状态线特性分析,采用以e-(p/p_a)ξ平面内的线性关系描述其等向压缩线. 通过对比分析两种不同压缩线函数 与临界状态线函数之间的关系提出更适合描述砂土在等向压缩下的参考压缩线,并给出了基于参考压缩线的等向硬化规律. 建议了适用于 描述砂土剪切特性的屈服面函数,并给出利用等向压缩和等p路径确定屈服面形状参数μ的方法. 将不同应力比对应的压缩线作为砂土状态参量参考线,以获取潜在强度M_f与特征状态应力比M_c,进而描述砂土压缩与剪切特性;基于等向压缩与等p路径建立了当前应力比与状态参量参考线之间的相关关系,从而实现了砂土状 态参量参考线由参考压缩线向临界状态线平稳过渡. 建立的砂土本构模型共11个参数,均能够通过常规土工试验或经验获取. 基于模型预测与Toyoura砂的等向压缩、三轴不排水剪切试验及排水剪切试验的对比结果,本文建立的砂土本构模型很好地描述了Toyoura 砂在不同孔隙比和不同压力下的压缩与剪切特性.      

Mechanical behaviour of the upper layers of a sandy loam soil under shear loading

The mechanical behaviour of the upper layers of a sandy loam soil was studied under standard triaxial compression and direct shear box tests. Variations of soil material properties were investigated at four different initial dry bulk densities of 1410, 1520, 1610 and 1670 kg/m³. Soil deformation and volume change under the triaxial compression loading were also studied at these bulk densities. Results from the two tests showed increases in the soil mechanical properties with the initial dry bulk density. The internal friction angle values measured with the triaxial compression apparatus exceeded those measured with the direct shear box. In contrast, the soil cohesion values measured with the direct shear box exceeded those measured with the triaxial compression apparatus. Under the triaxial compression test, the loose soil samples underwent contraction and volume reduction, whereas the dense samples swelled and failure cracks appeared clearly at various planes. The soil contraction for the former case characterizes the occurrence of soil compaction, whereas the cracks propagation and volume increase in the latter case characterizes the breaking up and loosening of soil during tillage operations. For the loose and moderately compacted states, the engineering Poisson's ratio increased with the axial strain until loading was completed. It also increased at the compacted and very compacted states until reaching given loading stages, after which its value started to decrease. This shifting in the engineering Poisson's ratio during loading may provide another identification of the moment of soil failure occurrence, in addition to that of the maximum shear stress.     

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.     

Simulation and calibration of granules using the discrete element method

The discrete element method (DEM), developed by Cundall and Strack (1979) to solve geomechanical problems, is used to simulate the mechanical behavior of granules. According to the DEM, an individual granule can be modeled as a realistic mechanical system consisting of primary particles bonded by interaction forces.Granulometric properties of the model material, zeolite 4A, have been measured to determine their macro properties. To investigate the compression behavior, a compression test was performed using a strength tester on single granules between two pistons. A modeled granule consisting of more than 22,000 primary particles was generated. The micro properties of the modeled granule have been precisely set to allow its macro properties to be equivalent to the macro properties of zeolite 4A granules. To calibrate the mechanical properties, diametrical compression was simulated using two rigid walls stressed at a constant stressing velocity. The force–displacement curve of the modeled granule at compression has been calibrated by the experimental curve of zeolite 4A.     

Static and fatigue behavior of metal cylinders with enhanced surface property

Effects of surface property on the static and fatigue strength of metal cylinders is determined by analysis of the redistribution of stress and energy density. Because the mechanical constraints on the material elements in a surface layer differ so drastically from those in the interior or bulk, changes of mechanical properties in a thin surface layer can significantly alter the global behavior of a metal part. Desirable gain in mechanical strength due to surface treatments such as laser spraying, ion implantation, etc. is now well known. Developed in this work is a method for evaluating the enhancement of the static and fatigue life of metals when their mechanical properties in a thin surface layer of material are improved. An extensive increase in the surface layer thickness does not lead to appreciable changes of mechanical strength.     

Neighborhood Relationships of Widely Distributed and Irregularly Shaped Particles in Partially Dewatered Filter Cakes

A more thorough understanding of the properties of bulk material structures in solid–liquid separation processes is essential to understand better and optimize industrially established processes, such as cake filtration, whose process outcome is mainly dependent on the properties of the bulk material structure. Here, changes of bulk properties like porosity and permeability can originate from local variations in particle size, especially for non-spherical particles. In this study, we mix self-similar fractions of crushed, irregularly shaped Al₂O₃ particles (20 to 90 µm and 55 to 300 µm) to bimodal distributions. These mixtures vary in volume fraction of fines (0, 20, 30, 40, 50, 60 and 100 vol.%). The self-similarity of both systems serves the improved parameter correlation in the case of multimodal distributed particle systems. We use nondestructive 3D X-ray microscopy to capture the filter cake microstructure directly after mechanical dewatering, whereby we give particular attention to packing structure and particle–particle relationships (porosity, coordination number, particle size and corresponding hydraulic isolated liquid areas). Our results reveal widely varying distributions of local porosity and particle contact points. An average coordination number (here 5.84 to 6.04) is no longer a sufficient measure to describe the significant bulk porosity variation (in our case, 40 and 49%). Therefore, the explanation of the correlation is provided on a discrete particle level. While individual particles?<?90 µm had only two or three contacts, others?>?100 µm took up to 25. Due to this higher local coordination number, the liquid load of corresponding particles (liquid volume/particle volume) after mechanical dewatering increases from 0.48 to 1.47.

     

ICE自由面台阶靶数据处理的直接计算方法

斜波压缩台阶靶实验中,不同厚度界面粒子速度历史与材料压缩特性参数存在联系。然而利用普遍采用的数据处理方法无法直接获得该联系。本文中借助特征线理论在建立上述关联的基础上,实现未知EOS下斜波压缩流场的直接计算过程。经数值计算表明,该方法不仅在无强度效应数据处理中能够准确计算理论值,而且在含有强度效应数据处理中也能够较好地逼近理论值,并可在真实实验数据处理中获得与文献符合较好的结果。该研究可为探索具有较完整理论的强度效应数据处理方法提供新途径。     

颗粒增强复合材料压缩行为的位错动力学模拟

颗粒增强铜基复合材料因具有极高的强度和弹性模量, 优异的导电、导热性能和抗磨损能力, 被广泛应用于航天航空、汽车、电子工业等领域. 第二相强化是其主要的强化方式, 其通过合金中弥散的微粒阻碍位错运动, 可有效提高金属材料的力学性能, 提高其服役安全. 针对该问题本文采用三维离散位错动力学(three-dimensional discrete dislocation dynamics, 3D-DDD)方法, 对微尺度颗粒增强铜基复合材料进行了微柱压缩模拟, 分析了位错与第二相颗粒交互作用对材料力学响应的影响, 揭示第二相颗粒强化的微观机理. 本研究将第二相颗粒视为位错不可穿透的球形微粒, 采用位错绕过机制模拟颗粒与位错的交互作用过程. 通过调控滑移面相对于第二相颗粒中心的距离发现: 屈服应力和应变硬化率均随距离的增大而减小. 研究也发现Schmid因子越高的滑移系, 屈服应力越低, 后续应变硬化率越低. 多位错与颗粒交互作用的模拟发现, 同一滑移面中位错间的反应和不同滑移系中位错的交互作用可能是导致屈服应力和应变硬化率降低的关键.      

胶粒改性陶粒混凝土力学性能试验研究

采用两种改性方案配制胶粒陶粒混凝土,试验研究规定龄期下改性陶粒混凝土的抗折、抗压和劈裂抗拉强度,并分析不同的养护环境对它的力学性能影响.试验结果表明,对胶粒进行改性处理后,能显著改善胶粒与水泥基材料的粘结性能,提高胶粒陶粒混凝土的力学性能;在不同的养护条件下,胶粒陶粒混凝土力学性能的整体变化规律是一致的,养护条件对胶粒陶粒混凝土力学性能的影响不是很大.     

漂珠颗粒材料静动态力学性能与破碎机理研究

为考察脆性空心颗粒材料冲击载荷下的力学特性,以具有不同粒径分布的粉煤灰漂珠为研究对象,对其静动态力学性能进行实验研究。通过限制颗粒材料压缩应变为50%,分析颗粒破碎率和破碎机理与材料宏观应变率效应的关系。结果表明:（1）不同粒径的漂珠颗粒材料在动态压缩下较准静态压缩下颗粒材料的强度均有明显的增强,在0.001和150 s?1大小颗粒的强度分别提高200%和195%,在150和300 s?1大小颗粒的强度分别提高39%和51.5%,在300和800 s?1大小颗粒的强度并未发生明显的变化;（2）在相同加载速度下粒径较小的颗粒比大粒径颗粒的强度和吸能效果分别提高35%～40%和35%～48%;（3）对破碎后颗粒粒径分布曲线分析可知,随着加载速度的增加,大小颗粒的破碎率和破碎程度都会增大,且在相同加载速度下大颗粒的破碎率较小颗粒的破碎率高;（4）Hardin破碎势分析表明,单位输入能量下颗粒的相对破碎势随冲击速度增大而减小,动态冲击下用于颗粒破碎的能量利用率降低,从而导致材料在相同压缩量下产生更高的能量耗散和应力水平,即表现为宏观的应变率效应。     

一种预测颗粒增强复合材料界面力学性能的新方法

界面在颗粒增强复合材料中起到传递载荷的关键作用, 界面性能对复合材料整体力学行为产生重要影响. 然而由于复合材料内部结构较为复杂, 颗粒与基体间的界面强度和界面断裂韧性难以确定, 尤其是法向与切向界面强度的分别预测缺乏有效方法. 本文以氧化锆颗粒增强聚二甲基硅氧烷(PDMS)复合材料为研究对象, 提出一种预测颗粒增强复合材料界面力学性能的新方法. 首先, 实验获得纯PDMS基体材料及单颗粒填充PDMS试样的单轴拉伸应力$\!-\!$应变曲线, 标定出PDMS基体材料的单轴拉伸超弹性本构关系; 其次, 建立与单颗粒填充试样一致的有限元模型, 选择特定的黏结区模型描述界面力学行为, 通过样品不同阶段拉伸力学响应的实验与数值结果对比, 分别给出颗粒与基体界面的法向强度、切向强度及界面断裂韧性; 进一步应用标定的界面力学参数, 开展不同尺寸及不同数目颗粒填充试样的实验与数值结果比较, 验证界面性能预测结果的合理性. 本文提出的界面力学性能预测方法简便、易操作、精度高, 对定量预测颗粒增强复合材料的力学性能具有一定帮助, 亦对定量预测纤维增强复合材料的界面性能具有一定参考意义.