Differential Equations Of Motion Research Articles

Dynamic modeling and analysis of viscoelastic hard-magnetic soft actuators with thermal effects

Hard-magnetic soft materials (HMSMs) are a class of magnetoactive smart polymers capable of sustaining a high residual magnetic flux density and can undergo large actuation strains under an external magnetic excitation. Due to these exceptional characteristics, soft actuators based on HMSMs hold enormous potential for remote-controlled applications. The temperature and viscoelasticity considerably influence the performance of these materials during the operation. This work aims to develop an analytical framework for modeling the dynamic behavior of a planar hard-magnetic soft actuators (HMSA) considering temperature and viscoelastic effects. The constitutive behavior of the viscoelastic HMSA is described by employing an incompressible neo-Hookean model in conjunction with a Zener rheological model and the Rayleigh dissipation function. The dynamic governing differential equations of motion are derived by utilizing the non-conservative form of the Euler–Lagrange equation. This study delves into the collective influence of temperature and viscoelastic properties on the stability, periodicity, and resonance characteristics of nonlinear vibrations exhibited by HMSM-based planar actuator subjected to dynamic magnetic loading, presenting the findings through time-history responses, Poincaré maps, and phase-plane plots. The presented results can help in the efficient and robust design of HMSM-based actuators and can also serve as an initial step toward the development of advanced actuators exposed to dynamic loading under variable temperatures for diverse applications in the fields of engineering and medicine.

SUPPRESSION OF SOMMERFELD EFFECT IN POWER TRANSMISSION SYSTEM EMPLOYING CARDAN SHAFT THROUGH PHASE ANGLE ARRANGEMENT

In a power transmission system with a cardan shaft, the Sommerfeld effect occurs, which is characterized by speed capture and release at the resonance range. Suppression of the Sommerfeld effect is critical for smooth and reliable operation. This study aims to suppress the Sommerfeld effect in a transmission system by compensating the phase angle between the two universal joints installed in the cardan shaft. The differential equations of motion representing the dynamics of the system are derived using the Lagrange equation. The responses are simulated numerically using the Runge–Kutta algorithm for scenarios with constant and gradually varying input torque. To suppress the Sommerfeld effect, the phase angle is set to 25%, 50%, 75% and 100% of the maximum twist angle observed in the subcritical speed range of the in-phase configuration. With the phase angle of 25%, the Sommerfeld effect is damped, where the output speed only deviates by 5% from the estimated value for both input torque scenarios. It is shown how the change of the phase angle attenuates the Sommerfeld effect and the system vibrations, which should be considered in the development and practical implementation.

Formation and development of vortices in power and propulsion systems

A differential equation of motion is proposed that reflects the dynamics of the vortex growth, which makes it possible to establish the boundaries of the existence of vortex flows and determine the vorticity field. Keywords:turbulence, impulse, vortex, Navier—Stokes equation, mass velocity, vorticity swgeorgy@gmail.com, toolgapodymova@yandex.ru

Semiclassical resolution of the black hole singularity inspired in the minimal uncertainty approach

We propose a new lapse function that simplifies the Hamiltonian constraint, describing the interior of the black hole in terms of the Ashtekar-Barbero variables, into a more straightforward form. The new Hamiltonian leads to different equations of motion than those found in the literature, but through a suitable transformation between temporal parameters, it is found that such a choice leads us to the classical solutions of the Schwarzschild metric, still preserving the physical singularity. In order to resolve this singularity, and inspired by the minimal uncertainty approach, we modify the classical algebra between the dynamic variables of the model, imposing an effective dynamics within the black hole. As a consequence, one of the dynamic variables, denoted by pb, acquires a minimum value at the singularity t=0, and on the other hand, the variable related to the radius of the 2-sphere, pc, leads to the resolution of the classical singularity of the black hole by replacing it with a bounce that connects the interior of the black hole with the interior of the white hole. This bounce occurs in the Planck-scale region, where a new event horizon manifests. Upon crossing this horizon, the nature of the interval changes from spatial to temporal outside the white hole.

Free flexural vibrations of standard wide-flange H-beams with consideration of the shear effect

The subject of the paper is a simply supported standard wide-flange H-beam. Cross sections of this beam is analytically described as a three-layer structure. The shear effect in its succesive layers is taking into account with consideration of the classical shear stress formula called Zhuravsky shear stress. Based on Hamilton’s principle, two differential equations of motion are obtained. These equations are analytically solved and the fundamental natural frequency of flexural vibration for this beam is derived. Examplary calculations are carried out for selected five I-beams.

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.

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.

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.

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

AbstractAn 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

AbstractA 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.