层状场址自由场动力分析等价线性法研究.docVIP

层状场址自由场动力分析等价线性法研究 　　摘要：基于等价线性法和一维剪切波传播理论，开发了等价线性法频域计算的C++程序。为验证自编C++程序的有效性，与在工程中广泛使用的SHAKE91程序就某一核电成层土地基进行算例对比，发现两种程序计算所得加速度峰值随高程的变化吻合较好，地震波时程向位移频谱转换时采用不同方法会对计算结果产生影响。 　　关键词：土层地震反应分析；等价线性频域方法；地震波时程处理；C++等价线性程序 　　中图分类号：TU352 文献标识码：A 文章编号：1000-0666（2016）01-0126-05 　　0 引言 　　土体是非线性很强的材料，承受地震荷载时，土的应力应变关系形成封闭的滞回圈。而各个滞回圈的形状、大小和方位都是变化的。为了简化地震作用下土的地震响应分析，Idriss和Seed（1968）最早提出了等价线性方法。由于该方法概念清晰，计算效率高，且计算精度能满足工程需要，目前被工程上广泛应用，已被我国《工程场地地震安全性评价技术规范》（GB17741-2005）推荐使用。 　　在与其他理论结合方面，吴再光等（1988）将等价线性法与随机理论相结合，解决了一维水平成层地基的地震反应分析问题；吴再光等（1989）将该方法推广到二维，提出了与随机理论相结合的概率平均等价线性法，随后，吴再光等（1992）又提出了水平成层地基非平稳随机地震反应分析方法，将概率平均等价线性法和时域模态分析相结合。为提高等价线性法的精度，针对其在强震作用下得到的地面反应与实测结果相比较小，且高频段的频响放大倍率较低的问题，学者们进行了一系列的研究：Yoshida等（2005）将频率分为3段，在每个频率段根据其特点采用不同的求等效剪应变抗解公式；沈建文等（2013）提出了用震级和距离参数修正土体的等效剪应变，使等价线性法更加完善；王伟等（2013）提出了刚度与阻尼频率相关的等效线性化方法，成功地改善了传统方法压制土体地震动高频成分的缺点，使结果更趋合理。 　　在实际工程应用中，学者们将等价线性法推广到不同的工程领域：徐学燕等（2003）将等价线性法应用于冻土和软土地基中，黄雨等（2005）利用Edushake程序对上海软土进行了分析，得到了该场地的频谱特性规律和地震加速度反应；刘海笑和王世水（2006）将等价线性法成功运用于结构-波浪-海床耦合系统的动力响应分析中。 　　本文通过等价线性法和一维波动传播理论，基于C++平台开发了等价线性法C++程序。对比了地震加速度时程到位移频谱的不同转换过程对结果的影响，并与SHAKE91进行了算例对比，验证了笔者开发程序的准确性。 　　1 等价线性法迭代过程 　　土是一种非线性材料，其主要参数剪切模量G和阻尼比D都是剪应变γ的函数，如图1所示。地震动越强，土的非线性特征越明显。在强震作用下，由于土体的非线性效应，叠加原理在动力求解的过程中不再适用。等价线性法在频域范围内用线性分析方法近似获得非线性地震响应，即假定在整个地震过程中土层的剪切模量G和阻尼比D均为常数，根据初始的G和D完成一个动力计算过程后，得出一个等效剪应变。根据新的剪应变插值出新的G和D，重复相同的计算过程，直到结果收敛误差达到允许范围为止，这就是等价线性法的非线性迭代过程。 　　2 基于Visual C++平台的软件开发与应用 　　Visual C++平台具有较高的计算和执行效率，且能与MATLAB结合进行数据库调用。为加深对等价线性法的理解，基于Visual C++开发平台，结合一维剪切谐波传播理论、傅立叶变换理论和露头基岩和下卧基岩地震运动的比例关系编制了等价线性法频域计算程序。 　　2.1 一维剪切谐波传播理论 　　图2显示了一维等效线性法场地剪切波传播反应分析模型，竖直向上的水平剪切波在一维成层土体中传播，则质点的水平向变化可表示为 　　2.2 算法描述与程序实现 　　在很多工程情况下，场地土层的波动传播性质沿横向变化幅度比沿深度小，因此在求解土层地震反应分析时主要考虑土体性质沿深度的变化。通常假定地震输入为竖直向上入射的水平剪切波，将土层的地震反应问题简化为一维的波动问题。通过一维剪切谐波传播理论和等价线性法理论，编制了水平成层覆盖土场地地震响应分析的频域计算程序，其计算流程如图3所示。 　　该程序通过输入土层编号、厚度、密度、剪切模量、剪切波速和泊松比等参数信息，结合等效剪应变与剪切模量、阻尼比关系曲线，输入地震波并进行土层地震响应的迭代计算，直至收敛。在进行计算时，可以选择露头基岩、内部层两种地震动输入方式，以通过露头基岩和下卧基岩地震运动的比例关系理论，实现地震波任意层输入所得到的土层地震响应结果。其参数输入界面如图4所示。该程序可输出各土层的位移、速度、加速度、最大剪应