Site quality evaluation of broadleaved mixed forest in Xiaoxing’an Mountains of northeastern China based on volume potential productivity
-
摘要:目的
立地质量评价是森林经营规划和决策的重要基础。建立林分蓄积生长模型,引入林分蓄积潜在生产力,对小兴安岭阔叶混交林进行立地质量评价,旨在为该地区森林生长收获预测、森林经营规划以及可持续经营提供科学依据。
方法利用黑龙江省2010—2015年两期连续清查中的小兴安岭阔叶混交林数据,建立林分平均高和蓄积生长模型,划分立地等级,并计算各立地等级阔叶混交林林分蓄积潜在生产力。结合林分现实生产力及其与蓄积潜在生产力之间的差异,对当前林龄林分的实际状况和可提升空间进行评价。
结果(1)小兴安岭阔叶混交林林分蓄积生长模型拟合效果较好,林分基准年龄(30 a)时,各立地等级的蓄积潜在生产力范围为3.87 ~ 6.10 m3/(hm2·a);阔叶混交林林分平均可提升空间介于50.78% ~ 62.97%之间,且随着立地等级的降低而逐渐增大。(2)小兴安岭阔叶混交林林分蓄积潜在生产力随林龄的增大而逐渐减小,最终趋于平缓。其与现实生产力的差异随林龄的增大而逐渐减小。(3)各林龄林分蓄积潜在生产力和现实生产力均呈现随立地等级降低而下降的趋势,立地质量中等的林分占总样本量的77.24%,且各个立地等级的林分蓄积潜在生产力始终大于现实生产力。
结论林分蓄积潜在生产力可以作为小兴安岭阔叶混交林立地质量评价的指标。针对小兴安岭阔叶混交林,各立地等级处于发育初期的林分蓄积潜在生产力比近成熟林的更大,且林分可提升空间也更大;随着立地质量的降低,林分可提升空间逐步扩大。因此,应根据林分发育时期和立地条件,采取差异化抚育措施,且需注重立地质量中等的林分,优先进行中幼龄林抚育,提升林分的整体生产力。
Abstract:ObjectiveSite quality evaluation is an important foundation for forest management planning and decision-making. This article introduces the potential productivity of stand volume by establishing a stand volume growth model to evaluate the site quality of broadleaved mixed forest in the Xiaoxing’an Mountains of northeastern China, providing a scientific basis for predicting forest growth and harvest, forest management planning, and sustainable management in the region.
MethodUsing data from the two consecutive inventory periods of the Xiaoxing’an Mountains broadleaved mixed forest in Heilongjiang Province from 2010 to 2015, a model for average tree height and volume growth of the forest stand was established. Site levels were classified, and potential productivity of each site level for broadleaved mixed forest stands was calculated. We evaluated the actual condition and potential improvement space of the current aged stand based on the difference between actual productivity of the forest and its potential productivity for accumulation.
Result(1) The fitting effect of volume growth model of broadleaved mixed forest in the Xiaoxing’an Mountains was good, and the range of potential productivity at the baseline age of forest was between 3.87−6.10 m3/(ha·year). The average improvement space of broadleaved mixed forest ranged from 50.78% to 62.97%, and gradually increased with decrease of site level. (2) The potential productivity of broadleaved mixed forest in Xiaoxing’an Mountains gradually decreased with increase of stand age, and eventually tended to flatten out. The difference between realized productivity and the potential productivity gradually decreased with increase of stand age. (3) The potential productivity and realized productivity of each age group showed a decreasing trend with decrease of site grade. Forests with moderate site quality accounted for 77.24% of total sample size, and potential productivity of each site grade was always greater than the realized productivity.
ConclusionThe potential productivity of stand volume can serve as an indicator for site quality evaluation of broadleaved mixed forest in the Xiaoxing’an Mountains. For the broadleaved mixed forest in Xiaoxing’an Mountains, stands at the initial stage of development within various site classes exhibit greater potential productivity of stand volume compared with those nearing maturity, and the potential for improvement in stand productivity is also larger. As site quality decreases, the scope for stand improvement gradually expands. Therefore, differentiated management measures should be adopted based on stand development stage and site conditions, with emphasis on stands of medium site quality. Priority should be given to tending of middle-aged and young forests to enhance overall stand productivity.
-
我国水果生产和消费能力均排名世界第一,但多数果园种植在丘陵山地,果树修整枝、疏花疏果、果实套袋和采收等作业主要依靠人工架梯来完成,劳动强度大、作业效率低、安全风险高[1-3],采摘器等辅助人工摘果装备同样存在作业效率低的问题[4-5]。丘陵山地地形地貌变化大,普通车辆很难进入,为了克服地形的影响,农业车辆多采用调平的方式来应对[6-9]。因此,设计研发一款可调平的果园升降作业平台对于提高果园机械化水平和作业效率,降低工作风险具有重要意义。
国外果园机械研究起步较早,北美自20世纪50年代开始把高空作业平台应用到果园中,发展至今已具有较高水平[10];意大利N.P.SEYMOUR公司生产的Windegger Picking Platform系列可全轮转向,工作台两级升降、角度可调,适用性较广[11];韩国SUNGBOO公司生产的SB系列采用纯电力能源,工作台可调角度10°,可满足缓丘陵作业需求[12-13]。
国内对于果园升降机械的研究起步较晚,主要因为果园种植模式以家庭分散种植为主,规模小,不利于果园机械的发展[14]。近些年,果园的主要种植模式为矮砧和密植,为果园机械化发展提供了前提条件。刘西宁等[15]研发了国内第一台果园升降平台,但不具备调平功能,不适合在丘陵山地作业。崔志超等[16]针对温室果蔬机械化低,传统燃油作业平台污染大等问题,设计制造了一款小型电动履带作业平台,采用电能作为动力,绿色环保,使用远程和在线的双操作系统。樊桂菊等[17-18]以拖拉机为载体加装工作台来实现调平,通过融合卡尔曼滤波的模糊PID控制系统实现自动调平,最大载荷150 kg,举升高度1.5 m,但是工作平台较小,在大规模果园作业时效率偏低。刘大为等[19-20]以南方丘陵山区柑橘园为对象研制了小型果园升降平台,通过“二次调平”方法,由液压缸伸缩进行调平,使用履带底盘,转向半径小、稳定性高、结构紧凑、操作简单,但调平方式为人工手动调平,调平精度低,而且作业效率不高。邱威等[21]根据中国南方果园多丘陵山区及果园作业的特点研制了一款折臂升降可调平平台,通过理论计算,仿真分析及样机试验分析验证了在丘陵山区的倾翻稳定性,并且利用正交试验分析了极限静态侧翻角的影响因子权重占比,在作业操作及相关机械的设计中有一定的借鉴意义。
国内公司现有果园产品多,但均不具备调平功能,无法适应丘陵山地果园作业,形成了“有机不好用”的局面,丘陵果园综合机械化水平仅为5.75%[2]。因此,针对丘陵山地矮砧密植苹果园地形与种植模式特点,研究升降平台调平机理,设计一款可调平的果园自走式升降作业平台,以期满足丘陵山地作业需求。
1. 整机结构与原理
1.1 设计要求
丘陵果园因地形特点,其种植模式主要有依山就势的缓坡模式和“坡改梯”的梯田模式两种。在这类果园中作业和转场对农业车辆的动力底盘要求很高,应具备良好的通过性和爬坡能力,本机选用越障能力强、转弯半径小的履带底盘。经查阅文献[22]和实地观测,苹果园一般株距1.0 ~ 1.5 m,行距3.5 ~ 4.0 m,树高不超过3.5 m。结合实际生产过程,可调平平台具体设计要求如表1所示。
表 1 果园作业平台设计参数Table 1. Design parameters of orchard operating platform参数 Parameter 数值 Value 整机尺寸 Overall size 1 900 mm × 1 400 mm × 2 000 mm 整机质量 Overall mass/kg 1 200 最大载质量 Maximum load/kg 300 升降高度 Lifting height/mm 0 ~ 1 150 转弯半径 Turning radius/m 3 行驶速度 Drving speed/(km·h−1) 0 ~ 5 1.2 整机结构
整机结构主要包括护栏、工作台、俯仰调平机构、机架、侧倾调平机构、履带底盘、操作台、剪叉机构和动力集成(图1)。动力集成安装在机架上,由电机、电池、液压管道和液压阀组成,为平台的升降、调平和行走等动作提供动力,并通过操作台控制平台升降和行走。工作台两侧可扩展,以增加工作平台面积,方便果农采摘套袋等作业。以下对平台的升降和调平机理做详细阐述。
1.3 工作台升降及俯仰调平原理
1.3.1 整体结构设计
因剪叉式升降机构具有较高的稳定性和承载能力[23-24],同时考虑到缓丘陵果园地形特点(坡度 ≤ 15°),以及矮砧密植型苹果树高一般不超过3.5 m,结合GB10000—1988《中国成年人人体尺寸标准》[25]可知平台升降高度达2 m即可满足果园作业要求,故选择单级剪叉升降机构,并在此基础上改进为折臂剪叉式俯仰调平机构(图2)。升降液压缸PQ伸长,工作台随之上升,反之工作台下降;俯仰液压缸MN伸长,俯仰臂O1D绕O1点旋转,工作台DE仰起一定角度θ(如图2中红色部分),反之工作台俯下一定角度−θ,达到俯仰调平目的。该调平方式没有额外增加过多机构,节省空间、简单可靠。
图 2 升降及俯仰调平机构原理O、O1分别为升降、俯仰调平机构坐标系原点;DE为工作台,PQ为升降液压缸,MN为俯仰液压缸;OD和EH为剪叉臂,其中OB、BD、EB和BH长度相等(mm);A、C两点为升降液压缸上下两端与剪叉臂的连接点,且AP与EH垂直,CQ与OB垂直;O和E分别为剪叉臂与机架、工作台的铰接点;D、H两点为铰接滑轮,滑轮与机架、工作台的滑轨相切,B点是两剪叉臂铰接点;h为工作台升起高度(mm);G为工作台和被载物体的总重(N);α为剪叉臂与水平面夹角(°);β为MO1与O1B的夹角(°);φ为俯仰臂O1D绕O1点转动的角度(°);θ为工作台与水平面夹角,即平台俯仰角(°);F1为升降液压缸推力(N);F2为俯仰液压缸推力(N)。O is the origin of lifting coordinate system, and O1 is the origin of pitch leveling mechanism coordinate system. DE is the platform, PQ is lifting hydraulic cylinder, MN is pitch hydraulic cylinder; OD and EH are scissor arms, and the lengths of OB, BD, EB and BH are equal (mm). A is the connection point between the upper end of lifting hydraulic cylinder and scissor arm, C is connection point between the lower end of lifting hydraulic cylinder and scissor arm, and AP is perpendicular to EH, and CQ is perpendicular to OB. O is hinge point between scissor arm and frame, E is the hinge point of scissor arm and platform; D and H are hinged pulleys, the pulley is tangent to slide rail of the rack and platform, B is hinge point of two scissor arms. h is lifting height of platform (mm). G is total mass of platform and the loaded object (N). α is the angle between scissor arm and horizontal plane (°). β is the angle between MO1 and O1B (°). φ is the angle that tilt arm O1D rotates around O1 point (°). θ is platform pitch angle between platform and horizontal plane (°). F1 is the thrust of lifting hydraulic cylinder (N). F2 is thrust of pitch hydraulic cylinder (N).Figure 2. Principles of lifting and pitch leveling mechanism1.3.2 升降和俯仰液压缸推力分析
根据图2建立液压缸推力、俯仰液压缸位移量和平台角度关系式,以确定优化液压缸所需推力和平台尺寸。由图2可知相关量的几何关系:
{xPQ=−lCQsinα+lBCcosα−lABcosα+lAPsinαyPQ=lCQcosα+lBCsinα+lABsinα+lAPcosαlPQ=√xPQ2+yPQ2h=lEHsinα (1) 式中:xPQ、yPQ分别为升降液压缸PQ在OX、OY轴上的投影长度(mm);lCQ、lBC、lAB和lAP分别为升降机构部分CQ、BC、AB和AP的长度(mm);lPQ为升降液压缸PQ的长度(mm);lEH为剪叉臂EH的长度(mm);α为剪叉臂与水平面夹角(°)。
根据虚位移原理[26]:
F1dlPQdα=Gdhdα (2) 得到:
F1=G2lEHcosα√asin2α+bsin2α+ccos2α(a−c)sin2α+2bcos2α (3) 式中:F1为升降液压缸推力(N);
a=(lAB+lBC)2+(lAP−lCQ)2 ,b=2(lABlCQ+lBClAP) ,c=(lBC−lAB)2+(lAP+lCQ)2 。由式(3)可知:F1不仅与lEH、lPQ、lBC、lAB和lAP有关,还与平台升降高度α有关,并且随着α的减小,液压缸所需推力随之增大。在工作台降到最低位置时,液压缸所需推力最大。
同理,在图2中以X1O1Y1为坐标系,根据虚位移原理可知:
F2=GlO1Dcos(α+φ)√lO1M2+lO1N2+2lO1MlO1Ncos(φ−β)2lO1MlO1Nsin(β−φ) (4) 式中:F2为俯仰液压缸推力(N);
lO1D 、lO1M 和lO1N 分别为俯仰调平机构部分O1D、O1M、O1N的长度(mm);β为MO1与O1B的夹角(°);φ 为俯仰臂O1D绕O1点转动的角度(°)。由式(4)可知:俯仰液压缸所需推力和平台俯仰角及平台升降高度有关,平台高度越低,俯仰角越大,所需推力越大,在平台降到最低且俯仰液压缸完全伸出时所需推力最大。考虑到工作台不能与底部动力集成干涉以及对升降高度的要求,对参数值选取合适的值。详细尺寸数据见表2。
表 2 升降及俯仰调平机构尺寸Table 2. Dimensions of lifting and pitch leveling structures参数 Parameter 数值 Value lAB/mm 220 lBC/mm 702 lAP/mm 70 lCQ/mm 97 lEH/mm 1 722 lO1N/mm 411 lO1M/mm 163 lO1D/mm 561 β/(°) 88 注:lAB、lBC、lAP和lCQ分别为升降机构AB、BC、AP和CQ的长度;lEH为剪叉臂EH的长度;lO1N、lO1M和lO1D分别为俯仰调平机构部分O1N、O1M和O1D的长度;β为MO1与O1B的夹角。Notes:lAB, lBC, lAP and lCQ are length of lifting mechanism AB, BC, AP and CQ, respectively. lEH is the length of scissor arm EH. lO1N, lO1M and lO1D are length of pitch leveling mechanism O1N, O1M and O1D, respectively. β is angle between MO1 and O1B. 1.3.3 平台俯仰调平倾角
在XOY坐标中:
{xOD1=lOO1cosα+lO1D1cos(α+φ)yOD1=lOO1sinα+lO1D1sin(α+φ) (5) 式中:
xOD1 、yOD1 分别为剪叉臂OD1在OY、OX轴上的投影长度(mm),lOO1 和lO1D1 分别为剪叉臂OO1和O1D1的长度(mm)。因为
tanθ=yOD1−hxOD1 ,可以得到:θ=arctan(lOO1−lEH)sinα+lO1D1sin(α+φ)lOO1cosα+lO1D1cos(α+φ) (6) 其中:
φ=arccoslO1M2+lO1N2−(lMN−y)22lO1MlO1N (7) 式中:y为俯仰液压缸活塞杆位移量(mm),
lMN 为俯仰液压缸MN的长度(mm),θ为工作台与水平面夹角,即平台俯仰角(°)。由式(6)和式(7)可知折臂剪叉俯仰调平机构可调角度与升降液压缸、俯仰液压缸的位移量有关。当平台降到最低点位置时,俯仰液压缸不能缩短,只能伸长。即平台只能进行仰调平,不能进行俯调平,否则工作台会与机架干涉。随着工作台的升起,工作台既能仰也能俯,可调角度均为10°。但工作台升到最高点时仰调平角度降至6.5°。
1.4 侧倾调平原理
1.4.1 侧倾调平机构设计
在履带底盘和机架之间设计了侧倾调平机构,起到连接和调平作用,机构简图如图3。液压缸I1J1和IJ做相反运动,机架O3J1J绕O3点旋转,工作台侧倾一定角度θ1,达到侧倾调平目的。
图 3 侧倾调平机构原理I1J1、IJ为侧倾液压缸;α1为侧倾液压缸I1J1与水平面的夹角(°);β1为侧倾液压缸IJ与水平面的夹角(°);θ1为平台与水平面的夹角,即平台侧倾角(°);O2为原点,I、I1、O3为底盘上的固定点;O3J1J为机架;O3R1与侧倾液压缸I1J1垂直,O3R与侧倾液压缸IJ垂直;e为机架上方总重G1在X2轴平行方向上与O2的距离,(mm);F3、F4为侧倾液压缸的推力(N)。I1J1 and IJ are the roll hydraulic cylinders, α1 is angle between roll hydraulic cylinder I1J1 and the horizontal plane (°). β1 is angle between roll hydraulic cylinder IJ and horizontal plane (°). θ1 is the angle between platform and horizontal plane, i.e. the platform roll angle (°). O2 is the origin, I,I1 and O3 are the fixed points on chassis, O3J1J is frame. O3R1 is perpendicular to the roll cylinder I1J1, O3R is perpendicular to roll cylinder IJ. e is the distance between total mass G1 and O2 in the direction parallel to X2 axis. F3 and F4 are thrust of roll hydraulic cylinder (N).Figure 3. Principle of roll leveling mechanism1.4.2 侧倾液压缸推力分析
由图3可知:
{∠IO3R=β1−∠IO3O2lO3R=lIO3cos∠IO3R∠I1O3R1=α1−∠I1O3O2lO3R1=lI1O3cos∠I1O3R1F3lO3R+F4lO3R1=G1e (8) 式中:
lO3R 、lIO3 、lO3R1 和lI1O3 分别为侧倾调平机构部分O3R、IO3、O3R1和I1O3的长度(mm);F3、F4为侧倾液压缸的推力(N);∠IO3R为IO3与O3R的夹角(°);∠IO3O2为IO3与O3O2的夹角(°);∠I1O3R1为I1O3与O3R1的夹角(°);∠I1O3O2为I1O3与O3O2的夹角(°);G1为机架上方总重(N);e为G1在X2轴平行方向上与O2的距离(mm)。根据选用的启创液压科技公司HSG-40系列工程油缸得1.45F3 = F4,则:
F4=G10.6elIO3cos(β1−∠IO3O2)+lI1O3cos(α1−∠I1O3O2) (9) 从式(9)中可以看出G1的增加及因G1升高引起的e增加会导致液压缸推力的增加,在设计中对工作台零部件轻量化,降低重心,以减小液压缸所需推力,同时也能增加平台的侧倾稳定性。
1.4.3 平台侧倾调平倾角
由图3中关系可知:
θ1=arctanlIJsinβ1−lI1J1sinα1lI1J1cosα1+lI1I+lIJcosβ1 (10) 式中:
lIJ 、lI1J1 和lI1I 分别为IJ、I1J1和I1I的长度(mm);α1为侧倾液压缸I1J1与水平面的夹角(°);β1为侧倾液压缸IJ与水平面的夹角(°);θ1为平台与水平面的夹角,即平台侧倾角(°)从式(10)中可以看出侧倾液压缸的长度和安装位置是影响平台调平角度的主要因素。结合设计要求与结构特点选取合适的尺寸,侧倾角度达 ± 8°。具体尺寸数据见表3。
表 3 侧倾调平结构尺寸Table 3. Dimensions of roll leveling structures参数 Parameter 数值 Value α1/(°) 28 ~ 70.5 lIO3/mm 264 lIJ/mm 320 ~ 390 lI1I/mm 373 注:α1为侧倾液压缸I1J1与水平面的夹角;lIO3为侧倾调平机构部分IO3的长度;lIJ和lI1I分别为IJ和I1I的长度。Notes: α1 is the angle between roll hydraulic cylinder I1J1 and horizontal plane. lIO3 is length of tilting leveling mechanism IO3. lIJ and lI1I are the lengths of IJ and I1I, respectively. 2. 控制系统仿真
2.1 系统数学模型建立
升降、俯仰和侧倾3个控制回路组成整个系统,且具有相似的结构形式。因升降部分需手动控制,故只分析俯仰、侧倾控制系统。控制系统由控制器、放大器、电液比例阀、液压缸、倾角传感器等组成。
2.1.1 放大器传递函数
电路部分看作时间常数很小的放大器,可以作为比例环节[27],所以放大器传递函数为
G1(s)=I(s)Θ(s)=Ke1 (11) G2(s)=I(s)Θ(s)=Ke2 (12) 式中:I为电流(A),Θ为平台角度与目标角度差值(°),Ke1、Ke2分别为俯仰、侧倾部分放大器比例系数,G1(s)和G2(s)分别为俯仰、侧倾调平机构的放大器传递函数,s代表时间函数经拉普拉斯变换为复变函数(下同)。
2.1.2 电液比例阀数学模型
选用北京尚泰合力液压科技公司生产固有角频率69 rad/s的STPLS系列电液比例阀。建立电液比例阀的传递函数时,如果系统中存在固有角频率较大的动力元件选用二阶震荡环节[28-29]。所以电液比例阀的传递函数为:
G3(s)=Xv(s)I(s)=Kvs2ωv2+2ξvωvs+1 (13) 式中:ωv为比例阀固有角频率(rad/s),
ξv 为比例阀阻尼比,Xv为比例阀芯位移(m),Kv比例阀增益。2.1.3 液压缸模型建立
选用启创液压科技公司的HSG-40系列工程油缸,属于阀控非对称缸,其传递函数[30]为:
G4(s)=Xp(s)Xv(s)=KqAps(s2ωh2+2ξhωhs+1) (14) 式中:液压固有角频率
ωh=√4βeAp2Vtm (rad/s);液压阻尼比ξh=KceAp√βemVt ;m为活塞及其负载总质量(kg);Kq为流量增益系数;Ap为无杆腔活塞有效面积(m2);Kce为总流量压力系数;Vt为阀腔、液压缸和管道总容积(m3);βe 为有效体积弹性模量(Pa);Xp为液压缸活塞位移量(m)。2.1.4 调平机构数学模型建立
根据式(6)、(7)、(10)可知俯仰、侧倾调平机构中液压缸位移量和平台角度可以看做线性关系。所以调平机构的传递函数为:
G5(s)=θ(s)Xp(s)=Kl1 (15) G6(s)=θ(s)Xp(s)=Kl2 (16) 式中:G5(s)、G6(s)分别为俯仰、侧倾调平机构的传递函数,Kl1、Kl2分别为俯仰、侧倾部分液压缸位移角度系数。
2.2 系统稳定性分析
由经典控制理论[31]可知,要使系统输出稳定,幅值裕度和相位裕度须满足如下条件:
{Kg⩾ (17) 将确定的液压系统元件参数带入俯仰、侧倾系统开环传递函数,令开环增益K = 1,绘制其Bode图。此时,俯仰系统幅值裕度Kg = 39.3 dB,相位裕度γ = 88.9°;侧倾系统Kg = 38.1 dB,γ =88.8°。可见系统不满足稳定性要求。利用公式(18)对开环增益校正,取Kg = 10 dB,则俯仰系统K = 29.2,γ = 54.6°,侧倾系统K = 25.3,γ = 58.4°,此时系统满足稳定性要求。可反推放大器比例系数Ke1、Ke2以设计电路部分。
20\lg \! K = {K_{{\rm{g}}1}} - {K_{\rm{g}}} (18) 式中:Kg1为原系统幅值裕度,Kg为系统期望幅值裕度。
则俯仰系统的开环传递函数:
\begin{aligned} & {G_f}(s) = \frac{{\theta (s)}}{{\Theta (s)}} = \\ & \frac{{29.2}}{{8.4 \! \times \! {{10}^{ - 10}}{s^5} \! + \! 9.68 \times {{10}^{ - 8}}{s^4} \! + \! 2.16 \! \times \! {{10}^{ - 4}}{s^3} \! + \! 2 \! \times \! {{10}^{ - 2}}{s^2} \! + \! s}} \end{aligned} (19) 侧倾系统的开环传递函数:
\begin{aligned} & {G_{\rm{c}}}(s) = \frac{{\theta (s)}}{{\Theta (s)}} =\\ & \frac{{25.3}}{{5.25 \! \times \! {{10}^{ - 9}}{s^5} \! + \! 6.17 \! \times \! {{10}^{ - 7}}{s^4} \! + \! 2.46 \! \times \! {{10}^{ - 4}}{s^3} \! + \! 2.06 \! \times \! {{10}^{ - 2}}{s^2} \! + \! s}} \end{aligned} (20) 2.3 仿真分析
调平系统的控制流程是控制器根据反馈的工作台角度差值计算输出电流控制电液比例阀的开口开度,进而控制液压缸伸缩运动,使平台角度始终维持在0°附近。同时,安装在工作台上的倾角传感器把实时角度发送给控制器,形成闭环反馈控制。控制框图如图4所示。
采用增量式PID控制器,优点是当计算机故障时,系统仍能保持原值,平台不会发生危险情况[32]。计算公式为:
\begin{aligned} u(k) = \; & u(k - 1) + {K_{\rm{p}}}(e(k) - e(k - 1)) + {K_{\rm{i}}}e(k) +\\ & {K_{\rm{d}}}(e(k) - 2e(k - 1) + e(k - 2)) \end{aligned} (21) 式中:u(k)为k时刻控制器输出的电流大小,e(k)为k时刻的控制偏差;Kp、Ki和Kd分别为比例、积分和微分的调节系数。
在Simulink建立调平平台仿真模型,使用零阶保持器,采样时间0.02 s;通过试验法确定俯仰系统比例、积分和微分调节系数分别为0.65、0.04、0,侧倾系统比例、积分和微分调节系数分别为0.85、0.04、0。丘陵山地果园坡度变化平缓,因此干扰选择5°的阶跃信号模拟路面凹坑等突变;因最大工作台俯仰和侧倾调平角度分别为10°、8°,调平角速度0.07 rad/s;确定频率0.1 Hz,幅值为10和8的正弦信号模拟连续坡度变化,此时处于平台极限工作状态。仿真结果如图5、图6所示。
从图5、图6中可以看出:俯仰、侧倾控制系统在阶跃干扰信号下能使工作台很快回到水平位置,调平时间分别为1.6、2.1 s,超调量均为0;俯仰、侧倾控制系统在正弦干扰信号下能使工作台始终保持在0°附近,波动范围分别在0.15°、0.19°内。仿真验证了平台调平的可行性,并且在不同干扰下调平性能较好,工作台角度变化较平稳,满足调平平台设计要求。
3. 安全性仿真实验验证
为测试平台的安全性和工作台调平对倾翻坡度的影响,根据结构和参数在SolidWorks中建立调平平台三维模型,导入Adams并建立倾翻实验台(图7)。实验台以一定角速度转动,平台以纵向、斜向、横向这3种不同姿态,以0、220、440、660、880、1 100 mm不同升降高度和不同载质量情况下进行仿真。当履带所受支持力为0时,实验台倾角δ记为极限倾翻角度(°),仿真结果如图8。
从图8中可以看出:升降高度和载质量一定时,平台在纵向、斜向和横向这3种姿态时的极限倾翻角度逐渐减小;且3种姿态下,升降高度和载质量增加后,平台极限倾翻角度都会减小。有调平时纵向、斜向和横向3姿态下最小倾翻角度相较于无调平时(26.36°、22.35°、20.27°)分别增加了2.47%、19.20%和24.77%,最大倾翻角度相较于无调平时无明显变化。但不同姿态、升降高度和载质量下,平台调平后对其极限倾翻角度的影响又有所不同。比如纵向时,有调平时极限倾翻角度相对于无调平时在分界线以上略有减小,分界线以下略有增加。这是因为剪叉折臂式调平方式会举升工作台,使平台重心沿坡面上移;当升降高度和载质量增加时(超过分界线时),重心位移量在平行坡面方向的投影占比会大于其在垂直坡面方向的占比,故有利于提高平台纵向上的抗倾翻能力。而在斜向和横向姿态时,有调平时极限倾翻角度较无调平时有所增加;当升降高度和载质量两个因素中的一个因素确定时,平台极限倾翻角度增量均随另一因素增加而增大。因为斜向和横向姿态调平时,均需要侧倾调平机构参与其中,使得平台重心在水平面上的投影会靠近上坡方向移动,且平台纵向尺寸大、横向小,重心移动量在横向上的占比大于纵向上占比,从而提升了平台的抗侧倾能力。
总体来说,无调平时最大和最小极限倾翻角度分别为49.85°和20.27°;有调平时最大极限倾翻角度为48.78°,最小为25.29°。整体来讲,最大极限倾翻角度(纵向、升降高度0 mm、载质量0 kg时)降低了2.15%,最小极限倾翻角度(当横向、升降高度1 100 mm、载质量300 kg时)增加了24.77%,平台安全性明显提高。
4. 结 论
(1)根据丘陵山地果园机械化发展现状设计了一款基于剪叉折臂式的自动调平平台,结构简单可靠;通过公式计算结合设计要求验证调平可行性,确定了平台各个部分的尺寸参数、液压缸位移量和平台角度的关系,计算出液压缸所需推力,建立了三维模型。对控制系统数学建模,在Simulink中采用增量式PID控制器对控制系统进行仿真,结果表明此平台有较好的抗阶跃和正弦干扰信号的能力,在极限工作状态下能始终保持在0°附近。
(2)使用Adams对在不同状态、不同举升高度和不同载质量情况下的平台极限倾翻坡度进行仿真。在仿真条件下,果园高位作业平台在调平状态下极限倾翻坡度最小值为25.29°,相较于无调平时增加了24.77%。这说明本文提出的调平机构可有效提高平台的抗倾翻能力。设计的高位调平平台能够满足丘陵山区果园的使用需求。
-
图 3 阔叶混交林各立地等级蓄积潜在生产力、现实生产力和林分可提升空间
L.立地等级。 PMIT和RMIT分别为林分平均年龄为T时,林分蓄积潜在生产力和现实生产力,ΔMIT为林分平均年龄为T时,林分蓄积年生长量可提升空间。L, site class. PMIT and RMIT represent potential productivity and realized productivity of stand volume at mean stand age T, respectively. ΔMIT represents space for improving annual growth of stand volume at mean stand age T.
Figure 3. Realized productivity,potential productivity of stand volume and their difference by each site class in broadleaved mixed forest
表 1 小兴安岭阔叶混交林林分数据
Table 1 Stand data of broadleaved mixed forest in Xiaoxing’an Mountains
指标 Index 蓄积/(m3·hm−2)
Volume/(m3·ha−1)平均年龄/a
Mean age/year密度指数/(株·hm−2)
Density index/(plant·ha−1)平均树高
Mean tree height/m最大值 Max. value 247.11 100 1974.47 22.8 最小值 Min. value 11.65 15 277.59 6.1 平均值 Mean 112.79 44 1 051.42 13.7 标准差 SD 51.44 15 367.94 3.0 表 2 林分平均高与平均年龄模型
Table 2 Stand average height and average age model
模型 Model 表达式 Expression Korf {H=a\mathrm{e}}^{{-bT}^{-c}} Richard {H=a(1-{\mathrm{e}}^{-cT})}^{b} Hossfeld H=\dfrac{a}{1 + b{T}^{-c}} Logistic H=\dfrac{a}{1 + b{\mathrm{e}}^{-cT}} Gompertz H=a{\mathrm{e}}^{{-b\mathrm{e}}^{-cT}} 注:H. 林分平均高;T. 林分平均年龄;a、b和c为模型参数。Notes: H, average stand height; T, average stand age; a, b and c are model parameters. 表 3 哑变量设置
Table 3 Dummy variable setting
立地等级
Site class变量 Variable x y z v 1 0 0 0 0 2 0 0 0 1 3 0 0 1 0 4 0 1 0 0 5 1 0 0 0 注:x、y、z、v代表哑变量参数,均通过显著性检验(P < 0.05)。Note: x, y, z, and v represent dummy variable parameters, all of which pass significance test (P < 0.05) 表 4 阔叶混交林树高曲线拟合
Table 4 Fitting of tree height curve for broadleaved mixed forest
模型
Model拟合精度 Fitting accuracy R2 RMSE Korf 0.335 9 2.441 5 Richard 0.333 4 2.445 9 Hossfeld 0.328 4 2.455 2 Logistic 0.335 1 2.442 9 Gompertz 0.330 5 2.451 4 表 5 地位级指数曲线参数
Table 5 Parameters of site class index curves
表 6 阔叶混交林不同立地等级林分蓄积生长模型参数估计值
Table 6 Parameter estimates for growth model of volume by each site class for broadleaved mixed forest
立地等级
Site classaL b′ c′ d′ 1 880.20 0.001 7 1.44 0.71 2 803.71 3 720.70 4 637.10 5 558.80 注:aL、b′、c′、d′为生长模型参数。Notes: aL, b′, c′, d′ are growth model parameters. 表 7 基准年龄30年时林分蓄积潜在生产力
Table 7 Potential productivity of stand volume at basic age of 30 years
立地等级
Site class蓄积潜在生产力/(m3·hm−2 a−1)
Potential productivity of volume/ (m3·ha−1·year−1)蓄积量/(m3·hm−2)
Volume/ (m3·ha−1)样本量
Sample size1 6.10 294.79 47 2 5.56 269.17 148 3 4.99 241.37 152 4 4.41 213.37 114 5 3.87 187.15 75 -
[1] 孟宪宇. 测树学[M]. 3版. 北京: 中国林业出版社, 2006. Meng X Y. Forest mensuration[M]. 3rd ed. Beijing: China Forestry Publishing House, 2006.
[2] 雷相东, 唐守正, 符利勇, 等. 森林立地质量定量评价:理论 方法 应用[M]. 北京: 中国林业出版社, 2020. Lei X D, Tang S Z, Fu L Y, et al. Quantitative evaluation of forest site quality: theory, method, application[M]. Beijing: China Forestry Publishing House, 2020.
[3] 雷相东, 符利勇, 李海奎, 等. 基于林分潜在生长量的立地质量评价方法与应用[J]. 林业科学, 2018, 54(12): 116−126. Lei X D, Fu L Y, Li H K, et al. Methodology and applications of site quality assessment based on potential mean annual increment[J]. Scientia Silvae Sinicae, 2018, 54(12): 116−126.
[4] 朱光玉, 康立. 森林立地生产力评价指标与方法[J]. 西北林学院学报, 2016, 31(6): 275−281. Zhu G Y, Kang L. A review of forest site productivity evaluation indicators and methods[J]. Journal of Northwest Forestry University, 2016, 31(6): 275−281.
[5] Yang Z H. An NFI-based site quality evaluation of Chinese fir plantation[J]. Journal of Sustainable Forestry, 2020, 39(2): 137−152. doi: 10.1080/10549811.2019.1623051
[6] Skovsgaard J P, Vanclay J K. Forest site productivity: a review of the evolution of dendrometric concepts for even-aged stands[J]. Forestry: An International Journal of Forest Research, 2008, 81(1): 13−31. doi: 10.1093/forestry/cpm041
[7] Fu L Y, Sharma R P, Zhu G Y, et al. A basal area increment-based approach of site productivity evaluation for multi-aged and mixed forests[J/OL]. Forests, 2017, 8(4): 119[2024−05−12]. https://doi.org/10.3390/f8040119.
[8] 张会儒, 雷相东, 张春雨, 等. 森林质量评价及精准提升理论与技术研究[J]. 北京林业大学学报, 2019, 41(5): 1−18. Zhang H R, Lei X D, Zhang C Y, et al. Research on theory and technology of forest quality evaluation and precision improvement[J]. Journal of Beijing Forestry University, 2019, 41(5): 1−18.
[9] Costa R A V, Gastón A, Bloszies A S, et al. Nature vs. nurture: drivers of site productivity in loblolly pine (Pinus taeda L.) forests in the southeastern US[J]. Forest Ecology and Management, 2024, 572(8): 122334.
[10] 郭艳荣, 吴保国, 刘洋, 等. 立地质量评价研究进展[J]. 世界林业研究, 2012, 25(5): 47−52. Guo Y R, Wu B G, Liu Y, et al. Research progress of site quality evaluation[J]. World Forestry Research, 2012, 25(5): 47−52.
[11] Hemingway H, Kimsey M. Estimating forest productivity using site characteristics, multipoint measures, and a nonparametric approach[J]. Forest Science, 2020, 66(6): 645−652. doi: 10.1093/forsci/fxaa023
[12] Milner K S, Running S W, Coble D W. A biophysical soil-site model for estimating potential productivity of forested landscapes[J]. Canadian Journal of Forest Research, 1996, 26(7): 1174−1186.
[13] 吴菲. 森林立地分类及质量评价研究综述[J]. 林业科技情报, 2010, 42(1): 12−14. Wu F. Summarize of forest site classification and quality appraisement[J]. Forestry Science and Technology Information, 2010, 42(1): 12−14.
[14] 何静, 李新建, 朱晋梅, 等. 基于最粗优势木胸径生长的湖南栎类天然林立地质量评价模型[J]. 林业科学, 2022, 58(8): 89−98. He J, Li X J, Zhu J M, et al. Site quality evaluation model of natural Quercus forests in Hunan based on the growth of the thickest dominant tree diameter at breast height[J]. Scientia Silvae Sinicae, 2022, 58(8): 89−98.
[15] 黄国胜, 马炜, 王雪军, 等. 基于一类清查数据的福建省立地质量评价技术[J]. 北京林业大学学报, 2014, 36(3): 1−8. Huang G S, Ma W, Wang X J, et al. Forestland site quality evaluation of Fujian Province based on continuous forest inventory data[J]. Journal of Beijing Forestry University, 2014, 36(3): 1−8.
[16] 沈剑波, 王应宽, 雷相东, 等. 基于BP神经网络的广东省针阔混交异龄林立地质量评价[J]. 北京林业大学学报, 2019, 41(5): 38−47. Shen J B, Wang Y K, Lei X D, et al. Site quality evaluation of uneven-aged mixed coniferous and broadleaved stands in Guangdong Province of southern China based on BP neural network[J]. Journal of Beijing Forestry University, 2019, 41(5): 38−47.
[17] 徐罗, 亢新刚, 郭韦韦, 等. 天然云冷杉针阔混交林立地质量评价[J]. 北京林业大学学报, 2016, 38(5): 11−22. Xu L, Kang X G, Guo W W, et al. Site quality evaluation of natural spruce-fir and broadleaf mixed stands[J]. Journal of Beijing Forestry University, 2016, 38(5): 11−22.
[18] 雷相东, 符利勇, 李海奎, 等. 森林立地质量评价技术规程: LY/T 3315-2022[S]. 北京: 国家林业和草业局, 2020. Lei X D, Fu L Y, Li H K, et al. Technical regulation for forest site quality evaluation: LY/T 3315-2022[S]. Beijing: National Forestry and Grassland Administration, 2020.
[19] Facundo J O, Cecilia C, Matías G G, et al. Modeling potential site productivity for Austrocedrus chilensis trees in northern Patagonia (Argentina)[J/OL]. Forest Ecology and Management, 2022, 524(46): 120525[2024−04−23]. https://doi.org/10.1016/j.foreco.2022.120525.
[20] Raul V, Duncan J G, Pedro M M, et al. Site quality for mahogany (Swietenia macrophylla King) in natural forests in Quintana Roo[J]. Journal of Sustainable Forestry, 2022, 41(1): 33−50. doi: 10.1080/10549811.2020.1841004
[21] 沈剑波, 雷相东, 雷渊才, 等. 长白落叶松人工林地位指数及立地形的比较研究[J]. 北京林业大学学报, 2018, 40(6): 1−8. Shen J B, Lei X D, Lei Y C, et al. Comparison between site index and site form for site quality evaluation of Larix olgensis plantation[J]. Journal of Beijing Forestry University, 2018, 40(6): 1−8.
[22] 高若楠, 谢阳生, 雷相东, 等. 基于随机森林模型的天然林立地生产力预测研究[J]. 中南林业科技大学学报, 2019, 39(4): 39−46. Gao R N, Xie Y S, Lei X D, et al. Study on prediction of natural forest productivity based on random forest model[J]. Journal of Central South University of Forestry & Technology, 2019, 39(4): 39−46.
[23] 段光爽. 基于生长过程的吉林省林分类型分类与潜在生产力预估[D]. 北京: 中国林业科学研究院, 2018. Duan G S. Growth process-based clustering of forest stand types and their potential productivity estimation in Jilin Province[D]. Beijing: Chinese Academy of Forestry, 2018.
[24] 段光爽, 郑亚丽, 洪亮, 等. 基于潜在生产力的华北落叶松纯林和白桦山杨混交林立地质量评价[J]. 林业科学, 2022, 58(10): 1−9. Duan G S, Zheng Y L, Hong L, et al. A potential productivity-based approach of site quality evaluation for larch pure forest and birch-aspen mixed forest[J]. Scientia Silvae Sinicae, 2022, 58(10): 1−9.
[25] 刘丹. 基于分布适宜性和潜在生产力的定量适地适树研究[D]. 北京: 中国林业科学研究院, 2018. Liu D. Qunatitative assessment of matching trees to sites based on both distribution suitability and potential site productivity[D]. Beijing: Chinese Academy of Forestry, 2018.
[26] Roberto J P, Carola P, Javier J C, et al. Evaluation of the site form as a site productive indicator in temperate uneven-aged multispecies forests in Durango, Mexico[J]. Plants, 2022, 11(20): 2764. doi: 10.3390/plants11202764
[27] 唐守正, 郎奎建, 李海奎. 统计和生物数学模型计算: ForStat教程[M]. 北京: 科学出版社, 2009: 343−347. Tang S Z, Lang K J, Li H K. Computation for statistics and biomathematical models: ForStat course[M]. Beijing: Science Press, 2009: 343−347.
[28] 唐守正, 李希菲, 孟昭和. 林分生长模型研究的进展[J]. 林业科学研究, 1993, 6(6): 672−679. Tang S Z, Li X F, Meng Z H. The development of studies on stand growth models[J]. Forest Research, 1993, 6(6): 672−679.
[29] 唐守正, 李希菲. 用全林整体模型计算林分纯生长量的方法及精度分析[J]. 林业科学研究, 1995, 8(5): 471−476. Tang S Z, Li X F. Precision analysis on growth rates estimated by integrated whole stand growth and yield model[J]. Forest Research, 1995, 8(5): 471−476.
[30] 董灵波, 邵威威, 田栋元, 等. 基于林木分级的大兴安岭天然兴安落叶松树高曲线研究[J]. 北京林业大学学报, 2023, 45(5): 88−96. doi: 10.12171/j.1000-1522.20210513 Dong L B, Shao W W, Tian D Y, et al. Height curve of natural Larix gmelinii in the Daxing’anling Mountains of northeastern China based on forest classification[J]. Journal of Beijing Forestry University, 2023, 45(5): 88−96. doi: 10.12171/j.1000-1522.20210513
[31] 王鹤智. 东北林区林分生长动态模拟系统的研究[D]. 哈尔滨: 东北林业大学, 2012. Wang H Z. Dynamic simulating system for stand growth of forests in Northeast China[D]. Harbin: Northeast Forestry University, 2012.
[32] 刘鑫. 长白山天然云冷杉针阔混交林立地质量评价研究[D]. 北京: 北京林业大学, 2017. Liu X. Evaluation of site quality in the natural spruce-fir mixed forests of Changbai Mountains[D]. Beijing: Beijing Forestry University, 2017.
[33] 宋长江. 帽儿山实验林场天然次生林立地质量综合评价[D]. 哈尔滨: 东北林业大学, 2023. Song C J. Comprehensive evaluation of natural secondary forest site quality in Maoershan Experimental Forest Farm[D]. Harbin: Northeast Forestry University, 2023.
[34] Hua W P, Pan X, Zhu D H, et al. Developing growth and harvest prediction models for mixed coniferous and broad-leaved forests at different ages[J/OL]. Forests, 2023, 14(7): 1416[2024−05−12]. https://doi.org/10.3390/f14071416.