付 函,刘恩秀,杨 洲,2,杨靖恺,张汉尧,谢博伟,杜文迪,段洁利
(1 华南农业大学 工程学院,广东 广州 510642; 2 嘉应学院/广东省山区特色农业资源保护与精准利用重点实验室,广东 梅州 514015)
国家统计局数据[1]显示,香蕉是我国排第五位的大宗水果,2019 年种植面积达 3.303×105hm2,产量达 1.1656×107t。在收获时为了长途运输,必须采摘未经催熟而暂不能食用的青果,即青香蕉。青香蕉受外力碰压或摩擦易发生机械损伤。机械损伤是影响果品质量的重要因素。机械损伤既能直接降低果实的外观品质,影响消费者的购买欲望[2],更会为某些微生物的侵染提供通道[3],致使果实腐败交叉感染,给果农和商家造成巨大的经济损失。研究果实产生机械损伤敏感性的成因及其影响因素,探索有效的防护措施,对于确保果品质量具有重要的意义。
碰撞损伤是水果机械损伤的主要类型,其敏感性被定义为单位吸收(或碰撞)能量下的损伤体积[4]。由于碰撞损伤敏感性的定义考虑了果实损伤中的吸收能量和损伤体积2个重要指标,因此被广泛地用于评价水果发生冲击损伤的可能。目前,有学者已经研究了苹果、梨、桃、番茄、猕猴桃等多种类球形水果的碰撞损伤敏感性,并揭示了相关的影响因素。这些因素既包括果实固有的物理特性(如尺寸、形状、硬度、含水率),也包含外部因素(如下落高度、碰撞能量、冲击力、冲击速度)。例如,Van Zeebroeck等[4]研究番茄时发现果实的碰撞损伤敏感性随其成熟度的增加而增加;番茄的温度与曲率的交互效应对其损伤敏感性存在显著影响,曲率越大、温度越低则番茄越易损伤。苹果和桃子的碰撞损伤敏感性与其成熟度之间也有类似规律[5-6]。
除了损伤敏感性,损伤率[7]、损伤抵抗性、损伤指数[8]、跌落损伤指数[9]等指标也被用于评价水果吸收(抵抗)外力致损的可能。其中损伤抵抗性的定义与损伤敏感性的定义相反[10]。Schoorl等[11]发现,苹果的损伤体积与吸收能量呈强线性正相关,于是采用损伤体积与吸收能量之间简单线性回归模型的斜率表达果实的碰撞损伤抵抗性。Hadi等[12]为了对油棕榈果的损伤进行分类评价,提出了损伤指数的概念,根据果实的成熟度,确定了损伤指数与游离脂肪酸(Free fatty acid,FFA)含量之间的关系。
目前有关香蕉机械损伤的研究文献很少,特别是针对青香蕉碰撞损伤敏感性的研究鲜见报道。Banks等[13]测量了准静态压缩下青香蕉损伤的阈值力及其随加载时间变化的规律,发现香蕉损伤敏感性随其成熟度的增加而降低,且损伤更易被辨识,同时损伤深度逐渐变浅。Akkaravessapong等[14]发现香蕉贮藏的相对湿度对其碰撞损伤敏感性没有显著影响。Kajuna等[15]利用摆式冲击的方法对比研究了贮藏时间对青香蕉与芭蕉碰撞损伤敏感性的影响,发现贮藏初期芭蕉的损伤敏感性较低,但贮藏后期其敏感性快速升高并超过香蕉。Bugaud等[16]从物理和化学的角度考虑,通过冲击试验揭示了基因表型和采后贮藏条件对香蕉果皮损伤阈值的影响。Maia等[17]验证了机械应力诱导下‘普拉塔’香蕉在室温条件下的物理和生理代谢变化。Fernando等[18]测量并分析了香蕉在长途运输过程道路振动对其机械损伤程度的影响,用于优化减损包装设计。
香蕉的结构特征与其他水果均不相同,本文结合青香蕉果穗的结构特点,通过嵌套和单因素试验设计,揭示梳层节点及其嵌套的排层、冲击水平和果实表面区域对果实碰撞损伤敏感性的影响;通过混合正交试验设计,揭示各因素对碰撞损伤敏感性的影响显著性以及影响主次关系。自行设计和搭建摆动冲击试验台,开展果实的碰撞损伤试验。
香蕉是一种草本植物,一棵香蕉树只结一束果穗,每束果穗一般包含7层果梳,每层果梳又包含内外2排果指,每根果指的长度依据香蕉品种的不同在10~30 cm之间[19],香蕉果穗和果实的形态和结构组成如图1所示[19-21]。本文以广州本地广泛种植的巴西蕉AAA品种为例开展试验,在2020年4月1日(断蕾后第110天)采收果穗并运回实验室,在室温条件下通风保存至整体温度稳定。试验时,将果指保留果柄在靠近梳柄处单根取下,挑选位于蕉梳中部、颜色、大小和形状一致、无机械损伤、无病虫害的果指。将选出的果指按生长节点和内排、外排分组,用清水清洗后软布擦干。受光照、植株长势等环境因素的影响,同一束果穗不同梳层、同一梳层的内外排果指以及同一果指表面的不同区域,果指的生长状况均可能存在差异,从而影响其损伤敏感性的差异。为此,本试验重点研究梳层节点、排层和果指表面区域3个重要因素对果实碰撞损伤敏感性的影响。
图1 香蕉果穗和果实的形态和结构组成Fig.1 Morphology and structure of banana bunch and fruit
果蔬碰撞损伤试验的实施一般选用跌落或摆动冲击2种类型的试验装置。由于香蕉果指是条状,采用跌落试验的方式难以精准控制其表面的碰撞位置,并且果指易发生弯转,不便计算碰撞部位的吸收能量。为此,参考前人开发的试验装置[10,22],设计并搭建了基于摆动冲击试验的香蕉果指碰撞损伤试验台,如图2所示。试验台中部的摆臂长0.6 m,截面为20 mm×15 mm的实心铝合金杆;摆锤碰撞头为直径40 mm的木质半球;在摆臂铰链处安装了一台角度传感器,用于测量摆角;试验台的整体框架选用铝合金型材搭建。试验时,香蕉果指被固定在具有升降功能的支撑平台上。为了避免果指碰撞区域外损伤,提前打印了多种规格的仿形部件,以适应香蕉接触表面的不同形态,增大接触面积。摆臂碰撞后的速度通过高速相机测量。试验台构成的主要硬件包括:高速相机(2 000 fps,Mini UX50,Photron),数据采集卡(26.5 kHz,NI9234,National Instrument),卡槽(cDAQ-9174,National Instrument),角度传感器(精度 0.3%,MCJSI420A-360-38AC,天宇恒创),升降台(LZ80-2,汇科)。
图2 摆动冲击试验台示意图Fig.2 Schematic of the pendulum impact test device
碰撞前,将香蕉果指试样固定于仿形部件,调节升降台以保证果指碰撞区域与木锤接触时摆臂水平。试验时,用手通过棉线将摆臂拉起,定位在设置的水平,然后自由释放,使木锤冲击在选定区域,在木锤再次冲击前用手阻挡摆臂,避免二次撞击。高速相机置于试验台一侧,平视冲击试验台,记录整个冲击过程。
参考苹果[23]、梨[24]等水果敏感性的定义,本文采用单位吸收能量下的损伤体积表达香蕉的碰撞损伤敏感性,其值越大,果实越易损。
将损伤敏感性记为SB,cm3/J;损伤体积记为VB,吸收能量为ΔE,则:
在评估损伤时,首先观察碰撞区域果皮的表面色态,用手触摸确定损伤的确切位置,然后将损伤区域的中间依次横切和纵切,用相机记录损伤截面。损伤的典型示例如图3所示。
图3 香蕉果指不同截面损伤情况示例Fig.3 Illustration of bruised banana finger in different sections
用ImageJ1.8.0软件测量损伤部位的面积和长度,分别记为S和b。损伤面积S的测量步骤如图4所示,首先建立如图4a所示的图片标尺,其次将图片二值化(图4b),然后选择适宜亮度值划分损伤部位和健康部位(图4c),计算损伤部位的面积。损伤长度b测量方法与损伤面积相似。
图4 香蕉损伤面积测量方法示意图Fig.4 Schematic of the method for measuring bruise area on banana
在苹果[25]、梨[24]等多种类球形水果损伤程度的测量中,通常将损伤部位的形状近似为半球或半椭球形,利用切片法进行积分计算。然而,香蕉横截面显示其损伤部位为3瓣状,与其他水果的损伤形状不同。为此,修改切片法测量损伤体积,将损伤在XOY截面上的边界近似为椭圆形,在YOZ截面上近似为半椭圆形,损伤形状类似于图5所示。设损伤最深处Z(0,0)=h,根据椭圆公式可得YOZ平面的损伤边界公式近似为:
图5 用于计算损伤体积的变量示意图Fig.5 Schematic of the definition of the variables for calculation of bruise volume
其中,y∈(-b/2,b/2)。
沿长度方向任意位置横截面的损伤面积设为S(y),根据实际情况S(y)在XOZ截面上最大,即S(0)=Smax。可以得到:
忽略摆锤铰链处的摩擦损耗、空气阻力和试验台振动损耗,根据能量守恒定律,香蕉果指的吸收能量ΔE为摆锤冲击能量E1与其反弹能量E2之差,即:
设摆锤的转动惯量为J,kg·m2,木锤质量为m1,摆臂质量为m2,木锤球径为r,摆臂长度为l,木锤质心到转动中心距离为l1,那么转动惯量为:
设ω为摆锤的角速度,rad/s,则钟摆的转动动能(E)为:
在进行冲击试验时,钟摆的线速度通过高速拍照的方式测量。设高速相机摄像时相邻2张图片的拍摄间隔为ΔT,分别选出摆锤与香蕉果指临近接触和分离的相邻2张图片,测量摆锤冲击和反弹的位移X1和X2,用位移与间隔时间的比计算摆锤冲击和反弹的速度,分别记为V1和V2:
由于摆锤固定于摆杆上绕摆杆一端转动,为了便于计算冲击能量,将摆锤冲击和反弹的瞬时线速度转换为角速度,分别记为ω1和ω2,设摆锤到转动中心距离为l1,那么
果指被撞部位的吸收能量为:
1.3.1 嵌套试验设计 香蕉一般有7个梳层节点,每个节点有2排果指,由于每一梳香蕉生长节点对应内外2排果指,不同节点对应的内排或外排果指无法交叉。为了确定一级因素梳层节点(A)和二级因素排层[B(A)]2个因素对果指碰撞损伤敏感性的影响,梳层节点和排层采用二级因素嵌套设计。梳层节点设为7个水平,排层为内、外2排2个水平,试验因素和水平如表1所示。每个嵌套组合重复5次,整个试验在24 h内完成,碰撞完成的样本在22 ℃的室内条件下储存24 h,使损伤特征充分显现。
图6 香蕉果指表面被撞位置标记示例Fig.6 Illustration of the impact zones on a banana finger
表1 嵌套试验因素和水平Table 1 Nested test factors and levels
表2 正交试验因素和水平Table 2 Orthogonal test factors and levels
1.3.2 单因素试验设计 香蕉果指较长,不同位置的损伤敏感性可能存在显著差异。为了检验果指表面区域对损伤的敏感性,从果指表面区域等间距选3个水平,即近果柄端、果指中部和近果顶端,如图6所示。前人研究[24-25]发现冲击水平是一种对水果碰撞损伤敏感性影响显著的因素,为了选取合适的冲击能量,在预试验中对果指中部开展了10°、15°、20°、35°和 50°多种摆角下的系列碰撞,静置 24 h 后发现:当摆角低于15°时,果指上几乎没有损伤显现;当摆角为20°时,可在果指上观察到小于1 cm2的不明显机械损伤;当摆角为35°时,发现约1 cm2的较明显损伤;当摆角为 50°时,发现大于 1 cm2但小于2 cm2的明显损伤。青香蕉品质分级标准[21]指出,若果指存在机械伤则被定为优等品以下(一等品、二等品或劣等品)。根据预试验结果,并结合青香蕉品质分级标准,本文选择了20°、35°和50°3种摆角,对应的冲击能量分别为0.52、0.88和1.17 J,研究果指表面区域和冲击水平2个因素对果指碰撞损伤敏感性的影响,每种处理重复7次。整个试验在24 h内完成,碰撞完成的样本放置在22 ℃的室内条件下储存24 h,以使损伤特征充分显现。
1.3.3 正交试验设计 为进一步探究各因素对冲击损伤敏感性影响的显著性以及影响的主次关系,在嵌套试验和单因素试验基础上,开展了混合正交试验。根据嵌套试验的结果,将B(A)因素的14个水平分为4组(见“2.1”节),因素记为生长位置,编号B ′。正交试验因素水平表见表2,每种处理重复5次。
所有数据用 Office excel 2016 软件记录;用Origin 2018 软件绘图;统计分析时,先用 IBM SPSS Statistics 22.0 对数据的正态性和方差齐性进行检验,然后通过方差分析(ANOVA)进行显著性检验,选用Duncan’s法对同一因素各水平间的均值进行多重比较。数据结果用平均值±标准差表示。
一级因素梳层节点各水平下的损伤敏感性服从正态分布,首梳至尾梳各水平下P值依次为0.33、0.79、0.45、0.63、0.46、0.09,且方差齐性检验结果为P=0.10。二级因素排层(梳层节点)各水平下的损伤敏感性服从正态分布,首梳外排至尾梳内排各水平下 P 值依次为 0.39、0.19、0.54、0.46、0.63、0.09、0.67、0.58、0.31、0.25、0.25、0.53、0.85、0.89,且方差齐性检验结果为P=0.09。
香蕉果指生长位置嵌套试验的方差分析结果显示,一级因素梳层节点间和二级因素排层(梳层节点)间的F值分别为0.17和6.94,P值分别为0.98和0.00,表明一级因素梳层节点对果指碰撞损伤敏感性影响不显著;二级因素排层(梳层节点)间的果指碰撞损伤敏感性差异显著。
为了进一步比较7个梳层节点内外2排不同水平间果指碰撞损伤敏感性的均值,利用Duncan’s法对其进行多重比较,结果如表3所示。从表3中可以看出,果指损伤敏感性的均值为1.06~1.34;同一梳层内排果指的损伤敏感性的均值总比外排果指的高。根据7个梳层节点内外2排不同水平间果指损伤敏感性均值的多重比较结果,可以将一条果穗14层不同生长位置的香蕉果指分为4组处理:第1组包括首梳内排、第5梳内排、第4梳内排、第2梳内排、第3梳内排、第7梳外排、第7梳内排、第6梳内排、第5梳外排和第4梳外排;第2组包括第7梳外排、第7梳内排、第6梳内排、第5梳外排、第4梳外排和第3梳外排;第3组包括第7梳内排、第6梳内排、第5梳外排、第4梳外排、第3梳外排和第6梳外排;第4组包括第6梳内排、第5梳外排、第4梳外排、第3梳外排、第6梳外排、第2梳外排和首梳外排。这种差异可能与果指的弯曲程度有关。香蕉果指总体向上生长,外排果指受到内排果指的影响,只能先水平再向上生长,所以内排果指的曲率半径总比外排果指的大,图1a与图1b较好地展现了这一点。Abedi等[26]的研究结果表明,当冲击水果的冲击能量超过0.24 J时,果实接触的曲率半径越大,其损伤敏感性越高,这与本文的试验结果是一致的。
表3 不同生长位置香蕉果指碰撞损伤敏感性均值的多重比较结果1)Table 3 Multiple comparisons of average bruise susceptibility of banana finger at various growth positions
2.2.1 果指表面区域对碰撞损伤敏感性的影响为研究果指表面区域对损伤敏感性的影响,从果指的凸起部分沿其生长方向等间距选取3个区域,分别记为近果柄端、果指中部和近果顶端,用低水平(0.52 J)的冲击能量对相同生长位置的香蕉进行碰撞,果指表面区域单因素试验结果如图7所示。
图7 外排香蕉果指不同表面区域的损伤敏感性Fig.7 Bruise susceptibility of different surface positions of the banana finger in the outer row
果指表面区域3种水平下的损伤敏感性的正态性检验显著,近果柄端水平下P=0.07;果指中部水平下P=0.23;近果顶端水平下P=0.60。方差齐性检验显著,P=0.14。
近果柄端、果柄中部和近果顶端的损伤敏感性分别为0.48、0.50和0.59 cm3/J,其损伤敏感性沿着果指生长方向逐渐增加。近果顶端的损伤敏感性与其他2个区域差异显著。果指中部与近果柄端的损伤敏感性无显著差异。
统计香蕉果指的形态发现,香蕉果指从果柄至果指中部区域弯曲较大,从果指中部至果顶几乎不弯曲或弯曲较小,说明果指的曲率半径从近果柄端到近果顶端依次增加,图6也很好地证实了这一点。先前的研究发现,水果的曲率半径越小,在低水平冲击下越易损伤[27]。然而,Abedi等[26]发现当冲击水平超过0.24 J时,撞击区域的曲率半径越大,水果的损伤程度越大,其损伤敏感性反而越高,本文的试验结果也很好地证明了这一点,即曲率半径对水果的碰撞损伤敏感性扮演了双重角色。Bugaud等[16]在研究贮藏条件与基因型对香蕉损伤敏感性影响的研究中发现,香蕉损伤敏感性与其硬度负相关。为了揭示造成香蕉果指不同区域损伤敏感性显著差异的机理,我们将在今后的研究中进一步证实。
2.2.2 冲击水平对碰撞损伤敏感性的影响 为了研究冲击水平对香蕉果指碰撞损伤敏感性的影响,选择相同生长位置的香蕉,对香蕉果指的近果顶端施加 3 种水平的冲击:低水平 (0.52 J)、中水平 (0.88 J)和高水平(1.17 J)。3种冲击水平下的损伤敏感性正态性检验显著,低水平下P=0.83,中水平下P=0.27,高水平下P=0.43;方差齐性检验显著,P=0.91。
损伤敏感性随冲击水平变化的结果如图8所示。从图8中可以看出,香蕉果指碰撞损伤敏感性在3种冲击水平下分别为0.52、0.95和1.35 cm3/J,随冲击水平的增加而显著增加。这个结果与先前其他水果碰撞损伤敏感性的规律是一致的。Van Zeebroeck等[4]用钟摆冲击的方式研究了番茄碰撞损伤的敏感性,发现被撞部位吸收的能量随冲击水平的提升而增加,即果实损伤的敏感性也在增加。Du等[28]通过自由下落试验与有限元分析的方法研究猕猴桃的碰撞损伤敏感性时也有类似发现。所以,在香蕉采收或采后处理时需尽量降低果实冲击,以保证果实质量。
图8 不同冲击水平下香蕉内排果指近果顶端的损伤敏感性Fig.8 Bruise susceptibility near apex of the banana finger in the inner row under different impact levels
正交试验结果如表4所示,方差分析如表5所示。
根据表4极差分析可知,影响香蕉果指损伤敏感性的主次顺序为:冲击水平>生长位置>果指表面区域;表5的正交试验方差分析结果显示,冲击水平、生长位置、果指表面区域对香蕉果实损伤敏感性影响极显著(P<0.01)。
表4 正交试验方案及结果1)Table 4 The scheme and result of orthogonal test
表5 正交试验方差分析Table 5 Variance analysis of orthogonal test
1)嵌套试验结果表明:香蕉果指的碰撞损伤敏感性从首梳节点到尾梳节点无显著性差异,而同梳节点嵌套的内外2排果指的损伤敏感性差异显著,从首梳至尾梳内排果指的损伤敏感性逐渐增加,而外排果指的损伤敏感性与此相反;多重比较结果显示,根据显著性水平,可以将节点嵌套排层的14种水平分为4组处理。
2)单因素试验结果表明:果指表面区域和冲击水平对其碰撞损伤敏感性影响极显著,其中近果顶端的损伤敏感性明显高于果指中部和近果柄端,而果指中部和近果柄端的损伤敏感性无显著差异;损伤敏感性随冲击水平的升高而增加。
3)正交试验结果表明:影响损伤敏感性各因素的主次顺序为冲击水平>生长位置>果指表面区域。
本文揭示了适收期内青香蕉果穗的梳层节点、嵌套的排层、果指表面区域和冲击水平对果实碰撞损伤敏感性的影响,初步分析了相关的原因,但对造成显著差异机理层面的解释还存在不足,将在后续的研究中进一步证实。