Optimized Calibration Procedure of the Damage Parameters of 6082-T6 Sheets

The constitutive equation of AA6082-T6 sheets is investigated by an extended Rousselier damage model. Previous research is mainly comprised of single-pull test specimens, notched tensile specimens, tensile specimens with a hole, and shear specimens. To decrease the natural material errors, a new specimen is used. In this paper, a clinched joint is regarded as a specimen to calibrate the initial void volume fraction kω and the shear damage parameter f0 by using the orthogonal analysis method, which can reduce the simulation times and accuracy. It also reveals that the initial void volume fraction f0 affects the void volume fraction in the neck and the groove of the clinched joint, and the shear damage parameter kω affects just the void volume fraction f in the neck of the clinched joint. It checks the force-displacement curve, shape of the clinched joint, and the fracture location, and approves that these damage parameters can describe the deformation process, fracture location, and shape of the clinched joint.


Introduction
The request to reduce fuel consumption, emissions, and weight in the automotive industry has previously led to the search for a new joining technique and material capability. The mechanical durability of the joints depends on several parameters, including the specific joining technique, the nature of assembled materials, and the corrosion stability of the constituents. AA6082-T6 sheets are widely used in the automotive industry as a basis for lightweight alternatives to classical designs with mild and high-strength-low alloyed steels. Dadbakhsh [1] explored the AA6082-T6 sheet and studied the alloy strengthening because of pre-and post-ECAP (Equal channel angular pressing) aging treatment. It showed that aging before and after ECAP processing was an effective method for strengthening of the alloy. Chen [2] suggested that the stress-strain behavior of extruded AA6xxx and AA7xxx aluminum alloys in T6 temper should be studied at a wide range of strain rates. The AA6xxx alloys were found to exhibit an insignificant rate sensitivity in the stress-strain behavior.
Some papers focused on the deformation of the alloy sheet to discuss the application in some fields, such as the deep drawing, cold welding process, tensile tension, and mechanical connection. A great deal of the experimental research will greatly increase the cost of the production, and the selection of a suitable constitutive equation to obtain the parameters is impossible. The FE (finite element) simulation analysis is an essential method in studying the deformation process, as it can obtain some parameters without any experimentation. The success of the simulation model depends upon the capacity of the constitutive models and respective material parameters to accurately reproduce the T6 are also constructed by the MMC3, DF2012, and DF2014 criteria for the purpose of comparison. These criteria are implemented into numerical prediction of ductile fracture in shear, uniaxial tension, plane strain tension, and balanced biaxial tension.
These damage models were applied to describe the deform process of metal material, such as deep drawing, cold welding process, tensile tension, and mechanical connection. Commonly, it is difficult to obtain a uniform constitutive equation for all materials. At present, the modified GTN and modified Rousselier damage models also introduced the shear damage parameter to describe the failure behavior. However, previous studies cannot describe the fracture in the groove of the clinched joint. In Xu's [8] paper, a clinched joint should be divided into four regions, including the no deforming region, bending region, compress region, and tensile region. In order to obtain accurate damage parameters, the clinched joint was regarded as a specimen to calibrate the damage parameters.
The present study uses data from the uniaxial tensile specimen, and the clinched joint test is also used to identify the intrinsic uniaxial stress-strain behavior, the initial void volume fraction, and the shear damage parameter. The paper is organized and closely follows the steps in the calibration agreement. In Section 2, the determination of the hardening behavior from the tensile tests is made. In Section 3, the determination of the initial void volume fraction and the shear damage parameter from the clinched joint by the orthogonal test is made. Simultaneously, a validation of the shear damage parameter and the initial void volume fraction of the extended Rousselier damage model [9] is performed. In Section 4, the applicability of the calibrated constitutive model to the clinched joint failure mode id discussed as one illustration. Possible variations in the identification agreement for other ductile materials are also discussed.

