Utilization of Function Generation Synthesis on Biomimetics: A Case Study on Moray Eel Double Jaw Design

Throughout history, humans have observed living or non-living things in nature and then imitated them in relation to these observations. This is due to the fact that the energy found in nature is generally consumed at an optimal level in order for it to endure. Biomimetic inspiration in many designs and applications is widely displayed, including within the field of engineering. In this paper, we were inspired by the double set of jaws found in the moray eel, which gives this fish a huge advantage while hunting, with a mobile pharyngeal jaw that works together with its oral jaw in order to overcome ineffective suction capabilities. A procedure that mimics the hunting motion of the moray eel was utilized by considering the overall movement as a single degree of freedom with multiple outputs on account of the repeating motion that is required during hunting. This procedure includes structural and dimensioning synthesis, wherein the latter was utilized with analytic kinematic synthesis for each linkage transfer. The flexibilities in parameters were taken into account with a novel multiple iterative kinematic synthesis algorithm that resulted in various mechanisms with the same purpose. Among the excessive number of resultant mechanisms, the optimization was carried out by considering the highest torque transmission ratio at critical timings that were specified as bio-constraints. In the end, the kinematic movement validation was utilized.


Introduction
Scientists have been trying to understand biological systems from various points of view throughout history. From an engineering point of view, the mechanical and electrical analogies of biological systems in technical works are said to have been started, at least, as long ago as Leonardo Da Vinci, Galvani, and Volta. The term "biomimetic" was born out of a need for those within technical fields to be able to describe the combination of some areas of biology with that of science or engineering wherein this combination wished to mimic the different functionalities of living or non-living creatures in the environment, including humans, with advances in material science, technology, and an understanding of biological systems [1]. Throughout these advances, the designing processes of many biomimetic machines that benefit humans have been inspired by nature [2,3]. It is undeniable that things in nature tend to use their energy at an optimum in order to endure. Thus, there is a duality between engineering and nature, as the performance of biologically inspired robots has robustness and sustainability [4]. Scientists have recently introduced this design approach of combining bio-inspiration to engineering students by adding it to curriculums as a new way to increase technical challenges [4]. Engineering companies should also implement biomimetic design processes into their procedures by changing organizational structures that, in the end, may combine the fields of engineering and biology together [5]. Biomimetic robots are robots that imitate not only humans but also animals that will eventually adapt to the environment so that they can learn and react faster, while, at the same time, may have novel mechanisms and manipulator structures according to their needs [6,7]. Over the past decade, many robotic devices, whose designs are based upon certain biological species in terms of user functionality in practical applications, have been developed, such as: humanoid robots [8,9] snake-like robots [10,11], swimming robots [12][13][14], flying robots [15,16], and even plants [17]. Moreover, there are studies where the functionality of certain parts of a biological species are mimicked without considering appearance, such as [18] and Liu et al.'s study [19] which is based on the compliant spine mechanism of the leaf-feeding insect, Motschulsky. In this study, they designed a wheeled wall-climbing robot and analyzed its performance with experiments and, eventually, concluded that climbing performance was improved by the bio-inspired compliant spine mechanism.
In this paper, the generated motion of the proposed mechanism was inspired by the combination of two different sets of jaws that have evolved in some fishes. While one of these jaws is the oral jaw that collects food, the other is the pharyngeal jaw normally located in the pharynx of the fish that transports the food into the esophagus [20,21]. Among these fishes that have this kind of jaw, moray eels have certain advantages in terms of hunting with respect to the point of view of this study. Although moray eels do not have effective suction feeding abilities, this disadvantage is overcome by a well-developed and wellpositioned pharyngeal jaw. More so than other fishes with a pharyngeal jaw, these creatures are capable of pulling prey into the esophagus with the help of their additional jaw, which gives them a huge advantage during hunting. Mehta and Wainwright [22] observed the moray eels by using high-speed video in a feeding sequence and they observed that the pharyngeal jaws are used together with the oral jaw. Moreover, it was noted that the upper and lower teeth are very convenient for extreme prey transport conditions during retraction and protraction. The advantage of this system overall appears especially when large prey is the target. In [23], Triyonoputro et al. designed and manufactured a double-jaw hand mechanism that mimicked the functionality of the moray eel jaw system. They imitated three motions of the mobile pharyngeal jaw: opening and moving forward, biting the prey, and pulling the prey with the double jaw and with their outer and inner gripper design. This overall hand system has four degrees of freedom (DOF) in total, with one of these dedicated to making the two grippers mobile.
When considering combining mechanisms with biomimetic studies, it is crucial to have knowledge about designing kinematic motions according to resultant tasks, which is, in this study, the motion of two sets of jaws dependently. Motion study of machines and mechanisms may be achieved by two methods of analysis and synthesis [16]. Simply, they can be described as studies that deal with mechanical systems in order to understand or fulfill certain motion characteristics. Kinematic synthesis is a design process for the functional dimensioning of mechanisms according to the required task, which allows the system DOF to be kept lower than the task degree of freedom [24]. The decision and procedures of the kinematic synthesis methods are the most important steps of the mechanical design procedure in terms of dimensioning the parts. Moreover, error characteristics of the designed mechanism are obtained by kinematic synthesis. Basically, there are three types of analytic kinematic synthesis methods used to designate the output characteristic of the mechanism, which are: function generation synthesis, body guidance synthesis, and path generation synthesis [25].
Function generation synthesis attempts to design the construction parameters of the mechanisms in order to provide a specified function relationship between the input and the output motion [26,27]. Path generation synthesis [28] guides a coupler point along a specified path and body guidance synthesis [29] is used when the body of the mechanism guides through a number of specified positions. Since these methodologies are based on keeping the DOF smaller than that of the tasks, it is not possible to obtain an exact solution, therefore, different approximation techniques are applied, and, eventually, motion is achieved with finite error. In applying this kind of approximation, precision points are defined regarding which motion is fulfilled. Overall, the number of the precision points must be equal to the number of the unknown construction parameters used in the synthesis [24].
In order to conclude a dimensioned mechanism or manipulator to serve as taskspecific by using a function generation synthesis methodology, as mentioned above, it is a requirement to have an input-output relation in order to fulfill the motion. This relationship has quite the effect on the resultant mechanism. On the other hand, if the required output paths are predictable, at least for some precision points but not for the input motion or range, then it is possible to construct a multiple input-output relationship within the permitted range of inputs and corresponding outputs. With an excessive number of input-output ranges, multiple solutions should be handled. Advances in computing technologies allow them to be used in different areas, including optimization computing. Ouezdou et al. [30] used an iterative optimization technique in kinematic synthesis with a rigid body guidance problem with six DOF manipulators. The Distributed Optimization Algorithm (DOM) was used to obtain analytical expressions for the minimal value of each parameter, whereas the dimensional parameters of the manipulator were determined. Without considering the kinematic synthesis methodology, in Kivela's thesis [31], a generic method for finding an optimal solution by using the Levenberg-Marguardt method was developed, eventually ranking and choosing the resultant robot according to the selected performance measure made possible. Patel and Sobh [32] developed a methodology to determine optimal manipulator configurations based on task descriptions. The idea behind this methodology was to get all of the possibilities in dimensioning that could reach the required poses, among which the best kinematic performance could be selected within the user-defined constraints with the least power consumption.
The purpose of this work is to design a mechanism that mimics the motion of two jaw sets of the moray eel with a single actuator. The oral jaw is driven by an inverted slider-crank mechanism that is responsible for initially catching the prey. The second set of the jaw (the pharyngeal jaw) pulls the prey inside through the esophagus, which is composed of two different mechanisms; anti parallel four-bar mechanism as the biting part and an offset slider-crank mechanism as the translating part between the throat and out end of the mouth. Analytic kinematic synthesis methodologies were applied in order to design a planar-type of mechanism by using an algorithm called a synthesis algorithm in order to have an iterative synthesis procedure that is made possible due to the flexibility of the parameters. Among different results, optimization is carried out in order to reach the most suitable mechanism to deliver the highest torque at a specified moment from the actuator to the outputs that are obtained with virtual work methods in each step. Finally, the mechanism is designed based on kinematics in order to validate the procedure.

