基于LF1352森林运载车的履带板仿生设计及通过性

刘九庆， 李新源， 孙术发

（东北林业大学 a. 机电工程学院， b. 工程技术学院， 哈尔滨 150040）

1 仿北极熊熊爪结构的履带设计

1.1 履带板形状设计

自然界很多动物都拥有在复杂道路上行走的能力，这得益于它们特殊的足部结构[11].北极熊是一种生活在北极的世界最大食肉动物之一.北极熊拥有在冰雪路面上快速奔跑的能力，这也源于它们特殊的足部结构.北极熊的熊掌宽超过25 cm，爪子长度可达10 cm.当北极熊在冰雪路面上行走时，尖锐的脚趾会插入坚硬的积雪中，积雪在挤压区域内聚集，有利于增大北极熊脚掌的牵引力，北极熊脚印如图1所示.

图1 北极熊脚印图Fig.1 Image of polar bear footprint

通过对北极熊的熊爪进行研究，对履带板形态进行仿生设计，履带板花纹的三角形区域和两侧的啮合部分相当于北极熊熊爪的脚趾区域，使履带板轻易地扎入积雪来增大附着力.履刺中间的四边形空腔相当于北极熊的脚掌挤压部分，这种履带板在运动过程中可以限制积雪的流动，减少积雪流失，增大履带车辆在积雪路面的通过性.履带板通过销轴进行链接，履刺花纹如图2所示.

图2 履刺花纹Fig.2 Track grouser pattern

履带板节距计算公式为

（1）

式中，m为履带车的质量.履带板宽度为

B=（0.2～0.4）L

（2）

式中，L为履带板整体长度.

1.2 履带板驱动轮啮合部位设计

将履带板和驱动轮啮合部分设计在履带板的两侧，履带车在行进过程中，驱动轮可以限制履带板的水平自由度.啮合部分的高度过大会对车体产生影响，因此啮合部位的高度等同于履刺的深度，以增大履带板受到的剪切力.履带板啮合部分尺寸如图3所示（单位：mm）.

图3 挤压区域及整体尺寸简化示意图Fig.3 Simplified schematic diagram of extrusion area and overall size

2 力学计算

履带车辆在行驶过程中，履带板受到重力的作用发生下陷，履刺完全扎入积雪中，并发生水平方向位移.履带车辆在行驶过程中的牵引能力主要由履带板受到的牵引力决定，牵引力和土壤剪切力是一对相互作用力.通过研究履带车辆在行驶过程中受到的土壤剪切力可以得出履带车辆的通过性，积雪通常情况下被认为是一种特殊的土壤.履带车辆受到的土壤剪切力一方面取决于车辆的设计参数，另一方面取决于土壤的力学性能.由于不同土壤的力学性能差异较大，所以需要设计适用于积雪路面的履带板机构.由文献[12]及[13]可知，履带板在行驶过程中所受到的水平方向剪切力主要分为4个部分，分别为作用于履带板上的力F1，作用于履刺底部的力F2，作用于履刺和啮合部位上的力F3，作用于履带板侧面的力F4.履带板的受力如图4所示，履带板履刺结构参数如图5所示.

图4 履带板受力分析Fig.4 Force analysis of track shoes

图5 履刺结构参数Fig.5 Track grouser structure parameters

由Bekker沉积公式可知

（3）

式中：p为法向应力；kc为土壤粘聚变形模量；kφ为土壤内摩擦变形模量；z为土壤压缩变形量；n为土壤变形指数.当履带车辆在路面上行驶时，可以把Bekker沉积公式简化为

p=Kzn

（4）

K可由式（5）推出，即

W=KznS1+K（z+h）nS2

（5）

式中：S1为履带板底部的面积；S2为履刺底部的面积.作用于履带板底部和履刺底部的水平剪切力F1和F2计算公式为

（6）

（7）

式中：C为土壤粘聚力；φ为土壤内摩擦角；η为履带板剪切速度系数，假设取1.

当履带板运动时，履刺受到车辆的重力作用完全压入土里.车辆的垂直载荷由履刺和履带板共同承担.

