A new modified couple stress theory for anisotropic elasticity and microscale laminated Kirchhoff plate model

A new modified couple stress theory for anisotropic elasticity is proposed. This theory contains three material length scale parameters. Differing from the modified couple stress theory, the couple stress constitutive relationships are introduced for anisotropic elasticity, in which the curvature (rotation gradient) tensor is asymmetric and the couple stress moment tensor is symmetric. However, under isotropic case, this theory can be identical to modified couple stress theory proposed by Yang et al. (Int J Solids Struct 39:2731–2743, 2002). The differences and relations of standard, modified and new modified couple stress theories are given herein. A detailed variational formulation is provided for this theory by using the principle of minimum total potential energy. Based on the new modified couple stress theory, composite laminated Kirchhoff plate models are developed in which new anisotropic constitutive relationships are defined. The First model contains two material length scale parameters, one related to fiber and the other related to matrix. The curvature tensor in this model is asymmetric; however, the couple stress moment tensor is symmetric. Under isotropic case, this theory can be identical to the modified couple stress theory proposed by Yang et al. (Int J Solids Struct 39:2731–2743, 2002). The present model can be viewed as a simplified couple stress theory in engineering mechanics. Moreover, a more simplified model of couple stress theory including only one material length scale parameter for modeling the cross-ply laminated Kirchhoff plate is suggested. Numerical results show that the proposed laminated Kirchhoff plate model can capture the scale effects of microstructures.     

Statics and kinematics of discrete Cosserat-type granular materials

A theoretical framework is presented for the statics and kinematics of discrete Cosserat-type granular materials. In analogy to the force and moment equilibrium equations for particles, compatibility equations for closed loops are formulated in the two-dimensional case for relative displacements and relative rotations at contacts. By taking moments of the equilibrium equations, micromechanical expressions are obtained for the static quantities average Cauchy stress tensor and average couple stress tensor. In analogy, by taking moments of the compatibility equations, micromechanical expressions are obtained for the (infinitesimal) kinematic quantities average rotation gradient tensor and average Cosserat strain tensor in the two-dimensional case. Alternatively, these expressions for the average Cauchy stress tensor and the average couple stress tensor are obtained from considerations of the equivalence of the continuum force and couple traction vectors acting on a plane and the resultant of the discrete forces and couples acting on this plane. In analogy, the expressions for the average rotation gradient tensor and the average Cosserat strain tensor are obtained from considerations of the change of length and change of rotation of a line element in the two-dimensional case. It is shown that the average particle stress tensor is always symmetrical, contrary to the average stress tensor of an equivalent homogenized continuum. Finally, discrete analogues of the virtual work and complementary virtual work principles from continuum mechanics are derived.     

修正的偶应力线弹性理论及广义线弹性体的有限元方法

以含偶应力的弹性理论为基础,考虑小变形情况下变形体的平动变形和旋转变形,提出关于偶应力与曲率张量的线性本构关系,建立一般弹性体的线性模型。为满足有限单元C1连续性要求,考虑转角为独立变量,利用罚方法引入约束条件,构造一般弹性体的约束变分形式。应用8节点48个自由度的实体等参元,建立一般弹性体力学响应分析的有限元方程。对悬臂梁的静力和动力分析表明,一般弹性体模型较之经典弹性力学更适合结构分析;较之Timoshenko梁模型,一般弹性体模型能够计及结构尺度对结构动力特性和动力响应造成的显著影响。     

基于Hellinger-Reissner变分原理的应变梯度杂交元设计

从一般的偶应力理论出发，基于Hellinger-Reissner变分原理，通过对有限元 离散体系的位移试解引入非协调位移函数，得到了偶应力理论下有限元离散系统的能量相容 条件，并由此建立了应变梯度杂交元的应力函数优化条件. 根据该优化条件，构造了一 个C0类的平面4节点梯度杂交元，数值结果表明，该单元对可压缩和不可压缩状态的 梯度材料均可给出合理的数值结果，再现材料的尺度效应.     

A physico-mechanical approach to modeling of metal forming processes—part I: theoretical framework

