非局部弹性理论概述及在当代材料背景下的一些进展

局部作用原理在发展经典连续介质力学的本构关系中起着重要的作用，由此导出的简单物质理论得到了广泛的应用．然而，随着科技的发展，各种具有微结构的新材料不断涌现，理论和实验表明，非局部理论可以更好地刻画这些材料的宏观力学行为．本文简要介绍了一些传统的非局部弹性理论，包括Eringen 理论、Kunin 理论、Mindlin 理论；阐述了针对复合材料发展的，具有时间-空间非局部特征的Willis 方程、最新的时间-空间耦合非局部弹性动力学理论以及近场动力学理论．时间-空间非局部理论反映了复合材料宏观性能固有的非局部特征，而具有空间非局部特征的近场动力学理论便于处理具有不连续性的问题．最后，本文讨论了非局部理论的发展中值得关注的一些问题.     

Nonlocal stress field of interface dislocations

Summary The basic theory of nonlocal elasticity is stated with emphasis on the difference between the nonlocal theory and classical continuum mechanics. The concept of Nonlocal Interface Residual (NIR) is introduced in nonlocal theory. With the concept of NIR and the nonlocal constitutive equation, we calculate nonlocal stresses due to an edge dislocation on the interface of bi-materials. The nonlocal stress distribution along an interface is quite different from the classical one. Instead of the singularity in the dislocation core, nonlocal stress gives a finite value in the core. A maximum of the stress is also found near the dislocation core. Received 27 May 1997; accepted for publication 1 July 1997     

Crack growth in planar elastic fiber materials

Particularly attention is here given to crack growth in opening mode in fiber networks. Low- and high-density cellulose fiber materials are used in synchrotron X-ray microtomography tensile experiments to illustrate phenomena arising during crack growth. To capture the observed fundamental mechanisms, significantly different from classical continua, a mechanical model based on a strong nonlocal theory is applied in which an intrinsic length reflects a characteristic length of the microstructure. Nonlocal stress and strain tensor fields are estimated by analytical solutions on closed form to a modified inhomogeneous Helmholtz equation using LEFM crack-tip fields as source terms. Justified by experimental observations, physical requirements of finite stresses and strains at infinity and at the tip are applied to remove singularities. The near-tip nonlocal hoop stress field is used to estimate crack growth directions and sizes of fracture process zones. Experimental observations are shown to be qualitatively well in accordance with numerical predictions, which justifies the adopted approach.     

Cosserat modeling of cellular solidsModélisation des solides cellulaires selon Cosserat

Cellular solids inherit their macroscopic mechanical properties directly from the cellular microstructure. However, the characteristic material length scale is often not small compared to macroscopic dimensions, which limits the applicability of classical continuum-type constitutive models. Cosserat theory, however, offers a continuum framework that naturally features a length scale related to rotation gradients. In this paper a homogenization procedure is proposed that enables the derivation of macroscopic Cosserat constitutive equations based on the underlying microstructural morphology and material behavior. To cite this article: P.R. Onck, C. R. Mecanique 330 (2002) 717–722.     

A review on the application of modified continuum models in modeling and simulation of nanostructures

Analysis of the mechanical behavior of nanostructures has been very challenging. Surface energy and nonlocal elasticity of materials have been incorporated into the traditional continuum analysis to create modified continuum mechanics models. This paper reviews recent advancements in the applications of such modified continuum models in nanostructures such as nanotubes, nanowires, nanobeams,graphenes, and nanoplates. A variety of models for these nanostructures under static and dynamic loadings are mentioned and reviewed. Applications of surface energy and nonlocal elasticity in analysis of piezoelectric nanomaterials are also mentioned. This paper provides a comprehensive introduction of the development of this area and inspires further applications of modified continuum models in modeling nanomaterials and nanostructures.     

Evolution of bond fractures in a randomly distributed fiber network

Fracture in a planar randomly ordered fiber network subjected to approximately homogenous macroscopic stress and strain field is considered. A theory describing material degradation on a macroscopic scale is derived via Griffith’s energy balance for an internal fractured area in the network assuming the active fracture process on the microscopic level is fiber–fiber bond breakage. Attention is confined to a purely mechanical theory assuming isothermal processes and the theory relies on equations commonly used in theories of statistical physics. In the theory, a bond breaking driving force is stated to be equal to the elastic strain energy density of a non-fractured network. A debond fraction can be coupled to a linearly decrease of the network’s macroscopic stiffness. The rate of the fracture processes is determined by the network’s inherent properties (bond and fiber density, bond strength, etc.). During the loading process, until onset of localization, the bond breaks occur at randomly distributed locations spread over the fiber network and the theory estimate material degradation on a macroscopic level. When localization takes place, the fracture process changes from a two-dimensional randomly distributed process to a one-dimensional process and other theories have to be included to describe post-localization behavior. An approximately in-plane isotropic low-density paper is used in tensile experiments while monitoring acoustic emission activity to evaluate the theory. The experimentally obtained results support the theory surprisingly well.     

On the coupling between macroscopic material degradation and interfiber bond fracture in an idealized fiber network

