下载:
-
中国竹林面积大，竹材资源丰富。据统计，2022年中国竹林面积高达709.92 万hm2，但由于竹子主要分布在海拔200~1 200 m[1]，受竹林立地条件和林道建设滞后制约，将竹子从山上采伐运输出来十分困难。长期以来，人工砍伐竹材后，通过人力或畜力将竹材运至山下，生产方式落后，劳动强度大、效率低[2]。为降低劳动强度，有人把钢丝绳作为运送索道来实现快速运竹，但运送索道前期布置比较困难，且对悬索和跑车的使用寿命及安全要求较高[3]。当前，中国林业索道集材领域跑车主要分为无动力源和自带动力源两大类。无动力源跑车不具备自主起吊或装卸货物的能力，此类索道装卸货物需人力或其他外力，功能单一，效率不高。自带动力源跑车可起吊和自主装卸货物，但此类跑车机械结构复杂，质量大，且因蓄电池容量限制，跑车的续航能力不足，若在工作途中蓄电池电量耗尽，易导致运输过程发生货物掉落等安全事故[4]。如何将采伐的竹材高效安全、低成本运输下山成为亟待解决的技术瓶颈。
本研究对竹材索道集材跑车运输系统中悬索和跑车等关键部件进行分析，设计了一种结构简单、续航能力强、运输安全性较高的自充电索道跑车，绘制了跑车三维模型，对跑车关键部件进行了有限元分析，对跑车-索道运输系统进行了仿真试验。自充电跑车设计大幅提高跑车续航能力和安全性，为竹材运输下山装备开发提供一种新的思路。
-
根据在湖南、福建等省部分竹林地理环境的航拍调研结果，以及竹材质量和结构特点，初步确定竹材索道集材跑车运输系统基本设计参数：竹林垂直高度≤200 m，索道倾斜角≤35°，索道跨距200~500 m，承载索速度≤1.0 m·s−1 (可调)，牵引索与回空索速度≤4.0 m·s−1 (可调)，跑车质量200 kg，运输跑车单次载质量≤400 kg。
-
在现代林业索道集材运输领域中，对工作索的要求为：具有可靠的静态抗拉力和交变的动态承载力；具有良好的抗磨性，在跑车滑轮接触、转动、滑动、挤压作用下不易损耗；具有适宜的弹性，易于吸收冲击。本设计选用“6×19股纤维芯”钢丝绳作为承载索，根据跑车和吊运质量、索道倾斜角等主要设计参数初步确定承载索规格。承载索金属截面积FMB可由下式计算[5]：
$$ {F}_{{\rm{MB}}}=\frac{2Pn\mathrm{cos}\beta \left({l}_{0}+f\mathrm{sin}2\beta \right)}{8\left[\sigma \right] f {\mathrm{cos}}^{2}\beta -0.1\mathrm{ln}\left({l}_{0}+2f\mathrm{sin}2\beta \right)}\mathrm{。} $$ (1) 式(1)中：P为总载荷(N)；n为安全系数，对于林业索道承载索，n一般取2；l0为承载索有效长度(m)；β为索道倾斜角；f为有荷最大挠度，f =Sl0，S为有荷最大挠度系数，S一般取0.03；[σ]为承载索钢丝有效强度，[σ]=0.85σb，σb为极限强度，σb=2 160 MPa。总载荷：
$$ P=Q(1+a)。 $$ (2) 式(2)中：Q为跑车与吊运质量(kg)；a为冲击载荷系数，查表取a=0.3。将Q=6 000 N代入式(2)可得：P=7 800 N。将P=7 800 N，[σ]=1 836 MPa，最大倾斜角β=35°等参数代入式(1)得：FMB=88.96 mm2。
根据承载索金属截面积，可计算承载索直径(d)：
$$ d=\sqrt{\frac{1.27{F}_{{\rm{MB}}}}{{K}'_{{\rm{n}}}}}\mathrm{。} $$ (3) 式(3)中：${K}_{{\rm{n}}}'$为钢索结构坚固系数，带有麻芯的钢索${K}'_{{\rm{n}}}\approx 0.6$。将FMB=88.96 mm2代入式(3)得：d=13.7 mm，取承载索钢丝绳直径d=14.0 mm。承载索型号为：14 6×19 NFC 1550 B ZZ，其单位长度重力6.85 N·m−1，许可破断拉力112 000 N，横截面面积72.49 mm2，弹性模量105 MPa。
-
①行走机构。行走机构设计的关键是行走滑轮的选择，行走滑轮个数和各滑轮间距大小对承载索弯曲应力影响很大。通过应变片试验测试，行走滑轮作用下承载索弯曲应力分布[5]如图1所示。当跑车吊重在承载索上运行时，行走轮正下方弯曲应力最大，弯曲应力随载荷作用点距离变大而逐渐减小。图1A比图1B应力分布情况更加合理，虽然弯曲应力最大值不变，但其弯曲应力的变化曲线较为平缓。因此，任意加大跑车行走轮的间距，对承载索的受力非常不利。为延长承载索使用寿命，应把行走轮布置得相互靠近并使之具有相等的间距。通常行走轮纵向轮距应为行走轮轮径2倍左右。本设计跑车行走轮间距取700 mm。行走轮个数(N)可由最大弯曲应力(${\sigma }_{{\rm{w}}}$)求得[5]：
图 1 行走滑轮作用下承载索弯曲应力
Figure 1. Bending stress of bearing cable under the action of traveling pulley
$$ N=\frac{P}{{\sigma }_{{\rm{w}}}}\sqrt{\frac{{E}_{{\rm{k}}}}{T{F}_{\mathrm{M}\mathrm{B}}}}\mathrm{。} $$ (4) 式(4)中：P为总载荷(N)；σw为承载索最大弯曲应力(N·mm−2)；Ek为承载索弹性模量(MPa)；T为承载索最大张力(N)；FMB为钢索横截面面积(mm2)。承载索最大张力(T)可由承载索的最大破断拉力(TP)求得：
$$ T\leqslant \frac{{T}_{{\rm{P}}}}{K}\mathrm{。} $$ (5) 式(5)中：TP为承载索最大破断拉力，TP=φS，φ为折断系数，φ=0.82~0.90，此处取φ=0.90，S为钢丝绳破断拉力(N)；K为考虑冲击系数后承载索的安全系数，取K=2.5。将TP=100 800 N代入式(5)，得T=40 320 N，承载索最大拉应力$ {\sigma }_{{\rm{i}}}=\dfrac{T}{{F}_{{\rm{MB}}}}= $556 MPa。由式(4)可知：拉应力与承载索张力T成正比，而弯曲应力与承载索最大张力T成反比，考虑到承载索总应力不得超过承载索钢丝绳的有效强度，一般规定：
$$ \sigma_{{\rm{w}}}\leqslant (0.5\sim 1.0)\sigma_{{\rm{i}}}。 $$ (6) 本次设计取$ {\sigma }_{{\rm{w}}}=0.5{\sigma }_{{\rm{i}}}=278 $ MPa。将P=7 800 N，$ {\sigma }_{w}=278 $ MPa，T、$ {E}_{\mathrm{k}} $和$ {F}_{\mathrm{M}\mathrm{B}} $等参数代入式(4)得$ N=5.3 $，取$ N $为6个。②升降卷扬机构。升降卷扬机构指安装于跑车车架底板上对称布置的2台升降卷扬机，主要包括电机、齿轮减速机构、刹车装置、起重索、防脱挂钩等。通过控制卷扬机起重索的升降实现林间集材和集材点打捆竹材的提升运输。根据索道集材跑车运输系统设计要求，跑车升降卷扬机选用坤航起重机械有限公司生产的迈士顿小型升降卷扬机。卷扬机参数如表1。③自充电机构。自充电机构主要由链传动装置、发电机和蓄电池组成。跑车行走轮在承载索上转动时，前后2个行走轮通过与其同轴连接的链传动装置将机械能传递给发电机，发电机将转换的电能传输至蓄电池，实现跑车边工作边充电，提高跑车的续航能力。跑车整体结构如图2所示。
表 1 升降卷扬机参数
Table 1. Lifting winch parameters
额定起
重/kg电机功
率/W额定电
压/V起升速度/
(m·min−1)卷起长
度/m挂钩质
量/kg整机质
量/kg200~400 1600 220 13~20 1~50 16.7 29 
图 2 跑车整体结构示意图
Figure 2. Overall structure diagram of sports car
-
集材作业时，主机上的卷扬机控制牵引索与回空索运行，牵引跑车在承载索上移动至林间集材点上空，并通过控制牵引索与回空索实现跑车在集材点上空的定位。当跑车固定于林间集材点上空时，用遥控器控制跑车升降卷扬机放索，人工拽拉绳索至林间伐倒竹材处，将卡住竹材根部的卡具挂在绳索上，跑车卷扬机收索，将竹子牵引至林间集材点。在集材点进行竹材打捆，用遥控器控制跑车升降卷扬机收索，提升打捆竹材至限位高度。再控制主机上的牵引索与回空索运行，牵引跑车运输竹材至转运点。
-
由于跑车与重物的整体质量全部通过行走轮支撑在承载索上，行走轮受力较大。若在工作过程中失效，将造成安全事故，所以需对行走轮进行静力学分析[6]。行走轮采用的灰铸铁HT250制作，其硬度低于钢丝绳。HT250力学特性如表2所示。
表 2 HT250材料特性
Table 2. HT250 material properties
材料 泊松比 密度/
(kg·m−3)抗拉强度/
Pa弹性模量/
PaHT250 0.3 7.15×103 2.5×108 9×1010 利用ABAQUS自带的四面体网格法[7]对行走轮进行网格划分。图3为网格划分图，网格节点总数为247 062个，单元数为82 354个。
图 3 行走轮网格划分图
Figure 3. Grid division diagram of walking wheel
-
行走轮工作时，围绕承载索垂直方向转动，因此，设置约束时，释放X轴方向的转动，其余方向的移动与转动均固定。
跑车工作时，跑车与重物总载荷全部作用在行走轮轴上，行走轮轴再作用于行走轮。设行走轮轴施加的压力均匀作用于轮孔上。6个行走轮受总载荷为7 800 N，则作用于中心孔的均布压力$ {G}_{\mathrm{n}} $：
$$ {G}_{\mathrm{n}}=\frac{{P}{{'}}}{{\text{π}}Dh}=138\;004.25\;\mathrm{P}\mathrm{a}\mathrm{。} $$ (7) 式(7)中：P′为单个行走轮中心孔压力(N)；D为行走轮中心孔直径，为50 mm；h为行走轮轮宽，为60 mm。
-
跑车行走轮在施加约束和载荷后，经求解运算得到行走轮等效应力最小值为3.121×10−3 MPa，最大值为11.84 MPa；等效位移最小值为0 m，最大值为1.778×10−4 m。等效应力和等效位移云图如图4A和图4B。
图 4 等效应力云图(A)和等效位移云图(B)
Figure 4. Equivalent stress nephogram and equivalent displacement nephogram
由图4A可知：行走轮的最大应力出现在轮中心孔圈中间处，为11.84 MPa；最小应力出现在行走轮轮缘最外侧，行走轮压力沿2侧轮缘变化平稳。由图4B知：行走轮的最大位移发生在外侧车轮轮缘和轮中心孔圈中间处，变形量为1.778×10−4 m。由于灰铸铁HT250的许用应力为250.00 MPa，远大于工作时所受的最大应力，且行走轮最大位移量对其工作状态影响极小，因此强度和刚度均满足设计要求。
-
本设计中承载索钢丝绳是由6绳股捻绕而成，具有一定的挠度。进行动力学分析时，需先把钢丝绳简化。在应用ADAMS研究耦合振动时，还需要对钢丝绳进行柔性化处理。此处选择利用ABAQUS软件将简化后的钢丝绳模型进行网格划分，模态处理后保存为MNF过渡文件格式[8]，再将文件导入ADAMS软件中，获得钢丝绳的柔性体。
-
将SolidWorks建立的跑车和模态处理后的承载索三维模型导入ADAMS 2020应用软件，在ADAMS中给跑车-承载索平移驱动motion设置step函数控制其运行速度[9]，当跑车速度为2 m·s−1时，step函数为step(time, 0, 0, 2, 2)+step(time, 250, 0, 252, −2)，即前2 s跑车速度由0 m·s−1加速到2 m·s−1，2~250 s以2 m·s−1保持匀速运动，接近承载索端点时，由2 m·s−1减速到0 m·s−1。承载索钢丝绳长度为500 m，取250 m处钢丝绳为研究点，跑车速度设置2和4 m·s−1，跑车总质量取200、400、600 kg，索道倾斜角为20°，通过ADAMS运动仿真得到承载索的垂直位移时间历程如图5所示。
图 5 不同跑车总质量下承载索的垂直位移时间历程
Figure 5. Vertical displacement time history of the load-bearing cable under different total mass of sports cars
图5可以看出：随着速度增加，跑车走完全程完成1次竹材运输所需时间减少。承载索垂直位移随跑车速度的增加而增大，且通过图5A与图5B对比发现：速度为2.0 m·s−1时位移变化比4.0 m·s−1时密集，此时跑车在承载索上运行时对承载索施加的载荷频率变化快。当速度为4.0 m·s−1时，虽位移变化频率比2.0 m·s−1时小很多，但此时振幅剧增，位移最大可达4.39 m，容易导致跑车运输的货物掉落或承载索因挠度过大发生断裂。跑车运输工作过程中应尽量避免在此工况下工作。速度一定时，对比跑车在不同载荷下索道垂直位移变化发现，随着跑车载荷的增加，承载索垂直位移越大。载荷一定时，随着跑车速度增加，位移急剧增大。因此，载运量越大，速度越快，往返运输竹子时间越短，效率越高，但此时跑车与承载索的位移大，可能导致跑车运输货物掉落或索道断裂风险，安全性较低。为保证跑车运输竹材安全性和高效率，需要确定跑车优化运输条件。
图 6 交互作用对综合评分指标(y3)的影响
Figure 6. Influence of interaction on comprehensive evaluation index (y3)
-
跑车运输竹材的仿真分析结果表明：影响竹材运输效率和安全性的主要因素为运输速度和运输质量。此外，索道的倾斜角对运输安全性也具有较大影响[10]，但因索道运输时索道倾斜角条件主要取决于竹林场实际山体高度，通过增大或减小山顶与山底塔架的高度对倾斜角的影响不大。若要将塔架增减到影响倾斜角的高度，将增加安装运输成本，且随着索道支撑塔架的增高，运输安全性将受到极大影响。因此，竹材索道集材跑车运输系统主要应用于倾斜角为20°±4°的竹林地环境。综合考虑实际作业过程中运输速度、跑车与竹材总质量和索道倾斜角的变化对整个运输系统效率和安全性的影响，基于第3节中ADAMS构建的索道跑车运输仿真模型，开展Box-Behnken仿真试验以确定优化参数组合[11]。
-
以运输速度(x1)、跑车与竹材总质量(x2)、索道倾斜角(x3)为试验因素，以索道最大位移(y1)和运输效率(y2)为评价指标进行仿真试验，因素水平编码如表3所示。因试验存在2个评价指标，且2个指标之间存在一定的矛盾，即安全性越高，y1应尽可能小；而运输效率越高，则y2要偏大。此处采用多指标问题处理的综合评分法中排队评分法，对2个试验指标进行综合评分(表4)，将多指标综合为单指标进行分析[12]。
表 3 试验因素编码
Table 3. Test factor coding
编码 因素 运输速度(x1)/
(m·s−1)跑车与竹材总质量(x2)/
kg倾斜角(x3)/
(°)−1 2 200 16 0 3 400 20 1 4 600 24 $$ y_{2}=x_{2}/t。 $$ (8) 式(8)中：t为运输总时间(s)。
表 4 试验结果
Table 4. Test result
序号 运输速度
(x1)/(m·s−1)总质量
(x2)/kg倾斜角
(x3)/(°)最大位移
(y1)/m效率
(y2)/(kg·s−1)综合评分
(y3)序号 运输速度
(x1)/(m·s−1)总质量
(x2)/kg倾斜角
(x3)/(°)最大位移
(y1)/m效率
(y2)/(kg·s−1)综合评分
(y3)1 2 400 24 2.47 1.59 8.92 10 4 400 16 2.89 3.10 11.34 2 3 600 16 2.97 3.54 12.16 11 3 400 20 3.01 2.36 9.30 3 3 400 20 3.05 2.36 9.20 12 2 600 20 2.98 2.38 9.18 4 3 200 16 1.37 1.18 10.88 13 2 400 16 2.17 1.59 9.72 5 3 600 24 4.05 3.54 9.30 14 3 400 20 3.11 2.36 9.04 6 3 200 24 1.42 1.18 10.75 15 3 400 20 2.95 2.36 9.46 7 4 600 20 4.39 4.65 11.00 16 3 400 20 2.92 2.36 9.54 8 2 200 20 0.98 0.79 11.00 17 4 400 24 3.77 3.10 9.01 9 4 200 20 1.65 1.55 10.99 根据仿真试验结果，对于最大位移指标，第8号试验为0.98 m，位移量最小，赋予10分；第7号试验为4.39 m，位移量最大，赋予1分。对于效率指标，第7号试验为4.65 kg·s−1，效率最高，赋予10分；第8号试验为0.79 kg·s−1，效率最低，赋予1分。对其他试验号各指标值的得分，根据其与该指标优秀值的差异按比例打分，最后对每号试验各指标的分数相加即得综合评分。
-
仿真试验结果如表4所示。采用Design-Expert 13对表4数据进行回归分析[13]，得到评价跑车运输安全性与运输效率的综合评分指标(y3)关于运输速度(x1)、跑车与竹材总质量(x2)和索道倾斜角(x3)二次回归模型为：$ y_{3}=9.31+0.44x_{1}-0.247\;5x_{2}-0.765\;0x_{3}+0.457\;5x_{1}x_{2}-0.382\;5x_{1}x_{3}-0.682\;5x_{2}x_{3}+ 0.104\;7x_{1}^{2}+ 1.130\;0x_{2}^{2} +0.334\;8x_{3}^{2} $。回归方程的方差分析如表5所示。回归方程模型显著性检验F1=25.15，P＜0.001 0，表明回归方程极为显著，与试验值拟合程度好，模型合理。模型的失拟性检验F2=2.82，P=0.1708＞0.05，说明回归方程与实际误差所占比例很小，即试验误差小。由方差分析结果可看出：各因素对跑车运输安全性与运输效率综合影响由大到小依次为索道倾斜角、运输速度、跑车与竹材总质量。
表 5 综合评分指标(y3)方差分析
Table 5. Analysis of variance of comprehensive evaluation index (y3)
变异来源 平方和 自由度 均方 F P 变异来源 平方和 自由度 均方 F P 模型 16.180 0 9 1.800 0 25.150 0 0.000 2** x12 0.046 2 1 0.046 2 0.646 3 0.447 6 x1 1.550 0 1 1.550 0 21.670 0 0.002 3** x22 5.370 0 1 5.370 0 75.180 0 ＜0.000 1** x2 0.490 0 1 0.490 0 6.860 0 0.034 5* x32 0.471 8 1 0.471 8 6.600 0 0.031 1* x3 4.680 0 1 4.680 0 65.500 0 ＜0.000 1** 残差 0.500 4 7 0.071 5 x1x2 0.837 2 1 0.837 2 11.710 0 0.011 1* 失拟 0.339 9 3 0.113 3 2.820 0 0.170 8 x1x3 0.585 2 1 0.585 2 8.190 0 0.024 3* 误差 0.160 5 4 0.040 1 x2x3 1.860 0 1 1.860 0 26.070 0 0.001 4** 总和 16.680 0 16 说明：**表示极显著(P＜0.01)；*表示显著(P＜0.05)。 -
为探究各因素的交互 作用，利用Design-Expert 13建立运输速度、跑车与竹材总质量及索道倾斜角的交互因素对综合评分指标(y3)的响应面曲线，结果如图6所示。由图6A可知：运输速度一定时，随着跑车与竹材总质量增大，y3先减小后增大，最佳运输总质量为400~600 kg；总质量一定时，y3与运输速度呈正相关，最佳运输速度为3.0~4.0 m·s−1，运输速度和跑车与竹材总质量的交互作用显著。由图6B可知：运输速度一定时，y3与倾斜角呈负相关，随着倾斜角的增大而减小，最佳倾斜角为16.4°~20.8°；倾斜角一定时，y3与运输速度呈正相关，随运输速度的增大而增大，最佳运行速度为2.8~4.0 m·s−1，跑车运输速度与倾斜角的交互作用不显著。由图6C可知，倾斜角一定时，随着跑车与竹材总质量增大，y3先减小后增大，最佳运输总质量为420~580 kg；总质量一定时，y3与倾斜角呈负相关，随倾斜角的增大而减小，最佳倾斜角度为16.5°~20.4°，跑车与竹材总质量和倾斜角的交互作用显著。综合各交互因素响应曲面分析可得：当索道倾斜角为16.5°~20.4°、运输速度为3.0~4.0 m·s−1、跑车与竹材总质量为420~580 kg时，综合考虑的跑车运输安全性和运输效率效果较好。
-
为明确在跑车运输安全性和运输效率综合达到最优情况下的最佳参数组合，以y3值最大为优化目标，采用Design-Expert 13中的优化模块对3种因素在约束条件下进行求解，优化约束条件为：
$$ \mathrm{M}\mathrm{a}\mathrm{x}{y}_{3}\left({x}_{1}\mathrm{、}{x}_{2}\mathrm{、}{x}_{3}\right) ,{\rm{s}}.{\rm{t}}. \left\{\begin{array}{c}2.0\;{\rm{m}}\cdot {\mathrm{s}}^{-1}\leqslant {x}_{1}\leqslant 4.0\;{\rm{m}}\cdot {\mathrm{s}}^{-1}\\ 200\;{\rm{kg}}\leqslant {x}_{2}\leqslant 600\;{\rm{kg}}\\ 16\text{°}\leqslant {x}_{3}\leqslant 24\text{°}\end{array}\right. 。 $$ (9) 求解得最佳参数组合：运输速度为3.95 m·s−1，跑车与竹材总质量为576.67 kg，索道倾斜角为17.09°，此时y3为12.31，达到最大，跑车运输安全性和运输效率综合达到最佳。
-
针对山地竹材砍伐后运输下山困难，成本高，现行简易索式、轨道式等机械设备多由木材索式、轨道式改造而来，效果不稳定，实用性不高，安全得不到保障等问题，提出了一种专用于山地竹材运输的索道集材自充电跑车运输系统。通过实地调研，确定了竹材索道集材跑车运输系统基本设计参数，对悬索及跑车整体尺寸、功能进行设计，并对跑车关键受力部件进行强度和刚度校核。经有限元静力学分析可知：运输过程中跑车行走轮所受最大应力为11.84 MPa，远小于行走轮许用应力，最大变形量为1.778×10−4 m，对跑车正常工作影响极小，强度和刚度均满足设计要求。使用ABAQUS和ADAMS软件对跑车与索道组成的系统进行了振动耦合分析，探究了运输速度、跑车与竹材总质量对运输安全性和效率的影响。当跑车速度越快，吊运竹材质量越大时，效率越高，但此时安全性急剧下降，最大位移量可达4.39 m，易发生安全事故。
承载索倾斜角为16.5°~20.4°、运输速度为3.0~4.0 m·s−1、跑车与竹材总质量为420~580 kg时，跑车运输安全性和运输效率均较好；运输速度为3.95 m·s−1、跑车与竹材总质量为576.67 kg、索道倾斜角为17.09°时，运输安全性和效率综合达到最佳。
Design and analysis of bamboo cableway skidding transportation equipment sports car and suspension cable
-
摘要:
目的 针对山地竹子砍伐后运输下山困难，现行简易索式机械设备效果不稳定，运输跑车续航能力不足，安全得不到保障等问题，设计了一种竹材索道集材跑车运输系统，以实现竹材下山安全、高效运输。 方法 对竹材资源丰富的福建、湖南等竹林场进行调研，确定了竹材索道集材跑车运输系统基本设计参数，设计了悬索和跑车整体功能、尺寸及关键部件，对跑车主要承载部件进行了有限元分析，校核了跑车强度刚度；利用ADAMS软件对索道-跑车系统进行耦合振动分析，确定索道运输竹材过程的主要影响因素。选取运输速度、跑车与竹材总质量、索道倾斜角为考察因素，以跑车运输竹材安全性和效率综合评分为指标，开展Box-Behnken仿真试验，求解优化参数组合。 结果 二次回归正交旋转组合试验结果表明：三因素对跑车运输安全性和效率的影响从大到小依次为索道倾斜角、运输速度、跑车与竹材总质量。当运输速度为3.95 m·s−1、跑车与竹材总质量为576.67 kg，索道倾斜角为17.09°时，运输安全性和效率综合达到最佳。 结论 当承载索倾斜角为16.50°~20.40°、运输速度为3.00~4.00 m·s−1、跑车与竹材总质量为420.00~580.00 kg时，运输安全性和运输效率较好，可在竹林索道集材运输作业的同时充电蓄能、竹材安全高效运输下山。图6表5参13 Abstract:Objective In response to the difficulties in transporting bamboo down the mountain after logging, the unstable effect of current simple cable type mechanical equipment, the insufficient endurance of transportation sports cars, and the lack of safety assurance, this study aims to design a bamboo cableway skidding sports car transportation system, so as to achieve safe and efficient transportation of bamboo down the mountain. Method Through research on bamboo forest farms in Fujian, Hunan Provinces with abundant bamboo resources, the basic design parameters of the bamboo cableway skidding and sports car transportation system were determined. The overall functions, dimensions, and key components of the suspension cable and sports car were designed. The main load-bearing components of the sports car were analyzed by finite element analysis, and the strength and stiffness of the sports car were verified. The coupling vibration of cableway sports car system was analyzed by ADAMS software, and the main influencing factors in transporting bamboo by cableway were determined. Taking transportation speed, the total weight of sports cars and bamboo, and the inclination angle of the cableway as evaluation factors, and comprehensive evaluation of the safety and efficiency of sports car transportation of bamboo as indicators, Box-Behnken simulation experiments were conducted to optimize parameter combination. Result The results of the quadratic regression orthogonal rotation combination experiment indicated that the three factors which had significant impact on the safety and efficiency of sports car transportation in descending order were cableway inclination angle, transportation speed, and total weight of sports car and bamboo. When the transportation speed was 3.95 m·s−1, the total weight of the sports car and bamboo was 576.67 kg, and the inclination angle of the cableway was 17.09°, the comprehensive transportation safety and efficiency reached the best. Conclusion The transportation safety and efficiency are good when the inclination angle of the carrying cable is 16.50° − 20.40°, the transportation speed is 3.00 − 4.00 m·s−1, and the total weight of the sports car and bamboo is 420.00 − 580.00 kg. It can achieve charging and energy storage during the bamboo cableway skidding transportation operation, and ensure the safe and efficient transportation of bamboo down the mountain. [Ch, 6 fig. 5 tab. 13 ref.] -
枇杷黄毛虫Melanographia flexilineata，又名枇杷瘤蛾，属鳞翅目Lepidoptera灯蛾科Arctiidae[1]，是中国南方枇杷Eriobotrya japonica上最主要的害虫。除危害枇杷外，枇杷黄毛虫还危害梨Pyrus spp.，李Prunus spp.，合欢Albizzia julibrissin，紫薇Lagerstroemia india等果树和绿化树。浙江省杭州市余杭区塘栖镇从唐代起就以盛产枇杷而著称于世[2]。近年来，鲜果枇杷价格快速上升，对品质的要求日益提高，枇杷产业迅速发展。枇杷黄毛虫在浙江余抗塘栖1 a发生4代，幼虫危害枇杷嫩芽和嫩叶，发生多时也危害老叶、嫩茎表皮和花果，严重时全树叶片被吃尽，削弱树势[3-4]，影响产量。国内对枇杷黄毛虫的寄主分布、形态、生物学特征等有一定的报道[5-9]，但开展系统监测调查、对其种群动态的研究及应用统计学方法建立预测模型等方面鲜有报道。本研究采用测报灯诱集与田间调查相结合的方法，通过对塘栖枇杷黄毛虫成虫灯下逐日系统监测结合田间定期系统调查，选择不同时期虫口基数、气象资料（气温、相对湿度、降水量）作为预测因子，使用逐步回归法开展其发生期和发生量的预测预报模型研究，对保障枇杷产量和品质，促进枇杷产业发展具有重要意义。
1. 材料和方法
1.1 诱捕设备
自动虫情测报灯由河南省佳多科工贸有限公司生产，灯下接诱导漏斗和接虫袋。
1.2 虫情数据的取得
试验地设在余杭区塘栖枇杷研究所试验基地，枇杷品种为‘红娘子’‘Hongniangzi’和‘白沙’‘Baisha’，试验时期为2008-2013年。枇杷黄毛虫成虫采用频振式测报灯进行逐日诱集，每年开灯时间为4月1日至9月30日，历时约183 d，幼虫隔3 d左右田间调查1次，选取样树5株进行调查，树冠按东南西北划分4个方位，各个方位分上层、中层、下层等3层，将样树划分为12个资源单位，选取新梢枝条1枝·资源单位-1，记载各枝条上的幼虫数量。统计、记录每天诱捕的黄毛虫成虫数量及每周田间查见的黄毛虫幼虫，进行系统监测。
1.3 数据的处理和预报模型的建立
利用2008-2013年枇杷黄毛虫在余杭区塘栖镇监测的系统历史资料，将前5 a的诱虫数据、田间系统调查数据和气象资料用来建模，应用SPSS 17.0软件，采取逐步回归分析方法[10]，建立发生期和发生量的预测预报模型，最后1 a的资料用来检验。原始调查数据不作任何转换，所用气象资料由杭州市气象局提供。
1.4 预测模型的检验
将2014年的实况资料，应用唐启义等[11]提出的病虫测报应验程度判定模式进行验证，从而判断组建模型的可行性。发生期预报应验程度判定模式：
${{T}_{s}}={{100}_{e}}^{\frac{-\pi }{\ln t}{{\left( \frac{{{a}_{1}}-{{a}_{2}}}{\delta } \right)}^{2}}}。$
发生量预报应验程度判定模式：
${{D}_{s}}=100e\frac{-2a\pi }{\left( a+{{a}_{1}} \right)\ln t}\bullet \left( \frac{{{a}_{1}}-{{a}_{2}}}{\delta } \right)。$
其中：Ts和Ds分别为发生期和发生量判定模式的分值，Ts和Ds<40表明预报不准确，40≤Ts和Ds< 60表明预报比较准确，Ts和Ds≥60表明预报准确；a为预报对象常年平均值，a1为实测值，a2为预测值；δ为预报对象常年标准差，t 为自预报发出至实际发生时的期距（d）。
2. 结果与分析
2.1 田间种群数量特征
2.1.1 灯下虫情特征
2008-2013年枇杷黄毛虫灯下成虫诱集量监测结果（图 1）显示：2008年以来，枇杷黄毛虫种群动态特征和数量结构发生了一定的变化：一是越冬代成虫始见期年度间差异较大，其中2008年、2009年和2013年成虫始见期较早，均为4月中旬，以2013年最早（4月14日），2010-2012年相对较迟。二是主害代灯下成虫峰期时间、蛾量年度间有差异，总体均以第3代峰期诱蛾量最高，其次为第2代；峰期持续时间、蛾量受当年气候等条件的影响呈相应的变化，2008年和2009年于7月中旬出现全年诱蛾最高峰，2010-2013年诱蛾最高峰出现在7月下旬；主害期危害程度年度间差异大，2011年7月下旬诱蛾量596头·灯-1·旬-1，与历年同期比第3代成虫峰期蛾量明显偏高，占全年诱集总量的55.18%，2013年7月高峰期蛾量为87头·灯-1·旬-1，明显低于历年同期，分析其原因可能与当年夏季出现罕见高温干旱天气，不利于枇杷黄毛虫生长发育有关。有研究表明：气候较为暖湿的条件下成虫产卵量多、幼虫发生量大[6]，低于20 ℃和高于32 ℃的温度条件下对枇杷黄毛虫世代存活率不利[12]，7月和8月干旱对黄毛虫的发生有抑制作用[6]，也印证了这个可能性；第4代成虫诱集量普遍下降较快，可能与第3代幼虫发生期正值高温天气、成虫在高温下产卵少、卵孵化率下降使幼虫虫量不高有关。
图 1 2008-2013年枇杷黄毛虫灯下成虫种群数量特征Figure 1. Quantitative characteristics of adult Melanographia flexilineata during 2008-20132.1.2 田间幼虫种群动态特点
2008-2013年田间幼虫系统调查（图 2）结果表明：年度间枇杷黄毛虫均以第2代幼虫危害最重，发生期为6月中下旬至7月上旬，高峰期主要出现在6月下旬，少数年份持续至7月上旬危害重，第2代幼虫期正值枇杷成熟阶段，除叶、芽外，果实也受害，对枇杷生产影响最大，2008-2013年高峰期平均虫量140.0~188.9头·株-1·旬-1，占全年总量的32.57%~47.23%。第1代幼虫发生期为5月上旬至6月上旬，高峰期主要出现在5月中旬，少数年份推迟至5月下旬，峰期平均虫量30.6~95.0头·株-1·旬-1，以2008年最高，2013年虫量最低。第3代幼虫发生于7月下旬至8月中旬，此时正值1 a中气温最高、雨水相对较少的季节，特别是2013年为持续高温干旱的年份，成虫产卵少，卵存活率低，发生期总虫量及高峰期虫量均下降，年度间发生量差别较大，2008年高峰期虫量49.0头·株-1·旬-1，高于其他年份，而2013年由于持续高温干旱，虫量明显低于其他年份。第4代幼虫田间查见幼虫较少，与前几代次比虫量下降较快，9月中旬后，幼虫陆续化蛹越冬。
图 2 2008-2013年枇杷黄毛虫幼虫种群数量特征Figure 2. Quantitative characteristics of larvae Melanographia flexilineata during 2008-20132.2 发生期、发生量的预测
2.2.1 发生期、发生量预测模型的建立
枇杷黄毛虫主要以幼虫啃食幼芽、叶片、嫩茎表皮和果实等，对产量和品质影响大。为开展枇杷黄毛虫的适期防治，加强危害风险监测，以第1代至第3代幼虫发生高峰期（y1,y2,y3）和高峰期发生量（y4,y5,y6）为预报对象，初选前期灯下诱蛾量、田间幼虫密度等虫情基数因子及上年12月至当年1月、4月至9月每旬平均气温、平均相对湿度、降水量等相关性较大的气象因子共56个指标作为预报因子，采用逐步回归分析，建立预测预报模型。结果表明：枇杷黄毛虫幼虫发生高峰期、高峰期虫量与前期虫口基数、气象因子密切相关，各有15个因子入选预测模型（表 1和表 2）。发生期预测模型表明：前期田间幼虫虫口密度、降水量对发生高峰期影响最大，各有5个因子入选预测模型，前期灯下诱蛾量、平均相对湿度、气温等也有一定的影响。发生量预测模型表明：枇杷黄毛虫第1代至第3代幼虫旬高峰期虫量与前期田间幼虫密度、灯下诱蛾量呈显著正相关，气象因子中平均相对湿度和降水量分别有5个和4个因子入选预测模型，气温因子中5月下旬、6月上旬、7月上旬平均气温也作为预测因子入选发生量预测模型。
表 1 枇杷黄毛虫发生期预测因子历史数据Table 1. Predictive factors for occurrence period of Melanographia flexilineata年 份 x2 x5 x9 x10 x12 x20 x21 x22 x28 x31 x38 x40 x43 x46 x53 2008 15 77 91.7 15.8 27.5 28 8 3 9.4 66.6 157 22 508 84.7 60.2 2009 20 61.3 36.3 19.2 50.7 9 11 48 31.4 35.8 148 14 380 74.2 0 2010 2 69.5 37.1 12.7 70.6 61 40 9 24.9 24.7 43 15 5 84.2 173.9 2011 5 60.6 68 16.9 16.3 36 25 8 33.8 29.4 165.9 14 289 75 42 2012 19 64.8 148.5 17.8 30.1 25.5 10 9 36 105.3 188.9 9 210 83.3 37.8 2013 3 69.3 23 18.6 12 24.5 8.1 6 53.9 9.9 106.5 5 29 79.7 0 x2，x43分别为4月下旬、7月中旬诱蛾量（头·灯-1·旬-1）；x20，x21，x22，x38，x40分别为5月下旬、6月上旬、6月中旬、6月下旬、7月下旬幼虫量（头·株-1·旬-1）；x5，x46分别为上年12月、当年6月下旬相对湿度（%）；x10为4月中旬气温（℃）；x9，x12，x28，x31，x53分别为1月、4月中旬、5月中旬、5月下旬、7月中旬降水量（mm）。 表 2 枇杷黄毛虫发生量预测因子历史数据Table 2. Predictive factors for occurrence quantity of Melanographia flexilineata年 份 x5 x9 x11 x14 x17 x20 x27 x29 x31 x32 x40 x42 x48 x50 x53 2008 77 91.7 78.7 57.1 64.6 28 54.4 25.3 66.6 24.2 22 0 30.8 26.9 60.2 2009 61.3 36.3 75.8 53.2 49.8 9 60.6 22.6 35.8 24.8 14 328 29.1 58.1 0 2010 69.5 37.1 80.3 60.4 61.3 61 76 22.3 24.7 21.8 15 0 28.9 38.4 173.9 2011 60.6 68 60.1 57.9 63.7 36 52.5 20.8 29.4 24.2 14 0 30.8 2.4 42 2012 64.8 148.5 72.7 66.4 66.3 25.5 70.6 21.1 105.3 24.1 9 0 31.7 15.1 37.8 2013 69.3 23 49 63.7 70.4 24.5 70.1 25.4 9.9 22.9 5 0 31.6 1.5 0 x5，x11，x14，x17，x27分别为上年12月、当年4月中旬、4月下旬、5月上旬、5月中旬相对湿度（%）；x9，x31，x50，x53分别为1月、5月下旬、7月上旬、7月中旬降水量（mm）；x29，x32，x48分别为5月下旬、6月上旬、7月上旬气温（℃）；x20，x40分别为5月下旬、7月下旬幼虫量（头·株-1·旬－1）；x42为7月上旬诱蛾量（头·灯-1·旬－1）。 表 3 枇杷黄毛虫发生期、发生量预测模型Table 3. Forecasting model for occurrence quantity and period of Melanographia flexilineata模型编号 预测模型 R2 1 y1=19.146+0.007x2-0.067 x5 -0.034 x9-0.486 x10+0.028 x12 1.000** 2 y2=-90.562+0.190 x5+0.142 x9+1.119 x11-0.044 x14+0.749 x17 1.000** 3 y3=3.005+0.245 x20+0.005 x21+0.085 x22-0.086 x28+0.010 x31 1.000** 4 y4=389.386+0.238 x20-1.230 x27-7.494 x29 +0.539 x31-0.763 x32 1.000** 5 y5= 34.876-0.139 x38-0.028 x40-0.004 x43-0.048 x46-0.043 x53 1.000** 6 y6=-492.071+2.252 x40+0.010 x42+15.380 x48-0.525 x50+0.061 x53 1.000** y1为第1代幼虫发生高峰期（5月15日为1）；y2为第2代幼虫发生高峰期（6月21日为1）；y3为第3代幼虫发生高峰期（8月1日为1）；y4为第1代幼虫高峰期发生量（头）；y5为第2代幼虫高峰期发生量（头）；y6为第3代幼虫高峰期发生量（头）。 2.2.2 发生期、发生量预测模型的拟合率检验
将2013年的虫情、气象数据，应用1.4节的应验程度判定模式，分别对枇杷黄毛虫第1代至第3代幼虫的发生高峰期、高峰期发生量预测结果进行检验（表 4和表 5）。结果表明：用逐步回归法拟合的枇杷黄毛虫幼虫发生期、发生量模型，评分分值高，均大于99分，在准确的范围内。根据模拟结果，第1代、第2代和第3代枇杷黄毛虫幼虫发生高峰期模型拟合值分别为5月19日、6月25日、8月16日，实际分别为5月19日、6月25日、8月16日；第1代、第2代和第3代幼虫高峰期发生量模型拟合值分别为30.62，106.28和5.83头·株-1·旬-1，实际分别为31，107和6头·株-1·旬-1。模型拟合值接近实测值，拟合结果准确。
表 4 枇杷黄毛虫发生期预测模型验证Table 4. Verification of forecastin model for occurrence period of Melanographia flexilineata预测对象 a1 a2 t δ Ts 结论 第1代幼虫 5 5.04 10 2.11 99.95 准确 第2代幼虫 5 5.02 5 3.77 99.99 准确 第3代幼虫 16 15.99 15 5.81 100 准确 表 5 枇杷黄毛虫发生量预测模型验证Table 5. Verification of forecastin model for occurrence quantity of Melanographia flexilineata预测对象 a a1 a2 t δ Ts 结论 第1代幼虫 53.75 31 30.62 10 15.03 100 准确 第2代幼虫 152.55 107 106.28 5 26.27 99.98 准确 第3代幼虫 22.67 6 5.83 15 13.22 99.97 准确 3. 小结与讨论
3.1 枇杷黄毛虫种群特征及发生规律
本研究对枇杷黄毛虫进行了灯下诱蛾及田间幼虫系统调查。结果表明：枇杷黄毛虫在浙江余杭塘栖田间虫量与越冬期温湿度条件、主害代前期虫口基数及气候条件等密切相关，灯下诱蛾年度间总虫量2008年和2011年最高，2013年明显低于历年同期，与当年夏季持续高温干旱、不利于枇杷黄毛虫生长发育有关。田间幼虫量也呈现出相似的趋性，以2011年最高，2013年最低。关于枇杷黄毛虫在中国南方枇杷产区的每年发生代数、种群动态等，国内相关文献主要为福建[3, 5]、浙江[6-7, 12]等地的调查研究，在福建福州1 a发生5代[5]，上海地区1 a发生3代[13]。本研究表明枇杷黄毛虫在余杭塘栖1 a发生4代，与前人报道一致[6-7]，越冬代成虫始见期2008-2013年在4月中旬至下旬（4月14日至4月29日），2008年、2009年、2013年较早，2010-2012年相对较迟，与王恩等[6]的研究结果较一致（4月20日），比何富泉等[7]报道的（5月上旬）约早15~20 d。塘栖枇杷黄毛虫田间幼虫量年度间均以第2代危害最重，发生高峰期主要出现在6月下旬，少数年份推迟至7月上旬，主害代灯下成虫峰期时间、蛾量年度间有差异，总体均以第3代峰期诱蛾量最高，其次为第2代，第4代灯下成虫与田间幼虫量均下降较快，9月中旬后，幼虫陆续在树皮裂缝中、分枝处或树干基部和附近灌木上吐丝、结茧、化蛹越冬。
3.2 种群动态预测模型的影响因素
本研究对枇杷黄毛虫发生的相关因素进行了调查和统计分析，筛选出发生期和发生量各15个因子，建立了种群动态预测模型，为预测各代次枇杷黄毛虫的发生高峰、开展适期防治提供了依据。枇杷黄毛虫发生高峰期受前期田间幼虫虫口密度、降水量影响最大；田间幼虫高峰期发生量与前期虫口基数呈极显著正相关，是枇杷黄毛虫发生的最重要因素，平均相对湿度、降水量、气温等气象因子也有一定的相关性，春夏气候较为暖湿的年份，成虫产卵量多，第1代和第2代幼虫田间发生量大，气温过高则卵孵化率、幼虫存活率下降[12]，特别是2013年夏季高温干旱对黄毛虫的发生有抑制作用，而比较暖湿的季节有利于黄毛虫发生。本研究引入了多年来枇杷黄毛虫发生期每旬平均气温、相对湿度、降水量等气象因子及前期虫口基数进行建模，应用唐启义等[11]提出的应验程度判定模式进行拟合率检验，评分分值均在90以上，模型拟合值与实测值接近，拟合结果准确，对今后枇杷黄毛虫的适时、准确预报具有一定的实用价值，可供杭州地区及气象等环境条件总体较为相似的地区应用。
3.3 预测预报技术的完善与发展
建立枇杷黄毛虫种群动态模型的目的是为了预测预报更为适时准确、更有效地为生产服务，以减少因该虫危害造成的损失。枇杷黄毛虫年度间发生期和发生量受当年气候条件等因素影响，在每年发生盛期内又呈季节性消长，采用多年数据建立种群动态模型可实行较准确的预测。本模型主要依据2008年以来的数据资料，为更充分地显现其年度间的周期性规律，增强预报结果的准确性，笔者将逐年追加监测资料，使样本总数增多，校正模型，使预测精度不断提高，更好地服务于生产实际。
-
图 1 行走滑轮作用下承载索弯曲应力
Figure 1 Bending stress of bearing cable under the action of traveling pulley
图 4 等效应力云图(A)和等效位移云图(B)
Figure 4 Equivalent stress nephogram and equivalent displacement nephogram
图 5 不同跑车总质量下承载索的垂直位移时间历程
Figure 5 Vertical displacement time history of the load-bearing cable under different total mass of sports cars
图 6 交互作用对综合评分指标(y3)的影响
Figure 6 Influence of interaction on comprehensive evaluation index (y3)
表 1 升降卷扬机参数
Table 1. Lifting winch parameters
额定起
重/kg电机功
率/W额定电
压/V起升速度/
(m·min−1)卷起长
度/m挂钩质
量/kg整机质
量/kg200~400 1600 220 13~20 1~50 16.7 29 
下载: 导出CSV
表 2 HT250材料特性
Table 2. HT250 material properties
材料 泊松比 密度/
(kg·m−3)抗拉强度/
Pa弹性模量/
PaHT250 0.3 7.15×103 2.5×108 9×1010
下载: 导出CSV
表 3 试验因素编码
Table 3. Test factor coding
编码 因素 运输速度(x1)/
(m·s−1)跑车与竹材总质量(x2)/
kg倾斜角(x3)/
(°)−1 2 200 16 0 3 400 20 1 4 600 24
下载: 导出CSV
表 4 试验结果
Table 4. Test result
序号 运输速度
(x1)/(m·s−1)总质量
(x2)/kg倾斜角
(x3)/(°)最大位移
(y1)/m效率
(y2)/(kg·s−1)综合评分
(y3)序号 运输速度
(x1)/(m·s−1)总质量
(x2)/kg倾斜角
(x3)/(°)最大位移
(y1)/m效率
(y2)/(kg·s−1)综合评分
(y3)1 2 400 24 2.47 1.59 8.92 10 4 400 16 2.89 3.10 11.34 2 3 600 16 2.97 3.54 12.16 11 3 400 20 3.01 2.36 9.30 3 3 400 20 3.05 2.36 9.20 12 2 600 20 2.98 2.38 9.18 4 3 200 16 1.37 1.18 10.88 13 2 400 16 2.17 1.59 9.72 5 3 600 24 4.05 3.54 9.30 14 3 400 20 3.11 2.36 9.04 6 3 200 24 1.42 1.18 10.75 15 3 400 20 2.95 2.36 9.46 7 4 600 20 4.39 4.65 11.00 16 3 400 20 2.92 2.36 9.54 8 2 200 20 0.98 0.79 11.00 17 4 400 24 3.77 3.10 9.01 9 4 200 20 1.65 1.55 10.99
下载: 导出CSV
表 5 综合评分指标(y3)方差分析
Table 5. Analysis of variance of comprehensive evaluation index (y3)
变异来源 平方和 自由度 均方 F P 变异来源 平方和 自由度 均方 F P 模型 16.180 0 9 1.800 0 25.150 0 0.000 2** x12 0.046 2 1 0.046 2 0.646 3 0.447 6 x1 1.550 0 1 1.550 0 21.670 0 0.002 3** x22 5.370 0 1 5.370 0 75.180 0 ＜0.000 1** x2 0.490 0 1 0.490 0 6.860 0 0.034 5* x32 0.471 8 1 0.471 8 6.600 0 0.031 1* x3 4.680 0 1 4.680 0 65.500 0 ＜0.000 1** 残差 0.500 4 7 0.071 5 x1x2 0.837 2 1 0.837 2 11.710 0 0.011 1* 失拟 0.339 9 3 0.113 3 2.820 0 0.170 8 x1x3 0.585 2 1 0.585 2 8.190 0 0.024 3* 误差 0.160 5 4 0.040 1 x2x3 1.860 0 1 1.860 0 26.070 0 0.001 4** 总和 16.680 0 16 说明：**表示极显著(P＜0.01)；*表示显著(P＜0.05)。
下载: 导出CSV
-
[1] QI Shuhua, SONG Bing, LIU Chong, et al. Bamboo forest mapping in China using the dense Landsat 8 Image Archive and Google Earth Engine [J/OL]. Remote Sensing, 2022, 14: 762[2023-01-09]. doi: 10.3390/rs14030762. [2] 刘延鹤, 周建波, 傅万四, 等. 中国现代竹林场全程机械化经营模式探析[J]. 林业机械与木工设备, 2020, 48(11): 13 − 17. LIU Yanhe, ZHOU Jianbo, FU Wansi, et al. Analysis on the whole-process mechanized management mode of modern bamboo forest farms in China [J]. Forestry Machinery and Woodworking Equipment, 2020, 48(11): 13 − 17. [3] 费叶琦, 刘英, 陈甲伟, 等. 竹材采运机械发展探究[J]. 世界林业研究, 2021, 34(1): 31 − 36. FEI Yeqi, LIU Ying, CHEN Jiawei, et al. Development of bamboo cutting and transporting machinery [J]. World Forestry Research, 2021, 34(1): 31 − 36. [4] 吴立国, 郭克君, 苗振坤, 等. 索道集材跑车液压系统的设计与仿真[J]. 林业机械与木工设备, 2016, 44(11): 14 − 17, 23. WU Liguo, GUO Kejun, MIAO Zhenkun, et al. Design and simulation of the hydraulic system for a ropeway skidding carriage [J]. Forestry Machinery and Woodworking Equipment, 2016, 44(11): 14 − 17, 23. [5] 冯建祥, 罗才英. 悬索工程[M]. 厦门: 厦门大学出版社, 2010: 20. FENG Jianxiang, LUO Caiying. Suspension Engineering [M]. Xiamen: Xiamen University Press, 2010: 20. [6] 徐建超, 谭超, 赵亚东. 采煤机行走轮的参数建模与有限元分析[J]. 煤矿机械, 2011, 32(5): 91 − 92. XU Jianchao, TAN Chao, ZHAO Yadong. Parameters modeling and FEA of shearer wheel [J]. Coal Mine Machinery, 2011, 32(5): 91 − 92. [7] 亢方亮, 盛选禹. 基于ABAQUS的反应堆主泵抗震强度分析[J]. 机电工程技术, 2010, 39(5): 64 − 67, 113. KANG Fangliang, SHENG Xuanyu. Aseismatic strength analysis of main nuclear reactor pump based on ABAQUS [J]. Mechanical &Electrical Engineering Technology, 2010, 39(5): 64 − 67, 113. [8] 周成, 邵跃林. 基于ADAMS与ABAQUS的刚柔耦合动力学分析方法[J]. 机械制造与自动化, 2014, 43(1): 131 − 133. ZHOU Cheng, SHAO Yuelin. Dynamic analysis method of flexible multi-body system based on ADAMS and ABAQUS [J]. Machine Building &Automation, 2014, 43(1): 131 − 133. [9] 钟一. 山地运输索道多因素耦合的动力学特性研究[D]. 长沙: 中南林业科技大学, 2021. ZHONG Yi. Study on Dynamic Characteristics of Multi-factor Coupling of Mountain Transportation Ropeway [D]. Changsha: Central South University of Forestry and Technology, 2021. [10] KNAWA M, BRYJA D. Modeling problems of steeply inclined cableway subjected to moving load [J]. Proceedings in Applied Mathematics and Mechanics, 2009, 9(1): 263 − 264. [11] 舒彩霞, 曹士川, 廖宜涛, 等. 基于ADAMS的油菜割晒机顺向侧铺装置参数优化与试验[J]. 农业机械学报, 2022, 53(增刊2): 11 − 19, 38. SHU Caixia, CAO Shichuan, LIAO Yitao, et al. Parameter optimization and experiment of forward laying device for rape windrower based on ADAMS [J]. Transactions of the Chinese Society for Agricultural Machinery, 2022, 53(suppl 2): 11 − 19, 38. [12] 苑玉凤. 多指标正交试验分析[J]. 湖北汽车工业学院学报, 2005, 19(4): 53 − 56. YUAN Yufeng. Analysis of multi-target orthogonal experiment [J]. Journal of Hubei University of Automotive Technology, 2005, 19(4): 53 − 56. [13] 王相友, 张蒙, 李学强, 等. 辊式导流马铃薯定重装袋机设计与试验[J]. 农业机械学报, 2021, 52(6): 159 − 168. WANG Xiangyou, ZHANG Meng, LI Xueqiang, et al. Design and test of roller-guided potato fixed weight bagging machine [J]. Transactions of the Chinese Society for Agricultural Machinery, 2021, 52(6): 159 − 168. 期刊类型引用(2)
1. 吴燕君，洪文英，章忠梅，吴耀，缪强. 设施黄瓜白粉病流行动态与预测模型. 浙江农业学报. 2022(01): 104-111 . 
百度学术2. 郑建洪. 4种常用药剂对枇杷黄毛虫的田间防效研究. 现代农业科技. 2017(21): 116-117 . 百度学术其他类型引用(0)
-
-
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20230141
点击查看大图
计量
- 文章访问数: 692
- HTML全文浏览量: 111
- PDF下载量: 42
- 被引次数: 2
