三维空间土壤-触土部件相互作用SPH模型与试验验证  

Three-dimensional SPH model for soil-soil engaging components interaction

作　　者：任逸 方瀚祯 刘航 刘鹏 李守太 胡嫚[1] REN Yi;FANG Hanzhen;LIU Hang;LIU Peng;LI Shoutai;HU Man(College of Engineering and Technology,Southwest University,Chongqing 400715,China;Chenzhi Technology Co.,Ltd,Chongqing 401120,China)

机构地区：[1]西南大学工程技术学院,重庆400715 [2]辰致科技有限公司,重庆401120

出　　处：《农业工程学报》2025年第3期53-62,共10页Transactions of the Chinese Society of Agricultural Engineering

基　　金：重庆市技术创新与应用发展(专项重点项目)(cstc2021jscxgksbX0002)。

摘　　要：触土部件在耕作、翻土和深松等作业中发挥着关键作用。然而,传统的力学分析方法难以准确量化土壤微观运动形态及三维空间中土壤-触土部件相互作用,在一定程度上限制了耕作效果的评估和触土部件的设计优化。该文采用光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法,在三维空间中构建了土壤-触土部件相互作用模型,研究基于土壤弹塑性本构理论,提出了基于虚拟平面和罚函数的三维接触算法,用于处理复杂空间触土部件的接触问题。采用C++语言借助GPU(graphics processing unit)加速技术在CUDA(compute unified device architecture)并行计算平台上实现了算法的数值求解。为验证模型的有效性,首先通过砂土坍塌试验验证了三维空间土壤弹塑性本构模型,并设计制作了小型土壤切削试验装置进行砂土切削试验,对比分析了数值模拟和物理试验中砂土切削力的变化规律。该研究还讨论了不同耕作深度、耕作速度、触土部件入土角对触土部件切削力的影响,探讨了模型在农业实践中的适用性。结果表明,所提出的模型不仅能够准确描述土壤切削形态,在切削力数值预测方面的相对误差仅为7.02%,展现出较高的精度。研究结果可为土壤-触土部件相互作用模型的建立提供思路,为农业机械触土部件的设计优化提供了新的数值分析工具。Agricultural operations(such as tillage,soil turning,and deep loosening) significantly affect the soil structure,aeration,and nutrient distribution in crop production.The performance of soil-engaging components during operations can directly impact the agricultural machinery and the quality of soil preparation.Traditional mechanical analysis has been limited to capture the micro-movements of soil particles and the dynamic interactions in three-dimensional space.It is very necessary to accurately assess the tillage in the optimization of soil-engaging component.In particular,the SPH method can also offer a mesh-free approach.The complex behavior of soil can be expected to effectively simulate,including large deformations and interactions with mechanical components.In this study,the more accurate and comprehensive model was developed for the soilengaging component interactions using SPH method.A three-dimensional contact algorithm was also incorporated using soil elastoplastic constitutive theory and penalty functions.The elastoplastic constitutive model was accurately represented the behavior of soil under various stress conditions,including elastic deformation and plastic yielding.The contact algorithm included the virtual planes and penalty functions.There was the complex contact behavior between soil particles and soilengaging components.The algorithm was implemented in C++ on the CUDA(Compute Unified Device Architecture) parallel computing platform,thus leveraging the power of GPU(Graphics Processing Unit) acceleration.Numerical simulation was conducted using Visual Studio 2019.A series of experiments were also conducted to validate the accuracy of the improved model.The sand collapse experiment was firstly conducted to create a controlled environment,where the sand was allowed to collapse under its own weight.The diameter of collapse extension was theoretically calculated using empirical formulas.A comparison was also made on the simulation and the observed one.Theoretical diameter of collapse extension was 0

关 键 词：土壤切削 弹塑性本构模型 触土部件 三维接触算法 光滑粒子流体动力学 

分 类 号：S222[农业科学—农业机械化工程]