Size-dependent electromechanical properties in piezoelectric superlattices due to flexoelectric effect

Piezoelectric superlattice is a potential component for nanoelectromechanical systems.Due to the strong nonlocal effect such as flexoelectric effect at interfaces,classical piezoelectric theory is unable to accurately describe the electromechanical response of piezoelectric superlattice at nanoscale scale.Based on the previous nonlocal thermodynamics theory with flexoelectric effect Liu et al.(2016),the sizedependent electromechanical properties of piezoelectric superlattices made of BaTiO_3(BTO)and PbTiO_3(PTO)layers are investigated systematically in the present work.Giant strain gradient is found near the interface between BTO and PTO layers,which leads to the significant enhancement of polarization in the superlattice due to the flexoelectric effect.For the piezoelectric BTO–PTO superlattices with different unitcell sizes,the thickness of interface with nontrivial strain gradient is almost constant.The influence of strain gradient at the interface becomes significant when the size of superlattice decreases.As a result,a strong size dependence of electromechanical properties is predicted for the piezoelectric BTO–PTO superlattices.In particular,for the superlattices with a specific thickness ratio of BTO and PTO layers,the piezoelectric response can be several times larger than that of bulk structure.The present work demonstrates a practical way to design the piezoelectric superlattices with high piezoelectric coefficient by using the nonlocal effect at nanoscale.     

Limit analysis of viscoplastic thick-walled cylinder and spherical shell under internal pressure using a strain gradient plasticity theory

Plastic limit load of viscoplastic thick-walled cylinder and spherical shell subjected to internal pressure is investigated analytically using a strain gradient plasticity theory. As a result, the current solutions can capture the size effect at the micron scale. Numerical results show that the smaller the inner radius of the cylinder or spherical shell, the more significant the scale effects. Results also show that the size effect is more evident with increasing strain or strain-rate sensitivity index. The classical plastic-based solutions of the same problems are shown to be a special case of the present solution.     

Bauschinger and size effects in thin-film plasticity due to defect-energy of geometrical necessary dislocations

The Bauschinger and size effects in the thinfilm plasticity theory arising from the defect-energy of geometrically necessary dislocations (GNDs) are analytically investigated in this paper. Firstly, this defect-energy is deduced based on the elastic interactions of coupling dislocations (or pile-ups) moving on the closed neighboring slip plane. This energy is a quadratic function of the GNDs density, and includes an elastic interaction coefficient and an energetic length scale L. By incorporating it into the work- conjugate strain gradient plasticity theory of Gurtin, an energetic stress associated with this defect energy is obtained, which just plays the role of back stress in the kinematic hardening model. Then this back-stress hardening model is used to investigate the Bauschinger and size effects in the tension problem of single crystal Al films with passivation layers. The tension stress in the film shows a reverse dependence on the film thickness h. By comparing it with discrete-dislocation simulation results, the length scale L is determined, which is just several slip plane spacing, and accords well with our physical interpretation for the defect- energy. The Bauschinger effect after unloading is analyzed by combining this back-stress hardening model with a friction model. The effects of film thickness and pre-strain on the reversed plastic strain after unloading are quantified and qualitatively compared with experiment results.     

On the homogenization of metal matrix composites using strain gradient plasticity

The homogenized response of metal matrix composites（MMC） is studied using strain gradient plasticity.The material model employed is a rate independent formulation of energetic strain gradient plasticity at the micro scale and conventional rate independent plasticity at the macro scale. Free energy inside the micro structure is included due to the elastic strains and plastic strain gradients. A unit cell containing a circular elastic fiber is analyzed under macroscopic simple shear in addition to transverse and longitudinal loading. The analyses are carried out under generalized plane strain condition. Micro-macro homogenization is performed observing the Hill-Mandel energy condition,and overall loading is considered such that the homogenized higher order terms vanish. The results highlight the intrinsic size-effects as well as the effect of fiber volume fraction on the overall response curves, plastic strain distributions and homogenized yield surfaces under different loading conditions. It is concluded that composites with smaller reinforcement size have larger initial yield surfaces and furthermore,they exhibit more kinematic hardening.     

