Experimental Verification of Methods for Converting Acceleration Data in High-Rise Buildings into Displacement Data by Shaking Table Test

When diagnosing damage to high-rise buildings during earthquakes, it is necessary to measure the displacement of each story. However, with respect to accuracy and cost, it is most reasonable to convert acceleration into displacement. In this study, shake table testing was carried out to verify the conversion methods, converting the acceleration data measured in a high-rise building into velocity and displacement. In the shaking table test, the displacement of a 10-story model building under strong motion was measured using high-speed imaging devices. High-speed images were taken at 1000 frames per second, reflecting the dynamic behavior of the model building. Then, this displacement was compared with the displacement obtained by processing the acceleration data. This study applied three methods for correcting and converting acceleration into velocity and displacement. Method 1 used the transfer function, H2(ω), which reflects the dynamic characteristics of the system. The displacements converted by this method showed the lowest accuracy, because the transfer function depends on the dynamic characteristics of the structure. Method 2 used the cosine Fourier transform for baseline correction, and the discrete input data are calculated as the sum of the cosine functions. Method 3 used the least-squares fitting in the first step to remove the linear drift in the acceleration and applied the high-pass Butterworth filter. The displacements converted by Method 2 were the most reliable, and were close to the displacements measured in the shaking table test. However, the response of high-rise buildings is affected by lowand high-frequency noise. It is necessary to further investigate the limitations and applicability of the conversion methods for providing reliable displacement of the building.


Introduction
The raw ground motion signals recorded by seismologists are always preprocessed before any engineering and seismological analysis takes place.Strong-motion data processing has two main objectives in making the data useful for engineering analysis: (1) correction of the strong motion sensor itself (analogue or digital sensors can be used), and (2) reduction of random noise in the recorded signals [1].Different authors and agencies around the world have used a variety of steps in data processing.The three major organizations in the United States for seismic analysis, the USGS (US Geological Survey), PEER (Pacific Earthquake Engineering Research Centre), and CSMIP (California Strong-Motion Instrumentation Program), also use different signal processing techniques to process records [2].Other European countries have also reported the specific data processing steps that they have adopted [3].Mollova (2006) presented the application of digital filtering in the data processing of acceleration records from earthquakes, whereby Butterworth, Chebyshev, and Bessel filters with different orders were employed to eliminate noise.One dataset under investigation included accelerograms from three stations located in Turkey, Dinar, Izmit, and Kusadasi, withal of which use an analogue seismograph, SMA-1 [2].Jones et al., (2018) presented the Processing and Review Interface for Strong Motion data (PRISM) software, developed by the National Strong Motion Project (NSMP) of the U.S. Geological Survey (USGS) [4].Corrected acceleration, velocity, and displacement time series could be obtained by applying baseline correction methods and filtering to ensure that the spectral content of the resulting products was within the sensors' frequency response range, and exhibited a signal-to-noise ratio (SNR) significantly improved with respect to the expected background noise level [5,6].In the United States and Europe, much research has been done on strong-motion data processing with the aim of accurately evaluating ground motion [7][8][9][10].Seismometers are able to measure the vertical and horizontal components of earthquakes.Analogue filters can be used for the data processing of analogue seismometers, with Butterworth [11] and Chebyshev filters being the most typical analogue filters used.The resulting frequency characteristics are thus referred to as Butterworth and Chebyshev characteristics, respectively [12].
This signal processing can be carried out for the seismological analysis of ground motion by seismologists.However, to perform a diagnosis of the structure, an accelerometer is installed in high-rise buildings [13].In order to carry out the diagnosis and health monitoring of high-rise buildings, the raw data measured by the accelerometer must be analyzed by seismologists directly, without preprocessing.This is because noise is generated due to the dynamic characteristics of the structure, as well as artificial factors such as walking and machine operation.Noise needs to be removed from the measured acceleration data, and the displacement should be calculated with baseline correction.Then, the converted displacement can be used to evaluate the baseline drift of the high-rise building.Therefore, acceleration must be converted into velocity and displacement through filtering and baseline correction.
The causes and characteristics of noise and baseline drift are as follows.Low-frequency signals, including acceleration, velocity, and displacement waveforms, are as crucial as high-frequency signals for many areas of study in the fields of seismology and earthquake engineering.For example, velocity signals are associated with the kinetic energy of seismic waves, and low-frequency seismic waves contribute substantially to the responses of high-rise buildings, as well as huge structures such as dams and bridges.Unfortunately, these low-frequency seismic waves are easily contaminated by baseline errors.Most baseline errors in seismic signals are negligible with respect to acceleration, but become larger with respect to velocity and displacement waveforms.A displacement waveform based directly on the integration of the raw accelerograms, without applying any preprocessing, may deviate greatly from the baseline.Most of these baseline errors are much larger than the ground displacement itself, and as such, they are larger than the limits for the usefulness of low-frequency signals.Therefore, baseline correction is an essential step in extracting more low-frequency information from strong-motion accelerograms [14].
Figure 1 shows the results of the numerical calculation of velocity and displacement from acceleration without baseline correction or filtering, measured by means of the shaking table test.The figure shows that the converted displacement is unreliable.In this paper, the abbreviation AVD (Acceleration-Velocity-Displacement) is used to refer to the conversion of acceleration into velocity and displacement.Velocity can be obtained by integrating acceleration, while displacement can be obtained either through the integration of velocity or a double integration of acceleration.However, this basic principle of physics is not valid for disseminated data acquired from various strong-motion data providers [15].Chiu (2012) examined several examples, which included three sets of three-component waveforms disseminated by the Berkeley Digital Seismic Network (BDSN), the National Strong-Motion Program of the US Geological Survey (USGS), and the California Strong-Motion Instrumentation Program of the California Geological Survey (CGS).The displacement waveforms were different from the corresponding velocity waveforms, which showed a significant linear baseline drift.Therefore, the initial data for velocity and displacement were assumed to be zero in order to remove the linear baseline drift [15].This incompatibility appears not only in displacement data, but also in data derived from acceleration, such as in long-period response spectra [16].Pecknold and Riddell [17,18] also found that incompatible acceleration may lead to low-frequency distortion in the calculation of structural response, and this effect becomes more serious in nonlinear cases [18].The idea underlying signal processing is the removal of noise from strong-motion data in order to be able to extract valuable information, such as peak ground characteristics and response spectra.Although noise causes distortion of the true signal over a wide frequency range, its influence is more dramatic in the low-and high-frequency ranges, where the signal-to-noise ratio tends to be significantly lower.Nowadays, in the signal processing of strong motion, frequencies lower than 1 Hz and higher than 20 Hz are regarded as low-and high-frequency limits, respectively [19].Athanasiou et al. (2018) removed low-frequency noise using a simple baseline-fitting scheme.The signal-processing scheme they developed and used defines the duration of the main event, removes background noise, and adjusts the distorted baseline drift using a polynomial curve.
Appl.Sci.2019, 1, x FOR PEER REVIEW 3 of 25 baseline drift.Therefore, the initial data for velocity and displacement were assumed to be zero in order to remove the linear baseline drift [15].This incompatibility appears not only in displacement data, but also in data derived from acceleration, such as in long-period response spectra [16].
Pecknold and Riddell [17,18] also found that incompatible acceleration may lead to low-frequency distortion in the calculation of structural response, and this effect becomes more serious in nonlinear cases [18].The idea underlying signal processing is the removal of noise from strong-motion data in order to be able to extract valuable information, such as peak ground characteristics and response spectra.Although noise causes distortion of the true signal over a wide frequency range, its influence is more dramatic in the low-and high-frequency ranges, where the signal-to-noise ratio tends to be significantly lower.Nowadays, in the signal processing of strong motion, frequencies lower than 1 Hz and higher than 20 Hz are regarded as low-and high-frequency limits, respectively [19].Athanasiou et al. ( 2018) removed low-frequency noise using a simple baseline-fitting scheme.The signal-processing scheme they developed and used defines the duration of the main event, removes background noise, and adjusts the distorted baseline drift using a polynomial curve.Several previous studies have reported the error and baseline drift in earthquake accelerometer records; these can be caused by various factors arising from the vicinity of the location at which the earthquake accelerometer is installed [14][15][16][17][18][19].The frequency band is distributed over a wide range, from 0.001 to several tens of Hz.Generally, low-frequency noise below about 0.1 Hz is called a microseism, and high-frequency noise from 0.5 to several tens of Hz is referred to a microtremor [20].Generally, noise in the range of 0.1 to 0.5 Hz is caused by waves on the sea surface, and becomes weaker as the distance from the coast gets greater.Noise lower than 0.1 Hz is known to be the main cause of changes in micro-pressure [21,22].By contrast, high-frequency noise above 1 Hz is typically caused by transportation and industrial facilities, other human factors, and natural factors, such as wind and waves [23,24].In particular, high-rise buildings are greatly influenced by natural factors such as wind, as well as artificial factors due to mechanical vibration and walking load.
Previous research on strong motion has identified the causes of noise in the frequency domain of accelerometers [14].A variety of processing and filtering techniques have been developed to compensate for seismic accelerometers [15,25].However, not much research has been done on Several previous studies have reported the error and baseline drift in earthquake accelerometer records; these can be caused by various factors arising from the vicinity of the location at which the earthquake accelerometer is installed [14][15][16][17][18][19].The frequency band is distributed over a wide range, from 0.001 to several tens of Hz.Generally, low-frequency noise below about 0.1 Hz is called a microseism, and high-frequency noise from 0.5 to several tens of Hz is referred to a microtremor [20].Generally, noise in the range of 0.1 to 0.5 Hz is caused by waves on the sea surface, and becomes weaker as the distance from the coast gets greater.Noise lower than 0.1 Hz is known to be the main cause of changes in micro-pressure [21,22].By contrast, high-frequency noise above 1 Hz is typically caused by transportation and industrial facilities, other human factors, and natural factors, such as wind and waves [23,24].In particular, high-rise buildings are greatly influenced by natural factors such as wind, as well as artificial factors due to mechanical vibration and walking load.
Previous research on strong motion has identified the causes of noise in the frequency domain of accelerometers [14].A variety of processing and filtering techniques have been developed to compensate for seismic accelerometers [15,25].However, not much research has been done on correction and filtering methods for accelerometers installed in high-rise buildings.The data obtained from accelerometers installed in a free field depend on the characteristics of the ground, and the artificial influence is relatively small.However, artificial factors have a great influence on the noise and dynamic characteristics of high-rise buildings.Therefore, it is necessary to verify the method by which acceleration data of high-rise buildings is converted into displacement data.In particular, acceleration measured over long periods leads to large error and baseline drift in displacement in the diagnostic assessment of high-rise buildings.The displacement of high-rise buildings can be measured by GPS, but high accuracy is required for soundness assessment, and the corresponding GPS is very expensive.High-speed imaging can measure the displacement of a structure, but it is difficult to perform measurements at all times, and challenging to apply this technique to high-rise buildings [26].Therefore, it is most efficient to install an accelerometer in high-rise buildings and convert this data into displacement.
In this study, the shaking table test was performed to verify the conversion methods.This involves converting the data from an accelerogram installed in a tall building into velocity and displacement data.In the shaking table test, the displacement of a 10-story model building under strong motion was measured using high-speed imaging.High-speed images were taken at 1000 frames per second, capturing the dynamic behavior of the model building.Then, this displacement was compared with the displacement obtained by processing the acceleration data.This study applied three methods for correcting and converting acceleration to velocity and displacement.The characteristics and the applicability of the conversion methods for carrying out diagnosis of high-rise buildings are evaluated in this study.

Shaking Table Test
In general, the shaking table test is performed in order (1) to identify the dynamic characteristics of the target structure system, and (2) to obtain a validated analytical model.It can also be used (3) to verify the durability of the structure, and (4) to diagnose and maintain the integrity of the structure.Changes in the vibration characteristics of a structure often indicate a problem related to the structural integrity of the structure.For example, when a structure is subjected to a load and reaches fracture, a crack occurs.In this case, the natural frequency is significantly different from the case in which there is no crack.Therefore, if there is a significant change in the natural frequency, it can be predicted that the safety of the structure is likely to be a problem [27][28][29][30][31].
In this study, the shaking table test was performed in order to verify the method used to convert acceleration data into velocity and displacement.
Figure 2 shows the shaking table test of a high-rise building.The shaking table test consists of a shaking table, a scale model consisting of 10 stories above ground and three stories below, as well as four accelerometers and a high-speed camera.The displacements of the high-rise building were estimated by analyzing images taken by a high-speed camera.The acceleration data measured by the accelerometer were then converted into displacement and compared with the displacement by image analysis.Figure 2c shows a model of a shaking table and a high-rise building.The height of the model is 147.2 cm, the height above the ground surface (upper part) is 113.0 cm, and the depth below the ground surface (lower part) is 34.2 cm.The lower part is a square of 40 cm × 40 cm, and the upper part is a square of 20 cm × 20 cm.The accelerometer was installed at the base, 4th floor, 7th floor, and roof.The accelerometer used in the shaking table test was iLOG-MEMS-ACC, which is developed and sold by Smart C & S as shown in Figure 3. iLOG-MEMS-ACC is a wireless acceleration sensor system that rapidly sends (115.2 kbps) acceleration signals without noise, using Bluetooth.The sensor device consists of a 2-axis MEMS acceleration sensor, a 16 bit A/D convertor, an 8 bit micro-controller, internal memory (2.0 MB), a Bluetooth communication module, and a Li-ion recharge module, which allows self-acquisition of the acceleration signal, as well as decision-making according to an embedded algorithm.The detailed specifications of iLOG-MEMS-ACC are shown in Table 1.The accelerometer used in the shaking table test was iLOG-MEMS-ACC, which is developed and sold by Smart C & S as shown in Figure 3. iLOG-MEMS-ACC is a wireless acceleration sensor system that rapidly sends (115.2 kbps) acceleration signals without noise, using Bluetooth.The sensor device consists of a 2-axis MEMS acceleration sensor, a 16 bit A/D convertor, an 8 bit micro-controller, internal memory (2.0 MB), a Bluetooth communication module, and a Li-ion recharge module, which allows algorithm.The detailed specifications of iLOG-MEMS-ACC are shown in Table 1.

