Numerical Measurement of a Virtual Model for the Knee Prosthesis Geometry

: The human knee joint usually suffers progressive deterioration with time. The conventional cure of this issue is to replace it with an alternate knee by applying the prosthesis implant. The reason is that the process causes the abrasion of the different materials rather than just sliding or rolling. This study aims to develop the numerical measurement of the knee prosthesis’s geometry, which fulﬁls the mechanical requirements of the human knee. The MSC.ADAMS programme was applied to demonstrate the movement of the human knee joint in terms of rotation and ﬂexion. The changes between the condyles of the developed multibody of the prosthesis related to the ﬂexion angle ranging from 20–120 ◦ were investigated and presented. The boundary conditions were determined, and simulations performed using the ADAM’s programme. An average value of 0.7 was reached for the slip ration, with the maximum getting up to 0.79. An angle between 110–120 ◦ for the ﬂexion angle was obtained. It can be said that the application of the multibody model saves time as there is no involvement of the tibia and the femur as required for the knee prosthesis. More importantly, as the application of the test machine is omitted in our process, our model’s approximations to a human knee are carried out directly. Without cost, several measurements for the knee prosthesis could be made and repaired. The study results provide the necessary insight for future tests regarding the movement of the knee joint. the two prosthetic components together according to the objective function of ﬂexion and rotation of the knee prosthesis. The results showed a close match with the curve of the movement of the normal human knee.


Introduction
Currently, knee replacement surgery is performed on a large number of people throughout the world. Modern technologies would be required to contribute to this subject within the near future. This means that techniques used for virtual prototyping of mechanical systems must be extended to advanced biomechanical systems to meet the specifications and demands [1][2][3][4][5][6][7].
The most efficient process for correcting deformities in the leg and removing pains associated with the knee, and enabling sanctionative patients to resume traditional daily activities, is by Total knee replacement (TKR) surgery. Since the primary TKR surgery was performed in 1968, enhancements with different materials and methods have drastically increased its effectiveness [7][8][9][10][11].
The proper handling of degenerative abrasion in the knee joint is currently one of the essential orthopaedic problems. The appropriate solution for this issue is to remove and replace it with a prosthesis. Unfortunately, no perfect knee prosthesis is available to replace the original knee because of the knee kinematics' complexity.
Results for the kinematical design method of the knee prosthesis geometries by Balassa, 2019, compared to the normal human knee [12], shows that the basic principle of their method was to move the two prosthetic components together according to the objective function of flexion and rotation of the knee prosthesis. The results showed a close match with the curve of the movement of the normal human knee.

Materials and Methods
Generally, for the typical normal human knee flexion or extension, the patellofemoral joint's local kinematics can be classified as partial rolling and sliding. The associated movement is controlled by connecting the tibial surfaces, the femoral, and the ligaments.
The exact relationship of the sliding rolling phenomenon regarding the knee's typical functional arc has still not been ascertained. However, it has been used by the early works of [18,19]. Between the angles of 20-30 • of flexion, rolling appears to be more significant, while beyond these angles, the sliding angle becomes more prevailing [20,21].
The "Screw home mechanism", which occurs in the arc between 0 • and 20 • , is of tremendous interest for an anatomist, although it may have little benefit in the daily living activities and activities such as one-legged stance or normal stance [22,23].
In the current work, the maximum arc angle of 110 • was considered. This is because there is no increase in the patellofemoral forces above 120 • of flexion. Also, in areas like Asia, where it appears in part of their cultures, its relevance is insignificant [24,25].

Knee Joint Geometry
The tibia, patella, and femur were taken to be rigid bodies as the effect of deformation is not considered in this study. The femur and tibia geometry was based on the prosthesis prototypes developed by Balassa 2019, as shown in Figures 1 and 2. These designs are still being researched and refined. The designed samples of the prosthesis are shown in Figures 1 and 2. As the basis for our model in the ADAMS software, these initial prosthesis prototypes were transferred together with the 3D model into the ADAMS programme as Appl. Sci. 2021, 11, 2541 3 of 18 the initial framework for the prosthesis design. The obtained STL file is converted to an IGS file and finally into the PARASOLID format for application in the ADAMS program.

Knee Joint Geometry
The tibia, patella, and femur were taken to be rigid bodies as the effect of deformation is not considered in this study. The femur and tibia geometry was based on the prosthesis prototypes developed by Balassa 2019, as shown in Figures 1 and 2. These designs are still being researched and refined. The designed samples of the prosthesis are shown in Figures 1 and 2. As the basis for our model in the ADAMS software, these initial prosthesis prototypes were transferred together with the 3D model into the ADAMS programme as the initial framework for the prosthesis design. The obtained STL file is converted to an IGS file and finally into the PARASOLID format for application in the ADAMS program.

Ligaments and Muscles
The various ligaments and the muscle are designed as simple linear springs. Considering both the ligaments and the muscles only in arithmetic, the stiffness modulus is set to 130 N/mm and the damping modulus to 0.15 Ns/mm in the case of both springs.

Mu Static
The coefficient of friction is the standard describing the force of friction between 2 objects and the natural reaction found for the objects in question. This value is used to determine the friction force of an object [26,27] and when using other methods as shown in Equation (1), the coefficient of friction. (1) where : The frictional force : The coefficient of friction : The normal force.

Knee Joint Geometry
The tibia, patella, and femur were taken to be rigid bodies as the effect of deformation is not considered in this study. The femur and tibia geometry was based on the prosthesis prototypes developed by Balassa 2019, as shown in Figures 1 and 2. These designs are still being researched and refined. The designed samples of the prosthesis are shown in Figures 1 and 2. As the basis for our model in the ADAMS software, these initial prosthesis prototypes were transferred together with the 3D model into the ADAMS programme as the initial framework for the prosthesis design. The obtained STL file is converted to an IGS file and finally into the PARASOLID format for application in the ADAMS program.

Ligaments and Muscles
The various ligaments and the muscle are designed as simple linear springs. Considering both the ligaments and the muscles only in arithmetic, the stiffness modulus is set to 130 N/mm and the damping modulus to 0.15 Ns/mm in the case of both springs.

Mu Static
The coefficient of friction is the standard describing the force of friction between 2 objects and the natural reaction found for the objects in question. This value is used to determine the friction force of an object [26,27] and when using other methods as shown in Equation (1), the coefficient of friction. (1) where : The frictional force : The coefficient of friction : The normal force.

Ligaments and Muscles
The various ligaments and the muscle are designed as simple linear springs. Considering both the ligaments and the muscles only in arithmetic, the stiffness modulus is set to 130 N/mm and the damping modulus to 0.15 Ns/mm in the case of both springs.

Mu Static
The coefficient of friction is the standard describing the force of friction between 2 objects and the natural reaction found for the objects in question. This value is used to determine the friction force of an object [26,27] and when using other methods as shown in Equation (1), the coefficient of friction.
where F f : The frictional force µ : The coefficient of friction F n : The normal force. µ stands for different items. It can be the coefficient of Static friction µ s , and also the coefficient of dynamic friction.
The expression for the frictional force is shown in Equation (2): F f = The frictional force (in Newtons). µ = The static (µ s ) or kinetic (µ k ) frictional coefficient (dimensional) and N stands for the normal force (in Newtons). µ depends on the rules that cause friction. Its magnitude ranges from zero to one but sometimes bigger than one. When it is zero, it indicates the presence of friction between the elements; such is possible with superfluidity. All things other than that will have some friction when they touch each other; the force of friction is equal to the normal force at the value 1, and the coefficient of friction is not just restrained within the range of 0 and 1. The coefficient of friction whose value exceeds the virtual force, for example, silicone rubber, the value of Mu static frictional force chosen for the modelling is 0.4, as shown in Table 1. Table 1. Definition of the various functions to determine the points connecting the femoral or tibial surfaces.

D c (t)
The scalar vector function determines the position at a particular point for the contact between the bodies involved presented in the coordinate system ( Figure 5). where i = 1 represents the femur and tibia contact, and i = 2 the contact between the patella and the femur.
The velocity scalar function determines the tangent vector at a particular position of the contact point (P) of the femoral or tibial surfaces, respectively, as shown in Figure 5.
These functions determine the current position of the centre of mass and also determine the angular velocity of the femur (F) and tibia (T), as shown in ( Figure 5).

Mu Dynamic
Strong friction opposes the lateral sliding movement between two objects in contact; they are divided into two main groups, static friction occurs between stationary objects, and kinetic friction occurs between objects or surfaces in motion as it is called (dynamic friction) (Air, 1921).
Coulomb friction is applied to dry friction, and the expression for its determination is presented in Equation (3): where F f is the frictional force resulting from the collision of the two surfaces and the application of the forces to each other. It is in line with the surface and applied in the opposite of the force.
µ is the coefficient of friction; this reveals the relationship and the normal reaction that exists between two bodies. F n is the average force obtained on each surface, directed perpendicular to the surface. The value of the Mu dynamic frictional force chosen for the modelling is 0.3, as shown in Table 1.

Description of the Calculations Employed in the ADAMS Model
In this section, the kinematic calculations for the motion as a function of time and automatically solved by the ADAMS programme are presented.
Because the multibody is considered solid, the solid body kinematics is applicable. The ration of the sliding roll is determined only between the femur and the tibia. This is what was noticed, that the patella does not appear in the account or the figures. The following Equations (4)-(6) were used to determine the velocity at a point of the connecting femoral or tibial surfaces ( Figure 3).
From Equation (5), we obtain (7) and (8): Velocities relating to the femur and tibia at the point of contact have been determined in the absolute coordinate system ( Figure 4).
By substituting Equation (5) into (4), we obtain: From Equation (5), we obtain (7) and (8): Velocities relating to the femur and tibia at the point of contact have been determined in the absolute coordinate system ( Figure 4). Appl. Sci. 2021, 11, x FOR PEER REVIEW 6 of 18 At the unit vector , we multiply Equations (7) and (8): The tangential scalar components are justified only under the following condition [28]. (11) This indicates that the typical normal standard components of the femoral and tibia contact velocities should be the same. Otherwise, there would be a collision or separation of the two surfaces.
By combining scalar velocities over time, the length of the continuous arc with respect to the femur and tibia can be calculated as follows: The sliding scroll ratio can be entered back by specifying the arc lengths on both connected bodies: where, These rates show the differences increasing with the longest connected arc. The sliding rolling ration is the difference between the distance travelled (∆ ) on the tibia and the additional distance travelled (∆ ) on the femur bone over the At the unit vector P c1 (t), we multiply Equations (7) and (8): The tangential scalar components are justified only under the following condition [28].
This indicates that the typical normal standard components of the femoral and tibia contact velocities should be the same. Otherwise, there would be a collision or separation of the two surfaces.
By combining scalar velocities over time, the length of the continuous arc with respect to the femur and tibia can be calculated as follows: (13) The sliding scroll ratio can be entered back by specifying the arc lengths on both connected bodies: where, These rates show the differences increasing with the longest connected arc. The sliding rolling ration is the difference between the distance travelled (∆S tibiaN ) on the tibia and the additional distance travelled (∆S f emurN ) on the femur bone over the increased distance travelled (∆S tibiaN ) on the tibia. N shows the arc length during the connection.
Thus, we can conclude accurate calculations about the sliding and gradient features of motion. The sliding toll ratio that translates to zero means this is pure rolling. Regarding a ratio between 0-1, the movement turns into partial rolling and sliding.
Determining the sliding roll ratio as a function of the bending angle is better than determining it as a function of time. We will explain what we have done with this sliding rolling ratio. The flexion angle (γ) was derived, adding the femur and tibia's angular velocities around the X-axis, and was set to an initial degree of 20 • of squatting.
By defining a function, the time can be altered by the angle of flexion, and rolling can be plotted in correlation with the angle of flexion.

The Multibody Model
The multibody model was created by applying the following procedures • The general point motion was used to stabilize the distal femur, where all the coordinates are shown ( Figure 5). This enables the distal femur to make a transitional movement along the y-axis.

•
The cylindrical joint model was used to restrict the knee part to allow rotation around all axes ( Figure 5). This enables the shin bone to conduct a natural rotation.

•
We only considered the patellar tendon and the rectus femur in the numericalkinematical model. We create both of them as simple linear springs. According to the literature, the rectal femoral stiffness modulus was determined between 25 and 100 N/mm, according to the literature [29,30]. As an average value, we set it to 80 N/mm. With the stabilization factor set at 0.15 Ns/mm, for all the strings to prevent oscillations in the system, the patellar tendon was set to inextensible ( Figure 5). • According to Coulomb's law, contact restrictions are established concerning static and low dynamic friction coefficient (µs = 0.4 µd = 0.3) between the femur, tibia, and patella, similarly to real joints. The kinetic relationship between systemic forces, frictional forces (Fn, Fs), and flexion angle is analyzed using this constraint.

Boundary Conditions for the Simulation
After the geometrical model is obtained, the MSC.ADAMS program was used to build the multibody model. The following boundary conditions were applied to our model (prosthesis geometry): Figure 6 presents the flow chart of the applied steps of the multibody virtual model created in the ADAMS software After the validation and acceptance of the results, the application stage will be to print the femur and the tibia by using 3D printing as in Figure 7. Appl. Sci. 2021, 11, x FOR PEER REVIEW 8 of 18

Boundary Conditions for the Simulation
After the geometrical model is obtained, the MSC.ADAMS program was used to build the multibody model. The following boundary conditions were applied to our model (prosthesis geometry): Figure 6 presents the flow chart of the applied steps of the multibody virtual model created in the ADAMS software After the validation and acceptance of the results, the application stage will be to print the femur and the tibia by using 3D printing as in Figure 7.

Results and Discussion
This section presents the results and the accompanying discussions for the study. The movement of the knee in the sagittal plane (lateral movement), flexion position, and the last stage of extension, and in the medial rotation or transverse motion are presented in

Results and Discussion
This section presents the results and the accompanying discussions for the study. The movement of the knee in the sagittal plane (lateral movement), flexion position, and the last stage of extension, and in the medial rotation or transverse motion are presented in Figure 7. Femur, Tibia, and the guide.

Results and Discussion
This section presents the results and the accompanying discussions for the study. The movement of the knee in the sagittal plane (lateral movement), flexion position, and the last stage of extension, and in the medial rotation or transverse motion are presented in the flexion and rotation diagrams. Thus, showing the measured ranges of the flexion and extension and rotation of the knee prosthesis during the simulation process. These figures are, in essence, used to determine the point of sliding and rolling.

The Virtual Multibody Model
The ADAMS programme could compute the forces directly. At first, we saved it as PARASOLID, and we imported it into the MSC.ADAMS. The flexion angle was derived by combining the femur and tibia's angular velocities about the x-axis. This was done considering that the model was at 20 • for the sliding and rolling at the start of the movement. The angles were divided into three to tackle the three-dimensional movement. The results are summarized in Figure 8.
To be able to describe all the coordinates, we have restricted the distal femur by the general point motion, as shown in Figure 8. The knee model was restricted by a cylindrical joint, which allows the flexion process between a femur and tibia. Simple linear springs are designed as the boundary between the rectus femur and the patellar tendon [29,30].
According to Coulomb's law, the limitations of contact between the femur and the patella tibia are established for low static and dynamic friction coefficients (µs = 0.4 µd = 0.3), similar to human joints [31]. A force vector was created on the femur distal, as shown in Figure 9, and the value set at 800 N, while defining it by a step function (A, x0, h0, x1, h1).

The Virtual Multibody Model
The ADAMS programme could compute the forces directly. At first, we saved it as PARASOLID, and we imported it into the MSC.ADAMS. The flexion angle was derived by combining the femur and tibia's angular velocities about the x-axis. This was done considering that the model was at 20° for the sliding and rolling at the start of the movement. The angles were divided into three to tackle the three-dimensional movement. The results are summarized in Figure 8.
To be able to describe all the coordinates, we have restricted the distal femur by the general point motion, as shown in Figure 8. The knee model was restricted by a cylindrical joint, which allows the flexion process between a femur and tibia. Simple linear springs are designed as the boundary between the rectus femur and the patellar tendon [29,30].
According to Coulomb's law, the limitations of contact between the femur and the patella tibia are established for low static and dynamic friction coefficients (μs = 0.4 μd = 0.3), similar to human joints [31]. A force vector was created on the femur distal, as shown in Figure 9, and the value set at 800 N, while defining it by a step function (A, x0, h0, x1, h1).

Simulation of the Multibody Model at the Position of 25°
The simulated model is shown in Figure 9. The graph of rotation against flexion is illustrated in Figure 10. It was noticed that the angle of rotation varies linearly with respect to the flexion. It was observed that the sudden increase in the rotational angle was offset by an increase in the angle of curvature. A sharp rise in the flexion angle till 25° was seen beyond the rotational angle of 20°. The maximum elevation of flexion against the angle of rotation was found to be 25°.

Simulation of the Multibody Model at the Position of 50°
The simulated model and the graphical representation of the rotational angle with flexion are depicted in Figures 11 and 12, respectively. It was observed that there was an

Simulation of the Multibody Model at the Position of 25 •
The simulated model is shown in Figure 9. The graph of rotation against flexion is illustrated in Figure 10. It was noticed that the angle of rotation varies linearly with respect to the flexion. It was observed that the sudden increase in the rotational angle was offset by an increase in the angle of curvature. A sharp rise in the flexion angle till 25 • was seen beyond the rotational angle of 20 • . The maximum elevation of flexion against the angle of rotation was found to be 25 • .

Simulation of the Multibody Model at the Position of 25°
The simulated model is shown in Figure 9. The graph of rotation against flexion is illustrated in Figure 10. It was noticed that the angle of rotation varies linearly with respect to the flexion. It was observed that the sudden increase in the rotational angle was offset by an increase in the angle of curvature. A sharp rise in the flexion angle till 25° was seen beyond the rotational angle of 20°. The maximum elevation of flexion against the angle of rotation was found to be 25°.

Simulation of the Multibody Model at the Position of 50°
The simulated model and the graphical representation of the rotational angle with flexion are depicted in Figures 11 and 12, respectively. It was observed that there was an

Simulation of the Multibody Model at the Position of 50 •
The simulated model and the graphical representation of the rotational angle with flexion are depicted in Figures 11 and 12, respectively. It was observed that there was an increase in the degree of rotation (virtual axis) from 0 to 9.5 • , whereas the flexion (horizontal axis) varies until 50 • . The flexion angle's relative variation was five times bigger than the rotational angle, which indicates that the onset of sliding between the tibia and the femur in the knee occurs in the range of 20-30 • of flexion angle, which conforms with the movement range of the normal human knee.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 12 of 18 increase in the degree of rotation (virtual axis) from 0 to 9.5°, whereas the flexion (horizontal axis) varies until 50°. The flexion angle's relative variation was five times bigger than the rotational angle, which indicates that the onset of sliding between the tibia and the femur in the knee occurs in the range of 20-30° of flexion angle, which conforms with the movement range of the normal human knee.

Simulation of the Multibody Model at the Position of 100 Degrees
The pictorial representations of the model at the angles of 0° and 100° of flexion are shown in Figures 13 and 14. The relationship between the angles of rotation and flexion is illustrated in Figure 15. It is observed that the sudden increase in the rotational angle was offset by an increase in the angle of curvature. A sharp rise in the flexion angle till 35° was seen beyond the rotational angle of 20°, which indicates the onset of sliding between the tibia and the femur in the knee. Similarly, the flexion angle in the range of 20-30° originates the joint prone to rolling. In contrast, the rotational angle's stability for the flexion angle lies in the range of 30-110°. The increasing flexion angle indicates that the tendency of sliding is predominant.  increase in the degree of rotation (virtual axis) from 0 to 9.5°, whereas the flexion (horizontal axis) varies until 50°. The flexion angle's relative variation was five times bigger than the rotational angle, which indicates that the onset of sliding between the tibia and the femur in the knee occurs in the range of 20-30° of flexion angle, which conforms with the movement range of the normal human knee.

Simulation of the Multibody Model at the Position of 100 Degrees
The pictorial representations of the model at the angles of 0° and 100° of flexion are shown in Figures 13 and 14. The relationship between the angles of rotation and flexion is illustrated in Figure 15. It is observed that the sudden increase in the rotational angle was offset by an increase in the angle of curvature. A sharp rise in the flexion angle till 35° was seen beyond the rotational angle of 20°, which indicates the onset of sliding between the tibia and the femur in the knee. Similarly, the flexion angle in the range of 20-30° originates the joint prone to rolling. In contrast, the rotational angle's stability for the flexion angle lies in the range of 30-110°. The increasing flexion angle indicates that the tendency of sliding is predominant.

Simulation of the Multibody Model at the Position of 100 Degrees
The pictorial representations of the model at the angles of 0 • and 100 • of flexion are shown in Figures 13 and 14. The relationship between the angles of rotation and flexion is illustrated in Figure 15. It is observed that the sudden increase in the rotational angle was offset by an increase in the angle of curvature. A sharp rise in the flexion angle till 35 • was seen beyond the rotational angle of 20 • , which indicates the onset of sliding between the tibia and the femur in the knee. Similarly, the flexion angle in the range of 20-30 • originates the joint prone to rolling. In contrast, the rotational angle's stability for the flexion angle lies in the range of 30-110 • . The increasing flexion angle indicates that the tendency of sliding is predominant. Appl. Sci. 2021, 11, x FOR PEER REVIEW 13 of 18

Experimental Measurement Result for Test Machine of Szent Istvan University
The result shown in Figure 16 is essential as it is the basis for our numerical experimentation. The research team of Szent Istvan University developed several prosthesis design methods initiated by [12] using the test machine they developed, as shown in Figure 17. It was mentioned by Balassa, 2019, that the presented results in his study was the best with respect to the closeness to the real knee movement.

Experimental Measurement Result for Test Machine of Szent Istvan University
The result shown in Figure 16 is essential as it is the basis for our numerical experimentation. The research team of Szent Istvan University developed several prosthesis design methods initiated by [12] using the test machine they developed, as shown in Figure 17. It was mentioned by Balassa, 2019, that the presented results in his study was the best with respect to the closeness to the real knee movement. Figure 16. A graph of rotation against flexion for the knee prosthesis designed by Balassa and tested with their machine as compared with the regular movement of the human knee [12]. In order to validate the results from the numerical studies, we compared our virtual numerical model with the prostheses joints that have been tested using the test machine in Szent Istvan University. The average values ( Figure 18) are plotted together against the virtual numerical model, as shown in Figure 19. The close similarity between the two curves indicates that this virtual model can replace the Szent Istvan University test machine's measurement. In order to validate the results from the numerical studies, we compared our virtual numerical model with the prostheses joints that have been tested using the test machine in Szent Istvan University. The average values ( Figure 18) are plotted together against the virtual numerical model, as shown in Figure 19. The close similarity between the two curves indicates that this virtual model can replace the Szent Istvan University test machine's measurement.
In order to validate the results from the numerical studies, we compared our virtual numerical model with the prostheses joints that have been tested using the test machine in Szent Istvan University. The average values ( Figure 18) are plotted together against the virtual numerical model, as shown in Figure 19. The close similarity between the two curves indicates that this virtual model can replace the Szent Istvan University test machine's measurement. Figure 18. A graph of Rotation against Flexion for the prosthesis developed with a test machine [12].
It was also found that there was a rise in the angle of rotation as the flexion angle varied from 0 to 30°. A good agreement of value obtained from our model with other prosthetic joints is established in the flexion range of 30-110°.
We noticed the same results for the proposed model and other prostheses. In this case, we were able to create a new model that enables us to make multiple measurements. We could change the materials made of artificial joints, make them more accessible, and obtain faster results without many calculations. A comparison of our model with the test machine is presented in Table 2.

FLEXION (°)
Average measureme nts of knee prosthesis for balasa Figure 18. A graph of Rotation against Flexion for the prosthesis developed with a test machine [12].

Virtual Model Test Machine
It saves time and material Time-consuming and involves wastage of material It saves money because it is objective focused High cost incurred because it is a procedure of trial and error The lack of transition in flection and rotation is used to simplify the geometry and diversity in motion, allowing our new virtual model to be more realistic.
This process cannot be applied with the test machine, hence the inefficiency of the procedure.

Conclusions
It can be said that the application of the multibody model saves time as there is no involvement of the tibia and femur, as needed for the knee prosthesis. More importantly, as the application of the test machine is omitted in our process, our model's approximations to a human knee are carried out directly. Without cost, we can make several measurements for a person's knee size to be repaired.
The MSC.ADAMS programme was applied to determine the human knee joint's movement in terms of rotation and flexion. The relationship between these two processes was also described. The changes between the condyles of the developed multibody of the prosthesis are also investigated concerning the flexion angle ranging from 20-120°. The boundary conditions were determined, and simulations performed using the ADAM's programme. Three-dimensional geometry was applied in the new virtual model, taking into account the influence of the condyles and collateral. The multibody modelling was used to measure the degree of flexion and rotation of the knee concerning its position, like extension, flexion or rotation, and insert a spring between the tibia and femur while observing its effects on the performance of the knee.
A slip ration, which is higher than 0.45, was achieved as was the limit in literature. Applying our model, an average value of 0.7 was reached, with the maximum reaching up to 0.79 and obtaining an angle between 110-120° for the flexion angle. The generated virtual model was used to measure the knee prosthesis size before its creation. It was also found that there was a rise in the angle of rotation as the flexion angle varied from 0 to 30 • . A good agreement of value obtained from our model with other prosthetic joints is established in the flexion range of 30-110 • .
We noticed the same results for the proposed model and other prostheses. In this case, we were able to create a new model that enables us to make multiple measurements. We could change the materials made of artificial joints, make them more accessible, and obtain faster results without many calculations. A comparison of our model with the test machine is presented in Table 2. Table 2. Comparison between the virtual model procedure and the test machine.

Virtual Model Test Machine
It saves time and material Time-consuming and involves wastage of material It saves money because it is objective focused High cost incurred because it is a procedure of trial and error The lack of transition in flection and rotation is used to simplify the geometry and diversity in motion, allowing our new virtual model to be more realistic.
This process cannot be applied with the test machine, hence the inefficiency of the procedure.

Conclusions
It can be said that the application of the multibody model saves time as there is no involvement of the tibia and femur, as needed for the knee prosthesis. More importantly, as the application of the test machine is omitted in our process, our model's approximations to a human knee are carried out directly. Without cost, we can make several measurements for a person's knee size to be repaired.
The MSC.ADAMS programme was applied to determine the human knee joint's movement in terms of rotation and flexion. The relationship between these two processes was also described. The changes between the condyles of the developed multibody of the prosthesis are also investigated concerning the flexion angle ranging from 20-120 • . The boundary conditions were determined, and simulations performed using the ADAM's programme. Three-dimensional geometry was applied in the new virtual model, taking into account the influence of the condyles and collateral. The multibody modelling was used to measure the degree of flexion and rotation of the knee concerning its position, like extension, flexion or rotation, and insert a spring between the tibia and femur while observing its effects on the performance of the knee.
A slip ration, which is higher than 0.45, was achieved as was the limit in literature. Applying our model, an average value of 0.7 was reached, with the maximum reaching up to 0.79 and obtaining an angle between 110-120 • for the flexion angle. The generated virtual model was used to measure the knee prosthesis size before its creation.
Finally, this virtual model could be used to measure the knee prosthesis size before creating it because it saves time, money, and effort instead of using 3D printing technology, and CNC milling consumes our time, money, and effort.
In future work, we will try to create a virtual model for the ankle to obtain a new virtual model for the complete human leg with the knee and ankle and all required movements. Analysis of the anatomical angles, such as the different rotation as human full legs, abduction, and adduction, will be conducted.