简谐荷载作用下非饱和土地基中波阻板隔振效果研究.pdf
《简谐荷载作用下非饱和土地基中波阻板隔振效果研究.pdf》由会员分享,可在线阅读,更多相关《简谐荷载作用下非饱和土地基中波阻板隔振效果研究.pdf(10页珍藏版)》请在咨信网上搜索。
1、第 36 卷第 4 期2023 年 8 月振 动 工 程 学 报Journal of Vibration EngineeringVol.36 No.4Aug.2023简谐荷载作用下非饱和土地基中波阻板隔振效果研究张猛1,马强1,2(1.青海大学土木工程学院,青海 西宁 810016;2.青海省建筑节能材料与工程安全重点实验室,青海 西宁 810016)摘要:鉴于非饱和土地基振动问题的普遍性和复杂性,环境振动下非饱和土地基振动控制已成为土动力学的研究热点。基于单相弹性介质和非饱和多孔介质理论,对简谐荷载作用下非饱和土地基中设置单相固体波阻板的隔振效果进行了研究。考虑地表排水排气的边界条件,利用
2、Fourier积分变换和 Helmholtz矢量分解原理,建立了动荷载作用下地基动力响应的计算列式。分析了非饱和土地基中土体饱和度、荷载频率、波阻板的埋深、厚度以及弹性模量对其隔振性能的影响规律。结果表明:非饱和土地基中设置波阻板能够取得很好的隔振效果。地表位移幅值随饱和度和波阻板埋深的减小而显著降低,随荷载频率、波阻板的厚度和弹性模量的增大而明显减小。关键词:非饱和土地基;波阻板;隔振效果;动力响应中图分类号:TU435 文献标志码:A 文章编号:1004-4523(2023)04-1136-10 DOI:10.16385/ki.issn.1004-4523.2023.04.027引 言随着
3、城镇化建设和现代化工业的迅速发展,各种人工振动引起的振动污染愈加频繁,如轨道交通、爆破等工程活动工作时产生的振动对各种精密仪器的正常运作以及人们的生活环境和工作环境产生了诸多不利影响。因此,对人工振动引起的振动传播过程和地基的振动规律进行分析研究,从而找到有效的减振隔振措施具有重要的实际意义,隔振理论和方法成为土动力学研究的热点。在既有振源与保护区之间设立屏障来切断弹性波的传播路径,消耗振动能量,降低振动幅度,是目前国际上普遍采用的隔振措施。Woods1最早通过现场原位试验对远场被动隔振和近场主动隔振进行了研究,提出了屏障隔振设计的基本准则,给出了衡量屏障隔振效果的系数振幅衰减系数。之后国内外
4、的学者对各种隔振屏障的减隔振性能开展了大量的研究工作212。根据隔振屏障的形式,高广运12按几何特性进一步将屏障分为连续屏障(如空沟、填充沟等)和非连续屏障(如桩列和板桩等)。除了以上所述的屏障隔振外,另一种可供选择的隔振屏障是在振源或被保护结构下一定深度内埋置波阻板进行隔振。Chouw 等13最先分析了弹性地基中波阻板主动和被动隔振效果,并对填充沟和波阻板的被动隔振效果进行了对比分析,结果表明波阻板的隔振效果要优于填充沟。Takemiya 等14采用波阻板研究了基岩上单一土层中群桩基础激振时的隔振效果,结果表明波阻板是一种有效的隔振措施。周凤玺等15研究了弹性地基中含液饱和多孔波阻板的隔振性
5、能。马强等16建立了弹性地基中功能梯度波阻板的地基隔振体系,研究表明梯度波阻板能有效地降低振动的振幅。田抒平等17和高盟等18提出一种带孔波阻板填充 Duxseal的联合隔振方法(简称 DXWIB),通过试验研究发现较浅埋深的 DXWIB 屏障可以有效地减小地表竖向位移。高广运等1920针对饱和土地基模型,对轨道交通载荷作用下饱和土地基中波阻板的隔振性能进行了研究。随后基于改进的三维边界有限元模型,Gao 等21研究了饱和土体中波阻板的隔振效果,分析了土地基波阻板相互作用的问题。徐长节等22建立了饱和土地基中夹水混凝土复合屏障的计算模型,结果表明,入射波的入射角度及屏障的弹性模量对屏障的隔振效
6、果影响最为明显。综上可知,目前关于波阻板隔振屏障的研究大都限于单相弹性地基与饱和土地基中。然而,非饱和土地基才是自然界中土体更为普遍的存在状态。地基振动由于振源位置、振源类型以及地基物理力学性质的不同而产生不同的波场特性,其振动的传播过程和衰减规律也不相同。因此,研究非饱和土收稿日期:2022-01-20;修订日期:2022-03-12基金项目:国家自然科学基金资助项目(52168053,51978320);青海省自然科学基金青年基金资助项目(2021-ZJ-943Q)。第 4 期张 猛,等:简谐荷载作用下非饱和土地基中波阻板隔振效果研究地基中波阻板的减振隔振效果具有更普遍意义。本文基于非饱和
7、多孔介质的控制方程,建立了非饱和土地基中设置均质波阻板的数学模型,利用 Fourier积分变换,通过 Helmholtz矢量分解原理,推导获得了土体在 Fourier变换域中动力问题的位移、应力的通解。通过数值算例,研究了非饱和土地基中设置波阻板后的隔振效果,且与饱和土地基比较了隔振效果的区别,分析了非饱和土地基中饱和度、荷载频率、波阻板埋深、厚度以及波阻板弹性模量对地基隔振性能的影响规律。1控制方程1.1单向固体介质各向同性的线弹性单相连续固体介质的控制方程:(e+e)(ue)+e2ue=eu e(1)式中 e和e表示固体材料的 Lame弹性常数;e为固体材料密度;ue表示位移矢量。相应的应
8、力位移关系为:e=e(ue)I+e(ue+ue)I(2)式中 e表示单相弹性固体中的应力张量;I为单位矩阵。根据 Helmholtz 矢量分解原理,位移矢量ue可以用势函数表示为:ue=e+e(3)式中 e和e分别为单相固体的标量位势函数和矢量位势函数。将式(3)代入式(1)中,整理得到弹性固体波动方程:2e=1v2p2et2,2e=1v2s2et2(4)式中 vp=()e+2e/e,vs=e/e分别为纵波和横波的传播速度。1.2非饱和多孔介质考虑土体受到简谐荷载作用,在稳态情况下,所有变量都可写成:f=f*eit(5)为了方便,后面的推导过程略去星号。基于连续介质力学和 Bishop 有效应
9、力公式,徐明江等23提出的非饱和土的动力控制方程为:p2us+()p+p us-apl-a()1-pg=-2 sus-2 lul-2 gug (6a)-pl=bli(ul-us)-2lul(6b)-pg=bgi(ug-us)-2gug(6c)-pl=a11 us+a12 ul+a13 ug(6d)-pg=a21 us+a22 ul+a23 ug(6e)式中 a11 a23,A11 A24,bl,bg,a,等系数详见文献23;us,ul,ug,s=(1-n)s,l=nSrg,g=n(1-Sr)g分别为非饱和土固、液、气相的位移和相对密度;s,l,g分别为固、液、气相的密度;为角频率;n为孔隙率;
10、Sr为饱和度;p,p为非饱和土的 Lame 系 数;pl,pg分 别 为 孔 隙 水 压 力、孔 隙 气压力。非饱和土的应力应变关系为:ij=peij+2pij-ijap(7)式中 ij为非饱和土体介质的总应力分量(i,j=1,3),e=us为土骨架体应变,ij为克罗内克符号,ij为土骨架的应变,p=pl+(1-)pg为等效孔隙流体压力。将式(6d),(6e)代入式(6a),(6b)和(6c)整理后可得:p2us+B1()us+B2()ul+B3()ug=-2 sus-2 lul-2 gug (8a)a11()us+a12()ul+a13()ug=C1us+C2ul(8b)a21()us+a2
11、2()ul+a23()ug=C3us+C4ug(8c)式中 B1=p+p+aa11+a(1-)a21,B2=aa12+a(1-)a22,B3=aa13+a(1-)a23,C1=-bli,C2=bli-2l,C3=-bgi,C4=bgi-2g。根据 Helmholtz矢量分解原理,位移矢量us,ul,ug可以用势函数表示为:uz=z+z(9)式中 z=s,l,g;s,l,g分别为固、液、气相的标量位势函数;s,l,g分别为固、液、气相的矢量位势函数。将式(9)代入式(8)中进行散度和旋度运算,整理可得:()p+B12s+B22l+B32g=-2 ss-2 ll-2 gg(10a)a112s+a1
12、22l+a132g=C1s+C2l (10b)a212s+a222l+a232g=C3s+C4g (10c)p2s=-2 ss-2 ll-2 gg (10d)C1s+C2l=0(10e)1137振 动 工 程 学 报第 36 卷C3s+C4g=0(10f)2位移势函数的解考虑非饱和土地基的表面受到频率为,幅值为q0的竖向条形简谐荷载q=q0eit作用,如图 1所示。在距离非饱和土地基表面深度为 H处设置厚度为 hw的波阻板。波阻板将非饱和土地基分成和两部分,图1中给出了波的传播以及波幅的示意图。2.1单向固体介质的解对空间变量x1进行 Fourier变换:f(,x3)=-f(x1,x3)e-i
13、x1dx1(11)将式(5)和式(11)代入方程(4)中,经过 Fourier变换整理可得:d2edx23+2e e=0,d2edx23+2ee=0(12)式中 e=2/v2p-2,e=2/v2s-2,故可得到单相弹性固体介质中位移势函数的通解为:e=Atee-iex3+Areeiex3(13a)e=Btee-iex3+Breeiex3(13b)式中 Ate和Are分别为单相固体介质中透射 P 波和反射 P波的波幅,Bte和Bre分别为单相固体介质中透射S波和反射S波的波幅。2.2非饱和多孔介质的解将式(11)代入方程(10a),(10b)和(10c)中,经过 Fourier变换整理可得:(p
14、+B1)d2sdx23+b11 s+B2d2ldx23+b12 l+B3d2gdx23+b13 g=0(14a)a11d2sdx23+b21 s+a12d2ldx23+b22 l+a13d2gdx23+b23 g=0(14b)a21d2sdx23+b31 s+a22d2ldx23+b32 l+a23d2gdx23+b33 g=0(14c)式 中 b11=s2-(p+B1)2,b12=l2-B22,b13=g2-B32,b21=-C1-a112,b22=-C2-a122,b23=-a132,b31=-C3-a212,b32=-a222,b33=-C4-a232。设方程组(14)的解为:s l g
15、T=csclcgTexp(x3)(15)将式(15)代入式(14)得到线性方程组:2B3+b132B2+b122(p+B1)+b112a13+b232a12+b222a11+b212a23+b332a22+b322a21+b31cgclcs=0(16)式(16)有 非 零 解 的 条 件 为 系 数 矩 阵 行 列 式为 0,即:16+24+32+4=0(17)设式(17)的根为n(n=1,2,3),则n由下式给出:n=rn(Ren 0,n=1,2,3)(18)式中 n(n=1,2,3,4)见附录 A。可得到常微分方程组(14)的解为:s=n=13()Atpne-nx3+Arpnenx3(19
16、a)l=n=13lpn()Atpne-nx3+Arpnenx3(19b)g=n=13gpn()Atpne-nx3+Arpnenx3(19c)式中 lpn,gpn见附录 A;Atpn和Arpn(n=1,2,3)分别为非饱和土多孔介质中透射 P1波,P2波,P3波和反射P1波,P2波,P3波的波幅。将式(11)代入方程(10d),(10e)和(10f)中,经过 Fourier变换整理得:pdsdx23+d11s+d12l+d13g=0(20a)d21s+d22l=0(20b)d31s+d32l=0(20c)式 中 d11=s2-p2,d12=l2,d13=g2,d21=C1,d22=C2,d31=
17、C3,d32=C4。设方程组(20)的解为:slgT=dsdldgTexp(rx3)(21)将式(21)代入式(20)得到线性方程组:图 1 简谐荷载作用下非饱和土地基设置波阻板的波传播示意图Fig.1 Wave propagation diagram of wave impeding block in unsaturated soil under harmonic load1138第 4 期张 猛,等:简谐荷载作用下非饱和土地基中波阻板隔振效果研究 d13d12pr2+d110d22d21d330d31dgdlds=0(22)式(22)有 非 零 解 的 条 件 为 系 数 矩 阵 行 列
18、式为 0,即:5r2+6=0(23)式(23)的根为r,则r由下式给出:r=-6/5 (Rer 0)(24)式中 n(n=5,6)见附录 A。可得到常微分方程组(20)的解为:s=Btse-rx3+Brserx3(25a)l=ls(Btse-rx3+Brserx3)(25b)g=gs(Btse-rx3+Brserx3)(25c)式中 ls,gs见附录 A;Bts和Brs分别为非饱和土多孔介质中透射S波和反射S波的波幅。3非饱和土地基动力响应求解3.1地基动力响应在直角坐标系中,各位移分量可用位移势函数和表示为:u1=x1-x3,u3=x3+x1(26)将式(13)和(26)代入式(2)和(3)
19、中,将式(19),(25)和(26)代入式(6d),(6e),(7)和(9)中,可分别得到弹性介质和非饱和多孔介质在 Fourier 变换域中的位移和应力表达式。在0 x3 H的区域内:usI1=in=13(AItpne-nx3+AIrpnenx3)+r(BItse-rx3-BIrserx3)(27a)usI3=n=13-n(AItpne-nx3-AIrpnenx3)+i(BItse-rx3+BIrserx3)(27b)ulI3=n=13-nlpn()AItpne-nx3-AIrpnenx3+ils(BItse-rx3+BIrserx3)(27c)ugI3=n=13-ngpn(AItpne-n
20、x3-AIrpnenx3)+igs(BItse-rx3+BIrserx3)(27d)I13=2pin=13-n(AItpne-nx3-AIrpnenx3)-p(r2+2)(BItse-rx3+BIrserx3)(27e)I33=n=13n(AItpne-nx3+AIrpnenx3)-2irp(BItse-rx3-BIrserx3)(27f)plI=n=13(a11+a12lpn+a13gpn)(2-2n)()AItpne-nx3+AIrpnenx3(27g)pgI=n=13(a21+a22lpn+a23gpn)(2-2n)()AItpne-nx3+AIrpnenx3(27h)式中 n=a(a1
21、1+a12lpn+a13gpn)+a()1-()a21+a22lpn+a23gpn()2n-2+()p+2p2n-p2。在H x3 H+hw的区域内:ue1=i(Atee-iex3+Areeiex3)+ie(Btee-iex3-Breeiex3)(28a)ue3=-ie(Atee-iex3-Areeiex3)+i(Btee-iex3+Breeiex3)(28b)e13=2ee(Atee-iex3-Areeiex3)+e(2e-2)(Btee-iex3+Breeiex3)(28c)e33=-e2-2e(e+2e)(Atee-iex3+Areeiex3)+2ee(Btee-iex3-Breeiex
22、3)(28d)如果波阻板为非饱和土多孔介质,则此区域内的式(28)将退化为式(27)。在x3 H+hw的区域内:usII1=in=13(AIItpne-nx3)+rBIItse-rx3(29a)usII3=n=13(-nAIItpne-nx3)+iBIItse-rx3(29b)ulII3=n=13(-nlpnAIItpne-nx3)+ilsBIItse-rx3 (29c)ugII3=n=13(-ngpnAIItpne-nx3)+igsBIItse-rx3 (29d)II13=2pin=13(-nAIItpne-nx3)-p(r2+2)BIItse-rx3(29e)II33=n=13(nAIIt
23、pne-nx3)-2irpBIItse-rx3(29f)3.2边界条件及求解对于荷载作用在半平面表面,考虑地表透气透水的边界条件以及各层面处的连续条件:在x3=0处:I33=q0sin()ll,I13=0,plI=0,pgI=0(30)在x3=H处:I33=e33,I13=e13,usI3=ue3,usI1=ue1,ulI3=usI3,ugI3=usI3(31)在x3=H+hw处:e33=II33,e13=II13,ue3=usII3,ue1=usII1,1139振 动 工 程 学 报第 36 卷usII3=ulII3,usII3=ugII3(32)将式(27),(28)和(29)代入到边界条
24、件式(30),(31)和(32)中,可以得到如下的矩阵方程组:Tx=f(33)式中 矩阵 T 以及矢量 x 和f中的各元素详见附录 B。通过求解方程组(33),获得 x内各类波的波幅,结合式(27),(28)和(29)即可获得 Fourier变换域内非饱和土地基和波阻板中任意点的应力和位移响应。4数值算例对于非饱和土地基,由于饱和度的变化会引起土中剪切模量的改变,因此在本文中采用修正后的动剪切模量2324公式:p=s+2050lnS-2e-1+S-1etan (34)为了研究波阻板对非饱和土地基振动的控制效果,本文选取一组非饱和土地基的物理力学参数23如表 1 所示。选取均质波阻板的物理力学参
25、数如下:弹性模量Ee=6.5 108 Pa,泊松比v=0.3,密度e=2458 kg/m3,取荷载幅值q0=1 kPa,分布长度 l=1 m。由于被积函数表达式较为复杂,因此很难得到 Fourier逆变换的封闭形式解,本文采用 FFT方法完成 Fourier 逆变换,波数的离散点为 1024,空间计算区间为 100 m。为了分析非饱和土地基与饱和土地基中波阻板隔振效果的区别,考虑上覆土层厚度H=2 m,波阻板厚度 hw=2 m,波阻板弹性模量Ee=6.5 108 Pa,荷载频率=1 rad/s的情形下,本文将非饱和土地基退化到饱和土地基,非饱和土地基中考虑饱和度Sr=0.8,图 2,3 分别给
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 荷载 作用 饱和 土地 中波 阻板隔振 效果 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。