Skip to content
BY 4.0 license Open Access Published by De Gruyter September 2, 2020

Optimization of spinning parameters of 20/316L bimetal composite tube based on orthogonal test

  • Wenbin Xu EMAIL logo , Yingying Yang , Chenwei Dai and Jingang Xie

Abstract

The spin forming parameters of bimetallic composite pipe were optimized by FEM and experiment, and the torque and the residual contact pressure were analyzed during forming process. Adoption of The orthogonal test of four factors and five levels to gain the optimal spinning parameters (δ = 0.16mm, Ψ = 0.4mm, β = 3.0, f = 0.5) and the effects of factors on analysis indicators under the condition of maximum residual contact pressure and small torque. The simulation results showed that the residual contact pressure was nearly 1MPa higher and the torque was lower about 160N·m than unoptimized results. It proved that the orthogonal test method was effective for the optimization of spinning parameters of bimetal composite tube. The accuracy of numerical simulation value was verified by the drawing experiment of the composite tube.

1 Introduction

Bimetal composite tube uses carbon steel tube as the base to ensure its excellent comprehensive mechanical properties, and according to the requirement with special performance material as the liner to meet the special requirements [1, 2, 3, 4]. Bimetal composite tubes fully use the base and liner in the best physical and chemical properties. It has high strength but also excellent corrosion resistance, high-temperature resistance, low-temperature resistance, wear-resistance and other special performance [5, 6, 7]. As a result, bimetal composite tube can be widely used in petroleum, chemical and other fields.

The forming technology of metal composite pipe mainly includes hydroforming, drawing forming and explosion forming. Due to the dual effects of physical non-linearity and geometric nonlinearity, it is difficult to study on the theoretical of the bimetal composite tube [8]. With the development of computer technology, the finite element numerical simulation method has been widely used in metal plastic forming. X.F.Lu and X. Zheng simulated the drawing forming process of bimetal-lined pipe made of 20 carbon steel as the base tube and 316L stainless steel as the liner tube by means of the finite element analysis software ABAQUS, optimized the drawing process parameters of bimetal composite tube, and obtained the optimal drawing process parameters at the maximum residual contact pressure [9]. Drawing forming of composite pipe also has been reported in other literatures [10, 11]. They analyzed the stress distribution in the drawing process of the composite pipe and studied the influence of drawing parameters on the bonding strength of the composite pipe by orthogonal test and got a modified approach for drawing force estimation of bi-metallic composite pipe forming in the light of the results of FEM simulation. Z.Y. Bi and C.Wang analyzed the minimum expansion pressure and the relationship between the expansion pressure and residual contact stress for the L360QS/Incoloy825 nickel base alloy hydraulic expansion pipe. They pointed out that the minimum expansion pressure and the residual contact pressure Pc* can be obtained by using ABAQUS finite element method and verified by hydraulic bulge test [12]. Other authors also investigated the hydroforming of composite tube by experimental or finite element [13, 14, 15, 16]. They studied the influence factors of hydraulic formation failure of composite pipe by experiment, and examined the interaction of material properties and die geometry on the selection of hydroforming process parameters and discussed the effects of the strain-hardening exponent, initial tube wall thickness, and die corner radii on corner filling and thickness distribution of the hydroformed tube. B. Cao and Z.S.Zhang introduced ANSYS finite element software to analyze the composite dynamic process, and obtained the displacement of key nodes of the liner and the stress distribution of the composite tube after forming [17]. Yu yong et al.applied PVDF piezoelectric thin film sensor to measure the dynamic parameters of the explosive welding and studied the impact pressure of the aluminum layer on the steel layer in the aluminum-steel underwater explosive [18]. The explosive forming of metal composite tubes also has been reported in other literatures [19, 20, 21]. They studied the effect of the explosive mass reduction on the bonding interface of the pipes and resolved the backward method of charge and low energy efficiency of explosives by using honeycomb structure explosives and double sided explosive cladding. T. Nechaykina et al. showed the viability of manufacturing of multi-layer steel/vanadium alloy/steel composite tubes by numerical simulations and experiment, studied the interface bonding of vanadium alloy and stainless steel by using scanning electron microscopy and microhardness measurements, and obtained the process conditions suitable for producing a good quality product [22]. Pawel Kazanowski et al. investigated the influence of the initial billet geometry on the product geometry during bi-material rod extrusion, optimized initial inner material length within a bi-material billet and noted the complex interface structure between AA2014 and AA6063 by light optical microscopy [23]. However, if the composite tube formed by drawing encounters high temperature, it would tend to be stratified, and the tube would fail due to the release of residual contact pressure. The forming force of hydroforming is uniform and can be calculated, but the device is pretty complicated and high requirements for the sealing of the device. The composite tube forming by Explosive is difficult to accurately determine the amount of explosive, the interface junction between base tube and liner is easy to form waveform and scrap, and special production site is needed and the production process has certain risks. The spinning forming equipment is simple, and the forming quality of the composite tube is relatively good. However, few literatures reported the FEM for the spinning process of composite tubes and optimized its forming parameters. This research mainly aimed at establishing the numerical simulation model of bimetal composite tube formed by spinning, analyzing the equivalent plastic strain, applying the orthogonal experiment design method and range difference to get the significant parameter of bimetal composite tube, and using experiments to verify the accuracy of the simulation.

2 Finite element analysis of spinning forming of compound tube

2.1 Build finite element model

The size of the carbon steel base tube is Φ 89×4.5mm and 316L liner size is Φ 79.2×1 mm. The gap between base and liner tube is 0.4mm. The embedding quantity δ is 0.14 mm, orientation Angle β (the centerline of the roller and liner axis Angle) is 2.5, the feed rate f (axial feed roller speed/angular velocity) is 0.3. The three-dimensional simulation model of the bimetallic clad pipe is shown as the Figure 1.

Figure 1 3D model of bimetallic clad pipe formed by spinning
Figure 1

3D model of bimetallic clad pipe formed by spinning

2.2 Defining the load and boundary conditions

In the spinning process of the composite tube, the liner should be flanged and fixed on the base tube. Therefore, the displacement/Angle constraint mode is adopted when defining boundary conditions. U2 = 0 is defined at the end of the liner and base tube. The axial loading speed defined as v2 = 18.84mm/s, the rotating speed n = 600r/min and the amplitude curve of smooth step is introduced for loading.

2.3 Contact mode and mesh division

Applying the method of surface to surface contact between the outer wall of the liner and the inner wall of the base tube and assuming the inner wall of the base tube as the main contact surface and the outer wall of the liner as the slave contact surface simulated the forming process. The surface to surface contact method was also adopted between the roller and the inner wall of the liner inner wall, and the roller as the main contact surface and the inner wall of the liner as the slave contact surface. The finite slip was introduced in each contact pair and the motion contact as the mechanical constraint. The C3D8R hexahedron unit was employed to mesh the liner and base tube, the liner was divided into 94000 hexahedron units and 79000 hexahedron units for the base tube was divided. Because of the shear load, the grid distortion may occur in the liner and base tube and it seriously affects the analytical accuracy, misconvergence and even interrupts the program. The usage of the ALE adaptive grid control prevents the occurrence of this phenomenon.

3 Simulation results

3.1 Equivalent plastic strain analyses

The equivalent plastic strain (PEEQ) of the base tube and liner in the stable spinning stage of the composite tube was shown in Figure 2:

Figure 2 Plastic strain of base tube and lined tube
Figure 2

Plastic strain of base tube and lined tube

The material hardening model was used to calculate the accumulated plastic strain, and the PEEQ can be used to express the accumulated plastic strain. According to the shear failure criterion, the material failure may happen when the PEEQ reached above 1.0. Figure 2 (a) showed that the equivalent plastic strain of the base tube was almost close to zero. It can be considered that the base tube was always within the elastic range during the forming process, the theoretical calculation with multiple tube spinning forming and the actual production safety requirements were consistent. The liner would be plastic deformation because of its large deformation. As it can be seen from Figure 2 (b) of the equivalent plastic strain distribution of the liner, plastic deformation occurred in the liner, and the maximum PEEQ was 0.2279, which didn’t exceed 1.0, so the material didn’t fail due to large deformation.

3.2 Residual contact pressure and torque

Residual contact pressure is the main index for the quality of composite tube. Under the condition of the base tube being in plastic deformation and the liner being well, the greater the residual contact pressure between the liner and the base tube, the better the quality of the composite tube. Figure 3 showed the distribution of the average residual contact pressure obtained along the uniformly selected path along the inner wall of the base tube. It can be seen that the residual contact pressure of each path was inconsistent, and the average residual contact pressure of each path was 8.20MPa.Metal accumulation formed by uneven metal-flowing in the forming process would lead to uneven residual contact pressure. In practical production, physical defects, such as small concave anti-bulge on inner wall of base tube or outer wall of liner, would also lead to uneven distribution of contact pressure between liner and base tube.

Figure 3 Residual contact pressure between lined tube and base tube
Figure 3

Residual contact pressure between lined tube and base tube

The torque distribution of the composite tube in spinning process was shown in Figure 4. The forming process went through three stages: in the initial spinning stage, the spinning mold came into contacting with the liner, the metal-flowing started to increase, the contact pressure also increased gradually, and the torque increased with time. In the stable spinning stage, the roller was in full contact with the liner, and the metal-flowing and contact pressure tended to be stable, so the torque fluctuated in a small range and tended to be stable. Spinning unloading stage, spinning roller and liner separated gradually, contact pressure reduced and torque also reduced to zero. The average torque of the entire spinning process was 3943.69 N·m.

Figure 4 Torque
Figure 4

Torque

4 Orthogonal experimental analysis of compound tube spinning forming

In order to get high bonding strength (residual contact pressure) and low required torque and reduce the number of numerical simulation as much as possible, the following four factors were selected: A-embedding amount δ, B-liner and base tube gap Ψ, C-guiding Angle β, D-feed rate f, and 5 levels for each factor. The orthogonal table selected according to the horizontal number was L25 (56), and the orthogonal test scheme and orthogonal test factor level developed were shown in the table below. Finite element analysis software ABAQUS was used to conduct numerical simulation on the spinning process parameters of 25 groups in Table 1.

Table 1

Scheme of orthogonal experiment

Test numberABCDPrc/MPaTorque/N·m
10.100.303.00.32.933219.61
20.100.353.50.353.163419.32
30.100.404.00.44.563559.63
40.100.454.50.454.743681.15
50.100.505.00.55.203873.26
60.120.303.50.45.633546.46
70.120.354.00.456.533637.29
80.120.404.50.57.423773.03
90.120.455.00.37.094186.13
100.120.503.00.354.833548.25
110.140.304.00.58.333748.13
120.140.354.50.37.204105.91
130.140.405.00.358.204305.54
140.140.453.00.46.263715.14
150.140.503.50.457.393845.23
160.160.304.50.358.744173.54
170.160.355.00.49.134354.74
180.160.403.00.458.083860.67
190.160.453.50.58.474001.76
200.160.504.00.39.114267.37
210.180.305.00.459.274457.39
220.180.353.00.58.443939.09
230.180.403.50.37.734253.25
240.180.454.00.357.834373.62
250.180.504.50.47.514445.12

4.1 Range analysis

The range analysis method can estimate the influence of various factors on the test results. The larger the range R value is, the greater the influence of the corresponding factors will be. Therefore, the influence of each factor on the test evaluation index can be ensured by the size of the range R value, and the better level of each factor can be ensured according to the size of the k value, so as to determine the optimal level combination.