Bio-Constraints and Structural Synthesis
The number of fishes that live in the ocean is a mystery. Eels are ray fin fish that belong to the Anguilliformes order, and moray eels (Muraenidae) are one of the members of this eel family. There are approximately 200 species of moray eel types in the world [33] and their heights can be 60-80 cm to 2-3 m [34]. In Figure 1, the moray eel's biological classification tree [34] is shown. achieved with finite error. In applying this kind of approximation, precision points are defined regarding which motion is fulfilled. Overall, the number of the precision points must be equal to the number of the unknown construction parameters used in the synthesis [24]. In order to conclude a dimensioned mechanism or manipulator to serve as task-specific by using a function generation synthesis methodology, as mentioned above, it is a requirement to have an input-output relation in order to fulfill the motion. This relationship has quite the effect on the resultant mechanism. On the other hand, if the required output paths are predictable, at least for some precision points but not for the input motion or range, then it is possible to construct a multiple input-output relationship within the permitted range of inputs and corresponding outputs. With an excessive number of input-output ranges, multiple solutions should be handled. Advances in computing technologies allow them to be used in different areas, including optimization computing. Ouezdou et al. [30] used an iterative optimization technique in kinematic synthesis with a rigid body guidance problem with six DOF manipulators. The Distributed Optimization Algorithm (DOM) was used to obtain analytical expressions for the minimal value of each parameter, whereas the dimensional parameters of the manipulator were determined. Without considering the kinematic synthesis methodology, in Kivela's thesis [31], a generic method for finding an optimal solution by using the Levenberg-Marguardt method was developed, eventually ranking and choosing the resultant robot according to the selected performance measure made possible. Patel and Sobh [32] developed a methodology to determine optimal manipulator configurations based on task descriptions. The idea behind this methodology was to get all of the possibilities in dimensioning that could reach the required poses, among which the best kinematic performance could be selected within the user-defined constraints with the least power consumption.
The purpose of this work is to design a mechanism that mimics the motion of two jaw sets of the moray eel with a single actuator. The oral jaw is driven by an inverted slider-crank mechanism that is responsible for initially catching the prey. The second set of the jaw (the pharyngeal jaw) pulls the prey inside through the esophagus, which is composed of two different mechanisms; anti parallel four-bar mechanism as the biting part and an offset slider-crank mechanism as the translating part between the throat and out end of the mouth. Analytic kinematic synthesis methodologies were applied in order to design a planar-type of mechanism by using an algorithm called a synthesis algorithm in order to have an iterative synthesis procedure that is made possible due to the flexibility of the parameters. Among different results, optimization is carried out in order to reach the most suitable mechanism to deliver the highest torque at a specified moment from the actuator to the outputs that are obtained with virtual work methods in each step. Finally, the mechanism is designed based on kinematics in order to validate the procedure.

Bio-Constraints and Structural Synthesis
The number of fishes that live in the ocean is a mystery. Eels are ray fin fish that belong to the Anguilliformes order, and moray eels (Muraenidae) are one of the members of this eel family. There are approximately 200 species of moray eel types in the world [33] and their heights can be 60-80 cm to 2-3 m [34]. In Figure 1, the moray eel's biological classification tree [34] is shown. The moray eels' set of jaws is the main difference that makes this research team want to study them. These creatures have two separate jaws; an oral jaw and a pharyngeal jaw. The moray eels' set of jaws is the main difference that makes this research team want to study them. These creatures have two separate jaws; an oral jaw and a pharyngeal jaw. None of the jaw sets have become blunt due to the fact that they give a huge advantage to those creatures under water and are very useful and frequently used during hunting ( Figure 2). None of the jaw sets have become blunt due to the fact that they give a huge advantage t those creatures under water and are very useful and frequently used during hunting (Fig  ure 2). In terms of mechanism structures, the oral jaw can perform an opening and closin motion and the pharyngeal jaw can additionally slide inside the oral jaw as a whole apar from its opening and closing motion. The hunting motion of the moray eel's overall jaw structure can be treated as a repeatable motion when the prey is being hunted. The ap proximate motion of the moray eel during hunting can be visualized as shown in Figure 3. In the nature of the moray eel's hunting, firstly, they squeeze their prey with thei oral set of jaws and make the prey ineludible. Then, the pharyngeal jaw comes into actio from the throat to the outside of the mouth in order to snatch the prey to pull it into th esophagus. Before pulling the prey into the esophagus, the oral jaw should release th prey in order for it to be pulled easily. Thus, the pharyngeal jaw of moray eel has a specifi motion characteristic working together with the oral jaw, which is combined with bot the translation to get out of the oral jaw and to pull back the jaw together with the prey t the esophagus, and the rotation to bite and squeeze the prey. In order to mimic the motio characteristics of the moray eel on a mechanism that is dependent on each other, bio-con straints that define different poses of the pharyngeal and the oral jaws to be fulfilled i terms of kinematics for this work can be visualized in Table 1. Table 1. Approximation of the different poses of the oral jaw and the pharyngeal jaw during hun ing.

Oral Jaw
Pharyngeal Jaw

Natural Pose
Slider to Backward, Jaw Closed-Opening In terms of mechanism structures, the oral jaw can perform an opening and closing motion and the pharyngeal jaw can additionally slide inside the oral jaw as a whole apart from its opening and closing motion. The hunting motion of the moray eel's overall jaw structure can be treated as a repeatable motion when the prey is being hunted. The approximate motion of the moray eel during hunting can be visualized as shown in Figure 3. None of the jaw sets have become blunt due to the fact that they give a huge advantage to those creatures under water and are very useful and frequently used during hunting (Figure 2). In terms of mechanism structures, the oral jaw can perform an opening and closing motion and the pharyngeal jaw can additionally slide inside the oral jaw as a whole apart from its opening and closing motion. The hunting motion of the moray eel's overall jaw structure can be treated as a repeatable motion when the prey is being hunted. The approximate motion of the moray eel during hunting can be visualized as shown in Figure 3. In the nature of the moray eel's hunting, firstly, they squeeze their prey with their oral set of jaws and make the prey ineludible. Then, the pharyngeal jaw comes into action from the throat to the outside of the mouth in order to snatch the prey to pull it into the esophagus. Before pulling the prey into the esophagus, the oral jaw should release the prey in order for it to be pulled easily. Thus, the pharyngeal jaw of moray eel has a specific motion characteristic working together with the oral jaw, which is combined with both the translation to get out of the oral jaw and to pull back the jaw together with the prey to the esophagus, and the rotation to bite and squeeze the prey. In order to mimic the motion characteristics of the moray eel on a mechanism that is dependent on each other, bio-constraints that define different poses of the pharyngeal and the oral jaws to be fulfilled in terms of kinematics for this work can be visualized in Table 1. Table 1. Approximation of the different poses of the oral jaw and the pharyngeal jaw during hunting.

Oral Jaw
Pharyngeal Jaw

Natural Pose
Slider to Backward, Jaw Closed-Opening In the nature of the moray eel's hunting, firstly, they squeeze their prey with their oral set of jaws and make the prey ineludible. Then, the pharyngeal jaw comes into action from the throat to the outside of the mouth in order to snatch the prey to pull it into the esophagus. Before pulling the prey into the esophagus, the oral jaw should release the prey in order for it to be pulled easily. Thus, the pharyngeal jaw of moray eel has a specific motion characteristic working together with the oral jaw, which is combined with both the translation to get out of the oral jaw and to pull back the jaw together with the prey to the esophagus, and the rotation to bite and squeeze the prey. In order to mimic the motion characteristics of the moray eel on a mechanism that is dependent on each other, bio-constraints that define different poses of the pharyngeal and the oral jaws to be fulfilled in terms of kinematics for this work can be visualized in Table 1. By means of considerations in the hunting of the moray eel with a visualization in Table 1, the whole mechanism has specific position and angle requirements at each instant of time that can be approximated as shown in Figure 4. The graphic in Figure 4 shows the rotational angle and translational stroke of the pharyngeal and oral jaws of the moray eel with respect to each other through the time. The solid red line is the angle of the pharyngeal jaw, the solid blue line is the oral jaw, the dashed red line is the stroke of the pharyngeal jaw, and the dashed blue line is the stroke of the oral jaw. In this work, the ultimate goal is specified with these approximate movement characteristics and can be treated as bio-constraints of the motion.
While the only rotational actuator that is placed in a proper point inside the biomimetic mechanism is driven with a constant rotational velocity, the required bio-constraints of the motions can be sorted as follows: If the pharyngeal jaw reaches out of the mouth, the pharyngeal jaw should be in a closed position (the action of biting). At the same time, the oral jaw should be in an opening position when the prey is just released.

•
If the pharyngeal jaw travels from out of the mouth to the throat, the pharyngeal jaw should be in a closed position (carrying prey to the esophagus). At the same time, the oral jaw should be in a closing manner to the natural pose, which is a closed oral jaw.
Before the structural design of the mechanism, functional design constraints should also be defined, which have a huge impact on the chosen structure. As is seen from the biomimetic motion characteristics, the overall motion of both jaws is dependent on each other during hunting. By considering the whole motion of both jaws, the mechanism may be constructed with multiple degrees of freedom, however, depending on the motion characteristics during the motion, it is highly possible to face control difficulties. For the sake of simplicity, this overall motion can be reduced to a single DOF motion including many sets of mechanisms that are driven by a common source. This creates dependency, which results in task-based results, however, it can be adapted to other automation applications without considering any relative control difficulties. The increase in overall safety is also an advantage in such a design due to the motion constraints.
Avoiding an increase in mechanical complexity, another functional design constraint, defined as the common actuation source, is a rotary electric motor that drives the whole The graphic in Figure 4 shows the rotational angle and translational stroke of the pharyngeal and oral jaws of the moray eel with respect to each other through the time. The solid red line is the angle of the pharyngeal jaw, the solid blue line is the oral jaw, the dashed red line is the stroke of the pharyngeal jaw, and the dashed blue line is the stroke of the oral jaw. In this work, the ultimate goal is specified with these approximate movement characteristics and can be treated as bio-constraints of the motion.
While the only rotational actuator that is placed in a proper point inside the biomimetic mechanism is driven with a constant rotational velocity, the required bio-constraints of the motions can be sorted as follows: If the pharyngeal jaw reaches out of the mouth, the pharyngeal jaw should be in a closed position (the action of biting). At the same time, the oral jaw should be in an opening position when the prey is just released.

