Next Article in Journal
A Multi-Year Study on Rice Morphological Parameter Estimation with X-Band Polsar Data
Previous Article in Journal
Endoscopic Laser-Based 3D Imaging for Functional Voice Diagnostics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Statistical Model Updating Method of Beam Structures with Random Parameters under Static Load

School of Civil Engineering and Architecture, Wuhan University of Technology, Wuhan 430070, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2017, 7(6), 601; https://doi.org/10.3390/app7060601
Submission received: 8 April 2017 / Revised: 31 May 2017 / Accepted: 5 June 2017 / Published: 9 June 2017

Abstract

:
This paper presents a new statistical model updating method of beam structures with random parameters under static load. The new updating method considers structural parameters and measurement errors to be random. To reduce the unmeasured degrees of freedom in the finite element model, a static condensation technique is used in this method. A statistical model updating equation with respect to element updated factors is established afterwards. The element updated factors are expanded as random multivariate power series. Using a high-order perturbation technique, the statistical model updating equation can be solved to obtain the coefficients of the power series expansions of the element updated factors. The results of two numerical examples show that for the solution of the statistical model updating equation, the accuracy of the proposed method agrees with that of the Monte Carlo simulation method very well. The static responses obtained by the updated finite element model coincide with the measured results very well. Finally, a series of static load tests of the concrete beam are conducted to testify the effectiveness of the proposed method.

1. Introduction

In recent years, structural health monitoring and safety assessment have been a hot research topic which has attracted the attention of numerous researchers. In structural health monitoring, the damage identification of the structure is a key point. However, to realize the damage identification, an initial referenced model must be established. Due to the influence of various factors, e.g., environment, construction styles etc., there is a difference between the finite element simulation model and the actual structure. To reduce the impact of these factors, it is necessary to update the simulation model to make it coincide with the actual structure using measured data. In this regard, one can consider the structural responses (e.g., displacements, strains, and dynamic characteristics, etc.) as the indicators, and update the parameters and boundary conditions of the initial simulation model based on the measured data. This process is generally named as the model updating. In this way, an eligible referenced model can be obtained so as to further achieve the structural damage identification. Therefore, the model updating will provide a reliable basis for the structural safety assessment [1,2,3,4].
Much research has been conducted on the model updating method in accordance with the structural measured data. The measured data may be displacements, strains, or dynamic characteristics of the structures. Since the static test equipment required is relatively cheap and the static displacements or strains of structures can be measured economically and accurately, the static model updating methods have attracted the attention of many researchers. For instance, Sanayei et al. [5,6,7] presented the analytical methods to identify the structural element stiffness using the applied forces and measured displacements. Subsequently, the static displacements and static strains are used to successfully evaluate the unknown stiffness parameters of the structural components. Liu et al. [8] proposed a superelement-based virtual distortion method to improve the efficiency of the finite element model updating of large-scaled bridges by using static information. Ren et al. [9] presented a method for updating a finite-element model based on the measured static responses of structures with the aid of the response surface methodology. Li et al. [10] proposed a method based on the meta-model for updating the finite element model of bridges by using measured static and dynamic results. Sanayei et al. [11] collected strain data during the load test to calibrate a detailed baseline finite element model in an effort to represent the 3D system behavior of a bridge.
For the above research, the uncertainty of the measurement errors is seldom involved. Considering the measurement errors as random quantities, many statistical approaches have been developed to update the parameters of structures based on the measured data. These approaches include the stochastic perturbation methods [3,12,13,14,15,16], the Monte Carlo simulation methods [17,18], the Bayesian updating methods [19,20,21,22,23,24,25,26] and so on. For example, Jacquelin et al. [12] proposed a random matrix approach to derive the closed-form expressions for the mean matrix and the covariance matrix of the updated stiffness matrix by the perturbation technique. Husain et al. [13] considered the statistical properties of experimental data and updating parameters as random variables, and used the perturbation method to update the parameters. Combined with the sensitive method, Hua et al. use a Monte Carlo simulation method to solve the updating parameters [18]. Although the results from the Monte Carlo simulation methods are traditionally regarded as the exact solution of random problems, these methods generally require much greater computational efforts to obtain an accurate solution, and are impractical in the case of large-scale problems. Zárate et al. [23] studied the behaviour of the Bayesian updating framework when both static and dynamic data are used to update the model. However, all these studies mainly focus on using dynamic data to update structural models. Although some researchers have tried to use different random solution methods based on the finite element model to conduct statistical model updating, few investigations involve the direct use of random finite element methods to perform this work with static data. In light of the evidence, it is clear that further research is needed in this area.
In this paper, a new statistical model updating method of beam structures with random parameters under static load is proposed. This new model updating method considers structural parameters and measurement errors as random quantities. A static condensation technique is used to reduce the unmeasured degrees of freedom in the random finite element model of the beam. Then, a statistical model updating equation with respect to element updated factors is established. The element updated factors are expanded as random multivariate power series. Using the high-order perturbation technique, the statistical model updating equation can be solved to determine the coefficients of the power series expansions of the element updated factors. The results of two numerical examples show that for the solution of the statistical model updating equation, the accuracy of the proposed method agrees with that of the Monte Carlo simulation method very well. The static responses obtained by the updated finite element model coincide with the measured results very well. Finally, static load tests of a concrete beam are conducted to testify the effectiveness of the proposed method.

