Birkhoffian symplectic algorithms derived from Hamiltonian symplectic algorithms

In this paper, we focus on the construction of structure preserving algorithms for Birkhoffian systems, based on existing symplectic schemes for the Hamiltonian equations. The key of the method is to seek an invertible transformation which drives the Birkhoffian equations reduce to the Hamiltonian equations. When there exists such a transformation,applying the corresponding inverse map to symplectic discretization of the Hamiltonian equations, then resulting difference schemes are verified to be Birkhoffian symplectic for the original Birkhoffian equations. To illustrate the operation process of the method, we construct several desirable algorithms for the linear damped oscillator and the single pendulum with linear dissipation respectively. All of them exhibit excellent numerical behavior, especially in preserving conserved quantities.     

基于辛理论的载流碳纳米管能带分析

基于连续介质力学理论和辛弹性理论,将载流碳纳米管等效为铁木辛柯梁,采用哈密顿变分原理建立了载流碳纳米管的振动控制方程;引入对偶变量将振动控制方程从拉格朗日体系导入到哈密顿体系下;通过波传播方法分析了载流碳纳米管的能带结构;研究了流体密度、流速对载流碳纳米管能带结构的影响;同时计算了载流碳纳米管的散射矩阵. 研究发现:管内流速以及流体密度对剪切频率和弯曲频率有着非常重要的影响. 研究结果表明:载流碳纳米管的剪切频率和弯曲频率因流体的加入而减小,并随流速及流体密度的增大而减小;通过对数值结果的分析发现:载流碳纳米管由于管内流体、流速以及流体密度的作用,会使得载流碳纳米管变的更“软”. 其中,哈密顿体系下所得出的载流碳纳米管弯曲频率随管内流体密度的增加而变小,有别于在拉格朗日体系下非局部梁理论所得的结论. 同时,数值结果表明散射矩阵是酉矩阵,辛体系下的入射波功率流与反射波功率流相等,即功率流守恒,体现了辛弹性力学理论的优越性.      

陀螺系统随机振动分析的辛本征展开方法

探讨了受随机载荷作用下陀螺阻尼系统随机动力响应问题.虚拟激励法作为随机振动分析的一种高效、精确方法已经广泛应用于结构抗震、抗风等工程领域.在以单类物理变量描述的Lagrange(拉格朗日)体系框架下,振型分解方法已被有效应用于上述随机振动问题的模型自由度缩减.然而,对于陀螺系统的随机振动问题,由于陀螺效应的存在,基于Rayleigh商本征值的振型分解方法受到很大限制.对此,首先给出了陀螺系统辛本征值问题的一般形式.然后对于受平稳随机载荷激励的陀螺系统(无阻尼或有阻尼)引入虚拟激励法,基于辛本征空间展开推导了系统随机振动响应功率谱的求解列式;对于仅考虑陀螺效应的保守系统(无阻尼),该求解列式可以表述为一个显式表达式.在数值算例中,应用该文提出的方法分析了平稳随机载荷作用下一类阻尼陀螺系统的随机振动响应问题,通过与其它方法进行对比,验证了该方法的精确性和有效性.     

Dynamic modeling and simulation of deploying process for space solar power satellite receiver

To reveal some dynamic properties of the deploying process for the solar power satellite via an arbitrarily large phased array (SPS-ALPHA) solar receiver, the symplectic Runge-Kuttamethod is used to simulate the simplified model with the consideration of the Rayleigh damping effect. The system containing the Rayleigh damping can be separated and transformed into the equivalent nondamping system formally to insure the application condition of the symplectic Runge-Kutta method©First, the Lagrange equation with the Rayleigh damping governing the motion of the system is derived via the variational principle. Then, with some reasonable assumptions on the relations among the damping, mass, and stiffness matrices, the Rayleigh damping system is equivalently converted into the nondamping system formally, so that the symplectic Runge-Kutta method can be used to simulate the deploying process for the solar receiver. Finally, some numerical results of the symplectic Runge-Kutta method for the dynamic properties of the solar receiver are reported. The numerical results show that the proposed simplified model is valid for the deploying process for the SPS-ALPHA solar receiver, and the symplectic Runge-Kutta method can preserve the displacement constraints of the system well with excellent long-time numerical stability.     

Self-Trapping in Discrete Nonlinear Schrdinger Equation with Next-Nearest Neighbor Interaction

The dynamical self-trapping of an excitation propagating on one-dimensional of different sizes with nextnearest neighbor (NNN) interaction is studied by means of an explicit fourth order symplectic integrator. Using localized initial conditions, the time-averaged occupation probability of the initial site is investigated which is a function of the degree of nonlinearity and the linear coupling strengths. The self-trapping transition occurs at larger values of the nonlinearity parameter as the NNN coupling strength of the lattice increases for fixed size. Furthermore, given NNN coupling strength, the self-trapping properties for different sizes are considered which are some different from the case with general nearest neighbor (NN) interaction.     

