Control of complex turbulent dynamical systems involving strong nonlinearity and high degrees of internal instability is an important topic in practice. Different from traditional methods for controlling individual trajectories, controlling the statistical features of a turbulent system offers a more robust and efficient approach. Crude first-order linear response approximations were typically employed in previous works for statistical control with small initial perturbations. This paper aims to develop two new statistical control strategies for scenarios with more significant initial perturbations and stronger nonlinear responses, allowing the statistical control framework to be applied to a much wider range of problems. First, higher-order methods, incorporating the second-order terms, are developed to resolve the full control-forcing relation. The corresponding changes to recovering the forcing perturbation effectively improve the performance of the statistical control strategy. Second, a mean closure model for the mean response is developed, which is based on the explicit mean dynamics given by the underlying turbulent dynamical system. The dependence of the mean dynamics on higher-order moments is closed using linear response theory but for the response of the second-order moments to the forcing perturbation rather than the mean response directly. The performance of these methods is evaluated extensively on prototype nonlinear test models, which exhibit crucial turbulent features, including non-Gaussian statistics and regime switching with large initial perturbations. The numerical results illustrate the feasibility of different approaches due to their physical and statistical structures and provide detailed guidelines for choosing the most suitable method based on the model properties.
{"title":"Effective statistical control strategies for complex turbulent dynamical systems","authors":"Jeffrey Covington, Di Qi, Nan Chen","doi":"10.1098/rspa.2023.0546","DOIUrl":"https://doi.org/10.1098/rspa.2023.0546","url":null,"abstract":"Control of complex turbulent dynamical systems involving strong nonlinearity and high degrees of internal instability is an important topic in practice. Different from traditional methods for controlling individual trajectories, controlling the statistical features of a turbulent system offers a more robust and efficient approach. Crude first-order linear response approximations were typically employed in previous works for statistical control with small initial perturbations. This paper aims to develop two new statistical control strategies for scenarios with more significant initial perturbations and stronger nonlinear responses, allowing the statistical control framework to be applied to a much wider range of problems. First, higher-order methods, incorporating the second-order terms, are developed to resolve the full control-forcing relation. The corresponding changes to recovering the forcing perturbation effectively improve the performance of the statistical control strategy. Second, a mean closure model for the mean response is developed, which is based on the explicit mean dynamics given by the underlying turbulent dynamical system. The dependence of the mean dynamics on higher-order moments is closed using linear response theory but for the response of the second-order moments to the forcing perturbation rather than the mean response directly. The performance of these methods is evaluated extensively on prototype nonlinear test models, which exhibit crucial turbulent features, including non-Gaussian statistics and regime switching with large initial perturbations. The numerical results illustrate the feasibility of different approaches due to their physical and statistical structures and provide detailed guidelines for choosing the most suitable method based on the model properties.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":"67 3","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135515145","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Honeycomb materials frequently encounter hypergravity conditions in both aerospace and biological contexts during phases such as launch, reentry or under centrifugal motion. The significant body force engendered by hypergravity induces alterations in the microstructure of honeycomb materials, which in turn, influences their macroscopic mechanical behaviour. Leveraging the stiffness of the beam element as a pivotal variable, we successfully derived the equivalent moduli of the honeycomb material under hypergravity conditions. We further proposed the concept of a ‘hypergravity factor’, elucidating that the density of the base material, the dimensions of honeycomb cells and the magnitude of the hypergravity contribute to amplifying hypergravity effects. The results, numerically validated through finite-element simulations, could be reduced to the case that neglects body force. The critical buckling load of the honeycomb material under hypergravity can be assessed by setting the derived moduli to zero. In the presence of hypergravity, a honeycomb material undergoes a transition into a gradient material along the hypergravity direction, thereby exacerbating anisotropy. This phenomenon is theoretically expected to occur in virtually all porous materials. The analytical framework adopted, which employs beam stiffness as an intermediary variable, facilitates the extension of these results to honeycomb materials which encompass beam elements with functional gradients or varying cross-sectional morphologies.
{"title":"Equivalent in-plane elastic moduli of honeycomb materials under hypergravity conditions","authors":"Lei Wang, Guannan Wang, Chaofeng Lü","doi":"10.1098/rspa.2023.0638","DOIUrl":"https://doi.org/10.1098/rspa.2023.0638","url":null,"abstract":"Honeycomb materials frequently encounter hypergravity conditions in both aerospace and biological contexts during phases such as launch, reentry or under centrifugal motion. The significant body force engendered by hypergravity induces alterations in the microstructure of honeycomb materials, which in turn, influences their macroscopic mechanical behaviour. Leveraging the stiffness of the beam element as a pivotal variable, we successfully derived the equivalent moduli of the honeycomb material under hypergravity conditions. We further proposed the concept of a ‘hypergravity factor’, elucidating that the density of the base material, the dimensions of honeycomb cells and the magnitude of the hypergravity contribute to amplifying hypergravity effects. The results, numerically validated through finite-element simulations, could be reduced to the case that neglects body force. The critical buckling load of the honeycomb material under hypergravity can be assessed by setting the derived moduli to zero. In the presence of hypergravity, a honeycomb material undergoes a transition into a gradient material along the hypergravity direction, thereby exacerbating anisotropy. This phenomenon is theoretically expected to occur in virtually all porous materials. The analytical framework adopted, which employs beam stiffness as an intermediary variable, facilitates the extension of these results to honeycomb materials which encompass beam elements with functional gradients or varying cross-sectional morphologies.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":"39 2","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135764574","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Some of the fundamentals of quantum information science are described, including the concepts of quantum resources, quantum states and mixedness of states. The explanations are detailed and include a combination of basic facts with fully worked examples, and some more advanced topics. The principles of quantum information are illustrated with chemical examples drawn from singlet fission, photophysics of radicals, molecular excitons, electron transfer and so on. Suggestions for prospects and challenges for the field are discussed.
{"title":"A molecular perspective on quantum information","authors":"Gregory D. Scholes","doi":"10.1098/rspa.2023.0599","DOIUrl":"https://doi.org/10.1098/rspa.2023.0599","url":null,"abstract":"Some of the fundamentals of quantum information science are described, including the concepts of quantum resources, quantum states and mixedness of states. The explanations are detailed and include a combination of basic facts with fully worked examples, and some more advanced topics. The principles of quantum information are illustrated with chemical examples drawn from singlet fission, photophysics of radicals, molecular excitons, electron transfer and so on. Suggestions for prospects and challenges for the field are discussed.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":"67 7","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135515357","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The energy transfer process for acoustic propagation and scattering in inhomogeneous, lossy media is rigorously investigated. Two problems are considered: an arbitrarily composed cluster of lossy scatterers, excited by a single point source and a layered, lossy medium excited by an internal point-source distribution. Energy conservation laws and optical theorems are established. The connection of active intensity with the scattering cross section and the relation between reactive intensity and Lagrangian density are determined. Furthermore, the influence of the losses of the involved media to the kinetic energy is unveiled. Reductions to lossless media and near-zero frequency values are derived.
{"title":"Analysis of the energy transfer process for multiple scattering problems involving lossy media","authors":"Andreas Kalogeropoulos, Nikolaos L. Tsitsas","doi":"10.1098/rspa.2023.0513","DOIUrl":"https://doi.org/10.1098/rspa.2023.0513","url":null,"abstract":"The energy transfer process for acoustic propagation and scattering in inhomogeneous, lossy media is rigorously investigated. Two problems are considered: an arbitrarily composed cluster of lossy scatterers, excited by a single point source and a layered, lossy medium excited by an internal point-source distribution. Energy conservation laws and optical theorems are established. The connection of active intensity with the scattering cross section and the relation between reactive intensity and Lagrangian density are determined. Furthermore, the influence of the losses of the involved media to the kinetic energy is unveiled. Reductions to lossless media and near-zero frequency values are derived.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":"167 4","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135161363","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We derive a Sturm–Liouville system of equations for the exact calculation of the survival probability in first-passage time problems. This system is the one associated with the Wiener–Hopf integral equation obtained from the theory of random walks. The derived approach is an alternative to the existing literature and we tested it against direct calculations from both discrete- and continuous-time random walks in a manageable, but meaningful, example. Within this framework, the Sparre Andersen theorem results to be a boundary condition for the system.
{"title":"Sturm–Liouville systems for the survival probability in first-passage time problems","authors":"Marcus Dahlenburg, Gianni Pagnini","doi":"10.1098/rspa.2023.0485","DOIUrl":"https://doi.org/10.1098/rspa.2023.0485","url":null,"abstract":"We derive a Sturm–Liouville system of equations for the exact calculation of the survival probability in first-passage time problems. This system is the one associated with the Wiener–Hopf integral equation obtained from the theory of random walks. The derived approach is an alternative to the existing literature and we tested it against direct calculations from both discrete- and continuous-time random walks in a manageable, but meaningful, example. Within this framework, the Sparre Andersen theorem results to be a boundary condition for the system.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":"38 4","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135764575","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The crosslinking technique (CXL) is an effective low-risk therapeutic treatment of keratoconus and other ectatic disorders of the human cornea. The effect of corneal CXL is to increase the stiffness of the stroma to prevent the progression of the cornea distortion. Several clinical and experimental studies have shown that the stiffening effects predominantly localize on the anterior portion of the stroma and that the in-depth stiffening distribution is highly dependent on the duration of treatment. Yet, how the stiffening effects distribute through the cornea thickness as a function of the treatment duration is an open question. Here, we propose an analytical model of the stiffening profile due to CXL treatment as a function of the irradiation time. We consider linear and nonlinear variations of the crosslinking effects across the thickness and implement them into a finite element model of the porcine cornea. We present a time-dependent in-depth stiffening profile that allows us to predict the post-operative cornea response to physiological intraocular pressure for different irradiation times. We anticipate that this predictive model will support the development of patient specific three-dimensional models that will allow clinicians to design customized CXL treatment, thus enhancing treatment outcomes.
{"title":"A predictive model of UV-A-riboflavin crosslinking treatment on porcine corneas","authors":"Alessandra Bonfanti, Anna Pandolfi","doi":"10.1098/rspa.2023.0323","DOIUrl":"https://doi.org/10.1098/rspa.2023.0323","url":null,"abstract":"The crosslinking technique (CXL) is an effective low-risk therapeutic treatment of keratoconus and other ectatic disorders of the human cornea. The effect of corneal CXL is to increase the stiffness of the stroma to prevent the progression of the cornea distortion. Several clinical and experimental studies have shown that the stiffening effects predominantly localize on the anterior portion of the stroma and that the in-depth stiffening distribution is highly dependent on the duration of treatment. Yet, how the stiffening effects distribute through the cornea thickness as a function of the treatment duration is an open question. Here, we propose an analytical model of the stiffening profile due to CXL treatment as a function of the irradiation time. We consider linear and nonlinear variations of the crosslinking effects across the thickness and implement them into a finite element model of the porcine cornea. We present a time-dependent in-depth stiffening profile that allows us to predict the post-operative cornea response to physiological intraocular pressure for different irradiation times. We anticipate that this predictive model will support the development of patient specific three-dimensional models that will allow clinicians to design customized CXL treatment, thus enhancing treatment outcomes.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":"45 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135764547","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In Tabula XCI verso of Codex Atlanticus, Leonardo da Vinci presents an ingenious masonry structure composed of segments in the shape of inverted triangles. These are assembled by contact in a chain to obtain a lintel or jack arch, where they are pressed together by the thrust of the end constraints. Drawing inspiration from Leonardo’s sketches, we show that, by connecting the segments in pair through elastic tendons, this system represents a new type of flextegrity beam. In a classical flextegrity, the contact surfaces of the segments are curved conjugate profiles, imposing a pure rolling motion along properly designed pitch lines: the consequent elongation of the tendon dictates the constitutive response as a function of the relative rotation of the segments. Here, the contact is through plane surfaces, so that the kinematics, besides the relative rotation, is characterized by segmental shearing. This system is theoretically analysed and a continuum model is derived as the length of the segments becomes small. Comparisons with experiments on three-dimensional-printed prototypes confirm the theoretical findings and highlight the possible competition between rotational and sliding types of deformation. Apart from the historical value, this type of construction can be used in innovative structures or metamaterials.
{"title":"Shear and flexural deformations in flextegrity segmental beams inspired by Leonardo’s triangular masonry construction","authors":"Claudio Boni, Gianni Royer-Carfagni","doi":"10.1098/rspa.2023.0453","DOIUrl":"https://doi.org/10.1098/rspa.2023.0453","url":null,"abstract":"In Tabula XCI verso of Codex Atlanticus, Leonardo da Vinci presents an ingenious masonry structure composed of segments in the shape of inverted triangles. These are assembled by contact in a chain to obtain a lintel or jack arch, where they are pressed together by the thrust of the end constraints. Drawing inspiration from Leonardo’s sketches, we show that, by connecting the segments in pair through elastic tendons, this system represents a new type of flextegrity beam. In a classical flextegrity, the contact surfaces of the segments are curved conjugate profiles, imposing a pure rolling motion along properly designed pitch lines: the consequent elongation of the tendon dictates the constitutive response as a function of the relative rotation of the segments. Here, the contact is through plane surfaces, so that the kinematics, besides the relative rotation, is characterized by segmental shearing. This system is theoretically analysed and a continuum model is derived as the length of the segments becomes small. Comparisons with experiments on three-dimensional-printed prototypes confirm the theoretical findings and highlight the possible competition between rotational and sliding types of deformation. Apart from the historical value, this type of construction can be used in innovative structures or metamaterials.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":"37 4-5","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135765803","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Developable degree-4 (DD4) vertices have four facets and four creases and can be unfolded flat. The rigid-folding kinematics of DD4 vertices is rich in that it generally has two folding modes and can get stuck when two facets bind together. To study the full spectrum of the kinematics of DD4 vertices, parametric solutions for fold angles in terms of the cotangents of half-angles are derived from the opposite and adjacent fold angle relationships. It is shown that any two fold angles of a general DD4 vertex are related by the equation of a hyperbola. When the vertex has collinear creases or is flat-foldable, the pertinent hyperbola equations degenerate into linear relationships. Meanwhile, when DD4 vertices are classified into three categories according to Grashof’s criterion, both unique and binding folds can be readily located from the facet with the largest or smallest sector angle. The rigid-folding kinematics of typical vertices is then investigated. In addition to the flat state, the two folding modes can also be switched at the binding states if self-intersection is permitted. The results provide new formulae and clear geometric views on the rigid-folding kinematics of DD4 vertices, which are fundamental for constructing larger origami patterns.
{"title":"Parametric solutions to the kinematics of developable degree-4 rigid origami vertices","authors":"Yucai Hu, Changjun Zheng, Chuanxing Bi, Haiyi Liang","doi":"10.1098/rspa.2023.0319","DOIUrl":"https://doi.org/10.1098/rspa.2023.0319","url":null,"abstract":"Developable degree-4 (DD4) vertices have four facets and four creases and can be unfolded flat. The rigid-folding kinematics of DD4 vertices is rich in that it generally has two folding modes and can get stuck when two facets bind together. To study the full spectrum of the kinematics of DD4 vertices, parametric solutions for fold angles in terms of the cotangents of half-angles are derived from the opposite and adjacent fold angle relationships. It is shown that any two fold angles of a general DD4 vertex are related by the equation of a hyperbola. When the vertex has collinear creases or is flat-foldable, the pertinent hyperbola equations degenerate into linear relationships. Meanwhile, when DD4 vertices are classified into three categories according to Grashof’s criterion, both unique and binding folds can be readily located from the facet with the largest or smallest sector angle. The rigid-folding kinematics of typical vertices is then investigated. In addition to the flat state, the two folding modes can also be switched at the binding states if self-intersection is permitted. The results provide new formulae and clear geometric views on the rigid-folding kinematics of DD4 vertices, which are fundamental for constructing larger origami patterns.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":"70 4","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135515334","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We consider the release and subsequent gravity-driven spreading of a dense finite volume of fluid in an anisotropic porous medium bounded by an impermeable substrate. When the permeability in the vertical direction is much smaller than in the horizontal direction, as is the case in many real geological reservoirs, this restricts the spread of the current to a very thin layer near the impermeable base. Using a combination of asymptotic analysis and finite difference computations of Darcy flow, we show that there exist two distinct flow regimes. At early times, the bulk of the current descends slowly and uniformly, injecting fluid into thin finger-like regions near the base. At much later times, the current transitions to the classical gravity-driven solution and continues to spread with a self-similar shape. One interesting consequence is that the swept volume of the current grows differently depending on the anisotropy of the medium. This has important consequences for managing contaminant spills, where it is important to minimize the contacted volume of the aquifer, or during geological CO2 sequestration where a larger contacted volume results in more CO2 being stored.
{"title":"Anisotropy distorts the spreading of a fixed volume porous gravity current","authors":"Graham P. Benham","doi":"10.1098/rspa.2023.0271","DOIUrl":"https://doi.org/10.1098/rspa.2023.0271","url":null,"abstract":"We consider the release and subsequent gravity-driven spreading of a dense finite volume of fluid in an anisotropic porous medium bounded by an impermeable substrate. When the permeability in the vertical direction is much smaller than in the horizontal direction, as is the case in many real geological reservoirs, this restricts the spread of the current to a very thin layer near the impermeable base. Using a combination of asymptotic analysis and finite difference computations of Darcy flow, we show that there exist two distinct flow regimes. At early times, the bulk of the current descends slowly and uniformly, injecting fluid into thin finger-like regions near the base. At much later times, the current transitions to the classical gravity-driven solution and continues to spread with a self-similar shape. One interesting consequence is that the swept volume of the current grows differently depending on the anisotropy of the medium. This has important consequences for managing contaminant spills, where it is important to minimize the contacted volume of the aquifer, or during geological <mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\"> <mml:msub> <mml:mrow> <mml:mi mathvariant=\"normal\">CO</mml:mi> </mml:mrow> <mml:mn>2</mml:mn> </mml:msub> </mml:math> sequestration where a larger contacted volume results in more <mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\"> <mml:msub> <mml:mrow> <mml:mi mathvariant=\"normal\">CO</mml:mi> </mml:mrow> <mml:mn>2</mml:mn> </mml:msub> </mml:math> being stored.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":"44 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135764551","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Plane and axisymmetric problems of interfacial cracks in a power-law graded infinite medium are examined. It is shown that the models are governed by integral equations of Mellin’s convolution type whose kernels are expressed through the hypergeometric function. Exact solutions are derived by the method of orthogonal Jacobi polynomials in a series form and by the Wiener–Hopf method by quadratures. Mode-I stress intensity factors and the associated weight functions for power-law graded materials in the plane and axisymmetric cases are introduced and evaluated. It is shown that, although the displacement jump and the normal traction component have power singularities at the crack tip different from 1/2, the strain energy variation is proportional to the crack length change as in the case of homogeneous materials. A Griffith-type criterion of crack propagation in power-law graded materials is proposed and results of numerical tests are reported.
{"title":"Plane and axisymmetric problems of an interfacial crack in a power-law graded material","authors":"Y. A. Antipov","doi":"10.1098/rspa.2023.0630","DOIUrl":"https://doi.org/10.1098/rspa.2023.0630","url":null,"abstract":"Plane and axisymmetric problems of interfacial cracks in a power-law graded infinite medium are examined. It is shown that the models are governed by integral equations of Mellin’s convolution type whose kernels are expressed through the hypergeometric function. Exact solutions are derived by the method of orthogonal Jacobi polynomials in a series form and by the Wiener–Hopf method by quadratures. Mode-I stress intensity factors and the associated weight functions for power-law graded materials in the plane and axisymmetric cases are introduced and evaluated. It is shown that, although the displacement jump and the normal traction component have power singularities at the crack tip different from 1/2, the strain energy variation is proportional to the crack length change as in the case of homogeneous materials. A Griffith-type criterion of crack propagation in power-law graded materials is proposed and results of numerical tests are reported.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":"41 2","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135764566","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}