Finite-Element-Based Nonlinear Modal Reduction of a Rotating Beam with Large-Amplitude Motion

There has been a growing demand for nonlinear model updating procedures in the structural dynamics community as advanced vehicles have started to incorporate nonlinearity into their designs. Finite element (FE) model updating is difficult for a nonlinear structure for several reasons: there may be many unknown parameters in the FE model so that multiple solutions may exist, and it is very expensive to compute the structure’s nonlinear normal modes from a full FE model that are an excellent metric to use for nonlinear updating. A recent work showed that some of these challenges can be overcome by updating a nonlinear reduced order model (ROM) rather than the full FE model. The updated ROM can be used to compute response statistics, stresses and ultimately life of the structure. However, one drawback of this approach is that one does not gain physical insight into which parameters in the FE model were in error, and so it is difficult to transfer the lessons learned to future FE models. This work explores the feasibility of updating a FE model to correlate with a ROM with a set of known parameters. The target ROM parameters consist of the linear stiffness and the nonlinear stiffness coefficients which have been identified previously by updating the ROM to correlate with experimental measurements. In this process an optimization routine is setup in which the free parameters are those of the FE model, such as boundary stiffness springs, imperfections, pre-stress, etc. The optimization procedure is wrapped around a ROM creation algorithm that iteratively tunes the FE parameters, creates a ROM, and evaluates the ROM parameters with respect to the target ROM in order to minimize the objective function. This work will demonstrate the procedure on a numerical example of a curved beam in order to verify its effectiveness on the nonlinear model updating.

POD reduced-order unstructured mesh modeling applied to 2D and 3D fluid flow

In recent years, autonomously guided parachute systems have been developed with the goal of increasing the landing accuracy of airdrop operations. These precision airdrop systems utilize an onboard global positioning system (GPS) coupled with a guidance, navigation and control (GNC) system to continuously monitor and correct the flight path of the parachute system during descent. A critical component of the GNC system is an algorithm that predicts the response of the system to various aerodynamic forces in real-time that are typically based on a low-order multi-degree of freedom (DOF) dynamics model. Reduced order modeling is the process of representing a complex dynamic system with many DOF with a model having relatively few DOF. The goal of this study is to examine the use of large scale finite element models to construct and calibrate Reduced Order Models (ROM) for parachute systems. The resulting ROM would provide parachute designers with a simple tool to examine parachute behavior under various operating conditions. The use of numerical simulations to calibrate the ROM would help to eliminate the number of physical tests needed to characterize a parachute system’s properties. In this study, a 6 degree of freedom ROM is adopted. The system dynamic properties for the ROM and the resultant forces and moments acting on the ROM are obtained from finite element simulations. The 6 DOF model originally proposed by Cockrell and Doherr and later used by Dobrokhodov is adopted for the ROM. The finite element (FE) code, TENSION, is used to perform the large scale simulations and extract the ROM dynamic properties and resultant forces. Several examples are performed that verify the ROM procedure by comparing ROM predictions to FE simulations for various parachute system trajectories. The performance of the ROM is shown to be excellent for in-plane motion under an oscillating wind field. It is shown, however, that the ROM procedure deteriorates as the FE model is refined by adding more DOF. It is demonstrated that, for highly deformable FE models with many DOF, more sophisticated ROM procedures are needed.

Reduced-order models for the analysis of a vertical rod under parametric excitation

This paper presents a general methodology to compute nonlinear frequency responses of flat structures subjected to large amplitude transverse vibrations, within a finite element context. A reduced-order model (ROM)is obtained by an expansion onto the eigenmode basis of the associated linearized problem, including transverse and in-plane modes. The coefficients of the nonlinear terms of the ROM are computed thanks to a non-intrusive method, using any existing nonlinear finite element code. The direct comparison to analytical models of beams and plates proves that a lot of coefficients can be neglected and that the in-plane motion can be condensed to the transverse motion, thus giving generic rules to simplify theROM. Then, a continuation technique, based on an asymptotic numerical method and the harmonic balance method, is used to compute the frequency response in free (nonlinear mode computation) or harmonically forced vibrations. The whole procedure is tested on a straight beam, a clamped circular plate and a free perforated plate for which some nonlinear modes are computed, including internal resonances. The convergence with harmonic numbers and oscillators is investigated. It shows that keeping a few of them is sufficient in a range of displacements corresponding to the order of the structure’s thickness, with a complexity of the simulated nonlinear phenomena that increase very fast with the number of harmonics and oscillators.

Bilinear modal representations for reduced-order modeling of localized piecewise-linear oscillators

On the Proper Orthogonal Modes and Normal Modes of Continuous Vibration Systems

Torsional vibrations experienced by drill strings can be detrimental to drilling operations. With a goal of understanding torsional vibrations experienced by drill strings and determining means to attenuate undesired vibrations, the authors have studied the effect of adding a sinusoidal modulation to a constant rotation speed of a drill string. A combination of modeling, analysis, and experiments is used to explore the influence of this rotation input modulation on the system response. The drill string is modeled as a modified Jeffcott rotor, which is described by a system with three degrees of freedom. Considering the case of forward whirling of a rotor in continuous contact with a stator, the equations of motion are reduced to a single degree-of-freedom nonlinear oscillator describing the torsional motions. In order to understand the fast time scale and slow time scale components of the motion, the method of direct partitions of motions is used to determine an approximate response to the nonlinear oscillator. The obtained results of the analysis illustrate that with the sinusoidal modulation of the rotor drive speed, the equivalent torsion stiffness can be enhanced and the character of the friction force at the contact can be made smooth. The analyses helps bring forth the stabilizing influence of the added sinusoidal input to the rotor drive speed. Over the considered parameter ranges, the numerical results obtained with the full three degree-of-freedom model and the reduced single degree-of-freedom model are found to be in agreement with each other. Furthermore, the results from these models are found to compare well with those obtained by using the method of direct partition of motions. Experiments with a laboratory scale drill-string arrangement are to be carried out to validate the analytical and numerical findings and further explore the effectiveness of the drive speed modulation on the rotor dynamics. Introduction Slender rotating structures are used in many engineering applications. Drill strings are long rotating slender structures, which are used in drilling operations. A schematic of a rotary drill rig is shown in Fig. 1 (e.g., Liao, Balachandran, Karkoub, and Abdel-Magid, 2011). Drill strings experience different types of vibrations (axial, torsional, and lateral vibrations) that may lead to detrimental failures of a drilling system. Drill-string vibrations have attracted the attention of many researchers, and many models have been developed to understand them. Since drill strings have a large length-to-diameter ratio, typically, the first torsional natural frequency and first lateral natural frequency are close to each other. This frequency proximity and the nature of the system allows for coupling and energy transfer between the associated vibration modes. Here, as a step towards developing further understanding, a drill string is modeled as an extended planar Jeffcott rotor with gravity acting normal to the rotor. Due to the planar motions, no gyroscopic effects are considered. The extended model, which was earlier considered in the work of Vlajic, Liu, Karki, and Balachandran (2014), is shown in Fig. 2. The model is described with three degrees of freedom (DOF), namely, x and y to account for lateral motions, and theta to account for torsional motion. A large number of research efforts have focused on controlling drill-string vibrations by using different feedback control algorithms, which need measurements along the drill string. In this work, the authors propose a different control approach to mitigate whirling motions during continuous rotor-stator contact. This can be compared to a situation of a drill string being in continuous contact with the borehole. The approach is open loop in implementation and this approach is based on adding a high frequency input to the drive speed of the drill string. Given the open-loop nature, the method does not depend on any measurements along the drill string or rotary table, which could be attractive for a practical stand point. Although the current focus is on motions of a rotor with continuous stator contact, it is planned to study stick-slip and non-contact cases in the future. System Modeling, Studies, and Results As previously mentioned, the drill string is modelled as a Jeffcott rotor with three DOF. Proceeding along the lines of the group's prior work (Vlajic, Liu, Karki, and Balachandran, 2014), after some approximations for the continuous rotor-stator contact case, the 3 DOF system is reduced to a single second-order nonlinear differential equation governing the torsional motion. In Fig. 3, for a representative case with a constant drive speed, the torsional state histories obtained for the full model and the reduced model are plotted. It is noted that the final state of the motion is captured quite well. To further analyze the response during forward whirling, the reduced-order model was nondimensionalized and an approximate solution was obtained by using the method of direct partitions. This method allows one to separate slow scale motions from fast scale motions. As discussed, in the group's prior work, this method can be useful to examine whirling motions. The results obtained by directly integrating the full model as well as the reduced-order model have been compared with that obtained by using the method of direct partitions of motions. It is seen that the perturbation analyses is able to provide an approximate solution that compares well with the numerical results obtained through direct integrations. Next, this analysis is used to examine the influence of the addition of secondary frequency component. It is seen that the addition of the secondary frequency, a high-frequency one, helps smooth out the friction coefficient variation with respect to the relative speed at contact. This is found to have a beneficial effect on the motion and helps suppress large-amplitude torsional motions. Stability analysis was also conducted to examine this effect. Concluding Remarks A study has been conducted to examine the influence of sinusoidal modulation of rotor drive speed in a system, wherein the rotor experiences continuous stator contact. A single degree of-freedom reduced model is developed to study the torsion response and it is found that the addition of a high frequency input can be beneficial in attenuating rotor motions. This is believed to be useful for developing open-loop control schemes for attenuating forward whirling motions of drill strings. Results obtained from a lab scale experimental arrangement will be used to examine the effect of this drive speed modulation further and they will be reported in the conference presentation. Future studies could build on the current effort to examine cases of backward whirling. Acknowledgment: The authors would like to gratefully acknowledge the support received from the Qatar National Research Fund for NPRP Project 7-083-2-041, to pursue this collaborative work between the University of Maryland, College Park, MD, USA and Qatar University, Doha, Qatar.

