基于SPH算法的深松铲破坏土壤仿真模型_张志军.pdf
《基于SPH算法的深松铲破坏土壤仿真模型_张志军.pdf》由会员分享,可在线阅读,更多相关《基于SPH算法的深松铲破坏土壤仿真模型_张志军.pdf(5页珍藏版)》请在咨信网上搜索。
1、收稿日期:2021-04-08 修回日期:2021-05-28 第 40 卷 第 4 期计 算 机 仿 真2023 年 4 月 文章编号:1006-9348(2023)04-0290-05基于 SPH 算法的深松铲破坏土壤仿真模型张志军,高奕珏(常州大学环境与安全工程学院,江苏 常州 214500)摘要:传统深松铲破坏土壤仿真模型存在土壤间作用力分析不够精准,无法确定深松铲作业的最佳入土深度。为此,构建基于光滑粒子流体动力学方法(Smoothed Particle Hydrodynamics,SPH)算法的深松铲破坏土壤仿真模型。将深松铲破坏土壤过程分为剪切失效和切削失效,采用光滑粒子流体动力
2、学方法,通过差值函数描述全部质点存在的相互作用,获取最佳光滑长度值。引入有限元分析方法,利用点-面形式构建土壤分析的耦合模型,确定土壤失效准则。利用此模型,分别探究铲尖与铲柄破坏土壤、刀面破坏土壤时深松铲受到的阻力情况。根据前铲与后铲对土壤的破坏作用,计算前后两铲对土壤的破坏半径及土壤扰动体积,完成深松铲破坏土壤的仿真模型构建。仿真结果证明,利用上述模型能够得出入土深度对土壤破坏的轮廓变化规律,可确定最佳入土深度,有效提高深松铲作业效果。关键词:光滑粒子流体动力学算法;深松铲;土壤破坏;仿真模型中图分类号:S222 文献标识码:BSimulation Model of Soil Damage
3、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.There
4、fore,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 des
5、cribed 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 su
6、b-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 o
7、f 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 subs
8、oiling shovel operation.KEYWORDS:Smoothed particle hydrodynamics algorithm;Sub-soiling shovel;Soil failure;Simulation model1 引言土壤是影响农作物产量的关键因素。当前,我国某些地区土壤出现沙化现象,导致肥力下降,农作物减产。此外,长期使用化学除草方式,使土壤抗药性增强,生态环境恶化。近年来,土壤深松方法逐渐引入到农业生产中,土壤深松方法属于一种保护性耕作方法,不会改变植被与土粒结构。深松铲是完成该项技术的主要部件,它的主要能耗是克服土壤阻力做功。由于受阻力与功耗影响较大,
9、还需进一步改进优化,同时深松作业涉及范围广,通过分析深松铲切削过程对土壤造成的影响,利用数值模拟方式获得土壤变化情况,为深松铲优化提供依据。为实现深松铲破坏土壤的仿真,马跃进1等人提出基于离散元法的深松铲减阻效果仿真分析。在深松铲顶部设置092减阻降耗的凸圆刃,利用非线性黏结弹性塑性接触模型构建土壤仿真模型,通过插件将颗粒和深松铲接触作用力导出,完成深松铲耕作阻力仿真。该研究可获取耕作设备切削阻力,但没有考虑到不同设备会产生不同阻力,导致该模型使用范围受到限制;丁启朔2等人提出基于田间摄像的深松扰动行为和效应研究。利用土壤耕作原位综合测试系统完成深松实验,从五个不同方位录制深松视频,通过视频分
10、析土壤紧实度与失效机理以及深松效应影响。但是由于该方法模型的实验成本较高,导致该模型无法广泛应用。为此,本文利用 SPH 方法构建深松铲破坏土壤仿真模型。SPH 方法属于拉格朗日算法3,是解决高度非线性问题的有效方式。将此方法与有限元(Finite element method,FEM)相结合,最大程度体现出两种算法优势,更加精准地实现深松铲破坏土壤过程仿真。2 深松耕作模型构建在深松耕作时,梨底层被间隔式破坏,生成相间存在的构造,此种构造能够缓解土壤退化现象,有助于土壤保持肥力。还能提高气体交换能力,确保土壤充分吸收养分,使农作物增产。深松铲对土壤的破坏原理如图 1 所示。图 1 深松铲破坏
11、原理示意图可将整个切削过程当作土壤受到外界作用力被划分成不同形状的过程。有关研究表明,入土深度不同土壤的失效方式也不同5。田地可当作没有边界的土壤,但受到仿真模型限制,无法建立较大的土壤模型。本文建立的土壤模型长宽高分别表示为 0.5m、0.3m、1.5m。结合有关约束条件仿真实际状况,利用 SPH 方法分析深松铲破坏土壤的过程。3 基于耦合分析的深松铲破坏土壤仿真模型3.1 SPH+FEM 耦合分析模型SPH 的基础方程为守恒方程6,能够有效解决大变形问题,适用于分析土壤变形破坏等非线性问题。差值理论是SPH 的根本理论,利用差值函数描述全部质点发生的相互作用,获取某点估计值,将守恒定律变换
12、为积分形式,完成转换求和。在 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 为固定值,容易出现数值畸变现象,光
13、滑长度的最佳变化范围规定为C1h0 h C2h0(4)式中,h0代表原始光滑长度,C1与 C2表示缩放因子。且 C1=C2,h 值始终不变。但是 SPH 算法在边缘约束处理方面还需改进,本文将 SPH 与 FEM 相结合,当变形区域较大时利用SPH 方法,如果变形区域较小则利用 FEM 方法,确保两种算法发挥出最大优势,提高仿真模型精度。图 2 耦合结构示意图图 2 代表两种算法的耦合模型示意图。上部分为 SPH质点,下侧区域是 FEM 网格。利用惩罚函数达到两个区域力学参数传递的目的。网格边界部分与 SPH 质点相互接触,采用点-面形式完成耦合。耦合算法下,粒子密度若比单元密度大,则表明仿真
14、精度较高。为提高仿真模型精度,将 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)铲尖和铲柄同时破坏土壤深松铲通常为平移前进,铲尖破坏土壤时会产生土壤耕作阻力。因为铲尖与土壤之间存在的深度较大,所以前进过程中受到的力来自土壤垂直方向。铲柄在作业时会受到土壤水平方向的耕作阻力,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 SPH 算法 深松铲 破坏 土壤 仿真 模型 张志军
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。