Anisotropic hardening model based on non-associated flow rule and combined nonlinear kinematic hardening for sheet materials

 A material model for more effective analysis of plastic deformation of sheet materials is presented in this paper. The model is capable of considering the following aspects of plastic deformation behavior of sheet materials: the anisotropy in yielding stresses in different directions by using a quadratic yield function (based on Hill’s 1948 model and stress ratios), the anisotropy in work hardening by introducing non-constant flow stress hardening in different directions, the anisotropy in plastic strains in different directions by using a quadratic plastic potential function and non-associated flow rule (based on Hill’s 1948 model and plastic strain ratios, r-values), and finally some of the cyclic hardening phenomena such as Bauschinger’s effect and transient behavior for reverse loading by using a coupled nonlinear kinematic hardening (so-called Armstrong-Frederick-Chaboche model). Basic fundamentals of the plasticity of the model are presented in a general framework. Then, the model adjustment procedure is derived for the plasticity formulations. Also, a generic numerical stress integration procedure is developed based on backward-Euler method (so-called multistage return mapping algorithm). Different aspects of the model are verified for DP600 steel sheet. Results show that the new model is able to predict the sheet material behavior in both anisotropic hardening and cyclic hardening regimes more accurately. By featuring the above-mentioned facts in the presented constitutive model, it is expected that more accurate results can be obtained by implementing this model in computational simulations of sheet material forming processes. For instance, more precise results of springback prediction of the parts formed from highly anisotropic hardened materials or that of determining the forming limit diagrams is highly expected by using the developed material model.

Anisotropic characterization of crack growth in the tertiary flow of asphalt mixtures in compression

