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.

At least 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Development of Additively Manufactured Complex Tools for Autoclave Cure Composites

IACMI Project 4.9, Tooling for Composites with Washout Features Produced by Additive Manufacturing, assembled a team including the industry lead, Ability Composites, NREL and Colorado State University (CSU). Ability Composites had originally expressed interest in alternate methods of producing tooling for composite parts. In follow-up discussions, it became clear that one of the bigger tooling challenges revolved around small production volume composite parts that were tooled on washout material due to the complex geometry. To build an understanding of the potential, both from a technology and a cost perspective, for replacing conventional washout tooling with 3D printed thermoplastic tooling, a number of commercially available dissolvable FDM printing materials were evaluated, leading to tooling representative of commercial articles of interest to Ability Composites. Ultimately, Ability Composites was able to directly compare autoclave processed prepreg composite parts produced on conventional washout tooling to composite parts molded on 3D printed dissolvable tooling produced at CSU. Small, laboratory test specimens were developed to investigate the structural performance of the candidate materials under autoclave processing conditions, which were nominally 121 °C (250 °F) and 345 kPa (50 psi). In addition, several internal structural configurations (infills) were evaluated under autoclave conditions using model materials. The results of these tests indicated that two materials, Stratasys ST 130 and Infinite Materials Solutions Aquasys 180 (AQ 180), were the best candidates, given the specified autoclave processing conditions. ST-130 was slightly more robust than AQ-180; however, the AQ-180 was carried forward as it was dissolvable in water, not requiring the basic solution needed to dissolve ST-130. Based on the preliminary material and 3D printed structures evaluations, larger tools with a truncated square pyramid geometry were created to produce prepreg composite test articles for 3D printed dissolvable tool evaluation under standard autoclave fabrication conditions. Two tools were manufactured using ST-130 and one tool using traditional ceramic washout tooling media. The tools were evaluated for geometric fidelity and surface roughness changes before and after carbon fiber/epoxy prepreg composites were manufactured on the tooling. The autoclave processing did not impact the geometry significantly and was completed at 121 °C and 345 kPa, indicating satisfactory tool performance. The results from surface roughness testing of both the resulting composite and the associated tooling indicated that an adequate surface resulted without the need for a surface sealing step, as was required for the conventional washout tooling. Based on results of the truncated pyramid tests as a basis, ST-130, AQ-120 and AQ-180 materials were carried forward to the tool geometry of interest to Ability Composites. These hollow rectangular bent ducts, which were complex in nature and not extractable after cure, were used to understand the impacts of tool material and thickness. One ST-130 tool was produced as a partially solid part, with an enclosed 40% dense infill region to reduce weight and material use. This was the same approach evaluated in the truncated pyramid portion of the study. This tool was to be envelope vacuum bagged and directly compared to a monolithic tool of conventional washout material. The traditional monolithic ceramic tool was manufactured by Ability composites using CNC-based subtractive methods. An additional five dissolvable polymer tools, manufactured from ST-130, AQ-120, and AQ-180, using a hollow design were 3D printed and used to produce carbon fiber/epoxy composite evaluation articles. These hollow dissolvable tools were expected to be less influenced by the autoclave conditions as the wall was solid and vacuum bagging was inside and outside the tool. This alternative geometry was also evaluated as an option in techno-economic modeling. Print times were reduced from in excess of 3 days to under 30 hours, while surface quality and and tool integrity were substantially improved in the transition from the partially solid tool to the hollow tooling concept. Ability Composites produced autoclave-cured prepreg ducts on each of these tools. The autoclave conditions utilized were more severe than those of the initial trials, reaching temperatures of 160 °C and a pressure of 414 kPa. Under these conditions, the partially solid 3D printed tool with skin and 40% dense infill crushed significantly; however, the thicker ST-130 hollow tool showed good promise, deforming only slightly. The thinner hollow tool walls were unsuccessful as were the other materials. Overall, the hollow tool manufacturing process saved significant amounts of time and material in manufacturing as compared to the solid ducts and produced composite surface quality improvements compared to the traditional washout tooling. The TEM was developed to allow direct comparisons between conventional washout tool manufacturing processes and those developed at CSU. It also allowed for two separate 3D printed tool geometries to be analyzed and compared. In this case, the square bent duct tool geometry was determined to be representative of common washout tools. This geometry was compared with a scaled-up version of it to assess differences in the two manufacturing processes based on tool size. The model was developed to make use of user input in the form of geometry details, process steps, manufacturing parameters, bulk material costs, capital equipment costs, and general costs to calculate overall labor, material, capital equipment, and energy costs per manufactured tool for the conventional and additive manufacturing processes for the two representative geometries. It was also able to estimate step-by-step process times for the manufacturing process and geometries. Based on significant input from Ability Composites and CSU from their knowledge gained from hands-on manufacturing of the 3D printed bent duct tool geometry, costs and process times were calculated for the two manufacturing processes. Results showed that the additive manufacturing techniques developed at CSU can substantially reduce the costs of tool manufacturing by reducing labor times and material usage. This is because additive manufacturing is a relatively hands-off process and allows for the tool design to be optimized to reduce material usage. The disadvantage, however, is that process times for additive manufacturing are significantly longer. The three-dimensional (3D) printing process is slow if tight tolerances are required, but the analysis did show that print times could be reduced with the hollow tool geometry. Also, further advances in additive manufacturing could expedite the process. Costs and process times for the tool washout process were calculated separately. They showed that costs are relatively insignificant when compared to the overall tool manufacturing processes, but with increases in tool size, costs for the conventional manufacturing approach are larger than for additive manufacturing. Again, the washout process for conventional tools is very hands-on, whereas for additively manufactured tools the print medium is dissolved in an automated detergent bath at the sacrifice of process time. The analysis showed that optimizing the additively manufactured tools may also reduce washout times. Overall, Project 4.9 demonstrated that commercially available dissolvable 3D printing materials exist that can be used to produce dissolvable tooling capable of surviving prepreg composites fabrication under autoclave conditions of 121 °C (250 °F) and 345 kPa (50 psi). An alternative hollow dissolvable tool design was developed which was structurally superior to the initial concept and was cost and time effective versus conventional washout tooling. The 3D printed sacrificial tool required no added surface sealing steps prior to composite part layup and cure, offering a significant advantage over the porous conventional washout tooling.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Integral Twist Actuation of Helicopter Rotor Blades for Vibration Reduction

Active integral twist control for vibration reduction of helicopter rotors during forward flight is investigated. The twist deformation is obtained using embedded anisotropic piezocomposite actuators. An analytical framework is developed to examine integrally-twisted blades and their aeroelastic response during different flight conditions: frequency domain analysis for hover, and time domain analysis for forward flight. Both stem from the same three-dimensional electroelastic beam formulation with geometrical-exactness, and axe coupled with a finite-state dynamic inflow aerodynamics model. A prototype Active Twist Rotor blade was designed with this framework using Active Fiber Composites as the actuator. The ATR prototype blade was successfully tested under non-rotating conditions. Hover testing was conducted to evaluate structural integrity and dynamic response. In both conditions, a very good correlation was obtained against the analysis. Finally, a four-bladed ATR system is built and tested to demonstrate its concept in forward flight. This experiment was conducted at NASA Langley T~ansonic Dynamics Tunnel and represents the first-of-a-kind Mach-scaled fully-active-twist rotor system to undergo forward flight test. In parallel, the impact upon the fixed- and rotating-system loads is estimated by the analysis. While discrepancies are found in the amplitude of the loads under actuation, the predicted trend of load variation with respect to its control phase correlates well. It was also shown, both experimentally and numerically, that the ATR blade design has the potential for hub vibratory load reduction of up to 90% using individual blade control actuation. Using the numerical framework, system identification is performed to estimate the harmonic transfer functions. The linear time-periodic system can be represented by a linear time-invariant system under the three modes of blade actuation: collective, longitudinal cyclic, and lateral cyclic. A vibration minimizing controller is designed based on this result, which implements classical disturbance rejection algorithm with some modifications. The controller is simulated numerically, and more than 90% of the 4P hub vibratory load is eliminated. By accomplishing the experimental and analytical steps described in this thesis, the present concept is found to be a viable candidate for future generation low-vibration helicopters. Also, the analytical framework is shown to be very appropriate for exploring active blade designs, aeroelastic behavior prediction, and as simulation tool for closed-loop controllers.

Shin, SangJoon↗
DOE OSTI2022

Micromechanics of Fracture Propagation During Multistage Stress Relaxation and Creep in Brittle Rocks

Time-dependent rock deformation caused by the initiation and growth of fractures leads to the weakening of the rock mass. Understanding the fracturing mechanisms involved in the time-dependent behavior in brittle rocks is very important and to achieve this goal, a systematic series of three types of experiments was performed on double-flawed prismatic Barre granite specimens under unconfined compression. The first series aimed to identify the failure mechanism in the short-term failure mode under monotonic loading, the, second series involved multistage relaxation (constant strain) experiments to analyze the damage at different strain levels, and the third series explored the fracture propagation under multistage creep (constant load) experiments. The spatial and temporal evolution of cracking mechanisms were evaluated using the acoustic emission (AE) and two-dimensional digital image correlation (2D-DIC) techniques to observe the whole crack growth process as well as the accumulated inelastic strain at the specified region of interest. Results suggest that in the case of multistage creep experiments, the time to failure was less compared to the multistage relaxation, when loaded above the crack damage threshold (CD) estimated from the monotonic testing. Furthermore, the frequency magnitude distribution of the AE events generated in the three loading conditions followed the Gutenberg Richter model. A relatively lower b-value was obtained for the creep experiments, indicative of high energy AE events and faster crack growth. In addition, the AE and DIC results also revealed high evolution of tensile cracks at-different stages of creep and relaxation compared to shear and mixed-mode cracks.

58 GEOSCIENCES↗
DOE OSTI2022

Elucidating the challenges in extracting ultra-slow flame speeds in a closed vessel - A CH 2 F 2 microgravity case study using optical and pressure-rise data

Refrigerants with a low global warming potential (GWP) possess mild flammability. Hence, a fundamental understanding of their combustion characteristics is required to assess their fire-hazardous potential. The laminar burning velocity is one fundamental property to describe these under safety evaluation aspects. However, measuring laminar burning velocities for slow-burning components, such as low-GWP alternatives, is challenging. This is due to influencing effects such as buoyant flame deformation, stretch, radiation, and confinement, especially at the flammability limits. Here, in the present study, we investigated near-limit flames of the representative refrigerant difluoromethane (CH 2 F 2 ) with nitrogen-enriched oxidizer mixtures to elucidate the potentials and limitations at ultra-slow flame speeds of two widely used flame speed measurement methods in the closed vessel configuration: the optical flame speed measurement in the early quasi-isobaric regime and the flame speed determination from the pressure-rise during the later phase of the isochoric combustion process. Experiments were performed under terrestrial- and microgravity, conducted in a specially designed setup suitable for drop tower applications, using both methods in parallel. Recommendations for the extraction of flame speed data are derived from the non-buoyant reference investigations under microgravity and transferred to ultra-slow flames at terrestrial gravity. Near-limit flames under microgravity were found to be significantly affected by stretch effects for the optical method so that flame speed extrapolations to unstretched conditions are prone to errors. Stretch also affects the pressure method’s lower evaluation limit, and errors can be estimated by combining flame speed data from the pressure method and Markstein lengths from optical experiments.

42 ENGINEERING↗
DOE OSTI2020

A numerical study of the meso-structure variability in the compaction process of prepreg platelet molded composites

Compaction deformations in a thermoplastic composite system molded from chopped prepreg platelets were analyzed to obtain statistical distributions of local distortions that define the final consolidated meso-structure. This study provides an insight into and allows quantification of the intrinsic geometric variability of the composite meso-structure that controls the variability of the effective composite mechanical properties. A two-phase viscoelastic finite element model was utilized to treat the behavior of molten prepreg platelets, which were assumed to be incompressible and inextensible in the fiber direction during processing. A two-step simulation procedure with discretely modeled platelets and platelet-to-platelet contacts is proposed. The analysis begins with a drop simulation of the uncontrolled material deposition followed by a 3D finite element analysis of the preform compaction and bulk material flow. The proposed model provides a method for generating a composite meso-structure ab initio starting from a statistical sampling of the initial conditions for platelet deposition. The predicted platelet thickness distribution, intra-platelet fiber waviness and out-of-plane platelet undulations were found to be in good agreement with experimental measurements. An estimate of the volume of resin pockets was inferred from the compaction simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

$S$ Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

The orientation of $S$ Hmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed $S$ Hmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

58 GEOSCIENCES↗
DOE OSTI2022

Validation of material models for puncture of 7075-T651 aluminum plate

Plate puncture simulations are challenging computational tasks that require advanced material models including high strain rate and thermal-mechanical effects on both deformation and failure, plus finite element techniques capable of representing large deformations and material failure. The focus of this work is on the material issues, which require large sets of experiments, flexible material models and challenging calibration procedures. In this study, we consider the puncture of 12.7 mm thick, 7075-T651 aluminum alloy plates by a cylindrical punch with a hemispherical nose and diameter of 12.7 mm. The plasticity and ductile failure models were isotropic with calibration data obtained from uniaxial tension tests at different temperatures and strain rates plus quasi-static notched tension tests and shear-dominated tests described here. Sixteen puncture experiments were conducted to identify the threshold penetration energy, mode of puncture and punch acceleration during impact, The punch was mounted on a 139 kg mass and dropped on the plates with different impact speeds. Since the mass was the same in all tests, the quantity of interest was the impact speed. The axis and velocity of the punch were perpendicular to the plate surface. The mean threshold punch speed was 3.05 m/s, and the mode of failure was plugging by thermal-mechanical shear banding accompanied by scabbing fragments. Application of the material models in simulations of the tests yielded accurate estimates of the threshold puncture speed and of the mode of failure. Time histories of the punch acceleration compared well between simulation and test. Remarkably, the success of the simulations occurred in spite of even the smallest element used being larger than the width of the shear bands.

36 MATERIALS SCIENCE↗
NASA NTRS2018

Optical Performance Prediction of the Thirty Meter Telescope after Initial Alignment Using Optical Modeling

We present an estimate of the optical performance of the Thirty Meter Telescope (TMT) after execution of the full telescope alignment plan. The TMT alignment is performed by the Global Metrology System (GMS) and the Alignment and Phasing System (APS). The GMS first measures the locations of the telescope optics and instruments as a function of elevation angle. These initial measurements will be used to adjust the optics positions and build initial elevation look-up tables. Then the telescope is aligned using starlight as the input for the APS at multiple elevation angles. APS measurements are used to refine the telescope alignment to build elevation and temperature dependent look-up tables. Due to the number of degrees of freedom in the telescope (over 10,000), the ability of the primary mirror to correct aberrations on other optics, the tight optical performance requirements and the multiple instrument locations, it is challenging to develop, test and validate these alignment procedures. In this paper, we consider several GMS and APS operational scenarios. We apply the alignment procedures to the model-generated TMT, which consists of various quasi-static errors such as polishing errors, passive supports errors, thermal and gravity deformations and installation position errors. Using an integrated optical model and Monte-Carlo framework, we evaluate the TMT’s aligned states using optical performance metrics at multiple instrument and field of view locations. The optical performance metrics include the Normalized Point Source Sensitivity (PSSN), RMS wavefront error before and after Adaptive Optics (AO) correction, pupil position change, and plate scale distortion.

Rogers, John↗
NASA NTRS1994

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

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↗
NASA NTRS2012

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

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.↗
DOE OSTI2021

Energy based fracture initiation criterion for strain-crystallizing rubber-like materials with pre-existing cracks

Fracture prediction is indispensable for polymers, like rubbers, which have a broad range of applications mainly due to their high extensibility. The phenomenon known as strain-induced crystallization further contributes to the fracture toughness of certain rubbers. Here in this study, a criterion based on internal bond energy, incorporating the effects of crystallization, is proposed to predict fracture initiation in rubber-like materials with pre-existing cracks. First, a multi-scale mechanical model is developed for characterizing the behavior of rubber when subjected to both uniaxial and biaxial deformation states. At the microscale, both the amorphous and crystalline chain segments are modeled as elastic in order to consider the energy contribution by the molecular bond distortions. This internal energy is considered along with the entropic and crystalline free energy for each chain. In the chain model, the effects of loading condition and the relative orientation of a chain on its crystallinity are taken into account. At the macroscopic scale, an existing crystallinity distribution function is adapted and a mixed finite element formulation with an augmented Lagrangian multiplier is utilized to impose the incompressibility constraint. A non-affine maximal advance path constraint based homogenization model is utilized for bridging the two scales. Its potential to account for anisotropy in the stretched network compels the model to be preferable due to its physical significance, for the purpose of fracture modeling. The rigidity of the crystallites is accounted for by proposing a crystallite distortion energy in addition to the critical bond dissociation energy, for fracture initiation to occur. The model is validated by comparison with existing experimental results for both crystallizing and non-crystallizing rubbers. In addition to its potential to predict the material behavior when subjected to uniaxial and biaxial loading, the capability of the model to quantitatively estimate the effect of crystallization on fracture initiation is also verified.

36 MATERIALS SCIENCE↗
NASA NTRS2009

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↗
NASA NTRS2017

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.↗
DOE OSTI2023

Effects of anisotropic shale creep on the stress and permeability evolution of a geological nuclear waste repository

To ensure long-term safety and performance, geological nuclear waste repositories require low-permeability barriers such as bentonite buffers and/or shale host rock. Shale is not only known for its low permeability but also for its trend to undergo time-dependent deformation (i.e., creep), which could heal damage, but the effects of shale creep on the long-term performance of nuclear waste repositories have not been clearly understood. In particular, the anisotropic nature of shale (i.e., bedding) could have a significant effect on its creep behavior, and consequently, on the long-term performance of nuclear waste repositories. In this research, numerical simulations were carried out with the objective of showing the effects of anisotropic shale creep on the stress and permeability evolution of a generic geological nuclear waste repository in shale. The TOUGH-FLAC simulator, a thermo-hydromechanically (THM) coupled numerical code, was used for the simulations. To achieve the objective, comparisons were performed between the results of anisotropic shale creep simulations and those of different simulation cases, namely, no creep (i.e., elastic), isotropic creep, and long-term creep shale cases. Results of the comparisons show that the elastic and isotropic creep shale cases respectively led to the overestimation and underestimation of stress and permeability in the repository, whereas the long-term creep shale case, which accumulated greater creep in later periods than in earlier periods, helped to keep large shear and tensile stresses from developing while maintaining compressive spherical stress, resulting in consistently low permeability levels. These results indicate that performance assessments with elastic and isotropic creep formation models will provide the upper and lower bound estimates of stress and permeability, while more reasonable estimates will be given by an anisotropic creep formation model, and that shale with long-term creep characteristics will be beneficial in many aspects of the safety and performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1994

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.↗
DOE OSTI2022

Properties of the neutron star crust: Quantifying and correlating uncertainties with improved nuclear physics

A compressible liquid-drop model (CLDM) is used to correlate uncertainties associated with the properties of the neutron star (NS) crust with theoretical estimates of the uncertainties associated with the equation of state (EOS) of homogeneous neutron and nuclear matter. For the latter, we employ recent calculations based on Hamiltonians constructed using chiral effective-field theory (χEFT). Fits to experimental nuclear masses are employed to constrain the CLDM further, and we find that they disfavor some of the χEFT Hamiltonians. The CLDM allows us to study the complex interplay between bulk, surface, curvature, and Coulomb contributions, and their impact on the NS crust. It also reveals how the curvature energy alters the correlation between the surface energy and the bulk symmetry energy. Here our analysis quantifies how the uncertainties associated with the EOS of homogeneous matter implies significant uncertainties for the composition of the crust, its proton fraction, and the volume fraction occupied by nuclei. We find that the finite-size effects impact the crust composition but have a negligible effect on the net isospin asymmetry of matter. The isospin asymmetry is largely determined by the bulk properties and the isospin dependence of the surface energy. The most significant uncertainties associated with matter properties in the densest regions of the crust, the precise location of the crust-core transition, are found to be strongly correlated with uncertainties associated with the Hamiltonians. By adopting a unified model to describe the crust and the core of NSs, we tighten the correlation between their global properties such as their mass-radius relationship, moment of inertia, crust thickness, and tidal deformability with uncertainties associated with the nuclear Hamiltonians.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