基于多尺度模型解释固体中的挠曲电效应

挠曲电效应是一种跨尺度的多场耦合现象。当前的宏观挠曲电理论均是基于应变梯度局部破坏晶体反演对称这一微观机理对该现象进行唯象描述。该宏观理论与基于晶格动力学及密度泛函理论的微观挠曲电理论模型之间存在较大差异。难以将两者结合用以跨尺度地研究材料中的挠曲电效应。针对该现状,本文基于前人提出的原子场理论,建立了一种新的多尺度挠曲电模型。并在该多尺度模型框架下解释了应变梯度诱发极化的微观机理。一方面,与基于连续介质力学的唯象理论不同,本文从材料微结构演化的角度推导了原子位移与极化的关系。另一方面,与通过晶格波假设原子位移的微观理论不同,本文得到的极化表达式更加真实和广义地解释了挠曲电效应。其能够适用于材料边界存在机械力作用,材料内部存在缺陷等复杂的情况。本文所建立的多尺度挠曲电模型能够为后续多尺度挠曲电效应的研究提供一些思路。     

A higher-order nonlocal elasticity and strain gradient theory and its applications in wave propagation

In recent years there have been many papers that considered the effects of material length scales in the study of mechanics of solids at micro- and/or nano-scales. There are a number of approaches and, among them, one set of papers deals with Eringen's differential nonlocal model and another deals with the strain gradient theories. The modified couple stress theory, which also accounts for a material length scale, is a form of a strain gradient theory. The large body of literature that has come into existence in the last several years has created significant confusion among researchers about the length scales that these various theories contain. The present paper has the objective of establishing the fact that the length scales present in nonlocal elasticity and strain gradient theory describe two entirely different physical characteristics of materials and structures at nanoscale. By using two principle kernel functions, the paper further presents a theory with application examples which relates the classical nonlocal elasticity and strain gradient theory and it results in a higher-order nonlocal strain gradient theory. In this theory, a higher-order nonlocal strain gradient elasticity system which considers higher-order stress gradients and strain gradient nonlocality is proposed. It is based on the nonlocal effects of the strain field and first gradient strain field. This theory intends to generalize the classical nonlocal elasticity theory by introducing a higher-order strain tensor with nonlocality into the stored energy function. The theory is distinctive because the classical nonlocal stress theory does not include nonlocality of higher-order stresses while the common strain gradient theory only considers local higher-order strain gradients without nonlocal effects in a global sense. By establishing the constitutive relation within the thermodynamic framework, the governing equations of equilibrium and all boundary conditions are derived via the variational approach. Two additional kinds of parameters, the higher-order nonlocal parameters and the nonlocal gradient length coefficients are introduced to account for the size-dependent characteristics of nonlocal gradient materials at nanoscale. To illustrate its application values, the theory is applied for wave propagation in a nonlocal strain gradient system and the new dispersion relations derived are presented through examples for wave propagating in Euler–Bernoulli and Timoshenko nanobeams. The numerical results based on the new nonlocal strain gradient theory reveal some new findings with respect to lattice dynamics and wave propagation experiment that could not be matched by both the classical nonlocal stress model and the contemporary strain gradient theory. Thus, this higher-order nonlocal strain gradient model provides an explanation to some observations in the classical and nonlocal stress theories as well as the strain gradient theory in these aspects.     

高维非线性系统的全局分岔和混沌动力学研究

综述了Melnikov方法的发展历史, 从1963年苏联学者Melnikov提出该方法开始, 一直到目前广义Melnikov方法的提出和发展. Melnikov方法的发展历程可以概括为3 个阶段, 分别综述了每一个阶段Melnikov方法的扩展和应用, 论述了国内外在该方向的研究现状和所获得的主要结果, 指出了各种Melnikov方法之间的联系、存在的问题和不足. 为了对比两种研究高维非线性系统多脉冲混沌动力学的理论, 本文综述了另外一种全局摄动理论, 即能量相位法, 总结了该方法十几年来的发展历史以及国内外的理论研究成果和工程应用实例, 阐述了能量相位法发展的根源以及与Melnikov方法的内在联系, 比较了能量相位法和广义Melnikov方法两种理论研究对象的差别, 以及各自所存在的不足和问题. 简要论述了能量相位法和广义Melnikov方法的理论体系, 并利用广义Melnikov方法分析了四边简支矩形薄板的多脉冲混沌动力学, 数值模拟进一步验证了理论研究的结果. 最后, 详细综述了两种理论的缺点和不足, 说明今后全局摄动理论的发展方向.     

