2. 西华大学 能源与动力工程学院,四川 成都 610039
2. School of Energy and Power Engineering, Xihua University, Chengdu 610039, China
盾构法隧道的防水关系到隧道功能的正常发挥.盾构法隧道管片接缝间的密封主要通过埋设在接缝间的橡胶密封垫压缩后的弹性反力防水,一般由橡胶密封垫与混凝土管片间的接触应力来表示橡胶密封垫弹性反力,并要求表面接触应力大于水压.橡胶密封垫通常采用三元乙丙橡胶,该橡胶具有优良的耐老化特性.然而,橡胶具有与高分子材料类似的黏弹性特性,在长期受压状态下,会出现力学松弛现象[1].在实际工程中,管片接缝宽度可近似认为保持不变,因此橡胶密封垫主要体现出应力松弛的力学行为.应力松弛的结果造成表面接触应力的衰减,从而对隧道的防水能力产生直接的不利影响,有可能无法满足隧道100年使用寿命的设计目标.
以往对盾构法隧道防水材料耐久性的研究,大多数是以拉伸力学指标(如硬度、伸长率、拉伸强度等[2])来评价防水材料耐久性.从橡胶密封垫的防水机理来看,通过研究接触应力松弛行为的时变特性来评价橡胶作为防水材料的耐久性更有实用价值[3].
由于原材料配比、加工误差等原因造成合成橡胶材料参数的不确定性.由文献[4-6]的研究结果可知,材料的超弹性属性参数C10和C01的变异对橡胶密封垫表面接触应力的影响较大,尤其C10是影响抗水能力的主要因素,但研究结果只能反映橡胶密封垫的短期力学行为,对橡胶密封垫耐久性产生的影响鲜有文献讨论.力学松弛特性还涉及材料的黏弹性参数.
本文中采用超弹性-黏弹性本构模型,选择橡胶密封垫材料参数作为随机输入变量,以表面平均接触应力为分析目标,运用Ansys-PDS随机分析模块计算各随机参数对橡胶密封垫应力松弛特性的影响.结合应力松弛试验,分析应力松弛行为的时变特性规律.
1 橡胶应力松弛本构模型Ansys软件中采用广义Maxwell模型表示黏弹性行为,黏弹性模型是由一个非线性弹簧(见图 1的A部分)和n个由弹簧和黏壶串联组成的Maxwell模型并联而成(见图 1的B部分)[7].A部分体现了与时间无关的超弹性,本文中采用Mooney-Rivlin超弹性模型;B部分体现了与时间响应密切相关的非线性黏弹性,本文中采用Prony级数形式黏弹性本构模型,该模型能较准确地表达橡胶的松弛特性[8].模型中Gi(i=0, 1, …, n)为剪切模量,ηi(i=0, 1, …, n)为黏壶系数.
弹性橡胶密封垫在使用过程中始终处于压缩状态,本文中所分析的断面(见图 2)最大压缩率约为35%,在不利工况下实际压缩率通常小于35%.因此,采用两参数的Mooney-Rivlin超弹性模型,该模型几乎可以模拟所有橡胶的力学行为,一般适用于应变约为100%(拉伸)和30%(压缩)的情况.
应变能密度函数表达式如下所示:
$ W = {C_{10}}\left( {{I_1} - 3} \right) + {C_{01}}\left( {{I_2} - 3} \right) $ | (1) |
式中:C10和C01为材料力学参数;I1和I2分别为Cauchy-Green应变张量的第一和第二应变不变量.
1.2 黏弹性本构模型在黏弹性单元中,体积形变和剪切形变行为是解耦的,黏弹性本构方程可写成如下形式[9]:
$ \sigma=\int_{0}^{t_{0}} 2 G(t) \frac{\mathit{\rm{d}} e}{\mathit{\rm{d}} t} \mathit{\rm{d}} t+\boldsymbol{I} \int_{0}^{t_{0}} K(t) \frac{\mathit{\rm{d}} \Delta}{\mathit{\rm{d}} t} \mathit{\rm{d}} t $ | (2) |
式中:σ为Cauchy应力;G(t)为剪切松弛函数;K(t)为体积松弛函数;e为应变偏量部分;Δ为体积应变部分; t0为当前时间;I为单位张量.式(2)中等号右边第一部分表示剪切形变引起的应力,第二部分表示体积形变引起的应力.
对于可视为不可压缩材料的橡胶来说,体积应变为零,因此式(2)中等号右边第二部分可以忽略不计,则橡胶的黏弹性本构方程可简化为
$ \sigma = \int_0^{{t_0}} 2 G(t)\frac{{{\rm{d}}e}}{{{\rm{d}}t}}{\rm{d}}t $ | (3) |
本文中采用Prony级数形式来表达式(3)中的剪切松弛函数G(t),如下所示:
$ G(t) = G(\infty ) + \sum\limits_{i = 1}^n {{G_i}} {{\rm{e}}^{ - \frac{t}{{{r_i}}}}} $ | (4) |
当时间t趋于无穷大时,G ∞ =G0, G0为图 1中纯弹性元件的剪切模量.定义相对剪切模量
$ {\alpha _i} = \frac{{{G_i}}}{{{G_0}}} $ | (5) |
式(4)两边同除以G0,则有
$ \frac{{G(t)}}{{{G_0}}} = 1 + \sum\limits_{i = 1}^n {{\alpha _i}} {{\rm{e}}^{ - \frac{t}{{{\tau _i}}}}} $ | (6) |
式(6)中τi被称作松弛时间,定义式为
$ {\tau _i} = {\eta _i}/{G_i} $ | (7) |
在有限元计算中,αi、τi即为反映材料黏弹性特性的参数.
2 有限元模型及计算参数的验证 2.1 有限元模型盾构法隧道接头橡胶密封垫和管片沟槽尺寸如图 2所示.模拟橡胶密封垫压缩采用的单侧压缩有限元模型如图 3所示.压缩过程中,混凝土管片沟槽及上压板可视为刚体,而橡胶密封垫可视为柔性体,沟槽、上压板与橡胶的接触采用刚-柔接触.
超弹性-黏弹性体的橡胶采用plane183单元;对于接触单元,采用targe169单元和conta172单元.
2.2 计算参数的选取及验证超弹性材料参数C10=0.700 MPa, C01=0.035 MPa[7].根据文献[10],图 1中B部分的黏弹性模型个数取2,即n=2,对应的黏弹性参数α1=0.119 696,τ1=241.820 8,α2=0.475 895,τ2=4.370 416.
对如图 3所示模型利用有限元分析软件Ansys进行确定参数的应力松弛计算,提取橡胶密封垫表面平均接触应力
从图 4可见,计算结果与实测结果拟合度较好.计算所得橡胶密封垫表面平均接触应力
由图 4可知,在开始阶段应力下降较快,随着时间的增加,下降速度逐渐减小,应力曲线渐趋收敛, 即橡胶密封垫上表面接触应力渐趋稳定.
3 各随机参数的统计特征选定橡胶密封垫材料的超弹性参数C10、C01和黏弹性参数αi、τi(i=1, 2)为随机参数,不考虑各变量之间的相关性.各变量随机分布和特征值以及统计特征如表 1所示.
将随机参数信息输入Ansys软件的概率设计PDS模块中,可获得橡胶密封垫防水性能指标——表面接触应力的统计特征,作为随机参数灵敏度分析的基础.
4 随机参数的灵敏度分析灵敏度分析是评价因设计变量或参数的改变而引起结构响应特性变化率的方法,可以分析随机因素对结构响应影响程度的大小.
灵敏度分析的方法有很多,本文中采用Spearman秩相关系数[11]来间接表示参数的灵敏度.Spearman秩相关系数可用于度量变量之间联系的强弱,如果计算所得系数接近1或-1,就认为输入变量对输出变量影响显著,如果计算所得系数趋近0,就认为影响微弱.如果系数为正值,表明增加输入变量的值,输出变量的值也增加;如果系数为负值,表明增加输入变量的值,输出变量的值反而减小.
选取橡胶密封垫上、下表面的平均接触应力为分析目标,分别用
(1) 各随机参数对上、下表面平均接触应力
(2) 以
(3) 表 2中,C10和τ1对应的灵敏度为正,表示C10和τ1与表面接触压力正相关,而其余4个参数对应的灵敏度为负,表明这些参数与表面接触压力负相关.这一点也可以从如图 5所示的
取松弛时间为1年,此时相邻2次应力计算结果差值小于10-4 MPa,说明应力松弛曲线已基本趋于稳定.分别作出超弹性随机参数、黏弹性随机参数以及综合考虑超弹性-黏弹性随机参数的最大极限值、均值和最小极限值下的橡胶密封垫表面接触应力曲线(见图 6),以探究应力松弛条件下橡胶密封垫防水能力的变化规律.
由图 6a~f可以看出,橡胶密封垫表面接触应力整个松弛过程可大致分为3个阶段:急速衰减期,发生在受压后的第1天;缓慢衰减期,在随后的100天左右基本完成;最后进入稳定期,得到后期相邻2次应力计算结果差值均小于10-4 MPa.
在只考虑超弹性随机参数时(见图 6a、b),图中曲线从上到下依次为随机参数取最大极限值下的应力曲线、取均值下的应力曲线、取最小极限值下的应力曲线,说明该随机参数取值越大,应力松弛程度越小,同时再次印证了超弹性参数与表面接触应力,也即橡胶密封垫防水能力正相关.
正好相反,黏弹性参数越大,应力松弛程度越大.
在综合考虑超弹性-黏弹性随机参数时(见图 6e、f),2种参数的作用相互抵消,3条曲线趋于重合.
6 结论(1) 随机参数对上、下表面平均接触应力
(2) 增大C10,
(3) 橡胶密封垫表面接触应力整个松弛过程可大致分为3个阶段:急速衰减期,发生在受压后的第1天;缓慢衰减期,在随后的100天左右基本完成;最后,进入稳定期.同时发现,超弹性参数越大,应力松弛越小,黏弹性参数越大,应力松弛越大,该变化规律对评价橡胶密封垫作为防水材料的耐久性具有一定的工程意义.
[1] |
刘礼华, 熊威, 张宏志, 等. 不可压缩超弹性止水材料的黏弹性计算方法研究[J]. 固体力学学报, 2007, 28(1): 83 LIU Lihua, XIONG Wei, ZHANG Hongzhi, et al. Viscoelasticity calculation method for incompressible hyperelastic waterstop material[J]. Acta Mechanica Solida Sinica, 2007, 28(1): 83 DOI:10.3969/j.issn.0254-7805.2007.01.015 |
[2] |
STÖGBAUER H, WOLF A T. The influence of heat aging on one-part construction silicone sealants[J]. Construction and Building Materials, 1991, 5(1): 27 DOI:10.1016/0950-0618(91)90029-K |
[3] |
伍振志, 杨林德, 季倩倩, 等. 越江盾构隧道防水密封垫应力松弛试验研究[J]. 建筑材料学报, 2009, 12(5): 539 WU Zhenzhi, YANG Linde, JI Qianqian, et al. Experimental study on stress relaxation of waterproof gasket for shield tunnel under river[J]. Journal of Building Materials, 2009, 12(5): 539 DOI:10.3969/j.issn.1007-9629.2009.05.008 |
[4] |
雷震宇, 刘明. 盾构管片弹性橡胶密封垫防水性能灵敏度分析[J]. 同济大学学报:自然科学版, 2017, 45(3): 336 LEI Zhenyu, LIU Ming. Sensitivity analysis of waterproof performance of elastic rubber seals for shield segments[J]. Journal of Tongji University: Natural Science, 2017, 45(3): 336 |
[5] |
雷震宇, 刘明. 具有随机参数的盾构管片弹性密封垫防水性能优化分析[J]. 应用数学和力学, 2017, 38(8): 899 LEI Zhenyu, LIU Ming. Optimization analysis of waterproof performance of shield segment elastic seal with random parameters[J]. Applied Mathematics and Mechanics, 2017, 38(8): 899 |
[6] |
雷震宇. 盾构隧道管片橡胶密封垫的优化设计方法[J]. 地下空间与工程学报, 2010, 6(4): 770 LEI Zhenyu. An optimal design approach for the rubber sealing gasket in shield tunneling[J]. Chinese Journal of Underground Space and Engineering, 2010, 6(4): 770 |
[7] |
于海富, 李凡珠, 杨海波, 等. 有限变形下的橡胶材料非线性高弹-黏弹性本构模型[J]. 橡胶工业, 2017, 64(11): 645 YU Haifu, LI Fanzhu, YANG Haibo, et al. Nonlinear high-elastic-viscoelastic constitutive model of rubber material under finite deformation[J]. Rubber Industry, 2017, 64(11): 645 DOI:10.3969/j.issn.1000-890X.2017.11.001 |
[8] |
BERGSTRÖM J S, BOYCE M C. Constitutive modeling of the time-dependent and cyclic loading of elastomers and application to soft biological tissues[J]. Mechanics of Materials, 2001, 33(9): 523 DOI:10.1016/S0167-6636(01)00070-9 |
[9] |
刘晓东, 张玉磊, 郭为君. 硅橡胶减震器振动试验的有限元模拟[J]. 噪声与振动控制, 2010, 30(5): 173 LIU Xiaodong, ZHANG Yulei, GUO Weijun. Finite element simulation of vibration test of silicone rubber shock absorber[J]. Noise and Vibration Control, 2010, 30(5): 173 DOI:10.3969/j.issn.1006-1355.2010.05.041 |
[10] |
冯希金, 危银涛, 李志超, 等. 未硫化橡胶非线性黏弹性本构模型研究[J]. 工程力学, 2016, 33(7): 212 FENG Xijin, WEI Yintao, LI Zhichao, et al. Study on nonlinear viscoelastic constitutive model of unvulcanized rubber[J]. Engineering Mechanics, 2016, 33(7): 212 |
[11] |
王涛. 基于Spearman秩相关系数的红外弱小目标检测[J]. 科学技术与工程, 2017, 17(2): 234 WANG Tao. Infrared weak small target detection based on Spearman rank correlation coefficient[J]. Science Technology and Engineering, 2017, 17(2): 234 DOI:10.3969/j.issn.1671-1815.2017.02.042 |