Electrical Resistivity-Based Study of Self-Sensing Properties for Shape Memory Alloy-Actuated Artificial Muscle

Shape memory alloy (SMA) has great potential to develop light and compact artificial muscle (AM) due to its muscle-like high power-to-weight ratio, flexibility and silent operation properties. In this paper, SMA self-sensing properties are explored and modeled in depth to imitate the integrated muscle-like functions of actuating and self-sensing for SMA-AM based on the investigation of SMA electrical resistivity (ER). Firstly, an ER transformation kinetics model is proposed based on the simulation of SMA differential scanning calorimetry (DSC) curves. Then a series of thermal-electrical-mechanical experiments are carried out to verify the validity of the ER model, whereby the SMA-AM self-sensing function is well established under different stress conditions. Finally the self-sensing capability is further demonstrated by its application to a novel SMA-AM-actuated active ankle-foot orthosis (AAFO).


OPEN ACCESS
SMA-based AM not only can be used for actuating, such as in artificial hands [8,9] or bipeds [10], but also can be used for self-sensing, such as the flexures designed by Lan et al. [11,12] and the AM developed by Zhang et al. [2]. For decades, many scholars have conducted and are still conducting research on SMA self-sensing properties in order to imitate the integrated muscle-like functions of actuating and self-sensing for SMA-AM [2,[11][12][13][14]. The self-sensing function mainly refers to using SMA resistance as the feedback signal for SMA length. It has been recognized that resistance plays a significant role in SMA self-sensing utilization, while there has been little research that explores and develops the self-sensing models due to the complexity of electrical resistivity (ER) variations. SMA mainly has two phases, namely martensite (M) and austenite (A). When being heated, the SMA changes its shape from M to A. Subsequent cooling returns its shape from A to M under stress [15]. For a certain given SMA wire, the intermediate R-phase (R) will also appear between the transformation from A to M during cooling. The ER is not only dependent on the stress and temperature, but also on the volume fractions of M, A and R during transformation [16]. Furthermore, the intermediate R-phase has different ER properties comparing with M and A, which further increases the complexity of the self-sensing modeling.
Up to now, few self-sensing models have been proposed in the literature [2,11,14]. Lan et al. [11] developed a self-sensing model directly based on the SMA strain to resistance curves by polynomial fitting without considering the ER effect. Nová k et al. [14] proposed a micromechanical model for ER variations based on s detailed analysis of the factors influencing ER variations from a microscopic perspective, such as stress, temperature, martensite fraction, martensite texture and R-phase distortion. Reasonable agreement between the simulation and experimental result of resistivity responses was achieved. However, apart from stress and temperature, the other three factors as mentioned above are difficult to measure by real-time control methods, which is not conducive to feedback control applications. In addition, these three factors are closely related to stress and temperature, and the ER variations can actually be attributed to the direct and indirect effect of stress and/or temperature. Besides, though R-phase may complicate the ER variations, the ER curves of A→M transformation and A→R transformation have good similarity according to our experimental results. Furthermore, similar to the linearity between resistance and strain during A-M forward and reverse transformation, the resistance and strain also has good linearity during A-R forward and reverse transformation [16]. Thus it is probable to develop more practical and accurate self-sensing model based on ER variations from control perspective.
The motivation of this paper was to propose a self-sensing model for SMA-actuated AM to establish resistance to length (R-L) relation based on deeply exploring inherent ER properties of SMA wires. The SMA-AM presents integrated muscle-like functions of actuating and self-sensing. Finally, we propose a novel active ankle-foot orthosis (AAFO) actuated and sensed by the SMA-AM to demonstrate the self-sensing capability. Experimental results show that the SMA-AM is robust to external stress variations and accurate angle control using is achieved the self-sensing method.

Analysis of ER Properties of SMA Actuated AM
SMA wire can be used as the active contractile element for artificial muscle design due to its muscle-like properties. The designed SMA-actuated AM, as shown in Figure 1, presents multifunctional actuating, energy-storing and self-sensing properties [1,2]. To obtain an accurate self-sensing function, the resistance-length (R-L) relationship should be correctly established. However, the R-L relationship has not yet been entirely explored due to the complexity of electrical resistivity (ER) variations. For a certain given SMA wire, ER variations not only depend on the direct effects of stress and temperature, but also on stress and/or temperature induced phase transformation, martensite reorientation and R-phase distortion. The detailed relationships of the influencing factors on ER variations are shown in Figure 2. As with common metals, the ER of SMA changes linearly with temperature and stress at stable individual phases (black arrow) since SMA is a metal. However, the ER is also closely related to phase transformation induced by stress and/or temperature (red arrow), and martensite reorientation induced by applied stress (blue arrow), which is in contrast to the behavior of common metals. Besides, when the R-phase appears in SMA, the ER also changes linearly with R-phase distortion induced by temperature (green arrow). Consequently, ER variations are attributed to the direct and indirect effect of stress and/or temperature for a certain given SMA wire. In this study, equiatomic NiTi wires were investigated as the sample obtained from Dynalloy, Inc. (Tustin, CA, USA). When no transformation occurs, the temperature coefficient of resistivity (TCR) and stress coefficient of resistivity (SCR) can be used to denote the direct effect of temperature and stress on the ER. Furthermore, martensite reorientation is proportional to applied stress and R-phase distortion is inversely proportional to temperature [14]. Thus, positive SCR and negative TCR can be used to denote the effect of martensite reorientation and R-phase distortion on ER respectively, and all the coefficients gained from experiments are listed in Tables 1 and 2. Therefore, the following linear equation can be obtained to express ER variations of individual phases outside the phase transformation region: where subscript i represents A, M or R, T is temperature and  is applied stress on SMA wire. 0 0i  is ER of individual phases measured at temperature 0i T and zero stress. i  is TCR and i  is SCR and iT   is ER of individual phases at temperature T and stress  .  During phase transformation, ER variations depend mainly on the volume fractions of austenite and martensite phases, which in turn depend on the stress and/or temperature. The appearance of R-phase in SMA wire complicates the change of ER, but we find that the ER variations of A→M and A→R have good similarity. Thus we mainly focus on the study of ER variations of phase transformation without R-phase, which can be deduced reasonably to the situation with R-phase.

DSC Curve Analysis and ER Modeling
In order to study ER variations during transformation, DSC tests were performed to analyze the phase transformation of SMA wire under zero stress by employing a Netzsch DSC 200 F3 Maia ® Instrument (Burlington. MA, USA). Three identical thermal cycles were performed on the sample. For each cycle, the temperature cycle range was from −10 °C to 120 °C at a constant heating and cooling rate of dT/dt = 10 °C /min. Figure 3 shows the third cycle heat flow vs. temperature curves of the SMA wires. An endothermic and two exothermic peaks (one exothermic peak is apparent while the other is not) are present due to latent heat absorption and release. s M , p M and f M are martensite start, peak and finish temperature respectively during cooling. A s , A p and A f are austenite start, peak and finish temperature, respectively, during heating. R p is R-phase peak temperature, while R-phase start R s and finish temperature R f are not detectable by DSC. Here we temporarily ignore the impact of R-phase. Note that the temperature parameters mentioned above approximately satisfy the following relationships: In order to clarify thermal transformation behavior of SMA, we use a "Cosine curve" to simulate relationship between heat flow and temperature during transformation as shown in Figure 3 (black line). This approach relies on the notion that the gradual phase transformation may be represented by relatively symmetric DSC peaks. Here reverse transformation is taken as a case study and the expression for heat flow and temperature is given as follows: where T is temperature of SMA, h and 1 C are constants and () Furthermore, latent heat Q associated with phase transformation can be calculated by integrating Equation (3) between s A and f A . They satisfy the following expression: According to the analysis in the last section, the ER variations of SMA wires depend mainly on the volume fractions of austenite and martensite phases during transformation, which actually depend on the extent of phase transformation. In addition, the gradual change of latent heat absorption and release can be used to represent the extent of phase transformation [17], thus the following equation can be obtained: where k and 2 C are constants. Substituting Equations (3) and (4)  we can obtain the following relationship between ER and temperature during reverse transformation: where 0 and 0 are the start and end ER of reverse transformation at zero stress, which can be calculated by Equation (1). Likewise, we can obtain ER vs. temperature during forward transformation: where 0 and 0 are the start and end ER of forward transformation at zero stress, which can be calculated by Equation (1). Besides, When stress is applied to SMA during thermal cycle, SMA transformation temperature will vary linearly [18], and the ER will also vary correspondingly. Besides, the ER varies linearly with applied stress according to our experimental results. Similar results are shown in literature [19]. In order to correctly model the relationship of ER, temperature and applied stress during reverse transformation, the following expression is proposed: and A C indicates the influence of stress on reverse transformation temperatures.
A  is SCR of austenite. Likewise, we can obtain ER, temperature and applied stress relationship during forward transformation: M  is the SCR of martensite. When R-phase appears, forward transformation will be separated into two steps (A→R and A→M), and the following two expressions can be obtained based on the same analysis of A→M: where 0 and 0 are the start and end ER of A→R transformation at zero stress, which can also be calculated by Equation (1) In the next section, we will validate the efficiency and accuracy of the developed model through experiments.

Experimental Setup
To study ER characteristics during transformation, we conducted a series of thermal-electricalmechanical experiments to the as-received SMA sample, as shown in Figure 4. The sample with diameter of 0.15 mm and length of 117 mm is connected to different constant loads varying from 28 to 224 MPa and a helical bias spring with stiffness of 390 N/m. A load cell is connected to the SMA wire to measure the contraction force. A linear variable differential transformer (LVDT) with resolution of 0.01 mm measures the SMA wire length change. Besides, a thermocouple (Omega type K) is positioned at the center of the SMA wire to measure the temperature. The wire is sheltered from air draughts. Heating current generated from current amplifier is used to drive the SMA wire and the resistivity is computed according to the following equation [20]: where U is the voltage and I is the heating current across the SMA wire. l is the wire length. 0 V is the wire volume which can be considered as a constant. A computer equipped with a data acquisition card (NI USB-6211) is used to store and process all the measured data.

Heat Transfer Analysis of SMA Wire
The thermocouple is capable of detecting steady state temperature or slow response temperature, but it is hard to detect fast response temperatures, and it is also hard to detect the temperature changes induced by latent heat absorption/release during fast phase transformation. In addition, with a view towards real-time feedback control applications, it is more practical to control heating/cooling time rather than temperature of SMA wires. Therefore, heat transfer analysis based on lumped parameter method (LPM) and energy conservation principle is performed to investigate the temperature response of SMA wires. Then the temperature vs. time and ER vs. time relationship can be established.
The heat transfer models outside phase transformation and during phase transformation are written as Equations (13a) and (13b) respectively: where m is the mass and surf A is the surface area of SMA wire, V C is heat capacity, R is the average resistance of SMA wire. amb T is the reference temperature, taken to be identical to the ambient temperature of the environment, I is the electric current, and h is the convection coefficient. Note that radiative effects are lumped into the convective term. i Q is the latent heat associated with the phase transformation. The subscript i refers to M→A during heating and A→M during cooling. Latent heat release of A→R transformation is ignored because no obvious heat flow occurs, as shown in Figure 3. In addition,  refers to martensitic volume fraction, which may be taken to be solely dependent on temperature for purely thermal transformation [21], and the specific functional dependence can follow from a simulation of DSC curves of the heat flow vs. temperature for the SMA wires. Considering Equation (3) The thermal response is obtained by numerically integrating the above Equations (13a)-(14b), and the temperature vs. time curve during heating (I = 350 mA) and cooling (I = 50 mA) is shown in Figure 5b. Note that outside the transformation region, the temperature changes in exponential form (segments AB, CD and segments DE, FG), while during transformation the latent heat absorption and release cause the sudden change of the slope, as is shown in segments BC and EF, which denote M→A and A→M transformation, respectively. The simulation results in Figure 5b are consistent with experimental results in literature [21]. For simplicity, temperature vs. time relationship can be approximately considered linear during transformation according to the simulation results, and the following temperature vs. time functions are obtained: The corresponding experimental results of ER (red line) and strain (blue line) response of the SMA wire are shown in Figure 5a. During heating, the ER response is non-monotonic with a rapid increase, decrease then increase again. In segment BC, the SMA wire undergoes reverse transformation, where the ER decreases because the ER of austenite is smaller than that of martensite. In segments AB and CD, no transformation occurs, where the ER increase is due to heating. In addition, the strain changes slightly in segments AB and CD but decreases rapidly in segment BC due to reverse transformation. During cooling, the opposite characteristics of ER and strain response occur. In segment EF, the ER increases due to forward transformation while in segments DE and FG, the ER decrease is due to cooling. Churchill et al. [22] showed similar results to those presented in Figure 5a through an external heating method.  Figure 6a,c show the experimental results of strain and corresponding ER curves during heating at various constant stresses. With the increase of the stress, the start (denoted by 12 AA ) and end (denoted by 12 BB ) temperature of the reverse transformation increase correspondingly. In addition, the whole ER value also increases in proportion to the stress. Besides, the simulation results of ER during heating can be obtained by combining Equations (1) and (9) with Equation (15), as shown in Figure 7a. The simulations agree reasonably well with the experimental results. Figure 6b,d shows the strain and corresponding ER curves during cooling at various constant applied stresses. Comparing with Figure 6a, the strain curves shown in Figure 6b are quite different as stress varies, which is caused by R-phase [17]. In contrast to the A→M transformation which yields a large shape change (5%), the A→R transformation yields only a small shape change (less than 1%).  During cooling, As is shown in Figure 6d, the ER increase mainly denotes A→R transformation (from to ) and A→M transformation (from E 1 E 2 to F 1 F 2 ) when the stress is lower than 56 MPa and higher than 196 MPa. This is because all the transformation temperatures increase linearly with increasing applied stress, but R s exhibits a stress dependence substantially smaller than that of M s or A f [16], thus R-phase fraction decreases gradually with the increase of the stress. For stresses greater than a certain value, M s exceeds R s and only one step (A→M) appears during the forward transformation.

Experimental Results at Various Constant Stresses
A revised critical stress-temperature plot including R-phase is shown schematically in Figure 8 to illustrate this situation. If applied stress is higher than σ f , no R-phase occurs during cooling (arrow a). If stress is lower than s  , two separate steps (A→R→M) appear successively (arrow c). If stress is within the above range, two mixed steps (A→R+M) appear (arrow b). Here, three constant applied stresses (28 MPa, 168 MPa and 224 MPa) are implemented as case studies.
Combining Equations (1), (10), (11) with Equations (15), we derive numerical simulation of ER variations. Figure 7b shows the comparison between simulation and experimental results during cooling. When A→R and A→M occur, the simulation results are in good agreement with the experimental results. When A→R+M occurs, the model reveals a relatively large error because the ER variations represent the mixture effect of R-phase and martensite.

Self-Sensing Modeling of SMA Actuated AM
To achieve self-sensing function of the SMA-AM, we should establish resistance-length (R-L) relationships of SMA wires. The resistivity-strain ( − ) relationship can be established firstly. According to the previous discussion, ER variations of SMA wires depend mainly on the extent of transformation during phase transformation, which can be described by martensite fraction [15]. Besides, the whole resistivity value also varies with applied stresses according to our finding. Thus the ER and martensitic fraction satisfy the following relationships: where the superscript 0 denotes the ER at free stress. i  is the stress coefficient of ER, which can be considered as a constant. In addition, the following constitutive model is used in this study [18]: where E is Young's modulus,  is thermal coefficient of expansion, and  is the transformation tensor of SMA. Besides, L E    , and L  as a constant is the maximum residual strain. Thus, we can obtain the following − expression by substituting Equation (16) into Equation (17): residual strain of R-phase transformation, which is generally far smaller than L  .
The simulation and experimental results of − curves during M→A and A→M are shown in Figure 9. During M→A, the − curves have good linearity. The whole resistivity value changes linearly with the stress and the simulation results agree well with the experimental results. However, during A→ M, a large error occurs between simulation and experimental results at lower applied stresses. This is due to the appearance of R-phase during cooling. With the increase of the stress, R-phase fraction decreases and the linearity becomes better. The experimental resistance-length (R-L) curves are shown in Figure 10a. When the applied stress is large, R-phase is inhibited and the linearity of R-L curve is quite good and the hysteresis gap of heating and cooling path is also small. However, when the applied stress is small, the resistance-length curve has two different steps with very different slopes during cooling caused by R-phase, which is not convenient for self-sensing application. Here we are only concerned with the case without R-phase. Assuming the initial length of SMA wire is 0 l , then: Combining Equations (12), (18) and (19), we can gain the following R-L relationship: The values h R and c R are the resistance on the heating path and cooling path respectively. As the heating and cooling paths are very close to each other, we average the two paths as follows: Combining Equations (19) and (20), we can gain simulation results of R-L relationships under different applied stresses, which agree well with experimental results as shown in Figure 10b.  Figure 4. The R-L curves with and without bias spring are shown in Figure 11. The differences between the two curves are minimal during heating path (H arrow) and they have the same slope. However, during cooling process (C arrow), there is a big difference between the two curves, which is caused by SMA R-phase. Therefore, as long as the R-phase is inhibited, the two curve can be coincide well with each other. Thus, the self-sensing model (Equations (20) and (21)) can predict the R-L relationships both under contant and variable stresses, and the self-sensing can be achieved for the SMA-AM.
It should be noted that, the above self-sensing model is developed based on natural cooling. Different cooling methods (such as natural cooling and forced cooling) have different heat transfer rate, as shown in Subsection 5.1, which is represented by convection coefficient h c . However, the self-sensing is not affected by the cooling method. As it is developed from ER variations and the ER depends mainly on martensite fraction during transformation, which is independent of heat transfer rate [17].  Figure 11. Experimental results of R-L relationships under constant (without spring) and variable (with spring) loads.

Application to Active Ankle-Foot Orthosis (AAFO)
To further demonstrate the self-sensing capability of the SMA-AM, we build an Active Ankle-Foot Orthosis (AAFO) actuated and sensed by the SMA-AM and with a bias spring with stiffness of 1,144 N/m as shown in Figure 12. The SMA-AM can provide a maximum contractile force of 230 N. The maximum torque of the AAFO can be determined as 10.9 Nm. It should be noted that different patients require different AAFO assist torque, which can be achieved by adding or reducing the amount of parallel SMA wires. The bias spring serves two purposes: to provide enough restoring force for the AM and to simulate variable load condition. The orthosis combines two light-weight thermoplastic shells with articulated joints, which can provide supplemental torque assistance for weak anterior muscles of the lower leg during dorsiflexion. Besides, stable and continuous airflow is applied to cool the AMs through the cooling vessel provided by a mini pump. The current heating and data processing methods are similar to the previous experiments. Before experiments, an initial pre-transformation stress of σ = 112 MPa, is applied on the SMA-AM to inhibit the appearance of R-phase.   [11]. The control diagram is shown in Figure 13.
where l is the SMA-AM length,  is the angular position and 2 a is the moment arm of the AAFO.
The subscript 0 denotes the initial length and angle value. The resistance of the SMA wire is measured as the feedback signal and then converted to the ankle angle combining the self-sensing model (Equations (20) and (21)) with Equation (22). Furthermore, the actual angle is also used for the feedback control as a comparison. Figure 14a shows the results of sinusoidal tracking response experiments using the self-sensor and encoder respectively at 0.2 Hz. Note that, the "Reference" angle ref  is the desired sinusoidal curve produced by host PC. The "Encoder" angle is the actual curve sensed by an encoder, and the "Self-sensor" angle is the actual curve sensed by the self-sensing model. The actual angle with the two feedback methods can accurately follow the reference angle, which demonstrates the validity of the self-sensing method. When the sinusoidal frequency increases to 0.5 Hz as shown in Figure 14b, a small phase lag appears between the actual angle and the reference angle. But the two actual angle curves coincide well with each other, because both resistance and length change of SMA wires are closely related to SMA transformation, which further demonstrates the validity of the self-sensing model.

Conclusions
This paper proposes a self-sensing model to achieve the integrated muscle-like functions of actuating and self-sensing for the designed SMA-AM based on a deep investigation of SMA electrical resistivity (ER). An ER transformation kinetics model is firstly established based on SMA DSC tests. Then a series of thermal-electrical-mechanical experiments are carried out to verify the validity of the ER model, whereby the SMA-AM self-sensing function is well established under different stress conditions. A novel AAFO actuated and sensed by the SMA-AM is built to further demonstrate the self-sensing capability of the SMA-AM. It is shown that the AAFO can accurately follow the desired angle value within SMA bandwidth, which is also robust to external stress variations.