•
If the pharyngeal jaw travels from out of the mouth to the throat, the pharyngeal jaw should be in a closed position (carrying prey to the esophagus). At the same time, the oral jaw should be in a closing manner to the natural pose, which is a closed oral jaw.
Before the structural design of the mechanism, functional design constraints should also be defined, which have a huge impact on the chosen structure. As is seen from the biomimetic motion characteristics, the overall motion of both jaws is dependent on each other during hunting. By considering the whole motion of both jaws, the mechanism may be constructed with multiple degrees of freedom, however, depending on the motion characteristics during the motion, it is highly possible to face control difficulties. For the sake of simplicity, this overall motion can be reduced to a single DOF motion including many sets of mechanisms that are driven by a common source. This creates dependency, which results in task-based results, however, it can be adapted to other automation applications without considering any relative control difficulties. The increase in overall safety is also an advantage in such a design due to the motion constraints. Avoiding an increase in mechanical complexity, another functional design constraint, defined as the common actuation source, is a rotary electric motor that drives the whole movement from an intersection point so that there will be no additional motion transferring other than the structural mechanisms. The moray eel's overall hunting characteristics will be integrated without any change during the mechanism's design.
Considering biomimetic and functional design constraints, the structure of the mechanism is chosen to be linkage type with different movement abilities in order to achieve multiple motions with kinematic synthesis methodology, since task mobility is greater than the targeted mobility.
Considering the motion of the pharyngeal jaw that has both translation and rotation motion, the structure of the mechanism has to contain some asymmetry due to the bioconstraints. This asymmetry may be allowable by adding some offset to a slider-crank mechanism. For this reason, a mechanism that combines a slider-crank mechanism with an offset for sliding and a four-bar mechanism for rotating was chosen, with the possibility of two different types ( Figure 5).
Biomimetics 2022, 7, x FOR PEER REVIEW 7 of 32 movement from an intersection point so that there will be no additional motion transferring other than the structural mechanisms. The moray eel's overall hunting characteristics will be integrated without any change during the mechanism's design.
Considering biomimetic and functional design constraints, the structure of the mechanism is chosen to be linkage type with different movement abilities in order to achieve multiple motions with kinematic synthesis methodology, since task mobility is greater than the targeted mobility.
Considering the motion of the pharyngeal jaw that has both translation and rotation motion, the structure of the mechanism has to contain some asymmetry due to the bioconstraints. This asymmetry may be allowable by adding some offset to a slider-crank mechanism. For this reason, a mechanism that combines a slider-crank mechanism with an offset for sliding and a four-bar mechanism for rotating was chosen, with the possibility of two different types ( Figure 5).
(a) (b) Figure 5. Kinematic chain representation of the pharyngeal jaw: slider-crank mechanism with (a) a four-bar mechanism; (b) an anti-parallel four-bar mechanism.
As seen in the kinematic chain representation of the structure, the actuation source of the movement is placed in the moray eel's body in a grounded position. The coupler link of the slider-crank mechanism is actually the input source for the translating four-bar mechanism itself. By considering the bio-constraints of the overall design process, as the slider moves out of the throat (+x direction in Figure 5a), which is possible by rotating the input source in a clockwise direction, the output of the sliding four-bar mechanism should rotate in the direction of the biting (clockwise direction) in order to catch the prey. In the other direction of the slider mechanism (-x direction in Figure 5a), the output of the sliding four-bar mechanism should keep its position closed (biting) in order to not allow the prey to flee. As a result, the structure of the sliding four-bar mechanism should be changed to an anti-parallel type in order to achieve this movement constraint ( Figure 5b). With this configuration, it is possible to achieve a closing action during the slider moves out of the mouth.
For mechanical simplicity and to avoid motion transmission with additional linkages or belts, the motion source of the overall mechanism is placed in the common point where it is the final actuator location. Considering the motion of the oral jaw, it should only achieve the rotational motion when necessary. According to the bio-constraints, in a onetime cycle of a full rotation of the rotary actuator, the rotational motion of the oral jaw should be utilized two times, one of which is to catch the prey with the oral jaw and the other is to keep the position as a natural pose, which is visualized in Table 1. On the other hand, as the functional design constraint, the mechanism should be a single DOF. As a result, the mechanism chosen to achieve this overall motion is the inverted slider-crank mechanism, where the input is the prismatic joint as seen from the kinematic chain representation in Figure 6. However, since it was planned to have a single rotational actuator on the system to drive all of the motion, this prismatic joint motion as input is given to the As seen in the kinematic chain representation of the structure, the actuation source of the movement is placed in the moray eel's body in a grounded position. The coupler link of the slider-crank mechanism is actually the input source for the translating four-bar mechanism itself. By considering the bio-constraints of the overall design process, as the slider moves out of the throat (+x direction in Figure 5a), which is possible by rotating the input source in a clockwise direction, the output of the sliding four-bar mechanism should rotate in the direction of the biting (clockwise direction) in order to catch the prey. In the other direction of the slider mechanism (-x direction in Figure 5a), the output of the sliding four-bar mechanism should keep its position closed (biting) in order to not allow the prey to flee. As a result, the structure of the sliding four-bar mechanism should be changed to an anti-parallel type in order to achieve this movement constraint ( Figure 5b). With this configuration, it is possible to achieve a closing action during the slider moves out of the mouth.
For mechanical simplicity and to avoid motion transmission with additional linkages or belts, the motion source of the overall mechanism is placed in the common point where it is the final actuator location. Considering the motion of the oral jaw, it should only achieve the rotational motion when necessary. According to the bio-constraints, in a one-time cycle of a full rotation of the rotary actuator, the rotational motion of the oral jaw should be utilized two times, one of which is to catch the prey with the oral jaw and the other is to keep the position as a natural pose, which is visualized in Table 1. On the other hand, as the functional design constraint, the mechanism should be a single DOF. As a result, the mechanism chosen to achieve this overall motion is the inverted slider-crank mechanism, where the input is the prismatic joint as seen from the kinematic chain representation in Figure 6. However, since it was planned to have a single rotational actuator on the system to drive all of the motion, this prismatic joint motion as input is given to the system with a cam profile, whose profile is extracted according to the timings of the oral jaw.
Biomimetics 2022, 7, x FOR PEER REVIEW 8 of 32 system with a cam profile, whose profile is extracted according to the timings of the oral jaw. As a result, the structures of both jaw mechanisms are chosen as an offset slidercrank mechanism combined with an anti-parallel four-bar mechanism for the pharyngeal jaw, and an inverted slider-crank mechanism with the input of a prismatic joint that is driven with a cam profile ( Figure 7).  As a result, the structures of both jaw mechanisms are chosen as an offset slider-crank mechanism combined with an anti-parallel four-bar mechanism for the pharyngeal jaw, and an inverted slider-crank mechanism with the input of a prismatic joint that is driven with a cam profile ( Figure 7). system with a cam profile, whose profile is extracted according to the timings of the oral jaw. Figure 6. Kinematic chain representation of the oral jaw: inverted slider-crank mechanism.
As a result, the structures of both jaw mechanisms are chosen as an offset slidercrank mechanism combined with an anti-parallel four-bar mechanism for the pharyngeal jaw, and an inverted slider-crank mechanism with the input of a prismatic joint that is driven with a cam profile ( Figure 7).  system with a cam profile, whose profile is extracted according to the timings of the oral jaw. As a result, the structures of both jaw mechanisms are chosen as an offset slidercrank mechanism combined with an anti-parallel four-bar mechanism for the pharyngeal jaw, and an inverted slider-crank mechanism with the input of a prismatic joint that is driven with a cam profile ( Figure 7).  Due to the design goal, both mechanisms are combined together with the common actuator. As a result, the mechanical structure of the pharyngeal jaw is composed of one offset slider-crank mechanism and one anti-parallel four-bar mechanism that is placed on the slider and driven by it ( Figure 9). Due to the design goal, both mechanisms are combined together with the common actuator. As a result, the mechanical structure of the pharyngeal jaw is composed of one offset slider-crank mechanism and one anti-parallel four-bar mechanism that is placed on the slider and driven by it ( Figure 9). As in the structural synthesis part, the oral jaw is constructed as an inverted slidercrank mechanism and this mechanism with labeled parameters can be seen in Figure 10. After the structural design is utilized, the dimensioning of the links comes into action as the next step. The functions that define the input-output relations should be revealed and the analytic kinematic synthesis procedure should be carried out.

Kinematic Synthesis of the Pharyngeal Jaw
In the proposed design, the motion of the pharyngeal jaw is generated by the combination of two different mechanisms; the anti-parallel four-bar mechanism and the offset slider-crank mechanism. The main mechanism is formed by the serial connection of these two mechanisms. Input motion acts on the crank (b1) of the slider-crank mechanism, and, thus, translational motion is generated. The translational motion of the slider provides the translational motion of the four-bar mechanism which, in turn, provides the rotational As in the structural synthesis part, the oral jaw is constructed as an inverted slidercrank mechanism and this mechanism with labeled parameters can be seen in Figure 10. Due to the design goal, both mechanisms are combined together with the common actuator. As a result, the mechanical structure of the pharyngeal jaw is composed of one offset slider-crank mechanism and one anti-parallel four-bar mechanism that is placed on the slider and driven by it ( Figure 9). As in the structural synthesis part, the oral jaw is constructed as an inverted slidercrank mechanism and this mechanism with labeled parameters can be seen in Figure 10. After the structural design is utilized, the dimensioning of the links comes into action as the next step. The functions that define the input-output relations should be revealed and the analytic kinematic synthesis procedure should be carried out.

Kinematic Synthesis of the Pharyngeal Jaw
In the proposed design, the motion of the pharyngeal jaw is generated by the combination of two different mechanisms; the anti-parallel four-bar mechanism and the offset slider-crank mechanism. The main mechanism is formed by the serial connection of these two mechanisms. Input motion acts on the crank (b1) of the slider-crank mechanism, and, thus, translational motion is generated. The translational motion of the slider provides the translational motion of the four-bar mechanism which, in turn, provides the rotational After the structural design is utilized, the dimensioning of the links comes into action as the next step. The functions that define the input-output relations should be revealed and the analytic kinematic synthesis procedure should be carried out.

Kinematic Synthesis of the Pharyngeal Jaw
In the proposed design, the motion of the pharyngeal jaw is generated by the combination of two different mechanisms; the anti-parallel four-bar mechanism and the offset slider-crank mechanism. The main mechanism is formed by the serial connection of these two mechanisms. Input motion acts on the crank (b 1 ) of the slider-crank mechanism, and, thus, translational motion is generated. The translational motion of the slider provides the translational motion of the four-bar mechanism which, in turn, provides the rotational motion of the pharyngeal jaw. To ensure the constraint that the overall mechanism is actuated by a single actuator, the rotational motion input of the four-bar mechanism is provided by one of the links of the slider mechanism. This situation is ensured by the rigid connection of the slider's coupler link (c 1 ) and the four-bar link (a 2 ). Thus, the motion of the slider makes (a 2 ) rotate anti-clockwise so that the translated jaw link (c 2 ) tends to close the pharyngeal jaw (the biting action). In this regard, it was decided that the most suitable mechanism type utilizing kinematics and geometrical constraints for the rigid connection of these two links is the anti-parallel four-bar mechanism. In light of the given information, both mechanisms will be synthesized by analytical methods.

Kinematic Synthesis of the Four-Bar Mechanism and the Pharyngeal Jaw
In Section 2, the motion profile of the pharyngeal jaw was explained in detail. The motion of this jaw can be examined as a combination of pure rotation and pure translation. In this section, the four-bar mechanism which presents pure rotation of this jaw will be examined.
The rotational motion of the pharyngeal jaw is constrained by the natural characteristics of the moray itself. In other words, output angle (α) and its time steps are bio-constraints ( Figure 8a). In order to ensure these characteristics, a function must be described between the input angle (θ 4 ), which comes from the slider, and the output angle (α). As explained before, such necessity can be solved by function generation synthesis and in this part of the study, this method is used.
In this configuration of the four-bar mechanism, the distance between the grounded joints is assumed as a biomimetic constraint and defined by the actual size of the slider of the pharyngeal jaw. This quantity affects the construction parameters (unknown link lengths of the four-bar mechanism) linearly. Therefore, it can be concluded that this parameter (d 2 in Figure 8a) is a scale factor for the construction parameters and the pharyngeal jaw and it is a predefined and known parameter.
Function generation synthesis starts with the x and y components of closed-loop equations which can be visualized in Figure 8a and revealed in Equations (1) and (2): In order to obtain a relationship between α and θ 4 , θ 5 should be eliminated from the equations by adding the squares of Equations (3) and (4): The resultant equation is shaped as follows to form the objective function as in Equation (5): With the necessary arrangements, the objective function in the form of a polynomial is shaped with the construction parameters (unknown link lengths) and variables: Equation (6) can be represented as the polynomial form as in Equation (7), in which the parameters were as in Equations (8) and (9): where, and, In order to solve Equation (6), the values of terms θ 4 and α should be given as precision points. As there exist only three unknowns in Equation (7) (P 0 , P 1 , P 2 ) three equations, collected in Equation (10), are needed to reach the solution by three sets of precision points α i = f(θ 4i ), i = 1, 2, 3: With three precision points, it is possible to obtain the construction parameters. The matrix form of Equation (10) is presented in Equation (11): Construction parameters are as in Equation (12): As a result, in order to define a four-bar mechanism that needs to mimic the biting action of the pharyngeal jaw of a moray eel, input-output relations at three different points and the distance between ground points (d 2 ) must be given to this system of equations.

Kinematic Synthesis of the Offset Slider-Crank Mechanism of the Pharyngeal Jaw
While the four-bar mechanism of the pharyngeal jaw is present for the biting action, another task of this jaw set is to travel from the throat to the outer oral jaw or vice-versa. As mentioned before, the motion of the pharyngeal jaw can be examined with two separate motions: pure rotation and pure translation. Translational motion is carried out by an offset slider-crank mechanism. This mechanism not only achieves the traveling (along with s 1 in Figure 8b) but also carries and actuates the four-bar mechanism by the rotational motion of the rigidly connected c 1 and a 2 links (Figure 9). Here, two dead-end points are considered so as to synthesize the slider-crank mechanism ( Figure 11).
Biomimetics 2022, 7, x FOR PEER REVIEW 12 of 32 offset slider-crank mechanism. This mechanism not only achieves the traveling (along with s1 in Figure 8b) but also carries and actuates the four-bar mechanism by the rotational motion of the rigidly connected c1 and a2 links ( Figure 9). Here, two dead-end points are considered so as to synthesize the slider-crank mechanism ( Figure 11). Figure 11. The Offset slider-crank mechanism poses at two dead-end points.
From the trigonometric relations, construction parameters of the offset slider-crank mechanism are derived as following in Equation (13) Figure 11. The Offset slider-crank mechanism poses at two dead-end points.
From the trigonometric relations, construction parameters of the offset slider-crank mechanism are derived as following in Equation (13): In this synthesis, procedure ∆s 1-max , θ 4-1-1 , and θ 4-1-2 are given to system equations according to the requirements of the jaw mechanism, which will be discussed in Section 4. Eventually, these requirements define the shape and limit of the offset slider-crank mechanism. Then, for constructing the overall shape of the pharyngeal jaw, biomimetic and functional design constraints for a combination of both mechanisms are considered in numerical works.

The Kinematic Synthesis of the Oral Jaw
The motion of the oral jaw is described as a pure rotation since its duty is to open and close the oral jaw within pre-defined time instants in order to catch the prey ( Figure 10). The biomimetic motion is achieved by using the inverted slider-crank mechanism, where the sliding input is given to the system by a synthesized motion profile. Here, this mechanism is convenient to fulfill the requirements and increase the design variations in the work. Moreover, it is easy to place the input motion profile by means of a mechanical point of view.
Since the actuator on the system has a rotary output, a motion profile of the sliding motion is obtained by creating a motion profile with respect to requirements at each time instant.
The Kinematic Synthesis of the Inverted Slider-Crank Mechanism of the Oral Jaw The oral jaw is the outer part of the mechanism, where the moray eel catches the prey with only a pure rotation motion as seen in Figure 6. In the previous sections, bio-constraints were defined for the motion of the oral jaw.
The biting action of the oral jaw is based on a rotation around a hinge, where the jaw achieves a closing and opening action. However, this rotation should work together sequentially with the pharyngeal jaw. Therefore, in the structural synthesis study, it was concluded that transferring the continuous rotational input motion to the slider-crank mechanism, so as to obtain a predefined oral jaw motion, was bound to the cam profile. In its simplest way, the system is built as an inverted slider-crank mechanism as in Figure 10.
As seen in Figure 10, there are some parameters that affect the efficiency of the resultant mechanism, thus, these parameters have some flexibility. For instance, we cannot conclude that the angle φ should be some sharp value due to the fact that it might be needing a different value to achieve the higher force transmission for biting, so the parameters φ and β are given to the system equations in the logical boundaries, which will be discussed in Section 4.2. Here, the analytic kinematic synthesis methodology is carried out as if those parameters are known and the variables are s 0 and ϕ, which are the input stroke profile and the oral jaw angle, respectively.
Function generation kinematic synthesis starts with x and y components of closed-loop equations which can be visualized in Figure 10 and revealed in Equations (14) and (15): In order to obtain a relation between s 0 and ϕ, γ should be eliminated from the equations by adding the squares of Equations (16) and (17): With the necessary arrangements, objective function in the form of a polynomial is shaped with the construction parameters (unknown link lengths) and variables: Equation (18) can be represented as the polynomial form as in Equation (19), in which parameters were as in Equations (20) and (21): where, and, In order to solve Equation (18), values of terms ϕ and s 0 should be given as precision points. As there exist only two unknowns in Equation (19) (P 0 , P 1 ) two equations, collected in Equation (22), are needed to reach the solution by two sets of precision points ϕ i = f(s 0i ), i = 1, 2: With two precision points, it is possible to get the construction parameters. The matrix form of Equation (22) is presented in Equation (23): Construction parameters are as in Equation (24): As a result, in order to define an inverted slider-crank mechanism that needs to mimic the biting action of the oral jaw of a moray eel, the input-output relation at two different points must be given to the system of equations.

Parameter Optimization
The main goal of this design is to create a single DOF mechanism that mimics the movements of moray eel jaws during hunting, which has three outputs in the view of kinematic representation. In general, after utilizing the kinematic synthesis procedures, the output position or orientation is optimized with various different techniques in order to reach the most accurate mechanism according to need. However, in this overall mechanism, there is no strict position or orientation needs at any instant, instead, relatively greater intervals are required such as at the specific instant minimum possible value of the jaw angle. For this mechanism, rather than positions and orientations, torque delivering at the time of when biting occurs is more crucial. Thus, it is focused on maximizing torque delivering at specific biting instants.
If the kinematic motion of the overall mechanism is examined, expected velocities of the mechanism parts are relatively slow. Together with this reason, although there are several methods available to analyze the force/torque relationship between the input and the output, instead of making a complex dynamic analysis, virtual work methodology is utilized without including the effects of the masses and inertias of the small links, which is negligible compared to the resultant biting forces.
In order to utilize the virtual work methodology in overall mechanism, position and velocity output depending on the crank and the prismatic joint input should be determined at every instant, thus, kinematic analysis in terms of position and velocity is a must.

Kinematic Analysis
After applying kinematic synthesis procedures to each individual mechanisms, kinematic analysis should be applied for simulating motions, utilizing virtual work method, and checking if the resultant mechanisms fit the constraints at different input angles.

Kinematic Position Analysis of the Four-Bar Mechanism of the Pharyngeal Jaw
As in the synthesis case, position analysis also starts with the closed-loop equations as in Equations (1) and (2) and dismiss the angle θ 5 as in Equations (3) and (4), eventually shaped Equations (25) and (26): From a half tangent formula, the procedure continues from Equation (27) to Equation (32): With the following equations, we can determine the output angles at each input instants as: Time dependent varying values of input angle θ 4 affect the outputs of α and θ 5 as in Equations (33) and (34), which are used in a simulation process and in dynamic relations.

Kinematic Position Analysis of the Offset Slider-Crank Mechanism of the Pharyngeal Jaw
Position analysis starts with the closed loop equations as in Equations (35) and (36) and dismiss the angle θ 3 as in Equations (37) and (38) and continues with Equations (39) and (40): By letting discriminant as in Equation (41): Time dependent varying values of input angle θ 2 affect outputs of s 1 and θ 3 as in Equations (42) and (43), which are used in a simulation process and in dynamic relations.

Kinematic Position Analysis of Inverted Slider-Crank Mechanism of Oral Jaw
As in the synthesis case, position analysis also starts with the closed-loop equations as in Equations (14) and (15) and dismiss the angle γ as in Equations (16) and (17) From a half tangent formula, the procedure continues from Equations (46) to (51): With the following equations, we can determine the output angles at each input instants as: Time dependent varying values of input profile s 0 affect outputs of ϕ and γ as in Equations (52) and (53), which are used in a simulation process and in dynamic relations.

Kinematic Velocity Analysis of the Overall Pharyngeal Jaw Mechanism
In order to determine the velocities at each instant, kinematic closed-loop equations are used as follows: Loop 1 in Equations (54) and (55): (56) and (57):

Loop 2 in Equations
With the following Equation (69), we can determine the rotational velocity of the pharyngeal jaw at each input instant: Time dependent varying position and velocity values of input angle θ 2 affect rotational velocity of the pharyngeal jaw ( . α) in Equation (69), which will be used in dynamic analysis.

Kinematic Velocity Analysis of the Overall Oral Jaw Mechanism
In order to determine the velocities at each instant, kinematic closed-loop equations, Equations (70) and (71), were used as follows from Equations (72) to (74): .
Time dependent varying position and velocity values of input profile s 0 affect rotational velocity of the oral jaw ( . ϕ) in Equation (74), which will be used in dynamic analysis.

Input Sets of Mechanisms
Since the main concern here is to mimic the desired action of jaws, kinematic position error is not critical, therefore, it is not needed to be so precise in terms of position, and, therefore, the motion characteristic is more important. By concerning the anti-parallel four-bar mechanism that slides as a whole, which makes the biting action of the pharyngeal jaw, the output jaw angle should have specific values at certain moments such as when biting takes place or at a time when the jaw is fully opened. At the same time, the oral jaw and the slider should also be at predefined positions, which were given as bio-constraints in Section 2. However, corresponding input angles cannot be specified sharply in the first place, since input has the flexibility of being in a different range. For instance, the input range of the four-bar mechanism (θ 4 ) of the pharyngeal jaw can be 30 • as well as 90 • , which directly affects the slider-crank mechanism. Here, what eventually matters is specified output angle.
For these reasons, it was decided that we should create an algorithm by using a computing program in order to utilize different kinematic synthesis by changing the input ranges as well as the additional parameters. These additional parameters are α offset for the four-bar mechanism of the pharyngeal jaw (Figure 8a), θ 4-2 for the slider-crank mechanism of the pharyngeal jaw (Figure 9), and φ and β for the inverted slider-crank mechanism of the oral jaw ( Figure 10). Firstly, input sets of all mechanisms should be created in a logical manner.

Input Sets of the Four-Bar Mechanism of the Pharyngeal Jaw
By considering the structural synthesis of combined mechanisms, while the slidercrank mechanism part carries the whole pharyngeal jaw by translating it out to the end of the mouth, it has a motion on +x direction (towards the end of the mouth) as the angle between +x and link c 1 gets bigger ( Figure 9). As the jaw slides out to the end of the mouth, the biting action should take place, which means that the jaw angle (α-α offset ) should get smaller, close to zero degrees. However, the initial value of the input angle of the sliding anti-parallel four-bar mechanism (θ 4 ) can be extremely flexible, therefore, its range of motion is limited to between 270 • and 360 • . Moreover, the difference between the initial and the final angle of θ 4 is chosen between 90 • and 30 • . With this information, an array of set of initial θ 4 and an array of set of final θ 4 are created with a step size of 5 • . As mentioned in the kinematic synthesis part, the sliding four-bar mechanism needs three precision points to be constructed, where the remaining θ 4 is defined with linear spacing. A numerical example of the first 8 index of θ 4 array can be seen in Table 2. The output set of the sliding four-bar mechanism should not change whereas the input set changes. When the input is on the state of initial θ 4 , the pharyngeal jaw should be open, which is considered to be 45 • , with the jaw angle as an open state. For the next step, when the input is on the state of mean θ 4 , the pharyngeal jaw should be approximately half open, which is considered to be 20 • , with the jaw angle as a medium state. Finally, when the input is on the state of final θ 4 , the pharyngeal jaw should be fully closed (biting state), which is considered to be 3 • , with the jaw angle as a closed state. However, as seen from the mechanism synthesis (Figure 9), there is also a platform at the end of the considered jaw, which directly affects the mechanism synthesis results, thus, its shape is flexible as well. For this reason, the angle in the platform (α offset ) is also considered to be variable, therefore, multiple αoffset angle values are created with the offset from 60 • to 160 • with a step size of 1 • .
For each αoffset and θ 4 difference values, 1911 different input sets are created and all of them are ready to be used in the synthesis study of the four-bar mechanism on the slider, which is eventually used in the biting action.

Input Sets of the Slider-Crank Mechanism of the Pharyngeal Jaw
In order to size the overall mechanism that mimics the average moray eel, the slider travel is chosen as 140 mm, which directly carries the pharyngeal jaw mechanism from throat to the end of the mouth. As seen in Figure 9, there is a fixed offset angle (θ 4-2 ), which creates a platform on the mechanism. This fixed angle creates flexibility on the design, which should not be chosen randomly because of the effect on the overall system. Therefore, in the algorithm, another loop is carried out in which the slider offset angle (θ 4-2 ) is changed from 91 • to 225 • with a step size of 1 • . Together with the previous selections, the input set increases to 9164 different combinations.
For each usable sliding four-bar mechanism, which are determined in previous section, with initial θ 4 angle with specified θ 4-2 offset, analytical kinematic synthesis with two positions is carried out for offset slider-crank mechanism. As a result, if the construction parameters (link lengths of the offset slider-crank mechanism) is real and positive, it is possible to conclude that offset the slider-crank mechanism is suitable to work with the sliding four-bar mechanism together without any constraint.

Input Sets of the Inverted Slider-Crank Mechanism of the Oral Jaw
By considering Figure 10, input of the oral jaw is actually a prismatic joint (s 0 ), whose motion can be given to the system as an input depending on the crank position by using a cam profile that drives the instantaneous output jaw angle, which is rotated to open or close oral jaw itself at every instant in one full turn of the crank. Eventually, the cam profile is actually extracted as a need to give the required motion for the prismatic joint in the oral jaw. Here, the flexibilities are on the direction of the input profile (φ) and output jaw platform (β), thus, these parameters should be multiple in logical boundaries. Moreover, the range and initial value of the input profile s 0 affects the shape and the force transmission of the resultant mechanism, as in the case of the input angle of the sliding anti-parallel four-bar mechanism (Section 4.2.1).
In logical boundaries, angle φ is from −50 • to 50 • , while angle β should change from 60 • to 150 • with 5 • step sizes. The input profile should also be within predefined values, which are from 20 mm to 50 mm with an initial value of 20 mm. For the different values of each parameter, a kinematic synthesis procedure was applied.

Multiple Synthesis Algorithm
Having flexibility on the parameters is a challenge for the design procedure due to the fact that it affects resultant actions dramatically. For this reason, the input-output relations were created within the logical boundaries so that the multiple results were concluded. Among these, the best result was chosen to be the one that has the highest output to input torque ratio at critical moments, which is at the fully closed position in order to increase the biting force, so that the smallest actuator can serve as the best torque transformer. Moreover, it is important to note that there are some functional kinematic constraints that depend on the motion characteristic to be followed by the overall mechanism in order to mimic the motion.
Defining input sets as discussed in previous sections, the flow of the multiple synthesis algorithm is based on a predefined logic. In this regard, a novel algorithm which is based on the kinematic synthesis is created in order to define the input-output set and the starting positions. The created algorithm, which worked in a MATLAB environment, starts with defining the input-output and the jaw angle sets of the sliding four-bar mechanism according to the fundamental needs of the biting action of the pharyngeal jaw. The needed equations for the synthesis study of the sliding four-bar mechanism are derived in Equation (12). As an input set changes with a predefined logic, a loop for determining the construction parameters (four-bar mechanism link lengths) is executed. For each input set, analytically calculated resultant link lengths are compared to zero, since a negative link length is not possible. If all link lengths are larger than zero, then it is possible to conclude that the resultant four-bar mechanism is usable without considering main constraints and design goals. As a result of the multiple synthesis for each input sets, besides link lengths, initial and final input angles (θ 4 ) and a jaw offset angle (α offset ) at the action of biting (slider is at the positive dead end) that meet predefined criteria are saved in the memory to be used later (Figure 12). link length is not possible. If all link lengths are larger than zero, then it is po conclude that the resultant four-bar mechanism is usable without considering m straints and design goals. As a result of the multiple synthesis for each input sets link lengths, initial and final input angles (θ4) and a jaw offset angle (αoffset) at the biting (slider is at the positive dead end) that meet predefined criteria are sav memory to be used later ( Figure 12). The offset slider-crank mechanism also has a predefined fixed parameter to in the synthesizing procedure (θ4-2). However, according to the requirements in this parameter has flexibility as well. For this reason, an interval is defined, wh are used in synthesizing the offset slider-crank mechanism for each step. By co the successful outputs of the four-bar mechanisms and a flexible parameter, the procedure is applied for each step and the ones that meet predefined criteria are the memory as in the previous case ( Figure 13). The offset slider-crank mechanism also has a predefined fixed parameter to be used in the synthesizing procedure (θ 4-2 ). However, according to the requirements in output, this parameter has flexibility as well. For this reason, an interval is defined, where they are used in synthesizing the offset slider-crank mechanism for each step. By combining the successful outputs of the four-bar mechanisms and a flexible parameter, the synthesis procedure is applied for each step and the ones that meet predefined criteria are saved in the memory as in the previous case ( Figure 13). At this step, it is reached some successful four-bar mechanisms and offset slidercrank mechanisms, which meet their criteria separately. The combination of these mechanisms should meet newer combined criteria, which are the main requirements for mimicking the pharyngeal jaw. The whole system should achieve the specific motion functional constraints that are formed by following:

•
At the negative dead end of the slider-crank mechanism, the pharyngeal jaw should be open. In order to achieve the above specific motion combinations, motion analysis of each created mechanism was carried out for a full revolution of the main input crank. At specified instants, the angles were compared to where they should occur. The first instant is the moment when slider is on the middle of travel from positive end to negative end (pulling the prey action), at which the pharyngeal jaw angle should be smaller than 5°. This means that the sliding four bar should be close to the pharyngeal jaw while being pulled into the esophagus. The second instant is the moment when slider is again in the middle of travel from the negative end to the positive end (catching the prey action), at which the pharyngeal jaw angle should be greater than 25°. This means that the sliding four bar should open the pharyngeal jaw while being pushed to the mouth's end. This algorithm eliminates impractical mechanisms which are not needed to be optimized for the purpose of mimicking the moray eel's pharyngeal jaw (Figure 14). At this step, it is reached some successful four-bar mechanisms and offset slider-crank mechanisms, which meet their criteria separately. The combination of these mechanisms should meet newer combined criteria, which are the main requirements for mimicking the pharyngeal jaw. The whole system should achieve the specific motion functional constraints that are formed by following:

•
At the negative dead end of the slider-crank mechanism, the pharyngeal jaw should be open. In order to achieve the above specific motion combinations, motion analysis of each created mechanism was carried out for a full revolution of the main input crank. At specified instants, the angles were compared to where they should occur. The first instant is the moment when slider is on the middle of travel from positive end to negative end (pulling the prey action), at which the pharyngeal jaw angle should be smaller than 5 • . This means that the sliding four bar should be close to the pharyngeal jaw while being pulled into the esophagus. The second instant is the moment when slider is again in the middle of travel from the negative end to the positive end (catching the prey action), at which the pharyngeal jaw angle should be greater than 25 • . This means that the sliding four bar should open the pharyngeal jaw while being pushed to the mouth's end. This algorithm eliminates impractical mechanisms which are not needed to be optimized for the purpose of mimicking the moray eel's pharyngeal jaw (Figure 14).

Figure 14.
Algorithm for combining the four-bar mechanism and the offset slider-crank mechanism to get a functional pharyngeal jaw. The next step is to design the oral jaw by considering the design constraints. As discussed in the synthesis procedure, two different precision points are required in order to define the mechanism. Here, the two dead ends of the input slider profile are given as flexible and the corresponding outputs of the main jaw mechanism are fully opened and fully closed. Motion in between these two ends is extracted according to the one cycle movement in one full turn, where it is fulfilled by the cam profile.
After utilizing the kinematic synthesis procedure for the oral jaw mechanism, according to the functional timing of the overall movement cycle, medium movements are extracted as defined by the cam profile. Functional timing of the overall movement cycle was approximately predicted in Figure 4. The resultant mechanisms that meet the criteria with created cam profiles were saved in the memory along with their important parameters ( Figure 15). Figure 14. Algorithm for combining the four-bar mechanism and the offset slider-crank mechanism to get a functional pharyngeal jaw. The next step is to design the oral jaw by considering the design constraints. As discussed in the synthesis procedure, two different precision points are required in order to define the mechanism. Here, the two dead ends of the input slider profile are given as flexible and the corresponding outputs of the main jaw mechanism are fully opened and fully closed. Motion in between these two ends is extracted according to the one cycle movement in one full turn, where it is fulfilled by the cam profile.
After utilizing the kinematic synthesis procedure for the oral jaw mechanism, according to the functional timing of the overall movement cycle, medium movements are extracted as defined by the cam profile. Functional timing of the overall movement cycle was approximately predicted in Figure 4. The resultant mechanisms that meet the criteria with created cam profiles were saved in the memory along with their important parameters ( Figure 15). In order to have an overall mechanism that completely mimics the motion of the two jaws of the moray eel, successful mechanisms should achieve the defined criteria. Eventually, it is possible to reach many distinct mechanisms, where optimization is needed.

Optimization of Mechanisms
With an excessive amount of resultant possible mechanisms, the most effective one should be taken into account. The effectiveness depends on certain parameters which affect overall efficiency, quality, and usability. In this design, the first criteria is chosen as force/torque transmitting from the actuator to both jaws, the oral and the pharyngeal jaw, when the biting action occurs, separately. From the input link to the output jaws, there are various links and connections to transform the desired motions. As mentioned in earlier sections, relative to required forces in the outputs, all of the masses and inertias can be neglected in order to simplify the force/torque analysis. In this work, a virtual work method is used for calculating the relationship between the input and the output torques.

The Optimization of the Pharyngeal Jaw
For the pharyngeal jaw, the first thing one must consider is the torque transmission in the anti-parallel four-bar mechanism, which plays a vital role in biting prey and pulling it through the esophagus. For this reason, this biting torque should be as high as possible Figure 15. Algorithm for designing the inverse slider-crank mechanism at the oral jaw.
In order to have an overall mechanism that completely mimics the motion of the two jaws of the moray eel, successful mechanisms should achieve the defined criteria. Eventually, it is possible to reach many distinct mechanisms, where optimization is needed.

Optimization of Mechanisms
With an excessive amount of resultant possible mechanisms, the most effective one should be taken into account. The effectiveness depends on certain parameters which affect overall efficiency, quality, and usability. In this design, the first criteria is chosen as force/torque transmitting from the actuator to both jaws, the oral and the pharyngeal jaw, when the biting action occurs, separately. From the input link to the output jaws, there are various links and connections to transform the desired motions. As mentioned in earlier sections, relative to required forces in the outputs, all of the masses and inertias can be neglected in order to simplify the force/torque analysis. In this work, a virtual work method is used for calculating the relationship between the input and the output torques.

The Optimization of the Pharyngeal Jaw
For the pharyngeal jaw, the first thing one must consider is the torque transmission in the anti-parallel four-bar mechanism, which plays a vital role in biting prey and pulling it through the esophagus. For this reason, this biting torque should be as high as possible in order to avoid the prey fleeing. With a unit applied, the torque from the input actuator, and the torque at the pharyngeal jaw can be calculated as in Equations (75) and (76): By using the result from the velocity analysis of the pharyngeal jaw, kinematic Equation (76) is combined with Equation (69), that shapes Equation (77), in order to determine torque transmission: Since the main objective is to keep the τ pharyngeal jaw to a maximum at the time when the biting occurs, this equation should be applied to all resultant possible mechanisms with a unit of τ 0 .

Optimization of the Oral Jaw
For the oral jaw, the main actuator creates the required torque that is then transformed to a linear force along with the slider. With the neglected masses and inertias compared to the desirable forces and torques, the torque transmission from the actuator to the teeth can be determined with a virtual work method.
The force on the slider depends on the cam profile as seen in Figure 16 with the unit input torque (τ 0 ) and can be calculated with Equation (78): Biomimetics 2022, 7, x FOR PEER REVIEW 26 of 32 in order to avoid the prey fleeing. With a unit applied, the torque from the input actuator, and the torque at the pharyngeal jaw can be calculated as in Equations (75) Since the main objective is to keep the τpharyngeal jaw to a maximum at the time when the biting occurs, this equation should be applied to all resultant possible mechanisms with a unit of τ0.

Optimization of the Oral Jaw
For the oral jaw, the main actuator creates the required torque that is then transformed to a linear force along with the slider. With the neglected masses and inertias compared to the desirable forces and torques, the torque transmission from the actuator to the teeth can be determined with a virtual work method.
The force on the slider depends on the cam profile as seen in Figure 16 with the unit input torque (τ0) and can be calculated with Equation (78): where Γ can be determined with vector operations. Tangential force drives the slider linearly, eventually it results in squeezing the oral jaw of the moray eel as in Equations (79) and (80): The main objective is to keep the τ oral jaw to a maximum at the time that the biting occurs, this equation should be applied to all resultant possible mechanisms with a unit of τ 0 .

Results and Discussions
In order to reach the ultimate goal of this work, all of the needed algorithms and analytic equations are derived in the previous sections. Since the procedure consists of different methodologies, it can be summarized with its application flow as follows: • Define the input sets of the four-bar mechanism of the pharyngeal jaw (θ 4 ) • Define the input sets of the four-bar mechanism output platform angle (α offset ) • Define the four-bar mechanism scale according to the requirement for size (d 2 ) • Utilize the Equations (11) and (12) for each input set in order to get the construction parameters of the four-bar mechanism of the pharyngeal jaw (a 2 , b 2 , c 2 ) • Define the pharyngeal jaw's sliding travel (∆s 1 ) • Define the input sets of the slider offset angle (θ 4-2 ) • Utilize Equation (13) for each input set in order to get the construction parameters of the offset slider-crank mechanism of the pharyngeal jaw (a 1 , b 1 , c 1 ) • Combine both pharyngeal jaw mechanisms that follow the functional constraints and then identify the one with the highest torque transmitting ratio at the moment of the pharyngeal jaw biting • Define the input sets of the inverted slider-crank mechanism of the oral jaw (s 0 ) • Define the input sets of the direction of the slider profile (φ) • Define the input sets of the output oral jaw platform (β) • Define the inverted slider-crank mechanism scale according to the requirement for the size (e 0 ) • Utilize Equations (23) and (24) for each input set in order to get the construction parameters of the inverted slider-crank mechanism of the oral jaw (c 0 , d 0 , and the resultant s 0 interval) • Combine all of the possible functional oral jaw mechanisms with the specified pharyngeal jaw and then identify the one with the highest torque transmitting ratio at the moment of oral jaw biting After utilizing the procedure above, the resultant mechanism is determined in terms of millimeters and degrees as shown in Table 3. In terms of kinematics, the resultant mechanism gives an output curve in one full rotation of the main actuator crank as seen in Figure 17. In Figure 4, the desired biting sequence is illustrated. Timing in the orientations of the jaws in Figure 4 was predicted in order to ensure harmony between the two jaws. A comparison between Figures 4 and 7 show that there exists a difference in the angles of the jaws. The main reason for this difference is one of the design goals of the optimization study. The proposed mechanical structure was generated with both kinematic and dynamic constraints. The constraint, which causes profile shifting in the biting sequence, is the maximization of the force/torque rate. In order to maximize the transferred torque through the pharyngeal jaw, it is necessary to manipulate the derivative of the position equation of the pharyngeal jaw. Therefore, optimization of the torque delivery through the pharyngeal jaw is dependent on the velocity equations (Equation (76)). This relationship between the optimization study and the dynamic behavior of the pharyngeal jaw causes the difference in the biting sequence of this jaw. difference is one of the design goals of the optimization study. The proposed mechanical structure was generated with both kinematic and dynamic constraints. The constraint, which causes profile shifting in the biting sequence, is the maximization of the force/torque rate. In order to maximize the transferred torque through the pharyngeal jaw, it is necessary to manipulate the derivative of the position equation of the pharyngeal jaw. Therefore, optimization of the torque delivery through the pharyngeal jaw is dependent on the velocity equations (Equation (76)). This relationship between the optimization study and the dynamic behavior of the pharyngeal jaw causes the difference in the biting sequence of this jaw. The resultant mechanism is designed roughly in SolidWorks, as seen in Figure 18, in order to simulate the resultant movement as well as perform the needed dynamic analysis at critical instants with a different computer program. In a MATLAB environment, all of the optimization and kinematic synthesis procedures are utilized. The two critical biting instants are determined in terms of variables. These variables for the pharyngeal jaw are shown in Table 4. At biting instance, the torque ratio is revealed as being 7.28 by using Equation (77). At this critical instant, the kinematic The resultant mechanism is designed roughly in SolidWorks, as seen in Figure 18, in order to simulate the resultant movement as well as perform the needed dynamic analysis at critical instants with a different computer program. structure was generated with both kinematic and dynamic constraints. The constraint, which causes profile shifting in the biting sequence, is the maximization of the force/torque rate. In order to maximize the transferred torque through the pharyngeal jaw, it is necessary to manipulate the derivative of the position equation of the pharyngeal jaw. Therefore, optimization of the torque delivery through the pharyngeal jaw is dependent on the velocity equations (Equation (76)). This relationship between the optimization study and the dynamic behavior of the pharyngeal jaw causes the difference in the biting sequence of this jaw. The resultant mechanism is designed roughly in SolidWorks, as seen in Figure 18, in order to simulate the resultant movement as well as perform the needed dynamic analysis at critical instants with a different computer program. In a MATLAB environment, all of the optimization and kinematic synthesis procedures are utilized. The two critical biting instants are determined in terms of variables. These variables for the pharyngeal jaw are shown in Table 4. At biting instance, the torque ratio is revealed as being 7.28 by using Equation (77). At this critical instant, the kinematic In a MATLAB environment, all of the optimization and kinematic synthesis procedures are utilized. The two critical biting instants are determined in terms of variables. These variables for the pharyngeal jaw are shown in Table 4. At biting instance, the torque ratio is revealed as being 7.28 by using Equation (77). At this critical instant, the kinematic design in SolidWorks can also be seen in Figure 19a with a corresponding position of moray illustration in Figure 19b. design in SolidWorks can also be seen in Figure 19a with a corresponding position of moray illustration in Figure 19b.  Figure 19. The pharyngeal jaw at biting instant: (a) as the mechanism; (b) as the moray illustration.
In this section, the control of the result of the mathematical model is carried out by using the Multi Body Dynamics (MBD) simulation method. In this method, the simulation CAD model was created by using construction parameters and variables. The model does not include manufacturing details or components but geometrical details for analysis.
The MBD simulation model of the proposed mechanism was integrated into simulation environments with boundary conditions and initial conditions. These conditions ensure the exact data of the transferred torque through the jaw at that instant time. In order to obtain this data, the joint of the pharyngeal jaw was locked at the biting time and unit vector torque 1 [N/mm] was applied through the input link. As a result of the simulation, the reaction torque was measured on the joint of the jaw. The simulation results as a reaction torque are shown in Figure 20. Due to the applied torque magnitude, the output torque is equal to the torque conversion ratio, derived from Equation (77). It was determined at 7.48, which deviates from the analytic result by 2.75%. Similarly, other critical biting instants for the oral jaw were determined in terms of variables as well and can be seen in Table 5 with its SolidWorks kinematic design in Figure  21a with corresponding position of moray illustration in Figure 21b.  In this section, the control of the result of the mathematical model is carried out by using the Multi Body Dynamics (MBD) simulation method. In this method, the simulation CAD model was created by using construction parameters and variables. The model does not include manufacturing details or components but geometrical details for analysis.
The MBD simulation model of the proposed mechanism was integrated into simulation environments with boundary conditions and initial conditions. These conditions ensure the exact data of the transferred torque through the jaw at that instant time. In order to obtain this data, the joint of the pharyngeal jaw was locked at the biting time and unit vector torque 1 [N/mm] was applied through the input link. As a result of the simulation, the reaction torque was measured on the joint of the jaw. The simulation results as a reaction torque are shown in Figure 20. Due to the applied torque magnitude, the output torque is equal to the torque conversion ratio, derived from Equation (77). It was determined at 7.48, which deviates from the analytic result by 2.75%. moray illustration in Figure 19b.  Figure 19. The pharyngeal jaw at biting instant: (a) as the mechanism; (b) as the moray illustration.
In this section, the control of the result of the mathematical model is carried out by using the Multi Body Dynamics (MBD) simulation method. In this method, the simulation CAD model was created by using construction parameters and variables. The model does not include manufacturing details or components but geometrical details for analysis.
The MBD simulation model of the proposed mechanism was integrated into simulation environments with boundary conditions and initial conditions. These conditions ensure the exact data of the transferred torque through the jaw at that instant time. In order to obtain this data, the joint of the pharyngeal jaw was locked at the biting time and unit vector torque 1 [N/mm] was applied through the input link. As a result of the simulation, the reaction torque was measured on the joint of the jaw. The simulation results as a reaction torque are shown in Figure 20. Due to the applied torque magnitude, the output torque is equal to the torque conversion ratio, derived from Equation (77). It was determined at 7.48, which deviates from the analytic result by 2.75%. Similarly, other critical biting instants for the oral jaw were determined in terms of variables as well and can be seen in Table 5 with its SolidWorks kinematic design in Figure  21a with corresponding position of moray illustration in Figure 21b.  Similarly, other critical biting instants for the oral jaw were determined in terms of variables as well and can be seen in Table 5 with its SolidWorks kinematic design in Figure 21a with corresponding position of moray illustration in Figure 21b.  The powerful nature behind the biomimetic designs can inspire similar designs utilized procedure stems from a biomimetic design with the aim of reducing the task D This procedure can be applied to different kinds of tasks that are mechanically conne to each other. In this work, the optimization parameter was chosen as the torque tran ring ratio. By following similar multiple analytical kinematic synthesis for input-ou flexibility, not only is the torque transferring ratio chosen as the optimization param but it can also be a kinematic position or an orientation error, by using existing optim tion methodologies. The main idea here is to have a procedure that can be follow order to get a combined mechanism that gives multiple outputs with a single degr freedom. For similar repeating tasks with such designs, it is possible to reduce the exi control and timing complexity and increase the safety. By considering the simula with different tools, it is shown that the possibility of such designs can be adapted e and can verify the techniques that are used. This kind of procedure can be applied cially to automation systems that repeat the movements in which more DOF was alr used and where algorithms are utilized. With a quicker kinematic design and embe optimization, the overall cost can be reduced due to there being no need for more a tors, sensors, and complex control algorithms.