From Table 3, The range analysis of the four factors on the residual contact pressure was RA>RC>RD>RB, which meant that the embedding amount exerted the greatest influence on the residual contact pressure, followed by the leading Angle, the feed rate, and the liner and base tube clearance, respectively. The four factors of A, B, C and D respectively gained the maximum residual contact pressure at the level of A4, B3, C5 and D5, so the corresponding optimal spinning process parameter was A4B3C5D5.

Table 2

Factors and levels of orthogonal experiment

levelfactA δ/(mm)B Ψ/(mm)C β()D f
level 10.100.303.00.3
level 20.120.353.50.35
level 30.140.404.00.4
level 40.160.454.50.45
level 50.180.505.00.5
Table 3

Residual contact pressure of orthogonal results

Residual contact pressure Prc/MPa
knABCD
k14.126.986.116.81
k26.306.896.486.55
k37.487.207.276.62
k48.716.887.127.20
k58.166.817.787.57
R4.590.391.671.02
OptimalA4B3C5D5

It can be seen from table 4 that the range analysis of torque by the four factors was RA>RC>RB>RD. It was also basically consistent with the range analysis results of residual contact pressure. The embedding amount and steering Angle had a large influence on the residual contact pressure and torque, and the gap between liner and base tube and the feed rate f had a low influence on them. The minimum torque was obtained at the level of A1, B1, C1 and D5, so the optimal process parameter should be A1B1C1D5. Overall consideration of torque and the <CJ/T 192-2004 lined stainless-steel tube> standard for the requirement of the residual contact pressure, the optimal parameter combination was acquired at A4B3C1D5.

4.2 Comparison of optimized and unoptimized results

The optimized parameters were used for numerical simulation and compared with unoptimized results. The differences between them were shown in Figure 5 and Figure 6.

Figure 5 Residual contact pressure
Figure 5

Residual contact pressure

Figure 6 Torque
Figure 6

Torque

Table 4

Torque of orthogonal results

Torque/N·m
knABCD
k13550.593829.033656.554006.45
k23738.233850.423813.203964.05
k33943.993891.273917.213924.22
k44131.623991.564035.753896.35
k54293.693995.854235.413867.05
R743.1145.435788.86139.4
OptimalA1B1C1D5

From Figure 5 (b), it can be seen that the range of residual contact pressure was form 8.4 MPa to 10 MPa, which was smaller than unoptimized results and the quality of the tube was more stable. The average residual contact pressure was 9.19 MPa, nearly 1 MPa higher than unoptimized results. From the Figure 6 (b), the average torque was 3783.61N·m, 160N·m less than unoptimized results.

5 Experimental analyses

The preparation of the test in accordance with Table 5 and composite tube samples were made to verify the reliability of the simulation. In the middle of the samples, the base tube carbon steel tube cut off 10mm, and the complete lined stainless-steel tube was left. The two Φ 20 mm holes were drilled at both ends of specimens to conveniently clamp for pull-out test. The resulting entity diagram was shown in Figure 7.

Table 5

Spinning test

Pilot programδΨβf
10.100.353.50.35
20.140.453.00.40
3(Optimization)0.160.403.00.50
40.180.504.50.40
Figure 7 The entity graph of bimetallic clad tube specimen
Figure 7

The entity graph of bimetallic clad tube specimen

5.1 Pull-out test

The holes on both ends of the specimen were clamped by the jacket on the material testing machine. One end was fixed and the other end could be moving. The testing machine applied tension to the removable jacket and pulled it out according to the standard of <CJ/T 192-2004 lined stainless-steel tube>. The force corresponding to the displacement of the liner and the base tube under the action of tension was the binding force F at the interface of the liner and the base tube separation, from the formula Prc = F/πμDiH, where Di: inner diameter of the liner base tube of the composite tube. H: effective binding length of composite tube sample. Friction coefficient of liner and base tube, μ = 0.21.

5.2 Experimental scheme

According to the influence factor of residual contact pressure and the limitation of test conditions, four groups of tests were conducted, as shown in Table 5. The test parameters in the test table were selected according to the orthogonal test and test conditions.

5.3 Experimental results

From Figure 8, the residual contact pressure by simulating was not much different from the experimental value, with the mutual error of about 15%, under the acceptable range of engineering, because the simulated value was obtained under the assumption of the liner and base tube being ideal, and it would cause the simulated value larger than the test value. Moreover, the concave resistance, bending, uneven wall thickness and poor cleaning of liner and base tube during the actual forming of composite tube would reduce the residual contact stress.

Figure 8 The comparison of residual contact pressure
Figure 8

The comparison of residual contact pressure

6 Conclusions

Through numerical simulation and experimental study on spinning forming of bimetallic composite pipe, the following conclusions are obtained:

  1. The numerical simulation model of bimetal composite tube formed by spinning was established. Experimental value and simulation value were compared and it proved the feasibility of the numerical simulation method for spinning forming of the metal composite tube.

  2. The residual contact pressure and the equivalent plastic strain of liner and base tube and the required torque were obtained. It provided the technical basis for the practical production of the composite tube.

  3. The influence of the embedding amount on the residual contact pressure was the largest, followed by the guiding angle, the feeding rate f and the liner and base tube clearance.

  4. The optimal spinning forming technology parameters (δ = 0.16mm, = 0.4mm, β = 3.0, f = 0.5), the residual contact pressure nearly 1 MPa higher and the required torque reduced 160N*m than that unoptimized results.

Acknowledgement

This work was supported by the Fundamental Research Funds for Anhui vocational and technical college, Anhui Provincial Natural Science Foundation (No. KJ2018B0003).

References

[1] Scott DA, Wolgemuth GA, Aikin JA. Hydraulically expanded tube-to-tube sheet joints. J Press Vessel Technol. 1984;106(1):104–9.10.1115/1.3264297Search in Google Scholar

[2] Berski S, Dyja H, Maranda A, Nowaczewski J, Banaszek G. Analysis of quality of bimetallic rod after extrusion process. J Mater Process Technol. 2006;177(1-3):582–6.10.1016/j.jmatprotec.2006.04.107Search in Google Scholar

[3] Guo X, Tao J, Wang W, Li H, Wang C. Effects of the inner mould material on the aluminium-316L stainless steel explosive clad tube. Mater Des. 2013;49:116–22.10.1016/j.matdes.2013.02.001Search in Google Scholar

[4] Krishnamurthy R. Hot Metal Gas Forming of Magnesium Tubes. Compos. Part. International Conference on Physical&Numerical Simulation of Materials Processing, 2004,34(2): 15-17.Search in Google Scholar

[5] Khosravifard A, Ebrahimi R. Investigation of parameters affecting interface strength in Al/Cu clad bimetal rod extrusion process. Mater Des. 2010;31(1):493–9.10.1016/j.matdes.2009.06.026Search in Google Scholar

[6] Alcaraz JL, Gil-Sevillano J. An analysis of the extrusion of bimetallic tubes by numerical simulation. J.Mech.Sci. 2016;38(2):157– 73.10.1016/0020-7403(95)00044-XSearch in Google Scholar

[7] MahjabinTaskin, GinshiroUtsumi, YoshimineKato. Observation of ultrasonic signal and measure-ment of H2 concentration from the exterior of a metal pipe. J Hydrogen Energy. 2019, 44(41): 03-12 .10.1016/j.ijhydene.2019.06.159Search in Google Scholar

[8] Krajewski PE, Morales AT. Tribological Issues During Quick Plastic Forming. J Mater Eng Perform. 2004;13(6):701–6.10.1361/10599490421330Search in Google Scholar

[9] Lu XF, Zheng X. Optimization of drawing parameters for 20/316L bimetal-lined tube based on FEM simulation. J Chinese Journal of Nonferrous Metals. 2011;21(1):205–13.Search in Google Scholar