Dynamics of Charged Particles Around a Magnetized Black Hole with Pseudo-Newtonian Potential

The linear stability of equilibria of charged particles moving near a compact object with a dipole magnetic field and a pseudo-Newtonian potential is analyzed detailedly. An optimal fourth-order force gradient symplectic method, as a global symplectic integrator that can simultaneously solve both the equations of motion and the variational equations, is used to calculate fast Lyapunov indicators. In this way, dynamical structures are described, and parameter domains for causing chaos are found.     

Two-dimensional complete rational analysis of functionally graded beams within symplectic framework

Exact solutions for generally supported functionally graded plane beams are given within the framework of symplectic elasticity. The Young’s modulus is assumed to exponentially vary along the longitudinal direction while the Poisson’s ratio remains constant. The state equation with a shift-Hamiltonian operator matrix has been established in the previous work, which is limited to the Saint-Venant solution. Here, a complete rational analysis of the displacement and stress distributions in the beam is presented by exploring the eigensolutions that are usually covered up by the Saint-Venant principle. These solutions play a significant role in the local behavior of materials that is usually ignored in the conventional elasticity methods but possibly crucial to the material/structure failures. The analysis makes full use of the symplectic orthogonality of the eigensolutions. Two illustrative examples are presented to compare the displacement and stress results with those for homogenous materials, demonstrating the effects of material inhomogeneity.     

Correspondence between quantum-optical transform and classical-optical transform explored by developing Dirac’s symbolic method

By virtue of the new technique of performing integration over Dirac’s ket–bra operators, we explore quantum optical version of classical optical transformations such as optical Fresnel transform, Hankel transform, fractional Fourier transform, Wigner transform, wavelet transform and Fresnel–Hadmard combinatorial transform etc. In this way one may gain benefit for developing classical optics theory from the research in quantum optics, or vice-versa. We cannot only find some new quantum mechanical unitary operators which correspond to the known optical transformations, deriving a new theorem for calculating quantum tomogram of density operators, but also can reveal some new classical optical transformations. For examples, we find the generalized Fresnel operator (GFO) to correspond to the generalized Fresnel transform (GFT) in classical optics. We derive GFO’s normal product form and its canonical coherent state representation and find that GFO is the loyal representation of symplectic group multiplication rule. We show that GFT is just the transformation matrix element of GFO in the coordinate representation such that two successive GFTs is still a GFT. The ABCD rule of the Gaussian beam propagation is directly demonstrated in the context of quantum optics. Especially, the introduction of quantum mechanical entangled state representations opens up a new area in finding new classical optical transformations. The complex wavelet transform and the condition of mother wavelet are studied in the context of quantum optics too. Throughout our discussions, the coherent state, the entangled state representation of the two-mode squeezing operators and the technique of integration within an ordered product (IWOP) of operators are fully used. All these have confirmed Dirac’s assertion: “...for a quantum dynamic system that has a classical analogue, unitary transformation in the quantum theory is the analogue of contact transformation in the classical theory”.     

General techniques for constructing variational integrators

The numerical analysis of variational integrators relies on variational error analysis, which relates the order of accuracy of a variational integrator with the order of approximation of the exact discrete Lagrangian by a computable discrete Lagrangian. The exact discrete Lagrangian can either be characterized variationally, or in terms of Jacobi’s solution of the Hamilton-Jacobi equation. These two characterizations lead to the Galerkin and shooting constructions for discrete Lagrangians, which depend on a choice of a numerical quadrature formula, together with either a finite-dimensional function space or a one-step method. We prove that the properties of the quadrature formula, finite-dimensional function space, and underlying one-step method determine the order of accuracy and momentum-conservation properties of the associated variational integrators. We also illustrate these systematic methods for constructing variational integrators with numerical examples.     

大型不正定陀螺系统本征值问题

陀螺动力系统可以导入哈密顿辛几何体系,在哈密顿陀螺系统的辛子空间迭代法的基础上提出了一种能够有效计算大型不正定哈密顿函数的陀螺系统本征值问题的算法．利用陀螺矩阵既为哈密顿矩阵而本征值又是纯虚数或零的特点,将对应哈密顿函数为负的本征值分离开来,构造出对应哈密顿函数全为正的本征值问题,利用陀螺系统的辛子空间迭代法计算出正定哈密顿矩阵的本征值,从而解决了大型不正定陀螺系统的本征值问题,算例证明,本征解收敛得很快．