CN106202649A - 考虑软基固结和混凝土徐变的闸首施工仿真方法 - Google Patents

考虑软基固结和混凝土徐变的闸首施工仿真方法 Download PDF

Info

Publication number
CN106202649A
CN106202649A CN201610499037.2A CN201610499037A CN106202649A CN 106202649 A CN106202649 A CN 106202649A CN 201610499037 A CN201610499037 A CN 201610499037A CN 106202649 A CN106202649 A CN 106202649A
Authority
CN
China
Prior art keywords
concrete
lock head
creep
soil
ground
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
CN201610499037.2A
Other languages
English (en)
Other versions
CN106202649B (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN201610499037.2A priority Critical patent/CN106202649B/zh
Publication of CN106202649A publication Critical patent/CN106202649A/zh
Application granted granted Critical
Publication of CN106202649B publication Critical patent/CN106202649B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING; 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

Abstract

本发明公开了一种考虑软基固结和混凝土徐变的闸首施工仿真方法,步骤为:1、确定混凝土和软土地基的热学参数、混凝土的徐变计算参数、地基的固结参数、地基的弹塑性力学参数;2、建立闸首结构的有限元网格;3、确定计算的时间步骤,混凝土的初始温度,混凝土保温措施的等效放热系数;4、根据热学参数、时程分析步骤、初始条件和气温预测计算温度场;5、计算开挖后的软土地基的孔压和地应力分布;6、施加温度荷载、地基初始孔压、地基地应力和需要考虑的其他荷载进行应力计算。本发明方法综合考虑了混凝土的徐变特性和软土地基的弹塑性固结特性,能较完善的预测施工期闸首结构的应力。

Description

考虑软基固结和混凝土徐变的闸首施工仿真方法
技术领域
[0001] 本发明涉及水利工程混凝土结构设计中软基上闸首结构施工期综合考虑软基固 结和混凝土徐变的应力分析方法。
背景技术
[0002] 内河水运具有运量大、成本低、效益大、对环境污染小等特点,对促进流域经济发 展发挥了重要作用。船闸是水运工程中的主要通航建筑物。
[0003] 闸首混凝土方量虽然相对于其他水工大坝等建筑物而言是较小的,但仍需要作为 大体积混凝土结构进行考虑。尤其是底板,一般厚度达到2~3米,且大部分一次浇筑完成; 边墩方量相对较小,但由于设置了廊道,空箱等薄壁结构,使得闸首结构错综复杂。
[0004] 现有闸首计算方法主要以规范计算方法为主,即将上部结构、基础和地基各自独 立分开,将每一部分作为独立的结构利用材料力学的方法进行计算。
[0005] 现有研究分析,船闸的大体积混凝土结构开裂的主要是由于结构浇筑后混凝土的 温度应力造成的。
[0006] 软土地基的渗透性低,具有明显的固结特性,造成地基的沉降随时间是变化的,且 现有土力学已指出土体具有弹塑性特性,传统的弹性理论已不适用。
[0007] 地基的较大沉降和不均匀沉降会对上部结构的安全性和稳定性有着严重的影响, 尤其对于闸首底板此类大体积混凝土板结构。
发明内容
[0008] 本发明的目的在于针对现有软基上闸首结构设计的不足,综合考虑上部结构-基 础-地基系统整体性、大体积混凝土的温度场和徐变特性、软土地基弹塑性固结沉降特性, 提供一种多物理场整体分析的软基上闸首结构施工期的应力分析方法。
[0009] 本发明所述的基于有限单元法的软基上闸首结构施工期结构应力的分析方法,包 括如下步骤:
[0010] -种考虑软基固结和混凝土徐变的闸首施工仿真方法,包括以下步骤:
[0011] (1)若在前期设计过程中已进行足够的实验测试,则根据前期设计报告,确定船闸 闸首结构的混凝土和软土地基在计算过程中所需要的热学和力学计算参数;若在前期设计 过程中未进行足够的实验测试,则根据当地工程的参数经验值进行预计算,并通过施工前 期的部分监测数据进行反演校正;
[0012] (2)根据闸首结构和地基开挖的设计方案,使用有限元前处理软件建立闸首结构 与地基的三维几何模型,并进行网格划分,生成可用于有限元计算的有限元网格;
[0013] (3)根据闸首的施工方案,确定计算分析的时间步骤,同时估算混凝土饶筑的初始 温度和计算期间内的气温,同时根据混凝土表面的保温措施,计算闸首结构与空气接触面 部分的等效的对流换热系数;
[0014] (4)跟据上述所得的初始温度、对流边界条件、水化热和时程分析步骤,计算混凝 土和地基的温度场变化;
[0015] (5)根据地基开挖的方案,计算开挖后的地基的初始孔压和地应力的分布;
[0016] (6)将已计算得到的温度场作为温度荷载施加到闸首结构和地基上,并在地基中 应用步骤5中得到的初始孔压和地应力,并施加其他荷载(如重力),混凝土采用徐变模型, 软土地基采用弹塑性固结模型,进行多物理场计算。
[0017] 进一步的,步骤1中所述的计算参数包括:混凝土的热传导系数、混凝土的比热、混 凝土的密度、混凝土的水化热模型参数、混凝土的对流换热系数、混凝土的热膨胀系数、混 凝土的弹性模量、混凝土的泊松比、混凝土的徐变模型参数、土体的热传导系数、土体的比 热、土体的密度、土的对流换热系数、土的热膨胀系数、土体的渗透系数、土体的弹性模量和 土体的弹塑性模型参数。
[0018] 进一步的,步骤2中所述的闸首结构包括底板、廊道、空箱、切角和门库。
[0019] 进一步的,步骤2中所述的地基根据真实的工程地质勘测进行合理的分层。
[0020] 进一步的,步骤3中所述的气温,如果对于短期计算预测,则使用天气预报数据;如 果对于长期计算预测,则使用工程经验模型。
[0021] 进一步的,所述的混凝土的水化热模型参数,其中水化热模型采用双指数模型为:
Figure CN106202649AD00051
[0023] 其中,Wc为混凝土中水泥用量,Q~为水泥最终水化放热量,a和b为混凝土水化热实 验参数。
[0024] 进一步的,所述的混凝土的弹性模量和徐变模型参数,其中弹性模量随时间的变 化模型为:
Figure CN106202649AD00052
[0026]其中,为最终弹性模量,a和b为混凝土弹性模量的实验参数;
[0027]徐变模型为:
[0028] C( t,τ) =Ca( 1+9 · 2τ-0.45) [ l-e-0.3(t-τ) ]+Cb( 1+1 · 7τ-0.45) [ l-e-0.005(t-τ)]
[0029] 其中,CjPCb为实验参数,t为时间。
[0030]进一步的,所述的土体弹塑性模型参数,土体的弹塑性模型采用修正剑桥模型,其 屈服函数为:
Figure CN106202649AD00053
[0032] 其中,Μ为应力比,p为平均应力,q为广义剪切应力,ρο为预固结压力;在计算时所 需要的参数有:正常固结线斜率λ,回弹线斜率κ,应力比M,预固结压力po。
[0033] 本发明一种考虑软基固结和混凝土徐变的闸首施工仿真方法与现有技术相比,其 优点在于:
[0034]第一、将上部结构-基础-地基系统进行整理考虑,不再将它们独立分开计算;
[0035]第二、考虑了对早期混凝土影响较大的徐变特性;
[0036]第三、考虑了软土地基的弹塑性固结特性。
[0037]以上三点使得闸首施工期的计算较为完善的考虑到各种影响因素,使得计算结果 更为合理准确。
附图说明
[0038] 图1为本发明分析方法流程图;
[0039] 图2为某闸首和开挖后地基的有限元模型;
[0040] 图3为温度监测点和数据提取点;
[0041] 图4为日平均气温曲线;
[0042] 图5为底板处温度监测值与计算值对比;
[0043]图6为开挖后孔隙水压力平衡值;
[0044]图7为开挖后铅直应力平衡值;
[0045] 图8为弹性和弹塑性固结计算的S2点和M2点垂直位移的变化对比图;
[0046] 图9为弹性和弹塑性固结计算的S1、S2、M1和M2点处最大主应力变化对比图;
[0047] 图10为边底板浇筑完毕时地基孔隙水压力分布图。
具体实施方式
[0048] 以下结合一个实例对本发明作进一步详细说明。
[0049] 本实例中混凝土水化热采取公式-(丨十,.热学计算参数参见表1。
[0050] 表1热学参数
Figure CN106202649AD00061
[0052] 混凝土的弹性模量公式为£'(r)-£, j ,其中E~=34.25X 106kPa,a = 0.28,b =0.52 ;混凝土的泊松比为0.167 ;混凝土徐变度公式为C( t,τ) = Ca( 1+9.2τ-Q·45) [ l-e -0·3(t-τ) ] +Cb (1+1 · 7 τ-。.45)[卜e-。.005 (τ)],其中 Ca = 6 · 39 X 10-6/MPa,Cb = 1 · 444 X 10-5/MPa。
[0053] 本实例中软土地基根据工程地质勘测共分为9层,使用修正剑桥固结模型,土的泊 松比为0.3,其他计算所需要的参数值参见表2。
[0054] 表2 土体弹塑性固结参数
Figure CN106202649AD00071
[0057]本实例中的有切角闸首结构在横河向宽为26.9m,在顺河向长为28.6m,在铅直向 上的高度为11.9m(其中底板厚2.6m),闸首结构和开挖后地基的有限元模型如附图2所示, 其中的底板上的温度监测点(T1、T2和T3)和结果数据提取点(M1、M2、S1和S2)的示意图如附 图3所示。
[0058]施工方案如表3所示,施工期间的气温如附图4所示。
[0059] 表3施工步骤及时长
Figure CN106202649AD00072
[0061] 根据上述参数及时程步骤计算温度场,温度监测点T1、T2和T3的监测结果与计算 结果的对比如附图5所示。
[0062] 根据基坑开挖方案:分两级向下开挖,平台上坡比为1:2,中间平台宽度为5m,平台 下坡比为1:2.5。进行开挖后的孔隙水压力和地应力的平衡。附图6为开挖后地基的孔隙水 压力平衡值,附图7为开挖后地基的铅直应力平衡值。
[0063]将平衡后的孔隙水压力和地应力作为地基弹塑性固结计算的初始条件,将温度场 结果作为闸首结构和地基的温度荷载,同时施加结构和地基的重力荷载。进行位移场和应 力场的计算。
[0064]作为对比,同时进行了将地基作为弹性模型的计算。附图8为一段时间内的M2和S2 点处垂直位移弹性计算和弹塑性固结结算的结果对比,可见弹塑性固结计算的结果具有明 显的时间效应,并考虑到了施工过程中孔隙水压力变化的影响。
[0065]附图9为弹性和弹塑性固结计算的S2、S2、M1和M2点处最大主应力变化对比图。弹 性计算的值总体上偏大,超过规范中的强度值,通过弹塑性固结计算得到的结果更为合理, 更符合现场观测得到的结果。
[0066] 使用弹塑性固结地基计算可同步得到软基中孔隙水压力的分布,从而反映孔隙水 压力对上部结构的影响,如附图10为边底板浇筑完毕时地基的孔隙水压力分布图,由于软 土的较低渗透性,孔隙水压力在边底板底部增大,从而影响边底板底部的位移和地基的沉 降,随着孔隙水压力的消散,地基沉降会逐渐增大。
[0067] 本发明不局限于上述实施方式,对于本技术领域的普通技术人员来说,在不脱离 本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护 范围之内。本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (8)

1 · 一种考虑软基固结和混凝土徐变的闸首施工仿真方法,其特征在于,包括以下步骤: (1) 若在前期设计过程中已进行足够的实验测试,则根据前期设计报告,确定船阐闹首 结构的混凝土和软土地基在计算过程中所需要的热学和力学计算参数;若在前期设计过程 中未进行足够的实验测试,则根据当地工程的参数经验值进行预计算,并通过施工前期的 部分监测数据进行反演校正; 。 (2) 根据闸首结构和地基开挖的设计方案,使用有限元前处理软件建立闸首结构与地 基的二维几何模型,并进行网格划分,生成可用于有限元计算的有限元网格; (3) 根据闸首的施工方案,确定计算分析的时间步骤,同时估算混凝土浇筑的初始温度 和计算期间内的气温,同时根据混凝土表面的保温措施,计算闸首结构与空气接触面部分 的等效的对流换热系数; (4) 跟据上述所得的初始温度、对流边界条件、水化热和时程分析步骤,计算混凝土和 地基的温度场变化; (5) 根据地基开挖的方案,计算开挖后的地基的初始孔压和地应力的分布; (6) 将已计算得到的温度场作为温度荷载施加到闸首结构和地基上,并在地基中应用 步骤5中得到的初始孔压和地应力,并施加其他荷载,如重力;混凝土采用徐变模型,软土地 基采用弹塑性固结模型,进行多物理场计算。 2 ·如权利要求1所述考虑软基固结和混凝土徐变的闹首施工仿真方法,其特征在于,步 骤1中所述的计算参数包括:混凝土的热传导系数、混凝土的比热、混凝土的密度、混凝土的 水化热模型参数、混凝土的对流换热系数、混凝土的热膨胀系数、混凝土的弹性模量、混凝 土的泊松比、混凝土的徐变模型参数、土体的热传导系数、土体的比热、土体的密度、土的对 流换热系数、土的热膨胀系数、土体的渗透系数、土体的弹性模量和土体的弹塑性模型参 数。
3. 如权利要求1所述考虑软基固结和混凝土徐变的闸首施工仿真方法,其特征在于,步 骤2中所述的闸首结构包括底板、廊道、空箱、切角和门库。
4. 如权利要求1所述考虑软基固结和混凝土徐变的闸首施工仿真方法,其特征在于,步 骤2中所述的地基根据真实的工程地质勘测进行合理的分层。
5. 如权利要求1所述考虑软基固结和混凝土徐变的阐首施工仿真方法,其特征在于,步 骤3中所述的气温,如果对于短期计算预测,则使用天气预报数据;如果对于长期计算预测, 则使用工程经验模型。
6. 如权利要求2所述考虑软基固结和混凝土徐变的闸首施工仿真方法,其特征在于,所 述的混凝土的水化热模型参数,其中水化热模型采用双指数模型为: 其中,Wc为混凝土中水泥用量,Qco为水泥最终水化放热量,a和b为混凝土水化热实验参 数。
7. 如权利要求2所述考虑软基固结和混凝土徐变的阐首施工仿真方法,其特征在于,所 述的混凝土的弹性模量和徐变模型参数,其中弹性模量随时间的变化模型为: 茗(輕 其中,Eoo为最终弹性模量,a和b为混凝土弹性模量的实验参数; 徐变模型为: Ο(ί,τ)=〇3(1+9.2τ^0·45)[1-θ_0· 3(^τ)]+α(1+1.7τ^0·45)[ΐ- θ·°·005(ί·τ)] 其中,CdPCb为实验参数,t为时间。
8.如权利要求2所述考虑软基固结和混凝土徐变的闸首施工仿真方法,其特征在于,所 述的土体弹塑性模型参数,土体的弹塑性模型采用修正剑桥模型,其屈服函数为:
Figure CN106202649AC00031
其中,Μ为应力比,p为平均应力,q为广义剪切应力,pQ为预固结压力;在计算时所需要的 参数有:正常固结线斜率λ,回弹线斜率κ,应力,预固结压力p〇。
CN201610499037.2A 2016-06-29 2016-06-29 考虑软基固结和混凝土徐变的闸首施工仿真方法 Active CN106202649B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610499037.2A CN106202649B (zh) 2016-06-29 2016-06-29 考虑软基固结和混凝土徐变的闸首施工仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610499037.2A CN106202649B (zh) 2016-06-29 2016-06-29 考虑软基固结和混凝土徐变的闸首施工仿真方法

Publications (2)

Publication Number Publication Date
CN106202649A true CN106202649A (zh) 2016-12-07
CN106202649B CN106202649B (zh) 2018-09-25

Family

ID=57463299

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610499037.2A Active CN106202649B (zh) 2016-06-29 2016-06-29 考虑软基固结和混凝土徐变的闸首施工仿真方法

Country Status (1)

Country Link
CN (1) CN106202649B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106919786A (zh) * 2017-02-07 2017-07-04 三峡大学 一种混凝土分数阶徐变模型
CN107633146A (zh) * 2017-09-28 2018-01-26 国家电网公司 一种不同尺度模型间地应力计算结果高精度转换的方法
CN109460582A (zh) * 2018-10-16 2019-03-12 河海大学 基于粒子群算法的顺水流向变高程廊道的人字门闸首结构优化设计方法
CN109726485A (zh) * 2018-12-30 2019-05-07 中国水利水电科学研究院 混凝土坝点真实应力可靠度预测分析方法及装置
CN111504252A (zh) * 2020-04-23 2020-08-07 长江水利委员会长江科学院 一种长距离隧洞膨胀性围岩变形超前预测预报方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101260664A (zh) * 2008-04-22 2008-09-10 中船第九设计研究院工程有限公司 一种软土地基大型钢浮箱闸首的施工方法
JP2009211452A (ja) * 2008-03-05 2009-09-17 Tokyo Electric Power Co Inc:The コンクリート柱強度計算装置
EP2107481A1 (en) * 2006-11-29 2009-10-07 Airbus España, S.L. Thermal simulation methods and systems for analysing fire in objects
CN104657538A (zh) * 2015-01-04 2015-05-27 河海大学 一种多节点航道调度的仿真方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2107481A1 (en) * 2006-11-29 2009-10-07 Airbus España, S.L. Thermal simulation methods and systems for analysing fire in objects
JP2009211452A (ja) * 2008-03-05 2009-09-17 Tokyo Electric Power Co Inc:The コンクリート柱強度計算装置
CN101260664A (zh) * 2008-04-22 2008-09-10 中船第九设计研究院工程有限公司 一种软土地基大型钢浮箱闸首的施工方法
CN104657538A (zh) * 2015-01-04 2015-05-27 河海大学 一种多节点航道调度的仿真方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
苏超 等: "三峡永久船闸衬砌墙结构及结构锚杆受力仿真分析", 《河海大学学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106919786A (zh) * 2017-02-07 2017-07-04 三峡大学 一种混凝土分数阶徐变模型
CN106919786B (zh) * 2017-02-07 2019-04-09 三峡大学 一种混凝土分数阶徐变模型
CN107633146A (zh) * 2017-09-28 2018-01-26 国家电网公司 一种不同尺度模型间地应力计算结果高精度转换的方法
CN107633146B (zh) * 2017-09-28 2021-03-30 国家电网公司 一种不同尺度模型间地应力计算结果高精度转换的方法
CN109460582A (zh) * 2018-10-16 2019-03-12 河海大学 基于粒子群算法的顺水流向变高程廊道的人字门闸首结构优化设计方法
CN109726485A (zh) * 2018-12-30 2019-05-07 中国水利水电科学研究院 混凝土坝点真实应力可靠度预测分析方法及装置
CN111504252A (zh) * 2020-04-23 2020-08-07 长江水利委员会长江科学院 一种长距离隧洞膨胀性围岩变形超前预测预报方法
CN111504252B (zh) * 2020-04-23 2021-07-02 长江水利委员会长江科学院 一种长距离隧洞膨胀性围岩变形超前预测预报方法

Also Published As

Publication number Publication date
CN106202649B (zh) 2018-09-25

Similar Documents

Publication Publication Date Title
CN106202649B (zh) 考虑软基固结和混凝土徐变的闸首施工仿真方法
Bofang Thermal stresses and temperature control of mass concrete
CN104480962B (zh) 一种有限填土挡墙的土压力分布计算方法
Hahn et al. Breach morphology observations of embankment overtopping tests
CN106649931A (zh) 一种矩形顶管隧道施工工后地表沉降计算方法
Xu et al. Frost heave of irrigation canals in seasonal frozen regions
CN110263355A (zh) 复杂地基特高拱坝整体抗裂设计方法
CN103215914B (zh) 堤身沉降控制方法
Liferov et al. 4.5 3D Finite element analysis of laboratory punch tests on ice rubble
CN104196024B (zh) 基于锚索预应力损失与岩土体蠕变耦合的预应力损失计算方法
Ahayan A constitutive Model for natural Clays: From Laboratory Testing to Modelling of Offshore Monopiles
Zheng et al. Cause and Stability Analysis of Cracks in Concrete Slab of Rockfill Dam under High Temperature Difference Condition
Corns et al. Gravity dam design and analysis
Lian et al. Surficial stability analysis of soil slope under seepage based on a novel failure mode
Simonsen et al. Prediction of pavement response during freezing and thawing using finite element approach
Zeidan Design and analysis of concrete gravity dams
CN106702999A (zh) 道路复合地基静载荷试验预判工后沉降的方法
Herdy et al. Preliminary investigation of bearing capacity of layered soils by centrifugal modeling
Wang et al. Arch dam and skewback and dam abutment rock mass deformation under seepage field and stress field coupling
Wieland et al. Earthquake aspects of roller compacted concrete and concrete-face rockfill dams
Li et al. Bearing Properties of Large-Diameter Embedded Rock-Socketed Pile with Steel Tube in Frame Wharf
Du et al. Numerical Simulation of Thermal Coupling Process of Foundation Pit under Freezing and Thawing
Lei et al. Analysis of the Causes of Cracks in the Bottom Floor of the Underground Garage of the Hefei Government Affairs Center by using 3D Finite-Element Analysis.
CN109853598A (zh) 一种兼具防滑坡和生产道路的微型桩刚性连接结构及施工方法
Zhang et al. Coupling Analysis of Fracture Seepage and Rock Deformation in DDA Method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant