1、山东农业大学学报(自然科学版),2023,54(3):454-463VOL.54 NO.3 2023Journal of ShandongAgricultural University(Natural Science Edition)doi:10.3969/j.issn.1000-2324.2023.03.017基于基于 Fluent 核桃定向吸附装置设计与仿真试验核桃定向吸附装置设计与仿真试验刘梦迪,王应彪*,王周梅,张超宇,李九峰,周 丹西南林业大学 机械与交通学院,云南 昆明 650224摘摘 要要:为确保核桃果体在吸附摆放过程中姿态固定，基于核桃外形特征、力学特性设计一种核桃定向吸附装
2、置，建立一个核桃吸附过程的受力模型得出在吸附核桃过程中最低的气流速度 V=19.67 m/s、最小压强 P=55.87 kpa、最大压强 P=1.40 MPa。采用 Fluent 软件，对优化的吸口机构进行气流场仿真模拟；设计了吸口锥角、气孔尺寸、气孔间距3 因素 3 水平正交试验；试验结果表明吸口的最优参数为吸口锥角 60、气孔直径 2 mm、气孔间距 15 mm。为后续设计核桃定向吸附装置提供新的方法与理论。关键词关键词:核桃;定向吸附;机械设计中图中图法法分类号分类号:S229.1文献标识码文献标识码:A文章编号文章编号:1000-2324(2023)03-0454-10Design a
3、nd Simulation of the Oriented Adsorption Device toWalnut on FluentLIU Meng-di,WANGYing-biao*,WANGZhou-mei,ZHANG Chao-yu,LI Jiu-feng,ZHOUDanSchool of Machinery and Transportation/Southwest Forestry University,Kunming 650224,ChinaAbstract:In order to ensure that the attitude of the walnut fruit body i
4、s fixed in the adsorption and placement process,awalnut directional adsorption device is designed based on the appearance and mechanical characteristics of the walnut.Aforce model of the walnut adsorption process is established to obtain the The lowest airflow velocity V=19.67 m/s,minimumpressure P=
5、55.87 kPa,maximum pressure P=1.40 MPa.Fluent software is used to simulate the airflow field of the optimizedsuction mechanism;the horizontal orthogonal test of three factors is designed:the cone angle,pore size and pore spacing ofthe suction port is designed.The test results show that the optimal pa
6、rameters of the suction port are the cone angle of 60,the pore diameter of 2 mm and the hole spacing is 15 mm.It provides new methods and theories for the subsequent design ofwalnut directional adsorption device.Keywords:Walnut;oriented adsorption;machine design核桃是一种传统的营养食品，含锌、锰、铬等对人体有益的微量元素1。基于核桃的经
7、济价值，核桃在我国云南、新疆等地广泛种植，已成为西部地区乡村振兴重要的经济作物，目前，我国核桃种植面积和产量均位居世界前列2。核桃属于坚果，在食用之前需将其外壳去除，核桃破壳是核桃生产中的重要工序。由于核桃大小不同、形状不规则，且其外壳坚硬，在其破壳取仁过程中，容易造成核桃仁破碎而降低其品质。为保证核桃仁完整，长期以来，核桃破壳采用手工作业。这种生产方式，工作效率低、劳动强度大，随着我国农村剩余劳动力的减少，核桃机械破壳已成为核桃加工业发展的必然方向3-5。核桃在破壳加工前的预处理对后期核桃破壳效果和质量有很大的影响，若能通过科学的方法在破壳时对核桃的姿态进行定向处理，将很大程度提高对核桃破壳
8、的质量与效果。国内的学者们对核桃定向装置进行了各方面的研究，马秋成通过采用一种振动调姿定向破壳方法，在核桃破壳前对核桃的姿态进行定向处理，设计出定向破壳的装置，基本满足行业要求为后续核桃破壳机械的研发和改进提供参考6；刘海江通过采用差速带定向的方法设计出核桃破壳的定向装置，为云南漾濞核桃破壳方式提供了一种技术借鉴7；陈建在通过强迫振动下提出一种核桃纵径定向运输方法，为后续核桃定向破壳加工提供新方法8。王超通过一种赫兹接触理论计算出核桃破壳的最优接触面，提高收稿日期收稿日期:2022-10-25修回日期修回日期:2022-12-14基金项目基金项目:国家自然科学基金项目(52165038);云南
9、省教育厅科学研究基金项目(2022Y574);云南省农业联合专项(202101BD07001-062)第第 1 作者简介作者简介:刘梦迪(2000-),女,硕士研究生,研究方向为农业机械.E-mail:*通讯作者通讯作者:Author for correspondence.E-mail:第 3 期刘梦迪等:基于 fluent 核桃定向吸附装置设计与仿真试验455了核桃的破壳率9。国内对于气吸式、气盘式的吸附核桃定向破壳技术研究分析的少之又少，而在气吸式排种器上有相关研究。赵湛采用负压吸孔将种子吸住，当滚筒转至输送带上方合适位置时隔气板阻断负压源，种子在正压和自重作用下离开吸口，设计的气吸滚筒式
10、排种器为排种器的结构调整和工作参数优化提供了依据10；赵学观采用真空吸口，将定向好的玉米种子以一定方向摆放到纸带上，制成种子带，利用种子带播种装置将种子带铺放到土壤中，完成玉米的定向播种，构建了玉米种子吸附摆放试验台得出最优吸口参数为后期优化提供数据依据11。定向破壳装置提高核桃破壳的质量，从而提高核桃的整仁率和脱壳率，获得更高的产量。目前，对核桃姿态固定的机构还有所欠缺。本文以此为基础，在核桃姿态固定过程中，采用真空吸口，将核桃果体按一定方向摆放在挤压工位上进行破壳加工，为确保核桃果体在吸附摆放过程中姿态固定，设计一种核桃定向吸附装置。采用 Fluent 软件，对不同定向吸口机构进行气流场仿
11、真模拟分析比较各流场速度分布情况，得出最优定向吸口机构。分析吸口锥角、气孔直径、气孔间隙对核桃吸附的影响，确定影响吸口机构主要因素。经过正交试验得出影响吸口因素的最优参数。为后续设计核桃定向吸附装置提供新的方法与理论。1核桃物理参数的测量和计算核桃物理参数的测量和计算1.1试验材料和试验仪器试验材料和试验仪器本文选择云南大理种植的漾濞大泡核桃为试验对象。在试验中，随机选取 100 颗无明显残缺的核桃，对各物理参数进行测量：a 横径、b 纵径、c 棱径、尖端壳厚、缝合线壳厚、最凸处壳厚、底端壳厚如图 1。图图 1 核桃的基本尺寸核桃的基本尺寸Fig 1 Basic size of walnut每
12、个方向测量三次取平均值作为结果，对测量结果进行数据分析和统计。其结果如表1所示。表表 1 核桃的尺寸数据核桃的尺寸数据Table 1 Dimension data of walnut分布范围/mmDistribution range横径/mmTransverse diameter纵径/mmVertical diameter棱径/mmRidge diameter粒数百分比粒数百分比粒数百分比25,30）00003333%30,35）1818%006767%35,40）7676%4343%0040,45）66%5454%0045,50）0033%00最大值42.1548.2934.84最小值30.
13、5635.6928.39均值38.0640.4230.68三径均值36.39重量均值/g13.01壳厚/mm尖端缝合线最凸底端4.183.301.462.851.2核桃物理参数计算核桃物理参数计算在设计核桃吸附机构时，核桃的外形特征对吸附机构设计有很大影响，为准确的描述核桃的形456山东农业大学学报(自然科学版)第 54 卷状特征，引入参数球形率 Sp12。其中球形率的计算公式为：3/13/123/1366sssssssssssPltwlltwltwlS,(1)式中：Sp球形率；ls核桃纵径；ws核桃横径；ts核桃棱径。根据式13（1）计算核桃的球形率为 0.89，当 Sp0.85 时，为类球
14、形，当 Sp0.85 时，为椭球形14。则核桃为类球形。1.3核桃壳力学性能计算核桃壳力学性能计算核桃壳承受外界载荷时,当载荷增大到某一个值时，壳体会失去原来的形状,发生失稳破壳。壳体失稳的临界压力为：22132rhEPcr,(2)式中：Pcr临界压力；E弹性模量；泊松比；h核桃壳厚度；r核桃半径。通过式15（2）壳体临界压力的计算公式，得出对于同一种核桃，其中 E，r 一定时，核桃不同位置的壳厚不同失稳临界压力就不同，壳越厚临界压力就越大。取核桃的参数：E=0.18 GPa、泊松比=0.316由上面测得核桃直径、4 处壳厚结果带入式子（2）中计算，得 4 个部位的临界压力值分别为 Pcr1=
15、15.22 MPa、Pcr2=6.04 MPa、Pcr3=1.40 MPa、Pcr4=5.34 MPa。由于在吸附核桃的时候壳体不能发生失稳破裂所以选最小的临界压力值 Pcr3=1.40 MPa 作为吸附试验时真空泵不能超过的压力值。2装置整体结构与工作原理装置整体结构与工作原理核桃的外形呈现不规则状，常规的建模方法无法准确的还原真实特征，本文采用 3D 扫描的方法获取核桃模型，获取核桃模型的的步骤：先通过分区拍照，将照片导入软件，进行处理，合成 3D模型如图 2；工作时即沿负压方向吸附。核桃定向吸附装置与吸附性能息息相关且影响到核桃姿态固定和吸附的效果；在吸附过程中，核桃在负压的作用下由于气
16、流不均匀极易发生偏转，影响核桃的姿态固定。本文研究设计的定向吸附机构能够使核桃姿态固定更稳定，为后续设计核桃定向吸附装置提供新的方法与理论。图图 2 核桃核桃 3D 扫描模型扫描模型Fig.2 3D scanning model of walnut2.1整体结构设计整体结构设计该装置主要部件包括真空泵、真空压强表、真空管、管卡、定向吸头、定向凹槽、吸盘、机架等，其结构组成如图 3 所示。其中提供负压环境的装置为真空泵，真空泵通电与真空管连接定向吸头进行抽取是坏境呈现负压环境；所显示负压值可以通过真空压强表随时观察，控制压强值在一定第 3 期刘梦迪等:基于 fluent 核桃定向吸附装置设计与仿
17、真试验457水平。定向吸头与定向凹槽配合使用，来完成核桃定向吸附的动作；管卡起到固定真空管的作用；吸盘在机架的四个支架下面是因为在真空泵进行工作的时候会产生震动，此时吸盘吸住地面就减少震动，使得定向吸附过程更加平稳可靠。1.吸盘 2.机架 3.真空压强表 4.真空泵 5.管卡 6.真空管 7.定向吸头 8.定向凹槽图图 3 核桃定向吸附机构核桃定向吸附机构Fig.3 Walnut directional adsorption mechanism2.2工作原理工作原理核桃以吸附方向放置在图 3 定向凹槽中，此时核桃的姿势方向就被固定了，呈现核桃最凸面朝上姿态就开始进行吸附；工作时，打开真空泵，真
18、空泵通过连接真空管与定向吸头来制造负压环境，此时核桃在定向凹槽中，由于真空泵提供的负压超过吸附起核桃的最小压强，此时核桃被吸起，随着负压不断的增大，增大到一定程度最终核桃被完全吸附在定向吸头的凹槽里，此时就完成了核桃的定向吸附过程如图 4。图图 4 核桃定向吸附过程图核桃定向吸附过程图Fig.4 Walnut directional adsorption process diagram2.2.1 核桃工作过程受力分析 在吸附过程中，设核桃直径为 d，密度为m；类球体为中心对称的旋转体，故核桃所受的升力 H=0，其他力有气动阻力和浮力，其计算公式如下17。气动阻力与气流方向一样所以是向上，浮力是
19、向上，受力分析如图 5。图图 5 核桃的简化模型核桃的简化模型Fig.5 A simplified model of walnut458山东农业大学学报(自然科学版)第 54 卷气动阻力公式为：2228121dVCAVCFddd(3)浮力公式为：gdFa361(4)向下有重力，其计算公式为：gdGm361(5)核桃被吸起的条件是：向上合力大于等于向下合力，即：GFFad(6)其临界状态为：GFFad(7)将式（1）、式（2）、式（3）代入式（5）则有：gdgdVdCmd3322616181,(8)式中：Fd气动阻力；Cd阻力系数；流体密度；V自由流速度/ms-1；A物体投影面积/m；d核桃直径
20、/mm；Fa浮力/N；G重力/N。由式（6）得出核桃被吸起时临界速度 V 为：gdCVmd34，(9)SGP。(10)对于本课题的研究情况，取 Cd=0.48，m=470 kg/m3，=1.2 kg/m3，g=9.8 m/s2，d=36.39 mm 带入式（7）得 V=19.67 m/s，吸附的受力面积是类半球所以 S=1/2（4 r2）带入式（8）P=55.87 kPa。从核桃在悬浮状态时的受力平衡方程式中，只要气流速度大于平衡时的临界速度 V=19.67 m/s时，核桃就会受到一个向上的加速度，此时需要的最低压力值为 P=55.87 kPa，从而就能被吸头所吸附。2.3关键部件设计关键部件
21、设计2.3.1 吸头外形尺寸设计 吸附机构是核桃定向吸附装置的关键部件之一，为保证核桃能顺利实现姿态固定，定向吸头和定向凹槽都是以核桃的外形尺寸进行设计，考虑到在吸附过程中核桃会有轻微转动定向凹槽应比本身核桃尺寸 b 多出 12 mm，为使气流更加均匀、定向效果更好，设计的吸口下端直径比凹槽要大、吸口锥角要平缓、吸头高度的尺寸要合理。本研究的设计的吸附机构的吸孔结构、形状、尺寸使其定向吸附合理容易实现如图 6。图图 6 吸头尺寸图吸头尺寸图Fig.6 Suction head dimensional drawing吸头形状尺寸取决核桃的形状尺寸要求为：L12.5b(11)L21/2a(12)L
22、3b(12)(13)第 3 期刘梦迪等:基于 fluent 核桃定向吸附装置设计与仿真试验459L42b(14)r1090，(15)式中：L1吸头高度/mm；L2吸头凹槽高度/mm；L3吸头凹槽直径/mm；L4吸头下端直径/mm；r1吸头锥角/；根据测的核桃尺寸则 L1为 80 mm，L2为 40 mm，L3为 42 mm；L4为 100 mm；r1为 30。2.3.2 吸头凹槽优化 观察核桃外形发现，核桃外轮廓不太平整有沟壑等，为满足吸附定向核桃的要求，设计了四种不同的凹槽造型。查阅农业机械设计手册其中气孔采用传统的圆型气孔，且均匀分布在凹圆槽上如图 7a；在 a 的基础上对吸孔进行优化，改
23、变传统的圆型吸孔，采用梯田扇形交错的气孔，两梯田间距均匀如图 c；在 a 的基础上，考虑到核桃表面粗糙,对凹圆槽进行优化,增加两个十字交叉的抓夹形状来增大摩擦这样更有利于吸附核桃，也防止在吸附过程中滑落如图 b；在图 c的基础上对凹圆槽进行优化，梯田扇形交错的气孔与两个十字交叉的抓夹形状结合尺寸不变如图 d。图图 7 不同定向吸头不同定向吸头Fig.7 Different directional suction heads3吸头特性仿真分析吸头特性仿真分析3.1模型建立及试验条件模型建立及试验条件在 UG 中建立吸口模型，导入到 Fluent 中，进行 Geometry 填充实体，进行 Mes
24、h 网格划分，对吸口附近区域进行网格加密处理最后进行计算。试验条件：Fluent 软件在相对压力-60 kPa 条件下，以吸口下平面为 Z=0 为基准面，选取基准面以下 20 mm，即 Z=-20 截面为观察面；在仿真模拟时对吸口模型作如下假设：1.空气为不可压缩理想气体。2.室温为恒 25。3.吸孔壁面施加无滑移边界条件。3.2仿真试验仿真试验对图 7a、b 进行初步单因素仿真试验，吸口锥角 30气孔直径为 2 mm、气孔间距 10 mm 进行初步仿真比较。Z=-20 截面的速度云图如图 8a、b 所示。由图 8a 所示，吸口流场的气流速度中心处高且气流速度分布均匀，吸口流场的气流速度分布面
25、积大，中心边缘处气流速度低；由图 8b 所示，中心气流速度高分布不均，中心边缘气流速度低且气流速度扩散不均，气流速度分布面积小，这主要是由于气流在水平方向上由内侧向外侧聚集。图图 8 气流速度云图气流速度云图Fig.8 Air velocity nephogram460山东农业大学学报(自然科学版)第 54 卷对图 7c、d 进行单因素仿真试验，吸口锥角 30梯田扇形角度 90间距为 1 mm 进行初步仿真比较。速度云图如图 9a、b 所示，由图 9a 所示，吸口流场的气流速度中间为 0，边缘气流场速度一圈高且边缘速度不均匀；图 9b 所示，吸口流场的气流速度中间为 0，边缘气流场外缘一圈高且
26、不均匀，吸盘外缘其他气流速度为 0。图图 9 气流速度云图气流速度云图Fig.9 Air velocity nephogram对四个仿真结果进行分析，图 8a、b 比较，a 的气流速度比 b 的气流速度均匀且气流速度流场面积大；图 9a、b 气流场中间速度为 0，周围气流不均匀。综上，选择图 7a 为最优吸头机构。3.3单因素仿真试验单因素仿真试验对图 7a 机构进行单因素试验，核桃所受吸附力与气孔参数、内外压强差、气流速度等都有关系，以吸口锥角、气孔直径、气孔间距为因素，在相对压力-60 kPa 条件下，以吸口下平面为 Z=0 为基准面，选取基准面以下 20 mm，即 Z=-20 截面为观察
27、面，以吸孔锥角 30，气孔直径 2 mm，气孔间距为 10 mm 为参考条件，分别进行单因素对比试验18,19。Z=-20 mm 截面云图如图 10 所示，由图 10b 分析，吸口锥角为 60时中间气流速度与 30中间气流速度减小，气流分布面积减小且相对均匀，气流速度由中间向外部逐渐减小；由图 10c 分析，当气孔直径为 4 mm 时中间气流速度增大，气流分布面积增大且速度相对均匀，气流速度由中间向外部逐渐减小；由图 10d 分析，当气孔间距为 20 mm 时，中间气流速度变小，气流分布面积变小且气流速度不均匀。可见吸口锥角对气流速度影响不太大，吸口直径越大中间气流速度越高气流分布面积增大且均
28、匀；气孔间距越大，气流速度减小且气流速度不均匀，气流分布面积减小。a.参考组；b.吸口锥角 60；c.气孔直径 4 mm；d.气孔间距 20 mm图图 10 单因素试验气流速度云图单因素试验气流速度云图Fig.10 Single factor test airflow velocity nephogram第 3 期刘梦迪等:基于 fluent 核桃定向吸附装置设计与仿真试验4613.4试验因素与试验指标试验因素与试验指标由单因素试验确定了吸口锥角、气孔直径、气孔间距 3 个因素对气流场的影响，为了进一步获得各因素的最优值；考虑到各因素之间的交互作用，设计三因素三水平正交试验如下表 2。在相同的
29、参数下，流场各处的 X 轴速度和 Z 轴速度比值不同；在不同参数下，流场各处的 X 轴速度和 Z 轴速度比值也不同，吸口气压决定了气流速度，气流速度决定了吸口的吸核桃的性能，为了探究吸口最佳性能，以吸口下平面中心为原点，在吸口流场 Z=-20 mm 截面 X 轴设置 4 个监测点分别为 a（5，5）、b（5，-5）、C（-5，5）、d（-5，-5）并确定这 4 个监测点水平速度与垂直速度比值20,21的平均值 K 及平均速度 F 为指标。表表 2 试验因素及水平表试验因素及水平表Table 2 Test factor and level table水平Level因素 Factor吸口锥角 A/
30、Cone angle of suction port气孔直径 B/mmAir hole diameter气孔间距 C/mmAir hole gap1302102453153604203.4.1 正交试验结果分析 考虑到各因素的交互作用，选择正交表 L27 试验和方差分析见表 3-5。表表 3 试验结果分析计算表试验结果分析计算表Table 3 Analysis and calculation table of test results序号No.各因素水平 Factors and levelsKFAB(AB)1(AB)2C(AC)1(AC)2(BC)1(BC)21111111111-2.9552
31、21111222220.2103631111333330.5106941222111235.3360.551222222311.0593.361222333122.1363.471333111320.3266.88133322213-1.347391333333210.0370.7102123123112.7344.611212323122-1.7482.612212331233-0.5589132231123332.7258.914223123111-0.4504152231312225.3306.5162312123228.2192.3172312231333.5359.118231231
32、2111.4291.219313213211-0.2231.5203132213221.1433.121313232133-0.2442.222321313233-4.9207.6233213213110.6302.624321332122-0.1257.125332113222-0.6134.6263321213330.2204.3273321321110.7195.8k5.2-16.55.710.610.84.59.32.9T23k21.211.69.222.23.213.90.44.814.8k-3.412.47.3-0.59.2-1.718.18.95.3极差24.613.42.722
33、.77.415.617.74.511.9f5085.3 5180.6 4267.54261.72548.833063420.33091.53385.7T10 822.3f3328.2 3153.9 3268.13266.64388.1 3793.33862.83636.53472.5f2408.8 2487.8 3286.732943885.437233539.44094.43964.1极差2676.5 2692.8999.41206.51839.3487.34421002.9578.4462山东农业大学学报(自然科学版)第 54 卷由表 4 可知，吸口锥角对指标 K 影响极显著，气孔直径对指
34、标影响显著。分析这 3 个交互作用发现吸口锥角和吸孔直径影响显著，其他交互影响不显著。以 K 值为指标，以因素方差显著性确定因素间主次顺序为气孔锥角气孔直径气孔间距，则气孔锥角、气孔直径影响气流稳定性较大，气孔间距影响气流稳定性较小。K 值越小，越不容易发生偏转选取三个因素中的最优水平为 A3B1C2即当吸口锥角 60、气孔直径 2 mm、气孔间距 15 mm。气流稳定性最好。表表 4 K值方差分析值方差分析表表Table 4 K-value analysis of variance table方差来源Soruce of variation离差平方和Sun of saquares of dev
35、iations自由度Degree of freedom均方值Quadratic meanF 值显著性SignificanceA34.6217.33.1*B12.626.31.1-C3.421.70.3-AB39.249.81.8*AC34.248.61.6*BC10.242.60.5-误差44.185.5-总和138.426-由表 5 可知，吸口锥角对指标 F 影响极显著，气孔直径对指标影响显著。以 F 值为指标，以因素方差显著性确定因素间主次顺序为吸口锥角气孔直径气孔间距。综上所述，吸口锥角越大，吸口的锥形腔越平缓，气压分布均匀，水平方向不易发生偏移。气孔直径越大，通过气孔的气流就越多，气流
36、速度越大。气孔间距对速度影响较次，由于吸口为凸圆槽，气孔间距越小分布越密，增大了气流向中心汇聚的强度，造成了气流速度小且分布不均。表表 5 F值方差分析表值方差分析表Table 5 F-value analysis of variance table方差来源Soruce of variation离差平方和Sun of saquares of deviations自由度Degree of freedom均方值Quadratic meanF 值显著性SignificanceA410975.922054882.8*B4371252218562.52.9*C200823.42100411.71.3-A
37、B144019.843600454.8*AC27229.746807.40.09-BC77719.3419429.80.3-误差59759.8874680-总和1357652.926-4结结 论论本文采用理论分析与仿真相结合的方法，为防止核桃在吸附过程中由于负压太大而裂开，根据核桃壳的受力算出壳体临界压力 Pcr3=1.40 MPa，通过建立核桃在吸附过程的受力模型，得出了吸附起核桃临界气流速度为 V=19.67 m/s，吸附最低气压 P=55.87 kPa 和吸附最大气压 P=1.40 MPa 为核桃定向吸附特性试验提供理论依据。本文设计了一种核桃定向吸附机构，通过分析工作原理，以云南漾濞核
38、桃为研究对象，对吸头的各个参数进行了设计保证了定向吸附的效果。对优化的四种吸头进行 Fluent 仿真分析，最终确定了吸口采用图 7a 定向吸头机构。根据单因素仿真试验确定了影响气流速度的三个关键因素（吸口锥角、气孔直径、气孔间距）。由正交仿真试验发现,吸口锥角对速度之比 K 的均值影响极显著，吸口锥角、气孔直径对气流速度影响显著。根据正交试验表分析得出最优参数为吸口锥角 60、气孔直径 2 mm、气孔间距 15 mm。为后续搭建试验台提供理论依据。第 3 期刘梦迪等:基于 fluent 核桃定向吸附装置设计与仿真试验463参考文献参考文献1张 亭,杜 倩,李 勇.核桃的营养成分及其保健功能的
39、研究进展J.中国食物与营养,2018,24(7):64-692张永成,马佳乐,唐玉荣,等.我国核桃初加工现状与分析J.食品工业,2020,41(7):198-2023Liu MZ,Li CH,Cao CM,et al.Walnut fruit processing equipment:academic insights and perspectives J.FoodEngineering Reviews,2021:1-364刘甲振,耿爱军,栗晓宇,等.国内外核桃破壳技术与装备研究进展J.中国果树,2019(1):16-20,265Zhang H,Shen LY,Lan HP,et al.Mec
40、hanical properties and finite element analysis of walnut under different crackingparts J.International Journal ofAgriculturaland Biological Engineering,2018,11(6):81-886马秋成,陈 强,卢安舸,等.核桃定向破壳装置设计及试验J.江西农业大学学报,2022,44(2):473-4857刘海江,王应彪,罗思蓝,等.云南漾濞核桃定向破壳设备的设计与研究J.山东农业大学学报(自然科学版),2020,51(3):487-4948陈 健,王
41、应彪,刘海江,等.强迫振动下核桃定向输送方法及仿真试验J.山东农业大学学报(自然科学版),2021,52(6):1035-10419王 超,郑甲红,霍启新,等.基于赫兹接触理论的核桃破壳分析与试验研究J.西南农业学报,2020,33(12):2967-297410 赵 湛,李耀明,陈 进,等.气吸滚筒式排种器吸种过程的动力学分析J.农业工程学报,2011(7):112-11611 赵学观,徐丽明,王应彪,等.基于 Fluent 与高速摄影的玉米种子定向吸附研究J.农业机械学报,2014,45(10):103-10912 Geng J,Wang S.Gao Z,et al.Influence o
42、finitial orientations on potato conveyor trajectories by machine vision J.Computers and Electronics inAgriculture,2019,163(38):10483813 崔 涛,刘 佳,杨 丽,等.基于高速摄像的玉米种子滚动摩擦特性试验与仿真J.农业工程学报,2013,29(15):34-4114 刘文政,何 进,李洪文,等.基于离散元的微型马铃薯仿真参数标定J.农业机械学报,2018,49(5):125-135,14215 徐芝纶.弹性力学（上,下）M.北京:高等教育出版社,199016 刘
43、明政,李长河,张彦彬,等.柔性带剪切挤压核桃破壳机理分析与性能试验J.农业机械学报,2016,47(7):266-27317 庄 森,田盈辉,蒋 晶,等.真空气吸式播种器吸嘴流场的研究J.农机化研究,2007(11):83-8618 韩占忠,王 敬,兰小平.Fluent 流体工程仿真计算实例与应用M.北京:北京理工大学出版社,200419 于 勇,张俊明,姜连田.Fluent 入门与进阶教程M.北京:北京理工大学出版社,200820 王延耀,李建东,王东伟,等.气吸式精密排种器正交试验优化J.农业机械学报,2012,43(s1):54-5821 张 敏,吴崇友,张文毅.吸盘式水稻育秧播种器吸孔气流场仿真分析J.农业工程学报,2011,27(7):162-167