Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities

This study aims to investigate the flow and fracture behavior of aluminum alloy 6082-T6 (AA6082-T6) at different strain rates and triaxialities. Two groups of Charpy impact tests were carried out to further investigate its dynamic impact fracture property. A series of tensile tests and numerical simulations based on finite element analysis (FEA) were performed. Experimental data on smooth specimens under various strain rates ranging from 0.0001~3400 s-1 shows that AA6082-T6 is rather insensitive to strain rates in general. However, clear rate sensitivity was observed in the range of 0.001~1 s-1 while such a characteristic is counteracted by the adiabatic heating of specimens under high strain rates. A Johnson-Cook constitutive model was proposed based on tensile tests at different strain rates. In this study, the average stress triaxiality and equivalent plastic strain at facture obtained from numerical simulations were used for the calibration of J-C fracture model. Both of the J-C constitutive model and fracture model were employed in numerical simulations and the results was compared with experimental results. The calibrated J-C fracture model exhibits higher accuracy than the J-C fracture model obtained by the common method in predicting the fracture behavior of AA6082-T6. Finally, the Scanning Electron Microscope (SEM) of fractured specimens with different initial stress triaxialities were analyzed. The magnified fractographs indicate that high initial stress triaxiality likely results in dimple fracture.


Introduction
Aluminum alloys, with low density, high strength to weight ratio, good ductility and excellent corrosion resistance, are widely used in aviation, aerospace, automotive, machinery manufacturing etc. 6082-T6 aluminum alloy, as an Al-Mg-Si alloy, is increasingly applied to manufacture high-speed trains due to its sufficient plasticity for extrusion, relatively high strength, excellent weldability, formability and machinability. It is well known that aluminum alloys have some mechanical properties like the dependence of mechanical properties on the a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 strain rate, stress state and temperature. Spigarelli et al. [1] analyzed the peak stress dependence on strain rate and temperature based on a series of torsion tests in a wide range of temperatures and strain rates. With regard to the 6082 aluminum mechanical properties, several fruitful studies were conducted, such as Tranə et al. [2], who have evaluated the drawing efficiency of 6082 0 temper aluminum alloy for cartridge tubes. They applied the Piecewise-Linear-Plasticity constitutive model and the Cockcroft-Latham fracture criterion to the simulations. Two important conclusions were drawn in their work that the stress triaxiality ratio dominated the material failure and the constitutive model for 6082 0 temper aluminum alloy was an exponential law. For 6082 aluminum alloy in T6 temper (solution heat treatment and artificial aging), little rate sensitivity and moderated anisotropy was found by Chen et al. [3].
Considerable quantity of classic researches [4][5][6] have been done to study the relation between mechanical properties and the strain rate. Comparing with other constitutive relations, the model proposed by Johnson and Cook [6], which describes the material stress response at large strain, various strain rates and elevated temperatures, is widely used in impact problems. Also, they provided material parameters of 12 different materials, Hu et al. [7] determined those for AerMet 100 steel, Feng et al. [8] for AZ31B magnesium alloy, Trajkovski et al. [9] for armor steel PROTAC 500. Singh et al. [10] derived the J-C parameters based on the dynamic increase factors of strengths for CP800 steel. Some authors, including Børvik et al. [11], Hou and Wang [12], and Li et al. [13] proposed modified J-C constitutive models to predict material mechanical response.
Fracture ductility, as the ability of a material to bear plastic deformation without failure, is of great significance for the application of a material. Studies of effect of stress state on fracture ductility may be dated back to the research of Ludwik and Scheu [14]. They hold the view that the strength-strain curve could be given by testing tensile specimens with circumferential notches of different depths. The stress triaxiality is defined as the ratio of the hydrostatic stress and the equivalent stress to depict the stress state of materials. Tests carried out by Bridgman [15] showed that a reduced stress triaxiality could increase the fracture strain greatly. Similar conclusion was drawn by Hancock and Mackenzie [16] through a series of tensile tests on prenotched round steel specimens.
Several failure criterions have been used in commercial FEA codes by now such as LS-Dyna, Abaqus, and Ansys. The Johnson-Cook failure criterion [17] has been widely used to study the dependence of fracture strain on the strain rate, stress triaxiality, and temperature through both experiments and numerical simulations. Considering the anisotropy of the material, Johnson-Cook failure criterions of AA5083-H116 were proposed by Clausen et al. [18] along different directions of the rolled material, and the 0˚direction specimens exhibit an opposite fracture strain tendency for various stress triaxialities with others. Erice et al. [19] calibrated the J-C fracture criterion of FV535 steel by averaging Bridgman's analysis and the numerical simulation result. Knowing that only high stress triaxiality is considered by the J-C failure criterion, Bao and Wierzbicki [20] studied the fracture ductility in the entire range of stress triaxialities by carrying out a series of upsetting, shear, and tension on 2024-T351 aluminum alloy, and suggested a segmented failure criterion. However, those few published studies of the fracture criterion coupled with dynamic constitutive mechanical behavior cannot match the increasingly engineering application of AA6082-T6. The objective of this paper is to bridge the gap.
In this paper, quasi-static and dynamic tensile tests for AA6082-T6 were carried out. A constitutive model was suggested to describe the dynamic mechanical behavior. Coupled with the constitutive model, a combined experimental and simulation method was used to calibrate the J-C fracture model. Further, the tensile tests were simulated to validate the constitutive model and fracture model. Magnified fractographs of smooth and notched specimens were recorded by SEM. In the end, this paper discusses the mechanical properties of the investigated material under various strain rates and stress triaxialities.

Experiment and method
Experimental procedure and materials Quasi-static and low strain rate tests. Experiments conducted in this study consist of axial tensile tests and notched specimen tensile tests. The former groups were set to study the influence of strain rate on mechanical behavior of aluminum alloy AA6082-T6, and the last one is for the influence of stress triaxiality on ductility. Three repeated experiments in each case were conducted to ensure the test reliability.
The chemical compositions (in wt. %) of aluminum alloy AA6082-T6 used in this study are listed in Table 1. The original material was rolled sheet, and all the specimens for the tensile experiment were cut with the length direction of the specimens coincident with the rolling direction of the sheet. For quasi-static (0.0001 s -1 ) and low strain rate (0.001~1 s -1 ) tension tests, the rectangular cross section flat specimens were used with total length of 259 mm, minimum width of 20 mm and thickness of 20 mm, as shown in Fig 1. The tensile tests were performed by MTS electronic universal testing machine (shown in Fig 2.) at room temperature 20˚C except high rate (100~10 4 s -1 ) tensile tests. The specimens were clamped with a wedge-shaped clamp in tension experiments, as shown in Fig 2. The displacement was measured by an extensometer with the resolution of 0.5 μm.
High strain rate tensile tests. For high rate tensile tests, a split Hopkinson tension bar (SHTB) testing apparatus was used to characterize the dynamic properties of the investigated material. The SHTB (shown in Fig 3(a)) device mainly consists of a striker, an incident bar, a transmitted bar, strain gauge, digital storage oscilloscope and an ultrahigh dynamic strainometer. Fig 3(b) is the geometry of the flat SHTB-specimen. Fig 3(c) shows a specimen sandwiched between the incident bar and transmission bar. The pressure bars are made of LY12 cryogenic aluminum alloy. The incident bar is 3300 mm in length and the transmitted bar is 1500 mm. The diameter is 20 mm for both bars. The tubular striker is made of the same material as the pressure bars. When the striker is propelled by the gas gun impacting the incident bar, an elastic wave is formed and travels through the incident bar. Once the wave reaches the specimen, due to the interaction with the specimen, part of the wave is reflected back into the incident bar and the remaining part, the transmitted wave, passes through the transmission bar. In the process of wave propagation, the specimen is under high rate tension, and the strain histories of those waves (the incident, reflected and transmitted wave) are recorded by the strain gauges mounted on both bars and collected by Tektronix 4104 digital storage oscilloscope. Fig 4 shows the strain waves at a strain rate of 800 s -1 .
In the high strain rate tensile tests, the strain rate was controlled by changing the length of the striker and the air pressure of the gas gun. The striker length varies from 300 mm to 1050 mm. The tests were conducted at three different strain rates, i.e., 800 s -1 , 2700 s -1 and 3400 s -1 . Three repeated tests under the same loading condition were carried out to ensure the reliability of the results. According to the principles of one-dimensional elastic-wave propagation in slender bars, the engineering stress, the engineering strain and the strain rate can be calculated from the reflected and transmitted wave amplitude via the following equations [10] where σ e (t), ε e (t) and _ εðtÞ are the engineering stress, the engineering strain and the strain rate, respectively. ε tw (t) and ε rw (t) refer to the amplitude of the transmitted wave and the reflected wave. E, C 0 and A 0 are Young's modulus, the stress wave speed and cross section area of the incident bar, respectively. A S and L represent cross section area and length of the specimen, respectively. The true stress and true strain (respectively denoted as σ t and ε t ) are given by the following equations Notched specimen tensile tests Studies [14][15][16][17][18][19][20] have been performed in the past that proved fracture ductility depends markedly on the triaxiality of the stress state. The triaxiality is usually represented by the dimensionless stress triaxiality ratio η, which is defined as Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities where σ m is the hydrostatic stress, σ eq is the von Mises equivalent stress, and σ 1 , σ 2 , σ 3 represent three principal stresses. For an uniaxial stress state where σ 1 is the only non-zero component, stress triaxiality η = 1/3. Hancock and Mackenzie [16] carried out a series of tensile tests on round pre-notched steel specimens. The initial stress triaxiality was calculated according to Bridgman's [15] analysis in their study.
where, r is the radius of the minimum cross-section and R is the radius of the circumferential notch. As (Eq 4) shows, for round pre-notched specimens the value of stress triaxiality increases when the radius of the notch decreases. To investigate the relation between fracture strain and stress triaxiality of AA6082-T6, three types of flat specimens with different notch size (notch radius 10mm, 40mm, 90mm respectively) were machined ( Fig 5) for tensile tests. Because of the different cross-section geometry of round and flat specimen, Bridgman's analysis is not suitable here. But the magnitude of the notch radius has the same influence on the stress triaxiality. That is, the specimen with the smaller notch radius exhibits the higher value of stress triaxiality. For each type of pre-notched specimens, tensile tests under quasi-static tensile loading were performed. The experimental conditions were exactly the same as those for the previously stated tensile tests on smooth specimens.

