带有点缺陷的二维颗粒系统离散元模拟

首先用离散元方法研究了颗粒系统在各向同性挤压和纯剪切状态下粒子间力的分布情况,并与相同条件下的实验结果进行对比. 然后模拟了带有不同数目点缺陷的二维颗粒系统在各向同性挤压和纯剪切时粒子间力的分布情况,并与无缺陷的情况做了比较,发现了点缺陷对颗粒系统的影响规律. 关键词： 颗粒系统 离散元 点缺陷 力概率分布     

光弹性颗粒的线扫描测力方法

提出了基于对光弹颗粒的应力光图进行加载点至圆心的线扫描以测接触力大小的办法.在正交圆偏振光场中获取了光弹颗粒的应力光图的等差条纹图案后,通过数字图像的处理,获取颗粒模型边界,得到颗粒圆心,半径等信息.并从应力光图的平方灰度梯度分布图出发,提出一种能够自动检测颗粒的接触点位置信息的方法.由多点加载的二维圆盘颗粒的光弹性原理,结合接触点位置信息,在力范围内,模拟得到每一种的接触力情况的每一加载点至圆心连线上光强灰度变化.将其与实验光强结果对比,可反推出颗粒接触力的大小.     

玻璃-橡胶混合颗粒体系的弹性行为研究

废旧橡胶制品颗粒与砂土颗粒混合物作为建筑填充材料具有环保、轻质、减震效果好等特点.软硬组分的混合比例可以调制体系力学性能从而实现兼顾材料柔韧性与强度的需求,但细观层面上材料性能改变的原因尚不明确.本文主要研究玻璃-橡胶混合颗粒体系的弹性行为及其微观机制.利用飞行时间法测量混合材料等效动弹性模量,发现随着橡胶颗粒增加,体系逐渐从类玻璃刚性行为转变为类橡胶柔性行为.离散元模拟结果与实验结果类似.此外,模拟显示低橡胶颗粒占比样品内主要由玻璃颗粒构成主力链结构,而橡胶颗粒基本不参与强力链的构成.当橡胶颗粒占比较大时,玻璃颗粒和橡胶颗粒共同构成主力链网络结构,但颗粒间法向接触力分布相对更为均匀,可视为玻璃颗粒悬浮于橡胶颗粒中.基于上述结果,提出了改进的等效介质理论,用于描述混合颗粒体系的弹性行为.研究认为:橡胶颗粒占比较小时内部颗粒的变形相对均匀,材料近似满足等应变假设,视为并联弹簧模型;橡胶颗粒占比较大时混合材料近似满足等应力假设,视为串联弹簧模型.两种模型得到的结果与模拟结果一致.上述结果有利于从微观角度揭示混合颗粒材料弹性行为的变化机制.     

火灾烟颗粒分形模型和球形模型光散射的比较研究

对烟颗粒的光散射进行模拟计算是研究火灾烟颗粒光散射特性的重要手段,目前对于火灾烟颗粒光散射的数值计算多采用球形或椭球模型.实际上,火灾烟颗粒的形貌与球形和椭球均存在着显著差异.扫描电子显微镜图像表明,烟颗粒具有近似分形的结构.本文利用离散偶极近似方法计算了随机取向的火灾烟颗粒分形凝团以及同体积的球形颗粒的光散射Muller矩阵,并对两者的归一化Muller矩阵元素随散射角的分布进行了比较.研究表明:火灾烟颗粒分形模型和球形模型的归一化矩阵元素F₁₁(θ)/     

同材质颗粒不同填充密度的随机介质中光场的空间分布

用时域有限差分法研究了同一材质的颗粒在不同填充密度下的随机介质中光局域化问题. 依据随机介质激光的实验参数,模拟了颗粒填充密度不同的随机介质中光场的空间分布. 结果表明：当散射颗粒的散射平均自由程与波长相当时,随机介质中的光场分布呈现局域化的特征,而且随着颗粒填充密度增大,光场的局域化程度增强. 因此,在同样的抽运激励下,颗粒填充密度越大的随机增益介质越容易产生激光辐射. 数值模拟结果与实验定性符合. 关键词： 随机激光 时域有限差分法 颗粒填充密度 局域化     

颗粒尺寸分散度对颗粒体系力学和几何结构特性的影响

利用颗粒离散元方法, 研究了由2048个光滑颗粒组成的体系在各向同性压缩条件下, 颗粒尺寸分散度s对颗粒体系力学和几何结构特性的影响. 结果表明: 配位数、剪切模量、静态结构因子以及取向序关联函数等都随分散度的变化而变化, 而力的累积分布不受分散度的影响. 其中, 单分散体系(亦即s=0)的静态结构因子在低波数区域(亦即低k)遵从幂律标度S(k)∝0.2k^-4/3, 各分散度下取向序关联函数峰值符合e指数变化规律g₆(r)∝ae^-r/ξ₆, 其中序关联长度ξ₆随分散度s增大而减小. 关键词： 颗粒物质 尺寸分散度 结构因子 取向序关联     

受激布里渊散射光脉冲波形的研究

从实验和理论上研究了受激布里渊散射(SBS)光脉冲波形.在Nd∶YAG调Q激光器中实验研究了SBS光脉冲波形随抽运光参数、结构参数和介质参数的变化规律，并利用SBS理论模型进行了数值模拟，实验与数值模拟结果的变化趋势基本一致.结果表明，抽运光能量越小，透镜焦距越短，镜-池间距越长，介质声子寿命越长，SBS光脉冲波形就越不容易出现调制现象.分析和讨论了抽运光参数、结构参数和介质参数对SBS光脉冲波形的影响机理. 关键词： 受激布里渊散射 脉冲波形 抽运光参数 结构参数     

时域有限差分法利用聚焦光脉冲模拟介质球微粒所受单频光阱力

用三维时域有限差分法对光镊装置中介质球微粒所受光阱力情况进行模拟。根据Richards-Wolf矢量场衍射积分公式对消球差会聚透镜像空间中电磁场分布的表示,在总场空间中实现了对聚焦光脉冲的模拟。聚焦光脉冲与介质球微粒的相互作用,通过离散傅里叶变换提取出单频成分,利用计算光阱力的麦克斯韦应力张量公式,计算单频激光对介质球的光阱力。聚焦光脉冲的引入可以使一次计算得到多个频率下的计算结果。对相同的聚焦装置下介质球受不同频率入射光的光阱力情况进行了计算。计算结果表明使用线偏振光作光源时,大数值孔径聚焦透镜和短入射波长有利于介质球的横向操纵;沿光轴方向对介质球进行纵向操纵,需要大数值孔径的物镜和与之相适应的入射波长。     

QUICK格式在湍流旋流流动数值模拟中的应用

为研究不同精度的离散格式对湍流旋流流动数值模拟结果的影响,同时应用QUICK格式和混合格式对同轴射流旋流燃烧室内的湍流流动进行了数值模拟.在采用k-ε湍流模型的条件下,QUICK格式计算得到的燃烧室内气体轴向与切向速度及轴向脉动速度均方根值分布与实验数据符合较好,而混合格式给出的数值模拟结果则与实验有一定的偏差.     

标准模型中两类中性玻色子混合圈链图传播子的重整化及其e+e-→μ+μ-反应截面

采用电弱统一标准模型, 对光子γ和中间玻色子Z₀ 混合圈链图传播子的构架方式及其重整化问题作了详细分析与讨论, 并完成了有关解析计算, 获得了由参与电弱相互作用的各种混合圈构成的γ和Z₀的重整化链图传播子的解析计算结果. 此外, 我们还将此结果应用于高能物理中备受关注的一类轻子反应e⁺e^-→μ⁺μ^-中, 获得了在γ和Z₀所构架出的重整化混合圈链图传播中e⁺e^-→μ⁺μ^-反应截面的解析计算结果, 并将本文所获得的这一理论计算结果与实验观测值作了对比分析, 发现本文理论计算结果与实验观测值符合非常好, 并由此也获得了有关辐射修正的重要信息. 本文研究结果也可为探讨一般复杂传播子的理论研究与应用提供某些参考. 关键词： 标准模型 链图传播子 重整化 辐射修正     

