Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 199 records · Page 11

Finite Element Analysis of Geodesically Stiffened Cylindrical Composite Shells Using a Layerwise Theory

Layerwise finite element analyses of geodesically stiffened cylindrical shells are presented. The layerwise laminate theory of Reddy (LWTR) is developed and adapted to circular cylindrical shells. The Ritz variational method is used to develop an analytical approach for studying the buckling of simply supported geodesically stiffened shells with discrete stiffeners. This method utilizes a Lagrange multiplier technique to attach the stiffeners to the shell. The development of the layerwise shells couples a one-dimensional finite element through the thickness with a Navier solution that satisfies the boundary conditions. The buckling results from the Ritz discrete analytical method are compared with smeared buckling results and with NASA Testbed finite element results. The development of layerwise shell and beam finite elements is presented and these elements are used to perform the displacement field, stress, and first-ply failure analyses. The layerwise shell elements are used to model the shell skin and the layerwise beam elements are used to model the stiffeners. This arrangement allows the beam stiffeners to be assembled directly into the global stiffness matrix. A series of analytical studies are made to compare the response of geodesically stiffened shells as a function of loading, shell geometry, shell radii, shell laminate thickness, stiffener height, and geometric nonlinearity. Comparisons of the structural response of geodesically stiffened shells, axial and ring stiffened shells, and unstiffened shells are provided. In addition, interlaminar stress results near the stiffener intersection are presented. First-ply failure analyses for geodesically stiffened shells utilizing the Tsai-Wu failure criterion are presented for a few selected cases.

Gerhard, Craig Steven↗

Global-Local Finite Element Analysis for Thermo-Mechanical Stresses in Bonded Joints

An analysis of adhesively bonded joints using conventional finite elements does not capture the singular behavior of the stress field in regions where two or three dissimilar materials form a junction with or without free edges. However, these regions are characteristic of the bonded joints and are prone to failure initiation. This study presents a method to capture the singular stress field arising from the geometric and material discontinuities in bonded composites. It is achieved by coupling the local (conventional) elements with global (special) elements whose interpolation functions are constructed from the asymptotic solution.

Shkarayev, S.↗

Computerized symbolic manipulation in nonlinear finite element analysis

The potential of using computerized symbolic manipulation in the development of nonlinear finite elements is discussed. Three tasks which can be efficiently performed using computerized symbolic manipulation are identified: (1) generation of algebraic expressions for the stiffness coefficients of nonlinear finite elements, (2) generation of FORTRAN source code for numerical evaluation of stiffness coefficients, and (3) checking the correctness of the FORTRAN statements for the arrays of coefficients. The symbolic and algebraic manipulation system MACSYMA is used in the present study. Two sample MACSYMA programs are presented for the development of the nonlinear stiffness coefficients of two-dimensional, shear-flexible, doubly-curved deep shell elements. The first program is for displacement models and the second program is for mixed models with discontinuous stress-resultant fields at interelement boundaries.

Noor, A. K.↗

Parallel processing techniques for finite element analysis of nonlinear large truss structures

Methods were developed for parallel processing of finite element solutions of large truss structures. The parallel processing techniques were implemented in two stages, i.e., the repeated forming of the nonlinear global stiffness matrix and the solving of the global system of equations. The Sequent Balance 21000 parallel computer was employed to demonstrate the procedures and the speed-up.

Chien, L. S.↗

Finite Element Analysis of a NASA National Transonic Facility Wind Tunnel Balance

This paper presents the results of finite element analyses and correlation studies performed on a NASA National Transonic Facility (NTF) Wind Tunnel balance. In the past NASA has relied primarily on classical hand analyses, coupled with relatively large safety factors, for predicting maximum stresses in wind tunnel balances. Now, with the significant advancements in computer technology and sophistication of general purpose analysis codes, it is more reasonable to pursue finite element analyses of these balances. The correlation studies of the present analyses show very good agreement between the analyses and data measured with strain gages and therefore the studies give higher confidence for using finite element analyses to analyze and optimize balance designs in the future.

Lindell, Michael C.↗

Finite Element Analysis of a NASA National Transonic Facility Wide Tunnel Balance

This paper presents the results of finite element analyses and correlation studies performed on a NASA National Transonic Facility (NTF) Wind Tunnel balance. In the past NASA has relied primarily on classical hand analyses, coupled with relatively large safety factors, for predicting maximum stresses in wind tunnel balances. Now, with the significant advancements in computer technology and sophistication of general purpose analysis codes, it is more reasonable to pursue finite element analyses of these balances. The correlation studies of the present analyses show very good agreement between the analyses and data measured with strain gages and therefore the studies give higher confidence for using finite element analyses to analyze and optimize balance designs in the future.

