Blunt Impact Damage Formation on Composite Aircraft Structures

Similar documents
Impact Damage Formation on Composite Aircraft Structures

Disclaimer for FAA Research Publication

Managing Impact Risk for Composite Structures: Unifying Durability and Damage Tolerance Perspective.

(A) 0 (B) mv (C) 2mv (D) 2mv sin θ (E) 2mv cos θ

Designation: NWI DIF Horizontal Hail impact Standard

E 490 FE Exam Prep. Engineering Mechanics

Two Tier projects for students in ME 160 class

Development of Reliability-Based Damage Tolerant Structural Design Methodology

Impulse/Momentum And Its Conservation

Strength and Damage Tolerance of Adhesively Bonded Joints: Stress-Based and Fracture-Based Approaches

A Numerical Study on Prediction of BFS in Composite Structures under Ballistic Impact

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

Conservation of Momentum and Energy

PRINCIPLE OF LINEAR IMPULSE AND MOMENTUM FOR A SYSTEM OF PARTICLES AND CONSERVATION OF LINEAR MOMENTUM FOR A SYSTEM OF PARTICLES

ROBUST DESIGN OF A MECHANICAL LATCH

EF 151 Final Exam, Fall, 2011 Page 1 of 11

Part Two: Earlier Material

Phy211: General Physics I Lab page 1 of 5 PCC-Cascade

PRINCIPLE OF LINEAR IMPULSE AND MOMENTUM AND CONSERVATION OF LINEAR MOMENTUM FOR SYSTEMS OF PARTICLES

Composite Sandwich Structures with Honeycomb Core subjected to Impact

CEE 271: Applied Mechanics II, Dynamics Lecture 17: Ch.15, Sec.2 4

APPLICATIONS. CEE 271: Applied Mechanics II, Dynamics Lecture 17: Ch.15, Sec.4 7. IMPACT (Section 15.4) APPLICATIONS (continued) IMPACT READING QUIZ

EXPERIMENTAL AND NUMERICAL STUDY OF OBLIQUE IMPACT ON HELICOPTER BLADES INFLUENCE OF THE CURVATURE

Center of Mass & Linear Momentum

First Name: Last Name: Section: 1. March 26, 2008 Physics 207 EXAM 2

The following examination-style questions are suitable for assessing evidence of learning in Topic 1. [This diagram is not drawn to scale.

Review Chapter 1 and 2 [184 marks]

Nonlinear Buckling Prediction in ANSYS. August 2009

Part I Review Unit Review Name Momentum and Impulse

Outline. Collisions in 1- and 2-D. Energies from Binary Star Expt. Energy Plot. Energies with Linear Fit. Energy Plot

Impact. m k. Natural Period of Vibration τ. Static load Gray area Impact load t > 3 τ. Absorbing energy. Carrying loads

Stress Analysis on Bulkhead Model for BWB Heavy Lifter Passenger Aircraft

Probabilistic Approaches

Questions on the December Assessment are broken into three categories: (Both MC and FR type questions can be in the following forms):

RELEASED. Go to next page. 2. The graph shows the acceleration of a car over time.

PSI AP Physics I Work and Energy

Impulse simply refers to a change in momentum, and is usually caused by a change in velocity, as described by p = m v.

Physics 12 Final Exam Review Booklet # 1

is acting on a body of mass m = 3.0 kg and changes its velocity from an initial

Using Energy History Data to Obtain Load vs. Deflection Curves from Quasi-Static Abaqus/Explicit Analyses

PARAMETRIC STUDY OF THE STRUCTURAL CAPACITY OF REINFORCED CONCRETE CONTAINMENT SUBJECTED TO LARGE COMMERICAL AIRCRFAT IMPACT

1. A 1,160-kg car traveling initially with a speed of 25.0 m/s in an easterly direction crashes into the rear end of a

SOLUTION. will destroy the integrity of the work and Anis not permitted

LAWS OF MOTION Newtons laws of motion. (i) First law: Law of inertia. Every body continues to be in its state of rest or of uniform motion in a

Momentum & Energy Review Checklist

Elementary Physics October 8, 2007

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

Extra Circular Motion Questions

LINEAR MOMENTUM AND COLLISIONS

Impact Analysis of Frontal Car Bumper using Long Fibre Reinforced Thermoplastics

Model-Based Engineering and Cyber-Physical Systems

Chapter 15 Kinematics of a Particle: Impulse and Momentum. Lecture Notes for Section 15-2~3

AP Physics 1 Momentum and Impulse Practice Test Name

Name SOLUTION Student ID Score Speed of blocks is is decreasing. Part III. [25 points] Two blocks move on a frictionless

Sometimes (like on AP test) you will see the equation like this:

ConcepTest PowerPoints

Improving Yield for High Pin Count Wafer Probing Applications

Chapter 6 - Linear Momemtum and Collisions

Algebra Based Physics

UNIVERSITY OF SASKATCHEWAN GE MECHANICS III FINAL EXAM APRIL 18, 2011 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS

Use a BLOCK letter to answer each question: A, B, C, or D (not lower case such a b or script such as D)

Mechanics Topic B (Momentum) - 1 David Apsley

Twentieth SLAPT Physics Contest Southern Illinois University Edwardsville April 30, Mechanics Test

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

Name: Block: Date: / / IP 614 Momentum Conservation Lab-Collision

Exam 3--PHYS 101--F15

Announcements. There will still be a WebAssign due this Friday, the last before the midterm.

1 Forces. 2 Energy & Work. GS 104, Exam II Review

Kul Aircraft Structural Design (4 cr) Fatigue Analyses

p p I p p p I p I p p

Curriculum Map-- Kings School District- Physics

Probabilistic Assessment of a Stiffened Carbon Fibre Composite Panel Operating in its Postbuckled Region Authors: Correspondence:

Projectile Motion Horizontal Distance

Momentum is a property of moving matter. Momentum describes the tendency of objects to keep going in the same direction with the same speed.

2. (a) Explain different types of wing structures. (b) Explain the advantages and disadvantages of different materials used for aircraft

The graph shows how an external force applied to an object of mass 2.0 kg varies with time. The object is initially at rest.

Comparison of LS-DYNA and NISA in Solving Dynamic Pulse Buckling Problems in Laminated Composite Beams

A Change in Momentum. Demonstrating the conservation of momentum.

Physics 1501 Lecture 17

Exam II. Spring 2004 Serway & Jewett, Chapters Fill in the bubble for the correct answer on the answer sheet. next to the number.

Curriculum Map-- Kings School District Honors Physics

AERODROMES PANEL (AP) VISUAL AIDS WORKING GROUP (VAWG) FIFTH MEETING. Montréal, Canada, 25 to 27 June 2008

Per 9 10 Momentum_Presentation.notebook. January 20, Momentum.

Parametric Studies of Low Velocity Impact on E-glass/Epoxy using Ls-Dyna

PHYSICS 218 EXAM 2 Tuesday, October 26, 2010

Momentum Practice Problems

PHYSICS 231 INTRODUCTORY PHYSICS I

Momentum and Impulse Concept Tests

1. The diagram below shows the variation with time t of the velocity v of an object.

Impact and Fracture Mechanics Assessment of a Fused Silica Window

PHYS 154 Practice Test 3 Spring 2018

COLLEGE OF FOUNDATION AND GENERAL STUDIES PUTRAJAYA CAMPUS FINAL EXAMINATION TRIMESTER I 2012/2013

LECTURE 13- PROBLEMS. Chapter 1-9,13 Professor Noronha-Hostler Professor Montalvo

Academic Physics! Work and Momentum Summary! Name

Practice Test for Midterm Exam

A study on wheel force measurement using strain gauge equipped wheels

Contact pressure distribution in joints formed by V-band clamps Simon M Barrans 1,a, Goodarz Khodabakhshi 1,b and Qiang Xu 1,c

A Repeated Dynamic Impact Analysis for 7x7 Spacer Grids by using ABAQUS/ Standard and Explicit

N - W = 0. + F = m a ; N = W. Fs = 0.7W r. Ans. r = 9.32 m

PRINCIPLE OF LINEAR IMPULSE AND MOMENTUM

Transcription:

Blunt Impact Damage Formation on Composite Aircraft Structures FAA Supported Research Project Hyonny Kim, Associate Professor University of California, San Diego Hyonny Kim, UCSD 1

incidental contact / collisions with ground vehicles and equipment can result in blunt impact damage a significant source of damage to commercial aircraft extensive sub-surface damage is often not externally visible new all-composite fuselage transport aircraft are coming into service significantly more composite skin surface area is exposed to ground vehicles and equipment Background Basic tools are needed for characterizing blunt impact events to aid in prediction of damage formation and its effect on structural performance. Airbus A330 Horizontal Stabilizer Source: S. Waite, FAA Damage Tolerance and Maintenance Workshop, Chicago, July 19-21, 2006 Hyonny Kim, UCSD 2

Objectives Newly-starting research project has three objectives: identify blunt impact scenarios that are commonly occurring of major concern to airline maintenance organizations and aircraft manufacturers develop methodology for blunt impact threat characterization and modeling experimental identification of key phenomena and parameters governing blunt impact damage formation Expected outcome upon accomplishment of these objectives (along with proper dissemination of the results) are simple modeling tools that can aid in: assessing whether incident could have caused damage improving the resistance and tolerance of composite aircraft structures to damage from blunt impacts establishing practical design criteria for damage resistance and tolerance Hyonny Kim, UCSD 3

Plan of Work Three tasks addressing each of the objectives: Task 1. Identification of Common Blunt Impact Scenarios and Establishment of Partner Relationships initial phase of the project active communication with the following organizations airlines aircraft manufacturers OEM suppliers FAA and EASA identify blunt impact scenarios that are: most commonly occurring of major concern e.g., source of impact, location, and type of damage formed visits to maintenance depots document findings via photographs (if permitted) and written notes outcome of this task will be summarized in a written document/faa report Hyonny Kim, UCSD 4

Task 2. Methodology for Impact Threat Characterization focus on defining generalized methodology for describing, characterizing, and establishing blunt impact threats e.g., assign all aircraft-adjacent ground items with threat index problem is very complex due to many variables that are important impactor can be different types of ground vehicles or equipment, and various locations on these equipment (e.g., corner, long edge, or flat face) target can be the many locations of the aircraft exposed to contact with ground vehicle/equipment fuselage, nacelles, wing skins, control surfaces, etc. impacts can be at or near internal stiffeners, or away from them, thereby greatly affecting the local stiffness of the structure local stiffness governs the contact forces generated during impact incidence angle between impactor and composite panel surface plays major role in nature of contact force history How are F 1 and F 2 related to contact angle? Hyonny Kim, UCSD 5

Task 2 Contd. Consider: 5,500 kg vehicle at 0.5 m/s (1.1 mph) collides with an aircraft. How do we describe this event? If kinetic energy based, is the entire kinetic energy of 687 J (507 ft-lb) involved in the projectile-target interaction? How much of the aircraft structure is active in reacting to the collision? Alternatively, peak contact force has been shown to be a descriptive parameter for describing damage initiation (Jackson & Poe 1992, Schoeppner & Abrate 2000) BIG ISSUE: contact force is not readily known quantity need to estimate using blunt impact characterization model First-Order Blunt Impact Characterization Model - similar momentum-force based tools exist for automobile impact simulation SMAC, EDSMAC, CRASH, etc. k 1 and c 1 represent impactor-target interaction k 1 could be nonlinear k 2 and c 2 represent stiffness supporting target, or portion of target active in the problem represented by mass M 2 applied force F(t) allows additional forcing input after initial contact parameters measured experimentally or determined via numerical simulation expected outcome: user-friendly analytical tool that estimates peak contact force and max deflection assess whether failure threshold has been exceeded if so, what type of inspection is needed validation needed experiments proposed in Task 3 Hyonny Kim, UCSD 6

Task 2 Preliminary Results: Interrogation of Contact Angle Effect via FEA High-mass projectile 500 kg (1103 lb) 127 mm (5 in.) corner radius initial velocity 0.447 m/s (1.0 mph) to right; KE = 50 J no applied external force constrained to exhibit only horizontal motion Curved Composite Panel 6.35 mm thickness (0.25 in.) radius of curvature 3 m clamped b.c. at top and bottom oriented at 45 o and 10 o angle w.r.t. ground plane FE simulation conducted in ABAQUS/Explicit Contact Force History? Same Panel and B.C. s Rotated to Adjust Angle of Contact w.r.t. Direction of Motion Hyonny Kim, UCSD 7

High Mass Impact Simulation Results Observations: for lower contact angle, Increased contact duration for 45 : 94 ms for 10 : 376 ms Contact spread across more elongated area 45 Relative to Horizontal Additional model info: linear material behavior i.e., damage not modeled geometric nonlinearity contact defined as frictionless 10 Relative to Horizontal Hyonny Kim, UCSD 8

Contact Force History Total contact force vector sum of x- and y-direction force components acts in direction normal to panel surface (frictionless contact defined) peak force NOT dependent on panel orientation panel target has identical stiffness thus same maximum displacement (quasi-static like event) Longer duration pulse can be more damaging Notes: same peak force only for projectile having constrained motion if x-direction force applied behind projectile, peak contact forces will NOT be the same (trigonometric balance) e.g., driving force from vehicle wheel torque X-Direction Component Contact Force (N) Total Contact Force (lbf) 8000 6000 4000 2000 0 10000 8000 6000 4000 2000 0 x-direction force component 0 0.1 0.2 0.3 0.4 0.5 Time (sec) total force Panel Angle 45 10 Mesh-Related Numerical Artifacts Panel Angle 45 10 0 0.1 0.2 0.3 0.4 0.5 Hyonny Kim, UCSD 9 Time (sec)

Momentum Transfer Momentum of projectile imparts impulse to structure during impact event projectile initial momentum is 500 kg x 0.447 m/s = 223.5 kg-m/s (or N-s) total momentum change is 2X due to projectile bouncing off target and returning with equal but opposite velocity: 447 N-s Total impulse on structure computed by integration of total force over time (area under f vs. t curve) dependent on panel orientation for 45 : 2,480 N-s for 10 : 623 N-s (4X higher than 45 ) acts normal to panel surface Impulse found to scale by trigonometric relationship p θ mvo = 2 Sinθ Normal-Direction Momentum Transfer (N-s) 2500 2000 1500 1000 where θ is angle between panel surface 500 and direction of projectile motion (constrained) good match-up with FEA 0 500 0 Total Impulse Imparted (N-s) Panel Angle 45 10 0 0.1 0.2 0.3 0.4 0.5 Time (sec) 3000 2500 2000 1500 1000 447 N-s Normal Impact Rebound 0 10 20 30 40 50 60 70 80 90 Surface Angle Relative to Horiz. Hyonny Kim, UCSD 10 FEA Data ScalingEqn

Critical Damage Threshold Force Concept Initiation of damage can be described by damage threshold load (DTL)* can be determined from contact force vs. time and force vs. displacement plots DTL shown to scale with panel thickness as t 3/2 : (relationship by Davies and Zhang**) Impact Contact Force Histories* Damage Threshold Load* Refs: * ** Hyonny Kim, UCSD 11

Task 3. Key Phenomena and Parameters Governing Blunt Impact Damage experimentally identify key phenomena and parameters governing blunt impact damage correlation of global parameters to local parameters measured by test and predicted by the previously described simple models global parameters: ground vehicle mass, velocity, aircraft substructure mass local parameter: contact stiffness, peak contact force and deflection during impact event damage extent resulting from various impact force levels on a given component can be catalogued and then looked up when an incident occurs in service can aid in making rapid maintenance/repair decisions Example of userfriendly lookup table for complex process of high velocity ice impact, by H. Kim et al. 2003. Panel thickness vs ice ball diameter H/D is found to be key parameter describing damage initiation. α 16 14 12 10 8 6 4 2 α = ρ ice ( τ 2 ILS V ρ = 870 kg ice m 3 2 FTE G ILS τ ILS = 100 MPa G ILS = 5. 17 GPa ) FTE for Ice Impact onto Woven Carbon/Epoxy Fit to Data Safe Design Curve Panel Type 8552, 1.22 mm 8552, 1.82 mm 8552, 2.44 mm 977, 1.42 mm 977, 1.91 mm 977, 2.62 mm 0 0x10 0 Hyonny Kim, UCSD 0.00 0.02 0.04 0.06 0.08 0.10 0.12 12 H/D 3x10 4 2x10 4 1x10 4 2 V FTE (m/s) 2

Task 3. Contd. Planned Experiments lab-scale panel tests quasi-static indentation and low velocity impact panels having varying boundary support stiffness measure initial stiffness w/out damage measure critical force for damage initiation as function of indenter radius, contact stiffness, boundary support stiffness, etc. generalization of results to encompass wide range of parameters full-scale blunt impact tests impact tests on full-scale structures by actual ground vehicles/equipment tests conducted using unique large high-rate equipment at UCSD or on test sled validate models in Task 2 Hyonny Kim, UCSD 13

Industry/Agency Participation key component of this research activity is participation from industry/agency partners research focus and activities should be relevant to the user community requested industry/agency participation broken down into three levels: Level 1. Initial Guidance. Contribution to Task 1. Hyonny Kim will communicate and/or visit the facility to learn and see first hand about the companies experiences in this topic and document experiences. Level 2. Advisory. Provide initial guidance, as well as ongoing advisement, particularly in Task 2 methodology development, and in defining parameters for experiments to be conducted in Task 3. Level 3. Provide material support for Task 3. Make available test panels, substructure components, and/or impactors (e.g., a cargo cart) for full-scale blunt impact investigations. Hyonny Kim, UCSD 14

Request From Today s Audience feedback on proposed activities wish list from industry/agency perspective what should analytical tool look like? be capable of? what quantities/outputs are most important to you? willingness to participate in survey to be issued querying about damage types and their sources, etc. additional participants contact info for Hyonny Kim: hyonny@ucsd.edu Hyonny Kim, UCSD 15