两相邻建筑“结构土结构体系”的动力特性

潘旦光 豆丽萍

摘要：为了研究建筑群中相邻结构存在导致结构动力特性的变化，将上部结构简化为等效单自由度模型，采用刚性基础明置于均质土层上，系统地研究两个相同结构所构成三维结构土结构体系中，邻近结构存在对结构动力特性的影响。根据结构土结构相互作用下体系存在两阶频率相近相位相反振动模态的现象，提出了孪生频率的概念。随后，探讨了孪生频率随结构与土层的频率比、结构之间的相对距离、基础宽度与土层厚度比等因素的变化规律，以反映两结构通过土体耦合对结构动力特性的影响。在此基础上，进行了行波地震输入下结构土结构相互作用体系动力反应计算，结果表明，由于孪生频率的存在，易于使结构和场地产生拍的现象。

关键词：结构土结构相互作用；动力特性；相对距离；频率比；拍

中图分类号：TU3113文献标志码：A文章编号：16744764（2014）03009207

Dynamic Characteristics

of Structure睸oil睸tructure System for Two Neighbor Buildings

Pan Danguang1,Dou Liping1，2

(1.School of Civil and Environmental Engineering, University of Science and Technology Beijing, Beijing 100083, P.R.China;

2.Yinchuan ShiMao Real Estate Development Co. Ltd, Yinchuan 750002, P.R. China)

Abstract: In order to realize the variation of structural dynamic characterics due to neighbor structures in buildings, the surface structure is idealized as an equivalent single degree of freedom system with rigid base whose site consists of a single homogeneous layer. The dynamic characteristics of three瞕imensional structure瞫oil瞫tructure interaction system, including two identical structures, are investigated to identify the additional effects caused by the presence of a second structure. When comparing the data from the soil瞫tructure system with only one structure and the structure瞫oil瞫tructure system, two close natural frequencies with opposite phase modes are identified. Therefore, the term twin瞗requency is proposed. Next, the value of twin瞗requency varying with the frequency ratio of structure to soil is discussed. Meanwhile, the relative distance between two structures and the ratio of foundation width to soil thickness are covered. The data will be used to consider the coupling effects of soil on the dynamic characteristics of structure瞫oil瞫tructure systems. Finally, a numerical case of seismic response of the structure瞫oil瞫tructure system is calculated under traveling wave excitations. The numerical results indicate that the twin瞗requencies cause beatings of the structures as well as ground motion.

Key words: structure瞫oil瞫tructure interaction; dynamic characteristics; relative distance; frequency ratio; beating

為了充分地利用土地，建筑与建筑之间的间距往往很小。日本的Kobe、法国的Nice、美国的Los Angeles、中国的上海等城市部分或全部修建在软土场地上，且位于地震带上或地震带的边缘。在地震作用下，土与结构之间的相互作用(soil瞫tructure interaction, 简称SSI)是影响结构动力反应的关键因素之一［13］。局部区域内的建筑必将通过地基土对相邻结构的动力反应产生影响［45］，由此存在结构土结构相互作用(structure瞫oil瞫tructure interaction, 简称SSSI)问题。

Luco等［6］首先研究了SH波作用下SSSI问题。Wolf［7］对反应堆建筑、反应堆辅助室和燃料处理室组合结构进行地震反应分析。田彼得等［8］采用二维边界元法研究弹性半空上两个条形结构的动力阻抗和动力反应。姜忻良等［9］研究两结构的动力反应随相互之间距离的变化规律。窦立军［10］分析了高层建筑对相邻多层建筑的影响。Padron等［11］研究了桩基结构的SSSI。Clouteau等［12］的计算结果表明，与SSI体系相比SSSI体系可使结构顶层的反应降低。王淮峰等［13］研究了由于相邻结构存在情况下基地剪力的变化情况。

在一定区域范围内的结构通过土体耦合而成为一个有机的整体。结构之间距离和结构基础宽度等的变化相当于改变了系统的质量和刚度分布，由此导致结构动力特性发生相应的变化。而结构的动力特性是影响结构的动力反应最关键的因素之一。为此，本文将上部结构简化为等效单自由度体系，比较两相同结构组成的结构土结构相互作用体系与单一结构下土结构相互作用体系动力特性的变化情况，系统地研究结构土结构体系动力特性随结构间的相对距离、结构与土体相对刚度比和基础宽度的变化规律。

1体系的计算模型

上部结构简化为等效单自由度体系，基础采用明置于等厚度均质水平土层上的方形刚性基础，所构成的结构土结构相互作用体系如图1所示。

图1结构土结构相互作用计算简图

