Nonlinear System Stability and Behavioral Analysis for Effective Implementation of Artificial Lower Limb

System performance and efficiency depends on the stability criteria. The lower limb prosthetic model design requires some prerequisites such as hardware design functionality and compatibility of the building block materials. Effective implementation of mathematical model simulation symmetry towards the achievement of hardware design is the focus of the present work. Different postures of lower limb have been considered in this paper to be analyzed for artificial system design of lower limb movement. The generated polynomial equations of the sitting and standing positions of the normal limb are represented with overall system transfer function. The behavioral analysis of the lower limb model shows the nonlinear nature. The Euler-Lagrange method is utilized to describe the nonlinearity in the field of forward dynamics of the artificial system. The stability factor through phase portrait analysis is checked with respect to nonlinear system characteristics of the lower limb. The asymptotic stability has been achieved utilizing the most applicable Lyapunov method for nonlinear systems. The stability checking of the proposed artificial lower extremity is the newer approach needed to take decisions on output implementation in the system design.


Introduction
Generally, nature is nonlinear as the responses of the physical systems in the applied fields also show nonlinear behavior. Some electromechanical processes show nonlinear behavior. Nonlinear controller can compensate the unwanted effects on the system. An important research field for nonlinear control systems is to design controllers to deal with uncertainty, mainly due to the unavailability of parametric information of the models and external disturbances. Hard nonlinearities such as dead-zones, hysteresis and saturation do not permit linear approximation of real-world systems. The method of approximation deals with the effects of nonlinearity. This can cause inaccuracy in the system. The effects of unwanted nonlinearities in the system should be analyzed to improve the dynamic performance of any real-time system development.
The applicability of forward dynamics [1] in the nonlinear systems such as lower limb models has a great impact on motion analysis. The application of Euler-Lagrange [2,3] method has been presented in the field of robotics for the movement analysis. The mathematical explanation of the systems is the newer approach using the Euler-Lagrange method to show the artificial motion [4]. There are many interesting applications of system performance analysis through mathematical model simulation in nonlinear domains, such as walking robots, Degree of Freedom (DoF) manipulator [5] working as an artificial limb, etc. The analysis and design of nonlinear control system stability has been discussed to The Lagrangian approach delivers the immediate reduction in the characteristic size to solve many equations. This describes the equation of motion. The Lagrange formulation is used to avoid the acceleration finding as shown in Figure 2. The Lagrangian equation can be expressed as the difference between the kinetic energy and the potential energy as shown in the Equation (1).
where, ∅ = time derivative of partial derivative of L with respect to ∅ ∅ = partial derivative of L with respect to ∅ τ = power produced by the joints related to ∅, ∅ .
In the form of vector components, it can be written as shown in Equation (3).

∅
The Lagrangian approach delivers the immediate reduction in the characteristic size to solve many equations. This describes the equation of motion. The Lagrange formulation is used to avoid the acceleration finding as shown in Figure 2. The Lagrangian equation can be expressed as the difference between the kinetic energy and the potential energy as shown in the Equation (1).
where, Kinetic energy = K ∅, . ∅ and Potential energy = P(∅). The Lagrangian approach delivers the immediate reduction in the characteristic size to solve many equations. This describes the equation of motion. The Lagrange formulation is used to avoid the acceleration finding as shown in Figure 2. The Lagrangian equation can be expressed as the difference between the kinetic energy and the potential energy as shown in the Equation (1).
where, Kinetic energy = ∅, ∅ and Potential energy = (∅). Euler-Lagrange equations are suitably obtained for the nonlinear behavioral analysis of the artificial movement of the lower limbs.
where, ∅ = time derivative of partial derivative of L with respect to ∅ ∅ = partial derivative of L with respect to ∅ τ = power produced by the joints related to ∅, ∅ .
In the form of vector components, it can be written as shown in Equation (3). Euler-Lagrange equations are suitably obtained for the nonlinear behavioral analysis of the artificial movement of the lower limbs.
where, In the form of vector components, it can be written as shown in Equation (3).
The joint1 refers to the ankle joint and the joint2 refers to the knee joint in Figure 3. From the above diagram, as shown in the Figure 3, Equation (4) can be written for the position of joint2 (a 2 , b 2 ).
The joint1 refers to the ankle joint and the joint2 refers to the knee joint in Figure 3. From the above diagram, as shown in the Figure 3, Equation (4) can be written for the position of joint2 ( , ).
Similarly, it can be written for the position of end-effector ( , ) as shown in the Equation (6) = cos ∅ + cos(∅ + ∅ ) sin ∅ + sin(∅ + ∅ ) This equation is for the position and velocity of the model as shown in Equation (7).
where, = length of the limb from joint2 to end-effector. Now, the Kinetic Energy for link1 and link2 are shown in the Equations (8) and (9) respectively.
where, = mass of the joint2, = velocity of the joint2. The derivative of (a 1 , b 1 ) is expressed in Equation (5).
Similarly, it can be written for the position of end-effector (a e , b e ) as shown in the Equation (6) This equation is for the position and velocity of the model as shown in Equation (7).
where, x 2 = length of the limb from joint2 to end-effector. Now, the Kinetic Energy for link1 and link2 are shown in the Equations (8) and (9) respectively.
One component of Lagrangian is L comp as shown in the Equation (14). Therefore, The effect of L comp on the 2nd joint torque is as given below in Equation (15).
Now torque for first joint is τ 1 ..
Now torque for second joint is τ 2 ..
Some terms are not dependent on the joint acceleration and velocity. So, the torque can be expressed as shown in Equation (19). where,
The overall acceleration depends not only on the joint accelerations but also on the products of the joint velocities. Joint coordinates are non-inertial. These velocity product terms appear because joint coordinates are non-inertial.

Stability Analysis of the Non-Linear System
The dynamic behavior of the artificial lower extremity model shows nonlinearity. The postures such as sitting and standing conditions of the lower part of human body exhibits the nonlinear characteristics. The control system [31][32][33] involves the mathematical relationship related to the nonlinear properties such as equilibrium point, saturation, dead zone, hysteresis, describing function, etc. Nonlinear systems may not produce a specific solution for all times but gives a decision maker feature [34,35] of the analyzed system. Through phase portrait analysis and the Lyapunov method, the system stability [36,37] of the designed system definition is feasible. Thus, system design method for nonlinear control systems are application-specific [38]. The describing function is the well-established process to analyze the nonlinear function [39] with an approximation, and to determine the stability of a nonlinear system. It is basically the approximated extended version of frequency domain analysis. For a higher degree of nonlinearity, the describing function is used but no information about time domain such as steady state error is given by this. Phase plane analysis can be used, but the order of the system should be less than second order as it is a graphical system. Sometimes it is called a sinusoidal describing function.

Workflow of the System Stability Analysis
The workflow of the analysis for lower limb system [40,41] stability with respect to linear and nonlinear systems is shown in Figures 4-6 The overall acceleration depends not only on the joint accelerations but also on the products of the joint velocities. Joint coordinates are non-inertial. These velocity product terms appear because joint coordinates are non-inertial.

Stability Analysis of the Non-Linear System
The dynamic behavior of the artificial lower extremity model shows nonlinearity. The postures such as sitting and standing conditions of the lower part of human body exhibits the nonlinear characteristics. The control system [31][32][33] involves the mathematical relationship related to the nonlinear properties such as equilibrium point, saturation, dead zone, hysteresis, describing function, etc. Nonlinear systems may not produce a specific solution for all times but gives a decision maker feature [34,35] of the analyzed system. Through phase portrait analysis and the Lyapunov method, the system stability [36,37] of the designed system definition is feasible. Thus, system design method for nonlinear control systems are application-specific [38]. The describing function is the wellestablished process to analyze the nonlinear function [39] with an approximation, and to determine the stability of a nonlinear system. It is basically the approximated extended version of frequency domain analysis. For a higher degree of nonlinearity, the describing function is used but no information about time domain such as steady state error is given by this. Phase plane analysis can be used, but the order of the system should be less than second order as it is a graphical system. Sometimes it is called a sinusoidal describing function.

