Search NASA⌕ Search

SEARCH · Search NASA

Results for “layered solid element”

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 19 records

Progressive Damage Analysis of Bonded Composite Joints

The present work is related to the development and application of progressive damage modeling techniques to bonded joint technology. The joint designs studied in this work include a conventional composite splice joint and a NASA-patented durable redundant joint. Both designs involve honeycomb sandwich structures with carbon/epoxy facesheets joined using adhesively bonded doublers.Progressive damage modeling allows for the prediction of the initiation and evolution of damage within a structure. For structures that include multiple material systems, such as the joint designs under consideration, the number of potential failure mechanisms that must be accounted for drastically increases the complexity of the analyses. Potential failure mechanisms include fiber fracture, intraply matrix cracking, delamination, core crushing, adhesive failure, and their interactions. The bonded joints were modeled using highly parametric, explicitly solved finite element models, with damage modeling implemented via custom user-written subroutines. Each ply was discretely meshed using three-dimensional solid elements. Layers of cohesive elements were included between each ply to account for the possibility of delaminations and were used to model the adhesive layers forming the joint. Good correlation with experimental results was achieved both in terms of load-displacement history and the predicted failure mechanism(s).

Leone, Frank A., Jr.↗

A spectral element approach to wave motion in layered solids

A matrix methodology similar to that of the finite element method is developed for the analysis of stress waves in layered solids. Because the mass distribution is modeled exactly, the approach gives the exact frequency response of each layer. The fast Fourier transform and Fourier series are used for inversion to the time/space domain. The impact of a structured medium with multiple layers is used to demonstrate the method. Comparison with existing propagator and direct global matrix methods show the present approach to be computationally more efficient.

Rizzi, S. A.↗

A Navier-Stokes solver for high speed equilibrium flows and application to blunt bodies

This paper presents a finite element method for the solution of Navier-Stokes equations with the assumption of thermodynamic and chemical equilibrium. The method employs an upwind finite element technique with an implicit time-marching scheme for the solution, and uses an adaptively generated unstructured triangular mesh with several layers of quadrilateral elements near solid walls. The complexity associated with the assumption that the flow is in equilibrium is treated consistently, and the inviscid flux Jacobian matrices are derived. Several problems involving inviscid and viscous hypersonic flow past blunt are solved. Results are compared with other numerical results and experimental data, and generally good agreement is observed.

Prabhu, Ramadas K.↗

Postseismic rebound due to creep of the lower lithosphere and asthenosphere

Postseismic surface deformations are attributed to the inelastic flow of the subcrustal regions of the Earth following an earthquake. A multilayer representation of the Earth's rheological properties is used in conjunction with a finite element computational scheme to calculate time dependent displacements and strains subsequent to a strike slip earthquake. The deviatoric stress strain relations for the uppermost layer is assumed elastic. Lower layers are assumed to be, in order of increasing depth, a standard, linear, three element, viscoelastic solid; a linear viscoelastic fluid; and another elastic solid. Physically these layers correspond to the upper lithosphere, lower lithosphere, asthenosphere, and lower mantle, respectively. Elastic dilatational properties are assumed throughout. Appreciable postseismic displacements, possibly approaching meters, for large earthquakes, arise from viscoelastic relaxation following sudden coseismic slip. Furthermore, compared to the simpler case of an elastic lithosphere over a viscoelastic asthenosphere and the near fault postseismic shear, strain is increased, by a factor of two or more in some cases, by the presence of a viscoelastic lower lithosphere. Also, the duration of postseismic straining is increased by the relatively slow relaxation of this layer.

Cohen, S. C.↗

Mixed Element Type Unstructured Grid Generation for Viscous Flow Applications

A procedure is presented for efficient generation of high-quality unstructured grids suitable for CFD simulation of high Reynolds number viscous flow fields. Layers of anisotropic elements are generated by advancing along prescribed normals from solid boundaries. The points are generated such that either pentahedral or tetrahedral elements with an implied connectivity can be be directly recovered. As points are generated they are temporarily attached to a volume triangulation of the boundary points. This triangulation allows efficient local search algorithms to be used when checking merging layers, The existing advancing-front/local-reconnection procedure is used to generate isotropic elements outside of the anisotropic region. Results are presented for a variety of applications. The results demonstrate that high-quality anisotropic unstructured grids can be efficiently and consistently generated for complex configurations.

Marcum, David L.↗

Finite Element Modeling of the Buckling Response of Sandwich Panels

A comparative study of different modeling approaches for predicting sandwich panel buckling response is described. The study considers sandwich panels with anisotropic face sheets and a very thick core. Results from conventional analytical solutions for sandwich panel overall buckling and face-sheet-wrinkling type modes are compared with solutions obtained using different finite element modeling approaches. Finite element solutions are obtained using layered shell element models, with and without transverse shear flexibility, layered shell/solid element models, with shell elements for the face sheets and solid elements for the core, and sandwich models using a recently developed specialty sandwich element. Convergence characteristics of the shell/solid and sandwich element modeling approaches with respect to in-plane and through-the-thickness discretization, are demonstrated. Results of the study indicate that the specialty sandwich element provides an accurate and effective modeling approach for predicting both overall and localized sandwich panel buckling response. Furthermore, results indicate that anisotropy of the face sheets, along with the ratio of principle elastic moduli, affect the buckling response and these effects may not be represented accurately by analytical solutions. Modeling recommendations are also provided.

Rose, Cheryl A.↗

Progressive Failure Analysis Correlation with Notched Composite Laminate Test Data

Testing of small and intermediate laminate panels with notches was performed and Progressive Failure Analyses (PFA) models were developed for aid in test planning and for correlation with the test data. Two progressive damage failure models were included in the study: the commercially available Abaqus built-in damage model and COmplete STress Reduction (COSTR) damage model developed at NASA Langley Research Center (LaRC). The finite elements models used for the analysis were developed using shell elements. The pre-test PFA results obtained from the two damage models were compared to test data. Then a post-test PFA with updated material properties based on additional available material property data was executed with the COSTR damage model and compared to the test data. The panels tested in compression all exhibited less scatter in the failure load with self-similar failure behavior. The panels tested in tension exhibited more scatter in the failure load with a failure mode that included delamination of plies. However, a self-similar crack path was also observed as the overall failure mode for all panels tested. When comparing test data to the pre-test analysis results, both damage models with shell element models were considered adequate for predicting the behavior and failure load of both the small and intermediate panels when tested in compression. However, for the tension load cases, larger differences were observed between the test and analysis results. When comparing test data to the post-test analysis performed using the COSTR damage model, better correlation was observed. However, large discrepancies were still observed for the tension load cases. Consequently, a higher fidelity finite element model including solid elements for sub-laminates and the modeling of cohesive layers between sub-laminate layers is recommended for tension loading of notched composite laminates.

Damage Model,Progressive Failure Analysis,Open Hol↗

Postseismic deformation due to subcrustal viscoelastic relaxation following dip-slip earthquakes

The deformation of the Earth following a dip-slip earthquake is calculated using a three layer rheological model and finite element techniques. The three layers are an elastic upper lithosphere, a standard linear solid lower lithosphere, and a Maxwell viscoelastic asthenosphere-a model previously analyzed in the strike-clip case (Cohen, 1981, 1982). Attention is focused on the magnitude of the postseismic subsidence and the width of the subsidence zone that can develop due to the viscoelastic response to coseismic reverse slip. Detailed analysis for a fault extending from the surface to 15 km with a 45 deg dip reveals that postseismic subsidence is sensitive to the depth to the asthenosphere but is only weakly dependent on lower lithosphere depth. The greatest subsidence occurs when the elastic lithosphere is about 30 km thick and the asthenosphere lies just below this layer (asthenosphere depth = 2 times the fault depth). The extremum in the subsidence pattern occurs at about 5 km from the surface trace of the fault and lies over the slip plane. In a typical case after a time t = 30 tau (tau = Maxwell time) following the earthquake the subsidence at this point is 60% of the coseismic uplift. Unlike the horizontal deformation following a strike slip earthquake, significant vertical deformation due to asthenosphere flow persists for many times tau and the magnitude of the vertical deformation is not necessarily enhanced by having a partially relaxing lower lithosphere.

