Differential Equations Of Motion Research Articles

Nonlinear Dynamic Properties of Rigid–Elastic–Liquid Coupled Ball Bearings Considering External Load Excitation

In this study, a novel dynamic model of a ball bearing based on rigid–elastic–liquid (REL) coupling structure considering external load excitation is established, and the dynamic response of its outer ring under external load excitation is examined. First, the differential equation of the journal motion (i.e. the outer ring) is derived based on Newton’s law, and the oil film force induced by the REL coupling structure is determined by using the central finite difference method, planar bending theory of thin ring, and Simpson’s rule. Further, the fourth and fifth-order variable step length Runge–Kutta method is used to solve the kinematic differential equation. Second, the dynamic characteristics of the outer ring of the bearing under different speeds and elastic support structural parameters are examined. Finally, the feasibility of the theoretical model is verified through experiments. The results show that compared with the traditional squeeze film damper (SFD) system, the ball bearing with REL coupling structure is more stable, and the vibration damping effect is better under higher bearing speed, lower squirrel cage stiffness, and smaller elastic ring flexibility coefficient. The above findings provide useful insights into the dynamic characteristics of REL coupled ball bearings under unbalanced excitation, which can serve as a reference for the optimization design of such ball bearings.

Analysis of nonlinear wave propagation within architected materials consisting of nonlinear Timoshenko beam structural elements

In the present work, a full nonlinear Timoshenko beam employing nonlinear shape functions is developed. The extended Hamilton principle is employed for deriving the differential equations of motion and the associated boundary conditions. The general form of the boundary conditions is then utilized to determine the static solution of the beam motion. Using this solution for the deformation and rotation of the beam, the nonlinear shape functions of the beam are identified, which leads to the linear and nonlinear mass and stiffness matrices of the Timoshenko beam element. The nonlinear dispersion diagram incorporating the non-linear corrections is obtained using the Linstedt–Poincaré perturbation method. An analysis of the effect of internal transverse shear and bending on the nonlinear dispersion characteristics of wave propagation in two-dimensional periodic network materials made of nonlinear Timoshenko beams is done. The formulated theory shows that the percentage of correction factor of the nonlinear kinematics versus the linear dynamical behavior is inversely proportional to the frequency amplitude. The shear and extension modes are shown to have the higher effect in the non-linear correction term in comparison to the flexural mode.

Vibrating roller with compacted soil interaction modelling

Introduction. Vibrating rollers are the most common means of compacting soils in construction. The nature of stress development on the contact surface of the roller with the ground depends on the technical characteristics of the vibrating roller (the mass of the roller, the mass of the roller frame, the frequency and driving force of vibrations, the number and characteristics of the roller shock absorbers) and the properties of the soil.Materials and methods. Simulation of the interaction of a vibrating roller with compacted soil was carried out using a three-mass rheological model of the frame-roller-soil system. Differential equations of mass motion in contact and separation modes were solved numerically. To determine the numerical values of the loading time (increase in contact stresses from zero to the maximum value) and the unloading time (decrease in contact stresses from the maximum value to zero), as well as the maximum reaction force of the soil, a computational experiment was conducted on a rheological model. The mass of the vibrating roller module (the mass of the front axle) and the relative driving force were used as independent parameters of the vibrating roller. The coefficients of elastic and viscous resistance of the soil were chosen as independent parameters of the soil. The total number of combinations of factors was 192. The values of the time of loading and unloading of the soil, as well as the maximum strength of the soil reaction, were determined by oscillograms of changes in the strength of the soil reaction over time.Results. Using the STATISTICA program, regression equations, to calculate the numerical values of the loading and unloading time of the soil, as well as the maximum reaction force of the soil and the corresponding values of the reliability coefficients of the multiple approximation, were obtained.Discussion and conclusion. The rheological model reproduces the asymmetric nature of changes in contact stresses during soil compaction by a vibrating roller, observed in experimental stress oscillograms obtained during field experimental studies. The results obtained are important for calculating the depth of stress propagation in the ground and the distribution of stresses in the ground after the passage of a vibrating roller using a wave approach to describing stress propagation in the ground. In the future, it is advisable to conduct a computational experiment with an expanded list of independent parameters of the roller, including the oscillation frequency.

