Engineering Papers⌕ Search

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.

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

At least 199 records · Page 11

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↗

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↗

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.↗

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↗

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↗

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.↗

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↗

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↗

Aeroelastic Trim Drag Optimization of Mach 0.745 Transonic Truss-Braced Wing Aircraft with Variable-Camber Continuous Trailing-Edge Flap

This paper presents an aeroelastic trim drag optimization study of the Mach 0.745 Transonic Truss-Braced Wing (TTBW) aircraft with the Variable-Camber Continuous Trailing-Edge Flap (VCCTEF). An aero-structural analysis solver VSPAERO with transonic small disturbance, integral boundary-layer, and wing-strut interference corrections coupled to finite-element code BEAM3D is developed to provide a rapid aircraft aeroelastic performance evaluation. The VSPAERO model has been validated against wind tunnel test data. Aeroelastic trim drag optimization studies are conducted for a VCCTEF configuration with 6-spanwise sections. Three different flight conditions corresponding to Mach 0.745 are selected for the TTBW aircraft aeroelastic trim drag optimization at the design and off-design cruise lift coefficients. The preliminary optimization results show that the TTBW aircraft with the optimized VCCTEF deflection achieves a drag reduction of about 9.1 counts, 5.9 counts, and 11.8 counts corresponding to the lift coefficients 0.68, 0.73, and 0.78, respectively. A high-fidelity CFD solver FUN3D is used to validate the aeroelastic trim drag optimization.

Aeroelastic trim drag↗

Testing of a Composite Conical-Cylindrical Shell

Launch-vehicle shell structures, which can be comprised of both cylindrical and conical sections, are known to be susceptible to buckling due to their large radius-to-thickness ratios. The advancements in composite manufacturing and numerical methods have enabled designers to consider more nontraditional shapes, such as connecting the conical and cylindrical sections with a toroidal transition to create a single-piece conical-cylindrical shell. This single-piece construction eliminates the need for a heavy interface ring between sections and has the potential to save mass. To better understand the buckling behavior, a composite conical-cylindrical shell was designed, fabricated, and tested. Prior to test, a finite element model that included thickness variations and radial imperfections was created. The test article buckled elastically at 251.8 kN, approximately 8.8% higher than the predicted buckling load of 231.4 kN Continued research in conical-cylindrical structures has the potential to expand the design space for launch-vehicle structures and lead to improved designs and reduced weight.

Buckling↗

Test and Analysis of a Composite Conical-Cylindrical Shell

Launch-vehicle shell structures, which can be comprised of both cylindrical and conical sections, are known to be susceptible to buckling due to their large radius-to-thickness ratios. The advancements in composite manufacturing and numerical methods have enabled designers to consider more nontraditional shapes, such as connecting the conical and cylindrical sections with a toroidal transition to create a single-piece conical-cylindrical shell. This single-piece construction eliminates the need for a heavy interface ring between sections and has the potential to save mass. To better understand the buckling behavior, a composite conical-cylindrical shell was designed, fabricated, and tested. Prior to test, a finite element model that included thickness variations and radial imperfections was created. The test article buckled elastically at 251.8 kN, approximately 8.8% higher than the predicted buckling load of 231.4 kN Continued research in conical-cylindrical structures has the potential to expand the design space for launch-vehicle structures and lead to improved designs and reduced weight.

Buckling↗

Effect of Phenolic Matrix Microcracking on the Structural Response of a 3-D Woven Thermal Protection System

The effect of microcracking in the phenolic matrix of a three-dimensional woven thermal protection system (TPS) and the resulting material stiffness reduction was studied via a comparison of finite element analysis results from a linear analysis and an iterative linear analysis. A TPS is necessary to protect space vehicles from the aerodynamic heating of planetary entry. The Heatshield for Extreme Entry Environment Technology (HEEET) project has developed a TPS for use in high heat-flux and pressure missions. The material is a dual-layer continuous dry weave, which is then infiltrated with a low-density phenolic resin matrix to form a rigid ablator. The phenolic resin matrix does not have structural load transfer requirements, and testing has shown that the phenolic resin can fully satisfy thermal requirements when the matrix contains microcracks. Due to high stresses in the through-the-thickness direction of the material, phenolic microcracks may form in the matrix material, which would result in a reduction of stiffness. An exploratory study was conducted to determine if reduction in material stiffness would change the load paths and/or decrease the structural margins. A comparison was performed between a linear finite element analysis that did not take into account phenolic microcracking and an iterative linear finite element analysis that accounted for propagation of phenolic microcracking. Four subcases using varying assumptions were analyzed and the results indicate that the assumed strength at which the phenolic microcracking propagates was the critical parameter for determining the extent of microcracking in the phenolic matrix. Phenolic microcracking does not have an adverse effect on the structural response of the test article and is not a critical failure.

Langston, Sarah L.↗

Directions for computational mechanics in automotive crashworthiness

The automotive industry has used computational methods for crashworthiness since the early 1970's. These methods have ranged from simple lumped parameter models to full finite element models. The emergence of the full finite element models in the mid 1980's has significantly altered the research direction. However, there remains a need for both simple, rapid modeling methods and complex detailed methods. Some directions for continuing research are discussed.

Bennett, James A.↗

