Effects of short-term exposure to high temperature at larval stage on the development and fecundity of Assara inouei
-
摘要:目的
明确幼虫期短时高温对井上蛀果斑螟 Assara inouei 存活、生长发育及生殖的影响。
方法将井上蛀果斑螟初孵幼虫和老熟幼虫在32、35、38、41和44 ℃条件下，分别处理1、2和4 h，测定其存活率；将井上蛀果斑螟初孵幼虫和老熟幼虫在35、38、41 ℃条件下，分别处理1、2和4 h后，研究其发育历期、化蛹率、羽化率、雌雄虫寿命、单雌产卵量及子代卵孵化率等。
结果随着温度的升高和处理时间的延长，井上蛀果斑螟幼虫存活率逐渐减小。井上蛀果斑螟初孵幼虫在35 ℃处理下发育历期缩短，38和41 ℃处理下发育历期延长。老熟幼虫经短时高温处理后，随着温度的升高及处理时间的延长，化蛹率、蛹质量、蛹羽化率、产卵期、雌虫寿命及产卵量均逐渐减小，蛹历期逐渐延长。
结论井上蛀果斑螟初孵幼虫对短时高温的适应性高于老熟幼虫，短时高温对井上蛀果斑螟老熟幼虫的存活、生长发育及生殖有显著影响。本研究揭示了井上蛀果斑螟幼虫对高温的适应性，可为井上蛀果斑螟种群动态的预测预报提供科学依据。
Abstract:ObjectiveTo investigate the effects of short-term exposure to high temperature at larval stage on the survival rate, development and fecundity of Assara inouei.
MethodNeonate larvae or mature larvae of A. inouei were exposed to high temperatures of 35, 38, 41 and 45 ℃ for 1, 2 and 4 h, respectively, and the survival rate was measured. After neonate larvae or mature larvae of A. inouei were exposed to high temperatures of 35, 38 and 41 ℃ for 1, 2 and 4 h, respectively, biological parameters of developmental duration, pupation rate, emergence rate, male and female longevity, number of eggs laid per female and egg hatching rate were investigated.
ResultThe survival rate of neonate larvae decreased gradually with the increase of temperature and treatment time. The developmental duration of neonate larvae was shortened by 35 ℃ treatment, but the developmental durations were prolonged by 38 and 41 ℃ treatments. After mature larvae were exposed to high temperature for short term, the pupation rate, pupal weight, emergence rate, oviposition period, female longevity and number of eggs laid per female decreased gradually, and pupal duration prolonged gradually with the increase of temperature and treatment time.
ConclusionNeonate larvae have higher adaptability to short-term high temperature than mature larvae. Short-term exposure to high temperature has significant effects on the survival, development and fecundity of mature larvae of A. inouei. The results reveal the adaptability of A. inouei larvae to high temperature, and provide a scientific basis for predicting the population dynamics of A. inouei.
-
Keywords:
- Assara inouei /
- short-term high temperature /
- survival /
- development /
- reproduction
-
玉米已成为我国种植面积最大的粮食作物，为了满足玉米机械化种植高精度与高效率的要求，高速精量排种技术是玉米播种装备发展的重要方向[1-3]。
高速精量排种器是实现机械化高速精量播种的核心工作部件，其排种性能直接影响播种质量[4-7]。按工作原理可将排种器分为机械式[8-9]和气力式[10-12]两种。机械式精量排种器因其结构简单、可靠性好、制造和使用成本低等优点，在我国中小规模的农业生产中具有广泛的应用前景[13-17]。但随着排种器作业速度的提高，机械式精量排种器在高速作业时，难以实现稳定的充种和清种，从而导致漏播和重播现象严重，无法发挥其优点。
为提升机械式排种器在高速作业时的充种性能，Ma等[18]针对现有机械式排种器作业速度慢、对种子适应能力差的问题，设计了一种取种滑块主动夹取种子。Chen等[19]为提高勺轮式排种器的高速作业性能，利用离散元法(Discrete element method and multi-body dynamics, DEM)对种勺的结构参数进行了优化。都鑫等[20]设计一种扰种条增加充种区种群活跃度、降低种群内摩擦力，从而提高了机械式排种器的充种率和对速度的适应性。在高速清种性能提升研究方面，唐汉[21]在排种盘上设计一种振动凹凸面结构，利用机械振动与清种毛刷配合进行复合清种，提高清种效果。李玉环等[22]采用双侧清种方法，有效降低了漏清率，对高速条件下的过清率也有所改善。上述研究取得了一些研究成果，但是主要侧重于优化型孔、增设促充和改进清种结构来提高排种器的工作性能，对于高速作业条件下，种子无序堆积、充种姿态各异的状态并未改变，这将会影响排种器的充种与清种过程，从而造成机械式排种器在高速作业时作业性能下降。
为此，本文设计一种基于分级限姿导引侧充−径向突变辅助清种的双列限姿导引侧充式玉米高速精量排种器。该排种器可通过分离种群并规划种子运动轨迹的方法，有效降低种子无序堆积、充种姿态各异对充种过程的干扰，同时借助径向突变辅助清种结构，提升排种器在高速作业条件下的清种效率。根据种子的物理特性和充种阶段种子的填充姿态，结合对排种器工作过程的理论分析，进一步确定关键部件的基本结构参数。采用离散元法与多体动力学(Discrete element method and multi-body dynamics, DEM-MBD)耦合仿真试验和响应面优化方法对充种型孔的关键结构参数进行优化，并通过开展台架试验，验证仿真试验的可靠性和排种器高速作业时的作业速度适应性，以期为机械式高速精量排种器的研发提供参考。
1. 排种器结构与工作原理
双列限姿导引侧充式玉米高速精量排种器总体结构如图1所示，主要由前壳、左限姿板、轴承、轴端螺母、左排种盘、左护种盘、外壳、Π形入种管、导种盘、右护种盘、右排种盘、传动轴、右限姿板、后壳、型孔插件、压缩弹簧等部件组成。
![]()
图 1 排种器结构示意图1：前壳，2：左限姿板，3：轴承，4：轴端螺母，5：左排种盘，6：左护种盘，7：外壳，8：Π形入种管，9：导种盘，10：右护种盘，11：右排种盘，12：传动轴，13：右限姿板，14：后壳，15：型孔插件，16：导种槽，17：压缩弹簧，18：充种型孔。Figure 1. Structure diagrams of seed metering device1: Front shell, 2: Left posture limiting plate, 3: Bearing, 4: Shaft end nut, 5: Left seed metering plate, 6: Left seed protection plate, 7: Shell, 8: Π shaped seed intake tube, 9: Seed guiding tray, 10: Right seed protection plate, 11: Right seed metering plate, 12: Drive shaft, 13: Right posture limiting plate, 14: Back shell, 15: Orifice plug, 16: Seed guiding channel, 17: Compression spring, 18: Seed filling hole.如图1b所示，左、右排种盘分别装有18个均匀分布的型孔插件，通过压缩弹簧分别卡入左、右护种盘中，共同组成充种型孔；把左、右护种盘固定在外壳上后，传动轴依次穿过右排种盘、导种盘和左排种盘后用轴端螺母固定，通过轴承连接到前后壳体上；前后壳体上分别固定左、右限姿板，与排种盘和护种盘共同组成种子的分级限姿引导通道，从而实现分级限制引导种子姿态的目的。
种子从Π形入种管进入前、后两对称种室后，依次经过充种、清种、递种、携种、投种这5个阶段。充种阶段可以分为2个过程，在过程一中，种子经过限姿板、排种盘、型孔插件以及护种盘共同组成的限姿通道以合适的姿态到达预充种位置，完成分级限姿引导过程；在过程二中，随着排种盘的转动，预充种位置的种子侧充进入护种盘与型孔插件共同组成的充种型孔内，完成整个充种过程。在清种阶段中，种子随排种盘转动到清种区域时，在重力作用以及型孔插件的径向突变产生的摩擦力的辅助作用下，完成快速清种，清种完成后，充种型孔内仅余单粒种子，完成清种过程。在递种阶段中，单粒种子随排种盘转动到护种盘的递种口处时，在递种口处种子失去护种盘的支持力，种子将在重力的影响下运动至导种盘的导种槽内，完成递种过程。在携种阶段中，种子随着导种盘的转动运动至投种口，完成携种过程。在投种阶段中，外壳内壁面对种子的支持力消失，种子在重力的作用下从投种口排出，完成整个排种过程。
2. 关键部件设计
2.1 种子的形状特征及尺寸
种子的形状特征及尺寸是排种器关键零部件结构设计的重要依据。本文选用‘郑单958’玉米种子为研究对象，根据玉米种子的形状特征将宽厚比大于1.5的定义为扁平种子，小于或等于1.5的定义为圆粒种子。然后根据玉米种子尺寸大小分为4类：大圆粒、小圆粒、大扁平、小扁平，其中l0和l1分别为圆粒种子和扁平种子的长度，w0和w1分别为圆粒种子和扁平种子的宽度，t0和t1分别为圆粒种子和扁平种子的厚度。选取
1000 粒玉米种子，使用游标卡尺测量它们的三轴尺寸，测量结果以及2种种子的占比如表1所示。表 1 玉米种子三轴尺寸和分布占比Table 1. Three-axis size and distribution proportion of corn seeds种子类型
Seed type大小
Size长度
Length宽度
Width厚度
Thickness占比/%
Proportion最小值
Min最大值
Max平均值
Avg最小值
Min最大值
Max平均值
Avg最小值
Min最大值
Max平均值
Avg圆粒
Round大 10.57 11.68 11.05 8.55 9.71 9.07 7.51 8.49 7.73 10.6 小 9.52 10.49 9.92 7.31 8.47 7.89 6.16 7.44 6.87 11.6 平均 9.52 11.68 10.34 7.31 9.71 8.77 6.16 8.49 7.06 22.2 扁平
Flat大 11.53 12.82 12.13 8.82 9.98 9.42 5.06 5.84 5.41 48.5 小 10.46 11.42 10.99 7.62 8.79 8.27 4.17 5.00 4.55 29.3 平均 10.42 12.82 11.68 7.62 9.98 9.15 4.17 5.84 5.06 77.8 2.2 排种盘结构设计
排种盘直径决定型孔插件的数量及分布，对整机性能至关重要。《农业机械设计手册》[23]中，机械式排种器的种盘直径通常取140~260 mm，不宜过大或过小。排种盘直径过小时，排种盘可分布型孔插件数量减少，排种器高速作业时转速提升，不利于排种均匀、稳定；排种盘直径过大时，型孔插件数目可增多，有利于降低排种器的工作转速，提高排种稳定性，但直径过大也会相应提升制造成本、装配难度及作业功耗。故本文综合考虑后选取排种盘外径d0为220 mm，内径d1为180 mm。
为实现型孔插件的径向滑动，排种盘上设计型孔插件的滑动槽，型孔插件在压缩弹簧的弹力作用下在滑动槽内可以径向滑动。同时排种盘上设有搅种齿用于搅动种室内的种群，在增加种群流动性的同时改变种群中种子的姿态，提高充种效率(图2)。
![]()
图 2 排种盘装配示意图d0：排种盘外径，d1：排种盘内径。Figure 2. Assembly schematic diagram of seed metering plated0: Outer diameter of seed metering plate, d1: Inner diameter of seed metering plate.2.3 限姿板结构设计
限姿板是配合排种盘与型孔插件形成限姿通道，在种子进入充种型孔前预先调控种子姿态，使种子以理想的姿态到达护种盘与型孔插件共同组成的充种型孔前端的预充种位置。限姿板确保排种器在高速作业时种子能够快速地充入充种型孔，从而保证排种器在高速作业时的充种速度和充种效果的稳定，提升排种性能。限姿板的装配展开图与剖视图如图3所示。
![]()
图 3 限姿板结构示意图u：限姿通道竖直距离，v：限姿通道倾斜距离，e：限姿通道夹角。Figure 3. Structural diagram of posture limiting plateu: Vertical distance of posture limiting channel, v: Inclination distance of posture limiting channel, e: Angle of posture limiting channel.由于种群中扁平种子的占比较大，所以选取扁平种子为主要研究对象。根据测量的玉米种子的三轴尺寸特征，可以发现扁平种子的长度均大于10 mm，宽度和厚度均小于10 mm。因此设计限姿通道竖直距离u为10 mm，限姿通道倾斜距离v为10 mm，为使种子顺利滑至预充种位置，限姿通道夹角e应大于种子的摩擦角(28°)加90°，即118°。综合考虑，设计夹角e为135°，使限姿板与排种盘和护种盘配合来限制调控玉米种子的充前姿态。
2.3.1 分级限姿引导过程分析
根据种子的分级限姿引导过程，把种子受限区域根据限姿板、排种盘以及护种盘的结构特征分为4个区域(图4)。
![]()
图 4 玉米种子分级限姿引导过程分析Figure 4. Analysis of staged posture control and guiding process for corn seeds首先是一级限姿引导，种子从Π形进种管进入种室后到达区域①，在区域①的种子姿态大致分为6种，分别为：正面平躺、侧面平躺、侧面竖立、正面竖立、底面竖立和底面平躺，其中正面平躺和侧面平躺2个姿态的扁平种子由于长度方向上垂直于排种盘，而区域②限姿通道竖直距离u为10 mm，所以这2种姿态的种子无法进入到区域②之中，因此进入区域②的种子姿态大致分为4种，分别是：侧面竖立、正面竖立、底面竖立和底面平躺。其次是二级限姿引导，区域③限姿通道倾斜距离v为10 mm，区域②与区域③的限姿通道夹角e为135°，所以区域②的正面竖立姿态的玉米种子无法进入到区域③之中，其余3个姿态则沿限姿板斜面滑跃进入区域③，所以区域③的种子姿态大致分为4种，分别是侧面竖立、侧面斜躺、底面竖立和底面斜躺。最后是三级限姿引导，区域③的种子在重力、离心力和种群间作用力的综合作用下，到达区域④的预充种位置，区域④的种子姿态大致分为4种，分别是侧面竖立、侧面平躺、底面竖立和底面平躺，从而实现了调整玉米种子进入充种型孔前姿态的目的，使种子更快更稳定地进入充种型孔。
2.3.2 种子运动轨迹力学分析
根据对种子分级限姿引导过程的分析，将种子的运动轨迹分为轨迹一与轨迹二2种情况(图5a)。并对限姿通道的4个区域内种子的运动过程进行受力分析(图5b)。
![]()
图 5 种子限姿过程中的运动轨迹与受力分析Fa：限姿板对种子的支持力，fa：种子与限姿板之间的摩擦力，Fb：排种盘对种子的支持力，fb：种子与排种盘之间的摩擦力，Fc：种子之间的接触力，Fe：型孔插件与种子之间的接触力，Fd1：护种盘底部对种子的支持力，Fd2：护种盘后侧对种子的支持力，Fr：种子受到的离心力，Fa’：限姿板对种子的支持力在y轴上的分力，G：种子所受重力。Figure 5. Motion trajectory and force analysis in the process of seed posture limitingFa: Seed support force by posture limiting plate, fa: Friction between seed and posture limiting plate, Fb: Seed support force by seed metering plate, fb: Friction between seed and seed metering plate, Fc: Contact force between seeds, Fe: Contact force between orifice plug and seed, Fd1: Seed support force by the bottom of seed protection plate, Fd2: Seed support force by the posterior side of seed protection plate, Fr: Centrifugal force on seed, Fa’: y-axis component of seed support force by posture limiting plate, G: Gravity on seed.在区域①中，种子受力为限姿板对种子的支持力Fa、种子与排种盘之间的摩擦力fb、种子与限姿板之间的摩擦力fa、种子之间的接触力Fc，向种子质心取距可得：
$$ \sum {M = } 0.5{l_1}({F_{\mathrm{c}}} - {f_{\mathrm{b}}}) + \frac{{\sqrt 2 }}{2}\Delta d({F_{\mathrm{a}}} - {f_{\mathrm{a}}}){\text{，}} $$ (1) 式中，M为种子受力×力臂所产生的使种子旋转的力矩，Δd为种子与限姿板接触点的质心距，l1为扁平种子的平均长度。
种子在排种盘径向(y轴向)的合力
${F_y} $ 与轴向(x轴向)的合力${F_x} $ 满足：$$ \left\{ \begin{gathered} \sum {{F_y} = G + {F_{\mathrm{r}}} + {F_{\mathrm{c}}} + \frac{{\sqrt 2 }}{2}({f_{\mathrm{a}}} - {F_{\mathrm{a}}}) - {f_{\mathrm{b}}} \geqslant 0} \\ \sum {{F_x} = {F_{\mathrm{b}}} - \frac{{\sqrt 2 }}{2}({F_{\mathrm{a}}} + {f_{\mathrm{a}}}) = 0} \\ \end{gathered} \right. {\text{，}}$$ (2) 其中，
$$ \left\{ \begin{gathered} G = mg \\ {f_{\mathrm{a}}} = \mu {F_{\mathrm{a}}} \\ {f_{\mathrm{b}}} = \mu {F_{\mathrm{b}}} \\ {F_{\mathrm{r}}} = m{\omega ^2}r \\ \end{gathered} \right.{\text{，}} $$ (3) 式中，m为玉米种子的质量，g为重力加速度，取9.8 m/s2，μ为种子与排种器间滑动摩擦因数，ω为排种器角速度，G为种子所受重力，r为种子受到的离心力半径，Fr为种子受到的离心力，Fb为排种盘对种子的支持力。
当种子即将调整姿态时，质心距
$\Delta d $ 约为种子长度的0.5倍，因此取Δd为0.5l1，代入公式(1)并联立公式(2)、(3)，得到种子质心处的合力距大于0，因此种子将绕质心进行转动。完成姿态调控后，种子由区域①进入区域②。在区域②中，当长轴垂直于排种盘的种子完成调姿后，种子可以分为2种运动轨迹进行分析。
第2种是如图5轨迹一所示，种子姿态满足从区域②下落至区域③，此时种子在区域②中的Fy与Fx满足：
$$ \left\{ \begin{gathered} \sum {{F_y} = G + {F_r} + {F_c} \geqslant 0} \\ \sum {{F_x} = 0} \\ \end{gathered} \right.{\text{，}} $$ (4) 由式(4)可知，在合外力的作用下，此时种子将沿排种盘径向方向由区域②下落至区域③。
第2种是如图5轨迹二所示，种子姿态保持正面竖立，此时由于限姿板的角度e，部分大扁平种子将无法从区域②下落至区域③，此时种子受力为型孔插件与种子之间的接触力Fe、Fc，向种子质心取距可得
$$ \sum {M = } 0.5{l_1}{F_{\mathrm{c}}} + \Delta {d_2}{F_{\mathrm{e}}}\cos \varphi {\text{，}}$$ (5) 式中，Δd2 为种子与型孔插件接触点的质心距，φ 为型孔插件侧面角度。
种子的
${F_y} $ 与${F_x} $ 满足$$ \left\{ \begin{gathered} \sum {{F_y} = G + {F_{\mathrm{r}}} + {F_{\mathrm{e}}}\sin \varphi + {F_{\mathrm{a}}}' = 0} \\ \sum {{F_x} = 0.5{l_1}{F_{\mathrm{c}}} + \Delta {d_2}{F_{\mathrm{e}}}\cos \varphi \geqslant 0} \\ \end{gathered} \right. {\text{，}}$$ (6) 其中，
$$ {F_{\mathrm{a}}}' = \frac{{\sqrt 2 }}{2}{F_{\mathrm{a}}}{\text{，}} $$ (7) 式中，
${F_{\mathrm{a}}}' $ 为限姿板对种子的支持力在y轴上的分力。当种子即将调整姿态时，质心距
$\Delta d_2$ 约为种子长度的0.3倍，因此取Δd2为0.3l1，代入公式(5)并联立公式(3)、(6)、(7)，得到种子质心处的合力距大于0，因此种子将绕质心进行转动。完成姿态调控后，种子以底面平躺姿态平躺在型孔插件上方。此时种子的
${F_y} $ 与${F_x} $ 满足$$ \left\{ \begin{gathered} \sum {{F_y} = G + {F_{\mathrm{r}}} + {F_{\mathrm{c}}} - {F_{\mathrm{e}}} = 0} \\ \sum {{F_x} = 0} \\ \end{gathered} \right. {\text{，}}$$ (8) 由式(8)可知，此时种子在排种盘径向截面中的合外力为0。而随着排种盘的转动，种子失去型孔插件的支持力Fe后，由区域②进入区域③。
在区域③中，部分保持侧面竖立姿态的大扁平种子在刚到达时，种子受力为Fb、Fa、fa、Fc，向种子质心取距可得：
$$ \sum {M = } \frac{{\sqrt 2 }}{4}{l_1}({F_{\mathrm{c}}} + {F_{\mathrm{a}}} - {f_{\mathrm{a}}}) + \Delta {d_1}{F_{\mathrm{b}}}{\text{，}} $$ (9) 式中，Δd1为种子与排种盘接触点的质心距。
种子的
${F_y} $ 与${F_x} $ 满足：$$ \left\{ \begin{gathered} \sum {{F_y} = G + {F_{\mathrm{r}}} + \frac{{\sqrt 2 }}{2}({F_{\mathrm{c}}} - {f_{\mathrm{a}}} - {F_{\mathrm{a}}}) - {f_{\mathrm{b}}} = 0} \\ \sum {{F_x} = \frac{{\sqrt 2 }}{2}({F_{\mathrm{a}}} + {F_{\mathrm{c}}} - {f_{\mathrm{a}}}) - {F_{\mathrm{b}}} \geqslant 0} \\ \end{gathered} \right. {\text{。}}$$ (10) 当种子即将调整姿态时，质心距
$\Delta d_1 $ 约为种子长度的0.5倍，因此取Δd1为0.5l1，代入公式(9)并联立公式(3)、(10)，得到种子质心处的合力距大于0，因此种子将绕质心进行转动完成姿态调控。完成姿态调控后，区域③的种子在与限姿板接触面法线方向(y轴向)的合力
${F_y}' $ 和接触面方向(x轴向)的合力${F_x}' $ 满足：$$ \left\{ \begin{gathered} \sum {{F_y}' = \frac{{\sqrt 2 }}{2}(G + {F_{\mathrm{r}}}) - {F_{\mathrm{a}}} = 0} \\ \sum {{F_x}' = \frac{{\sqrt 2 }}{2}(G + {F_{\mathrm{r}}}) + {F_{\mathrm{c}}} - {f_{\mathrm{a}}} \geqslant 0} \\ \end{gathered} \right. {\text{，}}$$ (11) 由式(11)可知，完成姿态调控后的侧面竖立的大扁平种子与其余3种姿态的种子在合外力的作用下，沿种子在与限姿板接触面方向由区域③滑至区域④。
在区域④中，种子到达预充种位置后，种子在与护种盘接触面法线方向(y轴向)的合力
$F_y'' $ 和接触面方向(x轴向)的合力$F_x'' $ 满足：$$ \left\{ \begin{gathered} \sum {{F_y}'' = (G + {F_{\mathrm{r}}})\sin \alpha - {F_{d2}} = 0} \\ \sum {{F_x}'' = (G + {F_{\mathrm{r}}})\cos \alpha + {F_{\mathrm{c}}} - {F_{d1}} = 0} \\ \end{gathered} \right.{\text{，}} $$ (12) 式中，α为护种盘底部角度，Fd1为护种盘底部对种子的支持力，Fd2为护种盘后侧对种子的支持力。
由式(12)可知，种子到达预充种位置后，在排种盘径向截面中种子的合外力为0。随着排种盘的转动，种子准备进入充种型孔，完成充种过程。
2.4 充种型孔结构设计
2.4.1 充种型孔填充姿态分析
为了达到单粒精密播种的要求，排种器在工作时应能将有且仅有一粒种子填入充种型孔内。如果充种型孔足够大，种子在充种型孔内的姿势可以大致分为6种：正面竖立、侧面竖立、底面竖立、正面平躺、侧面平躺、底面平躺。此时这个充种型孔满足种子的各种填充姿态，但是充种型孔就会很大，充种时会将多粒种子侧充进入充种型孔中，这将严重降低排种器的性能，大大增加排种器作业时的重播率。因此，在设计充种型孔时，必须选择理想的充填姿态，以达到只充填一粒种子的目的。
2.4.2 侧充过程分析
充种型孔的结构与尺寸大小应依据充入型孔内的种子姿态来确定。当种子呈正面竖立、侧面竖立、正面平躺或侧面平躺时，它们的长度方向均垂直于型孔插件与护种盘构成的充种型孔的下底面和左侧面，这将导致充种型孔过大，无法确保在清种后孔内仅剩1粒种子。因此，选择保留底面竖立和底面平躺这2种姿态，并通过调整护种盘与型孔插件的结构参数来限制或调整其他4种姿态。
种子经过限姿板限制姿态之后到达区域④的预充种位置，预充种位置的种子会随着排种盘的转动进行姿态调整，然后侧充进入充种型孔内直至达到稳定状态。如图6所示，底面平躺姿态的种子随着排种盘的转动，直接充入充种型孔。侧面平躺姿态的种子随着排种盘的转动，当种子与型孔插件发生接触时，型孔插件会对种子施加一个推力F2。由于种子与型孔的接触点在种子的边缘区域，支持力的作用点与种子质心之间存在一定距离，此时产生力矩M2。在M2的作用下，种子姿态被调整为底面平躺姿态，使其能够顺利进入充种型孔，随后通过排种盘的转动被充入充种型孔。对于底面竖立和侧面竖立姿态的种子，它们受到种群间其他种子的一个作用力Fc和护种盘上斜面的支持力F3，此时会产生力矩M3，在M3的作用下，种子姿态分别被调整为底面斜躺和侧面斜躺2个姿态，最终随着排种盘的转动侧充进入充种型孔并达到稳定姿态。
![]()
图 6 种子由区域④侧充进入充种型孔的充种过程及稳定姿态F2：型孔插件对种子的推力，Fc：种子之间的接触力，F3：护种盘上斜面对种子的支持力，M2：F2产生的力矩，M3：Fc与F4产生的力矩。Figure 6. Seed filling process and stable posture from region ④ to the filling orifice through side-fillingF2: Thrust of the orifice plug on the seed, Fc: Contact force between seeds, F3: Supporting force of the inclined plane on seed protection plate to the seed, M2: Torque generated by F2, M3: Torque generated by Fc and F4.2.4.3 型孔结构参数设计
由于型孔插件与护种盘共同组成的充种型孔直接与种子接触，因此，这2个零件结构的形状和参数对排种器工作性能有着重要的影响。
型孔插件结构形状如图7c所示，护种盘结构形状如图7b所示，两者结构参数中与充种型孔相关的是型孔插件顶部宽度W、型孔插件底部角度α0、型孔插件底部长度L0、型孔插件前侧高度H0、型孔插件顶部角度β0、和型孔插件侧面角度φ，护种盘后侧高度H、护种盘底部长度L、护种盘底部角度α、护种盘顶部角度β。
![]()
图 7 护种盘与型孔插件的结构参数设计α：护种盘底部角度，β：护种盘顶部角度，c：护种板顶部厚度，H：护种盘后侧高度，L：护种盘底部长度，α0：型孔插件底部角度，β0：型孔插件顶部角度，L0：型孔插件底部长度，H0：型孔插件前侧高度，b：型孔插件后侧厚度，e0：型孔插件后侧角度，c0：型孔插件后侧高度，W：型孔插件顶部宽度，φ：型孔插件侧面角度。Figure 7. kinds of plate type and structure parameters of hole plug-in protection designα: Bottom angle of seed protection plate, β: Top angle of seed protection plate, c: Top thickness of seed protection plate, H: Rear height of seed protection plate, L: Bottom length of seed protection plate, α0: Bottom angle of orifice plug, β0: Top angle of orifice plug, L0: Bottom length of orifice plug, H0: Front height of orifice plug, b: Rear thickness of orifice plug, e0: Rear angle of orifice plug, c0: Rear height of orifice plug, W: Top width of orifice plug, φ: Side angle of orifice plug.其中φ的设计主要是为了减小初始清种角，同时确保不影响充种效果。综合考虑，选取φ=15°。为了确保充种型孔的充种效果，型孔插件和护种盘与充种型孔相关的结构参数之间需要满足下式：
$$ \left\{ \begin{gathered} \alpha = {\alpha _0} \\ \beta = {\beta _0} \\ L = {L_0} + 0.5 \\ H = {H_0} + 0.5 \\ \end{gathered} \right.{\text{，}} $$ (13) 式中，0.5为型孔插件与护种盘的配合间隙宽度。
其余结构参数分别为型孔插件后侧厚度b、型孔插件后侧角度e0、型孔插件后侧高度c0、护种板顶部厚度c。结合限姿板与排种盘的尺寸特征，综合考虑，取b为9.5 mm，e0为135°，c0为2 mm。为了保证后期台架试验材料制作的护种盘质量，综合考虑，取c为3 mm。
为了确保充种型孔不仅适用于扁平种子，还能有效充入圆粒种子，通过对比扁平种子与圆粒种子的形状特征得到护种盘顶部角度β需要满足下式：
$$ \frac{{{w_0} - {t_1}}}{{0.5{w_0}}} < \tan (180 - \beta ) $$ (14) 式中，w0为圆粒种子的平均宽度，t1为扁平种子的平均厚度。
根据公式(14)求得β<49.77°。β不仅影响充种阶段的充种效果，而且对清种阶段有影响，因此，综合考虑，取β=45°。
充种型孔的参数取值应依据充入型孔内的种子姿态决定，由图6的侧充过程分析结合扁平种子与圆粒种子的形状特征，对型孔插件与护种盘的结构参数进行设计。
如图8的扁平种子和圆粒种子与充种型孔的结构示意图所示，充种型孔同时满足充入扁平种子与圆粒种子条件时，参数H、L、α应满足下式：
![]()
图 8 扁平种子(a)和圆粒种子(b)与充种型孔的结构示意图α：护种盘底部角度，β：护种盘顶部角度，H：护种盘后侧高度，L：护种盘底部长度，w0：圆粒种子宽度，t1：扁平种子厚度，w1：扁平种子宽度。Figure 8. Structural diagram of seed-filling orifices with flat seed (a) and round seed(b)α: Bottom angle of seed protection plate, β: Top angle of seed protection plate, H: Rear height of seed protection plate, L: Bottom length of seed protection plate, w0: Round seed width, t1: Flat seed thickness, w1: Flat seed width.$$ \left\{ \begin{gathered} {t_{1\min }} \leqslant H \leqslant {t_{1\max }} \\ {w_{1\max }} \leqslant L \leqslant {l_{1\min }} \\ {w_0} \leqslant \frac{{0.5L}}{{\tan \beta }} + H - \frac{{0.5L}}{{\tan \alpha }} \\ \tan \alpha \leqslant {90^{\text{o}}} \\ \end{gathered} \right. {\text{，}}$$ (15) 式中，w1max 为扁平种子宽度最大值，t1min 为扁平种子厚度最小值，t1max 为扁平种子厚度最大值，l1min 为扁平种子长度最小值。
通过公式(15)得到L的范围为9.98~10.42 mm，H的范围为4.17~5.86 mm，α的范围为76.02°~90.00°。综合考虑，取L=10 mm，而由于H和α取值范围过大，为了提高充种型孔对种子的适应性，H和α的取值依据具体试验优化结果。
根据侧充过程分析结果得到扁平种子侧面斜躺与底面平躺充种型孔的结构示意图(图9)。为了保证排种盘转动时充种型孔能够顺利充种与清种阶段充种型孔内的种子不掉落，则型孔插件顶部宽度
$W $ 的优选范围应靠近种子的质心，即应满足下式：![]()
图 9 扁平种子侧面斜躺(a)和底面斜躺(b)与充种型孔的结构示意图W：型孔插件顶部宽度，φ：型孔插件侧面角度，w1：扁平种子宽度，l1：扁平种子长度。Figure 9. Structural diagram of seed filling orifices with flat seeds in lateral oblique(a) and basal oblique(b) posturesW: Top width of orifice plug, φ: Side angle of orifice plug, w1: Flat seed width, l1: Flat seed length.$$ 0.5{w_{1\min }} \leqslant W \leqslant 0.5{l_{1\max }}{\text{，}} $$ (16) 式中，w1min 为扁平种子宽度最小值，l1max 为扁平种子长度最大值。
通过公式(16)得到W的范围为3.81~6.41 mm。由于W取值范围过大，为了提高充种型孔对种子的适应性，提高清种效果，W的取值依据具体试验优化结果。
2.4.4 型孔插件数量分布
排种器在满足充种、清种、递种、携种和投种5个工作阶段的需求下，应尽量增多型孔插件的数量，以便降低排种盘工作转速并提高投种频率，进而提高排种器的作业速度。型孔插件数量分布相关参数如图10所示。
![]()
图 10 型孔插件数量分布相关参数r1：导种槽递种口位置半径，q：递种口位置宽度，k：相邻两型孔插件的距离，i：导种槽两侧厚度。Figure 10. Parameters related to the quantity distribution of the orifice plugsr1: Radius at seed delivery port of seed guiding channel, q: Width of seed delivery port, k: Distance between two adjacent orifices, i: Thickness at both sides of seed guiding channel.为确保种子能够及时到达预充种位置，相邻两型孔插件的距离
$k $ 应大于1.5倍的种子最大长度，同时为了保证种子递种时能够顺利进入导种盘的导种槽内，导种槽的递种口位置宽度$q $ 应大于种子最大长度，即型孔插件数量$z $ 需满足下式$$ \left\{ \begin{array}{l} \dfrac{{{\text{π}} {d_1}}}{{\textit{z}}} - {d_1}\arcsin \left( {\dfrac{{{W_{\max }}}}{{{d_1}}}} \right) > {d_1}\arcsin \left( {\dfrac{k}{{{d_1}}}} \right) \\ \dfrac{{2{\text{π}} {r_1}}}{{2{\textit{z}}}} - 2{r_1}\arcsin \left( {\dfrac{i}{{2{r_1}}}} \right) > 2{r_1}\arcsin \left( {\dfrac{q}{{2{r_1}}}} \right) \end{array} \right. $$ (17) 式中，r1 为导种槽递种口位置半径，Wmax 为型孔插件顶部宽度最大值，i 为导种槽两侧厚度。
取i =3 mm，r1=95 mm，根据式(17)求得z<18.86，综合考虑，取z=18。
2.5 护种盘径向突变结构设计
根据侧充过程分析，当充种状态稳定之后，充种型孔内种子姿态大致分为3种，而底面平躺姿态下充种型孔内的剩余空间可能会有部分较小的种子卡入，导致种子难以利用自重快速完成清种，所以需要考虑底面平躺姿态的清种问题。
因此，设计一个径向突变结构辅助种子自重清种，从而加快清种速度，以达到减小重播指数的目的。如图11所示，型孔插件经过径向突变结构时，在压缩弹簧的弹力作用下实现变径。保留种子受到型孔插件对种子的摩擦力f1，减小了对待清种子的压力。待清种子受到型孔插件对种子的摩擦力f2，同时受到重力G，此时会产生力矩M4，最后种子在力矩M4的辅助下沿护种盘斜面下滑，达到快速清种的目的。
![]()
图 11 径向突变辅助清种分析f1：型孔插件对保留种子的摩擦力，f2：型孔插件对待清种子的摩擦力，M4：f2与G产生的力矩，G：种子所受重力。Figure 11. Analysis of seed cleaning assisted by radial mutationf1: Frictional force of orifice plug on retained seed, f2: Frictional force of orifice plug on the seed to be cleared, M4: Torque generated by f2 and G, G: Gravity on the seed.为了保证甩出的型孔插件能够在及时复位，设计1个护种盘椭圆结构将型孔插件推回原位，如图12a所示，其中椭圆短轴a为护种盘圆周半径R，长轴b为护种盘圆周半径R加上径向突变宽度p。护种盘径向突变轮廓曲线如图12b所示。
![]()
图 12 护种盘径向突变轮廓曲线设计a：椭圆短轴，b：椭圆长轴，R：护种盘圆周半径，p：径向突变宽度，θ：径向突变起始角度。Figure 12. Design of radial abrupt contour curve of seed protection platea: Short axis of ellipse, b: Long axis of ellipse, R: Radius of seed protection plate, p: Radial mutation width, θ: Radial mutation starting angle.为了更好地辅助种子自重清种，避免影响递种过程，综合考虑，径向突变起始角度θ取65°，p取1.5 mm。
3. 充种型孔参数优化仿真试验
充种型孔的结构参数是影响排种质量的重要因素。为了确定充种型孔的最佳结构参数组合，设计了基于DEM-MBD耦合参数优化仿真试验。
3.1 仿真模型的建立
3.1.1 物理参数设置
为了保证仿真模拟结果的准确性，需要获取排种器和玉米种子的物理参数与接触参数。排种器与种子接触的关键部件的材料均为树脂。根据前期试验得到种子的泊松比、剪切模量和密度分别为0.4、1.37×108 Pa和
1197 kg·m−3。树脂的泊松比、剪切模量和密度分别为0.5、1.77×108 Pa和1180 kg·m−3。种子间的碰撞恢复系数、静摩擦系数和动摩擦系数分别为0.18、0.03和0.01。种子与树脂间的碰撞恢复系数、静摩擦系数和动摩擦系数分别为0.71、0.46和0.09。3.1.2 玉米种子模型建立
本研究中，DEM方法是基于EDEM 2021软件来完成的。为了使模拟更真实，选取‘郑单958’玉米杂交种作为建模对象，在EDEM中采用颗粒填充法将种子分别填充建模，如图13所示。
![]()
图 13 基于‘郑单958’玉米不同类型种子的填充模型Figure 13. Filling models based on ‘Zhengdan 958’ corn seeds with different types3.1.3 仿真模型建立
本研究中，DEM-MBD方法是基于EDEM 2021与RecurDyn 2023软件来完成的。首先将SolidWorks中将装配好的排种器三维模型转化为x-t格式，导入至RecurDyn中，利用RecurDyn 对机构添加joint、force、contact等仿真参数。参数设置完成后，利用RecurDyn 中External SPI功能模块下的EDEM接口模块将各部件以wall.文件形式导出。在EDEM中Geometries模块下的Import Geometry from RecurDyn 功能将RecurDyn中生成的wall.文件导入。然后在EDEM中建立颗粒工厂生成颗粒，颗粒工厂在1 s内生成全部所需玉米籽粒。仿真完成后，在菜单栏File的Export功能下的Simulation Deck导出生成一个input文件，勾选Set Simulation Time并设置成为0 s。然后打开input文件，设置Fixed Time Step为1×10−6 s，最后打开RecurDyn的Analysis功能模块下的Dyn/Kin，设置仿真时间为20 s，step为1 000进行耦合仿真。
为了能对种子排种时间间隔进行检测，在落种口处设置收集器，并通过EDEM后处理功能添加监测器，记录种子与前一颗种子的排出时间间隔，如图14所示。参照《单粒(精密)播种机试验方法》(GB/T 6973-2005)[24]，排种时间间隔在0.5~1.5倍理论时间间隔的为合格，排种时间间隔大于1.5倍理论时间间隔的为漏播，排种时间间隔小于0.5倍理论时间间隔的为漏播。根据计算，当播种株距为25 cm、作业速度为12 km/h时，理论排种时间间隔为0.075 s，则当时间间隔小于0.039 s时为重播，时间间隔在0.039~0.113 s之间的为合格，时间间隔大于0.113 s的为漏播。
![]()
图 14 排种器EDEM后处理1：种子填充模型，2：排种器模型，3：收集器，4：监测器，5：收集器内种子质量。Figure 14. EDEM post-processing for seed metering device1: Seed filling model, 2: Seed metering device model, 3: Collector, 4: Monitor, 5: Seed mass in the collector.3.2 参数优化仿真试验
3.2.1 试验因素与试验指标
根据前面的理论分析选定型孔插件顶部宽度W、护种盘后侧高度H、护种盘底部角度α为影响排种器工作性能的主要因素，以合格指数Y1、重播指数Y2和漏播指数Y3为试验指标，以251粒种子为1组，记录其排种时间间隔，共250个时间间隔，每组试验重复3次，取平均值；各个指标的计算公式为：
$$ \left\{ {\begin{array}{*{20}{c}} {\mathop Y\nolimits_1 = \dfrac{{{{{n}}_1}}}{N} \times 100{\text{%}} } \\ {\mathop Y\nolimits_2 = \dfrac{{{{{n}}_2}}}{N} \times 100{\text{%}} } \\ {\mathop Y\nolimits_3 = \dfrac{{{{{n}}_3}}}{N} \times 100{\text{%}} } \end{array}} \right. {\text{，}}$$ (18) 式中，n1为排种时间间隔在0.5~1.5倍理论时间间隔的排种数，n2为排种时间间隔大于1.5倍理论时间间隔的排种数，n3为排种时间间隔小于0.5倍理论时间间隔的排种数，N为理论排种数。
3.2.2 试验方案
根据充种型孔结构参数对排种影响的理论分析，选取试验因素的取值范围为：根据前面的理论分析选定W为3.81~6.41 mm、H为4.17~5.86 mm、α为76.02°~90.00°。针对上述关键因素开展三因素五水平二次正交旋转回归组合试验，试验编码表如表2 所示。
表 2 排种试验因素和水平Table 2. Factors and levels of seed metering experiment水平
Level因素1) Factor X1/mm X2/mm X3/(°) 1.682 6.41 5.86 90.00 1 5.88 5.52 87.17 0 5.11 5.02 83.01 −1 4.34 4.52 78.85 −1.682 3.81 4.17 76.02 1) X1：型孔插件顶部宽度，X2：护种盘后侧高度，X3：护种盘底部角度。
1) X1: Top width of orifice plug, X2: Rear height of seed protection plate, X3: Bottom angle of seed protection plate.3.2.3 试验结果与分析
基于EDEM和RecurDyn耦合仿真的二次正交旋转回归组合试验的试验方案与结果如表3 所示，利用Design-expert 13 软件对试验结果进行回归分析，试验因素对试验指标的显著性检测结果见表4。
表 3 排种试验方案及试验结果Table 3. Design and result of seed metering experiment试验编号
Test number因素1) Factor 评价指标 Evaluation index X1/mm X2/mm X3/(°) 合格指数/%
Qualified index重播指数/%
Replay index漏播指数/%
Leakage index1 −1 −1 −1 91.47 2.28 6.25 2 1 −1 −1 86.46 8.71 4.83 3 −1 1 −1 92.16 4.14 3.70 4 1 1 −1 85.48 12.06 2.46 5 −1 −1 1 88.73 8.45 2.82 6 1 −1 1 82.74 15.11 2.15 7 −1 1 1 83.27 15.41 1.32 8 1 1 1 78.56 20.21 1.23 9 −1.682 0 0 91.95 3.85 4.20 10 1.682 0 0 82.43 16.03 1.54 11 0 −1.682 0 86.26 8.48 5.26 12 0 1.682 0 82.58 15.63 1.79 13 0 0 −1.682 91.21 3.83 4.96 14 0 0 1.682 82.12 16.23 1.65 15 0 0 0 91.92 6.25 1.83 16 0 0 0 90.61 6.77 2.62 17 0 0 0 89.54 7.96 2.50 18 0 0 0 90.24 7.58 2.18 19 0 0 0 90.16 7.91 1.93 20 0 0 0 91.39 6.93 1.68 21 0 0 0 90.29 7.60 2.11 22 0 0 0 91.57 6.57 1.86 23 0 0 0 90.18 8.01 1.81 1) X1：型孔插件顶部宽度，X2：护种盘后侧高度，X3：护种盘底部角度。
1) X1: Top width of orifice plug, X2: Rear height of seed protection plate, X3: Bottom angle of seed protection plate.根据表4 可知，合格指数、重播指数以及漏播指数的回归模型均为极显著(P<0.01)，失拟项均不显著(P>0.05)，表明回归模型的拟合程度较高，且模型的R2 均大于0.92，表明模型可以解释评价指标超过92% 的变化。对于合格指数回归模型，X1、X2、X3、X2X3、X12、X22 和X32 均为极显著，其余项均不显著；对于重播指数回归模型，X1、X2、X3、X2X3、X12、X22 和X32 均为极显著，其余项均为不显著；对于漏播指数回归模型，X1、X2、X3、X22和X32均为极显著，X2X3、X12为显著，其余项均不显著。
表 4 模型方差分析Table 4. Variance analysis of model方差来源
Variation source合格指数 Qualified index 重播指数 Replay index 漏播指数 Leakage index 平方和 SS Df F P1) 平方和 SS Df F P1) 平方和 SS Df F P1) 模型 Model 363.36 9 81.25 < 0.0001 **502.41 9 125.11 < 0.0001 **42.65 9 37.14 < 0.0001 **X1 107.98 1 217.29 < 0.0001 **156.93 1 351.7 < 0.0001 **4.56 1 35.76 < 0.0001 **X2 19.03 1 38.29 < 0.0001 **62.84 1 140.83 < 0.0001 **12.71 1 99.64 < 0.0001 **X3 103.29 1 207.85 < 0.0001 **204.48 1 458.26 < 0.0001 **17.11 1 134.12 < 0.0001 **X1X2 0.02 1 0.04 0.8479 0.02 1 0.04 0.8478 0.07 1 0.57 0.4653 X1X3 0.12 1 0.25 0.6278 1.04 1 2.34 0.1501 0.45 1 3.54 0.0826 X2X3 10.93 1 21.99 0.0004 **5.87 1 13.14 0.0031 **0.78 1 6.12 0.0279 *X12 21.54 1 43.35 < 0.0001 **13.47 1 30.18 0.0001 **0.94 1 7.4 0.0175 *X22 73.02 1 146.94 < 0.0001 **44.23 1 99.12 < 0.0001 **3.59 1 28.13 0.0001 **X33 28.96 1 58.27 < 0.0001 **14.41 1 32.3 < 0.0001 **2.51 1 19.68 0.0007 **残差 Residual 6.46 13 5.8 13 1.66 13 失拟 Lack of fit 1.46 5 0.47 0.7909 2.27 5 1.03 0.4618 0.82 5 1.56 0.2751 误差 Error 5 8 3.53 8 0.84 8 总和 Total 369.82 22 508.21 22 44.30 22 1) **表示影响极显著 (P<0.01)；*表示影响显著 (P<0.05)。
1) ** indicates extremely significant influence (P<0.01); * indicates significant influence (P<0.05).剔除回归方程中影响不显著因素的系数，建立评价指标与试验因素的二次回归方程为
$$ \left\{ \begin{gathered} {Y_1} = 90.65 - 2.81{X_1} - 1.18{X_2} - 2.75{X_3} -\\ \quad\quad 1.17{X_2}{X_3} - 1.16X_1^2 - 2.14X_2^2 - 1.35X_3^2 \\ {Y_2} = 7.29 + 3.39{X_1} + 2.15{X_2} + 3.87{X_3} +\\ \quad\quad 0.86{X_2}{X_3} + 0.92X_1^2 + 1.67X_2^2 + 0.95X_3^2 \\ {Y_3} = 2.06 - 0.58{X_1} - 0.96{X_2} - 1.12{X_3} +\\ \quad\quad 0.31{X_2}{X_3} + 0.24X_1^2 + 0.48X_2^2 + 0.40X_3^2 \\ \end{gathered} \right.{\text{。}} $$ (19) 为了探究试验因素对各指标的影响，分析仿真过程中结构参数的变化导致重播和漏播的原因，利用Design expert 13 软件绘制试验因素交互作用显著的响应曲面，如图15所示。
![]()
图 15 因素间交互作用对排种器性能影响的响应曲面Figure 15. Response surface of influence of interaction between factors on performance of seed metering device由表4和图15a可知，护种盘后侧高度与护种盘底部角度的交互项X2X3对合格指数的影响极显著。当型孔插件顶部宽度一定时，随着护种盘底部角度和护种盘后侧高度的增大，从而种子进入充种型孔的难度减小，有利于种子的充填，所以合格指数增大，但随着充种型孔的一直增大，则容易发生2粒种子同时进入充种型孔，所以合格指数减小，从而导致合格指数先增大后减小。
由表4和图15b可知，护种盘后侧高度与护种盘底部角度的交互项X2X3对重播指数的影响极显著。当型孔插件顶部宽度一定时，随着护种盘底部角度和护种盘后侧高度的增大，充种型孔逐渐变大，则容易发生2粒种子同时进入充种型孔，从而导致重播指数一直增大。
由表4和图15c可知，护种盘后侧高度与护种盘底部角度的交互项X2X3对漏播指数的影响显著。当型孔插件顶部宽度一定时，随着护种盘底部角度和护种盘后侧高度的减小，充种型孔逐渐减小，则充种时种子进入充种型孔难度增加，从而导致漏播指数一直增大。
为得到该排种器在高速作业条件下的充种型孔最优结构参数组合，以合格指数最大、重播与漏播指数最小为优化目标，对试验指标进行优化求解，目标函数与约束条件为：
$$ \left\{ \begin{gathered} \max {Y_1}\left( {{X_1},{X_2},{X_3}} \right) \\ \min {Y_2}\left( {{X_1},{X_2},{X_3}} \right) \\ \min {Y_3}\left( {{X_1},{X_2},{X_3}} \right) \\ {\mathrm{s.t.}}\left\{ \begin{gathered} 3.81\;{\mathrm{mm}} \leqslant W \leqslant 6.41\;{\mathrm{mm}} \\ 4.17\;{\mathrm{mm}} \leqslant H \leqslant 5.86\;{\mathrm{mm}} \\ {76.02^ \circ } \leqslant \alpha \leqslant {90.00^ \circ } \\ \end{gathered} \right. \\ \end{gathered} \right. {\text{，}}$$ (20) 求解得出型孔插件顶部宽度W为4.65 mm、护种盘后侧高度H为5.10 mm、护种盘底部角度α为81.86°时，排种质量最优，此时排种合格指数、重播指数和漏播指数分别为92.41%、4.88%和2.71%。为验证优化结果的准确性，在相同条件下进行3 组仿真试验，结果取平均值，得到排种器的合格指数、重播指数和漏播指数分别为92.18%、5.13% 和2.69%，接近预测值，说明二次正交旋转回归组合试验的优化结果准确。
3.2.4 种子受力破碎风险分析
在排种器工作过程中，随着作业速度的提高，种子受到的部件载荷大小也会随之发生变化。因此，为了得到排出的种子所受最大载荷对种子破碎的影响，在W为4.65 mm、H为5.10 mm、α为81.86°的条件下，作业速度分别为8、10、12、14、16 km/h时，分别统计4类种子在排出过程中所受的最大合力，每类种子各取50粒，取平均值并绘制点线图，如图16所示。
![]()
图 16 4类种子所受合力与作业速度的关系Figure 16. Relationship between the resultant force on four types of seeds and the operation speed由图16可以看出：随着作业速度的提高，4类种子排出过程中所受合力总体均呈现上升趋势，且种子所受合力的最大值出现在作业速度16 km/h下的小扁平种子，为6.27 N。籽粒抗压碎力与籽粒抗剪切力是影响破碎率的2个重要指标。经过查阅文献可知，玉米品种‘郑单958’的抗压碎力为105 N，抗剪切力为71.2 N[25]，均显著高于种子所受合力的最大值。这表明双列限姿导引侧充式排种器在当前试验条件下，种子的破碎风险较低。
4. 台架性能验证试验
4.1 试验材料与装置
为了验证仿真试验的可靠性，开展台架性能验证试验，同时验证排种器在不同作业速度下的排种性能，试验台架如图17 所示。
![]()
图 17 排种器试验台架1：操控台，2：MeterMax ultra 性能检测器，3：排种器，4：光电传感器，5：变压器，6：收集器，7：充种型孔装配，8：径向突变，9：护种盘，10：排种盘装配，11：排种盘，12：滑动槽。Figure 17. Test bench of seed metering device1: Console, 2: MeterMax ultra performance detector, 3: Seed metering device, 4: Photoelectric sensor, 5: Transformer, 6: Collector, 7: Seed filling orifice assembly, 8: Radial mutation, 9: Seed protection plate, 10: Seed metering plate assembly, 11: Seed metering plate, 12: Slide groove.试验选用未分级的‘郑单958’玉米杂交种作为研究对象，千粒质量351 g，含水率12.5%。试验地点为云南师范大学能源与环境科学学院D-02实验室。试验使用的排种器主要零部件通过光固化3D 打印技术加工而成，材料为光敏树脂，其余轴类零部件使用45号钢机加工；试验检测装置使用MeterMax ultra排种器性能检测仪(美国P&C Ag Solutions)，得到排种器的排种合格指数、重播指数和漏播指数。
4.2 作业速度适应性试验
为了体现本排种器在不同作业速度下的适应性，将作业速度分为3个速度范围，低速：8~10 km/h，中速：11~13 km/h，高速：14~16 km/h，在W为4.65 mm、H为5.10 mm、α为81.86°时，分别在3个作业范围下进行试验，每组重复5次取平均值，试验结果如表5所示。
表 5 作业速度适应性试验结果Table 5. Results of operation speed adaptability test速度
Speed作业速度/(km·h−1)
Working speed评价指标/%
Evaluation index合格指数
Qualified
index重播指数
Replay
index漏播指数
Leakage
index低
Low8 92.44 4.37 3.19 9 93.05 4.21 2.74 10 92.53 4.65 2.82 中
Moderate11 92.09 4.90 3.01 12 91.45 5.05 3.50 13 91.52 5.17 3.31 高
High14 91.04 5.52 3.44 15 90.11 6.26 3.63 16 90.57 6.01 3.42 由作业速度适应性试验结果可知，在作业速度12 km/h 时，排种合格指数为91.45%，与仿真试验的结果基本吻合，两者相对误差为0.79%，产生误差的主要原因有两点：一是仿真模拟的玉米籽粒模型与真实籽粒形状存在差异，二是排种器虚拟样机模型为刚体，与台架试验的排种器存在差异。在作业速度8~11 m/h时，排种器的合格指数均能达到92%以上，重播指数均低于5%，漏播指数均低于3.5%，在作业速度12~16 m/h时，排种器的合格指数均能达到90%以上，重播指数均低于7%，漏播指数均低于4%，满足国家标准的播种要求，因此本排种器对作业速度有较强的适应性。
4.3 排种器性能对比试验
基于上述仿真试验与台架试验的结果，双列限姿导引侧充式玉米高速精量排种器在W为4.65 mm、H为5.10 mm、α为81.86°的条件下，与勺轮式排种器分别开展作业速度为8、10、12、14、16 km/h的排种性能对比试验。在上文作业速度适应性台架试验方法的基础上，添加每组试验用种的破损率Y4统计结果，综合评价该排种器的高速作业性能。
由排种性能对比试验结果(表6)可知，相较传统的勺轮式排种器，当作业速度为8 km/h时，双列限姿导引侧充式排种器的作业效果与勺轮式排种器差异较小，排种合格指数提升2.10个百分点，种子破损率降低0.03个百分点，但随作业速度提升差异逐渐增大，在作业速度16 km/h下排种合格指数较勺轮式排种器提升11.66个百分点，种子破损率降低0.21个百分点，证明基于分级限姿导引侧充−径向突变辅助清种的双列限姿导引侧充式排种器可有效确保排种器高速作业时充、清种效果可靠，对比传统的勺轮式排种器对高速作业的适应性更强。
表 6 排种器性能对比试验结果Table 6. Results of comparative test on the performance of seed planter类型
Type作业速度/(km·h−1)
Working speed评价指标/%
Evaluation index合格指数
Qualified index重播指数
Replay index漏播指数
Leakage index破损率
Breakage rate双列限姿导引侧充式排种器
Double-row seed metering device with guided
side-filling and posture control8 92.44 4.37 3.19 0.13 10 92.53 4.65 2.82 0.19 12 91.45 5.05 3.50 0.25 14 91.04 5.52 3.44 0.29 16 90.57 6.01 3.42 0.36 勺轮式排种器
Spoon-wheel seed metering device8 90.34 6.65 3.01 0.16 10 88.12 7.05 4.83 0.25 12 84.53 8.67 6.80 0.33 14 81.58 10.33 8.09 0.41 16 78.91 11.06 10.03 0.57 5. 结论
1)本文针对现有的机械式玉米排种器高速作业性能较差的问题，设计了一种基于分级限姿导引侧充−径向突变辅助清种的双列限姿导引侧充式玉米高速精量排种器。通过理论分析和计算，结合种子形状特征尺寸与种子填充姿态对关键部件参数进行了设计。
2)基于EDEM和RecurDyn的耦合仿真的正交试验结果表明，排种器在型孔插件顶部宽度W为4.65 mm、护种盘后侧高度H为5.10 mm、护种盘底部角度α为81.86°的条件下，排种器的性能达到最优，排种器的合格指数、重播指数和漏播指数分别为92.18%、5.13% 和2.69%。
3)台架性能试验的结果表明，在作业速度12 km/h 时，排种合格指数为91.45%，与仿真试验的相对误差为0.79%，仿真优化试验结果可靠；在作业速度8~16 km/h时，排种器的合格指数均能达到90%以上，重播指数均低于7%，漏播指数均低于4%，满足玉米精量播种的技术要求；双列限姿导引侧充式排种器的高速作业性能优于传统的勺轮式排种器。证明采用限姿板分级限制引导种子姿态和径向突变结构辅助清种能有效提升其高速作业时的精量播种能力，使其具备良好的高速适应性及稳定性。
-
表 1 幼虫期短时高温处理对井上蛀果斑螟幼虫存活率的影响1)
Table 1 Effects of short-term exposure of Assara inouei larvae to high temperature on their survival rates
虫期
Developmental staget(处理)/h
Treatment time存活率/% Survival rate 32 ℃ 35 ℃ 38 ℃ 41 ℃ 44 ℃ 初孵幼虫
Neonate larvae1 100.00aA 100.00aA 92.62aA 52.14aB 0.00aC 2 100.00aA 98.56aA 89.82aB 20.39bC 0.00aD 4 100.00aA 88.32bB 72.14bC 12.14cD 0.00aE 老熟幼虫
Mature larvae1 100.00aA 100.00aA 100.00aA 82.25aB 11.80aC 2 100.00aA 100.00aA 100.00aA 67.62bB 7.51aC 4 100.00aA 100.00aA 91.57bB 45.27cC 0.00bD 1) 相同虫期同列数据后不同小写字母表示不同处理时间存在显著差异，同行数据后不同大写字母表示不同温度间存在显著差异(P<0.05，Duncan’s法)
1) Different lowercase letters within the same column at the same larval stage indicate significant difference among different treatment time, and different uppercase letters within the same row indicate significant difference among different treatment temperature (P<0.05, Duncan’s test)
下载: 导出CSV
表 2 幼虫期短时高温处理对井上蛀果斑螟初孵幼虫发育历期及存活率的影响1)
Table 2 Effects of short-term exposure of neonate larvae to high temperature on the larval duration and survival rate of Assara inouei
t(处理)/h
Treatment time发育历期/d Larval duration 存活率/% Survival rate 35 ℃ 38 ℃ 41 ℃ 35 ℃ 38 ℃ 41 ℃ 0(CK) 17.56±1.97aA 17.56±1.97cA 17.56±1.97bA 68.24±2.68aA 68.24±2.68aA 68.24±2.68aA 1 16.32±1.52aB 19.16±2.43bA 20.19±1.09abA 64.38±4.51aA 40.06±2.47bB 33.35±1.02bC 2 16.08±1.34aB 21.25±2.30aA 19.83±2.12bA 53.74±2.34bA 37.83±2.39bcB 26.77±2.11cC 4 15.69±1.21aB 20.84±1.23abA 21.74±1.34aA 48.86±1.78cA 34.51±0.76cB 22.24±0.64cC 1) 表中数据为平均值 ± 标准差，对同一参数，同列数据后不同小写字母表示不同处理时间存在显著差异，同行数据后不同大写字母表示不同温度间存在显著差异(P<0.05，Duncan’s法)
1) Datum in the table is mean ± SE, for the same parameter, different lowercase letters within the same column indicate significant difference among different treatment time, and different uppercase letters within the same row indicate significant difference among different treatment temperature(P<0.05, Duncan’s test)
下载: 导出CSV
表 3 井上蛀果斑螟幼虫期短时高温处理对化蛹率的影响1)
Table 3 Effect of short-term exposure of larvae to high temperature on the pupation rate of Assara inouei
处理 Treatment 化蛹率/% Pupation rate θ/℃
t/h
初孵幼虫
Neonate larvae老熟幼虫
Mature larvae25 (CK) 87.72±4.95a 87.72±4.95a 35 1 89.35±4.56a 86.49±5.43a 2 85.94±5.31a 83.15±3.64ab 4 87.13±6.37a 78.89±5.20bc 38 1 88.32±5.44a 74.72±6.31c 2 86.09±6.32a 75.40±5.96c 4 85.41±2.56a 67.67±3.67d 41 1 86.38±4.27a 70.44±4.19d 2 88.84±5.19a 58.98±6.04e 4 85.63±3.38a 50.53±4.77f 1) 表中数据为平均值 ± 标准差，同列数据后不同小字母表示差异显著(P<0.05，Duncan’s法)
1) Datum in the table is mean ± SE, different lowercase letters within the same column indicate significant difference (P<0.05, Duncan’s test)
下载: 导出CSV
表 4 井上蛀果斑幼虫期短时高温处理对蛹质量、蛹历期及羽化率的影响1)
Table 4 Effects of short-term exposure of larvae to high temperature on the pupal weight, pupa duration and emergence rate of Assara inouei
处理 Treatment 蛹质量/mg Pupal weight 蛹历期/d Pupa duration 羽化率/% Emergence rate θ/℃
t/h
初孵幼虫
Neonate larvae老熟幼虫
Mature larvae初孵幼虫
Neonate larvae老熟幼虫
Mature larvae初孵幼虫
Neonate larvae老熟幼虫
Mature larvae25 (CK) 15.65±0.67a 15.65±0.67a 7.76±0.55a 7.76±0.55b 94.48±2.35a 94.48±2.35a 35 1 14.67±0.72a 14.43±1.21ab 6.82±0.48a 7.20±1.58b 93.32±4.57a 91.74±5.31ab 2 15.91±0.81a 15.07±0.94a 7.74±0.61a 7.69±1.57b 94.74±4.32a 90.53±3.67ab 4 13.89±0.76a 13.91±0.78ab 7.63±0.39a 7.65±1.49b 91.90±5.38a 88.57±4.45b 38 1 14.26±1.13a 14.60±080ab 7.72±0.60a 8.47±2.54ab 93.37±3.48a 83.79±5.70c 2 14.75±0.94a 13.87±0.75ab 8.25±0.42a 8.78±2.83ab 93.59±5.31a 84.31±2.61c 4 13.79±0.98a 12.67±0.54bc 8.41±0.79a 9.31±2.70a 92.86±3.35a 77.62±3.70d 41 1 14.61±0.76a 13.14±0.62bc 7.95±0.95a 9.04±2.63a 92.24±4.06a 79.85±3.66cd 2 13.94±0.88a 12.85±0.83bc 8.04±0.72a 9.23±2.71a 92.73±4.58a 66.76±4.22e 4 13.85±1.01a 11.28±0.87c 7.79±0.86a 9.18±1.48a 91.31±2.34a 68.39±3.38e 1) 表中数据为平均值 ± 标准差，同列数据后不同小字母表示差异显著(P<0.05，Duncan’s法)
1) Datum in the table is mean ± SE, different lowercase letters within the same column indicate significant difference (P<0.05, Duncan’s test)
下载: 导出CSV
表 5 井上蛀果斑幼虫期短时高温处理对成虫产卵期及寿命的影响1)
Table 5 Effects of short-term exposure of larvae to high temperature on the oviposition period and longevity of Assara inouei adult
处理
Treatment产卵前期/d
Pre-oviposition period产卵期/d
Oviposition period雌虫寿命/d
Female longevity雄虫寿命/d
Male longevityθ/℃
t/h
初孵幼虫
Neonate
larvae老熟幼虫
Mature
larvae初孵幼虫
Neonate
larvae老熟幼虫
Mature
larvae初孵幼虫
Neonate
larvae老熟幼虫
Mature
larvae初孵幼虫
Neonate
larvae老熟幼虫
Mature
larvae25 (CK) 1.71±0.02a 1.71±0.02a 8.69±0.74a 8.69±0.74 a 14.47±1.73a 14.47±1.73a 12.59±1.46a 12.59±1.46a 35 1 1.68±0.03a 1.62±0.04a 7.31±0.85a 8.31±0.85a 13.38±1.18a 15.14±2.14a 11.87±1.33a 12.37±1.36a 2 1.84±0.05a 2.10±0.02a 7.55±1.09a 8.14±1.09a 15.11±1.36a 13.81±1.59ab 12.69±1.78a 11.76±1.80a 4 1.56±0.03a 1.74±0.04a 6.57±0.67a 7.63±0.67abc 14.28±1.55a 13.44±1.37ab 11.77±0.89a 13.04±1.93a 38 1 2.05±0.03a 1.60±0.00a 6.86±0.98a 6.48±0.98abc 13.71±1.42a 15.06±2.00a 12.28±1.74a 12.68±1.38a 2 1.79±0.02a 1.77±0.05a 7.25±0.54a 5.59±0.54c 13.69±1.97a 14.87±1.42a 11.12±1.62a 12.31±1.55a 4 2.04±0.05a 2.08±0.04a 7.84±0.71a 6.25±0.71bc 14.23±2.18a 11.90±1.63bc 10.94±0.90a 11.84±1.63a 41 1 1.73±0.06a 1.69±0.03a 6.63±0.66a 6.17±0.66bc 13.56±1.92a 12.78±1.81b 12.50±1.34a 12.00±1.49a 2 1.80±0.04a 1.84±0.02a 7.39±0.82a 5.80±0.82c 13.89±1.49a 11.06±1.76bc 11.53±1.12a 11.64±1.07a 4 1.89±0.07a 1.95±0.06a 6.76±0.47a 5.21±0.47c 14.04±1.85a 9.57±2.05c 11.67±1.37a 12.17±1.22a 1) 表中数据为平均值 ± 标准差，同列数据后不同小字母表示差异显著(P<0.05，Duncan’s法)
1) Datum in the table is mean ± SE, different lowercase letters within the same column indicate significant difference (P<0.05, Duncan’s test)
下载: 导出CSV
表 6 井上蛀果斑幼虫期短时高温处理对成虫平均单雌产卵量及卵孵化率的影响1)
Table 6 Effects of short-term exposure of larvae to high temperature on the number of eggs laid per female and egg hatching rate of Assara inouei adult
处理
Treatment平均单雌产卵量/粒
Number of eggs laid per female卵孵化率/%
Egg hatching rateθ/℃
t/h
初孵幼虫
Neonate larvae老熟幼虫
Mature larvae初孵幼虫
Neonate larvae老熟幼虫
Mature larvae25 (CK) 89.47±5.26a 89.47±5.26a 57.89±2.65a 57.89±3.65a 35 1 91.87±4.56a 84.45±4.57ab 60.14±1.54a 58.20±1.58a 2 88.69±4.34a 86.60±3.65ab 54.27±2.17a 64.69±2.57a 4 83.79±3.90a 85.09±3.32ab 56.63±2.23a 54.25±1.49a 38 1 90.28±4.96a 90.15±5.94a 59.72±3.75a 60.61±3.54a 2 91.12±4.62a 87.42±4.06a 62.25±2.69a 57.31±2.83a 4 86.34±3.98a 80.62±3.41ab 55.18±2.58a 63.74±2.70a 41 1 84.57±5.34a 76.88±2.94b 55.33±1.63a 54.24±2.63a 2 86.53±3.12a 75.09±4.51b 62.87±2.54a 59.23±2.71a 4 88.67±4.37a 64.30±5.37c 53.62±1.06a 58.48±1.48a 1) 表中数据为平均值 ± 标准差，同列数据后不同小字母表示差异显著(P<0.05，Duncan’s法)
1) Datum in the table is mean ± SE, different lowercase letters within the same column indicate significant difference (P<0.05, Duncan’s test)
下载: 导出CSV
-
[1] 何超, 沈登荣, 尹立红, 等. 不同石榴品种对井上蛀果斑螟生长发育和繁殖的影响[J]. 应用生态学报, 2017, 28(3): 935-940. [2] PAEK M K, CHA J Y, BAE Y S. Two new records of phycitine moths (Lepidoptera, Pyralidae) in Korea[J]. Korean J Entomol, 2002, 32(2): 95-98. doi: 10.1111/j.1748-5967.2002.tb00017.x
[3] DU Y L, LI H H, WANG S X. A taxonomic study on the genus Assara walker from China (Lepidoptera: Pyralidae: Phycitinae)[J]. Acta Zootaxonomica Sinica, 2002, 27(1): 8-19.
[4] 白玲玲, 张祖兵, 杨仕生, 等. 云南石榴新纪录害虫井上蛀果斑螟的形态学及种群动态特征[J]. 云南农业大学学报(自然科学), 2005, 20(2): 183-187. [5] 秦卓, 蒋智林, 李正跃, 等. 石榴不同部位对井上蛀果螟(Assara inoue)产卵行为的影响[J]. 云南大学学报(自然科学), 2008, 30(S1): 123-126. [6] 邵淑霞, 李春艳, 杨仕生, 等. 井上蛀果斑螟幼虫和蛹的形态描述及生物学特性[J]. 云南农业大学学报(自然科学), 2008, 23(1): 22-24. [7] 何超, 沈登荣, 尹立红, 等. 井上蛀果斑螟成虫饲养密度对其寿命及繁殖力的影响[J]. 西北农业学报, 2017, 26(6): 950-955. doi: 10.7606/j.issn.1004-1389.2017.06.020 [8] CUI Y D, DU Y Z, LU M X, et al. Antioxidant responses of Chilo suppressalis (Lepidoptera: Pyralidae) larvae exposed to thermal stress[J]. J Therm Biol, 2011, 36(5): 292-297. doi: 10.1016/j.jtherbio.2011.04.003
[9] MIRONIDIS G K, SAVOPOULOU-SOULTANI M. Effects of heat shock on survival and reproduction of Helicoverpa armigera (Lepidoptera: Noctuidae) adults[J]. J Therm Biol , 2010, 35(2): 59-69. doi: 10.1016/j.jtherbio.2009.11.001
[10] 蒋丰泽, 郑灵燕, 郭技星, 等. 温度对昆虫繁殖力的影响及其生理生化机制[J]. 环境昆虫学报, 2015, 37(3): 189-199. [11] RUKKE B A, AAK A, EDGAR K S. Mortality, temporary sterilization, and maternal effects of sublethal heat in bed bugs[J]. PLoS One, 2015, 10(5): e0127555. doi: 10.1371/journal.pone.0127555
[12] ZHAO M T, WANG Y, ZHOU Z S, et al. Effects of periodically repeated heat events on reproduction and ovary development of Agasicles hygrophila (Coleoptera: Chrysomelidae)[J]. J Econ Entomol, 2016, 109(4): 1586-1594. doi: 10.1093/jee/tow093
[13] 梁菲菲, 王振营, 何康来, 等. 幼虫期短时高温暴露对二点委夜蛾存活和繁殖的影响[J]. 昆虫学报, 2016, 59(9): 991-996. [14] 李定旭, 雷喜红, 徐艳彩, 等. 短时高温对桃小食心虫生长发育与繁殖的影响[J]. 昆虫学报, 2014, 57(2): 218-225. [15] 范秀娟, 陈丹, 孙志娟, 等. 短时高温对烟蚜生长发育、繁殖和取食行为的影响[J]. 昆虫学报, 2014, 57(10): 1188-1197. [16] 鲍晓文, 房爱省, 仵均祥, 等. 短期高温对梨小食心虫成虫生殖及寿命的影响[J]. 果树学报, 2019, 36(4): 486-492. [17] 姜姗, 李帅, 张彬, 等. 极端高温对西花蓟马存活、繁殖特性及体内海藻糖、山梨醇含量的影响[J]. 中国农业科学, 2016, 49(12): 2310-2321. doi: 10.3864/j.issn.0578-1752.2016.12.007 [18] 何超, 沈登荣, 尹立红, 等. 井上蛀果斑螟生物学特性研究[J]. 应用昆虫学报, 2017, 54(2): 292-297. [19] GENG S, JUNG C. Effect of temperature on longevity and fecundity of Phyllonorycter ringoniella (Lepidoptera: Gracillariidae) and its oviposition model[J]. J Asia-Pac Entomol, 2017, 20(4): 1294-1300. doi: 10.1016/j.aspen.2017.09.012
[20] 李国平, 封洪强, 黄博, 等. 短时高温暴露对绿盲蝽和中黑盲蝽存活及生殖的影响[J]. 生态学报, 2017, 37(11): 3939-3945. [21] 杨帅, 赵冰梅, 李广云, 等. 短时高温暴露对土耳其斯坦叶螨和截形叶螨的影响[J]. 昆虫学报, 2013, 56(3): 276-285. [22] 全玉东, 何康来, 王振营, 等. 短时极端高温对亚洲玉米螟卵、初孵幼虫和成虫的影响[J]. 植物保护学报, 2015, 42(6): 985-990. [23] MAHROOF R, ZHU K Y, SUBRAMANYAM B. Change in expression of heat shock proteins in Tribolium castaneum (Herbst) (Coleopteran: Tenebrionidae) in relation to developmental stage, exposure time, and temperature[J]. Ann Entomol Soc Am, 2005, 98(1): 100-107. doi: 10.1603/0013-8746(2005)098[0100:CIEOHS]2.0.CO;2
[24] 马罡, 马春森. 气候变化下极端高温对昆虫种群影响的研究进展[J]. 中国科学: 生命科学, 2016, 46(5): 556-564. [25] ZHU G D, XUE M, LUO Y, et al. Effects of short-term heat shock and physiological responses to heat stress in two Bradysia adults, Bradysia odoriphaga and Bradysia difformis[J]. Sci Rep, 2017, 7: 13381. doi: 10.1038/s41598-017-13560-4
[26] ZHANG W, CHANG X Q, HOFFMANN A, et al. Impact of hot events at different developmental stages of a moth: The closer to adult stage, the less reproductive output[J]. Sci Rep, 2015, 5: 10436. doi: 10.1038/srep10436
[27] PIYAPHONGKUL J, PRITCHARD J, BALE J. Heat stress impedes development and lowers fecundity of the brown planthopper Nilaparvata lugens (Stål)[J]. PLoS One, 2012, 7(10): e47413. doi: 10.1371/journal.pone.0047413
计量
- 文章访问数: 936
- HTML全文浏览量: 8
- PDF下载量: 753
