Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear structural 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 1,333 records · Page 74

Simulating the non-monotonic strain response of nanoporous multiferroic composites under electric field control

In this work, we simulate and analyze the mechanical response of a class of multiferroic materials consisting of a templated porous nanostructure made out of cobalt ferrite (CFO) partially filled by atomic layer deposition (ALD) with a ferroelectric phase of lead zirconate titanate (PZT). The strain in the device is measured when an electric field is applied for varying ALD thicknesses, displaying a non-monotonic dependence with a maximum strain achieved for a coating thickness of 3 nm. To understand this behavior, we apply finite element modeling to the smallest repeatable unit of the nanoporous template and simulate the mechanical response as a function of PZT coating thickness. We find that this non-monotonic response is caused by the interplay between two driving forces opposing one another. First, increased porosity works toward increasing the strain due to a reduced system stiffness. Second, decreased porosity involves a larger mass fraction of PZT, which drives the electro-mechanical response of the structure, thus leading to a larger strain. The balance between these two driving forces is controlled by the shear coupling at the CFO/PZT interface and the effective PZT cross section along the direction of the applied electric field. Here, our numerical results show that considering a nonlinear piezoelectric response for PZT leads to an improved agreement with the experimental data, consistent with ex situ poling of the nanostructure prior to magnetic measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Parametric, Data-Driven, Non-Intrusive Reduced-Order Model Framework for Crystal Plasticity Simulations of Voids

The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).

Anisotropy↗

Modal Survey Test of the SOTV 2X3 Meter Off-Axis Inflatable Concentrator

NASA's Marshall Space Flight Center has had several projects involving inflatable space structures. Projects in solar thermal propulsion have had the most involvement, primarily inflatable concentrators. A flight project called Shooting Star Experiment initiated the first detailed design, analysis and testing effort involving an inflatable concentrator that supported a Fresnel lens. The lens was to concentrate the sun's rays to provide an extremely large heat transfer for an experimental solar propulsion engine. Since the conclusion of this experiment, research and development activities for solar propulsion at Marshall Space Flight Center have continued both in the solar propulsion engine technology as well as inflatable space structures. Experience gained in conducting modal survey tests of inflatable structures for the Shooting Star Experiment has been used by dynamic test engineers at Marshall Space Flight Center to conduct a modal survey test on a Solar Orbital Transfer Vehicle (SOTV) off-axis inflatable concentrator. This paper describes how both previously learned test methods and new test methods that address the unique test requirements for inflatable structures were used. Effects of the inherent nonlinear response of the inflatable concentrator on test methods and test results are noted as well. Nine analytical mode shapes were successfully correlated to test mode shapes. The paper concludes with several "lessons learned" applicable to future dynamics testing and shows how Marshall Space Flight Center has utilized traditional and new methods for modal survey testing of inflatable space structures.

Engberg, Robert C.↗

Giant nonlinear optical responses from photon-avalanching nanoparticles

Avalanche phenomena use steeply nonlinear dynamics to generate disproportionately large responses from small perturbations, and are found in a multitude of events and materials. Photon avalanching enables technologies such as optical phase-conjugate imaging, infrared quantum counting and efficient upconverted lasing. However, the photon-avalanching mechanism underlying these optical applications has been observed only in bulk materials and aggregates, limiting its utility and impact. In this study we report the realization of photon avalanching at room temperature in single nanostructures-small, Tm 3+ -doped upconverting nanocrystals-and demonstrate their use in super-resolution imaging in near-infrared spectral windows of maximal biological transparency. Avalanching nanoparticles (ANPs) can be pumped by continuous-wave lasers, and exhibit all of the defining features of photon avalanching, including clear excitation-power thresholds, exceptionally long rise time at threshold, and a dominant excited-state absorption that is more than 10,000 times larger than ground-state absorption. Beyond the avalanching threshold, ANP emission scales nonlinearly with the 26th power of the pump intensity, owing to induced positive optical feedback in each nanocrystal. This enables the experimental realization of photon-avalanche single-beam super-resolution imaging with sub-70-nanometre spatial resolution, achieved by using only simple scanning confocal microscopy and without any computational analysis. Pairing their steep nonlinearity with existing super-resolution techniques and computational methods, ANPs enable imaging with higher resolution and at excitation intensities about 100 times lower than other probes. The low photon-avalanching threshold and excellent photostability of ANPs also suggest their utility in a diverse array of applications, including sub-wavelength imaging and optical and environmental sensing.

36 MATERIALS SCIENCE↗

Nonparametric identification of a class of nonlinear close-coupled dynamic systems

A nonparametric identification technique for the identification of close coupled dynamic systems with arbitrary memoryless nonlinearities is presented. The method utilizes noisy recorded data (acceleration, velocity and displacement) to identify the restoring forces in the system. The masses in the system are assumed to be known (or fairly well estimated from the design drawings). The restoring forces are expanded in a series of orthogonal polnomials and the coefficients of these polynomial expansions are obtained by using least square fit method. A particularly simple and computationally efficient method is proposed for dealing with separable restoring forces. The identified results are found to be relatively insensitive to measurement noise. An analysis of the effects of measurement noise on the quality of the estimates is given. The computations are shown to be relatively quick (when compared say to the Wiener identification method) and the core storage required relatively small, making the method suitable for onboard identification of large space structures.

Udwadia, F. E.↗

Effects of Accident Thermal Loading on In-Plane Shear Behavior of Steel-Plate Composite Walls

Structural walls in safety-related nuclear facilities are required to be designed for seismic and accident thermal (due to postulated high-energy pipe break events) loading combination. Current U.S. and international codes provide limited guidance for analysis and design of walls for this loading combination. This paper describes the experimental results and observations from tests conducted on a laboratory-scale (1:4 to 1:5) test unit representing steel-plate composite (SC) walls subjected to combined in-plane (seismic) and accident thermal loading. The test unit was subjected to surface temperatures of up to 450°F in combination with cyclic in-plane loading. Results of similar experiments recently conducted in Japan are also summarized (with surface temperatures up to 570°F). Surface heating combined with the low thermal conductivity and high specific heat of concrete resulted in nonlinear thermal gradients through the thickness of the specimens. These nonlinear thermal gradients and the associated self- or internal restraint led to extensive concrete cracking. This concrete cracking reduced the initial and secant stiffness of the specimens. The initial stiffness of the heated specimens reduced to 30 to 40% of the initial stiffness of the control (unheated) specimen. The secant stiffness of the heated specimens reduced up to 50% of the secant stiffness of the control (unheated) specimen. However, the in-plane shear strength of the heated SC specimens was still approximately 10 to 30% greater than the nominal in-plane shear strength calculated, for the limit state of steel plate von Mises yielding, using AISC N690 equations and measured material properties. Evaluation of the experimental results and observations suggests that the in-plane shear strength of SC walls subjected to typical accident surface temperatures (up to 570°F) can be estimated conservatively using the current provisions of AISC N690. The stiffness for accident thermal loading combinations can be considered to reduce from cracked composite stiffness at ambient temperature to fully cracked—that is, steel-only stiffness as the surface temperature increases up to typical accident value.

Construction & Building Technology↗

Model-size reduction for the analysis of symmetric structures with asymmetric boundary conditions

A simple computational procedure is presented for reducing the size of the analysis model for a symmetric structure with asymmetric boundary conditions to that of the corresponding structure with symmetric boundary conditions. The procedure is based on approximating the asymmetric response of the structure by a linear combination of symmetric and antisymmetric global approximation vectors (or modes). The key elements of the procedure are (1) restructuring the governing finite-element equations to delineate the contributions to the symmetric and antisymmetric components of the asymmetric response, (2) successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite-element method is first used to generate a few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The effectiveness of the computational procedure is demonstrated by means of numerical examples of linear static problems of shells, and its potential for solving nonlinear problems is discussed.

Noor, Ahmed K.↗

Review of probabilistic analysis of dynamic response of systems with random parameters

The various methods that have been studied in the past to allow probabilistic analysis of dynamic response for systems with random parameters are reviewed. Dynamic response may have been obtained deterministically if the variations about the nominal values were small; however, for space structures which require precise pointing, the variations about the nominal values of the structural details and of the environmental conditions are too large to be considered as negligible. These uncertainties are accounted for in terms of probability distributions about their nominal values. The quantities of concern for describing the response of the structure includes displacements, velocities, and the distributions of natural frequencies. The exact statistical characterization of the response would yield joint probability distributions for the response variables. Since the random quantities will appear as coefficients, determining the exact distributions will be difficult at best. Thus, certain approximations will have to be made. A number of techniques that are available are discussed, even in the nonlinear case. The methods that are described were: (1) Liouville's equation; (2) perturbation methods; (3) mean square approximate systems; and (4) nonlinear systems with approximation by linear systems.

Kozin, F.↗

Orion Parachute Riser Cutter Development

This paper presents the tests and analytical approach used on the development of a steel riser cutter for the CEV Parachute Assembly System (CPAS) used on the Orion crew module. Figure 1 shows the riser cutter and the steel riser bundle which consists of six individual cables. Due to the highly compressed schedule, initial unavailability of the riser material and the Orion Forward Bay mechanical constraints, JSC primarily relied on a combination of internal ballistics analysis and LS-DYNA simulation for this project. Various one dimensional internal ballistics codes that use standard equation of state and conservation of energy have commonly used in the development of CAD devices for initial first order estimates and as an enhancement to the test program. While these codes are very accurate for propellant performance prediction, they usually lack a fully defined kinematic model for dynamic predictions. A simple piston device can easily and accurately be modeled using an equation of motion. However, the accuracy of analytical models is greatly reduced on more complicated devices with complex external loads, nonlinear trajectories or unique unlocking features. A 3D finite element model of CAD device with all critical features included can vastly improve the analytical ballistic predictions when it is used as a supplement to the ballistic code. During this project, LS-DYNA structural 3D model was used to predict the riser resisting load that was needed for the ballistic code. A Lagrangian model with eroding elements shown in Figure 2 was used for the blade, steel riser and the anvil. The riser material failure strain was fine tuned by matching the dent depth on the anvil with the actual test data. LS-DYNA model was also utilized to optimize the blade tip design for the most efficient cut. In parallel, the propellant type and the amount were determined by using CADPROG internal ballistics code. Initial test results showed a good match with LS-DYNA and CADPROG simulations. Final paper will present a detailed roadmap from initial ballistic modeling and LS-DYNA simulation to the performance testing. Blade shape optimization study will also be presented.

Oguz, Sirri↗

Adaptive optics and optical structures; Proceedings of the Meeting, European Congress on Optics, 3rd, The Hague, Netherlands, Mar. 12-14, 1990

The present conference on adaptive optics (AO) and optical structures addresses AO systems and controls, AO components, nonlinear optics applications to AO, astronomical applications of AO, large telescopes and optical alignment, as well as the wavefront control experiment for the use of AO in beam propagation. Specific references are made to applications of electromagnetic theory in optics, theoretical studies of system performance and design parameters, Hartmann-Shack wavefront sensing, the use of ray-based techniques in cophasing segmented mirrors, the use of a phase-conjugating mirror for real-time phase visualization, and the absolute instability of oppositely directed waves with respect to a high-reflectivity phase-conjugate mirror. Also addressed are automatic control systems, precision segmented reflectors, AO system creation, the VLT's 8.2-m primary mirrors, an optical 12-m telescope, alignment optimization via the Talbot effect, and a combination of interferometry and ray-tracing analysis.

Tyson, Robert K.↗

Modifying PASVART to solve singular nonlinear 2-point boundary problems

To study the buckling and post-buckling behavior of shells and various other structures, one must solve a nonlinear 2-point boundary problem. Since closed-form analytic solutions for such problems are virtually nonexistent, numerical approximations are inevitable. This makes the availability of accurate and reliable software indispensable. In a series of papers Lentini and Pereyra, expanding on the work of Keller, developed PASVART: an adaptive finite difference solver for nonlinear 2-point boundary problems. While the program does produce extremely accurate solutions with great efficiency, it is hindered by a major limitation. PASVART will only locate isolated solutions of the problem. In buckling problems, the solution set is not unique. It will contain singular or bifurcation points, where different branches of the solution set may intersect. Thus, PASVART is useless precisely when the problem becomes interesting. To resolve this deficiency we propose a modification of PASVART that will enable the user to perform a more complete bifurcation analysis. PASVART would be combined with the Thurston bifurcation solution: as adaptation of Newton's method that was motivated by the work of Koiter 3 are reinterpreted in terms of an iterative computational method by Thurston.

Fulton, James P.↗

Partial Wave Analysis of Coupled Photonic Structures

The very high quality factors sustained by microcavity optical resonators are relevant to applications in wavelength filtering, routing, switching, modulation, and multiplexing/demultiplexing. Increases in the density of photonic elements require that attention be paid to how electromagnetic (EM) coupling modifies their optical properties. This is especially true when cavity resonances are involved, in which case, their characteristics may be fundamentally altered. Understanding the optical properties of microcavities that are near or in contact with photonic elements---such as other microcavities, nanostructures, couplers, and substrates---can be expected to advance our understanding of the roles that these structures may play in VLSI photonics, biosensors and similar device technologies. Wc present results from recent theoretical studies of the effects of inter- and intracavity coupling on optical resonances in compound spherical particles. Concentrically stratified spheres and bispheres constituted from homogeneous and stratified spheres are subjects of this investigation. A new formulation is introduced for the absorption of light in an arbitrary layer of a multilayered sphere, which is based on multiple reflections of the spherical partial waves of the Lorenz-Mie solution for scattering by a sphere. Absorption efficiencies, which can be used to profile cavity resonances and to infer fluorescence yields or the onset of nonlinear optical processes in the microcavities, are presented. Splitting of resonances in these multisphere systems is paid particular attention, and consequences for photonic device development and possible performance enhancements through carefully designed architectures that exploit EM coupling are considered.

Fuller, Kirk A.↗

Aeroelastic stability of complete rotors with application to a teetering rotor in forward flight

The derivation of a set of nonlinear coupled flap-lag-torsion equations of motion for moderately large deflections of an elastic, two bladed teetering helicopter rotor in forward flight is concisely outlined. The following degrees of freedom are included in the mathematical model: rigid body flapping, rigid body lead lag, elastic bending in flap and lead-lag blade root torsion and shaft torsion. Quasi-steady aerodynamic loads are considered and the effects of reversed flow are included. The aeroelastic stability of the complete rotor is investigated using a linearized system of equations of motion. The equilibrium position about which the equations are linearized is obtained by considering the trim state of the helicopter, in true or simulated forward flight conditions. The sensitivity of the aeroelastic stability boundaries to interblade structural and mechanical coupling is illustrated by comparing the complete rotor stability boundaries with those obtained from a single blade analysis for a number of hover and forward flight cases.

Shamie, J.↗

On Certain Theoretical Developments Underlying the Hilbert-Huang Transform

One of the main traditional tools used in scientific and engineering data spectral analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as being linear and stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectral analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposed data, the HHT allows spectral analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real-value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a nearly orthogonal derived from the data (adaptive) basis. The IMFs can be further analyzed for spectrum content by using the classical Hilbert Transform. A new engineering spectral analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications pose additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs nearly orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the development of new HHT processing options, such as real-time and 2-D processing using Field Programmable Gate Array (FPGA) computational resources,

Kizhner, Semion↗

Cascading economic losses from port disruptions under capacity constrained multimodal freight networks

This study quantifies how throughput disruptions at major seaports cascade through capacity-constrained multimodal freight networks and interregional production systems. We couple an agent-based model (ABM) multimodal freight simulation that resolves rerouting, terminal queueing, and inventory drawdown under binding modal and facility capacities with a multiregional output loss input-output (MRIIM) model that propagates realized delivery shortfalls across regions and sectors. The framework is demonstrated for the Port of Los Angeles using Freight Analysis Framework flows and Bureau of Economic Analysis input-output accounts and is evaluated over a 52-week horizon under deterministic sector targeted shocks and stochastic disruption realizations with uncertain severity and duration. Results indicate nonlinear amplification: realized national losses concentrate in manufacturing and transportation/warehousing even when exogenous port shocks are dispersed, suggesting that congestion spillback and limited short-run substitution can dominate the initial shock allocation. We further evaluate a tabular reinforcement-learning (Q-learning) intervention layer that selects among a small set of implementable system level levers (truck-to-rail and truck-to-barge shift settings) without overriding shipper routing, finding that such interventions reduce total losses for moderate disruptions but yield diminishing returns once substitute modes approach capacity. By linking operational freight behavior to system wide impacts under uncertainty, the proposed ABM-MRIIM pipeline provides a reusable workflow for port disruption stress testing, identification of structurally critical sectors/corridors, and evaluation of resilience interventions under realistic capacity limits.

42 ENGINEERING↗

In Situ Observation of Ion Migration in a Ferroelectric Ionic Conductor Rb-KTP during Thermal Annealing

Ion exchange in Rb-doped KTiOPO 4 has facilitated significant advancements in ferroelectric domain engineering, yet understanding the underlying mechanisms remains in its infancy. We perform time-of-flight secondary ion mass spectrometry analysis on multiple periodically ion-exchanged and periodically poled Rb-doped KTiOPO 4 samples under different temperatures and annealing durations. The results are compared between annealing in air, which involved ex situ annealing before periodic poling, and vacuum annealing conducted in situ after periodic poling. The Rb + diffusion profile after periodic ion exchange forms a tooth-shaped pattern. We show that in situ annealing causes a surface pinning effect on the nonpolar face, limiting Rb + migration along the polar axis at the surface. Once the pinned layer is removed through milling, the underlying Rb + diffusion is distinctively different from the surface. Additionally, the rate of Rb + diffusion during in situ annealing is linear, while the periodic domain structures remain stable after annealing. These results contribute to understanding the ionic diffusion process in a ferroelectric ionic conductor and using ion exchange to tailor the linear and nonlinear properties of KTiOPO 4 .

36 MATERIALS SCIENCE↗

The Generation of Smooth High Speed Solar Wind From Plume-Interplume Mixing

Ulysses has shown that fast solar wind is extremely smooth, with a variance of less than 5%, in contrast to slow wind with a variance of approximately 30%. Now UVCS has produced the surprising result that the flow speed within coronal holes, the source of fast wind, is not at all smooth. Specifically, Giordano et al. (ApJ, v531, L79-L82, 2000) report that at 1.7 R(sub SUN) the interplume flow speed is typically more than twice the plume flow speed. Other less direct evidence supports this same result, with speeds from less than 300 to over 1000 km/s reported at approximately 5 R(sub SUN). This presents the paradox of how strongly differing plume and interplume flow speeds can exist near the Sun and be absent far from the Sun. The only answer is that plume and interplume material or momentum must be strongly mixed and that the mixing must occur mainly inside 0.3 AU to be consistent with Helios observations of smooth fast wind. Pressure balance structures (PBSs) and He abundance anomalies (Reisenfeld, et al., GRL, v26(13), 1805-1808, 1999) have been identified as interplanetary remnants of plumes, implying momentum mixing is the dominant coronal process. One possible source for plume/interplume momentum mixing is MHD Kelvin-Helmholtz (KH) shear instabilities occurring on the velocity shear interfaces. The velocity shear is a source of free energy and KH fluctuations could, through nonlinear cascade, provide the forcing required for the plasma oscillations (Cranmer, ApJ, v532, 1197-1208, 2000) reported to exist in coronal holes. The physical properties in coronal holes are now sufficiently well known that we can show plume/interplume shear interfaces become unstable to the KH instability at 5 - 10 R(sub SUN). The KH dispersion relation can be used to analyze marginal stability, the most unstable wavelengths, and linear growth rates. Numerical simulations can be used to verify results from the linear analysis and study the nonlinear development of KH modes. Here we will describe the marginal stability criterion for the KH instability, how and where this condition is met in the corona, and the possible character of the resultant fluctuations.

Suess, Steve↗

On the Hilbert-Huang Transform Theoretical Developments

One of the main heritage tools used in scientific and engineering data spectrum analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as linearity, of being stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectrum analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposition data, the HHT allows spectrum analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a near orthogonal adaptive basis, a basis that is derived from the data. The IMFs can be further analyzed for spectrum interpretation by the classical Hilbert Transform. A new engineering spectrum analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications post additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs near orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the developments of new HHT processing options, such as real-time and 2-D processing using Field Programmable Array (FPGA) computational resources, enhanced HHT synthesis, and broaden the scope of HHT applications for signal processing.

Kizhner, Semion↗