A combined physico-mechanical approach to research and modeling of forming processes for metals with predictable properties is developed. The constitutive equations describing large plastic deformations under complex loading are based on both plastic flow theory and continuum damage mechanics. The model which is developed in order to study strongly plastically deformed materials represents their mechanical behavior by taking micro-structural damage induced by strain micro-defects into account. The symmetric second-rank order tensor of damage is applied for the estimation of the material damage connected with volume, shape, and orientation of micro-defects. The definition offered for this tensor is physically motivated since its hydrostatic and deviatoric parts describe the evolution of damage connected with a change in volume and shape of micro-defects, respectively. Such a representation of damage kinetics allows us to use two integral measures for the calculation of damage in deformed materials. The first measure determines plastic dilatation related to an increase in void volume. A critical amount of plastic dilatation enables a quantitative assessment of the risk of fracture of the deformed metal. By means of an experimental analysis we can determine the function of plastic dilatation which depends on the strain accumulated by material particles under various stress and temperature-rate conditions of forming. The second measure accounts for the deviatoric strain of voids which is connected with a change in their shape. The critical deformation of ellipsoidal voids corresponds to their intense coalescence and to formation of large cavernous defects. These two damage measures are important for the prediction of the meso-structure quality of metalware produced by metal forming techniques. Experimental results of various previous investigations are used during modeling of the damage process.      

On non-linear Beltrami-Mitchell equations

Beltrami-Mitchell equations for non-linear elasticity theory are derived using the first Piola-Kirchhoff stress and the deformation gradient tensors as field variables so as to yield linear equilibrium and compatibility equations, respectively. In the derivation it is assumed that a strain energy density and, correspondingly, a complementary strain energy density exist, and satisfy the axiom of objectivity. Substitution for the deformation gradient in the compatibility equations yields non-linear differential equations in terms of the first Piola-Kirchhoff stress tensor which may be regarded as the Beltrami-Mitchell equations of non-linear elasticity. The equations are also derived for “semi-linear” isotropic elastic materials and the theory is illustrated by three simple examples.     

Nonlinear Eulerian thermoelasticity for anisotropic crystals

A complete continuum thermoelastic theory for large deformation of crystals of arbitrary symmetry is developed. The theory incorporates as a fundamental state variable in the thermodynamic potentials what is termed an Eulerian strain tensor (in material coordinates) constructed from the inverse of the deformation gradient. Thermodynamic identities and relationships among Eulerian and the usual Lagrangian material coefficients are derived, significantly extending previous literature that focused on materials with cubic or hexagonal symmetry and hydrostatic loading conditions. Analytical solutions for homogeneous deformations of ideal cubic crystals are studied over a prescribed range of elastic coefficients; stress states and intrinsic stability measures are compared. For realistic coefficients, Eulerian theory is shown to predict more physically realistic behavior than Lagrangian theory under large compression and shear. Analytical solutions for shock compression of anisotropic single crystals are derived for internal energy functions quartic in Lagrangian or Eulerian strain and linear in entropy; results are analyzed for quartz, sapphire, and diamond. When elastic constants of up to order four are included, both Lagrangian and Eulerian theories are capable of matching Hugoniot data. When only the second-order elastic constant is known, an alternative theory incorporating a mixed Eulerian–Lagrangian strain tensor provides a reasonable approximation of experimental data.     

Basic Issues Concerning Finite Strain Measures and Isotropic Stress-Deformation Relations

It is indicated that the commonly-used Rivlin–Ericksen representation formula for isotropic tensor functions exhibits some properties that might be undesirable for its reasonable and effective applications. Towards clarification and improvement, a set of three mutually orthogonal tensor generators is introduced to achieve an alternative representation formula for isotropic symmetric tensor-valued functions of a symmetric tensor. This representation formula enables us to express the unknown representative coefficients in terms of simple, explicit tensorial inner products of the argument tensor and the value tensor without involving their eigenvalues. In particular, the tensorial interpolation expressions thus obtained assume a unified form for the three different cases of coalescence of the eigenvalues of the argument tensor. Moreover, each summand in the alternative representation formula is shown to inherit the continuity and differentiability properties of the represented isotropic tensor function. These results are used to study some basic issues concerning finite strain measures and stress-deformation relations of isotropic materials, such as continuity and differentiability properties of the representation, determination of the representative coefficients in terms of experimental data for stress and deformation tensors, and computations of finite strain measures. This revised version was published online in June 2006 with corrections to the Cover Date.     

An asymmetric theory of nonlocal elasticity—Part 1. Quasicontinuum theory

In this article, an asymmetric theory of nonlocal elasticity is developed on the basis of three dimensional atomic lattice model, the Galileo invariance for constitutive equations and by use of Fourier transformation of generalized function and energy method. It is shown that nonlocal characteristic functions (or constitutive parameters of internal elastic energy) can be explicitly expressed in terms of interacting forces connecting atoms, and the general model of nonlocal theory with rotation effects is asymmetric. Both symmetric stress and anti-symmetric stress is a nonlocal function of strain and local rotation for anisotropic materials. For isotropic materials, symmetric stress is only a nonlocal function of strain, while antisymmetric stress is only a nonlocal function of local rotation.     

Curvature and torsion of material lines in chaotic flows

