Abstract
Electronic skin for robotic tactile sensing has been studied extensively over the past years, yet practical applications of electronic skin for the grasping state monitoring during robotic manipulation are still limited. In this study, we present the fabrication and implementation of electronic skin sensor arrays for the detection of unstable grasping. The piezoresistive sensor arrays have the advantages of facile fabrication, fast response, and high reliability. With the tactile data from the sensor array, we propose two quantitative indicators, correlation coefficient and wavelet coefficient, to identify grasping with variable forces and slippage. Those two indicators reflect both time and frequency domain characteristics in the contact forces from the sensor array and can be obtained without large amount of calculation. We demonstrate the utility of this method under various conditions, the results indicate grasping with variable forces, and slippage can be distinguished by this method. The flexible sensor arrays are adopted for tactile sensing on a bionic hand, and the effectiveness of this method in detecting various grasping states has been verified. The electronic skin sensor array and the grasping state monitoring method are promising for applications in robotic dexterous manipulation.
Introduction
Similar as the “sense of touch” to human, tactile sensing is critical for robots to understand and interact with real-world objects and environment.1–5 Tactile sensing becomes an increasingly important capability for robots as their workplace has been shifting from structured industrial sectors to complex and dynamically changing environment.6–10 For instance, robotic dexterous manipulation is a typical scenario where tactile sensation is critical as the touch condition is constantly evolving and the grasping state has to be continuously monitored.11–15 Unstable grasping may lead to catastrophic failure and damage to operated objects, robots, as well as personnel.
Driven by the demand for robotic tactile perception, the development of electronic skin has been attracting significant interest over the past decade.16–20 Inspired by human skin, electronic skin refers to flexible sensors that can detect force,21–25 strain,26–28 temperature, 29 etc. 30 And electronic skin with high sensitivity, spatial resolution, and that can cover large area have been reported.31–36 However, due to the complexity in sensor data acquisition and analysis, especially for sensor arrays with a large number of units, electronic skin sensors actually adopted in robotic manipulation for grasping state monitoring are rather limited.
Flexible sensors that can measure both normal and shear forces have been developed to predict slippage during grasping37–39 ; however, this approach requires the prior knowledge of friction coefficients of contact interfaces, which differ significantly for objects to be operated and can experience variations even on the surface of the same object.40–42 In contrast, slippage detection approaches based on machine learning43–45 and neural network46,47 have been proposed through the analysis of the time and frequency domain features of sensor data. However, those approaches can be computationally costly due to the large amount of data, and the applications are typically limited to the scenarios similar to the training process.
To tackle these challenges, herein we design and fabricate flexible electronic skin sensor arrays and develop a highly effective approach for the monitoring of grasping state based on the sensor data. The piezoresistive sensor array and the data acquisition method can be easily scaled up. Due to the flexibility, the sensor array can be attached to curvy surfaces of robotic end effectors. In the grasping state monitoring approach, we propose two numerical indices by extracting both the time and frequency domain characteristics of the sensor array signals. Unstable grasping states, such as grasping with variable forces and slippage, can be distinguished by the magnitude of those two indices. The effectiveness of this method is demonstrated during various grasping tasks of a bionic hand. Owing to facile fabrication, easy implementation, and low computational cost, the flexible sensor array and the grasping state monitoring method are promising for applications in robotic dexterous manipulation.
Methods and Results
Fabrication of electronic skin sensor array and characterizations
We fabricated the electronic skin sensor array through the development of a flexible and porous piezoresistive film, laying out column and row electrode lines, followed by encapsulation using polymeric thin film. The piezoresistive film is porous thermoplastic polyurethanes matrix (TPU) doped with conductive carbon black (CB) particles. A piece of square film of 25 × 25 × 0.6 mm size was produced for the fabrication of a 5 × 5 sensor array. The intersection of each horizontal and longitudinal electrode line is a piezoresistive sensor unit, as shown in Figure 1a. Details of the fabrication processes can be found in the Supporting Note S1.

