INFORMATION ABOUT PROJECT,
SUPPORTED BY RUSSIAN SCIENCE FOUNDATION

The information is prepared on the basis of data from the information-analytical system RSF, informative part is represented in the author's edition. All rights belong to the authors, the use or reprinting of materials is permitted only with the prior consent of the authors.

 

COMMON PART


Project Number22-11-00087

Project titleAdvanced theoretical models and methods of studying elastic behavior and properties of nanostructures

Project LeadGrekov Mikhail

AffiliationFederal State Budgetary Educational Institution of Higher Education "Saint-Petersburg State University",

Implementation period 2022 - 2024 

Research area 01 - MATHEMATICS, INFORMATICS, AND SYSTEM SCIENCES, 01-301 - Elasticity theory, conjugated models

KeywordsElasticity, surface energy, surface defects, nanosized heterogeneities, pentagonal crystals, stress concentration, nanoplate, stability, diffusion, boundary perturbation method, integral equation


 

PROJECT CONTENT


Annotation
The object of the project is to develop and construct advanced theoretical methods for researching models of deformable solids containing nanostructures such as pores, inclusions, inhomogeneous, growths, relief regions of a surface and interface, and bodies having at least one dimension of a nanometer size (nanoplate, nanostrip, nanowire). A surface relief and roughness observed in experiments is a source of a higher stress concentration, the estimation of which demands new methods of the solution of corresponding boundary problems. Particularly, it is relevant to the nanomaterials and nanostructures which have the volume of nanometer thickness superficial layers comparable to the volume of a whole body. One of the actual problems of nanomechanics is the precise behavior of thin-walled nanostructures under bending, buckling and vibration taking into account surface properties. The relaxation process increasing the free surface of pentagonal crystals due to the diffusion void growth is of a particular interest. In the project, it is supposed to apply special methods to the solution of boundary value problems, using certain surface elasticity models for investigating elastic fields in subsurface layers and material properties when nanodefects exist on a free surface and interface. It is planned to study more completely an influence both of the surface elastic properties, and residual surface stress, on processes distorting the shape of the surface. In most publications, this stress either is not considered or considered in simplified versions. The special part of the project is methods allowing to predict an arising and evolving a relief of a free surface and interface at the nanoscale due to surface diffusion induced by the temperature and force fields. It is planned to fill the gap existing in the works published on bending models. All known models do not consider a whole system of parameters defining properties of surface layers in solids (surface tension like a liquid, the conjunction condition in the transverse direction, bending properties of a surface). Within the framework of this project, it is planned to develop new theoretical approaches to analyze the kinetics aspects of void nucleation and growth in the pentagonal wires (whiskers) under residual inhomogeneous elastic stress via bulk and grain- boundary vacancy diffusion. The results of scheduled investigations in the project are important to creating optoelectronic devices, accelerometers, frequency meters, biosensors, the work of which is especially sensitive to distortion of an external surface and surface between heterogeneous nanostructures.

Expected results
The basic results of the project are the modification of an existing mathematical methods and development of new methods for analytical solutions of the problems at the nano- and microscale for deformable bodies containing nanostructures such as pores, inclusions, inhomogeneous, growths, relief regions of a surface and interface, and bodies having at least one dimension of a nanometer size (nanoplate, nanostrip, nanowire), obtaining these solutions and analyzing them using numerical calculations. Analytical solutions of some problems either have been obtained or have been reduced by the participants of this project to the possibility of a numerical realization, using simplified Gurtin-Merdok constitutive equations of surface elasticity. It is planned in the project to construct and analyze new solutions of these problems, based on the exact Gurtin-Merdok relations in the original statement considering all terms with the residual surface stresses. The level of proposed investigations is not under the world-wide one and as for using specially developed methods allowing to consider all surface stresses arising at an arbitrarily distorted internal and external surfaces of nanomaterials by studying an elastic and strength properties of nanostructures, diffusion process of forming surface geometry, stability of thin-walled nanostructures, construct theoretical models and analytical solutions of equations of the void evolution in pentagonal crystals, analyze kinetics of origin and growth of an internal void in the field of heterogeneous elastic stresses for the cases of volume and grain-boundary diffusion, it exceeds the world-wide level. Renovating and embedding the sufficiently general method of solving new boundary two-dimensional problems of elasticity theory for solids with surface and interface nanodefects allows obtaining analytical solutions for comprehensive investigations of material properties. By applying to the nanomechanical systems, this method can be used for more precise computations of stresses in different planes and cylindrical nanostructures which are widely applied in the modern optoelectronics. Such computations are required for developing theoretical models of formation of different defects in heterogeneous nanostructures. The results obtaining during an investigation of a nanostructure formation mechanism on a surface/interface of a solid and the elastic field in the vicinity of these structures in an implementation process of the project can be found applications at the developing new structures and a formation technology of semiconductor devices which contain nanosized elements and thermo-defensive coatings, cylindrical nanosized inclusions and voids. The scientific, practical and social significance of the research of the void evolution in pentagonal nanowires is a creation of an academic basis for stress relaxation control in production and exploitation of these nanostructures. The results are expected to be used as recommendations for hollow pentagonal wires synthesis. The results will correspond to the international scientific standards and will be published in leading international academic journals of materials science, material mechanics and physics. Nonlinear models of the bending of thin-walled nanostructures are an important tool for modeling stiffness properties at the nanoscale. The construction of such models, taking into account the surface energy of a various nature (deformation and with surface tension according to Gurtin-Murdoch, with moment stresses according to Steigmann-Ogden) and comparison their results, allow us to evaluate the suitability of various approaches for solving specific problems. The planned scientific results: • Theoretical approach for studying a formation mechanism of irregularities on the free and interface boundaries, taking into account surface elasticity and surface tension, based on the complete Gurtin-Murdoch model. Investigation of the surface/interface morphological stability depending on the physical and geometric parameters. • Theoretical models describing the void kinetics under residual stress via bulk and grain boundary vacancy diffusion. The time dependencies of the void radius, the critical size of the pentagonal nanowire, the void growth time and void optimal size will be determined. • Method of the solution of the problems on finding elastic fields in the vicinity of nanosized asperities at a free plane surface and interface, using precise constitutive relations of Gurtin-Murdoch surface elasticity. Influence of geometrical and physical parameters of a surface and volume materials on the stress state of the body. • The model of nearly circular cylindrical defect (chanal, inclusion) with surface elastic properties different from similar ones of the material volume. Universal method of defining the stress-strain state in the vicinity of this defect incorporating surface stress and surface tension. It is supposed to construct an algorithm for finding stresses and displacements in any-order approximation of the boundary perturbation method. In the framework of the problem, it is supposed to study the stress concentration and distribution of stresses on the defect boundary and near it. The estimation of influence of geometrical and physical parameters of the problem will be done. • Stress field in the vicinity of a circular nanosized cavity in a plate of a nanometer thickness based on the solution of the modified Kirsch problem by the exact Gurtin-Murdoch model. By using the solution of these problems and analyzing zones of compressive stresses near the hole, the instability problem of such a plate under the uniaxial tension will be constructed. The estimation of a critical force depending on geometrical and physical parameters of the model will be given. • Nonlinear bending equations for nanoplates considering surface stress by Steigmann-Ogden. Comparative analysis the solutions of test problems considering surface stresses by Steigmann-Ogden and Gurtin-Murdoch (cylindrical bending of an infinite nanosrtip: nonlinear bending, free transverse vibration, compressive buckling and postcritical deformation of a nanoplate). The results of the Project are supposed to publish in a top of Russian and international journals such as: – International Journal of Engineering Science, Q1, IF=8.843; – Acta Materialia, Q1, IF=8.203; – Scripta Materialia, Q1, IF=5.611; – European Journal of Mechanics, Q1, IF=2.183; – Journal of Materials Science, Q1, IF=4.220; – International Journal of Solids and Structures, Q1, IF=3.9; – Mechanics of Materials, Q1, IF=3.266; – Journal of Elasticity, Q1, IF=2.372; – Mechanics Research Communications, Q1, IF=2.254; – Acta Mechanica, Q1, IF=2.698; – Continuum Mechanics and Thermodynamics, Q1, IF=3.822; – Physical Mesomechanics, Q2, IF= 1.368; – ZAMM Zeitschrift für Angewandte Mathematik und Mechanik, Q2, IF=1.603; – Engineering Fracture Mechanics, Q1, IF=4,406.


 

REPORTS


Annotation of the results obtained in 2022
As a result of studying the phenomenon of the relief formation on the surface of a nanomaterial, a generalization of the Asaro-Tiller-Grinfeld model has been proposed that allows considering the variation of the surface energy due to material deformation. The new method resulted in the analysis of the effect of the elastic constants, which characterize the elastic behavior of the surface layer, on the process of reorganization of surface atoms leading to the formation of nanoscale surface structures. The main driving force of such transformations is a variation of the chemical potential. Due to the reduced stability of the surface atomic layers, it leads to the diffusion of atoms along the surface from regions with a high value of the chemical potential to the regions with a low value. In our study, the change in the chemical potential is associated with the change in the stress field along the perturbed surface. To analyze the stresses, we developed the method based on the constitutive equations of the Gurtin-Murdoch model of surface elasticity in which the surface layer of a solid is considered as a separate substructure with mechanical properties different from the bulk material. In the frame of the research concerning the void formation and growth in multiply twinned particles the theoretical models of a void evolution in pentagonal wires and icosahedral particles are suggested to predict the pore kinetics in the dependence on the following parameters such as the elastic moduli of the particle material, the initial pore and particle radii, vacancy relaxation volume and specific surface energy. These models are based on the new analytical solutions of the boundary-value problems of vacancy diffusion induced by both the Gibbs-Thompson surface tension and the stress state of disclination defects (wedge disclination and Marks-Ioffe stereo-disclination) in hollow cylindrical and spherical bodies. The model of an elastic solid under the plane strain, containing the defect like a nearly circular cylindrical channel incorporating the surface stress and surface tension is constructed. Based on the boundary perturbation technique, the solution of the problem has been reduced to the singular integro-differential equation in the complex displacement for any order approximation. The stress fields obtained in the zero-order approximation by the original Gurtin-Murdoch surface elasticity model and its known simplified versions have been compared. The inaccuracy of each simplification is demonstrated. The same comparative analysis has been performed in the solution of the problem on an interaction of the periodic system of edge dislocations with the free plane boundary under plane strain incorporating surface stresses and surface tension. The stress fields at the boundary have been depicted depending on the distance between dislocations and the boundary using surface elasticity models under consideration. The stress field in the vicinity of a circular nanohole in a plate of a nanometer thickness obtained by the solution of the modified Kirsch problem with the original Gurtin-Murdoch model has been analyzed. The numerical results have been obtained in the form of graphic dependences of stresses on polar coordinates under the uniaxial tension of the plate. It is shown that allowing for the surface tension at the free surface of the hole and the faces of the plate essentially influence on the stress field in the vicinity of the nanohole and its surface. The Ritz method is developed for finding the critical buckling load of a stretched plate of a nanometer thickness with a circular hole, in which the obtained stress field is realized, taking into account the surface tension at the surface of the hole and on the faces of the plate in the transverse direction. The program was written in the Matlab environment. Using this program, the analysis of the results obtained was carried out depending on the parameters of the problem.

 

Publications

1. Grekov M.A., Vakaeva A.B., Müller W.H. Stress field around cylindrical nanopore by varies models of surface elasticity Continuum Mechanics and Thermodynamics, - (year - 2022) https://doi.org/10.1007/s00161-022-01168-7

2. Khramov A.S., Krasnitckii S.A., Smirnov A.M. , Gutkin M. Yu. The void evolution kinetics driven by residual stress in icosahedral particles Materials Physics and Mechanics, - (year - 2022)

3. Kostyrko S., Grekov M., Altenbach H. Stress distribution at the wavy surface of a solid incorporating surface stresses and surface tension Advances in Solid and Fracture Mechanics. Advanced Structured Materials., 180, 151-166 (year - 2022) https://doi.org/10.1007/978-3-031-18393-5_10

4. Vakaeva A., Grekov M. Stress field around cylindrical pore by various surface elasticity models European Conference on Fracture 2022. ESIS Publishing House, Madeira, Portugal, ECF23. Book of Abstracts, p. 521 (year - 2022)


Annotation of the results obtained in 2023
In this study, a mathematical model of a heteroepitaxial material has been developed to allow for the influence of geometric and physical-mechanical properties of an interfacial layer on the distribution of elastic stresses along the nanostructured interface between two bulk materials under plane strain conditions in accordance with the Gurtin–Murdoch model of surface/interfacial elasticity. Based on the boundary perturbation technique, we derived a solution to the boundary value problem, resulting in the obtaining of asymptotic representations for the stress and strain tensor components at any desired level of approximation. The stress state near the nanostructured interphase of an isotropic heteroepitaxial system has been studied using the first-order approximation. The effects of the perturbation wavelength, interfacial stiffness, relative stiffness of bulk phases, residual interfacial and misfit stresses have been analyzed. A comparative assessment has been conducted, juxtaposing the obtained solution with the outcomes derived from the simplified Gurtin–Murdoch model in earlier studies. A further development of the approach to studying the formation of surface irregularities in the stressed solids under the influence of diffusion processes was proposed for an arbitrary surface relief. The analysis involves examining the stress distribution along the surface of a complex-shape which forms due to the loss of the morphological stability predicted by the extended model. To consider the impact of the film thickness and elastic properties of the substrate on the formation of the surface relief in a heteroepitaxial film coating under the influence of the surface diffusion, the approach based on the utilization of the superposition method was proposed. This method enables to determine the stressed state of the film near the perturbed free surface through a linear combination of previously obtained solutions for a homogeneous solid and a heterogeneous layered material. Applying the established methodology, we investigated how the geometric and physical-mechanical parameters of the film-substrate system impact the formation of the nanosized relief on the film surface. The theoretical research concerning the void evolution in pentagonal nanowires and icosahedral particles due to relaxation of the residual stress has been proceeded in the second stage of the project. The residual stress induced by the five-fold symmetry axes (restricted by the classical crystallography) can be described within the disclination concept. During particle growth, the elastic energy stored by the residual stress is increased until it relaxes through the void formation. The proposed energetic model of void evolution in pentagonal particles (multiply twinned particles) allows anticipating void evolution scenarios in dependence of initial void radius and material parameters (elastic constants, surface tensile and disclination strength). Besides, the critical and optimal parameters of voids in pentagonal particles are determined. Continuous models of the nonlinear buckling of nanorod and nanoplate by the Steigmann-Ogden theory of surface elasticity have been constructed. Based on the Euler–Bernoulli kinematic hypothesis and Timoshenko hypothesis, the constitutive relations for the internal forces and moments under the nonlinear buckling of the nanorod have been derived. New effective models and rigidities were introduced for the description of these relations. The dynamic equations for the nonlinear buckling of the nanorod have been derived as well. Two static problems on buckling of a nanofiber under the dead weight and free (from both sides) and rigid (from one side - cantilever) supports have been solved. The size effect has been revealed, i.e., the surface stresses increase essentially the nanofiber rigid with decreasing the nanofiber thickness to the nanometer size. Comparative analysis of applying two beam theories considered and Gurtin-Murdoch and Steigmann-Ogden surface elasticity models has been carried out. It has been discovered that the Timoshenko beam theory and Steigmann-Ogden surface elasticity soften the size effect. The solution of the problem on the stress-strain state of an elastic body containing a nearly circular hole of a nanometer size under the plane strain conditions has been first obtained. The solution is based on the complete Gurtin-Murdoch surface elasticity model and the boundary perturbation method. The algorithm of the solution is adopted for any approximation and any continually differential function describing the boundary of the hole. Obtained results allow estimating influence of an error of the shape deviation from the circular one on the stress state and improving accuracy as compared with simplified surface elasticity models. The stress field in the Kirsch problem for a plate of a nanometer thickness with a circular hole incorporating surface stresses at the boundary of the hole and on the faces of the plate has been investigated based on the complete Gurtin-Murdoch surface elasticity model. As follows from the analytical solution, the elastic stress field depends on the thickness. The difference between modified and classic (not allowing for the surface effects) problems decreases with the distance from the hole boundary. Besides, numerical studies show that the stress field approaches the Kirsch solution when the hole radius increases. The solution of the problem on an interaction of an edge dislocations row with the flat surface of an elastic body incorporating the stretching and bending rigidity of the surface material according to the Steigmann-Ogden surface elasticity model has been constructed. The comparison of the stress fields at the surface, obtained by the Steigmann-Ogden, Gurtin-Murdoch and all simplified models, using actively in the literature, has been carried out. It was found that maximum values of longitudinal stresses decrease, and normal ones increase due to allowing for the bending rigidity of the surface. An influence of the surface bending rigidity on the image force acting on each dislocation has been studied as well. The results obtained by the Gurtin-Murdoch and Steigmann-Ogden models show that allowing for the surface bending rigidity leads to the decrease of the image force regardless of the Burgers vector direction and the relative distance between neighboring dislocations. It was established that the solution can be used in the case of an isolated dislocation if the distance between neighboring dislocations is two orders of magnitude greater than the distance of the dislocations to the surface.

 

Publications

1. Bochkarev A.O. Buckling of a nano-rod with taken into account of surface effect Zeitschrift für Angewandte Mathematik und Mechanik, e202300738 (year - 2023) https://doi.org/10.1002/zamm.202300738

2. Grekov M.A. Surface-dislocation interaction by various models of surface elasticity International Journal of Engineering Science, - (year - 2024)

3. Khramov A.S., Krasnitckii S.A., Smirnov A.M., Gutkin M.Yu. A Kinetic Model of the Stress-Induced Void Evolution in Pentagonal Whiskers and Rods Materials Physics and Mechanics (Механика и физика материалов), - (year - 2023)

4. Vakaeva A.B., Grekov M.A. Modified Kirsch problem incorporating surface stresses under plane stress coupled2023, Vol. Simulation methods for coupled problems, 42 (year - 2023) https://doi.org/10.23967/c.coupled.2023.038

5. - Программа для расчета напряженного состояния наноразмерных топологических дефектов на поверхности твердого тела и оценки их морфологической устойчивости (STS-NSS) -, 2023610171 - номер свидетельства (2022685603 - номер заявки) (year - )

6. - Программа для расчета напряженного состояния и оценки морфологической устойчивости наноразмерных топологических дефектов на поверхности ультратонкого пленочного покрытия (STS-UTF) -, 2023683568 - номер свидетельства (2023682660 - номер заявки) (year - )