静态颗粒体系力分布的系综理论.pdfVIP

第七届全国液体和软物质物理学术会议摘要 OA09 静态颗粒体系力分布的系综理论 张兴刚，胡林 贵州大学理学院，贵阳550025 近年来，J．H．Snoeijer等人用力网系综理论研究力分布问题，取得了许多与实 验相吻合的结果【¨l。由于力网系综理论基于等概率的假设，涉及到复杂奇异函数 的积分，所以前人主要采用MonteCarlo的方法进行数值模拟。为了对力的非单调 分布的机理有更深入的认识，需要提出能解析计算的模型。另外，也需要给出关于 奇异函数积分的新方法。 本人提出标量力网系综模型来研究力分布问题，模型中主要运用力网系综的理 论，借鉴了标量q模型4l的思想，只考虑颗粒体系在一个方向上的力平衡，将实际 堆积简化为图(1)所示的接触网。其中每条竖线或代表一个颗粒，由斜线∥相连 的两条竖线表示这两个颗粒相接触。不考虑外场，给出体系力平衡条件、边界条 件及非负约束后，可写出力向量的几率密度： 碟(雷固2丽1万(F一善g勉+1)。‘取(8(92是_i-乞2Hn-I—c箍∥(珐+-一乞2·n一一磙) (1) 艿(如一靠。一c黔11)8(gh一露一靠+11))】．【(兀秒(g：))(n口(√))】 m．n=l ，．J*1 而某个力醮的分布： J尸MNLg-，否) (2) p(／)=[：秘艿(／一繇n 对于c?的分布也一样。为了得到某个力的几率分布，需要求出(2)式中关于Dirac．6 函数、Heaviside．0函数的积分。这个问题的一般情形已得到解决，本人构造了几个 关于0函数的定理，利用它们可以直接计算任意情况下(2)式的积分问题。分析表明： (2)式积分的结果是若干个带有0函数的多项式的求和；计算过程及结果的复杂 程度随体系规模增大而急剧增加。直接计算的方法不便于对极限分布问题进行讨 43 第七届全国液体和软物质物理学术会议摘要 论。为此，提出了Laplace变换法计算(2)式的积分，解决了小M值的情况下力分 布的问题，得到1xN、2xN、3xN体系中力分布的严格解和相应极限分布函数的形 式。 m 图l 颗粒体系接触网：竖线代表颗粒，斜线表示颗粒问有接触 根据解析计算的结果得到以下结论：体系层数大于1时力分布都是非单调变 化的，原因是体系中层与层之间存在关联；极限分布一般为多项式与指数函数的 乘积，力分布在小于平均力时呈多项式增加，大于平均力时呈指数减小，接近平 均力处取最大值，这与相关实验和计算机模拟的结论一致‘5捌。 参考文献： 【l】Jacco J．H．Vlugt,et H．SnoeUer，ThUs a1．Phys．Rev．Lett，2004，92(5)：054302(4)． l 【2】JaccoH．Snoeijer,ThijsJ．H．Vlugt,eta1．Phys．Rev．E，2004，70(6)：06306(16)． Ellenbroek，et 【3】J．H．SnoeUer，W．G a1．Phys．Rev．Lett．2006，96：098001． 【4】s．N 996，53(5)：4673-4685． 079-1082． 【5】T．S．Majmudar,R．EBehringer．Nature，2005，435(23)：1 【6】孙其诚，王光谦．物理学报，2008，57