月球着陆器着陆过程的DEM-FEM耦合分析

实现月球着陆器的安全软着陆是宇航员、搭载设备、仪器安全以及着陆器后续正常工作的重要保证。本文采用离散元-有限元方法建立月壤与着陆器相互作用的耦合模型，其中月壤采用具有粘接作用的球体离散元单元，着陆器采用壳单元和梁单元组合的有限元模型，与缓冲垫连接的支撑腿采用可压缩弹簧模型，以实现着陆器自身的缓冲功能；为分析着陆于不同月表的过程，建立平坦、斜坡月球表面；针对不同坡度的斜坡月表，研究不同着陆模式下着陆器与月壤的相互作用，分析冲击力峰值的大小与作用时间的关系，并从能量方面讨论和解释大坡度下冲击力峰值小的原因。以上研究为月球着陆器的安全着陆分析提供有益的参考。     

不同载荷频率下有砟铁路道床累积沉降特性的离散元分析

为研究有砟道床在高频载荷作用下的累积沉降规律,建立工程尺度下枕木-有砟道床的离散元计算模型。考虑道砟的非规则真实几何形态,利用不可破碎的簇颗粒模型近似构造。在10 Hz～100 Hz载荷频率区间内,对比分析枕木位移以及道砟颗粒在垂向和横向上的运动分布规律,研究分析道床在循环正弦载荷作用下的动力特性。数值结果表明,在载荷频率不超过60 Hz时,道床累积沉降随着载荷频率的提高而逐渐增大;当载荷频率超过60 Hz时,道床累积沉降随着载荷频率的提高而急剧增大,并在载荷频率为80 Hz时达到最大值,而后道床累积沉降量随着频率的提高反而逐渐减小。此外,分析对比不同载荷频率下道床速度场分布及典型位置道砟的加速度等物理量,发现80 Hz的载荷频率下道床离散元模型的振动最强烈。     

海冰与自升式海洋平台相互作用GPU离散元模拟

在海冰与自升式海洋平台结构的相互作用过程中,冰载荷是影响平台结构振动响应和疲劳寿命的重要因素. 采用具有粘接-破碎效应的离散元模型,可对海冰与自升式海洋平台结构作用中的海冰破碎特征及相应冰载荷进行数值分析. 针对自升式海洋平台的多桩腿结构特性及其冰载荷离散元分析的大规模计算需求,建立了基于GPU 的并行算法并开发了相应的计算程序. 为实现离散元分析的高效计算,采用网格排序方法创建单元邻居列表,以快速确定海冰单元间及其与平台结构间的接触模式和作用力. 此外,还发展了球体单元与圆柱形结构在不同接触形式下的计算模型. 为检验该离散元模型的有效性,对渤海锥体海洋平台结构的作用过程进行了计算,并与现场实测冰力数据进行了对比验证. 在此基础上对多桩腿自升式平台结构的冰载荷进行了离散元分析,获得了海冰的破坏特性,确定了不同桩腿上的冰力时程. 该模型可进一步应用于不同类型海洋结构的冰载荷分析,为冰区海洋平台的结构设计和现役平台结构的疲劳分析提供参考依据.      

碎冰区海冰与船舶结构相互作用的离散元分析

采用离散元模型对碎冰区海冰与船舶结构的相互作用进行了数值模拟。碎冰由三维圆盘单元构成,并考虑了其在海流作用下的浮力、拖曳力、附加质量;采用单元组合的方式构造了船体的简易模型。通过海冰与船体单元间的接触判断和接触力计算确定海冰与船体结构之间的相互作用。采用以上离散单元模型对不同冰速、冰厚、冰块尺寸、密集度条件下海冰与船体的动力作用过程进行了数值计算,并对比分析了不同海冰参数下冰对船体作用力的幅值。结果表明：冰对船体的作用力随冰厚、流速、密集度、冰块尺寸的增加而增大。本文工作可为冰区船舶的安全运行和结构设计提供参考。     

水底管道抛石加固过程的DEM-SPH耦合分析

