邓肯-张模型公式推导高土
高等土力学-土的本构关系

第二章 土的本构关系
2.1 概述
仁者乐山 智者乐水
体积力 面力 静(动) 力平衡
应力
本构方程
位移
几何 相容
应变
本构关系在应力应变分析中的作用
第二章 土的本构关系
2.1 概述
传统土力 学分析方法
变形问题 (地基沉降量)
稳定问题 (边坡稳定性)
仁者乐山 智者乐水
• 弹性理论计算应力 • 压缩试验测定变形参数 • 弹性理论+经验公式计算变形
第二章 土的本构关系
本章内容提要
土的变形特性
土的非线性弹性模型 • 邓肯张EB和E模型
土的弹塑性模型
• 剑桥模型(CamClay) • Lade-Duncan模型 • 清华弹塑性模型 • 沈珠江双屈服面模型
第二章 土的本构关系
p108页 – 109页 第 14,18,19,33题
第二章 土的本构关系
O
B
A
v
3
1
平均应力 p
等向 固结线
B
土样总
A
体剪胀区
O
近似弹性
剪切起点
剪缩剪胀 分界点
V0
体应变 v
第二章 土的本构关系
2.3 土的应力变形特性
弹塑性、滞回圈、卸载体缩
400 q
仁者乐山 智者乐水
200
滞回圈
0
2
1
v
p
4
6
e
卸载 体缩
8 1 (%)
单调与循环加载的三轴试验曲线 (承德中密砂)
第二章 土的本构关系
土力学中应力符号规定
第二章 土的本构关系
2.2 应力和应变 – 应力
高等土力学(李广信)2-5章部分习题答案

2-1.什么叫材料的本构关系?在上述的本构关系中,土的强度和应力-应变有什么联系? 答:材料的本构关系是反映材料的力学性质的数学表达式,表现形式一般为应力-应变-强度-时间的关系,也成为本构定律,本构方程。
土的强度是土受力变形发展的一个阶段,即在微小的应力增量作用下,土单元会发生无限大或不可控制的应变增量,它实际上是土的本构关系的一个组成部分。
2-7什么是加工硬化?什么是加工软化?请绘出他们的典型的应力应变关系曲线。
答:加工硬化也称应变硬化,是指材料的应力随应变增加而增加,弹增加速率越来越慢,最后趋于稳定。
加工软化也称应变软化,指材料的应力在开始时随着应变增加而增加,达到一个峰值后,应力随应变增加而下降,最后也趋于稳定。
加工硬化与加工软化的应力应变关系曲线如右图。
2-8什么的是土的压硬性?什么是土的剪胀性?答:土的变形模量随着围压提高而提高的现象,称为土的压硬性。
土的剪胀性指土体在剪切时产生体积膨胀或收缩的特性。
2-9简述土的应力应变关系的特性及其影响因素。
答:土是岩石风化形成的碎散矿物颗粒的集合体,通常是固、液、气三相体。
其应力应变关系十分复杂,主要特性有非线性,弹塑性,剪胀性及各向异性。
主要的影响因素是应力水平,应力路径和应力历史。
2-10定性画出在高围压(MPa 303<σ)和低围压(KPa 1003=σ)下密砂三轴试验的v εεσσ--)(131-应力应变关系曲线。
答:如右图。
横坐标为1ε,竖坐标正半轴为)(31σσ-,竖坐标负半轴为v ε。
2-13粘土和砂土的各向异性是由于什么原因?什么是诱发各向异性?答:粘土和砂土的各向异性是由于其在沉积过程中,长宽比大于1的针、片、棒状颗粒在重力作用下倾向于长边沿水平方向排列而处于稳定的状态。
同时在随后的固结过程中,上覆土体重力产生的竖向应力与水平土压力大小不等,这种不等向固结也造成了土的各向异性。
诱发各向异性是指土颗粒受到一定的应力发生应变后,其空间位置将发生变化,从而造成土的空间结构的改变,这种结构的改变将影响土进一步加载的应力应变关系,并且使之不同于初始加载时的应力应变关系。
黄土邓肯-张模型有限元计算参数的试验

高 江 平1#李 ! 芳6
#1F长安大学 特殊地区公路工程教育部重点实验室%陕西 西安 713324& 6F长安大学 建筑工程学院%陕西 西安 713321$
摘 ! 要 !为 配 合 有 限 元 计 算 分 析 的 需 要 #从 工 地 现 场 取 样 #进 行 了 黄 土 的 室 内 三 轴 试 验 #并 对 试 验 结 果进行了分析#得到了邓肯 张非线性模型有限元计算的5个参数$黄土的有限元参数随着试样 的 压实度和试验时围压的不同而变化$在含水量保持不变的情况下#邓肯 张 C"# 模型的 5 个参数 中#参数 X%*%<.%M%’%8 均随压实度的提高而增大&D%O 随 X 的变化规律不明 显$ 在 含 水 量 和 压 实度保持不变的情况下进行三轴试验#土样的!$1d$8".%!$1d$8"&EP%C4%<. 均 随$8 的 增 大 而 增 大# #4 随$8 的增大而减小#O 随$8 变化的规律不明显$ 关键词!道路工程&邓肯 张模型&黄土&参数&试验 中 图 分 类 号 !S414F1&[S48! ! ! 文 献 标 识 码 !B! ! !
) T1Y3"强度恰好发挥*)&1Y3"则土体发生塑流#
) 值 愈 大"塑 流 变 形 愈 大"但 该 值 不 超 过 &$1 U
$8’&EP(&$1 U$8’.#
可 见"邓肯 张C (#模型具有5个参数"即)X!
*!<. !M!’!8!D!O#
4!试验简介
土 样 取 自 甘 肃 某 工 地 "按 现 行 +公 路 土 工 试 验 规 程,$8%规定的方法"经颗粒分析!液塑限试验"确 定 该
不同应力路径下的邓肯_张模型模量公式