Gradient elasticity in statics and dynamics: An overview of formulations,length scale identification procedures,finite element implementations and new results

In this paper, we discuss various formats of gradient elasticity and their performance in static and dynamic applications. Gradient elasticity theories provide extensions of the classical equations of elasticity with additional higher-order spatial derivatives of strains, stresses and/or accelerations. We focus on the versatile class of gradient elasticity theories whereby the higher-order terms are the Laplacian of the corresponding lower-order terms. One of the challenges of formulating gradient elasticity theories is to keep the number of additional constitutive parameters to a minimum. We start with discussing the general Mindlin theory, that in its most general form has 903 constitutive elastic parameters but which were reduced by Mindlin to three independent material length scales. Further simplifications are often possible. In particular, the Aifantis theory has only one additional parameter in statics and opens up a whole new field of analytical and numerical solution procedures. We also address how this can be extended to dynamics. An overview of length scale identification and quantification procedures is given. Finite element implementations of the most commonly used versions of gradient elasticity are discussed together with the variationally consistent boundary conditions. Details are provided for particular formats of gradient elasticity that can be implemented with simple, linear finite element shape functions. New numerical results show the removal of singularities in statics and dynamics, as well as the size-dependent mechanical response predicted by gradient elasticity.     

A large strain plasticity model for anisotropic materials — composite material application

In this work a generalized anisotropic model in large strains based on the classical isotropic plasticity theory is presented. The anisotropic theory is based on the concept of mapped tensors from the anisotropic real space to the isotropic fictitious one. In classical orthotropy theories it is necessary to use a special constitutive law for each material. The proposed theory is a generalization of classical theories and allows the use of models and algorithms developed for isotropic materials. It is based on establishing a one-to-one relationship between the behavior of an anisotropic real material and that of an isotropic fictitious one. Therefore, the problem is solved in the isotropic fictious space and the results are transported to the real field. This theory is applied to simulate the behavior of each material in the composite. The whole behavior of the composite is modeled by incorporating the anisotropic model within a model based on a modified mixing theory.     

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.     

Extended irreversible thermodynamics and its relation with other continuum approaches

We present an introduction to extended irreversible thermodynamics (EIT) as applied to polymer solutions in the presence of shear flow and of diffusion flux. We discuss with special attention the definition of chemical potential in non-equilibrium situations and its use in the analysis of shear-induced phase transitions. In the second part, we compare EIT with other contemporary continuum approaches: theories with internal variables, the GENERIC approach, and the matrix model. All these theories share an emphasis on the relations between dynamics and thermodynamics at a deeper level than in the classical theory, but each of them has some peculiar advantage in the analysis of some specific aspects of physical problems.     

About Clapeyron's Theorem in Linear Elasticity

We examine some elementary interpretations of the classical theorem of Clapeyron in linear elasticity theory. As we show, a straightforward application of this theorem in the purely mechanical setting leads to an apparent paradox which can be resolved by referring either to dynamics or to thermodynamics. These richer theories play an essential part in understanding the physical significance of this theorem.     

Deformation of a Peridynamic Bar

The deformation of an infinite bar subjected to a self-equilibrated load distribution is investigated using the peridynamic formulation of elasticity theory. The peridynamic theory differs from the classical theory and other nonlocal theories in that it does not involve spatial derivatives of the displacement field. The bar problem is formulated as a linear Fredholm integral equation and solved using Fourier transform methods. The solution is shown to exhibit, in general, features that are not found in the classical result. Among these are decaying oscillations in the displacement field and progressively weakening discontinuities that propagate outside of the loading region. These features, when present, are guaranteed to decay provided that the wave speeds are real. This leads to a one-dimensional version of St. Venant's principle for peridynamic materials that ensures the increasing smoothness of the displacement field remotely from the loading region. The peridynamic result converges to the classical result in the limit of short-range forces. An example gives the solution to the concentrated load problem, and hence provides the Green's function for general loading problems.     

Studying thin film damping in a micro-beam resonator based on non-classical theories

In this paper, a mathematical model is presented for studying thin film damping of the surrounding fluid in an in-plane oscillating micro-beam resonator. The proposed model for this study is made up of a clamped-clamped micro-beam bound between two fixed layers. The micro-gap between the micro-beam and fixed layers is filled with air. As classical theories are not properly capable of pre-dicting the size dependence behaviors of the micro-beam, and also behavior of micro-scale fluid media, hence in the presented model, equation of motion governing longitudinal displacement of the micro-beam has been extracted based on non-local elasticity theory. Furthermore, the fluid field has been modeled based on micro-polar theory. These coupled equations have been simplified using Newton-Laplace and continuity equations. After transforming to non-dimensional form and linearizing, the equations have been discretized and solved simultaneously using a Galerkin-based reduced order model. Considering slip boundary conditions and applying a complex frequency approach, the equivalent damping ratio and quality factor of the micro-beam resonator have been obtained. The obtained values for the quality factor have been compared to those based on classical theories. We have shown that applying non-classical theories underestimate the values of the quality factor obtained based on classical theo-ries. The effects of geometrical parameters of the micro-beam and micro-scale fluid field on the quality factor of the res-onator have also been investigated.     

Dynamical theory of hyperelastic rods

A rod is regarded as a one-dimensional mathematical model of a three-dimensional body. The exact field equations governing the motion of a hyperelastic rod are derived from the general three-dimensional theory. Then, by a suitable restriction on the number of displacement variables, a hierarchy of approximating theories is established. Because such theories are generated by a kinematic hypothesis, a precise, quantitative idea of the nature of the simplifying assumptions is furnished. An analysis of the structure of these approximating theories yields three distinct approaches by which they may be interpreted. Finally, constraints and their connection with other approximate theories are investigated. In particular, classical nonlinear theories and theories for planar motion are developed in this context of constrained theories.     

湍流被动标量研究的最新进展

被动标量场的统计性质, 在湍流理论以及湍流燃烧、污染物防治等工程领域都有非常重要的意义. 最近十几年, 大量的实验测量、数值模拟和理论分析结果表明, 标量场具有很多自身独特的性质, 甚至有些性质并不依赖于速度场本身. 这就促使人们对传统经典理论进行重新认识、修正或者提出新的理论来取代. 本文对标量场的各向异性、标量和标量耗散率的概率密度函数(PDF)、标量场的空间结构及演化过程、还有标量小尺度混合模型等几个方面进行综述.      

Buckling analysis of cylindrical thin-shells using strain gradient elasticity theory

The stability problem of cylindrical shells is addressed using higher-order continuum theories in a generalized framework. The length-scale effect which becomes prominent at microscale can be included in the continuum theory using gradient-based nonlocal theories such as the strain gradient elasticity theories. In this work, expressions for critical buckling stress under uniaxial compression are derived using an energy approach. The results are compared with the classical continuum theory, which can be obtained by setting the length-scale parameters to zero. A special case is obtained by setting two length scale parameters to zero. Thus, it is shown that both the couple stress theory and classical continuum theory forms a special case of the strain gradient theory. The effect of various parameters such as the shell-radius, shell-length, and length-scale parameters on the buckling stress are investigated. The dimensions and constants corresponding to that of a carbon nanotube, where the length-scale effect becomes prominent, is considered for this investigation.     

Singularity-free dislocation dynamics with strain gradient elasticity

The singular nature of the elastic fields produced by dislocations presents conceptual challenges and computational difficulties in the implementation of discrete dislocation-based models of plasticity. In the context of classical elasticity, attempts to regularize the elastic fields of discrete dislocations encounter intrinsic difficulties. On the other hand, in gradient elasticity, the issue of singularity can be removed at the outset and smooth elastic fields of dislocations are available. In this work we consider theoretical and numerical aspects of the non-singular theory of discrete dislocation loops in gradient elasticity of Helmholtz type, with interest in its applications to three dimensional dislocation dynamics (DD) simulations. The gradient solution is developed and compared to its singular and non-singular counterparts in classical elasticity using the unified framework of eigenstrain theory. The fundamental equations of curved dislocation theory are given as non-singular line integrals suitable for numerical implementation using fast one-dimensional quadrature. These include expressions for the interaction energy between two dislocation loops and the line integral form of the generalized solid angle associated with dislocations having a spread core. The single characteristic length scale of Helmholtz elasticity is determined from independent molecular statics (MS) calculations. The gradient solution is implemented numerically within our variational formulation of DD, with several examples illustrating the viability of the non-singular solution. The displacement field around a dislocation loop is shown to be smooth, and the loop self-energy non-divergent, as expected from atomic configurations of crystalline materials. The loop nucleation energy barrier and its dependence on the applied shear stress are computed and shown to be in good agreement with atomistic calculations. DD simulations of Lomer–Cottrell junctions in Al show that the strength of the junction and its configuration are easily obtained, without ad-hoc regularization of the singular fields. Numerical convergence studies related to the implementation of the non-singular theory in DD are presented.     

On canonical bending relationships for plates

In recent years, a series of papers have appeared on algebraic relationships between the solutions (e.g., deflections, buckling loads and frequencies) of a given higher-order plate theory and the classical plate theory. The bending relationships, for example, can be used to generate the transverse deflection of a plate according to the particular higher-order theory from the known deflection of the same problem according to the classical plate theory. In the present study relationships between the bending solutions of several higher-order plate theories and the classical plate theory are obtained in a canonical form (i.e., one set of relationships contain several theories and they can be specialized to a specific theory by assigning values to the constants appearing in the relationships). Numerical examples of bending solutions for rectangular plates with various boundary conditions are presented to show how the relations can be used to determine the deflections and bending moments for various theories. The relationships are validated by comparing the numerical results obtained using the relationships for the Mindlin plate theory against those computed using the ABAQUS finite element program.     

Finite element solutions for plane strain mode I crack with strain gradient effects

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom ω_i are introduced in addition to the conventional three translational degrees of freedom u_i. ω_i is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l₁. Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale.     

Forced oscillations and stability analysis of a nonlinear micro-rotating shaft incorporating a non-classical theory

In this paper, the stability and bifurcation analysis of symmetrical and asymmetrical micro-rotating shafts are investigated when the rotational speed is in the vicinity of the critical speed. With the help of Hamilton’s principle, nonlinear equations of motion are derived based on non-classical theories such as the strain gradient theory. In the dynamic modeling, the geometric nonlinearities due to strains, and strain gradients are considered. The bifurcations and steady state solution are compared between the classical theory and the non-classical theories. It is observed that using a non-classical theory has considerable effect in the steady-state response and bifurcations of the system. As a result, under the classical theory, the symmetrical shaft becomes completely stable in the least damping coefficient, while the asymmetrical shaft becomes completely stable in the highest damping coefficient. Under the modified strain gradient theory, the symmetrical shaft becomes completely stable in the least total eccentricity, and under the classical theory the asymmetrical shaft becomes completely stable in the highest total eccentricity. Also, it is shown that by increasing the ratio of the radius of gyration per length scale parameter, the results of the non-classical theory approach those of the classical theory.     

Analysis of micropolar elastic beams

In this paper, a linear theory for the analysis of beams based on the micropolar continuum mechanics is developed. Power series expansions for the axial displacement and micro-rotation fields are assumed. The governing equations are derived by integrating the momentum and moment of momentum equations in the micropolar continuum theory. Body couples and couple stresses can be supported in this theory. After some simplifications, this theory can be reduced to the well-known Timoshenko and Euler–Bernoulli beam theories. The nature of flexural and longitudinal waves in the infinite length micropolar beam has been investigated. This theory predicts the existence of micro-rotational waves which are not present in any of the known beam theories based on the classical continuum mechanics. Also, the deformation of a cantilever beam with transverse concentrated tip loading has been studied. The pattern of deflection of the beam is similar to the classical beam theories, but couple stress and micro-rotation show an oscillatory behavior along the beam for various loadings.     

An elementary molecular-statistical basis for the Mooney and Rivlin-Saunders theories of rubber elasticity

By relaxing the assumption that the end-to-end vectors of molecules transform as macroscopic material line elements, we arrive at a generalization of the molecular-statistical theory of rubber elasticity. This generalization includes as special cases continuum-mechanical theories proposed by Mooney and by Rivlin and Saunders as improvements upon the classical neo-Hookean theory.     

高雷诺数流动理论、算法和应用的若干研究进展

在黏性流体力学的历史发展中,Navier-Stokes (NS)方程组的 简化理论、相应算法和应用构成了不同历史时期的学科前沿、核心内 容的应用热点。以此为线索,简要评述经典边界层、多层(三层)边界 层、干扰边界层、扩散抛物化(DP) NS方程诸理论、相应算法和应用的 若干研究进展;诸理论之间以及他们与实验的关系;简化湍流计算的 一点注释以及物理分析和数值模拟相结合的一些问题。