淤积泥砂对垂直地运动作历时刚性坝面动压力的影响研究.docx
- 文档编号:27516666
- 上传时间:2023-07-02
- 格式:DOCX
- 页数:11
- 大小:19.90KB
淤积泥砂对垂直地运动作历时刚性坝面动压力的影响研究.docx
《淤积泥砂对垂直地运动作历时刚性坝面动压力的影响研究.docx》由会员分享,可在线阅读,更多相关《淤积泥砂对垂直地运动作历时刚性坝面动压力的影响研究.docx(11页珍藏版)》请在冰豆网上搜索。
淤积泥砂对垂直地运动作历时刚性坝面动压力的影响研究
淤积泥砂对垂直地运动作历时刚性坝面动压力的影响研究
摘要:
作者应历时域显式有限元模型分析了可紧缩库水条件下淤积泥砂对垂直地运动作历时刚性坝面上动压力的影响。
文中首先通过与解析解的比较,验证了计算模式的正确性。
随后的分析中研究了作为液固两相介质处置的淤砂层的饱和度、厚度、渗透系数和孔隙率对坝面动压力的影响,和作为液固两相介质、单相固体、单相重流体介质处置的淤砂层模型间的计算结果不同。
计算结果表明:
饱和淤砂层使得库水自振频率增加,而非饱和淤砂层使得库水自振频率明显减小。
关键词:
淤积泥砂垂直地运动坝面动压力
自从Westergaard[1]究了简谐水平地运动作用下二维刚性直立坝面的动水压力反映以来,许多学者为合理肯定地震时坝面动水压力的研究作出了贡[2]献。
但绝大多数的研究功效中都没有考虑库底淤积泥砂层的影响。
Fenves和Chopra[3]将淤砂层和地基的作用模拟为一种库底吸收边界条件,研究了表示库底吸收特性的压力波反射系数α对混凝土重力坝频域反映函数的影响,结果表明库底吸收边界条件对坝体频域反映函数有明显的影响,但文中没有给出压力波反射系数α的具体取值方式。
Medina等人[4]将库底淤积砂层模拟为线性粘弹性介质,研究了淤砂层对坝体反映的影响。
傅作新和陆瑞明[5]假设库底为成层弹性介质,成立了一般库底的吸收边界条件。
通常库底淤积泥砂的孔隙率很大,很显然将淤砂层模拟为包括固体骨架和孔隙水的液固两相介质应该是更为合理。
Cheng[6]将淤砂层模拟为多孔弹性介质,求得了垂直地震动时刚性坝面动水压力散布的解析解(一维模型),数值结果表明:
对于饱和孔隙水的情况,库水和淤砂的彼此作用可以忽略;但对于孔隙水略微不饱和的情况,即即是很薄的淤砂层能够明显改变库水反映曲线,随着孔隙水可紧缩性增加,共振峰对应频率减小,且峰值增加。
Bougaha和Tassoulas[7,8]将淤砂层模拟为两相多孔弹性介质,用超单元法研究了淤砂层对混凝土重力坝地震反映的影响。
Chen和Hung[9]将坝体、地基假定为刚性,淤砂层模拟为流体饱和的两相多孔介质,利用隐式有限差分法结合快速泊松解法(fastPoissionsolver)分析了在恒定地面加速度的作用下,淤砂层对坝面动力反映的影响。
Dominguez等人[10]利用边界元法分析了淤砂层对混凝土重力坝地震反映的影响。
杜修力、王进廷和通过室内振动台实验和理论方式,分析了垂直地运动作历时,淤砂层对刚性坝面动压力的影响。
上述研究都表明库底淤积泥砂对坝面动压力有着明显影响,淤砂层的存在有可能增加坝面的动压力。
本文基于单相固体、单相流体、液固两相介质的时域显式有限元波动分析方式[12~15],研究了可紧缩库水条件下淤积泥砂对垂直地运动时作用于刚性坝面上的动压力的影响。
研究中分析了作为液固两相介质处置的淤积泥砂的饱和度、厚度、渗透系数和孔隙率的影响,还比较了作为液固两相介质、单相固体、单相重流体处置的淤积泥砂层模型间的计算结果不同。
文中频域结果为时域计算结果直接进行傅氏变换后的表述形式。
1分析模型
分析模型如图1所示,坝体和地基为刚性,上游坝面直立;库水和淤砂层总深度H=dw+ds=150m,其中dw为库水深度,ds为淤砂层厚度,dw、ds沿上游方向维持不变;上游人工边界距坝面距离L=500m,人工边界采用多次透射人工边界的修正公式[16]。
波动方程及边界条件
库水波动方程在小扰动情况下,粘性可紧缩流体波动方程为[17]:
(1)
式中:
、、别离表示流体位移、速度和加速度向量;ρw为流体密度;Kw为流体体积模量;ζ、η别离为流体第一,第二粘性系数。
取ζ=0、η=0,式
(1)即为库水(理想可紧缩流体)波动方程。
液固两相介质波动方程本文采用Biot[18]液固两相介质的波动方程,如下:
(2)
(3)
式中:
u为固体骨架位移向量;U为平均孔隙液体位移向量;e=Δ·u;ε=Δ·U;G、A与固相骨架Lame常数μ、λ有关;R为维持多孔介质的整体积不变时,压入必然体积的液体时所需要的力;Q表示液相体积转变与固相体积转变之间的耦合关系;ρ11=(1-n)ρs+ρa;ρ12=-ρa;ρ22=nρl+ρa;耗散常数b=n2ρlg/k;其中ρs、ρl、ρa别离为固相密度、液相密度、附加视质量(很难测定,一般取ρa=0),n为孔隙率,g为重力加速度,k为达西渗透系数;Biot弹性常数G、A、R、Q由实验测定[19]。
当假定骨架颗粒不可紧缩时,R=nKl,Q=(1-n)Kl,G=μ,A=λ+Q2/R,其中Kl为孔隙液体积模量,对于饱和介质Kl=Kw,对于非饱和介质孔隙液体体积模量由下式肯定
(4)
其中:
S表示饱和度;p0为绝对孔隙静压力。
边界条件
(1)库水-固体边界①库水法向位移与固体边界法向位移持续。
②切向力为零。
(2)库水-液固两相介质交壤面①两相介质孔隙水压强和库水压强平衡。
②两相介质总法向力与库水压力平衡。
③两相介质骨架切向应力为零。
④两相介质法向总位移与库水法向位移持续。
(3)液固两相介质-刚体交壤面①两相介质液相、固相法向位移与刚性边界法位移持续。
②两相介质液相、固相切向应力为零。
显式有限元计算格式液固两相介质、库水和其交壤面波动分析均采历时域显式有限元法。
由于篇幅有限,本文再也不介绍,可参见文献[12~15]。
模型参数
库水模型参数库水为理想流体,密度ρw=1000kg/m3,体积模量Kw=×109Pa,为了消除多次透射人工边界的高频振荡,库水中引入很小的不影响结果的瑞利(Rayleigh)阻尼,阻尼系数αw=0,βw=。
淤砂层模型参数
(1)液固两相介质模型。
本文取固体骨架密度ρs=2640kg/m3,孔隙水密度ρl=1000kg/m3,并假定固体骨架颗粒不可紧缩,其它参数取值见表1。
表1淤砂层(液固两相介质)材料参数(SI)
组别
n
k
b(×106)
S(%)
Kl(×109)
λ(×107)
G(μ)(×107)
A(×108)
R(×109)
Q(×109)
1
2
3
4
5
6
7
100
100
100
100
100
(2)单相固体介质模型。
假定液固两相介质渗透系数为零,则固相和液相位移、速度、加速度向量相等,把液固两相介质Biot波动方程式
(2)和式(3)相加,则可得单相固体介质动力方程。
由此可得液固两相介质和单相固体介质参数的转换关系为μ′=G,λ′=A+2Q+R,其中λ′、μ′为单相固体介质的拉梅常数。
表1中液固两相介质第1组和第6组参数对应的单相固体介质参数见表2。
(3)单相重流体模型。
液固两相介质参数与含有固体颗粒的重流体参数转换关系
(5)
(6)
其中:
ρ'l为重流体密度;K′l为重流体体积模量;Es为固体骨架颗粒体积模量,若假定固体骨架颗粒不可紧缩,则Es为无穷大。
表1中液固两相介质第1组和第6组参数对应的重流体参数见表2。
表2与液固两相介质相对应的重流体和固体介质参数(SI)
组别
液固两相介质
含有固体颗粒的重流体
单相固体介质
密度
体积模量Kl
密度
λ′
μ′
1
2
表1中第1组参数
表1中第6组参数
1656
1656
×109
×108
1656
1656
×109
×108
×106
×106
2计算模型的验证
对于本文计算模式的验证可以直接通过与Cheng[6]给出的解析解的比较进行,但Cheng给出的计算公式较复杂,作者基于位移势函数也取得了该问题的解析解,图2为将10m厚的淤砂层别离模拟为液固两相介质(表1中的第一组数据)和单相流体(表2中第一组数据)时显式有限元模式计算取得的坝面动压力与解析解的比较,其中竖向轴表示受到简谐地面加速度g·eiωt作历时作用于上游坝面的库水和淤砂动压力积分值的模g·|F(ω)|和
作用于坝面静水压力积分值ρwgH2的比值,水平轴为无量纲频率ω/ωl,其中库水第一自振频率ω1=πcw/2H,cw为库水波速。
从图中可以看出,显式有限元计算结果与解析解大体吻合,表明了本文计算模式的正确性。
(a)液固两相介质
(b)单相流体
图2显示有限元计算结果与解析解的比较3淤积泥砂层对坝面动压力影响
淤砂层饱和度对坝面动压力影响图3(a、b)表示淤砂层厚度ds别离为10m和100m,淤砂层不同饱和度(对应于表1中第1组、第6组和第7组参数)对坝面动压力的影响。
从图中可以看出,在垂直地面加速度作历时,淤砂层的饱和度对坝面动压力有着明显的影响。
随着淤砂层饱和度的减小,自振频率逐渐减小。
对于厚
淤砂层这种趋势加倍明显。
(a)ds=10m
(b)ds=100m
图3饱和度对坝面动压力影响
淤砂层厚度对坝面动压力影响图4(a、b)表示淤砂层别离取表1中第1组和第6组参数时,淤砂层厚度的转变对坝面动压力的影响。
从图中可以看出,当淤砂层饱和时(图4(a)),随着淤砂层厚度的增加,库水的第一、二自振频率向高频方向移动,而当淤砂层非饱和时(图4(b)),随着淤砂层的增加,库水的第一自振频率及高阶自振频率向低频方向移动。
这是由于淤砂层饱和时,淤砂层中PI波波速1425m/s稍大于库水波速1414m/s,随着淤砂层厚度增加,库水和淤砂层中的平均纵波波速增大,从而致使系统的自振频率增高;而非饱和淤砂层中的PI波波速527m/s(S=%)小于库水波速1414m/s,随着淤砂层的增加,库水和淤砂层中的平均纵波波速减小,系统的自振频率也相应减小。
(a)饱和度S=100%
(b)饱和度S=%
图4淤沙层厚度对坝面动压力影响
淤砂层渗透系数对坝面动压力影响(饱和淤砂层)图5(a、b)表示淤砂层厚度ds别离为10m和100m时,不同淤砂层渗透系数(对应于表1中第2组、第1组和第3组参数)对坝面动压力的影响。
从图中可以看出,对于垂直地面加速度作用,渗透系数对自振峰值频率影响可以忽略。
(a)ds=10m
(b)ds=100m
图5渗透系数对坝面动压力影响
淤砂层孔隙率对坝面动压力影响(饱和淤砂层)图6(a、b)表示淤砂层厚度ds别离为10m和100m时,不同孔隙率(对应于表1中第1组、第4组和第5组参数)对坝面动压力的影响。
从图中可以看出随着孔隙率的减小库水和淤砂系统的自振频率逐渐增加。
其原因在于随着孔隙率的减小PI波波速逐渐增加,别离为1425m/s(n=、1490m/s(n=和1594m/s(n=,因此对于一样深度的淤砂层,孔隙率小时PI波波速快,库水和淤砂系统自振频率大;而孔隙率大时,PI波波速小,库水和淤砂系统自振频率小。
(a)ds=10m
(b)ds=100m
图6孔隙率对坝面动压力影响
4将淤砂层模拟为液固两相介质、重流体介质、单相固体介质时的比较
图7(a、b、c、d)为将淤砂层别离模拟为液固两相介质、重流体介质、单相固体介质时,坝面动压力的比较图。
从图中可以看出,将淤砂层模拟为单相固体介质与将淤砂层模拟为液固两相饱和介质完全吻合。
这是因为通过上述液固两相介质与单相固体介质参数间的换算关系,使得单相固体介质的P波和S波波速与液固两相介质的PI波和S波波速大体一致,而液固两相介质的PⅡ波衰减很快,对坝面动压力影响很小,因此二者完全吻合。
将淤砂层模拟为重流体介质时,只能传播P波,不能传播S波,因此将淤砂层别离模拟为重流体和液固两相介质时,坝面动压力共振频率稍有不同,重流体介质趋向于减小自振频率。
这说明本地面运动致使淤砂层液化后,可能使得坝面动压力峰值对应的频率略有减小。
5结论
本文研究了在刚性坝、刚性地基条件下,垂直地运动作历时,淤砂层对坝面动压力的影响,比较了作为液固两相介质、单相固体、单相重流体处置的淤积泥砂层模型间的计算结果。
通过研究,针对本文的计算参数,可以得出如下一些结论:
(1)淤砂层饱和度对坝面动压力共振频率有重要的影响,对于必然厚度的淤砂层,随着淤砂层饱和度的减小,共振频率慢慢减小。
(2)淤砂层厚度可以改变坝面动压力共振频率,对于饱和淤砂层,随着淤砂层厚度的增加,共振频率增加;而对于非饱和淤砂层,随着淤砂层厚度的增加,共振频率减小。
(3)淤砂层渗透系数对坝面动压力的共振频率影响可以忽略。
(4)随着淤砂层孔隙率的减小,坝面动压力共振频率增加。
(5)将淤砂层作为单相固体介质处置与作为液固两相介质处置的不同可以忽略,但将淤砂作为重流体处置时,坝面动水压力共振频率略有减小,也就是说,当淤砂层液化时,坝面动压力共振频率趋向于减小。
图7不同淤沙层模型计算结果比较
参考文献:
[1]WestergaardHM.Waterpressuresondamsduringearthquakes[J].TransactionsofASCE,1933,98:
418-433.
[2]杜修力,王进廷.动水压力及其对坝体地震反映影响的研究进展[J].水利学报,2001,(7):
13-21.
[3]FenvesG,ChopraAK.Effectsofreservoirbottomabsorptionanddam-water-foundationrockinteractiononfrequencyresponsefunctionsforconcretegravitydams[J].EarthquakeEngineeringandStructuralDynamics,1985,13
(1):
13-31.
[4]MedinaF,DominguezJ,TassoulasJL.Responseofdamstoearthq
uakeinculdingeffectsofsediment[J].JournalofStructuralEngineering,ASCE,1990,116
(1):
3108-3121.
[5]傅作新,陆瑞明.水库的库底条件和挡水结构的动水压力[J].水利学报,1987,(5):
28-34.
[6]ChengAHD.Effectsofsedimentonearthquakeinducedreservoirhydrodynamicresponse[J].JournalofEngineeringMechanics,ASCE,1986,112(7):
654-663.
[7]BougachaS,TassoulasJL.SeismicresponseofgravitydamsⅠ:
Modelingofsediments[J].JournalofEngineeringMechanics,ASCE,1991,117(8):
1826-1837.
[8]BougachaS,TassoulasJL.SeismicresponseofgravitydamsⅡ:
Effectsofsediments[J].JournalofEngineeringMechanics,ASCE,1991,117(8):
1839-1950.
[9]ChenBF,HungTK.Hydrodynamicpressureofwaterandsedimentonrigiddam[J].JournalofEngineeringMechanics,ASCE,1993,119(7):
1411-1433.
[10]DominguezJ,GallegoR,JaponBR.Effectsofporoussedimentsonseismicresponseofconcretegravitydams[J].JournalofEngineeringMechanics,ASCE,1997,123(4):
302-311.
[11]DuXiuli,WangJinting,HungTinkan.EffectsofSedimentontheDynamicPressureofWaterandSedimentondams[J].ChineseScienceBulletin,2001,46(7):
521-524.
[12]杜修力,王进廷.阻尼弹性结构动力计算的显式差分法[J].工程力学,2000,17(5):
37-43.
[13]王进廷,杜修力,赵成刚.液固两相饱和介质波动分析的一种显式有限元法[J].岩石力学与工程学报,2002,21(8):
1199-1204.
[14]王进廷,杜修力,赵成刚.分层分介质中波动传播分析的显式有限元法[J].岩石力学与工程学报,2003,22
(1).
[15]王进廷.高混凝土坝-可紧缩库水-淤砂-地基系统地震反映分析研究[D].北京:
中国水利水电科学研究院,2001.
[16]李小军.非线性场地地震反映分析方式的研究[D].哈尔滨:
国家地震局工程力学研究所,1993.
[17]王维佳.粘性可紧缩流体中的波[J].地球物理学报,1984,27
(1):
84-93.
[18]BiotMA.Theoryofpropagationofelasticwavesinafluidsaturatedporoussolid[J].JournaloftheAcousticalSocietyofAmerica,1965,28,168-191.
[19]BiotMA,WillisDG.Theelasticcoefficientsofthetheoryofconsolidation[J].JournalofAppliedMechanics,1957,24:
594-601.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 淤积 垂直 运动 历时 刚性 坝面动 压力 影响 研究