摘
要:实际工程中,土体可能处于轴向卸载、侧向加载、侧向卸载应力路径下,而邓肯–张模型是在轴向加载条件
下得到的,这限制了它的适用范围。本文模拟邓肯—张模型思路推导了不同应力路径下的切线模量公式,从而使模量 公式系列化,扩大了邓肯–张模型适用范围;根据不同应力路径试验的结果,发现公式推导中使用的假设可以得到验 证。 关键词:应力路径;邓肯–张模型;模量公式;三轴试验 中图分类号:TU431;X705 yindeshun @。 文献标识码:A 文章编号:1000–4548(2007)09–1380–06 作者简介:殷德顺(1972– ),男,山东人,讲师,主要从事岩土力学基础理论和岩土数值分析方面的研究。E-mail:
与邓肯–张模型类似,同样假设初始切线模量 Ei (Ei 为 2(σ rc − σ r ) – ε a − 2ε r 曲线上原点处的斜率) 随 轴 向 固 结 压 力 σ ac 而 变 化 , 且 在 双 对 数 纸 上 点 绘 lg( Ei / pa ) 和 lg(σ ac / pa ) 的关系,可用直线来拟合,设 直线的截距为 k、斜率为 n。
0
引
言
20世纪以来,伴随着土木水利工程的发展,高层 建筑、地下工程以及水利设施的数量迅速增多,这些 工程的建造使承担基础作用的土体处于不同荷载的作 用下。 很多学者对于土体应力–应变关系受应力路径影 响的问题进行了研究,而且这一方面的研究积累也非 常多。Lambe[1]首先提出了应力路径方法。曾国熙等[2] 的研究发现,软黏土的应力–应变关系不但具有非线 性特点,而且受应力路径的影响,正常固结饱和黏土 的应力–应变曲线可以按双曲线拟合,并可归一化; 他还将这个结论应用于软黏土的深基坑开挖非线性有 限元中。刘国彬、侯学渊[3]对上海地区有代表性软土 的卸荷模量的变化规律研究发现,软土的应力–应变 关系与应力路径密切相关,其卸荷应力–应变关系为 双曲线形式,并得到了切线卸荷模量的计算公式。应 宏伟等[4]对杭州、 上海饱和软黏土采用排水分析研究, 得到了同时考虑应力路径和应力历史影响的、用有效
有限元法的变分原理及其在土石坝设计中的应用

