Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Two-dimensional material”

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 235 records · Page 13

Edge delamination of composite laminates subject to combined tension and torsional loading

Delamination is a common failure mode of laminated composite materials. Edge delamination is important since it results in reduced stiffness and strength of the laminate. The tension/torsion load condition is of particular significance to the structural integrity of composite helicopter rotor systems. Material coupons can easily be tested under this type of loading in servo-hydraulic tension/torsion test stands using techniques very similar to those used for the Edge Delamination Tensile Test (EDT) delamination specimen. Edge delamination of specimens loaded in tension was successfully analyzed by several investigators using both classical laminate theory and quasi-three dimensional (Q3D) finite element techniques. The former analysis technique can be used to predict the total strain energy release rate, while the latter technique enables the calculation of the mixed-mode strain energy release rates. The Q3D analysis is very efficient since it produces a three-dimensional solution to a two-dimensional domain. A computer program was developed which generates PATRAN commands to generate the finite element model. PATRAN is a pre- and post-processor which is commonly used with a variety of finite element programs such as MCS/NASTRAN. The program creates a sufficiently dense mesh at the delamination crack tips to support a mixed-mode fracture mechanics analysis. The program creates a coarse mesh in those regions where the gradients in the stress field are low (away from the delamination regions). A transition mesh is defined between these regions. This program is capable of generating a mesh for an arbitrarily oriented matrix crack. This program significantly reduces the modeling time required to generate these finite element meshes, thus providing a realistic tool with which to investigate the tension torsion problem.

Hooper, Steven J.↗

MBE HgCdTe heterostructure detectors

HgCdTe has been the mainstay for medium (3 to 5 micron) and long (10 to 14 micron) wavelength infrared detectors in recent years. Conventional growth and processing techniques are continuing to improve the material. However, the additional ability to tailor composition and placement of doped layers on the tens of angstroms scale using molecular beam epitaxy (MBE) provides the opportunity for new device physics and concepts to be utilized. MBE-based device structures to be discussed here can be grouped into two categories: tailored conventional structures and quantum structures. The tailored conventional structures are improvements on familiar devices, but make use of the ability to create layers of varying composition, and thus band gap, at will. The heterostructure junction can be positioned independently of doping p-n junctions. This allows the small band gap region in which the absorption occurs to be separated from a larger band gap region in which the electric field is large and where unwanted tunneling can occur. Data from hybrid MBE/liquid phase epitaxy (LPE)/bulk structures are given. Quantum structures include the HgTe-CdTe superlattice, in which the band gap and transport can be controlled by alternating thin layers (tens of angstroms thick) of HgTe and CdTe. The superlattice has been shown to exhibit behavior which is non-alloy like, including very high hole mobilities, two-dimensional structure in the absorption coefficient, resonant tunneling, and anisotropic transport.

Schulman, Joel N.↗

Imaging of molecular hydrogen and oxygen by single and two-photon fluorescence using laser and flashlamp sources

Two flow visualization techniques, i.e., simultaneous two-dimensional fluorescence imaging of H2 and O2 in a diffusion flame, and quasi-linear fluorescence imaging of O2, are presented. The first uses an injection-locked argon-fluoride excimer laser and a partial overlap of a two-photon ground state absorption in H2 with a single photon absorption from a vibrational level in O2. The second uses a simple, high-intensity ultraviolet flashlamp which provides a flux of photons in the 180-195 nm range, sufficient to produce a quasi-one-dimensional fluorescence image of hot/room temperature oxygen. Both techniques do not require that a seed material be introduced into the flow, they can image major flow constituents, and provide an instantaneous snapshot of the flow.

Diskin, Glenn S.↗

Displaying Data As Movies

NMSB Movie computer program displays large sets of data (more than million individual values). Presentation dynamic, rapidly displaying sequential image "frames" in main "movie" window. Any sequence of two-dimensional sets of data scaled between 0 and 255 (1-byte resolution) displayed as movie. Time- or slice-wise progression of data illustrated. Originally written to present data from three-dimensional ultrasonic scans of damaged aerospace composite materials, illustrates data acquired by thermal-analysis systems measuring rates of heating and cooling of various materials. Developed on Macintosh IIx computer with 8-bit color display adapter and 8 megabytes of memory using Symantec Corporation's Think C, version 4.0.

Moore, Judith G.↗

Three-dimensional elastic-plastic finite-element analyses of constraint variations in cracked bodies

Three-dimensional elastic-plastic (small-strain) finite-element analyses were used to study the stresses, deformations, and constraint variations around a straight-through crack in finite-thickness plates for an elastic-perfectly plastic material under monotonic and cyclic loading. Middle-crack tension specimens were analyzed for thicknesses ranging from 1.25 to 20 mm with various crack lengths. Three local constraint parameters, related to the normal, tangential, and hydrostatic stresses, showed similar variations along the crack front for a given thickness and applied stress level. Numerical analyses indicated that cyclic stress history and crack growth reduced the local constraint parameters in the interior of a plate, especially at high applied stress levels. A global constraint factor alpha(sub g) was defined to simulate three-dimensional effects in two-dimensional crack analyses. The global constraint factor was calculated as an average through-the-thickness value over the crack-front plastic region. Values of alpha(sub g) were found to be nearly independent of crack length and were related to the stress-intensity factor for a given thickness.

Newman, J. C., Jr.↗

Design optimization of space structures

The topology-shape-size optimization of space structures is investigated through Kikuchi's homogenization method. The method starts from a 'design domain block,' which is a region of space into which the structure is to materialize. This domain is initially filled with a finite element mesh, typically regular. Force and displacement boundary conditions corresponding to applied loads and supports are applied at specific points in the domain. An optimal structure is to be 'carved out' of the design under two conditions: (1) a cost function is to be minimized, and (2) equality or inequality constraints are to be satisfied. The 'carving' process is accomplished by letting microstructure holes develop and grow in elements during the optimization process. These holes have a rectangular shape in two dimensions and a cubical shape in three dimensions, and may also rotate with respect to the reference axes. The properties of the perforated element are obtained through an homogenization procedure. Once a hole reaches the volume of the element, that element effectively disappears. The project has two phases. In the first phase the method was implemented as the combination of two computer programs: a finite element module, and an optimization driver. In the second part, focus is on the application of this technique to planetary structures. The finite element part of the method was programmed for the two-dimensional case using four-node quadrilateral elements to cover the design domain. An element homogenization technique different from that of Kikuchi and coworkers was implemented. The optimization driver is based on an augmented Lagrangian optimizer, with the volume constraint treated as a Courant penalty function. The optimizer has to be especially tuned to this type of optimization because the number of design variables can reach into the thousands. The driver is presently under development.

Felippa, Carlos↗

Interlaminar stress analysis at the skin-stiffener interface of a grid-stiffened composite panel

Results are presented of an analysis of the interlaminar stress components at the skin-wrap, skin-core, and wrap-core interfaces for an advanced-concept stiffened panel using a global/local finite element procedure. The procedure consists of a global model of two-dimensional shell elements that is used to design a grid-stiffened panel with blade-type stiffeners, a local model of three-dimensional solid elements that is used to compute interlaminar stress components, and a scheme devised to assign displacement boundary conditions for the local model that are derived from displacement and rotation data of a few nodes of the global model. The interlaminar stresses for this panel design were found to be well below typical tensile normal and shearing strengths of a graphite-epoxy material.

Wiggenraad, J. F. M.↗

Laplace's equation and the Dirichlet-Neumann map in multiply connected domains

A variety of problems in material science and fluid dynamics require the solution of Laplace's equation in multiply connected domains. Integral equation methods are natural candidates for such problems, since they discretize the boundary alone, require no special effort for free boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains in two space dimensions, regardless of shape. Current integral equation methods for the Dirichlet problem, however, require the solution of M independent problems of dimension N, where M is the number of boundary components and N is the total number of points in the discretization. In this paper, we present a new boundary integral equation approach, valid for both interior and exterior problems, which requires the solution of a single linear system of dimension N + M. We solve this system by making use of an iterative method (GMRES) combined with the last multipole method for the rapid calculation of the necessary matrix vector products. For a two-dimensional system with 200 components and 100 points on each boundary, we gain a speedup of a factor of 100 from the new analytic formulation and a factor of 50 from the fast multipole method. The resulting scheme brings large scale calculations in extremely complex domains within practical reach.

Greenbaum, A.↗

Strain energy of thermally stressed multilayer panels and its sensitivity coefficients

A study is made of the effect of lamination and material parameters of thermally stressed multilayer composite panels on the strain energy and its sensitivity coefficients obtained by the three-dimensional thermoelasticity model. The sensitivity coefficients measure the sensitivity of the strain energy to variations in the different lamination and material parameters of the panel. The strain energy and sensitivity coefficients obtained by the three-dimensional model are used as the basis for assessing the accuracy of the corresponding quantities obtained by two-dimensional first-order shear deformation theory. Numerical results are presented for antisymmetrically laminated angle-ply panels subjected to linear temperature variation through the thickness. The numerical results demonstrate the accuracy of the strain energy and its sensitivity coefficients obtained by the first-order shear deformation theory, for the range of parameters considered, and the importance of the particular terms in the strain energy, which vanish for single-layer panels.

Noor, Ahmed K.↗

Modal element method for scattering and absorbing of sound by two-dimensional bodies

The modal element method for acoustic scattering from a 2-D body is presented. The body may be acoustically soft (absorbing) or hard (reflecting). The infinite computational region is divided into two subdomains - the bounded finite element domain, which is characterized by complicated geometry and/or variable material properties, and the surrounding unbounded homogeneous domain. The acoustic pressure field is represented approximately in the finite element domain by a finite element solution, and is represented analytically by an eigenfunction expansion in the homogeneous domain. The two solutions are coupled by the continuity of pressure and velocity across the interface between the two subdomains. Also, for hard bodies, a compact modal ring grid system is introduced for which computing requirements are drastically reduced. Analysis for 2-D scattering from solid and coated (acoustically treated) bodies is presented, and several simple numerical examples are discussed. In addition, criteria are presented for determining the number of modes to accurately resolve the scattered pressure field from a solid cylinder as a function of the frequency of the incoming wave and the radius of the cylinder.

Baumeister, K. J.↗

Theoretical investigation of flutter of two-dimensional flat panels with one surface exposed to supersonic potential flow

A Rayleigh type analysis involving chosen modes of the panel as degrees of freedom is used to treat the flutter of a two-dimensional flat panel supported at its leading and trailing edges and subjected to a middle-plane tensile force. The panel has a supersonic stream passing over its upper surface and still air below. The aerodynamic forces due to the supersonic stream are obtained from the theory for linearized two-dimensional unsteady flow and the forces due to the still air are obtained from acoustical theory. In order to study the effect of increasing the number of modes in the analysis, two and then four modes are employed. The modes used are the first four natural modes of the panel in a vacuum with no tensile force acting. The analysis includes these variables: Mach number, structural damping, tensile force, density of the still air, and edge fixity (clamped and pinned). For certain combinations of these variables, stability boundaries are obtained which can be used to determine the panel thickness required to prevent flutter for any panel material and altitude.

Nelson, Herbert C↗

Failure prediction of thin beryllium sheets used in spacecraft structures

The primary objective of this study is to develop a method for prediction of failure of thin beryllium sheets that undergo complex states of stress. Major components of the research include experimental evaluation of strength parameters for cross-rolled beryllium sheet, application of the Tsai-Wu failure criterion to plate bending problems, development of a high order failure criterion, application of the new criterion to a variety of structures, and incorporation of both failure criteria into a finite element code. A Tsai-Wu failure model for SR-200 sheet material is developed from available tensile data, experiments carried out by NASA on two circular plates, and compression and off-axis experiments performed in this study. The failure surface obtained from the resulting criterion forms an ellipsoid. By supplementing experimental data used in the the two-dimensional criterion and modifying previously suggested failure criteria, a multi-dimensional failure surface is proposed for thin beryllium structures. The new criterion for orthotropic material is represented by a failure surface in six-dimensional stress space. In order to determine coefficients of the governing equation, a number of uniaxial, biaxial, and triaxial experiments are required. Details of these experiments and a complementary ultrasonic investigation are described in detail. Finally, validity of the criterion and newly determined mechanical properties is established through experiments on structures composed of SR200 sheet material. These experiments include a plate-plug arrangement under a complex state of stress and a series of plates with an out-of-plane central point load. Both criteria have been incorporated into a general purpose finite element analysis code. Numerical simulation incrementally applied loads to a structural component that is being designed and checks each nodal point in the model for exceedance of a failure criterion. If stresses at all locations do not exceed the failure criterion, the load is increased and the process is repeated. Failure results for the plate-plug and clamped plate tests are accurate to within 2 percent.

Roschke, Paul N.↗

Development of a portable multispectral thermal infrared camera

The purpose of this research and development effort was to design and build a prototype instrument designated the 'Thermal Infrared Multispectral Camera' (TIRC). The Phase 2 effort was a continuation of the Phase 1 feasibility study and preliminary design for such an instrument. The completed instrument designated AA465 has application in the field of geologic remote sensing and exploration. The AA465 Thermal Infrared Camera (TIRC) System is a field-portable multispectral thermal infrared camera operating over the 8.0 - 13.0 micron wavelength range. Its primary function is to acquire two-dimensional thermal infrared images of user-selected scenes. Thermal infrared energy emitted by the scene is collected, dispersed into ten 0.5 micron wide channels, and then measured and recorded by the AA465 System. This multispectral information is presented in real time on a color display to be used by the operator to identify spectral and spatial variations in the scenes emissivity and/or irradiance. This fundamental instrument capability has a wide variety of commercial and research applications. While ideally suited for two-man operation in the field, the AA465 System can be transported and operated effectively by a single user. Functionally, the instrument operates as if it were a single exposure camera. System measurement sensitivity requirements dictate relatively long (several minutes) instrument exposure times. As such, the instrument is not suited for recording time-variant information. The AA465 was fabricated, assembled, tested, and documented during this Phase 2 work period. The detailed design and fabrication of the instrument was performed during the period of June 1989 to July 1990. The software development effort and instrument integration/test extended from July 1990 to February 1991. Software development included an operator interface/menu structure, instrument internal control functions, DSP image processing code, and a display algorithm coding program. The instrument was delivered to NASA in March 1991. Potential commercial and research uses for this instrument are in its primary application as a field geologists exploration tool. Other applications have been suggested but not investigated in depth. These are measurements of process control in commercial materials processing and quality control functions which require information on surface heterogeneity.

Osterwisch, Frederick G.↗

Variations of a global constraint factor in cracked bodies under tension and bending loads

Elastic-plastic finite-element analyses were used to calculate stresses and displacements around a crack in finite-thickness plates for an elastic-perfectly plastic material. Middle- and edge-crack specimens were analyzed under tension and bending loads. Specimens were 1.25 to 20 mm thick with various widths and crack lengths. A global constraint factor alpha(sub g), an averaged normal-stress to flow-stress ratio over the plastic region, was defined to simulate three-dimensional (3D) effects in two-dimensional (2D) models. For crack lengths and uncracked ligament lengths greater than four times the thickness, the global constraint factor was found to be nearly a unique function of a normalized stress-intensity factor (related to plastic-zone size to thickness ratio) from small- to large-scale yielding conditions for various specimen types and thickness. For crack length-to-thickness ratios less than four, the global constraint factor was specimen type, crack length and thickness dependent. Using a 2D strip-yield model and the global constraint factors, plastic-zone sizes and crack-tip displacements agreed reasonably well with the 3D analyses. For a thin sheet aluminum alloy, the critical crack-tip-opening angle during stable tearing was found to be independent of specimen type and crack length for crack length-to-thickness ratios greater than 4.

Newman, J. C., Jr.↗

Optimal experimental designs for the estimation of thermal properties of composite materials

Reliable estimation of thermal properties is extremely important in the utilization of new advanced materials, such as composite materials. The accuracy of these estimates can be increased if the experiments are designed carefully. The objectives of this study are to design optimal experiments to be used in the prediction of these thermal properties and to then utilize these designs in the development of an estimation procedure to determine the effective thermal properties (thermal conductivity and volumetric heat capacity). The experiments were optimized by choosing experimental parameters that maximize the temperature derivatives with respect to all of the unknown thermal properties. This procedure has the effect of minimizing the confidence intervals of the resulting thermal property estimates. Both one-dimensional and two-dimensional experimental designs were optimized. A heat flux boundary condition is required in both analyses for the simultaneous estimation of the thermal properties. For the one-dimensional experiment, the parameters optimized were the heating time of the applied heat flux, the temperature sensor location, and the experimental time. In addition to these parameters, the optimal location of the heat flux was also determined for the two-dimensional experiments. Utilizing the optimal one-dimensional experiment, the effective thermal conductivity perpendicular to the fibers and the effective volumetric heat capacity were then estimated for an IM7-Bismaleimide composite material. The estimation procedure used is based on the minimization of a least squares function which incorporates both calculated and measured temperatures and allows for the parameters to be estimated simultaneously.

Scott, Elaine P.↗

On the dynamics of a shock-bubble interaction

We present a detailed numerical study of the interaction of a weak shock wave with an isolated cylindrical gas inhomogenity. Such interactions have been studied experimentally in an attempt to elucidate the mechanisms whereby shock waves propagating through random media enhance mixing. Our study concentrates on the early phases of the interaction process which are dominated by repeated refractions of acoustic fronts at the bubble interface. Specifically, we have reproduced two of the experiments performed by Haas and Sturtevant : M(sub s) = 1.22 planar shock wave, moving through air, impinges on a cylindrical bubble which contains either helium or Refrigerant 22. These flows are modelled using the two-dimensional, compressible Euler equations for a two component fluid (air-helium or air-Refrigerant 22). Although simulations of shock wave phenomena are now fairly commonplace, they are mostly restricted to single component flows. Unfortunately, multi-component extensions of successful single component schemes often suffer from spurious oscillations which are generated at material interfaces. Here we avoid such problems by employing a novel, nonconservative shock-capturing scheme. In addition, we have utilized a sophisticated adaptive mesh refinement algorithm which enables extremely high resolution simulations to be performed relatively cheaply. Thus we have been able to reproduce numerically all the intricate mechanisms that were observed experimentally (e.g., transitions from regular to irregular refraction, cusp formation and shock wave focusing, multi-shock and Mach shock structures, jet formation, etc.), and we can now present an updated description for the dynamics of a shock-bubble interaction.

Quirk, James J.↗

Effect of the fiber-matrix interphase on the transverse tensile strength of the unidirectional composite material

A simple one-dimensional fiber-matrix interphase model has been developed and analytical results obtained correlated well with available experimental data. It was found that by including the interphase between the fiber and matrix in the model, much better local stress results were obtained than with the model without the interphase. A more sophisticated two-dimensional micromechanical model, which included the interphase properties was also developed. Both one-dimensional and two-dimensional models were used to study the effect of the interphase properties on the local stresses at the fiber, interphase and matrix. From this study, it was found that interphase modulus and thickness have significant influence on the transverse tensile strength and mode of failure in fiber reinforced composites.

Tsai, H. C.↗

Simulations of Titan's Brightness by a Two-Dimensional Haze Model

We have used a 2-D microphysics model to study the effects of atmospheric motions on the albedo of Titan's thick haze layer. We compare our results to the observed variations of Titan's brightness with season and latitude. We use two wind fields. The first is an analytic solution to the simplified equations of motion. The second is based on the preliminary results of a Titan GCM. Seasonally varying wind fields, of typical horizontal velocity 1 cm/ s at optical depth unity, are capable of producing the observed change in geometric albedo of about 10% over the Titan year. Neither of the two wind fields can adequately reproduce the latitudinal distribution of the contrast seen by Voyager. At visible wavelengths, where only haze opacity is radiatively important, upwelling produces darkening by increasing the particle size at optical depth unity. This is due to the suspension of larger particles as well as the removal of smaller particles from the top of the atmosphere. From our results it can be inferred that the circulation on Titan just prior to the time of Voyager was upwelling in the northern hemisphere and downwelling in the southern hemisphere with broad regions of uniform wind separated by a sharp transition at 10 deg. S. At UV wavelengths and at 0.89 microns the albedo is determined by the competing effects of the gas and the haze material. Gas is bright in the UV and dark at 0.89 microns. Haze transport at high altitudes controls the UV albedo and transport at low altitude controls the 0.89 microns albedo. Comparisons between the hemispheric contrast at UV and IR wavelengths can be diagnostic of the vertical structure of the wind field on Titan. Infrared albedo features, which are in the lower atmosphere and persist for many rotations could reflect the haze distribution set up by a seasonally varying circulation rather than surface features as assumed in some observational studies.

Hutzell, W.↗