桩土动力分析中的有限元分析.docx

《桩土动力分析中的有限元分析.docx》由会员分享，可在线阅读，更多相关《桩土动力分析中的有限元分析.docx（5页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、桩土动力分析中的有限元分析0 土一前言要真实反映土结构的动力相互作用，重要因素是要充分考虑土与结构材料接触界面的变形 和动力相互作用。土与结构材料接触界面上通常存在较大的剪应力，这是两种材料的弹性 模量相差很大，界面两侧材料变形不一致引起的，而且由于土体即使在小变形下也表现出 明显的非线性，导致桩传给土体的能量有部分被耗散掉，反过来也是如此，这正好体现了 土一结动力相互作用的非线性性质。目前，在土一结构相互作用方面，研究包括滑移和分 离在内的相对运动的影响已受到一定的重视，为了模拟土-结构材料界面的性状，人们提 出了一些本构模型，但主要集中于静力分析，也没有在试验室用试验方法测定本构参数1接触

2、问题的有限元单元ANSYS软件含有丰富的接触单元模型，有点一点接触、点一面接触、面一面接触等下面简单分析ANSYS接触单元的计算过程。接触问题的关键在于接触体间的相互关系，此 关系又可分为在接触前后的法向关系与切向关系。在法向关系方面上，必须实现以下两点: （1）接触力的传递。（2）两接触面间没有穿透。在切向上要考虑两接触面间的摩擦力的 作用。ANSYS通过两种算法来实现此切向与法向接触关系：即罚函数法和拉格朗日乘子法。1.1接触刚度的确定通过引入人为定义的罚参数来实现接触约束条件，在反复迭代过程中，近似地实现接触面 的互不侵入条件。在ANSYS软件中，罚函数是以接触刚度的形式表现出来，通过接

3、触刚度 在接触力与接触面间的穿透值（接触位移）间建立力与位移的线性关系，然后合并到整个 系统的平衡方程中求解。可见接触刚度越大，则穿透就越小，理论上讲当接触刚度为无穷 大时，可以实现完全的接触状态，使穿透值等于零。但是显而易见，在程序计算中，接触 刚度不可能为无穷大（否则会出现病态矩阵），穿透也就不可能真实达到零，而只能是个 接近于零的有限值。当采用此法来考虑土一桩接触问题，由于计算是采用迭代过程，为保 证收敛性要求采用较小的接触刚度，但这样会带来较大的穿透值，产生较大的计算误差， 这是罚函数法无法解决的矛盾。因此ANSYS在接触算法中还引入了拉格朗日乘子法。2. 2基于anasas的接触单元

4、拉格朗日乘子法与罚函数法不同，不是采用力与位移的关系来求接触力，而是通过把接触 力作为拉格朗日乘子，与接触单元的位移场和约束共同建立拉格朗日函数，通过用梯度法 、拟牛顿法、广义逆等特殊方法进行求解，从而获得位移场和接触力场。由于接触力场和 位移场采用分别插值，所以此方法不要求在单元划分时考虑节点配对问题。但由于拉格朗 日乘子场即接触力场的引入，增大了方程组的尺度和求解难度，给计算带来困难。正因为如此，ANSYS联合两个方法的优点又形成了罚函数+拉格朗日乘子方法，这种方法可 以有更多更灵活的控制，可以更快的实现一个需要的穿透极限。但是不管采用哪种方法， 接触单元的计算都是以机械接触理论为基础，不

5、考虑接触过程中的能量损耗。由于桩一土 动力相互作用过程，土体表现明显的非线性，会消耗桩一土接触过程的能量，所以如采用 ANSYS软件中的接触单元分析桩一土接触问题会带来很大的误差。目前很多岩土研究者都 比较青睐于用ANSYS软件中的接触单元来考虑土-结构的接触分析，本文认为ANSYS中的 接触单元用于土一结构相互作用问题的适用性有待于进一步验证，目前可以有条件地使用。2单元刚度矩阵Goodman单元基本原理，Goodman等(1968)在线弹性假定下，界面单元的应力与相对位移之间的关系可用下式表示：式中T相应的单元刚度矩阵可以表示为其中：式中t为单元在第三维方向的尺寸；L为单元的长度。目前有关

6、描述界面应力和相对位移 之间的关系的界面本构关系的研究还很不成熟，对界面的性状尚难作细致的模拟3修改的软件设计桩一土动力相互作用过程中，一般认为能量损耗主要是由于桩周土体产生非线性变化而造 成的。根据Winkler地基梁模型，能量损耗可以利用阻尼来模拟式中有了 Goodman单元的刚度矩阵和阻尼矩阵后，即可将其组合到系统整体刚度矩阵和阻尼矩 阵中去，本文有限元的计算结果是利用Winlison-。法求解得到，此方法很常用，在本文 中毋需介绍。4子单元的确定为了验证Goodman接触单元加上阻尼后对计算结果的影响。本文假设某一桩一土一结构系 统(见图3),四层混凝土框架结构建于一软弱场地上，结构层

7、高均为3 m,假设上部结 构处于弹性阶段工作，钢筋混凝土的弹性模量取为31X106 kN/m(1) 土体计算参数。以假设的基岩上的4层土作为计算模型，计算参数见表1。(2)程序简介。本文编制了土一桩一结构二维有限元计算程序，桩和结构被简化为梁单 元梁单元采用四结点八自由度模型(含两个虚自由度)，梁单元两端分别有三自由度(两 个平动自由度和一转动自由度)和一虚自由度，虚自由度的设置是为了利于与土单元的四 节点八自由度统一，虚自由度上没有力和位移出现。土体划分为四边形单元，土单元与桩 单元之间加上修改的Goodman接触单元，人工边界采用粘弹性边界，其弹簧刚度和阻尼的 计算参考文献(3)输入的地震

8、动时程：Taft地震，见图4,震动持时为60 s,峰值加速度为0.15g, 为更好地分析Goodman接触单元中阻尼项对上部结构的影响，本文采用了对同一地震动， 通过调整加速度峰值的方法，得到不同强度的地震动输入，除了原有的0.15g,加速度峰 值分别调整到0. 2g和0. 3g两种工况，分别进行计算，研究了 Goodman单元在考虑阻尼的 情况下，对上部结构地震反应的影响。表2反应了土一桩一结构系统在同一地震不同的加 速度峰值作用下，接触单元考虑与不考虑阻尼的上部结构顶层处最大弯矩、最大剪力、最 大水平位移的比较，得出三种工况中上部结构在接触单元有阻尼时比无阻尼时最大反应都 表现减少的趋势。

9、上部结构顶层柱的最大弯矩、最大剪力、最大水平位移、峰值加速度在 加速度峰值0. 15g的情况下,有阻尼时，其值比无阻尼时分别减小7%, 4. 3% 11.6%、9%, 在加速度峰值0.2g的情况下，有阻尼时其值比无阻尼时分别减少了 8. 4%、7.1%, 12.78% ，11.4%,在加速度峰值0.3g的情况下，有阻尼时，其值比无阻尼时分别减小10虬12.3% 、15.13%、12.5%。而且不仅结构顶层的地震反应出现有阻尼时比无阻尼时减少，结构其 它层的地震反应也相应出现减少。从分析比较中还看出，在上调输入加速度的幅值，使用 考虑阻尼的Goodman接触单元后，上部结构的最大地震反应减小的幅

10、度更大一些。这与理 论上的解释也是相符的，由于输入地震动强度的增加，桩一土间的动力相互作用更加明显, 土的非线性发展更加充分，塑性耗能更加强烈，表现在土给桩和上部结构的能量要减小更 多，导致上部结构最大地震反应减小的幅度也相应增大。但笔者在计算中发现，当考虑接 触单元的阻尼项时，并不是结构各层的地震反应总是随地震加速度峰值的增加而减少，由 于桩土之间动力相互作用的问题是相当复杂，当地震加速度峰值增加到一定程度时，土一 桩一结构系统都处于严重的非线性状态，本文的一些假定不再适用，进一步的分析，将另 文报道。不过在中小地震作用下，本文提出的Goodman接触单元的修改模型，能够反应桩 土之间动力相

11、互作用过程中的能量损耗。5关于桩一结论与建议本文研究分析了桩一土一结构体系动力相互作用中的接触问题，指出了采用ANSYS软件中 接触单元分析桩一土接触问题存在的缺陷，并在工程研究中常用的Goodman接触单元的基 础上，提出一种不仅能够考虑桩一土动力相互作用中力的传递过程，又能考虑桩一土动力 相互作用过程中的能量损耗的改进的Goodman接触单元模型。通过分析比较得出以下结论:通过对ANSYS软件中接触单元的工作机理进行分析比较，发现ANSYS软件中的接触单元的 计算只反应了接触单元的力和位移的关系，不考虑接触单元的能量损耗，由于在桩一土相 互作用问题中，土体会表现明显的非线性并消耗一定的能量

12、。所以认为这种只考虑力和位移的相互关系的机械接触理论的接触单元，运用于桩一土接触界面分析中会带来一定的误 差。本文还分析了工程接触问题中常用的Goodman接触单元工作机理，同样存在着只考虑接触 单元与相邻单元的力的传递，不考虑动力相互作用中能量的损耗问题，所以用于分析土一 结构动力接触问题时，需要进一步修正。本文在Goodman单元的基础上，保留了此单元切向刚度和法向刚度的计算式，同时加上了 阻尼矩阵来表征土-结构动力相互作用中的能量损耗。阻尼矩阵的选取参照了 Winkler地 基梁中弹阻模型来表示桩一土的相互作用，认为接触单元中的阻尼必须能够反应桩周土的 非线性性质，又必须考虑到接触单元的

13、接触刚度影响，提出了阻尼矩阵的计算公式。编制了桩一土一结构动力相互作用的有限元计算程序，并分别计算了 Goodman单元有阻尼 和无阻尼的两种情况，发现在中小地震作用下，分析当接触单元含有阻尼时，上部结构的 最大地震反应会有所减小，可见采用原有的Goodman单元计算在一定程度上夸大了上部结 构的最大地震反应。通过分析输入地震动加速度峰值为0. 15g, 0. 2g和0. 3g的三种工况, 得出随着输入加速度幅值的增加，上部结构的最大反应减小的幅度会更大一些。存在问题和对今后研究的建议：到目前为止，对动力荷载下的桩一土界面性状的实验研究工作开展比较少，也没有合适的 桩一土界面的动力本构模型可用，这在一定程度上限制了土一结构动力相互作用研究和实 际应用，因此需要加强这方面的工作，为找到合适的动力接触刚度提供理论依据。本文对动力荷载下桩土界面存在的问题进行了一些分析，建议了一种Goodman修改模型， 并对比较简单的结构进行了一些比较分析，获得了以上初步结论，但其适用范围和计算精 度还有待于实践中进一步验证，并在此基础上发展更合理的，适合复杂结构实际应用的三 维分析方法。