n the fourth order tensor valued function of the stress in return map algorithm

针对各向同性材料,基于一组相互正交的基张量,建立了一套有 效的相关运算方法. 基张量中的两个分别是归一化的二阶单位张量和偏应力张量,另一个则 使用应力的各向同性二阶张量值函数经过归一化构造所得,三者共主轴. 根据张量函数表示 定理,本构方程和返回映射算法中所涉及到的应力的二阶、四阶张量值函数及其逆都由这组 基所表示. 推演结果表明：这些张量之间的运算,表现为对应系数矩阵之间的简单 关系. 其中,四阶张量求逆归结为对应的3\times3系数矩阵求逆,它对二阶张量的变换 则表现为该矩阵对3times 1列阵的变换. 最后,对这些变换关系应用于返回映 射算法的迭代格式进行了相关讨论.     

使用主值空间表示的各向同性塑性本构方程

针对各向同性材料,在内变量为标量的假定下,应用张量函数表示定理给出了其塑性应变增量的不变性表示.它的3个不可约基张量取决于应力张量、相互正交且共主轴.建立3个基张量构成的张量子空间与三维主值空间的对应关系,将共主轴的张量采用笛卡尔坐标系中的矢量描述,矢量在不同坐标系下的分量均为张量的一组不可约不变量.定义塑性应变增量对应的矢量为内变量增量,使用张量函数表示理论得到,内变量演化方程除取决于应力对应的矢量和内变量本身外,还取决于应力增量在张量子空间中的投影,该投影就是应力对应矢量的增量,因此,本构方程归结为确定主值空间中矢量之间的关系.最后表明,三维主值空间与张量子空间中的流动法则是等价的.     

非线性横观各向同性弹性材料的本构方程及其势函数

研究了非线性Green弹性材料弹性张量独立分量,归纳推导出横观各向同性Green弹性材料、各向同性非线性弹性材料独立的弹性常数个数.从张量函数出发,用含有高阶弹性张量的张量多项式,推导出四阶非线性横观各向同性,各向同性材料Green弹性材料本构方程及其势函数.并将本构方程及其势甬数用张量不变量,标量不变量表示.证明了方程是完备的,不可约的,满足张量函数表示定理.     

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

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

基于中间构形的大变形弹塑性模型

在大变形弹塑性本构理论中,一个基本的问题是弹性变形和塑性变形的分解.通常采用两种分解方式,一是将变形率(或应变率)加法分解为弹性和塑性两部分,其中,弹性变形率与Kirchhoff应力的客观率通过弹性张量联系起来构成所谓的次弹性模型,而塑性变形率与Kirchhoff应力使用流动法则建立联系;另一种是基于中间构形将变形梯度进行乘法分解,它假定通过虚拟的卸载过程得到一个无应力的中间构形,建立所谓超弹性-塑性模型.研究了基于变形梯度乘法分解并且基于中间构形的大变形弹塑性模型所具有的若干性质,包括:在不同的构形上,塑性旋率的存在性、背应力的对称性、塑性变形率与屈服面的正交性以及它们之间的关系.首先,使用张量函数表示理论,建立了各向同性函数的若干特殊性质,并导出了张量的张量值函数在中间构形到当前构形之间进行前推后拉的简单关系式.然后,基于这些特殊性质和关系式,从热力学定律出发,建立模型在不同构形上的数学表达,包括客观率表示的率形式和连续切向刚度等,从而获得模型所具有的若干性质.最后,将模型与4种其他模型进行了比较分析.      

基于中间构形的大变形弹塑性模型

在大变形弹塑性本构理论中,一个基本的问题是弹性变形和塑性变形的分解.通常采用两种分解方式,一是将变形率(或应变率)加法分解为弹性和塑性两部分,其中,弹性变形率与Kirchhoff应力的客观率通过弹性张量联系起来构成所谓的次弹性模型,而塑性变形率与Kirchhoff应力使用流动法则建立联系;另一种是基于中间构形将变形梯度进行乘法分解,它假定通过虚拟的卸载过程得到一个无应力的中间构形,建立所谓超弹性–塑性模型.研究了基于变形梯度乘法分解并且基于中间构形的大变形弹塑性模型所具有的若干性质,包括:在不同的构形上,塑性旋率的存在性、背应力的对称性、塑性变形率与屈服面的正交性以及它们之间的关系.首先,使用张量函数表示理论,建立了各向同性函数的若干特殊性质,并导出了张量的张量值函数在中间构形到当前构形之间进行前推后拉的简单关系式.然后,基于这些特殊性质和关系式,从热力学定律出发,建立模型在不同构形上的数学表达,包括客观率表示的率形式和连续切向刚度等,从而获得模型所具有的若干性质.最后,将模型与4种其他模型进行了比较分析.     

薄板分析的线性基梯度光滑伽辽金无网格法

