Wind Load and Structural Parameters Estimation from Incomplete Measurements

.e extended minimum variance unbiased estimation approach can be used for joint state/parameter/input estimation based on the measured structural responses. However, it is necessary to measure the structural displacement and acceleration responses at each story for the simultaneous identification of structural parameters and unknown wind load. A novel method of identifying structural state, parameters, and unknown wind load from incomplete measurements is proposed..e estimation is performed in a modal extended minimum variance unbiased manner, based on incomplete measurements of wind-induced structural displacement and acceleration responses. .e feasibility and accuracy of the proposed method are numerically validated by identifying the wind load and structural parameters on a ten-story shear building structure with incomplete measurements. .e effects of crucial factors, including sampling duration and the number of measurements, are discussed. Furthermore, the practical application of the developed inverse method is evaluated based on wind tunnel testing results of a 234m tall building structure. .e results indicate that the structural state, parameters, and unknown wind load can be identified accurately using the proposed approach.


Introduction
Wind load is one of the main loads during the design stage of tall buildings [1,2].Studies show that the wind load on a tall building varies depending on the terrain condition, the shape of the building, and the surrounding buildings [3].In most design processes, wind load is calculated according to the design code [4,5].However, it is difficult to calculate the time histories of wind load based on design code as the wind load design code is determined based on statistical information.Wind tunnel tests and computational fluid dynamic simulations are currently used to determine the timevarying wind load on a given structure [6,7].However, neither of the methods can exactly reproduce the incident turbulence and the characteristic of the surrounding buildings.Field measurement is regarded as the most accurate way of obtaining the time-varying wind load on tall buildings.However, due to the limitation of the wind load measurement technique, real-time measurement of wind load is difficult by field measurement.Compared to wind load measurement, the measurement of wind-induced displacement and acceleration responses is easier to achieve and more accurate.
In recent years, a lot of force estimation methods have been developed [8][9][10].Law et al. identified wind load on a 50 m guyed mast from structural responses by regularization of the identification equation [11].Lu and Law proposed a method for identifying unknown load using sensitivity measures of the dynamic response with respect to the input load [12].Liu and Shepard developed a dynamic force identification approach based on the enhanced least squares approach in the frequency domain [13].Ma et al. presented a Kalman filter method with a recursive estimator to determine the unknown excitation [14,15].e aforementioned approaches require that the unknown load has been generally point load acting at a specific DOF.However, wind loads acting on a building structure vary with space and time, and the methods in previous works cannot directly be used for wind load identification of building structures.
To address the aforementioned issue, Hwang et al. proposed the Kalman filtering approach in modal domain to estimate modal loads on a structure using limited measured response [16,17].Zhi et al. developed the Kalman filtering method and proper orthogonal decomposition technique to estimate the modal wind load on tall buildings [18,19].In 2007, Gillijns and Moor proposed an approach for joint input-state estimation in discrete-time dynamic system based on minimum variance unbiased solution [20].Lourens et al. applied the method to estimate structural responses and unknown inputs in both numerical and experimental studies [21].
is method requires no assumption or prior knowledge about the unknown inputs, and it can be used for wind load identification in the physical domain [22].Unfortunately, the aforementioned wind load estimation approaches assumed that the structural parameters are known a priori.However, for actual building structures, the structural parameters are generally determined based on the finite element model and it is difficult to calculate structural parameters exactly due to the material aging and damage in the structures.
To solve this problem, Wan et al. proposed a method called EGDF which is an extension of the unbiased minimum variance estimation for coupled state/input/parameter identification for nonlinear systems in state space [23].Song developed the joint input-state estimation technique for joint input-state-parameter estimation based on the unscented minimum variance unbiased (UMVU) estimation method [24].However, both EGDF and UMVU methods require that it is necessary to measure the acceleration responses at the locations where unknown inputs are applied, i.e., complete acceleration measurements at all DOFs are requested for simultaneous estimation of wind load and unknown structural parameters.
A novel method of the modal extended unbiased minimum variance estimation for joint state/parameter/wind load estimation from incomplete measurements is proposed in this study.e proposed method is able to simultaneously estimate the wind loads and structural parameters without using complete acceleration measurements.Moreover, the data fusion of acceleration responses and interstory displacement responses is used to prevent the drifts of the identified displacements responses and wind loads.e content of this paper is organized as follows.In Section 2, the proposed method of the modal extended unbiased minimum variance estimator is derived and the necessary mathematical proofs are given.In Section 3, a numerical validation by identifying wind loads and structural parameters from incomplete measurements on a ten-story shear building structure is addressed.Moreover, the effects of key factors, including sampling duration and the number of measurements, are discussed.In Section 4, a synchronous multipressure scanning system wind tunnel test on a 234 m tall building structure is carried out to validate the proposed method.In order to simplify the calculation, an equivalent model is carried out for experiment validation.Finally, in Section 5, the discussion and conclusion are summarized.

