Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deformation estimation”

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.

411 records · Page 23

Detection of atmospheric pressure loading using very long baseline interferometry measurements

Loading of the Earth by the temporal redistribution of global atmospheric mass is likely to displace the positions of geodetic monuments by tens of millimeters both vertically and horizontally. Estimates of these displacements are determined by convolving National Meteorological Center (NMC) global values of atmospheric surface pressure with Farrell's elastic Green's functions. An analysis of the distances between radio telescopes determined by very long baseline interferometry (VLBI) between 1984 and 1992 reveals that in many of the cases studied there is a significant contribution to baseline length change due to atmospheric pressure loading. Our analysis covers intersite distances of between 1000 and 10,000 km and is restricted to those baselines measured more than 100 times. Accounting for the load effects (after first removing a best fit slope) reduces the weighted root-mean-square (WRMS) scatter of the baseline length residuals on 11 of the 22 baselines investigated. The slight degradation observed in the WRMS scatter on the remaining baselines is largely consistent with the expected statistical fluctuations when a small correction is applied to a data set having a much larger random noise. The results from all baselines are consistent with approximately 60% of the computed pressure contribution being present in the VLBI length determinations. Site dependent coefficients determined by fitting local pressure to the theoretical radial displacement are found to reproduce the deformation caused by the regional pressure to within 25% for most inland sites. The coefficients are less reliable at near coastal and island stations.

Vandam, T. M.↗

Beyond Basic Drag in Interplanetary CME Modeling: Effects of Solar Wind Pileup and High-Speed Streams

Coronal mass ejections (CMEs) cause severe space weather effects throughout our solar system. As a fast CME propagates through interplanetary space, it accumulates solar wind materials at its front. This pileup of materials, or CME-driven sheath, can be important in determining the geoeffectiveness of a CME. We take an existing arrival time model that includes expansion and deformation of the CME flux rope (ANTEATR; Kay & Gopalswamy, 2018, https://doi.org/10.1029/2018JA025780; Kay & Nieves-Chinchilla, 2021a, https://doi.org/10.1029/2020JA028911) and add a pileup procedure (PUP) as a physics-based approach to modeling the CME-driven sheath. ANTEATR-PUP solves the Rankine-Hugoniot equations for an oblique shock to determine the shock speed and sheath density, magnetic field, and temperature. The extra sheath mass affects the background drag calculation. Additionally, ANTEATR can now use any 1D profile for the background solar wind as opposed to the simple empirical models it previously relied upon. We present initial results from ANTEATR-PUP and compare with previous ANTEATR findings. Using results from an MHD simulation, we explore the effects of interactions with a static high-speed stream (HSS) on the CME's and sheath's interplanetary evolution. The drag forces essentially disappear while a CME remains within the HSS, but reappear stronger once the CME exits. The HSS-CME interaction produces the largest changes in the CME and sheath properties at 1 au when it occurs either close to the Sun near the inner simulation boundary at 0.1 au or right before the CME reaches 1 au. We estimate that these changes could significantly affect the geoeffectiveness.

C Kay↗

CREME96 and Related Error Rate Prediction Methods

Predicting the rate of occurrence of single event effects (SEEs) in space requires knowledge of the radiation environment and the response of electronic devices to that environment. Several analytical models have been developed over the past 36 years to predict SEE rates. The first error rate calculations were performed by Binder, Smith and Holman. Bradford and Pickel and Blandford, in their CRIER (Cosmic-Ray-Induced-Error-Rate) analysis code introduced the basic Rectangular ParallelePiped (RPP) method for error rate calculations. For the radiation environment at the part, both made use of the Cosmic Ray LET (Linear Energy Transfer) spectra calculated by Heinrich for various absorber Depths. A more detailed model for the space radiation environment within spacecraft was developed by Adams and co-workers. This model, together with a reformulation of the RPP method published by Pickel and Blandford, was used to create the CR ME (Cosmic Ray Effects on Micro-Electronics) code. About the same time Shapiro wrote the CRUP (Cosmic Ray Upset Program) based on the RPP method published by Bradford. It was the first code to specifically take into account charge collection from outside the depletion region due to deformation of the electric field caused by the incident cosmic ray. Other early rate prediction methods and codes include the Single Event Figure of Merit, NOVICE, the Space Radiation code and the effective flux method of Binder which is the basis of the SEFA (Scott Effective Flux Approximation) model. By the early 1990s it was becoming clear that CREME and the other early models needed Revision. This revision, CREME96, was completed and released as a WWW-based tool, one of the first of its kind. The revisions in CREME96 included improved environmental models and improved models for calculating single event effects. The need for a revision of CREME also stimulated the development of the CHIME (CRRES/SPACERAD Heavy Ion Model of the Environment) and MACREE (Modeling and Analysis of Cosmic Ray Effects in Electronics). The Single Event Figure of Merit method was also revised to use the solar minimum galactic cosmic ray spectrum and extended to circular orbits down to 200 km at any inclination. More recently a series of commercial codes was developed by TRAD (Test & Radiations) which includes the OMERE code which calculates single event effects. There are other error rate prediction methods which use Monte Carlo techniques. In this chapter the analytic methods for estimating the environment within spacecraft will be discussed.

Adams, James H., Jr.↗

Hardness of CaF2 and BaF2 solid lubricants at 25 to 670 deg C

Plastic deformation is a prominent factor in determining the lubricating value of solid lubricants. Little information is available and its direct measurement is difficult so hardness, which is an indirect measure of this property was determined for fluoride solid lubricant compositions. The Vickers hardness of BaF2 and CaF2 single crystals was measured up to 670 C in a vacuum. The orientation of the BaF2 was near the (013) plane and the CaF2 was about 16 degrees from the degrees from the (1'11) plane. The BaF2 has a hardness of 83 kg/sq mm at the 25 C and 9 at the 600 C. The CaF2 is 170 at 25 C and 13 at 670 C. The decrease in hardness in the temperature range of 25 to 100 C is very rapid and amounts to 40% for both materials. Melts of BaF2 and CaF2 were made in a platinum crucible in ambient air with compositions of 50 to 100 wt% BaF2. The Vickers hardness of these polycrystalline binary compositions at 25 C increased with increasing CaF2 reaching a maximum of 150 kn/sq mm near the eutectic. The polycrystalline CaF2 was 14% softer than that of the single crystal surface and BsF2 was 30% harder than the single crystal surface. It is estimated that the brittle to ductile transition temperature for CaF2 and BaF2 is less than 100 C for the conditions present in the hardness tester.

Deadmore, Daniel L.↗

A Refined Zigzag Beam Theory for Composite and Sandwich Beams

A new refined theory for laminated composite and sandwich beams that contains the kinematics of the Timoshenko Beam Theory as a proper baseline subset is presented. This variationally consistent theory is derived from the virtual work principle and employs a novel piecewise linear zigzag function that provides a more realistic representation of the deformation states of transverse-shear flexible beams than other similar theories. This new zigzag function is unique in that it vanishes at the top and bottom bounding surfaces of a beam. The formulation does not enforce continuity of the transverse shear stress across the beam s cross-section, yet is robust. Two major shortcomings that are inherent in the previous zigzag theories, shear-force inconsistency and difficulties in simulating clamped boundary conditions, and that have greatly limited the utility of these previous theories are discussed in detail. An approach that has successfully resolved these shortcomings is presented herein. Exact solutions for simply supported and cantilevered beams subjected to static loads are derived and the improved modelling capability of the new zigzag beam theory is demonstrated. In particular, extensive results for thick beams with highly heterogeneous material lay-ups are discussed and compared with corresponding results obtained from elasticity solutions, two other zigzag theories, and high-fidelity finite element analyses. Comparisons with the baseline Timoshenko Beam Theory are also presented. The comparisons clearly show the improved accuracy of the new, refined zigzag theory presented herein over similar existing theories. This new theory can be readily extended to plate and shell structures, and should be useful for obtaining relatively low-cost, accurate estimates of structural response needed to design an important class of high-performance aerospace structures.

Tessler, Alexander↗

Experimental Study of Hypersonic Inflatable Aerodynamic Decelerator (HIAD) Aeroshell with Axisymmetric Surface Deflection Patterns

A wind tunnel test program was conducted to obtain aeroheating environment data on Hypersonic Inflatable Aerodynamic Decelerator aeroshells with flexible thermal protection systems. Data were obtained on a set of rigid wind tunnel models with surface deflection patterns of various heights that simulated a range of potential in-flight aeroshell deformations. Wind tunnel testing was conducted at Mach 6 at unit Reynolds numbers from 2.1 × 10(exp 6)/ft to 8.3 × 10(exp 6)/ft and angles of attack from 0 deg to 18 deg. Boundary-layer transition onset and global surface heating distribution measurements were performed using phosphor thermography and flow field images were obtained through schlieren photography. Surface deflections were found to both promote early transition of the boundary layer and to augment heating levels for both laminar and turbulent flows. A complimentary computational flow field study was also performed to provide heating predictions for comparison with the measurements as well as boundary layer flow field properties for use in correlating the data. Correlations of the wind tunnel data were developed to predict deflection effects on boundary layer transition and surface heating and were applied to both the wind tunnel test conditions and to the trajectory of NASA's successful IRVE-3 flight test. In general, the correlations produced at least qualitative agreement with the wind tunnel data, although the heating levels were underpredicted for some of the larger surface deflections. For the flight conditions, the correlations suggested that peak heating levels on the leeward side conical flank of the IRVE-3 vehicle may have exceeded those at nose for times late in the trajectory after the peak heating time point. However, the flight estimates were based on a conservative assumption of surface deflection magnitude (i.e., larger) than likely was produced in flight.

Hollis, Brian R.↗

The deep structure of Venusian plateau highlands

Magellan gravity data confirm that several of the large, tectonically deformed, plateau-like highlands on venus are shallowly compensated, most likely by crustal thickness variations. Apparent depths of isostatic compensation, computed in the spatial domain, range from 30 to 50 km for Alpha, Tellus, Ovda, and Thetis Regiones. Using a two-layer model for isostatic compensation, Alpha, Tellus, and Ovda are best represented as nearly completely compensated in crust that is regionally 20-40 km thick around these highlands, with little contribution from deeper mantle sources. In contrast to these three areas, a stronger regional gravity high associated with Thetis requires a significant upper mantle component to compensation. This is evident in the spectral admittance as a pronounced deep, long-wavelength anomaly. In the two-layer isostatic model, a broad, deeply compensated upland underlies a shallowly compensated central block of Thetis. If this deep component is interpreted as a thermal anomaly, the loci of maximum upwelling agree well with sites of recent extension. The plateau highlands are thus physiographically and isostatically equivalent to terrestrial continents, though probably not compositionally. They also share the record of a long tectonic history. The large regional gravity anomaly of Thetis indicates that active mantle proceses continue even beneath some areas (tessera) thought to be a relic of a former geological regime. The excellent agreement of modeled crustal thicknesses around Alpha, Tellus, and Ovda Regiones suggests that 20-40 km is a representative global value for the plains. Such a crust is thicker than previously estimated and about twice as thick as the expected thickness of crust produced at venusian spreading centers

Grimm, Robert E.↗

Palinspastic reconstruction of southeastern California and southwestern Arizona for the middle Miocene

A paleogeographic reconstruction of southeastern California and southwestern Arizona at 10 Ma was made based on available geologic and geophysical data. Clockwise rotation of 39 deg was reconstructed in the eastern Transverse Ranges, consistent with paleomagnetic data from late Miocene volcanic rocks, and with slip estimates for left-lateral faults within the eastern Transverse Ranges and NW-trending right lateral faults in the Mojave Desert. This domain of rotated rocks is bounded by the Pinto Mountain fault on the north. In the absence of evidence for rotation of the San Bernardino Mountains or for significant right slip faults within the San Bernardino Mountains, the model requires that the late Miocene Pinto Mountain fault become a thrust fault gaining displacement to the west. The Squaw Peak thrust system of Meisling and Weldon may be a western continuation of this fault system. The Sheep Hole fault bounds the rotating domain on the east. East of this fault an array of NW-trending right slip faults and south-trending extensional transfer zones has produced a basin and range physiography while accumulating up to 14 km of right slip. This maximum is significantly less than the 37.5 km of right slip required in this region by a recent reconstruction of the central Mojave Desert. Geologic relations along the southern boundary of the rotating domain are poorly known, but this boundary is interpreted to involve a series of curved strike slip faults and non-coaxial extension, bounded on the southeast by the Mammoth Wash and related faults in the eastern Chocolate Mountains. Available constraints on timing suggest that Quaternary movement on the Pinto Mountain and nearby faults is unrelated to the rotation of the eastern Transverse Ranges, and was preceded by a hiatus during part of Pliocene time which followed the deformation producing the rotation. The reconstructed Clemens Well fault in the Orocopia Mountains, proposed as a major early Miocene strand of the San Andreas fault, projects eastward towards Arizona, where early Miocene rocks and structures are continuous across its trace. The model predicts a 14 deg clockwise rotation and 55 km extension along the present trace of the San Andreas fault during late Miocene and early Pliocene time. Palinspastic reconstructions of the San Andreas system based on this proposed reconstruction may be significantly modified from current models.

Richard, Stephen M.↗

Development, Characterization, and Validation of Elevated-Temperature Constitutive Models: Deformation and Damage

This report provides a brief review of experimentally observed hereditary and nonhereditary material behavior along with background information on standard as well as advanced internal state variable constitutive modeling at elevated temperature. A description of exploratory, characterization, and validation testing is presented along with a detailed outline of what constitutes “sufficient” data content (i.e., quality and quantity) for developing or enhancing, characterizing, and validating a particular sophisticated nonlinear time- and history-dependent (hereditary) class of constitutive models known as GVIPS (generalized viscoplasticity with potential structure). The tests described are necessary to reveal a material’s behavior in the reversible (or viscoelastic) and irreversible (or viscoplastic) regimes, for the identification of both deformation and damage model parameters. Results presented are primarily for metallic materials. In all cases, both uniaxial and multiaxial tests are described, and the linkage between the specific tests and the parameters within the model that can be characterized from the results of these tests are also discussed. Discussion is also provided relative to the role information management must play relative to material data collection, analysis, maintenance, and dissemination. The need for such an information system is particularly important as both analyst and designer move toward utilizations of sophisticated, nonlinear time- and history dependent (hereditary) constitutive models. Lastly, the concept of state space and its utility in understanding and establishing constitutive models is addressed throughout. The intent behind this document is to help both the modeler and experimentalist understand each other’s specific points of view and provide guidance for both model development and characterization. Emphasis has been placed on providing the mechanician with information regarding how tests are performed and what issues to be aware of when interpreting results. It is hoped that experimentalists will take away a new perspective on the types of information that modelers and analysts are looking for from them.

Constitutive Modeling↗

The Evolution of Friction Stir Welding Theory at Marshall Space Flight Center

From 1995 to the present the friction stir welding (FSW) process has been under study at Marshall Space Flight Center (MSFC). This is an account of the progressive emergence of a set of conceptual tools beginning with the discovery of the shear surface, wiping metal transfer, and the invention of a kinematic model and making possible a treatment of both metallurgical structure formation and process dynamics in friction stir welding from a unified point of view. It is generally observed that the bulk of the deformation of weld metal around the FSW pin takes place in a very narrow, almost discontinuous zone with high deformation rates characteristic of metal cutting. By 1999 it was realized that this zone could be treated as a shear surface like that in simple metal cutting models. At the shear surface the seam is drawn out and compressed and pressure and flow conditions determine whether or not a sound weld is produced. The discovery of the shear surface was followed by the synthesis of a simple 3- flow kinematic model of the FSW process. Relative to the tool the flow components are: (1) an approaching translational flow at weld speed V, (2) a rotating cylindrical plug flow with the angular velocity of the tool , and (3) a relatively slow ring vortex flow (like a smoke ring) encircling the tool and driven by shoulder scrolls and pin threads. The rotating plug flow picks up an element of weld metal, rotates it around with the tool, and deposits it behind the tool ( wiping metal transfer ); it forms plan section loops in tracers cut through by the tool. Radially inward flow from the ring vortex component retains metal longer in the rotating plug and outward flow expels metal earlier; this interaction forms the looping weld seam trace and the tongue and groove bimetallic weld contour. The radial components of the translational and ring vortex flows introduce parent metal intrusions into the small grained nugget material close to the tool shoulder; if this feature is pronounced, nugget collapse may result. Certain weld features, in particular internal banding seen in transverse section as onion rings and associated surface ridges called tool marks , have long implied an oscillation flow component, but have only recently been attributed in the literature to tool eccentricity. Rotating plug shape, typically a hollow cylinder flared at the end where it sticks to the shoulder, varies as pressure distribution on the tool determines where sticking occurs. Simplified power input estimates balanced against heat loss estimates give reasonable temperature estimates, explain why the power requirement changes hardly at all over a wide range of RPM s, and yield isotherms that seem to fall along boundaries of parameter windows of operation.

Nunes, Arthur C.↗

Venus tectonic styles and crustal differentiation

Two of the most important constraints are known from Pioneer Venus data: the lack of a system of spreading rises, indicating distributed deformation rather than plate tectonics; and the high gravity/topography ratio, indicating the absence of an asthenosphere. In addition, the high depth/diameter ratios of craters on Venus indicate that Venus probably has no more crust than Earth. The problems of the character of tectonics and crustal formation and recycling are closely coupled. Venus appears to lack a recycling mechanism as effective as subduction, but may also have a low rate of crustal differentiation because of a mantle convection pattern that is more distributed, less concentrated, than Earth's. Distributed convection, coupled with the nonlinear dependence of volcanism on heat flow, would lead to much less magmatism, despite only moderately less heat flow, compared to Earth. The plausible reason for this difference in convective style is the absence of water in the upper mantle of Venus. We have applied finite element modeling to problems of the interaction of mantle convection and crust on Venus. The main emphasis has been on the tectonic evolution of Ishtar Terra, as the consequence of convergent mantle flow. The early stage evolution is primarily mechanical, with crust being piled up on the down-stream side. Then the downflow migrates away from the center. In the later stages, after more than 100 m.y., thermal effects develop due to the insulating influence of the thickened crust. An important feature of this modeling is the entrainment of some crustal material in downflows. An important general theme in both convergent and divergent flows is that of mixing vs. stratification. Models of multicomponent solid-state flow obtain that lower-density crustal material can be entrained and recycled, provided that the ration of low-density to high-density material is small enough (as in subducted slabs on Earth). The same considerations should apply in upflows; a small percent of partial melt may be carried along with its matrix and never escape to the surface. Models that assume melt automatically rising to the crust and no entrainment or other mechanism of recycling lower-density material obtain oscillatory behavior, because it takes a long time for heat to build up enough to overcome a Mg-rich low-density residuum. However, these models develop much thicker crust than consistent with estimates from crater depth/diameter ratios.

Kaula, W. M.↗

Overcoming Dynamic Disturbances in Imaging Systems

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, cost, and risk of space-borne optical systems in the presence of dynamic disturbances. More robust systems almost certainly will be a result as well. Many future space-based and ground-based optical systems will employ optical control systems to enhance imaging performance. The goal of the optical control subsystem is to determine the wavefront aberrations and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable mirror(s), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic control. These control subsystems may become elaborate indeed. But robust system performance will require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or atmospheric sensors for dynamic imaging of individual, severe storms. The technology developed here could have a substantial impact on the development of new systems in the Space Science Enterprise; such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optomechanical structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging of planets or the sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small targets as well and could benefit from the technology developed here. We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should work for all but the lowest light levels and highest frequency dynamic environments. This technique may serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk operation. For a given object, a short exposure image is "frozen" on the focal plane in the presence of the environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to-frame linear shift invariance. Therefore, although the Point Spread Function is varying from frame to frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limit of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. The magnitude of the complex object spectrum is straightforward to determine by well-established approaches. The phase of the complex object spectrum is decomposed into two parts. One is a single-valued function determined by the divergence of the optical phase gradient. The other is a multi-valued function determined by the circulation of the optical phase gradient-"hidden phase." Finite difference equations are developed for the phase. The novelty of this approach is captured in the inclusion of this "hidden phase." This technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure frames while simultaneously estimating the phase as a function of time from a set of exposures.

Young, Eric W.↗

Overcoming Dynamic Disturbances in Imaging Systems

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, cost, and risk of space-borne optical systems in the presence of dynamic disturbances. More robust systems almost certainly will be a result as well. Many future space-based and ground-based optical systems will employ optical control systems to enhance imaging performance. The goal of the optical control subsystem is to determine the wavefront aberrations and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable mirror(s), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic control. These control subsystems may become elaborate indeed. But robust system performance will require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or atmospheric sensors for dynamic imaging of individual, severe storms. The technology developed here could have a substantial impact on the development of new systems in the Space Science Enterprise; such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optormechanical structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging of planets or the sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small targets as well and could benefit from the technology developed here. We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should work for all but the lowest light levels and highest frequency dynamic environments. This technique may serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk operation. For a given object, a short exposure image is "frozen" on the focal plane in the presence of the environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to-frame linear shift invariance. Therefore, although the Point Spread Function is varying from frame to frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limit of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. The magnitude of the complex object spectrum is straightforward to determine by well-established approaches. The phase of the complex object spectrum is decomposed into two parts. One is a single-valued function determined by the divergence of the optical phase gradient. The other is a multi-valued function determined by, the circulation of the optical phase gradient-"hidden phase." Finite difference equations are developed for the phase. The novelty of this approach is captured in the inclusion of this "hidden phase." This technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure frames while simultaneously estimating the phase as a function of time from a set of exposures.

Young, Eric W.↗

Analysis and Design of Variable Stiffness Composite Cylinders

An investigation of the possible performance improvements of thin circular cylindrical shells through the use of the variable stiffness concept is presented. The variable stiffness concept implies that the stiffness parameters change spatially throughout the structure. This situation is achieved mainly through the use of curvilinear fibers within a fiber-reinforced composite laminate, though the possibility of thickness variations and discrete stiffening elements is also allowed. These three mechanisms are incorporated into the constitutive laws for thin shells through the use of Classical Lamination Theory. The existence of stiffness variation within the structure warrants a formulation of the static equilibrium equations from the most basic principles. The governing equations include sufficient detail to correctly model several types of nonlinearity, including the formation of a nonlinear shell boundary layer as well as the Brazier effect due to nonlinear bending of long cylinders. Stress analysis and initial buckling estimates are formulated for a general variable stiffness cylinder. Results and comparisons for several simplifications of these highly complex governing equations are presented so that the ensuing numerical solutions are considered reliable and efficient enough for in-depth optimization studies. Four distinct cases of loading and stiffness variation are chosen to investigate possible areas of improvement that the variable stiffness concept may offer over traditional constant stiffness and/or stiffened structures. The initial investigation deals with the simplest solution for cylindrical shells in which all quantities are constant around the circumference of the cylinder. This axisymmetric case includes a stiffness variation exclusively in the axial direction, and the only pertinent loading scenarios include constant loads of axial compression, pressure, and torsion. The results for these cases indicate that little improvement over traditional laminates exists through the use of curvilinear fibers, mainly due to the presence of a weak link area within the stiffness variation that limits the ultimate load that the structure can withstand. Rigorous optimization studies reveal that even though slight increases in the critical loads can be produced for designs with an arbitrary variation of the fiber orientation angle, the improvements are not significant when compared to traditional design techniques that utilize ring stiffeners and frames. The second problem that is studied involves arbitrary loading of a cylinder with a stiffness variation that changes only in the circumferential direction. The end effects of the cylinder are ignored, so that the problem takes the form of an analysis of a cross-section for a short cylinder segment. Various load cases including axial compression, pressure, torsion, bending, and transverse shear forces are investigated. It is found that the most significant improvements in load-carrying capability exist for cases which involve loads that also vary around the circumference of the shell, namely bending and shear forces. The stiffness variation of the optimal designs contribute to the increased performance in two ways: lowering the stresses in the critical areas through redistribution of the stresses; and providing a relatively stiff region that alters the buckling behavior of the structure. These results lead to an in-depth optimization study involving weight optimization of a fuselage structure subjected to typical design constraints. Comparisons of the curvilinear fiber format to traditional stiffened structures constructed of isotropic and composite materials are included. It is found that standard variable stiffness designs are quite comparable in terms of weight and load-carrying capability yet offer the added advantage of tailorability of distinct regions of the structure that experience drastically different loading conditions. The last two problems presented in this work involve the nonlinear phenomenon of long tubes under bending. Though this scenario is not as applicable to fuselage structures as the previous problems, the mechanisms that produce the nonlinear effect are ideally suited to be controlled by the variable stiffness concept. This is due to the fact that the dominating influence for long cylinders under bending is the ovalization of the cross-section, which is governed mainly by the stiffness parameters of the cylindrical shell. Possible improvement of the critical buckling moments for these structures is investigated using either a circumferential or axial stiffness variation. For the circumferential case involving infinite length cylinders, it is found that slight improvements can be observed by designing structures that resist the cross-sectional deformation yet do not detract from the buckling resistance at the critical location. The results also indicate that buckling behavior is extremely dependent on cylinder length. This effect is most easily seen in the solution of finite length cylinders under bending that contain an axial stiffness variation. For these structures, the only mechanism that exhibits improved response are those that effectively shorten the length of the cylinder, thus reducing the cross-sectional deformation due to the forced restraint at the ends. It was found that the use of curvilinear fibers was not able to achieve this effect in sufficient degree to resist the deformation, but that ring stiffeners produced the desired response admirably. Thus, it is shown that the variable stiffness concept is most effective at improving the bending response of long cylinders through the use of a circumferential stiffness variation.

Tatting, Brian F.↗

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