Next Article in Journal
Comprehensive Performance Evaluation of an Unmanned Excavator Based on Signal Stability Extraction
Previous Article in Journal
Image-Based Tool Characterization and DEM Simulation of Abrasive Brushing Processes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Optimization of a Robot Dosing Device for Aliquoting of Biological Samples Based on Genetic Algorithms

1
Research Institute Robotics and Control Systems, Belgorod State Technological University n.a. V.G. Shukhov, Kostyukova st., 46, 308012 Belgorod, Russia
2
Department of Mechanical, Energy and Management Engineering, University of Calabria, 87036 Rende, Italy
*
Author to whom correspondence should be addressed.
Submission received: 19 January 2024 / Revised: 26 February 2024 / Accepted: 27 February 2024 / Published: 1 March 2024
(This article belongs to the Special Issue Digital Twins in Robotics)

Abstract

:
Aliquoting of biological samples refers to the process of dividing a larger biological sample into smaller, representative portions known as aliquots. This procedure is commonly employed in laboratories, especially in fields like molecular biology, genetics, and clinical research. Currently, manual dosing devices are commonplace in laboratories, but they demand a significant amount of time for their manual operation. The automated dosing devices available are integrated into narrowly focused aliquoting systems and lack versatility as manipulator equipment. Addressing this limitation, a novel technical solution is proposed in this paper for a modular dosing device compatible with robotic manipulators. The paper introduces and details a mathematical model, optimizes its parameters, and constructs a detailed 3D model using the NX environment to demonstrate the engineering feasibility of our concept. It further outlines the development of a three-dimensional dynamic simulation model for the dosing device, comparing analytical calculations with simulation results. The construction of a dosing device prototype is discussed, followed by a comprehensive experimental validation.

1. Introduction

The last few years have shown a sharp increase in demand for intelligent robots capable of performing complex tasks without human intervention in almost every industry. The biomedicine industry is no exception. Compliance of biological samples with high international standards for laboratory and clinical research can be achieved only by robotization, reducing the duration of preparation and storage of initial samples at room temperature, minimizing the number of manual labor operations, and improving the accuracy of the considered characteristics of samples [1]. The need to reduce the number of production processes involving manual labor is also important to prevent contamination of equipment and work surfaces with infectious agents and, consequently, minimize the risk of infection of personnel.
The need for dosing devices, which are installed on automated manipulators and analyzers, is steadily growing. In modern research and clinical laboratories, sample preparation is becoming the most critical stage in experimental processes, especially when it comes to highly effective methods. The article [2] discusses various pipetting robots used for sequencing the human genome. In ref. [3], a robot has been developed with the ability to place from 1 to 10 pipettes, which allows users without programming skills to easily automate a wide range of experiments. In ref. [4], an automated pipette printed on a 3D printer is presented, designed for conducting chemical experiments in deserted laboratories.
There is a wide range of automated dosing or pipetting systems on the market today. In the article [5], the pipetting technology is considered, and the performance of working with liquids is monitored. In the research [6], a new immunoassay analyzer based on pipetting using a removable pipette tip with a magnetic ring was developed. This is in demand in pharmaceutical and biological laboratories. There is an invention [7], which represents a pipetting device that is highly flexible in use and has an integrated and compact design of independent pipettes, which allows aliquoting of several cells at once. In some areas, the key characteristic of a dosing device is not speed but accuracy [8,9,10]. These articles present automatic devices for dosing nanoliter droplets with high accuracy. The liquid is aspirated using a micro-pump and an electromagnetic valve. In ref. [11], a capillary chip was developed that automates aliquoting of the sample and eliminates the need for laboratory micropipettes due to the integration of CAC into the chip. The article [12] discusses an innovative method of passive aliquoting of samples based on pressure equilibrium, which provides accurate self-correcting aliquoting without any valves. There are tasks of liquid dosing where it is not individual approach or accuracy that is important, but high speed and precision of the aliquoting operation. Thus, in ref. [13], a new concept of automated sample processing for future proteomic and clinical studies is outlined. To save and archive samples of clinical patients, test tube systems were used with simultaneous dosing of tablets in the 96 and 384 cell formats, while the authors note that the mass character of this approach made it possible to complete the task of 11,000 aliquots of samples by 99.8%, while the quality of the samples remained at the same level. To improve the quality of samples and reduce errors, there are devices that perform automatic dosing or dispensing of the drug in accordance with the task [14]. In their investigation, the authors managed to reduce the overall frequency of unintended dosing errors from 0.015% to 0.002% using automatic dosing. Automated dosing devices find their field of application in other fields other than medicine and healthcare. Thus, the authors in [15] used an analyzer pipetting biological fluid samples to determine the molar ratio of magnesium and calcium (Mg/Ca) in foraminifera calcite, which is widely used to reconstruct past seawater temperatures.
At the moment, manual dosing devices are common, work that takes 60–70% of the time of an ordinary employee, for example, in laboratories of microbiology and virology. The design of a dosing device (DD) with an automatic tip removal mechanism is known [16]. This device is intended for manual use; however, the tip change mechanism may also be of interest for automated systems. The tip is replaced by a movable cone, which translates so that the tip contacts the fixed sleeve and detaches from the cone. In the article [17], a method for adapting a manual dosing device for installation on a manipulator using a specially made housing is proposed. The manual dosing device specified in the work takes only a specifically configured volume of liquid, does not have a flexible volume setting, and is often not universal. Due to the lack of manipulator equipment designed for aliquot dosing, there is a need to develop new technical solutions for modular dosing devices. Thus, the purpose of this work is to develop a new universal dosing device (DD). The systematic design of the proposed concept is carried out using the latest digital design tools, as well as multi-objective optimization methods based on genetic algorithms and particle swarm optimization. The main contributions are listed below:
(1)
A new design of a dosing device for aliquoting biological material has been proposed. The device is self-sufficient and can be installed on any manipulator thanks to a simple control system and connection. The motor input is around 1.2 W with a 5 V and 240 mA power supply. The gearbox used allows you to perform dosing operations, such as fixing/resetting the tip, when using the same engine. These technical solutions make the device lightweight and versatile. Cleaning of the device is not required, as the liquid is drawn only into the tip and does not come into contact with the internal pneumatic system. This device has no available analogues among the manipulator equipment and is in demand for automated dosing.
(2)
An approach to designing a dosing device using multi-objective optimization methods is presented. The dosing device was described in detail using well-known mathematical expressions, on the basis of which a set of Pareto was constructed. The set contains optimal configurations of geometric parameters, based on which the most preferred one is selected. This method allows you to both scale the resulting model and use it for different source data.
This document is organized as follows: in Section 2, the requirements for the parameters and design of the dosing device are formulated, the DD design scheme is proposed and described, and the mechanisms that need to be calculated are highlighted; in Section 3, a mathematical model of the dosing device is presented, the tip fixation mechanism is calculated, the pneumatic part is calculated, and the algorithm for determining the required engine torque is described; in Section 4, the dependences of the geometric parameters of the device and the performance of DD are presented, as well as the optimization of these parameters using a genetic algorithm; Section 5 describes the development of a 3D dynamic simulation DD model and also contains a comparison of the results of the analytical calculation and the results of the simulation; Section 6 is devoted to the finite element analysis of the DD rod under dynamic conditions, which was destroyed during the first tests of the layout, as a result of which a more thorough analysis was required; Section 7 describes the design of the prototype DD and experimental research with the results provided; and Section 8 formulates a conclusion on the work performed.

