ABAQUS在排水固结地基处理中数值模拟分析.docVIP

ABAQUS在排水固结地基处理中数值模拟分析 　　摘要：本文通过ABAQUS有限元软件模拟分析砂井的排水固结模型，并与理论计算结果对照，两者在计算结果、沉降发展规律等方面具有一致性，表明Abaqus有限元软件可以准确的反映砂井处理软土的过程及结果，可作为排水固结地基处理的设计辅助工具。 　　关键词：ABAQUS 排水固结 数值模拟 　　 　　中图分类号： TU47 文献标识码： A 文章编号： 　　一、前言 　　ABAQUS是国际上功能强大的大型通用有限元软件，其拥有能够真实反映土体性状的本构模型，能够进行有效应力和孔压的计算，具有强大的接触面处理功能来模拟土与结构之间的滑移等现象，具备处理填土或开挖等岩土工程中的特定问题的能力，可以灵活、准确的建立初始应力状态，对岩土工程有很强的适应性。本文结合工程案例，通过对比分析、验证得出abaqus在排水固结地基处理中的数值模拟具有较高的准确度，可用于实际的设计、检算。 　　二、排水固结地基处理原理及理论计算 　　1、排水固结地基处理原理 　　排水固结法是处理软黏土地基的有效方法之一。该法的基本原理是对天然地基、或先在地基中设置砂井、塑料排水带等竖向排水井，然后利用建筑物本身重量分级逐渐加载，或在建筑物建造以前，在场地先行加载预压，使得饱和软黏土地基在荷载作用下，土体中的孔隙水被排出，孔隙体积减小，地基发生固结变形，同时随着超静水压力逐渐消散，有效应力逐渐提高，地基土的强度得到增长。 　　2、排水固结地基处理中地基固结度理论计算 　　固结度的计算是竖井地基设计中一个很重要的内容。因为知道各级荷载下不同时间的固结度，就可以推算地基土强度的增长，从而可进行各级荷载下地基的稳定性分析，并确定相应的加载计划。已知固结度，就可推算加荷预压期间地基的沉降量，以便确定预压的期限。 　　由于砂井地基固结问题求解是土力学中的一个经典问题，已经有一些理论解答，因此本文采用Hansbo理论解结果与Abaqus解对照分析。 　　其中，Hansbo理论解给出的z深度处的平均固结度计算公式为： 　　 　　考虑到时间因子，可上述公式可简化为 　　式中： ??渗透系数（㎡/s）； 　　，，； 　　为砂井影响半径，由于砂井理论计算模型为轴对称模型，故，为砂井间距。 　　三、Abaqus有限元中固结计算简介 　　Abaqus程序以Biot理论为基础进行固结有限元计算，通过导入参数化语言来进行计算。在固结计算中，对于土体单元，一般采用孔压单元，其与常规力学分析中的连续单元的大体构造类似，最大区别在于具备了孔压自由度，可以施加孔压荷载，用于模拟饱和、非饱和土地。而对于砂井单元，由于不考虑井料刚度，除渗透系数外，其他参数均与其所在土层相同。 　　在初始条件设置上，对于土体这种空隙介质，为了正确定义初始状态，必须要给出初始孔隙比、初始孔压和初始有效应力的正确分布。这是应为Abaqus在计算前通过迭代来建立与荷载和边界条件相对应的平衡状态，并作为后续分析的初始状态，如果给出的初始应力状态偏离平衡条件过大，可能造成迭代的失败或迭代建立平衡后的位移多大从而造成分析失败或者分析结果与实际相差过大。 　　而对于时间步长的选择，一般采用自动时间步长，这是因为在固结计算的后期，孔压的变化较小，如果采用固定步长，就会造成不必要的计算时间浪费。 　　四、算例简述及计算结果 　　假设地基处理采用砂井加固，砂井为正方形布置，间距为=2m，长度为8m，砂井半径=0.1m，渗透系数为=10-2m/s。地基土厚度L=8m，土体渗透系数为=10-8m/s。为简化起见，将土体和砂井材料取为弹性模型，弹性模量=104kPa，不考虑土体的水平方向变形，泊松比=0。地基土底面不排水，顶部完全透水，土体及砂井表面作用有瞬时施加的荷载P0=20 kPa。 　　1、Hansbo理论解 　　由上述算例模型，可以求得砂井影响半径为， ，，=10-5㎡/s，由此可以推算出时间因子与土体固结度的相互对应关系如表1所列。 　　表1 Hansbo理论解时间因子与固结度 　　 　　 　　2、Abaqus有限元分析 　　在Abaqus有限元分析过程中，取模型底部固定三个方向的约束，不排水；四周仅能容许竖向位移，同样不排水；表面位移自由，且是自由排水面。荷载瞬时施加，然后根据土体表面各个时刻的竖向位移来求得时间因子与土体固结度的关系。模型中土体与砂井均采用C3D8P单元，土体与砂井的初始孔隙比均为1，设置时间总长为10天。 　　 　　图1 固结结束时模型中央剖面处孔压的分布 　　由图1可以看出，越靠近砂井处及地表处的孔压越小，也证实了砂井加快排水的特性。 　　 　　图2 固结结束时沉降等值线云图 　　由,2可以看出