Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear response”

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 469 records · Page 26

Experimental Calibration of Flow Rule Coefficients for LS-Dyna Mat_213

Composite materials have become a popular choice in aerospace engineering applications due to their low weight, high specific stiffness, high specific energy absorption, and economic synthesis. As use of composite materials has proliferated, the need for accurate simulation under extreme loading conditions is of increasing importance. *MAT_213 is a tabulated plasticity and damage material model in LS-DYNA designed for such use cases. To account for the nonlinear material response, a plasticity algorithm is used, which employs a non-associative flow rule with coefficients requiring experimental calibration. This work outlines an experimental methodology for calibrating these coefficients using the macroscopic stress-strain response and microscopic in-plane strain measurements from coupon scale experiments. The model material used for the development and demonstration of this methodology is a plain weave woven fabric with T300 carbon fiber tows in the warp direction and Kevlar© 49 fiber tows in the weft direction, and an Epon© 828/Epikure© 3223 thermoset epoxy matrix. Coupon-scale tensile testing performed at multiple fiber angles (carbon fibers oriented 0º, 22.5º, 45º, 67.5º, and 90º relative to the loading axis) provides data for multiple calibration processes. The measured plastic strain increments in the principal material coordinate are used to set limits for ratios between flow rule coefficients. The identified ranges are used to bound the search for a set of coefficients, leading to convergence of different coupon responses in the effective stress – effective plastic strain space. Flow rule coefficients found by this methodology are then used as (i) input for LS-DYNA simulations of uniaxial tension coupon tests at varying fiber angles, and (ii) comparison of predicted values for the plastic Poisson’s ratio in global coordinates with the experimentally measured one for all off-axis orientations. These results are quantitatively and qualitatively compared to macroscopic stress-strain curves and experimental strain components. The conformity of the experimentally measured quantities with those derived from the simulations provide a means for assessing the quality of calibrated coefficients.

Polymer Matrix Composites↗

Identifying approximate linear models for simple nonlinear systems

This paper addresses the identification (realization) of approximate linear models from response data for certain nonlinear dynamic systems. Response characteristics for several typical nonlinear joints are analyzed mathematically and represented by series expansions. The parameters of the series expansion are then compared with the modal parameters of a linear model identified by the Eigensystem Realization Algorithm. The agreement of the identified model and the analytically derived representation is excellent for the cases studied. Also laboratory data from a model which exhibited stiffening behavior was analyzed using the Eigensystem Realization algorithm and Fast Fourier Transform. The laboratory experiment demonstrated the ability of the technique to recover the model characteristics using real data.

Horta, L. G.↗

Response of symmetric rectangular composite laminates with nonlinear damping subjected to acoustic loading

Effects of both nonlinear damping and large-deflection are included in the analysis in an attempt to explain the experimental phenomena of aircraft panels excited at high sound pressure levels; that is the broadening and flattening out of the strain response peaks and the increasing of the response frequency. A nonlinear damping model is considered in the analysis using a single-mode approach. Root-mean-square (RMS) maximum deflection, RMS strains, spectral density functions of strain, and equivalent linear frequency are obtained for clamped and simply supported rectangular symmetric laminated composite plates. It is demonstrated that nonlinear damping is the one which causes the broadening of the response peaks, and it has great influence on RMS panel deflection, RMS strains, and frequency.

Mei, C.↗

Linear and nonlinear aspects of snow albedo feedbacks in atmospheric models

Namias' hypothesis, that anomalous snowcover on the eastern side of the North American continent can generate an anomalous east coast low pressure system and an anomalous inland high pressure system, is consistent with the time-averaged anomalous response from a nonlinear, primitive equation channel model with an idealized, flat land-sea arrangement. An attempt to understand and describe this anomalous response in the nonlinear model as a linear response to anomalous diabatic heating was largely unsuccessful, primarily because the anomalous eddy fluxes were also important. This unsuccessful attempt to describe the nonlinear model's time averages by linear theory then motivated several comparisons between linear and nonlinear severely truncated quasi-geostrophic models. It was also found in these models that the eddy fluxes were extremely important for forcing or dissipating the stationary eddies.

Roads, J. O.↗

Dynamic analysis of nonlinear rotor-housing systems

