某高速公路路基abaqus有限元分析报告.docVIP

《高等土力学》课程考察任务书 问题描述：选取某高速公路路基断面，宽度为12m。路基土层共分为三层，每层厚度均为4m，第一层土为砂土，其弹性模量E=（1+0.01×学号后两位数字）×kPa，第二层为黏土，其弹性模量E=（5+0.01×学号后两位数字）×kPa，第三层为砂土，其弹性模量E=(2+0.01×学号后两位数字)×kPa。三层土的其他有关基本物理特性参数由学生自己查找相关文献资料来确定。假设在高速公路路基断面上作用了两个对称的均布荷载q，宽度为0.5m，荷载作用边缘离断面中轴线的最短距离为2m。并且，均布荷载为一个呈正弦函数规律变化的动荷载，其幅值为（100+学号后两位数字）kPa，频率由学生自己确定。假设在均布荷载作用期间，土的基本物理特性参数均为恒量，不随时间改变。 基本要求：学生采用现成的大型商业软件，如ABAQUS,ANSYS,PLAXIS以及FLAC等。针对上述案例，建立自己的数值模型，模型可以是2D的，也可以是3D的。模型的本构关系要求在MC模型、Cam(clay模型、修正的Cam(clay模型、DP模型中进行选择。通过数值模拟，学生可以对路基断面中轴线上各点进行动力响应分析。要求在中轴线上选取5个关键点，分别为路面以下2m，4m，6m，8m，10m处。提取这五个点的位移响应曲线、速度响应曲线、应力(应变关系曲线。最后，提交一份完整的计算书，同时附上数值模拟计算的源程序。 某高速公路路基ABAQUS有限元分析报告 一、问题描述 选取某高速公路路基断面，宽度为12m，路基土层共分为三层，每层厚度均为4m，选用ABAQUS建立2D模型进行有限元分析。 （1）第一层土为砂土，其弹性模量E=1.21×107kPa，泊松比v=0.3，密度为1600kg/m3，本构关系采用摩尔-库伦模型，取其摩擦角为30°，膨胀角为30°，粘聚力为0； （2）第二层土为黏土，其弹性模量E=5.21×105kPa，泊松比v=0.4，密度为2000kg/m3，本构关系采用摩尔-库伦模型，取其摩擦角为20°，膨胀角为20°，粘聚力为45kPa； （3）第三层土为砂土，其弹性模量E=2.21×107kPa，泊松比v=0.3，密度为1600kg/m3，本构关系采用摩尔-库伦模型，取其摩擦角为30°，膨胀角为30°，粘聚力为0。 在高速公路路基断面上作用了两个对称的均布荷载q，宽度为0.5m，荷载作用边缘离断面中轴线的最短距离为2m。并且，均布荷载为一个按正弦函数规律变化的动荷载，其表达式为：，单位为kPa，频率ω=20s-1，其中，为了使该动荷载一直为压荷载，故加入了一个初始幅值121kPa。假设在均布荷载作用期间，土的基本物理特性参数均为恒量，不随时间改变。 通过数值模拟，我们可以对路基断面中轴线上的各点进行动力响应分析，在中轴线上选取5个关键点，分别为路面以下2m、4m、6m、8m和10m处，提取这五个点的位移响应曲线、速度响应曲线、应力(应变关系曲线。 二、ABAQUS建模过程 该三维高速公路路基的空间问题可简化为平面应变问题，因此，首先创建一个12m×12m的平面矩形部件，通过“拆分面：草图”工具将其按土层情况拆分成三层，如图1所示。再对其分别定义材料属性。 图1 创建部件 需要定义的材料属性包括：密度、弹性（弹性模量和泊松比）以及摩尔-库伦塑性，其中摩尔-库伦塑性中需要定义的参量分别为：摩擦角、膨胀角、粘聚力以及塑性应变。因为在ABAQUS的摩尔-库伦模型中要求粘聚力必须大于0而砂土的粘聚力为0，所以在此将砂土的粘聚力取一个趋近于0的较小值为1e-10kPa。对于塑性应变的定义方法，若不指定塑性应变，则ABAQUS认为其粘聚力保持不变，即为理想线弹塑性模型，本文采用此办法进行定义。其他材料属性均按问题描述中的参量一一定义完成。 装配好模型后，设置显式动态分析步，时间长度设为30，打开几何非线性开关，采用自动类型的增量，按此定义方式，模型分析增量步步长将控制在1.5左右。因为最终分析结果需要提取指定点的位移响应曲线、速度响应曲线、应力(应变关系曲线，所以按照此要求定义相应的场变量输出。 接下来，定义边界条件和载荷。点取本模型底边，设置成固定铰支约束；点取本模型左右两边，限制其水平方向的位移。因载荷按正弦规律变化，首先须创建周期型幅值曲线。 周期型幅值曲线用傅立叶级数表示，其表达式如下： 其中，A0为初始幅值，N为傅立叶级数项的个数，ω为圆频率，t0为初始时刻，An为cos项系数，Bn为sin项系数。 本模型正弦荷载幅值曲线设A0=121，N=1，ω=20s-1，t0=0，An=0，Bn=1。 选定荷载作用位置设线荷载大小为121kPa，方向竖直向下，选择上一步创建好的幅值曲线，如图2所示。 图2 施加正弦荷载 得到本模型示意图如图3所示