月壤双轴试验的剪切带离散元数值分析

土体的破坏问题一般都是从其剪切带入手进行研究的。针对真实月壤所处的环境（无水、低重力场、低气压等）和内摩擦角较大的特点,采用最近提出的1种考虑粒间抗转动作用与粒间范德华力2个因素的月壤散粒体力学接触模型,用离散单元法模拟了柔性边界条件下月球环境（含范德华力）与地球环境（不含范德华力）2种试样的双轴压缩试验,通过试样局部变形、速度场、孔隙比及转动场的变化情况研究了2种试样剪切带形成与发展。研究结果表明,剪切带的形成是试样内部应变局部化的结果,同时也伴随着试样内部颗粒相对转动的局部化,月面环境对试样的破坏形式与性状（剪切带的倾角及厚度）有显著的影响,探月工程中必须考虑到月面环境对月壤力学性质的影响。     

抗转动对颗粒材料组构特性的影响研究

颗粒材料大多由不规则形状的颗粒组成,如砂土、谷物等,抵抗转动是不规则形状颗粒的固有特性。已有研究表明,颗粒抗转动特性对其宏观力学特性有显著影响。因此,在颗粒材料的细观数值模拟中或采用非圆颗粒,或在圆颗粒离散元模拟中采用考虑抗转动的接触模型。采用有限元-离散元耦合方法（FDEM）和离散元方法（DEM）分别对椭球形状颗粒和具有抗转动能力的圆球颗粒进行三轴剪切数值模拟,指出了采用抗转动接触模型考虑颗粒形状影响的局限性,并基于颗粒的局部排布结构揭示了形状影响的细观来源。峰值内摩擦角和剪胀均随着转动摩擦系数和形状偏离圆球程度而单调增加,但颗粒形状对它们的影响呈现出明显的收敛趋势。细观组构分析也表明,虽然颗粒形状和转动摩擦都能显著增强组构各向异性,但是组构各向异性的演化模式有明显的区别。造成以上结果的差异在于其抵抗转动的影响机制不同。转动摩擦是通过限制颗粒转动,增强了颗粒间的稳定承载能力,而非圆颗粒是通过咬合作用形成稳定的局部排列结构。由于椭球颗粒腹部比端部能够传递更大的接触力,颗粒受剪切后发生转动,颗粒长轴倾向于正交大主应力方向,呈现交错排列,颗粒间相互锁定。     

月壤中的氦-3

除了极少数非常陡峭的山脉与撞击坑和火山通道的峭壁外，整个月球表面几乎都被一层厚度不等的月尘、岩屑和岩块的混合堆积物即月壤所覆盖。由于月球无大气层等特殊环境，太阳光长驱直入，太阳风粒子直接注入到月壤细小颗粒上使月壤中富含稀有气体等太阳风粒子组分。本文在系统阐述月壤的形成过程与形成机制的基础上，分析了月壤中稀有气体的来源及其浓度与月壤的成熟度、月壤颗粒大小、月壤矿物组成和化学成分的相关关系，进而利用已有的探测数据和分析结果，对月壤中氦-3资源的开发利用前景进行了初步评估。     

用于月面车辆力学试验的模拟月壤研究

研究月面环境下的车辆地面力学,对保证月球探测车辆的正常工作具有重要意义,而用于月面车辆力学试验的模拟月壤则影响其结果的准确度与可信度。用于月面车辆力学试验的模拟月壤,以月壤样本与JSC-1模拟月壤为参考标准,以吉林辉南县火山灰为主要原料、赤铁矿砂为辅料,通过调整试样的粒径分布与赤铁矿砂含量,使模拟月壤的比重、内摩擦角、凝聚力、承压特性、粒径分布、颗粒形态、矿物组成等参数接近JSC-1模拟月壤,并在月壤样本的参数变化范围之内,表明该模拟月壤可用于月面车辆力学试验研究。     