Modal Extended Minimum Variance Unbiased Estimation
2.1.System Model.e equation of motion of an n-degreesof-freedom building structure can be written as follows: where M, C, and K are the n × n structural mass, damping, and stiffness matrices, respectively; F(t) is the n × 1 wind load vector; € x(t), _ x(t), and x(t) are the n × 1 wind-induced acceleration, velocity, and displacement vector, respectively.
Based on modal coordinate transformation theory [25], the structural displacement vector can be obtained as follows: where Φ is the n × n mass normalized modal shape matrix and Y(t) is the n × 1 modal displacement vector.By premultiplying Φ T and using equation (2), equation ( 1) can be written as follows: where is the modal wind load vector.For proportionally damped system, the following can be obtained: where ω i is the i-th undamping natural frequency and ξ i is the i-th modal damping ratio.Equation (3) can be rewritten as follows: In general, due to the limitation of the number and location of the sensors, a reduced order representation of equation ( 6) is given: where € Y m (t), _ Y m (t), and Y m (t) are the first m order modal acceleration, modal velocity, and modal displacement vector, respectively.
). e augmented state vector consists of the modal displacement, modal velocity, and the unknown structural parameters: where θ p � [k 1 , k 2 , . . ., k n , α, β] T denotes the p × 1 unknown structural parameters vector.k i is the i-th structural stiffness coefficient.α and β are Rayleigh damping coefficients.Assuming that the unknown structural parameters are time invariant, the first-order differential equation of equation ( 8) can be obtained as follows: 2 Shock and Vibration where g(Z(t), d m (t)) is a nonlinear function consisting of the state vector, modal wind load vector, and time t.Denote that t k � k × Δt with Δt as the sampling interval and define  Z k|k and  d k as the estimated values of Z k and d k at time t k , respectively.Considering the process noise, the linearized expression of equation ( 9) can be expressed as follows: where w k is the process noise vector with a zero mean and covariance matrix Furthermore, the following is obtained [26]: into the left side of equation ( 9), the state space equation can be obtained: Only the partial wind-induced displacement and acceleration responses are measured.e measurement vector is expressed as follows: where x d denotes d × 1 wind-induced displacement responses and € x a denotes a × 1 wind-induced acceleration responses.
Using modal coordinate transformation theory, the measurement equation at the k-th time step can be expressed as follows: where T denotes the interstory displacements, D d is the d × n mapping matrix associated with the DOFs of the measured displacement, and D a is the a × n mapping matrix associated with the DOFs of the measured acceleration.
Furthermore, the following is obtained: Considering the measurement noise, the linearized measurement equation can be expressed as follows: where v k is the measurement noise vector at the k-th time step with zero mean and covariance Furthermore, the following is obtained: Shock and Vibration where (18)

Eigenvalue and Eigenvector
Sensitivity.e changes in the eigenvalues and eigenvectors of the system due to changes in system parameters are used to calculate matrix G Z  k|k and H Z k .For the case θ � [k 1 , k 2 , . . ., k n , α, β] T , the structural stiffness coefficient and damping coefficient are chosen as the unknown parameters.e eigenvalue problem of proportionally damped systems can be solved according to [27] as follows: where λ l � ω 2 l is the l-th eigenvalue, φ l is the l-th eigenvector, and θ j is the j-th unknown parameter.Assuming that m modes are used for joint state-parameter-input estimation for an n-degrees-of-freedom, the eigenvector derivative can be calculated based on Wang's method [28], as shown in the following equation: where us, the derivative of eigenvalue matrix Λ subject to structural parameter θ j can then be obtained based on equation (19): e sensitively of eigenvector matrix Φ to the parameter θ j can be given according to equation (20): e derivative of damping matrix Γ subject to structural parameter θ j can then be calculated as follows: For the case θ � [ω 1 , ω 2 , . . ., ω m , ξ 1 , ξ 2 , . . ., ξ m ] T , the natural frequency ω and damping ratio ξ are chosen as unknown parameters.
e eigenvalue and eigenvector sensitivity can be calculated based on [29] as follows: According to equation ( 4), the following equation can be obtained: in which c l � 2ω l ξ l is the l-th diagonal element in the Γ matrix.en, the derivative of the damping matrix subject to modal parameter θ j can then be calculated as follows: e derivative of eigenvalue matrix Λ subject to modal parameter θ j can then be obtained based on equation (25): e sensitively of eigenvector matrix Φ to the modal parameter θ j can be given: According to equations ( 10) and ( 31), the error of the predicted state estimate  Z k+1|k can be calculated as follows: where the coefficient matrices are the estimation errors of state and modal wind load at time t � k × Δt, respectively.
e covariance matrix P Z k+1|k related to the predicted state estimate  Z k+1|k can then be expressed as follows: where and

Modal Wind Load Estimation. Defining the innovation
), according to equation ( 16), the following equation can be obtained: where the error e k+1 is given by the following equation: As  Z k+1|k is unbiased, it follows from equation ( 35) that E[e k+1 ] � 0, and consequently according to equation (34), ] can be obtained.Assume that the form of the estimated modal wind load  d k+1 is as follows: erefore, E( According to equation ( 35), the covariance of the error e k+1 can be obtained as follows: Generally, P e k+1 ≠ cI, where c is a positive real number.is indicates that equation (34) does not satisfy homoscedasticity.erefore, the estimator given in equation ( 36) is not the minimum variance estimation of modal wind load according to the Gauss-Markov theorem [30].
To obtain the unbiased minimum variance estimation of modal wind load d k+1 , the optimal value of matrix S k+1 in equation ( 36) should be determined.Assume that the covariance matrix P e k+1 in equation ( 37) is positive definite (i.e., P e k+1 ≻ 0), an invertible matrix Ψ k+1 satisfying Ψ k+1 Ψ T k+1 � P e k+1 can be found.By premultiplying Ψ −1 k+1 to equation (34), one can obtain the following equation: Now the covariance E[(Ψ −1 k+1 e k+1 )(Ψ −1 k+1 e k+1 ) T ] � I, which satisfies homoscedasticity.Under the assumption that Ψ −1 k+1 H k+1 has full column rank, the unbiased minimum variance estimation of d k+1 can then be obtained based on the Gauss-Markov theorem [30] as follows: erefore, the optimal S k+1 is obtained: e error of the estimated modal wind load  d k+1 can be given based on equations ( 36) and (39): According to equation (41), the covariance matrix P d k+1 related to the estimated modal wind load  d k+1 is calculated as follows: where L k+1 is the gain matrix.e error ε Z k+1|k+1 of the updated state estimate  Z k+1|k+1 can be calculated according to equations ( 12) and (43): Based on Equations ( 44) and (45), the covariance matrix P Z k+1|k+1 related to the updated state estimate  Z k+1|k+1 can be obtained as follows: To obtain the unbiased minimum variance estimation of state Z k+1 , the optimal value of the gain matrix L k+1 should be determined.Based on the Lagrange multipliers method Shock and Vibration [31], the optimal gain matrix L k+1 can be calculated by minimizing the trace of P Z k+1|k+1 under the unbiased condition shown in equation (45): where By substituting equation ( 47) into (43), the unbiased minimum variance estimation of state Z k+1 can be calculated as follows: Similarly, substituting equation ( 47) into (46), the covariance matrix P Z k+1|k+1 related to  Z k+1|k+1 can be expressed as follows: Based on equations ( 41) and ( 44), the covariance matrices P Zd k+1 and P dZ k+1 can be obtained as follows: Now the estimated displacement response  x k+1 , velocity response  _ x k+1 , structural parameters  θ k+1 , and wind load  F k+1 at time t � (k + 1) × Δt can be calculated as follows: where  Y m,k+1 ,  _ Y m,k+1 , and  θ p,k+1 are the estimation of modal displacement response, modal velocity response, and structural parameters, respectively, which are obtained from the state estimate  Z k+1|k+1 .

Numerical Simulation
To verify the feasibility and accuracy of the proposed method, a ten-story shear building structure under wind load is considered.e mass coefficient of each floor is m � 10 5 kg, and the stiffness coefficient of each floor is k � 2.45 × 10 8 N/m e damping is assumed to be Rayleigh damping which is calculated as C � αM + βK with proportional coefficients of α � 0.557 and β � 0.0034.e corresponding damping ratio for the first two modes of vibration is approximate 5%.
e fluctuating wind speed is numerically simulated based on the autoregressive model method.
e power spectral density is Davenport spectral.
e vertical wind profile is taken as the power profile with an exponent of a � 0.22 and a reference height of z G � 10 m according to the Chinese National Load Code [32].e mean wind speed at the reference height is 10 m/s. Figure 1 shows the simulated fluctuating wind speed on the fifth and tenth floors.Figure 2 shows the comparison of the power spectral density between the simulated fluctuating wind speed and the Davenport spectral on the fifth and tenth floors.Figure 2 shows that the simulated power spectral density matches very well with the Davenport spectral.e wind load acting on the building structure is calculated according to [22].e air density is assumed to be ρ � 1.29 kg/m 3 .e drag coefficient is set to 1.3, and the orthogonal exposed wind area of each floor is set to 50 m 2 .

Joint State/Parameter/Wind Load Estimation Based on
Incomplete Measurements.In this section, the interstory displacement and acceleration taken as the "measurements" are calculated based on the Newmark-β method and superimposed with 2% root mean square (RMS) white noise.Only seven sets of measurements including seven accelerations and seven interstory displacements are used for the estimation.e locations of the measurements are listed in Table 1.e unknown structural parameters are structural stiffness coefficients on each floor and the two Rayleigh damping coefficients.
e initial value of the augmented state vector is e covariance matrix of process noise is Q � 10 −2 , and the covariance matrix of measurement noise is R � 10 −1 .Figure 3 shows the comparison of the estimated structural displacement responses with exact values on the fifth and tenth floors for time and frequency domains.Figure 4 shows the comparison of the estimated structural velocity responses with exact values on the fifth and tenth floors for time and frequency domains.e estimated curves coincide with the exact curves.
is indicates that the proposed method is capable of identifying structural responses.Furthermore, two obvious peaks at 1.185 Hz and 3.527 Hz can be obtained in Figures 3 and 4. e two peaks correspond to the first two translational natural frequency of the building.
Moreover, based on the proposed method, the unknown wind load is identified using incomplete measurements.e identified wind load and the relative errors between the identified wind load and the exact one in percentage on the fifth and tenth floors of the numerical model are plotted in Figure 5. Figure 5 shows that the iteration process cannot converge immediately, but after approximately 8 s, the relative errors are well converged to less than 5%.Table 2 shows the mean errors and the RMS errors between the identified wind loads and the exact ones on each floor.From Table 2, the maximum mean errors and RMS errors of the identified wind loads are 4.12% and 4.56%, respectively.is means that the identified wind loads match the exact ones very well after iterative convergence.
Figure 6 shows the estimated results of the fifth and tenth structural stiffness coefficients.It indicates that the stiffness 6 Shock and Vibration coefficients can converge to the real value within 2 s. e estimation results of Rayleigh damping coefficients α and β are plotted in Figure 7. Compared to the stiffness coefficient, the damping coefficient converges slightly slower and can converge to a real value within approximately 5 s.Table 3 lists the estimated values and errors of the structural parameters.e maximum error for the stiffness coefficient is 4.30%.e estimated errors for the damping coefficients α and β are 4.25% and 4.24%, respectively.is indicates that the proposed method is capable of identifying the structural parameters with incomplete measurements.

Effect of Sampling Duration.
To investigate the stability of the proposed method, the effect of the sampling duration on wind load and structural parameter estimation is discussed.e number and location of measurements are the same as in Table 1.e total duration of the measurements is set to 1 s, 10 s, 30 s, and 60 s.Table 4 shows the RMS errors of the estimated wind load on each floor under different sampling durations.e estimated results of the structural parameters are given in Table 5.From Table 4, the maximum RMS errors of wind load under 1 s, 10 s, 30 s, and 60 s are 14.23%, 7.68%, 4.58%, and 4.55%, respectively.From Table 5, the maximum errors of the structural parameters under 1 s, 10 s, 30 s, and 60 s are 32.78%,6.99%, 4.66%, and 4.25%, respectively.is indicates that the estimation errors decrease as the time duration increases.Table 5 shows that the estimation errors of the wind load and structural parameters are all less than 5% when the sampling duration reaches 30 s. Furthermore, comparing the estimated results between 30 s and 60 s shows that the estimation accuracy increases slightly.erefore, considering the computational costs, the sampling duration of 30 s is adequate for this numerical simulation.  is section discusses the influence of the number of measurements for joint state/parameter/wind load estimation.e considered set of measurements is four to nine, and each set of measurements includes interstory displacement and acceleration responses.Table 6 shows the number of measurements and the corresponding locations.
e time histories of the estimated wind loads on the fifth and tenth floors under different numbers of measurements are plotted in Figure 8. Figure 8 shows that the variation trends of the identified wind loads are identical to the exact ones.However, there is an obvious error between the estimated value and the exact one when the number of measurements is less than six.
e RMS errors of the identified wind load on each floor under different numbers of measurements are listed in Table 7. Table 7 shows that as the number of measurements increases, the RMS error of the estimated wind load decreases.e maximum error under six, five, and four sets of measurements is 5.33%, 7.08%, and 10.78%, respectively.is indicates that the estimation error increases significantly when the number of measurements is less than six.
e time histories of the estimated errors of the fifth and tenth structural stiffness coefficients for different numbers of measurements are shown in Figure 9. Figure 9 shows that at the beginning of the iteration are some fluctuations.However, after approximately 2 s, the iteration converges to the true value.Figure 10 plots the estimated errors of damping coefficients α and β. Figure 10 shows that both α and β have a large fluctuation at the beginning of the iteration, and they converge to the true value after approximately 10 s under more than seven measurements.However, the estimation results fluctuate around the true value when the number of measurements is less than seven.is indicates that damping coefficients are sensitive to the number of measurements and fluctuate around the true value if the number of measurements is insufficient.Table 8 gives the estimation results of the structural parameters with the number of measurements (i.e., from 4 to 9). e maximum errors of the estimated parameters under four to nine  Shock and Vibration measurements are 13.01%, 11.76%, 9.83%, 4.30%, 3.71%, and 3.63%, respectively.It indicates that as the number of measurements increases, the estimation accuracy of the structural parameters increases.e maximum errors of the structural stiffness coefficient and the damping coefficients α and β under six sets of measurements are 5.98%, 9.83%, and 5.33%, respectively.However, the estimation errors of each structural parameter are less than 5% under seven sets of measurements.e results indicate that, for this particular example, at least seven sets of measurements should be used.

Validation with Wind Tunnel Tests
In order to verify the effectiveness of the proposed method in practical engineering application, a tall building with 58 floors is considered as an example in this study.e size of the building is 234 m (height) × 39 m (width) × 39 m (length).e synchronous multipressure scanning system wind tunnel test for simultaneous estimation of state, unknown structural parameters, and wind load was conducted in wind environment wind tunnel laboratory at Harbin Institute of Technology (Shenzhen), China.According to Chinese National Load Code (GB 50009-2012) [32], the site around the test building is taken as terrain C.
e wind profile is regarded as power law with an exponent of 0.22.
e length scale of this wind tunnel test is 1 : 300. Figure 11 shows the simulated profiles of mean wind speed and turbulence intensity.
e model has the same length scale with that of wind field simulation, i.e., 1 : 300.e size of the model is 780 mm (height) × 130 mm (width) × 130 mm (length), as shown in Figure 12.In this wind tunnel test, wind direction was defined as an angle θ d from 0 °to 360 °with increments of 15 °, as shown in Figure 13.e total number of the pressure taps installed on the building was 441 with 98 on each side surface and 49 on top surface.e pressure data were collected using an microelectronic pressure scanner system made by PSI.
e data sampling frequency was 330 Hz, and the sampling length was 120 s. e wind speed at the top of the model with a 10-year return period is 39.6 m/s based on Chinese National Load Code [32].

Equivalent Model of the Tall Building.
For building structures with large degrees-of-freedom, it is complicated or even impossible to estimate the wind load and structural parameters strictly according to the prototype structure.Generally, a simplified equivalent model with less degreesof-freedom is needed in practical analysis.e equivalent model requires that the vibration responses of the two systems are same.
is means that the work-energy between the prototype structure and the equivalent structure is equal.Based on the method proposed in [33], the equivalent mass and equivalent stiffness can be obtained as follows: in which m(z) is the mass distribution of the prototype structure, k(z) is the stiffness distribution of the prototype structure, x(z, t) is the displacement of the prototype structure, M e is the equivalent mass, and K e is the equivalent stiffness.Based on the mode decomposition method, the displacement can be represented as follows: where Y j (t) is the j-th modal displacement and φ j (z) is the j-th mode shape vector.For tall building structures, the first mode is taken as the main mode of vibration.en, the equivalent mass and equivalent stiffness can be calculated as follows: Table 9 lists the energy contribution κ of the first five modal acceleration responses, modal velocity responses, and modal displacement responses, which are calculated according to the following equation: where σ i is the mean squared value of the responses calculated based on proper orthogonal decomposition [34].It can be seen that energy contribution of the first five modes is over 99%.erefore, the simplified equivalent model is chosen as a shear building structure with five degrees-offreedom.Assume that the simplified model quality is concentrated on the 12th, 24th, 36th, 48th, and 58th floors corresponding to the height of 50 m, 98 m, 146 m, 194 m, and 234 m, respectively.e equivalent mass coefficients and stiffness coefficients are calculated based on Equations ( 55) and (56), respectively, which are shown in Table 10.Table 11 gives the first five natural frequencies of the prototype structure and the equivalent model.It can be found that the relative errors of the first three natural frequencies are under 5%.Figure 14 shows the comparison of the first two mode shapes between the prototype structure and the equivalent structure.It can be found that the mode shapes of the equivalent structure are consistent with that of the prototype structure.e damping of this structure is considered as Rayleigh damping with damping ratio corresponding to the first two modes of vibration being 5%.

Wind Load and Structural Parameters Estimation.
e wind-induced responses of the tall building cannot be measured directly from the pressure measurement on a rigid model in the wind tunnel test.Hence, wind-induced responses including displacement, velocity, and acceleration responses are calculated based on the wind tunnel test results and the structural dynamic properties of the prototype building structure using the Newmark-β method.In order to
Figures 15 and 16 show the estimation results of structural displacement responses and velocity responses on the 24th and 58th floors, respectively.It can be found that the time histories of the estimated responses match very well with the exact responses.Meanwhile, the power spectrum of the estimated responses is found to agree well with those of the exact responses except slight differences in high frequency components (>0.8 Hz) of the responses on the 24th floor.is is mainly because (1) there no measurement information on the 24th floor, and the estimation results on the 24th floor is obtained from modal transformation and (2) the errors of the modal information between the equivalent model and prototype model are large in higher modes (i.e., the fourth and fifth modes).

Shock and Vibration
Displayed in Figure 17 is the comparison of the time histories of the equivalent wind load between the exact values and estimated values on the 24th floor and 58th floor.It can be observed that the estimated curves agree very well with the exact ones.e estimation errors of the equivalent wind load are listed in Table 12.From Table 12, the maximum mean error and RMS error of the identified wind load are 4.95% and 5.89%.Figure 18 shows the estimation results of the equivalent stiffness coefficients.It indicates that the stiffness coefficients can converge to the real value within 5 s. e estimation results of damping coefficients are shown in Figure 19, and it can be found that the damping  13. e errors of the stiffness coefficient are less than 5% except that on the 12th floor.e errors of the damping coefficient α and β are 3.8% and 2.11%, respectively.e comparison results indicate that the inverse method presented in this study is applicable to the practical engineering.

Conclusion
In this paper, a time domain joint state/parameter/wind load estimation method from incomplete measurements is proposed based on modal extended unbiased minimum variance estimation.e recursive procedure includes four parts: time update, modal wind load estimation, measurement update, and coordinate transformation.
e measurement responses include the interstory displacement and    e estimation quality of the proposed method is numerically validated by simultaneously identifying the structural state, parameters, and wind load of a ten-story shear building structure from incomplete measurements.e effects of crucial factors, including sampling duration and the number of measurements on the convergence and accuracy of the proposed method, are discussed.Moreover, a synchronous multipressure scanning system wind tunnel test on a 234 m tall building structure is used to demonstrate the proposed approach for real structure.A simplified equivalent model with five degrees-of-freedom is derived for experimental validation.Results indicate that the proposed method shows much potential as an alternative way for simultaneously identify the unknown structural parameters, wind load, and wind-induced responses from incomplete measurements.

Figure 1 :Figure 2 :
Figure 1: Simulated fluctuating wind speed on the (a) fifth floor and (b) tenth floor.

Figure 3 :
Figure 3: Comparison of the structural displacement responses for time and frequency domains.Displacement vs time on the (a) fifth floor and (b) tenth floor.Frequency vs power spectrum on the (c) fifth floor and (d) tenth floor.

Figure 4 :
Figure 4: Comparison of the structural velocity responses for time and frequency domains.Time vs velocity on the (a) fifth floor and (b) tenth floor.Frequency vs power spectrum on the (c) fifth floor and (d) tenth floor.

Figure 5 :
Figure 5: Comparison of the wind load time histories.Wind load on the (a) fifth floor and (b) tenth floor.Wind load errors on the (c) fifth floor and (d) tenth floor.

Figure 6 :Figure 7 :
Figure 6: Estimation results of unknown stiffness coefficients on the (a) fifth floor and (b) tenth floor.

Figure 8 :
Figure 8: Time histories of the estimated wind load for different numbers of measurements on the (a) fifth floor and (b) tenth floor.

Figure 9 :
Figure 9: Estimated stiffness coefficients for different numbers of measurements on the (a) fifth floor and (b) tenth floor.

Figure 11 :
Figure 11: Mean wind speed and turbulence intensity profiles: (a) mean wind speed; (b) turbulence intensity.

Figure 14 :Figure 15 :
Figure 14: Comparison of the mode shape: (a) the first mode shape; (b) the second mode shape.

Figure 16 :
Figure 16: Comparison of the structural velocity responses in direction X under wind direction of 180 °.Time history on the (a) 24th floor and (b) 58th floor.Power spectrum on the (c) 24th floor and (d) 58th floor.

Figure 17 :
Figure 17: Comparison of the equivalent wind load time histories in direction X under wind direction of 180 °on the (a) 24th floor and (b) 58th floor.

Figure 18 :
Figure 18: Estimation results of the equivalent stiffness coefficients.

Table 1 :
Number of measurements and corresponding locations.

Table 2 :
Mean errors and RMS errors for wind load estimation.

Table 3 :
Estimation values and errors of structural parameters.

Table 5 :
Structural parameter estimation results under different sampling durations.

Table 4 :
RMS errors of the estimated wind load under different sampling durations.

Table 6 :
Number of measurements and corresponding locations.

Table 7 :
RMS errors of wind load for different numbers of measurements.

Table 8 :
Estimation of structural parameters under different numbers of measurements.

Table 9 :
Energy contribution ratio of the first five modal responses.

Table 10 :
Equivalent mass and stiffness of the equivalent model.

Table 11 :
Comparison of the first five natural frequencies between the prototype model and the equivalent model.

Table 12 :
Mean errors and RMS errors for wind load estimation.

Table 13 :
Estimation values and errors of structural parameters.