An appropriate choice of contact forces for anthropomorphic robotic grasping devices is important for achieving a balanced grasp. Too little applied force may cause an object to slip or be dropped, and too much applied force may cause damage to delicate objects. Prior methods of grasping force optimization (GFO) in the literature can be difficult to compare due to variability in the parameters, such as the type of grasping device, the object grasped, and the contact model, among other factors. Additionally, methods are typically tested on a very small number of scenarios and may not be as robust in other settings. This paper presents a detailed analysis of three optimization approaches based on the literature, comparing them on the basis of accuracy and computational efficiency. Numerical examples are provided for three types of grasp commonly performed by the human hand (cylindrical grasp, tip grasp, and tripod grasp) using both soft finger (SF) contact and hard finger (HF) contact friction models. For each method and grasping example, an external force is applied to the object in eighteen different directions to provide a more complete picture of the methods' performance. Contact points between the hand and the object are predetermined (given). A comparison of the results showed that the nonlinear and linear matrix inequality (LMI) approaches perform best in terms of accuracy, while the computational efficiency of the linear method is stronger unless the number of contact points and segments becomes too large. In this case, the nonlinear method performs more quickly. Future work will extend the problem of GFO to real-time implementation, and a related work (briefly addressed here) examines the sensitivity of optimization methods to variability in the contact locations.

## Introduction

As robotic grasping systems become more versatile and complex, there is a growing need for robust and adaptable methods allowing these grasping devices to interact with objects. Within the scope of this problem, one important consideration is the development of automatic grasping methods resulting in adaptable and consistently stable grasps. The problem of automatic grasping has been of interest to researchers for several decades. Such considerations are applicable for a wide variety of robotic devices, including industrial robots, humanoid robots, robotic hands, and upper limb prostheses. Additionally, many robotic systems have gained the ability to interact with the environment using multiple contacts (either legs, as in a walking machine, fingers, or multiple limbs), and the principles involved in developing a grasping system are also applicable for other closed kinematic chains. Because a robotic hand's fundamental purpose is to grasp and manipulate objects, a proper choice of grasp is crucial to success in fulfilling tasks.

It is well-recognized that automatic grasping is a complex problem requiring many simplifications for reasonable implementation. Generally, automatic grasping involves the combination of several challenging tasks, including the synthesis of the optimal grasping contact points, load sharing, grasp control, grasp stability, force closure, obstacle avoidance, and task planning, among others.

While all components of the grasp planning and execution process are necessary for effective grasp, this work focuses exclusively on the optimization of grasping forces applied by the hand on the object. Effective force computation is a necessary and important consideration for the manipulation of an object by a multifingered robotic hand; further, it is critical to be able to identify appropriate grasping forces in real time. While the human hand is equipped with a complex and advanced kinematic and sensory structure which allows for natural control of contact forces and manipulation when grasping objects, the problem of choosing appropriate grasping forces for anthropomorphic robotic hands is less straightforward. Grasping force optimization (GFO) is complicated by a variety of factors, including the nonlinearity of constraints and objective functions, the relatively large number of variables and constraints, and the need to compute an optimal solution in real time. Further, some methods are acceptable for real-time computation for a small number of contact points but become unacceptable when the problem becomes more complex. Additionally, the uncertainty associated with contact locations, coefficients of friction, and the perception of object shapes due to occlusion or sensor calibration, among other factors, add complexity to the problem.

The purpose of this research is to perform a more detailed analysis of several possible grasping force optimization approaches as an extension of the work presented in Ref. [1]. The benefits and drawbacks unique to each method are compared to determine which approaches are most appropriate for determining contact forces. Namely, this work includes details regarding the tradeoff between accuracy and computational efficiency and associated factors (such as the number of segments used to approximate the friction constraint) to inform choices about which optimization approach is best. An anthropomorphic, five-fingered model of the hand is used to test the optimization approaches on a series of numerical examples representing grasps commonly performed by the human hand. The optimization approaches used include linear and nonlinear methods as well as the formulation of the friction constraints using linear matrix inequalities (LMIs).

One of the primary barriers to progress in GFO research has been the nonlinearity of the friction constraint represented by Coulomb's law. The most common early approach to GFO was to linearize the problem and solve for the optimal grasping forces using linear programming (LP) methods. While LP methods can be advantageous in terms of computational efficiency, they may provide grasping forces which are either too conservative or which violate the nonlinear friction constraint. Accuracy of LP methods also depends partially on how many segments are used in the linear approximation of the friction constraint. Adding more segments improves the accuracy of the solution but sacrifices computational efficiency; the computational time can become significantly larger for even a moderate number of segments. The most common LP approaches to solving the GFO problem are the simplex method [2–4] and interior point methods [5,6]. Interior point methods typically obtain solutions more quickly than the simplex method for large problems; however, they also require an initial feasible starting point.

Progress in GFO has typically been driven by breakthroughs in problem formulation, where the GFO for multifingered hands was first explored in greater depth by Kerr and Roth [7]. Buss et al. [8] discovered that the friction force constraints and force balance constraints may be represented by the positive definiteness of a certain matrix subject to linear constraints. This approach has been used by many researchers since its development and is beneficial because it retains the nonlinearity of the friction constraints without increasing the computation time unreasonably. Next, Han et al. [9] cast the friction cone constraints as LMIs and formulated the GFO problem as a convex optimization problem. The method presented by Buss et al. [8] also develops appropriate LMIs for both hard finger (HF) and soft finger (SF) friction constraints, in contrast to previous approaches. The combination of LMIs [10] and convex optimization has since been implemented by other researchers [11].

Another popular approach is to split algorithms into online and offline phases [12,13]. The computation time required during online implementation is greatly reduced by solving as much of the problem as possible offline. It is also possible to find a suboptimal but still acceptable solution to reduce computation time. Nonlinear approaches have also been studied because of their increased accuracy over LP methods. A central focus of research on nonlinear approaches is to make them fast enough for online applications. Among the most suitable approaches for real-time control and implementation are neural networks (NN). One approach of this type was developed by Xia et al. [14], where the method is able to converge to an optimal solution with any initial point. Other approaches have been developed for generic problems similar to GFO [15,16], but very few neural networks have been developed to solve the GFO problem specifically.

There are also multiple methods of decomposing contact forces for the improvement of GFO [17,18]. One such method separates the contact force into two components: the manipulation force, which contributes directly to balancing external forces applied to the object, and the internal force, which does not contribute to balancing object forces but which may be adjusted to meet other constraints [19]. Approaches have also been developed to assess the quality of a given grasp. A review of grasp quality measures is provided by Suarez et al. [20].

Currently available GFO approaches can be difficult to compare due to variable factors in the optimization approaches and in the numerical examples (e.g., object type, number of contact points, number of fingers, and contact model). In the previous work [1], the authors presented a basic examination of three methods of grasping force optimization: a linear method, a nonlinear method, and a linear matrix inequality method. As an extension of the work in Ref. [1], this paper presents a more thorough examination of three methods of GFO for anthropomorphic hands and for some grasps which have been identified as common in everyday life. Each method is tested for three different grasping scenarios, where the external force is applied to the object in eighteen different directions. The information provided presents a more complete picture about how each GFO method performs in a variety of settings.

This paper is organized as follows: First, the generic problem formulation for GFO is described. Next, the development of each specific optimization method from Ref. [1] is given. The three numerical examples are shown, along with the directions of applied forces for each example, and the sensitivity of methods with respect to variability in contact locations is briefly addressed. Results and discussion are provided for each force optimization method. Finally, concluding remarks are provided, and future work is discussed.

## Grasping Force Optimization Problem Formulation

The goal of grasping force optimization is to find a set of contact forces between the hand and the object which (1) balances the object weight and other external forces applied to the object; (2) minimizes the force applied to the object by the hand contacts; and (3) meets friction, torque limit, and joint limit constraints. For the following series of equations, parameter dimensions are provided in Table 1 for a hand with *n _{c}* contacts between the hand and the object. In this paper, both HF and SF contact models are used. Table 1 provides dimensions of parameters for both contact models. Figure 1 depicts a visual representation the grasping problem, showing two arbitrary fingers

*j*and

*m*, where ${N}$ represents an inertial frame, ${B}$ represents the coordinate frame at the center of mass of the object, and the contact frames between the hand and the object for $i=1,\u2026,nc$ are represented by ${Ci}$. In ${Ci}$,

*n*is an axis normal to the object, and $ti$ and $oi$ are the two axes tangential to the object. For finger

_{i}*i*with joints 1 through

*k*$qi=[qi,1\cdots qi,ki]T$ is the vector of joint angles, and $\tau i=[\tau i,1\cdots \tau i,ki]T$ is the vector of joint torques and where

_{i},*k*total joints for the hand, and the total vectors for joint angles and joint torques are $q=[q1T\cdots qncT]T$ and $\tau =[\tau 1T\cdots \tau ncT]T$, respectively (note that in Fig. 1, there are

_{1}+…+k_{nc}= k*k*joints associated with finger

_{j}*j*and

*k*joints associated with finger

_{m}*m*). In general, a subscript

*i,n*denotes a quantity associated with the normal axis in ${Ci}$. Similarly, subscripts

*i,t*and

*i,o*denote quantities associated with the tangential axes in ${Ci}$, and a subscript

*i,m*denotes a quantity associated with the moment about the normal axis in ${Ci}$.

The generic optimization formulation is, therefore, written as shown in Eq. (9), where *f* is the objective function to be minimized, and $\u2111$ is used to denote the generalized friction constraint. For this problem, contact points between the hand and the object are predetermined, and appropriate joint angles are found using inverse kinematics as described in Ref. [21]. Because joint angles are known, it is unnecessary to include the joint limit constraints at this stage. The specific formulation of the problem needs adaptations depending on the quantity to be minimized and the type of optimization used. Adaptations of the problem formulation are described in the following section:

find: $\lambda =[\lambda 1T,\u2026,\lambda nT]T$

minimize: *f*

## Optimization Approaches

### Manipulation and Internal Forces.

*G*are the internal forces, $\lambda int$, and can be adapted to satisfy the remaining constraints shown in Eq. (9). Manipulation forces can be solved for directly; thus, the problem can be simplified by solving for the manipulation force and using $\gamma $ as the design variable vector and internal force as the quantity to be minimized. The dimension of $\gamma $ is dependent on the size of the null space basis

### Problem Formulation With Nonlinear Constraints.

The problem formulation which allows for nonlinear friction constraints is shown in Eq. (12), where the objective function *f* represents the square summation of the internal forces. Because the force balance is already satisfied by direct calculation of $\lambda man$, the force balance constraint can be eliminated. The general representation of the friction constraint is $\lambda \u2208\u2111$; that is, the contact wrench must lie within the friction boundaries. For the HF contact, friction constraints are represented directly by Eq. (6). For the SF contact, the friction constraints are represented by Eq. (7). The optimization problem is solved using the *fmincon* function in MATLAB

find: $\gamma $

minimize: $f=\lambda intT\lambda int$

### Problem Formulation With LMIs.

The second problem formulation expresses friction constraints as LMIs as presented in Ref. [4]. Detailed proofs of equivalence to the nonlinear constraints are contained in that paper; the general problem formulation using LMIs is described here.

*P*for HF contact may be represented by LMIs as shown in Eq. (15), and

_{i}*P*for SF contact may be represented by LMIs as in Eq. (16). Each matrix

_{i}*P*may then be combined into one block diagonal matrix,

_{i}*P*, where

*P*matrices are contained along the diagonal as shown in Eq. (17). Equation (18) shows the resulting friction constraint, where 0 is a matrix of all 0 entries with the same dimensions as $P(\lambda )$

_{i}The new problem formulation using LMIs is shown in Eq. (19). Again, the objective function minimizes the square summation of the internal forces. Formulating the constraints both nonlinearly and as LMIs while using the same objective function will allow for a direct comparison of the computational efficiency gained through LMI constraints. Again, the problem is solved using the *fmincon* function in MATLAB

find: $\gamma $

minimize: $f=\lambda intT\lambda int$

### Linear Programming Method.

The simplex method of linear programming has previously been used as a common method for grasping force optimization [2]. In general, the simplex method works by generating a set of all basic feasible solutions, one of which is guaranteed to be an optimum, and systematically searching these solutions to find the optimum. Computational efficiency of the simplex method is sacrificed when the number of basic feasible solutions to be inspected becomes very large. By contrast, interior point methods, some of which have been implemented in grasping force optimization, are computationally efficient but require an initial feasible solution to start the algorithm. In the previous work, several LP algorithms were tested to determine which method performed the most quickly [1]. Based on this comparison, the LP problem in this paper is solved using the *linprog* function in MATLAB with the interior point legacy algorithm.

find: $\gamma $

minimize: $f=cT\gamma $

## Numerical Examples

Kinematics for the hand simulation were developed using the Denavit–Hartenberg method [22]. The dimensions of the palm were chosen to match the average male human hand [23], and the links of phalanges, along with upper and lower limits for each joint angle, were determined based on Ref. [24]. In total, the length of the hand is 18.9 cm and its length is 8.4 cm. Figure 3 shows a general view of the hand with a palm fixed coordinate system attached to the center of the palm which is used to define the direction of force applied to the object.

Three types of grasp which have been identified as important and frequently used in daily life are chosen for numerical examples, tip grasp, tripod grasp, and cylinder grasp [25], as shown in Fig. 4. Tip grasp and tripod grasps make contact only at the fingertips with $nc=2$ and $nc=3$, respectively, whereas cylinder grasp allows each phalanx and the palm to make contact with the object at one point (*n _{c}* = 15).

Varying parameters were assigned to each grasping example to test each optimization method under different conditions. All coefficients of friction were estimated based on ranges measured between the fingertips and various materials in Ref. [26]. For the tip grasp, the object is a cube with side lengths of 2.1 cm. The object mass is set to 0.1 kg with a small coefficient of friction of $\mu i=0.6$ for all contacts. Force is applied to the center of mass of the cube in one of 18 different directions, shown in Fig. 5 and detailed in Table 2 (in Table 2, a + indicates a force applied in the positive direction, and a– indicates a force applied in the negative direction). For the tripod grasp, the object is a cylinder with radius 2.1 cm and length 2.8 cm. The object mass is set to 1 kg with a coefficient of friction of $\mu i=2$ for all contacts. For the cylinder grasp, the object is a cylinder with radius 3.6 cm and length 12.6 cm. The object mass is set to 10 kg with a coefficient of friction of $\mu i=1$ for all contacts. The initial guesses for all cases are $\gamma =0$. Because the linear problem formulation is convex, its solution will provide a global optimal, whereas the nonlinear and LMI problems will yield local optimal solutions.

## Sensitivity Analysis

In this section, a brief overview is provided of a probabilistic sensitivity analysis which was conducted to determine whether the proposed optimization methods remain robust when uncertainty in the positions of the contact locations is considered; the details of the analysis, along with complete results, are given in Ref. [27].

## Results

### Grasping Force Optimization.

For each numerical example, the following comparisons were made: (1) difference in computation time between the HF contact and the SF contact, (2) difference in computation time among the three optimization methods (HF contact), (3) difference in computation time for the linear method with number of segments $ns=4,\u200910,\u200920$, (4) difference in calculated optimal forces between the HF and the SF contact, (5) difference in calculated optimal forces among the three optimization methods (HF contact), and (6) difference in calculated optimal forces for the linear method with $ns=4,\u200910,\u200920$.

For each direction in which force was applied, 100 examples were run with varying amounts of force in the given direction. Average values and standard deviations for computation time for all optimization methods are reported in Table 3.

Several patterns emerge from the data shown in Table 3. First, the nonlinear method consistently solves the problem more quickly than the LMI method. For tip grasp, when there are only two contact points with the object, the difference is small but becomes much larger as the number of contact points increases. The linear method also solves the problem more quickly than the nonlinear and LMI methods in most cases. Exceptions occur for the tripod and cylinder grasps, where the nonlinear method solves more quickly than the linear method using $ns=20$. In all cases, the computation time grows as more contact points are considered, and it grows substantially when whole hand grasp is considered; many of the computation times shown for the cylinder grasp are too slow for practical application. Similarly, in the case of whole hand grasp, the computation time grows quickly as the number of segments used to linearize the friction constraint becomes larger. The difference is so significant that the nonlinear method actually outperforms the linear method when the number of segments becomes too large. Finally, methods which use the HF contact model consistently solve more quickly than the SF model. So, if the optimal forces calculated for the HF and SF contact are similar, it would make sense to use the HF contact model since the computation time is smaller.

Tables 4–7 show some sample results of the forces computed for the tip grasp and cylinder grasp, considering both the HF and the SF contact. Results are presented in Newtons in terms of $\lambda $, the grasping force in the contact frames. In general, for the HF contact, $\lambda $ takes the form $\lambda =[\lambda 1,n\lambda 1,t\lambda 1,o\cdots \lambda nc,n\lambda nc,t\lambda nc,o]T$. For the tip grasp, the object is contacted in two places ($nc=2$). Similarly, results for the SF contact are shown for two methods (the approximation of SF contact for the linear case is not presented here) where $\lambda $ takes the form $\lambda =[\lambda 1,n\lambda 1,t\lambda 1,o\lambda 1,m\cdots \lambda nc,n\lambda nc,t\lambda nc,o\lambda nc,m]T$. In Table 6, only the rows with the significantly large differences between the nonlinear and linear cases are included; all other results are eliminated for brevity, as there are 15 contacts between the hand and the object for cylinder grasp. Similarly, in Table 7, only the rows with the largest differences between the nonlinear and LMI cases are shown. In Tables 6 and 7, the rows with the largest differences between the HF contact and the SF contact for the nonlinear case are shown in bold.

For all optimization methods presented in Tables 4 and 5, the forces meet the expected constraints (normal force must be non-negative, normal and tangential forces must fall within friction constraint boundaries). For tip grasp, the nonlinear and LMI results are identical in both the HF and SF contact cases. The normal forces for the linear methods differ from those for nonlinear and LMI methods but become acceptably close as the number of segments increases from 4 to 20. It is also notable that in the SF case, the additional term which represents the moment about the normal does not contribute much to the applied forces.

For the HF contact, the nonlinear and LMI results remain very close, while for SF contact, certain results differ much more significantly. In general, there is good agreement between the nonlinear and LMI methods about which phalanges of the fingers need to apply force (in most cases, the phalanges for which the forces are very close to zero are the same), but there are differences in the magnitudes of the forces, which typically differ by 2 or 3 N. For the HF contact, once again, there is a good agreement between nonlinear and linear methods about which phalanges need to apply forces, but the differences in magnitudes are in the 4–5 N range. Most notably, the force applied by the middle finger, distal phalanx is substantially different between the nonlinear and linear cases. For most contact points, the results for the three linear cases do not change much as the number of segments changes. For this reason, differences between the linear and nonlinear cases may be attributed to differences in the objective function rather than insufficient segments to approximate the friction constraint.

Tables 8 and 9 show the percent differences among optimization methods for tip and tripod grasps and HF contact, where representative percent differences are shown considering all force directions. Table 8 shows very small differences in the optimal forces calculated using the nonlinear and LMI methods for HF contact, and the differences between the nonlinear and linear methods become smaller as the number of segments increases. When using 20 segments to approximate the friction constraint, all differences are below 1%.

For tripod grasp, the results produced by the nonlinear and LMI methods are very small. In general, the percent difference between nonlinear and linear methods also decreases as more segments are used to approximate the friction cone. The differences are larger for tripod grasp than for tip grasp but are still always below 15%. However, it is also notable that for the middle finger contact point, the percent difference actually increases as more segments are used.

### Probabilistic and Sensitivity Analysis.

A brief summary of the major results of the probabilistic and sensitivity analysis is presented in this section, where complete details of the results may be found in Ref. [27]. Table 10 shows the probabilities of failure for each method and grasp under investigation. The main observation which can be drawn from these results is that the nonlinear and LMI optimization methods remain significantly more robust to variability in contact locations than the linear method; even despite increasing the number of segments used to approximate the friction constraint, the linear method remains far from acceptable for use in grasping force optimization.

## Discussion

In general, the nonlinear method seems to have the best compromise of computation speed and accuracy. Although the linear method performs more quickly for a small number of segments or a small number of fingers, differences occur between the nonlinear and linear methods due to an insufficient number of segments or variability in the objective functions. Also, for a large number of contacts and a large number of segments, the nonlinear method has better accuracy and better computation time compared to the linear method. For whole hand grasp, significantly large differences were also observed between the HF and the SF contact. Despite these differences, it is likely that HF contact is still a better approach due to the difference in computation time between the two methods.

A number of possible challenges exist when implementing GFO methods to control a real robotic hand. One significant factor affecting the control of an artificial hand would be the sensitivity of the presented methods to changes or errors in the contact location. The authors have conducted a detailed examination of this problem, presented in Ref. [27] and briefly discussed in this work, consisting of Monte Carlo simulations to identify both the probability of failure (where “failure” indicates that the constraints of the problem are violated, and the grasp cannot be maintained) and sensitivity levels when variability in the contact locations is considered. Two major results emerged from this analysis. First, the nonlinear and LMI methods both remain robust for all grasps when contact location variability is introduced (probabilities of failure, at a maximum, are 2% and are often significantly less), whereas the linear method, in general, is likely to yield higher probabilities of failure. Second, the probability of failure is likely to be higher, regardless of the optimization method, when the grasp becomes less restrictive, i.e., the highest probabilities of failure occur for tip grasp when the number of contact points between the hand and the object is lowest. The results presented in Ref. [27] complement and strengthen the results of this paper by showing that two of the optimization methods remain viable when variability is introduced; further, combining the results of the analysis presented in this work and the results presented in Ref. [27], the nonlinear method emerges clearly as the most appropriate approach since it requires less computation time than the LMI method. Additional details on the sensitivity of the methods presented here to changes in the contact location may be found in Ref. [27].

Other factors which may affect the control of the real robotic hand include the sensitivity of sensors to measure forces and contact points as well as the requirements for these sensors. It is crucial to identify these and other factors which may affect the implementation of GFO methods. The validation of these results using experimental methods is, therefore, a necessary subsequent step in this work.

Additionally, the numerical examples presented in this paper deal only with static forces. Realistically, an artificial hand must be capable of adapting to variable loads. Future work will include an extension of the methods presented here (focusing primarily on the nonlinear, HF contact method) to find solutions when the external force applied to the object varies. Other methods similar to this formulation will also be explored to improve computational efficiency for online implementation. The nonlinear, HF contact method will also be compared to a neural network to determine which approach is more appropriate. A suboptimal approach or an approach which splits the problem into offline and online components may be beneficial for obtaining solutions more quickly. It would also be desirable to allow the hand to change contact locations (i.e., manipulate an object) after contact has been established. Future work will also explore methods of computing optimal grasping forces combined with object manipulation.

## Conclusion

In this paper, a detailed examination of the three methods to predict optimal grasping forces in Ref. [1] was presented for three common types of grasps performed by anthropomorphic artificial hands. A probabilistic and sensitivity analysis was carried out. In general, the nonlinear method using the HF contact performed with the best combination of accuracy and computation speed, although the linear method was often more efficient except when the number of segments and contacts between the hand and the object became too large. However, in conjunction with the sensitivity analysis considering variability in the contact locations conducted in Ref. [28], the nonlinear method emerges clearly as the most appropriate optimization method since the linear method, in general, does not remain robust to changes in the contact locations. Future work will include extending the methods presented in this paper to cases for which the external load applied to the object varies and comparing them to a neural network. Additionally, methods of combining object manipulation with optimal force calculation will be explored.

## Funding Data

CH Foundation (2012).

NSF Graduate Research Fellowship (1356111).