岩石单颗粒压缩破碎的数值模拟分析

岩石颗粒破碎是影响粒状材料剪切强度和变形的最主要因素, 岩石颗粒破碎并不是想象的那么难, 像花岗岩颗粒有时在很小的压应力作用下就可以破碎。岩石单颗粒破碎的物理试验结果常常很离散, 完成大量单颗粒破碎的物理试验费时费力不现实, 采用离散单元法(Discrete element method, DEM)PFC软件模拟单颗粒压缩破碎试验, 既能克服单颗粒破碎物理试验的缺陷, 又能解决单颗粒破碎物理试验工作量大的难题, 是研究单颗粒破碎的理想选择。基于DEM的软件PFC2D, 将粒径为0.075～0.1245mm的基本粒子捆绑成不同粒径的单颗粒, 模拟岩石单颗粒压缩破碎试验, 观察颗粒破碎演化过程, 统计单颗粒破碎强度。计算单颗粒压缩破碎后颗粒分布的分维, 验证单颗粒破碎强度的分形模型和单颗粒破碎强度的尺寸效应。文中引用玄武岩单颗粒破碎试验结果, 与单颗粒破碎的离散单元模拟结果进行比较, 验证单颗粒破碎强度的尺寸效应和修正的Weibull理论的离散单元模拟结果。     

Numerical study of inter‐particle bond failure by 3D discrete element method

The cohesive‐frictional nature of cementitious geomaterials raises great interest in the discrete element method (DEM) simulation of their mechanical behavior, where a proper bond failure criterion is usually required. In this paper, the failure of bond material between two spheres was investigated numerically using DEM that can easily reproduce the failure process of brittle material. In the DEM simulations, a bonded‐grain system (composed of two particles and bond material in between) was discretized as a cylindrical assembly of very fine particles connecting two large end spheres. Then, the bonded‐grain system was subjected to compression/tension, shear, rolling and torsion loadings and their combinations until overall failure (peak state) was reached. Bonded‐grain systems with various sizes were employed to investigate bond geometry effects. The numerical results show that the compression strength is highly affected by bond geometry, with the tensile strength being dependent to a lesser degree. The shear, rolling and torsion strengths are all normal force dependent; i.e., with an increase in the normal force, these strengths first increase at a declining rate and then start to decrease upon the normal force exceeding a critical value. The combined actions of shear force, rolling moment and torque lead to a spherical failure envelope in a normalized loading space. The fitted bond geometry factors and bond failure envelopes obtained numerically in this three‐dimensional study are qualitatively consistent with those in previous two‐dimensional experiments. The obtained bond failure criterion can be incorporated into a future bond contact model. Copyright © 2015 John Wiley & Sons, Ltd.     

Effects of grain morphology on critical state: a computational analysis

We introduce a new DEM scheme (LS-DEM) that takes advantage of level sets to enable the inclusion of real grain shapes into a classical discrete element method. Then, LS-DEM is validated and calibrated with respect to real experimental results. Finally, we exploit part of LS-DEM potentiality by using it to study the dependency of critical state (CS) parameters such as critical state line (CSL) slope \(\lambda \), CSL intercept \(\varGamma \), and CS friction angle \(\varPhi _{\mathrm{CS}}\) on the grain’s morphology, i.e., sphericity, roundness, and regularity. This study is carried out in three steps. First, LS-DEM is used to capture and simulate the shape of five different two-dimensional cross sections of real grains, which have been previously classified according to the aforementioned morphological features. Second, the same LS-DEM simulations are carried out for idealized/simplified grains, which are morphologically equivalent to their real counterparts. Third, the results of real and idealized grains are compared, so the effect of “imperfections” on real particles is isolated. Finally, trends for the CS parameters (CSP) dependency on sphericity, roundness, and regularity are obtained as well as analyzed. The main observations and remarks connecting particle’s morphology, particle’s idealization, and CSP are summarized in a table that is attempted to help in keeping a general picture of the analysis, results, and corresponding implications.     

Mechanical characteristics of a lunar regolith simulant at low confining pressure

Manned lunar exploration has recently attracted renewed interest. This includes the NASA Constellation program to return humans to the Moon by 2020, the ESA Aurora program which may use the Moon as a way station to prepare for major interplanetary exploration by 2025, and the PRC program to send a human to the Moon by 2030 and build a temporary manned lunar base by 2040. One of the problems demanding a solution is the stresses on the mechanical characteristics of the lunar regolith under the microgravity environment. The gravity on the Moon is about 1/6 that on Earth. The regolith is subject to very low confining stresses under a microgravity environment and the mechanical properties can change correspondingly. Because of the limited amount of lunar regolith brought back to Earth by the Apollo missions, a lunar regolith simulant was developed using silicon carbide to investigate the properties of the lunar regolith. Based on triaxial tests, this study analyzed the mechanical properties of the lunar regolith simulant at low stresses including the shear strength, peak strength and dilatation angle. The research results provide useful information on lunar regolith characteristics for astronauts returning to the Moon and for building a temporary manned lunar base.     

基于FDTD的月壤分层雷达探测正演模拟

“嫦娥三号”巡视器上搭载的测月雷达通过500 MHz频率的天线,对巡视路线上的月壳浅层结构进行剖面式观测,能探测到月球地底下30 m深土壤层的结构。为了更好地分析测月雷达的探测结果,文中对月壤分层进行正演模拟,月岩层的介电常数采用与“嫦娥三号”登月点位置最近的Apollo 15的电性参数近似代替,月壤的介电常数是一个和密度相关的值且随深度的变化而变化。根据以上参数建立模型,基于时域有限差分（FDTD）原理,利用GprMax对月壤分层结构进行模拟,建立的模型包括月表介电常数随深度变化二层模型和月岩层存在下界面三层模型,并对其波形特征进行分析,找到月壤层次划分的一般规律,对“嫦娥三号”测月雷达数据的分析提供理论依据。     

堆石料单粒强度尺寸效应的颗粒流模拟方法研究

围绕堆石料单粒强度尺寸效应的颗粒流模拟方法展开研究。首先,基于FISH二次开发建立了堆石料的随机不规则单粒模型,充分考虑堆石料的形状特征和破碎现象;然后,建立了堆石料单粒强度尺寸效应的等效模拟方法,以单粒强度随其粒径的变化规律为基础,推导了堆石料模型中细观黏结强度与堆石料等效粒径的负指数经验公式;其次,基于建立的数值模型对堆石料的室内单粒压缩试验进行仿真模拟,验证数值模型的正确性和合理性,并对较大粒径堆石料的单粒强度进行模拟预测,突出数值试验的优势;最后,基于建立的数值模型对相同粒径不同形状特征堆石料的单粒强度分布特征进行模拟研究。研究结果表明：（1）堆石料内部缺陷含量和尺寸随粒径增加对其单粒强度所产生的尺寸效应,可通过堆石料模型中细观强度参数随粒径折减进行等效模拟;（2）形状特征对堆石料的破裂机制具有重要影响,方形颗粒为压剪破裂,单粒强度较高,而随机不规则颗粒和圆形颗粒为拉剪或劈裂,单粒强度相对较低;（3）拉剪或劈裂条件下,堆石料形状越不规则,其单粒强度的离散程度越高,反之则离散程度越低。相关研究成果可为进一步研究荷载作用下堆石体内各粒径段堆石料的破碎量奠定基础,从而更加真实地反映堆石体的级配演化规律。     

A micromechanical study of rolling and sliding contacts in assemblies of oval granules

The evolution of the microstructure of an assembly of cohesionless granular materials with associated pores, which carry the overall applied stresses through frictional contacts is a complex phenomenon. The macroscopic flow of such materials take place by the virtue of the relative rolling and sliding of the grains on the micro‐scale. A new discrete element method for biaxial compression simulations of random assemblies of oval particles with mixed sizes is introduced. During the course of deformation, the new positions of the grains are determined by employing the static equilibrium equations. A key aspect of the method is that, it is formulated for ellipse cross‐sectional particles, hence desirable inherent anisotropies are possible. A robust algorithm for the determination of the contact points between neighbouring grains is given. Employing the present methodology, many aspects of the behaviour of two‐dimensional assemblies of oval cross‐sectional rods have been successfully addressed. The effects of initial void ratio, interparticle friction angle, aspect ratio, and bedding angle on the rolling and sliding contacts are examined. The distribution of normals to the rolling and sliding contacts have different patterns and are concentrated along directions, which are approximately perpendicular to one another. On the other hand, the distribution of all contact normals (combined rolling and sliding) are close to that of rolling contacts, which confirm that rolling is the dominant mechanism. This phenomenon becomes more pronounced for higher intergranular friction angle. Characteristics of the rolling and sliding contacts are also discussed in the context of the force angle, which is the inclination of contact force with respect to the contact normal. Copyright © 2003 John Wiley & Sons, Ltd.     

On the stress-force-fabric equation in triaxial compressions: Some insights into the triaxial strength

The strength of granular materials during triaxial compression is investigated via a grain scale analysis in this paper. A 3D Discrete Element Method (DEM) program provides the triaxial strength data and helps to validate the micromechanical analysis. Some standard methods in statistics are employed first to quantitatively examine the assumptions made when deriving the stress-force-fabric (SFF) equation. After careful validation, a more concise format for the SFF equation is proposed for triaxial compressions. With this SFF equation, the strength is found to be jointly contributed by the magnitudes of the contact force anisotropy and fabric anisotropy. The influence of the initial void ratio, confining pressure and loading direction on the development of contact force anisotropy and fabric anisotropy is examined and presented. With similar techniques, the “force” term in the SFF equation is further decoupled, and an equation is obtained such that it explicitly links the contact force term with the friction coefficient between grains, a tensor defined as a statistic of the normal contact forces and a tensor defined as a statistic of the mobilisation status of contacts. Based on this equation, another equation regarding the stress ratio of granular assembly is obtained, and it clearly indicates two sources that contribute to the phenomenological friction nature of granular assembly. These two sources are caused by the contact force at the grain scale. The first is the anisotropy of the average normal contact forces, and the second is the mobilisation of contacts.     

月壤的物理和机械性质

月壤是在O2、水、风和生命活动都不存在的情况下，由陨石和微陨石撞击、宇宙射线和太阳风轰击、月表温差导致岩石热胀冷缩破碎等因素的共同作用下形成的。月壤独特的形成过程，加上独特的月表环境，使月壤在粒度分布、颗粒形态、颗粒比重、孔隙比和孔隙率、电性和电磁性质、压缩性、抗剪性、承载力等方面均与地球土壤存在较大差异，这些参数的平均值和最佳估计值，可以作为月表机械设计和操作、宇航员装备设计、月球着陆场选址的主要依据，对月球资源开发和利用以及月球基地建设具有极其重要的意义。     

简化接触模型的月壤离散元数值分析

根据月壤其颗粒级配可归类于粉质砂土。针对真实月壤所处的环境（无水、低重力场、低气压等）,将Perko等2001年提出的月壤颗粒间的范德华力植入离散元分析软件PFC2D中,模拟了刚性边界下加入该模型与未加该模型试样的双轴压缩试验,研究了颗粒间范德华力对试样的宏观力学特性与微观颗粒接触的影响。结果表明,颗粒间的范德华力对试样的抗剪强度、体应变以及颗粒平均配位数都有显著的影响     

Noble Gases in the Lunar Regolith

The most fundamental character of lunar soil is its high concentrations of solar-wind-implanted elements, and the concentrations and behavior of the noble gases He, Ne, Ar, and Xe, which provide unique and extensive information about a broad range of fundamental problems. In this paper, the authors studied the forming mechanism of lunar regolith, and proposed that most of the noble gases in lunar regolith come from the solar wind. Meteoroid bombardment controls the maturity of lunar soil, with the degree of maturation decreasing with grain size; the concentrations of the noble gases would be of slight variation with the depth of lunar soil but tend to decrease with grain size. In addition, the concentrations of noble gases in lunar soil also show a close relationship with its mineral and chemical compositions. The utilization prospects of the noble gas ^3He in lunar regolith will be further discussed.     

CAS-1模拟月壤

模拟月壤是与月球月壤具有相似的矿物组成、化学成分和物理力学性质的地球物质，是月球样品的地球化学复制品。长白山龙岗火山群金龙顶子火山喷发的四海火山渣具有与阿波罗14号采集的月球样品相似的化学和矿物组成，并含有20％～40％的玻璃物质。以四海火山渣为初始物质，研制成功CAS-1模拟月壤，并测量了CAS-1模拟月壤的主量和微量元素组成、矿物组成、密度、颗粒形态、粒度分布、抗剪性和复介电常数等参数。结果表明，CAS-1模拟月壤与Apollo 14号采集的月球样品具有相似的化学成分、矿物组成和物理力学性质，是一种理想的低钛玄武岩质模拟月壤。     

月球探测中月面热环境影响的研究现状

在月球探测中,多光谱、热红外、被动微波辐射等探测技术被广泛应用于月表物质组分和物理特性的探测,也积累了大量探测数据。月球太阳辐照、月球表面温度、地球反照和内部热流等月面热环境的变化,改变了月表物质反射率、热发射率以及其他电磁学等基本性质,制约了探测数据的准确解译;同时,大幅度的太阳辐射强度和月面温度变化也直接威胁月面探测中巡视探测器和宇航员的安全。但是,目前月球表面热环境对探测活动的影响认识还比较欠缺,月面热环境的时空变化规律认识还不够充分,在实验研究不足的情况下对各种探测方式的影响缺乏系统的理解。结合月球探测的发展,进一步立足实验手段和探测结果,通过开展不同地形条件下月面太阳辐射和物理温度的时间变化和空间分布规律研究、探测数据与月面热环境参数时空匹配问题研究、建设具备类似物质组成和月面热环境特征的实验场地以及开展系统的热环境影响模拟实验研究,全面认识月面热环境参数的时空变化规律,探讨月面热环境对不同探测方式的影响将是月面热环境研究的重要内容。     

控制应力路径散体材料真三轴试验强度特征的离散元模拟与细观机制分析

通过一系列真三轴离散元数值试验,模拟了不同应力路径下的等b试验中散体材料的强度特征。根据模拟结果详细地分析了三维应力条件下中主应力和应力路径对散体材料峰值强度的影响,研究了峰值摩擦角、峰值应力比的变化规律,并根据真应力的概念和组构张量的演化结果分析了散体材料的强度成因。研究表明,在不同类型的数值试验中峰值偏应力随b参数的变化规律不同,但采用初始围压归一化后的应力-应变曲线规律一致。峰值强度线的斜率只与b值有关而与应力路径无关,且随着b值的增加,峰值应力比qf /pf逐渐减小,数值模拟结果与室内试验结果吻合较好;随着应变的发展,数值试样的组构也随之发生变化,产生了明显的应力诱发各向异性;散体的强度为颗粒摩擦及材料各向异性共同作用的结果;理论上,组构比-应力比坐标系中破坏点位置仅取决于颗粒摩擦角 ,而数值模拟结果与理论值的差异源于颗粒间咬合和滚动摩擦的影响,其影响与颗粒表面摩擦系数有关,也受空间应力状态的影响。