研究发现，马铃薯70%的块茎损伤出现在收获环节^[1]，机械损伤严重影响了马铃薯的品质、产量和经济效益^[2-5]。随着马铃薯收获机械化、加工自动化的快速发展，马铃薯由碰撞导致的损伤问题亟待解决。有研究根据马铃薯果肉的变色来度量其损伤程度，提出了损伤的评价标准^[6-8]。吴永根等^[9]将经冲击的果实在自然环境中放置24 h，损伤的果肉会呈现褐色，从而获得损伤体积、损伤深度等试验指标。桑永英等^[10]根据淀粉遇碘变色原理，对碰撞后的果实涂抹碘，根据颜色深浅划分损伤程度。康璟等^[11]提出通过多因素试验研究马铃薯机械损伤的建议，如从挖掘深度、下落高度等角度分析马铃薯的机械损伤，从应力应变及其临界极限值方面研究马铃薯的机械损伤等。王咏梅等^[12]研究了马铃薯收获中的机械损伤，分析了机械损伤产生的过程和机理，建议对马铃薯动态力学特性进行深入研究。李辉等^[13]通过对比国内外的马铃薯机械损伤研究，鼓励多从实际作业中发现问题、解决问题。魏忠彩等^[14]系统研究了收获、清选和分级对马铃薯机械损伤的影响。冯斌等^[15]通过试验及参数的数学关系获得马铃薯块茎碰撞恢复系数，并分析跌落高度、碰撞材料等跌落过程的参数与块茎碰撞恢复系数的关系。吕金庆等^[16]研究了跌落高度、升运链倾角和线速度对马铃薯损伤综合指数和伤薯率的影响。Ito等^[17]研究了跌落高度和碰撞材料对马铃薯表面损伤的影响。Dwelle等^[18]研究了温度、放置时间、气体压力对马铃薯内部黑斑瘀伤形成的影响。Thomson等^[19]研究探讨了马铃薯块茎的特性及影响马铃薯表面损伤的处理技术。Nikara等^[20]采用扫描电镜(Scanning electron microscope, SEM)分析了马铃薯碰撞后组织的微观力学变化。卢立新等^[21]通过冲击力传感器研究苹果跌落的冲击力学特性，获得加速度与时间、冲击力与变形量、应力与应变以及跌落冲击弹性恢复系数等动态关系。李晓娟等^[22]通过悬摆式碰撞试验，采用加速度传感器，得到苹果碰撞过程的加速度−时间曲线及数学表达式。Brusewitz等^[23]使用力学传感器研究苹果的跌落，同时测得最大撞击力、撞击持续时间、到达最大撞击力所需要的时间、冲量等冲击特性的参数。刘治震^[24]通过冲击力传感器对马铃薯进行跌落碰撞试验，得到马铃薯的临界损伤高度。卢琦^[25]通过平膜盒测力传感器获得马铃薯的瞬时冲击力和冲击应力；通过与马铃薯的坚实度比较，得到马铃薯的损伤情况。胡奔^[26]通过力学传感器获得马铃薯跌落碰撞的冲击力，经过转换研究得到马铃薯损伤的极限冲击应力。冯斌等^[27]通过压力传感器测量马铃薯的冲击特性参数与损伤综合指数的关系。但是传统的测量方法不精确，普通的力学传感器只能测定马铃薯与其他材料发生碰撞时冲击力的大小，无法测定碰撞时的应力分布情况。对马铃薯与各作业部件的接触力学分布特性进行研究，可为马铃薯的机械化收割和自动化减损提供参考，具有重要的意义。

有学者用超声波技术^[28]或TekScan5051柔性薄膜网格压力传感器^[29]测量了苹果准静态下的压缩接触应力分布。美国和日本通过微胶囊颗粒缓释控制技术研制出了压力感应胶片，实现了动载下接触应力的精确测量，并成功应用于医学^[30-31]。采用感压胶片对苹果挤压和碰撞的接触应力进行分析，并获得损伤面积和体积^[32-33]。采用感压胶片对香梨跌落碰撞的接触应力分布进行研究，确定接触应力分布和香梨果实损伤间的关系^[34-35]。

