CN114169053A - 落石最大冲击力计算方法 - Google Patents

落石最大冲击力计算方法 Download PDF

Info

Publication number
CN114169053A
CN114169053A CN202111471102.8A CN202111471102A CN114169053A CN 114169053 A CN114169053 A CN 114169053A CN 202111471102 A CN202111471102 A CN 202111471102A CN 114169053 A CN114169053 A CN 114169053A
Authority
CN
China
Prior art keywords
impact force
poisson
ratio
elastic modulus
sphere
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.)
Granted
Application number
CN202111471102.8A
Other languages
English (en)
Other versions
CN114169053B (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.)
China Railway First Survey and Design Institute Group Ltd
Original Assignee
China Railway First Survey and Design Institute Group Ltd
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 China Railway First Survey and Design Institute Group Ltd filed Critical China Railway First Survey and Design Institute Group Ltd
Priority to CN202111471102.8A priority Critical patent/CN114169053B/zh
Publication of CN114169053A publication Critical patent/CN114169053A/zh
Application granted granted Critical
Publication of CN114169053B publication Critical patent/CN114169053B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Devices Affording Protection Of Roads Or Walls For Sound Insulation (AREA)

Abstract

本发明涉及一种落石最大冲击力计算方法。现有落石冲击力计算均计算的是落石平均冲击力,而不是最大冲击力,因而导致在工程应用中计算结果偏小。本发明将落石视为球体,获取球体弹性模量、拦石墙弹性模量、球体泊松比、拦石墙泊松比,计算与冲击体各自泊松比、弹性模量相关的参数;基于所述与冲击体各自泊松比、弹性模量相关的参数,计算落石最大冲击力。本发明着重计算危岩落石的最大冲击力,不同于现有技术的计算获得落石平均冲击力,最大冲击力高于平均冲击力,如此计算出最高可能值,对工程前期设计更具有安全意义。

Description

