高级检索

深部煤炭流态化开采装备自主行走机构多缸推进同步控制

鲍久圣, 李玥锋, 周恒, 阴妍, 赵少迪, 王忠宾, 葛世荣

鲍久圣,李玥锋,周 恒,等. 深部煤炭流态化开采装备自主行走机构多缸推进同步控制[J]. 煤炭科学技术,2024,52(2):264−278

. DOI: 10.12438/cst.2023-1513
引用本文:

鲍久圣,李玥锋,周 恒,等. 深部煤炭流态化开采装备自主行走机构多缸推进同步控制[J]. 煤炭科学技术,2024,52(2):264−278

. DOI: 10.12438/cst.2023-1513

BAO Jiusheng,LI Yuefeng,ZHOU Heng,et al. Synchronization control of multi-cylinder propulsion for autonomous walking of deep underground fluidized coal mining equipments[J]. Coal Science and Technology,2024,52(2):264−278

. DOI: 10.12438/cst.2023-1513
Citation:

BAO Jiusheng,LI Yuefeng,ZHOU Heng,et al. Synchronization control of multi-cylinder propulsion for autonomous walking of deep underground fluidized coal mining equipments[J]. Coal Science and Technology,2024,52(2):264−278

. DOI: 10.12438/cst.2023-1513

深部煤炭流态化开采装备自主行走机构多缸推进同步控制

基金项目: 

中国矿业大学“双一流”学科拓新平台建设专项资助项目(2019XKPT03);江苏高校优势学科建设工程资助项目(PAPD)

详细信息
    作者简介:

    鲍久圣: (1979—),男,安徽桐城人,教授,博士生导师,博士。E-mail:cumtbjs@cumt.edu.cn

  • 中图分类号: TD42

Synchronization control of multi-cylinder propulsion for autonomous walking of deep underground fluidized coal mining equipments

Funds: 