We adopted different weight ratios of CB and TPU, from 0.01:1 to 0.4:1, to fabricate the piezoresistive films and tested their force-resistance response. It can be seen from Supplementary Figure S1 that the ratio of percolation thresholds to achieve conductivity in this type of material combination is 0.03:1. The resistance drops rapidly when the ratio of CB exceeds the threshold. The measurement sensitivity of a piezoresistive sensor is defined as
Sodium chloride (NaCl) particles were added to the mixed powder to render the film porous and flexible. After dissolution of the NaCl particles, the porous structure results in lower modulus and larger deformation, and thus, larger change of resistance will occur under the same loading condition, leading to higher measurement sensitivity and high flexibility, as shown in Figure 1b and c. We adopted 1:1 weight ratio of NaCl and TPU powder to achieve the balance between structure integrity and measurement sensitivity. The conductivity of the composites composed of different weight ratios of CB and TPU is shown in Figure 1d. The best measurement performance of the film appears when the ratio of CB, NaCl, TPU is 0.3:1:1, which enables highly linear response and rather high sensitivity of 0.13/N, as shown in Figure 1e and Supplementary Figure S1.
The performances of the piezoresistive sensor were characterized using a unit in the 5 × 5 sensor array. The resistance of the unit varies linearly with applied force. The original resistance is 1400 Ω and drops to 495 Ω when the applied force is 5 N (Fig. 1e). The force-resistance curves of the units in the 5 × 5 sensor array are similar and individually calibrated before applications. The sensor also maintains high stability and repeatability as the resistance keeps almost constant after 2000 load-release cycles as shown in Supplementary Figure S2. Although we only demonstrate the 5 × 5 sensor array, the sensor array can be easily scaled up to large areas with large number of units. It is noted that sensor array can be obtained through facile fabrication schemes, and the flexibility allows easy attachment to even curvy surfaces.
As illustrated in Figure 2a, the equipotential circuit approach is adopted to collect the force data applied on the array. The details of process of data collection are described in Supporting Note S2. It is noted that only m + n wires are needed to cover m × n sensor units using the equipotential method, while 2 m × n connecting wires will be needed through traditional method to connect each electrode, which will significantly increase the complexity of data processing circuit. Figure 2b shows the pressure distribution images when a triangle and a pentagon wood block loaded by an iron cube was placed on the 5 × 5 array, respectively. It is observed that the force map is consistent with object shape with low level of cross talk. And the total forces from the sensor units are 2.241 and 2.303 N, respectively; both agree with the applied force of 2.254 N by the wood piece and iron cube (total weight 230 g), indicating the accuracy of the measurements. Moreover, the cross talk between sensing units in the array is also minimum, as shown in Supplementary Figure S4.

Figure 2c presents the force data on five sensor units when a wood cube was pushed on the sensor array with variable forces without sliding, and sliding against the array. The experiment details are shown in Supplementary Figure S3. It is clear that there is difference in contact forces for the two cases. Typically, large variation in contact forces can occur during both variable forces without sliding and sliding against array; however, it is difficult to distinguish those two conditions based on the magnitude of force data.
Method for monitoring robotic grasping state
Robotic manipulation can be very challenging as various conditions can occur in highly dynamic environment when the manipulator performs grasping tasks.48,49 In a stable grasping process, the force at each contact point remains stable, although tiny fluctuations may happen. The contact force may change as the manipulator follows the command instructions to open or close, and the change can be predicted as it occurs when the operation is executed. In this study, the unstable grasping state considered includes two types of unexpected conditions. The first is that the grasping is subject to external force interferences but the contact between the robot and grasped object is still sustained, and the other case is that the force applied is insufficient to provide enough friction to maintain contact, and slippage occurs at the contact point. Those two unstable grasping events need to be detected in real time and avoided for the safety of operation.
As shown in Figure 2d, when the electronic skin sensor array is adopted, the unstable grasping state typically causes short-term variation on the contact forces. To identify those unstable grasping state, it is necessary to propose quantitative indicators which can comprehensively represent the data changes of all sensor units and directly correlate to grasping state. In this study, we come up with two indicators which reflect both time and frequency domain changes in the contact forces from the sensor array. The complex change in contact forces during unstable state is manifested in the drastic change of the center of force (CoF) and the pressure distribution. In time domain analysis, we apply the calculation of correlation coefficient of the pressure distribution vector as the indicator for monitoring grasping state. In time-frequency domain analysis, we calculate CoF and process it with one-dimensional discrete wavelet transformation (DWT) to obtain the second indicator. The force data from the 5 × 5 array is adopted to illustrate the derivation of those two indicators.
Analysis in time domain
Correlation coefficient of the pressure distribution vector, which represents real-time changes in contact force distribution, is calculated from distribution vectors of the sensor array at two adjacent detection moments. Arrange each row of the 5 × 5 array sensor data collected at t0 in order and convert it into a vector Fdis,0 with 25 elements for correlation calculation. Fdis,0 represents the ratio of the pressure Fi,0 at each unit to the total force Fsum, 0 on the sensor array at moment t0:
Fdis,1 is calculated in the same manner at t1 (t1 = t0 + 100 ms), and the correlation coefficient is calculated based on the two vectors:
where N is the total length of data collected in the 400 ms period. Repeat the above calculation after 100 ms at t2 (t2 = t1 + 100 ms) to obtain Sumcoe,2 and then calculate the relative change:
It is noted that a time segment of 400 ms is adopted to calculate
Analysis in time-frequency domain
It has been known that high frequency signals are generated when unstable grasping such as slippage occurs during robotic manipulation, which may be utilized as the indicator of the grasping state.
50
Fourier transformation is widely used for the analysis of frequency features. However, the time characteristics of the signal is lost in Fourier transformation, that is, the time of occurrence of unstable grasping cannot be obtained using Fourier transformation. In this study, we use DWT to obtain both the temporal and frequency features in the sensor signal. Each sensor unit in the sensor arrays has a virtual coordinate based on the sampling order. We define CoF as the virtual center of the applied forces:
where Fi is the measured force of each unit, Fsum is the total force, xi and yi are the coordinates of each unit, and Xc and Yc are the coordinates of the force center along two dimensional orthogonal directions. Daubechies wavelet is adopted to perform DWT on the CoF data sequence in the period [t1−400 ms, t1], and the approximate coefficient dw1,L and the detail coefficient dw1,H are calculated by Mallat algorithm
51
:
where g and h are the low-pass filter and high-pass filter of Daubechies wavelet function, K is the length of the filter, and N is the total length of the data to be decomposed. The absolute value of detail coefficient represents the magnitude of high frequency component at the moment. During unstable grasping state, there will be drastic change in the CoF and thus large increase in dw1,H. We adopt Vcenter, the square of dw1,H, to enlarge the fluctuation:
and the sum of all the Vcenter values calculated during the [t1−400 ms, t1] period by this part of the data is defined as Sumdw,1;
Repeat the above calculation at t2 (t2 = t1 + 100 ms) to obtain Sumdw,2 and then calculate the relative change:
The relative amplitude change, Ratiodw and Ratiocoe, is used to detect the unstable grasping state. For the stable state, there is little force difference between data at the two periods analyzed; therefore, these two ratios are close to 1. However, in the unstable grasping state, the two ratios will far exceed 1. To focus on the difference between the unstable state and the stable cases, the above two ratios are processed as described below to achieve state indicators Rdw and Rcoe:
R dw and Rcoe are both close to zero during stable grasp and increase during unsteady grasp. It is noted that Rcoe represents the spatial change of grasping force distribution, while Rdw also captures the variation of applied forces in frequency domain during unstable grasping. In addition, the time of large change in Rdw corresponds to the occurrence time of grasping state change due to the adoption of wavelet analysis. Therefore, the combination of Rdw and Rcoe represents a comprehensive approach for unstable state detection as together they take account of the variation in both time and frequency domain.
Demonstration of the algorithm in detecting unstable states
We set up several scenarios to test the effectiveness of the algorithm in the detection of unstable contact between the sensor array and objects. Briefly, a cuboid wood block (size 20 × 20 × 20 mm) was placed on the 5 × 5 sensor array, and then the applied force on the sensor array by the block was adjusted in the range of 4–10 N. It is noted that the block covers at least 15 sensor units on the array. The block was pulled or pushed by a linear motor (LinMot HF01-23) at several constant speeds to create sliding with the array. The contact force and sliding speed were set up at different values for the following cases:
Case 1: Applied force varies from 6 to 10 N without sliding; Case 2:Applied force 8 N, slide once with speed 10 mm/s; Case 3:Applied force 8 N, slide once with speed 3 mm/s; Case 4: Applied force 4 N, slide once with speed 10 mm/s; Case 5: Applied force 4 N, slide once with speed 3 mm/s; Case 6: Applied force 6 N, slide four times back and forth with speed 5 mm/s.
It can be seen from Figure 3a that when the sensor array is subjected to variable pressure but no slippage occurs, the correlation coefficient has relatively large variation due to the large change in magnitude of the applied forces. The amplitude of the wavelet detail coefficient is high, leading to the rather high amplitude of

Data calculated from the sensor array at different scenarios.
In case 2–6, corresponding to Figure 3b–f, it is clear that slippage occurs, and the time of occurrence can be identified from the two indicators (9.4 s for case 2, for instance). Generally, there are two major peaks in the figures of coe,
There are four sliding peaks in the Figures 3f (case 6), indicating that the algorithm remains effective even in the presence of multiple sliding. Based on the results in Figure 3, it is found that variable force without slippage (case 1) leads to relatively lower magnitude of Rdw and Rcoe, and when slippage occurs both the ratios exceed 2.5. This threshold may be used to distinguish variable force condition and slippage.
Detection of unstable grasping by a bionic hand
We further demonstrate the applications of the electronic skin sensor array and the utility of the unstable state detection approach in the manipulation tasks by a bionic hand. To measure the grasping forces, the sensor array was fabricated to fit the size of the hand. The grasping experiment setup is shown in Figure 4a. The system consists of a UR5 robotic arm (Universal Robot Co., LTD), and a bionic hand is installed as the end effector. The bionic hand has five fingers with six degree of freedom for manipulation tasks. During grasping, the areas in contact with objects are mainly from the fingertips of the index, middle, ring, and little fingers and the palm. Therefore, four 1 × 5 sensor arrays were arranged on the four fingertips as shown in Figure 4a and Supplementary Figure S5 to measure the grasping forces, forming a 4 × 5 sensor array. Due to the flexibility, the sensor array could be readily attached to the fingertips. The bionic hand was programmed to go through six actions during the manipulation of a glass of water (Fig. 4b):

Approach the glass;
Grab the glass;
Pick up the glass;
Pour a little water into the glass;
Move the glass horizontally back and forth;
Pour more water until the glass slips.
The details of data analysis are shown in Figure 4c–f and Supplementary Figure S6. The variation of Rdw and Rcoe is shown in Figure 4g. Supplementary Video S1 shows the whole process.
Before grasping, the sensor array is not in contact with the glass in Action 1 and the two ratios maintain 0. A peak occurs in both Rdw and Rcoe at Action 2 due to the force applied when grabbing the glass. The applied forces are not changed during the pickup of the glass; thus, the ratios are 0 at Action 3, indicating a stable state. There is a low peak of Rcoe at Action 4 due to the flow-in of the small amount of water. At Action 5, the horizontal motion of the hand leads to acceleration and deceleration of the glass, resulting in relatively large variation in the handling force (Fig. 4c). The sensor array can effectively detect the fluctuation as there is notable increase in both Rdw and Rcoe. As more water is gradually poured into the glass, the handling force is not enough to overcome gravity at Action 6 and slippage occurs, which leads to a dramatic decrease in handling force and significantly high peaks in both Rdw and Rcoe.
In this manipulation task, the evolution of the two indicators is in line with expectations, verifying the accuracy and effectiveness of detection algorithm. Consistent with the results shown in Figure 3, Action 5 leads to variable handling forces with no relative sliding; thus, the increase in the magnitude of Rdw and Rcoe is rather limited, while slippage causes dramatic peaks in those two indicators in Action 6.
Conclusion
In summary, we have presented the fabrication and implementation of electronic skin sensor arrays for the detection of unstable grasping during robotic manipulation. The piezoresistive sensor arrays have the advantages of facile fabrication, fast response, and high reliability, while the equipotential circuit approach enables highly efficient data acquisition. The sensor array and the data acquisition method can be easily scaled up for a large number of sensing units. To reduce the complexity in tactile data analysis, we propose two quantitative indicators, correlation coefficient and wavelet coefficient, to identify unstable grasping states such as grasping with variable forces and slippage. Those two indicators reflect both time and frequency domain changes in the contact forces from the sensor array and can be obtained without large amount of calculation. With a low-cost microprocessor for signal analysis, the response time for unstable state monitoring can be as short as 100 ms in this study, which can be compared favorably with previous studies using data-intensive deep neural networks or sophisticated data acquisition systems, as shown in Supplementary Table S1.
We demonstrate the utility of this method under various unstable grasping conditions, and the results indicate that variable grasping force and slippage can be distinguished by this method. It is noted that the sensor array can be facilely fabricated; both the sensor array and data acquisition can be easily scaled up. Due to the low computation cost, the grasping state detection method demonstrates satisfactory performances with a low-cost microprocessor and no advanced electronic components. With dedicated data acquisition and analysis system, such as the embedded system, the efficiency of detection can be significantly improved. The flexible sensor arrays are adopted for tactile sensing on a bionic hand, and the effectiveness of this method in detecting various grasping states has been verified during manipulation tasks. The flexible electronic skin and the method for detecting unstable grasping have great potential for practical applications in robotic manipulation.
Footnotes
Acknowledgment
The authors acknowledge the support from Flexible Electronics Research Center of HUST for providing experiment facility.
Author Disclosure Statement
No competing financial interests exist.
Funding Information
We acknowledge the support from National Natural Science Foundation of China (nos. U2013213, 51820105008, 92048302), the Technology Innovation Project of Hubei Province of China under grant 2019AEA171.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