图中B、L、H、D分别表示基础宽度、两基础间净距、土层厚度以及基础边界距人工边界的距离。将土体和结构进行有限元划分后，体系的无阻尼自由振动方程可表示为：

［M］{}+［K］{u}=0(1)

式中，［M］和［K］分别表示质量矩阵和刚度矩阵，{}和{u}分别为节点的加速度和位移向量。本文采用Ansys有限元软件的Block Lanczos法进行模态分析。建模时，上部结构的等效单自由度体系采用带集中质量的梁进行模拟。梁采用beam188单元，集中质量为mass21单元。刚性基础采用shell63单元，土体采用solid45单元，基础与土体无相对滑移。

〖=D(〗潘旦光，等：两相邻建筑“结构土结构体系”的动力特性〖=〗土层实际为一半无限体，当采取用有限元法进行计算时，需要截取一定的范围进行计算，由此带来人工边界的问题［1415］。下面基于SSI体系的有限元模型讨论侧向人工边界范围的选取问题。计算中侧向人工边界为自由边界，底部采用固定边界。SSI体系基频随D/B变化的计算结果如表1所示。从计算结果可以看出，当D/B大于8后，结构基频的计算数值趋于稳定。因此，下面计算中侧向人工边界范围取为D/B=8，以消除侧向人工边界对结构基频的影响。计算中土体的参数如下：弹性模量为20×107 Pa，泊松比为03，密度为1 600 kg/m3，厚度为8 m。刚性基础宽度为8 m，忽略基础的质量。上部结构在刚性基础下的频率为1363 rad/s。此时，结构和土层的频率比=1，基础宽度与土层厚度比=1。表1不同D/B下SSI体系的基频

D/B481216基频/Hz1420 11420 31420 31420 3

从上部结构的计算模型看，刚性基础上的质量弹簧系统可以看作是采用筏基的建筑结构。设筏基的质量为mf，上部结构的质量为m。则不同mf/m下结构的基频如表2所示。表中的误差为采用无质量刚性基础近似计算所得频率的相对误差。显然随着基础质量的增加，采用无质量刚性基础SSI体系计算频率的误差将增加。实际工程中，筏基结构的mf/m一般为5%，此时，采用无质量基础的计算误差小于05%。因此，为简化计算参数，在后面的计算分析中将基础简化为无质量的刚性基础。

表2不同mf/m下SSI体系的基频

mf/m0250%5%750%10%基频/Hz1420 31417 61414 91412 21409 4误差/%00190 50381 70573 60773 4

图2为不同距离下结构土结构相互作用体系自振频率和相应土结构相互作用体系自振频率的计算结果。从计算结果可以看出，考虑土结构相互作用后，结构体系的频率小于刚性基础下结构的频率，这一点与其它相关文献的理论和试验结果相同［16］。但是，对于单自由度下单一结构的土结构相互作用体系而言，结构任一方向水平振动的主模态只有一阶，而与之对应的结构土结构相互作用体系的模态将是两阶频率相近的模态。其中一阶模态为两个结构同向运动，另一阶模态为两个结构相向运动，如图3所示。

图2SSSI体系频率随两结构相对距离的变化

(=100，=100)

圖3SSSI系统模态

结构土结构体系中两阶频率相近模态是两个相邻结构通过土体耦合而产生的特殊的动力特征，具有对称性相反、频率相近的特点。为描述方便，将这样的两阶频率称为“孪生频率”。并定义孪生频率相对变化量为：

Δt=(ωt2－ωt1)/ωsi(2)

式中ωt2和ωt1分别为孪生频率的上限和下限，ωsi为与孪生频率对应的SSI体系的频率。

从图2可以看出：1)结构土结构体系的孪生频率分别位于相应土结构相互作用体系频率的两侧；2)当两个结构相距越远，ωt2和ωt1逐渐收敛到ωsi，这表明随着两个结构彼此之间距离的增加，两结构通过土体的耦合效应逐渐减少，当两者之间的距离达到一定程度后，这种相互作用可以忽略不计。而Δt综合反映了结构土结构相互作用体系相对于土结构体系频率的变化情况，可作为描述结构之间相互影响的特征参数。

下面主要分析两相同结构所组成的结构土结构相互作用系统的频率变化情况。对结构体系频率的影响因素采用以下无量纲参数进行描述：结构和土层的频率比：=ωst/ωso；两结构相对距离：=L/B；宽厚比：=B/H。其中，ωst表示刚性基础下结构的频率，ωso表示仅均质土层的基频。宽厚比是基础宽度与土层厚度比，间接反映上部结构的细长比。

2土结构相互作用体系的频率

结构土结构相互作用对结构动力特性的影响是在土结构相互作用体系动力特性的基础上，由于邻近构筑物存在而发生的变化。为更好的理解孪生频率的变化规律，下面首先分析土结构相互作用体系频率ωsi随频率比和宽厚比的变化情况。计算中土层参数同前，保持不变。图4为改变结构刚度和结构基础宽度情况下土结构相互作用体系频率ωsi的变化规律。

图4土结构相互作用体系的频率

正如所料，增加结构的频率比将减少ωsi/ωst。这表明结构越刚，土层越软，土结构相互作用影响越强。→0，ωsi/ωst→1，土体相当于刚性基础，此时，不存在土结构相互作用。从结构频率的绝对量的大小看，>4以后，ωsi的值基本保持不变，这表明，当频率比大于4后，在其它条件不变的情况下，等效单自由度的频率取决于土层的频率，此时，增加上部结构的刚度对于土结构相互作用体系的动力特性没有影响。

宽厚比反映了结构基础的大小。对于特别细长的结构(如：=0125)，等效单自由度结构的频率将远小于结构自身的频率。此时，存在很明显土结构相互作用是由基础尺寸太小而使土对结构的约束动刚度很小引起的。在其它条件不变条件下，越大，结构的等效单自由度的频率越高，当>25以后，不同基础宽度下的ωsi的值基本相同。这表明，当宽厚比大于25后，单一均匀土层的动刚度对结构动力特性的影响基本不变。

因此，在下面分析结构土结构相互作用体系孪生频率相对变化量Δt时，频率比和宽厚比的计算范围分别取为［02,5］和［0125,4］。

3孪生频率相对变化量的影响因素

图5为不同宽厚比下孪生频率相对变化量Δt随结构间相对距离的变化情况。对于任意宽厚比，随着相对距离增加，孪生频率相对变化量Δt均逐渐减小并趋于零。这是因为相邻结构相互作用是通过地基土来实现的，结构的运动迫使周边的土层发生相应的振动，由此导致邻近结构的振动。两结构间的距离越近，相互作用影响越大。当建筑之间的距离足够大时，相互作用影响很小，可以忽略不计。虽然很大范围内，相邻结构之间的动力特性都存在相互影响，但是模态分析主要作用是为时程分析做准备。在时程分析时，当距离大于一定值后，由于结构和土体阻尼的存在，而使邻近结构导致结构动力特性的改变以及结构有效输入改变对动力反应的影响可以忽略不计。综合其他文献的研究成果［17］，可以发现，当>25后，孪生频率相对变化量Δt小于5%，此时，可以认为结构基本不受相邻结构动力特性的影响。

图5Δt随相对距离的变化

图6Δt随频率比的变化 (L=05)

从频率比的角度看，总的趋势是Δt随频率比的增加而增大，但并非单调变化。为更清楚的反映频率比的影响，图6显示=05情况下，Δt随频率比的变化情况。由图可知：1)当<07时，Δt<5%。这是因为结构的刚度小于土层刚度时，相当于结构修建于较硬的场地上，此时，土结构相互作用本身就很小，因此，通过土体对相邻构筑物的影响就更小，可以忽略不计。2)当>4时，由于土结构相互作用体系的频率取决于土层的频率，相应地Δt的也趋向于一个稳定的值。因此，对于>4的结构体系，结构间的相互影响可以用=4的计算结果近似表达。3)相邻结构对结构动力特性影响最大的区域通常在∈［07,25］，Δt值甚至超过15%。但是其峰值通常在>1时出现，这表明当结构刚度略大于土层刚度时，结构土结构的相互作用影响最大。

从宽厚比的角度看，当<05时，Δt<5%。结合土结构相互作用下频率的变化规律(图4)可知，对于宽厚比小的结构，由于结构本身存在强烈的土结构相互作用，相邻结构对结构动力特性的影响与之相比是一个小量。而当>25时，基础的面积很大，相邻结构通过地基对结构的影响范围将仅局限于基础很小的范围内，此时结构土结构的相互作用影响也很小。因此，对于长细比过大或过小的结构，相邻结构的影响都可以忽略不计。当结构的基础宽度05<<25时，需要考虑由相邻结构存在对结构动力特性的影响。图7表示=20时，Δt随宽厚比的变化情况。对于相邻结构有影响的<25范围内，Δt随宽厚比的变化并非单调函数，在区间∈［1,15］有峰值的存在。

图7Δt随基础宽度的变化(=200)

4行波地震输入下SSSI体系的反应

SSSI体系可作为分析城市环境下结构动力反应的基础。1985年墨西哥地震时，人们首先观察到城市条件下场地地震波中出现异常的长持时和拍的现象［17］，并基于结构土结构相互作用模型的地震反应分析再现了拍的情况。由前面的分析可知，结构土结构相互作用体系存在两阶频率相近但相位相反的模态，这表明建筑结构群下存在引起拍频振动的条件。为此下面以简谐荷载为例，说明结构土结构体系在地震作用下拍的现象。

地震时从震源释放出来的能量以地震波的形式在场地内传播，因而不同点的地震动必然存在差异。在当前对基岩地震波的分布规律尚缺乏深入认识的情况下，文中采用行进地震波假定，即土层地震输入为行波输入［18］。在有限元模型中，如果入射起始点的地震波位移时程为ub0(t)时，那么在地震波传播方向上，与入射点相距xi处的水平方向位移时程可以表示为：

ubi(t)=ub0(t－xi/c)(3)

式中c为地震波的行进速度。

算例分析两相同单自由度体系明置于一均质土层的SSSI体系，如图8所示。同时，在有限元模型中仅保留结构1即形成SSI体系。有限元模型中的参数为：土的计算范围为296 m×80 m（长×宽），厚度为16 m，土的弹性模量为80×107 Pa，泊松比为03，密度为1 600 kg/m3，刚性基础尺寸为16 m×16 m，两结构之间的距离为24 m。SSI体系的基频为166 Hz，SSSI体系与之相对应的孪生频率分别为164 Hz和168 Hz。

上部结构单元阻尼矩阵［C］est和土单元阻尼矩阵［C］eso分别为：

［C］est=ηstω［K］est，［C］eso=ηsoω［K］eso(4)

式中［K］est和［K］eso分别为结构和土的单元刚度矩阵；ηst和ηso分别为结构和土的阻尼系数。在算例分析中，ηst=001，ηso=005；ω为激振频率。

图8结构土结构体系有限元模型

为定性研究行波输入情况下，SSSI体系和SSI体系结构动力反应的差别，假定输入地震波为从左向右传播的简谐波，振幅为1 cm，激振频率取略小于基频的144 Hz，以突出孪生频率对拍的影响。计算中简谐波的行进速度假定为200 m/s，激振时间为40 s。

图9结构1位移响应

图10地表A点位移响应

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=048 s，其余计算参数同前。Δt/T=0表示行进波速为无穷大，即地震输入为一致输入。两相同结构所组成的SSSI体系中，孪生频率分别为反对称模态和正对称模态，在一致输入作用下无法激起正对称模态的振动，因此，不会因为孪生频率而产生拍，图11的动力反应也反映了这个现象。在行波输入下，可以同时激发正对称模态和反对称模态，因此，结构1自由振动阶段都出现了拍的现象。由此表明，当结构间距离导致不同结构物之间的输入地震波有明显相位差时，易于产生拍的现象。

图11不同相位差下SSSI体系中结构1位移响应

5结语

采用有限元方法研究了相邻结构对结构土结构体系动力特性的影响。系统地讨论了SSSI体系结构动力特性随两结构之间相对距离、频率比以及结构基础宽度与土层厚度比的变化规律。根据计算结果可以得出以下结论：

1)根据两相同结构下结构土结构体系中存在两阶频率相近，相位相反的模态，提出了“孪生频率”的概念。在此基础上，进一步提出“孪生频率相对变化量”指标用以定量描述两相邻结构的动力相互作用。

2)两结构间的相对距离大于25倍基础宽度时，相邻结构对结构动力特性的影响小于5%，可以忽略两结构间的动力相互作用。

3)结构与土体的频率比<07时，孪生频率相对变化量小于5%。因此，修建于场地刚度较大的结构，可忽略相邻结构动力相互作用的影响。

4)对于长细比过大(<05)或过小(>25)的结构，结构的动力特性由单一结构下土结构相互作用体系控制，因此，可以忽略相邻结构间的相互影响。

5)当两结构之间相对距离<25，结构与土体的频率比>07，基础与土层的宽厚比05<<25时，相邻结构之间存在较强的相互作用影响，由此引起的孪生频率相对变化量甚至超過15%。

6)当结构间距离导致不同结构物之间的输入地震波有明显相位差时，结构土结构相互作用体系中存在孪生频率而使结构和邻近场地易于产生拍的现象。

参考文献：

［1］Paul S J. Design aids for simplified nonlinear soil瞫tructure interaction analyses ［J］. Engineering Structures, 2012, 34: 572580.

［2］ Farshidianfar A, Soheili S. Ant colony optimization of tuned mass dampers for earthquake oscillations of high瞨ise structures including soil瞫tructure interaction ［J］. Soil Dynamics and Earthquake Engineering, 2013, 51:1422.

［3］ Jarernprasert S, Bazan瞆urita E, Bielak J. Seismic soil瞫tructure interaction response of inelastic structures [J]. Soil Dynamics and Earthquake Engineering, 2013, 47: 132143.

［4］ Ghergu M, Ionescu I R. Structure瞫oil瞫tructure coupling in seismic excitation and ‘city effect ［J］. International Journal of Engineering Science, 2009, 47(3): 342354.

［5］ Lou M, Wang H, Chen X, et al. Structure瞫oil瞫tructure interaction: literature review ［J］. Soil Dynamics and Earthquake Engineering, 2011, 31(12): 17241731.

［6］ Luco J E, Contesse L. Dynamic structure瞫oil瞫tructure interaction ［J］. Bulletin of the Seismological Society of America, 1973, 63(4): 12891303.

［7］ Wolf J P. Dynamic soil瞫tructure interaction ［M］. Englewood Cliffs, N. J.: Prentice睭all, 1985.

［8］田彼得，俞载道. 结构土结构相互作用体系的动力分析［J］. 同济大学学报:自然科学版, 1987, 15(2): 1627

Tian B D, Yu Z D. Dynamic analysis of structure瞫oil瞫tructure interaction ［J］. Journal of Tongji University: Natural Science, 1987, 15(2): 1627.

［9］姜忻良，严宗达，武志金. 相邻结构地基土相互作用的分支模态实用研究 ［J］. 地震工程与工程振动, 1998, 18(3): 914

Jiang X L, Yan Z D, Wu Z J. Research on practical method of branched mode for neighboring structures foundation soil interaction system ［J］. Journal of Earthquake Engineering and Engineering Vibration, 1998, 18(3): 813.

［10］窦立军，杨柏坡. 高层建筑与相邻多层建筑间的动力相互作用［J］. 地震工程与工程振动, 2000, 20(3): 1521

Dou L J, Yang B P. Dynamic interaction of tall building with neighboring multi瞫toried building ［J］. Journal of Earthquake Engineering and Engineering Vibration, 2000, 20(3): 1521.

［11］ Padron L A, Aznarez J J, Maeso O. Dynamic structure瞫oil瞫tructure interaction between nearby piled buildings under seismic excitation by BEM睩EM model ［J］. Soil Dynamics and Earthquake Engineering, 2009, 29(6):10841096.

［12］ Clouteau D, Broc D, Devesa G, et al. Calculation methods of Structure睸oil睸tructure Interaction (3SI) for embedded buildings: application to NUPEC tests ［J］. Soil Dynamics and Earthquake Engineering, 2012, 32 (1) 129142.

［13］王淮峰，楼梦麟，陈希,等. 建筑群结构土结构相互作用的影响参数研究［J］. 同济大学学报:自然科学版, 2013, 41(4): 510514

Wang H F, Lou ML, Chen X, et al. Parametric study on structure瞫oil瞫tructure interaction of high瞨ise buildings ［J］. Journal of Tongji University: Natural Science, 2013, 41(4): 510514.

［14］潘旦光，楼梦麟，董聪. 一致输入作用下土层的地震反应分析［J］. 计算力学学报, 2005, 22(5): 562567

Pan D G, Lou M L, Dong C. Seismic response analysis of soil layer under uniform excitation［J］. Chinese Journal of Computational Mechanics, 2005, 22(5):562567.

［15］李煌. 地基与高层建筑结构动力相互作用分析［D］. 南宁：广西大学，2006

［16］李辉，賴明，白绍良. 土结动力相互作用研究综述(Ⅰ)——研究的历史、现状与展望［J］. 重庆建筑大学学报, 1999, 21(4): 112116

Li H, Lai M, Bai S L. A Summary of research on soil瞫tructure dynamic interaction (I) 勃瞙istory, status in quo and prospects［J］. Journal of Chongqing Jianzhu University, 1999, 21(4):112116.

［17］ Groby J P, Wirgin A. Seismic motion in urban sites consisting of blocks in welded contact with a soft layer overlying a hard half瞫pace ［J］. Geophysical Journal International, 2008, 172(2): 725758.

［18］潘旦光，楼梦麟，董聪. 土层地震行波反应分析中侧向人工边界的影响［J］.岩土工程学报,2005，27(3): 308312

Pan D G, Lou M L, Dong C. Effect of vertical artifical boundary on seismic response of soil layer under traveling wave excitations ［J］. Chinese Journal of Geotechnical Engineering, 2005, 27(3): 308312.