SEARCH · Engineering Papers

Results for “continuous finite 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Time-domain finite elements in optimal control with application to launch-vehicle guidance

A time-domain finite element method is developed for optimal control problems. The theory derived is general enough to handle a large class of problems including optimal control problems that are continuous in the states and controls, problems with discontinuities in the states and/or system equations, problems with control inequality constraints, problems with state inequality constraints, or problems involving any combination of the above. The theory is developed in such a way that no numerical quadrature is necessary regardless of the degree of nonlinearity in the equations. Also, the same shape functions may be employed for every problem because all strong boundary conditions are transformed into natural or weak boundary conditions. In addition, the resulting nonlinear algebraic equations are very sparse. Use of sparse matrix solvers allows for the rapid and accurate solution of very difficult optimization problems. The formulation is applied to launch-vehicle trajectory optimization problems, and results show that real-time optimal guidance is realizable with this method. Finally, a general problem solving environment is created for solving a large class of optimal control problems. The algorithm uses both FORTRAN and a symbolic computation program to solve problems with a minimum of user interaction. The use of symbolic computation eliminates the need for user-written subroutines which greatly reduces the setup time for solving problems.

Bless, Robert R.↗

NASA NTRS1992

Contact thermal shock test of ceramics

A novel quantitative thermal shock test of ceramics is described. The technique employs contact between a metal-cooling rod and hot disk-shaped specimen. In contrast with traditional techniques, the well-defined thermal boundary condition allows for accurate analyses of heat transfer, stress, and fracture. Uniform equibiaxial tensile stresses are induced in the center of the test specimen. Transient specimen temperature and acoustic emission are monitored continuously during the thermal stress cycle. The technique is demonstrated with soda-lime glass specimens. Experimental results are compared with theoretical predictions based on a finite-element method thermal stress analysis combined with a statistical model of fracture. Material strength parameters are determined using concentric ring flexure tests. Good agreement is found between experimental results and theoretical predictions of failure probability as a function of time and initial specimen temperature.

Rogers, W. P.↗

NASA NTRS2003

Methods Developed by the Tools for Engine Diagnostics Task to Monitor and Predict Rotor Damage in Real Time

Tools for Engine Diagnostics is a major task in the Propulsion System Health Management area of the Single Aircraft Accident Prevention project under NASA s Aviation Safety Program. The major goal of the Aviation Safety Program is to reduce fatal aircraft accidents by 80 percent within 10 years and by 90 percent within 25 years. The goal of the Propulsion System Health Management area is to eliminate propulsion system malfunctions as a primary or contributing factor to the cause of aircraft accidents. The purpose of Tools for Engine Diagnostics, a 2-yr-old task, is to establish and improve tools for engine diagnostics and prognostics that measure the deformation and damage of rotating engine components at the ground level and that perform intermittent or continuous monitoring on the engine wing. In this work, nondestructive-evaluation- (NDE-) based technology is combined with model-dependent disk spin experimental simulation systems, like finite element modeling (FEM) and modal norms, to monitor and predict rotor damage in real time. Fracture mechanics time-dependent fatigue crack growth and damage-mechanics-based life estimation are being developed, and their potential use investigated. In addition, wireless eddy current and advanced acoustics are being developed for on-wing and just-in-time NDE engine inspection to provide deeper access and higher sensitivity to extend on-wing capabilities and improve inspection readiness. In the long run, these methods could establish a base for prognostic sensing while an engine is running, without any overt actions, like inspections. This damage-detection strategy includes experimentally acquired vibration-, eddy-current- and capacitance-based displacement measurements and analytically computed FEM-, modal norms-, and conventional rotordynamics-based models of well-defined damages and critical mass imbalances in rotating disks and rotors.

Baaklini, George Y.↗

DOE OSTI2023

Modular machine learning-based elastoplasticity: Generalization in the context of limited data

The development of highly accurate constitutive models for materials that undergo path-dependent processes continues to be a complex challenge in computational solid mechanics. Challenges arise both in considering the appropriate model assumptions and from the viewpoint of data availability, verification, and validation. Recently, data-driven modeling approaches have been proposed that aim to establish stress-evolution laws that avoid user-chosen functional forms by relying on machine learning representations and algorithms. However, these approaches not only require a significant amount of data but also need data that probes the full stress space with a variety of complex loading paths. Furthermore, they rarely enforce all necessary thermodynamic principles as hard constraints. Hence, they are in particular not suitable for low-data or limited-data regimes, where the first arises from the cost of obtaining the data and the latter from the experimental limitations of obtaining labeled data, which is commonly the case in engineering applications. In this work, we discuss a hybrid framework that can work on a variable amount of data by relying on the modularity of the elastoplasticity formulation where each component of the model can be chosen to be either a classical phenomenological or a data-driven model depending on the amount of available information and the complexity of the response. The method is tested on synthetic uniaxial data coming from simulations as well as cyclic experimental data for structural materials. The discovered material models are found to not only interpolate well but also allow for accurate extrapolation in a thermodynamically consistent manner far outside the domain of the training data. This ability to extrapolate from limited data was the main reason for the early and continued success of phenomenological models and the main shortcoming in machine learning-enabled constitutive modeling approaches. Training aspects and details of the implementation of these models into Finite Element simulations are discussed and analyzed.

42 ENGINEERING↗

NASA NTRS1978

Recent advances in thermal-structural analysis and design

A continuing effort to develop a comprehensive capability for thermal-structural analysis and automated design (sizing) is reviewed. A principal role in the activity is played by the finite element program SPAR which contains both an efficient structural and thermal analysis capability. The benefit of having thermal and structural analyses in the same finite element program is illustrated by the application of SPAR to design calculations for the National Transonic Facility - a cryogenic wind tunnel under construction at Langley. Some experience with large-scale thermal structural analysis problems, particularly the space shuttle orbiter, has led to the identification of some analysis needs. Those needs include automated model generation and data output for lumped parameter thermal analysis, faster solution methods for nonlinear transient heat transfer, automated interpolation of temperature data from a thermal finite element model to a dissimilar structural finite model, and automated techniques to identify the times at which the critical combinations of transient heating and loads occur on a structure.

Adelman, H. M.↗

NASA NTRS1996

A Thermostructural Analysis of a Diboride Composite Leading Edge

In an effort to support the design of zirconium diboride composite leading edges for hypersonic vehicles, a finite element model (FEM) of a prototype leading edge was created and finite element analysis (FEA) was employed to assess its thermal and structural response to aerothermal boundary conditions. Unidirectional material properties for the structural components of the leading edge, a continuous fiber reinforced diboride composite, were computed with COSTAR. These properties agree well with those experimentally measured. To verify the analytical approach taken with COSMOS/M, an independent FEA of one of the leading edge assembly components was also done with COSTAR. Good agreement was obtained between the two codes. Both showed that a unidirectional lay-up had the best margin of safety for a simple loading case. Both located the maximum stress in the same region and ply. The magnitudes agreed within 4 percent. Trajectory based aerothermal heating was then applied to the leading edge assembly FEM created with COSMOS/M to determine steady state temperature response, displacement, stresses, and contact forces due to thermal expansion and thermal strains. Results show that the leading edge stagnation line temperature reached 4700 F. The maximum computed failure index for the laminated composite components peaks at 4.2, and is located at the bolt flange in layer 2 of the side bracket. The temperature gradient in the tip causes a compressive stress of 279 ksi along its width and substantial tensile stresses within its depth.

Kowalski, Tom↗

NASA NTRS1998

Failure of Laminated Composites at Thickness Discontinuities: An Experimental and Analytical Study

Failure initiation of laminated composites at a thickness discontinuity is studied experimentally with the aid of an optical microscope under combined loading of tension, transverse shear and bending, making use of three- and four-point bending arrangements. Because transverse shear produced relatively small effects in failure initiation results are presented as tension-bending interactions. Two loading frames for three- and four-point bending were designed to apply moment and tension simultaneously to produce failure by generating a ply crack; this initiation was evaluated by finite element analysis using ABAQUS. For cross-plies bounding the interface at the base of the step it is found that a maximum strain criterion applied to the continuing (long) ply describes the failure initiation. Ultimate failure resulted at loads on the order of 25 to 35% higher than those at failure initiation.

