Abstract
Material length scale considerably affects the mechanical properties of microcantilever components. Recently, cantilever-plate-like structures have been commonly used, whereas the lack of studies on their size effects constrains the design, testing, and application of these structures. We have studied the size-dependent dynamic behavior of a cantilever plate based on a modified couple stress theory and the differential quadrature method in this note. The numerical solutions of microcantilever plate equation involving the size effect have been presented. We have also analyzed the bending and vibration of the microcantilever plates considering the size effect and discussed the dependence of the size effect on their geometric dimensions. The results have shown that (1) the mechanical characteristics of the cantilever plate show obvious size effects; as a result, the bending deflection of a microcantilever plate reduces whereas the natural frequency increases effectively and (2) for the plates with the same material, the size effect becomes more obvious when the plates are thinner.
The miniaturized cantilever device is one of the fastest developing and most efficient types of micro-nano devices; it is widely used in the chemistry, physical, and biological fields [1β3]. In these applications, the dimensions of the cantilever are typically on the order of microns or even sub-microns. And the performance of cantilever device depends on the accurate extract of static bending and resonant frequency. Using test measurements and theoretical studies, scholars have recently proved that the mechanical properties of miniaturized structures critically depend on material length-scale. The size effect shown in the mechanical properties considerably affects the static bending and resonant frequency shift of cantilever devices [4β6]. Therefore, studying the size-dependent behavior of cantilever devices is important.
Couple stress theory is an effective method for studying the size effect [7, 8]. By rotating the gradient, the theory introduces material length-scale parameters to the constitutive equation. Lam et al. [9] proposed a modified couple stress model in which only one material length parameter is needed to capture the size effect, and the model has attracted many researchers in the past years. Kong et al. [10] used this theory to study the size effect of the Euler-Bernoulli beam natural frequency. Ma et al. [11] then studied the Timoshenko size effect that accounts for axial deformation and the Poisson effect.
Most widely known miniaturized devices are typical cantilever plate structures. Compared with a cantilever beam, a cantilever plate with a large surface area has higher sensitivity. However, current studies focus only on micro-beam structures, and few studies on the size effect of cantilever plates have been reported. Tsiatas [12] developed a new Kirchhoff plate model for the static analysis of isotropic microplates with arbitrary shape based on a modified couple stress theory. Yin et al. [13] proposed a non-classical Kirchhoff plate model for the dynamic analysis of microscale. However, they have not studied the cantilever plate.
In this letter, we use the modified couple stress theory and the differential quadrature (DQ) method to study the size effect of a microcantilever plate. The characteristic length effect on static bending deflection and natural frequency of plate under uniform loads is analyzed. Furthermore, we discuss the dependence of the size effect on the geometric dimensions of the plate.
We study a rectangular Kirchhoff cantilever plate and assuming that the initial state of the microplate is flat and thin. The plate is subjected to a uniform distributed transverse load and hence undergoes a three-dimensional deformation state. The represents the deflection of the direction.
Based on the modified couple stress theory, the governing equation of the static bending of a microcantilever plate is [12] where is the classical bending rigidity of the plate and is the contribution of rotation gradients to the bending rigidity. The two parameters can be expressed as where is Youngβs modulus, is Poissonβs ratio, and denotes the material length-scale parameter representing the size effect, usually in the micro-nano order.
On the basis of (2), the governing equation for the vibration of the plate under uniform loads can be expressed as where is the plate density and is the thickness of microcantilever.
The equation based on the modified couple stress theory contains only one material length-scale parameter. Therefore, the size effect can be well characterized. Meanwhile, (2) shows that when is zero, the equation degenerates to the classical equations of the Kirchhoff plate. It should be noted, in this letter we focus on the size effect and therefore ignored surface effect and induce residual stress.
For solving the equation of microplate, the analytical method is well-known cumbersome, and with the change of boundary and loading conditions, the difficulty of solving a sharp increase even cannot solve. In commonly used finite difference method and finite element method, the solving regional needs to be divided into fairly large number of grid or node, leading to large amount of calculation. However, as an effective numerical approach to solving boundary value problems, the DQ method can get enough precision of the calculation results using a few grid points, and this is the most significant advantage in comparison to other methods. We obtain the numerical solutions of static bending and vibration equations using this method.
By setting , , where and are the length and width of the plate, respectively. The coordinate dimensionless form of (1) can be expressed as
In the DQ method, the value of function partial derivative at one point can be approximately to weighted summation of function value of all the nodes, and the -order derivative of function at node is expressed as [14]: where is the -order weighting coefficient and is the number of nodes.
By substituting (5) into (4), we can obtain the microplate static bending equation in DQ form: where , . and are the number of nodes along the and directions of the plate. When (6) is expressed in matrix form, the equation becomes where is the weight coefficient matrix, is the nodal displacement array, and is the load vector.
The boundary conditions of the cantilever plate are
We use the equation substitution method to processing boundary conditions. Based on this method, the DQ equations in (8) substituted the nodes equation on ; in (6). The correctional weight coefficient matrix is obtained. The boundary conditions are introduced into (7).
We can see that and are known matrixes in (7), and then we can obtain the nodal displacement array under uniform loads. Finally, we can obtain the full displacement fields using a higher-order Lagrange interpolation.
Similarly, by transforming and arranging (3), the cantilever vibration equation in DQ form becomes where is the mode function and is the angular frequency. The matrix form of (8) is
According to the approach in (7), the natural frequency of the cantilever vibration can be obtained. Because the uniform load does not affect the natural frequency of the plate, the obtained value of the natural frequency is the natural frequency of the free vibration of the plate.
We present the numerical solution of the mechanical equations of a cantilever plate by the DQ method, then we expect to bring value to further analyses of the size effects of cantilever plates.
We set both length and width as 50βΞΌm, thickness as 5βΞΌm, modulus of elasticity as 50βGpa, density as 2700βkg/m3, and uniform load as 1βN/ΞΌm. The nodes were chosen for the grid point in solving regional, and node distributions are taken based on well-accepted Gauss-Chebyshev-Lobatto points [14].
In order to verify the effectiveness of the DQ method, we present the comparison results of analytical and DQ solution for deflection on βΞΌm, as shown in Table 1. We can see that, compared with analytical solutions, the error of DQ solution is less than 1% at different position . It shows the satisfactory accuracy with a few grid points and full performances the characteristics of little computation quantity and high precision.
This letter is the first to present an investigation of the size effect of a cantilever plate by numerical simulation. The following issues are analyzed: (1) the changes in the cantilever deflection and natural frequency of the plate with certain geometric dimensions caused by the size effect and (2) the influence of the geometric dimensions of a plate on the strength of the size effect.
Figure 1 is obtained using (7), when Poissonβs ratio as 0.3, a thickness as 5βΞΌm, at different length-scale parameters, the bending deflection of the cantilever plate along the direction (βΞΌm). Here, has values of 0, 1.5, 3, and 4.5βΞΌm. Figure 1 shows that when the length-scale parameter increases, the deflection of the plate considerably decreases at the same location. For βΞΌm, when the length-scale parameter increases from 0 to 4.5βΞΌm, the bending deflection decreases from 0.15 to 0.03βΞΌmβa reduction of 80%. This reduction indicates an obvious size effect. We emphasize that when the length-scale parameter is zero (), the bending deflection is the classical elasticity solution. Therefore, the data show that, in considering the influence of the size effect, the deflection of the plate is substantially reduced.
Table 2 is obtained using (9), when having Poissonβs ratio as 0.3, a thickness as 5βΞΌm, at different length-scale parameters, the first four orders of the natural frequencies of the cantilever plate; has values of 0, 1.5, 3, and 4.5βΞΌm. Table 1 shows that when the length-scale parameters increase, the natural frequency of each order mode considerably increases. For the first-order mode, when the length-scale parameters increase from 0 to 4.5βΞΌm, the natural frequency increases from 14 to 30βMHzβan increase of 114%. Similarly, in considering the influence of the size effect, the natural frequency of the plate considerably increases.
Equation (2) can explain the phenomenon above, that is, when the length-scale parameters increases, the bending stiffness of the cantilever bending plate increases. Therefore, the plate bending deflection decreases and the natural frequency increases. In summary, when the length-scale parameters and plate thickness are in the same order of magnitude, the size effect cannot be disregarded. In addition, with the increase in length-scale parameters, the size effect becomes increasingly apparent. These results are in accordance with the size effect of microcantilever beams.
To study the dependence of the size effect on the geometric dimensions of the cantilever plate, we set the definitions of the relative change in deflection as and the relative change in frequency as , where and are respectively the bending deflection and natural frequency that do not consider the size effect (), and and are, respectively, the bending deflection and natural frequency that consider the size effect. By this definition, for the certain length-scale parameters, the influence of the size effect on mechanical properties is more substantial when relative change is greater.
The size effects of each point on the plate under the same conditions are the same. Thus, a fixed point (βΞΌm, βΞΌm) is chosen. Figure 2 and Table 3 show the same material (with βΞΌm), and the relative changes in deflection and natural frequency with thickness and Poissonβs ratios. The thickness is between 1 to 5βΞΌm, and Poissonβs ratios are 0.25, 0.3, and 0.35.
We can see that, the size effect on mechanical properties weakens with the increase in thickness. For a Poissonβs ratio of 0.3, when the thickness of the cantilever plate increases from 1 to 5βΞΌm, the relative change in deflection is reduced from 0.98 to 0.79 and the relative natural frequency decreases from 9.44 to 1.28. These results indicate a considerably reduced size effect. In addition, for the cantilever plate of a certain thickness, the relative changes in deflection and natural frequency increase with the decrease in Poissonβs ratio.
Based on the modified couple stress theory and the differential quadrature (DQ) method, we studied the size-dependent dynamic behavior of a rectangular cantilever plate. We provide numerical solutions to the static bending deflection and natural frequency and analyzed static bending deflection and the natural frequency vibration under the size effect and discussed the dependence of the size effect on the geometric dimensions of the plate. The results show that the DQ solution seems quite well to approach analytical solution and performances the satisfactory accuracy with a few grid points. When the length-scale parameters and the thickness are in the same magnitude, the size effect cannot be disregarded. Under the influence of the size effect, the bending deflection decreases and the natural frequency increases. For the cantilever plate with the same material, the size effect on the mechanical properties gradually increases when thickness decreases. This study provided data on the size effect of cantilever plates and serves as basis for future research.
Acknowledgment
This work is supported by the National Natural Science Foundation of China Project no. 10902031 and China Postdoctoral Science Foundation Project no. 20090451001.
Supplementary Materials
In the supplementary material we give a detailed derivation process for governing equation of micro cantilever plates based on the modified Cosserat theory. Meanwhile, we provide detailed solving steps of governing equation using the differential quadrature method (DQM).