薄板冲压成型中板料起皱的临界应力分析和预测

针对薄板冲压成型中起皱这一常见的材料失效形式，运用板料压缩失稳理论，提出虚拟刚度变量的概念和板料产生起皱的临界应力计算方法，并借助于计算机仿真技术和有限元计算方法，搜寻并计算得到反映板料各处发生起皱难易程度的临界因子，生成起皱云图来观察材料各部分的稳定状态，预测可能出现起皱的部位，为修改冲压工艺和修模提供依据，通过分析起皱云图得出的结论与实际情况非常接近。     

On the wrinkling of a pre-stressed annular thin film in tension

Asymptotic properties of the neutral stability curves for a linear boundary eigenvalue problem which models the wrinkling instability of an annular thin film in tension are considered. The film is subjected to imposed radial displacement fields on its inner and outer boundaries and, when these loads are sufficiently large, the film is susceptible to wrinkling. The critical values at which this onset occurs are dictated by the solution of a fourth-order ordinary differential eigensystem whose eigenvalue λ is a function of μ(?1), a quantity inversely proportional to the non-dimensional bending stiffness of the film, and n, the number of half-waves of the wrinkling pattern that sets in around the annular domain. Previously, Coman and Haughton [2006. Localised wrinkling instabilities in radially stretched annular thin films. Acta Mech. 185, 179-200] employed the compound matrix method together with a WKB technique to characterise the form of λ(μ,n) which essentially is related to a turning point in a reduced eigenproblem. The asymptotic analysis conducted therein pertained to the case when this turning point was not too close to the inner edge of the annulus. However, in the thin film limit μ→∞, the wrinkling load and the preferred instability mode are given by a modified eigenvalue problem that involves a turning point asymptotically close to the inner rim. Here WKB and boundary-layer asymptotic methods are used to examine these issues and comparisons with direct numerical simulations made.     

Experiment and evaluation of wrinkling strain in a corner tensioned square membrane

Prediction of wrinkling characteristics is strongly correlated with the strain perpendicular to wrinkling direc- tion. In this paper, the strain field of wrinkled membrane is tested by VIC-3D system based on the digital image correlation technique. Experimental results are validated by the tension wrinkling simulation. The experimental strain perpendicular to wrinkling direction is analyzed in depth. The wrinkling strain of a square wrinkled membrane under corner tension is extracted from experimental strain perpendicular to wrinkling direction. A quantitative characterization format of the experimental wrinkling strain is proposed. A modified prediction method of wrinkling amplitude is presented based on the experimental wrinkling strain. The re- sults show that the precision of modified prediction model has improved 13.2% compared with the classical prediction model. The results reveal that the modified model can give an accurate prediction of the wrinkling amplitude.     

Path-dependence of the forming limit stresses in a sheet metal

The effect of changing strain paths on the forming limit stresses of sheet metals is investigated using the Marciniak–Kuczyński model and a phenomenological plasticity model with non-normality effects [Kuroda, M., Tvergaard, V., 2001. A phenomenological plasticity model with non-normality effects representing observations in crystal plasticity. J. Mech. Phys. Solids 49, 1239–1263]. Forming limits are simulated for linear stress paths and two types of combined loading: a combined loading consisting of two linear stress paths in which unloading is included between the first and second loadings (combined loading A), and combined loading in which the strain path is abruptly changed without unloading (combined loading B). The forming limit stresses calculated for combined loading A agree well with those calculated for the linear stress paths, while the forming limit curves in strain space depend strongly on the strain paths. The forming limit stresses calculated for the combined loading B do not, however, coincide with those calculated for the linear stress paths. The strain-path dependence of the forming limit stress is discussed in detail by observing the strain localization process.     

Anisotropic finite elastoplasticity with nonlinear kinematic and isotropic hardening and application to sheet metal forming

The paper discusses the derivation and the numerical implementation of a finite strain material model for plastic anisotropy and nonlinear kinematic and isotropic hardening. The model is derived from a thermodynamic framework and is based on the multiplicative split of the deformation gradient in the context of hyperelasticity. The kinematic hardening component represents a continuum extension of the classical rheological model of Armstrong–Frederick kinematic hardening. Introducing the so-called structure tensors as additional tensor-valued arguments, plastic anisotropy can be modelled by representing the yield surface and the plastic flow rule as functions of the structure tensors. The evolution equations are integrated by a new form of the exponential map that preserves plastic incompressibility and uses the spectral decomposition to evaluate the exponential tensor functions in closed form. Finally, the applicability of the model is demonstrated by means of simulations of several deep drawing processes and comparisons with experiments.     

Surface wrinkling of mucosa induced by volumetric growth: Theory, simulation and experiment

Mechanics of living tissues focusing on the relationships between growth, morphology and function is not only of theoretical interest but can also be useful for diagnosis of certain diseases. In this paper, we model the surface wrinkling morphology of mucosa, the moist tissue that commonly lines organs and cavities throughout the body, induced by either physiological or pathological volumetric growth. A theoretical framework of finite deformation is adopted to analyze the deformation of a cylindrical cavity covered by mucosal and submucosal layers. It is shown that compressive residual stresses induced by the confined growth of mucosa can destabilize the tissue into various surface wrinkling patterns. A linear stability analysis of the critical condition and characteristic buckling patterns indicates that the wrinkling mode is sensitive to the thicknesses of the mucosal and submucosal layers, as well as the properties of the tissues. The thinner the mucosal layer and the lower its elastic modulus, the shorter the buckling wavelength. A series of finite element simulations are performed to validate the theoretical predictions and to study local wrinkling or non-uniform patterns associated with inhomogeneous growth. Our postbuckling analysis shows that the surface pattern may evolve towards a period-doubling morphology due to continuous growth of mucosa or submucosa beyond the critical state. Finally, the theoretical predictions and numerical simulations are compared to experimental observations.     

Deformation localization and shear band fracture in strong anisotropy sheet tension

The tensile deformation localization and the shear band fracture behaviors of sheet metals with strong anisotropy are numerically simulated by using Updating Lagrange finite element method, Quasi-flow plastic constitutive theory^[1] and B-L planar anisotropy yield criterion^[2]. Simulated results are compared with experimental ones. Very good consistence is obtained between numerical and experimental results. The relationship between the anisotropy coefficientR and the shear band angle θ is found. The project supported by the National Natural Science Foundation of China and the Excellent Youth Teacher Foundation of the State Education Commission of China     

Marciniak–Kuczynski type modelling of the effect of Through-Thickness Shear on the forming limits of sheet metal

The Marciniak–Kuczynski (MK) forming limit model is extended in order to predict localized necking in sheet metal forming operations in which Through-Thickness Shear (TTS), also known as out-of-plane shear, occurs. An example of such a forming operation is Single Point Incremental Forming. The Forming Limit Diagram (FLD) of a purely plastic, isotropic hardening material with von Mises yield locus is discussed, for monotonic deformation paths that include TTS. If TTS is present in the plane containing the critical groove direction in the MK model, it is seen that formability is increased for all in-plane strain modes, except equibiaxial stretching. The increase in formability due to TTS is explained through a detailed study of some selected deformation modes. The underlying mechanism is a change of the stress mode in the groove that results in a delay of the onset of localized necking.     

Failure in notched tension bars due to high-temperature creep: Interaction between nucleation controlled cavity growth and continuum cavity growth

For axi-symmetrically notched tension bars [Dyson, B.F., Loveday, M.S., 1981, Creep Fracture in Nimonic 80A under Tri-axial Tensile Stressing, In: Ponter A.R.S., Hayhurst, D.R. (Eds.), Creep in Structures, Springer-Verlag, Berlin, pp. 406–421] show two types of damage propagation are shown: for low stress, failure propagates from the outside notch surface to the centre-line; and for high stress, failure propagates from the centre-line to the outside notch surface. The objectives of the paper are to: identify the physics of the processes controlling global failure modes; and, describe the global behaviour using physics-based constitutive equations.Two sets of constitutive equations are used to model the softening which takes place in tertiary creep of Nimonic 80A at 750 °C. Softening by multiplication of mobile dislocations is firstly combined, for low stress, with softening due to nucleation controlled creep constrained cavity growth; and secondly combined, for high stress, with softening due to continuum void growth. The Continuum Damage Mechanics, CDM, Finite Element Solver DAMAGE XX has been used to study notch creep fracture. Low stress notch behaviour is accurately predicted provided that the constitutive equations take account of the effect of stress level on creep ductility. High stress notch behaviour is accurately predicted from a normalized inverse cavity spacing d/2? = 6, and an initial normalized cavity radius r_hi/? = 3.16 × 10^?3, where 2? is the cavity spacing, and d is the grain size; however, the constants in the strain rate equation required recalibration against high stress notch data. A void nucleation mechanism is postulated for high stress behaviour which involves decohesion where slip bands intersect second phase grain boundary particles. Both equation sets accurately predict experimentally observed global failure modes.     

多元固相FAE的离散与爆轰的协调研究

一次起爆技术使FAE的两个主要研究方向,即燃料选择与云爆技术,更加紧密地结合在一起,导致对燃料的离散与爆轰的自协调要求更加严格。根据燃料的离散与爆轰的自协调的各制约因素对燃料组分进行选择,再改变燃料组分的质量配比以反馈研究燃料的离散与爆轰的协调性。实验证明所选择的混合燃料的离散与爆轰的协调性决定了超压峰值、正压作用时间、火球覆盖体积等参数,故可根据实际需要通过改变燃料组分配比来调节其协调性以研制不同特点的一次起爆型FAE。另外,基于VB.net平台自编的图像批处理程序,对由高速摄影拍摄到的FAE离散爆轰过程的图像进行处理,效率高、精度高,而且能和Excel、Origin等数据处理软件对接。     

Initiation and propagation of localized deformation in elasto-plastic strips under uniaxial tension

This paper deals with the evolution of inhomogeneous deformation in shape memory alloy strips and mild steel strips under uniaxial tension. New experiments on NiTi strips, which initially are in an austenitic phase, show that at a critical stress level martensite nucleates in sharp bands inclined at 55° to the axis of loading. Under prescribed end displacement martensite subsequently spreads either by steady-state propagation of inclined transition fronts or via a criss-cross pattern of finger-like features. Similar events have been reported in the literature regarding the evolution of Lüders bands in fine grained steel strips and wires. The similarity of macroscopic events, despite the different mechanisms of instability at the micro-level, prompted us to approximate the material behavior as a finitely deforming elasto-plastic solid with a trilinear up-down-up nominal stress-strain response. Two such stress-strain responses were used in finite element simulations of strip tension tests. In the first the true stress-strain response maintains its stability and in the second the intermediate branch has a negative slope. While both material models produced inhomogeneous deformations with features similar to those of the experiments, the larger initiation peak associated with the second gave results which closely resembled specific experiments. The numerical simulations confirmed that the evolution of events seen in experiments on SMAs and mild steels is strongly influenced by overall geometric (structural) effects. Furthermore, the success of this simple continuum constitutive model strongly suggests that continuum level events remain dominant players in such fine grained materials.     

Identification and characterization of delamination in polymer coated metal sheet

In this paper, a mixed numerical-experimental approach is adopted to quantitatively investigate and characterize delamination in polymer coated steel. The integral bi-material system is analysed and special attention is given to the constitutive modelling of the polymer coating, the interface and the determination of all involved parameters. An extended cohesive zone model for large displacements is proposed, allowing for a mode-dependent behaviour in large deformations. Finally, peel tests are used to characterize the interface, whereby the interfacial properties are determined through an inverse parameter identification procedure.     

Response of moving load due to irregularity in slightly compressible,finitely deformed elastic media

The present paper has been framed to study the stresses produced on the rough surface of a slightly compressible, finitely deformed half space due to a normal moving load. The surface of the medium is irregular with parabolic type of irregularity. The perturbation method is applied to find the displacement field. The normal and shear stresses have been obtained in closed form and discussed numerically by means of figures. It has been observed that the shear stress developed at different depths below the surface depends on the irregularity depth, frictional coefficient and irregularity factor of the rough surface of the medium. Also, surface plots have been drawn to analyze the combined variation of non-dimensional stresses and irregularity factor against depth.     

Micro/macro-crack growth due to creep–fatigue dependency on time–temperature material behavior

The implicit character of micro-structural degradation is determined by specifying the time history of crack growth caused by creep–fatigue interaction at high temperature. A dual scale micro/macro-equivalent crack growth model is used to illustrate the underlying principle of multiscaling which can be applied equally well to nano/micro. A series of dual scale models can be connected to formulate triple or quadruple scale models. Temperature and time-dependent thermo-mechanical material properties are developed to dictate the design time history of creep–fatigue cracking that can serve as the master curve for health monitoring.In contrast to the conventional procedure of problem/solution approach by specifying the time- and temperature-dependent material properties as a priori, the desired solution is then defined for a class of anticipated loadings. A scheme for matching the loading history with the damage evolution is then obtained. The results depend on the initial crack size and the extent of creep in proportion to fatigue damage. The path dependent nature of damage is demonstrated by showing the range of the pertinent parameters that control the final destruction of the material. A possible scenario of 20 yr of life span for the 38Cr2Mo2VA ultra-high strength steel is used to develop the evolution of the micro-structural degradation. Three micro/macro-parameters μ^*, d^* and σ^* are used to exhibit the time-dependent variation of the material, geometry and load effects. They are necessary to reflect the scale transitory behavior of creep–fatigue damage. Once the algorithm is developed, the material can be tailor made to match the behavior. That is a different life span of the same material would alter the time behavior of μ^*, d^* and σ^* and hence the micro-structural degradation history. The one-to-one correspondence of the material micro-structure degradation history with that of damage by cracking is the essence of path dependency. Numerical results and graphs are obtained to demonstrate how the inherently implicit material micro-structure parameters can be evaluated from the uniaxial bulk material properties at the macroscopic scale.The combined behavior of creep and fatigue can be exhibited by specifying the parameter ξ with reference to the initial defect size a₀. Large ξ (0.90 and 0.85) gives critical crack size a_cr = 11–14 mm (at t < 20 yr) for a₀ about 1.3 mm. For small ξ (0.05 and 0.15), there results critical a_cr = 6–7 mm (at t < 20 yr) for a₀ about 0.7–0.8 mm. The initial crack is estimated to increase its length by an order of magnitude before triggering global to the instability. This also applies ξ ≈ 0.5 where creep interacts severely with fatigue. Fine tuning of a_cr and a₀ can be made to meet the condition oft = 20 yr.Trade off among load, material and geometric parameters are quantified such that the optimum conditions can be determined for the desired life qualified by the initial–final defect sizes. The scenario assumed in this work is indicative of the capability of the methodology. The initial–final defect sizes can be varied by re-designing the time–temperature material specifications. To reiterate, the uniqueness of solution requires the end result to match with the initial conditions for a given problem. This basic requirement has been accomplished by the dual scale micro/macro-crack growth model for creep and fatigue.     

Flow of a micropolar fluid due to a porous stretching sheet and heat transfer

The present work investigates the micropolar fluid flow due to a permeable stretching sheet and the resulting heat transfer. Unlike the existing numerical works on the flow phenomenon in the literature, the prime interest here is to analytically work out shape of the solutions and identify whether they are unique. Indeed, unique solutions are detected and presented in the exact formulas for the associated boundary layer equations. Temperature field influenced by the microrotation is also mathematically resolved in the cases of constant wall temperature, constant heat flux and Newtonian heating. To discover the salient physical features of many mechanisms acting on the considered problem, it is adequate to have the analytical velocity and temperature fields and also closed-form skin friction/couple stress/heat transfer coefficients, all as given in the current paper. For instance, the practically significant rate of heat transfer is represented by a single formula valid for all three temperature cases.     

Analytical and numerical solutions for axisymmetric flow of nanofluid due to non-linearly stretching sheet

This article reports the laminar axisymmetric flow of nanofluid over a non-linearly stretching sheet. The model used for nanofluid contains the simultaneous effects of Brownian motion and thermophoretic diffusion of nanoparticles. The recently proposed boundary condition is considered which requires the mass flux of nanoparticles at the wall to be zero. Analytic solutions of the arising boundary value problem are obtained by optimal homotopy analysis method. Moreover the numerical solutions are computed by Keller–Box method. Both the solutions are found in excellent agreement. The behavior of Brownian motion on the fluid temperature and wall heat transfer rate is insignificant. Further the nanoparticle volume fraction distribution is found to be negative near the vicinity of the stretching sheet.     

A constitutive model dealing with damage due to cavity growth and the Mullins effect in rubber-like materials under triaxial loading

In this work, we attempted to describe the evolution of damage in rubber-like materials due to the Mullins effect and the cavity growth process. To this end we introduced two distinct internal variables into the constitutive laws; the first one essentially describes the Mullins damage and the second describes the cavity growth. The Mullins effect was considered as a continuous type of damage that can be modelled within the continuum damage theory. The cavity growth, being discontinuous at the microscopic scale, was also modelled by a continuous variable after a homogenization procedure. These analyses allow the establishment of a compressible constitutive law describing the strain-softening phenomena for rubber-like materials. In order to identify the material parameters and to verify the efficiency of the model, we carried out experimental studies involving uniaxial, biaxial, and hydrostatic tensions under monotonic and cyclic loading. Comparison between the model-predicted results and the experimental data shows that the present model can efficiently describe both the Mullins damage and the porosity evolution of rubber-like materials under triaxial monotonic or cyclic loading with a satisfactory accuracy. The proposed concept is simple and easy to apply to engineering calculations.     

Orthotropic yield criteria for description of the anisotropy in tension and compression of sheet metals

In this paper, yield functions describing the anisotropic behavior of textured metals are proposed. These yield functions are extensions to orthotropy of the isotropic yield function proposed by Cazacu et al. (Cazacu, O., Plunkett, B., Barlat, F., 2006. Orthotropic yield criterion for hexagonal close packed metals. Int. J. Plasticity 22, 1171–1194). Anisotropy is introduced using linear transformations of the stress deviator. It is shown that the proposed anisotropic yield functions represent with great accuracy both the tensile and compressive anisotropy in yield stresses and r-values of materials with hcp crystal structure and of metal sheets with cubic crystal structure. Furthermore, it is demonstrated that the proposed formulations can describe very accurately the anisotropic behavior of metal sheets whose tensile and compressive stresses are equal.     

A three-dimensional model of anisotropic hardening in metals and its application to the analysis of sheet metal formability

The hardening model proposed by Z. Mróz based on the uniaxial fatigue behavior of many metals is adopted to derive an incremental constitutive equation for general three-dimensional problems. This constitutive law is then employed in the analysis of metal forming problems to assess the influence of loading cycles, of the types involved in standard forming processes, on the ultimate formability of sheet metals. The predicted forming limit curves differ quantitatively from results obtained via an isotropie hardening model and differ qualitatively from those obtained via a kinematic model. Also investigated are the effects of such loading cycles on material response to simple tensile loading, which is often used to characterize a material. Significant differences between the present model and the other two models considered are observed in such characterizers of simple tensile behavior as the stress-strain curve, the anisotropy parameter and the uniform elongation. These differences suggest a rather simple experiment to identify the proper material model to be used in analyses of problems which involve loading cycles. Comparisons with some experimental results reveal that the employment of an anisotropic hardening model, such as the generalized Mróz model derived herein, is indeed crucial in accurately predicting material response to complicated loading histories.