Geometrically nonlinear vibration of toroidal sandwich shells with auxetic honeycomb core under periodic/impulsive pressure

The main objective of this paper is to evaluate the dynamic behaviour of sandwich toroidal shells composed of a core with negative Poisson ratio (auxetic honeycomb material) and nanocomposite enriched coating layers under impulsive/periodic pressures considering the geometrical nonlinearity. The set of nonlinear equilibrium equations based on the first-order shear deformation theory are and the resultant nonlinear differential motion equations are approximately solved by the Galerkin method and 4th-order Runge–Kutta solution, respectively. The achieved responses show that the amplitude of outward deformations is intensively increased when the proposed sandwich shells are subjected to rectangular periodic and impulsive loadings which shows the destructive effects of this type dynamic loading. The performed parametric studies shows that these destructive effects can be controlled by changing the geometrical properties of the honeycomb cells and the volume fractions of the nanocomposites in the coating layers. It is notable that in the previous studies the nonlinear dynamic response of toroidal shells subjected to various periodic and impulsive pressures was not investigated. However, in the present work, the mentioned dynamic responses were achieved and investigated without spending time-consuming computational. In fact, the novelty of this work is to develop such a robust and efficient tool.

Calculation of Rayleigh Damping Coefficients for a Transient Structural Analysis

Direct numerical integration of differential equations of motion is widely used by engineers to describe the behavior of structures under dynamic loading. The method entails directly integrating the motion equations over time. In the direct method, the damping matrix is formed as a linear combination of mass and stiffness matrices multiplied by the Rayleigh damping coefficients α and β, respectively. The Rayleigh damping coefficients have a significant effect on the response of building structures under dynamic loading. Therefore, the design values of the damping coefficients α and β have crucial importance to ensure accurate and reliable results in a dynamic analysis. The paper presents a time domain analysis for a building subjected to seismic excitations using the modal superposition and direct integration methods. The direct method considers the damping properties of building structures by Rayleigh damping coefficients obtained using various approaches. The building's response to seismic load is compared by response spectra. The authors proposed the least conservative approach for calculating the Rayleigh damping coefficients for analyzing a building in the time domain.

A Novel Approach to Linear and Nonlinear Time-History Analysis of Structures: Gauss–Lobatto–Hermite 4-Point (GLH-4P) Method

An efficient time integration scheme has been devised for conducting time-history dynamic analysis of single-degree-of-freedom (SDOF) structural systems. The developed scheme belongs to Runge-Kuta family. The primary innovation of the new formulation lies in the constructive cooperation between the Hermite interpolation and the Gauss–Lobatto integration, linked to address the second-order ordinary differential equation of motion (DEOM). To put the idea into practice, several Hermite interpolators are generated to exploit the information of the internal points of Gauss-Lobatto integrator. The 4-point Gauss–Lobatto formula is employed to numerically integrate acceleration and velocity functions. This novel methodology has been designated the Gauss–Lobatto–Hermite 4-point (GLH-4P) method. GLH-4P furnishes a generally applicable implicit algorithm for conducting both linear and nonlinear analyses of structures subjected to seismic excitation. The advantage of the proposed algorithm over the conventional techniques is its higher accuracy level. Besides, it offers better stability compared with the conventional methods such as the Newmark-β and Wilson-θ methods. Unlike the other methods, GLH-4P can handle high-frequency systems without reducing the step size. Numerical examples reveal the efficacy of the GLH-4P method when juxtaposed against existing methodologies.

Analysis of fractional differential equation to compare the dynamic behavior of a single degree of freedom system models using a shaking table

A mathematical model to describe the dynamic behavior of a single degree of freedom system, SDFS, model in springs and damps assembles in series and parallel models in a laboratory, is presented. Kelvin–Voigt fractional derivative rheological models, KVFRdm, was used to solve the fractional differential equations of motion by numerical methods.

Nonlinear forced vibration of functionally graded hybrid three-phase nanocomposite toroidal shell segments reinforced by carbon nanotubes (CNTs) and graphene nanoplatelets (GPLs)

This study proposes a unique three-phase functionally graded (FG) hybrid nanocomposite material reinforcing the toroidal shell segments to investigate the nonlinear forced vibration using Reddy's higher-order shear deformation theory, von Karman geometrical nonlinearity and Stein-McElman's assumption along with the Hamilton principle. The examined toroidal shell's material consists of polymeric resin, carbon nanotubes (CNTs), and graphene platelets (GPLs) as nanocomposites. The material properties have been derived based on a modified Halpin-Tsai micromechanical model. Four distribution patterns have been studied: uniform (UD), FG-X, FG-O, and FG-V. The continuous distribution of GPLs and CNTs leads to inhomogeneous position-dependent properties throughout the shell thickness. The obtained differential equations of motion have been reduced to ordinary equations using Galerkin's technique. A multi-scales method (MSM) is used to estimate a closed-form solution representing the frequency-amplitude relation, and the state space representation is used along with numerical fourth-order Runge Kutta method (RK4) to obtain the nonlinear dynamic response of the toroidal shell. The accuracy of the current results obtained has been verified by comparing them with the relevant literature and numerical results using (RK4). In addition, the influence of both GPLs and CNTs weight fractions, nanofillers’ distribution types, longitudinal and circumferential wave numbers, elastic foundation parameters, static axial compression load, transverse excitation load, damping ratio, geometrical characteristics of the toroidal shell on the dimensionless natural frequency, nonlinear primary resonance (frequency-amplitude curve), and nonlinear dynamic response are carefully studied. The results found that combining (GPLs & CNTs) into the shell's matrix improves performance, especially in concave shells compared to convex ones. The FG-X pattern reduces peak resonance amplitude and improves natural frequency. A reduction of 12.3 % is observed for convex shells, and 21.5 % for concave shells compared to FG-O.

Approximate Noether theorem and its inverse for nonlinear dynamical systems with approximate nonstandard Lagrangian

The approximate Noether theorem and its inverse theorem for the nonlinear dynamical systems with approximate exponential Lagrangian and approximate power-law Lagrangian are investigated. For each case, the approximate differential equations of motion for the nonlinear dynamical systems with approximate nonstandard Lagrangian are established, the generalized Noether identities are given. The relationship between the approximate Noether symmetries and approximate conserved quantities for the system with approximate nonstandard Lagrangian are established, and the approximate Noether theorems and their inverse theorems are obtained. Two examples are given to illustrate the application of the results.

THE SIMULATION MODEL ALGORITHM TO ASSESS ACTIONS OF UNMANNED COMBAT AERIAL VEHICLES

The increase of combat capabilities of UCAVs brought them to the main range of means to defeat the enemy in both defensive and offensive operations. Therefore, the relevance of issues related to the justification of rational approaches to combat operations planning with the use of unmanned combat aerial vehicles makes it necessary to develop methods to evaluate such operations. This research direction has an important applied value. The article proposes an approach to create a simulation model to evaluate unmanned aerial vehicles actions. The research subject is the planning process of the unmanned combat aerial vehicles use. The work aim is to justify approaches to create ways to evaluate the actions of unmanned combat aerial vehicles based on the simulation model of the researched process. The research uses the discrete-time simulation method. The proposed simulation model essentially simulates two processes: the UAV flight process and the process of the enemy targets detecting. As a rule, the aircrafts flight in space is described by differential equations of motion. This approach requires the implementation of a continuous flow of model time, which leads to significant computational difficulties when building the model. Taking into account the necessary accuracy in the simulation model being created, it is permissible to proceed to a discrete change of the model time with some step comparable to the analyzed processes. This model links together such parameters as the barrage zone scope, the target detection intensity, the number of unmanned combat aerial vehicles, their weapons, and speed. When these parameters are varied, it becomes possible to establish the necessary dependencies between them. This can serve as the basis to develop rational ways of using unmanned combat aerial vehicles under the given conditions. The proposed simulation model can be used in the decision-making support system for combat planning with the use of unmanned combat aerial vehicles.