High-Speed Camera Shooting and Image Analysis
Table 2 and Figure 4 show the equipment used for high-speed camera shooting and image analysis, including a high-speed camera, targets, and an image analysis program.To estimate the displacements of super high-rise and complex facilities at the time of excitation, the high-speed camera was used for the image capture at 1000 frames per second and image analysis; (x, y, z) is the coordinate value of the tracker attached to the structure using a TEMA (Telecommunication Engineering and Manufacturing Association) image analysis program.The tracked tracker comprises the reference coordinate system, and then relative displacement is measured.

High-Speed Camera Shooting and Image Analysis
Table 2 and Figure 4 show the equipment used for high-speed camera shooting and image analysis, including a high-speed camera, targets, and an image analysis program.To estimate the displacements of super high-rise and complex facilities at the time of excitation, the high-speed camera was used for the image capture at 1000 frames per second and image analysis; (x, y, z) is the coordinate value of the tracker attached to the structure using a TEMA (Telecommunication Engineering and Manufacturing Association) image analysis program.The tracked tracker comprises the reference coordinate system, and then relative displacement is measured.As a device capable of shooting over 1000 frames per second, it captures moments that cannot be observed by the human eye and analyzes the behavior of the structure using these records.Table 2. Equipment and features used in high-speed camera shooting and image analysis.

Equipment Name
Model Name (Manufacturer) Special Feature Amount High-speed camera

Phantom V641 (Vision research)
As a device capable of shooting over 1000 frames per second, it captures moments that cannot be observed by the human eye and analyzes the behavior of the structure using these records.The test procedure for the image analysis and the high-speed camera setting procedure is described here.First, quadrant targets were attached to the corners of each layer in order to measure the three-axis displacement of each layer using a vibrating tester.The size of the target was chosen to be sufficiently large to identify and analyze each frame for a high-speed camera shooting.To measure three-dimensional coordinates using a high-speed camera image, two high-speed cameras should be used to start shooting at the same point in time.The target to be observed (target attached to each layer) should be captured on each camera image.The position of the camera should be adjusted so that the target does not deviate out of the image while the excited structure is moving.The camera angle is adjusted to be between 45 degrees and 180 degrees.
To correct the distortion caused by the lens curvature, which is a problem for all camera lenses, it is necessary to take certain calibration values (parameter values for distortion correction) by shooting the calibration board several times.A high-speed camera with a shooting speed of over 1000 frames per second has a limited shooting time, depending on resolution and frames per second, which depend on the camera's performance.Since the high-speed camera used in this test requires 50 s of shooting per experiment, the resolution was set at 2560 × 1600 pixels.
The basic concept of image analysis is to obtain a tracking shot of the trackers (targets) in each image taken by the high-speed camera, separating each frame using the algorithm included in the image analysis program, and converting the frames into data for easy analysis before finally reporting it.The data obtained by the image analysis program include the distance, speed, acceleration, relative displacement, angle, and angular velocity of the tracker.The test procedure for the image analysis and the high-speed camera setting procedure is described here.First, quadrant targets were attached to the corners of each layer in order to measure the three-axis displacement of each layer using a vibrating tester.The size of the target was chosen to be sufficiently large to identify and analyze each frame for a high-speed camera shooting.To measure three-dimensional coordinates using a high-speed camera image, two high-speed cameras should be used to start shooting at the same point in time.The target to be observed (target attached to each layer) should be captured on each camera image.The position of the camera should be adjusted so that the target does not deviate out of the image while the excited structure is moving.The camera angle is adjusted to be between 45 degrees and 180 degrees.
To correct the distortion caused by the lens curvature, which is a problem for all camera lenses, it is necessary to take certain calibration values (parameter values for distortion correction) by shooting the calibration board several times.A high-speed camera with a shooting speed of over 1000 frames per second has a limited shooting time, depending on resolution and frames per second, which depend on the camera's performance.Since the high-speed camera used in this test requires 50 s of shooting per experiment, the resolution was set at 2560 × 1600 pixels.
The basic concept of image analysis is to obtain a tracking shot of the trackers (targets) in each image taken by the high-speed camera, separating each frame using the algorithm included in the image analysis program, and converting the frames into data for easy analysis before finally reporting it.The data obtained by the image analysis program include the distance, speed, acceleration, relative displacement, angle, and angular velocity of the tracker.

AVD Conversion Method
Method 1 is called the Penzien filter, and it uses DFT and H 2 (ω) [32].Method 2 uses the baseline correction and cosine correction filter [33].Method 3 uses a least-square fitting and 4th-order Butterworth filter [14].The Penzien filter in Method 1 is a filtering process using the Fourier transform Appl.Sci.2019, 9, 1653 8 of 25 and transfer function, and it removes low-frequency noise by means of a high-pass filter.Method 2 computes the discrete data of the input seismic waves as the sum of the cosine functions by means of a cosine Fourier transform.In Method 3, the least-squares fitting implemented in the first step removes the linear drift in the acceleration.In this study, the long-term measurement of the acceleration of the high-rise building incorporates the error and drift of the baseline.This is the cause of some of the error in the measuring equipment, as well as the background noise.
At this point, it is instructive to once again solve for the steady-state harmonic response using the exponential form of a solution.Consider the general case of harmonic loading expressed in exponential form: ..
where φ is an arbitrary phase angle in the harmonic loading function, is loading frequency, and ω is natural frequency.In dealing with completely general harmonic loads, especially for the case of periodic loading where the earthquake is expressed as a series of harmonic terms, it is essential to define the input phase angle for each harmonic; however, this is usually most conveniently accomplished by expressing the input in complex number form rather than on the basis of amplitude and phase angle.In this section, only a single harmonic loading term will be considered; therefore, its phase angle is arbitrarily taken to be zero for the sake of simplicity.Therefore, it need not be included in the loading expression [32].
The specific solutions of Equation ( 1) and its first and second-time derivatives are . ..
where G is a complex constant.To evaluate these, substitute Equations ( 2)-(4) into Equation (1), cancel out the common quantity, exp(i t), in each term, substitute k ω 2 for m, and β for ω , and solve for G, which yields Equation ( 5) is rearranged as Equation (7).
By substituting Equation (11) into Equation ( 10), the displacement in the time domain can be calculated using Equation (12).
Here, H 2 (ω) is a correction filter in the frequency domain.
Since it is difficult to set the reference point (fixed point) in order to measure the temporal change, it is necessary to measure the temporal change of the acceleration, and then to calculate the displacement by integrating it twice.However, when the acceleration, including the measurement error, is integrated twice and integrated into the displacement, the displacement includes the amplified additional error, due to the two integral transformations.
In this study, numerical integration, such as Fourier transform, has been used to convert acceleration into velocity and displacement.The Fourier transform is more simply a numerical integration; it is the decomposition of a periodic loading into a harmonic series.Here, H 2 (ω) is a correction filter in the frequency domain [32].
Figure 5 shows the flowchart for the conversion of the acceleration data in Method 1.The acceleration data measured in the time domain should be transformed by the Fourier transform function in order to convert them into the frequency domain.
In addition, A f (ω) is calculated again in order to correct the data in the frequency domain; H 2 (ω) is the correction filter multiplied to obtain the corrected Fourier transformed function A c (ω).
Then, velocity and displacement are calculated in the frequency domain.
Finally, to convert them into velocity and displacement in the time domain, the inverse Fourier transform function is used.
In this process, the content of the original signal is modulated on the basis of the transform function, reflecting the system characteristics.After being passed through the high-pass filter used at this juncture, the low-frequency components are removed.This solution, using the time domain method, solves the equation of motion in the time domain directly.In contrast, solutions that use the frequency domain method are obtained by converting the equation of motion in the time domain into the frequency domain.Generally, by using a Fourier transform to convert the equation of motion into the frequency domain, a solution can be found more efficiently than by finding it in the time domain.Moreover, by applying a discrete Fourier transform (DFT) using a fast Fourier transform (FFT) algorithm with a Fourier transform, the integrals used in the Fourier transform can be easily and quickly processed.This solution, using the time domain method, solves the equation of motion in the time domain directly.In contrast, solutions that use the frequency domain method are obtained by converting the equation of motion in the time domain into the frequency domain.Generally, by using a Fourier transform to convert the equation of motion into the frequency domain, a solution can be found more efficiently than by finding it in the time domain.Moreover, by applying a discrete Fourier transform (DFT) using a fast Fourier transform (FFT) algorithm with a Fourier transform, the integrals used in the Fourier transform can be easily and quickly processed.
The second method is time domain analysis.Figure 6 shows the flowchart for the conversion of acceleration data using Method 2. The load magnitude after the first load and the end of the measured load are assumed to be zero.There is zero padding.First, the acceleration ( ) in the time domain is prepared and integrated so as to obtain the velocity ( ) and the displacement ( ).Next, a baseline correction of the initial acceleration ( ) in the first-order linear equation is performed. and  in Figure 6 are the coefficients of the first-order linear equation for the baseline correction.Then, the corrected acceleration () is integrated in order to calculate the new velocity, () , and the displacement, () .Finally, the displacement () is cosine-corrected.Using the cosine Fourier series transform, the discrete input data is calculated as the sum of the cosine functions.The second method is time domain analysis.Figure 6 shows the flowchart for the conversion of acceleration data using Method 2. The load magnitude after the first load and the end of the measured load are assumed to be zero.There is zero padding.First, the acceleration (a 0 ) in the time domain is prepared and integrated so as to obtain the velocity (v 0 ) and the displacement (d 0 ).Next, a baseline correction of the initial acceleration (a 0 ) in the first-order linear equation is performed.c 0 and c 1 in Figure 6 are the coefficients of the first-order linear equation for the baseline correction.
Then, the corrected acceleration (a(t)) is integrated in order to calculate the new velocity, (v(t)), and the displacement, (d(t)).Finally, the displacement (d(t)) is cosine-corrected.Using the cosine Fourier series transform, the discrete input data is calculated as the sum of the cosine functions.This solution, using the time domain method, solves the equation of motion in the time domain directly.In contrast, solutions that use the frequency domain method are obtained by converting the equation of motion in the time domain into the frequency domain.Generally, by using a Fourier transform to convert the equation of motion into the frequency domain, a solution can be found more efficiently than by finding it in the time domain.Moreover, by applying a discrete Fourier transform (DFT) using a fast Fourier transform (FFT) algorithm with a Fourier transform, the integrals used in the Fourier transform can be easily and quickly processed.
The second method is time domain analysis.Figure 6 shows the flowchart for the conversion of acceleration data using Method 2. The load magnitude after the first load and the end of the measured load are assumed to be zero.There is zero padding.First, the acceleration ( ) in the time domain is prepared and integrated so as to obtain the velocity ( ) and the displacement ( ).Next, a baseline correction of the initial acceleration ( ) in the first-order linear equation is performed. and  in Figure 6 are the coefficients of the first-order linear equation for the baseline correction.Then, the corrected acceleration () is integrated in order to calculate the new velocity, () , and the displacement, () .Finally, the displacement () is cosine-corrected.Using the cosine Fourier series transform, the discrete input data is calculated as the sum of the cosine functions.The third method uses the least-square fitting and Butterworth filters for digital strong-motion data, as proposed by Chiu [14].Figure 7 shows the flowchart for the conversion of acceleration data using Method 3, wehreby a 0 (t) is the initial acceleration, a(t) is the corrected acceleration, ω is natural frequency, ω c is cutoff frequency, and G 0 is a frequency response at zero frequency.The least-squares fitting in the first step removes the linear trend of the acceleration.The long-term measurement of the acceleration of a high-rise building incorporates the errors and drift of the baseline.These are caused by errors in the measuring equipment, along with background error.The correction method for their removal consists of the following three steps.First, the baseline drift of the acceleration is removed using the least-squares method, and the acceleration in the time domain is converted into the frequency domain by Fourier transform.Next, low-frequency error is removed using the fourth-order Butterworth filter.The Butterworth filter does not generate a ripple in the passband; it has better performance than the Bessel filter, and the bandpass filter can simultaneously remove low-frequency and high-frequency noise.Because of these advantages, passband Butterworth filters are used in the PEER (Pacific Earthquake Engineering Research Centre) strong-motion database [34].The Butterworth filter is also used in the Strong Motion data (PRISM) software developed by the US Geological Survey (USGS) [4].The accelerometer installed in the high-rise building is an analogue system, and a typical analogue Butterworth filter is used.The last step required is to remove the linear trend that arises from having a nonzero initial velocity.When the strong-motion instrument records seismic data, those data may come from the effect of the background noise.To ensure that the strong motion records start from zero, the records should be processed by means of baseline initialization [35].Furthermore, motions with very long periods must be sacrificed, and a combination of baseline correction and low-cut filtering can be used to produce motions that are accurate representations of the true ground motions obtained by many earthquake modelers [36].The third method uses the least-square fitting and Butterworth filters for digital strong-motion data, as proposed by Chiu [14].Figure 7 shows the flowchart for the conversion of acceleration data using Method 3, wehreby  () is the initial acceleration, () is the corrected acceleration, ω is natural frequency,  is cutoff frequency, and  is a frequency response at zero frequency.The least-squares fitting in the first step removes the linear trend of the acceleration.The long-term measurement of the acceleration of a high-rise building incorporates the errors and drift of the baseline.These are caused by errors in the measuring equipment, along with background error.The correction method for their removal consists of the following three steps.First, the baseline drift of the acceleration is removed using the least-squares method, and the acceleration in the time domain is converted into the frequency domain by Fourier transform.Next, low-frequency error is removed using the fourth-order Butterworth filter.The Butterworth filter does not generate a ripple in the passband; it has better performance than the Bessel filter, and the bandpass filter can simultaneously remove low-frequency and high-frequency noise.Because of these advantages, passband Butterworth filters are used in the PEER (Pacific Earthquake Engineering Research Centre) strongmotion database [34].The Butterworth filter is also used in the Strong Motion data (PRISM) software developed by the US Geological Survey (USGS) [4].The accelerometer installed in the high-rise building is an analogue system, and a typical analogue Butterworth filter is used.The last step required is to remove the linear trend that arises from having a nonzero initial velocity.When the strong-motion instrument records seismic data, those data may come from the effect of the background noise.To ensure that the strong motion records start from zero, the records should be processed by means of baseline initialization [35].Furthermore, motions with very long periods must be sacrificed, and a combination of baseline correction and low-cut filtering can be used to produce motions that are accurate representations of the true ground motions obtained by many earthquake modelers [36].

Test Conditions
Figure 8 shows the excitation force of the shaking table test.In the test, the accelerations from three different earthquake-the Manjil earthquake, the Superstition Hills earthquake, and the Kobe earthquake-were used.These three sets of acceleration data were selected from the 22 earthquakes presented in FEMA P-695, taking place between 1974 and 1999 [37].The Manjil earthquake occurred in 1990 in the city of Manjil in northern Iran.The shock had a moment magnitude of 7.4.The earthquake in Superstition Hills in California occurred in 1987, and had a magnitude of 6.6.The Kobe earthquake occurred in 1995.It measured 6.9 on the moment magnitude scale.These three earthquakes had in common that they occurred in the lower part of the metropolitan area.Because

Test Conditions
Figure 8 shows the excitation force of the shaking table test.In the test, the accelerations from three different earthquake-the Manjil earthquake, the Superstition Hills earthquake, and the Kobe earthquake-were used.These three sets of acceleration data were selected from the 22 earthquakes presented in FEMA P-695, taking place between 1974 and 1999 [37].The Manjil earthquake occurred in 1990 in the city of Manjil in northern Iran.The shock had a moment magnitude of 7.4.The earthquake in Superstition Hills in California occurred in 1987, and had a magnitude of 6.6.The Kobe earthquake occurred in 1995.It measured 6.9 on the moment magnitude scale.These three earthquakes had in common that they occurred in the lower part of the metropolitan area.Because skyscrapers are located in big cities, the acceleration data were selected from earthquakes that caused a great damage to large cities.
The accelerations and displacements of the structures generated during the excitation were measured by both accelerometer and high-speed imaging.The peak of the Manjil earthquake in Figure 8a is 0.4272 g, and the duration is 53.50 s.The peak of the Superstition Hills earthquake in Figure 8b is 0.45755 g, and the duration is 22.29 s.The peak of the Kobe earthquake in Figure 8c is 0.48557 g, and the duration is 40.95 s.
Appl.Sci.2019, 1, x FOR PEER REVIEW 12 of 25 skyscrapers are located in big cities, the acceleration data were selected from earthquakes that caused a great damage to large cities.The accelerations and displacements of the structures generated during the excitation were measured by both accelerometer and high-speed imaging.The peak of the Manjil earthquake in Figure 8a is 0.4272 g, and the duration is 53.50 s.The peak of the Superstition Hills earthquake in Figure 8b is 0.45755 g, and the duration is 22.29 s.The peak of the Kobe earthquake in Figure 8c is 0.48557 g, and the duration is 40.95 s.

Experimental Response with Application of Strong Motion
Figure 9 shows the structural response and roof displacement for the Manjil earthquake, as measured during the shaking table test.The highest acceleration occurred at the top layer, and the displacement of the top layer, determined by high-speed imaging and analysis, was up to 141.4 mm.The resulting acceleration of the measured structure was entered into the code, converted to displacement, and compared with the measured roof displacement.

Experimental Response with Application of Strong Motion
Figure 9 shows the structural response and roof displacement for the Manjil earthquake, as measured during the shaking table test.The highest acceleration occurred at the top layer, and the displacement of the top layer, determined by high-speed imaging and analysis, was up to 141.4 mm.The resulting acceleration of the measured structure was entered into the code, converted to displacement, and compared with the measured roof displacement.Figure 10 shows the response of the Superstition Hills earthquake and the roof displacement measured during the shaking table test.The highest acceleration was generated at the top layer, like in the Manjil earthquake, and the maximum displacement was 142.4 mm, as measured through highspeed imaging and analysis.The peak excitation of the Manjil earthquake was 0.4272 g, and the peak excitation of the Superstition Hills earthquake was 0.45755 g.The maximum displacement of the Superstition Hills earthquake was almost equal to the displacement caused by the Superstition Hills earthquake.
(a) (b) Figure 10 shows the response of the Superstition Hills earthquake and the roof displacement measured during the shaking table test.The highest acceleration was generated at the top layer, like in the Manjil earthquake, and the maximum displacement was 142.4 mm, as measured through high-speed imaging and analysis.The peak excitation of the Manjil earthquake was 0.4272 g, and the peak excitation of the Superstition Hills earthquake was 0.45755 g.The maximum displacement of the Superstition Hills earthquake was almost equal to the displacement caused by the Superstition Hills earthquake. Figure 10 shows the response of the Superstition Hills earthquake and the roof displacement measured during the shaking table test.The highest acceleration was generated at the top layer, like in the Manjil earthquake, and the maximum displacement was 142.4 mm, as measured through highspeed imaging and analysis.The peak excitation of the Manjil earthquake was 0.4272 g, and the peak excitation of the Superstition Hills earthquake was 0.45755 g.The maximum displacement of the Superstition Hills earthquake was almost equal to the displacement caused by the Superstition Hills earthquake.
(a) (b)  Figure 11 shows the structural response and top layer displacement for the Kobe earthquake, as measured during the shaking table test.The highest acceleration occurred at the roof, like in the previously examined Manjil and Superstition Hills earthquakes, and the displacement of the uppermost layer as determined by high-speed imaging and analysis was up to 99.4 mm. Figure 11 shows the structural response and top layer displacement for the Kobe earthquake, as measured during the shaking table test.The highest acceleration occurred at the roof, like in the previously examined Manjil and Superstition Hills earthquakes, and the displacement of the uppermost layer as determined by high-speed imaging and analysis was up to 99.4 mm. Figure 11 shows the structural response and top layer displacement for the Kobe earthquake, as measured during the shaking table test.The highest acceleration occurred at the roof, like in the previously examined Manjil and Superstition Hills earthquakes, and the displacement of the uppermost layer as determined by high-speed imaging and analysis was up to 99.4 mm.

