 Research Article
 Open Access
 Published:
Characterization of variablesensitivity force sensor using stiffness change of shapememory polymer based on temperature
ROBOMECH Journal volume 8, Article number: 24 (2021)
Abstract
In the present study, we propose a variablesensitivity force sensor using a shapememory polymer (SMP), the stiffness of which varies according to the temperature. Since the measurement range and sensitivity can be changed, it is not necessary to replace the force sensor to match the measurement target. Shapememory polymers are often described as twophase structures comprising a lowertemperature “glassy” hard phase and a highertemperature “rubbery” soft phase. The relationship between the applied force and the deformation of the SMP changes depending on the temperature. The proposed sensor consists of strain gauges bonded to an SMP bending beam and senses the applied force by measuring the strain. Therefore, the force measurement range and the sensitivity can be changed according to the temperature. In our previous study, we found that a sensor with one strain gauge and a steel plate had a small error and a large sensitivity range. Therefore, in the present study, we miniaturize this type of sensor. Moreover, in order to describe the viscoelastic behavior more accurately, we propose a transfer function using a generalized Maxwell model. We verify the proposed model experimentally and estimated the parameters by system identification. In addition, we realize miniaturization of the sensor and achieve the same performance as in our previous study. It is shown that the proposed transfer function can capture the viscoelastic behavior of the proposed SMP sensor quite well.
Introduction
Force sensors have been applied to various fields and are required to measure wider load ranges. One example in industry is a manufacturing system that has the flexibility to cope with various kinds of smallquantity production referred to as a flexible manufacturing system. Moreover, in rapidly aging societies, robotic technology has been applied to various fields, including industrial fields as well as nursing and welfare fields [1]. In these applications, widerange force sensors that can obtain load information can measure multiple biosignals, such as heart rate, respiration cycle, and weight transitions [2]. Most force sensors transform the mechanical deformation of the detection area under an applied force into a change in resistance, capacitance, or reflectance that can be measured using electric signals. For example, some force sensors consist of strain gauges bonded to a bending beam. However, with this approach, it is difficult to change the measurement range or sensitivity of a sensor, both of which depend on the material used, the type of strain gauge, and the measurement method. The deformation range depends on the sensor material, and it is difficult to change these specifications after the sensor is produced. For this reason, we previously developed a force sensor using a shapememory polymer (SMP), the measurement range and sensitivity of which can be changed [3, 4].
Shapememory polymers [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19] are increasingly being investigated as smart materials and are used in various fields. Shapememory polymers change their modulus around the glass transition temperature (T_{g}), and are often described as twophase structures comprising a lowertemperature “glassy” hard phase and a highertemperature “rubbery” soft phase. The hard and soft phases are characterized by two different elastic moduli: an elastic modulus for the lowertemperature, higherstiffness glassy plateau and an elastic modulus for the highertemperature, lowerstiffness rubbery plateau. The reversible change in the elastic modulus between the glassy and rubbery states of SMPs can be as high as several 100fold.
Since the stiffness of the SMP can be changed according to the temperature, the measurable force range determined based on the above strain range can also be changed. Moreover, even if the strain resolution is the same, the force resolution can be changed in a similar manner. In this way, the measurement range and sensitivity of the force sensor can be changed according to the temperature. Generally, environmental stability is a very important requirement for sensors. For example, special compensating elements are often incorporated either directly into sensors or into signal conditioning circuits in order to compensate for temperature errors [20]. Note that, inverting the above idea, the proposed sensor uses the temperaturedependent changes positively.
In previous studies [3, 4], we made several prototypes of this sensor by attaching a strain gauge to an SMP sheet with an embedded electrical heating wire and evaluated their basic characteristics. Through experiments with these prototypes, which use the stiffness change of the SMP based on the temperature, we showed that the measurement range and sensitivity can be changed without replacing the actual sensor [3]. On the other hand, the changes in measurement range and sensitivity (ranging from a 100fold to a thousandfold) depend on the Young’s modulus of the SMP and are not adjustable. However, the change may be too large for some applications. Therefore, we also affixed a thin steel plate (see “Basic concept of force sensor” section for additional details). Furthermore, we could reduce the influence of the difference in elastic modulus between the strain gauges and the SMP sheet. This made it possible to reduce the discrepancy between the theoretical values and the measured values. Moreover, SMP force sensors with either one or two strain gauges and steel plates were fabricated, and their accuracy and sensitivity were investigated under the same conditions [4]. Experiments using the prototypes demonstrated that a sensor with one steel plate had a small error and a large sensitivity range, although the dimensions of the sensors were not optimized.
Therefore, in the present study, we miniaturized this type of sensor for practical use. A prototype of this sensor was made by attaching a strain gauge to an SMP sheet with an embedded electrical heating wire, and we evaluated the basic characteristics of the prototype sensor.
In a previous study, we proposed a theoretical formula that took the viscosity of the SMP into consideration, which made it possible to reduce the effect of stress relaxation [4]. However, since the first derivative of the measurement values was used to estimate the force (see “Basic concept of force sensor” and “Generalized Maxwell model” sections for additional details), the measurement errors could be large. Therefore, in the present study, we propose a transfer function method using a generalized Maxwell model. We verified the proposed model experimentally and estimated the parameters by system identification.
Basic concept of force sensor
As shown in Fig. 1, the force sensor was fabricated by bonding an SMP sheet and a steel plate and attaching a strain gauge to the steel. If the viscosity of the SMP and the embedded wire are negligible, by assuming a composite beam consisting of SMP and steel sheets, then the strain on the strain gauge (ε) can be expressed as follows [3, 4, 21]:
where x is the distance between the strain gauge and the position at which the force is applied, E_{p} and E_{s} are the elastic moduli of the SMP and steel, respectively, and I_{p}' and I_{s}' are the area moments of inertia for the SMP and steel, respectively, about the neutral axis of the composite beam. Here, I_{p}' and I_{s}' are expressed as follows:
where b is the width of the beam, h_{p} and h_{s} are the thicknesses of the SMP and steel plates, respectively, and h_{1} is the distance between the neutral axis and the SMP surface and is expressed as follows:
As described in “Introduction” section, E_{p} can be changed according to the temperature. Therefore, the relationship between ε and W given in Eq. (1) can also be changed. Therefore, the change in the measurement range based on temperature can be modified. Moreover, as shown in Eq. (1), by changing the thickness of the steel plate, the measurement range and the sensitivity of the sensor can be modified.
In order to understand the mechanisms producing such unique properties and to design products including SMPs, various mathematical models have been proposed [5, 6, 12,13,14,15,16,17,18,19, 22]. To investigate relaxation processes in polymers, combinations of elements, including springs and dashpots, are widely used for modeling under isothermal conditions [12,13,14,15,16,17,18, 22]. For example, Tobushi et al. proposed a linear constitutive model by modifying a threeelement viscoelastic model combining two springs and a dashpot to represent the deformation characteristics of the SMP [12]. Similarly, in our previous study [4], considering the viscosity of the SMP, we assumed that the relationship between W and ε is given by
where L, M, and N are constants.
Miniaturization of sensor
Design of miniaturized sensor
The dimensions of the prototype SMP sensors used in our previous study [4] and the present study are shown in Table 1. Here, l' is the length of the sensor. The dimensions of the prototype sensor constructed in the present study are similar to those of the commercial force sensor (LVS2KA (T < T_{g}), LVS200GA (T > T_{g}), Kyowa Electronic Instruments Co., Ltd.) used in the experiment. The volume of the proposed sensor was reduced by 87% compared with that in our previous study. In our previous study, we applied a deformation of 5 mm to the tip of the sensor. In the present study, we determined b and h_{p} so that the reaction force W is the rated capacity of a commercial load cell (20 N) when the deformation is 1 mm below T_{g}. The ratio of the sensor length (14 mm) to the deformation (1 mm) is approximately the same as that in the previous study, because the large creep strain of the SMP below T_{g} is not recovered [7]. The thickness of the steel plate (h_{s} = 0.07 mm) is similar to that in our previous study [4]. The relationship between the applied force W and the deflection y is expressed as follows [21]:
where l (= 11 mm) is the distance between the fixed end and the position at which the force is applied. In the present study, we chose a polyurethane SMP (SMP Technologies Inc., MP4510, T_{g} = 45 °C). The fundamental characteristics of this material taken from the product catalogue are listed in Table 2. Substituting E_{p} = 1,350 MPa (T < T_{g}) and E_{s} = 193 GPa into this equation, we obtain W = 27 N.
Prototype
In this study, we prepared two prototype sensors (samples i and ii) to check the reproducibility of the sensor fabrication process. The prototype SMP force sensor is shown in Fig. 2. In the present study, we prepared an SMP sheet with an embedded electrical heating wire in a manner similar to that described in our previous studies [4, 8]. The shape and dimensions of the heating wire (nichrome, outer diameter: 0.26 mm, electrical resistivity: 108 × 10^{−6} Ω cm) are shown in Fig. 3. The underlined lengths in the figure were smaller than those in our previous study.
We bonded the SMP sheet and steel plate (SUS304H, thickness: 0.07 mm) using an adhesive (PPX, Cemedine Co., Ltd.). We attached one strain gauge (KFGS2120C116 L1M2R, Kyowa Electronic Instruments Co., Ltd.) to the steel plate and measured the strain on the surface of the steel plate. The distance between the strain gauge position and the position at which the force was applied (x) was 6 mm. We used a cyanoacrylate adhesive (CC36, Kyowa Electronic Instruments Co. Ltd., operating temperature range: − 30 to 100 °C). The size and electrical resistivity of the prototype sensors are listed in Table 3. Since we fabricated the sensors manually, there were several manufacturing errors. In future studies, it will be necessary to minimize these errors by using a jig and a manufacturing machine.
A schematic diagram of the temperature control system is shown in Fig. 4a. The heating wire was connected to a power supply (PE181.3AT, KENWOOD) with a stable direct current, and a voltage of 5 V was supplied to the prototype sensor. With a thermocouple attached to the surface of the SMP sheet (red circle in Fig. 2), we heated the sheet and maintained a temperature of 70 °C (above T_{g}), as in our previous study [4]. Note that we used only two temperatures above and below T_{g} (70 °C and room temperature, respectively), and maintained these temperatures in order to minimize creep and thermal expansion. A pulse width modulation (PWM) signal was controlled by software developed in LabVIEW (National Instruments Co.). When the SMP was heated from room temperature to 70 °C, we set the duty ratio to 100% and then modified it to compensate for the heat loss and maintain the SMP temperature at 70 °C. The changes in temperature and strain with time for samples i and ii are shown in Fig. 4b, c, respectively. After keeping the samples at room temperature for 100 s, the temperature was increased and maintained at 70 °C by PWM control. The strain below and above 70 °C was almost constant although the values were very different due to thermal expansion of the SMP sheet. The difference in the strain above T_{g} for the two sensors is attributed to the manufacturing errors shown in Table 3. The fluctuation of the strain at 70 °C is thought to be due to the viscoelasticity of the SMP sheet, which caused a change in its shape. In the experiments described in the “Experiments” section, we evaluated the prototypes after the temperature became almost constant.
We captured a thermogram of the heated SMP sheet shown in Fig. 5a using an infrared thermal camera (NEC Avio Infrared Technologies Co., Ltd., F30W), and the results are shown in Fig. 5b. The entire sheet was heated uniformly to approximately 70 °C. One reason for the temperature uniformity is the reduced distance between the heating wire segments. Since the tape covered the surface to attach the thermocouple, the center of the SMP sheet is shown in blue.
We created a base, a presser plate, and a cover using a 3D printer (Fig. 6). The prototype force sensor (samples i and ii) in Fig. 2 was fixed between the presser plate and the base, and was used in the experiments described in “Experiments” section. As shown in Fig. 6b), the total sensor dimensions including these parts are almost same as those for a commercial load cell (LVSA, Kyowa Electronic Instruments Co., Ltd.).
Experiments
The experimental apparatus is shown in Fig. 7. The applied force was measured at temperatures above and below T_{g}. The experiments below T_{g} were performed at room temperature. The relationship between the strain and the force applied using an indenter connected to the load cell was then evaluated. The indenter was placed in contact with the steel plate of the prototype sensor in order to prevent SMP surface deformation. The load cell and the sensor were attached to a manual stage and an automatic stage (OSMS2085, Sigma Koki Co., Ltd.), respectively. The prototype sensor was automatically displaced using the automatic stage. The strain gauge was connected to a PC through a bridge box (DB120A, Kyowa Electronic Instruments Co., Ltd.) and a strain amplifier (DPM711B, Kyowa Electronic Instruments Co., Ltd.). The load cell was also connected to the PC through a strain amplifier. The sampling frequency was 1 kHz. We resampled the obtained 1kHz signal at 100 Hz using the resample() function in MATLAB.
Experiment 1
We performed two types of experiments in order to characterize the proposed sensor. We first applied a random force to the prototype sensor in order to estimate the optimum transfer function and compare the proposed sensor with that used in our previous study [4]. The sensor was deformed as follows:
Step 1: The sensor was held motionless in the unloaded state (just before touching).
Step 2: After the unloaded state, the sensor was moved in the direction of the blue arrow in Fig. 7 and was brought into contact with the load cell to apply a deformation of 1 mm to the tip of the sensor.
Step 3: Leaving the tip deformed, the sensor was held motionless.
Step 4: The sensor was returned to the initial position.
Steps 1 through 4 were repeated. The number of repetitions was 5, which was larger than in our previous study [4]. For system identification, the input should be “persistently exciting”, i.e., it should contain many distinct frequencies [23]. Therefore, in the present study, considering the potential applications of our sensor (i.e., a wide range of inputs), we randomly set the velocity in Steps 2 and 4 from 0.5 to 5 mm/s, and the rest time in Steps 1 and 3 from 0 to 20 s. For each condition, the measurements were conducted six times. We evaluated two prototype sensors (samples 1 and 2). We prepared samples 1 and 2 by attaching samples i and ii, respectively, between the presser and base plates, as shown in Fig. 6b.
Experiment 2
We then evaluated the dynamic response of the proposed sensor to a step deformation. Similarly to experiment 1, the sensor was deformed as follows:
Step 1: The sensor was held motionless in the unloaded state (just before touching) for 10 s.
Step 2: After the unloaded state, the sensor was moved in the direction of the blue arrow in Fig. 7 (5 mm/s) and was brought into contact with the load cell in order to apply a deformation to the tip of the sensor.
Step 3: Leaving the tip deformed, the sensor was held motionless for 300 s.
We set the deformation in Step 2 to 0.25, 0.5, 0.75, or 1 mm in order to check whether the force increases with increasing deformation. For each condition, the measurements were conducted three times. We evaluated two prototype sensors (samples 1’ and 2’). After experiment 1, we removed the sensors from the presser and base plates, and reassembled samples 1’ and 2’ manually, as shown in Fig. 6b. Therefore, l in Fig. 2 was not exactly the same for samples 1 (2) and 1’ (2’).
Transfer function
Generalized Maxwell model
Here, \(\dot{\varepsilon }\) and ε are measured by the strain gauges, so \(\dot{W}\) is determined using Eq. (5) and the value of W at each time step. If W = 0 at t = 0, the value of \(\dot{W}\) determined from Eq. (5) can be substituted into Eq. (7) in order to sequentially calculate W:
In our previous study [4], after performing this procedure multiple times, the force measured using the load cell was compared with W obtained based on the strain measured using the strain gauge, and the error was calculated. Then, W was determined by combining the theoretical equations in “Basic concept of force sensor” section. In MATLAB, the leastsquares method was then used to determine the optimum values of L, M, and N in Eq. (5).
In the present study, in order to describe the viscoelastic behavior more accurately, we derived a transfer function using the generalized Maxwell model, as shown in Fig. 8 [13]. Westbrook developed a generalized Maxwell model to capture the shapememory effect using two sets of nonequilibrium branches for two fundamentally different modes of relaxation: the glassy mode and the Rouse modes [13]. Since the temperature of our sensor is fixed above and below T_{g}, we neglected the thermal expansion of the SMP. Then, using the transfer function, we calculated the force from the strain.
The derivation process is shown below. As shown in Fig. 8, the forces applied to each element \({f}_{0}\), \({f}_{i}\) (\(i=1\cdots n\)) are expressed as follows:
where A is the crosssectional area, u_{id} and u_{is} are the displacements of the dashpot and the spring, respectively, \({\eta }_{i}\) is the viscosity, E_{0} and \({E}_{i}\) are elastic moduli, L is the total length of the proposed model, and \({L}_{1}\) and \({L}_{2}\) are the length of the dashpot and the spring, respectively. The total applied force (f) and the displacement (u) of the model are expressed as follows:
Using Eqs. (8) and (11), we have
By Laplace transformation of Eqs. (9) and (12), we obtain
where
After Laplace transformation of Eq. (10), by substituting Eqs. (13) and (14), we have
Using Eq. (17), the transfer function G_{1}(s) for the generalized Maxwell model, shown in Fig. 8, with an input U(s) and an output F(s), is expressed as follows:
As shown in Eq. (18), the numbers of poles and zeros are n. In the present study, we assumed that the transfer function G(s) with an input ε(s) and an output W(s) had a structure similar to G_{1}(s) and n poles and n zeros. Namely, G(s) is expressed as follows:
Since the temperature of the proposed sensor is fixed above or below T_{g}, we assumed that \({\alpha }_{i}\)(i = 0 \(\cdots n\)) and \({\beta }_{i}\) (i = 0 \(\cdots n\)) are constants in the rubbery and glassy states. Note that by Laplace transformation of Eq. (5), G(s) with n = 1 can be obtained. We determined n using the experimental results.
Calculation of force
Using the results of experiments 1 and 2 shown in “Experiments” section, we first estimated the transfer function shown in Eq. (19) and calculated the fit ratio [24]. The fit ratio (hereinafter denoted as FIT) is defined as
where y(k) and \(\widehat{y}(k)\) are the measured output and the simulated output, respectively, at time k, \(\overline{y }\) is the average value of y, and m is the number of samples. When y(k) is identical to \(\widehat{y}(k)\) for all k = 1, …, m, FIT becomes 100. The procedure is as follows:

