Hysteresis Curve Fitting Optimization of Magnetic Controlled Shape Memory Alloy Actuator

As a new actuating material, magnetic controlled shape memory alloys (MSMAs) have excellent characteristics such as a large output strain, fast response, and high energy density. These excellent characteristics are very attractive for precision positioning systems. However, the availability of MSMAs in practical precision positioning is poor, caused by weak repeatability under a certain stimulus. This problem results from the error of a large magnetic hysteresis in an external magnetic field. A suitable hysteresis modelling method can reduce the error and improve the accuracy of the MSMA actuator. After analyzing the original hysteresis modelling methods, three kinds of hysteresis modelling methods are proposed: least squares method, back propagation (BP) artificial neural network, and BP artificial neural network based on genetic algorithms. Comparing the accuracy and convergence rate of three kinds of hysteresis modelling methods, the results show that the convergence rate of least squares method is the fastest, and the convergence accuracy of BP artificial neural networks based on genetic algorithms is the highest.


Introduction
During the end of last century, Ullakko et al. [1] found and proposed a magnetic controlled shape memory alloy (MSMA).Then, Murray et al. [2,3] conducted research into the mechanism and characteristics of MSMA.When a sample of MSMA is placed in a suitable external magnetic field, it can produce a deformation as high as 10%.Moreover, compared with traditional shape memory alloys, MSMA has characteristics of larger strain and better energy efficiency.These excellent characteristics are very attractive for precision positioning systems [4].Since then, some scholars and organizations have designed many kinds of MSMA actuators.However, there are still relatively few applications for MSMA actuators in precision positioning.The deficiency in practical precision positioning is caused by the weak repeatability under a certain stimulus.This is the result of the error caused by a large magnetic hysteresis [5,6] in an external magnetic field.To reveal the hysteresis characteristic of MSMA [7], related experiments with NiMnGa (which is the typical MSMA material) are carried out.
According to the research and analysis of Spanish scholars, suitable hysteresis modelling methods can reduce errors and improve the accuracy of MSMA actuators [8].Thus, it is very important to reference and establish hysteresis modelling methods.
(2) Jiles-Atherton model: The advantage of the Jiles-Atherton model [12,13] is its excellent convergence, while the disadvantage is that it is difficult to achieve fast modelling due to numerous undetermined parameters.
(3) Maxwell model: Maxwell model [14,15] can improve the convergence rate by reducing the number of subsystems, but the accuracy of the model is reduced accordingly.
In addition, the proportional-integral-derivative (PID) model and other models recently became research hotspots [16][17][18].To date, none of the modelling methods developed is perfect.By analyzing the original hysteresis modelling methods, three kinds of hysteresis modelling methods are proposed for fitting the hysteresis curve: least squares method [19], back propagation (BP) artificial neural network and BP artificial neural network based on genetic algorithm.Among these methods, there are special improvements.First, BP artificial neural network that adopts the Levenberg-Marquard algorithm is a novel method and has some advantages: (a) because the Levenberg-Marquard algorithm is a type of least squares method, it can it can be directly used for fitting the hysteresis curve; (b) mean square error of the Levenberg-Marquard algorithm [20] can be used to test the rate of convergence and error size; therefore, it can be used to test the reliability and accuracy of the BP artificial neural network.Moreover, the basic principle of BP artificial neural network based on a genetic algorithm [21,22] is as follows: through the selection, crossover and mutation operation of genetic algorithms, optimize the weights and thresholds of the neural network, and then return the optimal weights and the optimal thresholds to the BP artificial neural network.By comparing the precision of fitting the hysteresis curve, the best hysteresis modelling method will be selected.

Sample and Device of the Experiment
The sample of MSMA in the experiment is an NiMnGa alloy purchased from the Goodfellow company (Huntingdon, UK).The overall dimensions of the sample are 2 × 3 × 15 mm 3 , and its appearance is shown in Figure 1.Moreover, the atomic percentages of the Ni element, Mn element, and Ga element are 50%, 28%, and 22%, respectively.
Actuators 2016, 5, 25 2 of 13 (1) Preisach model: Preisach model [9][10][11] can fit a complex hysteresis curve.However, for some materials, the fitting error is large.(2) Jiles-Atherton model: The advantage of the Jiles-Atherton model [12,13] is its excellent convergence, while the disadvantage is that it is difficult to achieve fast modelling due to numerous undetermined parameters.(3) Maxwell model: Maxwell model [14,15] can improve the convergence rate by reducing the number of subsystems, but the accuracy of the model is reduced accordingly.
In addition, the proportional-integral-derivative (PID) model and other models recently became research hotspots [16][17][18].To date, none of the modelling methods developed is perfect.By analyzing the original hysteresis modelling methods, three kinds of hysteresis modelling methods are proposed for fitting the hysteresis curve: least squares method [19], back propagation (BP) artificial neural network and BP artificial neural network based on genetic algorithm.Among these methods, there are special improvements.First, BP artificial neural network that adopts the Levenberg-Marquard algorithm is a novel method and has some advantages: (a) because the Levenberg-Marquard algorithm is a type of least squares method, it can it can be directly used for fitting the hysteresis curve; (b) mean square error of the Levenberg-Marquard algorithm [20] can be used to test the rate of convergence and error size; therefore, it can be used to test the reliability and accuracy of the BP artificial neural network.Moreover, the basic principle of BP artificial neural network based on a genetic algorithm [21,22] is as follows: through the selection, crossover and mutation operation of genetic algorithms, optimize the weights and thresholds of the neural network, and then return the optimal weights and the optimal thresholds to the BP artificial neural network.By comparing the precision of fitting the hysteresis curve, the best hysteresis modelling method will be selected.

Sample and Device of the Experiment
The sample of MSMA in the experiment is an NiMnGa alloy purchased from the Goodfellow company (Huntingdon, UK).The overall dimensions of the sample are 2 × 3 × 15 mm , and its appearance is shown in Figure 1.Moreover, the atomic percentages of the Ni element, Mn element, and Ga element are 50%, 28%, and 22%, respectively.This experiment adopts a magnetic force coupling thermal loading measuring system (shown in Figure 2) to study multifield coupling properties of force and magnetic field of the MSMA under variable temperature magnetics.During the experiment, the load is provided by a hydraulic system and is measured by a force sensor.Moreover, the load ranges from 0 to 2.5 MPa.The external control magnetic field is provided by an electric current: the magnetic field range is 0-1.5 T, and the current is proportional to the magnetic field.The operating temperature is regulated by a circulating oil bath temperature loading system, and the working temperature range is 16-35 °C.When the magnetic This experiment adopts a magnetic force coupling thermal loading measuring system (shown in Figure 2) to study multifield coupling properties of force and magnetic field of the MSMA under variable temperature magnetics.During the experiment, the load is provided by a hydraulic system and is measured by a force sensor.Moreover, the load ranges from 0 to 2.5 MPa.The external control magnetic field is provided by an electric current: the magnetic field range is 0-1.5 T, and the current is proportional to the magnetic field.The operating temperature is regulated by a circulating oil bath temperature loading system, and the working temperature range is 16-35 • C. When the magnetic field is perpendicular or parallel to the preload stress, axial and transverse deformation of the MSMA sample will be produced.These will be measured with a resistance strain gauge.The basic characteristics of the resistance strain gauge are as follows: (1) The strain gauge is a series of BX (BX strain gauge refers to phenolic foil type strain gauge) whose properties include the following: the entire structure is sealed, stable performance, good flexibility, and applicability to the general accuracy of the sensor.The strain limit is 1.5%, and usage temperature range is −30 The strain gauge sampling frequency is 4 Hz.
properties include the following: the entire structure is sealed, stable performance, good flexibility, and applicability to the general accuracy of the sensor.The strain limit is 1.5%, and usage temperature range is −30 °C-+80 °C.
(2) The strain gauge sampling frequency is 4 Hz.
Experimental steps are as follows: (1) Turn off the power, place the sample in the cup, and fix the cup on the hydraulic loading device intermediate.(2) The hydraulic drive loading device is directly placed in the working range.At the same time, the magnetic field is set from 0 to 1.5 T. (3) Turn on the power, and set up the pre load of the sample.(4) Turn on the power supply for the temperature loading system, set the temperature parameters, and check whether the outer circulation system makes good contact.When everything is acceptable, turn on the heating power supply and the oil pump power supply, and regulate the flow rate of silicone oil circulation to prevent the silicone oil spilling.

