Modeling and Stochastic Model Updating of Bolt-Jointed Structure

Bolt-jointed structure is widely used in engineering fields. +e dynamic characteristics of bolt-jointed structure are complex, and there is a variety of uncertainties in the jointed structure. In this study, modeling and updating of a typical bolt-jointed structure are investigated. In modeling terms, three-dimensional brick elements are used to represent the substructures, and thin-layer elements with virtual material properties are employed to represent the joint interface. Modal tests and experimental modal analysis of substructures and built-up structure are performed. A hierarchical model updating strategy based on Bayesian inference is applied to identify the unknown parameters in the substructures model and those in the overall model. Radial basis function (RBF)models are used as surrogates of time-consuming finite element model with high resolution to avoid the enormous computational cost. +e results indicate that the updated model can reproduce modal frequencies used in updating and can predict those not used in the updating procedure.


Introduction
Modeling and simulation (M&S) is extensively used in the engineering community to represent actual structures, and it is the third pillar of science [1].A reasonable finite element model (FEM) makes a variety of model-based activities, such as response analysis and prediction, structural health monitoring, and damage identification, possible.As a result, a precise FEM is of great significance to characterize the structural behavior.In the reality, engineering structures are composed of several relatively simple substructures that are connected together by bolt joints which results in the overall structure to be a complex system.e precision of outputs obtained from FEM is greatly dependent on the model form and model parameters what makes the finite element modeling become a challenging task.In the bolt-jointed structure, the modeling of substructures and joint interfaces are the two main aspects that affect the model performance.As a rule, substructures can be modeled by conventional elements, such as shell and brick elements, or in some other multiscale ways.However, joint interfaces are usually represented by some equivalent elements, including but not limited to node-to-node elements, zero-thickness elements, and thin-layer elements as well as refined elements with high fidelity.Bograd et al. gave a summary of different approaches for modeling the dynamics of mechanical joints in composition structures [2].Sandia National Laboratories released a handbook on the dynamics of jointed structures which includes modeling, experiments, and suggestions for future works [3].Mayer and Gaul used segment-to-segment contact elements to define a contact interface in an assembled structure [4].Luan et al. proposed a bilinear spring dynamical model and used it to analyze a pipe structure with bolted flange joints [5].Zhao et al. employed the fractal theory to obtain the stiffness and damping of the contact interface in a bolted assembly [6].Li et al. applied the six-parameter Iwan model to describe the nonlinear mechanisms of bolt joints [7].Wang et al. investigated the equivalent parameters of mechanical characteristics of bolt joints in a complex mechanical system by thin-layer elements [8].Xiang et al. proposed an improved spring method and used it in the modeling of the composite bolt-jointed structures [9].
However, the parameters in the incipient FEM built on the basis of the designed blue print are usually determined according to design handbooks or analyst's experience.Apparently, responses, calculated utilizing such a model, often differ from those measured from a real-world structure.Model updating is a powerful approach to minimize the difference between the model analysis results and the corresponding measurements.Mottershead and Friswell gave a summary of the model updating on the structural dynamics and published a treatise of this field [10,11].In the last two decades, a variety of works on modeling and updating of bolt joints have been put forward and implemented.Ahmadian et al. carried out model updating of joints in the AWE-MACE system utilizing a sensitivity method within the linear range [12].Zhai et al. constructed a parametric model of bolt joints of aeroengine stator system and updated the model with multicharacteristic responses [13].Cunha et al. applied the model updating method to identify the stiffness of bolt joints in civil construction structural systems [14].Adel et al. presented a simple and applicable model to predict the dynamic behavior of bolt joints and updated the model using modal test data [15].An in-depth understanding of model updating in structural dynamics can be acquired from Reference [16].
In the early stage, model updating was performed under the deterministic framework; that is, model parameters were estimated utilizing a single set of measurement.However, uncertainties are irremissible in both model form and model parameters.ese uncertainties are usually introduced by simplification of geometric features, assumption of constitutive models, and representation of joint interfaces.Ibrahima and Pettit presented an overview on the nonlinear characteristics and uncertainties in bolt joints pertaining to their dynamics [17].For the concern of multisource uncertainties in various aspects of structure dynamics, the stochastic model updating which refers to determine the range or probability density function of model parameters utilizing several experimental tests was put forward and promoted in recent years.A variety of methods were proposed for stochastic model updating, including the perturbation method, interval analysis, Monte Carlo simulation, fuzzy arithmetic, and Bayesian inference.Khodaparast et al. developed the model updating equations using the sensitivity method and estimated the statistical moments of updating parameters utilizing perturbation methods [18].Jiang et al. identified the statistical characteristics of model parameters using vibration test data and first-order perturbation method [19].Khodaparast et al. investigated the Kriging predictor metamodel and presented an interval model updating approach based on the Kriging predictor [20].Mares et al. explained the theory of stochastic model updating using the Monte Carlo simulation and applied the theory to a benchmark structure [21].Haag et al. presented a method for the validation of models based on the inverse fuzzy arithmetic [22].Wan and Ren.combined the Bayesian inference and Gaussian process model in the stochastic model updating and illustrated the proposed framework utilizing a simulated plate and a real-world bridge [23].Goller et al. elaborated the fundamental rule of Bayesian model updating and validated the method utilizing numerical examples [24].
However, in the most stochastic model updating activities, iterative computations or huge amount of sampling simulations are extremely time-consuming, especially for complex FEMs with high fidelity.For the sake of improving computational efficiency, surrogate models and component mode synthesis methods are investigated in the issues of structural dynamics.Rui et al. developed a stochastic model updating method by combining the stochastic response surface method with the Monte Carlo inverse error propagation [25].Fang et al. proposed a simple and cost-efficient method for stochastic model updating utilizing the response surface method and Monte Carlo simulation [26].Dohnal et al. addressed the joint uncertainties in a finite element model based on component mode synthesis and the stochastic reduced basis method [27].Yin et al. detected a defect in a structure with bolt joints using the model reduction method and system mode concept under the Bayesian framework [28].
e remainder of this manuscript is organized as follows: in Section 2, the basic theory of stochastic model updating based on the Bayesian inference, surrogate modeling and evaluation, basic theory of thin-layer elements are briefly introduced.Section 3 is devoted to the real bolt lapped structure, including finite element modeling, surrogate modeling, and model updating of substructures and the overall structure.e work ends with conclusions in Section 4.

Methodology
Modeling and simulation of bolt joints has been an active research subject in mechanical engineering and some other engineering fields.Model updating is a synthesis approach that combines numerical simulations and experimental tests.With the popularization of model updating methodologies and awareness of uncertainties, stochastic model updating has been widely used.Details about methodologies used in this article will be introduced in the following parts.

Stochastic Model Updating Based on Bayesian Inference.
Stochastic model updating based on the Bayesian inference is a process to evaluate the probability distribution functions (PDFs) by combining the existing knowledge about updating parameters with experimental observations [23].In essence, the Bayesian inference approach updates the prior PDFs of updating parameters utilizing experimental data and resulting in posterior PDFs.
Given that θ is the vector of updating parameters and D is the information matrix of experimental observations, the updated PDFs of updating parameters are governed by the Bayesian inference theorem: where y(θ) is model predictions, y is the mean vector, and σ is the covariance matrix of experimental measurements, respectively.Actually, models of engineering structures are usually complex and unlikely to be solved analytically.As a result, the obtainment of the posterior PDFs of updating parameters utilizing direct multidimensional integration is impractical.Markov chain Monte Carlo (MCMC) sampling is a promising approach to generate sampling points according to a specific rule for rejection, and it can avoid multidimensional integration in the Bayesian inference effectively.In essence, the MCMC sampling is a dynamical Monte Carlo simulation in which the distribution of sampling points varies with the process of simulation.e success of MCMC sampling is greatly dependent on the degree of proximity between the proposal PDFs of updating parameters with the target ones.
e Metropolis-Hastings (MH) sampling algorithm is one of the most widely used MCMC sampling methods.However, the efficiency of MH sampling algorithm is affected by model dimensions.For a model with low dimensions, a quick and perfect convergence can be achieved, but the result is not contrary for a higher dimensional model.For the problem of model updating presented in this article, a revised MH algorithm called as the Delayed Rejection Adaptive Metropolis-Hastings (DRAM) sampling algorithm is proposed [29].e DRAM is a combination of delayed rejection method and adaptive MH algorithm.
e main advantages of DRAM algorithm are as follows: (1) adaptive MH algorithm uses the ergodic information of a chain to adjust the proposal distribution, and it is a global strategy and (2) delayed rejection increases the efficiency of sampling based on the rejected proposal only within each time step, which means that it is a local strategy.

Surrogate Modeling and Evaluation.
e engineering community demonstrated a growing interest in applying surrogate models to replace a complex FEM for practical engineering problems.Surrogate modeling is a synthesis technology of experimental design and approximation method.Experimental design generates a set of design variables and calculates model responses utilizing the FEM.Approximation method builds the relationship between model responses and design variables.A surrogate model with high precision used in engineering problems can avoid a huge amount of computation effectively.e common experimental design methods are orthogonal design, uniform design, and Latin hypercube sampling (LHS).However, the selection of experimental design method should be cautious because a set of inappropriate sampling points can bring about low accuracy in the surrogate models and even to a worse result with the failure to construct these models [30].LHS is a stratified sampling process that can reflect the distribution of variables over the whole design space.Hence, the LHS method is employed in this paper to generate sampling points for the construction of surrogate models.
e approximation methods including but not limited to the polynomial response surface method (PRSM) [26], Kriging [20], support vector machine (SVM) [31], Gaussian process (GP) model [32], and radial basis function (RBF) model [26] are commonly used.A RBF model is constructed based on the basis function, which is symmetric and centered at each sampling point.A brief depiction of RBF model is given as follows: where f(x) is the approximation of objective function, N s is the total number of sampling points, λ i is the weight coefficient which is to be determined by the process of approximation, and ϕ(d i ) is the basis function of Euclidean distance d i .e Euclidean distance d i � ‖x − x i ‖ is a normal value, which represents the distance from the design variables x to the normalized design variables x i at the ith training point.Different forms of basis function used in constructing a RBF model are listed in Table 1.
e symbol c in the mathematical expression of different basis function, named as the smoothing coefficient, is a vital factor that determines the precision of RBF model.e precondition of using a surrogate model to replace the FEM is that the surrogate model shall have sufficient precision over the design space.
erefore, the accuracy evaluation of surrogate model is necessary after it is constructed.In addition, the accuracy evaluation of surrogate model at points used in constructing the model and those resampled from the design space should be conducted simultaneously.Two assessment criteria, known as coefficient of determination and root mean square error, are employed in this paper [33]: Shock and Vibration 3 where N t is the number of points used in the accuracy evaluation, Y j (x) is the jth value calculated from the surrogate model, Y Rj is the corresponding value obtained from FEM, and Y RJ is the mean value of sampling values.e closer the R 2 is to 1 and the closer the RMSE is to 0, the more accurate the surrogate model is.

Basic eory of in-Layer
Element.Modeling bolt joints utilizing thin-layer elements is an effective way to investigate dynamical characteristics of this kind of structures.e idea of thin-layer element was first proposed by Desai, and then thin-layer elements were applied to model the soil-structure interface as well as rock joints [4,34].e thin-layer element is an isoparametric brick element with very small but finite thickness.A brick element with length l, width w, and thickness t, as shown in Figure 1, is taken as the example to elaborate the formulation of thin-layer element.
e virtual work of a thin-layer element can be expressed as where V l is the volume of the element in the local coordinate system, σ is the element stress, ε is the element strain, u n is the nodal displacement, and K is the stiffness matrix of element formulated as follows where B is the transformation matrix and D is the constitutive matrix.e integration of formula (6), in the local coordinate system xyz, can hardly be solved.However, in the natural coordinate system ξηζ, in which the origin is located at the center of the cube element, coordinates of element nodes are transformed, and the stiffness matrix K can be calculated as where J is the Jacobian matrix that represents the relationship of partial derivatives between local and natural coordinates.
Apparently, in the natural coordinate system, numerical integration is simplified.As shown in Figure 1, the Jacobian matrix J is derived and simplified as follows: Substituting formula ( 8) into ( 7) and adopting Gaussian integral method, the stiffness matrix K can be numerically solved where w ξ m , w η n , and w ζ p are weight coefficients of numerical integration.
In the thin-layer element, the thickness of element t is very small and close to zero.According to Desai et al., a constitutive relationship with zero Poisson's ratio results in the decoupling between the normal and tangential behavior of the thin-layer element and in-plane strains and stresses become ignorable.at is, in the coordinate system xyz, if the z-direction is defined as the normal of the contact region, tangential strain components ε x , ε y , and c xy and tangential stress components σ x , σ y , and τ xy are all approximately equal to zero.As a result, the constitutive relationship of the thin-layer element is simplified as follows: where σ z is the normal stress in the z-direction and τ yz and τ zx are the shear stresses in the yz and zx directions, respectively.ε z is the normal strain in the z direction, and c yz and c zx are the shear strains in the yz and zx directions, respectively.E z is the normal stiffness, and G yz and G zx are tangential stiffnesses.Equation ( 10) is derived on the basis of linear elastic constitutive model.In order to describe the energy dissipation in the joint interface, damping coefficient should be introduced in the modeling phase.However, unlike the modeling of stiffness and mass, there is still a lack of suitable

Basis function Mathematical expression
Shock and Vibration damping model for the joint interface [35].In the engineering practice, constant hysteretic model, which assumes frequency-and amplitude-independent damping properties, is often used in the modeling of joint interface [36]. is model is optional in some commercial finite element codes.While the damping is introduced into the model, a complex stiffness matrix of thin-layer element that contains normal and tangential stiffness is obtained.e performances of thin-layer element are dependent on the virtual material properties and element thickness.To determine a reasonable thickness of thin-layer element, the max aspect ratio is employed defining as follows: where l is the length, w is the width, and t is the thickness of thin-layer element, respectively.However, up to now, no standard is available for selecting the maximum aspect ratio of thin-layer elements.Desai et al. suggested a value less than 100, while Pande and Sharma consider that a value larger than 1000 will eliminate numerical errors significantly [37,38].All these suggestions have their own special conditions.In a finite element analysis, the maximum aspect ratio should be selected according to the special condition and the software used [4].
In general, structures excited by small loading are considered to work in the closed state which means that there is no slippage between components of the overall structures, and the frictional force is large enough to contend with the applied forces.As a result, behaviors of mechanical joints are perceived as linear or can be linearized.

Experiments and Application
In assembled structures, uncertainties are inevitable in both substructures and bolt joints.As a typical type of assembled structure, bolt lapped plate shears, many common issues of assembled structures are widely used in various engineering fields.For the convenience and without loss of generality, in this paper, modeling and hierarchical stochastic updating of a bolt lapped plate are investigated.

Model Description.
e bolt lapped plate is composed of two plates that are assembled together by two bolts, as shown in Figure 2. e width and thickness of the plates are both 120 mm and 2.76 mm, and their lengths are 450 mm and 550 mm, respectively.

Modeling of Bolt Lapped Structure.
In order to separate the uncertainties and to make the modeling convenient, the assembled bolt lapped plate is divided into two substructures, named as plate 1, plate 2, and two bolt joints.e two plates are molded utilizing isoparametric hexahedral elements, and the bolt joints are modeled by thin-layer elements.in-layer elements are connected with elements of plates by node superposition in the contact region.A lumped point element is created in the location of sensor to consider the weight of sensor and its base.e overall FEM of bolt lapped plate is shown as Figure 3.
e performances of structures are determined by material and connection properties.In the bolt lapped structure, isotropic material properties with nominal values of steel are assigned to substructures and orthogonal anisotropy material properties with virtual parameters are assigned to thin-layer elements.As shown in Figure 2, the length of the contact zone is 40 mm and the diameter of bolts is 10 mm. e stiffness and damping in the contact zone differ from those in the rest zone of the plates.In order to obtain precise models of substructures, as shown in Figure 3, material properties of elements in the contact zone and those in the rest region on the corresponding plates are defined separately.In addition, the acting force of bolt connection converges in the area close to the bolt and much smaller in the rest area of the contact region.Hence, thin-layer elements in the FEM are divided into three parts, two of them include elements in the range of two radiuses of each bolts, denoted as Bolt 1 and Bolt 2, and another part contains the rest thin-layer elements in the contact region, denoted as Contact.To account for the difference of the connection states, all of the three parts are assigned with distinct virtual material properties.e initial material properties of substructures and thin-layer elements are listed in Tables 2 and  3, respectively.
Apparently, inherent attributes and structural dynamic responses calculated by the initial FEM built based on nominal values or analyst's experience may differ from those of real structures, more or less.For the further study based on FEM, it is necessary to update variable parameters in the model to narrow the gap between simulations and experimental tests.

Modal Tests and Data Analysis.
A hierarchical test scheme is adopted in the modal tests; that is, substructures are tested firstly and then the overall plates are tested.e single-input single-output (SISO) scheme is adopted during the test process.e test plane is designed according to the Tutorial Guideline VDI 3830: Damping of Materials and Members [39].All test specimens are hanged horizontally by two rubber ropes on a rigid bracket to simulate the free-free boundary condition.Locations of accelerometers are selected elaborately to ensure that sufficient orders of modes can be identified from tests.Moreover, uncertainties of modal characteristics are unavoidable in substructures and assembled structures.In order to consider the uncertainties of modal characteristics, modal tests of specimens with the same sizes are conducted, and accelerometers are kept at their locations without movement during the test to eliminate the influences of sensors.
In the first stage, modal tests of substructures are conducted.Uncertainties of substructures are introduced by the variance of specimens with the same nominal sizes and material properties, but they are always taken from different production batches.In this paper, ten specimens of each kind of substructures are tested.For the case to eliminate the system error, every specimen is tested for three times in the Shock and Vibration same condition, and the results are averaged as the final modal parameters.
e experimental modal frequencies, modal damping and their statistics characteristics of substructures, are shown in Tables 4 and 5, respectively.
In the second stage, modal tests of overall structure are conducted.
e bolt lapped structure is assembled with plates that are randomly selected from the set of substructures tested in the first stage, and the torques used to connect substructures vary from 10 Nm to 30 Nm.In this stage, the bolt lapped structures are assembled for fifty times, and modal tests are conducted on every corresponding specimen.As is same to the substructures, every assembly of bolt lapped structure is tested for three times, and the test results are averaged as the final modal parameters.e test site of the overall structure is shown as Figure 4, and the identified modal frequencies and their statistic characteristics are shown in Table 6.

Surrogate Modeling of Structures.
In the model updating process, the iterative analysis of FEM with different groups of model parameters is usually timeconsuming and sometimes even impossible to be completed.Surrogate model is a promising alternative of FEM that represents the relationship between model outputs and design variables.In substructures, the elastic modulus, densities, and damping coefficients are selected as design variables, and the six modal frequencies and damping shown in Tables 4 and 5 are selected as objectives for the construction of RBF models.e elastic modulus and density parameters vary from 0.95 to 1.05 times of the corresponding values shown in Table 2, and the damping  e RBF model of first bending mode of plate 1 is shown as Figure 5.
In the overall structure, normal stiffness, tangential stiffness, and damping coefficients are selected as design variables, and six modal frequencies and the corresponding modal damping shown in Table 6 are selected as objectives of the construction of surrogate models.e design variables vary from 0.1 to 100 times of the values shown in Table 3. e assessment criteria, R 2 , vary from 0.9904 to 0.9996, and   Shock and Vibration RMSE varies from 1.53 × 10 −4 to 1.80 × 10 −3 , which indicates that the constructed RBF models of the bolt lapped plate are also fairly accurate.

Results of Model Updating.
Model updating is completed by revising model parameters to minimize the gap between FEM analysis results and experimental test results.Given the omnipresent uncertainties, all of the model updating tasks are completed under the uncertainty framework hierarchically.Details of updating schemes and results are presented as follows.
In the first stage, elastic modulus, densities, and damping coefficients of substructures are selected as updating parameters.During the process of stochastic model updating, the parameters are assumed to obey the prior multidimensional normal distribution and parameter boundaries are the same as those used in constructing surrogate models.Totally, 10000 sampling points are generated using the DRAM algorithm.Based on the experimental tests and constructed RBF models, updating parameters are calibrated.Statistical characteristic values, that is, mean values (Mean) and standard variances (SD) of updating parameters after updating are shown in Table 7.
After completing the process of stochastic model updating of substructures, a set of 10000 sampling points in the 99% confidence interval of updating parameters is generated to validate the accuracy of updated results.e statistical characteristics of experimental and updated modal frequencies of plate 1 and plate 2 are listed in Tables 4 and 5, respectively.
As shown in Table 4, the maximum error of modal frequency mean value is 0.26%, and that of standard variance value is 18.09%.e absolute mean error of modal frequency mean values is 0.16%, and that of modal frequency standard variances is 11.84%.In addition, the maximum error of modal damping mean value is 6.34%, and that of standard variance value is 22.03%.e absolute mean error of modal damping mean values is 2.83%, and that of modal frequency standard variances is 14.60%.
Meanwhile, as shown in Table 5, the maximum error of modal frequency mean value is −1.13%, and that of standard variance value is 18.82%.e absolute mean error of modal frequency mean values is 0.36%, and that of modal frequency standard variances is 10.47%.In addition, the maximum error of modal damping mean value is −5.91%, and that of standard variance value is −26.09%.e absolute mean error of modal damping mean values is 3.27%, and that of modal frequency standard variances is 19.47%.
In the second stage, stochastic model updating of the overall structure based on Bayesian inference is conducted.After the work done in the first stage, the statistic characteristics of material properties in the substructures are determined.During the process of updating for the overall structure, updated parameters of substructures are substituted into the model of overall structure and virtual material properties of thin-layer elements are assigned as the updating parameters.As mentioned above, thin-layer elements are divided into three parts in the modeling procedure, that is, Bolt 1, Bolt 2, and Contact region.All of these parts are mutually independent, and they are assigned with individual parameters.During the updating process, one normal stiffness, two tangential stiffness, and damping coefficient of every part thin-layer elements are assigned as updating parameters; that is, totally twelve parameters are allowed to vary until the stable convergence.For the convenience of listing the results of model updating, normal stiffness of Bolt 1, Bolt 2, and Contact region is assigned as E 1 , E 2 and E 3 , tangential stiffness is assigned as G yz1 , G zx1 , G yz2 , G zx2 , G yz3 , and G zx3 , and the corresponding damping coefficients are assigned as η 1 , η 2 , and η 3 , respectively.During the process of stochastic model updating, parameters are assumed to obey the prior multidimensional normal distribution, and variable spaces of updating parameters are the same as those used in constructing surrogate models.e advanced MCMC method DRAM algorithm with 10000 evaluations is adopted to determine the posterior PDFs of updating parameters.Statistic characteristics of updating parameters are listed in Table 8, and the comparison between experimental and updated results is shown in Table 6.Shock and Vibration As shown in Table 6, the maximum error of modal frequency mean value is 1.75%, and that of standard variance value is 35.99%.e absolute mean error of modal frequency mean value is 0.71%, and that of standard variance value is 18.83%.However, when it comes to damping, the maximum error of modal frequency mean value is 35.47%, and that of standard variance value is −29.93%.e absolute mean error of modal frequency mean value is 14.29%, and that of standard variance value is 20.20%.e results shown in  indicate that favorable agreements between experimental results and updated results of modal frequencies and modal damping are obtained by stochastic model updating based on the Bayesian inference.
Modal frequencies and modal damping calculated by sampling points generated according to the initial parameters and updated parameters, as well as the experimental results are shown in Figure 6.
As shown in Figure 6, the modal frequencies and damping between the results calculated by the initial parameters and the updated results are obvious.However, the updated results are closer to the experimental ones and the differences.
e purpose of model updating is to obtain a FEM with a high precision which can predict some structural characteristics or responses that are not used in the updating process.Modal frequencies and modal damping of third torsion and sixth bending are used to validate the prediction capability of the updated FEM.Mean value and standard variance of modal frequencies and modal damping used in validation are listed in Table 9.
As shown in Table 9, the statistics characteristics of modal frequency used to validate the precision of FEM are very close to the experimental ones.e maximum error of modal frequency mean value is 1.01%, and that of standard variance value is −10.71%.e absolute mean error of modal frequency mean value is 0.83%, and that of standard variance value is 5.85%.Meanwhile, the maximum error of modal damping mean value is −22.77%, and that of standard variance value is −26.91%.e absolute mean error of modal damping mean value is 16.38%, and that of standard variance value is 21.29%.
Moreover, a rough estimation of the time consumed in stochastic model updating of the overall bolt lapped structure with 10000 iterations based on FEM directly is about 4.6 days that means that it has a huge amount of computation.However, time consumed in the process of stochastic modal updating with the same times of iterations based on RBF models of the overall bolt lapped structure is about 4 minutes.e comparison of consumed time used in the model updating indicates that model updating based on surrogate models is a desirable way to save the computational cost.

Conclusion
In this paper, the issue of hierarchical finite element modeling and stochastic model updating of a bolt lapped plate is presented.e whole procedure is implemented in a progressive manner.In the part of modeling, isoparametric hexahedron elements and thin-layer elements with virtual material properties are adopted to represent substructures and bolt joint interfaces, respectively.in-layer elements make it possible to assign different material properties that represent the variance of interface stiffness and damping.In the part of stochastic model updating, material properties of substructures are calibrated firstly, and then virtual material properties that represent the interface stiffness and damping

Figure 3 :Figure 2 :
Figure 3: (a) FEM of the overall bolt lapped structure.(b) Partial enlarged detail of contact region.(c) Locations of bolts.(d) Partition of material properties on substructures.

Figure 4 :
Figure 4: Test site of overall plate.

Figure 5 :
Figure 5: RBF model of first bending mode of plate 1.

Figure 6 :
Figure 6: Comparison of modal frequencies and modal damping.

Table 1 :
Different forms of basis function for RBF model.

Table 2 :
Material properties of plates.

Table 3 :
Material properties of thin-layer elements.
parameters vary from 0.1 to 20 times of the corresponding values shown in Table2.e assessment criteria, R 2 , are all close to 1 and RMSE varies from 2.26 × 10 −6 to 2.20 × 10 −5 .evalues of assessment criteria indicate that the constructed RBF models of substructures are highly accurate.

Table 4 :
Experimental and updated results of plate 1.

Table 5 :
Experimental and updated results of plate 2.

Table 6 :
Experimental and updated results of overall structure.

Table 7 :
Updated parameters of substructures.

Table 8 :
Updated parameters of overall substructure.However, stochastic model updating based on the Bayesian inference requires a large amount of sampling points, and the direct model analysis at these sampling points is time-consuming.To solve this problem, surrogate models named as RBF models and an advanced MCMC algorithm named as DRAM algorithm are employed in this paper.e modal frequencies and modal damping calculated from the updated model show a great agreement with the corresponding experimental results.As a consequence, the strategy combining surrogate models and Bayesian inference in the stochastic model updating is a commendable way to improve the precision of FEM and enhancing computational efficiency in engineering practices.

Table 9 :
Results of validating modal frequencies.