Bending behaviors of large diameter thin-walled CP-Ti tube in rotary draw bending

Similar documents
Effect of Mandrel, Its Clearance and Pressure Die on Tube Bending Process via Rotary Draw Bending Method

A study of forming pressure in the tube-hydroforming process

An Analytical Model for Long Tube Hydroforming in a Square Cross-Section Die Considering Anisotropic Effects of the Material

Optimization of blank dimensions to reduce springback in the flexforming process

INVERSE METHOD FOR FLOW STRESS PARAMETERS IDENTIFICATION OF TUBE BULGE HYDROFORMING CONSIDERING ANISOTROPY

An alternative design method for the double-layer combined die using autofrettage theory

THE OPTIMIZATION OF SEVERAL TECHNOLOGICAL PARAMETERS OF ULTRASONIC VIBRATION AIDED FLOW-FORMING

For ASME Committee use only.

Prediction of Elastic Constants on 3D Four-directional Braided

Inverse method for flow stress parameters identification of tube bulge hydroforming considering anisotropy

Influence of impact velocity on transition time for V-notched Charpy specimen*

Analysis of axisymmetric cup forming of metal foil and micro hydroforming process

STANDARD SAMPLE. Reduced section " Diameter. Diameter. 2" Gauge length. Radius

Burst pressure estimation of reworked nozzle weld on spherical domes

Prediction of geometric dimensions for cold forgings using the finite element method

Carrying Capacity of Pressure Vessels under Hydrostatic Pressure

Expansion of circular tubes by rigid tubes as impact energy absorbers: experimental and theoretical investigation

The objective of this experiment is to investigate the behavior of steel specimen under a tensile test and to determine it's properties.

The Ultimate Load-Carrying Capacity of a Thin-Walled Shuttle Cylinder Structure with Cracks under Eccentric Compressive Force

Effect of Forming Conditions on Bending Characteristics in Push-Through Bending of Aluminum Extrusion Sections

Development and analysis for the large-scale tee-forming process

The numerical simulation research of an Ultra-Light Photovoltaic Cell multilayer composite structure

INFLUENCE OF A WELDED PIPE WHIP RESTRAINT ON THE CRITICAL CRACK SIZE IN A 90 BEND

CHAPTER 6 MECHANICAL PROPERTIES OF METALS PROBLEM SOLUTIONS

Advanced Structural Analysis EGF Cylinders Under Pressure

Simulation of the effect of DIE Radius on Deep Drawing Process

Plane Strain Test for Metal Sheet Characterization

The Nonlinear Dynamic Behavior Analysis and Research for Sheet Steel Stamping Forming Problems Based on Finite Element Method

4. Objectives of Research work

SPRING-BACK PREDICTION FOR STAMPINGS FROM THE THIN STAINLESS SHEETS

Strength Study of Spiral Flexure Spring of Stirling Cryocooler

G1RT-CT D. EXAMPLES F. GUTIÉRREZ-SOLANA S. CICERO J.A. ALVAREZ R. LACALLE W P 6: TRAINING & EDUCATION

End forming of thin-walled tubes

Effects of Forming Conditions of Roll Offset Method on Sectional Shape at the Corner of Square Steel Pipe +

Simplified Method for Mechanical Analysis of Safety Class 1 Piping

A novel technique of friction and material property measurement by tip test in cold forging

Analysis of Tube End Formability of AA 2024 Tubes Using FEM

Modelling of ductile failure in metal forming

3D Finite Element Analysis of Drilling of Ti-6Al-4V Alloy

Neural Network Prediction of Nonlinear Elastic Unloading for High Strength Steel

Failure analysis of serial pinned joints in composite materials

Finite element simulation of residual stresses in laser heating

Mechanical Properties of Polymer Rubber Materials Based on a New Constitutive Model

Predicting the shear strength of round clinched joint

Laboratory 4 Bending Test of Materials

Outline. Tensile-Test Specimen and Machine. Stress-Strain Curve. Review of Mechanical Properties. Mechanical Behaviour

Improvement of Low Strain Pile Integrity Test

Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure

Mechanics of Materials MENG 270 Fall 2003 Exam 3 Time allowed: 90min. Q.1(a) Q.1 (b) Q.2 Q.3 Q.4 Total

TOUGHNESS OF PLASTICALLY-DEFORMING ASYMMETRIC JOINTS. Ford Research Laboratory, Ford Motor Company, Dearborn, MI 48121, U.S.A. 1.

Effect of various stress ratio parameters on cold upset forging of irregular shaped billets using white grease as lubricant

THE HYDRAULIC BULGE TEST AND ITS IMPORTANCE FOR THE VEGTER YIELD CRITERION

Kuo-Long LEE, Wen-Fung PAN and Chien-Min HSU

7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment

Finite element analysis of relationship between tightening torque and initial load of bolted connections

12/8/2009. Prof. A.K.M.B. Rashid Department of MME BUET, Dhaka

A Study on the Tube of Integral Propeller Shaft for the Rear-wheel Drive Automobile Using Carbon Composite Fiber

Stress-Strain Behavior

Grid supports design for dual-cooled fuel rods

Influence of the filament winding process variables on the mechanical behavior of a composite pressure vessel

Experimental investigation on monotonic performance of steel curved knee braces for weld-free beam-to-column connections

Influence of friction in material characterization in microindentation measurement

N = N A Pb A Pb. = ln N Q v kt. = kt ln v N

Centrifuge Shaking Table Tests and FEM Analyses of RC Pile Foundation and Underground Structure

ME Final Exam. PROBLEM NO. 4 Part A (2 points max.) M (x) y. z (neutral axis) beam cross-sec+on. 20 kip ft. 0.2 ft. 10 ft. 0.1 ft.

Effect Of The In-Situ Stress Field On Casing Failure *

Study of Pile Interval of Landslide Restraint Piles by Centrifuge Test and FEM Analysis

Numerical Simulation and Experimental Study of Electromagnetic Forming

Reference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity",

Arch. Metall. Mater. 62 (2017), 3,

Natural Boundary Element Method for Stress Field in Rock Surrounding a Roadway with Weak Local Support

Investigation into the Self-loosening Trend of Bolt Joints on the Tower Crane ShengChun Wang1,a, PeiWei Ni1,b, Ye Zhang1,c and MingXiao Dong1,2,d

After lecture 16 you should be able to

OPTIMIZATION OF MACHINING PARAMETERS USING DESIRABILITY FUNCTION ANALYSIS AND ANOVA FOR THERMO-MECHANICAL FORM DRILLING

INFLUENCE OF PROCESSING PARAMETERS ON FORMING QUALITY OF NON-CIRCULAR SPINNING

Initial Stress Calculations

Research of concrete cracking propagation based on information entropy evolution

An integrated approach to the design of high performance carbon fibre reinforced risers - from micro to macro - scale

NUMERICAL SIMULATION ON THERMAL CONTACT RESISTANCE OF PRACTICAL SOLID SURFACES

Experimental Study and Numerical Simulation on Steel Plate Girders With Deep Section

6.4 A cylindrical specimen of a titanium alloy having an elastic modulus of 107 GPa ( psi) and

Characterizations of Aluminum Alloy Sheet Materials Numisheet 2005

Module 5: Theories of Failure

Chapter 5 Torsion STRUCTURAL MECHANICS: CE203. Notes are based on Mechanics of Materials: by R. C. Hibbeler, 7th Edition, Pearson

Determination of Poisson s Ratio of Rock Material by Changing Axial Stress and Unloading Lateral Stress Test

Finite Element Simulation with Coupled Thermo-Mechanical Analysis of Superplastic Dieless Tube Drawing Considering Strain Rate Sensitivity*

Chapter 7. Highlights:

Mechanics of Materials Primer

2 Experiment of GFRP bolt

Journal of NUCLEAR SCIENCE and TECHNOLOGY, Vol. 41, No. 7, p (July 2004)

FME461 Engineering Design II

CHAPTER 4: BENDING OF BEAMS

STRESS UPDATE ALGORITHM FOR NON-ASSOCIATED FLOW METAL PLASTICITY

Numerical simulation the bottom structures. grounding test by LS-DYNA

MATERIALS. Why do things break? Why are some materials stronger than others? Why is steel tough? Why is glass brittle?

Fluid structure interaction dynamic analysis of a mixed-flow waterjet pump

FEA A Guide to Good Practice. What to expect when you re expecting FEA A guide to good practice

Computational Analysis for Composites

MECHANICAL BEHAVIOUR OF DAMAGED HDPE GEOMEMBRANES

Nonlinear Analytical Model for Wire Strands

Transcription:

Bending behaviors of large diameter thin-walled CP-Ti tube in rotary draw bending Zhi-yong ZHANG, He YANG, Heng LI, Ning REN, Yu-li TIAN State Key Laboratory of Solidification Processing, School of Materials Science and Engineering, Northwestern Polytechnical University, Xi an 710072, China Received 16 June 2011; Accepted 15 September 2011 Abstract: The rotary draw bending(rdb) of large diameter thin-walled (LDTW) tube needs the strict cooperation of multiple processing parameters to avoid possible multiple defects. Due to the specific properties, the bending behaviors of LDTW commercial pure titanium (CP-Ti) tube are much more complex to achieve precise deformation. With the CP-Ti tube of 50.8 mm (out diameter, D) 0.508 mm (wall thickness, t) 101.6 mm (bending radius, R) as a representative component, the bending behaviors of LDTW CP-Ti tube under different processing conditions were investigated. With experiments and analytical analysis, the bending characteristics of the CP-Ti tube were identified. Then, based on the orthogonal experimental design, a series of three-dimensional FE models of RDB for LDTW CP-Ti tube were established, and the effects of processing parameters on the bending behaviors were numerically investigated. The results show that: 1) Wrinkling is the primary behavior for the LDTW CP-Ti tube in RDB, and the larger the difference between the maximum wall thickening and the maximum wall thinning degrees of the bent tube, the larger the wrinkling tendency; 2) The bending behaviors of the LDTW CP-Ti tube are very sensitive to the some processing parameters, and the wrinkling is significantly affected by the mandrel shank diameter, while the wall thinning is remarkably affected by the clearance between wiper die and tube, mandrel shank diameter. The qualified bent tube with the wall thinning of 11.43%, the cross-section distortion of 2.69% and the wrinkling height less than 2% is then obtained. Key words: CP-Ti tube; large diameter; rotary draw bending; bending behaviors; processing parameters 1 Introduction Due to high specific strength, excellent corrosion resistance, good formability and weldability, the commercial pure titanium (CP-Ti) has attracted more and more applications to manufacture lightweight components in various industries such as aviation and aerospace [1]. To satisfy the flight performance requirements of aeroplanes, CP-Ti bent tubes have been widely used as hydraulic and fuel piping [2]. Among many bending methods, the rotary-draw-bending (RDB) is the most universal technology for the thin-walled tube bending [3] (as shown in Fig. 1) because of its versatility, relatively higher forming quality and efficiency. However, the RDB of large diameter thin-walled (LDTW) tube is a tri-nonlinear plastic forming process with multi-factor interactive effects, which needs precise coordination of various dies and strictly controlling of processing parameters. When the processing parameters are selected unsuitably, it will easily lead to excessive wall thinning, wrinkling and severe cross-sectional distortion, especially for the tube with large size factor (ratio of out diameter to wall thickness, D/t). Moreover, compared with the Al-alloy tubes, the undesired defects of CP-Ti tube may more easily occur during the RDB due to the larger flow stress and lower hardening exponent, and the bending behaviors may be much more sensitive to the processing parameters. Thus, to improve the forming quality, it is of great importance to study the bending behaviors of the LDTW CP-Ti tube in RDB and obtain the effects of the processing parameters. Up to now, many efforts have been done on the effects of processing parameters on the bending behaviors in thin-walled tube bending. WU et al [4] experimentally studied the influences of temperature, bending velocity and original grain size on the bendability of wrought magnesium alloy AM30 tubes Foundation item: Project (50905144) supported by the National Natural Science Foundation of China; Project (09-10) supported by State Key Laboratory of Materials Processing and Die & Mould Technology, China; Project (JC201028) supported by NPU Foundation for Fundamental Research, China; Project supported by the Fund of the State Key Laboratory of Solidification Processing in NWPU; Project (2011JQ6004) supported by the Natural Science Basic Research Plan in Shaanxi Province and the 111 Project (B08040), China Corresponding author: He YANG; Tel: +86-29-88495632; Fax: +86-29-88495632; E-mail: yanghe@nwpu.edu.cn

