Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite element analysis”

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 379 records · Page 21

Parallel processing in finite element structural analysis

A brief review is made of the fundamental concepts and basic issues of parallel processing. Discussion focuses on parallel numerical algorithms, performance evaluation of machines and algorithms, and parallelism in finite element computations. A computational strategy is proposed for maximizing the degree of parallelism at different levels of the finite element analysis process including: 1) formulation level (through the use of mixed finite element models); 2) analysis level (through additive decomposition of the different arrays in the governing equations into the contributions to a symmetrized response plus correction terms); 3) numerical algorithm level (through the use of operator splitting techniques and application of iterative processes); and 4) implementation level (through the effective combination of vectorization, multitasking and microtasking, whenever available).

Noor, Ahmed K.↗

Interactive modeling using the PATRAN-G program

PATRAN-G is an interactive finite element generating and graphics display program. The program is currently executable on the VAX or PRIME computers and is compatible with 19 commercially available graphics terminals. The PATRAN program can be used in two modes: generation of a finite element model, or postprocessing of output data from a finite element analysis program. The finite element generating mode is divided into two phases. Phase I is the geometry construction. During this phase, the geometric boundaries of the model are defined and the material properties are input. Mirror imaging, translation and scaling of lines, surfaces, and volumes, and rotation of regions permit the geometry to be described efficiently. Points may be defined in alternate rectangular, cylindrical, or spherical coordinate systems. Also, PATRAN-G can calculate the intersection of lines or surfaces. This capability allows a user to easily model, for example, the intersection of a cylinder which cuts through a sphere at an angle. Once the geometry of the model has been described in terms of points, lines, surfaces, and volumes. Phase II is used to generate the computer finite element model.

Stalnaker, W. A.↗

Calculation of skin-stiffener interface stresses in stiffened composite panels

A method for computing the skin-stiffener interface stresses in stiffened composite panels is developed. Both geometrically linear and nonlinear analyses are considered. Particular attention is given to the flange termination region where stresses are expected to exhibit unbounded characteristics. The method is based on a finite-element analysis and an elasticity solution. The finite-element analysis is standard, while the elasticity solution is based on an eigenvalue expansion of the stress functions. The eigenvalue expansion is assumed to be valid in the local flange termination region and is coupled with the finite-element analysis using collocation of stresses on the local region boundaries. Accuracy and convergence of the local elasticity solution are assessed using a geometrically linear analysis. Using this analysis procedure, the influence of geometric nonlinearities and stiffener parameters on the skin-stiffener interface stresses is evaluated.

Cohen, David↗

Finite element stress analysis of idealized composite damage zones

A quasi three dimensional finite element stress analysis of idealized damage zones in composite laminates is presented. The damage zones consist of a long centered groove or cutout extending one or two layers in depth from both top and bottom surfaces of a thin composite laminate. Elastic results are presented for compressive loading of four and eight layer laminates. It is shown that a boundary layer exists near the cutout edge similar to that previously shown to exist along free edges. The cutout is shown to produce significant interlaminar stresses in the interior of the laminate away from free cutout edges. The interlaminar stresses are also shown to contribute to failure which is defined using the Tsai-Wu failure criteria.

Obrien, D.↗

A finite element thermal analysis procedure for several temperature-dependent parameters

A finite-element thermal analysis procedure for elements with several temperature-dependent thermal parameters is presented. The procedure, based on an application of the Newton-Raphson iteration technique, is formulated by resolving element matrices into component matrices, one component for each thermal parameter. Component conductance matrices are evaluated by assuming constant thermal parameters within an element and are computed once per unit thermal parameter. Significant savings in computer time result from the unit thermal parameter concept. The solution procedure applied to a convectively cooled structure with significantly varying thermal parameters converged in four iterations.

Thornton, E. A.↗

Efficient 3-D finite element failure analysis of compression loaded angle-ply plates with holes

Finite element stress analysis and the tensor polynomial failure criterion predict that failure always initiates at the interface between layers on the hole edge for notched angle-ply laminates loaded in compression. The angular location of initial failure is a function of the fiber orientation in the laminate. The dominant stress components initiating failure are shear. It is shown that approximate symmetry can be used to reduce the computer resources required for the case of unaxial loading.

Burns, S. W.↗

Delaminations in composite plates under impact loads

