Sage Journals: Discover world-class research

Abstract

This paper presents a new class of flexure hinges, namely, conic-V-shaped flexure hinges (CFHs), which can be used as a generalized model for flexure hinges with profiles such as parabolic-V-shape, elliptical-V-shape, and hyperbolic-V-shape. Compliance and precision equations for the CFHs were derived as a set of nonlinear equations using Castigliano’s second theorem. The parameters of the nonlinear equations inputted to the compliance and precision matrices were based on the generalized equations used for conic curves in polar coordinates. Furthermore, the compliance equations were verified by means of finite element analysis and experiments. The errors in the finite element and experimental results were within 10% and 8% compared to the analytical results, respectively. Finally, the effects of dimensional parameters on the analytical model could be effectively analyzed by numerical simulations and comparisons.

Keywords

Flexure hinge compliance finite element analysis numerical simulations

Introduction

Owing to small size, high sensitivity, monolithic manufacturing process, zero backlash and friction, lack of lubrication, smooth motion, and virtually-infinite resolution of flexure hinges, they are widely used in compliant mechanisms for micro/nano positioning, ultra-precise machining and manufacturing, and automated biological cell manipulation; other applications are in micro-positioning stages, scanning electron microscopy, X-ray lithography, and piezoelectric actuators.^1–6

Flexure hinges with various cutout profiles have been proposed and studied.^7–12 In 1965, Paros and Weisbord¹³ proposed circular flexure hinges and derived simplified design equations, including a complete theoretical analysis. Smith et al.¹⁴ derived a compliance matrix for corner-filleted and shallow-notch elliptical flexure hinges and constructed a compliance calculation model. Chen et al.¹⁵ introduced the centrifugal angle as the integral variable to derive a compliance equation for elliptical flexure hinges. Lobontiu et al.¹⁶ deduced compliance equations for parabolic and hyperbolic flexure hinges. Wang et al.¹⁷ proposed a novel class of exponent-sine-shaped flexure hinges with asymmetric structures and higher motion accuracy. Li et al.¹⁸ derived closed-form compliance equations for a power-function-shaped flexure hinge and investigated its motion accuracy. Lobontiu et al.¹⁹ presented an analytical approach to corner-filleted flexure hinges. Yong et al.²⁰ analyzed the comparison of various compliance/stiffness equations of circular flexure hinges with FEA results. Hsiao et al.²¹ propose a novel flexure hinge with 3° of freedom in which the motions are restricted in the same plane. Liu et al.²² designed and analyzed a new kind of flexure hinge obtained by using a topology optimization approach, namely, a quasi-V-shaped flexure hinge.

In the design of a flexure hinge, two main aspects need to be considered: the profile type and the parameter design.^23,24 A flexure hinge is typically designed by selecting a certain type of profile based on the designer’s experience. The corresponding design equations are used to design the parameters. The two aspects are theoretically independent of each other, and there is no clear basis for selecting the profile type. Hence, generalized models for flexure hinges with various cutout profiles have been proposed. For example, Vallance et al.²⁵ presented a unified geometric model for designing flexure hinges; the model is suitable for representing an entire class of conic-section hinges using quadratic rational Bézier curves. Chen et al.²⁶ incorporated elliptical arc, parabolic, and hyperbolic profiles in a generalized conic flexure hinge model and deduced all the elements in the compliance matrix. Tian et al.^27,28 presented filleted V-shaped flexure hinges (FVFHs), with profiles ranging from leaf-shaped to circular, and derived closed-form compliance equations and motion accuracy equations.

This paper proposes a new class of flexure hinges called conic-V-shaped flexure hinges (CFHs), in order to develop a generalized flexure hinge model and to improve the design of flexure hinges. The generalized flexure hinge model encompasses elliptical-V-shaped, parabolic-V-shaped, and hyperbolic-V-shaped. Figure 1 shows the 3D geometry of a CFH. The geometric parameters include the cutout length l, minimum thickness t, cross-sectional width w, cutout depth c, and the angle φ between the tangent to the V-shaped curve and the x axis. The V-shaped curves are tangent to the conic curves. Compliance and precision equations were deduced for this class of CFHs for the geometric optimization design in polar coordinates. Stress equations were also derived for the stress concentration. Finally, the finite element analysis and experiments were conducted to verify the compliance and precision equations.

Figure 1.

Schematic of a conic-V-shaped flexure hinge.

The generalized equation in polar coordinates

Figure 2 shows a conic curve along with its directrix and focus O_F. In polar coordinates, the conic curve with focus O_F as the origin can be expressed as

ρ (ϑ) = \frac{ep}{1 + e \cos ϑ}

(1)

where p (>0) is the distance between the corresponding directrix and the focus, and e (>0) is the eccentricity, which can be expressed as

e = \frac{ρ (ϑ)}{n (ϑ)}

(2)

When e > 1, the conic curve is a hyperbola; when e = 1, it is a parabola; and when e < 1, it is an ellipse.

Figure 2.

Conic curve represented in polar coordinates.

As shown in Figure 3(b), we assume that the entire notch of the CFH ranges from −θ _n = −π/2 to θ _n = π/2 (θ_n is the maximum polar angle of the notch). The conic curve ranges from −θ_m to θ_m (θ_m is the maximum polar angle of the conic curve). φ is the angle between the tangent to the V-shaped curve and the x axis. In this generalized model, e, p, φ, and t are the design parameters for the hinge cutouts. The equation for the tangents can be expressed as

r (ϑ) = \frac{ep (\cos θ_{m} + \sin θ_{m} \tan φ)}{(1 + ecos θ_{m}) (\cos θ + \sin θ \tan φ)}

(3)

where θ_m can be calculated as

ϑ_{m} = \arcsin (\frac{\sqrt{ta n^{2} φ (1 - e^{2}) + 1} + e}{ta n^{2} φ + 1} \tan φ)

(4)

The notch length can be calculated as

l = 2 ρ (\frac{π}{2}) = 2 ep

(5)

Figure 3 shows the relationships among three different types of flexure hinges. For different values of θ_m, φ, and e in a certain range, different cutout profiles of the flexure hinges can be obtained, such as elliptical-V-shape, parabolic-V-shape, hyperbolic-V-shape, leaf-shape, elliptical, parabolic, and hyperbolic. When φ = 0, the hinge is a leaf flexure hinge. When θ_m = π/2, the hinge is a conic flexure hinge (elliptical arc flexure hinge, parabolic flexure hinge, or hyperbolic flexure hinge). When π/2 > θ_m > 0, the hinge is a CFH (elliptical-V-shaped flexure hinge, parabolic-V-shaped flexure hinge, or hyperbolic-V-shaped flexure hinge). φ_m is the angle between the tangent and the x axis when θ_m = θ_n = π/2 which can be calculated as:

φ_{m} = \arctan \frac{1}{e}

(6)

Figure 3.

Relationships of three flexure hinges: (a) θ_m = 0, φ = 0, (b) 0 < θ_m < θ_n = π/2, 0 < φ < φ_m, and (c) θ_m = θ_n = π/2, φ = φ_m.

Compliance equations for conic-V-shaped flexure hinges

The notch area of the CFH is divided into countless cuboids vertically, and the height of a small bar dθ at position θ can be expressed as

h (ϑ) = {\begin{matrix} p (m - \frac{2 ecos ϑ}{1 + ecos ϑ}) = p g_{1} (ϑ) & - ϑ_{m} < ϑ < ϑ_{m} \\ p (n + \frac{2 esin (sign (ϑ) ϑ) \tan φ}{1 + ecos ϑ}) = p g_{2} (ϑ) & ϑ_{m} < | ϑ | < \frac{π}{2} \end{matrix}

(7)

where

m = \frac{2 e}{1 + e} + \frac{t}{p}

(8)

n = \frac{2 e}{1 + e} + \frac{t}{p} - \frac{2 e (ta n^{2} φ + 1)}{\sqrt{ta n^{2} φ (1 - e^{2}) + 1} + e}

(9)

As shown in Figure 3, the relationship between x and θ can be expressed as

x = \frac{l}{2} + ρ (ϑ) \sin (ϑ)

(10)

Differentiating the above equation yields:

dx = [ρ' (ϑ) \sin (ϑ) + ρ (ϑ) \cos (ϑ)] = \frac{ep (e + \cos ϑ)}{{(1 + ecos ϑ)}^{2}} d ϑ

(11)

Compliance matrix

When the minimum thickness of a flexure hinge t is assumed to be much lower than the notch length l, the hinge can be considered a fixed free beam with 6° of freedom subjected to bending, axial loading, shearing, and torsion, as shown in Figure 1. The force acting on the hinge can be defined by

F = {[F_{x}, F_{y}, F_{z}, M_{x}, M_{y}, M_{z}]}^{T}

(12)

and the deformations of the hinge can be defined by

X = {[Δ_{x}, Δ_{y}, Δ_{z}, α_{x}, α_{y}, α_{z}]}^{T}

(13)

Accordingly, the following equation can be obtained

X = C_{h} F

(14)

where C_h is the compliance matrix of the flexure hinge, which can be expressed as

C_{h} = {\begin{matrix} Δ_{x} / F_{x} & 0 & 0 & 0 & 0 & 0 \\ 0 & Δ_{y} / F_{y} & 0 & 0 & 0 & Δ_{y} / M_{z} \\ 0 & 0 & Δ_{z} / F_{z} & 0 & Δ_{z} / M_{y} & 0 \\ 0 & 0 & 0 & α_{x} / M_{x} & 0 & 0 \\ 0 & 0 & α_{y} / F_{z} & 0 & α_{y} / M_{y} & 0 \\ 0 & α_{z} / F_{y} & 0 & 0 & 0 & α_{z} / M_{z} \end{matrix}}

(15)

The compliance matrix is symmetric (α_z/F_y = Δ_y/M_z and α_y/F_z = Δ_z/M_y). Each factor in the compliance matrix is derived from the following beam theory.

Under an external load P, the elastic strain energy V_ε of a flexure hinge can be used to compute the corresponding deformation Δ based on Castiglione’s second theorem:

Δ = \frac{\partial V_{ε}}{\partial P}

(16)

The total strain energy is defined as follows:

V_{ε} = \frac{1}{2} \int_{- a}^{a} \frac{{F_{N}}^{2}}{EA} ds + \frac{1}{2} \int_{- a}^{a} \frac{M^{2}}{EI} ds + \frac{1}{2} \int_{- a}^{a} \frac{k {F_{S}}^{2}}{GA} ds

(17)

where E is the Young’s modulus, G is the shear modulus, M is the bending moment, F_N is the force along the x axis, F_S is the shearing force, I is the cross-sectional moment of area, and A is the cross-sectional area. The shearing coefficient k is equal to 6/5 for beams with rectangular cross sections.²⁹

Compliance equations

Angular compliance about the z axis

The z axis is the input axis of the flexure hinge, and the angular compliance corresponding to the torque M_z can be expressed as

\frac{α_{z}}{M_{z}} = \int_{- \frac{π}{2}}^{\frac{π}{2}} \frac{\frac{ep (e + \cos ϑ)}{{(1 + ecos ϑ)}^{2}}}{E I_{z} (ϑ)} d ϑ = \frac{24 e}{Ew p^{2}} I_{1}

(18)

where I_z(θ) represents the moment of inertia at the cross-sectional area about the neutral axis at position θ.

The effect of F_y on α_z is equal to its equivalent moment F_y[l/2 − ep sinθ/(1 + e cosθ)]. Therefore, the angular compliance corresponding to the force F_y can be calculated as

\frac{α_{z}}{F_{y}} = \int_{- \frac{π}{2}}^{\frac{π}{2}} \frac{(\frac{l}{2} - \frac{epsin ϑ}{1 + ecos ϑ}) \frac{ep (e + \cos ϑ)}{{(1 + ecos ϑ)}^{2}}}{E [w h^{3} (ϑ)] / 12} d θ = \frac{24 e^{2}}{Ewp} I_{1}

(19)

Angular compliance about the y axis

The angular compliance of the y axis corresponding to the torque M_y can be expressed as

\begin{matrix} \frac{α_{y}}{M_{y}} = \int_{- \frac{π}{2}}^{\frac{π}{2}} \frac{\frac{ep (\cos ϑ + e)}{{(1 + ecos ϑ)}^{2}}}{E I_{y} (ϑ)} d ϑ \\ = \int_{- \frac{π}{2}}^{- ϑ_{m}} \frac{\frac{ep (\cos ϑ + e)}{{(1 + ecos ϑ)}^{2}}}{E I_{y} (ϑ)} d ϑ + \int_{- ϑ_{m}}^{ϑ_{m}} \frac{\frac{ep (\cos ϑ + e)}{{(1 + ecos ϑ)}^{2}}}{E I_{y} (ϑ)} d ϑ + \int_{ϑ_{m}}^{\frac{π}{2}} \frac{\frac{ep (\cos ϑ + e)}{{(1 + ecos ϑ)}^{2}}}{E I_{y} (ϑ)} d ϑ \\ = \frac{24 e}{E w^{3}} I_{2} \end{matrix}

(20)

The effect of F_z on α_y is equal to its equivalent moment −F_z[l/2 − ep sinθ/(1 + e cosθ)]. Therefore, the angular compliance corresponding to the force F_z can be calculated as

\frac{α_{y}}{F_{z}} = \int_{- \frac{π}{2}}^{\frac{π}{2}} \frac{- (\frac{l}{2} - \frac{epsin ϑ}{1 + ecos ϑ}) \frac{ep (\cos ϑ + e)}{{(1 + ecos ϑ)}^{2}}}{E [w^{3} h (ϑ)] / 12} d ϑ = - \frac{24 e^{2} p}{E w^{3}} I_{2}

(21)

Linear compliance along the z axis

Since the compliance matrix is symmetric, the linear compliance of the torque M_y along the z axis can be expressed as

\frac{Δ z}{M_{y}} = \frac{α_{y}}{F_{z}} = - \frac{24 e^{2} p}{E w^{3}} I_{2}

(22)

The linear compliance of the force F_z along the z axis can be divided into two parts: a bending part (represented by Δ^b_z/F_z) and a shearing part (represented by $Δ_{z}^{s}$ /F_z). The shearing part can be calculated as

\frac{Δ_{z}^{s}}{F_{z}} = \int_{0}^{2 a} \frac{k}{Gwh (x)} dx = \frac{kep}{Gw} \int_{- \frac{π}{2}}^{\frac{π}{2}} \frac{e + \cos ϑ}{h (ϑ) {(1 + ecos ϑ)}^{2}} d ϑ = \frac{2 ke}{Gw} I_{2}

(23)

The bending part can be calculated as

\begin{matrix} \frac{Δ_{z}^{b}}{F_{z}} = - \int_{- \frac{π}{2}}^{\frac{π}{2}} [\int_{- \frac{π}{2}}^{ϑ} \frac{- (\frac{l}{2} - \frac{epsin γ}{1 + ecos γ}) \frac{ep (\cos γ + e)}{{(1 + ecos γ)}^{2}}}{E [w^{3} h (γ) / 12]} d γ] d (\frac{epsin ϑ}{1 + ecos ϑ}) \\ = \frac{24 e^{3} p^{2}}{E w^{3}} I_{2} + \frac{24 e^{3} p^{2}}{E w^{3}} I_{3} \end{matrix}

(24)

Therefore, the linear compliance along the z axis due to F_z can be expressed as

\frac{Δ_{z}}{F_{z}} = (\frac{2 ke}{Gw} + \frac{24 e^{3} p^{2}}{E w^{3}}) I_{2} + \frac{24 e^{3} p^{2}}{E w^{3}} I_{3}

(25)

Linear compliance along the y axis

Considering the symmetry in the compliance matrix, we can express the linear compliance of the torque M_z along the y axis as

\frac{Δ_{y}}{M_{z}} = \frac{α_{z}}{F_{y}} = \frac{24 e^{2}}{Ewp} I_{1}

(26)

Based on the derivation of Δ_z/M_y, the linear compliance along the y axis due to F_y can be expressed as

\frac{Δ_{y}}{F_{y}} = \frac{24 e^{3}}{E w^{3}} I_{1} + \frac{24 e^{3}}{Ew} I_{4} + \frac{2 ke}{Gw} I_{2}

(27)

Linear compliance along the x axis

The linear compliance along the x axis can be expressed as

\frac{Δ_{x}}{F_{x}} = \int_{- \frac{π}{2}}^{\frac{π}{2}} \frac{\frac{ep (e + \cos ϑ)}{{(1 + ecos ϑ)}^{2}}}{Ewh (ϑ)} d ϑ = \frac{2 e}{Ew} I_{2}

(28)

Angular compliance about the x axis

Since each infinitesimal element in the flexure hinge model can be considered a rectangular beam, according to the approximate torsion equation given in Hearn’s book,³⁰ the angular compliance about the x axis can be expressed as

\frac{α_{x}}{M_{x}} = \int_{- \frac{π}{2}}^{\frac{π}{2}} \frac{\frac{7 w^{3} h (ϑ) + 7 w h^{3} (ϑ)}{2} \frac{ep (e + \cos ϑ)}{{(1 + ecos ϑ)}^{2}}}{G w^{4} h^{4} (ϑ)} d ϑ = \frac{7 e}{Gw p^{2}} I_{1} + \frac{7 e}{G w^{3}} I_{2}

(29)

Precision equations for conic-V-shaped flexure hinges

To facilitate the analysis of the transmission precision of CFHs, the offset of the center point of the hinge profile (point S in Figure 3) is often used to characterize its transmission precision. By defining the offset of S as

O = {[ε_{x}, ε_{y}, ε_{z}]}^{T}

(30)

we can obtain the relationship between the offset and the force vector as

O = T_{h} F

(31)

where T_h is the precision matrix of the flexure hinge, which can be expressed as

T_{h} = {\begin{matrix} ε_{x} / F_{x} & 0 & 0 & 0 & 0 & 0 \\ 0 & ε_{y} / F_{y} & 0 & 0 & ε_{y} / M_{z} & 0 \\ 0 & 0 & ε_{z} / F_{z} & ε_{z} / M_{y} & 0 & 0 \end{matrix}}

(32)

The elements in the precision matrix are called precision factors.

In the precision matrix, the precision factors along the x axis can be given by

\frac{ε_{x}}{F_{x}} = \frac{e}{2 Ew} I_{2}

(33)

The precision factor corresponding to M_z along the y axis can be calculated as

\begin{matrix} \frac{ε_{y}}{M_{z}} = {[\frac{epsin ϑ}{1 + ecos ϑ} \int_{- \frac{π}{2}}^{ϑ} \frac{ep (\cos γ + e)}{Ew h^{3} (γ) {(1 + ecos γ)}^{2} / 12} d γ]}_{- \frac{π}{2}}^{0} - \int_{- \frac{π}{2}}^{0} \frac{12 e^{2} p^{2} \sin ϑ (\cos ϑ + e)}{Ew h^{3} (ϑ) {(1 + ecos ϑ)}^{3}} d ϑ \\ = \frac{12 e^{2}}{Ewp} I_{5} \end{matrix}

(34)

and the precision factor corresponding to the force F_y can be obtained using

\begin{matrix} \frac{ε_{y}}{F_{y}} = {[\frac{epsin ϑ}{1 + ecos ϑ} \int_{- \frac{π}{2}}^{ϑ} \frac{(\frac{l}{2} - \frac{epsin γ}{1 + ecos γ}) \frac{ep (\cos γ + e)}{{(1 + ecos γ)}^{2}}}{Ew h^{3} (γ) / 12} d γ]}_{- \frac{π}{2}}^{0} \\ - \int_{- \frac{π}{2}}^{0} \frac{l e^{2} p^{2} \sin ϑ (\cos ϑ + e)}{Ew h^{3} (ϑ) {(1 + ecos ϑ)}^{3} / 6} d ϑ \\ + \int_{- \frac{π}{2}}^{0} \frac{e^{3} p^{3} si n^{2} ϑ (\cos ϑ + e)}{Ew h^{3} (ϑ) {(1 + ecos ϑ)}^{4} / 12} d ϑ + \frac{kep}{Gw} \int_{- \frac{π}{2}}^{0} \frac{\cos ϑ + e}{h (ϑ) {(1 + ecos ϑ)}^{2}} d ϑ \\ = \frac{12 e^{3}}{Ew} I_{5} + \frac{12 e^{3}}{Ew} I_{4} + \frac{ke}{Gw} I_{2} \end{matrix}

(35)

The precision factor corresponding to M_z along the z axis can be calculated as

\begin{matrix} \frac{ε_{z}}{M_{y}} = - {[\frac{epsin ϑ}{1 + ecos ϑ} \int_{- \frac{π}{2}}^{ϑ} \frac{ep (\cos γ + e)}{E w^{3} h (γ) {(1 + ecos γ)}^{2} / 12} d γ]}_{- \frac{π}{2}}^{0} \\ + \int_{- \frac{π}{2}}^{0} \frac{12 e^{2} p^{2} \sin ϑ (\cos ϑ + e)}{E w^{3} h (ϑ) {(1 + ecos ϑ)}^{3}} d ϑ \\ = - \frac{12 e^{2} p}{E w^{3}} I_{6} \end{matrix}

(36)

the precision factor corresponding to force F_z is

\begin{matrix} \frac{ε_{z}}{F_{z}} = {[\frac{epsin ϑ}{1 + ecos ϑ} \int_{- \frac{π}{2}}^{ϑ} \frac{(\frac{l}{2} - \frac{epsin ϑ}{1 + ecos ϑ}) \frac{ep (\cos γ + e)}{{(1 + ecos γ)}^{2}}}{E w^{3} h (γ) / 12} d γ]}_{- \frac{π}{2}}^{0} \\ - \int_{- \frac{π}{2}}^{0} \frac{l e^{2} p^{2} \sin ϑ (\cos ϑ + e)}{E w^{3} h (ϑ) {(1 + ecos ϑ)}^{3} / 6} d ϑ \\ + \int_{- \frac{π}{2}}^{0} \frac{e^{3} p^{3} si n^{2} ϑ (\cos ϑ + e)}{E w^{3} h (ϑ) {(1 + ecos ϑ)}^{4} / 12} d ϑ + \frac{kep}{Gw} \int_{- \frac{π}{2}}^{0} \frac{\cos ϑ + e}{h (ϑ) {(1 + ecos ϑ)}^{2}} d ϑ \\ = \frac{12 e^{3} p^{2}}{E w^{3}} I_{6} + \frac{12 e^{3} p^{2}}{E w^{3}} I_{3} + \frac{ke}{Gw} I_{2} \end{matrix}

(37)

The calculations of I₁–I₆ are shown in Appendix 2.

Verification of closed-form equations

Finite element verification

The compliance equations were verified using the finite element method. ANSYS was employed to calculate the compliance factors of several flexure hinge designs under static loading. We first designed a 3D finite element model in SolidWorks and generated hexahedral grids using the HyperMesh software meshing model. We then established quality point mass, created MPC units, set up the related attributes, and finally performed a grid file import analysis in ANSYS, where each hinge model was fixed at one end and loaded at the other end. Figure 4 shows one of the CFH models. Considering the accuracy and efficiency of the calculation, we set the number of grids to no less than 0.4 million.

Figure 4.

Finite element model of a conic-V-shaped flexure hinge.

Table 1 lists the physical and geometrical parameters of the CFHs. The width (w) of the flexure hinges was 10 mm; the material was spring steel; the Young modulus was 2.07 × 10¹¹ N/m²; the shear modulus was 8.1 × 10¹⁰ N/m². The calculation results, listed in Table 2, were obtained by programming the above closed-form equation.

Table 1.

Design examples of conic-V-shaped flexure hinges.

Example	Type	E	p (mm)	t (mm)	φ (deg)
1	Elliptical-V-shaped	0.4	12.1	0.5	30
2	Elliptical-V-shaped	0.7	6.44	0.5	30
3	Elliptical-V-shaped	0.4	12.1	1	30
4	Elliptical-V-shaped	0.7	6.44	1	30
5	Parabolic-V-shaped	1	6	0.5	30
6	Parabolic-V-shaped	1	8	0.5	30
7	Parabolic-V-shaped	1	6	1	30
8	Parabolic-V-shaped	1	8	1	30
9	Hyperbolic-V-shaped	2	3	0.5	20
10	Hyperbolic-V-shaped	2.5	3	0.5	20
11	Hyperbolic-V-shaped	2	3	1	20
12	Hyperbolic-V-shaped	2.5	3	1	20

Table 2.

Comparison of compliance factors between theoretical results (denoted by C) and finite element results (denoted by F).

Example	α_z/M_z (rad/N·m)	α_y/M_y (rad/N·m)	Δ_z/M_y(×10⁻⁶ 1/N)	Δ_z/F_z (m/N)	Δ_y/M_z(×10⁻⁴ 1/N)	Δ_y/F_y (m/N)	Δ_x/F_x (m/N)	α_x/M_x(rad/N·m)
1(C)	8.427 × 10⁻²	4.530 × 10⁻⁴	–2.192	2.362 × 10⁻⁸	4.079	2.049 × 10⁻⁶	3.775 × 10⁻⁹	6.315 × 10⁻²
1(F)	8.771 × 10⁻²	4.774 × 10⁻⁴	–2.353	2.432 × 10⁻⁸	3.881	2.233 × 10⁻⁶	3.904 × 10⁻⁹	6.700 × 10⁻²
Error (%)	–4.08	–5.39	–7.31	–2.98	4.85	–8.96	8.61	–6.09
2(C)	8.158 × 10⁻²	4.357 × 10⁻⁴	–1.964	2.120 × 10⁻⁸	3.678	1.727 × 10⁻⁶	3.631 × 10⁻⁹	6.113 × 10⁻²
2(F)	7.350 × 10⁻²	4.601 × 10⁻⁴	–2.073	1.980 × 10⁻⁸	3.804	1.780 × 10⁻⁶	3.526 × 10⁻⁹	5.662 × 10⁻²
Error (%)	9.91	–5.59	–5.53	6.61	–3.44	–3.07	2.89	7.38
3(C)	1.482 × 10⁻²	2.916 × 10⁻⁴	–1.411	1.548 × 10⁻⁸	0.717	3.757 × 10⁻⁷	2.430 × 10⁻⁹	1.126 × 10⁻²
3(F)	1.460 × 10⁻²	3.119 × 10⁻⁴	–1.432	1.460 × 10⁻⁸	0.706	3.676 × 10⁻⁷	2.461 × 10⁻⁹	1.174 × 10⁻²
Error (%)	–0.74	–6.92	–2.94	–6.24	–5.89	4.48	–7.36	–6.15
4(C)	1.436 × 10⁻²	2.795 × 10⁻⁴	–1.260	1.383 × 10⁻⁸	0.065	3.182 × 10⁻⁷	2.329 × 10⁻⁹	1.091 × 10⁻²
4(F)	1.420 × 10⁻²	2.988 × 10⁻⁴	–1.273	1.268 × 10⁻⁸	0.064	3.062 × 10⁻⁷	2.461 × 10⁻⁹	1.107 × 10⁻²
Error (%)	1.14	–6.92	–1.05	7.96	1.42	3.78	–5.69	–1.45
5(C)	9.464 × 10⁻²	5.225 × 10⁻⁴	–3.135	3.467 × 10⁻⁸	5.679	3.513 × 10⁻⁶	4.354 × 10⁻⁹	7.093 × 10⁻²
5(F)	9.254 × 10⁻²	5.743 × 10⁻⁴	–3.396	3.427 × 10⁻⁸	5.913	3.667 × 10⁻⁶	4.781 × 10⁻⁹	7.110 × 10⁻²
Error (%)	2.22	–9.91	–8.33	1.13	–4.13	–4.38	–9.79	–0.23
6(C)	1.093 × 10⁻¹	6.196 × 10⁻⁴	–4.957	6.040 × 10⁻⁸	8.743	7.154 × 10⁻⁶	5.164 × 10⁻⁹	8.192 × 10⁻²
6(F)	1.074 × 10⁻¹	6.000 × 10⁻⁴	–4.899	5.493 × 10⁻⁸	8.705	7.299 × 10⁻⁶	5.099 × 10⁻⁹	7.465 × 10⁻²
Error (%)	1.77	3.17	1.17	9.06	0.43	–2.03	1.25	8.88
7(C)	1.673 × 10⁻²	3.397 × 10⁻⁴	–2.038	2.301 × 10⁻⁸	1.004	6.422 × 10⁻⁷	2.831 × 10⁻⁹	1.272 × 10⁻²
7(F)	1.710 × 10⁻²	3.585 × 10⁻⁴	–2.199	2.241 × 10⁻⁸	1.023	6.780 × 10⁻⁷	2.866 × 10⁻⁹	1.235 × 10⁻²
Error (%)	–2.24	–5.53	–7.89	2.60	–1.91	–5.56	–1.24	2.95
8(C)	1.932 × 10⁻²	4.077 × 10⁻⁴	–3.262	4.071 × 10⁻⁸	1.545	1.296 × 10⁻⁶	3.398 × 10⁻⁹	1.470 × 10⁻²
8(F)	1.940 × 10⁻²	4.404 × 10⁻⁴	–3.451	3.701 × 10⁻⁸	1.592	1.305 × 10⁻⁶	3.674 × 10⁻⁹	1.451 × 10⁻²
Error (%)	–0.42	–8.01	–5.80	9.11	–3.00	–0.69	–8.14	1.31
9(C)	9.737 × 10⁻²	5.560 × 10⁻⁴	–3.336	3.729 × 10⁻⁸	5.842	3.634 × 10⁻⁶	4.634 × 10⁻⁹	7.299 × 10⁻²
9(F)	9.254 × 10⁻²	5.492 × 10⁻⁴	–3.396	3.796 × 10⁻⁸	5.850	3.621 × 10⁻⁶	4.543 × 10⁻⁹	7.110 × 10⁻²
Error (%)	4.96	1.22	–1.80	–1.82	–0.13	0.38	1.96	2.59
10(C)	1.098 × 10⁻¹	6.466 × 10⁻⁴	–4.849	5.830 × 10⁻⁸	8.239	6.366 × 10⁻⁶	5.388 × 10⁻⁹	8.236 × 10⁻²
10(F)	9.900 × 10⁻¹	6.491 × 10⁻⁴	–4.575	6.045 × 10⁻⁸	8.311	6.580 × 10⁻⁶	5.396 × 10⁻⁹	8.084 × 10⁻²
Error (%)	9.88	–0.39	5.66	–3.69	–0.88	–3.35	–0.14	1.85
11(C)	1.758 × 10⁻²	3.626 × 10⁻⁴	–2.176	2.484 × 10⁻⁸	1.055	6.835 × 10⁻⁷	3.022 × 10⁻⁹	1.337 × 10⁻²
11(F)	1.709 × 10⁻²	3.585 × 10⁻⁴	–2.205	2.294 × 10⁻⁸	1.029	6.813 × 10⁻⁷	2.866 × 10⁻⁹	1.323 × 10⁻²
Error (%)	2.76	1.15	–1.33	7.64	2.47	0.32	5.17	1.03
12(C)	1.993 × 10⁻²	4.263 × 10⁻⁴	–3.197	3.937 × 10⁻⁸	1.495	1.195 × 10⁻⁶	3.552 × 10⁻⁹	1.517 × 10⁻²
12(F)	1.820 × 10⁻²	4.325 × 10⁻⁴	–3.451	3.974 × 10⁻⁸	1.452	1.274 × 10⁻⁶	3.447 × 10⁻⁹	1.495 × 10⁻²
Error (%)	8.68	–1.47	–7.95	–0.92	2.86	–6.56	2.97	1.50

Owing to the symmetry in the compliance matrix, only Δ_x/F_x, Δ_y/F_y, Δ_z/F_z, α_y/M_y, α_z/M_z, α_x/M_x, Δ_y/M_z, and Δ_z/M_y were verified. The data in Table 2 show that the error in the theoretical calculation and finite element analysis is within 10%, thus confirming the correctness and rationality of the theoretical model.

Experimental verification

To investigate the performance of the proposed generalized model, experimental tests were carried out on four CFH samples. Figure 5 shows the diagram of a sample. In Figure 6, one of the four samples, made of spring steel with an elastic modulus E = 2.07 × 10¹¹ Pa and a shear modulus G = 8.1 × 10¹⁰ N/m², was machined using wire electrical discharge machining with a low-traveling speed technology.

Figure 5.

Diagram of a flexure hinge sample.

Figure 6.

Flexure hinge sample (e = 0.4, p = 12.1 mm, w = 10 mm, t = 1 mm, φ = π/6).

Figure 7 shows the experimental setup, including a laser displacement sensor (KEYENCE^® LK-G155H with a 0.5 μm resolution), a power transformer, an optical platform, two support rods, four flexure hinge samples, and some standard weights (ranging from 0.01 to 0.5 kg). The support rods were mounted on the optical platform. The sensor head was fixed on one support rod using a holder. One end of the hinge sample was fixed on the other support rod using a holder, and the other end of the hinge sample was loaded with weights. As shown in Figure 7(a), a force F_x is applied, including a moment M_z = F_xL. When the displacement of the right bar is measured, the compliance factor Δ_y/M_z can be calculated. Because the deformation is small enough, the angle displacement can be considered α_z ≈ tanα_z and the compliance factor α_z/M_z also can be calculated. In Figure 7(b), the displacement of the right bar can be measured under the force F_y. Thus, the compliance factor Δ_y/F_y can be obtained. The flexure hinges are loaded with ten different loads. The force is evaluated by placing masses ranging from 50 to 500 g. The mass difference between adjacent weights is 50 grams. The compliance factors are calculated as average values. The experimental results are shown in the following table. Table 3 lists the experimental results. The results show that the experimental measurements are within an error of 8% compared to the analytical results.

Figure 7.

Experimental setups for the evaluations of: (a) α_z/M_z and Δ_y/M_z and (b) Δ_y/F_y.

Table 3.

Comparison of compliance factors between experimental results (denoted by X) and theoretical results (denoted by C).

Example	α_z/M_z (×10⁻² rad/Nm)			Δ_y/M_z (×10⁻⁴ 1/N)			Δ_y/F_y (×10⁻⁶ m/N)
Example	X	C	Error (%)	X	C	Error (%)	X	C	Error (%)
2	8.553	8.158	4.62	3.856	3.678	4.62	1.780	1.727	2.96
3	1.599	1.482	7.33	0.774	0.717	7.33	0.394	0.376	4.68
7	1.744	1.673	4.12	1.047	1.004	4.12	0.674	0.642	4.74
11	1.888	1.758	6.90	1.133	1.055	6.90	0.728	0.684	6.11

The following are some of the potential factors that may have caused an uncertainty in the experiments: (1) The difference between the designed samples and processed samples, (2) The accuracy of the measurement system, (3) The difference between the marked value and the actual value of weights.

Numerical simulations and comparisons

To investigate the effects of dimensional parameters on the analytical model, numerical simulations and comparisons are essential. In the generalized model for CFHs, the design parameters are e, p, w, t, and φ. The cutout profile of a CFH is determined by e, p, and φ, whereas the structural characteristics of a flexure hinge are determined by w and t. To investigate the influences of the profile, w, t, and φ were set to 10 mm, 0.5 mm, and π/6, respectively, while both e and p were varied from 0.1 to 1.9. Figure 8 shows the variations in the compliance factors with respect to e and p. Next, we set e, p, w, and t as 0.6, 5.33 mm, 10 mm, and 0.5 mm, respectively, while φ was varied from 0 to 50 deg. Figure 9 shows the variations in the compliance factors with respect to φ. To investigate the effects of the structural characteristics, e, p, and φ were set to 0.6, 5.33 mm, and π/6, respectively, while w and t were varied from 6 to 15 mm and 0.1 to 1 mm, respectively. Figure 10 shows the variations in the compliance factors with respect to w and t.

Figure 8.

Variations in the compliance factors of CFHs with respect to e and p: (a) Δ_x/F_x, (b) Δ_y/F_y, (c) Δ_y/M_z, (d) α_z/M_z, (e) α_y/M_y, (f) Δ_z/M_y, (g) Δ_z/F_z, and (h) α_x/M_x.

Figure 9.

Variations in the compliance factors of CFHs with respect to φ: (a) Δ_x/F_x, (b) Δ_y/F_y, (c) Δ_y/M_z, (d) α_z/M_z, (e) α_y/M_y, (f) Δ_z/M_y, (g) Δ_z/F_z, and (h) α_x/M_x.

Figure 10.

Variations in the compliance factors of CFHs with respect to w and t: (a) Δ_x/F_x, (b) Δ_y/F_y, (c) Δ_y/M_z, (d) α_z/M_z, (e) α_y/M_y, (f) Δ_z/M_y, (g) Δ_z/F_z, and (h) α_x/M_x.

The compliance factors nonlinearly increase with increasing eccentricity e and increasing distance p between the corresponding directrix and the focus, as shown in Figure 8. Moreover, the variation sensitivities of Δ_x/F_x, α_z/M_z, α_y/M_y, and α_x/M_x to e and p decrease, while the variation sensitivities of Δ_y/F_y, Δ_y/M_z, Δ_z/M_y, and Δ_z/F_z to e and p increase. The eccentricity e and the distance p between the corresponding directrix and the focus have largely the same influence on the performances of the conic-V-shaped flexure hinges.

Figure 9 shows that the compliance factors Δ_x/F_x, Δ_y/F_y, Δ_y/M_z, α_z/M_z, α_y/M_y, Δ_z/M_y, Δ_z/F_z, and α_x/M_x nonlinearly decrease with increasing angle φ, while their variation sensitivities to the angle φ decrease.

As shown in Figure 10, the compliance factors nonlinearly decrease with increasing cross-sectional width w and increasing minimum thickness t, while their variation sensitivities to cross-sectional width w and minimum thickness t decrease. The parameter t has a more significant influence on the compliance factors Δ_y/F_y, Δ_y/M_z, α_z/M_z, and α_x/M_x than the cross-sectional width w.

Conclusion

This paper presents a new class of flexure hinges, namely conic-V-shaped flexure hinges, which can be used as a generalized flexure hinge model for elliptical-V-shaped, parabolic-V-shaped, and hyperbolic-V-shaped. The following conclusions can be summarized:

The compliance and precision equations for the CFHs were derived using Castigliano’s second theorem. The factors in the compliance equations were verified by finite element analysis and experimental measurements. The errors in the finite element analysis and experimental results were within 10% and 8% compared to the analytical results, respectively. This demonstrates the practicality of the compliance and precision equations derived for CFHs.

Based on the analytical results, numerical simulations and comparisons are made to investigate the effects of dimensional parameters on the analytical model. The results showed that the compliance factors increased with the increase of the eccentricity e and the distance between the corresponding directrix and the focus p, but decreased with the increase of the angle φ, the minimum thickness t, and the cross-sectional width w.

Footnotes

Appendix 1

Appendix 2

I₁–I₆ are used to calculated factors in the compliance and precision matrixes as shown in sections 3 and 4. And I₁-I₆ can be expressed as

(2.1)

\begin{matrix} I_{1} = \int_{0}^{ϑ_{m}} \frac{e + \cos ϑ}{{g_{1}}^{3} (ϑ) {(1 + ecos ϑ)}^{2}} d ϑ + \int_{ϑ_{m}}^{\frac{π}{2}} \frac{e + \cos ϑ}{{g_{2}}^{3} (ϑ) {(1 + ecos ϑ)}^{2}} d ϑ \\ = - \frac{6 e (e^{2} (- 2 + m) - m) arctanh (\frac{(e (- 2 + m) - m) \tan (ϑ_{m} / 2)}{\sqrt{e^{2} {(- 2 + m)}^{2} - m^{2}}})}{{(e^{2} {(- 2 + m)}^{2} - m^{2})}^{5 / 2}} + \frac{(e^{4} {(- 2 + m)}^{2} (- 1 + m) + m^{3} + e^{2} m (2 + 5 m - 2 m^{2})) \sin (ϑ_{m})}{{(- e^{2} {(- 2 + m)}^{2} + m^{2})}^{2} {(m + e (- 2 + m) \cos (ϑ_{m}))}^{2}} \\ + \frac{e (e^{4} {(- 2 + m)}^{3} + (- 1 + m) m^{2} - e^{2} (8 + 2 m - 7 m^{2} + 2 m^{3})) \cos (ϑ_{m}) \sin (ϑ_{m})}{{(- e^{2} {(- 2 + m)}^{2} + m^{2})}^{2} {(m + e (- 2 + m) \cos (ϑ_{m}))}^{2}} \\ + \frac{2 co s^{3} φ (encos φ - (1 - e^{2}) \sin φ)}{e (4 + n^{2} + (- 4 + n^{2}) \cos 2 φ) {(ncos φ + 2 esin φ)}^{2}} \\ - \frac{2 co s^{3} φ (en (1 + ecos ϑ_{m}) \cos φ \sin ϑ_{m} - (1 + 2 ecos ϑ_{m} + e^{2} \cos 2 ϑ_{m}) \sin φ)}{e (4 + n^{2} + (- 4 + n^{2}) \cos 2 φ) {(n (1 + ecos ϑ_{m}) \cos φ + 2 esin ϑ_{m} \sin φ)}^{2}} \end{matrix}

(2.2)

\begin{matrix} I_{2} = \int_{0}^{ϑ_{m}} \frac{e + \cos ϑ}{g_{1} (ϑ) {(1 + ecos ϑ)}^{2}} d ϑ + \int_{ϑ_{m}}^{\frac{π}{2}} \frac{e + \cos ϑ}{g_{2} (ϑ) {(1 + ecos ϑ)}^{2}} d ϑ \\ = \frac{\sqrt{e^{2} - 1} \sqrt{e^{2} {(m - 2)}^{2} - m^{2}} arctanh (\frac{(e - 1) \tan (ϑ_{m} / 2)}{\sqrt{e^{2} - 1}})}{e \sqrt{e^{2} {(m - 2)}^{2} - m^{2}}} \\ + \frac{(- e^{2} (m - 2) + m) arctanh (\frac{(e (m - 2) - m) \tan (ϑ_{m} / 2)}{\sqrt{e^{2} {(m - 2)}^{2} - m^{2}}})}{e \sqrt{e^{2} {(m - 2)}^{2} - m^{2}}} \\ + \frac{\cot φ \log (((en + 2 e) \cos (π / 2 - φ) + 2 ncos φ) (1 + ecos ϑ_{m}))}{2 e} \\ + \frac{\cot φ \log (en - 2 e) \cos (π / 2 + φ) (1 + ecos ϑ_{m})}{2 e} \\ - \frac{\cot φ \log ((en + 2 e) \cos (ϑ_{m} - φ) + 2 ncos φ + (en - 2 e) \cos (ϑ_{m} + φ))}{2 e} \end{matrix}

(2.3)

\begin{matrix} I_{3} = \int_{0}^{ϑ_{m}} \frac{si n^{2} ϑ (\cos ϑ + e)}{g_{1} (ϑ) {(1 + ecos ϑ)}^{4}} d ϑ + \int_{ϑ_{m}}^{\frac{π}{2}} \frac{si n^{2} ϑ (\cos ϑ + e)}{g_{2} (ϑ) {(1 + ecos ϑ)}^{4}} d ϑ \\ = \frac{(e^{4} {(m - 2)}^{2} - 2 e^{2} {(m - 1)}^{2} + m^{2}) arctanh (\frac{(e - 1) \tan (ϑ_{m} / 2)}{\sqrt{e^{2} - 1}})}{4 e^{3} \sqrt{e^{2} - 1}} \\ - \frac{(e^{4} {(m - 2)}^{3} + m^{3} - 2 e^{2} m (2 - 3 m + m^{2})) arctanh (\frac{(e (m - 2) - m) \tan (ϑ_{m} / 2)}{\sqrt{e^{2} {(m - 2)}^{2} - m^{2}}})}{4 e^{3} \sqrt{e^{2} {(m - 2)}^{2} - m^{2}}} \\ + \frac{e (e^{2} - 1) \sin ϑ_{m}}{4 e^{3} {(1 + ecos ϑ_{m})}^{2}} + \frac{e (1 - e^{2} (m - 2) + m) \sin ϑ_{m}}{4 e^{3} (1 + ecos ϑ_{m})} \\ + \frac{\cot φ (2 (e^{2} + 1) + n^{2} co t^{2} φ \log (2 (ncos φ + 2 esin φ)) - 2 encot φ)}{8 e^{3}} \\ - \frac{\frac{2 \cot φ (e^{2} + 1 + 2 ecos ϑ_{m})}{{(1 + ecos ϑ_{m})}^{2}} - n^{2} co t^{3} φ \log (1 + ecos ϑ_{m})}{8 e^{3}} \\ + \frac{n^{2} co t^{3} φ \log (2 (n (1 + ecos ϑ_{m}) \cos φ + 2 esin φ \sin ϑ_{m}))}{8 e^{3}} + \frac{enco t^{2} φ \sin ϑ_{m}}{4 e^{3} (1 + ecos ϑ_{m})} \end{matrix}

(2.4)

\begin{matrix} I_{4} = \int_{0}^{ϑ_{m}} \frac{si n^{2} ϑ (\cos ϑ + e)}{{g_{1}}^{3} (ϑ) {(1 + ecos ϑ)}^{4}} d ϑ + \int_{ϑ_{m}}^{\frac{π}{2}} \frac{si n^{2} ϑ (\cos ϑ + e)}{{g_{2}}^{3} (ϑ) {(1 + ecos ϑ)}^{4}} d ϑ \\ = \frac{{(e^{2} - 1)}^{3 / 2} arctanh (\frac{(e - 1) \tan (ϑ_{m} / 2)}{\sqrt{e^{2} - 1}})}{4 e^{3}} - \frac{(e^{6} {(m - 2)}^{3} - m^{3}) arctanh (\frac{(e (m - 2) - m) \tan (ϑ_{m} / 2)}{\sqrt{e^{2} {(m - 2)}^{2} - m^{2}}})}{4 e^{3} {(e^{2} {(m - 2)}^{2} - m^{2})}^{3 / 2}} \\ - \frac{(3 e^{2} m (2 - 2 m + m^{2}) - 3 e^{4} (6 m - 4 - 4 m^{2} + m^{3})) arctanh (\frac{(e (m - 2) - m) \tan (ϑ_{m} / 2)}{\sqrt{e^{2} {(m - 2)}^{2} - m^{2}}})}{4 e^{3} {(e^{2} {(m - 2)}^{2} - m^{2})}^{3 / 2}} \\ + \frac{co t^{3} φ \log (2 (ncos φ + 2 esin φ))}{8 e^{3}} - \frac{co t^{2} φ \cos φ \cos (3 φ) (e n^{3} - 8 en)}{8 e^{3} (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ + 2 esin φ)}^{2}} \\ - \frac{co t^{2} φ \cos φ (\sin (3 φ) (3 e^{2} n^{2} - 16 e^{2} + n^{2}) + \cos φ (8 en + 3 e n^{3}) + \sin φ (48 e^{2} + 3 e^{2} n^{2} + n^{2}))}{8 e^{3} (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ + 2 esin φ)}^{2}} \\ - \frac{co t^{2} φ (- \cot φ \log (1 + ecos ϑ_{m}) + \cot φ \log (2 (n (1 + ecos ϑ_{m}) \cos φ + 2 esin ϑ_{m} \sin φ)))}{8 e^{3}} \\ + \frac{co t^{2} φ \cos φ (2 en (n^{2} - 2 n - 8) \sin (ϑ_{m} - 3 φ) + e^{2} (n^{3} + 2 n^{2} - 8 n - 16) \sin (2 ϑ_{m} - 3 φ))}{32 e^{3} (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ (1 + ecos ϑ_{m}) + 2 esin φ \sin ϑ_{m})}^{2}} \\ + \frac{co t^{2} φ \cos φ (16 en - 4 e n^{2} + 6 e n^{3}) \sin (ϑ_{m} - φ)}{32 e^{3} (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ (1 + ecos ϑ_{m}) + 2 esin φ \sin ϑ_{m})}^{2}} \\ + \frac{co t^{2} φ \cos φ ((96 e^{2} + 4 n^{2} + 8 e^{2} n^{2}) \sin φ + (4 n^{2} - 32 e^{2} + 8 e^{2} n^{2}) \sin (3 φ))}{32 e^{3} (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ (1 + ecos ϑ_{m}) + 2 esin φ \sin ϑ_{m})}^{2}} \\ + \frac{co t^{2} φ \cos φ ((16 en + 4 e n^{2} + 6 e n^{3}) \sin (ϑ_{m} + φ) + (48 e^{2} + 8 e^{2} n + 2 e^{2} n^{2} + 3 e^{2} n^{3}) \sin (2 ϑ_{m} - φ))}{32 e^{3} (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ (1 + ecos ϑ_{m}) + 2 esin φ \sin ϑ_{m})}^{2}} \\ + \frac{co t^{2} φ \cos φ ((8 e^{2} n - 48 e^{2} - 2 e^{2} n^{2} + 3 e^{2} n^{3}) \sin (2 ϑ_{m} + φ))}{32 e^{3} (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ (1 + ecos ϑ_{m}) + 2 esin φ \sin ϑ_{m})}^{2}} \\ + \frac{co t^{2} φ \cos φ ((4 e n^{2} - 16 en + 2 e n^{3}) \sin (ϑ_{m} + 3 φ) + (16 e^{2} - 8 e^{2} n - 2 e^{2} n^{2} + e^{2} n^{3}) \sin (2 ϑ_{m} + 3 φ))}{32 e^{3} (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ (1 + ecos ϑ_{m}) + 2 esin φ \sin ϑ_{m})}^{2}} \end{matrix}

(2.5)

\begin{matrix} I_{5} = \int_{0}^{ϑ_{m}} \frac{\sin ϑ (\cos ϑ + e)}{{g_{1}}^{3} (ϑ) {(1 + ecos ϑ)}^{3}} d ϑ + \int_{ϑ_{m}}^{\frac{π}{2}} \frac{\sin ϑ (\cos ϑ + e)}{{g_{2}}^{3} (ϑ) {(1 + ecos ϑ)}^{3}} d ϑ \\ = \frac{(e^{2} + 2 ecos ϑ_{m}) (m - 2) + m}{2 e^{2} {(m - 2)}^{2} {(m + e (m - 2) \cos ϑ_{m})}^{2}} - \frac{(e^{2} + 2 e) (m - 2) + m}{2 e^{2} {(m - 2)}^{2} {(m + e (m - 2))}^{2}} \\ + \frac{co s^{2} φ (8 e^{2} + n^{2} + e^{2} n^{2} + 4 ensin (2 φ) + (n^{2} - 8 e^{2} + e^{2} n) \cos (2 φ))}{2 e^{2} n (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ + 2 esin φ)}^{2}} \\ - \frac{co s^{2} φ (4 e^{2} + n^{2} + e^{2} n^{2} + 2 e n^{2} \cos ϑ_{m} - 4 e^{2} \cos (2 ϑ_{m}) + (2 en + e n^{2}) \cos (ϑ_{m} - 2 φ))}{2 e^{2} n (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ (1 + ecos ϑ_{m}) + 2 esin φ \sin ϑ_{m})}^{2}} \\ - \frac{co s^{2} φ ((2 e^{2} + e^{2} n) \cos (2 ϑ_{m} - 2 φ) + (n^{2} - 4 e^{2} + e^{2} n^{2}) \cos (2 φ))}{2 e^{2} n (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ (1 + ecos ϑ_{m}) + 2 esin φ \sin ϑ_{m})}^{2}} \\ - \frac{co s^{2} φ ((2 e^{2} - e^{2} n) \cos (2 ϑ_{m} + 2 φ) + (e n^{2} - 2 en) \cos (ϑ_{m} + 2 φ))}{2 e^{2} n (4 + n^{2} + (n^{2} - 4) \cos (2 φ)) {(ncos φ (1 + ecos ϑ_{m}) + 2 esin φ \sin ϑ_{m})}^{2}} \end{matrix}

(2.6)

\begin{matrix} I_{6} = \int_{0}^{ϑ_{m}} \frac{\sin ϑ (\cos ϑ + e)}{g_{1} (ϑ) {(1 + ecos ϑ)}^{3}} d ϑ + \int_{ϑ_{m}}^{\frac{π}{2}} \frac{\sin ϑ (\cos ϑ + e)}{g_{2} (ϑ) {(1 + ecos ϑ)}^{3}} d ϑ \\ = \frac{\frac{2 (e^{2} - 1)}{1 + ecos ϑ_{m}} + (e^{2} (m - 2) - m) \log (1 + ecos ϑ_{m})}{4 e^{2}} \\ - \frac{\frac{2 (e^{2} - 1)}{1 + e} + (e^{2} (m - 2) - m) \log (1 + e) + (m - e^{2} (m - 2)) \log (m + e (m - 2))}{4 e^{2}} \\ + \frac{\cot φ (2 e - ncos φ \log (2 (ncos φ + 2 esin φ)))}{4 e^{2}} \\ + \frac{(m - e^{2} (m - 2)) \log (m + e (m - 2) \cos ϑ_{m}) - \cot φ (ncot φ \log (1 + ecos ϑ_{m}))}{4 e^{2}} \\ + \frac{\cot φ (ncot φ \log (2 (n (1 + ecos ϑ_{m}) \cos φ + 2 esin ϑ_{m} \sin φ)) - \frac{2 esin ϑ_{m}}{1 + ecos ϑ_{m}})}{4 e^{2}} \end{matrix}

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (51705378 and 51875418), Hubei Provincial Natural Science Foundation of China (2018CFB757), Postdoctoral Sustentation Fund of China (2017M622532), Youth Fund Project for the State Key Laboratory of Refractories and Metallurgy (2018QN15).

ORCID iD

Puliang Yu

Author biographies

Jianyi Kong is a Professor in Wuhan University of Science and Technology, researching in mechanical design and optimization, and intelligent design and manufacturing.

Zhao Huang is a Postgraduate Student in Wuhan University of Science and Technology, researching in numerical simulation of flexure hinges.

Xiaodong Xian is a Postgraduate Student in Wuhan University of Science and Technology, researching in numerical simulation and experiment of flexure hinges.

Yingrui Wang is a Postgraduate Student in Wuhan University of Science and Technology, researching in mechanical design of flexure hinges.

Puliang Yu is an Assistant Professor in Wuhan University of Science and Technology. His main research interests in intelligent design and manufacturing, fluid mechanics and flow control.

References

Zhang

Guo

, et al. Design, modeling, and analysis of hybrid flexure hinges. Mech Mach Theory 2019; 131: 300–316.

Zhang

Fang

, et al. Multi-objective optimization of flexure hinge mechanism considering thermal–mechanical coupling deformation and natural frequency. Adv Mech Eng 2017; 137: 012301.

Zhang

Cai

, et al. Modeling and analysis of conical-shaped notch flexure hinges based on NURBS. Mech Mach Theory 2018; 128: 560–568.

Salinic

Nikolic

A new pseudo-rigid-body model approach for modeling the quasi-static response of planar flexure-hinge mechanisms. Mech Mach Theory 2018; 124: 150–161.

Lai

Zhu

ZN.

Design, modeling and testing of a novel flexure-based displacement amplification mechanism. Sens Actuators A Phys 2017; 266: 122–129.

Chen

Liu

Elliptical-arc-fillet flexure hinges: toward a generalized model for commonly used flexure hinges. J Mech Des 2011; 133: 081002.

Design, modeling, and analysis of a novel microgripper based on flexure hinges. Adv Mech Eng 2014; 2014: 1–11.

Jiang

, et al. Design and analysis of a 2-degree-of-freedom flexure-based micro-motion stage. Adv Mech Eng 2016; 8: 1–13.

Yao

Wang

, et al. Analysis of the displacement amplification ratio of bridge-type flexure hinge. Sens Actuators A Phys 2006; 132: 730–736.

10.

Verotti

Crescenzi

Balucani

, et al. MEMS-based conjugate surfaces flexure hinge. J Mech Des 2015; 137(1): 012301.

11.

Liu

Zhang

Fatikow

Design and analysis of a multi-notched flexure hinge for compliant mechanisms. Precis Eng 2017; 48: 292–304.

12.

Jouaneh

Yang

Modeling of flexure-hinge type lever mechanisms. Precis Eng 2003; 27: 407–418.

13.

Paros

Weisbord

How to design flexure hinges. Mach Des 1965; 37: 151–156.

14.

Smith

Badami

Dale

, et al. Elliptical flexure hinges. Rev Sci Instrum 1997; 68: 1474–1483.

15.

Chen

Shao

Huang

A new generalized model for elliptical arc flexure hinges. Rev Sci Instrum 2008; 79: 095103.

16.

Lobontiu

Paine

Garcia

, et al. Design of symmetric conic-section flexure hinges based on closed-form compliance equations. Mech Mach Theory 2002; 37: 477–498.

17.

Wang

Zhou

Zhu

Development of a novel sort of exponent-sine-shaped flexure hinges. Rev Sci Instrum 2013; 84: 095008.

18.

Pan

Closed-form compliance equations for power-function-shaped flexure hinge based on unit-load method. Precis Eng 2013; 37: 135–145.

19.

Lobontiu

Paine

JSN

Garcia

, et al. Corner-filleted flexure hinges. J Mech Des 2001; 123: 346–352.

20.

Yong

Handley

DC.

Review of circular flexure hinge design equations and derivation of empirical formulations. Precis Eng 2008; 32: 63–70.

21.

Hsiao

Lin

TW.

Analysis of a novel flexure hinge with three degrees of freedom. Rev Sci Instrum 2001; 72: 1565–1573.

22.

Liu

Zhang

Fatikow

Design and analysis of a high-accuracy flexure hinge. Rev Sci Instrum 2016; 87: 770–704.

23.

Kim

Kwak

YK.

Development and optimization of 3-D bridge-type hinge mechanisms. Sens Actuators A Phys 2004; 116:530–538.

24.

Ryu

Lee

Gweon

, et al. Inverse kinematic modeling of a coupled flexure hinge mechanism. Mechatronics 1999; 9: 657–674.

25.

Vallance

Haghighian

Marsh

ER.

A unified geometric model for designing elastic pivots. Precis Eng 2008; 32: 278–288.

26.

Chen

Liu

Gao

, et al. A generalized model for conic flexure hinges. Rev Sci Instrum 2009; 80: 055106.

27.

Tian

Shirinzadeh

Zhang

Closed-form compliance equations of filleted V-shaped flexure hinges for compliant mechanism design. Precis Eng 2010; 34: 408–418.

28.

Tian

Shirinzadeh

Zhang

, et al. Three flexure hinges for compliant mechanism designs based on dimensionless graph analysis. Precis Eng 2010; 34: 92–100.

29.

Spearing

SM.

Materials issues in microelectromechanical systems (MEMS). Acta Mater 2000; 48: 179–196.

30.

Hearn

EJ.

Mechanics of materials 2: an introduction to the mechanics of elastic and plastic deformation of solids and structural materials. Boston: Butterworth-Heinemann, 1997, p.143.

Generalized model for conic-V-shaped flexure hinges

Abstract

Keywords

Introduction

The generalized equation in polar coordinates

Compliance equations for conic-V-shaped flexure hinges

Compliance matrix

Compliance equations

Angular compliance about the z axis

Angular compliance about the y axis

Linear compliance along the z axis

Linear compliance along the y axis

Linear compliance along the x axis

Angular compliance about the x axis

Precision equations for conic-V-shaped flexure hinges

Verification of closed-form equations

Finite element verification

Experimental verification

Numerical simulations and comparisons

Conclusion

Footnotes

Appendix 1

Appendix 2

Declaration of conflicting interests

Funding

ORCID iD

Author biographies

References