Search NASA⌕ Search

SEARCH · Search NASA

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 631 records · Page 35

Structural Model Tuning Capability in an Object-Oriented Multidisciplinary Design, Analysis, and Optimization Tool

Updating the finite element model using measured data is a challenging problem in the area of structural dynamics. The model updating process requires not only satisfactory correlations between analytical and experimental results, but also the retention of dynamic properties of structures. Accurate rigid body dynamics are important for flight control system design and aeroelastic trim analysis. Minimizing the difference between analytical and experimental results is a type of optimization problem. In this research, a multidisciplinary design, analysis, and optimization (MDAO) tool is introduced to optimize the objective function and constraints such that the mass properties, the natural frequencies, and the mode shapes are matched to the target data as well as the mass matrix being orthogonalized.

Lung, Shun-fat↗

Aeroelastic Modeling of X-56A Stiff-Wing Configuration Flight Test Data

Aeroelastic stability and control derivatives for the X-56A Multi-Utility Technology Testbed (MUTT), in the stiff-wing configuration, were estimated from flight test data using the output-error method. Practical aspects of the analysis are discussed. The orthogonal phase-optimized multisine inputs provided excellent data information for aeroelastic modeling. Consistent parameter estimates were determined using output error in both the frequency and time domains. The frequency domain analysis converged faster and was less sensitive to starting values for the model parameters, which was useful for determining the aeroelastic model structure and obtaining starting values for the time domain analysis. Including a modal description of the structure from a finite element model reduced the complexity of the estimation problem and improved the modeling results. Effects of reducing the model order on the short period stability and control derivatives were investigated.

Grauer, Jared A.↗

Modeling and simulation challenges for softgoods in space deployable structures

This paper examines modeling and simulation challenges for deployable spacecraft structures with softgood components. Previous studies highlighted the numerical challenges due to the mixture of high stiffness and low stiffness contact in this class of structures. This paper examines the use of the commercial LS-DYNA nonlinear finite element code on a set of relevant benchmark problems.

Mobrem, Mehran↗

A shock capturing application of the finite element method

The paper is concerned with the development of finite element algorithms for the solution of viscous compressible flow problems with possible embedded shocks and recirculation regions. As a first step, the calculation in Cartesian coordinates of uniform flow on a rectangular region which encounters an embedded oblique shock with known turning angle is considered. A code is developed which is then used for computation of the boattail plume simulator problem in cylindrical coordinates. Oblique shock calculations are performed and the results are compared with a known finite difference solution.

Cooke, C. H.↗

A finite-element alternating method for two-dimensional Mode-1 crack configurations

A finite-element alternating method is presented for 2-D Mode-1 crack problems. An analytical solution for an arbitrary polynomial normal pressure distribution applied to the crack faces is obtained and used as the basic solution in the method. The method is applied to several crack problems to study its efficiency and the results are compared to accurate stress-intensity factor solutions in the literature. The method gave reasonably accurate stress-intensity factors and crack opening displacements with minimal computing effort. Because the method must model only the uncracked body, finite-element models with many degrees of freedom are not warranted and therefore, the method has been implemented on personal computers.

Raju, I. S.↗

Dynamics of a 4x6-Meter Thin Film Elliptical Inflated Membrane for Space Applications

Dynamic characterization of a thin film inflatable elliptical structure is described in detail. A two-step finite element modeling approach in MSC/NASTRAN is utilized, consisting of (1) a nonlinear static pressurization procedure used to obtain the updated stiffness matrix, and (2) a modal "restart" eigen solution that uses the modified stiffness matrix. Unique problems encountered in modeling of this large Hexameter lightweight inflatable arc identified, including considerable difficulty in obtaining convergence in the nonlinear finite element pressurization solution. It was found that the extremely thin polyimide film material (.001 in or 1 mil) presents tremendous problems in obtaining a converged solution when internal pressure loading is applied. Approaches utilized to overcome these difficulties are described. Comparison of finite element predictions for frequency and mode shapes of the inflated structure with closed-form solutions for a flat pre-tensioned membrane indicate reasonable agreement.

