颗粒介质弹性的弛豫

颗粒介质是复杂的多体相互作用体系, 其弹性源自内部的力链结构, 弹性能量处在亚稳态, 具有复杂的弛豫行为. 在常规作用下, 颗粒介质往往呈现明显的弹性弛豫. 应力松弛是应变恒定时应力的衰减现象, 弹性弛豫是应力松弛的主要原因. 在前期工作基础上, 从弹性势能面和双颗粒温度热力学角度分析了弹性弛豫的机理, 量化了弹性应力演化不可逆过程; 基于双颗粒温度热力学计算得到了弹性能、颗粒温度和应力的演化, 其中应力松弛的计算结果与实验结果基本一致, 讨论了颗粒温度初值和输运系数的影响. 指出, 开展力链结构及其动力学研究是揭示宏观弹性弛豫机理的关键.     

随机反应扩散格气模型上钙动力学的粗粒化模拟

发展了一种基于随机格气模型的粗粒化方法,该方法能有效模拟内质网表面钙动力学信息. 首先将相邻的微观节点合并成粗粒化节点,再根据局域平均场近似推导出粗粒化反应速率,然后执行粗粒化动力学蒙特卡洛模拟. 发现粗粒化动力学蒙特卡洛模拟结果和微观模拟结果非常吻合. 有趣的是,存在一个最佳的粗粒化比m,使得粗粒化模拟与微观模拟的相变点偏差最小. 固定m,发现临界点随体系尺度增加而单调增加,而且相变点的偏差与体系尺度存在一个标度关系.此外,该粗粒化方法大大地加快了蒙特卡洛模拟速率,并且与微观模拟直接相关. 该方法可以广泛用来研究体系尺度效应,而节省大量计算时间.     

基于微观结构计算晶体宏观弹性应变场

基于晶体的微观结构提出了一种计算晶体宏观弹性应变场的方法,将该方法应用于金属晶体Fe,得到的应变场与连续介质理论预测结果吻合地很好,从而证明了该方法的准确性.以该方法为基础,进一步考虑无限小应变与有限应变的差值,发现应变差值在晶体剧烈变形区域内较强,对金属晶体Cu位错应变场的计算表明这种应变差异可以用来描述位错的运动状态和芯域结构.该方法可为评估晶体缺陷附近应变状态提供依据.     

基于微观结构计算晶体宏观弹性应变场

基于晶体的微观结构提出了一种计算晶体宏观弹性应变场的方法,将该方法应用于金属晶体Fe,得到的应变场与连续介质理论预测结果吻合地很好,从而证明了该方法的准确性.以该方法为基础,进一步考虑无限小应变与有限应变的差值,发现应变差值在晶体剧烈变形区域内较强,对金属晶体Cu位错应变场的计算表明这种应变差异可以用来描述位错的运动状态和芯域结构.该方法可为评估晶体缺陷附近应变状态提供依据.     

复杂网络谱粗粒化方法的改进算法

大规模网络的同步问题是网络科学的重要研究课题之一.粗粒化方法提供了一种将大规模网络转化为小规模网络,同时又能较好地保持原始网络的拓扑性质或动态特性的研究途径,其中比较有代表性的谱粗粒化方法能较好地保持初始网络的同步能力.然而,谱粗粒化方法在实际计算中计算量大、对实际大规模网络可执行性差.本文提出一种改进的谱粗粒化算法,能大幅减少计算量,同时获得更好的谱粗粒化效果.通过理论分析和大量的数值仿真实验验证了所提改进算法的粗粒化效果和计算量都明显优于原谱粗粒化方法.     

不同生长环境下砷化镓纳米颗粒的应变场模拟

对于埋嵌在薄膜材料中的纳米颗粒,在其生长过程中总是不可避免地伴随着应变场的产生,而这种应变场的分布能反映纳米颗粒的结构变化,纳米颗粒结构与它的物理性能有重要的关系.研究埋嵌在不同薄膜材料中的纳米颗粒生长过程中的应变场分布对于调控纳米颗粒的物理性能有着重要的意义.本文利用有限元算法分别计算了埋嵌在非晶氧化铝薄膜和非晶二氧化硅薄膜材料中的砷化镓纳米颗粒生长过程中的应变场分布.砷化镓纳米颗粒在以上两薄膜材料生长过程中都受到非均匀偏应变作用.对于埋嵌在氧化铝薄膜中的砷化镓纳米颗粒,其生长过程中,纳米颗粒内部受到的应变大于纳米颗粒表面受到的应变;而对于埋嵌在二氧化硅薄膜中的砷化镓纳米颗粒,纳米颗粒内部受到的应变小于纳米颗粒表面受到的应变.选择砷化镓纳米颗粒生长的薄膜材料可以调控纳米颗粒生长过程中的应变场分布,从而进一步调控纳米颗粒的晶格结构和形貌及其物理性能.     

惯性颗粒所见被动标量统计特性的直接数值模拟研究