有限元法的变分原理及其在土石坝设计中的应用有限元法是采用直接法计算变分问题的重要方法,在土木工程计算领域的分析软件如ANSYS、Workbench、Autobank等均以变分法为理论基础。
本文将就有限元法的变分原理作一简单梳理,并采用Autobank软件建模分析某土石坝的渗流场及应力变形,计算结果表明大坝应力变形符合工程实际,计算分析对大坝设计工作起到了指导作用。
标签:有限元;变分法;Autobank;土石坝设计;应力变形分析引言随着坝工技术的发展,土石坝建设高度越来越高,其应力和变形计算越来越关系到大坝安全。
因此,结构计算分析将会在土石坝的设计和科学研究中发挥越来越重要的作用。
有限元法的理论基础为变分法,变分法历史悠久,是近代发展起来的一门重要数学分支,在工程技术及科学研究中有着广泛的应用。
变分法起源于泛函的极值问题,其关键定理是欧拉-拉格朗日方程。
Autobank软件应力变形分析模块是以变分法为理论基础开发的一款有限元分析软件,提供线弹性模型、非线性模型(如邓肯E-B、E-μ模型)等,在水利工程设计中有着广泛的应用。
1、有限元法简介目前在水利工程结构分析领域常用的数值计算方法有:有限差分法FDM、有限元法FEM、边界元法BEM、离散元法DEM等,其中有限元法是应用最广泛的方法。
有限元法是以变分原理为基础发展起来的,是一种高效的数值计算方法。
工程计算和科学研究领域,常常需要求解各类常微分方程(组)、偏微分方程(组),而许多微分方程(组)的解析解很难得到,甚至无法求出。
使用有限元法将微分方程离散化后,编制计算机程序辅助求解,是一种可行且高效的方法。
2、有限元法的变分原理2.1 泛函及其极值设有泛函的极值问题:研究泛函在某函数类中的极值问题即变分问题,例如最小曲面问题、悬链线问题、边坡稳定最小安全系数的滑弧问题、重力坝的最优断面问题等。
研究泛函极值的方法即变分法。
直接法是求解泛函极值的近似方法,对于无法求解解析解的变分问题及工程计算,有着及其重要的作用。
邓肯-张模型参数变化对计算结果的影响

