SEARCH · Search NASA
Results for “Finite Element Modeling”
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.
Analysis of Cavern and Well Stability at the Bryan Mound SPR Site Using the M–D Salt Creep Model
This report presents computational analyses that simulate the structural response of caverns at the Strategic Petroleum Reserve Bryan Mound site. The cavern field comprises 20 caverns. Five caverns (1, 2, 4, and 5; 3 was later plugged and abandoned) were acquired from industry and have unusual shapes and a history dating back to 1946. The other 16 caverns (101-116) were leached according to SPR standards in the mid-1980s and have tall cylindrical shapes. The history of the caverns and their shapes are simulated in a 3-D geomechanics model of the site that predicts deformations, strains, and stresses. Historical wellhead pressures are used to calculate cavern pressures up through July 2016. Because of the extent of heterogeneous creep behavior observed throughout the Bryan Mound site, a set of cavern-specific creep coefficients was developed to produce better matches with measured cavern closure and surface subsidence. For this new implementation of the model, there are two significant advances: the use of the multimechanism deformation (M-D) salt creep model to evaluate both steady-state and transient salt creep; and the creation of finite element mesh geometries for the caverns that nearly exactly match the geometries obtained through sonar measurements. The results of the finite element model are interpreted to provide information on the current and future status of subsidence, well integrity, cavern stability, and drawdown availability.
Vibration analysis of a 1/15 scale dynamic model of a space shuttle configuration
The natural frequencies and mode shapes of a 1/15 scale space shuttle dynamics model are analytically determined. The model, a parallel beam type structure with delta wings, is dynamically representative of the stiffness and mass properties of an early space shuttle design. Important characteristics of the model are elastic interfaces with adjustable spring rates. Normal mode computations are made using the finite element modeling technique as implemented in the NASTRAN (NASA Structural Analysis) computer program. The feasibility of neglecting elastic deformations in the lower modes was investigated using a rigid body model. Using NASTRAN, natural frequencies and mode shapes were first calculated for the booster fuselage, orbiter fuselage, and both delta wings in a free-free condition. Next, the fuselages were connected for various spring rates. Then the wings were attached to each fuselage, and the booster and orbiter were analyzed as separate airplanes. Finally, the two airplanes were elastically joined, and the complete model was analyzed.
Space shuttle lift-off dynamic model
Previously developed dynamic models for the calculation of lift-off dynamic response of the space shuttle vehicle can handle only response of the vehicle with eight holddown arms attached or with all holddown arms detached. The new model developed in the referenced report takes into account the transition period between holddown and lift-off by giving the model the ability to vary holddown point separation as a function of vehicle flexible body and rigid body motion. This report documents a study made to verify the new model's capability to simulate vehicle response at lift-off. To do this, a finite element model of the skirt is made and coupled to the free-free modes of the vehicle and cantilevered modes calculated and compared with the previously developed model. Results indicate that the new model will be able to predict accurate vehicle loads.
Generic model of a large flexible space structure for control concept evaluation
The objective of the paper is to propose a generic model on the shape and attitude control of a representative space structure. Reduction of the finite element model is based on modal identities and modal cost analysis. The model may help to: (a) suggest improved methods for designing reduced order controllers for large elastic spacecraft, (b) enhance coherence between the frequency-response and time-response multivariable applications, (c) examine computational efficiency of the multivariable control algorithms, and (d) suggest other desirable theoretical and numerical research areas.
Finite element analysis of the end notched flexure specimen for measuring Mode II fracture toughness
The paper presents a finite element analysis of the end-notched flexure (ENF) test specimen for Mode II interlaminar fracture testing of composite materials. Virtual crack closure and compliance techniques employed to calculate strain energy release rates from linear elastic two-dimensional analysis indicate that the ENF specimen is a pure Mode II fracture test within the constraints of small deflection theory. Furthermore, the ENF fracture specimen is shown to be relatively insensitive to process-induced cracks, offset from the laminate midplane. Frictional effects are investigated by including the contact problem in the finite element model. A parametric study investigating the influence of delamination length, span, thickness, and material properties assessed the accuracy of beam theory expressions for compliance and strain energy release rate, GII. Finite element results indicate that data reduction schemes based upon beam theory underestimate GII by approximately 20-40 percent for typical unidirectional graphite fiber composite test specimen geometries. Consequently, an improved data reduction scheme is proposed.
Analysis and test of a 16-foot radial rib reflector developmental model
Analytical and experimental modal tests were performed to determine the vibrational characteristics of a 16-foot diameter radial rib reflector model. Single rib analyses and experimental tests provided preliminary information relating to the reflector. A finite element model predicted mode shapes and frequencies of the reflector. The analyses correlated well with the experimental tests, verifying the modeling method used. The results indicate that five related, characteristic mode shapes form a group. The frequencies of the modes are determined by the relative phase of the radial ribs.
Mechanical Characterization and Micromechanical Modeling of Woven Carbon/Copper Composites
The present investigation examines the in-plane mechanical behavior of a particular woven metal matrix composite (MMC); 8-harness (8H) satin carbon/copper (C/Cu). This is accomplished via mechanical testing as well as micromechanical modeling. While the literature is replete with experimental and modeling efforts for woven and braided polymer matrix composites, little work has been done on woven and braided MMC's. Thus, the development and understanding of woven MMC's is at an early stage. 8H satin C/Cu owes its existence to the high thermal conductivity of copper and low density and thermal expansion of carbon fibers. It is a candidate material for high heat flux applications, such as space power radiator panels. The experimental portion of this investigation consists of monotonic and cyclic tension, compression, and Iosipescu shear tests, as well as combined tension-compression tests. Tests were performed on composite specimens with three copper matrix alloy types: pure Cu, Cu-0.5 weight percent Ti (Cu-Ti), and Cu-0.7 weight percent Cr (Cu-Cr). The small alloying additions are present to promote fiber/matrix interfacial bonding. The analytical modeling effort utilizes an approach in which a local micromechanical model is embedded in a global micromechanical model. This approach differs from previously developed analytical models for woven composites in that a true repeating unit cell is analyzed. However, unlike finite element modeling of woven composites, the geometry is sufficiently idealized to allow efficient geometric discretization and efficient execution.
Computational Models of the Eye and their Applications in Long Duration Space Flight
Astronauts are exposed to cephalad fluid shift, increased carbon dioxide levels and other environmental factors during space flight. As a result of these conditions, it is believed that they are at risk of developing increased intracranial pressure (ICP) and intraocular pressure (IOP), which in turn may cause papilledema and other disorders of the eye that can lead to temporary or permanent changes in vision. However, the mechanisms behind this risk are not fully understood. Ground analog and flight studies pose challenges because there are limited non-invasive methods that can be used to study the eye and intracranial space. Therefore it is proposed that computational models can be applied to help address this gap by providing a low cost method for studying the effects of IOP, ICP and various properties of the eye on these diseases. The information presented by the authors provides a summary of several models found in literature that could potentially be augmented and applied to inform research. Specifically, finite element models of the optic nerve head, sclera and other structures of the eye can be readily adapted as potential building blocks. These models may also be integrated with a brain/cerebrospinal fluid (CSF) model which will take into account the interaction between the CSF fluid and its pressure on the optic nerve. This integration can enable the study of the effects of microgravity on the interaction between the vasculature system and CSF system and can determine the effects of these changes on the optic nerve, and in turn the eye. Ultimately, it can help pinpoint the influences of long-term exposure to microgravity on vision and inform the future research into countermeasure development. In addition to spaceflight, these models can provide deeper understanding of the mechanisms of glaucoma, papilledema and other eye disorders observed in terrestrial conditions.
Shape Memory Alloy Material Maturation (SMA)
The Shape Memory Alloy (SMA) project will evaluate, process, and refine 3 SMA materials that have the potential to provide the properties needed to enable various technologies for human exploration missions on the Moon and Mars. The work will focus on the development and identification of a new SMA alloy targeted at enabling compliant rover tires capable of the extreme demands (high loads, faster drive speeds, extreme cold temperatures and exceedingly long mission durations). The necessary work will include the development of key material testing capabilities that will enable material properties development from 40K – 400K for all properties needed to qualify space flight hardware (Tensile, Fatigue/Durability, etc.). The SMA wire properties will be input into a database and model, ultimately to disseminate material property design data (for 1 alloy only) as parameters for finite element modeling of a tire structure.
Exploiting symmetries in the modeling and analysis of tires
A simple and efficient computational strategy for reducing both the size of a tire model and the cost of the analysis of tires in the presence of symmetry-breaking conditions (unsymmetry in the tire material, geometry, or loading) is presented. The strategy is based on approximating the unsymmetric response of the tire with a linear combination of symmetric and antisymmetric global approximation vectors (or modes). Details are presented for the three main elements of the computational strategy, which include: use of special three-field mixed finite-element models, use of operator splitting, and substantial reduction in the number of degrees of freedom. The proposed computational stategy is applied to three quasi-symmetric problems of tires: linear analysis of anisotropic tires, through use of semianalytic finite elements, nonlinear analysis of anisotropic tires through use of two-dimensional shell finite elements, and nonlinear analysis of orthotropic tires subjected to unsymmetric loading. Three basic types of symmetry (and their combinations) exhibited by the tire response are identified.
Modeling of the surface static displacements and fault plane slip for the 1979 Imperial Valley earthquake
Three-dimensional finite element modeling techniques are used to synthesize geodetic and seismological results for 1979 Imperial Valley earthquake. The strategy pursued consists of two principal steps. In the first step, the seismologically-derived coseismic fault slip is taken as a function of position in the fault plane and is applied directly to the three-dimensional dislocation model. In the second step, a physical model of stresses and constitutive parameters is perturbed so as to reproduce the observed fault slip. Hence, the principal features of the coseismic slip pattern are explained by a stress-driven fault model in which: (1) a spatially unresolved asperity is found equivalent to a stress drop of 18 MPa averaged over an area of 15 sq km, and (2) driving stress is essentially absent on the fault segment overlapping the 1940 earthquake rupture zone.
Uncertainty Models for the Hybrid Parametric Variation Method of Uncertainty Quantification; Analysis
There is some level of uncertainty in every finite element model (FEM), which flows to a level of uncertainty in predicted results. The purpose of uncertainty quantification (UQ) is to provide statistical bounds on prediction accuracy based on model uncertainty. This is distinct from model updating, which attempts to modify models to improve their accuracy. UQ does not improve the accuracy of models, but accepts that the models are inherently inaccurate and attempts to quantify the impact of that inaccuracy on predicted results. Previously, an alternate method for UQ, called the Hybrid Parametric Variation (HPV) method, was applied to Space Launch System (SLS) Hurty/Craig-Bampton (HCB) components to predict system-level statistics for launch vehicle attitude control transfer functions and core stage section loads due to buffet. The HPV method combines a parametric variation of the HCB fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method that randomly varies the HCB mass and stiffness matrices as Wishart random matrix distributions using random matrix theory (RMT). Alternatively, the most common method for modeling uncertainty in the structural dynamics community is a parametric approach, which varies physical parameters in the model. However, there are several disadvantages associated with the parametric method. Determining a reduced set of parameters that have a significant impact on the system response can be time consuming, and the selected parameter probability distributions are rarely reliably known. Therefore, in practice, the parameters are surrogates for the actual errors, and the link to parameter uncertainty is unknown. Another major drawback is that the uncertainty that can be represented is limited to the form of the nominal FEM. It is the experience of the authors that based on numerous aerospace programs, almost all FEM errors are in form rather than parameter values. This hypothesis is supported by the observation of the authors that it is almost never possible to ‘tune’ a FEM to match modal test results by only modifying model parameters. Model-form uncertainty cannot be directly represented by FEM input parameters nor included in a parametric approach. However, model-form uncertainty can be modeled using RMT, where a probability distribution is developed for the matrix ensemble of interest. The major advantage of the NPV method is that it covers errors in model form. The HPV method anchors uncertainty at the HCB component level to component modal test results by matching the HCB and test modes based on mode descriptions or other methods, and then applying differing levels of frequency variation. The specific variations depend on the confidence to which a component FEM has been validated through modal testing. The NPV method is layered on the frequency variation to match modal test self-orthogonality and cross-orthogonality (XO) results. Once the component uncertainty models are identified, they are assembled, and the uncertainty is propagated to the system level using a Monte Carlo (MC) analysis approach that generates statistics for system-level predictions This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The second example is an application to a component that is design specific to the SLS. Based on this work and other assessments, the HPV method provides another tool to the toolset used for complex system UQ analysis. From experience gathered to date using the HPV method, additional design specific applications must be investigated to provide further confidence in the validity of the HPV method of UQ analysis.
Numerical modelling of electrostatic and electromagnetic positioners for microgravity containerless processing
Numerical, three-dimensional, finite-difference and finite-element models of electrostatic and electromagnetic positioners are developed. These models allow a better understanding of the effects of the positioner on the sample during processing in microgravity and provide an engineering tool for optimizing the design of these positioners.
Development and Application of Benchmark Examples for Mixed-Mode I/II Quasi-Static Delamination Propagation Predictions
The development of benchmark examples for quasi-static delamination propagation prediction is presented and demonstrated for a commercial code. The examples are based on finite element models of the Mixed-Mode Bending (MMB) specimen. The examples are independent of the analysis software used and allow the assessment of the automated delamination propagation prediction capability in commercial finite element codes based on the virtual crack closure technique (VCCT). First, quasi-static benchmark examples were created for the specimen. Second, starting from an initially straight front, the delamination was allowed to propagate under quasi-static loading. Third, the load-displacement relationship from a propagation analysis and the benchmark results were compared, and good agreement could be achieved by selecting the appropriate input parameters. Good agreement between the results obtained from the automated propagation analysis and the benchmark results could be achieved by selecting input parameters that had previously been determined during analyses of mode I Double Cantilever Beam and mode II End Notched Flexure specimens. The benchmarking procedure proved valuable by highlighting the issues associated with choosing the input parameters of the particular implementation. Overall the results are encouraging, but further assessment for mixed-mode delamination fatigue onset and growth is required.
Two-Dimensional Nonlinear Finite Element Analysis of CMC Microstructures
A research program has been developed to quantify the effects of the microstructure of a woven ceramic matrix composite and its variability on the effective properties and response of the material. In order to characterize and quantify the variations in the microstructure of a five harness satin weave, chemical vapor infiltrated (CVI) SiC/SiC composite material, specimens were serially sectioned and polished to capture images that detailed the fiber tows, matrix, and porosity. Open source quantitative image analysis tools were then used to isolate the constituents, from which two dimensional finite element models were generated which approximated the actual specimen section geometry. A simplified elastic-plastic model, wherein all stress above yield is redistributed to lower stress regions, is used to approximate the progressive damage behavior for each of the composite constituents. Finite element analyses under in-plane tensile loading were performed to examine how the variability in the local microstructure affected the macroscopic stress-strain response of the material as well as the local initiation and progression of damage. The macroscopic stress-strain response appeared to be minimally affected by the variation in local microstructure, but the locations where damage initiated and propagated appeared to be linked to specific aspects of the local microstructure.
JWST ISIM Distortion Analysis Challenge
Very tight distortion requirements are imposed on the JWST's ISM structure due to the sensitivity of the telescope's mirror segment and science instrument positioning. The ISIM structure is a three dimensional truss with asymmetric gusseting and metal fittings. One of the primary challenges for ISIM's analysis team is predicting the thermal distortion of the structure both from the bulk cooldown from ambient to cryo, and the smaller temperature changes within the cryogenic operating environment. As a first cut to estimate thermal distortions, a finite element model of bar elements was created. Elements representing joint areas and metal fittings use effective properties that match the behavior of the stack-up of the composite tube, gusset and adhesive under mechanical and thermal loads. These properties were derived by matching tip deflections of a solid model simplified T-joint. Because of the structure s asymmetric gusseting, this effective property model is a first attempt at predicting rotations that cannot be captured with a smeared CTE approach. In addition to the finite element analysis, several first order calculations have been performed to gauge the feasibility of the material design. Because of the stringent thermal distortion requirements at cryogenic temperatures, a composite tube material with near zero or negative CTE is required. A preliminary hand analysis of the contribution of the various components along the distortion path between FGS and the other instruments, neglecting second order effects were examined. A plot of bounding tube longitudinal and transverse CTEs for thermal stability requirements was generated to help determine the feasibility of meeting these requirements. This analysis is a work in progress en route to a large degree of freedom hi-fidelity FEA model for distortion analysis. Methods of model reduction, such as superelements, are currently being investigated.
Mixed models and reduction method for dynamic analysis of anisotropic shells
A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.