Results for “continuous finite element” — Page 8

NASA NTRS2008

Simplified Shear Solution for Determination of the Shear Stress Distribution in a Composite Panel from the Applied Shear Resultant

The simplified shear solution method is presented for approximating the through-thickness shear stress distribution within a composite laminate or panel based on laminated beam theory. The method does not consider the solution of a particular boundary value problem; rather it requires only knowledge of the global shear loading, geometry, and material properties of the laminate or panel. It is thus analogous to lamination theory in that ply level stresses can be efficiently determined from global load resultants (as determined, for instance, by finite element analysis) at a given location in a structure and used to evaluate the margin of safety on a ply by ply basis. The simplified shear solution stress distribution is zero at free surfaces, continuous at ply boundaries, and integrates to the applied shear load. Comparisons to existing theories are made for a variety of laminates, and design examples are provided illustrating the use of the method for determining through-thickness shear stress margins in several types of composite panels and in the context of a finite element structural analysis.

Bednarcyk, Brett A.↗

NASA NTRS2015

Simulation of Delamination-Migration and Core Crushing in a CFRP Sandwich Structure

Following the onset of damage caused by an impact load on a composite laminate structure, delaminations often form propagating outwards from the point of impact and in some cases can migrate via matrix cracks between plies as they grow. The goal of the present study is to develop an accurate finite element modeling technique for simulation of the delamination-migration phenomena in laminate impact damage processes. An experiment was devised where, under a quasi-static indentation load, an embedded delamination in the facesheet of a laminate sandwich specimen migrates via a transverse matrix crack and then continues to grow on a new ply interface. The quasistatic nature of the indentation results in structural behavior equivalent to that seen in low-velocity impact and also allows for highly detailed real time damage characterization. Several finite element damage simulation methods were investigated. Comparing the experimental results with those of the different models reveals certain modeling features that are important to include in a numerical simulation of delamination-migration and some that may be neglected.

McElroy, M.↗

NASA NTRS2003

Simulating the Effect of Space Vehicle Environments on Directional Solidification of a Binary Alloy

Space microgravity missions are designed to provide a microgravity environment for scientific experiments, but these missions cannot provide a perfect environment, due to vibrations caused by crew activity, on-board experiments, support systems (pumps, fans, etc.), periodic orbital maneuvers, and water dumps. Therefore, it is necessary to predict the impact of these vibrations on space experiments, prior to performing them. Simulations were conducted to study the effect of the vibrations on the directional solidification of a dendritic alloy. Finite element ca!cu!attie?ls were dme with a simd2titcr based on a continuum model of dendritic solidification, using the Fractional Step Method (FSM). The FSM splits the solution of the momentum equation into two steps: the viscous intermediate step, which does not enforce continuity; and the inviscid projection step, which calculates the pressure and enforces continuity. The FSM provides significant computational benefits for predicting flows in a directionally solidified alloy, compared to other methods presently employed, because of the efficiency gains in the uncoupled solution of velocity and pressure. finite differences, arises when the interdendritic liquid reaches the eutectic temperature and concentration. When a node reaches eutectic temperature, it is assumed that the solidification of the eutectic liquid continues at constant temperature until all the eutectic is solidified. With this approach, solidification is not achieved continuously across an element; rather, the element is not considered solidified until the eutectic isotherm overtakes the top nodes. For microgravity simulations, where the convection is driven by shrinkage, it introduces large variations in the fluid velocity. When the eutectic isotherm reaches a node, all the eutectic must be solidified in a short period, causing an abrupt increase in velocity. To overcome this difficulty, we employed a scheme to numerically predict a more accurate value for the rate of change of fraction of liquid as the liquid in an element solidifies. The new method enables us to contrast results of simulations in which the alloy is subjected to no gravity or a steady-state acceleration versus simulations when the alloy is subjected to vibration disturbances; therefore, the effect of vibration disturbances can be assessed more accurately. To assess the impact of these vibration-perturbations, transient accelerometer data from a space shuttle mission are used as inputs for the simulation model. These on-orbit acceleration data were obtained from the Microgravity Science Division at Glenn Research Center (GRC- MSD) and are applied to the buoyancy term of the momentum equation in a simulation of a Pb-5.8 wt. % Sb alloy that solidifies in a thermal gradient of 4000 K/m and a translation velocity of 3 p d s . Figure 2 shows the vertical velocity of a node that begins in the all-liquid region and subsequently solidifies; the vibrations are applied at 5000 seconds in this simulation. An important difficulty, common to all solidification models based on finite elements or 2 The magnitudes of the velocity oscillations that are vibration-induced are very small and acceptable. The biggest concern is whether the concentration of the liquid near the dendrite tips is distorted because of the vibration-induced perturbations. Results for this case show no concentration oscillations present in the all-liquid region.

Westra, D. G.↗

NASA NTRS1979

Recent developments in finite element analysis for transonic airfoils

The prediction of aerodynamic forces in the transonic regime generally requires a flow field calculation to solve the governing non-linear mixed elliptic-hyperbolic partial differential equations. Finite difference techniques were developed to the point that design and analysis application are routine, and continual improvements are being made by various research groups. The principal limitation in extending finite difference methods to complex three-dimensional geometries is the construction of a suitable mesh system. Finite element techniques are attractive since their application to other problems have permitted irregular mesh elements to be employed. The purpose of this paper is to review the recent developments in the application of finite element methods to transonic flow problems and to report some recent results.

Hafez, M. M.↗

NASA NTRS1992

Global/local analysis of laminated composite plates using variable kinematic finite elements

A finite element modeling methodology is developed for the hierarchical, global/local analysis of laminated composite plates. The method incorporates a new variable kinematics, displacement-based, finite element that is developed using a multiple assumed displacement field approach. The variable kinematic elements provide a great degree of flexibility in defining the transverse (through thickness) variation of the assumed displacement field. The resulting finite element model permits different subregions of the computational domain to be described by different mathematical models. Enforcing displacement continuity along subregion boundaries requires only the specification of certain homogeneous essential boundary conditions, thus avoiding the inconvenience of multi-point constraints, penalty function methods, or special transition elements.

Robbins, D. H., Jr.↗

NASA NTRS1992

Stability of generally stiffened anisotropic noncircular cylinders

Continuous filament grid-stiffened structure is a stiffening concept that combines structural efficiency and damage tolerance. However, finite element design of such structures against buckling is expensive due to the complexities of the structure. An analytical model of such a structure is developed using a penalty method (artificial springs) with a first order shear deformation theory (FSDT). The buckling analysis under combined loadings is done using energy method with a penalty/Rayleigh-Ritz technique. The penalty/Rayleigh-Ritz approach is computationally less demanding when compared to the finite element solution and mesh generation. Apart from the published research works on buckling of stiffened plates and shells by finite element and finite strips, research works on buckling of stiffened plates and shells utilize three different approaches; smeared, column, and discrete approaches. The discrete approach considers the discrete effects of the stiffeners in the buckling behavior by modeling stiffeners as line of bending (EI) and torsion (GJ) stiffnesses on panel skin. Some local deformations are lost when stiffeners are modeled as (EI) and (GJ) stiffeners. This approach becomes difficult in the case of plate stiffened in more than two directions. Most of the work done using the discrete approach involved the Classical Plate Theory (CLPT) rather than the FSDT. We report on our formulation of a discrete approach coupled with a penalty formulation and FSDT.

Sobh, Nahil Atef↗

NASA NTRS1993

A multigrid solver for the semiconductor equations

We present a multigrid solver for the exponential fitting method. The solver is applied to the current continuity equations of semiconductor device simulation in two dimensions. The exponential fitting method is based on a mixed finite element discretization using the lowest-order Raviart-Thomas triangular element. This discretization method yields a good approximation of front layers and guarantees current conservation. The corresponding stiffness matrix is an M-matrix. 'Standard' multigrid solvers, however, cannot be applied to the resulting system, as this is dominated by an unsymmetric part, which is due to the presence of strong convection in part of the domain. To overcome this difficulty, we explore the connection between Raviart-Thomas mixed methods and the nonconforming Crouzeix-Raviart finite element discretization. In this way we can construct nonstandard prolongation and restriction operators using easily computable weighted L(exp 2)-projections based on suitable quadrature rules and the upwind effects of the discretization. The resulting multigrid algorithm shows very good results, even for real-world problems and for locally refined grids.

Bachmann, Bernhard↗

NASA NTRS1995

Damage-mitigating control of a reusable rocket engine for high performance and extended life

The goal of damage mitigating control in reusable rocket engines is to achieve high performance with increased durability of mechanical structures such that functional lives of the critical components are increased. The major benefit is an increase in structural durability with no significant loss of performance. This report investigates the feasibility of damage mitigating control of reusable rocket engines. Phenomenological models of creep and thermo-mechanical fatigue damage have been formulated in the state-variable setting such that these models can be combined with the plant model of a reusable rocket engine, such as the Space Shuttle Main Engine (SSME), for synthesizing an optimal control policy. Specifically, a creep damage model of the main thrust chamber wall is analytically derived based on the theories of sandwich beam and viscoplasticity. This model characterizes progressive bulging-out and incremental thinning of the coolant channel ligament leading to its eventual failure by tensile rupture. The objective is to generate a closed form solution of the wall thin-out phenomenon in real time where the ligament geometry is continuously updated to account for the resulting deformation. The results are in agreement with those obtained from the finite element analyses and experimental observation for both Oxygen Free High Conductivity (OFHC) copper and a copper-zerconium-silver alloy called NARloy-Z. Due to its computational efficiency, this damage model is suitable for on-line applications of life prediction and damage mitigating control, and also permits parametric studies for off-line synthesis of damage mitigating control systems. The results are presented to demonstrate the potential of life extension of reusable rocket engines via damage mitigating control. The control system has also been simulated on a testbed to observe how the damage at different critical points can be traded off without any significant loss of engine performance. The research work reported here is built upon concepts derived from the disciplines of Controls, Thermo-fluids, Structures, and Materials. The concept of damage mitigation, as presented in this report, is not restricted to control of rocket engines. It can be applied to any system where structural durability is an important issue.

Ray, Asok↗

NASA NTRS2009

A Theoretical Investigation of Composite Overwrapped Pressure Vessel (COPV) Mechanics Applied to NASA Full Scale Tests

A theoretical investigation of the factors controlling the stress rupture life of the National Aeronautics and Space Administration's (NASA) composite overwrapped pressure vessels (COPVs) continues. Kevlar (DuPont) fiber overwrapped tanks are of particular concern due to their long usage and the poorly understood stress rupture process in Kevlar filaments. Existing long term data show that the rupture process is a function of stress, temperature and time. However due to the presence of a load sharing liner, the manufacturing induced residual stresses and the complex mechanical response, the state of actual fiber stress in flight hardware and test articles is not clearly known. This paper is a companion to a previously reported experimental investigation and develops a theoretical framework necessary to design full-scale pathfinder experiments and accurately interpret the experimentally observed deformation and failure mechanisms leading up to static burst in COPVs. The fundamental mechanical response of COPVs is described using linear elasticity and thin shell theory and discussed in comparison to existing experimental observations. These comparisons reveal discrepancies between physical data and the current analytical results and suggest that the vessel s residual stress state and the spatial stress distribution as a function of pressure may be completely different from predictions based upon existing linear elastic analyses. The 3D elasticity of transversely isotropic spherical shells demonstrates that an overly compliant transverse stiffness relative to membrane stiffness can account for some of this by shifting a thin shell problem well into the realm of thick shell response. The use of calibration procedures are demonstrated as calibrated thin shell model results and finite element results are shown to be in good agreement with the experimental results. The successes reported here have lead to continuing work with full scale testing of larger NASA COPV hardware.

Thesken, John C.↗

NASA NTRS2006

A Theoretical Investigation of Composite Overwrapped Pressure Vessel (COPV) Mechanics Applied to NASA Full Scale Tests

A theoretical investigation of the factors controlling the stress rupture life of the National Aeronautics and Space Agency's (NASA) composite overwrapped pressure vessels (COPVs) continues. Kevlar(TradeMark) fiber overwrapped tanks are of particular concern due to their long usage and the poorly understood stress rupture process in Kevlar(TradeMark) filaments. Existing long term data show that the rupture process is a function of stress, temperature and time. However, due to the presence of a load sharing liner, the manufacturing induced residual stresses and the complex mechanical response, the state of actual fiber stress in flight hardware and test articles is not clearly known. This paper is a companion to the experimental investigation reported in [1] and develops a theoretical framework necessary to design full-scale pathfinder experiments and accurately interpret the experimentally observed deformation and failure mechanisms leading up to static burst in COPVs. The fundamental mechanical response of COPVs is described using linear elasticity and thin shell theory and discussed in comparison to existing experimental observations. These comparisons reveal discrepancies between physical data and the current analytical results and suggest that the vessel's residual stress state and the spatial stress distribution as a function of pressure may be completely different from predictions based upon existing linear elastic analyses. The 3D elasticity of transversely isotropic spherical shells demonstrates that an overly compliant transverse stiffness relative to membrane stiffness can account for some of this by shifting a thin shell problem well into the realm of thick shell response. The use of calibration procedures are demonstrated as calibrated thin shell model results and finite element results are shown to be in good agreement with the experimental results. The successes reported here have lead to continuing work with full scale testing of larger NASA COPV hardware.

Greene, N.↗

NASA NTRS1996

A Systematic Solution Approach for Neutron Transport Problems in Diffuse Regimes

A systematic solution approach for the neutron transport equation, based on a least-squares finite-element discretization, is presented. This approach includes the theory for the existence and uniqueness of the analytical as well as of the discrete solution, bounds for the discretization error, and guidance for the development of an efficient multigrid solver for the resulting discrete problem. To guarantee the accuracy of the discrete solution for diffusive regimes, a scaling transformation is applied to the transport operator prior to the discretization. The key result is the proof of the V-ellipticity and continuity of the scaled least-squares bilinear form with constants that are independent of the total cross section and the absorption cross section. For a variety of least-squares finite-element discretizations this leads to error bounds that remain valid in diffusive regimes. Moreover, for problems in slab geometry a full multigrid solver is presented with V(1, 1)-cycle convergence rates approximately equal to 0.1, independent of the size of the total cross section and the absorption cross section.

Manteuffel, T. A.↗

NASA NTRS2000

Advances in Fatigue and Fracture Mechanics Analyses for Metallic Aircraft Structures

This paper reviews some of the advances that have been made in stress analyses of cracked aircraft components, in the understanding of the fatigue and fatigue-crack growth process, and in the prediction of residual strength of complex aircraft structures with widespread fatigue damage. Finite-element analyses of cracked metallic structures are now used to determine accurate stress-intensity factors for cracks at structural details. Observations of small-crack behavior at open and rivet-loaded holes and the development of small-crack theory has lead to the prediction of stress-life behavior for components with stress concentrations under aircraft spectrum loading. Fatigue-crack growth under simulated aircraft spectra can now be predicted with the crack-closure concept. Residual strength of cracked panels with severe out-of-plane deformations (buckling) in the presence of stiffeners and multiple-site damage can be predicted with advanced elastic-plastic finite-element analyses and the critical crack-tip-opening angle (CTOA) fracture criterion. These advances are helping to assure continued safety of aircraft structures.

Newman, J. C., Jr.↗

NASA NTRS1999

Advances in Fatigue and Fracture Mechanics Analyses for Aircraft Structures

This paper reviews some of the advances that have been made in stress analyses of cracked aircraft components, in the understanding of the fatigue and fatigue-crack growth process, and in the prediction of residual strength of complex aircraft structures with widespread fatigue damage. Finite-element analyses of cracked structures are now used to determine accurate stress-intensity factors for cracks at structural details. Observations of small-crack behavior at open and rivet-loaded holes and the development of small-crack theory has lead to the prediction of stress-life behavior for components with stress concentrations under aircraft spectrum loading. Fatigue-crack growth under simulated aircraft spectra can now be predicted with the crack-closure concept. Residual strength of cracked panels with severe out-of-plane deformations (buckling) in the presence of stiffeners and multiple-site damage can be predicted with advanced elastic-plastic finite-element analyses and the critical crack-tip-opening angle (CTOA) fracture criterion. These advances are helping to assure continued safety of aircraft structures.

Newman, J. C., Jr.↗

NASA NTRSDate not supplied

Local Analysis-Test Correlation Of Tow-Steered Composite Shells With Small Cutouts

The prebuckling and postbuckling behavior of two composite tow-steered shells with small cutouts is assessed using nonlinear dynamic finite element analyses and compared in detail with experimental measurements. The cylindrical shells were manufactured without cutouts using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from ±10 degrees on the axially stiff crown and keel, to ±45 degrees on the shear-stiff sides. The first shell with overlaps has laminate thickness variations on the crown and keel that result from application of all 24 tows during each pass of the fiber placement system. The second shell without overlaps uses the fiber placement system’s tow drop/add capability to achieve a more uniform shell wall thickness. An unreinforced cutout representing a passenger door on a commercial aircraft fuselage is machined into the side of each of the two shells. These shells with cutouts were tested in axial compression and buckled elastically in previous work. Detailed nonlinear finite element analysis results are compared with their corresponding measured local load-displacement and load-strain responses in prebuckling, at global buckling, and into a stable postbuckled state. Test data from displacement transducers, strain gages, and digital image correlation are extracted at the centers of the crown and keel, and at the middles of the top and left edges of the cutout. The overall agreement between these measured and analytical responses is excellent in prebuckling through global buckling, and very good from global buckling through postbuckling. As such, the excellent correlation observed here increases confidence in applying tow-steered composites in operational vehicles.

Composites↗

NASA NTRSDate not supplied

Local Analysis-Test Correlation Of Tow-Steered Composite Shells With Small Cutouts

The prebuckling and postbuckling behavior of two composite tow-steered shells with small cutouts is assessed using nonlinear dynamic finite element analyses and compared in detail with experimental measurements. The cylindrical shells were manufactured without cutouts using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from ±10 degrees on the axially stiff crown and keel, to ±45 degrees on the shear-stiff sides. The first shell with overlaps has laminate thickness variations on the crown and keel that result from application of all 24 tows during each pass of the fiber placement system. The second shell without overlaps uses the fiber placement system’s tow drop/add capability to achieve a more uniform shell wall thickness. An unreinforced cutout representing a passenger door on a commercial aircraft fuselage is machined into the side of each of the two shells. These shells with cutouts were tested in axial compression and buckled elastically in previous work. Detailed nonlinear finite element analysis results are compared with their corresponding measured local load-displacement and load-strain responses in prebuckling, at global buckling, and into a stable postbuckled state. Test data from displacement transducers, strain gages, and digital image correlation are extracted at the centers of the crown and keel, and at the middles of the top and left edges of the cutout. The overall agreement between these measured and analytical responses is excellent in prebuckling through global buckling, and very good from global buckling through postbuckling. As such, the excellent correlation observed here increases confidence in applying tow-steered composites in operational vehicles.

Composites↗

NASA NTRSDate not supplied

Simplified Applied Loads of a Transonic Truss-Braced Wing Derived from Computational Fluid Dynamics

The transonic truss-braced wing (TTBW) concept is under consideration by the NASA Advanced Air Transport Technology (AATT) project due to its potential to reduce fuel consumption of transport aircraft. To properly evaluate the TTBW concept, detailed structural modeling and subsequent test planning is underway, which requires a simplified set of discrete applied loads that simulate the distributed aerodynamic pressures. Given the unique architecture of the TTBW and its intended operation in the transonic regime, computational fluid dynamics (CFD) is the ideal approach to employ to ensure the aerodynamics of this configuration are sufficiently captured. To convert aerodynamic results into input for structural models, application of a data transfer tool, Discrete Data Transfer Between Dissimilar Meshes (DDTBDM), is presented to resolve aerodynamic pressures from CFD into single load vectors at hardpoints along both the wing and the strut. Evaluation of mapping the aerodynamic forces to hardpoints at rib and spar intersections in a structural model is performed, and the results are also compared with a linear static aeroelastic analysis. The structural responses indicate that distributing forces along only the wing and strut forward spars is insufficient for simulating the distributed aerodynamic loads. Adding moments along the forward spars improved the global response of the wing, but the point moments created unrealistic local rotations within the structure. When distributing loads along both the forward and aft spars of the wing and strut, the global response was less sensitive to the inclusion of point moments, which indicated the arrangement of point forces alone was sufficient input for structural design work and testing and provided a more continuous loading profile.

Finite Element↗

NASA NTRSDate not supplied

Simplified Applied Loads of a Transonic Truss-Braced Wing Derived from Computational Fluid Dynamic

The transonic truss-braced wing (TTBW) concept is under consideration by the NASA Advanced Air Transport Technology (AATT) project due to its potential to reduce fuel consumption of transport aircraft. To properly evaluate the TTBW concept, detailed structural modeling and subsequent test planning is underway, which requires a simplified set of discrete applied loads that simulate the distributed aerodynamic pressures. Given the unique architecture of the TTBW and its intended operation in the transonic regime, computational fluid dynamics (CFD) is the ideal approach to employ to ensure the aerodynamics of this configuration are sufficiently captured. To convert aerodynamic results into input for structural models, application of a data transfer tool, Discrete Data Transfer Between Dissimilar Meshes (DDTBDM), is presented to resolve aerodynamic pressures from CFD into single load vectors at hardpoints along both the wing and the strut. Evaluation of mapping the aerodynamic forces to hardpoints at rib and spar intersections in a structural model is performed, and the results are also compared with a linear static aeroelastic analysis. The structural responses indicate that distributing forces along only the wing and strut forward spars is insufficient for simulating the distributed aerodynamic loads. Adding moments along the forward spars improved the global response of the wing, but the point moments created unrealistic local rotations within the structure. When distributing loads along both the forward and aft spars of the wing and strut, the global response was less sensitive to the inclusion of point moments, which indicated the arrangement of point forces alone was sufficient input for structural design work and testing and provided a more continuous loading profile.

Finite Element↗

NASA NTRS1981

Finite element analysis of inviscid subsonic boattail flow

A finite element code for analysis of inviscid subsonic flows over arbitrary nonlifting planar or axisymmetric bodies is described. The code solves a novel primitive variable formulation of the coupled irrotationality and compressible continuity equations. Results for flow over a cylinder, a sphere, and a NACA 0012 airfoil verify the code. Computed subcritical flows over an axisymmetric boattailed afterbody compare well with finite difference results and experimental data. Interative coupling with an integral turbulent boundary layer code shows strong viscous effects on the inviscid flow. Improvements in code efficiency and extensions to transonic flows are discussed.

Chima, R. V.↗