2. Requirements for the Parameters of the DD and the Proposed Design

Working in biological laboratories involves aliquoting the biological fluids of various patients. Aliquot dosing takes place using replaceable disposable tips, which must be changed after each dosage. At the moment, dosing with the help of manual dosing devices prevails, where the operator does all the work. Automated dosing devices, in the vast majority, are part of a robotic system, have a high price, and cannot be used in other robotic systems. To solve this problem, a modular dosing device that can be installed on any manipulator and is capable of performing all necessary operations is proposed. The design of a robotic system (RS) for aliquoting biological fluid with a modular dosing device is shown in Figure 1.
The proposed robotic system (RS) includes housing 1, in which there is a sterile work area 2 with clean tips 3, a tripod with small volume tubes 4, and a test tube with biological fluid 5. The manipulator 6 is fixed on the base 7 and positions the dosing device 8. The biological fluid is divided into two fractions, blood serum 9 and clot 10. Due to the algorithms of the technical vision system 11, the boundary between fractions and the coordinates of the blood serum is determined, after which the manipulator, using a dosing device, performs sampling into the tip and dosing of serum into a tripod with small volume tubes.
The following requirements for the dosing device were compiled to implement automated dosing in the process of blood aliquotation:
  • The operations that the dosing device (DD) must perform are fixing the pipette, taking biological fluid, draining the liquid in equal parts, and dumping the tip.
  • The DD must be able to work with standard laboratory tips, the manufacture of which is regulated by the ISO 13485 CE 100UL certificate [18].
  • The parameters of the tip are most often used for blood aliquoting: volume up to 1 mL, diameter of the larger hole D t i p = 7   m m , diameter of the smaller hole 0.7 mm, height 60 mm, material of manufacture plastic, without filter element.
  • The weight of the DD should be minimized, which will allow for low-load work. DD tests are planned to be carried out on robots available in the laboratory, with a load capacity of 0.5, 3 and 6 kg. Accordingly, the maximum mass of DD is selected to be no more than 0.5 kg.
  • The ability to dose liquid with a density of up to 1110 kg/m3, which corresponds to the maximum blood density.
  • The possibility of manufacturing a dosing device using 3D printing.
To meet these requirements, consider the following technical solution DD for performing operations: automatic fixing of the tip, liquid dosing, discharge of the tip; see Figure 2. DD includes the housing assembly 1, in which the motor 2 is installed, connected to the drive gear 3, due to the motor coupling 4, which ensures the rotation of the external gear wheel 5 and cam 6. When the cam is rotated, the rod 7 is displaced, which deforms the elastic O-ring 8. In the initial state, the diameter of the sealing ring is equal to D1 and is smaller than the diameter of the tip 9, and in the deformable state, D2 is larger than the diameter of the tip 9, which allows for the capture and sealing of the tip 9. The rod is connected to the cylindrical cam by means of a lever 10. The gears 11 are rigidly fixed to the corresponding screw of the piston 12, each of which rotates around its axis and are is connected to the piston 13. The pistons are translationally moved in the chamber 14 of the pneumatic cylinders, which provides a discharge in the vacuum line 15 and allows for the intake of liquid into the sealed tip 9.
Since the rubber seal must be placed in the tip, it has the following limitations: the outer diameter in the initial and deformable state is D 1 < D t i p < D 2 mm, and the inner diameter of the seal, D3, must be larger than the diameter of the rod. Based on the standard ISO 3601 [19] solutions available, the most suitable is a rubber seal with an outer diameter of 6.5 mm and a cross-section diameter from 0.1 to 3 mm. The parameters of the rubber seal D1 = 6.5 mm and d1 = 1 mm are selected, which will allow the sealing ring to be installed on a rod with a diameter of 1 mm.
In the described mechanism, the motor drives the pneumatic system as well as the rod, which deforms the O-ring. The DD should have a small mass and compact dimensions, which means that it is necessary to design considering the minimization of dimensions as well as to prevent the installation of an engine with excessive power and weight. Thus, it is necessary to determine the minimum torque to ensure the operation of all DD mechanisms. Let us determine the required torque for the engine.

3. Mathematical Model of the Dosing Device

