一种分析车轮沙石路面行驶性能的离散元路面更替法
【专利摘要】本发明提供一种分析车轮沙石路面行驶性能的离散元路面更替法，其包括：在给定空间范围内生成离散单元集，并使其在自重作用下达到稳定状态，然后将该离散单元集备份，该离散单元集即为路面样本块；将备份好的路面样本块复制两份，按照前后顺序组成初始路面，初始路面长度等于路面样本块的2倍；令车轮初始位置在第一个路面样本块的中间，且由第一个路面样本块向第二个路面样本块行驶，当车轮行驶距离等于路面样本块长度时，将第一个路面样本块删除，同时再在车轮行驶方向的前面铺设同样的路面样本块。本发明可以实现车轮行驶性能分析中沙石路面的离散元模拟，且将模型规模始终控制在个离散单元内，从而使离散元分析具有工程应用价值。
【专利说明】一种分析车轮沙石路面行驶性能的离散元路面更替法【技术领域】
[0001]本发明涉及一种大规模沙石路面的离散元分析法，具体为一种分析车轮沙石路面行驶性能的离散元路面更替法。
【背景技术】
[0002]开发沙漠中的丰富资源，探索月球、火星等未知领域，复杂地形下的军事作战，日益紧迫的风沙治理及效率低下的农业生产等都涉及到车轮与松散路面间的相互作用。沙石路面是典型的松散路面，沙石颗粒的非连续特性及其自身的转动、大位移运动等都会导致路面破坏，从而对车轮的行驶性能产生很大影响。因此，深入地研究车轮与沙石类松软路面间的相互作用原理，对加速提高我国越野车辆开发水平，促进民用工业、国防事业、农业机械化等的发展具有重要意义。
[0003]在该领域的研究方法中，传统的实验方法存在周期长、成本高、过程复杂等缺点，而且在观察沙石颗粒的流动趋势等方面存在一定的局限性，不能有效地分析沙石颗粒间的细观参数对车轮行驶行的影响。近年来，随着计算机技术的飞速发展，仿真方法被越来越多的应用于该领域的研究。有限元方法作为一种传统的仿真方法，被广泛应用于该领域的研究，例如米用粘弹性材料(T.Hiroma, S.ffanjii, T.Kataoka, etal.(1998) Stressanalysis using FEM on stress distribution under a wheel considering frictionwith adhesion between a wheel and soil.Journal of Terramechanics, 34(4):225-233)或弹塑性材料(J.P.Hambleton, A.Drescher.(2009)Modeling wheel-1nducedrutting in soils:rolling.Journal of Terramechanics, 46(2):35-47 和 KM.Xia, YMYang.(2012) Three-dimensional finite element modeling of tire/groundinteraction.1nternational Journal for Numerical and Analytical Methods inGeomechanics, 36(4):498-516等)模拟松软沙石路面,分析车轮行驶过程的下陷和路面的变形。但是有限元方法不能准确地描述沙石路面的非连续介质特性和沙石颗粒流动性，同时也不能合理模拟复杂的轮胎花纹与沙石颗粒间的相互作用。
[0004]Cundall等人于1971提出的离散元方法适用于研究颗粒介质的力学行为(PA.Cundall.(1971)A computer model for simulating progressive large movements inblocky systems.S-ymposium ISRM, Proc, 2:129-136)。该方法能够真实地再现车轮行驶过程中，沙石路面的破坏及颗粒飞溅、流动等现象，可以有效弥补传统有限元方法的不足，准确地描述沙石路面的非连续介质特性，是地面力学分析中最有前景的分析工具。但是，其在获得散体单元间的作用力时，效率随着单元数量规模的增大急剧降低。
[0005] 然而在实际工程问题中，测试车轮沙石路面行驶性能时，通常需要行驶较长距离，如果对该长度范围内的沙石颗粒采用离散单元逐一建模，则分析时间过长，仿真分析也失去实际工程意义。
【发明内容】
[0006]本发明提出一种分析车轮沙石路面行驶性能的离散元路面更替法，通过控制分析过程中离散单元的数量，实现对大规模沙石路面的离散元模拟，具体技术方案如下。
[0007]—种分析车轮沙石路面行驶性能的离散元路面更替法，其包括:
[0008]采用离散单元对沙石路面进行描述:在给定空间范围内生成离散单元集并使这些离散单元在自重作用下压实达到稳定状态，然后将该离散单元集中的离散单元编号、坐标、半径、速度和加速度信息存储、备份，该压实后的离散单元集即为路面样本块；
[0009]由两份所述路面样本块排列组成初始路面:将备份好的路面样本块复制两份，按照前后顺序组成初始路面，初始路面长度等于路面样本块的2倍；
[0010]令车轮初始位置在第一个路面样本块的中间，且车轮由第一个路面样本块向第二个路面样本块行驶，当车轮行驶距离等于路面样本块长度时，将第一个路面样本块删除，同时再在车轮行驶方向的前面铺设同样的路面样本块。
[0011]上述分析车轮沙石路面行驶性能的离散元路面更替法中，路面样本块及初始路面的生成方法为:首先在给定区域内生成离散单元集，该离散单元集中的离散单元坐标在该区域内随机分布，离散单元半径在给定半径范围内随机取值；其次，在该区域周围及底部添加刚性墙约束，并添加重力场使离散单元在自重下达到稳定状态，存储离散单元的编号、坐标、半径、速度和加速度信息，将达到稳定状态后的离散单元集备份、存储为B ;然后，将离散单元集B复制两份，分别记为Atl与A1,并将离散单元集Atl与A1按前后顺序排列，即组成初始路面。
[0012]上述分析车轮沙石路面行驶性能的离散元路面更替法中，初始路面的区域范围和单元数目等于离散单元集B的区域范围及单元数目的2倍。
[0013]上述分析车轮沙石路面行驶性能的离散元路面更替法中，初始路面组成后，在初始路面区域的两侧和底面添加无限长的刚性墙约束，以约束车轮行驶过程中离散单元在这些方向的运动。
[0014]上述分析车轮沙石路面行驶性能的离散元路面更替法中，当车轮行驶距离等于路面样本长度，即车轮到达离散单元集A1所在区域范围的中间位置时，将离散单元集Atl删除，同时，再次将离散单元集B复制一份，记为A2,并将离散单元集A2添加在A1前端，则A1和A2组成第一次更替后的路面，在车轮行驶方向不断重复此过程，实现离散元路面的更替。
[0015]上述分析车轮沙石路面行驶性能的离散元路面更替法中，路面更替完成后，沙石路面始终有两份路面样本块组成，且其区域范围和单元数目仍等于离散单元集B的区域范围和单元数目的2倍。
[0016]与现有技术相比，本发明具有如下优点和技术效果:本发明可以实现车轮行驶性能分析中沙石路面的离散元模拟，且将模型规模始终控制在2N个离散单元内，从而使离散元分析具有工程应用价值。该发明的成功实施对研究越野车辆沙石路面行驶性能具有重要意义。
【专利附图】
【附图说明】
[0017]图1表示执行本发明的装置结构图。
[0018]图2a表示本发明的执行流程图。
[0019]图2b表示流程图中的组成初始路面过程C的详细流程图。[0020]图2c表示流程图中的离散元路面更替过程G的详细流程图。
[0021]图3表示初始离散单元集自重压实后的状态及其所在的区域范围。
[0022]图4表示由离散单元集Aci和A1组成的初始路面。
[0023]图5将车轮置于离散单元集Atl所在区域的中间位置。
[0024]图6表示车轮行驶至离散单元集A1所在区域中点时Atl中的离散单元对车轮行驶行为的影响可以忽略且离散单元集A1的前端没有发生破坏，可以执行路面更替。
[0025]图7表示离散元路面的更替过程，即删除O?375_范围内的离散单元集Atl,同时在750?1125mm区域添加离散单元集A2。
【具体实施方式】
[0026]以下结合相关二维示意图对本发明的实施作进一步说明，但本发明实施不限于此。
[0027]执行本发明的装置结构图如图1所示，包括输入设备、输出设备、内存设备、外部存储设备及中央处理器。其中，输入设备如鼠标、键盘等，用于输入初始模型；输出设备用于显示分析结果，例如显示器、打印机等；内存设备(RAM/R0M)主要是用于暂时存储分析结果便于以后分析时的调用；外存储设备主要是将输出结果保存便于后处理中进行数据分析；中央处理器(CPU)集中处理数据读入、分析、更新等，执行本发明。
[0028]本发明的执行流程图如图2a所示，主要包含以下步骤:初始离散单元模型信息的读入；离散单元的自重压实及备份；调用备份离散单元集并复制两份，按顺序排列组成初始路面(图2b);将车轮放置于初始路面上；给车轮加载角速度ω，使其沿初始路面排列方向行驶；满足更替条件时执行离散元路面更替(图2c);输出并存储结果。
[0029]首先建立离散单元集的初始模型:在给定空间区域内(空间范围为:X方向O?a，Y方向O?b, Z方向O?c,其中Z方向为竖直方向，a, b, c的大小根据车轮尺寸和路面特性等具体参数确定)随机生成给定半径范围的离散单元集，单元数目记为N，单元编号为I?N。随后，在该区域四周边界及底部边界添加刚性墙约束。接着，添加重力场，使这些离散单元在自身重力作用下达到稳定状态。存储达到稳定状态后的离散单元的编号、坐标、半径、速度和加速度信息，将此刻的离散单元集记为B，也即路面样本块，其空间区域范围变为O?a, Y:0?b,Z:0?c’，其中C，< C。自重压实过程中离散单元间的接触判断采用C_grid方法(R.John, EP Williams, B.Cook.(2004) A contact algorithm for partitioning Narbitrary sized objects.Engineering computations, 21 (2/3/4):235_348)，离散单兀接触判断搜索范围为X:0?a，Y:0?b，Z:0?c，即自重压实之前离散单元集所处的区域。
[0030]其次，组成初始路面:将离散单元集B复制两份，分别记为Atl和4。然后，沿X方向将离散单元集Ap A1按前后顺序排列，组成初始路面，具体实现方法为:离散单元集Atl保持不变，将离散单元集A1中的所有离散单元的X坐标值都增加a、离散单元的编号修改为N+1?2N。最终生成的初始路面的大小为2aXbXc’，区域范围为X:0?2a，Y:0?b，Z:O?c’，离散单元数目为2N，离散单元的编号为I?2N。
[0031]然后，将车轮置于离散单元集Atl所在区域的中间位置，给车轮加载角速度ω，使车轮沿X正方向(Atl指向A1方向)行驶。同时，修改离散单元接触判断的搜索范围(改为X:O?2a，Y:0?b，Z:0?c’)，此时模型中离散单元数目为2N。[0032]离散元路面更替:当车轮行驶至离散单元集A1所在区域的中点位置时，将对车轮行驶性能没有影响的离散单元集Atl (其包含的离散单元的单元编号为I?N)删除，同时，在A1的前端添加一个新的离散单元集A2,由离散单元集A1和A2组成新的沙石路面，完成第一次路面更替。实现方法为:将离散单元集Atl中的离散单元信息置零(即删除单元)，同时，将离散单元集B复制一份，记为A2,并将A2中的离散单元的X坐标值增加2a，单元编号记为I?N，实现A” A2在X方向的顺序排列。此时，沙石路面的大小仍为2aXbXc’，且单元总数仍为2N，路面区域范围更改为(X:a?3a，Y:0?b，Z:0?c’)。随后，将离散单元接触判断的搜索范围修改为(X:a?3a，Y:0?b，Z:0?c’)。沿车轮行驶方向重复以上过程，实现车轮沙石路面行驶性能分析时的离散元路面更替。
[0033]以下再结合具体的二维实例对本发明的具体实施步骤作说明:
[0034]1、首先在区域范围X:0?375mm，Y:0?60mm内，随机生成半径范围为I?2mm的离散单元集，单元数目为1350个，存储于外部存储装置中。
[0035]2、将上述离散单元集作为初始输入，通过输入设备读取后，在离散单元集所在区域的四周及底部加刚性墙约束，同时添加重力场，使离散单元在自重作用下达到稳定状态，将此刻的离散单元的坐标、编号、半径、速度和加速度信息存储于内部存储设备中并备份，将稳定后的离散单元集记为B，如图3所示。离散单元的编号为I?1350，接触判断搜索范围为X:0?375臟，Y:0?55mm。
[0036]3、读取步骤2中备份的离散单元集B，并复制两份，其中，第一份保持不变，记为A0 ;将第二份离散单元集内的离散单元X坐标都增加375_，并将其离散单元的编号修改为1351?2700，并记为Ap此时，Atl与A1沿X方向按照前后顺序组成初始路面，结果如图4所示。其中离散单元集Aci的区域范围为X:0?375mm, Y:0?55mm,所包含的离散单元的单元编号为I?1350 ;离散单元集A1的区域范围为X:375?750mm，Y:0?55mm，所包含的离散单元的单元编号修改为1351?2700。则初始路面的区域范围为X:0?750mm，Y:0?55_，包含的离散单元总数为2700，单元编号为I?2700，将这些单元信息存储于内部存储设备中。
[0037]4、将直径为120mm的车轮放置于离散单元集Atl所在区域的中点位置，车轮的初始位置如图5所示，车轮质心坐标O为(187.5，115)。给车轮加载角速度ω =0.5rad/s，车轮开始沿X轴正方向行驶。在此过程中，离散单元接触判断搜索范围为X:0?750mm，Y:0?5 5mm ο
[0038]5、当车轮沿X方向行驶375_，车轮质心O到达离散单元集A1所在区域的中间位置，如图6所示(图中O’为初始时刻车轮质心位置)。此时，离散单元集Atl中的离散单元对车轮行驶性的能影响可以忽略。同时，离散单元集A1中的单元，其左端(图6中X坐标为750_处)没有发生明显的大变形破坏。此时，进行离散元路面更替:删除离散单元集Atl中的离散单元，同时，将离散单元集B复制一份，并将其中的离散单元X坐标增加750_，单元编号记为I?1350，然后将该离散单元集记为A2。此时离散单元集A1与A2组成新的沙石路面，完成第一次路面更替，结果如图7所示。其中，新的沙石路面所在的区域范围变为X:375?1125mm，Y:0?55mm，单元总数仍为2700。
[0039]沿车轮行驶方向重复以上步骤，不断更新接触判断搜索区域，在区域范围X:k*375?(k+2)*375，Y:0?55mm内，对2700个离散单元进行求解(k为更替次数)，并将分析结果输出、存储于外部存储装置中。从而实现对车轮在沙石路面行驶性能的仿真分析。每次更替时，新添加的离散单元集中的离散单元的编号及单元坐标更新方式为:当k为奇数时，离散单元编号取I~1350，任意单元i的坐标(夂1A4 I的更新方法为:
【权利要求】
1.一种分析车轮沙石路面行驶性能的离散元路面更替法，其特征在于包括: 采用离散单元对沙石路面进行描述:在给定空间范围内生成离散单元集并使这些离散单元在自重作用下压实达到稳定状态，然后将该离散单元集中的离散单元编号、坐标、半径、速度和加速度信息存储、备份，该压实后的离散单元集即为路面样本块； 由两份所述路面样本块排列组成初始路面:将备份好的路面样本块复制两份，按照前后顺序组成初始路面，初始路面长度等于路面样本块的2倍； 令车轮初始位置在第一个路面样本块的中间，且车轮由第一个路面样本块向第二个路面样本块行驶，当车轮行驶距离等于路面样本块长度时，将第一个路面样本块删除，同时再在车轮行驶方向的前面铺设同样的路面样本块。
2.根据权利要求1所述的一种分析车轮沙石路面行驶性能的离散元路面更替法，其特征在于，路面样本块及初始路面的生成方法为:首先在给定区域内生成离散单元集，该离散单元集中的离散单元坐标在该区域内随机分布，离散单元半径在给定半径范围内随机取值；其次，在该区域周围及底部添加刚性墙约束，并添加重力场使离散单元在自重下达到稳定状态，存储离散单元的编号、坐标、半径、速度和加速度信息，将达到稳定状态后的离散单元集备份、存储为B ;然后，将离散单元集B复制两份，分别记为Atl与A1,并将离散单元集Atl与A1按前后顺序排列，即组成初始路面。
3.根据权利要求1或2所述的一种分析车轮沙石路面行驶性能的离散元路面更替法，其特征在于初始路面的区域范围和单元数目等于离散单元集B的区域范围及单元数目的2倍。
4.根据权利要求1或2所述的一种分析车轮沙石路面行驶性能的离散元路面更替法，其特征在于初始路面组成后，在初始路面区域的两侧和底面添加无限长的刚性墙约束，以约束车轮行驶过程中离散单元在这些方向的运动。
5.根据权利要求1或2所述的一种分析车轮沙石路面行驶性能的离散元路面更替法，其特征在于当车轮行驶距离等于路面样本长度，即车轮到达离散单元集A1所在区域范围的中间位置时，将离散单元集Atl删除，同时，再次将离散单元集B复制一份，记为A2,并将离散单元集A2添加在A1前端，则A1和A2组成第一次更替后的路面，在车轮行驶方向不断重复此过程，实现离散元路面的更替。
6.根据权利要求1或5所述的一种分析车轮沙石路面行驶性能的离散元路面更替法，其特征在于路面更替完成后，沙石路面始终由两份路面样本块组成，且其区域范围和单元数目仍等于离散单元集B的区域范围和单元数目的2倍。
【文档编号】G06F17/50GK103955560SQ201410134643
【公开日】2014年7月30日 申请日期:2014年4月3日 优先权日:2014年4月3日
【发明者】臧孟炎, 赵春来 申请人:华南理工大学