无限单元在地基动力反应中的应用研究.pdfVIP

无限单元在地基动力反应中的应用 岳庆霞1 史晓军1 冯明国2 (1同济大学建筑工程系上海200092；2山东科技大学设计研究院山东泰安271019) 摘要 在地下结构或土体的动力分析中，无限地基的模拟是其中的一个关键问题。本文对无限地基模 拟的人工边界方法进行了若干理论探讨，并借助于有限元软件ABAQUS中对各类人工边界进行了算例 模拟，最后将无限单元、粘性边界和修正的粘弹性边界的数值计算结果同理论解进行了比较。结果表 明，相对于其他人工边界而言无限单元的效果较好，且使用简单，是无限地基的动力模拟的有效方法。 关键词 人工边界；粘性边界；粘弹性边界；无限单元 1 引言 在对土体或地下结构进行有限元分析的过程中，必须对实际上近于无限的计算域用某一人为 边界截断，才能在有限的区域内划分有限个单元进行计算。这样就不可避免地给计算结果引入截 断误差。特别是对于动力分析问题，应力波传到截断边界时将发生反射，常常会给计算结果带来 不能允许的误差。这就使得动力分析过程中边界处理问题尤为重要，也吸引了很多学者进行这方 面的研究。 最初和最简单的人工边界是远置人工边界。1968，Altermanll】首次将无限域用一个很大的有 限域来模拟，这种远置人工边界的方法成为现在常用的检验边界精度的一个比较精确的方法。一 般认为【2】，当人工边界与计算结构的距离为结构特征尺寸5-10倍时，上述远置人工边界可以取 boundary condition)。粘性边界利用在边界上设置法向和切向阻尼器的方法吸收反射波的能量，从而减少边 界反射波引入的误差。粘性边界因为使用简单，因此直到现在都应用较多。后来的研究者又提出 了粘弹性边界的概念，即在粘性边界的基础上并联一个弹簧，并建立了若干修正模型[4】。 随着单向波动理论的发展，各国学者开始从波动方程的角度的对人工边界条件进行研究。 Lindman首先从偏微分方程构造单侧波动方程的角度来近似模拟人工边界外行波场，并依据这一 边界(Paraxial Approximation 标志性成果。1984年，我国学者廖振鹏【5】在直接模拟单侧波动共同运动学特征的基础上，提出了 多次透射人工边界(MTF)。这种边界条件，以边界内节点的位移来推求边界点的位移，使得人 工边界变成一种真正的可动边界，但是在实际应用过程中，这种边界给出的计算结果有时会出现 低频漂移和高频失稳现象，并且控制参数的选取也带有一定的经验性。 下面仅对粘性边界、修正粘弹性边界及有限元软件ABAQUS自带的无限单元进行简要的说 明，并对各种边界的计算结果进行对比研究。结果表明，ABAQUS自带的无限单元使用简单，且 具有较高的效率和计算精度。 ·--——356·-—— 2各种边界条件的简要说明 ，2。1粘性边界口1 condition)。粘性边界 boundary 1969年，Lysmer和Kuhlemeyer提出了粘性边界条件(viscous 最初的想法是把无限域分成内部域和外部域，从内部域向外部域行进的波应该全部通过，而没有 反射，因此，可以把外部域看作是一个波动吸收器或是一个无反射边界。基于这种思想，Lysmer 等提出了边界条件用以下公式表示： 翌 (1) f=pK舀 式中，盯，f分别为边界的法向和切向应力；P为介质密度；圪，K分别为P波和s波波速；也zi 分别为法向和切向速度。 从以上我们可以看出，粘性边界的实质就是在人工边界上施加法向和切向的阻尼器。 2．2修正的粘弹-眭边界Ⅲ① 杜修力和赵密利用平面波和远场散射波经验叠加来反映外行波传播，并在推导过程中考虑了 多角度透射的影响，然后引入无限介质的线弹性本构关系建立了一种应力人工边界条件。对粘弹 性边界的粘性系数和弹性系数各提出一个经验系数进行了修正。其中对二维问题，法向和切向的 弹簧和阻尼系数分别为：