Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Element Method (FEM)”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Simulating Fatigue Crack Growth in Spiral Bevel Pinion

This project investigates computational modeling of fatigue crack growth in spiral bevel gears. Current work is a continuation of the previous efforts made to use the Boundary Element Method (BEM) to simulate tooth-bending fatigue failure in spiral bevel gears. This report summarizes new results predicting crack trajectory and fatigue life for a spiral bevel pinion using the Finite Element Method (FEM). Predicting crack trajectories is important in determining the failure mode of a gear. Cracks propagating through the rim may result in catastrophic failure, whereas the gear may remain intact if one tooth fails and this may allow for early detection of failure. Being able to predict crack trajectories is insightful for the designer. However, predicting growth of three-dimensional arbitrary cracks is complicated due to the difficulty of creating three-dimensional models, the computing power required, and absence of closed- form solutions of the problem. Another focus of this project was performing three-dimensional contact analysis of a spiral bevel gear set incorporating cracks. These analyses were significant in determining the influence of change of tooth flexibility due to crack growth on the magnitude and location of contact loads. This is an important concern since change in contact loads might lead to differences in SIFs and therefore result in alteration of the crack trajectory. Contact analyses performed in this report showed the expected trend of decreasing tooth loads carried by the cracked tooth with increasing crack length. Decrease in tooth loads lead to differences between SIFs extracted from finite element contact analysis and finite element analysis with Hertz contact loads. This effect became more pronounced as the crack grew.

Ural, Ani↗
NASA NTRSDate not supplied

Developing Procedures to Implement Geometric Imperfections Beyond Right Circular Cylindrical Shells in Finite Element Method Models

Analysis of aerospace structures is frequently conducted using nominal dimensions and frequently assumes ideal conditions in loading, contact, constraints, et cetera. Off-nominal dimensions and nonideal conditions, however, are present in all structures. These are the result of widely ranging causes from coefficient of thermal expansion mismatches, manufacturing tooling anomalies, to assembly procedures that inadvertently alter the structure. Specifically, geometric imperfections can have potentially significant influence on the response of a structural test article observed in an experiment versus the response given by a numerical simulation. The Python Tool for Implementing Geometric Imperfections in Reduced Structures (Py_TIGIRS) was previously presented as a set of Python scripts to calculate and implement as-manufactured geometric midsurface and thickness imperfections into finite element method (FEM) shell models of nominally right circular cylinders. By taking advantage of the simple shape of a right circular cylinder, interpolations of the measured data points were able to be performed along directions that aligned to the cylindrical coordinate system axes of the entire structure. By taking advantage of the shell representation of the real structure as opposed to modeling using a continuum representation, the thickness variation was able to be implemented by shell section definitions instead of having to modify the position of multiple nodes in the thickness direction. Py_TIGIRS is a useful tool that established a procedural example on how to implement geometric imperfections in right circular cylindrical shell structures. Three new procedures, each expanded from concepts established in Py_TIGIRS, are proposed for various test-article designs and are intended to broaden the range of structures that can be modeled with measured geometric imperfections in the structural analysis community. Each test-article design introduces new challenges to successfully implement geometric imperfections into a FEM model. The first test-article design consists of a carbon fiber reinforced polymer square plate with a hat-shaped stiffener co-cured on one side. This test-article design was for a novel seven-point bend test that was also previously presented. Manufacturing and cure-cycle imperfections are observed using digital image correlation (DIC) techniques. As thermal expansion coefficient mismatches between the plate and stiffener materials were anticipated, a thermal analysis study with continuum shell and solid elements was conducted to capture the global shape observed prior to testing. The second test-article design is of a similar hat-stiffened plate configuration, but with a side length ratio near 3:1 with elongation in the stiffener direction. The test article was used to characterize the response to uniaxial compressive loading in the direction of the stiffener. Due to differing manufacturing steps, a thermal analysis like the one developed for the seven-point bend configuration was unable to mimic the observed geometric imperfections. Instead, a strategy based on applying deformations directly to the structure during analysis was developed for continuum shell and solid element representation of a stiffened panel.

