基于深度学习和细观力学的颗粒材料本构关系研究

颗粒材料的本构关系对岩土工程等众多领域至关重要. 不同于传统的唯象本构理论, 本文基于机器学习模型探索了一种细观力学理论指导下的数据驱动型颗粒材料本构关系预测方法. 根据Vogit均质化假设, 建立了小应变条件下颗粒材料应力?应变解析关系, 此关系唯一地确定了一组与颗粒材料本构行为相关的细观组构变量. 这些变量与反应颗粒材料宏观性质的主应力和主应变信息通过一系列离散元三轴压缩数值试验获得. 考虑到细观组构变量为内变量, 不能直接作为本构模型的输入. 本文基于有向图方法将颗粒材料微观结构信息隐式地包含在应力?应变的预测当中, 并采用门控循环单元(GRU)循环神经网络作为基础深度学习模型描述有向图中结点之间的映射关系. 通过将有向图从目标节点沿源节点展开, 整个应力?应变预测模型可由两个神经网络分别训练并组装而成. 将训练后的深度学习模型在全新的数据集上进行测试, 结果表明该训练策略能有效捕捉到颗粒材料在常规三轴任意加卸载, 等中主应力系数b的真三轴加载, 和等平均有效应力p的真三轴加卸载等复杂多轴加载工况下的应力?应变响应关系, 模型具有良好的内插和外推预测能力. 考虑到深度学习模型捕捉颗粒材料力学响应的能力及其开放式学习的特点, 充分结合数据驱动方法和理论本构模型可能是颗粒材料本构研究的一个重要方向.      

基于深度强化学习算法的颗粒材料应力?应变关系数据驱动模拟研究

颗粒材料的宏观力学行为受颗粒组分等材料参数, 孔隙率、配位数等状态参数的影响, 同时又具备复杂的加载路径和加载历史相关性, 建立包含多个内变量以及各变量间相互关联的颗粒材料本构模型是一个重要的科学难题. 不同于传统的基于屈服面、流动法则和硬化函数框架下的唯象本构模型, 本文基于颗粒物质力学的研究基础, 以颗粒材料平均孔隙率、细观组构参数和弹性刚度参数作为内变量, 结合深度学习方法建立以有向图表征的数据本构模型. 有向图中以不同的链接网络表示不同的内变量信息流动方向, 各个内变量间的映射关系采用循环神经网络来建立, 将各个神经网络相互组合, 形成包含不同内变量且具有不同预测能力的本构模型. 该本构模型的建立过程等价于在众多可能的内变量链接关系空间中寻找最能描述实际材料宏观应力应变行为的优化问题. 因此, 可将有向图本构模型的建立过程看作“马尔可夫决策过程”, 采用深度强化学习算法构建有向图的内变量链接组合优化过程, 具体采用AlphaGo Zero算法自动寻找最优的颗粒材料数据驱动本构模型建模路径. 研究结果表明, 采用有向图和深度强化学习算法可建立起完全依靠“数据驱动”的颗粒材料应力?应变关系. 此外, 本方法提供了一种将不同理论模型从数据角度统一起来, 且基于人工智能算法发展更优模型的研究思路, 可为相似问题的研究提供借鉴.      

船-冰作用中海冰典型破坏模式的离散元数值模拟

在破冰船破冰过程中,冰排主要表现为挤压与弯曲两种破坏模式.本文基于黏结离散单元法对船-冰作用中的这两种典型破坏模式进行数值模拟.海冰离散元数值试样采用随机排布方式生成,采用单轴压缩试验与三点弯曲试验相结合的方式标定模型中的细观参数.将船-冰碰撞中的挤压和弯曲作用方式简化为直立或倾斜 平板与海冰的作用模式,构建挤压与弯曲破坏的海冰离散元数值试样,分析了挤压破坏模式中不同加载方位以及冰厚和加载速率对破坏模式的影响,以及弯曲破坏模式中不同冰排夹角以及冰厚和加载速率对破坏模式的影响.计算结果表明,离散元海冰数值模型可以很好地对冰排挤压与弯曲破坏现象进行模拟,可揭示冰排在不同破坏模式下的破坏机理.     

风沙影响下铁路道碴变形模量的离散元数值分析

针对风沙区有碴铁路道床的结构特性,对其在细沙贯入下的有效变形模量研究有助于理解风沙影响下的道床动力特性。采用离散单元模型对道碴碎石和细沙颗粒进行数值建模,并对不同含沙率γ下的有效变形模量进行了数值分析。结果表明,在低含沙率下(γ<30%)沙石混合体的变形模量基本保持不变,在高含沙率下(γ>30%),变形模量随含沙率的增加呈线性降低,以上变形模量随含沙量的变化规律与试验结果相一致。基于离散单元模型的数值模拟,在细观尺度上对沙石混合体的力链强度、空间分布及配位数进行了分析,揭示了沙石混合体有效变形模量随含沙量变化的内在机理。本文工作对风沙区有碴铁路道床力学行为的研究具有一定的借鉴意义,有助于促进风沙影响下有碴铁路道床的结构设计和沙害治理。     

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

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