SEARCH · Engineering Papers

Results for “Failure mode and effect analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Experimental observations and finite element analysis of the initiation of fiber microbuckling in notched composite laminates

An understanding was developed of the factors that determine the semi-circular edge-notched compressive strength and the associated failure mode(s) were identified of thermoplastic composite laminates with multidirectional stacking sequences. The experimental observations and the detailed literature review suggest at least four factors that affected the determination of the strain levels at which fiber microbuckling initiates and thus, partially control the composite's compression strength. The dependent variables studied are the compressive strength of a reduced gage section compression specimen and the compression strength of a compression specimen with two semi-circular edge notches (no opposite free edges) centered along the gage section. In this research, specimens containing two semi-circular edge notches (no opposite free edges) were loaded in compression at a relatively slow rate to provide more stable development of fiber microbuckling damage. The results indicate that the local constraints (free surfaces, supporting ply orientation, and resin-rich regions) significantly affect the strain level for the initiation of in-plane fiber microbuckling. Preliminary results at an elevated temperature, 77 C, showed the shear stress yield strength of the resin was reduced and consequently, the resistance to fiber microbuckling was also reduced. The finite element analysis of the perfectly straight fiber problem indicates that the free surface effect causes a 10 percent reduction in the critical buckling strain. However, the experimentally measured reduction for fibers with an initial fiber curvature, was 35 percent.

Guynn, E. Gail↗

NASA NTRS2020

Stochastic Simulation of Mudcrack Damage Formation in an Environmental Barrier Coating

The FEAMAC/CARES program, which integrates finite element analysis (FEA) with the MAC/GMC (Micromechanics Analysis Code with Generalized Method of Cells) and the CARES/Life (Ceramics Analysis and Reliability Evaluation of Structures / Life Prediction) programs, was used to simulate the formation of mudcracks during the cooling of a multilayered environmental barrier coating (EBC) deposited on a silicon carbide substrate. FEAMAC/CARES combines the MAC/GMC multiscale micromechanics analysis capability (primarily developed for composite materials) with the CARES/Life probabilistic multiaxial failure criteria (developed for brittle ceramic materials) and Abaqus (Dassault Systèmes) FEA. In this report, elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio elements enabled the formation of crack boundaries, leading to development of mudcrack-patterned damage. Finite element models of a disk-shaped three-dimensional specimen and a twodimensional model of a through-the-thickness cross section subjected to progressive cooling from 1,300 °C to an ambient temperature of 23 °C were made. Mudcrack damage in the coating resulted from the buildup of residual tensile stresses between the individual material constituents because of thermal expansion mismatches between coating layers and the substrate. A two-parameter Weibull distribution characterized the coating layer stochastic strength response and allowed the effect of the Weibull modulus on the formation of damage and crack segmentation lengths to be studied. The spontaneous initiation of cracking and crack coalescence resulted in progressively smaller mudcrack cells as cooling progressed, consistent with a fractal-behaved fracture pattern. Other failure modes such as delamination, and possibly spallation, could also be reproduced. The physical basis assumed and the heuristic approach employed, which involves a simple stochastic cellular automaton methodology to approximate the crack growth process, are described. The results ultimately show that a selforganizing mudcrack formation can derive from a Weibull distribution that is used to describe the stochastic strength response of the bulk brittle ceramic material layers of an EBC.

Nemeth, Noel N.↗

NASA NTRS2005

Hubble Space Telescope Star Tracker ad Two-Gyro Control Law Design, Implementation, and On-Orbit Performance

The Hubble Space Telescope (HST) normally requires three gyroscopes for three-axis rate control. The loss of the Space Shuttle Columbia on STS-107 resulted in the cancellation of a shuttle-based HST Servicing Mission 4. Therefore, HST must operate using the on-board hardware until an alternate means of servicing can be accomplished. The probability of gyro failure indicates that fewer than three gyros will be operable before any servicing mission can be performe& To mitigate this, and to extend the HST life expectancy, a rate estimation and control algorithm was developed that requires two gyros to measure rate about two axes, with the remaining axis rate estimated using one of three alternate sensors. Three-axis magnetometers (MSS) are used for coarse rate estimation during large maneuvers and during occultations of other sensors. Fixed-Head Star Trackers (FHSTs) are used for rate estimation during safe mode recovery and during transition to science operations. Fine rate estimation during science operations is performed using the Fine Guidance Sensors (FGSs). The FHST mode (T2G) relies on star vectors as measured by the FHSTs to estimate vehicle rate about the axis not measured by the gyros. Since the FHSTs were not designed to estimate body rate, this method involves a unique set of problems that had to be overcome in the final design, such as the effect of FHST break tracks and moving targets on rate estimation. The solutions to these problems, as well as a detailed description of the design and implementation of the rate estimation are presented Also included are the time domain and frequency domain analysis of the T2G control law. A high fidelity HST simulator (HSTSIM) was used to verify T2G performance prior to on-orbit use. Results of these simulations are also presented. Finally, analysis of actual T2G on-orbit test results is presented for design validation.

VanArsdall, John C.↗

NASA NTRS2015

Face Sheet/Core Disbond Growth in Honeycomb Sandwich Panels Subjected to Ground-Air-Ground Pressurization and In-Plane Loading

Typical damage modes in light honeycomb sandwich structures include face sheet/core disbonding and core fracture, both of which can pose a threat to the structural integrity of a component. These damage modes are of particular interest to aviation certification authorities since several in-service occurrences, such as rudder structural failure and other control surface malfunctions, have been attributed to face sheet/core disbonding. Extensive studies have shown that face sheet/core disbonding and core fracture can lead to damage propagation caused by internal pressure changes in the core. The increasing use of composite sandwich construction in aircraft applications makes it vitally important to understand the effect of ground-air-ground (GAG) cycles and conditions such as maneuver and gust loads on face sheet/core disbonding. The objective of the present study was to use a fracture mechanics based approach developed earlier to evaluate the loading at the disbond front caused by ground-air-ground pressurization and in-plane loading. A honeycomb sandwich panel containing a circular disbond at one face sheet/core interface was modeled with three-dimensional (3D) solid finite elements. The disbond was modeled as a discrete discontinuity and the strain energy release rate along the disbond front was computed using the Virtual Crack Closure Technique (VCCT). Special attention was paid to the pressure-deformation coupling which can decrease the pressure load within the disbonded sandwich section significantly when the structure is highly deformed. The commercial finite element analysis software, Abaqus/Standard, was used for the analyses. The recursive pressure-deformation coupling problem was solved by representing the entrapped air in the honeycomb cells as filled cavities in Abaqus/Standard. The results show that disbond size, face sheet thickness and core thickness are important parameters that determine crack tip loading at the disbond front. Further, the pressure-deformation coupling was found to have an important load decreasing effect [6]. In this paper, a detailed problem description is provided first. Second, the analysis methodology is presented. The fracture mechanics approach used is described and the specifics of the finite element model, including the fluid-filled cavities, are introduced. Third, the initial model verification and validation are discussed. Fourth, the findings from a closely related earlier study [6] are summarized. These findings provided the basis for the current investigation. Fifth, an aircraft ascent scenario from 0 to 12192 m (0 to 40000 ft) is considered and the resulting crack tip loading at the disbond front is determined. In-plane loading to simulate maneuvers and gust conditions are also considered. Sixth, the results are shown for a curved panel, which was used to simulate potential fuselage applications. Finally, a brief summary of observations is presented and recommendations for improvement are provided.

Chen, Zhi M.↗

NASA NTRS1961

A Study of Longitudinal Control Problems at Low and Negative Damping and Stability with Emphasis on Effects of Motion Cues

As part of a general investigation to determine the effects of simulator motions on pilot opinion and task performance over a wide range of vehicle longitudinal dynamics, a cooperative NASA-AMAL program was conducted on the centrifuge at Johnsville, Pennsylvania. The test parameters and measurements for this program duplicated those of earlier studies made at Ames Research Center with a variable-stability airplane and with a pitch-roll chair flight simulator. Particular emphasis was placed on the minimum basic damping and stability the pilots would accept and on the minimum dynamics they considered controllable in the event of stability-augmentation system failure. Results of the centrifuge-simulator program indicated that small positive damping was required by the pilots over most of the frequency range covered for configurations rated acceptable for emergency conditions only (e.g., failure of a pitch damper). It was shown that the pilot's tolerance for unstable dynamics was dependent primarily on the value of damping. For configurations rated acceptable for emergency operation only, the allowable instability and damping corresponded to a divergence time to double amplitude of about 1 second. Comparisons were made of centrifuge, pitch-chair and fixed-cockpit simulator tests with flight tests. Pilot ratings indicated that the effects of incomplete or spurious motion cues provided by these three modes of simulation were important only for high-frequency, lightly damped dynamics or unstable, moderately damped dynamics. The pitch- chair simulation, which provided accurate angular-acceleration cues to the pilot, compared most favorably with flight. For the centrifuge simulation, which furnished accurate normal accelerations but spurious pitching and longitudinal accelerations, there was a deterioration of pilots' opinion relative to flight results. Results of simulator studies with an analog pilot replacing the human pilot illustrated the adaptive capability of human pilots in coping with the wide range of vehicle dynamics and the control problems covered in this study. It was shown that pilot-response characteristics, deduced by the analog-pilot method, could be related to pilot opinion. Possible application of these results for predicting flight-control problems was illustrated by means of an example control-problem analysis. The results of a brief evaluation of a pencil-type side-arm controller in the centrifuge showed a considerable improvement in the pilots' ability to cope with high-frequency, low-damping dynamics, compared to results obtained with the center stick. This improvement with the pencil controller was attributed primarily to a marked reduction in the adverse effects of large and exaggerated pitching and longitudinal accelerations on pilot control precision.

Sadoff, Melvin↗

NASA NTRS2016

A Generic Modeling Process to Support Functional Fault Model Development

Functional fault models (FFMs) are qualitative representations of a system's failure space that are used to provide a diagnostic of the modeled system. An FFM simulates the failure effect propagation paths within a system between failure modes and observation points. These models contain a significant amount of information about the system including the design, operation and off nominal behavior. The development and verification of the models can be costly in both time and resources. In addition, models depicting similar components can be distinct, both in appearance and function, when created individually, because there are numerous ways of representing the failure space within each component. Generic application of FFMs has the advantages of software code reuse: reduction of time and resources in both development and verification, and a standard set of component models from which future system models can be generated with common appearance and diagnostic performance. This paper outlines the motivation to develop a generic modeling process for FFMs at the component level and the effort to implement that process through modeling conventions and a software tool. The implementation of this generic modeling process within a fault isolation demonstration for NASA's Advanced Ground System Maintenance (AGSM) Integrated Health Management (IHM) project is presented and the impact discussed.

Design Analysis↗

NASA NTRS1999

Lamb wave characterization of the effects of long-term thermal-mechanical aging on composite stiffness

Lamb waves offer a promising method of evaluating damage in composite materials. The Lamb wave velocity is directly related to the material parameters, so an effective tool exists to monitor damage in composites by measuring the velocity of these waves. The Lamb Wave Imager (LWI) uses a pulse/receive technique that excites an antisymmetric Lamb mode and measures the time-of-flight over a wide frequency range. Given the material density and plate thickness, the bending and out-of-plane shear stiffnesses are calculated from a reconstruction of the dispersion curve. In this study, the time-of-flight as well as the elastic stiffnesses D11, D22, A44, and A55 for composite samples which have undergone combined thermal and mechanical aging are obtained. The samples examined include a baseline specimen with 0 cycles, specimens which have been aged 2350 and 3530 cycles at high strain levels, and one specimen aged 3530 cycles at low strain levels.

Aircraft↗

NASA NTRS2007

Analysis of Composite Panel-Stiffener Debonding Using a Shell/3D Modeling Technique

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used successfully primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities, however, requires the successful demonstration of the methodology on the structural level. For this purpose, a panel was selected that is reinforced with stiffeners. Shear loading causes the panel to buckle, and the resulting out-of-plane deformations initiate skin/stiffener separation at the location of an embedded defect. A small section of the stiffener foot, web and noodle as well as the panel skin in the vicinity of the delamination front were modeled with a local 3D solid model. Across the width of the stiffener foot, the mixedmode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with a mixed-mode failure criterion of the graphite/epoxy material. Computed failure indices were compared to corresponding results where the entire web was modeled with shell elements and only a small section of the stiffener foot and panel were modeled locally with solid elements. Including the stiffener web in the local 3D solid model increased the computed failure index. Further including the noodle and transition radius in the local 3D solid model changed the local distribution across the width. The magnitude of the failure index decreased with increasing transition radius and noodle area. For the transition radii modeled, the material properties used for the noodle area had a negligible effect on the results. The results of this study are intended to be used as a guide for conducting finite element and fracture mechanics analyses of delamination and debonding in complex structures such as integrally stiffened panels.

Krueger, Ronald↗

NASA NTRS2008

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗

NASA NTRS2010

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Allgood, Daniel C.↗

NASA NTRS1993

Numerical and experimental investigation of the bending response of thin-walled composite cylinders

A numerical and experimental investigation of the bending behavior of six eight-ply graphite-epoxy circular cylinders is presented. Bending is induced by applying a known end-rotation to each end of the cylinders, analogous to a beam in bending. The cylinders have a nominal radius of 6 inches, a length-to-radius ratio of 2 and 5, and a radius-to-thickness ratio of approximately 160. A (+/- 45/0/90)S quasi-isotropic layup and two orthotropic layups, (+/- 45/0 sub 2)S and (+/- 45/90 sub 2)S, are studied. A geometrically nonlinear special-purpose analysis, based on Donnell's nonlinear shell equations, is developed to study the prebuckling responses and gain insight into the effects of non-ideal boundary conditions and initial geometric imperfections. A geometrically nonlinear finite element analysis is utilized to compare with the prebuckling solutions of the special-purpose analysis and to study the buckling and post buckling responses of both geometrically perfect and imperfect cylinders. The imperfect cylinder geometries are represented by an analytical approximation of the measured shape imperfections. Extensive experimental data are obtained from quasi-static tests of the cylinders using a test fixture specifically designed for the present investigation. A description of the test fixture is included. The experimental data are compared to predictions for both perfect and imperfect cylinder geometries. Prebuckling results are presented in the form of displacement and strain profiles. Buckling end-rotations, moments, and strains are reported, and predicted mode shapes are presented. Observed and predicted moment vs. end-rotation relations, deflection patterns, and strain profiles are illustrated for the post buckling responses. It is found that a geometrically nonlinear boundary layer behavior characterizes the prebuckling responses. The boundary layer behavior is sensitive to laminate orthotropy, cylinder geometry, initial geometric imperfections, applied end-rotation, and non-ideal boundary conditions. Buckling end-rotations, strains, and moments are influenced by laminate orthotropy and initial geometric imperfections. Measured buckling results correlate well with predictions for the geometrically imperfect specimens. The postbuckling analyses predict equilibrium paths with a number of scallop-shaped branches that correspond to unique deflection patterns. The observed postbuckling deflection patterns and measured strain profiles show striking similarities to the predictions in some cases. Ultimate failure of the cylinders is attributed to an interlaminar shear failure mode along the nodal lines of the postbuckling deflection patterns.

Fuchs, J. P.↗

NASA NTRS2011

Max Launch Abort System (MLAS) Pad Abort Test Vehicle (PATV) II Attitude Control System (ACS) Integration and Pressurization Subsystem Dynamic Random Vibration Analysis

In order to mitigate catastrophic failures on future generation space vehicles, engineers at the National Aeronautics and Space Administration have begun to integrate a novel crew abort systems that could pull a crew module away in case of an emergency at the launch pad or during ascent. The Max Launch Abort System (MLAS) is a recent test vehicle that was designed as an alternative to the baseline Orion Launch Abort System (LAS) to demonstrate the performance of a "tower-less" LAS configuration under abort conditions. The MLAS II test vehicle will execute a propulsive coast stabilization maneuver during abort to control the vehicles trajectory and thrust. To accomplish this, the spacecraft will integrate an Attitude Control System (ACS) with eight hypergolic monomethyl hydrazine liquid propulsion engines that are capable of operating in a quick pulsing mode. Two main elements of the ACS include a propellant distribution subsystem and a pressurization subsystem to regulate the flow of pressurized gas to the propellant tanks and the engines. The CAD assembly of the Attitude Control System (ACS) was configured and integrated into the Launch Abort Vehicle (LAV) design. A dynamic random vibration analysis was conducted on the Main Propulsion System (MPS) helium pressurization panels to assess the response of the panel and its components under increased gravitational acceleration loads during flight. The results indicated that the panels fundamental and natural frequencies were farther from the maximum Acceleration Spectral Density (ASD) vibrations which were in the range of 150-300 Hz. These values will direct how the components will be packaged in the vehicle to reduce the effects high gravitational loads.

Ekrami, Yasamin↗

NASA NTRS2019

Demisability of GFRP and CFRP Components of Reentering Orbital Debris: Phase I Test Results

Observations of surviving reentry debris on the ground and research performed by Hyperschall Technologie Göttingen (HTG) [1] indicated that significantly more glass fiber-reinforced polymer (GFRP) and carbon fiber-reinforced polymer (CFRP) components survive reentry than current models predict. NASA’s Orbital Debris Program Office conducted a series of tests to evaluate the accuracy of material demise models for reentering orbital debris used in NASA’s Object Reentry Survival Analysis Tool (ORSAT) and Debris Assessment Software (DAS). Testing is planned in a multi-phase series to allow for quick quantification of results as well as refinement of methods resulting from lessons learned during early phases. The Phase 1 tests discussed here validated ORSAT models for homogeneous metals, provided an efficient quantification of composite material demisability properties like mass loss rate and overall time to demise, and identified potential failure modes, which are currently not well understood. Phase 2 tests will be used to further understand mass loss rates and modes of both thermal and mechanical failure in composite materials. The authors exposed 95 samples of aluminum, CFRP, Kevlar fiber-reinforced polymer, GFRP, and sheets of G10 fiberglass to conditions approximating the reentry environment using an inductively coupled plasma (ICP) torch facility. The cylindrical CFRP samples were exposed to the atmospheric pressure plasma, at both the end and the midpoint, to investigate the difference in demisability between parts with exposed edges, like panels, and parts with no edges, such as carbonoverwrapped pressure vessels (COPVs). In a non-oxidative environment, no composite materials demised within the 5-minute test time. In the oxidative, elevated heat flux environment, CFRP samples demised between 210 s and 270 s. For the first 100 s of insertion time, most of the mass loss was due to pyrolysis of resin, creating an approximately bi-linear mass-loss rate curve with time. In a non-oxidative environment, carbon filaments were observed to unravel from some of the CFRP end-burned samples; however, this effect did not seem to affect the overall time to demise for the samples significantly. These results indicate that both GFRP and CFRP components survive reentry with significantly more remaining mass than current models predict.

Greene, Benton R.↗

NASA NTRS2006

Nonlinear Spring Finite Elements for Predicting Mode I-Dominated Delamination Growth in Laminated Structure with Through-Thickness reinforcement

One particular concern of polymer matrix composite laminates is the relatively low resistance to delamination cracking, in particular when the dominant type of failure is mode I opening. One method proposed for alleviating this problem involves the insertion pultruded carbon pins through the laminate thickness. The pins, known as z-pins, are inserted into the prepreg laminate using an ultrasonic hammer prior to the curing process, resulting in a field of pins embedded normal to the laminate plane as illustrated in Figure. 1. Pin diameters range between 0.28-mm to 0.5-mm and standard areal densities range from 0.5% to 4%. The z-pins are provided by the manufacturer, Aztex(Registered TradeMark) , in a low-density foam preform, which acts to stabilize orientation of the pins during the insertion process [1-3]. Typical pin materials include boron and carbon fibers embedded in a polymer matrix. A number of methods have been developed for predicting delamination growth in laminates reinforced with z-pins. During a study on the effect of z-pin reinforcement on mode I delamination resistance, finite element analyses of z-pin reinforced double cantilever beam (DCB) specimens were performed by Cartie and Partridge [4]. The z-pin bridging stresses were modeled by applying equivalent forces at the pin locations. Single z-pin pull-out tests were performed to characterize the traction law of the pins under mode I loading conditions. Analytical solutions for delamination growth in z-pin reinforced DCB specimens were independently derived by Robinson and Das [5] and Ratcliffe and O'Brien [6]. In the former case, pin bridging stresses were modeled using a distributed load and in the latter example the bridging stresses were discretely modeled by way of grounded springs. Additionally, Robinson and Das developed a data reduction strategy for calculating mode I fracture toughness, G(sub Ic), from a z-pin reinforced DCB specimen test [5]. In both cases a traction law similar to that adopted by Cartie and Partridge was used to represent z-pin failure under mode I loading conditions. In the current work spring elements available in most commercial finite element codes were used to model z-pins. The traction law used in previous analyses [4-6] was employed to represent z-pin damage. This method is intended for and is limited to simulating z-pins in composite laminate structure containing mode I-dominated delamination cracking. The current technique differs from previous analyses in that spring finite elements (available in commercial codes) are employed for simulating zpins, reducing the complexity of the analysis construction process. Furthermore, the analysis method can be applied to general structure that experiences mode I-dominated delamination cracking, in contrast to existing analytical solutions that are only applicable to coupon DCB specimens.

Ratcliffe, James G.↗

NASA NTRS2012

