We present the results of investigation of flexible stiffness of multilayer beams with geometric and/or stiffness asymmetry. An algorithm is proposed for the evaluation of coordinates of the geometric and stiffness centers and the analysis of flexural stiffness in any direction and its extreme values in multilayer beams with asymmetry of any kind. We study the kinetics of coordinates of the geometric and stiffness centers and flexural stiffness depending on the behavior of the geometric parameters and the ratio of the moduli of elasticity of layers and the cross-sectional shape of a multilayer beam.Introduction. Composite systems with layered structures are used in various fields of engineering because they combine high characteristics of strength, stiffness, and lightness of structural elements. The works [1-6] are devoted to the investigation of multilayer structural elements, mainly with one or two axes of symmetry. The influence of various factors on the regularities of the behavior of strength and stiffness in symmetric multilayer rods and beams is studied in [7,8]. In all cases, it is necessary to determine the location of the neutral layer. The methods used for its determination in the presence of one axis of symmetry are discussed in [1,7,9]. However, in the general case of geometric and stiffness asymmetry, these methods are quite complicated for application. As shown in [7], the location of the neutral layer and, hence, the stiffness of the beam are affected by various factors. Therefore, to study the kinetics of flexural stiffness in multilayer beams with geometric and stiffness asymmetry, it is necessary to determine the coordinates of the geometric and stiffness centers.The contemporary technologies used in manufacturing the products made of layered composite materials enable one to get structural elements of any given shape [10]. As the most extensively used structural elements of rolled steel, one can mention angles, channels, H-beams, and other sections. However, the application of these sections in the case of layered composite materials is restricted by the asymmetry of structural elements and the absence of procedures aimed at their numerical analysis and experimental investigation.The aim of the present work is to develop a mathematical procedure for the analysis of flexural stiffness in the process of elastic deformation in any direction with determination of its extreme values in multilayer beams with asymmetry of any kind and investigation of the kinetics of coordinates of the geometric and stiffness centers and flexural stiffness depending on the behavior of the ratio of moduli of elasticity and the geometric parameters of layers and the cross-sectional shape of multilayer beams.