The deformation of a material line as it evolves in a chaotic flow is considered. Of particular interest are the curvature and torsion as a function of arclength or material coordinate along the line and their near singular behavior that develops as the deformation proceeds. Recent work [Thiffeault 2004. Physica D 198, 169–181] has established an interesting connection between local line stretch and curvature for regions of high curvature. The present effort confirms this connection and develops a similar relation between torsion and curvature for regions of high absolute torsion. In addition it is shown that certain properties of the strain tensor play a key role in the location and structure of the near singular behavior in curvature.     

Couple stresses in crystalline solids: origins from plastic slip gradients, dislocation core distortions, and three-body interatomic potentials

In the presence of plastic slip gradients, compatibility requires gradients in elastic rotation and stretch tensors. In a crystal lattice the gradient in elastic rotation can be related to bond angle changes at cores of so-called geometrically necessary dislocations. The corresponding continuum strain energy density can be obtained from an interatomic potential that includes two- and three-body terms. The three-body terms induce restoring moments that lead to a couple stress tensor in the continuum limit. The resulting stress and couple stress jointly satisfy a balance law. Boundary conditions are obtained upon stress, couple stress, strain and strain gradient tensors. This higher-order continuum theory was formulated by Toupin (Arch. Ration. Mech. Anal. 11 (1962) 385). Toupin's theory has been extended in this work to incorporate constitutive relations for the stress and couple stress under multiplicative elastoplasticity. The higher-order continuum theory is exploited to solve a boundary value problem of relevance to single crystal and polycrystalline nano-devices. It is demonstrated that certain slip-dominated deformation mechanisms increase the compliance of nanostructures in bending-dominated situations. The significance of these ideas in the context of continuum plasticity models is also dwelt upon.     

The influence of particle fluctuations on the average rotation in an idealized granular material

Granular materials are a simple example of a Cosserat continuum in that the average particle rotations may differ from the rotation of the average deformation. In the absence of couple stress, this difference insures that the stress is symmetric. This has been shown in theories that assume that the displacement at particle contacts is given by the average deformation and spin. Here, we indicate how the difference between the average rotation of the particles and the average rotation of the deformation can be determined when fluctuations in particle displacements and rotations satisfy local force and moment equilibria in a random aggregate of identical spheres. The predictions based on this model are in better agreement with numerical simulation than that given by the simple average strain theory.     

偶应力理论下拉力型锚固体的受力特性

锚固体的受力特征及其影响因素是锚固体设计的重要依据,直接影响锚固效果。传统的经典弹性理论没有考虑应变梯度的影响。偶应力理论引进弯曲曲率,考虑了弯曲效应对介质变形特性的影响。基于偶应力理论,建立了平面应变问题的有限元计算模型,研究锚固体锚固段界面上的剪应力分布、锚固体轴力分布、偶应力的尺度效应以及弹性模量和围压对锚固力的影响,并将偶应力理论的计算结果和经典弹性理论的计算结果进行了比较。结果表明,在偶应力理论下,锚固体锚固段界面的剪应力有所减小,特别是峰值处的剪应力减小明显;岩土的弹性模量越大,锚固界面局部剪应力越大;锚固力随着围压的增大而增大,偶应力尺度效应明显。     

近场动力学框架下的键基对应模型

传统键基近场动力学模型存在泊松比限制的问题,为了解决这一问题发展了态基近场动力学模型。其中非常规态的近场动力学模型通过定义非局部的变形梯度将近场力和传统应力关联起来,方便使用传统本构,但是态基近场动力学计算效率低于键基近场动力学。结合态基模型和键基模型的优势,提出键基对应模型,定义了基于键的变形梯度,参考连续介质力学中变形梯度的极分解过程,将键的变形分为转动部分和伸长部分。从而进一步定义了应变,通过物理方程求应力,进而计算键传递的近场力。键基对应模型解决了键基近场动力学的泊松比限制问题,也不需要进行近场动力学微观材料常数的计算。数值算例和理论推导证明了键变形梯度定义以及近场力计算方式的正确性。     

A finite strain kinematic hardening constitutive model based on Hencky strain: General framework, solution algorithm and application to shape memory alloys

The logarithmic or Hencky strain measure is a favored measure of strain due to its remarkable properties in large deformation problems. Compared with other strain measures, e.g., the commonly used Green-Lagrange measure, logarithmic strain is a more physical measure of strain. In this paper, we present a Hencky-based phenomenological finite strain kinematic hardening, non-associated constitutive model, developed within the framework of irreversible thermodynamics with internal variables. The derivation is based on the multiplicative decomposition of the deformation gradient into elastic and inelastic parts, and on the use of the isotropic property of the Helmholtz strain energy function. We also use the fact that the corotational rate of the Eulerian Hencky strain associated with the so-called logarithmic spin is equal to the strain rate tensor (symmetric part of the velocity gradient tensor). Satisfying the second law of thermodynamics in the Clausius-Duhem inequality form, we derive a thermodynamically-consistent constitutive model in a Lagrangian form. In comparison with the available finite strain models in which the unsymmetric Mandel stress appears in the equations, the proposed constitutive model includes only symmetric variables. Introducing a logarithmic mapping, we also present an appropriate form of the proposed constitutive equations in the time-discrete frame. We then apply the developed constitutive model to shape memory alloys and propose a well-defined, non-singular definition for model variables. In addition, we present a nucleation-completion condition in constructing the solution algorithm. We finally solve several boundary value problems to demonstrate the proposed model features as well as the numerical counterpart capabilities.     

A robust hyper-elastic beam model under bi-axial normal-shear loadings

In this paper, a simple and robust constitutive model is proposed to simulate mechanical behaviors of hyper-elastic materials under bi-axial normal-shear loadings in the finite strain regime. The Mooney–Rivlin strain energy function is adopted to develop a two-dimensional (2D) normal-shear constitutive model within the framework of continuum mechanics. A motion field is first proposed for combined normal and shear deformations. The deformation gradient of the proposed field is calculated and then substituted into right Cauchy–Green deformation tensor. Constitutive equations are then derived for normal and shear deformations. They are two explicit coupled equations with high-level polynomial non-linearity. In order to examine capabilities of the developed hyper-elastic model, uniaxial tensile responses and non-linear stability behaviors of moderately thick straight and curved beams undergoing normal axial and transverse shear deformations are simulated and compared with experiments. Fused deposition modeling technique as a 3D printing technology is implemented to fabricate hyper-elastic beam structures from soft poly-lactic acid filaments. The printed specimens are tested under tensile/compressive in-plane and compressive out-of-plane forces. A finite element formulation along with the Newton–Raphson and Riks techniques is also developed to trace non-linear equilibrium path of beam structures in large defamation regimes. It is shown that the model is capable of predicting non-linear equilibrium characteristics of hyper-elastic straight and curved beams. It is found that the modeling of shear deformation and finite strain is essential toward an accurate prediction of the non-linear equilibrium responses of moderately thick hyper-elastic beams. Due to simplicity and accuracy, the model can serve in the future studies dealing with the analysis of hyper-elastic structures in which two normal and shear stress components are dominant.     

Modeling deformation behavior of polymers with viscoplasticity theory based on overstress

The nonlinear strain rate sensitivity, multiple creep and recovery behavior of polyphenylene oxide (PPO), which were explored through strain rate-controlled experiments at ambient temperature by Khan [The deformation behavior of solid polymers and modeling with the viscoplasticity theory based overstress, Ph.D. Thesis, Rensselaer Polytechnic Institute, New York], are modeled using the modified viscoplasticity theory based on overstress (VBO). In addition, VBO used by Krempl and Ho [An overstress model for solid polymer deformation behavior applied to Nylon 66, ASTM STP 1357, 2000, p. 118] and the classical VBO are used to demonstrate the improved modeling capabilities of VBO for solid polymer deformation. The unified model (VBO) has two tensor valued state variables, the equilibrium and kinematic stresses and two scalar valued states variables, drag and isotropic stresses. The simulations include monotonic loading and unloading at various strain rates, multiple creep and recovery at zero stress. Since creep behavior has been found to be profoundly influenced by the level of the stress, the tests are performed at different stresses above and below the yield point. Numerical results are compared to experimental data. It is shown that nonlinear rate sensitivity, nonlinear unloading, creep and recovery at zero stress can be reproduced using the modified viscoplasticity theory based on overstress.     

各向同性率无关材料本构关系的不变性表示

在内变量理论的框架下,针对各向同性率无关材料,使用张量函数表示理论建立了塑性应变全量及增量本构关系的最一般的张量不变性表示. 它们均由3个完备不可约的基张量组合构成,这3个基张量分别是应力的零次幂、一次幂和二次幂. 因此得出,塑性应变、塑性应变增量与应力三者共主轴. 通过对基张量的正交化,给出了本构关系式在主应力空间中的几何解释. 进一步,全量(或增量)本构关系中3个组合因子被表达为应力、塑性应变(或塑性应变增量)的不变量的函数. 当塑性应变(或塑性应变增量)的3个不变量之间满足一定关系时,所给出的本构关系将退化为经典的形变理论(或塑性势理论).最后,还讨论它与奇异屈服面理论的关系,当满足一定条件时,两者是一致的.