摘要
为解决目前笼养种鸭饲喂系统供料量不精确、饲喂过程中个体饲喂量可控程度较低导致的不利于种鸭育种选育研究的问题,设计了一种叶片式颗粒饲料饲喂器。分析饲喂器充料排料原理以及核心工作部件参数,以饲料通过间距、排料叶片螺旋升角、叶片转速为试验因素,并以单位时间饲喂器饲料排量的变异系数为评价指标建立三因素三水平正交试验,并结合实际养殖需求,分析各因素的影响大小以及其最优组合。仿真结果显示,各参数对排料均匀性的影响大小依次为饲料通过间距、排料叶片螺旋升角、叶片转速;当饲料通过间距为9.0 mm、叶片螺旋升角为71.5°、叶片转速为58.4 r/min时,饲喂器的排料变异系数有最小值18.66%。在最优参数组合下,饲喂器排料性能台架试验结果显示,饲喂器排料变异系数为18.74%,排料速度为109.26 g/s,排料预设量与实际量的最大相对误差为3.07%,满足种鸭供料需求。
中国是世界上最大的蛋鸭养殖和消费国,现有的蛋鸭存栏量超过2亿,占全球的90%以
孟庆军
本研究设计了一种叶片式种鸭饲料精量饲喂系统,以饲喂器的结构设计为核心研究对象,使用离散元仿真等工具,探索种鸭圆柱型颗粒饲料的充料机制,优化饲喂器主要结构参数和工作参数来提高叶片式饲料饲喂器的工作性能,为实现笼养家禽的精准化养殖提供参考。
种鸭精准饲喂系统的整体结构如

图1 饲喂系统结构示意图
Fig.1 Schematic diagram of the feeding system structure
1.饲料Feed;2.料塔Feed tower;3.支架Stand;4.饲料导管Feed conduit;5.笼具Cage;6.种鸭Breeding ducks;7.个体检测模块Individual detection module;8.采食余量模块Feed residual module;9.饲料槽Feed troug;10.条形码Bar code;11.扫码模块Code scanning module;12.鸭蛋Duck eggs;13.产蛋计数模块Egg counting module;14.饲喂器Feeder.
总饲料供应系统启动后,饲喂平台开始从笼具一端沿着笼具排列方向向另一端移动,在饲喂平台移动过程中:扫码模块读取条形码内容,在系统的预饲喂量数据库获取此条码对应笼具所设置的饲喂量,不同的饲喂量设置对应控制核心输出不同的PWM波信号时间,系统控制核心通过输出PWM波信号时间来控制饲喂器中步进电机的转动时间从而匹配不同的饲喂量,同时开启当前笼位的饲喂控制和信息采集工作使能;采食余量模块进行工作,记录采食数据;个体检测模块对个体信息进行识别,在识别到正常鸭个体时开启饲喂模块使能,控制电机的启动;产蛋计数模块进行鸭蛋识别和计数统计;OLED显示屏实时显示当前笼位各项工作数据。随着饲喂平台的继续移动,当扫描模块读取到下个条码识别时,系统完成当前1个笼位的信息采集和饲喂工作,同时开启下1个笼位对应的控制和数据采集工作。在读取到笼具列末端的终止信号后,饲喂系统采集工作停止,所有数据通过串口上传到上位机中。其中,残余的饲料的处理有2种方式,对应2种不同的采食量记录模式:(1)需要保留残余饲料,则不进行饲料处理,便在饲喂量的数据管理中,记录本次饲料供应量为残余量与添加量之和;(2)残余饲料需要立即取出,进行人工取出则饲料供应量为本次饲料添加量。上位机通过阿里云平台将数据上传并且生成数据报表,方便其他工作人员获取本批次的养殖信息,以便进行品种选育。
1)饲喂器排料工作原理。该系统可以通过扫描模块获取当前笼位的条形码数据并与饲喂预设量数据匹配,控制饲喂器上的步进电机工作,以实现饲喂量的精确控制。在此过程中饲喂器的排料性能决定了该系统中饲料饲喂量的控制精
(1) |

图2 饲喂器结构图
Fig.2 Structure diagram of the feeder
A:饲喂器结构Structure of feeder;B:饲喂器结构剖视图Section view of feeder structure;1.入料导管Feeding conduit;2.饲喂器外壳Shell of feede;3.步进电机Stepper motor;4.步进电机支架Stepper motor bracket;5.安装底板Mounting base plate;6.联轴器Coupling;7.排料扇叶Work fan blades;8.轴承Bearing;9.轴承座Housing for bearing.
控制信号由STM32F103单片机开发板输出PWM脉冲数来控制步进角的速度。饲料通过入料导管进入饲喂器,在重力和相互作用力下完成充料,通过排料叶片的旋转作用从出口排出,落入笼具的食槽中。
2)饲喂器主要参数设计。该饲喂器采用叶片螺旋输送结构,饲喂器的理论排量公式如下所示:
(2) |
(3) |
(4) |
式(
在饲喂养殖过程中,不同生长阶段、不同品种个体以及不同的培育模式,对饲喂量的需求都不同,需要装置在对应条件下实现不同排量的精准供给和快速调
1)主要仿真参数的确定。试验中使用的种鸭饲料外形为类圆柱体的固体压缩饲料,其平均直径为4.6 mm,长度为8~13 mm。实际饲喂过程中,饲料需要通过厂房外的螺旋输送机输送后到达料箱,再输送到饲喂器终端,但是在这2个输送阶段中饲料会发生部分断裂。因此,为还原真实饲喂过程,在对饲槽中的饲料长度进行采样计数统计后,在仿真中将饲料颗粒归纳设置为3种长度状态:1号短颗粒4.4 mm、2号中颗粒7.8 mm、3号长颗粒12.4 mm,在仿真中将所占总颗粒质量比分别对应设置为34.5%、22.9%、42.6%。仿真软件中颗粒模型如

