Wood densities of six hardwood tree species in Eastern Guangdong and Influencing factors
-
摘要:目的
分析不同影响因子对木材密度的影响，为木材密度的良种选育和碳汇计量提供理论和数据支撑。
方法在康禾自然保护区样地中，通过群落分析筛选出6个优势树种，测定其木材密度(包括生材密度和基本木材密度)，通过多种统计软件分析不同树种、不同径阶和不同坡向对木材密度的影响，以及树木性状因子与木材密度的相关关系。
结果木材密度在树种间的差异性极显著，生材密度、基本木材密度和含水率平均值最大的分别是米槠Castanopsis carlesii 、红锥C. hystrix 和米槠，其数值分别为1.10 g·cm-3、0.53 g·cm-3和122.84%。木材密度的径阶差异性分析结果显示，除了栲C. fargesii以外，其他树种的木材密度在径阶间变化并不显著。木材密度在坡向上的差异性分析结果表明，6个树种的木材密度在坡向间的差异均不显著。非参数相关分析结果表明，胸径、树高与生材密度的相关性强，而与基本木材密度的相关性弱。
结论红锥含水率低，基本木材密度高，具备优质木材的条件。木材密度随着径阶的增加有上升的趋势，但在坡向上没有规律性变化。生材密度、基本木材密度与胸径、树高的关系并不一致，可能是因为生材密度受到含水率增加的影响。
Abstract:ObjectiveThis study aims to provide the theoretical and data supports for the breeding of wood density and carbon sink accounting by analyzing the effects of different impact factors on wood density.
MethodSix tree species were selected through the analysis of communities in sample plots of Kanghe nature reserve, and their wood densities were measured. This study analyzed the effects of different species, different diameter grades and different slope aspects on wood density, and the correlation between tree traits and wood density using a variety of statistical softwares.
ResultThere were significant differences in wood density among different tree species. The highest averages of green densities (1.10 g·cm-3), basic densities (0.53 g·cm-3) and the moisture contents (122.84%) were found in Castanopsis carlesii, C. hystrix and C. carlesii respectively. There was no significant difference in wood density among different diameter grades for all species except C. fargesii. There was no significant difference in wood density among different slope aspects for all species. The nonparametric Spearman's rank correlation analysis showed that both diameter at breast height and tree height were significantly correlated with green density but not with basic density.
ConclusionC. hystrix produces high quality wood with low moisture content and high basic density. Wood density has a tendency to rise with the increase of diameter grade, but has no regular change with the change of slope aspect. The relationship between diameter at breast height, tree height and green density, and the relationship between diameter at breast height, tree height and basic density are different probably because green density is affected by the moisture content.
-
Keywords:
- wood density /
- green density /
- basic density /
- diameter grade /
- slope aspect /
- hardwood tree /
- influencing factor
-
设施大棚种植业的发展能够有效保证蔬菜、瓜果的淡季供应，具有巨大的社会效益和经济效益[1]。中国设施农业机械化发展水平相对滞后，综合机械化水平仅为33%[2-3]。施肥是设施大棚种植作业中最重要的环节之一，施撒有机肥能够有效改善土质、减少温室气体排放、改良作物品质、提高作物产量[4-8]。但是，目前大棚中有机肥的撒施仍然主要依靠人工，劳动强度大、作业效率低[9]，要实现高效高质的设施农业生产，必须依靠机械施肥，排肥装置是施撒有机肥的关键部件[10]。现阶段，国内外对于有机肥施肥机械的研究很多[11-13]，但专用于设施大棚的施肥机械鲜见介绍。为了提高我国设施大棚施肥机械化水平，有必要设计大棚有机肥排肥装置。目前用于有机肥施撒的主要排肥装置及原理为：离心圆盘式，利用圆盘的离心力抛撒肥料[14]；桨叶式，通过桨片高速旋转将肥料抛撒到田间[15]；锤片式，通过螺旋绞龙运送肥料至卸料口，再用高速运动的锤片将肥料撕裂、粉碎，自上而下抛出[16]；抛撒转轴式，通过排料口挡板调节落料量，落料通过旋转转轴抛撒而出[17]。以上排肥装置都只适用于大型农场或大田[18]，施撒效率高，排肥原理都是利用排肥器高速旋转的离心力惯性抛撒有机肥，肥料经排肥器排出过程中会有大量泄漏，导致排肥量不稳定。
刮板输送装置对散体物料良好的输送特性保证了散体物料运输、排施过程的稳定性和定量可控性[19-21]，但将刮板输送装置应用于有机肥施撒并不多见。本研究综合考虑刮板输送装置的机构特性和输送原理，设计一种刮板输送式有机肥排肥装置，对其关键结构参数进行设计和分析，搭建排肥器性能试验平台，对关键参数进行优化，以期为整机的设计提供参考依据。
1. 整体结构与工作原理
1.1 大棚有机肥施撒农艺要求
以苏南地区草莓大棚为例，该地区草莓种植大棚多为8332型大棚(图1)，棚高3.2 m，棚宽8 m，长 80 m，大棚门宽1.2 m，高度为1.75 m，大棚内作业空间有限，一般采用人工作业[22-23]。为了促进草莓开花结果，一般在草莓种植之前，在整个大棚内进行基肥施撒作业，施撒有机肥量为1.50~2.25 kg·m−2，有机肥种类为腐熟的肥类有机肥。施肥作业先将有机肥撒于地表，后用旋耕机将有机肥旋入土壤。
1.2 排肥装置整体结构
刮板输送式排肥装置主要由承料箱、挡肥板升降装置、刮板输送式排肥器、传动装置和控制器组成。承料箱由箱体和挡肥板组成，其中挡肥板安装在箱体两边的滑道中，可以上下升降移动。挡肥板升降装置主要由齿轮轴、齿轮和齿条组成，齿条焊接在挡肥板上，依靠齿轮旋转运动带动齿条上下移动。刮板输送式排肥器由排肥轴、链轮、链条、链板、刮板等组成，刮板用沉头螺钉安装在链板上，链板与链板之间采用铰接。传动装置主要由步进电机、联轴器、三相异步电机和减速机组成，步进电机通过联轴器连接齿轮轴，减速机侧面轴孔连接三相异步电机轴，正面轴孔连接排肥轴。控制器包括步进电机控制器、三相异步电机变频器以及编码器，编码器通过联轴器安装在排肥轴上，实时测量排肥轴转速。整体结构如图2所示。
![]()
图 2 排肥装置结构图1：步进电机；2：步进电机安装座；3：联轴器1；4：齿条；5：齿轮；6：齿轮轴；7：挡肥板；8：箱体；9：轴承座1；10：三相异步电机；11：减速机；12：三相异步电机安装座；13：支架；14：排肥轴；15：轴承座2；16：链轮；17：链条；18：链板；19：刮板；20：挡片；21：联轴器2；22：编码器Figure 2. Structure diagram of fertilizer discharging device1: Stepping motor; 2: Stepping motor mounting seat; 3: Coupling 1; 4: Rack and pinion; 5: Gear; 6: Gear shaft; 7: Fertilizer plate; 8: Box body; 9: Bearing seat 1; 10: AC motor; 11: Reducer; 12: AC motor mounting seat; 13: Bracket; 14: Sprocket shaft; 15: Bearing seat 2; 16: Sprocket; 17: Chain; 18: Chain plate; 19: Scraper; 20: Baffle; 21: Coupling 2; 22: Encoder1.3 排肥装置工作原理
步进电机通过联轴器将驱动力传输到齿轮轴，固定在齿轮轴上的齿轮带动齿条升降，把电机的旋转运动转化为直线运动从而驱动挡肥板升降，精准控制排肥装置开度。三相异步电机的动力通过减速机驱动排肥轴，带动链条、链板和刮板向排肥口处运动，刮板可把料箱底部的物料刮走，上方物料在自身重力作用下，克服摩擦力和内聚力而下落，继续受刮板作用力，使物料连续不断地稳定排出[24]。实施排肥作业时，分别调节步进电机控制器和三相异步电机变频器旋钮，设置对应排肥器开度和排肥轴转速。根据设定的开度值，步进电机转动相应的步距角，调节挡肥板开度；根据设定的转速值，三相异步电机带动排肥轴实施排肥作业。
2. 排肥器主要参数
2.1 排肥量的确定
排肥轴的单转排肥量是衡量排肥装置排肥性能的重要指标[25]，排肥器单转排肥量的稳定性决定了整个排肥装置的排肥稳定性。不考虑轴承座和支架两侧壁摩擦力影响，其排肥轴理论单转排肥量用下式计算：
$$ q = {\text{π}} dHb\left( {\frac{a}{{a + {\textit{z}}}}} \right)\rho, $$ (1) 式中，q表示输送链轮旋转一周排出的肥料，kg·r−1；H为排肥器的开度，m；b为排肥刮板的长度，m；a为刮板间距，m；z为刮板厚度，m；d为输送链轮直径，m；ρ为有机肥的密度，kg·m−3。
2.2 刮板间距与高度
有机肥在承料箱中运动时，呈连续的均质流输送，整个肥料流分为上下两层，下层有机肥高h1，刮板上方的有机肥高h2，任取1个刮板间距的有机肥为1个单元体，参考文献[26]，对刮板上层的有机肥进行受力分析，结果如图3所示。
![]()
图 3 单元体有机肥的受力分析a：刮板间距；h1：刮板间有机肥高度；h2：刮板上方的有机肥高度；h：输送有机肥整体高度；b：排肥刮板的长度；FS：有机肥与承料箱侧壁的摩擦力；Fr：有机肥体与承料箱底板的外摩擦力；FX：有机肥对承料箱壁的侧压力；FA：后刮板对单元体有机肥的推力；FB：前刮板对单元体有机肥的压力；G：上层有机肥重力；A：后刮板；B：前刮板；v：刮板前进速度Figure 3. Force analysis of unit organic fertilizera: Scraper spacing; h1: Height of organic fertilizer between scrapers; h2: Height of organic fertilizer above the scraper; H: Overall height of transported organic fertilizer; b: Legth of the scraper; FS: Friction between organic fertilizer and the side wall of the carrier box; Fr: External friction between the organic fertilizer body and the bottom plate of the carrier box; FX: Lateral pressure of organic fertilizer on the wall of the carrier box; FA: Thrust of the rear scraper on unit organic fertilizer; FB: Pressure of the front scraper on unit organic fertilizer; G: Upper organic fertilizer gravity; A: Back scraper; B: Front scraper; v: Scraper forward speed由散体力学理论可得，由于上层有机肥对承料箱壁有侧向压力作用，其大小随有机肥深度呈线性增加，并作三角形分布，其值为该点深度处的正压力与侧压系数之积，则承料箱每侧边单位面积的平均侧向力
$({F_{\rm{X}}} )$ 为：$$ {F_{\rm{X}}} = \frac{1}{2}\frac{g}{{ba}}n = \frac{1}{2}\rho g{h_2}{\tan ^2}\left( {\frac{{\text{π}}}{4} - \frac{\theta }{2}} \right), $$ (2) 式中，n为侧压系数；g为重力加速度，m·s−2；θ为有机肥的休止角，(°)。
有机肥对承料箱侧壁产生一个作用在1/3h2高度处的外摩擦力(FS），为：
$$ {F_{\rm{S}}} = \frac{{h_2^2a{f_1}}}{2}\rho g{\tan ^2}\left( {\frac{{\text{π}}}{4} - \frac{\theta }{2}} \right), $$ (3) 式中，f1为有机肥与承料箱的外摩擦系数。
如果假设在两刮板垂直平面之间的上层肥料流为一固定单元体，摩擦力FS使单元体承受一个弯矩，弯矩在刮板A 的上表面产生压应力，在刮板B的上表面产生一个大小相等、方向相反的拉应力，此应力将与上层有机肥在刮板A、B表面上产生的压应力叠加，使A压应力增加，B压应力减小，此时B内摩擦力产生的剪应力为上下层有机肥接触面上各处内摩擦力产生的剪应力的最小值。则B内摩擦力产生的剪应力(τB)为：
$$ {\tau _{\rm{B}}} = f\left( {{h_2}\rho g - \frac{{4{F_{\rm{S}}}{h_2}}}{{b{a^2}}}} \right), $$ (4) 式中，f为内摩擦系数。
只有使上、下层有机肥接触面上各处的内摩擦力产生的剪应力均大于有机肥和承料箱两侧外摩擦力作用于整个接触面上均匀分布的剪应力，才能保证有机肥能够整体运动，在没有挡肥板的情况下，倾仓而出，根据这一条件得：
$$ {\tau _{\rm{B}}} \geqslant \tau = \frac{{2{F_{\rm{S}}}}}{{ba}}, $$ (5) 式中，τ为有机肥和承料箱两侧外摩擦力作用于整个接触面上均匀分布的剪应力，N·m−2。
取式(4)的τB和式(3)的FS代入式(5)，得到h2的二次方程式：
$$ {h_2} \leqslant \dfrac{1}{{4f}}\left( {\sqrt {\dfrac{{8ab{f^2} + {a^2}{f_1}{{\tan }^2}\Bigg(\dfrac{{\text{π}}}{4} - \dfrac{\theta }{2}\Bigg)}}{{{f_1}{{\tan }^2}\Bigg(\dfrac{{\text{π}}}{4} - \dfrac{\theta }{2}\Bigg)}}} - a} \right)。 $$ (6) 经过试验测得有机肥休止角θ=47°，内摩擦系数f=tanθ=1.07，根据文献[27]，有机肥与承料箱的外摩擦系数f1=0.35，b为刮板长度，取b为1 m。刮板高度通过试验取为0.01 m，可以达到施肥量的要求。承料箱中有机肥的最大输送高度依据排肥装置最大开度设定并预留一定的承肥空间，承料箱高度取整为0.2 m。由式(6)可得h2和a的函数关系如图4所示。
![]()
图 4 肥料高度和刮板间距的关系Figure 4. The relationship between fertilizer height and scraper spacing由式(6)可得刮板间距的下极限值：
$$ {a_{\min }} = \dfrac{{2f{f_1}h_2^2{{\tan }^2}\left( {\dfrac{{\text{π}}}{4} - \dfrac{\theta }{2}} \right)}}{{bf - {h_2}{f_1}{{\tan }^2}\left( {\dfrac{{\text{π}}}{4} - \dfrac{\theta }{2}} \right)}}, $$ (7) 式中，
${a_{\min }} $ 是刮板间距(a)的最小值，当承料箱中有机肥的料层高度为最大h2时取到。对于一定的有机肥和承料箱宽度，a不能随意增大，因为a增大，将导致上层有机肥内部产生滑移，使有机肥上表层滞留在刮板后方，影响排肥稳定性。为避免在承料箱中形成滑移面，a必须有最大极限值。假设两刮板之间整体肥料为一固定单元体，不考虑挡肥板的阻力，将肥料体即将产生滑移面时的情况作为受力平衡情况，得：
$$ {F_{\rm{A}}} = {F_{\rm{r}}} + 2{F_{\rm{S}}} + {F_{\rm{B}}}, $$ (8) 式中，有机肥体与承料箱底板的外摩擦力(Fr)为：
$$ {F_{\rm{r}}} = b\left( {{h_1} + {h_2}} \right)a\rho g{f_1}。 $$ (9) 由于上层有机肥和承料箱侧壁外摩擦力Fs对上层肥料产生一个弯矩，该弯矩在前后刮板上表面分别产生拉、压应力，则前、后刮板上方的垂直压应力分别为
$$ {\sigma _{{\rm{A}}1}} = {h_2}\rho g + \frac{{4{F_{\rm{S}}}{h_2}}}{{ba}}, $$ (10) $$ {\sigma _{{\rm{B}}1}} = {h_2}\rho g - \frac{{4{F_{\rm{S}}}{h_2}}}{{ba}}, $$ (11) 式中，σA1
为后刮板上方的垂直压应力，N·m−2；σB1 为前刮板上方的垂直压应力，N·m−2。 刮板下方的垂直应力不受该弯矩的影响，则有：
$$ {\sigma _{{\rm{A}}2}} = {\sigma _{{\rm{B}}2}} = \left( {{h_1} + {h_2}} \right)\rho g , $$ (12) 式中，σA2为后刮板下方的垂直压应力，N·m−2；σB2为前刮板下方的垂直压应力，N·m−2。
用刮板上、下方的平均侧压力代替刮板对于单元体的作用力(FA、FB)，则FA、FB分别为：
$$ {F_{\rm{A}}} = \frac{1}{2}\left( {{\sigma _{{\rm{A}}1}} + {\sigma _{{\rm{A}}2}}} \right)\frac{1}{n}{h_1}b = \left[ {\left( {{h_2} + \frac{{{h_1}}}{2}} \right)\rho g + \frac{{2{F_{\rm{S}}}{h_2}}}{{ba}}} \right]\frac{1}{n}{h_1}b, $$ (13) $$ {F_{\rm{B}}} = \frac{1}{2}\left( {{\sigma _{{\rm{B}}1}} + {\sigma _{{\rm{B}}2}}} \right)n{h_1}b =\left[ {\left( {{h_2} + \frac{{{h_1}}}{2}} \right)\rho g - \frac{{2{F_{\rm{S}}}{h_2}}}{{ba}}} \right]n{h_1}b。 $$ (14) 将式(3)、(9)、(13)、(14)的各计算值代入式(8)，化简后得刮板间距上限值的二次方程式：
$$ \begin{split} &\left( {{h_2} + \frac{{{h_1}}}{2}} \right)\left( {\frac{1}{n} - n} \right){h_1}b + \frac{{h_2^3{f_1}}}{a}\left( {1 + {n^2}} \right){h_1} -\\ &\left[ {b\left( {{h_1} + {h_2}} \right) + h_2^2n} \right]{f_1}a = {\text{0}}, \end{split} $$ (15) 令
$$ B = \left( {{h_2} + \dfrac{{{h_1}}}{2}} \right)\left( {\dfrac{1}{n} - n} \right){h_1}b， $$ (16) $$ C = \left( {1 + {n^2}} \right){h_1}h_2^3{f_1}， $$ (17) $$ A = \left[ {b\left( {{h_1} + {h_2}} \right) + h_2^2n} \right]{f_1}， $$ (18) 代入式(15)得：
$$ B + \dfrac{C}{a} - Aa = 0， $$ (19) $$ {a_{\max }} = \frac{{B + \sqrt {{B^2} + 4AC} }}{{2A}}。 $$ (20) 为保证有机肥的稳定排施，a的合理范围应为amin~amax。依据式(7)、(20)得amin=0.004 m，amax=0.166 m。为了获得较大的有机肥承载量，由图4中h2和a的函数关系，刮板间距应该选取尽可能大的值，为了便于加工，a取值为0.15 m。
2.3 排肥器转速
根据刮板输送式排肥装置的排肥原理，由式(1)可以计算出刮板输送式排肥装置排施有机肥的理论单转排肥量。依据“1.1”所述苏南地区草莓大棚施撒基肥习惯用量，假设施肥机械的行驶速度为0.5~1.2 m·s−1，为了满足当地基肥习惯量，刮板输送式排肥器的转速应该为16.5~59.5 r·min−1。
2.4 排肥器开度
排肥器开度的下限值取为刮板的高度，上限值为最大单转排肥量所对应的排肥器开度。单位面积内的排肥量(Q)可由下式计算：
$$ Q = qn\frac{s}{{bv}}, $$ (21) 式中，q为排肥器单转的排肥量，kg·r−1；n为排肥器转速，r·min−1；s为施肥面积，m2；b为施肥幅宽(即刮板长度)，m；v为施肥机具前进速度，m·s−1。
根据“1.1”所述最大基肥习惯用量可得单位面积内的最大施肥量Q为2.25 kg，施肥机具前进速度为1.2 m·s−1；转速n为16.5 r·min−1时，所对应的挡肥板开度为35 mm，则排肥器开度调节为10~35 mm范围内即可满足排肥量的要求。
在确定h1、a、b、z、d、ρ这些参数后，排肥装置的排肥量只与排肥轴转速(n)和开度(H)有关。
3. 排肥器工作性能试验
排肥器开度大小会影响出料高度，进而影响排肥过程中肥料的下落，影响排肥量的稳定性；排肥轴转速影响出料的初始速度和排肥装置的振动，从而影响有机肥的落肥效果，造成排肥量稳定性的变化。排肥器开度、排肥轴转速不仅影响排肥稳定性，也是影响排肥量的重要因素，为探讨排肥器工作性能，优化其设计参数，搭建试验台并进行排肥量和排肥稳定性试验。
3.1 材料与评价指标
3.1.1 试验材料
试验使用的有机肥为长春市新阳光有机肥料有限公司生产的粉末状有机肥，主要以鸡粪为主，其中氮、磷、钾总养分质量分数为9.9%，有机质质量分数为75%。
试验装置如图5所示，采用步进电机(86BYG250H)通过步进电机控制器(KH-01)调节排肥器开度；采用三相异步电机(YS-903-4)连接减速机(型号为NMRV63，减速比为1∶20)，通过变频器(AE200-2-1PH-3.0G)调节排肥轴转速，并用编码器(LPD3806-400BM-G5-24C)实时读取排肥轴转速。
3.1.2 排肥性能评价指标
试验以排肥量稳定性变异系数作为评价标准，计算方法如下：
$$ x = \frac{{\displaystyle\sum\nolimits_{i = 1}^n {{x_i}} }}{n} , $$ (22) $$ S = \sqrt {\frac{{\displaystyle\sum\nolimits_{i = 1}^n {{{\left( {{x_i} - x} \right)}^2}} }}{{n - 1}}} , $$ (23) $$ V = \frac{S}{x} \times 100 , $$ (24) 式中，xi为各次排肥量，kg；x为各次排肥量的平均值，kg；S为标准差，kg；V为排肥稳定性变异系数，%；n为测定次数。
3.2 转速与排肥量相关性试验
排肥轴转速和排肥器开度影响排肥量，为考察各因素对排肥量的影响，基于前文的分析，在转速15~60 r·min−1和开度10~35 mm范围内进行试验，试验时根据排肥器开度与步进电机步距角对应关系，设置开度分别为10、15、20、25、30、35 mm，调节变频器旋钮，设置转速分别为15、24、33、42、51、60 r·min−1，排肥装置在设定参数下运行10 s，在排肥口处用塑料盆接肥，结束后采用精度±5 g的电子秤称量并记录排肥量，每个转速和开度重复3次试验，试验结果见图6。
![]()
图 6 不同开度水平下转速与排肥流量的关系Figure 6. The relationship between rotational speed and fertilizer discharging flow under different opening degrees由图6可得，排肥量与排肥轴转速、排肥器开度密切相关，排肥量随着排肥器开度和排肥轴转速的增大而增加，且同一排肥器开度下排肥轴转速与排肥量呈良好的线性关系(整体决定系数≥0.9857)，不同的转速匹配不同的排肥器开度可实现不同排肥量的需求。因此可以通过实时控制排肥轴转速或同时控制排肥器开度实现变量施肥。
3.3 响应面试验
3.3.1 试验设计
通过分析排肥轴转速和排肥器开度对刮板输送式排肥装置的排肥效果，对排肥装置的结果参数和运行参数进行优化。选用二次通用旋转组合试验方法，分析排肥轴转速和排肥器开度对排肥稳定性的影响，以排肥稳定性变异系数为评价指标，试验因素及水平编码见表1。
表 1 试验因素及水平编码表Table 1. Coding of test factors and levels编码值
Code value因素 Factor 开度/mm
Opening (x1)转速/(r·min−1)
Rotational speed (x2)γ(1.414) 35.0 59.5 1 31.3 53.0 0 22.5 38.0 −1 13.7 23.0 −γ(−1.414) 10.0 16.5 3.3.2 试验方案与结果分析
应用Design-Expert软件对试验结果进行回归分析，以确定2个试验因素下排肥稳定性变异系数的变化规律，试验方案及结果如表2所示。以排肥器开度和排肥轴转速为试验因素，以排肥稳定性变异系数为试验指标，运用Design-Expert软件对试验结果进行方差分析，方差分析结果见表3。
表 2 正交试验方案及结果1)Table 2. Orthogonal test plan and results试验编号
Test
number试验因素
Experimental factor排肥稳定性变异系数/%
Coefficient of variation of
fertilization stabilityx1 x2 1 22.5 38.0 2.47 2 31.3 53.0 1.89 3 35.0 38.0 1.85 4 22.5 38.0 2.53 5 31.3 23.0 3.21 6 22.5 16.5 3.38 7 22.5 59.5 2.47 8 13.7 23.0 2.62 9 22.5 38.0 2.31 10 13.7 53.0 2.59 11 22.5 38.0 2.45 12 10.0 38.0 2.43 13 22.5 38.0 2.31 1) x1、x2分别表示排肥器开度和排肥轴转速
1)x1、x2 indicated opening degree of fertilizer discharging device and rotational speed of the fertilizer discharging shaft表 3 方差分析表1)Table 3. Variation analysis table变异来源
Variance
source平方和
Sum of
squares自由度
Degree of
freedomF P 模型 Model 2.080 5 25.70 0.0002 x1 0.110 1 6.71 0.0359 x2 0.960 1 59.10 0.0001 x1x2 0.420 1 25.73 0.0014 x12 0.098 1 6.05 0.0434 x22 0.430 1 26.68 0.0013 残差 Residual 0.110 1 失拟项 Lack of fit 0.074 3 2.48 0.2000 纯误差 Pure error 0.040 4 总和 Sum 2.190 12 1) x1、x2分别表示排肥器开度和排肥轴转速
1)x1、x2 indicated opening degree of fertilizer discharging device and rotational speed of the fertilizer discharging shaft由表3可知，模型的显著性检验F=25.70，P<0.01，决定系数R2=0.9483，表示回归模型极显著，且在试验范围内拟合程度好。由变异系数模型可知，x2、x1x2、x22对方程影响极显著(P<0.01)，x1、x12对方程影响显著(P<0.05)，得到变异系数的回归方程为：
$$ \begin{split} V =& 3.41 - 0.12{x_1} - 0.34{x_2} - 0.32{x_1}{x_2} -\\ &0.12x_1^2 + 0.24x_2^2 {\text{。}} \end{split} $$ (25) 3.3.3 试验因素对指标的影响
排肥器开度和排肥轴转速对排肥稳定性变异系数的影响如图7所示，由图7a可知，在排肥器开度较大时变异系数随着排肥轴转速减小而增大，在排肥器开度较小时变异系数随着排肥轴转速先增大后减小。由图7b可知排肥轴转速对排肥稳定性变异系数的影响比排肥器开度大，该结论与方差分析得出的结论相吻合。
![]()
图 7 试验因素对变异系数的影响Figure 7. Influence of experimental factor on the coefficient of variation3.3.4 试验结果目标优化与试验验证
为寻求各因素最优组合，以排肥稳定性变异系数最小为原则，应用Design-Expert软件对排肥稳定性系数优化求解，得排肥器开度30.42 mm，排肥轴转速51.5 r·min−1，预测变异系数最优值为1.84%。在自制试验台上进行排肥稳定性变异系数和排肥量试验，重复5次，试验结果如表4所示。变异系数的相对误差最大值为4.89%，最小值为0.54%，平均相对误差仅为2.95%，相对误差均在5%范围内，为可接受的误差范围，满足实际试验要求。
表 4 最优条件验证试验结果Table 4. Verification test result of optimal condition试验编号
Test number排肥量/g
Fertilizer discharging volume排肥稳定性变异系数/%
Coefficient of variation of
fertilization stability变异系数相对误差/%
Relative error of coefficient
of variation1 2 3 均值
Mean1 60.22 61.56 60.59 61.46 1.93 4.89 2 60.36 62.64 61.15 61.45 1.89 2.71 3 62.34 60.15 61.49 61.33 1.80 2.17 4 61.06 62.40 63.23 62.23 1.76 4.43 5 61.54 63.75 62.15 62.48 1.83 0.54 均值 Mean 61.10 62.10 61.72 61.79 1.84 2.95 3.3.5 排肥稳定性对比验证试验
为验证排肥器排肥稳定性，与目前果园施撒有机肥应用最多的螺旋式排肥器[28-30]进行试验对比。果园内施肥作业空间有限，基肥的施撒量也同样巨大[31]。试验时设定两者的工作转速都是51.5 r·min−1，试验时间10 s，进行5组试验，每组试验重复3次取平均值，统计排肥量和试验的稳定性变异系数，试验结果如表5所示。本研究设计的刮板输送式排肥器排肥稳定性变异系数的均值为1.84%，而螺旋式排肥器排肥稳定性变异系数达6.12%，相较于螺旋式排肥器有效提高了有机肥的排肥稳定性，满足设计要求。
表 5 排肥稳定性对比验证试验结果Table 5. Comparative verification test results of fertilization stability试验编号
Test number排肥量/g
Fertilizer discharging volume变异系数/%
Coefficient of variation1 2 3 均值
Mean1 8.32 8.74 7.95 8.34 4.74 2 8.65 8.15 7.94 8.25 4.42 3 7.62 8.20 9.05 8.29 8.71 4 7.65 7.94 8.46 8.02 5.09 5 7.72 8.95 8.11 8.26 7.62 均值 Mean 7.99 8.40 8.30 8.23 6.12 4. 结论
1)设计了一种刮板输送式有机肥排肥装置，通过受力分析确定排肥装置的相关参数，保证排肥装置具有较好的排肥性能。
2)试验研究表明，在试验设定的转速和开度范围内，不同开度条件下，转速和排肥量均呈良好的线性关系，整体的决定系数大于0.9857，表明可通过实时调节开度和转速实现排肥量的调节。
3)通过排肥稳定性二次通用旋转组合试验的二次回归分析和响应面分析，表明排肥轴转速对排肥稳定性变异系数的影响大于排肥器开度；在排肥器开度为30.42 mm、排肥轴转速为51.5 r·min−1时排肥稳定性变异系数最优，为1.84%，验证试验与优化结果基本吻合，对比目前常用的螺旋式有机肥排肥器，设计的刮板输送式有机肥排肥器有效地提高了有机肥的排肥稳定性。
4)初步试验表明，该排肥装置设计合理，满足大棚内有机肥施撒的农艺要求，能够实现有机肥的定量准确施撒。
-
![]()
图 1 优势种序位曲线
米槠、木荷、栲、黄樟、红锥和华润楠：在A图中的序位依次为1、2、4、8、9和11；在B图中的序位依次为1、2、9、20、24和19；在C图中的序位依次为2、1、5、8、21和6；在D图中的序位依次为1、2、3、8、5和10。
Figure 1. The sequence curves of dominant species
表 2 6个树种木材密度在径阶间的差异1)
Table 2 Differences in wood density among diameter grades for six tree species
下载: 导出CSV
表 3 6个树种木材密度在坡向间的差异性1)
Table 3 Differences in wood density among slope aspects for six tree species
下载: 导出CSV
-
[1] BERGMAN R, CAI Z, CARLL C G. Wood handbook: Wood as an engineering material[M]. Madison: United States Department of Agriculture, Forest Service, Forest Products Laboratory, 2010: 1-2.
[2] WU R. Microstructural study of sanded and polished wood by replication[J]. Wood Sci Technol, 1998, 32(4): 247-260. doi: 10.1007/BF00702893
[3] GJERDRUM P, EIKENES B. A model for spatial wood density gradients in Norway spruce stems and stochastic between-stem dissimilarities for basic and dry density[J]. Wood Sci Technol, 2014, 48(1): 71-84. doi: 10.1007/s00226-013-0586-z
[4] MACHADO J S, LOUZADA J L, SANTOS A J A, et al. Variation of wood density and mechanical properties of blackwood (Acacia melanoxylon R. Br.)[J]. Mater Design, 2014, 56: 975-980. doi: 10.1016/j.matdes.2013.12.016
[5] MOLTEBERG D, HOIBO O. Modelling of wood density and fibre dimensions in mature Norway spruce[J]. Can J Forest Res, 2007, 37(8): 1373-1389. doi: 10.1139/X06-296
[6] POORTER L, WRIGHT S J, PAZ H, et al. Are functional traits good predictors of demographic rates? Evidence from five neotropical forests[J]. Ecology, 2008, 89(7): 1908-1920. doi: 10.1890/07-0207.1
[7] 罗云建. 华北落叶松人工林生物量碳计量参数研究[D]. 北京: 中国林业科学研究院, 2007: 52-53. http://cdmd.cnki.com.cn/Article/CDMD-82201-2007173541.htm [8] 中国林业科学研究院木材工业研究所.中国主要树种的木材物理力学性质[M].北京:中国林业出版社, 1982. [9] 侯振宏. 中国林业活动碳源汇及其潜力研究[D]. 北京: 中国林业科学研究院, 2010: 78-79. http://cdmd.cnki.com.cn/Article/CDMD-82201-2010264599.htm [10] 朱江玲, 石岳, 方乐祺, 等.中国主要树种木材物理力学属性的地理格局及其环境控制[J].中国科学(生命科学), 2015, 45(1): 56-67. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=JCXK201501006&dbname=CJFD&dbcode=CJFQ [11] 徐金梅, 吕建雄, 鲍甫成, 等.祁连山青海云杉木材密度对气候变化的响应[J].北京林业大学学报, 2011, 33(5): 115-121. http://d.old.wanfangdata.com.cn/Periodical/bjlydxxb201105021 [12] 刘福涛. 立地对马占相思人工林生长和材性的影响[D]. 福州: 福建农林大学, 2007: 25-29. http://cdmd.cnki.com.cn/article/cdmd-10389-2007135795.htm [13] PATINO S, LLOYD J, PAIVA R, et al. Branch xylem density variations across the Amazon Basin[J]. Biogeosciences, 2009, 6(4): 545-568. doi: 10.5194/bg-6-545-2009
[14] WIEMANN M C, WILLIAMSON G B. Geographic variation in wood specific gravity: Effects of latitude, temperature, and precipitation[J]. Wood Fiber Sci, 2002, 34(1): 96-107. http://cn.bing.com/academic/profile?id=521c8238fc0c20590a2546f032524c68&encoded=0&v=paper_preview&mkt=zh-cn
[15] 鲍甫成, 江泽慧.中国主要人工林树种木材性质[M].北京:中国林业出版社, 1998. [16] 成俊卿.木材学[M].北京:中国林业出版社, 1985. [17] 徐明锋, 胡砚秋, 李文斌, 等.土壤养分对亚热带天然林物种分布的影响[J].中南林业科技大学学报, 2014, 34(9): 91-97. doi: 10.3969/j.issn.1673-923X.2014.09.018 [18] 潘秋荣, 梁佰华, 陈小芸.广东东源康禾省级自然保护区野生药用植物资源调查[J].广东林业科技, 2012, 28(4): 25-30. doi: 10.3969/j.issn.1006-4427.2012.04.006 [19] 林中大.广东康禾自然保护区植物区系研究[J].广西林业科学, 2007, 36(2): 75-77. doi: 10.3969/j.issn.1006-1126.2007.02.003 [20] 胡传双, 苏志尧, 云虹, 等.流体静力称衡法测定生长锥取样木材密度的改进方法[J].华南农业大学学报, 2010, 31(3): 105-108. doi: 10.3969/j.issn.1001-411X.2010.03.026 [21] 蒋燚, 李志辉, 朱积余, 等.红锥家系木材密度等物理性状的遗传及变异性分析[J].中南林业科技大学学报, 2012, 32(11): 9-13. http://d.old.wanfangdata.com.cn/Periodical/znlxyxb201211002 [22] 刘青华, 张蕊, 金国庆, 等.马尾松年轮宽度和木材基本密度的种源变异及早期选择[J].林业科学, 2010, 46(5): 49-54. http://d.old.wanfangdata.com.cn/Periodical/lykx201005008 [23] NGUYEN H, FIRN J, LAMB D, et al. Wood density: A tool to find complementary species for the design of mixed species plantations[J]. Forest Ecol Manag, 2014, 334: 106-113. doi: 10.1016/j.foreco.2014.08.022
[24] 王秀花, 陈柳英, 马丽珍, 等. 7年生木荷生长和木材基本密度地理遗传变异及种源选择[J].林业科学研究, 2011, 24(3): 307-313. http://d.old.wanfangdata.com.cn/Periodical/lykxyj201103005
计量
- 文章访问数:
- HTML全文浏览量: 0
- PDF下载量:
