SEARCH · Search NASA
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.
An incremental, nonlinear, finite element program capable of tracing the damage that occurs in a notched laminated composite, before ultimate failure takes place, is described. Computer predictions for the ultimate failure load and failure modes are compared to experimental results for two laminates containing holes and loaded in uniaxial tension. The comparison is quite favorable.
A finite element hydroelastic analysis formulation is developed on the basis of Toupin's complementary variational principle. Emphasis is placed on the special case of an incompressible fluid model which is applicable to propellant tank hydroelastic analysis. A concise fluid inertia representation results from the assumption of incompressibility and the hydroelastic equations reduce to a simplified form associated with the structure alone. The efficiency of the incompressible hydroelastic formulation in unhanced for both fluid and structure by introduction of harmonic reduction as an alternative to Guyan reduction. The theoretical developments are implemented in the NASTRAN Program and the technique is verified and demonstrated as an efficient and accurate approach with a series of illustrative problems including the 1/8 scale space shuttle external tank.
A finite-element model of a 500-km diameter lunar basin flooded with 10 km of basalt was used to represent the central portion of a multi-ring basin undergoing elastic deformation as a result of basin subsidence. Gravitational forces due to the weight of the basalt were set as the initial forces in the model; elastic constants computed from seismic velocities of lunar samples provided strength specifications. In addition, the crust was constrained by permitting no displacements at a depth of 60 km. Because of the strength of the lithosphere, surface subsidence amounted to only 20 m. This magnitude of elastic deformation does not equal that found on the lunar surface, but the orientation of principal stress and the position of circumferential compressional stress are compatible with the evolution of inner-ring and radial ridge systems by gravitational-induced downdropping.
The finite element method for fluid dynamics was used to develop a two dimensional mesh generation scheme. The method consists of shearing and conformal maps with upper and lower surfaces handled independently to allow sharp leading edges. The method also generates meshes of triangular or quadrilateral elements.
A finite-element method for the solution of nonlinear hyperbolic systems of equations, such as those encountered in non-self-adjoint problems of transient phenomena in convection-diffusion or in the mixed representation of wave problems, is developed and demonstrated. The problem is rewritten in moving coordinates and reinterpolated to the original mesh by a Taylor expansion prior to a standard Galerkin spatial discretization, and it is shown that this procedure is equivalent to the time-discretization approach of Donea (1984). Numerical results for sample problems are presented graphically, including such shallow-water problems as the breaking of a dam, the shoaling of a wave, and the outflow of a river; compressible flows such as the isothermal flow in a nozzle and the Riemann shock-tube problem; and the two-dimensional scalar-advection, nonlinear-shallow-water, and Euler equations.
A finite element remeshing approach that makes use of quadrilateral and triangular elements is described. The approach uses the solution on a previous mesh to create a new mesh. Meshes are completely unstructured with highly refined elements in regions of steep gradients and larger elements where gradients are smaller. Studies of convergence rates for heat conduction problems with exact solutions show that for problems with highly localized solution variations, the remeshing approach gives smaller solution errors with fewer unknowns than refinement of uniform, structured meshes.
A finite-element formulation is developed for the free vibration of rectangular plates which are under the influence of moderately large stress-free initial deflections and large thermal deflections. The von Karman nonlinear strain-displacement relations are used to account for the thermal deflections. The plates are thin, isotropic, and Hookean in nature. The temperature imposed on the plate is assumed to be constant through the thickness of the plate. Uniform and sinusoidal temperature distributions are studied. The material properties of the plates are temperature-dependent due to the relatively high temperatures imposed on the plates.
A finite-element model of the Navier-Stokes equations is used for numerical simulation of flow past two normal flat plates arranged side by side at Reynolds number 80 and 160. The results from this simulation indicate that when the gap between the plates is twice the width of a single plate, the individual wakes of the plates behave independently, with the antiphase vortex shedding being dominant. At smaller gap sizes, the in-phase vortex shedding, with strong wake interaction, is favored. The gap flow in those cases becomes biased, with one of the wakes engulfing the other. The direction of the biased flow was found to be switching at irregular intervals, with the time histories of the indicative flow parameters and their power spectra resembling those of a chaotic system.
A 3D finite element (FE) approach was developed and implemented for computation of global magnetic fields in a 14.3 kVA modified Lundell alternator. The essence of the new method is the combined use of magnetic vector and scalar potential formulations in 3D FEs. This approach makes it practical, using state of the art supercomputer resources, to globally analyze magnetic fields and operating performances of rotating machines which have truly 3D magnetic flux patterns. The 3D FE-computed fields and machine inductances as well as various machine performance simulations of the 14.3 kVA machine are presented in this paper and its two companion papers.
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.
Computer-oriented finite-element techniques for extracting the required information about the influence of cracks on stress fields are reviewed based on the knowledge that cracks induce stress singularities. Attention is focused on modeling behavior at a crack tip. Also discussed are crack tip analysis, utility of crack tip finite elements, and newer developments for cases where the singularity is not known a priori, i.e. the strength of the singularity is not established in advance of the computation. Available results are briefly summarized.
Three-dimensional finite-element analyses are used to assess the accuracy of simplified composite micromechanics equations (SME) for hygral, thermal, and mechanical properties of unidirectional composites with orthotropic fibers. The properties predicted by the SME are in reasonably good agreement with those predicted by the three-dimensional finite-element analyses. This correlation demonstrates that the SME can be used with confidence in predicting the hygral, thermal, and mechanical behavior of unidirectional fiber composites.
An effective approach in the finite element analysis of the stress field at the traction free boundary of a solid continuum was studied. Conventional displacement and assumed stress finite elements were used in the determination of stress concentrations around circular and elliptical holes. Specialized hybrid elements were then developed to improve the satisfaction of prescribed traction boundary conditions. Results of the stress analysis indicated that finite elements which exactly satisfy the free stress boundary conditions are the most accurate and efficient in such problems. A general approach for hybrid finite elements which incorporate traction free boundaries of arbitrary geometry was formulated.
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 virtual 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 finte-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.
Finite element substructuring strategies are presented to obtain numerical solutions for three typical problems of interest to the composites community: (1) impact and toughness characterization of composites using Charpy's impact test specimen; (2) free-edge stress analysis of composite laminates; and (3) fracture toughness predictions of composites for individual and combined fracture of modes I, II, and III. The key issue common to these problems is the presence of singular or near singular stress fields. The regions prone to see steep stress gradients are substructured with progressively refined meshes to study the local response simultaneously with the global response. The results from the select examples indicate that finite element substructuring methods are computationally effective for composite singularity mechanics.
Finite element substructuring strategies are presented to obtain numerical solutions for three typical problems of interest to the composites community: (1) impact and toughness characterization of composites using Charpy's impact test specimen; (2) free-edge stress analysis of composite laminates; and (3) fracture toughness predictions of composites for individual and combined fracture of modes I, II, and III. The key issue common to these problems is the presence of singular or near singular stress fields. The regions prone to see steep stress gradients are substructured with progressively refined meshes to study the local response simultaneously with the global response. The results from the select examples indicate that finite element substructuring methods are computationally effective for composite singularity mechanics.
Electromagnetic finite elements based on a variational principle that uses the electromagnetic four-potential as primary variable are derived. This choice is used to construct elements suitable for downstream coupling with mechanical and thermal finite elements for the analysis of electromagnetic/mechanical systems that involve superconductors. The main advantages of the four-potential as a basis for finite element formulation are: (1) the number of degrees of freedom per node remains modest as the problem dimensionality increases, (2) jump discontinuities on interfaces are naturally accommodated, and (3) statics as well as dynamics may be treated without any a priori approximations. The new elements are tested on an axisymmetric problem under steady-state forcing conditions. The results are in excellent agreement with analytical solutions.
Electromagnetic finite elements based on a variational principle that uses the electromagnetic four-potential as a primary variable are derived. This choice is used to construct elements suitable for downstream coupling with mechanical and thermal finite elements for the analysis of electromagnetic/mechanical systems that involve superconductors. The main advantages of the four-potential as a basis for finite element formulation are that the number of degrees of freedom per node remains modest as the problem dimensionally increases, that jump discontinuities on interfaces are naturally accommodated, and that statics as well as dynamics may be treated without any a priori approximations. The new elements are tested on an axisymmetric problem under steady state forcing conditions. The results are in excellent agreement with analytical solutions.