Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

A finite element study of the EIDI system

This paper presents a method for modeling the structural dynamics of an Electro-Impulse De-Icing System, using finite element analyses procedures. A guideline for building a representative finite element model is discussed. Modeling was done initially using four noded cubic elements, four noded isoparametric plate elements and eight noded isoparametric shell elements. Due to the size of the problem and due to the underestimation of shear stress results when compared to previous analytical work an approximate model was created to predict possible areas of shedding of ice. There appears to be good agreement with the test data provided by The Boeing Commercial Airplane Company. Thus these initial results of this method were found to be encouraging. Additional analytical work and comparison with experiment is needed in order to completely evaluate this approach.

Khatkhate, A. A.↗
DOE OSTI2017

Finite-element modelling of multilayer X-ray optics

Multilayer optical elements for hard X-rays are an attractive alternative to crystals whenever high photon flux and moderate energy resolution are required. Prediction of the temperature, strain and stress distribution in the multilayer optics is essential in designing the cooling scheme and optimizing geometrical parameters for multilayer optics. The finite-element analysis (FEA) model of the multilayer optics is a well established tool for doing so. Multilayers used in X-ray optics typically consist of hundreds of periods of two types of materials. The thickness of one period is a few nanometers. Most multilayers are coated on silicon substrates of typical size 60 mm × 60 mm × 100–300 mm. The high aspect ratio between the size of the optics and the thickness of the multilayer (10 7 ) can lead to a huge number of elements for the finite-element model. For instance, meshing by the size of the layers will require more than 10 16 elements, which is an impossible task for present-day computers. Conversely, meshing by the size of the substrate will produce a too high element shape ratio (element geometry width/height > 10 6 ), which causes low solution accuracy; and the number of elements is still very large (10 6 ). In this work, by use of ANSYS layer-functioned elements, a thermal-structural FEA model has been implemented for multilayer X-ray optics. The possible number of layers that can be computed by presently available computers is increased considerably.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Finite Element Modeling of Orthogonal Machining of Brittle Materials Using an Embedded Cohesive Element Mesh

Machining of brittle materials is common in the manufacturing industry, but few modeling techniques are available to predict materials’ behavior in response to the cutting tool. The paper presents a fracture-based finite element model, named embedded cohesive zone–finite element method (ECZ–FEM). In ECZ–FEM, a network of cohesive zone (CZ) elements are embedded in the material body with regular elements to capture multiple randomized cracks during a cutting process. The CZ element is defined by the fracture energy and a scaling factor to control material ductility and chip behavior. The model is validated by an experimental study in terms of chip formation and cutting force with two different brittle materials and depths of cut. The results show that ECZ–FEM can capture various chip forms, such as dusty debris, irregular chips, and unstable crack propagation seen in the experimental cases. For the cutting force, the model can predict the relative difference among the experimental cases, but the force value is higher by 30–50%. The ECZ–FEM has demonstrated the feasibility of brittle cutting simulation with some limitations applied.

Takabi, Behrouz↗
NASA NTRS1992

Finite element simulations involving simultaneous multiple interface fronts in phase change problems

The present paper describes the simulation of phase change problems involving simultaneous multiple interface fronts employing the finite element method. Much of the past investigations employing finite elements have been restricted to primarily a single phase change situation. The existence of more than one phase, that is, the presence of multiple phase fronts poses certain challenges and further complications. However, the results provide a very interesting thermal behavior for this class of problems. In this paper, attention is focused on fixed grid methods and the trapezoidal family of one-step methods using the enthalpy formulations. Illustrative examples which handle simultaneous multiple fronts in phase change problems are presented.

Ouyang, Tianhong↗
NASA NTRS1982

Stability analysis of flexible wind turbine blades using finite element method

Static vibration and flutter analysis of a straight elastic axis blade was performed based on a finite element method solution. The total potential energy functional was formulated according to linear beam theory. The inertia and aerodynamic loads were formulated according to the blade absolute acceleration and absolute velocity vectors. In vibration analysis, the direction of motion of the blade during the first out-of-lane and first in-plane modes was examined; numerical results involve NASA/DOE Mod-0, McCauley propeller, north wind turbine and flat plate behavior. In flutter analysis, comparison cases were examined involving several references. Vibration analysis of a nonstraight elastic axis blade based on a finite element method solution was performed in a similar manner with the straight elastic axis blade, since it was recognized that a curved blade can be approximated by an assembly of a sufficient number of straight blade elements at different inclinations with respect to common system of axes. Numerical results involve comparison between the behavior of a straight and a curved cantilever beam during the lowest two in-plane and out-of-plane modes.

