Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Elastic Shape Analysis”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The physics potential of a reactor neutrino experiment with Skipper CCDs: Measuring the weak mixing angle

We analyze in detail the physics potential of an experiment like the one recently proposed by the vIOLETA collaboration: a kilogram-scale Skipper CCD detector deployed 12 meters away from a commercial nuclear reactor core. This experiment would be able to detect coherent elastic neutrino nucleus scattering from reactor neutrinos, capitalizing on the exceptionally low ionization energy threshold of Skipper CCDs. To estimate the physics reach, we elect the measurement of the weak mixing angle as a case study. We choose a realistic benchmark experimental setup and perform variations on this benchmark to understand the role of quenching factor and its systematic uncertainties, background rate and spectral shape, total exposure, and reactor antineutrino flux uncertainty. We take full advantage of the reactor flux measurement of the Daya Bay collaboration to perform a data driven analysis which is, up to a certain extent, independent of the theoretical un- certainties on the reactor antineutrino flux. We show that, under reasonable assumptions, this experimental setup may provide a competitive measurement of the weak mixing angle at few MeV scale with neutrino-nucleus scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1991

Advanced development of the boundary element method for elastic and inelastic thermal stress analysis

The focus of this dissertation is on advanced development of the boundary element method for elastic and inelastic thermal stress analysis. New formulations for the treatment of body forces and nonlinear effects are derived. These formulations, which are based on particular integral theory, eliminate the need for volume integrals or extra surface integrals to account for these effects. The formulations are presented for axisymmetric, two and three dimensional analysis. Also in this dissertation, two dimensional and axisymmetric formulations for elastic and inelastic, inhomogeneous stress analysis are introduced. The derivatives account for inhomogeneities due to spatially dependent material parameters, and thermally induced inhomogeneities. The nonlinear formulation of the present work are based on an incremental initial stress approach. Two inelastic solutions algorithms are implemented: an iterative; and a variable stiffness type approach. The Von Mises yield criterion with variable hardening and the associated flow rules are adopted in these algorithms. All formulations are implemented in a general purpose, multi-region computer code with the capability of local definition of boundary conditions. Quadratic, isoparametric shape functions are used to model the geometry and field variables of the boundary (and domain) of the problem. The multi-region implementation permits a body to be modeled in substructured parts, thus dramatically reducing the cost of analysis. Furthermore, it allows a body consisting of regions of different (homogeneous) material to be studied. To test the program, results obtained for simple test cases are checked against their analytic solutions. Thereafter, a range of problems of practical interest are analyzed. In addition to displacement and traction loads, problems with body forces due to self-weight, centrifugal, and thermal loads are considered.

Henry, Donald P., Jr.↗
NASA NTRSDate not supplied

Multiscale Mesh Adaptation for Transonic Aeroelastic Flutter Problems

This work applies multiscale mesh adaptation with refine to reduce spatial discretization error of aeroelastic computational fluid dynamics (CFD) simulations. Benchmark flutter models, such as the pitch and plunge NACA64A-010 airfoil and the benchmark supercritical wing, are studied with both a linearized frequency-domain solver and time-marching CFD coupled to a modal structural solver in FUN3D. The undeformed NASA Common Research Model (CRM), an aeroelastic jig shape variant of the CRM, is also studied with the linearized frequency-domain approach. For these cases, the adaptation process converges to comparable flutter predictions to hand-generated meshes but with smaller node counts. However the additional disciplines of the linearized frequency-domain analysis, the mesh deformation, and the unsteady finite-volume solver create robustness challenges that need to be addressed before it can be applied as a fully automated process for complex transonic aeroelastic problems. In particular, negative volumes are observed to be an issue for FUN3D’s linear elasticity mesh deformation solver when moving the adapted meshes.

Aeroelasticity↗
NASA NTRS1974

Transverse vibrations of non-uniform beams

The free vibrations of elastic beams with nonuniform characteristics are investigated theoretically by a new method. The new method is seen to combine the advantages of a finite element approach and of a Rayleigh-Ritz analysis. Comparison with the known analytical results for uniform beams shows good convergence of the method for natural frequencies and modes. For internal shear forces and bending moments, the rate of convergence is less rapid. Results from experiments conducted with a cantilevered helicopter blade with strong nonuniformities and also from alternative theoretical methods, indicate that the theory adequately predicts natural frequencies and mode shapes. General guidelines for efficient use of the method are presented.

Klein, L.↗
NASA NTRS1988

Time-domain transient elastodynamic analysis of 3-D solids by BEM

The BEM algorithm developed by Banerjee and Ahmad (1985) and Banerjee et al. (1986) for the analysis of three-dimensional time-dependent elastodynamics problems is extended and refined to improve the accuracy of the internal stress predictions. A quadratic spatial shape function and a linear temporal shape function are employed, and the complete numerical implementation (as part of BEST3D) is outlined. Applications of the improved BEM to typical problems are presented graphically and briefly characterized. Included are spherical cavities subjected to sudden radial expansion, triangular or rectangular pulses of radial pressure, or ramp loading; a circular loaded area on a half space; and a flexible square plate foundation on an elastic half space.

Ahmad, S.↗
NASA NTRSDate not supplied

Porous Microstructure Analysis (PuMA) software

The Porous Microstructure Analysis (PuMA) software was developed to provide a robust and efficient framework for computing material properties based on their microstructures. The development was motivated by advancements in X-ray microtomography, an imaging technology that can resolve the structure of a material at a sub-micron scale, in 3D and even in 4D (over time). PuMA provides the capability of computing a comprehensive spectrum of properties, from the most fundamental geometric features of a microstructure, to advanced anisotropic thermo-elastic properties. In addition, the software can generate artificial microstructures, ranging from simple analytical shapes to complex fibrous woven and non- woven geometries, which can be used in performance optimization studies. This presentation will highlight many of the capabilities of the recent open-source release.

microtomography↗
NASA NTRS1987

Damage growth in composite laminates with interleaves

The influence of placing interleaves between fiber reinforced plies in multilayered composite laminates is investigated. The geometry of the composite is idealized as two dimensional, isotropic, linearly elastic media made of a damaged layer bonded between two half planes and separated by thin interleaves of low extensional and shear moduli. The damage in the layer is taken in the form of a symmetric crack perpendicular to the interface and may extend up to the interface. The case of an H-shaped crack in the form of a broken layer with delamination along the interface is also analyzed. The interleaves are modeled as distributed shear and tension springs. Fourier integral transform techniques are used to develop solutions in terms of singular integral equations. An asymptotic analysis of the integral equations based on Muskhelishvili's techniques reveals logarithmically singular axial stresses in the half plane at the crack tips for the broken layer. For the H shaped crack, similar singularities are found to exist in the axial stresses at the interface crack tips in the layer and the half plane. The solution of the equations is found numerically for the stresses and displacements by using the Hadamard's concept of direct differentiation of Cauchy integrals as well as Gaussian integration techniques.

Goree, James G.↗
DOE OSTI2020

Comparing event generator predictions and ab-initio calculations of $\nu$-$^{12}$C neutral current quasi-elastic scattering at 1 GeV

The measurement of neutrino oscillations and exotic physics searches are important parts of the physics program in the near future, with new state-of-the-art experiments planned within the next decade. Future and modern experiments in these fields will make use of nuclear targets. Event Generators (EGs) are software used in the analysis of neutrino oscillation experiments. EGs use to predict kinematic observables for a range of neutrino energies. These simulations may lack physics captured by more rigorous theoretical calculations. This work compares EG performance to nuclear theory calculations by comparing observables generated in the two frameworks. We provide a common set of definitions between theory and experiment and assess the physics contained in EG simulations. Neutral current quasi-elastic (NCQE) scattering events for neutrinos and anti-neutrinos on a $^{12}$C target are simulated with a specific EG, NEUT, used by the T2K experiment for its analysis. The simulated cross sections are compared to analytic calculations from nuclear theory within the factorization scheme. We compare the NEUT implementation of two different models on nuclear spectral functions: the Relativistic Fermi Gas (RFG) and the correlated basis spectral function (CBF) to analytic calculations of the same models in the factorization scheme. For both nuclear physics models, we compare the appearance of features in the distributions relevant to experimental analyses. Qualitatively, the shape of the simulated distribution is similar to the one obtained through theory calculations; however, there are some discrepancies between the theory calculations and the NEUT simulation. While the EG simulations and analytic calculations with the same model of nuclear dynamics show similar overall features, there are still differences between the two.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1981

The failure strengths of perfect diamond crystals

Finite elasticity analysis is extended to the 110 direction, where off axis strain symmetry is not present, and the third order elastic data are obtained for diamond. The compressive yield strengths of perfect diamond crystals loaded in the 100, 110, and 111 directions are predicted to be 2.2, 5.6, and 2.8 Mbars, respectively, while the corresponding tensile fracture strengths are 1.0, 0.5, and 0.5 Mbars. From these results and from Hertz theory it is predicted that ring fracture of spherically tipped diamonds pressed against a flat will occur at pressures of 1.8-1.9 Mbars, substantially below the yield pressure (above 3 Mbars). Modification of the tip shape leads to a predicted increase in the pressure at which fracture occurs.

Whitlock, J.↗
DOE OSTI2020

Software-Defined Ultrasonic Communication System Based on Time-reversal Signal Processing

Ultrasonic communication through solid channels suffers from channel dispersion, attenuation, and multipath effect. For example, in elbow shape pipes or complex structures, the elastic waves are reflected and reverberated in multipath resulting in a major dispersion of the signals. In solid channels such as flat bar, multipath propagation is the main reason for bit-wavelet spreading and intersymbol interference (ISI). The conventional modulation method like OOK (on-off-keying) or PSK (phase-shift-keying) fails to provide a useful Signal-to-Noise Ratio (SNR). Time-reversal (TR) signal processing method is a promising technique that can compress the spreading of the bit-wavelet and strengthen the peak signal energy on the receiver. In this study, Finite Element Analysis (FEA) using ABAQUS software tools is used to characterize the bit-wavelet spreading effect caused by dispersive communication channels. Also, this study presents the impact of TR ultrasonic communication as a function of excitation pulse width, channel length, and carrier frequency through a flat-bar channel. A Software-Defined System-on-Chip (SD-SoC) platform is examined as a viable solution for the low-complexity and energy-efficient ultrasonic communication system. The proposed TR based SD-SoC ultrasonic communication system can make full use of signal dispersion for peak performance by minimizing the ISI within a dispersive and reverberant metal channel.

FEM simulation↗
NASA NTRSDate not supplied

Test and Analysis of a Composite Conical-Cylindrical Shell

Launch-vehicle shell structures, which can be comprised of both cylindrical and conical sections, are known to be susceptible to buckling due to their large radius-to-thickness ratios. The advancements in composite manufacturing and numerical methods have enabled designers to consider more nontraditional shapes, such as connecting the conical and cylindrical sections with a toroidal transition to create a single-piece conical-cylindrical shell. This single-piece construction eliminates the need for a heavy interface ring between sections and has the potential to save mass. To better understand the buckling behavior, a composite conical-cylindrical shell was designed, fabricated, and tested. Prior to test, a finite element model that included thickness variations and radial imperfections was created. The test article buckled elastically at 251.8 kN, approximately 8.8% higher than the predicted buckling load of 231.4 kN Continued research in conical-cylindrical structures has the potential to expand the design space for launch-vehicle structures and lead to improved designs and reduced weight.

Buckling↗
DOE OSTI2021

Stress-induced phase transitions in nanoscale CuInP 2 S 6

Using the Landau-Devonshire approach and available experimental results, we constructed multiwell thermodynamic potential of the layered ferroelectric CuInP 2 S 6 (CIPS). The analysis of temperature dependences of the dielectric permittivity and lattice constants for different applied pressures unexpectedly reveals the critically important role of a nonlinear electrostriction in this material. With the nonlinear electrostriction included we calculated the temperature and pressure phase diagrams and spontaneous polarization of a bulk CIPS, within the assumed range of applicable temperatures and applied pressures. Using the developed thermodynamic potential, we revealed the strain-induced phase transitions in thin epitaxial СIPS films, as well as the stress-induced phase transitions in СIPS nanoparticles, the shape of which varies from prolate needles to oblate disks. Furthermore, w also revealed the strong influence of a mismatch strain, elastic stress, and shape anisotropy on the phase diagrams and polar properties of a nanoscale CIPS, and derived analytical expressions allowing for elastic control of the nanoscale CIPS polar properties. Hence obtained results can be of particular interest for the strain engineering of nanoscale layered ferroelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1981

An empirical stress-intensity factor equation for the surface crack

This paper presents an empirical stress-intensity factor equation for a surface crack as a function of parametric angle, crack depth, crack length, plate thickness and plate width for tension and bending loads. The stress-intensity factors used to develop the equation were obtained from a previous three-dimensional, finite-element analysis of semielliptical surface cracks in finite elastic plates subjected to tension or bending loads. A wide range of configuration parameters was included in the equation. The ratios of crack length to plate thickness and the ratios of crack depth to crack length ranged from 0 to 1.0. The effects of plate width on stress-intensity variations along the crack front were also included. The equation was used to predict patterns of surface-crack growth under tension or bending fatigue loads. The equation was also used to correlate surface-crack fracture data for a brittle epoxy material within + or - 10 percent for a wide range of crack shapes and crack sizes.

Newman, J. C., Jr.↗
NASA NTRS1987

Finite element models of wire rope for vibration analysis

The usefulness of wire rope in shock and vibration isolation is briefly reviewed and its modeling for the purpose of vibration analysis is addressed. A model of a nominally straight segment of wire rope is described in which the rope structure is represented by a maiden, or central, strand of wire with one (or more) strand(s) wrapped around it in a helix (helices). The individual strands are modeled using finite elements and MSC NASTRAN. Small linear segments of each wire are modeled mathematically by dividing them lengthwise into triangular prisms representing each prism by a solid NASTRAN element. To model pretensioning and allow for extraction of internal force information from the NASTRAN model, the wound strands are connected to the maiden strand and each other using spring (scalar elastic) elements. Mode shapes for a length of wire rope with one and fixed to a moving base and the other attached to a point mass, are presented. The use of the NASTRAN derived mode shapes to approximate internal normal forces in equations of motion for vibration analyses is considered.

Cochran, J. E., Jr.↗
NASA NTRS1987

Mechanisms of rolling contact spalling

The results of a study aimed at analyzing the mechanical material interactions responsible for rolling contact spalling of the 440 C steel, high pressure oxygen turbopump bearings are presented. A coupled temperature displacement finite element analysis of the effects of friction heating under the contact is presented. The contact is modelled as a stationary, heat generating, 2 dimensional indent in an elastic perfectly plastic half-space with heat fluxes up to 8.6 x 10000 KW/m sq comparable to those generated in the bearing. Local temperatures in excess of 1000 C are treated. The calculations reveal high levels of residual tension after the contact is unloaded and cools. Efforts to promote Mode 2/Mode 3 fatigue crack growth under cyclic torsion in hardened 440 C steel are described. Spalls produced on 440 C steel by a 3 ball/rod rolling contact testing machine were studied with scanning microscopy. The shapes of the cyclic, stress strain hysteresis loops displayed by hardened 440 C steel in cyclic torsion at room temperature are defined for the plastic strain amplitudes encountered in rolling/sliding contact. Results of these analyses are discussed in detail.

Kumar, A. M.↗
NASA NTRS2000

A Strip-Yield Model for Predicting the Growth of Part-Through Cracks Under Cyclic Loading

Flaws exist in aircraft structures due to manufacturing operations and material defects. Under variable amplitude cyclic loading, these flaws grow as part-through cracks reducing the residual strength of structural components. To meet damage tolerant design requirements, accurate flaw growth predictions are needed which account for continual changes in crack shape as well as crack growth retardation and acceleration. Predicting the growth of part-through cracks under cyclic loading using an innovative and computationally efficient model is the focus of the research summarized in this report. In this research effort, a slice synthesis methodology was developed and used to construct a modified strip-yield model for the part-through semi-elliptical surface flaw, enabling prediction of plasticity-induced closure along the crack front and subsequent fatigue crack growth under constant amplitude and variable amplitude loading. While modeling the plasticity-induced closure in a part-through flaw may be performed using three dimensional elastic-plastic finite element analysis, this type of effort is impractical from an engineering perspective. A modified strip-yield model similar to that used in FASTRAN for part-through flaws is a much needed engineering design tool, particularly when computational resources are limited.

Daniewicz, S. R.↗
DOE OSTI2025

Sustainable Shape Memory Elastomers with Reduced Melt Viscosity and Enhanced Stiffness

Melt reactive processing of lignin with nitrile rubber is a promising approach to synthesizing shape memory materials. The strong intramolecular interactions in lignin macromolecular structures, caused by π–π stacking in aromatic rings and hydrogen bonding, often result in large phase separation or low miscibility with rubbers. In this study, we investigated the chemical and molecular characteristics, as well as the stiffness and complex viscosity, of modified kraft lignin melt-reacted with an acrylonitrile/butadiene copolymer containing 41% acrylonitrile (NBR41). To enhance the macromolecular compatibility of kraft lignin with NBR41, kraft lignin was cross-linked with poly(propylene glycol) diglycidyl ether (PPDE) and trimethylolpropane triglycidyl ether (TTE), both rich in epoxy reactive groups capable of forming chemical bonds with hydroxyl and carboxyl groups. Here, our findings demonstrate that the modification of kraft lignin with PPDE and TTE resulted in significantly increased stiffness of the composites. The elastic modulus of NBR41-Kraft lignin-PPDE and NBR41-Kraft lignin-TTE increased by 82 and 162%, respectively. Both the yield strength and Young’s modulus of these two samples showed dramatic improvements. Specifically, the yield strength and Young’s modulus of NBR41-Kraft lignin-TTE increased nearly 4 and 3-fold, respectively, compared to the control sample. Interestingly, despite significant improvements in mechanical properties, the viscosity of NBR41-Kraft lignin-PPDE was substantially lower than that of the control sample. At 210 °C and an angular frequency of 1 rad/s, the complex viscosity of NBR41-Kraft lignin was approximately 100.25 ± 4.77 kPa·s, while that of NBR41-Kraft lignin-PPDE was significantly lower at 56 ± 0.93 kPa·s. These findings were validated through Fourier transform infrared spectroscopy, scanning electron microscopy, dynamic mechanical analysis, thermal characterization, rheological tests, and quasi-elastic neutron scattering techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

Preliminary report on tests of tensile specimens with a part-through surface notch for a filament wound graphite/epoxy material

The behavior of tensile coupons with surface notches of various semi-elliptical shapes were evaluated for specimens obtained from a filament wound graphite/epoxy cylinder. The quasi-static test results, in some instances, are inadequate for defining complete trend curves and the interpretive analysis is considered to be preliminary. Specimens with very shallow notches were observed to be notch insensitive and the unnotched strength from these specimens was determined to be 54.97 Ksi. The failure strain of the laminate was found to be 1.328%. Specimens with deeper notches were sensitive to notch depth, notch aspect ratio, and specimen width. Using the unnotched strength of 54.97 Ksi and Poe's general toughness parameter, the fracture toughness was estimated to be 27.2 Ksi square root of In. Isotropic linear elastic fracture mechanics together with the estimated fracture toughness correctly predicted the influence of notch depth, aspect ratio, and specimen finitewidth.

Harris, C. E.↗