In-Situ Nondestructive Evaluation of Kevlar(Registered Trademark)and Carbon Fiber Reinforced Composite Micromechanics for Improved Composite Overwrapped Pressure Vessel Health Monitoring

NASA has been faced with recertification and life extension issues for epoxy-impregnated Kevlar 49 (K/Ep) and carbon (C/Ep) composite overwrapped pressure vessels (COPVs) used in various systems on the Space Shuttle and International Space Station, respectively. Each COPV has varying criticality, damage and repair histories, time at pressure, and pressure cycles. COPVs are of particular concern due to the insidious and catastrophic burst-before-leak failure mode caused by stress rupture (SR) of the composite overwrap. SR life has been defined [1] as the minimum time during which the composite maintains structural integrity considering the combined effects of stress level(s), time at stress level(s), and associated environment. SR has none of the features of predictability associated with metal pressure vessels, such as crack geometry, growth rate and size, or other features that lend themselves to nondestructive evaluation (NDE). In essence, the variability or surprise factor associated with SR cannot be eliminated. C/Ep COPVs are also susceptible to impact damage that can lead to reduced burst pressure even when the amount of damage to the COPV is below the visual detection threshold [2], thus necessitating implementation of a mechanical damage control plan [1]. Last, COPVs can also fail prematurely due to material or design noncompliance. In each case (SR, impact or noncompliance), out-of-family behavior is expected leading to a higher probability of failure at a given stress, hence, greater uncertainty in performance. For these reasons, NASA has been actively engaged in research to develop NDE methods that can be used during post-manufacture qualification, in-service inspection, and in-situ structural health monitoring. Acoustic emission (AE) is one of the more promising NDE techniques for detecting and monitoring, in real-time, the strain energy release and corresponding stress-wave propagation produced by actively growing flaws and defects in composite materials [3,4,5,6,7,8]. To gain further insight into the mechanisms responsible for composite rupture, broadband modal acoustic emission analysis was used. Also, since AE data reduction proved to be very time consuming, specialized data reduction software was written to automate the process.

Waller, Jess↗

NASA NTRS1994

Thin tailored composite wing for civil tiltrotor

The tiltrotor aircraft is a flight vehicle which combines the efficient low speed (i.e., take-off, landing, and hover) characteristics of a helicopter with the efficient cruise speed of a turboprop airplane. A well-known example of such vehicle is the Bell-Boeing V-22 Osprey. The high cruise speed and range constraints placed on the civil tiltrotor require a relatively thin wing to increase the drag-divergence Mach number which translates into lower compressibility drag. It is required to reduce the wing maximum thickness-to-chord ratio t/c from 23% (i.e., V-22 wing) to 18%. While a reduction in wing thickness results in improved aerodynamic efficiency, it has an adverse effect on the wing structure and it tends to reduce structural stiffness. If ignored, the reduction in wing stiffness leads to susceptibility to aeroelastic and dynamic instabilities which may consequently cause a catastrophic failure. By taking advantage of the directional stiffness characteristics of composite materials the wing structure may be tailored to have the necessary stiffness, at a lower thickness, while keeping the weight low. The goal of this study is to design a wing structure for minimum weight subject to structural, dynamic and aeroelastic constraints. The structural constraints are in terms of strength and buckling allowables. The dynamic constraints are in terms of wing natural frequencies in vertical and horizontal bending and torsion. The aeroelastic constraints are in terms of frequency placement of the wing structure relative to those of the rotor system. The wing-rotor-pylon aeroelastic and dynamic interactions are limited in this design study by holding the cruise speed, rotor-pylon system, and wing geometric attributes fixed. To assure that the wing-rotor stability margins are maintained a more rigorous analysis based on a detailed model of the rotor system will need to ensue following the design study. The skin-stringer-rib type architecture is used for the wing-box structure. The design variables include upper and lower skin ply thicknesses and orientation angles, spar and rib web thicknesses and cap areas, and stringer cross-sectional areas. These design variables will allow the maximum tailoring of the structure to meet the design requirements most efficiently. Initial dynamic analysis has been conducted using MSC/NASTRAN to determine the baseline wing's frequencies and mode shapes. For the design study we intend to use the finite-element based code called WIDOWAC (Wing Design Optimization With Aeroeastic Constraints) that was developed at NASA Langley in early 1970's for airplane wing structural analysis and preliminary design. Currently, the focus is on modification and validation of this code which will be used for the civil tiltrotor design efforts.

Rais-Rohani, Masoud↗

NASA NTRS2020

Evaluation of Fatigue Damage Accumulation Functions for Delamination Initiation and Propagation

The present report follows on the cohesive fatigue damage model methodology proposed in NASA-TP-2018-219838. In that report, an empirical function describing the incremental damage due to cyclic loading was used to calculate fatigue damage within a cohesive formulation. The form of the function was developed such that, when integrated at a constant stress amplitude from no damage to failure, it produces a life versus load response that is consistent with an S-N diagram. Therefore, the parameters of the damage model could be obtained by fitting the model predictions to an S-N diagram. The finite element analyses performed demonstrate that the cohesive fatigue accumulation function provides a link between the S-N diagram that describes crack initiation, and the Paris law that characterizes the rate of crack propagation. However, when the model was proposed, it was not known whether the form of the damage accumulation function associated with a desired S-N diagram is unique and, if not, if the link between S-N and the Paris law is unique and independent of the fatigue function selected. In the effort described herein, several alternative forms of the damage function that reproduce the desired features of S-N diagrams were found and evaluated. The effects of each of these functions on the predicted parameters of the Paris law and the propagation threshold are discussed. The results indicate that the predicted exponent m of the Paris law is indeed independent of the damage accumulation function. However, different functions predict different values for the pre-factor C of the Paris law. Therefore, the proper damage accumulation function must be selected by comparison with experiments. One of the new damage accumulation functions proposed herein was found to be particularly useful for analysis because of the ease with which the model parameters can be determined with a minimal amount of experimental information. The effectiveness of the proposed methodology and damage function was demonstrated by conducting analyses of a double cantilever beam test, a mixed-mode bending test, and a three-point bending test of a skin/doubler specimen. The results indicate that the same set of model parameters can provide accurate predictions of the rate of fatigue crack propagation for a variety of material interfaces, mode mixities, load levels, and stress ratios.

Cohesive elements↗

NASA NTRS2015

Investigation of Spiral Bevel Gear Condition Indicator Validation via AC-29-2C Combining Test Rig Damage Progression Data with Fielded Rotorcraft Data

This is the final of three reports published on the results of this project. In the first report, results were presented on nineteen tests performed in the NASA Glenn Spiral Bevel Gear Fatigue Test Rig on spiral bevel gear sets designed to simulate helicopter fielded failures. In the second report, fielded helicopter HUMS data from forty helicopters were processed with the same techniques that were applied to spiral bevel rig test data. Twenty of the forty helicopters experienced damage to the spiral bevel gears, while the other twenty helicopters had no known anomalies within the time frame of the datasets. In this report, results from the rig and helicopter data analysis will be compared for differences and similarities in condition indicator (CI) response. Observations and findings using sub-scale rig failure progression tests to validate helicopter gear condition indicators will be presented. In the helicopter, gear health monitoring data was measured when damage occurred and after the gear sets were replaced at two helicopter regimes. For the helicopters or tails, data was taken in the flat pitch ground 101 rotor speed (FPG101) regime. For nine tails, data was also taken at 120 knots true airspeed (120KTA) regime. In the test rig, gear sets were tested until damage initiated and progressed while gear health monitoring data and operational parameters were measured and tooth damage progression documented. For the rig tests, the gear speed was maintained at 3500RPM, a one hour run-in was performed at 4000 in-lb gear torque, than the torque was increased to 8000 in-lbs. The HUMS gear condition indicator data evaluated included Figure of Merit 4 (FM4), Root Mean Square (RMS) or Diagnostic Algorithm 1(DA1), + 3 Sideband Index (SI3) and + 1 Sideband Index (SI1). These were selected based on their sensitivity in detecting contact fatigue damage modes from analytical, experimental and historical helicopter data. For this report, the helicopter dataset was reduced to fourteen tails and the test rig data set was reduced to eight tested gear sets. The damage modes compared were separated into three cases. For case one, both the gear and pinion showed signs of contact fatigue or scuffing damage. For case two, only the pinion showed signs of contact fatigue damage or scuffing. Case three was limited to the gear tests when scuffing occurred immediately after the gear run-in. Results of this investigation highlighted the importance of understanding the complete monitored systems, for both the helicopter and test rig, before interpreting health monitoring data. Further work is required to better define these two systems that include better state awareness of the fielded systems, new sensing technologies, new experimental methods or models that quantify the effect of system design on CI response and new methods for setting thresholds that take into consideration the variance of each system.

spiral bevel gears↗