为了充分地利用土地,建筑与建筑之间的间距往往很小。日本的Kobe、法国的Nice、美国的Los Angeles、中国的上海等城市部分或全部修建在软土场地上,且位于地震带上或地震带的边缘。在地震作用下,土与结构之间的相互作用(soil-structure interaction, 简称SSI)是影响结构动力反应的关键因素之一[1-3]。局部区域内的建筑必将通过地基土对相邻结构的动力反应产生影响[4-5],由此存在结构-土-结构相互作用(structure-soil-structure interaction, 简称SSSI)问题。
Luco等[6]首先研究了SH波作用下SSSI问题。Wolf[7]对反应堆建筑、反应堆辅助室和燃料处理室组合结构进行地震反应分析。田彼得等[8]采用二维边界元法研究弹性半空上两个条形结构的动力阻抗和动力反应。姜忻良等[9]研究两结构的动力反应随相互之间距离的变化规律。窦立军[10]分析了高层建筑对相邻多层建筑的影响。Padron等[11]研究了桩基结构的SSSI。Clouteau等[12]的计算结果表明,与SSI体系相比SSSI体系可使结构顶层的反应降低。王淮峰等[13]研究了由于相邻结构存在情况下基地剪力的变化情况。
在一定区域范围内的结构通过土体耦合而成为一个有机的整体。结构之间距离和结构基础宽度等的变化相当于改变了系统的质量和刚度分布,由此导致结构动力特性发生相应的变化。而结构的动力特性是影响结构的动力反应最关键的因素之一。为此,本文将上部结构简化为等效单自由度体系,比较两相同结构组成的结构-土-结构相互作用体系与单一结构下土-结构相互作用体系动力特性的变化情况,系统地研究结构-土-结构体系动力特性随结构间的相对距离、结构与土体相对刚度比和基础宽度的变化规律。
上部结构简化为等效单自由度体系,基础采用明置于等厚度均质水平土层上的方形刚性基础,所构成的结构-土-结构相互作用体系如图 1所示。
图中B、L、H、D分别表示基础宽度、两基础间净距、土层厚度以及基础边界距人工边界的距离。将土体和结构进行有限元划分后,体系的无阻尼自由振动方程可表示为:
式中,[M]和[K]分别表示质量矩阵和刚度矩阵,{${\ddot u}$}和{u}分别为节点的加速度和位移向量。本文采用Ansys有限元软件的Block Lanczos法进行模态分析。建模时,上部结构的等效单自由度体系采用带集中质量的梁进行模拟。梁采用beam188单元,集中质量为mass21单元。刚性基础采用shell63单元,土体采用solid45单元,基础与土体无相对滑移。
土层实际为一半无限体,当采取用有限元法进行计算时,需要截取一定的范围进行计算,由此带来人工边界的问题[14-15]。下面基于SSI体系的有限元模型讨论侧向人工边界范围的选取问题。计算中侧向人工边界为自由边界,底部采用固定边界。SSI体系基频随D/B变化的计算结果如表 1所示。从计算结果可以看出,当D/B大于8后,结构基频的计算数值趋于稳定。因此,下面计算中侧向人工边界范围取为D/B=8,以消除侧向人工边界对结构基频的影响。计算中土体的参数如下:弹性模量为2.0×107 Pa,泊松比为0.3,密度为1 600 kg/m3,厚度为8 m。刚性基础宽度为8 m,忽略基础的质量。上部结构在刚性基础下的频率为13.63 rad/s。此时,结构和土层的频率比$\bar \omega $=1,基础宽度与土层厚度比$\bar B$=1。
从上部结构的计算模型看,刚性基础上的质量弹簧系统可以看作是采用筏基的建筑结构。设筏基的质量为mf,上部结构的质量为m。则不同mf/m下结构的基频如表 2所示。表中的误差为采用无质量刚性基础近似计算所得频率的相对误差。显然随着基础质量的增加,采用无质量刚性基础SSI体系计算频率的误差将增加。实际工程中,筏基结构的mf/m一般为5%,此时,采用无质量基础的计算误差小于0.5%。因此,为简化计算参数,在后面的计算分析中将基础简化为无质量的刚性基础。
图 2为不同距离下结构-土-结构相互作用体系自振频率和相应土-结构相互作用体系自振频率的计算结果。从计算结果可以看出,考虑土-结构相互作用后,结构体系的频率小于刚性基础下结构的频率,这一点与其它相关文献的理论和试验结果相同[16]。但是,对于单自由度下单一结构的土-结构相互作用体系而言,结构任一方向水平振动的主模态只有一阶,而与之对应的结构-土-结构相互作用体系的模态将是两阶频率相近的模态。其中一阶模态为两个结构同向运动,另一阶模态为两个结构相向运动,如图 3所示。
结构-土-结构体系中两阶频率相近模态是两个相邻结构通过土体耦合而产生的特殊的动力特征,具有对称性相反、频率相近的特点。为描述方便,将这样的两阶频率称为“孪生频率”。并定义孪生频率相对变化量为:
式中ωt2和ωt1分别为孪生频率的上限和下限,ωsi为与孪生频率对应的SSI体系的频率。
从图 2可以看出:1)结构-土-结构体系的孪生频率分别位于相应土-结构相互作用体系频率的两侧;2)当两个结构相距越远,ωt2和ωt1逐渐收敛到ωsi,这表明随着两个结构彼此之间距离的增加,两结构通过土体的耦合效应逐渐减少,当两者之间的距离达到一定程度后,这种相互作用可以忽略不计。而$\Delta {{\bar \omega }_{\rm{t}}}$综合反映了结构-土-结构相互作用体系相对于土-结构体系频率的变化情况,可作为描述结构之间相互影响的特征参数。
下面主要分析两相同结构所组成的结构-土-结构相互作用系统的频率变化情况。对结构体系频率的影响因素采用以下无量纲参数进行描述:结构和土层的频率比:${{\bar \omega }_{\rm{t}}} = {\omega _{{\rm{st}}}}/{\omega _{{\rm{so}}}}$;两结构相对距离:$\bar L = L/B$;宽厚比:$\bar B = B/H$。其中,ωst表示刚性基础下结构的频率,ωso表示仅均质土层的基频。宽厚比是基础宽度与土层厚度比,间接反映上部结构的细长比。
结构-土-结构相互作用对结构动力特性的影响是在土-结构相互作用体系动力特性的基础上,由于邻近构筑物存在而发生的变化。为更好的理解孪生频率的变化规律,下面首先分析土-结构相互作用体系频率ωsi随频率比${\bar \omega }$和宽厚比${\bar B}$的变化情况。计算中土层参数同前,保持不变。图 4为改变结构刚度和结构基础宽度情况下土-结构相互作用体系频率ωsi的变化规律。
正如所料,增加结构的频率比${\bar \omega }$将减少ωsi/ωst。这表明结构越刚,土层越软,土-结构相互作用影响越强。${\bar \omega }$→0,ωsi/ωst→1,土体相当于刚性基础,此时,不存在土-结构相互作用。从结构频率的绝对量的大小看,${\bar \omega }$ > 4以后,ωsi的值基本保持不变,这表明,当频率比大于4后,在其它条件不变的情况下,等效单自由度的频率取决于土层的频率,此时,增加上部结构的刚度对于土-结构相互作用体系的动力特性没有影响。
宽厚比反映了结构基础的大小。对于特别细长的结构(如:${\bar B}$=0.125),等效单自由度结构的频率将远小于结构自身的频率。此时,存在很明显土-结构相互作用是由基础尺寸太小而使土对结构的约束动刚度很小引起的。在其它条件不变条件下,${\bar B}$越大,结构的等效单自由度的频率越高,当${\bar B}$ > 2.5以后,不同基础宽度下的ωsi的值基本相同。这表明,当宽厚比大于2.5后,单一均匀土层的动刚度对结构动力特性的影响基本不变。
因此,在下面分析结构-土-结构相互作用体系孪生频率相对变化量$\Delta {{\bar \omega }_{\rm{t}}}$时,频率比${\bar \omega }$和宽厚比${\bar B}$的计算范围分别取为[0.2, 5]和[0.125, 4]。
图 5为不同宽厚比下孪生频率相对变化量$\Delta {{\bar \omega }_{\rm{t}}}$随结构间相对距离${\bar L}$的变化情况。对于任意宽厚比${\bar B}$,随着相对距离${\bar L}$增加,孪生频率相对变化量$\Delta {{\bar \omega }_{\rm{t}}}$均逐渐减小并趋于零。这是因为相邻结构相互作用是通过地基土来实现的,结构的运动迫使周边的土层发生相应的振动,由此导致邻近结构的振动。两结构间的距离越近,相互作用影响越大。当建筑之间的距离足够大时,相互作用影响很小,可以忽略不计。虽然很大范围内,相邻结构之间的动力特性都存在相互影响,但是模态分析主要作用是为时程分析做准备。在时程分析时,当距离大于一定值后,由于结构和土体阻尼的存在,而使邻近结构导致结构动力特性的改变以及结构有效输入改变对动力反应的影响可以忽略不计。综合其他文献的研究成果[17],可以发现,当${\bar L}$ > 2.5后,孪生频率相对变化量$\Delta {{\bar \omega }_{\rm{t}}}$小于5%,此时,可以认为结构基本不受相邻结构动力特性的影响。
从频率比的角度看,总的趋势是$\Delta {{\bar \omega }_{\rm{t}}}$随频率比的增加而增大,但并非单调变化。为更清楚的反映频率比的影响,图 6显示${\bar L}$=0.5情况下,$\Delta {{\bar \omega }_{\rm{t}}}$随频率比的变化情况。由图可知:1)当${\bar \omega }$ > 0.7时,$\Delta {{\bar \omega }_{\rm{t}}}$ < 5%。这是因为结构的刚度小于土层刚度时,相当于结构修建于较硬的场地上,此时,土-结构相互作用本身就很小,因此,通过土体对相邻构筑物的影响就更小,可以忽略不计。2)当${\bar \omega }$ > 4时,由于土-结构相互作用体系的频率取决于土层的频率,相应地$\Delta {{\bar \omega }_{\rm{t}}}$的也趋向于一个稳定的值。因此,对于${\bar \omega }$ > 4的结构体系,结构间的相互影响可以用${\bar \omega }$=4的计算结果近似表达。3)相邻结构对结构动力特性影响最大的区域通常在${\bar \omega }$∈[0.7, 2.5],$\Delta {{\bar \omega }_{\rm{t}}}$值甚至超过15%。但是其峰值通常在${\bar \omega }$ > 1时出现,这表明当结构刚度略大于土层刚度时,结构-土-结构的相互作用影响最大。
从宽厚比的角度看,当${\bar B}$ < 0.5时,$\Delta {{\bar \omega }_{\rm{t}}}$ < 5%。结合土-结构相互作用下频率的变化规律(图 4)可知,对于宽厚比小的结构,由于结构本身存在强烈的土-结构相互作用,相邻结构对结构动力特性的影响与之相比是一个小量。而当${\bar B}$ > 2.5时,基础的面积很大,相邻结构通过地基对结构的影响范围将仅局限于基础很小的范围内,此时结构-土-结构的相互作用影响也很小。因此,对于长细比过大或过小的结构,相邻结构的影响都可以忽略不计。当结构的基础宽度0.5 < ${\bar B}$ < 2.5时,需要考虑由相邻结构存在对结构动力特性的影响。图 7表示${\bar \omega }$=2.0时,$\Delta {{\bar \omega }_{\rm{t}}}$随宽厚比的变化情况。对于相邻结构有影响的${\bar L}$ < 2.5范围内,$\Delta {{\bar \omega }_{\rm{t}}}$随宽厚比的变化并非单调函数,在区间${\bar B}$∈[1, 1.5]有峰值的存在。
SSSI体系可作为分析城市环境下结构动力反应的基础。1985年墨西哥地震时,人们首先观察到城市条件下场地地震波中出现异常的长持时和拍的现象[17],并基于结构-土-结构相互作用模型的地震反应分析再现了拍的情况。由前面的分析可知,结构-土-结构相互作用体系存在两阶频率相近但相位相反的模态,这表明建筑结构群下存在引起拍频振动的条件。为此下面以简谐荷载为例,说明结构-土-结构体系在地震作用下拍的现象。
地震时从震源释放出来的能量以地震波的形式在场地内传播,因而不同点的地震动必然存在差异。在当前对基岩地震波的分布规律尚缺乏深入认识的情况下,文中采用行进地震波假定,即土层地震输入为行波输入[18]。在有限元模型中,如果入射起始点的地震波位移时程为ub0(t)时,那么在地震波传播方向上,与入射点相距xi处的水平方向位移时程可以表示为:
式中c为地震波的行进速度。
算例分析两相同单自由度体系明置于一均质土层的SSSI体系,如图 8所示。同时,在有限元模型中仅保留结构1即形成SSI体系。有限元模型中的参数为:土的计算范围为296 m×80 m(长×宽),厚度为16 m,土的弹性模量为8.0×107 Pa,泊松比为0.3,密度为1 600 kg/m3,刚性基础尺寸为16 m×16 m,两结构之间的距离为24 m。SSI体系的基频为1.66 Hz,SSSI体系与之相对应的孪生频率分别为1.64 Hz和1.68 Hz。
上部结构单元阻尼矩阵[C]ste和土单元阻尼矩阵[C]soe分别为:
式中[K]ste和[K]soe分别为结构和土的单元刚度矩阵;ηst和ηso分别为结构和土的阻尼系数。在算例分析中,ηst=0.01,ηso=0.05;ω为激振频率。
为定性研究行波输入情况下,SSSI体系和SSI体系结构动力反应的差别,假定输入地震波为从左向右传播的简谐波,振幅为1 cm,激振频率取略小于基频的1.44 Hz,以突出孪生频率对拍的影响。计算中简谐波的行进速度假定为200 m/s,激振时间为40 s。
SSSI体系和SSI体系结构1的谐振反应如图 9所示。图 10为地表A点的位移反应时程。由图可知:1)对于该简谐荷载和行波波速下,相邻结构对振动初期的瞬态振动有较大影响,对稳态振动影响很小。2)在地震波通过场地后,SSSI体系自由振动的结构和场地都出现拍的现象,但是SSI体系没有拍的现象。与SSI体系的位移时程相比,相邻结构使土层表面自由振动衰减速度降低和持时延长。
在行波输入下,不同激振点输入地震波有时间差(相位差)的存在。地震波的行进速度与下覆土层的剪切波速和入射角度有关,为此比较了Δt/T分别等于0、1/8、1/4、1/2和1情况下,SSSI体系中结构1的动力反应的差别,如图 11所示。图中Δt表示地震波通过相应于两结构水平投影位置的时间差,T表示简谐荷载的周期。计算中取T=0.48 s,其余计算参数同前。Δt/T=0表示行进波速为无穷大,即地震输入为一致输入。两相同结构所组成的SSSI体系中,孪生频率分别为反对称模态和正对称模态,在一致输入作用下无法激起正对称模态的振动,因此,不会因为孪生频率而产生拍,图 11的动力反应也反映了这个现象。在行波输入下,可以同时激发正对称模态和反对称模态,因此,结构1自由振动阶段都出现了拍的现象。由此表明,当结构间距离导致不同结构物之间的输入地震波有明显相位差时,易于产生拍的现象。
采用有限元方法研究了相邻结构对结构-土-结构体系动力特性的影响。系统地讨论了SSSI体系结构动力特性随两结构之间相对距离、频率比以及结构基础宽度与土层厚度比的变化规律。根据计算结果可以得出以下结论:
1) 根据两相同结构下结构-土-结构体系中存在两阶频率相近,相位相反的模态,提出了“孪生频率”的概念。在此基础上,进一步提出“孪生频率相对变化量”指标用以定量描述两相邻结构的动力相互作用。
2) 两结构间的相对距离大于2.5倍基础宽度时,相邻结构对结构动力特性的影响小于5%,可以忽略两结构间的动力相互作用。
3) 结构与土体的频率比${\bar \omega }$ < 0.7时,孪生频率相对变化量小于5%。因此,修建于场地刚度较大的结构,可忽略相邻结构动力相互作用的影响。
4) 对于长细比过大(${\bar B}$ < 0.5)或过小(${\bar B}$ > 2.5)的结构,结构的动力特性由单一结构下土-结构相互作用体系控制,因此,可以忽略相邻结构间的相互影响。
5) 当两结构之间相对距离${\bar L}$ < 2.5,结构与土体的频率比${\bar \omega }$ > 0.7,基础与土层的宽厚比0.5 < ${\bar B}$ < 2.5时,相邻结构之间存在较强的相互作用影响,由此引起的孪生频率相对变化量甚至超过15%。
6) 当结构间距离导致不同结构物之间的输入地震波有明显相位差时,结构-土-结构相互作用体系中存在孪生频率而使结构和邻近场地易于产生拍的现象。