Special Funding Project for the Construction of "Double First-Class" Discipline Innovation Platform of China University of Mining and Technology (2019XKPT03); Priority Discipline Construction Project of Jiangsu Higher Education Institutions (PAPD)

  • 摘要:

    我国浅部煤炭正逐渐开采殆尽,向地球深部开发资源已成为必然趋势和国家需求。针对2 000 m以深的深部煤炭开采难题,设计了一种适用于深部煤炭流态化开采装备的自主行走机构,并重点对自主行走机构中液压推进系统的多缸同步控制难题进行了研究。首先,基于流态化开采工艺原理及装备组成,设计了一种增阻迈步式自主行走机构,可实现采掘、转化和输出等多舱体的分段式自主行走;其次,针对液压推进系统的多缸同步控制要求,分析对比了主从控制、相邻交叉耦合控制、偏差耦合控制、均值耦合控制4种控制策略以及比例、积分、微分控制(PID)算法、自抗扰控制器(ADRC)的优缺点,分别开展了在均匀负载、突变负载、时变负载3种工况下的控制性能仿真试验;再次,采用雷达图测评法对不同控制策略下的同步控制性能进行综合评价,最终选定基于ADRC的均值耦合控制方法为最佳同步控制策略;最后,研制自主行走机构试验台并开展了多液压缸同步控制试验,试验结果表明:当采用基于ADRC的均值耦合同步控制策略时,4个液压缸在不同工况下的最大同步误差均能保持在±5 mm以内,且具有优异的鲁棒性,可满足流态化开采装备自主行走机构的多缸同步推进要求。

    Abstract:

    The shallow coal in China is gradually being exhausted, and mining resources from deep parts of the Earth has become an inevitable trend and national demand. Faced with the difficulty of deep coal mining beyond 2000 meters, an autonomous walking mechanism suitable for fluidized mining equipment was designed, and the focus was on the problem of multi-cylinder synchronous control in the hydraulic propulsion system of this mechanism. Firstly, based on the principle of fluidized mining technology and equipment composition, a resistance increasing step type autonomous walking mechanism was designed, which can achieve segmented autonomous walking of multiple compartments such as excavation, transformation, and output. Secondly, in response to the requirements of multi-cylinder synchronous control in hydraulic propulsion systems, four control strategies, namely master-slave control, adjacent cross-coupling control, deviation coupling control, and mean coupling control, as well as the advantages and disadvantages of proportional-integral-derivative control(PID) algorithm and automatic disturbance rejection controller (ADRC), were analyzed and compared. Control performance simulation tests were conducted under three working conditions: uniform load, sudden load, and time-varying load. Once again, the radar chart evaluation method was used to comprehensively evaluate the synchronization control performance under different control strategies, and the mean coupled control method based on ADRC was selected as the optimal synchronization control strategy. Finally, an autonomous walking mechanism test bench was developed and multiple hydraulic cylinder synchronization control experiments were conducted. The test results showed that when using the mean coupled synchronization control strategy based on ADRC, the maximum synchronization error of the four hydraulic cylinders under different working conditions can be maintained within ± 5 mm, and they have excellent robustness. It can meet the requirements of multi-cylinder synchronous propulsion for the autonomous walking mechanism of fluidized mining equipment.

  • 2022年我国煤炭产量为45.6亿t,比2021年增长10.5%[1] ,仍占据较大比例。然而,经过几十年大规模开采,露天矿以及浅煤层的煤炭资源已经被逐步开采殆尽,浅部资源正趋于枯竭[23]。研究表明,在现有技术条件下井工开采的极限深度约为1 500 m[4],因此突破煤炭开采方式和开采极限深度,向深部要资源已成为必然选择和重大需求[5]。2016年以来,谢和平院士团队[67]创新提出了深地资源流态化开采科学构想,该技术以现有全断面硬岩掘进机技术为先导,采用多功能盾构式结构的设备,实现煤炭采选充、电气热一体化开采和原位转化。流态化开采设备由采掘舱、流态化转化舱、产能输出舱和充填舱等多部分组成[8],多舱协同工作是实现煤炭原位开采、分选、转化、填充和运输的技术前提。然而,由于流态化开采系统舱体繁多、体积庞大,如何解决多舱体装备系统在深地条件下的自主行走和协同推进是目前亟待解决的关键基础问题之一。

    深地流态化开采系统借鉴敞开式全断面硬岩隧道掘进机(TBM)掘进技术[9],采用多缸协同作业的方式推进,为了实现盾构式多舱开采装备在深部地质条件下的直线掘进这一目标,需对各舱体间的多个液压缸同步推进控制问题进行研究。目前,国内外学者在此领域已开展了一些研究,例如:YANG等[10]提出一种基于PID算法的压力和速度双闭环耦合控制的盾构机液压缸同步控制策略,通过试验表明两个液压缸同步误差保持在3 mm以内;薛召等[11]提出了一种双缸同步控制策略,采用PID结合模糊单神经元的方式,并进行仿真和试验,结果表明双缸同步误差在±3.3 mm以内;徐莉萍等[12]提出了一种偏差耦合同步控制结合PID算法的控制方法,并通过仿真进行验证,结果表明双缸之间的同步误差小于1.5 mm;XU等[13]提出了一种基于径向基函数(RBF)网络自适应滑模控制的主从同步控制策略,提高了三液压缸同步动作的效果;SANATI等[14]提出了一种负载力和同步误差反馈的交叉耦合同步方法,使整个系统的非线性及不确定因素得到了很好的抑制,同时提高了系统的控制精度;卢子帅等[15]将相邻交叉耦合与自抗扰控制相结合,仿真结果表明该策略同步控制效果良好且加快了系统的响应速度。然而,现有研究多集中于双缸同步控制,且多采用传统的主从控制和偏差耦合控制等控制方法,对于三缸及以上的多液压缸同步控制,其结构和控制系统都较为复杂,同步控制难度大,目前相关研究较少,针对深地复杂环境下的多缸同步控制研究则未开展。

    深地岩体在“三高一扰动”的影响下,会出现大变形、强流变等特征,并且在地球内动力的作用下不同尺度的岩体岩石在不同时间和不同空间发生稳态/非稳态的演化[16],形成复杂的地质格局,如大尺度的均匀特质的岩石层、大尺度岩土突变的岩石层、小尺度非均匀特质的地质层。这些复杂的地质特征会导致开采设备在不同位置受到不同的地质约束,使得多缸同步控制系统的不同部分受到不同的阻力,产生振动、位移突变等现象,进而影响整个设备的运动平衡。鉴于此,针对深部流态化开采装备多缸同步推进控制需求,分析了多种不同的多缸闭环同步控制方法的优缺点,提出了一种均值耦合同步控制方法结合自抗扰控制器的多缸同步控制策略,并通过仿真和试验验证了其同步性和稳定性。

    深部煤炭流态化开采系统的结构原理如图1所示,通过多个舱体的协同作业实现流态化开采与原位转化,舱体间采用液压缸实现推进功能,并在舱与舱之间设置有柔性伸缩保护罩,防止煤层塌方卡死设备。

    图  1  深地流态化开采系统示意
    Figure  1.  Schematic diagram of deep ground fluidized mining system

    针对深部流态化开采系统的盾构式多舱开采装备结构特点,设计了一种增阻迈步式自主行走机构[17]。由图1可知,最少3个舱体即可实现装备的自主行走功能,所以选取采掘舱、流态转化舱及产能输出舱3个具有特色的舱体作为设计研究对象,如图2所示。自主行走机构是由采掘舱、转化舱、输出舱、前/后推进液压缸和增阻系统组成的,每个舱体之间采用液压缸铰接,通过液压缸之间的协同推进实现设备前行、转向等功能。借鉴敞开式TBM行进原理,在转化舱和输出舱上方布置有增阻液系统,采用增阻液压缸伸出顶板代替“撑靴”作用,实现舱体固定功能。

    图  2  增阻迈步式自主行走机构结构
    Figure  2.  Structural diagram of increased resistance walking equipment for autonomous walking

    自主行走机构的推进液压缸与增阻液压缸协同工作推动各舱增阻迈步式前进,整个迈步过程原理如图3所示,①转化舱上布置的增阻液压缸活塞顶出抵在煤层上,转化舱与煤层锁定;②转化舱上的4个前推进液压缸活塞顶出,推动采掘舱前行进行煤炭开采,实现头部采掘舱行进;③转化舱增阻液压缸活塞缩回,解除转化舱和煤层的锁定,尾部舱的增阻液压缸活塞顶出抵在煤层上,尾部舱与煤层锁定;④ 4个前推进液压缸活塞缩回,4个后推进液压缸活塞顶出,通过“前拉后顶”实现转化舱行进;⑤尾部舱增阻液压缸活塞缩回,解除尾部舱与煤层的锁定,转化舱上的增阻液压缸活塞再次顶出抵在煤层上,转化舱再次与煤层锁定;⑥后推进液压缸缩回,拉动尾部舱行进。

    图  3  增阻迈步式自主行走原理
    Figure  3.  Principle of increasing resistance and stepping autonomous walking

    为了解决深部煤炭流态化开采装备自主行走机构在推进过程中的多液压缸同步控制问题,对推进液压系统建模仿真,对比了所提出的控制策略与3种常用控制策略的综合性能,并对其鲁棒性进行优化。

    液压推进系统如图4所示,在采掘舱和转化舱之间布置4个液压缸作为前推进装置,在转化舱和输出舱布置4个液压缸为后推进装置。

    图  4  液压推进系统原理
    1—过滤器;2—恒压变量泵;3—安全阀;4—比例调速阀;5—三位四通换向阀;6—液压缸;7—压力传感器;8—位移传感器
    Figure  4.  Schematic diagram of hydraulic propulsion system

    由恒压变量泵2向前后推进系统的8个液压缸提供恒定的液压压力,采用比例调速阀4控制液压缸活塞伸出缩回的速度,通过换向阀5控制液压缸运动的方向。将位移传感器8测量的液压缸6伸缩时的位移信号进行微分,得到速度信号并将其反馈至控制器以实现液压系统的闭环控制。

    由于前、后推进系统的液压油路完全相同,因此选择前推进系统进行仿真,同时因为换向阀、过滤器等原件对液压系统影响较小,所以对前推进系统的建模进行一定的简化处理,所建立的模型如图5所示。

    图  5  推进系统仿真模型
    1—恒压源;2—调速阀;3—推进液压缸;4—速度传感器;5—模拟负载
    Figure  5.  Propulsion system simulation model

    当前多液压缸同步控制方法主要包括主从同步控制、相邻交叉耦合同步控制、偏差耦合同步控制等,并对以上3种控制方法及提出的均值耦合同步控制方法进行分析。

    1)主从同步控制。根据从动元件跟随不同主动元件耦合关系的不同,可将主从同步控制分成2种形式,如图6a所示为形式①,选择其中某个液压缸为主动元件,剩余液压缸为从动元件,从动液压缸跟随主动液压缸进行动作。如图6b所示为形式②,后一个液压缸跟随前一个液压缸的运动进行动作。主从控制通过在增加各控制通道之间的耦合关系实现液压缸同步控制,然而主从控制是从动液压缸跟随主动液压缸进行动作,因此存在从动液压缸跟随滞后时延的现象。

    图  6  主从同步控制
    Figure  6.  Master-slave synchronous control

    2)相邻交叉耦合同步控制。相邻交叉耦合液压缸同步控制是使相邻的两个液压缸之间的动作存在耦合关系,其结构如图7所示。图8为相邻交叉耦合误差补偿器结构,将相邻通道2的输出信号和相邻通道4的输出信号与通道1的输出信号进行比较,将差值送至补偿器并进行相加,作为反馈信号输入动作原件。然而非相邻之间的通道存在滞后现象,此外,相邻交叉耦合误差补偿器需要3个输入信号,对于具有n个液压缸的相邻交叉耦合同步控制系统,则需要3n个控制器。

    图  7  相邻交叉耦合同步控制
    Figure  7.  Adjacent cross-coupling synchronous control
    图  8  相邻交叉耦合误差补偿器结构
    Figure  8.  Adjacent cross-coupling error compensator structure

    3)偏差耦合同步控制。偏差耦合同步控制方法如图9所示,偏差耦合同步控制是在并行同步控制的基础上加入位置误差补偿器,如图10所示为偏差耦合误差补偿器结构,控制系统中每一个缸的位移信号均与通道1的输出信号进行比较,将比较后的差值传送至误差补偿器,然后将所有补偿器的结果累加,形成通道1液压缸的补偿器输出信号,反馈至控制器实现对通道1的液压缸的偏差耦合同步控制。

    图  9  偏差耦合同步控制
    Figure  9.  Deviation coupling synchronous control
    图  10  偏差耦合误差补偿器结构
    Figure  10.  Structure of deviation coupling error compensator

    误差补偿器考虑所有液压缸的运动状态,每一个液压缸位移信号都传送至误差补偿器,因此当一个通道的液压缸受到扰动导致活塞速度发生变化时,剩余液压缸立即进行跟随,该控制方法同步效果较好。然而偏差耦合系统需要n2个控制器。

    4)均值耦合控制方法。以上3种多缸同步闭环控制系统中,主从控制系统简单,然而存在从动件滞后的现象;相邻偏差耦合同步控制系统在相邻液压缸之间的同步效果较好,然而在非相邻的液压缸之间存在时延滞后现象;偏差耦合同步控制系统的液压缸同步效果较好,然而其存在运算量大,结构复杂,需要多个控制器等缺点。针对上述问题,提出了一种均值耦合同步控制方法,如图11所示。

    图  11  多缸均值耦合同步控制
    Figure  11.  Multi-cylinder mean coupled synchronous control

    均值耦合控制是一种用于实现多通道液压系统同步控制的方法。通过计算所有液压缸位移信号进行平均计算,求得平均值,将平均值作为每个通道的同步误差补偿器的输入信号。该方法综合分析所有液压缸的运行状况,将两舱体之间4个液压缸的移动量的平均值作为参考对象:

    $$ X\mathrm{_m}=\frac{1}{4}\sum\limits_{i=1}^4X_i $$ (1)

    式中:Xm为参考模型输出位移;Xi为第i个液压缸的输出位移。

    比例阀输入信号为

    $$ u_i=G\left[R\left(t\right)-V_i\right]+T\left(X_{\mathrm{m}}-X_i\right) $$ (2)

    式中:ui为输入信号;G为速度环控制器输出增益;R(t)为输入信号;Vi为液压缸的输出速度;T为位置环同步误差补偿器输出增益。

    从式(2)中可明显得出均值耦合同步控制算法整合了全部通道的液压缸的运行情况,每个通道中只需n个位置控制器和n个速度控制器,相较于上述2种控制方法,大幅减少了同步控制方法所需的控制器数量,同时均值耦合同步控制方法在具备偏差耦合同步控制高精度的基础上,又克服了相邻交叉耦合控制中同步时延滞后的难题。

    对主从控制、相邻交叉耦合控制、偏差耦合控制以及均值耦合控制4种同步控制方法进行仿真分析,对比其同步控制性能的优劣,将PID算法融入均值耦合控制方法中,得到如图12所示模型。

    图  12  多缸同步推进控制仿真模型
    1—恒压源;2—比例调速阀;3—推进液压缸;4—速度、位移传感器;5—负载模拟函数;6—位置耦合器;7—同步误差补偿器;8—控制器;9—阈值控制器;10—目标速度信号
    Figure  12.  Simulation model for multi cylinder synchronous propulsion control

    1)仿真工况设计。如图13所示,按照地质特点的不同,将地下煤层工况分成3类[1820],即流态化开采设备会遇到3类不同工况,第1类:均匀工况,设备在地质均质的工况下开采前行,所有液压缸承受相同负载,但随着设备前进速度加大,设备驱动力也随之变大。第2类:突变工况,设备在掘进过程中突然碰到形状不规则的硬岩,4个液压缸承受不同时刻的变负载。第3类:时变工况,设备在掘进过程中遇到复杂的地质特征,煤层中存在不同形貌、不同硬度、分布不均匀的硬岩,各个液压缸承受的外负载在时刻发生变化并在一定范围内波动,外部负载的时刻变化使液压阀的阀芯受到振动的影响,阀门开度难以控制,从而引起流量的波动,导致液压缸同步困难。

    图  13  不同地质工况
    Figure  13.  Different geological conditions

    因目前尚无2 000 m以深的煤层采掘数据实例,以锦屏水电站引水隧道工程作为依据[21],计算行走机构在不同工况下的驱动力,设置3种工况的加载曲线如图14所示。图14a为在稳定地质环境下推进速度与驱动力关系,拟定仿真工况0~5 s推进速度为1 mm/s,5~10 s推进速度均匀增加,10 s时推进速度达到2 mm/s。图14b为突变工况时的模拟负载曲线,前2 s内4个液压缸承受2 000 kN均匀外部载荷并保持不变,第2秒时采掘舱在液压缸推进下,刀盘的局部位置切割到复杂地质的硬岩,导致一个液压缸承受的载荷突变为3 000 kN,另外3个液压缸外部载荷保持不变,第4秒时随着推进系统继续向前推进,采掘舱刀盘切割硬岩的区域变大,第2个液压缸承受的外部载荷也突变到3 000 kN,随后的第6秒及第8秒,第3、第4油缸所支撑的刀盘区域承受的外部载荷突变为3 000 kN,并保持不变;图14c展示了随机选取4种在一定范围内波动的不同的白噪声,并乘上增益系数k形成的时变载荷,并将4种不同的时变负载分别导入4个液压缸的仿真模型,用以模拟设备的时变负载工况,时变负载模型如图15所示。

    图  14  不同地质工况加载曲线
    Figure  14.  Loading curves for different geological conditions
    图  15  时变负载模型
    Figure  15.  Time-varying load model

    同时拟用各缸的最大位移同步误差作为评判标准:

    $$ {e_{\max }} = \max \{ {X_i}\} - \min \{ {X_i}\} $$ (3)

    式中,emax为最大位移同步误差。

    2)均匀负载工况仿真。均匀负载工况下,由图16可以看出当主从控制在刚启动时,因该控制方法存在滞后性,最大位移同步误差呈上升趋势,1.5 s时趋向稳定,误差保持在3.0×10−5 m左右,然而5 s后随着目标速度的增加,主从控制的最大位移同步误差也开始逐渐上升,这是由于主从控制不存在耦合关系,当目标速度发生变化时,各个液压缸之间发生滞后现象从而产生误差。其他3种控制方法明显优于主从控制,并且在算法调整下,最大位移同步误差趋近于0 m。

    图  16  均匀负载工况液压缸同步仿真结果
    Figure  16.  Simulation results of hydraulic cylinder synchronization under uniform load conditions

    3)突变负载工况仿真。突变负载工况下,由图17可以看出主从控制的最大位移同步误差在2 、4、6、8 s时由于各液压缸承受的外负载突变而产生较大波动,在其余3种控制方法中,相邻交叉耦合同步控制下的最大位移同步误差大于偏差耦合控制和均值耦合控制,这是由于其非相邻的液压缸之间耦合存在滞后造成的,且3种控制方法最大同步误差均保持在2.5×10−5 m,具有较好的同步性能。

    图  17  突变负载工况液压缸同步仿真结果
    Figure  17.  Simulation results of hydraulic cylinder synchronization under sudden load conditions

    4)时变负载工况仿真。时变负载工况下,由图18可以看出主从控制方法的最大位移同步误差在白噪声的时变负载下出现无规则的发散现象,抗干扰能力较弱。另外3种同步控制方法最大位移同步误差趋于稳定保持在5×10−5 m以内,其中均值耦合控制方法最大同步误差保持在小范围内波动,且优于另外两种耦合控制方法,证明该同步控制方法具有更强的抗干扰能力,其同步性更加稳定。

    图  18  时变负载液压缸同步仿真结果
    Figure  18.  Synchronous simulation results of time-varying load hydraulic cylinder

    对控制方法的评估不能仅限于精度,而应全面考虑真实工况中的运算速度、鲁棒性和成本等因素。为了更准确评估不同控制方法性能的优劣,提出了一种综合考虑各因素的评估方法—雷达图面积测评法。评分标准见表1,该方法结合了4个指标,分别是精度、成本、速度和鲁棒性。这些指标通过最大位移同步误差、控制器数量、仿真时长及在时变负载下的仿真结果来体现。

    表  1  控制策略评分标准
    Table  1.  Control strategy evaluation criteria
    精度速度鲁棒性成本评分
    3
    2
    1
    下载: 导出CSV 
    | 显示表格

    绘制如图19所示雷达图,其中四边形的面积越大表明控制策略的综合性能越优异。雷达图中主从控制、相邻交叉耦合、偏差耦合、均值耦合面积依次为12、16、10、20。其中均值耦合方法的面积最大,各项性能指标较为突出,所以选择均值耦合控制方法应用到流态化开采装备推进系统中。然而,从雷达图中发现4种控制方法鲁棒性指标较差,这是因为经典PID算法抗干扰能力不强造成的,因此采用现代控制理论中的自抗扰控制器提升控制方法的鲁棒性性能指标。

    图  19  4种同步控制方法性能雷达图
    Figure  19.  Performance radar chart of four synchronous control methods

    前文同步仿真的控制器是基于PID算法设计的,而PID算法属于经典控制理论,是线性控制器,无法对系统内部的速度、加速度、非线性因素及系统外部扰动等进行调节控制,采用线性系统近似非线性系统导致控制了4种同步控制方法的精度下降。而自抗扰控制器(ADRC)属于现代控制理论的控制器,其能够估算和补偿液压同步控制系统中的非线性因素[22],有效应对系统的不确定性和变化,使得ADRC在面对参数变化、外部扰动和模型不确定性等情况下依然能够保持较好的控制性能,并且与传统的控制设计方法不同,自抗扰控制器不要求精确的数学模型,通过将复杂系统转化为简单的积分串联型系统,使得控制结构更简化,提高了控制器的可实现性和可调节性,综上自抗扰控制器拥有强大的鲁棒性、抗干扰能力和快速响应性能。

    首先对推进系统进行分析,确定系统阶数。线性化的调速阀流量方程为

    $$ {Q_{\mathrm{L}}} = {k_{\mathrm{q}}}{x_{\mathrm{v}}} - {k_{\mathrm{c}}}{p_{\mathrm{L}}} $$ (4)

    式中,QL为调速阀流量;kq为流量增益;xv为阀芯位移;kc为压力系数;pL为压力腔压力。

    液压缸流量连续性方程[23]

    $$ Q_{\mathrm{L}}=A_{\mathrm{p}}\frac{\mathrm{d}x_{\mathrm{p}}}{\mathrm{d}t}+\frac{V_{\mathrm{t}}}{4\beta_{\mathrm{e}}}\frac{\mathrm{d}p_{\mathrm{L}}}{\mathrm{d}t}+C_{\mathrm{tp}}p_{\mathrm{L}} $$ (5)

    式中,AP为压力腔面积;xp为活塞位移;Vt为压缩容积;βe为弹性系数;Ctp为泄露系数。

    液压缸平衡方程:

    $$ A_{\mathrm{p}}p_{\mathrm{L}}=M_{\mathrm{t}}\frac{\mathrm{d}^2x_{\mathrm{p}}}{\mathrm{d}t^2}+B_{\mathrm{p}}\frac{\mathrm{d}x_{\mathrm{p}}}{\mathrm{d}t}+k_{\mathrm{e}}x_{\mathrm{p}}+F_{\mathrm{L}} $$ (6)

    式中,Mt为总质量;Bp为黏性阻力系数;ke为负载弹性刚度;FL为负载力。

    联立上述式(4)—式(6)3个方程并对其进行拉氏变换得到液压系统的传递函数:

    $$ G\left( s \right) = \dfrac{{{k_{\mathrm{q}}}{x_{\mathrm{v}}} - \dfrac{{{k_{{\mathrm{ce}}}}}}{{{A_{\mathrm{p}}}}}\left( {\dfrac{{{V_t}}}{{4{\beta _{\mathrm{e}}}{k_{{\mathrm{ce}}}}}}s + 1} \right){F_{\mathrm{L}}}}}{{\dfrac{{{V_t}{M_t}}}{{4{\beta _{\mathrm{e}}}{A_{\mathrm{p}}}}}{s^3} + \left( {\dfrac{{{V_t}{B_{\mathrm{p}}}}}{{4{\beta _{\mathrm{e}}}{A_{\mathrm{p}}}}} + \dfrac{{{K_{{\mathrm{ce}}}}{M_t}}}{{{A_{\mathrm{p}}}}}} \right){s^2} + \left( {\dfrac{{{K_{{\mathrm{ce}}}}{B_{\mathrm{p}}}}}{{{A_{\mathrm{p}}}}} + {A_{\mathrm{p}}}} \right)s}} $$ (7)

    式中,kce为总压力系数。

    上式为三阶系统,其结构复杂,系数繁多,运算量大,经典控制理论中的控制器设计难以满足需求。而自抗干扰控制器可将三阶系统中的三阶项和其他扰动转化为一个总扰动,最终将上式复杂的三阶线性系统转化为一阶线性系统[2425]

    将式(7)以动态方程的形式改写为如下表达式:

    $$ {\dot x_{\mathrm{p}}} = - {\eta _1}{\dddot x_{\mathrm{p}}} - {\eta _2}{\ddot x_{\mathrm{p}}} + \omega + bu $$ (8)

    式中:$ \dot{x} $p、$ \ddot{x} $p、${\dddot x_{\mathrm{p}}} $为系统阶数;–η1, –η2为系统参数;ω为扰动;b为控制矩阵;u为控制器输出的控制信号。

    将系统中的高阶项及未建模的扰动部分视为系统的总扰动:

    $$ f\left( \cdot \right) = - {\eta _1}{\dddot x_{\mathrm{p}}} - {\eta _2}{\ddot x_{\mathrm{p}}} + \omega $$ (9)

    则:

    $$ \dot y = f\left( \cdot \right) + bu $$ (10)

    因此简化后的一阶系统控制器的控制结构如图20所示。

    图  20  一阶系统自抗扰控制结构
    Figure  20.  Control structure of ADRC

    非线性扩张状态观测器:

    $$ \left\{ \begin{gathered} {e_1} = {{\textit{z}}_1} - y \\ {{\dot {\textit{z}}}_1} = {{\textit{z}}_2} - {\beta _{01}}{e_1} + {b_0}u \\ {{\dot {\textit{z}}}_2} = - {\beta _{02}}{\left| {{e_1}} \right|^{\tfrac{1}{2}}}{\rm{sign}}({e_1}) \\ \end{gathered} \right. $$ (11)

    式中,z1, z2分别为v1v2的观测信号;z1,z2z1z2的一阶导数;β01β02为待整定系数;b0为控制矩阵;e1为观测误差;y为输出信号。

    设计非线性反馈控制律:

    $$ \left\{\begin{gathered}e_1=v_1-{{\textit{z}}}_1 \\ u_0=\beta_1{\mathrm{f}}\mathrm{al}\left(e,\alpha,\delta\right) \\ u=u_0-\frac{{{\textit{z}}}_2}{b_0} \\ {\mathrm{f}}\mathrm{al}\left(e,a,\delta\right)=\left\{\begin{gathered}\frac{e}{\delta^{1-\alpha}},\left|e\right|\leqslant\delta \\ \left|e\right|^{\alpha}\mathrm{sign}\left(e\right),\left|e\right| > \delta \\ \end{gathered}\right. \\ \end{gathered}\right. $$ (12)

    式中,u0为反馈控制率的控制信号;β1为待整定系数;e为fal的函数变量;δ为影响滤波效果的常数;a为0~1之间的常数。

    图21所示为在Simulink仿真软件中建立的多液压缸的同步控制器模型,其核心是采用了前文所述的均值耦合控制方法和经典控制理论中的自抗扰控制器,并同时建立速度环和位置环进行反馈,以实现液压缸的同步控制。该完整的控制器模型由S-function接口、位置耦合器、自抗扰控制器及同步误差补偿器组成。自抗扰控制器控制速度环,以降低系统受到的干扰,从而实现对液压缸活塞的速度同步控制。PI控制器具有积分环节,可以积累误差并进行补偿,从而提高系统的抗干扰能力,因此位置环采用PI控制器。整定参数见表2

    图  21  Simulink多缸同步控制器模型
    Figure  21.  Simulink multi-cylinder synchronous controller model
    表  2  ADRC控制器参数
    Table  2.  ADRC controller parameters
    待整定参数 b0 β01 β02 β1
    数值 10 5000 −200 10000
    下载: 导出CSV 
    | 显示表格

    图22所示为3种工况下对ADRC控制器和PID控制器的多液压缸同步控制仿真结果。由图22a图22b可以看出ADRC控制器的仿真结果明显优于PID控制器的同步控制效果。如图22c所示,PID控制器受时变负载的影响较大,其最大位移同步误差从最初的0.03 mm扩大至0.05 mm,并且最大位移同步误差仍逐渐增大,而ADRC控制器受时变负载影响较小,最大位移同步误差始终保持在0~0.04 mm范围内波动。通过对3种工况下的仿真结果表明,ADRC控制器比PID控制器具有更好的鲁棒性和抗干扰能力。

    图  22  联合仿真结果
    Figure  22.  Co-simulation results

    上文通过3种不同工况下的仿真,确定了均值耦合同步控制为4种同步控制方法中综合性能最为优异的方法,再对比PID和ADRC算法下均值耦合同步控制的结果,得出基于ADRC算法的均值耦合同步控制策略优于PID算法。为验证所提出自主行走机构和均值耦合同步控制方法结合自抗扰控制器的控制策略的可行性,研制了深部流态化开采装备自主行走机构试验台,并开展了多缸同步控制试验。

    图23a所示为自主行走机构试验台,主要是由采掘舱、转化舱、输出舱3个舱体组成,各舱体间布置4个液压缸,通过液压缸推进舱体前行。采掘舱前布置有负载液压缸,通过调节两个液压缸的压力来模拟3种不同工况。试验台测控软件调试界面如图23b所示。

    图  23  自主行走机构试验台
    Figure  23.  Autonomous walking mechanism test bench

    试验过程中需要对运动元件在线监测,采用图24所示的压力传感器和位移传感器实时采集试验中的压力、位移等信号,并整合到上位机进行显示和存储。传感器参数见表3

    图  24  传感器选型
    Figure  24.  Sensors selection
    表  3  传感器参数
    Table  3.  Sensor parameters
    传感器名称 型号 测量范围 精度 响应时间/
    ms
    工作温度/
    压力传感器 美控
    MIK–P300
    0~25 MPa 0.25%F.S 10 −20~85
    拉杆式位移
    传感器
    Accuiacy
    KTC–500
    0~500 mm 0.05 mm −60~150
    下载: 导出CSV 
    | 显示表格

    由于需进行不同工况负载下的同步试验,因此选取采掘舱与中间舱之间的液压缸进行同步控制试验。

    1)均匀负载下的速度阶跃响应试验。调节负载液压缸压力为零,将试验台空载启动,首先设置推进液压缸以10 mm/s的速度启动,待设备运行稳定后设定液压缸推进速度为20 mm/s。试验结果如图25所示,由图25a可知控制系统在两次不同阶跃信号下响应迅速,速度波动均控制在5 mm/s左右。由图25b可知第1次速度阶跃时最大位移同步误差在2.5 mm以下,第2次速度阶跃时最大位移同步误差在3 mm以下,试验证明均匀负载的速度阶跃情况下控制系统同步性能较好。

    图  25  速度阶跃响应试验结果
    Figure  25.  Speed step response test results

    2)均匀负载下的不同速度同步试验。推进系统在均匀负载时不同速度下最大位移同步误差如图26所示。观察图26曲线发现,当初始速度为5 mm/s时,系统的最大位移同步误差在1.4 mm/s上下波动,当初始速度为10 mm/s时,系统的最大位移同步误差在2.5 mm/s以内波动,当初始速度为15 mm/s时,系统的最大位移同步误差在2.6 mm/s上下波动,当初始速度为20 mm/s时,系统稳定后的最大位移同步误差在3 mm/s上下波动。由此得出推论,多缸推进系统的最大位移同步误差和液压缸活塞的速度有关,给定的液压缸速度越快,最大位移同步误差越大。

    图  26  不同速度下最大位移同步误差曲线
    Figure  26.  Maximum displacement synchronization error curve at different speeds

    1)突变负载同步响应试验。试验时首先调节左右两个负载液压缸压力为1 MPa,然后待液压缸稳定后,在第15秒将左右负载液压缸的压力调至2 MPa,获得试验曲线如图27所示。观察图27a可知,当负载液压缸压力调节至2 MPa时,4个液压缸的速度发生突降,从10 mm/s左右降至8 mm/s左右,7 s后速度恢复正常;观察图27b可知,当负载发生变化时,4个液压缸的最大位移同步误差曲线从2 mm突增至4 mm,随后曲线又恢复正常,保持在2 mm左右波动。

    图  27  突变负载试验结果
    Figure  27.  Sudden load test results

    2)时变负载工况下推进液压缸同步试验。为模拟时变负载工况,在左右两个负载液压缸压力为2 MPa的基础上,引入频率为0.1,振幅为1 MPa的正弦信号,获得如图28所示的试验曲线。由图28a可知,4个推进液压缸速度随正弦信号变化在8 mm至12 mm之间发生波动;由图28b可知多缸最大位移同步误差波动范围在4 mm以下。与突变负载相比,时变负载下的推进系统多液压缸速度和多液的压缸最大位移同步误差较大,然而其最大位移同步误差仍保持在一定范围内,具有较好的鲁棒性。

    图  28  时变载荷下多缸同步试验结果
    Figure  28.  Result of multi-cylinder synchronous test under time-varying load

    通过开展上述试验,分别从均匀负载下速度阶跃变化、推进系统不同速度以及突变负载和时变负载工况等多个方面,检测验证了基于ADRC的均值耦合同步控制策略的优越性,试验结果表明该控制策略下最大同步位移误差均在5 mm以下。

    1)通过开展3种工况下的仿真试验,构建4种不同控制方法的雷达图面积评价模型,结果表明均值耦合控制方法的综合性能最优,但鲁棒性不强。

    2)通过对比仿真分析发现,基于自抗扰控制的均值耦合控制策略比基于PID控制的均值耦合控制策略鲁棒性更高。

    3)通过自主研制自主行走机构试验台,验证了基于自抗扰控制的均值耦合控制策略的控制效果,结果表明所设计液压系统和提出的同步控制策略具有优异的同步性与鲁棒性。

  • 图  1   深地流态化开采系统示意

    Figure  1.   Schematic diagram of deep ground fluidized mining system

    图  2   增阻迈步式自主行走机构结构

    Figure  2.   Structural diagram of increased resistance walking equipment for autonomous walking

    图  3   增阻迈步式自主行走原理

    Figure  3.   Principle of increasing resistance and stepping autonomous walking

    图  4   液压推进系统原理

    1—过滤器;2—恒压变量泵;3—安全阀;4—比例调速阀;5—三位四通换向阀;6—液压缸;7—压力传感器;8—位移传感器

    Figure  4.   Schematic diagram of hydraulic propulsion system

    图  5   推进系统仿真模型

    1—恒压源;2—调速阀;3—推进液压缸;4—速度传感器;5—模拟负载

    Figure  5.   Propulsion system simulation model

    图  6   主从同步控制

    Figure  6.   Master-slave synchronous control

    图  7   相邻交叉耦合同步控制

    Figure  7.   Adjacent cross-coupling synchronous control

    图  8   相邻交叉耦合误差补偿器结构

    Figure  8.   Adjacent cross-coupling error compensator structure

    图  9   偏差耦合同步控制

    Figure  9.   Deviation coupling synchronous control

    图  10   偏差耦合误差补偿器结构

    Figure  10.   Structure of deviation coupling error compensator

    图  11   多缸均值耦合同步控制

    Figure  11.   Multi-cylinder mean coupled synchronous control

    图  12   多缸同步推进控制仿真模型

    1—恒压源;2—比例调速阀;3—推进液压缸;4—速度、位移传感器;5—负载模拟函数;6—位置耦合器;7—同步误差补偿器;8—控制器;9—阈值控制器;10—目标速度信号

    Figure  12.   Simulation model for multi cylinder synchronous propulsion control

    图  13   不同地质工况

    Figure  13.   Different geological conditions

    图  14   不同地质工况加载曲线

    Figure  14.   Loading curves for different geological conditions

    图  15   时变负载模型

    Figure  15.   Time-varying load model

    图  16   均匀负载工况液压缸同步仿真结果

    Figure  16.   Simulation results of hydraulic cylinder synchronization under uniform load conditions

    图  17   突变负载工况液压缸同步仿真结果

    Figure  17.   Simulation results of hydraulic cylinder synchronization under sudden load conditions

    图  18   时变负载液压缸同步仿真结果

    Figure  18.   Synchronous simulation results of time-varying load hydraulic cylinder

    图  19   4种同步控制方法性能雷达图

    Figure  19.   Performance radar chart of four synchronous control methods

    图  20   一阶系统自抗扰控制结构

    Figure  20.   Control structure of ADRC

    图  21   Simulink多缸同步控制器模型

    Figure  21.   Simulink multi-cylinder synchronous controller model

    图  22   联合仿真结果

    Figure  22.   Co-simulation results

    图  23   自主行走机构试验台

    Figure  23.   Autonomous walking mechanism test bench

    图  24   传感器选型

    Figure  24.   Sensors selection

    图  25   速度阶跃响应试验结果

    Figure  25.   Speed step response test results

    图  26   不同速度下最大位移同步误差曲线

    Figure  26.   Maximum displacement synchronization error curve at different speeds

    图  27   突变负载试验结果

    Figure  27.   Sudden load test results

    图  28   时变载荷下多缸同步试验结果

    Figure  28.   Result of multi-cylinder synchronous test under time-varying load

    表  1   控制策略评分标准

    Table  1   Control strategy evaluation criteria

    精度速度鲁棒性成本评分
    3
    2
    1
    下载: 导出CSV

    表  2   ADRC控制器参数

    Table  2   ADRC controller parameters

    待整定参数 b0 β01 β02 β1
    数值 10 5000 −200 10000
    下载: 导出CSV

    表  3   传感器参数

    Table  3   Sensor parameters

    传感器名称 型号 测量范围 精度 响应时间/
    ms
    工作温度/
    压力传感器 美控
    MIK–P300
    0~25 MPa 0.25%F.S 10 −20~85
    拉杆式位移
    传感器
    Accuiacy
    KTC–500
    0~500 mm 0.05 mm −60~150
    下载: 导出CSV
  • [1] 国家统计局. 中华人民共和国2022年国民经济和社会发展统计公报[N]. 人民日报,2023–03–01.
    [2] 秦容军. 我国煤炭开采现状及政策研究[J]. 煤炭经济研究,2019,39(1):57−61.

    QIN Rongjun. Research on current status and policy of coal mining in China[J]. Coal Economic Research,2019,39(1):57−61.

    [3] 何 花,李 鑫,郝成亮. 我国煤炭深部开采战略研究[J]. 煤炭经济研究,2021,41(9):15−51.

    HE Hua,LI Xin,HAO Chengliang. Research on deep coal mining strategy in China[J]. Coal Economic Research,2021,41( 9):15−51.

    [4] 谢和平,鞠 杨,高明忠,等. 煤炭深部原位流态化开采的理论与技术体系[J]. 煤炭学报,2018,43(5):1210−1219.

    XIE Heping,JU Yang,GAO Mingzhong,et al. Theories and technologies for in-situ fluidized mining of deep underground coal resources[J]. Journal of China Coal Society,2018,43( 5):1210−1219.

    [5] 谢和平,高 峰,鞠 杨,等. 深地科学领域的若干颠覆性技术构想和研究方向[J]. 工程科学与技术,2017,49(1):1−8.

    XIE Heping,GAO Feng,JU Yang,et al. Novel idea and disruptive technologies for the exploration and research of deep earth[J]. Advanced Engineering Sciences,2017,49(1):1−8.

    [6] 谢和平,王金华,王国法,等. 煤炭革命新理念与煤炭科技发展构想[J]. 煤炭学报,2017,43(5):1187−1197.

    XIE Heping,WANG Jinhua,WANG Guofa,et al. New ideas of coal revolution and layout of coal science and technology development[J]. Journal of China Coal Society,2017,43(5):1187−1197.

    [7]

    XIE Heping,JU Yang,REN ShiHua,et al. Theoretical and technological exploration of deep in situ fluidized coal mining[J]. Frontiers in Energy,2019,13(4):603−611. doi: 10.1007/s11708-019-0643-x

    [8] 鲍久圣,周恒,王忠宾,等. 一种深部煤炭流态化开采原位爆燃–蒸汽联合发电装置[P]. 中国:CN113250773B,2022–05–24.

    BAO Jiusheng,ZHOU Heng,WANG Zhongbin,et al. A kind of in-situ deflagration-steam combined power generation device for deep coal fluidized mining[P]. China:CN113250773B,2022–05–24.

    [9] 郑瑞良,左 浩,高伟民. 煤矿深部大断面岩巷TBM掘进围岩支护研究[J]. 煤炭科学技术,2020,48(S2):137−142.

    ZUO Zhengliang,ZUO Hao,GAO Weimin. Study on surrounding rock support of TBM driving in deep large section rock roadway of coal mine[J]. Coal Science and Technology,2020,48(S2):137−142.

    [10]

    YANG Huayong,SHI Hu,GONG Guofeng,et al. Electro-hydraulic proportional control of thrust system for shield tunneling machine[J]. Automation in Construction,2009,18(7):950−956. doi: 10.1016/j.autcon.2009.04.005

    [11] 薛 召,贾文昂,李展尚. 基于模糊单神经元PID的比例同步控制研究[J]. 机床与液压,2022,50(11):172−177. doi: 10.3969/j.issn.1001-3881.2022.11.030

    XUE Zhao,JIA Wenang,LI Zhangshang. Research on proportional synchronous control based on fuzzy single neuron PID[J]. Machine Tool & Hydraulics,2022,50(11):172−177. doi: 10.3969/j.issn.1001-3881.2022.11.030

    [12] 徐莉萍,张宇晖,周小磊,等. 隧道掘进机洞内拆装设备多缸同步控制研究[J]. 河南科技大学学报(自然科学版),2023,44(5):31−36,48.

    XU Liping,ZHANG Yuhui,ZHOU Xiaolei,et al. Research on multi cylinder synchronous control of tunnel boring machine tunnel disassembly and assembly equipment[J]. Journal of Henan University of Science and Technology(Natural Science),2023,44(5):31−36,48.

    [13]

    XU Cong,XU Xiao,LIU Zhao,et al. Research on multi-cylinder synchronous control system of multi-directional forging hydraulic press[J]. Journal of Physics:Conference Series,2022,2338(1).

    [14]

    BIJAN Sanati,NAOMI Kato. Trajectory synchronization of multi-cylinder electro-hydraulic lift system with huge load[J]. Journal of Fluids and Structures,2013,14(2):228−234.

    [15] 卢子帅,蔡 伟,赵静一,等. 液压系统群的相邻交叉耦合同步自抗扰控制[J]. 液压与气动,2020(8):29−34. doi: 10.11832/j.issn.1000-4858.2020.08.005

    LU Zishuai,CAI wei,ZHAO Jing yi,et al. Active disturbance rejection control for adjacent cross coupling synchronization of hydraulic system group[J]. Chinese Hydraulics and Pneumatics,2020(8):29−34. doi: 10.11832/j.issn.1000-4858.2020.08.005

    [16] 谢和平,张 茹,邓建辉,等. 基于“深地–地表”联动的深地科学与地灾防控技术体系初探[J]. 工程科学与技术,2021,53(4):1−12.

    XIE Heping,ZHANG Ru,DENG Jianhui,et al. A preliminary study on the technical system of deep earth science and geo disaster prevention-control based on the "deep earth-surface” linkage strategy[J]. Advanced Engineering Sciences,2021,53(4):1−12.

    [17] 鲍久圣,赵少迪,王忠宾,等. 一种深部煤炭流态化开采装备自主行走机构[P]. 中国:202110715377.5,2021–06–27.

    BAO Jiusheng,ZHAO Shaodi,WANG Zhongbin,et al. A autonomous walking mechanism for deep coal fluidized mining equipment[P]. China:202110715377.5,2021–06–27.

    [18] 周 恒. 深部煤炭流态化开采装备自主行走机构设计及控制策略研究[D]. 徐州:中国矿业大学,2022.

    ZHOU Heng. Research on design and control strategy of autonomous walking mechanism for fluidized mining equipment of deep underground coal resources[D]. Xuzhou:China University of Mining and Technology,2022.

    [19] 彭 左. TBM岩土载荷特征识别及姿态控制研究[D]. 杭州:浙江大学,2017.

    PENG Zuo. Research on excavation characteristic identification and attitude control of tunnel boring machine[D]. Hangzhou:Zhejiang University,2017.

    [20] 王利明,杨振兴,周建军,等. 不同围岩地质条件下盾构掘进效能评价方法[J]. 地下空间与工程学报,2019,15(1):240−246.

    WANG Liming,YANG Zhenxing,ZHOU Jianjun,et al. Evaluation method of shield tunneling efficiency in different Surrounding rock[J]. Chinese Journal of Underground Space and Engineering,2019,15(1):240−246.

    [21] 刘 宁,张春生,褚卫江,等. 深埋破碎岩体TBM掘进风险评估与应对[J]. 岩石力学与工程学报,2015,34(S1):3348−3355.

    LIU Ning,ZHANG Chunsheng,CHU Weijiang,et al. Risk assessment and response of TBM excavation in deep buried fractured rock mass[J]. Chinese Journal of Rock Mechanics and Engineering,2015,34(S1):3348−3355.

    [22] 陈 鑫,李 东,曾 冰. 基于改进差分进化算法的液压自抗扰控制器设计[J]. 机床与液压,2021,49(23):92−97. doi: 10.3969/j.issn.1001-3881.2021.23.018

    CHEN Xin,LI Dong,ZENG Bing. Design of hydraulic active disturbance rejection controller based on differential evolution algorithm[J]. Machine Tool & Hydraulics,2021,49(23):92−97. doi: 10.3969/j.issn.1001-3881.2021.23.018

    [23] 赵 爽,史红伟,刘小松. 三阶线性自抗扰控制器的液压伺服流量控制[J]. 液压与气压,2021,45(5):149−156.

    ZHAO Shuang,SHI Hongwei,LIU Xiaosong. Hydraulic servo flow control based on third order LADRC[J]. Chinese Hydraulics & Pneumatics,2021,45(5):149−156.

    [24] 司国雷,神英淇,王嘉磊,等. 电液位置伺服系统自抗扰控制研究[J]. 液压与气动,2020(12):14−21. doi: 10.11832/j.issn.1000-4858.2020.12.003

    SI Guolei,SHEN Yingqi,WANG Jialei,et al. Active disturbance rejection control of electro-hydraulic position servo system[J]. Chinese Hydraulics and Pneumatics,2020(12):14−21. doi: 10.11832/j.issn.1000-4858.2020.12.003

    [25] 郝建伟,鲍久圣,葛世荣,等. 带式输送机永磁驱动系统自抗扰同步控制策略[J]. 电机与控制应用,2021,48(9):27−35.

    HAO Jianwei,BAO Jiusheng,GE Shirong,et al. Active disturbance rejection synchronous control strategy of permanent magnet motor drive system for belt conveyor[J]. Electric Machines and Control Application,2021,48(9):27−35.

  • 期刊类型引用(2)

    1. 陈义,史宏宝,张世锐,刘晓龙. 近距离煤层下分层巷道变形破坏机理及控制技术. 科技和产业. 2024(17): 318-322 . 百度学术
    2. 付文俊,贾振刚,罗应奇. 立井施工智能化技术现状与展望. 煤炭工程. 2024(11): 68-75 . 百度学术

    其他类型引用(2)

图(28)  /  表(3)
计量
  • 文章访问数:  91
  • HTML全文浏览量:  19
  • PDF下载量:  30
  • 被引次数: 4
出版历程
  • 收稿日期:  2023-11-19
  • 网络出版日期:  2024-02-23
  • 刊出日期:  2024-02-22

目录

/

返回文章
返回