Workflow of the System Stability Analysis
The workflow of the analysis for lower limb system [40,41] stability with respect to linear and nonlinear systems is shown in Figure 4, Figures 5 and 6, respectively.    The polynomial equations generated from two different postures of the lower limb, such as sitting and standing, are considered as the primary conditions of the artificial system stability analysis. The transfer function is generated from the polynomials to achieve the mathematical model [42] of the mentioned system. The closed loop transfer function is generated with unity feedback to gather information about the controller, performance and robustness parameters of the system, which shows stability as shown in Figure 4.
The characteristics equation is generated from the polynomials of the system at two different positions of the lower limb. The roots are determined using a computational technique. In the nonlinear domain, using Eigen vectors and the generated matrix, the phase portrait analysis has been incorporated. The system shows asymptotic stability in the nonlinear domain for both the roots of the characteristics equation as shown in Figure 5. For better and stronger decision making on the stability of the system, the Lyapunov stability [43] analysis method has been performed with the equilibrium point detection. This method also shows asymptotic stability where the system is stable and roots are in a complex domain as shown in Figure 6.

Results and Discussions
In this present research work, a 3 DoF lower limb has been analyzed for the artificial electromechanical system development. The conditions of movement of the human body such as sitting and standing are the main concern in this discussion. Non applicability of the superposition and homogeneity property gives a greater emphasis on the analysis and design of the control systems in the nonlinear domain [15]. The controllability and observability [44] cannot be determined simply based on rank tests. The concept of system stability [45] and checking using Phase portrait and Lyapunov stability [17] has been incorporated to make it a complete evaluation. Polynomials for the two positions of the lower limb (sitting and standing) as shown in Figures 7 and 8 are presented below in Equations (23) and (24), respectively [5]:  The polynomial equations generated from two different postures of the lower limb, such as sitting and standing, are considered as the primary conditions of the artificial system stability analysis. The transfer function is generated from the polynomials to achieve the mathematical model [42] of the mentioned system. The closed loop transfer function is generated with unity feedback to gather information about the controller, performance and robustness parameters of the system, which shows stability as shown in Figure 4.
The characteristics equation is generated from the polynomials of the system at two different positions of the lower limb. The roots are determined using a computational technique. In the nonlinear domain, using Eigen vectors and the generated matrix, the phase portrait analysis has been incorporated. The system shows asymptotic stability in the nonlinear domain for both the roots of the characteristics equation as shown in Figure 5. For better and stronger decision making on the stability of the system, the Lyapunov stability [43] analysis method has been performed with the equilibrium point detection. This method also shows asymptotic stability where the system is stable and roots are in a complex domain as shown in Figure 6.

Results and Discussions
In this present research work, a 3 DoF lower limb has been analyzed for the artificial electromechanical system development. The conditions of movement of the human body such as sitting and standing are the main concern in this discussion. Non applicability of the superposition and homogeneity property gives a greater emphasis on the analysis and design of the control systems in the nonlinear domain [15]. The controllability and observability [44] cannot be determined simply based on rank tests. The concept of system stability [45] and checking using Phase portrait and Lyapunov stability [17] has been incorporated to make it a complete evaluation. Polynomials for the two positions of the lower limb (sitting and standing) as shown in Figures 7 and 8 are presented below in Equations (23) and (24), respectively [5]: The polynomial equations generated from two different postures of the lower limb, such as sitting and standing, are considered as the primary conditions of the artificial system stability analysis. The transfer function is generated from the polynomials to achieve the mathematical model [42] of the mentioned system. The closed loop transfer function is generated with unity feedback to gather information about the controller, performance and robustness parameters of the system, which shows stability as shown in Figure 4.
The characteristics equation is generated from the polynomials of the system at two different positions of the lower limb. The roots are determined using a computational technique. In the nonlinear domain, using Eigen vectors and the generated matrix, the phase portrait analysis has been incorporated. The system shows asymptotic stability in the nonlinear domain for both the roots of the characteristics equation as shown in Figure 5. For better and stronger decision making on the stability of the system, the Lyapunov stability [43] analysis method has been performed with the equilibrium point detection. This method also shows asymptotic stability where the system is stable and roots are in a complex domain as shown in Figure 6.