The mechanism of fixing the tip is the most loaded, since deformation of the O-ring to achieve a rigid fixation of the tip can occur with considerable effort. To determine the torque required for the operation of the tip fixation mechanism, we will conduct a kinetostatic analysis.
When the rod moves, the rubber seal must change the diameter; see Figure 3a. For uniform deformation, the rod shape is made in the form of a cone, which, when moving forward, causes an increase in the radius R (deformation) of the O-ring. Let us make up the task condition. Assume that the O-ring is made of a rubber cord with a length of l 0 and a mass of m. The ring is put on a smooth cone with an angle of 3α at the top; see Figure 3b,c. Let us determine the tension force T arising in the ring. To do this, select an infinitesimal element of the cord with length Δ l and mass Δ m , which are determined by the angle β. Gravity Δ m g acts vertically downwards on it, the resulting elastic force F e l a s t is directed from it to the center of the circle, and the reaction force of the support N is normal to the surface of the cone; see Figure 3.
Under the condition of the equilibrium of the system, the projection of forces on the X and Y axes is determined by the following expression:
F e l a s t x N cos α = 0 N sin α m g = 0 ,
where F e l a s t x = m g tan α .
Let us assume that at small angles, sin β / 2 = β . The angle β can be determined from the following expression: β = l R , and R = l 0 + x 2 π , where x is the relative change in length. Let us write down Hooke’s law in relative form:
x = T l 0 E S
Then, we express and define the tension force T taking into account the assumption ∆ m l = m l 0 ; then, we define T as follows:
T = m g E S E S 2 π tan α m g
where E is the Young’s modulus for the O-ring material [20], S is the cross-section of the cord. Let us express the friction force, provided that the rod moves inside the O-ring without rotation. Then the expression takes the following form:
F B = 2 π μ m g E S E S 2 π tan α m g
where μ is the coefficient of friction between the material of the sealing ring and the material from which the rod is made, rubber and steel are used in the calculation. The design scheme of the rod lifting mechanism is shown in Figure 4a. The mechanism works as follows: the lever CA is fixed on a rotary joint A and has a slider B, which in turn is fixed on a rod that deforms the O-ring. The diagram of the stem with an O-ring is shown in Figure 4b. The lever CA rotates at an angle γ due to sliding along the profile of the cylindrical cam at point C, which leads to a vertical displacement of the rod.
To determine the required torque on the cylindrical cam to lift the rod to a given height h 1 , we determine the height of the profile h 2 and the force F C , taking into account the difference in the lengths of the shoulder l 1 and the total length of the lever l 2 using the following expression:
sin α = h 2 l 2 sin α = h 1 l 1 ,   t h e r e f o r e   h 2 = h 1 · l 2 l 1 ,   F B = F C · l 2 l 1 ,
In order to determine the dependence of the force F C and the torque M1 required to achieve it, let us assume the surface of a cylindrical cam in the form of an expanded inclined surface on which the slider C slides; see Figure 5.
To turn the cam, it is necessary to apply torque M1, according to the formula M 1 = M 1 . f + M P . Where M 1 . f is the torque of friction force on the support end, M P is the torque of force on the cam profile M 1 . f = R C F C   f , and f is the reduced coefficient of friction.
According to the well-known theorem of mechanics, which considers friction forces, the slider is in equilibrium if the resultant Fn of the system of external forces is deviated from the normal n − n by the friction angle ϕ. In our case, the external ones are the axial force F C and the circumferential force F t = M 1 / R C . From the force triangle in Figure 5b, we can find the unknown force F t and determine the desired torque:
M P = F C R C tan φ + ψ ,
where ψ is the angle of lift of the cam helical surface, φ is the angle of friction of the cam helical surface. The moment of rotation of the gear wheel is expressed and defined as follows:
M 1 = F C R C f + F C R C f + L h 2 1 f L h 2 ,
Let us determine the force to be applied to the piston to lift the liquid column in the pipette to the required height, considering atmospheric pressure; see Figure 6.
The required torque of the screw for the movement of the pneumatic cylinder rod is determined similarly to the expression (6), then
M 2 = ( F 1 + F 3 F 2 ) tan ( φ + ψ P ) R z 3
where ψ P is the angle of the thread on the piston rod ψ P = tan 1 S p π d 2 , w h e r e   S p , is the pitch of the screw, d 2 is the diameter of the threaded rod, D P is the diameter of the piston, and F 1 is the force acting on the piston from the threaded rod, as a result of atmospheric pressure F 1 = P a π ( D P 2 d 2 2 ) 4 . Where P a is normal atmospheric pressure, F 2 is the force acting on the piston from the side of the vacuum line, which is determined by the following expression F 2 = ( P a ρ g H ) π D P 2 4 . Where ρ is the density of a column of liquid that rises to a height H, F 3 is the friction force in the piston seals F 3 = f p π D P l p P 0 , where f p —the coefficient of friction accepted for rubber cuffs; l p —width of the contact surface of the seal; P 0 —pressure in the pneumatic cylinder.
The required engine torque is defined as the sum of the rotational torques of the corresponding gears, taking into account their gear ratio.
M 3 = M 1 z 3 z 1 + 2 M 2 z 3 z 2
where z i is the number of cogs of the corresponding gears. Considering the fact that tan ( tan 1 f + tan 1 π d 2 S p ) = f + π d 2 S p 1 f π d 2 S p , the torque M 3 is equal to
M 3 = F C R C f + F C R C f + L h 2 1 f L h 2 + 2 ( F 1 + F 3 F 2 ) R z 3 tan ( φ + ψ P ) z 2 z 3 ,
Expression (9) allows you to determine the maximum torque of the engine, considering all of the design parameters. However, some of them may vary within certain limits. So, the diameter and height of the piston stroke can change, leaving the torque unchanged. To determine the most preferred configuration, we will optimize the geometric parameters of the DD.

4. Optimization of Geometric Parameters of the Mechanism

First, we obtain simple expressions that define the relationship between the geometric parameters of the device and its performance. We assume that, considering the gear modulus m = 2.5 , the linear displacement of the piston per turn of the engine l = z 1 z 3 S p , where z 1 ,   z 3 is the number of cogs of the gears and, S p is the pitch of the screw.
To take the required volume of liquid V with a piston diameter D P , the stroke of the rod l is required, and this dependence allows you to obtain the required geometric parameters of the mechanism for taking a given volume of liquid equal to 1000 mm3. Let us take into account only the working stroke of the mechanism n θ without taking into account the fixation V = π D P 2 4 l · n θ , and substitute the required number of cogs,
z 1 = V 4 z 3 S p π D P 2 n ,
where n is defined as n θ = 1 θ 360 ° . The volume of the tip is known as V = 1000   m m 3 , while the metering device must be compact, respectively. For the smallest gear wheel, the minimum number of cogs z 3 = 6 is selected, and the diameter of the piston D P may vary. The number of required cogs z 1 depends only on the thread pitch S p of the rod.
The dividing diameter of the gear wheel D z 1 is determined by the following:
D z 1 = 4 m V z 3 S p π D P 2 n + 2 m .
Thus, it is possible to determine the maximum dimensions of the gears and draw a conclusion about the dimensions of the metering device; however, a number of parameters may be different.
To determine the values of the design parameters, we use multi-objective optimization methods. To do this, we will find many solutions with different design parameters for the mechanism under given constraints. We will find out from the whole set those that will have the minimum values: the torque M2, the outer diameter of the case D o u t , and the height of the dosing device H o u t .
The formulation of the optimization problem involves the selection of criteria and constraints on the parameters. Let us consider the optimization criteria:
  • The torque of the engine (motor) tends to a minimum of M → min. Using Expression (9) we obtain the dependence of the required torque on the motor shaft on the design parameters of the mechanism to ensure rotation throughout the operation.
  • The outer diameter of the housing D o u t m i n , which is determined by the following formula:
    D o u t = D z 1 + 2 D b + H z + 2 c
    where D b is the diameter of the bearings of the outer gear, adopted by 6 mm since when 3D printing a groove for bearings of a smaller diameter, it was not possible to obtain a high-quality surface without using supports; therefore, c is the thickness of the walls of the housing, adopted by 1 mm for ease of calculation.
  • The height of the dosing device H o u t m i n , which consists of the length of the piston rod and its chamber, the lengths of which are assumed to be equal, as well as the thickness of the gears and the walls of the housing,
    H o u t = 2 l + H z + 3 c
    where H z , the thickness of the gear wheel adopted is 6 mm, as the minimum for the assembly of the bearing and as the minimum for assembling a bearing, which has an internal cage with a toothed profile and will also reduce the dimensions.
Let us consider the optimization parameters: the thread pitch on the piston rod S p , the diameter of the piston D P , and the angle of rotation θ of the external gear wheel to raise the rod to a predetermined height. Optimization limitations include the following: the thread pitch on the piston rod S p should not exceed 4 mm since the dosing accuracy will drop significantly as well as the required torque, which will increase many times; the minimum thread pitch should not be less than 1 mm since 3D printing will be difficult; and the movement of the rod per engine revolution is very small; therefore, 1 < S p < 4 mm. The diameter of the piston cannot be larger than the radius of the outer gear, as this will lead to an increase in the dimensions of the DD, and a very small piston size will require a large stroke of the piston rod, which will also increase the dimensions; then, we will limit the diameter of the piston within d 2 < D P < D z 1 2 . The angle of rotation θ of the external gear wheel to raise the rod to a given height can vary from 0° to 360° degrees; however, at 0°, the angle of elevation of the plane along which the lever slides will be vertical, and at 360° degrees, the angle is smooth, but the tip is fixed until it is completely rotated, which is unacceptable. Therefore, we limit the fixation of the tip to a quarter turn of 0° < θ < 90°.
To reduce the number of criteria and thereby build a Pareto set more efficiently, we will perform the convolution of criteria 2 (maximum diameter D o u t ) and 3 (maximum height of the device H o u t ) into one criterion, which is the volume V o u t occupied by the excavating device.
V = π D o u t 2 4 H o u t
As optimization algorithms, we use evolutionary algorithms that are widely used to solve similar problems [21,22,23]. One of the solutions is the “gamultiobj” module, which is included in the MATLAB 2019 software package and allows the use of a genetic algorithm to solve the problem of multi-objective optimization, including for shapes and surfaces, as shown in a number of studies [24,25,26]. To ensure the reliability of the result, two modifications of the multi-objective optimization (MO) algorithm developed earlier by the authors [27] were used, one of which is based on the parallel GWO algorithm (MO-GWO), and the second is based on the parallel PSO algorithm (MO-PSO).
Let us set the initial data for a computational experiment. The values of the constants used in the simulation are shown in Table 1. The parameters D t i p and H are taken in accordance with the parameters of standard tips; the friction coefficients μ, f p , and f , as well as the Young’s modulus E are selected in accordance with the selected materials of manufacture of all elements; the angle at the tip of the cone α° is selected experimentally; and the number of cogs z 2 is selected as the largest possible for the selected design.
Optimization parameters ranges: S P , ∈ [0.0001; 0.0025], D P ∈ [ D P ; 0.5], θ ° ∈ [0.00001; π/4].
The algorithm parameters were selected considering an approximately equal number of total iterations for the correct comparison of algorithms. Parameters of the NSGA-II algorithm of the «gamultiobj» module: the number of individuals in the population H G A = 10,000, the number of generations W G A = 100. The total number of iterations of calculating the criterion function as a result is
N G A = H G A W G A = 10 6
The parameters of the MO-GWO algorithm are as follows: the number of stages, the number of grid nodes for each criterion for sorting the importance coefficients n G W O = 16 = 16, the number of stages for each grid node S = 2, the number of individuals in the population at each stage H G W O = 390, and the number of generations at each stage W G W O = 5. So, the total number of iterations is the number of iterations.
N G W O = H G W O W G W O S n G W O 2 + N I = 998400 + N I
where N I is the number of iterations to calculate the average value of the coefficients, not exceeding 1600 iterations.
The optimization of the MONO and MOGWAI algorithms is performed using a software module in the C++ programming language. Parallel computing is implemented using the OpenMP library. The visualization of the results was performed using a Python script (Matplotlib and JSON libraries). A computer with a six-core processor and a frequency of 3.79 GHz was used for the simulation.
The calculation time and the proportion of approximation points in the resulting Pareto front obtained by each of the algorithms are used as metrics of the effectiveness of the algorithms. The optimization was performed 10 times. Table 2 shows the number of Pareto points and the proportion obtained using each of the algorithms. The average calculation time for the NSGA-II algorithm was 43 s, for MO-GWO it was 9 s, and for MO-PSO it was 3 s.
Figure 7 shows an example of the Pareto set obtained by each of the algorithms for experiment 4. The horizontal axis corresponds to the required torque, and the vertical axis corresponds to the volume of the DD.
Figure 8 shows a comparison of each of the Pareto sets and the resulting Pareto set.
As can be seen from the graphs, the NSGA-II algorithm made it possible to obtain a more uniform distribution of the Pareto set, in particular in its central zone. Nevertheless, the resulting Pareto set had an average of 57.9% of the points obtained by the MO-GWO algorithm compared with 33.1% by the NSGA-II algorithm. Due to the significantly higher performance of the MO-GWO and MO-PSO algorithms in comparison with the NSGA-II algorithm, additional computational experiments were performed to compare the efficiency of the algorithms with an equal calculation time. The NSGA-II algorithm’s operating time, which is 43 s, was taken as the total calculation time. To bring the time to an approximately equal value, the parameters of the MO-GWO and MO-PSO algorithms were changed as follows: n G W O = 25 , H G W O = 500   n P S O = 40 , and H P S O = 400 .
With equal computation time (Table 3), the MO-PSO algorithm made it possible to obtain an average share of 54.8% of the points of the resulting Pareto set, MO-GWO achieved 35.8%, and NSGA-II achieved only 9.4%. However, the MO-PSO algorithm, even with equal computation time, did not allow us to obtain a uniform distribution of points in the central part of the Pareto set (Figure 9).
The linear section of the Pareto front corresponds to the parameters of the DD, in which the radius of the outer gear is always larger than the diameter of the piston and is a key factor in determining the volume, while the torque is the lowest as a result of the highest gear ratio; however, there is a curved section on the graph where the diameter of the piston becomes the decisive value for the volume of the DD; in this case, the required torque can increase many times.
The torque of the selected stepper motor is 0.0441 Nm based on the obtained Pareto set, the closest value is 0.0406 Nm, for which the geometric parameters are accepted. The calculation results and the values of the calculated parameters that are not given in the description are presented in Table 4.

5. Building a Digital Twin

Figure 10 shows a detailed 3D model of the calculated dosing device, which was built using the CAD/CAE NX system. The layout of the gears is shown in Figure 10a. The drive gear 3 is coupled to an external gear 5, through which the torque is transmitted to the gears 11 connected to the screw. The cam profile is located on the upper side of the larger gear wheel; see Figure 10b. The use of 3D printing allows you to combine several kinematic functions into one part. The lever based on the cam profile has a folded shape, but its size in the projection exactly corresponds to the parameters obtained during the calculation. The edge of the housing, which is in contact with the tip, is made in a separate part and has a taper, which allows you to center the tip. An analysis of the assemblability of the entire product was carried out, as a result of which an assembly sequence was developed, and complex parts were divided into several parts for 3D printing.
The resulting, detailed 3D model is used for simulation and finite element analysis. The NX Nastran Advanced Kinematics module was used for this purpose. Each of the movable assembly elements is designated as a “link”, in the menu of which the main mass and inertial parameters are determined. The conjugated “links” form “nodes” and “connections”, the menu of which allows you to adjust degrees of freedom, friction forces, limits, and laws of motion, which are described in more detail in a number of works [28,29,30,31]. The coupling of the lever with the cam surface is realized using the “3D Contact” command, with stiffness parameters of 10,000 N/mm and friction parameters selected in accordance with the analytical calculation. The rubber seal in the simulation model was replaced with a spring with stiffness parameters of 0.0079 Nm and deformation or elongation of 2 mm. The force acting on the piston as a result of creating increased or decreased pressure is defined as a force vector perpendicular to the piston surface, the value of which is the difference between F1 and F2 of the mathematical model. As a result, a dynamic simulation model of the dosing device was obtained, and the results were obtained. Figure 11 shows a graph of the required torque DD, the peak value of which is reached at 0.54 s of operation of the mechanism and is equal to 0.03884 Nm. At this torque, the lever of the metering device overcomes the boundary between the rounded and rectilinear sections of the pushing profile. Further rotation of the mechanism requires less torque since the pushing profile changes the position of the lever evenly along the rectilinear inclined section of the surface of the pushing profile.
Figure 12 shows a graph of the movement of the rod, which is equal to 1.96 mm, which corresponds to the analytically specified ratio of lifting the lever to turning the cam. At the same time, the greatest torque ac of the engine is observed at the very beginning of the movement of the rod.

6. Finite Element Analysis

Compression of the rubber seal can occur at a sufficiently high load. There are CAE (computer aided engineering) analysis methods that are implemented in various software products [32,33]. To test the strength of the ABS plastic rod, a finite element model was built, which is deformed at various natural oscillations from 0.004 Hz to 567,060 Hz; see Figure 13. Thus, the load at one or another node of the CAE grid becomes known during various movements.
The “malleable body” simulation mode is activated. Using the “flexible connection” command, we replace the rod link with its calculated CAE file with the “sim2” extension. Thus, depending on the acting load on the «stem_sim2», the values of its deformation and stresses will be obtained (Figure 14).
The value of the limiting stress for ABS plastic is known as the tensile yield strength 24–40 MPa and the tensile strength 21–50 MPa. Thus, the rod has a much higher margin of safety. This calculation does not consider the heterogeneity of the material that will result from printing on a 3D printer. The presented calculation was performed due to the fact that during the manufacture of the first prototype, the rod was destroyed in the area of the washer installation, after which the model was changed, calculated, and manufactured with the appropriate cross-section; see Figure 15a. The plastic rod withstood the load, but was subjected to excessive wear, including in the area of the O-ring installation. Also, the rod was often subjected to shocks as a result of the erroneous positioning of the robot, and therefore it was made of St 37-2 steel; see Figure 15b.

7. Building a Prototype and Testing

A design of a dosing device using 3D printing has been developed, see Figure 16, for the construction of which relevant scientific works in the field of additive technologies were used [34,35]. A 28BYJ-48.5 V stepper motor is used as a drive, which is connected to a frequency converter.
The DD layout consists of a composite housing, including base 1, cover 2, and top 3. An engine 4 is installed in the housing, connected to the drive gear 5. An engine 4 is installed in the housing, connected to the drive gear 5. The driving gear 5 provides rotation of the external gear 6, which has a pushing profile at the end. When the outer gear is rotated by ten degrees or more, the rod 7 is displaced, which deforms the elastic O-ring 8. In the initial state, the diameter of the sealing ring 8 is equal to D1 and is smaller than the diameter of the fixed tip, and in the deformable D2, it is larger than the diameter of the tip, which allows for gripping and sealing. The rod 7 is coupled to an external gear wheel by means of a washer 9 and a lever 10, the use of which increases the clamping force. The lever 10 has an upper position in the range from 11 to 350 degrees of rotation of the external gear 6. In this range, each gear-screw 11 rotates, which is a single part consisting of a gear to which the screw is perpendicular and coaxially fixed. The forward movement of the pinion screw along the axis of rotation is limited by the stopper 12. This allows, as a result of rotation, to move the piston 13, on which the seal 14 is installed. Each of the pistons 13 translationally moves in the chamber of the pneumatic cylinder 15. The movement of the pistons in the chamber of the pneumatic cylinder provides a discharge in the vacuum line 16, which allows for the intake of liquid into the tip.
The motor rotation is controlled using an Arduino UNO, which is connected to the control panel of the AUBA i5 robot. The power supplied by the robot is controlled by two switches for various tests; see Figure 17a,b.
The control program of the AUBA i5 robot was developed to conduct the experiment. As a result, the robot grabs the tip, positioning it over a test tube with a liquid emulating a biological one. Liquid is taken in the amount of three equal fractions of 2 mL each; the dosing device is positioned above a small volume tube; three fractions are alternately drained into a small volume tube; and the tip is discharged. The experiment has been repeated many times.

8. Results and Discussion

The torque of the selected stepper motor is 0.0441 Nm, which is about 8% higher than the required operating time of the metering device, as a result of which, after 15 min of operation, the engine was significantly heated. It was possible to reduce the load on the engine by installing an O-ring made of softer rubber, which is shown in Figure 15b, mounted on a rod made of metal.
In 19 cases out of 20, the robot coped with fixing the tip and dosing the liquid; however, inaccurate dosing occurred, and contamination of the equipment was recorded (liquid got on the work surface). It was revealed that in the case of sequential immersion and liquid intake, partial liquid intake occurs due to the compression of gases in the vacuum line. Therefore, the intake of the correct volume of liquid is feasible in the following cases: When the tip is immersed to a level of no more than the volume being withdrawn (Figure 18a), it is necessary to synchronize the speed of liquid intake and the movement of the manipulator. It is also possible to correct the error at the software level, considering the proportion of self-collecting liquid.
During the intake of liquid, a micro drop occurs (Figure 18b), which can lead to contamination of the equipment. This can be avoided if, after taking the liquid, a small amount of air is taken, in which case a drop forms inside the tip and is unable to get to the work surface. The volume of air intake is calculated by the number of turns of the engine and is considered when dosing. There are studies on microdosing [36], which describe the modernization of the tip surface to eliminate the formation of micro drops on the surface.
However, our current research work still has some limitations. On the one hand, the resulting prototype is more of an experimental sample and only proves the correctness of the research methods. The metering device must be tested in real conditions during long hours of operation and ensure the required reliability. The question of fixing the measurement error of the dosed liquid also remained open. It is necessary to research feedback methods for dosing small volumes of liquid in more detail. On the other hand, the theoretical analysis and design of the dosing device in this research are carried out within the framework of the formulated requirements for a specialized task. In particular, for dosing biological fluid with a single replaceable tip, which is significantly slower than known laboratory dosing units with multiple tips, this is only the first step in our research. In future work, the dosing device will be studied as part of the aliquoting system of biological material, and a more advanced control system will be presented.

9. Conclusions

In conclusion, this paper presents a comprehensive analysis culminating in the development of a scalable mathematical model for the systematic design of a dosing device (DD). Multi-objective optimization using MO-GWO, MO-PSO, and NSGA-II algorithms was executed, revealing superior performance metrics for MO-GWO and MO-PSO in generating a larger number of optimal Pareto front values. Notably, NSGA-II demonstrated a more uniform distribution of points, constructing unique Pareto zones. The Pareto set analysis facilitated the selection of torque and volume parameters, showcasing superiority over existing stepper motors. The determined engine torque for the metering device stands at 0.0406 Nm. Subsequent 3D modeling and simulation, utilizing geometric parameters obtained during calculations, exhibited a small discrepancy of 4.4% from the analytical results, affirming the model’s reliability. Conducting a finite element analysis rectified initial oversights in the dynamic conditions of the DD rod. The manufactured ABS plastic prototype, addressing joint tightness during O-ring deformation, further integrated hardware successfully. Experimental tests, emulating biological fluid scenarios, validated the effectiveness of the redesigned DD, demonstrating its capability to fix, pipette, and eject the tip using a single motor. In essence, the proposed tip-fixing method emerged as efficient, ensuring simultaneous alignment, sealing, and dependable fixation. This research contributes valuable insights into the optimization of dosing devices, offering a promising avenue for enhanced precision and functionality in laboratory applications of automatized aliquotation.

Author Contributions

Conceptualization, L.R. and G.C.; methodology, L.R. and A.V.; software, D.M.; validation, L.R. and G.C.; formal analysis, D.M.; investigation, A.V.; resources, L.R.; data curation, A.V.; writing—original draft preparation, A.V.; writing—review and editing, G.C. and L.R.; visualization, D.M.; supervision, L.R. and G.C.; project administration, G.C.; funding acquisition, L.R. and G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the state assignment of Ministry of Science and Higher Education of the Russian Federation under Grant FZWN-2020-0017.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The work was carried out using equipment of High Technology Center at BSTU named after V.G. Shukhov.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. National Cancer Institute. NCI Best Practices for Biospecimen Resources. Biorepositories and Biospecimen Research Branch National Institutes of Health Table of Contents. National Cancer Institute. 2016. Available online: https://www.juntadeandalucia.es/salud/biobanco/sites/default/files/users/user29/2016-ncibestpractices.pdf (accessed on 20 December 2023).
  2. Houriiyah, T.; San, J.; Jennifer, G.; Tulio, O. Unlocking the efficiency of genomics laboratories with robotic liquid-handling. BMC Genom. 2020, 21, 729. [Google Scholar]
  3. Philip, D.; Tobias, K.; Geethika, A.; Nouraiz, A.; Yang, Z.; Florin, S.; Arne, W.; Daniel, S.; Shunsuke, K.; Dirk, L.; et al. Open-source personal pipetting robots with live-cell incubation and microscopy compatibility. Nat. Commun. 2022, 13, 2999. [Google Scholar]
  4. Naruki, Y.; Kourosh, D.; Mohammad, V.; Animesh, G.; Alán, A. Digital pipette: Open hardware for liquid transfer in self-driving laboratories. Digit. Discov. 2023, 2, 1745–1751. [Google Scholar]
  5. Thurow, K.; Junginger, S. Devices and Systems for Laboratory Automation; Wiley-VCH: Chennai, India, 2022; p. 498. [Google Scholar]
  6. Noh, J.Y.; Yoon, S.W.; Kim, Y.; Lo, T.V.; Ahn, M.J.; Jung, M.C.; Le, T.B.; Na, W.; Song, D.; Le, V.P.; et al. Pipetting-based immunoassay for point-of-care testing: Application for detection of the influenza A virus. Sci. Rep. 2019, 9, 16661. [Google Scholar] [CrossRef]
  7. Haack, C.; Engel, T.; Holenstein, T. Pipetting Device, Modular Pipetting Unit, Pipetting System and Method for Pipetting of Fluid Samples. U.S. Patent US8900878 B2, 2 December 2014. [Google Scholar]
  8. Xin, Y.; Xiaojie, W.; Baoqing, L.; Jiaru, C. A high-precision automated liquid pipetting device with an interchangeable tip. Rev. Sci. Instrum. 2023, 94, 094102. [Google Scholar]
  9. Choi, I.H.; Kim, H.; Lee, S.; Baek, S.; Kim, J. Plug-in nanoliter pneumatic liquid dispenser with nozzle design flexibility. Biomicrofluidics 2023, 9, 064102. [Google Scholar] [CrossRef]
  10. Mao, Y.; Pan, Y.; Li, X.; Li, B.; Chu, J.; Pan, T. High-Precision Digital Droplet Pipetting Enabled by Plug-and-Play Microfluidic Pipetting Chip. Lab Chip 2018, 18, 2720–2729. [Google Scholar] [CrossRef]
  11. Parandakh, A.; Ymbern, O.; Jogia, W.; Renault, J.; Ng, A.; Juncker, D. 3D-printed capillaric ELISA-on-a-chip with aliquoting. Lab. Chip. 2023, 23, 1547–1560. [Google Scholar] [CrossRef] [PubMed]
  12. Berasaluce, A.; Elizalde, J.; Castro, V.; Dominguez, E.; Berganzo, J.; Antoñana, M.; Ruano-López, J.M. Self corrected aliquoting method by overpressure for qPCR multiplexing. In Proceedings of the 19th International Conference on Miniaturized Systems for Chemistry and Life Sciences (µTAS 2015), Gyeongju, Republic of Korea, 25–29 October 2015. [Google Scholar]
  13. Malm, J.; Végvári, A.; Rezeli, M.; Upton, P.; Danmyr, P.; Nilsson, R.; Steinfelder, E.; Marko-Varga, G. Large scale biobanking of blood—The importance of high density sample processing procedures. J. Proteom. 2012, 76, 116–124. [Google Scholar] [CrossRef] [PubMed]
  14. Takase, T.; Masumoto, N.; Shibatani, N.; Matsuoka, Y.; Tanaka, F.; Hirabatake, M.; Kashiwagi, H.; Nishioka, I.; Ikesue, H.; Hashida, T.; et al. Evaluating the safety and efficiency of robotic dispensing systems. J. Pharm. Health Care Sci. 2022, 8, 24. [Google Scholar] [CrossRef] [PubMed]
  15. Hintikka, V. Pipette with a Tip Removal Mechanism, a Method for Removing a Tip, and a Method for Pipetting. European Patent EP3112026B1, 20 February 2008. [Google Scholar]
  16. Johnstone, H.J.H.; Steinke, S.; Kuhnert, H.; Bickert, T.; Pälike, H.; Mohtadi, M. Automated cleaning of foraminifera shells before Mg/Ca analysis using a pipette robot. Geochem. Geophys. Geosystems 2016, 17, 3502–3511. [Google Scholar] [CrossRef]
  17. Zhang, J.; Wan, W.; Tanaka, N.; Fujita, M.; Harada, K. Integrating a Manual Pipette into a Collaborative Robot Manipulator for Flexible Liquid Dispensing. J. Submiss. Reveiew 2022, in press. [Google Scholar]
  18. ISO 13485:2016; Medical Devices. ISO: Geneva, Switzerland, 2016. Available online: https://www.iso.org/standard/59752.html (accessed on 29 February 2024).
  19. ISO 3601-1:2012; Fluid Power Systems. O-rings. Part 1: Inside Diameters, Cross-Sections, Tolerances and Designation Codes. ISO: Geneva, Switzerland, 2012. Available online: https://cdn.standards.iteh.ai/samples/58043/f1be17df93fc436ca4e10e4cc09bf030/ISO-3601-1-2012.pdf (accessed on 29 February 2024).
  20. Żuk, D.; Abramczyk, N.; Panasiuk, K. Analysis of Strength Parameters at Tension of Epoxy-Glass Composites with Rubber Recyclate Addition. J. KONBiN 2022, 52, 131–150. [Google Scholar] [CrossRef]
  21. Steuer, R. Multiple Criteria Optimization: Theory, Computation, and Application; Robert E. Krieger Publishing: Malabar, FL, USA; New York, NY, USA, 1989; pp. 154–196. [Google Scholar]
  22. Simon, D. Evolutionary Optimization Algorithms: Biologically-Inspired and Population-Based Approaches to Computer Intelligence; John Wiley & Sons: Toronto, ON, Canada, 2013; p. 784. [Google Scholar]
  23. Shiakolas, P.S.; Koladiya, D.; Kebrle, J. Optimum robot design based on task specifications using evolutionary techniques kinematic, dynamic, and structural constraints. Machines 2023, 11, 814. [Google Scholar]
  24. Dimo, B.; Michael, E.; Boris, N.; Robin, P. Introduction to Many-Criteria Optimization and Decision Analysis, Many-Criteria Optimization and Decision Analysis; Springer: Cham, Switzerland, 2023; pp. 3–28. [Google Scholar]
  25. Lu, J.; Su, X.; Zhong, J.; Hu, G. Multi-objective shape optimization of developable Bézier-like surfaces using non-dominated sorting genetic algorithm. Mech. Ind. 2023, 24, 38. [Google Scholar] [CrossRef]
  26. Duong, V.; Cuong, N. Expert opinion-based multi objective optimization: An application in plasma coating technology. EUREKA Phys. Eng. 2022, 6, 175–184. [Google Scholar] [CrossRef]
  27. Rybak, L.; Malyshev, D.; Voloshkin, A. Algorithm for Multi-criteria Optimization of Robot Parameters for Fruit Harvesting Based on Evolutionary Methods, Taking into Account the Importance of Criteria. In OPTIMA 2023: Advances in Optimization and Applications; Springer Nature: Cham, Switzerland, 2023; Volume 14, pp. 241–255. [Google Scholar]
  28. Han, K. Dynamics modeling and control of a delta high-speed parallel robot. J. Korean Soc. Manuf. Process Eng. 2014, 13, 90–97. [Google Scholar]
  29. Williams, R.L., II. The Delta Parallel Robot: Kinematics Solutions. Available online: https://people.ohio.edu/williams/html/PDF/DeltaKin.pdf (accessed on 29 December 2023).
  30. Khalapyan, S.; Rybak, L.; Gaponenko, E.; Carbone, G. Motion Control of 6-DOF Relative Manipulation Device. Mech. Des. Robot. 2021, 103, 217–225. [Google Scholar]
  31. Voloshin, A.; Rybak, L.; Saitova, V.; Nozdracheva, A.; Gaponenko, E. Creation of a three-dimensional kinematic model of the Delta manipulator using the computer-aided design system in NX. Robot. Tech. Cybern. 2023, 11, 20–29. [Google Scholar] [CrossRef]
  32. Shuvra, D. An Introduction to Finite Element Analysis Using Matlab Tools; Springer: Cham, Switzerland; New York, NY, USA, 2023; Volume 1, p. 219. [Google Scholar]
  33. Prashant, A.; Amit, S. Finite Element Analysis and Failure Mechanisms of Porous Biomaterial Architecture for Prosthetic Device, Reliability and Risk Assessment in Engineering; Springer: Singapore, 2020; pp. 479–489. [Google Scholar]
  34. Sandeep, C. Improving Strength of 3D Print Using RSM. Int. J. Eng. Appl. Sci. (IJEAS) 2022, 9, 6–8. [Google Scholar]
  35. Bhatti, S.S.; Singh, J. 3D printing of biomaterials for biomedical applications: A review. Int. J. Interact. Des. Manuf. (IJIDeM) 2023, 1–12. [Google Scholar] [CrossRef]
  36. Tang, X.; Tian, Y.; Tian, X.; Li, W.; Han, X.; Kong, T.; Wang, L. Design of multi-scale textured surfaces for unconventional liquid harnessing. Mater. Today 2020, 43, 62–83. [Google Scholar] [CrossRef]
Figure 1. Design of a robotic system for aliquoting biological fluid ((I)—the housing of the system; (II)—sterile work area; (III)—tripod tips; 1—clean zone; 2—tips; 3 –tripod; 4—tripod with small volume tubes; 5—test tube with biological fluid; 6—manipulator; 7—base; 8—dosing device; 9—blood serum; 10—clot; 11—vision system).
Figure 1. Design of a robotic system for aliquoting biological fluid ((I)—the housing of the system; (II)—sterile work area; (III)—tripod tips; 1—clean zone; 2—tips; 3 –tripod; 4—tripod with small volume tubes; 5—test tube with biological fluid; 6—manipulator; 7—base; 8—dosing device; 9—blood serum; 10—clot; 11—vision system).
Machines 12 00172 g001
Figure 2. Diagram of the dosing device. (1—housing assembly; 2—motor coupling; 3—drive gear; 4—motor coupling; 5—external gear wheel; 6—tapper; 7—rod; 8—O-ring; 9—tip; 10—lever; 11—gears; 12—piston rod; 13—piston; 14—chamber; 15—vacuum line).
Figure 2. Diagram of the dosing device. (1—housing assembly; 2—motor coupling; 3—drive gear; 4—motor coupling; 5—external gear wheel; 6—tapper; 7—rod; 8—O-ring; 9—tip; 10—lever; 11—gears; 12—piston rod; 13—piston; 14—chamber; 15—vacuum line).
Machines 12 00172 g002
Figure 3. Calculation scheme of the rubber seal: (a) 3D model of the rod and the rubber seal; (b) calculation scheme in projection (x, y, z); (c) calculation scheme in projection (x, y).
Figure 3. Calculation scheme of the rubber seal: (a) 3D model of the rod and the rubber seal; (b) calculation scheme in projection (x, y, z); (c) calculation scheme in projection (x, y).
Machines 12 00172 g003
Figure 4. The scheme of fixing the tip: (a)—the design scheme of the mechanism; (b)—the scheme of deformation of the ring.
Figure 4. The scheme of fixing the tip: (a)—the design scheme of the mechanism; (b)—the scheme of deformation of the ring.
Machines 12 00172 g004
Figure 5. Calculation diagram: (a) the lever movement along the cam profile; (b) vector triangle.
Figure 5. Calculation diagram: (a) the lever movement along the cam profile; (b) vector triangle.
Machines 12 00172 g005
Figure 6. Design diagram of a pneumatic cylinder.
Figure 6. Design diagram of a pneumatic cylinder.
Machines 12 00172 g006
Figure 7. Pareto Set using an algorithm: (a) NSGA-II, (b) MO-GWO, (c) MO-PSO.
Figure 7. Pareto Set using an algorithm: (a) NSGA-II, (b) MO-GWO, (c) MO-PSO.
Machines 12 00172 g007
Figure 8. The resulting Pareto set: (a) points obtained by all algorithms; (b) points filtered by the Pareto optimality condition.
Figure 8. The resulting Pareto set: (a) points obtained by all algorithms; (b) points filtered by the Pareto optimality condition.
Machines 12 00172 g008
Figure 9. The Pareto set obtained using the MO-PSO at an equivalent time.
Figure 9. The Pareto set obtained using the MO-PSO at an equivalent time.
Machines 12 00172 g009
Figure 10. The mechanism of the dosing device: (a) detailed 3D model of the DD; (b) detailed 3D model of the external gear with a pushing profile (3—drive gear; 5—external gear wheel with pushing profile; 10—lever; 11—gear connected to a screw).
Figure 10. The mechanism of the dosing device: (a) detailed 3D model of the DD; (b) detailed 3D model of the external gear with a pushing profile (3—drive gear; 5—external gear wheel with pushing profile; 10—lever; 11—gear connected to a screw).
Machines 12 00172 g010
Figure 11. A graph of the torque required for the operation of the dosing device.
Figure 11. A graph of the torque required for the operation of the dosing device.
Machines 12 00172 g011
Figure 12. The graph of the rod movement.
Figure 12. The graph of the rod movement.
Machines 12 00172 g012
Figure 13. Rod stress at an oscillation frequency of 0.003 Hz.
Figure 13. Rod stress at an oscillation frequency of 0.003 Hz.
Machines 12 00172 g013
Figure 14. Rod stress during the dynamic movement of the metering device.
Figure 14. Rod stress during the dynamic movement of the metering device.
Machines 12 00172 g014
Figure 15. Manufactured rods: (a) rod manufactured before and after structural calculation with different cross–section and internal diameter; (b) all manufactured rods, including metal.
Figure 15. Manufactured rods: (a) rod manufactured before and after structural calculation with different cross–section and internal diameter; (b) all manufactured rods, including metal.
Machines 12 00172 g015
Figure 16. Design of the dosing device (1—base; 2—cover; 3—upper part of the body; 4—engine; 5—drive gear; 6—external gear wheel with pushing profile; 7—stem; 8—O–ring; 9—washer; 10—lever, 11—gear connected to a screw; 12—stopper; 13—piston; 14—seal; 15—pneumatic cylinder chamber; 16—vacuum line).
Figure 16. Design of the dosing device (1—base; 2—cover; 3—upper part of the body; 4—engine; 5—drive gear; 6—external gear wheel with pushing profile; 7—stem; 8—O–ring; 9—washer; 10—lever, 11—gear connected to a screw; 12—stopper; 13—piston; 14—seal; 15—pneumatic cylinder chamber; 16—vacuum line).
Machines 12 00172 g016
Figure 17. The prototype of the dosing device: (a) fixing the tip; (b) liquid dosing.
Figure 17. The prototype of the dosing device: (a) fixing the tip; (b) liquid dosing.
Machines 12 00172 g017
Figure 18. Examples of liquid intake into the tip: (a) an example of immersion of the tip to a different depth; (b) liquid intake with and without a micro drop.
Figure 18. Examples of liquid intake into the tip: (a) an example of immersion of the tip to a different depth; (b) liquid intake with and without a micro drop.
Machines 12 00172 g018
Table 1. Calculation results and values of calculated parameters.
Table 1. Calculation results and values of calculated parameters.
μα° d 2 , m D t i p , mm h 2 , m f p l p , m g , m/sE, Pa
0.8280.005870.00560.3420.0029.8 5 · 10 6
S, mm2 f m , kg r c , m z 1 z 2 z 3 N , m D z 1 ,   m
3.140.20.00250.032261060.060.0452
Table 2. Comparison of algorithms with an equal number of iterations.
Table 2. Comparison of algorithms with an equal number of iterations.
NSGA-IIMO-GWOMO-PSO
№ ExperimentThe Total Number of Points of the Pareto SetQuantityFraction, %QuantityFraction, %QuantityFraction, %
16358281344.2240837.9113717.9
28499271431.9505759.57288.6
38519274132.2490857.687010.2
48891270630.4534760.18389.4
58863270330.5572764.64334.9
68495272132.0501159.07639.0
78146283234.8499561.33193.9
88758273131.2514858.887910.0
98368273832.7498959.66417.7
108668272931.5527860.96617.6
Average8356.52742.833.14886.857.9726.98.9
Table 3. Comparison of algorithms with equal calculation time.
Table 3. Comparison of algorithms with equal calculation time.
NSGA-IIMO-GWOMO-PSO
№ ExperimentThe Total Number of Points of the Pareto SetQuantityFraction, %QuantityFraction, %QuantityFraction, %
121,816260311.910,23746.9897641.1
234,31825857.5886225.822,87166.6
327,82325959.3975835.115,47055.6
433,14526217.910,41431.420,11060.7
530,55526088.5930430.518,64361.0
622,994258711.310,48445.6992343.2
724,270260410.7978440.311,88249.0
825,883260810.110,10939.113,16650.9
931,63326098.210,91534.518,10957.2
1031,38825988.3915329.219,63762.6
Average28,382.52601.89.4990235.815,878.754.8
Table 4. Values obtained as a result of optimization.
Table 4. Values obtained as a result of optimization.
D P , m S P , mm θ , ° V ,   m m 3 M 2 ,   N · m
0.0120.982962,9390.0406
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rybak, L.; Carbone, G.; Malyshev, D.; Voloshkin, A. Design and Optimization of a Robot Dosing Device for Aliquoting of Biological Samples Based on Genetic Algorithms. Machines 2024, 12, 172. https://0-doi-org.brum.beds.ac.uk/10.3390/machines12030172

AMA Style

Rybak L, Carbone G, Malyshev D, Voloshkin A. Design and Optimization of a Robot Dosing Device for Aliquoting of Biological Samples Based on Genetic Algorithms. Machines. 2024; 12(3):172. https://0-doi-org.brum.beds.ac.uk/10.3390/machines12030172

Chicago/Turabian Style

Rybak, Larisa, Giuseppe Carbone, Dmitry Malyshev, and Artem Voloshkin. 2024. "Design and Optimization of a Robot Dosing Device for Aliquoting of Biological Samples Based on Genetic Algorithms" Machines 12, no. 3: 172. https://0-doi-org.brum.beds.ac.uk/10.3390/machines12030172

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop