基于ABAQUS路堤地震动力响应研究.docVIP

基于ABAQUS路堤地震动力响应研究摘要： 路堑路堤是我国西部山区常见的路基形式。该文结合工程实例，采用有限元数值分析软件ABAQUS对位于国道213DK732+560处坡度分别为1：1.5和1：1.75的半挖半填路堤进行了分析。本文首先选取5.12汶川地震中卧龙台站所记录的加速度时程，将有限元分析结果与Newmark法计算结果进行比较发现，对于坡度为1：1.5和1：1.75的路堤，有限元数值分析计算出的结果与PGA放大系数为1.25倍的Newmark法计算结果比较吻合，PGA放大系数取1.25倍时有更高的准确度。 Abstract： Cutting embankment is a common embankment form in China’s western mountains. Combined with an engineering example， the article uses the finite element analysis software ABAQUS to analyze the embankment located at National Road 213DK732+560. Firstly， the acceleration time history recorded at the Wolong stations is selected， by comparing the results from the finite element analysis， results show that the PGA amplification 1.25times is equal to and greater. If this coefficient is used to analyze embankments， the reliability of engineering design can be improved. 关键词： ABAQUAS；地震响应；路堤；PGA；Newmark Key words： ABAQUAS；Earthquake Response；Embankment；PGA；Newmark 中图分类号：U41 文献标识码：A 文章编号：1006—4311（2012）27—0077—02 0 引言 目前，在工程中常利用现有规范对高路堤边坡进行稳定性设计，大多采用拟静力法，即把地震力看作静力施加在路堤边坡上。这种方法有很大的缺陷性由于地震荷载的随机性，没有考虑地震波对结构稳定性的影响所具有的时程累计性，所以仅仅考虑峰值地面加速度PGA （Peak ground accelearation）对结构的稳定性影响是不够的。本文利用通用大型有限元软件ABAQUS对国道213DK732+560处半挖半填路堤边坡进行动力分析，并将分析结果与Newmark法进行对比分析，找出相关规律并得到了一些有价值的结论。 1 计算模型的确定和参数的选取 1.1 计算模型的确定 选取国道213DK732+560处半挖半填路堤为原型建模，公路设计等级为二级，路基所处的位置位于山腰上，路堤高20米，路基宽度8.5米，行车道宽度7米，土肩宽度2×0.75米。分别建立坡度为1：1.5和1：1.75的三维有限元模型进行计算对比分析，工况和有限元模型见图1。 1.2 计算参数的选择 1.2.1 材料的物理力学参数 1.2.2 材料的阻尼 本文将采用瑞利阻尼，以下将主要介绍瑞利阻尼。 1.2.3 粘弹性边界 本文中边界条件的选取在静力分析中主要采用固定边界，而在动力分析中采用的是人工边界中的粘弹性边界条件。 KBT=αT■；CBT=ρcs（1） KBN=αN■；CBN=ρcp（2） 式中：KBN、KBT分别为弹簧法向与切向刚度；CBN、CBT分别为阻尼器法向与切向的阻尼系数；R为波源至人工边界点的距离；cs和cp分别为S波和P波波速；G为介质剪切模量：ρ为介质质量密度；修正系数aT与aN取0.6和1.33。 1.2.4 地震波的输入 本文选择5.12汶川地震中卧龙台站所记录的加速度时程，为了节约计算时间，将汶川波持时控制在40s，保证其峰值加速度在此时间段内。为了能够尽量反映现实情况，本文考虑水平和竖直地震加速度对坡体的影响，汶川波水平和竖直地震加速度时程曲线如下图所示： 2 计算结果分析 2.1 有限元计算结果 通过ABAQUS数值分析计算得到路堤路面竖直方向随时间变化的沉降曲线，见图4。从图中可以看出，当路堤坡度为1：1.5时的永久沉降量为24cm左右，当路堤坡度为1：1.75时的永久沉降量为19cm左右，通过与现场破坏图1对比可以看出，对于半挖半填的路堤