CN117094253B - 推拼同步下基于推力均匀性的油缸推力分配方法及系统 - Google Patents
推拼同步下基于推力均匀性的油缸推力分配方法及系统 Download PDFInfo
- Publication number
- CN117094253B CN117094253B CN202311360805.2A CN202311360805A CN117094253B CN 117094253 B CN117094253 B CN 117094253B CN 202311360805 A CN202311360805 A CN 202311360805A CN 117094253 B CN117094253 B CN 117094253B
- Authority
- CN
- China
- Prior art keywords
- thrust
- oil cylinder
- partition
- oil
- cylinder
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 title claims description 29
- 238000005192 partition Methods 0.000 claims abstract description 94
- 238000005457 optimization Methods 0.000 claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 238000010276 construction Methods 0.000 claims abstract description 9
- 239000003921 oil Substances 0.000 claims description 252
- 230000005641 tunneling Effects 0.000 claims description 19
- 230000000903 blocking effect Effects 0.000 claims description 15
- 239000002689 soil Substances 0.000 claims description 9
- 239000010721 machine oil Substances 0.000 claims description 3
- 230000035515 penetration Effects 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 230000001360 synchronised effect Effects 0.000 abstract description 15
- 238000005516 engineering process Methods 0.000 abstract description 7
- 230000008602 contraction Effects 0.000 abstract description 3
- 239000013598 vector Substances 0.000 description 7
- 238000010835 comparative analysis Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000010720 hydraulic oil Substances 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21D—SHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
- E21D9/00—Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
- E21D9/06—Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining
- E21D9/0607—Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining the shield being provided with devices for lining the tunnel, e.g. shuttering
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21D—SHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
- E21D11/00—Lining tunnels, galleries or other underground cavities, e.g. large underground chambers; Linings therefor; Making such linings in situ, e.g. by assembling
- E21D11/40—Devices or apparatus specially adapted for handling or placing units of linings or supporting units for tunnels or galleries
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Evolutionary Computation (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mathematical Optimization (AREA)
- Geochemistry & Mineralogy (AREA)
- Mathematical Analysis (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Architecture (AREA)
- Algebra (AREA)
- Structural Engineering (AREA)
- Civil Engineering (AREA)
- Environmental & Geological Engineering (AREA)
- Computational Mathematics (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
本发明公开了推拼同步下基于推力均匀性的油缸推力分配算法及系统,属于盾构隧道施工技术和装备领域。基于管片分块的油缸分区,每个分区内至少设置一组油缸,每组油缸组内均设置至少一根油缸。构建油缸分区等效模型,油缸分区等效模型为合力和等效半径的显式函数;建立油缸组推力的约束优化模型,约束优化模型的目标函数为未回缩油缸的方差最小值,对约束优化模型采用优化算法得到合力的最优解,基于合力的最优解得到所需的单根油缸的推力。本发明将作用在每块管片上的油缸视为一个油缸分区,每个分区的油缸推力大小相同和伸缩同步。计算模型符合油缸分区独立控制的理念,可减少未知量的个数,提高了计算效率的同时不影响计算精度。
Description
技术领域
本发明属于盾构隧道施工技术和装置的技术领域,特别是涉及推拼同步下基于推力均匀性的油缸推力分配算法及系统。
背景技术
常规盾构施工采用推拼交替模式,即先掘进一环管片的宽度后停机拼装管片,待整环管片拼装完毕后,再进行下一循环。常规推拼模式的推拼时间可划分为掘进时间和拼装时间,两者时间大体相当,从而导致拼装效率较低。
为了克服常规推拼模式掘进效率较低的问题,日本首先提出了“管片同步拼装施工”的理念,即掘进的同时进行管片的拼装,相较于常规拼装模式,推拼同步模式可以节省近一半的时间。近年来,随着一大批长大隧道的规划与修建,推拼同步的技术优势充分体现。推拼同步技术经过不断地改进发展,油缸推力矢量控制技术成为一种可实施的技术方案。
在推拼同步过程中,边掘进边拼装管片,部分油缸需要缩回,因此盾构推力的矢量控制技术是推拼同步系统的关键,即通过算法控制油缸收缩前后总推力的大小和方向不变。现有的推拼同步推力矢量控制技术主要有三种:一种方式是缩回一组油缸用于管片拼装,再缩回一组与之对称的油缸(也称为同调)以克服偏载,其余油缸增大推力以实现合力点和总推力不变;另一种方式是使得未缩回的油缸形成梯度化的推进力以补偿缩回油缸的推进力;第三种方式是使得各油缸推力的变化值最小。上述三种方法均会导致油缸的性能无法充分发挥。
发明内容
本发明为解决上述背景技术中的问题,提供了推拼同步下基于推力均匀性的油缸推力分配算法及系统。
本发明采用以下技术方案:推拼同步下基于推力均匀性的油缸推力分配算法,至少包括以下步骤:每环管片由若干管片分块组成,于盾体内配置多组推进油缸组;推进时,每个管片分块上至少支撑有一组油缸;构建基于管片分块的油缸分区等效模型,所述管片分块油缸分区等效模型为将支撑于同一管片分块上的所有油缸视为一个油缸分区,分区内油缸推力的合力大小为Fj,等效半径为dj;建立油缸组推力的约束优化模型,所述约束优化模型的目标函数为使未缩回油缸推力的方差达到最小值,所述约束优化模型的约束条件包括:目标总推力大小约束函数、推力偏心量约束函数和油缸组推力范围约束函数;对所述约束优化模型采用优化算法得到各油缸分区的合力Fj的最优解,基于合力Fj的最优解得到所需的单根油缸推力fj。
在进一步的实施例中,根据需求设计和施工要求,每环管片由n块管片分块组成,每块管片分块j均配置有kj组油缸组,每组油缸组中含有根油缸;其中,k≥1,/>1≤j≤n;/>为管片分块j内油缸组kj中的油缸编号。
在进一步的实施例中,所述管片分块的油缸分区等效模型的构建流程如下:
取n块管片分块中的任意一块,以管片的圆心为坐标系原点(x0,y0),以选定的管片分块的对称轴为y’轴,建立Ox’y’局部坐标系:通过等效公式(1)将管片分块j上的所有油缸组划分为同一分区,将同一分区内的油缸组进行等效为一根油缸;
式中,为油缸组k中的油缸编号,/>表示管片分块j中第k组油缸中编号为/>的油缸推力;x’j为合力Fj在管片分块j的等效作用点的x轴坐标值,y’j为合力Fj在管片分块j的等效作用点的y’轴坐标值,θj为合力Fj的等效作用点与坐标原点的连线与Ox’轴的夹角;
其中,所述x’轴坐标值x’j的计算公式如下:
y’轴坐标值y’j的计算公式如下:式中,r0表示油缸推力点到坐标系原点(x0,y0)的距离,/>表示管片分块j中第k组油缸中编号为/>的油缸推力点与坐标系原点(x0,y0)的连线与坐标横轴的夹角。
在进一步的实施例中,还包括以下步骤:定义同一分区内的油缸推力相同,油缸在管片上沿对称轴对称分布,则所述等效公式(1)简化为:
其中,fj为管片分块j中任意一个油缸的推力;
对应的,x’j的取值简化为0,y’j的取值通过以下公式计算得到:
在进一步的实施例中,所述约束优化模型的目标函数具体如下:
式中,g为优化目标函数,即拼装时剩余n-1个管片分块对应等效油缸推力的方差;/>为剩余n-1个等效油缸的推力平均值,/>
在进一步的实施例中,所述总推力大小约束函数如下:
其中,F0为盾构掘进所需的目标总推力值;所述目标总推力值F0的计算公式如下:
F0=F阻+F侧压+F土压+F摩+F牵引,其中,F阻为刀盘贯入阻力,F侧压为作用在刀盘面板上的地层侧向压力,F土压为刀盘开口处的土仓压力,F摩为盾构推进过程中盾壳与地层间的摩擦力,F牵引为后配套台车牵引力。
在进一步的实施例中,建立以管片的圆心为坐标系原点水平方向为x轴,竖直方向为y轴的整体坐标系Oxy。所述推力偏心量约束函数包括:x轴方向偏心约束函数和y轴方向偏心约束函数。
所述x轴方向偏心约束函数表示为:
其中,αj表示管片分块j的油缸组的合力Fj作用点与x轴坐标之间的夹角,ex0表示对x轴的目标偏心距。
所述y轴方向偏心约束函数表示为:
其中,ey0表示对y轴的目标偏心距。
在进一步的实施例中,油缸组推力范围约束函数具体如下:
Fj,min≤Fj≤Fj,max;其中,Fj,min为保证安全拼装,第j根等效油缸所需最小推力,采用以下公式计算得到:
其中Gj表示第j块管片的自重,μ1为油缸撑靴与管片之间的摩擦系数;μ2为相邻两环管片之间的摩擦系数;Fj,max为第j根等效油缸最大推力,采用以下公式计算得到:/>其中,fmax表示单根油缸的最大推力。
在进一步的实施例中,所述对x轴的目标偏心距ex0的计算公式如下:
其中MRx为绕x轴的纠偏力矩,其大小与围岩情况、盾构姿态有关。
所述对y轴的目标偏心距ey0的计算公式如下:
其中MRy为绕y轴的纠偏力矩,其大小与围岩情况、盾构姿态有关。
推拼同步下基于推力均匀性的油缸推力分配系统,用于实现如上所述的分配算法;所述系统包括:设置在盾体内的硬件单元和连接于所述硬件单元的软件单元。
其中,所述硬件单元包括:若干油缸组,沿环向等间距分布在所述盾体内的油缸支撑架上;每组油缸由若干根油缸通过油缸撑靴联系在一起,所述油缸撑靴用于掘进时支撑在管片的横截面上。
在进一步的实施例中,所述油缸与油缸撑靴之间通过球铰铰接。
在进一步的实施例中,所述软件单元用于控制所述油缸组内单个油缸的推力大小,包括:
第一模块,根据管片分块,对盾构机油缸进行编号和分区,构建基于管片分块的油缸分区等效模型,所述管片分块油缸分区等效模型为将支撑于同一管片分块上的所有油缸视为一个油缸分区,分区内油缸推力的合力大小为Fj,作用等效半径为dj。
第二模块,建立油缸组推力的约束优化模型,所述约束优化模型的目标函数为使未回缩油缸推力的方差达到最小值,所述约束优化模型的约束条件包括:目标总推力大小约束函数、推力偏心量约束函数和油缸组推力范围约束函数。
第三模块,根据设定的目标总推力值和目标偏心距,对所述约束优化模型采用优化算法得到合力Fj的最优解,基于合力Fj的最优解得到所需的单根油缸推力fj,并将推力分配给相应的油缸。
本发明的有益效果:本发明提出了一种基于管片分块的油缸分区的计算模型,即将作用在每块管片上的油缸视为一个油缸分区,每个分区的油缸推力大小相同和伸缩同步。本计算模型首先符合油缸分区独立控制的理念,同时可减少未知量的个数,提高了计算效率的同时不影响计算精度。
与此同时,在保证推力矢量相同的条件下,尽可能地使各油缸的推力均匀,增加了油缸的寿命。
附图说明
图1为推拼同步模式下油缸压力均匀分配方法的流程图。
图2为管片分块及油缸分区示意图。
图3为单块管片的油缸组等效简化模型图。
图4为基于管片分块的油缸分区推力计算原理图。
图5为油缸推力合力示意图。
图6(a)为实施例1中目标总推力为100MN的同步推拼模式和非同步推拼模式下油缸推力对比分析图。
图6(b)为实施例1中目标总推力为150MN的同步推拼模式和非同步推拼模式下油缸推力对比分析图。
图6(c)为实施例1中目标总推力为211.1MN的同步推拼模式和非同步推拼模式下油缸推力对比分析图。
图6(d)为实施例1中目标总推力为232MN的同步推拼模式和非同步推拼模式下油缸推力对比分析图。
图7为实施例2中盾构机推进系统的三维结构示意图。
图8为某项目中盾构油缸推力作用点位置详图。
图1至图8中的各标注为:管片分块1、油缸撑靴2、油缸3、球铰4、盾体5、第一油缸分区1-1、第二油缸分区1-2、第三油缸分区1-3、第四油缸分区1-4、第五油缸分区1-5、第六油缸分区1-6、第七油缸分区1-7、第八油缸分区1-8、第九油缸分区1-9、第十油缸分区1-10。
具体实施方式
下面结合说明书附图和实施例对本发明做进一步的描述。
实施例1
如图1所示,推拼同步下基于推力均匀性的油缸推力分配算法,包括以下步骤:
步骤一、每环管片由若干个管片分块1组成,于盾体内配置多组推进油缸组;推进时,每个管片分块1上均分配有至少一组油缸组;
步骤二、构建基于管片分块的油缸分区等效模型,所述管片分块油缸分区等效模型为将支撑于同一管片分块上的所有油缸视为一个油缸分区,分区内油缸推力的合力大小为Fj,作用等效半径为dj;
步骤三、建立油缸组推力的约束优化模型,所述约束优化模型的目标函数为未回缩油缸推力的方差最小值,所述约束优化模型的约束条件包括:目标总推力大小约束函数、推力偏心量约束函数和油缸组推力范围约束函数;
步骤四、对所述约束优化模型采用优化算法得到合力Fj的最优解,基于合力Fj的最优解得到所需的单根油缸推力fj。
在进一步的实施例中,根据需求采用等分模式将管片均分成n块管片分块,每块管片分块j均配置有kj组油缸组,每组油缸组中含有根油缸;其中,k≥1,/>1≤j≤n;/>为管片分块j内油缸组kj中的油缸编号。
结合图2,每环管片由n块管片分块组成,即n=10。每个管片分块上均支撑有两组油缸组,即kj=k1=k2…=k10=2。根据需求,每组油缸组中的油缸数量可为相同或者不同:对于一般情况,每组油缸组中的油缸数量应相等,如图2所示,换言之,每组油缸组中均设置有3根油缸。对于管片非均匀分块的情况,若采用上述油缸均匀分配方案,将导致油缸撑靴骑压管片纵缝时,每块管片上支撑的油缸组以及油缸组内油缸数量可不相同。举例,10块管片分块中的第一块管片分块上的两组油缸组分别为油缸组11和油缸组21,其中,油缸组11中的油缸数量为3,油缸组21中的油缸数量为2,则表示为:/>和/>同一管片分块上的不同油缸组内的油缸数量不同且分别为三个和两个。
以每组油缸组中的油缸数量相同为例,推进时支撑在同一块管片分块上的油缸推力相等,假设不考虑拼装误差以及管片位移引起的推力偏转则油缸推力方向垂直于管片横截面。
结合图3和图4,所述管片分块油缸分区等效模型的构建流程如下:取n块管片分块中任意一块,以管片的圆心为坐标系原点(x0,y0),以选定的管片分块的对称轴为y’轴,建立Ox’y’局部坐标系:通过等效公式(1)将管片分块j上的所有油缸组划分为同一分区,将同一分区内的油缸组进行等效;
式中,为油缸组k中的油缸编号,/>表示管片分块j中第k组油缸中编号为/>的油缸推力;x’j为合力Fj在管片分块j的等效作用点的x’轴坐标值,y’j为合力Fj在管片分块j的等效作用点的y’轴坐标值;θj为合力Fj的等效作用点与坐标原点的连线与Ox’轴的夹角。
其中,所述x’轴坐标值x’j的计算公式如下:
y’轴坐标值y’j的计算公式如下:式中,r0表示油缸推力点到坐标系原点(x0,y0)的距离,/>表示管片分块j中第k组油缸中编号为/>的油缸推力点与坐标系原点(x0,y0)的连线与坐标横轴的夹角。
在执行步骤二时,将作用在每块管片分块上的油缸组视为一个油缸分区,每个分区内的油缸推力大小相同和伸缩同步。因此,本发明中提出的基于管片分块的油缸分区等效模型符合油缸分区独立控制的理念,同时可减少未知量的个数。本实施例中,n=10,即将管片等分成10块管片分块,共有60根油缸,每个油缸组设置有3根油缸,每块管片上支撑着2组油缸,可将支撑在同一块管片上的两组油缸划分同一油缸分区,得到以下油缸分区:第一油缸分区1-1、第二油缸分区1-2、第三油缸分区1-3、第四油缸分区1-4、第五油缸分区1-5、第六油缸分区1-6、第七油缸分区1-7、第八油缸分区1-8、第九油缸分区1-9、第十油缸分区1-10。
基于管片分块油缸分区等效模型,同一管片分块上的油缸推力的大小和方向均相同,将作用在每块管片上的油缸等效成一跟大油缸,等效原则为简化前后力的作用点和大小相同。图4中的F1、F2、…、F10则分别为第一油缸分区、第二油缸分区、…、第十管片分块上等效后的油缸组合力。对应的,d1、d2、…、d10则分别为第一管片分块、第二管片分块、…、第十管片分块上作用点的等效半径。图中,α3表示第三管片分块的油缸组的合力F3作用点与x轴坐标之间的夹角,α5表示第五油缸的油缸组的合力F5作用点与x轴坐标之间的夹角。
需要说明的是,专利CN114135300A曾提出各油缸推力变化最小的思想,但是基于全局考虑上述分配方法会加剧不同油缸的推力不均匀性,同时其计算模型将每根油缸推力视为优化变量,将大大增加计算工作量,降低计算效率。专利CN111810174A曾提出了存在四分区油缸推力计算模型,该模型将所有油缸等效为四组油缸,对于推拼同步模式不利于油缸推力的精确分配与盾构姿态的精确控制。
相比之下,本实施例提出的管片分块油缸分区等效模型中仅存在十个未知量,在提高了计算效率的同时不影响计算精度。
为了便于盾构姿态的调整以及防止管片在推拼过程中出现偏压,油缸一般对称支撑于管片上,由于每块管片上的油缸为同一分区,同一分区的油缸推力相同,所述等效公式(1)简化为:
其中,fj为管片分块j中任意一个油缸的推力。
对应的,x’j的取值简化为0,y’j的取值通过以下公式计算得到:
拼装某块管片分块时,支撑在对应管片分块上的kj组油缸组缩回,则对应的kj组油缸组缩回的推力为零。举例,若第二块管片分块(j=2)为待拼装管片,则第二块管片分块(j=2)上的k2组油缸组缩回后油缸推力为零。为了保证油缸收缩前后推力的矢量保持不变,应将缺失的相应的油缸顶推力分配给其他(n-1)kj组油缸,分配的原则为保证油缸缩回前后,推力矢量保持不变,即保持推力的大小和偏心量不变。
以图1为例,在拼装最顶部管片时,需将支撑在第二块管片分块(j=2)的两组油缸(6根)缩回,此时两组油缸的推力为零。为控制盾构掘进姿态,要保证油缸缩回前后,推力合力的大小和方向保持不变,因此,油缸收缩后,油缸推力的重分配是推拼同步系统的关键。
本实施例在保证推力矢量相同的条件下,尽可能地使各油缸的推力均匀,故以各油缸推力的方差最小为优化目标,即所述约束优化模型的目标函数具体如下:式中,g为优化目标函数,即推进时剩余n-1个管片分块对应油缸组的方差;/>为剩余n-1个k组油缸组的推力平均值,/>F0为盾构掘进所需的目标总推力值。
步骤三中的总推力大小约束函数如下:
其中,所述目标总推力值F0的计算公式如下:
F0=F阻+F侧压+F土压+F摩+F牵引,其中,F阻为刀盘贯入阻力,F侧压为作用在刀盘面板上的地层侧向压力,F土压为刀盘开口处的土仓压力,F摩为盾构推进过程中盾壳与地层间的摩擦力,F牵引为后配套台车牵引力。
推力偏心量约束函数包括:x轴方向偏心约束函数和y轴方向偏心约束函数;所述x轴方向偏心约束函数表示为:
其中,αj表示管片分块j的油缸组的合力Fj作用点与x轴坐标之间的夹角,ex0表示x轴方向目标偏心距。
所述y轴方向偏心约束函数表示为:
其中,ey0表示y轴方向目标偏心距。
在本实施例中,所述对x轴的目标偏心距ex0的计算公式如下:其中MRx为绕x轴的纠偏力矩,其大小与围岩情况、盾构姿态有关。
所述对y轴的目标偏心距ey0的计算公式如下:其中MRy为绕y轴的纠偏力矩,其大小与围岩情况、盾构姿态有关。
在进一步的实施例中,油缸组推力范围约束函数具体如下:
Fj,min≤Fj≤Fj,max;其中,Fj,min为管片分块内油缸最小推力,采用以下公式计算得到:
其中Gj表示第j块管片的自重,μ1为油缸撑靴与管片之间的摩擦系数;μ2为管片与管片之间的摩擦系数;Fj,max为管片分块内油缸最大推力,采用以下公式计算得到:/>其中,fmax表示单根油缸的最大推力。举例说明:某项目中,隧道长13.1km,盾构长度11.8km,拟采用同步推拼技术进行掘进。隧道管片外径14.3m,环宽2m,厚度600mm,采用10等分分块形式,每块管片的角度36°。推进系统包含沿管片均匀分布的20组三联油缸,每块管片上支撑2组油缸,每组油缸由三根组成,共60根。油缸参数:直径420/280mm,油缸行程4000mm,最大油压35MPa,故单根油缸最大推力fmax=4849kN。油缸支撑点位置如图8所示,r0=14m,每块管片支撑有2组油缸,两组油缸之间的夹角为9°,每组油缸内有3根油缸,油缸推力作用点的夹角为4.5°,故图3中的θj 1,1=76.5°,同理可得到其他油缸推力作用点的夹角。
代入到公式(1)中,可得到:
基于管片分块内油缸最大推力的计算公式,得到
Fj,max=6×4849=29094kN。
建立10等分油缸推力计算模型,分别对油缸分区进行编号。根据2.2与3.1.1可知,d1=…=dn=…=d10=13.80m,α1=18°,α2=54°,同理可得α3~α10。
如图6(a)~图6(d)给出了偏心量为0,总推力为100MN、150MN、211.1MN和232MN时等效油缸的推力。1-10油缸组收回后,为了保证偏心距为0,缺失的推力由同侧油缸分担,因此同侧油缸推力的推力要大于对侧油缸的推力。为保证推拼过程总推力大小不变,各个油缸均有不同比例的增大。在油缸分区推力达到极限推力之前,各个油缸的推力均成比例增加,其中第一油缸分区1-1和第九油缸分区1-9增加最大,增加22%。第五油缸分区1-5油缸推力增加最小,分别为2.2%。当第一油缸分区1-1和第九油缸分区1-9推力达到极限推力后,两侧油缸推力依次增加,对侧分区油缸推力(第五油缸分区1-5)减小,当两侧分区油缸推力达到最大,对侧油缸为0时,总推力达到最大。
实施例2
本实施例公开了推拼同步模式下油缸3的推力均匀分配系统,用于实现实施例1所述的推拼同步模式下油缸3的推力均匀分配方法。如图7所示,包括:设置在盾体5内的硬件单元和连接于所述硬件单元的软件单元。
其中,硬件单元包括:若干组推进油缸组,沿环向等间距分布在所述盾体5内的油缸支撑架上;每3根油缸通过油缸撑靴2联系在一起组成一个油缸组,所述油缸撑靴2用于掘进时支撑在管片的横截面上。在本实施例中,油缸组实际为20组,每组油缸组中均设置有三根油缸3,位于同一油缸组中的油缸3通过球铰4连接于同一个油缸撑靴2。为保证油缸3的稳定性和管片横截面受压均匀,一般将2~3根油缸3划分为一组,每组油缸3两端分别设置一个油缸撑靴2,油缸撑靴2在掘进时支撑在管片的横截面上。
为了更好的控制每块管片上的油缸组,以及每组油缸组内的油缸推力,本实施例中的软件单元被设置为控制所述油缸组内单个油缸的推力大小,包括:
第一模块,根据管片分块,对盾构机油缸进行编号和分区,构建基于管片分块的油缸分区等效模型,所述基于管片分块的油缸分区等效模型为将支撑于同一管片分块上的所有油缸视为一个油缸分区,分区内油缸推力的合力大小为Fj,作用等效半径为dj;
第二模块,建立油缸组推力的约束优化模型,所述约束优化模型的目标函数为使未回缩油缸推力的方差达到最小值,所述约束优化模型的约束条件包括:总推力大小约束函数、推力偏心量约束函数和油缸组推力范围约束函数;
第三模块,根据设定的目标总推力值和目标偏心距,对所述约束优化模型采用优化算法得到合力Fj的最优解,基于合力Fj的最优解得到所需的单根油缸推力fj,并将推力分配给相应的油缸。
Claims (8)
1.推拼同步下基于推力均匀性的油缸推力分配方法,其特征在于,至少包括以下步骤:
每环管片由若干管片分块组成,于盾体内配置多组推进油缸组;推进时,每个管片分块上至少支撑有一组油缸;构建基于管片分块的油缸分区等效模型,所述管片分块的油缸分区等效模型为将支撑于同一管片分块上的所有油缸视为一个油缸分区,分区内油缸推力的合力大小为Fj,等效半径为dj;
建立油缸组推力的约束优化模型,所述约束优化模型的目标函数为使未缩回油缸推力的方差达到最小值,所述约束优化模型的约束条件包括:总推力大小约束函数、推力偏心量约束函数和油缸组推力范围约束函数;
对所述约束优化模型采用优化算法得到合力Fj的最优解,基于合力Fj的最优解得到所需的单根油缸推力fj;
所述管片分块的油缸分区等效模型的构建流程如下:
取n块管片分块中的任意一块,以管片的圆心为坐标系原点(x0,y0),以选定的管片分块的对称轴为y’轴,建立Ox’y’局部坐标系:通过等效公式(1)将管片分块j上的所有油缸组划分为同一分区,将同一分区内的所有油缸进行等效为一根油缸;
式中,为油缸组kj中的油缸编号,/>表示管片分块j中第k组油缸中编号为/>的油缸推力;x’j为合力Fj在管片分块j的等效作用点的x’轴坐标值,y’j为合力Fj在管片分块j的等效作用点的y’轴坐标值,θj为合力Fj的等效作用点与坐标原点的连线与Ox’轴的夹角;
其中,所述x’轴坐标值x’j的计算公式如下:
y’轴坐标值y’j的计算公式如下:式中,r0表示油缸推力点到坐标系原点(x0,y0)的距离,/>表示管片分块j中第k组油缸中编号为/>的油缸推力点与坐标系原点(x0,y0)的连线与坐标横轴的夹角;
所述总推力大小约束函数如下:
其中,F0为盾构掘进所需的目标总推力值;所述目标总推力值F0的计算公式如下:
F0=F阻+F侧压+F土压+F摩+F牵引,其中,F阻为刀盘贯入阻力,F侧压为作用在刀盘面板上的地层侧向压力,F土压为刀盘开口处的土仓压力,F摩为盾构推进过程中盾壳与地层间的摩擦力,F牵引为后配套台车牵引力;
建立以管片的圆心为坐标系原点O(x0,y0),水平方向为x轴,竖直方向为y轴的整体坐标系Oxy;所述推力偏心量约束函数包括:x轴方向偏心约束函数和y轴方向偏心约束函数;
所述x轴方向偏心约束函数表示为:
其中,αj表示管片分块j的油缸组的合力Fj作用点与x轴坐标之间的夹角,ex0表示对x轴的目标偏心距;
所述y轴方向偏心约束函数表示为:
其中,ey0表示对y轴的目标偏心距;
油缸组推力范围约束函数具体如下:
Fj,min≤Fj≤Fj,max;其中,Fj,min为保证安全拼装,第j根等效油缸所需最小推力,采用以下公式计算得到:
其中Gj表示第j块管片的自重,μ1为油缸撑靴与管片之间的摩擦系数;μ2为相邻两环管片之间的摩擦系数;
Fj,max为第j根等效油缸的最大推力,采用以下公式计算得到:其中,fmax表示单根油缸的最大推力。
2.根据权利要求1所述的推拼同步下基于推力均匀性的油缸推力分配方法,其特征在于,根据需求,每环管片由n块管片分块组成,每块管片分块j均配置有kj组油缸组,每组油缸组中含有根油缸;其中,k≥1,/>1≤j≤n;/>为管片分块j内油缸组kj中的油缸编号。
3.根据权利要求1所述的推拼同步下基于推力均匀性的油缸推力分配方法,其特征在于,还包括以下步骤:定义同一分区内的油缸推力相同,油缸在管片上沿对称轴对称分布,则所述等效公式(1)简化为:
其中,fj为管片分块j中任意一个油缸的推力;
对应的,x’j的取值简化为0,y’j的取值通过以下公式计算得到:
4.根据权利要求1所述的推拼同步下基于推力均匀性的油缸推力分配方法,其特征在于,所述约束优化模型的目标函数具体如下:
式中,g为优化目标函数,即拼装时剩余n-1个管片分块对应等效油缸推力的方差;/>为剩余n-1个等效油缸的推力平均值,/>F0为盾构掘进所需的目标总推力值。
5.根据权利要求1所述的推拼同步下基于推力均匀性的油缸推力分配方法,其特征在于,
所述对x轴的目标偏心距ex0的计算公式如下:其中MRx为绕x轴的纠偏力矩;
所述对y轴的目标偏心距ey0的计算公式如下:其中MRy为绕y轴的纠偏力矩。
6.推拼同步下基于推力均匀性的油缸推力分配系统,用于实现如权利要求1至5中任意一项所述的分配方法;其特征在于,所述系统包括:设置在盾体内的硬件单元和连接于所述硬件单元的软件单元;
其中,所述硬件单元包括:若干油缸组,沿环向等间距分布在所述盾体内的油缸支撑架上;每组油缸由若干根油缸通过油缸撑靴联系在一起,所述油缸撑靴用于掘进时支撑在管片的横截面上。
7.根据权利要求6所述的推拼同步下基于推力均匀性的油缸推力分配系统,其特征在于,所述油缸与油缸撑靴之间通过球铰铰接。
8.根据权利要求6所述的推拼同步下基于推力均匀性的油缸推力分配系统,其特征在于,
所述软件单元用于控制所述油缸组内单个油缸的推力大小,包括:
第一模块,根据管片分块,对盾构机油缸进行编号和分区,构建基于管片分块的油缸分区等效模型,所述管片分块油缸分区等效模型为将支撑于同一管片分块上的所有油缸视为一个油缸分区,分区内油缸推力的合力大小为Fj,作用等效半径为dj;
第二模块,建立油缸组推力的约束优化模型,所述约束优化模型的目标函数为使未回缩油缸推力的方差达到最小值,所述约束优化模型的约束条件包括:总推力大小约束函数、推力偏心量约束函数和油缸组推力范围约束函数;
第三模块,根据设定的目标总推力值和目标偏心距,对所述约束优化模型采用优化算法得到合力Fj的最优解,基于合力Fj的最优解得到所需的单根油缸推力fj,并将推力分配给相应的油缸。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311360805.2A CN117094253B (zh) | 2023-10-20 | 2023-10-20 | 推拼同步下基于推力均匀性的油缸推力分配方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311360805.2A CN117094253B (zh) | 2023-10-20 | 2023-10-20 | 推拼同步下基于推力均匀性的油缸推力分配方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117094253A CN117094253A (zh) | 2023-11-21 |
CN117094253B true CN117094253B (zh) | 2024-01-23 |
Family
ID=88780270
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311360805.2A Active CN117094253B (zh) | 2023-10-20 | 2023-10-20 | 推拼同步下基于推力均匀性的油缸推力分配方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117094253B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117892044B (zh) * | 2024-03-18 | 2024-05-24 | 北京城建设计发展集团股份有限公司 | 盾构推进系统在不同位姿状态下产生三向末端力能力的评估方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201659082U (zh) * | 2010-03-04 | 2010-12-01 | 西安捷盛电子技术有限责任公司 | 机械式振动台四轴惯性式激振器偏心距同步驱动系统 |
WO2012080241A1 (en) * | 2010-12-14 | 2012-06-21 | Bluewater Energy Services B.V. | Thrust allocation in dynamic positioning systems |
WO2015180196A1 (zh) * | 2014-05-27 | 2015-12-03 | 中国空间技术研究院 | 一种基于多学科优化的卫星总体方案确定系统及实现方法 |
CN111810174A (zh) * | 2020-07-24 | 2020-10-23 | 上海隧道工程有限公司 | 推拼同步模式下盾构推进系统顶力分配计算方法 |
CN112832789A (zh) * | 2021-03-09 | 2021-05-25 | 上海隧道工程有限公司 | 盾构机推拼同步推进系统的阀控方法及系统 |
CN114135300A (zh) * | 2021-11-22 | 2022-03-04 | 中铁工程装备集团有限公司 | 盾构机同步推拼系统推进力分配方法、装置、设备和介质 |
CN114635698A (zh) * | 2022-03-22 | 2022-06-17 | 上海隧道工程有限公司 | 盾构推拼同步中基于同块同压的压力控制方法及系统 |
CN115478864A (zh) * | 2022-09-05 | 2022-12-16 | 北京交通大学 | 一种推拼同步模式下盾构推进系统推进力重分配计算方法 |
CN116856943A (zh) * | 2023-07-12 | 2023-10-10 | 中交疏浚技术装备国家工程研究中心有限公司 | 基于共面约束的四分区盾构推力及油压智能均衡分配方法 |
CN116856944A (zh) * | 2023-07-12 | 2023-10-10 | 中交疏浚技术装备国家工程研究中心有限公司 | 基于力矩补偿的盾构机智能掘进油压均衡分配方法 |
CN116882129A (zh) * | 2023-06-05 | 2023-10-13 | 华中科技大学 | 恒定/非恒定力学约束下盾构推进系统推力分配计算方法 |
-
2023
- 2023-10-20 CN CN202311360805.2A patent/CN117094253B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201659082U (zh) * | 2010-03-04 | 2010-12-01 | 西安捷盛电子技术有限责任公司 | 机械式振动台四轴惯性式激振器偏心距同步驱动系统 |
WO2012080241A1 (en) * | 2010-12-14 | 2012-06-21 | Bluewater Energy Services B.V. | Thrust allocation in dynamic positioning systems |
WO2015180196A1 (zh) * | 2014-05-27 | 2015-12-03 | 中国空间技术研究院 | 一种基于多学科优化的卫星总体方案确定系统及实现方法 |
CN111810174A (zh) * | 2020-07-24 | 2020-10-23 | 上海隧道工程有限公司 | 推拼同步模式下盾构推进系统顶力分配计算方法 |
CN112832789A (zh) * | 2021-03-09 | 2021-05-25 | 上海隧道工程有限公司 | 盾构机推拼同步推进系统的阀控方法及系统 |
CN114135300A (zh) * | 2021-11-22 | 2022-03-04 | 中铁工程装备集团有限公司 | 盾构机同步推拼系统推进力分配方法、装置、设备和介质 |
CN114635698A (zh) * | 2022-03-22 | 2022-06-17 | 上海隧道工程有限公司 | 盾构推拼同步中基于同块同压的压力控制方法及系统 |
CN115478864A (zh) * | 2022-09-05 | 2022-12-16 | 北京交通大学 | 一种推拼同步模式下盾构推进系统推进力重分配计算方法 |
CN116882129A (zh) * | 2023-06-05 | 2023-10-13 | 华中科技大学 | 恒定/非恒定力学约束下盾构推进系统推力分配计算方法 |
CN116856943A (zh) * | 2023-07-12 | 2023-10-10 | 中交疏浚技术装备国家工程研究中心有限公司 | 基于共面约束的四分区盾构推力及油压智能均衡分配方法 |
CN116856944A (zh) * | 2023-07-12 | 2023-10-10 | 中交疏浚技术装备国家工程研究中心有限公司 | 基于力矩补偿的盾构机智能掘进油压均衡分配方法 |
Non-Patent Citations (2)
Title |
---|
不同水平偏心距下群桩内力变化规律试验研究;陈仁朋;王书行;孔令刚;陈云敏;顾明;袁廉华;;岩石力学与工程学报(03);全文 * |
盾构推进系统分区性能分析的等效机构建模;邓颖聪;郭为忠;高峰;;机械工程学报(13);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN117094253A (zh) | 2023-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN117094253B (zh) | 推拼同步下基于推力均匀性的油缸推力分配方法及系统 | |
JP5584174B2 (ja) | シールド掘進機における推進ジャッキの推力設定方法 | |
CN111810174B (zh) | 推拼同步模式下盾构推进系统顶力分配计算方法 | |
CN117150826B (zh) | 适用于同步推拼系统的油缸推力精确分配方法及系统 | |
CN109578000B (zh) | 一种用于盾构推进系统的等比级数布局方法 | |
CN113431592B (zh) | 一种盾构控制方法及其系统 | |
CN113738390B (zh) | 一种矩形顶管机的分时迁移空间姿态微调纠偏方法 | |
CN116882129B (zh) | 恒定/非恒定力学约束下盾构推进系统推力分配计算方法 | |
CN106547942A (zh) | 一种外倾顺层节理岩体被动侧向岩石压力的计算方法 | |
CN110185456B (zh) | 一种盾构推进系统快速实现抗偏载的方法 | |
CN115478864A (zh) | 一种推拼同步模式下盾构推进系统推进力重分配计算方法 | |
CN106285711A (zh) | 一种基于光纤陀螺的顶管纠偏控制系统及其方法 | |
MY137228A (en) | Parallel-kinematical machine with an active measuring system | |
CN116856943A (zh) | 基于共面约束的四分区盾构推力及油压智能均衡分配方法 | |
CN116856944A (zh) | 基于力矩补偿的盾构机智能掘进油压均衡分配方法 | |
CN104732020A (zh) | 千吨级全地面起重机用六排滚子转盘轴承的校核方法 | |
CN115788459A (zh) | 基于掘进拼装同步的盾构机推进系统及推力分配方法 | |
JP3104955B2 (ja) | 同時掘進シールド工法 | |
JP2941671B2 (ja) | 中折れ型トンネル掘削機の方向制御方法 | |
CN113177270B (zh) | 一种盾构姿态引起的地层变形实时计算方法 | |
CN110273689B (zh) | 一种实现矩形盾构抗偏载的推进液压缸调节方法 | |
CN107435543A (zh) | Ⅳ级围岩超大断面三台阶五步施工方法 | |
CN116263109A (zh) | 一种盾构机管片同步拼装控制方法及盾构机 | |
JP3314366B2 (ja) | シールド掘進機における推進ジャッキの圧力設定方法 | |
CN220203896U (zh) | 一种用于36分度管片盾构机推进油缸的分区结构 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |