Search NASA⌕ Search

SEARCH · Search NASA

Results for “anisotropic elasticity”

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 145 records · Page 8

Nonlinear Elastic Effects on the Energy Flux Deviation of Ultrasonic Waves in GR/EP Composites

In isotropic materials, the direction of the energy flux (energy per unit time per unit area) of an ultrasonic plane wave is always along the same direction as the normal to the wave front. In anisotropic materials, however, this is true only along symmetry directions. Along other directions, the energy flux of the wave deviates from the intended direction of propagation. This phenomenon is known as energy flux deviation and is illustrated. The direction of the energy flux is dependent on the elastic coefficients of the material. This effect has been demonstrated in many anisotropic crystalline materials. In transparent quartz crystals, Schlieren photographs have been obtained which allow visualization of the ultrasonic waves and the energy flux deviation. The energy flux deviation in graphite/epoxy (gr/ep) composite materials can be quite large because of their high anisotropy. The flux deviation angle has been calculated for unidirectional gr/ep composites as a function of both fiber orientation and fiber volume content. Experimental measurements have also been made in unidirectional composites. It has been further demonstrated that changes in composite materials which alter the elastic properties such as moisture absorption by the matrix or fiber degradation, can be detected nondestructively by measurements of the energy flux shift. In this research, the effects of nonlinear elasticity on energy flux deviation in unidirectional gr/ep composites were studied. Because of elastic nonlinearity, the angle of the energy flux deviation was shown to be a function of applied stress. This shift in flux deviation was modeled using acoustoelastic theory and the previously measured second and third order elastic stiffness coefficients for T300/5208 gr/ep. Two conditions of applied uniaxial stress were considered. In the first case, the direction of applied uniaxial stress was along the fiber axis (x3) while in the second case it was perpendicular to the fiber axis along the laminate stacking direction (x1).

Prosser, William H.↗

Subsurface Stress Fields In Single Crystal (Anisotropic) Contacts

Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to conditions of high temperature, triaxial steady and alternating stresses, fretting stresses in the blade attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades are also subjected to extreme variations in temperature during start-up and shutdown transients. The most prevalent HCF failure modes observed in these blades during operation include crystallographic crack initiation/propagation on octahedral planes, and noncrystallographic initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation at the blade leading edge tip, blade attachment regions, and damper contact locations have been documented. Understanding crack initiation/propagation under mixed-mode loading conditions is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine blades. This paper presents analytical and numerical techniques for evaluating two and three dimensional subsurface stress fields in anisotropic contacts. The subsurface stress results are required for evaluating contact fatigue life at damper contacts and dovetail attachment regions in single crystal nickel-base superalloy turbine blades. An analytical procedure is , presented, for evaluating the subsurface stresses in the elastic half-space, using a complex potential method outlined by Lekhnitskii. Numerical results are presented for cylindrical and spherical anisotropic contacts, using finite element analysis. Effects of crystal orientation on stress response and fatigue life are examined.

Arakere, Nagaraj K.↗

Three-dimensional finite-element analysis of layered composite plates

Results are presented for an investigation of the three-dimensional, geometrically nonlinear, finite-element analysis of the bending of laminated anisotropic composite plates. The individual laminae are treated as homogeneous, transversely isotropic, and linearly elastic. A fully three-dimensional isoparametric finite element with eight modes (i.e., linear element) and 24 degrees of freedom (three displacement components per node) is used. The numerical results obtained using this linear analysis are compared with the exact solutions given in Pagano (1969, 1970). It is found that the results of the linear analysis converge to the exact solution as the mesh is refined.

Putcha, N. S.↗

A three-dimensional nonlinear analysis of cross-ply rectangular composite plates

The results of a three-dimensional, geometrically nonlinear, finite-element analysis of the bending of cross-ply laminated anisotropic composite plates are presented. Individual laminae are assumed to be homogeneous, orthotropic and linearly elastic. A fully three-dimensional isoparametric finite element with eight nodes (i.e. linear element) and 24 degrees of freedom (three displacement components per node) is used to model the laminated plate. The finite element results of the linear analysis are found to agree very well with the exact solutions of cross-ply laminated rectangular plates under sinusoidal loading. The finite element results of the three-dimensional, geometrically nonlinear analysis are compared with those obtained by using a shear deformable, geometrically nonlinear, plate theory. It is found that the deflections predicted by the shear deformable plate theory are in fair agreement with those predicted by three-dimensional elasticity theory; however stresses were found to be not in good agreement.

Kuppusamy, T.↗

Subsurface Stress Fields in FCC Single Crystal Anisotropic Contacts

Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to conditions of high temperature, triaxial steady and alternating stresses, fretting stresses in the blade attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades are also subjected to extreme variations in temperature during start-up and shutdown transients. The most prevalent high cycle fatigue (HCF) failure modes observed in these blades during operation include crystallographic crack initiation/propagation on octahedral planes, and non-crystallographic initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation at the blade leading edge tip, blade attachment regions, and damper contact locations have been documented. Understanding crack initiation/propagation under mixed-mode loading conditions is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine blades. This paper presents analytical and numerical techniques for evaluating two and three dimensional subsurface stress fields in anisotropic contacts. The subsurface stress results are required for evaluating contact fatigue life at damper contacts and dovetail attachment regions in single crystal nickel-base superalloy turbine blades. An analytical procedure is presented for evaluating the subsurface stresses in the elastic half-space, based on the adaptation of a stress function method outlined by Lekhnitskii. Numerical results are presented for cylindrical and spherical anisotropic contacts, using finite element analysis (FEA). Effects of crystal orientation on stress response and fatigue life are examined. Obtaining accurate subsurface stress results for anisotropic single crystal contact problems require extremely refined three-dimensional (3-D) finite element grids, especially in the edge of contact region. Obtaining resolved shear stresses (RSS) on the principal slip planes also involves considerable post-processing work. For these reasons it is very advantageous to develop analytical solution schemes for subsurface stresses, whenever possible.

Arakere, Nagaraj K.↗

Computation of Effective Mechanical Properties and Mechanical Erosion Modeling of TPS Materials

The goal of this presentation is to provide a general overview of the multi-scale modeling formulation to determine if there is additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. This modeling process is performed at different scales by leveraging two computational frameworks developed at NASA: the Porous Microstructure Analysis (PuMA) software, and the Porous material Analysis Toolbox based on OpenFOAM (PATO). PuMA specializes in computing effective macro-scale material properties by performing material response simulations on 3D digital micro-scale representations of porous micro-structures. The modeling of TPS materials at the micro-scale is essential to understand how they behave as part of a heat shield assembly. The first part of the presentation will detail the implementation of PuMA’s cell-centered finite volume elasticity solver, which allows the computation of macro-scale effective mechanical properties of heterogeneous and anisotropic materials such as fibrous and woven TPS composites. These homogenized mechanical properties are used by PATO’s mechanical erosion model to predict the TPS material’s recession at a larger scale. This work will also provide some examples of multi-scale analysis from the fiber level up to the unit cell. The second part of the presentation will focus on the macro-scale approach to determine if erosion at the heat shield’s surface occurs due to mechanical and thermal loads experienced during atmospheric entry. To accomplish this, a solid mechanics module was integrated within PATO enabling it to model the potential mechanical erosion in three steps: first, after obtaining the effective mechanical properties with PuMA, the implemented stress analysis solver computes the stress and the displacement fields for the TPS material using the wall shear stress tensor, computed using a CFD solver, as boundary conditions; then, regions on the surface where the stress meets the failure criteria are identified; finally, the failed material is removed and the mesh is redistributed accordingly. The outcome is a model capable of predicting the total recession in the material due to surface chemistry and mechanical erosion.

Mechanical Properties↗

Cell model and elastic moduli of disordered solids - Low temperature limit

The cell theory has been previously employed to compute the equation of state of a disordered condensed system. It is now generalized to include anisotropic stresses. The condition of affine deformation is adopted, transforming an orginally spherical into an ellipsoidal cell. With a Lennard-Jones n-m potential between nonbonded centers, the formal expression for the deformational free energy is derived. It is to be evaluated in the limit of the linear elastic range. Since the bulk modulus in this limit is already known, it is convenient to consider a uniaxial deformation. To begin with, restrictions are made to the low-temperature limit in the absence of entropy contributions. Young's modulus and Poisson's ratio then follow.

Peng, S. T. J.↗

Wave propagation in anisotropic medium due to an oscillatory point source with application to unidirectional composites

The far-field displacements in an infinite transversely isotropic elastic medium subjected to an oscillatory concentrated force are derived. The concepts of velocity surface, slowness surface and wave surface are used to describe the geometry of the wave propagation process. It is shown that the decay of the wave amplitudes depends not only on the distance from the source (as in isotropic media) but also depends on the direction of the point of interest from the source. As an example, the displacement field is computed for a laboratory fabricated unidirectional fiberglass epoxy composite. The solution for the displacements is expressed as an amplitude distribution and is presented in polar diagrams. This analysis has potential usefulness in the acoustic emission (AE) and ultrasonic nondestructive evaluation of composite materials. For example, the transient localized disturbances which are generally associated with AE sources can be modeled via this analysis. In which case, knowledge of the displacement field which arrives at a receiving transducer allows inferences regarding the strength and orientation of the source, and consequently perhaps the degree of damage within the composite.

Williams, J. H., Jr.↗

Composites with periodic microstructure

For an elastic body containing periodically distributed inhomogeneities, a general procedure is developed for estimating the overall properties of the composite in terms of several infinite series which, for the isotropic matrix (but anisotropic inclusions), depend only on the geometry of the inhomogeneities and hence can be calculated once and for all for each geometry. These infinite series are obtained and tabulated for ellipsoidal inhomogeneities, and the results are used to estimate the overall elastic moduli of composites which contain spherical or ellipsoidal voids or elastic inclusions.

Iwakuma, T.↗

Three-dimensional solutions for the free vibrations and buckling of thermally stressed multilayered angle-ply composite plates

Analytic three-dimensional elasticity solutions are developed for the free vibration and buckling of thermally stressed rectangular multilayered angle-ply anisotropic plates which are assumed to have an antisymmetric lamination with respect to the middle plane. Sensitivity derivatives are evaluated and used to investigate the sensitivity of the vibration and buckling responses to variations in the different lamination and material parameters of the plate. A Duhamel-Neumann-type constitutive model is used, and the material properties are assumed to be independent of temperature. Numerical results are presented, showing the effects of variations in the material characteristics and fiber orientation of different layers, as well as the effect of initial thermal deformation on the vibrational and buckling responses of the plate.

Noor, Ahmed K.↗

Anisotropic Plasticity of BN Nanotubes

Plastic collapse of compressed BN nanotubes are investigated and compared with carbon nanotubes of similar nature. Using a generalized tight-binding molecular dynamics (GTBMD) method for system containing B, N and C atoms we compute stiffness and plastic collapse of BN and C nanotubes under axial compression. For small compressional strain, BN nanotubes are found to be about 92% as stiff as similar C nanotubes. Due to BN bond buckling effect, however, the elastic limit of BN nanotubes is found to be more than C nanotubes. A route to plasticity is explored in which we find that at elastic limit the accumulated strain is released by a local plastic deformation of the nanotube. The mechanism of strain release and the resulting plastic deformation, however. are anisotropic in nature. The strain is released preferentially towards N as leading edge of a buckled BN bond and the tube, compressed at both ends, plastically collapses preferentially towards one end. Details of the anisotropic plasticity and prospective applications will be discussed in this presentation.

Madhu, Menon↗

An improved boundary force method for analyzing cracked anisotropic materials

The Boundary Force Method (BFM), a form of indirect boundary element method, is used to analyze composite laminates with cracks. The BFM uses the orthotropic elasticity solution for a concentrated horizontal and vertical force and a moment applied at a point in a cracked, infinite sheet as the fundamental solution. The necessary stress functions for this fundamental solution were formulated using the complex variable theory of orthotropic elasticity. The current method is an improvement over a previous method using only forces and no moment. The improved method was verified by comparing it to accepted solutions for a finite-width, center-crack specimen subjected to uniaxial tension. Four graphite/epoxy laminates were used: (0 + or - 45/90)sub s, (0), (+ or - 45)sub s, and (+ or - 30)sub s. The BFM results agreed well with accepted solutions. Convergence studies showed that with the addition of the moment in the fundamental solution, the number of boundary elements required for a converged solution was significantly reduced. Parametric studies were done for two configurations for which no orthotropic solutions are currently available; a single edge crack and an inclined single edge crack.

Tan, Paul W.↗

An improved boundary force method for analysing cracked anisotropic materials

The Boundary Force Method (BFM), a form of indirect boundary element method, is used to analyze composite laminates with cracks. The BFM uses the orthotropic elasticity solution for a concentrated horizontal and vertical force and a moment applied at a point in a cracked, infinite sheet as the fundamental solution. The necessary stress functions for this fundamental solution were formulated using the complex variables theory of orthotropic elasticity. The current method is an improvement over a previous method using only forces and no moment. The improved method was verified by comparing it to accepted solutions for a finite-width, center-crack specimen subjected to uniaxial tension. Four graphite/epoxy laminates were used: (0 + or - 45/90)sub s, (0), (+ or - 45)sub s, and (+ or - 30)sub s. The BFM results agreed well with accepted solutions. Convergence studies showed that with the addition of the moment in the fundamental solution, the number of boundary elements required for a converged solution was significantly reduced. Parametric studies were done for two configurations for which no orthotropic solutions are currently available; a single edge crack and an inclined single edge crack.

Tan, P. W.↗

The penetration of soft electrons into the ionosphere

Calculations are presented of energy spectra and angular and spatial distributions of electron fluxes in the ionosphere resulting from precipitation of monoenergetic (E = 25, 50, and 100 eV) electrons. The incident electrons are assumed to be isotropic over the downward direction. It is found that the resulting steady-state electron fluxes above ca. 300 km are highly anisotropic and that the pitch-angle distribution is energy dependent. About 15% of the incident electrons are backscattered elastically to the protonosphere. A much larger number of electrons escapes after they have deposited part of their energy in the atmosphere. The mean energy of the escaping electrons is about half that of the incident electrons. About 50% of the incident energy is absorbed in the atmosphere, the remainder being returned to the protonosphere. The rate of energy absorption is a maximum at heights between 300 and 400 km. Most of the energy is absorbed in ionization and excitation of atomic oxygen; an appreciable amount is absorbed as heat by the ambient electron gas. Altitude profiles are presented for the rates of ionization, excitation, and electron heating caused by soft-electron precipitation.

Mantas, G. P.↗

Procedure for implementation of temperature-dependent mechanical property capability in the Engineering Analysis Language (EAL) system

A procedure is presented to allow the use of temperature dependent mechanical properties in the Engineering Analysis Language (EAL) System for solid structural elements. This is accomplished by including a modular runstream in the main EAL runstream. The procedure is applicable for models with multiple materials and with anisotropic properties, and can easily be incorporated into an existing EAL runstream. The procedure (which is applicable for EAL elastic solid elements) is described in detail, followed by a description of the validation of the routine. A listing of the EAL runstream used to validate the procedure is included in the Appendix.

Glass, David E.↗

Three-Dimensional Stress Fields and Slip Systems for Single Crystal Superalloy Notched Specimens

Single crystal superalloys have become increasingly popular for turbine blade and vane applications due to their high strength, and creep and fatigue resistance at elevated temperatures. The crystallographic orientation of a single crystal material greatly affects its material properties, including elastic modulus, shear modulus, and ductility. These directional properties, along with the type of loading and temperature, dictate an anisotropic response in the yield strength, creep resistance, creep rupture ductility, fatigue resistance, etc. A significant amount of research has been conducted to determine the material properties in the <001> orientation, yet the material properties deviating from the <001> orientation have not been assessed for all cases. Based on the desired application and design criteria, a crystal orientation is selected to yield the maximum properties. Currently, single crystal manufacturing is able to control the primary crystallographic orientation within 15 of the target orientation, which is an acceptable deviation to meet both performance and cost guidelines; the secondary orientation is rarely specified. A common experiment is the standard load-controlled tensile test, in which specimens with different orientations can be loaded to observe the material response. The deformation behavior of single-crystal materials under tension and compression is known to be a function of not only material orientation, but also of varying microdeformation (i.e. dislocation) mechanisms. The underlying dislocation motion causes deformation via slip, and affects the activation of specific slip systems based on load and orientation. The slip can be analyzed by observing the visible traces left on the surface of the specimen from the slip activity within the single crystal material. The goal of this thesis was to predict the slip systems activated in three-dimensional stress fields of a notched tensile specimen, as a function of crystal orientation, using finite element analysis without addressing microstructural deformation mechanisms that govern their activation. Out of three orientations tested, the specimen with a [110] load orientation and a [001] growth direction had the lowest maximum resolved shear stress; this specimen orientation appears to be the best design candidate for a tensile application.

Magnan, Shannon M.↗