有初始几何缺陷陷的壳体弹塑性大变形分析

本文在等温小变形弹塑性内时本构方程偏量形式的基础上，导出了适用于大位移、大动、小应变分析的弹塑性内时本构方程，进一步推导出了带有初始几何缺陷的非线性弹塑性问题的有限元方程，可用于分析缺陷对结构非线性塑性反应的影响，也可用于带缺陷的非线性问题求解及稳定性分析。     

有初始几何缺陷的壳体弹塑性大变形分析

在等温小变形弹塑性内时本构方程偏量形式的基础上，导出了适用于大位移、小应变分析的弹塑性内时本构方程。并导出了带有初始几何缺陷的非线性弹塑性问题的有限元方程。文中给出的算例表明本方法是可行有效的。     

内蕴时间塑性理论及其新进展(续)

三、内蕴时间塑性理论的一些基本概念和本构方程本节首先分别从经典塑性理论与内时理论的基本观点出发,对Drucker在其奠基性论文[46]中用过的一维力学模型的塑性响应特性进行分析比较,以便用简洁的方式阐述内时理论的某些最基本的概念,并说明经典模型可作为内时模型的一种理想化情况而得到。接着评述广义时间在固体力学中引入的概况和笔者最近提出的耗散型材料本构方程的形式不变性定律。然后根据这一定律首先得到了Valanis(1971)提出的小变形小变温下内时弹塑性本构方程的显式。接着研究了塑性应变偏张量的欧几里得模作为内时测度定义时的内时本构方程的特性,特别是讨论了它与经典塑性理论之间的关系。然后引入了含弱奇异性的内时本构方程。最后讨论了Valanis与笔者最近提出的新型弹塑性内时本构方程。      

内时理论在板料成形失稳中的应用

应用Valanis提出的内时本构方程，研究了板料成形的拉伸失稳问题，推导出单向和双向拉伸应力状态下的内时本构方程，据此分析了分散性失稳和集中性失稳。该文推导出应用于拉伸失稳分析时内时理论的近似表达式，它对应于经典塑性理论解，同时给出了内时理论的完整迭代数值解。结果表明内时理论具有很好的适用性。     

内时本构方程材料常数简易确定方法

提出确定内时塑性本构方程材料常数的一种简便直观的方法．在获得材料单向拉伸曲线之后,使用ＭｉｃｒｏｓｏｆｔＥｘｃｅｌ的规划求解功能快速确定弹性模量和内时本构常数,取得了令人满意结果,对推进内时本构理论在工程界的应用具有重要意义．     

弹塑性屈曲问题的一种新的分析方法

本文利用内时本构模型提出了分析计算弹塑性压杆屈曲问题的一种新方法,建立了适应于整个长细比范围内的稳定屈曲统一公式.文中方法适应性强,分析计算过程明了.对铅合金柱的分析计算表明,利用文中分析方法可以得到有效合理的较为精确的屈曲结果.     

钢方管截面柱考虑损伤的滞回性能分析

结构钢本构关系的精度直接影响分析结果的可靠度.根据能量等效性假设、热力学第二定律,推出损伤材料的弹性本构方程;采用混合强化准则,考虑Bauschinger效应、屈服平台、硬化(软化)效应及损伤和损伤演化影响,建立了的结构钢弹塑性各向异性损伤本构关系.结合构建的本构关系,采用八节点超参数壳体单元,推导了用U.L.格式及Cauchy应力描述的板壳双重非线性有限元方程,并编制了计算程序.利用U.L.格式的壳体大挠度双重非线性有限元分析方法,对钢方管截面短柱进行面内拉压循环荷载作用下的滞回性能分析.     

裂尖大应变细观断裂研究

本文用反映空穴形核成长的Gurson本构方程来描述裂尖区域材料在大应变情形下的力学特性，并进一步考虑了空穴演变对材料杨氏模量的影响。文中用上述本构方程分别结合弹塑性大应变有限元方法对平面应变I型裂纹问题作了计算，分析了裂尖应力分布、裂尖形状变化和裂尖空穴演变过程，并与用Prandtl-Reuss本构方程教育处的结果作了比较。     

大变形中摩擦接触问题的数值模拟及应用

大变形的摩擦接触是复杂的非线性问题 ,本文介绍了一种处理摩擦接触问题的数值方法。采用接触单元技术模拟接触界面 ,基于弹塑性理论形式的非经典Coulomb摩擦定律及罚函数方法建立了摩擦接触的增量本构关系。结合大变形的增量分析格式给出了积分摩擦接触本构方程的回映方法。这种处理摩擦接触问题的方法计算简单、使用方便。给出的计算实例及应用实例说明了方法的精度与稳定性     

一种新的岩石非线性黏弹塑性蠕变模型研究

为了克服传统元件组合模型不能描述岩石蠕变过程中非线性特征的缺陷，首先根据加速蠕变阶段的应变和应变率随蠕变时间急剧增大的特点，建立黏塑性应变与蠕变时间的指数函数关系并提出非线性黏塑性体．将该非线性黏塑性体与广义Burgers蠕变模型串联，建立可以描述岩石全蠕变过程的非线性黏弹塑性蠕变模型，根据叠加原理得到一维应力状态下的轴向蠕变方程．然后基于塑性力学理论指出岩石三维蠕变本构方程建立过程中的不足之处，并给出非线性黏弹塑性蠕变模型合理的三维蠕变方程．最后采用不同应力水平下砂岩轴向蠕变试验对模型合理性进行验证，结果表明：拟合曲线与试验曲线吻合度较高，所建蠕变模型能够很好地描述砂岩在不同应力水平下的蠕变变形规律，尤其对加速蠕变阶段的非线性特征描述效果很好，验证了模型的合理性．     

Constitutive equations for describing the elastoplastic deformation of elements of a body along small-curvature paths in view of the stress mode

Equations relating the components of the stress and strain tensors (constitutive equations) are formulated in terms of Euler coordinates. The equations describe the finite elastoplastic deformation of an isotropic body along paths of small curvature. It is assumed that the stress deviator is coaxial with the plastic-strain differential deviator. The relationships between the first and second invariants of the stress and strain tensors in the case of complex elastoplastic deformation of the body’s elements are determined from base tests on tubular specimens loaded along rectilinear paths for several values of the stress mode angle. Methods for specification of these relationships are proposed. The assumptions adopted to derive the constitutive equations are validated experimentally __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 4, pp. 62–72, April 2006.     

Axisymmetric Thermoelastoplastic Stress State of Isotropic Solids of Revolution Under Impulsive Loading

Consideration is given to the solution of a dynamic problem for a solid of revolution with an arbitrary meridional section under impulsive thermomechanical loading inducing elastoplastic strains. The theory of small elastoplastic deformations is used. The constitutive equations are linearized by the variable-parameter method. The unloading process is described by a linear law. The solution technique involves the finite-element approximation in spatial coordinates and the finite-difference representation of time derivatives. Based on the principle of linear summation, recurrent relations are derived for successive evaluation of nodal displacements by an explicit scheme in time. Solution for cylinders and disks are presented to illustrate the influence of elastoplastic deformations on wave processes     

An elastoplastic framework for granular materials becoming cohesive through mechanical densification. Part II – the formulation of elastoplastic coupling at large strain

The two key phenomena occurring in the process of ceramic powder compaction are the progressive gain in cohesion and the increase of elastic stiffness, both related to the development of plastic deformation. The latter effect is an example of ‘elastoplastic coupling’, in which the plastic flow affects the elastic properties of the material, and has been so far considered only within the framework of small strain assumption (mainly to describe elastic degradation in rock-like materials), so that it remains completely unexplored for large strain. Therefore, a new finite strain generalization of elastoplastic coupling theory is given to describe the mechanical behaviour of materials evolving from a granular to a dense state.The correct account of elastoplastic coupling and of the specific characteristics of materials evolving from a loose to a dense state (for instance, nonlinear – or linear – dependence of the elastic part of the deformation on the forming pressure in the granular – or dense – state) makes the use of existing large strain formulations awkward, if even possible. Therefore, first, we have resorted to a very general setting allowing general transformations between work-conjugate stress and strain measures; second, we have introduced the multiplicative decomposition of the deformation gradient and, third, employing isotropy and hyperelasticity of elastic response, we have obtained a relation between the Biot stress and its ‘total’ and ‘plastic’ work-conjugate strain measure. This is a key result, since it allows an immediate achievement of the rate elastoplastic constitutive equations. Knowing the general form of these equations, all the specific laws governing the behaviour of ceramic powders are finally introduced as generalizations of the small strain counterparts given in Part I of this paper.     

Unified stress update algorithms for the numerical simulation of large deformation elasto-plastic and elasto-viscoplastic processes

This paper is concerned with unified stress update algorithms for elastoplastic and elasto-viscoplastic constitutive equations for metals submitted to large deformations. We present here a newly developed time integration algorithm which is, in the case of J2 flow theory material behavior, an extension to the viscoplastic range of the classical radial return algorithm for plasticity. The resulting unified implicit algorithm is both efficient and very inexpensive. Moreover, if there is no viscosity effect (rate-independent material) the presented algorithm degenerates exactly into the classical radial return algorithm for plasticity.     

Development of the large increment method for elastic perfectly plastic analysis of plane frame structures under monotonic loading

The displacement-based finite element method dominates current practice for material nonlinear analysis of structures. However, there are several characteristics that may limit the effectiveness of this approach. In particular, for elastoplastic analysis, the displacement method relies upon a step-by-step incremental approach stemming from flow theory and also requires significant mesh refinement to resolve behavior in plastic zones. This leads to computational inefficiencies that, in turn, encourage the reconsideration of force-based approaches for elastoplastic problems.One of these force algorithms that has been recently developed is the large increment method. The main advantage of the flexibility-based large increment method (LIM) over the displacement method is that it separates the global equilibrium and compatibility equations from the local constitutive relations. Consequently, LIM can reach the solution in one large increment or in a few large steps, thus, avoiding the development of cumulative errors. This paper discusses the extension of the large increment methodology for the nonlinear analysis of plane frame structures controlled by an elastic, perfectly plastic material model. The discussion focuses on the power of LIM to handle these nonlinear problems, especially when plastic hinges form in the frame and ultimately as the structure approaches the collapse stage. Illustrative planar frame examples are presented and the results are compared with those obtained from a standard displacement method.     

Pressure dependency,strength-differential effect,and plastic volume expansion in metals

Within the framework of a purely mechanical rate-type theory of finitely deforming elastic-plastic materials, a special set of constitutive equations is introduced. These equations include a dependence both on mean normal stress and on plastic strain, and involve non-normality of plastic strain-rate. They accommodate both the strength-differential effect and plastic volume expansion. Values of the material coefficients in the constitutive equations are approximately determined with the help of published experimental data for AISI 4330 steel. The plastic strain-rate is shown to have a large non-normal component, resulting in a small plastic volume expansion. The theory is in reasonable agreement with experiment for plastic strains up to about 2% in magnitude. Further experiments are suggested so that the constitutive coefficients can be completely and accurately determined.     

各向异性本构关系在板料成形数值模拟中的应用

对几种能表达面内各向异性的屈服准则Hill、Barlat-Lian、Barlat进行了比较。以弹性变形服从各向同性广义虎克定律的情况下，给出了基于张量算法推导的弹塑性本构关系的一般表达式，并由此导出了相应屈服准则的弹塑性本构关系的显式表达。借助ABAQUS软件本构模块用户子程序接口，分别实现了这些屈服准则在ABAQUS的嵌入。以模拟方形盒的拉延过程为例，分析了不同的屈服准则在板料成形过程数值模拟中的应用。模拟结果表明，基于弹塑性本构关系一般表达所列出的相应屈服准则的显式表达式是正确的；在采用壳元来模拟板料成形时，采用Barlat准则的模拟结果和采用Barlat-Lian准则的结果差别不大。     

解Drucker-Prager塑性问题的二阶锥互补法

基于经典弹塑性理论中多数屈服准则具有凸锥数学结构的事实,将在大规模计算中更具潜力的锥规划法引入弹塑性分析。考虑到弹塑性流动理论有关联与非关联之分,本文提出利用锥型互补法求解弹塑性问题。具体以Drucker-Prager弹塑性模型为例,首先利用最大塑性功耗散原理和圆锥对偶理论等工具,建立了弹塑性本构方程的等价二阶锥互补模型;然后,基于参变量变分原理和有限元技术,建立了弹塑性增量分析的二阶锥线性互补模型;最后,利用一类半光滑Newton算法求解。数值算例表明了本文方法的有效性。     

Unified elastoplastic finite difference and its application

Two elastoplastic constitutive models based on the unified strength theory (UST) are established and implemented in an explicit finite difference code, fast Lagrangian analysis of continua (FLAC/FLAC3D), which includes an associated/non-associated flow rule, strain-hardening/softening, and solutions of singularities. Those two constitutive models are appropriate for metallic and strength-different (SD) materials, respectively. Two verification examples are used to compare the computation results and test data using the two-dimensional finite difference code FLAC and the finite element code ANSYS, and the two constitutive models proposed in this paper are verified. Two application examples, the large deformation of a prismatic bar and the strain-softening behavior of soft rock under a complex stress state, are analyzed using the three-dimensional code FLAC3D. The two new elastoplastic constitutive models proposed in this paper can be used in bearing capacity evaluation or stability analysis of structures built of metallic or SD materials. The effect of the intermediate principal stress on metallic or SD material structures under complex stress states, including large deformation, three-dimensional and non-association problems, can be analyzed easily using the two constitutive models proposed in this paper.