Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element problems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,045 records · Page 58

Prediction of sound fields in cavities using boundary element methods

Two boundary element formulations of acoustical behavior exist; the Direct Boundary Element Methods (DBEM) based on the Helmholtz Integral Equations and the Indirect Boundary Element Methods (IBEM) based on Huygen's Principle. In this investigation, both methods are implemented utilizing a simple linear superparametric element. In addition the IBEM is studied using a quadratic isoparametric element. The accuracy and relative efficiency of the various techniques are examined. In order to properly model aircraft interior cavities the additional capability to model wall treatments and internal point sources is added to the methods. The procedures are verified for several well-understood cavity problems. The relative merits of each boundary element method and the finite element method are examined.

Bernhard, R. J.↗

Numerical simulation of electrophoresis separation processes

A new Petrov-Galerkin finite element formulation has been proposed for transient convection-diffusion problems. Most Petrov-Galerkin formulations take into account the spatial discretization, and the weighting functions so developed give satisfactory solutions for steady state problems. Though these schemes can be used for transient problems, there is scope for improvement. The schemes proposed here, which consider temporal as well as spatial discretization, provide improved solutions. Electrophoresis, which involves the motion of charged entities under the influence of an applied electric field, is governed by equations similiar to those encountered in fluid flow problems, i.e., transient convection-diffusion equations. Test problems are solved in electrophoresis and fluid flow. The results obtained are satisfactory. It is also expected that these schemes, suitably adapted, will improve the numerical solutions of the compressible Euler and the Navier-Stokes equations.

Ganjoo, D. K.↗

A survey of mixed finite element methods

This paper is an introduction to and an overview of mixed finite element methods. It discusses the mixed formulation of certain basic problems in elasticity and hydrodynamics. It also discusses special techniques for solving the discrete problem.

Brezzi, F.↗

Optimal cure cycle design of a resin-fiber composite laminate

A unified computed aided design method was studied for the cure cycle design that incorporates an optimal design technique with the analytical model of a composite cure process. The preliminary results of using this proposed method for optimal cure cycle design are reported and discussed. The cure process of interest is the compression molding of a polyester which is described by a diffusion reaction system. The finite element method is employed to convert the initial boundary value problem into a set of first order differential equations which are solved simultaneously by the DE program. The equations for thermal design sensitivities are derived by using the direct differentiation method and are solved by the DE program. A recursive quadratic programming algorithm with an active set strategy called a linearization method is used to optimally design the cure cycle, subjected to the given design performance requirements. The difficulty of casting the cure cycle design process into a proper mathematical form is recognized. Various optimal design problems are formulated to address theses aspects. The optimal solutions of these formulations are compared and discussed.

Hou, Jean W.↗

Dynamic response of laminated composite plates using a three-dimensional hybrid-stress finite-element formulation

A method of analysis of dynamic response of laminated composite plates is presented. The analysis is carried by using a hybrid-stress finite element numerical technique. By using this approach, the response of simply supported laminated plates subjected to sinusoidal loading are investigated. For the solution of the finite element equations of motion of free vibrations and dynamic response problems, two effective methods of solution, the space iteration method and the Newmark direct integration method are used. These two methods are discussed here.

Liou, W. J.↗

Computing Deformations Of Rubbery Materials

Better use made of experimental data in finite-element computations. New formulation of constitutive equations of rubbery, nonlinearly elastic material enables finite-element analysis of boundary-value stress-and-strain problems involving arbitrary shapes and loads. In development of formulation, principal stretches used as arguments of strain-energy-density function.

Peng, Steven T. J.↗

Elevated temperature crack growth

A study was performed to examine the applicability of path-independent (P-I) integrals to crack growth problems in hot section components of gas turbine aircraft engines. Alloy 718 was used and the experimental parameters included combined temperature and strain cycling, thermal gradients, elastic-plastic strain levels, and mean strains. A literature review was conducted of proposed P-I integrals, and those capable of analyzing hot section component problems were selected and programmed into the postprocessor of a finite element code. Detailed elastic-plastic finite element analyses were conducted to simulate crack growth and crack closure of the test specimen, and to evaluate the P-I integrals. It was shown that the selected P-I integrals are very effective for predicting crack growth for isothermal conditions.

Kim, K. S.↗

Probabilistic Structures Analysis Methods (PSAM) for select space propulsion system components

The basic formulation for probabilistic finite element analysis is described and demonstrated on a few sample problems. This formulation is based on iterative perturbation that uses the factorized stiffness on the unperturbed system as the iteration preconditioner for obtaining the solution to the perturbed problem. This approach eliminates the need to compute, store and manipulate explicit partial derivatives of the element matrices and force vector, which not only reduces memory usage considerably, but also greatly simplifies the coding and validation tasks. All aspects for the proposed formulation were combined in a demonstration problem using a simplified model of a curved turbine blade discretized with 48 shell elements, and having random pressure and temperature fields with partial correlation, random uniform thickness, and random stiffness at the root.

Source record↗

Nonlinear static and dynamic finite element analysis of an eccentrically loaded graphite-epoxy beam

The Dynamic Crash Analysis of Structures (DYCAT) and NIKE3D nonlinear finite element codes were used to model the static and implulsive response of an eccentrically loaded graphite-epoxy beam. A 48-ply unidirectional composite beam was tested under an eccentric axial compressive load until failure. This loading configuration was chosen to highlight the capabilities of two finite element codes for modeling a highly nonlinear, large deflection structural problem which has an exact solution. These codes are currently used to perform dynamic analyses of aircraft structures under impact loads to study crashworthiness and energy absorbing capabilities. Both beam and plate element models were developed to compare with the experimental data using the DYCAST and NIKE3D codes.

Fasanella, Edwin L.↗

Parallel aeroelastic computations for wing and wing-body configurations

The objective of this research is to develop computationally efficient methods for solving fluid-structural interaction problems by directly coupling finite difference Euler/Navier-Stokes equations for fluids and finite element dynamics equations for structures on parallel computers. This capability will significantly impact many aerospace projects of national importance such as Advanced Subsonic Civil Transport (ASCT), where the structural stability margin becomes very critical at the transonic region. This research effort will have direct impact on the High Performance Computing and Communication (HPCC) Program of NASA in the area of parallel computing.

Byun, Chansup↗

Modeling Composite Laminate Crushing for Crash Analysis

Crash modeling of composite structures remains limited in application and has not been effectively demonstrated as a predictive tool. While the global response of composite structures may be well modeled, when composite structures act as energy-absorbing members through direct laminate crushing the modeling accuracy is greatly reduced. The most efficient composite energy absorbing structures, in terms of energy absorbed per unit mass, are those that absorb energy through a complex progressive crushing response in which fiber and matrix fractures on a small scale dominate the behavior. Such failure modes simultaneously include delamination of plies, failure of the matrix to produce fiber bundles, and subsequent failure of fiber bundles either in bending or in shear. In addition, the response may include the significant action of friction, both internally (between delaminated plies or fiber bundles) or externally (between the laminate and the crushing surface). A figure shows the crushing damage observed in a fiberglass composite tube specimen, illustrating the complexity of the response. To achieve a finite element model of such complex behavior is an extremely challenging problem. A practical crushing model based on detailed modeling of the physical mechanisms of crushing behavior is not expected in the foreseeable future. The present research describes attempts to model composite crushing behavior using a novel hybrid modeling procedure. Experimental testing is done is support of the modeling efforts, and a test specimen is developed to provide data for validating laminate crushing models.

Fleming, David C.↗

A Capability to Generate Physics-based Mass Estimating Relationships for Conceptual Space Vehicle Design

This paper is written in support of the on-going research into conceptual space vehicle design conducted at the Space Systems Design Laboratory (SSDL) at the Georgia Institute of Technology. Research at the SSDL follows a sequence of a number of the traditional aerospace disciplines. The sequence of disciplines and interrelationship among them is shown in the Design Structure Matrix (DSM). The discipline of Weights and Sizing occupies a central location in the design of a new space vehicle. Weights and Sizing interact, either in a feed forward or feed back manner, with every other discipline in the DSM. Because of this principle location, accuracy in Weights and Sizing is integral to producing an accurate model of a space vehicle concept. Instead of using conceptual level techniques, a simplified Finite Element Analysis (FEA) technique is described as applied to the problem of the Liquid Oxygen (LOX) tank bending loads applied to the forward Liquid Hydrogen (LH2) tank of the Georgia Tech Air Breathing Launch Vehicle (ABLV).

Olds, John R.↗

GVIPS Models and Software

Two reports discuss, respectively, (1) the generalized viscoplasticity with potential structure (GVIPS) class of mathematical models and (2) the Constitutive Material Parameter Estimator (COMPARE) computer program. GVIPS models are constructed within a thermodynamics- and potential-based theoretical framework, wherein one uses internal state variables and derives constitutive equations for both the reversible (elastic) and the irreversible (viscoplastic) behaviors of materials. Because of the underlying potential structure, GVIPS models not only capture a variety of material behaviors but also are very computationally efficient. COMPARE comprises (1) an analysis core and (2) a C++-language subprogram that implements a Windows-based graphical user interface (GUI) for controlling the core. The GUI relieves the user of the sometimes tedious task of preparing data for the analysis core, freeing the user to concentrate on the task of fitting experimental data and ultimately obtaining a set of material parameters. The analysis core consists of three modules: one for GVIPS material models, an analysis module containing a specialized finite-element solution algorithm, and an optimization module. COMPARE solves the problem of finding GVIPS material parameters in the manner of a design-optimization problem in which the parameters are the design variables.

Arnold, Steven M.↗

Coupled finite-difference/finite-element approach for wing-body aeroelasticity

Computational methods using finite-difference approaches for fluids and finite-element approaches for structures have individually advanced to solve almost full-aircraft configurations. However, coupled approaches to solve fluid/structural interaction problems are still in their early stages of development, particularly for complex geometries using complete equations such as the Euler/Navier-Stokes equations. Earlier work demonstrated the success of coupling finite-difference and finite-element methods for simple wing configurations using the Euler/Navier-Stokes equations. In this paper, the same approach is extended for general wing-body configurations. The structural properties are represented by beam-type finite elements. The flow is modeled using the Euler/Navier-Stokes equations. A general procedure to fully couple structural finite-element boundary conditions with fluid finite-difference boundary conditions is developed for wing-body configurations. Computations are made using moving grids that adapt to wing-body structural deformations. Results are illustrated for a typical wing-body configuration.

Guruswamy, Guru P.↗

Practical implementation of an accurate method for multilevel design sensitivity analysis

Solution techniques for handling large scale engineering optimization problems are reviewed. Potentials for practical applications as well as their limited capabilities are discussed. A new solution algorithm for design sensitivity is proposed. The algorithm is based upon the multilevel substructuring concept to be coupled with the adjoint method of sensitivity analysis. There are no approximations involved in the present algorithm except the usual approximations introduced due to the discretization of the finite element model. Results from the six- and thirty-bar planar truss problems show that the proposed multilevel scheme for sensitivity analysis is more effective (in terms of computer incore memory and the total CPU time) than a conventional (one level) scheme even on small problems. The new algorithm is expected to perform better for larger problems and its applications on the new generation of computer hardwares with 'parallel processing' capability is very promising.

Nguyen, Duc T.↗

A finite element procedure for radio-frequency sheath–plasma interactions based on a sheath impedance model

A finite element code that solves self-consistent radio-frequency (RF) sheath-plasma interaction problems is improved by incorporating a generalized sheath boundary condition in the macroscopic solution scheme. This sheath boundary condition makes use of a complex sheath impedance including both the sheath capacitance and resistance, which enables evaluation of not only the RF voltage across the sheath but also the power dissipation in the sheath. The newly developed finite element procedure is applied to cases where the background magnetic field is perpendicular to the sheath surface in one- and two-dimensional domains filled by uniform low- and high-density plasmas. The numerical results are compared with those obtained by employing the previous capacitive sheath model at a typical frequency for ion cyclotron heating used in fusion experiments. It is shown that for sheaths on the order of 100 V in a high-density plasma, localized RF power deposition can reach a level which causes material damage. It is also shown that the sheath-plasma wave resonances predicted by the capacitive sheath model do not occur when parameters are such that the generalized sheath impedance model substantially modifies the capacitive character of the sheath. Here, possible explanations for the difference in the maximum RF sheath voltage depending on the plasma density are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of the shifted boundary method for the Stokes problem

The analysis of stability and accuracy of the shifted boundary method is developed for the Stokes flow equations. The key feature of the shifted boundary method, an embedded finite element method, is the shifting of the location where boundary conditions are applied from the true to a surrogate boundary, and an appropriate modification (shifting) of the value of the boundary conditions. An inf–sup condition is proved for the variational formulation associated to the shifted boundary method and we derive, by way of Strang’s second lemma, an optimal error estimate in the natural SBM norm. Here we also derive an L 2 -error estimate for the velocity field, by means of an extension of the Aubin–Nitsche approach to embedded, non-consistent, mixed finite element methods.

97 MATHEMATICS AND COMPUTING↗

Cost Considerations in Nonlinear Finite-Element Computing

Conference paper discusses computational requirements for finiteelement analysis using quasi-linear approach to nonlinear problems. Paper evaluates computational efficiency of different computer architecturtural types in terms of relative cost and computing time.

Utku, S.↗