Archives
January 2025
Categories
All
|
Back to Blog
Nami Mogharabin and Amin Ghadami
J. Comput. Nonlinear Dynam. Feb 2025, 20(2): 021006 https://doi.org/10.1115/1.4067297 Characterizing the nonlinear behavior of dynamical systems near the stability boundary is a critical step toward understanding, designing, and controlling systems prone to stability concerns. Traditional methods for bifurcation analysis in both experimental systems and large-dimensional models are often hindered either by the absence of an accurate model or by the analytical complexity involved. This paper presents a novel approach that combines the theoretical frameworks of nonlinear reduced-order modeling and stability analysis with advanced machine learning techniques to perform bifurcation analysis in dynamical systems. By focusing on a low-dimensional nonlinear invariant manifold, this work proposes a data-driven methodology that simplifies the process of bifurcation analysis in dynamical systems. The core of our approach lies in utilizing carefully designed neural networks to identify nonlinear transformations that map observation space into reduced manifold coordinates in its normal form where bifurcation analysis can be performed. The unique integration of analytical and data-driven approaches in the proposed method enables the learning of these transformations and the performance of bifurcation analysis with a limited number of trajectories. Therefore, this approach improves bifurcation analysis in model-less experimental systems and cost-sensitive high-fidelity simulations. The effectiveness of this approach is demonstrated across several examples.
0 Comments
Read More
Back to Blog
Complex Modal Synthesis Method for Viscoelastic Flexible Multibody System Described by ANCF1/24/2025 Zuqing Yu, Zhuo Liu, Yu Wang, and Qinglong Tian
J. Comput. Nonlinear Dynam. Mar 2025, 20(3): 031004 https://doi.org/10.1115/1.4067522 The viscoelastic dynamic model of flexible multibody coupled with large rotation and deformation can be described by the absolute nodal coordinate formulation (ANCF). However, with the increase of degrees-of-freedom, the computational cost of viscoelastic multibody systems will be very high. In addition, for nonproportionally viscoelastic flexible multibody systems, the orthogonality and superposition of complex modes only exist in the state space. In this investigation, a systematical procedure of model reduction method for viscoelastic flexible multibody systems described by ANCF is proposed based on the complex modal synthesis method. First, the whole motion process of the system is divided into a series of quasi-static equilibrium configurations. Then the dynamic equation is locally linearized based on the Taylor expansion to obtain the constant tangent stiffness matrix and damping matrix. The initial modes and modal coordinates need to be updated for each subinterval. A modal selection criterion based on the energy judgment is proposed to ensure the energy conservation and accuracy by the minimum number of truncations. Finally, three numerical examples are carried out as verification. Simulation results indicate that the method proposed procedure reduces the system scale and improves the computational efficiency under the premise of ensuring the simulation accuracy.
Back to Blog
N'Gbo N'Gbo, Changpin Li, and Min Cai
J. Comput. Nonlinear Dynam. Mar 2025, 20(3): 031002. https://doi.org/10.1115/1.4067471 This article focuses on investigating fractional Lyapunov exponents for generalized ψψ-fractional differential systems. By employing a new and more suitable definition, we derive an expression for the fractional Lyapunov exponents using the inverse of the Mittag-Leffler function, which depends on the kernel, weight, and order of the considered fractional derivative. We also provide an upper bound for the obtained fractional Lyapunov exponents that is tighter than the one available in existing literature. Finally, experiments conducted on a hyperchaotic 5D system and the well-known Lorenz system serve to illustrate and verify our main results.
Back to Blog
Nonlinear Dynamic Analysis Framework for Slender Structures Using the Modal Rotation Method12/10/2024 Yoshitaka Shizuno, Shuonan Dong, Ryo Kuzuno, Taiki Okada, Shugo Kawashima, Kanjuro Makihara, and Keisuke Otsuka
J. Comput. Nonlinear Dynam. March 2025, 20(2): 021002. https://doi.org/10.1115/1.4067201 Owing to their low induced drag, high-aspect-ratio wings are often applied to aircraft, particularly high-altitude long-endurance (HALE) aircraft. An analytical method that considers geometrical nonlinearity is necessary for the analysis of high-aspect-ratio wings as they tend to undergo large deformations. Nonlinear shell/plate or solid finite element methods are widely used for the static analysis of wing strength. However, an increase in the number of elements drastically increases the computational costs owing to the complexity of wing shapes. The modal rotation method (MRM) can avoid this additional expense by analyzing large deformations based on modes and stiffness matrices obtained from any linear or linearized model. However, MRM has only been formulated as a static analysis method. In this study, a novel modal-based dynamic analysis framework, referred to as dynamic MRM (DMRM), is developed to analyze slender cantilever structures. This paper proposes a method to discretize dynamics by capitalizing on the fact that MRM considers geometrical nonlinearity based on deformed shapes. The proposed method targets slender structures with small strains and large displacements and considers geometrical nonlinearity, but not material nonlinearity. Additionally, a formulation method for the work performed by a follower force is proposed. The energy stored in the structure agreed with the work performed by an external force in each performed simulation. DMRM achieved a 95% reduction in the calculation time compared with a nonlinear plate finite element method in a performed simulation.
Back to Blog
Jamal-Odysseas Maaita and Dimitrios Prousalis
J. Comput. Nonlinear Dynam. January 2025, 20(1): 011008. https://doi.org/10.1115/1.4067010 A nonregular oscillation is not enough to define a system as chaotic. A more in-depth investigation is required to prove the existence of chaotic behavior, which is challenging. Although many scientists use the Lyapunov Characteristic Exponents to detect chaos, it is not the only method. Several scientists have introduced different methods that utilize various properties of dynamical systems. Hidden Attractors may be chaotic or regular. The fact that they have small basins of attraction introduces difficulties in locating and characterizing them. The paper presents four different chaotic indicators based on the evolution of the deviation vectors: the maximal Lyapunov Exponent, the Lyapunov Characteristic Exponents, the Fast Lyapunov Index (FLI), and the Small Alignment Index. It includes their properties and the advantages and disadvantages of each method. Also, it includes the algorithms to calculate them and their implementation in Python. The paper closes with a comparison between the four indices applied to a system with hidden attractors.
Back to Blog
Praveen Kumar
J. Comput. Nonlinear Dynam. December 2024, 19(12): 121010. https://doi.org/10.1115/1.4066728 Nonlinear intermodal coupling based on internal resonances in MEMS resonators has advanced significantly over the past two decades for various real-world applications. In this study, we demonstrate the existence of various three-mode combination internal resonances between the first five flexural modes of electrostatically actuated flexible–flexible beams and dynamic modal interaction between three modes via internal resonance. We first calculate the natural frequencies of the beam as a function of the stiffnesses of the transverse and rotational springs of the flexible supports, utilizing both analytical formulation and finite element analysis (FEA). Following this, we identify six combination internal resonances among the first five modes and use applied DC voltage to validate the exactness of one commensurable internal resonance condition (ω2=ω5−ω4). Subsequently, we studied a detailed forced vibration analysis corresponding to this resonance condition by solving the five-mode coupled governing equations through numerical time integration and the method of multiple scales. The results compellingly exhibit three-mode intermodal coupling among the second, fourth, and fifth modes as a function of excitation amplitude and frequency. Alongside this, intriguing nonlinear phenomena such as threshold behavior, saturation phenomena, and autoparametric instability are observed. Finally, this paper provides a systematic methodology for investigating three-mode combination internal resonances and related nonlinear dynamics, offering significant insights that could be used in observing phonon or mechanical lasing phenomena in MEMS resonators.
Back to Blog
Mohammed F. Daqaq
J. Comput. Nonlinear Dynam. January 2025, 20(1): 014501. https://doi.org/10.1115/1.4066659 The symmetric monostable Duffing oscillator exhibits a superharmonic resonance of order three when excited harmonically at an excitation frequency that is one third its linear natural frequency. In this letter, it is shown that a certain class of periodic excitations can inherently quench the superharmonic resonance of order three. The Fourier series expansion of such excitations yields a harmonic component at the natural frequency whose magnitude can be properly tuned to completely quench the effect of the superharmonic component. Based on this understanding, the parameters of a piecewise periodic function and the modulus of the cosine Jacobi elliptic function are intentionally designed to passively suppress the superharmonic resonance. Such periodic functions can be used to replace single-frequency harmonic excitations whenever the effects of the superharmonic resonance are to be passively mitigated.
Back to Blog
Jean M. Souza, Luciana Loureiro S. Monteiro, Marcelo A. Savi
J. Comput. Nonlinear Dynam. December 2024, 19(12): 121005. https://doi.org/10.1115/1.4066469 Magnetic shape memory alloys (MSMAs) constitute a class of smart materials capable of exhibiting large magnetic field induced strain (MFIS) when subjected to magnetomechanical loadings. Two distinct mechanisms are responsible for the induced strain: martensitic variant reorientation and phase transformation. The martensitic reorientation is the most explored mechanism presenting the advantage to potential provide high-frequency actuation since it does not rely on phase transformation cycles. Despite its capabilities and potential dynamical applications, the dynamical behavior of MSMAs is not extensively explored in the literature that is usually focused on quasi-static behavior. Thereby, the objective of this work is to analyze the nonlinear dynamics of MSMAs. In this regard, an MSMA nonlinear oscillator is investigated, exploiting the system response under different bias magnetic field levels and actuation frequencies. A phenomenological model is employed to describe the MSMA magnetomechanical behavior. Numerical simulations are carried out using the operator split technique together with an iterative process and the fourth-order Runge–Kutta method. Results show that the application of a bias magnetic field can reduce the mean displacement of the system, increasing the oscillation amplitude. Furthermore, the period of oscillation can be modified, even achieving complex behaviors, including chaos. The potential use of MSMAs to dynamical systems is explored showing the possibility to provide adaptive behaviors.
Back to Blog
R. Ju, S. M. Yang, H. Ren, W. Fan, R. C. Ni, and P. Gu
J. Comput. Nonlinear Dynam. Dec 2024, 19(12): 121001 https://doi.org/10.1115/1.4066221 Steady-state rotary periodic responses of mechanisms lead to stress cycling in flexible structures or connecting joints, which in turn can result in structural fatigue. A general approach is developed to study rotary periodic solutions of rigid and flexible mechanisms with large spatial rotations based on the incremental harmonic balance (IHB) method. The challenge in analyzing such dynamic systems emanates from the noncommutativity of the spatial rotation and the nonsuperposition nature of the rotational coordinates. The generally used rotational coordinates, such as Euler angles, cannot be expanded into Fourier series, which prevents direct usage of the IHB method. To overcome the problem, the natural coordinates method and absolute nodal coordinate formulation (ANCF) are used herein for the dynamic modeling of the rigid and flexible bodies, respectively. The absolute positions and gradients are used as generalized coordinates, and rotational coordinates are naturally avoided. Equations of motions of the system are differential-algebraic equations (DAEs), and they are solved by the IHB method to obtain the steady-state rotary periodic solutions. The effectiveness of the proposed approach is verified by the simulation of rigid and flexible examples with spatial rotations. The approach is general and robust, and it has the potential to be further extended for other extensive multibody dynamic systems.
Back to Blog
Xu Dai, József Kövecses, and Marek Teichmann
J. Comput. Nonlinear Dynam. Nov 2024, 19(11): 111005 https://doi.org/10.1115/1.4066329 Contact simulation is essential in modeling mechanical systems. The contact models require accurate geometric information, which is determined through collision detection methods. When the mechanical system includes flexible bodies such as structural components, the dynamic formulation and collision detection can be more challenging, as the geometric boundaries of such components keep changing during the simulation. The floating frame of reference (FFR) formulation is suitable for flexible systems with small deformation. In this work, a stable and efficient dynamic simulation method is introduced for flexible systems with contact based on the FFR formulation. In addition, a curve-based collision detection method is proposed, which is more consistent with the dynamic formulation and more efficient than common existing collision detection methods. Case studies of flexible beams and multibody systems are employed to demonstrate the performance of the proposed dynamic simulation and collision detection methods.
Back to Blog
Paolo Magri, Marco Gadola, Daniel Chindamo, and Giulia Sandrini ASME. J. Comput. Nonlinear Dyn. Oct 2024, 19(10): 101006 https://doi.org/10.1115/1.4066092 Suspension and steering design play a major role in ensuring the correct dynamic behaviour of road vehicles. Passenger cars are especially demanding from this point of view: NVH and ride comfort requirements often collide with active safety-related requirements such as road holding in steady-state conditions and stability in transients. Driving pleasure is also important for market success, therefore accurate steering feedback and predictable handling properties are additional priorities.
Since flexible bushings are used as interface between the suspension arms and the chassis, extra degrees of freedom make the design process a complex task. While the use of a multibody software is common practice in the industry, a dedicated computational tool can be more practical and straightforward, especially when undertaking the design of a new suspension concept ground-up. The paper presents a computational methodology for the design of an independent suspension with the associated kinematic and compliance attributes. Typical elastokinematic properties like toe, camber, wheelbase, and track variations vs tyre forces and moments can be computed by means of a dedicated software tool. A sort of validation was performed either by means of a comparison with a MathWorks Simscape® Multibody based model. Finally, a sensitivity analysis is given as an example. Computationally, the method proposed is intuitively based on the equilibrium equations. The nonlinear equations are then solved with Newton-Raphson algorithm. The method can be also optimized for computational efficiency and is thoroughly described so that the reader can easily replicate it in the desired programming environment.
Back to Blog
Sankalp Tiwari, Junaidvali Shaik, and C. P. Vyasarayani J. Comput. Nonlinear Dynam. Feb 2024, 19(2): 021002 https://doi.org/10.1115/1.4064251 Delayed systems are those in which the present dynamics is governed by what happened in the past. They are encountered in manufacturing, biology, population dynamics, control systems, etc. Determining stability of such systems is an important and difficult problem. In the existing works, stability is determined by assuming the governing differential equation. However, the equation may not be known or difficult to obtain. Unlike existing works, our method determines the linear stability of a delayed system using its response to a few known inputs. In particular, our method does not require or assume the differential equation governing that system. The only system information we use is its largest delay time, and the only assumption we make about the underlying equation is that its coefficients are either constant or time-periodic. Our approach involves giving the first few functions of an orthonormal polynomial basis as input and measuring/computing the corresponding responses to generate a state transition matrix, whose largest eigenvalue determines the stability. We demonstrate our method's correctness, efficacy, and convergence by studying four candidate DDEs with differing features. Importantly, we show that our approach is robust to noise in measurement, thereby establishing its suitability for practical applications.
Back to Blog
Milena Petrini, Lucio Demeio and Stefano Lenci ASME. J. Comput. Nonlinear Dyn. October 2023, 18(10): 101004 https://doi.org/10.1115/1.4063034 In this work we examine the nonlinear dynamics of an inverted pendulum between lateral rebounding barriers. We continue the numerical investigation started in (Demeio et al., 2006, “Response Scenario and Non-Smooth Features in the Nonlinear Dynamics of an Impacting Inverted Pendulum”, ASME J. Comput. Nonlin. Dyn., 1(1), pp. 56-64) by adding the contribution of the second harmonic in the external forcing term. We investigate the behavior of the periodic attractors by bifurcation diagrams with respect to each amplitude and by behavior charts of single attractors in the amplitude parameters plane for fixed frequency. We study the effects of the second harmonic term on the existence domain of each attractor, on local bifurcations and on the changes in the basins of attraction. The behavior of some robust chaotic attractor is also considered. In the evolution of the periodic attractors we have observed that the addition of the second harmonic generates a rich variety of behaviors, such as loss of stability and formation of isolas of periodic orbits. In the case of chaotic attractors, we have studied one attractor at high frequency, ω = 18, and one at low frequency, ω = 3. In the high frequency case we detect a transition from a scattered to a confined attractor, whereas at the lower frequency the chaotic attractor is present over a wide range of the second harmonic’s amplitude. Finally, we extend the investigation of the chaotic attractors by bifurcation diagrams with respect to the frequency.
Back to Blog
Tomoyuki Suzuki, Kenji Hirohata, Yasutaka Ito, Takehiro Hato, Akira Kano J. Comput. Nonlinear Dynam. November 2023, 18(11): 111001. https://doi.org/10.1115/1.4063224 This paper proposes a sparse modeling method for automatically creating a surrogate model for nonlinear time-variant systems from a small number of time-series data. The proposed method is an improvement over a method for sparse identification of nonlinear dynamical systems first proposed in 2016, for application to temperature prediction simulations. The form of the thermal model is constrained by the physical model, and we use three novel machine-learning methods to efficiently estimate the model parameters. We verify the proposed method’s effectiveness using time-series data obtained by thermo-fluid analysis of a power module mounted on a comb-shaped heat sink. The proposed method has potential applications in a wide range of fields where the concept of equivalent circuits is applicable. Because the proposed method requires few data, has high extrapolation accuracy, and is easily interpreted, we expect that design parameters can be fine-tuned and actual loads considered, and that condition-based maintenance can be realized through real-time simulations.
Back to Blog
Eduardo Okabe, Victor Paiva, Luis Silva-Teixeira, Jaime Izuka J. Comput. Nonlinear Dynam. Oct 2023, 18(10): 104501 https://doi.org/10.1115/1.4063222 The industry and the scientific community have shown interest in SCARA (Selective Compliance Assembly Robot Arm) robots due to their high accuracy. In this paper, a two-link SCARA has its end-effector pulled by three cables that generate a triangular-shaped workspace. Moving the end-effector in this region is a relatively straightforward task, but placing the end-effector outside it requires a nonlinear dynamic model and a state-of-the-art controller. To address this problem in a simpler, more efficient and innovative manner, the equations of motion are derived and three reinforcement learning algorithms are employed: Proximal Policy Optimization (the same used by the chatbot ChatGPT), Soft Actor-Critic and Twin Delayed Deep Deterministic Policy Gradient. Three targets outside the triangular workspace are considered and the trained networks have their results compared in terms of displacement error, velocity and standard deviation. The Twin Delayed Deep Deterministic Policy Gradient provides creative trajectories, the Soft Actor-Critic presents better solutions for two out of three targets, while the Proximal Policy Algorithm appears to be the most consistent considering all targets under analysis.
|