Rehabilitation Robotics 2013View this Special Issue
Research Article | Open Access
Soumya Kanti Manna, Subhasis Bhaumik, "A Bioinspired 10 DOF Wearable Powered Arm Exoskeleton for Rehabilitation", Journal of Robotics, vol. 2013, Article ID 741359, 15 pages, 2013. https://doi.org/10.1155/2013/741359
A Bioinspired 10 DOF Wearable Powered Arm Exoskeleton for Rehabilitation
The developed exoskeleton device (Exorn) has ten degrees of freedom to control joints starting from shoulder griddle to wrist to provide better redundancy, portability, and flexibility to the human arm motion. A 3D conceptual model is being designed to make the system wearable by human arm. All the joints are simple revolute joints with desired motion limit. A Simulink model of the human arm is being developed with proper mass and length to determine proper torque required for actuating those joints. Forward kinematics of the whole system has been formulated for getting desired dexterous workspace. A proper and simple Graphical User Interface (GUI) and the required embedded system have been designed for providing physiotherapy lessons to the patients. In the literature review it has been found that researchers have generally ignored the motion of shoulder griddle. Here we have implemented those motions in our design. It has also been found that people have taken elbow pronation and supination motion as a part of shoulder internal and external rotation though both motions are quite different. A predefined resolved motion rate control structure with independent joint control is used so that all movements can be controlled in a predefined way.
It is seen that a major stroke is viewed by more than half of those at risk as being worse than death. Paralysis is caused due to complete loss of muscle function. In both cases patients have stiff muscles which restrict them from any physical movement of the affected part. Even the patients suffering from spinal cord injury or several nerve diseases may also lose their muscle strength gradually. According to Krakauer  the degree of improvement at 6 months is best predicted by the motor deficit at 1 month despite standard rehabilitative interventions in the ensuing 5 months. It is already proved that if they are under the process of rehabilitation for several months after stroke, their active range of motion as well as muscle strength can increase significantly. The training includes all the orthopedic and neurological lessons so that it is effective to the human muscle treatment. The rehabilitation training is generally performed by a physiotherapist, but the duration of the training is not adequate due to the fatigue of the therapist. It is observed that a physiotherapist can perform the training for 8–10 hours a day and can provide service to 2 to 3 patients at a very high cost. They may even omit certain exercises which are essential to the patients. The person also does not get repetitive and accurate movement sessions in case of a manual interaction. Considering time, cost, repeatability, reproducibility, and accuracy, robotically assisted rehabilitation process is far better than human assisted training process. Several assistive devices like ARMin, MEDARM, and SMA are already being developed for the rehabilitation process. But each one of them has several advantages and disadvantages. The purpose of Exorn is to provide maximum degree of freedom so that the device becomes as redundant as human hand and can provide better rehabilitation process in future.
2. The Literature Review
Until now two types of design are being proposed. First one is an exoskeleton structure connected to a base platform from where all the actuators are being controlled taking the platform as a reference. Second one is a portable structure wearable by a human. First one is a bulky system like ARMin , MEDARM , and Sarcos Master Arm  where the motion is stable, but it cannot be effective as an assistive device for the paralyzed people. Also there is a problem of misalignment of the system with the human hand at the time of joint movement due to which there is a possibility of injury . In some design the number of degree of freedom is restricted to certain directions . It will omit certain necessary movements which are very much needed for the proper rehabilitation process including solder griddle movement. Hydraulic and pneumatic actuators/devices like NEUROExos  have a big cylinder and pump connected with it to provide actuation signal. This kind of system is impossible to relocate with the human movement to make system portable.
Mistry et al. have already carried out experiments on joint space force field trajectories  to estimate kinematics and dynamics of the human arm with several trails in case of 3D movement of the human arm. Parallel manipulators have a problem of high stiffness and range of motion with several restrictions. In case of a serial manipulator if the offset between actuators is less then it acts like a joint with several DOFs. Yang et al. control puma manipulator using arm exoskeleton because control system is quite the same for both the cases . Pneumatic and hydraulic system are used for actuation. electroactive polymer (EAP) has been used for haptic device. Ball et al. have developed actuation system using cables which in reverse increases the stiffness factor by 10 to 12 . So it is expected that motor driven system should be used to make the exoskeleton robot simple. Researchers are also trying to control the actuation through EMG sensors using fuzzy algorithm . Generally impedance control is usually implemented in which the Cartesian forces applied at the joint are converted into joint torque commands using the Jacobian  matrix. The main disadvantage of impedance control is that good force replication requires compensation of the natural dynamics of the exoskeleton, such as gravity loading and drives frictions. The Exoskeleton Arm-Master  and the L-EXOS Exoskeleton  are classic examples of exoskeletons that use this approach. An alternative approach of joint control is called “admittance” control which is primarily used to control manipulators used as large reach haptic devices. However, it has the major drawback of instability for high admittance (low impedance), which is the opposite of impedance control . The Sensor Arm  is an example of an exoskeleton implementing this approach.
3. Mechanical Model
The whole structure is made of Perspex which is lighter in weight and much hardy than any other similar material. High torque DC geared motor is used for controlling the joint movements of shoulder griddle, glenohumeral joint, and elbow joint for gross positioning. For the motion of wrist joint small BLDC servo motor is being used for fine positioning. Gear drives are made of Nylon 101 material. As there are some joint limitations because of mechanical constraints, it cannot give the same motion as human joint range as shown in Table 1. In few locations suspensors are used to sustain any jerk and to incorporate the variation of link. Human joint of two or more degrees of freedom is being replaced by serial manipulator with several actuators. The offset between those actuators is less so that it can act like a single joint actuator with several degrees of freedom. The axis of rotation of human joint has been considered for placement of the actuators. In this model, the length of upper arm and forearm supporting link can be varied according to the different human arm length. The CAD model of Exorn and manufactured model are shown in Figures 1 and 2. The other parts of exoskeleton structure and its manufactured models are shown in Figures 3, 4, 5, 6, 7, 8, 9, and 10. From those given figures the design of any particular joint model can be analyzed.
Shoulder griddle is one of the vital joints with two degrees of freedom. Actuator 1 is used for giving the elevation and depression motion in vertical direction and actuator 2 is used for scapular abduction and adduction in horizontal direction. For elevation and depression screw nut mechanism is being used. Glenohumeral joint is a ball-socket joint having three degrees of freedom. The actuators required for giving the motion to this joint are shown in Figure 5. Actuators 3 to actuator 5 are DC motors used for three different rotations. Actuator 6 is a BLDC motor used for shoulder internal and external rotation. One hollow cylindrical spur gear meshes with a small pinion gear which is driven by actuator 6. The length between shoulder joint and elbow is changeable. Actuator 7 is a DC motor used for controlling elbow flexion and extension as shown in Figure 7. The same kind of technique is applied for forearm pronation and supination as that happened in the case of shoulder internal and external rotation. A suspensor is attached to upper portion of elbow joint to make it jerk-free and incorporate the change in arm length. Wrist joint has two degrees of freedom in two planes. Actuators 9 and 10 are used for providing those motions. Here we can see that all the actuators are dedicated to give the motion in all directions according to the different movements of human arm as shown in Figure 11.
Weight of the Prototype: 5 kg.
Specification of the DC Geared Motor Weight—500 gm Speed—10 rpm Torque—100 kg-cm at 4 ampere (Max).
Specification of BLDC Motor Company-Faulhaber Weight—500 gm Speed—10 rpm Torque—100 kg-cm at 4 ampere (Max).
Gears are made up of Nylon 101. Structures are used for twisting moment at human elbow and wrist with the help of gear arrangement which provide better flexibility and strength.
Specification of Nylon 101 Tensile strength: 79289709 N/m2 Mass density: 1150 Kg/m3 Elastic modulus: 1000000000 N/m2 Yield strength: 60000000 N/m2 Turn ratio at the elbow joint: 35 : 8 Turn ratio at wrist joint: 15 : 4.
At the shoulder joint bush is used to provide the support to the joint motor axis. It is made up of brass to reduce the wear and tear between perspex and motor spindle. One cylindrical rod of the same material moves inside the bush and provides the support to the other end of structure at the time of joint flexion and extension.
The whole body of Exorn is made up of perspex. It is hardy and elastic. It gives less weight in respect to aluminum. Perplex sheet is cut into different pieces and attached with each other through industrial gum and nut (M4 thread). Three different thickness perspex sheets have been used for the structure. Higher width is used at the arm structure to give better support and the lower one is used for hand structure: Type 1—Width 5 mm Type 2—Width 10 mm Type 3—Width 20 mm.
Specification of Perspex (TM) GS Acrylic Cast Sheet Tensile strength: 75000000 N/m2 Density: 1190 Kg/m3.
Size of the Exoskeleton. The size of the exoskeleton (approximately 96.73 cm) has been taken as per the standard length of a human arm and it can be varied from person to person: Length of shoulder griddle: Length of upper arm: Length of forearm: Length of hand: .
Here is defined as the height of the human arm and if we consider the height of a standard human as 170 cm, the length will be approximately 96.73 cm. Obviously the length of human arm is different from person to person. But as we have made a physical system, we have taken the standard length. But the system has a provision so that it is wearable to all the people of variable arm length which has been discussed in the portion of mechanical model.
4. Kinematic Analysis of Exorn
Kinematic model of Exorn is designed according to the anthropometric data of standard human arm. The mechanical structure is such that it can align with human arm and can provide the desirable joint motion. Some mechanical barrier is located to restrict the joint from hyperextension. Here the forward kinematics of the whole system is discussed which provides the end point position and orientation due to several joint motions. It is basically acting as a series manipulator. The kinematic relation of a manipulator is represented by .
is the dimensional vector specifying the task velocity of end effector. is dimensional vector representing the joint velocity. is the by Jacobian matrix. For a redundant manipulator is always greater than . Here the kinematic redundancy is 3.
4.1. D-H Parameter of the 3D Model
To do kinematic analysis of the whole system it is essential to derive D-H parameter of the whole system. The coordinate system of the whole structure is shown in Figure 12. The coordinate structure of the exoskeleton system is based on 3D conceptual model of Exorn. D-H parameters of the structure are shown in Table 2. It is essential to calculate the position and orientation of the end point in 3D space with respect to the reference coordinate.
4.2. Transformation Matrix Generation
For getting the position and orientation of the end point of Exorn system homogeneous transformation matrix is considered. Each homogeneous transformation matrix was calculated for the position and orientation of present frame with respect to the previous frame. It is a matrix. The column is position of the end point and is the orientation of end point in 3D space. The transformation matrix is: The homogeneous matrix is being explained as For calculating the homogeneous matrix due to 10 degrees of freedom the 10 homogeneous matrixes are been multiplied sequentially. The final matrix is which is represented as The calculation can be made easy by dividing the total multiplication into two parts such that . As it is a 10 DOF system. Finally matrix of the end point of Exorn is represented with respect to reference and its upper 3 values of last column are position vector of end point related to world frame.
5. Simulation Model of Human Arm
A human arm model is being made in Simulink to determine the required torque to move a joint. The length of upper arm, lower arm, and hand is considered according to the anthropometric data of standard human arm. The inertia matrix of each segment is determined according to the mass of that segment and length of the centre of gravity from the proximal end as per Table 3. In Simulink, it is possible to provide motions in terms of velocity to joint actuator and from the joint sensor torque level change is noted. The Simulink model of human arm is shown in Figure 13. Maximum amount of torque is needed to move those joints which are in vertical direction. Here only joint motion related to shoulder, elbow, and wrist flexion and extension is considered. Next those internal and external rotations are being taken care of. In exoskeleton system those motions cannot be controlled from the fixed end due to its twisting type. Therefore the end effector should be rotated with respect to fixed end.
The Simulink blocks used for making the human arm model are shown in Figures 14 and 15. A joint actuator block is connected to provide the actuation signal and joint sensor block is attached to measure the computed torque. In order to make a joint with several DOFs, several bodies of zero mass and zero length with actuator block of different axis of rotations are connected. Here in Figure 15 body 1 and body 2 have zero length and zero mass. Only body 3 is the actual body. It is under every human arm joint. Body sensor is attached at the end point to get points in 3D space at the time of motion.
6. Embedded System
A PC based control system is being developed to provide particular orthopedic lesson. First of all command signal from PC is transmitted to the master controller. Master controller activates those joints which are needed to achieve the desired motion. It sends the command bits in a package through serial communication. Enable signal is being sent to all the local controllers, but only few are being activated as per command with the help of multiprocessor communication technique. All the motor controllers are connected in a daisy chain network that means and pins of master controller are connected to all and pins of each controller, respectively. Each controller has a particular node as an activating signal. When the master controller sends the activation signal, the particular motor controller with that particular node is activated.
The interfacing circuit between microcontroller and motor driver is shown in Figures 16(a) and 16(b). Driver IC has 4 input pins and 2 enable pins. Here all the inputs pins and enable pins are connecting to port 0 and port 1. Here all the enable pins are sorted and connected through one port pin. In this case all the enable pins become activated at one time. So the motion of motor can be controlled by changing driver inputs. Those inputs are amplified from 5 V to 12 V to run the motor in clockwise and anticlockwise manner (Figure 17). As shown in the figure one dedicated 7805 IC (Voltage regulator) is connected to microcontroller for 5 V power supply. A serial IC MAX232 is connected to the microcontroller to convert the RS232’s signals to TTL voltage levels that is acceptable to 8051 microcontroller and pins.
A bidirectional H bridgeDC motor control circuit is shown in Figures 17(a) and 17(b). The circuit is based on the IC L298 from ST Microelectronics L298. It is a dual full bridge driver that has a wide operating voltage range and can handle load currents up to 3A. The IC also features low saturation voltage and over temperature protection. In the circuit diodes D1 to D8 are protection diodes. The state of the motor depends on the logic level of the pins 5, 6, and 7 and 10, 11, and 12 and it is described in the table shown below the circuit diagram.
The diodes D1 to D4 and D5 to D8 provide a safer path for the back e.m.f from the motor to dissipate and thus it protects the corresponding bipolar transistors from damage. The input pins from motor driver are being connected to the microcontroller port pins. When the enable pins of the driver pin are active, then the input signal of the motor driver is amplified to 12 volt and sent to the motor. According to the configuration of port pins, motor moves in clockwise or anticlockwise direction.
In case of BLDC motor controller, every controller has its particular node for distinguishing one from the other. It is possible to control the speed of the motor either by varying the voltage supplied to it or by sending the serial command to it. The speed of the motor is controlled by sending a string of command to the BLDC motor controller serially. Also there is an encoder inside the BLDC motor for getting the angle value simultaneously in a feedback loop.
All the supply voltages are being bypassed through those junction circuits. Actually all voltage sources are being sorted. A 7812 IC voltage regulator is used to convert the 24 V into 12 V DC supply. When the diode connected in junction circuit is in forward bias, it sends the power supply to all driver circuits, but in the opposite case current from driver circuit cannot be sent to the controller circuit. The flow chart (Figure 18) is showing the logic behind the algorithm.
7. Graphical User Interface Model of Exorn v.1
A proper Graphical User Interface (Figure 19) has been developed for controlling all the motions serially through PC based system. It is being designed in Microsoft visual studio platform. The whole GUI model is divided into four divisions. First one is the configuration of serial communication where user has to select the particular COM port and baud rate. It has two buttons for opening and closing of COM port. After opening of COM port, any command can be sent to the serial port. There are three divisions regarding the motion of human joint for proper treatment. Here only one control can be selected out of three options. Selection of one option will freeze the others. The flow chart behind the GUI model is shown in Figure 20.
7.1. Orthopedic Lessons
There are some movements which are generally prescribed by different doctors for stroke patients for different problems. All the motions are stored in memory. After selecting some motion, the patient presses the start button and a command will be sent through serially for taking proper action. It will be updated with time according to orthopedic doctor advice.
7.2. Isolateral Exercises
In isolateral exercises there are some predefined exercises which are standard motions of all joints of human arm. Patients can follow those exercises so that their muscle rigidity becomes less. After each selection a command is being sent serially to master controller. Master controller always takes care of proper actions.
7.3. Angle Oriented Control
In this configuration each joint is controlled in clockwise or anticlockwise manner with different speed. More than one joint can be governed through it. When user selects joint motions, then a frame of command bit serially sends with a start and stop bit to distinguish between different commands. User has to write the required speed beside the box of direction of rotation. In case two motions are contradictory, then the controller does not allow that motion. It is already preprogrammed in the controller.
8. Predefined Control System
It is a predefined control (Figure 21). Already a vast schedule of exercise modules is to be loaded in the memory. When someone wants to perform a task, it is being sent to the controller through a manual switch or PC. Master controller will generate a reference trajectory for that particular task. Then it is being sent to the local controller which takes care of that particular trajectory. Always there is feedback loop connected between the system and the master controller. With the help of resolved motion rate control (Whitney 72) approach each joint can be controlled in an independent way. According to forward kinematics the kinematic model can be expressed using the equation Outside the singularities the inverse kinematics can be expressed as There is a relation between task space and joint space. Here the current position at some angle can be expressed as The desired location in 3D space is defined as . So error between current and desired position can be expressed as The next will be expressed as . In this way value can be increased or decreased in an iterative way so that the desired location can be reached. The error between desired location and current location can be compared using forward kinematics in a feedback loop.
In case of finding inverse kinematics there is a problem of singularity. So we are imposing damped least square method for calculating the inverse kinematics. The damped least squares solutions are intended to ensure a well-conditioned matrix for the solution algorithm while limiting the size of the solution vector . This is done by adding a diagonal matrix, , to the matrix . The approximate solution of damped least squares method can be expressed as , where can be expressed as the equation given below: When this equation is used then the solution satisfies where denotes the Euclidean norm. For the damped least squares solution will become pseudoinverse solution. Then the Jacobian inverse matrix will be , where become expressed as .
In the damped least squares method, the size of error in the task space is weighed against the size of the resulting solution. For a given the size of the solution vector is decreased by increasing . However, this is done at the expense of using an approximate solution. As increases so does the size of the error in the task space. A large has the other benefit of ensuring a well-conditioned matrix for inversion. It has been shown that the condition number, , of the matrix , is . Here are the singular values of the Jacobian matrix. It can be seen in this equation that is made arbitrarily close to 1 by increasing , thus ensuring a well-conditioned matrix for inversion even in a singular configuration, where . Simply stated, if one is willing to give up the exactness of the solution then the size of the joint rate can be reduced and a well-conditioned matrix for inversion can be ensured by increasing .
9. Simulation Results of Torque Level of General Human Arm
Maximum amount of torque is required to actuate human joints in vertical direction. Here in Figure 25 the torque needed to twist the elbow with respect to shoulder joint is shown. Figure 26 is showing the same kind of motion in wrist joint with respect to the elbow. Figures 22, 23, 24, 25, and 26 will show torque variation in different joints in vertical direction because of proper actuation signal. From those below graphs max. and min. torque for joint actuator can be determined because it is very much essential to select proper motor.
10. Workspace Calculation and Motion due to Forward Kinematics
For calculating the workspace of the robotic system, a MATLAB code of nested for loop of all DOFs is being taken with the joint range. It will plot all the reachable points of arm exoskeleton in 3D space. Figure 27 is showing the outer region covered by human hand. Figures 28, 29, 30, and 31 are showing the motion of exoskeleton model due to several joint movements.
From the above discussion it may be concluded that the developed system is useful for rehabilitation. Comparisons between human joint’s range and that of Exorn are already tabulated in Table 1 and the ranges are quite the same for both the cases. If the portable mechanical structure is attached to the human arm, it can be used as an assistive device for the paralyzed people. The 3D model of exoskeleton is already being developed. The model has been fabricated. With the help of the GUI all the motions can be controlled through serial communication. The developed prototype is under working condition. We would like to build the real system so that patient can be benefitted.
It is a video of the whole exoskeleton system. It is an experimental model that was built in the university laboratory. All the motions are controlled from PC through serial communication using the designed GUI (Graphical User Interface). The whole system can be controlled in three modes. These are as follows:
(i) Orthopaedic lessons are basically predefined exercises which are previously loaded into PC memory. These kinds of motions have been designed according to the lessons provided by orthopaedic doctor after any stroke or neuromuscular diseases. Operator has to select proper exercise for his/her treatment.
(ii) In the second mode, there are different isolateral exercises. Those are standard human motions having fixed range of movement. But it may differ based on the patient’s constraint.
(iii) Last one is the angle oriented control where we can control each actuator in clockwise or anticlockwise direction with different speed.
In this video we are basically controlling the human arm in 2nd and 3rd mode. First of all we are selecting the proper communication port and baud rate to establish the serial communication between PC and the controlling board of the exoskeleton system. By pressing the “open” and “close” button we can open the communication port and close it respectively.
To get internal and external rotation of glenohumeral joint, we have to provide motion in actuator no. 6. For that purpose, we are selecting the motor 6 following by its type of motion-clockwise or anticlockwise and different speed. After that, as soon as the ‘start’ key is pressed, the actuator will rotate the following part of human arm and stop the system by pressing the stop key. The speed can be changed by giving proper speed in the provided space for the speed of the motor in ‘angle oriented control’ mode. BLDC servo motors are used for those kinds of motions.
In the next portion, we are selecting elbow pronation and supination in ‘Isolateral Exercises’ mode. Here the range of motion is fixed because it is standard in range. For this kind of motion we are using the DC geared motor. We can also control the motion in slow, medium and fast range.
- J. W. Krakauer, “Arm function after stroke: from physiology to recovery,” Seminars in Neurology, vol. 25, no. 4, pp. 384–395, 2005.
- T. Nef and R. Riener, “ARMin—design of a novel arm rehabilitation robot,” in Proceedings of the IEEE 9th International Conference on Rehabilitation Robotics (ICORR '05), pp. 57–60, July 2005.
- S. J. Ball, I. E. Brown, and S. H. Scott, “MEDARM: a rehabilitation robot with 5DOF at the shoulder complex,” in Proceedings of the IEEE/ASME International Conference on Advanced Intelligent Mechatronics (AIM '07), September 2007.
- M. Mistry, P. Mohajerian, and S. Schaal, “Arm movement experiments with joint space force fields using an exoskeleton robot,” in Proceedings of the IEEE 9th International Conference on Rehabilitation Robotics (ICORR '05), pp. 408–413, July 2005.
- A. Schiele and F. C. T. van der Helm, “Kinematic design to improve ergonomics in human machine interaction,” IEEE Transactions on Neural System and Rehabilitation, vol. 14, no. 4, pp. 1534–4320, 2006.
- S. J. Ball, I. E. Brown, and S. H. Scott, “A planar 3DOF robotic exoskeleton for rehabilitation and assessment,” in Proceedings of the 29th Annual International Conference of IEEE-EMBS, Engineering in Medicine and Biology Society (EMBC '07), pp. 4024–4027, August 2007.
- I. Sardellitti, E. Cattin, S. Roccella et al., “Description, characterization and assessment of a bio-inspired shoulder joint-first link robot for neuro-robotic applications,” in Proceedings of the 1st IEEE/RAS-EMBS International Conference on Biomedical Robotics and Biomechatronics (BioRob '06), pp. 112–117, February 2006.
- Source-winter, Biomechanics and Motor Control of Human Motion, ISBN-978-0470398180, 2010.
- C.-J. Yang, B. Niu, J.-F. Zhang, and Y. Chen, “Different structure based control system of the PUMA manipulator with an arm exoskeleton,” in Proceedings of the IEEE Conference on Robotics, Automation and Mechatronics, pp. 572–577, December 2004.
- K. Kiguchi, T. Tanaka, K. Watanabe, and T. Fukuda, “Exoskeleton for human upper-limb motion support,” in Proceedings of the IEEE International Conference on Robotics and Automation, pp. 2206–2211, September 2003.
- T. H. Massie and K. J. Salisbury, “PHANToM haptic interface: a device for probing virtual objects,” in Proceedings of the International Mechanical Engineering Congress and Exposition, pp. 295–299, November 1994.
- C. Carignan and K. Cleary, “Closed-loop force control for haptic simulation of virtual environments,” The Electronic Journal of Haptics Research, vol. 2, no. 2, pp. 1–14, 2000.
- A. Frisoli, F. Rocchi, S. Marcheschi, A. Dettori, F. Salsedo, and M. Bergamasco, “A new force-feedback arm exoskeleton for haptic interaction in virtual environments,” in Proceedings of the First Joint Eurohaptics Conference and Symposium on Haptic Interfaces for Virtual Environment and Teleoperator Systems, pp. 195–201, 2005.
- D. Lawrence, “Impedance control stability properties in common implementations,” in Proceedings of the IEEE International Conference on Robotics and Automation, pp. 1185–1190, 1988.
- A. Nakai, Y. Kunii, H. Hashimoto, and F. Harashima, “Arm type haptic human interface: sensor arm,” in Proceedings of the International Conference on Artificial Reality and Tele-Existence (ICAT), pp. 77–84, Tokyo, Japan, December 1997.
Copyright © 2013 Soumya Kanti Manna and Subhasis Bhaumik. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.