晶粒数量对多晶集合体初始各向异性的影响

Taylor类多晶晶体粘塑性模型被用于研究晶粒数量对随机分布多晶体拉伸塑性各向异性的影响。分别沿包含不同晶粒数量的多晶集合体的各方向进行单向拉伸数值模拟实验,得到多晶集合体各方向在一定等效应变下的等效应力,并用云图和等高线表示在多晶体的参考球面上。定义了描述多晶集合体各向异性程度的参考指标。讨论了三种确定晶体随机取向的方法。计算结果表明：晶粒数量有限的多晶集合体的应力应变响应仍有一定的各向异性,且随着晶粒数量增多,多晶集合体的各向异性程度降低;就所包含晶粒数相同的多晶集合体来说,在确定晶粒随机取向时,选取不同的方法对它的各向异性程度也有一定的影响。     

Influence of domain switching zones on the fracture toughness of ferroelectrics

The paper presents a fracture model for ferroelectric materials taking into account the hysteretic domain switching processes near to the tip of a macroscopic crack. The model is based on the balance of energy supplied by the driving forces, on the one hand, and the total of energies either dissipated by domain switching, stored in the crack wake region or consumed by the formation of new fracture surface, on the other hand. An internal variable theory describes the nonlinear coupled electromechanical material response within the framework of a three-dimensional continuum model. For simplicity, the complex orientation distribution function of domains in a polycrystalline ceramic is approximated by only six representative space orientations. The theory predicts certain dimensionless material parameter combinations which govern the change of fracture toughness under the application of different mechanical and electrical loadings. A comparison with data available in the literature for barium titanate ceramics yields a reasonable coincidence.     

Fast calculation of average Taylor factors and Mandel spins for all possible strain modes

The average Taylor factor of the crystallites of a polycrystalline material and the average spins (‘spin’ here means rate of rotation) of the crystal lattices (Mandel spin) are functions of the strain mode to which the material is subjected. These functions can be very helpful when constructing material models to be used by elastic–plastic FE simulations of forming processes in which the anisotropy due to texture, back stress etc. is taken into account. These functions can be calculated from the crystallographic texture. It is shown in the present paper, that a calculation done for all possible strain modes can be organised in a very efficient way. This is made possible by the Fourier series expansion method used in quantitative texture analysis. The newest calculation method which uses this methodology is 400 times faster than one which would not, and this with sufficient accuracy. This makes it possible to incorporate calculations of complete yield loci in FE simulations without prohibitive increase in computer time.     

Evaluation of identification methods for YLD2004-18p

Four calibration methods have been evaluated for the linear transformation-based anisotropic yield function YLD2004-18p (Barlat, F., Aretz, H., Yoon, J.W., Karabin, M.E., Brem, J.C., Dick, R.E., 2005. Linear transformation-based anisotropic yield functions. Int. J. Plasticity 21, 1009–1039) and the aluminium alloy AA5083-H116. The different parameter identifications are based on least squares fits to combinations of uniaxial tensile tests in seven directions with respect to the rolling direction, compression (upsetting) tests in the normal direction and stress states found using the full-constraint (FC) Taylor model for 690 evenly distributed strain paths. An elastic–plastic constitutive model based on YLD2004-18p has been implemented in a non-linear finite element code and used in finite element simulations of plane-strain tension tests, shear tests and upsetting tests. The experimental results as well as the Taylor model predictions can be satisfactorily reproduced by the considered yield function. However, the lacking ability of the Taylor model to quantitatively reproduce the experiments calls for more advanced crystal plasticity models.     

Non-linear vibration of rectangular reticulated shallow shell structures

ＮＯＮ－ＬＩＮＥＡＲＶＩＢＲＡＴＩＯＮＯＦＲＥＣＴＡＮＧＵＬＡＲＲＥＴＩＣＵＬＡＴＥＤＳＨＡＬＬＯＷＳＨＥＬＬＳＴＲＵＣＴＵＲＥＳＮｉｅＧｕｏ－ｈｕａ（聂国华）（ＤｅｐａｒｔｍｅｎｔｏｆＥｎｇｉｎｅｅｒｉｎｇＭｅｃｈａｎｉｃｓ,ＴｏｎｇｊｉＵｎｉｖｅ...     

NUMERICAL SIMULATION OF NON-PLANAR PLASTIC ANISOTROPY OF COLD ROLLING POLYCRYSTALLINE MATERIALS

The plastic anisotropy of sheet metal is usually caused by preferred orientation of grains, developed by mechanical deformation and thermal treatment. In the present study, a Taylor-like polycrystal model suggested by Asaro and Needleman is applied to investigate the evolution of the anisotropic behavior of a face centered cubic (FCC) polycrystalline metal, which is considered having {111} (110) slip systems, by stretching it along an arbitrary direction after it has undergonea plane-strata compression that rationally simulates the cold rolling process of FCC polycrystalline pure aluminium. By using the Taylor-like polycrystal model, pole figures are obtained to describe the texture development of polycrystalline aggregate after plane-strain compression, and then the plastic anisotropy of polycrystalline aggregate is evaluated by stretching the polycrystalline aggregate in different direction in term of yield stress. According to the results, the contours of longitudinal flow stress in three-dimensional orientation space are given and analyzed. Experiment results similar to the prediction of planar anisotropy can be found inthe literature written by Takahashi et al. that in directly show the correctness of the prediction of non-planar plastic anisotropy by this analysis.     

On the r-value of textured sheet metals

The r-value of a sheet metal is a measure of plastic anisotropy frequently used for prediction of performance in deep-drawing. It has also figured prominently in the literature for validation of theories where the predicted angular dependence of r is compared with the measured dependence. As plastic anisotropy in sheet metals is caused mainly by the preferred orientations of grains within the polycrystalline metal, it is natural to ask how r would depend on the orientation distribution function (ODF) w which defines the crystallographic texture of the polycrystal. In this paper a general formula relating r to w is derived for textured sheet metals whose plastic flow behavior is governed by a plastic potential f(σ, w), the anisotropic part of which depends linearly on the texture coefficients; here σ denotes the deviator of the Cauchy stress. Specific forms of this formula for orthorhombic sheets of cubic and of hexagonal metals are explicitly given.     

OPERATOR–SPLITTING COMPUTATION OF TURBULENT FLOW IN AN AXISYMMETRIC 180° NARROWING BEND USING SEVERAL k—ϵ MODELS AND WALL FUNCTIONS

This paper describes a finite element implementation of an operator-splitting algorithm for solving transient/steady turbulent flows and presents solutions for the turbulent flow in an axisymmetric 180° narrowing bend, a benchmark problem dealt with at the 1994 WUA-CFD annual meeting. Three k–ϵ based models are used: the standard linear k–ϵ model, a non-linear k–ϵ model and an RNG k–ϵ model. Flow separation after the bend, as observed in the experiment, is predicted by the RNG model and by both the linear and non-linear k–ϵε models with van Driest mixing length wall functions. Good agreement with experimental data of pressure distribution on bending walls is obtained by the present numerical simulation. Results show that there is very little difference between the linear and non-linear k–ϵε models in terms of predicted velocity fields and that the non-linearities mainly affect the distribution of turbulent normal stress and pressure, in analogy to the effect of second-order viscoelastic fluid models on laminar flow. Both the linear and non-linear k–ϵε models fail to predict any flow separation if logarithmic wall functions are used.     

The evolution of crystalline stresses of a polycrystalline metal during cyclic loading

The presence of a positive average applied stress during cyclic uniaxial loading leads to a reduction in fatigue life of metallic parts. The metals are typically polycrystalline, with stresses varying from crystal to crystal due to differences in lattice orientation and slip system strength. Simulations enable us to better understand how polycrystals behave under cyclic loading and how the changing stress over many cycles influences fatigue life. Specifically, uniaxial cyclic simulations of pre-strained HY100 steel were conducted using an elastic viscoplastic continuum slip model employing a Taylor hypothesis. Stress-controlled loading conditions were employed to mimic fatigue tests on cold-bent bar specimens for three different load levels. The macroscopic axial strains and the crystal axial stresses were monitored during the cycles. The stress–strain response for the first cycle was used to determine the load input for the material point simulations. The peak values of crystal axial stress were found to evolve continuously with the number of loading cycles. It was found that the stress change in a crystal is influenced not only by its own orientation but also by the orientations of the other crystals in the aggregate. Furthermore, the distribution of crystal stresses after thousands of cycles at a lower stress amplitude closely resembled the distribution after tens of cycles at a larger stress amplitude.     

Self-similar problems of elastic contact for non-convex punches

Self-similar problems of contact for non-convex punches are considered. The non-convexity of the punch shapes introduces differences from the traditional self-similar contact problems when punch profiles are convex and their shapes are described by homogeneous functions. First, three-dimensional Hertz type contact problems are considered for non-convex punches whose shapes are described by parametric-homogeneous functions. Examples of such functions are numerous including both fractal Weierstrass type functions and smooth log-periodic sine functions. It is shown that the region of contact in the problems is discrete and the solutions obey a non-classical self-similar law. Then the solution to a particular case of the contact problem for an isotropic linear elastic half-space when the surface roughness is described by a log-periodic function, is studied numerically, i.e. the contact problem for rough punches is studied as a Hertz type contact problem without employing additional assumptions of the multi-asperity approach. To obtain the solution, the method of non-linear boundary integral equations is developed. The problem is solved only on the fundamental domain for the parameter of self-similarity because solutions for other values of the parameter can be obtained by renormalization of this solution. It is shown that the problem has some features of chaotic systems, namely the global character of the solution is independent of fine distinctions between parametric-homogeneous functions describing roughness, while the stress field of the problem is sensitive to small perturbations of the punch shape.     

Anti-plane shear of an arbitrary oriented crack in a functionally graded strip bonded with two dissimilar half-planes

An internal crack located within a functionally graded material (FGM) strip bonded with two dissimilar half-planes and under an anti-plane load is considered. The crack is oriented in an arbitrary direction. The material properties of strip are assumed to vary exponentially in the thickness direction and two half-planes are assumed to be isotropic. Governing differential equations are derived and to reduce the difficulty of the problem dealing with solution of a system of singular integral equations Fourier integral transform is employed. Semi closed form solution for the stress distribution in the medium is obtained and mode III stress intensity factor (SIF), at the crack tip is calculated and its validity was verified. Finally, the effects of nonhomogeneous material parameter and crack orientation on the stress intensity factor are studied.     

ORIENTATION DISTRIBUTION FUNCTIONS FOR MICROSTRUCTURES OF HETEROGENEOUS MATERIALS (Ⅰ) ─ DIRECTIONAL DISTRIBUTION FUNCTIONS AND IRREDUCIBLE TENSORS

In this two-part paper, a thorough investigation is made on Fourier expansions with irreducible tensorial coefficients for orientation distribution functions (ODFs) and crystal orientation distribution functions (CODFs), which are scalar functions defined on the unit sphere and the rotation group, respectively. Recently it has been becoming clearer and clearer that concepts of ODF and CODF play a dominant role in various micromechanically-based approaches to mechanical and physical properties of heterogeneous materials. The theory of group representations shows that a square integrable ODF can be expanded as an absolutely convergent Fourier series of spherical harmonics and these spherical harmonics can further be expressed in terms of irreducible tensors. The fundamental importance of such irreducible tensorial coefficients is that they characterize the macroscopic or overall effect of the orientation distribution of the size, shape, phase, position of the material constitutions and defects. In Part (Ⅰ), the investigation about the irreducible tensorial Fourier expansions of ODFs defined on the N-dimensional (N-D) unit sphere is carried out. Attention is particularly paid to constructing simple expressions for 2- and 3-D irreducible tensors of any orders in accordance with the convenience of arriving at their restricted forms imposed by various point-group (the synonym of subgroup of the full orthogonal group) symmetries. In the continued work -Part (Ⅱ), the explicit expression for the irreducible tensorial expansions of CODFs is established. The restricted forms of irreducible tensors and irreducible tensorial Fourier expansions of ODFs and CODFs imposed by various point-group symmetries are derived.     

Harmonic analysis in rheological property measurement

 Sudden change in the linear stress-strain relationship has long been known to produce interesting discontinuities in the flow properties of polymer melts, concentrated particle suspensions, and gels as the shear rate amplitude is increased. The non-linear effect produces stress harmonics, which contain information about the material properties that is otherwise difficult to obtain. This note describes a generalized procedure for extracting yield stress material functions from the discrete Fourier transform of the stress signal. Received: 23 April 2001 Accepted: 27 July 2001     

Reynolds’ dream?

For many non-linear phenomena, it is necessary to solve infinite systems of equations for correlation functions with a wide range of parameters. This paper can be seen as a first step in addressing this problem. Correlation functions related to the ?³ and to the ?⁴ field theories are described by means of generating Fock space vectors constructed with the help of Cuntz algebra. The equations obtained are easily transformed and general solutions are constructed. Various expansions of these solutions are developed using the explicitly constructed right inverse operators related to linear and non-linear parts of the theory. Based on the idea of information loss and using the language of classical mechanics, a solution to the closure problem for correlation functions is proposed. The method described in this paper can be used to obtain approximated correlation functions with strong and weak non-linearity.     

Limit analysis for a general class of yield conditions

The paper describes a generalisation of the programming method described by Ponter and Carter (1997) for the evaluation of optimal upper bounds on the limit load of a body composed of a rigid/perfectly plastic material. The method is based upon similar principles to the `Elastic Compensation' method which has been used for design calculations for some years but re-interpreted as a non-linear programming method. A sufficient condition for convergence is derived which relates properties of the yield surface to those of the linear solutions solved at each iteration. The method is demonstrated through an application to a Drucker–Prager yield condition in terms of the Von Mises effective stress and the hydrostatic pressure. Implementation is shown to be possible using the user routines in a commercial finite element code, ABAQUS. The examples chosen indicate that stable convergent solutions may be obtained. There are, however, limits to the application of the method if isotropic linear solutions are used for an isotropic yield surface. In an accompanying paper (Ponter and Engelhardt, 2000) the method is extended to shakedown and related problems.     

Solutions of some simple boundary value problems within the context of a new class of elastic materials

Some simple boundary value problems are studied, for a new class of elastic materials, wherein deformations are expressed as non-linear functions of the stresses. Problems involving ‘homogeneous’ stress distributions and one-dimensional stress distributions are considered. For such problems, deformations are calculated corresponding to the assumed stress distributions. In some of the situations, it is found that non-unique solutions are possible. Interestingly, non-monotonic response of the deformation is possible corresponding to monotonic increase in loading. For a subclass of models, the strain-stress relationship leads to a pronounced strain-gradient concentration domain in the body in that the strains increase tremendously with the stress for small range of the stress (or put differently, the gradient of the strain with respect to the stress is very large in a narrow domain), and they remain practically constant as the stress increases further. Most importantly, we find that for a large subclass of the models considered, the strain remains bounded as the stresses become arbitrarily large, an impossibility in the case of the classical linearized elastic model. This last result has relevance to important problems in which singularities in stresses develop, such as fracture mechanics and other problems involving the application of concentrated loads.     

弹塑性结构安定下限分析的无网格局部

将基于Voronoi结构的无网格局部Petrov-Galerkin法与减缩基技术相结合,建立了一种安定下限分析的新方法．为了克服移动最小二乘近似难以准确施加本质边界条件的缺点,采用了自然邻近插值构造试函数.通过引入基准载荷域上载荷角点的概念,消除了安定下限分析中由时间参数所引起的求解困难．利用减缩基技术,将安定分析问题化为一系列未知变量较少的非线性规划子问题．在每个非线性规划子问题中,自平衡应力场由一组带有待定系数的自平衡应力场基矢量的线性组合进行模拟,而这些自平衡应力场基矢量可应用弹塑性增量分析中的平衡迭代结果得到．算例结果证明了提出的分析方法的有效性．      

The Facet method: A hierarchical multilevel modelling scheme for anisotropic convex plastic potentials

An entirely new analytical expression describing plastic anisotropy is presented. It is designed to be used in combination with multilevel models. It makes use of the theory of dual plastic potentials, which is shortly revisited. An analysis of convexity is presented. Note that the new method is not optimal when not used in combination with a multilevel model; other methods are better suited for identification on the basis of mechanical test data. Compared with already existing methods which work with multilevel models, the new method has the following advantages: (i) it is automatically convex anywhere in the six-dimensional stress or strain rate space; (ii) it can be used for materials with a stress differential effect, such as hcp metals or pre-strained cubic materials; (iii) its identification procedure is such that not only the Taylor theory, but also more advanced theories, such as the Alamel-model or self-consistent models, can be used to identify the parameters; (iv) an analytical expression of the plastic potential can be obtained in both strain rate and stress space, which is an important advantage when implementing the model in finite element codes for metal forming. Equipotential surfaces in strain rate space and corresponding yield loci obtained by the new method for four materials (one ferrite single crystal, one aluminium alloy and two types of steel) are presented and discussed.     

基于晶体塑性理论研究铝材料高压高应变率下的强度特性

为了了解金属材料在极端加载下复杂动态响应过程中的多种机制和效应,重点针对Al材料在高压、高应变率加载下的塑性变形机制,在经典晶体塑性模型的基础上,对其中的非线性弹性、位错动力学和硬化形式进行改进,建立适用于高压、高应变率加载下的热弹-黏塑性晶体塑性模型。该模型可以较好地描述单晶铝和多晶铝材料屈服强度随压力的变化过程,相比宏观模型,用该模型还获得了多晶Al材料在冲击加载下的织构演化规律,揭示了织构择优取向行为和压力的关系。     

Fourier-transform rheology

Oscillatory shear of polymeric liquids in the non-linear regime generates higher harmonic contributions in the shear stress response. These non-linear contributions are analyzed in Fourier space with respect to the different frequencies and intensities. Simulated and experimental Fourier rheology spectra for atactic poly(propylene) melts are shown. Received: 28 January 1998 Accepted: 26 May 1998