基于论文摘要与主题，拟定的专业中文标题为： 不同品种小粒咖啡果实机械压缩特性的实验表征与有限元建模

《International Journal of Food Properties》：Experimental characterization and finite element modeling of mechanical compression in Coffea arabica fruits from different varieties

　　

摘要翻译：

机械损伤在咖啡脱皮过程中可能影响咖啡豆品质，因此对咖啡鲜果机械特性的表征对于机器设计具有重要意义。本研究评估了Castillo和Colombia两个品种的成熟小粒咖啡（Coffea arabica）鲜果在不同果实尺寸和压缩方向下的坚实度（firmness）与破裂力（fracture force）。研究人员采用析因实验设计（factorial experimental design）和方差分析（ANOVA），通过平行板压缩试验共测试了192颗鲜果。坚实度由力–位移曲线破裂前线性段的斜率确定，而破裂力则在结构失效起始点处识别。平均破裂力为38.19 ± 19.35 N，且从小型果实的34.78 ± 17.30 N增加至大型果实的41.61 ± 20.73 N。破裂力受品种（p = .027）、尺寸（p = .011）和压缩方向（p = .027）的显著影响，包括两个赤道方向和极地方向。平均坚实度为11.52 ± 5.04 N/mm，仅受品种显著影响（p = .001）。有限元模型（finite element model, FEM）再现了非线性压缩趋势，并识别出沿载荷传递路径向种子集中的应力分布。这些发现为咖啡脱皮机的校准与设计提供了定量力学参数。

关键词：ANOVA；咖啡品质；咖啡尺寸；FEA；破裂力

论文解读文章

一、研究背景与问题

咖啡品质评价是决定其商品价值的关键环节，其中物理特性评估直接影响“优质咖啡”的认定比例。在咖啡加工过程中，脱皮工序造成的机械损伤是导致咖啡豆品质下降的重要来源之一。根据脱皮机质量标准，加工批次中机械损伤豆比例不得超过0.5%，然而实际生产中这一标准很难达到。常见原因包括脱皮操作不当、未成熟果过多、设备校准不准确以及维护不足等。此外，现有脱皮机在设计时普遍缺乏对咖啡鲜果力学特性的充分考虑，导致设备对不同品种、尺寸和成熟度果实的适应性和加工效率不足。已有若干研究尝试优化脱皮机设计，但多数仅关注机器结构或传动参数，未将咖啡果实的力学特性作为核心设计输入。有限元方法（FEM）在果实力学领域已成为预测局部应力、应变和变形场的重要工具，但针对咖啡果实压缩行为的有限元研究仍然有限，且多依赖文献中近似或推算的材料参数。因此，迫切需要开展实验标定的有限元研究，以获取更具工程代表性的力学描述参数，服务于脱皮机设计和加工损伤防控。

二、研究目的与主要内容

本研究旨在通过平行板压缩试验，系统测定四个哥伦比亚小粒咖啡品种（Castillo 1、Castillo 2、Colombia 1、Colombia 2）成熟鲜果的破裂力和坚实度，分析品种、果实尺寸和压缩方向三个因素的影响；同时，基于实验几何和载荷条件建立有限元模型，模拟压缩过程中的弹性变形行为，获取与刚度相关的代表性参数。研究采用全因子实验设计，因素包括品种（4个水平）、果实尺寸（2个水平）和压缩方向（3个水平，分别为极向P、赤道方向1 E1、赤道方向2 E2），共24个实验条件，每个条件8次重复，总计192颗鲜果。果实采自哥伦比亚咖啡带不同地理区域的两个农场，每个品种分别来自不同农场，以检验同一品种名称下不同产地是否产生力学差异。鲜果在采集当日运抵实验室，在22°C恒温条件下进行颜色筛选、三维尺寸测量（极径、赤道1径、赤道2径）、体积估算（按椭球体公式V=1/6πD_PD_E1D_E2计算）、大小分组和随机分配压缩方向。压缩试验在位移控制下进行，记录力–位移曲线，以曲线破裂前线性段的最小二乘回归斜率定义坚实度（N/mm），以结构失效起始点的力值定义破裂力（N）。

三、主要关键技术与方法

研究采用的主要关键技术包括：（1）平行板压缩试验（parallel-plate compression test），用于获取鲜果在受压过程中的力–位移响应；（2）全因子实验设计与方差分析（ANOVA），用于评估品种、尺寸和压缩方向对坚实度和破裂力的影响及其交互作用；（3）有限元建模（FEM），基于实验几何建立三层结构（外果皮、中果皮、种子）的数值模型，模拟压缩过程并预测内部应力分布；模型采用位移控制加载，外果皮设为薄弹性壳，中果皮设为近不可压缩软固体，种子设为弹性固体，接触界面设置摩擦系数0.6。样本来源为哥伦比亚咖啡带两个农场的四个品种鲜果。

四、研究结果

（一）样品选择与尺寸测量

