SEARCH · Search NASA
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.
Regional Source-Type Discrimination Using Nonlinear Alignment Algorithms
The discrimination problem in seismology aims to accurately classify different underground source types based on local, regional, and/or teleseismic observations of ground motion. Typical discriminant approaches are rooted in fundamental, physics-based differences in radiation pattern or wave excitation, which can be frequency-dependent and may not make use of the full waveform. In this article, we explore whether phase and amplitude distances derived from dynamic time warping (DTW) and elastic shape analysis (ESA) can inform event discrimination. We demonstrate the ability to distinguish underground point sources using synthetic waveforms calculated for a 1D Earth model and various source mechanisms. We then apply the method to recorded data from events in the Korean Peninsula, which includes declared nuclear explosions, a collapse event, and naturally occurring earthquakes. Phase and amplitude distances derived from DTW and ESA are then used to classify the event types via dendrogram and k-nearest-neighbor clustering analyses. Using information from the full waveform, we show how different underground sources can be distinguished at regional distances. We highlight the potential of these nonlinear alignment algorithms for discrimination and comment on ways we can extend the framework presented here.
Assessing Dynamic Time Warping Techniques for Discriminating Seismic Sources at Local and Regional Distances
Effective monitoring of seismic explosions and hazard assessment relies heavily on the accurate discrimination of underground seismic sources. This study investigates the application of novel nonlinear alignment techniques, specifically Dynamic Time Warping (DTW), for event-type discrimination at regional and local distances. Building on prior research that used DTW and Elastic Shape Analysis (ESA) in discrimination at regional distances, we evaluate the performance of recently developed variants of DTW, including a method that employs Pearson cross-correlation as a measure of warping distance and a time distortion coefficient that quantifies the type and degree of time distortion between signals. By analyzing observational datasets that include different source types, we assess the performance of these approaches for realistic monitoring scenarios. Specifically, we consider a dataset recorded at regional distances in the Korean Peninsula and a local-distance subset from the Unconstrained Utah Event Bulletin catalog to evaluate DTW-based discrimination across multiple distance scales. Additionally, we introduce the maximum cross-correlations of warped waveforms as a similarity metric for event classification. Through hierarchical cluster analysis and dendrogram interpretation, we present our findings, highlighting the strengths and limitations of these techniques in seismic event classification.
Explainable Machine Learning for Functional Data
Black-box machine learning models are recognized as useful tools for prediction applications, but the algorithmic complexity of some models causes interpretation challenges. Explainability methods have been proposed to provide insight into these models, but there is little research focused on supervised modeling with functional data inputs. We argue that, especially in applications of high consequence, it is important to explicitly model the functional dependence in a black-box analysis to not obscure or misrepresent patterns in explanations. As such, we propose the V ariable importance E xplainable E lastic S hape A nalysis (VEESA) pipeline for training supervised machine learning models with functional inputs. The pipeline is an analysis process that includes the data preprocessing, modeling, and post-hoc explanations. The preprocessing is done using elastic functional principal components analysis, which accounts for vertical and horizontal variability in functional data and, ultimately, allows for explanations in the original data space that identify the important functional variability without bias due to correlated variables. Here, we demonstrate the pipeline on two high-consequence applications: explosives classification for national security and inkjet printer identification in forensic science. The applications exhibit the VEESA pipeline’s ability to provide an understanding of the characteristics of the functional data useful for prediction. Code for implementing the pipeline is available in the veesa R package (and supplemental python code).
Shape design sensitivity analysis and optimization of three dimensional elastic solids using geometric modeling and automatic regridding
An automatic regridding method and a three dimensional shape design parameterization technique were constructed and integrated into a unified theory of shape design sensitivity analysis. An algorithm was developed for general shape design sensitivity analysis of three dimensional eleastic solids. Numerical implementation of this shape design sensitivity analysis method was carried out using the finite element code ANSYS. The unified theory of shape design sensitivity analysis uses the material derivative of continuum mechanics with a design velocity field that represents shape change effects over the structural design. Automatic regridding methods were developed by generating a domain velocity field with boundary displacement method. Shape design parameterization for three dimensional surface design problems was illustrated using a Bezier surface with boundary perturbations that depend linearly on the perturbation of design parameters. A linearization method of optimization, LINRM, was used to obtain optimum shapes. Three examples from different engineering disciplines were investigated to demonstrate the accuracy and versatility of this shape design sensitivity analysis method.
Three-dimensional finite strip analysis of elastic solids
Finite element shape functions were combined with beam eigenfunctions to formulate three-dimensional finite strips. Because of the orthogonality of the beam functions, three-dimensional problems were reduced to a series of two-dimensional problems, often with stiffness matrices of very narrow bandwidth. Isoparametric and high order finite element shape functions were used in the formulation of the 3-D finite strips. Numerical examples such as the static and free vibration analyses of simply supported thick plates are presented. Results were compared with existing solutions and good agreements were obtained in all cases.
Sensitivity analysis of a wing aeroelastic response
The present shape-sensitivity analysis of wing aeroelastic response proceeds from aeroelastic response sensitivities obtained on the basis of the aerodynamic performance valid for high aspect ratio wings in subsonic, subcritical flow. Attention is given to the shape sensitivity of various static aeroelastic responses; the formulation is general, and assumes that, for a given shape and elastic deformation, the aerodynamic analysis will furnish the distribution of the pressure and the pressure sensitivity derivatives with respect to the shape parameters of interest. Wing displacements are obtained by means of an iterative scheme.
Investigating Some Technical Issues on Cohesive Zone Modeling of Fracture
This study investigates some technical issues related to the use of cohesive zone models (CZMs) in modeling fracture processes. These issues include: why cohesive laws of different shapes can produce similar fracture predictions; under what conditions CZM predictions have a high degree of agreement with linear elastic fracture mechanics (LEFM) analysis results; when the shape of cohesive laws becomes important in the fracture predictions; and why the opening profile along the cohesive zone length needs to be accurately predicted. Two cohesive models were used in this study to address these technical issues. They are the linear softening cohesive model and the Dugdale perfectly plastic cohesive model. Each cohesive model constitutes five cohesive laws of different maximum tractions. All cohesive laws have the same cohesive work rate (CWR) which is defined by the area under the traction-separation curve. The effects of the maximum traction on the cohesive zone length and the critical remote applied stress are investigated for both models. For a CZM to predict a fracture load similar to that obtained by an LEFM analysis, the cohesive zone length needs to be much smaller than the crack length, which reflects the small scale yielding condition requirement for LEFM analysis to be valid. For large-scale cohesive zone cases, the predicted critical remote applied stresses depend on the shape of cohesive models used and can significantly deviate from LEFM results. Furthermore, this study also reveals the importance of accurately predicting the cohesive zone profile in determining the critical remote applied load.
An experimental-theoretical study of free vibrations of plates on elastic point supports
A theoretical and experimental study is made to investigate the effect on plate vibrations of varying the stiffness of corner elastic point supports. A theoretical model is developed using a Rayleigh-Ritz analysis which approximates the plate mode shapes as products of free-free beam modes. The elastic point supports are modelled both as massless translational springs, and springs with tip masses. The tip masses are included to better represent the experimental supports. An experiment is constructed using the bending stiffness of horizontal beams to support a square plate at its four corners. The stiffness of these supports can be varied over such a range that the plate fundamental frequency is lowered to 40% of the rigid support frequency. The variation with support stiffness of the frequencies of the first eight plate modes is measured, and compared with the theoretical results. The plate mode shapes for rigid supports are analyzed using holographic interferometry. There is excellent agreement between the theoretical and experimental results, except for high plate modes where the theoretical model is demonstrated to be inadequate.
A life comparison of tube and channel cooling passages for thrust chambers
The life analysis used to compare copper tubes and milled copper channels for rocket engine cooling passages is described. Copper tubes were chosen as a possible replacement for the existing milled copper channel configuration because (1) they offer increased surface area for additional enthalpy extraction; (2) they have ideal pressure vessel characteristics; and (3) the shape of the tube is believed to allow free expansion, thus accommodating the strain resulting from thermal expansion. The analysis was a two-dimensional elastic-plastic comparison, using a finite element method, to illustrate that, under the same thermal and pressure loading, the compliant shape of the tube increases the life of the chamber. The analysis indicates that for a hot-gas-side-wall temperature of 100 F the critical strain decreases from 1.25 percent in the channel to 0.94 percent in the tube. Since the life of rocket thrust chambers is most often limited by cyclic strain or strain range, this decrease corresponds to an expected tube life which is nearly twice the channel life.
Stress analysis of spacecraft parachutes using finite elements and large deformation theory
Spacecraft parachute stress analysis, using finite elements with nonlinear elastic properties to obtain shape and load distribution
BEST3D user's manual: Boundary Element Solution Technology, 3-Dimensional Version 3.0
The theoretical basis and programming strategy utilized in the construction of the computer program BEST3D (boundary element solution technology - three dimensional) and detailed input instructions are provided for the use of the program. An extensive set of test cases and sample problems is included in the manual and is also available for distribution with the program. The BEST3D program was developed under the 3-D Inelastic Analysis Methods for Hot Section Components contract (NAS3-23697). The overall objective of this program was the development of new computer programs allowing more accurate and efficient three-dimensional thermal and stress analysis of hot section components, i.e., combustor liners, turbine blades, and turbine vanes. The BEST3D program allows both linear and nonlinear analysis of static and quasi-static elastic problems and transient dynamic analysis for elastic problems. Calculation of elastic natural frequencies and mode shapes is also provided.
Improvements to the SADSAM computer program for aeroelasticity analysis
A computer program for calculating the elastic element forces and strain energy distribution for normal mode shapes is presented. The subjects discussed are: (1) revised aerodynamic force representation, (2) flap aerodynamics, (3) chordwise motions and aerodynamic forces, (4) wing-tail interaction, and (5) gust analysis. An aircraft configuration used in sample gust analysis is included.
Approximating the stress field within the unit cell of a fabric reinforced composite using replacement elements
This report concerns the prediction of the elastic moduli and the internal stresses within the unit cell of a fabric reinforced composite. In the proposed analysis no restrictions or assumptions are necessary concerning yarn or tow cross-sectional shapes or paths through the unit cell but the unit cell itself must be a right hexagonal parallelepiped. All the unit cell dimensions are assumed to be small with respect to the thickness of the composite structure that it models. The finite element analysis of a unit cell is usually complicated by the mesh generation problems and the non-standard, adjacent-cell boundary conditions. This analysis avoids these problems through the use of preprogrammed boundary conditions and replacement materials (or elements). With replacement elements it is not necessary to match all the constitutional material interfaces with finite element boundaries. Simple brick-shaped elements can be used to model the unit cell structure. The analysis predicts the elastic constants and the average stresses within each constituent material of each brick element. The application and results of this analysis are demonstrated through several example problems which include a number of composite microstructures.
Orientation effects on the measurement and analysis of critical CTOA in an aluminum alloy sheet
Fracture tests were conducted on 76.2mm wide, 2.3mm thick middle crack tension (M(T)) specimens machined from 2024-T3 aluminum sheet. The specimens were tested on the T-L orientation and comparisons were made to similar tests conducted in the L-T orientation. Measurement of critical crack tip opening angle (CTOA), applied stress, and crack front shape were made as a function of crack extension. A two-dimensional, elastic-plastic finite element analysis was used to simulate the fracture behavior for both orientations. The results indicate that the T-L orientation had a 10 percent lower stress at fracture than similar tests conducted in the L-T orientation. Correspondingly, the critical CTOA in the T-L tests reached a constant value of 4.7 degrees after 2-3mm of crack extension and the L-T tests reached a value of 6 degrees. The fracture surfaces of the T-L specimens were observed to remain flat, while those of the L-T specimens transitioned to a 45 degree slant fracture after about 2-3mm of crack extension. The tunneling behavior of the two orientations also differed; the T-L specimens reached a deeply tunneled stabilized crack front shape while, the L-T specimens were observed to have only a small amount of tunneling once the crack began to grow on the 45 degree slant. The two-dimensional, elastic-plastic finite element analysis was able to simulate the fracture behavior for both the T-L and L-T orientations.
Mesh Deformation Based on Fully Stressed Design: The Method and Two-Dimensional Examples
Mesh deformation in response to redefined boundary geometry is a frequently encountered task in shape optimization and analysis of fluid-structure interaction. We propose a simple and concise method for deforming meshes defined with three-node triangular or four-node tetrahedral elements. The mesh deformation method is suitable for large boundary movement. The approach requires two consecutive linear elastic finite-element analyses of an isotropic continuum using a prescribed displacement at the mesh boundaries. The first analysis is performed with homogeneous elastic property and the second with inhomogeneous elastic property. The fully stressed design is employed with a vanishing Poisson s ratio and a proposed form of equivalent strain (modified Tresca equivalent strain) to calculate, from the strain result of the first analysis, the element-specific Young s modulus for the second analysis. The theoretical aspect of the proposed method, its convenient numerical implementation using a typical linear elastic finite-element code in conjunction with very minor extra coding for data processing, and results for examples of large deformation of two-dimensional meshes are presented in this paper. KEY WORDS: Mesh deformation, shape optimization, fluid-structure interaction, fully stressed design, finite-element analysis, linear elasticity, strain failure, equivalent strain, Tresca failure criterion
Vibration analysis of a 1/15 scale dynamic model of a space shuttle configuration
The natural frequencies and mode shapes of a 1/15 scale space shuttle dynamics model are analytically determined. The model, a parallel beam type structure with delta wings, is dynamically representative of the stiffness and mass properties of an early space shuttle design. Important characteristics of the model are elastic interfaces with adjustable spring rates. Normal mode computations are made using the finite element modeling technique as implemented in the NASTRAN (NASA Structural Analysis) computer program. The feasibility of neglecting elastic deformations in the lower modes was investigated using a rigid body model. Using NASTRAN, natural frequencies and mode shapes were first calculated for the booster fuselage, orbiter fuselage, and both delta wings in a free-free condition. Next, the fuselages were connected for various spring rates. Then the wings were attached to each fuselage, and the booster and orbiter were analyzed as separate airplanes. Finally, the two airplanes were elastically joined, and the complete model was analyzed.
Application of elastic and elastic-plastic fracture mechanics methods to surface flaws
Fuel tanks that are a part of the External Tank assembly for the Space Shuttle are made of relatively thin 2219-T87 aluminum plate. These tanks contain about 917 m of fusion weld seam, all of which is nondestructively inspected for flaws and all those found are repaired. The tanks are subsequently proof-tested to a pressure that is sufficiently severe to cause weld metal yielding in a few local regions of the weld seam. The work undertaken in the present project was to develop a capability to predict flaw growth from undetected surface flaws that are assumed to be located in the highly stressed regions. The technical challenge was to develop R-curve prediction capability for surface cracks in specimens that contain the flaws of unusual sizes and shapes deemed to be of interest. The test techniques developed and the elastic-plastic analysis concepts adopted are presented. The flaws of interest were quite small surface cracks that were narrow-deep ellipses that served to exacerbate the technical difficulties involved.