Nonlinear analysis methods are developed which will enable the reliable prediction of the dynamic behavior of the space shuttle main engine (SSME) turbopumps in the presence of bearing clearances and other local nonlinearities. A computationally efficient convolution method, based on discretized Duhamel and transition matrix integral formulations, is developed for the transient analysis. In the formulation, the coupling forces due to the nonlinearities are treated as external forces acting on the coupled subsystems. Iteration is utilized to determine their magnitudes at each time increment. The method is applied to a nonlinear generic model of the high pressure oxygen turbopump (HPOTP). As compared to the fourth order Runge-Kutta numerical integration methods, the convolution approach proved to be more accurate and more highly efficient. For determining the nonlinear, steady-state periodic responses, an incremental harmonic balance method was also developed. The method was successfully used to determine dominantly harmonic and subharmonic responses fo the HPOTP generic model with bearing clearances. A reduction method similar to the impedance formulation utilized with linear systems is used to reduce the housing-rotor models to their coordinates at the bearing clearances. Recommendations are included for further development of the method, for extending the analysis to aperiodic and chaotic regimes and for conducting critical parameteric studies of the nonlinear response of the current SSME turbopumps.

Noah, Sherif T.↗

Design Optimization of Stiffened Panels with Postbuckling Constraints

The funding provided by the grant is used to complete the final stages of the development of a geometrically nonlinear analysis and design capability for the static response of compressively loaded prismatic plate structures. The analysis is based on the nonlinear finite strip method and is applicable for structures, such as stiffened panels or box columns, that can be modeled as assemblages of finite length plate strips. In an effort to reduce the computational cost of the nonlinear finite strip method, thus making it suitable for use in the design optimization environment, reduced basis techniques as described in various references by Noor are used in conjunction with the finite strip method. In addition, an efficient scheme for tracing the nonlinear equilibrium paths through highly nonlinear response curves was implemented. The new scheme, which is referred to as the normal flow algorithm, is based on homotopy methods, and is capable of negotiating highly nonlinear limit point instabilities with a reduced computational cost compared to the popular Ricks/Wempner and Chrisfield algorithms.

Guerdal, Zafer↗

High-Fidelity Nonlinear Analysis of Compression-Loaded Composite Shells

The results of an experimental and analytical study of the effects of initial imperfections on the buckling and postbuckling response of unstiffened thin-walled compression-loaded graphite-epoxy cylindrical shells are presented. The shells considered in the study have four different shell-wall laminates and two different shell-radius-to-thickness ratios. The shell-wall laminates include two different orthotropic laminates and two different quasi-isotropic laminates. The shell-radius-to-thickness ratios include shell-radius-to-thickness ratios equal to 100 and 200. The results identify the effects of traditional and nontraditional initial imperfections on the nonlinear response characteristics and buckling loads of the shells. The traditional imperfections include the geometric shell-wall mid-surface imperfections that are commonly discussed in the literature on thin shell buckling. The nontraditional imperfections include shell-wall thickness variations, local shell-wall ply-gaps associated with the fabrication process, shell-end geometric imperfections, nonuniform applied end loads, and variations in the boundary conditions including the effects of elastic boundary conditions. A high-fidelity nonlinear shell analysis procedure that accurately accounts for the effects of these traditional and nontraditional imperfections on the nonlinear response characteristics and buckling loads of the shells is described. The analysis procedure includes a nonlinear static analysis that predicts the stable response characteristics of the shells, and a nonlinear transient analysis that predicts the unstable response characteristics. The results of a local shell-wall stress analysis used to estimate failure stresses are also described.

Hilburger, Mark W.↗

Computer analysis of ring stiffened shells of revolution

The equations and method of solution for a series of five compatible computer programs for structural analysis of axisymmetric shell structures are presented. These programs, designated as the SRA programs, apply to a common structural model but analyze different modes of structural response. They are: (1) linear asymmetric static response (SRA 100), (2) buckling of linearized asymmetric equilibrium states (SRA 101), (3) nonlinear axisymmetric static response (SRA 200), (4) buckling of nonlinear axisymmetric equilibrium states(SRA 201), and (5) vibrations about nonlinear axisymmetric equilibrium state (SRA 300).

Cohen, G. A.↗

Response of convectively controlled burning to nonlinear disturbances.

A process where energy release rate is proportional to the square root of the Reynolds number of a burning particle or surface is used to demonstrate some distinctive dynamic properties of a disturbed convectively controlled burning process. The instantaneous energy release rates caused by assumed periodic flow field disturbances are numerically evaluated and examined. Correlation coefficients which express the energy released in-phase with the pressure disturbance are evaluated and analytical solutions for these coefficients are derived. The response of the process is shown to be highly sensitive to the harmonic distortion of the disturbance. For some disturbances the energy released in-phase with the pressure disturbance is an order of magnitude larger than that for linear (sinusoidal) disturbances. The results show that harmonic distortion increases the coupling between the burning process and the flow field. This amplifying effect of harmonic distortion is suppressed when a disturbance becomes steep-fronted; therefore, steepening can act to limit the equilibrium amplitudes in disturbed systems. Amplification is also suppressed when the velocity disturbance is phase shifted with respect to the pressure disturbance or when the process is exposed to a high steady velocity.

Heidmann, M. F.↗

Capabilities of electrodynamic shakers when used for mechanical shock testing

The results of a research task to investigate the capabilities of electrodynamic vibrators (shakers) to perform mechanical shock tests are presented. The simulation method employed was that of developing a transient whose shock response spectrum matched the desired shock response spectrum. Areas investigated included the maximum amplitude capabilities of the shaker systems, the ability to control the shape of the resultant shock response spectrum, the response levels induced at frequencies outside the controlled bandwidth, and the nonlinearities in structural response induced by a change in test level.

Keegan, W. B.↗

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar M Rao↗

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar Rao↗

A Nonlinear, Human-Centered Approach to Motion Cueing with a Neurocomputing Solver

This paper discusses the continuation of research into the development of new motion cueing algorithms first reported in 1999. In this earlier work, two viable approaches to motion cueing were identified: the coordinated adaptive washout algorithm or 'adaptive algorithm', and the 'optimal algorithm'. In this study, a novel approach to motion cueing is discussed that would combine features of both algorithms. The new algorithm is formulated as a linear optimal control problem, incorporating improved vestibular models and an integrated visual-vestibular motion perception model previously reported. A control law is generated from the motion platform states, resulting in a set of nonlinear cueing filters. The time-varying control law requires the matrix Riccati equation to be solved in real time. Therefore, in order to meet the real time requirement, a neurocomputing approach is used to solve this computationally challenging problem. Single degree-of-freedom responses for the nonlinear algorithm were generated and compared to the adaptive and optimal algorithms. Results for the heave mode show the nonlinear algorithm producing a motion cue with a time-varying washout, sustaining small cues for a longer duration and washing out larger cues more quickly. The addition of the optokinetic influence from the integrated perception model was shown to improve the response to a surge input, producing a specific force response with no steady-state washout. Improved cues are also observed for responses to a sway input. Yaw mode responses reveal that the nonlinear algorithm improves the motion cues by reducing the magnitude of negative cues. The effectiveness of the nonlinear algorithm as compared to the adaptive and linear optimal algorithms will be evaluated on a motion platform, the NASA Langley Research Center Visual Motion Simulator (VMS), and ultimately the Cockpit Motion Facility (CMF) with a series of pilot controlled maneuvers. A proposed experimental procedure is discussed. The results of this evaluation will be used to assess motion cueing performance.

Telban, Robert J.↗

Nonlinear structural analysis for fiber-reinforced superalloy turbine blades

A computational capability for predicting the nonlinear thermomechanical structural response of fiber-reinforced superalloy (FRS) turbine blades is described. This capability is embedded in a special purpose computer code (COBSTRAN) developed at the NASA Lewis Research Center. Special features of this computational capability include accounting for: fiber/matrix reaction, nonlinear and anisotropic material behavior, complex stress distribution due to local and global heterogeneity, and residual stresses due to initial fabrication and/or inelastic behavior during subsequent missions. Numerical results are presented from analyses of a hypothetical FRS turbine blade subjected to a fabrication process and subsequent mission cycle. The results demonstrate the capabilities of this computational tool to; predict local stress/strain response and capture trends of local nonlinear and anisotropic material behavior, relate the effects of this local behavior to the global response of a multilayered fiber-composite turbine blade, and trace material history from fabrication through successive missions.

Hopkins, D. A.↗

Simulation of multicorrelated random processes using the FFT algorithm

A technique for the digital simulation of multicorrelated Gaussian random processes is described. This technique is based upon generating discrete frequency functions which correspond to the Fourier transform of the desired random processes, and then using the fast Fourier transform (FFT) algorithm to obtain the actual random processes. The main advantage of this method of simulation over other methods is computation time; it appears to be more than an order of magnitude faster than present methods of simulation. One of the main uses of multicorrelated simulated random processes is in solving nonlinear random vibration problems by numerical integration of the governing differential equations. The response of a nonlinear string to a distributed noise input is presented as an example.

Wittig, L. E.↗