落石最大冲击力计算方法
技术领域
本发明涉及落石拦挡结构设计技术领域,具体涉及一种落石最大冲击力计算方法。
背景技术
进入21世纪以来,随着国内高速公路和高速铁路的大力修建,落石发生在国内公路、铁路的频率增加,会对公路、铁路等基础设施、设备和车辆及人员造成重大损失,落石灾害逐渐被国内学者提出来,落石理论的研究受到越来越多学者的重视。
我国《铁路隧道设计规范》明洞和棚洞设计荷载相关条文中对落石冲击力没有给出相关计算方法,仅规定通过现场调查或有关计算验证确定。《公路路基设计规范》、《铁路路基设计规范》和铁路行业《路基手册》等均没有给出落石冲击力计算方法。国内落石冲击力计算方法可参照的仅有《公路隧道设计细则》明洞设计章节和铁路行业《隧道手册》,以及杨其新、关宝树教授建立的基于室内试验所得的冲击力计算方法。细则方法、隧道手册方法、杨其新和关宝树提出的算法,计算的是落石平均冲击力,而不是最大冲击力,因而导致在工程应用中计算结果偏小。国内其他学者研究成果中冲击力计算都考虑缓冲或为平均冲击力。
发明内容
本发明的目的是提供一种落石最大冲击力计算方法,以至少解决现有技术仅计算落石平均冲击力、而不是最大冲击力导致的计算结果偏小的问题。
为了达到上述目的,本发明所采用的技术方案为:
落石最大冲击力计算方法,其特征在于:
所述方法包括以下步骤:
将落石视为球体,获取球体弹性模量、拦石墙弹性模量、球体泊松比、拦石墙泊松比,计算与冲击体各自泊松比、弹性模量相关的参数;
Figure BDA0003392263250000021
式中:
E——与冲击体各自泊松比、弹性模量相关的参数,Pa;
E1——球体弹性模量,Pa;
E2——拦石墙弹性模量,Pa;
μ1——球体泊松比;
μ2——拦石墙泊松比;
基于所述与冲击体各自泊松比、弹性模量相关的参数,计算落石最大冲击力;
Figure BDA0003392263250000022
式中:
Pmax——落石最大冲击力,N;
R1——球体半径,m;
m1——球体的质量,kg;
v1——球体的初始速度,m/s。
进一步地,所述球体弹性模量和球体泊松比根据岩石力学性质指标的经验数据取值,或通过岩石物理力学性质试验获得。
进一步地,所述拦石墙弹性模量和拦石墙泊松比根据混凝土结构设计规范取值。
进一步地,所述球体的初始速度通过RocFall软件计算获得。
与现有技术相比,本发明的有益效果如下:
1、本发明考虑到落石对刚性拦石墙冲击状态的受力复杂性,基于拦石墙主要受落石冲击力,在落石最大冲击力计算中将落石简化为球体,将拦石墙简化为平面体,计算过程更明确简单,计算结果准确可靠。
2、长期以来,落石最大冲击力的计算一直缺乏适用性的理论计算公式作为设计指导,本发明方案的提出对工程前期设计的落石拦挡结构受力情况分析意义重大。
3、本发明着重计算落石的最大冲击力,不同于现有技术的计算获得落石平均冲击力,最大冲击力高于平均冲击力,如此计算出最高可能值,对工程前期设计更具有安全意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他实施例的附图。
图1是实施例中把边坡坡型导入到RocFall软件中的示意图。
图2是实施例中输入边坡法向、切向碰撞恢复系数和粗糙度标准偏差的示意图。
图3是实施例中输入落石的原始速度和质量的示意图。
图4为实施例中计算落石在拦石墙位置运动速度的示意图。
具体实施方式
为了便于理解本发明,下面将参照相关附图对本发明进行更全面的描述。附图中给出了本发明的较佳实施方式。但是,本发明可以以许多不同的形式来实现,并不限于本文所描述的实施方式。相反地,提供这些实施方式的目的是使对本发明的公开内容理解的更加透彻全面。
应注意到,相似的标号和字母表示类似项,因此,一旦某一项在一个实施例中被定义,则在随后的实施例中不需要对其进行进一步定义和解释。
本发明提供了一种落石最大冲击力计算方法,具体来说,该方法将落石视为球体,将拦石墙视为被撞击面,根据球体半径、球体质量、球体初始速度、球体泊松比、球体弹性模量、以及拦石墙泊松比、拦石墙弹性模量等关键因素,计算出落石的最大冲击力。所述方法具体包括以下计算过程:
将落石视为岩石球体,将拦石墙视为混凝土被撞击面,获取球体弹性模量、拦石墙弹性模量、球体泊松比、拦石墙泊松比,计算与冲击体各自泊松比、弹性模量相关的参数;
Figure BDA0003392263250000041
式中:
E——与冲击体各自泊松比、弹性模量相关的参数,Pa;
E1——球体弹性模量,Pa;
E2——拦石墙弹性模量,Pa;
μ1——球体泊松比;
μ2——拦石墙泊松比;
基于所述与冲击体各自泊松比、弹性模量相关的参数,计算落石最大冲击力;
Figure BDA0003392263250000042
式中:
Pmax——落石最大冲击力,N;
R1——球体半径,m;
m1——球体的质量,kg;
v1——球体的初始速度,m/s。
所述球体弹性模量和球体泊松比根据岩石力学性质指标(如《工程地质手册》的岩石力学性质指标)的经验数据取值,或通过岩石物理力学性质试验(如岩石物理力学性质试验规程(DZT 0276.19-2015)第19部分岩石单轴压缩变形试验)获得。
所述拦石墙弹性模量和拦石墙泊松比根据混凝土结构设计规范(如混凝土结构设计规范(GB 50010-2010)4.1.5条)取值。
所述球体的初始速度通过RocFall软件计算获得。
应注意到,上述方法中,所述球体弹性模量、球体泊松比、拦石墙弹性模量、拦石墙泊松比以及球体的初始速度等参数,均基于现有成熟的公开方法和应用软件获得,本发明对方法的改进并不在于这些参数的获取过程,任何基于现有技术或公知方法获得上述参数的手段均在本申请的保护范围内。
实施例:
(1)将落石简化为球体并假定简化为球体的落石的半径R1为0.5m,将拦石墙简化为混凝土平面;
(2)获取球体弹性模量E1和球体泊松比μ1,根据《工程地质手册》中岩石力学性质指标的经验数据取值,经验数据取值如表1:
表1岩石力学性质指标的经验数据取值
Figure BDA0003392263250000051
以花岗岩为例,密度取表1中花岗岩密度的平均值,为2680kg/m3,球体弹性模量取表1中的静模量平均值,即E1取3.5×104MPa,球体泊松比表1中的泊松比平均值,即μ1取0.26。
(3)获取拦石墙弹性模量E2和拦石墙泊松比μ2,根据混凝土结构设计规范(GB50010-2010)4.1.5条,选取拦石墙弹性模量E2如表2:
表2混凝土弹性模量(×104N/mm2)
Figure BDA0003392263250000061
以C35混凝土为例,拦石墙弹性模量E2取3.15×104MPa,拦石墙泊松比μ2即各等级混凝土泊松比按0.2计算(根据混凝土结构设计规范(GB 50010-2010)4.1.5条的规定和建议)。
(4)根据球体弹性模量E1、球体泊松比μ1、拦石墙弹性模量E2和拦石墙泊松比μ2,计算与冲击体各自泊松比、弹性模量相关的参数E:
Figure BDA0003392263250000062
求得E=1.75×1010Pa。
(5)计算球体的质量m1,假定简化为球体的落石的半径为0.5m,其密度为2680kg/m3,则求得球体的质量m1为1403.2kg。
(6)获取球体的初始速度v1,假定落石半径为0.5m,通过RocFall软件对落石的冲击速度进行计算,计算得到落石在拦石墙位置运动速度为4m/s,即v1为4m/s。
RocFall软件的具体操作过程为:
把边坡坡型导入到RocFall软件中(图1),输入边坡法向、切向碰撞恢复系数和粗糙度标准偏差(图2),以及落石的原始速度和质量(图3),让软件开始计算,可以得到落石的冲击速度曲线(图4)。根据冲击速度曲线可以得到落石在边坡不同位置的冲击速度,可为落石防护措施的设计提供参数。
其中法向和切向碰撞恢复系数可根据《铁路特殊路基设计规范》(TB10035-2018)条文说明表12.2.6(表3)对其进行选取。
表3法向和切向碰撞恢复系数
Figure BDA0003392263250000071
注:坡面被灌木丛覆盖,Rt减少0.1。
软件使用手册指出粗糙度标准偏差在0~5之间,一般不会超过5,软件计算可以选择输入粗糙度标准偏差的值为平均值2.5。
以上给出了RocFall软件的具体计算过程,实施例中的参数值和附图中涉及的输入参数值均是为了具体描述实施例的过程,不应理解为对参数的唯一限定。
(7)计算落石最大冲击力;
Figure BDA0003392263250000072
带入步骤(1)-(6)中的参数得到,Pmax=5.7×103kN。
上述实施例描述了将本发明的方法应用于落石撞击交通工程构筑物的过程,撞击体为动态球体,被撞击体为静态平面。本发明的方法还可应用于其他受力情况,如,撞击体和被撞击体均为动态,以及被撞击体也为球体的情况。还有,本发明的方法除运用到落石最大冲击力计算外,还可运用到滚石撞击最大冲击力计算。另外,本发明的方法还可运用到其他领域,如轮船靠岸碰撞、卫星回收落地碰撞等涉及碰撞的领域。基于本发明的方法,运用任何公知常识做改变的方案都在本发明的保护范围内。
以上应用了具体个例对本发明进行阐述,只是用于帮助理解本发明,并不用以限制本发明。对于本发明所属技术领域的技术人员,依据本发明的思想,还可以做出若干简单推演、变形或替换。

