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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Aeroelastic Modeling of Elastically Shaped Aircraft Concept via Wing Shaping Control for Drag Reduction

Lightweight aircraft design has received a considerable attention in recent years as a means for improving cruise efficiency. Reducing aircraft weight results in lower lift requirements which directly translate into lower drag, hence reduced engine thrust requirements during cruise. The use of lightweight materials such as advanced composite materials has been adopted by airframe manufacturers in current and future aircraft. Modern lightweight materials can provide less structural rigidity while maintaining load-carrying capacity. As structural flexibility increases, aeroelastic interactions with aerodynamic forces and moments become an increasingly important consideration in aircraft design and aerodynamic performance. Furthermore, aeroelastic interactions with flight dynamics can result in issues with vehicle stability and control. Abstract This paper describes a recent aeroelastic modeling effort for an elastically shaped aircraft concept (ESAC). The aircraft model is based on the rigid-body generic transport model (GTM) originally developed at NASA Langley Research Center. The ESAC distinguishes itself from the GTM in that it is equipped with highly flexible wing structures as a weight reduction design feature. More significantly, the wings are outfitted with a novel control effector concept called variable camber continuous trailing edge (VCCTE) flap system for active control of wing aeroelastic deflections to optimize the local angle of attack of wing sections for improved aerodynamic efficiency through cruise drag reduction and lift enhancement during take-off and landing. The VCCTE flap is a multi-functional and aerodynamically efficient device capable of achieving high lift-to-drag ratios. The flap system is comprised of three chordwise segments that form the variable camber feature of the flap and multiple spanwise segments that form a piecewise continuous trailing edge. By configuring the flap camber and trailing edge shape, drag reduction could be achieved. Moreover, some parts of the flap system can be made to have a high frequency response for roll control, gust load alleviation, and aeroservoelastic (ASE) modal suppression control. Abstract The aeroelastic model of the ESAC is based on one-dimensional structural dynamic theory that captures the aeroelastic deformation of a wing structure in a combined motion that involves flapwise bending, chordwise bending, and torsion. The model includes the effect of aircraft propulsion due to wing flexibility which causes the propulsive forces and moments to couple with the wing elastic motion. Engine mass is also accounted in the model. A fuel management model is developed to describe the wing mass change due to fuel usage in the main tank and wing tanks during cruise. Abstract The model computes both static and dynamic responses of the wing structures. The static aeroelastic deflections are used to estimate the effect of wing flexibility on induced drag and the potential drag reduction by the VCCTE flap system. A flutter analysis is conducted to estimate the flutter speed boundary. Gust load alleviation via adaptive control has been recently investigated to address flexibility of aircraft structures. A multi-objective flight control approach is presented for drag reduction control. The approach is based on an optimal control framework using a multi-objective cost function. Future studies will demonstrate the potential benefits of the approach.

Aeroelasticity↗
DOE OSTI2021

Morphology and Dynamics of Bottlebrush Polymers

This thesis makes significant advances to the understanding of bottlebrush polymers. While bottlebrushes have received much attention due to the recent discovery of their unprecedented properties, including supersoftness, ultra-low viscosity, and hyperelasticity, this thesis is the first fundamental investigation at the molecular level that comprises structure and dynamics. Neutron scattering experiments, detailed within, reveal spherical or cylindrical shapes, instead of a random coil conformation. Another highlight is the analysis of the fast dynamics at the sub nm-length scale. The combination of three neutron spectrometers and the development of a new analysis technique enabled the calculation of the mean-square displacement over seven orders of magnitude in time scale. This unprecedented result can be applied to a broad class of samples, including polymers and other materials. The thesis is accessible to scientists from other fields, provides the reader with easily understandable guidelines for applying this analysis to other materials, and has the potential to make a significant impact on the analysis of neutron scattering data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Sensitivity Analysis of Numerical Modeling Input Parameters on Wind Turbine Loads in Deterministic Transient Load Cases

Aero-hydro-elastic-servo numerical models used to design and analyze wind turbines are based on thousands of variable input parameters that dictate the inflow, aerodynamic, structural, and control characteristics of the system as well as sea state, hydrodynamic, and mooring characteristics for fixed-bottom and floating offshore wind turbines. Each of these parameters has some level of uncertainty, which can significantly impact the predicted loads. Understanding the uncertainty in the inputs is critical to understanding the uncertainty in the outputs. This work demonstrates a screening technique to identify which parameters ultimate loads are most sensitive to so that more focus can be given to quantifying the possible range of those parameters. This technique has been demonstrated previously for different turbine and load case types and is extended here for a floating offshore wind turbine in design load cases with transient events both in the inflow and operations. Each load case features a deterministic gust, including variations in wind speed, direction, and shear. Load cases are considered with an operating turbine as well as with prescribed fault, startup, and shutdown procedures. The study found that key input parameters with a large impact on loads include the length of the gust, the magnitude of direction change and speed in the gust, the initial wind speed, and the shape of the gust profile.

17 WIND ENERGY↗
DOE OSTI2022

Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals

In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.

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

The method of lines in analyzing solids containing cracks

A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.

Gyekenyesi, John P.↗
DOE OSTI2022

Damage Tolerance Assessment of Laser Clad Repairs of Coarse Grain Ti-6Al-4V

Background: Laser cladding offers a promising approach to repairing aviation components, with the ability to restore the geometric shape and the static and fatigue strength. Objective: The damage tolerance of a coarse grain Ti-6Al-4 V material is assessed via three sets of samples: (i) baseline, (ii) laser clad repaired, and (iii) laser cladding followed by heat treatment. Method: The specimens were subjected to fatigue loading using a marker band technique, in which cracks initiated from pre-defined laser slots and the small crack growth rates were measured by quantitative fractography. Additional characterization was performed on these samples, including hardness measurements, microstructure analysis, and multimodal residual stress analyses. Results: The region of the specimen undergoing laser cladding exhibited higher hardness values and a compressive residual stress, based on the slitting method. The samples subjected to heat treatment after cladding displayed negligible differences in grain size and minimal residual stress relaxation. Conclusions: The results demonstrated significant variability, with data demonstrating that cracks in the laser clad region (exhibiting the largest compressive residual stress) grew faster than cracks in the substrate region of the sample or in the baseline cases. Finally, these counterintuitive results are rationalized with additional analyses from synchrotron-based energy dispersive X-ray diffraction and high energy X-ray diffraction, demonstrating significant spatial variability in the residual elastic strains in the millimeter-sized prior β grains within the material.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Fully Coupled Aero-Thermochemical-Elastic Simulations of an Eroding Graphite Nozzle

A multiphysics simulation capability has been developed that incorporates mutual interactions between aerodynamics, structural response from aero/thermal loading, ablation/pyrolysis, heating, and surface-to-surface radiation to perform high-fidelity, fully coupled aerothermoelastic ablation simulations, which to date had been unattainable. The multiphysics framework couples CHAR (a 3-D implicit charring ablator solver), Loci/CHEM (a computational fluid dynamics solver for high-speed chemically reacting flows), and Abaqus (a nonlinear structural dynamics solver) to create a fully coupled aerothermoelastic charring ablative solver. The solvers are tightly coupled in a fully integrated fashion to resolve the effects of the ablation pyrolysis and charring process and chemistry products upon the flow field, the changes in surface geometry due to recession upon the flow field, and thermal-structural analysis of the body from the induced aerodynamic heating from the flow field. The multiphysics framework was successfully demonstrated on a solid rocket motor graphite nozzle erosion application. Comparisons were made with available experimental data that measured the throat erosion during the motor firing. The erosion data is well characterized, as the test rig was equipped with a windowed nozzle section for real-time X-ray radiography diagnostics of the instantaneous throat variations for deducing the instantaneous erosion rates. The nozzle initially undergoes a nozzle contraction due to thermal expansion before ablation effects are able to widen the throat. A series of parameters studies were conducted using the coupled simulation capability to determine the sensitivity of the nozzle erosion to different parameters. The parameter studies included the shape of the nozzle throat (flat versus rounded), the material properties, the effect of the choice of turbulence model, and the inclusion or exclusion of the mechanical thermal expansion. Overall, the predicted results match the experiment very well, and the predictions were able to bound the data within acceptable limits.

Blades, E. L.↗
NASA NTRS2000

Mechanical properties of a biodegradable bone regeneration scaffold

Poly (Propylene Fumarate) (PPF), a novel, bulk erosion, biodegradable polymer, has been shown to have osteoconductive effects in vivo when used as a bone regeneration scaffold (Peter, S. J., Suggs, L. J., Yaszemski, M. J., Engel, P. S., and Mikos, A. J., 1999, J. Biomater. Sci. Polym. Ed., 10, pp. 363-373). The material properties of the polymer allow it to be injected into irregularly shaped voids in vivo and provide mechanical stability as well as function as a bone regeneration scaffold. We fabricated a series of biomaterial composites, comprised of varying quantities of PPF, NaCl and beta-tricalcium phosphate (beta-TCP), into the shape of right circular cylinders and tested the mechanical properties in four-point bending and compression. The mean modulus of elasticity in compression (Ec) was 1204.2 MPa (SD 32.2) and the mean modulus of elasticity in bending (Eb) was 1274.7 MPa (SD 125.7). All of the moduli were on the order of magnitude of trabecular bone. Changing the level of NaCl from 20 to 40 percent, by mass, did not decrease Ec and Eb significantly, but did decrease bending and compressive strength significantly. Increasing the beta-TCP from 0.25 g/g PPF to 0.5 g/g PPF increased all of the measured mechanical properties of PPF/NVP composites. These results indicate that this biodegradable polymer composite is an attractive candidate for use as a replacement scaffold for trabecular bone.

NASA Discipline Cell Biology↗
NASA NTRS1997

Mechanical continuity and reversible chromosome disassembly within intact genomes removed from living cells

Chromatin is thought to be structurally discontinuous because it is packaged into morphologically distinct chromosomes that appear physically isolated from one another in metaphase preparations used for cytogenetic studies. However, analysis of chromosome positioning and movement suggest that different chromosomes often behave as if they were physically connected in interphase as well as mitosis. To address this paradox directly, we used a microsurgical technique to physically remove nucleoplasm or chromosomes from living cells under isotonic conditions. Using this approach, we found that pulling a single nucleolus or chromosome out from interphase or mitotic cells resulted in sequential removal of the remaining nucleoli and chromosomes, interconnected by a continuous elastic thread. Enzymatic treatments of interphase nucleoplasm and chromosome chains held under tension revealed that mechanical continuity within the chromatin was mediated by elements sensitive to DNase or micrococcal nuclease, but not RNases, formamide at high temperature, or proteases. In contrast, mechanical coupling between mitotic chromosomes and the surrounding cytoplasm appeared to be mediated by gelsolin-sensitive microfilaments. Furthermore, when ion concentrations were raised and lowered, both the chromosomes and the interconnecting strands underwent multiple rounds of decondensation and recondensation. As a result of these dynamic structural alterations, the mitotic chains also became sensitive to disruption by restriction enzymes. Ion-induced chromosome decondensation could be blocked by treatment with DNA binding dyes, agents that reduce protein disulfide linkages within nuclear matrix, or an antibody directed against histones. Fully decondensed chromatin strands also could be induced to recondense into chromosomes with pre-existing size, shape, number, and position by adding anti-histone antibodies. Conversely, removal of histones by proteolysis or heparin treatment produced chromosome decondensation which could be reversed by addition of histone H1, but not histones H2b or H3. These data suggest that DNA, its associated protein scaffolds, and surrounding cytoskeletal networks function as a structurally-unified system. Mechanical coupling within the nucleoplasm may coordinate dynamic alterations in chromatin structure, guide chromosome movement, and ensure fidelity of mitosis.

NASA Discipline Cell Biology↗
NASA NTRS1988

Stiffnesses of a solid-rocket motor from an ambient vibration survey

Experience with many spacecraft configurations boosted by a variety of launch vehicles indicates that the maximum loads experienced throughout most of the structure are inertial in origin. These loads arise from the dynamic elastic response of the flight vehicle to the transient disturbances of launch and flight, and are highly dependent on the dynamic characteristics of both the spacecraft and the launch vehicle. It has proved to be most advantageous, in the analysis of this critical dependency of loads upon vehicle dynamic properties, to establish a mathematical model in terms of normal mode characteristics. In this way, the vibration behavior of an elastomechanical structure (or substructure) can be described by means of the so-called modal or natural degrees of freedom. The conduct of a mode survey test and the use of a suitably test-verified model in loads analyses is essential to the flight worthiness certification process of space systems. The desirability of such tests is confirmed by the fact that, almost invariably, significant deficiencies in the analytical models are revealed by the results. Therefore, this experimental program was undertaken to determine those properties of a solid-propellant rocket motor (SRM) which are required to characterize a dynamic model. Random ambient-excited accelerations were measured at a series of stations along the motor for the purpose of identifying the motor beam-like stiffnesses in bending, shear, and torsion. From a system identification point of view, it is significant that stiffness properties of a subsystem (the motor) are determined from modes of the full system (motor/stand configuration) using mode shape data of the subsystem only. This contrasts with traditional system identification approaches which rely upon complete system mode shapes.

Rubin, S.↗
DOE OSTI2020

Regression models using shapes of functions as predictors

Functional variables are often used as predictors in regression problems. A commonly used parametric approach, called scalar-on-function regression, uses the $\mathbb L^2$ inner product to map functional predictors into scalar responses. This method can perform poorly when predictor functions contain undesired phase variability, causing phases to have disproportionately large influence on the response variable. One past solution has been to perform phase–amplitude separation (as a pre-processing step) and then use only the amplitudes in the regression model. In this paper, we propose a more integrated approach, termed elastic functional regression model (EFRM), where phase-separation is performed inside the regression model, rather than as a pre-processing step. This approach generalizes the notion of phase in functional data, and is based on the norm-preserving time warping of predictors. Due to its invariance properties, this representation provides robustness to predictor phase variability and results in improved predictions of the response variable over traditional models. We demonstrate this framework using a number of datasets involving gait signals, NMR data, and stock market prices.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

A multiscale model was developed for predicting the thermoelastic behavior of semi-crystalline thermoplastic materials for composite aerospace applications. At the highest scale containing the semi-crystalline spherulite in an amorphous matrix, the generalized method of cells, or high fidelity method of cells, was used to perform the homogenization calculations to obtain the effective properties. Models were developed assuming a cubic, or spherical shape, for the spherulite to understand if the morphology of the spherulite affects the effective thermoelastic properties. The generalized method of cells was used to model at the repeating unit cells at the subscales of the microstructure including the lamellae stacks and granular crystal blocks. The scales are integrated using the multiscale micromechanics method in the NASA Multiscale Analysis Tool. Data from molecular dynamics simulations were used as inputs for the amorphous and crystalline constituents. Convergence studies were performed to determine the best level of discretization for the repeating unit cell at the highest scale. Effective Young’s modulus, shear modulus, Poisson’s ratio, coefficient of thermal expansion, and thermal conductivity were predicted for polyether ether ketone and polyether ketone ketone, and very good agreement between the model utilizing the cubic spherulite and the experimental data, where available, was observed for polyether ketone ketone. Normalization of the data for the bulk polyether ketone ketone, against amorphous data, improved the predictions as compared to experimental data. Overall, the high fidelity method of cells predicted a stiffer response then the generalized method of cells as the crystallinity was increased. The shape of the spherulite had a minimal effect on the predicted bulk properties of the polymers.

thermoplastics↗
NASA NTRSDate not supplied

Nonlinear Buckling Analysis of Tow-Steered Composite Cylinders with Cutouts

The buckling and postbuckling behavior of two composite tow-steered shells with cutouts of different sizes is assessed using nonlinear finite element (FE) analysis and compared to experimental measurements. The cylindrical shells are manufactured using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from ±10 degrees on the crown and keel to ±45 degrees on the sides. One shell features thickness variations due to tow overlaps that result from application of all 24 tows during each pass of the fiber placement system. The second shell uses the system’s tow drop/add capability to achieve a more uniform wall thickness without overlaps. Unreinforced cutouts of two different sizes—the first smaller cutout representing a passenger door on a commercial aircraft and the second larger cutout a cargo door—were machined into each of the two cylinders resulting in a total of four test cases. These cylinders were tested in axial compression and buckled elastically in previous work and are now analyzed using nonlinear FE models to compare bifurcation buckling loads as well as the load-displacement response in the prebuckling and postbuckling regimes. For all four shells analyzed, the prebuckling stiffness, buckling load, and deformation mode sequence throughout the loading-unloading cycle is accurately reproduced by the models. In particular, the shells first buckle locally around the cutouts in a stable (super-critical) manner with only a slight decrease in axial stiffness, which occurs due to the favorable load redistribution facilitated by tow steering. The shells then buckle globally in an unstable (sub-critical) manner with diamond-shaped buckles forming to the left and right of the cutouts. The buckling load of all shells with cutouts is at least 82% of the buckling load of the pristine shells without cutouts. Overall, the ability to sustain local buckling phenomena, and the relatively small reductions in global buckling load compared to pristine shells without cutouts, demonstrates the great potential of using tow steering to mitigate the adverse effects of cutouts in axially-compressed shell structures.

composites↗
NASA NTRSDate not supplied

Nonlinear Buckling Analysis of Tow-Steered Composite Cylinders with Cutouts

The buckling and postbuckling behavior of two composite tow-steered shells with cutouts of different sizes is assessed using nonlinear finite element (FE) analysis and compared to experimental measurements. The cylindrical shells are manufactured using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from 10 degrees on the crown and keel to 45 degrees on the sides. One shell features thickness variations due to tow overlaps that result from application of all 24 tows during each pass of the fiber placement system. The second shell uses the system’s tow drop/add capability to achieve a more uniform wall thickness without overlaps. Unreinforced cutouts of two different sizes—the first smaller cutout representing a passenger door on a commercial aircraft and the second larger cutout a cargo door—were machined into each of the two cylinders resulting in a total of four test cases. These cylinders were tested in axial compression and buckled elastically in previous work and are now analyzed using nonlinear FE models to compare bifurcation buckling loads as well as the load-displacement response in the prebuckling and postbuckling regimes. For all four shells analyzed, the prebuckling stiffness, buckling load, and deformation mode sequence throughout the loading-unloading cycle is accurately reproduced by the models. In particular, the shells first buckle locally around the cutouts in a stable (super-critical) manner with only a slight decrease in axial stiffness, which occurs due to the favorable load redistribution facilitated by tow steering. The shells then buckle globally in an unstable (sub-critical) manner with diamond-shaped buckles forming to the left and right of the cutouts. The buckling load of all shells with cutouts is at least 82% of the buckling load of the pristine shells without cutouts. Overall, the ability to sustain local buckling phenomena, and the relatively small reductions in global buckling load compared to pristine shells without cutouts, demonstrates the great potential of using tow steering to mitigate the adverse effects of cutouts in axially-compressed shell structures.

Buckling↗
DOE OSTI2021

Seismic waves in medium with poroelastic/elastic interfaces: a two-dimensional P-SV finite-difference modelling

SUMMARY We present a new methodology of the finite-difference (FD) modelling of seismic wave propagation in a strongly heterogeneous medium composed of poroelastic (P) and (strictly) elastic (E) parts. The medium can include P/P, P/E and E/E material interfaces of arbitrary shapes. The poroelastic part can be with (i) zero resistive friction, (ii) non-zero constant resistive friction or (iii) JKD model of the frequency-dependent permeability and resistive friction. Our FD scheme is capable of subcell resolution: a material interface can have an arbitrary position in the spatial grid. The scheme keeps computational efficiency of the scheme for a smoothly and weakly heterogeneous medium (medium without material interfaces). Numerical tests against independent analytical, semi-analytical and spectral-element methods prove the efficiency and accuracy of our FD modelling. In numerical examples, we indicate effect of the P/E interfaces for the poroelastic medium with a constant resistive friction and medium with the JKD model of the frequency-dependent permeability and resistive friction. We address the 2-D P-SV problem. The approach can be readily extended to the 3-D problem.

58 GEOSCIENCES↗
NASA NTRS2009

Subspace Iteration Method for Complex Eigenvalue Problems with Nonsymmetric Matrices in Aeroelastic System

Modern airplane design is a multidisciplinary task which combines several disciplines such as structures, aerodynamics, flight controls, and sometimes heat transfer. Historically, analytical and experimental investigations concerning the interaction of the elastic airframe with aerodynamic and in retia loads have been conducted during the design phase to determine the existence of aeroelastic instabilities, so called flutter .With the advent and increased usage of flight control systems, there is also a likelihood of instabilities caused by the interaction of the flight control system and the aeroelastic response of the airplane, known as aeroservoelastic instabilities. An in -house code MPASES (Ref. 1), modified from PASES (Ref. 2), is a general purpose digital computer program for the analysis of the closed-loop stability problem. This program used subroutines given in the International Mathematical and Statistical Library (IMSL) (Ref. 3) to compute all of the real and/or complex conjugate pairs of eigenvalues of the Hessenberg matrix. For high fidelity configuration, these aeroelastic system matrices are large and compute all eigenvalues will be time consuming. A subspace iteration method (Ref. 4) for complex eigenvalues problems with nonsymmetric matrices has been formulated and incorporated into the modified program for aeroservoelastic stability (MPASES code). Subspace iteration method only solve for the lowest p eigenvalues and corresponding eigenvectors for aeroelastic and aeroservoelastic analysis. In general, the selection of p is ranging from 10 for wing flutter analysis to 50 for an entire aircraft flutter analysis. The application of this newly incorporated code is an experiment known as the Aerostructures Test Wing (ATW) which was designed by the National Aeronautic and Space Administration (NASA) Dryden Flight Research Center, Edwards, California to research aeroelastic instabilities. Specifically, this experiment was used to study an instability known as flutter. ATW was a small-scale airplane wing comprised of an airfoil and wing tip boom. This wing was formulated based on a NACA-65A004 airfoil shape with a 3.28 aspect ratio. The wing had a span of 18 inch with root chord length of 13.2 inch and tip chord length of 8.7 inch. The total area of this wing was 197 square inch. The wing tip boom was a 1 inch diameter hollow tube of length 21.5 inch. The total weight of the wing was 2.66 lbs.

Pak, Chan-gi↗