Linear constitutive model for mechano-sorptive creep in paper

The creep of paper is accelerated by moisture cycling. This effect is known as mechano-sorptive creep. It is assumed that this is an effect of transient stresses produced during moisture content changes in combination with non-linear creep behaviour of the fibres. The stresses produced by the moisture content changes are often much larger than the applied mechanical loads. If this is the case, the mechanical loads are only a perturbation to the internal stress state, and it will appear as if the mechano-sorptive creep is linear in stress. It is possible to take advantage of this feature. In the present report the pure moisture problem is first solved. The mechanical load is then treated as a perturbation of the solution to the moisture problem. Using this strategy, it is possible to linearize a non-linear network model for mechano-sorptive creep and to formulate a continuum model. As a result, the number of variables in the model is reduced. This is a significant improvement as it will be possible to use the linearized model to describe the material in a finite element program and solve problems with complicated geometries.     

A phenomenological model for anisotropic creep in a multipass weld metal

Creep strength of welded joints can be estimated by continuum damage mechanics. In this case constitutive equations are required for different constituents of the welded joint: the weld metal, the heat-affected zone, and the parent material. The objective of this paper is to model the anisotropic creep behavior in a weld metal produced by multipass welding. To explain the origins of anisotropic creep, a mechanical model for a binary structure composed of fine-grained and coarse-grained constituents with different creep properties is introduced. The results illustrate the basic features of the stress redistribution and damage growth in the constituents of the weld metal and agree qualitatively with experimental observations. The structural analysis of a welded joint requires a model of creep for the weld metal under multiaxial stress states. For this purpose the engineering creep theory based on the creep potential hypothesis, the flow rule, and assumption of transverse isotropy is applied. The outcome is a coordinate-free equation for secondary creep formulated in terms of the Norton–Bailey–Odqvist creep potential and three invariants of the stress tensor. The material constants are identified according to the experimental data presented in the literature.     

An anisotropic rotary diffusion model for fiber orientation in short- and long-fiber thermoplastics

The Folgar–Tucker model, which is widely-used to predict fiber orientation in injection-molded composites, accounts for fiber–fiber interactions using isotropic rotary diffusion. However, this model does not match all aspects of experimental fiber orientation data, especially for composites with long discontinuous fibers. This paper develops a fiber orientation model that incorporates anisotropic rotary diffusion. From kinetic theory we derive the evolution equation for the second-order orientation tensor, correcting some errors in earlier treatments. The diffusivity is assumed to depend on a second-order space tensor, which is taken to be a function of the orientation state and the rate of deformation. Model parameters are selected by matching the experimental steady-state orientation in simple shear flow, and by requiring stable steady states and physically realizable solutions. Also, concentrated fiber suspensions align more slowly with respect to strain than models based on Jeffery's equation, and we incorporate this behavior in an objective way. The final model is suitable for use in mold filling and other flow simulations, and it gives improved predictions of fiber orientation for injection molded long-fiber composites.     

An intermediate model method for obtaining a discrete relaxation spectrum from creep data

A convenient method is described for obtaining a discrete stress relaxation spectrum from linear viscoelastic creep data by means of a three-stage process. In stage one, a discrete retardation spectrum is fitted to the creep data using a least squares procedure, subject to the constraint that the discrete spectrum must be a specified order of polynomial function of the retardation time. In stage two, the resulting generalised Voigt model is solved numerically for an imposed step in strain, to determine the stress relaxation modulus function of time. In stage three, a discrete relaxation spectrum is fitted to the calculated stress relaxation modulus function. Although three stages are involved instead of the usual two, the procedure has been found to have certain practical advantages. These advantages make it suitable for the generation of relaxation spectra needed in viscoelastic stress analyses of solids, for example by the finite element method. In order to illustrate the proposed procedure it is applied to both artificial data and experimental creep data for poly(methyl methacrylate) at 70°C and at the glass transition.     

An anisotropic micromechanics model for predicting the rafting direction in Ni-based single crystal superalloys

An anisotropic micromechanics model based on the equivalent inclusion method is developed to investigate the rafting direction of Ni-based single crystal superalloys. The micromechanical model considers actual cubic structure and orthogonal anisotropy properties. The von Mises stress, elastic strain energy density, and hydrostatic pressure in different inclusions of micromechanical model are calculated when applying a tensile or compressive loading along the [001] direction. The calculated results can successfully predict the rafting direction for alloys exhibiting a positive or a negative mismatch, which are in agreement with pervious experimental and theoretical studies. Moreover, the elastic constant differences and mismatch degree of the matrix and precipitate phases and their influences on the rafting direction are carefully discussed.     

A hardening model for anisotropic materials

Results from proportionate-loading experiments under plane-stress states illustrate the existence of a field of uniform-hardening potentials for the yield and creep deformation behavior of isotropic, slightly anisotropic, aelotropic and orthotropic polycrystalline materials in the initially strain-free condition. For two different plane-stress states, it is shown that a linear functional relationship holds between the plastic-strain increment ratio and the stress ratio in these materials and that, consequently, the field is adequately modeled by a uniform-hardening anisotropic function that is quadratic in the components of deviatoric stress. Anisotropic plane-stress yield functions are formulated for any stage in the deformation process by combining the uniform-hardening function with the kinematic-hardening rule. The resulting functions, which correspond to rigid translations of initial yield loci according to Ziegler's rule, provide good agreement with experimental observations on a marked Bauschinger effect and an absence of cross hardening.     

Multi-mechanism anisotropic model for granular materials

By representing the assembly by a simplified column model, a constitutive theory, referred to as sliding–rolling theory, was recently developed for a two-dimensional assembly of rods subjected to biaxial loading, and then extended to a three-dimensional assembly of spheres subjected to triaxial (equibiaxial) loading. The sliding–rolling theory provides a framework for developing a phenomenological constitutive law for granular materials, which is the objective of the present work. The sliding–rolling theory provides information concerning yield and flow directions during radial and non-radial loading. In addition, the theory provides information on the role of fabric anisotropy on the stress–strain behavior and critical state shear strength. In the present paper, a multi-axial phenomenological model is developed within the sliding–rolling framework by utilizing the concepts of critical state, classical elasto-plasticity and bounding surface. The resulting theory involves two yield surfaces and falls within the definition of the multi-mechanism models. Computational issues concerning the solution uniqueness for stress states at the corner of yield surfaces are addressed. The effect of initial and induced fabric anisotropy on the constitutive behavior is incorporated. It is shown that the model is capable of simulating the effect of anisotropy, and the behavior of loose and dense sands under drained and undrained loading.     

An elasto-plastic constitutive model with plastic strain rate potentials for anisotropic cubic metals

An elasto-plastic constitutive model with the plastic strain rate potential was developed for finite element analysis. In the model, isotropic-kinematic hardening was incorporated under the plane stress condition for anisotropic sheet cubic metal forming analysis. The formulation is general enough for any homogeneous plastic strain rate potential (with the first-order homogeneous effective strain rate) but the plastic strain rate potential Srp2004-18p was considered here. Attention was focused on the development of the elasto-plastic transition criterion and the effective stress update algorithm. Also, to assure the quadratic convergence rate in Newton’s method, the elasto-plastic tangent modulus was analytically derived. Accuracy and convergence of the stress update algorithm were assessed by the iso-error maps, whereas stability of the algorithm was confirmed by analytical procedure. Validations were performed for the examples of the circular cup drawing, 2D draw-bending and unconstrained cylindrical bending tests, utilizing aluminum sheet alloys.     

An anisotropic elastic-plastic constitutive model for single and polycrystalline metals. I-theoretical developments

An anisotropic elastic-plastic constitutive model for single and polycrystalline metals is proposed. The anisotropic hardening of single crystals, at first, is discussed with the viewpoint of yield surface and a new formulation of it is proposed. Then, a model for the anisotropic hardening of polycrystals is suggested by increasing the number of slip systems and incorporating the interaction of all slip systems. The interaction of grains through grain boundaries is shown to be similar to, and incorporated into, the interaction of slip systems in grains. The numerical predictions and their comparisons with experiments will follow in Part II of this paper.     

The spherical vessel with anisotropic creep properties considering large strains

The effect of material anisotropy on creep of pressurized thick-walled spherical vessel has been discussed considering the large strain theory. It is found that the creep strain varies with varying anisotropy of the material. The results obtained for the anisotropic case have been compared with those obtained for the isotropic case and it is observed that the stress and strain distribution in the wall of the vessel is strikingly different for the two cases.     

A dislocation density based material model to simulate the anisotropic creep behavior of single-phase and two-phase single crystals

The primary and secondary creep behavior of single crystals is observed by a material model using evolution equations for dislocation densities on individual slip systems. An interaction matrix defines the mutual influence of dislocation densities on different glide systems. Face-centered cubic (fcc), body-centered cubic (bcc) and hexagonal closed packed (hcp) lattice structures have been investigated. The material model is implemented in a finite element method to analyze the orientation dependent creep behavior of two-phase single crystals. Three finite element models are introduced to simulate creep of a γ′ strengthened nickel base superalloy in 〈1 0 0〉, 〈1 1 0〉 and 〈1 1 1〉 directions. This approach allows to examine the influence of crystal slip and cuboidal microstructure on the deformation process.     

A boundary model for the thermal creep problem

A new integral boundary condition model for determining the particle velocity distribution function in a Knudsen layer, which makes it possible to improve the accuracy of the moment approximations, is proposed. The results for various models are compared with the exact numerical solution of the problem.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.6, pp. 160–163, November–December, 1993.This research was made possible by financial support from the Exploratory Research Program, US EPA and by a Fulbright grant for Professor I. N. Ivchenko.