基于SPH算法的深松铲破坏土壤仿真模型_张志军.pdf

收稿日期:2021-04-08 修回日期:2021-05-28 第 40 卷 第 4 期计 算 机 仿 真2023 年 4 月 文章编号:1006-9348(2023)04-0290-05基于 SPH 算法的深松铲破坏土壤仿真模型张志军,高奕珏(常州大学环境与安全工程学院,江苏 常州 214500)摘要:传统深松铲破坏土壤仿真模型存在土壤间作用力分析不够精准,无法确定深松铲作业的最佳入土深度。为此,构建基于光滑粒子流体动力学方法(Smoothed Particle Hydrodynamics,SPH)算法的深松铲破坏土壤仿真模型。将深松铲破坏土壤过程分为剪切失效和切削失效,采用光滑粒子流体动力学方法,通过差值函数描述全部质点存在的相互作用,获取最佳光滑长度值。引入有限元分析方法,利用点-面形式构建土壤分析的耦合模型,确定土壤失效准则。利用此模型,分别探究铲尖与铲柄破坏土壤、刀面破坏土壤时深松铲受到的阻力情况。根据前铲与后铲对土壤的破坏作用,计算前后两铲对土壤的破坏半径及土壤扰动体积,完成深松铲破坏土壤的仿真模型构建。仿真结果证明,利用上述模型能够得出入土深度对土壤破坏的轮廓变化规律,可确定最佳入土深度,有效提高深松铲作业效果。关键词:光滑粒子流体动力学算法;深松铲;土壤破坏;仿真模型中图分类号:S222 文献标识码:BSimulation Model of Soil Damage by Sub-soilingshovel Based on SPH AlgorithmZHANG Zhi-jun,GAO Yi-yu(School of Environmental&Safety Engineering,Changzhou University,Changzhou Jiangsu 214500,China)ABSTRACT:Traditionally,the simulation model cannot accurately analyze the forces between soils and determinethe best depth into soil.Therefore,a simulation model of soil destroyed by sub-soiling shovel was built based onsmoothed particle hydrodynamics(SPH)algorithm.The process of soil failure includes shear failure and cutting fail-ure.According to the smoothed particle hydrodynamics algorithm,the interaction between particles was described bythe difference function,so that the best smooth length was obtained.The finite element analysis was introduced.Andthe coupling model of soil analysis was constructed by combining point and plane to determine the criterion of soilfailure.Based on this model,we researched the resistance of sub-soiling shovel when the tine and handle of shoveldestroyed the soil and the tool face destroyed the soil.According to the damaging effect of front shovel and back shov-el on soil,we calculated the radius and volume of the soil destroyed by front shovel and back shovel.Finally,thesimulation model of the soil destroyed by sub-soiling shove was built.Simulation results prove that the above modelcan be used to obtain the contour change rule of soil damage caused by the depth of entering and leaving the soil,de-termine the best depth of entering the soil,and effectively improve the effect of subsoiling shovel operation.KEYWORDS:Smoothed particle hydrodynamics algorithm;Sub-soiling shovel;Soil failure;Simulation model1 引言土壤是影响农作物产量的关键因素。当前,我国某些地区土壤出现沙化现象,导致肥力下降,农作物减产。此外,长期使用化学除草方式,使土壤抗药性增强,生态环境恶化。近年来,土壤深松方法逐渐引入到农业生产中,土壤深松方法属于一种保护性耕作方法,不会改变植被与土粒结构。深松铲是完成该项技术的主要部件,它的主要能耗是克服土壤阻力做功。由于受阻力与功耗影响较大,还需进一步改进优化,同时深松作业涉及范围广,通过分析深松铲切削过程对土壤造成的影响,利用数值模拟方式获得土壤变化情况,为深松铲优化提供依据。为实现深松铲破坏土壤的仿真,马跃进1等人提出基于离散元法的深松铲减阻效果仿真分析。在深松铲顶部设置092减阻降耗的凸圆刃,利用非线性黏结弹性塑性接触模型构建土壤仿真模型,通过插件将颗粒和深松铲接触作用力导出,完成深松铲耕作阻力仿真。该研究可获取耕作设备切削阻力,但没有考虑到不同设备会产生不同阻力,导致该模型使用范围受到限制;丁启朔2等人提出基于田间摄像的深松扰动行为和效应研究。利用土壤耕作原位综合测试系统完成深松实验,从五个不同方位录制深松视频,通过视频分析土壤紧实度与失效机理以及深松效应影响。但是由于该方法模型的实验成本较高,导致该模型无法广泛应用。为此,本文利用 SPH 方法构建深松铲破坏土壤仿真模型。SPH 方法属于拉格朗日算法3,是解决高度非线性问题的有效方式。将此方法与有限元(Finite element method,FEM)相结合,最大程度体现出两种算法优势,更加精准地实现深松铲破坏土壤过程仿真。2 深松耕作模型构建在深松耕作时,梨底层被间隔式破坏,生成相间存在的构造,此种构造能够缓解土壤退化现象,有助于土壤保持肥力。还能提高气体交换能力,确保土壤充分吸收养分,使农作物增产。深松铲对土壤的破坏原理如图 1 所示。图 1 深松铲破坏原理示意图可将整个切削过程当作土壤受到外界作用力被划分成不同形状的过程。有关研究表明,入土深度不同土壤的失效方式也不同5。田地可当作没有边界的土壤,但受到仿真模型限制,无法建立较大的土壤模型。本文建立的土壤模型长宽高分别表示为 0.5m、0.3m、1.5m。结合有关约束条件仿真实际状况,利用 SPH 方法分析深松铲破坏土壤的过程。3 基于耦合分析的深松铲破坏土壤仿真模型3.1 SPH+FEM 耦合分析模型SPH 的基础方程为守恒方程6,能够有效解决大变形问题,适用于分析土壤变形破坏等非线性问题。差值理论是SPH 的根本理论,利用差值函数描述全部质点发生的相互作用,获取某点估计值,将守恒定律变换为积分形式,完成转换求和。在 SPH 算法中,质点近似函数公式如下(x)dx=f(y)W(x-y,h)dy(1)式中,W 代表核函数,h 表示光滑长度,该值随时间与空间的变化而变化。将核函数 W 通过辅助函数 表示,获取尖峰函数 W(x,h)W(x,h)=(x)h(x)d(2)式中,d 代表空间维数。该算法中普遍使用的光滑核表示为(u)=C 1-32u2+34u3|u|114(2-u)31|u|20|u|2|(3)式中,C 为常量,取决于空间维数多少。u 为光滑核系数。h 的不断变化会对仿真精度产生重要影响,它可以确保存在足够多的质点,使质点连续有效。若 h 为固定值,容易出现数值畸变现象,光滑长度的最佳变化范围规定为C1h0 h C2h0(4)式中,h0代表原始光滑长度,C1与 C2表示缩放因子。且 C1=C2,h 值始终不变。但是 SPH 算法在边缘约束处理方面还需改进,本文将 SPH 与 FEM 相结合,当变形区域较大时利用SPH 方法,如果变形区域较小则利用 FEM 方法,确保两种算法发挥出最大优势,提高仿真模型精度。图 2 耦合结构示意图图 2 代表两种算法的耦合模型示意图。上部分为 SPH质点,下侧区域是 FEM 网格。利用惩罚函数达到两个区域力学参数传递的目的。网格边界部分与 SPH 质点相互接触,采用点-面形式完成耦合。耦合算法下,粒子密度若比单元密度大,则表明仿真精度较高。为提高仿真模型精度,将 0.2 倍与 2.0 粒子间距分别当作光滑长度极小、极大值。因此失效准则7可定义为|fn|fn,fail|m1+fsfs,fail|m2 1(5)式中,fn与 fs分别代表法向力与剪切力,fn,fail与 fs,fail分别表示法向与剪切力指数。1923.2 深松铲破坏土壤仿真模型构建3.2.1 深松铲受力分析利用上述构建的耦合模型分析深松铲作业时,产生的阻力情况。1)铲尖和铲柄同时破坏土壤深松铲通常为平移前进,铲尖破坏土壤时会产生土壤耕作阻力。因为铲尖与土壤之间存在的深度较大,所以前进过程中受到的力来自土壤垂直方向。铲柄在作业时会受到土壤水平方向的耕作阻力,假设两种阻力分别表示为 k1与 k2,铲尖与铲柄的长度通过 l1、l2描述,则同时生成的阻力 F1的计算公式如下F1=k1l1+k2l2(6)2)刀面破坏土壤深松铲工作时,刀面对土壤产生挤压作用,土壤与刀面之间的力包括法向压力 FN与摩擦阻力 Ff=FN,这两种力在水平方向中的投影为阻力 F2F2=FNsin +Ffcos =FN(sin +cos)(7)式中,代表深松铲入土角度,为深松铲和土壤之间的摩擦系数8。已有研究显示,当耕地深度在 400-500mm 区间内,土壤的正向压力情况近似表示为 0=90kN/m2,因此土壤的正向压力9计算公式如下=0(8)式中,代表土壤有关系数。假定深松铲的工作区域面积表示为 S,法向压力 FN的计算公式如下FN=S0(9)则得出法相压力 FN与摩擦阻力 Ff存在如下三角函数关系Ff=FNtan =FN(10)式中,代表土壤内摩擦角。将式(9)与式(10)代入到式(7)中,重新获得刀面破坏土壤生成的阻力 F2的计算公式F2=S0(sin +tan cos)(11)3.2.2 深松铲对土壤的破坏过程获取深松铲与土壤之间的摩擦力后,对深松铲破坏土壤的过程完成仿真建模。破坏过程可分为前铲对土壤的切削与破碎、后铲对土壤的切削与破碎两个阶段。1)前铲破坏过程前铲的入土深度小,先与土壤发生接触,如果铲尖表层光滑,没有磨损,这时铲尖作用在无限大土体上,形成作用力P 及夹角。在此过程中,土壤受到的压缩力 r呈现径向分布,而垂直方向上表现为自重应力 sz。土壤产生的抗剪性能持续增大,当出现最大值,铲尖对土壤生成削切作用,整个土块失效破坏。发生形变的土块变换成有限体积,这时仅受到来自没有失效土壤的作用力10。前铲逐渐前行,铲尖使部分土壤提升。在持续作用力下导致上升的土壤与土体之间形成作用应力。在不同力的相互作用下,升起的土壤表层与内部出现褶皱现象,生成拉伸应