软黏土广泛分布在我国沿海、沿江等经济发达地区,其具有含水率高、孔隙比大、压缩性强、黏滞性等特点,这使得在软黏土地基上修筑结构物时,其沉降问题尤其是长期沉降问题就显得尤为突出。软黏土的变形过程与时间密切相关,这既包括超静孔隙水压力消散、土体体积压密的固结过程,还包括由于其黏滞特性所引起的蠕变变形,两者互相影响,共同决定着土体的长期沉降。因此,如何基于土体变形的时间效应,建立起可以反映软黏土蠕变特性的本构模型并实现蠕变-固结的耦合作用,是工程界和学术界共同探讨的热点问题。本研究使用室内试验、理论分析及数值模拟等方法对软黏土的固结、蠕变特性及控制方程的建立进行了综合研究,主要研究内容及成果如下:(1)开展了一维K0侧限固结试验、三轴固结不排水剪切试验、三轴固结不排水蠕变试验,分析了软黏土变形的时效性特征和变形机理,并为进一步弹黏塑性本构模型的建立提供了试验基础。试验成果表明:1)土体的次固结系数随着加载时间的增长及预压作用而减小;且次固结系数随荷载的增加而增大,但其增加幅度依赖于土体的固结状态,次固结系数随荷载变化曲线会在与pc有关的某一范围处出现转折点。次固结系数与压缩指数之间呈较好的线性关系,其比值在0.028~0.049之间。2)当上部荷载较小时,土体的三轴固结不排水蠕变变形随时间会逐渐减弱并最终趋于稳定;而当荷载达到一临界强度时,土体的蠕变变形会在加载一段时间后随时间快速发展,土体发生剪切蠕变破坏。3)在应变控制模式下,土体的不排水剪切强度会随着剪切应变速率的提高而增加,且其率相关性可以很好地用速率参数表征。(2)基于一维K0侧限固结试验的成果,在Bjerrum等时间线理论的基础上验证了将蠕变变形当作与时间相关塑性变形这一假设的合理性,并分别从蠕变效应及应变率效应两个角度推导了一维K0侧限条件下的非线性弹黏塑性模型,给出了模型参数的物理意义及获取方法。将所建立的非线性弹黏塑性模型引入一维Terzaghi固结理论,建立了一维K0侧限条件下的蠕变固结耦合控制方程,并通过Crank Nicolson有限差分格式对控制方程进行了求解。通过对一维K0侧限固结试验的模拟及与Imai和Leroueil的试验结果对比,验证了所建模型的合理性。计算结果表明,随着土样高度的增加、土体黏滞特性的增强,本文模型与假设A计算所得的主固结完成时的变形量之间的差值也会越来越大。(3)基于改进的GJY型K0固结仪开展了一维柔性等比侧限固结试验,分析了土体在允许一定侧向变形时的压缩特性,并将试验结果与K0侧限条件下的成果进行了对比分析。结果表明,柔性等比侧限条件下的压缩曲线、压缩指数及次固结系数的发展规律均与K0侧限条件下的一致,但由于侧向变形的产生,其压缩指数约为K0侧限条件下的1.08~1.16倍,次固结系数约为K0侧限条件下的1.1~1.3倍。在柔性等比侧限固结试验的基础上,根据达西定律、有效应力原理及连续性原理推导出柔性等比侧限条件下的一维固结方程,且与建立的非线性弹黏塑性模型进行耦合,得到柔性等比侧限条件下的蠕变固结耦合控制方程,并给出了模型参数的意义及确定方法。通过对柔性等比侧限压缩试验结果的模拟,验证了模型的合理性,且模拟结果表明,柔性等比侧限条件下土样厚度、次固结系数对主固结完成时应变值的影响均较K0侧限条件下的大,本文对单层地基在柔性等比侧限条件下计算所得的主固结完成时的应变值约为K0侧限条件下的1.17倍。(4)为拓宽蠕变固结共同作用的路径,本文还基于室内蠕变试验成果建立了含分数阶导数的Merchant经验模型,并对一维K0侧限固结试验及三轴固结不排水蠕变试验成果进行了拟合,给出了模型参数的取值范围。该模型可以回避蠕变固结相互作用的复杂机理,从宏观上预测实际工程的长期变形。计算结果表明,与传统的Merchant模型相比,分数阶导数Merchant模型能更好地描述软黏土的蠕变变形,尤其是对于模拟试验曲线拐点处的变形具有明显优势,且能更好地反映土体变形随时间发展的非线性化程度。模型参数分析表明,在各级荷载作用下阶数α的变化范围最小,而黏滞性系数η1的变化范围最大,但从精度上总的来说可以满足工程计算的要求。(5)基于Perzyna过应力理论,引入参考面的概念建立了一般应力条件下的三维非线性弹黏塑性模型,并考虑了初始各向异性的影响。通过将所建立的本构模型引入至Biot固结理论,实现了三维条件下的蠕变固结耦合,并推导出了应力增量及增量形式的雅克比矩阵的表达式,编制了Abaqus有限元软件中的UMAT子程序,对软基上堆载问题进行了模拟。计算结果表明,由于土体的蠕变作用,即使外荷载保持不变,当排水条件较差时,孔隙水压力随时间也会有一定幅度的增加,并对地基的变形和长期变形产生实际影响。