Asphalt mixtures exhibit primary, secondary, and tertiary stages in sequence during a rutting deterioration. Many field asphalt pavements are still in service even when the asphalt layer is in the tertiary stage, and rehabilitation is not performed until a significant amount of rutting accompanied by numerous macrocracks is observed. The objective of this study was to provide a mechanistic method to model the anisotropic cracking of the asphalt mixtures in compression during the tertiary stage of rutting. Laboratory tests including nondestructive and destructive tests were performed to obtain the viscoelastic and viscofracture properties of the asphalt mixtures. Each of the measured axial and radial total strains in the destructive tests were decomposed into elastic, plastic, viscoelastic, viscoplastic, and viscofracture strains using the pseudostrain method in an extended elastic-viscoelastic correspondence principle. The viscofracture strains are caused by the crack growth, which is primarily signaled by the increase of phase angle in the tertiary flow. The viscofracture properties are characterized using the anisotropic damage densities (i.e., the ratio of the lost area caused by cracks to the original total area in orthogonal directions). Using the decomposed axial and radial viscofracture strains, the axial and radial damage densities were determined by using a dissipated pseudostrain energy balance principle and a geometric analysis of the cracks, respectively. Anisotropic pseudo J-integral Paris' laws in terms of damage densities were used to characterize the evolution of the cracks in compression. The material constants in the Paris' law are determined and found to be highly correlated. These tests, analysis, and modeling were performed on different asphalt mixtures with two binders, two air void contents, and three aging periods. Consistent results were obtained; for instance, a stiffer asphalt mixture is demonstrated to have a higher modulus, a lower phase angle, a greater flow number, and a larger n1 value (exponent of Paris' law). The calculation of the orientation of cracks demonstrates that the asphalt mixture with 4% air voids has a brittle fracture and a splitting crack mode, whereas the asphalt mixture with 7% air voids tends to have a ductile fracture and a diagonal sliding crack mode. Cracks of the asphalt mixtures in compression are inclined to propagate along the direction of the external compressive load. © 2014 American Society of Civil Engineers.

Anisotropic damage mechanics as a novel approach to improve pre-and post-failure borehole stability analysis

Anisotropic damage distribution and evolution have a profound effect on borehole stress concentrations. Damage evolution is an irreversible process that is not adequately described within classical equilibrium thermodynamics. Therefore, we propose a constitutive model, based on non-equilibrium thermodynamics, that accounts for anisotropic damage distribution, anisotropic damage threshold and anisotropic damage evolution. We implemented this constitutive model numerically, using the finite element method, to calculate stress–strain curves and borehole stresses. The resulting stress–strain curves are distinctively different from linear elastic-brittle and linear elastic-ideal plastic constitutive models and realistically model experimental responses of brittle rocks. We show that the onset of damage evolution leads to an inhomogeneous redistribution of material properties and stresses along the borehole wall. The classical linear elastic-brittle approach to borehole stability analysis systematically overestimates the stress concentrations on the borehole wall, because dissipative strain-softening is underestimated. The proposed damage mechanics approach explicitly models dissipative behaviour and leads to non-conservative mud window estimations. Furthermore, anisotropic rocks with preferential planes of failure, like shales, can be addressed with our model.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter = 9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 10(8). Under these conditions the convection is turbulent, and the time-averaged velocity at any point is `zero'. The Reynolds number based on the Taylor microscale, Re-lambda, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6-7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as (RaSc1/2)-Sc-1/2, and the Reynolds number would scale as (RaSc-1/2)-Sc-1/2. The velocity and the flux measurements appear to be consistent with the Ra-1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are Much higher compared to what would be obtained in Rayleigh-Benard (R-B) convection for similar density differences.

Depressed weir with two unequal sheetpiles in anisotropic soil

An analytical solution is presented, making use of the Schwartz-Christoffel transformation, for determining the seepage characteristics for the problem of flow under a weir having two unequal sheetpiles at the ends and embedded in an anisotropic porous medium of finite thickness. Results for several particular cases of simple hydraulic structures can be obtained from the general solution presented. Numerical results in nondimensional form have been given for quantity of seepage and exit gradient distribution for various conditions in the equivalent transformed isotropic section and, by making use of the physical parameters in the actual anisotropic plane and the set of transformation relations given, these quantities (seepage loss, exit gradient) can be interpreted in the actual anisotropic physical plane.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter=9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 108. Under these conditions the convection is turbulent, and the time-averaged velocity at any point is ‘zero’. The Reynolds number based on the Taylor microscale, Reλ, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6–7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as Ra1/2Sc1/2, and the Reynolds number would scale as Ra1/2Sc−1/2. The velocity and the flux measurements appear to be consistent with the Ra1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are much higher compared to what would be obtained in Rayleigh–Bénard (R–B) convection for similar density differences.

Dynamic contact angle beating from drops impacting onto solid surfaces exhibiting anisotropic wetting

This paper reports an experimental investigation of low Weber number water drops impacting onto solid surfaces exhibiting anisotropic wetting. The wetting anisotropy is created by patterning the solid surfaces with unidirectional parallel grooves. Temporal measurements of impacting drop parameters such as drop base contact diameter, apparent contact angle of drop, and drop height at the center are obtained from high-speed video recordings of drop impacts. The study shows that the impact of low Weber number water drops on the grooved surface exhibits beating phenomenon in the temporal variations of the dynamic contact angle anisotropy and drop height at the center of the impacting drop. It is observed that the beating phenomenon of impacting drop parameters is caused by the frequency difference between the dynamic contact angle oscillations of impacting drop liquid oriented perpendicular and parallel to the direction of grooves on the grooved surface. The primary trigger for the phenomenon is the existence of non-axisymmetric drop flow on the grooved surface featuring pinned and free motions of drop liquid in the directions perpendicular and parallel to the grooves, respectively. The beat frequency is almost independent of the impact drop Weber number. Further experimental measurements with solid surfaces of different groove textures show that the grooved surface with larger wetting anisotropy may be expected to show a dominant beating phenomenon. The phenomenon is gradually damped out with time and is fully unrecognizable at higher drop impact Weber numbers. (C) 2011 Elsevier B.V. All rights reserved.

Tensile flow and work hardening behavior of hot cross-rolled AA7010 aluminum alloy sheets

The present work describes the tensile flow and work hardening behavior of a high strength 7010 aluminum alloy by constitutive relations. The alloy has been hot rolled by three different cross-rolling schedules. Room temperature tensile properties have been evaluated as a function of tensile axis orientation in the as-hot rolled as well as peak aged conditions. It is found that both the Ludwigson and a generalized Voce-Bergstrom relation adequately describe the tensile flow behavior of the present alloy in all conditions compared to the Hollomon relation. The variation in the Ludwigson fitting parameter could be correlated well with the microstructural features and anisotropic contribution of strengthening precipitates in the as-rolled and peak aged conditions, respectively. The hardening rate and the saturation stress of the first Voce-Bergstrom parameter, on the other hand, depend mainly on the crystallographic texture of the specimens. It is further shown that for the peak aged specimens the uniform elongation (epsilon(u)) derived from the Ludwigson relation matches well with the measured epsilon(u) irrespective of processing and loading directions. However, the Ludwigson fit overestimates the epsilon(u) in case of the as-rolled specimens. The Hollomon fit, on the other hand, predicts well the measured epsilon(u), of the as-rolled specimens but severely underestimates the epsilon(u), for the peak aged specimens. Contrarily, both the relations significantly overestimate the UTS of the as-rolled and the peak aged specimens. The Voce-Bergstrom parameters define the slope of e Theta-sigma plots in the stage-III regime when the specimens show a classical linear decrease in hardening rate in stage-III. Further analysis of work hardening behavior throws some light on the effect of texture on the dislocation storage and dynamic recovery.

Bioconvection in an anisotropic scattering suspension of phototactic algae

Phototaxis is a directed swimming response dependent upon the light intensity sensed by microorganisms. Positive phototaxis denotes motion directed towards the source of light and negative phototaxis is motion directed away from it. In this paper, we investigate the onset of bioconvection in a suspension of anisotropic scattering phototactic algae illuminated by collimated radiation at the top. The basic state of the system is defined by the zero fluid flow and the up and down swimming, caused by the positive and negative phototaxis, is balanced by the diffusion. A comprehensive numerical study of the linear stability is presented with particular emphasis on the forward scattering effect. The onset of bioconvection occurs either via a stationary mode or an oscillatory mode. The transition from a stationary mode to an oscillatory mode or vice versa has been observed as the anisotropic coefficient is varied for certain parameter values. (C) 2012 Elsevier Masson SAS. All rights reserved.

Instabilities in a fluid overlying an inclined anisotropic and inhomogeneous porous layer

In this paper, linear stability analysis on a Newtonian fluid film flowing under the effect of gravity over an inclined porous medium saturated with the same fluid in isothermal condition is carried out. The focus is placed on the effect of the anisotropic and inhomogeneous variations in the permeability of the porous medium on the shear mode and surface mode instabilities. The fluid-porous system is modelled by a coupled two-dimensional Navier-Stokes/Darcy problem. The perturbation equations are solved numerically using the Chebyshev collocation method. Detailed stability characteristics as a function of the depth ratio (the ratio of the depth of the fluid layer to that of the porous layer), the anisotropic parameter (the ratio of the permeability in the direction of the basic flow to that in the direction transverse to the basic flow) and the inhomogeneity functions are presented.

Anisotropic Deformation Of Polystyrene Particles By Mev Au Ion Irradiation

MeV An irradiation leads to a shape change of polystyrene (PS) and SiO2 particles from spherical to ellipsoidal, with an aspect ratio that can be precisely controlled by the ion fluence. Sub-micrometer PS and SiO2 particles were deposited on copper substrates and irradiated with Au ions at 230 K, using an ion energy and fluence ranging from 2 to 10 MeV and 1 x 10(14) ions/cm(2) to 1 x 10(15) ions/cm(2). The mechanisms of anisotropic deformation of PS and SiO2 particles are different because of their distinct physical and chemical properties. At the start of irradiation, the volume of PS particles decrease, then the aspect ratio increases with fluence, whereas for SiO2 particles the volume remains constant. (C) 2008 Elsevier B.V. All rights reserved.

Triple-region structure for turbulent-flow in a square duct - a finite-element approach

The Reynolds-averaged Navier-Stokes equations for describing the turbulent flow in a straight square duct are formulated with two different turbulence models. The governing equations are then expanded as a multi-deck structure in a plane perpendicular to the streamwise direction, with each deck characterized by its dominant physical forces as commonly carried out in analytical work using triple-deck expansion. The resulting equations are numerically integrated using higher polynomial (H-P) finite element technique for each cross-sectional plane to be followed by finite difference representation in the streamwise direction until a fully developed state is reached. The computed results using the two different turbulence models show fair agreement with each other, and concur with the vast body of available experimental data. There is also general agreement between our results and the recent numerical works anisotropic k-epsilon turbulence model.

Secondary flow near an undulatory surface induced by wall blocking effect

Prandtl's secondary mean motions of the second kind near an undulating surface were explained in terms of turbulent blocking effect and kinematic boundary conditions at the surface, and its order of magnitude was estimated. Isotropic turbulence is distorted by the undulating surface of wavelength λ and amplitude h with a low slope, so that h « λ. The prime mechanism for generating the mean flow is that the far-field Isotropic turbulence is distorted by the non-local blocking effect of the surface to become anisotropic axisymmetric turbulence near the surface with principal axis that is not aligned with the local curvature of the undulation. Then the local analysis can be applied and the mechanism is similar to the mean flow generation mechanism for homogeneous axisymmetric turbulence over a planer surface, i.e. gradients of the Reynolds stress caused by the turbulent blocking effect generate the mean motions. The results from this simple analysis are consistent with previous exact analysis in which the effects of curvature are strictly taken into account. The results also qualitatively agree with flow visualization over an undulating surface in a mixing-box.

Reynolds number influence on turbulent blocking effects of a rigid plane boundary (3rd report, direct numerical simulation on wall blocking effects for anisotropic turbulence)

The Reynolds number influence on turbulent blocking effects by a rigid plane boundary is studied using direct numerical simulation (DNS). A new forcing method proposed in the second report using Townsend's "simple model eddies" for DNS was extended to generate axisymmetric anisotropic turbulence. A force field is obtained in real space by sprinkling many space-filling "simple model eddies" whose centers are randomly but uniformly distributed in space. The axes of rotation are controlled in this study to generate axisymmetric anisotropic turbulence. The method is applied to a shear-free turbulent boundary layer over a rigid plane boundary and the blocking effects for anisotropic turbulence are investigated. The results show that stationary axisymmetric anisotropic turbulence is generated using the present method. Turbulence intensities near the wall showed good agreements with the rapid distortion theory (RDT) for small t (t ≪ TL), where TL. is the eddy turnover time. The splat effect (i. e. turbulence intensities of the components parallel to the surface are amplified) occurs near the boundary and the viscous effect attenuates the splat effect at the quasi steady state at low Reynolds number as for Isotropic turbulence. Prandtl's secondary flow of the second kind does not occur for low Reynolds number flows, which qualitatively agrees with previous observetion in a mixing-box.