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 181 records · Page 10

Time-domain Green's Function Method for three-dimensional nonlinear subsonic flows

The Green's Function Method for linearized 3D unsteady potential flow (embedded in the computer code SOUSSA P) is extended to include the time-domain analysis as well as the nonlinear term retained in the transonic small disturbance equation. The differential-delay equations in time, as obtained by applying the Green's Function Method (in a generalized sense) and the finite-element technique to the transonic equation, are solved directly in the time domain. Comparisons are made with both linearized frequency-domain calculations and existing nonlinear results.

Tseng, K.↗

Biaxial load effects in fracture mechanics

It is found that the standard expressions for elastic stress and displacement in the crack-tip region (i.e., the so-called singular solution) cannot be considered to be approximations that are acceptable in a completely general sense. This conclusion is best illustrated by the instance of a biaxially loaded infinite sheet with a flat horizontal central crack, where the effect of load applied parallel to the plane of the crack appears entirely in the second terms of the series representations for local stresses and displacements. An elastoplastic finite-element analysis of the same biaxially loaded finite specimen geometry shows that the global energy release rate, the J-integral, the plastic stress and strain intensity factors (in the sense of Hilton and Hutchinson), and the size of the crack border region plastic yield, all have pronounced biaxial load dependence.

Liebowitz, H.↗

Computerized structural sizing at NASA Langley Research Center

Programs at the NASA Langley Research Center associated with the development of computerized structural sizing technology are reviewed. Particular attention is given to (1) lightweight columns for space structure applications, (2) stiffened composite panels for aerospace structures, (3) thermal structures for high-speed aircraft and space vehicles, (4) structural sizing methodology for finite-element structural models, (5) the sizing of large complex structural systems in multidisciplinary environments. Improvements to computational efficiency are noted with reference to a reduced number of sizing variables, a reduced number of constraints, and improved sizing algorithms.

Stroud, W. J.↗

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

Dynamic simulation and stability analysis

Dynamic Interaction Simulation of Controls and Structure (DISCOS) program was developed for dynamic simulation and stability analysis of passive and actively controlled spacecraft. Program is written in FORTRAN IV for batch execution and requires access to finite-element structures program as NASTRAN for flexible-body input data.

Frisch, H. P.↗

Analysis of fatigue damage in composites

Finite-element heat-transfer analysis determines sites of potential failure in composite materials. Method is sensitive to matrix damage and fiber disbonding that occurs long before actual fiber breakage.

Whitcomb, J. D.↗

Nonlinear structural analysis

Development of computer programs for nonlinear structural analyses has progressed from special application programs to large, generalized programs. AGGIE I program is moderately-sized finite-element program that was developed specifically for nonlinear structural analysis. It is based on two- and three-dimensional isoparametric solid elements. AGGIE I accounts for nonlinearities due to large displacements, large strains and non-linear material behavior.

Haisler, W. E.↗

An integrated computer procedure for sizing composite airframe structures

A computerized algorithm to generate cross-sectional dimensions and fiber orientations for composite airframe structures is described, and its application in a wing structural synthesis is established. The algorithm unifies computations of aeroelastic loads, stresses, and deflections, as well as optimal structural sizing and fiber orientations in an open-ended system of integrated computer programs. A finite-element analysis and a mathematical-optimization technique are discussed.

Sobieszczanski-Sobieski, J.↗

User's guide to computer programs JET 5A and CIVM-JET 5B to calculate the large elastic-plastic dynamically-induced deformations of multilayer partial and/or complete structural rings

These structural ring deflections lie essentially in one plane and, hence, are called two-dimensional (2-d). The structural rings may be complete or partial; the former may be regarded as representing a fragment containment ring while the latter may be viewed as a 2-d fragment-deflector structure. These two types of rings may be either free or supported in various ways (pinned-fixed, locally clamped, elastic-foundation supported, mounting-bracket supported, etc.). The initial geometry of each ring may be circular or arbitrarily curved; uniform-thickness or variable-thickness rings may be analyzed. Strain-hardening and strain-rate effects of initially-isotropic material are taken into account. An approximate analysis utilizing kinetic energy and momentum conservation relations is used to predict the after-impact velocities of each fragment and of the impact-affected region of the ring; this procedure is termed the collision-imparted velocity method (CIVM) and is used in the CIVM-JET 5 B program. This imparted-velocity information is used in conjunction with a finite-element structural response computation code to predict the transient, large-deflection, elastic-plastic responses of the ring. Similarly, the equations of motion of each fragment are solved in small steps in time. Provisions are made in the CIVM-JET 5B code to analyze structural ring response to impact attack by from 1 to 3 fragments, each with its own size, mass, translational velocity components, and rotational velocity. The effects of friction between each fragment and the impacted ring are included.

Wu, R. W. H.↗

Analytical and scale model research aimed at improved hangglider design

Research consisted of a theoretical analysis which attempts to predict aerodynamic characteristics using lifting surface theory and finite-element structural analysis as well as an experimental investigation using 1/5 scale elastically similar models in the NASA Ames 2m x 3m (7' x 10') wind tunnel. Experimental data were compared with theoretical results in the development of a computer program which may be used in the design and evaluation of ultralight gliders.

Kroo, I.↗

Crack tip finite element analyses

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.

Swedlow, J. L.↗

Effect of crack size on the natural frequencies of a cracked plate

Results are presented for a finite element modal analysis of a centrally cracked rectangular plate made of linearly elastic material. The objective is to assess the effect of crack size on the natural vibration frequencies of the cracked plate. Only the in-plane vibration modes are studied. The results presented are finite-element mesh-size dependent. Namely, shrinking the mesh size, especially in the crack tip region, would change the magnitudes of the natural frequencies, and the trends would be the same as shown in the present note.

Ko, W. L.↗

Simplified and refined finite element approaches for determining stresses and internal forces in geometrically nonlinear structural analysis

Two methods for determining stresses and internal forces in geometrically nonlinear structural analysis are presented. The simplified approach uses the mid-deformed structural position to evaluate strains when rigid body rotation is present. The important feature of this approach is that it can easily be used with a general-purpose finite-element computer program. The refined approach uses element intrinsic or corotational coordinates and a geometric transformation to determine element strains from joint displacements. Results are presented which demonstrate the capabilities of these potentially useful approaches for geometrically nonlinear structural analysis.

Robinson, J. C.↗

Stress-intensity factors for a wide range of semi-elliptical surface cracks in finite-thickness plates

Surface cracks are among the more common flaws in aircraft and pressure vessel components. Several calculations of stress-intensity factors for semi-elliptical surface cracks subjected to tension have appeared in the literature. However, some of these solutions are in disagreement by 50-100%. In this paper, stress-intensity factors for shallow and deep semi-elliptical surface cracks in plates subjected to tension are presented. To verify the accuracy of the three-dimensional finite-element models employed, convergence was studied by varying the number of degrees of freedom in the models from 1500 to 6900. The 6900 degrees of freedom used here were more than twice the number used in previously reported solutions. Also, the stress-intensity variations in the boundary-layer region at the intersection of the crack with the free surface were investigated.

Raju, I. S.↗

Influence of two-dimensional hygrothermal gradients on interlaminar stresses near free edges

Interlaminar stresses are determined for mechanical loading, uniform hygrothermal loading, and gradient moisture loading through implementation of a finite-element computer code. Nonuniform two-dimensional hygroscopic gradients are obtained from a finite-difference solution of the diffusion equation. It is shown that hygroscopic-induced stresses can be larger than those resulting from mechanical and thermal loading, and that the distribution of the interlaminar normal stress may be changed significantly in the presence of a two-dimensional moisture gradient in the boundary layer of a composite laminate.

Farley, G. L.↗

Analysis of surface cracks in finite plates under tension or bending loads

Stress-intensity factors calculated with a three-dimensional, finite-element analysis for shallow and deep semielliptical surface cracks in finite elastic isotropic plates subjected to tension or bending loads are presented. A wide range of configuration parameters was investigated. The ratio of crack depth to plate thickness ranged from 0.2 to 0.8 and the ratio of crack depth to crack length ranged from 0.2 to 2.0. The effects of plate width on stress-intensity variations along the crack front was also investigated. A wide-range equation for stress-intensity factors along the crack front as a function of crack depth, crack length, plate thickness, and plate width was developed for tension and bending loads. The equation was used to predict patterns of surface-crack growth under tension or bending fatigue loads. A modified form of the equation was also used to correlate surface-crack fracture data for a brittle epoxy material within + or - 10 percent for a wide range of crack shapes and crack sizes.

Newman, J. C., Jr.↗

Thermal elastohydrodynamic lubrication of spur gears

An analysis and computer program called TELSGE were developed to predict the variations of dynamic load, surface temperature, and lubricant film thickness along the contacting path during the engagement of a pair of involute spur gears. The analysis of dynamic load includes the effect of gear inertia, the effect of load sharing of adjacent teeth, and the effect of variable tooth stiffness which are obtained by a finite-element method. Results obtained from TELSGE for the dynamic load distributions along the contacting path for various speeds of a pair of test gears show patterns similar to that observed experimentally. Effects of damping ratio, contact ratio, tip relief, and tooth error on the dynamic load were examined. In addition, two dimensionless charts are included for predicting the maximum equilibrium surface temperature, which can be used to estimate directly the lubricant film thickness based on well established EHD analysis.

Wang, K. L.↗

The Effects of Variable Mass and Geometry, Pretwist, Shear Deformation and Rotatory Inertia on the Resonant Frequencies of Intact Long Bones: A Finite Element Model Analysis

The influence of pretwist, nonuniformities in mass and flexural stiffness, rotatory inertia and shear deformation on the natural frequencies of intact bones is evaluated by means of a linear elastic, finite-element model which has been programmed for solution on the digital computer. Theoretical results are compared to the results on the forced vibration of intact canine radii obtained experimentally by Thompson. Surprisingly, inclusion of fairly large pretwist angles (from -14 to 12 deg for one specimen) had little affect on the first three frequencies of transverse vibration in either the cranial or lateral directions. Inclusion of shear deformation reduced the third-mode frequency in the stiffest (lateral) direction by about six percent, otherwise shear deformation played a minor role in determining natural frequencies. Similarly. rotatory inertia had negligible influence up to the third natural frequency. The predominant influence on the first three natural frequencies of transverse vibration could be attributed to the variations in mass and flexural stiffness along the length of the test specimens. Different effective moduli of elasticity are required to yield correct absolute values for the frequencies which correspond to experimental findings. thus implying the presence of some inhomogeneities in material properties around the bone cross-section and/or along its length.

Young, Donald R.↗