Parameter plane analysis of forced nonlinear oscillations.
Parameter plane analysis of forced oscillations and jump resonance phenomena in nonlinear systems with periodic forced signals
SEARCH · Search NASA
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.
Parameter plane analysis of forced oscillations and jump resonance phenomena in nonlinear systems with periodic forced signals
Parameter plane analysis of forced oscillations and jump resonance phenomena in nonlinear systems
Numerical solution of nonlinear equilibrium problems of structures by means of Newton-Raphson type iterations is reviewed. Each step of the iteration is shown to correspond to the solution of a linear problem, therefore the feasibility of the finite element method for nonlinear analysis is established. Organization and flow of data for various types of digital computers, such as single-processor/single-level memory, single-processor/two-level-memory, vector-processor/two-level-memory, and parallel-processors, with and without sub-structuring (i.e. partitioning) are given. The effect of the relative costs of computation, memory and data transfer on substructuring is shown. The idea of assigning comparable size substructures to parallel processors is exploited. Under Cholesky type factorization schemes, the efficiency of parallel processing is shown to decrease due to the occasional shared data, just as that due to the shared facilities.
Minimum mass designs are obtained for a simple box beam structure subject to bending, torque and combined bending/torque load cases. These designs are obtained subject to point strain and linear buckling constraints. The present work differs from previous efforts in that special attention is payed to including the effects of subcomponent panel interaction in the optimal design process. Two different approaches are used to impose the buckling constraints. When the global approach is used, buckling constraints are imposed on the global structure via a linear eigenvalue analysis. This approach allows the subcomponent panels to interact in a realistic manner. The results obtained using this approach are compared to results obtained using a traditional, less expensive approach, called the local approach. When the local approach is used, in-plane loads are extracted from the global model and used to impose buckling constraints on each subcomponent panel individually. In the global cases, it is found that there can be significant interaction between skin, spar, and rib design variables. This coupling is weak or nonexistent in the local designs. It is determined that weight savings of up to 7% may be obtained by using the global approach instead of the local approach to design these structures. Several of the designs obtained using the linear buckling analysis are subjected to a geometrically nonlinear analysis. For the designs which were subjected to bending loads, the innermost rib panel begins to collapse at less than half the intended design load and in a mode different from that predicted by linear analysis. The discrepancy between the predicted linear and nonlinear responses is attributed to the effects of the nonlinear rib crushing load, and the parameter which controls this rib collapse failure mode is shown to be the rib thickness. The rib collapse failure mode may be avoided by increasing the rib thickness above the value obtained from the (linear analysis based) optimizer. It is concluded that it would be necessary to include geometric nonlinearities in the design optimization process if the true optimum in this case were to be found.
We develop a theory for the experimentally observed formation of multiple quantum dots (QDs) in strained-layer heteroepitaxy based on surface morphological stability analysis of a coherently strained epitaxial thin film on a crystalline substrate. Using a fully nonlinear model of surface morphological evolution that accounts for a wetting potential contribution to the epitaxial film's free energy as well as surface diffusional anisotropy, we demonstrate the formation of multiple QD patterns in self-consistent dynamical simulations of the evolution of the epitaxial film surface perturbed from its planar state. The simulation predictions are supported by weakly nonlinear analysis of the epitaxial film surface morphological stability. We find that, in addition to the Stranski-Krastanow instability, long-wavelength perturbations from the planar film surface morphology can trigger a nonlinear instability, resulting in the splitting of a single QD into multiple QDs of smaller sizes, and predict the critical wavelength of the film surface perturbation for the onset of the nonlinear tip-splitting instability. The theory provides a fundamental interpretation for the observations of “QD pairs” or “double QDs” and other multiple QDs reported in experimental studies of epitaxial growth of semiconductor strained layers and sets the stage for precise engineering of tunable-size nanoscale surface features in strained-layer heteroepitaxy by exploiting film surface nonlinear, pattern forming phenomena.
The various types of local-global finite-element problems point out the need to develop a new generation of software. First, this new software needs to have a complete analysis capability, encompassing linear and nonlinear analysis of 1-, 2-, and 3-dimensional finite-element models, as well as mixed dimensional models. The software must be capable of treating static and dynamic (vibration and transient response) problems, including the stability effects of initial stress, and the software should be able to treat both elastic and elasto-plastic materials. The software should carry a set of optional diagnostics to assist the program user during model generation in order to help avoid obvious structural modeling errors. In addition, the program software should be well documented so the user has a complete technical reference for each type of element contained in the program library, including information on such topics as the type of numerical integration, use of underintegration, and inclusion of incompatible modes, etc. Some packaged information should also be available to assist the user in building mixed-dimensional models. An important advancement in finite-element software should be in the development of program modularity, so that the user can select from a menu various basic operations in matrix structural analysis.
A new type of self-similar solution of ideal magnetohydrodynamics (MHD) in the nonlinear stage of the undular model (k parallel to B) of the magnetic buoyancy instability (the ballooning instability in fusion plasma physics or the Parker instability in astrophysics) is found through MHD simulation and theory. The linear theory developed agrees well with the simulation in the early (linear) stage. The nonlinear stages of the instability in the simulation show the self-similar evolution. One of the solutions obtained from the nonlinear analysis has the characteristics of nonlinear instability in Lagrangian coordinates; the fluid velocity and the Alfven speed on each magnetic loop increase exponentially with time, because the loop is evacuated by the field-aligned motion of matter resulting from gravitational acceleration. In the later stage of the nonlinear evolution, the solution property changes from exponential to power-law time dependence. The latter corresponds to a force-free expansion solution. The later saturation of the velocity increment is also discussed.
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.
Most published work relating to dynamically loaded journal bearings are directed to determining the minimum film thickness from the predicted journal trajectories. These do not give any information about the subsynchronous whirl stability of journal bearing systems since they do not consider the equations of motion. It is, however, necessary to know whether the bearing system operation is stable or not under such an operating condition. The stability characteristics of the system are analyzed. A linearized perturbation theory about the equilibrium point can predict the threshold of stability; however it does not indicate postwhirl orbit detail. The linearized method may indicate that a bearing is unstable for a given operating condition whereas the nonlinear analysis may indicate that it forms a stable limit cycle. For this reason, a nonlinear transient analysis of a rigid rotor supported on oil journal bearings under: (1) a unidirectional constant load, (2) a unidirectional periodic load, and (3) variable rotating load are performed. The hydrodynamic forces are calculated after solving the time-dependent Reynolds equation by a finite difference method with a successive overrelaxation scheme. Using these forces, equations of motion are solved by the fourth-order Runge-Kutta method to predict the transient behavior of the rotor. With the aid of a high-speed digital computer and graphics, the journal trajectories are obtained for several different operating conditions.
The focus of this research is on the development of analysis and sensitivity analysis equations for nonlinear, transient heat transfer problems modeled by p-version, time discontinuous finite element approximation. The resulting matrix equation of the state equation is simply in the form ofA(x)x = c, representing a single step, time marching scheme. The Newton-Raphson's method is used to solve the nonlinear equation. Examples are first provided to demonstrate the accuracy characteristics of the resultant finite element approximation. A direct differentiation approach is then used to compute the thermal sensitivities of a nonlinear heat transfer problem. The report shows that only minimal coding effort is required to enhance the analysis code with the sensitivity analysis capability.
The prediction of the ultimate load carrying capability for compressively loaded shell structures is a challenging nonlinear analysis problem. Selected areas of finite element technology research and nonlinear solution technology are assessed. Herein, a finite element analysis procedure is applied to four shell collapse problems which have been used by computational structural mechanics researchers in the past. This assessment will focus on a number of different shell element formulations and on different approaches used to account for geometric nonlinearities. The results presented confirm that these aspects of nonlinear shell analysis can have a significant effect on the predicted nonlinear structural response. All analyses were performed using the CSM Testbed software system which allowed a convenient assessment of different element formulations with a consistent approach to solving the discretized nonlinear equations.
The existence of standard forms of finite element equations for nonlinear analysis is observed. These are noted to provide a sufficiently special circumstance to admit the construction of asymptotic solution techniques. Although efforts in this direction have been made only recently, a number of distinct approaches to asymptotic solution techniques for finite element analysis have been put forward and are examined herein. Consideration is given to the differences and similarities of the approaches and the resulting methods of solution. Special characteristics regarding applicability, accuracy and efficiency are delineated.
The solar magnetic activity is governed by a complex dynamo mechanism and exhibits a nonlinear dissipation behavior in nature. The chaotic and fractal properties of solar time series are of great importance to understanding the solar dynamo actions, especially with regard to the nonlinear dynamo theories. In the present work, several nonlinear analysis approaches are proposed to investigate the nonlinear dynamical behavior of the polar faculae and sunspot activity for the time interval from 1951 August to 1998 December. The following prominent results are found: (1) both the high- and the low-latitude solar activity are governed by a three-dimensional chaotic attractor, and the chaotic behavior of polar faculae is the most complex, followed by that of the sunspot areas, and then the sunspot numbers; (2) both the high- and low-latitude solar activity exhibit a high degree of persistent behavior, and their fractal nature is due to such long-range correlation; (3) the solar magnetic activity cycle is predictable in nature, but the high-accuracy prediction should only be done for short- to mid-term due to its intrinsically dynamical complexity. With the help of the Babcock–Leighton dynamo model, we suggest that the nonlinear coupling of the polar magnetic fields with strong active-region fields exhibits a complex manner, causing the statistical similarities and differences between the polar faculae and the sunspot-related indicators.
A new computational approach for transient nonlinear thermal analysis of structures is proposed. It is a hybrid approach which combines the modeling versatility of contemporary finite elements in conjunction with transform methods and classical Bubnov-Galerkin schemes. The present study is limited to nonlinearities due to temperature-dependent thermophysical properties. Numerical test cases attest to the basic capabilities and therein validate the transfinite element approach by means of comparisons with conventional finite element schemes and/or available solutions.
Although various adaptive control techniques have been around for a long time and in spite of successful proofs of stability and even successful demonstrations of performance, the eventual use of adaptive control methodologies in practical real world systems has met a rather strong resistance from practitioners and has remained limited. Apparently, it is difficult to guarantee or even understand the conditions that can guarantee stable operations of adaptive control systems under realistic operational environments. Besides, it is difficult to measure the robustness of adaptive control system stability and allow it to be compared with the common and widely used measure of phase margin and gain margin that is utilized by present, mainly LTI, controllers. Furthermore, customary stability analysis methods seem to imply that the mere stability of adaptive systems may be adversely affected by any tiny deviation from the pretty idealistic and assumably required stability conditions. This paper first revisits the fundamental qualities of customary direct adaptive control methodologies, in particular the classical Model Reference Adaptive Control, and shows that some of their basic drawbacks have been addressed and eliminated within the so-called Simple Adaptive Control methodology. Moreover, recent developments in the stability analysis methods of nonlinear systems show that prior conditions that were customarily assumed to be needed for stability are only apparent and can be eliminated. As a result, sufficient conditions that guarantee stability are clearly stated and lead to similarly clear proofs of stability. As many real-world applications show, once robust stability of the adaptive systems can be guaranteed, the added value of using Add-On Adaptive Control along with classical Control design techniques is pushing the desired performance beyond any previous limits.
Soil-structure interaction (SSI) analysis in the nuclear industry is currently performed using linear codes that function in the frequency domain. There is a consensus that these frequency-domain codes give reasonably accurate results for low-intensity ground motions that result in almost linear response. For higher intensity ground motions, which may result in nonlinear response in the soil, structure or at the vicinity of the foundation, the adequacy of frequency-domain codes is unproven. Nonlinear analysis, which is only possible in the time domain, is theoretically more appropriate in such cases. These methods are available but are rarely used due to the large computational requirements and a lack of experience with analysts and regulators. This paper presents an assessment of the linear frequency-domain code, SASSI, which is widely used in the nuclear industry, and the time-domain commercial finite-element code, LS-DYNA, for SSI analysis. The assessment involves benchmarking the SSI analysis procedure in LS-DYNA against SASSI for linearly elastic models. After affirming that SASSI and LS-DYNA result in almost identical responses for these models, they are used to perform nonlinear SSI analyses of two structures founded on soft soil. An examination of the results shows that, in spite of using identical material properties, the predictions of frequency- and time-domain codes are significantly different in the presence of nonlinear behavior such as gapping and sliding of the foundation.
Techniques used in simulating an electrical circuit with nonlinear elements for use in computer-aided circuit analysis programs are described. Elements of the circuit include capacitors, resistors, inductors, transistors, diodes, and voltage and current sources (constant or time varying). Simulation features are discussed for dc, ac, and/or transient circuit analysis. Calculations are based on the model approach of formulating the circuit equations. A particular solution of transient analysis for nonlinear storage elements is described.
An explanation for the origin and number of clumps along the equatorial ring of Supernova 1987A has eluded decades of research. Our linear analysis and hydrodynamic simulations of the expanding ring prior to the supernova reveal that it is subject to the Crow instability between vortex cores. Here, the dominant wave number is remarkably consistent with the number of clumps, suggesting that the Crow instability stimulates clump formation. Although the present analysis focuses on linear fluid flow, future nonlinear analysis and the incorporation of additional stellar physics may further elucidate the remnant structure and the evolution of the progenitor and other stars.