CN115575004B - 落石冲击埋地管道的冲击性能计算方法、设备和存储介质 - Google Patents

落石冲击埋地管道的冲击性能计算方法、设备和存储介质 Download PDF

Info

Publication number
CN115575004B
CN115575004B CN202210990958.4A CN202210990958A CN115575004B CN 115575004 B CN115575004 B CN 115575004B CN 202210990958 A CN202210990958 A CN 202210990958A CN 115575004 B CN115575004 B CN 115575004B
Authority
CN
China
Prior art keywords
impact
falling
stone
vertical
maximum
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
CN202210990958.4A
Other languages
English (en)
Other versions
CN115575004A (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 CN202210990958.4A priority Critical patent/CN115575004B/zh
Publication of CN115575004A publication Critical patent/CN115575004A/zh
Application granted granted Critical
Publication of CN115575004B publication Critical patent/CN115575004B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • 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
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Probability & Statistics with Applications (AREA)
  • Devices Affording Protection Of Roads Or Walls For Sound Insulation (AREA)

Abstract

本发明公开了落石冲击埋地管道的冲击性能计算方法、设备和存储介质，以解决现有技术中对斜向落石冲击力计算不准的问题。落石冲击埋地管道的冲击性能计算方法，包括以下步骤：获取落石参数和地质参数，所述落石参数包括落石的球体等效半径、落石质量、落石冲击埋地管道的对地面入射角度与入射速度，所述地质参数包括管道的覆土层厚度、泊松比、土层剪切波速和土层密度；将所述落石参数和地质参数输入至落石冲击计算模型中，获取所述落石冲击计算模型输出的冲击性能数据；所述的冲击性能数据至少包括落石的最大冲击力和冲击动力放大系数，所述冲击动力放大系数为落石最大冲击力与落石重力的比值。本发明的算法概念明确，计算简单，计算参数获取方便。

Description

落石冲击埋地管道的冲击性能计算方法、设备和存储介质
技术领域
本发明涉及落石冲击埋地管道的技术领域，具体而言，涉及落石冲击埋地管道的冲击性能计算方法、设备和存储介质。
背景技术
崩塌落石作为山区常见的一种地质灾害，在我国西部山区尤为频发。当管道沿山体坡脚铺设时，陡崖上的危岩就成为威胁管道安全的主要地质灾害隐患之一。危岩在风化崩解、人为扰动、地震动力、暴雨影响等外在因素作用下可能脱离母岩而启动，在短时间内获得巨大能量，可对埋地管道造成严重的冲击危害，甚者能直接击穿管道上方的覆土层至管道破损。
在实际的管道工程建设中，常会遇到高陡崖壁上的危岩地形。由于危岩运动轨迹难以预测，导致拦石墙、防护网等防治工程结构有时难以发挥实际作用效果，而通常可在管道上方设置一定厚度的覆土缓冲层，以其缓冲耗散能力来减缓落石冲击对管道所造成的危害。因此，确定合理的管道的覆土层厚度是危岩落石区域管道防冲击的关键参数之一，其中，控制性设计环节在于如何合理确定落石冲击力。
对落石冲击力的计算，我国公路路基规范、铁路隧道手册等曾建议了算法。然而，现有算法在实际工程中常存在冲击力计算不准而造成防护结构失效的问题，尤其是落石斜向冲击地面的情况。
发明内容
本发明的主要目的在于提供落石冲击埋地管道的冲击性能计算方法、设备和存储介质，以解决现有技术中对斜向落石冲击力计算不准的问题。
为了实现上述目的，根据本发明的第一个方面，提供了落石冲击埋地管道的冲击性能计算方法，技术方案如下：
落石冲击埋地管道的冲击性能计算方法，包括以下步骤：
获取落石参数和地质参数，所述落石参数包括落石的球体等效半径、落石质量、落石冲击埋地管道的对地面入射角度与入射速度，所述地质参数包括管道的覆土层厚度、泊松比、土层剪切波速和土层密度；
将所述落石参数和地质参数输入至落石冲击计算模型中，获取所述落石冲击计算模型输出的冲击性能数据；
所述的冲击性能数据至少包括落石的最大冲击力和冲击动力放大系数，所述冲击动力放大系数为落石最大冲击力与落石重力的比值。
为了实现上述目的，根据本发明的第二个方面，还提供了落石冲击埋地管道的冲击性能计算设备，技术方案如下：
落石冲击埋地管道的冲击性能计算设备，包括处理器、通信接口、存储器和通信总线，处理器、通信接口和存储器通过通信总线完成相互间的通信；存储器，用于存放计算机程序；处理器，用于执行存储器上所存放的程序时，实现上述第一方面所述的落石冲击埋地管道的冲击性能计算方法。
为了实现上述目的，根据本发明的第三个方面，还提供了计算机可读存储介质，技术方案如下：
计算机可读存储介质，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现上述第一方面所述的落石冲击埋地管道的冲击性能计算方法。
首先，本发明所求得的是落石最大冲击力，而非既有方法的平均冲击力，这不同于公路路基规范、铁路隧道手册及既有文献的相关算法。其次，本发明的算法概念明确，计算简单，计算参数获取方便。本发明在实际工程中的应用便捷，求解得到的冲击动力放大系数，进而可分析埋地管道所受冲击力以及验算管体强度，从而为埋地管道埋深的合理且经济设计提供技术支撑，达到埋地管道落石冲击灾害的合理防治与节约工程成本的目标。
下面结合附图和具体实施方式对本发明做进一步的说明。本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
构成本发明的一部分的附图用来辅助对本发明的理解，附图中所提供的内容及其在本发明中有关的说明可用于解释本发明，但不构成对本发明的不当限定。在附图中：
图1为本发明的冲击性能计算方法中斜向冲击的示意图。
图2为本发明的冲击性能计算方法中竖向冲击下的阻尼-弹簧等效模型的示意图。
图3为本发明的冲击性能计算方法中水平向冲击下的等效弹簧模型的示意图。
图4为本发明的冲击性能计算方法中水平向冲击下的等效弹塑性模型的示意图。
图5为本发明的冲击性能计算方法中恢复系数和触地瞬时速度与落石初速度之比的关系图。
图6为本发明的冲击性能计算方法与既有计算方法的计算结果对比图。
图7为本发明的冲击性能计算方法与现有试验方法的计算结果对比图。
上述附图中的有关标记为：
100-落石，200-管道。
具体实施方式
下面结合附图对本发明进行清楚、完整的说明。本领域普通技术人员在基于这些说明的情况下将能够实现本发明。在结合附图对本发明进行说明前，需要特别指出的是：
本发明中在包括下述说明在内的各部分中所提供的技术方案和技术特征，在不冲突的情况下，这些技术方案和技术特征可以相互组合。
此外，下述说明中涉及到的本发明的实施例通常仅是本发明一部分的实施例，而不是全部的实施例。因此，基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
关于本发明中术语和单位。本发明的说明书和权利要求书及有关的部分中的术语“包括”、“具有”以及它们的任何变形，意图在于覆盖不排他的包含。
本发明的落石冲击埋地管道的冲击性能计算方法的具体实施方式为包括以下步骤：
获取落石参数和地质参数，所述落石参数包括落石球体等效半径、落石质量、落石冲击埋地管道的对地面入射角度与入射速度，所述地质参数包括管道的覆土层厚度、泊松比、土层剪切波速和土层密度；
将所述落石参数和地质参数输入至落石冲击计算模型中，获取所述落石冲击计算模型输出的冲击性能数据；如图1所示，将落石100对埋地管道200的斜向冲击分解为竖向冲击与水平向冲击，因此，所述落石冲击计算模型又包括竖向冲击计算模型和水平向冲击计算模型。
所述的冲击性能数据包括落石的最大侵入深度、最大冲击力和冲击动力放大系数。
以下对本发明的落石冲击计算模型进行详细说明。
1.竖向冲击计算模型
落石100在竖向以一定的速度冲击土层，然后在土层表面做自由振动，因此，将土层简化为如图2所示的阻尼-弹簧等效模型。
如图2所示，建立竖向冲击的动力平衡方程为式1。
式1：
式中，左边两项为落石100与地面瞬间接触时所产生的作用力，右边两项为土层对落石100的阻力与变形阻力，整理后求得通解为式2.
式2：
根据位移边界条件为u(t)|t＝0＝0，速度边界条件为以此求解得竖向冲击地面时的位移时程方程即为式3，对式3求二阶导数得到竖向冲击地面时的落石加速度时程方程即为式4。
式3：
式4：
为分析落石100在地层中侵入特征，需要求得落石最大侵入深度，对此可令位移时程方程的一阶导数为零，即可得到竖向冲击时的落石最大侵入深度所对应的时间的计算式为式5：
式5：
从而，将tA代入式3即可求得落石最大侵入深度。
2.水平向冲击计算模型
水平向冲击计算模型可以根据土层性质进一步分为等效弹簧模型和等效弹塑性模型。当土层工程性质较好时，即水平向恢复系数较大(优选为大于0.4)时，可采用等效弹簧模型，反之，可采用等效弹塑性模型。
2.1等效弹簧模型
图3为本发明的冲击性能计算方法中水平向冲击下的等效弹簧模型的示意图。
如图3所示，在等效弹簧模型中，落石100在地层中沿水平向运动过程中受到一个等效弹簧的变形阻力，动力平衡方程为式6。
式6：
根据位移边界条件为v(t)|t＝0＝0，速度边界条件为求得式6所对应的位移时程方程为式7，再次对式7求二阶导数即可得到在水平向的落石加速度时程方程为式8。
式7：
式8：
2.2等效弹塑性模型
图4为本发明的冲击性能计算方法中水平向冲击下的等效弹塑性模型的示意图。
如图4所示，在等效弹塑性模型中，落石水平向运动过程中可认为受到一个等效阻尼和等效弹簧作用力，动力平衡方程为：
式9：
根据位移边界条件为v(t)|t＝0＝0，速度边界条件为求得式9所对应的位移时方程为式10，再次对式10求二阶导数即可得到水平向的落石加速度时程方程为式11。
式10：
式11：
3.最大冲击力和冲击动力放大系数
将上述所求得竖向加速度时程方程与水平向加速度时程方程采用矢量叠加即可求得斜向冲击时的落石加速度a计算式为式12：
式12：
令式12的一阶导数为零，可得到斜向冲击时的最大加速度对应的时间tmax及最大加速度amax，分别为式13和式14。
式13：
式14：
由此，最大冲击力Fmax可以采用以下两种计算方式。
3.1第一种计算方式
首先，将tmax分别代入式4、8或11，求得竖向与水平向的最大加速度，然后代入根据牛顿第二定律F＝m1a(F为或)，求出竖向冲击的最大冲击力为式15，水平向等效弹簧模型的最大冲击力为式16，水平向等效弹塑性模型的最大冲击力为式17。
式15：
式16：
式17：
然后，根据式18即可计算得到最大冲击力。
式18：
3.2第二种计算方式
首先根据式14求解出落石最大加速度amax，然后代入牛顿第二定律Fmax＝m1amax，即可求出Fmax。
3.3冲击动力放大系数的计算
所述冲击动力放大系数为最大冲击力与落石重力的比值，按照式19计算。
式19：
可以按照式20和式21分别求得竖向放大系数Kv和水平向放大系数Kh。
式20：
式21：
在式1～式21中：K为落石冲击动力放大系数；Fmax为最大冲击力；m1为落石质量；g为重力加速度；为最大冲击力的竖向分力，为最大冲击力的水平向分力；e为自然指数，取值为2.71828；λ为土层的竖向特征参数；λh＝为土层的水平向特征参数；t为冲击作用时间；tA为最大侵入深度相应时间；tmax为斜向冲击最大加速度所对应的时间；为落石100与土层接触后在竖向的瞬时速度，其与竖向的初速度关系可由图5所示的恢复系数Rc和触地瞬时速度v2与落石初速度v1之比v2/v1的关系图确定；为落石100与土层接触后在水平向的瞬时速度，其与在水平向的接触土层前的初速度关系可由图5所示的恢复系数Rc和触地瞬时速度v2与落石初速度v1之比v2/v1的关系图确定；ωd为竖向有阻尼圆频率；ωn为竖向无阻尼圆频率；为水平向有阻尼圆频率；u(t)为竖向冲击的位移；ü(t)为竖向的加速度，是t的函数，令t＝tmax求解得到ü(tmax)；v(t)为水平向冲击的位移；为水平向的加速度，是t的函数，令t＝tmax求解得到a为斜向冲击的落石加速度；amax为斜向冲击的落石最大加速度；kz为土层竖向弹簧系数；kh为土层水平向弹簧系数，其值为kz的0.6～0.8倍；cz为土层的竖向阻尼系数；ch为土层的水平向阻尼系数，其值为cz的0.6～0.8倍；A、B为方程通解中的待定系数，根据具体问题的初始与边界条件确定。
在上述符号中：
其中，r0为落石100的球体等效半径；G、ρ、μ分别为冲击土层的剪切模量、密度和泊松比；vs为土层剪切波速；β为土层竖向阻尼比；为土层水平向无阻尼圆频率；βh为土层水平向阻尼比。
以下通过具体的应用例来说明本发明的有益效果。
应用例
地面下土层中有一埋地管道200，落石仅沿竖向冲击地面，从而对埋地管道200产生影响作用。落石的初速度取为其它相关参数如表1所示。
(一)采用本发明的计算方法和既有方法包括路基规范法[1]、隧道手册方法[2](管道埋深为1.5m、土层弹性模量为35MPa)、杨其新-关宝树法[3]和日本道路公团法[4]对最大冲击力进行计算，计算结果如图6所示。
[1]交通部第二公路勘察设计院.公路路基设计规范(JTJ13-95)[S].北京:人民交通出版社，1995。
[2]铁道部专业设计院，铁路工程设计技术手册:隧道[M].北京:中国铁道出版社，1995。
[3]杨其新，关宝树.落石冲击力计算方法的试验研究[J].铁道学报，1996,18(1):101-106。
[4]Compilation by Japan RoadAssociation.Handbook for rockfallmeasures.Tokyo,2000.(in Japanese)。
表1
由图6可见，随落石质量的增加，最大冲击力呈现不断增大的趋势，但是不同的计算方法结果差异较大，其中，本发明的计算方法的计算结果最大，日本道路公团法[4]其次，隧道手册法[2]最小，路基规范法[1]与杨其新-关宝树法[3]居中。
分析以上结果差异较大的主要原因为：路基规范法[1]、隧道手册法[2]及杨其新-关宝树法[3]计算的落石冲击力实质上均为落石平均冲击力，而非落石的最大冲击力，这也是实际中导致防护结构失效的原因之一。日本道路公团法[4]是基于落石现场试验实测冲击力结果所拟合的经验公式，所反映的冲击力即为最大冲击力，它的计算结果远大于理论方法计算所得的平均冲击力，但相较于以往的试验得到的冲击力[4]，日本道路公团法[4]仍偏小30％～60％。与日本道路公团法[4]相比较，本发明的计算方法计算结果约大37％。因此，就以往试验结果，综合比较而言，本发明的计算方法更加符合实际。
(二)在第一部分的基础上，计算质量为2000kg的落石100的最大冲击力随落距的变化值，并将其与Pichler等[7]所做的现场试验值进行对比，结果如图7所示。
[7]Pichler B.,Hellmich Ch.,Mang H.A.Impact ofrocks onto gravel Designand evaluation of experiments[J].International Journal ofImpact Engineering,2005,31(5):559-578.
从图7可以看出，本发明的计算方法的计算结果介于试验值的范围内，符合实际。
(三)在第一部分的基础上，对以下几种具体情况分别进行计算，并给出相应的冲击动力放大系数计算结果表。
所述典型情况包括：根据土层的剪切波速Vs不同，将土层分为三个等级，即：Vs＝150m/s、Vs＝300m/s和Vs＝450m/s的工程性质较差土、中等土和良好土。竖向采用等效弹塑性模型，考虑到水平向恢复系数较大，故水平向采用等效弹簧模型。
在每种情况中，将落石100对地面的入射角度、入射速度(即为落石初速度v1)、落石球体等效半径r0作为可变参数，以便工程应用，下面直接给出竖向与水平向对应的冲击动力放大系数。
(1)工程性质较差土(Vs＝150m/s)
根据土层性质，取竖向水平向落石密度ρ1取为2500kg/m3，土层密度ρ取为1600kg/m3，泊松比μ取为0.35，入射角θ为90°时冲击动力放大系数K大小等于竖向放大系数Kv，计算结果如表2～7所示。
表2为θ＝90°，时的冲击动力放大系数计算表
落石的球体等效半径 落石质量 最大冲击力时间 竖向放大系数
r0(m) m1(kg) tmax(s) Kv
0.1 10.47 0.00145 377.6
0.2 83.73 0.00290 188.6
0.3 282.60 0.00437 125.5
0.4 669.87 0.00584 94.1
0.5 1308.33 0.00732 75.1
表3为θ＝90°，时的冲击动力放大系数计算表
落石的球体等效半径 落石质量 最大冲击力时间 竖向放大系数
r0(m) m1(kg) tmax(s) Kv
0.1 10.47 0.00144 755.8
0.2 83.73 0.00289 377.6
0.3 282.60 0.00434 251.6
0.4 669.87 0.00580 188.6
0.5 1308.33 0.00726 150.8
表4为θ＝90°，时的冲击动力放大系数计算表
落石的球体等效半径 落石质量 最大冲击力时间 竖向放大系数
r0(m) m1(kg) tmax(s) Kv
0.1 10.47 0.00144 1133.9
0.2 83.73 0.00289 566.7
0.3 282.60 0.00434 377.6
0.4 669.87 0.00579 283.2
0.5 1308.33 0.00725 226.5
表5为θ＝60°，v1＝10m/s时的冲击动力放大系数计算表
表6为θ＝60°，v1＝20m/s时的冲击动力放大系数计算表
表7为θ＝60°，v1＝30m/s时的冲击动力放大系数计算表
(2)工程性质中等土(Vs＝300m/s)
根据土层性质，取竖向水平向落石密度ρ1取为2500kg/m3，土层密度ρ取为1800kg/m3，泊松比μ取为0.30，入射角为90°时冲击动力放大系数K大小等于竖向放大系数Kv，计算结果如表8～13所示。
表8为θ＝90°，时的冲击动力放大系数计算表
表9为θ＝90°，时的冲击动力放大系数计算表
表10为θ＝90°，时的冲击动力放大系数计算表
表11为θ＝60°，v1＝10m/s时的冲击动力放大系数计算表
表12为θ＝60°，v1＝20m/s时的冲击动力放大系数计算表
表13为θ＝60°，v1＝30m/s时的冲击动力放大系数计算表
(3)工程性质良好土(Vs＝450m/s)
根据土层性质，取竖向水平向落石密度ρ1取为2500kg/m3，土层密度ρ取为2000kg/m3，泊松比μ取为0.25，入射角为90°时冲击动力放大系数K大小等于竖向放大系数Kv，计算结果如表14～19所示。
表14为θ＝90°，时的冲击动力放大系数计算表
表15为θ＝90°，时的冲击动力放大系数计算表
表16为θ＝90°，时的冲击动力放大系数计算表
表17为θ＝60°，v1＝10m/s时的冲击动力放大系数计算表
表18为θ＝60°，v1＝20m/s时的冲击动力放大系数计算表
表19为θ＝60°，v1＝30m/s时的冲击动力放大系数计算表
由此，在实际工程中，既可以根据落石对地面的入射速度、落石球体等效半径等参数查询表2～表19而直接获取落石的冲击性能数据，也可以将参数输入到落石冲击计算模型中，计算获取所述落石冲击计算模型输出的冲击性能数据。
综上可知，本发明具有以下优点：
(1)本发明考虑了落石最为常见的斜向冲击问题，并将斜向冲击问题转化为竖向冲击与水平向冲击，与公路路基规范、铁路隧道手册及既有文献的相关算法有根本区别，解决了非垂直冲击作用时冲击力的计算问题。
(2)本发明以碰撞理论为基础，并基于竖向与水平向不同的恢复系数，对地层土体做出合理简化，考虑了土层的塑性变形，更符合土层变形特性。
(3)不同于公路路基规范、铁路隧道手册及既有文献的相关算法，本发明方法所求得的是落石的最大冲击力，而非既有方法的平均冲击力。同时，本发明方法明确给出竖向与水平向的位移时程方程与加速度时程方程，可以直接计算出落石冲击过程中的最大位移及最大加速度，从而求出落石实际的最大冲击力以及最大侵入深度。
(4)本发明中的算法概念明确，计算简单，计算参数获取方便，实际工程中应用便捷。
本发明的落石冲击埋地管道的冲击性能计算设备的实施例为包括处理器、通信接口、存储器和通信总线，处理器、通信接口和存储器通过通信总线完成相互间的通信；所述存储器用于存放计算机程序；所述处理器用于执行存储器上所存放的程序时，实现上述的落石冲击埋地管道的冲击性能计算方法。
本发明的计算机可读存储介质的实施例为存储有一个或者多个程序，一个或者多个程序可被一个或者多个处理器执行，以实现上述的落石冲击埋地管道的冲击性能计算方法。
在上述实施例中，可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件实现时，可以全部或部分地以计算机程序产品的形式实现。计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行计算机程序指令时，全部或部分地产生按照本发明实施例的流程或功能。计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。计算机指令可以存储在计算机可读存储介质中，或者从一个计算机可读存储介质向另一个计算机可读存储介质传输，例如，计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，DVD)、或者半导体介质(例如固态硬盘Solid State Disk(SSD))等。
以上对本发明的有关内容进行了说明。本领域普通技术人员在基于这些说明的情况下将能够实现本发明。基于本发明的上述内容，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都应当属于本发明保护的范围。

