圆盘DDA及无网格法耦合DDA理论应用

《圆盘DDA及无网格法耦合DDA理论应用》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、国内图书分类号：O346.1密级：公开国际图书分类号：531西南交通大学研究生学位论文圆盘DDA及无网格法耦合DDA理论应用年级二〇一五级姓名葛发乐申请学位级别硕士专业力学指导老师喻勇教授二零一八年五月ClassifiedIndex:O346.1U.D.C:531SouthwestJiaotongUniversityMasterDegreeThesisDisc-basedDDAandtheApplicationofMeshlessMethodandDDACoupledModelGrade：2015Candidate：GeFaleAcademicDegreeAppliedfor

2、：MasterDegreeSpeciality：MechanicsSupervisor：Prof.YuYongMay,2018西南交通大学学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权西南交通大学可以将本论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复印手段保存和汇编本学位论文。本学位论文属于1．保密□，在年解密后适用本授权书；2．不保密□，使用本授权书。（请在以上方框内打“√”）学位论文作者签名：指导老师签名：日期：日期：西南交通

3、大学硕士学位论文主要工作（贡献）声明本人在学位论文中所做的主要工作或贡献如下：（1）本文介绍了刚性圆盘和边接触的处理方法，使用罚函数法推导了接触子矩阵，并针运动方向的不同推导了摩擦力子矩阵。针对原DDA中判断圆盘与凸角、凹角接触时出现多加接触弹簧的问题，提出了新的圆盘与凸角、凹角的接触判断，把角平分线法应用到接触判断中去，编写程序验证角平分线法的准确性与有效性。通过分析圆盘在斜面运动时分别与凸角、凹角接触时系统能量计算值与理论值的误差，验证了角平分线法的实用性。（2）针对使用非刚性圆盘离散的悬臂梁模型，在相同约束下分析圆盘尺寸对计算结果精度的影响。当圆盘半径由0.1m减小到0

4、.015m时，圆盘个数由原来的360个增加到16000个，应力计算精度大大提高。但圆盘数的增加，使得平衡方程维数的大大增加，计算时长骤增。（3）对DDA计算结果使用点插值法和移动最小二乘法进行插值后处理，结果证明插值结果依赖于圆盘尺寸的选取。圆盘尺寸过大或过小均会影响插值结果，圆盘尺寸过大，插值点过少，无法得到更高精度的插值解过，圆盘尺寸过小，将会使插值点失去特性而过于光滑，亦会增大误差。算例证明，在本文使用的悬臂梁离散化模型中，圆盘半径为0.03~0.02m时插值效果最好。（4）使用无网格位移模式改进DDA方法，推导了总体平衡方程和各个子矩阵。研究加载速率对计算结果的影响，

5、并绘制悬臂梁的应力云图。推导了块体接触矩阵。本人郑重声明：所呈交的学位论文，是在导师指导下独立进行研究工作所得的成果。除文中已经注明引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写过的研究成果。对本文的研究做出贡献的个人和集体，均已在文中作了明确说明。本人完全了解违反上述声明所引起的一切法律责任将由本人承担。学位论文作者签名：日期：西南交通大学硕士研究生学位论文I摘要原DDA理论在计算离散体运动学速度和位移时较准确，但对于块体间接触或边界接触，使用罚函数法会因为参数选取没有准确的标准而使计算结果因人而异，以致在模拟块体系统运动时块体间的总能量计算不准确，在有摩擦力的

6、情况下误差更大。本文使用较准确的块体位移公式修正原DDA中因块体转动而造成的块体面积虚假膨胀，使用角平分线法修正圆盘与凸角、凹角的接触，并编写MATLAB程序进行验证。研究了圆盘离散悬臂梁模型，探讨了不同圆盘半径对自由端受荷载的悬臂梁计算精度的影响，使用点插值法和移动最小二乘法对DDA计算结果进行了后处理，并推导了相关公式。使用无网格位移模式改进DDA方法，推导了系统总体方程和各个子矩阵，研究了加载速率对计算结果的影响，并编写程序绘制了自由端受荷载的悬臂梁的应力云图。主要结论如下：（1）本文介绍了刚性圆盘和边接触的处理方法，使用罚函数法推导了接触子矩阵，并针对不同的运动方向推

7、导了摩擦力子矩阵。针对原DDA中接触判断不准确的问题，提出了新的圆盘与凸角、凹角接触判断，将角平分线法应用于接触判断，并编写程序验证了角平分线法的准确性与有效性。通过分析圆盘在斜面运动过程中分别与凸角、凹角接触时系统能量计算值与理论值的误差，验证了角平分线法的实用性。（2）针对使用非刚性圆盘离散的悬臂梁模型，在相同约束下分析圆盘尺寸对计算结果精度的影响。当圆盘半径由0.1m减小到0.015m时，圆盘个数由原来的360个增加到16000个，应力计算精度大大提高。但圆盘数的增加，使得平衡方程维数大大增加，计