Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element 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.

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

At least 739 records · Page 41

Elastic-plastic crack analysis using a global-local approach on a parallel computer

A global-local finite-element analysis procedure is employed to solve elastoplastic crack problems. This procedure involves a global analysis using a coarse mesh and the subsequent locally refined analysis of a number of subregions. Parallel computation is involved in the local analyses. Numerical examples are performed on a MIMD parallel computer to demonstrate the time-saving capability of the proposed procedure.

Sun, C. T.↗

Structural Analysis of Test Flight Vehicles for Application of Multifunctional Energy Storage System

Under the NASA Aeronautics Research Mission Directorate (ARMD) Convergent Aeronautical Solutions (CAS) project, NASA Glenn Research Center has been leading Multifunctional Structures for High Energy Lightweight Load-bearing Storage (M-SHELLS) research efforts. The technology of integrating load-carrying structures with electrical energy storage capacity has the potential to reduce the overall weight of future electric aircraft. The proposed project goals were to develop M-SHELLS in the form of honeycomb coupons and subcomponents, integrate them into the structure, and conduct low-risk flight-tests onboard a remotely piloted small aircraft. Experimental M-SHELLS energy-storing coupons were fabricated and tested in the laboratory for their electrical and mechanical properties. In this report, finite element model development and structural analyses of two small test aircraft candidates are presented. The finite element analysis of the initial two-spar wing is described for strain, deflection, and weight estimation. After a test aircraft Tempest was acquired, a load-deflection test of the wing was conducted. A finite element model of the Tempest was then developed based on the test aircraft dimensions and construction detail. The component weight analyses from the finite element model and test measurements were correlated. Structural analysis results with multifunctional energy storage panels in the fuselage of the test vehicle are presented. Although the flight test was cancelled because of programmatic reasons and time constraints, the structural analysis results indicate that the mid-fuselage floor composite panel could provide structural integrity with minimal weight penalty while supplying electrical energy. To explore potential future applications of the multifunctional structure, analyses of the NASA X-57 Maxwell electric aircraft and a NASA N+3 Technology Conventional Configuration (N3CC) fuselage are presented. Secondary aluminum structures in the fuselage sub-floor and cargo area were partially replaced with reinforced five-layer composite panels with M-SHELLS honeycomb core. The N3CC fuselage weight reduction associated with each design without risking structural integrity are described. The structural analysis and weight estimation with the application of composite M-SHELLS panels to the N3CC fuselage indicate a 3.2% reduction in the fuselage structural weight, prior to accounting for the additional weight of core material required to complete the energy storage functionality.

Mukhopadhyay, Vivek↗

NASTRAN analysis for the Airmass Sunburst model 'C' Ultralight Aircraft

The purpose of this project was to create a three dimensional NASTRAN model of the Airmass Sunburst Ultralight comparable to one made for finite element analysis. A two dimensional sample problem will be calculated by hand and by NASTRAN to make sure that NASTRAN finds similar results. A three dimensional model, similar to the one analyzed by the finite element program, will be run on NASTRAN. A comparison will be done between the NASTRAN results and the finite element program results. This study will deal mainly with the aerodynamic loads on the wing and surrounding support structure at an attack angle of 10 degrees.

Verbestel, John↗

A formulation of rotor-airframe coupling for design analysis of vibrations of helicopter airframes

A linear formulation of rotor airframe coupling intended for vibration analysis in airframe structural design is presented. The airframe is represented by a finite element analysis model; the rotor is represented by a general set of linear differential equations with periodic coefficients; and the connections between the rotor and airframe are specified through general linear equations of constraint. Coupling equations are applied to the rotor and airframe equations to produce one set of linear differential equations governing vibrations of the combined rotor airframe system. These equations are solved by the harmonic balance method for the system steady state vibrations. A feature of the solution process is the representation of the airframe in terms of forced responses calculated at the rotor harmonics of interest. A method based on matrix partitioning is worked out for quick recalculations of vibrations in design studies when only relatively few airframe members are varied. All relations are presented in forms suitable for direct computer implementation.

Kvaternik, R. G.↗

Improvements to a method for the geometrically nonlinear analysis of compressively loaded stiffened composite panels

The NLPAN computer code uses a finite-strip approach to the analysis of thin-walled prismatic composite structures such as stiffened panels. The code can model in-plane axial loading, transverse pressure loading, and constant through-the-thickness thermal loading, and can account for shape imperfections. The NLPAN code represents an attempt to extend the buckling analysis of the VIPASA computer code into the geometrically nonlinear regime. Buckling mode shapes generated using VIPASA are used in NLPAN as global functions for representing displacements in the nonlinear regime. While the NLPAN analysis is approximate in nature, it is computationally economical in comparison with finite-element analysis, and is thus suitable for use in preliminary design and design optimization. A comprehensive description of the theoretical approach of NLPAN is provided. A discussion of some operational considerations for the NLPAN code is included. NLPAN is applied to several test problems in order to demonstrate new program capabilities, and to assess the accuracy of the code in modeling various types of loading and response. User instructions for the NLPAN computer program are provided, including a detailed description of the input requirements and example input files for two stiffened-panel configurations.

Stoll, Frederick↗

Stress Analysis and Permeability Testing of Cryogenic Composite Feed Line

For the next generation Single-Stage-To-Orbit (SSTO) Reusable Launch Vehicle (RLV), the use of advanced composite materials is highly desirable and critical to the success of the mission. NASA Marshall Space Flight Center (MSFC) has been working with the aerospace industry for many years to develop and demonstrate the cryogenic composite propellant tanks and feed lines technologies. A 50.8-mm diameter composite feed line for the Clipper Graham (DCY.A) was developed and tested. The purpose of the program is to demonstrate the LH2 permeability, composite to composite and metal joints, as well as composite flange interface of the composite feed line. Stress analysis and permeability testing have been performed on this article. Recently, a larger composite feed line design is being investigated and developed at MSFC for potential use in future RLV. The diameter of the feed line is 203 mm and the overall length is approximately 2.2 meters. This one piece unlined feed line consists of three straight tubular sections joined by two 90 degree elbows. The material chosen is IM7/977-3 prepreg fabric. The lay-up pattern is [0/90, plus or minus 45]s and is built up to 18 plies to the flanges at both ends. A preliminary stress analysis has been conducted to identify potential critical stresses and to develop the finite element analysis (FEA) capability of composite feed lines. As expected, the critical stresses occurred at the rims of some flange holes and the onset of the tapered tubular sections. Further analysis is required to determine the loads, flange deflection, vibration, and combined maximum loads. Two permeability-testing apparatuses were also designed for both flat panel specimens and curved feed line sections after impact damage. A larger permeant gas exposed area is required to accurately determine the effect of impact damage on the permeability of the feed line materials. The flat panel tester was fabricated and assembled. Three test coupons were made of graphite/epoxy plain weave prepreg and unidirectional prepreg. The flow rate was calculated by measuring the change of manometer fluid height per unit time. The permeability of the coupons is then calculated according to Darcy's Law. The pressure increase due to the head rise in the manometer was not considered for the initial estimation. The result of the unidirectional composite coupon agreed with that obtained from the DC-XA feed line testing of a similar material. Further improvement on the design may be required to improve the accuracy and efficiency of permeability testing.

Chu, Tsuchin Philip↗

Recent finite element studies in plasticity and fracture mechanics

The paper reviews recent work on fundamentals of elastic-plastic finite-element analysis and its applications to the mechanics of crack opening and growth in ductile solids. The presentation begins with a precise formulation of incremental equilibrium equations and their finite-element forms in a manner valid for deformations of arbitrary magnitude. Special features of computational procedures are outlined for accuracy in view of the near-incompressibility of elastic-plastic response. Applications to crack mechanics include the analysis of large plastic deformations at a progressively opening crack tip, the determination of J integral values and of limitations to J characterizations of the intensity of the crack tip field, and the determination of crack tip fields in stable crack growth.

Rice, J. R.↗

Finite element methods of analysis for 3D inviscid compressible flows

The applicants have developed a finite element based approach for the solution of three-dimensional compressible flows. The procedure enables flow solutions to be obtained on tetrahedral discretizations of computational domains of complex form. A further development was the incorporation of a solution adaptive mesh strategy in which the adaptivity is achieved by complete remeshing of the solution domain. During the previous year, the applicants were working with the Advanced Aerodynamics Concepts Branch at NASA Ames Research Center with an implementation of the basic meshing and solution procedure. The objective of the work to be performed over this twelve month period was the transfer of the adaptive mesh technology and also the undertaking of basic research into alternative flow algorithms for the Euler equations on unstructured meshes.

Peraire, Jaime↗

A finite element random response analysis of a complex panel with fluid-conveying pipes

A finite element approach is developed for a complex panel with fluid-conveying pipes undergoing large deflections subjected to random loadings. The influence of fluid velocity on the random response is investigated. The root-mean-square (rms) deflections and frequencies for different sound spectrum level value are studied. All four edges of the panel are considered to have the same conditions and restrictions. Either simply supported of clamped boundary condition with respect to the transverse deflection of the panel is considered. The prediction of fatigue life is then based on obtained rms stress. This analytical investigation will help to broaden the basic understanding of the role of fluid flow within structures subjected to acoustic loading.

Chiang, C. K.↗

A simple element for multilayer beams in NASTRAN thermal stress analysis

In the application of NASTRAN, structural members are usually represented by bar elements with multipoint constraint cards to enforce the interface conditions. While this is a very powerful method in principle, it was found that in practice the process for specification of constraints became tedious and error prone, unless the geometry was simple and the number of grid points low. An alternative approach was found within the framework of the NASTRAN program. This approach made use of the idea that a thermal distortion in a multilayer beam may be similar to a homogeneous beam with a thermal gradient across the cross section. The exact mathematical derivation for the equivalent beam, and all the necessary formulae for the equivalent parameters in NASTRAN analysis are presented. Some numerical examples illustrate the simplicity and ease of this approach for finite element analysis.

Chen, W. T.↗

Nonlinear analysis for the response and failure of compression-loaded angle-ply laminates with a hole

The objective was to determine the effect of nonlinear material behavior on the response and failure of unnotched and notched angle-ply laminates under uniaxial compressive loading. The endochronic theory was chosen as the constitutive theory to model the AS4/3502 graphite-epoxy material system. Three-dimensional finite element analysis incorporating the endochronic theory was used to determine the stresses and strains in the laminates. An incremental/iterative initial strain algorithm was used in the finite element program. To increase computational efficiency, a 180 deg rotational symmetry relationship was utilized and the finite element program was vectorized to run on a supercomputer. Laminate response was compared to experimentation revealing excellent agreement for both the unnotched and notched angle-ply laminates. Predicted stresses in the region of the hole were examined and are presented, comparing linear elastic analysis to the inelastic endochronic theory analysis. A failure analysis of the unnotched and notched laminates was performed using the quadratic tensor polynomial. Predicted fracture loads compared well with experimentation for the unnotched laminates, but were very conservative in comparison with experiments for the notched laminates.

Mathison, Steven R.↗

Analysis and Testing of a Metallic Repair Applicable to Pressurized Composite Aircraft Structure

Development of repair technology is vital to the long-term application of new structural concepts on aircraft structure. The design, analysis, and testing of a repair concept applicable to a stiffened composite panel based on the Pultruded Rod Stitched Efficient Unitized Structure was recently completed. The damage scenario considered was a mid-bay to mid-bay saw-cut with a severed stiffener, flange, and skin. A bolted metallic repair was selected so that it could be easily applied in the operational environment. The present work describes results obtained from tension and pressure panel tests conducted to validate both the repair concept and finite element analysis techniques used in the design effort. Simulation and experimental strain and displacement results show good correlation, indicating that the finite element modeling techniques applied in the effort are an appropriate compromise between required fidelity and computational effort. Static tests under tension and pressure loadings proved that the proposed repair concept is capable of sustaining load levels that are higher than those resulting from the current working stress allowables. Furthermore, the pressure repair panel was subjected to 55,000 pressure load cycles to verify that the design can withstand a life cycle representative for a transport category aircraft. These findings enable upward revision of the stress allowables that had been kept at an overly-conservative level due to concerns associated with repairability of the panels. This conclusion enables more weight efficient structural designs utilizing the composite concept under investigation.

Przekop, Adam↗

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. The transient finite-element analysis is done using TACS, and the cell-discharge battery model uses OpenMDAO and dymos. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗

Delamination growth analysis in quasi-isotropic laminates under loads simulating low-velocity impact

A geometrically nonlinear finite-element analysis was developed to calculate the strain energy released by delamination plates during impact loading. Only the first mode of deformation, which is equivalent to static deflection, was treated. Both the impact loading and delamination in the plate were assumed to be axisymmetric. The strain energy release rate in peeling, G sub I, and shear sliding, G sub II, modes were calculated using the fracture mechanics crack closure technique. Energy release rates for various delamination sizes and locations and for various plate configurations and materials were compared. The analysis indicated that shear sliding (G sub II) was the primary mode of delamination growth. The analysis also indicated that the midplane (maximum transverse shear stress plane) delamination was more critical and would grow before any other delamination of the same size near the midplane region. The delamination growth rate was higher (neutrally stable) for a low toughness (brittle) matrix and slower (stable) for high toughness matrix. The energy release rate in the peeling mode, G sub I, for a near-surface delamination can be as high as 0.5G sub II and can contribute significantly to the delamination growth.

Shivakumar, K. N.↗

Delamination growth analysis in quasi-isotropic laminates under loads simulating low-velocity impact

A geometrically nonlinear finite-element analysis has been developed to calculate the strain energy released by delaminating plates during impact loading. Only the first mode of deformation, which is equivalent to static deflection, was treated. Both the impact loading and delamination in the plate were assumed to be axisymmetric. The strain energy release rate in peeling, GI, and shear sliding, GII, modes were calculated using the fracture mechanics crack closure technique. Energy release rates for various delamination sizes and locations and for various plate configurations and materials were compared. The analysis indicated that shear sliding was the primary mode of delamination growth. The analysis also indicated that the midplane (maximum transverse shear stress plane) delamination was more critical and would grow first before any other delamination of the same size near the midplane region. The delamination growth rate was higher (neutrally stable) for a low toughness (brittle) matrix and slower (stable) for high toughness matrix. The energy release rate in the peeling mode, GI, for a near-surface delamination can be as high as 0.5GII, and can contribute significantly to the delamination growth.

Shivakumar, K. N.↗

Analysis of notched metal matrix composites under tension loading

Presented are techniques based on 3-D finite-element analysis for the examination of continuous fiber reinforced metal matrix composites. Examples are shown for specific metal matrix composites such as boron/aluminum and silicon-carbide/aluminum. Specimen stress-strain behavior and stress at first fiber failure were predicted for boron/aluminum laminates containing circular holes and crack-like slits. The predictions compared very well for (+ or - 45) sub 2s laminates. Mesh configuration was shown to have an effect on the calculation of stresses local to the notch. The presence of thin interface layers of matrix material had a significant influence on the slit-tip stress state, causing sharper stress gradients near the notch. Interface layers reduced the slit-tip fibers stresses in a (+ or - 45) sub s silicon-carbide/aluminum laminate but increased them in a (0/90) sub s laminate.

Bigelow, C. A.↗

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.↗

Stress analysis method for clearance-fit joints with bearing-bypass loads

Within a multi-fastener joint, fastener holes may be subjected to the combined effects of bearing loads and loads that bypass the hole to be reacted elsewhere in the joint. The analysis of a joint subjected to search combined bearing and bypass loads is complicated by the usual clearance between the hole and the fastener. A simple analysis method for such clearance-fit joints subjected to bearing-bypass loading has been developed in the present study. It uses an inverse formulation with a linear elastic finite-element analysis. Conditions along the bolt-hole contact arc are specified by displacement constraint equations. The present method is simple to apply and can be implemented with most general purpose finite-element programs since it does not use complicated iterative-incremental procedures. The method was used to study the effects of bearing-bypass loading on bolt-hole contact angles and local stresses. In this study, a rigid, frictionless bolt was used with a plate having the properties of a quasi-isotropic graphite/epoxy laminate. Results showed that the contact angle as well as the peak stresses around the hole and their locations were strongly influenced by the ratio of bearing and bypass loads. For single contact, tension and compression bearing-bypass loading had opposite effects on the contact angle. For some compressive bearing-bypass loads, the hole tended to close on the fastener leading to dual contact. It was shown that dual contact reduces the stress concentration at the fastener and would, therefore, increase joint strength in compression. The results illustrate the general importance of accounting for bolt-hole clearance and contact to accurately compute local bolt-hole stresses for combined bearings and bypass loading.

Naik, R. A.↗