Search NASA⌕ Search

SEARCH · Search NASA

Results for “method development”

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 649 records · Page 36

Solar wind disturbances caused by solar flares - Equatorial plane

The propagation of solar wind disturbances caused by single, double and six successive flares in the dipolar and quadrupolar patterns of the interplanetary magnetic field (IMF) and the associated solar wind flow is studied. This study is based on a kinematic and empirical method developed by Hakamada and Akasofu (1982). Each flare is characterized by six parameters (such as the highest speed flow, its extent and duration). The successive IMF patterns in the equatorial plane of the heliosphere during a time span of 0.5-60 days after flares are presented for a variety of flares. The solar wind speed and IMF magnitude are given as a function of distance along a radial line fixed in space and also as a function of time at several points fixed in space (simulating approximately space probe observations). Some of the results are qualitatively compared with recent space probe observations, demonstrating fair similarity with the observed time profiles of solar wind speed variations over a wide range of both distances (0-10 AU) and time spans (60 days). The method provides a first-order construction, temporal and spatial, of flare-induced shocks and their multiple interactions with each other, as well as with the corotating interaction regions.

Akasofu, S.-I.↗

Thermal verification testing of commercial printed-circuit boards for spaceflight

A method developed to verify commercial printed-circuit boards for a Shuttle orbital flight is discussed. The test sequence is based on early fault detection, desire to test the final assembly, and integration with other verification testing. A component thermal screening test is performed first to force flaws in design, workmanship, parts, processes, and materials into observable failures. Temperature definition and vibration tests are performed next. Final assembly testing is performed to simulate the Shuttle flight. An abbreviated thermal screening test is performed as a check after the vibration test, and then a complete thermal operational test is performed. The final assembly test finishes up with a burn-in of 100 h of trouble-free operation. Verification is successful when all components and final assemblies have passed each test. This method was very successful in verifying that commercial printed-circuit boards will survive in the Shuttle environment.

Foster, William M., II↗

Production of a long-term global water vapor and liquid water data set using ultra-fast methods to assimilate multi-satellite and radiosonde observations

There is a well-documented requirement for a comprehensive and accurate global moisture data set to assist many important studies in atmospheric science. Currently, atmospheric water vapor measurements are made from a variety of sources including radiosondes, aircraft and surface observations, and in recent years, by various satellite instruments. Creating a global data set from a single measuring system produces results that are useful and accurate only in specific situations and/or areas. Therefore, an accurate global moisture data set has been derived from a combination of these measurement systems. Under a NASA peer-reviewed contract, STC-METSAT produced two 5-yr (1988-1992) global data sets. One is the total column (integrated) water vapor data set and the other, a global layered water vapor data set using a combination of radiosonde observations, Television and Infrared Observation Satellite (TIROS) Operational Satellite (TOVS), and Special Sensor Microwave/Imager (SSM/I) data sets. STC-METSAT also produced a companion, global, integrated liquid water data set. The complete data set (all three products) has been named NVAP, an anachronym for NASA Water Vapor Project. STC-METSAT developed methods to process the data at a daily time scale and 1 x 1 deg spatial resolution.

Vonderhaar, Thomas H.↗

Simulation-Based Height of Burst Map for Asteroid Airburst Damage Prediction

Entry and breakup models predict that airburst in the Earth's atmosphere is likely for asteroids up to approximately 200 meters in diameter. Objects of this size can deposit over 250 megatons of energy into the atmosphere. Fast-running ground damage prediction codes for such events rely heavily upon methods developed from nuclear weapons research to estimate the damage potential for an airburst at altitude. (Collins, 2005; Mathias, 2017; Hills and Goda, 1993). In particular, these tools rely upon the powerful yield scaling laws developed for point-source blasts that are used in conjunction with a Height of Burst (HOB) map to predict ground damage for an airburst of a specific energy at a given altitude. While this approach works extremely well for yields as large as tens of megatons, it becomes less accurate as yields increase to the hundreds of megatons potentially released by larger airburst events. This study revisits the assumptions underlying this approach and shows how atmospheric buoyancy becomes important as yield increases beyond a few megatons. We then use large-scale three-dimensional simulations to construct numerically generated height of burst maps that are appropriate at the higher energy levels associated with the entry of asteroids with diameters of hundreds of meters. These numerically generated HOB maps can then be incorporated into engineering methods for damage prediction, significantly improving their accuracy for asteroids with diameters greater than 80-100 m.

Burst Map↗

Co-leveraging Scientific Advances in Space Biology and Astrobiology Towards Achieving NASA’s Life Science Objectives

Executive Summary: Distinct lines of scientific inquiry drives the separation of NASA’s fundamental life science research into Space Biology and Astrobiology. This division developed as a way to place life scientists alongside experts in the physical constraints that define the acclimation, adaptation and evolution of biology systems relevant to their respective subjects. For astrobiology, integration with disciplines such as geology, geochemistry, astronomy, planetary science, etc., enables a comprehensive assessment of the physical environment and its co-evolution with biological processes. Space Biology’s co-location with Physical Sciences places life science researchers adjacent to experts in the physical phenomena associated with microgravity and spaceflight, enabling an understanding of how the spaceflight environment affects biological systems. Despite this separation, aspects of both disciplines have converged on a similar, fundamental objective: to describe and understand the dynamics of complex living communities in the contexts of their physical environments. While the environmental systems and timescales are dramatically different, continuing to motivate the separation into distinct fields, similarities in the underlying objective present opportunities to find efficiencies, reduce overlap, and minimize duplication of effort. Space Biology and Astrobiology share a common need to understand microbial physiology in extreme environments – whether the ‘built’ spaceflight environment or the natural environments in which many astrobiology studies are conducted. In particular, open questions in each discipline require the development of quantitative frameworks, applicable at the ecosystem level, that support predictive capabilities for environments where observations are sparse. Additionally, both disciplines have a need to prepare, detect, and analyze the (potential) biological signal in complex samples-often in a completely autonomous fashion. The next decade will see NASA Space Biology moving to understand and describe the effects of the beyond low-earth orbit (BLEO) spaceflight environment on living systems. This new direction will dramatically reduce the opportunities for ground-based analysis of space-flown samples, driving space biology investigations towards fully autonomous experiments and missions. At the same time, astrobiology life detection missions aimed at detecting biosignatures on Mars and icy moons in the outer solar system could benefit from fully automated sample processing and analysis. There are opportunities to leverage instrument and method development between both disciplines within the context of these BLEO missions.

Astrobiology↗

Energy Management of Manned Boost-Glide Vehicles: A Historical Perspective

As flight progressed from propellers to jets to rockets, the propulsive energy grew exponentially. With the development of rocket-only boosted vehicles, energy management of these boost-gliders became a distinct requirement for the unpowered return to base, alternate landing site, or water-parachute landing, starting with the X-series rocket aircraft and terminating with the present-day Shuttle. The problem presented here consists of: speed (kinetic energy) - altitude (potential energy) - steep glide angles created by low lift-to-drag ratios (L/D) - distance to landing site - and the bothersome effects of the atmospheric characteristics varying with altitude. The primary discussion regards post-boost, stabilized glides; however, the effects of centrifugal and geopotential acceleration are discussed as well. The aircraft and spacecraft discussed here are the X-1, X-2, X-15, and the Shuttle; and to a lesser, comparative extent, Mercury, Gemini, Apollo, and lifting bodies. The footprints, landfalls, and methods developed for energy management are also described. The essential tools required for energy management - simulator planning, instrumentation, radar, telemetry, extended land or water range, Mission Control Center (with specialist controllers), and emergency alternate landing sites - were first established through development of early concepts and were then validated by research flight tests.

Day, Richard E.↗

Detecting Small-Scale Topographic Changes and Relict Geomorphic Features on Barrier Islands Using SAR

The shapes and elevations of barrier islands may change dramatically over a short period of time during a storm. Coastal scientists and engineers, however, are currently unable to measure these changes occurring over an entire barrier island at once. This three-year project, which is funded by NASA and jointly conducted by the Bureau of Economic Geology and the Center for Space Research at The University of Texas at Austin, is designed to overcome this problem by developing the use of interferometry from airborne synthetic aperture radar (AIRSAR) to measure coastal topography and to detect storm-induced changes in topography. Surrogate measures of topography observed in multiband, fully polarimetric AIRSAR (This type of data are now referred to as POLSAR data.) are also being investigated. Digital elevation models (DEM) of Galveston Island and Bolivar Peninsula, Texas obtained with Topographic SAR (TOPSAR) are compared with measurements by Global Positioning System (GPS) ground surveys and electronic total station surveys. In addition to topographic mapping, this project is evaluating the use of POLSAR to detect old features such as storm scarps, storm channels, former tidal inlets, and beach ridges that have been obscured by vegetation, erosion, deposition, and artificial filling. We have also expanded the work from the original proposal to include the mapping of coastal wetland vegetation and depositional environments. Methods developed during this project will provide coastal geologists with an unprecedented tool for monitoring and understanding barrier island systems. This understanding will improve overall coastal management policies and will help reduce the effects of natural and man-induced coastal hazards. This report summarizes our accomplishments during the second year of the study. Also included is a discussion of our planned activities for year 3 and a revised budget.

Gibeaut, James C.↗

Navier-Stokes simulations of unsteady transonic flow phenomena

Numerical simulations of two classes of unsteady flows are obtained via the Navier-Stokes equations: a blast-wave/target interaction problem class and a transonic cavity flow problem class. The method developed for the viscous blast-wave/target interaction problem assumes a laminar, perfect gas implemented in a structured finite-volume framework. The approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the full viscous terms are computed by central differencing. Comparisons of unsteady numerical, analytical, and experimental results are made in two- and three-dimensions for Couette flows, a starting shock-tunnel, and a shock-tube blockage study. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing of the predominantly inviscid flows. Viscous effects were increasingly significant at large post-interaction times. While the blast-wave/target interaction problem benefits from high-resolution methods applied to the Euler terms, the transonic cavity flow problem requires the use of an efficient scheme implemented in a geometrically flexible overset mesh environment. Hence, the Reynolds averaged Navier-Stokes equations implemented in a diagonal form are applied to the cavity flow class of problems. Comparisons between numerical and experimental results are made in two-dimensions for free shear layers and both rectangular and quieted cavities, and in three-dimensions for Stratospheric Observatory For Infrared Astronomy (SOFIA) geometries. The acoustic behavior of the rectangular and three-dimensional cavity flows compare well with experiment in terms of frequency, magnitude, and quieting trends. However, there is a more rapid decrease in computed acoustic energy with frequency than observed experimentally owing to numerical dissipation. In addition, optical phase distortion due to the time-varying density field is modelled using geometrical constructs. The computed optical distortion trends compare with the experimentally inferred result, but underpredicts the fluctuating phase difference magnitude.

Atwood, C. A.↗

Analysis of a novel non-contacting waveguide backshort

A new non-contacting waveguide backshort has been developed for millimeter and submillimeter wave frequencies. The design consists of a metal bar with rectangular or circular holes cut into it, which is covered with a dielectric (mylar) layer to form a snug fit with the walls of a waveguide. Hole geometries are adjusted to obtain a periodic variation of the guide impedance on the correct length scale, in order to produce efficient reflection of RF power. It is a mechanically rugged design which can be easily fabricated for frequencies from 1 to 1000 GHz and is thus a sound alternative to the miniaturization of conventional non-contacting shorts. To aid in high-frequency design, a rigorous full-wave analysis has been completed, which will allow variations of the size, number and spacing of the holes to be easily analyzed. This paper will review the backshort design and the method developed for theoretical characterization, followed by a comparison of the experimental and numerical results. Low frequency models operating from 4-6 GHz are shown to demonstrate return loss of greater than -0.2 dB over a 33 percent bandwidth. The theory is in good agreement with measured data.

Weller, T. M.↗

Optimization of Neutral Atom Imagers

The interactions between plasma structures and neutral atom populations in interplanetary space can be effectively studied with energetic neutral atom imagers. For neutral atoms with energies less than 1 keV, the most efficient detection method that preserves direction and energy information is conversion to negative ions on surfaces. We have examined a variety of surface materials and conversion geometries in order to identify the factors that determine conversion efficiency. For chemically and physically stable surfaces smoothness is of primary importance while properties such as work function have no obvious correlation to conversion efficiency. For the noble metals, tungsten, silicon, and graphite with comparable smoothness, conversion efficiency varies by a factor of two to three. We have also examined the way in which surface conversion efficiency varies with the angle of incidence of the neutral atom and have found that the highest efficiencies are obtained at angles of incidence greater then 80deg. The conversion efficiency of silicon, tungsten and graphite were examined most closely and the energy dependent variation of conversion efficiency measured over a range of incident angles. We have also developed methods for micromachining silicon in order to reduce the volume to surface area over that of a single flat surface and have been able to reduce volume to surface area ratios by up to a factor of 60. With smooth micro-machined surfaces of the optimum geometry, conversion efficiencies can be increased by an order of magnitude over instruments like LENA on the IMAGE spacecraft without increase the instruments mass or volume.

Shappirio, M.↗

An approximate Riemann solver for magnetohydrodynamics (that works in more than one dimension)

An approximate Riemann solver is developed for the governing equations of ideal magnetohydrodynamics (MHD). The Riemann solver has an eight-wave structure, where seven of the waves are those used in previous work on upwind schemes for MHD, and the eighth wave is related to the divergence of the magnetic field. The structure of the eighth wave is not immediately obvious from the governing equations as they are usually written, but arises from a modification of the equations that is presented in this paper. The addition of the eighth wave allows multidimensional MHD problems to be solved without the use of staggered grids or a projection scheme, one or the other of which was necessary in previous work on upwind schemes for MHD. A test problem made up of a shock tube with rotated initial conditions is solved to show that the two-dimensional code yields answers consistent with the one-dimensional methods developed previously.

Powell, Kenneth G.↗

Nondestructive Evaluation for the Space Shuttle's Wing Leading Edge

The loss of the Space Shuttle Columbia highlighted concerns about the integrity of the Shuttle's thermal protection system, which includes Reinforced Carbon-Carbon (RCC) on the leading edge. This led NASA to investigate nondestructive evaluation (NDE) methods for certifying the integrity of the Shuttle's wing leading edge. That investigation was performed simultaneously with a large study conducted to understand the impact damage caused by errant debris. Among the many advanced NDE methods investigated for applicability to the RCC material, advanced digital radiography, high resolution computed tomography, thermography, ultrasound, acoustic emission and eddy current systems have demonstrated the maturity and success for application to the Shuttle RCC panels. For the purposes of evaluating the RCC panels while they are installed on the orbiters, thermographic detection incorporating principal component analysis (PCA) and eddy current array scanning systems demonstrated the ability to measure the RCC panels from one side only and to detect several flaw types of concern. These systems were field tested at Kennedy Space Center (KSC) and at several locations where impact testing was being conducted. Another advanced method that NASA has been investigating is an automated acoustic based detection system. Such a system would be based in part on methods developed over the years for acoustic emission testing. Impact sensing has been demonstrated through numerous impact tests on both reinforced carbon-carbon (RCC) leading edge materials as well as Shuttle tile materials on representative aluminum wing structures. A variety of impact materials and conditions have been evaluated including foam, ice, and ablator materials at ascent velocities as well as simulated hypervelocity micrometeoroid and orbital debris impacts. These tests have successfully demonstrated the capability to detect and localize impact events on Shuttle's wing structures. A first generation impact sensing system has been designed for the next Shuttle flight and is undergoing final evaluation for deployment on the Shuttle's first return to flight. This system will employ wireless accelerometer sensors that were qualified for other applications on previous Shuttle flights. These sensors will be deployed on the wing's leading edge to detect impacts on the RCC leading edge panels. The application of these methods will help to insure the continued integrity of the Shuttle wing's leading edge system as the Shuttle flights resume and until their retirement.

Eric I Madaras↗

CloudSat-Constrained Cloud Ice Water Path and Cloud Top Height Retrievals from MHS 157 and 183.3 GHz Radiances

Ice water path (IWP) and cloud top height (ht) are two of the key variables in determining cloud radiative and thermodynamical properties in climate models. Large uncertainty remains among IWP measurements from satellite sensors, in large part due to the assumptions made for cloud microphysics in these retrievals. In this study, we develop a fast algorithm to retrieve IWP from the 157, 183.3+/-3 and 190.3 GHz radiances of the Microwave Humidity Sounder (MHS) such that the MHS cloud ice retrieval is consistent with CloudSat IWP measurements. This retrieval is obtained by constraining the empirical forward models between collocated and coincident measurements of CloudSat IWP and MHS cloud-induced radiance depression (Tcir) at these channels. The empirical forward model is represented by a lookup table (LUT) of Tcir-IWP relationships as a function of ht and the frequency channel.With ht simultaneously retrieved, the IWP is found to be more accurate. The useful range of the MHS IWP retrieval is between 0.5 and 10 kg/sq m, and agrees well with CloudSat in terms of the normalized probability density function (PDF). Compared to the empirical model, current operational radiative transfer models (RTMs) still have significant uncertainties in characterizing the observed Tcir-IWP relationships. Therefore, the empirical LUT method developed here remains an effective approach to retrieving ice cloud properties from the MHS-like microwave channels.

cloud top height↗

Chemical Reactivity of In-Situ Lunar Dust for Biotoxicity Assessment

How does the chemical reactivity of in-situ lunar dust compare to Apollo samples currently stored in curation facilities here on Earth? Essential investigations of this question will help us to further mitigate exploration risks for future human explorers on the Moon and will also provide critical information for astrobiologists and space biologists using the Moon for scientific inquiry. Apollo 14 dust biotoxicity studies, carried out by the NASA Lunar Airborne Dust Toxicity Assessment Group (LADTAG), included numerous cellular and animal experiments. Intratracheal instillation and inhalation studies in rats both showed Apollo 14 dust to be intermediate in toxicity compared to low-tox titanium dusts and high-tox quartz dusts of similar particle sizes. The collective results were used in models to establish a safe exposure limit for astronauts. Although LADTAG took extensive steps to preserve what chemical reactivity may still have existed in the samples, it is simply unknown if they possessed true in-situ chemical reactivity or if that reactivity has decayed. Initial gas loss on collection and other alterations, and even intermittent exposure to Earth-normal conditions during subsequent decades of handling, obscure a forensic reconstruction of the initial state. Because a mineral dust’s chemical reactivity influences its biotoxicity, researchers have developed methods to “activate” lunar dust and simulants. Past studies that modeled impact processes and radiation in the lunar environment suggest that in-situ lunar dust is likely to be more chemically reactive than Earth-exposed samples. Because of these results, in-situ measurements are warranted. Since the lunar surface is heterogeneous, dust biotoxicity is expected to vary from site to site due to particle size, mineralogy, physical characteristics, degree of space weathering, and chemical reactivity. This circumstance dictates dust assessments at a suite of lunar sites enabled by CLPS opportunities. Dose, location, and duration of particle exposure will also affect biological responses. In-situ chemical reactivity measurements can inform cross-cutting collaborative research campaigns such as astrobiology studies examining regolith interactions with organisms and its ability to preserve chemical and structural biomarkers, as well as space biology investigations that examine regolith-microbe interactions relating to life support systems, plant growth, biomining, and development of regolith biocomposites.

Jon C Rask↗

Aeroacoustic sensitivity analysis and optimal aeroacoustic design of turbomachinery blades

During the first year of the project, we have developed a theoretical analysis - and wrote a computer code based on this analysis - to compute the sensitivity of unsteady aerodynamic loads acting on airfoils in cascades due to small changes in airfoil geometry. The steady and unsteady flow though a cascade of airfoils is computed using the full potential equation. Once the nominal solutions have been computed, one computes the sensitivity. The analysis takes advantage of the fact that LU decomposition is used to compute the nominal steady and unsteady flow fields. If the LU factors are saved, then the computer time required to compute the sensitivity of both the steady and unsteady flows to changes in airfoil geometry is quite small. The results to date are quite encouraging, and may be summarized as follows: (1) The sensitivity procedure has been validated by comparing the results obtained by 'finite difference' techniques, that is, computing the flow using the nominal flow solver for two slightly different airfoils and differencing the results. The 'analytic' solution computed using the method developed under this grant and the finite difference results are found to be in almost perfect agreement. (2) The present sensitivity analysis is computationally much more efficient than finite difference techniques. We found that using a 129 by 33 node computational grid, the present sensitivity analysis can compute the steady flow sensitivity about ten times more efficiently that the finite difference approach. For the unsteady flow problem, the present sensitivity analysis is about two and one-half times as fast as the finite difference approach. We expect that the relative efficiencies will be even larger for the finer grids which will be used to compute high frequency aeroacoustic solutions. Computational results show that the sensitivity analysis is valid for small to moderate sized design perturbations. (3) We found that the sensitivity analysis provided important insight into how airfoils should be modified to improve their aeroelastic stability. Using this insight, we redesigned an aeroelastically unstable cascade to produce a stable cascade.

Hall, Kenneth C.↗

Development of a Liner Design Methodology and Relevant Results of Acoustic Suppression in the Farfield for Mixer-Ejector Nozzles

We have developed a process to predict noise field interior to the ejector and in the farfield for any liner design for a mixer-ejector of arbitrary scale factor. However, a number of assumptions, not verified for the current application, utilized in this process, introduce uncertainties in the final result, especially, on a quantitative basis. The normal impedance model for bulk with perforated facesheet is based on homogeneous foam materials of low resistivity. The impact of flow conditions for HSCT application as well as the impact of perforated facesheet on predicted impedance is not properly accounted. Based on the measured normal impedance for deeper bulk samples (i.e., 2.0 in.) the predicted reactance is much higher compared to the data at frequencies above 2 kHz for T-foam and 200 ppi SiC. The resistance is under predicted at lower frequencies (below 4 kHz) for these samples. Thus, the use of such predicted data in acoustic suppression is likely to introduce inaccuracies. It should be noted that the impedance prediction methods developed recently under liner technology program are not utilized in the studies described in this report due to the program closeout. Acoustic suppression prediction is based on the uniform flow and temperature conditions in a two-sided treated constant area rectangular duct. In addition, assumptions of equal energy per mode noise field and interaction of all frequencies with the treated surface for the entire ejector length may not be accurate. While, the use of acoustic transfer factor minimizes the inaccuracies associated with the prediction for a known test case, the assumption of the same factor for other liner designs and with different linear scale factor ejectors seems to be very optimistic. As illustrated in appendix D that the predicted noise suppression for LSM-1 is lower compared to the measured data is an indication of the above argument. However, the process seems to be more reliable when used for the same scale models for different liner designs as demonstrated for Gen. 1 mixer-ejectors.

Salikuddin, M.↗

Designing Agent Utilities for Coordinated, Scalable and Robust Multi-Agent Systems

Coordinating the behavior of a large number of agents to achieve a system level goal poses unique design challenges. In particular, problems of scaling (number of agents in the thousands to tens of thousands), observability (agents have limited sensing capabilities), and robustness (the agents are unreliable) make it impossible to simply apply methods developed for small multi-agent systems composed of reliable agents. To address these problems, we present an approach based on deriving agent goals that are aligned with the overall system goal, and can be computed using information readily available to the agents. Then, each agent uses a simple reinforcement learning algorithm to pursue its own goals. Because of the way in which those goals are derived, there is no need to use difficult to scale external mechanisms to force collaboration or coordination among the agents, or to ensure that agents actively attempt to appropriate the tasks of agents that suffered failures. To present these results in a concrete setting, we focus on the problem of finding the sub-set of a set of imperfect devices that results in the best aggregate device. This is a large distributed agent coordination problem where each agent (e.g., device) needs to determine whether to be part of the aggregate device. Our results show that the approach proposed in this work provides improvements of over an order of magnitude over both traditional search methods and traditional multi-agent methods. Furthermore, the results show that even in extreme cases of agent failures (i.e., half the agents failed midway through the simulation) the system's performance degrades gracefully and still outperforms a failure-free and centralized search algorithm. The results also show that the gains increase as the size of the system (e.g., number of agents) increases. This latter result is particularly encouraging and suggests that this method is ideally suited for domains where the number of agents is currently in the thousands and will reach tens or hundreds of thousands in the near future.

Tumer, Kagan↗

Mechanically fastened composite laminates subjected to combined bearing-bypass and shear loading

Bolts and rivets provide a means of load transfer in the construction of aircraft. However, they give rise to stress concentrations and are often the source and location of static and fatigue failures. Furthermore, fastener holes are prone to cracks during take-off and landing. These cracks present the most common origin of structural failures in aircraft. Therefore, accurate determination of the contact stresses associated with such loaded holes in mechanically fastened joints is essential to reliable strength evaluation and failure prediction. As the laminate is subjected to loading, the contact region, whose extent is not known, develops between the fastener and the hole boundary through this contact region, which consists of slip and no-slip zones due to friction. The presence of the unknown contact stress distribution over the contact region between the pin and the composite laminate, material anisotropy, friction between the pin and the laminate, pin-hole clearance, combined bearing-bypass and shear loading, and finite geometry of the laminate result in a complex non-linear problem. In the case of bearing-bypass loading in compression, this non-linear problem is further complicated by the presence of dual contact regions. Previous research concerning the analysis of mechanical joints subjected to combined bearing-bypass and shear loading is non-existent. In the case of bearing-bypass loading only, except for the study conducted by Naik and Crews (1991), others employed the concept of superposition which is not valid for this non-linear problem. Naik and Crews applied a linear finite element analysis with conditions along the pin-hole contact region specified as displacement constraint equations. The major shortcoming of this method is that the variation of the contract region as a function of the applied load should be known a priori. Also, their analysis is limited to symmetric geometry and material systems, and frictionless boundary conditions. Since the contact stress distribution and the contact region are not known a priori, they did not directly impose the boundary conditions appropriate for modelling the contact and on-contact regions between the fastener and the hole. Furthermore, finite element analysis is not suitable for iterative design calculations for optimizing laminate construction in the presence of fasteners under complex loading conditions. In this study, the solution method developed by Madenci and Ileri (1992a,b) has been extended to determine the contact stresses in mechanical joints under combined bearing-bypass and shear loading, and bearing-bypass loading in compression resulting in dual contact regions.

Madenci, Erdogan↗