Conversion Process of Each Method
Three different conversion methods were used for this research, as already described in "2.3.AVD Conversion Method".The techniques used in each method are shown in Table 3.

Method
Filter Technique The acceleration response measured three times on the basis of the applied force is converted into velocity and displacement using the Penzien filter in Figure 12

Conversion Process of Each Method
Three different conversion methods were used for this research, as already described in "2.3.AVD Conversion Method".The techniques used in each method are shown in Table 3.
Table 3. Techniques for AVD conversion methods.

Method
Filter Technique The acceleration response measured three times on the basis of the applied force is converted into velocity and displacement using the Penzien filter in Figures 12-14.A description of the conversion steps and the graphs of the results are provided below.
First, the applied load is discretized.The duration of the acceleration (t d ) input for the shaking table test was 53.50 s for the Manjil earthquake, 22.29 s for the Superstition Hills earthquake, and 40.95 s for the Kobe earthquake.The measured time interval (∆t) was 0.05 s.The Manjil earthquake had 10784, the Superstition Hills earthquake had 4764, and the Kobe earthquake had 9784 data points.Therefore, FFT was applied, and the analysis time (T 0 ) was selected in consideration of the period component of the response and the load, with it duration being longer than that of the excited load.In the case of the Manjil earthquake,∆t was 0.005 s, and the number of data points (N) was 10784.The analysis time (T 0 ) was then calculated using Equation (20).Here, the size of the first load and the magnitude of the load after the end of the load measurement were assumed to be zero, and zero padding was performed.The acceleration response measured three times on the basis of the applied force is converted into velocity and displacement using the Penzien filter in Figure 12       First, the applied load is discretized.The duration of the acceleration ( ) input for the shaking table test was 53.50 s for the Manjil earthquake, 22.29 s for the Superstition Hills earthquake, and 40.95 s for the Kobe earthquake.The measured time interval (∆) was 0.05 s.The Manjil earthquake had 10784, the Superstition Hills earthquake had 4764, and the Kobe earthquake had 9784 data points.Therefore, FFT was applied, and the analysis time ( ) was selected in consideration of the period component of the response and the load, with it duration being longer than that of the excited load.Next, the discrete Fourier transform of the applied load was performed.When the applied load in the time domain is converted into the frequency domain using the discrete Fourier transform (DFT), the discrete Fourier transform F( f ) is calculated.The F( f ) of the applied load of the Manjil earthquake is shown in Figure 12a, that of the Superstition Hills earthquake is shown in Figure 13a, and that of the Kobe earthquake is shown in Figure 14a.In the discrete Fourier transform F( f ) of the excited load of the Manjil earthquake, the sampling frequency ( f s ), the maximum frequency ( f max ), and the frequency interval (∆ f ) were obtained from Equations ( 21) to (23).
Next, the complex frequency response function (or compliance) H( f ) was calculated using Equation (12).The complex frequency response function H( f ) of the Manjil earthquake is shown in Figure 12b, that of the Superstition Hills earthquake is shown in Figure 13b, and that of the Kobe earthquake is shown in Figure 14b.The responses of the displacement (U d ), velocity (U v ), and acceleration (U a ) in the frequency domain were obtained using Equations ( 24) to (26) by multiplying the complex frequency response function and the load by the frequency in the frequency domain.
The displacement (U d ), velocity (U v ), and acceleration (U a ) responses in the frequency domain were obtained using the discrete Fourier inverse transform of the response (u), velocity .u , and acceleration .. u in the time domain.The displacement (u) in the time domain was obtained using the discrete Fourier inverse transformation of Equation ( 27), the displacement (U d ) in the frequency domain was calculated using Equation (24).The displacement (u) in the time domain of the Manjil earthquake is shown in Figure 12e, that of the Superstition Hills earthquake in Figure 13e, and that of the Kobe earthquake in Figure 14e.
The velocity .
u in the time domain is determined by the discrete Fourier inverse transformation, as in Equation ( 28), and the velocity (U v ) by the frequency domain estimated using Equation (25).The displacement (u) in the time domain of the Manjil earthquake is shown in Figure 12d, that of the Superstition Hills earthquake in Figure 13d, and that of the Kobe earthquake in Figure 14d. .
u in the time domain was obtained using the discrete Fourier inverse transformation in Equation ( 29) with respect to acceleration (U a ), and in the frequency domain this was calculated as per Equation (26).The acceleration .. u in the time domain of the Manjil earthquake is shown in Figure 12c, that of the Superstition Hills earthquake in Figure 13c, and that of the Kobe earthquake in Figure 14c Here, as mentioned in the discrete Fourier transform of the load, before the discrete Fourier inverse transformation, the response in the frequency domain needs to be made a symmetrical complex conjugate with respect to f max on the right and left sides [33].
The damping ratios of high-rise buildings are related to their shape, dimensions, room size, etc., and these were determined to be 3.5~21.7%,based on room size [38].The average value, 10%, was selected for the Method 1; however, the results of Method 1 were not much affected by the damping ratios.In addition, the baseline drifts were not corrected in Method 1; therefore, the initial value of the displacement graphs did not start at zero point, meaning that zero padding is required.

