Abstract
In this article, we present a new contact resonance atomic force microscopy-based method utilizing a square, plate-like microsensor to accurately estimate viscoelastic sample properties. A theoretical derivation, based on Rayleigh–Ritz method and on an “unconventional” generalized eigenvalue problem, is presented and a numerical experiment is devised to verify the method. We present an updated sensitivity criterion that allows users, given a set of measured in-contact eigenfrequencies and modal damping ratios, to select the best eigenfrequency for accurate data estimation. The verification results are then presented and discussed. Results show that the proposed method performs extremely well in the identification of viscoelastic properties over broad ranges of nondimensional sample stiffness and damping values.
1 Introduction
Contact resonance atomic force microscopy (CR-AFM), originally pioneered by Rabe et al. [1], uses the measured freely vibrating and in-contact resonance frequencies of a coupled microsensor-sample system to estimate material properties of the sample. Contact resonance (CR) spectroscopy AFM techniques are an ideal candidate for nanoscale quantitative mechanical property extraction due to their avoidance of nonlinear interactions present in other AFM measurement techniques. CR-AFM operates primarily in the linear net-repulsive regime of the tip-sample interaction force curve. CR-AFM has been used to measure both elastic and viscoelastic sample properties [2–6]. Recently, Aureli et al. [7] showed that by using a square, plate-like cantilevered microsensor, as opposed to the traditional rectangular, narrow beam AFM cantilevered sensor, several advantages arise. These advantages include an increase of sensing modes available in a given frequency bandwidth and the ability to optimally place the sensing tip within the plate domain to maximize measurement sensitivity.
In this study, we expand the method presented in Ref. [7], which was formulated to test purely elastic samples, to include viscoelastic materials. With this addition, the aforementioned measurement benefits of the proposed sensor can be extended to include the measurement of many viscoelastic materials of interest, including, for example, polymer composites and biological materials. In addition, with the inclusion of damping in the system, it becomes possible to characterize and account for hydrodynamic forces present in liquid measurement environments, such as those encountered when conducting AFM measurements of live biological samples ex vivo.
To incorporate viscoelastic effects, Hamilton’s least action principle is used to create a variational formulation of the problem. The strain energy of the plate is calculated following the assumptions of the Kirchhoff–Love plate equations. The sample contact, originally modeled in Ref. [7] as a linear spring element, is replaced with a Kelvin–Voigt element. We include an additional term in the variational formulation of the method that accounts for the nonconservative work done by the new viscous damper added to the model. A Rayleigh–Ritz method is used to transform the problem into a generalized eigenvalue problem (EVP) with complex eigenvalues. A nontraditional interpretation of the eigenvalue problem leads to an efficient method for the estimation of the sample viscoelasticity. Specifically, given the measured in-contact natural frequencies and modal damping ratios of the coupled sensor-sample system, the solution of the generalized EVP yields estimates of the viscoelastic material properties, including stiffness and damping.
Following the derivation of the theoretical model, a numerical experiment, using the finite element method, is created to verify the method. Given a set of measured in-contact eigenmodes, a criterion is proposed to select which of these frequencies should be used to produce the most accurate material property estimates. In addition, we discuss methods to extract the model input data from experimentally measured frequency response functions of the system. Finally, the verification results are presented over a broad range of parameters of interest, and future research directions of the work are discussed.
2 Theoretical Framework
2.1 Governing Equations.
In this study, we extend the theoretical formalism developed in our previous paper [7], which is concerned with purely elastic response of the sample, to incorporate the effect of viscoelasticity in the sample mechanical behavior. As in Ref. [7], we will consider the free vibrations of a thin rectangular cantilevered Kirchhoff–Love plate as illustrated in Fig. 1. The plate is clamped at its edge at x = 0 and free at the edges x = Lx and y = ±Ly/2, see Fig. 1. In-plane dimensions of the plate are its length Lx and width Ly, and the plate thickness is indicated with h. We assume that the mechanical constitutive behavior of the plate is linear elastic, isotropic and homogeneous, with Young’s modulus E and Poisson’s ratio , and we denote its mass density per unit volume with .
The plate is in contact through an attachment point of coordinates (xs, ys) with a spring-damper element, which models the Kelvin–Voigt viscoelastic behavior of a sample being analyzed via contact resonance AFM. For simplicity, both the elastic and the viscous behaviors are assumed linear and fully described by a single constant each, that is, the sample stiffness ks and the sample damping coefficient cs. We restrict our study to small oscillations, so that the sample is assumed to react to displacements of the attachment point in the z-direction only.
2.2 The Rayleigh–Ritz Eigenvalue Problem.
By following our original derivation in Ref. [7], we now develop a Rayleigh–Ritz method [11,12] to estimate the complex eigenvalues of Eq. (4). We assume the following ansatz for the motion of the plate , where W(x, y) indicates the complex displacement amplitude and is a separation constant. This important difference from Ref. [7] is necessary to capture both frequency and damping ratio for the separable solutions of Eq. (4). Note, in fact, that the system studied in Ref. [7] is undamped so it was therein assumed that .
2.3 Estimation of the Sample Viscoelasticity.
While Eq. (14) is typically used to calculate natural frequencies and modal damping ratios as described earlier, in this work, the model developed will be used for the identification of the sample stiffness parameter α and viscosity parameter in Eq. (3). Specifically, following in spirit the practice of CR-AFM measurements in which sample properties are estimated from a set of known frequency measurements, here we are concerned with the estimation of α and when , , and are known, along with Ω and a set of measured , with i = 1, …M ≪ N. As in Ref. [7], Ω can, in principle, be identified from measurements of the free unsprung vibrations of the plate.
2.4 Some Implementation Remarks.
Construction of the Φi(x, y) functions as the Cartesian product of “modes” and , their numbering scheme, and further properties are discussed in Ref. [7]. The integrals in Eqs. (11) and (13) are approximated by using Gauss–Legendre quadrature points [17]. The complete details of the implementation are extensively discussed in Ref. [7]. We only desire to mention here that the proposed method results in a reduced-order model for the plate system, including only a relatively small number of modes (with much smaller matrix sizes than a comparably accurate finite element model) capable of accurate prediction of the dynamic characteristics of the system. As an example, when 20 modes are extracted from an N = (15 × 15) size model, calculation run-time is less than 2 s on a typical desktop personal computer, with approximately 75% of the time required to calculate the system matrices.
3 Numerical Experiments
3.1 Finite Element Analysis Simulations.
In lieu of CR-AFM experimental measurements, we thoroughly test our method for sample viscoelasticity identification via numerical experiments. Our numerical campaign is centered around a square plate-like cantilevered microscale sensor of dimensions 225 × 225 × 3 μm3, similar to one of the prototype plate sensors considered in Ref. [7]. The plate material is silicon, with Young’s modulus 169 GPa, Poisson’s ratio , and density . The plate is connected with a linear viscoelastic element located at (xs, ys) = (0.9Lx, 0.4Ly). This location is close to, but not coincident with, the optimal location determined in Ref. [7] for maximum sensitivity for a square plate. While this location is chosen for its favorable sensitivity properties, the choice has to be regarded as purely illustrative, as the application of the method would be unchanged if a different attachment point were to be selected. The viscoelastic element can be assigned directly a stiffness ks and a viscous damping coefficient cs. A model for finite element analysis (FEA) is created for the system using the commercial software package ansys mechanical apdl v.17.0. Shell elements “SHELL63” are utilized to create the structure with a mapped quadrilateral mesh. The element size is selected as Lx/100 to ensure accuracy in the numerical results. Membrane effects are not considered. One combination element “COMBIN14” is used to model the viscoelastic sample. A comprehensive convergence study for the model was already presented in our previous work [7].
For each analyzed case, a known stiffness and damping are assigned to the linear viscoelastic element modeling, and the resulting “in-contact” natural frequencies and modal damping ratios of the system are then calculated via the software “MODAL” analysis solution routine. Because of the explicit presence of the damping in the viscoelastic element, the “DAMP” (damped system) solution method is selected. The output of this analysis consists of a set of complex eigenvalues, which essentially play the role of the in Eq. (14). Thus, the output from the FEA is used in our algorithm described in Sec. 2 for estimating α and . After the modal analysis is completed, we use the FEA model to conduct the identification of the system frequency response function (FRF) via the software “HARMIC" analysis solution routine. Specifically, we study the driving point FRF by applying a harmonic unit force F at the attachment point (xs, ys) and determining the harmonic displacement response p at the same point (xs, ys) over the range of frequencies f ∈ (0, fmax] of interest. This choice is representative and, certainly, not the only possible one. Further study is necessary to understand the potential of nondriving point FRFs in the context of viscoelasticity identification.
For the numerical study, we choose a measurement bandwidth of fmax = 2 MHz to represent frequency limitations present in real experimental systems, similar to Ref. [7]. We use 1000 equally spaced sampling points in the band 0–2 MHz to obtain a fine representation of the system FRF with complex displacement data every 2 kHz. Then, in the neighborhood of each of the resonance peaks identified from the FRF (and coincident with the output of the previous modal analysis), we conduct a second scan with another 100 frequency points to obtain a very fine FRF of the system in such a way as to obtain at least ten sample points in the 3dB bandwidth of each resonance peak, especially for lightly damped modes.
Note that, in our previous analysis in Ref. [7], we have shown how the parameter Ω can be directly estimated from the eigenvalues of the unsprung, freely vibrating system from the FEA simulations. This procedure is analogous to the method used in traditional CR studies and, rather than trying to identify the individual geometric and material parameters separately, we choose to robustly identify the important governing parameter Ω. Therefore, in the remainder of the article, the value of Ω ≈ 1.5079 × 105 s−1, calculated from nominal geometric and material parameters of the plate system, will be presumed known.
3.2 Eigenvalue Extraction From the Frequency Response Function.
In a real experimental measurement, obtained, for example, by contact resonance atomic force microscopy, the frequency response function of the system is the primary measurand. In the proposed method, we require the complex eigenvalues of the in-contact system to generate material property estimates. Here, we show that the eigenvalues of the system can be effectively extracted from measured FRFs using a single-mode approximation, provided that the modes are not too closely spaced. Although, in the following, the procedure will be illustrated on a numeric dataset, it could be implemented seamlessly on actual experimental data.
As stated earlier, the free response of the plate is given by . Assuming that the modes of the system response are well spaced, we approximate each mode as a single-degree-of-freedom (SDOF) oscillator, whose free response is given by + , where C± are constants determined by the initial conditions, is the SDOF natural frequency, and is the damping ratio [13].

(Top) Sample FRF for α = 100, , and . (Bottom) Comparison between the eigenvalues generated by FEA through modal analysis versus FRF extraction. Close-up view of the first mode is shown. The complex eigenvalues can be accurately estimated from the FRF of the system.
Table 1 presents the extraction results for several combinations of assigned α and . In general, our results demonstrate that, as long as the mode is adequately sampled (>10 sample points in the 3 dB bandwidth), the eigenvalues can be reliably calculated. Having shown the applicability of the extraction of eigenvalues from the system FRF, in the rest of this study, we will directly use the eigenvalues calculated through the FEA modal analysis routine as an input for our estimation procedure.
Eigenvalues extracted from FRFs for several combinations of α and for various eigenmodes and comparison with results from FEA modal analysis
FRF fit | FEA value | % Difference | ||||||
---|---|---|---|---|---|---|---|---|
α | Mode | |||||||
0.01 | 1.20 × 10−03 | 4 | 0.65757 | 1.36 × 10−04 | 0.65757 | 1.36 × 10−04 | 0.00 | −0.06 |
0.1 | 1.20 × 10−03 | 7 | 1.4743 | 2.10 × 10−05 | 1.4743 | 2.11 × 10−05 | 0.00 | 0.42 |
0.1 | 1.20 × 10−02 | 3 | 0.51682 | 1.84 × 10−04 | 0.51682 | 1.83 × 10−04 | 0.00 | −0.21 |
1 | 1.20 × 10−03 | 2 | 0.22845 | 2.21 × 10−04 | 0.22845 | 2.21 × 10−04 | 0.00 | 0.00 |
1 | 1.20 × 10−02 | 1 | 0.10389 | 7.32 × 10−04 | 0.10386 | 7.36 × 10−04 | 0.03 | 0.47 |
10 | 1.20 × 10−03 | 6 | 0.80441 | 2.68 × 10−04 | 0.80441 | 2.68 × 10−04 | 0.00 | 0.00 |
10 | 1.20 × 10−02 | 7 | 1.3343 | 8.03 × 10−04 | 1.3343 | 8.04 × 10−04 | 0.00 | 0.12 |
100 | 1.20 × 10−02 | 5 | 1.1723 | 6.01 × 10−04 | 1.1723 | 6.01 × 10−04 | 0.00 | −0.02 |
100 | 1.20 × 10+01 | 1 | 0.14219 | 4.04 × 10−04 | 0.14219 | 4.04 × 10−04 | 0.00 | 0.00 |
1000 | 1.20 × 10+00 | 5 | 1.2413 | 3.16 × 10−04 | 1.2413 | 3.16 × 10−04 | 0.00 | 0.01 |
FRF fit | FEA value | % Difference | ||||||
---|---|---|---|---|---|---|---|---|
α | Mode | |||||||
0.01 | 1.20 × 10−03 | 4 | 0.65757 | 1.36 × 10−04 | 0.65757 | 1.36 × 10−04 | 0.00 | −0.06 |
0.1 | 1.20 × 10−03 | 7 | 1.4743 | 2.10 × 10−05 | 1.4743 | 2.11 × 10−05 | 0.00 | 0.42 |
0.1 | 1.20 × 10−02 | 3 | 0.51682 | 1.84 × 10−04 | 0.51682 | 1.83 × 10−04 | 0.00 | −0.21 |
1 | 1.20 × 10−03 | 2 | 0.22845 | 2.21 × 10−04 | 0.22845 | 2.21 × 10−04 | 0.00 | 0.00 |
1 | 1.20 × 10−02 | 1 | 0.10389 | 7.32 × 10−04 | 0.10386 | 7.36 × 10−04 | 0.03 | 0.47 |
10 | 1.20 × 10−03 | 6 | 0.80441 | 2.68 × 10−04 | 0.80441 | 2.68 × 10−04 | 0.00 | 0.00 |
10 | 1.20 × 10−02 | 7 | 1.3343 | 8.03 × 10−04 | 1.3343 | 8.04 × 10−04 | 0.00 | 0.12 |
100 | 1.20 × 10−02 | 5 | 1.1723 | 6.01 × 10−04 | 1.1723 | 6.01 × 10−04 | 0.00 | −0.02 |
100 | 1.20 × 10+01 | 1 | 0.14219 | 4.04 × 10−04 | 0.14219 | 4.04 × 10−04 | 0.00 | 0.00 |
1000 | 1.20 × 10+00 | 5 | 1.2413 | 3.16 × 10−04 | 1.2413 | 3.16 × 10−04 | 0.00 | 0.01 |
4 Results and Discussion
The results of the estimations, using the FEA calculated eigenvalues and Eqs. (16) and (17), are shown in Table 2, which reports the percent difference between the estimated values and assigned values in the two-dimensional (α, ) space explored.
Predicted and versus assigned values α and , with percent differences (%D)
α | Mode | %D α | %D | |||
---|---|---|---|---|---|---|
0.01 | 1.20 × 10−06 | 1 | 9.71 × 10−03 | 1.20 × 10−06 | −2.9 | 0.0 |
0.01 | 1.20 × 10−05 | 1 | 9.71 × 10−03 | 1.20 × 10−05 | −2.9 | 0.0 |
0.01 | 1.20 × 10−04 | 1 | 9.71 × 10−03 | 1.20 × 10−04 | −2.9 | 0.0 |
0.01 | 1.20 × 10−03 | 1 | 9.71 × 10−03 | 1.20 × 10−03 | −2.9 | 0.0 |
0.1 | 1.20 × 10−05 | 1 | 9.97 × 10−02 | 1.20 × 10−05 | −0.3 | 0.0 |
0.1 | 1.20 × 10−04 | 1 | 9.97 × 10−02 | 1.20 × 10−04 | −0.3 | 0.0 |
0.1 | 1.20 × 10−03 | 1 | 9.97 × 10−02 | 1.20 × 10−03 | −0.3 | 0.0 |
0.1 | 1.20 × 10−02 | 1 | 9.97 × 10−02 | 1.20 × 10−02 | −0.3 | 0.0 |
1 | 1.20 × 10−04 | 2 | 9.98 × 10−01 | 1.20 × 10−04 | −0.2 | 0.0 |
1 | 1.20 × 10−03 | 2 | 9.98 × 10−01 | 1.20 × 10−03 | −0.2 | 0.0 |
1 | 1.20 × 10−02 | 2 | 9.98 × 10−01 | 1.20 × 10−02 | −0.2 | 0.0 |
1 | 1.20 × 10−01 | 2 | 9.99 × 10−01 | 1.20 × 10−01 | −0.1 | 0.0 |
10 | 1.20 × 10−03 | 2 | 9.98 × 10+00 | 1.20 × 10−03 | −0.2 | −0.4 |
10 | 1.20 × 10−02 | 2 | 9.98 × 10+00 | 1.20 × 10−02 | −0.2 | −0.4 |
10 | 1.20 × 10−01 | 2 | 9.98 × 10+00 | 1.20 × 10−01 | −0.2 | −0.4 |
10 | 1.20 × 10+00 | 4 | 1.01 × 10+01 | 1.21 × 10+00 | 0.9 | 0.5 |
100 | 1.20 × 10−02 | 5 | 9.87 × 10+01 | 1.17 × 10−02 | −1.3 | −2.6 |
100 | 1.20 × 10−01 | 5 | 9.87 × 10+01 | 1.17 × 10−01 | −1.3 | −2.6 |
100 | 1.20 × 10+00 | 5 | 9.93 × 10+01 | 1.17 × 10+00 | −0.7 | −2.5 |
100 | 1.20 × 10+01 | 1 | 9.86 × 10+01 | 1.13 × 10+01 | −1.4 | −5.5 |
1000 | 1.20 × 10−01 | 2 | 8.47 × 10+02 | 8.62 × 10−02 | −15.3 | −28.1 |
1000 | 1.20 × 10+00 | 2 | 8.47 × 10+02 | 8.61 × 10−01 | −15.3 | −28.2 |
1000 | 1.20 × 10+01 | 2 | 8.55 × 10+02 | 8.64 × 10+00 | −14.5 | −28.0 |
1000 | 1.20 × 10+02 | 2 | 1.40 × 10+03 | 7.61 × 10+01 | 40.0 | −36.6 |
α | Mode | %D α | %D | |||
---|---|---|---|---|---|---|
0.01 | 1.20 × 10−06 | 1 | 9.71 × 10−03 | 1.20 × 10−06 | −2.9 | 0.0 |
0.01 | 1.20 × 10−05 | 1 | 9.71 × 10−03 | 1.20 × 10−05 | −2.9 | 0.0 |
0.01 | 1.20 × 10−04 | 1 | 9.71 × 10−03 | 1.20 × 10−04 | −2.9 | 0.0 |
0.01 | 1.20 × 10−03 | 1 | 9.71 × 10−03 | 1.20 × 10−03 | −2.9 | 0.0 |
0.1 | 1.20 × 10−05 | 1 | 9.97 × 10−02 | 1.20 × 10−05 | −0.3 | 0.0 |
0.1 | 1.20 × 10−04 | 1 | 9.97 × 10−02 | 1.20 × 10−04 | −0.3 | 0.0 |
0.1 | 1.20 × 10−03 | 1 | 9.97 × 10−02 | 1.20 × 10−03 | −0.3 | 0.0 |
0.1 | 1.20 × 10−02 | 1 | 9.97 × 10−02 | 1.20 × 10−02 | −0.3 | 0.0 |
1 | 1.20 × 10−04 | 2 | 9.98 × 10−01 | 1.20 × 10−04 | −0.2 | 0.0 |
1 | 1.20 × 10−03 | 2 | 9.98 × 10−01 | 1.20 × 10−03 | −0.2 | 0.0 |
1 | 1.20 × 10−02 | 2 | 9.98 × 10−01 | 1.20 × 10−02 | −0.2 | 0.0 |
1 | 1.20 × 10−01 | 2 | 9.99 × 10−01 | 1.20 × 10−01 | −0.1 | 0.0 |
10 | 1.20 × 10−03 | 2 | 9.98 × 10+00 | 1.20 × 10−03 | −0.2 | −0.4 |
10 | 1.20 × 10−02 | 2 | 9.98 × 10+00 | 1.20 × 10−02 | −0.2 | −0.4 |
10 | 1.20 × 10−01 | 2 | 9.98 × 10+00 | 1.20 × 10−01 | −0.2 | −0.4 |
10 | 1.20 × 10+00 | 4 | 1.01 × 10+01 | 1.21 × 10+00 | 0.9 | 0.5 |
100 | 1.20 × 10−02 | 5 | 9.87 × 10+01 | 1.17 × 10−02 | −1.3 | −2.6 |
100 | 1.20 × 10−01 | 5 | 9.87 × 10+01 | 1.17 × 10−01 | −1.3 | −2.6 |
100 | 1.20 × 10+00 | 5 | 9.93 × 10+01 | 1.17 × 10+00 | −0.7 | −2.5 |
100 | 1.20 × 10+01 | 1 | 9.86 × 10+01 | 1.13 × 10+01 | −1.4 | −5.5 |
1000 | 1.20 × 10−01 | 2 | 8.47 × 10+02 | 8.62 × 10−02 | −15.3 | −28.1 |
1000 | 1.20 × 10+00 | 2 | 8.47 × 10+02 | 8.61 × 10−01 | −15.3 | −28.2 |
1000 | 1.20 × 10+01 | 2 | 8.55 × 10+02 | 8.64 × 10+00 | −14.5 | −28.0 |
1000 | 1.20 × 10+02 | 2 | 1.40 × 10+03 | 7.61 × 10+01 | 40.0 | −36.6 |
Note: The eigenmode used for the prediction, selected via the maximum sensitivity criterion in Eq. (18), is indicated.
Table 2 also indicates the specific mode used for the estimation, selected by the maximum sensitivity criterion given in Eq. (18). We see that the error of the estimation is below for values of α ≤ 100 for all combinations of . For α ≥ 1000, the error in the estimation begins to increase, becoming successively larger with increasing . Interestingly, the mode number used does not increase monotonically when using the sensitivity criterion in Eq. (18). For example, for α = 100, , the most sensitive mode is mode 1, while for α = 1000, the criterion in Eq. (18) selects mode 2.
The selection of the form in Eq. (18) for the sensitivity is inspired by the results in Ref. [7], in which, however, the eigenvalues are purely imaginary numbers since the system is undamped. Thus, it is of interest to evaluate whether other sensitivity parameters, leveraging the complex nature of the present eigenvalues, may offer a better estimation performance. To this aim, Table 3 presents the prediction results for α and using the mode with maximum sensitivity using . We see that by using this sensitivity criterion, the predictions are less accurate compared with those using Eq. (18) in nearly all cases. In particular, the α estimates for assigned values of α ≤ 0.1 are extremely poor. On the other hand, the α estimates for assigned values of 1 ≤ α ≤ 100 are comparable to estimates using Eq. (18), except for the α = 100, case. In general, estimates perform well for assigned values of α ≤ 100. For large assigned values of α, the predictions for both α and diverge. In addition, and very interestingly, the mode selected by this criterion is much different than those selected by Eq. (18), possibly because of how the natural frequency of higher modes enters in the proposed definition of sensitivity.
Predicted and versus assigned values and , with percent differences (%D), using the eigenmode (indicated) with maximum sensitivity
α | Mode | %D α | %D | |||
---|---|---|---|---|---|---|
0.01 | 1.20 × 10−06 | 4 | −2.08 × 10−01 | 1.20 × 10−06 | 0.2 | |
0.01 | 1.20 × 10−05 | 4 | −2.08 × 10−01 | 1.20 × 10−05 | 0.2 | |
0.01 | 1.20 × 10−04 | 4 | −2.08 × 10−01 | 1.20 × 10−04 | 0.2 | |
0.01 | 1.20 × 10−03 | 4 | −2.08 × 10−01 | 1.20 × 10−03 | 0.2 | |
0.1 | 1.20 × 10−05 | 4 | −1.18 × 10−01 | 1.20 × 10−05 | 0.2 | |
0.1 | 1.20 × 10−04 | 4 | −1.18 × 10−01 | 1.20 × 10−04 | 0.2 | |
0.1 | 1.20 × 10−03 | 4 | −1.18 × 10−01 | 1.20 × 10−03 | 0.2 | |
0.1 | 1.20 × 10−02 | 4 | −1.18 × 10−01 | 1.20 × 10−02 | 0.2 | |
1 | 1.20 × 10−04 | 2 | 9.98 × 10−01 | 1.20 × 10−04 | 0.0 | |
1 | 1.20 × 10−03 | 2 | 9.98 × 10−01 | 1.20 × 10−03 | 0.0 | |
1 | 1.20 × 10−02 | 2 | 9.98 × 10−01 | 1.20 × 10−02 | 0.0 | |
1 | 1.20 × 10−01 | 2 | 9.99 × 10−01 | 1.20 × 10−01 | 0.0 | |
10 | 1.20 × 10−03 | 2 | 9.98 × 10+00 | 1.20 × 10−03 | ||
10 | 1.20 × 10−02 | 2 | 9.98 × 10+00 | 1.20 × 10−02 | ||
10 | 1.20 × 10−01 | 2 | 9.98 × 10+00 | 1.20 × 10−01 | ||
10 | 1.20 × 10+00 | 4 | 1.01 × 10+01 | 1.21 × 10+00 | 0.9 | 0.5 |
100 | 1.20 × 10−02 | 5 | 9.87 × 10+01 | 1.17 × 10−02 | ||
100 | 1.20 × 10−01 | 5 | 9.87 × 10+01 | 1.17 × 10−01 | ||
100 | 1.20 × 10+00 | 7 | 1.06 × 10+02 | 8.75 × 10−01 | 6.3 | |
100 | 1.20 × 10+01 | 6 | 4.51 × 10+01 | 1.17 × 10+01 | ||
1000 | 1.20 × 10−01 | 4 | 5.01 × 10+02 | 3.02 × 10−02 | ||
1000 | 1.20 × 10+00 | 4 | 5.01 × 10+02 | 3.02 × 10−01 | ||
1000 | 1.20 × 10+01 | 4 | 5.22 × 10+02 | 3.00 × 10+00 | ||
1000 | 1.20 × 10+02 | 8 | 2.79 × 10+03 | 1.05 × 10+02 | 178.8 |
α | Mode | %D α | %D | |||
---|---|---|---|---|---|---|
0.01 | 1.20 × 10−06 | 4 | −2.08 × 10−01 | 1.20 × 10−06 | 0.2 | |
0.01 | 1.20 × 10−05 | 4 | −2.08 × 10−01 | 1.20 × 10−05 | 0.2 | |
0.01 | 1.20 × 10−04 | 4 | −2.08 × 10−01 | 1.20 × 10−04 | 0.2 | |
0.01 | 1.20 × 10−03 | 4 | −2.08 × 10−01 | 1.20 × 10−03 | 0.2 | |
0.1 | 1.20 × 10−05 | 4 | −1.18 × 10−01 | 1.20 × 10−05 | 0.2 | |
0.1 | 1.20 × 10−04 | 4 | −1.18 × 10−01 | 1.20 × 10−04 | 0.2 | |
0.1 | 1.20 × 10−03 | 4 | −1.18 × 10−01 | 1.20 × 10−03 | 0.2 | |
0.1 | 1.20 × 10−02 | 4 | −1.18 × 10−01 | 1.20 × 10−02 | 0.2 | |
1 | 1.20 × 10−04 | 2 | 9.98 × 10−01 | 1.20 × 10−04 | 0.0 | |
1 | 1.20 × 10−03 | 2 | 9.98 × 10−01 | 1.20 × 10−03 | 0.0 | |
1 | 1.20 × 10−02 | 2 | 9.98 × 10−01 | 1.20 × 10−02 | 0.0 | |
1 | 1.20 × 10−01 | 2 | 9.99 × 10−01 | 1.20 × 10−01 | 0.0 | |
10 | 1.20 × 10−03 | 2 | 9.98 × 10+00 | 1.20 × 10−03 | ||
10 | 1.20 × 10−02 | 2 | 9.98 × 10+00 | 1.20 × 10−02 | ||
10 | 1.20 × 10−01 | 2 | 9.98 × 10+00 | 1.20 × 10−01 | ||
10 | 1.20 × 10+00 | 4 | 1.01 × 10+01 | 1.21 × 10+00 | 0.9 | 0.5 |
100 | 1.20 × 10−02 | 5 | 9.87 × 10+01 | 1.17 × 10−02 | ||
100 | 1.20 × 10−01 | 5 | 9.87 × 10+01 | 1.17 × 10−01 | ||
100 | 1.20 × 10+00 | 7 | 1.06 × 10+02 | 8.75 × 10−01 | 6.3 | |
100 | 1.20 × 10+01 | 6 | 4.51 × 10+01 | 1.17 × 10+01 | ||
1000 | 1.20 × 10−01 | 4 | 5.01 × 10+02 | 3.02 × 10−02 | ||
1000 | 1.20 × 10+00 | 4 | 5.01 × 10+02 | 3.02 × 10−01 | ||
1000 | 1.20 × 10+01 | 4 | 5.22 × 10+02 | 3.00 × 10+00 | ||
1000 | 1.20 × 10+02 | 8 | 2.79 × 10+03 | 1.05 × 10+02 | 178.8 |
Table 4 presents the prediction results for α and using the mode with maximum sensitivity using . We see that, by using this criterion, the predictions are very similar to those using Eq. (18), except for cases of assigned values of α ≤ 0.01. As before, the modes selected by this criterion differ from the two previous criteria used. On the basis of these comparisons, thus, we recommend that the sensitivity in Eq. (18) be used to isolate the “best” prediction of the α and values from an array of available modal data.
Predicted and versus assigned values and , with percent differences (%D), using the eigenmode (indicated) with maximum sensitivity
α | Mode | %D α | %D | |||
---|---|---|---|---|---|---|
0.01 | 1.20 × 10−06 | 2 | 8.72 × 10−03 | 1.20 × 10−06 | −12.8 | 0.0 |
0.01 | 1.20 × 10−05 | 2 | 8.72 × 10−03 | 1.20 × 10−05 | −12.8 | 0.0 |
0.01 | 1.20 × 10−04 | 2 | 8.72 × 10−03 | 1.20 × 10−04 | −12.8 | 0.0 |
0.01 | 1.20 × 10−03 | 2 | 8.73 × 10−03 | 1.20 × 10−03 | −12.7 | 0.0 |
0.1 | 1.20 × 10−05 | 1 | 9.97 × 10−02 | 1.20 × 10−05 | −0.3 | 0.0 |
0.1 | 1.20 × 10−04 | 1 | 9.97 × 10−02 | 1.20 × 10−04 | −0.3 | 0.0 |
0.1 | 1.20 × 10−03 | 1 | 9.97 × 10−02 | 1.20 × 10−03 | −0.3 | 0.0 |
0.1 | 1.20 × 10−02 | 2 | 9.87 × 10−02 | 1.20 × 10−02 | −1.3 | 0.0 |
1 | 1.20 × 10−04 | 2 | 9.98 × 10−01 | 1.20 × 10−04 | −0.2 | 0.0 |
1 | 1.20 × 10−03 | 2 | 9.98 × 10−01 | 1.20 × 10−03 | −0.2 | 0.0 |
1 | 1.20 × 10−02 | 2 | 9.98 × 10−01 | 1.20 × 10−02 | −0.2 | 0.0 |
1 | 1.20 × 10−01 | 2 | 9.99 × 10−01 | 1.20 × 10−01 | −0.1 | 0.0 |
10 | 1.20 × 10−03 | 2 | 9.98 × 10+00 | 1.20 × 10−03 | −0.2 | −0.4 |
10 | 1.20 × 10−02 | 2 | 9.98 × 10+00 | 1.20 × 10−02 | −0.2 | −0.4 |
10 | 1.20 × 10−01 | 2 | 9.98 × 10+00 | 1.20 × 10−01 | −0.2 | −0.4 |
10 | 1.20 × 10+00 | 4 | 1.01 × 10+01 | 1.21 × 10+00 | 0.9 | 0.5 |
100 | 1.20 × 10−02 | 5 | 9.87 × 10+01 | 1.17 × 10−02 | −1.3 | −2.6 |
100 | 1.20 × 10−01 | 5 | 9.87 × 10+01 | 1.17 × 10−01 | −1.3 | −2.6 |
100 | 1.20 × 10+00 | 5 | 9.93 × 10+01 | 1.17 × 10+00 | −0.7 | −2.5 |
100 | 1.20 × 10+01 | 2 | 1.08 × 10+02 | 1.16 × 10+01 | 7.5 | −3.6 |
1000 | 1.20 × 10−01 | 1 | 7.48 × 10+02 | 6.71 × 10−02 | −25.2 | −44.1 |
1000 | 1.20 × 10+00 | 1 | 7.48 × 10+02 | 6.71 × 10−01 | −25.2 | −44.1 |
1000 | 1.20 × 10+01 | 1 | 7.46 × 10+02 | 6.66 × 10+00 | −25.4 | −44.5 |
1000 | 1.20 × 10+02 | 2 | 1.40 × 10+03 | 7.61 × 10+01 | 40.0 | −36.6 |
α | Mode | %D α | %D | |||
---|---|---|---|---|---|---|
0.01 | 1.20 × 10−06 | 2 | 8.72 × 10−03 | 1.20 × 10−06 | −12.8 | 0.0 |
0.01 | 1.20 × 10−05 | 2 | 8.72 × 10−03 | 1.20 × 10−05 | −12.8 | 0.0 |
0.01 | 1.20 × 10−04 | 2 | 8.72 × 10−03 | 1.20 × 10−04 | −12.8 | 0.0 |
0.01 | 1.20 × 10−03 | 2 | 8.73 × 10−03 | 1.20 × 10−03 | −12.7 | 0.0 |
0.1 | 1.20 × 10−05 | 1 | 9.97 × 10−02 | 1.20 × 10−05 | −0.3 | 0.0 |
0.1 | 1.20 × 10−04 | 1 | 9.97 × 10−02 | 1.20 × 10−04 | −0.3 | 0.0 |
0.1 | 1.20 × 10−03 | 1 | 9.97 × 10−02 | 1.20 × 10−03 | −0.3 | 0.0 |
0.1 | 1.20 × 10−02 | 2 | 9.87 × 10−02 | 1.20 × 10−02 | −1.3 | 0.0 |
1 | 1.20 × 10−04 | 2 | 9.98 × 10−01 | 1.20 × 10−04 | −0.2 | 0.0 |
1 | 1.20 × 10−03 | 2 | 9.98 × 10−01 | 1.20 × 10−03 | −0.2 | 0.0 |
1 | 1.20 × 10−02 | 2 | 9.98 × 10−01 | 1.20 × 10−02 | −0.2 | 0.0 |
1 | 1.20 × 10−01 | 2 | 9.99 × 10−01 | 1.20 × 10−01 | −0.1 | 0.0 |
10 | 1.20 × 10−03 | 2 | 9.98 × 10+00 | 1.20 × 10−03 | −0.2 | −0.4 |
10 | 1.20 × 10−02 | 2 | 9.98 × 10+00 | 1.20 × 10−02 | −0.2 | −0.4 |
10 | 1.20 × 10−01 | 2 | 9.98 × 10+00 | 1.20 × 10−01 | −0.2 | −0.4 |
10 | 1.20 × 10+00 | 4 | 1.01 × 10+01 | 1.21 × 10+00 | 0.9 | 0.5 |
100 | 1.20 × 10−02 | 5 | 9.87 × 10+01 | 1.17 × 10−02 | −1.3 | −2.6 |
100 | 1.20 × 10−01 | 5 | 9.87 × 10+01 | 1.17 × 10−01 | −1.3 | −2.6 |
100 | 1.20 × 10+00 | 5 | 9.93 × 10+01 | 1.17 × 10+00 | −0.7 | −2.5 |
100 | 1.20 × 10+01 | 2 | 1.08 × 10+02 | 1.16 × 10+01 | 7.5 | −3.6 |
1000 | 1.20 × 10−01 | 1 | 7.48 × 10+02 | 6.71 × 10−02 | −25.2 | −44.1 |
1000 | 1.20 × 10+00 | 1 | 7.48 × 10+02 | 6.71 × 10−01 | −25.2 | −44.1 |
1000 | 1.20 × 10+01 | 1 | 7.46 × 10+02 | 6.66 × 10+00 | −25.4 | −44.5 |
1000 | 1.20 × 10+02 | 2 | 1.40 × 10+03 | 7.61 × 10+01 | 40.0 | −36.6 |
5 Conclusions
In this study, we have presented a new method, to be used in contact resonance atomic force microscopy, to accurately determine the viscoelastic properties of a sample using a plate-like microsensor. The theoretical formulation, based on the Rayleigh–Ritz method and on an “unconventional” generalized eigenvalue problem, was presented and a numerical experiment was devised to validate the method capability to identify viscoelastic parameters from modal data. The proposed method performed extremely well for nondimensional sample stiffnesses in the range 0.01 ≤ α ≤ 100 and nondimensional damping values in the range . For values of α ≥ 1000, the prediction error became increasingly larger with successively larger values of .
This study opens up several avenues for future research. First, future work should include validation of the method against experimental results. While we anticipate that the method will perform well in a practical implementation, additional study must be performed to determine the validity of the sensitivity metric used in Eq. (18) and whether better metrics exist to accurately select the modal data used in the predictions, given experimental test data rather than numerical surrogates. In addition, in the future work, we would like to address: the optimal placement of the sensing tip given damping, the effect of large material damping, and the implications of fluid–structure interactions in liquid environments and their effects on the sensor estimates.
Acknowledgment
This material is based in part upon work supported by the National Science Foundation under grant CMMI-1660448 (to R. T.) and grant CMMI-1847513 (to M. A.).
Conflict of Interest
There are no conflicts of interest.