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 1,135 records · Page 63

Fracture mechanics analysis for various fiber/matrix interface loadings

Fiber/matrix (F/M) cracking was analyzed to provide better understanding and guidance in developing F/M interface fracture toughness tests. Two configurations, corresponding to F/M cracking at a broken fiber and at the free edge, were investigated. The effects of mechanical loading, thermal cooldown, and friction were investigated. Each configuration was analyzed for two loadings: longitudinal and normal to the fiber. A nonlinear finite element analysis was performed to model friction and slip at the F/M interface. A new procedure for fitting a square-root singularity to calculated stresses was developed to determine stress intensity factors (K sub I and K sub II) for a bimaterial interface crack. For the case of F/M cracking at a broken fiber with longitudinal loading, crack tip conditions were strongly influenced by interface friction. As a result, an F/M interface toughness test based on this case was not recommended because nonlinear data analysis methods would be required. For the free edge crack configuration, both mechanical and thermal loading caused crack opening, thereby avoiding frictional effects. An F/M interface toughness test based on this configuration would provide data for K(sub I)/K(sub II) ratios of about 0.7 and 1.6 for fiber and radial normal loading, respectively. However, thermal effects must be accounted for in the data analysis.

Naik, R. A.↗

Fracture mechanics analysis for various fiber/matrix interface loadings

Fiber/matrix (F/M) cracking was analyzed to provide better understanding and guidance in developing F/M interface fracture toughness tests. Two configurations, corresponding to F/M cracking at a broken fiber and at the free edge, were investigated. The effects of mechanical loading, thermal cooldown, and friction were investigated. Each configuration was analyzed for two loadings: longitudinal and normal to the fiber. A nonlinear finite element analysis was performed to model friction and slip at the F/M interface. A new procedure for fitting a square-root singularity to calculated stresses was developed to determine stress intensity factors (K sub I and K sub II) for a bimaterial interface crack. For the case of F/M cracking at a broken fiber with longitudinal loading, crack tip conditions were strongly influenced by interface friction. As a result, an F/M interface toughness test based on this case was not recommended because nonlinear data analysis methods would be required. For the free edge crack configuration, both mechanical and thermal loading caused crack opening, theory avoiding fractional effects. A F/M interface toughness test based on this configuration would provide data for K(sub I/K(sub II) ratios of about 0.7 and 1.6 for fiber and radial normal loading, respectively. However, thermal effects must be accounted for in the data analysis.

Naik, Rajiv A.↗

Analysis of the shearout failure mode in composite bolted joints

A semi-empirical shearout strength model has been formulated for the analysis of composite bolted joints with allowance for the effects of joint geometry. The model employs a polynomial stress function in conjunction with a point stress failure criterion to predict strength as a function of fastener size, edge distance, and half spacing. The stress function is obtained by two-dimensional plane-stress finite element analysis using quadrilateral elements with orthotropic material properties. Comparison of experimentally determined shearout strength data with model predicted failures has substantiated the accuracy of the model.

Wilson, D. W.↗

Finite-Element Modeling For Structural Analysis

Report presents study of finite-element mathematical modeling as used in analyzing stresses and strains at joints between thin, shell-like components (e.g., ducts) and thicker components (e.g., flanges or engine blocks). First approach uses global/local model to evaluate system. Provides correct total response and correct representation of stresses away from any discontinuities. Second approach involves development of special transition finite elements to model transitions between shells and thicker structural components.

Min, J. B.↗

Ceramics Analysis and Reliability Evaluation of Structures (CARES). Users and programmers manual

This manual describes how to use the Ceramics Analysis and Reliability Evaluation of Structures (CARES) computer program. The primary function of the code is to calculate the fast fracture reliability or failure probability of macroscopically isotropic ceramic components. These components may be subjected to complex thermomechanical loadings, such as those found in heat engine applications. The program uses results from MSC/NASTRAN or ANSYS finite element analysis programs to evaluate component reliability due to inherent surface and/or volume type flaws. CARES utilizes the Batdorf model and the two-parameter Weibull cumulative distribution function to describe the effect of multiaxial stress states on material strength. The principle of independent action (PIA) and the Weibull normal stress averaging models are also included. Weibull material strength parameters, the Batdorf crack density coefficient, and other related statistical quantities are estimated from four-point bend bar or unifrom uniaxial tensile specimen fracture strength data. Parameter estimation can be performed for single or multiple failure modes by using the least-square analysis or the maximum likelihood method. Kolmogorov-Smirnov and Anderson-Darling goodness-of-fit tests, ninety percent confidence intervals on the Weibull parameters, and Kanofsky-Srinivasan ninety percent confidence band values are also provided. The probabilistic fast-fracture theories used in CARES, along with the input and output for CARES, are described. Example problems to demonstrate various feature of the program are also included. This manual describes the MSC/NASTRAN version of the CARES program.

Nemeth, Noel N.↗

Modeling the Effective Elasticity of Anisotropic Porous Materials

The development and optimization of composite materials designed for thermal protection of NASA’s spacecraft require understanding their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the microscale on realistic representations of these composites. The digital microstructures are generated either synthetically or through X-ray micro-computed tomography reconstructions. One of the main challenges in the prediction of the structural response of heatshield materials is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the microscale. These are driven by the fiber layout within the microstructure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled using a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA) that was presented at the 10th Ablation Workshop. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases, and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Elasticity↗

Elastic-plastic analysis of AS4/PEEK composite laminate using a one-parameter plasticity model

A one-parameter plasticity model was shown to adequately describe the plastic deformation of AS4/PEEK (APC-2) unidirectional thermoplastic composite. This model was verified further for unidirectional and laminated composite panels with and without a hole. The elastic-plastic stress-strain relations of coupon specimens were measured and compared with those predicted by the finite element analysis using the one-parameter plasticity model. The results show that the one-parameter plasticity model is suitable for the analysis of elastic-plastic deformation of AS4/PEEK composite laminates.

Sun, C. T.↗

The Inclusion of In-Plane Stresses in Delamination Criteria

A study of delamination failure was conducted with emphasis on delamination criteria. Evidence is presented which supports the inclusion of the in-plane stresses in addition to the interlaminar stress terms in delamination criteria. The delamination is characterized as the failure of a resin rich region in between ply sets. The entire six component stress state in this resin layer is calculated through a finite element analysis, averaged over a dimension of 1.75 ply thicknesses, and used in a Modified von Mises Delamination Criterion. This criterion builds onto previous criteria by including all six stress components in the interply resin layer. The MVMDC shows good correlation to experimental data. The results show that the treatment of delamination as the failure of a finite interply resin layer is a valid method and that the MVMDC, considering the full stress state, accurately indicates delamination for different laminate families.

Fenske, Matthew T.↗

Multidisciplinary Tool for Systems Analysis of Planetary Entry, Descent, and Landing

Systems analysis of a planetary entry (SAPE), descent, and landing (EDL) is a multidisciplinary activity in nature. SAPE improves the performance of the systems analysis team by automating and streamlining the process, and this improvement can reduce the errors that stem from manual data transfer among discipline experts. SAPE is a multidisciplinary tool for systems analysis of planetary EDL for Venus, Earth, Mars, Jupiter, Saturn, Uranus, Neptune, and Titan. It performs EDL systems analysis for any planet, operates cross-platform (i.e., Windows, Mac, and Linux operating systems), uses existing software components and open-source software to avoid software licensing issues, performs low-fidelity systems analysis in one hour on a computer that is comparable to an average laptop, and keeps discipline experts in the analysis loop. SAPE uses Python, a platform-independent, open-source language, for integration and for the user interface. Development has relied heavily on the object-oriented programming capabilities that are available in Python. Modules are provided to interface with commercial and government off-the-shelf software components (e.g., thermal protection systems and finite-element analysis). SAPE currently includes the following analysis modules: geometry, trajectory, aerodynamics, aerothermal, thermal protection system, and interface for structural sizing.

Samareh, Jamshid A.↗

Flutter Analysis for Turbomachinery Using Volterra Series

The objective of this paper is to describe an accurate and efficient reduced order modeling method for aeroelastic (AE) analysis and for determining the flutter boundary. Without losing accuracy, we develop a reduced order model based on the Volterra series to achieve significant savings in computational cost. The aerodynamic force is provided by a high-fidelity solution from the Reynolds-averaged Navier-Stokes (RANS) equations; the structural mode shapes are determined from the finite element analysis. The fluid-structure coupling is then modeled by the state-space formulation with the structural displacement as input and the aerodynamic force as output, which in turn acts as an external force to the aeroelastic displacement equation for providing the structural deformation. NASA's rotor 67 blade is used to study its aeroelastic characteristics under the designated operating condition. First, the CFD results are validated against measured data available for the steady state condition. Then, the accuracy of the developed reduced order model is compared with the full-order solutions. Finally the aeroelastic solutions of the blade are computed and a flutter boundary is identified, suggesting that the rotor, with the material property chosen for the study, is structurally stable at the operating condition, free of encountering flutter.

computational fluid dynamics↗

Static and dynamic structural-sensitivity derivative calculations in the finite-element-based Engineering Analysis Language (EAL) system

The implementation of static and dynamic structural-sensitivity derivative calculations in a general purpose, finite-element computer program denoted the Engineering Analysis Language (EAL) System is described. Derivatives are calculated with respect to structural parameters, specifically, member sectional properties including thicknesses, cross-sectional areas, and moments of inertia. Derivatives are obtained for displacements, stresses, vibration frequencies and mode shapes, and buckling loads and mode shapes. Three methods for calculating derivatives are implemented (analytical, semianalytical, and finite differences), and comparisons of computer time and accuracy are made. Results are presented for four examples: a swept wing, a box beam, a stiffened cylinder with a cutout, and a space radiometer-antenna truss.

Camarda, C. J.↗

Compliance calibration of specimens used in the R-curve practice

The compliance calibrations for the compact (CS) and crack-line-wedge-loaded (CLWL) specimens have been determined by experimental measurements and by boundary-collocation analysis. The CS and CLWL specimen configurations were modeled more accurately than those used in previous analytical investigations. Polynomial expressions for the compliance at various stations along the crack line for CS and CLWL specimens are presented. The compliance calibrations for the center-crack tension (CCT) specimen have been determined theoretically by boundary-collocation and finite-element analysis. The calculated compliance values for the CCT specimen are compared with values obtained from the Irwin-Westergaard expression and from a modification to the Irwin-Westergaard expression proposed by Eftis and Liebowitz. The Eftis-Liebowitz expression was found to be in good agreement (plus or minus 2 percent) with both analyses for crack aspect ratios up to 0.8 and for gage half-span to specimen width ratios up to 0.5.

Mccabe, D. E.↗

Progressive Damage Analysis of Post-Buckled Stiffened Panels under Static Compressive Loading

The validation of progressive damage analysis methods requires experimental data that captures the initiation of damage, its evolution with continued loading, and the morphology of damage throughout its evolution. In support of the NASA Advanced Composites Project, several validation specimens of increasing complexity have been tested and characterized in order to provide the data required to assess the capabilities of several progressive damage analysis methods for fiber-reinforced composite materials. The analytical capabilities of interest for this validation program include networks of interacting matrix cracks and delaminations in post-buckled stiffened structures. A series of single stringer compression specimens was tested and characterized at NASA Langley Research Center. The experimental data from these tests are compared with progressive damage finite element analysis results. Test/analysis comparisons are made in terms of load/displacement histories, damage morphology, and damage size versus load metrics.

Frank A Leone↗

Residual Stresses Induced due to Curing of the Bulk Matrix in a Simplified Three-Dimensional (3D) Woven Repeating Unit Cell

It has been observed that 3D woven materials are highly susceptible to processes induced defects including tow misalignment, intratow and intertow cracking and/or voids. Adequate process modeling of 3D woven Polymer Matrix Composites (PMCs) is necessary to predict and estimate the effects of the manufacturing process on these defects. A preliminary step towards this is achieved with the commercial Finite Element Analysis (FEA) tool Abaqus complemented with written user subroutines to account for the effect of shrinkage and thermoelastic properties evolution as a function of the degree of cure of the bulk matrix phase in a 3D woven composite repeating unit cell (RUC). As a pathfinder, the process modeling framework is demonstrated on a simplified 3D woven RUC of an AS4/RIMR 135 system. Curing simulations are performed on the 3D woven RUC to observe the processing induced residual stresses. The results of the analysis show that high stress concentrations, present within the 3D intertow matrix bordering the binder tow, coincide with X-Ray computed tomography (X-Ray-CT) scan data depicting cracking in the same location

3D Woven↗

Fitness-for-Service Analysis of Reactor Components under Flexible Load-Following Operating Conditions

Conventional power-generation plants, including nuclear plants, have been traditionally designed to provide a steady baseload energy capacity, optimizing output efficiency while minimizing variable costs. However, the growing adoption of large-scale renewable energy-generation systems, which rely on intermittent sources such as solar and wind, has introduced more variability into the energy supply in interconnected electricity grids. As a result, the next generation of power plants needs to operate in what is known as the load-following mode, requiring flexible adjustments in electricity production to align with the energy demand on the grid. This transition from the steady baseload operation to load-following operating conditions can significantly increase the number of times various plant components are exposed to transient stresses. This increased thermo-mechanical cycling can lead to accelerated material degradation, thereby elevating the risk of premature failure of a component. It becomes imperative to conduct a comprehensive analysis of fatigue, creep-fatigue, and stress corrosion cracking life, to assess the resilience of the various engineering components under these flexible load-following operating conditions. This study aims to develop a comprehensive numerical model of a light-water reactor pressure vessel (RPV) to investigate its degradation under various operating scenarios. This coupled thermo-mechanical finite element analysis evaluated the stress response of the RPV caused by considering fluctuations in thermal and mechanical loads caused by the varying pressure and temperature occurring during the load-following operation. Critical locations on the RPV are subsequently identified based on the stress response. The stress intensity factors for the postulated flaws at those locations are then calculated, followed by an evaluation of the reactor's life in accordance with the ASME Boiler and Pressure Vessel Code Section XI. This comprehensive life assessment covers a number of transients expected during the flexible load-following operation, providing invaluable insights into the RPV's structural integrity. Moreover, the development methodology can be adapted to other reactor components, as well as components of conventional power stations that are affected by varying operating conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermal-Hydraulic Analysis of the AFIP-7 Irradiation Test in the Advanced Test Reactor - Oxide Growth Prediction and Correlation

Knowing the thickness of the oxide layer on the surface of aluminum fuel cladding is vitally important for predicting fuel temperature due to the low thermal conductivity of the oxide layer. Several correlation models for predicting oxide growth can be found in the literature. In previous research, the correlations were combined with heat transfer simulations in Abaqus, a finite element analysis code, to forecast the oxide growth. However, this approach requires heat transfer coefficients for modeling heat exchanges with the external flow field, and such coefficients were obtained through empirical equations. Since different empirical equations yield varying heat transfer coefficients, the cladding temperature and predicted oxide thickness both carry a high degree of uncertainty. Here, this research develops a new approach that integrates the fluid flow, fluid and solid heat transfer, and oxide growth correlation(s) into a single computational fluid dynamics model. We demonstrate this approach’s ability to predict oxide development on the AFIP-7 plates during two Advanced Test Reactor (ATR) irradiation cycles. The projected oxide thickness falls within the experimental measurements taken during post irradiation examination.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Mechanically fastened joints in woven fabric composites

Strength analysis for composite bolted joints involves the mating of a stress analysis with an appropriate mode specific failure criterion for each of the primary failure modes. The stress analysis and failure criteria are independent of each other and can be manipulated separately in order to optimize the strength analysis package formed by their coupling. Material properties tests were conducted on rubber toughened graphite-epoxy material to measure the basic strength and stiffness in the warp and fill directions and in shear. Test matrices are summarized for investigations of laminate configuration, stacking sequence, fastener diameter, edge distance, fastener half spacing, laminate thickness, and fastener torque. A three dimensional finite element analysis computer program was written and failure criteria for net tension, shearout, and bearing were determined.

Wilson, D. W.↗

Design of ceramic components with the NASA/CARES computer program

The ceramics analysis and reliability evaluation of structures (CARES) computer program is described. The primary function of the code is to calculate the fast-fracture reliability or failure probability of macro-scopically isotropic ceramic components. These components may be subjected to complex thermomechanical loadings, such as those found in heat engine applications. CARES uses results from MSC/NASTRAN or ANSYS finite-element analysis programs to evaluate how inherent surface and/or volume type flaws component reliability. CARES utilizes the Batdorf model and the two-parameter Weibull cumulative distribution function to describe the effects of multiaxial stress states on material strength. The principle of independent action (PIA) and the Weibull normal stress averaging models are also included. Weibull material strength parameters, the Batdorf crack density coefficient, and other related statistical quantities are estimated from four-point bend bar or uniform uniaxial tensile specimen fracture strength data. Parameter estimation can be performed for a single or multiple failure modes by using a least-squares analysis or a maximum likelihood method. Kolmogorov-Smirnov and Anderson-Darling goodness-to-fit-tests, 90 percent confidence intervals on the Weibull parameters, and Kanofsky-Srinivasan 90 percent confidence band values are also provided. Examples are provided to illustrate the various features of CARES.

Nemeth, Noel N.↗