图3 饲料颗粒模型
Fig.3 Feed pellet model
A:短颗粒 Short particles;B:中颗粒 Medium particles;C:长颗粒 Long particles.
饲喂器与饲料颗粒接触部分的材料为光敏树脂。在饲喂器工作过程中,物体之间光滑无粘附,因此设置仿真设计中颗粒与颗粒和颗粒与饲喂器之间接触模型为 Hertz-Mindlin 无滑移接触模型。通过查阅资料和相关试验测定
对象Subject | 属性Causality | 仿真参数 Simulation parameters |
---|---|---|
饲料颗粒 Feed pellets | 泊松比 Poisson’s ratio | 0.32 |
剪切模量 Pashear modulus |
7.57×1 | |
真密度/(kg/ True density |
1.07×1 | |
短饲料颗粒 Short feed pellets | 长度/mm Length | 4.20 |
中饲料颗粒 Medium feed pellets | 长度/mm Length | 7.80 |
长饲料颗粒 Long feed pellets | 长度/mm Length | 12.60 |
9400树脂/ 旋转扇叶 9400 resin/rotary fan blade | 泊松比Poisson’s ratio | 0.41 |
剪切模量/Pa Shear modulus |
2.20×1 | |
密度/(kg/ |
1.13×1 | |
颗粒-颗粒 Pellet-pellet |
碰撞恢复系数 Crash recovery factor | 0.24 |
静摩擦因素 Coefficient of static friction | 0.30 | |
动摩擦因数 Coefficient of kinetic friction | 0.10 | |
颗粒-扇叶 Pellet-fan |
碰撞恢复系数 Crash recovery factor | 0.32 |
静摩擦因素 Coefficient of static friction | 0.21 | |
动摩擦因数 Coefficient of kinetic friction | 0.05 |
2)仿真试验的设计。上述饲喂器结构使用SOLIDWORKS进行三维建模,完成建模、装配和调整后,保存为STEP格式文件导入EDEM中。为了提高仿真效率,较长的饲槽模型在仿真中使用了圆筒模型代替。最后所得的饲喂器模型结构和仿真过程示意图如

图4 仿真模型图
Fig.4 Simulation model diagram
A:饲喂器三维模型 3D model of feeder;B:饲喂器仿真工作图Feeder simulation working diagram.
为了实时记录排料仿真过程中饲料的排量,在出料口设置了颗粒通过质量监测站(

图5 试验中设置的饲料质量检测站
Fig.5 Feed quality sensor in simulation test
通过EDEM中进行模拟仿真试验,仿真中将在饲料导管设置颗粒生成工厂,饲料以150 g/s的质量流进入饲喂器内腔,以确保在仿真过程中对饲喂器饲料供应充足。仿真过程中与颗粒有接触的模型设置振幅0.2 mm、频率10 Hz、方向垂直地面、持续时间0.1 s、间隔时间0.4 s的正弦运动,以模拟饲喂过程中饲喂平台主电机运动时产生的震动。
为了准确还原饲喂器排料的过程,仿真试验流程分成4个阶段(

图6 饲喂器工作示意图
Fig.6 Diagram of the feeder in working condition
A:自然充料阶段 Natural seeding phase;B:排料充料阶段 Discharge seeding phase;C:饲喂工作阶段 Feeding work phase;D:监测等待阶段 Monitor the waiting phase.
育种饲喂中饲料每次平均供应量为95 g,故饲喂器结构中饲料通过间距设计为6 mm、排料叶片数为6片、叶片有效工作长度为20 mm、叶片螺旋升角为75°,叶片转速为60 r/min,工作持续1 s,20组平均饲喂量为96.5 g。在饲喂器排料性能检测试验阶段,发现叶片到内腔底部的饲料通过间距、排料叶片的螺旋升角和排料叶片的角速度等因素对排料速率、排料质量波动变化等排料结果影响较大。为探索结构参数对饲喂器排料功能的影响,本研究以上述3个因素为研究对象,探究3个因素对排料均匀性的影响。
饲料通过间距为饲喂器内腔中排料叶片最低处到饲喂器内腔底部最低的距离,如

图7 饲料通过间距示意图
Fig.7 Schematic diagram of the spacing
叶片螺旋升角示意图如
(5) |

图8 叶片螺旋升角示意图
Fig.8 Schematic diagram of the spiral rising angle
螺旋升角范围为0°~90°,升角越大,排料叶片接近直型叶片。随着排料叶片螺旋升角从90°开始减小,饲料颗粒受到的叶片轴向方向的力会从零开始增加,同时饲料受径向力开始减小。轴向力方向为大拇指所指的方向(根据右手螺旋法则,当四指与叶片旋转方向相同时)。轴向力使得饲料颗粒有垂直向内的运动趋势,使得当螺旋升角过小时,饲料颗粒有在饲喂器内壁聚集的可能性,从而影响饲喂器的排料性能。因此,在本研究中设计螺旋升角应当有最小值。试验中使用单片机向步进电机发送不同的脉冲频率信号以控制叶片转速,步进电机的转速公式为:
(6) |
为考察上述3个因素对变异系数的影响,设置三因素三水平正交试验。在预试验阶段,发现饲料通过间距和叶片转速对饲料排料影响效果明显,超出一定范围会使得饲喂器系统在短时间完不成指定饲喂量供应或者使得单位时间饲喂量过大不利于系统通过时间调节饲喂量。因此,选取饲料通过间距为2、10 mm作为水平最值,原转速的80%和120%作为水平最值。当叶片螺旋升角度为60°时,相邻2块叶片边缘的连线与旋转轴平行,使得轴向截面上,叶片占据饲料的体积均匀分布,能缓解传统直型扇叶的脉冲效应,故选取叶片螺旋升角60°与传统直型叶片螺旋升角90°作为水平最值。试验因素编码表如
水平 Level | X1饲料通过间距/mm Feed through spacing | X2叶片螺旋 升角/(°) Blade spiral angle of rise | X3叶片转速/ (r/min) Blade speed |
---|---|---|---|
-1 | 2 | 60 | 48 |
0 | 6 | 75 | 60 |
1 | 10 | 90 | 72 |
试验编号 Test No. | X1 | X2 | X3 | Y1/% |
---|---|---|---|---|
1 | -1 | -1 | 0 | 21.591 9 |
2 | 1 | -1 | 0 | 20.330 7 |
3 | -1 | 1 | 0 | 33.814 6 |
4 | 1 | 1 | 0 | 20.937 5 |
5 | -1 | 0 | -1 | 31.614 7 |
6 | 1 | 0 | -1 | 19.016 9 |
7 | -1 | 0 | 1 | 23.942 7 |
8 | 1 | 0 | 1 | 20.625 4 |
9 | 0 | -1 | -1 | 21.533 1 |
10 | 0 | 1 | -1 | 28.397 5 |
11 | 0 | -1 | 1 | 20.017 6 |
12 | 0 | 1 | 1 | 22.772 1 |
13 | 0 | 0 | 0 | 19.603 3 |
14 | 0 | 0 | 0 | 19.494 8 |
15 | 0 | 0 | 0 | 22.167 2 |
16 | 0 | 0 | 0 | 19.571 7 |
17 | 0 | 0 | 0 | 20.603 3 |
注: Note:X1:饲料通过间距Feed through spacing;X2:叶片螺旋升角Blade spiral rise angle;X3:叶片转速Blade speed;Y1:变异系数Coefficient of variation.
将
方差来源 Source of variance | 均方 Mean square | 自由度 Degrees of freedom | 平方和 Square sum | F | P |
---|---|---|---|---|---|
模型 Model | 0.030 1 | 9 | 0.003 3 | 33.07 |
<0.000 |
X1 | 0.011 3 | 1 | 0.011 3 | 111.78 |
<0.000 |
X2 | 0.006 3 | 1 | 0.006 3 | 62.37 |
<0.000 |
X3 | 0.002 2 | 1 | 0.002 2 | 21.58 |
0.002 |
X1X2 | 0.003 4 | 1 | 0.003 4 | 33.40 |
0.000 |
X1X3 | 0.002 2 | 1 | 0.002 2 | 21.32 |
0.002 |
X2X3 | 0.000 4 | 1 | 0.000 4 | 4.18 | 0.080 2 |
X1² | 0.002 1 | 1 | 0.002 1 | 21.11 |
0.002 |
X2² | 0.001 1 | 1 | 0.001 1 | 11.08 |
0.012 |
X3² | 0.000 7 | 1 | 0.000 7 | 6.63 |
0.036 |
残差 Residual | 0.000 7 | 7 | 0.000 1 | ||
失拟 Misfitting | 0.000 2 | 3 | 0.000 1 | 0.465 0 | 0.722 3 |
误差 Error | 0.000 5 | 4 | 0.000 1 | ||
总和 Total | 0.030 8 | 16 |
注: “**”表示该项极显著(P<0.01);“*”表示该项显著(P<0.05)。Note:“**” indicates that the item is highly significant(P<0.01);“*” indicates that the item is significant(P<0.05).
由

图9 交互作用对变异系数影响的响应面
Fig.9 Response surface of the influence of interaction on coefficient of variation
去除不显著项,得到单位时间排量变异系数二次多项式回归模型如
Y1=20.29-3.76X1+2.81X2-1.65X3-2.90X1X2+
2.32X1X3 +2.25X | (7) |
设计的饲喂系统通过控制排料时间来控制饲料排量,饲喂时间中位数设置为1 s。根据饲喂供料需求,需要饲喂器1 s的饲料排量在60~120 g。重复变异系数函数求解过程,根据方差分析和拟合结果得到饲料通过间距、叶片螺旋升角、叶片转速与排量呈一次关系,拟合方程模型显著,方程拟合性好,拟合函数如下所示。
Y2=4.937+0.768X1+0.120X2+0.863X3 | (8) |
为探索该结构下饲喂器最佳结构参数,根据拟合方程,以饲料通过间距、叶片螺旋升角、叶片转速为自变量,以变异系数最小为优化目标,并结合实际饲喂排量需求条件建立数学模型。
(9) |
使用Design Expert中Numerical Optimization 模块进行目标参数优化求解Y1,得到当饲料通过间距为9.0 mm、叶片螺旋升角为71.5°、叶片转速为58.4 r/min时,变异系数有最小值18.66%,此时每0.05 s平均饲料排量为5.5 g,满足生产需求。为了验证仿真试验的准确性,进行饲喂器排种试验,

图10 验证试验台架
Fig.10 Verification test bench
1.电源Power supply;2.饲料导管Feed conduit;3.饲喂器 Feeders;4.步进电机Stepper motor;5.步进电机控制器 Stepper motor controller;6.单片机Microcontroller;7.台架 Racks;8.烧杯 Beaker;9.高精度电子秤High precision electronic scale.
使用STM32控制步进电机转动,排料转动时间设置为0.05 s,转动停止间隔为2 s,排料转动和停止重复进行20次,使总排料时间为1 s。通过电子秤测量记录烧杯内饲料质量变化,得到每次排料的质量,计算落入烧杯中饲料的平均质量、质量变异系数、饲喂总质量以及实际质量与预测质量之间的相对误差。同时进行变量排料试验,增加相同条件下总饲喂时间为分别0.75、1.25 s的验证试验,总时间不同的每组试验进行10次,取10次结果的平均值。试验结果如
总饲喂时间/s Feeding time | 单位时间平均饲喂质量/g Average feeding quality | 质量变异 系数/% Quality coefficient of variation | 总质量/g Total mass | 相对误差/% Relative error |
---|---|---|---|---|
0.75 | 5.61 | 18.88 | 84.16 | 2.01 |
1.00 | 5.46 | 18.74 | 109.26 | 0.67 |
1.25 | 5.67 | 18.67 | 141.73 | 3.07 |
本研究设计一套笼养种鸭饲喂监控信息采集系统,该系统能通过对笼养种鸭的身份识别后对种鸭个体实现精准定量饲喂,同时能采集种鸭的健康信息以及采食量和产蛋量等重要的养殖育种信息,提高了育种研究的智能化精确化工作水平。
本研究设计了一种基于单片机控制、以步进电机为控制终端的叶片式精量饲喂器,着重对饲喂器的结构设计和参数优化进行了研究,通过正交试验分析了饲料通过间距、叶片螺旋升角、叶片转速对饲喂器排料均匀性的影响,结果表明:当饲料通过间距为9.0 mm、叶片螺旋升角为71.5°、叶片转速为58.4 r/min时,系统能获取最优工作性能参数,此时饲喂变异系数为18.66%,排料速率为110 g/s,各因素对变异系数的影响主次顺序为饲料通过间距、叶片螺旋升角、叶片转速。通过台架试验对预测结果进行试验验证,得到饲喂变异系数为18.74%,饲喂质量流为109.26 g/s,多组试验中预设量与实际量的最大相对误差值为3.07%,试验结果显示,该结构下饲喂器排料速率稳定,在1.5 s的饲喂时间内,饲喂器能实现单次0~160 g内可控的饲喂量供应,能满足种鸭选育研究所需要的供料要求。
研究结果表明,本研究设计的种鸭选育饲喂系统可满足种鸭饲喂过程可控、精量的供料需求,下一步考虑对饲喂器结构中叶片数目、叶片直径等其他设计参数进行研究,通过优化其他参数来进一步提升饲喂器的饲喂性能。
参考文献 References
赵静,张红凤.我国家禽养殖现状与发展趋势研究[J].现代畜牧科技,2020(12):21-22.ZHAO J,ZHANG H F.Study on the present situation and development trend of poultry breeding in China[J].Modern animal husbandry science & technology,2020(12):21-22(in Chinese). [百度学术]
邓见文,宁中华.笼养和福利散养模式下蛋鸡相关性能比较分析[J].中国家禽,2022,44(9):74-79.DENG J W,NING Z H.Comparative analysis on performances of laying hens under cage and welfare free range modes[J].China poultry,2022,44(9):74-79(in Chinese with English abstract). [百度学术]
匡伟,郗正林,伍冠锁,等.蛋鸭笼养技术实践与研究进展[J].家禽科学,2020(5):54-56.KUANG W,XI Z L,WU G S,et al.Practice and improvement of cage rearing technique for laying duck[J].Poultry science,2020(5):54-56(in Chinese with English abstract). [百度学术]
陈斌丹,张科杰,李柳萌,等.蛋鸭笼养现状与发展趋势[J].浙江农业科学,2022,63(8):1868-1870.CHEN B D, ZHANG K J ,LI L M,et al.Present status and development trend of egg duck cage rearing[J].Journal of Zhejiang agricultural sciences,2022,63(8):1868-1870(in Chinese with English abstract). [百度学术]
王玉泉.蛋鸭层叠笼养设备及控制系统的设计与研究[D].武汉:华中农业大学,2020.WANG Y Q.Design and research of cascading cage-rearing equipment and control system for laying duck[D].Wuhan:Huazhong Agricultural University,2020(in Chinese with English abstract). [百度学术]
王开云,钟日开,罗土玉,等.鸭只生产性能测定系统设计和试验[J].黑龙江畜牧兽医,2022(24):55-61.WANG K Y,ZHONG R K,LUO T Y,et al.Design and experiment of production performance testing system for ducks[J].Heilongjiang animal science and veterinary medicine,2022(24):55-61(in Chinese with English abstract). [百度学术]
孟庆军,邱俊柱.基于离散元的自走式家禽饲料喂料机设计及仿真[J].林业机械与木工设备,2022,50(3):48-51.MENG Q J,QIU J Z.Design and simulation of self-propelled poultry feeder based on discrete element[J].Forestry machinery & woodworking equipment,2022,50(3):48-51(in Chinese with English abstract). [百度学术]
戚江涛,蒙贺伟,坎杂,等.基于EDEM的双螺旋奶牛饲喂装置给料性能分析与试验[J].农业工程学报,2017,33(24):65-71.QI J T,MENG H W,KAN Z,et al.Analysis and test of feeding performance of dual-spiral cow feeding device based on EDEM[J].Transactions of the CSAE,2017,33(24):65-71(in Chinese with English abstract). [百度学术]
聂传斌,孙卫明,张燕军,等.基于DEM仿真的种鹅精准饲喂装置设计优选研究[J].中国农机化学报,2022,43(6):67-75.NIE C B,SUN W M,ZHANG Y J,et al.Design and optimization of breeding goose precision feeding device based on DEM simulation[J].Journal of Chinese agricultural mechanization,2022,43(6):67-75(in Chinese with English abstract). [百度学术]
杜华兴,李久熙,李丽华.基于EDEM的平养育雏喂料装置设计与试验[J].中国农机化学报,2022,43(2):60-66.DU H X,LI J X,LI L H.Design and experiment of a feeding device for flat rearing chicks based on EDEM[J].Journal of Chinese agricultural mechanization,2022,43(2):60-66(in Chinese with English abstract). [百度学术]
赵思琪,丁为民,赵三琴,等.基于EDEM-Fluent的气动式鱼塘投饲机性能优化[J].农业机械学报,2019,50(9):130-139.ZHAO S Q,DING W M,ZHAO S Q,et al.Performance analysis and optimization of pneumatic fishpond feeder based on EDEM-fluent coupling[J].Transactions of the CSAM,2019,50(9):130-139(in Chinese with English abstract). [百度学术]
刘朝贤.猪用饲喂器拨片式供料装置设计与试验研究[D].哈尔滨:东北农业大学,2021.LIU C X.Design and experiment of the paddle-unloading feeding equipment for pig-feeder[D].Harbin:Northeast Agricultural University,2021(in Chinese with English abstract). [百度学术]
吴艳,皮劲松,梁振华,等.农湖2号蛋鸭配套系选育研究[J].中国畜牧杂志,2023,59(3):120-122.WU Y,PI J S,LIANG Z H,et al.Study on breeding of Longhu No.2 laying duck mating line[J].Chinese journal of animal science,2023,59(3):120-122(in Chinese). [百度学术]
牛智有,孔宪锐,沈柏胜,等.颗粒饲料破损离散元仿真参数标定[J].农业机械学报,2022,53(7):132-140.NIU Z Y,KONG X R,SHEN B S,et al.Parameters calibration of discrete element simulation for pellet feed attrition[J].Transactions of the CSAM,2022,53(7):132-140(in Chinese with English abstract). [百度学术]
郭三琴,万鹏,杨俊,等.基于离散元的鱼饲料仿真参数标定与试验[J].华中农业大学学报,2022,41(6):286-295.GUO S Q,WAN P,YANG J,et al.Discrete element based calibration and test of simulation parameters for fish feed[J].Journal of Huazhong Agricultural University,2022,41(6):286-295(in Chinese with English abstract). [百度学术]
赖庆辉,贾广鑫,苏微,等.基于DEM-MBD耦合的链勺式人参精密排种器研究[J].农业机械学报,2022,53(3):91-104.LAI Q H,JIA G X,SU W,et al.Design and test of chain-spoon type precision seed-metering device for ginseng based on DEM-MBD coupling[J].Transactions of the CSAM,2022,53(3):91-104(in Chinese with English abstract). [百度学术]
韩丹丹,何彬,周毅,等.气吸式密植精量排种器的设计与试验[J].华中农业大学学报,2023,42(1):237-247.HAN D D,HE B,ZHOU Y,et al.Design and experiment of air suction type dense planting precision seed metering device[J].Journal of Huazhong Agricultural University,2023,42(1):237-247(in Chinese with English abstract). [百度学术]
高全伟.蛋鸭不同品种及其杂交后代的生长性能及限饲对开产日龄的影响[D].南京:南京农业大学,2016.GAO Q W.The ducks of different varieties and their hybrid offspring growth performance and feed restriction on the influence of day age[D].Nanjing:Nanjing Agricultural University,2016(in Chinese with English abstract). [百度学术]