Enhancing low-orbit vibration energy harvesting by a tri-stable piezoelectric energy harvester with an innovative dynamic amplifier

Piezoelectric energy harvesting faces a primary challenge in effectively capturing low-orbit vibration energy across a broad frequency range. In this paper, we present a tri-stable piezoelectric energy harvester that incorporates a dynamic amplifier (TPEH + DM), specifically designed for efficient collection of low-orbit vibration energy. The TPEH + DM comprises a piezoelectric cantilever beam connected to an innovative dynamic amplifier at its restrained end, which enhances both the rotational and lateral displacement of the piezoelectric cantilever beam simultaneously. The governing coupled differential equations of motion for the system is derived based on the Lagrange equation, and analytical expressions for its steady-state response are obtained using the multi-scale method. The influence of factors such as the mass and the stiffness ratio of the dynamic amplifier on the steady-state dynamic output characteristics of the system is investigated using the theoretical analysis and numerical simulation. The results indicate that TPEH + DM exhibits significantly improved energy harvesting performance compared to TPEH under low-orbit external excitations. The bandwidth of inter-well motion and the TPEH + DM power output may be further increased by suitably modifying the relative stiffness between the cantilever beam and the dynamic amplifier. In addition, we analyze the time-domain behavior of the system’s output voltage using the ode45 solver under various external excitation frequencies and intensities. The results demonstrate that with appropriate adjustments to the mass of the tip magnet and the stiffness ratio of the dynamic amplifier, the proposed TPEH + DM system can harvest energy efficiently across a broad frequency range, even under low-orbit excitations.

Large-amplitude nonlinear free vibrational behavior of the rotating rectangular plate reinforced by functionally graded carbon nanotubes

The present study discusses the nonlinear free vibrational behavior of the rotating rectangular plate reinforced by functionally graded carbon nanotubes (FG-CNTs). Mechanical properties of the plate have been continuously changed along the thickness of the plate by considering four different distribution patterns of CNTs. Based on the first-order shear deformation theory (FSDT) and von Ḱarḿan’s geometrical nonlinearity, the governing equations of the plate are derived using Hamilton’s principle. Galerkin discretization technique is employed to convert the partial differential equations of motion to ordinary differential equations. Afterward, the nonlinear time-dependent equations of the plate were analytically solved using the multiple time scales method. Eventually, the effect of CNT parameters, the geometrical ratios and rotational speed on the nonlinear frequency of the plate have been studied. For verification purposes, the present numerical results have been compared with those available in the literature with evident agreement. The present research may give benchmark solutions and some guidelines for designing FG-CNTs rotating plates.

Joint Resonance Analysis in Multiple Modes of Soft Ferromagnetic Rectangular Thin Plate

In this article, the nonlinear principal and internal resonance properties of a soft ferromagnetic rectangular thin plate are investigated in a magnetic field environment. The nonlinear partial differential equation of motion of a soft ferromagnetic rectangular thin plate is derived under the effect of homogeneous simple harmonic excitation. The system of nonlinear differential equations with multiple degrees of freedom is established by the assumed one-sided fixed trilateral simply support condition using the Galerkin’s method. The system of nonlinear differential equations is solved by the multiscale method to obtain the response of two modes under the simple harmonic force at the principal and internal resonance. The numerical results of the system response show that when the frequency of the simple harmonic force is close to one of the modes (first-order or second-order mode) causing it to resonate, the other mode will also resonate internally. The magnetic field can have an inhibiting effect on the resonant response of the system and also affect the kinematic state of the system. The internal resonance provides a mechanism for transferring energy from a high mode to a lower mode.

A unified impact response analysis for FG-CNTRC and 3D-GrF sandwich composite sector- rectangular coupled plates

