基于广义应变梯度理论的纳米梁挠曲电效应研究

挠曲电效应是应变梯度与电极化的耦合,它存在于所有的电介质材料中。在纳米电介质结构的挠曲电效应研究中,应变梯度弹性对挠曲电响应的影响一直以来被低估甚至被忽略了。根据广义应变梯度理论,应变梯度弹性中独立的尺度参数只有三个,而文献中所采用的一个或两个尺度参数的应变梯度理论只是它的简化形式。基于该理论,论文建立了考虑广义应变梯度弹性的三维电介质结构的理论模型,并以一维纳米梁为例研究了其弯曲问题的挠曲电响应及其能量俘获特性。结果表明,纳米梁的挠曲电响应存在尺寸效应,并且弹性应变梯度会影响结构挠曲电的尺寸效应,特别是当结构的特征尺寸低于尺度参数时。论文的工作为更进一步理解纳米尺度下的挠曲电机理和能量俘获特性提供理论基础和设计依据。     

应用CBR确定面心立方金属的表面弹性参量

首先从能量变分出发基于同时考虑应变梯度效应和表面效应的跨尺度力学理论, 推导出表面能和表面弹性本构等基本关系, 然后基于简单的准连续Cauchy-Born法则(CBR)建立一种确定表面能密度以及表面弹性参量的方法.进一步以面心立方(face-centre-cubic,FCC)金属为例, 系统地获得了常用FCC金属表面弹性参量的数值, 结果与他人应用分子动力学计算得到的结果相吻合.     

Experiments and theory in strain gradient elasticity

Conventional strain-based mechanics theory does not account for contributions from strain gradients. Failure to include strain gradient contributions can lead to underestimates of stresses and size-dependent behaviors in small-scale structures. In this paper, a new set of higher-order metrics is developed to characterize strain gradient behaviors. This set enables the application of the higher-order equilibrium conditions to strain gradient elasticity theory and reduces the number of independent elastic length scale parameters from five to three. On the basis of this new strain gradient theory, a strain gradient elastic bending theory for plane-strain beams is developed. Solutions for cantilever bending with a moment and line force applied at the free end are constructed based on the new higher-order bending theory. In classical bending theory, the normalized bending rigidity is independent of the length and thickness of the beam. In the solutions developed from the higher-order bending theory, the normalized higher-order bending rigidity has a new dependence on the thickness of the beam and on a higher-order bending parameter, b_h. To determine the significance of the size dependence, we fabricated micron-sized beams and conducted bending tests using a nanoindenter. We found that the normalized beam rigidity exhibited an inverse squared dependence on the beam's thickness as predicted by the strain gradient elastic bending theory, and that the higher-order bending parameter, b_h, is on the micron-scale. Potential errors from the experiments, model and fabrication were estimated and determined to be small relative to the observed increase in beam's bending rigidity. The present results indicate that the elastic strain gradient effect is significant in elastic deformation of small-scale structures.     

A micromechanics-based strain gradient damage model for fracture prediction of brittle materials – Part I: Homogenization methodology and constitutive relations

In this paper, we first describe a homogenization methodology with the aim of establishing strain gradient constitutive relations for heterogeneous materials. The methodology presented in this work includes two main steps. The first one is the construction of the average strain-energy density for a well-chosen RVE by using a homogenization technique. The second one is the transformation of the obtained average strain-energy density to that for the continuum. An important characteristic of this method is its self-consistency with respect to the choice of the RVE: the strain gradient constitutive law built by using the present method is independent of the size and the form of the RVE. In the frame of this homogenization procedure, we have constructed a strain gradient constitutive relation for a two-dimensional elastic material with many microcracks by adopting the self-consistent scheme. It was shown that the effective behavior of cracked solids depends not only on the crack density but also on the average crack size with which the strain gradient is associated. The proposed constitutive relation provides a starting point for the development of an evolution law of damage including strain gradient effect, which will be presented in the second part of this work.     

Static, stability and dynamic analysis of gradient elastic flexural Kirchhoff plates

The governing equation of motion of gradient elastic flexural Kirchhoff plates, including the effect of in-plane constant forces on bending, is explicitly derived. This is accomplished by appropriately combining the equations of flexural motion in terms of moments, shear and in-plane forces, the moment–stress relations and the stress–strain equations of a simple strain gradient elastic theory with just one constant (the internal length squared), in addition to the two classical elastic moduli. The resulting partial differential equation in terms of the lateral deflection of the plate is of the sixth order instead of the fourth, which is the case for the classical elastic case. Three boundary value problems dealing with static, stability and dynamic analysis of a rectangular simply supported all-around gradient elastic flexural plate are solved analytically. Non-classical boundary conditions, in additional to the classical ones, have to be utilized. An assessment of the effect of the gradient coefficient on the static or dynamic response of the plate, its buckling load and natural frequencies is also made by comparing the gradient type of solutions against the classical ones.     