Material model
Johnson-Cook constitutive relation. The Johnson-Cook constitutive model, which is an empirical viscoplastic constitutive model, can well depict the work hardening, strain rate hardening and the thermal softening effect of metal materials. The J-C constitutive relation is shown as follows, where σ is the stress; ε is the equivalent plastic strain; _ ε Ã is the dimensionless plastic strain rate given by the equation of _ ε Ã ¼ _ ε= _ ε 0 in which _ ε 0 is the reference strain rate (0.0001 s -1 ); T Ã is the dimensionless temperature expressed as Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities PLOS ONE | https://doi.org/10.1371/journal.pone.0181983 July 31, 2017 environment, reference and material melt temperature respectively; A is the yield stress; B is the strain hardening coefficient; n is the strain hardening exponent; C is the rate sensitivity coefficient. During high strain rate loading specimens are reasonable to be assumed in adiabatic conditions. In this case the rate of temperature change is calculated through where ρ and C p and are the density and heat capacity of material. χ is the conversion factor of work into heat. The temperature rise can be estimated using the following equation where ε e is the strain corresponding to the maximum stress. ρ and C p are assumed as constants, while χ is taken as 0.9 for metals [11]. The equation can be rewritten as In this study all experiments were conducted at room temperature. For tensile tests at quasistatic and low strain rate, the thermal softening effect is neglected for simplicity. Therefore, the constitutive relation can be simplified as Here, the quasi-static strain rate is taken as 0.0001 s -1 . For quasi-static experiments, the value of _ ε Ã can naturally be taken as 1. Therefore, the constitutive model can be further simplified as where A is the true yield stress of the material at a strain rate of 0.0001 s -1 which is equal to the stress corresponding to 0.2% offset strain. By placing constant A on the left side of (Eq 10), the following equation can be obtained (Eq 11) can be used to fit the quasi-static experimental data in plastic deformation stage using the least square method. The values for each B and n can be directly obtained as parameters a and b respectively. At room temperature, the value of C can be obtained by fitting the strainstress data of different strain rate experiments. The constitutive model can be rewritten as where σ 0 is the yield stress under strain rate 0.0001 s -1 ; σ i is the yield stress under different strain rates; _ ε 0 is the reference strain rate (taken 0.0001 s -1 as mentioned earlier); _ ε represents strain rate. (Eq 12) can be converted to the following equation With this method, the constitutive relation of aluminum alloy 6082-T6 without thermal softening effect can be established effectively.
Johnson-Cook fracture model. Johnson and Cook proposed a fracture model in which the equivalent plastic fracture strain ε f is dependent on strain rate and temperature in addition to stress triaxiality. The general expression for the equivalent plastic strain at fracture given by J-C failure model is The constants D 1 , D 2 , D 3 are determined from the quasi-static tests on smooth and notched axis symmetric specimens and D 4 from various strain rate tests of smooth specimens. In the present study, all tests were conducted at room temperature, temperature influence on fracture strain is thus neglected. The J-C fracture model expression without temperature effect then is simplified as (Eq 15) where only four parameters need to be determined.
Calibration of the Johnson-Cook fracture model. During tensile tests, the strain field is not homogeneously distributed throughout the whole specimen and the maximum strain occurs on the fracture surface. Constant volume of the specimen and homogenous strain over the fracture surface are assumed and then the fracture strain [19] could be obtained through where A represent minimum sectional area on the specimen. a and b are the width and the height of specimens at minimum sectional area, respectively. Both of them are obtained from measurements. The subscript "0" and "f" represent before and after the test. Previous works [7][8][9][10][11][12][13] have calculated the stress triaxiality of round specimens with Bridgman's formulation ((Eq 4)) which is a valid method to obtain the stress triaxiality over the minimum cross section accurately. Actually, the stress triaxiality of the specimen is not a constant value in the plastic strain process [18,19]. The stress triaxiality changes with the geometry evolution of a specimen. What's more, the stress triaxiality distributes not uniformly across the minimum sectional area. Thus, stress triaxiality is a complicated function of equivalent strain and space distribution. In the fracture model parameter fitting process, the computing method for determining the values of specimen stress triaxiality η corresponding to the fracture strain plays a decisive role to the fitting result. However, Bridgman's formulation is not suitable for flat specimens either. For these reasons, a different technique is required to obtain a reasonable value of stress triaxiality. Bao and Wierzbicki [19] used a new concept of average stress triaxiality η av . The average stress triaxiality is interpreted by the integration of the stress triaxiality with respect to the equivalent plastic strain as shown in (Eq 17).
where t f is the time to fracture; ε p,eq is the equivalent plastic strain. The average stress triaxiality is the summation of stress triaxiality multiplied by the incremental equivalent plastic strain at t f . In this work, a method combining experimental data and numerical simulations was applied to obtain the stress triaxiality history along the entire experiment. The procedure is as follows: Step 1: Compute numerical simulations of the tested specimens with a suitable non-linear finite element code. Confirm that the engineering stress-strain curves obtained from FE simulation and tests are coincident to each other. It should be noted that a fracture criterion is not used.
Step 2: Identify the engineering failure strain by comparing numerical and experimental curves. Find the time step of specimen fracture.
Step 3: Obtain stress triaxiality and equivalent plastic strain histories in the most critical elements, which may be defined as those for which the values of stress triaxiality and equivalent plastic strain are most unfavorable. One should note that there is an assumption for the original J-C calibration [17], stress triaxiality and equivalent plastic strain distribute homogeneously across the minimum cross sectional area.
Step 4: Obtain the equivalent plastic strain at fracture for the most critical elements.
Step 5: Once the history of stress triaxiality will be obtained, the initial stress triaxiality can be read. Calculate the average stress triaxiality with (Eq 17). Then, two sets of stress triaxiality, i.e. the initial and average, are obtained.
Step 6: Plot the stress triaxiality versus the equivalent plastic strain at fracture. Formulate the relation between the stress triaxiality and the equivalent plastic strain at fracture. By using the above methodology, two sets of points in the stress triaxiality versus equivalent plastic strain space can be obtained. One of the sets is the lower fracture limit. And the other one is the proposed calibration point set. In fact, according to previous works [18][19][20][21], with the necking developing the stress triaxiality in the minimum section area is increasing. And high stress triaxiality contributes to early fracture. Therefore, the initial stress triaxiality is supposed to correspond to the lower fracture limit while the other set corresponds to the average stress triaxiality, which is closer to the actual situation. It should be emphasized that in the stress triaxiality versus plastic strain at fracture space, the initial stress triaxiality corresponds to the equivalent plastic strain at fracture calculated with (Eq 16) which is considered to be the average value of plastic strain at fracture along the minimum cross section. The reason of choosing the most unfavorable elements for obtaining the stress triaxiality and equivalent plastic strain histories, instead of the average stress triaxiality, is that fracture begins once some element has failed.