Conclusion and Future Works
In this work, a mechanism that mimics the double jaw apparatus of the mora was designed and optimized with a lower DOF than task DOF. A mechanism with a s actuator point drives the resultant mechanism and forms the kinematic movements. I kinematic synthesis procedure, the bio-constraints were defined and, depending on t constraints, some construction parameters were forced into logical intervals. A struc design procedure was utilized in order to reach the ultimate goal and the dimensio part was carried out with multiple kinematic synthesis algorithms due to the flexibil terms of inputs and constraints. The procedure started with the pharyngeal jaw, w slides and rotates simultaneously. Depending on the resultant construction paramete the pharyngeal jaw, a similar approach was utilized for the oral jaw. As a result o algorithm, many different mechanisms were created, among which the most torqu livering device at critical instants was chosen. The procedure was explained in deta cluding bio-constraints, structural synthesis, kinematic synthesis, and optimization. I end, the resultant mechanism was roughly designed in SolidWorks in order to visu the movements.
In future work, the mechanical design will be improved and the manufacturing cess will take place. The plan is to carry out the validation procedure in terms of kin and kinematics. A kinematic synthesis and optimization procedure is handled with a ware that is developed in a MATLAB environment. However, since there are man rameters that should be specified by the user, it is not easy to hard-code using the sy The powerful nature behind the biomimetic designs can inspire similar designs. The utilized procedure stems from a biomimetic design with the aim of reducing the task DOF. This procedure can be applied to different kinds of tasks that are mechanically connected to each other. In this work, the optimization parameter was chosen as the torque transferring ratio. By following similar multiple analytical kinematic synthesis for input-output flexibility, not only is the torque transferring ratio chosen as the optimization parameter but it can also be a kinematic position or an orientation error, by using existing optimization methodologies. The main idea here is to have a procedure that can be followed in order to get a combined mechanism that gives multiple outputs with a single degree of freedom. For similar repeating tasks with such designs, it is possible to reduce the existing control and timing complexity and increase the safety. By considering the simulations with different tools, it is shown that the possibility of such designs can be adapted easily and can verify the techniques that are used. This kind of procedure can be applied especially to automation systems that repeat the movements in which more DOF was already used and where algorithms are utilized. With a quicker kinematic design and embedded optimization, the overall cost can be reduced due to there being no need for more actuators, sensors, and complex control algorithms.

Conclusions and Future Works
In this work, a mechanism that mimics the double jaw apparatus of the moray eel was designed and optimized with a lower DOF than task DOF. A mechanism with a single actuator point drives the resultant mechanism and forms the kinematic movements. In the kinematic synthesis procedure, the bio-constraints were defined and, depending on those constraints, some construction parameters were forced into logical intervals. A structural design procedure was utilized in order to reach the ultimate goal and the dimensioning part was carried out with multiple kinematic synthesis algorithms due to the flexibility in terms of inputs and constraints. The procedure started with the pharyngeal jaw, which slides and rotates simultaneously. Depending on the resultant construction parameters of the pharyngeal jaw, a similar approach was utilized for the oral jaw. As a result of the algorithm, many different mechanisms were created, among which the most torque delivering device at critical instants was chosen. The procedure was explained in detail including bio-constraints, structural synthesis, kinematic synthesis, and optimization. In the end, the resultant mechanism was roughly designed in SolidWorks in order to visualize the movements.
In future work, the mechanical design will be improved and the manufacturing process will take place. The plan is to carry out the validation procedure in terms of kinetics and kinematics. A kinematic synthesis and optimization procedure is handled with a software that is developed in a MATLAB environment. However, since there are many parameters that should be specified by the user, it is not easy to hard-code using the synthesis and optimization program once the different requirements have arisen. Since the system can be used in a different automation process or in a similar project, a user-friendly graphical user interface is planned to be utilized so that the user can alter the parameters easily. Besides the torque transmitting option, the kinematic movement optimization algorithms will also be applied to this interface so that the user can decide which parameter would be optimized according to their needs. Moreover, a visual simulation will be added to this interface so that the user will be able to see the results easily, apart from just seeing the numerical values. In the end, this kind of design process could be utilized by anybody who needs such a design.
Author Contributions: M.K.: conceptualization, methodology, software, validation, investigation, data curation, visualization, and writing-original draft preparation; M.V.Y.: conceptualization, methodology, software, validation, investigation, data curation, visualization, and writing-original draft preparation; E.A.: conceptualization, methodology, software, validation, investigation, data curation, visualization, and writing-original draft preparation; F.C.C.: conceptualization, methodology, formal analysis, writing-original draft preparation, writing-review and editing, supervision, and project administration; E.G.: conceptualization, methodology, formal analysis, writing-original draft preparation, writing-review and editing, supervision, and project administration. All authors have read and agreed to the published version of the manuscript.