CN114706322B - 一种盾构机姿态自动控制仿真系统 - Google Patents

一种盾构机姿态自动控制仿真系统 Download PDF

Info

Publication number
CN114706322B
CN114706322B CN202210331663.6A CN202210331663A CN114706322B CN 114706322 B CN114706322 B CN 114706322B CN 202210331663 A CN202210331663 A CN 202210331663A CN 114706322 B CN114706322 B CN 114706322B
Authority
CN
China
Prior art keywords
shield
module
tail
soil
coordinate system
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
Application number
CN202210331663.6A
Other languages
English (en)
Other versions
CN114706322A (zh
Inventor
黄德青
宋晨健
秦娜
徐进
严一舟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202210331663.6A priority Critical patent/CN114706322B/zh
Publication of CN114706322A publication Critical patent/CN114706322A/zh
Application granted granted Critical
Publication of CN114706322B publication Critical patent/CN114706322B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B17/00—Systems involving the use of models or simulators of said systems
    • G05B17/02—Systems involving the use of models or simulators of said systems electric
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • 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

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Excavating Of Shafts Or Tunnels (AREA)

Abstract

本发明公开了一种盾构机姿态自动控制仿真系统，通过构建数字通信接口模块、推进动力学模块、盾体运动学模块、尾盾动力学模块、基于盾体‑土体动态作用模型的土体环境模块、虚拟导向系统模块以及盾构三维场景展示模块，组建盾构机姿态自动控制仿真系统；本发明构建基于盾体‑土体动态作用模型的土体环境模块，并结合各模块间的数据传输合理反映盾构机运动过程盾体所受阻力，阻力矩的情况；并通过数字通信接口模块预留外部的数据/控制接口，通过外部控制信号完善对基于盾构姿态自动控制的仿真测试，为盾构姿态自动控制算法的设计与验证提供了一个低成本，高效率，高可靠性的实现方案，并实时可视化展示盾构机的运动过程。

Description

一种盾构机姿态自动控制仿真系统
技术领域
本发明涉及盾构施工智能化建设技术领域，具体涉及一种盾构机姿态自动控制仿真系统。
背景技术
目前国内各大城市中，地铁建设已经成为市政交通建设的关键部分，现有的地铁隧道挖掘方案多采用盾构机来施工，然而由于城市地下具有错综复杂的建筑基础以及复杂的管线，施工面临的地质条件错综复杂，这就使得盾构机司机具备丰富的施工经验以及较高的专业素养，为了培养出专业的盾构司机，需要研发盾构机模拟驾驶的操作平台用于进行盾构驾驶教学；并用于模拟产生环境对盾构机产生的作用效果，使仿真更具有实际效果和意义。
发明内容
针对现有技术中的上述不足，本发明提供了一种盾构机姿态自动控制仿真系统，通过构建基于盾构机数学模型的盾构姿态自动控制仿真系统，解决了现有盾构姿态自动控制算法验证成本高，耗时长的问题，并通过构建基于盾体-土体动态作用模型的基于盾体-土体动态作用模型的土体环境模块，结合各模块间的数据传输合理反映盾构机运动过程盾体所受阻力，阻力矩的情况。
为了达到上述发明目的，本发明采用的技术方案为：
一方面，一种盾构机姿态自动控制仿真系统，包括：数字通信接口模块、推进动力学模块、盾体运动学模块、尾盾动力学模块、基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块以及盾构三维场景展示模块；
其中，数字通信接口模块分别与推进动力学模块、尾盾动力学模块、盾体运动学模块、基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块以及盾构三维场景展示模块连接；用于接收并存储推进动力学模块、尾盾动力学模块、盾体运动学模块、基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块的输出数据；并为各模块提供基础数据支撑；
推进动力学模块与盾体运动学模块连接；用于根据数字通信接口模块传输的数据计算盾构机油缸的推进行程，并将推进行程传输到盾体运动学模块；
尾盾动力学模块与盾体运动学模块连接；用于接收数字通信接口模块传输数据，并根据传输数据计算铰链行程；并将铰链行程传输到盾体运动学模块；
盾体运动学模块与基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块连接；用于接收推进行程与铰链行程，并根据推进行程与铰链行程计算盾体实时位姿，并将盾体实时位姿传输给基于盾体-土体动态作用模型的土体环境模块与虚拟导向系统模块；
基于盾体-土体动态作用模型的土体环境模块，用于根据盾体实时位姿计算土体环境对推进油缸的推进阻力以及铰接拉力；
虚拟导向系统模块，用于盾体实时位姿计算盾体姿态偏差；
盾构三维场景展示模块，用于根据推数字通信接口模块的输出参数，进行三维建模，并对盾构机运动过程进行可视化展示。
本发明具有以下有益效果：
构建数字通信接口模块、推进动力学模块、盾体运动学模块、尾盾动力学模块、基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块以及盾构三维场景展示模块，组建盾构机姿态自动控制仿真系统；本发明基于盾体-土体间的动态作用构建基于盾体-土体动态作用模型的土体环境模块，并结合各模块间的数据传输合理反映盾构机运动过程盾体所受阻力，阻力矩的情况；并通过数字通信接口模块预留外部的数据/控制接口，通过外部控制信号完善对基于盾构姿态自动控制的仿真测试，为盾构姿态自动控制算法的设计与验证提供了一个低成本，高效率，高可靠性的实现方案，并实时可视化展示盾构机的运动过程。
附图说明
图1为本发明实施例提供的一种盾构机姿态自动控制仿真系统的结构示意图；
图2为本发明实施例中基于推进阻力的盾体-土体动态作用关系示意图；
图3为本发明实施例中基于铰链拉力的盾体-土体动态作用关系示意图；
图4为本发明实施例中盾体运动学模块的结构示意图；
图5为本发明实施例中前中盾运动学子模块的结构示意图；
图6为本发明实施例中尾盾运动学子模块的结构示意图；
图7为本发明实施例中基于盾体-土体动态作用模型的土体环境模块的结构示意图；
图8为本发明实施例中铰链拉力计算子模块的结构示意图；
图9为本发明实施例中推进阻力计算子模块的结构示意图；
图10为本发明实施例中盾构机三维可视化线性框图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
如图1所示，本发明实施例提供一种盾构机姿态自动控制仿真系统，包括：数字通信接口模块、推进动力学模块、盾体运动学模块、尾盾动力学模块、基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块以及盾构三维场景展示模块；
其中，数字通信接口模块分别与推进动力学模块、尾盾动力学模块、盾体运动学模块、基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块以及盾构三维场景展示模块连接；用于接收并存储推进动力学模块、尾盾动力学模块、盾体运动学模块、基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块的输出数据；并为各模块提供基础数据支撑；
本发明实施例中，数字通信接口模块为盾构机姿态自动控制仿真系统提供与外部程序的通信途径，外部程序通过向数字通信接口模块提供基础数据，来实现对盾构机的控制或仿真参数的设置，其中，外部程序为本发明实施例提供给的盾构机姿态自动控制仿真系统提供包括：推进启动、推进模式选择、整体推进速度给定百分比、压力给定百分比、推进速度给定百分比等数据，并存于数字通信接口模块的数据/指令存储区，其存储数据详如表1-4所示；其中各类型数据均对应一个固定的存储地址。
表1为盾构机系统存储区的控制输入数据
Figure BDA0003573289440000041
Figure BDA0003573289440000051
表2土体环境参数存储区的土地模型参数
土体模型参数类型 状态输出数据类型
地基基床系数 环境阻力
盾构机挤压土体所受到的粘滞系数 铰接拉力
等效的摩擦系数 预留的状态输出存储区
水土压压强相关参数
推进阻力系数
表3虚拟导向系统存储区的状态输出数据类型
状态输出数据类型
前盾趋向水平/垂直偏差
盾尾趋向水平/垂直偏差
盾首中心水平/垂直偏差
盾中中心水平/垂直偏差
盾尾中心水平/垂直偏差
当前推进里程
当前推进环号
前中盾位姿
尾盾位姿
表4盾构三维场景展示存储区
Figure BDA0003573289440000052
Figure BDA0003573289440000061
推进动力学模块与盾体运动学模块连接；用于根据数字通信接口模块传输的数据计算盾构机油缸的推进行程，并将推进行程传输到盾体运动学模块；
本发明实施例中，推进动力学模块根据数据/指令存储区中的压力指令和速度指令，以及基于盾体-土体动态作用模型的土体环境模块计算所得的推进阻力(即虚拟压力传感器测得的压力)，结合推进动力学方程，即可求解出虚拟盾构机油缸的推进行程，推进速度，最后将油缸的推进行程输入到盾体运动学模块中前中盾运动学方程，并保存到数据/指令存储区。
尾盾动力学模块与盾体运动学模块连接；用于接收数字通信接口模块传输数据，并根据传输数据计算铰链行程；并将铰链行程传输到盾体运动学模块；
本实施例中，尾盾动力学模块根据数据/指令存储区中的铰接状态输入信号(保持、释放、拖拉)，以及基于盾体-土体动态作用模型的土体环境模块计算所得的盾尾拉力，结合盾尾动力学方程，即可求解出实时铰接行程，最后将铰接行程输入到盾体运动学模型中尾盾运动学方程，并保存到数据/指令存储区。
盾体运动学模块与基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块连接；用于接收推进行程与铰链行程，并根据推进行程与铰链行程计算盾体实时位姿，并将盾体实时位姿传输给基于盾体-土体动态作用模型的土体环境模块与虚拟导向系统模块；
本发明实施例中，盾体运动学模块根据推进动力学模块输出的油缸推进行程，结合前中盾运动学方程，即可求解出虚拟盾构机前中盾盾体在城市坐标系下的实时位姿Pf；同样的，根据尾盾动力学模块输出的铰接行程，结合尾盾运动学方程，即可求解出虚拟盾构机尾盾在城市坐标系下的实时位姿Pe；前中盾位姿和尾盾位姿共同被称为盾体实时位姿，最后盾体实时位姿将被保存到数据/指令存储区，并被输入到虚拟导向系统模块用于计算盾体各个标志点相对于期望隧道中心轴线的姿态偏差，同时也被输入到基于盾体-土体动态作用模型的土体环境模块用于计算土体环境对推进油缸和铰接油缸分别产生的推进阻力与盾尾拉力。
基于盾体-土体动态作用模型的土体环境模块，用于根据盾体实时位姿计算土体环境对推进油缸的推进阻力以及铰接拉力；
本发明实施例中，基于盾体-土体动态作用模型的土体环境模块的作用是根据盾体运动学模块计算出的盾体实时位姿以及推进油缸的推进速度，结合盾体-土体动态作用模型，计算土体环境对推进油缸的推进阻力以及铰接拉力。
本发明中的盾体-土体动态作用模型在静态温克尔弹性地基模型的基础上，将原本弹簧地基模型改进为弹簧阻尼地基模型，从而有效反映盾体-土体之间的作用关系，其中温克尔弹性地基模型给出了地基反力与地基压缩量之间的关系，其将地基等效为多个独立且互不影响的弹簧，该模型适合于对盾体进行静态的受力分析；但是针对盾构推进时的动态情况，土体存在一定的流动性。
本实施例中，计算土体环境对推进油缸的推进阻力的原理，如图2所示，为水平方向和垂直方向转弯或纠偏过程中土体与盾体的作用力示意图，其中盾体轴向所受到的总的环境阻力为Fss，盾体受到的总的环境力矩为τsr，盾体旋转中心相对于前中盾坐标系的位置为CR，重心相对于前中盾坐标系的位置为CG；对于垂直方向，由于刀盘上下部水土压分布不均匀，压力差会对盾体产生一个阻力矩τsp，除此之外，盾构机自身的重力也会产生一个相应的力矩τG，结合盾体轴向所受到的总的环境阻力为Fss，盾体受到的总的环境力矩为τsr，结合力和力矩平衡方程，可以得到各组油缸所受到的环境阻力Fet；
盾尾通过铰接油缸与中盾相连，盾尾在摩擦力的作用下，会对铰接油缸产生相应的拉力，从而拉动铰接油缸伸出缸壳；因此本发明实施例中计算土体环境对推进油缸的铰接拉力的原理，如图3所示，为尾盾-环境-中盾的作用关系，其中Fess为环境对盾尾的拉力，包含盾尾与管片，土体之间的摩擦力等，τesr为环境对盾尾的扭矩，盾尾对应一个旋转中心CER；通过Fess和τesr即可计算出环境对铰接油缸拉力。
虚拟导向系统模块，用于根据盾体实时位姿计算盾体姿态偏差。
本发明实施例中，虚拟导向系统模块的作用是根据盾体运动学模块计算出的盾体实时位姿，以及仿真系统中建立的盾构施工任务信息(即盾构机基础信息)，通过几何关系的求解，计算得到与真实盾构导向系统物理含义一致的数据(包括盾首偏差，盾尾偏差等)，并将计算出的数据存储到数据/指令存储区，具体的输出数据类型如表3所示。
盾构三维场景展示模块，用于根据数字通信接口模块的输出参数，进行三维建模，并对盾构机运动过程进行可视化展示。
本发明实施例中，盾构三维场景展示模块的作用是对数据/指令存储区中虚拟盾构机盾体实时位姿，推进油缸的推进行程，推进速度，铰接油缸的行程，施工场景信息进行可视化展示。其技术原理如下：通过solidworks对盾构机进行三维建模，然后导出盾构机各个模块的stl文件，在开源机器人仿真环境中导入各个机构的三维模型并创建各个推进油缸的移动关节，从而实现虚拟盾构机的可视化；结合盾构机运动学模型与动力学模型得到的盾构机位姿与推进油缸的相关数据，并可实时的将盾构机的运动过程进行可视化展示。
优选地，推进动力学模块具体为：
接收数字通信接口模块传输的压力指令、速度指令及推进阻力，利用推进动力学方程计算盾构机油缸的推进行程，并将推进行程传输到盾体运动学模块；其推进动力学方程表示为：
Figure BDA0003573289440000091
其中，Mc为活塞盾体中质量矩阵，Bc为盾体与土体负载之间的粘性阻尼系数，Kc为活塞盾体中负载的刚度系数，Xc为盾构机油缸的实际行程，即为盾构机油缸的推进行程，FL为推进阻力，Ac为盾构机油缸中各分区总的油缸横截面，Pc为盾构机油缸的油压。
其中，推进行程的计算式表示为：
Figure BDA0003573289440000092
其中，Ts表示单步计算的时间周期，Xc[t]表示当前t时刻的推进行程，Vc[t-Ts]表示t-Ts时刻油缸的伸出速度，Ac[t-Ts]表示t-Ts时刻油缸伸出加速度，Xc[t-Ts]为t-Ts时刻推进油缸行程。
优选地，尾盾动力学模块具体为：
接收数字通信接口模块传输的铰链状态输入信号及盾尾拉力，并结合尾盾动力学模型计算实时铰链行程；并将铰链行程传输到盾体运动学模块。
其中，尾盾动力学模型计算式表示为：
Figure BDA0003573289440000101
其中，Xj为铰接油缸的行程，Fsj为铰接油缸的拉力，Mj为盾尾等效到铰接油缸的等效质量矩阵，Bj为铰接油缸的阻尼系数，Kj为铰接油缸的弹簧系数，Fe为盾尾拉力，Xj为铰接油缸的位移向量，当铰接油缸为保持模式时，Xj为铰接油缸的实际行程；
本发明实施例中，铰接行程表示为：
Figure BDA0003573289440000102
其中，Ts表示单步计算的时间周期，Xj[t]表示当前t时刻的铰接行程，Vc[t-Ts]表示t-Ts时刻油缸的伸出速度，Aj[t-Ts]表示t-Ts油缸伸出加速度。
如图4所示，优选地，盾体运动学模块具体：
包括前中盾运动学子模块与尾盾运动学子模块；
其中，前中盾运动学子模块，用于根据推进动力学模块输出的盾构机油缸的推进行程计算前中盾实时位姿；
如图5所示，优选地，前中盾运动学子模块具体包括：
盾体刀盘中心位置向量计算单元，用于根据推进动力学模块输出的盾构机油缸的推进行程计算当前环号对应的管片坐标系下盾体刀盘中心位置向量，其计算式表示为：
Figure BDA0003573289440000111
其中，BNPc为当前环号对应的管片坐标系下盾体刀盘中心位置向量，h为前中盾盾体的长度，xc为上下左右四组推进油缸的平均行程，θ、ψ分别为盾体相对于当前环号对应的管片坐标系的水平偏转角和垂直偏转角，满足：
Figure BDA0003573289440000112
xci,＝1,…,4分别为右下左上四组油缸的行程,lAC,lBD分别为左右两组油缸距离和上下两组油缸的距离；sin(.)为正弦函数，cos(.)为余弦函数；
位移向量计算单元，用于根据当前环境坐标系下盾体刀盘中心位置向量计算刀盘中心在城市坐标系下的位移向量，其计算式表示为：
Figure BDA0003573289440000113
其中，WPc为刀盘中心在城市坐标系下的位移向量，
Figure BDA0003573289440000114
为当前环境坐标系相对于城市坐标系的旋转矩阵，WPBN为当前环号对应的管片坐标系原点在城市坐标系位置向量；
第一旋转参数计算单元，用于计算盾体相对于城市坐标系的旋转参数，其计算式表示为：
Figure BDA0003573289440000115
其中，
Figure BDA0003573289440000116
为盾体相对于城市坐标系的旋转参数，/>
Figure BDA0003573289440000117
为当前环号对应的管片坐标系相对于城市坐标系的旋转矩阵，/>
Figure BDA0003573289440000121
为盾体相对于当前环号对应的管片坐标系的实时姿态，满足：/>
Figure BDA0003573289440000122
旋转四元数计算单元，用于根据旋转参数构建盾体相对于城市坐标系的旋转四元数，其计算式表示为：
Figure BDA0003573289440000123
其中，Wqf为旋转四元数，dcm2quat(.)为旋转矩阵转四元数函数；
前中盾实时位姿计算单元，用于根据旋转四元数与位移向量计算前中盾实时位姿，其计算式表示为：
Figure BDA0003573289440000124
其中，Pf为前中盾实时位姿。
尾盾运动学子模块，用于根据尾盾动力学模块输出的实时铰链行程计算尾盾在城市坐标系下的实时位姿。
如图6所示，优选地，尾盾运动学子模块具体包括：
第一位置向量计算单元，用于计算盾尾在前中盾坐标系下的位置向量，其计算式表示为：
Figure BDA0003573289440000125
其中，CPE为盾尾在前中盾坐标系下的位置向量，xji,i＝1,…,4分别表示右下左上四组铰接油缸的行程，θj1,ψj1分别为尾盾相对于前中盾的水平偏转角和垂直偏转角，he为尾盾盾体的长度；
第二位置向量计算单元，用于根据盾尾在前中盾坐标系下的位置向量计算盾尾中心相对于城市坐标系的位置向量，其计算式表示为：
Figure BDA0003573289440000131
其中，WPE为盾尾中心相对于城市坐标系的位置向量，WPc为刀盘中心在城市坐标系下的位移向量，
Figure BDA0003573289440000132
为前中盾相对于城市坐标系的旋转矩阵，CPE为盾尾在前中盾坐标系下的位置向量；
实时位姿计算单元，用于计算尾盾相对于前中盾坐标系的实时位姿，其计算式表示为：
Figure BDA0003573289440000133
其中，
Figure BDA0003573289440000134
为尾盾相对于前中盾坐标系的实时位姿，θj1,ψj1分别表示尾盾相对于前中盾的水平偏转角和垂直偏转角，满足：/>
Figure BDA0003573289440000135
Figure BDA0003573289440000136
其中xji,i＝1,…,4分别表示右下左上四组铰接油缸的行程,lACj,lBDj分别表示左右两组铰接油缸距离和上下两组铰接油缸的距离
第二旋转参数计算单元，用于根据尾盾相对于前中盾坐标系的实时位姿计算尾盾相对于城市坐标系的旋转参数，其计算式表示为：
Figure BDA0003573289440000141
其中，
Figure BDA0003573289440000142
为尾盾相对于城市坐标系的旋转参数；
第三旋转参数计算单元，用于根据尾盾相对于城市坐标系的旋转参数计算城市坐标系到盾尾坐标系的旋转参数，其计算式表示为：
Figure BDA0003573289440000143
其中，Wqe为城市坐标系到盾尾坐标系的旋转参数，dcm2quat(.)为旋转矩阵转四元数函数；
尾盾实时位姿计算单元，用于根据城市坐标系到盾尾坐标系的旋转参数与盾尾中心相对于城市坐标系的位置向量计算尾盾在城市坐标系下的实时位姿，其计算式表示为：
Figure BDA0003573289440000144
其中，Pe为尾盾在城市坐标系下的实时位姿。
如图7所示，优选地，基于盾体-土体动态作用模型的土体环境模块具体包括：
铰链拉力计算子模块，用于根据盾体运动学模块输出的盾体实时位姿结合盾体-土体动态作用模型计算铰接拉力；
如图8所示，优选地，铰链拉力计算子模块具体包括：
土体抗力力矩计算单元，用于利用盾体-土地动态作用模型根据盾体运动学模块输出的盾体实时位姿，计算盾尾的土体抗力力矩，其计算式表示为：
Figure BDA0003573289440000151
其中，τesr为盾尾的土体抗力力矩，se2为盾尾头部中点里程，Zj(s),
Figure BDA0003573289440000152
分别为盾尾压缩土体量和压缩土体量的变化量，K为地基基床系数矩阵，D为盾构机挤压土体所受到的粘滞系数，max(.)为最值函数；
轴向阻力计算单元，用于根据推进油缸的推进速度计算盾尾受到的轴向阻力，其计算式表示为：
Figure BDA0003573289440000153
其中，Fess为盾尾受到的轴向阻力，N为总的油缸分组个数，Vt为推进油缸的推进速度，μ为摩擦系数，F2为尾盾侧边所受到的水土压力；
环境拉力计算单元，用于根据盾尾的土体抗力力矩构建第二扭矩平衡方程，计算各组铰接油缸所受环境拉力，其第二扭矩平衡方程表示为：
Figure BDA0003573289440000154
其中，Fej1为尾盾所受力矩折算到各组铰接油缸上面的拉力值，作为各组铰接油缸所受环境拉力，Lej1为各组铰接油缸对应的力臂；
铰链拉力计算单元，用于根据各组铰接油缸所受环境拉力与轴向阻力计算铰链拉力，其计算式表示为：
Fej＝Fej1+Fess
其中，Fej为铰链拉力。
推进阻力计算子模块，用于根据盾体运动学模块输出的盾体实时位姿结合盾体-土体动态作用模型计算推进阻力；
如图9所示，优选地，推进阻力计算子模块具体包括：
刀盘中心历史位置轨迹计算单元，用于构建刀盘中心历史轨迹的插值函数，并结合历史盾体实时位姿计算刀盘行进轨迹点，得到前中盾坐标系下刀盘中心历史位置轨迹，其中前中盾坐标系下刀盘中心历史位置轨迹的计算式表示为：
Figure BDA0003573289440000161
其中，CXct(s)为前中盾坐标系下刀盘历史位置轨迹，即表示盾构机各个方向投影面的面积，se1为盾尾尾部中点里程，Xct(s)为插值后的刀盘位置函数，
Figure BDA0003573289440000162
为盾体相对于世界坐标系的旋转矩阵，(.)-1为矩阵求逆函数；
本发明实施例中，插值后的刀盘位置函数Xct(s)表示为：
Figure BDA0003573289440000163
其中，xct(s)为刀盘中心的X坐标值关于刀盘里程的函数，yct(s)为刀盘中心的Y坐标值关于刀盘里程的函数，zct(s)为刀盘中心的Z坐标值关于刀盘里程的函数，
Figure BDA0003573289440000165
为当前刀盘里程。
土体压缩量计算单元，用于根据前中盾坐标系下刀盘中心历史位置轨迹计算盾体里程处对土体压缩量，其中盾体里程处对土体的压缩量的计算式表示为：
Figure BDA0003573289440000164
其中，Z(s)为盾体里程S处对土体的压缩量，z1(s)为水平方向上土体压缩量，z2(s)为垂直方向上土体压缩量，
Figure BDA0003573289440000171
分别为前中盾坐标系下历史位置轨迹/>
Figure BDA0003573289440000172
在水平方向与垂直方向上坐标值，R为前中盾的外壳半径，sgn(.)为阶跃函数，arccos(.)为反余弦函数；
压缩变化量计算单元，用于根据盾体里程对土体的压缩量计算压缩变化量，计算式表示为：
Figure BDA0003573289440000173
其中，
Figure BDA0003573289440000174
为压缩变化量；
盾体-土地动态系统单元，用于利用温克尔弹性地基原理，根据压缩量及压缩变化量计算推进过程中受到侧边土体抗力产生的力矩，得到侧边抗力力矩，其计算式表示为：
Figure BDA0003573289440000175
其中，τsr为推进过程中受到侧边土体抗力产生的力矩，作为侧边抗力力矩，C(.)为单位转换矩阵，即为盾构机各个方向上的投影面的面积，K为地基基床系数矩阵，D为盾构机挤压土体所受到的粘滞系数；s0为中盾尾部中心的里程，sc为刀盘中心的里程；L(s)为前中盾上里程为s的各个受力点相对于等效旋转中心的力臂；
本发明实施例中，温克尔弹性地基模型，具体为通过：假设地基表面任一点的压力强度与该点的沉降成正比，即p＝ks；式中p为地基表面某点单位面积上的压力；s为相应点的竖向位移；k为地基反力系数，又称基床系数。
重力扭矩计算单元，用于根据预设条件计算重力扭矩，其计算式表示为：
τG＝Gs·LG
其中，τG为重力扭矩，Gs为前中盾与刀盘等效的重力中心。LG为盾体重心相对于转动中心的力臂；
本发明实施例中，在假设盾体的中心位于CG处，且盾体的转动中心和盾体的重心不重合的条件下计算计算重力扭矩。
等效阻力计算单元，用于根据重力扭矩与侧边抗力力矩构建第一扭矩平衡方程，得到总环境扭矩对各推进油缸的等效阻力，其第一扭矩平衡方程表示为：
Figure BDA0003573289440000181
其中，Fet1,i为盾体所受阻力矩折算到第i组推进油缸上面的阻力值，作为总环境扭矩对推进油缸的等效阻力，Let1,i为第i组推进油缸对应的力臂，N为推进油缸的分组个数；
盾体轴向所受阻力计算单元，用于计算盾体轴向所受阻力，其计算式表示为：
Fss＝Fe1+Fe2+Fe3
Figure BDA0003573289440000182
Fe2＝∫ΔSPrdS/N
Figure BDA0003573289440000183
其中，Fss为盾体轴向所受阻力，Fe1为推进过程中受到的摩擦力，包含牵引力，μ为摩擦系数，F1为盾体所受到水土压等压力的总和，Fe2为静止水土压对各组油缸产生的轴向作用力，ΔS为刀盘面面积微元，Pr为刀盘面所受水土压的压强函数，Fe3代表在盾体推进过程到盘面所受前方土体产生的动态阻力，Vt为推进油缸的推进速度，α为阻力系数；
推进阻力计算单元，用于利用盾体轴向所受阻力及总环境扭矩对各推进油缸的等效阻力，计算推进阻力，其计算式表示为：
Fet＝Fss+Fet1
其中，Fet为推进阻力，Fet1为总环境扭矩对各推进油缸的等效阻力。
本发明实施例中，虚拟导向系统模块的计算过程包括以下步骤：
A1、利用线性插值法根据隧道期望中心轴线的勘测坐标，计算隧道期望中心轴线的插值函数；
A2、利用隧道期望中心轴线的插值函数分别计算当前刀盘里程与盾尾里程；计算式表示为：
Figure BDA0003573289440000191
Figure BDA0003573289440000192
其中，sc,se分别为当前刀盘里程与盾尾里程，WpE为盾尾中心位置，WpC为刀盘中心位置，Xdta(s)为隧道期望中心轴线的插值函数，
Figure BDA0003573289440000195
为最小值计算函数，||.||为2范数；
A3、利用隧道期望中心轴线的插值函数分别前中盾坐标系、盾尾坐标系下期望轴线的函数值，计算式表示为：
Figure BDA0003573289440000193
Figure BDA0003573289440000194
其中，CXcdta(s)、EXedta(s)分别为前中盾坐标系、盾尾坐标系下期望轴线的函数值，
Figure BDA0003573289440000201
为盾体相对于世界坐标系的旋转矩阵，/>
Figure BDA0003573289440000202
为尾盾相对于城市坐标系的旋转矩阵，(.)-1为矩阵求逆函数；
A4、根据当前刀盘里程与盾尾里程与前中盾坐标系、盾尾坐标系下期望轴线的函数值计算盾体姿态偏差。
其中盾体姿态偏差包括盾首偏差，盾首趋向偏差，盾尾偏差，盾尾趋向偏差，其计算式分别表示为：
Figure BDA0003573289440000203
Figure BDA0003573289440000204
Figure BDA0003573289440000205
Figure BDA0003573289440000206
其中，Ec、Etc、Ee及Ete分别表示为盾首偏差、盾首趋向偏差、盾尾偏差、盾尾趋向偏差，
Figure BDA0003573289440000207
为偏差转换矩阵，满足/>
Figure BDA0003573289440000208
本发明实施例中可视化场景的线性图如图10所示，其中只包含盾体姿态控制相关的可视化组件，包括：刀盘、前中盾、铰接油缸、推进油缸、尾盾、管片以及期望隧道中心轴线。在仿真系统运行时，推进油缸，铰接油缸会对每一时刻数学模型计算出的行程同步显示，除此之外，盾体整体的位姿的变化情况也会同步的通过该三维场景进行展示；本发明实施例为盾构姿态自动控制算法的设计与验证提供了一个低成本，高效率，高可靠性的实现方案。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

Claims (9)

1.一种盾构机姿态自动控制仿真系统，其特征在于，包括：数字通信接口模块、推进动力学模块、盾体运动学模块、尾盾动力学模块、基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块以及盾构三维场景展示模块；
其中，数字通信接口模块分别与推进动力学模块、尾盾动力学模块、盾体运动学模块、基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块以及盾构三维场景展示模块连接；用于接收并存储推进动力学模块、尾盾动力学模块、盾体运动学模块、基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块的输出数据；并为各模块提供基础数据支撑；
推进动力学模块与盾体运动学模块连接；用于根据数字通信接口模块传输的数据计算盾构机油缸的推进行程，并将推进行程传输到盾体运动学模块；
尾盾动力学模块与盾体运动学模块连接；用于接收数字通信接口模块传输数据，并根据传输数据计算铰链行程；并将铰链行程传输到盾体运动学模块；
盾体运动学模块与基于盾体-土体动态作用模型的土体环境模块、虚拟导向系统模块连接；用于接收推进行程与铰链行程，并根据推进行程与铰链行程计算盾体实时位姿，并将盾体实时位姿传输给基于盾体-土体动态作用模型的土体环境模块与虚拟导向系统模块；
基于盾体-土体动态作用模型的土体环境模块，用于根据盾体实时位姿计算土体环境对推进油缸的推进阻力以及铰接拉力；
虚拟导向系统模块，用于根据盾体实时位姿计算盾体姿态偏差；
盾构三维场景展示模块，用于根据数字通信接口模块的输出参数，进行三维建模，并对盾构机运动过程进行可视化展示。
2.根据权利要求1所述的盾构机姿态自动控制仿真系统，其特征在于，推进动力学模块具体为：
接收数字通信接口模块传输的压力指令、速度指令及推进阻力，利用推进动力学方程计算盾构机油缸的推进行程，并将推进行程传输到盾体运动学模块；其推进动力学方程表示为：
Figure QLYQS_1
其中，Mc为活塞盾体中质量矩阵，Bc为盾体与土体负载之间的粘性阻尼系数，Kc为活塞盾体中负载的刚度系数，Xc为盾构机油缸的实际行程，即为盾构机油缸的推进行程，FL为推进阻力，Ac为盾构机油缸中各分区总的油缸横截面，Pc为盾构机油缸的油压。
3.根据权利要求1所述的盾构机姿态自动控制仿真系统，其特征在于，尾盾动力学模块具体为：
接收数字通信接口模块传输的铰链状态输入信号及盾尾拉力，并结合尾盾动力学模型计算实时铰链行程；并将铰链行程传输到盾体运动学模块；
其中，尾盾动力学模型计算式表示为：
Figure QLYQS_2
其中，Xj为铰接油缸的行程，Fsj为铰接油缸的拉力，Mj为盾尾等效到铰接油缸的等效质量矩阵，Bj为铰接油缸的阻尼系数，Kj为铰接油缸的弹簧系数，Fe为盾尾拉力。
4.根据权利要求1所述的盾构机姿态自动控制仿真系统，其特征在于，盾体运动学模块具体：
包括前中盾运动学子模块与尾盾运动学子模块；
其中，前中盾运动学子模块，用于根据推进动力学模块输出的盾构机油缸的推进行程计算前中盾实时位姿；
尾盾运动学子模块，用于根据尾盾动力学模块输出的实时铰链行程计算尾盾在城市坐标系下的实时位姿。
5.根据权利要求4所述的盾构机姿态自动控制仿真系统，其特征在于，前中盾运动学子模块具体包括：
盾体刀盘中心位置向量计算单元，用于根据推进动力学模块输出的盾构机油缸的推进行程计算当前环号对应的管片坐标系盾体刀盘中心位置向量，其计算式表示为：
Figure QLYQS_3
其中，BNPc为当前环号对应的管片坐标系盾体刀盘中心位置向量，h为前中盾盾体的长度，
Figure QLYQS_4
为上下左右四组推进油缸的平均行程，θ、ψ分别为盾体相对于当前环号对应的管片坐标系的水平偏转角和垂直偏转角，sin(.)为正弦函数，cos(.)为余弦函数；
位移向量计算单元，用于根据当前环境坐标系下盾体刀盘中心位置向量计算刀盘中心在城市坐标系下的位移向量，其计算式表示为：
Figure QLYQS_5
其中，WPc为刀盘中心在城市坐标系下的位移向量，
Figure QLYQS_6
为当前环境坐标系相对于城市坐标系的旋转矩阵，WPBN为当前环号对应的管片坐标系原点在城市坐标系位置向量；
第一旋转参数计算单元，用于计算盾体相对于城市坐标系的旋转参数，其计算式表示为：
Figure QLYQS_7
其中，
Figure QLYQS_8
为盾体相对于城市坐标系的旋转参数，/>
Figure QLYQS_9
为当前环号对应的管片坐标系下相对于城市坐标系的旋转矩阵，/>
Figure QLYQS_10
为盾体相对于当前环号对应的管片坐标系下实时姿态；
旋转四元数计算单元，用于根据旋转参数构建盾体相对于城市坐标系的旋转四元数；
前中盾实时位姿计算单元，用于根据旋转四元数与位移向量计算前中盾实时位姿，其计算式表示为：
Figure QLYQS_11
其中，Pf为前中盾实时位姿。
6.根据权利要求4所述的盾构机姿态自动控制仿真系统，其特征在于，尾盾运动学子模块具体包括：
第一位置向量计算单元，用于计算盾尾在前中盾坐标系下的位置向量；
第二位置向量计算单元，用于根据盾尾在前中盾坐标系下的位置向量计算盾尾中心相对于城市坐标系的位置向量，其计算式表示为：
Figure QLYQS_12
其中，WPE为盾尾中心相对于城市坐标系的位置向量，WPc为刀盘中心在城市坐标系下的位移向量，
Figure QLYQS_13
为前中盾相对于城市坐标系的旋转矩阵，CPE为盾尾在前中盾坐标系下的位置向量；
实时位姿计算单元，用于计算尾盾相对于前中盾坐标系的实时位姿；
第二旋转参数计算单元，用于根据尾盾相对于前中盾坐标系的实时位姿计算尾盾相对于城市坐标系的旋转参数，其计算式表示为：
Figure QLYQS_14
其中，
Figure QLYQS_15
为尾盾相对于城市坐标系的旋转参数；
第三旋转参数计算单元，用于根据尾盾相对于城市坐标系的旋转参数计算城市坐标系到盾尾坐标系的旋转参数；
尾盾实时位姿计算单元，用于根据城市坐标系到盾尾坐标系的旋转参数与盾尾中心相对于城市坐标系的位置向量计算尾盾在城市坐标系下的实时位姿，其计算式表示为：
Figure QLYQS_16
其中，Pe为尾盾在城市坐标系下的实时位姿。
7.根据权利要求1所述的盾构机姿态自动控制仿真系统，其特征在于，基于盾体-土体动态作用模型的土体环境模块具体包括：
铰链拉力计算子模块，用于根据盾体运动学模块输出的盾体实时位姿结合盾体-土体动态作用模型计算铰接拉力；
推进阻力计算子模块，用于根据盾体运动学模块输出的盾体实时位姿结合盾体-土体动态作用模型计算推进阻力。
8.根据权利要求7所述的盾构机姿态自动控制仿真系统，其特征在于，推进阻力计算子模块具体包括：
刀盘中心历史位置轨迹计算单元，用于构建刀盘中心历史轨迹的插值函数，并结合历史盾体实时位姿计算刀盘行进轨迹点，得到前中盾坐标系下刀盘中心历史位置轨迹，其中前中盾坐标系下刀盘中心历史位置轨迹的计算式表示为：
Figure QLYQS_17
其中，CXct(s)为前中盾坐标系下刀盘历史位置轨迹，Xct(s)为插值后的刀盘位置函数，
Figure QLYQS_18
为盾体相对于世界坐标系的旋转矩阵，(.)-1为矩阵求逆运算函数；
土体压缩量计算单元，用于根据前中盾坐标系下刀盘中心历史位置轨迹计算盾体里程处对土体压缩量，其中盾体里程处对土体的压缩量的计算式表示为：
Figure QLYQS_19
其中，Z(s)为盾体里程S处对土体的压缩量，z1(s)为水平方向上土体压缩量，z2(s)为垂直方向上土体压缩量，Cxct(s)，Cyct(s)分别为前中盾坐标系下历史位置轨迹CXct(s)在水平方向与垂直方向上坐标值，R为前中盾的外壳半径，sgn(.)为阶跃函数，arccos(.)为反余弦函数；
压缩变化量计算单元，用于根据盾体里程对土体的压缩量计算压缩变化量，计算式表示为：
Figure QLYQS_20
其中，
Figure QLYQS_21
为压缩变化量；
盾体-土地动态系统单元，用于利用温克尔弹性地基原理，根据压缩量及压缩变化量计算推进过程中受到侧边土体抗力产生的力矩，得到侧边抗力力矩，其计算式表示为：
Figure QLYQS_22
其中，τsr为推进过程中受到侧边土体抗力产生的力矩，作为侧边抗力力矩，C(.)为单位转换矩阵，K为地基基床系数矩阵，D为盾构机挤压土体所受到的粘滞系数；s0为中盾尾部中心的里程，sc为刀盘中心的里程；L(s)为前中盾上里程为s的各个受力点相对于等效旋转中心的力臂；
重力扭矩计算单元，用于根据预设条件计算重力扭矩，其计算式表示为：
τG＝Gs·LG
其中，τG为重力扭矩，Gs为前中盾与刀盘等效的重力中心，LG为盾体重心相对于转动中心的力臂；
等效阻力计算单元，用于根据重力扭矩与侧边抗力力矩构建第一扭矩平衡方程，得到总环境扭矩对各推进油缸的等效阻力，其第一扭矩平衡方程表示为：
Figure QLYQS_23
其中，Fet1,i为盾体所受阻力矩折算到第i组推进油缸上面的阻力值，作为总环境扭矩对推进油缸的等效阻力，Let1,i为第i组推进油缸对应的力臂，N为推进油缸的分组个数；
盾体轴向所受阻力计算单元，用于计算盾体轴向所受阻力，其计算式表示为：
Fss＝Fe1+Fe2+Fe3
Figure QLYQS_24
Fe2＝∫ΔSPrdS/N
Figure QLYQS_25
其中，Fss为盾体轴向所受阻力，Fe1为推进过程中受到的摩擦力，μ为摩擦系数，F1为盾体所受到水土压等压力的总和，Fe2为静止水土压对各组油缸产生的轴向作用力，ΔS为刀盘面面积微元，Pr为刀盘面所受水土压的压强函数，Fe3代表在盾体推进过程到盘面所受前方土体产生的动态阻力，Vt为推进油缸的推进速度，α为阻力系数；
推进阻力计算单元，用于利用盾体轴向所受阻力及总环境扭矩对各推进油缸的等效阻力，计算推进阻力，其计算式表示为：
Fet＝Fss+Fet1
其中，Fet为推进阻力，Fet1为总环境扭矩对各推进油缸的等效阻力。
9.根据权利要求7所述的盾构机姿态自动控制仿真系统，其特征在于，铰链拉力计算子模块具体包括：
土体抗力力矩计算单元，用于利用盾体-土地动态作用模型根据盾体运动学模块输出的盾体实时位姿，计算盾尾的土体抗力力矩，其计算式表示为：
Figure QLYQS_26
其中，τesr为盾尾的土体抗力力矩，se1为盾尾尾部中点里程，se2为盾尾头部中点里程，Zj(s),
Figure QLYQS_27
分别为盾尾压缩土体量和压缩土体量的变化量，K为地基基床系数矩阵，D为盾构机挤压土体所受到的粘滞系数，max(.)为最值函数；
轴向阻力计算单元，用于根据推进油缸的推进速度计算盾尾受到的轴向阻力，其计算式表示为：
Figure QLYQS_28
其中，Fess为盾尾受到的轴向阻力，N为总的油缸分组个数，Vt为推进油缸的推进速度，μ为摩擦系数，F2为尾盾侧边所受到的水土压力；
环境拉力计算单元，用于根据盾尾的土体抗力力矩构建第二扭矩平衡方程，计算各组铰接油缸所受环境拉力，其第二扭矩平衡方程表示为：
Figure QLYQS_29
其中，Fej1为尾盾所受力矩折算到各组铰接油缸上面的拉力值，作为各组铰接油缸所受环境拉力，Lej1为各组铰接油缸对应的力臂；
铰链拉力计算单元，用于根据各组铰接油缸所受环境拉力与轴向阻力计算铰链拉力，其计算式表示为：
Fej＝Fej1+Fess
其中，Fej为铰链拉力。
CN202210331663.6A 2022-03-30 2022-03-30 一种盾构机姿态自动控制仿真系统 Active CN114706322B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210331663.6A CN114706322B (zh) 2022-03-30 2022-03-30 一种盾构机姿态自动控制仿真系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210331663.6A CN114706322B (zh) 2022-03-30 2022-03-30 一种盾构机姿态自动控制仿真系统

Publications (2)

Publication Number Publication Date
CN114706322A CN114706322A (zh) 2022-07-05
CN114706322B true CN114706322B (zh) 2023-06-23

Family

ID=82170415

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210331663.6A Active CN114706322B (zh) 2022-03-30 2022-03-30 一种盾构机姿态自动控制仿真系统

Country Status (1)

Country Link
CN (1) CN114706322B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115562072B (zh) * 2022-10-13 2024-04-26 中煤科工集团上海有限公司 一种带式输送机动态分析系统
CN115982515B (zh) * 2023-01-05 2023-09-29 西南交通大学 一种盾构机姿态控制参数最佳取值的获取方法
CN117552796B (zh) * 2024-01-11 2024-04-16 北京城建设计发展集团股份有限公司 盾构推进系统油缸伸缩速度控制方法、装置、设备及介质

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101211185A (zh) * 2006-12-31 2008-07-02 中铁隧道股份有限公司 盾构机控制系统检测试验台
CN102034005A (zh) * 2010-12-17 2011-04-27 中铁十四局集团有限公司 盾构隧道施工盾构机姿态仿真检测系统
KR20130060478A (ko) * 2011-11-30 2013-06-10 호서대학교 산학협력단 터널 보어링 머신의 굴착 모사 시험 장치
CN103714214A (zh) * 2013-12-30 2014-04-09 上海市机械施工集团有限公司 盾构机组合仿真设备和仿真方法
CN103745648A (zh) * 2013-12-30 2014-04-23 上海市机械施工集团有限公司 盾构机双模型仿真设备与方法
CN106354029A (zh) * 2016-11-10 2017-01-25 西安电子科技大学 盾构机三维虚拟仿真系统及仿真方法
CN106481344A (zh) * 2016-10-25 2017-03-08 华中科技大学 一种土压平衡盾构机隧道掘进参数智能控制方法
CN108166985A (zh) * 2017-12-06 2018-06-15 浙江大学 一种基于驱动轴扭矩估计的盾构机刀盘驱动系统的控制方法
KR20180083471A (ko) * 2017-01-12 2018-07-23 현대건설주식회사 토사지반 tbm 성능 예측을 위한 굴진성능 평가 시험장치 및 시험방법
CN110208485A (zh) * 2019-05-22 2019-09-06 西南交通大学 一种用于泥水盾构泥膜形成及渗透破坏试验的试验装置
CN110486029A (zh) * 2019-09-09 2019-11-22 辽宁石油化工大学 盾构机土压平衡控制的优化方法及装置
CN110633525A (zh) * 2019-09-12 2019-12-31 辽宁石油化工大学 土压平衡盾构机系统能耗的优化方法及装置
CN111179677A (zh) * 2019-12-25 2020-05-19 中交天和机械设备制造有限公司 一种盾构机操作工培训考评系统
CN114019795A (zh) * 2021-10-15 2022-02-08 中铁高新工业股份有限公司 一种基于强化学习的盾构掘进纠偏智能决策方法

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101211185A (zh) * 2006-12-31 2008-07-02 中铁隧道股份有限公司 盾构机控制系统检测试验台
CN102034005A (zh) * 2010-12-17 2011-04-27 中铁十四局集团有限公司 盾构隧道施工盾构机姿态仿真检测系统
KR20130060478A (ko) * 2011-11-30 2013-06-10 호서대학교 산학협력단 터널 보어링 머신의 굴착 모사 시험 장치
CN103714214A (zh) * 2013-12-30 2014-04-09 上海市机械施工集团有限公司 盾构机组合仿真设备和仿真方法
CN103745648A (zh) * 2013-12-30 2014-04-23 上海市机械施工集团有限公司 盾构机双模型仿真设备与方法
CN106481344A (zh) * 2016-10-25 2017-03-08 华中科技大学 一种土压平衡盾构机隧道掘进参数智能控制方法
CN106354029A (zh) * 2016-11-10 2017-01-25 西安电子科技大学 盾构机三维虚拟仿真系统及仿真方法
KR20180083471A (ko) * 2017-01-12 2018-07-23 현대건설주식회사 토사지반 tbm 성능 예측을 위한 굴진성능 평가 시험장치 및 시험방법
CN108166985A (zh) * 2017-12-06 2018-06-15 浙江大学 一种基于驱动轴扭矩估计的盾构机刀盘驱动系统的控制方法
CN110208485A (zh) * 2019-05-22 2019-09-06 西南交通大学 一种用于泥水盾构泥膜形成及渗透破坏试验的试验装置
CN110486029A (zh) * 2019-09-09 2019-11-22 辽宁石油化工大学 盾构机土压平衡控制的优化方法及装置
CN110633525A (zh) * 2019-09-12 2019-12-31 辽宁石油化工大学 土压平衡盾构机系统能耗的优化方法及装置
CN111179677A (zh) * 2019-12-25 2020-05-19 中交天和机械设备制造有限公司 一种盾构机操作工培训考评系统
CN114019795A (zh) * 2021-10-15 2022-02-08 中铁高新工业股份有限公司 一种基于强化学习的盾构掘进纠偏智能决策方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
盾构掘进模拟试验系统及相关技术研究;朱北斗;《中国优秀硕士学位论文全文数据库工程科技II辑》(第7期);全文 *
盾构推进系统的机构-环境一体化建模与推进分区研究;郭文韬;《中国博士学位论文全文数据库工程科技II辑》(第8期);全文 *

Also Published As

Publication number Publication date
CN114706322A (zh) 2022-07-05

Similar Documents

Publication Publication Date Title
CN114706322B (zh) 一种盾构机姿态自动控制仿真系统
CN103745648B (zh) 盾构机双模型仿真设备与方法
CN103303433B (zh) 一种船舶性能虚拟测试系统
CN102426425B (zh) 一种汽车abs虚拟现实仿真系统
CN104484522A (zh) 一种基于现实场景的机器人模拟演练系统的构建方法
CN102561395B (zh) 一种面向沉管隧道抗震设计的三维精细化建模方法
CN101944299B (zh) 全断面掘进机虚拟施工系统
CN103714214A (zh) 盾构机组合仿真设备和仿真方法
CN106493735A (zh) 存在外界扰动的柔性机械臂扰动观测控制方法
CN108205516A (zh) 一种基于条件传递的海洋结构流固耦合动力响应算法
CN101609615A (zh) 一种实时半实物仿真系统及其数据采集方法
CN111797459A (zh) 一种有砟轨道-桥梁动力学耦合模型的构建方法
CN105759632B (zh) 一种水下机器人伺服阀控制液压推进器动力学仿真系统及其仿真方法
CN104200052A (zh) 液压驱动六足机器人仿真系统及仿真方法
CN110261141B (zh) 一种悬浮隧道车-隧动力耦合混合模拟测试方法与装置
CN110263359B (zh) 一种悬浮隧道管道性能混合模拟测试方法与装置
CN104318601A (zh) 一种流体环境下人体运动仿真方法
CN102213646B (zh) 铺管船托管架动力子结构实验装置及实验方法
CN204101270U (zh) 液罐车辆液-固双向耦合实时模拟试验台
CN111179677A (zh) 一种盾构机操作工培训考评系统
CN115310339A (zh) 基于物质点法的具有表面张力效应的固液耦合模拟方法
CN104915481A (zh) 基于虚拟样机建模和周期性规划的球形电机协同控制
Jia et al. A virtual laboratory on fluid mechanics
CN113867175A (zh) 轨道交通模型创建方法、装置、计算机设备和存储介质
RU150806U1 (ru) Динамический тренажер вождения автомобиля

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