30.08.2013 Views

Finite Element Modeling of Crushing Behaviour of Thin Tubes with ...

Finite Element Modeling of Crushing Behaviour of Thin Tubes with ...

Finite Element Modeling of Crushing Behaviour of Thin Tubes with ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

13 th International Conference on<br />

AEROSPACE SCIENCES & AVIATION TECHNOLOGY,<br />

ASAT- 13, May 26 – 28, 2009, E-Mail: asat@mtc.edu.eg<br />

Military Technical College, Kobry Elkobbah, Cairo, Egypt<br />

Tel : +(202) 24025292 – 24036138, Fax: +(202) 22621908<br />

<strong>Finite</strong> <strong>Element</strong> <strong>Modeling</strong> <strong>of</strong> <strong>Crushing</strong> <strong>Behaviour</strong><br />

<strong>of</strong> <strong>Thin</strong> <strong>Tubes</strong> <strong>with</strong> Various Cross-Sections<br />

M. M. Younes *<br />

1/19<br />

Paper: ASAT-13-ST-34<br />

Abstract: Typical steel square plates <strong>of</strong> 300mm side length and 1mm thickness were used to<br />

simulate six tubes <strong>with</strong> different cross-sections. These cross-sections were selected symmetric<br />

from round-shapes as circle and ellipse and polygon-shapes as triangle, square, pentagon and<br />

hexagon. The effect <strong>of</strong> varying configurations <strong>of</strong> tube cross-section shape on the deformation<br />

response, collapse mode and energy absorption characteristics <strong>of</strong> tubes under quasi-static<br />

axial compression have been studied numerically. The commercial finite element package<br />

ABAQUS/Explicit was used in the present analysis. The axial load-displacement results<br />

accompanied by the fold formation <strong>of</strong> various tubes were investigated and compared <strong>with</strong><br />

published experimental work. Variation <strong>of</strong> the initial peak load and the mean crushing force<br />

<strong>with</strong> the tube side breadth and the fold depth were carefully examined. Based on the finite<br />

element simulation results, empirical formulas for the absorbed energy and the fold depth<br />

were developed as a function <strong>of</strong> the side breadth-to-thickness ratio <strong>of</strong> the tubes.<br />

Keywords: <strong>Finite</strong> element; Initial peak load; Mean crushing force; Fold formation; Axial<br />

loading; Absorbed energy; <strong>Tubes</strong><br />

Introduction<br />

Numerical methods are now extensively applied in engineering due to the advances in<br />

computing. Of all the numerical methods, the finite element method is the most popular<br />

convenient approach, because it is easy to implement for all kinds <strong>of</strong> boundary and loading<br />

conditions and it can be used for the analysis <strong>of</strong> large complex structures. It is worth<br />

mentioning that experimental measurements are considered a powerful data required to<br />

approve most <strong>of</strong> mathematical models. Paik et al. [1] performed a series <strong>of</strong> experiments on the<br />

quasi-static crushing <strong>of</strong> stiffened square tubes and brought forward an empirical formula for<br />

the mean crushing load. The experimental results showed that the mean crushing load<br />

changed <strong>with</strong> the variation <strong>of</strong> the stiffener directions. Zhang and Suzuki [2] analyzed the<br />

effects <strong>of</strong> longitudinal and transverse stiffeners on the quasi-static crushing <strong>of</strong> stiffened square<br />

tubes by using the non-linear finite element program LS-DYNA. The obtained numerical data<br />

was used in prediction the mean crushing load and the equivalent thickness <strong>of</strong> the stiffened<br />

tubes. Recent numerical tools have allowed that a range <strong>of</strong> quasi-static problems can be<br />

analyzed effectively. Gupta et al. [3] applied the finite element method to simulate the axial<br />

compression <strong>of</strong> the tubes <strong>of</strong> round cross-sections. An axi-symmetric model was used to obtain<br />

the deformed shapes and the crushing load at different stages <strong>of</strong> loading. Gupta and Gupta [4]<br />

have carried out experimental studies <strong>of</strong> axial compression <strong>of</strong> thin walled aluminum and mild<br />

steel tubes <strong>of</strong> varying diameter to thickness ratios. The effect <strong>of</strong> annealing and presence <strong>of</strong> a<br />

* Egyptian Armed Forces


8/19<br />

Paper: ASAT-13-ST-34<br />

hole at mid height <strong>of</strong> the tube on the energy absorption characteristics and collapse mode<br />

transitions was examined. Results revealed that a diamond deformation mode was generated<br />

in all the tube specimens. A great number <strong>of</strong> studies have been carried out on the axial<br />

crushing <strong>of</strong> thin-walled tubes. Wierzbicki and Abramowicz [5] developed a simplified model<br />

based on rigid plastic assumptions and obtained the mean crushing force for square tubes.<br />

Abramowicz and Jones [6,7] conducted a series <strong>of</strong> experiments, and modified the collapse<br />

models by taking account <strong>of</strong> strain-rate sensitivity and effective crushing distance. Grzebieta<br />

[8] proposed a method for determining the load history between a peak and a minimum during<br />

an oscillation <strong>of</strong> the load-compression curve <strong>of</strong> round tubes. Gupta and Velmurugan [9]<br />

studied experimentally the internal/external folding <strong>of</strong> round tubes. The folding parameter<br />

(ratio <strong>of</strong> the inside to the total fold length) was determined experimentally and employed in a<br />

proposed analysis. Wierzbicki et al. [10] studied the axi-symmetric mode <strong>of</strong> deformation <strong>of</strong><br />

round tubes by considering partly internal and partly external folding. It was shown that the<br />

load compression curve is dependent on the folding parameter, while the mean collapse load<br />

and the folding length are independent <strong>of</strong> it. Closed form solutions were obtained for the<br />

instantaneous and mean crushing force, the effective crush distance and the length <strong>of</strong> the local<br />

folding wave. Singace et al. [11-13] gave an analysis <strong>of</strong> the axi-symmetric and multi-lobe or<br />

diamond mode <strong>of</strong> deformation to determine the eccentricity factor and crushing load. It was<br />

concluded that measured values <strong>of</strong> the eccentricity factor and the critical folding angles<br />

obtained for tubes <strong>of</strong> different materials and geometric ratios are independent <strong>of</strong> the tubes<br />

material and geometric ratios. There is a considerable amount <strong>of</strong> published data on the<br />

response <strong>of</strong> composite tubes to axial crushing [14]. Many <strong>of</strong> these studies utilize circular<br />

cross-section tubular specimens to determine the energy absorption capability <strong>of</strong> the material.<br />

Farley [15] studied the effect <strong>of</strong> specimen geometry on the energy absorption capability <strong>of</strong><br />

composite materials. He found that, changes in section lay-up that lead to an increase in<br />

modulus lead to higher crush strengths and energy absorption. Mamalis et al. [16-18]<br />

analyzed the collapse behavior and deformation mechanism <strong>of</strong> thin-walled non-circular<br />

composite tube, thin-walled composite conical shell, and braced elliptical tube. They<br />

contributed <strong>with</strong> a valuable data <strong>of</strong> crashworthiness <strong>of</strong> composite structures.<br />

The aims <strong>of</strong> this study are to obtain numerical data on the crushing <strong>of</strong> different cross-section<br />

tubes, and to develop simple empirical expressions for predicting the absorbed energy and the<br />

fold formation <strong>of</strong> steel tubes. A series <strong>of</strong> finite element calculations was carried out on six<br />

tube models crushed axially in a quasi-static condition by using ABAQUS/Explicit. The<br />

cross-section shapes <strong>of</strong> the six tubes were circle, ellipse, triangle, square, pentagon and<br />

hexagon. Effect <strong>of</strong> the side breadth and the generated fold depth on the initiated peak load and<br />

the mean crushing force <strong>of</strong> different tubes were carefully examined.<br />

<strong>Finite</strong> <strong>Element</strong> Model<br />

The finite element simulation <strong>of</strong> the case described herein was carried out using the<br />

commercial code ABAQUS/Explicit version 6.4. Six tubes <strong>with</strong> different geometrical crosssections<br />

were modeled in the present study. All the tube models were developed <strong>with</strong> same<br />

lengths and equal cross-section perimeters <strong>of</strong> 300 mm each. A uniform thickness <strong>of</strong> 1 mm is<br />

considered for all the wall tubes. The difference between the whole tubes was only the shape<br />

<strong>of</strong> the cross-sections. The latter were selected from uniform geometries like a circle, an<br />

ellipse, a triangle, a square, a pentagon and a hexagon. Detailed description <strong>of</strong> the tube crosssections<br />

is illustrated in Fig. (1). For comparing reasons; the material used in the simulation<br />

was provided from the actual material properties <strong>of</strong> steel experimentally measured by Paik et<br />

al. [1]. The mechanical properties <strong>of</strong> the steel material were assigned into the ABAQUS input


9/19<br />

Paper: ASAT-13-ST-34<br />

file. A list <strong>of</strong> numerical values <strong>of</strong> the basic material properties required for the tube models is<br />

given in Fig. (2). All the tube models used in the present finite element simulation were<br />

generated by using the element S4R. This element is a three-dimensional doubly curved fournode<br />

shell element. Each node has three-displacement and three-rotation degrees <strong>of</strong> freedom.<br />

Moreover, this element is considered a general-purpose shell element where it allows for large<br />

strains as load increases. In order to predict the overall response accurately, the mesh <strong>of</strong> finite<br />

elements was fine and uniform <strong>with</strong> equal number <strong>of</strong> elements along the length <strong>of</strong> the tube<br />

and through its circumferential direction. On the other side, all the tube models were crushed<br />

axially between upper and lower rigid parallel plates. These two parallel plates were<br />

simulated by using three-dimensional four-node rigid elements R3D4. The tubes were crushed<br />

by pushing down the upper rigid plate. The latter was fixed in all degrees <strong>of</strong> freedom except<br />

the vertical downward displacement where the loading was applied. However, the lower rigid<br />

plate was stationary by constraining its whole degrees <strong>of</strong> freedom. All tube models rested free<br />

on the lower rigid plate however the applied force was attained as a reaction created from the<br />

pushing <strong>of</strong> the upper plate. The quasi-static loading condition was achieved by moving the<br />

upper plate slowly downward over a sufficiently long time. It is noted that the tubes were not<br />

completely free due to the effect <strong>of</strong> friction between the tube models and the crushing parallel<br />

plates. The progressive deformation shape was continuously monitored for each tube and the<br />

corresponding force-axial displacement curve was depicted.<br />

Experimental Validation<br />

The validation <strong>of</strong> the finite element tube model was made by direct comparison <strong>with</strong> the<br />

experimental results and the collapse observations <strong>of</strong> the square tube US-2 subjected to quasistatic<br />

axial compression and examined by Paik et al. [1]. The structural geometry <strong>of</strong> the<br />

specimen and the properties <strong>of</strong> the used material are shown in Fig. (2). All other testing<br />

details used in the experimental works were necessarily simulated in finite element model.<br />

The progressive collapse <strong>of</strong> the finite element tube model and the examined specimen at<br />

various stages <strong>of</strong> the compression are compared in Fig. (3). Visual examination <strong>of</strong> the pictures<br />

showed good agreements between the numerical simulations and the experimental outcomes.<br />

Comparing the predicted load-displacement behaviour <strong>with</strong> the experimental results in [1] and<br />

the obtained calculations in [2] revealed good correlation between curves at the same<br />

displacement values as shown in Fig. (4).<br />

Table (1) lists the comparison <strong>of</strong> results between the previously published work and the<br />

present finite element analysis.<br />

Investigation <strong>of</strong> <strong>Crushing</strong> <strong>Behaviour</strong><br />

The present numerical simulation provides an opportunity to predict the fold formation <strong>of</strong><br />

various tubes. Figures (5-10) display the collapse mechanisms accompanied by the<br />

corresponding axial- and energy-displacement curves <strong>of</strong> different tube models. Folding modes<br />

were symmetric about the principal axes <strong>of</strong> all the tubes cross-sections. Moreover, an axisymmetric<br />

mode was noticed during the deformation <strong>of</strong> the circular tube. The entire tube<br />

models initiated folding from the top end in contact <strong>with</strong> the movable upper plate. However,<br />

the circular tube began to fold from the bottom end in contact <strong>with</strong> the fixed lower plate. It is<br />

worth mentioning that, the folding process <strong>of</strong> polygonal tubes was initiated by random interior<br />

and exterior rotating <strong>of</strong> the upper edges <strong>of</strong> sidewalls. <strong>Tubes</strong> <strong>of</strong> square and hexagonal crosssections<br />

created symmetric inward and outward folds on each two opposite sidewalls. Folds<br />

<strong>of</strong> the triangular tube initiated by curving its upper edges completely inward. However, the


10/19<br />

Paper: ASAT-13-ST-34<br />

initial inward folding <strong>of</strong> the pentagonal tube was appeared in only two non-successive<br />

sidewalls. After complete constituting <strong>of</strong> the initial folds, similar sequential folds were<br />

regularly created and the crushing was propagated along the sidewalls <strong>of</strong> all polygonal tubes.<br />

The folds generated in all tube models were in contact and adjoining each other except for<br />

that ones <strong>of</strong> the triangular tube. This may be attributed to the acute inner angles between the<br />

sidewalls which increase the strength <strong>of</strong> the triangular tube and allow the rigid body<br />

behaviour to appear early before occurring a complete adjoining <strong>of</strong> folds. On the other hand,<br />

the crushing behaviour <strong>of</strong> the elliptical tube was quite complex and totally different from the<br />

deformation characteristics <strong>of</strong> the other tube models. The upper cross-section <strong>of</strong> the elliptical<br />

tube was initially stretched in the direction <strong>of</strong> maximum diameter as shown in Fig. (6).<br />

Afterwards, the crushing was propagated suddenly allover a region bounded by the upper<br />

stretched cross-section and another similar one located nearly at the middle height <strong>of</strong> the tube.<br />

Progressive compression generated similar crushing characteristics on the rest <strong>of</strong> the elliptical<br />

tube model. When the walls <strong>of</strong> the tubes stopped folding and behaved as a rigid body, the<br />

number <strong>of</strong> folds was observed maximum <strong>with</strong> fourteen folds in the circular tube and<br />

minimum <strong>with</strong> three folds in the triangular tube.<br />

Load- and Energy-Displacement Characteristics<br />

Figures (5-10) indicate the initiation and end <strong>of</strong> the typical axial load- and the energydisplacement<br />

curves <strong>of</strong> various tube models. The absorbed energy was obtained by integrating<br />

the area under the axial load-displacement curve until the folding <strong>of</strong> the walls ended.<br />

Afterwards, the mean crushing load was calculated by dividing the absorbed energy by the<br />

total displacement. Maintaining all parameters <strong>of</strong> the finite element model, a series <strong>of</strong> the<br />

load-displacement curves were obtained from different shape tube models. The general loaddisplacement<br />

characteristics <strong>of</strong> all tubes can be described by a rapid rise <strong>of</strong> the force initiated<br />

due to the elastic compression <strong>of</strong> the tubes. After that a rapid drop <strong>of</strong> the force was observed<br />

when one end <strong>of</strong> the tube models started folding. A steep rising <strong>of</strong> the force was again noticed<br />

after the first folding <strong>of</strong> walls ended and the second folding <strong>of</strong> the adjacent walls started<br />

generation. A drop <strong>of</strong> the force occurred again when a new folding started, and the force rose<br />

if the walls came into contact. This response was repeated until the folding ended and the<br />

models behaved as a rigid body. A uniform response was noted in all tube curves except for<br />

the elliptical tube, the response was irregular. It is worth mentioning that all the numerical<br />

results <strong>of</strong> the energy-displacement relations were nearly linear all over the crushing intervals.<br />

Comparing the force-displacement curves seen in Figs. (5-10) can help analyze the behaviours<br />

<strong>of</strong> different crushing mechanisms. It was clearly seen in Fig. (5) that the higher number <strong>of</strong><br />

local peaks along the force-displacement curve reflected the great number <strong>of</strong> folds created in<br />

the circular tube. The difference between the upper and lower bounds <strong>of</strong> the load values was<br />

significantly large in the square and triangular tube curves and quite narrow in the hexagonal<br />

tube curve. For the pentagonal tube, the difference between upper and lower bounds <strong>of</strong> the<br />

load values started large then damped gradually <strong>with</strong> the compression progress until reaching<br />

a stable load as shown in Fig.(8). On the other side, the general crushing response <strong>of</strong> the<br />

elliptical tube was totally different where the values <strong>of</strong> the load and consequently the<br />

absorbed energy were pronouncedly high and the load peaks were scattered in a distinct<br />

manner all over the curve. This may be attributed to the disturbance occurred during the<br />

folding process <strong>of</strong> this tube shape. Figure (11) shows that the absorbed energy <strong>of</strong> the finite<br />

element tube model decreases <strong>with</strong> an increase in the tube side breadth-to-thickness ratio up<br />

to the tube <strong>of</strong> four sides. Then, the absorbed energy starts increasing at the triangular tube.<br />

The relationship between the absorbed energy and the side breadth-to-thickness ratio <strong>of</strong> the<br />

tubes can be empirically expressed as:


11/19<br />

Paper: ASAT-13-ST-34<br />

E 5.19 1.48cos[(0.05 Sb/ t)<br />

0.6]<br />

(1)<br />

where;<br />

E = Absorbed energy<br />

Sb = Side breadth,<br />

T = Tube thickness<br />

Initial Peak Load and Mean <strong>Crushing</strong> Force<br />

Figure (12) depicts the nonlinear relation between the increasing <strong>of</strong> the number <strong>of</strong> sides and<br />

the decreasing <strong>of</strong> the side breadth <strong>of</strong> a polygonal tube <strong>of</strong> 300mm perimeter. Variations <strong>of</strong> the<br />

initial peak load and mean crushing force <strong>with</strong> different number <strong>of</strong> sides (3-6) and various<br />

side breadths (50-100 mm) are shown in Figs. (13,14) for the polygonal tubes. Generally, it<br />

was observed that the initial peak load highly increased <strong>with</strong> the increasing <strong>of</strong> the number <strong>of</strong><br />

sides and greatly decreased <strong>with</strong> the increasing <strong>of</strong> the side breadth. Moreover, the initial peak<br />

loads <strong>of</strong> the triangular and square tubes had closed magnitudes while values <strong>of</strong> the initial peak<br />

loads <strong>of</strong> the pentagonal and hexagonal tubes were sharply far <strong>of</strong>f. On the other side, the mean<br />

crushing load was changed in a different manner <strong>with</strong> the number <strong>of</strong> sides and the side<br />

breadth <strong>of</strong> the tubes. The mean crushing load greatly decreased <strong>with</strong> the increasing <strong>of</strong> both the<br />

number <strong>of</strong> sides and the side breadth until reaching a minimum value at the square tube.<br />

Subsequently, the mean crushing load began to increase <strong>with</strong> the increasing <strong>of</strong> the number <strong>of</strong><br />

sides and the side breadth <strong>of</strong> the tubes. It is worth mentioning that the initial peak load and<br />

mean crushing force <strong>of</strong> the round (circular and elliptical) tubes were greatly higher than that<br />

obtained from the polygonal tubes as compared in Table (2).<br />

The Fold Depth and the Side Breadth <strong>of</strong> <strong>Tubes</strong><br />

The fold depth was investigated through a longitudinal path <strong>of</strong> the tube. This path was located<br />

at the middle <strong>of</strong> one arbitrary face for all the sided tubes. The fold depth <strong>of</strong> the sided tubes<br />

was compared <strong>with</strong> that obtained from the circular tube overall the tube span as shown in Fig.<br />

(15). The bottom end <strong>of</strong> all tubes was located at the origin <strong>of</strong> the horizontal axis, however the<br />

inward and/or outward fold pr<strong>of</strong>iles were recorded on the vertical axis <strong>of</strong> Fig. (15). It can be<br />

noticed that the folds <strong>of</strong> the circular tube have the lowest depth (approximately 8 mm) and<br />

bulged mainly outward. The pr<strong>of</strong>iles <strong>of</strong> all the folds are not so different except for the<br />

pentagonal tube where its fold pr<strong>of</strong>ile was not symmetric about the middle path. On the other<br />

side and apart from the triangular tube, the fold depth was seen maximum in the square tube<br />

and it generally decreased <strong>with</strong> increasing the number <strong>of</strong> sides <strong>of</strong> the polygonal tubes. In other<br />

words, starting from a polygonal tube <strong>with</strong> four sides, the initiated fold depth is decreased<br />

<strong>with</strong> the decreasing <strong>of</strong> the tube side breadth as depicted in Fig. (16). The relationship between<br />

the fold depth and the side breadth-to-thickness ratio <strong>of</strong> the tubes can be empirically<br />

expressed as:<br />

where;<br />

Fd = Fold depth,<br />

Sb = Side breadth,<br />

t = Tube thickness<br />

Fd 8.302exp[0.021( Sb / t)]<br />

(2)


12/19<br />

Paper: ASAT-13-ST-34<br />

Figure (16) shows that the results <strong>of</strong> the fold depth formula coincide <strong>with</strong> the numerical data.<br />

This relation is valid for the polygonal tubes <strong>with</strong> number <strong>of</strong> sides greater than or equal four<br />

sides.<br />

Conclusion<br />

The non-linear finite element model simulations provide much more detailed information than<br />

other approaches. These are especially efficient in local folding and internal contact <strong>of</strong><br />

structural elements. In this paper, a study on the quasi-static crushing collapse <strong>of</strong> multi-sided<br />

and round tubes was performed using a non-linear finite element program ABAQUS/Explicit.<br />

Based on the numerical data, simple empirical formulas were developed to predict the<br />

absorbed energy and the fold depth as a function <strong>of</strong> the side breadth-to-thickness ratio <strong>of</strong> the<br />

tubes. From the present study, the following conclusion can be drawn:<br />

- All the tube models were crushed in symmetric concertina modes.<br />

- The greatest number <strong>of</strong> folds was observed fourteen in the circular cylinder and minimum<br />

<strong>with</strong> three folds in the triangular tube.<br />

- The folds <strong>of</strong> the cylindrical tube were compact <strong>with</strong> short depth while the triangular tube<br />

terminated folding <strong>with</strong>out contact between folds.<br />

- The initial peak load increases <strong>with</strong> the increasing <strong>of</strong> the no. <strong>of</strong> sides <strong>of</strong> the tube and<br />

decreases <strong>with</strong> the increasing <strong>of</strong> the side breadth.<br />

- The greatest energy absorption capacity was obtained by the circular cylinder while the<br />

square tube absorbed minimum energy.<br />

- The mean crushing load recorded lowest value in case <strong>of</strong> the square tube however the<br />

highest magnitude was noticed in the triangular tube.<br />

References<br />

[1] J. K. Paik, J. Y. Chung, and M. S. Chun, “On quasi-static crushing <strong>of</strong> a stiffened square<br />

tube”, J. Ship Res., Vol. 40, No. 3, pp. 258–267, 1996.<br />

[2] A. Zhang, and K. Suzuki, “A study on the effect <strong>of</strong> stiffeners on quasi-static crushing <strong>of</strong><br />

stiffened square tube <strong>with</strong> non-linear finite element method”, Int. J. <strong>of</strong> Imp. Eng., Vol.<br />

34, pp. 544–555, 2007.<br />

[3] Gupta N. K., Sekhon G. S. and Gupta P. K., “A Study <strong>of</strong> Formation in Axi-symmetric<br />

Axial Collapse <strong>of</strong> Round <strong>Tubes</strong>”, Int. J. Impact. Eng., Vol. 27, pp. 87-117, 2002.<br />

[4] Gupta N. K. and Gupta S. K., “Effect <strong>of</strong> Annealing, Size and Cutouts on Axial Collapse<br />

<strong>Behaviour</strong> <strong>of</strong> Circular <strong>Tubes</strong>”, Int. J. <strong>of</strong> Mech. Sci., Vol. 35, pp. 519-613, 1993.<br />

[5] Wierzbicki T. and Abramowicz W., “On the <strong>Crushing</strong> Mechanics <strong>of</strong> <strong>Thin</strong>-Walled<br />

Structures”, J. Appl. Mech., Vol. 50, No. 4a, pp. 727–34, 1983.<br />

[6] Abramowicz W. and Jones N., “Dynamic Axial <strong>Crushing</strong> <strong>of</strong> Square <strong>Tubes</strong>”, Int. J.<br />

Impact Eng., Vol. 2, No. 2, pp. 179–208, 1984.<br />

[7] Abramowicz W. and Jones N., “Dynamic Progressive Buckling <strong>of</strong> Circular and Square<br />

<strong>Tubes</strong>”, Int. J. Impact Eng., Vol. 4, No. 4, pp. 243–70, 1986.<br />

[8] Grzebeita R. H., “An Alternate Method for Determining the <strong>Behaviour</strong> <strong>of</strong> Round Stocky<br />

<strong>Tubes</strong> Subjected to an Axial Crush Load”, <strong>Thin</strong> Wall Struct., Vol., 9, pp. 61-89, 1990.<br />

[9] Gupta N. K., Velmurugan R., “Consideration <strong>of</strong> Internal Folding and non Symmetric<br />

Fold Formation in Axi-symmetric Axial Collapse <strong>of</strong> Round <strong>Tubes</strong>”, Int. J. Solids<br />

Struct., Vol. 34, No. 20, pp. 2611-40, 1997.


13/19<br />

Paper: ASAT-13-ST-34<br />

[10] Wierzbicki T., Bhat S. U., Abramowicz W. and Brodikin D., “A Two Folding <strong>Element</strong>s<br />

Model <strong>of</strong> Progressive <strong>Crushing</strong> <strong>of</strong> <strong>Tubes</strong>”, Int. J. Solids Struct., Vol. 29, No. 24, pp.<br />

3269-88, 1992.<br />

[11] Singace A. A. and EI-Sobky H., “Further Experimental Investigation on the Eccentricity<br />

Factor in the Progressive <strong>Crushing</strong> <strong>of</strong> <strong>Tubes</strong>”, Int. J. Solids Struct., Vol. 33, No. 24, pp.<br />

3517-38, 1996.<br />

[12] Singace A. A, EI-Sobky H. and Reddy T. Y., “On the Eccentricity Factor in the<br />

Progressive <strong>Crushing</strong> <strong>of</strong> <strong>Tubes</strong>”, Int. J. Solids Struct., Vol. 32, No. 24, pp. 3589-602,<br />

1995.<br />

[13] Singace A. A., “Axial <strong>Crushing</strong> Analysis <strong>of</strong> <strong>Tubes</strong> Deforming in the Multi-Lobe<br />

Mode”, Int. J. Mech. Sci., Vol. 41, pp. 865-90, 1999.<br />

[14] Hull D., “A Unified Approach to Progressive <strong>Crushing</strong> <strong>of</strong> Fiber-Reinforced Composite<br />

<strong>Tubes</strong>”, Comp. Sci. Technol., Vol. 40, pp. 377-421, 1991.<br />

[15] Farley G. L., “Effects <strong>of</strong> Specimen Geometry on the Energy Absorption Capability <strong>of</strong><br />

Composite Materials”, J. Comp. Mater., Vol. 20, pp. 390-400, 1986.<br />

[16] Mamalis A. G., Manolakos D. E., Demosthenous G. A. and Ioannidis M. B., “The<br />

Deformation Mechanism <strong>of</strong> <strong>Thin</strong>-Walled Non-circular Composite <strong>Tubes</strong> Subjected to<br />

Bending”, Comp. Struct., Vol. 30, pp. 131-146, 1995.<br />

[17] Mamalis A. G., Manolakos D. E., Demosthenous G. A. and Ioannidis M. B., “Analysis<br />

<strong>of</strong> Failure Mechanisms Observed in Axial Collapse <strong>of</strong> <strong>Thin</strong>-Walled Circular Fiberglass<br />

Composite <strong>Tubes</strong>”, <strong>Thin</strong>-Walled Struct., Vol. 24, pp. 335-352, 1996.<br />

[18] Mamalis A. G., Manolakos D. E., Demosthenous G. A. and Ioannidis M. B.,<br />

“Analytical <strong>Modeling</strong> <strong>of</strong> the Static and Dynamic Axial Collapse <strong>of</strong> <strong>Thin</strong>-Walled<br />

Fiberglass Composite Conical Shells”, Int. J. Impact Eng., Vol. 19, No. 5-6, pp. 477-<br />

492, 1997.<br />

75 mm 100 mm 95.49 mm<br />

60 mm 50 100 mm<br />

Fig. (1) Different shapes <strong>of</strong> uniform cross-sections used<br />

in the F.E. analysis <strong>of</strong> a 300 mm perimeter tube.<br />

50 mm


Dimensions <strong>of</strong><br />

the test specimen<br />

L = 450 mm<br />

b = 100 mm<br />

t = 2.8 mm<br />

L<br />

b<br />

t<br />

Fig. (2) Configuration <strong>of</strong> the test specimen US-2 in Ref. [1].<br />

Fig. (3) Gradual crushing <strong>of</strong> the square tube US-2 [1] using the<br />

present ABAQUS solution.<br />

Fig. (4) Comparison between the crushing curves <strong>of</strong> the<br />

square tube US-2 [1].<br />

14/19<br />

Paper: ASAT-13-ST-34<br />

Material properties required for<br />

finite element model<br />

Density = 7860 kg/m 3<br />

Young’s modulus = 205.8 GPa<br />

Poisson’s ratio = 0.3<br />

Yield strength = 310 MPa<br />

Ultimate stress = 363.3 MPa<br />

Fracture strain = 0.354<br />

1<br />

2


Load(kN)<br />

15/19<br />

Paper: ASAT-13-ST-34<br />

Fig. (5) Gradual collapse and crushing results <strong>of</strong> the cylindrical tube.<br />

Load (kN)<br />

Fig. (6) Gradual collapse and crushing results <strong>of</strong> the elliptical tube.<br />

Load (kN)<br />

Fig. (7) Gradual collapse and crushing results <strong>of</strong> the hexagonal tube.<br />

100<br />

80<br />

60<br />

40<br />

20<br />

0<br />

450<br />

360<br />

270<br />

180<br />

90<br />

0<br />

60<br />

48<br />

36<br />

24<br />

12<br />

0<br />

Load<br />

Energy<br />

0 50 100 150 200 250<br />

Displacement (mm)<br />

Load<br />

0 50 100 150 200 250<br />

Displacement (mm)<br />

Load<br />

Energy<br />

Energy<br />

0 50 100 150 200 250<br />

Displacement (mm)<br />

7<br />

5.6<br />

4.2<br />

2.8<br />

1.4<br />

0<br />

45<br />

36<br />

27<br />

18<br />

9<br />

0<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

Energy (kJ)<br />

Energy (kJ)<br />

Energy (kJ)


Load (kN)<br />

45<br />

36<br />

27<br />

18<br />

9<br />

0<br />

Paper: ASAT-13-ST-34<br />

0 50 100 150 200 250<br />

Displacement (mm)<br />

Fig. (8) Gradual collapse and crushing results <strong>of</strong> the pentagonal tube.<br />

Load (kN)<br />

Fig. (9) Gradual collapse and crushing results <strong>of</strong> the square tube.<br />

Load (kN)<br />

40<br />

32<br />

24<br />

16<br />

40<br />

32<br />

24<br />

16<br />

8<br />

0<br />

8<br />

0<br />

Energy<br />

Load<br />

0 50 100 150 200 250<br />

Displacement (mm)<br />

Load<br />

Energy<br />

Energy<br />

Load<br />

0 50 100 150 200 250<br />

Displacement (mm)<br />

Fig. (10) Gradual collapse and crushing results <strong>of</strong> the triangular tube.<br />

16/19<br />

4<br />

0<br />

4.5<br />

3.6<br />

2.7<br />

1.8<br />

0.9<br />

0<br />

3.2<br />

2.4<br />

1.6<br />

0.8<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

Energy (kJ)<br />

Energy (kJ)<br />

Energy (kJ)


Fig. (11) Variation <strong>of</strong> the absorbed energy<br />

<strong>of</strong> polygonal tubes <strong>with</strong> the side<br />

breadth-to-thickness ratio.<br />

Initial peak load<br />

Mean crushing load<br />

55<br />

50<br />

45<br />

40<br />

35<br />

30<br />

2 3 4 5 6 7<br />

No. <strong>of</strong> sides<br />

27<br />

24<br />

21<br />

18<br />

15<br />

Side breadth<br />

17/19<br />

Paper: ASAT-13-ST-34<br />

0<br />

0 20 40 60<br />

No. <strong>of</strong> sides<br />

Fig. (13) Variation <strong>of</strong> the initial peak load <strong>with</strong> the number <strong>of</strong> sides<br />

(left) and <strong>with</strong> the side breadth (right) <strong>of</strong> the polygonal tubes.<br />

12<br />

2 3 4 5 6 7<br />

No. <strong>of</strong> sides<br />

120<br />

100<br />

80<br />

60<br />

40<br />

20<br />

Fig. (12) Variation <strong>of</strong> the side breadth <strong>of</strong><br />

a polygonal tube <strong>of</strong> 300mm perimeter <strong>with</strong><br />

the number <strong>of</strong> sides.<br />

Initial peak load<br />

Mean crushing Load<br />

55<br />

50<br />

45<br />

40<br />

35<br />

30<br />

40 60 80 100<br />

side breadth (mm)<br />

27<br />

24<br />

21<br />

18<br />

15<br />

12<br />

40 60 80 100 120<br />

Side breadth (mm)<br />

Fig. (14) Variation <strong>of</strong> the mean crushing load <strong>with</strong> the number <strong>of</strong><br />

sides (left) and <strong>with</strong> the side breadth (right) <strong>of</strong> the polygonal tubes


Fig. (15) Variation <strong>of</strong> the fold formation along axial path at<br />

the mid side breadth <strong>of</strong> tubes.<br />

Fig. (16) Variation <strong>of</strong> the fold depth <strong>with</strong> the side<br />

breadth-to-thickness ratio <strong>of</strong> the tube models.<br />

18/19<br />

Paper: ASAT-13-ST-34


Table (1) Comparison between the presented results <strong>of</strong> F.E.M.,<br />

experimental results in Ref.[1] and numerical solution in Ref.[2]<br />

Parameter Ref.[1] Ref.[2] F.E.M.<br />

Peak load (kN) 337.12 357.23 348.25<br />

Mean load (kN) 100.98 101.32 102.78<br />

Table (2) Effect <strong>of</strong> variation <strong>of</strong> the tube cross-section on<br />

the obtained finite element results<br />

19/19<br />

Paper: ASAT-13-ST-34<br />

Cross-section<br />

Max. load<br />

(kN)<br />

Mean load<br />

(kN)<br />

Max. deformation<br />

(mm)<br />

Absorbed energy<br />

(kJ)<br />

Triangular 34.07 23.53 191 4.41<br />

Squarer 34.43 16.84 234.73 3.75<br />

Pentagonal 38.64 19.47 217 4.22<br />

Hexagonal 52.59 22.76 218.92 4.90<br />

Circular 92.92 37.82 235.39 6.42<br />

Elliptical 268.98 162.71 233.29 40.59

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!