Data-driven modeling of a four-dimensional stochastic projectile system

The dynamical modeling of projectile systems with sufficient accuracy is of great difficulty due to high-dimensional space and various perturbations. With the rapid development of data science and scientific tools of measurement recently, there are numerous data-driven methods devoted to discovering governing laws from data. In this work, a data-driven method is employed to perform the modeling of the projectile based on the Kramers-Moyal formulas. More specifically, the four-dimensional projectile system is assumed as an Itô stochastic differential equation. Then the least square method and sparse learning are applied to identify the drift coefficient and diffusion matrix from sample path data, which agree well with the real system. The effectiveness of the data-driven method demonstrates that it will become a powerful tool in extracting governing equations and predicting complex dynamical behaviors of the projectile.     

Energy diffusion controlled reaction rate of reacting particle driven by broad-band noise

The energy diffusion controlled reaction rate of a reacting particle with linear weak damping and broad-band noise excitation is studied by using the stochastic averaging method. First, the stochastic averaging method for strongly nonlinear oscillators under broad-band noise excitation using generalized harmonic functions is briefly introduced. Then, the reaction rate of the classical Kramers' reacting model with linear weak damping and broad-band noise excitation is investigated by using the stochastic averaging method. The averaged It? stochastic differential equation describing the energy diffusion and the Pontryagin equation governing the mean first-passage time (MFPT) are established. The energy diffusion controlled reaction rate is obtained as the inverse of the MFPT by solving the Pontryagin equation. The results of two special cases of broad-band noises, i.e. the harmonic noise and the exponentially corrected noise, are discussed in details. It is demonstrated that the general expression of reaction rate derived by the authors can be reduced to the classical ones via linear approximation and high potential barrier approximation. The good agreement with the results of the Monte Carlo simulation verifies that the reaction rate can be well predicted using the stochastic averaging method.     

Coupled magneto-electro-mechanical lumped parameter model for a novel vibration-based magneto-electro-elastic energy harvesting systems

The objective of this paper is to present a coupled magneto-electro-mechanical (MEM) lumped parameter model for the response of the proposed magneto-electro-elastic (MEE) energy harvesting systems under base excitation. The proposed model can be used to create self-powering systems, which are not limited to a finite battery energy. As a novel approach, the MEE composites are used instead of the conventional piezoelectric materials in order to enhance the harvested electrical power. The considered structure consists of a MEE layer deposited on a layer of non-MEE material, in the framework of unimorph cantilever bars (longitudinal displacement) and beams (transverse displacement). To use the generated electrical potential, two electrodes are connected to the top and bottom surfaces of the MEE layer. Additionally, a stationary external coil is wrapped around the vibrating structure to induce a voltage in the coil by the magnetic field generated in the MEE layer. In order to simplify the design procedure of the proposed energy harvester and obtain closed form solutions, a lumped parameter model is prepared. As a first step in modeling process, the governing constitutive equations, Gauss's and Faraday's laws, are used to derive the coupled MEM differential equations. The derived equations are then solved analytically to obtain the dynamic behavior and the harvested voltages and powers of the proposed energy harvesting systems. Finally, the influences of the parameters that affect the performance of the MEE energy harvesters such as excitation frequency, external resistive loads and number of coil turns are discussed in detail. The results clearly show the benefit of the coil circuit implementation, whereby significant increases in the total useful harvested power as much as 38% and 36% are obtained for the beam and bar systems, respectively.     

M-shaped asymmetric nonlinear oscillator for broadband vibration energy harvesting: Harmonic balance analysis and experimental validation

Over the past few years, nonlinear oscillators have been given growing attention due to their ability to enhance the performance of energy harvesting devices by increasing the frequency bandwidth. Duffing oscillators are a type of nonlinear oscillator characterized by a symmetric hardening or softening cubic restoring force. In order to realize the cubic nonlinearity in a cantilever at reasonable excitation levels, often an external magnetic field or mechanical load is imposed, since the inherent geometric nonlinearity would otherwise require impractically high excitation levels to be pronounced. As an alternative to magnetoelastic structures and other complex forms of symmetric Duffing oscillators, an M-shaped nonlinear bent beam with clamped end conditions is presented and investigated for bandwidth enhancement under base excitation. The proposed M-shaped oscillator made of spring steel is very easy to fabricate as it does not require extra discrete components to assemble, and furthermore, its asymmetric nonlinear behavior can be pronounced yielding broadband behavior under low excitation levels. For a prototype configuration, linear and nonlinear system parameters extracted from experiments are used to develop a lumped-parameter mathematical model. Quadratic damping is included in the model to account for nonlinear dissipative effects. A multi-term harmonic balance solution is obtained to study the effects of higher harmonics and a constant term. A single-term closed-form frequency response equation is also extracted and compared with the multi-term harmonic balance solution. It is observed that the single-term solution overestimates the frequency of upper saddle-node bifurcation point and underestimates the response magnitude in the large response branch. Multi-term solutions can be as accurate as time-domain solutions, with the advantage of significantly reduced computation time. Overall, substantial bandwidth enhancement with increasing base excitation is validated experimentally, analytically, and numerically. As compared to the 3 dB bandwidth of the corresponding linear system with the same linear damping ratio, the M-shaped oscillator offers 3200, 5600, and 8900 percent bandwidth enhancement at the root-mean-square base excitation levels of 0.03g, 0.05g, and 0.07g, respectively. The M-shaped configuration can easily be exploited in piezoelectric and electromagnetic energy harvesting as well as their hybrid combinations due to the existence of both large strain and kinetic energy regions. A demonstrative case study is given for electromagnetic energy harvesting, revealing the importance of higher harmonics and the need for multi-term harmonic balance analysis for predicting the electrical power output accurately.     

