# The sensitivity analysis of the translation and the rotation angle of the first-order mode shape of the joints in frame structures

## Yi Chen Yuan1, Lin Li2, Hongping Zhu3

1, 2, 3School of Civil Engineering and Mechanics, Huazhong University of Science and Technology, Wuhan 430074, China

2Corresponding author

Vibroengineering PROCEDIA, Vol. 5, 2015, p. 281-286.
Accepted 4 September 2015; published 18 September 2015

Copyright © 2015 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Views 57
Abstract.

When damage occurs, there are changes in the structural stiffness. This causes changes in the vibrational information, such as translation and rotation angle of the joints in the structure. This paper presents a sensitivity study of translation and slope of the first-order mode shape of the joints in the frame structure, and which will contribute to structural damage identification further. Starting from the equation of natural vibration of frame structure, and according to the characteristics of the stiffness matrix, derive sensitivity coefficient expression of translation and slope of first-order mode shape of the joints in the structure, and obtain the first-order relative sensitivity coefficient of the first story is always smaller than zero when single column of the first story damage occurs. The numerical analysis of the four stories three spans frame shows consistent results with the formula derivation.

Keywords: sensitivity analysis, translation, rotation angle.

#### 1. Introduction

The modal parameters of the structure (such as frequency, modal shape and modal damping) are the function of structural physical parameters (such as mass, stiffness and damping), so the change of structural physical parameters will lead to the changes of vibration modal parameters of structure. In order to correctly establish the relationship between the dynamic properties of structure and damage of structure, it is necessary to conduct the damage sensitivity analysis of the dynamic characteristics.

Zhao and DeWolf  derived sensitivity coefficients of natural frequency, vibration mode, and the modal flexibility, which was applied to the spring-mass model. It was concluded that the modal flexibility is more sensitive to the damage. Zhu et al. [2, 3] used optimization method to identify the damage of periodic structure based on sensitivity analysis of modal parameters of periodic structures. Li et al.  derived sensitivity coefficient expressions of structure natural frequency, mode shapes and mode shape slops based on the characteristics of shear frame structural stiffness matrix. Yang  presented a mixed sensitivity method that combines eigenvalue and flexibility sensitivity to identify structural damage. Chen  analyzed the effectiveness of damage localization method based on the change in the value of the slope of the first mode. Yam et al.  studied the sensitivity of static and dynamic parameters to damage occurring in plate-like structures and proposed corresponding damage indices.

The first-order mode shape slope is a good index of damage detection for shear frames, which can be used to detect the damage of story stiffness . As for the general frame, there is not only the translation but also rotation angle. Thus the damage of the beam and column will cause the changes of the slope and rotation angle.

In this paper sensitivity study of the translation and slope of the first-order mode shape of first story is conducted for further damage detection of general frame structures.

#### 2. Sensitivity analyses

In the sensitivity analysis, the sensitivity coefficient of the structure modal parameter is divided into absolute sensitivity coefficient and relative sensitivity coefficient. Absolute sensitivity coefficient is defined as the derivative of structure modal parameters with respect to structural physical parameters, and relative sensitivity coefficient is defined as the ratio of the rate of change of structure modal parameter to the rate of change of structural physical parameter.

#### 2.1. Sensitivity analyses of translation of mode shapes

The equilibrium equation for undamped structural vibration is:

(1)
$\left(K-{\omega }^{2}M\right){\phi }_{i}=0,$

where $K$ and $M$ are the stiffness and mass matrices, respectively, ${\omega }_{i}$ is the $i$th circular frequency, and ${\phi }_{i}$ is the $i$th mode shape corresponding to ${\omega }_{i}$, which is normalized to be the unit-mass mode shape, i.e.:

(2)
${\phi }_{i}^{T}M{\phi }_{i}=1,$

the derivative of Eq. (1) with respect to any beam or column stiffness, ${k}_{t}$, is:

(3)
$\left(\frac{\partial K}{\partial {k}_{t}}-2{\omega }_{i}\frac{\partial {\omega }_{i}}{\partial {k}_{t}}M\right){\phi }_{i}+\left(K-{\omega }_{i}^{2}M\right)\frac{\partial {\phi }_{i}}{\partial {k}_{t}}=0.$