2. New Statistical Model Updating Method

2.1. Initial Equilibrium Equations

Consider one N degrees of freedom of beam structure under the external static load. The static equilibrium equation of the beam structure at initial simulation state is written as:
K a x a = F
where K a is an N × N dimensional stiffness matrix, x a is an N × 1 dimensional displacement response vector of the structure under N × 1 dimensional static load vector F .
Considering the translational and rotational degrees of freedom of the structure, Equation (1) can be rewritten as:
[ K a t t K a t θ K a θ t K a θ θ ] { x a t x a θ } = { F a t 0 }
where K a t t and K a θ t are respectively sub-matrices, which are related to the translational degrees of freedom of the initial stiffness matrix; K a , K a t θ and K a θ θ are sub-matrices of K a with respect to the rotation degrees of freedom, which are not measured in this paper. x a t and x a θ are displacement vectors related to the translational and rotation degrees of freedom of the initial structure respectively. F a t is the concentrated vertical force vector.
According to the static condensation method, which is also used in literature [5] by Sanayei et al., the rotational degrees of freedom of the stiffness matrix in the static equilibrium Equation (2) can be eliminated. In this way, Equation (2) can be condensed as:
K a t x a t = F a t
where K a t = K a t t K a t θ K a θ θ 1 K a θ t .
Of course, in the same manner, other unmeasured displacements can be removed.

2.2. Statistical Model Updating Equation

