Abstract
The present work presents a model for nonlocal and viscoelastic Euler-Bernoulli beams and aspects of its calibration are addressed. The nonlocal feature of the model is described by the nonlocal elasticity theory proposed by Eringen and its viscoelastic behaviour is modelled by means of internal variables. Parametric analyses are performed to determine the impact of the nonlocal and viscoelastic parameters on the modal properties of the system. Inverse analyses are performed under the Bayesian framework and samples of the posterior density function are obtained by means of the Delayed Rejection Adaptive Metropolis (DRAM). The inverse analyses consider a nonlocal viscoelastic beam model with one internal variable and they address three aspects, namely: the impact of a misspecification of the beam diameter, the impact of modelling the beam diameter as an unknown but uninteresting model parameter and the model calibration when synthetic experimental data comes from a model containing two internal variables. The model parameters were chosen such that the system resembles a Single-Walled Carbon Nanotube (SWCNT).
1. Introduction
Due to the superiority of their mechanical, electrical, thermal, chemical and optical properties, the carbon-based nanomaterials such as the Carbon Nanotubes (CNTs) and Graphene Sheets (GSs) are of great interest in science and technology and they find applications in nanoelectronics [Citation1,Citation2], nanosensors [Citation3,Citation4] and nanostructured composites [Citation5,Citation6]. However, since experimental analyses at the nanometre scale are extremely difficult to carry out, mathematical models for nanocomponents play an important role on the analysis and design of nanoelectromechanical systems.
Two approaches may be used to build models for nanostructures: atomistic modelling and continuum mechanics based modelling [Citation7–9]. In the atomistic approach the structure is modelled considering the individual motion of each atom [Citation10] and it has been reported to provide accurate predictions. On the other hand, the atomistic approach tends to be computationally prohibitive. Regarding the continuum mechanics approach, it is relatively simple and less computationally intense than the atomistic one. On the other hand, the classical (local) continuum elasticity theory fails on accurately modelling nanoscale structures due to the size dependence effects, which demands the use of nonlocal elasticity theories such as the one proposed by Eringen [Citation11,Citation12]. Mechanical models based on the nonlocal elasticity theory of Eringen have been proposed in the literature for modelling nanomaterials and nanoelectromechanical systems (NEMS) [Citation13–23]. Inherent to these models is the nonlocal parameter , where is a material parameter and a is an internal characteristic length of the system, usually adopted as the length of the carbon-carbon bond. Although the nonlocal parameter is of fundamental importance in the theory of nonlocal elasticity, there is no consensus on its determination. As the experimental analyses at the nanometre scale are extremely difficult to carry out, most of these models are calibrated based on numerical results provided by atomistic approaches, specially from molecular dynamics simulations. Eringen [Citation12] proposed for the material parameter the value by matching the dispersion curves of plane waves with those of Born-Karman model of atomic lattice dynamics and the value when considering the Rayleigh surface waves. Regarding the task of estimating the parameter , one finds works based on nonlocal Euler-Bernoulli beam model [Citation24], nonlocal Timoshenko beam model [Citation25] and nonlocal shell models [Citation26–28]. The task of simultaneously estimating the nonlocal parameter and Young's modulus of one dimensional atomic chains and armchair Single-Walled Carbon Nanotube SWCNTs is approached in [Citation29].
Since the dynamic behaviour of a nanostructure may be significantly affected by damping or viscoelastic effects, an accurate modelling of these dissipative mechanisms is of prime importance [Citation30–33]. At the nanoscale, the damping may be generated by external forces, such as magnetic forces [Citation34], the interaction between the nanostructure and the substrate, humidity or thermal effects [Citation35]. Regarding the description of dissipative mechanisms, generalized viscous and viscoelastic damping models are found in the literature within the context of nonlocal elasticity theory [Citation36–39]. Nonlocality and viscoelastic models are approached, for instance, in the works by Lei [Citation36,Citation37], Ansari et al. [Citation38] and Cajić et al. [Citation39]. Nonlocal viscoelastic beam models embedded in viscoelastic mediums are addressed by Zhang et al. [Citation40–42]. The aforementioned works on nonlocal viscoelastic damping models provide valuable insights about the dynamic behaviour of nanomaterials through parametric studies of the effects of nonlocal, viscoelastic and other parameters on their modal properties or time domain responses. However, it should be highlighted that none of them addresses the problem of estimating the nonlocal and/or viscoelastic parameters.
The present work is concerned with modelling nonlocal viscoelastic behaviour in Euler-Bernoulli beams and presenting aspects of the corresponding model calibration under the formalism of Bayesian inference [Citation43]. The nonlocal model is described by the nonlocal elasticity theory proposed by Eringen [Citation12] combined with the theory of internal variables for viscoelasticity [Citation44–46]. The posterior probability density function (PDF) of the model parameters is sampled using the Delayed Rejection Adaptive Metropolis (DRAM) algorithm [Citation47]. The inverse analyses consider a nonlocal viscoelastic beam model with one internal variable and they address three aspects, namely: the impact of a misspecification of the beam diameter, the impact of modelling the beam diameter as an unknown but uninteresting model parameter and the model calibration when synthetic experimental data comes from a model containing two internal variables. From the authors best knowledge, there exists no work in the literature that addresses the simultaneous estimation of the nonlocal parameter and viscoelastic parameters using the Bayesian framework.
2. Mathematical modelling
2.1. Nonlocal elasticity theory
According to the nonlocal elasticity theory of Eringen [Citation12], the stress at a given reference point in the body depends on the strain field at every point in the elastic domain. The nonlocal stress tensor is given by (1) (1) where is the classical local stress tensor, the kernel function is the nonlocal modulus (attenuation function), is the Euclidean distance and ϕ is a material constant that depends on internal and external characteristic lengths, such as lattice distance and wave length, respectively. Usually, one has , where is a parameter used for calibrating the nonlocal model by means of data obtained from experiments, a represents an internal characteristic length and l an external one.
The nonlocal kernel proposed by Eringen [Citation12] is given by (2) (2) where is the modified Bessel function.
In order to obtain a differential form for the nonlocal constitutive relation in Equation (Equation1(1) (1) ), the nonlocal kernel is considered as a Green's function of a linear differential operator as follows [Citation12] (3) (3) where is the linear differential operator and δ is the Dirac delta distribution. Applying the linear differential operator in the nonlocal constitutive relation in Equation (Equation1(1) (1) ) and considering Equation (Equation3(3) (3) ), one has (4) (4) Considering 1D apllications, one may combine Equation (Equation4(4) (4) ) and the nonlocal modulus in Equation (Equation2(2) (2) ) to obtain the nonlocal constitutive relation in differential form as follows (5) (5) where is the Laplacian operator.
2.2. Nonlocal viscoelastic model based on internal variables
A nonlocal viscoelastic constitutive relation can be obtained by combining the nonlocal elasticity theory with a local constitutive viscoelastic relation [Citation36]. For the one-dimensional case, the integral form of the local constitutive relation of a linear and homogeneous viscoelastic material is given by [Citation48] (6) (6) where is the extensional relaxation kernel of the viscoelastic material, τ and ε are, respectively, the longitudinal local stress and strain [Citation49]. A one-dimensional nonlocal viscoelastic constitutive relation is provided by Equations (Equation5(5) (5) ) and (Equation6(6) (6) ) as follows (7) (7) In the present work, the relaxation kernel is adopted as follows (8) (8) where is the relaxed or static modulus, N is the number of internal variables which in turn are in charge of describing viscoelastic dissipation mechanisms and and correspond to the the relaxation amplitude/modulus and to the inverse of relaxation time associated to the ith internal variable as shown next.
The combination of the Laplace transform of Equations (Equation7(7) (7) ) and (Equation8(8) (8) ) leads to (9) (9) where denotes the Laplace transform of the variable . The ith internal variable is defined as follows (10) (10) Finally, combining Equations (Equation9(9) (9) ) and (Equation10(10) (10) ) and applying the inverse Laplace transform leads to the differential form of the constitutive relation of a nonlocal viscoelastic material (11) (11) where the second equation in (Equation11(11) (11) ) is the evolution equation of the ith internal variable.
2.3. Nonlocal viscoelastic Euler-Bernoulli beam model
In the Euler-Bernoulli beam model [Citation50], the resultant bending moment is given by (12) (12) where z is the coordinate along the thickness of the beam, such that the neutral plane is given by z = 0, is the axial stress and A is the cross sectional area.
The governing equation of the Euler-Bernoulli beam model casts as follows (13) (13) where ρ is the mass density, w is the transverse deflection of the midline of the beam, C is the viscous damping coefficient and f is the distributed transverse load per unit length.
Applying the nonlocal operator in Equation (Equation12(12) (12) ) and considering the first equation in Equation (Equation11(11) (11) ) leads to (14) (14) Considering the kinematics description assumed for Euler-Bernoulli beams and the hypotheses of small strains and small displacements, the axial strain is given as (15) (15) As for the internal variables , inspired by Equation (Equation10(10) (10) ), one may consider its definition linked to an auxiliary internal field as follows (16) (16) where the governing equation for is given as (17) (17) Substituting Equations (Equation15(15) (15) ) and (Equation16(16) (16) ) into Equation (Equation14(14) (14) ) leads to (18) (18) where I is the cross-sectional area moment of inertia about the y axis.
Considering Equations (Equation13(13) (13) ), (Equation17(17) (17) ) and (Equation18(18) (18) ), the differential equations that govern the dynamics of the nonlocal viscoelastic Euler-Bernoulli beam cast as follows (19) (19) Equation (Equation19(19) (19) ) is a generalization of the Euler-Bernoulli beam model. For instance, the equation of motion of a local and undamped Euler-Bernoulli beam may be obtained by taking , C = 0 and and the equation of an undamped nonlocal Euler-Bernoulli beam is obtained by taking C = 0 and , .
2.4. Finite element model of a nonlocal and viscoelastic Euler-Bernoulli beam
In this work, the finite element model of the nonlocal viscoelastic Euler-Bernoulli beam is derived by means of the weighted residuals approach [Citation19,Citation51]. Starting from Equation (Equation19(19) (19) ), the weighted-integral equations within a finite element of the beam may be written as (20) (20) where denotes the length of the element, represents the local coordinate and is the Galerkin's weight function [Citation51].
Considering Hermite interpolation functions for the transverse displacement and for the corresponding auxiliary internal field, within a finite element of the beam, one has (21) (21) where is the matrix comprised of the interpolation functions, is the vector of generalized nodal displacements of the element and is the vector describing the variables of ith internal variable.
The governing equations of a beam element is given as follows (22) (22) where , and are, respectively, the mass, viscous damping and relaxed stiffness matrices of a finite element and is the vector of generalized nodal forces. The assembly of element matrices allows one to build the finite element model of the nonlocal viscoelastic beam. Finally, the vectors and matrices in Equation (Equation22(22) (22) ) are defined in Appendix 1.
3. Inverse problem
In the Bayesian framework for inverse problems, model parameters are modelled as random variables, uncertainties associated to a random variable are described by means of its probability density function (PDF) and the problem solution becomes the determination of the posterior PDF of model parameters given a set of measured data [Citation43], namely .
Assuming an additive error observation model one may write (23) (23) where is a random vector in charge of representing measurement errors and is an operator used to provide model predictions. As , and are random variables, one may obtain a relation among them by means of the Bayes rule as follows (24) (24) where is the posterior PDF of model parameters, is the likelihood function and is the prior density, which in turn represents the information one has about the parameters before considering any information about the experimental data . Concerning , it represents the distribution of measured data and is given by (25) (25) Assuming that follows a zero mean Gaussian density with covariance matrix formally represented by , and that and are independent random variables, the likelihood function is given (26) (26) Equation (Equation24(24) (24) ) presents the full probabilistic model for obtaining information about model parameters . Nevertheless, Equation (Equation24(24) (24) ) is analytically intractable in general. Therefore, in order to obtain information about the posterior one may resort to sampling based techniques, such as Markov Chain Monte Carlo (MCMC) methods. The Delayed Rejection Adaptive Metropolis algorithm is used to simulate the posterior density in the present work.
3.1. Delayed rejection adaptive metropolis (DRAM)
The DRAM algorithm, proposed by Haario et al. [Citation47], with the aim of improving the efficiency of the Metropolis-Hastings (MH) algorithm, combines two strategies: the Adaptive Metropolis (AM) [Citation52] and the Delayed Rejection (DR) [Citation53]. Assuming the target density that one wishes to simulate is the posterior density , the algorithm is briefly described next.
In the standard MH algorithm, given a current state of the Markov Chain, a candidate is generated, from an auxiliary PDF , and it is accepted with probability given as follows (27) (27) Upon acceptance, the candidate state is the new state of the chain, . Otherwise, the current state is retained, .
In the DR algorithm, upon rejection of the candidate , instead of retaining the current state, a second candidate is proposed from the auxiliary PDF which explicitly depends on the current state and the previously proposed and rejected candidate . The second candidate is accepted with probability given as follows (28) (28) The process of generating new candidates, built on the current state and the previously proposed and rejected candidates, within a step of the Markov Chain, may be performed for a fixed or random number of stages. At the ith stage, the candidate is accepted with probability (29) (29) It should be emphasized that, since the acceptance probability of each stage of the DR algorithm is defined so that the reversibility with respect to the target density is preserved, the delaying rejection process can be interrupted at any stage [Citation47].
In the present work, the auxiliary PDF of the first stage of the DR is derived according to the AM strategy. More specifically, considering the nth state of the chain, the auxiliary density of the first stage of the DR is a multivariate Gaussian centred at the current state with covariance matrix , whose adaptation is based on the history of the chain, as shown in what follows. In higher stages, the auxiliary PDFs are also Gaussians whose covariance matrices are simply scaled versions of the covariance at the first stage.
According to the AM strategy, the auxiliary Gaussian density function is centred at the current state of the chain and its covariance matrix is given by (30) (30) where is the length of the initial non-adaptation period, is an arbitrary strictly positive definite matrix, representing the covariance matrix of the auxiliary density during the non-adaptation period, is a parameter that depends only on the dimension of the domain of the posterior, is the covariance matrix computed from the previous states of the chain, λ is an arbitrarily small positive constant and is the dimensional identity matrix.
Considering a sequence of states , the empirical covariance matrix may be written as (31) (31) where represents the mean value, given by (32) (32) Considering Equations (Equation30(30) (30) ) and (Equation31(31) (31) ), for , one has the recursive formula for the covariance matrix given by Equation (Equation33(33) (33) ) that yields a reduction in the computational burden of the AM algorithm [Citation47]. (33) (33) Therefore, the DRAM algorithm is such that, at the first stage of the DR strategy, the auxiliary PDF is Gaussian, centred at the current state of the Markov chain and with covariance matrix , i.e. . At higher stages, the auxiliary density is where is a shrinkage factor.
4. Numerical results
This section presents some forward and inverse analyses with the nonlocal viscoelastic beam model described by Equation (Equation19(19) (19) ).
The geometrical and physical characteristics of the system are chosen such that they resemble the ones of a Single-Walled Carbon Nanotube (SWCNT) as shown in Table . Further, it should be remarked that, in the analyses that follows, the parameters adopted are in agreement with the ones presented in the literature [Citation14,Citation36].
The numerical analyses are made considering the FE model of the beam which was spatially discretized into 20 finite elements. As for the internal variables used to describe dissipative behaviour, all the analyses considered a single internal variable. Therefore, for the sake of simplicity, from now on, the relaxation modulus and the inverse of relaxation time shown in Equation (Equation8(8) (8) ) will be simply termed Δ and Ω, respectively.
Finally, it should be highlighted that the natural frequencies of nanobeams are within the high-frequency range. As a basis of reference for the next analyses, a classical (local) and elastic Euler-Bernoulli beam model was initially built considering the properties shown in Table and assuming and in Equation (Equation22(22) (22) ). Table presents the first three undamped natural frequencies of this local and elastic beam.
The effects of the nonlocality and viscoelasticity on modal properties of the system will be shown in Subsection 4.1.
4.1. Nonlocal and viscoelastic effects on the modal parameters of the system
The effects of the nonlocal and viscoelastic parameters on the first three undamped natural frequencies and damping ratios of the system are addressed in this subsection. The analysis is carried out considering the basic parameters in Table along with the nonlocal and viscoelastic parameters presented in Table .
The first three undamped natural frequencies and damping ratios of the nonlocal viscoelastic Euler-Bernoulli beam model of the SWCNT are computed from the poles of the system. The modal parameters are computed from the rth complex pole of the system by means of the relations: (34) (34) (35) (35) where is the notation used for the complex conjugate of a variable q. Appendix 2 shows how the poles of the system are obtained from the finite element model matrices.
Figures present the undamped natural frequencies and damping ratios considering the geometrical and physical parameters shown in Table along with the nonlocal and viscoelastic parameter ranges presented in Table . In these figures, the nonlocal undamped natural frequencies () are normalized by the undamped natural frequencies () shown in Table , which are provided by the corresponding local and elastic beam model for r = 1, 2, 3.
From Figures to , one may observe that the nonlocal parameter significantly affects the undamped natural frequencies of the system. An increase in the nonlocal parameter results in a decrease in the natural frequencies. Besides, the higher the mode, the greater the effect of on the undamped natural frequency, showing that higher frequencies are more sensitive to small scale effects than the fundamental frequency. On the other hand, the effect of the nonlocal parameter on the damping ratios of the first three vibration modes is rather small. The damping ratios of the first two vibration modes remain practically constant with respect to .
Regarding the dependence of the modal properties with the viscoelastic parameter Ω, from Figures to , one may observe that the undamped natural frequencies present a transition region when the viscoelastic parameter Ω is near the value of the corresponding natural frequency of a given mode [Citation36]. Within the transition region, a decrease in the undamped natural frequencies is observed with an increase in the viscoelastic parameter Ω. Outside the transition region, the undamped natural frequencies present extremely low sensitivity with respect to the viscoelastic parameter Ω. On the other hand, Ω has a very significant effect on the damping ratios. One may observe that the damping ratios present a very pronounced ridge when Ω is in the vicinity of the corresponding natural frequency. Besides, the damping ratios monotonically decrease when Ω deviates from the value of the corresponding natural frequency. Small damping ratios are observed when Ω is far from the pronounced ridges.
4.2. Bayesian parameter estimation
This subsection is devoted to present some numerical analyses concerning the calibration of the nonlocal viscoelastic Euler-Bernoulli beam model of a simply supported SWCNT. Here, a model containing N = 1 internal variable is used to provide synthetic data for the inverse analyses.
The inverse problem is formulated considering the first three undamped natural frequencies and damping ratios of the system. The generalized response vector is defined as follows (36) (36) where and are, respectively, the undamped natural frequency and damping ratio of the rth vibration mode of the system.
Synthetic experimental modal data were computed considering the corrupting effects of noise as follows (37) (37) where is the reference generalized response vector, obtained considering the finite element model of the system and is a zero-mean Gaussian random vector whose covariance matrix is , i.e. .
The modal parameters that compose the reference response vector , presented in Table , were obtained considering the finite element model of the system and the parameters presented in Tables and . Besides, the measurement error covariance matrix was assumed as follows (38) (38) where and are diagonal covariance matrices containing the variances of the measurement errors of the undamped natural frequencies and damping ratios, respectively. A relative standard deviation of 0.3% was adopted for the frequency measurement errors and a relative standard deviation of 3% was adopted for the measurement errors in the damping ratios. The standard deviations are relative to the corresponding reference values for the undamped natural frequencies and damping ratios presented in Table . It is worth mentioning that the relative values adopted for the standard deviations in the present manuscript are in agreement with experimental values found in the literature [Citation54,Citation55]. These assumptions yield GHz and .
It should be highlighted that, in practice, the modal parameters of the system may be extracted from measured time domain data by means of the Eigensystem Realization Algorithm (ERA) [Citation56], the Time Domain Decomposition method (TDD) [Citation57] or the Subspace System Identification (SSI) [Citation58], to cite a few.
Firstly, six cases were addressed in the numerical simulations, as summarized in Table . In Cases 1A, 1B and 1C, the nonlocal and the viscoelastic parameters were considered as unknowns, yielding the vector of unknown parameters . Therefore, for cases 1A, 1B and 1C it is assumed that the diameter of the model, used along the inverse analyses, has a constant value and it is equal to a nominal value . In particular, for case 1A one considers the nominal diameter for the inverse analysis exactly equal to the diameter d of the model that was used to generate the synthetic data, as shown in Table . On the other hand, for cases 1B and 1C one may analyze the impact of the misspecification of the system diameter value on the results provided by the inverse analyses. In Cases 2A, 2B and 2C, the diameter of the system was assumed as unknown and, therefore, it was also considered as a component of the vector of unknown model parameters, i.e. . In this case, as shown next, the prior densities for the diameter d were assumed as Gaussians centred at a nominal diameter value . Different values for the nominal diameter of the SWCNT distinguish Cases 2A, 2B and 2C, as pinpointed in Table . Although the diameter d belongs to the vector of unknown parameters , it should be highlighted that it consists of an uninteresting parameter.
In the numerical results that follow, samples of the posterior probability density function of the unknown parameters were drawn by the DRAM algorithm. A total of states were drawn, where the first were considered as burn-in states. Two stages were considered. During the non-adaptation period states, the covariance matrix of the auxiliary probability density function at the first stage was (provided later, according to the considered case). The adaptation of the auxiliary density of the first stage started at the state and the computation of the covariance matrix did not take into account the first states of the burn-in period. At the second stage of the DRAM algorithm, which represents the DR step, the auxiliary density was obtained by scaling the auxiliary density of the first stage by a shrinking factor , as described in Subsection 3.1. The prior probability density functions and the initial states of the Markov chains are also provided later, according to the considered case.
4.2.1. Cases 1A, 1B and 1C
Numerical results regarding Cases 1A, 1B and 1C described in Table are presented here. In these cases, the vector of unknown parameters is given by .
The covariance matrix of the non-adaptation period of the DRAM was adopted as , with nm, and rad/s. During the adaptation period, the parameter was adopted, Equation (Equation30(30) (30) ).
As for the prior density, it was assumed that all model parameters are independent random variables whose marginal densities are: , and . Here, is used to denote a uniform PDF defined within the interval .
Figures present the Markov chains and corresponding marginal posterior PDFs of the model parameters for Cases 1A, 1B and 1C.
From Figures to , one may observe that the chains converged to their stationary distributions and presented good mixing. In all cases considered in the present work, the convergence and mixing of the chains were verified by inspection of the cumulative mean values and autocorrelation functions of the unknown parameters, not shown here for the sake of space. Figure (a,b) show that the misspecification of the system diameter yielded extremely slight variations in the estimates of the nonlocal parameter , whose marginal PDFs presented approximately the same dispersion, in all cases, and all of them encompass the exact value of this parameter. Figures – pinpoint relatively biased solutions for the viscoelastic parameters Δ and Ω due to the adoption of a diameter , for the inverse model, that is different from the one used to generate the synthetic experimental data. Although the PDFs of these parameters presented approximately the same dispersion, only those obtained from Case 1A encompass the corresponding exact values of these parameters.
Table presents the exact value of each estimated parameter, some corresponding sample statistical properties and the maximum a posteriori (MAP) estimate, obtained for Cases 1A, 1B and 1C. The considered sample statistics are: sample mean (μ), sample standard deviation (σ), coefficient of variation () and 95% credibility interval (). The MAP estimate was obtained considering the Nelder-Mead optimization method [Citation59].
Relative errors around 18% and , between the sample mean estimates and exact parameter values, for Δ and Ω, respectively, were obtained in Case 1B, when an error of in the system diameter was considered. Relative erros around and for Δ and Ω, respectively, were obtained in Case 1C, where an error of in the system diameter was considered. Besides, the marginal PDFs of these parameters presented the greater coefficients of variation and the credibility intervals did not encompass their corresponding exact values.
The marginal posterior distributions and the scatter plots of the model parameters, for Cases 1A, 1B and 1C, are presented in Figure , where one may observe a high correlation between the nonlocal parameter and the viscoelastic parameter Δ. This correlation may be explained by the fact that, according to Figures –, an increase in yields a decrease in the values of the natural frequencies, whereas, on the other hand, the viscoelastic parameter Δ is related to the stiffness of the system, as can be observed in Equation (Equation22(22) (22) ), such that an increase in its value yields also an increase in the natural frequencies of the system. Therefore, in sampling the posterior probability density function of the unknown parameters, greater (smaller) values of may be compensated by greater (smaller) values of Δ, and vice versa, as can be seen by the positive correlation between them in Figure .
4.2.2. Cases 2A, 2B and 2C – addressing misspecifications of the SWCNT diameter
In the results that follows, the diameter of the system is assumed uncertain and the vector of unknown parameters is defined as , according to Table .
The covariance matrix of the non-adaptation period of the DRAM was adopted as , with nm, , rad/s and nm. During the adaptation period, the parameter was adopted, Equation (Equation30(30) (30) ).
The same prior PDFs were adopted for the nonlocal and viscoelastic parameters: , and . Regarding the prior for d, it was assumed as .
Figures , and present the Markov chains and marginal PDFs for the interesting parameters for cases 2A, 2B and 2C. It is interesting to note that the level of uncertainty provided by the marginal posteriors is approximately the same as the one provided by the marginal posteriors in cases 1A, 1B and 1C, as shown in Figure .
As for the viscoelastic parameters , Figures and indicate a larger level of uncertainty for them when compared to the one provided by the posteriors in cases 1A, 1B and 1C, as shown in Figures and . Finally, Figure presents the chains and marginal densities for the uninteresting parameter d.
Figure depicts the marginal posterior distributions, along with the scatter plots of the model parameters, for Cases 2A, 2B and 2C. High correlations between the parameters d, Δ and Ω may be clearly observed in Figure . The correlation between the parameters d and Δ may be explained as follows. The diameter d of the nanotube, similar to the viscoelastic parameter Δ, is related to the stiffness of the system. The higher the value of these parameters, the greater the stiffness and, consequently, the greater the natural frequencies of the system. Therefore, in sampling the posterior probability density function of the parameters, greater (smaller) values of d may be compensated by smaller (greater) values of Δ, and vice versa, yielding a negative correlation between them.
Table presents some statical properties from cases 2A, 2B and 2C. An important aspect concerning the information presented in Table and in the marginal PDFs shown in Figures , and is that, despite the fact that the level of uncertainty associated to the blue model parameters has increased, when compared to cases 1A, 1B and 1C, the true parameter values are within the feasible values provided by the posterior solutions of Cases 2A, 2B and 2C.
From these analyses one may observe that slight levels of misspecification in the system diameter may lead to significant biased estimates for the viscoelastic parameters , as shown in Figures and , when one does not consider the diameter as an uncertain model parameter. At the same time, when assuming the diameter as an unknown but uninteresting model parameter, this biased characteristic could be properly alleviated.
It is worth mentioning the fact that, in this subsection, the model used for inversion analyses is able to describe the entire physical mechanisms associated to the model that provided the set of synthetic data. Nevertheless, in practice one will likely face situations where the model used for inversion is not able to completely describe the damping/dissipative mechanisms present in the system. This issue will be approached in Subsection 4.3, shown next using a preliminary approach.
4.3. Further analysis
In order to simulate a more realistic situation, in the following analysis, the nonlocal viscoelastic beam model with N = 1 internal variable is calibrated in the presence of both measurement and modelling errors. The modelling error here stems from the incomplete modelling of the damping mechanisms of the system. More specifically, one considers a model with internal variables to provide the reference data for computing the synthetic experimental data , whereas a forward model with N = 1 internal variable is used for the inverse analysis.
Given the fact that the model used to derive the synthetic data and the one used in the inverse analysis are different, one may rewrite the observation model, relating observations and model predictions, as follows [Citation60] (39) (39) where is an accurate forward operator able to completely describe the physical mechanisms associated to the system under analysis, such that discrepancies between measured data and model predictions could be completely described by measurement errors . The accurate forward operator in Equation (Equation39(39) (39) ) requires the specification of its own set of parameters . As for the right side of Equation (Equation39(39) (39) ), it explicitly states that if one decides to write the observation model using a less accurate predictive model , one should take into account the modelling errors, defined as . Besides, according to Equation (Equation39(39) (39) ), one has that the difference between measured data and predictions provided by a less accurate operator depends on both measurement errors and modelling errors . One aspect that should be emphasized is the fact that the modelling errors depend on the structure of the model and also on the structure of the accurate model .
In the present case, as the models and belong to different model classes, inasmuch as is built with internal variables and is built with N = 1 internal variable, a simplifying hypothesis was adopted here which proved amenable for the present application. It was assumed that the effects of the measurement errors and modelling errors could be globally described by a single random variable , as shown in Equation (Equation39(39) (39) ), and that it follows a zero mean Gaussian density function, i.e. , where is unknown. Further, for the inverse analysis it is assumed that and the global error are independent random variables, such that the likelihood function may be written as (40) (40) Finally, given the fact that the vector contains information of measured natural frequencies and damping ratios, one decided for a parameterization of the covariance matrix as follows (41) (41) where β is an unknown parameter and the matrix is the covariance matrix of the measurement errors shown in Equation (Equation38(38) (38) ).
In the present case (Case 3), the modal parameters that compose the reference response vector were obtained considering a model built with internal variables, whose parameters are presented in Tables and . The reference modal parameters for Case 3 are presented in Table . As in the previous cases, the experimental response vector was computed considering a relative standard deviation of 0.3% for the frequency errors and a relative standard deviation of 3% for the errors in the damping ratios. The standard deviations are relative to the corresponding reference values for the undamped natural frequencies and damping ratios presented in Table , yielding GHz and .
The nonlocal and viscoelastic parameters , along with the nanotube diameter d and the parameter β, related to the covariance matrix of the global error vector , Equation (Equation41(41) (41) ), were considered as unknowns in Case 3, i.e. the vector of unknown parameters is , as indicated in Table .
Concerning the prior probability density function, it was built assuming the unknown parameters as independent random variables whose marginal densities are: , , , and , where stands for the inverse gamma distribution, with shape parameter a and scale parameter b. As for the sampling of the posterior probability density function, it was considered samples where the first ones were considered as burn-in states.
Figure depicts the marginal posterior distributions, along with the scatter plots of the model parameters, obtained for Case 3. Regarding the marginal posterior distributions shown in Figure , one may highlight the fact that comparisons between feasible values and reference ones may be meaningless, except for the nonlocal parameter . This is due to the fact that the nonlocal parameter is the only constitutive model parameter with the same physical meaning in the model with N = 1 internal variable, used for the inverse analysis, and in the model with internal variables, used to generate the synthetic data.
Table presents the statistical properties of the model parameters provided by the inverse analysis. From Table , one may observe that the reference value of the nonlocal parameter is not within the 95% credibility interval. However, the reference value of the nonlocal parameter is within the feasible region provided by the posterior marginal density , as can be seen from Figure . Regarding the scaling parameter β, which is used in Equation (Equation41(41) (41) ) to characterize the covariance of the random errors in Equation (Equation39(39) (39) ), one may observe that its marginal posterior distribution is centred around 4. This fact indicates that the support of the likelihood function had to expand in order to cope with the discrepancies associated to the modelling error that stemmed from the fact that the model used in the inverse analysis is not able to completely describe the dissipation mechanisms associated to the system that provided the synthetic data.
As a final comment concerning Case 3, one may highlight that the adopted approach could be used when the user cannot properly select an accurate model to completely describe the dissipation mechanisms of the system under analysis. In this perspective, this simple approach could be used as an alternative to build an equivalent model for predictive analyses, when considering the simplifying hypothesis of a zero mean Gaussian global error in Equation (Equation39(39) (39) ). From this point of view, the raw examination of the samples may possibly not be useful. Nevertheless, the analysis of model predictions provided by the equivalent model, for diverse conditions, when may prove useful. In this sense, as an example, one decided to analyze the predictions provided by the calibrated model when the quantity of interest is the Frequency Response Function (FRF) of the nonlocal viscoelastic beam, i.e. , where ω is the angular frequency. The FRF relating an excitation at position 0.45 nm with the displacement of the system at 0.85 nm is considered here. Figure depicts the FRF provided by the reference model with internal variables in blue line. The model with N = 1 internal variable is used to provide model predictions , using samples (after discarding the brun-in states) from the posterior probability density obtained for Case 3. Realizations of model predictions are shown in grey lines and the black line indicates its expected values; moreover, dashed lines indicate the 99% stochastic envelope provided by the model predictions. From Figure one may observe that the stochastic envelope provided by the model with N = 1 internal variable encompasses the FRF provided by the reference model with internal variables until the fourth mode of the system, around 1250 GHz.
Concerning the calibration of a model containing N = 1 internal variable, assuming the conventional observation model, i.e. in Equation (Equation41(41) (41) ), although it is not shown here one states that: (i) its model predictions present quite low levels of uncertainty when compared to the ones shown in Figure , and (ii) it fails to properly describe the second and fourth modes of the system under analysis.
Finally, one should notice that the choice for this simplifying approach in the presence of modelling errors is an attempt to build an equivalent predictive model for the system under analysis. More rigorous and comprehensive approaches could be considered aiming at obtaining an accurate calibrated model and its proper parameter estimates. This requires that the user has candidate feasible models that are believed to accurately describe the physical mechanisms inherent to the system under analysis. In this sense, as an example, one may cite the approach proposed by Muto and Beck [Citation61], where a set of feasible models of different classes are assessed and the probability of each model , given a set of measured data , is computed taking into account a balance between its predictive capabilities and its complexity. This approach would be interesting in case the user has, for example, model classes describing different dissipations mechanisms or even different numbers of internal variables, but he is unable, at a first sight, to specify which model would be the most appropriate one to describe the physical phenomena of the system under analysis.
5. Conclusions
The present work presented a model for nonlocal and viscoelastic Euler-Bernoulli beams and some aspects of the calibration of its nonlocal viscoleastic model were addressed. The inverse analyses was performed based on the Bayesian framework. The main contributions of the present work are:
to present an approach to simultaneously estimate the nonlocal and viscoelastic constitutive model parameters;
to assess the impact of a misspecification of the beam diameter d on the estimation of the constitutive parameters and its compensation when considering d as an unknown but uninteresting parameter;
to build an equivalent model to describe more complex damping mechanisms than the one inherent to the model with N = 1 internal variable, based on the hypothesis that the model errors could be described by means of a zero mean Gaussian global error with unknown covariance matrix.
to show that the equivalent predictive model is somehow able to describe the response of the system under analysis when the global error model is calibrated.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Additional information
Funding
References
- Il'ina MV, Il'in OI, Blinov YF, et al. Piezoelectric response of multi-walled carbon nanotubes. Materials. 2018;11(4):638.
- Ageev O, Blinov YF, Il'ina MV, et al. Resistive switching of vertically aligned carbon nanotubes for advanced nanoelectronic devices. In: Tiwari A, Mishra YK, Kobayashi H, Turner APF, editors. Intelligent nanomaterials. Basel: John Wiley & Sons; 2017. p. 361–394.
- Fakher M, Rahmanian S, Hosseini-Hashemi S. On the carbon nanotube mass nanosensor by integral form of nonlocal elasticity. Int J Mech Sci. 2019;150:445–457.
- Karličić D, Kozić P, Adhikari S, et al. Nonlocal mass-nanosensor model based on the damped vibration of single-layer graphene sheet influenced by in-plane magnetic field. Int J Mech Sci. 2015;96–97:132–142.
- Gibson RF, Ayorinde EO, Wen YF. Vibrations of carbon nanotubes and their composites: a review. Compos Sci Technol. 2007;67(1):1–28.
- Kinloch IA, Suhr J, Lou J, et al. Composites with carbon nanotubes and graphene: an outlook. Science. 2018;362(6414):547–553.
- Arash B, Wang Q. A review on the application of nonlocal elastic models in modeling of carbon nanotubes and graphenes. Comput Mater Sci. 2012;51(1):303–313.
- Eltaher MA, Khater ME, Emam SA. A review on nonlocal elastic models for bending, buckling, vibrations, and wave propagation of nanoscale beams. Appl Math Model. 2016;40(5):4109–4128.
- Wang Y, LI F. Dynamical properties of nanotubes with nonlocal continuum theory: a review. Sci China Phys Mech Astron. 2012;55(7):1210–1224.
- Büyüköztürk O, Buehler MJ, Lau D, et al. Structural solution using molecular dynamics: fundamentals and a case study of epoxy-silica interface. Int J Solids Struct. 2011;48(14-15):2131–2140.
- Adhikari S, Gilchrist D, Murmu T, et al. Nonlocal normal modes in nanoscale dynamical systems. Mech Syst Signal Process. 2015;60:583–603.
- Eringen AC. On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. J Appl Phys. 1983;54(9):4703–4710.
- Reddy JN. Nonlocal theories for bending, buckling and vibration of beams. Int J Eng Sci. 2007;45:288–307.
- Wang Q, Varadan VK. Vibration of carbon nanotubes studied using nonlocal continuum mechanics. Smart Mater Struct. 2006;15:659–666.
- Wang CM, Zhang YY, He XQ. Vibration of nonlocal Timoshenko beams. Nanotechnology. 2007;18(10):105401.
- Wang Q, Varadan VK. Stability analysis of carbon nanotubes via continuum models. Smart Mater Struct. 2005;14(1):281.
- Eltaher MA, Alshorbagy AE, Mahmoud FF. Vibration analysis of Euler-Bernoulli nanobeams by using finite element method. Appl Math Model. 2013;37(7):4787–4797.
- Eltaher MA, Emam SA, Mahmoud FF. Free vibration analysis of functionally graded size-dependent nanobeams. Appl Math Comput. 2012;218(14):7406–7420.
- Phadikar JK, Pradhan SC. Variational formulation and finite element analysis for nonlocal elastic nanobeams and nanoplates. Comput Mater Sci. 2010;49(3):492–499.
- Lee H-L, Hsu J-C, Chang W-J. Frequency shift of carbon-nanotube-based mass sensor using nonlocal elasticity theory. Nanoscale Res Lett. 2010;5(11):1774.
- Li X-F, Tang G-J, Shen Z-B, et al. Resonance frequency and mass identification of zeptogram-scale nanosensor based on the nonlocal beam theory. Ultrasonics. 2015;55:75–84.
- De Rosa MA, Lippiello M, Martin HD. Free vibrations of a cantilevered SWCNT with distributed mass in the presence of nonlocal effect. Sci World J. 2015;2015.
- Aydogdu M, Filiz S. Modeling carbon nanotube-based mass sensors using axial vibration and nonlocal elasticity. Physica E. 2011;43(6):1229–1234.
- Murmu T, Adhikari S. Nonlocal frequency analysis of nanoscale biosensors. Sens Actuators A. 2012;173(1):41–48.
- Duan WH, Wang CM, Zhang YY. Calibration of nonlocal scaling effect parameter for free vibration of carbon nanotubes by molecular dynamics. J Appl Phys. 2007;101(2):024305.
- Liang Y, Han Q. Prediction of nonlocal scale parameter for carbon nanotubes. Sci China Phys Mech Astron. 2012;55(9):1670–1678.
- Ansari R, Rouhi H. Analytical treatment of the free vibration of single-walled carbon nanotubes based on the nonlocal Flugge shell theory. J Eng Mater Technol. 2012;134(1):011008.
- Arash B, Ansari R. Evaluation of nonlocal parameter in the vibrations of single-walled carbon nanotubes with initial strain. Physica E. 2010;42(8):2058–2064.
- Tuna M, Kırca M. Unification of Eringen's nonlocal parameter through an optimization-based approach. Mech Adv Mater Struct. 2019. DOI:https://doi.org/10.1080/15376494.2019.1601312
- Ke C, Espinosa HDNanoelectromechanical systems and modeling. In: Handbook of theoretical and computational nanotechnology. American Scientific Publishers. Vol. 1. 2005. p. 1–38.
- Calleja M, Kosaka PM, San Paulo Á. Challenges for nanomechanical sensors in biological detection. Nanoscale. 2012;4(16):4925–4938.
- Tamayo J, Kosaka PM, Ruz JJ, et al. Biosensors based on nanomechanical systems. Chem Soc Rev. 2013;42(3):1287–1311.
- Lim T-C, editor. Nanosensors: theory and applications in industry, healthcare and defense. Boca Raton: CRC Press; 2016.
- Lee JL, Lin CB. The magnetic viscous damping effect on the natural frequency of a beam plate subject to an in-plane magnetic field. J Appl Mech. 2010;77(1):011014.
- Chen C, Ma M, Zhe Liu J, et al. Viscous damping of nanobeam resonators: humidity, thermal noise, and a paddling effect. J Appl Phys. 2011;110(3):034320.
- Lei Y, Murmu T, Adhikari S, et al. Dynamic characteristics of damped viscoelastic nonlocal Euler-Bernoulli beams. Eur J Mech A Solids. 2013;42:125–136.
- Lei Y, Adhikari S, Friswell MI. Vibration of nonlocal Kelvin-Voigt viscoelastic damped Timoshenko beams. Int J Eng Sci. 2013;66:1–13.
- Ansari R, Oskouie MF, Sadeghi F, et al. Free vibration of fractional viscoelastic Timoshenko nanobeams using the nonlocal elasticity theory. Physica E. 2015;74:318–327.
- Cajić M, Karličić D, Lazarević M. Nonlocal vibration of a fractional order viscoelastic nanobeam with attached nanoparticle. Theor Appl Mech. 2015;42(3):167–190.
- Zhang DP, Lei Y, Wang CY, et al. Vibration analysis of viscoelastic single-walled carbon nanotubes resting on a viscoelastic foundation. J Mech Sci Technol. 2017;31(1):87–98.
- Zhang DP, Lei Y, Shen ZB. Vibration analysis of horn-shaped single-walled carbon nanotubes embedded in viscoelastic medium under a longitudinal magnetic field. Int J Mech Sci. 2016;118:219–230.
- Zhang DP, Lei Y, Shen ZB. Free transverse vibration of double-walled carbon nanotubes embedded in viscoelastic medium. Acta Mech. 2016;227(12):3657–3670.
- Kaipio J, Somersalo E. Statistical and computational inverse problems. New York: Springer Science & Business Media; 2006.
- Castello DA, Rochinha FA, Roitman N, et al. Constitutive parameter estimation of a viscoelastic model with internal variables. Mech Syst Signal Process. 2008;22(8):1840–1857.
- Borges FCL, Castello DA, Magluta C, et al. An experimental assessment of internal variables constitutive models for viscoelastic materials. Mech Syst Signal Process. 2015;50:27–40.
- Vasques CMA, Moreira RAS, Rodrigues JD. Viscoelastic damping technologies-part I: modeling and finite element implementation. J Adv Res Mech Eng. 2010;1(2).
- Haario H, Laine M, Mira A, et al. DRAM: efficient adaptive MCMC. Stat Comput. 2006;16(4):339–354.
- Christensen R. Theory of viscoelasticity. New York: Elsevier; 2012.
- Wineman AS, Kumbakonam RR. Mechanical response of polymers: an introduction. Cambridge: Cambridge University Press; 2000.
- Inman DJ. Engineering vibration. Upper Saddle River, NJ: Prentice Hall; 2001.
- Reddy JN. An introduction to the finite element method. New York: McGraw-Hill; 1993.
- Haario H, Saksman E, Tamminen J. An adaptive metropolis algorithm. Bernoulli. 2001;7(2):223–242.
- Mira A. On metropolis-Hastings algorithms with delayed rejection. Metron. 2001;59(3-4):231–241.
- Mellinger P, Döhler M, Mevel L. Variance estimation of modal parameters from output-only and input/output subspace-based system identification. J Sound Vib. 2016;379:1–27.
- Souza M, Castello D, Roitman N, et al. Impact of damping models in damage identification. Shock Vib. 2019;2019:Article ID 4652328.
- Juang JN, Pappa RS. An eigensystem realization algorithm for modal parameter identification and model reduction. J Guid Control Dyn. 1985;8(5):620–627.
- Byeong HK, Stubbs N, Park T. A new method to extract modal parameters using output-only responses. J Sound Vib. 2005;282:215–230.
- Overschee PV, De Moor BL. Subspace identification for linear systems, theory, implementation, applications. Dordrecht: Kluwer Academic Publishers; 1996.
- Nelder JA, Mead R. A simplex method for function minimization. Comput J. 1965;7(4):308–313.
- Kolehmainen V, Tarvainen T, Arridge SR, et al. Marginalization of unisteresting distributed parameters in inverse problems – application to diffuse optical tomography. Int J Uncertain Quantif. 2011;1(1):1–7.
- Muto M, Beck JL. Bayesian updating and model class selection for hysteretic structural models using stochastic simulation. J Vib Control. 2008;14(1–2):7–34.
Appendices
Appendix 1.
Finite element matrices and vectors
The vector of generalized element nodal displacements and the vector of internal variables are, respectively, defined as (A1) (A1) (A2) (A2) where and are transverse displacements and and are rotations at the element nodes.
The matrix with the interpolation functions is given by (A3) (A3) The element mass matrix is defined as (A4) (A4) The element viscous damping matrix is defined as (A5) (A5) The relaxed stiffness matrix is defined as (A6) (A6) The vector of generalized nodal forces is defined as (A7) (A7)
Appendix 2.
State-space model and system poles
Taking into account the elemental equations in (Equation22(22) (22) ), the equations of motion of a nonlocal viscoelastic Euler-Bernoulli beam can be written as (A8) (A8) Defining the augmented displacement vector as (A9) (A9) the coupled equations in (EquationA8(A8) (A8) ) may be expressed in the compact form [Citation46] (A10) (A10) where (A11) (A11) with (A12) (A12) Defining the augmented state vector (A13) (A13) the second-order equation of motion in Equation (EquationA10(A10) (A10) ) can be written in the first-order state-space form (A14) (A14) where the augmented state matrix and the augmented input influence matrix are given, respectively, by (A15) (A15) Considering Equation (EquationA14(A14) (A14) ), the poles of the system are defined as the roots of the characteristic equation (A16) (A16)