三维抗拔桩土模型无剪应力（带命令流文件） [复制链接]

1、模型描述
    三维竖直水泥土抗拔桩土模型，按1/2对称建模，如下图所示。水泥土桩体内埋设了钢绞线（无分散板），钢绞线分上部自由段和下部锚固段，桩体周围为匀质土体，土体本构模型为D-P模型。为了简化分析，上部自由段钢绞线没有建，直接在下部锚固段顶作用了上拔位移UZ=100mm。同时假定钢绞线-水泥土桩体间、桩-土接触处完全粘结（通过共用节点），没有考虑初始地应力，假设竖向重力和上拔位移同时作用，命令流如附件。      
                                  

图片:抗拔桩土模型.jpg

2、模型存在的问题-桩周土体无剪应力和塑性变形
    该模型计算可以收敛，土体、桩体和钢绞线的竖向位移和应力也还正常，如下图所示。存在的问题是，不管桩顶上拔位移多大，桩周土体始终无剪应力和塑性变形，我认为这是有问题的，主要原因有二，一是按照力的平衡原理，本抗拔桩即为纯摩擦桩，桩周土体肯定是有剪应力的，二是把本模型的材料改为完全弹性的话，桩周土体却又有剪应力了！！！！这个问题困扰了我好久，一直得不到解决，现上传到本坛子，望有兴趣的学友解决下。
                                                

图片:正常的土体竖向位移.jpg

                                                                 正常的土体竖向位移
                                                

图片:1.png

                                                              不正常的土体剪应力（为0）
                           三维抗拔桩土模型命令流文件.txt (6 K) 下载次数:26                                                                                                        

   上述问题困扰了我很久，我也想了很多方法，其中，我认为没有剪应力的原因可能是单元类型不同造成的，因此进行了下面更加简单的模型试算，经过下面的试算可得出，土体没有剪应力的原因可能是单元不同造成的。详见“http://”blog.sina.com.cn/s/blog_9f5beef70102wqhh.html。需手工删除上引号。

我是大评委，评分给态度