Experimental Results
In order to achieve good experimental results, the experimental conditions must meet the following requirements.First, magnetic field intensity H is parallel to the external load.The operating temperature is 16 °C, and the size of external load is 3.9 kg.In the experimental conditions, the key original experimental data of magnetic field intensity H and magnetic induction intensity B are shown in Table 1.Experimental steps are as follows: (1) Turn off the power, place the sample in the cup, and fix the cup on the hydraulic loading device intermediate.( 2) The hydraulic drive loading device is directly placed in the working range.At the same time, the magnetic field is set from 0 to 1.5 T. (3) Turn on the power, and set up the pre load of the sample.(4) Turn on the power supply for the temperature loading system, set the temperature parameters, and check whether the outer circulation system makes good contact.When everything is acceptable, turn on the heating power supply and the oil pump power supply, and regulate the flow rate of silicone oil circulation to prevent the silicone oil spilling.( 5) Start the test.The corresponding deformation of MSMA is measured.(6) After unloading, view and save the experimental data.

Experimental Results
In order to achieve good experimental results, the experimental conditions must meet the following requirements.First, magnetic field intensity H is parallel to the external load.The operating temperature is 16 • C, and the size of external load is 3.9 kg.In the experimental conditions, the key original experimental data of magnetic field intensity H and magnetic induction intensity B are shown in Table 1.

Least Squares Method
As we know, MATLAB software (MATLAB 2013, MathWorks Company, Natick, MA, USA) is the most common software numerical calculation and graphics processing.In this chapter, we apply the least squares method with MATLAB software to fit the original experimental data.
Compared with some other least squares methods, we adopt the quadratic polynomial method for the following reasons.First, the quadratic polynomial method is one of the most common least squares methods.Second, the calculation time for the quadratic polynomial method is the shortest of the least squares methods.In this paper, we assume a fitting equation as follows: The programming procedure to fit the original experimental data is as follows: clear a = load("data.txt");x = a(:, 1); y = a(:, 2); P = polyfit(x, y, 4) Application results of the programming procedure are shown as follows: Through a simulation analysis using MATLAB software, the fitting curve and the original experimental data curve are shown in Figure 3.
Through the original experimental data curve in Figure 3, magnetic induction intensity B is maintained at approximately −3.4 T when magnetic field intensity H ranges from −0.25 T to −0.54 T. Shown as curve 1 in Figure 3, magnetic induction intensity B gradually increases from −3.4 T to 6.8 T when magnetic field intensity H ranges from −0.25 T to 0.4 T.However, shown as curve 2 in Figure 3, magnetic induction intensity B gradually decreases from 6.8 T to −3.4 T when magnetic field intensity H ranges from 0.4 T to −2.5 T. Obviously, curve 1 and curve 2 do not coincide.Therefore, as shown in Figure 3, the relation curve between magnetic field intensity H and magnetic induction intensity B is the hysteresis curve of the MSMA.With regard to Equation ( 2), because the number of undetermined parameters during the hysteresis curve fitting is lower, the fitting speed of the least squares method is faster.A wealth of information can be found in the statistics figures.For example, with regard to Figures 3 and 4, the fitting curve using the least squares method can show the trend of the hysteresis curve as 0.33% of the largest fitting error rate.More importantly, the largest fitting error rate of 0.33% is quite accurate and can improve the accuracy of the MSMA actuator.In order to further improve the fitting accuracy of the fitting curve, this paper adopts a BP artificial neural network to carry out the fitting.

BP Artificial Neural Network
This paper adopts a three-layer BP artificial neural network.The network is divided into an input layer, hidden layer, and output layer.The input and output layers are, respectively, magnetic field intensity H and magnetic induction intensity B. The number of input layer nodes and output layer nodes are both 1.Like the least squares method, we apply a BP artificial neural network with MATLAB software to fit the original experimental data.We apply equation calculations combined with programming algorithms to establish a magnetic hysteresis model that is faster and more accurate.With regard to Equation ( 2), because the number of undetermined parameters during the hysteresis curve fitting is lower, the fitting speed of the least squares method is faster.A wealth of information can be found in the statistics figures.For example, with regard to Figures 3 and 4, the fitting curve using the least squares method can show the trend of the hysteresis curve as 0.33% of the largest fitting error rate.With regard to Equation ( 2), because the number of undetermined parameters during the hysteresis curve fitting is lower, the fitting speed of the least squares method is faster.A wealth of information can be found in the statistics figures.For example, with regard to Figures 3 and 4, the fitting curve using the least squares method can show the trend of the hysteresis curve as 0.33% of the largest fitting error rate.More importantly, the largest fitting error rate of 0.33% is quite accurate and can improve the accuracy of the MSMA actuator.In order to further improve the fitting accuracy of the fitting curve, this paper adopts a BP artificial neural network to carry out the fitting.

BP Artificial Neural Network
This paper adopts a three-layer BP artificial neural network.The network is divided into an input layer, hidden layer, and output layer.The input and output layers are, respectively, magnetic field intensity H and magnetic induction intensity B. The number of input layer nodes and output layer nodes are both 1.Like the least squares method, we apply a BP artificial neural network with MATLAB software to fit the original experimental data.We apply equation calculations combined with programming algorithms to establish a magnetic hysteresis model that is faster and more accurate.More importantly, the largest fitting error rate of 0.33% is quite accurate and can improve the accuracy of the MSMA actuator.In order to further improve the fitting accuracy of the fitting curve, this paper adopts a BP artificial neural network to carry out the fitting.

