Nội dung

TẠP CHÍ KHOA HỌC CÔNG NGHỆ GIAO THÔNG VẬN TẢI, SỐ 30-11/2018 57 DEFORMATION OF THIN-WALLED CIRCULAR TUBE SUBJECTED TO IMPACT THEREE-POINT BENDING BY USING NUMERICAL SIMULATION NGHIÊN CỨU BIẾN DẠNG CỦA ỐNG TRÒN THÀNH MỎNG CHỊU TẢI VA ĐẬP UỐN BA ĐIỂM BẰNG MÔ PHỎNG SỐ Ly Hung Anh Ho Chi Minh City University of Technology - VNU-HCM lyhunganh@hcmut.edu.vn Abstract: Crashworthiness is one of the most important criteria in the design of piping systems, suspension pipes in particular or energy absorbers in general. The objective of this paper is to study the deformation of thin-walled tube subjected to impact three-point bending using numerical simulation. Results are agreed very well with theoretical and experimental results. Based on the finite element modeling, deformation of thin-walled tube is presented when the diameter and spacing of two supporters change. Limited ratio of diameter and thickness is clarified to prevent overall bending in designing circular tube. Keywords: Crashworthiness, circular tube, impact, three-point bending, simulation. Classification number: 2.4 Tóm tắt: An toàn khi va chạm luôn được xem là một trong những tiêu chí quan trọng của thiết kế các hệ thống ống dẫn, ống treo nói riêng hay những thiết bị hấp thụ năng lượng nói chung. Nội dung chính của bài báo này là nghiên cứu ứng xử biến dạng của ống tròn thành mỏng chịu tải va đập uốn ba điểm bằng phương pháp mô phỏng số. Kết quả mô phỏng bằng phương pháp phần tử hữu hạn đúng so với tính toán bằng lý thuyết và cả thực nghiệm. Biến dạng của ống tròn thành mỏng được trình bày với các giá trị khác nhau của đường kính và khoảng cách của hai gối đỡ. Tỷ số giới hạn giữa đường kính ống và bề dày được tìm thấy để tránh hiện tượng uốn toàn cục khi ống tròn chịu tải va đập ngang. Từ khóa: Va chạm, ống tròn, uốn ba điểm, mô phỏng. Chỉ số phân loại: 2.4 1. Introduction Crashworthiness is always one of the most important criteria in design. Crashworthiness is defined a deformation in controlled manners without failure of structure. In this point of view, behaviour of thin - walled circular tube subjected to bending impact load by numerical method using LS - DYNA is presented in this paper. Crushing force and displacement at impact position are considered. In this paper, behaviour of circular tube is analyzed based on three-point bending theory that is developed by Wierzbicky [1]. The energy from impactor is absorbed entirely by the formation of hinge lines that cause deformation on the tube. The deformation of the tube is divided into distinct zones of compression and tension. Energy equilibrium equation is applied to obtain solutions for specific parameters such as mean crushing force, instantaneo us crushing force as well as bending moment. (1) Eext = Eint and E ext = E int E ext = P × δ and E int = ∑ E l With: Eext : Energy of impactor; Eint : Energy absorbed of tube; El : Energy absorbed by folds; P : Instantaneous force; δ : Displacement. Fig.1. Theoretical deformation of circular tube [2],[3]. (2) 58 Journal of Transportation Science and Technology, Vol 30, Nov 2018 Integrate equation (7) over α to find the During impact, the total energy from the impactor is absorbed by the displacement of value of mean force P m the hinge lines as shown in figure 1 and H  Pm M 0 ( 0.5π R + 476 H +0.61 H 2 + 2.47 R 2  (8) displacement of hinge lines are listed in table = R  1. Energy absorbed of tube, Eint , is calculated And then the bending moment M(θ) by sum of the energy generated by the   54  movement of hinge lines.  R  1  = M (θ ) M 0 1.76  1    + 3.15 Rate of energy dissipation for each hinge θ     4  t   line: 1  1  (3) E1 = M 0π Rα (9) + 2 R + 1.36t ( Rt ) 4    θ  (4) E 2 = M 0π Rθ 2. Finite element model 2 2 The finite element model shown in figure (π R ) + H 2 γ = E 3 M 0 (5) 2 is developed based on the three-point 4 bending model. The dimension of the circular (6) tube is the same as the sample in Mamalis [4] E int = E1 + E 2 + 4 E 3 to provide a comparison with empirical result. In which: Table 1. Displacement of hinge lines. - Outer diameter D = 30 mm; Hinge Rotational Lengths - Thickness t = 1.4 mmp; Line rates - Length L = 200 mm; Line 1 α πR - Impactor and supporter radius are 5 mm; Line 2 πR θ - Distance between 2 supporters L sup = 160 mm. 2 (π R ) 2 γ Line 3 Numerical simulation is performed with +H 4 constraint on 6 degrees of freedom for the 2 Instantaneous crushing force P (α ) , stoppers and 5 degrees of freedom for the impactor, which is only allowed to move mean crushing force Pm and bending moment vertically. Belytschko- Tsay four-node shell M (θ ) are obtained by substituting the above element are used with size of element to be 2.5 x 2.5 mm for whole specimen. Material equations into equation (1) and (2). properties (stainless steel 316CW) of the tubes π  πR P (α ) M0  = + are chosen to match with experiment in  2 H α 1.63 Mamalis [4] are presented in table 2. With: 2 2  π H  1  σ y : Yield stress; (7) +  + H2    4 1.63 R α      σ u : Ultimate stress; With: E: Young’s modulus; M 0 = σ 0t 2 / 4 : Yield moment per unit ν : Poison ration. length, σ 0 = 0.92σ u : Flow stress of tube’s material, H: Half length, α : Folding angle. Fig.2. Model of three-point bending. TẠP CHÍ KHOA HỌC CÔNG NGHỆ GIAO THÔNG VẬN TẢI, SỐ 30-11/2018 59 Table 2. Material properties of stainless steel 316CW [4], [5]. σ y (MPa) σ u (MPa) E (GPa) ν 420 950 207 0.3 3. Simulation results In bending problem, crushing force and displacement at impact loaded position on the tube need to be understand well. Results of numerical simulation are compared with experimental results in Mamalis [4] and analytical results in Yucheng Liu [2]. Fig. 4. Comparison of instantaneous force between simulation and experiment. Fig.3. Deformation in simulation and in experiment. Qualitative comparison of deformation in numerical simulation and experiment is presented in figure 3. Numerical simulation of instantaneous crushing force, which is shown in figure 4, is comparable fairy well with experimental results in Mamalis [4]. Several value of friction coefficients between impactor and tube, between tube and two supporters are plotted in figure 4 to clarify the effect of friction to numerical results. It is proved that this effect is not much if the displacement at impact loaded position on the tube rather small. Mean crushing force obtained from numerical simulation shown in figure 5 has good agreement with analytica l result. Therefore, it can be convinced that the numerical model in this study is reliable. According to simulation as well as experiment, deformation process of circular tube consists of three sequential phases. In phase 1, the tube is dented at contacting surface without bending the bottom surface. In consequence, the force rises up. In phase 2, stress increases and approaches to yield stress while denting goes on and the tube starts to bend. At this time, the force is almost saturated at the peak. After that, the tube is totally bent and the force drops rapidly corresponding to phase 3. Fig.5. Comparison of crushing force between simulation and analytics. 4. Deformation of cross-section at impacted position In industry, circular tube is widely used in piping systems. Hence, the tube section is an important factor need to be considered during impact. In this section, the pipe cross section is investigated by changing the pipe diameter parameters based on the standard dimens io ns in application. Fig.6. Location of calculation node. Finite element model has the same dimensions as the previous section. The impactor is applied mass of 40 kg and velocity of 10 m/s to get the initial kinetic energy. The thickness of 2 mm is maintained for all of circular tubes which have different diameters. The results are expressed by the displaceme nt 60 Journal of Transportation Science and Technology, Vol 30, Nov 2018 of the impactor, Node 1 (the node at impact position) and Node 2 (the node at the bottom of the tube as shown in figure 6. Some typical results of deformation behaviour with different values of D/t and distance between two supporters are shown in figure 7. When the diameter increases and the thickness is kept constant, the tube is less dented. The deformation process is divided into 3 phases as mentioned in part 3. The bending phase does not happen as the diameter increases because the energy of the load has been absorbed entirely by the first denting phase. At that time, the displacement is only due to denting behaviour without bending. When D/t is about 135, displacement starts to increase. That value is considered the limit in designing circular tube to prevent overall bending. Besides, L sup = 2000÷4500 mm are all calculated similarly to the case of L sup = 2000 mm and each of them has its own limited D/t value, which tents to rise following the increase of L sup . Fig.7. Displacement study with different D/t ratio and L sup . Figure 8 provides a good estimatio n capability for designing supporters for piping system in industry, which helps reduce time for preliminary design. Fig.8. Limited value of D/t with distance between two supporters. 5. Conclusion This paper presents bending behavior of thin-walled circular tube used widely in industry. Simulation results have good agreement with theoretical and empirica l results. Tubes with increasing diameter and constant thickness are less dented; however, this is only true for a limited D/t value. Given the distances of the two stoppers, the limited value of the D/t is identified. These results support more detailed understanding of three-point bending behavior and contribute practically for industry using circular piping system  Acknowledgement This work was supported in part by Ho Chi Minh City University of Technology – VNU-HCM under Grant T-KTGT-2017-58. References [1] Wierzbicky T., Recke L., Abramowicz T. G. W., Huang J., Stress profiles in thin-walled primatic TẠP CHÍ KHOA HỌC CÔNG NGHỆ GIAO THÔNG VẬN TẢI, SỐ 30-11/2018 columns subjected to crush loading-II bending, Computers & Structures, 51 (1994) 625-641. [2] Liu Y., Day M., Bending collapse of thin-walled circular tubes and computational application, International Journal of Crashworthiness, 46 (2008) 442-450. [3] Liu Y., Day M., Bending collapse of thin-walled beams with channel cross-section, International Journal of Crashworthiness, 11 (2006) 251-261. [4] Mamalis A. G., Manolakos D. E., Ioannidis M. B., 61 Kostazos P. K., Bending of cylindrical steel tubes: numerical modelling, International Journal of Crashworthiness, 11 (2006) 37-47. [5] Kamaya M., Stress–strain curve estimation procedures for stainless steels, Engineerin g Fracture Mechanics, 127 (2014) 194-21 Ngày nhận bài: 8/10/2018 Ngày chuyển phản biện: 12/10/2018 Ngày hoàn thành sửa bài: 1/11/2018 Ngày chấp nhận đăng: 9/11/2018