本文采用直接数值模拟方法,对惯性颗粒所见均匀各向同性湍流中具有平均标量梯度的被动标量场统计特性进行了研究。模拟结果表明:与颗粒所见流体速度的自相关特性不同,颗粒所见标量的自相关特性随颗粒惯性的增加而单调减少;颗粒所见标量脉动能随颗粒惯性的增加先减少再增大,在St≈1.0的临界颗粒附近达到最小值,而颗粒所见标量耗散率随颗粒惯性的变化行为与颗粒所见标量脉动能的变化行为相反;数值模拟的结果进一步揭示,在St≈1.0时颗粒所见流体标量脉动能和耗散率的极值是因为St≈1.0的临界颗粒聚集于低涡量、高应变区域和标量场在高应变区域形成强耗散的片状结构所致。     

颗粒介质的结构及热力学

颗粒介质具有远程无序和近程有序的结构, 是产生动力学不均匀性(dynamical heterogeneity) 和复杂不可逆过程的根源. 本文分析了颗粒介质的结构特征、变形和能量耗散之间的内在关联, 讨论了颗粒介质的弹性, 提出了流变应变增量、耦合应变增量和弹性应变增量的应变增量分解方式. 沿用非平衡热力学框架, 引入表征运动无序的动理学颗粒温度T^k和表征弹性应力涨落的构型温度T^c, 作为非平衡态变量, 建立了双颗粒温度热力学(two-granular-temperature thermodynamics, TGT理论), 注重分析了不可逆过程中的热力学力和流, 并与著名的砂土内变量热力学进行了对比.     

相场法模拟弹性场对沉淀相变组织演化及相平衡成分的影响

采用相场法模拟沉淀组织在不均匀共格本征弹性场及外加约束弹性应变作用下的演化历程，探讨了不均匀共格本征弹性场及外加约束弹性应变对共格沉淀组织演化的影响规律，并分析了弹性场在沉淀相变中对基体与沉淀相的平衡成分的影响. 关键词： 相场法 沉淀组织 弹性场 相平衡成分     

基于链化分析的磁流变液剪切屈服应力模型

基于磁性颗粒在磁场作用下的链化分析和统计分析方法,建立了磁流变液的宏观屈服剪应力的分析模型.模型考虑了磁感应强度、颗粒尺寸、体积分数、剪应变率以及饱和磁化强度等因素,能描述不同剪切应变率下剪切应力的变化及在高应变率下可能出现的剪切稀化效应.分析了不同因素对剪切屈服应力的影响,讨论了提高磁流变液剪切屈服应力的途径.分析表明,该模型的计算结果能较好地描述有关实验现象,并可用于高性能磁流变液的设计分析.     

Softening induced instability of a stretched cohesive granular layer

We report on a cellular pattern which spontaneously forms at the surface of a thin layer of a cohesive granular material submitted to in-plane stretching. We present a simple model in which the mechanism responsible of the instability is the "strain softening" exhibited by humid granular materials above a typical strain. Our analysis indicates that such a type of instability should be observed in any system presenting a negative stress sensitivity to strain perturbations.     

Discrete element simulation of localized deformation in stochastic distributed granular materials

The deformation in granular material under loading conditions is a problem of great interest currently. In this paper, the micro-mechanism of the localized deformations in stochastically distributed granular materials is investigated based on the modified distinct element method under the plane strain conditions, and the influences of the confining pressure, the initial void ratio and the friction coefficient on the localized deformation and the stability of granular materials are also studied. It is concluded, based on the numerical simulation testing, that two crossed shear sliding planes may occur inside the granular assembly, and deformation patterns vary with the increasing of transverse strain. These conclusions are in good agreement with the present experimental results. By tangential velocity profiles along the direction normal to the two shear sliding planes, it can be found that there are two different shear deformation patterns: one is the fluid-like shear mode and the other is the solid-like shear mode. At last, the influences of various material parameters or factors on localized deformation features and patterns of granular materials are discussed in detail. Supported by the Key Project of the National Natural Science Foundation of China (Grant No. 10532040)     

Rheology and contact lifetimes in dense granular flows

We study the rheology and distribution of interparticle contact lifetimes for gravity-driven, dense granular flows of noncohesive particles down an inclined plane using large-scale, three dimensional, granular dynamics simulations. Rather than observing a large number of long-lived contacts as might be expected for dense flows, brief binary collisions predominate. In the hard-particle limit, the rheology conforms to Bagnold scaling, where the shear stress is quadratic in the strain rate. As the particles are made softer, however, we find significant deviations from Bagnold rheology; the material flows more like a viscous fluid. We attribute this change in the collective rheology of the material to subtle changes in the contact lifetime distribution involving the increasing lifetime and number of the long-lived contacts in the softer particle systems.     

Protocol dependence of mechanical properties in granular systems

We study the protocol dependence of the mechanical properties of granular media by means of computer simulations. We control a protocol of realizing disk packings in a systematic manner. In 2D, by keeping material properties of the constituents identical, we carry out compaction with various strain rates. The disk packings exhibit the strain rate dependence of the critical packing fraction above which the pressure becomes non-zero. The observed behavior contrasts with the well-studied jamming transitions for frictionless disk packings. We also observe that the elastic moduli of the disk packings depend on the strain rate logarithmically. Our results suggest that there exists a time-dependent state variable to describe macroscopic material properties of disk packings, which depend on its protocol.     

