Introduction. The sediment transportation is one of the major processes that define the magnitude and back surface changing rate for water bodies. The most used prognostic studies in this field are based on the mathematical models that allow reducing, and in some cases, eliminating the expensive and often environmentally burdensome experiments. Spatially one-dimensional models are mainly used to predict changes in bottom topography. For actual coastal systems with irregular coastal configuration, the flow vector is generally not orthogonal to the tangent line for the coastline at each of its points. It also may not coincide with the wind stress vector. Therefore, to solve lots of practically important problems associated with the prediction of the dynamics of the back surface of water basins, it is necessary to use spatially two-dimensional models of sediment transportation and effective numerical methods of their implementation. Materials and Methods. The authors (A.I. Sukhinov, A.E. Chistyakov, E.A. Protsenko, and V.V. Sidoryakina) have earlier proposed a spatially two-dimensional model of sediment transport that satisfies the basic conservation laws (of material balance and momentum) which is a quasilinear equation of parabolic type. The linear difference schemes are constructed and studied; model and some practically important tasks are solved. However, the theoretical study on the proximity of solutions for the original nonlinear initial-boundary value problem and the linearized continuous task - on which basis a discrete model (difference scheme) was built - remained in the shadow. The study of the linearized problem correctness and the determination of sufficient conditions for positivity of solutions are of special interest because only positive solutions to this sediment transport problem have physical value within the framework of the considered models. Research Results. The investigated nonlinear two-dimensional model of sediment transport in the coastal zone of shallow reservoirs takes account of the following physically significant factors and parameters: soil porosity; critical value of the tangent stress at which load transport starts; turbulent interaction; dynamically varying of the bottom geometry; wind currents; and bottom friction. Linearization is carried out on the time grid - nonlinear coefficients of the parabolic equation are taken for the previous time grid step. Next, a chain of tasks connected by the initial data - final solutions of the linearized mixed Cauchy problems chain on a uniform time grid is constructed, and thus, the linearization for the initial 2D nonlinear model is carried out in large. Earlier, the authors have proved the existence and uniqueness of the solution to a linear tasks chain. Prior estimate of the proximity of the linearized problem chain solution to the initial non-linear task solution has been also obtained. The conditions of its solution positivity and their convergence to the nonlinear sediment transport problem are defined in the norm of the Hilbert space L₁ with the rate O(τ) where τ is a time step. Discussion and Conclusions. The obtained research results of the spatially two-dimensional nonlinear sediment transport model can be used for predicting the nonlinear hydrodynamic processes, improving their accuracy and reliability due to the availability of new accounting functionality of physically important factors, including the specification of the boundary conditions.

Introduction. The investigation of elastic wave penetration through periodic gratings is an important problem in the fields of the ultrasonic quantitative evaluation of materials, sound propagation and electromagnetic waveguides with diaphragms. In practice, analytical results can be obtained under the assumption of low frequency, with a weak interaction regime where some approximated results can be established in an analytical form. Materials and Methods . In the previous papers, the 3-D normal penetration of a wave through a plane screen with an infinite doubly periodic system of cracks with a low frequency assumption and the 2-D normal penetration of a wave through a couple of the systems with an infinite periodic array of cracks are studied. The present study objective is to generalize the results obtained before extending the explicit analytical expressions of he considered coefficients for the system of parallel plane screens based on the context of the in-plane problem for the wave propagation through elastic solids with a periodic array of cracks. Research Results . The present work continues to study the 2-D problem for three such parallel arrays which form a doubly-periodic system. The investigation is devoted to the derivation of the analytic expressions for coefficients of reflection and transmission when a longitudinal plane wave falls on the system of three identical two-dimensional gratings. In the one-mode range, the approximation of the problem is reduced to a system of the hypersingular integral equations whose solution gives these coefficients and an explicit representation of the wave field inside the structure. Discussion and Conclusions. The desired control of the acoustic filtering in the considered grating can be arranged by the appropriate choice of the crack’s length, respective frequency interval, and finally - by the distance between two vertical arrays containing periodic systems of cracks.

Introduction. The powder antifriction materials applied in the oil and gas industry are considered. Special attention is given to the materials with the effect of self-lubricity. The use of such materials in a frictional unit of one of the most widespread types of the oil and gas well service equipment - a high-lift plunger pump - is offered. The pump construction diagram with the description of the modernized unit is provided; the problems arising under the production and exploitation of the connecting-rod pin bearing group are described. Materials and Methods. Upon the experience in applying bearing bushes from powder antifriction materials in mechanical engineering, the iron graphite porous material is selected. The design change considering the self-lubricating materials behavior and allowing the best usage of their properties is made. The bushing-to-shaft gap, one of the most important parameters of the plain bearer performance, is considered. The fit providing the fullest seating of the bushing to the shaft is selected. That improves the heat removal, promotes the even load distribution over the frictional area, and so enhances the operating capacity of the bearing unit. Research Results. A change option of the frictional unit design of a connecting-rod small power end with the use of the self-lubricated bearing from the powder antifriction iron-based material is offered. The proposed bearing meets the loading conditions of the considered unit. The design is developed due to the condition of self-lubricity that allows simplifying the production of a connecting rod obviating the manufacturing of oilholes in the connecting rod body. The analysis of the offered modernization is carried out with the application of 3D modeling. The constructed models permit to visually estimate the simplification of the entire connecting rod design, the production of its constituent parts, as well as its assembling. Discussion and Conclusions . The obtained models allow drawing a conclusion on the feasibility of the bearing design change for the sliding friction unit of the connecting-rod end of the high-lift plunger pump using the powder antifriction material to replace the existing one. The application of this bearing with the effect of self-lubricity will allow a substantial simplification of the pin bearing unit design.

Introduction. Any metal-cutting machine has errors of the executive elements movement depending on its geometric accuracy and state. These errors change the forming motions indirectly, through their transformation by the dynamic cutting system with account for the elastic deformations of the tool and the workpiece. Materials and Methods. Laws of the relation between the executive elements and forming paths for the longitudinal cutting case are considered. All deformation displacements are counted off from the independent coordinate system anchored to the carrier system of the machine which is considered to be absolutely rigid and non-deforming. Research Results. A mathematical transformation model consisting of two elastic-deformation subsystems on the part of the tool and workpiece that interact through the dynamic link generated by the cutting process is proposed. This model allows for the errors of the executive elements movement. They are specified by the kinematic disturbances and spindle wavering, and a workpiece in the space. Discussion and Conclusions. Therefore, as opposed to the known studies, it is shown how the mathematical formulation and parameters of the dynamic link generated under the cutting process are changed depending on the kinematic and other types of disturbances. Two types of disturbances are fully considered: those inducing speed variations of the tool motion relative to the workpiece in the direction to the caliper mobility, and disturbances in the form of oscillating displacements in the orthogonal direction to the feed velocity attitude. These disturbances are generated as radial and axial spindle wavering, and as speed variations of the longitudinal feed and variations of the caliper position relative to the ideal rotation axis of the part. It is shown, that the law of the disturbance transformation into the forming movements trajectories depends on the frequency content of the disturbances. The examples are provided, and the conclusions about the transformation properties are made. In particular, it is shown that the speed variations in the direction to the longitudinal feed with the frequency equal to or multiple of the spindle rotation frequency in no way image into the forming movement trajectories.

Introduction. The work objective is the field research of vibrowave processes in the utilization technology of assembly units (AU) of the agricultural machinery. The reducer from SK-4 discarded combine harvester is used for testing. Materials and Methods. A reducer component complex of different sizes and configurations, various materials, varying oil-consumption pollution is considered. After the exterior check of these test pieces, their weighing was performed. The parts information was registered in the related tables. After weighing and checking, the following stage - vibration washing and clearing operation - was carried out. Various components of cleaning the reducer parts surfaces were used for that. They differed both in composition, and in appearance. The following components: PT15х15, porcelain (mm d=6-8.) and steel (d=5-7mm.) full-spheres, were selected and applied. These experiments were performed in the laboratory of the Engineering Technology Department of Don State Technical University. UVG40 vibratory cleaning washers (chamber volume is 40 l) and UVG4x10 (4 chambers per 10 l each) are installed in this laboratory. Both chambers were activated as not only small-sized pieces, but also large-size ones were involved. Cleaning was conducted in 3 steps. The authors consider this quantity optimum for cleaning the rundown reducer parts. Each period was equal to 90 minutes of the chamber operation with the pieces and an abrasive cleaning component. 2% soda ash solution was added to the operation environment. It helped to soften the oil-consumption pollution for its fastest removal from the reducer parts surfaces. After every step, pieces were removed from the cleaning chambers, washed and dried. Further on, each piece was weighed and compared with the basic weight; its degree of the purification of both external and internal part surfaces was studied. Research Results. On completing the cleaning operation, a comparative analysis on the cleaning degree of the surface of each component was carried out. The most effective component was selected. Discussion and Conclusions. Within the framework of the study, the most effective components of cleaning the reducer surfaces are selected.

Introduction. Measuring the linear and angular components of small displacements of the monitor object surfaces is considered. A high-precision non-contact device based on the advanced laser technologies and methods of optical interferometry, namely a double-pass laser interferometer with integrated branches, is developed to perform these measurements. The calculated and experimental validation of the device is offered. Thus, the quality and condition of the materials and products can be controlled by means of the acoustic non-destructive methods. Materials and Methods. New software and mathematical models describing optical fields of the interference patterns generated by the interferometer under consideration are offered. Various types of beam splitters and cases of polarization are used under the numerical simulation of the intensity distributions in the interference patterns. Research Results. New mathematical models and software for the numerical simulation of the optical radiation intensity fields in the interference patterns considering both linear and angular components of small displacement of the surface of the test object are developed. A new interference technique of measuring small linear and angular displacements is developed and theoretically grounded. It allows combining the meter measuring capabilities for both small linear and angular displacements of the monitor object surfaces in a single tool. New optical interference methods and means for the noncontact measurements of small linear and angular movements of the monitor object surfaces that implement the method capabilities and extend the functionality of the known indicators of small displacements are offered. Discussion and Conclusions. The obtained results can be used for the high-precision measurements of small linear and angular displacements of surfaces of the objects under control. Areas of application are the following: experimental research; assessment and diagnostics of the structural materials samples; investigations of fast wave processes in the layered structures of complex shape made of the anisotropic composite materials. The methods developed in this study are applicable in such industries as mechanical engineering, shipbuilding, aircraft engineering, instrument making, power engineering, and etc.

Introduction. Usage of the biocomposites is a vital component of the environment protection policy. According to the classification of the petrochemical polymers, biopolymers can be divided into thermal and plastic heat-strengthening ones. The former can be modified for reheating, and the latter cannot - due to the calcium carbonate solidification. Thermal bioplastics are developed mainly for the packaging industry, and they can be used as a binder in the biocomposite materials. Materials and Methods. Sawdust wood fibers are selected as an additive to produce reinforced rubber plastics (polyester fiber) for this study. The work objective is to determine the optimum ratio of the sawdust volume and size. In the course of the experiments, various sawdust mixtures are tested through varying the combination of their size, volume and applied loads. Research Results. The obtained results show a growth of the composite material deformation when increasing the proportion of sawdust of any volumes studied under the pressure testing, as well as a decrease in strength and tensile strain when increasing the proportion of sawdust of medium and low volumes. The strength and elongation of the composite grows with increasing the sawdust fraction volume for more than 120 mm³. Stresses increase when the size of wood sawdust is 20 mm³ and mixing ratio is 25%, therewith the material elasticity decreases. Discussion and Conclusions. In bending tests, the reduction in strength and increase in tension strain when mixing a sample with sawdust of 5 mm³, gain in strength and elongation at increasing their volume up to 25 mm³, and then the reduction of strength and tension strain at increasing their volume to 120 mm³ are observed.

Introduction . One of the widely accepted methods of the precision metal turning provided a fairly complex shape of the workpiece is the method of multiple repetitions of turning. The principle of this method is cutting an extra thin layer of the work material. At that, the required operating depth is obtained through the successive passes of the cutting tool along the workpiece. An essential fault of this method of complex parts processing is low efficiency due to the long processing time. Materials and Methods . The improvement of efficiency is possible through the introduction of new modern methods of the treatment process control. Modern approaches to the analysis and synthesis of the control systems of complex operating procedures based on the development of the synergetic concept of management could be used as such methods. A special feature of this approach is the identification and control of the inherent feedbacks stabilizing the treatment process. On the basis of the analysis conducted, the control that is a new nonlinear degenerative feedback is synthesized. Its task is to assure adequate quality and efficiency of the processing under the asymptotic stability of the required attractors. Research Results . A fundamentally new mathematical model considering non-linear processes of the force changes is developed to simulate the system and allowance for the effect of the engine parameters. The analysis of the dynamics of the metal turning control system with the evaluation of the process quality is carried out. A synergistic approach to the construction of the optimal production control system considering the interconnectivity of the actuators providing the workpiece rotation and supply is used. The analysis of the synthesized system dynamics with the comparative evaluation of the quality growth is conducted. Discussion and Conclusions . The synthesized synergetic control law changes fundamentally the properties of the metal turning control system. A comparative analysis of the obtained dependences allows observing a substantial increase in the processing accuracy. At that, a synthesized system control asymptotically approaches the deviation of the allowance to the value specified as one of the constants on the control channel of the feed velocity.

Introduction. There is a performance improvement problem of the existing wind-driven power plants (WPP). It is worthwhile investigating nonclassical types of the WPP to solve it. The most promising are the WPP which use special guide static elements to create vortex structures. The latter contributes to the increased aerodynamic efficiency of the WPP. Materials and Methods . Ansys Workbench and Ansys Fluent software packages are used to calculate the design aerodynamics of the wind turbines. Research Results. On the analysis of the particular publications, the viability of the vortex-type WPP studies in which the additional stator structures are used to increase their energy efficiency, is proved. The aerodynamics analysis of the WPP with “rotor-in-socket” assembly type is carried out using the continuum dynamics computational methods. An advanced WPP design is obtained by the aerodynamic optimization technique. Discussion and Conclusions. The obtained results show that the construction of the vortex-type turbine on the basis of the favourable interference concept of rotor and stator is promising. The specific form of WPP with a vertical axis of rotation, axisymmetric stator and socket of special forms is analyzed. The methods of the aerodynamic optimization of the described WPP design are considered. On the basis of this methodology, an advanced socket form and its assembly with the rotor which allows reaching an average calculated value of the favourable air moment that exceeds more than twofold the corresponding value for the initial geometry of the WPP, are obtained.

Introduction. Currently, wheel tread surfacing is implemented on special design wheel lathes. This method provides the desired technological strength of the wheel treads, but does not provide hard facing and, consequently, does not increase wear resistance. Ball-rod hardening (BRH) has great advantages as a process, especially for composite geometric surfaces; and besides, it enables the necessary compressive stresses in the surface of the hardenable products. However, it has a significant disadvantage in terms of the increased levels of the acoustic pressure in the operators’ working area. Materials and Methods . The experimental studies under BRH of the wheelpair assembly elements are carried out by several process steps. Wheel tread surfacing is performed on a special design wheel lathe on which ball-rod reinforcers are installed instead of a cutter. Hardening of the mounting surfaces of the wheelset axles is performed on a special axial lathe on which, in a similar manner, ball-rod reinforcers are installed instead of a cutter. Hardening of the wheel mounting surface is performed on a boring-and-turning mill. The octave sound pressure levels are measured on all of these lathes at idle and under the hardening process implementation. It allows determining the contribution of the sound radiation of the acoustic subsystem “hardenable product - ball-rod reinforcer” into the soundfield generated by a common lathe loudspeaker system in the operators’ working area. Research Results . Due to the high intensity of the radiant sound energy, the equipment is installed in separate facilities commensurable with the overall machine dimensions that consequently contributes to the increased noise activity. Therefore, the research data are of current interest for the enterprises of mechanical engineering, and they have both scholary interest and practical applicability. Discussion and Conclusions. It should be noted that the data obtained on the regular engine lathes, i.e. on round work pieces, are a particular case and are not applicable to the assembly of such complex geometrics as wheelpairs of rolling stock and various types of cranes. The implemented measurements show that under BRH, sound pressure levels exceed health standards by 5 to 15 dB of the high-frequency part of 500-8000 Hz spectrum, i.e. in the range in which the maximum permissible levels should not exceed 78-69 dB.

Introduction. Research results of the part surface treatment by a ball-and-rod reinforcer (multicontact shock-vibrating tool) through the surface plastic deformation method are provided. The considered process can be used for hardening both flat and curved surfaces, as well as for generating compressing residual voltages, pocket smoothing, regular microrelief fitting on friction pairs, and for joint weld machining. When carrying out this research, the manufacturing processability is identified, the regularities of selecting design data of the device are determined. The treatment process is percussive, and it is attended by heavy acoustic radiation. It is found that noise under processing far exceeds the allowable range. Materials and Methods. The investigation was conducted on a universal milling machine. Long flat parts were exposed to hardening. The parts were installed directly on the machine table or in special devices. A line noise source was accepted as a model. Research Results. Dependences are obtained for determining the sound pressure, natural frequencies of the part oscillations, and noise spectra. A set of noise abatement procedures under processing is worked out. A considerable amount of the experimental research confirming adequacy of the obtained theoretical models is undertaken. Discussion and Conclusions. When carrying out the studies, it is established that efficiency in noise abatement cannot be obtained by the rational selection of the vibration-absorbing material. In this regard, the use of passive noise protection of the machine site in the form of a flat polycarbonate-glass screen is recommended. This will allow reducing noise in the work zone to the standard values. Design data of the offered acoustic screen are calculated.

Introduction. Discrepancy of the mathematical model structure to the observed process is one of the causes for loss in the control system operation. It determines the urgency of the problem of structural synthesis of the adaptive tracking systems which is considered in the extremum formulation. Materials and Methods. The combined maximum principle methodology is used to solve the synthesis problem in the quasideterministic approximation. In this case, the minimum condition of the objective functional which is the result of the convolution of the quality criterion and action integral is investigated. This allows building a state filter model with the structural adaptation as a solution to the inverse dynamic problem. Practical implementation of the results received through this approach in the digital tracking systems sets a problem of selecting a finite-dimensional approximation scheme for the obtained model which is associated with the determination of the synthesizing function approximation form and the use of the general solution to the homogeneous differential equation. Research Results. Embodiment of the digital tracking systems based on the finite-dimensional approximation with a new structure of the transfer matrix and covariance matrix of the disturbance intensity vector in discrete time is considered. Discussion and Conclusions. The mathematical modeling results show that a new method of the tracking systems synthesis in discrete time under regular perturbances improves accuracy of the estimators operation compared with the traditional approaches.

Introduction. The authors’ method of the objective evaluation of the maturity level of the organization for obtaining sustainable success on the basis of self-assessment is described. The technique is based on the use of expert information and mathematical apparatus of the fuzzy sets theory. The substantiation of the chosen approach is presented. It allows promptly analyze the organization state and apply full expert information. It is proposed to use expert assessments of maturity levels for each of the activities according to GOST R ISO 9004-2010. A brief description of the fuzzification stages content, composition and defuzzification which are the essence of the fuzzy modeling process to obtain reliable estimates is given. Materials and Methods . As follows from the analysis, a generalized model of the “evaluation of the maturity level of the organization” domain is developed. The input and output linguistic variables corresponding to the maturity main criteria in accordance with GOST R ISO 9004-2010 are simulated; and the base term sets are defined. Membership functions of all the linguistic variables are developed according to the estimates of four experts for which purpose the typical trapezoidal functions are used. An expanded form of the fuzzy logic inference for the considered system of knowledge is presented. Taking into account the calculated significance coefficients of the basic maturity criteria, a knowledge base is created. It includes 729 production rules and serves as a basis for the logical deduction of the evaluation of the organization maturity level. Research Results. The barycenter method is used to calculate the precise value of the maturity level. It is implemented in Matlab with the help of Fuzzy Logic Toolbox application package. Response surfaces showing the relationship of the linguistic variables are constructed. The fuzzy inference of the evaluation of the organization maturity level using the proposed techniques is exemplified. Discussion and Conclusions. The presented methodology using a mathematical apparatus of the fuzzy sets theory allows an objective estimation of the organization maturity level for sustained success. It is designed for prompt analysis of the organization on the basis of the self-assessment performed by experts.

Introduction. Boolean functions of several variables are of paramount importance in the coding theory and cryptography. The compositions of these functions are used in a set of the symmetric cryptosystems; therewith, some error-control codes, such as Reed-Muller codes, Kerdock codes, can be defined; as well as some new decoders operating beyond half of the code distance can be constructed. The task of restoring a Boolean function from its derivative which is called a Boolean function integration problem is considered. A Boolean function being restored, the vector towards which the derivative is calculated is supposed unknown. Materials and Methods . The results are obtained on the basis of the following methodology: theory of Boolean functions, theory of finite fields and polynomial rings, linear algebra. The space of Boolean functions is considered a certain isomorphic factor-ring that allows reducing the task to finding solutions to a polynomial set of equations of a special form. The constructed isomorphism enables to check whether the integration problem is decidable, and also to offer a new method of its solution. Research Results . The algorithm of searching preimage by the full enumeration method is formally constructed; and its algorithmic complexity is calculated. The theorem of necessary and sufficient conditions for the existence of an arbitrary Boolean function preimage regarded as the directional derivative value is proved. The provided proofs are constructive. On the basis of the established facts, the algorithms of checking the preimage existence for the specified Boolean function and of building the preimage are developed. In the proposed version, the algorithm forms only one of the possible preimages under the condition of its existence. The proposed algorithm of the preimage generation is significantly efficient from the standpoint of the algorithmic complexity compared to the full enumeration method. Time estimates of the complexity of the basic formal algorithms developed for solving the formulated problems are given. The comparison of their operation complexity to the algorithm of Boolean functions integration complexity by the complete enumeration method is described. Discussion and Conclusions. The research performed can be useful for special sections of the coding theory and cryptography where Boolean functions of several variables are used.

Introduction. Schedules of classes and exams are crucial components in the training process organization. Their quality directly affects the successful mastering of learning materials and ensures physical and psychological health of the educational process participants. The influence of educational systems features on the problem setting and solving methods under timetabling is considered. The emphasis is upon the timetabling problems in the higher education systems. Materials and Methods. The analytical study results of a number of the Russian and foreign sources on methods and algorithms of the classes and exams timetabling automation are described. These techniques and algorithms are compared from the standpoint of modern system-analysis methods. Schedule is considered as part of the educational process support system, and, in its turn, has features of the system combining such objects as students, teachers, disciplines, and classrooms. The application of the system analysis methods allows allocating essential features of the implemented systems of timetabling, classifying and evaluating them. Research Results. It is found that the tasking of classes scheduling and exams timetabling are practically equivalent and include the timetable quality criterion, soft and hard constraints. Russian sources contain detailed studies of various tasks of the classes timetabling. At the same time, no papers devoted to developing methods and algorithms for the examination schedules are found among these sources. However, foreign sources are largely devoted to the examination schedules. Graph theory methods, heuristic and hybrid optimization algorithms, linear integer programming techniques, system analysis methods, and other modern techniques of discrete mathematics are applied to solve the problems in hand. Discussion and conclusions. In general, in one and the same institution, the dimension, as well as the solution to the examination scheduling problem is less complicated than the classes scheduling. In connection with the expansion of the individualization of the learning process in the Russian educational system, the improvement of the research methods and algorithms relevance for constructing exams schedule is expected.

Introduction. The research area of “natural calculation” is now widely used for the solution to optimization NP-complete problems including combinatorial tasks of cryptanalysis. A quick overview of the publications devoted to the application of the natural (bioinspired) methods for cryptanalysis is provided. The main work objective is to investigate a possibility of applying bee colony algorithms to the realization of block cipher cryptanalysis. Materials and Methods . The known bee colony techniques belonging to a relatively new class of the bioinspired optimization methods that simulate the processes occurring in wildlife are applied to solve this optimization problem. The description and the block diagram of the bee colony algorithm for the solution to a cryptanalysis task are provided; basic operations performed in parallel at the global level are noted. In the following, a set of independent operators allowing for the concurrent execution is defined. For this purpose, information-logical flowgraphs of the algorithm with the input control and information links are built, as well as matrices of succession, logical incompatibility, and independence are formed. This matrix of independence allows the definition of a set of algorithm operators admitting parallel execution. At that, the dimensionality of the maximal internally stable sets defines the maximum number of the processors used for the algorithm implementation. Research Results . Theoretical estimates of time complexity of the bee colony algorithm are given as the key data. Besides, the problem solution is provided: to find the required smallest number of processors of the homogeneous parallel computing systems with distributed memory, and a uniform plan for the implementation of operators for them, for the cryptanalysis algorithm based on the constructed information-logical graph data-logical graph, and for the preset time. The assessment of the wanted smallest number of processors for the cryptanalysis algorithm implementation, and the evaluation of the total time of the algorithm realization are given. Discussion and Conclusions. The basic research results are: the development of the bee colony algorithm used for the cryptanalysis task solution; the description of its flowchart and the principal par-allel executed stages; the construction of a matrix of independence; the evaluation of the number of processors for the algorithm imple-mentation. It should be noted (and it was observed in the previous works) that the distinctive feature of applying the bioinspired meth-ods of cryptanalysis is the applicability of the encryption-decryption algorithm as a criterion function for the evaluation of the key ac-ceptability defined by the bioinspired method operations. Thus, it can be affirmed that when using the bioinspired methods, the secret key definition process depends not so much on the complexity of the encryption transformations, as on the bioinspired method itself which should provide a sufficient diversity of the key generation