Without losing the generality, it is assumed here that the mass of the structural simulation model is invariable compared with the actual structure [27], and the discrepancy between the simulation model and the actual beam structure is mainly due to the variation of the structural stiffness, which is represented as:
Δ K = i = 1 n α i K i
where n is the number of the structural elements; α i is the element updated factor of the ith element of the beam structure, which is the variation ratio of a structural parameter such as the bending rigidity of the beam; K i is an N × N dimensional expanded matrix of the ith element stiffness matrix, where all elements in expanded parts are zero.
In regard to the stiffness matrix of the initial simulation model, the updated stiffness matrix, K m can be expressed as follows:
K m = K a + i = 1 n α i K i
Similar to the initial simulation model, the updated stiffness matrix can be condensed as:
K m t = K m t t K m t θ K m θ θ 1 K m θ t
where K m t t , K m t θ , K m θ t and K m θ θ are sub-matrices of K m , which are related to the translational and rotational degrees of freedom of the updated finite element model respectively like the situation at the initial simulation state.
Equations (5) and (6) can be rewritten as:
K m t = ( K a t t + i = 1 n α i K i t t ) ( K a t θ + i = 1 n α i K i t θ ) ( K a θ θ + i = 1 n α i K i θ θ ) 1 ( K a θ t + i = 1 n α i K i θ t )
where K i t t , K i t θ , K i θ t and K i θ θ are sub-matrices of the expanded matrix K i , which corresponds to the sub-matrices of the updated stiffness matrix.
Taking the first-order partial derivatives of K m t with respect to α i , there is:
K m t α i | α i = 0 = K i t t K i t θ K a θ θ 1 K a θ t K a t θ K a θ θ 1 K i θ t + K a t θ K a θ θ 1 K i θ θ K a θ θ 1 K a θ t
It is worth pointing out that in the sensitive methods as presented by Mottershead et al. in the literature [28,29], the calculation steps similar to the above are used to obtain the partial derivatives of the mass matrix and stiffness matrix with respect to updating parameters based on the first-order Taylor expansion.
Since the applied load is the same for the initial model and the updated model, the product of the updated stiffness matrix K m t and the measured displacement vector x m t should equal to the concentrated force vector F a t . Note that this principle is different from that presented in literature [5] by Sanayei et al., where a displacement equation error function is defined and minimized to determine the updating parameters. Using Equation (3), the following equation can be obtained:
K a t x a t = ( K m t | α = 0 + i = 1 n α i K m t α i | α i = 0 ) x m t
where K m t | α = 0 = K a t the translational part of the stiffness matrix of the initial simulation model.
Letting K t i = K m t α i | α i = 0 and Δ x t = x a t x m t , Equation (9) can be rewritten as:
i = 1 n α i K t i x m t = K a t Δ x t
Actually, Equation (10) is the model updating equation with respect to the element updated factors α i (i = 1, …, n). If the measured displacements or x m t contain the measurement errors, and these errors are considered as random quantities, the element updated factors will also be random. In this regard, Equation (10) is named as the statistical model updating equation in this paper.

2.3. Solution of Statistical Model Updating Equation

Assuming that the measurement errors of the beam structure under static load are random, the measured displacement vector of the actual structure can be expressed as:
x m t = x m t 0 + i = 1 l β i x m t i
where l is the number of random factors affecting the measurement, and equals to 1 given that the measurement errors are completely dependent. x m t 0 is the mean of the measured displacement vector, β i is the independent random variable, and x m t i is the adjoint vector of the random variable β i .
Using the multivariable power series to expand the element updated factors α i (i = 1, …, n), one can have
α i = α i 0 + j = 1 l α i j β j + j = 1 l k = 1 j α i j k β j β k + j = 1 l k = 1 j p = 1 k α i j k p β j β k β p + ...
where α i 0 , α i j , α i j k and α i j k p are respectively the zero-order to third-order of unknown expansion coefficients of the multivariable power series (12).
Substituting Equation (12) for Equation (10) leads to:
K a t ( x a t ( x m t 0 + i = 1 l β i x m t i ) ) = i = 1 n ( α i 0 + j = 1 l α i j β j + j = 1 l k = 1 j α i j k β j β k + ) K t i ( x m t 0 + i = 1 l β i x m t i )
Obviously, Equation (13) is a random algebraic equation with unknown expansion coefficients α i 0 , α i j , α i j k , etc. A high-order perturbation technique, which has been successfully used to determine the random frequencies of a large span cable-stayed bridge in the literature [27], is suggested to solve Equation (13). Using the high-order perturbation technique, the expansion coefficients α i 0 , α i j , α i j k and α i j k p , which correspond to the power polynomials 1, β j , β j β k and β j β k β p respectively, can be recursively determined. At first, for the zero-order polynomial term, there is:
α i 0 = A 0 1 ( K a t x a t K a t x m t 0 )
where A 0 = [ K t 1 x m t 0 , K t 2 x m t 0 , ... , K t n x m t 0 ] .
Next, for the first-order power polynomial β j , one has:
α i j = A 0 1 ( K a t x m t i + A 1 α i 0 )
where A 1 = [ K t 1 x m t j , K t 2 x m t j , ... , K t n x m t j ] .
Then, considering the second-order power polynomial, β j β k results in:
α i j k = A 0 1 A 2 α i j
where A 2 = [ K t 1 x m t k , K t 2 x m t k , ... , K t n x m t k ] .
For the third-order power polynomial β j β k β p , the expansion coefficient α i j k p is obtained as:
α i j k p = A 0 1 A 3 α i j k
where A 3 = [ K t 1 x m t p , K t 2 x m t p , ... , K t n x m t p ] .
In the same way, the fourth- to the higher-order of expansion coefficients can be determined. Taking into account the balance between the accuracy and the efficiency, up to the fourth-order of coefficients are used in this paper. Further, by the power polynomial expansion solved, the statistics of the element updated factors can be obtained easily. Note that many random methods, such as the Monte Carlo simulation methods, the spectral stochastic finite element methods (SSFEM), etc., can be used to solve the statistical model updating Equation (10). However, although SSFEMs have exhibited relatively high accuracy for many elastic mechanical problems [30], the calculation effectiveness or CPU time of the methods is not satisfactory for complex problems due to the exponential increase of the number of polynomial terms. Compared with these random methods, the proposed method is very efficient only with a little loss of the accuracy in some cases; it is therefore very suitable for solving large-scale engineering problems. For convenience, this paper denotes the proposed statistical model updating method (SMUM) as HP-SMUM, where HP means high-order perturbation.
Since the multivariate power series expansion has been obtained, Equation (1) can be rewritten as:
[ K a + i = 1 n ( α i 0 + j = 1 l α i j β j + j = 1 l k = 1 j α i j k β j β k + ) K i ] x a = F a
Using the high-order perturbation technique or other generally used methods, such as the spectral random finite element methods, the Monte Carlo simulation methods and so on, the statistical vertical displacements of the updated model can be obtained by Equation (18). Of course, facing the large-scale model updating problems, the high-order perturbation technique is still a very competitive solution method.

3. Numerical Examples

3.1. A Simply Supported Beam

Consider a simply supported beam with a rectangular section as shown in Figure 1. The length of the beam l = 6.0 m , the elastic modulus E = 2.8 × 10 10 Pa , the sectional area A = 0.3 m × 0.8 m , and the sectional inertia moment I = 0.0128 m 4 . The simply supported beam is divided into six elements with seven nodes, and each node includes the vertical deflection and rotation. A concentrated load P = 200 kN is applied in the middle of the beam. The rotational degrees of freedom at all nodes are eliminated using the static condensation method.
It is assumed that for the elastic modulus of the 4th element of the beam, a 10% reduction happens, which is due to the discreteness of the elastic modulus of concrete. According to the deterministic mechanical parameters given, the mean values of the assumed measurement deflections are simulated in this case. Considering the uncertainty of the measurement deflections, it is assumed that the coefficients of variation of the measured displacements at all nodes are 0.01 according to our experimental results which are 0.01~0.03, and the displacements are of Beta distributions.
By the simulated measurement data, the proposed HP-SMUM and the Monte Carlo simulation method with 100,000 samples are used to solve the statistical model updating equation to modify the original simulation model. The statistical results of the element updated factors are plotted in Figure 2. From Figure 2, it is found that the means and standard deviations of the element updated factors from the proposed method are in good agreement with those of the Monte Carlo simulation method. Furthermore, the mean value of the element updated factor of the 4th element indicates that the updated result consists of the assumed case where the 10% reduction of the elastic modulus of the 4th element is produced. Especially, from Figure 2c, it is observed that the probability density functions of the 4th element selected from the two calculation methods coincide with each other very well. These results illustrate that the proposed HP-SMUM is a very accurate solution for the statistical model updating equation. Meanwhile, by using the Monte Carlo simulation method to solve Equation (18), the means and standard deviations of the node deflections of the updated beam structure are obtained and shown in Figure 3. From Figure 3, it is found that the first two moments of deflections and the probability density function (pdf) of the deflection in the mid-span determined by the updated finite element model match with the assumed measurement values very well. For the mean of deflection, actual calculated relative errors are less than 1.38%. On the other hand, for the standard deviation of deflection, the relative errors at all nodes are close to zero. These results demonstrate that the proposed method has very good updated effectiveness.