Results and Discussions
In this present research work, a 3 DoF lower limb has been analyzed for the artificial electro-mechanical system development. The conditions of movement of the human body such as sitting and standing are the main concern in this discussion. Non applicability of the superposition and homogeneity property gives a greater emphasis on the analysis and design of the control systems in the nonlinear domain [15]. The controllability and observability [44] cannot be determined simply based on rank tests. The concept of system stability [45] and checking using Phase portrait and Lyapunov stability [17] has been incorporated to make it a complete evaluation. Polynomials for the two positions of the lower limb (sitting and standing) as shown in Figures 7 and 8 (23) and (24), the open loop transfer function is obtained, as given below in the Equation (25) The output graphical presentation of the step response of the open loop transfer function of the system has been shown in the Figure 9.   (23) and (24), the open loop transfer function is obtained, as given below in the Equation (25) The output graphical presentation of the step response of the open loop transfer function of the system has been shown in the Figure 9. See Appendix A. f 2 (t) = t 2 +19t + 18 See Appendix A. From Equations (23) and (24), the open loop transfer function is obtained, as given below in the Equation (25) The output graphical presentation of the step response of the open loop transfer function of the system has been shown in the Figure 9.
The closed loop transfer function is given below in Equation (26) The output graphical presentation of the closed loop transfer function of the system has been shown in Figure 10. The PID (Proportional-Integral-Derivative) controller has been used to tune the closed loop transfer function of the lower limb system. The controller parameters of the closed loop transfer function of the system are shown in Table 1 and the performance and robustness parameters of the closed loop transfer function of the system are shown in Table 2. It is seen that the system is underdamped in Figure 10. The status related to the stability analysis of the mentioned system is prominent from the graphical presentation in Figure 10. The red line is representing the lower limb closed loop transfer function presented in Equation (26) in this paper. Comparing the reference signal represented through the blue line in Figure 10 and the lower limb system output presented in Equation (26), it is observed that the red line, the lower limb system output signal, has reached the stable condition with respect to the blue line at the settling point of the graph. The representations of the red and the blue lines are shown in the legends in Figure 10. The closed loop transfer function is given below in Equation (26) derived from Equation (25) The output graphical presentation of the closed loop transfer function of the system has been shown in Figure 10. The PID (Proportional-Integral-Derivative) controller has been used to tune the closed loop transfer function of the lower limb system. The controller parameters of the closed loop transfer function of the system are shown in Table 1 and the performance and robustness parameters of the closed loop transfer function of the system are shown in Table 2. It is seen that the system is underdamped in Figure 10. The status related to the stability analysis of the mentioned system is prominent from the graphical presentation in Figure 10. The red line is representing the lower limb closed loop transfer function presented in Equation (26) in this paper. Comparing the reference signal represented through the blue line in Figure 10 and the lower limb system output presented in Equation (26), it is observed that the red line, the lower limb system output signal, has reached the stable condition with respect to the blue line at the settling point of the graph. The representations of the red and the blue lines are shown in the legends in Figure 10.    (26)  The output graphical presentation of the closed loop transfer function of the system has been shown in Figure 10. The PID (Proportional-Integral-Derivative) controller has been used to tune the closed loop transfer function of the lower limb system. The controller parameters of the closed loop transfer function of the system are shown in Table 1 and the performance and robustness parameters of the closed loop transfer function of the system are shown in Table 2. It is seen that the system is underdamped in Figure 10. The status related to the stability analysis of the mentioned system is prominent from the graphical presentation in Figure 10. The red line is representing the lower limb closed loop transfer function presented in Equation (26) in this paper. Comparing the reference signal represented through the blue line in Figure 10 and the lower limb system output presented in Equation (26), it is observed that the red line, the lower limb system output signal, has reached the stable condition with respect to the blue line at the settling point of the graph. The representations of the red and the blue lines are shown in the legends in Figure 10.   In Figure 10, Tables 1 and 2 are representing the tuned graphical analysis of the lower limb transfer function presented in Equation (26), using the traditional PID controller. The output performance is not very satisfactory here. So, an advanced controller, the "Control System Designer", is preferred in order to show the better-tuned output presented in Figure 11 of the lower limb transfer function.

The closed loop transfer function is given below in Equation
In Figure 10, Tables 1 and 2 are representing the tuned graphical analysis of the lower limb transfer function presented in Equation (26), using the traditional PID controller. The output performance is not very satisfactory here. So, an advanced controller, the "Control System Designer", is preferred in order to show the better-tuned output presented in Figure 11 of the lower limb transfer function. Figure 11. Graphical representation of the closed loop transfer function of the system using an advanced control system designer. Therefore, the decision making regarding the system stability is possible observing the whole performance of the advanced controller output. In Table 3, the related characteristic parameters are shown. This is clearly shown in that the peak value and the settling time is reduced to a certain level. Modern technology requires advanced control laws to meet the design specifications for humanoid limb design, thus highlighting the position of nonlinear control systems with recent advancements in the artificial lower limb design. The application area in system design and analysis means that many technologies are developed based on this advanced nonlinear control system stability analysis method. The non-linear system stability analysis is processed as discussed below.
The generated matrix from Equations (23) and (24) is as given below: Then where, λ = eigen value, I = Identity matrix. Suppose X = A − λI Now, det(X) = λ 2 − 20λ − 6.6354 (29) det(X) is called the characteristic equation of the matrix A. Since the matrix A − λI is singular, det(X) = 0. Therefore, The roots of the equation are λ = 10 ± 10.3264. Therefore, the two roots are λ 1 = 20.3264 and λ 2 = −0.3264. For λ 1 = 20.3264 from matrix X, the obtained value is given below: Using the Gauss elimination method, the obtained simplified equation is Or, For λ 2 = −0.3264 from matrix X, obtained value is given below: Using the Gauss elimination method, the obtained simplified equation is Now from Equation (27)  The produced graphical phase portrait [10] plot is given below in Figure 11. the first phase plot is going outward or away from the origin of the plane and the second phase plot is coming towards the origin of the plane as shown in Figure 12. These conditions of phase portrait plot proves an asymptotic stability. The phase plot lines are presented through red lines in Figure 12.
Now from Equation (27) and using the roots = 20.3264 and = −0.3264. The produced graphical phase portrait [10] plot is given below in Figure 11. This output shows asymptotic stability because for Eigen vectors the first phase plot is going outward or away from the origin of the plane and the second phase plot is coming towards the origin of the plane as shown in Figure 12. These conditions of phase portrait plot proves an asymptotic stability. The phase plot lines are presented through red lines in Figure 12.  The theory of Lyapunov stability is a standard theory and one of the most important mathematical tools in the analysis of non-linear systems. Now, for further stability checking, the Lyapunov stability method is strongly recommended for nonlinear systems. The followed steps are given below for stability checking: Step 1: Representation of polynomials in characteristics equation form.
Step 2: Determination of the equilibrium point of the system.
Step 3: Generation of Jacobian matrix to implement in the characteristics equation.
Step 4: Determination of roots from the characteristics equation.
Step 5: Verifying the stability conditions from the roots.
Step 6: Declaration of Asymptotic stability condition. Now the process of Lyapunov stability method analysis is as shown below: From Equation (23) it can be written as From Equation (24) it can be written as After applying a derivative process in the Equations (36) and (37) respectively, the obtained Equations are as given below: 2y + 25.6354 = 0 (38) 2y Now, putting in values from Equation (35), the obtained Equations are as given below: After applying the derivative process again in the Equations (40) and (41) respectively, the obtained Equations are as given below: 2 . .
x 1 = 0 or it can be written as Now, to find the equilibrium point, it is similar as .
x 2 = 0 or it can be written as Therefore, the equilibrium point is at (0, 0) from Equations (43) and (44). Now, the Jacobian matrix can be formed using the formula as given below in the Equation (45): Or, λ 2 + λ + 1 = 0 Therefore, the roots are λ 1 , λ 2 = − 1 2 ± j 1 √ 3 . The abovementioned system has roots containing real and imaginary parts. Therefore, the system can be declared as asymptotically stable. The system performance efficiency is checked through this method.
This method is involved to ensure stability of the electro-mechanical system in the field of behavioral quality improvement of the designed artificial model. According to the above results shown through the mathematical models, the nonlinearity nature present in the lower limb system is a dead zone combined with a saturation condition. Figure 13 represents the lower limb system characteristics of overall nonlinearity. through this method.
This method is involved to ensure stability of the electro-mechanical system in the field of behavioral quality improvement of the designed artificial model. According to the above results shown through the mathematical models, the nonlinearity nature present in the lower limb system is a dead zone combined with a saturation condition. Figure 13 represents the lower limb system characteristics of overall nonlinearity. The nonlinear simulation of the lower limb system mentioned in Equation (26) is simulated using the dead zone and saturation functions to show satisfactory output in this paper for the overall simulation of the system. The graphical presentations are presented in Figures 14 and 15, respectively.  The nonlinear simulation of the lower limb system mentioned in Equation (26) is simulated using the dead zone and saturation functions to show satisfactory output in this paper for the overall simulation of the system. The graphical presentations are presented in Figures 14 and 15, respectively. This method is involved to ensure stability of the electro-mechanical system in the field of behavioral quality improvement of the designed artificial model. According to the above results shown through the mathematical models, the nonlinearity nature present in the lower limb system is a dead zone combined with a saturation condition. Figure 13 represents the lower limb system characteristics of overall nonlinearity. The nonlinear simulation of the lower limb system mentioned in Equation (26) is simulated using the dead zone and saturation functions to show satisfactory output in this paper for the overall simulation of the system. The graphical presentations are presented in Figures 14 and 15, respectively.  The graphical scenario depicts the lower limb system performance regarding nonlinearity. In Figure 14, initially the signal has started to show the activities but after a certain time it has reached the dead zone or the "no action" condition. In Figure 15, initially the signal has started to show output but after a certain period it has reached to saturation condition, being stable at a fixed level. These dead zone and saturation conditions are the specified nonlinear characteristics of the lower limb function only. The graphical scenario depicts the lower limb system performance regarding nonlinearity. In Figure 14, initially the signal has started to show the activities but after a certain time it has reached the dead zone or the "no action" condition. In Figure 15, initially the signal has started to show output but after a certain period it has reached to saturation condition, being stable at a fixed level. These dead zone and saturation conditions are the specified nonlinear characteristics of the lower limb function only.
The artificial limb design needs to be developed on the basis of knowledge-based system development. This process demands some expertise in the domain of the embedded system. This also requires the intelligent incorporation of an input-output database for the trained functional implementation of humanoid movement with an artificial lower extremity. Modern approaches for nonlinear control systems can be represented as an intelligent control. Various methods can be utilized to develop an automated system using neural networks, machine learning, and fuzzy logic, etc.

Conclusions
The nonlinear control systems have gained a challenging position in the technological development of the biomedical application field. In this recent work, the forward dynamics of the robotic system with the DoF concept has been explained using the Euler-Lagrange method. The idea of the mathematical model representation has been demonstrated with the polynomial equation generation from the experimental setup of the lower limb model. The comparative analysis has been incorporated for the system performance analysis using the traditional PID controller and an advanced control system designer. In future, there are lots of opportunities to use a more advanced application of controllers. The system performance stability is checked in two ways. It is found to be asymptotically stable using phase portrait analysis and the most important Lyapunov stability method. The lower limb shows nonlinear behavior. The system remains in the dead zone initially and after a certain period it starts operating. Then it reaches to the saturation level. The overall nonlinear behavior of the lower limb system has been shown here. The role of the describing function is the most important factor to determine the system nature selection. In future, the control system simulation approach can be implemented to get a more informative system identification view of the artificial system design.