BP Artificial Neural Network
This paper adopts a three-layer BP artificial neural network.The network is divided into an input layer, hidden layer, and output layer.The input and output layers are, respectively, magnetic field intensity H and magnetic induction intensity B. The number of input layer nodes and output layer nodes are both 1.Like the least squares method, we apply a BP artificial neural network with MATLAB software to fit the original experimental data.We apply equation calculations combined with programming algorithms to establish a magnetic hysteresis model that is faster and more accurate.
Through various analyses, the number of hidden layer nodes is determined as In Equation ( 3), a and b are, respectively, the number of input layer nodes and output layer nodes.In addition, c ranges from 1 to 10.When c is 3, the number of hidden layer nodes is 5. Therefore, the structure of the artificial neural network is 1-5-1, as shown in Figure 5. ( In Equation ( 3), a and b are, respectively, the number of input layer nodes and output layer nodes.In addition, c ranges from 1 to 10.When c is 3, the number of hidden layer nodes is 5. Therefore, the structure of the artificial neural network is 1-5-1, as shown in Figure 5.The activation function of neurons g(x) in the hidden layer is a bipolar compressed function: Correspondingly, the activation function of neurons h(x) in the output layer is a symmetric saturated linear satlins function: In this paper, the original BP artificial neural network that adopts the Levenberg-Marquard algorithm has some advantages: (a) because the Levenberg-Marquard algorithm is a type of least squares method, it can be directly used for fitting the hysteresis curve; and the (b) mean square error of the Levenberg-Marquard algorithm can be used to test the rate of convergence and error size; therefore, it can be used to test the reliability and accuracy of the BP artificial neural network.
The equation for mean square error in the Levenberg-Marquard algorithm is In Equation ( 6), is original experimental data, and is expected data.If mean square error is important, equally so is weight adjustment ∆ : where is the unit matrix, is the user-defined learning rate, ( ) is the error vector, and ( ) is the jacobian reciprocal matrix of the error to weights, as shown in Equation ( 8): Moreover, the equation of threshold adjustment ∆ is shown in Equation ( 9): The activation function of neurons g(x) in the hidden layer is a bipolar compressed function: Correspondingly, the activation function of neurons h(x) in the output layer is a symmetric saturated linear satlins function: In this paper, the original BP artificial neural network that adopts the Levenberg-Marquard algorithm has some advantages: (a) because the Levenberg-Marquard algorithm is a type of least squares method, it can be directly used for fitting the hysteresis curve; and the (b) mean square error of the Levenberg-Marquard algorithm can be used to test the rate of convergence and error size; therefore, it can be used to test the reliability and accuracy of the BP artificial neural network.
The equation for mean square error E mse in the Levenberg-Marquard algorithm is In Equation ( 6), orig i is original experimental data, and exp i is expected data.If mean square error E mse is important, equally so is weight adjustment ∆w: where I is the unit matrix, µ 1 is the user-defined learning rate, e (w) is the error vector, and J (W) is the jacobian reciprocal matrix of the error to weights, as shown in Equation (8): Moreover, the equation of threshold adjustment ∆θ is shown in Equation ( 9): where µ 2 is the user-defined learning rate, x is the input vector, and f θ T x is the output of the neurons.Through an analysis of the activation function, weight adjustment, and threshold adjustment, this paper uses MATLAB software to realize a hysteresis curve fitting based on a BP artificial neural network.
The fitting curve and the original experimental data curve are shown in Figure 6.
Actuators 2016, 5, 25 7 of 13 where is the user-defined learning rate, is the input vector, and ( ) is the output of the neurons.
Through an analysis of the activation function, weight adjustment, and threshold adjustment, this paper uses MATLAB software to realize a hysteresis curve fitting based on a BP artificial neural network.
The fitting curve and the original experimental data curve are shown in Figure 6.As shown in Figure 7, the fitting error rate between the fitting curve and the original experimental data curve has large fluctuations.It was found that the largest fitting error rate for the original BP artificial neural network (0.07%) is much smaller than the largest fitting error rate through the least squares method (0.33%).To further analyse the fitting effect using an original BP artificial neural network, the mean square error is introduced.As shown in Figure 8, the mean square error sharply decreases during 0-10 iterations.After 10 iterations, the mean square error slows until it stops changing.Finally, the mean square error reaches its minimum of 9.9413 × 10 at the 100th iteration.Therefore, As shown in Figure 7, the fitting error rate between the fitting curve and the original experimental data curve has large fluctuations.It was found that the largest fitting error rate for the original BP artificial neural network (0.07%) is much smaller than the largest fitting error rate through the least squares method (0.33%).
Actuators 2016, 5, 25 7 of 13 where is the user-defined learning rate, is the input vector, and ( ) is the output of the neurons.
Through an analysis of the activation function, weight adjustment, and threshold adjustment, this paper uses MATLAB software to realize a hysteresis curve fitting based on a BP artificial neural network.
The fitting curve and the original experimental data curve are shown in Figure 6.As shown in Figure 7, the fitting error rate between the fitting curve and the original experimental data curve has large fluctuations.It was found that the largest fitting error rate for the original BP artificial neural network (0.07%) is much smaller than the largest fitting error rate through the least squares method (0.33%).To further analyse the fitting effect using an original BP artificial neural network, the mean square error is introduced.As shown in Figure 8, the mean square error sharply decreases during 0-10 iterations.After 10 iterations, the mean square error slows until it stops changing.Finally, the mean square error reaches its minimum of 9.9413 × 10 at the 100th iteration.Therefore, To further analyse the fitting effect using an original BP artificial neural network, the mean square error is introduced.As shown in Figure 8, the mean square error sharply decreases during 0-10 iterations.After 10 iterations, the mean square error slows until it stops changing.Finally, the mean square error reaches its minimum of 9.9413 × 10 −5 at the 100th iteration.Therefore, compared with the least squares method, the BP artificial neural network based on a genetic algorithm can decrease the fitting error and improve the accuracy of the MSMA actuator.In addition, its convergence rate is relatively slow.

BP Artificial Neural Network Based on Genetic Algorithm
Based on an analysis of Figures 7 and 8, the convergence at the 100th iteration is slow, and 0.07% of the largest fitting-error rate is not small.
In order to accelerate the convergence rate of the original BP artificial neural network, this paper uses a genetic algorithm to improve BP artificial neural network.As shown in Figure 9, the weight and threshold values of the BP artificial neural network are optimized through selection, crossover, and mutation of the genetic algorithm.

BP Artificial Neural Network Based on Genetic Algorithm
Based on an analysis of Figures 7 and 8, the convergence at the 100th iteration is slow, and 0.07% of the largest fitting-error rate is not small.
In order to accelerate the convergence rate of the original BP artificial neural network, this paper uses a genetic algorithm to improve BP artificial neural network.As shown in Figure 9, the weight and threshold values of the BP artificial neural network are optimized through selection, crossover, and mutation of the genetic algorithm.
In this paper, the initial parameters of genetic algorithm are as follows: the population size is 20, crossover probability is 0.1, and mutation probability is 0.0001.
Through selection, crossover and mutation of genetic algorithm, the optimized weight and threshold value are as follows: The following require special attention: W 1 in Equation ( 10) is the weight between the input layer and the hidden layer, W 2 in Equation ( 11) is the weight between the hidden layer and the output layer, θ 1 in Equation ( 11) is the threshold value between the input layer and the hidden layer, and θ 2 in Equation ( 12) is the threshold value between the hidden layer and the output layer.
Based on an analysis of Figures 7 and 8, the convergence at the 100th iteration is slow, and 0.07% of the largest fitting-error rate is not small.
In order to accelerate the convergence rate of the original BP artificial neural network, this paper uses a genetic algorithm to improve BP artificial neural network.As shown in Figure 9, the weight and threshold values of the BP artificial neural network are optimized through selection, crossover, and mutation of the genetic algorithm.In this paper, the initial parameters of genetic algorithm are as follows: the population size is 20, crossover probability is 0.1, and mutation probability is 0.0001.As shown in Figure 10, it is found that 0.03% of the largest fitting error rate through an improved BP artificial neural network is much smaller than 0.33% of the largest fitting error rate through an original BP artificial neural network.The following require special attention: in Equation ( 10) is the weight between the input layer and the hidden layer, in Equation ( 11) is the weight between the hidden layer and the output layer, in Equation (11) is the threshold value between the input layer and the hidden layer, and in Equation ( 12) is the threshold value between the hidden layer and the output layer.As shown in Figure 10, it is found that 0.03% of the largest fitting error rate through an improved BP artificial neural network is much smaller than 0.33% of the largest fitting error rate through an original BP artificial neural network.To further show the advantages of the BP artificial neural network based on a genetic algorithm, as shown in Figure 11, the mean square error of the BP artificial neural network based on a genetic algorithm reaches its minimum of 4.9804 × 10 −5 at the 16th iteration.Compared with the mean-square-error curve of the original BP artificial neural network in Figure 8, a BP artificial neural network based on a genetic algorithm not only accelerates the convergence rate but also decreases the fitting error.As shown in Figure 10, it is found that 0.03% of the largest fitting error rate through an improved BP artificial neural network is much smaller than 0.33% of the largest fitting error rate through an original BP artificial neural network.Therefore, it can improve the accuracy of the MSMA actuator and then enhance the utilization rate of the MSMA actuator in precision positioning.

Validation
To validate the reliability and accuracy of the BP neural network based on a genetic algorithm, this paper makes use of additional original experimental data in experimental conditions.The experimental conditions are as follows: magnetic field intensity H is parallel to the external load; operating temperature is 16 • C; the size of the external load is 4.9 kg; and the range of the magnetic field intensity H is 0-1 T.
To test the fitting effect, through a simulation analysis using MATLAB software, the fitting curve and the original experimental data curve are shown in Figure 12.As preliminarily seen in Figure 12, the fitting accuracy of the BP neural network based on a genetic algorithm is excellent.To further show the advantages of the BP artificial neural network based on a genetic algorithm, as shown in Figure 11, the mean square error of the BP artificial neural network based on a genetic algorithm reaches its minimum of 4.9804 × 10 at the 16th iteration.Compared with the meansquare-error curve of the original BP artificial neural network in Figure 8, a BP artificial neural network based on a genetic algorithm not only accelerates the convergence rate but also decreases the fitting error.
Therefore, it can improve the accuracy of the MSMA actuator and then enhance the utilization rate of the MSMA actuator in precision positioning.

Validation
To validate the reliability and accuracy of the BP neural network based on a genetic algorithm, this paper makes use of additional original experimental data in experimental conditions.The experimental conditions are as follows: magnetic field intensity H is parallel to the external load; operating temperature is 16 °C; the size of the external load is 4.9 kg; and the range of the magnetic field intensity H is 0-1 T.
To test the fitting effect, through a simulation analysis using MATLAB software, the fitting curve and the original experimental data curve are shown in Figure 12.As preliminarily seen in Figure 12, the fitting accuracy of the BP neural network based on a genetic algorithm is excellent.As shown in Figure 13, the largest fitting error rate of BP artificial neural network is 0.06%, and the largest fitting error rate of the BP artificial neural network based on a genetic algorithm is 0.025%.Thus, a BP artificial neural network based on a genetic algorithm can significantly reduce the fitting error.As shown in Figure 13, the largest fitting error rate of BP artificial neural network is 0.06%, and the largest fitting error rate of the BP artificial neural network based on a genetic algorithm is 0.025%.Thus, a BP artificial neural network based on a genetic algorithm can significantly reduce the fitting error.To further show the advantage of the BP artificial neural network based on a genetic algorithm, as shown in Figure 14, the mean square error of the BP artificial neural network based on a genetic algorithm reaches its minimum of

Conclusions
Through the above experiments and simulation analysis, the conclusions are as follows: (1) The fitting accuracy rate is quite high and can improve the accuracy of MSMA actuators.
Because there are fewer undetermined parameters during the hysteresis curve fitting, the fitting speed of least squares method is fast.(2) Compared with BP artificial neural networks, BP artificial neural networks based on genetic algorithms can accelerate the convergence rate and decrease the fitting error.Therefore, they can improve the accuracy of MSMA actuators and enhance the utilization rate of To further show the advantage of the BP artificial neural network based on a genetic algorithm, as shown in Figure 14, the mean square error of the BP artificial neural network based on a genetic algorithm reaches its minimum of 5.0384 × 10 −5 at the 19th iteration.There is no doubt that a BP artificial neural network based on a genetic algorithm can accelerate the convergence rate and decrease the fitting error.To further show the advantage of the BP artificial neural network based on a genetic algorithm, as shown in Figure 14, the mean square error of the BP artificial neural network based on a genetic algorithm reaches its minimum of

Conclusions
Through the above experiments and simulation analysis, the conclusions are as follows: (1) The fitting accuracy rate is quite high and can improve the accuracy of MSMA actuators.

Conclusions
Through the above experiments and simulation analysis, the conclusions are as follows: (1) The fitting accuracy rate is quite high and can improve the accuracy of MSMA actuators.Because there are fewer undetermined parameters during the hysteresis curve fitting, the fitting speed of least squares method is fast.
(2) Compared with BP artificial neural networks, BP artificial neural networks based on genetic algorithms can accelerate the convergence rate and decrease the fitting error.Therefore, they can improve the accuracy of MSMA actuators and enhance the utilization rate of MSMA actuators in precision positioning.

Figure 1 .
Figure 1.Sample of magnetic controlled shape memory alloy (MSMA) for the experiment.

Figure 1 .
Figure 1.Sample of magnetic controlled shape memory alloy (MSMA) for the experiment.

( 5 )
Start the test.The corresponding deformation of MSMA is measured.(6) After unloading, view and save the experimental data.

Figure 2 .
Figure 2. Experimental setup: (1) power; (2) loading device; (3) power supply for temperature loading system; (4) computer with data acquisition card, which collects input magnetic field intensity H and output magnetic induction intensity B; and (5) control device for magnetic field.

Figure 2 .
Figure 2. Experimental setup: (1) power; (2) loading device; (3) power supply for temperature loading system; (4) computer with data acquisition card, which collects input magnetic field intensity H and output magnetic induction intensity B; and (5) control device for magnetic field.

Figure 3 .
Figure 3. Fitting curve and original experimental data curve."Processed 1" represents the fitting curve through least squares method."Original data" represents the original experimental data.

Figure 4 .
Figure 4.The fitting error rate curve of the least squares method.

Figure 3 .
Figure 3. Fitting curve and original experimental data curve."Processed 1" represents the fitting curve through least squares method."Original data" represents the original experimental data.

Figure 3 .
Figure 3. Fitting curve and original experimental data curve."Processed 1" represents the fitting curve through least squares method."Original data" represents the original experimental data.

Figure 4 .
Figure 4.The fitting error rate curve of the least squares method.

Figure 4 .
Figure 4.The fitting error rate curve of the least squares method.

Figure 5 .
Figure 5. Structure of the back propagation (BP) artificial neural network.

Figure 5 .
Figure 5. Structure of the back propagation (BP) artificial neural network.

Figure 6 .
Figure 6.Fitting curve and original experimental data curve."Processed 2" represents the fitting curve through the original BP artificial neural network."Original data" represents the original experimental data curve.

Figure 7 .
Figure 7. Fitting-error-rate curve for the original BP artificial neural network.

Figure 6 .
Figure 6.Fitting curve and original experimental data curve."Processed 2" represents the fitting curve through the original BP artificial neural network."Original data" represents the original experimental data curve.

Figure 6 .
Figure 6.Fitting curve and original experimental data curve."Processed 2" represents the fitting curve through the original BP artificial neural network."Original data" represents the original experimental data curve.

Figure 7 .
Figure 7. Fitting-error-rate curve for the original BP artificial neural network.

Figure 7 .
Figure 7. Fitting-error-rate curve for the original BP artificial neural network.

Figure 9 .
Figure 9. Flowchart of the improved original BP artificial neural network.

Figure 9 .
Figure 9. Flowchart of the improved original BP artificial neural network.

Figure 9 .
Figure 9. Flowchart of the improved original BP artificial neural network.

Figure 10 .
Figure 10.Fitting-error-rate curve."Original" and "improved" represent the fitting-error-rate curve of the BP artificial neural network and BP artificial neural network based on genetic algorithms, respectively.

Figure 10 .
Figure 10.Fitting-error-rate curve."Original" and "improved" represent the fitting-error-rate curve of the BP artificial neural network and BP artificial neural network based on genetic algorithms, respectively.

Figure 10 .
Figure 10.Fitting-error-rate curve."Original" and "improved" represent the fitting-error-rate curve of the BP artificial neural network and BP artificial neural network based on genetic algorithms, respectively.

Figure 11 .
Figure 11.Mean square error curve of the BP artificial neural network based on genetic algorithms.Figure 11.Mean square error curve of the BP artificial neural network based on genetic algorithms.

Figure 11 .
Figure 11.Mean square error curve of the BP artificial neural network based on genetic algorithms.Figure 11.Mean square error curve of the BP artificial neural network based on genetic algorithms.

Figure 12 .
Figure 12.Fitting curve and the original experimental data curve."Processed 4" represents the fitting curve through the BP neural network based on a genetic algorithm."Original data" represents the original experimental data.

Figure 12 .
Figure 12.Fitting curve and the original experimental data curve."Processed 4" represents the fitting curve through the BP neural network based on a genetic algorithm."Original data" represents the original experimental data.

Actuators 2016, 5 , 25 11 of 13 Figure 13 .
Figure 13.Fitting error rate curve of the BP artificial neural network and the BP artificial neural network based on genetic algorithms.

Figure 14 .
Figure 14.The mean square error curve of the BP artificial neural network based on a genetic algorithm.

5 5 .
0384 10 − × at the 19th iteration.There is no doubt that a BP artificial neural network based on a genetic algorithm can accelerate the convergence rate and decrease the fitting error.

Figure 13 .
Figure 13.Fitting error rate curve of the BP artificial neural network and the BP artificial neural network based on genetic algorithms.

Actuators 2016, 5 , 25 11 of 13 Figure 13 .
Figure 13.Fitting error rate curve of the BP artificial neural network and the BP artificial neural network based on genetic algorithms.

Figure 14 .
Figure 14.The mean square error curve of the BP artificial neural network based on a genetic algorithm.

5 5 .
0384 10 − × at the 19th iteration.There is no doubt that a BP artificial neural network based on a genetic algorithm can accelerate the convergence rate and decrease the fitting error.

Figure 14 .
Figure 14.The mean square error curve of the BP artificial neural network based on a genetic algorithm.

Table 1 .
Key original experimental data.