高等土力学一课程考查任务书解析

《高等土力学》课程考察任务书

问题描述：选取某高速公路路基断面，宽度为12m。路基土层共分为三层，每层厚度均为4m，第一层土为砂土，其弹性模量E=（1+0.01×学号后两位数字）×107kpa，第二层为黏土，其弹性模量E=（5+0.01×学号后两位数字）×105kpa，第三层为砂土，其弹性模量E=(2+0.01×学号后两位数字)×107kpa。三层土的其他有关基本物理特性参数由学生自己查找相关文献资料来确定。假设在高速公路路基断面上作用了两个对称的均布荷载q，宽度为0.5m，荷载作用边缘离断面中轴线的最短距离为2m。并且，均布荷载为一个呈正弦函数规律变化的动荷载，其幅值为（100+学号后两位数字）kpa，频率由学生自己确定。假设在均布荷载作用期间，土的基本物理特性参数均为恒量，不随时间改变。

基本要求：学生采用现成的大型商业软件，如ABAQUS,ANSYS,PLAXIS以及FLAC 等。针对上述案例，建立自己的数值模型，模型可以是2D的，也可以是3D的。模型的本构关系要求在MC模型、Cam-clay模型、修正的Cam-clay模型。DP模型中进行选择。通过数值模拟，学生可以对路基断面中轴线上各点进行动力响应分析。要求在中轴线上选取5个关键点，分别为路面以下2m，4m，6m，8m，10m 处。提取这五个点的位移响应曲线、速度响应曲线、应力-应变关系曲线。最后，提交一份完整的计算书，同时附上数值模拟计算的源程序。

解：

1、假设高速公路为单向的三车道第一层砂土刚好在地平线上，第二层为黏土刚

好地平线以下，边坡为1：1.5。断面为12m，则边坡边缘刚好为24m。

2、利用圣维南原理，取x与y方向各为路基断面的6倍，因此计算模型可简化

为平面应变问题。假设断面承受的外力不随z轴变化，位移和应变都发生在自身平面内。

3、因为荷载随时间的变化而变化，可采用瞬态分析。

土体单元选取plane183单元来模拟。Plane83为2D八节点单元，用于模拟实体

结构。该单元可以作为平面单元（平面应力、平面应变和广义平面应变），每个节点具有2个自由度，即沿着节点坐标系x和y方向的平动位移。可利用混合模式模拟不可压缩的弹塑性行为或完全不可压缩材料的超弹性行为。其特点能够有效的反应土体的各种性状。

4、结构各层材料参数

类别厚度（m）弹性模量E

（KPa）泊松比v 密度D 粘聚力c

内摩擦

角膨胀角

砂土 4 1.39E+07 0.3 1600 0 33 20 黏土 4 5.39E+05 0.35 1800 27 37 29 砂土 4 2.39E+07 0.3 1600 0 35 20 硬质黏土（假设）48 4.00E+07 0.2 2100 27.4 27 8

139000sin(/2)

force t

=⨯

5、施加的正弦荷载

6、建立的模型

7、网格划分

8、竖直方向上的位移云图

9、中轴线-2、-4、-6、-8、-10处的位移响应曲线

结论：由竖向位移响应曲线可以得到，路基断面在受正弦荷载情况下，整体由于受到重力影响存在重力加速度，刚开始受力的时候，力逐渐增大又逐渐减小，因此力为正值时是所加速度逐渐增大又逐渐减小的向下加速运动，位移曲线为向上凸的抛物线，当力为负值时，力所对应的加速度减去重力加速度，整体加速度可能向上也可能向下，但由于正值力对应下，使整体具有向下的速度和惯性力，由结果曲线可知，加速度任然为向下，只是加速度逐渐减小，因此总体位移向下。路基越下面，位移越大。

10、中轴线-2、-4、-6、-8、-10处的速度响应曲线

结论：由达朗贝尔定理，能量原理可知，力做的功转化成动能、势能和应变能，整体做加速度逐渐增大再减小，再减小再增大的循环变化，因此，总体的速度方向不变，为向下的增大的，只是向下增大的快慢，受外力f的影响。

11、中轴线-2、-4、-6、-8、-10m处的应变曲线

结论：由弹塑性原理可知应变为位移的导数：y y y u u t t

v y

t y y

ε∂∂∂∂=

=•

=•∂∂∂∂可知应变随着速度和时间，还有竖向位移的变化而变化。由图可知应变随着深度的增加而增加，证明了土力学，等值线下凹的现象 12、

中轴线-2、-4、-6、-8、-10m 处的应力曲线

应力曲线结果显示为0，具体还没弄明白。

13、中轴线-2、-4、-6、-8、-10m处的应力应变曲线

结论：变形在弹性范围内，没有塑性变形。

14、命令流

FINISH

/CLEAR

/CONFIG,NRES,1000 !设置子步结果限值

/nopr !菜单过滤设置

keyw,pr_struc,1 !保留结构分析部分菜单

/units,SI !单位设置

/FILNAME, The subgrade cross-section analysis

/TITLE, The subgrade cross-section analysis

/PREP7

ET,1,plane183 !定义单元类型Keyopt,1,3,2 !求解类型为平面应变问题!土基4m砂土

MP,EX,1,1.39e10

MP,PRXY,1,0.3

MP,DENS,1,1600

tb,dp,1

tbdata,,0,33,20

!底基层4m的黏土

MP,EX,2,5.39E8

MP,PRXY,2,0.35

MP,DENS,2,18000

tb,dp,2

tbdata,,27400,37,29