固体电介质中的挠曲电效应

挠曲电效应通常描述为非均匀变形如应变梯度引起的电极化或者电场梯度引起的变形.应变梯度能够局部破坏晶体的反演对称从而在材料中诱导电极化,因此挠曲电效应是固体电介质材料中普遍存在的一种力电耦合效应.应变梯度和电场梯度均随材料尺寸的减小而迅速增大,在宏观尺度通常被忽略的挠曲电效应在微纳尺度反而起着非常重要的作用,会显著影响材...     

An asymptotic strain gradient Reissner-Mindlin plate model

In this paper we derive a strain gradient plate model from the three-dimensional equations of strain gradient linearized elasticity. The deduction is based on the asymptotic analysis with respect of a small real parameter being the thickness of the elastic body we consider. The body is constituted by a second gradient isotropic linearly elastic material. The obtained model is recognized as a strain gradient Reissner-Mindlin plate model. We also provide a mathematical justification of the obtained plate model by means of a variational weak convergence result.     

一种基于应变梯度监测Ⅰ型裂纹的方法

基于线弹性断裂力学中I型裂纹的欧文解答,解析推导了在单向拉伸作用下无限大平板中I型裂纹尖端应变梯度场,建立了应变梯度与裂纹扩展之间的关联;基于挠曲电效应建立了电极化强度与应变梯度之间的力电耦合关系,提出了一种利用应变梯度传感器监测I型裂纹的方法,获知裂纹尖端坐标和裂纹扩展长度.本研究拟为应用应变梯度传感器对工程结构中裂纹扩展的实时监测提供初步的理论依据及方法.挠曲电感应技术在结构健康监测领域前景广阔.     

Suppressed plastic deformation at blunt crack-tips due to strain gradient effects

Large deformation gradients occur near a crack-tip and strain gradient dependent crack-tip deformation and stress fields are expected. Nevertheless, for material length scales much smaller than the scale of the deformation gradients, a conventional elastic–plastic solution is obtained. On the other hand, for significant large material length scales, a conventional elastic solution is obtained. This transition in behaviour is investigated based on a finite strain version of the Fleck–Hutchinson strain gradient plasticity model from 2001. The predictions show that for a wide range of material parameters, the transition from the conventional elastic–plastic to the elastic solution occurs for length scales ranging from 0.001 times the size of the plastic zone to a length scale of the same order of magnitude as the plastic zone.     

Derivation of Mindlin’s first and second strain gradient elastic theory via simple lattice and continuum models

Mindlin, in his celebrated papers of Arch. Rat. Mech. Anal. 16, 51–78, 1964 and Int. J. Solids Struct. 1, 417–438, 1965, proposed two enhanced strain gradient elastic theories to describe linear elastic behavior of isotropic materials with micro-structural effects. Since then, many works dealing with strain gradient elastic theories, derived either from lattice models or homogenization approaches, have appeared in the literature. Although elegant, none of them reproduces entirely the equation of motion as well as the classical and non-classical boundary conditions appearing in Mindlin theory, in terms of the considered lattice or continuum unit cell. Furthermore, no lattice or continuum models that confirm the second gradient elastic theory of Mindlin have been reported in the literature. The present work demonstrates two simple one dimensional models that conclude to first and second strain gradient elastic theories being identical to the corresponding ones proposed by Mindlin. The first is based on the standard continualization of the equation of motion taken for a sequence of mass-spring lattices, while the second one exploits average processes valid in continuum mechanics. Furthermore, Mindlin developed his theory by adding new terms in the expressions of potential and kinetic energy and introducing intrinsic micro-structural parameter without however providing explicit expressions that correlate micro-structure with macro-structure. This is accomplished in the present work where in both models the derived internal length scale parameters are correlated to the size of the considered unit cell.     

On bending of strain gradient elastic micro-plates

Bending of strain gradient elastic thin plates is studied, adopting Kirchhoff’s theory of plates. Simple linear strain gradient elastic theory with surface energy is employed. The governing plate equation with its boundary conditions are derived through a variational method. It turns out that new terms are introduced, indicating the importance of the cross-section area in bending of thin plates. Those terms are missing from the existing strain gradient plate theories; however, they strongly increase the stiffness of the thin plate.     

Strain gradient solution for a finite-domain Eshelby-type anti-plane strain inclusion problem

A solution for the finite-domain Eshelby-type inclusion problem of a finite elastic body containing an anti-plane strain inclusion of arbitrary cross-sectional shape prescribed with a uniform eigenstrain and a uniform eigenstrain gradient is derived in a general form using a simplified strain gradient elasticity theory (SSGET). The formulation is facilitated by an extended Betti’s reciprocal theorem and an extended Somigliana’s identity based on the SSGET and suitable for anti-plane strain problems. The disturbed displacement field is obtained in terms of the SSGET-based Green’s function for an infinite anti-plane strain elastic body. The solution reduces to that of the infinite-domain anti-plane strain inclusion problem when the boundary effect is not considered. The problem of a circular cylindrical inclusion embedded concentrically in a finite cylindrical elastic matrix undergoing anti-plane strain deformations is analytically solved by applying the general solution, with the Eshelby tensor and its average over the circular cross section of the inclusion obtained in closed forms. This Eshelby tensor, being dependent on the position, inclusion size, matrix size, and a material length scale parameter, captures the inclusion size and boundary effects, unlike existing ones. It reduces to the classical linear elasticity-based Eshelby tensor for the circular cylindrical inclusion in an infinite matrix if both the strain gradient and boundary effects are suppressed. Numerical results quantitatively show that the inclusion size effect can be quite large when the inclusion is small and that the boundary effect can dominate when the inclusion volume fraction is high. However, the inclusion size effect is diminishing with the increase of the inclusion size, and the boundary effect is vanishing as the inclusion volume fraction becomes sufficiently low.     

Bending and buckling of thin strain gradient elastic beams

Bending of strain gradient elastic thin beams is studied adopting Bernoulli-Euler principle. Simple linear strain gradient elastic theory with surface energy is employed. The governing beam equations with its boundary conditions are derived through a variational method. It turns out that new terms are introduced, indicating the importance of the cross-section area in bending of thin beams. Those terms are missing from the existing strain gradient beam theories. Those terms increase highly the stiffness of the thin beam. The buckling problem of the thin beams is also discussed.     

Wedge and twist disclinations in second strain gradient elasticity

The aim of this paper is to study disclinations in the framework of a second strain gradient elasticity theory. This second strain gradient elasticity has been proposed based on the first and second gradients of the strain tensor by Lazar et al. [Lazar, M., Maugin, G.A., Aifantis, E.C., 2006. Dislocations in second strain gradient elasticity. Int. J. Solids Struct. 43, 1787–1817]. Such a theory is an extension of the first strain gradient elasticity [Lazar, M., Maugin, G.A., 2005. Nonsingular stress and strain fields of dislocations and disclinations in first strain gradient elasticity. Int. J. Eng. Sci. 43, 1157–1184] with triple stress. By means of the stress function method, the exact analytical solutions for stress and strain fields of straight disclinations in an infinitely extended linear isotropic medium have been found. An important result is that the force stress, double stress and triple stress produced by wedge and twist disclinations are nonsingular. Meanwhile, the corresponding elastic strain and its gradients are also nonsingular. Analytical results indicate that the second strain gradient theory has the capacity of eliminating all unphysical singularities of physical fields.     

Microstructural effects in non linear stratified composites

In this paper, we analyze the microstructural effects on non linear elastic and periodic composites within the framework of asymptotic homogenization. We assume that the constitutive laws of the individual constituents derive from strain potentials. The microstructural effects are incorporated by considering the higher order terms, which come from the asymptotic series expansion. The complete solution at any order requires the resolution of a chain of cell problems in which the source terms depend on the solution at the lower order. The influence of these terms on the macroscopic response of the non linear composite is evaluated in the particular case of a stratified microstructure. The analytic solutions of the cell problems at the first and second order are provided for arbitrary local strain–stress laws which derive from potentials. As classically, the non-linear dependence on the applied macroscopic strain is retrieved for the solution at the first order. It is proved that the second order term in the expansion series also exhibits a non linear dependence with the macroscopic strain but linearly depends on the gradient of macroscopic strain. As a consequence, the macroscopic potential obtained by homogenization is a quadratic function of the macroscopic strain gradient when the expansion series is truncated at the second order. This model generalizes the well known first strain gradient elasticity theory to the case of non linear elastic material. The influence of the non local correctors on the macroscopic potential is investigated in the case of power law elasticity under macroscopic plane strain or antiplane conditions.     

基于应变率阶跃测试的晶体铜压入应变率敏感性研究

纳米压入测试可以原位获取材料的诸多力学性能,包括弹性模量,硬度,屈服应力,应变率敏感指数等。本文利用应变率阶跃测试技术对多晶铜试样的应变率敏感性进行测试分析,硬度-位移曲线表明压头下方所存在的变形梯度对各阶跃应变率下的硬度值存在明显影响;采用基于晶体细观机制的塑性应变梯度理论对压入变形梯度效应予以修正,比较了修正与未修正数据所得的应变率敏感指数,在有效剔除压入变形梯度影响的基础上,应变率阶跃测试可实现单次压入下材料应变率敏感性的测试表征。     

Nonlinear strain gradient elastic thin shallow shells

The governing equilibrium equations for strain gradient elastic thin shallow shells are derived, considering nonlinear strains and linear constitutive strain gradient elastic relations. Adopting Kirchhoff’s theory of thin shallow structures, the equilibrium equations, along with the boundary conditions, are formulated through a variational procedure. It turns out that new terms are introduced, indicating the importance of the cross-section area in bending of thin plates. Those terms are missing from the existing strain gradient shallow thin shell theories. Those terms highly increase the stiffness of the structures. When the curvature of the shallow shell becomes zero, the governing equilibrium for the plates is derived.     

纳米圆轴在弹性介质中的扭转振动分析

基于非局部应变梯度理论,考虑周围弹性介质的影响,研究纳米圆轴的扭转自由振动。首先通过Hamilton原理推导纳米圆轴扭转振动的控制方程及边界条件,接着采用微分求积法得到控制方程及三类边界条件的离散形式,最后由数值计算结果分析扭转振动特性。讨论了两个小尺度参数和弹性介质刚度的变化对振动频率的影响,并通过小尺度参数比对振动频率的影响分析两个尺度参数的耦合作用。研究结果表明,扭转自由振动频率随非局部参数增加而减小,随应变梯度尺度参数、弹性介质刚度增加而增大;当非局部参数大于应变梯度尺度参数时,小尺度效应体现为非局部效应,相反则体现为应变梯度效应。     

A framework for numerical integration of crystal elasto-plastic constitutive equations compatible with explicit finite element codes

In this paper we summarize the elements of a numerical integration scheme for elasto-plastic response of single crystals. This is intended to be compatible with large-scale explicit finite element codes and therefore can be used for problems involving multiple crystals and also overall behavior of polycrystalline materials. The steps described here are general for anisotropic elastic and plastic response of crystals. The crystallographic axes of the lattice are explicitly stored and updated at each time step. A plastic predictor–elastic corrector scheme is used to calculate the plastic strain rates on all active slip systems based on a rate-dependent physics-based constitutive model without the need of further auxiliary assumptions. Finally we present the results of numerous calculations using a physics-based rate- and temperature-dependent model of copper and the effect of elastic unloading, elastic crystal anisotropy, and deformation-induced lattice rotation are emphasized.     

Strain gradient solution for a finite-domain Eshelby-type plane strain inclusion problem and Eshelby’s tensor for a cylindrical inclusion in a finite elastic matrix

A solution for the finite-domain Eshelby-type inclusion problem of a finite elastic body containing a plane strain inclusion prescribed with a uniform eigenstrain and a uniform eigenstrain gradient is derived in a general form using a simplified strain gradient elasticity theory (SSGET). The formulation is facilitated by an extended Betti’s reciprocal theorem and an extended Somigliana’s identity based on the SSGET and suitable for plane strain problems. The disturbed displacement field is obtained in terms of the SSGET-based Green’s function for an infinite plane strain elastic body, which differs from that in earlier studies using the three-dimensional Green’s function. The solution reduces to that of the infinite-domain inclusion problem when the boundary effect is suppressed. The problem of a cylindrical inclusion embedded concentrically in a finite plane strain cylindrical elastic matrix of an enhanced continuum is analytically solved for the first time by applying the general solution, with the Eshelby tensor and its average over the circular cross section of the inclusion obtained in closed forms. This Eshelby tensor, being dependent on the position, inclusion size, matrix size, and a material length scale parameter, captures the inclusion size and boundary effects, unlike existing ones. It reduces to the classical elasticity-based Eshelby tensor for the cylindrical inclusion in an infinite matrix if both the strain gradient and boundary effects are not considered. Numerical results quantitatively show that the inclusion size effect can be quite large when the inclusion is very small and that the boundary effect can dominate when the inclusion volume fraction is very high. However, the inclusion size effect is diminishing with the increase of the inclusion size, and the boundary effect is vanishing as the inclusion volume fraction becomes sufficiently low.