402 Zhi-yong ZHANG, et al/progress in Natural Science: Materials International 21(2011) 401 412 Fig. 1 Schematic diagram of RDB and typical stress/strain distributions of tube materials with the size factor of 10.5. YANG et al [5] studied the influences of bending radius and wall thickness on the prebending of the steel tube with the size factor of 24 during hydroforming. LI et al [6] numerically investigated the influences of mandrel parameters (mandrel diameter, mandrel extension, number of balls, thickness of balls, etc) on stress distribution and wrinkling of steel and Al-alloy tube with the size factor of 38. LI et al [7] experimentally studied the effects of processing parameters on the wall thinning and the cross section distortion for LDTW Al-alloy tubes with the D/t of 47 in RDB. YANG et al [8] obtained the effects of clearances and friction between different dies and tube on the maximum wrinkling factor of Al-alloy thin-walled tube bending with the size factors from 38 to 100. Most of these studies were carried out on the Al-alloy and steel tubes, and most of the size factors of the tubes were under 100. The results identified the key parameters affecting the bending behaviors of thin-walled tubes, including bending angle (θ), push assistant level (v p =V p /V, where V p is the pushing speed of the pressure die, and V is the speed of the bend die), mandrel shank diameter (d), mandrel extension length (e), number of balls (N), ball diameter (d 0 ), thickness of ball (k), clearance between pressure die and tube (C p ), clearance between wiper die and tube (C w ). These results provide useful knowledge to improve the bending quality of the LDTW CP-Ti tube. However, it is time consuming to test all the variables by single factor analysis method. By the range analysis and analysis of variance(ava), the orthogonal experimental design [9] has been one efficient way to qualitatively analyze the correlations among the relevant factors at different levels. Due to the multi-die constraints and multi-defect inherent to the LDTW CP-Ti tube during RDB, it is difficult to precisely analyze this process only analytically or experimentally. Based on the orthogonal experimental design, finite element (FE) simulation is an effective way to deal with the complex forming problems combined with the experimental and analytical methods [10]. Consequently, the tube specification is 50.8 mm (outer diameter, D) 0.508 mm (wall thickness, t) CP-Ti tube with large size factor D/t of 100 and bending radius of 2D. Based on the experiments and analytical analysis, the bending behaviors of LDTW CP-Ti tube are described and the corresponding objective functions are identified in terms of primary defects. Then, using the orthogonal experiment method, a series of threedimensional (3D) FE modeling and simulations are performed under the ABAQUS/Explicit platform. Finally, the effects of the processing parameters on the bending behaviors of LDTW CP-Ti tube are investigated and the significance sequences are lined via range analysis and AVA technique. 2 Experimental and analytical description of bending characteristics 2.1 Bending characteristics by experiments To understand the bending behaviors of the LDTW CP-Ti tube in RDB, the specification of 50.8 mm (D) 0.508 mm (t) CP-Ti tube with the D/t of 100 is selected for the physical experiments with the W27YPC-63 hydraulic tube bender. The bending dies are designed based on Ref. [2] and the mandrel parameters are chosen according to Ref. [6]. The experimental conditions are listed in Table 1. The experiments show that the wrinkling is the 402