Casiano, Matthew J.↗

Recent advances in thermostructural finite element analysis

A comparison study is described which has indicated that the finite element and lumped parameter methods are comparable for steady-state conduction/convection analyses. The finite element model was advantageous because of the ease of model verification with computer graphics. Some problems encountered in constructing compatible thermal and structural models have indicated the need for further methodology developments to automate the transition between thermal and structural analyses. The need for more efficient finite element nonlinear transient algorithms has been cited. With further development in the thermal analysis area, the finite element method offers high potential for an integrated thermal/structural analysis capability.

Thornton, E. A.↗

A Unified Development of Basis Reduction Methods for Rotor Blade Analysis

The axial foreshortening effect plays a key role in rotor blade dynamics, but approximating it accurately in reduced basis models has long posed a difficult problem for analysts. Recently, though, several methods have been shown to be effective in obtaining accurate,reduced basis models for rotor blades. These methods are the axial elongation method,the mixed finite element method, and the nonlinear normal mode method. The main objective of this paper is to demonstrate the close relationships among these methods, which are seemingly disparate at first glance. First, the difficulties inherent in obtaining reduced basis models of rotor blades are illustrated by examining the modal reduction accuracy of several blade analysis formulations. It is shown that classical, displacement-based finite elements are ill-suited for rotor blade analysis because they can't accurately represent the axial strain in modal space, and that this problem may be solved by employing the axial force as a variable in the analysis. It is shown that the mixed finite element method is a convenient means for accomplishing this, and the derivation of a mixed finite element for rotor blade analysis is outlined. A shortcoming of the mixed finite element method is that is that it increases the number of variables in the analysis. It is demonstrated that this problem may be rectified by solving for the axial displacements in terms of the axial forces and the bending displacements. Effectively, this procedure constitutes a generalization of the widely used axial elongation method to blades of arbitrary topology. The procedure is developed first for a single element, and then extended to an arbitrary assemblage of elements of arbitrary type. Finally, it is shown that the generalized axial elongation method is essentially an approximate solution for an invariant manifold that can be used as the basis for a nonlinear normal mode.

Ruzicka, Gene C.↗

Face Sheet/Core Disbond Growth in Honeycomb Sandwich Panels Subjected to Ground-Air-Ground Pressurization and In-Plane Loading

Typical damage modes in light honeycomb sandwich structures include face sheet/core disbonding and core fracture, both of which can pose a threat to the structural integrity of a component. These damage modes are of particular interest to aviation certification authorities since several in-service occurrences, such as rudder structural failure and other control surface malfunctions, have been attributed to face sheet/core disbonding. Extensive studies have shown that face sheet/core disbonding and core fracture can lead to damage propagation caused by internal pressure changes in the core. The increasing use of composite sandwich construction in aircraft applications makes it vitally important to understand the effect of ground-air-ground (GAG) cycles and conditions such as maneuver and gust loads on face sheet/core disbonding. The objective of the present study was to use a fracture mechanics based approach developed earlier to evaluate the loading at the disbond front caused by ground-air-ground pressurization and in-plane loading. A honeycomb sandwich panel containing a circular disbond at one face sheet/core interface was modeled with three-dimensional (3D) solid finite elements. The disbond was modeled as a discrete discontinuity and the strain energy release rate along the disbond front was computed using the Virtual Crack Closure Technique (VCCT). Special attention was paid to the pressure-deformation coupling which can decrease the pressure load within the disbonded sandwich section significantly when the structure is highly deformed. The commercial finite element analysis software, Abaqus/Standard, was used for the analyses. The recursive pressure-deformation coupling problem was solved by representing the entrapped air in the honeycomb cells as filled cavities in Abaqus/Standard. The results show that disbond size, face sheet thickness and core thickness are important parameters that determine crack tip loading at the disbond front. Further, the pressure-deformation coupling was found to have an important load decreasing effect [6]. In this paper, a detailed problem description is provided first. Second, the analysis methodology is presented. The fracture mechanics approach used is described and the specifics of the finite element model, including the fluid-filled cavities, are introduced. Third, the initial model verification and validation are discussed. Fourth, the findings from a closely related earlier study [6] are summarized. These findings provided the basis for the current investigation. Fifth, an aircraft ascent scenario from 0 to 12192 m (0 to 40000 ft) is considered and the resulting crack tip loading at the disbond front is determined. In-plane loading to simulate maneuvers and gust conditions are also considered. Sixth, the results are shown for a curved panel, which was used to simulate potential fuselage applications. Finally, a brief summary of observations is presented and recommendations for improvement are provided.

