粘湿土壤力学行为的离散元模拟.pdfVIP

粘湿土壤力学行为的离散元模拟1 张锐，李建桥 吉林大学地面机械仿生技术教育部重点实验室(130025) E-mail: zhangrui@ 摘 要：本文系统分析了粘湿土壤的微观力学结构，并指出了水分对粘湿土壤力学行为的影 响。根据土壤颗粒的基本力学性质以及水分对粘湿土壤颗粒的作用，建立了粘湿土壤离散元 力学模型，并使用离散元专用软件PFC2D®对推土板与粘湿土壤之间相互作用的动态行为进行 了模拟。模拟结果与真实粘湿土壤在外力作用下的情况十分相似。 关键词：粘湿土壤 液桥 离散单元法 1. 概述 粘湿土壤是耕地和施工中经常遇到的土壤类型，分析粘湿土壤与触土部件之间的力学相 互作用以及外力作用下粘湿土壤内部的动态行为，对于设计和优化减粘降阻的地面机械具有 十分重要的意义。粘湿土壤内部含有水分，水分的存在使土壤颗粒之间具有粘附性，这使粘 湿土壤的力学结构与普通土壤相比更加复杂。只是通过土壤运动过程的宏观现象建立离散元 仿真模型，不足以描述粘湿土壤的力学行为。充分考虑粘湿土壤的微观结构，会使土壤离散 元动态仿真结果更加符合实际情况。 离散单元法(Discrete Element Method, DEM)是Cundall P.A.教授在1971年基于分子动力 学原理首次提出的一种不连续数值模拟方法[1] 。这种方法克服了连续机理模型(Continuum Mechanics Method, CMM) 中忽略单元个体性质，过分依赖高度简化的、规定性质的本构方程 的不足，在分析具有离散性质的物质上表现出了极大的优越性。作为多相混合体的土壤本身 即是离散结构，被触土部件切断或分离的土壤更是不连续的，所以离散单元法是一种分析象 土壤耕作这样的条件下，土壤巨大非连续变形的理想方法[2] 。 2. 水分对粘湿土壤微观力学结构的影响 根据饱和程度的不同，颗粒集合体内空隙之间的液体以摆动、链索、毛细或浸渍状态存 在[3] 。粘湿土壤的水分含量介于干性土壤和饱和土壤之间，水在粘湿土壤中以非均匀、非连 续形式分布，所以水在粘湿土壤内部土壤颗粒空隙之间以摆动状态存在。在摆动状态下，土 壤颗粒之间的水分以不连续的液桥分布，如图 1 所示。Fisher R.A.认为颗粒之间的液桥可以 近似认为由两个相互垂直的弧形曲面组成，两土壤颗粒之间的液桥模型如图 2 所示[4] 。粘湿 土壤颗粒空隙之间不连续的液桥会产生 “静态”毛细管力和 “动态”粘附力。毛细管力产生 于液桥静水压力差∆P和液桥表面张力γ；动态粘附力包括空隙水分径向压流产生的法向分量 和空隙水分剪切流产生的切向分量[5] 。根据图 2 所示的几何关系，可以求得凹面半径ρ ，凸 1 面半径ρ2分别为： 1本课题得到国家自然科学基金 ；教育部博士点基金 （1999018504）；教育部高等院校骨干教 师基金（[2000]65-1）资助。 - 1 - s R + (1−cosφ ) i i i ρ1 ρ R sinφ −ρ [1−sin(φ +θ)] i=1 或 2 （1） 2 i i 1 i cos(φ +θ) i 式中 R ——两颗粒半径；s ——两颗粒到液桥颈部的距离；θ——颗粒的水接触角；φ—— i i