循环载荷下非凝聚土本构方程的数值计算.pdfVIP

!三矍查堂!堕型 !!!!生 ：!!!： 循环载荷下非凝聚土本 构方程的数值计算 章根德 韩贝传 (中国科学院力学研究所，北京l∞o鲫)(中国科学院地质研究所．北京l∞029) 摘要循环荷藏下非凝聚土弹-塑性本构方程的鼓值计算包含了许多离散化．这将导 致偏离屈服面的河题，弹塑性加截增量所预言的应力状态不蔷在属晨面上．这样，偏 离将积累起来，所以在每一步加载增量后进行修正，以确保应力仍落在屈服面上就显 得十分重要丁．本文考察了五种计算应力点信移的方j圭．将数值结果和实验资料进行 丁比较．发现沿弹塑性流动方向反向投影的惨正方{蠹给逝比较精确的预言，如用其他 方法对屈服面的偏移进行修正将产生较明显的误差．所以。这种修正方法可以在弹． 塑性分析的数值计算中广泛地使用． 关键词本构方程。数值计算，修正方法 1引言 解决大型的、复杂的土-结构相互作用问题，如水坝、核电站、海洋工程结构受到地 震、爆炸冲击、应力波作用的问题，都需要先进的计算模型．饲如典型的海洋平台，自重就 达105吨，而引起的倾覆力矩可达106吨米．在这样的载荷条件下，如再借助于常规的经 验公式来设计将会出现许多不可估量的问题．于是，常常采用先进的数值计算方法，如有 限元方法等，来代替常规的承载力公式进行设计．面这些效值方法中，关键就是对土介质 本构响应的近似模拟．这就导致了人们对土本构模型研究的极大关注，使最近20年来土本 构模型的研究取得了较大的进展． 在岩土介质的弹一塑性分析中，常采用小步长进行数值积分．弹塑性本构矩阵【B，J实 际上是随着应力的变化而连续变化的。但在计算中，对每一小步长总假设【G，】是个常数． 所以，尽管采用小步长，但加载增量结束后，应力状态往往会不在屈服面上．特别是当屈 服面随着载荷的增加而运动时。这种偏离就会更加明显．这就要求我们在计算的每一步都 要通过修正将应力拉回到屈服面上来． 可以采用不同的方法来修正{“q．经过比较表明沿着塑性流动方向反向投影进行修正 有较高的精度．所以本文所给出的弹塑性应力一应变关系的数值计算采用了沿着塑性流动 方向反向投影修正的方法． 2循环荷载下非凝聚土的本构方程 文献【61提出了一个描述循环载荷下非凝聚土本构特性的模型．分析常规土力学试验 所得到的应力．应变曲线，如图1所示．可以发现，在加载过程中，当土介质达到某种应 力状态后。应力一应变曲线呈现出非线性。即人们所常说的达到了屈服．随后便是继续屈 章根德．男．1942年lO月出生，研究员 国家自然科学基金资助项耳． ·338- 《工程力学，增刊 1999年 变硬化过程．所以，在加载过程中，土介质的屈服现象都可以与屈服面发生联系．屈服面可 以表示为： F=J2D—02并(1一nl开一2)(1一n3耳)一}(1) 式中，oI’铆，03为材料常数， J1=以l， J2。=；＆％ ^。=；最j岛·t 重新加戴 sT=jInl啦 对于卸载与重新加载。实验表明 许多种土的应力·应变性状几乎是线 弹性的，并且与再次加载过程所处的 应力水平无关【7_…．我们实验所采用 田l加载-卸曩过程的近戗裹示 的福建标准砂、常规三轴压缩实验所 得的结果也是这样。如图5所示． 文献【6】提出的模拟循环荷载下非凝聚土的本构模型具有传统的增量理论的结构．在三 维笛卡尔坐标系中，模型的数学表达可以写成： dE-，=dE0+出3 (2) as；=壶如帆t+嘉(aa玎一扣帆t) (3)