作者于2018年1月3日在重庆市涪陵区渝东南农业科学研究院“涪杂2号”试验田中，采用5点取样法统一选取9月上旬播种育苗、生态健康、外观无损伤、无病害的20颗青菜头制作压缩试件，其中10颗进行轴向压缩试验，另外10颗进行径向压缩试验.

参照木材抗压强度试验方法，结合青菜头实际尺寸，将试件制成长、宽、高都为10 mm的正方体，试件尺寸公差为±0.5 mm.试件如图 1所示.

本试验采用WDW3100型微机控制电子式万能试验机，由于农业物料具有粘弹性和小力大变形的特点，因此必须采用连续加载将力和变形真实地记录下来，该万能试验机配备全数字测量控制系统，作业时能够实现计算机自动控制和数据自动采集^[15].

试验机最大测试力有1，10，100 kN自选，有效测力范围：0.4%~100%(F.S)(内外不分档)，力分辨率：1/240 000(全程分辨力不变)，力示值相对误差为±0.5%，试验速度范围为0.005~500 mm/min(无级调速，任意设置)，速度示值相对误差为±0.5%，可依据应力、应变、位移3种闭环控制方式求出最大力、抗拉压强度、剪切强度及弹性模量等参数.

如图 2所示，将试件放置在压缩工作台中心，操作控制手柄将移动横梁移动到距试件顶端约5 mm处，选择1 kN测力传感器、压缩速度1 mm/min、断裂百分比10%、采样频率20点/s后开始压缩.由于本次试验仅需要得到试件的最大抗压强度及弹性模量，故当力-位移曲线出现急速上升趋势时即试件被压实后，手动停止试验并输出试验结果，缩短单次试验时间.

由图 3所示，压缩试件发生劈裂，与木材抗压试验中最常见的试件破坏形状相似^[16]，结合力-位移输出曲线分析可知，加工的青菜头压缩试件在压缩试验中表现良好，达到压缩破坏标准，实现预期试验目标.

如图 4所示为压缩试验曲线，曲线表明压缩试验中试件经过了4个阶段：

第一阶段为万能试验机压头靠近试件阶段(0~1.3 mm)，为了便于试件摆放在球面滑动支座及出于保护机器考虑，将压头调节至试件上方一定高度后进行试验.故从图中可以看到，0~1.3 mm阶段的力为试验机压头与机器内部作用力，呈非线性趋势.

第二阶段为试件弹性变形阶段(1.3~2.4 mm)，随着压缩力的增加，从图中弹性段起点至弹性段终点区间内，青菜头试件的压缩力与位移呈现出近似线性变化的关系，表明青菜头试件弹性变形阶段符合胡克定理.

第三阶段为试件粘弹性变形阶段(2.4~3.5 mm)，青菜头与金属的性质不同，不仅表现为各向异性，而且具有粘弹性.一定温度下，大分子变形要达到与外力相适应的平衡状态，中间有一个松弛过程，由于松弛过程的存在，使青菜头试件的变形往往要落后于应力的变化，产生滞后效应，此阶段压缩力的增长速度比位移增长速度快.

第四阶段为试件屈服阶段(3.5~5.3 mm)，此时压缩力保持在一定范围内波动，位移迅速增加.当位移到达完全压实点后，由于试件被完全破坏，压缩力在短时间内迅速增大.

青菜头缩短茎压缩试验中的表现与金属材料压缩试验表现不同，说明了青菜头缩短茎这一材料具有一定的特殊性.

通过青菜头缩短茎的压缩试验，可以得到在缩短茎切割部位轴向、径向方向的一系列力学参数，如抗压强度、屈服强度、弹性模量等，为切割青菜头的数值模拟分析提供支撑数据.作者对青菜头缩短茎切割部位轴向、径向各做10组试验，得到了两个方向上的抗压强度、屈服强度、弹性模量，如表 1所示.

弹性模量是材料重要的力学性能参数.在压缩试验力-位移曲线中，可以根据公式

求解材料的弹性模量，其中E为材料弹性模量，F为试件受到的压力，A为试件的截面积，l为试件被压缩长度，万能试验机根据以上3个参数，自动输出材料的弹性模量.

分析可得(表 2)，径向抗压强度最大值为2.31 MPa，最小值为1.71 MPa，平均值为1.95 MPa，标准差为0.22，变异系数为11%.径向弹性模量最大值为5.89 MPa，最小值为3.87 MPa，平均值为4.91 MPa，标准差为0.66，变异系数为13%.

轴向抗压强度最大值为2.34 MPa，最小值为1.09 MPa，平均值为1.65 MPa，标准差为0.38，变异系数为23%，轴向弹性模量最大值为8.95 MPa，最小值为3.18 MPa，标准差为1.59，变异系数为34%.

对比轴向、径向得到的力学参数可知，轴向与径向通过压缩试验得到的弹性模量有差异，这是由于青菜头缩短茎自身的材料结构特点决定的.观察发现青菜头缩短茎外圈存在纤维，纤维本身的数量、强度以及连接形式等要素决定了青菜头缩短茎轴向、径向方向的力学性质差异.

本文采用斜面法测定青菜头缩短茎相对65 Mn弹簧钢的静摩擦系数^[17].由图 5所示，当θ由0缓慢增加到一定值时，试件在斜面上开始匀速滑动，此时的试件受力如公式(2)所示，合外力为0，化简得到静摩擦系数μ_s关于倾角θ的表达式，即倾角的正切值，由公式(3)所示.

根据斜面法测定静摩擦系数的原理，作者设计了一个摩擦系数测定斜面仪，如图 6所示.它主要由合页、限位螺钉、标尺、PP塑料板、丝杆、加厚螺母组成.试验时，如图 7所示，将65 Mn弹簧钢板固定在斜面PP塑料板上，青菜头缩短茎放在钢板上，通过旋转丝杆让斜面慢慢抬高，直至试件发生滑动时，使用手机自带角度测量软件，测量斜面倾角θ值即可.

本文测定了“涪杂2号”青菜头缩短茎对于65 Mn弹簧钢的静摩擦系数，共准备了5个生长情况相同的缩短茎试件，为了减小随机误差，每个试件测试3次，然后求其平均值，静摩擦系数测量平均值为0.58.

本文中动摩擦系数测试平台分为两部分(图 8)，一部分为摩擦系数测定斜面仪，另一部分为高速摄像分析系统，摩擦系数测定斜面仪已在前文中详细介绍.

由图 8所示，高速摄像分析系统硬件由高速摄像机、LED补光灯、移动工作站构成，软件由高速摄像机控制软件界面(phantom camera control，简称PCC)构成，该软件具有设置拍摄分辨率、拍摄速率、触发位置等参数的功能，同时也能对拍摄图像处理和存储.作者将使用该软件拍摄试件在摩擦系数测定斜面仪上滑动的过程，分析计算试件在65 Mn弹簧钢板上的动摩擦系数.

动摩擦系数测定采用与静摩擦系数测定方法相同的斜面法，通过高速摄像计算出试件在摩擦系数测定斜面仪上滑动的距离和时间，进而计算出试件滑动的加速度，根据斜面上试件受力平衡方程可计算出试件在65 Mn弹簧钢板上的动摩擦系数，具体原理如下：

如图 9所示，在高速摄像机控制软件界面设定拍摄频率为每秒1 000张图片，故每一帧画面试件滑动时间为1 ms，通过计算出起始位置与终止位置之间帧数差值，可以得到试件滑动的时间，同时通过坐标纸读出试件滑动的距离.由公式(4)及公式(5)可以计算出试件在斜面上的加速度.

S为试件滑动的距离，单位m；t为试件滑动的时间，单位s；a为试件在斜面上的加速度，单位m/s².

通过对斜面滑动的试件进行受力分析得到方程(6)，化简得到关于试件加速度、摩擦系数测定斜面仪倾角的动摩擦系数公式(7)，将相关参数带入后可得试件在65 Mn弹簧钢板上的动摩擦系数μ_k.

公式中m为试件质量，单位g；θ为摩擦系数测定斜面仪倾角；μ_k为试件动摩擦系数；g为重力加速度，单位m/s².

由静摩擦系数的测定可知，试件在斜面上滑动的临界角度大约为30°，故将制作好的试件放在摩擦系数测定仪斜面，以40°、45°、50°倾角测定其动摩擦系数，本次试验共测定5个试件.

为了避免65 Mn弹簧钢板表面污渍对测试结果的影响，试验时应将试件放在相对光泽、干净的位置，对于试件滑动影像后处理阶段，尽量选择短行程来分析计算，避免钢板表面粗糙度不同对试验结果造成影响.

将表 3中相关参数代入公式(7)中，计算出青菜头缩短茎在65 Mn弹簧钢板上的动摩擦系数.

自制试验平台测定的动摩擦系数最大值为0.509 6，最小值为0.426 3，平均值为0.468 1，标准差为0.027 7，变异系数为5.9%.试验结果看出，动摩擦系数变异系数较低，试验结果离散程度较低，根据材料摩擦性质可知，材料的动摩擦系数与材料质量、接触面积、斜面角度等因素无关，验证了基于高速摄像测定青菜头动摩擦系数方法的可行性.

选取重庆地区青菜头主要种植品种“涪杂2号”为研究对象，通过试验测定青菜头缩短茎的抗压强度、弹性模量、动摩擦系数、静摩擦系数等基础物理力学参数，为青菜头收获机械仿真优化提供数据支撑，为青菜头收获和加工机器系统的研发提供基础理论依据.

1) 青菜头缩短茎径向抗压强度在2.31~1.71 MPa范围内，径向弹性模量在5.89~3.87 MPa范围内；轴向抗压强度在2.34~1.09 MPa范围内，轴向弹性模量在8.95~3.18 MPa范围内.

2) 基于斜面法采用高速摄像测定青菜头缩短茎对于65 Mn钢的动摩擦系数平均值为0.468 1，同时测定其静摩擦系数平均值为0.577 4.