# Motion Path Design for Specific Muscle Training Using Neural Network

Authors: Kenta Itokazu; Takanori Miyoshi; Kazuhiko Terashima
Journal: Journal of Robotics
 (2013)

---

## Abstract

Specific muscle training is expected to be used for efficient rehabilitation and care prevention. In this paper, we propose algorithms for designing a motion path capable of strengthening specific muscles. By using the proposed algorithms, it is possible to design a motion path maximizing the activity of an agonist muscle and minimizing that of other muscles. For training, the load is applied by using a 2-link arm. EMG signal is measured during a training experiment, and the degree of muscular revitalization is evaluated by the amplitude of EMG signal. Finally, the effectiveness of the proposed approach is demonstrated through experiments.


---

## Body

## 1. Introduction

In recent years, in the context of the emergence of population ageing as a social issue in developed countries, the importance of regaining muscle strength for care prevention has become increasingly apparent. Prolonged immobility induces muscle weakness, which affects activities of daily living (ADLs) directly.Much research is being done on rehabilitation robotics that is pertinent to strength training. Lum et al. indicated that, compared with conventional therapy techniques, robot-assisted training is more efficient for improving muscle strength and path-following capability [1]. For lower limb rehabilitation, Akdoğan and Adli developed a therapeutic exercise robot that enables rehabilitation for spinal cord injury in diverse ways, including both isotonically and isometrically [2]. Such research aims to regain a muscle strength of an entire arm or leg. Therefore, the robots are unable to apply a load to specific muscles. However, the degree of muscle weakness differs according to each muscle. Thus, the application of a load to specific muscles that require strengthening is expected to lead to more efficient and safer training.In the authors' previous research, a method of estimating muscle force or level of muscle activation was proposed [3]. Though the objective of the research was to strengthen a muscle by isometric exercise, an isotonic exercise is more effective for ADLs training. The purpose of the present work is to develop an algorithm for designing a motion path capable of strengthening specific muscles for isotonic exercise.The method proposed in previous research does not consider the coordinated motion of an antagonistic muscle. However, when doing an exercise, an antagonistic muscle works to increase the stiffness of each joint, such as a shoulder or an elbow. A method of estimating muscle activity should consider the coordinated motion. In this paper, a neural network is used since it is suited to estimation of nonlinear data such as muscle activity. First, the neural network is trained by a backpropagation algorithm. A training data differs according to each subject person. The neural network is able to estimate the level of muscle activation. Secondly, we design the optimal motion path by using a multiobjective optimization method for the obtained neural network. The objective of optimization is to maximize the activity of an agonist muscle and minimize the activity of other muscles. The motion paths obtained by an optimization method for the learned neural network model are evaluated through experiments, and the validity of the proposed approach is then demonstrated.


## 2. System Description

### 2.1. System Structure

Figure1 shows overall look of the 2-link arm that is used as a training device. There are DC motors and rotary encoders at each joint, and 6-axis force sensor at the end-point. An emergency stop switch is mounted. The device enables training of upper limbs by horizontal motion. Figure 2 shows appearance of an experiment. The EMG sensors are attached to the subject's arm.Figure 1
2-link arm.Figure 2
Appearance of experiment.


### 2.2. Control System

An end-point velocity of the device is controlled by force feedback. To apply a load to muscles, virtual mass and viscosity are set at the end-point (see Figure3).Figure 3
Virtual mass and viscosity.A motion equation of the end-point is described by(1)m·a+D·v=f,
where m is virtual mass, D is virtual viscosity coefficient, f=(fx,fy)T is force measured by force sensor, and a=(x¨,y¨)T and v=(x˙,y˙)T are acceleration and velocity of the end-point, respectively.The control input for motors is given by(2)vi+1=a·dt+vi,Θ˙i+1=Ji+1-1vi+1,J=[l1S1+l2S12l2S12l1C1+l2C12l2C12],s=KΘ˙i+1,
where dt is a sampling time set as 1 [msec], Θ˙=(θ˙1,θ˙2)T is angular velocity of each joint, and S1, C1, S12, and C12 are sinθ1, cosθ1, sin(θ1+θ2), and cos(θ1+θ2), respectively. Furthermore, s=(s1,s2)T is input signal, K is conversion coefficient, and l1 and l2 are length of upper arm and lower arm that are set as 0.26 and 0.3 [m], respectively.


## 2.1. System Structure

Figure1 shows overall look of the 2-link arm that is used as a training device. There are DC motors and rotary encoders at each joint, and 6-axis force sensor at the end-point. An emergency stop switch is mounted. The device enables training of upper limbs by horizontal motion. Figure 2 shows appearance of an experiment. The EMG sensors are attached to the subject's arm.Figure 1
2-link arm.Figure 2
Appearance of experiment.


## 2.2. Control System

An end-point velocity of the device is controlled by force feedback. To apply a load to muscles, virtual mass and viscosity are set at the end-point (see Figure3).Figure 3
Virtual mass and viscosity.A motion equation of the end-point is described by(1)m·a+D·v=f,
where m is virtual mass, D is virtual viscosity coefficient, f=(fx,fy)T is force measured by force sensor, and a=(x¨,y¨)T and v=(x˙,y˙)T are acceleration and velocity of the end-point, respectively.The control input for motors is given by(2)vi+1=a·dt+vi,Θ˙i+1=Ji+1-1vi+1,J=[l1S1+l2S12l2S12l1C1+l2C12l2C12],s=KΘ˙i+1,
where dt is a sampling time set as 1 [msec], Θ˙=(θ˙1,θ˙2)T is angular velocity of each joint, and S1, C1, S12, and C12 are sinθ1, cosθ1, sin(θ1+θ2), and cos(θ1+θ2), respectively. Furthermore, s=(s1,s2)T is input signal, K is conversion coefficient, and l1 and l2 are length of upper arm and lower arm that are set as 0.26 and 0.3 [m], respectively.


## 3. Neural Network

### 3.1. Outline

In this study, three-layer Artificial Neural Network (ANN) consisting of an input layer, a hidden layer, and an output layer is used, as shown in Figure4. The inputs of ANN are joint torque (T1, T2) and joint angle (θ1, θ2), as shown in Figure 5. And, the outputs are level of muscle activation described as α^1,…,α^5.α1,…,α5 are measured values.Figure 4
Neural network.Figure 5
Upper limb model.The level of muscle activationα1,…,α5 is defined as
(3)αi(t)=EMGi(t)maxEMGi(i=1~5),
where EMGi(t) is electromyographic (EMG) signal that is measured during an experiment. maxEMGi is maximum EMG signal for each muscle. During an experiment, EMG signals of muscles shown in Table 1 are measured.Table 1
Measured muscles.
Name
M
1
Pectoralis major
M
2
Latissimus dorsi
M
3
Brachioradialis
M
4
Triceps brachii Lateral head
M
5
Biceps brachiiThe number of hidden units is determined by comparing root-mean-square errors (RMSEs) of estimation results. RMSEs are calculated by(4)RMSE=1N∑i=1N‍(α^i-αi)2.Figure6 shows RMSEs of estimation results obtained using ANN that is trained by training data-sets as described in Section 3.2. Even where there are more than 10 units, RMSE is virtually unchanged. Therefore, we set the number as 10.Figure 6
RMSE for various number of units in hidden layer.The outputs of ANN are calculated by(5)Sj(1)=∑n=1N=4‍wn,j(1)xn(1)(j=1~10),yj(1)=sigmoid(Sj(1)+bj(1)),Si(2)=∑n=1N=10‍wn,i(2)yn(1)(i=1~5),αi=sigmoid(Si(2)+bi(2)),sigmoid(z)=11+e-z,
where x is input (joint torque and angle), α is output (level of muscle activation), w is connection weight, and b is bias.


### 3.2. Training of Neural Network

In order to determine connection weights and biases, training of ANN is required. The training data sets are obtained through exploratory experiments for each subject. The data is measured during linear motion on horizontal plane as shown in Figure7. During the experiments, the motion paths, force, and position of the end-point are displayed on a monitor. The displayed force and end-point position are updated in a real time. And then, subjects follow the paths at a constant force (30 [N]). The ANN is trained by a backpropagation algorithm.Figure 7
Motion path for exploratory experiment.


## 3.1. Outline

In this study, three-layer Artificial Neural Network (ANN) consisting of an input layer, a hidden layer, and an output layer is used, as shown in Figure4. The inputs of ANN are joint torque (T1, T2) and joint angle (θ1, θ2), as shown in Figure 5. And, the outputs are level of muscle activation described as α^1,…,α^5.α1,…,α5 are measured values.Figure 4
Neural network.Figure 5
Upper limb model.The level of muscle activationα1,…,α5 is defined as
(3)αi(t)=EMGi(t)maxEMGi(i=1~5),
where EMGi(t) is electromyographic (EMG) signal that is measured during an experiment. maxEMGi is maximum EMG signal for each muscle. During an experiment, EMG signals of muscles shown in Table 1 are measured.Table 1
Measured muscles.
Name
M
1
Pectoralis major
M
2
Latissimus dorsi
M
3
Brachioradialis
M
4
Triceps brachii Lateral head
M
5
Biceps brachiiThe number of hidden units is determined by comparing root-mean-square errors (RMSEs) of estimation results. RMSEs are calculated by(4)RMSE=1N∑i=1N‍(α^i-αi)2.Figure6 shows RMSEs of estimation results obtained using ANN that is trained by training data-sets as described in Section 3.2. Even where there are more than 10 units, RMSE is virtually unchanged. Therefore, we set the number as 10.Figure 6
RMSE for various number of units in hidden layer.The outputs of ANN are calculated by(5)Sj(1)=∑n=1N=4‍wn,j(1)xn(1)(j=1~10),yj(1)=sigmoid(Sj(1)+bj(1)),Si(2)=∑n=1N=10‍wn,i(2)yn(1)(i=1~5),αi=sigmoid(Si(2)+bi(2)),sigmoid(z)=11+e-z,
where x is input (joint torque and angle), α is output (level of muscle activation), w is connection weight, and b is bias.