Kamoulakos, A.↗
NASA NTRS2001

The Use of Non-Standard Devices in Finite Element Analysis

A general mathematical description of the response behavior of thin-skin pneumatic envelopes and many other membrane and cable structures produces under-constrained systems that pose severe difficulties to analysis. These systems are mobile, and the general mathematical description exposes the mobility. Yet the response behavior of special under-constrained structures under special loadings can be accurately predicted using a constrained mathematical description. The static response behavior of systems that are infinitesimally mobile, such as a non-slack membrane subtended from a rigid or elastic boundary frame, can be easily analyzed using such general mathematical description as afforded by the non-linear, finite element method using an implicit solution scheme if the incremental uploading is guided through a suitable path. Similarly, if such structures are assembled with structural lack of fit that provides suitable self-stress, then dynamic response behavior can be predicted by the non-linear, finite element method and an implicit solution scheme. An explicit solution scheme is available for evolution problems. Such scheme can be used via the method of dynamic relaxation to obtain the solution to a static problem. In some sense, pneumatic envelopes and many other compliant structures can be said to have destiny under a specified loading system. What that means to the analyst is that what happens on the evolution path of the solution is irrelevant as long as equilibrium is achieved at destiny under full load and that the equilibrium is stable in the vicinity of that load. The purpose of this paper is to alert practitioners to the fact that non-standard procedures in finite element analysis are useful and can be legitimate although they burden their users with the requirement to use special caution. Some interesting findings that are useful to the US Scientific Balloon Program and that could not be obtained without non-standard techniques are presented.

Schur, Willi W.↗
NASA NTRS1999

Reflections of AE Waves in Finite Plates: Finite Element Modeling and Experimental Measurements

The capability of a three-dimensional dynamic finite element method for predicting far-field acoustic emission signals in thin plates of finite lateral extent, including their reflections from the plate edges, was investigated. A lead break (Hsu-Neilsen) source to simulate AE was modeled and used in the experimental measurements. For the thin plate studied, the signals were primarily composed of the lowest order symmetric (S0) and antisymmetric (A0) Lamb modes. Experimental waveforms were detected with an absolutely calibrated, wideband, conical element transducer. The conditions of lead fractures both on the surface of the plate as well as on the edge of the plate were investigated. Surface lead breaks preferentially generate the A0 mode while edge lead breaks generate the S0 mode. Reflections of developed plate waves from both normal and oblique incidence angles were evaluated. Particularly interesting for the case of the lead break on the plate edge were S0 waves produced by the interaction of a Rayleigh wave with the plate corner and by a bulk shear wave mode converting at the side edge. The Rayleigh wave, in this case, propagated along the specimen edge. For all cases considered, the experimental measurements were in good agreement with the predictions of the finite element model.

Prosser, W. H.↗
NASA NTRS1992

Improved accuracy for finite element structural analysis via a new integrated force method

A comparative study was carried out to determine the accuracy of finite element analyses based on the stiffness method, a mixed method, and the new integrated force and dual integrated force methods. The numerical results were obtained with the following software: MSC/NASTRAN and ASKA for the stiffness method; an MHOST implementation method for the mixed method; and GIFT for the integrated force methods. The results indicate that on an overall basis, the stiffness and mixed methods present some limitations. The stiffness method generally requires a large number of elements in the model to achieve acceptable accuracy. The MHOST method tends to achieve a higher degree of accuracy for course models than does the stiffness method implemented by MSC/NASTRAN and ASKA. The two integrated force methods, which bestow simultaneous emphasis on stress equilibrium and strain compatibility, yield accurate solutions with fewer elements in a model. The full potential of these new integrated force methods remains largely unexploited, and they hold the promise of spawning new finite element structural analysis tools.

Patnaik, Surya N.↗
DOE OSTI2022

A New Perspective of Post-Weld Baking Effect on Al-Steel Resistance Spot Weld Properties through Machine Learning and Finite Element Modeling