Method 1 (proposed in the present study): We estimated G(s) (n = 1–5) using the tfest() function in MATLAB. For the estimation, we imposed the condition that α_{i} and β_{i} are positive because the values of L, A, \({E}_{0}\), \({B}_{i}\), and \({C}_{i}\) in Eq. (18) are positive. The initial values of α_{i} and β_{i} were set to 1 or NaN. The use of NaN indicates unknown coefficients. First, regarding each experimental data, we determined α_{i} and β_{i} with maximum FIT value as the optima. Second, we calculated the average values of α_{i} and β_{i} below and above T_{g}, and determined G(s) with n = 1–5. We then calculated FIT using the compare() function to validate the determined G(s). Comparing FIT for each condition, we determined the optimal value of n.
Furthermore, using the method in our previous study [4] (namely Eqs. (5) and (7), method 2) and the results of experiment 1, we estimated W, calculated FIT, and compared the obtained results with those measured by Method 1. The procedure is as follows.

Method 2 (used in our previous study): Using Eqs. (5) and (7), we calculated L, M, and N. In MATLAB, the leastsquares method was used to determine the optimum values of L, M, and N in Eq. (5). Here, L, M, and N were changed at intervals of 10 N s^{−1}, 10 Ns, and 0.1 s, respectively. Based on these values, we also calculated FIT.
Results and discussion
Identification of transfer function (experiment 1, method 1)
Typical transitions of the force below and above T_{g} are shown in Fig. 9. As shown in Fig. 9, the two measured forces are almost identical. Moreover, by considering the viscosity of the SMP, the estimated force can reproduce the stress relaxation phenomenon. The maximum force and strain for different conditions are shown in Fig. 10. Although we applied a similar displacement below and above T_{g}, the measured force range is significantly different. Based on the above results, it was shown that the miniaturized sensor achieved the same performance as in our previous study. On the other hand, the maximum strain below and above T_{g} is similar. The measured force below T_{g} is smaller than expected (27 N, see “Design of miniaturized sensor” section for additional details) and both the maximum force and strain are different for the two prototype sensors. These differences may be attributed to errors during sensor manufacture, as shown in Table 3.
On the other hand, when the sensor was returned to the initial position, the force measured by the load cell became zero although the estimated values were not zero (arrow in Fig. 9b). One reason is that the SMP sheet could not recover to the initial shape quickly because of its viscosity, and the indenter of the load cell could not contact the prototype sensor. Above T_{g}, a fluctuation of the estimated force can be seen. The transitions of the measured strain are shown in Fig. 11. Similarly to Fig. 9, fluctuations can also be seen, and are attributed to electrical noise in the heating wire.
Using Method 1 (n = 1–5), we calculated FIT. The mean ± the standard deviation below and above T_{g} are shown in Fig. 12a, b, respectively. When n = 4 and 5, by imposing the condition that the initial values of α_{i} and β_{i} are both 1, an error occurred during the calculation in MATLAB. When n = 5, below T_{g}, there was a case when FIT was negative value, and the average value was small, as shown in Fig. 12a. The values of FIT below T_{g} are larger than those above T_{g}. One reason would be the electrical noise shown in Figs. 9 and 11. As shown in Fig. 12, except for n = 5, as the order of the transfer function became large (n increased), the FIT values also increased.
We then calculated the average values of the coefficients in Eq. (19) below and above T_{g}, and determined G(s) with n = 1–5. Using the same experimental results, we calculated the mean and the standard deviation of FIT below and above T_{g} for each n. The calculated FIT values are shown in Fig. 13; for most conditions, the largest value was obtained for n = 3. Therefore, a transfer function with n = 3 is considered to be optimal.
Bode diagrams for the transfer function (n = 3) obtained using Method 1 below and above T_{g} are shown in Fig. 14. It can be seen that the gain has a significant temperature dependence. The phases above T_{g} were larger than those below T_{g}, which can be attributed to the viscosity of the SMP. There were several differences between the gains and the phases for the two sensors, which may be attributed to the manufacturing errors shown in Table 3.
Comparison with our previous studies (experiment 1, comparison between Methods 1 and 2)
The FIT values determined using Method 2 are shown in Fig. 15. As shown in Figs. 12 and 15, the maximum values of FIT determined using Method 1 were larger than those determined using Method 2. Using the proposed transfer function model, the FIT values became larger than those in our previous studies.
Bode diagrams of the transfer function obtained using Method 2 below and above T_{g} are shown in Fig. 16. The gains in Fig. 16 are similar to those in Fig. 14. However, above T_{g,} the phases in Fig. 16 are smaller than those in Fig. 14. This difference would cause a decrease in FIT.
Step deformation response (Experiment 2)
Using Method 1 (n = 1–5) and the results of experiment 2, we calculated FIT for different conditions. The mean ± the standard deviation below and above T_{g} are shown in Fig. 17a, b, respectively. When n = 4 and 5, by imposing the condition that the initial values of α_{i} and β_{i} are both 1, an error occurred during the calculation in MATLAB. When n = 4, there were two cases when FIT was negative value, and its average value was small, as shown in Fig. 17. Similarly to experiment 1, Fig. 17 shows excellent agreement between the model estimations and the experimental data, although the deformations of the SMP sensors were not the same. As shown in Figs. 12 and 17, the values of FIT in experiment 2 were smaller than those in experiment 1. One reason would be the changes in the temperature and the sensor characteristics over time, as shown in Fig. 4b, c, because experiment 2 (more than 310 s) was longer than experiment 1 (about 104 s). In the future, we will attempt to minimize these effects, for example by using improved temperature control, to determine the optimal value of n.
We then calculated the average values of the coefficients in Eq. (19) below and above T_{g}, and determined G(s) for n = 1–5. Using the same experimental results, we calculated the mean and the standard deviation of FIT below and above T_{g} for each n. The calculated values of FIT are shown in Fig. 18, and are seen to be much lower than those in Fig. 17. The main reason is that the coefficients in Eq. (19) are different for each deformation. Therefore, for practical applications of the proposed sensor, it would be necessary to set the optimum coefficients according to the operating conditions.
Conclusion
We have developed a variablesensitivity force sensor using an SMP sheet with an embedded electrical heating wire. In the present study, we miniaturized this type of sensor while referencing the dimensions and rated capacity of a commercial load cell. The volume was decreased by 87% compared with that in our previous study. The entire sheet of the prototype sensor was heated uniformly to approximately 70 °C.
Moreover, we proposed a transfer function using a generalized Maxwell model. Using identification experimental results, we determined the numbers of poles and zeros and compared the FIT value between our previous and present studies. Models were introduced and were validated experimentally, and there was excellent agreement between the model estimations and the experimental data. A transfer function with n = 3 was found to be optimal. Using the proposed model, the FIT value became larger than in our previous studies.
Availability of data and materials
Not applicable.
References
 1.
Mukai T, Hirano S, Nakashima H et al (2011) Manipulation Using Tactile Information for a NursingCare Assistant Robot in WholeBody Contact with the Object. Trans JAPAN Soc Mech Eng Ser C 77:3794–3807. https://doi.org/10.1299/kikaic.77.3794
 2.
Murozaki Y, Sakuma S, Arai F (2017) Improvement of the Measurement Range and Temperature Characteristics of a Load Sensor Using a Quartz Crystal Resonator with All Crystal Layer Components. Sensors 17:1067. https://doi.org/10.3390/s17051067
 3.
Takashima K, Kamizono H, Takenaka M, Mukai T (2017) Force sensor utilizing stiffness change of shapememory polymer based on temperature. ROBOMECH J 4:17. https://doi.org/10.1186/s4064801700862
 4.
Takashima K, Onoda R, Mukai T (2018) Evaluation of Error and Sensitivity for Force Sensor Using ShapeMemory Polymer. In: 2018 57th Annual Conference of the Society of Instrument and Control Engineers of Japan, SICE 2018. Institute of Electrical and Electronics Engineers Inc., pp 1433–1438
 5.
Barot G, Rao IJ (2006) Constitutive modeling of the mechanics associated with crystallizable shape memory polymers. Zeitschrift fur Angew Math und Phys 57:652–681. https://doi.org/10.1007/s0003300500096
 6.
