地震动力响应问题方法研究进展..docx

地震动力响应问题方法研究进展随着科学技术的不断发展，国内外学者和研究人员对边坡问题的认识也不断完善，特别是近几十年来，国内外对土石坝地震动力响应问题的研究取得了比较丰硕的成果，而且关于边坡地震动力响应问题的研究方法也逐渐完善。目前，最常采用的研究方法是拟静力法、Newmark 滑块位移法和动力有限法。静力计算的基础上，将地震作用简化为一个惯性力系，将其所产生的惯性力假定为一个恒定的静力，并将其作用在边坡潜在的不稳定滑体上，然后根据极限平衡理论，求出边坡的抗震安全系数，其核心是设计地震加速度的确定问题；随着对边坡动力问题认识的不断深入，最早把坝坡认为是变形体的是 Mononobe HA 等人，并开始了以变形体的观点来探究土质边坡的动力反应问题，并首次提出了剪切楔法模型的概念，发明了边坡地震反应分析的新方法-剪切楔法；随着边坡地震响应分析方法不断发展，到 20 世纪中期，Newmark 通过假定滑移面的方法来确定变形体的屈服加速度值，并采用动力分析手段来判定是否产生滑移，并估算其永久位移，然后根据潜在变形来评价土坝坝坡的动力稳定性，这就是非常有名的有限滑块位移法；周健、徐志英等发展了基于粘弹性本构关系的动力有限单元法，黄建梁等借用 Sarma 法进行了地震稳定性的动态理论分析，在同时考虑水平和竖向地震动基础上，给出了坡体临界加速度计算公式，建立了根据水平和铅直地震加速度时程估计坡体失稳的加速度、速度和位移时程的方法，解决了地震加速度时程的确定问题、地震过程中坡体抗滑强度的衰减问题和孔隙的动态响应问题及坡体稳定性的评价问题；薄景山建立了计算土质边坡地震反应及评价其动力稳定性的数值分析模型；我国学者王思敬较早的研究了岩体边坡的动力问题，通过振动模拟试验探索并建立了边坡块体运动的动力微分方程，通过数值积分求得块体滑动的动力学特征，即块体运动加速度和块体相对基岩的运动加速度、运动速度和位移曲线多种动力响应分析方法，大致可分为 3 种(贾俊)[10]：a.解析方法；b.物理模拟方法；c.数值模拟方法。l)解析方法边坡根据临空面的数目可以分为单面坡和双面坡(比如坝坡)。双面坡有两个自由面，它的动力反应三量(速度、加速度和应力)的分布规律可以采用解析的方法—剪切楔法来获得。一维剪切楔法是 1936 年由 Mononobe 等提出的，随后大量的文献对该法进行了改进，并把它推广到三维情形。对于顺层岩质边坡，大多只有一个自由面，因此其基本属于单面坡，解析的方法对单面坡是无能为力的。(2)物理方法物理模拟是科学研究的重要手段，能较为合理的揭示事物的本质，但是由于受实验材料、实验设备以及实验技术等限制，物理模拟无法保证模型与原型的真正相似。同时，物理模拟会存在尺寸效应问题，要研究边坡在整个剖面上的动力响应规律，在动力作用过程中必须在边坡体内布置大量的监测点。这对于物理模拟来说是非常困难的。(3)数值模拟随着计算机技术的发展，数值模拟己逐渐成为边坡动力问题研究中最重要的方法之一，它不但可以有效的弥补物理模拟的不足，而且计算结果明了直观，又可以节省大量的时间和金钱。现在常用的数值分析方法有：有限单元法、有限差分法、离散单元法等等。最早研究边坡动力问题的是土力学领域，当时是为了解决土石坝和堤坝在地震作用下的稳定性问题。在早期的关于土坝地震设计方法中人们往往错误地假定坝坡是绝对刚性体，因而在分析中采用拟静力法（pseudo-static approach）[11]。Leshchinsky 等[12]采用拟静力法来评价简单边坡的稳定，用数值方法计算了潜在滑动面上的正应力分布，用此正应力确定了满足所有极限平衡方程的最小安全系数，且提出了简单边坡地震稳定评估的设计表，此表在非地震条件下与 Taylor表相同。Ling 等[13]将拟静力法用于沿节理面滑动的岩体地震稳定性分析中，进行了地震稳定分析和永久位移计算。Siyahi 等[14]在正常固结土边坡地震稳定性分析中采用了拟静力法，采用参数分析确定了不同剪切强度的安全系数并考虑了剪切强度降低的影响。实践证明，用拟静力法设计往往低估含易液化土坡破坏的可能性，而对无液化可能的边坡，则往往高估其破坏的可能性。通常，地震动特性用峰值、频谱和持时三要素来描述，拟静力法的根本缺陷是未能考虑地震动的频谱特性和持时的影响。随着地震反应分析方法取得实质性进展，20 世纪 60～70 年代更多的注意力集中到发展可靠的步骤和标准来评价强震期间土坝的稳定性和安全性。拟静力法由于本身的缺陷显然已经无法胜任。为了解决这一问题，N M Newmark 注意到无论什么时候只要作用在潜在滑体上的惯性力超过了它的屈服阻力，滑动便发生；当惯性力改变方向时，滑动停止甚至向回滑动。Newmark 于 1965 年以屈服加速度 ay概念为基础提出了有限滑动位移的