3.2. A Two-Span Continuous Beam

Consider a two-span continuous beam as shown in Figure 4. The mechanical parameters of the continuous beam are the same as those of the simply supported beam in Section 3.1. The concentrated load P = 200 kN is applied in the middle of the left span of the continuous beam. Here, the 10% reduction of the elastic modulus of the 3rd and 4th element is assumed. Additionally, it is assumed that the coefficients of variation of the measured displacements at all nodes are 0.01, and they are of Beta distributions.
By the simulated measurement data, the proposed HP-SMUM and the Monte Carlo simulation method with 100,000 samples are used to obtain the statistics of the element updated factors, and the results are plotted in Figure 5. From Figure 5, it is found that the means and standard deviations of the element updated factors from the proposed HP-SMUM are in good agreement with those of the Monte Carlo simulation method. Furthermore, the means of the 3rd and 4th element shown in Figure 5a testify that a 10% reduction of the element stiffness has happened.
However, the probability density function of the 4th element selected from the proposed method is little different from the results of the Monte Carlo simulation method. This may be attributed to the non-symmetry of the continuous beam. Meanwhile, using the obtained element updated factors to solve Equation (10), the statistics of deflections of the updated continuous beam can be determined and shown in Figure 6 which shows that the first two moments of deflections by the updated finite element model match with the assumed measurement values very well. For the mean of deflection, actual calculated relative errors at all nodes are less than 0.54%. On the other hand, for the standard deviation of deflection, the relative errors are less than 7%. In addition, the pdf of the deflection in the mid-span by the updated model hardly deviates from the measured result.

4. Model Updating Test

A simply supported concrete beam under static load is shown in Figure 7a. The section of the reinforced concrete beam is 150 mm × 250 mm, and the length of the span is 1900 mm. For the distributive girder of static load, the length of the span is 950 mm. The grade of concrete strength is C25. This test was implemented according to a test standard as illustrated in literature [31].
The tested simply supported beam is divided into eight elements with nine nodes, as shown in Figure 7b. The concentrated load P is 15 kN. The deflections of seven nodes are measured several times so that the statistic values of the deflections are obtained. The statistics of the measured deflections at seven nodes are listed in Table 1 where the node deflections are assumed to be of Beta distributions. Based on the measured data, the proposed HP-SMUM is applied to update the concrete beam.
The means and standard deviations of the element updated factors are plotted in Figure 8 which shows that some degree of damage exists in the middle of the beam, and this result confirms the true case where small cracking happens in the concrete beam under the load.
Further, using the obtained element updated factors, the deflections of the concrete beam are calculated by means of Equation (18), and the results are shown in Figure 9. From Figure 9, it is found that the statistics of calculated deflections of the updated concrete beam coincide with the actual measured data very well. This result testifies the effectiveness of the proposed statistical model updating method once again.

5. Conclusions

A new statistical model updating method has been developed for modifying the beam structures with random parameters under static load. The new updating method considered the uncertainty of structural parameters and measurement errors. A static condensation technique is used to reduce the unmeasured degrees of freedom in the finite element model. A statistical model updating equation with respect to element updated factors has been established afterwards. Using the high-order perturbation technique, the statistical model updating equation has been solved to obtain the statistics of the element updated factors. The results of two numerical examples show that the accuracy of the proposed method is very good. The static responses obtained by the updated finite element model coincide with the measured results very well. Further, the updated results of the concrete beam based on static load tests testified the effectiveness of the proposed method.

Acknowledgments

The authors would like to acknowledge the financial support by the National Natural Science Foundation of China (Project No. 51378407 and 51578431).

Author Contributions

