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 1,171 records · Page 65

Float Zone Experiments in Space

The objective of this work has been to evaluate whether or not Marangoni flow could be suppressed in molten metals by the presence of very thin oxide films. Experimental work has been carried out on molten Sn under UHV conditions. A disk floating zone arrangement was developed to allow in situ Auger examination of molten surfaces. An electron energy loss technique was developed which allows detection of continuous tin oxide films of 6 A or greater. Experiments were planned to detect the effects of oxide formation upon Marangoni flow by measuring: (1) temperature profiles, (2) solid liquid interface shapes, (3) macrosegregation, and (4) the onset of oscillatory Marangoni flow by detecting oscillating temperature variations. Work on (4) showed that oscillatory temperature variations of frequency or = 10 Hz were not present in the disk float zone geometry under conditions of Ma = 4300 with an oxide free molten surface. The disk float zone geometry was modeled with a finite element analysis and temperature and velocity profiles were determined.

Verhoeven, J. D.↗

Mixed-mode cyclic debonding of adhesively bonded composite joints

A combined experimental-analytical investigation to characterize the cyclic failure mechanism of a simple composite-to-composite bonded joint is conducted. The cracked lap shear (CLS) specimens of graphite/epoxy adherend bonded with EC-3445 adhesive are tested under combined mode 1 and 2 loading. In all specimens tested, fatigue failure occurs in the form of cyclic debonding. The cyclic debond growth rates are measured. The finite element analysis is employed to compute the mode 1, mode 2, and total strain energy release rates (i.e., GI, GII, and GT). A wide range of mixed-mode loading, i.e., GI/GII ranging from 0.03 to 0.38, is obtained. The total strain energy release rate, G sub T, appeared to be the driving parameter for cyclic debonding in the tested composite bonded system.

Rezaizadeh, M. A.↗

Interlaminar fracture toughness: Three-dimensional finite element modeling for end-notch and mixed-mode flexure

A computational procedure is described for evaluating End-Notch-Flexure (ENF) and Mixed-Mode-Flexure (MMF) interlaminar fracture toughness in unidirectional fiber composites. The procedure consists of a three-dimensional finite element analysis in conjunction with the strain energy release rate concept and with composite micromechanics. The procedure is used to analyze select cases of ENF and MMF. The strain energy release rate predicted by this procedure is in good agreement with limited experimental data. The procedure is used to identify significant parameters associated with interlaminar fracture toughness. It is also used to determine the critical strain energy release rate and its attendant crack length in ENF and/or MMF. This computational procedure has considerable versatility/generality and provides extensive information about interlaminar fracture toughness in fiber composites.

Murthy, P. L. N.↗

Application of mathematical optimization procedures to a structural model of a large finite-element wing

An optimization system is described that consists of the finite-element analysis program known as Engineering Analysis Language (EAL), the CONMIN optimization program, and a piecewise linear analysis. The optimization analysis system (OPTEAL) is used to minimize the mass of a transport wing subject to stress constraints for two static load conditions. The optimization procedure and the optimization analysis system (OPTEAL) are described and the optiomized results with fully stressed design results compared.

Walsh, J. L.↗

Experimental and analytical generic space station dynamic models

A dynamic model used for verification of analytical and experimental methods is documented. The model consists of five substructures to simulate the multibody, low frequency nature of large space structures. Design considerations which led to a fundamental vibration frequency of less than one Hz are described. Finite element analysis used to predict the vibration modes and frequencies of the experimental model is presented. In addition, modeling of cable suspension effects using prestressed vibration analysis is described. Details of the expermental and analytical models are included to permit replication of the study. Results of the modal vibration tests and analysis are presented in a separate document.

Belvin, W. K.↗

Combined bearing and bypass loading on a graphite/epoxy laminate

A combined experimental and analytical study was conducted to determine the behavior of a graphite/epoxy laminate subjected to combined bearing and bypass loading. Single-fastener quasi-isotropic specimens were loaded at various bearing-bypass ratios until damage was produced at the fastener hole. Damage-onset strengths and damage modes were then analyzed using local hole-boundary stresses calculated by a finite-element analysis. The tension data showed the expected linear interaction for combined bearing and bypass loading with damage developing in the net-section tension mode. However, the compression bearing-bypass strengths showed an unexpected interaction involving the bearing mode. Compressive bypass loads reduced the bearing strength by decreasing the bolt-hole contact arc and thus increasing the severity of the bearing loads. The bearing stresses at the hole boundary were not accurately estimated by superposition of the stress components for separate bearing and bypass loading. However, superposition produced reasonably accurate estimates for tangential stresses especially near the specimen net-section.

Crews, J. H., Jr.↗

A comparison of NASTRAN (Cosmic) and experimental results for the vibration of thick open cylindrical cantilevered shells

The natural frequencies and associated mode shapes for three thick open cantilevered cylindrical shells were determined both numerically and experimentally. The shells ranged in size from moderately to very thick with length to thickness ratios of 16, 8 and 5.6, the independent dimension being the shell thickness. The shell geometry is characterized by a circumferential angle of the 142 degrees and a ratio of length to inner radii arc length near 1.0. The finite element analysis was performed using NASTRAN's (COSMIC) triangular plate bending element CTRIA2, which includes membrane effects. The experimental results were obtained through holographic interferometry which enables one to determine the resonant frequencies as well as mode shapes from photographs of time-averaged holograms.

Mindle, W. L.↗

Computational simulation of progressive fracture in fiber composites

Computational methods for simulating and predicting progressive fracture in fiber composite structures are presented. These methods are integrated into a computer code of modular form. The modules include composite mechanics, finite element analysis, and fracture criteria. The code is used to computationally simulate progressive fracture in composite laminates with and without defects. The simulation tracks the fracture progression in terms of modes initiating fracture, damage growth, and imminent global (catastrophic) laminate fracture.

Chamis, C. C.↗

Multigrid methods in structural mechanics

Although the application of multigrid methods to the equations of elasticity has been suggested, few such applications have been reported in the literature. In the present work, multigrid techniques are applied to the finite element analysis of a simply supported Bernoulli-Euler beam, and various aspects of the multigrid algorithm are studied and explained in detail. In this study, six grid levels were used to model half the beam. With linear prolongation and sequential ordering, the multigrid algorithm yielded results which were of machine accuracy with work equivalent to 200 standard Gauss-Seidel iterations on the fine grid. Also with linear prolongation and sequential ordering, the V(1,n) cycle with n greater than 2 yielded better convergence rates than the V(n,1) cycle. The restriction and prolongation operators were derived based on energy principles. Conserving energy during the inter-grid transfers required that the prolongation operator be the transpose of the restriction operator, and led to improved convergence rates. With energy-conserving prolongation and sequential ordering, the multigrid algorithm yielded results of machine accuracy with a work equivalent to 45 Gauss-Seidel iterations on the fine grid. The red-black ordering of relaxations yielded solutions of machine accuracy in a single V(1,1) cycle, which required work equivalent to about 4 iterations on the finest grid level.

Raju, I. S.↗

Large Deformation Dynamic Bending of Composite Beams

Studies were conducted on the large deformation response of composite beams subjected to a dynamic axial load. The beams were loaded with a moderate eccentricity to promote bending. The study was primarily experimental but some finite element results were obtained. Both the deformation and the failure of the beams were of interest. The static response of the beams was also studied to determine potential differences between the static and dynamic failure. Twelve different laminate types were tested. The beams were loaded dynamically with a gravity driven impactor traveling at 19.6 ft/sec and quasi-static tests were conducted on identical beams in a displacement controlled manner. For laminates of practical interest, the failure modes under static and dynamic loadings were identical. Failure in most of the laminate types occurred in a single event involving 40% to 50% of the plies. However, failure in laminates with 30 deg or 15 deg off-axis plies occured in several events. All laminates exhibited bimodular elastic properties. Using empirically determined flexural properties, a finite element analysis was reasonably accurate in predicting the static and dynamic deformation response.

Derian, E. J.↗

The use of finite element methods for determining the cross-sectional temperature distribution in heat pipes

A model that is currently used to predict the priming and performance limitations of a monogroove heat pipe is expanded to include the boiling limitation and the cross-sectional temperature distribution as determined from a multidimensional finite element analysis technique. The improved model is verified experimentally and shown to accurately predict the cross-sectional temperature distribution when the heat flux distribution is known. The model provides a way to estimate the level at which nucleate boiling and the associated dryout of the capillary wick occurs.

Peterson, G. P.↗

Longitudinal compressive failure modes in fiber composites End attachment effects on IITRI type test specimens

The end-attachment effects on longitudinal compressive strength of IITRI type specimen unidirectional fiber composites are formally assessed using finite-element analysis (FEA) in conjunction with composite mechanics. Sixteen different cases were analyzed to evaluate end-attachment effects (such as degree of misalignment, type of misalignment, progressive end-tab debonding, and specimen thickness) on stress distribution, peak stresses, buckling loads, and buckling mode shapes. The results obtained from the FEA and comparisons with fractured specimens show that eccentricities induce bending-type stresses which peak near the end-tabs and cause flexural type fracture. Also, guidelines are included for placing back-to-back strain gages to measure the presence/absence of possible end-attachment and eccentricity effects.

Chamis, C. C.↗

Characterization of Mode I and mixed-mode delamination growth in T300/5208 graphite/epoxy

The roles played by Mode I and Mode II strain-energy release rates (G-I and G-II, respectively) in inducing delamination growth under static and fatigue loading were investigated, using T300/5208 graphite/epoxy specimens. Double cantilever beam (DCB) specimens and cracked lap shear (CLS) specimens were used for pure Mode I and mixed-mode tests, respectively. Fatigue-induced delamination growth was characterized by constant-amplitude fatigue tests at a minimum to maximum cyclic load ratio of 0.05 and a frequency of 10 Hz. During the tests, the maximum and minimum strain-energy release rates (Gmax, Gmin) and the delamination growth rate (da/dN) were monitored. Static tests on mixed-mode CLS specimens measured the total strain-energy release rate, which was broken into G-I and G-II components using finite-element analysis. A power-law relationship between da/dN and G-Imax, and da/dN and Gmax were obtained from fatigue test results on DCB and CLS specimens, respectively. The power law for a pure Mode II delamination was derived from CLS results by subtracting the contribution due to G-I.

Ramkumar, R. L.↗

Postbuckling behavior of selected graphite-epoxy cylindrical panels loaded in axial compression

Results of an experimental and analytical study of the postbuckling behavior of selected graphite-epoxy cylindrical panels loaded in axial compression are presented. The postbuckling response and failure characteristics of the panels are described. The postbuckling response of each specimen is typical of axially-compressed cylindrical shells and curved panels in that a severe reduction in load occurs at buckling. Failure of all panels initiated near regions with severe local bending gradients. Analytical results from a nonlinear general shell finite element analysis computer code correlate well with typical experimental results up to buckling. Measured initial geometric imperfections were included in the postbuckling analysis. Analytically-determined stress distributions in the postbuckling response were used with failure criteria to identify the load level and the location of first-ply failure.

Knight, N. F., Jr.↗

An analytical treatment of discretely varying constraints and inertial properties in multi-body dynamics

A computerized procedure is presented for calculating instantaneous velocity changes due to discrete imposition of constraints, impact or discrete inertial changes. Such discrete changes are impulsive in nature, involve wave propagation and therefore usually require detailed modeling for accurate prediction. It is demonstrated that when modeling is coarse, physically unrealistic results can be produced. The present procedure, which concentrates on finite element analysis, reduces the need for detailed modeling by assuming that only those velocities in the neighborhood of grid points involved in the constraint or those rigidly connected to grid points involved in the constraint need be considered in calculating velocity changes. The procedure utilizes Newton's second law and as such conserves momentum where valid. Sample problems are provided.

Housner, J. M.↗

Dynamic delamination crack propagation in a graphite/epoxy laminate

The dynamic delamination crack propagation behavior during ballistic tests of (90/0)5s T-300/934 graphite/epoxy laminates with embedded interfacial cracks was investigated using high speed photography. The impact on the beam-like specimen was produced with a silicon rubber ball, and the crack propagation speeds and the threshold impact velocities required to initiate dynamic crack propagation were determined for several crack positions. The results suggest that the mode of crack propagation depends on the specimen geometry as well as the loading condition. A simplified finite element analysis of the experimental data obtained from one of the midplane-cracked specimens was used to estimate the critical strain energy release rate, which may determine the onset of unstable crack propagation.

Sun, C. T.↗

Some observations on the behavior of the Langley model rotor blade

The design of the model rotor and the comparative study of coupled beam theory and the finite element analysis performed earlier at the Aerostructures Directorate by Robert Hodges and Mark Nixon is examined. Attention is focused upon two matters: (1) an examination of the small discrepancies between twist angle predictions under pure torque and radial loading, and (2) an assessment of nonclassical effects in bending behavior. The primary objective is understanding, particularly with regard to cause and effect relationships. Understanding, together with the simple, affordable nature of the coupled beam analysis, provides a sound basis for design.

Rehfield, L. W.↗

STAEBL: Structural tailoring of engine blades, phase 2

The Structural Tailoring of Engine Blades (STAEBL) program was initiated at NASA Lewis Research Center in 1980 to introduce optimal structural tailoring into the design process for aircraft gas turbine engine blades. The standard procedure for blade design is highly iterative with the engineer directly providing most of the decisions that control the design process. The goal of the STAEBL program has been to develop an automated approach to generate structurally optimal blade designs. The program has evolved as a three-phase effort with the developmental work being performed contractually by Pratt & Whitney Aircraft. Phase 1 was intended as a proof of concept in which two fan blades were structurally tailored to meet a full set of structural design constraints while minimizing DOC+I (direct operating cost plus interest) for a representative aircraft. This phase was successfully completed and was reported in reference 1 and 2. Phase 2 has recently been completed and is the basis for this discussion. During this phase, three tasks were accomplished: (1) a nonproprietary structural tailoring computer code was developed; (2) a dedicated approximate finite-element analysis was developed; and (3) an approximate large-deflection analysis was developed to assess local foreign object damage. Phase 3 is just beginning and is designed to incorporated aerodynamic analyses directly into the structural tailoring system in order to relax current geometric constraints.

Hirschbein, M. S.↗