薄板问题的控制方程为四阶微分方程,因而当采用伽辽金法进行分析时,形函数需要满足C$^{1}$连续性要求,且至少使用二次基函数才能保证方法的收敛性.无网格形函数虽然易于满足C$^{1}$连续性要求,但由于不是多项式,其二阶导数的计算较为复杂耗时,同时也对刚度矩阵的数值积分提出了更高的要求.本文提出了一种薄板分析的线性基梯度光滑伽辽金无网格法,该方法的基础是线性基无网格形函数的光滑梯度.在梯度光滑构造的理论框架内,无网格形函数的二阶光滑梯度可以表示为形函数一阶梯度的线性组合,因而可以提高形函数二阶梯度的计算效率.分析表明,线性基无网格形函数的光滑梯度不仅满足其固有的线性梯度一致性条件,还满足本属于二次基函数对应的额外高阶一致性条件,因此能够恰当地运用到薄板结构的伽辽金分析.此外,插值误差分析也很好地验证了线性基无网格光滑梯度的收敛特性.算例结果进一步表明,线性基梯度光滑伽辽金无网格法的收敛率与传统二次基伽辽金无网格法相当,但精度更高,同时刚度矩阵所需的高斯积分点数明显减少.      

关于张量函数表示理论的标量不变量的讨论

发现文献[1,2]提出的张量函数表示理论中的完备而不可约的不变量不是互相完全独立的．分别对一个任意二阶张量和两个对称二阶张量的标量不变量进行了计算,证明前者只有六个不变量是独立的,后者只有九个是独立的．     

薄板分析的线性基梯度光滑伽辽金无网格法

薄板问题的控制方程为四阶微分方程,因而当采用伽辽金法进行分析时,形函数需要满足C~1连续性要求,且至少使用二次基函数才能保证方法的收敛性.无网格形函数虽然易于满足C~1连续性要求,但由于不是多项式,其二阶导数的计算较为复杂耗时,同时也对刚度矩阵的数值积分提出了更高的要求.本文提出了一种薄板分析的线性基梯度光滑伽辽金无网格法,该方法的基础是线性基无网格形函数的光滑梯度.在梯度光滑构造的理论框架内,无网格形函数的二阶光滑梯度可以表示为形函数一阶梯度的线性组合,因而可以提高形函数二阶梯度的计算效率.分析表明,线性基无网格形函数的光滑梯度不仅满足其固有的线性梯度一致性条件,还满足本属于二次基函数对应的额外高阶一致性条件,因此能够恰当地运用到薄板结构的伽辽金分析.此外,插值误差分析也很好地验证了线性基无网格光滑梯度的收敛特性.算例结果进一步表明,线性基梯度光滑伽辽金无网格法的收敛率与传统二次基伽辽金无网格法相当,但精度更高,同时刚度矩阵所需的高斯积分点数明显减少.     

面向应力约束的独立连续映射方法

大多数已有的拓扑优化研究为系统刚度最大化设计,尤其以体积比约束下的静态柔顺度最小化问题为典型.从工程角度出发,结构强度设计至关重要.以往的应力研究表明,应力约束拓扑优化存在着奇异性、约束数目庞大、高度非线性特性等诸多数值困难.为了实现应力约束下的拓扑优化设计,采用归一化p范数应力指标以减少单元应力约束数目.遵循独立连续映射建模方式,引入密度变量的倒变量函数作为设计变量.推导了应力约束函数和体积目标函数对设计变量的敏度,并基于一阶和二阶泰勒近似得到各自的显式表达式.通过构造的系列二次规划子问题,原拓扑优化问题采用序列二次规划算法高效求解.二维数值算例考察了结构刚度和强度设计结果的异同,以及不同应力约束上限值对应力约束拓扑优化结果的影响.通过提出方法与传统变密度法结果的比较,说明提出的独立连续映射方法在应力约束下具有可行性和有效性.优化结果也表明了考虑应力约束的连续体拓扑优化具有必要性.      

General invariant representations of the constitutive equations for isotropic nonlinearly elastic materials

This paper develops general invariant representations of the constitutive equations for isotropic nonlinearly elastic materials. Different sets of mutually orthogonal unit tensor bases are constructed from the strain argument tensor by using the representation theorem and corresponding irreducible invariants are defined. Their relations and geometrical interpretations are established in three dimensional principal space. It is shown that the constitutive law linking the stress and strain tensors is revealed to be a simple relationship between two vectors in the principal space. Relative to two different sets of the basis tensors, the constitutive equations are transformed according to the transformation rule of vectors. When a potential function is assumed to exist, the vector associated with the stress tensor is expressed in terms of its gradient with respect to the vector associated with the strain tensor. The Hill’s stability condition is shown to be that the scalar product of the increment of those two vectors must be positive. When potential function exists, it becomes to be that the 3 × 3 constitutive matrix derived from its second order derivative with respect to the vector associated with the strain must be positive definite. By decomposing the second order symmetric tensor space into the direct sum of a coaxial tensor subspace and another one orthogonal to it, the closed form representations for the fourth order tangent operator and its inversion are derived in an extremely simple way.     

弹性大变形问题中的应力状态描述、奇异性问题和余能原理

对弹性大变形理论中的3方面问题进行了综述.首先,对各种应变度量的共轭应力进行综述.大变形问题引起的应力状态描述的复杂性引起了许多学者的兴趣,对这个问题的研究也促进了大变形弹性理论的发展.在各种特定问题中,人们提出了不同的应力张量来描述应力状态,如Caucby应力张量、第一类和第一二类Piola-Kirchhoff应力张...     

Objective Tensor Rates and Applications in Formulation of Hyperelastic Relations

Following Ogden, a class of objective (Lagrangian and Eulerian) tensors is identified among the second-rank tensors characterizing continuum deformation, but a more general definition of objectivity than that used by Ogden is introduced. Time rates of tensors are determined using convective rates. Sufficient conditions of objectivity are obtained for convective rates of objective tensors. Objective convective rates of strain tensors are used to introduce pairs of symmetric stress and strain tensors conjugate in a generalized sense. The classical definitions of conjugate Lagrangian (after Hill) and Eulerian (after Xiao et al.) stress and strain tensors are particular cases of the definition of conjugacy of stress and strain tensors in the generalized sense used in the present paper. Pairs of objective stress and strain tensors conjugate in the generalized sense are used to formulate constitutive relations for a hyperelastic medium. A family of objective generalized strain tensors is introduced, which is broader than Hill’s family of strain tensors. The basic forms of the hyperelastic constitutive relations are obtained with the aid of pairs of Lagrangian stress and strain tensors conjugate after Hill (the strain tensors in these pairs belong to the family of generalized strain tensors). A method is presented for generating reduced forms of the constitutive relations with the aid of pairs of Lagrangian and Eulerian stress and strain tensors conjugate in the generalized sense which are obtained from pairs of Lagrangian tensors conjugate after Hill by mapping tensor fields on one configuration of a deformable body to tensor fields on another configuration.      

Application of corotational rates of the logarithmic strain in constitutive modeling of hardening materials at finite deformations

In this paper a finite deformation constitutive model for rigid plastic hardening materials based on the logarithmic strain tensor is introduced. The flow rule of this constitutive model relates the corotational rate of the logarithmic strain tensor to the difference of the deviatoric Cauchy stress and the back stress tensors. The evolution equation for the kinematic hardening of this model relates the corotational rate of the back stress tensor to the corotational rate of the logarithmic strain tensor. Using Jaumann, Green–Naghdi, Eulerian and logarithmic corotational rates in the proposed constitutive model, stress–strain responses and subsequent yield surfaces are determined for rigid plastic kinematic and isotropic hardening materials in the simple shear problem at finite deformations.     

Macro–micro relations in granular mechanics

In granular mechanics, macroscopic approaches treat a granular material as an equivalent continuum at macro-scale, and study its constitutive relationship between macro-quantities, such as stresses and strains. On the other hand, microscopic approaches consider a granular material as an assembly of individual particles interacting with each other at micro-scale (i.e., particle-scale), and the physical quantities under study are forces and displacements. This paper aims at linking up the findings from these two scales and to establish the macro–micro relations in granular mechanics.Three aspects of the macro–micro relations are investigated. They are about the internal structure, the stress tensor and the strain tensor. The internal structure is described with geometrical systems at the particle scale. Micro-structural definitions of the stress and strain tensors are derived, which link the macro-stress tensor with the contact forces, and the macro-strain tensor with the relative displacements at contact. In addition to a brief review of the past research work on these topics, further generalizations are made in this paper. In particular, the two cell systems proposed by Li and Li (2009), namely the solid cell system and the void cell system, are introduced and used for the derivation of the macro-structural expressions. The stress expression is derived based on Newton’s second law of motion. The result is valid for both static and dynamic cases. The strain expression is derived based on the compatibility requirement. And the expression is valid for any tessellation subdividing the granular assembly into polyhedral elements.The homogenization for deriving a macroscopic constitutive relationship from microscopic behaviour is discussed. Attention is placed on the macroscopic quantification of the internal structure in terms of a second rank tensor, known as the fabric tensor. Existing definitions of the fabric tensors have been reviewed. The correlations among different fabric tensors and their relations with the stress–strain behaviours have been investigated.     

Flows of constant stretch history for polymeric materials with power-law distributions of relaxation times

It is herein shown that for separable integral constitutive equations with power-law distributions of relaxation times, the streamlines in creeping flow are independent of flow rate.For planar flows of constant stretch history, the stress tensor is the sum of three terms, one proportional to the rate-of-deformation tensor, one to the square of this tensor, and the other to the Jaumann derivative of the rate-of-deformation tensor. The three tensors are the same as occur in the Criminale-Ericksen-Filbey Equation, but the coefficients of these tensors depend not only on the second invariant of the strain rate, but also on another invariant which is a measure of flow strength. With the power-law distribution of relaxation times, each coefficient is equal to the second invariant of the strain rate tensor raised to a power, times a function that depends only on strength of the flow. Axisymmetric flows of constant stretch history are more complicated than the planar flows, because three instead of two nonzero normal components appear in the velocity gradient tensor. For homogeneous axisymmetric flows of constant stretch history, the stress tensor is given by the sum of the same three terms. The coefficients of these terms again depend on the flow strength parameter, but in general the dependences are not the same as in planar flow.