And based on vibration theory, the mode shapes of a structure are independent of each other, and thus a complete set of mode shapes can be used to constitute the sensitivity coefficients of the $i$th mode shape. Therefore, the absolute sensitivity coefficients of the mode shapes can be defined as follows:

(4)
$\frac{\partial {\phi }_{i}}{\partial {k}_{t}}=\sum _{r=1}^{N}{\beta }_{r}{\phi }_{r},$

both sides of Eq. (3) are premultiplied by the transposed vector of the $r$th mode shape, ${\phi }_{r}^{T}$.Making use of the orthogonal properties of the mode shapes, Eq. (3) can be rewritten as:

(5)

because ${\phi }_{r}^{T}K={\phi }_{r}^{T}{\omega }_{r}^{2}M$, Eq. (5) is simplified:

(6)

For the first story column, substituting Eq. (4) into Eq. (6) and making use of the orthogonal properties of the mode shapes leads to:

(7)

where $n$ and $m$ are the number of stories and the number of spans respectively, $s$ is the $s$th column in the first story ,which is no more than $\left(m+1\right)$. $E$ and ${l}_{1}$ are respectively the elastic modulus and the length of the column. And when $r=i$:

(8)
${\beta }_{r}=0.$

Thus, the first-order absolute sensitivity coefficient of the mode shapes can be summarized as follows:

(9)
${S}_{it}^{M}=\frac{\partial {\phi }_{i}}{\partial {k}_{t}}=\sum _{r=1}^{N}{\beta }_{r}{\phi }_{r},$

where and the first-order relative sensitivity coefficient of the mode shapes can be summarized as follows:

(10)
${\stackrel{-}{S}}_{it}^{M}=\underset{\Delta {k}_{t}\to 0}{lim}\frac{\frac{\Delta {\phi }_{i}}{{\phi }_{i}}}{\frac{\Delta {k}_{t}}{{k}_{t}}}\approx \frac{\partial {\phi }_{i}}{\partial {k}_{t}}\cdot \frac{{k}_{t}}{{\phi }_{i}}.$

#### 2.2. Sensitivity analyses of mode shape slopes

For general frames, the mode shape slope of the $j$th story in the $i$th mode is defined as:

