# Neuroendocrine-Based Cooperative Intelligent Control System for Multiobjective Integrated Control of a Parallel Manipulator

Authors: Chongbin Guo; Kuangrong Hao; Yongsheng Ding
Journal: Mathematical Problems in Engineering
 (2012)

---

## Abstract

This paper presents a novel multiloop and Multi-objective cooperative intelligent control system (MMCICS) used to improve the performance of position, velocity and acceleration integrated control on a complex multichannel plant. Based on regulation mechanism of the neuroendocrine system (NES), a bioinspired motion control approach has been used in the MMCICS which includes four cooperative units. The planning unit outputs the desired signals. The selection unit chooses the real-time dominant control mode. The coordination unit uses the velocity Jacobian matrix to regulate the cooperative control signals. The execution unit achieves the ultimate task based on sub-channel controllers with the proposed hormone regulation self-adaptive Modules (HRSMs). Parameter tuning is given to facilitate the MMCICS implementation. The MMCICS is applied to an actual 2-DOF redundant parallel manipulator where the feasibility of the new control system is demonstrated. The MMCICS keeps its subchannels interacting harmoniously and systematically. Therefore, the plant has fast response, smooth velocity, accurate position, strong self-adaptability, and high stability. The HRSM improves the control performance of the local controllers and the global system as well, especially for manipulators running at high velocities and accelerations.


---

## Body

## 1. Introduction

With the development of the high-standard manufacturing requirement, plants become more complex while controlled by several of subchannels [1, 2]. Usually, the different sub-channels have different characteristics and control requirements. Therefore, the sub-channels have to interact harmoniously and systematically to achieve multiobjective integrated control [3, 4]. In that manner, the plants can have a quick start and stop, a fine uniform movement, an accurate destination, and a strong self-adaptability with stability [5]. In general, position based control cannot keep uniform velocity, while velocity-based control cannot satisfy accurate position requirement [6, 7]. It’s also challenging to achieve acceleration control directly [8]. Some bio-intelligent control algorithms can overcome mathematical model problem of complex plants and have better control performances with physiological regulation to achieve multiobjective control [1, 9].Neuroendocrine system (NES) is a major homeostatic system in human body and has some outstanding multiobjective cooperative modulation mechanisms. Being a multiloop feedback mechanism, NES can still regulate the functions of several organs and glands with high self-adaptability and stability, by means of regulating their hormone secretions synchronously [10]. Some researchers have presented several models for modulation mechanism [11], feedback control [12], and hormone release [13] of NES. Based on such mechanisms, some novel artificial neuroendocrine systems (ANES) have been developed and applied to the complex control field. Neal and Timmis [14] proposed the first artificial endocrine system (AES) which includes secretion, regulation and control of hormones. The theory is applied to design a useful emotional mechanism for robot control. Vargas et al. [15] has extended the previous work of literature [14], studied the interactions between the nervous and endocrine systems and provided a comprehensive methodology to design a novel AES for autonomous robot navigation. Córdova and Cañete [16] discussed in conceptual terms the feasibility of designing an ANES in robots and to reflect upon the bionic issues highly associated with complex automatons.To achieve multi-objective cooperative control, some recent work concentrates on how to use multi-loop and multi-objective regulation mechanism of the NES to design some novel control structures and systems. Stear [11] summarized all hormone regulation processes and described a series of control structures. Liu et al. [17] designed a NES-based two-level structure controller, which can not only achieve accurate control but adjust control parameters in real time as well. Ding and Liu et al. [18] developed a bio-inspired decoupling controller from the bi regulation principle of growth hormone in NES. Tang et al. [19] presents an NES-inspired approach for adaptive manufacturing control system. Based on NES, Guo et al. [20] proposes a position-velocity cooperative intelligent controller for motor motion. Compared to conventional control system, these novel control systems always have better simplicity, practicality, stability, and adaptability. These approaches provide some new ideas to multi-objective integrated control field and have good results in simulation. Nevertheless, no experiment has been done on actual plants, especially for multi-objective cooperative control of the position, velocity and acceleration of different parts of the controlled plant.In this paper, a novel multi-loop and multi-objective cooperative intelligent control system (MMCICS) based on regulation mechanism of NES is proposed. Inherited from NES, the MMCICS consists of four subunits: Planning unit regulates position, velocity and acceleration signals based on ultralong loop feedback. Selection unit is a soft switcher to smoothly select dominant motion control signal based on long loop feedback. As short loop feedback, coordination unit is responsible for processing and transmitting coordination signals to several sub-channels in execution unit. The execution unit is an integrity whose sub-channels interact harmoniously and systematically based on ultra-short loop feedback. Each channel has a proposed hormone regulation self-adaptive module (HRSM) which identifies control error and regulates control parameters in real time. The control performance of the proposed MMCICS is verified by an actual 2-DOF redundant parallel manipulator. The experimental results demonstrate that, through regulation mechanism of the MMCICS, the multiobjective integrated control task can be achieved easily while the stability, accuracy, adaptability, and response rate of the plant is improved by proposed HRSM.The main contribution of this paper lies in that it generalizes the characteristics of the NES for regulation, and then reveals the similarity between the NES and a motion control system where the coordination of position, velocity, and acceleration are implemented by the cooperation of different subchannels of the plant. Furthermore, based on the regulation characteristics of the NES, a bioinspired motion control approach is provided, it has been used in MMCICS design. According to our knowledge, this is the first time that the MMCICS based on biological NES is proposed and especially applied to an actual manipulator. The proposed approach is practical and easy to implement, which provides a new efficient method for the intelligent control of complex systems.The remainder of this paper is arranged as follows. In Section2, the regulation mechanism of the NES is described while a corresponding bio-inspired motion control approach is presented. In Section 3, the detailed design of the MMCICS is elaborated including system structure, control algorithms, and parameters tuning methods. The experimental results are given to verify the effectiveness of the proposed control system in Section 4. Finally, the work is summarized in Section 5.


## 2. Regulation Mechanism and Bioinspired Motion Control Approach

### 2.1. Regulation Mechanism of Neuroendocrine System