Zhi-yong ZHANG, et al/progress in Natural Science: Materials International 21(2011) 401 412 403 Table 1 Experimental conditions for CP-Ti tube in RDB Forming parameter Value Diameter of mandrel/mm 49.64 Ball diameter/mm 49.54 Mandrel extension length/mm 0-6 Number of balls 3 10 Push assistant speed/(mm s 1 ) 3.5 Bending speed/(( ) s 1 ) 2.0 Final bending angle/rad π/2 Bending radius/mm 101.6 Pressure die length/mm 450 Wiper die length/mm 203 Mandrel shank length/mm 203 Pitch of balls/mm 2.3 Clamp die length/mm 203 Ball thickness/mm 12 primary behavior of the CP-Ti tube, which easily occurs near the clamp die, as shown in Fig. 2(a). One reason is that the large size factor leads to small anti-wrinkling ability and large wrinkling possibility. Another reason is the specific material properties of the CP-Ti tube at room temperature. The CP-Ti tubes exhibit remarkable anisotropy and are hard to deform due to the hexagonal close-packed crystal structure at room temperature [11], which makes the material flow hard during the bending process, and it will easily lead to the instability at the inner side of the tube, viz. wrinkling. Besides, because of the large flow stress of the CP-Ti tube, the wrinkling is difficult to be smoothed by the flexible balls. Since the waves of wrinkling are small, they may insert into the space between the first flexible ball and the mandrel shank or that of the flexible balls. When it happens, the joint of flexible balls may crack (as shown Fig. 3). On the other hand, due to the large size factor D/t of 100, the tube cross section will be distorted when clamped, which may lead to the crack of the joint when too many balls are used (about more than 6). The crack of the joint will lead to serious defect such as severe cross-section distortion (Fig. 2(b)). This phenomenon is rarely found in Al-alloy tube bending [6, 7]. So, compared with the Al-alloy and steel tubes, the CP-Ti tube bending needs even more strict cooperations of the processing parameters to avoid the wrinkling. Consequently, the preventions of wrinkling and crack of the joint are the main problems that should be solved in RDB of the LDTW CP-Ti tubes. When the number of flexible balls is less than 6, the crack of the joint due to the large number of balls can be avoided. On Fig. 2 Wrinkling and serious cross-section distortion: (a) Wrinkling; (b) Serious cross-section distortion due to crack of joint Fig. 3 Crack of joint this basis, the cross-section distortion can be reduced efficiently by increasing the number of balls [6]. Thus, the wrinkling is the key defect to be prevented. Due to the interactive effects of wall thinning and wrinkling, it is hard to make both of them in control. To achieve the stable bending of the LDTW CP-Ti tubes, it is necessary to study the effects of the processing parameters on wall thinning and wrinkling. 2.2 Analytical description and target functions In order to quantitively analyze the effects of the processing parameters on the wall thinning and wrinkling, the target functions of wall thinning and wrinkling are deduced according to the geometrical features in RDB 403

404 Zhi-yong ZHANG, et al/progress in Natural Science: Materials International 21(2011) 401 412 and deformation plasticity theory. To simplify the analysis, some assumptions are proposed as follows: 1) The proportional loading and uniform distributions of stress/strain are assumed. The tangent stress σ φ, hoop stress σ θ and normal stress σ t are regarded as the principal stresses, and corresponding strains are the principal strains. 2) As shown in Fig. 1, the normal stress σ t can be neglected compared with σ φ and σ θ. For the large D/t of 100, the bending process can be regarded as plane stress problem. 3) The elastic deformation is regarded as zero, and the constitutive relationships of tube materials accord with the Hill anisotropic yield model and Swift strain n hardening function σ = K ( ε + b), where σ and ε refer to the equivalent stress and strain, K is the strength coefficient, and b is the material constant. As shown in Fig. 4, according to the geometrical features of bending tube, the strain field can be approximately calculated by [12] (1 + rr ) σ 1+ r σϕ =..( εϕ + εθ) 1+ 2r ε r (1 + rr ) σ 1+ r σθ =..( εϕ + εθ) 1+ 2r ε r (2) where r refers to normal anisotropic exponent, and the equivalent strain ε can be calculated by Hill anisotropic yield model: 1+ r 2 2 2r ε = εϕ + εθ + εϕεθ (3) 1 + 2r 1 + r By substituting Eq. (3) into Eq. (2), the stress state can be obtained. When the point is at the outside of the tube, ε φ 0, Eq. (2) gives 1+ r 1 σϕ = LK( )( Mεϕ + b) r 4 1+ r n σθ = LK(1 )( Mεϕ + b) 4r n (4) (1 + r) r 1 17 where + r r L = > 0 and 1+ 2r 1 2 16 2(1 ) + r + r 1 εϕ = ln[ R /( Rd E0 )] εθ = 1/4εϕ εt = 3/4εϕ where R = Rd r cosθ, θ [0, π], r is the average radius of the tube cross-section, θ [0,π / 2] refers to inside tube and θ [π / 2, π] refers to outside tube. The neutral layer (NL)/neutral axis (NA) shift value E 0 = 2 0.84r / R d [13], and the hoop strain is obtained by the experiments [14]. Then, on the basis of the deformation plasticity theory, the stress state can be derived as (1) 1+ r 17 r M = > 0. 1+ 2r 16 2(1 + r) When the point is at the inside of the tube, ε φ 0, the equation becomes 1+ r 1 σϕ = LK( )( Mεϕ b) r 4 1+ r n σθ = LK(1 )( Mεϕ b) 4r n (5) Thus, the key deformation behaviors in terms of wall thinning and wrinkling instability can be discussed analytically. Fig. 4 Schematic diagrams of longitudinal (a) and transverse (b) cross-sections of bent tube geometry [10] 404

