skip to main content


Title: Model calibration of locally nonlinear dynamical systems: Extended constitutive relation error with multi-harmonic coefficients
Purpose This paper aims to present an approach for calibrating the numerical models of dynamical systems that have spatially localized nonlinear components. The approach implements the extended constitutive relation error (ECRE) method using multi-harmonic coefficients and is conceived to separate the errors in the representation of the global, linear and local, nonlinear components of the dynamical system through a two-step process. Design/methodology/approach The first step focuses on the system’s predominantly linear dynamic response under a low magnitude periodic excitation. In this step, the discrepancy between measured and predicted multi-harmonic coefficients is calculated in terms of residual energy. This residual energy is in turn used to spatially locate errors in the model, through which one can identify the erroneous model inputs which govern the linear behavior that need to be calibrated. The second step involves measuring the system’s nonlinear dynamic response under a high magnitude periodic excitation. In this step, the response measurements under both low and high magnitude excitation are used to iteratively calibrate the identified linear and nonlinear input parameters. Findings When model error is present in both linear and nonlinear components, the proposed iterative combined multi-harmonic balance method (MHB)-ECRE calibration approach has shown superiority to the conventional MHB-ECRE method, while providing more reliable calibration results of the nonlinear parameter with less dependency on a priori knowledge of the associated linear system. Originality/value This two-step process is advantageous as it reduces the confounding effects of the uncertain model parameters associated with the linear and locally nonlinear components of the system.  more » « less
Award ID(s):
1633608
NSF-PAR ID:
10341182
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Engineering Computations
Volume:
36
Issue:
2
ISSN:
0264-4401
Page Range / eLocation ID:
466 to 490
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. A second‐order accurate, linear numerical method is analyzed for the Landau–Lifshitz equation with large damping parameters. This equation describes the dynamics of magnetization, with a non‐convexity constraint of unit length of the magnetization. The numerical method is based on the second‐order backward differentiation formula in time, combined with an implicit treatment for the linear diffusion term from the harmonic mapping part and explicit extrapolation for the nonlinear terms. Afterward, a projection step is applied to normalize the numerical solution at a point‐wise level. This numerical scheme has shown extensive advantages in the practical computations for the physical model with large damping parameters, which comes from the fact that only a linear system with constant coefficients (independent of both time and the updated magnetization) needs to be solved at each time step, and has greatly improved the numerical efficiency. Meanwhile, a theoretical analysis for this linear numerical scheme has not been available. In this paper, we provide a rigorous error estimate of the numerical scheme, in the discrete norm, under suitable regularity assumptions and reasonable ratio between the time step size and the spatial mesh size. In particular, the projection operation is nonlinear, and a stability estimate for the projection step turns out to be highly challenging. Such a stability estimate is derived in details, which will play an essential role in the convergence analysis for the numerical scheme, if the damping parameter is greater than 3.

     
    more » « less
  2. null (Ed.)
    In the present article, we follow up our recent work on the experimental assessment of two data-driven nonlinear system identification methodologies. The first methodology constructs a single nonlinear-mode model from periodic vibration data obtained under phase-controlled harmonic excitation. The second methodology constructs a state-space model with polynomial nonlinear terms from vibration data obtained under uncontrolled broadband random excitation. The conclusions drawn from our previous work (experimental) were limited by uncertainties inherent to the specimen, instrumentation, and signal processing. To avoid these uncertainties in the present work, we pursued a completely numerical approach based on synthetic measurement data obtained from simulated experiments. Three benchmarks are considered, which feature geometric, unilateral contact, and dry friction nonlinearity, respectively. As in our previous work, we assessed the prediction accuracy of the identified models with a focus on the regime near a particular resonance. This way, we confirmed our findings on the strengths and weaknesses of the two methodologies and derive several new findings: First, the state-space method struggles even for polynomial nonlinearities if the training data is chaotic. Second, the polynomial state-space models can reach high accuracy only in a rather limited range of vibration levels for systems with non-polynomial nonlinearities. Such cases demonstrate the sensitivity to training data inherent in the method, as model errors are inevitable here. Third, although the excitation does not perfectly isolate the nonlinear mode (exciter-structure interaction, uncontrolled higher harmonics, local instead of distributed excitation), the modal properties are identified with high accuracy. 
    more » « less
  3. Schmidt, Dirk ; Schreiber, Laura ; Vernet, Elise (Ed.)
    Uncorrected residual wavefront errors limit the ultimate performance of adaptive optics (AO) systems. We present different contributing factors and techniques to estimate and compensate these wavefront errors in the Keck natural guide star (NGS) AO systems. The error terms include low order static and semi-static aberrations from multiple sources, periodic and random segment piston errors, single-segment low order aberrations, wavefront sensor aliasing, vibrations, calibration drifts, and AO-to-telescope offload related errors. We present the design of a new AO subsystem, a residual wavefront controller (rWFC) to monitor the performance of the AO control loops and the image quality of the AO science instruments and apply the necessary changes to the telescope and AO parameters to minimize the residual wavefront errors. The distributed system consists of components at the telescope, AO bench and the science instruments. A few components of this system are already tested as on-demand standalone tools and will be integrated into a high-level graphical user interface (GUI) to operate the system. The software tool will periodically collect AO telemetry data, perform control loop parameter optimization and update AO parameters such as loop gains, centroid gain, etc. In addition, the system will analyze the science data at the end of each exposure and estimate telescope/AO performance when a bright point source is available in the science field. The benefits of reducing or eliminating the residual wavefront errors have broad implications for optical astronomy. Testing these techniques on a segmented telescope will be extremely useful to the teams developing high contrast AO systems for all extremely large telescopes and future segmented space telescopes. 
    more » « less
  4. The chief objective of this paper is to explore energy transfer mechanism between the sub-systems that are coupled by a nonlinear elastic path. In the proposed model (via a minimal order, two degree of freedom system), both sub-systems are defined as damped harmonic oscillators with linear springs and dampers. The first sub-system is attached to the ground on one side but connected to the second sub-system on the other side. In addition, linear elastic and dissipative characteristics of both oscillators are assumed to be identical, and a harmonic force excitation is applied only on the mass element of second oscillator. The nonlinear spring (placed in between the two sub-systems) is assumed to exhibit cubic, hardening type nonlinearity. First, the governing equations of the two degree of freedom system with a nonlinear elastic path are obtained. Second, the nonlinear differential equations are solved with a semi-analytical (multi-term harmonic balance) method, and nonlinear frequency responses of the system are calculated for different path coupling cases. As such, the nonlinear path stiffness is gradually increased so that the stiffness ratio of nonlinear element to the linear element is 0.01, 0.05, 0.1, 0.5 and 1.0 while the absolute value of linear spring stiffness is kept intact. In all solutions, it is observed that the frequency response curves at the vicinity of resonant frequencies bend towards higher frequencies as expected due to the hardening effect. However, at moderate or higher levels of path coupling (say 0.1, 0.5 and 1.0), additional branches emerge in the frequency response curves but only at the first resonant frequency. This is due to higher displacement amplitudes at the first resonant frequency as compared to the second one. Even though the oscillators move in-phase around the first natural frequency, high amplitudes increase the contribution of the stored potential energy in the nonlinear spring to the total mechanical energy. The out-of-phase motion around the second natural frequency cannot significantly contribute due to very low motion amplitudes. Finally, the governing equations are numerically solved for the same levels of nonlinearity, and the motion responses of both sub-systems are calculated. Both in-phase and out-of-phase motion responses are successfully shown in numerical solutions, and phase portraits of the system are generated in order to illustrate its nonlinear dynamics. In conclusion, a better understanding of the effect of nonlinear elastic path on two damped harmonic oscillators is gained. 
    more » « less
  5. null (Ed.)
    Loss of operation or devastating damage to buildings and industrial structures, as well as equipment housed in them, has been observed due to earthquake-induced vibrations. A common source of operational downtime is due to the performance reduction of vital equipment, which are sensitive to the total transmitted acceleration. A well-known method of protecting such equipment is seismic isolation of the equipment itself (or a group of equipment), as opposed to the entire structure due to the lower cost of implementation. The first objective of this dissertation is assessing a rolling isolation system (RIS) based on existing design guidelines for telecommunications equipment. A discrepancy is observed between the required response spectrum (RRS) and the one and only accelerogram recommended in the guideline. Several filters are developed to generate synthetic accelerograms that are compatible with the RRS. The generated accelerograms are used for probabilistic assessment of a RIS that is acceptable per the guideline. This assessment reveals large failure probability due to displacement demands in excess of the displacement capacity of the RIS. When the displacement demands on an isolation system are in excess of its capacity, impacts result in spikes in transmitted acceleration. Therefore, the second objective of this dissertation is to design impact prevention/mitigation mechanisms. A dual-mode system is proposed where the behavior changes when the displacement exceeds a predefined threshold. A new piecewise optimal control approach is developed and applied to find the best possible mechanism for the region beyond the threshold. By utilizing the designed curves obtained from the proposed optimal control procedure, a Kelvin-Voigt device is tuned for illustrative purposes. On the other hand, the preference for protecting equipment decreases as the earthquake intensity increases. In extreme seismic loading, the response mitigation of the primary structure (i.e., life safety and collapse prevention) is of greater concern than protecting isolated equipment. Therefore, the third objective of this dissertation is to develop an innovative dual-mode system that can behave as equipment isolation under low to moderate seismic loading and passively transition to behave as a vibration absorber for the primary structure under extreme seismic loading. To reduce the computational cost of simulating a large linear elastic structure with nonlinear attachments (i.e., equipment isolation with cubic hardening nonlinearity), a reduced order modeling method is introduced that can capture the behavior of such nonlinear coupled systems. The method is applied to study the feasibility of dual-mode vibration isolation/absorber. To this end, nonlinear transmissibility curves for the roof displacement and isolated mass total acceleration are developed from the steady-state responses of dual-mode systems using the harmonic balanced method. The final objective of this dissertation is to extend the reduced order modeling method developed for linear elastic structure with nonlinear attachment to inelastic structures (without attachments). The new inelastic model condensation (IMC) method uses the modal properties of the full structural model (in the elastic range) to construct a linear reduced order model in conjunction with a hysteresis model to capture the hysteretic inter-story restoring forces. The parameters of these hysteretic forces are easily tuned, in order to fit the inelastic behavior of the condensed structure to that of the full model under a variety of simple loading scenarios. The fidelity of structural models condensed in this way is demonstrated via simulation for different ground motion intensities on three different building structures with various heights. The simplicity, accuracy, and efficiency of this approach could significantly alleviate the computational burden of performance-based earthquake engineering. 
    more » « less