Generalizing J 2 flow theory: Fundamental issues in strain gradient plasticity

It has not been a simple matter to obtain a sound extension of the classical J₂ flow theory of plasticity that in- corporates a dependence on plastic strain gradients and that is capable of capturing size-dependent behaviour of metals at the micron scale. Two classes of basic extensions of clas- sical J₂ theory have been proposed: one with increments in higher order stresses related to increments of strain gradi- ents and the other characterized by the higher order stresses themselves expressed in terms of increments of strain gra- dients. The theories proposed by Muhlhaus and Aifantis in 1991 and Fleck and Hutchinson in 2001 are in the first class, and, as formulated, these do not always satisfy ther- modynamic requirements on plastic dissipation. On the other hand, theories of the second class proposed by Gudmundson in 2004 and Gurtin and Anand in 2009 have the physical deficiency that the higher order stress quantities can change discontinuously for bodies subject to arbitrarily small load changes. The present paper lays out this background to the quest for a sound phenomenological extension of the rate- independent J₂ flow theory of plasticity to include a de- pendence on gradients of plastic strain. A modification of the Fleck-Hutchinson formulation that ensures its thermo- dynamic integrity is presented and contrasted with a compa- rable formulation of the second class where in the higher or- der stresses are expressed in terms of the plastic strain rate. Both versions are constructed to reduce to the classical J₂ flow theory of plasticity when the gradients can be neglected and to coincide with the simpler and more readily formulated J₂ deformation theory of gradient plasticity for deformation histories characterized by proportional straining.     

SIMULATIONS OF MECHANICAL BEHAVIOR OF POLYCRYSTALLINE COPPER WITH NANO-TWINS

Mechanical behavior and microstructure evolution of polycrystalline copper with nano-twins were investigated in the present work by finite element simulations. The fracture of grain boundaries are described by a cohesive interface constitutive model based on the strain gradient plasticity theory. A systematic study of the strength and ductility for different grain sizes and twin lamellae distributions is performed. The results show that the material strength and ductility strongly depend on the grain size and the distribution of twin lamellae microstructures in the polycrystalline copper.     

A unified plasticity methodology for rate-and temperature-sensitive alloys exhibiting a non-linear kinematic hardening behavior

In this paper, a novel unified plasticity methodology is proposed to allow the coupling of rate-and temperature-sensitivity of engineering alloys as well as the non-linear kinematic hardening behavior often observed during cyclic loading. The proposed methodology is general in the sense that an arbitrary constitutive model may be chosen for the viscoplastic part, as well as the cyclic part. We adapt our model with a physically-motivated viscoplasticity flow rule and a nonlinear kinematic hardening model. In contrast with other unified plasticity models, the simplified theory involves few material parameters that can be readily calibrated from standard mechanical tests. The capabilities of the proposed theory are demonstrated for a hot rolled annealed 304 L stainless steel supplied by Vimetal Peckover. The model is tested with stress–strain curves obtained from standard tensile and cyclic uniaxial tests at various strain amplitudes and strain-rates, and good accuracy of the response is obtained for strains up to 15%, within a temperature range of 293–673 K. We note that the cyclic plasticity model in our adapted theory can be readily enhanced with ratchetting, mean stress relaxation, strain amplitude history, Masing effects or other complex capabilities.     

Spinodal surface instability of soft elastic thin films

When the thicknesses of thin films reduce to microns or even nanometers, surface energy and surface interaction often play a significant role in their deformation behavior and surface morphology. The spinodal surface instability induced by the van der Waals force in a soft elastic thin film perfectly bonded to a rigid substrate is investigated theoretically using the bifurcation theory of elastic structures. The analytical solution is derived for the critical condition of spinodal surface morphology instability by accounting for the competition of the van der Waals interaction energy, elastic strain energy and surface energy. Detailed examinations on the effect of surface energy, thickness and elastic properties of the film show that the characteristic wavelength of the deformation bifurcation mode depends on the film thickness via an exponential relation, with the power index in the range from 0.749 to 1.0. The theoretical solution has a good agreement with relevant experiment results.     

Hierarchical structure observation and nanoindentation size effect characterization for a limnetic shell

In the present research, hierarchical structure observation and mechanical property characterization for a type of biomaterial are carried out. The investigated biomaterial is Hyriopsis cumingii, a typical limnetic shell, which consists of two different structural layers, a prismatic "pillar"structure and a nacreous "brick and mortar" structure. The prismatic layer looks like a "pillar forest" with variationsection pillars sized on the order of several tens of microns.The nacreous material looks like a "brick wall" with bricks sized on the order of several microns. Both pillars and bricks are composed of nanoparticles. The mechanical properties of the hierarchical biomaterial are measured by using the nanoindentation test. Hardness and modulus are measured for both the nacre layer and the prismatic layer, respectively.The nanoindentation size effects for the hierarchical structural materials are investigated experimentally. The results show that the prismatic nanostructured material has a higher stiffness and hardness than the nacre nanostructured material.In addition, the nanoindentation size effects for the hierarchical structural materials are described theoretically, by using the trans-scale mechanics theory considering both strain gradient effect and the surface/interface effect. The modeling results are consistent with experimental ones.     

Finite versus small strain discrete dislocation analysis of cantilever bending of single crystals

Plastic size effects in single crystals are investi-gated by using finite strain and small strain discrete dislo-cation plasticity to analyse the response of cantilever beam specimens. Crystals with both one and two active slip sys-tems are analysed, as well as specimens with different beam aspect ratios. Over the range of specimen sizes analysed here, the bending stress versus applied tip displacement response has a strong hardening plastic component. This hardening rate increases with decreasing specimen size. The hardening rates are slightly lower when the finite strain discrete disloca-tion plasticity (DDP) formulation is employed as curving of the slip planes is accounted for in the finite strain formulation. This relaxes the back-stresses in the dislocation pile-ups and thereby reduces the hardening rate. Our calculations show that in line with the pure bending case, the bending stress in cantilever bending displays a plastic size dependence. How-ever, unlike pure bending, the bending flow strength of the larger aspect ratio cantilever beams is appreciably smaller. This is attributed to the fact that for the same applied bend-ing stress, longer beams have lower shear forces acting upon them and this results in a lower density of statistically stored dislocations.     

Wedge indentation of thin films modelled by strain gradient plasticity

A plane strain study of wedge indentation of a thin film on a substrate is performed. The film is modelled with the strain gradient plasticity theory by Gudmundson [Gudmundson, P., 2004. A unified treatment of strain gradient plasticity. Journal of the Mechanics and Physics of Solids 52, 1379–1406] and analysed using finite element simulations. Several trends that have been experimentally observed elsewhere are captured in the predictions of the mechanical behaviour of the thin film. Such trends include increased hardness at shallow depths due to gradient effects as well as increased hardness at larger depths due to the influence of the substrate. In between, a plateau is found which is observed to scale linearly with the material length scale parameter. It is shown that the degree of hardening of the material has a strong influence on the substrate effect, where a high hardening modulus gives a larger impact on this effect. Furthermore, pile-up deformation dominated by plasticity at small values of the internal length scale parameter is turned into sink-in deformation where plasticity is suppressed for larger values of the length scale parameter. Finally, it is demonstrated that the effect of substrate compliance has a significant effect on the hardness predictions if the effective stiffness of the substrate is of the same order as the stiffness of the film.     

Indentation of a hard film on a soft substrate: Strain gradient hardening effects

The strain gradient work hardening is important in micro-indentation of bulk metals and thin metallic films, though the indentation of thin films may display very different behavior from that of bulk metals. We use the conventional theory of mechanism-based strain gradient plasticity (CMSG) to study the indentation of a hard tungsten film on soft aluminum substrate, and find good agreement with experiments. The effect of friction stress (intrinsic lattice resistance), which is important in body-center-cubic tungsten, is accounted for. We also extend CMSG to a finite deformation theory since the indentation depth in experiments can be as large as the film thickness. Contrary to indentation of bulk metals or soft metallic films on hard substrate, the micro-indentation hardness of a hard tungsten film on soft aluminum substrate decreases monotonically with the increasing depth of indentation, and it never approaches a constant (macroscopic hardness). It is also shown that the strain gradient effect in the soft aluminum substrate is insignificant, but that in the hard tungsten thin film is important in shallow indentation. The strain gradient effect in tungsten, however, disappears rapidly as the indentation depth increases because the intrinsic material length in tungsten is rather small.     

Size effects in the conical indentation of an elasto-plastic solid

The size effect in conical indentation of an elasto-plastic solid is predicted via the Fleck and Willis formulation of strain gradient plasticity (Fleck, N.A. and Willis, J.R., 2009, A mathematical basis for strain gradient plasticity theory. Part II: tensorial plastic multiplier, J. Mech. Phys. Solids, 57, 1045–1057). The rate-dependent formulation is implemented numerically and the full-field indentation problem is analyzed via finite element calculations, for both ideally plastic behavior and dissipative hardening. The isotropic strain-gradient theory involves three material length scales, and the relative significance of these length scales upon the degree of size effect is assessed. Indentation maps are generated to summarize the sensitivity of indentation hardness to indent size, indenter geometry and material properties (such as yield strain and strain hardening index). The finite element model is also used to evaluate the pertinence of the Johnson cavity expansion model and of the Nix–Gao model, which have been extensively used to predict size effects in indentation hardness.     

Small scale,grain size and substrate effects in nano-indentation experiment of film–substrate systems

The mechanical properties of film–substrate systems have been investigated through nano-indentation experiments in our former paper (Chen, S.H., Liu, L., Wang, T.C., 2005. Investigation of the mechanical properties of thin films by nano-indentation, considering the effects of thickness and different coating–substrate combinations. Surf. Coat. Technol., 191, 25–32), in which Al–Glass with three different film thicknesses are adopted and it is found that the relation between the hardness H and normalized indentation depth h/t, where t denotes the film thickness, exhibits three different regimes: (i) the hardness decreases obviously with increasing indentation depth; (ii) then, the hardness keeps an almost constant value in the range of 0.1–0.7 of the normalized indentation depth h/t; (iii) after that, the hardness increases with increasing indentation depth. In this paper, the indentation image is further investigated and finite element method is used to analyze the nano-indentation phenomena with both classical plasticity and strain gradient plasticity theories. Not only the case with an ideal sharp indenter tip but also that with a round one is considered in both theories. Finally, we find that the classical plasticity theory can not predict the experimental results, even considering the indenter tip curvature. However, the strain gradient plasticity theory can describe the experimental data very well not only at a shallow indentation depth but also at a deep depth. Strain gradient and substrate effects are proved to coexist in film–substrate nano-indentation experiments.     

The finite deformation theory of Taylor-based nonlocal plasticity

Recent experiments have shown that metallic materials display significant size effect at the micron and sub-micron scales. This has motivated the development of strain gradient plasticity theories, which usually involve extra boundary conditions and possibly higher-order governing equations. We propose a finite deformation theory of nonlocal plasticity based on the Taylor dislocation model. The theory falls into Rice's theoretical framework of internal variables [J Mech Phys Solids 19 (1971) 433], and it does not require any extra boundary conditions. We apply the theory to study the micro-indentation hardness experiments, and it agrees very well with the experimental data over a wide range of indentation depth.     

Modeling the interfacial effect on the yield strength and flow stress of thin metal films on substrates

It is shown in this paper that interfacial effects have a profound impact on the scale-dependent yield strength and strain hardening rates (flow stress) of metallic thin films on elastic substrates. This is achieved by developing a higher-order strain gradient plasticity theory based on the principle of virtual power and the laws of thermodynamics. This theory enforces microscopic boundary conditions at interfaces which relate a microtraction stress to the interfacial energy at the interface. It is shown that the film bulk length scale controls the size effect if a rigid interface is assumed whereas the interfacial length scale dominates if a compliant interface is assumed.     

Particulate size effects in the particle-reinforced metal-matrix composites

The influences of particle size on the mechanical properties of the particulate metal matrix composite are obviously displayed in the experimental observations. However, the phenomenon can not be predicted directly using the conventional elastic-plastic theory. It is because that no length scale parameters are involved in the conventional theory. In the present research, using the strain gradient plasticity theory, a systematic research of the particle size effect in the particulate metal matrix composite is carried out. The roles of many composite factors, such as: the particle size, the Young's modulus of the particle, the particle aspect ratio and volume fraction, as well as the plastic strain hardening exponent of the matrix material, are studied in detail. In order to obtain a general understanding for the composite behavior, two kinds of particle shapes, ellipsoid and cylinder, are considered to check the strength dependence of the smooth or non-smooth particle surface. Finally, the prediction results will be applied to the several experiments about the ceramic particle-reinforced metal-matrix composites. The material length scale parameter is predicted. The project supported by the National Natural Science Foundation of China (19891180, 19925211) and by the Chinese Academy of Sciences (KJ951-1-201) and “Bai Ren” plan     

An expanding cavity model incorporating strain-hardening and indentation size effects

An expanding cavity model (ECM) for determining indentation hardness of elastic strain-hardening plastic materials is developed. The derivation is based on a strain gradient plasticity solution for an internally pressurized thick-walled spherical shell of an elastic power-law hardening material. Closed-form formulas are provided for both conical and spherical indentations. The indentation radius enters these formulas with its own dimensional identity, unlike that in classical plasticity based ECMs where indentation geometrical parameters appear only in non-dimensional forms. As a result, the newly developed ECM can capture the indentation size effect. The formulas explicitly show that indentation hardness depends on Young’s modulus, yield stress, strain-hardening exponent and strain gradient coefficient of the indented material as well as on the geometry of the indenter. The new model reduces to existing classical plasticity based ECMs (including Johnson’s ECM for elastic–perfectly plastic materials) when the strain gradient effect is not considered. The numerical results obtained using the newly developed model reveal that the hardness is indeed indentation size dependent when the indentation radius is very small: the smaller the indentation, the larger the hardness. Also, the indentation hardness is seen to increase with the Young’s modulus and strain-hardening level of the indented material for both conical and spherical indentations. The strain-hardening effect on the hardness is observed to be significant for materials having strong strain-hardening characteristics. In addition, it is found that the indentation hardness increases with decreasing cone angle of the conical indenter or decreasing radius of the spherical indenter. These trends agree with existing experimental observations and model predictions.     

颗粒增强复合材料的界面开裂与尺度效应

采用基于Huang等提出的塑性应变梯度传统理论发展的有限元方法,模拟了颗粒增强金属基复合材料的界面开裂与颗粒尺度效应.分别针对考虑颗粒与基体间界面开裂和不开裂两种情况进行分析,并将考虑界面开裂的模拟结果与实验结果进行比较,证明了模型的有效性,同时也获得应变梯度理论中所包含的材料特征尺度参量的取值.     

A conventional theory of mechanism-based strain gradient plasticity

There exist two frameworks of strain gradient plasticity theories to model size effects observed at the micron and sub-micron scales in experiments. The first framework involves the higher-order stress and therefore requires extra boundary conditions, such as the theory of mechanism-based strain gradient (MSG) plasticity [J Mech Phys Solids 47 (1999) 1239; J Mech Phys Solids 48 (2000) 99; J Mater Res 15 (2000) 1786] established from the Taylor dislocation model. The other framework does not involve the higher-order stress, and the strain gradient effect come into play via the incremental plastic moduli. A conventional theory of mechanism-based strain gradient plasticity is established in this paper. It is also based on the Taylor dislocation model, but it does not involve the higher-order stress and therefore falls into the second strain gradient plasticity framework that preserves the structure of conventional plasticity theories. The plastic strain gradient appears only in the constitutive model, and the equilibrium equations and boundary conditions are the same as the conventional continuum theories. It is shown that the difference between this theory and the higher-order MSG plasticity theory based on the same dislocation model is only significant within a thin boundary layer of the solid.