Cohen, S. C.↗

Postseismic deformation due to subcrustal viscoelastic relaxation following dip-slip earthquakes

The deformation of the earth following a dip-slip earthquake is calculated using a three layer rheological model and finite element techniques. The three layers are an elastic upper lithosphere, a standard linear solid lower lithosphere, and a Maxwell viscoelastic asthenosphere - a model previously analyzed in the strike-clip case (Cohen, 1981, 1982). Attention is focused on the magnitude of the postseismic subsidence and the width of the subsidence zone that can develop due to the viscoelastic response to coseismic reverse slip. Detailed analysis for a fault extending from the surface to 15 km with a 45 deg dip reveals that postseismic subsidence is sensitive to the depth to the asthenosphere but is only weakly dependent on lower lithosphere depth. The greatest subsidence occurs when the elastic lithosphere is about 30 km thick and the asthenosphere lies just below this layer (asthenosphere depth = 2 times the fault depth). The extremum in the subsidence pattern occurs at about 5 km from the surface trace of the fault and lies over the slip plane. In a typical case after a time t = 30 tau (tau = Maxwell time) following the earthquake, the subsidence at this point is 60 percent of the coseismic uplift. Unlike the horizontal deformation following a strike slip earthquake, significant vertical deformation due to asthenosphere flow persists for many times tau and the magnitude of the vertical deformation is not necessarily enhanced by having a partially relaxing lower lithosphere. Previously announced in STAR as N83-13683

Cohen, S. C.↗

Sputtering as a Technique for Applying Tribological Coatings

Friction and wear induced mechanical failures may be controlled to extend the life of tribological components through the interposition of selected solid materials between contacting surfaces. Thin solid films of soft and hard materials are appropriate to lower friction and enhance the wear resistance of precision tribo-elements. Thin film coating technologies that have been developed may be used to interpose thin solid layers between contacting surfaces to reduce friction or to extend wear life. Film deposition technologies developed include chemical vapor deposition (CVD) and physical vapor deposition (PVD). Sputter coating is a PVD process for thin film deposition. It allows close control of coating thickness and composition and is well-suited for the deposition of tribological coatings. In this paper, the utility of magnetron sputtering for the deposition of tribological coatings is examined. Experimental results obtained when coated bodies are tested under a number of tribological operating conditions are presented. It is shown that when film adhesion problems are resolved, sputtering can be an effective means of depositing tribological coatings.

S Ramalingam↗

Unstructured Cartesian/prismatic grid generation for complex geometries

The generation of a hybrid grid system for discretizing complex three dimensional (3D) geometries is described. The primary grid system is an unstructured Cartesian grid automatically generated using recursive cell subdivision. This grid system is sufficient for computing Euler solutions about extremely complex 3D geometries. A secondary grid system, using triangular-prismatic elements, may be added for resolving the boundary layer region of viscous flows near surfaces of solid bodies. This paper describes the grid generation processes used to generate each grid type. Several example grids are shown, demonstrating the ability of the method to discretize complex geometries, with very little pre-processing required by the user.

Karman, Steve L., Jr.↗

Simulation of Dynamic Crushing Behavior of Simulated Carbon/Epoxy Stanchions Using MAT58 and MAT213 in LS-DYNA

Carbon fiber-reinforced polymer (CFRP) composites have shown promise as a material for structures designed to absorb energy in crush-style impact loading. In complementing the lightweight and tailorable characteristics of CFRP structures, the benefits of impact energy absorption are particularly interesting in aerospace vehicles that require excellent crashworthiness performance. However, simulating the behavior of CFRP structures in dynamic, crushing-style impact loading is challenging because of the many damage and failure modes that are essential to capture in the model. For the present investigation, CFRP stanchions were tested using a crash sled experimental test rig. The stanchions were designed by the Composite Materials Handbook-17 (CMH-17) Crashworthiness Working Group for the purpose of comparing experimental crash sled tests to corresponding simulations of the tests. The stanchions are C-channel shaped and represent a geometry common in the interior of aerospace vehicle structures. Explicit simulations in LSDYNA were performed using a well-established composite material model (MAT58) and a next-generation material model (MAT213). Fully integrated shell elements were utilized instead of 3D solid elements to limit computation time. In all models, each of the sixteen plies were represented with individual layers of shell elements with tiebreak contact between each adjacent ply pair to simulate interlaminar fracture. Simulating the crushing of the stanchions occurred in two phases. First, the material models were calibrated using flat specimens that were manufactured with the same layup as the stanchions. While the original goal was calibration of material-related properties, meshdependent behavior was observed in simulations with either material model, and an unstructured mesh was selected to remediate undesirable mesh-dependent failure modes. Additionally, for both the MAT58 and MAT213 models, it was found that either the crush force or the failure mode could be modeled accurately, but no set of parameters could be identified to attain both results in the same model. Once satisfactory calibration was achieved, the same material parameters were applied to the stanchion crush simulations. The stanchion simulations showed that MAT213 more accurately predicted the experimentally determined crush force, and both material models predicted key aspects of the experimentally observed failure modes.

Composite Materials↗

Opaque Mineral Assemblages at Chondrule Boundaries in the Vigarano CV Chondrite: Evidence for Gas-Solid Reactions Following Chondrule Formation

Recent studies of opaque minerals in primitive ordinary chondrites suggest that metal grains exposed at chondrule boundaries were corroded when volatile elements recondensed after the transient heating event responsible for chondrule formation. Metal grains at chondrule boundaries in the Bishunpur (LL3.1) chondrite are rimmed by troilite and fayalite. If these layers formed by gas solid reaction, then the composition of the corrosion products can provide information on the chondrule formation environment. Given the broad similarities among chondrules from different chondrite groups, similar scale layers should occur on chondrules in other primitive meteorite groups. Here I report on metal grains at chondrule boundaries in Vigarano (CV3).

Lauretta, Dante S.↗

Finite Element Models for Electron Beam Freeform Fabrication Process

Electron beam freeform fabrication (EBF3) is a member of an emerging class of direct manufacturing processes known as solid freeform fabrication (SFF); another member of the class is the laser deposition process. Successful application of the EBF3 process requires precise control of a number of process parameters such as the EB power, speed, and metal feed rate in order to ensure thermal management; good fusion between the substrate and the first layer and between successive layers; minimize part distortion and residual stresses; and control the microstructure of the finished product. This is the only effort thus far that has addressed computer simulation of the EBF3 process. The models developed in this effort can assist in reducing the number of trials in the laboratory or on the shop floor while making high-quality parts. With some modifications, their use can be further extended to the simulation of laser, TIG (tungsten inert gas), and other deposition processes. A solid mechanics-based finite element code, ABAQUS, was chosen as the primary engine in developing these models whereas a computational fluid dynamics (CFD) code, Fluent, was used in a support role. Several innovative concepts were developed, some of which are highlighted below. These concepts were implemented in a number of new computer models either in the form of stand-alone programs or as user subroutines for ABAQUS and Fluent codes. A database of thermo-physical, mechanical, fluid, and metallurgical properties of stainless steel 304 was developed. Computing models for Gaussian and raster modes of the electron beam heat input were developed. Also, new schemes were devised to account for the heat sink effect during the deposition process. These innovations, and others, lead to improved models for thermal management and prediction of transient/residual stresses and distortions. Two approaches for the prediction of microstructure were pursued. The first was an empirical approach involving the computation of thermal gradient, solidification rate, and velocity (G,R,V) coupled with the use of a solidification map that should be known a priori. The second approach relies completely on computer simulation. For this purpose a criterion for the prediction of morphology was proposed, which was combined with three alternative models for the prediction of microstructure; one based on solidification kinetics, the second on phase diagram, and the third on differential scanning calorimetry data. The last was found to be the simplest and the most versatile; it can be used with multicomponent alloys and rapid solidification without any additional difficulty. For the purpose of (limited) experimental validation, finite element models developed in this effort were applied to three different shapes made of stainless steel 304 material, designed expressly for this effort with an increasing level of complexity. These finite element models require large computation time, especially when applied to deposits with multiple adjacent beads and layers. This problem can be overcome, to some extent, by the use of fast, multi-core computers. Also, due to their numerical nature coupled with the fact that solid mechanics- based models are being used to represent the material behavior in liquid and vapor phases as well, the models have some inherent approximations that become more pronounced when dealing with multi-bead and multi-layer deposits.

Chandra, Umesh↗

Secondary flow spanwise deviation model for the stators of NASA middle compressor stages

A model of the spanwise variation of deviation for stator blades is presented. Deviation is defined as the difference between the passage mean flow angle and the metal angle at the outlet of a blade element of an axial compressor stage. The variation of deviation is taken as the difference above or below that predicted by blade element, (i.e., two-dimensional) theory at any spanwise location. The variation of deviation is dependent upon the blade camber, solidity and inlet boundary layer thickness at the hub or tip end-wall, and the blade channel aspect ratio. If these parameters are known or can be calculated, the model provides a reasonable approximation of the spanwise variation of deviation for most compressor middle stage stators operating at subsonic inlet Mach numbers.

Roberts, W. B.↗

The Liquid Annular Reactor System (LARS) propulsion

A concept for very high specific impulse (greater than 2000 seconds) direct nuclear propulsion is described. The concept, termed the liquid annular reactor system (LARS), uses liquid nuclear fuel elements to heat hydrogen propellant to very high temperatures (approximately 6000 K). Operating pressure is moderate (approximately 10 atm), with the result that the outlet hydrogen is virtually 100 percent dissociated to monatomic H. The molten fuel is contained in a solid container of its own material, which is rotated to stabilize the liquid layer by centripetal force. LARS reactor designs are described, together with neutronic and thermal-hydraulic analyses. Power levels are on the order of 200 megawatts. Typically, LARS designs use seven rotating fuel elements, are beryllium moderated, and have critical radii of approximately 100 cm (core L/D approximately equal to 1.5).

Powell, James↗

Development of an improved model for runback water on aircraft surfaces

A computer simulation for 'running wet' and evaporative aircraft anti-icing systems is developed. The model is based on the analysis of the liquid water film which forms in the regions of direct impingement and, then, breaks up near the impingement limits into rivulets. The wetness factor distribution resulting from the film breakup and the rivulet configuration on the surface are predicted using a stability analysis theory and the laws of mass energy conservation. The solid structure is modeled as a multiple layer wall. The anti-icing system modeled is of the thermal type utilizing hot air and/or electrical heating elements embedded within the wall layers. Experimental observations revealing some of the basic physics of the water flow on the surface are presented. Detailed qualitative documentation of the tests are given. Several numerical examples are considered, and the effect of some of the involved parameters on the system performance are investigated.

Al-Khalil, Kamel M.↗

Computational models for the nonlinear analysis of reinforced concrete plates

A finite element computational model for the nonlinear analysis of reinforced concrete solid, stiffened and cellular plates is briefly outlined. Typically, Mindlin elements are used to model the plates whereas eccentric Timoshenko elements are adopted to represent the beams. The layering technique, common in the analysis of reinforced concrete flexural systems, is incorporated in the model. The proposed model provides an inexpensive and reasonably accurate approach which can be extended for use with voided plates.

Hinton, E.↗