用粒状体力学方法的研究沙的理论负功区.pdfVIP

· · 333 用粒状体力学方法研究沙的理论负功区 武建勋1，高尔新1，张海涛2，琚素芳2，黄建鸿2 中国矿业大学北京控区资源系．北京100083 摘要 粒状体力学方法为砂和砾建立了一种理想化的模型．可以将眇和砾看作很多小颗粒的集合体，并铂：计算机 中模拟的力学表现．此法特．男嵫含于研究离散体性质特别突：lj的情况，例如颗粒的飞散或砂土的后屈服特抓。奉 义用粒状九素法对Hill稳定条件进行了较深入的研究，发现非关联流动则在很多场合成为粒状体理想本构关系的 很好的近似．健只是接近而不是完全遵守非关联流动则。所以理论J-t#-关联流动则战该发生的Hill失稳现象实际 I．并不会发生。本文关于粒状体力学方法的理论根据提沁了“几率场模拟”的解释． 关键词：粒状元素泫，粒状体．本构关系．数值模拟，非关联流动则 1导言 粒状体材料指由很多小颗粒组成的集合体，可以包括沙，土，岩币^ 粮食．散装煤炭，洗农粉，食盐等。粒状体不是连续体，未必均匀，成分构成 复杂．试验误差人，影响冈素多而且缺乏成熟的理论。粒状体力学方法从 微观机理出发提供了一种理想化的模型，便于作某些深入研究。 Rowe往提出沙的”锯齿模型”时曾用十几个圆球的集台作剪切试验研 究“剪涨”现象，计算机出现后删数值化的圆盘或圆球的集合模拟颗粒的集 魁～l网盘组成的粒状体 合就更方便。Cundall(1971)提出离散单元法(DEM)。岸野111(1990) NNA,J 的运动，易控制应变而应力控制不佳，就产生上述差别。GEM法可以用于动力学和运动学，不要求粒子 保持接触，可以计算粒子的ls散平“碰撞，所以应用更7‘泛。 粒：状体力学方法可以记录颗粒运动和接触的一切细+饥可以随 本，可以随意调1y试验条什．这些都是传统的土力学试验难于做到 就是一个{良难州十力学试验实现的本构关系研究的例子。 2．粒状体力学一般方法 图一1所示平面颗粒集合由395个大小不等的圆盘形的粒子组成 一1．0cm．，接触点性质可概略表示如图．2，其中摩擦力满足．|=孳仑法!il ；P，}≤p。，tg巾 (1) r)。≥O (2) 这里P。P。分别为接触点的法向和切向的接触力，法向力受乐为豇二 l武堪撕 l㈨H．I}t：乍】，5．T学媾{：．救授．结构．岩上 9i．I，：}…’l-． 1．岛尔似 男．T学溥l：．救授．结构．爆破 1q#1．……’l ：帐淞涛 粥．硕1．研究!h绌拗 2 f城羝嚣 1977．i…，1址．1ij{Ij研宄!l÷．结构 1970．1：l：’I ：c^缝鸿 粥．碗1．研究瓠绻树 ·334· △LI．=k。P。 (3) Au，=ktP， (4) 摩擦角q’=25。。在平衡状态内部粒子在各接触点的作刚力下平衡，周边粒 子除其他粒子接触力还受到外界的作用力，同一边界粒子所受各力也应该 平衡。颗粒集合所受外界力的综合可以定义出该粒状体内的平均应力，根 据在连续体情况下的等价性经常定义为 oi．i一--(1／V)∑(xifj) (5) 这里V是面积(二维)或体积(三维)，指所有周边粒子中心连线所围住的 而积或体积，x；为周边粒子的坐标，fJ为周边粒子所受外力的合力的分撼， ∑对所有周边粒子求和，i，jl，2(二维)或i。j--1，2，3(三维)。 颗粒集合受到外力后每个粒子会发生变位，由这些变位的综合可以 定义出该粒状体内的平均应变，根据在连续体情况下的等