基于LS―DYNA深松铲阻力仿真.docVIP

基于LS―DYNA深松铲阻力仿真 　　摘要： 为降低耕作阻力，揭示深松铲与土壤之间的关系特性，根据深松铲切削土壤的工作特点，利用LS-DYNA模拟分析深松铲切削土壤的过程，获得切削土壤的应力变化规律。结果表明，当深松铲以0.49 m/s的初速度切削、深松深度为250 mm时，单个深松铲的最大切削阻力为2 230 N，土壤在1.75 s发生崩裂，达到深松效果。通过试验测试，单个深松铲受力约为2 332.5 N，与仿真误差为4.5%，验证了仿真的合理性，仿真数据结果可以为实际深松过程提供参考。 　　关键词： 土壤；切削；有限元法；深松铲；LS-DYNA模拟分析；阻力仿真 　　中图分类号： S222.19 文献标志码： A 　　文章编号：1002-1302（2017）22-0257-03 　　深松铲是深松的核心部件，通过深松可以改善土壤三相比，从而促进植物的生长[1]。然而，深松过程中有特别大的阻力。因此，须要探明在深松土壤切削过程中阻力的变化规律，达到减少阻力的目的。 　　在以往关于土壤切削的研究中，主要采用试验方法，但试验具有受外界环境因素影响大且经费多、周期长等缺点[2]。近50年来，随着数值分析模拟技术的发展，为研究土壤切削提供了一种新的方法。数学模拟切削土壤方法主要包括传统分析方法、有限单元法（finite element method，简称FEM）、离散单元法（distinct element method，简称DEM）[3]。近年来，多利用有限元法研究土壤切削问题，通过模拟仿真，可以方便清楚地观察到土壤和工作部件的情况，例如土壤的运移状况、破坏位置和工件部件各部分的受力状态。这些模拟仿真数据为耕作部件的设计优化提供了理论依据。 　　针对土壤非线性不均匀的特性，本研究采用动力分析方法探讨深松铲对土壤的切削过程，通过显式非线性动力分析软件ANSYS/LS_DYNA来进行深松铲切削土壤的数值模拟，建立深松铲深松工作过程的有限元仿真模型，对深松过程进行三维数值仿真研究。 　　1 有限元建模 　　1.1 土壤模型 　　ANSYS/LS-DYNA 12.1程序材料库中193号材料（MAT_DRUCKER_PRAGER）模型是一种弹塑性本构模型。由于这种材料本构模型能准确地体现土壤特性，简单实用而且能够得到准确的结果，所以被广泛地用于模拟岩土材料中，具体参数取值见表1[4]。Mohr-Coulomb屈服条件通常被使用在土壤切削仿真中，但其有不可忽略的缺点：它的屈服面在平面上表现为特殊的等边六边形（具有棱角奇异性的不等角），这大大增加了数值计算的困难。与Mohr-Coulomb屈服条件相比，Drucker-Prager屈服准则不仅计算简单，而且在选择恰当材料常数的前提下，就可以寻找到与之匹配的Coulomb模型。基于以上比较，本研究采用改进扩展后的Drucker-Prager准则[5]，即： 　　式中：C为偏应力第3应力不变量；I为单轴拉伸屈服应力与单轴压缩屈服应力比值；F1为凝聚力；为材料摩擦角；F2为Mises等效应力；F3为平均压应力；σ1、σ2、σ3为3个方向的主应力；t为应力。 　　选定土壤的实体模型为800×800×800的立方体，由于模型比较简单，所以直接在ANSYS中生成。土壤选用solid164单元，这种单元常用于三维显示结构实体，由8个节点构成而且支持所有许可的非线性特性，对土壤采用智能网格划分，选择土壤模型的12条边，设置网格数为20，单元形状采用六面体8个节点，划分为9 261个节点，8 000个单元[6]。 　　1.2 深松铲模型 　　ANSYS在构建复杂实体模型方面比较薄弱，所以本研究利用SolidWorks构建深松铲模型并通过接口导入到ANSYS中[7]。深松铲采用 65Mn材料，其密度为7.8×103 kg/m3，弹性模量为2.10万MPa，泊松比为0.3。由于深松铲是一个整体，所以采用计算精度较高的四面体单元solid45单元类型，采用智能自由网格划分方法对深松铲模型进行网格划分，设定尺寸级别为5级，为减小应力集中的影响，在深松铲铲刃部分进行网格细化，深松铲划分为12 724个节点， 19 046 个单元。在建立土壤-深松铲数值模拟模型时，对模型进行网格划分后如图1所示。 　　1.3 创建单元 　　为了给后面定义接触类型和施加边界条件作基准，要先创建单元（简称PART）。PARTS是1个单元集合体，它是由相同单元类型、材料号和实常数组成的。若要生成PARTS，必须要求选择所有实体。本案例生成的PARTS结果如图2所示，其中1代表深松铲模型，2代表土壤模型。网格划分数量相当是能够计算的必要条件，从图2可以看出，深松铲模型和土壤模型网格数量相当，这为后续计算提供了可能。