Claims (4)

1.落石最大冲击力计算方法,其特征在于:
所述方法包括以下步骤:
将落石视为球体,获取球体弹性模量、拦石墙弹性模量、球体泊松比、拦石墙泊松比,计算与冲击体各自泊松比、弹性模量相关的参数;
Figure FDA0003392263240000011
式中:
E——与冲击体各自泊松比、弹性模量相关的参数,Pa;
E1——球体弹性模量,Pa;
E2——拦石墙弹性模量,Pa;
μ1——球体泊松比;
μ2——拦石墙泊松比;
基于所述与冲击体各自泊松比、弹性模量相关的参数,计算落石最大冲击力;
Figure FDA0003392263240000012
式中:
Pmax——落石最大冲击力,N;
R1——球体半径,m;
m1——球体的质量,kg;
v1——球体的初始速度,m/s。
2.根据权利要求1所述的方法,其特征在于:
所述球体弹性模量和球体泊松比根据岩石力学性质指标的经验数据取值,或通过岩石物理力学性质试验获得。
3.根据权利要求1所述的方法,其特征在于:
所述拦石墙弹性模量和拦石墙泊松比根据混凝土结构设计规范取值。
4.根据权利要求1所述的方法,其特征在于:
所述球体的初始速度通过RocFall软件计算获得。
CN202111471102.8A 2021-12-03 2021-12-03 落石最大冲击力计算方法 Active CN114169053B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111471102.8A CN114169053B (zh) 2021-12-03 2021-12-03 落石最大冲击力计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111471102.8A CN114169053B (zh) 2021-12-03 2021-12-03 落石最大冲击力计算方法

Publications (2)

Publication Number Publication Date
CN114169053A true CN114169053A (zh) 2022-03-11
CN114169053B CN114169053B (zh) 2024-04-30

Family

ID=80482997

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111471102.8A Active CN114169053B (zh) 2021-12-03 2021-12-03 落石最大冲击力计算方法

Country Status (1)

Country Link
CN (1) CN114169053B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115575004A (zh) * 2022-08-18 2023-01-06 西南交通大学 落石冲击埋地管道的冲击性能计算方法、设备和存储介质
CN116955915A (zh) * 2023-07-21 2023-10-27 中南大学 一种测算落石碰撞恢复系数的方法及测算装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005163532A (ja) * 2003-11-12 2005-06-23 Nihon Samicon Co Ltd 落石防護構造物の設計方法と落石防護構造物
CN104794352A (zh) * 2015-04-24 2015-07-22 成都理工大学 一种滚石的冲击力计算方法
CN105590030A (zh) * 2016-01-28 2016-05-18 成都理工大学 水下滚石冲击力的计算方法
CN112818532A (zh) * 2021-01-27 2021-05-18 中国地质大学(北京) 崩塌落石地质灾害中钢筋混凝土桩板墙的动态工程响应测算方法
CN112924305A (zh) * 2021-01-27 2021-06-08 中国地质大学(北京) 获取被动桩板拦石墙在崩塌落石地质灾害中动态响应数据的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005163532A (ja) * 2003-11-12 2005-06-23 Nihon Samicon Co Ltd 落石防護構造物の設計方法と落石防護構造物
CN104794352A (zh) * 2015-04-24 2015-07-22 成都理工大学 一种滚石的冲击力计算方法
CN105590030A (zh) * 2016-01-28 2016-05-18 成都理工大学 水下滚石冲击力的计算方法
CN112818532A (zh) * 2021-01-27 2021-05-18 中国地质大学(北京) 崩塌落石地质灾害中钢筋混凝土桩板墙的动态工程响应测算方法
CN112924305A (zh) * 2021-01-27 2021-06-08 中国地质大学(北京) 获取被动桩板拦石墙在崩塌落石地质灾害中动态响应数据的方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115575004A (zh) * 2022-08-18 2023-01-06 西南交通大学 落石冲击埋地管道的冲击性能计算方法、设备和存储介质
CN116955915A (zh) * 2023-07-21 2023-10-27 中南大学 一种测算落石碰撞恢复系数的方法及测算装置
CN116955915B (zh) * 2023-07-21 2024-03-08 中南大学 一种测算落石碰撞恢复系数的方法及测算装置

Also Published As

Publication number Publication date
CN114169053B (zh) 2024-04-30

Similar Documents

Publication Publication Date Title
CN114169053A (zh) 落石最大冲击力计算方法
CN109657339B (zh) 一种铁道车辆坡道运行综合性能的评估方法
Taherinezhad et al. A review of behaviour of prestressed concrete sleepers
CN109670217B (zh) 隧道气动载荷影响因素的权重分析方法及计算机系统
CN104458173B (zh) 钢框架结构突变损伤的识别方法及系统
Recuero et al. Finite-element analysis of unsupported sleepers using three-dimensional wheel–rail contact formulation
CN111859580B (zh) 一种铁路线路线型动态分析与设计方法
Fu et al. Statistical analysis of the causes of bridge collapse in China
CN112818532A (zh) 崩塌落石地质灾害中钢筋混凝土桩板墙的动态工程响应测算方法
CN111325482A (zh) 一种大跨度地下洞库的围岩快速分级方法
CN107478803A (zh) 一种隧道悬臂掘进机施工适应性分级方法
CN110261578B (zh) 一种考虑结构面粗糙度的裂隙岩体稳定性分析系统
Benedetti et al. Evaluation of the seismic response of masonry buildings based on energy functions
CN117217562B (zh) 一种隧道涌水淹没-抽排水后的结构安全评价方法及系统
CN112924305B (zh) 获取被动桩板拦石墙在崩塌落石地质灾害中动态响应数据的方法
CN110688698B (zh) 一种基于随机森林算法的围岩大变形智能评估方法
CN110032757B (zh) 一种强夯施工振动对周边建筑物影响安全距离的计算方法
Zhou et al. Numerical simulation of a cable-stayed bridge subjected to ship collision
CN114970276A (zh) 一种基于bp神经网络的沥青路面潜在破坏参数计算方法
Xin et al. Seismic damage to non-fault-crossing and fault-crossing tunnels: Comparative study of the 2008 Wenchuan earthquake (Mw 7.9) and the 2022 Menyuan earthquake (Mw 6.7)
CN110298095B (zh) 一种门洞形裂隙岩体的破碎机理数值分析方法
CN112329206A (zh) 一种隧道施工非爆破区判定方法
CN113652980A (zh) 多级耗能缓冲系统、设计方法及拦石墙结构损伤监测方法
Vijayakumar et al. The High-Speed Railway Bridges Under Vehicle Moving Load And Near Fault Seismic Ground Motions–Review.
Kawecki et al. Consideration of criteria of vibration comfort of people in diagnosis and design of buildings

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