Fluctuations of fluctuation-induced casimir-like forces

The force experienced by objects embedded in a correlated medium undergoing thermal fluctuations-the so-called fluctuation-induced force-is actually itself a fluctuating quantity. Using a scalar field model, we compute the corresponding probability distribution and show that it is a Gaussian centered on the well-known Casimir force, with a nonuniversal standard deviation that can be typically as large as the mean force itself. The relevance of these results to the experimental measurement of fluctuation-induced forces in soft condensed matter is discussed, as well as the influence of the finite temporal resolution of the measuring apparatus.     

A new paradigm for the molecular basis of rubber elasticity

The molecular basis for rubber elasticity is arguably the oldest and one of the most important questions in the field of polymer physics. The theoretical investigation of rubber elasticity began in earnest almost a century ago with the development of analytic thermodynamic models, based on simple, highly-symmetric configurations of so-called Gaussian chains, i.e. polymer chains that obey Markov statistics. Numerous theories have been proposed over the past 90 years based on the ansatz that the elastic force for individual network chains arises from the entropy change associated with the distribution of end-to-end distances of a free polymer chain. There are serious conceptual objections to this assumption and others, such as the assumption that all network nodes undergo a simple volume-preserving linear motion and that all of the network chains have the same length. Recently, a new paradigm for elasticity in rubber networks has been proposed that is based on mechanisms that originate at the molecular level. Using conventional statistical mechanics analyses, Quantum Chemistry, and Molecular Dynamics simulations, the fundamental entropic and enthalpic chain extension forces for polyisoprene (natural rubber) have been determined, along with estimates for the basic force constants. Concurrently, the complex morphology of natural rubber networks (the joint probability density distributions that relate the chain end-to-end distance to its contour length) has also been captured in a numerical model (EPnet). When molecular chain forces are merged with the network structure in this model, it is possible to study the mechanical response to tensile and compressive strains of a representative volume element of a polymer network. As strain is imposed on a network, pathways of connected taut chains, that completely span the network along strain axis, emerge. Although these chains represent only a few percent of the total, they account for nearly all of the elastic stress at high strain. Here we provide a brief review of previous elasticity theories and their deficiencies, and present a new paradigm with an emphasis on experimental comparisons.     

The effect of electrostatic force on the evolution of sand saltation cloud

In a wind-blown sand layer, it has been found that wind transport of particles is always associated with separation of electric charge. This electrification in turn produces some electrostatic forces in addition to the gravitational and fluid friction forces that affect the movement of saltating sand particles, further, the wind-blown sand saltation. To evaluate this effect quantitatively, this paper presents a simulation of evolution of wind-blown sand grains after the electrostatic forces exerted on the grains are taken into account in the wind feedback mechanism of wind-blown saltation. That is, the coupling interaction between the wind flow and the saltating sand particles is employed in the simulation to the non-stationary wind and sand flows when considering fluid drag, gravitation, and a kind of electrostatic force generated from a distribution of electric field changing with time in the evolution process of the sand saltation. On the basis of the proposed simulation model, a numerical program is given to perform the simulation of this dynamic process and some characteristic quantities, e.g., duration of the system to reach the steady state, and curves of the saltating grain number, grain transport rate, mass-flux profile, and wind profile varying with time during the non-stationary evolution are displayed. The obtained numerical results exhibit that the electrostatic force is closely related to the average charge-to-mass ratio of sand particles and has obvious influence on these characteristic quantities. The obtained results also show that the duration of the system to reach the steady state, the sand transport rate and the mass flux profile coincide well with experimental results by Shao and Raupach (1992) when the average charge-to-mass ratio of sand particles is 60 μC/kg for the sand particles with average diameter of 0.25 mm. When the average charge-to-mass ratios of sand particles are taken as some other certain values, the calculation results still show that the mass flux profiles are well in agreement with the experimental data by Rasmussen and Mikkelsen (1998) for another category of sand particles, which tell us that the electrostatic force is one of main factors that have to be considered in the research of mechanism of wind-blown sand saltation.     

基于离散元方法的颗粒材料缓冲性能及影响因素分析

在冲击荷载作用下, 颗粒材料通过颗粒间的摩擦及非弹性碰撞可有效进行能量耗散实现缓冲作用. 本文采用离散元方法对冲击载荷下颗粒材料的缓冲过程进行数值分析, 研究不同厚度下颗粒材料的缓冲性能. 计算结果表明: 颗粒层厚度H是影响颗粒材料缓冲性能的关键因素, 并存在一个临界厚度H_c. 当H < H_c时, 冲击力随H的增加而降低; 当H > H_c时, 冲击力对H的变化不敏感并趋于稳定值. 此外, 在不同颗粒摩擦系数和初始密集度下对缓冲过程的离散元分析表明, 光滑和疏松颗粒材料具有更好的缓冲性能. 最后, 对颗粒材料在冲击过程中的力链结构和底板的压力分布进行了讨论, 以揭示颗粒材料缓冲性能的内在机理.     

Nondestructive measurement of thermal contact resistance for the power vertical double-diffused metal-oxide-semiconductor

To obtain thermal contact resistance(TCR) between the vertical double-diffused metal-oxide-semiconductor(VDMOS) and the heat sink, we derived the relationship between the total thermal resistance and the contact force imposed on the VDMOS. The total thermal resistance from the chip to the heat sink is measured under different contact forces, and the TCR can be extracted nondestructively from the derived relationship. Finally, the experimental results are compared with the simulation results.     

Determination of tip–sample interaction forces from measured dynamic force spectroscopy curves

The forces between a sharp tip and a sample are characteristic for different sample materials. A new method for quantifying the elastic tip–sample interaction forces from measured frequency vs. distance curves is presented. The dynamic force–spectroscopy curves investigated were obtained by dynamic force microscopy under ultrahigh vacuum (UHV) conditions for large vibration amplitudes with commercial levers/tips. The full non-linear force–distance relationship is deduced via a numerical algorithm, where the equation of motion describing the oscillation of the tip is solved explicitly. The elastic force distance dependence can be determined by fitting the results of a computer simulation to experimental frequency vs. distance data. The obtained force–distance curves can be compared quantitatively with theoretical models.     

Unfolding proteins under external forces: a solvable model under the self-consistent pair contact probability approximation

We extend a model of Micheletti et al. [Phys. Rev. Lett. 87, 088102 (2001)]] used to study protein con-formations to the case in which there is an external force field. Under the self-consistent pair contact probability approximation, this residue-level resolution model can still be solved under pulling forces. We implement the algorithm using heterogeneous parameters and study the force-induced unfolding of a helical segment from the protein transformylase and of the beta-stranded domains from the protein titin. The results are qualitatively consistent with the results from more expensive, atomistic dynamics simulation. Despite the mean-field-like approach, we observed a sharp and cooperative unfolding transition.     

Modeling force transmission in granular materials

The probability density function of contact forces in granular materials has been extensively studied and modeled as an outstanding signature of granular microstructure. Arguing that particle environments play a fundamental role in force transmission, we analyze the effects of steric constraints with respect to force balance condition and show that each force may be considered as resulting from a balance between lower and larger forces in proportions that mainly depend on steric effects. This idea leads to a general model that predicts an analytical expression of force density with a single free parameter. This expression fits well our simulation data and generically predicts the exponential fall-off of strong forces, a small peak below the mean force and the non-zero probability of vanishingly small forces.     

Two different kinds of time delays in a stochastic system

The steady state properties of a noise-driven bistable system are investigated when there are two different kinds of time delays existed in the deterministic and fluctuating forces respectively. Using the approximation of the probability density approach, the delayed Fokker-Planck equation is obtained. The stationary probability distribution (SPD) and the variance of the system are derived. It is found that the time delay τ in the deterministic force can reduce the fluctuations while the time delay β in the fluctuating force can enhance the fluctuations. Numerical simulations are presented and are in good agreement with the approximate theoretical results.