Renormalized diffusion coefficients for electron-hole plasmas in crossed fields

The diffusion and drift of an excess plasma in a semiconductor is described with magnetohydrodynamic two-fluid equations including the fluctuating electric field produced by the equilibrium plasma in the sample. Using the weak coupling limit an equation of motion for the mean density of the excess plasma is established with renormalized drift and diffusion coefficients. With the aid of the fluctuation dissipation theorem these coefficients are expressed in terms of the dielectric function and discussed in detail for stable systems. The renormalized diffusion coefficient differs from the bare one by an additional term with thet ^–3/2-long time dependence. It is shown that this term in addition represents an anomalous diffusion rate proportionalB ^–1 which overweights the classical ambipolar diffusion for sufficiently strong fields, but decreases with increasing external electric field. The results are compared with experimental data.     

Microscopic transport theory of heavy-ion collisions

Drift and diffusion coefficients for the variables of mass fragmentation and excitation energy are studied for deeply inelastic collisions. The transport coefficients are obtained in closed form as functions of the parameters of the interaction matrix elements between nucleonic states and as functions of the binding energy of the intermediate rotating quasimolecular configuration. Drift and diffusion coefficients for excitation and mass transfer are related (dissipation fluctuation theorem). In good approximation these relations take the simple form of Einstein's relation between the mobility and the diffusion coefficient of a particle in a medium. For the total mass numbers 100, 250 and 500 the results are discussed in detail. The transport coefficients are compared with experimental results. Within the uncertainties of determination from experimental data, the drift and diffusion coefficients are well described by two previously adjusted parameters of the mean interaction matrix elements. Consequences for the production of superheavy elements in deeply inelastic collisions of U on U or Cf are discussed.     

Investigations of Viscous Dissipation in Stagnation Point Flow Past a Stretchable Riga Wall: Modern Analysis of Heat Transport

This article scrutinizes the features of viscous dissipation in the stagnation point flow past through a linearly stretched Riga wall by implementing Cattaneo-Christov heat flux model. Viscous dissipation is carried out in Cattaneo-Christov diffusion analysis for the first time in this letter. As a result of Cattaneo-Christov model, some extra terms of viscous dissipation are appeared in the energy equation. These extra terms of viscous dissipation are missing in the literature. On the utilization of suitable transformations, the equations governing the problem are reduced under the boundary layer approximation into the non-linear and dimensionless ordinary differential equations. Convergent approach is utilized to solve the dimensionless governing equations. The solution thus acquired is used to highlight the effects of emerging parameters on velocity distribution and fluid's temperature through the graphs. Features of the drag force (or skin friction co-efficient) are graphically interpreted. It is noticed that the presence of modified Hartman number helps to reduce the fluid's temperature but enhances the velocity profile. Further an enlargement in the value of thermal time relaxation parameter helps to decrease the temperature distribution.     

Equivalent damping and frequency change for linear and nonlinear hybrid vibrational energy harvesting systems

A unified approximation method is derived to illustrate the effect of electro-mechanical coupling on vibration-based energy harvesting systems caused by variations in damping ratio and excitation frequency of the mechanical subsystem. Vibrational energy harvesters are electro-mechanical systems that generate power from the ambient oscillations. Typically vibration-based energy harvesters employ a mechanical subsystem tuned to resonate with ambient oscillations. The piezoelectric or electromagnetic coupling mechanisms utilized in energy harvesters, transfers some energy from the mechanical subsystem and converts it to an electric energy. Recently the focus of energy harvesting community has shifted toward nonlinear energy harvesters that are less sensitive to the frequency of ambient vibrations. We consider the general class of hybrid energy harvesters that use both piezoelectric and electromagnetic energy harvesting mechanisms. Through using perturbation methods for low amplitude oscillations and numerical integration for large amplitude vibrations we establish a unified approximation method for linear, softly nonlinear, and bi-stable nonlinear energy harvesters. The method quantifies equivalent changes in damping and excitation frequency of the mechanical subsystem that resembles the backward coupling from energy harvesting. We investigate a novel nonlinear hybrid energy harvester as a case study of the proposed method. The approximation method is accurate, provides an intuitive explanation for backward coupling effects and in some cases reduces the computational efforts by an order of magnitude.     

Collective Excitations and Nonlinear Dynamics of 1D BEC with Two- and Three-body Interactions in Anharmonic Traps

The collective excitations of low-dimensional Bose-Einstein condensates with two- and three-body interactions in anharmonic potentials are investigated. Using the standard variational approach, the governing equations of motions for the low-energy excitations are obtained by solving time-dependent Gross-Pitaevskii-Ginzburg equation, and the excitation spectrums are calculated in small amplitude limit. The frequency shift and nonlinear mode coupling induced by the anharmonic distortion (adding cubic, quartic, or quintic term to a harmonic trap) are studied.     

扰动扩散方程初值问题的近似对称约化

本文利用近似广义条件对称方法研究一类带有源项的非线性扩散方程的初值问题.给出所研究方程的分类并将偏微分方程的初值问题约化为常微分方程的初值问题,通过求解约化后的常微分方程组可得相对应偏微分方程初值问题的近似解.     

带有分数阶阻尼的压电能量采集系统相干共振

研究了含分数阶阻尼的双稳态能量采集系统的相干共振. 建立了带有分数阶阻尼的轴向受压梁压电能量采集系统动力学模型. 对于分数阶方程, 采用Euler-Maruyama-Leipnik方法进行求解, 计算了不同阻尼阶数下的能量采集系统的信噪比、响应均值、跃迁数目等统计物理量. 结果表明: 此压电能量采集系统在随机激励下可以实现相干共振, 阻尼阶数对相干共振的临界噪声强度和相干共振幅值有很大影响. 关键词： 分数阶阻尼 随机激励 能量采集系统 相干共振     

Coupled out of plane vibrations of spiral beams for micro-scale applications

An analytical method is proposed to calculate the natural frequencies and the corresponding mode shape functions of an Archimedean spiral beam. The deflection of the beam is due to both bending and torsion, which makes the problem coupled in nature. The governing partial differential equations and the boundary conditions are derived using Hamilton’s principle. Two factors make the vibrations of spirals different from oscillations of constant radius arcs. The first is the presence of terms with derivatives of the radius in the governing equations of spirals and the second is the fact that variations of radius of the beam causes the coefficients of the differential equations to be variable. It is demonstrated, using perturbation techniques that the derivative of the radius terms have negligible effect on structure’s dynamics. The spiral is then approximated with many merging constant-radius curved sections joined together to approximate the slow change of radius along the spiral. The equations of motion are formulated in non-dimensional form and the effect of all the key parameters on natural frequencies is presented. Non-dimensional curves are used to summarize the results for clarity. We also solve the governing equations using Rayleigh’s approximate method. The fundamental frequency results of the exact and Rayleigh’s method are in close agreement. This to some extent verifies the exact solutions. The results show that the vibration of spirals is mostly torsional which complicates using the spiral beam as a host for a sensor or energy harvesting device.     

Dynamic stability of a slender beam under horizontal–vertical excitations

The dynamic stability of a vertically standing cantilevered beam simultaneously excited in both horizontal and vertical directions at its base is studied theoretically. The beam is assumed to be an inextensible Euler–Bernoulli beam. The governing equation of motion is derived using Hamilton's principle and has a nonlinear elastic term and a nonlinear inertia term. A forced horizontal external term is added to the parametrically excited system. Applying Galerkin's method for the first bending mode, the forced Mathieu–Duffing equation is derived. The frequency response is obtained by the harmonic balance method, and its stability is investigated using the phase plane method. Excitation frequencies in the horizontal and vertical directions are taken as 1:2, from which we can investigate the influence of the forced response under horizontal excitation on the parametric instability region under vertical excitation. Three criteria for the instability boundary are proposed. The influences of nonlinearities and damping of the beam on the frequency response and parametric instability region are also investigated. The present analytical results for instability boundaries are compared with those of experiments carried out by one of the authors.     

Vibration confinement and energy harvesting in flexible structures using collocated absorbers and piezoelectric devices

