Next Article in Journal
Porous Mullite Ceramic Modification with Nano-WO3
Previous Article in Journal
Ultrasensitive Detection of Malachite Green Isothiocyanate Using Nanoporous Gold as SERS Substrate
Previous Article in Special Issue
Multiscale Simulation of Composite Structures: Damage Assessment, Mechanical Analysis and Prediction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Anisotropic Coefficients in the Non-Principal Axis Directions of Tubular Materials Using Hole Bulging Test

1
School of Mechanical Engineering, Dalian University of Technology, Dalian 116024, China
2
State Key Laboratory of High-Performance Precision Manufacturing, Dalian 116024, China
*
Author to whom correspondence should be addressed.
Materials 2023, 16(13), 4629; https://doi.org/10.3390/ma16134629
Submission received: 12 May 2023 / Revised: 23 June 2023 / Accepted: 25 June 2023 / Published: 27 June 2023

Abstract

:
We propose an experimental method to identify anisotropic coefficients in non-principal axis directions of thin-walled tubes. The method involves extracting specimens from the parent tubes and machining a hole in the axial center. The specimens are then inserted into a tube without a hole. The inner diameter of the specimen is theoretically equal to the outer diameter of the inner tube. The double-layer tube undergoes free bulging under internal pressure in our self-developed experimental equipment, with the hole on the specimen expanding simultaneously. The stress states around the hole are uniaxial, and the hole deformation can reflect the anisotropic plastic flow characteristics of the tube. Furthermore, based on the information obtained from the proposed experimental method, a hybrid numerical–experimental method was used to identify the anisotropic coefficients of tubes. Through FE simulations, the relationships between the thickness, stress, and strain states around the hole, the hole shape, and anisotropic coefficients of non-principal axis directions are revealed, and the factors that affect the hole deformation are analyzed. Finally, the hole bulging experiments and FE simulations of AA6061-O extruded tube were conducted, and modeled with Hill48 and calibrated by uniaxial tensile and hoop tensile tests. Its in-plane anisotropy coefficients in any direction are given for the first time which first increase and then decrease from 0° to 90°, reaching a maximum of 1.13 in 60° and a minimum of 0.69 in 0°. This work can provide the key experimental data for establishing an accurate anisotropic plastic constitutive model of thin-walled tubes.

1. Introduction

Tube hydroforming is an advanced forming technology for the production of hollow and variable cross-section complex parts which are widely used in the aerospace and auto industries [1,2,3,4]. However, tube hydroforming results are greatly affected by the loading path, and extremely prone to cracking or wrinkling [5,6,7]. The shapes of the formed parts are generally complex, and thin-walled tubes exhibit significant plastic anisotropy [8,9,10], making it difficult to accurately predict the reasonable loading paths through theoretical analysis and experiments [11]. Therefore, FE simulation is usually used to determine and optimize process parameters before experiments [12]. However, the accuracy of FE simulation largely depends on the accuracy of the constitutive model used to describe the material properties including the plastic flow, yielding, and hardening characteristics [13]. Hence, constructing an accurate plastic constitutive model for thin-walled tubes is crucial for achieving accurate FE simulation of tube hydroforming. The anisotropy coefficients are usually used to calibrate the plastic potential and yield function, which are important for accurately reflecting the plastic deformation. Therefore, the anisotropy coefficients of the thin-walled tube in any direction are very important for constructing its accurate plastic constitutive model [14].
However, it is impossible to measure the anisotropy coefficients in any direction through uniaxial tensile tests such as sheet metal for thin-walled tubes due to its closed structure along the hoop direction [15]. At present, only the anisotropy coefficients along the axial and hoop directions can be measured by the uniaxial tensile test and the hoop tensile test [16,17,18], respectively. Some researchers performed uniaxial tensile tests by flattening specimens which are cut in different directions of the tube [19]. However, this method would produce pre-strain or even rupture during the flattening process, which causes the test results to be inaccurate.
In order to obtain the anisotropy coefficients of thin-walled tubes in the non-principal axis directions, some scholars first calibrate the anisotropic plastic constitutive model of thin-walled tubes through some experimental data of tension-tension [20,21,22,23] or even tension-compression [24,25,26] stress states, and then reverse determine the anisotropy coefficient through the flow rule. However, the two kinds of experimental data can only be used to calibrate the coefficients related to normal stress in the plastic constitutive model, while the coefficients related to shear stress in the plastic constitutive model cannot be determined. Due to lack of experimental data to calibrate the coefficients related to shear stress components in the plastic constitutive model of tubes, some researchers [26,27] assume them to be isotropic. Then the anisotropic coefficients of tubular materials in any direction are inversely calculated by the flow rule. For example, Kuwabara et al. [14] assumed that the undetermined coefficients α 7 and α 8 in the Yld2000 yield criteria are 1. However, the yield criteria based on this assumption cannot correctly describe the anisotropic plastic flow characteristics in the non-principal axial directions, so the anisotropic coefficients of thin-walled tubes in any non-principal axis directions determined by this method are not accurate.
To address the problem, Zhang et al. [25] proposed an experimental method for measuring the axial pure shear yield stress of thin-walled tubes. Based on the assumption that the coefficients α 7 and α 8 of the Yld2000 yield criteria are equal, α 7 and α 8 are determined by pure shear yield stress. Further, the anisotropy coefficients in any direction of thin-walled tubes can be determined inversely by the associated flow rule. However, the axial pure shear yield stress can only reflect the anisotropic yield characteristics of thin-walled tubes, and cannot accurately reflect the anisotropic plastic flow characteristics in the non-principal axis directions. The anisotropy coefficients of thin-walled tubes in the non-principal axis directions determined by this method are still inaccurate.
In summary, the anisotropy coefficients of thin-walled tubes in the non-principal axis directions are the key data for calibrating the plastic constitutive models. However, so far, there are few studies about it that have been publicly reported.
In this paper, in order to identify anisotropic coefficients in non-principal axis directions of thin-walled tubes, a new experimental method of the hole bulging test (HBT) is proposed, which is inspired by the hole expansion tests of sheets. The deformation characteristics at the hole periphery can reflect the anisotropic plastic flow characteristics of tubes. Based on this experiment, the anisotropy coefficients of thin-walled tubes in any non-principal axis direction are determined by a hybrid numerical–experimental method. The FE simulations and experiments are used to analyze the deformation characteristics and influencing factors of the hole periphery during the experiment, and the feasibility of the method is verified. Furthermore, the hole bulging experiments of aluminum alloy AA6061-O extruded tube are conducted, and its in-plane anisotropy coefficients in the non-principal axis directions are successfully obtained, which will provide the key experimental data for establishing an accurate anisotropic plastic constitutive model.

2. Hole Bulging Test (HBT) and Hybrid Numerical–Experimental Method

2.1. Test Procedure of Hole Bulging Test

The hole bulging test (HBT) is proposed to determine the anisotropy coefficients of tubular materials, the idea of which is inspired by the hole expansion test (HET) of sheet metal. HET is a widely used method to measure the forming properties of the edge quality [28,29,30,31]. In recent years, more than the formability evaluation, the HET is increasingly being used to validate the anisotropic plastic constitutive model of materials, including plastic anisotropy [30,31,32] due to the following reasons: (1) larger deformation can be achieved than a standard uniaxial tensile test due to compatibility with the surrounding materials; (2) the material is deformed under various stress states between uniaxial tension, plane strain tension, and equal biaxial tension from the hole edge and radially inland [29]. Furthermore, the deformation of the hole during the HET has been used to determine the anisotropy coefficients of sheet metal in any direction [33].
For the HBT, the specimen is a tube with a hole in the axial center. Due to the fact that the tube is a circumferentially closed structure, how to deform the hole is completely different from how to expand the hole in the sheet metal. To solve this problem, the tubular specimen is inserted into a tube without a hole (inner tube), and the inner diameter of the tubular specimen is theoretically equal to the outer diameter of the inner tube. The double-layer tube undergoes free bulging under the internal pressure applied to the inner tube, and the hole on the specimen expands simultaneously. The schematic of the HBT is shown in Figure 1. The edge of the hole on the tubular specimen is free, so that it can deform freely and the stress state is uniaxial, which is similar to that of the hole expansion test of sheet metal but completely different from the stress state of the single-layer tube free bulging test (see Figure 2). The stress direction of any point at the hole periphery is along the tangent direction. Each tangent line can represent a direction in the tube plane (see Figure 2b). All directions in the tube plane are considered to be involved in the deformation of the hole periphery. From the stress–strain analysis of Section 4.2, the plastic deformation occurs in almost all directions of the hole periphery. Therefore, the deformation characteristics of the hole (such as the hole shape, the thickness of the hole periphery, etc.) can reflect the anisotropic plastic flow characteristics in any direction in the tube plane, which can be used to determine the anisotropy coefficients in the tube plane.

2.2. Identification of Anisotropy Coefficients by Hybrid Numerical–Experimental Method

The HET can directly obtain the anisotropy coefficients of the sheet metal in any direction by measuring the strain [33]. However, due to the particularity of the tube structure, the anisotropy coefficients of the tube cannot be directly obtained by the same method. Therefore, the hybrid numerical–experimental method shown in Figure 3 is proposed to determine the in-plane anisotropy coefficients of tubes.
The method is based on the results of the deformation characteristics of the hole in the HBT, and the simulation analysis is carried out by changing the undetermined parameters of the constitutive model. When the error between the results of FE simulation and experiments is reduced to an allowable range, the parameters of the constitutive model used in the simulation are the actual values of the material. The parameters of the constitutive model are generally calibrated by the in-plane anisotropy coefficient r θ , the subscript θ represents the angle with respect to the axial direction. The anisotropy coefficients in the principal axis directions, r0 and r90, can be measured by experiments, while the anisotropy coefficients in other directions such as r45 cannot be obtained by experiments which are the undetermined parameters of the constitutive model and need to be determined through the hybrid numerical–experimental method, as shown in Figure 3. After the anisotropy coefficients of these specific non-principal directions are determined, the constitutive model can be determined, and further, the anisotropy coefficients of any other non-principal directions can be determined by the constitutive model.
The deformation characteristics of the hole in the HBT include the hole shape and the thickness at the hole periphery. The hole shape can be characterized by the long axis of the deformed hole in this study (hereinafter referred to as the aperture).
The accurate determination of the anisotropy coefficients needs a reasonable cost function. The cost function is defined as the root mean square error (RMSE) between the results of FE simulations and experiments:
η ¯ = 1 N i = 1 N ( R i e x p R i f e m R i e x p ) 2
where N is the number of experiments, R i e x p and R i f e m are the values of the iterative metric in the results of experiment and simulation, respectively.
The number of non-principal anisotropy coefficients to be determined is related to the constitutive model selected. When there is only one coefficient to be determined, a reasonable RMSE threshold Δ η is selected. Firstly, the bisection iteration is carried out to quickly reduce the iteration interval. When η ¯ Δ η terminates the iteration, the r θ is the optimal result. If the iteration step size Δ r θ < 0.05 and the RMSE η ¯ > Δ η appears in the iteration process of the bisection iteration, the fixed step search is carried out in the interval. The search step size of the fixed step search is Δ r θ = 0.01 . If the iteration error increases as the iterative step size increases or decreases, it indicates that the iteration error has reached the minimum value. Therefore, the r θ at this point can be considered the optimal result. When the number of non-principal anisotropy coefficients to be determined is more than one, the grid search method can be used to determine the optimal solution of each non-principal anisotropy coefficient.

2.3. Experiment Equipment of Hole Bulging Test

The home-developed THF.HIT-160/110-A tube mechanical properties test equipment was used for the HBT, as shown in Figure 4. The equipment is composed of mechanical structure, high-pressure system, hydraulic system, control system, and data analysis system. The free bulging test of the tube can be conducted on it, through which the mechanical properties of tubes such as stress–strain component curves and hardening curves can be obtained. A maximum axial force of 50 t and a maximum working pressure of 160 MPa can be provided by the equipment.

2.4. Advantages and Disadvantages of the HBT

The advantages of the HBT include the following: (1) The specimen is a tube with a hole, which is directly extracted from the parent tube without changing the curvature shape of the tested tube and introducing additional strain hardening; (2) The edge of the hole is a free boundary and is subjected to uniaxial stress state, stress state is simple; (3) The experimental process only requires the application of internal pressure, which is simple and does not require a complex feedback control system; (4) The experiment can be conducted on free bulging test equipment without the need for additional development of new devices.
Although this experiment has the above advantages, it also has some disadvantages, such as: (1) the influencing factors during the experimental process are unclear, (2) the feature quantities of hole edges during deformation are difficult to be measured, and so on. Therefore, extensive research is needed.

3. Material, and Anisotropy Coefficients in the Principal Axis Directions

3.1. Material

The material of this study is aluminum alloy AA6061-O extruded thin-walled tubes, which have obvious in-plane anisotropy and are widely used in aerospace, automotive, and other fields due to their excellent mechanical properties. The initial outer diameter of the tube is 60 mm and the initial thickness is 2 mm.

3.2. Anisotropic Yield Criterion