Zhi-yong ZHANG, et al/progress in Natural Science: Materials International 21(2011) 401 412 405 1) As shown in Fig. 4, the wall thickness changing degree of tube, t, can be calculated by the normal strain: t t Δt = = expε t 1 (6) t where t is initial tube wall thickness, and t' the wall thickness after bending. The wall thinning degree of tube is expressed as t t Δ t = = ( 1 expε t ) (7) t 2) As described in Ref. [15], the difference between the maximum tangent compressive stress and the corresponding tangent tensile stress I w is proposed to represent the wrinkling possibility from the stress point of view: I = σ σ (8) w max,c max,t where σ max,c is the maximum tangent compressive stress at the inner side of tube; and σ max,t is the maximum tangent compressive stress at outer side of tube. The larger I w means the larger tube wrinkling tendency. For stable forming process, I w follows the function: I I = σ σ (9) w ϕ,c ϕ,t where I is the analytical difference between the maximum tangent compressive stress and the corresponding tangent tensile stress, while, σ φ,c is the analytical tangent compressive stress at the extrados and σ φ,t is the analytical tangent tensile stress at the intrados obtained by Eqs. (4) and (5). It can be seen that both σ φ,c and σ φ,t are the functions of the tangent strain ε φ, and they have the same changing tendency with ε φ, so, it is reasonable to use the difference between the maximum tangent compressive strain and the corresponding tangent tensile strain to represent the wrinkling possibility. By Eqs. (1) and (6), the maximum tangent compressive strain and the corresponding tangent tensile strain can be replaced by the maximum wall thickening degree t max,c and the maximum wall thinning degree t max,t, so Eq. (9) can be changed to I = Δt Δt I = Δt ϕ Δt (10) w max,c max,t,c ϕ,t where I w is the difference between the maximum wall thickening degree and the maximum wall thinning degree of the bent tube, while I is the analytical one, and Δ t ϕ, c is the analytical maximum wall thickening degree and Δ t ϕ, t is the analytical maximum wall thinning degree. When I w I, it means that the tube has the wrinkling tendency or the wrinkling has occurred, and the larger the value of I, the larger the wrinkling tendency or the more serious the wrinkling. 2.3 Evaluation of analytical description To examine the precision of analytical results, in terms of maximum wall thickening degree Δ t max, c, maximum wall thinning degree Δ t max, t and I w, the analytical results are compared with experimental ones. The LDTW CP-Ti tube of 50.8 mm (D) 0.508 mm (t) is conducted. Tables 1 and 2 give the forming conditions and the material properties, and the number of balls is 4, while the mandrel extension length is 3 mm. The Δ t ϕ, c, Δ t ϕ,t and I are 18.96%, 18.57% and 0.39%, respectively, and Fig. 5 shows the comparison results. It can be seen that Δ t ϕ, c is near the experimental Δ t max, c at the stable forming part, and the experimental Δ t max, c is much larger near the initial bending section, where the minor wrinkling happens. But the Δ t ϕ, t is much larger than the experimental results Δ t max, t. The main reasons of the discrepancy are: 1) The analytical analysis does not consider the constraints of the dies, the pushing effect of the pressure die, etc, and they may have positive effects to reduce the wall thinning; 2) The cross-section of the bent tube has a little distortion, and it is unable to avoid in the physical experiments, while cross-section is kept invariant in the analytical analysis. Because of the error of Δ t ϕ, t, the I is smaller than the experimental I w. But the changing trend of I w can reflect the wrinkling tendency, which is helpful to target the wrinkling tendency. Thus, I w can be used as the target value to express the wrinkling tendency for further analysis. 3 Orthogonal experimental design method and FE modeling of RDB 3.1 Orthogonal experimental design The orthogonal experimental design is an efficient and convenient method to identify the effects of processing parameters on wall thinning and wrinkling. Table 2 Mechanical properties of CP-Ti tube Material parameter Value Elastic modulus, E/GPa 90.325 Fracture elongation, δ/% 30.014 Initial yield stress, σ 0.2 /MPa 351.428 Ultimate strength, σ b /MPa 492.032 Strength coefficient, K/MPa 756.038 Hardening exponent, n 0.149 Material constant, b 0.001 93 Normal anisotropic exponent, r 1.50 Density, ρ/(kg m 3 ) 4 540 405

406 Zhi-yong ZHANG, et al/progress in Natural Science: Materials International 21(2011) 401 412 Table 3 Factors and their levels in orthogonal experimental design Code Factor Level 1 Level 2 Level 3 Level 4 A θ/rad π/3 π/2 2π/3 π B v p /% 90 100 110 120 C d/mm 49.60 49.64 49.68 49.72 D E/mm 0 2 4 6 E N 3 4 5 6 F d 0 /mm 49.50 49.54 49.58 49.62 G K/mm 10 12 14 16 H C p /mm 0 0.02 0.04 0.06 I C w /mm 0 0.02 0.04 0.06 wall thinning. Based on the method described by YANG et al [8] and the orthogonal experimental design plan, a series of 3D-FE models are established under FE platform ABAQUS/Explicit to simulate the RDB of LDTW CP-Ti tube with multiple dies (as shown in Fig. 6). The basic geometry conditions are the same as the experimental ones (as listed in Table 1). The element size of the tube is 1.5 mm 1.5 mm to obtain the tradeoff between numerical accuracy and stability, and to neglect the dynamic effects of explicit FE formulation. The mass scaling factor is set as 6 000. Fig. 5 Comparisons of analytical and experimental results: (a) Δ tmax ; (b) I w The wall thinning degree Δ t max, t and the wrinkling tendency I w are regarded as the target values to mark wall thinning and wrinkling, respectively. They are described as follows: t te Δ tmax,t = (11) t where t is the initial tube wall thickness, t e is the wall thickness at the extrados of the bent tube. I = Δt Δt (12) w max,c max,t Nine relevant processing parameters with four levels for each are investigated, as listed in Table 3. Based on the influence factors and corresponding levels, the standard orthogonal array L 32 (4 9 ) is arranged. 3.2 FE modeling of RDB with LDTW CP-Ti tube The FE simulation can quantitatively describe the influences of processing parameters on wrinkling and Fig. 6 3D FE model for RDB of large diameter thin-walled CP-Ti tube with multiple dies The material studied in this work is CP-Ti welded tube (SAE AMS4941E, 2008) with the outside diameter of 50.8 mm and wall thickness of 0.508 mm. The weld seam is placed at the neutral position of tube during the RDB. Since the plastic strain at this place is minor in bending, the effects of weld and heat affect zone on the deformation behaviors of the CP-Ti tube can be neglected. The properties of the parent metal are presented in Table 2. The Coulomb friction model is used to describe the tangent behavior between different contact interfaces. Table 4 gives the friction coefficients about different contact pairs. The double precision is used for bending calculation to ensure the accuracy. 406

Zhi-yong ZHANG, et al/progress in Natural Science: Materials International 21(2011) 401 412 407 Table 4 Coefficients on various contact interfaces Contact interface Friction coefficient Tube-wiper die 0.10 Tube-pressure die 0.30 Tube-clamp die 0.6(rough) Tube-bend die 0.25 Tube-mandrel 0.10 Tube-balls 0.10 3.3 Evaluation of 3D FE model To validate the FE model, the experiments are conducted by W27YPC-63 hydraulic tube bender. The experimental conditions are the same as the Section 2.3. Figure 7 shows the comparison results in terms of the maximum wall thinning degree t max,t and the maximum wall thickening degree t max,c. It can be seen that, t max,c of the FE simulation fluctuates near the initial bending section where the minor wrinkling occurs, which is in good accordance with the experimental result. The simulated t max,t and t max,c are 12.7% and 24.8%, while those of experimental results are 11.21% and 26.73%, and the relative error is 13.29% and 7.22%, respectively. Thus, the results show that the FE model is reliable. Fig. 7 Comparisons of FE simulated and experimental results: (a) Maximum wall thinning degree t max,t ; (b) Maximum wall thickening degree t max,c 4 Results and discussion 4.1 Results of FE simulation After the simulation of the RDB arrayed by orthogonal experimental design, the simulated results are obtained by postprocessing, which are listed in Table 5. It is found that the bending behaviors of the bent tubes present four typical patterns, viz. wrinkling, wrinkling tendency, severe cross-sectional distortion and good forming (as shown in Fig. 8). Table 5 Results of FE simulation No. t max,t I w Bending status 1 14.55 3.29 Severe distortion 2 12.17 10.21 Wrinkling tendency 3 12.36 5.30 Good 4 14.37 2.82 Good 5 17.40 5.65 Good 6 12.87 12.18 Wrinkling tendency 7 11.63 1.24 Good 8 27.54 21.11 Severe distortion 9 16.92 0.81 Good 10 13.06 10.90 Wrinkling tendency 11 18.69 8.90 Good 12 15.82 5.79 Severe distortion 13 17.89 2.78 Wrinkling tendency 14 14.39 17.26 Wrinkling 15 13.90 2.52 Good 16 13.43 2.00 Good 17 15.92 0.14 Wrinkling tendency 18 12.87 4.03 Good 19 12.70 8.13 Wrinkling tendency 20 12.63 18.83 Wrinkling 21 17.22 8.47 Good 22 17.15 5.83 Good 23 13.75 17.18 Wrinkling 24 14.47 24.06 Wrinkling 25 17.93 9.23 Wrinkling 26 25.43 19.32 Severe distortion 27 18.40 36.8 Wrinkling 28 11.72 2.66 Good 29 17.42 0.17 Good 30 12.45 0.50 Good 31 14.86 26.14 Wrinkling 32 14.03 10.5 Wrinkling tendency 407

408 Zhi-yong ZHANG, et al/progress in Natural Science: Materials International 21(2011) 401 412 Fig. 8 Typical patterns of simulated results: (a) Wrinkling; (b) Wrinkling tendency; (c) Severe distortion cross-sectional caused by over thinning; (d) Good forming As listed in Table 5, t max,t ranges from 11.63% to 27.54%, and I w ranges from 19.32% to 26.14%. The wrinkling is the most possible defect during the LDTW CP-Ti tube RDB. When I w is near the analytical result I (0.39%), the forming quality is good. But when I w is too large (>8 or so) or too small (< 19 or so), the wrinkling and the severe cross-sectional distortion may occur, respectively. When N<3, there is a larger chance that the severe cross-sectional distortion occurs, shown as No. 1 and 12. The over thinning may occur while d is 49.72 mm, which may cause severe cross-sectional distortion, such as No. 8 and 26. When v p <100%, because of the negative effect of the pressure die, the wrinkling is much easier to occur, shown as No. 13 and 17. So, for the stable bending process, N should be more than 3 and v p should not be smaller than 100%. These results will be helpful to the experimental research of LDTW CP-Ti tube in RDB. To extract the significance sequence of the factors, the simulated results are analyzed by range analysis. Then, the analysis of variance (AVA) technique is applied to determine the statistical significance of the investigated processing parameters. 4.2 Range analysis The range value (R) for each factor is the difference between the maximal and the minimal value of the specific levels. It can be used to clarify the significance level of each factor on the target values. During the range analysis, the average of eight values for each level K is calculated, and then R is gotten, as listed in Table 6. Table 6 Range value for t max,t and I w Factor R t R I θ 3.80 3.95 v p 2.37 9.62 d 4.18 24.89 e 1.60 1.95 N 2.43 3.08 d 0 1.05 5.20 k 3.67 6.63 C p 1.52 7.75 C w 4.28 10.31 R t is range value for t max,t, and R I is range value for I w. 408

Zhi-yong ZHANG, et al/progress in Natural Science: Materials International 21(2011) 401 412 409 Based on the results of range analysis, the significance sequence of all the investigated factors from the largest to the smallest is lined. For t max,t, the order is as follows: C w, d, θ, k, N, v p, e, C p, d 0. And for I w, the order is: d, C w, v p, C p, k, d 0, θ, N, e. Compared with the RDB of the Al-alloy tubes, the order is similar for I w ' [6, 8], but it is different for t max,t [7]. To further understand the effects of each influencing factor, the influencing tendencies of each investigated factor on wall thinning and wrinkling are clarified, as shown in Fig. 9. From Fig. 9, it can be seen that: values of t max,t are all smaller and most values of I w are larger than the analytical ones, and the reason may be the limitation of the analytical description analyzed before. The effects of the processing parameters on t max,t and I w are complicated and they have interactive effects with each other. t max,t and I w are both sensitive to d, and t max,t increases and I w decreases sharply with the increasing of d. With the increase of θ and k, both t max,t and I w have a trend to increase. With C w and v p increasing, there is the tendency that t max,t decreases and I w increases. These results are similar to those of the Al-alloy tubes [6, 8]. However, the effect of mandrel shank diameter d on the bending behaviors of the LDTW CP-Ti tube is much more remarkable, meanwhile, the mandrel extension length e and the ball diameter d 0 have little effects on the wall thinning, which is different from the results of Al-alloy tube [7]. In order to better understand the influence mechanism of d, e and d 0, the stress distribution under different bending conditions has been revealed using the FE simulation. Figure 10 shows the history curve of the maximum tangent stress with different bending conditions. From Figs. 10(a) and (b), it can be seen that the maximum tangent compressive stress rises by 40 MPa with the small decrease of d of 0.04mm. Due to the small hardening exponent of the CP-Ti tube, the larger compressive stress should cause larger increasing of the inner wall thickness, which implies larger possibility of the wrinkling. From Figs. 10(a), (c) and (d), it is found that the distribution of the maximum tangent compressive stress and the tensile stress are almost unchanged with the changes of e and d 0. So, due to the large flow stress of the CP-Ti tube, the increased backward dragging force caused by increasing e and d 0 slightly affects the strain during the RDB, and Fig. 9 Influencing tendency of factors on t max,t ((a) and (b)) and I w ((c) and (d)) 409

410 Zhi-yong ZHANG, et al/progress in Natural Science: Materials International 21(2011) 401 412 Fig. 10 History curves of maximum tangent stress with different bending conditions: (a) Basic condition (d=49.68 mm, e=2 mm, d 0 =49.50 mm); (b) d=49.64 mm, e=2 mm, d 0 =49.50 mm; (c) d=49.68 mm, e=6mm, d 0 =49.50 mm; (d) d=49.68 mm, e=2 mm, d 0 =49.62 mm the effects of e and d 0 are limited to the wall thinning and wrinkling. Considering both the effects of the processing parameters on the t max,t and I w, the optimal combination is A 2 B 2 C 3 D 2 E 4 F 1 G 1 H 1 I 2. In order to verify the optimum conditions, corresponding confirmation experiments are conducted, and the result shows that t max,t is 12.71% and I w is 0.10%. At these optimum conditions, t max,t and I w are both at a low level compared with the results obtained from orthogonal array, meanwhile, I w is very close to the analytical result I. 4.3 Analysis of variance (AVA) In the orthogonal experimental design, AVA is used to analyze the results of the orthogonal experiments and to determine how much variation that each factor has contributed to the target values. The AVA on t max,t and I w is carried out and the results are given in Tables 7 and 8. The calculated ratio F>F 0.01 (3,4) indicates the factor with high significance on the results, F 0.05 (3,4)< F<F 0.01 (3,4) indicates significant, F 0.1 (3,4)<F<F 0.05 (3,4) indicates lightly significant, and F<F 0.1 (3,4) indicates statistically nonsignificant. From the calculated results, it can be seen that C w, d and θ have significant influences on t max,t, and it is statistically significant at 95% confidence limit. k shows lightly significant influences and the other factors have statistical nonsignificance on t max,t. The wrinkling( I w ) is significantly affected by d, while the other factors have little effects on the wrinkling of LDTW CP-Ti tube in RDB. Thus, d and C w are the most important processing parameters that affect the bending behaviors, and they should be considered preferentially during the RDB of LDTW CP-Ti tube. Following the above analysis, by mainly adjusting mandrel shank diameter (d), the clearance between wiper die and tube (C w ) and the push assistant level (v p ), the qualified bent CP-Ti tubes are finally achieved. The maximum wall thinning of the qualified part is 11.43%, while the maximum cross-section distortion is 2.69%, and the wrinkling height is less than 2% (as shown in Fig. 11). 410

Zhi-yong ZHANG, et al/progress in Natural Science: Materials International 21(2011) 401 412 411 Table 7 AVA for t max,t Code Source Sum of squares Degree of freedom Mean square F value Critical value Significance A θ 70.160 3 23.387 7.26 F 0.01 (3,4)= * * B v p 24.880 3 8.293 2.57 16.7 C d 81.618 3 27.206 8.44 * * D e 11.642 3 3.881 1.20 E N 37.351 3 12.450 3.86 F 0.05 (3,4)= 6.59 F d 0 4.522 3 1.507 0.468 G k 58.796 3 19.599 6.08 * H C p 11.315 3 3.772 0.23 F 0.1 (3,4)= 4.19 I C w 84.976 3 28.325 8.79 * * Error 16.114 4 3.223 Total 401.372 9 31 ** Significant; * Lightly significant; Nonsignificant. Table 8 AVA for I w Code Source Sum of squares Degree of freedom Mean square F value Critical value Significance A θ 81.176 3 27.059 0.327 B v p 401.461 3 133.820 1.615 F 0.01 (3,4)= 16.7 C d 2663.554 3 887.851 10.717 * * D e 23.453 3 7.818 0.094 E N 42.684 3 14.228 0.172 F 0.05 (3,4)= 6.59 F d 0 178.258 3 59.419 0.717 G k 183.354 3 61.118 0.738 H C p 242.734 3 80.911 0.244 F 0.1 (3,4)= 4.19 I C w 524.211 3 174.737 2.109 Error 331.371 4 82.843 Total 3 741.842 31 ** Significant; Nonsignificant. 5 Conclusions Fig. 11 Qualified bent LDTW CP-Ti tubes achieved by experiments 1) Due to the large size factors and specific material properties, the bending quality of the LDTW CP-Ti tube is very sensitive to the processing parameters and the tube is much more difficult to bend. The wrinkling is the primary behavior of the CP-Ti tube during the physical experiments. The crack of the joint easily occurs when the wrinkling has happened or more flexible balls are used, and it is rarely found in Al-alloy tube bending process. 2) It is found that the wrinkling tendency can be represented by the difference between the maximum wall thickening degree and the maximum wall thinning degree of the bent tube, and it is very sensitive to the 411

412 Zhi-yong ZHANG, et al/progress in Natural Science: Materials International 21(2011) 401 412 mandrel shank diameter d during the RDB of the LDTW CP-Ti tube. The significance sequence of the investigated processing parameters from the largest to the smallest for wrinkling is: d, C w, v p, C p, k, d 0, θ, N, e. 3) In RDB process of the LDTW CP-Ti tube, the wall thinning is significantly affected by the mandrel shank diameter d and the clearance between wiper die and tube C w. But, the effects of the mandrel extension length e and the ball diameter d 0 are limited, which is different with the Al-alloy tube. The significance sequence of the investigated factors from the largest to the smallest for wall thinning is: C w, d, θ, k, N, v p, e, C p, d 0. 4) By adjusting the processing parameters based on the above analysis, qualified bending of LDTW CP-Ti tube (50.8 mm (D) 0.508 mm (t)) with the bending radius 101.6mm is obtained with the wall thinning of 11.43%, the cross-section distortion of 2.69% and the wrinkling height less than 2%. The results may be helpful to better understand the influencing role of processing parameters in RDB and improve the forming quality of LDTW CP-Ti bent tubes. References [1] BAPTISTA C A R P, BARBOZA M J R, ADIB A M L, et al. High temperature cyclic pressurization of titanium ducts for use in aircraft pneumatic systems [J]. Mater & Des, 2009, 30(5): 1503 1510. [2] SCHWEIHOLD D J. Titanium tube bending for aerospace [J]. Aircr Eng Aerosp Tec, 1989: 61(12): 8 29. [3] YANG H, SUN Z C, LIN Y, et al. Advanced plastic processing technology and research progress on tube forming [J]. J Plast Eng, 2001, 8(2): 86 88. (in Chinese). [4] WU W, ZHANG P, ZENG X, et al. Bendability of the wrought magnesium alloy AM30tubes using a rotary draw bender [J]. Mater Sci Eng A, 2008, 48(1 2): 596 601. [5] YANG J,JEON B, OH S. The tube bending technology of a hydroforming processes for an automotive part [J]. J Mater Process Technol 2001, 111: 175 181. [6] LI H, YANG H, ZHAN M, et al. Role of mandrel in NC precision bending process of thin-walled tube [J]. Int J Mach Tools Manuf, 2007, 47(7 8): 1164 1175. [7] LI C, YANG H, ZHAN M, et al. Effects of process parameters on numerical control bending process for large diameter thin-walled aluminum alloy tubes [J]. Trans Nonferrous Met Soc China, 2009, 19: 668 673. [8] YANG H, YAN J, ZHAN M, et al. 3D numerical study on wrinkling characteristics in NC bending of aluminum alloy thin-walled tubes with large diameters under multi-die constraints [J]. Comput Mater Sci, 2009, 45: 1052 1067. [9] LIU W Q. Design of experiments [M]. Beijing: Tsinghua University Press, 2005: 64 76. (in Chinese) [10] ZHAN M, YANG H, JIANG Z Q, et al. 3D FEM analysis of forming parameters on cone spinning based on orthogonal experimental design method [C]// The 9th International Conference on Technology of Plasticity. Gyeongju, Korea, 2008: 374 388. [11] LUTJERING G, WILLIAMS J C. Titanium (2nd edition) [M]. Heidelberg: Springer, 2007: 15 22. [12] LI H, YANG H, YAN J, et al. Numerical study on deformation behaviors of thin-walled tube NC bending with large diameter and small bending radius [J]. Comput Mater Sci, 2009, 45: 921 934. [13] TANG N C. Plastic-deformation analysis in tube bending [J]. Int J Pressure Piping, 2000, 77: 751 759. [14] KOU Y L. Experimental research on NC bending with small bending radius of aluminum alloy thin-walled tube [D]. Xi an: Northwestern Polytechnical University, 2007. (in Chinese) [15] YANG H, LI H, ZHAN M. Friction role in bending behaviors of thin-walled tube in rotary-draw-bending under small bending radii [J]. J Mater Process Technol, 2010, 210: 2273 2284. 412