Next Article in Journal
Design and Fabrication of an Integrated Hollow Concave Cilium MEMS Cardiac Sound Sensor
Next Article in Special Issue
Relationship between the Young’s Moduli of Whole Microcapsules and Their Shell Material Established by Micromanipulation Measurements Based on Diametric Compression between Two Parallel Surfaces and Numerical Modelling
Previous Article in Journal
Miniature Manipulator Design and Cartesian Control for Minimally Invasive Transluminal Endoscopic Surgery
Previous Article in Special Issue
A Substructure Condensed Approach for Kinetostatic Modeling of Compliant Mechanisms with Complex Topology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Direct Kinetostatic Analysis of a Gripper with Curved Flexures

by
Alessandro Cammarata
1,*,
Pietro Davide Maddio
1,
Rosario Sinatra
1 and
Nicola Pio Belfiore
2
1
Department of Civil Engineering and Architecture, University of Catania, Viale Andrea Doria, 6, 95123 Catania, Italy
2
Department of Industrial, Electronic and Mechanical Engineering, University of Roma Tre, Via Vito Volterra, 62, 00146 Rome, Italy
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(12), 2172; https://doi.org/10.3390/mi13122172
Submission received: 3 November 2022 / Revised: 22 November 2022 / Accepted: 3 December 2022 / Published: 8 December 2022
(This article belongs to the Special Issue Flexible Micromanipulators and Micromanipulation)

Abstract

:
Micro-electro-mechanical-systems (MEMS) extensively employed planar mechanisms with elastic curved beams. However, using a curved circular beam as a flexure hinge, in most cases, needs a more sophisticated kinetostatic model than the conventional planar flexures. An elastic curved beam generally allows its outer sections to experience full plane mobility with three degrees of freedom, making complex non-linear models necessary to predict their behavior. This paper describes the direct kinetostatic analysis of a planar gripper with an elastic curved beam is described and then solved by calculating the tangent stiffness matrix in closed form. Two simplified models and different contributions to derive their tangent stiffness matrices are considered. Then, the Newton–Raphson iterative method solves the non-linear direct kinetostatic problem. The technique, which appears particularly useful for real-time applications, is finally applied to a case study consisting of a four-bar linkage gripper with elastic curved beam joints that can be used in real-time grasping operations at the microscale.

1. Introduction

Compliant mechanisms have been adopted in several applications for centuries, because of several well-known advantages, such as the absence of sliding friction, backlash, and significant wear, with the advantage of requiring minimum effort assembly. The first complete and systematic description of their characteristics appeared in the literature about thirty years ago, around 1994, when Midha and Howell introduced a classification [1] that has been used until today. Another significant step forward in developing the compliant mechanisms was taken a couple of years later, when the Pseudo-Rigid-Body equivalent Model (PRBM) was introduced to evaluate the elasticity of compliant mechanisms with significant deflection capabilities [2]. More materials concerning the compliant mechanisms were presented in 2001 [3].
Compliant mechanisms are nowadays used in many fields [4], as, for instance, for the development of aerospace [5] and biomedical [6,7] devices, compliant bistable mechanisms [8,9,10,11], grasping and releasing micro-objects devices [12], precision engineering [13], MEMS [14,15,16,17,18,19], polishing and deburring [20,21,22], lab-on-chip micro and nanosystems [23], and automotive devices [24,25].
The actual configuration of a compliant mechanism depends not only on the applied forces and torques but also on its geometric characteristics, with kinematic and mechanical coupling and non-linearity problems that especially arise in case of large deformations.
A comprehensive survey on many recent techniques for modeling the kinetostatic and dynamic behavior of flexure-based compliant mechanisms has been recently presented [26].
Kinetostatic models of complex plane compliant mechanisms have been developed in both micro and macro scale devices by using a wide variety of different linear and non-linear methods, such as, for only representative example, those based on: Loop-closure equations and the static equilibrium conditions for multi-loops compliant mechanisms [27]; chained beam constraint model and geometric parameter optimization, specially conceived for translational motion [28]; a combination of a beam constraint model, load equilibrium conditions, and geometric compatibility equations, specially conceived for 3-PPR compliant parallel mechanisms [29]; a kinetostatic modeling approach that integrates the screw theory with the energy method, with consequent avoidance of the problem of finding a solution to the equilibrium equations of nodal forces and the possibility of taking into account the parasitic deformations in space [30]; an extension of the chained beam constraint model specially revisited to analyze flexible beams of effective variable length [31]; a mathematical formulation of the compliance matrix method, combined with the inverse kinematic, specially introduced for modeling the flexure-based parallel compliant mechanisms with multiple actuation forces [32]; the adoption of three representations of multiple segments 2D beam models, namely, beam constraint model, linear Euler–Bernoulli beam, and PRBM [33]; the adoption of a new two-colored digraph representation of planar flexure-based compliant mechanisms for the automatic generation of the kinetostatic equations [34]; the introduction of virtual flexure hinges, link-flexure incidence matrices, and path matrices to generate automatically the formulation of the kinetostatic equations [35].
A more specific contribution has been dedicated in 2013 [36] to the solution of the problem of inverse kinetostatic analysis of a compliant four-bar linkage with flexible circular joints and pseudo-rigid bodies. This problem was attacked by extensively applying the theory of curved beams to the flexible parts, which gave rise to the closed-form symbolic expression of the compliance matrix, and by applying the static balance equations to both the elastic and pseudo-rigid parts.
The present investigation explores the opposite problem of direct kinetostatic analysis of a planar gripper with circular flexures. Two possible models based on the static balance of flexures are provided. The first linear model considers the static equilibrium in the undeformed configuration, while the second considers the balance in the deformed configuration. Both models simplify the fully non-linear model by exploiting a constant stiffness matrix of the undeformed curved beam element. These models allow us to find the tangent stiffness matrix in closed form as the sum of different contributions. Furthermore, dividing the tangent stiffness matrix into its contributions allows for evaluating each term’s importance and setting strategies to speed up convergence. Furthermore, through a validation process of the fully non-linear model results performed on a case study, it will be possible to ascertain how the simplifications still provide accurate values in almost the entire mechanism’s range of motion. As known, the tangent stiffness matrix is the heart of an iterative solving method. It is the basis of many implicit integrators widely used for the study of flexible mechanisms, such as the generalized α -method [37] or the HHT-method [38].
The main target of this article is to create two simplified models:
  • Solving the problem of the direct kinetostatic analysis of planar grippers with curved beams;
  • Being reliable in terms of motion accuracy and actuation forces;
  • Being computationally efficient to extend the formulation for real-time applications.
Any Finite Element Analysis (FEA) or Multibody Dynamics Simulation (MBDS) package is very reliable for solving any general problem in kinetostatic analysis numerically. Despite this, the availability of a ready-to-use independent algorithm to solve the direct kinetostatic problem gives rise to the possibility of implementing it in any real-time applications. Nevertheless, the MBDS Adams software has been used herein for validation purposes.
The paper is divided into the following sections. Section 2 gives the fundamentals of the curved beam model. Section 3 outlines the kinetostatic analysis. Two simplified linear and partial non-linear models are developed, and their tangent stiffness matrices are obtained in closed form. Section 4 includes a detailed case study description. Section 5 compares and validates the two models and gives essential insights into convergence and computational burden. Finally, Section 6 gives the concluding remarks.

2. The Adopted Curved Beam Model

Flexures employed in this context are curved beams. It has been demonstrated that curved beams can provide large rotations while maintaining small errors in terms of displacements of its center [39], as it is typical for classic revolute pairs. This feature is important to guarantee finite rotations of the end-effector in monolithic structures such as MEMS-based grippers. Furthermore, a linear model is capable of faithfully reproducing the displacements and in-plane rotation of the curved beam tip up to rotations of approximately ± 20 . This feature has the considerable advantage of using a constant stiffness matrix, as will be recalled below.
In the following, the curved beam compliance matrix, and its inverse stiffness matrix, will be recalled from [36]. Let us consider a curved beam with a circular profile of radius r f and beam characteristic angle θ f , as displayed in Figure 1. First, let us consider the generalized displacement array ψ f = [ ξ ^ f T , ϕ f ] T containing the displacement ξ ^ f and the rotation angle ϕ f of the end section due to the deformation. Then, introducing the generalized wrench array w f = [ F f T , M f ] T containing the force vector and the torque applied to the end section, the compliance matrix C f , derived in [36], follows from
ψ f = C f w f
and depends only on the geometric and structural parameters of the curved beam, i.e.,
C f = 1 E I r f 3 4 ( 6 θ f + s ( 2 θ f ) 8 s ( θ f ) ) r f 3 2 ( c 2 ( θ f ) 2 c ( θ f ) + 1 ) r f 2 ( θ f s ( θ f ) r f 3 4 ( 2 θ f s ( 2 θ f ) ) r f 2 ( 1 c ( θ f ) ) ( sym ) r f θ f
where E is Young’s modulus and I is the area moment of inertia, assumed both constant for the circular profile. In the following sections, the inverse of the compliance matrix C f , i.e., the stiffness matrix K f will be employed to write the kinetostatic equations of planar mechanisms with curved beams. Furthermore, the stiffness matrix will be expressed in its locale frame S ^ f attached to the end-section of the curved beam in the undeformed configuration, as shown in Figure 1.

3. Kinetostatic Analysis

Hereafter, all vectors denoted with the hat will refer to the undeformed configuration, while the same vectors will indicate the deformed configuration without the hat. Position vectors of local frames as well as rotation matrices of frames describing the orientation of bodies in the undeformed configuration are constant.
A curved beam links two components, as it happens for the two bodies displayed in Figure 2. First, consider the undeformed system composed of two rigid bodies, identified by the reference frames S ^ i and S ^ j , and by the flexure f ^ . The vectors r ^ i and r ^ j denote the positions of the body reference frame origins with respect to the fixed reference frame Σ . In contrast, the vectors s ^ i f and s ^ j f , respectively, indicate the distance vectors going from the body-reference frame origins to the attachment points of the curved beam to the bodies.
In the undeformed configuration, the position vector p ^ f going from the attachment point on body i to that on body j is obtained through the following expression
p ^ f = r ^ j + s ^ j f r ^ i s ^ i f
Then, consider a generic configuration in which the bodies undergo finite displacements and rotations, and the flexure is deformed. For the assumption of rigid bodies, it follows that s ^ i f ( S ^ i ) s i f ( S i ) s ¯ i f and s ^ j f ( S ^ j ) s j f ( S j ) s ¯ j f where the superscript denotes the reference frame in which the vector is expressed. In the previous expressions, s ¯ i f and s ¯ j f have been introduced to simplify the notation. Then, the following closure equation stands,
r i + A i s ¯ i f + p f r j A j s ¯ j f = 0
where A i and A j are the rotation matrices mapping S i and S j into Σ and p f is the distance vector between the two flexure extremities in the deformed configuration.
Let us introduce the deformation vector x f containing the deformations of the flexure due to the displaced configuration described in Section 2. As known from the continuum mechanics, this vector can be represented using either the material or the spatial description of motion. In the following, only the material description is implemented. Therefore, expressing x f in the frame S ^ i f of Figure 2, it follows
x f = A ^ i A i T p f p ^ f
where p f has been pulled back to the undeformed configuration as required in the material description of motion.
Then, as recalled in Section 2, the circular flexure model requires x f to be expressed in the frame S ^ f j instead of S ^ i f , therefore
x f ( S ^ f j ) = A ^ f j T A ^ j T x f
where A ^ f j is the constant rotation matrix mapping S ^ f j into S ^ j . The rotation angle ϕ f due to the flexure deformation reads
ϕ f = θ j θ i θ ^ j + θ ^ i Δ θ i j Δ θ ^ i j
where θ and θ ^ , respectively, are the rotation angles of the bodies in the spatial and material configurations and Δ θ , Δ θ ^ denote the corresponding relative rotation angles. The generalized displacement array of the curved beam f, already introduced in Section 2, becomes
ψ f ( S ^ f j ) = x f ( S ^ f j ) ϕ f

3.1. Jacobian of the Deformation Vector

Since the direct kinetostatic analysis will be solved using an iterative procedure, the variation of the generalized displacement array must be calculated. Using Equation (5), the variation δ x f is
δ x f = A ^ i A ¯ i T p f δ θ i + A ^ i A i T ( δ r j δ r i + A ¯ j s ¯ j f δ θ j A ¯ i s ¯ i f δ θ i )
where A ¯ = A / θ while δ r , δ θ are the variations of the body coordinates in the deformed configuration. The variation δ x f is evaluated in the reference frame Σ but can be easily expressed in S ^ f j remembering that A ^ f j and A ^ j in Equation (6) are constant, i.e.,
δ x f ( S ^ f j ) = A ^ f j T A ^ j T δ x f
Considering the angles, the variation of Equation (7), leads to
δ ϕ f = δ θ j δ θ i δ Δ θ i j
Finally, the variations can be combined to form the variation of the generalized deformation vector δ ψ f . The latter satisfies the following expression
δ ψ f = J ˘ f δ q i j
where δ q i j = [ δ q i T δ q j T ] T is a 6-dimensional vector containing the variations of the body coordinates of bodies i and j and J ˘ f is the ( 3 × 6 ) Jacobian matrix, defined as
J ˘ f = J ˘ f i | J ˘ f j = A ^ i A i T A ^ i ( A ¯ i T p f 1 ˜ s ¯ i f ) 0 T 1 | A ^ i A i T A ^ i A i T A ¯ j s ¯ j f 0 T 1
In deriving J ˘ f , the property A i T A ¯ i = 1 ˜ has been employed, being
1 ˜ = 0 1 1 0
a particular skew-symmetric matrix used to define the cross-product in the planar case.

3.2. Kinetostatic Equations

The kinetostatic equations of the system require the static equilibrium of a curved beam. Consider the layout of Figure 3 showing the static balance of a curved beam connecting the bodies i and j. The deformation of the beam yields force and moment applied on the section S ^ f j that must be equilibrated at section S ^ i f . A first simplified model, hereafter referred to as the linear model, performs the balance in the undeformed configuration and leads to the following expressions
F i f ( S ^ i f ) M i f = A f 0 ( A f T d f ( S ^ i f ) ) T 1 ˜ 1 T f F f j ( S ^ f j ) M f j
where A f is the rotation matrix mapping S ^ f j to S ^ i f and d f is the distance vector between the two sections, respectively, defined as
A f = cos ( θ f ) sin ( θ f ) sin ( θ f ) cos ( θ f ) , d f ( S ^ i f ) = r f sin ( θ f ) r f ( 1 cos ( θ f ) )
In Section 2, the stiffness matrix of the curved beam has been derived considering the undeformed configuration, meaning that the stiffness model is linear and cannot capture the geometrical non-linearity coming from the change of configuration during the beam deformation. Despite this, the balance in the deformed configuration can be modified including the tip displacement due to deformation. Referring to Figure 3, the tip displacement can be included in deriving the moment M i f at the first section S ^ i f , therefore modifying the previous Equation (15) into
F i f ( S ^ i f ) M i f = A f 0 ( A f T d f ( S ^ i f ) + x f ( S ^ f j ) ) T 1 ˜ 1 T f F f j ( S ^ f j ) M f j
where F f j and M f j are referred to the deformed configuration. It is noteworthy that this partial non-linear model is not the geometrically exact fully non-linear model of the curved beam since the forces and moment at section S ^ f j are still obtained using a linear stiffness model for the curved beam.
In the following, either the linear or the partial non-linear model will be included to derive the kinetostatic equations of a planar mechanism. Let us consider the body i in its deformed configuration, as displayed in Figure 4. The flexures have been removed and replaced with their reaction forces and torques where the minus signs come from Newton’s third law. The static balance of body i requires that the following system be satisfied
F i F 1 i F i 2 = 0
M i M 1 i M i 2 + s i 1 T 1 ˜ F 1 i + s i 2 T 1 ˜ F i 2 = 0
where F i and M i are the external force and torque applied to body i, respectively. From the balance Equations (15) and (17), the forces and moments coming from the flexures have been expressed in the undeformed configuration and are now turned into the deformed one. From Figure 4, it can be found that A 1 i F 1 i ( S 1 i ) A ^ 1 i F ^ 1 i ( S ^ 1 i ) and A i 2 F i 2 ( S i 2 ) A ^ i 2 F ^ i 2 ( S ^ i 2 ) , hence it follows that
F i A i A ^ 1 i F ^ 1 i ( S ^ 1 i ) A i A ^ i 2 F ^ i 2 ( S ^ i 2 ) = 0
M i M 1 i M i 2 + s i 1 T 1 ˜ F 1 i + s i 2 T 1 ˜ F i 2 = 0
Considering the frame invariance of the scalar equation of moments, the final system reads
F i A i A ^ 1 i F ^ 1 i ( S ^ 1 i ) A i A ^ i 2 F ^ i 2 ( S ^ i 2 ) = 0
M i M ^ 1 i M ^ i 2 + s ^ i 1 T 1 ˜ A ^ 1 i F ^ 1 i ( S ^ 1 i ) + s ^ i 2 T 1 ˜ A ^ i 2 F ^ i 2 ( S ^ i 2 ) = 0
or in matrix form
F i M i A i A ^ 1 i 0 s ^ i 1 T 1 ˜ A ^ 1 i 1 F ^ 1 i ( S ^ 1 i ) M ^ 1 i A i A ^ i 2 0 s ^ i 2 T 1 ˜ A ^ i 2 1 F ^ i 2 ( S ^ i 2 ) M ^ i 2 = 0
Then, denoting with
N 1 i = A i A ^ 1 i 0 s ^ i 1 T 1 ˜ A ^ 1 i 1 , N i 2 = A i A ^ i 2 0 s ^ i 2 T 1 ˜ A ^ i 2 1
and with R i the residual vector of body i, the final kinetostatic model is
R i w i + N 1 i K 1 ( S ^ 1 i ) ψ 1 ( S ^ 1 i ) + N i 2 T 2 K 2 ( S ^ 2 j ) ψ 2 ( S ^ 2 j ) = 0
where w i = [ F i T , M i ] T is the generalized force vector or wrench acting on body i, K 1 and K 2 are the stiffness matrices of the curved beams connected to the body, and T 2 is the matrix defined in Equation (17). If the latter is replaced by the matrix T 2 of the linear model in Equation (15), the kinetostatic model turns into
R i w i + N 1 i K 1 ( S ^ 1 i ) ψ 1 ( S ^ 1 i ) + N i 2 T 2 K 2 ( S ^ 2 j ) ψ 2 ( S ^ 2 j ) = 0
The kinetostatic equations of a complex multibody system are derived by assembling the residual vectors of all rigid bodies. The final system can be cast in the form
R ( w , q ) = 0 ,
where R indicates the global residual vector, w is the global wrench of external forces and torques, and q is the global vector of generalized coordinates. The elastostatics model offers two types of analyses: the inverse and the direct kinetostatic analysis. The deformed configuration is the input, and the global wrench is the output in the inverse analysis, as has been already described in [36]. The solution of the inverse kinetostatic analysis is straightforward and does not require an iterative procedure. In the direct kinetostatic analysis, the forces and moments applied to the system are known, while the final configuration of the deformed mechanism is sought. This highly non-linear problem can be solved using an iterative procedure such as the Newton–Raphson method described in Algorithm 1.
Algorithm 1 Newton–Raphson iterative method
1:
ϵ given threshold to terminate iterations
2:
k = 1 iteration number
3:
q ( k ) = q ^ set the solution guess value
4:
procedureNewton–Raphson iterative method( q ( k ) , q ^ , w )
5:
     R ( w , q ( k ) ) calculate the residuals as in Equation (24)
6:
     K T ( q ( k ) ) calculate the tangent stiffness matrix as in Equation (26)
7:
     q ( k + 1 ) = q ( k ) K T ( k ) 1 R ( w , q ( k ) ) update the solution
8:
    if  q ( k + 1 ) q ( k ) < ϵ  then
9:
         q = q ( k + 1 ) DKP solution
10:
        exit procedure
11:
    else
12:
         k k + 1
13:
        goto step 4
14:
    end if
15:
end procedure

3.3. Tangent Stiffness Matrix Determination

Suppose that the external forces and moments are fixed in space. Then, considering the residual of the partial non-linear model of Equation (23), the tangent stiffness matrix is
K T i = R i q K T i I + K T i I I + K T i I I I
where
K T i I = N 1 i K 1 ( S ^ 1 i ) ψ 1 ( S ^ 1 i ) q + N i 2 T 2 K 2 ( S ^ 2 j ) ψ 2 ( S ^ 2 j ) q
K T i I I = N 1 i q K 1 ( S ^ 1 i ) ψ 1 ( S ^ 1 i ) + N i 2 q T 2 K 2 ( S ^ 2 j ) ψ 2 ( S ^ 2 j )
K T i I I I = N i 2 T 2 q K 2 ( S ^ 2 j ) ψ 2 ( S ^ 2 j )
If the residual of the linear model of Equation (24) is used instead, the tangent stiffness matrix turns into
K T i = R i q K T i I + K T i I I
The Appendix A reports the expressions for the terms of K T i .

4. Case Study: The Four-Bar Linkage

In this section, a compliant gripper with curved beams, shown in Figure 5a, is studied. The monolithic structure of the mechanism can be reduced to two in-parallel four-bar linkages, as revealed in Figure 5b. For symmetry along the vertical axis, only half a mechanism will be analyzed, i.e., the right side of Figure 5b. The layout of Figure 6 has been plotted following the notation presented in Section 3. Body 1 is the frame, here considered fixed, while bodies 2, 3, and 4 are moving rigid bodies. Considering grippers with CSFH joints, one link is a part of the monolithic structure enclosed between two circular beams. In an MEMS device, the entire monolithic structure can deform. However, the assumption of rigid links coupled to flexible circular beams has been verified using FE models. It is fully justified since the links are at least one order of magnitude stiffer than the circular beam flexures.
First, vectors p ^ f can be calculated knowing the undeformed configuration, i.e., q ^ , therefore
p ^ a = r ^ 2 r ^ 1 + A ^ 2 s ^ 2 a A ^ 1 s ^ 1 a
p ^ b = r ^ 3 r ^ 2 + A ^ 3 s ^ 3 b A ^ 2 s ^ 2 b
p ^ c = r ^ 4 r ^ 3 + A ^ 4 s ^ 4 c A ^ 3 s ^ 3 c
p ^ d = r ^ 1 r ^ 4 + A ^ 1 s ^ 1 d A ^ 4 s ^ 4 d
Similarly, the vectors p f of the flexures in the deformed configuration, i.e., q , are
p a = r 2 r 1 + A 2 s 2 a A 1 s 1 a
p b = r 3 r 2 + A 3 s 3 b A 2 s 2 b
p c = r 4 r 3 + A 4 s 4 c A 3 s 3 c
p d = r 1 r 4 + A 1 s 1 d A 4 s 4 d
Since body 1 is the frame, r 1 = r ^ 1 and A 1 = A ^ 1 . Furthermore, if the frame S 1 of body 1 is coincident with Σ , it follows that r 1 = 0 and A 1 = 1 . Here, these matrices are written for the sake of completeness.
Notice that q could be one of the iterative solutions q ( k ) employed in the Newton–Raphson algorithm. The deformation vectors x f in the material description and expressed in the local frames of the undeformed flexures are
x a ( S ^ a 2 ) = A ^ a 2 T A ^ 2 T ( A ^ 1 A 1 T p a p ^ a )
x b ( S ^ b 3 ) = A ^ b 3 T A ^ 3 T ( A ^ 2 A 2 T p b p ^ b )
x c ( S ^ c 4 ) = A ^ c 4 T A ^ 4 T ( A ^ 3 A 3 T p c p ^ c )
x d ( S ^ d 1 ) = A ^ d 1 T A ^ 1 T ( A ^ 4 A 4 T p d p ^ d )
The angular deformation ϕ f is obtained as
ϕ a = θ 2 θ 1 θ ^ 2 + θ ^ 1
ϕ b = θ 3 θ 2 θ ^ 3 + θ ^ 2
ϕ c = θ 4 θ 3 θ ^ 4 + θ ^ 3
ϕ d = θ 1 θ 4 θ ^ 1 + θ ^ 4
The expressions (31) and (32) allows for determining the flexure generalized deformations ψ a ( S ^ a 2 ) , ψ b ( S ^ b 3 ) , ψ c ( S ^ c 4 ) , and ψ d ( S ^ d 1 ) .
The transformation matrices T f of Equation (17) are defined as
T a = A a 0 ( A a T d a ( S ^ 1 a ) + x a ( S ^ a 2 ) ) T 1 ˜ 0
T b = A b 0 ( A b T d b ( S ^ 2 b ) + x b ( S ^ b 3 ) ) T 1 ˜ 0
T c = A c 0 ( A c T d c ( S ^ 3 c ) + x c ( S ^ c 4 ) ) T 1 ˜ 0
T d = A d 0 ( A d T d d ( S ^ 4 d ) + x d ( S ^ d 1 ) ) T 1 ˜ 0
where A f and d f can be found for each curved beam using Equation (16). Expressions similar to T f , not reported for brevity, can be written to determine T f of Equation (15).
Then, the matrices N of Equation (22) are
N d 1 = A 1 A ^ d 1 0 s ^ 1 d T 1 ˜ A ^ d 1 1 , N 1 a = A 1 A ^ 1 a 0 s ^ 1 a T 1 ˜ A ^ 1 a 1
N a 2 = A 2 A ^ a 2 0 s ^ 2 a T 1 ˜ A ^ a 2 1 , N 2 b = A 2 A ^ 2 b 0 s ^ 2 b T 1 ˜ A ^ 2 b 1
N b 3 = A 3 A ^ b 3 0 s ^ 3 b T 1 ˜ A ^ b 3 1 , N 3 c = A 3 A ^ 3 c 0 s ^ 3 c T 1 ˜ A ^ 3 c 1
N c 4 = A 4 A ^ c 4 0 s ^ 4 c T 1 ˜ A ^ c 4 1 , N 4 d = A 4 A ^ 4 d 0 s ^ 4 d T 1 ˜ A ^ 4 d 1
Setting the external wrenches w i applied at the mass centers G i of the rigid bodies, the four residual vectors R i , i = 1 , , 4 , are
R 1 w 1 + N d 1 K d ( S ^ d 1 ) ψ d ( S ^ d 1 ) + N 1 a T a K a ( S ^ a 2 ) ψ a ( S ^ a 2 )
R 2 w 2 + N a 2 K a ( S ^ a 2 ) ψ a ( S ^ a 2 ) + N 2 b T b K b ( S ^ b 3 ) ψ b ( S ^ b 3 )
R 3 w 3 + N b 3 K b ( S ^ b 3 ) ψ b ( S ^ b 3 ) + N 3 c T c K c ( S ^ c 4 ) ψ c ( S ^ c 4 )
R 4 w 4 + N c 4 K c ( S ^ c 4 ) ψ c ( S ^ c 4 ) + N 4 d T d K d ( S ^ d 1 ) ψ d ( S ^ d 1 )
The residuals are employed to form a system of 12 non-linear kinetostatic equations, i.e.,
R [ R 1 T , R 2 T , R 3 T , R 4 T ] T = 0
that must be solved using an iterative procedure. To calculate the tangent stiffness matrix necessary to apply the Newton–Raphson algorithm, let us define the Jacobians J ˘ f ( S ^ f j ) , i.e.,
J ˘ a ( S ^ a 2 ) [   J ˘ a 1 ( S ^ a 2 )   |   J ˘ a 2 ( S ^ a 2 )   ] = A ^ a 2 T A ^ 2 T 0 0 T 1 J ˘ a
J ˘ b ( S ^ b 3 ) [   J ˘ b 2 ( S ^ b 3 )   |   J ˘ b 3 ( S ^ b 3 )   ] = [ A ^ b 3 T A ^ 3 T 0 0 T 1 ] J ˘ b
J ˘ c ( S ^ c 4 ) [   J ˘ c 3 ( S ^ c 4 )   |   J ˘ c 4 ( S ^ c 4 )   ] = [ A ^ c 4 T A ^ 4 T 0 0 T 1 ] J ˘ c
J ˘ d ( S ^ d 1 ) [   J ˘ d 4 ( S ^ d 1 )   |   J ˘ d 1 ( S ^ d 1 )   ] = [ A ^ d 1 T A ^ 1 T 0 0 T 1 ] J ˘ d
where
J ˘ a = A ^ 1 A 1 T A ^ 1 ( A ¯ 1 T p a 1 ˜ s ¯ 1 a ) 0 T 1 | A ^ 1 A 1 T A ^ 1 A 1 T A ¯ 2 s ¯ 2 a 0 T 1
J ˘ b = A ^ 2 A 2 T A ^ 2 ( A ¯ 2 T p b 1 ˜ s ¯ 2 b ) 0 T 1 | A ^ 2 A 2 T A ^ 2 A 2 T A ¯ 3 s ¯ 3 b 0 T 1
J ˘ c = A ^ 3 A 3 T A ^ 1 ( A ¯ 3 T p c 1 ˜ s ¯ 3 c ) 0 T 1 | A ^ 3 A 3 T A ^ 3 A 3 T A ¯ 4 s ¯ 4 c 0 T 1
J ˘ d = A ^ 4 A 4 T A ^ 4 ( A ¯ 4 T p d 1 ˜ s ¯ 4 d ) 0 T 1 | A ^ 4 A 4 T A ^ 4 A 4 T A ¯ 1 s ¯ 1 d 0 T 1
The Jacobians J ˘ f ( S ^ f j ) can be mapped using Boolean matrices to the final dimension of the system, i.e.,
J a ( S ^ a 2 ) = [   J ˘ a 1 ( S ^ a 2 )   |   J ˘ a 2 ( S ^ a 2 )   |   O   |   O ]
J b ( S ^ b 3 ) = [   O   |   J ˘ b 2 ( S ^ b 3 )   |   J ˘ b 3 ( S ^ b 3 )   |   O   ]
J c ( S ^ c 4 ) = [   O   |   O   |   J ˘ c 3 ( S ^ c 4 )   |   J ˘ c 4 ( S ^ c 4 )   ]
J d ( S ^ d 1 ) = [   J ˘ d 1 ( S ^ d 1 )   |   O   |   O   |   J ˘ d 4 ( S ^ d 1 )   ]
Then, following the expression (A2) of K T i I reported in the Appendix, it yields
K T 1 I = N d 1 K d ( S ^ d 1 ) J d ( S ^ d 1 ) + N 1 a T a K a ( S ^ a 2 ) J a ( S ^ a 2 )
K T 2 I = N a 2 K a ( S ^ a 2 ) J a ( S ^ a 2 ) + N 2 b T b K b ( S ^ b 3 ) J b ( S ^ b 3 )
K T 3 I = N b 3 K b ( S ^ b 3 ) J b ( S ^ b 3 ) + N 3 c T c K c ( S ^ c 4 ) J c ( S ^ c 4 )
K T 4 I = N c 4 K c ( S ^ c 4 ) J c ( S ^ c 4 ) + N 4 d T d K d ( S ^ d 1 ) J d ( S ^ d 1 )
The final expression for the first part of the tangent stiffness matrix is
K T I = K T 1 K T 2 I K T 3 I K T 4 I
The expressions for K T i I I can be obtained starting from z i in Equation (A7), i.e.,
z 1 = G d 1 K d ( S ^ d 1 ) ψ d ( S ^ d 1 ) + G 1 a T a K a ( S ^ a 2 ) ψ a ( S ^ a 2 )
z 2 = G a 2 K a ( S ^ a 2 ) ψ a ( S ^ a 2 ) + G 2 b T b K b ( S ^ b 3 ) ψ b ( S ^ b 3 )
z 3 = G b 3 K b ( S ^ b 3 ) ψ b ( S ^ b 3 ) + G 3 c T c K c ( S ^ c 3 ) ψ c ( S ^ c 4 )
z 4 = G c 4 K c ( S ^ c 4 ) ψ c ( S ^ c 4 ) + G 2 b T d K d ( S ^ d 4 ) ψ d ( S ^ d 4 )
where the matrices G of Equation (A4) are defined as
G d 1 = A ¯ 1 A ^ d 1 0 0 T 1 , G 1 a = A ¯ 1 A ^ 1 a 0 0 T 1
G a 2 = A ¯ 2 A ^ a 2 0 0 T 1 , G 2 b = A ¯ 2 A ^ 2 b 0 0 T 1
G b 3 = A ¯ 3 A ^ b 3 0 0 T 1 , G 3 c = A ¯ 3 A ^ 3 c 0 0 T 1
G c 4 = A ¯ 4 A ^ c 4 0 0 T 1 , G 4 d = A ¯ 4 A ^ 4 d 0 0 T 1
Therefore, K T I I becomes
K T I I = 0 0 z 1 0 0 0 0 0 0 0 0 0 | 0 0 0 0 0 z 2 0 0 0 0 0 0 | 0 0 0 0 0 0 0 0 z 3 0 0 0 | 0 0 0 0 0 0 0 0 0 0 0 z 4
Finally, the third part of the tangent stiffness matrix K T I I I is derived through the 6-dimensional vectors v f of Equation (A11), i.e.,
v a = J ˘ a x T A ^ 2 A ^ a 2 1 ˜ F a ( S ^ a 2 ) J a x ( S ^ a 2 ) T 1 ˜ F a ( S ^ a 2 )
v b = J ˘ b x T A ^ 3 A ^ b 3 1 ˜ F b ( S ^ b 3 ) J b x ( S ^ b 3 ) T 1 ˜ F b ( S ^ b 3 )
v c = J ˘ c x T A ^ 4 A ^ c 4 1 ˜ F c ( S ^ c 4 ) J c x ( S ^ c 4 ) T 1 ˜ F c ( S ^ c 4 )
v d = J ˘ d x T A ^ 1 A ^ d 1 1 ˜ F d ( S ^ d 1 ) J d x ( S ^ d 1 ) T 1 ˜ F d ( S ^ d 1 )
where F f ( S ^ f j ) is obtained taking the force vector from the flexure wrench w f ( S ^ f j ) = K f ( S ^ f j ) ψ f ( S ^ f j ) . The block-matrices J ˘ f x , or J f x ( S ^ f j ) , are derived taking only the first two rows, i.e., the block of x f , from the corresponding Jacobian matrices. Using the equations of (A12), the matrices V f = V fi |Vfj can be obtained and therefore, the matrix K T I I I becomes
K T I I I = N 1 a V a 1 O O N 4 d V d 1 | N 1 a V a 2 N 2 b V b 2 O O | O N 2 b V b 3 N 3 c V c 3 O | O O N 3 c V c 4 N 4 d V d 4
The final expression for the tangent stiffness matrix is K T = K T I + K T I I + K T I I I . If the DOFs of the first body, i.e., the fixed frame, are removed by imposing fixed boundary conditions, the final form of K T will be a ( 9 × 9 ) matrix with the following pattern
K T = | | | | | |

5. Numerical Application

Referring to Figure 6, body 1 is fixed while body 2 is actuated through a vertical force applied at its center of mass. Finally, the end-effector is attached to body 3. Following the layout of Figure 6, all geometric and structural parameters necessary for direct kinetostatic analysis of the case study are reported in Table 1.

5.1. Comparison and Validation

First, let us consider the linear flexure model expressed through Equation (24) of Section 3. Considering an initial actuation force F 2 y = 60 ( μ N), the latter is increased until the x-coordinate of the end-effector becomes zero, i.e., the gripper clamp is completely closed. The gripper deforms as displayed on the left side of Figure 7, wherein the limit cases and the undeformed configuration are reported. Then, let us consider the partial non-linear flexure model expressed through the Equation (23) of Section 3. Performing the same simulations, the workspace becomes that of the right side of Figure 7. The values of F 2 y for which the clamp is closed, respectively, are F 2 y = 63 ( μ N) for the linear model and F 2 y = 77 ( μ N) for the partial non-linear model.
The two models have been compared to a model implemented using the commercial multibody software MSC Adams © . The model includes rigid links and flexible circular beams. The latter are combined to the links’ endpoints through fixed connections. The model has been developed to be comparable with the Matlab model. The only difference pertains to the flexures since the curved beams have been modeled using a two-dimensional geometrical non-linear representation for beam-like structures. Compared to the previous Matlab models, this representation is fully non-linear, as the stiffness matrix of each flexure is updated during deformation. The actuation force has been applied through step functions, and two simulations have been carried out to achieve convergency; one for the forward movement of closing and the other for the backward movement of opening. Figure 8 shows the two simulations and the trajectories accomplished by the end-effector point. As for Figure 7, the deformed configuration has also been plotted to understand the deformation process better.
The three models have been compared in Figure 9 in terms of the end-effector trajectory, also referred to as the mechanism’s workspace (left subplot) and actuation forces (right subplot).
Two materials have been modeled for the curved beams: silicon with Young modulus E y = 100 (GPa) and nylon with Young modulus E y = 3.84 (GPa). Considering the same rectangular cross-section, whose dimensions are reported in Table 1, the silicon bending stiffness is E I = 26.0417 ( μ N μ m 2 ) while the nylon bending stiffness is E I = 1 ( μ N μ m 2 ). The two materials have different properties in terms of elasticity. Silicon has a brittle behavior, while nylon has an anisotropic hyperelastic or visco-hyperelastic behavior. The two bending stiffness values should be seen as extreme cases to test the proposed method. With this premise in mind, the two materials will be assumed to have both isotropic linear behavior, while the different degrees of non-linearity of the models will only concern the geometric stiffness.
The first row of plots in Figure 9 pertains to the silicon while the second one is the nylon. First, let us consider silicon. Observing the top-left subplot, the three arc-shaped trajectories of the end-effector reveal relevant differences only in the final part of the path. The influence of the fully non-linear flexures becomes more evident in the opening movement, where the trajectories become more distant. Compared to the linear model, the top-right subplot reveals that the flexure non-linearities introduce a stiffening effect, and the actuation force required to produce the same displacement grows. It can be observed that the partial non-linear model is stiffer than the fully non-linear model in the forward closing movement and softer in the opening movement. The three models have equal stiffness only at the undeformed configuration where the actuation force is zero.
Then, let us consider nylon. Observing the bottom-left subplot, the end-effector trajectory follows a trend similar to the previous case. The bottom-right plot reveals differences in force range, as could be expected considering the lower bending stiffness of the nylon. Now, the non-linearities are more pronounced, and three inflection points appear for the fully non-linear flexures, which are totally absent in the remaining cases. Despite this, the plot is similar to the simplified cases.
Excluding the limit points of the workspace during the opening movement, the simplified models provide excellent results. Added to this is that the proper workspace is usually limited by other constraints such as the electrical interfaces or the maximum stress in the material, thus making the three models closer than they might appear in Figure 9.
For example, the Conjugate Surfaces Flexure Hinges (CSFH) employed in MEMS micro-grippers have a rotation range limited to ±20 to prevent the silicon from breaking [39,40]. This range is displayed in the opening–closing movement of Figure 9. However, this range is further limited to about ±2 by other phenomena coming from the electrostatic actuation such as sticking-friction anomalies, the pull-in or the impossibility to generating high actuation forces [41].

5.2. Shape Optimization

The simplified models have been employed to perform the cross-section optimization of the curved beams, as shown in Figure 10. These plots can be used in various ways; for example, knowing the maximum actuation force the electrical interface can produce, the section parameters can be chosen to cope with this value. Another example could be related to the choice of the section parameters based on the maximum allowable stress of the material or its fatigue limit. Likewise, the optimization could affect other structural parameters or mechanism lengths.

5.3. Tangent Stiffness Matrix

Since the tangent stiffness matrix is the key element of the direct kinetostatic analysis, in the following, a detailed analysis of the tangent stiffness matrix and its role in the Newton–Raphson algorithm convergence is detailed. As already described in Section 4, after imposing fixed boundary conditions on body 1, the tangent stiffness matrix turns into a 9 × 9 symmetric matrix. Referring to the partial non-linear model, K T has the expression reported in Table 2 in the undeformed configuration.
Let us consider the mechanism in the final deformed configuration obtained by applying F 2 y = 50 ( μ N). In this case, the expression of K T is reported in Table 3, while the percentage difference between the deformed and undeformed case is provided in Table 4. It can be observed that relevant differences appear during the deformation process, especially in the off-diagonal components. Furthermore, the matrix K T is no longer symmetrical in the deformed configuration.
Since K T is composed of three terms, it is legitimate to ask what the contribution of each term is. As reported in Table 5, the first term K T I is the closest to the final expression of K T .
The second term K T I I is reported in Table 6. Remembering the expression for z i in Equation (A7) and the form of K T I I in Equation (A8), observing Table 6, it can be found that at the equilibrium, i.e., if and only if the residuals are zero, the following expression stands
z i = 1 ˜ 0 0 T 0 , w i F i y + F i x 0
For the case study, only body 2, and therefore z 2 , has components different from zero.
Finally, the third term K T I I I is reported in Table 7. It can be noticed that only the components of the inner moments, i.e., due to the flexures, are different from zero. This feature comes from the particular form of V f i , or V f j , in Equation (A12).
To better focus on the importance of the tangent stiffness matrix in achieving solution convergence, let us keep the tangent stiffness matrix constant and equal to that obtained in the undeformed configuration for the entire simulation, i.e., K T = K T ( q ^ ) . The results of Figure 11 reveal that the number of iterations necessary to achieve convergence grows exponentially as the input load increases and no longer converges beyond F 2 y = 76 ( μ N).
The same simulation has been repeated by updating the tangent stiffness matrix following four strategies based on:
  • The complete matrix K T starting each simulation from the undeformed solution;
  • The first term K T I starting each simulation from the undeformed solution;
  • The complete matrix K T starting each simulation from the previous converged solution;
  • The first term K T I starting each simulation from the previous converged solution.
It can be seen that the approaches using the previous converged solution reduce the number of iterations. Similarly, using the complete matrix instead of the first-term approximated matrix results in fewer iterations. The results are displayed in Figure 12.
It is interesting to observe how these trends translate into a computational burden. The reduced number of iterations provided by the strategies based on the previous converged solution translates directly into savings in computation time. The two strategies based on the first term of the tangent stiffness matrix lead to a higher number of iterations but, simultaneously, require a smaller number of variables to be determined and save on calculation times, as shown in Figure 13. The peaks observed in Figure 13 come from memory allocation and other inner processes of Matlab during the first computation. On the other hand, it can be observed from Figure 12 that a higher number of iterations do not correspond to these CPU times. The results have been obtained using an HP workstation equipped with an Intel Xeon CPU @3.20GHz with 32 GB of RAM.
This section is concluded by giving some insights into the computational time obtained using Adams. In Table 8, the CPU times obtained in the opening–closing movement for the Adams fully non-linear model, the Matlab linear model, and the Matlab partial non-linear model are compared. Adams simulations have been performed by disabling the graphic display. The comparison has been carried out for both silicon and nylon. As can be observed, the CPU time of the simplified methods is from 30 to 500 times faster than Adams. It is noteworthy that the simulation time for nylon is ten times faster than silicon for the Adams fully non-linear model.

6. Conclusions and Discussion

The tangent stiffness matrix has been used as a conceptual base to solve the direct kinetostatic problem of planar grippers with curved beams. Two models have been presented to cope with flexure deformations. The first linear model considers the flexure equilibrium in the initial undeformed configuration, while the second partial non-linear model considers the equilibrium in the deformed configuration. Both methods do not include a fully non-linear geometric description of the curved beam flexure whose stiffness matrix is kept constant and equal to that obtained in the undeformed state. The tangent stiffness matrix has been divided into sub-parts to facilitate both the theoretical treatment and the numerical implementation. The linear model led to two sub-parts, while the partial non-linear model introduced a further third sub-part. Both models were tested and compared with a fully non-linear model obtained using the commercial software MSC Adams. The results proved to be in good agreement on most of the mechanism’s workspace, except for the extreme areas wherein the geometric non-linear effects become relevant. The same case study was used to show the method’s potential; for example, in conducting a shape optimization of the flexure cross-section. Finally, the importance of each term of the tangent stiffness matrix in the convergence process was detailed in terms of the number of iterations required to achieve convergence and computational load.
From what has been outlined, the proposed method offers various advantages:
1.
The results of Figure 13 suggest a possible extension to real-time applications of micro and nano-grippers. It is known that the control often requires simplified models to be executed quickly by the control unit. Often these models are obtained by linearizing the equilibrium equations around one or more operating points. Using models with reduced complexity would allow more efficient control strategies such as control in the operating space, inverse dynamics control, pre-calculated torque control. Furthermore, the closed form helps creating more efficient reduced order models [42,43,44].
2.
The tangent stiffness matrix is obtained in closed form. This feature prevents the use of numeric differentiation, making the convergence process of the direct kinetostatic solution more robust. Furthermore, splitting the expression of K T allowed for identifying its most basic terms and calibrating the compromise between the number of iterations and calculation time. The calculation times are considerably reduced by using only the first term of the tangent stiffness matrix and recalculating it at each iteration of the Newton–Raphson algorithm described in Algorithm 1.
3.
The tangent stiffness matrix can be employed to develop a dynamics model to study vibrations. The tangent stiffness matrix is the core of implicit time integration methods primarily employed in flexible multibody dynamics [45]. Shape optimization takes further advantage of the closed form of K T opening scenarios to gradient-based constrained optimization problems based on the kinetostatic analysis.
4.
Both the two simplified models employ curved beams modeled by a constant stiffness matrix. Despite this, the curved beams guarantee finite displacements/rotations in the mechanism, allowing for the expansion of the reachable workspace. The model remains reliable for most of the mechanism’s workspace. The results are accurate in the functional area except for the limit zones of the workspace in which physical constraints usually prevent motion. When the maximum rotations of the curved beams exceed ± 20 the constant stiffness hypothesis can no longer capture the geometric nonlinearities, and the results deviate from the actual case.
5.
Although the proposed method is valid only for planar cases, it can be extended to other compliant mechanisms with constant stiffness flexures without changing the mathematical background.

Author Contributions

Conceptualization, A.C. and N.P.B.; methodology, A.C. and N.P.B.; software, A.C.; validation, A.C.; formal analysis, A.C. and N.P.B.; investigation, P.D.M.; resources, R.S.; data curation, P.D.M.; writing—original draft preparation, A.C.; writing—review and editing, R.S.; visualization, A.C.; supervision, R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The first term K T i I of the tangent stiffness matrix in Equation (26) can be derived by mapping the Jacobian of Equation (13) to the whole n-dimensional vector of generalized coordinates q , therefore obtaining the new ( 3 × n ) extended Jacobian J f
J f = O J ˘ f i J ˘ f j O
where O is the 3 × 3 zero matrix. By substituting Equation (A1) into Equation (27a), the ( 3 × n ) matrix K T i I becomes
K T i I = N 1 i K 1 ( S ^ 1 i ) J 1 ( S ^ 1 i ) + N i 2 T 2 K 2 ( S ^ 2 j ) J 2 ( S ^ 2 j )
with J 1 and J 2 referred to the curved beams connected to body i.
To find the second term K T i I I , let us consider the variation of matrices N 1 i and N i 2 of Equation (22), i.e.,
δ N 1 i = δ θ i G 1 i , δ N i 2 = δ θ i G i 2
where
G 1 i = A ¯ i A ^ 1 i 0 0 T 0 , G i 2 = A ¯ i A ^ i 2 0 0 T 0
From this expression, it follows that
δ N 1 i K 1 ( S ^ 1 i ) ψ 1 ( S ^ 1 i ) + δ N i 2 T 2 K 2 ( S ^ 2 j ) ψ 2 ( S ^ 2 j ) δ θ i z i
where z i is a 3-dimensional vector defined as
z i = N 1 i θ i K 1 ( S ^ 1 i ) ψ 1 ( S ^ 1 i ) + N i 2 θ i T 2 K 2 ( S ^ 2 j ) ψ 2 ( S ^ 2 j )
Using Equation (A3), z i becomes
z i = G 1 i K 1 ( S ^ 1 i ) ψ 1 ( S ^ 1 i ) + G i 2 T 2 K 2 ( S ^ 2 j ) ψ 2 ( S ^ 2 j )
Comparing Equation (A5) with (27b), the ( 3 × n ) matrix K T i I I reads
K T i I I = [ O | | 0   0   z i | | O ]
where z i is mapped in the column corresponding to the angle θ i while all the other entries are zero.
Finally, starting from the expression of T f in Equation (17) and calculating its variation provides
δ T f = O 0 ( δ x f ( S ^ f 2 ) ) T 1 ˜ 0 O 0 ( A ^ f j T A ^ j T J ˘ f x δ q i j ) T 1 ˜ 0
where Equation (10) has been employed to pass from the global to the local system, and J ˘ f x is the block matrix obtained considering the first two rows of J ˘ f pertaining x f only. Multiplying δ T f by the wrench w f j S ^ f j = K f S ^ f j ψ f S ^ f j at section S ^ f j , the following expression is obtained
δ T f w f j ( S ^ f j ) = 0 v f T δ q i j = O 2 , 6 v f T δ q i j = V f δ q i j
where V f is a (3×6) matrix while the 6-dimensional vector v f = [ v f i T , v f j T ] T is defined as
v f = J ˘ f x T A ^ j A ^ f j 1 ˜ F f S ^ f j
being F f ( S ^ f j ) the force composing the wrench w f ( S ^ f j ) = [ F f ( S ^ f j ) T , M f ] T . The matrix V f can also be written as
Vf= V f i | V f j , V f i = 0T0TvfiT,Vfj=0T0TvfjT
Using these expressions, the ( 3 × n ) matrix K T i I I I becomes
K T i I I I = O | | N i 2 V 2 i | | N i 2 V 2 j | | O
in which the matrices are mapped in correspondence with the dofs of the body i and j, respectively.

References

  1. Midha, A.; Norton, T.W.; Howell, L.L. On the nomenclature, classification and abstractions of compliant mechanisms. J. Mech. Des. 1994, 116, 270–279. [Google Scholar] [CrossRef]
  2. Howell, L.L.; Midha, A.; Norton, T.W. Evaluation of Equivalent Spring Stiffness for Use in a Pseudo-Rigid-Body Model of Large-Deflection Compliant Mechanisms. J. Mech. Des. Trans. ASME 1996, 118, 126–131. [Google Scholar] [CrossRef]
  3. Howell, L.L. Compliant Mechanisms; John Wiley & Sons: Hoboken, NJ, USA, 2001; p. 459. [Google Scholar]
  4. Shuib, S.; Ridzwan, M.I.Z.; Kadarman, A.H. Methodology of compliant mechanisms and its current developments in applications: A review. Am. J. Appl. Sci. 2007, 4, 160–167. [Google Scholar] [CrossRef]
  5. Mathew, B.C.; Bharatpatil, V.; Anilchamoli.; Raikwar, M.; Negi, M.S.; Singh, H. Compliant Mechanism and Origami usage in Aerospace and Space Application. In IOP Conference Series: Earth and Environmental Science; IOP Publishing: Bristol, UK, 2021; Volume 775. [Google Scholar]
  6. Kota, S.; Lu, K.J.; Kreiner, Z.; Trease, B.; Arenas, J.; Geiger, J. Design and application of compliant mechanisms for surgical tools. J. Biomech. Eng. 2005, 127, 981–989. [Google Scholar] [CrossRef] [PubMed]
  7. Peterlík, I.; Nouicer, M.; Duriez, C.; Cotin, S.; Kheddar, A. Constraint-based haptic rendering of multirate compliant mechanisms. IEEE Trans. Haptics 2011, 4, 175–187. [Google Scholar] [CrossRef]
  8. Ma, F.; Chen, G. Influence of non-ideal fixed-end constraints on kinetostatic behaviors of compliant bistable mechanisms. Mech. Mach. Theory 2019, 133, 267–277. [Google Scholar] [CrossRef]
  9. Li, G.; Chen, G. A function for characterizing complete kinetostatic behaviors of compliant bistable mechanisms. Mech. Sci. 2014, 5, 67–78. [Google Scholar] [CrossRef]
  10. Han, Q.; Huang, X.; Shao, X. Nonlinear kinetostatic modeling of double-tensural fully-compliant bistable mechanisms. Int. J. Non-Linear Mech. 2017, 93, 41–46. [Google Scholar] [CrossRef]
  11. Chen, G.; Ma, F. Kinetostatic Modeling of Fully Compliant Bistable Mechanisms Using Timoshenko Beam Constraint Model. J. Mech. Des. Trans. ASME 2015, 137, 022301. [Google Scholar] [CrossRef]
  12. Vurchio, F.; Ursi, P.; Buzzin, A.; Veroli, A.; Scorza, A.; Verotti, M.; Sciuto, S.A.; Belfiore, N.P. Grasping and releasing agarose micro beads in water drops. Micromachines 2019, 10, 436. [Google Scholar] [CrossRef]
  13. Hartmann, L.; Zentner, L. A compliant mechanism as rocker arm with spring capability for precision engineering applications. Mech. Mach. Sci. 2015, 30, 1–7. [Google Scholar] [CrossRef]
  14. Kemeny, D.C.; Howell, L.L.; Magleby, S.P. Using compliant mechanisms to improve manufacturability in MEMS. In Proceedings of the ASME Design Engineering Technical Conference, Montreal, ON, Canada, 29 September–2 October 2002; Volume 3, pp. 247–254. [Google Scholar]
  15. Belfiore, N.P.; Bagolini, A.; Rossi, A.; Bocchetta, G.; Vurchio, F.; Crescenzi, R.; Scorza, A.; Bellutti, P.; Sciuto, S.A. Design, fabrication, testing and simulation of a rotary double comb drives actuated microgripper. Micromachines 2021, 12, 1263. [Google Scholar] [CrossRef]
  16. Kota, S. Design of compliant mechanisms: Applications to MEMS. In Proceedings of the SPIE—The International Society for Optical Engineering, San Jose, CA, USA, 28–29 January 1999; Volume 3673, pp. 45–54. [Google Scholar]
  17. Kota, S.; Joo, J.; Li, Z.; Rodgers, S.; Sniegowski, J. Design of compliant mechanisms: Applications to MEMS. Analog. Integr. Circuits Signal Process. 2001, 29, 7–15. [Google Scholar] [CrossRef]
  18. Tolou, N.; Henneken, V.A.; Herder, J.L. Statically balanced compliant micro mechanisms (sb-mems): Concepts and simulation. In Proceedings of the ASME Design Engineering Technical Conference, Montreal, QC, Canada, 15–18 August 2010; Volume 2, pp. 447–454. [Google Scholar]
  19. Xi, X.; Clancy, T.; Wu, X.; Sun, Y.; Liu, X. A MEMS XY-stage integrating compliant mechanism for nanopositioning at sub-nanometer resolution. J. Micromech. Microeng. 2016, 26, 025014. [Google Scholar] [CrossRef]
  20. Guo, T.; Li, J.; Wang, H. Application of compliant mechanism on the polishing robot. Appl. Mech. Mater. 2011, 63–64, 593–596. [Google Scholar] [CrossRef]
  21. Zhan, J.; Zhang, M. Study on edge-forming mechanism for aspheric surface compliant polishing. Zhongguo Jixie Gongcheng/China Mech. Eng. 2013, 24, 1587–1591. [Google Scholar] [CrossRef]
  22. Zhu, W.; Liu, J.; Li, H.; Gu, K. Design and analysis of a compliant polishing manipulator with tensegrity-based parallel mechanism. Aust. J. Mech. Eng. 2021, 19, 414–422. [Google Scholar] [CrossRef]
  23. Verotti, M.; Bagolini, A.; Bellutti, P.; Belfiore, N.P. Design and validation of a single-SOI-wafer 4-DOF crawling microgripper. Micromachines 2019, 10, 376. [Google Scholar] [CrossRef] [PubMed]
  24. Pantano, A.; Minore, G.; Carollo, G. Design of a wiper as compliant mechanisms with a monolithic layout. SAE Int. J. Passeng. Cars—Mech. Syst. 2020, 13, 173–188. [Google Scholar] [CrossRef]
  25. Shinde, S.; Bandi, P.; Detwiler, D.; Tovar, A. Structural Optimization of Thin-Walled Tubular Structures for Progressive Buckling Using Compliant Mechanism Approach. SAE Int. J. Passeng. Cars—Mech. Syst. 2013, 6, 109–120. [Google Scholar] [CrossRef]
  26. Ling, M.; Howell, L.L.; Cao, J.; Chen, G. Kinetostatic and dynamic modeling of flexure-based compliant mechanisms: A survey. Appl. Mech. Rev. 2020, 72, 030802. [Google Scholar] [CrossRef]
  27. Chi, I.T.; Chang, P.L.; Tran, N.; Wang, D.A. Kinetostatic Modeling of Planar Compliant Mechanisms with Flexible Beams, Linear Sliders, Multinary Rigid Links, and Multiple Loops. J. Mech. Robot. 2021, 13, 051002. [Google Scholar] [CrossRef]
  28. Kong, K.; Chen, G.; Hao, G. Kinetostatic modeling and optimization of a novel horizontal-displacement compliant mechanism. J. Mech. Robot. 2019, 11, 064502. [Google Scholar] [CrossRef]
  29. Zhang, C.; Yu, H.; Yang, M.; Chen, S.; Yang, G. Nonlinear kinetostatic modeling and analysis of a large range 3-PPR planar compliant parallel mechanism. Precis. Eng. 2022, 74, 264–277. [Google Scholar] [CrossRef]
  30. Wu, S.; Shao, Z.; Su, H.; Fu, H. An energy-based approach for kinetostatic modeling of general compliant mechanisms. Mech. Mach. Theory 2019, 142, 103588. [Google Scholar] [CrossRef]
  31. Ma, F.; Chen, G. Kinetostatic modeling and characterization of compliant mechanisms containing flexible beams of variable effective length. Mech. Mach. Theory 2020, 147, 103770. [Google Scholar] [CrossRef]
  32. Arredondo-Soto, M.; Cuan-Urquizo, E.; Gómez-Espinosa, A. The compliance matrix method for the kinetostatic analysis of flexure-based compliant parallel mechanisms: Conventions and general force–displacement cases. Mech. Mach. Theory 2022, 168, 104583. [Google Scholar] [CrossRef]
  33. Turkkan, O.; Su, H.J. A general and efficient multiple segment method for kinetostatic analysis of planar compliant mechanisms. Mech. Mach. Theory 2017, 112, 205–217. [Google Scholar] [CrossRef]
  34. Li, J.; Liu, X.; Chen, G. Automatically generating kinetostatic model for planar flexure-based compliant mechanisms. J. Xidian Univ. 2016, 43, 74–79. [Google Scholar]
  35. Li, J.; Chen, G. A general approach for generating kinetostatic models for planar flexure-based compliant mechanisms using matrix representation. Mech. Mach. Theory 2018, 129, 131–147. [Google Scholar] [CrossRef]
  36. Belfiore, N.P.; Simeone, P. Inverse kinetostatic analysis of compliant four-bar linkages. Mech. Mach. Theory 2013, 69, 350–372. [Google Scholar] [CrossRef]
  37. Chung, J.; Hulbert, G. A time integration algorithm for structural dynamics with improved numerical dissipation: The generalized-α method. J. Appl. Mech. 1993, 60, 371–375. [Google Scholar] [CrossRef]
  38. Hilber, H.M.; Hughes, T.J.; Taylor, R.L. Improved numerical dissipation for time integration algorithms in structural dynamics. Earthq. Eng. Struct. Dyn. 1977, 5, 283–292. [Google Scholar] [CrossRef]
  39. Verotti, M.; Crescenzi, R.; Balucani, M.; Belfiore, N.P. MEMS-based conjugate surfaces flexure hinge. J. Mech. Des. 2015, 137, 012301. [Google Scholar] [CrossRef]
  40. Bagolini, A.; Ronchin, S.; Bellutti, P.; Chistè, M.; Verotti, M.; Belfiore, N.P. Fabrication of novel MEMS microgrippers by deep reactive ion etching with metal hard mask. J. Microelectromech. Syst. 2017, 26, 926–934. [Google Scholar] [CrossRef]
  41. Buzzin, A.; Cupo, S.; Giovine, E.; de Cesare, G.; Belfiore, N.P. Compliant nano-pliers as a biomedical tool at the nanoscale: Design, simulation and fabrication. Micromachines 2020, 11, 1087. [Google Scholar] [CrossRef]
  42. Cammarata, A.; Sinatra, R.; Maddìo, P.D. Static condensation method for the reduced dynamic modeling of mechanisms and structures. Arch. Appl. Mech. 2019, 89, 2033–2051. [Google Scholar] [CrossRef]
  43. Cammarata, A.; Sinatra, R.; Maddìo, P.D. A system-based reduction method for spatial deformable multibody systems using global flexible modes. J. Sound Vib. 2021, 504, 116118. [Google Scholar] [CrossRef]
  44. Cammarata, A.; Sinatra, R.; Maddìo, P.D. Interface reduction in flexible multibody systems using the Floating Frame of Reference Formulation. J. Sound Vib. 2022, 523, 116720. [Google Scholar] [CrossRef]
  45. Bauchau, O.A. Flexible Multibody Dynamics; Springer: Dordrecht, The Netherlands; Heidelberg, Germany; London, UK; New York, NY, USA, 2011; Volume 176. [Google Scholar]
Figure 1. Curved beam in its initial (dashed line) and deformed configuration (solid line).
Figure 1. Curved beam in its initial (dashed line) and deformed configuration (solid line).
Micromachines 13 02172 g001
Figure 2. Deformation of a curved beam due to the relative motions of the bodies connected to its extremities.
Figure 2. Deformation of a curved beam due to the relative motions of the bodies connected to its extremities.
Micromachines 13 02172 g002
Figure 3. Static balance of a curved beam. Dashed line for the undeformed beam and a solid line for the deformed beam.
Figure 3. Static balance of a curved beam. Dashed line for the undeformed beam and a solid line for the deformed beam.
Micromachines 13 02172 g003
Figure 4. Kinetostatic balance of a rigid-body.
Figure 4. Kinetostatic balance of a rigid-body.
Micromachines 13 02172 g004
Figure 5. MEMS-based gripper with curved beams. (a) CAD layout. (b) Matlab model: (red) deformed configuration, (blue) undeformed configuration.
Figure 5. MEMS-based gripper with curved beams. (a) CAD layout. (b) Matlab model: (red) deformed configuration, (blue) undeformed configuration.
Micromachines 13 02172 g005
Figure 6. Layout of a four-bar linkage with curved beams.
Figure 6. Layout of a four-bar linkage with curved beams.
Micromachines 13 02172 g006
Figure 7. Workspace of the gripper and end-effector trajectory for the linear and partial non-linear models during the opening/closing maneuver: (red) deformed configuration, (blue) undeformed configuration, (black) end-effector trajectory. Units in millimeters.
Figure 7. Workspace of the gripper and end-effector trajectory for the linear and partial non-linear models during the opening/closing maneuver: (red) deformed configuration, (blue) undeformed configuration, (black) end-effector trajectory. Units in millimeters.
Micromachines 13 02172 g007
Figure 8. Adams results: (left) forward movement; (right) backward movement. The deformed gripper is in red, while the undeformed configuration is represented in blue color. The trajectory of the end-effector for the two movements is displayed in black color.
Figure 8. Adams results: (left) forward movement; (right) backward movement. The deformed gripper is in red, while the undeformed configuration is represented in blue color. The trajectory of the end-effector for the two movements is displayed in black color.
Micromachines 13 02172 g008
Figure 9. Model comparison in terms of workspace and actuation forces. Left subplot: end-effector trajectory in an opening–closing movement; right subplot: actuation force vs. x-coordinate of the end-effector during the opening–closing movement.
Figure 9. Model comparison in terms of workspace and actuation forces. Left subplot: end-effector trajectory in an opening–closing movement; right subplot: actuation force vs. x-coordinate of the end-effector during the opening–closing movement.
Micromachines 13 02172 g009
Figure 10. Cross-section optimization of the flexures.
Figure 10. Cross-section optimization of the flexures.
Micromachines 13 02172 g010
Figure 11. Number of iterations to achieve convergence of the direct kinetostatic analysis varying the input force. The simulation employs the constant tangent stiffness matrix obtained in the undeformed configuration.
Figure 11. Number of iterations to achieve convergence of the direct kinetostatic analysis varying the input force. The simulation employs the constant tangent stiffness matrix obtained in the undeformed configuration.
Micromachines 13 02172 g011
Figure 12. Number of iterations to achieve convergency of the direct kinetostatic analysis varying the input force. The simulation compares four strategies to upload the tangent stiffness matrix: considering the complete matrix K T or only the first term K T I starting each simulation from the undeformed solution, considering the complete matrix K T or only the first term K T I starting each simulation from the previous converged solution.
Figure 12. Number of iterations to achieve convergency of the direct kinetostatic analysis varying the input force. The simulation compares four strategies to upload the tangent stiffness matrix: considering the complete matrix K T or only the first term K T I starting each simulation from the undeformed solution, considering the complete matrix K T or only the first term K T I starting each simulation from the previous converged solution.
Micromachines 13 02172 g012
Figure 13. Computation time of the four strategies of Figure 12.
Figure 13. Computation time of the four strategies of Figure 12.
Micromachines 13 02172 g013
Table 1. Geometric and structural parameters of the case study.
Table 1. Geometric and structural parameters of the case study.
Four-Bar Mechanism
mass center G 2 [1.5000, 0.2000](mm)
mass center G 3 [1.8000, 1.5000](mm)
mass center G 4 [1.0000, 1.5000](mm)
hinge center O a [0.5000, 0.0000](mm)
hinge center O b [1.8478, 0.7654](mm)
hinge center O c [1.1928, 1.9000](mm)
hinge center O d [0.5000, 1.5000](mm)
local base vector s ¯ 1 a [0.6848, 0.0765](mm)
local base vector s ¯ 2 a [−0.8413, −0.3218](mm)
local base vector s ¯ 2 b [0.2139, 0.7140](mm)
local base vector s ¯ 3 b [0.0408, −0.5348](mm)
local base vector s ¯ 3 c [−0.4140, 0.3482](mm)
local base vector s ¯ 3 E E [−1.3571, 1.6991](mm)
local base vector s ¯ 4 c [0.3342, 0.2586](mm)
local base vector s ¯ 4 d [−0.6000, −0.1732](mm)
local base vector s ¯ 1 d [0.3714, 1.6532](mm)
curved beam
beam radius r f , ( f [ a , b , c , d ] ) 0.2(mm)
beam characteristic angle [ θ a , θ b , θ c , θ d ] [300, 320, 330, 250]( )
Young modulus E y 100(GPa)
cross-section base b25( μ m)
cross-section height h5( μ m)
Table 2. Tangent stiffness matrix of the partial non-linear model in the undeformed configuration.
Table 2. Tangent stiffness matrix of the partial non-linear model in the undeformed configuration.
2.40 × 10 3 1.84 × 10 2 −5.42 × 10 2 −1.32 × 10 3 −1.11 × 10 2 −1.01 × 10 3
1.84 × 10 2 2.71 × 10 3 −1.34 × 10 3 −1.11 × 10 2 −1.09 × 10 3 −1.46 × 10 2
−5.42 × 10 2 −1.34 × 10 3 2.29 × 10 3 6.76 × 10 2 −3.29 × 10 2 4.70 × 10 2
−1.32 × 10 3 −1.11 × 10 2 6.76 × 10 2 2.41 × 10 3 2.06 × 10 2 4.98 × 10 2 −1.10 × 10 3 −9.53 × 10 1 4.31 × 10 2
−1.11 × 10 2 −1.09 × 10 3 −3.29 × 10 2 2.06 × 10 2 2.30 × 10 3 −6.44 × 10 2 −9.53 × 10 1 −1.21 × 10 3 −1.75 × 10 2
−1.01 × 10 3 −1.46 × 10 2 4.70 × 10 2 4.98 × 10 2 −6.44 × 10 2 1.51 × 10 3 5.07 × 10 2 7.89 × 10 2 −9.01 × 10 1
−1.10 × 10 3 −9.53 × 10 1 5.07 × 10 2 2.34 × 10 3 −5.58 × 10 1 −3.75 × 10 2
−9.53 × 10 1 −1.21 × 10 3 7.89 × 10 2 −5.58 × 10 1 4.16 × 10 3 −1.08 × 10 3
4.31 × 10 2 −1.75 × 10 2 −9.01 × 10 1 −3.75 × 10 2 −1.08 × 10 3 7.93 × 10 2
Table 3. Tangent stiffness matrix of the partial non-linear model in the final deformed configuration obtained applying F 2 y = 50 ( μ N).
Table 3. Tangent stiffness matrix of the partial non-linear model in the final deformed configuration obtained applying F 2 y = 50 ( μ N).
2.31 × 10 3 −9.15 × 10 1 −1.66 × 10 2 −1.27 × 10 3 −1.45 × 10 2 −9.67 × 10 2
4.22 × 10 2 2.75 × 10 3 −1.20 × 10 3 −1.45 × 10 2 −1.14 × 10 3 −2.65 × 10 2
−5.90 × 10 2 −1.45 × 10 3 2.20 × 10 3 7.24 × 10 2 −2.25 × 10 2 4.81 × 10 2
−1.28 × 10 3 −2.79 × 10 2 6.74 × 10 2 2.36 × 10 3 3.65 × 10 2 5.15 × 10 2 −1.08 × 10 3 −8.63 × 10 1 5.32 × 10 2
6.01 × 10 0 −1.11 × 10 3 −3.28 × 10 2 8.03 × 10 1 2.33 × 10 3 −6.07 × 10 2 −8.63 × 10 1 −1.22 × 10 3 −2.72 × 10 1
−9.59 × 10 2 −3.34 × 10 2 4.46 × 10 2 4.64 × 10 2 −4.42 × 10 2 1.46 × 10 3 4.95 × 10 2 7.76 × 10 2 −2.18 × 10 2
−1.02 × 10 3 2.52 × 10 2 2.50 × 10 2 2.23 × 10 3 −1.41 × 10 3 −1.90 × 10 2
−3.80 × 10 2 −1.20 × 10 3 8.58 × 10 2 6.65 × 10 2 3.92 × 10 3 −1.12 × 10 3
4.43 × 10 2 −1.43 × 10 2 −1.12 × 10 2 −4.97 × 10 2 −9.28 × 10 2 6.76 × 10 2
Table 4. Percentage difference of the tangent stiffness matrices of the partial non-linear model between the final deformed configuration of Table 3 and the undeformed configuration of Table 2.
Table 4. Percentage difference of the tangent stiffness matrices of the partial non-linear model between the final deformed configuration of Table 3 and the undeformed configuration of Table 2.
−3.5−149.8−69.4−3.831.4−3.8
129.61.2−10.431.44.582.0
8.87.8−4.17.1−31.62.3
−3.1152.5−0.4−2.377.63.5−1.2−9.423.4
−105.42.2−0.2−61.01.6−5.6−9.41.1−84.4
−4.6129.3−5.1−6.7−31.4−3.9−2.5−1.7141.6
−7.2−364.2−50.7−4.42424.5−49.4
299.1−0.88.7−1290.9−5.83.7
2.7−18.223.932.7−14.1−14.7
Table 5. The first term K T I in the deformed configuration obtained applying F 2 y = 50 ( μ N).
Table 5. The first term K T I in the deformed configuration obtained applying F 2 y = 50 ( μ N).
2.31 × 10 3 −9.15 × 10 1 −1.16 × 10 2 −1.27 × 10 3 −1.45 × 10 2 −9.67 × 10 2
4.22 × 10 2 2.75 × 10 3 −1.20 × 10 3 −1.45 × 10 2 −1.14 × 10 3 −2.65 × 10 2
−5.74 × 10 2 −1.42 × 10 3 2.19 × 10 3 7.08 × 10 2 −2.46 × 10 2 4.71 × 10 2
−1.28 × 10 3 −2.79 × 10 2 6.74 × 10 2 2.36 × 10 3 3.65 × 10 2 5.15 × 10 2 −1.08 × 10 3 −8.63 × 10 1 5.32 × 10 2
6.01 × 10 0 −1.11 × 10 3 −3.28 × 10 2 8.03 × 10 1 2.33 × 10 3 −6.07 × 10 2 −8.63 × 10 1 −1.22 × 10 3 −2.72 × 10 1
−9.59 × 10 2 −3.34 × 10 2 4.46 × 10 2 4.83 × 10 2 −4.22 × 10 2 1.44 × 10 3 4.76 × 10 2 7.56 × 10 2 −2.15 × 10 2
−1.02 × 10 3 2.52 × 10 2 2.50 × 10 2 2.23 × 10 3 −1.41 × 10 3 −1.90 × 10 2
−3.80 × 10 2 −1.20 × 10 3 8.58 × 10 2 6.65 × 10 2 3.92 × 10 3 −1.12 × 10 3
4.43 × 10 2 −1.43 × 10 2 −1.12 × 10 2 −4.77 × 10 2 −9.10 × 10 2 6.62 × 10 2
Table 6. Second term K T I I in the deformed configuration obtained applying F 2 y = 50 ( μ N).
Table 6. Second term K T I I in the deformed configuration obtained applying F 2 y = 50 ( μ N).
00−5.00 × 10 1 000
00−6.25 × 10 12 000
000000
00000−6.75 × 10 13 000
000004.73 × 10 13 000
000000000
00000−1.95 × 10 11
000007.21 × 10 11
000000
Table 7. Third term K T I I I in the deformed configuration obtained applying F 2 y = 50 ( μ N).
Table 7. Third term K T I I I in the deformed configuration obtained applying F 2 y = 50 ( μ N).
000000
000000
−1.62 × 10 1 −2.14 × 10 1 8.03 × 10 0 1.62 × 10 1 −2.14 × 10 1 −1.03 × 10 1
000000000
000000000
000−1.86 × 10 1 −1.94 × 10 1 1.30 × 10 1 1.86 × 10 1 1.94 × 10 1 −2.29 × 10 0
000000
000000
000−1.99 × 10 1 −1.80 × 10 1 1.45 × 10 1
Table 8. CPU-time comparison in seconds for the opening–closing movement.
Table 8. CPU-time comparison in seconds for the opening–closing movement.
Silicon
openingclosing
Adams fully non-linear model41.0050.00
Matlab linear model0.1500.100
Matlab partial non-linear model0.1500.150
Nylon
openingclosing
Adams Fully non-linear model4.5004.500
Matlab linear model0.1350.100
Matlab partial non-linear model0.1450.135
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cammarata, A.; Maddio, P.D.; Sinatra, R.; Belfiore, N.P. Direct Kinetostatic Analysis of a Gripper with Curved Flexures. Micromachines 2022, 13, 2172. https://doi.org/10.3390/mi13122172

AMA Style

Cammarata A, Maddio PD, Sinatra R, Belfiore NP. Direct Kinetostatic Analysis of a Gripper with Curved Flexures. Micromachines. 2022; 13(12):2172. https://doi.org/10.3390/mi13122172

Chicago/Turabian Style

Cammarata, Alessandro, Pietro Davide Maddio, Rosario Sinatra, and Nicola Pio Belfiore. 2022. "Direct Kinetostatic Analysis of a Gripper with Curved Flexures" Micromachines 13, no. 12: 2172. https://doi.org/10.3390/mi13122172

APA Style

Cammarata, A., Maddio, P. D., Sinatra, R., & Belfiore, N. P. (2022). Direct Kinetostatic Analysis of a Gripper with Curved Flexures. Micromachines, 13(12), 2172. https://doi.org/10.3390/mi13122172

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop