Differential Equations Of Motion Research Articles

Free vibration Analysis of Beams with Functional Gradient Materials with Cracks

The article presents the development of a dynamic model for a functionally graded material (FGM) beam incorporating cracks. Initially, it assumes an exponential distribution of material properties along the thickness direction of the beam and simulates the opening crack using a zero-mass rotational spring model. This approach enables the calculation of bending stiffness and local flexibility at the cracked section. Subsequently, drawing upon Timoshenko beam theory and von Kármán geometric nonlinear theory, the study formulates the energy equation of the beam and establishes the partial differential control equations for the cracked FGM using Hamilton's principle. The method of separation of variables is employed to discretize the partial differential motion equations into ordinary differential motion equations. Beam functions serve as mode functions, whose unknown coefficients are determined according to the boundary and continuity conditions, thereby yielding the natural frequencies and mode shapes of the cracked FGM beam. Numerical analysis is conducted to evaluate the impact of boundary conditions, the relative position of cracks, and the length-to-thickness ratio on the natural frequencies of the cracked FGM beam.

Analysis of bending vibrations of a three-layered pre-twisted sandwich beam with an exact dynamic stiffness matrix

This paper presents the novel development of the dynamic stiffness matrix for a three-layered symmetric pre-twisted sandwich beam (PTSB), aiming to investigate its free vibration characteristics. The outer layers of the beam are modeled using the Euler–Bernoulli theory, while the core is assumed to deform solely in shear. The boundary conditions and governing partial differential equations of motion are derived based on Hamilton's principle. By applying harmonic variations of displacements, the governing equations of motion are expressed as a tenth-order equation, which is solved to obtain the desired dynamic stiffness matrix. To compute the natural frequencies of in-plane and out-of-plane free vibration for both uniform and PTSBs, the Wittrick–Williams algorithm is employed. The computed frequencies are compared with the results obtained by other authors as well as those obtained from ABAQUS simulations. Various vibration modes of uniform and twisted sandwich beams are plotted and thoroughly discussed. Interestingly, contrary to straight symmetric sandwich beams, the results indicate that flexural displacements in pre-twisted symmetric sandwich beams exhibit coupling in two planes. Additionally, although there are minor changes in vibration frequencies, the mode shapes undergo significant transformations as the pre-twist angle is altered.

Noether’s Theorem of Herglotz Type for Fractional Lagrange System with Nonholonomic Constraints

This research aims to investigate the Noether symmetry and conserved quantity for the fractional Lagrange system with nonholonomic constraints, which are based on the Herglotz principle. Firstly, the fractional-order Herglotz principle is given, and the Herglotz-type fractional-order differential equations of motion for the fractional Lagrange system with nonholonomic constraints are derived. Secondly, by introducing infinitesimal generating functions of space and time, the Noether symmetry of the Herglotz type is defined, along with its criteria, and the conserved quantity of the Herglotz type is given. Finally, to demonstrate how to use this method, two examples are provided.

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.