Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixed 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1966

Wind-tunnel boundary interference for V/STOL TESTING

The wake skew angle used in applying the theory of NASA TR R-124 to data correction should be such that the angular deflection of the wake vorticity from the horizontal is one-half that calculated from momentum theory at the lifting element. This usage is in contrast to that of the original paper which used the angle of the mass flow. Because of large-scale recirculation effects, there is a finite lower limit to the test speed at which reliable and correctable data can be obtained in closed wind tunnels. Although a zero-correction wind tunnel for V/STOL testing has not yet been achieved, it is shown that the use of suitably mixed wind-tunnel boundaries can alleviate boundary effects on V/STOL data.

V/STOL AIRCRAFT↗
NASA NTRS1985

A first-order Green's function approach to supersonic oscillatory flow: A mixed analytic and numeric treatment

A frequency domain Green's Function Method for unsteady supersonic potential flow around complex aircraft configurations is presented. The focus is on the supersonic range wherein the linear potential flow assumption is valid. In this range the effects of the nonlinear terms in the unsteady supersonic compressible velocity potential equation are negligible and therefore these terms will be omitted. The Green's function method is employed in order to convert the potential flow differential equation into an integral one. This integral equation is then discretized, through standard finite element technique, to yield a linear algebraic system of equations relating the unknown potential to its prescribed co-normalwash (boundary condition) on the surface of the aircraft. The arbitrary complex aircraft configuration (e.g., finite-thickness wing, wing-body-tail) is discretized into hyperboloidal (twisted quadrilateral) panels. The potential and co-normalwash are assumed to vary linearly within each panel. The long range goal is to develop a comprehensive theory for unsteady supersonic potential aerodynamic which is capable of yielding accurate results even in the low supersonic (i.e., high transonic) range.

Freedman, M. I.↗
NASA NTRS1990

A Poisson equation formulation for pressure calculations in penalty finite element models for viscous incompressible flows

The calculation of pressures when the penalty-function approximation is used in finite-element solutions of laminar incompressible flows is addressed. A Poisson equation for the pressure is formulated that involves third derivatives of the velocity field. The second derivatives appearing in the weak formulation of the Poisson equation are calculated from the C0 velocity approximation using a least-squares method. The present scheme is shown to be efficient, free of spurious oscillations, and accurate. Examples of applications are given and compared with results obtained using mixed formulations.

Sohn, J. L.↗
NASA NTRS1998

Fracture Mechanics Analysis of Stitched Stiffener-Skin Debonding

An analysis based on plate finite elements and the virtual crack closure technique has been implemented to study the effect of stitching on mode I and mode II strain energy release rates for debond configurations. The stitches were modeled as discrete nonlinear fastener elements with a compliance determined by experiment. The axial and shear behavior of the stitches was considered, however, the two compliances and failure loads were assumed to be independent. Both a double cantilever beam (mode I) and a mixed mode skin-stiffener debond configuration were studied. In the double cantilever beam configurations, G(sub I) began to decrease once the debond had grown beyond the first row of stitches and was reduced to zero for long debonds. In the mixed-mode skin-stiffener configurations, G(sub I) showed a similar behavior as in the double cantilever beam configurations, however, G(sub u), continued to increase with increasing debond length.

Glaessgen, E. H.↗
NASA NTRS2014

A Numerical and Experimental Study of Damage Growth in a Composite Laminate

The present study has three goals: (1) perform an experiment where a simple laminate damage process can be characterized in high detail; (2) evaluate the performance of existing commercially available laminate damage simulation tools by modeling the experiment; (3) observe and understand the underlying physics of damage in a composite honeycomb sandwich structure subjected to low-velocity impact. A quasi-static indentation experiment has been devised to provide detailed information about a simple mixed-mode damage growth process. The test specimens consist of an aluminum honeycomb core with a cross-ply laminate facesheet supported on a stiff uniform surface. When the sample is subjected to an indentation load, the honeycomb core provides support to the facesheet resulting in a gradual and stable damage growth process in the skin. This enables real time observation as a matrix crack forms, propagates through a ply, and then causes a delamination. Finite element analyses were conducted in ABAQUS/Explicit(TradeMark) 6.13 that used continuum and cohesive modeling techniques to simulate facesheet damage and a geometric and material nonlinear model to simulate core crushing. The high fidelity of the experimental data allows a detailed investigation and discussion of the accuracy of each numerical modeling approach.

McElroy, Mark↗
NASA NTRS2017

Computationally Efficient Beam Elements for Accurate Stresses in Sandwich Laminates and Laminated Composites with Delaminations

Laminated composites are prone to delamination failure due to the lack of reinforcement through the thickness. Therefore, during the design process the initiation and propagation of delaminations should be accounted for as early as possible. This paper presents computationally efficient nine degree-of-freedom (dof) and eight-dof shear locking-free beam elements using the mixed form of the refined zigzag theory (RZT(sup(m))). The corresponding nine-dof and eight-dof elements use the anisoparametric and constrained anisoparametric interpolation schemes, respectively, to eliminate shear locking in slender beams. The advantage of the present element over previous RZT beam elements is that no post-processing is required to accurately model the transverse shear stress while maintaining the computational efficiency of a low-order beam element. Comparisons with high-fidelity finite element models and three-dimensional elasticity solutions show that the elements can robustly and accurately predict the displacement field, axial stress and transverse shear stress through the thickness of a sandwich beam or a composite laminate with an embedded delamination. In fact, the accuracy and computational efficiency of predicting stresses in laminates with embedded delaminations make the present elements attractive choices for RZT-based delamination initiation and propagation methodologies available in the literature.

R. M. J. Groh↗
NASA NTRS2014

TPSAS-NF1676L-19616-DND

The present study has three goals: (1) perform an experiment where a simple laminate damage process can be characterized in high detail; (2) evaluate the performance of existing commercially available laminate damage simulation tools by modeling the experiment; (3) observe and understand the underlying physics of damage in a composite honeycomb sandwich structure subjected to low-velocity impact. A quasi-static indentation experiment has been devised to provide detailed information about a simple mixed-mode damage growth process. The test specimens consist of an aluminum honeycomb core with a cross-ply laminate facesheet supported on a stiff uniform surface. When the sample is subjected to an indentation load, the honeycomb core provides support to the facesheet resulting in a gradual and stable damage growth process in the skin. This enables real time observation as a matrix crack forms, propagates through a ply, and then causes a delamination. Finite element analyses were conducted in ABAQUS/Explicit(TradeMark) 6.13 that used continuum and cohesive modeling techniques to simulate facesheet damage and a geometric and material nonlinear model to simulate core crushing. The high fidelity of the experimental data allows a detailed investigation and discussion of the accuracy of each numerical modeling approach.

Mark McElroy↗
NASA NTRS1980

ACCESS 3. Approximation concepts code for efficient structural synthesis: User's guide

A user's guide is presented for ACCESS-3, a research oriented program which combines dual methods and a collection of approximation concepts to achieve excellent efficiency in structural synthesis. The finite element method is used for structural analysis and dual algorithms of mathematical programming are applied in the design optimization procedure. This program retains all of the ACCESS-2 capabilities and the data preparation formats are fully compatible. Four distinct optimizer options were added: interior point penalty function method (NEWSUMT); second order primal projection method (PRIMAL2); second order Newton-type dual method (DUAL2); and first order gradient projection-type dual method (DUAL1). A pure discrete and mixed continuous-discrete design variable capability, and zero order approximation of the stress constraints are also included.

Fleury, C.↗
NASA NTRS1981

International Conference on Numerical Methods in Fluid Dynamics, 7th, Stanford University, Stanford and Moffett Field, CA, June 23-27, 1980, Proceedings

Topics discussed include polygon transformations in fluid mechanics, computation of three-dimensional horseshoe vortex flow using the Navier-Stokes equations, an improved surface velocity method for transonic finite-volume solutions, transonic flow calculations with higher order finite elements, the numerical calculation of transonic axial turbomachinery flows, and the simultaneous solutions of inviscid flow and boundary layer at transonic speeds. Also considered are analytical solutions for the reflection of unsteady shock waves and relevant numerical tests, reformulation of the method of characteristics for multidimensional flows, direct numerical simulations of turbulent shear flows, the stability and separation of freely interacting boundary layers, computational models of convective motions at fluid interfaces, viscous transonic flow over airfoils, and mixed spectral/finite difference approximations for slightly viscous flows.

Reynolds, W. C.↗
NASA NTRS1988

Effect of vertical magnetic field on convection and segregation in vertical Bridgman crystal growth

A previous finite-element analysis of vertical Bridgman growth for dilute and nondilute alloys is extended to include the effect of a vertically-aligned magnetic field in the limit of zero magnetic Reynolds number. Calculations are presented for growth of a dilute gallium-germanium alloy in a vertically stabilized Bridgman-Stockbarger system and in a furnace with a uniform temperature gradient imposed along the ampoule. Steady cellular convection driven by radial temperature gradients causes good axial and radial mixing in both systems without a magnetic field. A weak magnetic field decreases the intensity of convection and the effectiveness of solute mixing. The radial nonuniformity is greatest for an intermediate field strength. Stronger fields suppress flow recirculation completely, and lead to uniform solute segregation across the crystal and to diffusion-controlled axial segregation.

Kim, Do Hyun↗
NASA NTRS1981

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

Analysis of X-Ray Microradiographs of Al-Au Interface Quench Profile using Modeling of Solidification Including Double-Diffusion and Convection in the Melt

Experimental data on Al-0.8Au horizontal solidification of a 1 mm thick specimen in a BN crucible shows the effect of growth rate on the solidification interface shape. For translation rates below 0.5 micron/s the interface maintains a plain and flat shape. When the translation rate is 3 to 5 micron/s or more, the interface appearance changes to two planar zones, with the zone closer to the bottom having higher inclination. The interface shapes were measured by first quenching in place during growth. X-ray microscopy shows the interface shape within the quenched sample by viewing through the side of the specimen. In order to provide theoretical explanation of the phenomena, numerical modeling was undertaken using finite element code FIDAP. Double diffusion convection in Al-0.8Au melt and crystal-melt interface curvature during directional solidification was analyzed numerically. Actual thermophysical properties of Al-0.8Au including the binary Al-Au phase diagram were used. Although convection in the sample is weak, for the slower translation rate convection and diffusion is sufficient for the redistribution of initial compositional stratification caused by gravity. When translation rate is raised, neither convection nor diffusion can provide proper mixing so that solidification temperatures differ significantly near the bottom within the bulk of the sample. As a result, the solid-liquid interface appears to have two planar zones with different inclination.

Bune, Andris V.↗
NASA NTRS1983

Advances and trends in structural and solid mechanics; Proceedings of the Symposium, Washington, DC, October 4-7, 1982

The mechanics of materials and material characterization are considered, taking into account micromechanics, the behavior of steel structures at elevated temperatures, and an anisotropic plasticity model for inelastic multiaxial cyclic deformation. Other topics explored are related to advances and trends in finite element technology, classical analytical techniques and their computer implementation, interactive computing and computational strategies for nonlinear problems, advances and trends in numerical analysis, database management systems and CAD/CAM, space structures and vehicle crashworthiness, beams, plates and fibrous composite structures, design-oriented analysis, artificial intelligence and optimization, contact problems, random waves, and lifetime prediction. Earthquake-resistant structures and other advanced structural applications are also discussed, giving attention to cumulative damage in steel structures subjected to earthquake ground motions, and a mixed domain analysis of nuclear containment structures using impulse functions.

Noor, A. K.↗
NASA NTRSDate not supplied

Impact of Thermoplastic Composites: Testing and Modeling

Thermoplastic Composites (TPCs) are being increasingly considered for aerospace applications given their faster manufacturing cycles and lower cost. The NASA High-Rate Composite Aircraft Manufacturing (HiCAM) project aims to evaluate and mature thermoplastic manufacturing technologies to achieve a 4X – 6X increase in the production rate of Next Generation Single Aisle (NGSA) commercial aircraft. Within the purview of the NASA HiCAM project, this work investigates the response of thermoplastic composites under High Energy Dynamic Impact (HEDI) conditions. HEDI tests were conducted on panels fabricated with a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK) material system. A blunt metal projectile was chosen as the impactor which impacted the T700/LMPAEK panels over a range of impact velocities imparted by a gas gun test setup. As expected, lower velocity impacts caused the projectile to rebound, whereas higher velocity impacts resulted in the projectile penetrating the test panels. For the test cases in which the projectile rebounded, the predominant damage modes recorded with Ultrasonic (UT) scans included interlaminar delaminations which exhibited a “rotating fan” like structure when viewed in the through-thickness direction. Additionally, the dynamic deflection of the center point of the back face of the test panels and the velocity of the projectile was also measured during the test. Finite element models were developed to model a test case in which the projectile rebounded. These models were developed using LS-DYNA® wherein the primary objective was to evaluate the ability of the material model MAT299 to capture the dynamic response and damage modes in the T700/LMPAEK test panels. MAT299 is a Deformation Gradient Decomposition (DGD) based Continuum Damage Mechanics (CDM) material model. The interlaminar delaminations were modeled using the cohesive contact formulation available within LS-DYNA, wherein the mixed-mode fracture is captured via the Benzegaggh-Kenane (B-K) mode-dependent fracture energy interpolation law. The above-mentioned modeling methods have mostly been applied to thermoset composites and adapting them to TPCs involves addressing challenges associated with appropriately representing material behavior. Therefore, the proposed paper would discuss the systematic approach taken to adapt thermoset composites modeling practices for applications to TPCs, while appropriately addressing material behavior. The proposed paper would show that the predicted back face center point deflection correlated reasonably well with experimental data and the peak deflection was predicted to be within 5% of the experimental measurements. The projectile rebound velocity, while predicted to be higher than the experimental measurement, was within reasonable bounds. Additionally, it will be shown that the predicted delamination shapes were in reasonable agreement with experimental data.

Composite Materials↗
NASA NTRS1985

Delamination and debonding of materials

The general topics consist of stress analysis, mechanical behavior, and fractography/NDI of composite laminates. Papers are presented on a dynamic hybrid finite-element analysis for interfacial cracks in composites, energy release rate during delamination crack growth in composite laminates, matrix deformation and fracture in graphite-reinforced epoxies, and the role of delamination and damage development on the strength of thick notched laminates. In addition, consideration is given to a new ply model for interlaminar stress analysis, a fracture mechanics approach for designing adhesively bonded joints, the analysis of local delaminations and their influence on composite laminate behavior, and moisture and temperature effects on the mixed-mode delamination fracture of unidirectional graphite/epoxy.

Johnson, W. S.↗
NASA NTRS1999

A Viscoelastic Hybrid Shell Finite Element

An elastic large displacement thick-shell hybrid finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at he element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses in included in the mixed variational functional. Nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to numerically simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

Johnson, Arthur↗
NASA NTRS1987

Interlaminar shear fracture toughness and fatigue thresholds for composite materials

Static and cyclic end notched flexure tests were conducted on a graphite epoxy, a glass epoxy, and graphite thermoplastic to determine their interlaminar shear fracture toughness and fatigue thresholds for delamination in terms of limiting values of the mode II strain energy release rate, G-II, for delamination growth. The influence of precracking and data reduction schemes are discussed. Finite element analysis indicated that the beam theory calculation for G-II with the transverse shear contribution included was reasonably accurate over the entire range of crack lengths. Cyclic loading significantly reduced the critical G-II for delamination. A threshold value of the maximum cyclic G-II below which no delamination occurred after one million cycles was identified for each material. Also, residual static toughness tests were conducted on glass epoxy specimens that had undergone one million cycles without delamination. A linear mixed-mode delamination criteria was used to characterize the static toughness of several composite materials; however, a total G threshold criterion appears to characterize the fatigue delamination durability of composite materials with a wide range of static toughness.

Obrien, T. Kevin↗
NASA NTRS1989

Interlaminar shear fracture toughness and fatigue thresholds for composite materials

Static and cyclic end notched flexure tests were conducted on a graphite epoxy, a glass epoxy, and graphite thermoplastic to determine their interlaminar shear fracture toughness and fatigue thresholds for delamination in terms of limiting values of the mode II strain energy release rate, G-II, for delamination growth. The influence of precracking and data reduction schemes are discussed. Finite element analysis indicated that the beam theory calculation for G-II with the transverse shear contribution included was reasonably accurate over the entire range of crack lengths. Cyclic loading significantly reduced the critical G-II for delamination. A threshold value of the maximum cyclic G-II below which no delamination occurred after one million cycles was identified for each material. Also, residual static toughness tests were conducted on glass epoxy specimens that had undergone one million cycles without delamination. A linear mixed-mode delamination criteria was used to characterize the static toughness of several composite materials; however, a total G threshold criterion appears to characterize the fatigue delamination durability of composite materials with a wide range of static toughness.

O'Brien, T. Kevin↗