IMPACT STRENGTH AND RESPONSE BEHAVIOR OF CFRP GUARDER BELT FOR SIDE COLLISION OF AUTOMOBILES

Similar documents
COMPARISON OF EXPERIMENTAL RESULTS WITH FEM ONES OF RECTANGULAR CFRP TUBES FOR FRONT SIDE MEMBERS OF AUTOMOBILES

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

IDE 110 Mechanics of Materials Spring 2006 Final Examination FOR GRADING ONLY

Impact Analysis of Frontal Car Bumper using Long Fibre Reinforced Thermoplastics

3 Nihon University Izumicho, Narashino, Chiba, , JAPAN

ANALYTICAL PENDULUM METHOD USED TO PREDICT THE ROLLOVER BEHAVIOR OF A BODY STRUCTURE

BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS

Failure analysis of serial pinned joints in composite materials

The University of Melbourne Engineering Mechanics

Module III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina

Lab Exercise #5: Tension and Bending with Strain Gages

Experimentally Calibrating Cohesive Zone Models for Structural Automotive Adhesives

Modeling non-isothermal thermoforming of fabricreinforced thermoplastic composites

Numerical simulation of the crash behaviour of braided composites

Crash Simulation and Analysis of a Car Body Using ANSYS LS- DYNA

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3

Mechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002

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

Calibration and Experimental Validation of LS-DYNA Composite Material Models by Multi Objective Optimization Techniques

Samantha Ramirez, MSE. Stress. The intensity of the internal force acting on a specific plane (area) passing through a point. F 2

IMECE CRASHWORTHINESS OF AIRCRAFT COMPOSITES STRUCTURES

Computational Analysis for Composites

NUMERICAL SIMULATION OF DAMAGE IN THERMOPLASTIC COMPOSITE MATERIALS

An orthotropic damage model for crash simulation of composites

Comparison of Ply-wise Stress-Strain results for graphite/epoxy laminated plate subjected to in-plane normal loads using CLT and ANSYS ACP PrepPost

Special edition paper

Effect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test

INCREASING RUPTURE PREDICTABILITY FOR ALUMINUM

MECHANICAL FAILURE OF A COMPOSITE HELICOPTER STRUCTURE UNDER STATIC LOADING

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.

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

NONLINEAR STATIC SIMULATION OF AUTOMOTIVE BUMPER OF A PASSENGER CAR IN LOW SPEED IMPACT CRASH SIMULATION

Quasi-static simulation approaches on rollover impact of a bus structure

Structural Metals Lab 1.2. Torsion Testing of Structural Metals. Standards ASTM E143: Shear Modulus at Room Temperature

Stress-Strain Behavior

International Journal of Scientific and Research Publications, Volume 5, Issue 8, August ISSN

THE EFFECT OF GEOMETRY ON FATIGUE LIFE FOR BELLOWS

Purpose of this Guide: To thoroughly prepare students for the exact types of problems that will be on Exam 3.

SCALING EFFECTS IN THE LOW VELOCITY IMPACT RESPONSE OF FIBRE METAL

Crashworthiness of composite structures: Experiment and Simulation

An investigation of the mechanical behaviour of carbon epoxy cross ply cruciform specimens under biaxial loading

Finite Element Simulation of Bar-Plate Friction Welded Joints Steel Product Subjected to Impact Loading

MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING

Mechanical Properties of Materials

PREDICTION OF FAILURE BEHAVIOUR OF COMPOSITE LATTICE STRUCTURE UNDER COMPRESSIVE LOAD

Failure Analysis of Unidirectional Composite Pinned- Joints

FLEXURAL RESPONSE OF FIBER RENFORCED PLASTIC DECKS USING HIGHER-ORDER SHEAR DEFORMABLE PLATE THEORY

External Pressure... Thermal Expansion in un-restrained pipeline... The critical (buckling) pressure is calculated as follows:

EXPERIMENTAL AND NUMERICAL STUDY OF THE ENERGY ABSORPTION CAPACITY OF PULTRUDED COMPOSITE TUBES

COMPARISON OF NUMERICAL SIMULATION AND EXPERIMENT OF A FLEXIBLE COMPOSITE CONNECTING ROD

Stress Analysis Lecture 3 ME 276 Spring Dr./ Ahmed Mohamed Nagib Elmekawy

N = Shear stress / Shear strain

[5] Stress and Strain

Mechanical Behavior of Circular Composite Springs with Extended Flat Contact Surfaces

Analysis Of Naca 2412 For Automobile Rear Spoiler Using Composite Material *

SOME RESEARCH ON FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS

OPTIMAL FIBER PLACEMENT INCLUDING EFFECTS OF EMBROIDERY

ME 243. Mechanics of Solids

Experimental Approach to Determine the Stress at a Section of Semi Circular Curved Beam Subjected to Out-Of-Plane Load Using Strain Rosette

Effect of Deformation Mode of Cylindrical Tubes with Corrugated Surface Dimensional Errors

Parameter Design of High Speed Coupling for 6 MW Wind Turbine Considering Torsional Vibration

BTECH MECHANICAL PRINCIPLES AND APPLICATIONS. Level 3 Unit 5

Progressive Damage of GFRP Composite Plate Under Ballistic Impact: Experimental and Numerical Study

SHEAR CAPACITY OF REINFORCED CONCRETE COLUMNS RETROFITTED WITH VERY FLEXIBLE FIBER REINFORCED POLYMER WITH VERY LOW YOUNG S MODULUS

Experimental Test and Computer Simulation Research on Rollover Impact of a Bus Structure

A ROLLOVER TEST OF BUS BODY SECTIONS USING ANSYS

RATE-DEPENDENT OFF-AXIS COMPRESSIVE STRENGTH OF A UNIDIRECTIONAL CARBON/EPOXY LAMINATE AT HIGH TEMPERATURE

Strength of Material. Shear Strain. Dr. Attaullah Shah

Discrete Element Modelling of a Reinforced Concrete Structure

KINK BAND FORMATION OF FIBER REINFORCED POLYMER (FRP)

Static and Time Dependent Failure of Fibre Reinforced Elastomeric Components. Salim Mirza Element Materials Technology Hitchin, UK

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011

Study on the Energy Absorption and the Local Buckling of the Expansion Tubes

Name :. Roll No. :... Invigilator s Signature :.. CS/B.TECH (CE-NEW)/SEM-3/CE-301/ SOLID MECHANICS

Introduction to Engineering Materials ENGR2000. Dr. Coates

Impact Simulation of Extreme Wind Generated Missiles on Radioactive Waste Storage Facilities

Crash and Impact Simulation of Composite Structures by Using CAE Process Chain

DESIGNING COMPOSITE STRUCTURES FOR IMPACT PERFORMANCE WHAT CAN WE LEARN FROM THE AEROSPACE INDUSTRY?

EDEXCEL NATIONAL CERTIFICATE/DIPLOMA SCIENCE FOR TECHNICIANS OUTCOME 1 - STATIC AND DYNAMIC FORCES TUTORIAL 3 STRESS AND STRAIN

Analysis of Composite Pressure Vessels

NUMERICAL AND EXPERIMENTAL STUDY OF FAILURE IN STEEL BEAMS UNDER IMPACT CONDITIONS

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

NUMERICAL MODELLING OF COMPOSITE PIN- JOINTS AND EXPERIMENTAL VALIDATION

Solution: T, A1, A2, A3, L1, L2, L3, E1, E2, E3, P are known Five equations in five unknowns, F1, F2, F3, ua and va

PREDICTION OF OUT-OF-PLANE FAILURE MODES IN CFRP

Design of Beams (Unit - 8)

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

FINITE ELEMENT AND EXPERIMENTAL STUDY OF NOVEL CONCEPT OF 3D FIBRE CELL STRUCTURE

Open-hole compressive strength prediction of CFRP composite laminates

NUMERICAL MODELING OF GAS LEAKAGE THROUGH DAMAGED COMPOSITE LAMINATES

4.MECHANICAL PROPERTIES OF MATERIALS

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

Sabah Shawkat Cabinet of Structural Engineering Walls carrying vertical loads should be designed as columns. Basically walls are designed in

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

MECHANICAL PROPERTIES OF LIGHTWEIGHT COMPOSITES REINFORCED WITH MICRO GLASS BALLOONS

Evaluation of in-plane orthotropic elastic constants of paper and paperboard

THREE DIMENSIONAL STRESS ANALYSIS OF THE T BOLT JOINT

ME 202 STRENGTH OF MATERIALS SPRING 2014 HOMEWORK 4 SOLUTIONS

2. Rigid bar ABC supports a weight of W = 50 kn. Bar ABC is pinned at A and supported at B by rod (1). What is the axial force in rod (1)?

2012 MECHANICS OF SOLIDS

Transcription:

IMPACT STRENGTH AND RESPONSE BEHAVIOR OF CFRP GUARDER BELT FOR SIDE COLLISION OF AUTOMOBILES AOKI Y.*, BEN G.**, KIM H. S.* *College of Science & Technology, Nihon University **College of Industrial Technology, Nihon University 741 Narashinodai, Funabashi, Chiba 748501 JAPAN aoki.yoshio@nihonu.ac.jp SUMMARY In this paper, the CFRP guarder belt equipped in the automotive door is developed and examined by an experiment and a numerical analysis for replacing the conventional steel door guarder beam. As the experimental relation of impact load to displacement for CFRP guarder belt agreed well with that of numerical result, the numerical method developed here is quite useful for estimating impact behaviours of CFRP guarder belt. Keywords: Impact behaviour, CFRP, Absorbed energy, Automobile, Collision safety INTRODUCTION It is well known that CO emissions, which are one of the greenhouse gases emitted from passenger vehicles such as automobile, are major cause of global warming. In the automobile industry, to reduce CO emissions, it is well known that the most effective method is to produce the fuel efficient automobile. To increase the fuel efficiency of the automobile, the most effective approach is to reduce the automobile weight by using lightweight material such as composite materials. Carbon fiber reinforced plastics (CFRP) have been widely used in aerospace, industrial goods and other application fields because of their high specific strength and high specific modulus compared with metal. This means that the CFRP contributes to lighten automobiles greatly. Otherwise, the safety of automobiles is also very important and the collision safety of the automobile has been evaluated by full flap frontal crash, offset frontal crash and side impact tests. In the frontal crash test, it is possible to absorb the energy by largely deforming the front part and the rear. With an increasing an interest in the lightening of the automobile and in the securing the safety of passengers, many researches for them have been performed [19]. However, in the side impact test, it is hard to absorb the energy similarly, because of being very narrow for the survival space of passengers. In the inside of the door, a reinforcement member as shown in Figure 1, namely door guarder beam made of steel has been installed to absorb impact energy and its deformation is limited to about 150mm. In this study, the CFRP door guarder belt is developed for the purpose of designing impact energy absorption members under side collision as shown in Figure. A drop weight impact tests are carried out to investigate impact response behaviours and impact energy absorption characteristics of the CFRP door guarder belt. Also, a finite element

(FE) model was developed by using the ANSYS to simulate the impact response behaviour and the absorbed energy of the CFRP door guarder belt under impact loading. Rotary pin Tension Impactor Tension Tension Tension CFRP Belt Rotary pin Figure 1 Conventional door guarder beam Figure Energy absorption mechanism of CFRP guarder belt EXPERIMENT Specimen fabrication Figure shows the schematic diagram of the impact energy absorption by the CFRP guarder belt. The impact energy is effectively absorbed by installing the CFRP belt between two free fulcrums for the rotation and changing the vertical impact load of the falling weight to the tensile load. In order to prevent fracture in the support edge of specimen due to the stress concentration, a beltshaped specimen was adopted. Thin CFRP belt specimens were manufactured from unidirectional prepregs (T800H/Epoxy) by using the sheet winding method and its thickness, width and length were 0.3mm, 50mm and 164mm, respectively. Figure 3 shows an actual CFRP guarder belt. Tower drop weight impact test In order to evaluate the capacity of crash energy absorption and to show the micro and macro fracture behavior of the CFRP guarder belt, the large size of drop tower facility for the impact test was designed as shown in Figure 3. The CFRP guarder belt was received the impact load generated by a free drop weight of 100kg from 1m height. Therefore, the impact speed was approximately 55km/h just before the impact. The shape of impactor was a half cylinder having 100mm radius and 00mm width. The Road cell Strain gage (a) CFRP guarder belt specimen (b) mounted specimen

(c) Impactor (d) measured location of the longitudinal strain Impactor Mass=100kg Figure 3 CFRP guarder belt specimen and Tower drop impact test setup impact load of the specimen was measured by a load cell installed behind the rotary pin. In order to investigate the fracture mechanism of the CFRP guarder belt, a high speed camera was chosen. And the dynamic strain of the specimen from collision to fracture was measured by a strain gauge. The strain gauge stuck on three places at the center of the specimen as a collision point (T1), near the rotary pin (T3) and middle point of both (T) as shown in Figure 3. The CFRP guarder belt specimen supported both ends in the rotary pin of a diameter of 40mm. Figure 4 shows the measured horizontal load and the longitudinal strain of the center of the specimen variation after impact. The horizontal impact load of CFRP belt specimen increase nonlinearly, and the impact load reached to the maximum and then instantly dropped to almost zero as shown in the left side of Figure 4. Within the time of 5.0~5.4msec, the impact load recovered to 38~40kN and become to zero at the displacement of about 85mm. The longitudinal strain of the specimen at the impact location also increases nonlinearly from a moment of impact, and becomes the maximum value just before fracture of the CFRP belt specimen. Figure 5 shows the observed fracture modes and the fracture location of CFRP guarder belt specimen. In the experiment, the observed fracture mode was fiber breakage in the entire width of CFRP belt specimen due to the tensile load acted on the whole the

specimen. And almost specimen fractured at the location of T3, that is larger tensile stress occurred around the rotary pin is supposed. 45 40 : Horizontal load 16000 14000 : Longitudinal strain (T1) 35 30 L 5 o 0 a d 15 (kn) 10 S 1000 t 10000 r a 8000 i 6000 n (µε) 4000 5 000 0 0.0 1.0.0 3.0 4.0 5.0 6.0 Time(msec) 0 0.0 1.0.0 3.0 4.0 5.0 6.0 Time(msec) Figure 4 Horizontal load and Strain of the center of the specimen variation after impact Fracture Point Figure 5 Fracture mode of CFRP in the impact test for the CFRP guarder belt Figure 6 shows the comparison of experimental loaddisplacement curve between the vertical impact and the impact with offset angle of 15. The horizontal load of the impact with offset angle of 15 increases nonlinearly until the displacement of 80mm, and it s almost same tendency for the vertical impact case. The maximum impact load for offset impact case is 4% smaller than the vertical impact case. Figure6 Comparison of experimental loaddisplacement curves between the vertical impact and the impact with offset angle of 15

NUMERICAL ANALYSIS To simulate the impact response behavior and absorbed energy of the CFRP guarder belt under impact loading, a finite element model was developed by using the commercial FEM solver (ANSYS). In this FE analysis, the laminated shell element (SHELL_63) and the solid element (SOLID_45) were used in the CFRP guarder belts and the impactor and rotary pins, respectively. Details of the finite element model are shown in Figure 7. The contact element between the impactor and the upper surface of CFRP guarder belt and between the rotary pins and the inner surface of the belt was Contact element 173 with the Target element 170 of the friction coefficients of 0.5. And, we adopted 1/ models with a symmetric condition to reduce a calculation time. Table 1 shows the material properties of CFRP guarder belt used for the analysis. Prediction of the Fracture time of CFRP guarder belt Next, the fracture time of the CFRP guarder belt was predicted from the results of the numerical analysis by using the fracture criterion of composites. In case of the unidirectional reinforced CFRP, the material property in the transverse direction is same as that of thickness direction, therefore TsaiHill fracture criterion was used in this analysis. σ L FL σ σ L T FL σ T + FT τ + F LT LT = 1 (1) Where F L, F T and F LT are the strength of the CFRP belt along the fiber direction, that of the transverse direction and that of the inplane shear, respectively. The F L being much bigger than the stress of the transverse direction σ T for thin CFRP belt specimen, so that equation (1) can be rewritten as below σ L FL σ T + FT τ + F LT LT = 1 () Substituting the occurred stress and the strength in fiber direction σ L, those of transverse σ T and those of membrane shear τ LT in the equation, it can be consider the fracture of CFRP belt specimen when the value of the lefthand side become 1. Calculation of the absorbed energy for CFRP guarder belt The absorbed energy of CFRP guarder belt can be calculated by an area of the curve of an impact load P and a displacement of the impact location of the specimen δ until its fracture as ( δ) = P( δ) E dδ (3)

CFRP Impactor, Belt Coner Pin 3 D Structual Elastic Structural Element Shell Solid 3 D Target Segment Contact Element 3 D Surface to Surface Contact Figure 7 Details of the finite element model Table 1 Material properties of CFRP guarder belt E x [GPa] E y [GPa] E z [GPa] G xy [GPa] G yz [GPa] F x [GPa] F y [GPa] Fxy[GPa] ν Density[kg/m 3 ] CFRP Belt 130 8 8 4 4.45 0.07 0.098 0.34 1600 Impactor, Coner Pin 06 80 0.3 8000 RESULTS AND DISCUSSIONS Figure 8 shows the experimental and the numerical relation of impact load to displacement of the CFRP guarder belt for vertical impact. The experimental and numerical impact load increase with increasing the displacement of CFRP belt specimen nonlinearly, then, they become the maximum value just before fracture. The experimental and the numerical relation of impact load to displacement for CFRP guarder belt agreed well generally. Next, Table shows the predicted fracture time was calculated by equation (), the location of fracture in the specimen and the absorbed energy obtained by FE analysis are compared with the experimental ones. The fracture time and the location of fracture are almost exactly predicted and the numerical absorbed energy is also close to the experimental ones. The specific absorbed energy divided by the individual weights of the CFRP guarder belt was compared with the conventional steel impact beam shown in Table 3. The conventional impact beam consists of steel pipe and has diameter of 35mm, thickness of mm and length of 945mm. Though experimental energy absorption of CFRP guarder belt is smaller than it of the conventional impact beam, the specific energy absorption of the CFRP belt is 30 times bigger than it of the impact beam. Therefore the energy absorption of CFRP belt bigger than the impact beam is possible if the width and thickness of the CFRP belt is increased.

Load(KN) The experimental strain variations at three locations in the specimen compare with the numerical ones as shown in Figure 9. The strain of each locations increase with time transit nonlinearly. The measured maximum impact strain at each location is almost same, though the numerical strain variation at the location of T3 is larger than experimental ones. This is because the numerical result wasn t considered a little looseness between the CFRP guarder belt and the rotary pin. Both results of the strain variations at the locations of T1 and T agreed well. From the comparison of both results, the numerical method developed here is quite useful for estimating impact behavior of the CFRP guarder belt. In the case of side collision, there is the collision from diagonal direction actually. Therefore, the impact response behavior for the side collision with offset angle of CFRP guarder belt was examined. The experiment of the side collision with offset angle was performed to incline the supported base of the specimen so that the impactor hit the specimen diagonally as shown in Figure 10. The drop weight impact test with offset angle of 15 was carried out in this experiment, then the difference in impact response behavior of the side collision was examined. 45 40 35 ; Vertical impact(exp.) ; Vertical impact(fem) 30 5 0 15 10 5 0 0 0 40 60 80 Displacement(mm) Figure 8 Comparison of experimental and predicted loaddisplacement curve of CFRP guarder belt for vertical impact Table Fracture time and absorbed energy of CFRP guarder belt for vertical impact Experiment FEM Fracture time (msec) 5.63 5.7 Location of the fracture T3 T3 Absorbed energy (J) 1118 1010 Table 3 Comparison of Absorbed energy for CFRP guarder belt and steel impact beam Absorbed energy (J) Specific absorbed energy (kj/kg) CFRP guarder belt 1118 37.14 Steel impact beam 1900 1.43

Figure 9 Comparisons of experimental and predicted the strain variation of CFRP guarder belt for vertical impact Figure 11 shows the experimental and the numerical relation of impact load to displacement of the CFRP guarder belt for impact with offset angle of 15. Though the maximum impact load by numerical result is slightly larger than experimental ones, both relation of impact load to displacement for CFRP guarder belt agreed well generally. Table 4 shows the predicted fracture time, the location of fracture in the specimen and the absorbed energy obtained by FE analysis are compared with the experimental ones

Figure 10 The drop weight impact test with offset angle of 15 Figure 11 Comparison of experimental and predicted loaddisplacement curve for impact with offset angle of 15 Table 4 Experimental fracture time and absorbed energy of CFRP guarder belt for impact with offset angle of 15 Experiment F.E.M Fracture time (msec) 5.4 5.93 Location of the fracture T3 T3 Absorbed energy (J) 865 973 for the impact with offset angle of 15. The fracture time and the absorbed energy are almost exactly predicted. The fracture time for the offset impact test is slightly longer than that for vertical impact, but the absorbed energy is almost same for the both case. The experimental strain variations at three locations in the specimen compare with the numerical ones as shown in Figure 1. The strain of each locations increase with time transit nonlinearly. Both results of the strain variations at the locations of T1 and T agreed well. CONCLUSIONS The CFRP guarder belt was developed for the purpose of designing impact energy absorption member under side collision. The drop weight impact tests were carried out and the impact response behavior and the absorbed energy of the CFRP door guarder belt under impact loading were examined by using the numerical analysis and the experimental results. From these results, we could be concluded as below. 1. The CFRP guarder belt absorbed crash energy along the entire length of it and tension stress is applied on both the upper and lower side of belt.. From the comparison of FEM results with the experimental ones for the specimens, the proposed numerical method by ANSYS code was supposed to be useful for analyzing the CFRP guarder belt. 3. The fracture time of the CFRP guarder belt can be predicted by using FEM and fracture criteria of composite materials and agreed well with the experimental fracture time. 4. The impact response behavior and impact strength of the CFRP guarder belt were obtained by the tower drop weight impact test.

The CFRP guarder belt contributes to lightweighting and safety improvement of safety of the car body greatly is supposed. Figure 1 Comparisons of experimental and predicted the strain variation for impact with offset angle of 15 ACKNOWLEDGEMENTS This study was conducted as part of the Japanese National Project "R&D of Carbon FiberReinforced Composite Materials to Reduce Automobile Weight" supported by NEDO (New Energy and Industrial Technology Development Organization). The authors acknowledge the assistance of Toray Industries, Inc. who supplied the materials for these test specimens. References 1. G. Ben, T. Uzawa, H.S. Kim, Y. Aoki et. al., J. JSME Series A, 70 (694), 84, 004 (in Japanese). Ari G. Caliskan., Proceedings of IMECE 00, November17, 00 3. H. S. Kim, G. Ben, Y. Aoki, A. Shikada., Proceedings of the 5th JapanKorea Joint Symposium on Composite Materials, October, 005 4. Y. Aoki and G. Ben., IATSS Review, 9(4), 1, 005 (in Japanese)

5. G. Ben, Y. Furuta, T. Uzawa., Proceedings of the 14th Conference on the JSME Computational Mechanics Division, November, 003 (in Japanese) 6. Yoshio Aoki, Goichi Ben, Yuka Iizuka, Key Engineering Materials, 334335, 197, 006 7. Yoshio Aoki, Goichi Ben, Hyoung Soo Kim, Akihisa Tabata, Proceedings of conference of the Society for the advancement of material and process engineering, L66, May, 008 8. Hyoung Soo Kim, Goichi Ben, Yoshio Aoki, Proceedings of the USJapan Conference on Composite Materials 008, AUT3, June, 008 9. Yoshio Aoki, Goichi Ben, Hyoung Soo Kim, Proceedings of Internetional Crashworthiness Conference 008, 008080, July, 008Reference 1