## 3.2. Training of Neural Network

In order to determine connection weights and biases, training of ANN is required. The training data sets are obtained through exploratory experiments for each subject. The data is measured during linear motion on horizontal plane as shown in Figure7. During the experiments, the motion paths, force, and position of the end-point are displayed on a monitor. The displayed force and end-point position are updated in a real time. And then, subjects follow the paths at a constant force (30 [N]). The ANN is trained by a backpropagation algorithm.Figure 7
Motion path for exploratory experiment.


## 4. Motion Path Design

### 4.1. Outline

To design motion paths, we proposed two algorithms(i)
algorithm for designing a motion path that always has the same initial posture (Algorithm 1),(ii)
algorithm for designing a motion path that passes through the most effective point (Algorithm 2).These algorithms use a multiobjective problem. When designing a path using Algorithm 1, the initial value of multiobjective problem is the same for each time. Therefore, the designed paths for each muscles (M1~M5) have the same initial position. The motion path designed by Algorithm 1 does not necessarily pass through the most effective point that is capable of maximizing the activity of an agonist muscle and minimizing the activity of other muscles. By contrast, Algorithm 2 designs a motion path based on the most effective point.


### 4.2. Algorithm 1

The algorithm uses ANN and multiobjective optimization to design motion paths. Figure8 shows a flowchart of the algorithm.Figure 8
Flowchart of Algorithm 1.First, define an initial posture asθ1=π/4 [rad] and θ2=π/2 [rad]. Next, solve an optimization problem as shown below to satisfy both:
(6)maxθ000αn(θ)minθ000α-(θ),σ2(θ)subjectto0≤θ<2π,
where θ is direction of force, αn is level of muscle activation for target muscle, and α- and σ2 are average and variance of nontarget muscles. αn, α- and σ2 are calculated by using ANN. The aim of this optimization problem is to maximize the activity of an agonist muscle and minimize the activity of other muscles. In order to solve the optimization problem, sequential quadratic programming method is used. The method is also used to solve (8) and (9). To calculate the velocity from end-point force, (1) is simplified by setting m=0. By solving the optimization problem, we obtain the optimal direction of force. Thirdly, calculate the position of the end-point. The position is calculated by
(7)f=[FxFy]=[|F|cos(θ)|F|sin(θ)],v=[VxVy]=[FxDFyD],Θ˙=[θ˙1θ˙2]=J-1v,Θ=[θ1θ2]=∫‍Θ˙dt,T=[T1T2]=JTf,J=[l1S1+l2S12l2S12l1C1+l2C12l2C12],P=[xy]=[-l1C1-l2C12l1S1+l2S12],
where F is end-point force, |F| is absolute value of end-point force that is set as 30 [N], V is end-point velocity, P is position of end-point, Θ˙ and Θ are angular velocity and angle of joint, and T is joint torque. And dt is sampling time that is set as 1 [msec], D is virtual viscosity coefficient that is set as 900 [Ns/m], and l1 and l2 are length of upper arm and lower arm that are set as 0.26 and 0.3 [m].By repeating the procedure until the distance between end-point position and origin exceeds 0.2 [m] or level of muscle activation of target muscle is less than 0.5, a motion path is designed. Figure10 shows an example of a designed motion path. P1~P5 express motion paths of strength training that target each muscle (e.g., P1 targets M1). In this algorithm, these paths have the same initial position.M
1
~
M
5 have different function of arm movement. The shown paths are designed to enhance muscle activation of target muscle. Therefore, direction of the paths depend on the function of each muscle. For example, P1 shows adduction movement of shoulder, because M1 has function of shoulder adduction.


### 4.3. Algorithm 2

In Algorithm 2, most effective point that is capable of maximizing the activity of an agonist muscle and minimizing the activity of other muscles is calculated before designing a motion path to use in (9) as an initial value.A posture to take the most effective point is given by(8)maxθ,θ1,θ2000000αn(θ,θ1,θ2)minθ,θ1,θ2000000α-(θ,θ1,θ2),σ2(θ,θ1,θ2)subjectto0000≤θ<2π00.00000vv0-π6≤θ1≤135π1800000000vvv00≤θ2≤145π180.These variables are the same as for (6). When solving (8), initial posture is set as θ1=π/4 [rad] and θ2=π/2 [rad].An end-point force is set as 30 [N] in (8). However, after the beginning of a real experiment, the end-point force is lower than 30 [N] for a moment. Therefore, using the value obtained by (8) as initial posture of designed path is unsuitable. Thus, the most effective point obtained by (8) is passed through in mid-course of designed path in order to fit a situation of (8).In this algorithm, we design two paths based on the posture obtained by (8). And then, the two paths are integrated to one path that passes through the most effective point in mid-course.First, the value obtained by (8) is defined as the initial value of (9). Secondly, an optimization problem is solved as follows:
(9)maxθαn(θ)minθα-(θ),σ2(θ)subjecttolb≤θ<ub,
where lb and ub are lower bound and upper bound that are set as lb=θ-π/2 and ub=θ+π/2. θ used in lb and ub is constant value obtained by (8). Thirdly, a path is designed by the flow described in Figure 9 using (7). The end condition is the same as Algorithm 1. Next, the values of lb and ub are switched, and the value obtained by (8) is set as initial value. And then, another path is designed by the same flow. To integrate the two paths into one path, these paths must not have opposite directions. Thus, in this case, the force is reversed when it is used for calculating joint torques that are inputs of ANN in order to align the direction of the paths. Finally, the two paths are integrated to one path. Figure 11 shows an example of a designed motion path. As described in Section 4.2, these paths depend on the function of each muscle.Figure 9
Flowchart of Algorithm 2.Figure 10
One example of motion path (Algorithm 1).Figure 11
One example of motion path (Algorithm 2).


## 4.1. Outline

To design motion paths, we proposed two algorithms(i)
algorithm for designing a motion path that always has the same initial posture (Algorithm 1),(ii)
algorithm for designing a motion path that passes through the most effective point (Algorithm 2).These algorithms use a multiobjective problem. When designing a path using Algorithm 1, the initial value of multiobjective problem is the same for each time. Therefore, the designed paths for each muscles (M1~M5) have the same initial position. The motion path designed by Algorithm 1 does not necessarily pass through the most effective point that is capable of maximizing the activity of an agonist muscle and minimizing the activity of other muscles. By contrast, Algorithm 2 designs a motion path based on the most effective point.


## 4.2. Algorithm 1

The algorithm uses ANN and multiobjective optimization to design motion paths. Figure8 shows a flowchart of the algorithm.Figure 8
Flowchart of Algorithm 1.First, define an initial posture asθ1=π/4 [rad] and θ2=π/2 [rad]. Next, solve an optimization problem as shown below to satisfy both:
(6)maxθ000αn(θ)minθ000α-(θ),σ2(θ)subjectto0≤θ<2π,
where θ is direction of force, αn is level of muscle activation for target muscle, and α- and σ2 are average and variance of nontarget muscles. αn, α- and σ2 are calculated by using ANN. The aim of this optimization problem is to maximize the activity of an agonist muscle and minimize the activity of other muscles. In order to solve the optimization problem, sequential quadratic programming method is used. The method is also used to solve (8) and (9). To calculate the velocity from end-point force, (1) is simplified by setting m=0. By solving the optimization problem, we obtain the optimal direction of force. Thirdly, calculate the position of the end-point. The position is calculated by
(7)f=[FxFy]=[|F|cos(θ)|F|sin(θ)],v=[VxVy]=[FxDFyD],Θ˙=[θ˙1θ˙2]=J-1v,Θ=[θ1θ2]=∫‍Θ˙dt,T=[T1T2]=JTf,J=[l1S1+l2S12l2S12l1C1+l2C12l2C12],P=[xy]=[-l1C1-l2C12l1S1+l2S12],
where F is end-point force, |F| is absolute value of end-point force that is set as 30 [N], V is end-point velocity, P is position of end-point, Θ˙ and Θ are angular velocity and angle of joint, and T is joint torque. And dt is sampling time that is set as 1 [msec], D is virtual viscosity coefficient that is set as 900 [Ns/m], and l1 and l2 are length of upper arm and lower arm that are set as 0.26 and 0.3 [m].By repeating the procedure until the distance between end-point position and origin exceeds 0.2 [m] or level of muscle activation of target muscle is less than 0.5, a motion path is designed. Figure10 shows an example of a designed motion path. P1~P5 express motion paths of strength training that target each muscle (e.g., P1 targets M1). In this algorithm, these paths have the same initial position.M
1
~
M
5 have different function of arm movement. The shown paths are designed to enhance muscle activation of target muscle. Therefore, direction of the paths depend on the function of each muscle. For example, P1 shows adduction movement of shoulder, because M1 has function of shoulder adduction.


## 4.3. Algorithm 2

In Algorithm 2, most effective point that is capable of maximizing the activity of an agonist muscle and minimizing the activity of other muscles is calculated before designing a motion path to use in (9) as an initial value.A posture to take the most effective point is given by(8)maxθ,θ1,θ2000000αn(θ,θ1,θ2)minθ,θ1,θ2000000α-(θ,θ1,θ2),σ2(θ,θ1,θ2)subjectto0000≤θ<2π00.00000vv0-π6≤θ1≤135π1800000000vvv00≤θ2≤145π180.These variables are the same as for (6). When solving (8), initial posture is set as θ1=π/4 [rad] and θ2=π/2 [rad].An end-point force is set as 30 [N] in (8). However, after the beginning of a real experiment, the end-point force is lower than 30 [N] for a moment. Therefore, using the value obtained by (8) as initial posture of designed path is unsuitable. Thus, the most effective point obtained by (8) is passed through in mid-course of designed path in order to fit a situation of (8).In this algorithm, we design two paths based on the posture obtained by (8). And then, the two paths are integrated to one path that passes through the most effective point in mid-course.First, the value obtained by (8) is defined as the initial value of (9). Secondly, an optimization problem is solved as follows:
(9)maxθαn(θ)minθα-(θ),σ2(θ)subjecttolb≤θ<ub,
where lb and ub are lower bound and upper bound that are set as lb=θ-π/2 and ub=θ+π/2. θ used in lb and ub is constant value obtained by (8). Thirdly, a path is designed by the flow described in Figure 9 using (7). The end condition is the same as Algorithm 1. Next, the values of lb and ub are switched, and the value obtained by (8) is set as initial value. And then, another path is designed by the same flow. To integrate the two paths into one path, these paths must not have opposite directions. Thus, in this case, the force is reversed when it is used for calculating joint torques that are inputs of ANN in order to align the direction of the paths. Finally, the two paths are integrated to one path. Figure 11 shows an example of a designed motion path. As described in Section 4.2, these paths depend on the function of each muscle.Figure 9
Flowchart of Algorithm 2.Figure 10
One example of motion path (Algorithm 1).Figure 11
One example of motion path (Algorithm 2).


## 5. Results of Experiments

### 5.1. Outline

The effectiveness of the proposed approach is evaluated through experiments. An outline of the experiments is as follows:(i)
comparative experiments of Algorithms 1 and 2.(ii)
comparative experiments of personalized path and nonpersonalized paths.Table2 shows body dimension data of subjects. All subjects are male.Table 2
Body dimension.
Sub.
Age
Height [cm]
Weight [kg]
Length of lower arm [cm]
Length of upper arm [cm]
A
23
161
48
23
28
B
23
167
60
25
31
C
22
170
51
25
29
D
21
163
52
25
31
E
24
171
60
26
32
F
23
168
54
25
27
G
24
170
62
24
26
H
24
174
94
27
28
I
24
170
62
24
27
J
22
179
69
26
29


### 5.2. Comparative Experiments of Algorithms 1 and 2 (Experiment 1)

#### 5.2.1. Experiment Objective

The objective of the experiment is to evaluate the results of each paths designed by Algorithms 1 and 2. In this experiment, ANN is trained by own data of each subject person.


#### 5.2.2. Results

The experimental procedure is the same as the exploratory experiments described in Section3.2. During the experiments, the designed motion paths, force, and position of end-point are displayed on a monitor.The results of the experiments are shown in Figures12, 13, 14, 15, and 16. For example, Figure 12 shows the result of P1 that is capable of strengthening M1. In Figures 12–16, “Measured” means measured value during the experiment. And “Predicted” means output of ANN calculated before the experiment. The shown values are average of the results of ten subjects.Figure 12
Level of muscle activation (Target:M1).Figure 13
Level of muscle activation (Target:M2).Figure 14
Level of muscle activation (Target:M3).Figure 15
Level of muscle activation (Target:M4).Figure 16
Level of muscle activation (Target:M5).There is a difference between the predicted value and the measured values. It is a consequence of redundancy of muscles. It is not necessarily always the same level of activation even if doing the same motion. However, it is possible to predict which muscle is the agonist muscle in most cases. Therefore, ANN is capable of designing a motion path that achieve the objective of this research.


### 5.3. Comparative Experiments of Personalized Path and Other Paths (Experiment 2)

#### 5.3.1. Experiment Objective

The objective of the experiment is to evaluate an influence exerted by personal experimental data. The motion paths are obtained using ANN. And ANN has to be trained by experimental data before using for path design. Thus, there are personalized path and nonpersonalized path. The personalized path is designed using ANN trained by own experimental data. And the nonpersonalized path uses ANN trained by experimental data of other subjects. The objective of the experiment is to compare a personalized path and other paths. In this experiment, the paths are designed by Algorithm 2.


#### 5.3.2. Results

Figures17, 18, 19, 20, 21, 22, 23, 24, 25, and 26 show experimental results of each subject. The level of muscle activation of the targeted muscle is shown in the upper portion of the graph, and the average of muscle activation of nontargeted muscles is shown in the lower portion. The number of motion path is equal to the number of subject. Therefore, nine paths are nonpersonalized path for each subject. The values of “Nonpersonalized” shown in Figures 17–26 are average of experimental results using the nine paths.Level of muscle activation (Sub. A).
(a)
(b)Level of muscle activation (Sub. B).
(a)
(b)Level of muscle activation (Sub. C).
(a)
(b)Level of muscle activation (Sub. D).
(a)
(b)Level of muscle activation (Sub. E).
(a)
(b)Level of muscle activation (Sub. F).
(a)
(b)Level of muscle activation (Sub. G).
(a)
(b)Level of muscle activation (Sub. H).
(a)
(b)Level of muscle activation (Sub. I).
(a)
(b)Level of muscle activation (Sub. J).
(a)
(b)


### 5.4. Discussion

To evaluate the result of Experiment1 (Section 5.2) quantitatively, we calculated the improvement rate of muscle activation. The value is given by
(10)Rallα=1S1N∑i=1S=10‍∑j=1N=5‍αa,jiαb,ji-1,Rallα-=1S1N∑i=1S=10‍∑j=1N=5‍α-a,jiα-b,ji-1,RMjα=1S∑i=1S=10‍αa,jiαb,ji-1,RMjα-=1S∑i=1S=10‍α-a,jiα-b,ji-1,
where α is level of muscle activation of target muscle, α- is average of muscle activation of nontargeted muscle, a is Algorithm 2, and b is Algorithm 1. i means a number of subjects, j is number of muscles. Rα is improvement rate for level of muscle activation of targeted muscle and Rα- is improvement rate for average of muscle activation of nontargeted muscle, respectively.The results of the calculation is shown in Table3. As shown in Table 3, RM1α is small value. The results indicate that there is little difference between Algorithms 1 and 2 when designing a path for M1. By contrast, RM5α shows remarkable improvement compared to other results. Furthermore, RM5α- shows that average of muscle activation of nontargeted muscle is decreased drastically. The results show that Algorithm 2 is more effective than Algorithm 1 for M5. The difference between M5 and others is classification of muscle. M5 is biarticular muscle. On the other hand, others are monoarticular muscle. It would appear that the difference affects experimental results.Table 3
Improvement rate (Experiment1).
Improvement rate [%]
All
M
1
M
2
M
3
M
4
M
5
R
α
20.9
−0.28
12.3
22.8
11.2
58.4
R
α
-
−15.6
−6.52
−14.5
−15.5
−2.12
−39.5Next, in order to evaluate the result of Experiment 2 (Section5.3), we calculated the ratio of muscle activation by (10). In this case, αa, αb, α-a, and α-b are replaced with αp, αn, α-p, and α-n. p means personalized path, and n means nonpersonalized path. The results of the calculation are shown in Table 4. If all improvement rates were 0%, the result shows that it would make no difference between personalized path and nonpersonalized path. However, all value of Rα is positive and Rα- is negative. Thus, the results indicate that personalized path is more effective than nonpersonalized path in all cases. Therefore, the training of ANN by own experimental data is needed before path design.Table 4
Improvement rate (Experiment2).
Improvement rate [%]
All
M
1
M
2
M
3
M
4
M
5
R
α
41.4
27.9
44.2
47.7
21.2
65.7
R
α
-
−20.1
−20.1
−30.1
−14.4
−13.8
−22.1


## 5.1. Outline

The effectiveness of the proposed approach is evaluated through experiments. An outline of the experiments is as follows:(i)
comparative experiments of Algorithms 1 and 2.(ii)
comparative experiments of personalized path and nonpersonalized paths.Table2 shows body dimension data of subjects. All subjects are male.Table 2
Body dimension.
Sub.
Age
Height [cm]
Weight [kg]
Length of lower arm [cm]
Length of upper arm [cm]
A
23
161
48
23
28
B
23
167
60
25
31
C
22
170
51
25
29
D
21
163
52
25
31
E
24
171
60
26
32
F
23
168
54
25
27
G
24
170
62
24
26
H
24
174
94
27
28
I
24
170
62
24
27
J
22
179
69
26
29


