Stiffness and Damping Simultaneous Identification Based on Extrapolation of Unrecorded Response to Ground Motion

A new method is proposed for stiffness-damping simultaneous identification in a building structure when responses are recorded at only two stories, i.e. the base and the second floor. Uniqueness of stiffness and damping is guaranteed if the transfer function of the absolute acceleration (or displacement) at the second floor relative to the base is given. Stiffness and damping distributions can then be obtained by minimizing the error between two transfer functions, i.e. the ratio of records at the above-mentioned two floors in the frequency domain and the corresponding theoretical value of the system with assumed stiffness and damping. The stiffness and damping for the theoretical transfer function are evaluated by using identification functions through the extrapolation of unrecorded story responses. The validity of the proposed method is investigated through numerical examples. and use of recurrence equations and to overcome the difficulty in application to actually recorded data. The validity and effectiveness of the proposed methods are demonstrated through numerical examples using theoretical data in the frequency domain and time-history data. Identification Method using Limit Manipulation Identification of Stiffness and Damping Coefficient The methods developed in [12,15] are explained here (Section 2) for the development of new theories in Sections 3 and 4. Consider an N-story shear building model, as shown in Figure 1, with viscous damping. The node and element numbers are defined from the top. The jth node and the jth element from the top are called the “node j” and the “element j”, respectively. Let mj and kj denote the mass of the node j and the story stiffness of the element j and let cj be the viscous damping coefficient in the element j. The frequency-domain equations of motion for this shear building model subjected to the horizontal base acceleration ( )  z t can be expressed by ( ) ( ) ( ) ( ) ω ω ω ω = Z A U F (1) ( ) ω U and ( ) ω Z are the Fourier-transforms of the nodal absolute horizontal displacements ( ) t u and the base displacement ( ) z t , respectively. The vector ( ) ω F in Eq.(1) indicates ( ) {0 0 i } ω ω = +  T N N c k F and the symbol ( ) ⋅ T means the transpose of a vector. The matrix ( ) ω A in Eq.(1) is defined by 2 ( ) i ω ω ω = − + + A M C K (2) Citation: Kojima K, Kuwabara M, Yoshitomi S, Fujita K, Takewaki I (2014) Stiffness and Damping Simultaneous Identification Based on Extrapolation of Unrecorded Response to Ground Motion. J Archit Eng Tech 3: 117. doi:10.4172/2168-9717.1000117

In general, responses need to be measured at every floor above the target story for identification to evaluate the story shear force required in the stiffness-damping evaluation. Therefore, in order to identify physical parameters at the lowest story or at every story, responses at all stories are necessary. It is rather impractical to measure responses at every story from the viewpoint of instrumentation and data management, especially in high-rise buildings. However, limited location of measurement causes problems of non-existence or non-uniqueness of solution [12,13]. To avoid this difficulty, several approaches have been proposed. One approach is to alleviate the degree of ill-posedness using techniques such as the regularization technique [14]. Another approach is to certify the uniqueness of solution for the limited location of measurement.
In the latter approach, theories for unique identification of local stiffness and viscous damping coefficients at a specified target story are proposed based on the local response records at the floors just above and below a target story [12,15]. Furthermore, some methods for elimination of noise effects have been proposed [16,17] based on these theories. In these papers, the mathematical limit manipulation toward an infinite frequency or zero frequency has been introduced and simultaneous identification of local stiffness and damping can be performed in a unified manner and only an FFT technique is needed.
As for the identification of physical parameters at whole stories using limited local records, it has been shown that unique identification can be performed using both the limit manipulation toward an infinite frequency and recurrence equations when records at the base and the second floor are obtained [12]. In this paper, two new SI methods are proposed when the responses at the base and the second floor are obtained. These methods enable one to avoid the limit manipulation

Identification of Stiffness and Damping Coefficient
The methods developed in [12,15] are explained here (Section 2) for the development of new theories in Sections 3 and 4.
Consider an N-story shear building model, as shown in Figure 1, with viscous damping. The node and element numbers are defined from the top. The jth node and the jth element from the top are called the "node j" and the "element j", respectively. Let m j and k j denote the mass of the node j and the story stiffness of the element j and let c j be the viscous damping coefficient in the element j. The frequency-domain equations of motion for this shear building model subjected to the horizontal base acceleration ( )  z t can be expressed by where, i indicates the imaginary unit.
The jth component ( ) It should be noted that the matrix for which ∆ j is defined has a tri-diagonal property only for = j N and ∆ N can be expressed rather compactly. This property is closely related to the fact that the story stiffness of the element N and its material damping ratio can be identified uniquely when the acceleration ( )  N u t of the node N is observed in addition to the base acceleration ( )  z t . Substitution of Eq.(5) into Eq.(4) provides the following relation.
After some manipulations [15], the following relations can be derived.
The stiffness can be obtained from Eq.(8) and the damping coefficient can be derived from Eq.(9) when the masses m j are given and the floor accelerations or displacements are recorded (see Figure 1a). In Eq.

Uniqueness of identification of structural parameters including unrecorded story
Uniqueness of stiffness and damping distributions is discussed here when there are some stories of unrecorded responses (see Figure 1b). The 2-story models with different sensor locations as shown in Figure 2 (a) SI method developed in [15] (b) SI method developed in [12] ground motion  are used and stiffness and damper distributions are derived analytically for both models.
Sensor location guaranteeing uniqueness of identification: As for the sensor location as shown in Figure 2a, i.e. sensors at the base and the second floor, it can be proved that a unique identification of stiffness and damping distributions is possible based on the recorded response data [12,15].
From Eq.(4), 1 ( ) ω U and 2 ( ) ω U can be expressed as Here, two models with different sets of stiffness and damping coefficients, i.e.
The left-hand side is equal to the right-hand side for every ω if and only if the following relations hold.
The following solution is then obtained which satisfies Eq.(17).
This means that, if the transfer function 2 ( ) ω H of the lowest story is given, there exists only one corresponding set of stiffness and damping coefficients. In this sensor location, the uniqueness of stiffness and damper coefficient distributions holds. This result holds also for the N-story structure [12]. In this paper, this sensor location is used later.
Sensor location without uniqueness: As for the sensor location shown in Figure 2b, i.e. sensors at the base and the top floor, it can be proved that the unique identification of stiffness and damping distributions is not possible from the recorded response data.
As in the previous section, two sets of stiffness and damping coefficients, i.e. 1 The left-hand side is equal to the right-hand side for every ω if and only if the following relations hold.
The following two solutions are obtained which satisfy the relations of Eq.(20).
where, 1 1 This means that, if the transfer function 1 ( ) ω H of the top story is given, there exist two corresponding sets of stiffness and damping coefficients. It is therefore concluded that, in this sensor location, the uniqueness of stiffness and damper coefficient distributions is not guaranteed. In general, if the transfer function 1 ( ) ω H of the top story is given in the N-story model, it can be proved that there exists N! stiffness and damping coefficient distributions.

Identification Method Based on Extrapolation of Unrecorded Data (Method 1)
In this section, a new identification method using extrapolation of unrecorded data is proposed when responses are recorded at the base and the second floor (see Figure 1b). This method is called "Method 1".

Extrapolation of unrecorded data
A method is proposed for generating unrecorded data based on extrapolation. Here, a ratio L of the interstory drift to that of the lowest story is introduced. The conceptual figure is shown in Figure 3. Let define the following ratio j L .

Identification method using extrapolated data
By using extrapolated response data in the previous section, unique identification of stiffness and damping coefficient distributions can be performed. The outline of the identification is shown in Figure 4. This method identifies the stiffness and damping coefficient distributions by determining ratios Ground motion  are assumed, all the displacement can be estimated from Eq.(26), and by using these extrapolated data, identification can be performed by using Eqs.(7)- (9). When all the stiffness and damping coefficients can be obtained by substituting these values into Eq.(1), the transfer function of the lowest story can be estimated.
In general, this estimated transfer function constructed from extrapolated data is different from the actual transfer function constructed from recorded data. According to the discussion in Section 2.2.1, these two functions are equal if and only if the true distributions of stiffness and damping coefficients are used. When actual recorded data are used, it is difficult to guarantee the equality of these two functions exactly. Therefore, let us define a function ( ) f L in Eq.(29) to estimate the error between these two transfer functions. Here, As for the 2-story model shown in Table 1 . When the extrapolation is performed using these exact values, identification functions are estimated as shown in Figure 7. From Figure 7, the limit values of identification functions at frequency 0    . Therefore, in the numerical examination below, it is assumed that k N and c N can be identified exactly.
The result of the identification for the 2-story model based on the "Method 1" is shown in Table 2. It can be observed that the identification based on "Method 1" has been performed effectively. Figure 9 shows the transfer functions of the lowest story using given data ( ) The above-mentioned 2-story model has proportional damping. The "Method 1" is then applied to a 5-story model with non-proportional damping shown in Table 3. The result of the identification is shown in Table 5 which indicates that the "Method 1" can identify the stiffness and damping coefficients of all the stories with high accuracy using only two given response data.        Figure 11 indicates an identification function in the lowest story of the model shown in Table  1. This figure was obtained by applying ( ) to Eq.(8). This figure shows that the identification function using time-history data converges to 0 at frequency 0 because of the unfavorable effect of noise. Therefore, some approximation methods may be necessary to estimate the limit value and it is difficult to use time-history data directly in "Method 1". In the next section, another identification method without the limit value is proposed.