The Hill48 yield criterion is used to characterize the plastic anisotropy of the material, due to this study only focusing on the in-plane anisotropic plastic flow characteristics of the tube. For the plane stress state, the Hill48 yield criterion is expressed as [34]:
2 f ( σ i j ) ( G + H ) σ 11 2 2 H σ 11 σ 22 + ( H + F ) σ 22 2 + 2 N σ 12 2 = 1
where F, G, H, and N are the undetermined coefficients of Hill48 yield criterion.
Substituting the equivalent stress σ ¯ , which is the flow stress of the axial tensile test, into Equation (2), the following equation can be obtained as:
G + H = 1 / σ ¯ 2
In order to accurately describe the in-plane anisotropic plastic flow characteristics of tubes, the in-plane anisotropy coefficients, including the axial anisotropy coefficient r0, diagonal anisotropy coefficient r45 and hoop anisotropy coefficient r90, are used to determine the undetermined coefficients of the Hill48 yield criterion, as shown in the following equations,
r 0 = H G ;         r 90 = H F ;         r 45 = N F + G 1 2
According to Equations (3) and (4), the coefficients can be expressed as:
{ G + H = 1 / σ ¯ 2 2 H = 2 r 0 σ ¯ 2 ( 1 + r 0 ) H + F = r 0 ( 1 + r 90 ) σ ¯ 2 r 90 ( 1 + r 0 ) 2 N = r 0 + r 90 σ ¯ 2 r 90 ( 1 + r 0 ) ( 2 r 45 + 1 )
Thus, the Hill48 yield criterion can be expressed as:
σ 11 2 2 r 0 1 + r 0 σ 11 σ 22 + r 0 ( 1 + r 90 ) r 90 ( 1 + r 0 ) σ 22 2 + r 0 + r 90 r 90 ( 1 + r 0 ) ( 2 r 45 + 1 ) σ 12 2 = σ ¯ 2
The anisotropic coefficients r0, r90 and hardening curve σ ¯ ( ε ¯ ) can be accurately determined by the uniaxial tensile tests and the hoop tensile test. While the anisotropic coefficient r45 will be determined by the hybrid numerical–experimental method as shown in Section 2.2.

3.3. Measurement of Anisotropy Coefficients in the Principal Axis Directions

The uniaxial tensile tests are conducted on the LE5105 100 KN electronic universal tensile testing machine to measure the axial anisotropy coefficient r0 of AA6061-O tube. Axial uniaxial tensile specimens are designed according to the Chinese Standard of GB/T 228.1-2010 [35] and cut along the axial orientation of the tube as shown in Figure 5 (every 90° in the hoop direction).
The strain is measured in real-time by a 3D digital image correlation (3D-DIC). The stress–strain curves measured using specimens cut from four locations are almost overlapping, so only one result is plotted in Figure 6a. As shown in Figure 6b, the average value of r0 is 0.69, which does not change significantly with the increase in strain.
The hoop anisotropy coefficient r90 of AA6061-O tubes is measured by the hoop tensile test, and the corresponding specimen and experimental device are shown in Figure 7. The length–width ratio of the gauge section of the specimen is 4:1 and other dimensions are designed according to the Chinese standard GB/T 228.1-2010, as shown in Figure 7a.
As shown in Figure 7b, the hoop tensile test device consists of two clevis grips and two motion blocks. The specimen is elongated and deformed by the separation of the upper and lower motion blocks. As shown in Figure 7c, compared with the traditional device with a central mandrel, the central mandrel and the lower moving block are designed as a single unit, which can reduce the influence of fixture jitter during the testing process. Keeping the gauge section at the 3 o’clock and 9 o’clock orientations of the fixture and lubricating the specimen and fixture with PVC films can minimize the influence of friction between the fixture and the specimen. The test is performed on the LE5105 100 KN electronic universal tensile testing machine with DIC for strain measurement in the gauge section. The specimens after the hoop tensile test are shown in Figure 8, and the measured r90 of AA6061-O aluminum alloy tubes is 1.05.

4. FE Simulation

4.1. FE Model

In this section, ABAQUS is used to simulate the HBT. As shown in Figure 9, the FE model includes the tubular specimen with a hole, the inner tube, the punches, and the clamping sleeves. The simulation process is consistent with the experiment, which includes two analysis steps: first, the displacements of the left and right punches are controlled to achieve the seal of the double-layer tube, and then the internal pressure of the double-layer tube is applied according to the actual loading path of the experiment. The punches and the clamping sleeves are set as rigid bodies. The tubular specimen with a hole and the inner tube are set to be deformable. The elastic modulus, hardening curve, and other parameters are obtained by the uniaxial tensile tests. The contact between each surface is set as face-to-face contact, and Coulomb’s law is chosen. The friction coefficient between the inner tube and the specimen is 0.02, between the inner tube and the punch is 0.15, and other contact regions are 0.17.
As shown in Figure 10, the tubular specimen with a hole adopts a 4-node linear reduced integration shell element (S4R). Five integral points are set in the thickness direction, and the global element size of the tube is set to 2.0 mm. Local mesh refinement is performed around the hole, and 80 elements are evenly distributed around the hole to ensure simulation accuracy. The global size of the inner tube, the punches, and the sleeves grid are set to 3.0 mm.

4.2. Stress–Strain Analysis of the Hole Periphery

The stress–strain characteristics of the hole periphery in deformation are analyzed by FE simulation. When analyzing the stress and strain of each point on the edge of the hole, the coordinate system is established in the radial (r direction), tangent (ta direction), and thickness (t direction) of the point, which changes with the position of the analysis point. Each point of the hole periphery can be represented by the angle between the point and the tube axis direction.
Figure 11 shows the stress distribution at each point of the hole periphery under different bulging heights. The radial stress of each point at the hole periphery is almost 0 MPa at different bulging heights, while tangent stress is completely different from radial stress. From 0° to 90°, the tangential stress gradually transforms from tensile stress to compressive stress, with the maximum tensile stress at 0° and the maximum compressive stress at 90°. Therefore, it can be approximated that each point around the hole is only subjected to tangential stress and is in a uniaxial stress state. The direction of stress at each point around the hole can represent one specific orientation in the tube plane. In addition, the stress in almost all regions of the hole periphery exceeds the yield strength of the material, which indicates that plastic deformation occurs in almost all orientations around the hole, that is, all directions in the tube plane are considered to involve the deformation around the hole. Therefore, the deformation characteristics around the hole can reflect the differences in plastic flow in any orientation within the tube plane, and can be used to further determine the anisotropy coefficient of the tube.
Figure 12 shows the radial, tangent, and thickness strain contour diagrams around the hole. The tangent strain at the hole periphery is the largest, with the maximum tensile strain in the 0° and 180° orientations and the maximum compressive strain in the 90° and 270° orientations, which is consistent with the variation law of tangential stress shown in Figure 11. Due to the fact that the edge of the hole is in a tangential uniaxial stress state, the position with a large positive tangent strain corresponds to a large thickness reduction, the rupture around the hole is most likely to occur at 0° and 180° positions, and the smaller the r value, the earlier the rupture occurs. The position with the negative tangential strain exhibits thickening, and the smaller the r value, the more thickening occurs. Therefore, it can be concluded that the thickness variation around the hole is closely related to the in-plane anisotropy of the tubular specimen. The thickness characteristics around the hole can reflect the differences of plastic flow in any direction within the tube plane and can be used to further determine the anisotropy coefficient of the tube.
Figure 13 shows the strain distribution at the hole periphery with different bulging heights. It can be seen that the radial strain and the thickness strain are different in the same position, and the difference is especially significant at the positions with large deformation, such as 0° position and 180° position. As the bulging height increases, the strain difference between the radial and thickness directions of each point at the hole periphery is more obvious. Because the stress state around the hole is uniaxial, the difference between the radial and thickness strains can directly reflect the in-plane anisotropy of the tube. Therefore, the difference between the radial strain and thickness strain around the hole can be used to determine the in-plane anisotropy of the tube. When the bulging height is constant, the greater the difference between the two, indicating that the tube has more obvious in-plane anisotropy.
Through the stress–strain analysis of the hole periphery, the deformation characteristics around the hole, such as the distribution of thickness and the difference between the radial and thickness strains, can reflect the anisotropic plastic flow characteristics of tubes, which can be used to determine the in-plane anisotropy coefficients of tubes.

4.3. Deformation Characteristics of the Hole Periphery

Hole deformation characteristics mainly include aperture, thickness of hole periphery, etc. In this section, the relationship between the deformation characteristics around the hole and in-plane anisotropy coefficients, is analyzed by FE simulation. The analysis of the deformation characteristics at the hole periphery, which are highly sensitive to the anisotropy coefficients, is important for determining the anisotropy coefficients in the non-principal axis direction of the tube. Taking the Hill48 yield criterion as an example, the undetermined anisotropy coefficient is r45. The influence of r45 on the deformation characteristics of the hole periphery will be analyzed. The r45 is taken every 0.1 in the range of 0.30–1.70 for the simulation analysis. Figure 14 shows the relationship between the apertures and the bulging height with different r45. It can be seen that different r45 have obvious regularity on the influence of the apertures after bulging. The larger the r45, the larger the aperture at the same bulging height, and the aperture difference will increase with the increase in bulging height. The apertures for different r45 were extracted when the bulging height was 2.80 mm, as shown in Figure 15. At the same bulging height, the larger the r45, the larger the aperture after bulging. There is a one-to-one correspondence between r45 and the apertures, which indicates that it is feasible to determine r45 inversely by the apertures.
At 2.5 mm bulging height, the thickness around the hole was extracted, as shown in Figure 16. The thickness reduction is the largest at the orientations of 0° and 180°, reaching about 0.6 mm, while there is almost no thinning in the direction of 90° and 270° (note that initial thickness is 2 mm), and even a little thickening, which is consistent with the results of Section 4.2 hole periphery strain analysis.
As shown in Figure 17, the variations of thickness with bulging height in three orientations of 0°, 45°, and 90° under different r45 are analyzed. The thicknesses of 0°and 90° are more sensitive to the changes of r45 than that of 45°. The thickness of 0° is the most sensitive to r45. When r45 changes from 0.3 to 1.7 with the bulging height of 2.5 mm, the difference in thickness thinning at the position of 0° is 0.27 mm, and the difference rate reaches 13.5%. This is because throughout the entire HBT, the stresses near 45° and 90° positions are always very small which lead to a small deformation amount at both locations, so the differences in deformation that could be reflected are very small. However, the stress near 0° is always the maximum stress around the hole, and the deformation amount is large, which naturally better reflected the difference in deformation, as shown in Figure 11. Fortunately, the wall thickness at the 0° position can well reflect the changes in r45, so the wall thickness, especially at the 0° position, can be used to determine the anisotropy coefficient of the thin-walled tubular specimens in reverse.
In summary, both the aperture of the hole and the thickness around the hole (especially thickness in the 0° orientation) can reflect the in-plane anisotropy of tubes, and have a significant correlation with the undetermined anisotropy coefficient r45. Considering that the aperture is easy to measure accurately, the aperture is chosen as the iterative metric of the hole deformation used to determine the anisotropy coefficients in the non-principal axis directions in Section 2.2.

4.4. Influencing Factors of Hole Deformation

The strength relationship between the inner tube and the tubular specimen, the initial aperture, and the friction coefficient between the inner and outer tube on the hole deformation were analyzed by FE simulations. The optimal experimental conditions are determined to make the hole deformation characteristics significantly reflect the in-plane anisotropy of tubes and reduce the influence of external factors on the experimental results.

4.4.1. Strength Relationship between the Inner and Outer Tube

The inner and outer tube strength-matching relationship can be categorized as “high internal strength” (HIS), “equal strength” (ES), and “low internal strength” (LIS). As shown in Table 1, the ultimate bulging height is the maximum expansion height of the hole before rupturing at r45 = 1.5, the aperture is the major axis length of the hole under the ultimate bulging height in this study, and the deformation difference is the difference between the long axis after experiment and the initial aperture divided by the initial aperture. From Table 1, it can be seen that the ultimate bulging height, aperture, and deformation difference degree corresponding to HIS are the largest, which will be more beneficial to reflect the deformation difference caused by different r45 values. As shown in Figure 18, the overall hole size of HIS is the largest at the ultimate bulging height, and the difference in hole size between different anisotropic coefficients is the most significant. Thus, it is recommended to choose HIS for the HB experiments, and a 304 stainless steel tube was selected as the inner tube in this paper as it has higher strength than AA6061-O tube. In addition, 304 stainless steel tubes are generally isotropic, so it is modeled with the Mises yield criterion in the simulation.

4.4.2. Effect of Initial Aperture

The hole bulging tests with initial apertures of 2 mm, 3 mm, 4 mm, and 5 mm were simulated and analyzed, respectively, and the results are shown in Table 2. The ultimate bulging height in the table is the maximum bulging height before the hole ruptures when r45 = 1.5. The aperture and deformation difference degree are consistent with the definition in Section 4.4.1. A small initial aperture corresponds to a large ultimate impacting height and a large deformation difference degree.
Figure 19 shows the thickness distribution at the hole periphery with different initial apertures at the ultimate bulging height. The results show that the thickness distribution is not very different when the initial aperture changes. The thickness at a position of 0° is the smallest, and the thickness at a position of 90° is the largest.
From the above results, the smaller the initial aperture is, the larger the ultimate bulging height is, and the larger the hole deformation is. However, a small initial aperture will increase the difficulty of measurement, resulting in increased measurement error. Although the ultimate bulging height of the large aperture is smaller, the size of the aperture changes more greatly, which is easier to be measured. To balance this problem, preliminary HB experiments were carried out on 2 mm and 4 mm apertures, as shown in Figure 20. The experimental results are basically consistent with the simulation results. The ultimate bulging height of the tube with 2 mm aperture is 3.34 mm, while at 4 mm aperture it is 3.10 mm. However, when measuring the aperture after deformation, the 2 mm aperture tube is difficult to be measured and the measurement error is large, while more accurate measurement values can be obtained for the 4 mm aperture tube. Therefore, to obtain more accurate experimental results, tubes with 4 mm initial aperture are chosen to conduct HB experiments.

4.4.3. Friction between Inner and Outer Tubes

For different r45, the friction coefficients between the inner and outer tubes are set to 0.1, 0.2, and 0.3 for FE simulations, and the apertures at the bulging height of 2.75 mm are extracted, as shown in Table 3. It can be seen that the difference of apertures at the same bulging height is very small when r45 is the same and the friction coefficient increases from 0.1 to 0.3, with the maximum difference not exceeding 1%. Therefore, it can be concluded that the friction coefficient has little effect on the deformation of the hole and is not an important factor. During the HB experiments in this paper, the inner tube is wrapped by PVC films to fully reduce the friction between the two tubes.

5. Hole Bulging Experiment

5.1. Experimental Process

Before the HB experiments, the AA6061 tubes were annealed by holding the temperature at 405 °C for 3 h, cooled to 260 °C with the furnace, and then air-cooled. The total length of the tubular specimen with a hole and the length of the bulging zone used in this study are 160 mm and 70 mm, respectively. A hole with a diameter of 4 mm was drilled in the middle of the tubular specimen. The inner tube is a 304 stainless steel tube, the initial diameter, thickness, and total length of which are 55 mm, 1 mm, and 160 mm, respectively.
The measurement points are marked at the hole periphery as shown in Figure 21. The surface of the inner tube is wrapped around by PVC films. The inner tube is placed inside the tubular specimen with a hole. Then the double-layer tube is installed on the equipment shown in Figure 4. The double-layer tube is sealed and loaded according to the preset pressure loading path to bulge the tubular specimen with holes. After the experiments, the apertures, the thicknesses at the hole periphery, and the bulging zone contours were measured.

5.2. Experimental Results

Figure 22 is the tubular specimen with an initial hole diameter 4 mm after HB experiment with 3.31 mm bulging height. The hole ruptures at the positions of 0° and 180°, and local thickening occurs at the positions of 90° and 270°, which are consistent with the simulation results.
Seven groups of experiments under different bulging heights are conducted on the same batch of tubes. During the HB experiments, the bulging height is measured in real-time by a laser displacement sensor. The results are shown in Table 4 and Figure 23.

6. Results and Verification

6.1. Results of r-Value

When determining anisotropy coefficients in the non-principal axis directions using the method in Section 2.2, it is necessary to determine the iterative range of the undetermined anisotropy coefficients. In this study, the iterative range of the anisotropy coefficient r45 is determined by comparing the results of HB experiments and FE simulations. As shown in Figure 24, the bulging height–aperture curve obtained by the HB experiments is located between the FE simulation results of r45 = 0.9 and r45 = 1.1. Therefore, it can be concluded that the actual r45 of the tube is between 0.9 and 1.1. The iterative range of the r45 is determined to be [0.9, 1.1].
Through the identification process of anisotropic coefficients proposed in Section 2.2, the RMSE (see Equation (1)) threshold is Δ η = 0.5 % , and the iteration results are shown in Table 5. Finally, r45 = 1.06 is determined, and the RMSE between the results of experiments, and FE simulations of r45 = 1.06 is 0.47%. The difference in the bulging height–aperture curve is also the smallest, as shown in Figure 25, and the maximum error is 0.51%, which indicates that the proposed method is reasonable.

6.2. Verification

Other characteristics of the tubular specimen with a hole after the experiments can be used to verify the accuracy of the determined anisotropy coefficients, such as the thickness around the hole and the profile of the bulging zone. The r45 determined in Section 6.1 is substituted into the simulation model for FE simulation, and the simulation results of the hole periphery thickness and the contour of the bulging zone are compared with the experimental results.
Figure 26 shows the results of the thickness around the hole at the bulging height of 2.50 mm, obtained by the experiment, and FE simulation of r45 = 1.06. The thicknesses obtained by FE simulations are in good consistency with the experimental results, and the average error is only 1.25%. From the experimental results, the thickness reduction is the largest at the position of 0°, which is 0.65 mm. The thickness reduction decreases gradually from 0° to 90°, and there is a very small thickness thickening at the position of 90°.
As shown in Figure 27a, the whole bulging zone is divided into 15 sampling points. The profiles of the bulging zone obtained by the HB experiments at the bulging heights of 2.10 mm, 2.41 mm, and 2.74 mm are compared with the corresponding results of FE simulation, as shown in Figure 27b. The comparison results indicate that they are in good consistency, and the average error of the contour is the largest at the bulging height of 2.74 mm, which is 1.94%.
From the comparative results of the thickness around the hole and the profile of the bulging zone, the hole deformation characteristics of the HBT can reflect the anisotropic plastic flow characteristics of tubes. It is feasible to determine the non-principal axis directions anisotropy coefficient of tubes through the hybrid numerical–experimental method based on the HBT, and the anisotropy coefficients obtained can guarantee high accuracy.

6.3. Results of AA6061-O in-Plane Anisotropy

The Hill48 plastic potential equation of AA6061-O extruded aluminum alloy tubes can be obtained by substituting r0 = 0.67, r45 = 1.06, and r90 = 1.05 determined in this study into Equation (6).
σ 11 2 0.80 × σ 11 σ 22 + 0.78 × σ 22 2 + 3.09 × σ 12 2 = σ ¯ 2
According to the Drucker flow rule, the anisotropy coefficients of AA6061-O extruded aluminum alloy tubes in any direction θ is:
r θ = 1.45 × cos 4 θ + 0.95 × sin 4 θ + cos 2 2 θ + 1.87 × sin 2 2 θ 1.45 × cos 2 θ + 0.95 × sin 2 θ 1
The anisotropy coefficients of AA6061-O extruded aluminum alloy tubes in any non-principal axis directions are given for the first time, as shown in Figure 28. The in-plane anisotropy coefficients of AA6061-O extruded aluminum alloy tubes increase first and then decrease from 0° to 90°, and reach the maximum value of 1.13 in the orientation of 60°.

7. Conclusions

In this study, in order to determine the anisotropic coefficients in non-principal axis directions of thin-walled tubes, a new experiment was proposed, which was called the hole bulging test (abbreviated as HBT). The specimen of HBT is a tube with a hole in the axial center. The specimen is inserted into a tube without a hole. The hole on the specimen expands simultaneously as the double-layer tube undergoes free bulging under internal pressure. Furthermore, based on the information obtained from the proposed experimental method, a hybrid numerical–experimental method was used to identify the anisotropic coefficients of tubes. The thickness, stress, and strain state around the hole, the hole shape, as well as the factors that affect the hole deformation were analyzed through FE simulations and the aperture reflecting the hole shape was determined as the iterative metric of the hybrid numerical–experimental method. Uniaxial tensile tests and hoop tensile tests were performed to determine the tubular materials’ yield stresses, hardening curve, and anisotropy coefficients with respect to axial and hoop directions. The HB experiments and FE simulations of the aluminum alloy AA6061-O extruded tube were performed, and the undetermined anisotropic coefficient in the plastic constitutive model was identified by the hybrid numerical–experimental method and the accuracy of it was further verified. Finally, the in-plane anisotropic coefficients in any direction of AA6061-O extruded tube were determined through the plastic flow rule. The main conclusions are as follows:
  • In the HBT, the stress state around the hole is uniaxial, and the hole deformation is a comprehensive result of the deformations of various points around the hole, which can reflect the in-plane anisotropic plastic flow characteristics of the tube.
  • The aperture of the hole and thickness around the hole after deformation vary significantly with the variation of anisotropy coefficients in the non-principal axis directions of thin-walled tubes, and can be used to determine them.
  • For the HBT, a higher strength of the inner tube than the tubular specimen, and a smaller initial diameter of the circular hole are recommended. The friction coefficient between the double-layer tube has little effect on the hole deformation.
  • Compared with the experimental results, both the average errors of the thickness around the hole and the profile of the bulging zone obtained from the final iterative simulation analysis do not exceed 2%, verifying the feasibility of the proposed method in this paper.
  • The aluminum alloy AA6061-O extruded tube exhibits significant in-plane anisotropic plastic flow characteristics, and its in-plane anisotropic coefficients in any direction are given for the first time, which increase first and then decrease from 0° to 90°, reaching a maximum value of 1.13 at 60°, a minimum value of 0.69 at 0°.
The above conclusions reveal that the HBT is a promising test for accurate identification of the in-plane anisotropic coefficients of the non-principal axial directions of thin-walled tubes. This study only investigates the feasibility and application effect of the HBT. In future research, the feasibility of directly measuring the in-plane anisotropic coefficients of the non-principal axial directions of thin-walled tubes through the HBT will be explored.

Author Contributions

The research in this paper was performed as a collaboration between all the authors. Y.L. and Z.H. designed the project and research theme and provided scientific guides and discussions. Y.W., Y.S. and J.L. performed the uniaxial tensile test, hoop tensile tests and the hole bulging test. Y.L. and Z.H. solved the model building process. Y.L. and Y.W. wrote the manuscript. Z.H. and K.C. improved the manuscript. All authors analyzed the correction results. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (52075075) and supported by Liaoning Revitalization Talents Program (No. XLYC1802065) and supported by the Fundamental Research Funds for the Central Universities of China (DUT20RC(5)031 and DUT23RC(3)026).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

HBTHole bulging test, a new experiment proposed in this paper HET hole expansion test of sheet metal
r θ Anisotropy coefficient of the tube, θ is the direction relative to the axial direction R i e x p the values of the iterative metric in the results of experiment
R i f e m the values of the iterative metric in the results of simulation Δ η the iteration step
σ ¯ the equivalent stress, which is the yield stress of the axial tensile test σ i j Stress tensor
HISHigh internal strength, that is, the strength of the inner tube is higher than the strength of the tubular specimenESEqual strength, that is the strength of the inner tube is equal to the strength of the tubular specimen
LISLow internal strength, that is, the strength of the inner tube is lower than the strength of the tubular specimen

References

  1. Yuan, S.J.; Fan, X.B. Developments and Perspectives on the Precision Forming Processes for Ultra-Large Size Integrated Components. Int. J. Extrem. Manuf. 2019, 1, 22002. [Google Scholar] [CrossRef]
  2. Lang, L.H.; Liu, K.N.; Cai, G.; Yang, X.; Guo, C.; Bu, G. A Critical Review on Special Forming Processes and Associated Research for Lightweight Components Based on Sheet and Tube Materials. Manuf. Rev. 2014, 1, 9. [Google Scholar] [CrossRef]
  3. He, X.; Ding, J.; Habibi, M.; Safarpour, H.; Safarpourf, M. Non-polynomial Framework for Bending Responses of the Multi-scale Hybrid Laminated Nanocomposite Reinforced Circular/Annular Plate. Thin-Walled Struct. 2021, 166, 108019. [Google Scholar] [CrossRef]
  4. Liu, Z.; Su, S.; Xi, D.; Habibi, M. Vibrational Responses of a MHC Viscoelastic Thick Annular Plate in Thermal Environment Using GDQ Method. Mech. Based Des. Struct. Mach. 2020, 1, 108019. [Google Scholar] [CrossRef]
  5. Chen, X.; Yu, Z.; Hou, B.; Li, S.; Lin, Z. A Theoretical and Experimental Study on Forming Limit Diagram for a Seamed Tube Hydroforming. J. Mater. Process. Technol. 2011, 211, 2012–2021. [Google Scholar] [CrossRef]
  6. Hashemi, S.J.; Naeini, H.M.; Liaghat, G.; Tafti, R.A.; Rahmani, F. Forming Limit Diagram of Aluminum Aa6063 Tubes at High Temperatures by Bulge Tests. J. Mech. Sci. Technol. 2014, 28, 4745–4752. [Google Scholar] [CrossRef]
  7. Nikhare, C.; Weiss, M.; Hodgson, P.D. Fea Comparison of High and Low Pressure Tube Hydroforming of Trip Steel. Comput. Mater. Sci. 2009, 47, 146–152. [Google Scholar] [CrossRef]
  8. Christodoulou, N.; Turner, P.A.; Ho, E.T.C.; Chow, C.K.; Resta Levi, M. Anisotropy of Yielding in a Zr-2.5Nb Pressure Tube Material. Metall. Mater. Trans. A-Phys. Metall. Mater. Sci. 2000, 31, 409–420. [Google Scholar] [CrossRef]
  9. Matlock, D.K.; Cho, S.H.; Speer, J.G.; Al Jabr, H.M.; Zhang, P. Anisotropy of Mechanical Properties of Api-X70 Spiral Welded Pipe Steels. Mater. Sci. Forum 2013, 753, 538–541. [Google Scholar] [CrossRef]
  10. Yuan, S.J.; Xie, W.C.; Lin, Y.L.; He, Z.B. Analytical Model and Testing Method for Equivalent Stress-Strain Relation of Anisotropic Thin-Walled Steel Tube. J. Test. Eval. 2019, 47, 775–790. [Google Scholar] [CrossRef]
  11. Banabic, D. Sheet Metal Forming Processes: Constitutive Modelling and Numerical Simulation; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  12. Cai, Y.; Wang, X.S.; Yuan, S.J. Pre-Form Design for Hydro-Forming of Aluminum Alloy Automotive Cross Members. Int. J. Adv. Manuf. Technol. 2016, 86, 463–473. [Google Scholar] [CrossRef]
  13. Jansson, M.; Nilsson, L.; Simonsson, K. On Constitutive Modeling of Aluminum Alloys for Tube Hydroforming Applications. Int. J. Plast. 2005, 21, 1041–1058. [Google Scholar] [CrossRef]
  14. Kuwabara, T.; Yoshida, K.; Narihara, K.; Takahashi, S. Anisotropic Plastic Deformation of Extruded Aluminum Alloy Tube Under Axial Forces and Internal Pressure. Int. J. Plast. 2005, 21, 101–117. [Google Scholar] [CrossRef]
  15. Kuwabara, T. Advances in Experiments on Metal Sheets and Tubes in Support of Constitutive Modeling and Forming Simulations. Int. J. Plast. 2007, 23, 385–419. [Google Scholar] [CrossRef]
  16. Dick, C.P.; Korkolis, Y.P. Mechanics and Full-Field Deformation Study of the Ring Hoop Tension Test. Int. J. Solids Struct. 2014, 51, 3042–3057. [Google Scholar] [CrossRef] [Green Version]
  17. Dick, C.P.; Korkolis, Y.P. Anisotropy of Thin-Walled Tubes by a New Method of Combined Tension and Shear Loading. Int. J. Plast. 2015, 71, 87–112. [Google Scholar] [CrossRef] [Green Version]
  18. Kamerman, D.; Cappia, F.; Wheeler, K.; Petersen, P.; Rosvall, E.; Dabney, T.; Yeom, H.; Sridharan, K.; Ševeček, M.; Schulthess, J. Development of Axial and Ring Hoop Tension Testing Methods for Nuclear Fuel Cladding Tubes. Nucl. Mater. Energy 2022, 31, 101175. [Google Scholar] [CrossRef]
  19. Venkatsurya, P.K.C.; Jia, Z.; Misra, R.D.K.; Mulholland, M.D.; Manohar, M.; Hartmann, J.E. Understanding Mechanical Property Anisotropy in High Strength Niobium-Microalloyed Linepipe Steels. Mater. Sci. Eng. A Struct. Mater. Prop. Microstruct. Process. 2012, 556, 194–210. [Google Scholar] [CrossRef]
  20. Korkolis, Y.P.; Kyriakides, S. Path-Dependent Failure of Inflated Aluminum Tubes. Int. J. Plast. 2009, 25, 2059–2080. [Google Scholar] [CrossRef]
  21. Yoshida, K.; Kuwabara, T. Effect of Strain Hardening Behavior on Forming Limit Stresses of Steel Tube Subjected to Nonproportional Loading Paths. Int. J. Plast. 2007, 23, 1260–1284. [Google Scholar] [CrossRef]
  22. Korkolis, Y.; Kyriakides, S. Inflation and Burst of Anisotropic Aluminum Tubes for Hydroforming Applications. Int. J. Plast. 2008, 24, 509–543. [Google Scholar] [CrossRef]
  23. Kuwabara, T.; Sugawara, F. Multiaxial Tube Expansion Test Method for Measurement of Sheet Metal Deformation Behavior under Biaxial Tension for a Large Strain Range. Int. J. Plast. 2013, 45, 103–118. [Google Scholar] [CrossRef]
  24. He, Z.B.; Zhang, K.; Lin, Y.L.; Yuan, S.J. An Accurate Determination Method for Constitutive Model of Anisotropic Tubular Materials with Dic-Based Controlled Biaxial Tensile Test. Int. J. Mech. Sci. 2020, 181, 105715. [Google Scholar] [CrossRef]
  25. Zhang, K.; He, Z.B.; Zheng, K.; Yuan, S.J. Experimental Verification of Anisotropic Constitutive Models under Tension-Tension and Tension-Compression Stress States. Int. J. Mech. Sci. 2020, 178, 105618. [Google Scholar] [CrossRef]
  26. Wang, X.; Hu, W.; Huang, S.; Ding, R. Experimental Investigations on Extruded 6063 Aluminium Alloy Tubes under Complex Tension-Compression Stress States. Int. J. Solids Struct. 2019, 168, 123–137. [Google Scholar] [CrossRef]
  27. Nazari Tiji, S.A.; Park, T.; Asgharzadeh, A.; Kim, H.; Athale, M.; Kim, J.H.; Pourboghrat, F. Characterization of Yield Stress Surface and Strain-Rate Potential for Tubular Materials Using Multiaxial Tube Expansion Test Method. Int. J. Plast. 2020, 133, 102838. [Google Scholar] [CrossRef]
  28. Parmar, A.; Mellor, P.B. Plastic Expansion of a Circular Hole in Sheet Metal Subjected to Biaxial Tensile Stress. Int. J. Mech. Sci. 1978, 20, 707–720. [Google Scholar] [CrossRef]
  29. Choi, Y.; Ha, J.; Lee, M.G.; Korkolis, Y.P. Effect of Plastic Anisotropy and Portevin-Le Chatelier Bands on Hole-Expansion in Aa7075 Sheets in -T6 and -W Tempers. J. Mater. Process. Technol. 2021, 296, 117211. [Google Scholar] [CrossRef]
  30. Ha, J.; Coppieters, S.; Korkolis, Y.P. On the Expansion of a Circular Hole in an Orthotropic Elastoplastic Thin Sheet. Int. J. Mech. Sci. 2020, 182, 105706. [Google Scholar] [CrossRef]
  31. Barlat, K.F. Numerical Modeling for Accurate Prediction of Strain Localization in Hole Expansion of a Steel Sheet. Int. J. Solids Struct. 2019, 156–157, 107–118. [Google Scholar]
  32. Suzuki, T.; Okamura, K.; Capilla, G.; Hamasaki, H.; Yoshida, F. Effect of Anisotropy Evolution on Circular and Oval Hole Expansion Behavior of High-Strength Steel Sheets. Int. J. Mech. Sci. 2018, 146–147, 556–570. [Google Scholar] [CrossRef]
  33. Liu, W.; Iizuka, T. Trials to Evaluate Distribution of Lankford Value Using Hole Expansion Test. Procedia Manuf. 2018, 15, 1754–1761. [Google Scholar] [CrossRef]
  34. Hill, R. A Theory of the Yielding and Plastic Flow of Anisotropic Metals. Proc. R. Soc. Lond. 1948, 193, 281–297. [Google Scholar] [CrossRef] [Green Version]
  35. GB/T 228.1-2010; Metallic Materials-Tensile Testing-Part 1: Method of Test at Room Temperature. Standards Press of China: Beijing, China, 2011.
Figure 1. Schematic of the HBT.
Figure 1. Schematic of the HBT.
Materials 16 04629 g001
Figure 2. (a) Stress state of single-layer tube during free bulging, and (b) Stress state of the HBT.
Figure 2. (a) Stress state of single-layer tube during free bulging, and (b) Stress state of the HBT.
Materials 16 04629 g002
Figure 3. Flow-chart of hybrid numerical–experimental method.
Figure 3. Flow-chart of hybrid numerical–experimental method.
Materials 16 04629 g003
Figure 4. Tube free bulging test equipment.
Figure 4. Tube free bulging test equipment.
Materials 16 04629 g004
Figure 5. (a) Specimen sampling diagram, and (b) Specimen dimensions (unit mm).
Figure 5. (a) Specimen sampling diagram, and (b) Specimen dimensions (unit mm).
Materials 16 04629 g005
Figure 6. (a) True stress–strain curve, and (b) Evolution of the r-value (r0) with strain.
Figure 6. (a) True stress–strain curve, and (b) Evolution of the r-value (r0) with strain.
Materials 16 04629 g006
Figure 7. (a) Specimen dimensions (unit mm), (b) Clamps, and (c) Test device.
Figure 7. (a) Specimen dimensions (unit mm), (b) Clamps, and (c) Test device.
Materials 16 04629 g007
Figure 8. Specimens after hoop tensile tests: (a) Front view, and (b) Side view.
Figure 8. Specimens after hoop tensile tests: (a) Front view, and (b) Side view.
Materials 16 04629 g008
Figure 9. FE model for the HBT.
Figure 9. FE model for the HBT.
Materials 16 04629 g009
Figure 10. FE mesh: (a) the entire model, and (b) Local mesh near the hole.
Figure 10. FE mesh: (a) the entire model, and (b) Local mesh near the hole.
Materials 16 04629 g010
Figure 11. Stress distribution at the hole periphery for different bulging heights: (a) 1.0 mm, (b) 1.5 mm, (c) 2.0 mm, and (d) 2.5 mm.
Figure 11. Stress distribution at the hole periphery for different bulging heights: (a) 1.0 mm, (b) 1.5 mm, (c) 2.0 mm, and (d) 2.5 mm.
Materials 16 04629 g011aMaterials 16 04629 g011b
Figure 12. Strain contour diagrams of hole periphery with different bulging heights.
Figure 12. Strain contour diagrams of hole periphery with different bulging heights.
Materials 16 04629 g012
Figure 13. Strain distribution at the hole periphery for different bulging heights: (a) 1.0 mm, (b) 1.5 mm, (c) 2.0 mm, and (d) 2.5 mm.
Figure 13. Strain distribution at the hole periphery for different bulging heights: (a) 1.0 mm, (b) 1.5 mm, (c) 2.0 mm, and (d) 2.5 mm.
Materials 16 04629 g013
Figure 14. The relationship between aperture and bulging height for different r45: (a) Overall graph, and (b) Local zoomed-in graph.
Figure 14. The relationship between aperture and bulging height for different r45: (a) Overall graph, and (b) Local zoomed-in graph.
Materials 16 04629 g014
Figure 15. Relationship between aperture and r45 at the same bulging height.
Figure 15. Relationship between aperture and r45 at the same bulging height.
Materials 16 04629 g015
Figure 16. Thickness of hole periphery at 2.5 mm bulging height.
Figure 16. Thickness of hole periphery at 2.5 mm bulging height.
Materials 16 04629 g016
Figure 17. Variation of thickness with bulging height at different locations of hole periphery: (a) 0°, (b) 45°, and (c) 90°.
Figure 17. Variation of thickness with bulging height at different locations of hole periphery: (a) 0°, (b) 45°, and (c) 90°.
Materials 16 04629 g017
Figure 18. Effect of different inner and outer tube strength relationships on the shape of the hole after deformation.
Figure 18. Effect of different inner and outer tube strength relationships on the shape of the hole after deformation.
Materials 16 04629 g018
Figure 19. Thickness distribution at ultimate bulging height for different initial aperture.
Figure 19. Thickness distribution at ultimate bulging height for different initial aperture.
Materials 16 04629 g019
Figure 20. Tubular specimens after HB experiments with initial apertures of: (a) 2 mm, and (b) 4 mm.
Figure 20. Tubular specimens after HB experiments with initial apertures of: (a) 2 mm, and (b) 4 mm.
Materials 16 04629 g020
Figure 21. Schematic diagram of the tubular specimen with a hole before experiment.
Figure 21. Schematic diagram of the tubular specimen with a hole before experiment.
Materials 16 04629 g021
Figure 22. Experimental results of 3.31 mm bulging height.
Figure 22. Experimental results of 3.31 mm bulging height.
Materials 16 04629 g022
Figure 23. Test results of the HB experiments.
Figure 23. Test results of the HB experiments.
Materials 16 04629 g023
Figure 24. The bulging height–aperture curve of tubular specimen with a hole.
Figure 24. The bulging height–aperture curve of tubular specimen with a hole.
Materials 16 04629 g024
Figure 25. Comparison of aperture between HB experiments and simulation results for r45 = 1.06.
Figure 25. Comparison of aperture between HB experiments and simulation results for r45 = 1.06.
Materials 16 04629 g025
Figure 26. Comparison of thicknesses at the hole periphery from experiment and simulation with r45 = 1.06.
Figure 26. Comparison of thicknesses at the hole periphery from experiment and simulation with r45 = 1.06.
Materials 16 04629 g026
Figure 27. Comparison of experimental and simulated contour of the bulging zone for r45 = 1.06: (a) Distribution of sampling points, and (b) Comparison of the contour of the bulging zone.
Figure 27. Comparison of experimental and simulated contour of the bulging zone for r45 = 1.06: (a) Distribution of sampling points, and (b) Comparison of the contour of the bulging zone.
Materials 16 04629 g027
Figure 28. AA6061-O extruded aluminum alloy tubes in-plane anisotropy coefficients.
Figure 28. AA6061-O extruded aluminum alloy tubes in-plane anisotropy coefficients.
Materials 16 04629 g028
Table 1. Effect of strength relationship between inner and outer tube on hole deformation.
Table 1. Effect of strength relationship between inner and outer tube on hole deformation.
Strength RelationshipHISESLIS
Ultimate bulging height/mm3.203.112.95
Aperture/mm3.603.052.97
Deformation difference degree/%80.052.548.5
Table 2. Simulation results for different initial apertures.
Table 2. Simulation results for different initial apertures.
Initial Aperture/mm2.03.04.05.0
Ultimate bulging height/mm3.253.183.102.89
Aperture/mm3.665.306.757.80
Deformation difference degree/%83.076.768.856.0
Table 3. Aperture after deformation with different friction coefficients.
Table 3. Aperture after deformation with different friction coefficients.
Friction Coefficient0.10.20.3
r45 = 0.55.58 mm5.58 mm5.60 mm
r45 = 1.05.82 mm5.80 mm5.83 mm
r45 = 1.56.08 mm6.01 mm6.03 mm
Table 4. Test results of the HBT at different bulging heights.
Table 4. Test results of the HBT at different bulging heights.
Bulging Height/mmInternal Pressure/MPaAperture/mm
2.1128.075.50
2.4129.345.73
2.5129.685.84
2.7130.026.05
2.9430.326.35
3.1530.516.58
3.3130.70Rupture
Table 5. Iterative results.
Table 5. Iterative results.
Iterative MethodsNumber of IterationsIteration Step Size r45 RMSE   η ¯
Bisection method 0.91.47%
1.10.53%
10.11.00.68%
20.051.050.53%
Fixed step method3-10.011.060.47%
3-20.011.070.51%
3-30.011.080.60%
3-40.011.090.64%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lin, Y.; Wang, Y.; Su, Y.; Liu, J.; Chen, K.; He, Z. Identification of Anisotropic Coefficients in the Non-Principal Axis Directions of Tubular Materials Using Hole Bulging Test. Materials 2023, 16, 4629. https://doi.org/10.3390/ma16134629

AMA Style

Lin Y, Wang Y, Su Y, Liu J, Chen K, He Z. Identification of Anisotropic Coefficients in the Non-Principal Axis Directions of Tubular Materials Using Hole Bulging Test. Materials. 2023; 16(13):4629. https://doi.org/10.3390/ma16134629

Chicago/Turabian Style

Lin, Yanli, Yifan Wang, Yibo Su, Junpeng Liu, Kelin Chen, and Zhubin He. 2023. "Identification of Anisotropic Coefficients in the Non-Principal Axis Directions of Tubular Materials Using Hole Bulging Test" Materials 16, no. 13: 4629. https://doi.org/10.3390/ma16134629

APA Style

Lin, Y., Wang, Y., Su, Y., Liu, J., Chen, K., & He, Z. (2023). Identification of Anisotropic Coefficients in the Non-Principal Axis Directions of Tubular Materials Using Hole Bulging Test. Materials, 16(13), 4629. https://doi.org/10.3390/ma16134629

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop