
Citation: | WANG Yusuo, LYU Ningning, WANG Mingnian, TANG Jianhui, YANG Junxiang, ZHAO Zhuang, XIAO Peng, AN Bo. Method of Determining Design Value of Rockfall Block Size Based on the Return Period[J]. Journal of Southwest Jiaotong University, 2025, 60(2): 262-270. doi: 10.3969/j.issn.0258-2724.20240240 |
The rockfall events on transport lines in mountainous areas occur frequently and are sudden and random. Probability-based reliability design of protective structures is a necessary and effective means to reduce the hazard risk. However, the method of determining the design value of rockfall size, as a basic parameter of the design, is still unclear, making the reliability design fail to be carried out effectively. The prediction method of rockfall size and the return period was introduced based on the Poisson-generalized Pareto distribution (GPD) composite model. According to the rockfall record data of the Bai−Upper section of Baoji−Chengdu Railway, the GPD threshold was selected, and the impact of the missing maximum and minimum rockfall sizes (vmax and vmin) in the record on the threshold selection and prediction results was analyzed. Corresponding to the exceedance probability of the standard value of earthquake action, a method of determining the design value of rockfall size was proposed, which was based on the three-level protection standard featuring “the structure doesn’t break under small rockfalls, can be repaired under medium rockfalls, and doesn’t collapse under large rockfalls”. The results show that the missing maximum rockfall size in the record has a small impact on the threshold selection but a large impact on the prediction results of the rockfall size. The prediction value of rockfall size corresponding to the return period increases exponentially with the set maximum rockfall size, while the recorded minimum rockfall size has less impact on both the threshold selection and the prediction results. By taking the record data of the Bai−Upper section of Baoji−Chengdu Railway as an example, according to the selected reasonable threshold, different maximum rockfall size values are set. By using the method presented in this paper, the interval range of rockfall size in return periods of 50, 100, 1 000, and 10 000 years on the Bai−Upper section of Baoji−Chengdu Railway is predicted, and the standard value (design value) of rockfall size corresponding to the three-level protection standard in design reference periods of 50 and 100 years is calculated, which can provide a reference for similar engineering designs.
波形钢腹板箱梁因为能减轻结构自重、减小预应力损失、避免腹板斜裂缝而成为具有广阔应用前景和吸引力的新型桥梁结构[1-2]. 不同于混凝土箱梁,波形钢腹板箱梁用波形钢腹板替换混凝土腹板后,箱梁的抗扭刚度被显著削弱[2-3],因此其扭转和畸变效应也更加突出.
国内外学者对波形钢腹板箱梁的约束扭转效应开展了不少研究工作. 杨丙文等[4]分析了单箱单室波形钢腹板箱梁的约束扭转效应,但在悬臂板自由端出现了剪应力,与实际不符. 马磊等[5]推导了单箱双室波形钢腹板箱梁扭转控制微分方程,但仅分析了翘曲正应力. 邓文琴等[6]推导了单箱三室波形钢腹板悬臂梁扭转控制微分方程,并给出了翘曲正应力和剪应力计算式. Nie等[7-8]等在波形钢腹板箱梁的试验研究结果中靠近顶底板局部区域内出现了可观的正应变. 张元海等[9]指出分析波形钢腹板箱梁的约束扭转时应考虑顶底板对波形钢腹板的约束作用. 以上都是基于乌曼斯基第二理论(下文简称乌-Ⅱ理论)建立约束扭转控制微分方程并开展的理论分析,而乌-Ⅱ理论应用了剪应力的两种算式,出现了两种不一致的纵向翘曲位移[10]. 张元海等[11]指出约束扭转剪应力应按微元体平衡条件或根据自由扭矩和二次扭矩计算,不能直接按胡克定律计算. 此外,Kato等[12]采用矩阵法提出了波形钢腹板箱梁扭转效应计算理论. 李运生等[13]采用能量变分法推导了波形钢腹板曲线梁弯扭控制微分方程,并采用伽辽金法求得其解析解. 徐勋等[14]在广义坐标法的基础上基于混合变分原理建立了薄壁箱梁约束扭转分析新理论,与Reissner原理的结果是统一的,均考虑了全部次生剪应力对中面剪切变形的影响,计算精度更高. 鲍永方等[15]对适用于刚性扭转的乌曼斯基理论和Reissner原理进行比较,指出Reissner原理的计算精度高于乌曼斯基理论.
为了更加合理地分析波形钢腹板箱梁约束扭转效应,本文将位移和应力同时作为变分参量,应用Reissner原理推导出约束扭转控制微分方程,给出了不同于乌-Ⅱ理论的翘曲系数. 同时考虑波形钢腹板的褶皱效应,合理计算截面几何特性,推演出翘曲正应力和剪应力的计算公式. 结合波形钢腹板简支箱梁数值算例,研究了应力分布特征及截面几何参数变化对应力的影响规律.
图1为波形钢腹板的形状简图. 图中:a、c、d分别为直板段长、斜板段投影长、斜板段长;hw、tw分别为波高和板厚.
因组成波形钢腹板箱梁的各板件材料不同,约束扭转分析时需通过等效原则将材料统一. 本文根据波形钢腹板的力学性能先将其等效为正交异性平钢板,再等效为混凝土腹板. 将波形钢腹板等效为正交异性平钢板后的等效剪切模量[16]为
Gse=(a+c)/(a+d)Gs=ηGs, | (1) |
式中:Gs为钢板剪切模量;
波形钢腹板在箱梁的受力中主要承受剪力,由式(1)采用剪切模量比可计算波形钢腹板等效为混凝土腹板的厚度,波形钢腹板的等效厚度为
twe=Gse/Gctw=ηnGtw, | (2) |
式中:Gc为混凝土剪切模量;
在一个波长的波形钢腹板两端施加一对等值、反向的轴向力,采用卡氏第二定理计算出弯矩和轴力下的轴向位移,根据位移等效原则可求得波形钢腹板等效为平钢板的纵向弹性模量为
Ese=κEs, | (3) |
式中:Es为钢板弹性模量;
κ=a+c(a+c2/d)(tw/hw)2+(3a+d)(twhw)2. |
根据
Ece=κnE/(ηnG)Ec = λEc, | (4) |
式中:Ec为混凝土弹性模量;
由于tw/hw一般很小,κ、λ则更小,波形钢腹板表现出明显的褶皱效应.
根据薄壁箱梁的约束扭转分析理论,结合波形钢腹板的力学性能,对波形钢腹板箱梁约束扭转分析时作以下基本假定:
1) 刚性周边假设,横截面周边不变形;
2) 横截面上的翘曲正应力和剪应力沿壁厚均匀分布;
3) 混凝土顶底板与波形钢腹板间不产生界面滑移,波形钢腹板不发生剪切屈曲;
4) 波形钢腹板箱梁处于线性弹性工作状态.
分析波形钢腹板箱梁的约束扭转效应时需考虑二次剪切变形的影响,引入广义翘曲位移
u(x,s)=−β′(x)ˉω(s), |
v(x,s)=φ(x)ρ(s), |
式中:
由式(5)可得壁厚中面任一点的翘曲正应变和剪应变分别为
ε(x,s)=∂u/∂x=−β″ˉω, |
γ(x,s)=∂u/∂s+∂v/∂x=φ′ρ−β′ˉω′. |
根据刚性周边假设得各板件横向正应变
σ={−Eoβ″ˉω, 顶底板,−λEoβ″ˉω,等效腹板, | (7) |
式中:
引起翘曲正应力的内力定义为扭转双力矩
Bˉω=∫s1σˉωtds=−EoIˉωβ″, | (8) |
式中:
联立式(7)和式(8)可得翘曲正应力与扭转双力矩的关系式为
σ={Bˉωˉω/Iˉω, 顶底板,λBˉωˉω/Iˉω,等效腹板. | (9) |
约束扭转剪应力应按微元体平衡条件或根据自由扭矩和二次扭矩计算. 约束扭转总扭矩Mz分解为自由扭转扭矩M1和二次扭矩
闭口箱壁中面上任一点的自由扭转剪应力为
τ1=Msc/(Ωt)=M1ψ/(Idt), | (10) |
式中:Msc为闭口箱壁承担的自由扭矩;
根据波形钢腹板箱梁的抗扭特性需对抗扭惯性矩进行修正,考虑悬臂板的贡献,抗扭惯性矩[17]为
Id=Ω2/Ω2[bttt(1−ζ)+bbtb(1−ζ)+2bwtwe(1+ζ)][bttt(1−ζ)+bbtb(1−ζ)+2bwtwe(1+ζ)]+2bft3t3, |
式中:ζ为修正系数,对矩形波形钢腹板箱梁,
将式(9)代入各板件微元体平衡方程
{\tau _{2}} = {q_{\bar \omega }}/t = - ({M_{2}}{\bar S_{\bar \omega }})/({I_{\bar \omega }}t), | (11) |
式中:
各板壁厚中面上任一点的
{\bar{S}}_{\bar{\omega }}=\left\{\begin{array}{l}{S}_{\bar{\omega }}, \quad 悬臂板,\\ {S}_{\bar{\omega }}-{\displaystyle \oint ({S}_{\bar{\omega }}/t\text{)d}s}\Big/{\displaystyle \oint \;1/t\text{d}s},\quad 闭口箱壁,\end{array}\right. |
式中:
由式(10)、(11)可得各板壁厚中面上任一点的总剪应力为
{\tau }_{{\text{z}}}=\left\{\begin{array}{l}-{M}_{2}{\bar{S}}_{\bar{\omega }}/({I}_{\bar{\omega }}t) , \quad \ 悬臂板 , \\ {M}_{1}\psi/({I}_{{\rm{d}}}t)-{M}_{2}{\bar{S}}_{\bar{\omega }}/({I}_{\bar{\omega }}t) , \quad 闭口箱壁. \end{array}\right. | (12) |
基于Reissner原理建立波形钢腹板箱梁的约束扭转控制微分方程时,将位移和应力同时作为了变分参量,建立的能量泛函[18]为
\begin{split} &{\varPi_{\rm{R}}}({\boldsymbol{\sigma}} ,{\boldsymbol{u}}) = \int_V {\left({{\boldsymbol{\sigma}} ^{\rm{T}}}{\boldsymbol{\varepsilon}} ({\boldsymbol{u}}) - {\boldsymbol{V}}({\boldsymbol{\sigma}} ) - {\boldsymbol{Fu}} \right)} {\text{d}}V -\\ &\quad \int_l {{\boldsymbol{Pu}}} {\text{d}}x, \end{split} | (13) |
式中:
能量泛函中的应变能为
\int_V {{{\boldsymbol{\sigma}} ^{\rm{T}}}{\boldsymbol{\varepsilon}} ({\boldsymbol{u}})} {\text{d}}V = \int_V {({\sigma }\varepsilon + {\tau _{\text{z}}}\gamma )} {\text{d}}V. | (14) |
将式(6)、(9)和式(12)代入式(14),得
\int_V {{{\boldsymbol{\sigma}} ^{\rm{T}}}{\boldsymbol{\varepsilon}} } ({\boldsymbol{u}}){\text{d}}V = \int_0^l {( - {B_{\bar \omega }}\beta '' + {M_{\text{z}}}\varphi ' - {M_{2}}\beta '{\text{)d}}x} . | (15) |
余能函数为
\int_V {{\boldsymbol{V}}({\boldsymbol{\sigma}} )} {\text{d}}V = \int_V {\left(\dfrac{{\sigma ^2}}{{2E_{\rm{o}}}} + \dfrac{{\tau _{\text{z}}^2}}{{2{G_{\rm{c}}}}}\right)} {\text{d}}V. | (16) |
将式(9)和式(12)代入式(16),得
\int_V {{\boldsymbol{V}}({\boldsymbol{\sigma}})} {\text{d}}V = \int_0^l \left({\dfrac{{{B_{\bar \omega }^2}}}{{2E_{\rm{o}}{I_{\bar \omega }}}} + \dfrac{{M_{\rm{1}}^2}}{{2{G_{\rm{c}}}{I_{\rm{d}}}}} + } \dfrac{{\chi M_{2 }^2}}{{2{G_{\rm{c}}}{I_{\rm{d}}}}}\right){\text{d}}x , | (17) |
式中:
当箱梁上仅作用分布扭矩mt时,外力势能为
- \int_0^l {{\boldsymbol{Pu{}}\text{d}}x} = - \int_0^l {{m_{\rm{t}}}\varphi {\text{d}}x} . | (18) |
将式(15)、(17)和式(18)代入式(13),并对能量泛函数
\begin{split} &\delta {\varPi_{\rm{R}}} = \int_0^l \left[{\left( - \beta '' - \dfrac{{{B_{\bar \omega }}}}{{E_{\rm{o}}{I_{\bar \omega }}}}\right)\delta {B_{\bar \omega }} + \left(\varphi ' - \dfrac{{{M_{\rm{s}}}}}{{{G_{\rm{c}}}{I_{\rm{d}}}}}\right)\delta {M_{1}}} + \right.\\ &\quad \left.\left(\varphi ' - \beta ' - \dfrac{{\chi {M_{2 }}}}{{{G_{\rm{c}}}{I_{\rm{d}}}}}\right)\delta {M_{2}}{\text{ + }}\left({B'_{\bar \omega }} - {M_{2}}\right)\delta \beta' - \right. \\ &\quad \Bigg.({M'_1}+ {M'_{2}} + {m_{\rm{t}}})\delta \varphi \Bigg]{\text{d}}x . \\[-12pt] \end{split} | (19) |
由
1) 力与位移的关系条件
\left\{\begin{array} {l} - \beta '' = \dfrac{{{B_{\bar \omega }}}}{{E_{\rm{o}}{I_{\bar \omega }}}}, \\ \varphi ' = \dfrac{{{M_{\rm{1}}}}}{{{G_{\rm{c}}}{I_{\rm{d}}}}}, \\ \varphi ' - \beta ' = \dfrac{{\chi {M_{2}}}}{{{G_{\rm{c}}}{I_{\rm{d}}}}}. \end{array} \right. | (20) |
2) 力的平衡条件
\left\{ \begin{array}{l} {B'_{\bar \omega }} = {M_{2}}, \\ {M'_1} + {M'_{2}} = - {m_{\rm{t}}} . \end{array} \right. | (21) |
借助式(20)和式(21),并忽略mt的高阶微分,可得波形钢腹板箱梁约束扭转控制微分方程为
\varphi^ {(4)} - k_{\text{R}}^2\varphi '' = {\xi _{\text{R}}}{m_{\rm{t}}}/(E_{\rm{o}}{I_{\bar \omega }}), | (22) |
式中:
不同于Reissner原理,采用乌Ⅱ-理论推导波形钢腹板箱梁约束扭转控制微分方程时,首先要确定广义翘曲位移
q_{\text{z}} = {\tau _{\text{z}}}t = {G_{\text{c}}}t(\varphi '\rho - \beta '\bar \omega '). | (23) |
将式(23)代入静力平衡方程
\varphi ' - {\xi _{\text{w}}}\beta = {M_{\text{z}}}/({G_{\rm{c}}}{I_\rho }). | (24) |
式中:
由微元体平衡方程
由功能原理及自由扭转剪力流qs满足的静力平衡方程
M1和M2叠加得Mz,即
{M_{\text{z}}} = {G_{\rm{c}}}{I_{\text{d}}}\varphi ' - E_{\rm{o}}{I_{\bar \omega }}\beta ''' . | (25) |
将式(24)及
\varphi^ {(4)} - k_{\text{w}}^2\varphi '' = {\xi _{\text{w}}}{m_{\rm{t}}}/(E_{\rm{o}}{I_{\bar \omega }}) , | (26) |
式中:
两种方法推导的约束扭转控制微分方程表达式相同,且均可采用初参数法求解[19],但翘曲系数公式不同. 采用乌-Ⅱ理论推导约束扭转控制微分方程时应用了由胡克定律和微元体平衡条件得到的剪力流的两种算式及两种不一致的纵向翘曲位移[10],且未考虑全部次生剪应力对中面剪切变形的影响. 应用Reissner原理建立扭转控制微分方程时将位移和应力同时作为变分参量,并考虑了全部次生剪应力对中面剪切变形的影响,通过建立能量泛函并对其变分,得出力与位移的关系条件及力的平衡条件,与刚性周边假设下的力与位移的关系条件及力的平衡条件是相符的,因此该理论更准确.
算例采用计算跨径为40 m的简支波形钢腹板箱梁,梁端设置1.5 m厚混凝土横隔板,跨中作用集中扭矩荷载
图3为按本文方法和乌-Ⅱ理论计算的扭转内力分布曲线. 由图3可以看出:按不同方法计算的双力矩和二次扭矩在跨中截面有明显差异,按乌-Ⅱ理论计算的双力矩和二次扭矩分别比本文方法计算的约增大了51%和128%. 由于两种方法计算的翘曲系数差异较大导致扭转内力的差异较大,也说明不同方法计算的翘曲系数对扭转内力有显著影响.
为验证本文方法合理性,按本文方法、乌-Ⅱ理论计算横截面翘曲正应力和总剪应力,并与有限元解对比. 有限元解是采用ANSYS软件建立空间有限元模型求得,有限元模型如图4所示. 混凝土板和波形钢腹板分别采用SOLID45实体单元和SHELL63壳单元模拟,由于分析时假设波形钢腹板与混凝土顶底板之间无界面滑移[20],故在模型中壳单元和实体单元间用主从约束法作刚性连接处理;每侧梁端底板上设置2个支座,左侧梁端约束为3个坐标轴方向的平动位移Ux 、Uy和Uz,右侧梁端约束为横向和竖向平动位移Uy和Uz;集中扭矩荷载等效成刚性扭转荷载的力流形式,均匀施加于跨中截面闭合箱壁各单元的节点上,顶底板施加反向水平荷载
为避免有限元模型中集中扭矩施加位置的局部应力影响,根据扭转内力分布及有限元模型的应力分布,选取距跨中1.6 m的左截面正面(外法线方向与x轴一致的横截面为正面)为分析截面. 以该截面上的Ⅰ、Ⅱ、Ⅴ 3个计算点进行翘曲正应力对比,Ⅲ、Ⅳ、Ⅵ 3个计算点进行总剪应力对比,各计算点的应力结果见表1. 表1中的相对误差为本文方法解或乌-Ⅱ理论解与有限元解的差除以有限元解求得. 由表1可以看出,本文方法计算的应力与有限元解吻合更好. 但是由于刚性周边等假设,理论解与有限元解之间还是存在一定的误差.
名称 | 计算点 | 有限元 解/kPa | 本文方 法/kPa | 乌-Ⅱ 理论/kPa | 相对误差/% | |
本文 方法 | 乌-Ⅱ 理论 | |||||
{\sigma} | Ⅰ | 54.06 | 57.35 | 63.61 | 6.09 | 17.67 |
Ⅱ | −27.83 | −25.42 | −28.18 | −8.66 | 1.28 | |
Ⅴ | −75.83 | −81.24 | −90.10 | 7.13 | 18.82 | |
{\tau _{\text{z}}} | Ⅲ | 82.37 | 84.59 | 77.37 | 2.70 | −6.07 |
Ⅳ | 1670.20 | 1480.32 | 1354.03 | −11.37 | −18.93 | |
Ⅵ | 154.55 | 154.50 | 176.98 | −0.03 | 15.51 |
为了解波形钢腹板箱梁的应力在横截面上的分布规律,图5为跨中左截面翘曲正应力、二次剪应力、总剪应力的分布图,其中箭头方向为剪应力的方向. 由图5可以看出:点Ⅰ和点Ⅴ的翘曲正应力绝对值较大,而波形钢腹板几乎不产生翘曲正应力;点Ⅴ和点Ⅵ的二次剪应力绝对值较大,腹板的二次剪应力与自由扭转剪应力叠加后使得总剪应力减小且沿梁高呈均匀分布,而底板的二次剪应力与自由扭转剪应力叠加后使得总剪应力增大,悬臂板壁厚中面上不产生自由扭转剪应力,但存在二次剪应力.
为进一步探索截面几何参数变化对波形钢腹板箱梁约束扭转效应的影响,选取不同高宽比的梁为分析对象. 梁高h分别取值为2.15、4.30 m,即高宽比h/bt分别为0.5、1.0. 根据图2分别改变tw、bf,腹板厚度tw以2 mm为步长从8 mm增大至22 mm,悬臂板单侧宽度bf从0开始以0.215 m为步长增大至4.300 m,即悬臂板宽度与闭合箱壁宽度的比值(简称悬臂板宽度比,bf /bt)从0开始以0.05为步长增大至1.00.
为比较乌-Ⅱ理论与本文方法计算的翘曲系数的差异,图6绘制了乌-Ⅱ理论与本文方法计算的翘曲系数比值
由于点Ⅰ和点Ⅴ的翘曲正应力绝对值较大,图7绘制了跨中左截面点Ⅰ和点Ⅴ翘曲正应力随tw、bf /bt的变化曲线. 由图7可以看出:点Ⅰ和点Ⅴ的翘曲正应力绝对值均随腹板厚度的增大而减少,随悬臂板宽度比的增大先增大后减少,bf /bt=0.10时开始减小;当bf /bt>0.10时,悬臂板宽度增大可以减少翘曲正应力.
图8绘制了跨中左截面点 Ⅳ 和点 Ⅵ 的总剪应力随tw、bf/bt的变化曲线. 由图8可以看出:高宽比较小的梁的总剪应力明显大于高宽比较大的梁;对高宽比较大的梁,当tw=8 mm和bf/bt=0.30时,点Ⅵ总剪应力分别达到点 Ⅳ总剪应力的34.4%、30.2%,而混凝土与钢板的剪切模量比为0.17,底板总剪应力显然不容忽视,经分析,悬臂板根部的总剪应力同样不容忽视. 因此,在波形钢腹板箱梁设计时,应对顶底板进行精细化设计,防止斜裂缝的产生;随着悬臂板宽度的增大,当bf /bt>0.30时,总剪应力几乎不变;对高宽比较小的梁的总剪应力随腹板厚度增大而减小,而对高宽比较大的梁,腹板厚度对总剪应力的影响不显著.
1) Reissner原理将位移和应力同时作为变分参量,考虑约束扭转全部次生剪应力对中面剪切变形的影响来推导波形钢腹板箱梁约束扭转控制微分方程,理论更准确. 本文推导的约束扭转控制微分方程表达式与乌-Ⅱ理论的相同,但翘曲系数公式不同. 数值算例的结果表明,相对于乌-Ⅱ理论,本文方法的计算结果与有限元解吻合更好,验证了本文方法的合理性.
2) 按乌-Ⅱ理论与本文方法计算的翘曲系数的比值表征了两种理论计算结果的差异,翘曲系数比值偏离1越远,则计算结果的差异越大. 两种方法计算的翘曲系数受腹板厚度和悬臂板宽度的影响显著,翘曲系数比值最大可达到4.70.
3) 波形钢腹板箱梁发生约束扭转时,波形钢腹板几乎不承担翘曲正应力,主要承担剪应力,混凝土顶底板既承担翘曲正应力又承担剪应力. 底板总剪应力可达到腹板总剪应力的34.4%,而混凝土相对钢板的抗剪性能较差,应对顶底板予以重视,并对其精细化设计,防止主拉应力超限引起斜裂缝.
4) 波形钢腹板厚度增大可以减小约束扭转翘曲正应力;随着悬臂板宽度的增大,悬臂板宽度比大于0.10时,翘曲正应力减小,而当悬臂板宽度比大于0.30时,总剪应力几乎无变化.
[1] |
MARCHELLI M, COLTRINARI G, ALFARO DEGAN G, et al. Towards a procedure to manage safety on construction sites of rockfall protective measures[J]. Safety Science, 2023, 168: 106307.1-106307.8.
|
[2] |
中华人民共和国住房和城乡建设部. 铁路工程结构可靠性设计统一标准:GB 50216—2019[S]. 北京:中国计划出版社,2019.
|
[3] |
中华人民共和国交通运输部. 公路工程结构可靠性统一标准:JTG 2120—2020[S]. 北京:人民交通出版社,2020.
|
[4] |
国家铁路局. 铁路隧道设计规范:TB 10003—2016[S]. 北京:中国铁道出版社, 2017.
|
[5] |
中华人民共和国交通运输部. 公路隧道设计规范 第一册 土建工程:JTG 3370.1—2018[S]. 北京:人民交通出版社, 2019.
|
[6] |
中华人民共和国住房和城乡建设部. 建筑结构可靠性设计统一标准:GB 50068—2018[S]. 北京:中国建筑工业出版社.
|
[7] |
铁道第二勘测设计院. 铁路工程设计技术手册. 隧道[M]. 北京:中国铁道出版社,1999.
|
[8] |
曾廉. 明洞顶设计荷载的研究[J]. 铁路标准设计通讯,1974(7): 3-17.
|
[9] |
ILLEDITSCH M, PREH A. Determination of meaningful block sizes for rockfall modelling[J]. Natural Hazards, 2024, 120(6): 5685-5710. doi: 10.1007/s11069-024-06432-4
|
[10] |
STELZER G, BICHLER A. ONR 24810–A comprehensive guideline for building better rockfall protection structures[C]//64th Annual Highway Geology Symposium. New Hampshire: Highway Geology Symposium, 2013.
|
[11] |
MARCHELLI M, DE BIAGI V, PEILA D. Reliability-based design of rockfall passive systems height[J]. International Journal of Rock Mechanics and Mining Sciences, 2021, 139: 104664.1-104664.9.
|
[12] |
MARCHELLI M, DE BIAGI V, PEILA D. Reliability-based design of protection net fences: influence of rockfall uncertainties through a statistical analysis[J]. Geosciences, 2020, 10(8): 280.1-280.24.
|
[13] |
DE BIAGI V, MARCHELLI M, PEILA D. Reliability analysis and partial safety factors approach for rockfall protection structures[J]. Engineering Structures, 2020, 213: 110553.1-110553.1.
|
[14] |
Austrian Standards Institute. Technical protection against rockfall-terms and definitions, effects of actions, design, monitoring and maintenance: ONR 24810[S]. Vienna: [s.n.], 2017.
|
[15] |
MELZNER S, ROSSI M, GUZZETTI F. Impact of mapping strategies on rockfall frequency-size distributions[J]. Engineering Geology, 2020, 272: 105639.1-105639.11.
|
[16] |
LAIMER H J. Determination of rockfall design blocks in Upper Triassic limestones and Dolomites (Dachstein Formation, Northern Calcareous Alps)[J]. Bulletin of Engineering Geology and the Environment, 2020, 79(3): 1581-1590. doi: 10.1007/s10064-019-01640-w
|
[17] |
DE BIAGI V, LIA NAPOLI M, BARBERO M, et al. Estimation of the return period of rockfall blocks according to their size[J]. Natural Hazards and Earth System Sciences, 2017, 17(1): 103-113. doi: 10.5194/nhess-17-103-2017
|
[18] |
CURCEAC S, ATKINSON P M, MILNE A, et al. An evaluation of automated GPD threshold selection methods for hydrological extremes across different scales[J]. Journal of Hydrology, 2020, 585: 124845.1-124845.15.
|
[19] |
STORMY A. MATLAB[M]. Fifth edition. Boston: [s.n.], 2019.
|
[20] |
COLES S. An Introduction to Statistical Modeling of Extreme Values[M]. London: Springer, 2001.
|
[21] |
史道济. 实用极值统计方法[M]. 天津:天津科学技术出版社,2006.
|
[22] |
仇学艳,王超,秦崇仁. 阈值法在河海工程设计中的应用[J]. 水利学报,2001(8): 32-37. doi: 10.3321/j.issn:0559-9350.2001.08.006
QIU Xueyan, WANG Chao, QIN Chongren. The application of threshold procedure in design of river and coastal engineering[J]. Shuili Xuebao, 2001(8): 32-37. doi: 10.3321/j.issn:0559-9350.2001.08.006
|
[23] |
PHILIP J, DAVID R, JENNY W, et al. Uncertainties in return from extreme value analysis of peaks over threshold using the generalised Pareto distribution[J]. Ocean Engineering, 2021, 220: 107725.1-107725.17.
|
[24] |
THOMAS R K, ILARIA P. Use of peak over threshold data for flood frequency estimation: An application at the UK national scale.[J]. Jounal of Hydrology, 2023, 626(5): 130235.1-130235.31.
|
[25] |
中华人民共和国住房和城乡建设部. 建筑抗震设计规范GB 50011—2010(2016版)[S]. 北京:中国建筑工业出版社,2016.
|
名称 | 计算点 | 有限元 解/kPa | 本文方 法/kPa | 乌-Ⅱ 理论/kPa | 相对误差/% | |
本文 方法 | 乌-Ⅱ 理论 | |||||
{\sigma} | Ⅰ | 54.06 | 57.35 | 63.61 | 6.09 | 17.67 |
Ⅱ | −27.83 | −25.42 | −28.18 | −8.66 | 1.28 | |
Ⅴ | −75.83 | −81.24 | −90.10 | 7.13 | 18.82 | |
{\tau _{\text{z}}} | Ⅲ | 82.37 | 84.59 | 77.37 | 2.70 | −6.07 |
Ⅳ | 1670.20 | 1480.32 | 1354.03 | −11.37 | −18.93 | |
Ⅵ | 154.55 | 154.50 | 176.98 | −0.03 | 15.51 |