Stress-strain curves at different strain rates
The engineering stress-strain curves of AA6082-T6 under quasi tension are shown for 3 samples in Fig 6. In the beginning of tensile loading, the material shows linear elastic behavior within the elastic regime. When the stress reaches the yield limit (about 270 MPa) which means that the material comes into the plastic stage the variation of the stress tends to be much gentler till the stress meets ultimate strength. After yielding, further loading makes the material enter the plastic phase. A strain hardening effect can be seen in each case. Fig 7 shows true stress-strain curves of AA6082-T6 alloy at different strain rates under tension loading. The true stress and true strain are calculated using (Eq 2), which is valid up to necking. Curves in Fig 7 show that the ultimate strength appears in the final stage of deformation, which means that (Eq 2) is valid up to necking. The limit stresses are presented in Table 2. From a general view, no significant differences on the values of yield stresses and ultimate stresses can be observed between quasi-static tension tests and dynamic tests for this material during the plastic strain hardening phase. However, positive effects of strain rate on the yield stress can be clearly observed within the strain rate range from 0.001 s -1 to 1 s -1 . Experimental data for AA6082-T6 indicated that an increase in strain rate from 0.001 s -1 to 1 s -1 gives about a 2.5% increase in yield stress, while the ultimate strength rises from 351.51 MPa to 370.99 MPa with an increase of 5.5%. And it is worth noting that curves for quasi-static tension and strain rate 0.001 s -1 are overlapped by each other completely which means that the strain rate sensitivity is approximately zero in this range. True stress-strain curves from strain rate 0.001 s -1 to 1 s -1 are plotted in Fig 8. As we can see, experimental data of three repeated tests for each strain rate agrees very well, which excludes the effect of experimental spread. Thus, a conclusion can be drawn that AA6082-T6 exhibits a kind of relatively weak strain rate sensitivity within strain rate range from 0.001 s -1 to 1 s -1 while for high strain rate (above 800 s -1 ) strain rate have little influence on the material's plastic mechanical response. Based on the above point of view a strain rate related Johnson-Cook constitutive model is proposed in the sections below.

J-C constitutive model of AA6082-T6 alloy
The constants of the Johnson-Cook model for the investigated material are listed in Table 3. In (Eq 9), let the strain ε equal zero then the flow stresses were obtained at low strain rates and shown in Fig 9. The results calculated by the J-C model are in good agreement with the experimental values, and verify its effectiveness. In Fig 9, the stress vs. strain rate experimental data point at strain rate 0.001 s -1 is compared with the value of the J-C model as well, and the difference in stress is less than 2 MPa and small enough to neglect. Thus, for convenience the strain rate sensitivity range of the investigated material can be expanded to the range from 0.0001 s -1 Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities to 1 s -1 . The J-C model is used to predict the yield stresses under both low and high strain rates. The predicted results are listed in Table 4 and compared with the experimental results (Relative error ¼ s JÀ C y À s exp y s exp y Â 100%). It can be found that the yield stresses given by the J-C model at high strain rates are almost the same due to the small value of the rate sensitivity coefficient. And the three predicted values are all greater than the experimental results. This is because the loading process at a high strain rate increases the temperature and softens the specimens. The temperature increase can be computed with (Eq 6). The coefficient m can be determined with the least squares method. The parameters used in (Eq 6) and the value of m are listed in Table 5  Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities proposed J-C model at low strain rates. For high strain rate tension, the three predicted curves are almost completely overlapped due to the small value of rate sensitivity coefficient as well.
In general, the predicted curves can represent the actual plastic stress. In conclusion, the J-C model is able to describe the flow behavior of AA6082-T6.

Numerical simulations of smooth and notched specimens
In this paper, a non-linear finite element code, ABAQUS/Standard and ABAQUS/Explicit were both used for the numerical simulations under quasi-static and low strain rates respectively. Four smooth and notched specimen models were simulated. The J-C constitutive relation built above was used to define isotropic elastic-plastic material property. Other parameters of material property are referred to Table 5. It is notable that no failure criterion is introduced due to the calibration aim of simulation instead of validation.
Only the 1/4 of each of the specimens is modeled to reduce the computing time. The element type is the three dimensional hexahedral element with reduced integration (C3D8R). As with the experiments, axial velocity loads are applied at the end of the finite element models. The meshes of the finite element models are shown in Fig 12. Mesh refinement technique is adopted near the fracture zone. Fig 13 shows a typical finite element mesh of a notched specimen model with the notch radius of 10 mm. The mesh of the smooth specimen has 119120 elements while the specimens with notch radius R = 90 mm, R = 40 mm and R = 10 mm have 75400, 34000 and 54440 elements, respectively.  Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities The engineering stress-strain curves from numerical simulations agree well with those from experiments as shown in Fig 14. By simulating the processes of tensile tests, the spatial distribution of the stress triaxiality and the equivalent plastic strain from the beginning to the fracture can be obtained. Fig 15 shows stress triaxiality and equivalent plastic strain profiles of smooth and notched specimens along two perpendicular axes at the minimum cross section. Each of the profile is obtained from ABAQUS/Standard numerical simulations at the step time of fracture. The stress triaxiality values at the minimum cross section decrease from the center of the rectangular cross section to the edge. However, the stress triaxialities of the edges for different specimens are in a rough range from 0.33 to 0.4 where the difference is not as marked as the difference at central area. In the view of the space, the plastic deformation concentrates nearby the minimum section area, and the maximum value of the stress triaxiality is located in the center of the minimum section area almost throughout the whole plastic strain stage, as  expected. And the stress triaxiality values along two perpendicular axis show a small difference. The equivalent plastic strain shows the same tendency as the stress triaxiality. Nonetheless, as shown in Fig 15(b), the distribution of equivalent plastic strain along the transverse axis of 10 mm notch radius specimen exhibits an opposite behavior, which is opposite with the variation along the transverse axis too and may lead to an advanced failure. After all, due to the marked difference of the stress triaxiality level between the cross section center and the notch surface along the transverse axis for the 10 mm notch radius specimen, failure is supposed to begin at the center of the minimum cross section. In conclusion, the elements that failure begin with located in the center of the minimum cross sections of the specimens. The two sets of points in the equivalent plastic strain verse the stress triaxiality space are plotted in Fig 16. The first set of points given by a common method, in which the initial stress triaxiality (denoted as η in ) and the fracture strain calculated with (Eq 16) were adopted. All the fracture area ratios (A 0 /A f ) of minimum sectional area and the calculated fracture strain ε f eq are listed in Table 6 as well as the fracture strains from numerical simulations ε f À s eq , and the data  Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities are all averaged results of repeated tests. With the same height and width at the minimum cross section, a 0 and b 0 are both 20 mm for each specimen. And the second set is the averaged stress triaxiality versus fracture strain obtained from ABAQUS/Standard numerical simulation. Compared with the latter set points the points obtained by the common way incline to locate below, which implies a more conservative fracture criterion. Stress triaxiality histories of the critical elements are plotted as well. Curves in Fig 16 exhibit a generally increasing tendency with the increase of plastic strain for each type of the specimens. This tendency is more conspicuous for specimens with larger initial stress triaxialities.

Johnson-Cook fracture model constants identification
The construction of the J-C fracture model implies that five constants must be identified (see (Eq 10)). In present work, the effect of temperature is not taken into consideration. Thus, the constant in the temperature-dependent factor is neglected, then the formula, with four constants left to identify, is simplified as (Eq 9). The calibration of the J-C fracture model is performed using the data set of the average stress triaxiality versus equivalent plastic strain at fracture while the other J-C fracture model (the lower limit) is built using the data set of the initial average stress triaxiality versus equivalent plastic strain at fracture. Similar parameter fitting processes have been done by minimizing the residuals with success [5][6][7]. In the present work, the same method is adopted. The quasi-static tensile test data of smooth and notched specimens is used to fit the first three constants, D 1 , D 2 and D 3 . (Eq 9) Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities can be rewritten in such case just with the first factor: J-C fracture model parameters for the common method and calibration are listed in Table 7.

Dynamic mechanical response
Tests on aluminum alloy strain-rate sensitivity are reported by a large number of researchers. Many studies [3,22,23] have investigated the stress-strain behavior of AA6xxx and other aluminum alloys at a wide strain rate range. Compared with others AA6xxx aluminum alloys exhibit a rather low rate sensitivity, while AA7xxx alloys show moderate rate sensitivity. However, AA6005-T6 shows a rather strong rate sensitivity according to [24]. Actually, tests on AA6082-T6 aluminum alloy strain-rate sensitivity have been studied by Chen et al. [3]. In this paper, a series of tensile tests on AA6082-T6 smooth specimens in the strain rate range from quasi-static to 3400s -1 were conducted, and the stress-strain curves are plotted in Fig 7. A Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities similar result is obtained, that is, little rate sensitivity is observed on the mechanical response of AA6082. However, an additional phenomenon is observed that a slight but clear rate sensitivity exists in the strain rate range from quasi-static to 1s -1 , while the rate sensitivity is not found in dynamic tensile tests. The phenomenon is considered to be caused by the adiabatic heating of specimens at high strain rates. Such rate insensitivity or even negative sensitivity feature is an unfavorable property because it contributes to a faster fracture after onset of plastic deformation localization. Actually, plastic deformation is a process that concentrates in the localization zone which is a relatively small region. This concentration may increase the strain Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities rate locally. In general, increased strain rate implies high flow stress, stresses can thus redistribute in the necking zone. However, for a material featuring negative rate sensitivity, fracture may take place prematurely because the softening effect prevents the redistribution.

Validation of the J-C fracture model
In this paper, a calibration methodology is proposed to obtain a simplified J-C fracture model to describe the ductile fracture more accurately. Previous works have built J-C fracture models for various materials with success. A common method is used to obtain the data set of stress triaxialities and equivalent plastic strain at fracture that take into account the initial stress triaxiality (given by Bridgman's analysis (Eq 4)) and the equivalent plastic strain at fracture calculated with (Eq 14). However, this method may lead to a relatively conservative result due to the Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities underestimate of global stress triaxialities. For this reason, some calibrations for the J-C fracture model were proposed. Bao et al. [20] proposed to integrate the stress triaxiality of the entire process with respect to the equivalent plastic strain as shown in (Eq 15). Such a method is used by Choung et al. [25,26] to estimate the failure strain of EH36, and a good prediction was provided. Erice et al. [19] suggested an approach by averaging two data sets, the "lower" and the "upper" fracture limits. They took the data set from common method as the "lower"  Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities limit, and obtain the "upper" data set from numerical simulations where the stress triaxialities and the equivalent plastic fracture strains of the critical elements at the time step of fracture were used. Such calibration by merely averaging the two limits seems closer to reality than the common method, however, the accuracy may not be good enough. Here, the two methods were combined to calibrate the J-C fracture model. The average stress triaxiality is used and the equivalent plastic fracture strains are obtained from numerical simulations at the time step of fracture. Fig 17 shows the two models in the stress triaxiality versus equivalent plastic strain space. What is noteworthy is that the relation is given only in the region of high stress triaxialities due to the limitation of J-C fracture model [20]. The calibrated J-C fracture model was used in the numerical simulations of the quasi-static tensile tests for smooth and notched specimens. Also, the J-C fracture model was used in ABA-QUS/Explicit simulations for comparison. Fig 18 shows the fracture process of the smooth specimen in the numerical simulation. The fracture simulation gives a depiction of the stress triaxiality distribution and fracture process. When the equivalent plastic strain reaches the fracture strain elements are deleted automatically. It can be seen that failure started in the Table 7. J-C fracture model parameters for the common method and calibration.
Parameter  Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities to validate the calibration of the J-C fracture model. The common J-C fracture model shows a relatively "conservative" estimating on the ductility of AA6082-T6 alloy as expected. The underestimation reaches to an averaging extent of 13.50%. And the average error of the engineering fracture strain after calibration is 6.05% which is a great improvement on fracture strain estimation. The results give an evidence that the calibration methodology can improve the accuracy of the J-C fracture model effectively.

Charpy impact test
Railway collision accidents are rare, but once it happened the consequences could be catastrophic. In recent years, researchers such as Peng et al. [27] have investigated the collision performance of train energy absorbing structures. AA6082-T6, as a train body manufacturing material is supposed to be of good impact resistance. To further investigate its dynamic impact fracture property, two groups of Charpy impact tests are carried out. Two groups of specimens with the same size but different notch depth are machined. The dimension is 55×10×10 mm as shown in Fig 20 (unit: mm). The pendulum impact testing machine provides 150 J impact energy by lifting the pendulum to a height, then released to impact the specimen placed horizontally bellow the pendulum. During the impact, part of the mechanical energy is dissipated by specimen deformation and fracture. The height of pendulum after impact is recorded, and the energy difference before and after the impact can thus be calculated.
Two finite element models corresponding to the two tests are established based on ABA-QUS/Explicit to verify the effectiveness of the calibrated J-C fracture model. The material property is defined same as the models built above. Both of the common and calibrated J-C fracture model are used to simulate the tests. The models before and after the impact are shown in Fig 21. The deformation shape in simulation qualitatively agrees with the experiment results. On ther other hand, the Charpy impact energy obtained directly from experiment and that from the simulation are compared in Table 8. The energy absorption differences between Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities  simulation and experiment are 16.34% and 24.38% using the common J-C fracture model for specimens with the notch depth of 2 mm and 5 mm, respectively, which are much higher than them (4.91% and 11.14%) using the calibrated J-C fracture model. It indicates that fracture behavior of AA6082-T6 alloy can be more effectively described through the calibrated J-C failure model.

Fractography analysis
From the results of a series of tensile tests, it is found that AA6082-T6 exhibits an excellent ductility, with an elongation at fracture of 10% to 25%. A number of studies have been done on ductility damage and ductility fracture [28][29][30][31]. It is common view that the dimple fracture (experiencing the process of nucleation, growth, and coalescence of voids) and shear fracture (developed from shear band localization) are two primary fracture forms for ductile materials. Fig 22 shows the smooth and notched specimens after fracture at quasi-static strain rate. Typical cup-and-cone form fracture surfaces are shown. The central zone of the fracture surface shows ductile fracture full of bumps and hollows while the residuary area represents shear fracture with the surface along the maximum shear stress plane, at an angle of 45 degree to the tensile axis. SEM fractographs of the smooth specimen fracture surface are shown in Fig 23. Fig  23(a) shows a typical nucleation-growth-coalescence process of ductile fracture, where a lot of Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities voids and dimples can be observed. In the central zone, the stress triaxiality is supposed to be higher than at the area close to the specimen surface. When a ductile material is subjected to hydro static stress, micro-voids nucleate in the central zone giving priority to impurities and second-phase particles, then grow bigger with the development of deformation, at last coalescence of grown void leads to many micro cracks. When the micro cracks link up and expand to the area close to the specimen surface shear fracture takes place. Fig 23(b) shows shear fracture mode. SEM micrographs of the fracture surfaces of smooth and notched specimens tested at quasi-static strain rates are compared in Fig 24. It is found that dimples had various sizes and heights in the central zone of fracture surface. The reason is that in the nucleation stage, voids grow rapidly due to high stress triaxiality. High hydrostatic pressure likely results in large dimples which may accelerate the dimple fracture. And low hydrostatic pressure leads to shear band localization and further causes shear fracture.

Conclusions
In this paper, a series of tests were carried out to investigate the mechanical properties of AA6082-T6 under various strain rates and different stress states. Numerical simulations based on FEA were successfully performed. Primary conclusions can be drawn as follows: • From a general view, strain rate has little influence on flow stress. However, clear rate sensitivity was observed in the strain rate range of 0.001 s -1 to 1 s -1 . The softening effect of adiabatic heating counteracts the positive strain rate sensitivity at high strain rates. A Johnson-Cook constitutive model was thus proposed. Flow and fracture behavior of aluminum alloy 6082-T6 at different tensile strain rates and triaxialities • A calibration methodology for the J-C fracture model was proposed by using the average stress triaxiality and the equivalent plastic fracture strain of the critical elements (in the central zone of the minimum cross section).
• Numerical simulations based on FEA were performed using ABAQUS. By simulating the tensile process, it was proven that the stress triaxiality is highest along the fracture cross section in the central zone, where failure begins first. The calibrated J-C fracture model was proven to be valid in predicting the fracture strain of AA6082-T6 alloy.
• Ductile fracture is the main fracture form in tensile test for the present material. Dimple fracture appears in the center of the minimum cross section due to the high level of stress triaxiality and the remaining cross sectional area is sheared to fracture. This may be caused by the rapid growth of voids in the nucleated stage due to a high stress triaxiality.