本文采用Prescale感压胶片及配套的压力图像数字解析系统对马铃薯跌落碰撞的接触应力分布进行研究，精确测量马铃薯在动载下的接触应力分布特性，研究马铃薯在动载下的力学特性和损伤规律，以实现对其机械损伤的可靠性评估和预测。

1 材料与方法 1.1 材料

马铃薯试样为华中农业大学农业农村部马铃薯生物学与生物技术国家重点试验室试验田培育的‘中薯5号’，早熟品种，从出苗到收获约60~80 d，株高约50 cm，休眠期短。‘中薯5号’具有较好的增产性，薯块大而整齐、表面光滑、芽眼极浅、综合性状好，结薯集中，商品薯率高^[36-37]。

收获期为2019年6月，在确保马铃薯无虫害、无损伤后迅速冷藏。马铃薯的含水率(w)为77.36%，贮藏温度为−5~5 ℃，相对湿度为80%~90%。

试验材料为土块、橡胶、塑料和钢板。土块采集自马铃薯收获时的试验田，含水率(w)为22.58%；橡胶为工业橡胶中的丁晴橡胶；塑料的材料是ABS。钢板采用65Mn钢，为我国薯类收获机上升运链建议采用的材料^[38]。

高速摄像机（Pco.dimaxHD）为德国PCO公司生产。高速摄像机放置在距离试验台架水平2 m处，记录马铃薯跌落碰撞时的运动情况，选定的拍摄帧率为2 000帧/s，曝光时间为942 ns，像素分辨率为1 920 ×1 080 。

Fuji公司研发的感压胶片(Prescale^®双层特超低压LLLW型)测量范围为0.2~0.6 MPa，测试精度≤±10%，配套的压力图像数字解析系统V2.0(FPD-8010E型)。

EPSON平板扫描仪(V370)为爱普生公司生产。

马铃薯跌落碰撞测试系统示意图见图1。

1.2 方法 1.2.1 马铃薯跌落试验

马铃薯试样通过调节阀固定在跌落孔中，松开调节阀，马铃薯试样自由下落，保证其长轴扁平部位与不同的材料碰撞。反弹后在再次落地前接住马铃薯试样，保证马铃薯试样经过反弹达最大高度，也避免二次碰撞。高速摄影机记录全过程。

1.2.2 接触应力分布的测定

将感压胶片放置在碰撞接触材料表面，使感压胶片能够获取完整的碰撞接触部位。碰撞后的感压胶片，通过扫描仪进行平滑处理，获得接触应力分布的特征图，然后经相配套的压力图像数字解析系统V2.0读取特征图并进行数值解析，最后获得接触面积、碰撞冲击力、接触应力及其分布等参数。

1.2.3 马铃薯冲击压缩变形量的测定

通过马铃薯的冲击压缩变形量来衡量马铃薯的损伤。马铃薯的冲击压缩变形量越大，损伤越大。

冲击压缩变形量(D)表征马铃薯跌落碰撞时弹性形变的最大值，定义公式如下：

$ D=\frac{\Delta s}{d} ,$

(1)

式中， $\Delta {{s}}$ 是指形变量，马铃薯从碰撞开始到压缩到最大值（下一秒会发生反弹）(图2)；d是与接触材料发生碰撞时，马铃薯竖直方向的直径。

形变量 $\Delta s={\int }_{0}^{t}{v}_{t}{\rm d}t$ ， ${v}_{t}={v}_{0}+{\int }_{0}^{T}a{\rm d}T$ ，T为碰撞弹性变形时间，通过高速摄影技术获得，t=2T。

1.3 马铃薯试验方案设计

通过Design-expert设计了3因素3水平的组合正交试验方案，研究因素及组合正交试验设计方案的编码值和实际值的对应关系如表1所示。

跌落高度是根据马铃薯收获机振动筛与地面的高度进行筛选的。碰撞材料是根据马铃薯收获时可能发生碰撞的材料选取的。马铃薯质量是将收获后的马铃薯进行分级所得。最后通过分析较显著的因素，补充详细完整的因素水平，进行单因素试验。

2 结果与分析 2.1 方差分析及响应曲面分析

组合正交试验的方差分析见表2。跌落高度和碰撞材料的P<0.000 1，对冲击压缩变形量的影响极显著；马铃薯质量的P= 0.026 5<0.05，对冲击压缩变形量的影响是显著的。而两两因素间的交互作用影响不显著，因素的二次方影响也不显著。通过方差分析获得冲击压缩变形量(D)与3个因素的二次回归模型。

$ \begin{split} D= & 11.27-2.60A+3.31B-1.17C+0.45AB-0.25AC-\\ & 0.22BC+0.16{A}^{2}-1.22{B}^{2}+0.082{C}^{2}{\text{。}} \end{split} \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!$

(2)

失拟P=0.747 5>0.05，影响不显著，说明方差分析模型是可信的。

对跌落高度、碰撞材料和马铃薯质量的响应曲面分析如图3所示，二维的投影曲线反映因素间的关系，曲线在不同坐标轴具有不同疏密程度，疏密程度反映因素对冲击压缩变形量的不同贡献程度，曲线越密集说明对冲击压缩变形量的贡献越大。

由图3a、3c可看出，随着跌落高度增加，曲线斜率均较陡峭，说明跌落高度对冲击压缩变形量的影响较显著，冲击压缩变形量随着跌落高度的增加而增大。由二维的投影可知，跌落高度的曲线密集程度比马铃薯质量的高，比碰撞材料的低。由图3a、3b可看出，不同的碰撞材料中马铃薯的冲击压缩变形量差异较大，碰撞材料的曲线密集程度比跌落高度、马铃薯质量的高。由图3b、3c可看出，当马铃薯质量较小时，曲线较陡峭，说明此时马铃薯质量对冲击压缩变形量影响较显著。马铃薯质量的曲线密集程度比跌落高度、碰撞材料都低。综上可得各因素对冲击压缩变形量贡献的程度为碰撞材料>跌落高度>马铃薯质量。因此后续着重分析碰撞材料和跌落高度对马铃薯跌落碰撞损伤的影响。

2.2 马铃薯跌落碰撞的损伤规律

为减少马铃薯外形差异对试验结果的影响，保证试验因素以外的环境和条件一致，选取的马铃薯为100~150 g，以减小其质量对试验的误差影响。不同跌落高度下马铃薯与不同材料碰撞的冲击压缩变形量见图4。

由图4可知，随着跌落高度的增加，不同碰撞材料的冲击压缩变形量均呈上升趋势：土块和马铃薯的曲线交织上升，65Mn钢、塑料ABS和丁晴橡胶曲线上升趋势比较清晰。65Mn钢的冲击压缩变形量随着跌落高度增加而增加，在跌落高度>300 mm后，曲线呈近似直线上升趋势，冲击压缩变形量比其他材料的大。由于硬度比较大，65Mn钢碰撞时的弹性变形较小，使马铃薯的形变较大，因而65Mn钢的冲击压缩变形量比其他材料大。

塑料ABS的冲击压缩变形量曲线，首先是逐渐上升，在跌落高度>400 mm后趋于平缓。塑料ABS具有一定的弹性，对马铃薯的碰撞变形起到缓冲作用，因而塑料ABS的冲击压缩变形量处于中间水平，与土块和马铃薯的冲击压缩变形量相差不大。

马铃薯和土块的冲击压缩变形量曲线在跌落高度>400 mm后上升，在跌落高度>600 mm之后趋于平缓。马铃薯具有弹性缓冲性能，因此两两碰撞时，能够减少变形，而使冲击压缩变形量不大。土块本身具有一定的硬度，但由于收获时的土块含水量较高，能够对马铃薯的碰撞有一定的缓冲作用，加之土块具有较多的孔隙，跌落碰撞时的接触面积较大，从而减少了形变，因而土块的冲击压缩变形量处于中间水平。

丁晴橡胶的冲击压缩变形量在跌落高度为300 mm时最小，在跌落高度>400 mm后，曲线呈近似直线上升，在跌落高度为600~800 mm，冲击压缩变形量的差异很小。丁晴橡胶的冲击压缩变形量整体上比其他材料低。丁晴橡胶是缓冲材料，工业设计材料内部具有较多孔隙，能够对马铃薯起到缓冲减震的作用，使马铃薯的冲击变形较小，因而丁晴橡胶的冲击压缩变形量最小。