Wang ZD, Li DF, Xiong ZY, Chang RN (2009) Modeling thermomechanical behaviors of shape memory polymer. J Appl Polym Sci 113:651–656. https://doi.org/10.1002/app.29656
 7.
Tobushi H, Hara H, Yamada E, Hayashi S (1996) Thermomechanical properties in a thin film of shape memory polymer of polyurethane series. Smart Mater Struct 5:483–491. https://doi.org/10.1088/09641726/5/4/012
 8.
Takashima K, Sugitani K, Morimoto N et al (2014) Pneumatic artificial rubber muscle using shapememory polymer sheet with embedded electrical heating wire. Smart Mater Struct 23:125005. https://doi.org/10.1088/09641726/23/12/125005
 9.
Takashima K, Noritsugu T, Rossiter J et al (2012) Curved type pneumatic artificial rubber muscle using ShapeMemory Polymer. J Robot Mechatronics 24:472–479
 10.
Abavisani I, Rezaifar O, Kheyroddin A (2021) Multifunctional properties of shape memory materials in civil engineering applications: a stateoftheart review. J Build Eng. https://doi.org/10.1016/j.jobe.2021.102657
 11.
An Y, Okuzaki H (2020) Novel electroActive shape memory polymers for soft actuators. Jpn J Appl Phys 59:061002. https://doi.org/10.35848/13474065/ab8e08
 12.
Tobushi H, Hashimoto T, Hayashi S, Yamada E (1997) Thermomechanical Constitutive Modeling in Shape Memory Polymer of Polyurethane Series. J Intell Mater Syst Struct 8:711–718. https://doi.org/10.1177/1045389X9700800808
 13.
Westbrook KK, Kao PH, Castro F et al (2011) A 3D finite deformation constitutive model for amorphous shape memory polymers: A multibranch modeling approach for nonequilibrium relaxation processes. Mech Mater 43:853–869. https://doi.org/10.1016/j.mechmat.2011.09.004
 14.
Niwa Y, Ikeda T, Senba A (2014) Micromechanical model of shape memory polymer including temperaturestrain history dependence. Trans JSME (in Japanese). https://doi.org/10.1299/transjsme.2014smm0310
 15.
Leng J, Lan X, Liu Y, Du S (2011) Shapememory polymers and their composites: Stimulus methods and applications. Prog Mater Sci 56:1077–1135
 16.
Kim JH, Kang TJ, Yu WR (2010) Thermomechanical constitutive modeling of shape memory polyurethanes using a phenomenological approach. Int J Plast 26:204–218. https://doi.org/10.1016/j.ijplas.2009.06.006
 17.
Diani J, Liu Y, Gall K (2006) Finite strain 3D thermoviscoelastic constitutive model for shape memory polymers. Polym Eng Sci 46:486–492. https://doi.org/10.1002/pen.20497
 18.