Claims (6)

1.落石冲击埋地管道的冲击性能计算方法，包括以下步骤：
获取落石参数和地质参数，所述落石参数包括落石的球体等效半径、落石质量、落石冲击埋地管道的对地面入射角度与入射速度，所述地质参数包括管道的覆土层厚度、泊松比、土层剪切波速和土层密度；
将所述落石参数和地质参数输入至落石冲击计算模型中，获取所述落石冲击计算模型输出的冲击性能数据；
所述的冲击性能数据包括落石的最大冲击力和冲击动力放大系数，所述冲击动力放大系数为落石最大冲击力与落石重力的比值；
其中，将落石对埋地管道的斜向冲击分解为竖向冲击与水平向冲击；
所述落石冲击计算模型包括：
其中，K为落石冲击动力放大系数；Fmax为最大冲击力；m1为落石质量；g为重力加速度；为最大冲击力的竖向分力；为最大冲击力的水平向分力；amax为斜向冲击的落石最大加速度；
的计算模型包括：
的计算模型包括：
当管顶覆土层视为等效弹簧模型时：
当管顶覆土层视为等效弹塑性模型时：
其中，e为自然指数；λ为土层的竖向特征参数；tmax为斜向冲击的落石最大加速度所对应的时间；为落石与土层接触后在竖向的瞬时速度；ωd为竖向有阻尼圆频率；ωn为竖向无阻尼圆频率；为落石与土层接触后在水平向的瞬时速度；kh为土层水平向弹簧系数；为水平向有阻尼圆频率；λh＝为土层的水平向特征参数。
2.如权利要求1所述的计算方法，其特征在于：
落石最大加速度及其所对应的时间的计算模型包括：
其中，t为冲击作用时间；tmax为落石加速度a的一阶导数为零时的作用时间，a为斜向冲击时的落石加速度；为竖向的加速度，是t的函数，令t＝tmax求解得到 为水平向的加速度，是t的函数，令t＝tmax求解得到
3.如权利要求2所述的计算方法，其特征在于：
竖向的加速度的计算模型包括：
水平向的加速度的计算模型包括：
当管顶覆土层视为等效弹簧模型时：
当管顶覆土层视为等效弹塑性模型时：
4.如权利要求1所述的计算方法，其特征在于：所述的冲击性能数据还包括落石的最大侵入深度；所述落石冲击计算模型包括：
其中，tA为最大侵入深度相应时间；u(t)为落石竖向冲击的位移，是t的函数，令t＝tA求解得到的u(tA)即为最大侵入深度；e为自然指数；λ为土层的竖向特征参数；t为冲击作用时间；为竖向的落石与土层接触后的瞬时速度；kz为土层的竖向弹簧系数；m1为落石质量；g为重力加速度；ωd为竖向有阻尼圆频率。
5.落石冲击埋地管道的冲击性能计算设备，其特征在于，包括处理器、通信接口、存储器和通信总线，处理器、通信接口和存储器通过通信总线完成相互间的通信；
存储器，用于存放计算机程序；
处理器，用于执行存储器上所存放的程序时，实现权利要求1-4中任一所述的落石冲击埋地管道的冲击性能计算方法。
6.计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现权利要求1-4中任一所述的落石冲击埋地管道的冲击性能计算方法。
CN202210990958.4A 2022-08-18 2022-08-18 落石冲击埋地管道的冲击性能计算方法、设备和存储介质 Active CN115575004B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210990958.4A CN115575004B (zh) 2022-08-18 2022-08-18 落石冲击埋地管道的冲击性能计算方法、设备和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210990958.4A CN115575004B (zh) 2022-08-18 2022-08-18 落石冲击埋地管道的冲击性能计算方法、设备和存储介质

