A Novel and Inexpensive Approach for Force Sensing Based on FSR Piezocapacitance Aimed at Hysteresis Error Reduction

Force-sensing resistors (FSRs) are inexpensive alternatives to load cells. They are suitable for applications where noninvasive devices are needed to measure force, stress, or pressure. However, they have been proved to be hysteresis prone and offer nonrepeatable readings due to their highly voltage-dependent electrical resistance. A piezocapacitive effect has been found as an alternative phenomenon that is able to offer force-dependent readings of capacitance with less hysteresis error. Also, this capacitance is not dependent on voltage, which also improves repeatability in force measurements. Since measuring capacitance is more expensive than resistance, the least costly conditioning circuitry is desired. An inexpensive alternative using an LM555 that oscillates depending on capacitance is here presented. Hysteresis and repeatability errors have been reduced for a widespread-used force-sensing resistor brand.


Introduction
Conductive polymer composites (CPCs), also known as quantum tunneling composites (QTCs), are manufactured by randomly dispersing conductive particles within an insulating polymer matrix.They exhibit a considerable decrement of electrical resistance when subjected to mechanical stress [1].They are therefore of extensive use as force-sensing resistors (FSR) in robotic links, biomechanical diagnosis, touchpads, touchscreens, and industrial processes, which are demanding force-sensing accuracy with less space consumption at reduced cost [2].Load cells have been the most linear, accurate, repeatable, and temperature-independent devices that have successfully accomplished the task of converting force to voltage [3,4].Nevertheless, their cost and size have been a drawback when applications have budget and space constraints.FSRs dramatically reduce the space needs for operation as well as they are less invasive and can be easily deployed without impacting greatly the mechanical design at final application.However, they are not as linear and temperature independent as load cells, exhibit large amounts of hysteresis, when dynamic forces are exerted on them, and also drift, when the application is about static forces [5].Moreover, its accuracy depends on the possibility of performing a correct modelling since nonlinearities are present.Thus, research endeavors are aimed at deeply understanding FSRs' behavior and also at complementing their usage with appropriate electronic conditioning circuitry and data processing techniques that could compensate errors and nonlinearities [6][7][8][9].
Hysteresis hinders the applicability of some materials in actuators [10] and also in sensors.The overall performance of FSRs is affected when dynamic forces are measured.For example, considerable linearity can be observed at FSRs by loading forces, but unloading forces are associated to a high nonlinearity and a lack of repeatability when the corresponding deformation-stress curves are studied.A complementary higher-order model for the strain-stress characteristic under unloading forces is presented in this paper in order to explain electrical hysteresis using the sensor electric model [11].
Previous works [12][13][14][15][16] have demonstrated that FSRs based on CPCs also embrace a capacitive effect.Their capacitance can be measured and it varies when the CPC is subjected to mechanical stress.Thus, measuring capacitance instead of resistance is an alternative procedure worth to assess for force sensing since the theoretical approach yields in lower hysteresis errors when capacitance readings are considered, enhancing accuracy for plenty of applications [17,18].Nevertheless, measuring capacitance involves additional considerations that increase the complexity of the conditioning circuitry and, thus, the resulting costs.Thus, also a simple and cost-reduced approach for measuring capacitance is here presented.
This paper is organized as follows: Section 2 proposes an approach to mathematically explain the hysteresis phenomenon in CPCs.Section 3 combines the proposed approach of Section 2 with the electrical model of CPCs when subjected to mechanical stress.Section 4 explains and introduces the piezocapacitive effect and relates it to the previous discussion.Section 5 proposes a novel and inexpensive method to measure FSRs' capacitance.Section 6 compares the traditional with the proposed method to sense forces and exposes the experimental results.Conclusions are stated in Section 7.