综上，不同碰撞材料造成冲击压缩变形量的差异与其接触应力、冲击力、接触面积等有关。

2.3 马铃薯碰撞接触应力的分布特性 2.3.1 不同碰撞材料的接触应力分布特性

马铃薯与5种接触材料进行碰撞试验，通过感压胶片及其压力数字解析系统获得马铃薯碰撞的接触应力及其分布。马铃薯与各种材料在不同高度碰撞时的接触应力与接触面积的关系见图5，不同接触材料在不同高度碰撞后的应力分布如图6所示。

马铃薯与各种材料在不同高度碰撞时，接触应力≤0.50 MPa占主要的接触面积。接触应力为0.50~0.60 MPa应是应力的峰值，因为接触应力>0.60 MPa的接触面积非常小，约为0.01~1.00 mm²，几乎可以忽略不计。接触应力在≤0.50 MPa的接触面积最大，对马铃薯的损伤起主要贡献作用；而接触应力为0.50~0.60 MPa的接触面积很小，不是造成马铃薯损伤的主要应力。

当接触应力≤0.20 MPa时，不同跌落高度间的接触面积相差不大。当接触应力>0.20 MPa后，跌落高度为200和300 mm的接触面积相差很小，而跌落高度≥400 mm的接触面积比200和300 mm的接触面积大很多。说明随着跌落高度增加，接触应力增大，接触面积也增大，即高应力占主要的接触面积，是导致马铃薯碰撞损伤的主要原因。

由图6可以看出，马铃薯与65Mn钢碰撞，当跌落高度为200 mm时，低应力区域(绿色区域)占主要面积；当跌落高度≥300 mm后，低应力区域分布在边缘，高应力区域占主要面积。马铃薯与塑料ABS、马铃薯、土块和丁晴橡胶碰撞，当跌落高度≤300 mm时，低应力区域占主要面积；当跌落高度≥400 mm后，低应力区域分布在边缘，而高应力区域占主要面积。从图6可看出，马铃薯与土块碰撞时的接触应力分布的轮廓非常分散，存在较多孔隙，同时接触面积也比其他材料大；马铃薯与65Mn钢、塑料ABS、马铃薯、丁晴橡胶碰撞的轮廓均呈较规则的椭圆形。

以上分析说明，与65Mn钢碰撞，跌落高度≥300 mm后，马铃薯开始出现损伤；与塑料ABS、土块、马铃薯、丁晴橡胶碰撞，跌落高度≥400 mm后，马铃薯开始产生损伤。

与65Mn钢碰撞，跌落高度<300 mm时，≤0.20 MPa的区域(低应力区域)占据主要的接触面积 (74.57%)；在跌落高度≥300 mm后，>0.20 MPa的区域(高应力区域)为主要接触面积 (74.37%)。与塑料ABS、马铃薯、土块、丁晴橡胶碰撞时，马铃薯在跌落高度<400 mm时，≤0.20 MPa的区域(低应力区域)占主要的接触面积，分别为200 mm的74.37%、60.38%、61.73%和65.68%，300 mm的58.23%、76.98%、71.42%和55.41%；在跌落高度≥400 mm时，>0.20 MPa的区域(高应力区域)占主要接触面积，分别为83.52%，73.31%，66.22%和83.87%。因此可以将马铃薯受到0.20 MPa的应力作为其损伤的临界应力。

2.3.2 马铃薯接触面积和接触应力

从图7可以看出，随着跌落高度的增加，不同材料的接触面积随之上升，且接触面积(A_c)随着跌落高度(H)上升呈高度线性正相关，决定系数(R²)均大于0.95。

