Research ArticleEffect of Soil-Structure Interaction on Seismic Performance ofLong-Span Bridge Tested by Dynamic Substructuring Method
Zhenyun Tang, HuaMa, Jun Guo, and Zhenbao Li
The Key Laboratory of Urban Security and Disaster Engineering, Ministry of Education,Beijing University of Technology, Beijing 100124, China
Correspondence should be addressed to Hua Ma; [email protected]
Received 9 October 2016; Revised 13 December 2016; Accepted 19 January 2017; Published 8 March 2017
Academic Editor: Salvatore Strano
Copyright Β© 2017 Zhenyun Tang et al. This is an open access article distributed under the Creative Commons Attribution License,which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Because of the limitations of testing facilities and techniques, the seismic performance of soil-structure interaction (SSI) system canonly be tested in a quite small scale model in laboratory. Especially for long-span bridge, a smaller tested model is required whenSSI phenomenon is considered in the physical test. The scale effect resulting from the small scale model is always coupled withthe dynamic performance, so that the seismic performance of bridge considering SSI effect cannot be uncovered accurately by thetraditional testing method. This paper presented the implementation of real-time dynamic substructuring (RTDS), involving thecombined use of shake table array and computational engines for the seismic simulation of SSI. In RTDS system, the bridge withsoil-foundation system is divided into physical and numerical substructures, in which the bridge is seen as physical substructuresand the remaining part is seen as numerical substructures.The interface response between the physical and numerical substructuresis imposed by shake table and resulting reaction force is fed back to the computational engine. The unique aspect of the methodis to simulate the SSI systems subjected to multisupport excitation in terms of a larger physical model. The substructuring strategyand the control performance associated with the real-time substructuring testing for SSI were performed. And the influence of SSIon a long-span bridge was tested by this novel testing method.
1. Introduction
In the seismic analysis of a structure founded on ground, theground motion passes to the base of structure and then loadson structure. The response of the foundation system affectsthe response of the structure and vice versa, which is calleddynamical soil-structure interaction (SSI).Theoretical results[1, 2] indicate that SSI is sometimes beneficial and sometimesdetrimental to structural performance. Therefore, the effectof soil cannot be neglected, because SSI phenomenon isclosely related to its dynamic characteristics [3], especiallythe damping [4] of the whole system. To explore the effectof SSI on the seismic performance of engineering structures,the finite element methods [5] and theoretical analysis [6] areoften used. However, the uncertainties [7, 8] and boundarycondition [9] existing in SSI system have not been simulatedproperly by thesemethods.At present, experimental evidencerelating to SSI system is scarce. The SSI system is difficult tobe investigated by testing full-size specimen (including both
structure and soil-foundation systems) under earthquakeloading because of the size and power of the testing facil-ities, especially, for the structures with large spatial extent.In current testing of SSI, the soil-foundation systems arereplaced by laminar shear box with soil [10β12], and thestructures are usually scaled down to small size (often scaledto 1/30 or even smaller) or a simple cantilevered mass [10, 11],which results in the inevitable possibility of size effects. Inaddition, time must be scaled, accompanying with frequencyscaled up, in test to provide rapid shaking, often preventedby the frequency bandwidth of shaking table. SSI plays amore important role in the seismic performance of long-spanbridge [13]. However, the smaller specimen is required in thelaboratory experiment because of its dimension. As a result,simplified models are tested at reduced scales making theinterpretation and application of results challenging.
The method of real-time substructuring derived fromhardware-in-the-loop test [14] has been paid close attentionsince the first real-time substructuring test was reported by
HindawiShock and VibrationVolume 2017, Article ID 4358081, 12 pageshttps://doi.org/10.1155/2017/4358081
2 Shock and Vibration
Numerical modelEarthquake
loading Shaking tableInterface response
Physical model Interface response
Interface reaction
Figure 1: Generally schematic representation of RTDS using shaking table.
Nakashima et al. [15], which is a method of dynamicallytesting a structure without experimentally testing a physicalmodel of the entire system; instead the structure can besplit into two coupled parts, the region of particular interest,which is tested experimentally, and the remainder which istested numerically. It also received interest in SSI systems.To simulate the soil medium numerically in a shakingtable substructuring test where the superstructure is testedexperimentally, Konagai et al. [16, 17] used analogue electriccircuits to generate the transient response of both linearand nonlinear soils. Heath et al. [18] and Wang et al. [19]investigated an inverted substructuring system by modellingthe superstructure numerically. However, only very simplesoil model and superstructures were adopted in these works.To simulate the SSI phenomenon with bridge experimentallyusing RTDS, further validation is required involving morecomplex structures and more complex interactions.
In the research presented in this paper, a framework ofsubstructuring for SSI system with bridge was establishedfirstly; after that the control strategies for the real-timedynamic substructuring testing were analysed. Finally, theinfluence of SSI on a long-span bridge was tested through thisnovel testing method.
2. Modelling of RTDS Systemsfor SSI Simulation
2.1. RTDS Scheme Based on Shaking Table Array. The devel-opment of advanced computation and control techniquesresulted in the great progress of real-time substructure[20, 21]. The basis of such tests is to divide the overallsystem into experimental test specimen and the numericalmodel. Typically the experimental test specimen containsthe area of particular interest which behaves nonlinearly anduncertainly; the numerical model contains the remainderof the system which behaves linearly or can be modelledadequately. These two parts are tested in parallel and inreal-time with information relating to the interface betweenthe two parts being exchanged between them, allowing thetwo parts to emulate the whole system. Figure 1 shows thegenerally schematic representation of substructuring used toemulate dynamics of structures in earthquake engineering.Interface responses are passed from the numerical model toshaking table, which reproduces accurately the responses atthe interface of the test specimen.The interface reaction force
measured from test specimen then feed back to the numericalmodel.
Importantly, the most significant advantage of real-timesubstructuring is that the process is implemented as a time-stepping routine in real-time, which allows that not onlycan the nonlinear and uncertain behaviour of physical sub-structure be accurately tested, but also the time-dependentnonlinearity of numerical substructure can be adopted in thetesting. The test method mentioned above supplies a novelsolution for the SSI testing. As known, compared with soil-foundation system, the behaviour of structure is relativelyeasily modelled, especially, the nonlinearity including gap-ping between the foundations and soil and nonlinear stress-strain response of the soil, and so on. However, the realinterest in SSI is the response of structure under earthquakeloading; meanwhile, the complexity of structural form andmaterial, geometrical, and contacting nonlinearity of struc-ture are also not well described theoretically. Therefore, boththe two kinds of substructuring (soil-foundation system orstructure modelled numerically) for SSI need to be devel-oped. In this work, the soil-foundation system is treated asnumerical substructure and simulated in computer, while thebridge seems as the physical substructure and tested in thelaboratory. A diagrammatic representation of RTDS for abridge is shown in Figure 2. The emulated SSI system shownin Figure 2(a) can be tested by the RTDS strategy shown inFigure 2(b). Unlike other building structures, a single shakingtable is enough to carry out the test; the length of bridgeis quite large, even the scale model is often more than tenmeters long. To overcome this drawback, the shaking tablearray composed of multishaking tables was constructed andused to test the seismic performance of bridge [22].Therefore,the shaking table array is adopted to reproduce the interfaceresponse of each bridge pier and foundation.
2.2. Numerical Model of Soil-Foundation System. The groundmotion imposing on superstructure caused by earthquakeincludes three portions (as illustrated in Figure 3(a)): themovement of free-field soil resulting from earthquake (π’g),the deviation of the foundation movement resulting fromthe free-field soil motion (π’gs), and the deformation offoundation produced by the reaction force from superstruc-ture (inertia interaction: πΏ and π). This work focuses onthe interaction between soil and structure; consequently,the foundation movement is identical to the free-field soilmovement (π’gs = 0). In the seismic analysis for SSI, the
Shock and Vibration 3
(a) Whole soil-structure system under multisupport excitation
Reaction force: F, M
Physical substructure
Numerical substructure
Interface responses: οΏ½νΏ, οΏ½ν
kh, ch
kr, cr
(b) Substructured system:soil-foundation system represented as numerical substructure
Figure 2: Substructuring system for bridge subject to earthquake excitation.
Semi-infinite
FM
οΏ½νΏ
οΏ½ν
ug + ugs
ug
(a) Deformation of soil-pile-foundation
FoundationFd
Kb Cb
(b) Time-domain difference model
Figure 3: Numerical model of soil-pile-foundation system.
substructure method is cost-effective, comparing with thedirect integrity approach, in which the behaviour of the soil-foundation system (see in Figure 3(a)) was described as animpedance function (also shown in Figure 3(b)):
πΎh = πh (π) + πππh (π) , (1)
where πh is horizontal stiffness, πh is horizontal damping, πis circular frequency, and π is imaginary unit. However, theimpedance function is frequency-dependent, which cannotbe used to represent the nonlinear effect of dynamic interac-tion between soil and structure.
In order to overcome this disadvantage, two main meth-ods [23β25] were developed to transform the frequency-dependent impedance function to time-domain: lumpedparameter model approximately fitting to impedance func-tions using different mass-spring-damping systems and fastintegral algorithms (also referred to as recursivemodel) basedon the properties of Fourier inverse transform. Nevertheless,the two methods have their own defects: lumped parametermodel [23, 24] can only fit the actual impedance functions
in narrow frequency band accurately but cannot reflect theregular component, which corresponds to time-delay effect;meanwhile, recursive model [25] has intrinsically limitationat Nyquist frequency; namely, the imaginary part of filtermust be zero at Nyquist frequency.
So as to take the full advantage of the two models andmake up for their deficiency as well, a time-domain differencemodel (TDDM) for soil-foundation system was proposed byDu and Zhao [26], combining lumped parameter model andtime-domain recursive model, which can take into accountboth the singular and regular component of foundationimpedance comprehensively.Themechanicalmodel is shownin Figure 3(b), where
πΎπ = ππΎπ π1,πΆπ = ππΎπ π1,
πΉπ =π
βπ=1
πππ¦π (π‘ β πΞπ‘) +π
βπ=1
πππΉπ (π‘ β πΞπ‘)(2)
4 Shock and Vibration
Table 1: Parameters of the used TDD model.
π1 π2 π3 π4 π5 π6 π1 π11.198 β0.013 β0.415 β0.127 0.02 0.004 0.659 0.011π1 π2 π3 π4 π5 π6 β ββ0.681 β0.205 0.433 0.028 β0.029 β0.003 β β
with π being numbers of piles, πΎπ static stiffness of singlepile, π1 normalized stiffness, π1 normalized damping, πΉπ(π‘)pseudo-force, π¦π(π‘) displacement of foundation, Ξπ‘ integra-tion time step, and ππ, ππ the coefficients of pseudo-force anddisplacement. From (2) we can see that the determinationof parameters π1, π1, ππ, ππ is the key point for the time-domain differencemodel, the best values of which are derivedusing hybrid optimization based on the theory or numericalsolution of impedance function.
Herein, we assume that the bridge pier is constructed on a3Γ 3 pile group foundation, and the parameters are as follows:damping ratio of soil is 0.05, Poissonβs ratio of soil is 0.4, ratioof pile spacing and diameter is 5, YoungModulus ratio of pileto soil is 1000, mass density ratio of pile to soil is 1.42, andratio of pile length and diameter is 15. Normalized by singlepile static stiffness, the impedance function of the pile groupfoundation can be expressed as
πΎh (π0) = 9πsh [πh (π0) + ππh (π0)] , (3)
where π0 = ππ/Vπ is dimensionless frequency, πsh is singlepile horizontal static stiffness, Vπ is shear wave velocity of soil,π is pile diameter. Using the dynamic impedance functiondeveloped by Makris and Gazetas [27, 28] for the pile groupfoundation, the parameters of TDDM were estimated bymultiple regression analysis. The evaluated values of allparameters for this model are shown in Table 1. The com-parative results of the TDDM and Markrisβs analytical modelpresented in Figure 4 show that TDDMis suitable to representthe soil-foundation system in RTDS.
3. Control Strategies for RTDS
3.1. RTDS Controller Design. In commonly used RTDS sys-tem, displacement or force was used to be target signalto drive transfer system; however in some SSI systems,such as sandy soil, the effect of SSI on the accelerationof structure is noticeable, but the effect on displacementis inconspicuous. In that case, the acceleration of interfacecannot be reproduced well by using displacement control. Inthis work, the acceleration driving was attempted to be usedin the RTDS for SSI of this work. In this RTDS system,shaking table array was used as the transfer system toproduce the interface response of foundation, resulting fromthe dynamics of which, the response between physical andnumerical substructure is always asynchronous, even causingthe instability, especially for the low damping system [29].In standard earthquake testing, the shaking table is drivenby a conventional controller, such as PID or three-variablecontroller. These controllers cannot give a perfect control.Phase lag and magnitude error always exist. In order to
Real part (Makrisβs result)Imaginary part (Makrisβs result)
Real part (TDDM)Imaginary part (TDDM)
Kh/9k
sh
a0
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
0.0 0.2 0.4 0.6 0.8
Figure 4: Normalized dynamic stiffness of pile group foundationreported by Du and Zhao [34].
impose the interface response calculated from numericalsubstructure on the physical substructure through shakingtable synchronously, an additional controller is necessary tocancel the shaking table dynamics.
Based on TDDM presented in Section 2.2, the emulatedsystem shown in Figure 5(a) can be substructured as theRTDS system shown in Figure 5(b). TDDMmodel representsthe soil-foundation system.The bridge is installed on shakingtable array, which is used as the transfer system to simulate theinterface acceleration (ππ). As can be seen in Figure 5(b), theinterface acceleration resulted from external excitation (π’g)and reaction force of superstructure (π):
ππ = πΊπππ’g + πΊπππ. (4)
From Figure 5(b), the transfer function models for thedifferent parts can be determined.
πΊππ (π )
= π 2πππππ 2 + πππ β (π1π + π1 + β6π=1 πππ§βπ/ (1 + β6π=1 πππ§βπ))
,
πΊππ (π )
= π 2πππ 2 + πππ β (π1π + π1 + β6π=1 πππ§βπ/ (1 + β6π=1 πππ§βπ))
,
(5)
Shock and Vibration 5
Shaking tableactuator
Ground motion
Bridge
Soilbox
ug
(a) Emulated system
Inner-loopcontroller Shaking table
actuator
RTDScontroller
u
eKb
Cb
Mp
Mb
ap
an
an
Kp , Cp
β
+
f
Ξ£p
Ξ£Nug
(b) Substructured system
Figure 5: Framework of RTDS for SSI simulation.
where π§ = πππΞπ‘, Ξπ‘ is integration time step. Normally, theachieved value of interface acceleration is not ππ but ππ aftergoing through shaking table because of its dynamics. Thatleads to the needs of RTDS controller. As reference [30]mentioned, the shaking table can be approximated to a first-order system such that
πΊπ (π ) = ππ + π , (6)
where π is the product of the time constant and the propor-tional gain of shaking table. And then the reaction force isformulated as
π = πΊπππ. (7)
πΊπ represents the transfer function of the total shear force (π)of superstructure from excitation (ππ), which is determinedby the dynamic parameters. From Figure 5(b), πΊπ can beexpressed as
πΊπ =πΆππ + πΎπ
πππ 2 + πΆππ + πΎπ , (8)
where Mπ, Cπ, and Kπ are the mass, damping, and stiffnessof bridge, respectively.
To formulate the synthesis procedure of RTDS dynamicsand control, a framework of linear substructuring controller(LSC) [31] together with the minimal control synthesis witherror feedback (MCSEF) algorithm [32] was proposed. Tosimplify the LSC for multidegree of freedom system, a mod-ified linear substructuring controller (MLSC) together withMCSEF was developed by Guo et al. [33]. As illustrated inFigure 6, an emulated system is conceptually decomposedinto at least two substructures, Ξ£π and Ξ£π. For practical casesΞ£π represents the numerical model and Ξ£π represents the
+
+
+
+
MCSEF
ug
ap
an
u
e
β
β
Kr
Ke
G1
G2
G
Ξ£p
Ξ£N
Figure 6: MLSC-MCSEF controlled RTDS system.
physical model. Control and excitation signals are denoted byπ’ and π, respectively. Here, the RTDS dynamics are rep-resented by three generalized blocks, πΊ, constraint; πΊ1,excitation; and πΊ2, transfer system dynamics. πΊ1(π ) is therelationship between the displacements of the top of thefoundation (πππ) and the ground motion (π’g); πΊ(π ) is therelationship between the control signals (π’) and the accel-eration of the top of the foundation (πππ) from the reactionforce, including the TS dynamics πΊπ (π ), the reaction forcedynamics πΊπ(π ), and the numerical substructure dynamicsπΊππ(π ); πΊ2(π ) is the TS dynamics; thus
πΊ1 (π ) = πΊππ (π ) ,πΊ (π ) = πΊπ (π ) πΊπ (π ) πΊππ (π ) ,πΊ2 (π ) = πΊπ (π ) .
(9)
6 Shock and Vibration
With reference to [31], the error dynamics (π = ππ β ππ) of anLSC-controlled DSS can be written as
πΊerr (π ) = πΊ1 (π ) β [πΊ (π ) + πΊ2 (π )]πΎπ (π )1 + [πΊ (π ) + πΊ2 (π )]πΎπ . (10)
To synchronize the two outputs {ππ, ππ}, the substructuringerror is expected to equal zero, the numerator of (10) is setequal to zero, and then
πΎπ (π ) = πΊ1 (π )πΊ (π ) + πΊ2 (π ) . (11)
The control signal and adaptive gains of MCSEF in Figure 6are generated from (12), where {πΌ, π½} are adaptive weights;the ratio πΌ = 10π½, which has been shown to work wellempirically [29], is used here.The term π¦π is the output error,generated directly from π, according to π¦π(π‘) = πΆππ(π‘), whereπΆπ is selected to ensure a strictly positive real dynamic in thehyperstability proof for the MCSEF controller [29]; for first-order control,πΆπ is normally defined as [29]:πΆπ = 4/π‘π , whereπ‘π is the step response time of the implicit reference model.
π’ (π‘) = πΎπ (π ) οΏ½οΏ½g (π‘) + πΎπ (π ) π (π ) + πΎππ (π‘) οΏ½οΏ½g (π‘)+ πΎππ (π‘) π (π‘) ,
πΎππ (π‘) = πΌβ«π‘
0π¦π (π‘) ππ (π‘) ππ‘ + π½π¦π (π‘) ππ (π‘) ππ‘,
πΎππ (π‘) = πΌβ«π‘
0π¦π (π‘) ππ (π‘) ππ‘ + π½π¦π (π‘) ππ (π‘) ππ‘.
(12)
3.2. Validation of RTDSController Performance. For verifyingthe validity of the RTDS for SSI simulation, further tests ofa one-story steel frame including SSI considerations wereperformed. The parameters for the soil model are shown inTable 1, and πsh = 62600N/m.The parameters of the physicalsubstructure are as follows: ππ = 115 kg, πΆπ = 17.67N/ms,and πΎπ = 7.56 Γ 104N/m. The shaking table used here is ashake table array including nine shaking tables (shown in Fig-ure 7) fromBeijing University of Technology (BJUT), and thespecifications are shown in Table 2. The estimated parameterof the shaking table used in MLSC design was π = 16Hz. Theacceleration record of El Centro earthquake (NS componentof the 1940 El Centro, Station: 117 El Centro Array #9,Component 270β) was chosen as the ground motion. Exper-iment and simulation results are shown in Figure 8, whichdemonstrate that the RTDS method emulates dynamic soil-structure interaction effectively. That means the MLSC-MCSEF controller compensated the dynamics of shakingtable successfully.
4. The Effect of SSI on Seismic Performance ofa Long-Span Bridge
4.1. Test Design. Thedynamic response of a large-span bridgewith the considerations of SSI subjected to earthquake excita-tion was tested using RTDS method developed in this work.
Table 2: Specifications of shake table array in BJUT.
Table size 1m Γ 1m (π/π),B600 (π)Operation mode π/π/πMaximum specimen mass 5 tonVelocity Β±80 cm/sFrequency of operation 0.1βΌ50HzTable mass 1 tonNumber 9Displacement Β±7.5 cmAcceleration (full load) Β±1.5 g (π/π), 0.8 g (π)Control mode TVC
Figure 7: Shake table array in BJUT.
In this RTDS test, the specimen (see in Figure 9) was takenas physical substructure.The parameters of this specimen aremass 4500 kg, length 13.84m, and height 1.2m. The time-domain difference model described in Section 2 was chosenas numerical substructure. The parameters of soil model arelisted in Table 1. And the Young modulus of pile is 113.4Gpa,the density of soil is 2000 kg/m3, and diameter of pile is 0.4m.The photo of the specimen is shown in Figure 10. Herein, fourshaking tables are used to produce the interface response ofeach foundation under piper. The acceleration record of ElCentro (See in Figure 11(a)) and Wenchuan earthquake (Seein Figure 11(b)) was chosen as the ground motion.
Before conducting RTDS test, a conventional shakingtable test with white noise excitation was carried on toidentify the dynamic parameters of the physical substructure.The Fourier spectrums of the acceleration measured at pier2were shown in Figure 12. As can be seen, the frequencies of thefirst and second mode were 4.76Hz and 9.58Hz, and thefirst mode is the dominated mode. For the design of MLSCcontroller in this test, the physical substructure still can betreated as a single degree-of-freedom system. The dampingratio of the first mode obtained from Figure 12 is 2.6%.
4.2. Effect of SSI on the Response of the Long-Span Bridge. Toevaluate the effect of SSI on the dynamic response of long-span bridge, the soft soil-foundation system with the shearwave velocity Vπ = 100m/s, 200m/s, and 400m/s was testedusing the MLSC-MCSEF controlled RTDS. For comparison,
Shock and Vibration 7
0 1 2 3 4 5 6 7 8 9
0
1
2
Time (sec.)
TestingSimulation
β1
β2
Acc.(m
2/s)
Figure 8: Acceleration of physical substructure measured from the RTDS for one-story frame.
Table1 Table2 Table3 Table4
Pier1 Pier2 Pier3 Pier4
3.68m 6.48m 3.68m
1.2mug
Figure 9: Dimensions of specimen.
Figure 10: Photo of specimen.
the hard soil-foundation system was also tested using theconventional shaking table test without the consideration ofSSI.
In the experiment, the MLSC-MCSEF strategy wasadopted. Figure 13 displayed the desired and achieved accel-eration of shaking table in the form of synchronization sub-plots. It was seen that MLSC-MCSEF supplied an acceptableaccuracy for RTDS-SSI test.
The tested results of acceleration and strain responseexcited by Wenchuan earthquake were presented in Figures14(a) and 14(b), respectively. As can be seen from Figure 14,the SSI effect enlarged the dynamic response of bridge whenthe shear wave velocity of soil is 400m/s and 200m/s. Theamplification ratios of acceleration for Vπ = 400m/s and200m/s were around 45% and 30%, respectively, while thecorresponding value of strain was around 106% and 86%.However, when the shear wave velocity of soil is too low (e.g.,Vπ = 100m/s), the SSI effect reduced the dynamic response ofbridge in this case. The acceleration and strain were bothreduced to 50% of hard soil case (Vπ = infinity). The clearevidence can be obtained from the acceleration measured on
the shaking table shown in Figure 15. Compared with thehard soil case (Vπ = infinity), the excitation imposed on thebridge was changed significantly by soil. Generally, the highfrequency components of the seismic wave were filtered bysoil, while the low frequency components were amplified.Thecan be seen clearly from the time history (see Figure 15(a))and Fourier spectrum (see Figure 15(b)). In the hard soil case,the seismic wave used in this test has abundant componentsfrom 4Hz to 16Hz. After considering SSI, the componentsof the frequency more than 10Hz were lessened, and thecomponents of the frequency less than 8Hz were enhancedin the cases of Vπ = 400m/s and 200m/s. However, the com-ponents of the low and high frequency were both filtered inthe case of Vπ = 100m/s. In this case, the soil is too soft andequivalent to an isolator.When Vπ = 400m/s and 200m/s, theamplified components were close to the frequency of the firstmode of the bridge; hence, the acceleration and the strain ofthe bridge were enlarged.
Table 3 summarized themaximum strain at the bottom ofPier2 and Pier3 measured from El Centro and Wenchuanearthquake. It can be seen that Pier2 and Pier3 gave the
8 Shock and Vibration
0 10 20 30 40
0
0.05
0.1
0.15
Time (sec.)
Acc.
(g)
β0.05
β0.1
β0.15
β0.2
(a) El Centro
Time (sec.)
Acc.
(g)
0 10 20 30 40 50 60
0
0.05
0.1
0.15
β0.05
β0.1
β0.15
β0.2
(b) Wenchuan
Figure 11: Seismic wave used in the test.
0 5 10 15 200
0.5
1
1.5
2
2.5
3
3.5
Frequency (Hz)
Γ10β4
4.76Hz
9.58HzFour
ier s
pect
rum
(gΒ·s)
Figure 12: Flourier spectrum of the acceleration at the top of pier2.
0.05 0.15 0.25
0.05
0.15
0.25
β0.05
β0.15
β0.25β0.05β0.15β0.25
an (g)
a p(g
)
Figure 13: Control performance of the table below pier2.
Shock and Vibration 9
0
0.5
1
Acc.
(g)
0
0.5
1
Acc.
(g)
0
0.5
1
Acc.
(g)
0 5 10 15
0
0.5
1
Acc.
(g)
Time (sec.)
0 5 10 15Time (sec.)
0 5 10 15Time (sec.)
0 5 10 15Time (sec.)
οΏ½s = infinity
οΏ½s = 400 m/s
οΏ½s = 200 m/s
οΏ½s = 100 m/s
β0.5
β1
β0.5
β1
β0.5
β1
β0.5
β1
(a) Acceleration at the top of pier2
0
400
800
0
400
800
0
400
800
0 5 10 15
0
400
800
Time (sec.)
0 5 10 15Time (sec.)
0 5 10 15Time (sec.)
0 5 10 15Time (sec.)
οΏ½s = infinity
οΏ½s = 400 m/s
οΏ½s = 200 m/s
οΏ½s = 100 m/s
Stra
in (οΏ½ν
οΏ½ν)St
rain
(οΏ½νοΏ½ν)
Stra
in (οΏ½ν
οΏ½ν)St
rain
(οΏ½νοΏ½ν)
β400
β800
β400
β800
β400
β800
β400
β800
(b) Strain at the bottom of pier2
Figure 14: Dynamic response of bridge.
Table 3: Maximum strain at the bottom of Pier2 and Pier3.
Vπ (m/s) Infinity 400 200 100
El Centro Pier2 (ππ) 269 559 601 272Pier3 (ππ) 193 416 505 216
Wenchuan Pier2 (ππ) 345 710 643 185Pier3 (ππ) 289 725 513 118
similar conclusion. Whereas almost the same reductionratio was measured from El Centro excitation when Vπ =100m/s, the amplification ratio is smaller than Wenchuan
excitation when Vπ = 200m/s and 400m/s. This resultedfrom the frequency characteristics of the two excitations.TheFlourier spectrum plotted in Figure 16 demonstrated that ElCentro excitation (see Figure 16(a)) has different frequencycomponents with Wenchuan excitation (see Figure 16(b)). Itshowed that the characteristics of excitation determined theresponse of bridge together with the soil properties in SSIsystem.
5. Conclusions
The purpose of this work is to establish a potential test-ing method for the simulation of soil-structure interaction
10 Shock and Vibration
0 5 10 15 20Time (sec.)
0 5 10 15 20Time (sec.)
0 5 10 15 20Time (sec.)
0 5 10 15 20Time (sec.)
οΏ½s = infinity
οΏ½s = 400 m/s
οΏ½s = 200 m/s
οΏ½s = 100 m/s
0
0.15
0.3
β0.15
β0.3
Acc.
(g)
0
0.15
0.3
β0.15
β0.3
Acc.
(g)
0
0.15
0.3
β0.15
β0.3
Acc.
(g)
0
0.15
0.3
β0.15
β0.3
Acc.
(g)
(a) Time history
0 5 10 15 20 25Frequency (Hz)
0 5 10 15 20 25Frequency (Hz)
0 5 10 15 20 25Frequency (Hz)
0 5 10 15 20 25Frequency (Hz)
οΏ½s = infinity
οΏ½s = 400 m/s
οΏ½s = 200 m/s
οΏ½s = 100 m/s
0
1
2
Four
ier s
pec.
Γ10β3
0
1
2
Four
ier s
pec.
Γ10β3
0
1
2Fo
urie
r spe
c.Γ10β3
0
1
2
Four
ier s
pec.
Γ10β3
(b) Fourier spectrum
Figure 15: Acceleration response of the table below pier2.
based on real-time dynamic substructuring (RTDS), testingbridge experimentally and modelling the remainder partnumerically, which make the long-span bridge testing withconsideration of SSI possible. The βsize effectβ resulting fromthe scaled model in whole system tests because of the sizeand power of the testing facilities is weakened; meanwhilethe nonlinear dynamics of the SSI system can be investigatedusing this method. The feasibility of this method has beenverified experimentally.
The influence of SSI on a long-span bridge was testedthrough RTDS method. The SSI effect enlarged the dynamicresponse of bridge when the shear wave velocity of soil is400m/s and 200m/s, while the SSI effect reduced thedynamic response of bridge when the shear wave velocity of
soil is too low (e.g., Vπ = 100m/s). In other words, themediumsoft soil-foundation may be detrimental to the seismicperformance of bridge, and the soft soil may be beneficialto the seismic performance based on the test results of thiswork.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the support of theNational Science Foundation of China, Grant no. 51608016,
Shock and Vibration 11
0 5 10 15 20 250
0.5
1
1.5
2
2.5
3
3.5
Frequency (Hz)
Γ10β4Fo
urie
r spe
ctru
m (g
Β·s)
(a) El Centro
0 5 10 15 20 250
0.2
0.4
0.6
0.8
1
1.2
1.4
Frequency (Hz)
Γ10β4
Four
ier s
pect
rum
(gΒ·s)
(b) Wenchuan
Figure 16: Flourier spectrum of the seismic wave used in the test.
and the support of the Beijing Natural Science Foundation(8164050) in the pursuance of this work.
References
[1] J. P. Wolf, Dynamic Soil-Structure Interaction, Prentice-Hall,Englewood Cliffs, NJ, USA, 1985.
[2] G. Mylonakis and G. Gazetas, βSeismic soil-structure interac-tion: beneficial or detrimental?β Journal of Earthquake Engi-neering, vol. 4, no. 3, pp. 277β301, 2000.
[3] M.-L. Lou, W.-J. Wang, H.-C. Ma, and T. Zhu, βStudy on soil-pile-structure interaction system by shaking table model test,βJournal of Tongji University, vol. 29, no. 7, pp. 763β768, 2001.
[4] Z. Zhang, C. Cho, X. Lu, and M. Lou, βShaking table tests ofthe damping behavior of SSI systems,β China Civil EngineeringJournal, vol. 43, no. 2, pp. 100β104, 2010 (Chinese).
[5] A. Barcena and L. Esteva, βInfluence of dynamic soil-structureinteraction on the nonlinear response and seismic reliabilityof multistorey systems,β Earthquake Engineering & StructuralDynamics, vol. 36, no. 3, pp. 327β346, 2007.
[6] H. Matinmanesh and M. S. Asheghabadi, βSeismic analysis onsoil-structure interaction of buildings over sandy soil,β ProcediaEngineering, vol. 14, pp. 1737β1743, 2011.
[7] D. Gilles, G. McClure, and L. E. Chouinard, βUncertaintyin fundamental period estimates leads to inaccurate designseismic loads,β Canadian Journal of Civil Engineering, vol. 38,no. 8, pp. 870β880, 2011.
[8] F. Mirzaie, M. Mahsuli, and M. A. Ghannad, βProbabilisticanalysis of soil-structure interaction effects on the seismic per-formance of structures,β Earthquake Engineering & StructuralDynamics, 2016.
[9] X. Gao, X.-Z. Ling, L. Tang, and P.-J. Xu, βSoil-pile-bridgestructure interaction in liquefying ground using shake tabletesting,β Soil Dynamics and Earthquake Engineering, vol. 31, no.7, pp. 1009β1017, 2011.
[10] Y. Chen, X. Lu, P. Li, and B. Chen, βComparative study onthe dynamic soil-structure interaction systemwith various soils
by using shaking table model tests,β China Civil EngineeringJournal, vol. 39, no. 5, pp. 57β64, 2006.
[11] D. Pitilakis, M. Dietz, D. M. Wood, D. Clouteau, and A.Modaressi, βNumerical simulation of dynamic soil-structureinteraction in shaking table testing,β Soil Dynamics and Earth-quake Engineering, vol. 28, no. 6, pp. 453β467, 2008.
[12] K. T. Chau, C. Y. Shen, andX. Guo, βNonlinear seismic soil-pile-structure interactions: Shaking table tests and FEM analyses,βSoil Dynamics and Earthquake Engineering, vol. 29, no. 2, pp.300β310, 2009.
[13] L. Tang, X. Ling, P. Xu, X. Gao, and D. Wang, βShake table testof soil-pile groups-bridge structure interaction in liquefiableground,β Earthquake Engineering and Engineering Vibration,vol. 9, no. 1, pp. 39β50, 2010.
[14] A. Collina, A. Facchinetti, F. Fossati, and F. Resta, βHardwarein the loop test-rig for identification and control application onhigh speed pantographs,β Shock and Vibration, vol. 11, no. 3-4,pp. 445β456, 2004.
[15] M. Nakashima, H. Kato, and E. Takaoka, βDevelopment ofrealβtime pseudo dynamic testing,β Earthquake Engineering &Structural Dynamics, vol. 21, no. 1, pp. 79β92, 1992.
[16] K. Konagai and T. Nogami, βAnalog circuit to simulate dynamicsoil-structure interaction in shake table test,β Soil Dynamics andEarthquake Engineering, vol. 17, no. 5, pp. 279β287, 1998.
[17] K. Konagai and R. Ahsan, βSimulation of nonlinear soil-structure interaction on a shaking table,β Journal of EarthquakeEngineering, vol. 6, no. 1, pp. 31β51, 2002.
[18] A. Heath, A. P. Darby, and J. Bawcombe, βSubstructure testingfor dynamic soil-structure interaction,β in Proceedings of the2nd BritishGeotechnical Association International Conference onFoundations (ICOF β08), pp. 1213β1224,Watford, UK, June 2008.
[19] Q. Wang, J.-T. Wang, F. Jin, F.-D. Chi, and C.-H. Zhang,βReal-time dynamic hybrid testing for soil-structure interactionanalysis,β Soil Dynamics and Earthquake Engineering, vol. 31, no.12, pp. 1690β1702, 2011.
[20] M. I. Wallace, D. J. Wagg, and S. A. Neild, βAn adaptivepolynomial based forward prediction algorithm for multi-actuator real-time dynamic substructuring,β Proceedings of The
12 Shock and Vibration
Royal Society of London. Series A. Mathematical, Physical andEngineering Sciences, vol. 461, no. 2064, pp. 3807β3826, 2005.
[21] J. Y. Tu, P. Y. Lin, D. P. Stoten, and G. Li, βTesting of dynamicallysubstructured, base-isolated systems using adaptive controltechniques,β Earthquake Engineering and Structural Dynamics,vol. 39, no. 6, pp. 661β681, 2010.
[22] X. Li, D.-Y. Zhang, W.-M. Yan, Y.-J. Chen, and W.-C. Xie,βShake-table test for a typical curved bridge: wave passage andlocal site Effects,β Journal of Bridge Engineering, vol. 20, no. 2,pp. 1β14, 2015.
[23] J. Lysmer and F. E. Richart, βDynamic response of footings tovertical loading,β Journal of Engineering Mechanics, vol. 95, no.4, pp. 65β91, 1969.
[24] W.-H.Wu andW.-H. Lee, βSystematic lumped-parameter mod-els for foundations based on polynomial-fraction approxima-tion,β Earthquake Engineering and Structural Dynamics, vol. 31,no. 7, pp. 1383β1412, 2002.
[25] E. Safak, βTime-domain representation of frequency-depend-ent foundation impedance functions,β Soil Dynamics and Earth-quake Engineering, vol. 26, no. 1, pp. 65β70, 2006.
[26] X. L. Du and M. Zhao, βStability and identification for rationalapproximation of frequency response function of unboundedsoil,β Earthquake Engineering and Structural Dynamics, vol. 39,no. 2, pp. 165β186, 2010.
[27] N. Makris and G. Gazetas, βDynamic pile-soil-pile interaction.Part II: lateral and seismic response,β Earthquake Engineering &Structural Dynamics, vol. 21, no. 2, pp. 145β162, 1992.
[28] N. Makris and G. Gazetas, βDisplacement phase differences ina harmonically oscillating pile,β Geotechnique, vol. 43, no. 1, pp.135β150, 1993.
[29] D. P. Stoten, C. N. Lim, and S. A. Neild, βAssessment ofcontroller strategies for real-time dynamic substructuring ofa lightly damped system,β Proceedings of the Institution ofMechanical Engineers. Part I: Journal of Systems and ControlEngineering, vol. 221, no. 2, pp. 235β250, 2007.
[30] S. A. Neild, D. P. Stoten, D. Drury, and D. J. Wagg, βControlissues relating to real-time substructuring experiments using ashaking table,β Earthquake Engineering and Structural Dynam-ics, vol. 34, no. 9, pp. 1171β1192, 2005.
[31] D. P. Stoten and R. A. Hyde, βAdaptive control of dynamicallysubstructured systems: the single-input single-output case,βProceedings of the Institution of Mechanical Engineers. Part I:Journal of Systems and Control Engineering, vol. 220, no. 2, pp.63β79, 2006.
[32] D. P. Stoten, J. Y. Tu, and G. Li, βSynthesis and control ofgeneralized dynamically substructured systems,β Proceedings ofthe Institution ofMechanical Engineers. Part I: Journal of Systemsand Control Engineering, vol. 223, no. 3, pp. 371β392, 2009.
[33] J. Guo, Z. Tang, S. Chen, and Z. Li, βControl strategy forthe substructuring testing systems to simulate soil-structureinteraction,β Smart Structures and Systems, vol. 18, no. 6, pp.1169β1188, 2016.
[34] X. L. Du and J. F. Zhao, βTime domain analysis procedureaccounting for frequency dependent of foundation impedance,βin Proceedings of the 19th International Conference on StructuralMechanics in Reactor Technology, Toronto, Canada, August2007.
International Journal of
AerospaceEngineeringHindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
RoboticsJournal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Active and Passive Electronic Components
Control Scienceand Engineering
Journal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
International Journal of
RotatingMachinery
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Journal ofEngineeringVolume 2014
Submit your manuscripts athttps://www.hindawi.com
VLSI Design
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Shock and Vibration
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Civil EngineeringAdvances in
Acoustics and VibrationAdvances in
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Electrical and Computer Engineering
Journal of
Advances inOptoElectronics
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
The Scientific World JournalHindawi Publishing Corporation http://www.hindawi.com Volume 2014
SensorsJournal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Modelling & Simulation in EngineeringHindawi Publishing Corporation http://www.hindawi.com Volume 2014
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Chemical EngineeringInternational Journal of Antennas and
Propagation
International Journal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
Navigation and Observation
International Journal of
Hindawi Publishing Corporationhttp://www.hindawi.com Volume 2014
DistributedSensor Networks
International Journal of