## 5.2. Comparative Experiments of Algorithms 1 and 2 (Experiment 1)

### 5.2.1. Experiment Objective

The objective of the experiment is to evaluate the results of each paths designed by Algorithms 1 and 2. In this experiment, ANN is trained by own data of each subject person.


### 5.2.2. Results

The experimental procedure is the same as the exploratory experiments described in Section3.2. During the experiments, the designed motion paths, force, and position of end-point are displayed on a monitor.The results of the experiments are shown in Figures12, 13, 14, 15, and 16. For example, Figure 12 shows the result of P1 that is capable of strengthening M1. In Figures 12–16, “Measured” means measured value during the experiment. And “Predicted” means output of ANN calculated before the experiment. The shown values are average of the results of ten subjects.Figure 12
Level of muscle activation (Target:M1).Figure 13
Level of muscle activation (Target:M2).Figure 14
Level of muscle activation (Target:M3).Figure 15
Level of muscle activation (Target:M4).Figure 16
Level of muscle activation (Target:M5).There is a difference between the predicted value and the measured values. It is a consequence of redundancy of muscles. It is not necessarily always the same level of activation even if doing the same motion. However, it is possible to predict which muscle is the agonist muscle in most cases. Therefore, ANN is capable of designing a motion path that achieve the objective of this research.


## 5.2.1. Experiment Objective

The objective of the experiment is to evaluate the results of each paths designed by Algorithms 1 and 2. In this experiment, ANN is trained by own data of each subject person.


## 5.2.2. Results

The experimental procedure is the same as the exploratory experiments described in Section3.2. During the experiments, the designed motion paths, force, and position of end-point are displayed on a monitor.The results of the experiments are shown in Figures12, 13, 14, 15, and 16. For example, Figure 12 shows the result of P1 that is capable of strengthening M1. In Figures 12–16, “Measured” means measured value during the experiment. And “Predicted” means output of ANN calculated before the experiment. The shown values are average of the results of ten subjects.Figure 12
Level of muscle activation (Target:M1).Figure 13
Level of muscle activation (Target:M2).Figure 14
Level of muscle activation (Target:M3).Figure 15
Level of muscle activation (Target:M4).Figure 16
Level of muscle activation (Target:M5).There is a difference between the predicted value and the measured values. It is a consequence of redundancy of muscles. It is not necessarily always the same level of activation even if doing the same motion. However, it is possible to predict which muscle is the agonist muscle in most cases. Therefore, ANN is capable of designing a motion path that achieve the objective of this research.


## 5.3. Comparative Experiments of Personalized Path and Other Paths (Experiment 2)

### 5.3.1. Experiment Objective

The objective of the experiment is to evaluate an influence exerted by personal experimental data. The motion paths are obtained using ANN. And ANN has to be trained by experimental data before using for path design. Thus, there are personalized path and nonpersonalized path. The personalized path is designed using ANN trained by own experimental data. And the nonpersonalized path uses ANN trained by experimental data of other subjects. The objective of the experiment is to compare a personalized path and other paths. In this experiment, the paths are designed by Algorithm 2.


### 5.3.2. Results

Figures17, 18, 19, 20, 21, 22, 23, 24, 25, and 26 show experimental results of each subject. The level of muscle activation of the targeted muscle is shown in the upper portion of the graph, and the average of muscle activation of nontargeted muscles is shown in the lower portion. The number of motion path is equal to the number of subject. Therefore, nine paths are nonpersonalized path for each subject. The values of “Nonpersonalized” shown in Figures 17–26 are average of experimental results using the nine paths.Level of muscle activation (Sub. A).
(a)
(b)Level of muscle activation (Sub. B).
(a)
(b)Level of muscle activation (Sub. C).
(a)
(b)Level of muscle activation (Sub. D).
(a)
(b)Level of muscle activation (Sub. E).
(a)
(b)Level of muscle activation (Sub. F).
(a)
(b)Level of muscle activation (Sub. G).
(a)
(b)Level of muscle activation (Sub. H).
(a)
(b)Level of muscle activation (Sub. I).
(a)
(b)Level of muscle activation (Sub. J).
(a)
(b)


## 5.3.1. Experiment Objective

The objective of the experiment is to evaluate an influence exerted by personal experimental data. The motion paths are obtained using ANN. And ANN has to be trained by experimental data before using for path design. Thus, there are personalized path and nonpersonalized path. The personalized path is designed using ANN trained by own experimental data. And the nonpersonalized path uses ANN trained by experimental data of other subjects. The objective of the experiment is to compare a personalized path and other paths. In this experiment, the paths are designed by Algorithm 2.