当车辆行驶时，作用在履刺上的法向压力主要分为两个部分，一部分是由土的自重引起的法向压力，另一部分是由粘聚力C所导致的土压力.土压力随着深度的变化不断增加，作用在履刺上的土压力如图6所示.

图6 均布载荷作用下土壤破坏模型Fig.6 Soil failure model under uniform load

考虑到履带板有一定的接地压力，使得土壤具有一定的均布载荷.土壤作用在履刺上的水平推力为

（8）

式中：γ为土壤容重；κp为被动土压力系数；P为均布载荷；I为啮合区域长度.

土壤对倒三角倾斜部分的额外推力为

F32=chb3tanα

（9）

式中，α为履刺的拔摸角度.则作用在履带板上的水平推力F3为

F3=F31+14F32

（10）

作用在履刺和履带板两侧的力F4可以通过弹性力学理论进行计算，即

（11）

因此，履带板在行驶过程中所受到的合力为

FH=F1+F2+F3+F4

（12）

以LF1352森林运载车为例，表1为履带板结构参数和土壤力学性能表.

3 积雪和履带板的仿真实验分析

当履带板受到垂直载荷时，积雪处于弹性阶段，考虑到本文中的履带车辆载荷为26 kPa，雪的变形分为两个区域：第一个区域为弹性区域，随着压力逐渐接近雪的粘聚力，积雪发生弹性形变，该区域通常来说比较小；第二个区域为材料塑性变形区，如果雪有足够的深度，接地压力会一直传递到雪的底部.在计算过程中，积雪被假设为连续体，把积雪视为一种固体材料.为了研究履带板在雪的剪切过程中受到的剪切力，运用有限元软件建立积雪弹塑性本构模型.

表1 履带板结构参数和土壤力学性能Tab.1 Structure parameters of track shoes and mechanical properties of soil

3.1 单条履带板仿真分析

在仿真过程中将履带板的运动简化为沉陷和水平剪切两个部分.当履带板运动时，在履带板顶部添加载荷来模拟接地比压，再给履带板施加一个沿Y轴正方向的水平位移来模拟履带板的水平剪切运动.由于履带板在雪中运动是一个非线性、大位移的过程，可通过动态显式分析方法对其进行分析.动态显式方法通过在分析过程中不断重新划分网格来求解，可以不断改变网格质量，有效解决分析过程中的网格畸变问题.考虑到履带板的材料属性为钢，在运动过程中基本不变形，为了减少计算时间，将履带板设置为刚体，并对履带板与雪模型接触的部分进行网格加密以保证计算数值的准确性.履带板运动分析过程中将履带板模型简化为宽度为570 mm、履带板节距为155 mm、履刺高度为30 mm、拔模角度为60°的履带板，如图7所示.

图7 履带板模型Fig.7 Track shoes model

将雪的三维模型简化为弹性模型和Mohr-Coulomb塑性模型.雪模型的面积超过履带板两倍尺寸外的部分基本没有变形，因此积雪模型尺寸为长1 500 mm、宽1 500 mm、高500 mm，相关参数见表1.将履带板放置在雪模型上，限制雪四周的自由度，在第一个分析步给履带板添加接地比压让履带沉陷，在第二个分析步给履带添加Y轴正方向50 mm位移，让履带板在积雪模型上水平剪切，加密模型如图8所示.

图8 仿真实验网格图Fig.8 Simulation experiment meshing

输出变量是履带板在运动过程中所受到的水平方向反作用力之和，即履带板受到的剪切力.在仿真结束时，履带板前方的积雪不断堆积在一起，正下方形成一个空腔，如图9所示，与实际情况基本一致.

图9 积雪受力图Fig.9 Force diagram of snow

为了计算车辆在积雪路面上行驶的牵引力，选用压实积雪的力学性能参数，计算出单块履带板受到的剪切力计算值约为3 600 N.

为了验证有限元分析的可靠性，将仿真结果与公式计算结果进行对比分析，结果如图10所示.履带板受到的剪切力在较小位移时快速上升到达顶点，仿真剪切力要稍小于计算剪切力，可以认定仿真有效.

