桩土桩相互的作用有限元接触分析.docVIP

桩土桩相互的作用有限元接触分析 　　摘要:桩土体作为一个共同工作的系统,广泛存在于土木工程实践中,是典型的接触问题之一,对桩-土-桩相互作用的研究也是工程十分关心的,其中桩身摩阻力的分布更是关键所在。本文基于有限元数值分析方法软件对此进行了深入研究。 　　Abstract: Pile and soil which work together as a system are applied in the engineering practice of civil engineering, and this problem is also one of the typical contact topics. And in pile engineering it is always concerned about the distribution of friction resistance at the lateral pile.This paper discusses these issues in depth based on the numerieal method of finite element. 　　关键词:有限单元法;接触非线性;桩土相互作用;桩侧摩阻力 　　Key words: finite element method;non-linear contact;pile-soil interaction;friction resistance at the lateral pile 　　中图分类号:TU43 文献标识码:A文章编号:1006-4311(2010)11-0108-02 　　 　　0引言 　　桩土相互作用问题的实质是固体力学中不同介质的接触问题,具体表现为材料非线性、接触非线性等。目前,有限单元法是解决复杂空间结构静、动力问题、弹塑性问题最有效的数值方法之一。本文对桩土相互作用中接触问题进行分析时主要采用接触非线性有限元法,利用ABAQUS有限元软件进行研究。 　　1ABAQUS软件概述 　　ABAQUS是功能强大的有限元法软件[1,2],提供了广泛的功能且使用起来十分简明。对于非线性分析,ABAQUS能自动选择合适的荷载增量和收敛精度,且拥有十分丰富的、可模拟任意实际形状的单元库。 　　2ABAQUS桩土接触分析中需解决的问题 　　2.1 单元类型的选择在接触模拟中采用二阶单元会引起接触面上等效节点力的计算出现混淆,因此接触面两侧的单元一般不宜采用二阶单元,只能采用线性单元。 　　2.2 主从接触面的建立可以通过定义接触面(surface)来模拟接触问题,本文所涉及的桩土体之间的接触面主要有两类:①桩侧单元构成的柔性接触面(桩侧土体表面)或刚性接触面(桩表面);②桩底土体一般采用节点构成的接触面,选取桩底土体节点时,不包含己定义在柔性接触面上的节点。 　　在模拟过程中,接触方向总是主面的法线方向,从面上的节点不会穿越主面,但主面上的节点可以穿越从面。一般遵循以下原则:①应选择刚度较大的面作为主面,对于刚度相似的两个面,应选择网格较粗的面作为主面;②主面不能是由节点构成的面,并且必须是连续的;③如果接触面在发生接触的部位有很大的凹角或尖角,应该将其分别定义为两个面;④如果两个接触面之间的相对滑动小于接触面单元尺寸的20%,选用小滑动,否则选用有限滑动。有限滑动应打开几何非线性开关。 　　2.3 接触属性的定义接触表面间的相互作用包括两部分:一是垂直于接触面,另一是沿接触面的切向。①接触面法向性质。当两个面之间的间隙变为零,接触约束就起了作用,当接触压力变为零或负值时,接触面分离,约束就被撤出。这个行为称为“硬”接触。②接触面切向性质。接触面接触时,接触面间要传递切向力,需要考虑阻止面之间相对滑动趋势的摩擦力。库仑摩擦[3]是常用的描述接触面的相互作用的摩擦模型,用摩擦系数来表示接触面之间的摩擦特性。 　　3桩土作用分析中的几个关键问题 　　3.1 初始地应力场的计算在大多数岩土工程问题中,土体的初始应力场即为自重应力场,水平应力与竖向应力为[3]: 　　σz=γz(2) 　　σx=σy=k0σ0(3) 　　式中,k0为静止侧压力系数,k0=1-sin′或k0=μ/1-μ。 　　 初始应力场的设定可通过下列命令进行设定: 　　* initial conditions,type=stress,geostatic 　　Setname,stress1,coord1,stress2,coord2,k0 　　3.2 桩土之间接触摩擦系数的确定可根据经验公式估算出桩土界面的摩擦角δ,进而确定桩土间的摩擦系数。研究表明:对于粘性土取δ/′=0.6~0.7是比较合适的,下