砂土和粘性土的颗粒流模型及其试验研究.docxVIP

砂土和粘性土的颗粒流模型及其试验研究 0 基于颗粒流的土体微观观模拟 土壤结构可以被视为由单粒、收集粒、凝块等骨架单元组成的空间结构体系。单元的形状确定了力的传输性能和土壤变形性质的性能。通过连接方法确定了土壤的结构强度，排列方法确定了土壤的稳定性。这种纤细剖面结构的复杂性决定了土壤的复杂工程特征。因此，对于土壤的这种复杂工程特征的研究，可以从薄微观的力学方面开始。然而，这项研究非常复杂，很难建立宏观和微观的定量关系。传统的连续介质力学模型的宏观连续性假设可以从薄微观层面模拟土壤工程特征，并通过详细观测参数的研究来分析宏观力学行为。国外的dannin和appellee（1991.1992）使用理想的二维线衍射模型模拟了骨料中条带状结构，并研究了带的厚度、带内位移、间隙比、体变形、旋转等。多孔方程模型（1999）用于模拟砂岩和混合土壤之间的张力-反应关系，验证了周健等人（2000.2002）提出的砂模和粘土条件下的张力-反应关系曲线。 本文通过引入颗粒流理论和开发颗粒流数值模拟技术, 对土体应力应变关系和剪切带形成机理进行细观数值模拟, 将土体微细观结构与宏观力学反应联系起来, 对土体工程力学特性和剪切带形成与发展等渐进破坏过程有更深入的了解和发现. 1 接触部分的运动定律 颗粒流理论在整个计算循环过程中, 交替应用力-位移定律和牛顿运动定律.通过力-位移定律更新接触部分的接触力.通过运动定律, 更新颗粒与墙 (边界) 的位置, 构成颗粒之间的新接触. 1.1 接触点刚度计算 颗粒流理论通过力-位移定律把相互接触两部分的力与位移联系起来. 接触力Fi可以分解为切向与法向分量 Fi=Fni+Fsi(1)Fi=Fin+Fis(1) 式中:Fniin为法向分量;Fsiis为切向分量. 法向分量可以根据下式计算 Fni=KnUnni(2)Fin=ΚnUnni(2) 式中:Kn为接触点法向刚度;Un为接触“重叠”量;ni为接触面单位法向量. 而切向接触力以增量的形式计算 ΔFsi=?KsΔUsi(3)ΔUsi=VsiΔt(4)ΔFis=-ΚsΔUis(3)ΔUis=VisΔt(4) 式中:Ks为接触点切向刚度;ΔUsiis为计算时步内接触位移增量的切向分量;Vsiis为接触点速度的切向分量;Δt为计算时步. 式 (2) 和式 (3) 中法向接触刚度Kn和切向接触刚度Ks, 是根据相互接触颗粒的几何参数以及接触模量确定, 在具体计算时首先根据被模拟介质特性设定某一值, 然后通过试算逼近目标值的方法确定其值. 通过迭加求出切向接触力分量 Fsi←Fsi+ΔFsi(5)Fis←Fis+ΔFis(5) 调整由式 (2) 和式 (5) 确定的法向与切向接触力, 使其满足接触本构关系. 1.2 旋转运动方程 单个颗粒的运动是由作用于其上的合力和合力矩决定, 可以用颗粒内一点的线速度与颗粒的角速度来描述.运动方程由两组向量方程表示, 一组是合力与线性运动的关系, 另一组是表示合力矩与旋转运动的关系分别如式 (6) 与式 (7) 所示. Fi=m(x¨?gi)线性运动(6)Mi=H˙i?旋转运动(7)Fi=m(x¨-gi)线性运动(6)Μi=Η˙i?旋转运动(7) 式中:Fi为合力;m为颗粒总质量;gi为重力加速度;Mi为合力矩;H˙Η˙i为角动量. 2 接触结构模型 2.1 颗粒间重叠量为6.2动时fsis评分的判别条件 滑动模型在相互接触颗粒之间没有法向抗拉强度, 允许颗粒在其抗剪强度范围内发生滑动, 该模型适用于模拟颗粒间不存在粘结力的散体材料 (如砂土) . 滑动模型是通过两接触颗粒间最小摩擦系数μ定义, 若颗粒间重叠量Un小于或等于零, 则令法向和切向接触力等于零.颗粒之间发生滑动的判别条件为 Fsmax=μ|Fni|(8)Fmaxs=μ|Fin|(8) 若|Fsiis|Fsmaxmaxs, 则可发生滑动, 并在下一循环中Fsiis为 Fsi←Fsi(Fsmax/|Fsi|)(9)Fis←Fis(Fmaxs/|Fis|)(9) 通过这样的循环迭代, 直到|Fsiis|与Fsmaxmaxs非常逼近, 确定发生滑动时Fsi的临界值. 2.2 接触连接模型 颗粒流模型允许相互接触颗粒连接在一起, 有两种连接模型:即接触连接与平行连接模型.接触连接假设连接只发生在接触点很小范围内, 而平行连接发生在接触颗粒间有限范围内.接触连接只能传递力, 而平行连接同时能传递力矩.这种模型适用于模拟颗粒之间存在粘聚力的材料 (粘性土) . 接触连接可以想象为一对有恒定法向刚度与切向刚度的弹簧作用于颗粒接触点处, 并假设这些弹簧有一定的抗拉强度与抗剪强度.当接触连接存在时颗粒间没有滑动, 即切向接触力不满足式 (8) .接触连接模型当颗粒间重叠量Un0时, 允许出现