We propose an optimal design for supplementing flexible structures with a set of absorbers and piezoelectric devices for vibration confinement and energy harvesting. We assume that the original structure is sensitive to vibrations and that the absorbers are the elements where the vibration energy is confined and then harvested by means of piezoelectric devices. The design of the additional mechanical and electrical components is formulated as a dynamic optimization problem in which the objective function is the total energy of the uncontrolled structure. The locations, masses, stiffnesses, and damping coefficients of these absorbers and capacitances, load resistances, and electromechanical coupling coefficients are optimized to minimize the total energy of the structure. We use the Galerkin procedure to discretize the equations of motion that describe the coupled dynamics of the flexible structure and the added absorbers and harvesting devices. We develop a numerical code that determines the unknown parameters of a pre-specified set of absorbers and harvesting components. We input a set of initial values for these parameters, and the code updates them while minimizing the total energy in the uncontrolled structure. To illustrate the proposed design, we consider a simply supported beam with harmonic external excitations. Here, we consider two possible configurations for each of the additional piezoelectric devices, either embedded between the structure and the absorbers or between the ground and absorbers. We present simulations of the harvested power and associated voltage for each pair of collocated absorber and piezoelectric device. The simulated responses of the beam show that its energy is confined and harvested simultaneously.     

Nonlinear diffusion under a time dependent external force: q-maximum entropy solutions

Nonlinear diffusion equations provide useful models for a number of interesting phenomena, such as diffusion processes in porous media. We study here a family of nonlinear Fokker-Planck equations endowed both with a power-law nonlinear diffusion term and a drift term with a time dependent force linear in the spatial variable. We show that these partial differential equations exhibit exact time dependent particular solutions of the Tsallis maximum entropy (q-MaxEnt) form. These results constitute generalizations of previous ones recently discussed in the literature [C. Tsallis, D.J. Bukman, Phys. Rev. E 54, R2197 (1996)], concerning q-MaxEnt solutions to nonlinear Fokker-Planck equations with linear, time independent drift forces. We also show that the present formalism can be used to generate approximate q-MaxEnt solutions for nonlinear Fokker-Planck equations with time independent drift forces characterized by a general spatial dependence. Received 25 April 2001 and Received in final form 6 June 2001     

Elasticity and dislocations in quasicrystals with 18-fold symmetry

Elasticity problems of quasicrystals with 18-fold rotational symmetry are studied. Constitutive equations and governing equations are obtained. For static elasticity problems, the displacement vectors in two phason fields are expressed in terms of two pairs of associated harmonic functions or two analytic functions. For dynamic problems, the displacement vectors can be represented in terms of an auxiliary function satisfying a fourth-order partial differential equation. A general solution of phasons is given by the solution of two diffusion equations. Phason elastic fields induced by a dislocation in a quasicrystal with 18-fold symmetry are determined and exhibit an inverse singularity.     

时变偏微分方程的贝叶斯稀疏识别方法

在数据驱动的建模中,通过测量或模拟得到时空数据,我们发现基于拉普拉斯先验的贝叶斯稀疏识别方法能有效地恢复时变偏微分方程的稀疏系数.本文将贝叶斯稀疏识别方法运用于各种时变偏微分方程模型(KdV方程、Burgers方程、Kuramoto-Sivashinsky方程、反应-扩散方程、非线性薛定谔方程和纳维-斯托克斯方程)的方...     

轴向加速运动黏弹性梁受迫振动中的混沌动力学

研究了黏弹性轴向运动梁在外部激励和参数激励共同作用下横向振动的混沌非线性动力学行为. 引入有限支撑刚度, 并考虑黏弹性本构关系取物质导数, 同时计入由梁轴向加速度引起的沿径向变化的轴力, 建立轴向运动黏弹性梁横向非线性振动的偏微分-积分模型. 通过Galerkin截断方法研究了外部激励的频率和因速度简谐脉动引起的参数激励的频率在不可通约关系时轴向运动连续体的非线性动力学行为, 并对不同截断阶数的数值预测进行了对比. 基于对控制方程的Galerkin截断, 得到离散化的常微分方程组, 使用四阶Runge-Kutta方法求解. 基于此数值解, 运用非线性动力学时间序列分析方法, 通过Poincaré 映射, 观察到轴向运动梁随扰动速度幅值的倍周期分岔现象, 并比较了有无外部激励对倍周期分岔的影响. 分别在低速以及近临界高速运动状态下, 从相平面图、Poincaré 映射以及频谱分析的角度识别了系统中存在的准周期运动形态. 关键词： 轴向运动梁 非线性 混沌 分岔     

外电场作用下C6H6的分子结构及其特性

采用 B3P86方法计算了C6H6在外电场作用下分子系统的能量及分子轨道能级、谐振频率、CH键离解能和极化率等分子特性;通过跃迁矩阵元计算了各激发态的激发能,吸收及自发辐射系数,从而进一步分析了苯分子的外场效应.对于高分子化合物尤其是具有苯环结构的有机高分子化合物的外场作用特性的研究具有参考价值.