## 5.3.2. Results

Figures17, 18, 19, 20, 21, 22, 23, 24, 25, and 26 show experimental results of each subject. The level of muscle activation of the targeted muscle is shown in the upper portion of the graph, and the average of muscle activation of nontargeted muscles is shown in the lower portion. The number of motion path is equal to the number of subject. Therefore, nine paths are nonpersonalized path for each subject. The values of “Nonpersonalized” shown in Figures 17–26 are average of experimental results using the nine paths.Level of muscle activation (Sub. A).
(a)
(b)Level of muscle activation (Sub. B).
(a)
(b)Level of muscle activation (Sub. C).
(a)
(b)Level of muscle activation (Sub. D).
(a)
(b)Level of muscle activation (Sub. E).
(a)
(b)Level of muscle activation (Sub. F).
(a)
(b)Level of muscle activation (Sub. G).
(a)
(b)Level of muscle activation (Sub. H).
(a)
(b)Level of muscle activation (Sub. I).
(a)
(b)Level of muscle activation (Sub. J).
(a)
(b)


## 5.4. Discussion

To evaluate the result of Experiment1 (Section 5.2) quantitatively, we calculated the improvement rate of muscle activation. The value is given by
(10)Rallα=1S1N∑i=1S=10‍∑j=1N=5‍αa,jiαb,ji-1,Rallα-=1S1N∑i=1S=10‍∑j=1N=5‍α-a,jiα-b,ji-1,RMjα=1S∑i=1S=10‍αa,jiαb,ji-1,RMjα-=1S∑i=1S=10‍α-a,jiα-b,ji-1,
where α is level of muscle activation of target muscle, α- is average of muscle activation of nontargeted muscle, a is Algorithm 2, and b is Algorithm 1. i means a number of subjects, j is number of muscles. Rα is improvement rate for level of muscle activation of targeted muscle and Rα- is improvement rate for average of muscle activation of nontargeted muscle, respectively.The results of the calculation is shown in Table3. As shown in Table 3, RM1α is small value. The results indicate that there is little difference between Algorithms 1 and 2 when designing a path for M1. By contrast, RM5α shows remarkable improvement compared to other results. Furthermore, RM5α- shows that average of muscle activation of nontargeted muscle is decreased drastically. The results show that Algorithm 2 is more effective than Algorithm 1 for M5. The difference between M5 and others is classification of muscle. M5 is biarticular muscle. On the other hand, others are monoarticular muscle. It would appear that the difference affects experimental results.Table 3
Improvement rate (Experiment1).
Improvement rate [%]
All
M
1
M
2
M
3
M
4
M
5
R
α
20.9
−0.28
12.3
22.8
11.2
58.4
R
α
-
−15.6
−6.52
−14.5
−15.5
−2.12
−39.5Next, in order to evaluate the result of Experiment 2 (Section5.3), we calculated the ratio of muscle activation by (10). In this case, αa, αb, α-a, and α-b are replaced with αp, αn, α-p, and α-n. p means personalized path, and n means nonpersonalized path. The results of the calculation are shown in Table 4. If all improvement rates were 0%, the result shows that it would make no difference between personalized path and nonpersonalized path. However, all value of Rα is positive and Rα- is negative. Thus, the results indicate that personalized path is more effective than nonpersonalized path in all cases. Therefore, the training of ANN by own experimental data is needed before path design.Table 4
Improvement rate (Experiment2).
Improvement rate [%]
All
M
1
M
2
M
3
M
4
M
5
R
α
41.4
27.9
44.2
47.7
21.2
65.7
R
α
-
−20.1
−20.1
−30.1
−14.4
−13.8
−22.1


## 6. Conclusion

In this study, we developed two algorithms for designing a motion path capable of strengthening specific muscles. Algorithm 1 designs a path that always has the same initial posture, and Algorithm 2 designs a path that has free initial posture. As a result of the experiments, Algorithm 2 is found to be more effective than Algorithm 1 from the viewpoint of the objective of this study. The path designed by Algorithm 1 always has the same initial posture. Therefore, the initial condition is stricter than that for Algorithm 2. The experimental result shows that less strict constraint condition leads to more effective result.Furthermore, the results show that a personalized path produces a better outcome than nonpersonalized paths. This indicates that ANN trained by own experimental data is capable of designing a more effective motion path.The objective of this study is to strengthen upper limb muscles. However, ADLs depend on not only upper limb muscles but also lower limb muscles. In future work, we intended to apply the proposed approach to training of lower limb muscles in preliminary rehabilitation for walking.



---

*Source: 810909-2013-12-04.xml*