Finite element models can be used to model and predict the hysteresis and energy dissipation exhibited by nonlinear joints in structures. As a result of the nonlinearity, the frequency and damping of a mode is dependent on excitation amplitude, and when the modes remain uncoupled, quasi-static modal analysis has been shown to efficiently predict this behavior. However, in some cases the modes have been observed to couple such that the frequency and damping of one mode is dependent on the amplitude of other modes. To model the interactions between modes, one must integrate the dynamic equations in time, which is several orders of magnitude more expensive than quasi-static analysis. This work explores an alternative where quasi-static forces are applied in the shapes of two or more modes of vibration simultaneously, and the resulting load–displacement curves are used to deduce the effect of other modes on the effective frequency and damping of the mode in question. This methodology is demonstrated on a simple 2D cantilever beam structure with a single bolted joint which exhibits micro-slip nonlinearity over a range of vibration amplitudes. The predicted frequency and damping are compared with those extracted from a few expensive dynamic simulations of the structure, showing that the quasi-static approach produces reasonable albeit highly conservative bounds on the observed dynamics. This framework is also demonstrated on a 3D structure where dynamic simulations are infeasible.
Finite element models can be used to model and predict the hysteresis and energy dissipation exhibited by nonlinear joints in structures. As a result of the nonlinearity, the frequency and damping of a mode is dependent on excitation amplitude, and when the modes remain uncoupled, quasi-static modal analysis has been shown to efficiently predict this behavior. However, in some cases the modes have been observed to couple such that the frequency and damping of one mode is dependent on the amplitude of other modes. To model the interactions between modes, one must integrate the dynamic equations in time, which is several orders of magnitude more expensive than quasi-static analysis. This work explores an alternative where quasi-static forces are applied in the shapes of two or more modes of vibration simultaneously, and the resulting load–displacement curves are used to deduce the effect of other modes on the effective frequency and damping of the mode in question. This methodology is demonstrated on a simple 2D cantilever beam structure with a single bolted joint which exhibits micro-slip nonlinearity over a range of vibration amplitudes. The predicted frequency and damping are compared with those extracted from a few expensive dynamic simulations of the structure, showing that the quasi-static approach produces reasonable albeit highly conservative bounds on the observed dynamics. This framework is also demonstrated on a 3D structure where dynamic simulations are infeasible.
In the context of experimental vibration data, strain gauges can obtain linear and nonlinear dynamic measurements. However, measuring strain can be disincentivizing and expensive due to the complexity of data acquisition systems, lack of portability, and high costs. This research introduces the use of a low-cost efficient wireless intelligent sensor for strain (LEWIS-S) that is based on a portable-sensor-design platform that streamlines strain sensing. Additionally, the softening behavior of a cantilever beam with geometric and inertial nonlinearities is characterized by the LEWIS-S based on high force level inputs. Two experiments were performed on a nonlinear cantilever beam with measurements obtained by the LEWIS-S sensor and an accelerometer. First, a sine sweep test was performed through the fundamental resonance of the system, then a ring-down test was performed from a large initial static deformation. Good agreement was revealed in quantities of interest such as frequency response functions, the continuous wavelet transforms, and softening behavior in the backbone curves.
In this work, we study how a contact/impact nonlinearity interacts with a geometric cubic nonlinearity in an oscillator system. Specific focus is shown to the effects on bifurcation behavior and secondary resonances (i.e., super- and sub-harmonic resonances). The effects of the individual nonlinearities are first explored for comparison, and then the influences of the combined nonlinearities, varying one parameter at a time, are analyzed and discussed. Nonlinear characterization is then performed on an arbitrary system configuration to study super- and sub-harmonic resonances and grazing contacts or bifurcations. Both the cubic and contact nonlinearities cause a drop in amplitude and shift up in frequency for the primary resonance, and they activate high-amplitude subharmonic resonance regions. The nonlinearities seem to never destructively interfere. The contact nonlinearity generally affects the system's superharmonic resonance behavior more, particularly with regard to the occurrence of grazing contacts and the activation of many bifurcations in the system's response. The subharmonic resonance behavior is more strongly affected by the cubic nonlinearity and is prone to multistable behavior. Perturbation theory proved useful for determining when the cubic nonlinearity would be dominant compared to the contact nonlinearity. The limiting behaviors of the contact stiffness and freeplay gap size indicate the cubic nonlinearity is dominant overall. It is demonstrated that the presence of contact may result in the activation of several bifurcations. In addition, it is proved that the system's subharmonic resonance region is prone to multistable dynamical responses having distinct magnitudes.
Dynamical systems subject to intermittent contact are often modeled with piecewise-smooth contact forces. However, the discontinuous nature of the contact can cause inaccuracies in numerical results or failure in numerical solvers. Representing the piecewise contact force with a continuous and smooth function can mitigate these problems, but not all continuous representations may be appropriate for this use. In this work, five representations used by previous researchers (polynomial, rational polynomial, hyperbolic tangent, arctangent, and logarithm-arctangent functions) are studied to determine which ones most accurately capture nonlinear behaviors including super- and subharmonic resonances, multiple solutions, and chaos. The test case is a single-DOF forced Duffing oscillator with freeplay nonlinearity, solved using direct time integration. This work intends to expand on past studies by determining the limits of applicability for each representation and what numerical problems may occur.
Conference Proceedings of the Society for Experimental Mechanics Series
Saunders, Brian E.; Vasconcellos, Rui M.G.; Kuether, Robert J.; Abdelkefi, Abdessattar
Dynamical systems containing contact/impact between parts can be modeled as piecewise-smooth reduced-order models. The most common example is freeplay, which can manifest as a loose support, worn hinges, or backlash. Freeplay causes very complex, nonlinear responses in a system that range from isolated resonances to grazing bifurcations to chaos. This can be an issue because classical solution methods, such as direct time integration (e.g., Runge-Kutta) or harmonic balance methods, can fail to accurately detect some of the nonlinear behavior or fail to run altogether. To deal with this limitation, researchers often approximate piecewise freeplay terms in the equations of motion using continuous, fully smooth functions. While this strategy can be convenient, it may not always be appropriate for use. For example, past investigation on freeplay in an aeroelastic control surface showed that, compared to the exact piecewise representation, some approximations are not as effective at capturing freeplay behavior as other ones. Another potential issue is the effectiveness of continuous representations at capturing grazing contacts and grazing-type bifurcations. These can cause the system to transition to high-amplitude responses with frequent contact/impact and be particularly damaging. In this work, a bifurcation study is performed on a model of a forced Duffing oscillator with freeplay nonlinearity. Various representations are used to approximate the freeplay including polynomial, absolute value, and hyperbolic tangent representations. Bifurcation analysis results for each type are compared to results using the exact piecewise-smooth representation computed using MATLAB® Event Location. The effectiveness of each representation is compared and ranked in terms of numerical accuracy, ability to capture multiple response types, ability to predict chaos, and computation time.
The goal of this paper is to present a set of measurements from a benchmark structure containing two bolted joints to support future efforts to predict the damping due to the joints and to model nonlinear coupling between the first two elastic modes. Bolted joints introduce nonlinearities in structures, typically causing a softening in the natural frequency and an increase in damping because of frictional slip between the contact interfaces within the joint. These nonlinearities pose significant challenges when characterizing the response of the structure under a large range of load amplitudes, especially when the modal responses become coupled, causing the effective damping and natural frequency to not only depend on the excitation amplitude of the targeted mode, but also the relative amplitudes of other modes. In this work, two nominally identical benchmark structures, known in some prior works as the S4 beam, are tested to characterize their nonlinear properties for the first two elastic modes. Detailed surface measurements are presented and validated through finite element analysis and reveal distinct contact interactions between the two sets of beams. The free-free test structures are excited with an impact hammer and the transient response is analyzed to extract the damping and frequency backbone curves. A range of impact amplitudes and drive points are used to isolate a single mode or to excite both modes simultaneously. Differences in the nonlinear response correlate with the relative strength of the modes that are excited, allowing one to characterize mode coupling. Each of the beams shows different nonlinear properties for each mode, which is attributed to the different contact pressure distributions between the parts, although the mode coupling relationship is found to be consistent between the two. The test data key finding are presented in this paper and the supporting data is available on a public repository for interested researchers.
Conference Proceedings of the Society for Experimental Mechanics Series
Singh, Aabhas; Wielgus, Kayla M.; Dimino, Ignazio; Kuether, Robert J.; Allen, Matthew S.
Morphing wings have great potential to dramatically improve the efficiency of future generations of aircraft and to reduce noise and emissions. Among many camber morphing wing concepts, shape changing fingerlike mechanisms consist of components, such as torsion bars, bushings, bearings, and joints, all of which exhibit damping and stiffness nonlinearities that are dependent on excitation amplitude. These nonlinearities make the dynamic response difficult to model accurately with traditional simulation approaches. As a result, at high excitation levels, linear finite element models may be inaccurate, and a nonlinear modeling approach is required to capture the necessary physics. This work seeks to better understand the influence of nonlinearity on the effective damping and natural frequency of the morphing wing through the use of quasi-static modal analysis and model reduction techniques that employ multipoint constraints (i.e., spider elements). With over 500,000 elements and 39 frictional contact surfaces, this represents one of the most complicated models to which these methods have been applied to date. The results to date are summarized and lessons learned are highlighted.
Nonlinear force appropriation is an extension of its linear counterpart where sinusoidal excitation is applied to a structure with a modal shaker and phase quadrature is achieved between the excitation and response. While a standard practice in modal testing, modal shaker excitation has the potential to alter the dynamics of the structure under test. Previous studies have been conducted to address several concerns, but this work specifically focuses on a shaker-structure interaction phenomenon which arises during the force appropriation testing of a nonlinear structure. Under pure-tone sinusoidal forcing, a nonlinear structure may respond not only at the fundamental harmonic but also potentially at sub- or superharmonics, or it can even produce aperiodic and chaotic motion in certain cases. Shaker-structure interaction occurs when the response physically pushes back against the shaker attachment, producing non-fundamental harmonic content in the force measured by the load cell, even for pure tone voltage input to the shaker. This work develops a model to replicate these physics and investigates their influence on the response of a nonlinear normal mode of the structure. Experimental evidence is first provided that demonstrates the generation of harmonic content in the measured load cell force during a force appropriation test. This interaction is replicated by developing an electromechanical model of a modal shaker attached to a nonlinear, three-mass dynamical system. Several simulated experiments are conducted both with and without the shaker model in order to identify which effects are specifically due to the presence of the shaker. The results of these simulations are then compared to the undamped nonlinear normal modes of the structure under test to evaluate the influence of shaker-structure interaction on the identified system's dynamics.
A proper understanding of the complex physics associated with nonlinear dynamics can improve the accuracy of predictive engineering models and provide a foundation for understanding nonlinear response during environmental testing. Several researchers and studies have previously shown how localized nonlinearities can influence the global vibration modes of a system. This current work builds upon the study of a demonstration aluminum aircraft with a mock pylon with an intentionally designed, localized nonlinearity. In an effort to simplify the identification of the localized nonlinearity, previous work has developed a simplified experimental setup to collect experimental data for the isolated pylon mounted to a stiff fixture. This study builds on these test results by correlating a multi-degree-of-freedom model of the pylon to identify the appropriate model form and parameters of the nonlinear element. The experimentally measured backbone curves are correlated with a nonlinear Hurty/Craig-Bampton (HCB) reduced order model (ROM) using the calculated nonlinear normal modes (NNMs). Following the calibration, the nonlinear HCB ROM of the pylon is attached to a linear HCB ROM of the wing to predict the NNMs of the next-level wing-pylon assembly as a pre-test analysis to better understand the significance of the localized nonlinearity on the global modes of the wing structure.
Conference Proceedings of the Society for Experimental Mechanics Series
Saunders, Brian E.; Vasconcellos, Rui M.G.; Kuether, Robert J.; Abdelkefi, Abdessattar
Physical systems that are subject to intermittent contact/impact are often studied using piecewise-smooth models. Freeplay is a common type of piecewise-smooth system and has been studied extensively for gear systems (backlash) and aeroelastic systems (control surfaces like ailerons and rudders). These systems can experience complex nonlinear behavior including isolated resonance, chaos, and discontinuity-induced bifurcations. This behavior can lead to undesired damaging responses in the system. In this work, bifurcation analysis is performed for a forced Duffing oscillator with freeplay. The freeplay nonlinearity in this system is dependent on the contact stiffness, the size of the freeplay region, and the symmetry/asymmetry of the freeplay region with respect to the system’s equilibrium. Past work on this system has shown that a rich variety of nonlinear behaviors is present. Modern methods of nonlinear dynamics are used to characterize the transitions in system response including phase portraits, frequency spectra, and Poincaré maps. Different freeplay contact stiffnesses are studied including soft, medium, and hard in order to determine how the system response changes as the freeplay transitions from soft contact to near-impact. Particular focus is given to the effects of different initial conditions on the activation of secondary- and isolated-resonance responses. Preliminary results show isolated resonances to occur only for softer-contact cases, regions of superharmonic resonances are more prevalent for harder-contact cases, and more nonlinear behavior occurs for higher initial conditions.
This paper presents a set of tests on a bolted benchmark structure called the S4 beam with a focus on evaluating coupling between the first two modes due to nonlinearity. Bolted joints are of interest in dynamically loaded structures because frictional slipping at the contact interface can introduce amplitude-dependent nonlinearities into the system, where the frequency of the structure decreases, and the damping increases. The challenge to model this phenomenon is even more difficult if the modes of the structure become coupled, violating a common assumption of mode orthogonality. This work presents a detailed set of measurements in which the nonlinearities of a bolted structure are highly coupled for the first two modes. Two nominally identical bolted structures are excited using an impact hammer test. The nonlinear damping curves for each beam are calculated using the Hilbert transform. Although the two structures have different frequency and damping characteristics, the mode coupling relationship between the first two modes of the structures is shown to be consistent and significant. The data is intended as a challenge problem for interested researchers; all data from these tests are available upon request.
The 2020 Nonlinear Mechanics and Dynamics (NOMAD) Research Institute was successfully held from June 15 to July 30, 2020. NOMAD brings together participants with diverse technical backgrounds to work in small teams to cultivate new ideas and approaches in engineering mechanics and dynamics research. NOMAD provides an opportunity for researchers – especially early career researchers - to develop lasting collaborations that go beyond what can be established from the limited interactions at their institutions or at annual conferences. A total of 11 students participated in the seven-week long program held virtually due to the COVID-19 health pandemic. The students collaborated on one of four research projects that were developed by various mentors from Sandia National Laboratories, the University of New Mexico, and other academic and research institutions. In addition to the research activities, the students attended weekly technical seminars, various virtual tours, and socialized at virtual gatherings. At the end of the summer, the students gave a final technical presentation on their research findings. Many of the research discoveries made at NOMAD 2020 are published as proceedings at technical conferences and have direct alignment with the critical mission work performed at Sandia.
This work studies the different types of behavior and inaccuracies that can occur when contact is not adequately accounted for in a dynamical system with freeplay, as the strength of the contact stiffness increases. The MATLAB® ode45 time integration solver, with the built-in Event Location capability, is first validated using past experimental data from a forced Duffing oscillator with freeplay. Next, numerical results utilizing event location are compared to results neglecting event location in order to highlight possible numerical errors and effects on multistable dynamical responses. Inaccuracies tend to occur in two different ways. First, neglecting event location can affect the boundaries between basins of attraction. Second, neglecting event location has little effect on the behaviors of the attractor solutions themselves besides merely resembling poorly converged solutions. Errors are less pronounced at the limits of soft or hard contact stiffness. This study shows the importance of accurately solving piecewise-smooth systems and the existing correlation between the strength of the contact force and possible numerical inaccuracies.
Virtual prototyping in engineering design rely on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.
The industrial approach to nonlinearities in structural dynamics is still very conservative, particularly from an experimental point of view. A demo aluminum aircraft has been equipped with discrete nonlinear elements designed to replicate real-world engine pylon subassemblies to increase awareness on the effects of nonlinearities in design, and understand how these effects can be positively exploited, if properly understood. After some preliminary experiments aimed at understanding the coupled behavior of the aircraft-pylon mockup, it became clear that more in-depth numerical and experimental analyses are required on the pylon subassembly alone. For this paper, experimental data is collected to analyze the nonlinear dynamic behavior of the pylon, leading to better understanding of the subassembly once it connects to the aircraft. The pylon element has three main sources of nonlinearities: (1) geometric nonlinearities of the connecting beam, (2) contact as the beam presses into the tapered block surface and (3) friction in the bolted connections. Backbone curves are generated, which map the evolution of natural frequency and damping ratio with excitation amplitude. Using the experimental data, a low-order nonlinear model is identified to replicate the backbone characteristics and response of the pylon.
In this study, the response of a dynamical system with an increasingly stiff freeplay nonlinearity is investigated. Numerical simulations are carried out on a forced Duffing oscillator with low, intermediate, and hard levels of contact. This study is intended to explore the effects of contact stiffness on the system response as the contact approaches hard impact. The importance of accurately capturing points of contact in the system will also be examined. Preliminary results indicate significant shifts in the frequency-response curves of the system along with reduction in amplitudes. These effects are amplified if the contact is asymmetrically spaced in the system. Further, for soft contact, event detection is not always essential to result accuracy. However, for hard contact, event detection is important to ensure accurate results. Differences with/without event detection range from “poorly converged” appearance to slower transient decay to incorrectly predicting potentially dangerous subharmonic resonances.
Bolted joints often risk failure due to the loss of fastener preload when subjected to dynamic, multiaxial loads. This process is a complex problem that depends on multiple attributes such as loading direction, rate, contact within the threads and the interface, material properties, and many others. Current literature suggests that oscillatory shearing loads appear to be most detrimental to the loss of preload in threaded fasteners. To investigate the effect of less idealized loading conditions, an experimental setup employing a bolted c-beam structure is used to study loss of preload from various initial preloads during harmonic excitation near specific resonant frequencies of the structure. The preload force is measured using bolts equipped with internal strain gauges and the structure is excited at specific modes via sine dwell excitation with an electrodynamic shaker. The experiments were designed to measure loss of preload as a function of excitation duration and strength. A finite element model incorporating a fully-threaded joint is developed in parallel to investigate the effectiveness of each at measuring and predicting bolt loosening.
Structural dynamic finite element models typically use multipoint constraints (MPC) to condense the degrees of freedom (DOF) near bolted joints down to a single node, which can then be joined to neighboring structures with linear springs or nonlinear elements. Scalability becomes an issue when multiple joints are present in a system, because each requires its own model to capture the nonlinear behavior. While this increases the computational cost, the larger problem is that the parameters of the joint models are not known, and so one must solve a nonlinear model updating problem with potentially hundreds of unknown variables to fit the model to measurements. Furthermore, traditional MPC approaches are limited in how the flexibility of the interface is treated (i.e. with rigid bar elements the interface has no flexibility). To resolve this shortcoming, this work presents an alternative approach where the contact interface is reduced to a set of modal DOF which retain the flexibility of the interface and are capable of modeling multiple joints simultaneously. Specifically, system-level characteristic constraint (S-CC) reduction is used to reduce the motion at the contact interface to a small number of shapes. To capture the hysteresis and energy dissipation that is present during microslip of joints, a hysteretic element is applied to a small number of the S-CC Shapes. This method is compared against a traditional MPC method (using rigid bar elements) on a two-dimensional finite element model of a cantilever beam with a single joint near the free end. For all methods, a four-parameter Iwan element is applied to the interface DOF to capture how the amplitude dependent modal frequency and damping change with vibration amplitude.
Lacayo et al. (Mechanical Systems and Signal Processing, 118: 133–157, 2019) recently proposed a fast model updating approach for finite element models that include Iwan models to represent mechanical joints. The joints are defined by using RBE3 averaging constraints or RBAR rigid constraints to tie the contact surface nodes to a single node on each side, and these nodes are then connected with discrete Iwan elements to capture tangential frictional forces that contribute to the nonlinear behavior of the mechanical interfaces between bolted joints. Linear spring elements are used in the remaining directions to capture the joint stiffness. The finite element model is reduced using a Hurty/Craig-Bampton approach such that the physical interface nodes are preserved, and the Quasi-Static Modal Analysis approach is used to quickly predict the effective natural frequency and damping ratio as a function of vibration amplitude for each mode of interest. Model updating is then used to iteratively update the model such that it reproduces the correct natural frequency and damping at each amplitude level of interest. In this paper, Lacayo’s updating approach is applied to the S4 Beam (Singh et al., IMAC XXXVI, 2018) giving special attention to the size and type of the multi-point constraints used to connect the structures, and their effect on the linear and nonlinear modal characteristics.
Transient simulations of linear viscoelastically damped structures require excessive computational resources to directly integrate the full-order finite element model with time-stepping algorithms. Traditional modal reduction techniques are not directly applicable to these systems since viscoelastic materials depend on time and frequency. A more appropriate reduction basis is obtained from the nonlinear, complex eigenvalue problem, whose eigenvectors capture the appropriate kinematics and enable frequency-based mode selection; unfortunately, the computational cost is prohibitive for computing these modes from large-scale engineering models. To address this shortcoming, this work proposes a novel two-tier reduction procedure to reduce the upfront cost of solving the complex, nonlinear eigenvalue problem. The first reduction step reduces the full-order model with real mode shapes linearized about various centering frequencies to capture the kinematics over a full range of viscoelastic material behavior (glassy, rubbery, and glass-transition zones). This tier-one reduction preserves time-temperature superposition and allows the equations to depend parametrically on operating temperature. The second-level reduction then solves the complex, nonlinear eigenmode solutions in the tier-one reduced space about a fixed temperature to further reduce the equations-of-motion. The method is demonstrated on a cantilevered sandwich plate to showcase its accuracy and efficiency in comparison to full-order model predictions.
In computational structural dynamics problems, the ability to calibrate numerical models to physical test data often depends on determining the correct constraints within a structure with mechanical interfaces. These interfaces are defined as the locations within a built-up assembly where two or more disjointed structures are connected. In reality, the normal and tangential forces arising from friction and contact, respectively, are the only means of transferring loads between structures. In linear structural dynamics, a typical modeling approach is to linearize the interface using springs and dampers to connect the disjoint structures, then tune the coefficients to obtain sufficient accuracy between numerically predicted and experimentally measured results. This work explores the use of a numerical inverse method to predict the area of the contact patch located within a bolted interface by defining multi-point constraints. The presented model updating procedure assigns contact definitions (fully stuck, slipping, or no contact) in a finite element model of a jointed structure as a function of contact pressure computed from a nonlinear static analysis. The contact definitions are adjusted until the computed modes agree with experimental test data. The methodology is demonstrated on a C-shape beam system with two bolted interfaces, and the calibrated model predicts modal frequencies with <3% total error summed across the first six elastic modes.
Contact in structures with mechanical interfaces has the ability to significantly influence the system dynamics, such that the energy dissipation and resonant frequencies vary as a function of the response amplitude. Finite element analysis is commonly used to study the physics of such problems, particularly when examining the local behavior at the interfaces. These high fidelity, nonlinear models are computationally expensive to run with time-stepping solvers due to their large mesh densities at the interface, and because of the high expense required to update the tangent operators. Hurty/Craig-Bampton substructuring and interface reduction techniques are commonly utilized to reduce computation time for jointed structures. In the past, these methods have only been applied to substructures rigidly attached to one another, resulting in a linear model. The present work explores the performance of a particular interface reduction technique (system-level characteristic constraint modes) on a nonlinear model with node-to-node contact for a benchmark structure consisting of two c-shape beams bolted together at each end.
Krattiger, Dimitri; Wu, Long; Zacharczuk, Martin; Buck, Martin; Kuether, Robert J.; Allen, Matthew S.; Tiso, Paolo; Brake, Matthew R.W.
The Hurty/Craig-Bampton method in structural dynamics represents the interior dynamics of each subcomponent in a substructured system with a truncated set of normal modes and retains all of the physical degrees of freedom at the substructure interfaces. This makes the assembly of substructures into a reduced-order system model relatively simple, but means that the reduced-order assembly will have as many interface degrees of freedom as the full model. When the full-model mesh is highly refined, and/or when the system is divided into many subcomponents, this can lead to an unacceptably large system of equations of motion. To overcome this, interface reduction methods aim to reduce the size of the Hurty/Craig-Bampton model by reducing the number of interface degrees of freedom. This research presents a survey of interface reduction methods for Hurty/Craig-Bampton models, and proposes improvements and generalizations to some of the methods. Some of these interface reductions operate on the assembled system-level matrices while others perform reduction locally by considering the uncoupled substructures. The advantages and disadvantages of these methods are highlighted and assessed through comparisons of results obtained from a variety of representative linear FE models.
A reduced order modeling capability has been developed to reduce the computational burden associated with time-domain solutions of structural dynamic models with linear viscoelastic materials. The discretized equations-of-motion produce convolution integrals resulting in a linear system with nonviscous damping forces. The challenge associated with the reduction of nonviscously damped, linear systems is the selection and computation of the appropriate modal basis to perform modal projection. The system produces a nonlinear eigenvalue problem that is challenging to solve and requires use of specialized algorithms not readily available in commercial finite element packages. This SAND report summarizes the LDRD discoveries of a reduction scheme developed for monolithic finite element models and provides preliminary investigations to extensions of the method using component mode synthesis. In addition, this report provides a background overview of structural dynamic modeling of structures with linear viscoelastic materials, and provides an overview of a new code capability in Sierra Structural Dynamics to output the system level matrices computed on multiple processors.
The 2017 Nonlinear Mechanics and Dynamics (NOMAD) Research Institute was successfully held from June 19 to July 28, 2017. NOMAD seeks to bring together participants with diverse technical backgrounds to work in small teams to utilize an interactive approach to cultivate new ideas and approaches in engineering . NOMAD provides an opportunity for researchers - especially early career researchers - to develop lasting collaborations that go beyond what can be established from the limited interactions at their institutions or at annual conferences. A total of 17 students from around the world came to Albuquerque, New Mexico to participate in the six - week long program held at the University of New Mexico campus. The students collaborated on one of six research projects that were developed by various mentors from Sandia National Laboratories, academia, and other government laboratories. In addition to the research activities, the students attended weekly technical seminars, toured the National Museum of Nuclear Science & History, and socialized at various off - hour events including an Albuquerque Isotopes baseball game. At the end of the summer, the students gave a final technical presentation o n their research findings that was broadcast via Skype. Many of the research discoveries made at NOMAD are published as proceedings at technical conference s and have direct alignment with the critical mission work performed at Sandia.
The ability to model nonlinear structures subject to random excitation is of key importance in designing hypersonic aircraft and other advanced aerospace vehicles. When a structure is linear, superposition can be used to construct its response to a known spectrum in terms of its linear modes. Superposition does not hold for a nonlinear system, but several works have shown that a system's dynamics can still be understood qualitatively in terms of its nonlinear normal modes (NNMs). This work investigates the connection between a structure's undamped nonlinear normal modes and the spectrum of its response to high amplitude random forcing. Two examples are investigated: a spring-mass system and a clamped-clamped beam modeled within a geometrically nonlinear finite element package. In both cases, an intimate connection is observed between the smeared peaks in the response spectrum and the frequency-energy dependence of the nonlinear normal modes. In order to understand the role of coupling between the underlying linear modes, reduced order models with and without modal coupling terms are used to separate the effect of each NNM's backbone from the nonlinear couplings that give rise to internal resonances. In the cases shown here, uncoupled, single-degree-of-freedom nonlinear models are found to predict major features in the response with reasonable accuracy; a highly inexpensive approximation such as this could be useful in design and optimization studies. More importantly, the results show that a reduced order model can be expected to give accurate results only if it is also capable of accurately predicting the frequency-energy dependence of the nonlinear modes that are excited.
Wu, Long W.; Krattiger, Dimitri K.; Zacharczuk, Martin Z.; Buck, Martin B.; Kuether, Robert J.; Allen, Matthew S.; Brake, Matthew R.W.; Tiso, Paolo T.; Reuss, Pascal R.; Salles, Loic S.
Structural dynamics models with localized nonlinearities can be reduced using Hurty/Craig-Bampton component mode synthesis methods. The interior degrees-of-freedom of the linear subcomponents are reduced with a set of dynamic fixedinterface modes while the static constraint modes preserve the physical coordinates at which the nonlinear restoring forces are applied. For finite element models with a highly refined mesh at the boundary, a secondary modal analysis can be performed to reduce the interface down to a truncated set of local-level characteristic constraint modes. In this research, the cost savings and accuracy of the interface reduction technique are evaluated on a simple example problem involving two elastic blocks coming into contact.
Experimental dynamic substructuring is a means whereby a mathematical model for a substructure can be obtained experimentally and then coupled to a model for the rest of the assembly to predict the response. Recently, various methods have been proposed that use a transmission simulator to overcome sensitivity to measurement errors and to exercise the interface between the substructures; including the Craig-Bampton, Dual Craig-Bampton, and Craig-Mayes methods. This work compares the advantages and disadvantages of these reduced order modeling strategies for two dynamic substructuring problems. The methods are first used on an analytical beam model to validate the methodologies. Then they are used to obtain an experimental model for structure consisting of a cylinder with several components inside connected to the outside case by foam with uncertain properties. This represents an exceedingly difficult structure to model and so experimental substructuring could be an attractive way to obtain a model of the system.
Experimental dynamic substructuring is a means whereby a mathematical model for a substructure can be obtained experimentally and then coupled to a model for the rest of the assembly to predict the response. Recently, various methods have been proposed that use a transmission simulator to overcome sensitivity to measurement errors and to exercise the interface between the substructures; including the Craig-Bampton, Dual Craig-Bampton, and Craig-Mayes methods. This work compares the advantages and disadvantages of these reduced order modeling strategies for two dynamic substructuring problems. The methods are first used on an analytical beam model to validate the methodologies. Then they are used to obtain an experimental model for structure consisting of a cylinder with several components inside connected to the outside case by foam with uncertain properties. This represents an exceedingly difficult structure to model and so experimental substructuring could be an attractive way to obtain a model of the system.
The efficiency of a modal substructuring method depends on the component modes used to reduce each subcomponent model. Methods such as Craig-Bampton have been used extensively to reduce linear finite-element models with thousands or even millions of degrees of freedom down orders of magnitude while maintaining acceptable accuracy. A novel reduction method is proposed here for geometrically nonlinear finite-element models using the fixed-interface and constraint modes of the linearized system to reduce each subcomponent model. The geometric nonlinearity requires an additional cubic and quadratic polynomial function in the modal equations, and the nonlinear stiffness coefficients are determined by applying a series of static loads and using the finite-element code to compute the response. The geometrically nonlinear, reduced modal equations for each subcomponent are then coupled by satisfying compatibility and force equilibrium. This modal substructuring approach is an extension of the Craig-Bampton method and is readily applied to geometrically nonlinear models built directly within commercial finite-element packages. The efficiency of this new approach is demonstrated on two example problems: one that couples two geometrically nonlinear beams at a shared rotational degree of freedom, and another that couples an axial spring element to the axial degree of freedom of a geometrically nonlinear beam. The nonlinear normal modes of the assembled models are compared with those of a truth model to assess the accuracy of the novel modal substructuring approach.
Armand, Jason A.; Lacayo, Robert M.; Gross, Johann G.; Reuss, Pascal R.; Salles, Loic S.; Schwingshackl, Christoph W.; Brake, Matthew R.; Kuether, Robert J.
Broadband impact excitation in structural dynamics is a common technique used to detect and characterize nonlinearities in mechanical systems since it excites many frequencies of a structure at once and can be applied with a variety of boundary conditions. Non-stationary time signals from transient ring-down measurements require time-frequency analysis tools to observe variations in frequency and energy dissipation as the response evolves. This work uses the short-time Fourier transform to estimate the instantaneous frequency and damping ratio from either measured or simulated transient ring-down data. By combining the discrete Fourier transform with an expanding or contracting window function that moves along the time axis, the resulting spectrum is used to estimate the instantaneous frequencies, damping and complex Fourier coefficients. This method is demonstrated on a multi-degree-of-freedom beam with a cubic spring attachment, and investigates the amplitudefrequency dependence in connection to the undamped nonlinear normal modes. A second example shows the results from experiment ring-down response on a beam with a lap joint, and reveals how the system behaves as energy dissipates.
Several reduced-order modeling strategies have been developed to create low-order models of geometrically nonlinear structures from detailed finite element models, allowing one to compute the dynamic response of the structure at a dramatically reduced cost. But, the parameters of these reduced-order models are estimated by applying a series of static loads to the finite element model, and the quality of the reduced-order model can be highly sensitive to the amplitudes of the static load cases used and to the type/number of modes used in the basis. Our paper proposes to combine reduced-order modeling and numerical continuation to estimate the nonlinear normal modes of geometrically nonlinear finite element models. Not only does this make it possible to compute the nonlinear normal modes far more quickly than existing approaches, but the nonlinear normal modes are also shown to be an excellent metric by which the quality of the reduced-order model can be assessed. Hence, the second contribution of this work is to demonstrate how nonlinear normal modes can be used as a metric by which nonlinear reduced-order models can be compared. Moreover, various reduced-order models with hardening nonlinearities are compared for two different structures to demonstrate these concepts: a clamped–clamped beam model, and a more complicated finite element model of an exhaust panel cover.
The objective of the present study is to explore the connection between the nonlinear normal modes of an undamped and unforced nonlinear system and the isolated resonance curves that may appear in the damped response of the forced system. To this end, an energy balance technique is used to predict the amplitude of the harmonic forcing that is necessary to excite a specific nonlinear normal mode. A cantilever beam with a nonlinear spring at its tip serves to illustrate the developments. Furthermore, the practical implications of isolated resonance curves are also discussed by computing the beam response to sine sweep excitations of increasing amplitudes.