Research Article  Open Access
NghiaHuu Nguyen, MooYeon Lee, JiSoo Kim, DongYeon Lee, "Compliance Matrix of a SingleBent Leaf Flexure for a Modal Analysis", Shock and Vibration, vol. 2015, Article ID 672831, 10 pages, 2015. https://doi.org/10.1155/2015/672831
Compliance Matrix of a SingleBent Leaf Flexure for a Modal Analysis
Abstract
We present a compliance matrix for a singlebent leaf flexure (SBLF) that shows the relationships between the deformations and the sixaxis loads applied to the SBLF. Higherorder beam theory that considers the variable shear and warping effect is considered in bending. The partially restrained warping at the junction between elements is also considered in torsion. The total strain energy is calculated, and the complete compliance matrix is derived by using Castigliano’s second theorem. Sensitivity analyses over the compliance elements are performed and verified via finite element analysis (FEA). The results show that the derived compliance elements are in good agreement with FEA, with errors of less than 7.6%. We suggest that theoretical compliance elements considering variable shear and warping in bending and partially restrained warping in torsion give highly accurate design equations representing the compliant mechanism of the SBLF. The present work could be used in a modal analysis of a singlebent leaf flexure.
1. Introduction
Flexure guides are used in precision engineering, especially in nanoscanner devices, because of their smooth, elastic, and nofriction characteristics. Flexure guides provide short or moderate ranges of travel due to their elastic deformation characteristics [1]. New flexure guide which has large working range and nanoresolution is hot topic in the recent years. Several types of flexure guides have been described in the literature including hinge flexures [2–7] and leafspring flexures [8–13]. A singlebent leaf flexure, investigated in this study, is an Lshaped leafspring mechanism. The SLBF has remarkable advantages for flexure such as a moderate range of deflection and easily obtained uniform spring material. Many studies have focused on flexure guides. Formulas are derived for the stiffness of a parallel leafspring flexure in [8]. A displacement reduction mechanism based on torsion leafspring hinges was developed in [9]. An orthoplanar spring design that operates by raising or lowering its platform relative to the base with no rotation is presented in [10]. A planar positioning mechanism with three degrees of freedom using a flexure guide is presented in [11]. New flexures were designed and developed with a large range and compact dimensions, as described in [1, 12, 14].
The compliance matrix is often used to express the loaddeformation relationship for a linearly elastic system. It has been applied in many previous studies in the precision field [3, 4, 15–17] to derive the equations of motion for general flexure mechanisms. The compliance matrix for a prismatic beam is presented in [3], but the shear deformation and the warping effect in bending and torsion were not considered. The constant shear deformation is considered in bending, the warping is not included in torsion [4, 15]. The warping was considered in the torsion theory; however, it was applied to the beams with a fully fixed (restrained) end [18, 19]. The partially restrained warping in torsion was investigated in [20] and was applied to truck chassis frames. The general higherorder beam theory (HBT) was presented in [21–23], and Levinson [21] proposed an equation for the deflection of a beam that includes the shear and warping in bending analysis. Relationships between Levinson beam theory and classical EulerBernoulli theory are presented and applied to examples by Reddy et al. [22]. However, these models are usually applied to the cantilever beam or the simply supported beam models, not to the SBLF. The compliance matrix of a bent beam microsuspension was presented in [24]. The beam model is similar to the SBLF model; however, only the inplane deformations were derived and the shear deformations were not considered in the bending analysis.
In this study, we analyzed and derived a compliance matrix that expresses the relationship between the deformations and applied loads of the SBLF by using Castigliano’s second theorem. The HBT of Levinson is applied in the bending analysis. The partially restrained warping at the junction between elements is introduced, and the full and free warping is also investigated in the torsional analysis. The theoretical results are verified by finite element analysis (FEA) results. The good agreement between the results of the two methods validates the accuracy of the theoretical equations.
2. Generalized Modeling of the SBLF
Figure 1(a) shows a model of the SBLF that consists of two leaf flexure elements. One end is fixed, and other end has six components of loading: three forces (, , and ) and three moments (, , and ). Dimensions , and are the length, width, and thickness, respectively. Under the 6axis loads, axial, bending, shear, and torsional deformations will occur in both elements of the SBLF. Figure 1(b) shows one example of the SBLF in a planar nanoscanner, which includes four SBLFs symmetrically connected with the four corners of the square moving body to ensure smooth and nonparasitic error motion of the system. In this study, the deformations due to the six components of loading need to be analyzed and determined in the compliance matrix form.
(a)
(b)
2.1. Derivation of the Total Strain Energy
Castigliano’s second theorem was used to find the translational and rotational deformations of the SBLF. Thus, the deformations due to loads are defined by the partial derivative of the total strain energy (SE) with respect to the loads as follows:where is the total strain energy that is stored in the flexure, is the deflection due to the force , and is the rotation due to the moment . When the six components of loading are applied, the total SE stored in two elements of the SBLF can be classified as the extensional, shear, bending, and torsional deformation energies.
The extensional SE () in element 1 due to the axial force areIn element 2, the extensional SE () due to the axial force is given by where is the crosssectional area, and is the elastic modulus (Young’s modulus).
In the bending and shear analyses, the shear deformations due to the transverse loads are usually neglected in EulerBernoulli beam theory (EBT). Accordingly, the EBT can predict the deflection accurately only for a thin beam (ratio ). Timoshenko beam theory considers the additional constant shear deformation. Thus, TBT shows more precise analysis results than EBT. However, for thick beams (), the higherorder beam theory that considers the variable shear deformation and the warping in bending should be applied. In the HBT of Levinson, the variable shear deformation and warping in bending are considered and the deflection is described in [21]. Therefore, the shear SE () in element 1 due to the shear forces , is given bywhere is Poisson’s ratio. In element 2, the shear SE () due to the shear forces , is given byThe bending SE () in element 1 due to , , , and is given by In element 2, the bending SE () due to , , , , and is given by where and is the inertia moment about the , and axes, respectively.
A partially restrained warping torsion analysis at the joint between elements 1 and 2 is presented in [25]. The warping restraint factor () is adopted to represent the degree of warping restraint. In the present study, the torsional SE () was defined in element 1 due to moment as follows [25]:where In element 2, the torsional SE () due to moments and is as follows:wherewhere .
The torsion constant of a bar with a rectangular cross section can be determined by Equation in [26]. Therefore, Warping constant with respect to the shear center is defined by Equation in [27] as the warping moment of inertia: where is the warping function and is the torsionbending constant ().
, the warping restraint factor, is introduced to express the degree of restraint from the warping at the joint of elements 1 and 2. Equations (9) and (11) show that when , free warping occurs at end or (), and fully restrained warping occurs when .
The total SE () of a SBLF with sixaxis loads at free end was summarized as follows:where is the total SE of element 1, ; and is the total SE of element 2, .
From (1), the translational displacements due to forces , , and according to the , and axes, respectively, are given byThe rotational displacements due to moments , , and about the , , and axes, respectively, are given bywhere
2.2. Compliance Matrix of the SBLF
The relationship between the translational and rotational displacements and the applied forces and moments at the free end of the SBLF is formulated as follows:where , , [] is the compliance matrix or flexibility matrix, [] is the stiffness matrix, and are the force and translational displacement with respect to the axis, and and are the moment and rotational displacement about the axis, respectively. The compliance matrix [] includes 11 independent elements as follows:
3. FEA Verification
In this study, FEA was conducted by using ProMechanica commercial software (Wildfire 5, PTC Corp., MA, USA) to verify the results obtained by theory and to validate the reliability of the calculated results. The default values of the SBLF were chosen as follows: length mm, width mm, and thickness mm. The forces and moments applied to the SBLF are N, N, and Nmm. The warping restraint factor was chosen to be . The material selected in this simulation is aluminum Al 6061. In the parametric analysis, the sensitivity analysis was performed to determine the agreement between FEA and theory with the variation of length to 20 mm, width to 8 mm, and thickness to 1 mm, respectively. If the errors between the two methods are lower than 10%, the results can usually be accepted and used from the viewpoint of the engineering design.
Table 1 shows the comparison between the theory and FEA results at the default values of the SBLF. Equations (15)–(16) were applied to find the translational and rotational displacements by the theoretical method that considers the variable shear deformation in the bending and the warping in both the bending and torsion and the partially restrained warping at the joint of elements 1 and 2. The derived theoretical equations showed a good agreement with the FEA results; all of the errors of these simulations were less than 5%, whereas the results of without shear and warp effect showed higher errors, up to 7.4% as presented in Table 1. The results of inplane deformations , , and in this study were consistent with results of the equations in [24]. Although these results were showing the good agreement at the default values, the sensitive parameters analysis was also conducted to ensure the reliability of the derived equations.

In the research of the reference [28], we analyzed the displacement of a SBLF under transverse loading by the EulerBernoulli, Timoshenko, and thirdorder beam theory and the results showed that the thirdorder beam theory of Levinson which considers the shear and warp effect in bending gave the good results. In this study, the above beam theory was also used and the translational displacements , , and under , , and according to the variations of length , thickness , and width , respectively, are presented in Table 2. Figure 2 shows typical simulation results under force . The errors in these investigations are lower (7.6%). The sensitive analysis of is also presented in Figure 3. All the errors are lower (5%). These results reveal the accuracy of the derived theoretical equations for the translational displacements of SBLF.

(a)
(b)
(c)
(a)
(b)
(c)
The comparison results among the fully restrained, partially restrained, and free warping in torsion analysis of SBLF were investigated in [25] and the partially restrained warping with factor was applied in this torsion simulation. Table 3 shows the variation of rotational displacements , , and due to moments , , and with the variation of parameters , , and of the SBLF, respectively. It can be observed that the results of theory and FEA are in good agreement, with the errors lower than 6%. Figure 4 presents the typical variation of rotation according to , , and under . The maximum error is lower than 6%.

(a)
(b)
(c)
In summary, the relationship between the deformation and the applied loads for the SBLF is presented in compliance matrix form. The variable shear deformation and warping effect were considered, and partial restraint was introduced in the analysis. All results were verified by FEA, with strong agreement between the two methods. The errors were lower than 7.6%. These results demonstrate the high accuracy and reliability of the proposed theoretical equations.
4. Conclusion
In this study, a compliance matrix that expresses the relationship between the deformations and applied loads of the SBLF was analyzed and derived by using Castigliano’s second theorem. In bending analysis, higherorder beam theory was applied wherein variable shear deformation and the warping were considered in the calculated formulas for the shear forces. At the joint of two elements of the SBLF, the partially restrained warping in torsion was analysed with consideration of the warping restraint factor. The theoretical results were verified by FEA at both the default and the sensitive values. The results indicate that there is strong agreement between the two methods, with errors below 7.6%. This suggests the accuracy of the proposed theoretical equations and that they can be used in the precision machine design.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Authors’ Contribution
NghiaHuu Nguyen and MooYeon Lee equally contributed to this work.
Acknowledgment
This research was supported in part by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF2013R1A1A4A01009657).
References
 E.J. Park, J. Shim, D.Y. Lee, and J. Lee, “A doublebent planar leaf flexure guide for a nanoscanner,” Journal of the Korean Physical Society, vol. 57, no. 61, pp. 1581–1588, 2010. View at: Publisher Site  Google Scholar
 K. Kim, D. Ahn, and D. Gweon, “Optimal design of a 1rotational DOF flexure joint for a 3DOF Htype stage,” Mechatronics, vol. 22, no. 1, pp. 24–32, 2011. View at: Publisher Site  Google Scholar
 T. T. Y. Koseki, T. Tanikawa, N. Koyachi, and T. Arai, “Kinematic analysis of a translational 3d.o.f. microparallel mechanism using the matrix method,” Advanced Robotics, vol. 16, no. 3, pp. 251–264, 2002. View at: Publisher Site  Google Scholar
 J. H. Kim, S. H. Kim, and Y. K. Kwak, “Development and optimization of 3D bridgetype hinge mechanisms,” Sensors and Actuators A: Physical, vol. 116, no. 3, pp. 530–538, 2004. View at: Publisher Site  Google Scholar
 J. W. Ryu, D.G. Gweon, and K. S. Moon, “Optimal design of a flexure hinge based XYφ wafer stage,” Precision Engineering, vol. 21, no. 1, pp. 18–28, 1997. View at: Publisher Site  Google Scholar
 N. Lobontiu and E. Garcia, “Analytical model of displacement amplification and stiffness optimization for a class of flexurebased compliant mechanisms,” Computers & Structures, vol. 81, no. 32, pp. 2797–2810, 2003. View at: Publisher Site  Google Scholar
 Y. Li and Q. Xu, “A novel piezoactuated XY stage with parallel, decoupled, and stacked flexure structure for micro/nanopositioning,” IEEE Transactions on Industrial Electronics, vol. 58, no. 8, pp. 3601–3615, 2011. View at: Publisher Site  Google Scholar
 D. M. Brouwer, J. P. Meijaard, and J. B. Jonker, “Large deflection stiffness analysis of parallel prismatic leafspring flexures,” Precision Engineering, vol. 37, no. 3, pp. 505–521, 2013. View at: Publisher Site  Google Scholar
 M. Hayashi and M. Fukuda, “Generation of nanometer displacement using reduction mechanism consisting of torsional leaf spring hinges,” International Journal of Precision Engineering and Manufacturing, vol. 13, no. 5, pp. 679–684, 2012. View at: Publisher Site  Google Scholar
 J. J. Parise, L. L. Howell, and S. P. Magleby, “Orthoplanar linearmotion springs,” Mechanism and Machine Theory, vol. 36, no. 1112, pp. 1281–1299, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Fukada and K. Nishimura, “Nanometric positioning over a onemillimeter stroke using a flexure guide and electromagnetic linear motor,” International Journal of Precision Engineering and Manufacturing, vol. 8, pp. 49–53, 2007. View at: Google Scholar
 Q. S. Xu, “Design and development of a compact flexurebased XY precision positioning system with centimeter range,” IEEE Transactions on Industrial Electronics, vol. 61, no. 2, pp. 893–903, 2014. View at: Publisher Site  Google Scholar
 X.P. S. Su and H. S. Yang, “Design of compliant microleverage mechanisms,” Sensors and Actuators, A: Physical, vol. 87, no. 3, pp. 146–156, 2001. View at: Publisher Site  Google Scholar
 Y. K. Yong, S. S. Aphale, and S. O. R. Moheimani, “Design, identification, and control of a flexurebased XY stage for fast nanoscale positioning,” IEEE Transactions on Nanotechnology, vol. 8, no. 1, pp. 46–54, 2009. View at: Publisher Site  Google Scholar
 J.J. Kim, Y.M. Choi, D. Ahn, B. Hwang, D.G. Gweon, and J. Jeong, “A millimeterrange flexurebased nanopositioning stage using a selfguided displacement amplification mechanism,” Mechanism and Machine Theory, vol. 50, pp. 109–120, 2012. View at: Publisher Site  Google Scholar
 S. Xiao, Y. Li, and Q. Meng, “Mobility analysis of a 3PUU flexurebased manipulator based on screw theory and compliance matrix method,” International Journal of Precision Engineering and Manufacturing, vol. 14, no. 8, pp. 1345–1353, 2013. View at: Publisher Site  Google Scholar
 N. Lobontiu, “Compliancebased matrix method for modeling the quasistatic response of planar serial flexurehinge mechanisms,” Precision Engineering, vol. 38, no. 3, pp. 639–650, 2014. View at: Publisher Site  Google Scholar
 M. Kujawa, “Torsion of restrained thinwalled bars of open constraint bisymmetric crosssection,” Tast Quarterly, vol. 16, pp. 5–15, 2011. View at: Google Scholar
 E. J. Sapountzakis, “Bars under torsional loading: a generalized beam theory approach,” ISRN Civil Engineering, vol. 2013, Article ID 916581, 39 pages, 2013. View at: Publisher Site  Google Scholar
 A. H. AlHaKeem, Structural analysis of truck chassis frames under longitudinal loads considering bimoment effects [Ph.D. thesis], Cranfield Institute of Technology, Cranfield, UK, 1991.
 M. Levinson, “A new rectangular beam theory,” Journal of Sound and Vibration, vol. 74, no. 1, pp. 81–87, 1981. View at: Publisher Site  Google Scholar
 J. N. Reddy, C. M. Wang, G. T. Lim, and K. H. Ng, “Bending solutions of Levinson beams and plates in terms of the classical Theories,” International Journal of Solids and Structures, vol. 38, no. 2627, pp. 4701–4720, 2001. View at: Publisher Site  Google Scholar
 C. M. Wang, J. N. Reddy, and K. H. Lee, Shear Deformable Beams and Plates: Relationships with Classical Solutions, Elsevier Science, Oxford, UK, 2000. View at: MathSciNet
 N. Lobontiu and E. Garcia, Mechanics of Microelectromechanical Systems, Kluwer Academic, New York, NY, USA, 2005.
 N. H. Nguyen, B. D. Lim, and D. Y. Lee, “Torsional analysis of a singlebent leaf flexure,” Structural Engineering and Mechanics, vol. 54, no. 1, pp. 189–198, 2015. View at: Publisher Site  Google Scholar
 S. P. Timoshenko and J. N. Goodier, Theory of Elasticity, McGrawHill, 2nd edition, 1951. View at: MathSciNet
 W. D. Pilkey, Analysis and Design of Elastic Beams: Computational Methods, John Wiley & Sons, New York, NY, USA, 2002.
 N. H. Nguyen, B. D. Lim, and D. Y. Lee, “Displacement analysis of a singlebent leaf flexure under transverse load,” International Journal of Precision Engineering and Manufacturing, vol. 16, no. 4, pp. 749–754, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 NghiaHuu Nguyen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.