A generalised plate with kinematically independent thickness for modelling shapes of corneas affected by keratoconus before and after penetrating keratoplasty
Restricted accessResearch articleFirst published online July, 2026
A generalised plate with kinematically independent thickness for modelling shapes of corneas affected by keratoconus before and after penetrating keratoplasty
A two-dimensional (2D) reduced-order generalised continuum model within the framework of the three-dimensional (3D) deformations is deduced from a 3D Cauchy continuum model by imposing a micro-macro kinematical map, which is linear in the direction normal to the corneal surface. This kinematical assumption is plausible as the cornea thickness is much smaller than its diameter. We use the obtained 2D generalised continuum that incorporates a kinematically independent thickness to model the changes of shape induced in corneas: (1) by the changes of cornea mechanical properties whose aetiology can be found in the complex (and not completely understood yet) pathogenic process causing keratoconus, (2) by penetrating keratoplasty, and (3) degeneration of both patient residual corneal tissue and transplanted corneal tissue after transplant. We postulate that growth and regeneration phenomena occurring in the cornea shape it following the “elastic” solutions, which we have calculated. The preliminary obtained predictions seem to promise significant applicative developments and are in good qualitative agreement with experimental results: future investigations will need to improve the presented model by considering explicitly the remodelling phenomena and a more detailed analysis of the evolution of metabolically driven mechanical damage of corneal tissue and its visco-plasticity.
NottinghamJ. Practical observations on conical cornea, and on the short sight, and other defects of vision connected with it. London: J. Churchill, 1854.
2.
PandolfiAHolzapfelGA. Three-dimensional modeling and computational analysis of the human cornea considering distributed collagen fibril orientations. J Biomech Eng2008; 130(6): 061006.
3.
PandolfiAManganielloF. A model for the human cornea: constitutive formulation and numerical analysis. Biomech Model Mechanobiol2006; 5(4): 237–246.
dell’IsolaFVastaMGizziA, et al. 2D generalized continuum modeling cornea changes of shapes induced by pathogenetical process related mechanical damage. In: Mathematical modelling in biology and medicine, Arpino, Italy, 8–12 May, 2023.
6.
EpsteinA. Keratoconus and related disorders. North Shore Contact Lens2000; 87: 87.
7.
Romero-JiménezMSantodomingo-RubidoJWolffsohnJS. Keratoconus: a review. Cont Lens Anterior Eye2010; 33(4): 157–166.
8.
HashemiHHeydarianSHooshmandE, et al. The prevalence and risk factors for keratoconus: a systematic review and meta-analysis. Cornea2020; 39(2): 263–270.
9.
GrilloADi StefanoS. Comparison between different viewpoints on bulk growth mechanics. Math Mech Complex Syst2023; 11(2): 287–311.
10.
GrilloADi StefanoS. An a posteriori approach to the mechanics of volumetric growth. Math Mech Complex Syst2023; 11(1): 57–86.
11.
TepedinoM. The mechanical role of the periodontal ligament for developing mathematical models in orthodontics. Math Mech Complex Syst2023; 11(4): 525–539.
12.
CoganDG. Clinical physiology of the cornea. Bull Ophthal Soc Egypt1953; 46: 484–498.
13.
AltenbachJAltenbachHEremeyevVA. On generalized Cosserat-type theories of plates and shells: a short review and bibliography. Arch Appl Mech2010; 80(1): 73–92.
14.
AltenbachHEremeyevV, et al. Actual developments in the nonlinear shell theory—state of the art and new applications of the six-parameter shell theory. Shell Struct Theory Appl2014; 3: 3–12.
15.
EremeyevVPietraszkiewiczW, et al. Refined theories of plates and shells. Z Angew Math Mech2014; 94(1–2): 5.
16.
EremeyevVALebedevLPCloudMJ. The Rayleigh and Courant variational principles in the six-parameter shell theory. Math Mech Solids2015; 20(7): 806–822.
17.
PietraszkiewiczWEremeyevVKonopińskaV. Extended non-linear relations of elastic shells undergoing phase transitions. Z Angew Math Mech2007; 87(2): 150–159.
18.
AuffrayNdell’IsolaFEremeyevVA, et al. Analytical continuum mechanics à la Hamilton–Piola least action principle for second gradient continua and capillary fluids. Math Mech Solids2015; 20(4): 375–417.
19.
dell’IsolaFSeppecherPMadeoA. How contact interactions may depend on the shape of Cauchy cuts in Nth gradient continua: approach “à la D’Alembert”. Z Angew Math Phys2012; 63: 1119–1141.
LubardaVChenM. On the elastic moduli and compliances of transversely isotropic and orthotropic materials. J Mech Mater Struct2008; 3(1): 153–171.
22.
De BellisMLVastaMGizziA, et al. A numerical model of the human cornea accounting for the fiber-distributed collagen microstructure. Math Mech Solids. Epub ahead of print 8 November 2023. DOI: 10.1177/10812865231202024.
23.
De BellisMLVastaMGizziA, et al. Modeling the deterioration of the stiffness and of the collagen fibril distribution in a discrete model of the cornea microstructure. Int J Non Linear Mech2024; 163: 104736.
24.
PandolfiAGizziAVastaM. A microstructural model of cross-link interaction between collagen fibrils in the human cornea. Philos Trans R Soc A2019; 377(2144): 20180079.
25.
PandolfiADe BellisMLGizziA, et al. Modeling the degeneration of the collagen architecture in a microstructural model of the human cornea. Math Mech Solids2023; 28(1): 196–207.
26.
PyeDC. A clinical method for estimating the modulus of elasticity of the human cornea in vivo. PLoS ONE2020; 15(1): e0224824.
27.
ShihPJHuangCJHuangTH, et al. Estimation of the corneal Young’s modulus in vivo based on a fluid-filled spherical-shell model with Scheimpflug imaging. Journal of ophthalmology2017; 2017(1): 5410143.
28.
dell’IsolaFMaierGPeregoU, et al. (eds) The complete works of Gabrio Piola: Volume I. commented English translation. advanced structured materials. Cham: Springer, 2014.
29.
Dell’IsolaFKosinskiW. Deduction of thermodynamic balance laws for bidimensional nonmaterial directed continua modelling interphase layers. Arch Mech1993; 45(3): 333–359.
30.
BoutinCVivergeK. Generalized plate model for highly contrasted laminates. Eur J Mech A/Solids2016; 55: 149–166.
31.
JakabčinLSeppecherP. On periodic homogenization of highly contrasted elastic structures. J Mech Phys Solids2020; 144: 104104.
32.
VivergeKBoutinCSalletF. Model of highly contrasted plates versus experiments on laminated glass. Int J Solids Struct2016; 102: 238–258.
33.
KöryJStewartPSHillNA, et al. A discrete-to-continuum model for the human cornea with application to keratoconus. R Soc Open Sci2024; 11(7): 240265.
34.
MurgatroydHBembridgeJ. Intraocular pressure. Conti Edu Anaesth Critical Care Pain2008; 8(3): 100–103.
35.
RamierAEltonyAMChenY, et al. In vivo measurement of shear modulus of the human cornea using optical coherence elastography. Sci Report2020; 10(1): 17366.
36.
GrecoLCuomoMContrafattoL. Two new triangular G1-conforming finite elements with cubic edge rotation for the analysis of Kirchhoff plates. Comput Method Appl M2019; 356: 354–386.
37.
CuomoMGrecoL. An implicit strong G1-conforming formulation for the analysis of the Kirchhoff plate model. Continuum Mech Therm2020; 32(3): 621–645.
38.
GrecoLCuomoM. An implicit G1-conforming bi-cubic interpolation for the analysis of smooth and folded Kirchhoff–Love shell assemblies. Comput Method Appl M2021; 373: 113476.
39.
BattistaARosaLDell’ErbaR, et al. Numerical investigation of a particle system compared with first and second gradient continua: deformation and fracture phenomena. Math Mech Solids2017; 22(11): 2120–2134.
40.
Dell’ErbaRD’avanzoPRapisardaAC. A comparison between the finite element method and a kinematic model derived from robot swarms for first and second gradient continua. Continuum Mech Therm2023; 35(4): 1769–1786.
41.
LaudatoMManzariLGöranssonP, et al. Experimental analysis on metamaterials boundary layers by means of a pantographic structure under large deformations. Mech Res Commun2022; 125: 103990.
42.
BarchiesiECiallellaAGiorgioI. On boundary layers observed in some 1D second-gradient theories. In: BarchiesiECiallellaAGiorgioI (eds) Theoretical analyses, computations, and experiments of multiscale materials: a tribute to francesco dell’isola. Cham: Springer, 2022, pp. 359–376.
43.
EremeyevVALurieSASolyaevYO, et al. On the well posedness of static boundary value problem within the linear dilatational strain gradient elasticity. Z Angew Math Phys2020; 71: 1–16.
44.
TurcoEBarchiesiE. Kinematically triggered nonlinear vibrations of Hencky-type pantographic sheets. Math Mech Complex Syst2022; 9(3): 311–335.
45.
CiallellaA. Research perspective on multiphysics and multiscale materials: a paradigmatic case. Continuum Mech Therm2020; 32: 527–539.
46.
EugsterSR. Numerical analysis of nonlinear wave propagation in a pantographic sheet. Math Mech Complex Syst2022; 9(3): 293–310.
47.
YildizdagMESararBCSalvatoriA, et al. Analysis of transmission and reflection characteristics of linear plane waves in pantographic lattices. Z Angew Math Phys2023; 74(5): 178.
48.
ElsheikhAWangDBrownM, et al. Assessment of corneal biomechanical properties and their variation with age. Current eye research2007; 32(1): 11–19.
49.
MeekKMTuftSJHuangY, et al. Changes in collagen orientation and distribution in keratoconus corneas. Invest Ophthalmol Vis Sci2005; 46(6): 1948–1956.
50.
Romero-JiménezMSantodomingo-RubidoJWolffsohnJS. Keratoconus: a review. Cont Lens Anterior Eye2010; 33(4): 157–166.
51.
ShawAJCollinsMJDavisBA, et al. Eyelid pressure and contact with the ocular surface. Invest Ophthalmol Vis Sci2010; 51(4): 1911–1917.
52.
DoaneMG. Interactions of eyelids and tears in corneal wetting and the dynamics of the normal human eyeblink. American journal of ophthalmology1980; 89(4): 507–516.
53.
MathewsPBenbowACorcoranK, et al. 2022eye banking statistical report—executive summary. Eye Bank Corneal Transplant2023; 2(3): e1–e12.
54.
GurnaniBKaurK. Penetrating keratoplasty. Treasure Island, FL: StatPearls Publishing, 2023.
55.
AbaliBEMüllerWHEremeyevVA. Strain gradient elasticity with geometric nonlinearities and its computational evaluation. Mech Adv Mater Mod Process2015; 1: 1–11.
56.
EremeyevVAAltenbachH. On the direct approach in the theory of second gradient plates. In: AltenbachHMikhasevGI (eds) Shell and membrane theories in mechanics and biology: from macro-to nanoscale structures. Cham: Springer, 2014, pp. 147–154.
57.
AbaliBE. Energy based methods applied in mechanics by using the extended Noether’s formalism. Z Angew Math Mech2023; 103: e202300020.
58.
AbaliBEVazicBNewellP. Influence of microstructure on size effect for metamaterials applied in composite structures. Mech Res Commun2022; 122: 103877.
EremeyevVA. Strong ellipticity conditions and infinitesimal stability within nonlinear strain gradient elasticity. Mech Res Commun2021; 117: 103782.
61.
EremeyevVALazarM. Strong ellipticity within the Toupin–Mindlin first strain gradient elasticity theory. Mech Res Commun2022; 124: 103944.
62.
Dell’IsolaFSciarraGVidoliS. Generalized Hooke’s law for isotropic second gradient materials. Proc R Soc A Math Phys Eng Sci2009; 465(2107): 2177–2196.
63.
Abdoul-AnzizHSeppecherP. Strain gradient and generalized continua obtained by homogenizing frame lattices. Math Mech Complex Syst2018; 6(3): 213–250.
64.
FedeleR. Piola’s approach to the equilibrium problem for bodies with second gradient energies. Part I: first gradient theory and differential geometry. Continuum Mech Therm2022; 34(2): 445–474.
65.
FedeleR. Approach à la Piola for the equilibrium problem of bodies with second gradient energies. Part II: variational derivation of second gradient equations and their transport. Continuum Mech Therm2022; 34(5): 1087–1111.
66.
SugarJMacsaiMS. What causes keratoconus?Cornea2012; 31(6): 716–719.
67.
ZhangHScarcelliGRandlemanJB. Comparison of superior and inferior corneal biomechanics by polarization-sensitive Brillouin microscopy. Invest Ophth Vis Sci2024; 65(7): 2060–2060.
68.
DebourdeauEPlanellsGChamardC, et al. New keratoconus risk factors: a cross-sectional case—control study. Journal of Ophthalmology2022; 2022(1): 6605771.
69.
KrachmerJHFederRSBelinMW. Keratoconus and related noninflammatory corneal thinning disorders. Surv Ophthalmol1984; 28(4): 293–322.
70.
GiorgioI. A variational formulation for one-dimensional linear thermoviscoelasticity. Math Mech Complex Syst2022; 9(4): 397–412.
71.
CuomoMContrafattoLGrecoL. A cohesive interface model with degrading friction coefficient. Math Mech Complex Syst2024; 12(2): 113–133.
72.
AltenbachHEremeyevV. Thin-walled structural elements: classification, classical and advanced theories, new applications. Shelllike Struct Adv TheorAppl. Epub ahead of print 10 August 2016. DOI: 10.1007/978-3-319-42277-0_1.
73.
AltenbachHEremeyevV. On the effective stiffness of plates made of hyperelastic materials with initial stresses. Int J Non Linear Mech2010; 45(10): 976–981.
74.
AltenbachHEremeyevVANaumenkoK. On the use of the first order shear deformation plate theory for the analysis of three-layer plates with thin soft core layer. Z Angew Math Mech2015; 95(10): 1004–1011.
75.
FedeleR. Third-gradient continua: nonstandard equilibrium equations and selection of work conjugate variables. Math Mech Solids2022; 27(10): 2046–2072.
76.
FedeleR. Deformation-induced coupling of the generalized external actions in third-gradient materials. Z Angew Math Phys2022; 73(5): 218.
77.
De AngeloMBarchiesiEGiorgioI, et al. Numerical identification of constitutive parameters in reduced-order bi-dimensional models for pantographic structures: application to out-of-plane buckling. Arch Appl Mech2019; 89(7): 1333–1358.
78.
VazicBAbaliBEYangH, et al. Mechanical analysis of heterogeneous materials with higher-order parameters. Eng with Comput2022; 38(6): 5051–5067.
79.
AbaliBEKlunkerABarchiesiE, et al. A novel phase-field approach to brittle damage mechanics of gradient metamaterials combining action formalism and history variable. Z Angew Math Mech2021; 101(9): e202000289.
80.
PlacidiLBarchiesiEMisraA. A strain gradient variational approach to damage: a comparison with damage gradient models and numerical results. Math Mech Complex Syst2018; 6(2): 77–100.
81.
PlacidiLBarchiesiEMisraA, et al. Micromechanics-based elasto-plastic–damage energy formulation for strain gradient solids with granular microstructure. Continuum Mech Therm2021; 33(5): 2213–2241.
82.
PlacidiLMisraABarchiesiE. Two-dimensional strain gradient damage modeling: a variational approach. Z Angew Math Phys2018; 69: 1–19.
83.
CuomoM. Continuum damage model for strain gradient materials with applications to 1D examples. Continuum Mech Therm2019; 31(4): 969–987.
84.
NaumenkoKEremeyevVA. A non-linear direct peridynamics plate theory. Compos Struct2022; 279: 114728.
85.
PlacidiLGrecoLBucciS, et al. A second gradient formulation for a 2D fabric sheet with inextensible fibres. Z Angew Math Phys2016; 67: 1–24.
86.
PideriCSeppecherP. A second gradient material resulting from the homogenization of an heterogeneous linear elastic medium. Continuum Mech Therm1997; 9: 241–257.
87.
EremeyevVAdell’IsolaFBoutinC, et al. Linear pantographic sheets: existence and uniqueness of weak solutions. J Elasticity2018; 132: 175–196.
88.
TrinhDKJanickeRAuffrayN, et al. Evaluation of generalized continuum substitution models for heterogeneous materials. Int J Multiscale Com2012; 10(6).
89.
EremeyevVAFreidinABPavlyuchenkoVN, et al. Instability of hollow polymeric microspheres upon swelling. Dokl Phys2007; 52(1): 37–40.
90.
SoizeC. An overview on uncertainty quantification and probabilistic learning on manifolds in multiscale mechanics of materials. Math Mech Complex Syst2023; 11(1): 87–174.
91.
La ValleGCiallellaAFalsoneG. The effect of local random defects on the response of pantographic sheets. Math Mech Solids2022; 27(10): 2147–2169.
92.
MühlichUAbaliBEdell’IsolaF. Commented translation of Erwin Schrödinger’s paper ‘On the dynamics of elastically coupled point systems’(Zur Dynamik elastisch gekoppelter Punktsysteme). Math Mech Solids2020; 26(1): 1081286520942955.
93.
MandadapuKKAbaliBEPapadopoulosP. On the polar nature and invariance properties of a thermomechanical theory for continuum-on-continuum homogenization. Math Mech Solids2021; 26(11): 1581–1598.
EremeyevVAPietraszkiewiczW. Phase transitions in thermoelastic and thermoviscoelastic shells. Arch Mech2009; 61(1): 41–67.
96.
EremeyevVAKonopińska-ZmysłowskaV. On the correspondence between two-and three-dimensional Eshelby tensors. Continuum Mech Therm2019; 31: 1615–1625.
97.
AbaliBEBarchiesiE. Additive manufacturing introduced substructure and computational determination of metamaterials parameters by means of the asymptotic homogenization. Continuum Mech Therm2021; 33: 993–1009.
98.
VazicBAbaliBENewellP. Generalized thermo-mechanical framework for heterogeneous materials through asymptotic homogenization. Continuum Mech Therm2023; 35(1): 159–181.
99.
YangHAbaliBEMüllerWH, et al. Verification of asymptotic homogenization method developed for periodic architected materials in strain gradient continuum. Int J Solids Struct2022; 238: 111386.
100.
AydinGSararBCYildizdagME, et al. Investigating infill density and pattern effects in additive manufacturing by characterizing metamaterials along the strain-gradient theory. Math Mech Solids2022; 27: 2002–2016.
101.
SararBCYildizdagMEAbaliBE. Comparison of homogenization techniques in strain gradient elasticity for determining material parameters. In: AltenbachHBerezovskiAdell’IsolaF, et al. (eds) Sixty shades of generalized continua: dedicated to the 60th birthday of Prof. Victor A. Eremeyev. Cham: Springer, 2023, pp. 631–644.
102.
EremeyevVA. Two-and three-dimensional elastic networks with rigid junctions: modeling within the theory of micropolar shells and solids. Acta Mech2019; 230(11): 3875–3887.
EremeyevVAPietraszkiewiczW. Nonlinear resultant theory of shells accounting for thermodiffusion. Continuum Mech Therm2021; 33: 893–909.
105.
EremeyevVAPietraszkiewiczW. Local symmetry group in the general theory of elastic shells. J Elasticity2006; 85: 125–152.
106.
EremeyevVZubovL. On constitutive inequalities in nonlinear theory of elastic shells. Z Angew Math Mech2007; 87(2): 94–101.
107.
MaHGaoXLReddyJ. A non-classical Mindlin plate model based on a modified couple stress theory. Acta Mech2011; 220(1): 217–235.
108.
AbaliBEMüllerWHdell’IsolaF. Theory and computation of higher gradient elasticity theories based on action principles. Arch Appl Mech2017; 87(9): 1495–1510.
109.
ShekarchizadehNAbaliBEBersaniAM. A benchmark strain gradient elasticity solution in two-dimensions for verifying computational approaches by means of the finite element method. Math Mech Solids2022; 27(10): 2218–2238.
110.
PlacidiLEl SherbinyMGBaragattiP. Experimental investigation for the existence of frequency band gap in a microstructure model. Math Mech Complex Syst2022; 9(4): 413–421.
111.
BersaniAMCaressaPCiallellaA. Numerical evidence for the approximation of dissipative systems by gyroscopically coupled oscillator chains. Math Mech Complex Syst2022; 10(3): 265–278.
112.
Osorio-BlancoERBergueiroJAbaliBE, et al. Effect of core nanostructure on the thermo-mechanical properties of soft nanoparticles. Chem Mater2020; 32: 518–528.
113.
ZhangJYangHAbaliBE, et al. Dynamic mechanics-modulated hydrogels to regulate the differentiation of stem-cell spheroids in soft microniches and modeling of the nonlinear behavior. Small2019; 15: 1901920.
114.
WagnerHBarrJZadnikK, et al. Collaborative longitudinal evaluation of keratoconus (CLEK) study: methods and findings to date. Cont Lens Anterior Eye2007; 30(4): 223–232.
115.
FedeleRRakaBHildF, et al. Identification of adhesive properties in glare assemblies using digital image correlation. J Mech Phys Solids2009; 57(7): 1003–1016.
116.
FedeleRGalantucciLCianiA. Global 2D digital image correlation for motion estimation in a finite element framework: a variational formulation and a regularized, pyramidal, multi-grid implementation. Int J Numer Meth Eng2013; 96(12): 739–762.
117.
FedeleR. Simultaneous assessment of mechanical properties and boundary conditions based on digital image correlation. Exp Mech2015; 55: 139–153.
118.
ClaireDHildFRouxS. A finite element formulation to identify damage fields: the equilibrium gap method. Int J Numer Method Eng2004; 61(2): 189–208.
119.
SpagnuoloMYildizdagMEPinelliX, et al. Out-of-plane deformation reduction via inelastic hinges in fibrous metamaterials and simplified damage approach. Math Mech Solids2022; 27(6): 1011–1031.
120.
ValmalleMSmaniottoBSpagnuoloM, et al. Mesoscale DVC analyses and parameter calibration for pantographic block in 3-point flexure. Eur J Mech A/Solids2023; 101: 105063.
121.
ValmalleMVintacheASmaniottoB, et al. Local–global DVC analyses confirm theoretical predictions for deformation and damage onset in torsion of pantographic metamaterial. Mech Mater2022; 172: 104379.
122.
MaierGBocciarelliMBolzonG, et al. Inverse analyses in fracture mechanics. Int J Fracture2006; 138(1): 47–73.
123.
FedeleRPlacidiLFabbrocinoF. A review of inverse problems for generalized elastic media: formulations, experiments, synthesis. Continuum Mech Therm2024; 36: 1413–1453.
124.
SuttonMHildF. Recent advances and perspectives in digital image correlation. Exp Mech2015; 55(1): 1–8.
125.
HildFRouxS. Digital image correlation: from displacement measurement to identification of elastic properties–a review. Strain2006; 42(2): 69–80.
126.
GrédiacMHildF. Full-field measurements and identification in solid mechanics. Hoboken, NJ: John Wiley, 2012.
127.
GermainP. The method of virtual power in the mechanics of continuous media, I: second-gradient theory. Math Mech Complex Syst2020; 8(2): 153–190.
128.
AlibertJJSeppecherPDell’IsolaF. Truss modular beams with deformation energy depending on higher displacement gradients. Math Mech Solids2003; 8(1): 51–73.
129.
BîrsanMShiraniMSteigmannDJ. Convexity conditions for fiber-reinforced elastic shells. Math Mech Solids. Epub ahead of print August 20, 2024. DOI: 10.1177/10812865241261485.
130.
BîrsanMShiraniMSteigmannDJ. The coupled Legendre-Hadamard condition for fiber-reinforced materials: three-dimensional solids and two-dimensional shells. Continuum Mech Therm2025; 37(2): 26.
131.
ShiraniMLuoCSteigmannDJ. Cosserat elasticity of lattice shells with kinematically independent flexure and twist. Continuum Mech Therm2019; 31: 1087–1097.
132.
OliveMAuffrayN. Symmetry classes in piezoelectricity from second-order symmetries. Math Mech Complex Syst2021; 9(1): 77–105.
SpagnuoloMYildizdagMEAndreausU, et al. Are higher-gradient models also capable of predicting mechanical behavior in the case of wide-knit pantographic structures?Math Mech Solids2021; 26(1): 18–29.
135.
Erden YildizdagMPlacidiLTurcoE. Modeling and numerical investigation of damage behavior in pantographic layers using a hemivariational formulation adapted for a Hencky-type discrete model. Continuum Mech Therm2023; 35(4): 1481–1494.
136.
ValorosoNFedeleR. Characterization of a cohesive-zone model describing damage and de-cohesion at bonded interfaces. sensitivity analysis and mode-I parameter identification. Int J Solids Struct2010; 47(13): 1666–1677.
137.
PonomarevD. A generalised time-evolution model for contact problems with wear and its analysis. Math Mech Complex Syst2022; 10(3): 279–319.
