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 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS2013

Structural Assessment of Advanced Composite Tow-Steered Shells

The structural performance of two advanced composite tow-steered shells, manufactured using a fiber placement system, is assessed using both experimental and analytical methods. The fiber orientation angles vary continuously around the shell circumference from 10 degrees on the shell crown and keel, to 45 degrees on the shell sides. The two shells differ in that one shell has the full 24-tow course applied during each pass of the fiber placement system, while the second shell uses the fiber placement system s tow drop/add capability to achieve a more uniform shell wall thickness. The shells are tested in axial compression, and estimates of their prebuckling axial stiffnesses and bifurcation buckling loads are predicted using linear finite element analyses. These preliminary predictions compare well with the test results, with an average agreement of approximately 10 percent.

Wu, K. Chauncey↗

NASA NTRS1983

Improved transistor-controlled and commutated brushless DC motors for electric vehicle propulsion

The development, design, construction, and testing processes of two electronically (transistor) controlled and commutated permanent magnet brushless dc machine systems, for propulsion of electric vehicles are detailed. One machine system was designed and constructed using samarium cobalt for permanent magnets, which supply the rotor (field) excitation. Meanwhile, the other machine system was designed and constructed with strontium ferrite permanent magnets as the source of rotor (field) excitation. These machine systems were designed for continuous rated power output of 15 hp (11.2 kw), and a peak one minute rated power output of 35 hp (26.1 kw). Both power ratings are for a rated voltage of 115 volts dc, assuming a voltage drop in the source (battery) of about 5 volts. That is, an internal source voltage of 120 volts dc. Machine-power conditioner system computer-aided simulations were used extensively in the design process. These simulations relied heavily on the magnetic field analysis in these machines using the method of finite elements, as well as methods of modeling of the machine power conditioner system dynamic interaction. These simulation processes are detailed. Testing revealed that typical machine system efficiencies at 15 hp (11.2 kw) were about 88% and 84% for the samarium cobalt and strontium ferrite based machine systems, respectively. Both systems met the peak one minute rating of 35 hp.

Demerdash, N. A.↗

NASA NTRS1985

Float Zone Experiments in Space

The objective of this work has been to evaluate whether or not Marangoni flow could be suppressed in molten metals by the presence of very thin oxide films. Experimental work has been carried out on molten Sn under UHV conditions. A disk floating zone arrangement was developed to allow in situ Auger examination of molten surfaces. An electron energy loss technique was developed which allows detection of continuous tin oxide films of 6 A or greater. Experiments were planned to detect the effects of oxide formation upon Marangoni flow by measuring: (1) temperature profiles, (2) solid liquid interface shapes, (3) macrosegregation, and (4) the onset of oscillatory Marangoni flow by detecting oscillating temperature variations. Work on (4) showed that oscillatory temperature variations of frequency or = 10 Hz were not present in the disk float zone geometry under conditions of Ma = 4300 with an oxide free molten surface. The disk float zone geometry was modeled with a finite element analysis and temperature and velocity profiles were determined.

Verhoeven, J. D.↗

NASA NTRS1992

Adaptive explicit and implicit finite element methods for transient thermal analysis

The application of adaptive finite element methods to the solution of transient heat conduction problems in two dimensions is investigated. The computational domain is represented by an unstructured assembly of linear triangular elements and the mesh adaptation is achieved by local regeneration of the grid, using an error estimation procedure coupled to an automatic triangular mesh generator. Two alternative solution procedures are considered. In the first procedure, the solution is advanced by explicit timestepping, with domain decomposition being used to improve the computational efficiency of the method. In the second procedure, an algorithm for constructing continuous lines which pass only once through each node of the mesh is employed. The lines are used as the basis of a fully implicit method, in which the equation system is solved by line relaxation using a block tridiagonal equation solver. The numerical performance of the two procedures is compared for the analysis of a problem involving a moving heat source applied to a convectively cooled cylindrical leading edge.

Probert, E. J.↗

NASA NTRS1980

Calculation of residual principal stresses in CVD boron on carbon filaments

A three-dimensional finite element model of the chemical vapor deposition of boron on a carbon substrate (B/C) is developed. The model includes an expansion of the boron after deposition due to atomic rearrangement and includes creep of the boron and carbon. Curves are presented showing the variation of the principal residual stresses and the filament elongation with the parameters defining deposition strain and creep. The calculated results are compared with experimental axial residual stress and elongation measurements made on B/C filaments. For good agreement between calculated and experimental results, the deposited boron must continue to expand after deposition, and the build up of residual stresses must be limited by significant boron and carbon creep rates.

Behrendt, D. R.↗

NASA NTRS1989

Application and implementation of transient algorithms in computer programs

A brief introduction is given to the nonlinear finite element programs developed at Lawrence Livermore National Laboratory. The four programs are DYNA3D and DYNA2D, which are explicit hydrocodes, and NIKE3D and NIKE2D, which are implicit programs. The main emphasis is on DYNA3D with asides about the other programs. During the past year several new features were added to DYNA3D, and major improvements were made in the computational efficiency of the shell and beam elements. Most of these new features and improvements will eventually make their way into the other programs. The emphasis in the computational mechanics effort was always, and continues to be, efficiency. To get the most out of the supercomputers, all Crays, the programs were vectorized where possible. Several of the more interesting capabilities of DYNA3D will be described and the impact on efficiency will be discussed. Some of the recent work on NIKE3D and NIKE2D will also be presented. In the belief that a single example is worth a thousand equations, the theory is skipped entirely and the examples presented.

Benson, David J.↗

NASA NTRSDate not supplied

Using Distributed Fiber-optic Strain Sensing to Estimate Generalized Modal Coordinates from Flight-test Data

A method for estimating the generalized modal coordinates of an aircraft during flight has been developed. The Fiber-optic Sensing System (FOSS) offers an efficient and cost-effective method of measuring the strain at thousands of points along the wings. The estimation of modal coordinates was implemented as a two-step process. First, a maximum likelihood method is used to estimate the statistical properties of the sensors and generalized modal coordinates. Second, the strain mode shapes from the finite element model are used along with the statistical properties from the first step to estimate the generalized modal coordinates over time. Using simulated data from the X-56A Multi-Utility Technology Testbed (MUTT), different methods of modal coordinate estimation were compared to demonstrate the benefits and weaknesses of each. These were compared against the exact solution and an ordinary least squares (a more traditional method) estimate. Modal coordinate estimation methods were then applied to flight-test data from the X-56A aircraft to show that the method continues to work well with actual test data. The new estimation method provides insights unavailable from more classical approaches.

Jeffrey Ouellette↗

NASA NTRSDate not supplied

A Moving Discontinuous Galerkin Method with Interface Condition Enforcement for Reacting Hypersonic Flows

The necessity of enforcing conservation in computational elements or cells (element conservation) for discontinuous solutions is well understood and respected for solving conservation laws in computational fluid dynamics (CFD). In contrast, interface conservation, where the conservation across cell interfaces is enforced, is long ignored, and yet is also ruled and required by the underlying physics just like element conservation. Violation of the interface conservation across discontinuities is the root cause why an exact discontinuous solution can never be achieved in shock capturing methods. The interface conservation is examined and explored in this talk. Moving discontinuous Galerkin (MDG) finite element method with interface conservation enforcement (MDG-ICE) [1],[2] are then presented for solving compressible flow problems with discontinuities based on the observation that the interface conservation can only be satisfied, only when mesh interfaces are aligned with discontinuities. In the MDG-ICE formulation, both conservative quantities and grid geometry are considered as independent variables. A space-time DG formulation is used to solve the multi-material compressible Euler equations in the standard discontinuous solution space and the discrete grid geometry is solved using a variational formulation in a continuous space. A self-adaptive Levenberg-Marquardt method is utilized to solve the resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation. A number of numerical experiments for a variety of flow problems are conducted to assess the accuracy and performance of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to deliver the designed order of both h- and p-convergence even for discontinuous solutions, and detect all types of interfaces, via interface condition enforcement and satisfy, via grid movement, the compressible Euler equations and the associated interface condition.

Discontinuous Galekin↗

NASA NTRS2015

Static Aeroelastic Modeling of a Sub-Scale Wind Tunnel Model with Novel Flap Concept

This paper examines the static aeroelastic modeling of a flexible wind tunnel model equipped with a novel control effector known as the Variable Camber Continuous Trailing Edge Flap (VCCTEF) system. The wind tunnel model is an approximately 10 percent sub-scale version of a wing designed to be equipped on a commerical transport aircraft such as the full-scale NASA Generic Transport Model (GTM). The structure of the model is made highly flexible such that a 10 percent of semi-span wing tip deflection is expected at a design lift condition, and a representation of the VCCTEF concept is incorporated on the model. Static aeroelastic modeling is conducted by using a representative single beam structural finite-element model coupled to a vortex-lattice aerodynamic model. The resulting aeroelastic model of the flexible wind tunnel model is compared against experimental wind tunnel test results from the actual model tested at the University of Washington Aeronautical Laboratory (UWAL). Comparison cases are made using reference and VICON measurements of VCCTEF flap deflections from the wind tunnel test data, and percent errors between the lift curve parameters are presented. The results show up to a maximum percent of error of <10% with regards to lift curve slope values between the aeroelastic model and the UWAL test data. The agreement of lift curve slope values from the aeroelastic model and experimental test results serves as validation for the coupled vortex-lattice finite-element static aeroelastic model.

Aeroelasticity↗

NASA NTRS2018

A Higher-Order Space-Time Finite-Element Method for Moving-Body and Fluid-Structure Interaction Problems

We present a high-order finite-element method for moving body and fluid/structure interaction problems. Our solution strategy is based on a space-time discontinuous Galerkin (DG) spectral-element discretization which extends to arbitrary order of accuracy. The space-time DG discretization is a natural choice for moving body and fluid-structure interaction problems as moving surfaces are incorporated simply by considering curved space-time elements whose space-time faces align with the moving body. We present a discontinuous-Galerkin in time discretization for six-degree of motion modeling of rigid bodies, and a continuous-Galerkin discretization for equations of linear elasticity to generate curved space-time meshes. Numerical results for several simple 2D test cases are presented in order to verify the implementation of the different models. Finally we present a preliminary dynamic simulation of a parachute.

High-Order↗

NASA NTRS1980

The finite element duct eigenvalue problem - An improved formulation with Hermitian elements and no-flow condensation

Hermitian elements are used in a finite element solution for the eigenvalue problem in lined ducts with flow. These elements give significantly greater accuracy for reduced dimensionality when compared with Lagrangian elements. Spurious mode generation associated with the Lagrangian formulation is eliminated. A dramatic improvement in the ratio of the number of reliable eigenvalues to the total number of computed eigenvalues is effected by the use of a condensation scheme based on the no-flow eigenvectors. Results are presented for two dimensional and axisymmetric ducts. In the axisymmetric case good resolution is obtained even for high order, high frequency modes by the use of continuously graded meshes.

Astley, R. J.↗

NASA NTRS1994

Buckling and vibration analysis of laminated panels using VICONOPT

The analysis aspects of the 23,000-line FORTRAN program VIPASA with constraints and optimization (VICONOPT) are described. Overall stiffness matrices assembled from the earlier exact vibration and instability of plate assemblies including shear and anisotropy (VIPASA) flat plate stiffnesses are optionally coupled by Langrangian multipliers to find critical buckling loads, or natural frequencies of undamped vibration, of prismatic assemblies of anisotropic flat plates with arbitrarily located point supports or simple transverse supporting frames. The longitudunal continuity of typical wing and fuselage panels is closely approximated because the solutions are for the infinitely long structure obtained by repeating a bay and its supports longitudinally. Any longitudinally invariant in-plane plate stresses are permitted, and very rapid solutions are guaranteed by numerous refinements, including multilevel substructuring and a method for repetitive cross sections that is exact for regular polygons used to represent cylinders. Modal displacements and stresses in or between plies of laminated plates are calculated and plotted, with values being recovered at all nodes of substructures. Comparison with usual approximate finite-element methods confirms that, for comparably converged solutions, VICONOPT is typically between 100 and 10(exp 4) times faster.

Kennedy, David↗

NASA NTRS1999

Development of Fuzzy Logic and Soft Computing Methodologies

Our earlier research on computing with words (CW) has led to a new direction in fuzzy logic which points to a major enlargement of the role of natural languages in information processing, decision analysis and control. This direction is based on the methodology of computing with words and embodies a new theory which is referred to as the computational theory of perceptions (CTP). An important feature of this theory is that it can be added to any existing theory - especially to probability theory, decision analysis, and control - and enhance the ability of the theory to deal with real-world problems in which the decision-relevant information is a mixture of measurements and perceptions. The new direction is centered on an old concept - the concept of a perception - a concept which plays a central role in human cognition. The ability to reason with perceptions perceptions of time, distance, force, direction, shape, intent, likelihood, truth and other attributes of physical and mental objects - underlies the remarkable human capability to perform a wide variety of physical and mental tasks without any measurements and any computations. Everyday examples of such tasks are parking a car, driving in city traffic, cooking a meal, playing golf and summarizing a story. Perceptions are intrinsically imprecise. Imprecision of perceptions reflects the finite ability of sensory organs and ultimately, the brain, to resolve detail and store information. More concretely, perceptions are both fuzzy and granular, or, for short, f-granular. Perceptions are f-granular in the sense that: (a) the boundaries of perceived classes are not sharply defined; and (b) the elements of classes are grouped into granules, with a granule being a clump of elements drawn together by indistinguishability, similarity. proximity or functionality. F-granularity of perceptions may be viewed as a human way of achieving data compression. In large measure, scientific progress has been, and continues to be, driven by a quest to progress from perceptions to measurements. Pursuit of this aim has led to brilliant successes. But alongside the successes stand problems whose solutions are not in sight. Representative of such problems is the problem of automation of driving in city traffic. In this case, as in many others, what can be done with ease by humans - without any measurements and a computations - is an intractable task for machines.

Zadeh, L. A.↗

NASA NTRS2008

Integrated Force Method for Indeterminate Structures

Two methods of solving indeterminate structural-mechanics problems have been developed as products of research on the theory of strain compatibility. In these methods, stresses are considered to be the primary unknowns (in contrast to strains and displacements being considered as the primary unknowns in some prior methods). One of these methods, denoted the integrated force method (IFM), makes it possible to compute stresses, strains, and displacements with high fidelity by use of modest finite-element models that entail relatively small amounts of computation. The other method, denoted the completed Beltrami Mitchell formulation (CBMF), enables direct determination of stresses in an elastic continuum with general boundary conditions, without the need to first calculate displacements as in traditional methods. The equilibrium equation, the compatibility condition, and the material law are the three fundamental concepts of the theory of structures. For almost 150 years, it has been commonly supposed that the theory is complete. However, until now, the understanding of the compatibility condition remained incomplete, and the compatibility condition was confused with the continuity condition. Furthermore, the compatibility condition as applied to structures in its previous incomplete form was inconsistent with the strain formulation in elasticity.

Hopkins, Dale A.↗

NASA NTRS2014

A Multiscale Progressive Failure Modeling Methodology for Composites that Includes Fiber Strength Stochastics

A multiscale modeling methodology was developed for continuous fiber composites that incorporates a statistical distribution of fiber strengths into coupled multiscale micromechanics/finite element (FE) analyses. A modified two-parameter Weibull cumulative distribution function, which accounts for the effect of fiber length on the probability of failure, was used to characterize the statistical distribution of fiber strengths. A parametric study using the NASA Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC) was performed to assess the effect of variable fiber strengths on local composite failure within a repeating unit cell (RUC) and subsequent global failure. The NASA code FEAMAC and the ABAQUS finite element solver were used to analyze the progressive failure of a unidirectional SCS-6/TIMETAL 21S metal matrix composite tensile dogbone specimen at 650 degC. Multiscale progressive failure analyses were performed to quantify the effect of spatially varying fiber strengths on the RUC-averaged and global stress-strain responses and failure. The ultimate composite strengths and distribution of failure locations (predominately within the gage section) reasonably matched the experimentally observed failure behavior. The predicted composite failure behavior suggests that use of macroscale models that exploit global geometric symmetries are inappropriate for cases where the actual distribution of local fiber strengths displays no such symmetries. This issue has not received much attention in the literature. Moreover, the model discretization at a specific length scale can have a profound effect on the computational costs associated with multiscale simulations.models that yield accurate yet tractable results.

Generalized Method of Cells↗

NASA NTRS1980

Calculation of residual principal stresses in CVD boron on carbon filaments

A three-dimensional finite element model of the chemical vapor deposition (CVD) of boron on a carbon substrate (B/C) is developed. The model includes an expansion of the boron after deposition due to atomic rearrangement and includes creep of the boron and carbon. Curves are presented to show how the principal residual stresses and the filament elongation vary as the parameters defining deposition strain and creep are varied. The calculated results are compared with experimental axial residual stress and elongation measurements made on B/C filaments. This comparison requires that for good agreement between calculated and experimental results, the deposited boron must continue to expand after deposition, and that the build-up of residual stresses is limited by significant boron and carbon creep rates.

Behrendt, D. R.↗

NASA NTRS1993

Integrated analysis and design of thick composite structures for optimal passive damping characteristics

The development of novel composite mechanics for the analysis of damping in composite laminates and structures and the more significant results of this effort are summarized. Laminate mechanics based on piecewise continuous in-plane displacement fields are described that can represent both intralaminar stresses and interlaminar shear stresses and the associated effects on the stiffness and damping characteristics of a composite laminate. Among other features, the mechanics can accurately model the static and damped dynamic response of either thin or thick composite laminates, as well as, specialty laminates with embedded compliant damping layers. The discrete laminate damping theory is further incorporated into structural analysis methods. In this context, an exact semi-analytical method for the simulation of the damped dynamic response of composite plates was developed. A finite element based method and a specialty four-node plate element were also developed for the analysis of composite structures of variable shape and boundary conditions. Numerous evaluations and applications demonstrate the quality and superiority of the mechanics in predicting the damped dynamic characteristics of composite structures. Finally, additional development was focused on the development of optimal tailoring methods for the design of thick composite structures based on the developed analytical capability. Applications on composite plates illustrated the influence of composite mechanics in the optimal design of composites and the potential for significant deviations in the resultant designs when more simplified (classical) laminate theories are used.

Saravanos, D. A.↗

NASA NTRS1999

Open-Mode Debonding Analysis of Curved Sandwich Panels Subjected to Heating and Cryogenic Cooling on Opposite Faces

Increasing use of curved sandwich panels as aerospace structure components makes it vital to fully understand their thermostructural behavior and identify key factors affecting the open-mode debonding failure. Open-mode debonding analysis is performed on a family of curved honeycomb-core sandwich panels with different radii of curvature. The curved sandwich panels are either simply supported or clamped, and are subjected to uniform heating on the convex side and uniform cryogenic cooling on the concave side. The finite-element method was used to study the effects of panel curvature and boundary condition on the open-mode stress (radial tensile stress) and displacement fields in the curved sandwich panels. The critical stress point, where potential debonding failure could initiate, was found to be at the midspan (or outer span) of the inner bonding interface between the sandwich core and face sheet on the concave side, depending on the boundary condition and panel curvature. Open-mode stress increases with increasing panel curvature, reaching a maximum value at certain high curvature, and then decreases slightly as the panel curvature continues to increase and approach that of quarter circle. Changing the boundary condition from simply supported to clamped reduces the magnitudes of open-mode stresses and the associated sandwich core depth stretching.

Ko, William L.↗