Method 2
Figures 15-17 show the results of converting the accelerations of the Manjil, Superstition Hills, and Kobe earthquakes using Method 2. Method 2 computes the discrete data of the input seismic waves as the sum of the cosine functions using the cosine Fourier transform.The load magnitude after the first load and after the end of the measured load are assumed to be zero.There is zero padding.First, the acceleration (a 0 ) in the time domain is prepared and integrated so as to obtain the velocity (v 0 ) and the displacement (d 0 ).Baseline correction of the initial acceleration (a 0 ) was performed using the 1st-order equation.Then, the corrected acceleration (a) was integrated to calculate the new velocity (v) and the displacement (d).Finally, the displacement (d) was cosine corrected.Using the cosine Fourier series transform, the discrete input data were calculated as the sum of the cosine functions.
as to obtain the velocity ( ) and the displacement ( ).Baseline correction of the initial acceleration ( ) was performed using the 1st-order equation.Then, the corrected acceleration () was integrated to calculate the new velocity () and the displacement ().Finally, the displacement () was cosine corrected.Using the cosine Fourier series transform, the discrete input data were calculated as the sum of the cosine functions.

Method 3
Figure 18 Figure 19 Figure 20 show the results of converting the accelerations of the Manjil, Superstition Hills, and Kobe earthquakes using Method 3. The Butterworth filter is a type of signal processing filter designed to have a frequency response that is as flat as possible in the passband [11].It is also used to refer to a maximally flat magnitude filter.It was first described in 1930 by the British engineer and physicist, Stephen Butterworth, in his paper entitled "On the Theory of Filter Amplifiers".The filter is a 4th-order Butterworth filter design and a high-pass filter that passes through the 0.1  section.In addition,  is 1 2 of the sample frequency, and is also referred to as the cutoff frequency or Nyquist frequecy.

Method 3
Figures 18-20 show the results of converting the accelerations of the Manjil, Superstition Hills, and Kobe earthquakes using Method 3. The Butterworth filter is a type of signal processing filter designed to have a frequency response that is as flat as possible in the passband [11].It is also used to refer to a maximally flat magnitude filter.It was first described in 1930 by the British engineer and physicist, Stephen Butterworth, in his paper entitled "On the Theory of Filter Amplifiers".The filter is a 4th-order Butterworth filter design and a high-pass filter that passes through the 0.1 f c section.In addition, f c is 1  2 of the sample frequency, and is also referred to as the cutoff frequency or Nyquist frequecy.
The magnitude of the frequency response G(ω) of an n-order Butterworth low-pass filter is given in terms of the transfer function H(iω) as where n is the order of the filter, which is 4 in this study; ω c is the cutoff frequency, which in this study is 0.03; and G 0 is the frequency response at a frequency of zero.

Comparison of Acceleration Data Converted to Velocity and Displacement
Figure 21 compares the displacements at the roof of the tested model, estimated on the basis of high-speed imaging and analysis, with the displacements obtained based on conversion from acceleration data.Method 2 (baseline and cosine correction of the acceleration data) and Method 3 (least-square fitting and Butterworth filter) provided analyzed displacements that were closer to those measured than Method 1 (Penzien filter) did.Method 1 showed a substantial difference from the actual displacement.The differences between the analyzed and measured results occurred because Method 1 analyzed the raw data without zero padding, baseline correction or the least-squares method.The errors of displacement would have accumulated as a result of integrating the acceleration data twice, as shown in Figure 21.Also, H 2 (ω) is a value that reflects the dynamic characteristics of the structure, and can be calculated differently depending on the mass and stiffness of the structure.When it is difficult to accurately calculate the stiffness, the converted displacement differs from the actual value.The output of Method 3 depends on the degree of the Butterworth filter and the cutoff frequency.As shown by the National Disaster Management Institute (2012) [15], it is necessary to determine the cutoff frequency based on the number of samples per second.where n is the order of the filter, which is 4 in this study;  is the cutoff frequency, which in this study is 0.03; and  is the frequency response at a frequency of zero.

Comparison of Acceleration Data Converted to Velocity and Displacement
Figure 21 compares the displacements at the roof of the tested model, estimated on the basis of high-speed imaging and analysis, with the displacements obtained based on conversion from acceleration data.Method 2 (baseline and cosine correction of the acceleration data) and Method 3 (least-square fitting and Butterworth filter) provided analyzed displacements that were closer to those measured than Method 1 (Penzien filter) did.Method 1 showed a substantial difference from the actual displacement.The differences between the analyzed and measured results occurred because Method 1 analyzed the raw data without zero padding, baseline correction or the least-

Conclusions
In this study, the shake table test was performed to verify the conversion methods that were employed to convert acceleration data measured in a high-rise building into velocity and displacement data.Then, the displacement obtained by high-speed imaging during the shaking table test was compared with the displacement obtained by processing acceleration data.This study applied three methods for the correction and conversion of acceleration data into velocity and displacement data.The resulting characteristics and applicability of the conversion methods for evaluating the diagnosis of high-rise buildings were derived as follows.
Method 1 used the transfer function  (), which reflects the dynamic characteristics of the

Conclusions
In this study, the shake table test was performed to verify the conversion methods that were employed to convert acceleration data measured in a high-rise building into velocity and displacement data.Then, the displacement obtained by high-speed imaging during the shaking table test was compared with the displacement obtained by processing acceleration data.This study applied three methods for the correction and conversion of acceleration data into velocity and displacement data.The resulting characteristics and applicability of the conversion methods for evaluating the diagnosis of high-rise buildings were derived as follows.
Method 1 used the transfer function H 2 (ω), which reflects the dynamic characteristics of the system.Displacements converted using this method showed the lowest accuracy, because the transfer function depends on the dynamic characteristics of the structure.Method 1 did not use the zero padding and baseline correction methods.Therefore, even if discrete Fourier transformation was employed, there were a lot of differences between the analyzed and measured displacement values occurring based on a non-fixed starting point and the integration of acceleration.
Method 2 used the cosine Fourier series transform and baseline correction.Discrete input data was calculated as the sum of the cosine functions.Method 2 used the zero padding and baseline correction techniques in order to correct the initial value of the displacement.In this way, the accumulated displacement errors were corrected.Furthermore, even though the cosine correction filter functions in the time domain, the displacements were very close to the experimental values.
Method 3 used the least-squares fitting in the first step in order to remove the linear drift in the acceleration, and applied a high-pass Butterworth filter.
Method 3 did not used a zero padding; however, least-square fitting was used as the first filter.After that, the corrected acceleration data was passed through a Butterworth filter.The first filtered data was used to make the second filter efficient.Therefore, the analyzed data obtained using Method 3 are also very close to the experimental data.
For this study, there were a number of reasons for which a scale model was implemented: (1) The researchers had difficulty obtaining real data for a high-rise building; and (2) if the scale building model was to be used, the conditions and data related to strong-motion would be simplified.However, if the AVD methods for this study arweree applied to a real high-rise building, the accuracy would decrease.
Furthermore, the techniques adopted for these three methods could be separated and combined with each other.Based on Methods 1~3, the application and importance of each method is known.Based on this research, the analytical method for converting acceleration into displacement can be summarized by the following three steps: 1.
First step: A technique is required to establish a fixed point, so that the integration of acceleration data will not be subject to accumulated errors.2.
Second step: Methods for eliminating baseline drift are also taken into consideration, including baseline correction filters and least-square fitting.These could be referred to as a first filter.

3.
Third step: After these 2 steps, the noise contained in the real data needs to be eliminated.This step is the most important step.There are many potential second filters, including discrete Fourier transformation (Method 1), the cosine correction filter (Method 2), and the Butterworth filter (Method 3).
The displacements converted by Method 2 proved to be the most reliable when compared with displacements measured in the shaking table test.In addition, it is necessary to further investigate the limitations and applicability of the conversion methods for the purpose of estimating the reliable displacement of high-rise buildings.

Figure 1 .
Figure 1.Results of the numerical integration of velocity and displacement without baseline correction from the acceleration measured in the shaking table test.

Figure 1 .
Figure 1.Results of the numerical integration of velocity and displacement without baseline correction from the acceleration measured in the shaking table test; (a) Acceleration; (b) Velocity; (c) Displacement.

Figure 2 .
Figure 2. Shaking table test for verification of displacement conversion method of acceleration data; (a) Photo of experimental equipment; (b) 3D drawing; (c) 2D drawing.

Figure 2 .
Figure 2. Shaking table test for verification of displacement conversion method of acceleration data; (a) Photo of experimental equipment; (b) 3D drawing; (c) 2D drawing.
can analyze length, speed, acceleration, angle, angular velocity, etc., through a captured image or images 1 EA can analyze length, speed, acceleration, angle, angular velocity, etc., through a captured image or images 1 EAOne ea

25 Figure 5 .
Figure 5. Flowchart for the conversion of acceleration data using Method 1.

Figure 6 .
Figure 6.Flowchart for the conversion of acceleration data using Method 2.

Figure 5 .
Figure 5. Flowchart for the conversion of acceleration data using Method 1.

Figure 5 .
Figure 5. Flowchart for the conversion of acceleration data using Method 1.

Figure 6 .
Figure 6.Flowchart for the conversion of acceleration data using Method 2.

Figure 6 .
Figure 6.Flowchart for the conversion of acceleration data using Method 2.

Figure 7 .
Figure 7. Flowchart for conversion of acceleration data in Method 3.

Figure 7 .
Figure 7. Flowchart for conversion of acceleration data in Method 3.

Figure 9 .
Figure 9. Dynamic response and roof displacement of the model for the Manjil earthquake; (a) base acceleration, (b) 4th floor acceleration, (c) 7th floor acceleration, (d) roof acceleration, and (e) roof displacement.

Figure 9 .
Figure 9. Dynamic response and roof displacement of the model for the Manjil earthquake; (a) base acceleration, (b) 4th floor acceleration, (c) 7th floor acceleration, (d) roof acceleration, and (e) roof displacement.

Figure 10 .
Figure 10.Dynamic response and roof displacement of the model for the Superstition Hills earthquake; (a) base acceleration, (b) 4th floor acceleration, (c) 7th floor acceleration, (d) roof acceleration, and (e) roof displacement.

Figure 10 .
Figure 10.Dynamic response and roof displacement of the model for the Superstition Hills earthquake; (a) base acceleration, (b) 4th floor acceleration, (c) 7th floor acceleration, (d) roof acceleration, and (e) roof displacement.

Figure 10 .
Figure 10.Dynamic response and roof displacement of the model for the Superstition Hills earthquake; (a) base acceleration, (b) 4th floor acceleration, (c) 7th floor acceleration, (d) roof acceleration, and (e) roof displacement.

Figure 11 .
Figure 11.Dynamic response and roof displacement of the model for the Kobe earthquake; (a) base acceleration, (b) 4th floor acceleration, (c) 7th floor acceleration, (d) roof acceleration, and (e) roof displacement.

Figure 13 Figure 14
. A description of the conversion steps and the graphs of the results are provided below.

Figure 11 .
Figure 11.Dynamic response and roof displacement of the model for the Kobe earthquake; (a) base acceleration, (b) 4th floor acceleration, (c) 7th floor acceleration, (d) roof acceleration, and (e) roof displacement.

Figure 13 Figure 14
. A description of the conversion steps and the graphs of the results are provided below.

Figure 12 . 25 Figure 12 .
Figure 12. Analysis process for the 10th-story acceleration in the Manjil earthquake using Method 1: (a) discrete Fourier Transform of the applied load; (b) transfer function; (c) acceleration; (d) velocity; and (e) displacement.

Figure 13 .
Figure 13.Analysis process for the 10th-story acceleration in the Superstition Hills earthquake using Method 1: (a) discrete Fourier Transform of the applied load; (b) transfer function; (c) acceleration; (d) velocity; and (e) displacement.

Figure 13 .
Figure 13.Analysis process for the 10th-story acceleration in the Superstition Hills earthquake using Method 1: (a) discrete Fourier Transform of the applied load; (b) transfer function; (c) acceleration; (d) velocity; and (e) displacement.

Figure 13 .
Figure 13.Analysis process for the 10th-story acceleration in the Superstition Hills earthquake using Method 1: (a) discrete Fourier Transform of the applied load; (b) transfer function; (c) acceleration; (d) velocity; and (e) displacement.

Figure 14 .
Figure 14.Analysis process for the 10th-story acceleration in the Kobe earthquake using Method 1: (a) discrete Fourier Transform of applied load; (b) transfer function; (c) acceleration; (d) velocity; and (e) displacement.

Figure 14 .
Figure 14.Analysis process for the 10th-story acceleration in the Kobe earthquake using Method 1: (a) discrete Fourier Transform of applied load; (b) transfer function; (c) acceleration; (d) velocity; and (e) displacement. .
Appl.Sci.2019, 1, x FOR PEER REVIEW 20 of 25The magnitude of the frequency response G() of an n-order Butterworth low-pass filter is given in terms of the transfer function () as  () = |()
The shaking table test was conducted by KICT (Korea Institute of Civil Engineering and Building Technology) in May 2018 at Smart C & S (Control & Sensing) in Daejeon, Korea.The model of the high-rise building was made without scale because there was no actual reference structure.The input seismic waves were scaled considering the natural frequencies of the model, 1st = 2.8230 Hz, 2nd = 8.7980 Hz, and 3rd = 14.6990Hz.

Table 1 .
Detailed specification of acceleration (iLOG-MEMS-ACC) used in shaking table test.

Table 1 .
Detailed specification of acceleration (iLOG-MEMS-ACC) used in shaking table test.

Table 2 .
Equipment and features used in high-speed camera shooting and image analysis.

Table 3 .
Techniques for AVD conversion methods.