Constitutive Modelling of Superplastic AA-5083

Similar documents
EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

Module 3: Element Properties Lecture 1: Natural Coordinates

Finite Element Modelling of truss/cable structures

Lifetime prediction of EP and NBR rubber seal by thermos-viscoelastic model

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

DUE: WEDS FEB 21ST 2018

Inductance Calculation for Conductors of Arbitrary Shape

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST

CHAPTER 9 CONCLUSIONS

Numerical Transient Heat Conduction Experiment

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

NUMERICAL DIFFERENTIATION

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

Army Ants Tunneling for Classical Simulations

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Chapter - 2. Distribution System Power Flow Analysis

Kernel Methods and SVMs Extension

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

Week 9 Chapter 10 Section 1-5

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH

PLATE BENDING ELEMENTS

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Second Order Analysis

Experiment 1 Mass, volume and density

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

One-sided finite-difference approximations suitable for use with Richardson extrapolation

DESIGN OPTIMIZATION OF CFRP RECTANGULAR BOX SUBJECTED TO ARBITRARY LOADINGS

Modeling of Dynamic Systems

Numerical Nonlinear Analysis with the Boundary Element Method

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES

Linear Approximation with Regularization and Moving Least Squares

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Influence of longitudinal and transverse bulkheads on ship grounding resistance and damage size

In this section is given an overview of the common elasticity models.

Chapter 13: Multiple Regression

An Interactive Optimisation Tool for Allocation Problems

FUZZY FINITE ELEMENT METHOD

THE CURRENT BALANCE Physics 258/259

Torsion Stiffness of Thin-walled Steel Beams with Web Holes

Global Sensitivity. Tuesday 20 th February, 2018

Modal Strain Energy Decomposition Method for Damage Detection of an Offshore Structure Using Modal Testing Information

16 Circular Footing on a Semi-infinite Elastic Medium

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Three-dimensional eddy current analysis by the boundary element method using vector potential

The Finite Element Method

Implicit Integration Henyey Method

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE

x = , so that calculated

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE

Numerical Heat and Mass Transfer

Solution Thermodynamics

Supplementary Notes for Chapter 9 Mixture Thermodynamics

One Dimensional Axial Deformations

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

MODELLING OF ELASTO-STATICS OF POWER LINES BY NEW COMPOSITE BEAM FINITE ELEMENT Bratislava

RELIABILITY ASSESSMENT

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

The Minimum Universal Cost Flow in an Infeasible Flow Network

Chapter 11: Simple Linear Regression and Correlation

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

ˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m)

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Supplemental document

Canonical transformations

Note 10. Modeling and Simulation of Dynamic Systems

Linear Regression Analysis: Terminology and Notation

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

Inexact Newton Methods for Inverse Eigenvalue Problems

A Modified Neuber Method Avoiding Artefacts Under Random Loads

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Effect of loading frequency on the settlement of granular layer

9.2 Seismic Loads Using ASCE Standard 7-93

Pressure Measurements Laboratory

EEE 241: Linear Systems

CIVL 8/7117 Chapter 10 - Isoparametric Formulation 42/56

829. An adaptive method for inertia force identification in cantilever under moving mass

Identification of Instantaneous Modal Parameters of A Nonlinear Structure Via Amplitude-Dependent ARX Model

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Analytical and Numerical Analysis of Free Bulge Tube Hydroforming

Δ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions

THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES

2 Finite difference basics

Instituto Tecnológico de Aeronáutica FINITE ELEMENTS I. Class notes AE-245

Construction of Serendipity Shape Functions by Geometrical Probability

Appendix B. The Finite Difference Scheme

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

The Feynman path integral

Supplementary Materials for

Transcription:

TECHNISCHE MECHANIK, 3, -5, (01, 1-6 submtted: September 19, 011 Consttutve Modellng of Superplastc AA-5083 G. Gulano In ths study a fast procedure for determnng the constants of superplastc 5083 Al alloy at 73 K s analyzed. To evaluate the dsplacement and the thckness evolutons at the dome apex of the metal sheet bulge formng experments at constant pressure are performed. The dsplacement-tme curves at the constant pressures of 0.30, 0.35, 0.40, 0.45 and 0.50 MPa are obtaned by laser measurements whereas many bulgng tests at the same gas pressures stopped up at prefxed ntervals of tme were carred out to evaluate the thckness at the dome apex. Drect measurements of the thckness were obtaned usng a centesmal mcrometer. It s known that an equbaxal stress state s present at the dome apex of the metal sheet. In contrast, at other locatons, a nonequbaxal stress state prevals. Due to the nonequbaxal stress state, the thckest locaton occurs at the sheet perphery and the thnnest locaton at the dome apex. Therefore, by usng the fnte element method, t s possble to accurately analyse the superplastc formng process n the real stress state condtons. The materal behavour s modelled by the power law relatonshp between the effectve flow stress, the effectve stran and the effectve stran-rate. In the power law relatonshp the stran-rate senstvty ndex, m, s a constant, n and K are varables that are functons of the effectve stran-rate. Expermental results allowed to easly determne the value of m. Moreover, by means of some numercal smulatons t was possble to determne the values of n and K and assocate them wth specfc stran-rates. The results of comparson of the numercal smulatons of a bulge formng process wth the expermental tests have shown good agreement and they ndcate that the materal constants are relable. 1 Introducton Ths paper descrbes the characterzaton of AA5083 alumnum sheet by usng an nverse analyss technque n order to mnmze the dfference between numercal data and expermental measurements. Usng the commercal fnte element software MSC.Marc, the bulgng process at constant pressure s smulated. A rgd-plastc flow formulaton s appled to the superplastc formng analyss (Zenkewcz, 1984. The mesh s composed of fournode, soparametrc elements used for axsymmetrc applcatons. As ths element uses blnear nterpolaton functons, the strans tend to be constant throughout the element. The stffness of ths element s formed usng four-pont Gaussan ntegraton. The element has two coordnates n the global z- and r-drecton and two degrees of freedom for node (Msc.Marc, 005. An ntegraton scheme whch mposes a constant dlatatonal stran constrant on the element s used. The deformable body correspondng to the metal sheet s dscretzed wth two layers of 18 contnuum elements wth 4 nodes. The de s consdered as a perfectly rgd body. Constant pressure s appled as a dstrbuted load. Because of the symmetry of the geometry, the load and the constrant condtons, half of the cross-secton of the sheet metal s analysed. It s necessary to lock the movement of the nodes along the axs of symmetry n a drecton that s orthogonal of the axs of symmetry, to avod penetraton of the adjacent elements. Moreover, t s necessary to mpose constrant condtons on the perphery of the sheet n order to smulate the acton of a blank holder. Fgure 1a shows the shape of the de, the fnte element mesh of the sheet and the dstrbuted load over the edges of the sheet. The materal behavour s modelled by the power law relatonshp between the effectve flow stress, s, the effectve stran, e, and the effectve stran-rate, e &, by usng the followng equaton n m s = Ke e& (1 where the stran-rate senstvty ndex, m, s a constant, n and K are varables that are functons of the effectve stran-rate. The propertes of the materal and the pressure value, appled to the top part of the sheet metal, are ncluded usng a subroutne. 1

Materal and Formng Tests Alumnum alloy, AA5083 n the form of a 1mm thck sheet s characterzed by usng bulge formng tests at the constant pressures of 0.30, 0.35, 0.40, 0.45 and 0.50 MPa. The bulgng tests are carred out at temperature of 73 K (450 C. A crcular 79-mm-dameter specmen s nterposed between two heated steel des. The upper de has a crcular geometry wth an aperture radus of 30.0mm and a de entry radus of.0mm. The lower de works as a blankholder. From the lower de a pressure gas acts on one sde of the sheet forcng t to expand nto the upper de cavty. The ar pressure s obtaned by a compressor and t s regulated by a proportonal valve. A pressure transducer s nserted on the ar njecton lne to ensure that the pressure nsde the de s equal to the one set for the bulgng test. The specmen temperature s ndrectly controlled through the temperature of the heatng bands wrapped around the des and connected to an electrcal feeder. The temperature of the bands s measured and checked by a LabVIEW program. Fgure 1b shows the heatng bands wrapped around the des. Further detals about the formng system are reported n (Govnco, Gulano and Testa, 010. (a (b Fgure 1: a Schematc vew of the de and sheet mesh and b formng des for bulgng test The objects of the expermental actvty are to evaluate the dsplacement and the thckness evolutons at the dome apex of the metal sheet. The dsplacement-tme curves are obtaned by laser measurements whereas many bulgng tests at the same gas pressures stopped up at prefxed ntervals of tme were carred out to evaluate the thckness at the dome apex as showed n (Gulano, 010 for a magnesum-based alloy. Drect measurements of the thckness are obtaned usng a centesmal mcrometer. formng pressure (MPa formng tme to H=1 (s average stran-rate (s -1 0.30 545 4.1x10-4 0.35 173 6.13x10-4 0.40 16 8.83x10-4 0.45 896 1.3x10-3 0.50 671 1.66x10-3 Table 1: Expermental data Fgure shows the expermental trends of the dsplacement-tme curves for constant gas pressures of 0.30, 0.35, 0.40, 0.45 and 0.50 MPa, respectvely. Fgure 3 shows the dome apex thckness evoluton at the same pressures consdered n Fg.: t s possble to consder the relatonshp [1-(s/s 0 ] versus [H /(1+H ] as lnear. In ths fgure H s the normalzed polar dsplacement (defned as the rato of the dsplacement and the de radus, s 0 and s are the ntal and the current thckness at the dome apex. Table 1 summarzes data from the expermental tests. The average stran-rate s e t, where e s the effectve stran at the dome apex for H=1 and t s the formng tme.

Fgure : Dsplacement-tme curves n superplastc formng process at constant pressure 3 Characterzaton of the Superplastc 5083 Al Alloy The model proposed to descrbe the behavour of the 5083 Al alloy, relates the flow stress to the stran and the stran rate, accordng to the equaton (1. The materal constants were obtaned by usng the expermental bulgng tests at constant pressure and dfferent numercal smulatons. Under a gven constant pressure, the dsplacements at the dome apex are measured by a laser devce and recorded as a functon of tme. In order to obtan materal constants, the bulgng test data at fve dfferent pressure values were used: for each sngle pressure value, the tme step taken for the sheet to pass through the normalzed polar dsplacement H=0 to H=1 was measured. Fgure 3: Dome apex thckness evoluton n superplastc formng process at constant pressure 3

Fgure 4: Pressure-tme expermental data at H=1 Fgure 4 shows the pressure-tme expermental data at H=1. It s possble to note that -b t = ap ( where t s the formng tme at the polar dsplacement H=1, p s the pressure and a and b are constants. Therefore, the stran rate senstvty ndex, m, can be determned by adoptng the expresson presented n (Gulano, Franchtt, 007; Enkeev, Kruglov, 1995; Gulano, 009 ln( p1 p m = (3 ln t t ( 1 where t 1 and t are formng tmes at constant pressures p 1 and p, respectvely. In ths paper, from t-p expermental relatonshp eq. (, m may be defned from the followng expresson ( ln p m = - (4 ( lnt Then, from eq. (, one can obtan that m = 1 (5 b 0,00 0,010 n 0,000-8,0-7,6-7, -6,8-6,4-6,0-0,010 ln(stran-rate -0,00 Fgure 5: Relatonshp between n and the stran-rate 4

Moreover, by means of some numercal smulatons t was possble to determne the values of n and K and assocate them wth specfc stran-rates. The frst numercal phase requres the bulge formng smulatons at the same pressure consdered n the expermental actvty: m s well-known through eq. (5, to K s assgned an arbtrary value and n s made to vary n a sutable range. The numercal value of n s calculated mnmzng the functon Q by N ( ln ( lnt ( lnt N ( e ( e ( e Ê t - ˆ Ê - ˆ Q( n = Á + Á (6 Ë Ë where ( ln t N and ( e N represent the values of the normalzed tme and the equvalent stran obtaned from the numercal smulaton for a fxed value of n whereas ( ln t and ( e, nstead, are the normalzed tme and the equvalent stran obtaned expermentally. The parameter t or normalzed tme to be defned as follows t t H t (7 = =1 H =0.5 t H= 0. where t and H= 1 5 constant pressure-formng process. are, respectvely, the formng tmes taken to reach the confguratons H=1 and H=0.5 n a Fgure 5 shows that t s possble to consder the relatonshp = ln( e& n as lnear. The last phase of the procedure to determne the materal constants requres the value of the constant K to be determned: for each sngle pressure value, m s well-known through eq. (5, n s known by mnmzng the functon Q(n, and K s made to vary n a sutable range. The numercal value of K s calculated by mnmzng the functon F gven by where ( N ( th = 1 ( th = 1 ( t Ê ( Á - K = Á Ë H = 1 ˆ F (8 N H 1 t = and ( t H = 1 are the formng tmes necessary to obtan the same dome geometry at constant pressure (H=1 from the numercal smulaton and from the expermental tests, respectvely. Fgure 6 shows that a thrd-order polynomal functon can gve an adequate relatonshp between - ln( e& K. 19 188 K 184 180-8,0-7,6-7, -6,8-6,4-6,0 ln(stran-rate Fgure 6: Relatonshp between K and the stran-rate 4 FEM Results and Dscussons Usng the commercal fnte element software MSC.Marc, the bulgng process at the constant pressures of 0.30, 0.35, 0.40, 0.45 and 0.50 MPa has been smulated. The materal behavour s modelled by the power law relatonshp between the effectve flow stress, the effectve stran and the effectve stran-rate accordng to eq. (1. 5

In the stran rate nterval of 4.1x10-4 - 1.66x10-3 s -1 whch corresponds to the regon of the superplastc behavour for alumnum alloy AA5083, t s assumed n FEM calculatons that n = Aln ( e& + B (9 K = C ( [ ( ] [ ( ] 3 0 + C1 ln e & + C ln e& + C3 ln e& (10 where A, B, C are constants. The values for the constants are: A=-0.01, B=-0.148, C 0 =36.80, C 1 =5.816, C =11.317 and C 3 =0.684. For each sngle numercal smulaton, the dsplacement-tme curve and the [1-(s/s 0 ]-[H /(1+H ] curve have been obtaned. All the numercal analyses show that t s possble to consder the relatonshp [1-(s/s 0 ] versus [H /(1+H ] as lnear. In the range 0 H 1, the dstance between numercal calculatons and expermental measurements was evaluated. Ths dstance s expressed as t = tnum - t t (11 ( exp 100 ( s - s 100 % exp s num (1 % = exp sexp where t s the formng tme, s s the current thckness at the dome apex and the subscrpts num and exp ndcate the numercal and expermental values, respectvely. In all the bulgng tests, t% and s% were found to be less than 10%. 5 Conclusons Constant pressure superplastc bulgng tests of AA5083 alumnum sheet at 73 K are consdered n the present paper. The dsplacement-tme curve and the thckness evoluton at the dome apex of the metal sheet are evaluated. Commercal fnte element software s used to determne the constants of superplastc alloy. The results of comparson between the numercal smulatons of bulge formng process and the expermental tests show good agreement and they ndcate that the materal constants are relable. References Zenkewcz, O.C.: Numercal Analyss of Formng Processes. Wley, New York, (1984. MSC.Marc, User s Manual, Vol. B, (005. Govnco, G.; Gulano, G.; Testa, G.: Formng apparatus to nvestgate the effect of temperature on the superplastc behavor of alloys. NUMIFORM 010, Pohang, Republc of Korea, 13-17 June, (010. Gulano, G.: Modellng of superplastc formng of AZ31 magnesum alloy. AMPT 010, Pars, France, 4-7 October, (010. Gulano, G.; Franchtt, S.: On the evaluaton of superplastc characterstcs usng the fnte element method. Int. J. Mach. Tool Manufact., 47, (007, 471-476. Enkeev, F. U.; Kruglov, A.A.: An analyss of the superplastc formng of a thn crcular daphragm. Int. J. Mech. Sc., 37, (1995, 473-483. Gulano, G.: On the constants of the superplastc magnesum-based AZ31 alloy. ICMM, Dortmund, Germany, 15-17 September, (009. Address: Dr. Gllo Gulano, DMSAT, Unversty of Cassno, 03043 Cassno (FR ITALY. emal: gulano@uncas.t. 6