Lee, Sangwook↗

NASA NTRS2024

Design of an In-Slot Cooled Air-Core Flux-Focusing Permanent Magnet Synchronous Machine for Electric Aircraft Applications

High performance and reliable permanent magnet synchronous machines are potentially a key technology for enabling reduced emissions in future generations of sustainable aircraft. In this paper, a novel air-core flux focusing permanent magnet synchronous machine with carbon fiber shafting for aviation applications is presented. A low-fidelity analytical sizing tool is detailed for the motor topology. The sizing tool is used to explore design trades and create an example 100 kW flux-focusing motor for eVTOL applications. Higher fidelity finite element analysis was used to validate the predicted performance of the motor. The selected design is able to achieve roughly 97 percent efficiency and 6 kW/kg at a continuous operating power of 100 kW.

Jonathan Gutknecht↗

NASA NTRS1983

Research on the control of large space structures

The research effort on the control of large space structures at the University of Houston has concentrated on the mathematical theory of finite-element models; identification of the mass, damping, and stiffness matrix; assignment of damping to structures; and decoupling of structure dynamics. The objective of the work has been and will continue to be the development of efficient numerical algorithms for analysis, control, and identification of large space structures. The major consideration in the development of the algorithms has been the large number of equations that must be handled by the algorithm as well as sensitivity of the algorithms to numerical errors.

Denman, E. D.↗

NASA NTRS1987

Hierarchically partitioned nonlinear equation solvers

By partitioning solution space into a number of subspaces, a new multiply constrained partitioned Newton-Raphson nonlinear equation solver is developed. Specifically, for a given iteration, each of the various separate partitions are individually and simultaneously controlled. Due to the generality of the scheme, a hierarchy of partition levels can be employed. For finite-element-type applications, this includes the possibility of degree-of-freedom, nodal, elemental, geometric substructural, material and kinematically nonlinear group controls. It is noted that such partitioning can be continuously updated, depending on solution conditioning. In this context, convergence is ascertained at the individual partition level.

Padovan, Joseph↗

NASA NTRS1993

TransFit: Finite element analysis data fitting software

The Advanced X-Ray Astrophysics Facility (AXAF) mission support team has made extensive use of geometric ray tracing to analyze the performance of AXAF developmental and flight optics. One important aspect of this performance modeling is the incorporation of finite element analysis (FEA) data into the surface deformations of the optical elements. TransFit is software designed for the fitting of FEA data of Wolter I optical surface distortions with a continuous surface description which can then be used by SAO's analytic ray tracing software, currently OSAC (Optical Surface Analysis Code). The improved capabilities of Transfit over previous methods include bicubic spline fitting of FEA data to accommodate higher spatial frequency distortions, fitted data visualization for assessing the quality of fit, the ability to accommodate input data from three FEA codes plus other standard formats, and options for alignment of the model coordinate system with the ray trace coordinate system. TransFit uses the AnswerGarden graphical user interface (GUI) to edit input parameters and then access routines written in PV-WAVE, C, and FORTRAN to allow the user to interactively create, evaluate, and modify the fit. The topics covered include an introduction to TransFit: requirements, designs philosophy, and implementation; design specifics: modules, parameters, fitting algorithms, and data displays; a procedural example; verification of performance; future work; and appendices on online help and ray trace results of the verification section.

Freeman, Mark↗

DOE OSTI2025

Adaptation of Wallace's approach to the specific heat of elemental solids with significant intrinsic anharmonicity, particularly the light actinide metals

The quasiharmonic approximation is the most common method for modeling the specific heat of solids; however, it fails to capture the effects of intrinsic anharmonicity. In this study, we introduce the “elastic softening approximation,” an alternative approach to modeling intrinsic anharmonic effects on thermodynamic quantities, which is grounded in Wallace's thermodynamic framework that tracks entropy changes resulting from the continuous change (e.g., softening) of phonons as a function of temperature. A key finding of our study is a correlation between Poisson's ratio and the differential rate of phonon softening at finite frequencies, compared to lower frequencies relevant to elastic moduli measurements. We observe that elemental solids such as α -Be, diamond, Al, Cu, In, W, Au, and Pb, which span a wide range of Poisson's ratios and exhibit varying degrees of intrinsic anharmonicity, consistently follow this trend. When applied to α -U, α -Pu, and δ -Pu, our method reveals unusually large anharmonic phonon contributions at elevated temperatures across all three light actinide metals. These findings are attributed to the unique combination of enhanced covalency and softer elastic moduli inherent in the actinides, potentially influenced by their 5 f -electron bonding. Published by the American Physical Society 2025

36 MATERIALS SCIENCE↗

DOE OSTI2020

Comparison of Thermal Management Approaches for Integrated Traction Drives in Electric Vehicles

The continuous push to increase power densities of electric vehicle (EV) traction drive systems necessitates combining electric motor and power electronics into one unit. A single, compact traction drive unit with fewer interconnecting components also facilitates fast, automated assembly of electric vehicles, driving production costs down and enabling wider adoption of EVs. There are a number of challenges associated with the integration of power electronics with the electric machine, including thermal management of the combined traction drive system. However, one important benefit of integration from the thermal management system perspective is the potential for using a single fluid loop instead of two separate cooling systems for the electric machine and the power electronics/inverter. This paper reviews several integration approaches and, employing finite element analysis (FEA), compares thermal management solutions for the combined electric machine and power electronics systems. Namely, three different scenarios are modeled: (1) independent component (motor and power electronics) cooling, which is compared to the combined cooling system approach for (2) radially and (3) axially integrated power electronics modules into the motor enclosure. Temperature distributions for selected thermal loads and thermal resistances from the key heat-generating components to the cooling fluid are compared for each scenario.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

A least-squares finite element method based on the Helmholtz decomposition for hyperbolic balance laws

In this paper, a least-squares finite element method for scalar nonlinear hyperbolic balance laws is proposed and studied. The approach is based on a formulation that utilizes an appropriate Helmholtz decomposition of the flux vector and is related to the standard notion of a weak solution. This relationship, together with a corresponding connection to negative-norm least-squares, is described in detail. As a consequence, an important numerical conservation theorem is obtained, similar to the famous Lax–Wendroff theorem. The numerical conservation properties of the method in this paper do not fall precisely in the framework introduced by Lax and Wendroff, but they are similar in spirit as they guarantee that when L 2 convergence holds, the resulting approximations approach a weak solution to the hyperbolic problem. The least-squares functional is continuous and coercive in an H -1 -type norm, but not L 2 -coercive. Nevertheless, the L 2 convergence properties of the method are discussed. Convergence can be obtained either by an explicit regularization of the functional, that provides control of the L 2 norm, or by properly choosing the finite element spaces, providing implicit control of the L 2 norm. Numerical results for the inviscid Burgers equation with discontinuous source terms are shown, demonstrating the L 2 convergence of the obtained approximations to the physically admissible solution. The numerical method utilizes a least-squares functional, minimized on finite element spaces, and a Gauss–Newton technique with nested iteration. Finally, we believe that the linear systems encountered with this formulation are amenable to multigrid techniques and combining the method with adaptive mesh refinement would make this approach an efficient tool for solving balance laws (this is the focus of a future study).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Digital image correlation and infrared thermography data for seven unique geometries of 304L stainless steel

Material Testing 2.0 (MT2.0) is a paradigm that advocates for the use of rich, full-field data, such as from digital image correlation and infrared thermography, for material identification. By employing heterogeneous, multi-axial data in conjunction with sophisticated inverse calibration techniques such as finite element model updating and the virtual fields method, MT2.0 aims to reduce the number of specimens needed for material identification and to increase confidence in the calibration results. To support continued development, improvement, and validation of such inverse methods—specifically for rate-dependent, temperature-dependent, and anisotropic metal plasticity models—we provide here a thorough experimental data set for 304L stainless steel sheet metal. The data set includes full-field displacement, strain, and temperature data for seven unique specimen geometries tested at different strain rates and in different material orientations. Commensurate extensometer strain data from tensile dog bones is provided as well for comparison. We believe this complete data set will be a valuable contribution to the experimental and computational mechanics communities, supporting continued advances in material identification methods.

36 MATERIALS SCIENCE↗

NASA NTRS1989

Axisymmetric shapes and stability of charged drops in an external electric field

A highly conducting charged drop that is surrounded by a fluid insulator of another density can be levitated by suitably applying a uniform electric field. Axisymmetric equilibrium shapes and stability of the levitated drop are found by solving simultaneously the augmented Young-Laplace equation for surface shape and the Laplace equation for the elecric field, together with constraints of fixed drop volume, charge, and center of mass. The means are a method of subdomains, finite element basis functions, and Galerkin's method of weighted residuals, all facilitated by a large-scale computer. Shape families of fixed charge are treated systematically by first-order continuation. Previous analyses by Abbas et al. in 1967 and Abbas and Latham in 1969, in which the shapes of levitated drops are approximated as spheroids, are corrected. The new analysis shows that drops charged to less than the Rayleigh limit lose shape stability at turning points, with respect to external field strength, and that the instability seen in experiments of Doyle et al. in 1964 and others is not a bifurcation to a family of two-lobed shapes, but rather is a related imperfect bifurcation.

Basaran, O. A.↗

NASA NTRS1990

Dynamic characteristics of a Space Station Freedom Mars evolution reference configuration

One concept for a manned mission to Mars uses as a transportation mode an evolutionary version of Space Station Freedom (SSF), termed the Mars Evolution Reference Configuration (MERC). The MERC is configured by adding to SSF dual keels, an upper and lower boom, additional laboratory and habitation modules, increased power, and an assembly platform. A finite-element model of the MERC was formulated to investigate the expected low frequency modes, and its variations with the addition of large payload. A basic reboost procedure using near-continuous firing of reaction control system jets is described with the closed-loop attitude control system. The elastic dynamic behavior at critical points during a reboost of the MERC, both with and without the Mars piloted vehicle installed, is examined.

Ayers, J. Kirk↗

NASA NTRS1990

Dynamic Characteristics of Space Station Freedom Mars and Lunar Evolution Reference Configurations

One concept for a manned mission to Mars uses an evolutionary version of Space Station Freedom (SSF) as a transportation node. The station is modified by the addition of dual keels, an upper and lower boom, additional laboratory and habitation modules, increased power and an assembly platform. With these modifications the station is called the Mars Evolution Reference Configuration (MERC). The mass of the station is 65 percent greater than the mass of SSF and its moments of inertia through the mass center are greater by approximately a factor of four. Over a period of months, several flights from Earth to low-Earth-orbit carry the components of a manned Mars piloted vehicle (MPV) to the MERC where the vehicle is constructed on the assembly platform. After each flight the station is reboosted to an appropriate altitude, such that the orbit decay due to atmospheric drag forces lowers the spacecraft to the proper altitude at the appropriate time for rendezvous with the next assembly flight. When the assembly process is completed, the MPV, which has a mass of approximately 200,000 lbm, is situated on the evolutionary station. The mass increase of the MERC with MPV system over SSF is 112 percent and the moments of inertia about axes through the mass center increase by up to a factor of 12. When the MPV is assembled, inspected and verified, the mission is ready to proceed and the MPV is moved from the station to a staging area and mated with fueled trans-Mars injection stages for the flight to Mars. This presentation describes a finite element model of the MERC formulated to investigate the expected low frequency modes and its variation with the addition of a large payload. A basic reboost procedure using near-continuous firing of reaction control system jets is proposed with off-modulation of the jets used to control flight attitude. The reboost procedure is described with the closed-loop attitude control dictating jet on/off cycling based on feedback signals which contain both the rigid body rotation information and the elastic rotations local to the attitude sensor. The presentation contains a description of the dynamic response at critical points of the station during the reboost and concludes with results of a brief study of the dynamic characteristics of a Lunar transportation node configuration.

Ayers, J. Kirk↗