断层倾角对隧道震动响应影响分析论文_黄博,张丽

重庆中检工程质量检测有限公司重庆水务集团股份有限公司

摘要：为了对比不同断层倾角对衬砌地震动响应的影响差异，分别取倾角为45°、60°、90°三种工况，采用有限元软件ANSYS进行动力分析，研究地震响应作用下不同倾角断层对隧道结构的影响规律。

关键词：隧道工程；断层倾角；有限元方法；动力分析；

在西部地震区山岭隧道的修建中常常会遇到断层，在断层处围岩结构松散、破碎，特别是在地震中岩体表现得极不稳定性，使得隧道衬砌与围岩之间出现脱空现象，衬砌得不到有力的约束，最终导致破坏。所以本文将采用有限元动力时程分析方法研究在地震作用下不同断层倾角对衬砌影响。

1有限元模型建立要点

有限单元法就是把物体，假想地划分成有限个简单形状的单元，用这种有限个单元的集合体来替代原来的物体。而各个单元之间则靠节点相连，节点相当于一个铰链，传递着单元之间的相互作用力，即为节点力。作用在节点上的外力称为节点荷载。节点力与节点荷载不同，前者是内力，后者是外力。单元的划分和节点的选择，除了根据物体的特点、承受荷载的情况、计算精度的要求以及考虑计算机的容量等因素外，很大程度上是人为的。

目前以有限元法[1]为基础理论的很多通用商业软件用在各个主要领域，其中ANSYS是使用较频繁的一种，也是本文采用的主要分析手段。

1.1单元选取

本文计算所选用单元介绍：

① Shell63壳单元

Shell63壳单元具有弯曲及薄膜特性，与平面同方向或者法线同方向的负载皆可承受。单元形式有平面四边形和三角形两种，并具有x、y、z位移方向及绕x、y、z轴旋转方向的6个自由度，包括应力强大变形能力。可选择连续性剪切举证，用于大变形（有限的旋转）分析。其中需注意的是Shell63壳单元的面积不可以为0，通常单元连接号码不正确就会造成该现象。单元的厚度为0或线性厚度变化至单元的四点时为0接皆不许可。四点所构成的单元必须为一平面，然而仍然允许少量的非平面，但会造成少许歪曲。适度的歪曲结果在输出时会有警告信息。若歪曲严重则会导致严重错误，此时最好使用三角形单元。

② Solid45（3-D）实体单元

Solid45（3-D）实体单元由8个节点结合而成，每个节点具有x、y、z位移方向的3个自由度。单元具有塑性、蠕变、膨胀、应力强化、大变形和大应变的特征。与壳单元相似。单元不允许零体积，且不可有扭曲而造成两块体积分离，这两点通常是因为节点连接错误造成的。单元有8个节点，但棱柱及角锥单元也可以接受。

③ Combin14弹簧单元

Combin14是一种无质量单元，可应用于模拟一维、二维或三维空间在纵向或扭转的弹簧—阻尼效果。当考虑纵向弹簧—阻尼时，该单元承受单向受拉或受压，每个节点具有x、y、z位移方向的3个自由度，不考虑弯曲及扭转。可用于模拟模型中的边界条件。

1.2边界条件理论及选取

在土体结构动力相互作用问题中，由于振动或波动要向无限远处传播，实际研究对象是一个半无限空间。而有限元的计算对象只能是有界域，故此在用有限元法研究结构在地震激励下的振动时，只是将感兴趣的部分人为地切割出来进行离散化[2]。这个切割面即称为人工边界，实际上它在原连续介质中并不存在，因此要对计算边界作一定的近似处理。

本文在对模型设置边界条件时，地表采用自由边界，模型左右及底部均采用粘弹性人工边界。

1.3本构模型选取

目前还没有一种本构模型能够模拟各种情况下土的动力非线性特征，并具有实际应用中必要的简单性。因此针对具体的问题应该选择合理而简单的本构模型。本文采用了D—P本构模型[3]，认为土体是一种完全弹塑性的材料，这种假设目前被我国工程界广泛应用。

1.4地震波选取

本次计算采用在地震反应分析时常用的EL-Centro波加速度曲线作为地震输入。美国EI-Centro波南北向加速度记录，1984年发生于美国Califljnia，峰值加速度341.7cm/s2，时间间隔0.02秒，持续时间30s。受计算条件限制，截取包含其最大峰值的前20s共1000个记录作为输入。