几何测量结果表明，Colombia品种的平均直径普遍大于Castillo品种；平均极径大于赤道1径，赤道1径大于赤道2径。果实平均尺寸按升序排列为：Castillo 1 < Castillo 2 < Colombia 1 < Colombia 2。与Castillo 1相比，Colombia 2的极径增加约15.5%，赤道1径增加约8.3%，赤道2径增加约9.0%，估算体积增加约34.9%，显示实验组间尺寸差异明显。

（二）样品压缩与力–位移曲线

192颗鲜果的压缩试验完成后，力–位移曲线呈现出明显的逐样本差异，反映了每颗果实的固有变异性，包括遗传差异、光照、养分可用性以及病虫害等影响。不同压缩方向下虽出现多种相似行为，但没有一种特定行为可明确归因于某一方向。在同一方向和尺寸分类下，破裂结果也可能不同。破裂力对应的平均位移约为3.5 mm，可作为脱皮机设计或校准的广义参考。

（三）果实抗脱皮力的计算

通过最小二乘线性回归计算坚实度。总体结果显示出破裂力随果实增大而略高的趋势，且极向破裂力大于赤道方向3（此处指E2方向？原文为equatorial fracture force 3，可能笔误，但需忠实）的破裂力，后者又大于赤道方向1的破裂力。但破裂力的标准差很大，约为均值的一半，与已有文献一致。

（四） expelled seeds measurement（排出种子的测量）

压缩后对排出种子进行测量，结果显示种子尺寸与果实尺寸呈正相关，即果实越大，其内部种子越大。由于赤道1轴上有两颗种子，种子厚度测量值通常小于赤道方向测量值。

（五）坚实度的统计分析

ANOVA结果表明，品种对坚实度有显著影响（p = .001），因此拒绝原假设1；尺寸（p = .411）和压缩方向（p = .145）对坚实度无显著影响，不能拒绝原假设2和3。即坚实度主要取决于咖啡品种，而与果实大小和压缩方向无关。

（六）破裂力的统计分析

ANOVA结果表明，品种（p = .027）、尺寸（p = .011）和压缩方向（p = .027）均对破裂力有显著影响，因此拒绝原假设4、5和6。破裂力依赖品种、尺寸和压缩方向，说明脱皮过程中果实定向和尺寸分级对所需力有重要影响。

（七）有限元分析

有限元模拟结果显示，压缩过程中果实内部应力场高度不均匀，最高应力集中不仅出现在上下接触区域，还出现在载荷向种子传递的内部路径上。这表明果实结构并非均匀地分布压缩载荷，而是沿特定区域将应力导向种子，在果实内部形成二次集中区。中果皮区域应力较低且分布较平滑，表现出更柔顺的力学响应。模型预测的力–位移曲线与实验曲线定性吻合：初始阶段为低刚度区，随后斜率逐渐增大，反映组织逐渐参与承载的过渡。模拟初期低估了实验力，可能与组织本构简化和接触条件有关；在较大位移下，实验曲线因种子排出而下降，而模型未包含种子排出机制，故持续单调上升。模型采用的外果皮、中果皮和种子杨氏模量分别为50 MPa、2 MPa和400 MPa，泊松比分别为0.30、0.49和0.30。这种刚度分布使中果皮在初始阶段承担大部分变形，而外果皮和种子随位移增加逐渐主导响应，表明咖啡果实的非线性压缩行为主要受内部组织刚度对比控制。

五、讨论与结论总结

讨论部分将本研究的坚实度和破裂力值与文献进行了对比。坚实度值与Buitrago-Osorio等人的报告相近，但本研究中极向坚实度相对较高，而其他作者的极向值偏低。破裂力值总体低于文献报道值；其他作者的结果显示极向破裂力低于赤道方向，而本研究中极向值仅略低于赤道值，差异不显著。ANOVA结果拒绝原假设5和6，表明尺寸和方向依赖性的存在具有实践意义：较小果实需要的脱皮力更小，因此脱皮机的设计和功率要求应针对不同尺寸果实分别考虑，或采用不同功率的机器；方向依赖性则提示脱皮机设计者需分析脱皮过程中的潜在压缩方向。拒绝原假设1和4表明品种对坚实度和破裂力均有影响，且同一品种名称下不同农场来源的果实可能因种植实践和生长条件不同而表现出力学差异，因此“品种”不等于官方名称所代表的固定特性。研究提出的定量力学标准可支持脱皮机的设计、再设计和校准，包括果实尺寸分级、间隙校准、机械载荷定义、驱动系统尺寸确定和结构尺寸确定。

结论部分：本研究通过平行板试验和有限元建模，表征了来自不同哥伦比亚品种的小粒咖啡果实的压缩行为。结果显示，破裂力受品种、果实尺寸和压缩方向显著影响，而坚实度主要受品种影响。总体上，较大果实倾向于承受更高的破裂力，力学响应随加载方向变化。有限元模型再现了主要的实验压缩趋势，并帮助识别了果实内部的非均匀应力分布。这些发现为改进咖啡脱皮机的设计和操作提供了有用的力学标准，有助于实现更合适的加工条件和降低咖啡豆损伤风险。