High-Fidelity Modeling Approach for Coupling Thermal, Fluid, and Neutronic Analysis in an NTP Engine

Nuclear thermal propulsion (NTP), with its increased efficiency compared to traditional rocket engines, is a potentially enabling technology for space missions of higher complexity, such as delivering heavy machinery or crew to the moon or other planets. Thus, exploring modeling approaches that can effectively couple the thermal, fluid, and neutronic behaviors of NTP reactors is essential to mitigate cost and schedule challenges with hardware development. This paper details one such approach which utilizes Cardinal to unify reactor physics analysis performed using OpenMC, heat transfer analysis performed using the MOOSE framework, and thermal hydraulic analysis performed using the MOOSE Thermal Hydraulic Module. This approach is demonstrated with a proof of concept that focuses on a legacy reactor geometry from the historic NERVA program. This preliminary model found an overprediction in, when compared to the reference data, axial power peaking and the temperatures of both solids and fluids by roughly 20%. It is clear, however, that the model conserved a number of broad trends, and a few model refinements should reduce the observed error. This approach has definite merit and should enable the production of informative results for more modern NTP systems. This work is a first step in creating a transient reactor model that can be updated across a sequence of time steps based on continually fluctuating, time-dependent boundary conditions.

Nuclear Thermal Propulsion↗

Curie isotherm surfaces inferred from high-altitude magnetic anomaly data

Two-dimensional Curie depth models along two sections across the western United States are derived using an equivalent layer magnetization model derived from Magsat data, and the results are used to constrain finite element thermal models developed by Mayhew and Majer (1980). Regional heat flow variations predicted by the models compare favorably with those inferred from direct measurements. The methodology is applied to upward-continued aeromagnetic data, and the results are found to be in good agreement with previous Curie depth estimates and average measured heat flow.

Mayhew, M. A.↗

Mechanical characterization and modeling of non-linear deformation and fracture of a fiber reinforced metal matrix composite

The nonlinear anisotropic mechanical behavior of an aluminum alloy metal matrix composite reinforced with continuous alumina fibers was determined experimentally. The mechanical behavior of the composite were modeled by assuming that the composite has a periodical microstructure. The resulting unit cell problem was solved with the finite element method. Excellent agreement was found between theoretically predicted and measured stress-strain responses for various tensile and shear loadings. The stress-strain responses for transverse and inplane shear were found to be identical and this will provide a simplification of the constitutive equations for the composite. The composite has a very low ductility in transverse tension and a limited ductility in transverse shear that was correlated to high hydrostatic stresses that develop in the matrix. The shape of the initial yield surface was calculated and good agreement was found between the calculated shape and the experimentally determined shape.

Jansson, S.↗

A hybrid formulation of a component mode synthesis method

Component mode synthesis is a substructuring technique frequently employed in structural dynamics. In this method, a given structure is subdivided into components or substructures, each of which is analyzed independently for natural frequencies and for mode shapes. The substructure mode shapes are then assembled to give displacement shapes or load patterns of the original structure. An analytical justification of the basic concept is presented using spectral decompositions, and a variant substructuring approach where intersubstructure continuity is enforced in a weak form is derived. This leads to a hybrid formulation of the basic method which is particularly suitable for assembling heterogeneous substructures and analyzing nonconforming and incompatible finite element substructure models. For problems where both the basic and hybrid methods are applicable, the hybrid variant can be computationally more advantageous.

Farhat, Charbel↗

Shuttle Debris Impact Tool Assessment Using the Modern Design of Experiments

Computational tools have been developed to estimate thermal and mechanical reentry loads experienced by the Space Shuttle Orbiter as the result of cavities in the Thermal Protection System (TPS). Such cavities can be caused by impact from ice or insulating foam debris shed from the External Tank (ET) on liftoff. The reentry loads depend on cavity geometry and certain Shuttle state variables, among other factors. Certain simplifying assumptions have been made in the tool development about the cavity geometry variables. For example, the cavities are all modeled as shoeboxes , with rectangular cross-sections and planar walls. So an actual cavity is typically approximated with an idealized cavity described in terms of its length, width, and depth, as well as its entry angle, exit angle, and side angles (assumed to be the same for both sides). As part of a comprehensive assessment of the uncertainty in reentry loads estimated by the debris impact assessment tools, an effort has been initiated to quantify the component of the uncertainty that is due to imperfect geometry specifications for the debris impact cavities. The approach is to compute predicted loads for a set of geometry factor combinations sufficient to develop polynomial approximations to the complex, nonparametric underlying computational models. Such polynomial models are continuous and feature estimable, continuous derivatives, conditions that facilitate the propagation of independent variable errors. As an additional benefit, once the polynomial models have been developed, they require fewer computational resources to execute than the underlying finite element and computational fluid dynamics codes, and can generate reentry loads estimates in significantly less time. This provides a practical screening capability, in which a large number of debris impact cavities can be quickly classified either as harmless, or subject to additional analysis with the more comprehensive underlying computational tools. The polynomial models also provide useful insights into the sensitivity of reentry loads to various cavity geometry variables, and reveal complex interactions among those variables that indicate how the sensitivity of one variable depends on the level of one or more other variables. For example, the effect of cavity length on certain reentry loads depends on the depth of the cavity. Such interactions are clearly displayed in the polynomial response models.

DeLoach, R.↗