Wang Y, Wang J, Peng X (2021) Refinement of a 3D finite strain viscoelastic constitutive model for thermally induced shape memory polymers. Polym Test 96:107139. https://doi.org/10.1016/j.polymertesting.2021.107139
 19.
Liu Y, Gall K, Dunn ML et al (2006) Thermomechanics of shape memory polymers: Uniaxial experiments and constitutive modeling. Int J Plast 22:279–313. https://doi.org/10.1016/j.ijplas.2005.03.004
 20.
Fraden J (2003) Handbook of modern sensors: Physics, designs, and applications, 3rd edn. Springer, Berlin
 21.
Shibata T, Otani R, Komai K, Inoue T (1991) Zairyo Rikigaku No Kiso. Baifukan, Tokyo
 22.
Kawahara T, Tokuda K, Tanaka N, Kaneko M (2006) Noncontact impedance sensing Artif Life Robot 10:35–40. https://doi.org/10.1007/s1001500503637
 23.
Ljung L (1999) System Identification: Theory for the User, 2nd Edition  Pearson. Prentice Hall
 24.
Muroi H, Adachi S (2015) Model Validation Criteria for System Identification in Time Domain. IFACPapersOnLine 48:86–91. https://doi.org/10.1016/j.ifacol.2015.12.105
Acknowledgements
Not applicable.
Funding
The present study was supported by JSPS KAKENHI Grant Numbers JP17K06265 and JP20K04401.
Author information
Affiliations
Contributions
KT conceived and led the study, and wrote this paper as the corresponding author. JK developed the sensors, carried out all experiments, and analyzed the data. NK, KT and TM participated in the research design. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Takashima, K., Kobuchi, J., Kamamichi, N. et al. Characterization of variablesensitivity force sensor using stiffness change of shapememory polymer based on temperature. Robomech J 8, 24 (2021). https://doi.org/10.1186/s40648021002108
Received:
Accepted:
Published:
Keywords
 Shapememory polymer
 Force sensor
 Glass transition temperature
 Cantilever
 Strain gauge
 System identification
 Viscoelasticity