水底管道的抛石加固过程是典型的颗粒-流体耦合问题.采用DEM-SPH耦合方法模拟颗粒-流体系统,其中离散元方法(DEM)用于模拟落石,光滑粒子流体动力学方法(SPH)用于模拟流体.通过三维Voronoi切割算法生成不规则形状的多面体,并基于闵可夫斯基原理构造扩展多面体形态的落石单元.通过SPH的边界排斥力模型计算颗粒与流体间的作用力,从而建立DEM-SPH耦合方法.采用该方法模拟溃坝与楔形块入水的过程,将计算结果与试验结果及其他数值结果进行对比分析,分别验证了 SPH与DEM-SPH耦合方法的合理性.建立锥形结构模拟卸料斗和水底管道,采用DEM-SPH耦合方法模拟落石通过卸料斗入水并与水底管道相互作用的过程,确定了落石和水对管道的作用力,并分析了卸料斗静止与运动时的落石堆积情况.以上研究表明,DEM-SPH方法可有效模拟颗粒材料、水和工程结构的相互作用,可进一步应用于水下抛石过程的结构和参数设计.     

非均匀颗粒材料的类固-液相变行为及本构方程

以非均匀颗粒介质为研究对象，采用三维离散元方法对其在不同密集度和剪切速率下的动 力过程进行了数值模拟，分析了其在由瞬时接触的快速流动向持续接触的准静态流动的转变 过程及其行为特点. 通过对不同材料性质下相变过渡区内颗粒材料的宏观应力、接触时间数、 配位数、团聚颗粒数量、有效摩擦系数等参量的计算，更加全面地描述了非均匀颗粒材料在 类固-液相变过程中的基本特征. 基于以上数值计算结果，建立了一个适用于颗粒材料 类固态、类液态以及其相变过程的本构方程，并通过剪切室实验结果验证了它的合理性.     

岩土类颗粒物质宏-细观力学研究进展

岩土类颗粒物质在自然界、工程建设以及日常生活中普遍存在,其运动特性的研究在力学界已经开展了几十年.在近20年开展的一系列小尺寸物理实验中,颗粒物质表现出许多新奇现象,人们从物理角度开展了系统研究,在统计力学中,颗粒固体的流体动力学等理论研究以及实验检测技术等方面都取得突破性进展,深刻地揭示了颗粒材料的物理机制,促使力学...     

含石量对土石混合体压桩承载力影响的离散元分析

土石混合体是介于土体和岩体之间的一种非均质、非连续和非线性的特殊工程地质材料,其在压桩贯入过程中的承载力受含石量的影响非常显著。本文分别采用球形颗粒和非规则镶嵌组合颗粒模拟土体颗粒和块石,对不同含石量下压桩贯入过程进行离散元数值分析。计算结果表明,桩柱阻力及其波动规律在不同含石量下有很大的差别。高含石量下的阻力要明显大于低含石量下的阻力,且其波动性也更加明显。通过对土石混合体内部力链结构的细观分析,揭示了压桩贯入过程中承载力随含石量变化的内在机理。以上研究有助于分析土石混合体材料的宏观力学行为,深入研究其在复杂工程条件下的力学特性。     

考虑等效曲率的超二次曲面单元非线性接触模型

基于连续函数包络的超二次曲面单元可有效地描述自然界和工业生产中的非球体颗粒形态,并通过非线性迭代方法精确计算单元间的接触力.对于具有复杂几何形态的超二次曲面单元,线性接触模型不能准确地计算不同接触模式下的作用力.考虑超二次曲面单元相互作用时不同颗粒形状及表面曲率的影响,本文发展了相应的非线性黏弹性接触模型.该模型将不同接触模式下的法向刚度和黏滞力统一表述为单元间局部接触点处等效曲率半径的函数;切向接触作用则借鉴基于Mohr-Coulomb摩擦定律的球体单元非线性接触模型的计算方法.为检验超二次曲面单元接触模型的可靠性,对球形颗粒间的法向碰撞、椭球体颗粒间的斜冲击过程、圆柱体的静态堆积和椭球体的动态卸料过程进行离散元模拟,并与有限元数值结果及试验结果进行对比验证.计算表明,考虑接触点处等效曲率半径的超二次曲面非线性接触模型可准确地计算单元间的接触碰撞作用,并合理地反映非球形颗粒体系的运动规律.在此基础上进一步分析了不同长宽比和表面尖锐度对卸料过程中颗粒流动特性的影响,为非球形颗粒材料的流动特性分析提供了一种有效的离散元方法.     

转动驱动圆角立方体颗粒有序堆积的离散元模拟

宏观颗粒体系的有序堆积研究可以为热系统中微观粒子的自组装提供研究模型,也有助于工业生产中提高颗粒材料的填充率.实验发现圆筒中的圆角立方体颗粒受容器往复旋转剪切作用会实现有序堆积.为探究旋转圆筒中立方体颗粒有序堆积过程的内部结构演化过程和动力学机理,采用超二次曲面方程构造了圆角立方体颗粒,基于离散元方法对旋转剪切作用下圆...