图10 仿真数值与计算结果对比Fig.10 Comparison between simulation values and calculation results

3.2 多条履带板仿真分析

为了研究履带板受到的剪切力与履带板花纹数量的关系，对不同数量履带板进行仿真分析.选用履刺高度为30 mm、履带板节距为155 mm、拔模角度为60°的履带板数据，实验输出结果为履带板受到的全部反力.以两块履带板为例，结构如图11所示.

图11 双排履带板Fig.11 Double row track shoes

将多排履带板的数据进行整理，结果如图12所示.

图12 履带板数量剪切力折线图Fig.12 Number of track shoes-shear force line chart

由图12可知，随着履带板数量的增加，剪切力也随之增加，增幅约为60%.因此整车的受力为

（13）

4 履带板优化设计

4.1 履刺花纹挤压区域面积对剪切力的影响

履带板的设计参数对剪切力有很大影响，为了探究挤压区域面积对履带板剪切力的影响，选用履刺高度为30 mm、履刺节距为155 mm、履刺角度为60°的无啮合区域履带板模型.图13为单因素试验结果.由于增加履带板水平距离会令挤压区域变形，失去挤压能力，因此选择增加成对履带花纹之间的纵向距离来增大挤压区域.履刺角度β为

β=π-2α

（14）

图13 单因素实验结果Fig.13 Single factor test results

由图13可知，当履带花纹的垂直距离增加时，履带板受到的剪切力的变化幅度很小，挤压区域的面积大小对履带板受到的剪切力影响较小.

4.2 履带板节距、履刺高度、履刺角度对履带板剪切力的影响

为了研究不同履带板设计参数对于履带在积雪路面时受到的剪切力的影响，考察的因素有履带板的履刺高度、履带板节距、履刺角度，将剪切力作为考察指标.本实验采用Box-Behnken进行分析，将考察因素作为自变量，将剪切力作为响应值.确定的因素水平如表2～3所示.

表2 Box-Behnken实验因素水平Tab.2 Box-Behnken test factor levels

对实验进行方差分析，结果如表4所示.由方差分析可知，履刺高度、节距对剪切力的影响最为显著，履带板的拔模角度对履带板的剪切力影响较小.图14为Box-Behnken实验响应图.将显著性因素作为自变量，剪切力作为响应值，随着履刺节距和履刺高度的增加，履带的剪切力也不断增加，履带板受到的剪切力与履刺高度、履带板节距呈正相关.

表3 Box-Behnken实验结果Tab.3 Box-Behnken test results

图14 Box-Behnken实验响应Fig.14 Box-Behnken test response

4.3 履带板设计参数优化

考虑到履带板尺寸参数的增大会造成车辆行驶噪音大、动载冲击强、转向困难等问题，对履带板的设计参数进行优化，优化目标为在符合通过性的前提下尽量减小履带板的尺寸参数.优化对象参数选择履带板节距和履刺高度，约束条件选择LF1352森林运载车在林区的爬坡能力.履带板节距的最小值为145 mm，履刺高度的最小值为0 mm.

优化的目标函数为履带板体积的最小值，即

（15）

约束条件为：

1） LF1352森林运载车能通过45°的坡面角度，约束函数为

W（ξ+sin 45°）≤2F总

（16）

式中，ξ为履带板的阻力系数，本文取0.2.

2） 本文中垂直于履刺剪切面的力的计算公式基于Rankine被动土压力理论，即

（17）

优化结果如表5所示.

表5 优化结果Tab.5 Optimization results

5 结 论

2） 通过研究履带板因素对剪切力的影响，得出履刺高度和履带板节距对剪切力的影响较为显著，并呈正相关.剪切力随着履带板节距、履刺高度的增加而增大.拔模角度、挤压区域面积对剪切力影响较小.

3） 得出了履带车辆在积雪路面上整条履带板与单一履带板的受力关系，推导出了整条履带板的剪切力计算公式，并对履带板的结构尺寸进行了优化.