Calibration of the Hardening Behavior
The AA6082-T6 sheets were cut into tensile test specimens. The samples for the tensile test were 200 mm long and 25 mm wide, as shown in Figure 1. The mesh size in the densest area is 100 µm × 100 µm. The mechanical response of the undamaged materials was obtained from quasi-static uniaxial tensile tests on the specimen ( Figure 1) with elastic-plastic finite element computations. The test geometry and finite elements are shown in Figure 1. AA6082-T6 with a nominal chemical composition in Table 1 was the candidate for the micro-indentation test to inspect the fracture toughness. The micro-indentation specimens were carefully ground with sandpaper. The specimens were polished with 5th and 1th grit diamonds to a mirror finish. The specimen tests were conducted under quasi-static loading, whereas all simulations were performed with the dynamic code in ABAQUS/Explicit [9]. To minimize the inertial effects and efficiently simulate the quasi-static tests in the explicit code, a preliminary series of calculations with different fixed and applied loading rates was performed for each test configuration. At some loading rates, the simulations converge to a quasi-static limit as the rates decrease. The loading rate was then employed in all subsequent calculations. purpose of comparison. These criteria are implemented into numerical prediction of ductile fracture in shear, uniaxial tension, plane strain tension, and balanced biaxial tension. These damage models were applied to describe the deform process of metal material, such as deep drawing, cold welding process, tensile tension, and mechanical connection. Commonly, it is difficult to obtain a uniform constitutive equation for all materials. At present, the modified GTN and modified Rousselier damage models also introduced the shear damage parameter to describe the failure behavior. However, previous studies cannot describe the fracture in the groove of the clinched joint. In Xu's [8] paper, a clinched joint should be divided into four regions, including the no deforming region, bending region, compress region, and tensile region. In order to obtain accurate damage parameters, the clinched joint was regarded as a specimen to calibrate the damage parameters.
The present study uses data from the uniaxial tensile specimen, and the clinched joint test is also used to identify the intrinsic uniaxial stress-strain behavior, the initial void volume fraction, and the shear damage parameter. The paper is organized and closely follows the steps in the calibration agreement. In Section 2, the determination of the hardening behavior from the tensile tests is made. In Section 3, the determination of the initial void volume fraction and the shear damage parameter from the clinched joint by the orthogonal test is made. Simultaneously, a validation of the shear damage parameter and the initial void volume fraction of the extended Rousselier damage model [9] is performed. In Section 4, the applicability of the calibrated constitutive model to the clinched joint failure mode id discussed as one illustration. Possible variations in the identification agreement for other ductile materials are also discussed.

Calibration of the Hardening Behavior
The AA6082-T6 sheets were cut into tensile test specimens. The samples for the tensile test were 200 mm long and 25 mm wide, as shown in Figure 1. The mesh size in the densest area is 100 μm × 100 μm. The mechanical response of the undamaged materials was obtained from quasi-static uniaxial tensile tests on the specimen ( Figure 1) with elastic-plastic finite element computations. The test geometry and finite elements are shown in Figure 1. AA6082-T6 with a nominal chemical composition in Table 1 was the candidate for the micro-indentation test to inspect the fracture toughness. The micro-indentation specimens were carefully ground with sandpaper. The specimens were polished with 5th and 1th grit diamonds to a mirror finish. The specimen tests were conducted under quasi-static loading, whereas all simulations were performed with the dynamic code in ABAQUS/Explicit [9]. To minimize the inertial effects and efficiently simulate the quasi-static tests in the explicit code, a preliminary series of calculations with different fixed and applied loading rates was performed for each test configuration. At some loading rates, the simulations converge to a quasi-static limit as the rates decrease. The loading rate was then employed in all subsequent calculations.   In the paper by Luo [19], the experimental results show a strong dependency of the strain to fracture the material orientation with respect to the loading direction. The use of the isotropic modified Mohr-Coulomb stress state weighting function in this anisotropic fracture modeling framework provides accurate predictions of the onset of fracture for all thirteen fracture experiments. Material strain-rate dependence is ignored in the present computations. Each test was conducted five times, and the average values were calculated to eliminate the error.
According to the tensile test, the engineering stress and strain curve can be obtained. Simultaneously, the true stress and true strain can be obtained with the equations. In Figure 2, better hardening behavior was shown for a better long space; therefore, the stress strain curve for a long space is difficult to check with the tensile test. In this paper, the Voce equation is used to fit the performance of the AA6082-T6 sheet with a thickness of 2 mm. The simulation results and the experiment results of the stress-strain curves are shown in Figure 2. Basically, prior to necking, the true strain is given by ε T = ln(1 + ε N ), and the true stress is given by σ T = σ N (1 + ε N ). To deduce the true response in the post-necking regime, computations were performed with an assumed form of the stress-strain relation and matching the predicted force-displacement curves with those obtained experimentally. In the paper by Luo [19], the experimental results show a strong dependency of the strain to fracture the material orientation with respect to the loading direction. The use of the isotropic modified Mohr-Coulomb stress state weighting function in this anisotropic fracture modeling framework provides accurate predictions of the onset of fracture for all thirteen fracture experiments. Material strain-rate dependence is ignored in the present computations. Each test was conducted five times, and the average values were calculated to eliminate the error.
According to the tensile test, the engineering stress and strain curve can be obtained. Simultaneously, the true stress and true strain can be obtained with the equations. In Figure 2, better hardening behavior was shown for a better long space; therefore, the stress strain curve for a long space is difficult to check with the tensile test. In this paper, the Voce equation is used to fit the performance of the AA6082-T6 sheet with a thickness of 2 mm. The simulation results and the experiment results of the stress-strain curves are shown in Figure 2. Basically, prior to necking, the true strain is given by T N ε ε , and the true stress is given by (1 ) To deduce the true response in the post-necking regime, computations were performed with an assumed form of the stress-strain relation and matching the predicted force-displacement curves with those obtained experimentally. As shown in Figure 2, this hardening equation cannot accurately describe the hardening process of the AA6082-T6 sheet. The VOCE hardening equation can describe the part before the uniform deformation, and the part after the uniform deformation needs a new hardening equation to fit. Consequently, in this paper, the VOCE hardening behavior is modified. Before the uniform deformation, the hardening model VOCE was used. After the uniform deformation, the shear band was shown, and the normal hardening model cannot be used. However, the material sheet should still be hardened, and the hardening curve is unknown. Luo [20] discussed an optimized method comprehension for all experimental results and adjusted the hardening curve to save the accordance between the simulation and experimental results according to the displacement and the force curve. Sam [21] presented a method to obtain the hardening equation with DIC (Digital Image Correlation). According to Sam's method for cold and rolled interstitial-free steel sheet, we can compare the results of experiments and simulation by PNTE (Post-NeckingTensile Experiment) -SWIFT and PNTE-VOCE. The tensile process of PNTE-SWIFT hardening behavior is regarded as the hardening equation of the Al6082-T6 sheet with 2 mm thickness. Thus, an important subject As shown in Figure 2, this hardening equation cannot accurately describe the hardening process of the AA6082-T6 sheet. The VOCE hardening equation can describe the part before the uniform deformation, and the part after the uniform deformation needs a new hardening equation to fit. Consequently, in this paper, the VOCE hardening behavior is modified. Before the uniform deformation, the hardening model VOCE was used. After the uniform deformation, the shear band was shown, and the normal hardening model cannot be used. However, the material sheet should still be hardened, and the hardening curve is unknown. Coppieters et al. [20] presented a method to obtain the hardening equation with Digital Image Correlation (DIC). Luo [21] discussed an optimized method comprehension for all experimental results and adjusted the hardening curve to save the accordance between the simulation and experimental results according to the displacement and the force curve.
This study agreed with the paper by Xue [22]. The material was unstable when the strain value was 0.14. Actually, the maximum equivalent strain should be 0.1522 when the deformation of material reaches instability. Thus, the two hardening equations can be described in the AA6082-T6 sheet: the hardening behavior and the mixed hardening behavior: where A = 402.05, B = 292.56, C = 12.26, σ peak = 385.5326, and ε peak = 0.1522.
A preliminary estimate of the strain hardening exponent n was obtained by fitting. A series of finite element computations was performed to ascertain the full nominal tensile stress-strain curve with a range of n values. As shown in Figure 3, the results for n = 0.18 accurately replicate the experimental measurements up to the onset of rupture. When the strain is smaller 0.1522, the voce hardening equation was used. When the strain is larger than 0.1522, Xue's equation was used to describe the hardening behavior. In other words, Xue's equation will be used after instability. must be found to a suitable hardening behavior to take the place of PNTE-VOCE and PNTE-SWIFT. This study agreed with the paper by Xue [22]. The material was unstable when the strain value was 0.14. Actually, the maximum equivalent strain should be 0.1522 when the deformation of material reaches instability. Thus, the two hardening equations can be described in the AA6082-T6 sheet: the hardening behavior and the mixed hardening behavior: 402.05 (402.05 292.56) exp ( 12.26 ) 0.1522 A preliminary estimate of the strain hardening exponent n was obtained by fitting. A series of finite element computations was performed to ascertain the full nominal tensile stress-strain curve with a range of n values. As shown in Figure 3, the results for n = 0.18 accurately replicate the experimental measurements up to the onset of rupture. When the strain is smaller 0.1522, the voce hardening equation was used. When the strain is larger than 0.1522, Xue's equation was used to describe the hardening behavior. In other words, Xue's equation will be used after instability. All the simulation results show that the mixed hardening equation can accurately be described by the hardening process of the AA6082-T6 sheet with 2 mm when the hardening index n is equal to 0.18. Other than 0 f and w k , the basic parameters characterizing the constitutive model that are used in all simulation in this paper are the following: In Guo's paper [14], the yield potential is written as Equation (3), so it is a coupled constitutive equation in which the damage accumulation and hydrostatic stress are incorporated: where m σ is hydrostatic stress, All the simulation results show that the mixed hardening equation can accurately be described by the hardening process of the AA6082-T6 sheet with 2 mm when the hardening index n is equal to 0.18. Other than f 0 and k ω , the basic parameters characterizing the constitutive model that are used in all simulation in this paper are the following: In Guo's paper [14], the yield potential is written as Equation (3), so it is a coupled constitutive equation in which the damage accumulation and hydrostatic stress are incorporated: where σ m is hydrostatic stress, σ eq = (3σ d : σ d /2) 1/2 is the von Mises equivalent stress, ρ = (1 − f )/(1 − f 0 ) is the relative density, f is the damage variable or void volume fraction, f 0 is the initial is the equivalent plastic strain, D and σ 1 are Rousselier material constants: Here, k ω is the shear damage coefficient, which sets the magnitude of the void coalescence rate in shear deformation on the basis of Junhang Guo et al. [14]. And the invariant measure ω(σ) is given by: where J 3 is the third stress invariant of the deviatoric stress tensor σ d (σ d is the deviatoric stress tensor); σ 1 , σ 2 , and σ 3 are the principal stresses of the stress tensor σ and are assumed to be ordered as The non-dimensional metric in Equation (5) lies in the range 0 ≤ ω ≤ 1 to discriminate between axisymmetric and shear-dominated stress states. For all axisymmetric stress states, ω = 0. For all states comprised of a pure shear stress plus a hydrostatic contribution, ω = 1. With f 0 as the initial void volume fraction, the analytical solution can be derived as: Then for the shear stress state, the solution can be particularized with D = 2 (Junhang Guo [14]), and ω(σ) = 1 as: From Equation (8), it can be seen that the damage coefficients f 0 and k ω need be calibrated. In this paper, the damage coefficients of the material were identified by an inverse method. It described the clinching process to calibrate the extended Rousselier damage model.

3D Model
Based on the previous discussion, the use of the failure characteristics during the clinching process calibrated the shear damage parameter and the initial void volume fraction. In the first section, a new specimen was proposed to optimize the calibration of damage parameters. Figure 4 shows a 3D model of the experimental sample with the groove of the clinched joint and in the neck of the clinched joint. Figure 4a is the 3D model based on the paper by Xu [8]. The mesh was divided into three parts: the intensive, transition, and loose regions. The symmetric boundary condition is used and one-quarter geometry is modeled. To improve the computational efficiency, only the material in the clinched joint region is modeled using the modified Rousselier model. Outside this zone the plate is modeled using the standard von Mises plasticity model provided by ABAQUS/Explicit. The minimum element size is 100 µm × 100 µm × 100 µm. Adaptive meshing is employed to avoid element distortion in the large localized shear deformation in the simulation. Figure 4b shows the tools and the clinched joint, wherein the die is closed. Figure 4c shows a specimen where the fracture just happened in the groove of the clinched joint, but in the neck of the clinched joint. Figure 4d shows that the cross-section of the clinched joint, including all geometrical parameters of the clinched joint, is simultaneously used to find the obvious fracture that happened in the groove and neck of the clinched joint. clinched joint. Figure 4d shows that the cross-section of the clinched joint, including all geometrical parameters of the clinched joint, is simultaneously used to find the obvious fracture that happened in the groove and neck of the clinched joint. The fracture obviously happened in the groove of the clinched joint but not in the neck of the said joint in Figure 4. The tensile stress causes a crack in the groove of the clinched joint. In Xue's paper [7], the groove region belongs to the tensile deformation region. In the previous discussion, Xu [8] proposed that the shear stress and tensile stress cause the crack in the neck of the clinched joint. If the crack happens in the neck of the clinched joint, the shear damage parameter should be controlled. During the tensile test process, the reduced thickness of the sample (12%) is smaller than the thickness (84%) of the neck of the clinched joint. When the bottom thickness reached at 0.56 mm, the reduced ratio thickness is 86%. The crack happened in the groove of the clinched joint. Local tensile deformation leads to the crack that happens in the groove of the clinched joint. The local deformation region reached the yield limitation, and then the crack starts at this point. Therefore, the fracture can be controlled in the groove.

Orthogonal Analysis
To reduce the fitting time and evaluate the computation results, an effective design method is needed [23]. Several different factors and levels need to be considered when optimizing the parameters during the clinching joint. The purpose of the optimization of parameters during the clinched joint is to get the damage parameters of the Rousselier damage model. Orthogonal design is a kind of design method that is mainly used to study multiple factors and multiple levels. This design method is uniformly dispersed, neat, and comparable, making each design highly representative. The selection of the representative points from the comprehensive design can fully reflect the impact of different levels of each factor on the design result that takes less time. Oudjene and Ben-Ayed [24] used the Taguchi's experiment design method that is used to investigate the effects of geometry tools on the clinch joint resistance as well as on its shape. Chen [25] designed a plan to optimize the shape of the rivet to reduce the protrusion height and increase the strength of clinched joint. Wen [26] presented a response surface optimization of the clinching tools.
During the clinched process, the clinched joints with fracture zone in the groove of the clinched joint are obtained. The previous research mainly discussed how to get some qualified samples to avoid the failure samples. In the discussion, some ideas are discussed. A lot of experiments show that the surface stage will affect the quality of the clinched joint; thus, in our research, the grease is used to avoid the damage that comes from the surface stage. The friction between tools and sheets slightly influences the intensity of the punch force and the distribution of deformation between the upper and lower sheets. This tendency is even greater with respect to metal hardening.
Since the fracture just happened, in the groove of the clinched joint during the clinched process, it must be ensured that there will be no cracks in the neck of the clinched joint. In the paper by Xu [8], the clinched joint is divided into four parts, including no deformation zone, tensile deformation The fracture obviously happened in the groove of the clinched joint but not in the neck of the said joint in Figure 4. The tensile stress causes a crack in the groove of the clinched joint. In Xue's paper [7], the groove region belongs to the tensile deformation region. In the previous discussion, Xu [8] proposed that the shear stress and tensile stress cause the crack in the neck of the clinched joint. If the crack happens in the neck of the clinched joint, the shear damage parameter should be controlled. During the tensile test process, the reduced thickness of the sample (12%) is smaller than the thickness (84%) of the neck of the clinched joint. When the bottom thickness reached at 0.56 mm, the reduced ratio thickness is 86%. The crack happened in the groove of the clinched joint. Local tensile deformation leads to the crack that happens in the groove of the clinched joint. The local deformation region reached the yield limitation, and then the crack starts at this point. Therefore, the fracture can be controlled in the groove.

Orthogonal Analysis
To reduce the fitting time and evaluate the computation results, an effective design method is needed [23]. Several different factors and levels need to be considered when optimizing the parameters during the clinching joint. The purpose of the optimization of parameters during the clinched joint is to get the damage parameters of the Rousselier damage model. Orthogonal design is a kind of design method that is mainly used to study multiple factors and multiple levels. This design method is uniformly dispersed, neat, and comparable, making each design highly representative. The selection of the representative points from the comprehensive design can fully reflect the impact of different levels of each factor on the design result that takes less time. Oudjene and Ben-Ayed [24] used the Taguchi's experiment design method that is used to investigate the effects of geometry tools on the clinch joint resistance as well as on its shape. Chen [25] designed a plan to optimize the shape of the rivet to reduce the protrusion height and increase the strength of clinched joint. Wen [26] presented a response surface optimization of the clinching tools.
During the clinched process, the clinched joints with fracture zone in the groove of the clinched joint are obtained. The previous research mainly discussed how to get some qualified samples to avoid the failure samples. In the discussion, some ideas are discussed. A lot of experiments show that the surface stage will affect the quality of the clinched joint; thus, in our research, the grease is used to avoid the damage that comes from the surface stage. The friction between tools and sheets slightly influences the intensity of the punch force and the distribution of deformation between the upper and lower sheets. This tendency is even greater with respect to metal hardening.
Since the fracture just happened, in the groove of the clinched joint during the clinched process, it must be ensured that there will be no cracks in the neck of the clinched joint. In the paper by Xu [8], the clinched joint is divided into four parts, including no deformation zone, tensile deformation zone, shear deformation, and bending deformation zone. Zhao [9] presented that the shear damage will happen in the neck of the clinching joint. Therefore, the assumption is that the shear deformation region check the shear damage parameter k ω and the tensile deformation region and validate the initial void volume fraction f 0 has been made. In other words, the shear damage parameter k ω and the initial void volume fraction f 0 affect the quality of the clinched joint. The purpose functions are f 1 , f 2 , and F. The function f 1 expresses the void volume fraction when crack happens in the neck of the clinched joint, whereas f 2 expresses the void volume fraction when crack happens in the groove of the clinched joint. The function F expresses the deformed force from the experimental results when the bottom thickness of the clinched joint is 0.56 mm.
Several factors and levels have been investigated. In order to facilitate research, it need to construct a suitable orthogonal table [23]. In Section 2, it introduces the function of the void volume fraction and discusses the relationship between the shear damage parameter and the initial void volume fraction. It has explained that the shear damage parameter and the initial void volume fraction, they are independent parameters. In this paper, the relationship between shear damage parameter and the initial void volume fraction is ignored. The suitable level range will reduce the number of optimization design. Therefore, before the orthogonal analysis, it is necessary to predict the shear damage parameter and the initial void volume fraction, that is, use the specimen shown in Figure 1 to determine a rough range. Finally, six five-level factors are contained in the orthogonal array used for the optimization of parameters during the clinched process as presented in Table 2 (L25 (5 6 )). These factors are more helpful in getting the suitable damage parameters during the clinching process. Commonly, lubricant is used between punch and upper sheet, so µ 1 is equal to 0. The critical void volume fraction is equal to 0.2 [27]. In fact, four parameters, µ 2 , µ 3 , k ω , and f 0 need be calibrated by this method. µ 1 and f c are regarded as fixed constants. Then a group of 25 different parameter combinations should be simulated to find the optimal the damage parameters and the state parameters according to the orthogonal array as shown in Table 3A-F represent the friction coefficient between punch and the upper sheet µ 1 , the friction coefficient between the upper sheet and the under sheet µ 2 , the friction coefficient between the die and the under sheet µ 3 , the shear damage parameter k ω , the initial void volume fraction f 0 , the critical void volume fraction f c , respectively.  The extent of the effect in the void volume fraction during the clinching process should be investigate; in other words, a method should be developed on how to get the crack that happens in the groove of the clinched joint. The range analysis was made according to the results of the numerical simulation. The data of the range analysis are shown in Tables 4-6. The effect curves of the parameters are shown in Figures 5-7.  The most important factor that affects the fracture location in the neck of the clinching joint is shown in Table 4. K a expresses the sum of test indicators (a expresses level (1, 2, 3, 4, 5)), i expresses level (1,2,3,4,5), j expresses factor (A, B, C, D, E, F), k ij is equal to K ij divided by 5, rank is equal to the maximum (k ij ) minus the minimum (k ij ). Table 4 shows thatextreme difference analysis for the void volume fraction on the neck of the clinched joint. Factor D has an influence on the void volume fraction f on the neck of the clinched joint. D expresses the shear damage parameter k ω . The previous research shows that the shear stress leads to the failure on the neck of the clinched joint. Thus, the orthogonal analysis method can describe the trend. Table 4 shows that the shear damage parameter k ω has a great influence on the void volume on the neck of the clinched joint; simultaneously, the initial void volume fraction f 0 also influences the void volume fraction f on the neck of the clinched joint. The shear damage parameter k ω is the one of the most important factors on the neck of the clinched joint. Except on the shear damage parameter k ω and the initial void volume fraction f 0 , the shear damage parameter k ω also has an impact on other factors. Table 5 shows that the initial void volume f 0 has a signification influence on the void volume fraction f on the groove of the clinched joint. Table 6 shows that the coefficient friction µ 2 between the two sheets is one of the most important factors on the force F. The extent in which the factors affect the deform force F is shown in Figure 6, and the effect curves of the factors are also shown in Figure 6. The final optimized results cannot be obtained through extreme difference analysis because the orthogonal analysis' maximum or minimum value is not a subject value. The orthogonal analysis was used to get the influence trend of the void volume fraction on the different deformation regions (neck and grove regions). The initial void volume fraction f 0 and the shear damage parameter k ω also obviously affect the void volume fraction f on the neck of the clinched joint ( Figure 5 and Table 4). During the simulation process, the initial void volume fraction f 0 and the shear damage parameter k ω must be controlled to ensure that there will be an absence of cracks in the neck of the clinched joint. Apparently, the initial void volume fraction f 0 affected the void volume fraction f distribution in the groove of the clinched joint as shown in Figure 6 and Table 5. Simultaneously, the factor that influences the deform force is known.   Unlike other studies, in this paper, the orthogonal analysis method was used, but the subject function cannot be considered as the final result. The maximum or minimum value can evaluate the influence trend, but it cannot be considered as the final result. The aim of the paper is to discuss the calibration of the damage parameters for AA6082-T6 sheet with a thickness of 2 mm by the clinching process according to the fracture location. After 25 simulation results, it is easy to find that the simulation of No.10 is close to the experimental result. However, as shown in Figure 8a, no fracture   Unlike other studies, in this paper, the orthogonal analysis method was used, but the subject function cannot be considered as the final result. The maximum or minimum value can evaluate the influence trend, but it cannot be considered as the final result. The aim of the paper is to discuss the calibration of the damage parameters for AA6082-T6 sheet with a thickness of 2 mm by the clinching process according to the fracture location. After 25 simulation results, it is easy to find that the simulation of No.10 is close to the experimental result. However, as shown in Figure 8a, no fracture   Unlike other studies, in this paper, the orthogonal analysis method was used, but the subject function cannot be considered as the final result. The maximum or minimum value can evaluate the influence trend, but it cannot be considered as the final result. The aim of the paper is to discuss the calibration of the damage parameters for AA6082-T6 sheet with a thickness of 2 mm by the clinching process according to the fracture location. After 25 simulation results, it is easy to find that the Unlike other studies, in this paper, the orthogonal analysis method was used, but the subject function cannot be considered as the final result. The maximum or minimum value can evaluate the influence trend, but it cannot be considered as the final result. The aim of the paper is to discuss the calibration of the damage parameters for AA6082-T6 sheet with a thickness of 2 mm by the clinching process according to the fracture location. After 25 simulation results, it is easy to find that the simulation of No.10 is close to the experimental result. However, as shown in Figure 8a, no fracture happens in the neck and the groove of the clinched joint. Apparently, the void volume fraction f in the groove of the clinched joint is close to 0.2, whereas the void volume fraction f in the neck of the clinched joint is almost close to 0.1. The subject stage is the fracture in the groove of the clinched joint; thus, according to the above analysis, an investigation must be conducted if the initial void volume fraction meets the fracture location. The second optimized result was shown in Table 7. Compared with the two group parameters as shown in Figure 8, the fracture obviously happens in the groove of the clinched joint, but in the neck of the clinched joint, as shown in Figure 8b,there was no crack that was observed anywhere because the void volume fraction f is smaller than the critical void volume fraction c f when the initial void volume fraction 0 f is 0.009. Evolution of the damage at the neck of the clinched joint, shear damage parameter ω k , and the initial void volume fraction 0 f lead to the damage. In the groove of the clinched joint, the initial void volume fraction 0 f leads to the fracture. If the shear damage parameter ω k was not discussed, the shear stress in the neck of the clinched joint will be neglected. That means an important factor can be ignored, and then there is a greater difference between the simulation results and the experiment results. It was observed that the void volume fraction of tracked nodes normally ranged from 0.001 to 0.01. The shear damage parameter ω k normally ranged from 0.5 to 1. The simulation result for the damage parameter of the initial void volume fraction 0 f distribution on the clinched joint is shown in Figure 9. The maximum value of SDV16 (the void volume fraction) is shown in the groove of the clinched joint. Compared with the experimental results, the fracture that happens in the groove of the clinched joint is obviously noted. Thus, this model is known to describe the failure location; in other words, these damage parameters can accurately describe the performance of the AA6082 sheet with a thickness of 2 mm during the clinched process.  Compared with the two group parameters as shown in Figure 8, the fracture obviously happens in the groove of the clinched joint, but in the neck of the clinched joint, as shown in Figure 8b, there was no crack that was observed anywhere because the void volume fraction f is smaller than the critical void volume fraction f c when the initial void volume fraction f 0 is 0.009. Evolution of the damage at the neck of the clinched joint, shear damage parameter k ω , and the initial void volume fraction f 0 lead to the damage. In the groove of the clinched joint, the initial void volume fraction f 0 leads to the fracture. If the shear damage parameter k ω was not discussed, the shear stress in the neck of the clinched joint will be neglected. That means an important factor can be ignored, and then there is a greater difference between the simulation results and the experiment results. It was observed that the void volume fraction of tracked nodes normally ranged from 0.001 to 0.01. The shear damage parameter k ω normally ranged from 0.5 to 1.
The simulation result for the damage parameter of the initial void volume fraction f 0 distribution on the clinched joint is shown in Figure 9. The maximum value of SDV16 (the void volume fraction) is shown in the groove of the clinched joint. Compared with the experimental results, the fracture that happens in the groove of the clinched joint is obviously noted. Thus, this model is known to describe the failure location; in other words, these damage parameters can accurately describe the performance of the AA6082 sheet with a thickness of 2 mm during the clinched process.  Figure 10 shows that the displacement-force curve and the shape of the clinched joint is shown in Figure 10. The comparison experiment results and simulation results, including undamaged model and Rousselier damage model, are shown in Figure 10a. The undamaged model was used to describe the clinching process, and a larger difference between the simulation results and the experiment results is noted. When the Rousselier damage model was used to describe the clinching process, a good consistency was shown in Figure 10a. Therefore, this model is suitable for describing the clinching process and explaining the performance of the AA6082-T6 sheet with 2 mm thickness. In Figure 10b, the shape from simulation is known to be in accordance with the shape from the experiment. In summary, this model can describe the clinching process, including the force-displacement curve, the crack location, and the clinched joint shape. Thus, the orthogonal analysis method is also shown to obtain the suitable damage parameters, and the damage parameters can be obtained by the clinching process.

Concluding Remarks
The present work describes a methodology to identify the Rousselier damage model parameters with the clinched joint experimental-numerical approach. A clinched joint test program includes experiments for both tension-and shear-dominated loading. Based on the paper by Fan [8], the neck of the clinched joint can check the shear damage parameter, whereas the groove of the clinched joint can simultaneously validate the void volume fraction. The main conclusions of this work as follows: (1) The constructs demonstrate a segment-hardening behavior, and this hardening behavior can be described by the post-necking hardening phenomena. When the equivalent strain eq ε is  Figure 10 shows that the displacement-force curve and the shape of the clinched joint is shown in Figure 10. The comparison experiment results and simulation results, including undamaged model and Rousselier damage model, are shown in Figure 10a. The undamaged model was used to describe the clinching process, and a larger difference between the simulation results and the experiment results is noted. When the Rousselier damage model was used to describe the clinching process, a good consistency was shown in Figure 10a. Therefore, this model is suitable for describing the clinching process and explaining the performance of the AA6082-T6 sheet with 2 mm thickness. In Figure 10b, the shape from simulation is known to be in accordance with the shape from the experiment. In summary, this model can describe the clinching process, including the force-displacement curve, the crack location, and the clinched joint shape. Thus, the orthogonal analysis method is also shown to obtain the suitable damage parameters, and the damage parameters can be obtained by the clinching process.  Figure 10 shows that the displacement-force curve and the shape of the clinched joint is shown in Figure 10. The comparison experiment results and simulation results, including undamaged model and Rousselier damage model, are shown in Figure 10a. The undamaged model was used to describe the clinching process, and a larger difference between the simulation results and the experiment results is noted. When the Rousselier damage model was used to describe the clinching process, a good consistency was shown in Figure 10a. Therefore, this model is suitable for describing the clinching process and explaining the performance of the AA6082-T6 sheet with 2 mm thickness. In Figure 10b, the shape from simulation is known to be in accordance with the shape from the experiment. In summary, this model can describe the clinching process, including the force-displacement curve, the crack location, and the clinched joint shape. Thus, the orthogonal analysis method is also shown to obtain the suitable damage parameters, and the damage parameters can be obtained by the clinching process.

Concluding Remarks
The present work describes a methodology to identify the Rousselier damage model parameters with the clinched joint experimental-numerical approach. A clinched joint test program includes experiments for both tension-and shear-dominated loading. Based on the paper by Fan [8], the neck of the clinched joint can check the shear damage parameter, whereas the groove of the clinched joint can simultaneously validate the void volume fraction. The main conclusions of this work as follows: (1) The constructs demonstrate a segment-hardening behavior, and this hardening behavior can

Concluding Remarks
The present work describes a methodology to identify the Rousselier damage model parameters with the clinched joint experimental-numerical approach. A clinched joint test program includes experiments for both tension-and shear-dominated loading. Based on the paper by Xu et al. [8], the neck of the clinched joint can check the shear damage parameter, whereas the groove of the clinched joint can simultaneously validate the void volume fraction. The main conclusions of this work as follows: (1) The constructs demonstrate a segment-hardening behavior, and this hardening behavior can be described by the post-necking hardening phenomena. When the equivalent strain ε eq is less than 0.1522, the Voce function is used as the hardening equation. When the equivalent strain ε eq is more than 0.1522, the equation fits the hardening behavior. Finally, the hardening equation at room temperature for AA6082-T6 sheet is determined. (2) To calibrate the shear damage parameter and the initial void volume fraction simultaneously, a clinched joint is used. The reasons for the use of this sample are analyzed. Orthogonal analysis is used to calibrate the parameters during the clinching process and to build five of the six-level factors. According to the fracture that happens in the groove of the clinching joint, the void volume fraction f is compared with the critical void volume fracture f c to obtain the purpose function. The factors f 1 , f 2 and F are used to evaluate the results. Simultaneously, orthogonal analysis can be described by the sensitivity of the damage parameters on the different deformation regions. The initial void volume fraction f 0 affects the void volume fraction in the neck and the groove of the clinched joint, and the shear damage parameter k ω affects just the void volume fraction f in the neck of the clinched joint. The second optimization obtained the better damage parameters. (3) The good agreement of the numerical predictions with all experimentally-measured force-displacement curves validates the extended modified Rousselier model. Simultaneously, the experimental results on the force-displacement curve, the shape of the clinched joint, and the fracture location are compared with the simulation results. It proved that these damage parameters can describe the deformation process, fracture location, and shape of the clinched joint. These parameters can describe the performance of AA6082-T6 sheet with 2 mm thickness.