The root cause of post-weld baking on the mechanical performance of Al-steel dissimilar resistance spot welds (RSWs) has been determined by machine learning (ML) and finite element modeling (FEM) in this study. A deep neural network (DNN) model was constructed to associate the spot weld performance with the joint attributes, stacking materials, and other conditions, using a comprehensive experimental dataset. The DNN model positively identified that the post-weld baking reduces the joint performance, and the extent of degradation depends on the thickness of stacking materials. A three-dimensional finite element (FE) model was then used to investigate the root cause and the mechanism of the baking effect. It revealed that the formation of high thermal stresses during baking, from the mismatch of thermal expansion between steel and Al alloy, causes damage and cracking of the brittle intermetallic compound (IMC) formed at the interface of the weld nugget during welding. This in turn reduces the joint performance by promoting undesirable interfacial fracture when the welds were subjected to externally applied loads. The FEM model further revealed that increase in structural stiffness, because of increase in steel sheet thickness, reduces the thermal stresses at the interface caused by the thermal expansion mismatch and consequently lessens the detrimental effect of post-weld baking on the joint performance.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Discrete-ordinates finite-element method for atmospheric radiative transfer and remote sensing

The principal features of the discrete-ordinates finite-element method are reviewed, and the applicability of general-purpose discrete-ordinates codes to atmospheric radiative transfer and remote sensing problems is demonstrated. In particular, numerical results for typical problems arising in meteorology, climatology, and remote sensing are shown to be in good agreement with results from other methods and measurements. A sample two-dimensional calculation demonstrates that specific capabilities available in the discrete-ordinates code TWOTRAN can produce new results that are valuable in the characterization of atmospheric effects on remote sensing (e.g., the adjacency effect). The intrinsic limitations of the method are also considered, and it is concluded that the strengths of the discrete-ordinates finite-element method outweigh its weaknesses.

Gerstl, S. A. W.↗
NASA NTRS1998

A Posteriori Finite Element Bounds for Sensitivity Derivatives of Partial-Differential-Equation Outputs

We present a Neumann-subproblem a posteriori finite element procedure for the efficient and accurate calculation of rigorous, 'constant-free' upper and lower bounds for sensitivity derivatives of functionals of the solutions of partial differential equations. The design motivation for sensitivity derivative error control is discussed; the a posteriori finite element procedure is described; the asymptotic bounding properties and computational complexity of the method are summarized; and illustrative numerical results are presented.

Lewis, Robert Michael↗
NASA NTRS1996

Verification of Orthogrid Finite Element Modeling Techniques

The stress analysis of orthogrid structures, specifically with I-beam sections, is regularly performed using finite elements. Various modeling techniques are often used to simplify the modeling process but still adequately capture the actual hardware behavior. The accuracy of such 'Oshort cutso' is sometimes in question. This report compares three modeling techniques to actual test results from a loaded orthogrid panel. The finite element models include a beam, shell, and mixed beam and shell element model. Results show that the shell element model performs the best, but that the simpler beam and beam and shell element models provide reasonable to conservative results for a stress analysis. When deflection and stiffness is critical, it is important to capture the effect of the orthogrid nodes in the model.

Steeve, B. E.↗
DOE OSTI2017

Modal Substructuring of Geometrically Nonlinear Finite Element Models with Interface Reduction

Substructuring methods have been widely used in structural dynamics to divide large, complicated finite element models into smaller substructures. For linear systems, many methods have been developed to reduce the subcomponents down to a low order set of equations using a special set of component modes, and these are then assembled to approximate the dynamics of a large scale model. In this paper, a substructuring approach is developed for coupling geometrically nonlinear structures, where each subcomponent is drastically reduced to a low order set of nonlinear equations using a truncated set of fixedinterface and characteristic constraint modes. The method used to extract the coefficients of the nonlinear reduced order model (NLROM) is non-intrusive in that it does not require any modification to the commercial FEA code, but computes the NLROM from the results of several nonlinear static analyses. The NLROMs are then assembled to approximate the nonlinear differential equations of the global assembly. The method is demonstrated on the coupling of two geometrically nonlinear plates with simple supports at all edges. The plates are joined at a continuous interface through the rotational degrees-of-freedom (DOF), and the nonlinear normal modes (NNMs) of the assembled equations are computed to validate the models. The proposed substructuring approach reduces a 12,861 DOF nonlinear finite element model down to only 23 DOF, while still accurately reproducing the first three NNMs of the full order model.

42 ENGINEERING↗
NASA NTRS1982

Three-dimensional finite-element analysis of layered composite plates

Results are presented for an investigation of the three-dimensional, geometrically nonlinear, finite-element analysis of the bending of laminated anisotropic composite plates. The individual laminae are treated as homogeneous, transversely isotropic, and linearly elastic. A fully three-dimensional isoparametric finite element with eight modes (i.e., linear element) and 24 degrees of freedom (three displacement components per node) is used. The numerical results obtained using this linear analysis are compared with the exact solutions given in Pagano (1969, 1970). It is found that the results of the linear analysis converge to the exact solution as the mesh is refined.

Putcha, N. S.↗
NASA NTRS1993

Domain decomposition methods for nonconforming finite element spaces of Lagrange-type

In this article, we consider the application of three popular domain decomposition methods to Lagrange-type nonconforming finite element discretizations of scalar, self-adjoint, second order elliptic equations. The additive Schwarz method of Dryja and Widlund, the vertex space method of Smith, and the balancing method of Mandel applied to nonconforming elements are shown to converge at a rate no worse than their applications to the standard conforming piecewise linear Galerkin discretization. Essentially, the theory for the nonconforming elements is inherited from the existing theory for the conforming elements with only modest modification by constructing an isomorphism between the nonconforming finite element space and a space of continuous piecewise linear functions.

Cowsar, Lawrence C.↗
NASA NTRS2000

Finite Elements Analysis of a Composite Semi-Span Test Article With and Without Discrete Damage

AS&M Inc. performed finite element analysis, with and without discrete damage, of a composite semi-span test article that represents the Boeing 220-passenger transport aircraft composite semi-span test article. A NASTRAN bulk data file and drawings of the test mount fixtures and semi-span components were utilized to generate the baseline finite element model. In this model, the stringer blades are represented by shell elements, and the stringer flanges are combined with the skin. Numerous modeling modifications and discrete source damage scenarios were applied to the test article model throughout the course of the study. This report details the analysis method and results obtained from the composite semi-span study. Analyses were carried out for three load cases: Braked Roll, LOG Down-Bending and 2.5G Up-Bending. These analyses included linear and nonlinear static response, as well as linear and nonlinear buckling response. Results are presented in the form of stress and strain plots. factors of safety for failed elements, buckling loads and modes, deflection prediction tables and plots, and strainage prediction tables and plots. The collected results are presented within this report for comparison to test results.

Lovejoy, Andrew E.↗
DOE OSTI2018

Finite deformation cohesive polygonal finite elements for modeling pervasive fracture

In this work, we introduce a framework for modeling dynamic fracture problems using cohesive polygonal finite elements. Random polygonal meshes provide a robust, efficient method for generating an unbiased network of fracture surfaces. Further, these meshes have more facets per element than standard triangle or quadrilateral meshes, providing more possible facets per element to insert cohesive surfaces. This property of polygonal meshes is advantageous for the modeling of pervasive fracture. We use both Wachspress and maximum entropy shape functions to form a finite element basis over the polygons. Fracture surfaces are captured through dynamically inserted cohesive zone elements at facets between the polygons in the mesh. Contact is enforced through a penalty method that is applied to both closed cohesive surfaces and general interpenetration of two polygonal elements. Finally, several numerical examples are presented that illustrate the capabilities of the method and demonstrate convergence of solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Extended Finite Element Based Approach in Additive Manufacturing Modeling for Optimizing Highly Complex Manifold in Protonic Ceramic Electrochemical Cells

The objective of this project is to develop a novel numerical approach based on the extended finite element method (XFEM) for modeling moving boundaries of material deposited in additive manufacturing (AM) process with improved accuracy and reduced computational cost. One of the major challenges in simulating AM processes is that the boundaries of the computational domain need to change as material is deposited. Previous approaches have achieved this by activating new finite elements, but this requires a high level of mesh refinement to capture small movement of the boundary. XFEM allows the solution boundary to move independently of the mesh, permitting smooth representation of the evolution of the boundary as material is deposited. As a result, this new method will significantly improve the simulation accuracy while reducing the computational cost compared with existing approaches. Furthermore, this project will bring significant improvements for the design and AM process optimization in an iterative fashion among numerical simuation, parameter optimization, and experimental validation. Specifically, this project will support the protonic ceramic electrochemical cells (PCEC) stack development at INL by providing insights into the density, residual stress, thermal and mechanical properties of the PCEC manifold and interconnect, which is critical in enhancing the system lifetime and reducing the overall cost. Upon the success of proposed development and validation of the proposed approach, simulation will be applied to determine an optimal set of AM process parameters for the PCEC manifold production. This project will consolidate Idaho National Laboratory (INL)’s simulation capabilities provided by the MOOSE framework and the Valhalla AM simulation application, with encouraging expansion to emerging PCEC applications and AM technologies at INL and beyond.

97 MATHEMATICS AND COMPUTING↗