不同材料的接触面积由高到低依次为：土块、65Mn钢、丁晴橡胶、马铃薯、塑料ABS。从图7可以看出，土块的接触面积十分分散，是因为土块本身具有较多孔隙，碰撞时通过增大接触面积，降低了马铃薯的变形和损伤。马铃薯与65Mn钢碰撞时的接触面积较大，是由于碰撞时65Mn钢塑性变形很小，形变集中于马铃薯果肉，使马铃薯碰撞变形和接触面积较大。马铃薯与马铃薯、丁晴橡胶碰撞时的接触面积较大，且马铃薯和丁晴橡胶均具有缓冲减震的作用，因此接触面积较大能够减少对马铃薯的损伤。马铃薯与塑料ABS碰撞时接触面积较小，应力集中，是马铃薯碰撞损伤严重的主要原因。

由此可见，马铃薯接触面积在不同碰撞材料上变化的原理不同，导致马铃薯损伤程度的不同。

如图8所示，马铃薯与各材料碰撞的接触应力随着高度增加呈正态分布。65Mn钢板和塑料ABS的接触应力较大，范围为0.20~0.31 MPa，随着跌落高度增大略微起伏；而丁晴橡胶、马铃薯土块的接触应力明显小于钢板和塑料ABS，为0.18~0.26 MPa。

不同材料的接触应力具有显著差异，说明接触应力与马铃薯跌落的接触面积存在关联。

2.3.3 马铃薯碰撞接触应力分布与损伤的关系

马铃薯的接触应力（P）与接触面积（A_c）存在关联，而且接触应力与接触面积的乘积，即为碰撞冲击力（F），即：

$ {{F}} = {{P}} \times {{{A}}_{\rm{c}}}{\text{。}} $

(3)

马铃薯冲击压缩变形量与碰撞冲击力的线性拟合见图9，两者高度线性相关。冲击压缩变形量的回归方程见表3，除了土块(R²=0.963 4)，其他材料的R²均大于0.99。说明通过感压胶片测量的接触应力和接触面积能够准确地预测和评估马铃薯的损伤。

3 结论

本研究采用感压胶片及配套的压力图像数字解析系统分析马铃薯跌落碰撞的接触应力及其分布特性，获得马铃薯跌落碰撞的损伤预测模型。研究了碰撞材料和跌落高度对马铃薯碰撞损伤影响的数学关系，建立了相关模型。在方差分析中，马铃薯质量对其碰撞损伤同样具有显著影响，马铃薯质量对马铃薯碰撞损伤机理有待进一步的研究。

马铃薯组合正交试验的响应曲面分析表明，跌落高度、碰撞材料、马铃薯质量对马铃薯碰撞的冲击压缩变形量影响显著，且影响程度由大到小为碰撞材料>跌落高度>马铃薯质量。

跌落高度为200~800 mm时，接触应力≤0.50 MPa占主要的接触面积，对马铃薯的损伤起主要贡献作用。马铃薯接触应力0.50~0.60 MPa是接触应力峰值，因为马铃薯接触应力>0.60 MPa的接触面积约为0.01~1.00 mm²，可忽略不计。不同碰撞材料马铃薯损伤的跌落高度不一致，碰撞材料为65Mn钢时，在跌落高度≥300 mm后，马铃薯开始出现损伤，而塑料ABS、土块、马铃薯、丁晴橡胶在跌落高度≥400 mm后，马铃薯开始发生损伤。当碰撞材料的跌落高度<400 mm(65Mn钢<300 mm)时，≤0.20 MPa的区域占据主要的接触面积，当跌落高度≥400 mm(65Mn钢≥300 mm)时，(0.20，0.60] MPa的区域为主要接触面积。因此可以将0.20 MPa的接触应力作为马铃薯损伤的临界应力。跌落高度与接触面积呈高度线性正相关，R²均大于0.95。随着跌落高度增大，65Mn钢和塑料ABS的接触应力有略微起伏，为0.20~0.31 MPa；马铃薯、土块和丁晴橡胶的接触应力明显低于65Mn钢和塑料ABS，为0.18~0.26 MPa。冲击压缩变形量与接触应力和接触面积的乘积即碰撞冲击力呈高度线性相关，R²大于0.96，说明冲击压缩变形量与碰撞冲击力的回归方程模型能够准确地预测和评估马铃薯的损伤。

　致谢：感谢宋波涛教授和陈惠兰教授对试验给予的帮助和鼓励！