Chen, Zhi M.↗

How Certain Physical Considerations Impact Aerostructural Wing Optimization

Wing design optimization has been studied extensively and is of continued interest as optimization tools are developed and become more accessible. In each of these studies, certain assumptions and simplifications are made to make the design problem tractable. However, it is difficult to find systematic studies in which several considerations are added or removed one at a time to study how much impact they have. In this work, we examine how certain physical considerations (viscous drag, wave drag, thrust loads, and inertial relief from structural, fuel, and engine masses), impact the aerostructural optimization results for three distinct aircraft wings. The goal is to help develop a rough idea of how important these physical considerations are. We do this using gradient-based optimization and a multidisciplinary design optimization framework, OpenMDAO. We use the open-source tool OpenAeroStruct that couples a vortex lattice method to a finite element method. We establish a baseline aerostructural design optimization problem then perform a series of optimizations, each with one physical consideration removed from the baseline case. We find that depending on the size of the aircraft and flight conditions, the importance of some of these physical considerations varies considerably whereas the importance of others do not. Specifically, the optimal designs change radically without proper viscous and wave drag considerations and smaller aircraft with more distributed propulsion are more affected by the inclusion of engine loads.

Aeropropulsion↗

Meshfree Methods

Meshfree methods have undergone substantial development and have received much attention in the last two decades. This new family of numerical methods is designed to inherit the main advantages of the finite element method such as compact supports of shape functions and good approximation properties while, at the same time, overcome the main disadvantages of the finite element method caused by the mesh dependence. The meshfree methods share a common feature that no mesh is needed and shape functions are constructed from sets of points, thus eliminating the need for time consuming mesh generation. The most significant advantage of meshfree methods is the flexibility in customizing approximation functions for desired regularity and for capturing essential physics and features of the particular problems of interest. Adaptivity formulation and multiple-scale solution strategies also can be implemented with relative ease. It has become clear that the meshfree methods provide considerable advantages over the conventional finite element methods in solving problems involving moving discontinuities, evolving material interfaces, multiple-scale phenomena, large material distortion and structural deformation, and fracture and damage processes. This Chapter gives an overview of many classes of meshfree methods, with more detailed discussions on Smoothed Particle Hydrodynamics (SPH), the Reproducing Kernel Particle Method (RKPM), Peridynamics (PD), the Material Point Method (MPM), as well as their applications in various challenging engineering problems.2

Chen, Jiun-Shyan↗

Iterative Nonlinear Tikhonov Algorithm with Constraints for Electromagnetic Tomography

Low frequency electromagnetic tomography such as the capacitance tomography (ECT) has been proposed for monitoring and mass-gauging of gas-liquid two-phase system under microgravity condition in NASA's future long-term space missions. Due to the ill-posed inverse problem of ECT, images reconstructed using conventional linear algorithms often suffer from limitations such as low resolution and blurred edges. Hence, new efficient high resolution nonlinear imaging algorithms are needed for accurate two-phase imaging. The proposed Iterative Nonlinear Tikhonov Regularized Algorithm with Constraints (INTAC) is based on an efficient finite element method (FEM) forward model of quasi-static electromagnetic problem. It iteratively minimizes the discrepancy between FEM simulated and actual measured capacitances by adjusting the reconstructed image using the Tikhonov regularized method. More importantly, it enforces the known permittivity of two phases to the unknown pixels which exceed the reasonable range of permittivity in each iteration. This strategy does not only stabilize the converging process, but also produces sharper images. Simulations show that resolution improvement of over 2 times can be achieved by INTAC with respect to conventional approaches. Strategies to further improve spatial imaging resolution are suggested, as well as techniques to accelerate nonlinear forward model and thus increase the temporal resolution.

Xu, Feng↗

Analysis of delamination in fiber composite laminates out-of-plane under bending

Delamination in the form of cracking or separation between plies in an advanced fiber composite laminate is a problem of major concern. Both advanced analytical methods and advanced computational analyses are conducted to: (1) develop an asymptotic solution for a composite laminate subject to out-of-plane bending; (2) construct advanced singular finite elements in conjunction with the development of nonsingular elements for this bending problem; and (3) evaluate the delamination failure mechanics parameters and the subsequent modes of fracture. A parametric study was also conducted to evaluate the influences of various lamination parameters on the delaminated composites.

Wang, S. S.↗

Multi-Body Analysis of a Tiltrotor Configuration

The paper describes the aeroelastic analysis of a tiltrotor configuration. The 1/5 scale wind tunnel semispan model of the V-22 tiltrotor aircraft is considered. The analysis is performed by means of a multi-body code, based on an original formulation. The differential equilibrium problem is stated in terms of first order differential equations. The equilibrium equations of every rigid body are written, together with the definitions of the momenta. The bodies are connected by kinematic constraints, applied in form of Lagrangian multipliers. Deformable components are mainly modelled by means of beam elements, based on an original finite volume formulation. Multi-disciplinar problems can be solved by adding user-defined differential equations. In the presented analysis the equations related to the control of the swash-plate of the model are considered. Advantages of a multi-body aeroelastic code over existing comprehensive rotorcraft codes include the exact modelling of the kinematics of the hub, the detailed modelling of the flexibility of critical hub components, and the possibility to simulate steady flight conditions as well as wind-up and maneuvers. The simulations described in the paper include: 1) the analysis of the aeroelastic stability, with particular regard to the proprotor/pylon instability that is peculiar to tiltrotors, 2) the determination of the dynamic behavior of the system and of the loads due to typical maneuvers, with particular regard to the conversion from helicopter to airplane mode, and 3) the stress evaluation in critical components, such as the pitch links and the conversion downstop spring.

Ghiringhelli, G. L.↗

Effects of local vibrations on the dynamics of space truss structures

The paper discusses the influence of local member vibrations on the dynamics of repetitive space truss structures. Several focus problems wherein local member vibration modes are in the frequency range of the global truss modes are discussed. Special attention is given to defining methods that can be used to identify the global modes of a truss structure amidst many local modes. Significant interactions between the motions of local member vibrations and the global behavior are shown to occur in truss structures when: (1) the natural frequencies of the individual members for clamped-clamped boundary conditions are in the vicinity of the global truss frequency; and (2) the total mass of the individual members represents a large portion of the mass of the whole structure. The analysis is carried out with a structural analysis code which uses exact member theory. The modeling detail required using conventional finite element codes to adequately represent such a class of problems is examined. The paper concludes with some practical considerations for the design and dynamic testing of structures which might exhibit such behavior.

Warnaar, Dirk B.↗

A high order accurate finite element algorithm for high Reynolds number flow prediction

A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.

Baker, A. J.↗

Three-dimensional finite strip analysis of elastic solids

Finite element shape functions were combined with beam eigenfunctions to formulate three-dimensional finite strips. Because of the orthogonality of the beam functions, three-dimensional problems were reduced to a series of two-dimensional problems, often with stiffness matrices of very narrow bandwidth. Isoparametric and high order finite element shape functions were used in the formulation of the 3-D finite strips. Numerical examples such as the static and free vibration analyses of simply supported thick plates are presented. Results were compared with existing solutions and good agreements were obtained in all cases.

Cheung, M. S.↗