Study on Design of Non-Circular Gears for Speed Control of the Squid Belly Opening and Gutting Machine (SBOGM)

: Non-circular gears can maintain rotational motions of general gears and implement all varying rotational motions of the cam. They adjust the angular velocity of driven gear according to operating conditions and make precise changes in angular motion. The design of non-circular gears has not been sufﬁciently studied because of their particularity and complex design methods unlike spur gears. In the gutting section of the Squid Belly Opening and Gutting Machine (SBOGM), spur gears generate rotational impact due to constant angular velocities, causing noise and equipment damage; so, efﬁciency should be improved by varying sectional angular velocity. Therefore, we derived pitch curves by selecting angular velocity ratio considering operating environments, and the tooth proﬁle was designed by calculating module for each section according to radius through theorical analysis for precise expression of angular velocity ratio. To conﬁrm reliability of design, angular velocity ratio and structural safety of designed non-circular gears were veriﬁed using, commercial software, ‘DAFUL 2020 R1 (cid:48) .


Introduction
The Squid Belly Opening and Gutting Machine (SBOGM) is a machine that automates processes such as belly opening, gutting, and separation of the body and legs of squids, which were previously performed manually by people, and all processes of the machine are carried out through gears and chain transmission as shown in Figure 1. Among the processes, viscera removal is done by fixing the viscera of squids with the inhaler and then plucking them out. Since the inhaler always moves at the same speed, there is the problem that the viscera of squids may not be removed completely depending on the type and size of the squid, and thus only squids of a specific type and size can be processed. In addition, in the section where the inhaler rotates and returns to its original position for repetitive processes, the heavy weight of the inhaler causes impact during rotation, generating noise, and at a high velocity, there is a possibility of damage to the machine; so, it is not possible to increase the processing speed due to concerns about damage to the machine.
As a result, the quantity of squids processed per hour of the SBOGM is limited, resulting in reduced work efficiency and economy.
In order to solve the problem, it is necessary to control the speed of the gutting section, which was previously controlled by the spur gear and chain drive, and although we initially attempted to use an electric motor for speed control, we intend to use a noncircular gear instead of an electric motor because the use environment is a salty area near the sea and it is difficult to use an electric motor that has relatively higher failure rates than a mechanical motor. Jian [1] obtained the angular velocity graph of driven gear through the design of high-order elliptic gears and dynamics simulation using software MAPLE. Seong-Cheol Lee [2] applied the average module for an arbitrary number of teeth to all teeth and connected the tooth roots in a circular arc form in the design of non-circular gears according to the angular velocity ratio diagram. Zhou Kaihong [3] showed the theoretical curve analysis of the angular velocity ratio and the consistency of analysis results for non-circular gears with an oval pitch curve by using the analysis software ADAMS. Marius VASIE [4] considered the simulation of rolling in a 2D space and verified the geometric and kinematic calculations performed in the design stage. Sang-Hoon Choi [5] studied optimal ranges of pressure angle, ratio of tooth number and radius of curvature of an elliptical gear. Dong-Woo Kang [6] represented the pitch curve of an elliptical gear in the form of a biarc and suggested a method of creating a tooth profile with an involute curve on each arc of biarc curves. Hao Xu [7] studied a general form of gears known as non-circular gears that can transfer periodic motion with variable speed through their irregular shapes and eccentric rotation centers. Sang-Hoon Choi [8] suggested one of references in applying non-circular gear to industrial plant, and suitable range of application by pressure angle curvature and angle ratio. Mircea Niculescu [9] introduced a new kinematics of the billet reheating furnace unloading door, by using a noncircular gear train in the drive mechanism. Marius Vasie [10] showed about the meshing conditions of non-circular gears. Fangyan Zheng [11] showed synthesizes still another special kind of indexing mechanism, using planetary non-circular gear train covering three gear trains with different gear types. Zhang Jian [12] proposed a simple algorithm for designing pitch curves for non-circular gears based on variable gear ratio functions. A A Lyashkov [13] proposed a geometrical model of solution to the task of profiling a non-circular gear, the centroid of which consists of interconnected arcs. D Mundo [14] presented a concept of epicyclical gear train able to generate a variable gear ratio law. CRISTESCU Ana [15] designed of multispeed gears pitch curves was taken into account, as the first important step in the gear design procedure. Specific original algorithms were dedicated to multispeed gear pitch curve modeling, based on the definition of gear transmission ratio variation law. Ming-Feng Tsay [16] derived a mathematical model of the tooth profile of a non-circular gear based on the approximate pitch curve. SABAH KHAN [17] measured the transmission error to determine the noise and vibration of non-circular using 'ANSYS'. A Cristescu [18] investigated multi-speed gear generation procedures and gear bending conditions. However, previous studies used highly similar pitch curve design methods and dealt with the cases where the pitch curve is elliptical or relatively gentle. Hence, there has been a lack of research on module selection and tooth profile design considering the operating environment as well as the curvature of the pitch curve of non-circular gears in relation to the characteristics of non-circular gears with angular velocity ratio varying continuously from moment to moment unlike general spur gears. to all teeth and connected the tooth roots in a circular arc form in the design of non-circular gears according to the angular velocity ratio diagram. Zhou Kaihong [3] showed the theoretical curve analysis of the angular velocity ratio and the consistency of analysis results for non-circular gears with an oval pitch curve by using the analysis software ADAMS. Marius VASIE [4] considered the simulation of rolling in a 2D space and verified the geometric and kinematic calculations performed in the design stage. Sang-Hoon Choi [5] studied optimal ranges of pressure angle, ratio of tooth number and radius of curvature of an elliptical gear. Dong-Woo Kang [6] represented the pitch curve of an elliptical gear in the form of a biarc and suggested a method of creating a tooth profile with an involute curve on each arc of biarc curves. Hao Xu [7] studied a general form of gears known as non-circular gears that can transfer periodic motion with variable speed through their irregular shapes and eccentric rotation centers. Sang-Hoon Choi [8] suggested one of references in applying non-circular gear to industrial plant, and suitable range of application by pressure angle curvature and angle ratio. Mircea Niculescu [9] introduced a new kinematics of the billet reheating furnace unloading door, by using a noncircular gear train in the drive mechanism. Marius Vasie [10] showed about the meshing conditions of noncircular gears. Fangyan Zheng [11] showed synthesizes still another special kind of indexing mechanism, using planetary non-circular gear train covering three gear trains with different gear types. Zhang Jian [12] proposed a simple algorithm for designing pitch curves for non-circular gears based on variable gear ratio functions. A A Lyashkov [13] proposed a geometrical model of solution to the task of profiling a non-circular gear, the centroid of which consists of interconnected arcs. D Mundo [14] presented a concept of epicyclical gear train able to generate a variable gear ratio law. CRISTESCU Ana [15] designed of multispeed gears pitch curves was taken into account, as the first important step in the gear design procedure. Specific original algorithms were dedicated to multispeed gear pitch curve modeling, based on the definition of gear transmission ratio variation law. Ming-Feng Tsay [16] derived a mathematical model of the tooth profile of a noncircular gear based on the approximate pitch curve. SABAH KHAN [17] measured the transmission error to determine the noise and vibration of non-circular using 'ANSYS'. A Cristescu [18] investigated multi-speed gear generation procedures and gear bending conditions. However, previous studies used highly similar pitch curve design methods and dealt with the cases where the pitch curve is elliptical or relatively gentle. Hence, there has been a lack of research on module selection and tooth profile design considering the operating environment as well as the curvature of the pitch curve of non-circular gears in relation to the characteristics of non-circular gears with angular velocity ratio varying continuously from moment to moment unlike general spur gears. Therefore, in this study, we intend to design a non-circular gear that can be applied to a special use environment where the velocity should be decreased to reduce impact in the rotation section while the transportation speed of the inhaler should be rapidly increased to increase the force for viscera removal. The ratio of angular velocity was selected Therefore, in this study, we intend to design a non-circular gear that can be applied to a special use environment where the velocity should be decreased to reduce impact in the rotation section while the transportation speed of the inhaler should be rapidly increased to increase the force for viscera removal. The ratio of angular velocity was selected in consideration of the operating environment and expressed as a function for each section. The curvature of each section of the derived pitch curve was calculated, and modules were calculated for each section, based on the transfer load and radius of the pitch. The driving gear was designed by applying the involute tooth shape determined for each section, and the driven gear for the designed driving gear was generated using the driven gear generation program. Finally, the ratio of angular velocity and structural safety of the designed non-circular gear were verified by conducting a flexible body dynamics analysis using 'DAFUL' according to design procedures shown in Figure 2.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 3 of 14 in consideration of the operating environment and expressed as a function for each section. The curvature of each section of the derived pitch curve was calculated, and modules were calculated for each section, based on the transfer load and radius of the pitch. The driving gear was designed by applying the involute tooth shape determined for each section, and the driven gear for the designed driving gear was generated using the driven gear generation program. Finally, the ratio of angular velocity and structural safety of the designed non-circular gear were verified by conducting a flexible body dynamics analysis using 'DAFUL' according to design procedures shown in Figure 2.   In order to design the pitch curve, the angular velocity ratio of the viscera removal section needs to be selected. In the existing equipment, the section was operated through a spur gear at a constant angular velocity of 60 rpm, but to reduce rotational impact and improve work efficiency, the velocity was increased by 30% compared to that of the existing equipment in the section where the actual gutting process occurs ( 2 ~3 ), while the velocity was decreased by 30% compared to that of the existing equipment to reduce rotational impact in the section in which the inhaler rotates ( 4 ~5 , 7 ~8 ). The entire acceleration and deceleration sections are as shown in Figure 3.

Design of the Angular Velocity Ratio
In order to design the pitch curve, the angular velocity ratio of the viscera removal section needs to be selected. In the existing equipment, the section was operated through a spur gear at a constant angular velocity of 60 rpm, but to reduce rotational impact and improve work efficiency, the velocity was increased by 30% compared to that of the existing equipment in the section where the actual gutting process occurs (②~③), while the velocity was decreased by 30% compared to that of the existing equipment to reduce rotational impact in the section in which the inhaler rotates (④~⑤, ⑦~⑧). The entire acceleration and deceleration sections are as shown in Figure 3. The graph of the angular velocity ratio of the driven gear to the driving gear according to the selected acceleration and deceleration sections is shown in Figure 4. In the case of non-circular gears, the ratio of the driving gear to the driven gear must be designed with a natural number, and the rotation starting point of the driving gear and driven gear must always be constant. To satisfy these conditions, the areas of deceleration and acceleration sections in the angular velocity ratio graph of the driving gear must be identical.
Here, to express each section as a single function, a combination of functions for each section was selected as shown in Table 1 so that the average angular velocity ratio per period for each section would be 1. The average angular velocity ratio was set as 1, and 'h' was set as 0.3 to increase or decrease the angular velocity by 30% in the acceleration/deceleration section.  The graph of the angular velocity ratio of the driven gear to the driving gear according to the selected acceleration and deceleration sections is shown in Figure 4. In the case of non-circular gears, the ratio of the driving gear to the driven gear must be designed with a natural number, and the rotation starting point of the driving gear and driven gear must always be constant. To satisfy these conditions, the areas of deceleration and acceleration sections in the angular velocity ratio graph of the driving gear must be identical.

