Bammann-Chiesa-Johnson粘塑性本构模型材料参数的一种识别方法

粘塑性本构模型能否成功模拟金属高应变率大应变变形过程依赖于材料参数识别结果的好坏。由于BCJ模型考虑了应变率、温度与材料硬化之间的耦合效应以及应变率、温度历史效应,同时模型中包含了多个材料参数,因此很难通过试验直接识别模型的材料参数。本文针对BCJ模型中的耦合效应和历史效应,基于对模型中材料参数物理涵义的界定,给出了一种对材料参数解耦、分离并进行估计的方法,获得了模型材料参数估计公式,估计了材料参数的取值范围。在此基础上,编制了BCJ模型应力积分径向返回算法和粒子群优化算法的计算程序,应用重新设计了BCJ模型耦合效应和历史效应的反分析方法,在参数取值范围内对材料参数进行了优化识别。以OFHC Cu为例,应用提出的识别方法对BCJ模型的材料参数进行了识别,计算结果和试验结果符合较好。     

Constitutive modeling of the material behavior at elevated temperatures and related material testing

Constant strain rate tests, both in tension and in torsion, were performed on 2618-T61 aluminum alloy at 200°C. The results were compared with those from the creep tests. It was demonstrated that the macroscopically observed material behavior was essentially a consequence of the inherent time dependence of the inelastic strain. Based on the experimental results, the traditional approach and the unified approach for modeling the material behavior at elevated temperatures were compared. The material tests required for the development of the constitutive model and their limitations were also discussed.     

Some new extensions to multi-mechanism models for plastic and viscoplastic material behavior under small strains

Multi-mechanism models (MM models) have become an important tool for modeling complex material behavior. In particular, two-mechanism models are used. They are applied to model ratcheting in metal plasticity as well as steel behavior during phase transformations. We consider a small-deformation setting. The characteristic trait of multi-mechanism models is the additive decomposition of the inelastic (e.g., plastic or viscoplastic) strain into several parts. These parts are sometimes called mechanisms. In comparison with rheological models, the mechanisms can interact with each other. This leads to new properties and allows to describe important observable effects. Up to now, each mechanism has one kinematic internal variable. As a new feature, we develop multi-mechanism models (in series) with several kinematic variables for each mechanism as well as with several isotropic variables for each flow criterion. We describe this complex situation by three structural matrices which express the mutual relations between mechanisms, flow criteria, kinematic, and isotropic variables. The well-known Chaboche model with a unique inelastic strain and several kinematic variables represents a special case of these general multi-mechanism models. In this work, we also present a matrix-based approach for these new complex MM models. The presented models can form the basis for developing numerical algorithms for simulation and parameter identification.     

Thermo-viscoelastic and viscoplastic behavior of high-density polyethylene

Observations are reported on high-density polyethylene in uniaxial tensile tests with constant strain rates and relaxation tests at various temperatures ranging from 25 to 90 °C. A constitutive model is derived for the nonlinear viscoelastic and viscoplastic behavior of semi-crystalline polymers at three-dimensional deformations. Adjustable parameters in the stress–strain relations are found by fitting the experimental data. It is demonstrated that (i) the model correctly approximates the observations and (ii) material parameters are independent of strain rate and change consistently with temperature.     

Unsteady shear flow of a viscoplastic material

A viscoplastic, or yield-stress, liquid occupies the space between two infinite parallel plates. Initially the whole system is at rest. The lower plate is suddenly jerked into motion with given speed or shear stress, while the upper plate is kept fixed. The flow consists of two regions; (1) a lower sheared region bounded above by the yield surface, (2) an upper unyielded region bounded below by the yield surface. The yield surface propagates to the upper plate as time proceeds. We first consider the equivalent one plate problem of flow over a jerked plate, and find similarity solutions and small time asymptotic solutions for prescribed shear and speed cases respectively. These solutions are used as initial solutions for the two plate case. The motion of the yield surface and the time taken for the entire material to yield are investigated. The problems considered here are two dimensional representations of some control devices, for example the light duty clutch, which consists of two corotating, coaxial discs separated by a layer of electrorheological material. In this application it is useful to know the time taken for all the material to yield.     

The viscoelastic and viscoplastic behavior of low-density polyethylene

Three series of tensile tests with constant cross-head speeds (ranging from 5 to 200 mm/min), tensile relaxation tests (at strains from 0.03 to 0.09) and tensile creep tests (at stresses from 2.0 to 6.0 MPa) are performed on low-density polyethylene at room temperature. Constitutive equations are derived for the time-dependent response of semicrystalline polymers at isothermal deformation with small strains. A polymer is treated as an equivalent heterogeneous network of chains bridged by temporary junctions (entanglements, physical cross-links and lamellar blocks). The network is thought of as an ensemble of meso-regions linked with each other. The viscoelastic behavior of a polymer is modelled as thermally-induced rearrangement of strands (separation of active strands from temporary junctions and merging of dangling strands with the network). The viscoplastic response reflects mutual displacement of meso-domains driven by macro-strains. Stress–strain relations for uniaxial deformation are developed by using the laws of thermodynamics. The governing equations involve five material constants that are found by fitting the observations. Fair agreement is demonstrated between the experimental data and the results of numerical simulation. It is shown that observations in conventional creep tests reflect not only the viscoelastic, but also the viscoplastic behavior of an ensemble of meso-regions.     

A model for high temperature creep of single crystal superalloys based on nonlocal damage and viscoplastic material behavior

A model for high temperature creep of single crystal superalloys is developed, which includes constitutive laws for nonlocal damage and viscoplasticity. It is based on a variational formulation, employing potentials for free energy, and dissipation originating from plasticity and damage. Evolution equations for plastic strain and damage variables are derived from the well-established minimum principle for the dissipation potential. The model is capable of describing the different stages of creep in a unified way. Plastic deformation in superalloys incorporates the evolution of dislocation densities of the different phases present. It results in a time dependence of the creep rate in primary and secondary creep. Tertiary creep is taken into account by introducing local and nonlocal damage. Herein, the nonlocal one is included in order to model strain localization as well as to remove mesh dependence of finite element calculations. Numerical results and comparisons with experimental data of the single crystal superalloy LEK94 are shown.     

The viscoplastic behavior of hastelloy-X single crystal

A viscoplastic constitutive model for Hastelloy-X single crystal material is developed based on crystallographic slip theory. The constitutive model was constructed for use in a viscoplastic self-consistent model for isotropic Hastelloy-X polycrystalline material, which has been described in a recent publication. It is found that, by using the slip geometry known from the metallurgical literature, the anisotropic response can be accurately predicted. The model was verified by using tension and torsion data taken at 982°C (1800°F). The constitutive model used on each slip system is a simple unified visoplastic power law model in which weak latent interaction effects are taken into account. The drag stress evolution equations for the octahedral system are written in a hardening/recovery format in which both hardening and recovery depend on separate latent interaction effects between the octahedral crystallographic slip systems. The strain rate behavior of the single crystal material is well correlated by the constitutive model in uniaxial and torsion tests, but it is necessary to include latent information effects between the octahedral slip systems in order to obtain the best possible representation of biaxial cyclic strain rate behavior. Finally, it was observed that the single crystal exhibited dynamic strain aging at 871°C (1600°F). Similar dynamic strain aging occurs at 649°C (1200°F) in the polycrystalline version of the alloy.     

Hydraulic impact of “exponential” and nonlinearly viscoplastic media in pipes made of a viscoplastic material

Differential equations are derived and the hydraulic impact process for “exponential” and nonlinearly viscoplastic media in pipes made of a viscoelastic material is analyzed. Hydraulic impact problems for actual media in pipes has been repeatedly treated in the literature [1–4]. The hydraulic impact of a viscous and linearly viscoplastic media in pipes made of an elastic and viscoelastic material was studied in this work. It is well known [5] that many media in the region of low and moderate shear rates reveal a nonlinearity of the flow curve (oil, drilling fluids, polymer solutions and melts, loaded fuels, fuel mixtures, blood, etc.). It should be noted that flexible pipes made of natural materials (pipe boreholes made of polymer materials, membranes of blood vessels, etc.) are described by complicated rheological equations of state for viscoelastic media. Thus a calculation of the influence of nonlinearity of these media and of the viscoelastic properties of the pipe material on the hydraulic impact process is of theoretical and practical interest in many engineering problems.     

Observation and modeling of the anisotropic visco-hyperelastic behavior of a rubberlike material

The visco-hyperelastic behavior of a filled rubberlike material has been studied experimentally by large deformation cyclic uniaxial loadings, and an anisotropy induced by the Mullins effect has been demonstrated. By applying a generalized Maxwell model to a set of material directions, damage could be included in order to reproduce the stress softening due to the Mullins effect. This induces also an anisotropic mechanical response, and the model compares favorably with the experimental measures.     

Experimental investigation and modeling of non-monotonic creep behavior in polymers

The deformation behavior of two unfilled engineering thermoplastics, ultra high molecular weight polyethylene (UHMWPE) and polycarbonate (PC), has been investigated in creep test conditions. It has been found that a loading history (prior to the creep test) comprising of loading to a maximum stress or strain value followed by partial unloading to arrive at the target stress value can greatly modify the strain-time behavior. Under such a test protocol, while the expected increase in strain during creep (constant tensile load) is observed, at relatively low creep stresses specimens have also demonstrated a monotonic decrease in strain. In an intermediate stress range, specimens have demonstrated time dependent behavior comprising of a transition from decreasing to increasing strain during creep in tension. This paper presents experimental results to delineate these findings and explore the effect of prior strain rate on the qualitative and quantitative changes in the output (strain-time) behavior. Furthermore, modification of the viscoplasticity theory based on overstress (VBO) model into a double element configuration is introduced. These changes confer upon the model the ability to yield non-monotonic behavior in creep, and supporting simulation results have been included. These changes, therefore, allow the model to simulate strain rate sensitivity, creep, relaxation, and recovery behavior, but more importantly address the issue of non-monotonic changes in creep and relaxation when a loading history involves some degree of unloading.     

Development of viscoplastic constitutive equation through biaxial material testing

Internal-state variables have been used to represent the deformation historyin the recently proposed viscoplastic constitutive equations. In the current study, creep tests under nonproportional loadings were used to study the relative roles played by the internal-state varaibles in the constitutive equation by tracing the strain trajectory in strain space for a given stress trajectory in stress space. An experimental approach to studying the evolution rule for the tensorial state variable is also proposed. The experimental results on 2618-T61 aluminum alloy suggest that the scalarstate variable should play a much more dominant role than the tensorial state varaible in the constitutive modeling of 2618-T61 aluminum alloy.     

Characteristics of a viscoplastic material in the couette flow

A model governing a steady flow of a viscoplastic material between coaxial cylinders is proposed. Nonlinear velocity sensitivity typical of superplastic materials is taken into account. An algorithm of calculating the characteristics of the material is developed. The algorithm is based on the experimental data on moments and angular velocities of the rotating coaxial cylinders. The stability of the algorithm to errors in the initial data is estimated.     

The flow of viscoplastic material between two concentric spheres

The motion of a viscoplastic medium between two concentric spheres is considered upon rotation of one sphere with constant angular velocity. This problem is solved by an heuristic iterative method. The boundary of the stagnation zones is found and its specific shape is shown. The flow characteristics versus the parameter of the medium are obtained. Voronezh State Engineering Academy, Voronezh 394017. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 1, pp. 133–139, January–February, 1999.     

Mechanical characterisation of a viscoplastic material sensitive to hydrostatic pressure

The present paper deals with the characterisation of the static mechanical behaviour of an energetic material all along its lifespan. The material behaviour is viscoplastic, damageable and sensitive to hydrostatic pressure. For such materials, existing models have generally been developed in the framework of transient dynamic behaviour. These models are not suitable for a static study. Therefore a specific experimental protocol and an associated model are developed. Characterisation is derived from both uniaxial compressive, tensile tests and triaxial tests. Plastic behaviour is described by means of a parabolic yield criterion and a new hardening law. Non-associated plastic flow rule and isotropic damage complete the model. The performance of the model is illustrated through the simulation of various loading paths.     

Problems and trends in mathematical modeling

Main problems and trends in mathematical modeling — a new line of research into various processes and phenomena — are formulated. The status and future prospects are analyzed using as an example the mechanics of continuous media. Emphasis is on two stages of modeling — the selection of physicomathematical models of the mechanics of continuous media and numerical algorithms of solution.     

Time-dependent creep of a dual-phase viscoplastic material with lamellar structure

Exact solutions for the time-dependent creep behavior of a two-phase material with a lamellar microstructure are derived as a function of volume concentration and the properties of its constituents. Each phase is taken to be elastic–viscoplastic, exhibiting work-hardening characteristics. The derivation takes advantage of the condition of interfacial discontinuities over the interfaces, with a result given in a rate form for the general combined loading. Specific overall creep strains are presented along five distinctive loading directions for two kinds of viscoplastic composites: one involving an elastic and a viscoplastic phases and the other with dual viscoplastic phases. In addition to providing insightful information for the overall time-dependent creep, the exact nature of the results can also serve as a bench mark to test the accuracy of the approximate theories. In this light a secant-viscosity approach recently developed for a particle-reinforced solid (Li and Weng (1997). A secant-viscosity approach to the time-dependent creep of an elastic–viscoplastic composite. J. Mech. Phys. Solids, 45, 1069) is extended to a lamellar structure and the results are tested against these exact solutions. Comparison between the two indicates that the secant-viscosity concept is a sufficiently accurate one and it can be applied to composites with other types of microgeometries.     

Experimental analysis and modeling of biaxial mechanical behavior of woven composite reinforcements

This paper presents experimental studies on the mechanical behavior of fiber fabrics using a biaxial tensile device based on two deformable parallelograms. The cross-shaped specimens are well adapted to fabrics because of their lack of shear stiffness. Tension versus deformation curves, for different strain ratios, are obtained in the case of composite woven reinforcements used in aeronautic applications. It is shown that the tensile behavior of the fabric is strongly nonlinear due to the weaving undulations and the yarn contraction, and that the phenomenon is clearly biaxial. A constitutive model is described and identified from the experimental data. The essential role played by the yarn crushing will be pointed out.     

Parameter identification for improved viscoplastic model considering dynamic recrystallization

The successful application of viscoplastic model considering dynamic recrystallization depends on how well the parameters are identified. However, it is difficult to obtain satisfactory parameters using conventional parameter identification methods. The reasons are due to difficulties in obtaining homogeneous deformation, high complexity of physical process described by the model and large number of parameters. In this paper, the material parameters are identified by inverse analysis. Global information on objective function is firstly studied by an improved uniform random sampling method; secondly, a hybrid global optimization method, which combines the genetic algorithm, the Levenberg–Marquardt algorithm, the augmented Gauss–Newton algorithm and the flexible tolerance method, is constructed and an inverse analysis numerical procedure, which combines the proposed optimization method with the finite element analysis, is proposed; at last, a set of satisfactory material parameters for 26Cr2Ni4MoV is obtained by the proposed inverse analysis numerical procedure.