!"#$% ?
# 值的变化对最大位移及应力水平 " 的影响 " 1 变化率 8 ; > ?B@ > =B< = E =B< E ?BF " , 8 02 > &A B &C > &A B =? > &A B @< > &A B =? > &< B G? " , 变化率 8 ; > ?BC > <B@ = > <B@ > AB? " 变化率 8 ; E ?B& E =BD = > <BD > =BD
邓肯 ! 张模型在国内外广泛使用近 )’ 年, 大量的 试验成果表明, 由于取样制样、 试验仪器、 试验方法与 过程、 试验人员操作熟练程度、 整理分析资料等诸多因 素, 使其 * 个参数变化较大,! 值可成倍甚至成量级相
[# B )] , 用于计算所得结果的差别也较大。邓肯等人 差
曾对该模型的参数作了初步讨论, 对几种不同类型的
!
前
言
!
偏应力 ("# ! ") ) 不太大时, 就能达到较高应力水平 ’ , ( 从而使 " > 降低, 变形增大。 ’ F( G "# ! ") ) "# ! ") ) ?, 黏聚力 # (#) 黏聚力 # 的增减对水平位移 ’ 4 、 垂直位移 ’ H 、 应 力水平 ’ 的影响见图 # B ) 及表 #。
[-] 土给出了参数的范围, 并编制了图表 。这些图表变
化范围较大, 不同的取值对计算结果的影响没作进一 步讨论。 本文依据文献 [(] 分别增减 * 个参数, 用有限元法 考察对一个均质土坝的位移和应力水平的影响程度。 考察某一参数时, 其余 / 个参数不变, 即保持试验取值。
土体邓肯—张非线性弹性模型参数反演分析
土体邓肯—张非线性弹性模型参数反演分析《土体邓肯张非线性弹性模型参数反演分析》是一项重要的科学研究,在这项研究中,研究人员将利用非线性弹性模型的参数反演分析技术,来研究土体的张拉性能。
对于土体的张拉性能,其抗拉强度、塑性性能和力学特征都具有较强的不确定性。
为了克服这种不确定性,非线性弹性模型的参数反演分析技术提供了一种有效的手段,可以在试验过程中获取到准确的参数模型,并进一步深入研究土体的张拉性能。
一、非线性弹性模型的概述非线性弹性模型是一种利用实验测量的参数,结合物理模型的理论模型,用来研究弹性反应的内部结构和行为的模型。
它可以用来描述弹性物质的非线性物性,以及土体张拉时的力学特性。
典型的非线性弹性模型包括邓肯张模型、兰氏莫尔斯模型、HarrisYarwood模型、RiceVangenuchten模型等,其中邓肯张模型是最常用也是最具代表性的模型。
邓肯张模型是一种非线性弹性模型,其物理模型以弹性超塑性为基础,以土体张拉过程中的变形量、压实率及应力量的关系为参数,描述土体张拉过程中的行为特性。
二、参数反演分析技术参数反演分析技术是一种基于回归模型的参数估计的数学方法,可以利用与实验数据相关的模型参数,经过迭代优化,最终得到最佳匹配的参数模型。
非线性弹性模型的参数反演分析技术由实验中获取的非线性参数和迭代优化模型结合而成,可以对非线性弹性模型的参数进行更为准确的反演分析,实现对土体张拉特性和参数之间联系的准确描述。
三、土体邓肯张非线性弹性模型参数反演分析非线性弹性模型的参数反演分析,可以通过试验确定土体张拉过程中的力学特性,并反演出最佳的张拉参数,以深入了解土体的张拉性能。
本文以著名的邓肯张非线性弹性模型为例,通过实验过程,获取相关参数,建立非线性弹性模型,并以最小二乘法、Simplex法等为基础,进行参数反演分析。
经过迭代优化,最终获取到准确的非线性弹性模型参数,从而对土体的张拉性能更为准确的描述。
邓肯-张模型开发及其在面板坝计算中的应用
邓肯-张模型开发及其在面板坝计算
中的应用
邓肯-张模型是一种用于计算面板坝的模型,它是由美国土木
工程师邓肯和张在20世纪50年代提出的。
该模型基于坝面的水力学原理,以及坝面的结构特性,建立了一个简单的模型,用于计算面板坝的水力特性。
邓肯-张模型的基本原理是,坝面上的水流可以分解为两个部分:一个是水流的横向分布,另一个是水流的纵向分布。
横向分布的水流可以用一个简单的模型来描述,即水流的流量与坝面的宽度成正比,而纵向分布的水流可以用一个简单的模型来描述,即水流的流量与坝面的高度成正比。
根据邓肯-张模型,可以计算出面板坝的水力特性,包括水流
的流量、水流的速度、水流的压力等。
此外,该模型还可以用来计算面板坝的抗滑性能,以及面板坝的稳定性。
邓肯-张模型在面板坝计算中的应用非常广泛,它可以用来计
算面板坝的水力特性,以及面板坝的抗滑性能和稳定性。
此外,该模型还可以用来计算面板坝的抗滑性能,以及面板坝的稳定性。
此外,该模型还可以用来计算面板坝的抗滑性能,以及面板坝的稳定性。
总之,邓肯-张模型是一种简单而有效的模型,可以用来计算
面板坝的水力特性、抗滑性能和稳定性,在面板坝计算中有着广泛的应用。
高等土力学土的本构关系PPT课件
(3)应力历史
第32页/共86页
2.4 土的弹性模型
概述
1、线弹性模型
x y
1
E 1
E
x y
y
z
z
x
z
1 E
x
y
z
xy
21
E
xy
yz
21
E
yz
zx
21
E
zx
第33页/共86页
亦可表示成:
p Kv q 3G
式中:
K
E
31 2v
G
E
21
a 1 3
a
b a
常规三轴试验, d 2 d3 0
Et
d1 3
d1
a
a b12
试验起点,1 0 ,Et Ei
Ei
1 a
第39页/共86页
如果 1
1
3 ult
1 b
所以:
a 1 Ei
b
1
1
3
ult
定义破坏比 Rf
Rf
1 1
3f 3 ult
第40页/共86页
则有:
b
1
1 3 ult
x
y
xzy
yz
zx
T x , y , z , xy, yz, zx
第9页/共86页
2. 应力张量的主应力和应力不变量
l cos
m
cos
n cos
第10页/共86页
x0 y0 z 0
x
xyl
l
y
yxm m
zx zy
n n
0 0
zx l
yzm
土力学问题: 变形问题:弹性理论 强度(稳定)问题:极限平衡分析
二维非线性邓肯-张模型分析软基沉降变形规律
利用物 理 、 位移 、 几 何 以及空 隙压力 连续 性方程 反应 应 力 和位移 的关 联 , 并 且考 虑 了泊 松 比和 时 间 对路 基 固结 沉降 的影 响 , 在 理论上 是 比较完善 的 , 但 由于 其 复杂 性很 难给 出解 析解 , 目前 都 多采 用 数 值解 法 进 行计 算 。
坡 部分 对应 的软 基表 面上 , 而深度 较 大 (>1 5 m 以上 ) 时, 路堤 中心点 处的沉 降 变形最 大。 关键 词 : 软基 ; 沉 降 ;邓肯 一 张模 型 中图分类号 : U 4 1 6 . 1 文献标 识码 : A
0 前 言
目前 , 我 国公 路软 土 路 堤… 沉 降 分 析 主要 还 是
中填方 的施 工速 度 , 路堤高 2 m时 , 为 1 级加载 ; 路
3 4
湖 南
交 通 科 技
4 l卷
为讨论 不 同路 堤高 度 、 软 基厚度 的情况 下 , 软 基
1 0 9 5 0・ 9 O. 8 5
0 8
的沉降变形规律 , 分别计算 了路堤顶面宽度为 3 0 m
{ s ) = [ 】 { ) +[ 露】 { )
( 3 )
1 理 论 基 础关 系采 用 Ha r d i n—D me v i c h双 曲线 型本 构关 系 , 见式 ( 1 ) :
式中: [ ] 为通常的 单元刚 度矩阵; { 6 ) 为单元的
时, 软基厚 度 分别 为 6 、 1 2 、 1 5 、 2 0 m, 路堤 高 度为 2 、
5 、 7 I l l 时路 堤沉 降变 形 规 律 , 路 堤 的边 坡 坡 比均 为
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
邓肯-张模型是一个非线性本构模型,既然是一个本构模型,可想而之他反应的是应力与应变之间的关系。
说它是非线性的,那么反映应力应变关系的模量就不是一个常数E那么简单。
在介绍该模型之前,先要介绍一个概念,就是反映非线性关系的增量广义胡克定律: 1123()t
t t
v d d d d E E σεσσ=
-+ (1) 1963年,康纳(Kondner )根据大量土的三
轴试验的应力应变关系曲线,提出可以用双曲线拟合出一般土的三轴试验13()~a σσε-曲线,即:
13a
a
a b εσσε-=
+ (2)
其中,a 、b 为试验常数。
对于常规三轴压缩试验,1a εε=。
邓肯等人根据这一双曲线应力应变关系提出了一种目前被广泛的增量弹性模型,一般被称为邓肯-张(Duncan-Chang )模型。
在常规三轴压缩试验中,13a
a
a b εσσε-=+可以写成:
1113
a b εεσσ=+- (3)
将常规三轴压缩试验的结果按
11
13
~εεσσ-的关系进行整理,则二者近似成线性关系(见图
1)。
其中,a 为直线的截距;b 为直线的斜率。
在常规三轴压缩试验中,由于
230d d σσ==,所以切线模量为
ε1
/(σ1
-σ3
)
1
-σ3
)ult
图1
1113
~εεσσ-线性关系图
132
11()()
t d a
E d a b σσεε-=
=+ (4) 在试验的起始点,10ε=,t i E E =,则:
1i E a
=
,这表明a 表示的是在这个试验中的起始变形模量E i 的倒数。
如果1ε→∞,则: 131
()ult b
σσ-=
(5) 由此可以看出b 代表的是双曲线的渐近线所对应的极限偏差应力13()ult σσ-的倒数。
在土的试样中,如果应力应变曲线近似于双曲线关系,则往往是根据一定的应变值(如
115%ε=)来确定土的强度13()f σσ-,而不可
能在试验中使1ε无限大,求取13()ult σσ-;对于有峰值点的情况,取1313()()f σσσσ-=-峰,
这样1313()()f σσσσ--ult <。
定义破坏比R f 为:
1313()()f f R σσσσ-=-ult
(6)
而
13131
()()f f
R b σσσσ==
--ult (7) 将上式与1
i E a
=
代入 132
11()()
t d a
E d a b σσεε-==+ (8) 得到:
塑性力学读书报告
2
113111()t f
i i E R E E εσσ⎡⎤
⎢⎥⎢⎥=⎢⎥+⎢⎥
-⎣⎦
(9)
该式表示为应变1ε的函数,使用时不方便,可将t E 表示为应力的函数形式。
由式
1113
a b εεσσ=+- (10)
可以得到13113()
1()
a b σσεσσ-=
--,将该式代入
132
11()()t d a
E d a b σσεε-=
=+得到
22
1313131()1
[][]
1()
1()
t a
E ab a a b b σσσσσσ=
=
-+
----将式13131
()()f f
R b σσσσ=
=
--ult 和1i E a =代入上式得到:
2
1313[1]()t i i
f
E E R σσσσ-=-- (11)
根据莫尔-库仑强度准则,有 3132cos 2sin ()1sin f c ϕσϕ
σσϕ
+-=
- (12)
如果绘制lg(/)i a E P 与3lg(/)a P σ的关系图,可以发现二者近似呈直线关系,所以得式:
3
(
)n i a a
E KP P σ= (13)
其中,a P 为大气压(a P =101.4k a P ),量纲与3σ相同;K 、n 为试验常数,分别代表
lg(/)i a E P 与3lg(/)a P σ直线的截距和斜率。
将
3132cos 2sin ()1sin f c ϕσϕ
σσϕ
+-=
- (14)
和3
(
)n i a a
E KP P σ=
代入2
1313[1]()t i i
f
E E R σσσσ-=--,则得到:
1323
3()(1sin )(
)[1]2cos 2sin f n
t a a
R E KP P c σσϕσϕσϕ
--=-
+(15)
可见切线变形模量中包括5个材料常数K 、n 、ϕ、c 、f R 。
2 切线泊松比(poisson's ratio)
Duncan 等人根据一些试验资料,假定在常规三轴压缩试验中轴向应变1ε与侧向应变3ε-之间也存在双曲线关系 ()
3
13f D εεε-=
+- (16)
或者 ()3
331f D f D εεεε-
=+-=- (17)
从上式,试验得到的31
εε-与3ε-的关系可近似为直线关系,从而确定截距f 与斜率D 。
从式上式可见当30ε-→时,
310i f V εε⎛⎫-→== ⎪⎝⎭
i V 即为初始泊松比。
见图 (a)。
D 为
13εε-:关系渐近线的倒数,见图(b)。
试验表明
土的切线
王吉亮(2006631011)
泊松比i V 是与试验的围压3σ有关的。
它们画在单对数坐标中,可假设是一条直线,见图(c),这样:
()3f=G-Flg Pa i v σ= (18)
G 、F 为试验常数,其确定见图(c)。
将(16)式微分:
()()()
11322
111111i
t D f D f d V v d D D εεεεεε-+-=
==-- (19) 将1ε表达式代入式(19),则得到
()
()()()32
131333lg 11sin 12cos 2sin t n
f G F Pa v D R KPa Pa c σσσσσϕσϕσϕ-=
⎧⎫
⎪⎪
-⎪⎪
-⎨⎬--⎡⎤⎛⎫⎪⎪-⎢⎥ ⎪⎪⎪+⎝⎭⎣⎦⎩⎭
在切线泊松比式中又引入G 、F 、D 三个材料常数。
加上t E 中五个常数,共有八个常数。
根据弹性理论,
00.5t v <<。