CN117216862B - 一种基于纤维梁单元的三维盾构隧道动力分析模型 - Google Patents
一种基于纤维梁单元的三维盾构隧道动力分析模型 Download PDFInfo
- Publication number
- CN117216862B CN117216862B CN202311483822.5A CN202311483822A CN117216862B CN 117216862 B CN117216862 B CN 117216862B CN 202311483822 A CN202311483822 A CN 202311483822A CN 117216862 B CN117216862 B CN 117216862B
- Authority
- CN
- China
- Prior art keywords
- joint
- segment
- concrete
- strain
- gap
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 45
- 239000000835 fiber Substances 0.000 title claims abstract description 24
- 239000004567 concrete Substances 0.000 claims description 48
- 239000000463 material Substances 0.000 claims description 21
- 238000006073 displacement reaction Methods 0.000 claims description 18
- 238000005452 bending Methods 0.000 claims description 15
- 230000003993 interaction Effects 0.000 claims description 11
- 239000011241 protective layer Substances 0.000 claims description 9
- 238000010008 shearing Methods 0.000 claims description 7
- 239000010410 layer Substances 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 4
- 238000007906 compression Methods 0.000 claims description 4
- 229910001294 Reinforcing steel Inorganic materials 0.000 claims description 3
- 229910000831 Steel Inorganic materials 0.000 claims description 3
- 239000012792 core layer Substances 0.000 claims description 3
- 239000010959 steel Substances 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 18
- 238000004364 calculation method Methods 0.000 abstract description 16
- 238000012360 testing method Methods 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract 1
- 230000006399 behavior Effects 0.000 description 22
- 239000002689 soil Substances 0.000 description 11
- 230000008569 process Effects 0.000 description 8
- 238000013016 damping Methods 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012821 model calculation Methods 0.000 description 2
- 239000011150 reinforced concrete Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 241000272525 Anas platyrhynchos Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- 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
- Lining And Supports For Tunnels (AREA)
Abstract
本申请公开了一种基于纤维梁单元的三维盾构隧道动力分析模型,属于隧道地震安全韧性评估技术领域,具体包括:设置模型的单元和在分析模型中引入本构关系描述管片及接头的非线性力学行为。本申请可以准确模拟管片拼装方式,结构非线性力学行为,提高动力分析计算效率及精度;通过与试验结果对比,验证了模型的可行性及准确性;通过多个地震动分析实例,验证模型的适用性、稳定性、高效性、必要性。
Description
技术领域
本申请涉及隧道地震安全韧性评估技术领域,具体涉及一种基于纤维梁单元的三维盾构隧道动力分析模型。
背景技术
盾构隧道有助于通过掘进机高效、自动化地隧道施工,并已广泛应用于城市生命线基础设施,例如地铁系统和地下公路。因此,确保盾构隧道的地震安全对于维持城市的韧性至关重要。盾构隧道等地下基础设施的脆弱性极大地影响了城市的抗震能力。
如何准确评估过去几十年普遍建设的盾构隧道的地震易损性是一个实际问题。这个实际问题变得具有挑战性,因为盾构隧道是钢筋混凝土管片、螺栓、垫片和间隙等的组合体。每个部件的机械行为可能是高度非线性的,并且涉及与其他部件的复杂相互作用,如通过离心机、振动台和其他静载荷试验都观察到这一现象。例如,接缝的存在极大地改变了隧道衬砌的内力分布,也就是说,忽略它可能会导致隧道衬砌的抗震需求要么高估要么低估;接缝的开裂和破损标志着整体衬砌的破损程度,应计入隧道衬砌的性能指标或破损指标。因此,盾构隧道的理想模型应该能够捕获外部荷载下的上述高度非线性响应。
现今,在开发盾构隧道地震分析数值模型方面付出了巨大的努力,这些模型可以分为两类:整体方法和复杂方法。前一种方法通常将盾构隧道衬砌视为线弹性均匀梁或壳单元,其中接缝的贡献通过所谓的等效刚度来近似。 显然,这种方法可能能够再现盾构隧道的整体变形,但很难估计它们的损坏或破坏模式。一些研究人员试图通过引入独立的非线性弹簧来代表关节来完善整体方法。尽管如此,这些改进的方法无法捕捉装配部件机械状态的演变,例如螺栓的屈服、间隙的接触和挤压。后一种方法可以用实体单元和接触单元复杂地表示盾构隧道。此类模型可以很好地模拟盾构隧道的非线性响应,但包含巨大的自由度,导致计算成本非常高,尤其是在进行三维动力时程分析过程中,分析成本呈指数增长。此外大部分现有地震分析模型由于计算效率的限制,多数面向解决二维地震响应问题开发。然而盾构隧道作为一种拼装结构,其本身的力学行为不仅受到拼装方式的影响(错缝拼装或通缝拼装),同时在多向地震作用下结构损伤更明显。迄今为止,仍然缺乏一个既具有可以反映结构拼装方式与多向地震动影响的,同时可以保证计算高精度、高计算效率的盾构隧道数值分析模型。
发明内容
本申请实施例的目的是提供一种基于纤维梁单元的三维盾构隧道动力分析模型,其可以准确模拟管片拼装方式,结构非线性力学行为,提高动力分析计算效率及精度,从而可以解决背景技术中涉及的至少一个技术问题。
为了解决上述技术问题,本申请是这样实现的:
本申请实施例提供了一种基于纤维梁单元的三维盾构隧道动力分析模型,包括:
设置模型的单元,包括:基于OpenSees中的纤维梁单元模拟管片非线性轴弯相互作用行为,采用零长度单元模拟管片的纵向接头轴弯相互作用,采用刚臂单元及零长度单元模拟管片纵向延伸及环间接头轴弯相互作用,采用Section Aggregator命令在接头单元中耦合剪切弹簧本构模拟接头的剪切行为;
在分析模型中引入本构关系描述管片及接头的非线性力学行为,包括:
引入三种本构关系来描述纤维梁管片的非线性力学行为:保护层混凝土单轴本构、核心混凝土单轴本构及纵向钢筋单轴本构;
引入了五种本构关系描述接头的非线性力学行为:间隙-间隙接触、核心混凝土接触、垫片-垫片接触、螺栓及剪切弹簧的刚度。
可选的,核心层混凝土和保护层混凝土均采用OpenSees中的Concrete 02单轴材料进行模拟;钢筋采用Steel02 单轴材料进行模拟。
可选的,间隙-间隙接触,由完美弹塑性间隙的单轴材料来模拟,表示为:
;
;
;
式中,σ是压缩位移为u时材料中的应力;ugc1是间隙的宽度;ugc2和ugc3分别是间隙-间隙接触的屈服和极限位移;εc0为达到保护层混凝土抗压强度fc0时的应变;εcu为保护层混凝土极限压缩应变;lc是接头的标称长度,具体取决于接头的配置。
可选的,核心混凝土接触,采用弹塑性单轴材料模拟,表示为:
;
;
;
式中,uc1和uc2分别是核心混凝土接触的屈服和极限位移;为核心混凝土抗压强度/>时的应变;/>是核心混凝土的极限压应变;
核心混凝土的抗压强度由下式计算:
;
式中,hcore为核心混凝土厚度,h为节段接缝总厚度;
核心混凝土的抗压强度应变 和极限压缩应变/>分别由下式计算:
;
;
;
式中,σ2是约束应力。
可选的,垫片-垫片接触,将垫片的压缩行为简化为双线性模型,表示为:
;
式中,kg1和kg2分别为垫片的弹性模量和硬化模量;ug为垫片的屈服位移。
可选的,螺栓,由仅受拉的单轴材料Steel 02进行模拟,初始应力表示为fb0,屈服强度和极限强度分别表示为fb1和fb2,则屈服及极限位移由下式确定:
;
;
式中,lb是螺栓长度;εb1 和εb2 分别是屈服应变和极限应变。
可选的,剪切弹簧的刚度,由下式表示:
;
式中,Ec为管段的弹性模量;I为隧道管片横截面的惯性矩;L和l分别为段环的宽度和螺栓的长度;n为螺栓数量;Gb为一根螺栓的剪切模量;Sb为一根螺栓的横截面积;m为分段截面形状系数,矩形截面取m=1.5。
本申请具有如下有益效果:
1、可以准确模拟管片拼装方式,结构非线性力学行为,提高动力分析计算效率及精度;
2、通过与试验结果对比,验证了模型的可行性及准确性;
3、通过多个地震动分析实例,验证模型的适用性、稳定性、高效性、必要性。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图,其中:
图1为本申请三维纤维梁单元管片结构示意图;
图2为本申请提供的管片荷载-挠度曲线图;
图3为本申请提供的管片钢筋荷载-应变曲线图;
图4为本申请提供的直螺栓接头弯矩-张开量曲线图;
图5为本申请提供的直螺栓接头弯矩-转角曲线图;
图6为本申请提供的土-盾构隧道模型二维图;
图7为本申请提供的隧道内力分布图;
图8-图10为本申请提供的接头及螺栓变形时程曲线图;
图11为本申请提供的基于纤维梁单元的三维土-隧道模型图;
图12为本申请提供的不同时刻盾构隧道面内弯矩分布图;
图13为本申请提供的不同时刻盾构隧道位移轮廓图;
图14为本申请提供的所选地震动的5%阻尼比的加速度响应谱图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
本申请的说明书和权利要求书中的术语“第一”、“第二”等是用于区别类似的对象,而不用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便本申请的实施例能够以除了在这里图示或描述的那些以外的顺序实施,且“第一”、“第二”等所区分的对象通常为一类,并不限定对象的个数,例如第一对象可以是一个,也可以是多个。此外,说明书以及权利要求中“和/或”表示所连接对象的至少其中之一,字符“/”,一般表示前后关联对象是一种“或”的关系。
请参见图1所示,本申请实施例提供的一种基于纤维梁单元的三维盾构隧道动力分析模型,包括:
设置模型的单元,包括:基于OpenSees中的纤维梁单元模拟管片非线性轴弯相互作用行为,采用零长度单元模拟管片的纵向接头轴弯相互作用,采用刚臂单元及零长度单元模拟管片纵向延伸及环间接头轴弯相互作用,采用Section Aggregator命令在接头单元中耦合剪切弹簧本构模拟接头的剪切行为;
在分析模型中引入本构关系描述管片及接头的非线性力学行为,包括:
引入三种本构关系来描述纤维梁管片的非线性力学行为:保护层混凝土单轴本构、核心混凝土单轴本构及纵向钢筋单轴本构;
引入了五种本构关系描述接头的非线性力学行为:间隙-间隙接触、核心混凝土接触、垫片-垫片接触、螺栓及剪切弹簧的刚度。
具体的,核心层混凝土和保护层混凝土均采用OpenSees中的Concrete 02单轴材料进行模拟,其单调压缩行为由修正的 Kent&Park 模型表示,而其拉伸和循环行为遵循Yassin的模型。 核心混凝土的约束效应集中在一个称为约束系数K的参数中;钢筋采用Steel02单轴材料进行模拟。
需要说明的是,在有间隙的区域,混凝土不会接触,直到间隙闭合,因此,间隙-间隙接触,可以由一种称为完美弹塑性间隙的单轴材料来模拟,表示为:
;
;
;
式中,σ是压缩位移为u时材料中的应力;ugc1是间隙的宽度;ugc2和ugc3分别是间隙-间隙接触的屈服和极限位移;εc0为达到保护层混凝土抗压强度 fc0时的应变;εcu为保护层混凝土极限压缩应变;lc是接头的标称长度,具体取决于接头的配置。
核心混凝土接触,采用弹塑性单轴材料模拟,表示为:
;
;
;
式中,uc1和uc2分别是核心混凝土接触的屈服和极限位移;为核心混凝土抗压强度/>时的应变;/>是核心混凝土的极限压应变;
核心混凝土的抗压强度按照下式计算:
;
式中,hcore为核心混凝土厚度,h为节段接缝总厚度;
核心混凝土的抗压强度应变 和极限压缩应变/>按欧洲规范 2 确定,分别由下式计算:
;
;
;
式中,σ2是约束应力。
需要说明的是,垫片与垫片接触描述了分段接头内两个相邻垫片的接触和分离。这种接触只能传递压应力,因此,垫片-垫片接触,Shen将垫片的压缩行为简化为双线性模型,表示为:
;
式中,kg1和kg2分别为垫片的弹性模量和硬化模量;ug为垫片的屈服位移。
需要说明的是,螺栓由仅受拉的单轴材料Steel02进行模拟,螺栓由螺栓杆和螺母组成,只能传递拉应力。 这里通过仅受拉的单轴材料钢02进行模拟,初始应力表示为fb0,屈服强度和极限强度分别表示为fb1和fb2,则屈服及极限位移由下式确定:
;
;
式中,lb是螺栓长度;εb1 和 εb2 分别是屈服应变和极限应变。
剪切弹簧的刚度,由下式表示:
;
式中,Ec为管段的弹性模量;I为隧道管片横截面的惯性矩;L和l分别为段环的宽度和螺栓的长度;n为螺栓数量;Gb为一根螺栓的剪切模量;Sb为一根螺栓的横截面积;m为分段截面形状系数,矩形截面取m=1.5。
通过与多个足尺试验结果进行对比,验证了本模型的准确性及适用性;在多个试验验证过程中,模型的本构参数完全服从本申请所提出的本构关系及材料属性定义,具体对比结果如图2-图5所示,验证结果表明本模型不仅可以反映钢筋混凝土管片的非线性力学行为,还可以准确捕捉管片接头非线性行为。
本申请提供的模型,可以被应用来研究盾构隧道的二维及三维的地震响应,具体如下:
对于二维地震动模拟而言:
在OpenSees中对盾构隧道进行了一系列共计12次地震分析,以证明所提出的模型在抗震性能评估中的适用性。盾构隧道埋在9m的均质软粘土下,它包含六个段:一个关键段(F)、两个相邻段(L)和三个标准段(B)。相邻标准段的圆心角为67.5°。
盾构隧道的数值模型如图6所示。隧道衬砌通过六十四个基于位移的纤维梁单元进行模拟,而接缝则通过所提出的接缝模型进行建模。土壤域由 FourNodeQuad 单元表示,以模拟平面应变条件。 这些元件的厚度定义为 1.0 m,遵循下式规定:
;
式中,λwave,min是最小波长,采用PIMY本构模型模拟了土体的非线性行为。
在隧道-土壤界面处,滑移条件通过 ZeroLengthContactASDimplex 元件实现。该单元由 STKO开发,是一种遵循 Mohr-Coulomb 准则的基于惩罚的摩擦接触单元。惩罚刚度的值由下式定义:
;
式中,O是刚度矩阵的最大条目的阶数,库仑摩擦系数定义为 0.8。
在土域边界处,分别采用固定边界和人工边界进行重力分析和地震分析。 这里,采用ASDabsorbingBounda-ry单元来模拟人工边界条件。它由粘性缓冲器和剪切土柱组成。前者用于模拟粘性人工边界,而后者用于模拟自由场条件。粘性阻尼器在剪切和法向方向上的阻尼系数为cs=ρbvs,cp=ρbvp,其中ρb是土域密度,vs和vp是边界处速度的剪切和法向分量。总的来讲,基于纤维束的盾构隧道模型开发效率较高。借助OpenSees的预处理软件STKO,它在一个小时内开发完成。
地震分析包括四个步骤:
(1)固定模型边界,定义初始土体模型为弹性土体模型,然后进行初始重力分析;
(2)将土体模型转换为弹塑性模型并进行第二次重力分析;
(3) 移除步骤1中的所有固定边界约束,并用步骤(2)中的节点力替换它们;
(4)沿模型底部添加分布剪力τ= 2csv来模拟地震动,其中v为地震动速度,地震动具体信息如图14所示。
为了消除高频域中出现的寄生振荡,采用结合数值阻尼的TRBDF2积分器对系统进行求解。还采用自适应时间步长来增加收敛概率并加快分析速度。在此过程中,初始时间步长和最大时间步长定义为0.005s,如下式:
;
式中,Fc0是截止频率,Lw是截止频率下的波长。如果收敛失败,则时间步减半。如果在所需迭代次数之前实现收敛,则时间步长会增加1.5倍。
此外,采用3%的瑞利阻尼比来模拟土体非常小的剪应变水平下的能量耗散。质量和刚度比例阻尼系数按照下式进行定义:
;
;
式中,ωi为土层固有频率,n为大于模型基础输入运动基频与土层固有频率之比的最接近的奇数。
总的来说,基于纤维的隧道模型的地震分析计算效率很高。所有12次地震分析均使用顺序 Open-Sees 在配备 Intel(r) Xeon(r) Gold 6248r CPU @ 3.00GHZ 的个人计算机上进行。G12下的计算时间最短为1447秒,G3下的最长计算时间为5244秒,平均计算时间为3029.83秒。
以G1地震动计算结果为例,模型计算获得的隧道内力分布如图7所示,模型计算获得的管片接头变形及螺栓变形时程曲线如图8-图10所示。根据图7模拟结果发现隧道内力,包括轴力和弯矩,在接头周围会显著变化。例如,接头3(J3)处的静态弯矩是其右侧段的6.0倍,是其左侧段的1.3倍。这种变化可以归因于管片与接头的不同刚度,突出了接头在盾构隧道地震分析中的不可或缺的作用。
根据图8-图10模拟结果可得出以下结论,总体上来讲,接头的变形在地震期间会累积,并在地震后展现出显著的残余值。具体来说,接头2(J2)的张开和旋转受到最大影响,幅度分别为1.787毫米和-0.0077弧度。中国的盾构隧道设计标准建议限制接头的张开为2.000毫米。这表明地震后隧道存在严重的渗水风险。此外,接头1(J1)的螺栓应变最大,振幅和残余值分别为1065.4微应变和949.30微应变,分别为屈服值的44.40%和39.55%。总之,接头的潜在损伤进一步证明了在盾构隧道的地震分析中,接头及其组装构件的重要性。
根据上述相同的材料设置,单元选择,及地震动施加方法与边界条件设置方法,本申请还又针对三维地震分析开展系列研究。具体的三维隧道-地层相互作用数值分析模型如图9所示。图10记载了不同时刻盾构隧道弯矩分布图,接头处的内力存在明显突变,这与二维模型计算结果相互印证,强调管片接头在盾构隧道地震分析中的作用。图11记载了所示不同时刻盾构隧道位移轮廓图,计算结果表明,在地震过程中盾构隧道受剪切波的作用存在显著的水平向变形,变形后隧道呈“鸭蛋”状。
以上所有的分析,都表现出本申请提供的模型具有很高的建模效率,计算效率,计算精度,同时避免了收敛性及不稳定向问题,因此本申请提出的基于纤维梁单元的三为盾构隧道地震分析模型有望成为一种盾构隧道地震分析模型新范式。
本申请具有如下有益效果:
1、可以准确模拟管片拼装方式,结构非线性力学行为,提高动力分析计算效率及精度;
2、通过与试验结果对比,验证了模型的可行性及准确性;
3、通过多个地震动分析实例,验证模型的适用性、稳定性、高效性、必要性。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
此外,需要指出的是,本申请实施方式中的方法和系统的范围不限按示出或讨论的顺序来执行功能,还可包括根据所涉及的功能按基本同时的方式或按相反的顺序来执行功能,例如,可以按不同于所描述的次序来执行所描述的方法,并且还可以添加、省去、或组合各种步骤。另外,参照某些示例所描述的特征可在其他示例中被组合。
上面结合附图对本申请的实施例进行了描述,但是本申请并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本申请的启示下,在不脱离本申请宗旨和权利要求所保护的范围情况下,还可做出很多形式,均属于本申请的保护之内。
Claims (5)
1.一种基于纤维梁单元的三维盾构隧道动力分析模型,其特征在于,包括:
设置模型的单元,包括:基于OpenSees中的纤维梁单元模拟管片非线性轴弯相互作用行为,采用零长度单元模拟管片的纵向接头轴弯相互作用,采用刚臂单元及零长度单元模拟管片纵向延伸及环间接头轴弯相互作用,采用Section Aggregator命令在接头单元中耦合剪切弹簧本构模拟接头的剪切行为;
在分析模型中引入本构关系描述管片及接头的非线性力学行为,包括:
引入三种本构关系来描述纤维梁管片的非线性力学行为:保护层混凝土单轴本构、核心混凝土单轴本构及纵向钢筋单轴本构;
引入了五种本构关系描述接头的非线性力学行为:间隙-间隙接触、核心混凝土接触、垫片-垫片接触、螺栓及剪切弹簧的刚度,核心层混凝土和保护层混凝土均采用OpenSees中的Concrete 02单轴材料进行模拟;钢筋采用Steel02 单轴材料进行模拟;间隙-间隙接触,由完美弹塑性间隙的单轴材料来模拟,表示为:
;
;
;
式中,σ是压缩位移为u时材料中的应力;ugc1是间隙的宽度;ugc2和ugc3分别是间隙-间隙接触的屈服和极限位移;εc0为达到保护层混凝土抗压强度fc0时的应变;εcu为保护层混凝土极限压缩应变;lc是接头的标称长度,具体取决于接头的配置。
2.根据权利要求1所述的基于纤维梁单元的三维盾构隧道动力分析模型,其特征在于,核心混凝土接触,采用弹塑性单轴材料模拟,表示为:
;
;
;
式中,uc1和uc2分别是核心混凝土接触的屈服和极限位移;为核心混凝土抗压强度时的应变;/>是核心混凝土的极限压应变;
核心混凝土的抗压强度由下式计算:
;
式中,hcore为核心混凝土厚度,h为节段接缝总厚度;
核心混凝土的抗压强度应变 和极限压缩应变/>分别由下式计算:
;
;
;
式中,σ2是约束应力。
3.根据权利要求2所述的基于纤维梁单元的三维盾构隧道动力分析模型,其特征在于,垫片-垫片接触,将垫片的压缩行为简化为双线性模型,表示为:
;
式中,kg1和kg2分别为垫片的弹性模量和硬化模量;ug为垫片的屈服位移。
4.根据权利要求3所述的基于纤维梁单元的三维盾构隧道动力分析模型,其特征在于,螺栓,由仅受拉的单轴材料Steel 02进行模拟,初始应力表示为fb0,屈服强度和极限强度分别表示为fb1和fb2,则屈服及极限位移由下式确定:
;
;
式中,lb是螺栓长度;εb1 和εb2 分别是屈服应变和极限应变。
5.根据权利要求4所述的基于纤维梁单元的三维盾构隧道动力分析模型,其特征在于,剪切弹簧的刚度,由下式表示:
;
式中,Ec为管段的弹性模量;I为隧道管片横截面的惯性矩;L和l分别为段环的宽度和螺栓的长度;n为螺栓数量;Gb为一根螺栓的剪切模量;Sb为一根螺栓的横截面积;m为分段截面形状系数,矩形截面取m=1.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311483822.5A CN117216862B (zh) | 2023-11-09 | 2023-11-09 | 一种基于纤维梁单元的三维盾构隧道动力分析模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311483822.5A CN117216862B (zh) | 2023-11-09 | 2023-11-09 | 一种基于纤维梁单元的三维盾构隧道动力分析模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117216862A CN117216862A (zh) | 2023-12-12 |
CN117216862B true CN117216862B (zh) | 2024-02-02 |
Family
ID=89049706
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311483822.5A Active CN117216862B (zh) | 2023-11-09 | 2023-11-09 | 一种基于纤维梁单元的三维盾构隧道动力分析模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117216862B (zh) |
Citations (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20020021390A (ko) * | 2002-01-04 | 2002-03-20 | 김용성 | 다층지반의 지진응답해석을 위한 동적 점탄-점소성구성모델 및 이를 이용한 해석방법 |
JP2005076339A (ja) * | 2003-09-02 | 2005-03-24 | Shimizu Corp | 鋼コンクリート合成構造の構造解析方法及び同方法を用いて設計された合成構造。 |
WO2009116048A2 (en) * | 2008-03-20 | 2009-09-24 | Technion Research & Development Foundation Ltd. | A method for cosserat point element (cpe) modeling of nonlinear elastic materials |
CN101586463A (zh) * | 2009-06-18 | 2009-11-25 | 苏州市交通设计研究院有限责任公司 | 一种盾构隧道衬砌设计方法 |
JP2009301149A (ja) * | 2008-06-10 | 2009-12-24 | Kyushu Univ | 非線形の有限要素法による構造解析方法、及びプログラム、記録媒体、シミュレーション装置 |
CN107122536A (zh) * | 2017-04-20 | 2017-09-01 | 济南轨道交通集团有限公司 | 一种区间隧道抗震数值模拟方法 |
WO2019140794A1 (zh) * | 2018-01-16 | 2019-07-25 | 天津大学 | 一种基于量纲分析的盾构掘进总推力计算方法 |
CN111535830A (zh) * | 2020-07-09 | 2020-08-14 | 湖南大学 | 盾构隧道管片超高性能混凝土加固结构及其设计方法 |
CN111597617A (zh) * | 2020-05-18 | 2020-08-28 | 中铁第六勘察设计院集团有限公司 | 一种盾构隧道疲劳寿命预测方法 |
CN111950184A (zh) * | 2020-08-21 | 2020-11-17 | 上海市隧道工程轨道交通设计研究院 | 一种多环盾构衬砌结构的三维仿真计算方法 |
CN112329122A (zh) * | 2021-01-07 | 2021-02-05 | 湖南大学 | 一种旁侧基坑开挖引起的盾构隧道横向变形和内力确定方法 |
CN112632837A (zh) * | 2020-12-31 | 2021-04-09 | 华中科技大学 | 一种地下结构纵向抗震数值确定方法 |
CN113221399A (zh) * | 2021-04-13 | 2021-08-06 | 北京工业大学 | 一种铁路隧道装配式衬砌构件划分方案决策方法 |
CN113283142A (zh) * | 2021-05-31 | 2021-08-20 | 华中科技大学 | 一种基于解析解法分析盾构下穿对既有隧道影响的方法 |
CN113360985A (zh) * | 2021-06-10 | 2021-09-07 | 安徽省建筑科学研究设计院 | 一种新建隧道上穿引起既有盾构隧道纵向变形的高精度预测方法 |
CN113420482A (zh) * | 2021-06-24 | 2021-09-21 | 北京安捷工程咨询有限公司 | 一种基于结构内力监测值的管片荷载正交数值反演方法 |
CN114033405A (zh) * | 2021-09-14 | 2022-02-11 | 山东交通职业学院 | 一种高强度耐磨损的盾构机切割刀片 |
CN114662180A (zh) * | 2021-10-08 | 2022-06-24 | 温州大学 | 基于环间弱连接效应的盾构隧道纵向力学响应计算方法 |
CN114707382A (zh) * | 2022-04-07 | 2022-07-05 | 合肥工业大学 | 一种含弹性衬垫管片接头的抗弯刚度计算方法 |
CN115344935A (zh) * | 2022-09-14 | 2022-11-15 | 浙江大学 | 一种考虑浆液凝固特性的盾构隧道施工期上浮分析计算方法 |
CN115355056A (zh) * | 2022-07-22 | 2022-11-18 | 温州大学 | 市域铁路越江隧道复合衬砌结构受力性能自动化监测系统 |
US11556680B1 (en) * | 2022-01-27 | 2023-01-17 | Qingdao university of technology | Test method based on test system for five-dimensional space effects of ground surface settlement caused by dual-mode shield construction |
CN116305455A (zh) * | 2023-03-13 | 2023-06-23 | 广州大学 | 一种基于OpenSees的盾构隧道-土体动力相互作用模拟方法 |
CN116305424A (zh) * | 2023-02-17 | 2023-06-23 | 浙江大学 | 一种基于状态空间法的考虑错缝拼装盾构隧道开口响应计算方法 |
CN116542109A (zh) * | 2023-05-21 | 2023-08-04 | 南京工业大学 | 一种长大盾构隧道纵向地震反应动力子结构分析方法 |
CN116537820A (zh) * | 2023-04-07 | 2023-08-04 | 湖南大学 | 一种高耐久性叠合式盾构隧道管片结构及其设计方法 |
CN116776435A (zh) * | 2023-06-20 | 2023-09-19 | 同济大学 | 一种盾构隧道纵向多级损伤指标阈值确定方法 |
CN116822267A (zh) * | 2023-03-07 | 2023-09-29 | 北京工业大学 | 一种基于OPenSeesMP考虑局部非连续域的高效并行计算模拟方法 |
WO2023197829A1 (zh) * | 2022-04-14 | 2023-10-19 | 华南理工大学 | 可更换耗能钢框架节点的两阶段抗震设计方法 |
CN116976108A (zh) * | 2023-07-31 | 2023-10-31 | 同济大学 | 一种盾构隧道纵缝接头抗弯刚度反演计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130289953A1 (en) * | 2012-01-24 | 2013-10-31 | The University Of Akron | Self-optimizing, inverse analysis method for parameter identification of nonlinear material constitutive models |
-
2023
- 2023-11-09 CN CN202311483822.5A patent/CN117216862B/zh active Active
Patent Citations (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20020021390A (ko) * | 2002-01-04 | 2002-03-20 | 김용성 | 다층지반의 지진응답해석을 위한 동적 점탄-점소성구성모델 및 이를 이용한 해석방법 |
JP2005076339A (ja) * | 2003-09-02 | 2005-03-24 | Shimizu Corp | 鋼コンクリート合成構造の構造解析方法及び同方法を用いて設計された合成構造。 |
WO2009116048A2 (en) * | 2008-03-20 | 2009-09-24 | Technion Research & Development Foundation Ltd. | A method for cosserat point element (cpe) modeling of nonlinear elastic materials |
JP2009301149A (ja) * | 2008-06-10 | 2009-12-24 | Kyushu Univ | 非線形の有限要素法による構造解析方法、及びプログラム、記録媒体、シミュレーション装置 |
CN101586463A (zh) * | 2009-06-18 | 2009-11-25 | 苏州市交通设计研究院有限责任公司 | 一种盾构隧道衬砌设计方法 |
CN107122536A (zh) * | 2017-04-20 | 2017-09-01 | 济南轨道交通集团有限公司 | 一种区间隧道抗震数值模拟方法 |
WO2019140794A1 (zh) * | 2018-01-16 | 2019-07-25 | 天津大学 | 一种基于量纲分析的盾构掘进总推力计算方法 |
CN111597617A (zh) * | 2020-05-18 | 2020-08-28 | 中铁第六勘察设计院集团有限公司 | 一种盾构隧道疲劳寿命预测方法 |
CN111535830A (zh) * | 2020-07-09 | 2020-08-14 | 湖南大学 | 盾构隧道管片超高性能混凝土加固结构及其设计方法 |
CN111950184A (zh) * | 2020-08-21 | 2020-11-17 | 上海市隧道工程轨道交通设计研究院 | 一种多环盾构衬砌结构的三维仿真计算方法 |
CN112632837A (zh) * | 2020-12-31 | 2021-04-09 | 华中科技大学 | 一种地下结构纵向抗震数值确定方法 |
CN112329122A (zh) * | 2021-01-07 | 2021-02-05 | 湖南大学 | 一种旁侧基坑开挖引起的盾构隧道横向变形和内力确定方法 |
CN113221399A (zh) * | 2021-04-13 | 2021-08-06 | 北京工业大学 | 一种铁路隧道装配式衬砌构件划分方案决策方法 |
CN113283142A (zh) * | 2021-05-31 | 2021-08-20 | 华中科技大学 | 一种基于解析解法分析盾构下穿对既有隧道影响的方法 |
CN113360985A (zh) * | 2021-06-10 | 2021-09-07 | 安徽省建筑科学研究设计院 | 一种新建隧道上穿引起既有盾构隧道纵向变形的高精度预测方法 |
CN113420482A (zh) * | 2021-06-24 | 2021-09-21 | 北京安捷工程咨询有限公司 | 一种基于结构内力监测值的管片荷载正交数值反演方法 |
CN114033405A (zh) * | 2021-09-14 | 2022-02-11 | 山东交通职业学院 | 一种高强度耐磨损的盾构机切割刀片 |
CN114662180A (zh) * | 2021-10-08 | 2022-06-24 | 温州大学 | 基于环间弱连接效应的盾构隧道纵向力学响应计算方法 |
US11556680B1 (en) * | 2022-01-27 | 2023-01-17 | Qingdao university of technology | Test method based on test system for five-dimensional space effects of ground surface settlement caused by dual-mode shield construction |
CN114707382A (zh) * | 2022-04-07 | 2022-07-05 | 合肥工业大学 | 一种含弹性衬垫管片接头的抗弯刚度计算方法 |
WO2023197829A1 (zh) * | 2022-04-14 | 2023-10-19 | 华南理工大学 | 可更换耗能钢框架节点的两阶段抗震设计方法 |
CN115355056A (zh) * | 2022-07-22 | 2022-11-18 | 温州大学 | 市域铁路越江隧道复合衬砌结构受力性能自动化监测系统 |
CN115344935A (zh) * | 2022-09-14 | 2022-11-15 | 浙江大学 | 一种考虑浆液凝固特性的盾构隧道施工期上浮分析计算方法 |
CN116305424A (zh) * | 2023-02-17 | 2023-06-23 | 浙江大学 | 一种基于状态空间法的考虑错缝拼装盾构隧道开口响应计算方法 |
CN116822267A (zh) * | 2023-03-07 | 2023-09-29 | 北京工业大学 | 一种基于OPenSeesMP考虑局部非连续域的高效并行计算模拟方法 |
CN116305455A (zh) * | 2023-03-13 | 2023-06-23 | 广州大学 | 一种基于OpenSees的盾构隧道-土体动力相互作用模拟方法 |
CN116537820A (zh) * | 2023-04-07 | 2023-08-04 | 湖南大学 | 一种高耐久性叠合式盾构隧道管片结构及其设计方法 |
CN116542109A (zh) * | 2023-05-21 | 2023-08-04 | 南京工业大学 | 一种长大盾构隧道纵向地震反应动力子结构分析方法 |
CN116776435A (zh) * | 2023-06-20 | 2023-09-19 | 同济大学 | 一种盾构隧道纵向多级损伤指标阈值确定方法 |
CN116976108A (zh) * | 2023-07-31 | 2023-10-31 | 同济大学 | 一种盾构隧道纵缝接头抗弯刚度反演计算方法 |
Non-Patent Citations (15)
Title |
---|
Experimental study on the vibration damping of two parallel shield tunnels connected by an assembled transverse passage;Junchen Zhang, Qixiang Yan, Minghui Sun, Binjia Li, Wenyu Chen, Hang Chen;《Tunnelling and Underground Space Technology》;第107卷;103659 * |
剪切波作用下盾构隧道地震效应的拟静力分析方法研究;晏启祥;张煜;王春艳;耿萍;;工程力学(第05期);全文 * |
地震作用下盾构隧道衬砌结构动力响应分析;汪亦显;单生彪;任伟新;靖洪文;曹平;;土木工程学报(第S2期);全文 * |
地震荷载作用下大坂盾构隧道管片配筋校核;戴永浩;陈卫忠;谭贤君;黄胜;;岩土力学(第S2期);全文 * |
基于三维非连续接触模型的管片接头静动力特性分析;艾辉军;彭立敏;施成华;;岩土工程学报(第11期);全文 * |
基于层状土推覆方法的盾构隧道抗震性能分析;张业勤;邹炎;刘华北;;防灾减灾工程学报(第06期);全文 * |
岩石初始宏观变形模拟及微裂纹闭合应力确定方法;张超;曹文贵;徐赞;贺敏;;岩土力学(第04期);全文 * |
模拟岩石应变软化变形全过程的统计损伤本构模型研究;汪辉平;曹文贵;王江营;张超;;水文地质工程地质(第04期);全文 * |
泥炭质软土地区盾构隧道抗震性能研究;吕广;刘远才;许纯梅;包琪;;现代城市轨道交通(第06期);全文 * |
深部岩石统计损伤本构模型研究;曹文贵;戴笠;张超;;水文地质工程地质(第04期);全文 * |
盾构管片UHPC 加固技术及力学性能分析;陈仁朋,鲁 立,张 阳,吴怀娜;《工程力学》;第36卷(第11期);第41-50页 * |
盾构隧道地震响应及减震措施研究;孔戈;周健;;现代隧道技术(第S1期);全文 * |
盾构隧道地震响应的三维数值模拟方法及应用;丁峻宏;金先龙;郭毅之;曹伟飚;李根国;;岩石力学与工程学报(第07期);全文 * |
盾构隧道结构计算分析方法研究;何川;张景;封坤;;中国公路学报(第08期);全文 * |
考虑软硬物质双变形特征的脆性岩石损伤本构模型研究;张超;雷勇;曹文贵;;应用力学学报(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN117216862A (zh) | 2023-12-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Abedini et al. | Comparison of ALE, LBE and pressure time history methods to evaluate extreme loading effects in RC column | |
Esmaeily et al. | Behavior of reinforced concrete columns under variable axial loads: analysis | |
An et al. | Numerical simulation of size effect in shear strength of RC beams | |
Padgett et al. | Experimental response modification of a four‐span bridge retrofit with shape memory alloys | |
Abdelnaby et al. | Numerical modeling and analysis of RC frames subjected to multiple earthquakes | |
Zivaljic et al. | A combined finite-discrete element model for RC structures under dynamic loading | |
Zong et al. | Collapse failure of prestressed concrete continuous rigid-frame bridge under strong earthquake excitation: Testing and simulation | |
Tran et al. | Blast behaviour of precast segmental vs monolithic concrete beams prestressed with unbonded tendons: A numerical investigation | |
Mahmoudi Moazam et al. | Three-dimensional modelling for seismic assessment of plain concrete arch bridges | |
Sayed-Ahmed et al. | Nonlinear finite-element model of hollow masonry | |
Zhang et al. | Transverse seismic response of diamond-shaped pylon in cable-stayed bridge: Experiment and analysis | |
Stablon et al. | Influence of building process on stiffness: numerical analysis of a masonry vault including mortar joint shrinkage and crack re-closure effect | |
Anil et al. | Experimental and numerical study of rc columns under lateral low-velocity impact load | |
Scattarreggia et al. | A new Distinct Element meso‐model for simulating the rocking‐dominated seismic response of RC columns | |
Zhao et al. | An effective framework for performance evaluation of reinforced concrete beams under impact loadings | |
Gerolymos et al. | Macroelement modeling of piles in cohesive soil subjected to combined lateral and axial loading | |
Blasi et al. | In-plane and out-of-plane model for retrofitted infill walls in reinforced concrete framed buildings | |
CN117216862B (zh) | 一种基于纤维梁单元的三维盾构隧道动力分析模型 | |
Haido et al. | Simulation of dynamic response for steel fibrous concrete members using new material modeling | |
Mazars et al. | Numerical modelling for earthquake engineering: the case of lightly RC structural walls | |
Živaljić et al. | Numerical simulation of reinforced concrete structures under impact loading | |
Zhang et al. | Mechanical response of tunnels crossing active fault zones under normal or reverse faulting: A refined nonlinear analytical approximate solution | |
Wang et al. | Development of dimensionless PI diagram for curved SCS sandwich shell subjected to uniformly distributed blast pressure | |
Cui et al. | Numerical study on parameter impact on fundamental frequencies and dynamic characteristics of pre-stressed concrete beams | |
Yousefi et al. | A multi-fiber approach with directional stiffness matrix in reinforced concrete structures |
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 |