The numerical analysis performed here, using a finite element network model, provides a number of important results regarding the evolution of micro fractures in planar random fiber networks where the only active microscopic fracture mechanism is bond fracture. The fibers are randomly distributed in the network meaning that the network is considered having in-plane isotropic properties on the macroscopic scale. The network is loaded so that, in an average sense, homogenous macroscopic stress and strain fields are present.Several conclusions are drawn. It is found that the development of macroscopic material degradation follows an exponential two-parameter law, consisting of an onset parameter and a fracture rate parameter, justifying a previous theory derived by the authors. The fracture rate parameter is linearly related to the inverse of the bond density above a certain density limit (percolation) and increases with increasing slenderness ratio of the fibers when keeping the bond density at a constant level. The strain energies stored in interfiber bonds are exponentially distributed over the whole network. The numerical analysis reveals that there is a linear relation between the ratio of fractured and initial number of loaded bonds, and the network’s macroscopic material stiffness normalized with its pristine stiffness, confirming earlier findings based on experimental observations. At localization the analyzed theory looses its validity because the fracture process is no longer randomly distributed over the whole network. Localization coincides with location of peak load in force–displacement tensile tests.     

Elasto-plastic analysis of influences of bond deformability on the mechanical behavior of fiber networks

Micromechanical modeling of three-dimensional (3D) fiber networks is performed by reducing the web structure to fiber segments and fiber–fiber bonds to explore the influence of fiber–fiber bond deformability on fiber network’s elasto-plasticity. The fiber segment between every two adjacent bonds is described by a Timoshenko beam element, while fiber–fiber bonds are taken as a different two-node element due to the extremely large height-to-span ratio. This overcomes the rigid-bond assumption employed by most previous network models, providing the feasibility to build the relationship between bonding condition and global mechanical properties, resulting in a 3D network structure and able to accommodate curled fibers. Both fibers and bonds are assumed to be elasto-plastic and described by the bilinear kinematic hardening model. Deformation of the network with load increasing is simulated by the Newton–Raphson method. Numerical tensile tests show reasonable results and agree qualitatively with experiments. The influence of bond deformability on the mechanical behavior of the network is discussed in detail.     

微结构非局部效应的多重性模型及在微梁弯曲中的应用

基于Eringen非局部弹性理论,直接利用逐次逼近法推导了非局部应力场的精确表达,该精确的非局部应力可具体表示为一个无穷级数的形式. 然后以微梁的横向弯曲和纯弯曲变形为例,建立平衡方程并求解及分析了挠度受非局部效应的影响. 结果表明：根据所取非局部小尺度参数大小的不同,非局部微梁的弯曲挠度可低于也可以高于经典力学下的挠度,非局部效应的增大可提高亦可降低结构的抗弯刚度. 本文结果证明了Wang以及Lim等人分别提出的两种相反的非局部模型的各自正确性. 同时首次发现,弯曲挠度随着非局部效应的增大而上下波动且存在若干跳跃点,挠度是非局部小尺度参数的非单调函数,研究同时给出了一种确定材料非局部常数的建议途径.     

A Note on Fracture Models Incorporating Surface Elasticity

Several recent papers that describe an approach to developing a theory of brittle fracture incorporating crack-surface elasticity have asserted that such models predict bounded stresses and strains at crack-tips in two dimensions. It is proved in this short note that the assertion is incorrect. While incorporating surface elasticity into fracture models in the setting of the linearized theory of elasticity does remove the square-root stress/strain singularity at crack-tips seen in classical theories that do not take account of crack-surface interfacial effects, it does not predict bounded stresses and strains but rather diminishes the singularity to being logarithmic.     

近场动力学研究进展与岩石破裂过程模拟

近场动力学（Peridynamics，PD）作为一种新兴的非局部性理论，在非连续处不需要任何处理，能够很好表述模型从连续到非连续的过程。首先，在PD基本理论简介的基础上，系统回顾了PD的国内外研究现状。其次，采用键型PD理论对非均匀性的圆孔岩板单轴拉伸破裂过程进行了二维数值模拟，采用态型PD理论对单轴、常规三轴以及真三轴等不同压缩条件下的岩石破裂过程进行了三维数值模拟，并以加拿大Mine-by隧洞为例对现场岩体破裂过程进行了模拟，结果表明PD在岩石破裂过程模拟上具有较强适用性。最后，指出当前PD在岩石破裂过程模拟中存在的主要问题和未来值得开展的若干研究课题。     

具有非局部体力矩的非局部弹性理论

本文基于非局部连续统场论的公理系统,建立了具有非局部体力矩作用的非局部弹性理论,我们证明了,在非局部弹性固体中存在着非局部体力矩,非局部体力矩引起了应力的非对称和非局部体力矩是由材料中的共价键产生的。     

A mixed Von Mises distribution for modeling soft biological tissues with two distributed fiber properties

Constitutive modeling of biological tissues plays an important role in the understanding of tissue behavior and the development of synthetic materials for medical and bio-inspired applications. A structural continuum model that incorporates principal structural features of the tissue can potentially provide the link between microstructure and the macroscopic mechanical response of biological tissues. For most soft biological tissues, including arterial walls and skin tissue, the main load-carrying constituent is presumed to be the distributed collagen fibers embedded in a base matrix. It is believed that the organization of the collagen fibers gives rise to the anisotropy of the material. In this paper, a semi-structural constitutive model is proposed to account for planar fiber distributions with more than one distributed planar fiber property. Motivated by histology information of the wing membrane of the bat, a statistical treatment is formulated in this paper to capture the overall effect of the distribution of fiber cross-sectional area and the distribution of the number of fibers. This formulation is suitable for general cases when more than one fiber property varies spatially. Furthermore, this model is a two-dimensional specialization within the framework of a three-dimensional theory, which is different the formulation based on a fundamentally two-dimensional theory.     

The scattering of harmonic elastic anti-plane shear waves by two collinear cracks in anisotropic material plane by using the non-local theory

In this paper, the dynamic behavior of two collinear cracks in the anisotropic elasticity material plane subjected to the harmonic anti-plane shear waves is investigated by use of the nonlocal theory. To overcome the mathematical difficulties, a one-dimensional nonlocal kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress field near the crack tips. By use of the Fourier transform, the problem can be solved with the help of a pair of triple integral equations, in which the unknown variable is the displacement on the crack surfaces. To solve the triple integral equations, the displacement on the crack surfaces is expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present near crack tips. The nonlocal elasticity solutions yield a finite hoop stress at the crack tips, thus allowing us to using the maximum stress as a fracture criterion. The magnitude of the finite stress field not only depends on the crack length but also on the frequency of the incident waves and the lattice parameter of the materials.     

A homogenized microstructural approach to creep fracture

A two-dimensional nonlocal continuum model is proposed in this paper for creep damage in polycrystalline materials. Starting from previous micromechanical modeling, a heuristic homogenization approach is adopted to derive a theory for the macroscopic response. The model accounts for the main damage mechanisms (grain boundary sliding, nucleation, growth and coalescence of cavities along the grain boundaries) responsible for the creep fracture process. The resulting constitutive law takes into account the nonlocalities expressed through the gradients of the stresses and the damage variables.     

Analysis of wave propagation in micro/nanobeam-like structures: A size-dependent model

By incorporating the strain gradient elasticity into the classical Bernoulli-Euler beam and Timoshenko beam models, the size-dependent characteristics of wave propagation in micro/nanobeams is studied. The formulations of dispersion relation are explicitly derived for both strain gradient beam models, and presented for different material length scale parameters (MLSPs). For both phenomenological sizedependent beam models, the angular frequency, phase velocity and group velocity increase with increasing wave number. However, the velocity ratios approach different values for different beam models, indicating an interesting behavior of the asymptotic velocity ratio. The present theory is also compared with the nonlocal continuum beam models.     

Nonlocal theory solution of two collinear cracks in the functionally graded materials

In this paper, the interaction of two collinear cracks in functionally graded materials subjected to a uniform anti-plane shear loading is investigated by means of nonlocal theory. The traditional concepts of the nonlocal theory are extended to solve the fracture problem of functionally graded materials. To make the analysis tractable, it is assumed that the shear modulus varies exponentially with the coordinate vertical to the crack. By use of the Fourier transform, the problem can be solved with the help of a pair of triple integral equations, in which the unknown variable is the displacement on the crack surfaces. To solve the triple integral equations, the displacement on the crack surfaces is expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present near the crack tips. The nonlocal elastic solutions yield a finite hoop stress at the crack tip, thus allowing us to use the maximum stress as a fracture criterion in functionally graded materials. The magnitude of the finite stress field depends on the crack length, the distance between two cracks, the parameter describing the functionally graded materials and the lattice parameter of the materials.     

A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations

This study develops a gradient theory of single-crystal plasticity that accounts for geometrically necessary dislocations. The theory is based on classical crystalline kinematics; classical macroforces; microforces for each slip system consistent with a microforce balance; a mechanical version of the second law that includes, via the microforces, work performed during slip; a rate-independent constitutive theory that includes dependences on a tensorial measure of geometrically necessary dislocations. The microforce balances are shown to be equivalent to nonlocal yield conditions for the individual slip systems. The field equations consist of the yield conditions coupled to the standard macroscopic force balance; these are supplemented by classical macroscopic boundary conditions in conjunction with nonstandard boundary conditions associated with slip. As an aid to solution, a weak (virtual power) formulation of the nonlocal yield conditions is derived. To make contact with classical dislocation theory, the microstresses are shown to represent counterparts of the Peach-Koehler force on a single dislocation.     

Dynamic meshless method applied to nonlocal crack problems

The dynamic meshless methods for local and nonlocal field theories are formulated in this paper. Application to two crack problems is presented. The meshless method of local theory gives solution that is in good agreement with the classical analytical crack tip solution, while the nonlocal theory yields a solution without stress singularity at the crack tip. The numerical results also show the embedded nonlocal nature of meshless methods.