Deformation versus Stress Curves: A Proposed Complementary Approach for Unloading Stress
Due to the capability that CPCs or QTCs have to change their electrical resistance when subjected to mechanical stress, its use as force sensors has been considered for many applications in medicine and automation.A sandwiched CPC between two electrodes, known as a traditional sandwich element (TSE) [9], is shown in Figure 1.This is the basic configuration for the Tekscan FlexiForce A201-1 according to the registered patents that explain about the electrical conduction mechanism and the way they are manufactured [19,20].Rubbers, elastomers, and polydimethylsilicone (PDMS) are the preferred materials to compose the nonconductive matrix where the nanoparticles are randomly dispersed in.Nanoparticles, whose dimensions are within the range of 10 −8 and 10 −6 meters, are made of conductive materials such as copper or nickel.They allow electrical conduction between electrodes although thin films of the polymer matrix avoid direct physical contact among them [21][22][23].Let s be the average interparticle distance inside of the CPC, which changes with the applied axial stress due to polymer deformation, so it can be stated that s = s(σ), where σ is the mechanical stress (σ) over the sensor-sensitive area (SSA).Let the stress (σ) also be the ratio between the exerted force (F) and the area (A FSR ) where F is applied on.Loading forces are meant to reduce s from an unloaded initial nondeformed state (s 0 ) to a loaded lower deformed state (s f ) due to axial stress.Unloading forces are applied when the FSR is released until σ = 0.During unloading forces, the FSR recovers its nondeformed state and the average filler interparticle distance returns from s f to s 0 , but this happens through a very different path [24].As shown in Figure 1(a), a straight linear path is established form s 0 to s f as the FSR is loaded by a timeincreasing stress, but the way the interparticle separation recovers to its nonloaded state is clearly nonlinear.This hysteresis phenomenon affects also the electric behavior of the FSR since its electrical resistance depends on the instantaneous value of s [11].
Let ε be the strain ratio between the instantaneous interparticle separation s and s 0 .Equation (1) shows the linear relationship that can be established between deformation and strain in terms of s 0 for loading forces [25].
As expected, there is a linear relation between interparticle separation and stress for loading forces which depends on the average interparticle distance at σ = 0, s 0 , and the material stiffness, M, both as invariant parameters.Nevertheless, for unloading forces, this linear approach is no longer valid as suggested by the aforementioned reasons.The interparticle separation recovers to its nonloaded state through a different nonlinear path.Let (2) be the relationship between s and the strain during unloading forces.
where ε rec and s u are, respectively, the strain and the stress occurring when the CPC is being unloaded.It is worth to remark that σ/M ≠ ε rec and s l ≠ s u due to the hysteretic phenomenon.Equation ( 2) describes the path that s follows to recover to the nonloaded state, from s f to s 0 .It is equivalent to the analysis needed to derive (1), but considering a positive sign since s increments with the strain.When no recovery strain (ε rec ) exists, s = s f and σ = σ f .Therefore, the point in Figure 1(a) (σ f , s f ) is the one where loading stops and recovery path begins from.Equation (3) defines s f in terms of σ f which is the final and maximal stress applied when loading stops (and unloading begins).
The next step to define an approach for establishing a proper relationship between interparticle separation and stress with hysteresis is about defining ε rec .For unloading forces in the presence of hysteresis, it is straightforward to understand that (1) is not enough to represent a nonlinear relationship between strain and stress for two reasons: first, as the Young's modulus is constant for the considered material under any circumstances, a linear relation like (1) for unloading forces would yield no hysteresis at all and second, a different general linear relation considering applying a gain or an offset would be unable to describe a recovery state from s f to s 0 as Figure 1(a) shows.Thus, including an exponent to the expression (2) would ensure it to represent a nonlinear relation.Equation (4) combines the whole analysis discussed up to now.
where a and b are real nonzero positive constants to be defined in a way that (1) and (4) meet in s 0 and s f .The 2 Journal of Sensors constant k is the exponent whose effect in this approach will be explained later in this paper.To define a and b, it is needed to remark what (5) state for.
s l 0 = s u 0 = s 0 , After the respective algebraic operations needed from the consideration of ( 5), ( 6) is obtained.Thus, values for a and b can be determined: Equation (8) already considers that k is an even integer in order to keep real all parameters in the equations.Combining (7) and ( 8) in (4), an expression for s u (σ) can be found: Equation ( 9) is a nonlinear approach to represent deformation at unloading forces.It accomplishes with one of the observed characteristics of hysteresis: consideration of σ f , since the shape of the unloading path varies depending on the stress level reached by the exerted force.Simulations in Figure 2 show how (1) and ( 9) describe a hysteretic pattern when considered for loading and unloading forces, respectively.It comprises different values of σ f and k.
Once the proper relationships to represent deformation on stress with considerable hysteresis are established, the next step of the analysis is about its consideration in the Simmons' model [11].It is worth to remark that the proposed approaches for s l (σ) and s u (σ) (( 1) and ( 9) are not a full model of hysteresis in CPCs.The hysteresis phenomenon is widely more complex than the stated equations, since it comprises a memory effect as already demonstrated.Moreover, it has been observed that several cycles on loading and unloading forces yield in different hysteresis paths for every cycle.That is, if (1) and ( 9) are considered valid and accurate enough for a first cycle, with all likelihood, they will be no longer valid for the next force cycle due to the memory effect and also due to the creep effect present in CPCs [26].For a more accurate approach able to be considered as a model, the exponent k might be considered as a time-and stress-dependent function, but this extension of the approach is beyond the scope of this paper.More thorough modeling (variational formulation, nonlinear finite elements) has been carried out for electroactive polymers (EAPs) not used for sensing though, but as actuators.However, the material EAPs are similar to 3 Journal of Sensors CPCs and has the same electromechanical characteristics [27,28].Also, existing modelling alternatives have not been aimed at acquiring data with less hysteresis but at considering more complete models to compensate hysteresis [29].The proposed approach will be enough to assess the electrical hysteresis as a consequence of hysteresis in the interparticle separation versus stress characteristic.

Consideration of the Proposed Approach in the Simmons' Model for Force Sensing from Conductance Variation
The Simmons' model [11] consists a set of three equations that predict the current density (J) based on s(σ) and the applied voltage u in the CPC, excluding contact resistance between both electrodes and the CPC.Since CPCs are a combination of conductive particles separated by very thin layers of polymer, quantum tunneling effect happens.Electrons can flow through these thin layers which is a potential barrier even if their particle energy is not enough to classically overcome this barrier.Quantum tunneling consists in the existing probability of establishing a considerable amount that overcome a potential barrier, which is the solution of the Schrödinger equation [30].When tunneling occurs, electrical current can be established from particle to particle.Thus, several conduction paths appear between electrodes.As stated in [30,31], the narrower the potential barrier is, the higher the probability of the electrons to overcome the barrier through quantum tunneling is.Hence, FSR's resistance depends on this probability and this is the physical basis of the piezoresistive effect.A model of current density J considering s and u, which E mainly depends on, was developed by Simmons [11] considering the multiple quantum tunneling effects occurring in a CPC.This approach was moved forward for FSRs in [31], where (11) and ( 12) are taken into account.This analysis complements former approaches [6,8,9] where just (10) was embraced.where e is the electron's charge, m its mass, φ 0 is the height of the potential barrier, and the constant ℏ = h/2π, being h the Planck's constant.Equations ( 10) through (12) show two aspects to be assessed.The dependence on s(σ) and on u are as follows: (i) Dependence on s(σ): the shape of J when the FSR is subjected to increasing stress (loading) will be different to the shape when unloading since J u, s l σ ≠ J u, s u σ .This leads to understand that mechanical hysteresis provokes electrical hysteresis.
(ii) Dependence on u: let V be the constant voltage applied by the conditioning circuitry on the FSR.Let I be the current that flows through the FSR and R C the sum of the contact resistance that the interface of the CPC and electrodes provides to the whole resistance of the FSR (see Figure 3).It is straightforward to 8,32,33].Therefore, u is also a function of stress u = u σ and this makes the current density J highly dependent on stress.Hence, the characteristic J versus σ is also significantly affected by hysteresis as Figure 4(a) shows.This has been simulated considering the parameters shown in Table 1.
As well as the J versus σ one, the characteristic resistance per unit area, R A versus σ, where R A = u/J, also presents considerable hysteresis due to the already-studied hysteresis observed in the s versus σ characteristic.Figure 4(b) shows that the behavior of resistance per unit area is as linear as expected.This is desirable to measure pressure or force based on resistance.Moreover, Figure 4(a) also presents the desired linear increment of density current with stress.Current through the whole CPC in axial direction is the product of the effective conductive area A and density current, which also depend on stress.Thus, A = A σ .For comparison purposes, the hysteresis error will be defined in a manner that A(σ) will not affect calculations.Hence, hysteresis error can be defined in terms either of J or I, which is I = I u, s σ = A σ J u, s σ .The formulation for the hysteresis error for any physical quantity is defined by   12) was considered for both simulations carried out with MATLAB R2014.Sensor invariant parameters presented in [31] were considered.They are stated in Table 1.
where Q is the physical quantity to calculate its hysteresis error (HE), Q u is Q at unloading forces, Q l is Q at loading forces, σ d is the stress where the maximal difference between Q u and Q l occurs, and σ f is the already-defined final stress.Figure 5 depicts terms in (13).HE will be considered hereafter in terms of percentage.

Consideration of the Proposed Approach in the Capacitive Model
CPCs present the already described piezoresistive effect due to the quantum tunneling conduction that occurs when conductive particles are separated by thin insulating films.However, this is not the only electrical phenomenon that can be observed in CPCs and specifically in FSRs.Between two conductive particles, which are separated by an insulating material, also, a capacitive effect takes place.Confronted particles are able to accumulate electrical load and to establish an electrical field as shown in Figure 6.The usage of FSRs' based on piezocapacitive effect for several applications has been reported by [13][14][15], and strong evidence of capacitance existence in FSRs has been published in [12] where phase shifting between voltage and current was observed when the FSR was fed with AC power.The combined piezocapacitive effect with quantum tunneling conduction result in the R-C circuit is depicted in Figure 7. Resistances, originated from quantum tunneling conduction, can be combined in parallel if the tunneling conduction occurs between particles that are at the same potential.Then, the resulting resistances happen to be connected in series: vertically placed particles where quantum tunneling occurs successively as a conduction path that exists from electrode to electrode through the whole polymer matrix.These quantum tunneling resistances are connected in series.Thus, a total resistance of the FSR can be obtained but its calculation as a circuit is not straightforward.Resistance per unit area is the result of dividing the voltage in the CPC, u, by the current density J from ( 10) through (12).This yields an expression that depends on s(σ) and as a consequence of quantum tunneling, on u.Similarly, capacitances in the CPC can also be combined in series and in parallel in order to obtain an equivalent capacitance for the FSR.Equation ( 14) is a valid manner of calculating capacitance in piezoresistive sensors [34].For the sake of fair comparison and regarding to the fact that the behavior of the effective conductive area is beyond the scope of this paper, let C A be the capacitance per unit area of the CPC, defined in (15).
where ε 0 is the dielectric permittivity of vacuum and ε r is the relative dielectric permittivity of rubber, that is, the material polymer matrix is composed of.Thus, expressions for the resistance per unit area as well as for the capacitance per unit area are finally obtained in the set from ( 16) through (18) and in (15), respectively.

Journal of Sensors
It is worth to recall that ( 16) through (18) are the result of dividing u by ( 10) through (12).Hence, ( 16) yields in a voltage-independent expression, only valid when u is small.As long as u increases, its effect on the resistance is noticeable and ( 17) and ( 18) are voltage and stress dependent.Figure 8 shows how R A and C A behave when a loading-unloading pattern of forces is exerted on a FSR.For the simulation shown, the applied voltage V is larger than φ 0 /e, so ( 18) is valid to calculate R A .
Numerical values of hysteresis errors of ( 13) can be calculated based on simulations of Figure 8. Table 2 presents the values of HE for the simulations performed in Figures 9(a) and 9(b).As already mentioned, R A is highly dependent on u whereas C A does not depend on u as shown in ( 14) and (15).Therefore, Table 2 shows different HE for R A whereas the HE for C A remains unchanged regardless of the voltage level, u.
In all cases, the HE is lower for capacitance than the corresponding HE for the resistance in the simulation.The reason for this reduction of hysteresis error in capacitance readings relies on the fact that s(σ) is not in the argument of an exponential function neither in (14) nor in (15) since piezocapacitance is a nonquantum phenomenon, but a classic consequence of confronting two conductive plates with a dielectric in the middle.Conversely, s(σ) exists in the argument of exponential functions in ( 16) through ( 18) that describe R A , which exacerbates rheological hysteresis.
Capacitance's expressions are also voltage independent, which could improve repeatability towards changes in voltage (either in V or in u), which are common in FSR, as u = u σ due to the fact that u = V − R C I and R C = R C σ .Therefore, the proposal of measuring capacitance instead of resistance for force sensing seems attractive due to its likely hysteresis error reduction and also its lower repeatability error due to the simulations' results shown in Figure 8 and Table 2. Nevertheless, measuring capacitance is more complicated and expensive than measuring resistance, but the likely error reduction could pay off the cost increment in circuitry.The next section explains how a low-cost circuit used to measure capacitance was designed.

Capacitance Measurement with a LM555 in Modified Astable Operation
Measuring capacitance is not as straightforward as measuring resistance, which can be done with inexpensive and simple circuitry [35,36].Nevertheless, several strategies to measure capacitance have been successfully implemented with the purpose of force sensing as well as what Paredes-Madrid et al. and Damilano et al. have done [12,14].However, the most concerning drawback of those approaches was costs.Although it can be understood that more accurate force reading based on capacitance can payoff the boost in costs, it is also interesting to lower them as much as possible.
A LM555 single timer in astable operation was considered to be an alternative strategy to measure capacitance due to its low cost.Figure 9(a) shows the traditional astable configuration for the LM555.
In astable operation, the output of the depicted circuit in Figure 9(a) is a switched response whose frequency and duty cycle depends on the value of resistances R A and R B and on capacitance C1.Thus, the output frequency f can be related to force as the applied stress modifies sensor's capacitance.Equation (19) shows the relationship between frequency and capacitance in LM555 astable operation whose derivation can be consulted in [37].This analysis is performed considering the periods of charge and discharge of C1.
If the value of R is constant, frequency and capacitance are in an inverse-proportional relationship.Hence, measuring FSR piezocapacitance can be achieved by considering the circuit in Figure 9(a), where C1 is replaced by the FSR itself (capacitance, CFSR, and resistance, RFSR), as depicted in Figure 9(b).However, it must be recalled that the circuital model for the FSR is a resistance (RFSR) and a capacitance (CFSR) in parallel, which cannot be avoided.The method to propose is about measuring piezocapacitance aiming to reduce hysteresis error due its nonquantum electrical behavior.Neither piezoresistance nor piezoimpedance would lead to a less hysteretic response due its quantum resistive component that makes exponential relations appear.The new condition modifies the typical analysis done in [37] and, therefore, (19).Hence, the needed approach to derive the equivalent expression to (19) considering the modification shown in Figure 9(b) is presented next.

When Capacitor CFSR Charges.
According to the analysis performed in [37] but considering that the capacitance C1 is replaced by an FSR, the resulting circuit in the charge period is shown in Figure 10.Thus, the voltage in the capacitor CFSR (V CFSR t ) is derived using the general solution for the obtained RC circuit: that is, the voltage that CFSR reaches when it is fully charged; V CFSR 0 is V CFSR t = 0 , that is, the voltage across CFSR when the charging period has just begun; and R th is the Thèvenin equivalent resistance seen from CFSR.
As also done in the derivation of expression found in [37], V CFSR ∞ is the voltage divider seen on CFSR and V CFSR 0 = Vcc/3.It can be observed that R th is the parallel of 2R and RFSR.So, (20) can be rewritten as Let τ H be the time period that CFSR takes to reach 2/3 Vcc:  The equivalent analysis is followed with the resulting circuit shown in Figure 11.
Considering V CFSR 0 = 2/3 Vcc, V CFSR ∞ = 0, and R th the parallel of R and RFSR, (20) can be rewritten as Let τ L be the time period that CFSR takes to discharge up to Vcc/3 again: The output frequency of the circuit is The duty cycle (DC) of the output signal is Analyzing (23), restrictions for the values of R have to be set in order to ensure τ H to be nonnegative.Natural logarithm's argument has to be Equation ( 30) is the conclusion for the aforesaid restriction about the natural algorithm argument.This means that if (30) is not satisfied, the output of the circuit in Figure 9(b) will not oscillate.Let (31) be the relation between R and RFSR, where N ∈ 4, ∞ , as (30) demands.
The minimal value of the ratio τ H /τ L for N > 4 is 24.37.This means that the duty cycle of the output signal cannot be smaller than 24 37/25 37 = 0 96.An output signal with a duty cycle of at least 96% is unpractical for sensing.This is quite close to the point of no oscillation at all.Physical explanation behind this mathematical approach consists in the fact that τ L is too short due to the highly fast discharge of CFSR.In Figure 11, it can be seen that during the discharge period, CFSR discharges through the paralleled resistances R and RFSR, which is very low and demands considerable amount of current.This is a concerning pitfall for performing force sensing with this method that could be solved by modifying the circuit in Figure 9(b) including a series diode (D1) with R B as Figure 13 shows.
The existence of D1, considered as an ideal diode with no voltage drop, does not modify the analysis done in the charging period and it remains unchanged.Equation (21) though is still valid for the circuit in Figure 13.However, the analysis for the discharging period changes considerably because D1 disconnects R from the equivalent circuit.The existence of D1 modifies the circuit in Figure 11 in a manner that CFSR is discharged only through RFSR and no longer the parallel between R and RFSR, as Figure 14 shows.
Just with the exception of a new time constant, the analysis that led to (24) remains unchanged.Considering the new situation, ( 20) is rewritten as follows: Vcc e −t/ RFSR * CFSR 33 The new value of τ L is Equations ( 27), (28), and (32) change as follows with the inclusion of D1.As analyzed with (32), Figure 15 is the graph of (37), where there is no local minimum beyond N = 4, so almost any value of τ H /τ L can be found.Due to the fact that lim N→∞ τ H /τ L = 2, considering the definition of (37), the smallest duty cycle possible is 66.6% and any value above this can be achieved by selecting the proper value of R.This configuration offers a wider range for the duty cycle and allows to work away from the point where the output of the circuit stops oscillating.Finally, an offset capacitance (C_off) has been added in parallel to the FSR.Since the relationship between frequency and capacitance is inverse, if CFSR happens to be too low, the resulting frequency would be too high.With a minimal capacitance, a maximal frequency could be calculated and thus, a specific range of frequencies can be established.This is important for the logic digital circuit that will measure the output frequency.The resulting circuit is shown in Figure 16.The output is connected to a logic circuit based on Arduino, whose programmed algorithm decodes frequency and duty cycle and also multiplexes either C_off and the FSR itself, regarding the force range that it is intended to measure.It is worth to remark that the presence of RFSR in the circuit cannot be neither eliminated nor diminished since this is a fundamental characteristic of FSRs.A possible concern that the parallel RFSR may cause would be the possibility to bypass CFSR when the applied stress is high, which makes RFSR to be considerably reduced.However, it must be remarked that RFSR is both stress and voltage dependent, and therefore, in order to avoid bypassing CFSR, the sourcing voltage Vcc has to be set to a reduced value so that amid the large applied stress, RFSR remains large enough.Equations ( 16) through (18) show that the lower the applied voltage is, the higher the resistance becomes [38].From Figure 16, the largest voltage connected to the FSR is two-thirds of the sourcing voltage, and thus, RFSR is consequently held within acceptable limits.Given that Vcc = 5 V and considering the highest voltage across the FSR, 2/3Vcc = 3 33 V, there are no concerns about bypassing CFSR.
Once the conditioning circuit has been fully defined, the next testing process made was about applying different forces and registering the output frequencies for each applied force.Thus, a force versus frequency characteristic can be graphed and analyzed.A versatile workbench has been considered in order to gather the data for the forthcoming experiments.A stepper motor was accommodated to exert controlled forces on FSR specimens though a spring.Loading and unloading forces could have been applied in order to acquire the frequency data for several force levels.Detailed information about the workbench can be found in [31,36].
Figure 17(a) shows that a linear fitting for the graphed characteristic is not appropriate, which suggests that this approach is not to be the most adequate for the obtained data.Although data are approximately well distributed on both sides of the line graph, the fitting process yielded a poor goodness of fit (R 2 ) of 0.9432, what seems not attractive since this proposal is about improving error percentages which are below 11% (see Table 2).The reasons that make the linear 11 Journal of Sensors approach not suitable although it was theoretically expected are listed next: (i) Conductive area influence: it is worth to recall that either the FSR total resistance or its total capacitance is the product of R A or C A , respectively, with the effective conductive area, which is the gathering of all the formed conductive paths and depends highly on stress, but is worth to state that its influence reduces linearity to all relationships between physical quantities it multiplies [31].
(ii) RFSR influence: as (35) shows, even if R ≪ RFSR, the dependence on RFSR still exists.Equations ( 16) through (18) are worth to be evaluated in order to recall how nonlinear the relation between resistance and stress is.This remaining dependence on RFSR affects linearity considerably.
However, Figure 17(b) is a third-order polynomial approach that fits better the data.
It is worth to recall that not only that this endeavor of measuring capacitance with the proposed method is mainly for the sake of cost but also that it provides considerable flexibility: It allows determining the frequency range to consider with the proper selection of C_off and high-speed multiplexing.More accurate solutions, such as the ones presented here, offer a more linear characteristic to measure capacitance, but their main drawbacks are cost [12] and flexibility [14].Hence, a cost of inaccuracy is paid if this inexpensive conditioning circuit is considered, but the cost seem not to be too high when the aim of considering capacitance is to lower an error from theoretically less than 11% to an error below 2%.The third-order polynomial approach proposes to pay a cost of computational processing instead of the cost of inaccuracy, which seems to be more attractive considering that linearity error introduced may be greater than the reduction of error to expect with this proposed method.