(11)
${\phi \text{'}}_{ij}=\frac{{\phi }_{ij}-{\phi }_{i\left(j-1\right)}}{{h}_{j}},$

where ${h}_{j}$ is the height of the $j$th story. Taking the first derivatives of Eq. (11) with respect to ${k}_{t}$ leads to:

(12)
${S}_{ij,t}^{S}=\frac{\partial {\phi \text{'}}_{ij}}{\partial {k}_{t}}=\frac{1}{{h}_{j}}\left(\frac{\partial {\phi }_{ij}}{\partial {k}_{t}}-\frac{\partial {\phi }_{i\left(j-1\right)}}{\partial {k}_{t}}\right),$

where ${S}_{ij,t}^{S}$ is the first-order absolute sensitivity coefficient of the mode shape slope of the $j$th story in the $i$th mode with respect to the $t$th column or beam stiffness. And the first--order relative sensitivity coefficient:

(13)
${\stackrel{-}{S}}_{ij,t}^{S}=\underset{\Delta {k}_{t}\to 0}{lim}\frac{\frac{\Delta {\phi \text{'}}_{ij}}{{\phi \text{'}}_{ij}}}{\frac{\Delta {k}_{t}}{{k}_{t}}}\approx \frac{\partial {{\phi }^{\text{'}}}_{ij}}{\partial {k}_{t}}\cdot \frac{{k}_{t}}{{{\phi }^{\text{'}}}_{ij}}=\frac{{k}_{t}}{\partial {\phi }_{ij}-\partial {\phi }_{i\left(j-1\right)}}\left(\frac{\partial {\phi }_{ij}}{\partial {k}_{t}}-\frac{\partial {\phi }_{i\left(j-1\right)}}{\partial {k}_{t}}\right),$

when $j=1$:

(14)
${\stackrel{-}{S}}_{i1,t}^{S}\approx \frac{\partial {\phi }_{i1}}{\partial {k}_{t}}\cdot \frac{{k}_{t}}{{\phi }_{i1}},$

when $i=1$:

(15)
$\frac{\partial {\phi }_{11}}{\partial {k}_{t}}\cdot \frac{{k}_{t}}{{\phi }_{11}}=\frac{{k}_{t}}{{\phi }_{11}}\cdot \sum _{r=2}^{N}{\beta }_{r}{\phi }_{r1}=\frac{{k}_{t}}{{\phi }_{11}}\cdot \left({\beta }_{2}{\phi }_{21}+{\beta }_{3}{\phi }_{31}+{\beta }_{4}{\phi }_{41}+...\right).$

According to the analysis of the examples, we can get nearly accurate answer of $\frac{\partial {\phi }_{11}}{\partial {k}_{t}}$ when $r=n$. As for ${\beta }_{r}{\phi }_{r1}$, suppose:

(16)
${\beta }_{r}{\phi }_{r1}<0,$
(17)
${\beta }_{r}{\phi }_{r1}=\frac{2E\left(6{\phi }_{1,1}{\phi }_{r,1}+3{\phi }_{1,1}{\phi }_{r,n+s}{l}_{1}+3{\phi }_{r,1}{\phi }_{1,n+s}{l}_{1}+2{\phi }_{1,n+s}{\phi }_{r,n+s}{{l}_{1}}^{2}\right)}{{l}_{1}^{3}\left(-{\omega }_{r}^{2}+{\omega }_{1}^{2}\right)}\cdot {\phi }_{r,1}<0.$

Therefore:

(18)
$\frac{\left(3{\phi }_{1,1}{\phi }_{r,n+s}{l}_{1}+2{\phi }_{1,n+s}{\phi }_{r,n+s}{{l}_{1}}^{2}\right){\phi }_{r,1}}{6{\phi }_{1,1}{\phi }_{r,1}^{2}+3{\phi }_{r,1}^{2}{\phi }_{1,n+s}{l}_{1}}>-1,$

when ${\phi }_{r,n+t}$, ${\mathrm{\phi }}_{\mathrm{r},1}$ symbol are the same, inequality (18) is correct. And when ${\phi }_{r,n+t}$, ${\phi }_{r,1}$ have different symbol:

(19)
$\frac{\left(3{\phi }_{1,1}{l}_{1}+2{\phi }_{1,n+s}{{l}_{1}}^{2}\right)}{6{\phi }_{1,1}+3{\phi }_{1,n+s}{l}_{1}}<-\frac{{\phi }_{r,1}}{{\phi }_{r,n+s}}.$

According to the analysis of the examples, $|\frac{{\phi }_{1,1}}{{\phi }_{1,n+s}{l}_{1}}|$ is the minimum in all of $\left|\frac{{\phi }_{i,1}}{{\phi }_{i,n+s}{l}_{1}}\right|$ ($i=$ 1, 2, 3, ...). And making $-\frac{{\phi }_{r,1}}{{\phi }_{r,n+s}{l}_{1}}$ equals to $\frac{{\phi }_{1,1}}{{\phi }_{1,n+s}{l}_{1}}$, inequality (19) yields:

(20)
$\frac{\left(3\frac{{\phi }_{1,1}}{{\phi }_{1,n+s}{l}_{1}}+2\right)}{6\frac{{\phi }_{1,1}}{{\phi }_{1,n+s}{l}_{1}}+3}<\frac{{\phi }_{1,1}}{{\phi }_{1,n+s}{l}_{1}},$
(21)
$\frac{{\phi }_{1,1}}{{\phi }_{1,n+s}{l}_{1}}>\frac{\sqrt{3}}{3}\approx 0.577.$

Thus, inequality (16) is necessarily true when $\frac{{\phi }_{1,1}}{{\phi }_{1,n+s}{l}_{1}}>\frac{\sqrt{3}}{3}$. And according to the analysis of cantilever column ,taking concrete frame columns minimum section 300 mm×300 mm, and the corresponding $\frac{{\phi }_{1,1}}{{\phi }_{1,n+s}{l}_{1}}$ is 0.669, which is more than $\frac{\sqrt{3}}{3}$. So inequality (16) is true, that is the first-order relative sensitivity coefficient of the first story is always smaller than zero when the first story column is damaged.

#### 3. Numerical example

The proposed method is demonstrated using numerical example of three spans four stories frame, as shown in Fig. 1, which is referred to as model 1 hereinafter. All of the column and beam length is 3 m, and the cross section of the column is 500 mm×500 mm and cross section of the beam is 300 mm×500 mm. Other parameter information is in the Table 1.

Fig. 1. Three spans four stories frame Table 1. Parameter information of the frame

 Elastic Modulus Poisson’s ratio Density 2.486E+10 0.2 2400 kg/m3

When the frame is damaged, stiffness of the column or beam will reduce. In the study, the reducing of the stiffness is simulated by the reducing of moment of inertia of column or beam. The following discussing the sensitivity of the first-order mode shape slope of first story when the first story column is damaged.

Assumes that the damage degree is 20 %, and the translation of the first order mode shape before and after the damage is shown in Table 2, and the slope of the first order mode shape before and after the damage is shown in Table 3.

Table 2. Translation of the first order mode shape before and after the damage

 Story The translation before damage The translation after damage 1 0.001817 0.001849 2 0.004335 0.004349 3 0.006258 0.006252 4 0.007249 0.007231

Table 3. Slope of the first order mode shape before and after the damage

 Story The slope before damage The slope after damage 1 0.000606 0.000616 2 0.000839 0.000833 3 0.000641 0.000634 4 0.000330 0.000326

The Table 3 shows that when the first story column is damaged the slope of the first order mode shape of first story will increase, which also means that the first-order relative sensitivity coefficient of mode shape slope of first story is smaller than zero.

From the above numerical example we can obtain that the first-order relative sensitivity coefficient of the first story is smaller than zero when the first story column is damaged, and this conclusion is the same with the results of the formula derivation.

#### 4. Conclusions

Starting from the equation of natural vibration of frame structure, and according to the characteristics of the stiffness matrix, the paper has derived sensitivity coefficient expression of translation and slope of first-order mode shape of the joints in the structure. And some important observations can be pointed out:

1. When the single column of the first story damage occurs, the first-order relative sensitivity coefficient of the first story is always smaller than zero.

2. The damage of column of the first story can be detected by observing the changes of first-order mode shape slope of the first story.

#### Acknowledgements

The authors would like to express their most sincere thanks for the support of the National Natural Science Foundation of China (Grant No. 51278214).

1. Jun Zhao, De Wolf John T. Sensitivity study for vibrational parameters used in damage detection. Journal of Structural Engineering, Vol. 125, Issue 4, 1999, p. 410-416. [Search CrossRef]
2. Zhu Hongping, Li Lin, Wang Dansheng Damage assessment in periodic structures from measured natural frequencies by a sensitivity and transfer matrix-based method. Structural Engineering and Mechanics, Vol. 16, Issue 1, 2003, p. 17-34. [Search CrossRef]
3. Zhu Hongping, Xu Y. L. Damage detection of mono-coupled periodic structures based on sensitivity of modal parameters. Journal of Sound and Vibration, Vol. 285, Issue 1, 2005, p. 365-390. [Search CrossRef]
4. Li Lin, Zhu Hongping, Chen Huan Sensitivity analysis of shear frame damage modal parameters. Journal of Huazhong University of Science and Technology, Vol. 24, Issue 3, 2007, p. 35-39. [Search CrossRef]
5. Yang Q. W. A mixed sensitivity method for structural damage detection. Communications in Numerical Methods in Engineering, Vol. 25, Issue 4, 2009, p. 381-389. [Search CrossRef]
6. Chen Huan Study on Damage Identification of Frame Structures Based on the Slope of First Order Mode Shape. Huazhong University of Science and Technology, Wuhan, 2008. [Search CrossRef]
7. Yam L. H., Li Y. Y., Wong Wo Sensitivity studies of parameters for damage detection of plate-like structures using static and dynamic approaches. Engineering Structures, Vol. 24, Issue 11, 2002, p. 1465-1475. [Search CrossRef]
8. Zhu Hongping, Li Lin, He Xiaoqiao Damage detection method for shear buildings using the changes in the first mode shape slopes. Computers and Structures, Vol. 89, Issues 9-10, 2011, p. 733-743. [Search CrossRef]