Design of the Angular Velocity Ratio
In order to design the pitch curve, the angular velocity ratio of the viscera removal section needs to be selected. In the existing equipment, the section was operated through a spur gear at a constant angular velocity of 60 rpm, but to reduce rotational impact and improve work efficiency, the velocity was increased by 30% compared to that of the existing equipment in the section where the actual gutting process occurs (②~③), while the velocity was decreased by 30% compared to that of the existing equipment to reduce rotational impact in the section in which the inhaler rotates (④~⑤, ⑦~⑧). The entire acceleration and deceleration sections are as shown in Figure 3. The graph of the angular velocity ratio of the driven gear to the driving gear according to the selected acceleration and deceleration sections is shown in Figure 4. In the case of non-circular gears, the ratio of the driving gear to the driven gear must be designed with a natural number, and the rotation starting point of the driving gear and driven gear must always be constant. To satisfy these conditions, the areas of deceleration and acceleration sections in the angular velocity ratio graph of the driving gear must be identical.
Here, to express each section as a single function, a combination of functions for each section was selected as shown in Table 1 so that the average angular velocity ratio per period for each section would be 1. The average angular velocity ratio was set as 1, and 'h' was set as 0.3 to increase or decrease the angular velocity by 30% in the acceleration/deceleration section.  Here, to express each section as a single function, a combination of functions for each section was selected as shown in Table 1 so that the average angular velocity ratio per period for each section would be 1. The average angular velocity ratio was set as 1, and 'h' was set as 0.3 to increase or decrease the angular velocity by 30% in the acceleration/deceleration section.

Design of the Pitch Curve
In non-circular gears, the rotation angle (ϕ) of the driven gear can be obtained differently from the rotation angle (θ) of the driving gear with a constant angular velocity, and the rotation angle (ϕ) of the driven gear is expressed by Equation (2).
Here, d∅ dθ = f (θ) represents the angular velocity ratio for the rotation angle (∅) of the driven gear to the rotation angle (θ) of the driving gear.
As shown in Figure 5, if the distance (c) between the rotation axes of the driving and driven gears is constant as shown in Equation (3) and the gear ratio is 1:1, since the pitch curves are in rolling contact with each other, rolling distance l 1 of the pitch curve of the driving gear is equal to the rolling distance l 2 of the pitch curve of the driven gear. Therefore, an infinitesimal change in each rolling distance at an arbitrary pitch point can be expressed by Equation (4). where, Appl. Sci. 2021, 11, x FOR PEER REVIEW 5 of 14

Design of the Pitch Curve
In non-circular gears, the rotation angle (φ) of the driven gear can be obtained differently from the rotation angle (θ) of the driving gear with a constant angular velocity, and the rotation angle (φ) of the driven gear is expressed by Equation (2).
Here, ∅ = ( ) represents the angular velocity ratio for the rotation angle (∅) of the driven gear to the rotation angle ( ) of the driving gear. As shown in Figure 5, if the distance (c) between the rotation axes of the driving and driven gears is constant as shown in Equation (3) and the gear ratio is 1:1, since the pitch curves are in rolling contact with each other, rolling distance of the pitch curve of the driving gear is equal to the rolling distance of the pitch curve of the driven gear. Therefore, an infinitesimal change in each rolling distance at an arbitrary pitch point can be expressed by Equation (4).
where, = ∅ (4) Figure 5. Pitch curves of driving gear and driven gear. Figure 5. Pitch curves of driving gear and driven gear.
By combining Equations (3) and (4), r 1 and r 2 are derived as follows: The length of the pitch curve of each section is given by Equation (6).
The total length of the pitch curve is calculated by Equation (7).
In other words, if the angular velocity ratio f (θ) and the inter-axis distance 'c' are determined, the pitch curve can be determined from the derived r 1 (θ) and r 2 (θ), and its results were shown in Table 2. If f (θ) obtained in Table 1 is substituted into Equation (5), the coordinates of the pitch curves of the driving and driven gears can be obtained. In this study, the inter-axis distance was set as 100 mm according to the actual operating environment, and for data reduction for the ease of 2D and 3D modeling, the error of the pitch curve according to the number of coordinate points was examined. When using 3600 coordinates obtained in every case of θ = 0.1 • , the maximum distance error was generated as 1 × 10 −6 mm or less.

Tooth Profile Design
Unlike general spur gears that maintain the same velocity, non-circular gears have continuously changing angular velocity of the driven gear during one rotation of the driving gear, so it is difficult to design a non-circular gear with the same module. So, minimum module was determined through theoretical equation for each section, satisfied with structural safety and obtained precision of angular velocity described.

Calculation of Minimum Module
In order to determine module for each section, it is necessary to calculate the minimum module that can be used in the operating environment. In the case of non-circular gears, the distance from the center of rotation to the pitch point continuously changes during the rotation of gears which are engaged with each other. The largest changes in angular velocity occur in the section with the largest curvature, and in order to accurately represent the angular velocity, module was based on the position with the smallest radius in the driving gear, and a sufficient safety factor was given in consideration of the specificity of non-circular gears. For the transfer load (W), the minimum module is given by Equation (8).
The minimum module is obtained from Equation (8), where w is the transfer load, s is the safety factor, f v is the velocity factor, y is the tooth form factor, and σ b is the bending stress in gear tooth.
The torque load required in the field is 50 000 N·mm in the viscera removal section, and the transfer load W n is calculated by Equation (9). Each interval of driving gear is shown in Figure 6, and theoretically calculated module is shown in Table 3.
. . Natural number) (10) gears, the distance from the center of rotation to the pitch point continuously change during the rotation of gears which are engaged with each other. The largest changes i angular velocity occur in the section with the largest curvature, and in order to accuratel represent the angular velocity, module was based on the position with the smallest radiu in the driving gear, and a sufficient safety factor was given in consideration of the spec ficity of non-circular gears. For the transfer load (W), the minimum module is given b Equation (8).
The minimum module is obtained from Equation (8), where w is the transfer load, is the safety factor, is the velocity factor, y is the tooth form factor, and is the bend ing stress in gear tooth.
The torque load required in the field is 50′000 N • mm in the viscera removal section and the transfer load is calculated by Equation (9). Each interval of driving gear shown in Figure 6, and theoretically calculated module is shown in Table 3.  The minimum module obtained from Equation (8)    The minimum module obtained from Equation (8) was adjusted by Equation (10) because the number of teeth within each section must be a natural number.
Design procedures of non-circular gear are composed of 4 steps and are as follows.

1.
Creation of involute tooth according to the modules calculated (Figure 7a).

2.
Equal division of pitch curve at each section according to the modules adjusted (Figure 7b).

3.
Arrangement of centers (pitch points) of circle in the equal points by translation and rotation (Figure 7c).
3. Arrangement of centers (pitch points) of circle in the equal points by translation and rotation (Figure7c). 4. Connection of roots between designed teeth (Figure 7d). Proceeding in the above order, the final driving gear design is completed. The shape of the final designed driving gear is shown in Figure 8.

Design of the Tooth Profile of the Driven Gear
The design of the driven gear was carried out using the gear design program 'Gearify' to ensure that the driven gear is precisely engaged with the driving gear while it is rotating. The pitch curve of the driving gear and its final shape were imported into 'Gearify' as shown in Figure 9.
'Extra High Accuracy' was selected to derive a precise design of the driven gear, and its teeth were smoothly connected by setting 'Integration Accuracy' as high as shown in Figure 10.  Proceeding in the above order, the final driving gear design is completed. The shape of the final designed driving gear is shown in Figure 8.
of the final designed driving gear is shown in Figure 8.

Design of the Tooth Profile of the Driven Gear
The design of the driven gear was carried out using the gear design program 'Gearify' to ensure that the driven gear is precisely engaged with the driving gear while it is rotating. The pitch curve of the driving gear and its final shape were imported into 'Gearify' as shown in Figure 9.
'Extra High Accuracy' was selected to derive a precise design of the driven gear, and its teeth were smoothly connected by setting 'Integration Accuracy' as high as shown in Figure 10.

Design of the Tooth Profile of the Driven Gear
The design of the driven gear was carried out using the gear design program 'Gearify' to ensure that the driven gear is precisely engaged with the driving gear while it is rotating. The pitch curve of the driving gear and its final shape were imported into 'Gearify' as shown in Figure 9. Sci. 2021, 11, x FOR PEER REVIEW 9 of 14  If the designed pitch curve and gear curve of the driving gear are input into the program which is proceeded with, the driven gear is generated through 'Gearify' and is shown in Figure 11.

3D Modeling
Flexible body dynamic analysis was conducted to examine the angular velocities of designed driving and driven gears and their structural safety was confirmed using the dynamic analysis program 'DAFUL'; analytical modeling was carried out using 'Space-Claim'.
Tooth thickness (10 mm) was the same as the chain gear and the spur gear used in the actual field, and gears were divided into two parts (tooth part and body part) in order to reduce analysis time, as shown in Figure 12, offsetting pitch curve. 'Extra High Accuracy' was selected to derive a precise design of the driven gear, and its teeth were smoothly connected by setting 'Integration Accuracy' as high as shown in Figure 10.  If the designed pitch curve and gear curve of the driving gear are input into the program which is proceeded with, the driven gear is generated through 'Gearify' and is shown in Figure 11.

3D Modeling
Flexible body dynamic analysis was conducted to examine the angular velocities of designed driving and driven gears and their structural safety was confirmed using the dynamic analysis program 'DAFUL'; analytical modeling was carried out using 'Space-Claim'.
Tooth thickness (10 mm) was the same as the chain gear and the spur gear used in the actual field, and gears were divided into two parts (tooth part and body part) in order If the designed pitch curve and gear curve of the driving gear are input into the program which is proceeded with, the driven gear is generated through 'Gearify' and is shown in Figure 11.  If the designed pitch curve and gear curve of the driving gear are input into the program which is proceeded with, the driven gear is generated through 'Gearify' and is shown in Figure 11.

3D Modeling
Flexible body dynamic analysis was conducted to examine the angular velocities of designed driving and driven gears and their structural safety was confirmed using the dynamic analysis program 'DAFUL'; analytical modeling was carried out using 'Space-Claim'.
Tooth thickness (10 mm) was the same as the chain gear and the spur gear used in the actual field, and gears were divided into two parts (tooth part and body part) in order

3D Modeling
Flexible body dynamic analysis was conducted to examine the angular velocities of designed driving and driven gears and their structural safety was confirmed using the dynamic analysis program 'DAFUL'; analytical modeling was carried out using 'SpaceClaim'.
Tooth thickness (10 mm) was the same as the chain gear and the spur gear used in the actual field, and gears were divided into two parts (tooth part and body part) in order to reduce analysis time, as shown in Figure 12, offsetting pitch curve.

Boundary Conditions
Based on analytical modeling, the mesh for the tooth part was created with Workbench'. The mesh was controlled by edge sizing and the sweep method, an consisting of 231′612 elements and 1′226′260 nodes were generated through the of mesh quality as shown in Figure 13. The generated mesh data was exported a file' to transfer it to 'Daful'. After importing the generated 'Mesh file' from 'Daful' to conduct dynamic the 'Joint condition' was given so that the driving and driven gears could rotat the rotation axis. Also, considering the actual operating environment, the rotation of 60 rpm was applied for the driving gear and the torque load of 50′000 Nmm wa for the driven gear. For the convergence of analysis, as shown in Figure 14, r velocity and torque load were set so that they would gradually increase with tim ing the 'Step function', as shown in Figure 15.

Boundary Conditions
Based on analytical modeling, the mesh for the tooth part was created with 'Ansys-Workbench'. The mesh was controlled by edge sizing and the sweep method, and a mesh consisting of 231 612 elements and 1 226 260 nodes were generated through the analysis of mesh quality as shown in Figure 13. The generated mesh data was exported as an '.inp file' to transfer it to 'Daful'.

Boundary Conditions
Based on analytical modeling, the mesh for the tooth part was created with ' Workbench'. The mesh was controlled by edge sizing and the sweep method, and consisting of 231′612 elements and 1′226′260 nodes were generated through the a of mesh quality as shown in Figure 13. The generated mesh data was exported as a file' to transfer it to 'Daful'. After importing the generated 'Mesh file' from 'Daful' to conduct dynamic an the 'Joint condition' was given so that the driving and driven gears could rotate the rotation axis. Also, considering the actual operating environment, the rotationa of 60 rpm was applied for the driving gear and the torque load of 50′000 Nmm was a for the driven gear. For the convergence of analysis, as shown in Figure 14, rot velocity and torque load were set so that they would gradually increase with time ing the 'Step function', as shown in Figure 15. After importing the generated 'Mesh file' from 'Daful' to conduct dynamic analysis, the 'Joint condition' was given so that the driving and driven gears could rotate around the rotation axis. Also, considering the actual operating environment, the rotational speed of 60 rpm was applied for the driving gear and the torque load of 50 000 Nmm was applied for the driven gear. For the convergence of analysis, as shown in Figure 14, rotational velocity and torque load were set so that they would gradually increase with time by using the 'Step function', as shown in Figure 15.
the 'Joint condition' was given so that the driving and driven gears could rotate the rotation axis. Also, considering the actual operating environment, the rotation of 60 rpm was applied for the driving gear and the torque load of 50′000 Nmm was for the driven gear. For the convergence of analysis, as shown in Figure 14, ro velocity and torque load were set so that they would gradually increase with tim ing the 'Step function', as shown in Figure 15. In the case of contact conditions, as shown in Figure 16, the 'Friction Coeffici set as 0.08 according to the general friction conditions of gears. The total analysis t 3 s, and the 'Output Step' was set as 2000 to obtain detailed analysis results.

Analysis Results
The angular velocity ratio during one revolution of the driving gear was extracted through 'DAFUL' analysis, based on the time points when the velocity was sufficiently developed.
The results are shown in Figure 17, and it can be seen that analysis results of angular velocity ratio are identical to the graph of angular velocity ratio obtained by theoretical calculations. In the case of contact conditions, as shown in Figure 16, the 'Friction Coefficient' was set as 0.08 according to the general friction conditions of gears. The total analysis time was 3 s, and the 'Output Step' was set as 2000 to obtain detailed analysis results.

Analysis Results
The angular velocity ratio during one revolution of the driving gear was extracted through 'DAFUL' analysis, based on the time points when the velocity was sufficiently developed.
The results are shown in Figure 17, and it can be seen that analysis results of angular velocity ratio are identical to the graph of angular velocity ratio obtained by theoretical calculations.

Analysis Results
The angular velocity ratio during one revolution of the driving gear was extracted through 'DAFUL' analysis, based on the time points when the velocity was sufficiently developed.
The results are shown in Figure 17, and it can be seen that analysis results of angular velocity ratio are identical to the graph of angular velocity ratio obtained by theoretical calculations.

Analysis Results
The angular velocity ratio during one revolution of the driving gear was extracted through 'DAFUL' analysis, based on the time points when the velocity was sufficiently developed.
The results are shown in Figure 17, and it can be seen that analysis results of angular velocity ratio are identical to the graph of angular velocity ratio obtained by theoretical calculations. In order to verify the structural safety, Von-Mises stresses of driving and driven gears were extracted. Figure 18 show the stress distribution of each section during one revolution of the driving gear and driven gear, respectively. In the sections of corresponding to circular motion with a constant speed, Von-Mises stress is relatively constant, but in the sections where the angular velocity of the driven gear increases or decreases, stress is higher than in the constant velocity section. In order to verify the structural safety, Von-Mises stresses of driving and driven gears were extracted. Figure 18 show the stress distribution of each section during one revolution of the driving gear and driven gear, respectively. In the sections of corresponding to circular motion with a constant speed, Von-Mises stress is relatively constant, but in the sections where the angular velocity of the driven gear increases or decreases, stress is higher than in the constant velocity section.  Table 4 showed the maximum equivalent stress at each section of the driving and driven gears considering safety factor 1.2, and the maximum stress in Section (8) was 411.45 MPa, which was lower than the yield strength (SUS440C) as shown in Figure 19, i.e., 445 MPa so that structural evaluation of non-circular gear was considered to be safe.  Table 4 showed the maximum equivalent stress at each section of the driving and driven gears considering safety factor 1.2, and the maximum stress in Section (8) was 411.45 MPa, which was lower than the yield strength (SUS440C) as shown in Figure 19, i.e., 445 MPa so that structural evaluation of non-circular gear was considered to be safe.

Conclusions
In this study, we intend to design non-circular gears that can be applied to a special environment where the velocity should be decreased to reduce impact in the rotation section while the transportation speed of the inhaler should be rapidly increased to increase

Conclusions
In this study, we intend to design non-circular gears that can be applied to a special environment where the velocity should be decreased to reduce impact in the rotation section while the transportation speed of the inhaler should be rapidly increased to increase the force for viscera removal. And the ratio of angular velocity and structural safety of the designed non-circular gear were verified by conducting a flexible body dynamics analysis using 'DAFUL'. The results of the study are summarized as follows.
(1) To improve the efficiency of the processes of the existing Squid Belly Opening and Gutting Machine (SBOGM) and to reduce impact during rotation, non-circular gears alternative to the chain gear and spur gear were developed for the SBOGM. (2) Driving gear was designed according to the tooth design procedures suggested for non-circular gears, and the driven gear was generated through its program 'Gearify' based on the designed driving gear. (3) The angular velocity ratio during one revolution of the driving gear was extracted through 'DAFUL' analysis, whose results are identical to the graph of angular velocity ratio obtained by theoretical calculations. (4) Flexible body dynamics analysis was conducted using the commercial program 'DAFUL', and the maximum Von-Mises stress at the Section (8) was 411.45 MPa, which was lower than the yield strength of 445 MPa, so the structural safety was considered to be safe.