收稿日期 :2008‐10‐18 ;修改稿收到日期 :2009‐12‐16畅
基金项目 :国家自然科学基金(10872050)资助项目 畅
作者简介 :王燕昌(1956‐) ,男 ,教授 ;
郑 静
倡
(1982‐) ,女 ,硕士
(E‐mail :zhengjing_abc@ 163 .com) .
第27卷第6期
2010 年 12 月
计 算 力 学 学 报
Chinese Journal of Computational Mechanics
Vol .27 ,No .6
December 2010
文章编号 :1007‐4708(2010)06‐1011‐05
径 向 积 分 边 界 元 法 确 定 非 均 质 土 石 坝 渗 流 自 由 面
王燕昌
1
, 郑 静
倡 1
, 高效伟
2
(1 .宁夏大学 固体力学研究所 ,银川 750021 ;2 .大连理工大学 工业装备结构分析国家重点实验室 ,大连 116024)
摘 要 :浸润面位置的确定是无压渗流分析中非常重要的问题 ,本文将径向积分边界单元法应用到渗流问题中 ,
通过直接将区域积分转化为边界积分的技术 ,克服了传统边界元解决该类问题的缺陷 ,编写了迭代计算程序 ,并
用二 、三维算例证明了算法的有效性 。
关键词 :径向积分法 ;边界单元法 ;浸润线
中图分类号 :TV223 .4 ;O241 文献标识码 :A
1 引 言
土石坝是挡水建筑物 ,它与渗流同时并存 。 在
土石坝设计中 ,渗流分析是有效降低土石坝渗透破
坏 、布设防渗措施的重要方面 。 所以关于渗流问题
的数值计算具有重要的应用价值 。 土石坝渗流为
无压渗流 ,有浸润面 ,可将其视为稳定层流 ,满足达
西定律 。 浸润面的确定是无压渗流分析中非常重
要的问题
[1]
。 浸润线是介质中的渗流水平面(又称
浸润面)与其铅垂面的交线 ,是渗流的自由液面 。
对于无压渗流而言 ,渗流分析中重要的问题是作为
渗流边界一部分的渗流自由面事先是未知的 ,因而
使这类问题成为非线性问题 ,解决此类问题的方法
是迭代法 。
柴军瑞和杨德全等
[2‐4]
曾用边界元法研究各向
同性浸润线的位置 ,在处理变系数问题时 ,需要将
计算域划分为块 ,每块边界上渗透系数假设相同 。
本文研究各向同性变渗透系数下自由面的确定 ,不
需要将计算域划分为块 。
2 边界单元法
[5‐8]
2 .1 边界单元法基本方程
渗流问题稳定场的控制方程如下 :
磹 k 磹 H
+
Q
=
0 (1)
式中 k 是渗透系数 ,是坐标的函数 ,H 是水头 ,Q
是区域
Ω
内的汇源项 ,是一个已知项 。
边界条件为
H
|
Γ
1
=
H
v
|
Γ
2
= -
(k
抄 H
抄 x
n
x
+
k
抄 H
抄
y
n
y
) = v
(2)
式中 H 和 v 是边界
Γ
1
和
Γ
2
上的已知函数 ,v 为渗
流速度 ,且
Γ
=
Γ
1
+
Γ
2
。
引用高斯公式并通过相应的转化 ,方程(1) 可
表示为
[9]
∫
Ω
G
抄
抄 x
i
(k
抄 H
抄 x
i
)d
Ω
=
∫
Γ
Gk
抄 H
抄 x
i
n
i
d
Γ
-
∫
Γ
k H
抄G
抄x
i
n
i
d
Γ
+
∫
Ω
H
抄G
抄 x
i
d
Ω
+
∫
Ω
k H
抄
抄x
i
(
抄G
抄x
i
)d
Ω
(3)
式中 G
=
1
2
π
log(
1
r
) (对于二维问题)
1
4
π
r
(对于三维问题)
式中 i
=
1 为 x 方向的值 ,i
=
2 为
y
方向的值 ,
Γ
是区域上的边界 。
通过转化 ,式(3) 可以转化为解决变渗透系数
渗流问题边界积分方程的标准方程 :
珮
H (
y
) = -
∫
Γ
G (x ,
y
)v(x)d
Γ
(x) -
∫
Γ
抄G(x ,
y
)
抄n
珮
H (x)d
Γ
(x) +
∫
Ω
G (x ,
y
)Q(x)d
Ω
(x) +
∫
Ω
V (x ,
y
)
珮
H (x)d
Ω
(x) (4)
式中
珮
H (x) = k(x)H(x) ,k
~
(x) = log k(x)
V x ,
y
=
抄 k
~
(x)
抄x
i
抄G(x ,
y
)
抄x
i