Reduced order models (ROMs) are computational models whose dimension is significantly lower than those obtained through classical numerical discretizations (e.g., finite element, finite difference, finite volume, or spectral methods). Thus, ROMs have been used to accelerate numerical simulations of many query problems, e.g., uncertainty quantification, control, and shape optimization. Projection-based ROMs have been particularly successful in the numerical simulation of fluid flows. In this brief survey, we summarize some recent ROM developments for the quasi-geostrophic equations (QGE) (also known as the barotropic vorticity equations), which are a simplified model for geophysical flows in which rotation plays a central role, such as wind-driven ocean circulation in mid-latitude ocean basins. Since the QGE represent a practical compromise between efficient numerical simulations of ocean flows and accurate representations of large scale ocean dynamics, these equations have often been used in the testing of new numerical methods for ocean flows. ROMs have also been tested on the QGE for various settings in order to understand their potential in efficient numerical simulations of ocean flows. In this paper, we survey the ROMs developed for the QGE in order to understand their potential in efficient numerical simulations of more complex ocean flows: We explain how classical numerical methods for the QGE are used to generate the ROM basis functions, we outline the main steps in the construction of projection-based ROMs (with a particular focus on the under-resolved regime, when the closure problem needs to be addressed), we illustrate the ROMs in the numerical simulation of the QGE for various settings, and we present several potential future research avenues in the ROM exploration of the QGE and more complex models of geophysical flows.

A strongly coupled model reduction of vibro-acoustic interaction

Inverse viscoelastic material characterization using POD reduced-order modeling in acoustic–structure interaction

A discussion regarding reduced-order modelling of inclined elastic and immersed cables under support excitation

In this paper an efficient algorithm is presented for the development of compact and passive macro‐models of electromagnetic devices through the systematic reduction of the order of discrete models for these devices obtained through the use of finite elements. The proposed methodology is founded on a new finite element formulation that casts Maxwell's curl equations in a state‐space form. Such state‐space representations are very compatible with a class of robust model order reduction techniques based on Krylov subspaces. However, the advantage of this compatibility appears to be hindered by the fact that the state matrix of the discretized Maxwellian system is of dimension almost twice that obtained from the finite element approximation of the electromagnetic vector wave equation. It is shown in this paper that the apparent penalty in both memory and computation efficiency can be avoided by a proper selection of the finite element expansion functions used for the discretization of the electromagnetic fields. More specifically, it is shown that the proper selection of expansion functions renders the state‐space form of the Maxwellian system equivalent to the discrete problem obtained from the approximation of the vector wave equation using tangentially continuous vector finite elements. This equivalence is then used to effect Krylov‐based model order reduction directly from the finite element approximation of the vector wave equation. In particular, a passive model order reduction algorithm is used for this purpose. The proposed reduced order macro‐modelling algorithm is demonstrated through its application to a variety of microwave passive components. Copyright © 2002 John Wiley & Sons, Ltd.

Building on previous efforts dedicated to educing the flow motions associated with low-angle sound emission from turbulent jets ,2 , the current paper proposes a strategy for the reduced-order dynamical modelling of these motions. The technique comprises a dynamic Ansatz obtained by Galerkin projection of the Navier-Stokes equations onto a set of orthonormal spatial basis functions; however, rather than directly computing the coefficients of the resulting evolution equation, which is only possible when full flow information is available, the coefficients are here identified. The described approach, though validated here using partial flow information from a numerical simulation, is thus applicable in an experimental context. A key objective of such modelling is to obtain a dynamical model with as few degrees of freedom as possible, this being motivated by the need to perform real-time sensing and control. The number of degrees of freedom is dictated by the number of basis functions deemed necessary for accurate and robust reproduction of the dynamics, and this raises the question regarding what kind of basis functions are best suited to a given problem. The empirical eigenfunctions obtained by Proper Orthogonal Decomposition (hereafter POD) are frequently used, largely on account of the optimality property of this kind of decomposition. POD is used in this work; however, rather than applying it to the full flow solution, it is applied to the reduced complexity, sound-producing flow skeleton identified in the previous studies . These studies showed that the number of modes necessary to capture a given percentage of fluctuation energy of the acousticallyimportant coherent structures is an order of magnitude less than the number necessary to capture the same percentage of fluctuation energy of the turbulence. The POD modes used here can thus be considered as acoustically-optimised. A reduced-order model (twenty degrees of freedom) is constructed that reproduces the dynamics of the essential sound-producing flow motions. The approach combines the flow calibration technique of Perret et al. and Tikhonov regularisation such as proposed by Cordier et al. . The model successfully computes the evolution of the flow motions over a short time horizon (about six convective time units), and a further intrinsic stabilisation technique is then implemented in order to enable long-time prediction. Sound radiation is computed by means of the jittering wavepacket source Ansatz of Cavalieri et al., as per Kerherve et al., and good agreement is obtained with the sound field of the original simulation, validating the modelling approach with respect to sound source dynamics.

The large computational effort required for finite element analysis predictions of the dynamic response of thin aircraft panels in extreme aeroacoustic environments has led to significant research on the development of nonlinear reduced order models (ROMs). ROMs reduce the finite element model to a low-order system of nonlinear modal-like equations which can be integrated in the time domain. Although much more computationally efficient than finite element models, it may still be expensive to build and solve the ROM equations of motion for complex structures owing to the large number of nonlinear terms involved. To further reduce this computational cost, deterministic and stochastic partial linearization methodologies are introduced in this paper, in which the nonlinearity is retained only for the dominant modes of the basis. In the stochastic approach, the solution approximation error introduced by the partial linearization is recognized as an epistemic uncertainty, which is modeled by randomizing the linear stiffness terms of the non-dominant modes. The nonparametric stochastic approach is adopted for this modeling and its uncertainty level is estimated using maximum likelihood estimation. Both methods are verified on the ROM of a 96,000 degrees-of-freedom, 9-bay panel, where significant computational gains are observed for dynamic simulations. The response of the full nonlinear ROM in the verification study is within the 95% confidence bounds of the stochastic partially linearized model.