The NES mainly includes nervous system and endocrine system [15]. The nervous system is primarily responsible for receiving stimuli of environmental change and processing corresponding nerve impulse. The endocrine system can be viewed as a system of glands that works with the nervous system in regulating the activity of internal glands and coordinating the long-range response to external stimuli [21]. One of the most important interactions between them is regulated by means of their hormone secretions.A typical regulation mechanism of the neuroendocrine hormone can be generalized as follows [12, 13, 20, 22]: central nervous system detects the changes in the internal and external environments and transmits the nerve impulse as appropriate response to hypothalamus. Hypothalamus receives the nerve impulses and secretes relevant releasing hormone (RH), which stimulates pituitary to secrete tropic hormone (TH). Under the influence of pituitary’s TH, other glands (such as thyroid, adrenal, gonads, etc.) secrete corresponding hormones which regulate the situation of human physiological balance. There are massive of feedback loops in neuroendocrine system. Four types of typical feedbacks include ultra-short, short, long and ultra-long loop feedbacks [22, 23]. The ultrashort loop feedback means that the hormone released by a certain gland is directly fed back to its source and changes its status. In the short, long and ultralong loop feedbacks, the concentration of corresponding hormone is fed back to the pituitary, hypothalamus and central nervous system, respectively. Through the multiloop feedback mechanism, multihormone control is stable and easy to practice, as shown in Figure 1.Figure 1
Hormone regulation of the NES.


### 2.2. Regulation Characteristics and Bioinspired Motion Control Approach

The regulation characteristics of NES can be summarized as below: (1) the NES has several feedback loops and glands. Each feedback mechanism has its own function and different messages can be transferred among them so that the whole system has a multiobjective regulation mechanism of integrity. (2) Central nervous system is the foremost command center. (3) Hypothalamus is the medium between the nervous system and the endocrine system. (4) Pituitary has the ability to achieve multi-hormone coordinative control. (5) The different glands always have different hormone secretion scopes and different hormone secretion standards. But they have the similar regulation mechanism that can enhance identification and secretion precision within a certain range of stimulus [13].Therefore, corresponding to the motion control system, the central nervous system, the hypothalamus, the pituitary, and glands of NES can be regarded as the planning unit, the selection unit, the coordination unit, and the execution unit, respectively. In this scenario, the planning unit receives input signal and transmits the suitable motion planning signal to the selection unit. The selection unit processes the motion planning signal and chooses the dominant motion control signal. And then, the coordination unit converts dominant motion control signal to various coordination signals according to its performance characteristic. Various sub-channels in the execution unit receive their own coordination signal from the coordination unit and accomplish homologous task. Ultimately, the whole system could be controlled through the combined action of these sub-channels.


## 2.1. Regulation Mechanism of Neuroendocrine System

The NES mainly includes nervous system and endocrine system [15]. The nervous system is primarily responsible for receiving stimuli of environmental change and processing corresponding nerve impulse. The endocrine system can be viewed as a system of glands that works with the nervous system in regulating the activity of internal glands and coordinating the long-range response to external stimuli [21]. One of the most important interactions between them is regulated by means of their hormone secretions.A typical regulation mechanism of the neuroendocrine hormone can be generalized as follows [12, 13, 20, 22]: central nervous system detects the changes in the internal and external environments and transmits the nerve impulse as appropriate response to hypothalamus. Hypothalamus receives the nerve impulses and secretes relevant releasing hormone (RH), which stimulates pituitary to secrete tropic hormone (TH). Under the influence of pituitary’s TH, other glands (such as thyroid, adrenal, gonads, etc.) secrete corresponding hormones which regulate the situation of human physiological balance. There are massive of feedback loops in neuroendocrine system. Four types of typical feedbacks include ultra-short, short, long and ultra-long loop feedbacks [22, 23]. The ultrashort loop feedback means that the hormone released by a certain gland is directly fed back to its source and changes its status. In the short, long and ultralong loop feedbacks, the concentration of corresponding hormone is fed back to the pituitary, hypothalamus and central nervous system, respectively. Through the multiloop feedback mechanism, multihormone control is stable and easy to practice, as shown in Figure 1.Figure 1
Hormone regulation of the NES.


## 2.2. Regulation Characteristics and Bioinspired Motion Control Approach

The regulation characteristics of NES can be summarized as below: (1) the NES has several feedback loops and glands. Each feedback mechanism has its own function and different messages can be transferred among them so that the whole system has a multiobjective regulation mechanism of integrity. (2) Central nervous system is the foremost command center. (3) Hypothalamus is the medium between the nervous system and the endocrine system. (4) Pituitary has the ability to achieve multi-hormone coordinative control. (5) The different glands always have different hormone secretion scopes and different hormone secretion standards. But they have the similar regulation mechanism that can enhance identification and secretion precision within a certain range of stimulus [13].Therefore, corresponding to the motion control system, the central nervous system, the hypothalamus, the pituitary, and glands of NES can be regarded as the planning unit, the selection unit, the coordination unit, and the execution unit, respectively. In this scenario, the planning unit receives input signal and transmits the suitable motion planning signal to the selection unit. The selection unit processes the motion planning signal and chooses the dominant motion control signal. And then, the coordination unit converts dominant motion control signal to various coordination signals according to its performance characteristic. Various sub-channels in the execution unit receive their own coordination signal from the coordination unit and accomplish homologous task. Ultimately, the whole system could be controlled through the combined action of these sub-channels.


## 3. MMCICS Design Inspired from NES

### 3.1. MMCICS Structure Design

According to the bioinspired motion control approach, a novel multi-loop and multi-objective cooperative intelligent control system (MMCICS) is proposed to achieve intelligent coordination of position, velocity, and aceleration implemented by cooperation of several subchannels of plants, as shown in Figure2.Figure 2
The structure of MMCICS.


### 3.2. Units Design of MMCICS

#### 3.2.1. Planning Unit

The planning unit is primarily responsible for receiving and processing input signals of the positionPin(t), the velocity Vin(t), and the acceleration Ain(t), and transmitting the desired position Pout(t), the desired velocity Vout(t), and the breaking factor εbrake signals to the selection unit. The planning algorithm includes the automatic braking process and the cooperative planning process.(1)
Automatic braking process. The position error is defined as(3.1)eP1(t)=Pin(t)-P(t).
When it satisfies
(3.2)|eP1(t)|≤εbrake,
the input velocity signal is changed automatically to
(3.3)Vin(t)=0,
where
(3.4)εbrake=|V2(tbrake)2Ain(tbrake)|
is the braking factor, tbrake is the initial time of the automatic braking process. The actual position signal P(t) and the actual velocity signal V(t) are obtained via ultra-long feedback.(2)
Cooperative planning process. Since acceleration is hardly to be controlled directly, theVin(t) and the Ain(t) are regulated by the cooperative planning process while the Pin(t) is sent to the selection unit directly. Some typical planning methods have good results and have been used in practice for a long time. In order to test the control performance of the MMCICS more clearly, trapezoid curve method has been chosen in this paper. The algorithm can be described as
(3.5)Pout(t)=Pin(t),Vout(t)={Ain(t)⋅(t-tup)+V(tup),t∈TupVin(t),t∈(Tup∪Tdown)cV(tdown)-Ain(t)⋅(t-tdown),t∈Tdown,
where
(3.6)Tup={t∣Ain(t)⋅(t-tup)+Vout(tup)≤Vin(tup)},Tdown={t∣Vout(tdown)-Ain(t)⋅(t-tdown)≥Vin(tdown)},
where, tup and tdown is the initial time when Vin(tup)>Vout(tup) and Vin(tdown)<Vout(tdown), respectively.


#### 3.2.2. Selection Unit

The selection unit is designed as a switcher for the real-time dominant control mode. This unit receives the actual position feedback signal via long-loop feedback mechanism while the dominant motion control signal is transmitted to the coordination unit. Velocity-velocity control mode is on when the actual position is far from desired position while velocity control signal is sent to keep smooth movement. Position-velocity control mode takes over when the actual position is close to the desired position while position control signal is send to achieve accurate position. This rule for automatic switching is described as follows [6, 20]:
(3.7)strategy={velocity-velocity,r>rc,position-velocity,r≤rc,
where r is the distance between actual position and desired position, and rc is a switcher distance which is decided by current state of plant and switching strategy. To guarantee smooth switch, a simple conversion factor Kc is also designed in the selection unit. The control algorithm can be designed as follows:
(3.8)H(t)={Vout(t),|eP2(t)|>εbrake⋅ηswitcheP2(t)⋅Kc,|eP2(t)|≤εbrake⋅ηswitch,
where
(3.9)eP2(t)=Pout(t)-P(t),Kc=|Vout(tswitch)|εbrake⋅ηswitch,
where H(t) is the output of the selection unit, eP2(t) is the error signal between desired and actual position, 0%<ηswitch≤100% is a switching coefficient which decides switching position, Kc is the conversion factor, and tswitch is the initial time of the switching process.


#### 3.2.3. Coordination Unit

The coordination unit is a coordinator which sends cooperative control signals to each sub-channel of the plant. Many methods and mathematic models are suitable for this unit, the velocity Jacobian matrix is chosen in this paper due to the velocity control is our foremost object. In this scenario, all the input signals and output signals are regarded as the velocity signals whether the velocity-velocity control mode or the position-velocity control mode is selected. That output signals can be calculated by(3.10)[C1(t),C2(t),…,Cn(t)]T=J⋅H(t),
where Ci(t) is the ouput signal of the coordination unit to channel i, (i=1,2,…,n) of the execution unit, J is the velocity Jacobian matrix of the plant.


#### 3.2.4. Execution Unit

The execution unit, which includes a number of sub-channels, is the core and key unit of the MMCICS. To keep sub-channels interact harmoniously and systematically, the same control method and control structure have been applied to each channel. As shown in Figure3, each channel has its own independent control subsystem which includes a primary controller, a hormone regulation self-adaptive module (HRSM), and a controlled subpart of plant. There are two ultra-short loop feedbacks. One is that the actual velocity signal is fed back to the primary controller; the other is that the adjusted control parameters are fed back to the HRSM, which can improve the local and global control effectiveness.Figure 3
The structure of sub-channel.Some advanced controllers widely used in industry can be applied as primary controller. The controller can obey PID control algorithm, fuzzy control algorithm [24, 25], H-infinity control algorithm [26, 27], and so forth. Due to their simpledescription, high-dependability, and satisfactory performances, in the MMCICS, the control law of primary controller obeys the conventional PID control algorithm
(3.11)Oi(t)=Kpi0⋅ei(t)+Kii0⋅∫ei(t)dt+Kdi0⋅dei(t)dt,
where
(3.12)ei(t)=Ci(t)-vi(t)
is the error signal between the input signal Ci(t) and the actual velocity vi(t) of the part i,Oi(t) is the output of the primary controller, Kpi0, Kii0, and Kdi0 are the initial PID parameters.The HRSM is designed to improve primary controller self-adaptive performance. The regulation algorithm of HRSM is inspired from hormone regulation mechanism which includes identification and regulation processes.(1)
Identification. In NES, the gland can enhance identification and secretion precision within the working scope. However,  when the stimulate signal beyond the control scope, hormone secretion rate is at its high limit. Similarly, the control errorei(t) in HRSM can be regarded as the stimulate signal, and its identification approach follows the principle of the hormone secretion. Therefore, the absolute value of control error ei(t) is calculated at first and then mapped to the correspond ding regulation scope. Hormone identification error 0≤Ei(t)≤1 is defined as
(3.13)Ei(t)={|ei(t)|eimax⁡-eimin⁡,|ein(t)|<eimax⁡-eimin⁡,1,|ein(t)|≥eimax⁡-eimin⁡,
where eimax and eimin are the high and low limited error of the optimal working scope, respectively.(2)
Regulation. The hormone secretion rate in NES is always nonnegative and monotone, and its secretion regulation mechanism usually follows the Hill functions, the growth curve, and so forth [13, 21]. Based on the Sigmoid function, a hormone regulation factor is designed to regulate primary controller parameter as
(3.14)αij(t)=kij1+(kij-1)e-βij((Ei(t)/ηij)-1),
where j=p,i,d, 0%<ηij≤100% is the critical regulation coefficient, kij≥1 is the high limited regulation coefficient, 0<βij≤10 is the sensitivity regulation coefficient. These three coefficients joint control the function curve’s slope. Where ηij decides the critical point between the up- and down-regulation, as
(3.15)αij(t)<1,Ei(t)<ηij,αij(t)=1,Ei(t)=ηij,αij(t)>1,Ei(t)>ηij.
The kij decides the high limited value. Because if (Ei(t)/ηij)≫1, then e-βij((Ei(t)/ηij)-1)→0 that αij(t)→kij. Meanwhile, it also should be noted that if  kij=1, then αij(t)=1. The βij decides the response rate and has a major impact on the low limited value of αij(t). When βij is bigger, the αij(t) curve changes acutely and the low limit of αij(t) is lower; in contrast, the gentle changes results to higher low limit.Then primary controller parameter can be regulated by its control characteristic. In the PID control algorithm, when the control error is too big, the proportion gainKpi0 should decrease to weaken the control action, thus reduces the overshoot. In contrast, the proportion gain should increase to enhance control precision and eliminate control error quickly [20]. The correcting regulation of the integral coefficient Kii0 and the differential coefficient Kdi0 are similar to that of the proportion gain. Therefore, the parameter regulation algorithm of the PID controller is
(3.16)Kpi(t)=Kpi0/αip(t)Kii(t)=Kii0⋅αii(t)Kdi(t)=Kdi0αid(t).
where, when αip(t)>1, Kpi0 will be reduced; when αip(t)<1, Kpi0 will be increased; when αip(t)=1, Kpi0 will not be changed. Meanwhile, Kii0 and Kdi0 have similar regulation characteristics. The regulation principle of the HRSM satisfies the optimization task and then (3.11) will be changed to optimized control law
(3.17)Oi(t)=Kpi(t)⋅ei(t)+Kii(t)⋅∫ei(t)dt+Kdi(t)⋅dei(t)dt,
where Kpi(t), Kii(t), and Kdi(t) are optimized control parameters.


### 3.3. Parameters Tuning of MMCICS

(1)
Tune the primary controller parameter. First, only take the primary controller into action, and then tune the initial control parametersKpi0, Kii0, and Kdi0 approximately.
(2)
Determine the high and low limited hormone identification error. According to the response characteristics of the experimental results in step (1), determine the high limited erroreimax and low limited error eimin of the optimal working scope.
(3)
Tune the regulation coefficients of the hormone regulator. Take the execution unit into action, according to the response characteristic and overshoot of the experimental results, tune the critical regulation coefficientηij to decide critical working point of the hormone regulator. And then when control error ei(t) is too big, tune the high limited regulation coefficient kij to ensure a stable and faster movement of the plant with little or without overshoot. In contrast, tune the sensitivity regulation coefficient βij to ensure accuracy and stability.
(4)
Determine the switching coefficient. Take the MMCICS into action  and then determine the switching coefficientηswitch to ensure the control strategy switching smoothly.


## 3.1. MMCICS Structure Design

According to the bioinspired motion control approach, a novel multi-loop and multi-objective cooperative intelligent control system (MMCICS) is proposed to achieve intelligent coordination of position, velocity, and aceleration implemented by cooperation of several subchannels of plants, as shown in Figure2.Figure 2
The structure of MMCICS.


## 3.2. Units Design of MMCICS

### 3.2.1. Planning Unit

The planning unit is primarily responsible for receiving and processing input signals of the positionPin(t), the velocity Vin(t), and the acceleration Ain(t), and transmitting the desired position Pout(t), the desired velocity Vout(t), and the breaking factor εbrake signals to the selection unit. The planning algorithm includes the automatic braking process and the cooperative planning process.(1)
Automatic braking process. The position error is defined as(3.1)eP1(t)=Pin(t)-P(t).
When it satisfies
(3.2)|eP1(t)|≤εbrake,
the input velocity signal is changed automatically to
(3.3)Vin(t)=0,
where
(3.4)εbrake=|V2(tbrake)2Ain(tbrake)|
is the braking factor, tbrake is the initial time of the automatic braking process. The actual position signal P(t) and the actual velocity signal V(t) are obtained via ultra-long feedback.(2)
Cooperative planning process. Since acceleration is hardly to be controlled directly, theVin(t) and the Ain(t) are regulated by the cooperative planning process while the Pin(t) is sent to the selection unit directly. Some typical planning methods have good results and have been used in practice for a long time. In order to test the control performance of the MMCICS more clearly, trapezoid curve method has been chosen in this paper. The algorithm can be described as
(3.5)Pout(t)=Pin(t),Vout(t)={Ain(t)⋅(t-tup)+V(tup),t∈TupVin(t),t∈(Tup∪Tdown)cV(tdown)-Ain(t)⋅(t-tdown),t∈Tdown,
where
(3.6)Tup={t∣Ain(t)⋅(t-tup)+Vout(tup)≤Vin(tup)},Tdown={t∣Vout(tdown)-Ain(t)⋅(t-tdown)≥Vin(tdown)},
where, tup and tdown is the initial time when Vin(tup)>Vout(tup) and Vin(tdown)<Vout(tdown), respectively.


### 3.2.2. Selection Unit

The selection unit is designed as a switcher for the real-time dominant control mode. This unit receives the actual position feedback signal via long-loop feedback mechanism while the dominant motion control signal is transmitted to the coordination unit. Velocity-velocity control mode is on when the actual position is far from desired position while velocity control signal is sent to keep smooth movement. Position-velocity control mode takes over when the actual position is close to the desired position while position control signal is send to achieve accurate position. This rule for automatic switching is described as follows [6, 20]:
(3.7)strategy={velocity-velocity,r>rc,position-velocity,r≤rc,
where r is the distance between actual position and desired position, and rc is a switcher distance which is decided by current state of plant and switching strategy. To guarantee smooth switch, a simple conversion factor Kc is also designed in the selection unit. The control algorithm can be designed as follows:
(3.8)H(t)={Vout(t),|eP2(t)|>εbrake⋅ηswitcheP2(t)⋅Kc,|eP2(t)|≤εbrake⋅ηswitch,
where
(3.9)eP2(t)=Pout(t)-P(t),Kc=|Vout(tswitch)|εbrake⋅ηswitch,
where H(t) is the output of the selection unit, eP2(t) is the error signal between desired and actual position, 0%<ηswitch≤100% is a switching coefficient which decides switching position, Kc is the conversion factor, and tswitch is the initial time of the switching process.


### 3.2.3. Coordination Unit

The coordination unit is a coordinator which sends cooperative control signals to each sub-channel of the plant. Many methods and mathematic models are suitable for this unit, the velocity Jacobian matrix is chosen in this paper due to the velocity control is our foremost object. In this scenario, all the input signals and output signals are regarded as the velocity signals whether the velocity-velocity control mode or the position-velocity control mode is selected. That output signals can be calculated by(3.10)[C1(t),C2(t),…,Cn(t)]T=J⋅H(t),
where Ci(t) is the ouput signal of the coordination unit to channel i, (i=1,2,…,n) of the execution unit, J is the velocity Jacobian matrix of the plant.


### 3.2.4. Execution Unit

The execution unit, which includes a number of sub-channels, is the core and key unit of the MMCICS. To keep sub-channels interact harmoniously and systematically, the same control method and control structure have been applied to each channel. As shown in Figure3, each channel has its own independent control subsystem which includes a primary controller, a hormone regulation self-adaptive module (HRSM), and a controlled subpart of plant. There are two ultra-short loop feedbacks. One is that the actual velocity signal is fed back to the primary controller; the other is that the adjusted control parameters are fed back to the HRSM, which can improve the local and global control effectiveness.Figure 3
The structure of sub-channel.Some advanced controllers widely used in industry can be applied as primary controller. The controller can obey PID control algorithm, fuzzy control algorithm [24, 25], H-infinity control algorithm [26, 27], and so forth. Due to their simpledescription, high-dependability, and satisfactory performances, in the MMCICS, the control law of primary controller obeys the conventional PID control algorithm
(3.11)Oi(t)=Kpi0⋅ei(t)+Kii0⋅∫ei(t)dt+Kdi0⋅dei(t)dt,
where
(3.12)ei(t)=Ci(t)-vi(t)
is the error signal between the input signal Ci(t) and the actual velocity vi(t) of the part i,Oi(t) is the output of the primary controller, Kpi0, Kii0, and Kdi0 are the initial PID parameters.The HRSM is designed to improve primary controller self-adaptive performance. The regulation algorithm of HRSM is inspired from hormone regulation mechanism which includes identification and regulation processes.(1)
Identification. In NES, the gland can enhance identification and secretion precision within the working scope. However,  when the stimulate signal beyond the control scope, hormone secretion rate is at its high limit. Similarly, the control errorei(t) in HRSM can be regarded as the stimulate signal, and its identification approach follows the principle of the hormone secretion. Therefore, the absolute value of control error ei(t) is calculated at first and then mapped to the correspond ding regulation scope. Hormone identification error 0≤Ei(t)≤1 is defined as
(3.13)Ei(t)={|ei(t)|eimax⁡-eimin⁡,|ein(t)|<eimax⁡-eimin⁡,1,|ein(t)|≥eimax⁡-eimin⁡,
where eimax and eimin are the high and low limited error of the optimal working scope, respectively.(2)
Regulation. The hormone secretion rate in NES is always nonnegative and monotone, and its secretion regulation mechanism usually follows the Hill functions, the growth curve, and so forth [13, 21]. Based on the Sigmoid function, a hormone regulation factor is designed to regulate primary controller parameter as
(3.14)αij(t)=kij1+(kij-1)e-βij((Ei(t)/ηij)-1),
where j=p,i,d, 0%<ηij≤100% is the critical regulation coefficient, kij≥1 is the high limited regulation coefficient, 0<βij≤10 is the sensitivity regulation coefficient. These three coefficients joint control the function curve’s slope. Where ηij decides the critical point between the up- and down-regulation, as
(3.15)αij(t)<1,Ei(t)<ηij,αij(t)=1,Ei(t)=ηij,αij(t)>1,Ei(t)>ηij.
The kij decides the high limited value. Because if (Ei(t)/ηij)≫1, then e-βij((Ei(t)/ηij)-1)→0 that αij(t)→kij. Meanwhile, it also should be noted that if  kij=1, then αij(t)=1. The βij decides the response rate and has a major impact on the low limited value of αij(t). When βij is bigger, the αij(t) curve changes acutely and the low limit of αij(t) is lower; in contrast, the gentle changes results to higher low limit.Then primary controller parameter can be regulated by its control characteristic. In the PID control algorithm, when the control error is too big, the proportion gainKpi0 should decrease to weaken the control action, thus reduces the overshoot. In contrast, the proportion gain should increase to enhance control precision and eliminate control error quickly [20]. The correcting regulation of the integral coefficient Kii0 and the differential coefficient Kdi0 are similar to that of the proportion gain. Therefore, the parameter regulation algorithm of the PID controller is
(3.16)Kpi(t)=Kpi0/αip(t)Kii(t)=Kii0⋅αii(t)Kdi(t)=Kdi0αid(t).
where, when αip(t)>1, Kpi0 will be reduced; when αip(t)<1, Kpi0 will be increased; when αip(t)=1, Kpi0 will not be changed. Meanwhile, Kii0 and Kdi0 have similar regulation characteristics. The regulation principle of the HRSM satisfies the optimization task and then (3.11) will be changed to optimized control law
(3.17)Oi(t)=Kpi(t)⋅ei(t)+Kii(t)⋅∫ei(t)dt+Kdi(t)⋅dei(t)dt,
where Kpi(t), Kii(t), and Kdi(t) are optimized control parameters.


## 3.2.1. Planning Unit

The planning unit is primarily responsible for receiving and processing input signals of the positionPin(t), the velocity Vin(t), and the acceleration Ain(t), and transmitting the desired position Pout(t), the desired velocity Vout(t), and the breaking factor εbrake signals to the selection unit. The planning algorithm includes the automatic braking process and the cooperative planning process.(1)
Automatic braking process. The position error is defined as(3.1)eP1(t)=Pin(t)-P(t).
When it satisfies
(3.2)|eP1(t)|≤εbrake,
the input velocity signal is changed automatically to
(3.3)Vin(t)=0,
where
(3.4)εbrake=|V2(tbrake)2Ain(tbrake)|
is the braking factor, tbrake is the initial time of the automatic braking process. The actual position signal P(t) and the actual velocity signal V(t) are obtained via ultra-long feedback.(2)
Cooperative planning process. Since acceleration is hardly to be controlled directly, theVin(t) and the Ain(t) are regulated by the cooperative planning process while the Pin(t) is sent to the selection unit directly. Some typical planning methods have good results and have been used in practice for a long time. In order to test the control performance of the MMCICS more clearly, trapezoid curve method has been chosen in this paper. The algorithm can be described as
(3.5)Pout(t)=Pin(t),Vout(t)={Ain(t)⋅(t-tup)+V(tup),t∈TupVin(t),t∈(Tup∪Tdown)cV(tdown)-Ain(t)⋅(t-tdown),t∈Tdown,
where
(3.6)Tup={t∣Ain(t)⋅(t-tup)+Vout(tup)≤Vin(tup)},Tdown={t∣Vout(tdown)-Ain(t)⋅(t-tdown)≥Vin(tdown)},
where, tup and tdown is the initial time when Vin(tup)>Vout(tup) and Vin(tdown)<Vout(tdown), respectively.


## 3.2.2. Selection Unit

The selection unit is designed as a switcher for the real-time dominant control mode. This unit receives the actual position feedback signal via long-loop feedback mechanism while the dominant motion control signal is transmitted to the coordination unit. Velocity-velocity control mode is on when the actual position is far from desired position while velocity control signal is sent to keep smooth movement. Position-velocity control mode takes over when the actual position is close to the desired position while position control signal is send to achieve accurate position. This rule for automatic switching is described as follows [6, 20]:
(3.7)strategy={velocity-velocity,r>rc,position-velocity,r≤rc,
where r is the distance between actual position and desired position, and rc is a switcher distance which is decided by current state of plant and switching strategy. To guarantee smooth switch, a simple conversion factor Kc is also designed in the selection unit. The control algorithm can be designed as follows:
(3.8)H(t)={Vout(t),|eP2(t)|>εbrake⋅ηswitcheP2(t)⋅Kc,|eP2(t)|≤εbrake⋅ηswitch,
where
(3.9)eP2(t)=Pout(t)-P(t),Kc=|Vout(tswitch)|εbrake⋅ηswitch,
where H(t) is the output of the selection unit, eP2(t) is the error signal between desired and actual position, 0%<ηswitch≤100% is a switching coefficient which decides switching position, Kc is the conversion factor, and tswitch is the initial time of the switching process.


## 3.2.3. Coordination Unit

The coordination unit is a coordinator which sends cooperative control signals to each sub-channel of the plant. Many methods and mathematic models are suitable for this unit, the velocity Jacobian matrix is chosen in this paper due to the velocity control is our foremost object. In this scenario, all the input signals and output signals are regarded as the velocity signals whether the velocity-velocity control mode or the position-velocity control mode is selected. That output signals can be calculated by(3.10)[C1(t),C2(t),…,Cn(t)]T=J⋅H(t),
where Ci(t) is the ouput signal of the coordination unit to channel i, (i=1,2,…,n) of the execution unit, J is the velocity Jacobian matrix of the plant.


## 3.2.4. Execution Unit

The execution unit, which includes a number of sub-channels, is the core and key unit of the MMCICS. To keep sub-channels interact harmoniously and systematically, the same control method and control structure have been applied to each channel. As shown in Figure3, each channel has its own independent control subsystem which includes a primary controller, a hormone regulation self-adaptive module (HRSM), and a controlled subpart of plant. There are two ultra-short loop feedbacks. One is that the actual velocity signal is fed back to the primary controller; the other is that the adjusted control parameters are fed back to the HRSM, which can improve the local and global control effectiveness.Figure 3
The structure of sub-channel.Some advanced controllers widely used in industry can be applied as primary controller. The controller can obey PID control algorithm, fuzzy control algorithm [24, 25], H-infinity control algorithm [26, 27], and so forth. Due to their simpledescription, high-dependability, and satisfactory performances, in the MMCICS, the control law of primary controller obeys the conventional PID control algorithm
(3.11)Oi(t)=Kpi0⋅ei(t)+Kii0⋅∫ei(t)dt+Kdi0⋅dei(t)dt,
where
(3.12)ei(t)=Ci(t)-vi(t)
is the error signal between the input signal Ci(t) and the actual velocity vi(t) of the part i,Oi(t) is the output of the primary controller, Kpi0, Kii0, and Kdi0 are the initial PID parameters.The HRSM is designed to improve primary controller self-adaptive performance. The regulation algorithm of HRSM is inspired from hormone regulation mechanism which includes identification and regulation processes.(1)
Identification. In NES, the gland can enhance identification and secretion precision within the working scope. However,  when the stimulate signal beyond the control scope, hormone secretion rate is at its high limit. Similarly, the control errorei(t) in HRSM can be regarded as the stimulate signal, and its identification approach follows the principle of the hormone secretion. Therefore, the absolute value of control error ei(t) is calculated at first and then mapped to the correspond ding regulation scope. Hormone identification error 0≤Ei(t)≤1 is defined as
(3.13)Ei(t)={|ei(t)|eimax⁡-eimin⁡,|ein(t)|<eimax⁡-eimin⁡,1,|ein(t)|≥eimax⁡-eimin⁡,
where eimax and eimin are the high and low limited error of the optimal working scope, respectively.(2)
Regulation. The hormone secretion rate in NES is always nonnegative and monotone, and its secretion regulation mechanism usually follows the Hill functions, the growth curve, and so forth [13, 21]. Based on the Sigmoid function, a hormone regulation factor is designed to regulate primary controller parameter as
(3.14)αij(t)=kij1+(kij-1)e-βij((Ei(t)/ηij)-1),
where j=p,i,d, 0%<ηij≤100% is the critical regulation coefficient, kij≥1 is the high limited regulation coefficient, 0<βij≤10 is the sensitivity regulation coefficient. These three coefficients joint control the function curve’s slope. Where ηij decides the critical point between the up- and down-regulation, as
(3.15)αij(t)<1,Ei(t)<ηij,αij(t)=1,Ei(t)=ηij,αij(t)>1,Ei(t)>ηij.
The kij decides the high limited value. Because if (Ei(t)/ηij)≫1, then e-βij((Ei(t)/ηij)-1)→0 that αij(t)→kij. Meanwhile, it also should be noted that if  kij=1, then αij(t)=1. The βij decides the response rate and has a major impact on the low limited value of αij(t). When βij is bigger, the αij(t) curve changes acutely and the low limit of αij(t) is lower; in contrast, the gentle changes results to higher low limit.Then primary controller parameter can be regulated by its control characteristic. In the PID control algorithm, when the control error is too big, the proportion gainKpi0 should decrease to weaken the control action, thus reduces the overshoot. In contrast, the proportion gain should increase to enhance control precision and eliminate control error quickly [20]. The correcting regulation of the integral coefficient Kii0 and the differential coefficient Kdi0 are similar to that of the proportion gain. Therefore, the parameter regulation algorithm of the PID controller is
(3.16)Kpi(t)=Kpi0/αip(t)Kii(t)=Kii0⋅αii(t)Kdi(t)=Kdi0αid(t).
where, when αip(t)>1, Kpi0 will be reduced; when αip(t)<1, Kpi0 will be increased; when αip(t)=1, Kpi0 will not be changed. Meanwhile, Kii0 and Kdi0 have similar regulation characteristics. The regulation principle of the HRSM satisfies the optimization task and then (3.11) will be changed to optimized control law
(3.17)Oi(t)=Kpi(t)⋅ei(t)+Kii(t)⋅∫ei(t)dt+Kdi(t)⋅dei(t)dt,
where Kpi(t), Kii(t), and Kdi(t) are optimized control parameters.


## 3.3. Parameters Tuning of MMCICS

(1)
Tune the primary controller parameter. First, only take the primary controller into action, and then tune the initial control parametersKpi0, Kii0, and Kdi0 approximately.
(2)
Determine the high and low limited hormone identification error. According to the response characteristics of the experimental results in step (1), determine the high limited erroreimax and low limited error eimin of the optimal working scope.
(3)
Tune the regulation coefficients of the hormone regulator. Take the execution unit into action, according to the response characteristic and overshoot of the experimental results, tune the critical regulation coefficientηij to decide critical working point of the hormone regulator. And then when control error ei(t) is too big, tune the high limited regulation coefficient kij to ensure a stable and faster movement of the plant with little or without overshoot. In contrast, tune the sensitivity regulation coefficient βij to ensure accuracy and stability.
(4)
Determine the switching coefficient. Take the MMCICS into action  and then determine the switching coefficientηswitch to ensure the control strategy switching smoothly.


## 4. Experimental Results and Analysis

Some typical experimental results are provided in this section to explore two main experiments of proposed MMCICS. Firstly, the control results with and without HRSM are compared to find out whether HRSM yields better in subchannel experiment. Next more comprehensive experiments are performed to verify multiobject cooperative control performance of the MMCICS, and whether HRSM has better global control effect.As shown in Figure4, a 2-DOF redundant parallel manipulator (Googol Tech Ltd.’s GPM2002) [28, 29] is selected as the experiment platform due to its complex redundancy structure and multi-channel inputs. Three bases of the manipulator are equipped with three AC servo motors with harmonic gear drives. The coordinates of three bases are A1(0, 250), A2(433, 0), and A3(433, 500), and all the links have the same length l=244. The unit of coordinates and length is millimeter. Active joint angles are qa1, qa2 and qa3, and passive joint angles are qb1, qb2 and qb3. Position signals of the motors are measured with the absolute optical electrical encoders, and input voltage signals are controlled by a motion control board. All algorithms are implemented with Matlab/Simulink environment on an industrial controlling computer with a 2.8 GHz processor and 1024 MB memory. The real-time implementation is executed with the Real Time Workshop (RTW) of Matlab, and sampling period is 5 ms.Figure 4
The 2-DOF redundant parallel manipulator.Firstly, to verify the effectiveness of the proposed HRSM in the execution unit, we only take active joint 1 (baseA1) without loads and links into action. The control performance of the conventional PID controller and the PID controller with HRSM (HRSM-PID) are compared under the six different velocities of the servo motor 1, namely the motor of base A1. To make the contrast effect more clearly, the conventional PID parameters are designed as the same as the initial PID parameters in HRSM-PID controller, as shown in Table 1.Table 1
Parameter set.
Initial PID
Error factors
Hormone regulation factors
K
p
i
0
=
0.03
e
i
max
=
1 r/s
k
i
p
=
5, βip=1.2, ηip=20%
K
i
i
0
=
0.004
e
i
min
=
-
1 r/s
k
i
i
=
5,  βii=1.2, ηii=15%
K
d
i
0
=
0.005
k
i
d
=
3, βid=1,  ηid=10%Motor in sub-channel has different dynamic characteristics at different velocities but has similar results in the same parameter sets. Multiple experiments have the similar results, and a typical result is as shown in Figure5(a), when motor is running at low velocities, the steady-state errors are obvious due to load influence. The HRSM-PID controller achieves better stabilities, higher accuracies, slightly faster dynamic responses, and lower or no overshoots, compared with the conventional PID controller. Figure 5(b) shows that when running at high velocities, the motor has better motion performance and spends more times to achieve higher velocity. HRSM-PID controller achieves significantly faster dynamic responses compared with the PID controller. Figure 5(c) shows a typical output control signal O1(t) when input velocity step is 5. As the expected, when the error is too big, the HRSM decreases the output control signals to reduce the overshoot. In contrast, the output control signals are increased to enhance control precision and eliminate control error quickly. With such strong self- adaptability, the HRSM improves the dynamic performances. The detailed lower quartile, median, upper quartile, average, and variance of the 10 time’s results are shown in Table 2. Where, Vd is the desired velocity, ts is the settling time, σ is the overshoot, and |ess| is the absolutely value of steady-state error. The sub-channel experimental results show that based on hormone regulation mechanism, the HRSM owns strong self-adaptability that improves the response, accuracy, and stability of the subchannel.Table 2
Performance evaluation for subchannel experiment.
V
d (r/s)
1
5
10
20
40
50
HRSMPID
PID
HRSMPID
PID
HRSMPID
PID
HRSMPID
PID
HRSMPID
PID
Lower quartile
0.0475
0.1575
0.18
0.2275
0.25
0.32
0.3
0.3575
0.32
0.4275
0.3575
0.48
Median
0.06
0.175
0.2
0.245
0.265
0.33
0.32
0.375
0.33
0.43
0.38
0.5
t
s(s)
Upper quartile
0.0725
0.185
0.22
0.2525
0.285
0.3525
0.34
0.385
0.3525
0.4525
0.4
0.54
Average
0.061
0.175
0.201
0.246
0.272
0.335
0.324
0.37
0.336
0.439
0.379
0.506
Variance
0.000309
0.000445
0.000369
0.000484
0.000496
0.000445
0.000384
0.00052
0.000524
0.000589
0.000609
0.000964
Lower quartile
0.06
0.095
0.0375
0.1675
0.04
0.16
0.02
0.0575
0.0175
0.04
0.01
0.03
Median
0.08
0.12
0.04
0.185
0.055
0.19
0.04
0.075
0.02
0.045
0.02
0.04
σ(r/s)
Upper quartile
0.1
0.14
0.06
0.225
0.0775
0.22
0.05
0.1
0.04
0.06
0.025
0.0525
Average
0.079
0.12
0.049
0.195
0.11
0.192
0.038
0.077
0.024
0.049
0.021
0.045
Variance
0.000449
0.00064
0.000449
0.001165
0.027
0.000736
0.000176
0.000501
0.000124
0.000249
0.000109
0.000225
Lower quartile
0.05
0.08
0.08
0.215
0.06
0.195
0.04
0.08
0.03
0.06
0.02
0.0475
Median
0.06
0.105
0.1
0.245
0.08
0.22
0.05
0.1
0.04
0.065
0.035
0.06
|
ess
| (r/s)
Upper quartile
0.085
0.12
0.1275
0.2525
0.085
0.2575
0.0625
0.12
0.045
0.08
0.04
0.08
Average
0.068
0.102
0.103
0.237
0.078
0.225
0.052
0.099
0.04
0.066
0.033
0.063
Variance
0.000576
0.000656
0.000921
0.001021
0.000376
0.001705
0.000256
0.000609
0.00014
0.000144
0.000181
0.000321Contrast effect of the velocity control. (a) Low velocity control, (b) high velocity control. (c) Output control signal.
(a)
(b)
(c)To verify the multiobject cooperative control performance of the MMCICS, the end-effector of the redundant parallel manipulator is viewed as a controlled plant, and three active joints are viewed as three subchannels. The velocity Jacobian matrix between the end-effector and three active joints is(4.1)J=1l(cos⁡qb1sin⁡(qb1-qa1)sin⁡qb1sin⁡(qb1-qa1)cos⁡qb2sin⁡(qb2-qa2)sin⁡qb2sin⁡(qb2-qa2)cos⁡qb3sin⁡(qb3-qa3)sin⁡qb3sin⁡(qb3-qa3)).Due to the complex mechanism structure of the parallel manipulator with actuation redundancy, it is a typical nonlinear system and difficult to get the accurate dynamic and friction model [28, 29]. Although the manipulator has different dynamic characteristics in different positions, velocities, and accelerations, the proposed MMCICS can overcome accurate mathematical model problem. To verify the control performance of the MMCICS more thoroughly and whether the HRSM also achieves better control effectiveness in the proposed MMCICS, many different experiments were tested and have similar results. A representative contrast experimental result is shown in Figure 6, where the MMCICS without HRSM is chosen as contrast control system (CCS). The experiments are implemented with the same input signals and control parameters. The starting position, input goal position, and input acceleration are [216.5,250]T, [316.5,350]T, and [1500,1500]T, respectively. The input velocity signal is
(4.2)Vin(t)=[0,0]T,t≥0sVin(t)=[100,100]T,t≥0.1sVin(t)=[300,300]T,t≥3.5s.
The switching coefficient is ηswitch=[20%,20%]T in selection unit, and control parameters in channel 1, 2, and 3 are the same as in Table 1. Similarly, the parameters of CCS are the same as MMCICS.Multichannel control experimental results. (a)X-direction velocity. (b) Y-direction velocity. (c) X-direction position. (d) Y-direction position.
(a)
(b)
(c)
(d)As shown in Figures6(a) and 6(b), MMCICS achieves a faster response, better stability, and higher accuracy of velocity control compared with CCS. Especially, when the manipulator is running at high velocities, it is hard to achieve object velocity using CCS, due to complex plan structure and big load. However, MMCICS still maintains high performance as low velocity process. From the velocity response during the ascent, it’s easy to find that the MMCICS has more stable acceleration response than CCS does. That means, based on the HRSM, the cooperative planning algorithm in the planning unit can be implemented easier for acceleration control. Moreover, during control strategy switching, CCS always has a significant negative overshoot of the velocity in braking process. In contrast, MMCICS can stop quickly with little or no negative overshoot due to its strong adaptability. Compared with Figures 6(a) and 6(b), we can find that, due to uneven distribution of loads, the CCS performance in Y-direction is worse than X-direction. However, MMCICS can overcome this problem, since its local self-adaptability improves the global self-adaptability.Figures6(c) and 6(d) shows that due to faster, more stable, and accurate velocity response, the MMCICS can achieve better position accuracy compared with the CCS. In the braking process, because of its better adaptability when control strategy is switched from the velocity-velocity control to the position-velocity control, the MMCICS has a faster position response, which makes position stable with lower overshoot or no overshoot.Some compare results of the 10 time’s average absolute values are shown in Table3. The experimental results show that, with the planning algorithm in the planning unit, the soft switching algorithm in the selection unit, and the velocity cooperative control in the coordination unit, both MMCICS and CCS take advantages of position control and velocity control, and achieve cooperative control for position, velocity and acceleration. Particularly, with strong self-adaptability, faster response, and better stability of HRSM, control potentials of the MMCICS are exploited more thoroughly. The MMCICS achieves multi-objective cooperative intelligent control with higher performance even at high velocities and accelerations, for a nonlinear multi-input complex plant without accurate dynamics model.Table 3
Performance evaluation for comprehensive experiment.
Positon
Veloctiy
Final error
Settling time
Error
Braking overshoot
V
in
(
t
)
=
[
100
,
100
]
T
V
in
(
t
)
=
[
300
,
300
]
T
MMCICS
[
0.02
,
0.05
]
T
0.82 s
[
6.62
,
7.50
]
T
[
7.22
,
5.64
]
T
[
2.41
,
3.13
]
T
CCS
[
0.05
,
0.22
]
T
1.07 s
[
13.82
,
18.20
]
T
[
32.2
,
52.06
]
T
[
10.25
,
25.20
]
T


## 5. Conclusions

This work presents a bioinspired cooperative intelligent control system for position, velocity, and acceleration multi-objective integrated control of a parallel plant. The similarity between the NES and motion control system revealed, and a bio-inspired motion control approach is proposed. Under the context of such approach, the MMCICS with system structure, algorithm, and steps in parameter tuning is proposed to achieve multiobjective control. The experiments are carried out with a 2-DOF redundant parallel manipulator where the feasibility of the new control system is demonstrated. The contrast effect shows that the stability, accuracy, adaptability, response rate of the proposed MMCICS is superior to those of the conventional controllers. According to our knowledge, this is the first time that NES-based MMCICS and HRSM are proposed and used for an actual parallel manipulator. The proposed MMCICS can be implemented easily and provides a new and efficient method for multiobjective integrated control of complex multichannel systems. In future works, force and torque control will be considered to establish a more complete multi-objective control system. More rigorous and advanced algorithm and proof are required instead of the PID controller. Besides, parameter optimization, dynamics, and stability analysis can be conducted on MMCICS.



---

*Source: 467402-2012-09-12.xml*