Lindell, Michael C.↗

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.↗

Finite element analysis of instability related delamination growth

Instability-related growth of through-width delaminations in laminated coupons was investigated analytically and experimentally. In the analytical effort, a geometrically nonlinear finite-element stress analysis was used, and the numerical results were compared with exact solutions. In addition, measured lateral deflections of postbuckled through-width delaminations were compared with predicted deflections. Lateral deflections, stress distributions, and strain energy release rates were calculated for various delamination lengths, delamination depths, and applied loads. Delamination growth rates observed for unidirectional graphite/epoxy bonded to aluminum were qualitatively correlated with calculated strain-energy release rates and it was found that the Mode I component of energy release dominates the growth process.

Whitcomb, J. D.↗

A Procedure for Modeling Structural Component/Attachment Failure Using Transient Finite Element Analysis

Structures often comprise smaller substructures that are connected to each other or attached to the ground by a set of finite connections. Under static loading one or more of these connections may exceed allowable limits and be deemed to fail. Of particular interest is the structural response when a connection is severed (failed) while the structure is under static load. A transient failure analysis procedure was developed by which it is possible to examine the dynamic effects that result from introducing a discrete failure while a structure is under static load. The failure is introduced by replacing a connection load history by a time-dependent load set that removes the connection load at the time of failure. The subsequent transient response is examined to determine the importance of the dynamic effects by comparing the structural response with the appropriate allowables. Additionally, this procedure utilizes a standard finite element transient analysis that is readily available in most commercial software, permitting the study of dynamic failures without the need to purchase software specifically for this purpose. The procedure is developed and explained, demonstrated on a simple cantilever box example, and finally demonstrated on a real-world example, the American Airlines Flight 587 (AA587) vertical tail plane (VTP).

Lovejoy, Andrew E.↗

Global functions in global-local finite-element analysis of localized stresses in prismatic structures

An important consideration in the global local finite-element method (GLFEM) is the availability of global functions for the given problem. The role and mathematical requirements of these global functions in a GLFEM analysis of localized stress states in prismatic structures are discussed. A method is described for determining these global functions. Underlying this method are theorems due to Toupin and Knowles on strain energy decay rates, which are related to a quantitative expression of Saint-Venant's principle. It is mentioned that a mathematically complete set of global functions can be generated, so that any arbitrary interface condition between the finite element and global subregions can be represented. Convergence to the true behavior can be achieved with increasing global functions and finite-element degrees of freedom. Specific attention is devoted to mathematically two-dimensional and three-dimensional prismatic structures. Comments are offered on the GLFEM analysis of NASA flat panel with a discontinuous stiffener. Methods for determining global functions for other effects are also indicated, such as steady-state dynamics and bodies under initial stress.

Dong, Stanley B.↗

A fluid--structure interaction finite element analysis of pulsatile blood flow through a compliant stenotic artery

A new model is used to analyze the fully coupled problem of pulsatile blood flow through a compliant, axisymmetric stenotic artery using the finite element method. The model uses large displacement and large strain theory for the solid, and the full Navier-Stokes equations for the fluid. The effect of increasing area reduction on fluid dynamic and structural stresses is presented. Results show that pressure drop, peak wall shear stress, and maximum principal stress in the lesion all increase dramatically as the area reduction in the stenosis is increased from 51 to 89 percent. Further reductions in stenosis cross-sectional area, however, produce relatively little additional change in these parameters due to a concomitant reduction in flow rate caused by the losses in the constriction. Inner wall hoop stretch amplitude just distal to the stenosis also increases with increasing stenosis severity, as downstream pressures are reduced to a physiological minimum. The contraction of the artery distal to the stenosis generates a significant compressive stress on the downstream shoulder of the lesion. Dynamic narrowing of the stenosis is also seen, further augmenting area constriction at times of peak flow. Pressure drop results are found to compare well to an experimentally based theoretical curve, despite the assumption of laminar flow.

Non-NASA Center↗

ZIP3D: An elastic and elastic-plastic finite-element analysis program for cracked bodies

