Optimization of rice straw disassembly process parameters and fiber self-interweaving structure
-
摘要:
大尺寸高长径比纤维是影响非织造秸秆纤维生态毯质量的关键。机械揉搓是获得高长径比纤维的最佳处理技术工艺。该研究探究了揉丝机主轴转速和锤齿间隙对不同含水率秸秆纤维拆解揉丝的影响,其工艺参数为主轴转速、锤齿间隙及物料含水率,指标输出为秸秆丝化率与标定单位生产率,并对纤维自交织结构特性进行评价。结果表明,采用Box-Behnken响应面法(response surface method, RSM)开发了一个可确定输入和输出参数间的函数,2种评价指标模型均具有较高的可信度,决定系数分别为R2>0.95,R2>0.84。优化数学模型工艺参数:含水率41.2%,主轴转速2 498 r/min,锤齿间隙12.84 mm组合下秸秆丝化率达 96.93%,标定单位功率生产率为29.33 kg/(kW·h)。大尺寸高长径比秸秆纤维长度>70 mm和长宽比>80比例增幅分别达120.55%和16.01%,有效改善纤维自交织能力,所开发的RSM模型秸秆的实际揉搓拆解丝化率达96.27%,加工成本降低27.50元/t,可广泛应用于生产大尺寸纤维秸秆基高值农用制品。该研究为生物质秸秆的有效循环利用和高值化材料产业的绿色可持续发展提供重要依据和技术支撑。
Abstract:Rice straw is an important renewable resource. Under the circumstance of global resource shortage and energy shortage, it is of great significance to vigorously develop comprehensive utilization of straw. It is not only an important way to realize agricultural circular economy and low-carbon economy, but also involves sustainable development issues such as soil fertility, environmental security and efficient utilization of renewable resources in the entire agricultural ecosystem. Straw fiber ecological blanket has a good function of air permeability and water permeability, and has a significant positive effect on the ecological protection of soil surface microenvironment. Mechanical kneading is the key technology for the dismantling and optimization of agricultural waste straw fiber. In this study, the optimal process input is the spindle speed, hammer tooth clearance and material water content, and the output is the straw silk rate and calibrated unit productivity. An intelligent model was developed using Box-Behnken response surface method (RSM) to determine the functional relationship between input and output parameters. The three evaluation index models all have high reliability, and the determination coefficients are R2 > 0.95 and R2 > 0.84, respectively. The variance analysis of RSM regression model clarified that the influence of various factors on straw silk ratio was in the order of spindle speed > water content > hammer tooth clearance, and the linear relationship between spindle speed and straw silk ratio was increasing. The significant influence on the calibrated unit power productivity of the kneading machine is in the order of hammer tooth clearance > spindle speed > moisture content. When the hammer tooth clearance is constant, it decreases first and then increases with the spindle speed and straw moisture content. The results of factor interaction analysis showed that when the water content of straw was 40%, the spindle speed was 2 500 r/min, and the hammer tooth clearance was 13 mm, the silk ratio of rice straw was 96.75%. When the straw moisture content is 60%, the spindle speed is 2 300 r/min, and the hammer tooth clearance is 11.08 mm, the maximum calibrated unit power productivity is 29.96 kg/(kW·h). Mathematical model process parameters optimization: straw moisture content is 41.2%, spindle speed is
2498 r/min, hammer tooth clearance is 12.84 mm. Under the combination of these parameters, the rice straw silk ratio is 96.93%, and the calibrated unit power productivity is 29.33 kg/(kW·h). Under optimized conditions, the proportion of fibers with length of 40-70 mm and > 70 mm showed an increasing trend, especially the fiber with length of > 70 mm increased by 120.55%. The proportion of fiber length-to-width ratio > 50-80 and > 80 increased from 17.39% and 6.52% to 38.64% and 16.01%, respectively. Ultra depth of field microscopy showed that the self-interweaving formability of kneaded straw fiber was significantly improved. The optimized mechanical kneading cost of straw fiber can be reduced by 27.50 yuan/t. The mechanical kneading technology strategy of the RSM model developed in this paper is effective for improving the decomposition rate of biomass straw and reducing the consumption of biomass straw, and can be widely used in the production of high-value agricultural products based on large fiber bundles of straw. The mechanical kneading technique strategy of the RSM model developed in this paper is effective for improving the decomposing silk rate and reducing the consumption of biomass straw, and can be widely used in the production of large-size fiber bundle straw based high-value agricultural products. -
0. 引 言
东北地区苏打盐碱地是中国盐碱地的主要集中区,对盐碱土壤进行改良研究与实践,使之适宜农作物的生长发育,对提高粮食安全具有重要意义[1-2]。针对东北苏打盐碱地特征,课题组设计了一种具有叶脉状内通道的深松铲[3],可以实现自润减阻配合改良液层施,进而达到改良东北地区苏打盐碱土的目的。目前多孔分布管广泛应用于给排水、滴喷灌等领域[4]。为研究其出流特性,国内外众多学者开展了大量分析研究[5]。不同学者针对不同情况下对多孔管内流体分布的计算方程进行推导, RAVIKUMAR等[6]建立了分析方程来确定压力的变化系数,推导出计算流量变化系数的简单解析方程,用于求解数值问题中的流量变异系数;CHELLAM等[7]研究了具有1道和2道可渗透壁的均匀多孔管内不可压缩层流的特性,建立了四阶微分方程并得到了所有壁面雷诺数的解;刘文华等[8]结合水力学分析理论,对坡度、压力水头、孔距、管长4个参数对沿程压力分布的影响规律开展研究;杜涛等[9]建立多孔管出流控制体模型并对达西-韦斯巴赫公式进行完善,得出了多孔管沿程水头损失的计算公式;一些学者通过物理试验方法对多孔管的液体出流进行了相应研究,王亚朦等[10]根据达西-韦斯巴赫公式,结合理论分析与试验的方法,对不同坡度下多孔分布管内流体压力分布规律开展研究;鞠学良等[11]根据微灌毛管水力解析模型,提出了考虑适宜布置形式的均匀坡微灌毛管设计方案以提高灌水均匀度;杨宝中等[12]结合单孔出流试验结果,对多孔管土中出流计算公式进行分析研究,为多孔管地下灌溉提供参考;随着计算流体动力学的发展,数值模拟技术已成为研究单相流和结构优化的重要工具,钱卫忠等[13]使用有限元分析的方法分析多孔管孔径、管长、水压等因素对多孔管出流速度和流量的影响,并得到影响多孔管出流速度及出流量影响显著的因素;王鹏辉[14]提出了基于CFD的“流场分析设计方法”并对喷淋装置几何结构进行了设计优化;韩圆圆[15]通过数值模拟的方法分析了T型三通管内的气液流动特性及相分布规律。
综上,关于多孔管出流速度及出流量的研究多与水利灌溉工程相关,目前未见将其与深松部件结合的相关研究。针对叶脉状内通道的流场分布规律及影响因素不明确的问题,本文拟采用基于计算流体动力学的仿真方法结合物理试验对其进行研究。首先利用仿真的方法,开展单因素试验,优选出影响内通道分支出口流速的因素水平范围值,其次开展Box-Behnken试验,建立各影响因素与分支出口流速之间的二阶回归模型,进而得到最优参数组合。按最优参数组合条件对叶脉状内通道进行设计加工后进行相应的物理试验,得出物理试验值。通过仿真与物理试验结果的对比,检验叶脉状内通道结构优化设计的效果,从而达到对叶脉状内通道优化设计的目的,为利用叶脉状内通道多分支出口管结构在东北地区苏打盐碱地进行深松及化学改良耦合技术提供理论基础。
1. 数值计算模型与计算方法
1.1 数值计算模型
以课题组自研的叶脉状内通道自润滑减阻深松铲的内部结构——叶脉状内通道多分支出口管为研究对象。其深松铲结构由铲柄、弧形铲刃、凿型铲尖及叶脉状多分支出口内通道组成。深松铲铲体长度为1200 mm,其中铲柄长500 mm,弧形铲刃长600 mm,铲尖长100 mm,铲刃横截面为等腰三角形且刃口角度为42°。叶脉状内通道与深松铲铲柄及弧形铲刃的总长度尺寸一致,由主通道直管部分及弧形叶脉状多分支出口管两部分构成,其总长度为1100 mm,主通道直径10 mm,直管部分长度500 mm,弧形叶脉状多分支出口管部分是半径550 mm的圆弧,其长度为600 mm。其中每组分支出口由两个在主通道X轴线上呈42°夹角对称分布的分支出口组成,同时每组分支出口在Y轴上以相同间距在叶脉状内通道多分支出口弧形管段上均匀分布。其叶脉状内通道深松铲及内通道多分支出口管结构如图1。
图 1 叶脉状内通道深松铲与叶脉状内通道多分支出口管结构图1.主通道 2.深松铲体 3.深松铲刃口夹角42° 4.直管部分 5.弧形管部分 6.末端分支出口 7.分支出口夹角42° 8.首端分支出口 9.主通道入口Figure 1. Structure diagram of deep loosening shovel and a porous pipe flow channel model in leaf vein shaped inner channel1.Main channel 2.Deep loosening shovel body 3. Deep loosening shovel blade angle 42° 4. Straight pipe section 5. Curved tube section 6. End branch outlet 7. Branch outlet angle 42° 8. The head branch outlet 9. Main channel entrance1.2 数值模型网格划分
数值模拟计算中选取的叶脉状内通道多孔管的流体计算流域为主通道入口至叶脉状多分支出口处。采用ICEM中的O-Block技术进行流体计算域网格划分[16],可获得较规整的结构网格和较好的边界层网格。计算域采用适应度较强的非结构四面体网格进行划分,算法采用补丁适形法,采用美国机械工程协会(American Society of Mechanical Engineers,ASME)推荐的网格收敛指数 GCI 进行网格离散误差的估计,对网格进行无关性验证[17-19]。最终确定选取中等网格,此时的模型网格总数为198353个,节点总数为43799个,并检查网格问题与质量,结果显示平均网格质量、偏度和纵横比分别为0.837、0.228、1.85, 所生成网格的质量较高,符合本文计算要求。
1.3 数值计算方法及试验设计
1.3.1 数值计算方法
图2为多孔管计算域网格。液体入口位于直管部分顶端,出口位于弧形管部分上的各分支出口上,其余位置均密闭,壁面处设置为无滑移壁面,计算流体使用不可压缩的液态水,进口边界选择流速入口边界条件,对于出口边界条件,采用压力出口边界条件压力值为0,参考压力为标准大气压[20]。采用k-ε湍流模型,依靠计算流体力学求解软件Fluent进行求解[21-24]。在多孔管设置中给定主通道入口流速,经过软件模拟可得到多孔管在该主通道入口流速下达到稳定状态时的内通道各分支出口流速。
1.3.2 试验设计
文章中所讨论的叶脉状内通道分支多孔管,目的是在深松过程中,利用分支出口流出的液体改良剂在铲刃面上形成液膜达到自润滑减阻目的,同时利用流出的改良剂对盐碱地进行化学改良。这就需要每组分支出口的流速要达到一定要求,同时每组分支出口的流量也要达到一定要求,通过大量田间试验发现,叶脉状内通道的结构参数(孔径、孔间距、孔数)及运行参数(主通道入口流速)对分支出口的出流速度及液膜的形成起到决定性作用。因此本文选择分支出口的孔径、孔间距、孔数及主通道入口流速为研究因素。其中孔径即为分支出口的直径,孔间距即为相邻两组分支出口的距离,孔数即为分支出口的组数,主通道入口流速即为改良剂流入主通道的速度。
结合实际工况及多次预试验发现,分支出口6 mm以下的孔径所流出的液体在深松过程中不能够形成将整个深松铲铲面完全覆盖的液膜,即不能保证流量要求,因此选择6 mm以上的孔径进行讨论。分支出口孔间距距离太小及孔数过少则达不到分层减阻及化学改良液层施时所需效果,同时分支出口孔间距距离太大及孔数过多则使最末端分支出口流量太小。在确保各组分支出口的表面积流速大于1 m/s以及所流出的液体能够使铲刃面形成液膜的前提下,初步确定叶脉状内通道分支出口孔径、孔间距、孔数和主通道入口流速4个因素的水平取值范围,进而开展单因素仿真模拟试验,对因素水平进行初步优选,以确定各因素最佳工作水平范围,具体因素水平见表1。研究单因素对各组分支出口流速的影响时,其他因素取中间水平值。
表 1 叶脉状内通道多孔管数值模拟参数Table 1. Numerical simulation parameters of leaf vein shaped inner channel porous tube水平
Factors分支出口孔径
Branch outlet
aperture D/mm分支出口孔间距
Branch outlet
hole spacing
A/mm孔数
Hole number
B主通道入口流速
Main channel
inlet velocity
C/(m·s−1)1 6 80 3 3 2 7 90 4 4 3 8 100 5 5 4 9 110 6 6 5 10 120 7 7 内通道分支出口截面上的最大流速即为表面积最大流速,内通道分支出口截面上的平均流速即为平均流速,两种流速表达的含义不同。距主通道入口最近的一组分支出口即为首端出口,距主通道入口最远的一组分支出口即为末端出口;主通道入口流速在沿程上经过首端或末端分支出口的分流而减少的值即为首端或末端出口分流损失;主通道入口流速在沿程上经过各分支出口的分流而减少的值与主通道入口流速的比值即为流速的分流损失率,以此来量化表征主通道入口流速经过各分支出口被分流的多少,进而分析不同影响因素水平对流速的影响大小。即:首端出口分流损失等于主通道入口流速减去首端分支出口流速;末端出口分流损失等于主通道入口流速减去末端分支出口流速;首末端出口分流损失等于首端分支出口流速减去末端分支出口流速;首端出口分流损失率等于首端出口分流损失与主通道入口流速的比值;末端出口分流损失率等于末端出口分流损失与主通道入口流速的比值;首末端出口分流损失率等于首末端出口分流损失与首端分支出口流速的比值。
2. 结果与分析
2.1 叶脉状内通道多分支出口管沿程流速变化规律
图3为叶脉状内通道主通道流速剖面图。由图3发现叶脉状主通道上沿程流速变化趋势为先逐渐增大后逐渐减小。以叶脉状内通道多分支出口管主通道直径10 mm、内通道分支出口孔径6 mm、孔间距100 mm、孔数5、主通道入口端流速5 m/s时为例,叶脉状主通道上最大流速值为5.736 m/s,最小流速为0。沿程流速从入口处开始逐渐增加直至首组分支出口处上端达到沿程流速最大值,经过首端分支出口下端后开始逐渐减小,在达到末端分支出口后沿程流速逐渐减小至最小值。出现主通道内部沿程流速先逐渐增大后逐渐减小的原因,可能是在射流和惯性重力及离心力的影响作用下,在主通道内部产生沿程流速加速现象,使主通道内部沿程流速先逐渐增大到流速最大值[25];随后液体经过首组内通道分支出口后压力减小,液体开始逐级分流,沿程流速开始逐渐减小。
2.2 单因素试验结果
通过对初步确定的影响叶脉状内通道分支出口流速的因素进行单因素试验,根据试验数据分析各因素对分支出口流速的影响情况,进一步确定参数范围,试验结果如图4所示。
图 4 各因素下的内通道分支出口平均流速注:B1为首端分支出口;B2为第二分支出口;B3为第三分支出口;B4为第四分支出口;B6为第六分支出口;B5和B7为末端分支出口。Figure 4. Average flow velocity at the outlet of internal channel branches under various factorsNote:B1is the head branch outlet; B2 is the second branch outlet; B3 is the third branch outlet; B4 is the fourth branch outlet; B6 is the sixth branch outlet; B5 and B7 are the end branch outlet.2.2.1 分支出口孔径对流速的影响
根据图4a可知,随着内通道分支出口孔径的增加,同一分支出口下的流速逐渐减小,同一孔径下的各分支出口流速逐渐减小。当孔径由6 mm增加到10 mm时,表面积最大流速的首端出口分流损失率、末端出口分流损失率及首末端出口分流损失率均逐渐增加,当孔径为6 mm时,表面积最大流速的三种分流损失率分别为4.43%、38.52%和35.67%,均为最小。同时随着孔径的增加内通道分支出口平均流速的首端出口分流损失率、末端出口分流损失率及首末端出口分流损失率同样均逐渐增加,当孔径为6 mm时,平均流速的三种分流损失率分别为67.14%、58.32%和−26.83%,均为最小。孔径由6 mm增加到10 mm时,内通道分支出口表面积最大流速首端出口减小1.567 m/s,末端出口减小1.878 m/s;平均流速首端出口减小0.657 m/s,末端出口减小1.475 m/s。根据结果数据发现在主通道入口流速统一的情况下,内通道分支出口孔径越小,相邻两分支出口的出流速度差值变化越小,随着分支出口孔径的增加,同一孔径下相邻两分支出口出流速度差值逐渐增加,进而使各个分支出口平均流速减小。在孔径为6 mm时液体出流速度最快,因此液体流量也最大,能够达到自润滑减阻的目的。因此优选出叶脉状内通道多孔管模型内通道分支出口孔径为6 mm。
2.2.2 分支出口孔间距对流速的影响
根据图4b可知在内通道分支出口孔间距增加的情况下,内通道分支出口首端出口流速逐渐增大,末端出口流速逐渐减小,这是因为分支出口流速随着孔间距的增加在第二、三分支出口之后降幅开始加剧。同一孔间距下各分支出口流速整体上逐渐减小,只有当孔间距为80 mm时各分支出口流速先增加后减小。当孔间距为80 mm时,首端出口分流损失率同其他孔间距相比最大为26.30%;末端出口分流损失率及首末端出口分流损失率同其他孔间距相比最小,分别为60.22%和46.02%。当孔间距为120 mm时,首端出口分流损失率同比最小为15.46%;末端出口分流损失率及首末端出口分流损失率同其他孔间距相比最大,分别为68.58%和62.83%。同时随着孔间距的增加内通道分支出口平均流速的首端出口分流损失率、末端出口分流损失率及首末端出口分流损失率变化同表面积最大流速变化趋势一样。孔间距由80 mm增加到120 mm时,内通道分支出口表面积最大流速首端出口增加0.542 m/s,末端出口减小0.418 m/s;平均流速首端出口增加0.268 m/s,末端出口减小0.285 m/s。在主通道入口流速统一的情况下,内通道分支出口孔间距越小对各分支出口的首端出口出流速度影响变化较大,反而对末端出口出流速度影响变化较小,当孔间距由90 mm增加到110 mm时流速变化不大,这是叶脉状内通道分支出口这一结构造成的,但在对因素水平讨论优选的过程中是有必要的,因为在实际工作环境中分支出口的位置关系到改良剂液体对深松铲减阻的部位,以及利用叶脉状内通道这一结构进行化学改良过程中,达到化学改良剂分层施用这一目的。因此考虑综合实际试验因素,优选出叶脉状内通道多孔管模型内通道分支出口孔间距90、100、110 mm进行Box-Behnken仿真试验,进一步探究孔间距及其它因素水平交互下对分支出口流速的影响。
2.2.3 分支出口孔数对流速的影响
根据图4c可知在内通道分支出口孔数增加的情况下,内通道分支出口首末端流速均逐渐减小,同一孔数下各分支出口流速整体上是逐渐减小的,只有当孔数为3时各分支出口流速逐渐增加。孔数由3增加到7时,内通道分支出口表面积最大流速的首端出口分流损失率、末端出口分流损失率及首末端出口分流损失率均先增加后减小,当孔数为3时,三种分流损失率均为最小,分别为10.23%、21.49%和12.54%。当孔数为7时,此时首端出口分流损失率不是最大,分别为24.20%、67.44%和57.05%。同时随着孔数的增加内通道分支出口平均流速的首端出口分流损失率、末端出口分流损失率及首末端出口分流损失率变化同表面积最大流速变化趋势一样。孔数由3增加到7时,内通道分支出口表面积最大流速首端出口减少0.698 m/s,末端出口减小2.298 m/s;平均流速首端出口增加0.455 m/s,末端出口减小1.663 m/s。在主通道入口流速统一的情况下,内通道分支出口孔数越少对各分支孔的首末两端组出流速度影响变化越小,通过试验数据可以得出孔数在4~6时对流速影响变化幅度适中,同时为了考虑综合实际试验因素,优选出叶脉状内通道多孔管模型内通道分支出口孔数4、5、6进行Box-Behnken仿真试验,进一步探究孔数对分支出口流速的影响。
2.2.4 主通道入口流速对流速的影响
根据图4 d可知在内通道主通道入口流速增加的情况下,内通道各分支出口的流速是逐渐增大的,同一主通道入口流速下各分支出口流速先增加后逐渐减小。主通道入口流速由3 m/s增加到7 m/s时,内通道分支出口表面积最大流速的首端出口分流损失率逐渐增加,末端出口分流损失率逐渐减少,当主通道入口流速为3 m/s时,首端出口分流损失率同比最小为14.83%;末端出口分流损失率及首末端出口分流损失率同比最大为82.74%和79.74%。当主通道入口流速为7 m/s时,首端出口分流损失率同比最大为26.03%;末端出口分流损失率及首末端出口分流损失率同比最小为54.53%和38.53%。同时随着主通道入口流速的增加内通道分支出口平均流速的首端出口分流损失率先减少后增加,末端出口分流损失率及首末端出口分流损失率均减少。当主通道入口流速为7 m/s时,首端出口分流损失率为78.58%,末端出口分流损失率及首末端出口分流损失率74.38%和−19.58%,均为最小。在其他因素不变的情况下,主通道入口流速为5、6、7 m/s时首端出口分流损失率、末端出口分流损失率及首末端出口分流损失率均为中间值,因此优选出叶脉状内通道多孔管模型内通道分支出口主通道入口流速为5、6、7 m/s进行Box-Behnken仿真试验,进一步探究主通道入口流速对分支出口流速的影响。
3. Box-Behnken仿真试验
基于单因素仿真试验中优选出的各因素范围,设计二次回归正交试验,建立叶脉状内通道分支出口平均流速回归模型,以得到最优参数组合[26]。内通道分支出口孔径在单因素试验中已经优选出最优值为6 mm,其他因素水平如表2所示,Box-Behnken试验设计及结果如表3所示。
表 2 因素水平表Table 2. The factor level table水平
Level分支出口孔间距
Branch outlet
hole spacing
A/mm分支出口孔数
Holes number of branch
outlet B主通道入口流速
Main channel
inlet velocity
C/(m·s−1)−1 90 4 5 0 100 5 6 1 110 6 7 表 3 Box-Behnken试验设计及结果Table 3. Box Behnken experimental design and results序号No. X1 X2 X3 出口流速Outlet velocity/(m·s−1) 1 −1 −1 0 6.069 2 1 −1 0 5.896 3 −1 1 0 5.483 4 1 1 0 5.889 5 −1 0 −1 4.713 6 1 0 −1 4.961 7 −1 0 1 6.456 8 1 0 1 6.806 9 0 −1 −1 4.903 10 0 1 −1 4.769 11 0 −1 1 7.263 12 0 1 1 6.334 13 0 0 0 5.621 14 0 0 0 5.685 15 0 0 0 5.679 16 0 0 0 5.558 17 0 0 0 5.681 注:X1、X2、X3分别为分支出口孔间距、分支出口孔数、主通道入口流速的水平值。 Note: X1, X2, X3 are the levels of branch outlet hole spacing, branch outlet hole number, main channel inlet velocity respectively. Box-Behnken试验方差分析结果如表4所示,由表4的分析结果可知,X2、X3和X2X3对平均流速影响极其显著, X1、X1X2和X 22对平均流速影响显著;X1X3、X12和X32对平均流速影响不显著。模型的P<0.0001,说明回归模型高度显著;模型失拟项P>0.05,说明模型失拟性不显著,回归模型拟合程度高。回归方程决定系数R2=0.9911,校正决定系数R2=0.9796,与1非常接近,且变异系数CV=1.75%。根据统计学上F值越大、表明该因素对结果的影响越大[27-28],通过F值大小,可以判定各因素对平均流速影响的重要性,试验因素对平均流速影响从大到小为主通道入口流速、内通道分支出口孔数、内通道分支出口孔间距;综上表示该回归模型极其显著,能够可靠和真实的反映真实情况,可用于进一步的目标平均流速的预测分析。
表 4 Box-Behnken试验方差分析Table 4. Variance analysis of Box Behnken experiment方差源
Variance
source均方
Mean
square自由度
Degree of
freedom平方和
Sum of
squareF值
F-valueP值
P-value模型Model 7.830 9 0.870 86.21 <0.0001** X1 0.086 1 0.086 8.55 0.0222* X2 0.343 1 0.343 33.96 0.0006** X3 7.060 1 7.060 698.99 <0.0001** X1X2 0.084 1 0.084 8.30 0.0236* X1X3 0.003 1 0.003 0.26 0.6273 X2X3 0.158 1 0.158 15.65 0.0055** X1² 0.012 1 0.012 1.18 0.3141 X2² 0.078 1 0.078 7.75 0.0271* X3² 0.006 1 0.006 0.54 0.4849 残差
Residual0.071 7 0.010 失拟项
Lack of fit0.059 3 0.020 6.39 0.0525 纯误差
Pure error0.012 4 0.003 总和Sum 7.900 16 注: **表明影响极显著(P<0.01),*表明影响显著(P<0.05)。 Note: ** indicates that the impact is extremely significant (P<0.01), and * indicates that the impact is significant (P<0.05). 通过Design–Expert 11.0对Box-Behnken试验结果进行多元回归拟合,得到仿真试验出口流速二阶回归方程:
V=5.64+0.104X1−0.207X2+0.939X3+0.145X1X2+0.026X1X3−0.199X2X3+0.053X21+0.136X22+0.036X23 (1) 在Design–Expert 11.0软件的优化模块中,以分支出口平均流速最大值为优化目标值,将二阶回归方程式(1)中的不显著项剔除后进行优化求解,求解出最优参数组合为内通道分支出口孔间距110 mm、内通道分支出口孔数4、主通道入口流速7 m/s。
为验证仿真优化参数的可靠性,在最优参数条件下进行3次Fluent仿真模拟,得到分支出口平均流速为3.262、3.255、3.275 m/s,其平均值为3.264 m/s。
4. 试验验证
根据数值模拟选定的最优参数组合,加工出用于验证试验的叶脉状内通道多孔管,主通道入口端通过软管、阀门、流量计、变径接头与增压泵相连接构成试验通路,各分支出口由软管及量杯进行承接构成泄放通路,主通道入口流速通过控制流量计算获得,液体流动示意图及实际装置见图5。
图 5 叶脉状内通道多分支出口流动示意图及试验装置图1.叶脉状内通道深松铲 2.深松铲卡钳 3.流量计 4.变径接头 5.水箱 6.机架 7.软管 8.量杯 9.分支出口Figure 5. Schematic diagram and experimental setup diagram of multi branch outlet flow in leaf vein like internal channels1.Vein shaped inner channel deep loosening shovel 2. Deep loosening shovel caliper 3. Flowmeter 4. Reducer union 5. Water tank 6. Frame 7. Hose 8. Counting cup 9. Branch outlet因无法直接测定各分支出口流速,利用图5中量杯,通过软管与深松铲分支出口进行密闭连接,当深松铲达到稳定工作水平状态时,测量出单位时间内各个分支出口所流出液体量Q,根据液体流量计算式(2)计算各分支出口的平均流速[29-31]。
Q=SV (2) 式中Q为分支出口流量,L/min;S为分支出口面积,cm2;V为分支出口平均流速,m/s。
在最优参数组合下,经过3次重复物理试验计算出内通道分支出口平均流速值为2.988、3.102、2.824 m/s,其平均值为2.971 m/s。与仿真值平均误差为8.97%,验证了数值模拟的可靠性,同时检验了叶脉状内通道结构优化设计的效果。分析各分支出口的平均流速小于数值模拟,造成该数据异常的原因可能是实际多分支开孔处不平整,增加了叶脉状内通道多分支出口处的压力损失;另外,实际管路开孔直径的精度有限,而分支孔管出流量对孔径的变化较为敏感,孔径的变化都会对管的流量造成一定的影响。
5. 结 论
1)通过数值仿真分析得到叶脉状内通道多孔管主通道上的沿程流速变化规律为先增大后减小。沿程流速从入口处开始逐渐增加直至首端分支出口处上端达到沿程流速最大值,经过首端分支出口下端后开始逐渐减小,在达到最末组出口端后沿程流速逐渐减小至最小值。
2)采用单因素仿真分析对4个因素进行初步优选,以流速最大为前提,优选出内通道分支出口孔径为6 mm、孔间距为90、100、110 mm、孔数为4、5、6、主通道入口流速为5、6、7 m/s。
3)利用Design–Expert 11.0对仿真数据进行分析,求解得到最优参数组合为内通道分支出口孔间距为110 mm、孔数为4、主通道入口流速为7 m/s。对流速影响因素大小依次为主通道入口流速、内通道分支出口孔数、内通道分支出口孔间距。
4)在最优参数组合下进行仿真与物理试验,内通道分支出口平均流速分别为3.264 m/s和2.971 m/s,误差为8.97%。
-
表 1 试验因素及水平
Table 1 Factors and level of test
水平
Level因素 Factor 含水率
Moisture content X1/%主轴转速
Spindle speed X2/(r·min−1)锤齿间隙
Clearance between hammer
and gear X3/mm-1 40 2 300 10 0 50 2 400 13 1 60 2 500 16 表 2 Box-Behnken试验方案与结果
Table 2 Test plan and result of Box-Behnken
试验号
No.含水率
Moisture content
X1/%主轴转速
Spindle speed
X2/(r·min−1)锤齿间隙
Clearance between hammer and
gear X3/mm丝化率
Filamentation rate φ1/%标定单位
功率生产率
Calibration of
unit power
productivity φ2/(kg·(kW·h)−1)1 −1 0 0 96.23 25.00 2 1 0 0 94.68 25.53 3 0 −1 0 93.79 28.21 4 0 1 0 96.42 27.69 5 0 0 −1 93.28 26.95 6 0 0 1 95.95 25.28 7 −1 −1 0 94.16 27.07 8 −1 1 0 96.75 29.61 9 1 −1 0 92.18 30.15 10 1 1 0 95.79 28.57 11 −1 0 −1 95.28 27.31 12 −1 0 1 96.16 25.42 13 1 0 −1 91.26 26.59 14 1 0 1 93.52 25.17 15 0 −1 −1 90.28 27.61 16 0 −1 1 93.87 25.79 17 0 1 −1 94.25 28.08 表 3 回归模型方差分析
Table 3 Variance analysis regression model
指标
index方差来源
Source of variance平方和
Sum of squares自由度
Degree of freedom均方
Mean squareF P 丝化率
Filamentation rate φ1/%模型 52.95 9 5.88 14.95 0.0009 **X1 12.43 1 12.43 31.58 0.0008 **X2 16.05 1 16.05 40.78 0.0004 **X3 8.71 1 8.71 22.12 0.0022 **X1X2 0.2601 1 0.2601 0.6608 0.4431 X1X3 0.4761 1 0.4761 1.21 0.3078 X2X3 0.8657 1 0.8657 2.2 0.1817 X12 0.5203 1 0.5203 1.32 0.2880 X22 2.19 1 2.19 5.55 0.0496 *X32 5.86 1 5.86 14.87 0.0062 **残差 2.7555 7 0.3936 失拟 0.64 3 0.33 0.39 0.8426 误差 2.115 4 0.46 总和 55.71 16 标定单位功率生产率
Calibration of unit power
productivity φ2/
(kg·(kW·h)−1)模型 33 9 3.67 4.09 0.0384 *X1 0.2579 1 0.2579 0.2876 0.6084 X2 0.5671 1 0.5671 0.6323 0.4526 X3 4.16 1 4.16 4.63 0.0684 X1X2 4.24 1 4.24 4.72 0.0663 X1X3 0.057 1 0.057 0.0635 0.8083 X2X3 0.2522 1 0.2522 0.2812 0.6123 X12 1.25 1 1.25 1.4 0.2755 X22 13.06 1 13.06 14.56 0.0066 **X32 0.0241 1 0.0241 0.0269 0.8745 残差 6.28 7 0.8969 失拟 3.125 4 2.872 0.46 0.2572 误差 2.56 4 3.16 总和 39.28 16 注:P<0.01(极显著**);0.01≤P<0.05(显著*)。 Note: P<0.01 (highly significant **); 0.01≤P<0.05 (significant *). 表 4 机械揉搓1 000 kg秸秆加工成本分析
Table 4 Analysis of processing cost of mechanical kneading 1 000 kg straw
项目
Item丝化率
Filamentation rate/%电费
Electricity consumption
/(元·t−1)原料损耗率
Loss rate/%原料费
Raw material cost
/(元·t−1)水费
Water cost
/(元·t−1)人工费
Labor cost/(元·d−1)加工成本费用
Processing cost expense
/(元·t−1)优化工艺
Optimization process96.7 7.81 3.3 310 0.26 16.67 334.74 原有工艺
Original process90.0 9.50 10.0 333 0.51 19.23 362.24 -
[1] YANG J, CHANG J, CAO X, et al. Energy, environmental, and socioeconomic potential benefit of straw resources utilization under water availability limit condition[J]. Journal of Cleaner Production, 2023, 392,15: 136274.
[2] MOUSAVI-AVVAL S, SAHOO K, NEPAL P, et al. Environmental impacts and techno-economic assessments of biobased products: A review[J]. Renewable and Sustainable Energy Reviews, 2023, 180: 113302. doi: 10.1016/j.rser.2023.113302
[3] SIN X, XIE M, MAI L, et al. Biobased plastic: A plausible solution toward carbon neutrality in plastic industry[J]. Journal of Hazardous Materials, 2022, 435(5): 129037.
[4] WU Z, PENG K, ZHANG Y, et al. Lignocellulose dissociation with biological pretreatment towards the biochemical platform: A review[J]. Materials Today Bio, 2022, 16: 100445. doi: 10.1016/j.mtbio.2022.100445
[5] SUN E, ZHANG Y, YONG C, et al. Biological fermentation pretreatment accelerated the depolymerization of straw fiber and its mechanical properties as raw material for mulch film[J]. Journal of Cleaner Production, 2021, 284(15): 124688.
[6] 范肖东,韩善明,曲萍,等. 氯化胆碱/尿素预处理水稻秸秆及其解离纤维结构性能的表征[J]. 林产化学与工业,2023,43(4):9-16. FAN Xiaodong, HAN Shanming, QU Ping, et al. Characterization of choline chloride/urea pretreated rice straw and its dissociative fiber structure and properties[J]. Chemistry and Industry of Forest Products, 2023, 43(4): 9-16. (in Chinese with English abstract)
[7] 孙恩惠,黄红英,武国峰,等. 稻壳/大豆蛋白基黏合剂成型育苗钵性能评价及成因分析[J]. 农业环境科学学报,2015,34(6):1202-1209. SUN Enhui, HUANG Hongying, WU Guofeng, et al. Performance evaluation and analysis of nursery containers made from rice husk/soybean protein adhesives[J]. Journal of Agro-Environment Science, 2015, 34(6): 1202-1209. (in Chinese with English abstract)
[8] 曲萍,黄红英,孙恩惠,等. 一种高长径比秸秆纤维的制备方法及其在针刺成型秸秆毯中的应用:202211300323.3[P]. (2022-10-24). [9] ZHANG Y, WANG J, GONG S, et al. Straw mulching improves soil water content, increases flag leaf photosynthetic parameters and maintaines the yield of winter wheat with different irrigation amounts[J]. Agricultural Water Management, 2021, 249(30): 106809.
[10] 彭鵾,靳红梅,孙恩惠,等. 秸秆纤维毯覆盖对葡萄园土壤理化性质、杂草防控及果实品质的影响[J]. 农业资源与环境学报,2024,DOI: 10.13254/j.jare.202.0005. PENG Kun, JIN Hongmei, SUN Enhui, et al. Impact of straw fiber blanket mulching on soil properties, weed management, and grape quality in vineyards[J]. Journal of Agricultural Resources and Environment, 2024, DOI: 10.13254/j.jare.202.0005. (in Chinese with English abstract)
[11] ZHAO W, YU R, DONG W, et al. The influence of long carbon fiber and its orientation on the properties of three-dimensional needle-punched CF/PEEK composites[J]. Composites Science and Technology, 2021, 203(8): 108565.
[12] MENG X, FAN W, MAHARI W, et al. Production of three-dimensional fiber needle-punching composites from denim waste for utilization as furniture materials[J]. Journal of Cleaner Production, 2021, 281(25): 125321.
[13] SINGH C, MUKHOPADHYAY S, RENGASAMY R. A sustainable approach to oil spill cleanup by kapok and waste cotton needle punched nonwoven blends[J]. Industrial Crops and Products, 2023, 191: 115939. doi: 10.1016/j.indcrop.2022.115939
[14] 刘禹辰. 基于离散元法的玉米秸秆揉丝仿真与试验研究[D]. 兰州:甘肃农业大学,2023. LIU Yuchen. Simulation and Experimental Study of Corn Stalk Kneading Based on Discrete Element Method[D]. Lanzhou: Gansu Agricultural University, 2023. (in Chinese with English abstract)
[15] 张馨怡. 玉米秸秆滚压揉搓破碎机理与机构仿真研究[D]. 哈尔滨:哈尔滨工业大学,2020. ZHANG Xinyi. Study on Mechanism and Mechanism Simulation of Rolling, Rubbing and Crushing of Corn Stalks[D]. Harbin: Harbin Institute of Technology, 2020. (in Chinese with English abstract)
[16] 成翔. 秸秆饲料收获机捡拾揉碎装置设计与实验研究[D]. 呼和浩特:内蒙古农业大学,2023. CHENG Xiang. Design and Experiment of Picking and Crushing Device in Straw Feed Harvester[D]. Hohhot: Inner Mongolia Agricultural University, 2023. (in Chinese with English abstract)
[17] 李秀清. 基于离散元方法的锤片粉碎机玉米粉碎过程模拟及锤片性能试验研究[D]. 呼和浩特:内蒙古农业大学,2020. LI Xiuqing. The Simulation of Maize Crushing by Hammer Mill and the Test of Hammer Performance Based on Discrete Element Method[D]. Hohhot: Inner Mongolia Agricultural University, 2020. (in Chinese with English abstract)
[18] 曹丽英,杨左文,焦巍,等. 锤片式粉碎机噪声源识别及降噪方法[J]. 农业工程学报,2018,34(7):59-65. CAO Liying, YANG Zuowen, JIAO Wei, et al. Noise sources identification and noise reduction methods of hammer mill[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(7): 59-65. (in Chinese with English abstract)
[19] 翟之平,张龙,刘长增,等. 秸秆抛送装置外壳振动辐射噪声数值模拟与试验验证[J]. 农业工程学报,2017,33(16):71-79. ZHAI Zhiping, ZHANG Long, LIU Changzeng, et al. Numerical simulation and experimental validation of radiation noise from vibrating shell of stalk impeller blower[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(16): 71-79. (in Chinese with English abstract)
[20] 方梅,郁志宏,张文杰,等. 盘刀式铡草机粉碎物料运动过程分析与试验[J]. 农业工程学报,2021,37(7):76-84. FANG Mei, YU Zhihong, ZHANG Wenjie, et al. Analysis and experiments of the movement process for the shredded material of disc knife chaff cutter[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(7): 76-84. (in Chinese with English abstract)
[21] 钱义,王迪,张珏,等. 粉碎机异形筛片气流场数值模拟及试验研究[J]. 中国农业大学学报,2020,25(3):79-87. QIAN Yi, WANG Di, ZHANG Jue, et al. Numerical simulation and experimental research on the gas flow field of special-shaped sieve sheet of crusher[J]. Journal of China Agricultural University, 2020, 25(3): 79-87. (in Chinese with English abstract)
[22] 王海庆,玛纳拉苏仁,田海清,等. 锤片式粉碎机三角形筛片的设计及试验研究[J]. 农机化研究,2018,40(7):72-76. WANG Haiqing, MANALASUREN, TIAN Haiqing, et al. The design and experimental study on the triangular screen of hammer mill[J]. Journal of Agricultural Mechanization Research, 2018, 40(7): 72-76. (in Chinese with English abstract)
[23] 侯泽华,崔红梅,闫鹏,等. 9RS-60型揉碎机工作性能试验研究[J]. 农机化研究,2020,42(1):141-146. HOU Zehua, CUI Hongmei, YAN Peng, et al. The orthogonal test study on working performance of 9R-60 crushed mill[J]. Journal of Agricultural Mechanization Research, 2020, 42(1): 141-146. (in Chinese with English abstract)
[24] 李松毅. 玉米秸秆揉碎装置设计与性能试验研究[D]. 呼和浩特:内蒙古农业大学,2019. LI Songyi. Design and Performance Test of Corn Straw Crushing Device[D]. Hohhot: Inner Mongolia Agricultural University, 2019. (in Chinese with English abstract)
[25] 王铁军,王铁良,崔红光,等. 喂入调节式秸秆破包揉丝机设计与试验[J]. 农业机械学报,2021,52(6):148-158. WANG Tiejun, WANG Tieliang, CUI Hongguang, et al. Design and experiment of adjustable feeding straw bale-breaking and rubbing filament machine[J]. Transactions of the Chinese Society for Agricultural Machinery, 2021, 52(6): 148-158. (in Chinese with English abstract)
[26] MATSUMOTO J, OHSAKI S, NAKAMURA H, et al. Quantitative analysis of wet kneading states by a novel compression test[J]. Advanced Powder Technology, 2024, 35(2): 104336. doi: 10.1016/j.apt.2024.104336
[27] LIANG R, ZHANG B, ZHOU P, et al. Power consumption analysis of the multi-edge toothed device for shredding the residual film and impurity mixture[J]. Computers and Electronics in Agriculture, 2022, 196: 106898. doi: 10.1016/j.compag.2022.106898
[28] 麻乾,刘飞,赵满全,等. 揉碎机揉碎机理分析及锤片结构优化[J]. 农业工程学报,2016,32(S2):7-15. MA Qian, LIU Fei, ZHAO Manquan, et al. Working mechanism and structure optimization of hammer of rubbing machine[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(S2): 7-15. (in Chinese with English abstract)
[29] 李洪成. 锤片式粉碎机工作性能影响机理与优化研究[D]. 武汉:华中农业大学,2023. LI Hongcheng. Study on Influence Mechanism and Optimization of the Hammer Mill Grinding Performance[D]. Wuhan: Huazhong Agricultural University, 2023. (in Chinese with English abstract)
[30] 陈涛,衣淑娟,李衣菲,等. 铡切揉碎协同式牧草揉丝机设计与试验[J]. 农业机械学报,2024,55(2):149-159. CHEN Tao, YI Shujuan, LI Yifei, et al. Design and test of cutting and crushing cooperative silk kneading machine[J]. Transactions of the Chinese Society for Agricultural Machinery, 2024, 55(2): 149-159. (in Chinese with English abstract)
[31] CHEN X, JIANG H, WANG G, et al. Disposable bamboo fiber meal boxes characterized by efficient preparation, excellent performance, and the potential for beneficial degradation[J]. Journal of Cleaner Production, 2024, 434: 139973. doi: 10.1016/j.jclepro.2023.139973
-
期刊类型引用(1)
1. 郭佩杰,陈雨,陈宇翔,汪海英,肖爱萍,陈军,钟培林. 液态粪肥施肥机肥罐防浪板的设计与试验. 农业工程学报. 2025(02): 47-57 . 本站查看
其他类型引用(0)