This paper proposes a general method for analyzing the dynamic response characteristics of sector- rectangular coupled plates subjected to external impacts. This method combines the method of reverberation-ray matrix (MRRM) with the strip transfer function method, referred to as the strip method of reverberation-ray matrix (S-MRRM). The traditional MRRM has some limitations for the dynamic modeling of the sector structure and also requires the Lévy-type exact solution, leading to simply supported boundary conditions in certain direction of the structure. The present method can not only address the aforementioned limitations but also maintain the efficiency advantage of the wave method for solution. The investigated structure consists of sandwich panels composed of functionally graded carbon nanotube reinforced (FG-CNTRC) panels and the three-dimensional graphene foam (3D-GrF) core. The structure is discretized into several strip elements by the strip transfer function method, and the whole dynamic model of the structure is obtained through matrix assembly. The motion differential equation in matrix form is derived by using Hamilton principle, followed by obtaining the matrix-form wave number equation. Upon deriving the wave solution of the structure based on this equation, the transient and steady-state responses of the structure are obtained by employing MRRM. Convergence analysis is conducted on these calculation results, followed by comparison with FEM results to validate its prediction accuracy. Finally, through a comprehensive parametric study, the impacts of various factors such as CNTs distribution form, CNTs volume fraction, 3D-GrF pore distribution form, 3D-GrF porosity, and sandwich panel thickness ratio on both transient and steady-state responses of these structures were investigated, thereby providing guidance for engineering applications.

Gyroscopic effects of the spinning rotor-blades assembly on dynamic response of offshore wind turbines

An analytical solution for the gyroscopic effect of the spinning rotor-blades assembly on the dynamic response of offshore wind turbines (OWT) is presented. A continuous coupled model is rigorously developed to form the partial differential equations of fore-aft, side-side, and yaw motions. The gyroscopic moments caused by the angular momentum of the spinning rotor-blades assembly are formulated and handled into the three boundary conditions at the nacelle. The procedure for obtaining the operational natural frequencies of the structure including these coupled boundary conditions is developed. Furthermore, a response function for each fore-aft and yaw motions is obtained by solving the equations of motion under a wave load applied in only the fore-aft direction. Finally, the operational natural frequencies are calculated and compared to the idling ones for the considered example of OWT. The response of the structure under different values of the rotor-blades assembly's angular velocity is investigated. The proposed solution in this study unfolds a revelational procedure in capturing the gyroscopic effect in wind turbine industries which also can be a guideline for developing the floating OWT models.

Thermal effects on nonlinear vibration of nonlocal nanobeam embedded in nonlinear elastic medium

Nanotechnology has an impact on our lives in a many ways, from better medical treatments and more efficient energy sources to stronger and lighter materials and advanced electronics and this article presents the implementation of a perturbation method for the vibration analysis of simply supported and clamped–clamped Euler–Bernoulli nanobeams resting on nonlinear elastic foundations in thermal environment using nonlocal elasticity theory. Hamilton's principle is used to construct the differential equation of motion of a nanobeam in conjunction with appropriate boundary conditions. The equations of motion of the Euler–Bernoulli nanobeam are determined using nonlocal elasticity theory. It is shown how thermal loadings affect the vibration of the Euler–Bernoulli nanobeam. The multiple scale method, which is one of the perturbation method, is used to get an approximated solution for the presented system. The effects of temperature, Winkler, Pasternak and nonlinear foundation parameters on the vibration analysis of simply supported and clamped–clamped nanobeams are determined and results are given in tables and graphs.

Vibration Characteristics of a Functionally Graded Viscoelastic Fluid-Conveying Pipe with Initial Geometric Defects under Thermal–Magnetic Coupling Fields

In this study, the Kevin–Voigt viscoelastic constitutive relationship is used to investigate the vibration characteristics and stability of a functionally graded viscoelastic(FGV) fluid-conveying pipe with initial geometric defects under thermal–magnetic coupling fields. First, the nonlinear dimensionless differential equations of motion are derived by applying Timoshenko beam theory. Second, by solving the equilibrium position of the system, the nonlinear term in the differential equations of motion is approximated as the sum of the longitudinal displacement at the current time and longitudinal displacement relative to the position, and the equations are linearized. Third, these equations are discretized using the Galerkin method and are numerically solved under simply supported conditions. Finally, the effects of dimensionless temperature field parameters, dimensionless magnetic field parameters, thermal–magnetic coupling, initial geometric defect types, and the power-law exponent on the complex frequency of the pipe are examined. Results show that increasing the magnetic field intensity enhances the critical velocity of first-order mode instability, whereas a heightened temperature variation reduces the critical velocity of first-order diverge instability. Under thermal–magnetic fields, when the magnetic field intensity and temperature difference are simultaneously increased, their effects on the complex frequency can partially offset each other. Increasing the initial geometric defect amplitude increases the imaginary parts of the complex frequencies; however, for different types of initial geometric defect tubes, it exhibits the most distinct influence only on a certain order.

Linear and Nonlinear Dynamics Responses of an Axially Moving Laminated Composite Plate-Reinforced with Graphene Nanoplatelets

The subharmonic resonances of an axially moving graphene-reinforced laminated composite plate are studied based on the Galerkin and multiscale methods. Graphene nanoplatelets (GPLs) are added into matrix material which acts as the basic layer of the plate, and a graphene-reinforced nanocomposite plate is thus obtained. Different GPL distribution patterns in the laminated plate are considered. The Halpin–Tsai model is selected to predict the physical properties of the nanocomposite. Hamilton’s principle is utilized to conduct the dynamic modeling of the plate and the von Kármán deformation theory is used. The velocity is assumed to be a combination of constant and harmonically varied velocities. The natural frequencies of the linear system with constant velocity can be obtained using the eigenvalues of the coefficient matrix of the ordinary differential equations after the governing partial differential equations of motion are discretized through the Galerkin method. The instability regions of the linear system and the amplitude–frequency relations of the nonlinear system considering the harmonically varied velocity are obtained based on the multiscale analysis. The effect of GPL reinforcement on the subharmonic resonances of the linear and nonlinear systems is analyzed in detail.

Formulation and free vibration analysis of spinning Timoshenko composite beams using nonlinear geometrically exact theory

AbstractIn this study, a set of generalized nonlinear equations of motion for Timoshenko composite shafts is derived using the geometrically exact approach. These partial differential equations of motion are discretized by Galerkin method and the free vibration of an orthotropic spinning shaft with hinged‐hinged boundary conditions is investigated using the multiple‐scale method. We find that in the case of two mode discretization, composite couplings can have significant linear and nonlinear effects on the behavior of the shaft and that these effects were not predictable in the classical model. Also, in the investigating of composite layering we find that the effects of the fibers angle as well as the composite couplings can be noticeable, and these effects are greater in the case of the asymmetrical layering.

Structural design and dynamic characteristics analysis of braided composite two-stage gear transmission system

In order to realize the lightweight design of the transmission system, the braided composite material is applied to the two-stage gear transmission system. According to the structural characteristics of the two-stage gear reducer box, the whole box is designed to be assembled with the braided base and the box wall. Woven composite materials are applied to the web parts of mixed metal composite gears to realize the design goal of lightweight gears. Then, under the assumption of ignoring the influence of friction, bearings and other factors on the system, the dynamic model of the two-stage gear transmission system considering the box is established. By normalizing and dimensionless processing of the equations, the dimensionless differential equations of motion are obtained. The fourth-order Runge–Kutta method is used to analyze the relationship between the connection parameters and the dynamic characteristics of the system under the two working conditions of rigid and flexible connection between the composite base and the box wall.Through the analytical analysis of vibration displacement of two-stage gear reducer and box, the theoretical basis is found for the numerical analysis results. Finally, the dynamic characteristics of the transmission system are studied by vibration resonance analysis through high and low frequency interference. It is found that in a certain frequency range, with the decrease of the mass and moment of inertia of the transmission parts corresponding to the mixed metal composite gear, the amplitude-frequency characteristic Q of the lightweight gear and gearbox transmission system is slightly lower than that of the common gear and gearbox system, and the stability of the system is increased, and the dynamic characteristics of the system are improved.