An Energy Closure Criterion for Model Reduction of a Kicked Euler-Bernoulli Beam

MATLAB source code accompanying the forthcoming research article in the Journal of vibration and acoustics: Bhattacharyya Suparno, Cusumano, Joseph P., 2020, "An Energy Closure Criterion for Model Reduction of a Kicked Euler-Bernoulli Beam", Journal of vibration and acoustics. Contact: Bhattacharyya Suparno ( sxb1086@psu.edu), Penn State University. Abstract: Reduced order models (ROMs) can be simulated with lower computational cost while being more amenable to theoretical analysis. Here, we examine the performance of the proper orthogonal decomposition (POD), a data-driven model reduction technique. We show that the accuracy of ROMs obtained using POD depends on the type of data used and, more crucially, on the criterion used to select the number of proper orthogonal modes (POMs) used for the model. Simulations of a simply supported Euler-Bernoulli beam subjected to periodic impulsive loads are used to generate ROMs via POD, which are then simulated for comparison with the full system. We assess the accuracy of ROMs obtained using steady-state displacement, velocity, and strain fields, tuning the spatiotemporal localization of applied impulses to control the number of excited modes in, and hence the dimensionality of, the system's response. We show that conventional variance-based mode selection leads to inaccurate models for sufficiently impulsive loading, and that this poor performance is explained by the energy imbalance on the reduced subspace. Specifically, the subspace of POMs capturing a fixed amount (say, 99.9%) of the total variance underestimates the energy input and dissipated in the ROM, yielding inaccurate reduced-order simulations. This problem becomes more acute as the loading becomes more spatio-temporally localized (more impulsive). Thus, energy closure analysis provides an improved method for generating ROMs with energetics that properly reflect that of the full system, resulting in simulations that accurately represent the system's true behavior. Code description: These two MATLAB codes were used to generate Figs. 5c, 5d, and Figs. 7 of the paper, which show the energy balance error (\(e_{_W}\)) on the reduced subspace, and the associated modeling errors (\(\widehat{e}_d\) and \(\widehat{e}_v\)) against different values of the subspace dimension \(P\). This code can also be used to calculate these errors in various other situations, with different parameter values. Instructions: Download both codes Open "Energy_Balance_Errors_and_Modeling_Errors.app" in MATLAB Execute the app Enter the values of \(\epsilon\) (parameter controlling spatial localization of the forcing function; \(\epsilon \leq 1\)) \(\tau\) (parameter controlling temporal localization of the external forcing; \(\tau \leq T\)) \(T\) (time period of the periodic forcing function) \(c_m\)(material damping coefficent) \(c_v\) (viscous damping coefficent) field: (choose -1,0, or 1 to select the field to be used for generating the proper orthogonal modes; -1 for strain, 0 for velocity, and 1 for displacement field) \(s\) (location on the beam where the amplitude of the forcing is maximum; in the paper \(s=1/\sqrt{2}\)) Press "Run" Check out the file titled "App_interface.JPG" for a better understanding