[10] Zheng M, Yu J. A modified approach for drawing force estimation of bi-metallic composite pipe Forming. Compos. Part. ICMMBE. 2016;83:270–5.10.2991/icmmbe-16.2016.53Search in Google Scholar

[11] Yang XJ, Sun FH, Zhang ZM, et al. Optimization of drawing parameters for copper tubes with hollow sinking based on FEM simulation. J Nonferrous Metals. 2008;18(12):2245–52.Search in Google Scholar

[12] Bi ZY, Wang C. Finite element simulation and experimental study of L360QS/Incoloy825 nickel base alloy hydraulic tube. J. Plast Eng. 2016;23(1):131–5.Search in Google Scholar

[13] Olabi AG, Alaswad A. Experimental and finite element investigation of formability and failures in bi-layered tube hydroforming. Adv Eng Softw. 2011;42(10):815–20.10.1016/j.advengsoft.2011.05.022Search in Google Scholar

[14] Islam MD, Olabi AG, Hashmi MS. Feasibility of multi-layered tubular components forming by hydroforming and finite element simulation. J Mater Process Technol. 2006;174(1-3):394–8.10.1016/j.jmatprotec.2006.02.016Search in Google Scholar

[15] Hossein SM. Seyedkashi, Panahizdeh R. Process Analysis of Two-layered Tube Hydroformingwith Analytical and Experimental verification. J Mech Sci Technol. 2013;27(1):169–75.10.1007/s12206-012-1216-7Search in Google Scholar

[16] Kridli GT, Bao L,Mallick PK, Tian Y. Investigation of Thickness Variation and Corner Filling in Tube Hydroforming. J Mater Process Technol. 2003;133(3):287–96.10.1016/S0924-0136(02)01004-XSearch in Google Scholar

[17] Cao JB, Zhang ZS. Finite element analysis of underwater explosive forming of bimetal composite tube. Transactions of the China Welding Institution. 2018;39(2):61–5.Search in Google Scholar

[18] Yu Y, Ma HH, Zhao K. Study on Explosive Welding of Al-Steel Coaxial Pipes Underwater. Cent Eur J Energ Mater. 2017;14(1):251–6.10.22211/cejem/68696Search in Google Scholar

[19] Zamani E, Liaghat GH. Explosive welding of stainless steel–carbonsteel coaxial pipes. J Mater Sci. 2012;47(2):685–95.10.1007/s10853-011-5841-9Search in Google Scholar

[20] Sui GF, Li JS, Li HW, Sun F, Zhang T, Fu H. Investigation on the explosive welding mechanism of corrosion-resisting aluminium and stainless steel tubes through finite element simulation and experiments. Int J Miner Metall Mater. 2012;19(2):151–8.10.1007/s12613-012-0531-4Search in Google Scholar

[21] Miao GH, Ma HH, Shen ZW, Yu Y. Research on honeycomb structure explosives and double sided explosive cladding. Mater Des. 2014;63:538–43.10.1016/j.matdes.2014.06.050Search in Google Scholar

[22] Nechaykina T, Nikulin S, Rozhnov A, Molotnikov A, Zavodchikov S, Estrin Y. Proving the viability of manufacturing of multi-layer steel/vanadium alloy/steel composite tubes by numerical simulations and experiment. J Nucl Mater. 2018;503:178–90.10.1016/j.jnucmat.2018.02.043Search in Google Scholar

[23] Kazanowski P, Mario E, Wojciech Z. Bi-metal rod extrusion-process and product optimization. Mater Sci Eng. 2004;369(1-2):170–80.10.1016/j.msea.2003.11.002Search in Google Scholar

Received: 2019-07-18
Accepted: 2019-12-02
Published Online: 2020-09-02

© 2020 W. Xu et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 25.4.2024 from https://www.degruyter.com/document/doi/10.1515/secm-2020-0026/html
Scroll to top button