考虑孔隙比和渗透系数随土体当前应力变化深基坑降水开挖变形分析.docVIP

考虑孔隙比和渗透系数随土体当前应力变化深基坑降水开挖变形分析 　　摘 要：基于修正剑桥模型理论，推导了孔隙比e随土体当前应力变化的方程，同时比选了4组经典的描述渗透系数k随孔隙比变化的方程，选择了其中一组最佳的估算公式，编写ABAQUS用户子程序VOIDRI和USDFLD，以实现孔隙比和渗透系数随土体当前应力的变化。在此基础上，研究深基坑降水开挖所致的坑内外土体的变形、围护结构的变形及弯矩，得到以下结论：当考虑孔隙比随土体当前应力变化时，坑外地表沉降量、墙体的水平位移、地下连续墙的弯矩、坑底隆起量均大于孔隙比为定值时的情况；当考虑渗透系数随土体当前应力变化时，坑外地表沉降量、墙体的水平位移、地下连续墙的弯矩均小于渗透系数为定值时的情况，但同时考虑渗透系数和孔隙比变化情况时，其对坑底的隆起量的影响可以忽略不计。 　　关键词：孔隙比；渗透系数；用户子程序；深基坑 　　中图分类号：O319.56 文献标志码：A 文章编号：1674-4764（2015）02-0092-09 　　孔隙比和渗透系数的取值对土体的强度和变形产生影响较大。早期的岩土工作者在试验的基础上，考虑影响渗透系数的相关因素，提出了一系列针对不同土体e-k之间的关系表达式[1]。对于实际工程，若能掌握孔隙比和渗透系数随土体当前应力变化的规律，将使工程师能准确的计算如卸载、加载等不同工况条件下土体当前的孔隙比和渗透系数，为准确的预测土体的强度和变形奠定坚实的基础。为此，岩土工作者们在此方面做了不少有意义的研究。曾玲玲等[2]通过分析23组重塑黏土的固结渗透试验数据，提出土体压缩过程中当土体的应变大于20%时，渗透系数的对数与孔隙比的关系并非线性关系，而是呈非线性的指数关系。孙立强等[3]通过室内试验得出了吹填土孔隙比随固结应力变化的关系式及渗透系数与孔隙比的关系式，在此基础上建立了渗透系数与土体固结应力的变化关系。刘维正等[4]对太湖湖沼相粉质黏土原状样与不同工况下重塑样的渗透系数进行了测定，试验表明：对于同一土体，渗透系数主要受孔隙比的影响，而受土的结构性和应力历史的影响较小，并在此基础上建立了lg（1+e）-lgkv的渗透模型。邱长林等[5]通过理论分析和室内模型试验，对吹填土的真空预压固结的影响因素进行了研究，结果表明：相对传统计算方法而言，对于吹填土地基进行固结设计时，应当考虑孔隙比和渗透系数这两种重要因素的变化的影响，否则得到偏小的总体固结时间。 柯瀚等[6]以自制垃圾为研究对象，采用室内三轴渗透试验和常水头渗透试验，分析了孔隙比、降解龄期、有效固结应力等因素对自制垃圾的渗透系数的影响。 　　孔隙比和渗透系数是影响土体强度和土体固结快慢的重要参数，其取值对于基坑降水开挖这种复杂工况的工程有着重要的影响，继而影响基坑降水开挖所致的坑内外土体的变形及围护结构的变形。本文在前人研究的基础上，首先基于修正剑桥模型，推导出土体处于K0状态下孔隙比随土体当前应力变化的关系式，同时，比选了几组经典的渗透系数k随孔隙比变化的估算公式，选择其中一组最佳的估算公式，并据此编写了ABAQUS用户子程序VOIDRI和USDFLD，以实现孔隙比和渗透系数随土体当前应力的变化。在此基础上，分析了孔隙比和渗透系数随土体当前应力的变化而变化这种更符合实际的情况下，深基坑降水开挖对坑内外土体变形、围护结构的变形和内力的影响。 　　邵 羽，等：考虑孔隙比和渗透系数随土体当前应力变化的深基坑降水开挖变形分析 　　1 基于MCC的孔隙比和渗透系数随土体当前应力变化的关系 　　为了实现孔隙比随土体深度变化而变化，也即随土体当前应力的变化而变化，根据下图1中e-lnp′平面中的NCL线（各向等压固结线）、CSL线（临界状态线）与K0状态固结线之间的关系，并结合修正剑桥模型的原理推导出正常固结黏土K0状态下孔隙比e与有效平均主应力的关系式，其具体推导过程可参见文献[7]。 　　关于渗透系数与孔隙比之间的关系，不少研究者通过室内外试验提出了不同的经验关系式，其中运用较为普遍的主要为表1中4种关系式。 　　将上述所得的斜率Ck与e0的关系描绘在文献[15-16]中所统计的Champlain Clay、加拿大国内及一些国外粘土的Ck与e0的关系数据库中，如图2所示。从图2可以看出，Ck（M）和Ck（S）在数据较为集中的上方，Ck（K）处于数据集中的下方，其与数据库数据的相关性较差，而 Ck（T）与其他3式相比，其与数据库的相关性较好，所以本文拟采用Talay提出的渗透系数与孔隙比的关系式[14]，并结合前面推导的孔隙比随土体当前应力变化的关系式，可以得到下式： 　　式中e0为K0状态下的初始孔隙比，e1为各向等压固结下的初始孔隙比。 　　按式（6）编写子程序USDFLD，以实现渗透系数随