Publications (2)

Publication Number Publication Date
CN115575004A CN115575004A (zh) 2023-01-06
CN115575004B true CN115575004B (zh) 2024-07-09

Family

ID=84579302

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210990958.4A Active CN115575004B (zh) 2022-08-18 2022-08-18 落石冲击埋地管道的冲击性能计算方法、设备和存储介质

Country Status (1)

Country Link
CN (1) CN115575004B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103308663A (zh) * 2013-05-13 2013-09-18 西南交通大学 岩土工程滑坡模型试验箱

Family Cites Families (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101242491B1 (ko) * 2012-07-11 2013-03-12 한국지질자원연구원 토석류 내부충격력 측정장치
US9904579B2 (en) * 2013-03-15 2018-02-27 Advanced Elemental Technologies, Inc. Methods and systems for purposeful computing
CN104794352B (zh) * 2015-04-24 2017-11-10 成都理工大学 一种滚石的冲击力计算方法
DE102015119461B4 (de) * 2015-11-11 2022-11-17 Fujitsu Client Computing Limited Anordnung zum Festlegen einer Einschubbaugruppe in einem Einbaukäfig eines Computersystems sowie Verfahren zum Herstellen einer Anordnung
CN105590030A (zh) * 2016-01-28 2016-05-18 成都理工大学 水下滚石冲击力的计算方法
CN105740620B (zh) * 2016-01-28 2019-01-29 成都理工大学 泥石流冲击力的计算方法
CN106529198B (zh) * 2016-12-13 2019-01-29 中国科学院、水利部成都山地灾害与环境研究所 一种泥石流全过程数值模拟及数值计算方法
CN107239645B (zh) * 2017-08-10 2020-03-10 西南交通大学 落石冲击下无回填土拱形明洞结构概率可靠度设计方法
CN108733942B (zh) * 2018-05-28 2020-03-10 中国地质大学（武汉） 隧道围岩可靠性的评估方法及装置
CN108846172B (zh) * 2018-05-28 2019-08-30 中国地质大学（武汉） 泥石流冲击荷载函数生成方法及装置
CN108804785B (zh) * 2018-05-28 2019-07-09 中国地质大学（武汉） 一种冲击荷载计算方法及装置
CN109214064B (zh) * 2018-08-13 2022-09-23 长安大学 简支梁结构的水平层状围岩隧道临界开挖跨度计算方法
JP3218943U (ja) * 2018-09-07 2018-11-15 四川興邦工程勘察有限公司 落石衝撃力試験装置
CN110222369B (zh) * 2019-05-05 2022-11-22 西南交通大学 一种考虑回填缓冲层材料强化的落石冲击力计算方法
CN110334455B (zh) * 2019-07-11 2020-12-18 四川大学 一种水平式泥石流拦砂坝长度确定计算方法
CN110702302A (zh) * 2019-10-16 2020-01-17 西南交通大学 一种泥石流撞击力测量装置
CN111737871B (zh) * 2020-06-23 2024-02-20 重庆交通大学 一种结合岩土材料特性的落石轨迹三维预测分析方法
CN111852567A (zh) * 2020-07-10 2020-10-30 西南交通大学 一种爆破冲击作用下支护钢拱架动力响应评估方法
CN114169053B (zh) * 2021-12-03 2024-04-30 中铁第一勘察设计院集团有限公司 落石最大冲击力计算方法
CN114611195B (zh) * 2022-03-18 2024-06-04 西南交通大学 基于落石质量与下落高度的缓冲层-拱形明洞设计方法
CN114818294B (zh) * 2022-04-15 2023-05-02 成都理工大学 滚石冲击测量方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103308663A (zh) * 2013-05-13 2013-09-18 西南交通大学 岩土工程滑坡模型试验箱

Also Published As

Publication number Publication date
CN115575004A (zh) 2023-01-06

Similar Documents

Publication Publication Date Title
Huang et al. Seismic vulnerability of circular tunnels in soft soil deposits: The case of Shanghai metropolitan system
Tekie et al. Seismic fragility assessment of concrete gravity dams
Argyroudis et al. Analytical seismic fragility functions for highway and railway embankments and cuts
US11493656B1 (en) Full probability-based seismic risk analysis method for tunnel under fault dislocation
Haiyang et al. Seismic performance of underground subway station structure considering connection modes and diaphragm wall
Nie et al. Deformation characteristics and mechanism of the landslide in West Open-Pit Mine, Fushun, China
Shen et al. Seismic response of shield tunnel structure embedded in soil deposit with liquefiable interlayer
Suh et al. Subsidence hazard assessment at the Samcheok coalfield, South Korea: a case study using GIS
Dadfar et al. Vulnerability of buried energy pipelines subject to earthquake-triggered transverse landslides in permafrost thawing slopes
Xu et al. Dynamic response of buried gas pipeline under excavator loading: Experimental/numerical study
CN115575004B (zh) 落石冲击埋地管道的冲击性能计算方法、设备和存储介质
Li et al. Experimental research on seismic response of split‐type prefabricated utility tunnels through shaking table tests
Mantakas et al. Seismic response of a motorway bridge founded in an active landslide: a case study
Fan et al. Stochastic dynamic response analysis and probability evaluation of subway station considering subjected to stochastic earthquake excitation
Chen et al. 3D seismic response characteristics of a pile-mat-founded AP1000 nuclear-island building considering nonlinear hysteretic behavior of soil
Tavakoli Mehrjardi et al. Numerical modeling of buried steel pipe subjected to impact load
Sun et al. Numerical study on the use of soft material walls to enhance seismic performance of an existing tunnel
Xu et al. Buffering effect of overlying sand layer technology for dealing with rockfall disaster in tunnels and a case study
Xia et al. Stochastic method for predicting risk of slope failure subjected to unsaturated infiltration flow
Shi et al. A Dynamic Strain‐Rate‐Dependent Contact Model and Its Application in Hongshiyan Landslide
Park et al. Assessment of risk potential due to underground box structure installation employing ANN model and field experimental approaches
Cui et al. Structure strengthening method for enhancing seismic behavior of soft tunnel portal section
Garcia et al. Space-Time Pressure Distribution Applied to a Stiff Concrete Structure through a Protective Sand Layer from Full-Scale Experimental Rockfall Tests
Bong et al. Geostatistical properties of liquefaction vulnerability index for reclaimed land in South Korea
Nandi et al. Seismic stability assessment of rock slopes using limiting slope face concept

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