APPLICATION OF MULTIPLE REGRESSION ANALYSIS TO IDENTIFICATION TIME-VARYING LINEAR DYNAMIC SYSTEMS
Dynamic response characteristics of time-varying linear systems
SEARCH · Engineering Papers
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.
Dynamic response characteristics of time-varying linear systems
AS&M Inc. performed finite element analysis, with and without discrete damage, of a composite semi-span test article that represents the Boeing 220-passenger transport aircraft composite semi-span test article. A NASTRAN bulk data file and drawings of the test mount fixtures and semi-span components were utilized to generate the baseline finite element model. In this model, the stringer blades are represented by shell elements, and the stringer flanges are combined with the skin. Numerous modeling modifications and discrete source damage scenarios were applied to the test article model throughout the course of the study. This report details the analysis method and results obtained from the composite semi-span study. Analyses were carried out for three load cases: Braked Roll, LOG Down-Bending and 2.5G Up-Bending. These analyses included linear and nonlinear static response, as well as linear and nonlinear buckling response. Results are presented in the form of stress and strain plots. factors of safety for failed elements, buckling loads and modes, deflection prediction tables and plots, and strainage prediction tables and plots. The collected results are presented within this report for comparison to test results.
Soft ferromagnetic alloys are often utilized in electromagnetic applications due to their desirable magnetic properties. In support of these applications, the ferromagnetic alloys are also required to bear mechanical load under various loading and environmental conditions. In this study, a Fe–49Co–2V alloy was dynamically characterized in tension with a Kolsky tension bar and a Drop–Hopkinson bar at various strain rates and temperatures. Dynamic tensile stress–strain curves of the Fe–49Co–2V alloy were obtained at strain rates ranging from 40 to 230 s -1 and temperatures from - 100 to 100 °C. All dynamic tensile stress–strain curves exhibited an initial linear elastic response to an upper yield followed by Lüders band response and then a nearly linear work-hardening behavior. The yield strength of this material was found to be sensitive to both strain rate and temperature, whereas the hardening rate was independent of strain rate or temperature. The Fe–49Co–2V alloy exhibited a feature of brittle fracture in tension under dynamic loading with no necking being observed.
Interactions between a flexible beam and a fluid in a channel are of great relevance to biological hairy surfaces, aquatic vegetation, marine life (e.g., fish gills), and many industrial systems alike. While steady state response of a beam to such flows is fairly well-explored, their behavior in the transient regime is not fully understood. Here, a series of numerical simulations are performed to study the laminar Couette flow of an incompressible viscous fluid past an elastic beam in a two-dimensional channel. The flexible beam is perpendicular to the direction of flow, and its base is fixed to the stationary bottom of the channel. We measure the evolution of the Vogel exponent, drag reduction, and reconfiguration number during the transient and steady-state response of the fluid–structure system for different geometrical and physical properties. Our benchmark shows a good agreement between numerical and experimental observations. Our results show that the system's steady-state response at different bulk-fluid velocities can be reproduced by investigating the shear flow response during the transient regime. We define a new variable that characterizes the evolution of the local velocity profile in the proximity of the free end of the beam and use that to characterize the transient-regime response. The analysis yields insight into the competing effects of elasticity of the beam and non-linear flow response.
Both the Greenland and Antarctic ice sheets have been melting at an accelerating rate over recent decades. Meltwater from Greenland might be expected to initiate a climate response which is distinct, and perhaps different from, that associated with Antarctic meltwater. Which one might elicit a greater climate response, and what mechanisms are involved? To explore these questions, we apply “Climate Response Functions (CRFs)” to guide a series of meltwater perturbation experiments using a fully-coupled climate model. In both hemispheres, meltwater drives atmospheric cooling, sea-ice expansion, and strengthened Hadley and Ferrel cells. Greenland meltwater induces a slowdown of the Atlantic Meridional Overturning Circulation (AMOC) and a cooling of the subsurface ocean in the northern high-latitudes. Antarctic meltwater, instead, induces a slowdown of the Antarctic Bottom Water formation and a warming of the subsurface ocean around Antarctica. For melt-rates up to 2000 Gt yr −1 , the climate response is rather linear. However, as melt-rates increase to 5000 Gt yr −1 , the climate response becomes non-linear. Due to a collapsed AMOC, the climate response is super-linear at high Greenland melt-rates. Instead, the climate response is sub-linear at high Antarctic melt-rates, due to the halting of the northward expansion of Antarctic sea ice by warm surface waters. Finally, in the linear limit, we use CRFs and linear convolution theory to make projections of important climate parameters in response to meltwater scenarios, which suggest that Antarctic meltwater will become a major driver of climate change, dominating that of Greenland meltwater, as the current century proceeds.
A true airspeed sensor based on the precession of a vortex whistle for sensing airspeeds up to 321.9 km/hr (200 mph). In an attempt to model the complicated fluid mechanics of the vortex precession, three dimensional, inviscid, unsteady, incompressible fluid flow was studied by using the hydrodynamical linearized stability theory. The temporal stability approach was used to derive the relationship between the true airspeed and frequency response. The results show that the frequency response is linearly proportional to the airspeed. A computer program was developed to obtain the numerical solution. Computational results for various parameters were obtained. The designed sensor basically consisted of a vortex tube, a swirler, and a transducer system. A microphone converted the audible tone to an electronic frequency signal. Measurements for both the closed conduit tests and wind tunnel tests were recorded. For a specific flow rate or airspeed, larger exit swirler angles produced higher frequencies. For a smaller cross sectional area in the precessional flow region, the frequency was higher. It was observed that as the airspeed was increased the Strouhal number remained constant.
A summary of recent developments and results in the area of inelastic response of advanced composites obtained with the micromechanics model proposed by Aboudi is presented. The recent developments include investigation of linear viscoelastic response of polymer matrix composites and yielding and elastoplastic response of metal matrix composites. The predictive capability and accuracy of the model is demonstrated through comparison of the predicted results with experimental data and results of finite-element analyses.
Flight control system design and analysis for aircraft rely on mathematical models of the vehicle dynamics. In addition to a six degree of freedom nonlinear simulation, the X-29A flight controls group developed a set of programs that calculate linear perturbation models throughout the X-29A flight envelope. The models include the aerodynamics as well as flight control system dynamics and were used for stability, controllability, and handling qualities analysis. These linear models were compared to flight test results to help provide a safe flight envelope expansion. A description is given of the linear models at three flight conditions and two flight control system modes. The models are presented with a level of detail that would allow the reader to reproduce the linear results if desired. Comparison between the response of the linear model and flight measured responses are presented to demonstrate the strengths and weaknesses of the linear models' ability to predict flight dynamics.
High and low intensity dynamic environments experienced by a spacecraft during launch and on-orbit operations, respectively, induce structural loads and motions, which are difficult to reliably predict. Structural dynamics in low- and mid-frequency bands are sensitive to component interface uncertainty and non-linearity as evidenced in laboratory testing and flight operations. Analytical tools for prediction of linear system response are not necessarily adequate for reliable prediction of mid-frequency band dynamics and analysis of measured laboratory and flight data. A new MATLAB toolbox, designed to address the key challenges of mid-frequency band dynamics, is introduced in this paper. Finite-element models of major subassemblies are defined following rational frequency-wavelength guidelines. For computational efficiency, these subassemblies are described as linear, component mode models. The complete structural system model is composed of component mode subassemblies and linear or non-linear joint descriptions. Computation and display of structural dynamic responses are accomplished employing well-established, stable numerical methods, modern signal processing procedures and descriptive graphical tools. Parametric sensitivity and Monte-Carlo based system identification tools are used to reconcile models with experimental data and investigate the effects of uncertainties. Models and dynamic responses are exported for employment in applications, such as detailed structural integrity and mechanical-optical-control performance analyses.
The three-dimensional, linearized Euler analysis, LINFLUX, is being developed to provide a comprehensive and efficient unsteady aerodynamic scheme for predicting the aeroacoustic and aeroelastic responses of axial-flow turbomachinery blading. LINFLUX couples a near-field, implicit, wave-split, finite-volume solution to far-field acoustic eigensolutions, to predict the aerodynamic responses of a blade row to prescribed structural and aerodynamic excitations. It is applied herein to predict the acoustic responses of a fan exit guide vane (FEGV) to rotor wake excitations. The intent is to demonstrate and assess the LINFLUX analysis via application to realistic wake/blade-row interactions. Numerical results are given for the unsteady pressure responses of the FEGV, including the modal pressure responses at inlet and exit. In addition, predictions for the modal and total acoustic power levels at the FEGV exit are compared with measurements. The present results indicate that the LINFLUX analysis should be useful in the aeroacoustic design process, and for understanding the three-dimensional flow physics relevant to blade-row noise generation and propagation.
A shock spectrum procedure is developed to estimate the peak displacement response of linear flexible rotor-bearing systems sugjected to a step change in unbalance (i.e., a blade loss). A progressive and a retrograde response spectrum are established. These blade loss response spectra are expressed in a unique non-dimensional form and are functions of the modal damping ratio and the ratio of rotor spin speed to modal damped whirl speed. Modal decomposition using complex modes is utilized to make use of the unique feature of the spectra for the calculation of the peak blade loss displacement response of the rotor system. The procedure is applied to three example systems using several modal superposition strategies. The results of each are compared to true peak displacement obtained by a separate transient response program.
Computerized FORTRAN equations simulation of linear and nonlinear circuits for symbolic analysis and transient response
Due to the high computational cost associated with unsteady aeroelastic analysis, state-of-the-art aeroelastic optimizations based on computational fluid dynamics typically ignore critical constraints like flutter and aeroelastic gust response. The linearized frequency-domain method offers an approach for adding high-fidelity flutter constraints to multidisciplinary optimizations at relatively low cost compared to other unsteady computational fluid dynamics methods. In recent work, sensitivities have been implemented for the linearized frequency-domain method in FUN3D. In this work, the linearized frequency-domain method and associated sensitivities are applied to gradient-based aeroelastic optimization with flutter constraints based on computational fluid dynamics. An overview of the flutter constraint formulation and implementation is provided, and then two optimization problems are presented. The first optimization increases the flutter speed of a pitch and plunge airfoil at transonic conditions using the minimal amount of geometric changes. The second optimization minimizes the mass of the AGARD 445.6 wing subject to a flutter constraint.
Due to the high computational cost associated with unsteady aeroelastic analysis, state-of-the-art aeroelastic optimizations based on computational fluid dynamics typically ignore critical constraints like flutter and aeroelastic gust response. The linearized frequency-domain method offers an approach for adding high-fidelity flutter constraints to multidisciplinary optimizations at relatively low cost compared to other unsteady computational fluid dynamics methods. In recent work, sensitivities have been implemented for the linearized frequency-domain method in FUN3D. In this work, the linearized frequency-domain method and associated sensitivities are applied to gradient-based aeroelastic optimization with flutter constraints based on computational fluid dynamics. An overview of the flutter constraint formulation and implementation is provided, and then two optimization problems are presented. The first optimization increases the flutter speed of a pitch and plunge airfoil at transonic conditions using the minimal amount of geometric changes. The second optimization minimizes the mass of the AGARD 445.6 wing subject to a flutter constraint.
The gust response of a coupled hingeless rotor-fuselage system is studied in both hover and forward flight. Each rotor blade undergoes flap bending, lag bending, and torsional deflections. The blades are discretized into beam elements, each with fifteen nodal degrees of freedom. The fuselage is permitted three translational (vertical, longitudinal, and lateral) and two rotational (pitch and roll) degrees of freedom. The formulation considers a three-dimensional gust field wherein each gust velocity component can have an arbitrary variation in space and time. Aerodynamic loads are obtained using quasisteady strip theory. Wake-induced effects are introduced through dynamic inflow modeling. Dynamic stall and reverse flow effects are also included. Equations governing the rotor-fuselage gust response are linearized about the vehicle propulsive trim state and the blade steady-state defected positions, and solved by time integration. The effects of several parameters on the helicopter gust response are illustrated, including dynamic inflow, lag stiffness, forward speed, gust profile, gust penetration rate, and gust velocity direction.
A general method is presented for determining the dynamic torsional/axial response of linear structures composed of either tapered bars or shafts to transient excitations. The method consists of formulating and solving the dynamic problem in the Laplace transform domain by the finite element method and obtaining the response by a numerical inversion of the transformed solution. The derivation of the torsional and axial stiffness matrices is based on the exact solution of the transformed governing equation of motion, and it consequently leads to the exact solution of the problem. The solution permits treatment of the most practical cases of linear tapered bars and shafts, and employs modeling of structures with only one element per member which reduces the number of degrees of freedom involved. The effects of external viscous or internal viscoelastic damping are also taken into account.
The interactions of dark matter with Standard Model particles can be systematically studied in the language of effective field theories. We investigate dark matter interactions with Standard Model particles, including spin-dependent interactions, for direct detection experiments and demonstrate that, although the scattering rate generally depends on multiple types of material response functions, certain linear combinations of these material response functions vanish if the initial and final electronic states share the same Hamiltonian. We also find that several other response functions vanish in parity-symmetric materials, making these systems as simple as isotropic detectors in some respects. Finally, we present the scattering rate for an anisotropic, possibly chiral detector, for generic dark matter-electron spin interactions. These relations reduce the number of independent response functions needed, thereby simplifying the computational complexity for a broad class of dark matter models. Our results provide a complete and efficient toolkit for analyzing electron recoil signals in diverse detector materials.
Abstract The response to Antarctic sea ice loss within a coupled modeling framework is examined in comparison to the response to Arctic sea ice loss and within the context of general greenhouse warming. Sea ice loss responses are found to be linear (particularly in response to Antarctic or global sea ice loss) with respect to the degree of imposed perturbation and additive when perturbations are applied in hemispheres separately and concurrently. Globally, and in the tropical Pacific in particular, Antarctic sea ice loss plays a relatively larger role than Arctic sea ice loss in both the atmosphere and the ocean, within the parameters of our experiments. The pattern of response to Antarctic sea ice loss is also found to more closely resemble that of greenhouse warming, again particularly in the tropics. An extension to multiparameter pattern scaling is developed to include a scaling factor for Antarctic change in addition to those for tropical warming and Arctic sea ice loss. The decomposition is applied to the modeled response to Antarctic sea ice loss to break it down into component partial responses that scale with Antarctic, tropical, and Arctic changes. This reveals the aspects of the response that are directly related to Antarctic change, such as an equatorward intensification of tropical precipitation in the Northern Hemisphere, and those that are modified via the induced changes in the tropics and Arctic, such as Northern Hemisphere temperature change. With this, we hope to gain a deeper understanding of the role of each of these changes for the development of physical mechanisms of the response.