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

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

以Atluri应力和左伸长张量为共轭变量的有限弹性广义变分原理

本文提出了以Atluri应力张量和左伸长张量为共轭应力应变量的有限弹性变形广义变分原理。     

拉格朗日几何非线性混合应变元

本文应用由Ｓｉｍｏ和Ｒｉｆａｉ建议的混合假定附加应变途径，采用第二Ｐｉｏｌａ－Ｋｉｒｃｈｈｏｆｆ应力张量和Ｇｒｅｅｎ－Ｌａｇｒａｎｇｅ应变张量作能量共轭的应力应变度量，导出了Ｌａｇｒａｎｇｅ几何非线性下的胡海昌－Ｗａｓｈｉｚｕ三变量变分原理的Ｇａｌｅｒｋｉｎ形式以及相应的混合假定应变元公式。     

旋转湍流中雷诺应力模型压力应变关联项的修正研究

利用张量的不变量理论,推导得出传统雷诺应力模型中压力应变关联项模型应用于旋转湍流模拟中的一些基本问题,即在纯旋转条件下,传统模型所描述的初始各向异性的湍流中雷诺应力张量演化规律是一个无衰减振荡过程,而快速畸变理论推导结果显示,其演化应是一个阻尼振荡衰减的过程。以衰减雷诺应力为目的,构造出包含旋转率张量高阶量的关联项。然后,结合变形率张量的高阶项,将修正模型扩展至椭圆形流线类型流动。最后,将修正模型应用于轴向旋转圆管内湍流流场的模拟,并将结果与实测结果进行了对比。     

含微裂纹材料的损伤理论

本文从含微裂纹材料的变形能出发引出了裂纹的方位张量。在考虑裂纹受压闭合与滑动摩擦的基础上,给出了损伤张量、损伤应变及有效弹性常数。文中给出了损伤机构离散化的方法,并对方位密度给出了演化方程。最后给出一个单向拉压的应力应变关系例子,并揭示了裂纹扩展时的应力突跌现象。     

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

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

Lode角对应力张量的导数的基本性质

为计及岩土类材料塑性力学行为的中主应力影响或应力路径相关性,通常将应力张量Lode角/Lode数引入屈服函数与塑性势函数。由此在计算塑性应变增量时必然涉及Lode角/Lode数对应力的导数张量（记为 ）。然而,应力张量主值有重根时 的计算存在困难。本文给出了 的主值计算方法及谱分解表达式并详细讨论了张量 的基本性质。     

以Jaumann应力和右伸长张量为共轭变量的一对广义变分原理

本文提出了以Jaumann应力和右伸长张量作为共轭应力应变变量的一对非线性广义变分原理,其中一个是势能原理,另一个是余能原理。     

连续介质力学中变形梯度张量客观性表述的分歧

张量的客观性是连续介质力学中一个重要的概念,但现有文献对张量客观性的定义不一致导致有关变形梯度张量客观性的表述存在分歧.本文主要基于张量的逆及功共轭角度分析了不同客观性定义的差别,旨在加深对张量客观性,特别是对变形梯度等两点张量客观性的认识.     

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

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

The invariant representation of spins with applications in the theory of finite deformations

Based on the general solution given to a kind of linear tensor equations, the spin of a symmetric tensor is derived in an invariant form. The result is applied to find the spins of the left and the right stretch tensors and the relation among different rotation rate tensors has been discussed. According to work conjugacy, the relations between Cauchy stress and the stresses conjugate to Hill's generalized strains are obtained. Particularly, the logarithmic strain, its time rate and the conjugate stress have been discussed in detail. These results are important in modeling the constitutive relations for finite deformations in continuum mechanics. The project is supported by the National Natural Science Foundation of China and the Chinese Academy of Sciences (No. 87-52).     

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.      

A formulation of anisotropic continuum elastoplasticity at finite strains. Part I: Modelling

A constitutive model for anisotropic elastoplasticity at finite strains is developed together with its numerical implementation. An anisotropic elastic constitutive law is described in an invariant setting by use of structural tensors and the elastic strain measure C_e. The elastic strain tensor as well as the structural tensors are assumed to be invariant in relation to superimposed rigid body rotations. An anisotropic Hill-type yield criterion, described by a non-symmetric Eshelby-like stress tensor and further structural tensors, is developed, where use is made of representation theorems for functions with non-symmetric arguments. The model also considers non-linear isotropic hardening. Explicit results for the specific case of orthotropic anisotropy are given. The associative flow rule is employed and the features of the inelastic flow rule are discussed in full. It is shown that the classical definition of the plastic material spin is meaningless in conjunction with the present formulation. Instead, the study motivates an alternative definition, which is based on the demand that such a quantity must be dissipation-free, as the plastic material spin is in the case of isotropy. Equivalent spatial formulations are presented too. The full numerical treatment is considered in Part II.     

An elasto-plastic theory of dislocation and disclination fields

A linear theory of the elasto-plasticity of crystalline solids based on a continuous representation of crystal defects – dislocations and disclinations – is presented. The model accounts for the translational and rotational aspects of lattice incompatibility, respectively associated with the presence of dislocations and disclinations. The defects content relates to the incompatible plastic strain and curvature tensors. The stress state is described by using the conjugate variables to strain and curvature, i.e., the stress and couple-stress tensors. Defect motion is described by two transport equations. A dynamic interplay between dislocations and disclinations results from a disclination-induced source term in the transport of dislocations. Thermodynamic guidance provides the driving forces conjugate to dislocation and disclination velocity in a continuous context, as well as admissible constitutive relations for the latter. When dislocation and disclination velocity vanish, the model reduces to deWit’s elasto-static theory of crystal defects. It also reduces to Acharya’s linear elasto-plastic theory for dislocation fields when the disclination density is ignored. The theory is intended for use in instances where rotational defects matter, such as grain boundaries. To illustrate its applicability, a finite high-angle tilt boundary is modeled using a disclination dipole and its behavior under tensile loading normal to the boundary is shown.     

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

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

A non-ordinary state-based peridynamic method to model solid material deformation and fracture

In this paper, we develop a new non-ordinary state-based peridynamic method to solve transient dynamic solid mechanics problems. This new peridynamic method has advantages over the previously developed bond-based and ordinary state-based peridynamic methods in that its bonds are not restricted to central forces, nor is it restricted to a Poisson’s ratio of 1/4 as with the bond-based method. First, we obtain non-local nodal deformation gradients that are used to define nodal strain tensors. The deformation gradient tensors are used with the nodal strain tensors to obtain rate of deformation tensors in the deformed configuration. The polar decomposition of the deformation gradient tensors are then used to obtain the nodal rotation tensors which are used to rotate the rate of deformation tensors and previous Cauchy stress tensors into an unrotated configuration. These are then used with conventional Cauchy stress constitutive models in the unrotated state where the unrotated Cauchy stress rate is objective. We then obtain the unrotated Cauchy nodal stress tensors and rotate them back into the deformed configuration where they are used to define the forces in the nodal connecting bonds. As a first example we quasi-statically stretch a bar, hold it, and then rotate it ninety degrees to illustrate the methods finite rotation capabilities. Next, we verify our new method by comparing small strain results from a bar fixed at one end and subjected to an initial velocity gradient with results obtained from the corresponding one-dimensional small strain analytical solution. As a last example, we show the fracture capabilities of the method using both a notched and un-notched bar.     

Large strain plastic deformation by crystallographic slio in metals

Large strain plastic deformation of metallic materials is investigated. The rate of change of the resolved shear stress applied to a slip system in a grain is expressed as a linear function of the slip rates. For this purpose, using an elastic-plastic analysis of the inclusion problem, the Jaumann rate of change of the stress state applied to the grain, together with its total spin, are obtained as a function of its deformation rate in the rigid-plastic case. It is shown that the existing stress has a significant influence on the accommodation tensors involved in the solution.     

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.