读书笔记2..pptxVIP

读书笔记2郑泽松Particulate discrete element modelling——chapter2运用离散元法研究土体管涌机理——第4、5章管涌现象细观机理的模型试验与颗粒流数值模拟研究——第1章Particulate discrete element modellingCatherine o’SullivanChapter 2Tessellation 曲面细分离散颗粒颗粒系统颗粒间的接触初始样本几何构建useful系统三角网应力边界的应用应变计算材料的组构分析使用三角网将DEM颗粒代码与连续体力学结合起来是可能的表示液相Tessellation 曲面细分： 将一个空间分成一组子空间，这些子空间无重叠且完全充满原空间（无空隙） 各个Tessellation密切相关，每一个称为另一个的 “对偶dual”Delaunay triangulation复杂几何体的有 限元分析的网格生成useful一个n点集合P={P1，Pk,Pn}的一个三角网是由一组m三角形集合T={T1,TK，Tm}组成的，其内部不相交点集的一个delaunay三角网，特性： 点集内没有落在三角网内任意三角形外接圆内的点 顶点集的delaunay三角网是唯一的可推广至三维一个n点集合P={P1，Pk,Pn}的voronoi图是m多边形集合V={V 1，Vk，Vm}每个多边形Vk是围绕响应点Pk的每个多边形Vk包含一块区域或体积（a）十个随机点 （b） Delaunay三角网 （c） voronoi图舍入误差：以有限位的浮点数来表示实际数字所产生的误差在精确求解接触几何或时间积分法时，注意算法选择Chapter 2 particle motion离散元分析是动态或瞬时的分析考虑颗粒接触的系统内的动态相互作用非线性分析 几何非线性 颗粒集合体几何形状改变 接触时力位移的非线性通过直接考虑独立颗粒的动力平衡来引入DEM基本法则： 矩阵结构分析 有限元分析颗粒与 矩阵结构分析中的随机程度相似 有限元网格中的结点相似系统动力平衡方程M：重力矩阵C：阻尼矩阵U:位移（平移/转动位移） ：力的增量向量K:整体刚度矩阵DEM中的颗粒与有限元分析中的结点相似但是。DEM中颗粒可自由转动2维 三个自由度2平移 1旋转 3维 六个自由度 3平移 3旋转求“多点”系统动力平衡方程的两种方法 含蓄法implicit 明晰法explicit对系统所有颗粒中心点用单一向量u表示联合位移增量与有限元分析中用单一向量表示所有结点的位移相似在cundall和stract的明显单元法和分子动力学中，方程的全局系统的解法，是通过考虑独立颗粒的动力平衡而不是同时解决整个系统。这种方法避免了巨大全局刚度矩阵的创建和存储。一个质量mp的颗粒p的平移动力平衡一般方程：式中， ——当颗粒p与其他颗粒或边界有Nc,p次接触，由接触c引起的接触力 ——当颗粒p与其他颗粒或边界有Nnc,p次接触的非接触力（不饱和土 毛细力） ——液体作用在颗粒p上的相互作用力 ——颗粒p的重力 ——颗粒p的具体应用力（可能在应力控制膜中出现）动态旋转平衡方程 ——角速度向量Mpj——在总共Nmom个传递接触的动量里，第j个传递涉及 颗粒p接触力的动量接触点处接触力 切向力 正交向力假设应用力及扭矩已知，据上述两个公式求解颗粒的平移和旋转的加速度时间积分法time integration methods用中心差异时间积分法central-difference time integration approach渐进性求解下式：——通过颗粒中心，垂直分析平面的轴线的角速度角速度 接触力的切向分力 更新非球形颗粒边缘的位置 全部颗粒的旋转（作为颗粒在材料内定位的 重要指标）关键： 的选择 计算由泰勒展开式所成的截断误差与 成比例， 越小，如 1，n越大，误差越小在给定时间点产生误差E0，n次计算后的误差En,称整体误差En 的线性增长是不可避免的若En与局部误差local error的关系：En=CnE0,C为常数，则算法稳定；若En=C nE0 ,n1，则算法不稳定通过计算系统全部能量来监测数值模型稳定性若系统能量没有流动，则数值积分稳定不稳定的系统中，能量不守恒接触弹性应变能颗粒的动能对单一自由度系统，动力平衡方程：a=-kx, k为刚度对此简单系统，若用中心差异法，最大时间增量为：T为系统自由振动周期若时间增量 超过上述临界值，预测结果将很快变得没有物理意义，分析就是不稳定的使用考虑放大（增强）矩阵的线性稳定分析来计算稳定分析的临界时间增量增强矩阵A定义：若A的任何特征值ui超过1，则任何初始小误差都将无限增长，分