Geometric imperfections↗
NASA NTRS1977

Transmission in nonuniform ducts - A comparative evaluation of finite element and weighted residuals computational schemes

The Method of Weighted Residuals (MWR) and the Finite Element Method (FEM) are considered as computational schemes in the problem of acoustic transmission in nonuniform ducts. MWR is presented in an improved form which includes the interaction of acoustic modes (irrotational) and hydrodynamic modes (rotational). FEM is based on a weighted residuals formulation using eight noded isoparametric elements. Both are applicable to two-dimensional and axially symmetric problems. Calculations are made for several sample problems to demonstrate accuracy and relative efficiency. One test case has implications in the phenomenon of subsonic acoustic choking and it is found that a large transmission loss is not an automatic consequence of propagation against a high subsonic mean flow.

Eversman, W.↗
NASA NTRS1979

The application of finite element techniques to acoustic transmission in lined ducts with flow

The finite element method (FEM) is used to analyze the propagation of sound in two-dimensional nonuniform ducts carrying a compressible subsonic mean flow. Galerkin and residual least squares (RLS) methods with natural and forced boundary conditions are considered. The accuracy of FEM results for the eigenvalue and transmission problems is assessed by comparison with alternative numerical schemes for nonuniform ducts. The results presented and those from associated investigations indicate that modal coupling is a significant feature of the acoustic field, especially at high Mach numbers. A multimodal model therefore appears to be essential if any reliable conclusions are to be drawn in the context of turbofan inlet regions. Improvements to the eigenvalue scheme following the implementation of higher-order Hermitian elements indicate a similar modification for the transmission problem.

Astley, R. J.↗
NASA NTRS1981

Accuracy and stability of finite element schemes for the duct transmission problem

An investigation is conducted regarding the feasibility of approaches for improving the efficiency and stability of existing finite element method (FEM) schemes, taking into account both analytical and numerical studies. Of the four schemes considered for the 'steady' problem, the Hermitian Galerkin formulation appears to be the most efficient and therefore the most suitable scheme for futher full scale implementation. The Hermitian residual least squares (RLS) scheme although comparable in accuracy for the cases considered exhibits a slight tendency to cumulative errors. The performance of both the Lagrangian element schemes considered compares poorly with that of their Hermitian element counterparts. This is particularly true of the Lagrangian RLS scheme. The presence of internal oscillatory components is an inevitable consequence of all Galerkin schemes irrespective of element type.

Astley, R. J.↗
NASA NTRS1984

Development of large rotorcraft transmissions

The U.S. Army Heavy Lift Helicopter (HLH) represents a large rotorcraft which was developed by an American aerospace company. In the early 1970's with the HLH Advanced Technology Components (ATC) program, the development of large rotorcraft transmission and drive systems was started. Failures in the spiral bevel gearing were experienced in tests because the employed method of analysis had not considered the effect of rim bending. Consequently, new gears with strengthened rims were designed and fabricated. For a more accurate prediction of the load capacity of the gears, an extensive Finite Element Method (FEM) system was developed. The U.S. Army's XCH-62 HLH aft rotor transmission was finally successfully tested at full design torque and speed. A description of the test program is provided, and the analytical program is discussed. The analytical phase includes the development of a preprocessing program which aids in the review of calculated FEM stresses.

Samanich, N. E.↗
NASA NTRS1983

A finite element model for sound transmission through panels

A finite element method (FEM) is applied to predicting coupled frequencies and pressures within an acoustic cavity in order to characterize sound transmission through a panel. Structural equations of motion are defined and the FEM model is configured with four-noded plate elements, the lateral displacement and two slopes being the unknowns at every node. Each element then has 12 degrees of freedom (DOF) and the displacement variation is expressed by a 12-term nonconforming polynomial. A consistent mass matrix is used to represent the panel mass matrix and a wave equation governs the acoustic volume. Analysis of pressure and displacement over the panel yields a square coupling matrix, and an eigenanalysis leads to a solution of the forced vibration problem.

Ramakrishnan, J. V.↗
NASA NTRS1985

Analysis of Spiral Bevel Gearing

Program generalized for automated finite-element-method model generation. Tedious task of finite-element-method (FEM) model preparation greatly reduced by development of preprocessing program for automated model generation. Utility of program, including ability to rapidly evaluate parametric changes, demonstrated for spiral, helical, and spiral bevel gearing. Postprocessing program also developed, aids in review and interpretation of calculated FEM stresses.

Drago, R. J.↗
NASA NTRS1986

Survey of the status of finite element methods for partial differential equations

The finite element methods (FEM) have proved to be a powerful technique for the solution of boundary value problems associated with partial differential equations of either elliptic, parabolic, or hyperbolic type. They also have a good potential for utilization on parallel computers particularly in relation to the concept of domain decomposition. This report is intended as an introduction to the FEM for the nonspecialist. It contains a survey which is totally nonexhaustive, and it also contains as an illustration, a report on some new results concerning two specific applications, namely a free boundary fluid-structure interaction problem and the Euler equations for inviscid flows.

Temam, Roger↗
NASA NTRS1987

Theoretical analysis of HVAC duct hanger systems

Several methods are presented which, together, may be used in the analysis of duct hanger systems over a wide range of frequencies. The finite element method (FEM) and component mode synthesis (CMS) method are used for low- to mid-frequency range computations and have been shown to yield reasonably close results. The statistical energy analysis (SEA) method yields predictions which agree with the CMS results for the 800 to 1000 Hz range provided that a sufficient number of modes participate. The CMS approach has been shown to yield valuable insight into the mid-frequency range of the analysis. It has been demonstrated that it is possible to conduct an analysis of a duct/hanger system in a cost-effective way for a wide frequency range, using several methods which overlap for several frequency bands.

Miller, R. D.↗
NASA NTRS1987

An approximate semi-analytical method for prediction of interlaminar shear stresses in an arbitrarily laminated thick plate

An approximate semianalytical method for determination of interlaminar shear stress distribution through the thickness of an arbitrarily laminated thick plate has been presented. The method is based on the assumptions of transverse inextensibility and layerwise constant shear angle theory (LCST) and utilizes an assumed quadratic displacement potential energy based finite element method (FEM). Centroid of the triangular surface has been proved from a rigorous mathematical point of view (Aubin-Nitsche theory), to be the point of exceptional accuracy for the interlaminar shear stresses. Numerical results indicate close agreement with the available three-dimensional elasticity theory solutions. A comparison between the present theory and that due to an assumed stress hybrid FEM suggest that the (normal) traction-free-edge condition is not satisfied in the latter approach. Furthermore, the present paper is the first to present the results for interlaminar shear stresses in a two-layer thick square plate of balanced unsymmetric angle-ply construction. A comparison with the recently proposed Equilibrium Method (EM) indicates the superiority of the present method, because the latter assures faster convergence as well as simultaneous vanishing of the transverse shear stresses on both of the exposed surfaces of the laminate. Superiority of the present method over the EM, in the case of a symmetric laminate, is limited to faster convergence alone. It has also been demonstrated that the combination of the present method and the reduced (quadratic order) numerical integration scheme yields convergence of the interlaminar shear stresses almost as rapidly as that of the nodal displacements, in the case of a thin plate.

Chaudhuri, Reaz A.↗
NASA NTRS1990

Dynamic substructuring by the boundary flexibility vector method of component mode synthesis

Component mode synthesis (CMS) is a method of dynamic analysis, for structures having a large number of degrees of freedom (DOF). These structures often required lengthy computer CPU time and large computer memory resources, if solved directly by the finite-element method (FEM). In CMS, the structure is divided into independent components in which the DOF are defined by a set of generalized coordinates defined by displacement shapes. The number of the generalized coordinates are much less than the original number of physical DOF, in the component. The displacement shapes are used to transform the component property matrices and any applied external loads, to a reduced system of coordinates. Reduced system property matrices are assembled, and any type of dynamic analysis is carried out in the reduced coordinate system. Any obtained results are back transformed to the original component coordinate systems. In all conventional methods of CMS, the mode shapes used for components are dynamic mode shapes, supplemented by static deflected shapes. Historically, all the dynamic mode shapes used in conventional CMS are the natural modes (eigenvectors) of components. A new method of CMS, namely the boundary flexibility vector method of CMS, is presented. The method provides for the incorporation of a set of static Ritz vectors, referred to as boundary flexibility vectors, as a replacement and/or supplement to conventional eigenvectors, as displacement shapes for components. The generation of these vectors does not require the solution of a costly eigenvalue problem, as in the case of natural modes in conventional CMS, and hence a substantial saving in CPU time can be achieved. The boundary flexibility vectors are generated from flexibility (or stiffness) properties of components. The formulation presented is for both free and fixed-interface components, and for both the free and forced vibration problems. Free and forced vibration numerical examples are presented to verify the accuracy of the method and the saving in CPU time. Compared to conventional methods of CMS, the results indicate that by using the new method, more accurate results can be obtained with a substantial saving in CPU time.

Abdallah, Ayman Ahmed↗
NASA NTRS1989

Elevated temperature crack growth

Alloy 718 crack growth experiments were conducted to assess the ability of the selected path-independent (P-I) integrals to describe the elevated temperature crack growth behavior. These tests were performed on single edge notch (SEN) specimens under displacement control with multiple extensometers to monitor the specimen and crack mouth opening displacement (CMOD). The displacements in these tests were sufficiently high to induce bulk cyclic inelastic deformation of the specimen. Under these conditions, the linear elastic fracture mechanics (LEFM) parameter K does not correlate the crack growth data. The experimentally measured displacement gradients at the end of specimen gage length were used as the boundary conditions in elastic-plastic finite element method (FEM) analyses. These analyses were performed with a node release approach using CYANIDE, a GEAE FEM code, which included a gap element which is capable of efficiently simulating crack closure. Excellent correlation was obtained between the experimentally measured and predicted variation of stress and CMOD with crack length and the stress-CMOD loops for Alloy 718 tests conducted at 538 C. This confirmed the accuracy of the FEM crack growth simulation approach. The experimentally measured crack growth rate data correlated well the selected P-I integrals. These investigations have produced significant progress in developing P-I integrals as non-linear fracture mechanics parameters. The results suggest that this methodology has the potential of accurately describing elevated temperature crack growth behavior under the combined influence of thermal cycling and bulk elastic-inelastic deformation states.

Kim, K. S.↗
NASA NTRS1989

Probabilistic boundary element method

The purpose of the Probabilistic Structural Analysis Method (PSAM) project is to develop structural analysis capabilities for the design analysis of advanced space propulsion system hardware. The boundary element method (BEM) is used as the basis of the Probabilistic Advanced Analysis Methods (PADAM) which is discussed. The probabilistic BEM code (PBEM) is used to obtain the structural response and sensitivity results to a set of random variables. As such, PBEM performs analogous to other structural analysis codes such as finite elements in the PSAM system. For linear problems, unlike the finite element method (FEM), the BEM governing equations are written at the boundary of the body only, thus, the method eliminates the need to model the volume of the body. However, for general body force problems, a direct condensation of the governing equations to the boundary of the body is not possible and therefore volume modeling is generally required.

Cruse, T. A.↗
NASA NTRS1990

Static aeroelastic analysis of composite wing

A static aeroelastic analysis capability that can predict aerodynamic loads for the deformed shape of the composite wing has been developed. The finite element method (FEM) was used for composite plate structural analysis, and the linear vortex lattice method (VLM) was used for steady aerodynamic analysis. The final deformed shape of the wing due to applied forces is determined by iterative manner using FEM and VLM. FEM and VLM analysis are related by a surface spline interpolation procedure. The wing with Gr/Ep composite material has been investigated to see the wing deformation effect. Aerodynamic load change due to wing flexibility has been investigated. Also, the effect of fiber orientation and sweep angle on the deformation pattern and aerodynamic coefficients are examined. For a certain fiber orientation, the deflection and aerodynamic loading of the composite wing is very much reduced. The swept forward wing has more significant effect of wing flexibility on aerodynamic coefficient than the swept back wing does.

Lee, IN↗
NASA NTRS1992

Validation of finite element and boundary element methods for predicting structural vibration and radiated noise

Analytical and experimental validation of methods to predict structural vibration and radiated noise are presented. A rectangular box excited by a mechanical shaker was used as a vibrating structure. Combined finite element method (FEM) and boundary element method (BEM) models of the apparatus were used to predict the noise radiated from the box. The FEM was used to predict the vibration, and the surface vibration was used as input to the BEM to predict the sound intensity and sound power. Vibration predicted by the FEM model was validated by experimental modal analysis. Noise predicted by the BEM was validated by sound intensity measurements. Three types of results are presented for the total radiated sound power: (1) sound power predicted by the BEM modeling using vibration data measured on the surface of the box; (2) sound power predicted by the FEM/BEM model; and (3) sound power measured by a sound intensity scan. The sound power predicted from the BEM model using measured vibration data yields an excellent prediction of radiated noise. The sound power predicted by the combined FEM/BEM model also gives a good prediction of radiated noise except for a shift of the natural frequencies that are due to limitations in the FEM model.

Seybert, A. F.↗
NASA NTRS1992

Response of launch pad structures to random acoustic excitation

Two solutions (probabilistic and deterministic) for the random vibration problem are presented in this paper from the standpoint of their applicability to predict the response of ground structures subjected to acoustic loading during the launch of a Space Shuttle. Deficiencies of the probabilistic method, especially to predict response in the low-frequency regime, prompted the development of the deterministic analysis, which offers a valid alternative. Challenges associated with the implementation of these response solutions in a commercially available Finite Element Method (FEM) code are briefly addressed.

Margasahayam, Ravi↗
NASA NTRS1993

Application of symbolic and algebraic manipulation software in solving applied mechanics problems

As its name implies, symbolic and algebraic manipulation is an operational tool which not only can retain symbols throughout computations but also can express results in terms of symbols. This report starts with a history of symbolic and algebraic manipulators and a review of the literatures. With the help of selected examples, the capabilities of symbolic and algebraic manipulators are demonstrated. These applications to problems of applied mechanics are then presented. They are the application of automatic formulation to applied mechanics problems, application to a materially nonlinear problem (rigid-plastic ring compression) by finite element method (FEM) and application to plate problems by FEM. The advantages and difficulties, contributions, education, and perspectives of symbolic and algebraic manipulation are discussed. It is well known that there exist some fundamental difficulties in symbolic and algebraic manipulation, such as internal swelling and mathematical limitation. A remedy for these difficulties is proposed, and the three applications mentioned are solved successfully. For example, the closed from solution of stiffness matrix of four-node isoparametrical quadrilateral element for 2-D elasticity problem was not available before. Due to the work presented, the automatic construction of it becomes feasible. In addition, a new advantage of the application of symbolic and algebraic manipulation found is believed to be crucial in improving the efficiency of program execution in the future. This will substantially shorten the response time of a system. It is very significant for certain systems, such as missile and high speed aircraft systems, in which time plays an important role.

Tsai, Wen-Lang↗