Identification method without limit manipulation
The second approach for identification is proposed in this section when responses are recorded at the base and the second floor. This method is called "Method 2". This method does not use extrapolation and directly identify the physical parameters, i.e. stiffness and damping coefficients.
Let us introduce the following ratios X kj and X cj .
where, k N and c N can be identified by applying recorded data to Eqs.(7)- (9). X kN and X cN are equal to 1 and can be specified independently.
A detailed algorithm of the proposed method is shown in Figure  12. When the ratios are assumed, all the stiffness and damping coefficients can be obtained. By substituting these values into Eq.(1), the transfer function of the lowest story can be estimated. As discussed in Section 3.2, a function ( ) f X is defined by Eq.(32) to estimate the error between two transfer functions. Here, As for the "Method 2", the design variables and their initial values are set as and the ratios X kN , X cN of the lowest story are fixed to 1.0. In the numerical examination below, it is assumed that k N and c N can be identified exactly and the optimization is performed by the same procedure as in Section 3.3.  Using theoretical data Using simulated data   Figure 12: Flow chart of the proposed identification method without extrapolation.
The result of identification for the 2-story model (Table 1) is shown in Table 5. The result of identification for the 5-story model (Table 3) is shown in Table 6. These tables indicate that "Method 2" can identify stiffness and damping coefficients of all the stories based on only two given response data with a high accuracy comparable to "Method 1". Investigation using simulated time-history data: The validity of "Method 2" is investigated here through numerical examples using simulated time-history data. An input ( )  z t at the base is generated as a band-limited white noise. A time-history response analysis is performed for ( )  z t . ( )  N u t and ( )  z t are then used as known recorded data. It is assumed here that k N and c N can be identified exactly from the recorded data by using Eqs. The results of identification for the 2-story and 5-story models are shown in Table 7. "Method 2" has a better performance because of avoiding the limit manipulation. The accuracy of the 2-story model is higher than that of the 5-story model and the value of ( ) f X are greater than that in Tables 4,5 and 6 using theoretical data in the frequency domain. Figures 13 and 14 show the target transfer function and the identified transfer function at the lowest story corresponding to the 2-story and 5-story models, respectively. In both figures, the difference between two transfer functions increases in higher frequency range. To investigate this error, two transfer functions of the 2-story model estimated by a different approach is shown in Figure 15, i.e. the transfer function estimated using theoretical data in the frequency domain and the transfer function estimated using time-history data. Figure 15 indicates that the difference between two functions increase in higher frequency range. The validity of "Method 2" depends on the assumption that these two functions are considered to be equal with sufficient accuracy. By improving the performance in this accuracy problem, the identification based on "Method 2" is expected to be useful.

Conclusions
In order to overcome the difficulties in the past work [12] in which stiffness and damping coefficients are identified uniquely when     responses at the base and the second floor are obtained, more practical and stable physical-parameter SI methods called "Method 1" and "Method 2" have been proposed in this paper.
The advantageous features of these methods can be summarized as follows: i. These methods enable identification by minimizing the error between the target transfer function and an assumed transfer function. The target transfer function is estimated directly from recorded data. The assumed transfer functions are estimated in a different way in each method. This approach is based on the uniqueness of identification proved in the previous paper [12] and it enables one to avoid the process to solve the recurrence equation.
ii. "Method 1" uses the extrapolation of unrecorded responses in the frequency domain. By introducing the ratios among interstory drifts, unrecorded data (floor acceleration or displacement) can be extrapolated. Applying the previously proposed SI method with limit manipulation [15,16] to these extrapolated data, physical parameters can be obtained at all stories. Based on these assumed physical parameters, the transfer functions can be calculated. Then optimization is performed to minimize the error between the target transfer functions and assumed transfer functions with the ratios among interstory drifts as the design parameters.
iii. The validity of the proposed method, "Method 1", has been investigated numerically. When the theoretical data in the frequency domain is used, SI is performed with sufficiently high accuracy. However, when the time-history data are used, the limit manipulation is difficult to perform because of the effect of noise in the low frequency range.
iv. To overcome the difficulty in the limit manipulation, "Method 2" has been proposed. This method enables one to avoid the extrapolation of data and the limit manipulation by regarding directly the physical parameters as the design variables. Optimization has been introduced to minimize the error between the target transfer function and an assumed transfer function with the physical parameters as the design parameters.
v. The validity of the proposed method, "Method 2", has been investigated. When the theoretical data in the frequency domain is used, SI can be performed with the accuracy comparable to the "Method 1". Furthermore, when the time-history data are used, this method is able to perform SI very well. However, the deterioration of the accuracy is caused by the error between two transfer functions, i.e. the transfer function estimated using theoretical data in the frequency domain and the transfer function estimated using time-history data.