Compactivity and transmission of stress in granular materials

We outline a statistical-mechanical theory of granular materials. Stress propagation and force fluctuations in static granular media are still poorly understood. We develop the statistical-mechanical theory that delivers the fundamental equations of stress equilibrium. The formalism is based on the assumptions that grains are rigid, cohesionless, and that friction is perfect. Since grains are assumed perfectly rigid, no strain or displacement field can enter the equations for static equilibrium of the stress field. The complete system of equations for the stress tensor is derived from the equations of intergranular force and torque balance, given the geometric specification of the material. These new constitutive equations are indeed fundamental and are based on relations between various components of the stress tensor within the material, and depend on the topology of the granular packing. The problem of incorporating into the formalism the "no tensile forces" constraint is considered. The compactivity concept is reviewed. We discuss the relation between the concept of compactivity and the problem of stress transmission. (c) 1999 American Institute of Physics.     

基于$lt;i$gt;J$lt;/i$gt;积分的颗粒煤岩单轴压缩下裂纹扩展研究

颗粒煤岩是由众多离散的煤岩颗粒组成的固态多层次多结构物质,具有煤岩与颗粒物质的双重性质,其裂纹扩展规律可以从煤岩力学特性和颗粒物质多尺度特性进行研究. 首先,从能量角度对线弹性材料受压破坏,裂纹扩展产生原因进行了阐述,指出线弹性阶段裂纹的扩展动力源自应变能的释放. 然后,通过物理实验和数值试验从宏观和细观两方面对颗粒煤岩受压破裂过程中裂纹扩展做了进一步研究,结果表明:颗粒煤岩完全破裂后,底部会形成一个锥形堆,裂纹的扩展随着煤岩颗粒粒径的减小而减缓,部分裂纹扩展会出现突变点,且裂纹无光滑性;由于煤岩颗粒粒径等引起介质的非均匀性对裂纹扩展有重要的影响,均质度系数越大裂纹起裂时间越晚,声发射能量释放在裂纹扩展的轻度、中度和深度三个不同阶段逐渐变得频繁、剧烈. 研究结果将有利于进一步研究岩土类颗粒材料受压破裂过程的裂纹扩展规律. 关键词： J积分')" href="#">J积分 颗粒煤岩 单轴压缩 裂纹扩展     

Strain stiffening in random packings of entangled granular chains

Random packings of granular chains are presented as a model system to investigate the contribution of entanglements to strain stiffening. The chain packings are sheared in uniaxial compression experiments. For short chain lengths, these packings yield when the shear stress exceeds the scale of the confining pressure, similar to granular packings of unconnected particles. In contrast, packings of chains which are long enough to form loops exhibit strain stiffening, in which the effective stiffness of the material increases with strain, similar to many polymer materials. The latter packings can sustain stresses orders-of-magnitude greater than the confining pressure, and do not yield until the chain links break. X-ray tomography measurements reveal that the strain-stiffening packings contain system-spanning clusters of entangled chains.     

Refraction of shear zones in granular materials

We study strain localization in slow shear flow focusing on layered granular materials. A heretofore unknown effect is presented here. We show that shear zones are refracted at material interfaces in analogy with refraction of light beams in optics. This phenomenon can be obtained as a consequence of a recent variational model of shear zones. The predictions of the model are tested and confirmed by 3D discrete element simulations. We found that shear zones follow Snell's law of light refraction.     

Granular flows in a rotating drum: the scaling law between velocity and thickness of the flow

The flow of dry granular material in a half-filled rotating drum is studied. The thickness of the flowing zone is measured for several rotation speeds, drum sizes and beads sizes (size ratio between drum and beads ranging from 47 to 7400). Varying the rotation speed, a scaling law linking mean velocity vs. thickness of the flow, v∼h^m, is deduced for each couple (beads, drum). The obtained exponent m is not always equal to 1, the value previously reported for a drum in litterature, but varies with the geometry of the system. For small size ratios, exponents higher than 1 are obtained due to a saturation of the flowing zone thickness. The exponent of the power law decreases with the size ratio, leading to exponents lower than 1 for high size ratios. These exponents imply that the velocity gradient of a dry granular flow in a rotating drum is not constant. More fundamentally, these results show that the flow of a granular material in a rotating drum is very sensible to the geometry, and that the deduction of the “rheology” of a granular medium flowing in such a geometry is not obvious.     

Numerical method for constructing a continual deformation model equivalent to a specified discrete element model

A class of boundary value problems of solid mechanics and mechanics of granular media is identified in which uniform strain and strain rate distributions are possible. Boundary conditions appropriate to these problems are applied for statements of numerical experiments. An affine transformation method is proposed for constructing a continual model equivalent to a specified discrete element model. An example of using the method to construct a three-dimensional model of granular media is given.