没有合适的资源?快使用搜索试试~ 我知道了~
温馨提示
建立了能考虑塑性变形、温度、饱和度和孔隙比等对材料特性影响、土-水特性曲线并采用扩展的 Mo-hr-Coulomb屈服准则、以屈服接近度为安全指标的尾矿坝流固耦合有限元分析模型。分别运用上述方法和强度折减系数法对龙都尾矿坝干滩面为100 m时的饱和-非饱和渗流场和变形场的耦合响应特性进行了模拟,比较表明两者的结论吻合。但本文所采用的方法仅需一次计算就能获得坝体变形、浸润线等重要指标,以及考虑各种影响因素的屈服接近度分布,据此可直观全面地了解不同工况下坝体安全的全场信息,同时大大简化了分析过程和减少了计算工
资源推荐
资源详情
资源评论
第
42
卷
第
3
期
2010
年
6
月
西 安 建 筑 科 技 大 学 学 报
(自然科学版)
JXi′anUniv.ofArch.& Tech.
(
NaturalScienceEdition
)
Vol.42
No.3
Jun.2010
基于屈服接近度概念的尾矿坝静力稳定性分析
及其与强度折减系数法的对比
曹林卫
1
,彭向和
1
,杨春和
1
,
2
(
1.
重庆大学西南资源开发及环境灾害控制工程教育部重点实验室,重庆
400044
;
2.
中国科学院武汉岩土力学研究所,湖北 武汉
430071
)
摘
要:建立了能考虑塑性变形、温度、饱和度和孔隙比等对材料特性影响、土
水特性曲线并采用扩展的
Mo
hrCoulomb
屈服准则、以屈服接近度为安全指标的尾矿坝流固耦合有限元分析模型
.
分别运用上述方法和强
度折减系数法对龙都尾矿坝干滩面为
100m
时的饱和
非饱和渗流场和变形场的耦合响应特性进行了模拟,
比较表明两者的结论吻合
.
但本文所采用的方法仅需一次计算就能获得坝体变形、浸润线等重要指标,以及考
虑各种影响因素的屈服接近度分布,据此可直观全面地了解不同工况下坝体安全的全场信息,同 时 大 大 简 化
了分析过程和减少了计算工作量和人为因素的影响,有利于对尾矿坝的稳定性进行客观的评价
.
关键词:尾矿坝;稳定性;有限元分析;屈服接近度;强度折减系数法
中图分类号:
TV64
文献标识码:
A
文章编号:
10067930
(
2010
)
03040708
.
尾矿坝是以尾矿为主要材料建成的一种特殊的构筑物,坝坡经受长期多循环水力充填以及固结沉
降作用,坝体材料呈层状特征,其多相性与不均匀性使尾矿坝的分析较通常的边坡或土石坝更为复杂
.
尾矿库既是矿山生产的重要设施,又是重要的危险源
.
尤其是许多尾矿库的下游有村落、工厂学校、交通
设施等,其中相当部分属于一级保护设施
.
因此,尾矿库的稳定性不仅关系矿山自身的安全,而且涉及到
周边及其下游居民的生命财产安危
.
所以尾矿坝的稳定性分析十分重要
.
目前尾矿坝稳定性分析最常用的是极限平衡法、数值分析法,其 中极限 平衡法 原理简 单、且能 直接
提供有关整个坡体安全性的评估,在设计验算中应用较多,但此法以假设潜 在滑动 面为基 础,认 为 滑坡
面以上土体为理想刚塑性体,容易造成对实际情形的偏离且无法直观显示坝内各部分所处状态;随着计
算机和数值计算方法的迅速发展,数值分析法特别是有限元法的运用日益广泛,它采用弹塑性应力应变
关系进行变形和稳定分析,并提供坡体在 载荷作 用 下应力
应变状 态以及 相关安 全性指 标 的分 布,便于
及时地发现薄弱环节并进行维护
.
采用常规的有限元法进行尾矿坝稳定性分析不能提供坝体的整体安全系数
.
为此,近年来人们开展
了相关研究并取得一些成果
.
强度折减法
[
1
]
以有限元数值计算不收 敛或者 滑 面上节 点塑性 应变和 位移
发生突变作为土体破坏的依据,通过试算,得出坝体的安全系数
.
此法存在两点不足:(
1
)需多次试算,计
算成本较高;(
2
)失稳判据的选择尚未统一
[
2
]
.
此外,折减系数法改变 了材料 的承载 特 性,可能 造成坝 体
中应力分布的失真
.
对于岩土结构而言,合理的计算应该基于能考虑各种主要影响因素的应力与变形分布、能考虑各种
主要影响因素的本构模型和失效判据、简单易行的计算思路以及能直观反映材料和结构安全程度的指
标
.
本文采用有限元软件
ABAQUS
/
Standard
对尾矿坝进行饱和与非饱和条件下渗流场和变形场直接
耦合分析,采用扩展的
MohrCoulomb
屈服准则,并考虑饱和度 变化、等效塑性应变以及其他由用户预
收稿日期:
20081005
修改稿日期 :
20100418
基金项目:国家自然科学基金资助项目(
10872220
)
作者简介:曹林卫(
1981
),男,江苏如东人,博士研究 生,主要从事边坡稳定性方向的研究
.
先定义场变量对饱和非饱和土抗剪强度的影响;引入屈服接近度
η
[
3
]
作为材料屈服程度度量指标,结果
可直观全面地反映不同工况下坝体安全的全场信息
.
此方法可大大简化分析过程、减少计算工作量和避
免人为因素的影响,具有重要的应用价值
.
1
渗流场与变形场耦合的分析模型
1.1
基本假定
基于下列假设建立渗流场与变形场耦合的尾矿坝分析模型:
(
1
)尾矿是由固体尾矿颗粒与孔隙(水和空 气)所组 成的多 相介质,如果认为孔隙气压力足够小而
忽略其影响,则其力学特性决定于水和尾矿颗粒的特性及相互作用;
(
2
)假设尾矿材料初始各向同性,且颗粒和水都不可压缩,尾矿所表现出的体积压缩性取决于孔隙
比的变化;
(
3
)渗流满足达西定律
v
w
=-
k
sn
·
H
x
(
1
)
式中:
s
为饱和度;
n
、
k
和
H
为分别为孔隙率、渗透系数和压力水头,分别表示为
n
=
e
1+
e
,
k
=
y
(
s
)
F
(
e
),
H
=
z
+
u
w
g
ρ
w
(
2
)
式中:
e
为孔隙比;
v
w
为水相对于土体 骨架运动 的平均 速度(渗流 速度);
z
为某点 相对于 参考 面的 位 置
标高(位置水头)
.
u
w
和
ρ
w
分别为孔隙水压力和水的密度;
g
为重力加速度(沿
z
方向);
y
(
s
)和
F
(
e
)分别
用于描述
s
和
e
对渗透系数的影响,在分析中取
y
(
s
)
=
s
3
,
F
(
e
)
=
k
0
I
,
k
0
为不同尾矿层在饱和条件下的
初始渗透系数
.
计算中函数的定义通过用户子程序
USDFLD
实现;
(
4
)当孔隙水压力
u
w
≥
0
时,尾矿饱和度
s
=1.0
,否则其饱和度通过土
水特性曲线
u
w
~
s
确定
.
1.2
控制方程
以多孔介质理论为基础,采用
Galerkin
法,以节点位移和孔隙水压力为独力变量进行空间离散,得
到耦合变形场和渗流场的平衡方程
[
4
5
]
[
K
]{
Δδ
}
-
[
L
]{
Δ
p
}
=
{
F
}
-
{
I
} (
3
)
和渗流连续性方程 [
B
^
]
T
{
v
}
+
[
H
^
]
T
{
p
}
=
{
Q
} (
4
)
式中:[
K
]是通常的刚度矩阵;{
Δδ
}为节点位移增量;[
L
]为节点孔隙水压力所对应的节点力;{
Δ
p
}为
孔隙水压力增量;{
F
}和{
I
}分别为节点外载荷和增量迭代过程中上一增量步中的不平衡力;[
B
^
]为节点
变形率对应的体积改变率;{
v
}为位移{
δ
}对时间的导数;[
H
^
]为孔压变化对应的体积改变率;{
p
}为孔隙
水压力,{
Q
}为节点流量
.
1.3
扩展的
MohrCoulomb
屈服准则及屈服接近度
Fredlund
等(
1978
年)提出非饱和土的抗剪强度公式可用独立的应力状态变量来表达,在抗剪强度
公式中可采用三个应力状态变量中的两个,现采用应力状态 变量(
σ
n
-
u
w
)和 (
u
a
-
u
w
)来 表示非 饱和土
的抗剪强度公式
[
6
]
τ
f
=
c′
+
(
σ
n
-
u
w
)
tan
φ
′
+
(
u
a
-
u
w
)
tan
φ
″
(
5
)
其中:
τ
f
为非饱和土抗剪强度;(
σ
n
-
u
w
)和(
u
a
-
u
w
)分别 为 净法向 应力和 基质吸 力;
c′
(
ε
p
l
,
θ
,
f
α
)为有 效
粘聚力,
φ
′
(
θ
,
f
α
)为有效内摩擦角,
ε
p
l
是等效塑性应变,
θ
是温度,
f
α
,
α
=1
,
2
,…是用户预先定义场变量;
φ
″
为与基质吸力有关的内摩擦角,且
-
φ
′
≤
φ
″
≤
0
;
σ
n
为法向总应力;由式(
4
)可见非饱和土的抗剪强度
除了与
c′
、
φ
′
以及
σ
n
有关,还与基质吸 力(
u
a
-
u
w
)有关
.
当土体 趋于饱 和时,基质 吸力趋 于 零,则 式(
5
)
变为
804
西
安
建
筑
科
技
大
学
学
报(自然科学版)
第
42
卷
剩余7页未读,继续阅读
资源评论
weixin_38575536
- 粉丝: 3
- 资源: 926
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功