CN103593553B - 一种盾构隧道管片非均质等效梁单元模型结构计算方法 - Google Patents
一种盾构隧道管片非均质等效梁单元模型结构计算方法 Download PDFInfo
- Publication number
- CN103593553B CN103593553B CN201310532655.9A CN201310532655A CN103593553B CN 103593553 B CN103593553 B CN 103593553B CN 201310532655 A CN201310532655 A CN 201310532655A CN 103593553 B CN103593553 B CN 103593553B
- Authority
- CN
- China
- Prior art keywords
- lambda
- section
- equivalent
- max
- sigma
- 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
Abstract
本发明公开了一种盾构隧道管片非均质等效梁单元模型结构计算方法,方法包括步骤:(1)收集基本资料,建立惯用法模型并计算,提取各接头位置截面内力值;(2)计算等效截面力学参数;(3)建立非均质等效梁单元模型,并进行结构计算;(4)在次计算其截面等效参数;(5)比较步骤(3)输入的等效截面力学参数与步骤(4)输出的截面力学参数,是否满足误差精度要求?若不满足,则回到步骤(3),再次迭代计算,重复步骤(3)~步骤(5)。本发明采用的计算模型更接近真实,计算结果更加准确。
Description
技术领域:
本发明涉及盾构隧道管片结构的计算方法,特别是涉及一种盾构隧道管片非均质等效梁单元模型结构计算方法。
背景技术:
在盾构隧道设计施工以及运营过程中,如何进行盾构隧道管片结构的内力计算是当前的研究热点和难点。近几十年来,伴随城市轨道交通的大规模发展,针对盾构隧道管片结构的计算模型研究取得了长足的发展和进步,衍生出了包括匀质圆环模型、等效刚度圆环模型、自由铰圆环模型、梁—弹簧模型和梁—接头模型在内的多种计算结构模型。然而,由于结构本身的以及赋存环境的复杂性,当前应用在盾构隧道管片结构力学计算的几种模型中,仍存在一定的不足和缺陷,如:匀质圆环模型完全忽略了管片接头的影响,显然与实际情况不符,导致硬质地层内力计算结果偏大,造成浪费;软弱地层变形偏小,导致安全风险增大。而修正的惯用法模型本身仍是通过人为折减整个管片环接头刚度以及弯矩增大系数来体现,具有较大的随意性和不确定性;自由铰圆环模型将接头作为单铰结构,不仅没有反映出接头的截面传力性质,而且其本身是非静定体系,需借助隧道围岩的支撑作用才能进行静力求解,而隧道围岩的支撑力依然没有确切可靠的计算参数,故该模型也具有较大的局限性和不确定性;梁—弹簧和梁—接头模型能够较好地反映管片接头对整个环结构的刚度影响,但其弹簧系数的取得依然没有成熟的理论基础,大多是需要通过具体试验获取或工程经验来设置,其计算结果不可避免地存在一定的任意性。
中国发明专利公开号CN101364241A公开的盾构隧道衬砌连续—非均匀刚度模型结构计算方法,在处理管片接头等效刚度时,仅采用了一个通过接头转角反算得到的综合等效刚度系数,并没有分别考虑等效梁单元受压区高度、宽度以及材料弹性模量的等效性及其实际受力状态,事实上,这三者以及等效刚度与管片结构的实际受力状态息息相关,且三者之间并非为线性变化关系,因此采用综合的等效刚度系数在某种程度上仍不足以表征接头的实际受力状态和以反映接头的刚度影响。
可见,尽管广大科技工作者针对盾构隧道管片结构的内力计算做了大量的研究工作,取得了丰富的成果,但仍存在诸多难解之题,仍需进一步广泛研究。
发明内容:
本发明的目的在于克服现有技术的不足,提供一种新的计算模型来进行盾构隧道管片的结构计算,使得计算结果更加准确。
为了解决上述技术问题,本发明一种盾构隧道管片非均质等效梁单元模型结构计算方法,包括以下步骤:
(1)收集基本资料,建立惯用法模型并计算,提取各接头位置截面内力值,即弯矩M和轴力N;
(2)根据接头等效截面受力类型及内力值大小,计算等效截面力学参数,即混凝土等效受压区高度h′、等效受压区宽度b′及等效弹性模量E′e;
(3)建立非均质等效梁单元模型,将步骤(2)得到的等效截面力学参数作为初始值输入该模型,并进行结构计算;
(4)提取步骤(3)计算得到的等效梁单元截面内力值,再次计算其截面等效参数;
(5)比较步骤(3)输入的等效截面力学参数与步骤(4)输出的截面力学参数,看是否满足误差精度要求,若不满足,则回到步骤(3),再次迭代计算,重复步骤(3)~步骤(5);若满足,则输出最终计算结果。
进一步地,所述步骤(2)在不同受力类型条件下,等效力学参数计算公式如下:
①正弯矩小偏心受压
其中,λ1=nT0+N;
式中,b′、h′、E′e分别为等效梁单元截面的等效受压区宽度、高度、弹性模量及等效截面计算受压区高度;b、h、M、N、n、σc,max、σc,min、分别为原管片接头截面宽度、管片厚度、弯矩、轴力、单位宽度管片内横向螺栓根数、管片上缘混凝土应力、管片下缘混凝土应力;Ee为接头梁单元等效抗压弹性模量,采用Ee=EcEp(l1+l2)/(Ecl1+Epl2)计算,Ec、Ep分别为管片混凝土和止水橡胶条(传力衬垫)的弹性模量;T0为螺栓的初始预紧力;h0为螺栓轴线离管片內缘的距离;K为螺栓的抗拉刚度;
②正弯矩大偏心受压
式中,x′为等效梁单元截面的等效截面计算受压区高度;αcs为螺栓横截面积与管片内侧钢筋面积之比;Tb为考虑了初始预紧力后的原接头截面螺栓拉力,Tb=T0+(h-h0-x)σc,maxK/EeL;L为螺栓的有效计算长度;x为原管片接头截面计算受压区高度,按式Ax2+Bx+C=0——关于受压区高度x的一元二次方程求得的在[0,h]范围的解取值;σc,max按式
③负弯矩小偏心受压
其中,λ1=nT0+N;
④负弯矩大偏心受压
负弯矩大偏心受压条件下的等效截面高度和宽度的计算方法同正弯矩大偏心受压情况。
进一步地,负弯矩大偏心受压中的原接头截面混凝土受压区高度x和混凝土下缘最大应力σc,max分别按以下两种情况计算确定:
第一种情况:原接头截面混凝土受压区高度x由式计算确定,且x∈(h0,h],混凝土下缘最大应力σc,max由式计算确定,
其中, λ5=λ2h0;
第二种情况:原接头截面混凝土受压区高度x由式λ3x2-λ4x+λ5=0计算确定,且x∈[0,h0],求得x后,混凝土下缘最大应力σc,max由式 计算确定,
其中, λ3=bλ1; λ5=λ2h0。
与现有技术相比,本发明的有益效果是:考虑了管片接头及螺栓初始预紧力对结构体系的刚度影响,使得结构计算结果与实际情况更为吻合;分别针对等效梁单元的等效高度、宽度和弹性模量进行了推导分析,并给出了具体的计算公式,摒弃了采用单一综合等效刚度的做法,使得等效力学截面的受力特征更为明确、具体,更能反映接头及螺栓预紧力对接头刚度的影响。
附图说明:
图1为现有的盾构隧道管片惯用法结构计算模型。
图2为本发明正弯矩小偏心受压条件接头截面力学分析图。
图3为本发明接头截面等效弹性模量分析图。
图4为本发明正弯矩小偏心条件下等效梁单元截面受力分析图。
图5为本发明正弯矩大偏心受压条件接头截面力学分析图。
图6为本发明正弯矩大偏心条件下等效梁单元截面受力分析图。
图7为本发明负弯矩小偏心受压条件接头截面力学分析图。
图8为本发明负弯矩小偏心条件下等效梁单元截面受力分析图。
图9为本发明负弯矩大偏心受压条件接头截面力学分析图。
图10为本发明负弯矩大偏心受压条件接头截面力学分析图。
图11为本发明盾构隧道管片接头非均质等效梁单元模型。
图12为本发明盾构隧道管片结构非均质等效梁结构计算模型。
图13为本发明盾构隧道管片结构非均质等效梁模型计算流程图。
具体实施方式:
下面结合附图和具体实施方式对本发明作进一步描述:
(1)采用惯用法模型(图1)建立盾构隧道管片结构计算模型,按相关规范要求计算各荷载即模型边界,进行计算,提取各接头截面的内力值——弯矩M和轴力N。
(2)根据接头等效截面受力类型及内力值大小,计算等效截面力学参数,即混凝土等效受压区高度h′、等效受压区宽度b′及等效弹性模量E′e。不同受力类型条件下,等效力学参数计算公式如下:
①正弯矩小偏心受压(图2、图4)
其中,λ1=nT0+N;
式中,b′、h′、E′e分别为等效梁单元截面的等效受压区宽度、高度、弹性模量及等效截面计算受压区高度;b、h、M、N、n、σc,max、σc,min、分别为原管片接头截面宽度、管片厚度、弯矩、轴力、单位宽度管片内横向螺栓根数、管片上缘混凝土应力、管片下缘混凝土应力;Ee为接头梁单元等效抗压弹性模量(图3),采用Ee=EcEp(l1+l2)/(Ecl1+Epl2)计算,Ec、Ep分别为管片混凝土和止水橡胶条(传力衬垫)的弹性模量;T0为螺栓的初始预紧力;h0为螺栓轴线离管片內缘的距离;K为螺栓的抗拉刚度。上述参数符号意义同样适用于以下内容。
②正弯矩大偏心受压(图5、图6)
式中,x′为等效梁单元截面的等效截面计算受压区高度;αcs为螺栓横截面积与管片内侧钢筋面积之比;Tb为考虑了初始预紧力后的原接头截面螺栓拉力,Tb=T0+(h-h0-x)σc,maxK/EeL;L为螺栓的有效计算长度;x为原管片接头截面计算受压区高度,按式(3)——关于受压区高度x的一元二次方程求得的在[0,h]范围的解取值。σc,max按式(4)计算。
Ax2+Bx+C=0 (3)
③负弯矩小偏心受压(图7、图8)
其中,λ1=nT0+N;
④负弯矩大偏心受压(图9、图10)
负弯矩大偏心受压条件下的等效截面高度和宽度的计算方法同正弯矩大偏心受压情况。式(6)中原接头截面混凝土受压区高度x和混凝土下缘最大应力σc,max分别按以下两种情况计算确定。
第一种情况(图9a):原接头截面混凝土受压区高度x由式(7)计算确定,且x∈(h0,h],混凝土下缘最大应力σc,max由式(8)计算确定。
其中, λ5=λ2h0。
第二种情况(图9b):原接头截面混凝土受压区高度x由式(9)计算确定,且x∈[0,h0],求得x后,混凝土下缘最大应力σc,max由式(10)计算确定。
λ3x2-λ4x+λ5=0 (9)
其中, λ3=bλ1; λ5=λ2h0。
实际计算过程中为一试算过程,首先按第一种情况计算得到x值,若x∈(h0,h],即为第一种情况,若则再按第二种情况计算。
(3)将管片接头视为一受压区高度、宽度及弹性模量均不同于管片主截面的梁单元(图11),等效梁单元长度为L=2l1+l2,建立非均质等效梁单元模型,将步骤(2)得到的等效截面力学参数作为初始值输入该模型(图12),并进行结构计算。
(4)提取步骤(3)计算得到的等效梁单元截面内力值,再次计算其截面等效参数。
(5)比较步骤(3)输入的等效截面力学参数与步骤(4)输出的截面力学参数,是否满足误差精度要求?若不满足,则回到步骤(3),再次迭代计算,重复步骤(3)~步骤(5);若满足,则输出最终计算结果。
Claims (2)
1.一种盾构隧道管片非均质等效梁单元模型结构计算方法,其特征在于,它由以下步骤实现:
(1)收集基本资料,建立惯用法模型并计算,提取各接头位置截面内力值,即弯矩M和轴力N;
(2)根据接头等效截面受力类型及内力值大小,计算等效截面力学参数,即混凝土等效受压区高度h′、等效受压区宽度b′及等效弹性模量E′e;
(3)建立非均质等效梁单元模型,将步骤(2)得到的等效截面力学参数作为初始值输入该模型,并进行结构计算;
(4)提取步骤(3)计算得到的等效梁单元截面内力值,再次计算其截面等效参数;
(5)比较步骤(3)输入的等效截面力学参数与步骤(4)输出的截面力学参数,看是否满足误差精度要求,若不满足,则回到步骤(3),再次迭代计算,重复步骤(3)~步骤(5);若满足,则输出最终计算结果;
所述步骤(2)在不同受力类型条件下,等效力学参数计算公式如下:
①正弯矩小偏心受压
其中,λ1=nT0+N;
式中,b′、h′、E′e分别为等效梁单元截面的等效受压区宽度、高度、弹性模量及等效截面计算受压区高度;b、h、M、N、n、σc,max、σc,min、分别为原管片接头截面宽度、管片厚度、弯矩、轴力、单位宽度管片内横向螺栓根数、管片上缘混凝土应力、管片下缘混凝土应力;Ee为接头梁单元等效抗压弹性模量,T0为螺栓的初始预紧力;h0为螺栓轴线离管片內缘的距离;K为螺栓的抗拉刚度;
②正弯矩大偏心受压
式中,x′为等效梁单元截面的等效截面计算受压区高度;αcs为螺栓横截面积与管片内侧钢筋面积之比;Tb为考虑了初始预紧力后的原接头截面螺栓拉力,Tb=T0+(h-h0-x)σc, maxK/EeL;L为螺栓的有效计算长度;x为原管片接头截面计算受压区高度,按式Ax2+Bx+C=0——关于受压区高度x的一元二次方程求得的在[0,h]范围的解取值;σc,max按式
或计算其中,A=λ2λ5;B=λ3λ5-λ1λ6;C=λ4λ5-λ1λ7;
λ2=-b/6;
λ5=N+nT0;
③负弯矩小偏心受压
其中,λ1=nT0+N;
④负弯矩大偏心受压
负弯矩大偏心受压条件下的等效截面高度和宽度的计算方法同正弯矩大偏心受压情况。
2.根据权利要求1所述的一种盾构隧道管片非均质等效梁单元模型结构计算方法,其特征在于,所述负弯矩大偏心受压中的原接头截面混凝土受压区高度x和混凝土下缘最大应力σc,max分别按以下两种情况计算确定:
第一种情况:原接头截面混凝土受压区高度x由式计算确定,且x∈(h0,h],混凝土下缘最大应力σc,max由式计算确定,
其中,
λ5=λ2h0;
第二种情况:原接头截面混凝土受压区高度x由式λ3x2-λ4x+λ5=0计算确定,且x∈[0,h0],求得x后,混凝土下缘最大应力σc,max由式
计算确定,
其中,λ3=bλ1;
λ5=λ2h0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310532655.9A CN103593553B (zh) | 2013-11-01 | 2013-11-01 | 一种盾构隧道管片非均质等效梁单元模型结构计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310532655.9A CN103593553B (zh) | 2013-11-01 | 2013-11-01 | 一种盾构隧道管片非均质等效梁单元模型结构计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103593553A CN103593553A (zh) | 2014-02-19 |
CN103593553B true CN103593553B (zh) | 2017-01-25 |
Family
ID=50083691
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310532655.9A Active CN103593553B (zh) | 2013-11-01 | 2013-11-01 | 一种盾构隧道管片非均质等效梁单元模型结构计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103593553B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105178158B (zh) * | 2015-05-29 | 2017-03-01 | 广西交通科学研究院 | 大跨径混凝梁桥腹板横向钢筋配置方法 |
CN106529004B (zh) * | 2016-10-28 | 2019-06-04 | 中铁二院工程集团有限责任公司 | 一种隧道钢筋混凝土衬砌截面钢筋面积的设计方法 |
CN107301305B (zh) * | 2017-07-31 | 2021-03-09 | 中铁隧道局集团有限公司 | 一种实现管片三维受载的结构力学性能的建模方法 |
CN109682557B (zh) * | 2018-12-03 | 2020-11-10 | 北京市市政工程研究院 | 一种人行天桥栏杆结构承载能力试验评定方法 |
CN110909406B (zh) * | 2019-11-20 | 2022-08-30 | 西南交通大学 | 基于模拟接头模型的盾构隧道内力计算方法 |
CN113361169B (zh) * | 2021-06-10 | 2022-04-22 | 安徽省建筑科学研究设计院 | 一种地表突发堆载引起盾构隧道纵向变形的高效预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101364241A (zh) * | 2007-08-08 | 2009-02-11 | 同济大学 | 盾构隧道衬砌连续-非均匀刚度模型结构计算方法 |
CN101363323A (zh) * | 2007-08-08 | 2009-02-11 | 同济大学 | 盾构隧道预应力衬砌设计方法 |
CN101667220A (zh) * | 2009-09-22 | 2010-03-10 | 同济大学 | 盾构隧道接头构造的模拟方法 |
CN102306225A (zh) * | 2011-09-27 | 2012-01-04 | 上海大学 | 多线叠交隧道施工过程及对隧道变形影响数值模拟方法 |
-
2013
- 2013-11-01 CN CN201310532655.9A patent/CN103593553B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101364241A (zh) * | 2007-08-08 | 2009-02-11 | 同济大学 | 盾构隧道衬砌连续-非均匀刚度模型结构计算方法 |
CN101363323A (zh) * | 2007-08-08 | 2009-02-11 | 同济大学 | 盾构隧道预应力衬砌设计方法 |
CN101667220A (zh) * | 2009-09-22 | 2010-03-10 | 同济大学 | 盾构隧道接头构造的模拟方法 |
CN102306225A (zh) * | 2011-09-27 | 2012-01-04 | 上海大学 | 多线叠交隧道施工过程及对隧道变形影响数值模拟方法 |
Non-Patent Citations (2)
Title |
---|
盾构隧道非均质等效梁模型的建立与分析;高云龙等;《铁道学报》;20130715;第35卷(第7期);摘要,正文第1-3节 * |
迎坡条件下盾构隧道开挖面极限支护力计算与分析;雷明锋等;《岩土工程学报》;20100315;第32卷(第3期);第488-492页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103593553A (zh) | 2014-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103593553B (zh) | 一种盾构隧道管片非均质等效梁单元模型结构计算方法 | |
CN101364241B (zh) | 盾构隧道衬砌连续-非均匀刚度模型结构计算方法 | |
CN102592029B (zh) | 湿陷性黄土路基工后沉降的分析预测方法 | |
Zhang et al. | A novel displacement back analysis method considering the displacement loss for underground rock mass engineering | |
CN102446241B (zh) | 一种翼面结构刚度仿真方法 | |
Zhang et al. | A numerical study on choosing the best configuration of the blade for vertical axis wind turbines | |
CN104252576B (zh) | 一种地铁隧道岩土体等效计算参数反演方法 | |
Ho et al. | Numerical flow analysis for spillways | |
CN103399992A (zh) | 一种基于可靠寿命的结构耐久性优化设计方法 | |
CN103778298A (zh) | 改进的模拟多孔介质中二维水流运动的多尺度有限元方法 | |
Chakraborty et al. | An efficient algorithm for building locally refined hp–adaptive H-PCFE: Application to uncertainty quantification | |
CN104794356B (zh) | 一种高耸升船机塔柱结构抗风分析方法 | |
Wu et al. | An approach for deformation modulus mechanism of super-high arch dams | |
Zhou et al. | A novel method for the dynamic reliability analysis of slopes considering dependent random parameters via the direct probability integral method | |
Liu et al. | Cracking risk analysis of face slabs in concrete face rockfill dams during the operation period | |
Zhu et al. | A new dam reliability analysis considering fluid structure interaction | |
Tseng | Kinematic wave computation using an efficient implicit method | |
CN105631093A (zh) | 一种基于m-bswa多目标优化的机械结构设计方法 | |
Basile | A practical method for the non-linear analysis of piled rafts | |
Zhang et al. | Using softened contact relationship describing compressible membrane in FEA of spiral case structure | |
Kissami et al. | Numerical assessment of criteria for mesh adaptation in the finite volume solution of shallow water equations | |
Fu et al. | Slope stability analysis based on big data and convolutional neural network | |
Cai et al. | Numerical simulation of overtopping breach processes caused by failure of landslide dams | |
CN108229069B (zh) | 一种提高液氢储罐防旋转装置有限元分析精度的方法 | |
CN106874611A (zh) | 一种基于超体积迭代策略的含区间参数结构响应区间的分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |