地铁隧道水平冻结法施工冻结壁温度场影响参数分析精Word下载.docx
- 文档编号:21828972
- 上传时间:2023-02-01
- 格式:DOCX
- 页数:12
- 大小:511.43KB
地铁隧道水平冻结法施工冻结壁温度场影响参数分析精Word下载.docx
《地铁隧道水平冻结法施工冻结壁温度场影响参数分析精Word下载.docx》由会员分享,可在线阅读,更多相关《地铁隧道水平冻结法施工冻结壁温度场影响参数分析精Word下载.docx(12页珍藏版)》请在冰豆网上搜索。
A
Analysisofinfluencingparametersoftemperaturefieldwithinfreezing
wallofsubwaytunnelconstructionwithhorizontalgroundfreezing
BIGuiquan1,2,CHENGXingyan1,SHILei2,MAXijin1
(1.CollegeofEnergyandPowerEngineering,LanzhouUniv.ofTech.,Lanzhou730050,China;
2.StateKeyLaboratoryofFrozenSoil
Engineering,CAREERI,CAS,LanzhouGansu730000,China
Abstract:
Takingthemostcommonsiltyclayusedinsubwaytunnelconstructionwithhorizontalfreezingmethodasanexample,theeffectsofbrinetemperature,diameteroffrozenpipes,intervaloffrozenpipes,andmoisturecontentofsoillayeronthetemperaturefieldwithinfrozenwallforthefreeingconstructionweresystematicallyinvestigatedwithlargescalenumericalsoftwareANSYS.Theformulasofsensitivityofeveryabovementionedfactortothethicknessofthefrozenwallandfreezingtimewerepresented.Byanalyzingandcontrastingthesensitivitytoeveryfactor,thepatternofinfluenceofeveryfactoronthesensitivityoffrozenwallthicknessresultshowedthatindecreasingorder,thebrinetemperature,intervaloffrozenpipes,diameteroffrozenpipes,andmoisturecontentofsoillayerinfluencedthesensitivity.
Keywords:
frozenwall;
temperaturefield;
diameteroffrozenpipes;
intervaloffrozenpipes;
brinetemperature;
moisturecontentofsoillayer
地层人工冻结技术是利用人工制冷的方法,降低土体的温度使含水地层形成冻结体,以达到维护开挖面周围土体稳定、抵抗周围水土压力、防止地下水侵入的目的.近年来,人工冻结技术被越来越多的应用到复杂条件下的矿山、港口工程及水工工程中,特别是沿海城市的地铁隧道及其联络通道施工中.地铁隧道目前广泛采用的水平冻结法施工形成的冻结壁是抵御水土压力的临时地下结构物[1],要确定
收稿日期:
20080925
基金项目:
科学院知识创新工程重大项目(KZCX1SW04
作者简介:
毕贵权(1979,男,吉林九台人,讲师.冻结壁承载力并进行冻结壁厚度和其它参数设计,必须要对温度场进行分析.冻结壁温度场是一个相变、移动边界和有内热源、边界条件复杂的不稳定导热问题,其求解极其困难,至今还没有完善的理论解.通常采用的研究方法有半解析法、数值模拟法、物理模型试验法等.由于物理模型试验法具有费用高、耗时长等缺点,因此,采用计算机进行数值模拟具有非常好的经济价值和使用价值.数值法则在近年来又有了新的进展,尤其表现在相变的处理技术上.目前相变处理主要有显热容法和焓法[2],国内外许多的研究人员对该问题采用数值方法进行了研
第35卷第3期2009年6月
兰州理工大学学报
JournalofLanzhouUniversityofTechnology
Vol.35No.3
Jun.2009
究,并取得了一定的成果.徐士良介绍了利用ANSYS有限元软件对冻结壁温度场分布进行数值模拟的基本过程,包括模型的建立,问题的求解及结果分析等[3].丁德文等,任松保采用有限差分方法对人工冻结温度场进行分析[45].李磊等考虑了冻结管偏斜布置情况,利用有限元方法对上海市复兴东路隧道联络通道进行了三维温度场分析.汪仁和等根据冻结壁温度场物理模型试验的结果,建立非线性多维导热系数反分析的数学模型,采用选择法和有限元法,反分析冻结壁整体导热系数[6]
.汪东波通过模型试验及其有限元方法,对双排管冻结温度场分布规律进行研究;
并通过近似方法推导出双排管冻结壁内外侧厚度以及冻结壁交圈时间解析解[7]
.曹荣斌建立双排冻结管冻结条件下平面冻结壁温度场与空间冻结壁温度场数值模型,利用ANSYS有限元软件数值模拟平面冻结壁温度场的发展过程与分布规律,以及空间状态冻结壁温度场分布规律[8].目前尚未有人采用数值方法系统研究各人为可控制因素(盐水温度、冻结管直径、冻结管间距以及土层含水量对冻结壁温度场的量化影响规律,而这也这是工程设计人员所关心的问题.
针对各因素对冻结壁温度场的影响进行分析研究.在所有的影响因素中,盐水温度、冻结管直径、冻结管间距和土层含水量是主要的影响因素[9],它们的动态变化直接影响着冻结壁温度场的变化.以地铁水平冻结法施工中最常见的粉质粘土为例,应用数值软件ANSYS,采用热焓法考虑冰水相变作用,针对冻结管直径、冻结管间距、盐水温度、土层含水量4大因素对冻结壁温度场的影响进行系统研究,提出各因素对冻结壁厚度和冻结时间的灵敏度公式,并对各个因素影响下冻结壁厚度关于灵敏度进行分析,得出各因素对冻结壁温度场的影响规律.
1基于ANSYS的数值分析
对粉质粘土隧道冻结法施工中冻结壁温度场进行数值试验,试验中各参数采取了工程经验常用值;
在此种土质下对冻结壁温度场发展影响较大的冻结管直径、冻结管间距、盐水温度、土层含水量4个因素进行探讨.主要进行试验前的基本假定、数值的模拟及边界条件、参数的选取、相变的模拟.1.1基本假定
1视未冻土体、冻土体为均质、各向同性材料,沿冻结管水平方向无温度梯度.
2仅考虑土体冻结过程中土骨架和介质水的
热传导以及冰水相变作用.
3冻土的导热系数和比热容不随着温度变化.1.2数值试验方案及边界条件
选用单向地铁隧道为模拟对象,隧道埋深在地面以下20m.采用二维计算模型,x方向取-30m到30m;
y方向取-30m到20m,隧道轴心在原点处.应用大型计算软件ANSYS进行数值试验,选取PLANE77高阶单元.PLANE77是2维4节点热单元PLANE55的高阶版本,每个节点只有1个自由度温度,该8节点单元具有一致的温度形函数,可以较好地适应具有曲线边界的模型,适用于2维稳态或瞬态热分析.为了使计算有足够的精度且节点尽量少,划分单元时以隧道轴线为坐标原点,半径为5m的范围内均用三角形单元,以外用四边形单元.模型两侧为绝热边界,下边界有地热作用,热流密度154J/(hm2,上边界与大气相通,大气与地表的热交换系数为18W/(m2
.采用某地区的多年气温统计资料显示的地温变化曲线[10]
:
T=23.5+
7.8sin28760th+
2.边界条件及有限元网格划分
如图1所示.
图1数值试验模型Fig.1Modelofnumericaltest
1.3相变的模拟
土壤冻结过程为相变导热过程,相变潜热的释放对冻结壁温度场的发展有重要影响.通过ANSYS热分析中的热焓来考虑人工冻结过程中的冰水相变,考虑到相变是在很短的时间内放出或吸收大量的热,其计算如下.
T1=-30H1=0(1T2=-0.5H2=CfT2-T1(2T3=0H3=H2+Q(3T4=30H4=H3+CuT4-T3(4其中,T为温度,H为热焓,Cu、Cf分别为冻土和融土的比热,Q为相变热.
单位体积土中由于水的相态改变所放出和吸收的热量,单位J/m3
可按下式计算
[11]
Q=L
dW-Wu式中:
L为水的结晶或融化潜热,L=334.56kJ/kg;
d为土的干密度;
W为土中水总的质量分数;
Wu为
122
兰州理工大学学报第35卷
冻土中水的质量分数.
对粉质粘土Wu=K(TWp,Wp为塑限水的质量分数,K(T为温度修正系数
[12]
.
1.4数值试验参数
粉质粘土各参数如表1
[12,13]
表1粉质粘土的热学参数
Tab.1Parametersofthecalorificsofsiltyclay
/(kgm-3w/%Cu/(kJm-3-1Cf/(kJm-3-1u/(Jm-1
s-1-1f/(Jm-1
s-1
-1
6001103362.32070.10.631.32900703538.02354.50.751.141200352760.12032.50.861.451500
25
2822.9
2229.0
1.13
1.58
注:
Cu、Cf分别为融土和冻土的比热容,u、f分别为融土和冻土的导热系数.
2数值试验结果及分析
粉质粘土隧道冻结法施工中对冻结壁温度场进行数值试验,试验中各参数采取工程经验常用值.对此种土质影响冻结壁温度场发展较大的4大因素进行探讨.分析各参数对冻结壁厚度和冻结速率的发展规律的影响,以及灵敏度的分析.2.1参数影响
1冻结管直径的影响.
取各个因素的工程经验常用值,当其余参数保持不变时,通过改变冻结管直径数值来进行数值试验,分别为6、8、10、12、14cm5种方案.冻结管间距为0.8m,盐水温度为-26,水的质量分数为35%.因为冻结壁厚度和冻结速率均是时间的函数,所以对冻结160d内的冻结壁温度场进行分析.图2为不同冻结管直径下冻结壁的厚度随冻结时间的的变化图.随着冻结时间的增加冻结壁的厚度均是基本呈线性增大,同一时刻冻结管的直径越大冻结壁越厚.冻结壁的交圈厚度都在1m左右
图2不同冻结管直径下冻结壁厚度随冻结时间的变化Fig.2Variationofthicknessoffrozenwallwithfreezing
timeincaseofdifferentdiametersofthefrozenpiples
冻结壁交圈时间为10d到16d之间,第160d冻结
壁的厚度最小为4.9m,最大为5.73m.随着冻结时间的增长同一时刻不同冻结管直径下的冻结壁厚度的差值变大.
图3为不同冻结管直径下冻结壁的冻结速率随冻结时间的变化图.各个因素下的冻结速率发展趋势是一样的,均是随着时间的增大而减小,在冻结初期冻结速率很大,最大为98mm/d,最小为58mm/d,到90d后趋于平缓变化很小,140d后各直径下的冻结速率随着时间的增大有增大的趋势
图3不同冻结管直径下冻结壁的冻结速率随冻结时间
的变化
Fig.3Variationoffreezingrateoffrozenwallwith
freezingtimeincaseofdifferentdiametersofthefrozenpiples
2冻结管间距的影响.
取工程经验常用值,当其余参数保持不变时,通过改变冻结管间距来进行数值试验,参数为0.6、0.8、1.0、1.2m.冻结管直径为10cm,盐水温度为-26,水的质量分数为35%.
图4为不同冻结管间距下冻结壁的厚度随冻结时间的变化图,可见各条曲线随着时间的增加逐渐增大,冻结管间距越小同一时刻冻结壁就越厚,交圈时间也越早,交圈时间在20d左右.在160d时间距为0.6m的冻结壁厚度最大为5.8m,间距为1.2m的冻结壁的厚度为4.4m.
图5为不同间距下冻结壁的冻结速率随冻结时间的变化图,可见间距越小冻结速率就越大,冻结速率均随着时间的增加递减;
且起始的冻结速率降低的很快,间距为0.6m时的交圈速率最大,在16d时为81.5mm/d.随着时间的增加各冻结管间距下的冻结速率逐渐趋于平缓.
3盐水温度的影响.
取工程经验常用值,当其余参数保持不变时,通过改变盐水温度的数值来进行数值试验,分别为-14、-18、-22、-26、-30、-346种方案.冻结管直径为10cm,冻结管间距为0.8m,水的质量
123第3期毕贵权等:
地铁隧道水平冻结法施工冻结壁温度场影响参数分析
图4不同间距下冻结壁厚度随冻结时间的变化Fig.4Variationofthicknessoffrozenwallwithfreezing
timeincaseofdifferentintervalofthefrozenpi
ples
图5不同间距下冻结壁的冻结速率随冻结时间的变化Fig.5Variationoffreezingratefreezeoffrozenwall
withfreezingtimeincaseofdifferentintervaloffrozenpiples
分数为35%.
图6为不同盐水温度下冻结壁的厚度随冻结时间的变化图,可见随着时间的增加各个盐水温度下的冻结壁厚度是逐渐增大的,刚交圈时冻结壁的厚度相差很小,随着时间的增加各盐水温度下的冻结壁的厚度差别越来越大;
盐水温度越低同一时刻的冻结壁厚度越大,交圈时间越早;
同一时刻盐水温度每降低4,冻结壁的厚度增加的量值是不一样的,以88d为例,温度每降低4冻结壁厚度增加值分别为0.37、0.42、0.35、0.22、0.33
m.
图6不同盐水温度下冻结壁的厚度随冻结时间的变化Fig.6Variationofthicknessoffrozenwallwithfreezing
timeincaseofdifferentbrinetemperature
图7为不同盐水温度下冻结壁的冻结速率随
冻结时间的变化图,可见各个盐水温度下冻结速率均是先快速降低,随着时间增加降低速度慢慢变小,
最终趋于稳定.盐水温度越低交圈冻结速率就越大,达到稳定时的冻结速率也最大;
随着冻结时间的增加各盐水温度下的冻结速率的差值逐渐减小.
图7不同盐水温度下冻结壁的冻结速率随冻结时间的
变化
Fig.7Variationofratefreezingoffrozenwallwith
freezingtimeincaseofdifferentbrinetemperature
4水的质量分数的影响.
取工程经验常用值,当其余参数保持不变时,通
过改变水的质量分数的数值来进行数值试验,分别为25%、35%、70%、110%4种方案,冻结管直径为10cm,冻结管间距为0.8m,盐水温度为-26.
图8是不同水的质量分数下冻结壁的厚度随冻结时间的变化图,可见水的质量分数越大同一时刻冻结壁厚度就越小,水的质量分数增加相同的量时冻结壁厚度急剧减小.以124d为例,水的质量分数从25%增加到35%、70%、110%时冻结壁厚度的减少量分别为0.97、1.56、0.14m;
水的质量分数越大,同一时刻冻结壁厚度的变化也就越小;
由图知水的质量分数为25%时在140d之前冻结壁厚度增加很快,140d之后冻结壁厚度发展很慢;
水的质量分数为70%、110%时冻结壁厚度基本呈线性发展.
图8不同水的质量分数下冻结壁厚度随冻结时间的变化Fig.8Variationofthicknessoffrozenwallwithfreezing
timeincaseofdifferentmoisturecontent
图9是不同水的质量分数下冻结壁的冻结速
124
率随冻结时间的变化图.可见刚交圈时由于受交圈冻结速率的影响冻结速率都很大,随着时间的增加
慢慢趋于平稳;
水的质量分数越小冻结速率越大,水的质量分数25%时冻结速率先减小后趋于稳定,在140d后又继续减少,水的质量分数35%时冻结速率先减少后趋于平稳,而后有变大的趋势
图9不同水的质量分数下冻结壁的冻结速率随冻结时
间的变化
Fig.9Variationoffreezingrateoffrozenwallwith
freezin
gtimeincaseofdifferentmoisturecontent
2.2灵敏度分析
1灵敏度定义.
工程上常将函数F(X关于分量xi的一阶偏导数称为函数F(X关于xi的灵敏度函数,记作[14]
Sxi(X=F(X
xi
(5
实际工程中,多数结构状态函数(如位移、内力、
稳定安全系数等为关于参数的隐函数,无法直接由式(5求得关于某一参数的显式灵敏度函数Sxi(X,故通常用差分法作近似数值计算.如欲求灵敏度Sxi(X0,则先对参数x0i作微小扰动!
xi,然后按下式计算.
Sxi(X0
=F(x0
1,!
x0
i+!
xi,!
n-F(X0
i
(6
此灵敏度定义表征了由于计算参数的变化引起计算结果的变化大小.考虑到各个影响因素的取值范围和单位的不一致,从而可能使得灵敏度无法比较大小,所以这里试图引入相对变化量的概念来消除这一影响.因为!
xi与xi变化域(xn
i-x0
i的比值具有相同的物理含义,所以重新将灵敏度定义为
Sxi(X0=
F(x0
ini-x0
i(7
它说明了由于某一计算参数相对其变化域的单
位变化引起计算结果的变化大小.
2灵敏度分析.
由以上试验结果,得出各影响因素下灵敏度大小分析如图10.
图10显示了冻结50d时各参数对冻结壁厚度灵敏度示意图,可见冻结管直径对冻结壁厚度的灵敏度为0.35~0.60;
冻结管间距对冻结壁厚度的灵
图10冻结50d时各参数对冻结壁厚度灵敏度示意图
Fig.10Schematicdiagramofsensitivityofeveryparametertofrozenwallthicknessafter50day∀sfreezing
125第3期毕贵权等:
126兰州理工大学学报[2][3]第35卷敏度的绝对值为0.72~0.76;
盐水温度对冻结壁厚度的灵敏度的绝对值为0.70~1.70;
水的质量分数对冻结壁厚度的灵敏度的绝对值为0.15~2.60.所以各因素对冻结壁厚度灵敏度的影响由大到小依次为盐水温度、冻结管间距、冻结管直径、土层含水量.郭宽良,孔祥谦,陈善年.计算传热学[M].合肥:
中国科技大学出版社,1988.徐士良.ANSYS在冻结壁温度场分布分析中的运用[J].安徽建筑工业学院学报,2005,13(2:
3436.汤国璋,王星华.温度场控制在多年冻土隧道施工中的作用[J].岩土力学,2007,28(3:
455460.丁德文,傅连弟,庞荣庆.冻结壁变化的数学模型及其计算[J].科学通报.1982,27(14:
875879.任松保,郑正泉.冻结过程中相变界面移动及传热的计算机模拟[J].华中科技大学学报,2003,31(2:
6163.汪仁和,徐士良.冻结壁温度场模型试验及其导热系数反分析[J].安徽理工大学学报,2003,23(4:
1822.汪东波.双排管冻结温度场分布规律理论与试验研究[D].淮南:
安徽理工大学,2002.曹荣斌.深井冻结凿井法数值模拟与工程计算[D].淮南:
安徽理工大学,2002.马巍.中国地层土冻结技术研究的回顾与展望[J].冰川冻[3][4]3结论在试验中对粉质粘土的冻结壁温度场进行数值[5][6][7]分析,从冻结管直径、冻结管间距、盐水温度、土层含水量4个因素对冻结壁厚度和冻结速率的影响进行研究,通过对各因素下灵敏度进行对比分析得出各因素对冻结壁温度场的影响大小.1同一时刻各因素对冻结壁厚度的影响大小依次为盐水温度、冻结管间距、冻结管直径、土层含水量.2同一冻结壁厚度下各因素对冻结时间的影响大小依次为盐水温度、冻结管间距、冻结管直径、土层含水量.3土层含水量对冻结壁厚度影响不大,对冻结时间较为敏感.4各因素对冻结壁厚度灵敏度的影响由大到小依次为盐水温度、冻结管间距、冻结管直径、土层含水量.参考文献:
[1]汪仁和,曹荣斌.双排管冻结下冻结壁温度场形成特征的数值分析[J].冰川冻土,2002,24(2:
181184.[8][9][10]土,2001,23(3:
218224.陈方方.地下洞室弹塑性仿真反分析研究[D].西安:
西安理工大学,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地铁 隧道 水平 冻结 施工 温度场 影响 参数 分析
![提示](https://static.bdocx.com/images/bang_tan.gif)