Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Elements”

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 937 records · Page 52

Nuclear Materials Packaging, Transportation, and Systems Analysis Group Software Quality Assurance Plan: ANSYS Mechanical Finite Element Analysis Software Version 2023R1

ANSYS Inc. develops and markets engineering simulation software and services used in the aerospace, automotive, manufacturing, electronics, biomedical, energy, defense, and many other industries. ANSYS is dedicated to engineering simulation and is the world’s leading software provider. ANSYS was founded in 1970 and is headquartered in Canonsburg, Pennsylvania. ANSYS provides an engineering analysis tool combining structural, thermal, computational fluid dynamics, acoustic, and electromagnetic simulation capabilities. ANSYS has two main programs, which use the same solvers: (1) Mechanical APDL (ANSYS Design Parametric Language), a Fortran-based coding platform, and (2) ANSYS Workbench, which uses a graphical user interface to aid in finite element analysis implementation. This plan covers both APDL and Workbench. The ANSYS computer program is a large-scale, multipurpose finite element program that can be used to solve several classes of engineering analyses. The analysis capabilities of ANSYS include the ability to solve static and dynamic structural analyses, steady-state and transient heat transfer problems, mode-frequency and buckling eigenvalue problems, static or time-varying magnetic analyses, and various types of field and coupled-field applications. The program contains many special features that allow nonlinearities or secondary effects such as plasticity, large strain, hyperelasticity, creep, swelling, large deflections, contact, stress stiffening, temperature dependency, material anisotropy, and radiation to be included in the solution. As ANSYS has been developed, other special capabilities such as substructuring, submodeling, random vibration, kinetostatics, kinetodynamics, free convection fluid analysis, acoustics, magnetics, piezoelectrics, coupled-field analysis, and design optimization have been added to the program. These capabilities contribute further to making ANSYS a multipurpose analysis tool for varied engineering disciplines. The ANSYS program has been in commercial use for over 50 years and has been used extensively in the aerospace, automotive, construction, electronic, energy services, manufacturing, nuclear, plastics, oil, and steel industries. Additionally, many consulting firms and hundreds of universities have used ANSYS for analysis, research, and educational purposes. ANSYS is recognized worldwide as one of the most widely used and capable programs of its type. Ansys design analysis software is the first created within a quality system with ISO 9001 certification, the internationally accepted quality standard. Product development, testing, maintenance and support processes also meet the United States Nuclear Regulatory Commission's quality requirements, as they have for nearly four decades. The Quality Assurance Service Agreement is suitable for the customers working in the nuclear industry who need to meet specific federal regulations including 10CRF50 Appendix B and provisions of 10CFR21. ANSYS has retained its original International Organization for Standardization (ISO) 9001 accreditation certificate since1995-05-04, It’s current certificate is valid until 2027-05-29.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Implementing Geometric Surface Imperfections into Sandwich Composite Cylinder Finite Element Method Models

The buckling responses of certain cylindrical shell structures are extremely sensitive to geometric imperfections. The NASA Engineering and Safety Center (NESC) Shell Buckling Knockdown Factor Project (SBKF) is conducting research to develop analysis-based buckling design recommendations. Experiments are used to verify the analysis-based factors, but the sensitivity of the test articles to geometric imperfections requires implementing as-manufactured imperfections into high-fidelity finite element method (FEM) models. Geometry measurement methods such as structured light scanning are used for all geometric surface data used in this work. Common preprocessing and visualization steps used in SBKF are discussed, and steps of how surface scans are prepared for implementation into a finite element model is described. The Python Tool for Implementing Geometric Imperfections in Reduced Structures (Py_TIGIRS), written specifically for the use with SBKF, is briefly described and uses eight functions to extract, modify, and write geometric imperfections into Abaqus input files. Results of the preprocessing methods and results from Py_TIGIRS are provided and compared for Composite Test Articles (CTA) 8.2, 8.2B, and 8.3. Excellent agreement between the visualized scan data and the FEM-extracted geometry is demonstrated. A brief example of why geometric surface imperfections are significant in nonlinear numerical analyses for thin cylinders in axial compression is provided as motivation to use tools such as Py_TIGIRS. Future developments of Py_TIGIRS including expansion to structures of arbitrary geometry is planned.

Geometric imperfections↗

Implementing Geometric Surface Imperfections into Sandwich Composite Cylinder Finite Element Method Models

The buckling responses of certain cylindrical shell structures are extremely sensitive to geometric imperfections. The NASA Engineering and Safety Center (NESC) Shell Buckling Knockdown Factor Project (SBKF) is conducting research to develop analysis-based buckling design recommendations. Experiments are used to verify the analysis-based factors, but the sensitivity of the test articles to geometric imperfections requires implementing as-manufactured imperfections into high-fidelity finite element method (FEM) models. Geometry measurement methods such as structured light scanning are used for all geometric surface data used in this work. Common preprocessing and visualization steps used in SBKF are discussed, and steps of how surface scans are prepared for implementation into a finite element model is described. The Python Tool for Implementing Geometric Imperfections in Reduced Structures (Py_TIGIRS), written specifically for the use with SBKF, is briefly described and uses eight functions to extract, modify, and write geometric imperfections into Abaqus input files. Results of the preprocessing methods and results from Py_TIGIRS are provided and compared for Composite Test Articles (CTA) 8.2, 8.2B, and 8.3. Excellent agreement between the visualized scan data and the FEM-extracted geometry is demonstrated. A brief example of why geometric surface imperfections are significant in nonlinear numerical analyses for thin cylinders in axial compression is provided as motivation to use tools such as Py_TIGIRS. Future developments of Py_TIGIRS including expansion to structures of arbitrary geometry is planned.

Geometric imperfections↗

Implementing Geometric Surface Imperfections into Sandwich Composite Cylinder Finite Element Method Models

The buckling responses of certain cylindrical shell structures are extremely sensitive to geometric surface imperfections. The NASA Engineering and Safety Center (NESC) Shell Buckling Knockdown Factor Project (SBKF) is conducting research to develop analysis-based buckling design recommendations. Experiments are used to verify the analysis-based factors, but the sensitivity of the test articles to geometric imperfections requires implementing as-manufactured imperfections into high-fidelity finite element method models. Data collection methods such as structured light scanning are used for all geometric surface data used in this work. Common preprocessing and visualization steps used in SBKF are discussed, and steps on how surface scans are prepared for implementation into a finite element model is described. The Python Tool for Implementing Geometric Imperfections in Reduced Structures (Py_TIGIRS), written specifically for the use with SBKF, is briefly described and uses eight functions to extract, modify, and write geometric imperfections into Abaqus input files. Results of the pre-processing methods and results from Py_TIGIRS are provided and compared for Composite Test Article (CTA) 8.2B. Excellent agreement between the visualized scan data and the FEM-extracted geometry is demonstrated. A brief example of why geometric surface imperfections are significant in nonlinear numerical analyses for thin cylinders in axial compression is provided as motivation to use tools such as Py_TIGIRS. Future development of Py_TIGIRS including expansion to structures of arbitrary geometry is planned.

Sandwich structures↗

Finite element modeling of cracked bodies using the Bodner-Partom flow law

The Bodner-Partom flow law which models viscoplastic material behavior has been used to represent two nickel-base superalloys, Gatorized IN100 and Inconel 718 at elevated temperature. Procedures for the determination of the material parameters are presented along with a discussion of the physical significance of each parameter. The material model is then used in finite element computations to evaluate the response of cracked bodies to monotonic, sustained, or cyclic loading. Geometries investigated include the center cracked panel, the compact tension specimen, and the single cracked ring under tension. A Hybrid Experimental Numerical (HEN) procedure has been used to deduce crack growth rates from experimental displacement measurements which are input into finite element computations. The results of several studies conducted over the last several years are summarized.

Nicholas, T.↗

Performance Portable Graphics Processing Unit Acceleration of a High-Order Finite Element Multiphysics Application

The Lawrence Livermore National Laboratory (LLNL) will soon have in place the El Capitan exascale supercomputer, based on advanced micro devices (AMD) graphics processing units (GPUs). As part of a multiyear effort under the National Nuclear Security Administration (NNSA) Advanced Simulation and Computing (ASC) program, we have been developing marbl, a next generation, performance portable multiphysics application based on high-order finite elements. In previous years, we successfully ported the Arbitrary Lagrangian–Eulerian (ALE), multimaterial, compressible flow capabilities of marbl to nvidia GPUs as described in Vargas et al. Here, in this paper, we describe our ongoing effort in extending marbl's GPU capabilities with additional physics, including multigroup radiation diffusion and thermonuclear burn for high energy density physics (HEDP) and fusion modeling. We also describe how our portability abstraction approach based on the raja Portability Suite and the mfem finite element discretization library has enabled us to achieve high performance on AMD based GPUs with minimal effort in hardware-specific porting. Throughout this work, we highlight numerical and algorithmic developments that were required to achieve GPU performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Compatibility conditions of structural mechanics for finite element analysis

The equilibrium equations and the compatibility conditions are fundamental to the analyses of structures. However, anyone who undertakes even a cursory generic study of the compatibility conditions can discover, with little effort, that historically this facet of structural mechanics had not been adequately researched by the profession. Now the compatibility conditions (CC's) have been researched and are understood to a great extent. For finite element discretizations, the CC's are banded and can be divided into three distinct categories: (1) the interface CC's; (2) the cluster or field CC's; and (3) the external CC's. The generation of CC's requires the separating of a local region, then writing the deformation displacement relation (ddr) for the region, and finally, the eliminating of the displacements from the ddr. The procedure to generate all three types of CC's is presented and illustrated through examples of finite element models. The uniqueness of the CC's thus generated is shown.

Patnaik, Surya N.↗

Compatibility conditions of structural mechanics for finite element analysis

The equilibrium equations and the compatibility conditions are fundamental to the analyses of structures. However, anyone who undertakes even a cursory generic study of the compatibility conditions can discover, with little effort, that historically this facet of structural mechanics had not been adequately researched by the profession. Now the compatibility conditions (CC's) have been researched and are understood to a great extent. For finite element discretizations, the CC's are banded and can be divided into three distinct categories: (1) the interface CC's; (2) the cluster or field CC's; and (3) the external CC's. The generation of CC's requires the separating of a local region, then writing the deformation displacement relation (ddr) for the region, and finally, the eliminating of the displacements from the ddr. The procedure to generate all three types of CC's is presented and illustrated through examples of finite element models. The uniqueness of the CC's thus generated is shown.

Patnaik, S. N.↗

Integrated force method - Compatibility conditions of structural mechanics for finite element analysis

The equilibrium equations and the compatibility conditions are fundamental to the analyses of structures. However, anyone who undertakes even a cursory generic study of the compatibility conditions can discover, with little effort, that historically this facet of structural mechanics had not been adequately researched by the profession. Now the compatibility conditions (CC's) have been researched and are understood to a great extent. For finite element discretizations, the CC's are banded and can be divided into three distinct categories: (1) the interface CC's, (2) the cluster or field CC's, and (3) the external CC's. The generation of CC's requires the separating of a local region, then writing the deformation displacement relation (ddr) for the region, and finally, the eliminating of the displacements from the ddr. The procedure to generate all three types of CC's is presented and illustrated through examples of finite element models. The uniqueness of the CC's thus generated is shown. The utilization of CC's has resulted in the novel integrated force method (IFM). The solution that is obtained by the IFM converges with a significantly fewer number of elements, compared to the stiffness and the hybrid methods.

Patnaik, Surya N.↗

Finite element analysis of the impact of beam heating mode in molten salt corrosion experiments employing simultaneous ion irradiation

Finite element analysis was used to investigate the temperature and stress profiles that develop in 316L stainless steel membranes being irradiated using different proton beam conditions in contact with a molten salt environment. It was shown that in addition to a nonuniform irradiation profile, a focused 2 MeV proton beam leads to very strong temperature and stress gradients in the membrane, introducing highly localized driving forces that complicate and even compromise the integrity and reliability of the experimental results of corrosion studies. Here, the use of a focused beam in corrosion studies can create experimental artifacts that may misrepresent the true corrosion behavior. In contrast, the use of a rastered beam is shown to distribute the protons and resulting radiation damage uniformly across the membrane face, and more importantly, results in temperature and stress profiles that are not only very uniform but are of much lower magnitude. The use of a rastered beam during molten salt corrosion experiments is therefore recommended to achieve uniform damage rates, thereby reducing both gradients and magnitudes of the temperature and stress distributions.

43 - PARTICLE ACCELERATORS↗

Finite Element Modeling of the Posterior Eye in Microgravity

Microgravity experienced during spaceflight affects astronauts in various ways, including weakened muscles and loss of bone density. Recently, visual impairment and intracranial pressure (VIIP) syndrome has become a major concern for space missions lasting longer than 30 days. Astronauts suffering from VIIP syndrome have changes in ocular anatomical and visual impairment that persist after returning to earth. It is hypothesized that a cephalad fluid shift in microgravity may increase the intracranial pressure (ICP), which leads to an altered biomechanical environment of the posterior globe and optic nerve sheath (ONS).Currently, there is a lack of knowledge of how elevated ICP may lead to vision impairment and connective tissue changes in VIIP. Our goal was to develop a finite element model to simulate the acute effects of elevated ICP on the posterior eye and optic nerve sheath. We used a finite element (FE) analysis approach to understand the response of the lamina cribrosa and optic nerve to the elevations in ICP thought to occur in microgravity and to identify which tissue components have the greatest impact on strain experienced by optic nerve head tissues.

intraocular pressure↗

Hybrid-finite-element analysis of some nonlinear and 3-dimensional problems of engineering fracture mechanics

In this paper, efficient numerical methods for the analysis of crack-closure effects on fatigue-crack-growth-rates, in plane stress situations, and for the solution of stress-intensity factors for arbitrary shaped surface flaws in pressure vessels, are presented. For the former problem, an elastic-plastic finite element procedure valid for the case of finite deformation gradients is developed and crack growth is simulated by the translation of near-crack-tip elements with embedded plastic singularities. For the latter problem, an embedded-elastic-singularity hybrid finite element method, which leads to a direct evaluation of K-factors, is employed.

Atluri, S. N.↗

Efficient exascale discretizations: High-order finite element methods

Efficient exploitation of exascale architectures requires rethinking of the numerical algorithms used in many large-scale applications. These architectures favor algorithms that expose ultra fine-grain parallelism and maximize the ratio of floating point operations to energy intensive data movement. One of the few viable approaches to achieve high efficiency in the area of PDE discretizations on unstructured grids is to use matrix-free/partially assembled high-order finite element methods, since these methods can increase the accuracy and/or lower the computational time due to reduced data motion. In this paper we provide an overview of the research and development activities in the Center for Efficient Exascale Discretizations (CEED), a co-design center in the Exascale Computing Project that is focused on the development of next-generation discretization software and algorithms to enable a wide range of finite element applications to run efficiently on future hardware. CEED is a research partnership involving more than 30 computational scientists from two US national labs and five universities, including members of the Nek5000, MFEM, MAGMA and PETSc projects. We discuss the CEED co-design activities based on targeted benchmarks, miniapps and discretization libraries and our work on performance optimizations for large-scale GPU architectures. We also provide a broad overview of research and development activities in areas such as unstructured adaptive mesh refinement algorithms, matrix-free linear solvers, high-order data visualization, and list examples of collaborations with several ECP and external applications.

97 MATHEMATICS AND COMPUTING↗

Finite element structural redesign by large admissible perturbations

In structural redesign, two structural states are involved; the baseline (known) State S1 with unacceptable performance, and the objective (unknown) State S2 with given performance specifications. The difference between the two states in performance and design variables may be as high as 100 percent or more depending on the scale of the structure. A Perturbation Approach to Redesign (PAR) is presented to relate any two structural states S1 and S2 that are modeled by the same finite element model and represented by different values of the design variables. General perturbation equations are derived expressing implicitly the natural frequencies, dynamic modes, static deflections, static stresses, Euler buckling loads, and buckling modes of the objective S2 in terms of its performance specifications, and S1 data and Finite Element Analysis (FEA) results. Large Admissible Perturbation (LEAP) algorithms are implemented in code RESTRUCT to define the objective S2 incrementally without trial and error by postprocessing FEA results of S1 with no additional FEAs. Systematic numerical applications in redesign of a 10 element 48 degree of freedom (dof) beam, a 104 element 192 dof offshore tower, a 64 element 216 dof plate, and a 144 element 896 dof cylindrical shell show the accuracy, efficiency, and potential of PAR to find an objective state that may differ 100 percent from the baseline design.

Bernitsas, Michael M.↗

A Mixed Multi-Field Finite Element Formulation for Thermopiezoelectric Composite Shells

Analytical formulations are presented which account for the coupled mechanical, electrical, and thermal response of piezoelectric composite shell structures. A new mixed multi-field laminate theory is developed which combines "single layer" assumptions for the displacements along with layerwise fields for the electric potential and temperature. This laminate theory is formulated using curvilinear coordinates and is based on the principles of linear thermopiezoelectricity. The mechanics have the inherent capability to explicitly model both the active and sensory responses of piezoelectric composite shells in thermal environment. Finite element equations are derived and implemented for an eight-noded shell element. Numerical studies are conducted to investigate both the sensory and active responses of piezoelectric composite shell structures subjected to thermal loads. Results for a cantilevered plate with an attached piezoelectric layer are com- pared with corresponding results from a commercial finite element code and a previously developed program. Additional studies are conducted on a cylindrical shell with an attached piezoelectric layer to demonstrate capabilities to achieve thermal shape control on curved piezoelectric structures.

Lee, Ho-Jun↗

Design of Rigid-Flex PCB Robotics Leveraging Validated Finite Element Simulations

The use of rigid-flex printed circuit board (PCB) as primary structure has the potential to reduce the weight and volume of robotic systems. In the case of robotics for interplanetary exploration, these systems can leverage origamiinspired folding for increased mobility options and reduced storage volume. Folding rigid-flex PCB robotics can be constructed with rigid PCB connected by short Nomex fabric hinges coupled with flex PCB ribbon cables that permits enhanced system flexibility and energy dissipation to promote impact survivability. This paper presents a design methodology of rigidflex PCB systems with an emphasis on impact resistance. The design process considers solder joint adequacy, panel bending, and fracture using a finite element (FE) model. The proposed design methodology is developed using a case study with NASA JPL’s Pop-Up Folding Flat Explorer Robot (PUFFER). First, the finite-element (FE) modeling methodology is presented with consideration to both frequency and time-domain modeling applications, which include operational self-contact analysis and high impact scenarios. The time-domain impact modeling methodology utilizes hyperelastic material properties for the Nomex hinges. This modeling method is validated using image correlation of PUFFER drop tests. A flowchart is presented to guide users through a validated Abaqus modeling procedure for highly flexible rigid-flex systems. Next, a case study is presented in which PUFFER is subject to drop heights representative of falls into Lunar pits and then the design is refined for a more optimum impact performance. Finally, the results of the case study are used to inform a generalized design methodology for rigid-flex PCB robotics subject to high-impact loads with the considerations presented.

de la Croix, Jean-Pierre↗

Finite Element Modeling of Thermal Cycling Induced Microcracking in Carbon/Epoxy Triaxial Braided Composites

The microcrack distribution and mass change in PR520/T700s and 3502/T700s carbon/epoxy braided composites exposed to thermal cycling was evaluated experimentally. Acoustic emission was utilized to record the crack initiation and propagation under cyclic thermal loading between -55 C and 120 C. Transverse microcrack morphology was investigated using X-ray Computed Tomography. Different performance of two kinds of composites was discovered and analyzed. Based on the observations of microcrack formation, a meso-mechanical finite element model was developed to obtain the resultant mechanical properties. The simulation results exhibited a decrease in strength and stiffness with increasing crack density. Strength and stiffness reduction versus crack densities in different orientations were compared. The changes of global mechanical behavior in both axial and transverse loading conditions were studied. Keywords: Thermal cycles; Microcrack; Finite Element Model; Braided Composite

Zhang, Chao↗

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.

Mital, Subodh K.↗