Bin Huang, Zhifeng Wu presented the model updating method; Zhifeng Wu, Yejun Li, Wuchuan Pu completed the numerical examples, designed and performed the experiments; Zhifeng Wu and Bin Huang wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Banan, M.; Hjelmstad, K. Parameter estimation of structures from static response. I. Computational aspects. J. Struct. Eng. 1994, 120, 3243–3258. [Google Scholar] [CrossRef]
  2. Banan, M.; Hjelmstad, K. Parameter estimation of structures from static response. II. Numerical simulation studies. J. Struct. Eng. 1994, 120, 3259–3283. [Google Scholar] [CrossRef]
  3. Simoen, E.; Roeck, G.; Lombaert, G. Dealing with uncertainty in model updating for damage assessment: A review. Mech. Syst. Signal Process 2015, 56–57, 123–149. [Google Scholar] [CrossRef]
  4. Erdogan, Y.; Gul, M.; Catbas, F.; Bakir, P. Investigation of uncertainty changes in model outputs for finite-element model updating using structural health monitoring data. J. Struct. Eng. 2014, 140, 1–14. [Google Scholar] [CrossRef]
  5. Sanayei, M.; Imbaro, G.; Mcclain, J.; Brown, L. Structural model updating using experimental static measurements. J. Struct. Eng. 1997, 123, 792–798. [Google Scholar] [CrossRef]
  6. Sanayei, M.; Onipede, O. Damage assessment of structures using static test data. Aiaa J. 1991, 29, 1174–1179. [Google Scholar] [CrossRef]
  7. Sanayei, M.; Scampoli, S. Structural element stiffness identification from static test data. J. Eng. Mech. 1991, 117, 1021–1036. [Google Scholar] [CrossRef]
  8. Liu, Y.; Ma, J.; Nie, J.; Zhang, S. Virtual distortion method-based finite element model updating of bridges by using static deformation. J. Eng. Mech. 2017, 143, 1–9. [Google Scholar] [CrossRef]
  9. Ren, W.; Fang, S.; Deng, M. Response surface-based finite-element-model updating using structural static responses. J. Eng. Mech. 2011, 137, 248–257. [Google Scholar] [CrossRef]
  10. Li, Y.; Zhang, W.; Liu, Y. A method based on meta-model for updating the finite element model of bridges using the measured static and dynamic data. In Proceedings of the Twelfth COTA International Conference of Transportation Professionals, Beijing, China, 3–6 August 2012; pp. 3175–3186. [Google Scholar]
  11. Sanayei, M.; Phelps, J.; Sipple, J.; Bell, E.; Brenner, B. Instrumentation, nondestructive testing and finite-element model updating for bridge evaluation using strain measurements. J. Bridge Eng. 2011, 17, 130–138. [Google Scholar] [CrossRef]
  12. Jacquelin, E.; Adhikari, S.; Friswell, M.I. A second-moment approach for direct probabilistic model updating in structural dynamics. Mech. Syst. Signal Process 2012, 29, 262–283. [Google Scholar] [CrossRef]
  13. Husain, N.; Khodaparast, H.; Ouyang, H. Parameter selection and stochastic model updating using perturbation methods with parameter weighting matrix assignment. Mech. Syst. Signal Process 2012, 32, 135–152. [Google Scholar] [CrossRef]
  14. Hua, X.; Ni, Y.; Chen, Z.; Ko, J. An improved perturbation method for stochastic finite element model updating. Int. J. Numer. Meth. Eng. 2008, 73, 1845–1864. [Google Scholar] [CrossRef]
  15. Govers, Y.; Haddad, H.; Link, M.; Mottershead, J. A comparison of two stochastic model updating methods using the DLR AIRMOD test structure. Mech. Syst. Signal Process 2015, 52–53, 105–114. [Google Scholar] [CrossRef]
  16. Khodaparast, H.; Mottershead, J.; Friswell, M. Perturbation methods for the estimation of parameter variability in stochastic model updating. Mech. Syst. Signal Process 2008, 22, 1751–1773. [Google Scholar] [CrossRef]
  17. Boulkaibet, I.; Mthembu, L.; Marwala, T.; Friswell, M.; Adhikari, S. Finite element model updating using the shadow hybrid Monte Carlo technique. Mech. Syst. Signal Process 2015, 52–53, 115–132. [Google Scholar] [CrossRef]
  18. Hua, X.; Ni, Y.; Chen, Z.; He, X. Monte Carlo study of the effect of measurement noise in model updating with regularization. J. Eng. Mech. 2015, 138, 71–81. [Google Scholar] [CrossRef]
  19. Wan, H.; Ren, W. Stochastic model updating utilizing Bayesian approach and Gaussian process model. Mech. Syst. Signal Process 2015, 70–71, 245–268. [Google Scholar] [CrossRef]
  20. Gomes, H.; Broggi, M.; Patelli, E.; Mottershead, J. Model updating by uncertain parameter inference. Vulnerability Uncertain. Risk 2014, 9, 1523–1532. [Google Scholar]
  21. Beck, J.; Katafygiotis, L. Updating models and their uncertainties. I: Bayesian statistical framework. J. Eng. Mech. 1998, 124, 455–461. [Google Scholar] [CrossRef]
  22. Katafygiotis, L.; Beck, J. Updating models and their uncertainties. II: Model identifiability. J. Eng. Mech. 1998, 124, 463–467. [Google Scholar] [CrossRef]
  23. Zárate, B.; Caicedo, J.; Wieger, G.; Marulanda, J. Bayesian finite element model updating using static and dynamic data. In Conference Society for Experimental Mechanics; Springer: New York, NY, USA, 2011; Volume 3, pp. 395–402. [Google Scholar]
  24. Cheung, S.; Bansal, S. A new Gibbs sampling based algorithm for Bayesian model updating with incomplete complex modal data. Mech. Syst. Signal Process 2017, 92, 156–172. [Google Scholar] [CrossRef]
  25. Goller, B.; Beck, J.; Schuëller, G. Evidence-based identification of weighting factors in Bayesian model updating using modal data. J. Eng. Mech. 2012, 138, 430–440. [Google Scholar] [CrossRef]
  26. Beck, J.; Au, S. Bayesian updating of structural models and reliability using Markov chain Monte Carlo Simulation. J. Eng. Mech. 2002, 128, 380–391. [Google Scholar] [CrossRef]
  27. Huang, B.; Seresh, R.; Zhu, L. Statistical analysis of dynamic characteristic of large span cable-stayed bridge based on high order perturbation stochastic FEM. Adv. Struct. Eng. 2013, 16, 1499–1512. [Google Scholar]
  28. Mottershead, J.; Friswell, M. Model Updating In Structural Dynamics: A Survey. J. Sound Vib. 1993, 167, 347–375. [Google Scholar] [CrossRef]
  29. Mottershead, J.; Link, M.; Friswell, M. The sensitivity method in finite element model updating: A tutorial. Mech. Syst. Signal Process 2011, 25, 2275–2296. [Google Scholar] [CrossRef]
  30. Ghanem, R.; Spanos, P. Stochastic Finite Elements: A Spectral Approach, 1st ed.; Springer: New York, NY, USA, 1991. [Google Scholar]
  31. Standard for Test Method of Concrete Structures (GB/T50152-2012); China Building Industry Press: Beijing, China, 2012.
Figure 1. A simply supported beam.
Figure 1. A simply supported beam.
Applsci 07 00601 g001
Figure 2. The statistical results of the element updated factors of the simply supported beam from the proposed high-order perturbation statistical model updating method (HP-SMUM) and the Monte Carlo simulation method (MC). (a) Means of element updated factors; (b) Standard deviations of element updated factors; (c) Probability density functions of the 4th element updated factor.
Figure 2. The statistical results of the element updated factors of the simply supported beam from the proposed high-order perturbation statistical model updating method (HP-SMUM) and the Monte Carlo simulation method (MC). (a) Means of element updated factors; (b) Standard deviations of element updated factors; (c) Probability density functions of the 4th element updated factor.
Applsci 07 00601 g002
Figure 3. The statistics of deflections of the simply supported beam based on the updated finite element model. (a) Means of deflections; (b) Standard deviations of deflections; (c) Probability density functions of the deflection at the 4th node.
Figure 3. The statistics of deflections of the simply supported beam based on the updated finite element model. (a) Means of deflections; (b) Standard deviations of deflections; (c) Probability density functions of the deflection at the 4th node.
Applsci 07 00601 g003
Figure 4. A two-span continuous beam.
Figure 4. A two-span continuous beam.
Applsci 07 00601 g004
Figure 5. The statistical results of the element updated factors of the continuous beam from the proposed HP-SMUM and the Monte Carlo simulation method (MC). (a) Means of element updated factors; (b) Standard deviations of element updated factors; (c) Probability density functions of the 3rd element updated factor; (d) Probability density functions of the 4th element updated factor.
Figure 5. The statistical results of the element updated factors of the continuous beam from the proposed HP-SMUM and the Monte Carlo simulation method (MC). (a) Means of element updated factors; (b) Standard deviations of element updated factors; (c) Probability density functions of the 3rd element updated factor; (d) Probability density functions of the 4th element updated factor.
Applsci 07 00601 g005
Figure 6. The statistics of deflections of the continuous beam based on the updated finite element model. (a) Means of deflections; (b) Standard deviations of deflections; (c) Probability density functions of the deflection at the 4th node.
Figure 6. The statistics of deflections of the continuous beam based on the updated finite element model. (a) Means of deflections; (b) Standard deviations of deflections; (c) Probability density functions of the deflection at the 4th node.
Applsci 07 00601 g006
Figure 7. Bending test of a simply supported concrete beam. (a) Bending test; (b) Mechanical model of the simply supported beam.
Figure 7. Bending test of a simply supported concrete beam. (a) Bending test; (b) Mechanical model of the simply supported beam.
Applsci 07 00601 g007
Figure 8. The statistical results of the element updated factors of the tested beam from the proposed HP-SMUM. (a) Means of element updated factors; (b) Standard deviations of element updated factors.
Figure 8. The statistical results of the element updated factors of the tested beam from the proposed HP-SMUM. (a) Means of element updated factors; (b) Standard deviations of element updated factors.
Applsci 07 00601 g008
Figure 9. The statistics of deflections of the tested beam based on the updated finite element model. (a) Means of deflections; (b) Standard deviations of deflections; (c) Probability density functions of the deflection at the 5th node.
Figure 9. The statistics of deflections of the tested beam based on the updated finite element model. (a) Means of deflections; (b) Standard deviations of deflections; (c) Probability density functions of the deflection at the 5th node.
Applsci 07 00601 g009
Table 1. Statistics of measured deflections of seven nodes.
Table 1. Statistics of measured deflections of seven nodes.
Node2345678
Mean (mm)0.1920.3570.4540.4880.4550.3580.196
Standard deviation (×10−2 mm)0.5761.0711.3621.4641.3651.0740.588

Share and Cite

MDPI and ACS Style

Wu, Z.; Huang, B.; Li, Y.; Pu, W. A Statistical Model Updating Method of Beam Structures with Random Parameters under Static Load. Appl. Sci. 2017, 7, 601. https://doi.org/10.3390/app7060601

AMA Style

Wu Z, Huang B, Li Y, Pu W. A Statistical Model Updating Method of Beam Structures with Random Parameters under Static Load. Applied Sciences. 2017; 7(6):601. https://doi.org/10.3390/app7060601

Chicago/Turabian Style

Wu, Zhifeng, Bin Huang, Yejun Li, and Wuchuan Pu. 2017. "A Statistical Model Updating Method of Beam Structures with Random Parameters under Static Load" Applied Sciences 7, no. 6: 601. https://doi.org/10.3390/app7060601

APA Style

Wu, Z., Huang, B., Li, Y., & Pu, W. (2017). A Statistical Model Updating Method of Beam Structures with Random Parameters under Static Load. Applied Sciences, 7(6), 601. https://doi.org/10.3390/app7060601

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop