一种描述形状记忆合金拟弹性变形行为的本构关系

本文给出一种描述形状记忆合金拟弹性变形现象的本构关系,可用于多晶材料在一般应力状态下及单晶材料在单轴应力下的变形情况。该本构关系采用弹性应变与相变应变迭加形式,物理意义明显,形式简洁。对 Cu-Zn-Sn 合金及 Ti-Ni 合金材料的变形行为进行了模拟计算,结果与实验值有较好的吻合。     

一种描述形状记忆合金拟弹性变形行为的本构关系

本文给出一种描述形状记忆合金拟弹性变形现象的本构关系,可用于多晶材料在一般应力状态下及单晶材料在单轴应力下的变形情况。该本构关系采用弹性应变与相变应变迭加形式,物理意义明显,形式简洁。对 Cu-Zn-Sn 合金及 Ti-Ni 合金材料的变形行为进行了模拟计算,结果与实验值有较好的吻合。     

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

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

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

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

应变相关三维蠕变损伤和破坏的热力学基础

应用热力学定律阐明了非等温多维蠕变损伤和破坏的概念。从以应力、温度和损伤诸函数的时间积分表示的Gibbs自由能泛函出发,并借助热力学第二定律,导出了关于应变、熵和损伤的本构律,它们一般足以包括现有的各种理论。然后,选取更具体的Gibbs自由能形式,得到了应变相关损伤理论的控制本构律。此外还对耦合和准静态非耦合两种情况,讨论了应力、位移、温度和损伤的边值问题的表述。     

一种新的应力率定义

本文研究客观应力率的定义及表示,在放弃了客观应力率形式上的对称性以后定义了一类广义的非对称的客观应力率。从Cauchy应力原理出发,通过构造应力向量形式上和变形及变形率无关的客观导数,得到两种形式的应力率;它们表示作用在参考构形中单位物质面元应力向量的真实变化,这样就消除了应力率定义中的不确定性。用新应力率对次弹性体简单剪切进行了计算,得到的应力不振荡,且满足超弹性正应力-剪应力普适关系。     

生长对超弹性球壳变形和稳定性的影响

应用连续介质力学有限变形理论建立受内压作用不可压超弹性球壳大变形问题的力学模型, 且运用基于变形梯度张量极分解的弹性体积生长理论分析生长对不可压超弹性球壳变形和稳定性的影响. 通过对球壳变形与内压关系式的数值计算得到不同生长条件下球壳的变形曲线和应力分布曲线及由生长引起的残余变形和残余应力分布. 计算结果表明生长对球壳变形特性有明显的影响, 生长影响球壳可产生不稳定变形的临界壁厚和临界内压, 且在某些情况下可改变球壳的稳定性. 生长在球壳中可产生一定的残余变形和残余应力, 对球壳中的应力分布有一定的影响. 另外当生长的程度足够大时, 即便没有外力作用,球壳仅在生长引起的残余应力作用下也可产生不稳定变形.      

有机玻璃在高应变率下的损伤型非线性粘弹性本构关系及破坏准则

对两种有机玻璃高应变率下的大变形和破坏行为进行了实验研究,通过改进文献[2]本构关系的非线性弹性项并引入损伤参量,建立了一个适用于更大变形范围、能描述应力平台及本构失稳的损伤型非线性粘弹性本构方程。相应地,从临界损伤量概念出发,提出以应变和应变率为控制变量的破坏准则。不论是本构关系还是破坏准则,理论计算均与试验结果吻合良好。     

基于对数应力率大应变本构关系的参数标定研究

在所有率型弹塑性本构模型中,只有对数应力率对应的本构模型能够满足自适应准则.基于对数应力率,采用实心圆轴扭转实验,对大应变弹塑性本构模型中的参数标定问题进行了讨论.推导出了考虑Swift效应时端部自由实心圆轴扭转变形的变形率、对数旋率、Kirchhoff应力及Kirchhoff应力的对数应力率.对于等向强化大应变弹塑性本构关系,给出了由实心圆轴扭转实验标定的、基于Kirchhhoff应力对数应力率的本构关系中塑性刚度函数的表达式.分析了扭转圆轴的Swift效应对塑性刚度函数的影响.结果表明,实心圆轴扭转的轴向伸长变形和径向变形对基于对数应力率大应变本构关系中的塑性刚度函数都有影响.当不考虑Swift效应时,所得塑性刚度函数表达式与不考虑Swift效应时基于Jaumann应力率的塑性刚度函数表达式相同.     

硬质聚氨酯泡沫塑料本构关系的研究

介绍用大尺寸分离式Ｈｏｐｋｉｎｓｏｎ压杆对四种密度的硬质聚氨酯泡沫塑料进行高应变率实验，完整地给出了这种材料在１０３／ｓ高应变率下的包括弹性区、屈服区和致密区变形全过程的动态应力应变曲线，并提出了包括应力、应变、应变率和密度等参量的本构关系．     

Corotational rates in constitutive modeling of elastic-plastic deformation

The principal axes technique is used to develop a new hypoelastic constitutive model for an isotropic elastic solid in finite deformation. The new model is shown to produce solutions that are independent of the choice of objective stress rate. In addition, the new model is found to be equivalent to the isotropic finite elastic model; this is essential if both models describe the same material.

The new hypoelastic model is combined with an isotropic flow rule to form an elastic-plastic rate constitutive equation. Use of the principal axes technique ensures that the stress tensor is coaxial with the elastic stretch tensor and that solutions do not depend on the choice of objective stress rate. The flow rule of von Mises and a parabolic hardening law are used to provide an example of application of the new theory. A solution is obtained for the prescribed deformation of simple rectilinear shear of an isotropic elastic and isotropic elastic-plastic material.     

Constitutive relations for isotropic or kinematic hardening at finite elastic–plastic deformations

The rate-type constitutive relations of rate-independent metals with isotropic or kinematic hardening at finite elastic–plastic deformations were presented through a phenomenological approach. This approach includes the decomposition of finite deformation into elastic and plastic parts, which is different from both the elastic–plastic additive decomposition of deformation rate and Lee’s elastic–plastic multiplicative decomposition of deformation gradient. The objectivity of the constitutive relations was dealt with in integrating the constitutive equations. A new objective derivative of back stress was proposed for kinematic hardening. In addition, the loading criteria were discussed. Finally, the stress for simple shear elastic–plastic deformation was worked out.     

A model for compression-weakening materials and the elastic fields due to contractile cells

We construct a homogeneous, nonlinear elastic constitutive law that models aspects of the mechanical behavior of inhomogeneous fibrin networks. Fibers in such networks buckle when in compression. We model this as a loss of stiffness in compression in the stress–strain relations of the homogeneous constitutive model. Problems that model a contracting biological cell in a finite matrix are solved. It is found that matrix displacements and stresses induced by cell contraction decay slower (with distance from the cell) in a compression weakening material than linear elasticity would predict. This points toward a mechanism for long-range cell mechanosensing. In contrast, an expanding cell would induce displacements that decay faster than in a linear elastic matrix.     

弹塑性有限变形的拟流动理论

本文提出一种弹塑性有限变形的拟流动理论。该理论从正交性法则出发，通过引入“拟弹性模量”和模量衰减函数并改进应变率的弹塑性分解，实现了由有限变形Ｐｒａｎｄｔｌ－Ｒｅｕｓｓ流动理论（Ｊ２Ｆ）向基于非正交法则的率形式形变理论（Ｊ２Ｄ）的合理的光滑过渡；并适用于初始及后继各向异性变形分析。在特殊条件下，可退化为Ｊ２Ｆ、Ｊ２Ｄ理论以及由任意各向异性屈服函数描述的流动理论。将该理论用于韧性金属平面应力／应变拉伸失稳与变形局部化的有限元模拟，并与理论分析及实验结果相比较，表明了本文理论的正确性。     

A unified thermodynamic theory of elasticity and plasticity

A thermodynamics is developed for a unified theory of elasticity and plasticity in infinitesmal strain. The constitutive equations which relate stress and strain deviators are rate type differential equations. When they satisfy a Lipschitz condition, uniqueness for the initial value problem dictates that the stress and strain will be related through elastic relations. Failure of the Lipschitz condition occurs when a von Mises yield condition is achieved: Plastic yield then occurs and the deviator relations turn into the Prandtl-Reuss equations. The plastic yield solution is stable during loading and unstable during unloading. The requirement that the solution followed during unloading be stable dictates entry into an elastic regime. Appropriate thermodynamic functions are constructed. It then appears that stress deviator (not strain deviator) is a viable state variable, and the thermodynamic relations are constructed in terms of a Gibbs function. The energy balance leads to satisfaction of the Clausius-Duhem inequality (and thus the second law of thermodynamics) in an elastic regime because it is shown that in an elastic regime entropy production is caused only by heat flux. During yield, the proper method of differentiating yields entropy production terms in addition to those arising from heat flux. These terms are positive during loading, whence it is concluded that the requirement that a stable solution be followed leads to satisfaction of the Clausius-Duhem inequality during plastic as well as elastic behavior.     

Towards the solution of finite plane-strain problems for compressible elastic solids

A new approach to the solution of finite plane-strain problems for compressible Isotropie elastic solids is considered. The general problem is formulated in terms of a pair of deformation invariants different from those normally used, enabling the components of (nominal) stress to be expressed in terms of four functions, two of which are rotations associated with the deformation. Moreover, the inverse constitutive law can be written in a simple form involving the same two rotations, and this allows the problem to be formulated in a dual fashion.For particular choices of strain-energy function of the elastic material solutions are found in which the governing differential equations partially decouple, and the theory is then illustrated by simple examples. It is also shown how this part of the analysis is related to the work of F. John on harmonic materials.Detailed consideration is given to the problem of a circular cylindrical annulus whose inner surface is fixed and whose outer surface is subjected to a circular shear stress. We note, in particular, that material circles concentric with the annulus and near its surface decrease in radius whatever the form of constitutive law within the given class. Whether the volume of the material constituting the annulus increases or decreases depends on the form of law and the magnitude of the applied shear stress.     

The effect of rotation and initial stress on the propagation of waves in a transversely isotropic elastic solid

In this paper the equations governing small amplitude motions in a rotating transversely isotropic initially stressed elastic solid are derived, both for compressible and incompressible linearly elastic materials. The equations are first applied to study the effects of initial stress and rotation on the speed of homogeneous plane waves propagating in a configuration with uniform initial stress. The general forms of the constitutive law, stresses and the elasticity tensor are derived within the finite deformation context and then summarized for the considered transversely isotropic material with initial stress in terms of invariants, following which they are specialized for linear elastic response and, for an incompressible material, to the case of plane strain, which involves considerable simplification. The equations for two-dimensional motions in the considered plane are then applied to the study of Rayleigh waves in a rotating half-space with the initial stress parallel to its boundary and the preferred direction of transverse isotropy either parallel to or normal to the boundary within the sagittal plane. The secular equation governing the wave speed is then derived for a general strain–energy function in the plane strain specialization, which involves only two material parameters. The results are illustrated graphically, first by showing how the wave speed depends on the material parameters and the rotation without specifying the constitutive law and, second, for a simple material model to highlight the effects of the rotation and initial stress on the surface wave speed.     

复合材料层合板粘弹性固化残余应力分析

基于Ｓｃｈａｐｅｒｙ积分型粘弹性本构关系，推导了考虑横向剪切效应的复合材料层合板线性热粘弹性有限元分析列式，对层合板的粘弹性响应和加工成型过程中的残余应力进行了分析，给出一些有意义的结果     

The Peierls stress in non-local elasticity

The effect of non-locality on the Peierls stress of a dislocation, predicted within the framework of the Peierls-Nabarro model, is investigated. Both the integral formulation of non-local elasticity and the gradient elasticity model are considered. A modification of the non-local kernel of the integral formulation is proposed and its effect on the dislocation core shape and size, and on the Peierls stress are discussed. The new kernel is longer ranged and physically meaningful, improving therefore upon the existing Gaussian-like non-locality kernels. As in the original Peierls-Nabarro model, lattice trapping cannot be captured in the purely continuum non-local formulation and therefore, a semi-discrete framework is used. The constitutive law of the elastic continuum and that of the glide plane are considered both local and non-local in separate models. The major effect is obtained upon rendering non-local the constitutive law of the continuum, while non-locality in the rebound force law of the glide plane has a marginal effect. The Peierls stress is seen to increase with increasing the intrinsic length scale of the non-local formulation, while the core size decreases accordingly. The solution becomes unstable at intrinsic length scales larger than a critical value. Modifications of the rebound force law entail significant changes in the core configuration and critical stress. The discussion provides insight into the issue of internal length scale selection in non-local elasticity models.     

A crystal plasticity model that incorporates stresses and strains due to slip gradients

This work is concerned with incorporating the kinematic and stress effects of excess dislocations in a constitutive model for the elastoplastic behavior of crystalline materials. The foundation of the model is a three term multiplicative decomposition of the deformation gradient in which the two classical terms of plastic and elastic deformation are included along with an additional term for long range strain due to the collective effects of excess dislocations. The long range strain is obtained from an assumed density of Volterra edge dislocations and is directly related to gradients in slip. A new material parameter emerges which is the size the region about a continuum point that contributes to long range strains.Using Hookean elasticity, the stress at a point is linearly related to the sum of the elastic plus the long range strain fields. However, the driving force for slip is postulated to be due only to the elastic stress so that the long range stress is a back stress in the constitutive relationship for plastic deformation. A consistent balance of the total deformation rate with the three proposed mechanisms of deformation leads to a set of differential equations that can be solved for the elastic stress, rotation and pressure which then implicitly defines the material state and equilibrium stress. Results from the simulation of a tapered tensile specimen demonstrate that the constitutive model exhibits isotropic and kinematic type hardening effects as well as changes in the pattern of plastic deformation and necking when compared to a material without slip gradient effects.