ZIP3D is an elastic and an elastic-plastic finite element program to analyze cracks in three dimensional solids. The program may also be used to analyze uncracked bodies or multi-body problems involving contacting surfaces. For crack problems, the program has several unique features including the calculation of mixed-mode strain energy release rates using the three dimensional virtual crack closure technique, the calculation of the J integral using the equivalent domain integral method, the capability to extend the crack front under monotonic or cyclic loading, and the capability to close or open the crack surfaces during cyclic loading. The theories behind the various aspects of the program are explained briefly. Line-by-line data preparation is presented. Input data and results for an elastic analysis of a surface crack in a plate and for an elastic-plastic analysis of a single-edge-crack-tension specimen are also presented.

Shivakumar, K. N.↗

Efficient algorithms for use in probabilistic finite element analysis

This paper investigates the use of Fast Probability Integration (FPI) algorithms in a Finite Element environment. A method allowing the representation of correlated fields in terms of a vector of uncorrelated transformed variables, based on the spectral decomposition of the variance-covariance matrix is developed. The response of the deterministic model corresponding to selected perturbations of these uncorrelated variables is then obtained via a Newton-type iterative scheme. The results of the perturbed problems are used to construct a local representation of the model's behavior in the neighborhood of the deterministic state, which the FPI algorithm will use to estimate the reliability of the system. Although the proposed strategy has thus far only been applied to linear elastostatics, the extension of the method to a broader class of problems appears to be feasible.

Dias, J. B.↗

A finite element analysis of the freeze/thaw behavior of external artery heat pipes

A two-dimensional finite element model was used to determine the freeze/thaw characteristics of an external artery heat pipe. During startup, the working fluid, which was located in the liquid channel and the circumferential wall grooves, experienced a phase transformation from a solid to a liquid state. The transient heat conduction equations with moving interfacial conditions were solved using the appropriate initial boundary conditions. The modelling results include the cross-sectional temperature distribution and the interfacial or melt front position as a function of time. A fixed grid approach was adopted in the model for the phase-change process during thawing of frozen working fluid. The interfacial position between the liquid and solid regions was found by balancing the latent heat caused by interfacial movement with the heat addition or extraction at the related grid points.

Lu, X. J.↗

Finite-element-analysis model and preliminary ground testing of controls-structures interaction evolutionary model reflector

Results of two different nonlinear finite element analyses and preliminary test results for the final design of the Controls-Structures Interaction Evolutionary Model are presented. Load-deflection data bases are generalized from analysis and testing of the 16-foot diameter, dish shaped reflector. Natural frequencies and mode shapes are obtained from vibrational analysis. Experimental and analytical results show similar trends; however, future test hardware modifications and finite element model refinement would be necessary to obtain better correlation. The two nonlinear analysis procedures are both adequate techniques for the analysis of prestressed structures with complex geometries.

Reaves, Mercedes C.↗

Validation of finite element analysis strategy to investigate acoustic levitation in a two-axis acoustic levitator

A two-axis acoustic levitator can be used to generate a standing pressure wave capable of levitating solid and liquid particles at appropriate input conditions. This work proposes a simulation framework to investigate the two-axis levitation particle stability using a commercial, computational fluid dynamics software based on the harmonic solution to the acoustic wave equation. The simulation produced predictions of the standing wave that include a strong "+" shaped pattern of nodes and anti-nodes that are aligned with the levitator axes. To verify the simulation, a levitator was built and used to generate the standing wave. The field was probed with a microphone and a motorized-scanning system. After scaling the simulated pressure to the measured pressure, the magnitudes of the sound pressure level at corresponding high-pressure locations were different by no more than 5%. This is the first time a measurement of a two-axis levitator standing pressure wave has been presented and shown to verify simulations. Furthermore, as an additional verification, the authors consulted high speed camera measurements of a reference-levitator transducer, which was found to have a maximum peak-to-peak displacement of 50 +/- 5 mu m. The reference-levitator is known to levitate water at 160 dB. The system for this work was simulated to match the operation of the reference-levitator so that it produced sound pressure levels of 160 dB. This pressure was achieved when the transducer maximum peak-to-peak displacement was 50.8 mu m. The agreement between the two levitators' displacements provides good justification that the modeling approach presented here produces reliable results.

2 axis levitator↗

Finite element analysis of flexible, rotating blades

A reference guide that can be used when using the finite element method to approximate the static and dynamic behavior of flexible, rotating blades is given. Important parameters such as twist, sweep, camber, co-planar shell elements, centrifugal loads, and inertia properties are studied. Comparisons are made between NASTRAN elements through published benchmark tests. The main purpose is to summarize blade modeling strategies and to document capabilities and limitations (for flexible, rotating blades) of various NASTRAN elements.

Mcgee, Oliver G.↗