A method is presented for calculating the locations, shapes, and sizes of delaminations which occur in a fiber reinforced composite plate subjected to non-penetrating (low velocity) impact of a solid object. The plate may be simply supported, clamped, or free along its edges. A failure model of the delamination formation was developed. This model was then coupled with a finite element analysis. The model and the finite element analysis were then implemented by a computer code (IMPACT-ST) which can be used to estimate the damage initiation load and the locations, shapes, and sizes of the delaminations. Tests were performed measuring the geometries of the delaminations in graphite-epoxy, graphite-toughened epoxy, and graphite-PEEK plates impacted by a projectile with a spherical tip having masses ranging from 0.355 lbm to 0.963 lbm and velocities from 50 in/sec to 225 in/sec. The data were compared to the results of the model, and good agreements were found between the measured and the calculated delamination lengths and widths.

Finn, Scott R.↗

Delaminations in composite plates under transverse static or impact loads - A model

A method is presented for calculating the locations, shapes, and sizes of delaminations which occur in a fiber reinforced composite plate subjected to transverse static or dynamic (impact) loads. The plate may be simply supported, clamped, or free along its edges. A model of the delamination formation was developed. This model was then coupled with a finite element analysis. The model and the finite element analysis were implemented by a computer code which can be used to estimate the load at which damage is initiated as well as the locations, shapes, and sizes of the delaminations.

Finn, Scott R.↗

Finite element dynamic analysis on CDC STAR-100 computer

Computational algorithms are presented for the finite element dynamic analysis of structures on the CDC STAR-100 computer. The spatial behavior is described using higher-order finite elements. The temporal behavior is approximated by using either the central difference explicit scheme or Newmark's implicit scheme. In each case the analysis is broken up into a number of basic macro-operations. Discussion is focused on the organization of the computation and the mode of storage of different arrays to take advantage of the STAR pipeline capability. The potential of the proposed algorithms is discussed and CPU times are given for performing the different macro-operations for a shell modeled by higher order composite shallow shell elements having 80 degrees of freedom.

Noor, A. K.↗

Three-Dimensional Magnetic Analysis Technique Developed for Evaluating Stirling Convertor Linear Alternators

The Department of Energy, the Stirling Technology Company (STC), and the NASA Glenn Research Center are developing Stirling convertors for Stirling radioisotope generators to provide electrical power for future NASA deep space missions. STC is developing the 55-We technology demonstration convertor (TDC) under contract to the Department of Energy. The Department of Energy recently named Lockheed Martin as the system integration contractor for the Stirling radioisotope generator development project. Lockheed Martin will develop the Stirling radioisotope generator engineering unit and has contract options to develop the qualification unit and the first flight unit. Glenn s role includes an in-house project to provide convertor, component, and materials testing and evaluation in support of the overall power system development. As a part of this work, Glenn has established an in-house Stirling research laboratory for testing, analyzing, and evaluating Stirling machines. STC has built four 55-We convertors for NASA, and these are being tested at Glenn. A cross-sectional view of the 55-We TDC is shown in the figure. Of critical importance to the successful development of the Stirling convertor for space power applications is the development of a lightweight and highly efficient linear alternator. In support, Glenn has been developing finite element analysis and finite element method tools for performing various linear alternator thermal and electromagnetic analyses and evaluating design configurations. A three-dimensional magnetostatic finite element model of STC's 55-We TDC linear alternator was developed to evaluate the demagnetization fields affecting the alternator magnets. Since the actual linear alternator hardware is symmetric to the quarter section about the axis of motion, only a quarter section of the alternator was modeled. The components modeled included the mover laminations, the neodymium-iron-boron magnets, the stator laminations, and the copper coils. The three-dimensional magnetostatic model was then coupled with a circuit simulator model of the alternator load and convertor controller. The coupled model was then used to generate alternator terminal voltage and current predictions. The predicted voltage and current waveforms agreed well with the experimental data, which tended to validate the accuracy of the coupled model. The model was then used to generate predictions of the demagnetization fields acting on the alternator magnets for the alternator under load. The preliminary model predictions indicate that the highest potential for demagnetization is along the inside surface of the uncovered magnets. The demagnetization field for the uncovered magnets when the mover is positioned at the end of a stroke is higher than it is when the mover is at the position of maximum induced voltage or maximum alternator current. Assuming normal load conditions, the model predicted that the onset of demagnetization is most likely to occur for magnet temperatures above 101 C.

Geng, Steven M.↗

Finite element stress analysis of a notched coupon specimen for in-plane shear behaviour of composites

The results of a linear elastic, plane stress finite element investigation of the stress distribution in a double V-notched coupon specimen are presented for an isotropic material (steel) and five graphite/polyimide laminates: (0), (90), (0/90)s, (plus or minus 45)s and (0/90/plus or minus 45)s. Stress contours, stress profiles, and stress concentration factors are presented for specimens loaded through rigid fixtures. It is shown that the finite element analysis predicts a region of essentially uniform pure shear in the central portion of the specimen for all laminates considered. Stress concentrations vary considerably with laminate configuration. The influences of elastic fixtures and thermal stresses on the stress distribution in (0) laminates are also investigated. The (0) laminate is predicted to exhibit maximum pure shear in the center of the test section for both rigid and elastic fixtures. It is concluded that the specimen is a practicable candidate for use as a composite shear specimen.

Herakovich, C. T.↗

Study on prestressed concrete reactor vessel structures. II-5: Crack analysis by three dimensional finite elements method of 1/20 multicavity type PCRV subjected to internal pressure

A three-dimensional finite elements analysis is reported of the nonlinear behavior of PCRV subjected to internal pressure by comparing calculated results with test results. As the first stage, an analysis considering the nonlinearity of cracking in concrete was attempted. As a result, it is found possible to make an analysis up to three times the design pressure (50 kg/sqcm), and calculated results agree well with test results.

Source record↗

Strain-energy-release rate analysis of the end-notched flexure specimen using the finite-element method

Two-dimensional finite-element analysis of the end-notched flexure specimen was performed using 8-node isoparametric, parabolic elements to evaluate compliance and mode II strain energy release rates, G sub II. The G sub II values were computed using two different techniques: the virtural crack-closure technique (VCCT) and the rate of change of compliance with crack length (compliance derivative method). The analysis was performed for various crack-length-to-semi-span (a/L) ratios ranging from 0.2 to 0.9. Three material systems representing a wide range of material properties were analyzed. The compliance and strain energy release rates of the specimen calculated with the present finite-element analysis agree very well with beam theory equations including transverse shear. The G sub II values calculated using the compliance derivative method compared extremely well with those calculated using the VCCT. The G sub II values obtained by the compliance derivative method using the top or bottom beam deflections agreed closely with each other. The strain energy release rates from a plane-stress analysis were higher than the plane-strain values by only a small percentage, indicating that either assumption may be used in the analysis. The G sub II values for one material system calculated from the finite-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.

Salpekar, S. A.↗

A multi-microprocessor system for finite element structural analysis

During the last few years, advances in microprocessor technology have spurred a renewed interest in special-purpose computers. The microprocessor has become small, inexpensive, and powerful enough to be considered as a building block for special-purpose hardware. A description is presented of the architecture of a prototype 'finite element machine' currently being built. Attention is given to details regarding the finite element analysis problem, the arrangement of the processors as finite element nodes in the structural model, the influence of the architecture on the solution algorithm, interprocessor communication primitives, and the performance of the finite element machine.

Jordan, H. F.↗

Elimination of Gibbs' phenomena from error analysis of finite element results

This paper is one of a series on error analysis and correction of finite element solutions for plates and shells. The error analysis in the earlier papers used half-range double Fourier sine series for numerical harmonic analysis. The half-range formulas are simple to apply, but they can be inaccurate near the ends of the ranges of the independent variables. The Gibbs' phenomenon exhibited by half-range sine series in one independent variable has a two-dimensional analog; a classic example is the Navier solution in a double half-range sine series for the simply supported plate under a uniform load. A simple change of variables is introduced in the paper to improve the accuracy of the double Fourier sine series without adding complexity to the numerical analysis. The change of variables is applied to the problem of approximating a transverse load that is tabulated on a rectangular grid. A solution based on the change of variables is compared with results from the Navier solution for the simply supported plate problem and finite element results for the same problem.

Thurston, Gaylen A.↗

Wavelet and Multiresolution Analysis for Finite Element Networking Paradigms

This paper presents a final report on Wavelet and Multiresolution Analysis for Finite Element Networking Paradigms. The focus of this research is to derive and implement: 1) Wavelet based methodologies for the compression, transmission, decoding, and visualization of three dimensional finite element geometry and simulation data in a network environment; 2) methodologies for interactive algorithm monitoring and tracking in computational mechanics; and 3) Methodologies for interactive algorithm steering for the acceleration of large scale finite element simulations. Also included in this report are appendices describing the derivation of wavelet based Particle Image Velocity algorithms and reduced order input-output models for nonlinear systems by utilizing wavelet approximations.

Kurdila, Andrew J.↗

Simplified analytical procedures for representing material cyclic response

Requirements for increased durability of gas turbine hot section structural components have made it necessary to place greater emphasis on accurate structural analysis and life prediction. Linear finite-element analysis is generally sufficient for structural analysis applications. However, for structures in the hot part of the engine, nonlinear structural analysis may be required under certain conditions for the accurate prediction of the local stress-strain response. Nonlinear finite element analysis represents a costly effort which is generally incompatible with the iterative nature of the design process. The present investigation is, therefore, concerned with two simplified procedures for estimating the local hysteretic response produced by cyclic thermal loading. These procedures reduce the need for nonlinear finite-element analysis.

Moreno, V.↗