Experimental Error Comparison
6.1.Hysteresis Error.Four (04) Tekscan FlexiForce A201-1 sensors were considered to be the sample for the experimental procedure.Hysteresis errors were measured using the traditional method and also the proposed capacitive technique.Figure 18 shows the simple conditioning circuit to measure a physical quantity related to the FSR's resistance.In order to avoid working with negative voltage signals, Vcc is set to a negative value.Capacitance CFSR does not play any role since Vcc is direct current.The proposed method considers the circuit already studied and shown in Figure 16.It is worth to recall that the output measurable physical quantity related to force for the proposed method is frequency.Conversely, for the traditional method, the output physical quantity related to force that has to be considered is the output voltage of the amplifier.The relation between RFSR and the output voltage is given as follows: The circuit in Figure 16 was used to acquire force measurements based on capacitance.Loading and unloading forces were applied to four (04) sensors in order to calculate hysteresis error according to (13).The test bench described in [31] was employed again.It is worth to remark that these measurements were done to the same sensors with both methods.The obtained errors are presented in Table 3.The conditioning circuit in Figure 18 was powered up with Vcc = −5 V and its RF = 20 kΩ. Figure 19 shows the readings of frequency and voltage that were acquired with each method.
As expected, in Table 3, it is clear that considering the theoretical approach presented in Sections 2 and 3, the errors for the proposed method are considerably lower than the errors that the traditional method presents.Hysteresis error of the traditional method (HE TM ) is similar to the reported errors calculated in [28].However, according to simulated results in Table 2, two-digit errors were not expected.The following aspects have to be taken in account to explain hysteresis errors above the simulated results: (i) Effective conductive area: simulated results just embrace resistance per unit area as well as capacitance per unit area, not considering the effect of conductive areas.Experimental measurements rely on resistance and capacitance already affected by conductive areas.Thus, a comparison between data in Tables 2 and 3 has to consider that the effect of conductive areas is not taken into account in both analyses.As stated in [31], effective conductive area depends on the applied stress and it has been studied for loading forces.A model of the effective conductive area at unloading forces has not been proposed, but the data presented in this paper suggests that the effective conductive area presents also additional hysteresis, which boosts the total HE.
(ii) Test bench hysteresis: it is also true that the test bench (springs and stepper motor) provides hysteresis beyond the amount predicted at simulations, which may exaggerate the HE as Table 3 shows.However, since both HE TM and HE PM were calculated based on data acquired with the same test bench at the same circumstances, comparison may be carried out.
6.2.Repeatability Error.Repeatability error was also evaluated in this comparison through (39).Two manners of calculating errors were employed.Since repeatability error is about static loading, dead weights were considered instead of the stepper motor, which was appropriate for dynamic loading.The data to calculate errors in both ways is gathered according to the following process: (i) Loading the sensor with a 200 g dead weight to acquire the measurements of voltage or frequency, depending on the method to assess.Then, retire the mass pattern.
(ii) Repeat the previous step five times.A set of five values of voltage or frequency has been collected for the same load applied.
(iii) Repeat the two previous steps with another dead weight to acquire another data set in a different loading condition.
With the acquired data, RE 1 and RE 2 can be calculated: where X Max is the maximal value in the acquired set, X Min is the minimal one, X is the mean value for the same set, X is any value that belongs to the set, and n is the quantity of values that the set has.It is worth to remark that RE 2 is the ratio between the standard deviation of the set and its mean, in percentage.Both errors have been calculated for the same four sensors that the hysteresis error was calculated on.Also, two loading conditions were considered: 200 g and 400 g.In all cases and with all analyses done, repeatability error dramatically decreases with the proposed method in comparison with the traditional one.This was observed regardless of the loading condition (200 g or 400 g) or the error formula indicator ((39)) to consider for comparison.The errors in Table 4 outperform in all circumstances to those on Table 5.This was expectable due to the capacitance's nondependence on voltage.Since piezoresistance in CPCs is a quantum tunneling phenomenon, it is evidently dependent on voltage as theoretically predicted in [11] and demonstrated in [30].As seen in Table 2, different errors are exhibited, for readings based on conductance, depending on the considered voltage to apply.If voltage is not constant, resistance neither is.Therefore, repeatability is affected.This is worth to remark due to the widespread use of voltage dividers as conditioning circuit [8,29,39] that do not guarantee constant voltage on the FSR.Even when the considered circuit for this work (Figure 18) is not a voltage divider but a single OpAmp inverter circuit, the lack of repeatability still appears due to other causes.However, the circuit in Figure 18 is a better solution in comparison with voltage dividers, which is discouraged due to the aforementioned reasons.The proposed method that comprises the circuit in Figure 16 has outperformed even the OpAmp inverter regarding to repeatability error due to the independence of piezocapacitance on voltage.

Conclusions
CPCs exhibit both piezoresistive and piezocapacitive effects due to two factors: quantum tunneling and electrical field confirmation among conductive particles, respectively.Piezoresistive effect has been useful for force sensing due to sensitivity, low invasive profile, and low cost.However, hysteresis phenomenon is observed due to the polymeric behavior of its matrix and this is a concerning drawback worth to be lessen through seeking alternative phenomena in them or modifying the conditioning circuits to gain accuracy although cost increases.Load cells, FSR competitors, are enormously more expensive.Their cost ascends up to USD 500.00 while FSRs' cost roughly reaches not more than USD 25.00.This price gap permits to boost considerably the conditioning circuits' complexity.However, research aims to find out even cheaper alternatives to make them more attractive.Piezocapacitive effect is an alternative phenomenon with the advantage of not having exponential relations that increases hysteresis and therefore, considering piezocapacitance for force sensing is a valid alternative to reduce hysteresis error.Moreover, its nondependence on voltage also lessens repeatability as an additional advantage that contributes to accuracy.Although measuring capacitance is not that straightforward as measuring resistance, it is allowed to boost conditioning circuit complexity due to the aforementioned price gap in comparison to load cells.The proposed method in this paper is the least costly solution in comparison to those found in literature.Current endeavors aim to improve accuracy of force sensing methods to broaden FSR usage at industry.The main goal is to make them measure as good as load cells with considerable cost reduction.

Conflicts of Interest
The authors declare that they have no conflicts of interest.

Figure 1 :
Figure 1: (a) Graph of the interparticle separation (s) versus stress (σ) for a hysteretic conductive polymer composite.Deformation (i) decreases linearly as stress increases and (ii) increases nonlinearly when the CPC is released (unloaded).(b) Scheme of polymer matrix (light gray) with conductive nanoparticles (black) sandwiched between electrodes (dark gray) not subjected to stress, where s 0 depicts the natural average interparticle separation.(c) Scheme of the same polymer matrix subjected to stress (σ f ) showing that the interparticle separation has shrunk to s f as a consequence of stress.

Figure 2 :
Figure 2: (a) Simulated plot for s l (σ) and s u (σ), where different values of σ f are considered.(b) Simulated plot for s l (σ) and s u (σ), where different values of k are considered.MATLAB R2014 was employed to perform both simulations.

Figure 3 :
Figure 3: Electrical model for a FSR.

Figure 4 :
Figure4: Simulation of (a) current density (J) versus stress (σ) at loading forces (blue) and at unloading forces (green).(b) Resistance per unit area (R A ) versus stress (σ) for the same situations.Equation (12) was considered for both simulations carried out with MATLAB R2014.Sensor invariant parameters presented in[31] were considered.They are stated in Table1.

Figure 5 :
Figure 5: Terms to calculate HE regardless of the involved physical quantity.

Figure 7 :Figure 8 :
Figure 7: Resulting RC circuit in a FSR subjected to a potential V.

Figure 13 :Figure 14 :Figure 15 :
Figure 13: LM555 in astable operation with FSR replacing capacitance C1 and diode in series with R B .

Figure 16 :
Figure 16: Final conditioning circuit for force sensing based on capacitance and frequency.

Figure 17 :Figure 18 :
Figure 17: (a) Linear fitting for frequency versus force characteristic with R 2 = 0.9432.(b) Third-order polynomial fitting for frequency versus force characteristic with R 2 = 0.9974.

Figure 19 :
Figure 19: Experimental data taken from (a) the circuit in Figure 16 and (b) the circuit in Figure 18.

Table 1 :
Parameters for simulation in Figure4.
Figure 6: Diagram of FSR where possible electrical fields between particles are depicted.A couple of confronted conductive particles can behave as capacitor.Although not depicted, electrical field between electrodes and neighboring particles also exists.

Table 2 :
Hysteresis error for resistance and capacitance at several levels of u calculated from the simulations of Figure8.

Table 3 :
Hysteresis error for both the traditional and proposed methods.

Table 4 :
Repeatability error for the proposed method based on capacitance measurements.

Table 5 :
Repeatability error for the traditional method based on conductance measurements.