Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational structural dynamics”

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,531 records · Page 85

On The Computational Capabilities of Physical Systems: Relationship With Conventional Computer Science - Part 2

In the first of this pair of papers, it was proven that there cannot be a physical computer to which one can properly pose any and all computational tasks concerning the physical universe. It was then further proven that no physical computer C can correctly carry out all computational tasks that can be posed to C. As a particular example, this result means that no physical computer that can, for any physical system external to that computer, take the specification of that external system's state as input and then correctly predict its future state before that future state actually occurs; one cannot build a physical computer that can be assured of correctly "processing information faster than the universe does". These results do not rely on systems that are infinite, and/or non-classical, and/or obey chaotic dynamics. They also hold even if one uses an infinitely fast, infinitely dense computer, with computational powers greater than that of a Turing Machine. This generality is a direct consequence of the fact that a novel definition of computation - "physical computation" - is needed to address the issues considered in these papers, which concern real physical computers. While this novel definition does not fit into the traditional Chomsky hierarchy, the mathematical structure and impossibility results associated with it have parallels in the mathematics of the Chomsky hierarchy. This second paper of the pair presents a preliminary exploration of some of this mathematical structure. Analogues of Chomskian results concerning universal Turing Machines and the Halting theorem are derived, as are results concerning the (im)possibility of certain kinds of error-correcting codes. In addition, an analogue of algorithmic information complexity, "prediction complexity", is elaborated. A task-independent bound is derived on how much the prediction complexity of a computational task can differ for two different reference universal physical computers used to solve that task, a bound similar to the "encoding" bound governing how much the algorithm information complexity of a Turing machine calculation can differ for two reference universal Turing machines. Finally, it is proven that either the Hamiltonian of our universe proscribes a certain type of computation, or prediction complexity is unique (unlike algorithmic information complexity), in that there is one and only version of it that can be applicable throughout our universe.

Wolpert, David H.↗

A Time Integration Algorithm Based on the State Transition Matrix for Structures with Time Varying and Nonlinear Properties

A variable order method of integrating the structural dynamics equations that is based on the state transition matrix has been developed. The method has been evaluated for linear time variant and nonlinear systems of equations. When the time variation of the system can be modeled exactly by a polynomial it produces nearly exact solutions for a wide range of time step sizes. Solutions of a model nonlinear dynamic response exhibiting chaotic behavior have been computed. Accuracy of the method has been demonstrated by comparison with solutions obtained by established methods.

Bartels, Robert E.↗

Frequencies and Flutter Speed Estimation for Damaged Aircraft Wing Using Scaled Equivalent Plate Analysis

Equivalent plate analysis is often used to replace the computationally expensive finite element analysis in initial design stages or in conceptual design of aircraft wing structures. The equivalent plate model can also be used to design a wind tunnel model to match the stiffness characteristics of the wing box of a full-scale aircraft wing model while satisfying strength-based requirements An equivalent plate analysis technique is presented to predict the static and dynamic response of an aircraft wing with or without damage. First, a geometric scale factor and a dynamic pressure scale factor are defined to relate the stiffness, load and deformation of the equivalent plate to the aircraft wing. A procedure using an optimization technique is presented to create scaled equivalent plate models from the full scale aircraft wing using geometric and dynamic pressure scale factors. The scaled models are constructed by matching the stiffness of the scaled equivalent plate with the scaled aircraft wing stiffness. It is demonstrated that the scaled equivalent plate model can be used to predict the deformation of the aircraft wing accurately. Once the full equivalent plate geometry is obtained, any other scaled equivalent plate geometry can be obtained using the geometric scale factor. Next, an average frequency scale factor is defined as the average ratio of the frequencies of the aircraft wing to the frequencies of the full-scaled equivalent plate. The average frequency scale factor combined with the geometric scale factor is used to predict the frequency response of the aircraft wing from the scaled equivalent plate analysis. A procedure is outlined to estimate the frequency response and the flutter speed of an aircraft wing from the equivalent plate analysis using the frequency scale factor and geometric scale factor. The equivalent plate analysis is demonstrated using an aircraft wing without damage and another with damage. Both of the problems show that the scaled equivalent plate analysis can be successfully used to predict the frequencies and flutter speed of a typical aircraft wing.

Krishnamurthy, Thiagarajan↗

Nonlinear modeling and estimation of slew induced structural deformations

A model of the nonlinear dynamics of a deformable maneuvering multibody system is described, whereby elastic deformation are modeled by restoring forces and dissipative forces at point mass appendages. This model is brought into bilinear form. Estimation of deformations occasioned by rapid slewing maneuvers is carried out by a filter based on a globally equivalent linear model of the bilinear dynamics, and is shown to be an improvement over the extended Kalman filter. To further alleviate the computational burden, the estimated deformation state is propagated between observations by a low dimensional operator spline interpolator of bilinear system Volterra series, which is easily implemented.

Dwyer, T. A. W., III↗

Dynamic upper atmospheric force model on stabilized vehicles for a high-precision trajectory computer program

The upper atmosphere model draws heavily on the behavior of the earth's upper atmosphere which exhibits cyclic as well as irregular variations in density profile, temperature, pressure, and composition in unison with solar activities as deduced from the more recent land-based and satellite observations. The lift and drag model is designed specifically for inertially stabilized vehicles of the Mariner class, with possible extension to gravity gradient stabilized vehicles of the GEOS class. The model considers operation in the free molecular flow regimes with large Knudsen numbers. The vehicle is considered a composite structure with basic components having well-defined shapes, each with its own surface characteristics in terms of temperature, reflectivity, and accommodation of free stream molecules. The model takes into account both the calculation of precise aerodynamic force coefficients in terms of expansion of modified Bessel functions in speed ratios and angle of attack, and approximate force coefficients when the speed ratios approach infinity. Other considerations include specular and diffused reflectivity, shielding, and shadow effects.

Khatib, A. R.↗

Implicit flux-split Euler schemes for unsteady aerodynamic analysis involving unstructured dynamic meshes

Improved algorithms for the solution of the time-dependent Euler equations are presented for unsteady aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux-split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves an implicit time-integration shceme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. For example, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady and unsteady flow results are presented for the NACA 0012 airfoil to demonstrate applications of the new Euler solvers. The unsteady results were obtained for the airfoil pitching harmonically about the quarter chord. The resulting instantaneous pressure distributions and lift and moment coefficients during a cycle of motion compare well with experimental data. The paper presents a description of the Euler solvers along with results and comparisons which assess the capability.

Batina, John T.↗

Dynamic Analyses Including Joints Of Truss Structures

Method for mathematically modeling joints to assess influences of joints on dynamic response of truss structures developed in study. Only structures with low-frequency oscillations considered; only Coulomb friction and viscous damping included in analysis. Focus of effort to obtain finite-element mathematical models of joints exhibiting load-vs.-deflection behavior similar to measured load-vs.-deflection behavior of real joints. Experiments performed to determine stiffness and damping nonlinearities typical of joint hardware. Algorithm for computing coefficients of analytical joint models based on test data developed to enable study of linear and nonlinear effects of joints on global structural response. Besides intended application to large space structures, applications in nonaerospace community include ground-based antennas and earthquake-resistant steel-framed buildings.

Belvin, W. Keith↗

Analysis of Arc Welding Process in Space

This work is motivated by NASA plans to conduct welding experiments on the ISS. It is expected that deployment of welding and additive manufacturing technologies in the space environment has the potential to revolutionise how orbiting platforms are designed, manufactured, and assembled. However, the structure, composition and quality of a weld is extremely dependent on the environment and can be difficult to control in space. Shielding gases would also be tough to manage as gases behave differently in zero gravity and airless environments. Additional points of concern are related to the spatter and sparks dynamics in space. Therefore, there is a need for a more basic understanding of welding processes by computational modelling. To provide such an insight we developed state of art models of the ARC torch, droplet detachment and transfer, and the melt pool build up using magnetohydrodynamics approximation and level set method for two-phase liquid metal/-gas flow modelling. Two finite element models were built using 2D axisymmetric geometry in COMSOL Multiphysics®: (i) stationary model of the ARC torch and (ii) dynamical model of droplet detachment and transfer. Both models demonstrate reasonable agreement with earlier experimental observations and high sensitivity to the temperature dependence of the thermophysical parameters on the system materials. The models were used to provide physical insight into ARC welding in various environments and geometries. It was shown, in particular, that there is a significant probability of gas bubbles trapping in the meltpool and a possibility of unbounded wondering of sparks in the welding chamber in zero gravity. In addition, we estimated metal evaporation rate that may be hazardous in the confined environment of the ISS. These issues raise concern of quality of the weld and safety for ISS applications.

ARC welding↗

Development of a Probabilistic Dynamic Synthesis Method for the Analysis of Nondeterministic Structures

Accounting for the statistical geometric and material variability of structures in analysis has been a topic of considerable research for the last 30 years. The determination of quantifiable measures of statistical probability of a desired response variable, such as natural frequency, maximum displacement, or stress, to replace experience-based "safety factors" has been a primary goal of these studies. There are, however, several problems associated with their satisfactory application to realistic structures, such as bladed disks in turbomachinery. These include the accurate definition of the input random variables (rv's), the large size of the finite element models frequently used to simulate these structures, which makes even a single deterministic analysis expensive, and accurate generation of the cumulative distribution function (CDF) necessary to obtain the probability of the desired response variables. The research presented here applies a methodology called probabilistic dynamic synthesis (PDS) to solve these problems. The PDS method uses dynamic characteristics of substructures measured from modal test as the input rv's, rather than "primitive" rv's such as material or geometric uncertainties. These dynamic characteristics, which are the free-free eigenvalues, eigenvectors, and residual flexibility (RF), are readily measured and for many substructures, a reasonable sample set of these measurements can be obtained. The statistics for these rv's accurately account for the entire random character of the substructure. Using the RF method of component mode synthesis, these dynamic characteristics are used to generate reduced-size sample models of the substructures, which are then coupled to form system models. These sample models are used to obtain the CDF of the response variable by either applying Monte Carlo simulation or by generating data points for use in the response surface reliability method, which can perform the probabilistic analysis with an order of magnitude less computational effort. Both free- and forced-response analyses have been performed, and the results indicate that, while there is considerable room for improvement, the method produces usable and more representative solutions for the design of realistic structures with a substantial savings in computer time.

Brown, A. M.↗

Preliminary control/structure interaction study of coupled Space Station Freedom/Assembly Work Platform/orbiter

Results are presented from a preliminary control/structure interaction study of the Space Station, the Assembly Work Platform, and the STS orbiter dynamics coupled with the orbiter and station control systems. The first three Space Station assembly flight configurations and their finite element representations are illustrated. These configurations are compared in terms of control authority in each axis and propellant usage. The control systems design parameters during assembly are computed. Although the rigid body response was acceptable with the orbiter Primary Reaction Control System, the flexible body response showed large structural deflections and loads. It was found that severe control/structure interaction occurred if the stiffness of the Assembly Work Platform was equal to that of the station truss. Also, the response of the orbiter Vernier Reaction Control System to small changes in inertia properties is examined.

Singh, Sudeep K.↗

An Experimental Investigation on the Static Equilibria and Dynamics of Liquid Bridges

A liquid bridge is a volume of liquid held between two or more solid supports. In the case of small disk supports with a sharp edge, the contact line between the bridge and the support disk will be anchored along the edge of the disk. For these cases the solid presents a geometrical singularity and the contact angle is indeterminate within a given range. This dissertation presents research conducted on liquid bridges with anchored contact lines. The three major topics covered are: determining the role of support geometry on static equilibria, liquid bridge dynamical behavior, and forces exerted by a liquid bridge on a support structure. The work was primarily experimental and conducted in a "Plateau tank" that allowed for the simulation of equivalent low-gravity conditions. The main thrust of the experimental work involved the use of a high resolution optical measurement system for imaging the dynamic zone shape, measurement of the static and dynamic contact angles and non-invasive analysis of excited surface modes. The liquid bridge was manipulated by computer controlled linear actuators which allowed precise control over the physical characteristics of the bridge. Experiments have been carried out to locate a bifurcation point along the maximum volume axisymmetric stability margin. Below the critical slenderness the bifurcation from an axisymmetric to a stable nonaxisymmetric configuration is supercritical. However, above this critical slenderness, the bifurcation is subcritical. A series of experiments analyzed the effect on axisymmetric bridge stability by using support disks of different radii, The shape behavior as transition points were approached, as well as the limiting case of a vanishing support radius was investigated. Experiments were performed to determine the resonant frequencies of axisymmetric bridges subject to lateral vibrations. Anomolous results led to a series of experiments to characterize nonlinearities present in the dynamic bridge shape. Finally, an attempt was made to experimentally measure the force exerted by the bridge on the lower support disk. This was done through use of a force balance apparatus. Particular attention was paid to the behavior of the bridge as the minimum volume stability limit was approached.

Resnick, Andrew Howard↗

Health Monitoring of Composite Material Structures Using a Vibrometry Technique

Non-destructive evaluation (NDE) methods for quantifying and locating damage are essential for inspecting structures to ensure safety and reliability. Transmittance function monitoring is a potentially new NDE technique being tested as a tool to detect, quantify, and locate damage on flexible structures. The technique has a large spatial range that is practical for detecting damage on large composite material structures such as a reusable launch vehicle. The Transmittance Function (TF) theory is based on structural dynamics principles that define how vibration at one point in a structure is related to a force at another point. This relationship is called the Frequency Response Function (FRF). A Transmittance Function (TF) is derived as the ratio of FRFs, and can detect damage because the FRFs change due to damage. If one excitation is used for the testing, the force does not need to be measured to compute the TF. In the damage detection procedure, the structure is subjected to wide-band vibration and TFs are computed between different accelerometers to detect changes in the structure, presumably due to damage. In the first year of the project the TF method was tested on a bolted panel, a curved panel, and beams, all made of fiberglass. It was shown that damage could be detected using low frequency vibration, 250 to 1,250 Hz. The technique is sensitive to damage, but it requires storage of historical or pre-damage TFs for the healthy structure. This would become a large data storage requirement for large structures. Thus one objective for the second year of the project was to eliminate the need to store historical data. The second year report gives details of how storage of historical data was eliminated. Further results of testing panel structures are also given.

Schulz, Mark J.↗

Slot Nozzle Effects for Reduced Sonic Boom on a Generic Supersonic Wing Section

NASA has conducted research programs to reduce or eliminate the operational restrictions of supersonic aircraft over populated areas. Restrictions are due to the disturbance from the sonic boom, caused by the coalescence of shock waves formed off the aircraft. Results from two-dimensional computational fluid dynamic (CFD) analyses (performed on a baseline Mach 2.0 nozzle in a simulated Mach 2.2 flow) indicate that over-expanded and under-expanded operation of the nozzle has an effect on the N-wave boom signature. Analyses demonstrate the feasibility of reducing the magnitude of the sonic boom N-wave by controlling the nozzle plume interaction with the nozzle boat tail shock structure. This work was extended to study the impact of integrating a high aspect ratio exhaust nozzle or long slot nozzle on the trailing edge of a supersonic wing. The nozzle is operated in a highly under-expanded condition, creating a large exhaust plume and a shock at the trailing edge of the wing. This shock interacts with and suppresses the expansion wave caused by the wing, a major contributor to the sonic boom signature. The goal was to reduce the near field pressures caused by the expansion using a slot nozzle located at the wing trailing edge. Results from CFD analysis on a simulated wing cross-section and a slot nozzle indicate potential reductions in sonic boom signature compared to a baseline wing with no propulsion or trailing edge exhaust. Future studies could investigate if this effect could be useful on a supersonic aircraft for main propulsion, auxiliary propulsion, or flow control.

Castner, Raymond S.↗

Slot Nozzle Effects for Reduced Sonic Boom on a Generic Supersonic Wing Section

NASA has conducted research programs to reduce or eliminate the operational restrictions of supersonic aircraft over populated areas. Restrictions are due to the disturbance from the sonic boom, caused by the coalescence of shock waves formed off the aircraft. Results from two-dimensional computational fluid dynamic (CFD) analyses (performed on a baseline Mach 2.0 nozzle in a simulated Mach 2.2 flow) indicate that over-expanded and under-expanded operation of the nozzle has an effect on the N-wave boom signature. Analyses demonstrate the feasibility of reducing the magnitude of the sonic boom N-wave by controlling the nozzle plume interaction with the nozzle boat tail shock structure. This work was extended to study the impact of integrating a high aspect ratio exhaust nozzle or long slot nozzle on the trailing edge of a supersonic wing. The nozzle is operated in a highly under-expanded condition, creating a large exhaust plume and a shock at the trailing edge of the wing. This shock interacts with and suppresses the expansion wave caused by the wing, a major contributor to the sonic boom signature. The goal was to reduce the near field pressures caused by the expansion using a slot nozzle located at the wing trailing edge. Results from CFD analysis on a simulated wing cross-section and a slot nozzle indicate potential reductions in sonic boom signature compared to a baseline wing with no propulsion or trailing edge exhaust. Future studies could investigate if this effect could be useful on a supersonic aircraft for main propulsion, auxiliary propulsion, or flow control.

Caster, Raymond S.↗

Mechanically Tunable Radiative Cooling for Adaptive Thermal Control

Passive radiative cooling is currently thefrontier technology in renewable-energy research. Interms of extraterrestrial applications, radiative cooling is a critical component to the thermalmanagement system of a spacecraft, where the extreme environment of space can cause largetemperature variations that can break and damage equipment. For terrestrial applications,nocturnal or daytime radiative cooling is expected to lead to cost-effective passive heatmanagement without the need of inefficient and costly artificial refrigeration technologies.However, most currently available radiative cooling systems cannot be changed dynamically andradiate a constant static amount of thermal power. Dynamically tunable adaptive radiativecooling systems will be a critical development to prolong the lifetime of spacecraft or improvethe efficiency of terrestrial cooling systems. Here we propose stretchable radiative coolingdesigns that can be substantially tuned by using the simple physical mechanism of mechanical strain. When their structure is stretched, the radiated power is significantly reduced. We developa modeling method that can simulate mechanical stretching combined with electromagnetic response to compute the tunable thermal emission of these new adaptive radiative coolingsystems. The presented photonically engineered structures can be used as coatings to achieveefficient adaptive thermal control of various objects in a cost-effective and environmentallyfriendly way. The proposed designs are much simpler to be realized than others found in theliterature and the best design achieves a high thermal emission power with a tunable range on theorder of 132 W/m2.

Andrew Butler↗

Computed and Experimental Flutter/LCO Onset for the Boeing Truss-Braced Wing Wind-Tunnel Model

This paper presents high fidelity Navier-Stokes simulations of the Boeing Subsonic Ultra Green Aircraft Research truss-braced wing wind-tunnel model and compares the results to linear MSC. Nastran flutter analysis and preliminary data from a recent wind-tunnel test of that model at the NASA Langley Research Center Transonic Dynamics Tunnel. The simulated conditions under consideration are zero angle of attack, so that structural nonlinearity can be neglected. It is found that, for Mach number greater than 0.78, the linear flutter analysis predicts flutter onset dynamic pressure below the wind-tunnel test and that predicted by the Navier-Stokes analysis. Furthermore, the wind-tunnel test revealed that the majority of the high structural dynamics cases were wing limit cycle oscillation (LCO) rather than flutter. Most Navier-Stokes simulated cases were also LCO rather than hard flutter. There is dip in the wind-tunnel test flutter/LCO onset in the Mach 0.76-0.80 range. Conditions tested above that Mach number exhibited no aeroelastic instability at the dynamic pressures reached in the tunnel. The linear flutter analyses do not show a flutter/LCO dip. The Navier-Stokes simulations also do not reveal a dip; however, the flutter/LCO onset is at a significantly higher dynamic pressure at Mach 0.90 than at lower Mach numbers. The Navier-Stokes simulations indicate a mild LCO onset at Mach 0.82, then a more rapidly growing instability at Mach 0.86 and 0.90. Finally, the modeling issues and their solution related to the use of a beam and pod finite element model to generate the Navier-Stokes structure mode shapes are discussed.

Bartels, Robert E.↗

Computational Approaches to Vestibular Research

The Biocomputation Center at NASA Ames Research Center is dedicated to a union between computational, experimental and theoretical approaches to the study of neuroscience and of life sciences in general. The current emphasis is on computer reconstruction and visualization of vestibular macular architecture in three-dimensions (3-D), and on mathematical modeling and computer simulation of neural activity in the functioning system. Our methods are being used to interpret the influence of spaceflight on mammalian vestibular maculas in a model system, that of the adult Sprague-Dawley rat. More than twenty 3-D reconstructions of type I and type II hair cells and their afferents have been completed by digitization of contours traced from serial sections photographed in a transmission electron microscope. This labor-intensive method has now been replace d by a semiautomated method developed in the Biocomputation Center in which conventional photography is eliminated. All viewing, storage and manipulation of original data is done using Silicon Graphics workstations. Recent improvements to the software include a new mesh generation method for connecting contours. This method will permit the investigator to describe any surface, regardless of complexity, including highly branched structures such as are routinely found in neurons. This same mesh can be used for 3-D, finite volume simulation of synapse activation and voltage spread on neuronal surfaces visualized via the reconstruction process. These simulations help the investigator interpret the relationship between neuroarchitecture and physiology, and are of assistance in determining which experiments will best test theoretical interpretations. Data are also used to develop abstract, 3-D models that dynamically display neuronal activity ongoing in the system. Finally, the same data can be used to visualize the neural tissue in a virtual environment. Our exhibit will depict capabilities of our computational approaches and some of our findings from their application. For example, our research has demonstrated that maculas of adult mammals retain the property of synaptic plasticity. Ribbon synapses increase numerically and undergo changes in type and distribution (p<0.0001) in type II hair cells after exposure to microgravity for as few as nine days. The finding of macular synaptic plasticity is pertinent to the clinic, and may help explain some. balance disorders in humans. The software used in our investigations will be demonstrated for those interested in applying it in their own research.

Ross, Muriel D.↗

Control of linear uncertain systems utilizing mismatched state observers

The control of linear continuous dynamical systems is investigated as a problem of limited state feedback control. The equations which describe the structure of an observer are developed constrained to time-invarient systems. The optimal control problem is formulated, accounting for the uncertainty in the design parameters. Expressions for bounds on closed loop stability are also developed. The results indicate that very little uncertainty may be tolerated before divergence occurs in the recursive computation algorithms, and the derived stability bound yields extremely conservative estimates of regions of allowable parameter variations.

Goldstein, B.↗