开发排水板单元问题的交流

sharp_liu_73

费康的书《ABAQUS在岩土工程中的应用》里面有较为详细的自己开发排水板单元的内容。我按照里面叙述的方法和步骤做了一个例子，但是发现总沉降被大大加大了，不知是何原因？按理说排水板只能加速排水固结的过程，也就是说只是加速沉降，而不能加大总沉降。
      而且，在计算过程中出现以下的warnings：User subroutine uel will be used with the stavev array dimensioned to one since the variable parameter has not been specified. Consequently, defining more than one statev entry in subroutine uel will cause code execution errors.不知这个有什么样的影响？
      有朋友在做塑料排水板数值模拟吗？希望多交流一下。

sharp_liu_73

经过进一步的研究，发现采用排水板单元后不是最终总沉降被大大加大了，只是在平衡地应力那步产生了很大的沉降，平衡地应力产生的沉降最大值大致为0.8m，而后续荷载产生的沉降只有0.06m。在不添加排水板单元时（天然地基），平衡地应力产生的沉降最大值大致为0.004m，而后续荷载产生的沉降为0.221m。所以采用排水板单元并未加大总沉降。附图是同一个节点采用排水板单元和不采用排水板单元两种情况下的沉降时程曲线（时间单位是天）。但是问题是为何添加了排水板单元会造成平衡地应力那步产生很大的沉降呢？而且，采用排水板单元的情况下在荷载保持恒定的情况下，为何还会出现地基隆起（沉降量减小）的情况呢？附上计算的相关模型文件，供感兴趣的同仁研究。

sharp_liu_73

为了方便大家参考，现将费康的书《ABAQUS在岩土工程中的应用》里面开发排水板单元的相关内容也一起附上。

sharp_liu_73

经过与同门师兄弟中的高手探讨，发现了以上问题的原因。主要是排水板单元在平衡地应力过程中的排水作用大大加大了沉降量，而实际情况是，在计算平衡地应力时不应该将排水板单元加入（应该是均质地基在漫长的地质年代过程中形成的）。采用单元生死技术，在geostatic步将排水板单元移除，在geostatic后面的实际荷载步激活排水板单元，效果就比较令人满意了，平衡地应力的沉降只有几个毫米。

sharp_liu_73

以上的做法沉降基本是比较正常了，但是感觉孔压还是不太正常，正在进一步的摸索当中。请关心这个问题的同仁共同关注此问题。我会把自己得到的最新结果提出来与大家共同分享。

sharp_liu_73

为何添加排水板单元后，在荷载施加过程中及其以后会出现局部的负的孔压？奇怪！还在摸索中。希望高手参与讨论。

sharp_liu_73

基础底面（是排水边界条件）结点孔压消散最终值有两种情况，一种是+47kPa，一种是-47kPa。基本上排水板处的都是+47kPa，其它点的有正也有负的。时间单位是天。添加了排水板单元，地基中出现负的孔压，无法理解啊！期望高手出现点拨一下！

sharp_liu_73

问题已经解决。就是不能用总孔压，只能按照超孔压来计算。也就是说，平衡地应力时不能施加gravity荷载，只能施加体积力。

凯迪_拉克

很有用的帖子，本人也在模拟塑料排水板，看了此贴才恍然大悟用单元生死功能来平衡地应力，现在是地应力平衡了，但load步却计算不收敛，用的是修正剑桥，希望能向楼主探讨下

something2005

楼主很有研究精神，赞一个！

LHP1211

sharp_liu_73 发表于 2011-3-2 16:47
为何添加排水板单元后，在荷载施加过程中及其以后会出现局部的负的孔压？奇怪！还在摸索中。希望高手参与讨 ...

楼主，能不能加个qq交流下，我的排水板堆载预压孔压计算结果也不正常，在坡脚正下方出现一个核，很奇怪，可不可以指点一下qq346998592