Structure-preserving properties of three differential schemes for oscillator system

A numerical method for the Hamiltonian system is required to preserve some structure-preserving properties. The current structure-preserving method satisfies the requirements that a symplectic method can preserve the symplectic structure of a finite dimension Hamiltonian system, and a multi-symplectic method can preserve the multi-symplectic structure of an infinite dimension Hamiltonian system. In this paper,the structure-preserving properties of three differential schemes for an oscillator system are investigated in detail. Both the theoretical results and the numerical results show that the results obtained by the standard forward Euler scheme lost all the three geometric properties of the oscillator system, i.e., periodicity, boundedness, and total energy,the symplectic scheme can preserve the first two geometric properties of the oscillator system, and the St¨ormer-Verlet scheme can preserve the three geometric properties of the oscillator system well. In addition, the relative errors for the Hamiltonian function of the symplectic scheme increase with the increase in the step length, suggesting that the symplectic scheme possesses good structure-preserving properties only if the step length is small enough.     

STRESS SINGULARITY ANALYSIS OF MULTI-MATERIAL WEDGES UNDER ANTIPLANE DEFORMATION

With the help of the coordinate transformation technique, the symplectic dual solving system is developed for multi-material wedges under antiplane deformation. A virtue of present method is that the compatibility conditions at interfaces of a multi-material wedge are expressed directly by the dual variables, therefore the governing equation of eigenvalue can be derived easily even with the increase of the material number. Then, stress singularity on multi-material wedges under antiplane deformation is investigated, and some solutions can be presented to show the validity of the method. Simultaneously, an interesting phenomenon is found and proved strictly that one of the singularities of a special five-material wedge is independent of the crack direction.     

基于哈密顿动力系统新变分原理的保辛算法之一：变分原理和算法构造

提出了哈密顿动力系统的一个新变分原理,并基于此变分原理构造了四类保辛算法。通过新的变分原理定义修正作用量,然后将位移和动量采用拉格朗日多项式近似,并采用高斯积分对时间近似积分得到近似的修正作用量。在修正作用量的基础上,通过选择时间步两端不同的位移或动量作为独立变量,可构造四种不同类型的保辛算法。     

连续系统的界带分析方法

基于分析结构力学提出的界带分析方法,将子结构间的分界面延拓为有一定宽度的分界带/分界域,从而可以用于分析计算结构的非局部效应。界带分析方法首先在离散结构的分析计算中取得了成功,从而验证了该套理论算法的准确性。离散结构按界带宽度（影响域范围）划分子结构,因而限制了子结构区段积分计算的最小步长;而连续系统则要求可以实现任意步长的积分运算。通过引入步进的计算方法,界带分析方法可以实现任意步长的积分计算,进而可以解决连续系统的积分问题。通过数值算例验证了连续系统的界带分析方法的准确性和可行性,也为进一步研究该套计算方法在分析动力学中的应用打下基础。     

Estimation of goal functional error arising from iterative solution of Euler equations

Estimation of the error arising in the cost (goal) functional due to stopping the iterative process is considered for a steady problem solved by temporal relaxation. The functional error is calculated using an iteration residual along with related adjoint parameters. Numerical tests demonstrate the applicability of this approach for the steady 2D Euler equations.     

Hamilton体系下功能梯度梁的热冲击动力屈曲分析

在Hamilton体系下,基于Euler梁理论研究了功能梯度材料梁受热冲击载荷作用时的动力屈曲问题;将非均匀功能梯度复合材料的物性参数假设为厚度坐标的幂函数形式,采用Laplace变换法和幂级数法解析求得热冲击下功能梯度梁内的动态温度场:首先将功能梯度梁的屈曲问题归结为辛空间中系统的零本征值问题,梁的屈曲载荷与屈曲模态分别对应于Hamilton体系下的辛本征值和本征解问题,由分叉条件求得屈曲模态和屈曲热轴力,根据屈曲热轴力求解临界屈曲升温载荷。给出了热冲击载荷作用下一类非均匀梯度材料梁屈曲特性的辛方法研究过程,讨论了材料的梯度特性、结构几何参数和热冲击载荷参数对临界温度的影响。     

Hamiltonian Systems on Time Scales

Linear and nonlinear Hamiltonian systems are studied on time scales . We unify symplectic flow properties of discrete and continuous Hamiltonian systems. A chain rule which unifies discrete and continuous settings is presented for our so-called alpha derivatives on generalized time scales. This chain rule allows transformation of linear Hamiltonian systems on time scales under simultaneous change of independent and dependent variables, thus extending the change of dependent variables recently obtained by Do lý and Hilscher. We also give the Legendre transformation for nonlinear Euler–Lagrange equations on time scales to Hamiltonian systems on time scales.     

Non-Parallel Acoustic Receptivity of a Blasius Boundary Layer Using an Adjoint Approach

Localized and non-localized acoustic receptivity for a Blasius boundary layer is investigated using the adjoint Parabolized Stability Equations. The scattering of an acoustic wave onto a hump, a rectangular roughness or a wall steady blowing and suction is described. Comparisons with local approaches, triple deck theory, direct numerical simulations and experiments are successfully shown. Non-parallel effects are discussed. For comparable parameters, the non-localized receptivity problem produces amplitudes one order of magnitude larger than for the case of localized receptivity. This revised version was published online in July 2006 with corrections to the Cover Date.     

SYMPLECTIC SCHEMES OF VORTEX SYSTEM IN HALF PLAIN

l introductionThe vortex system is a system of POint voids, and a model Of incompressible inviscidnow inspired by the idea of an almOSt POtence flOW. The voracity in the now is concentratedin N-vortices (i. e., POints at which the vortloty field is singUlar) [4]. An ideal incompressible now can be approximated by the motion Of a ~ system which is not only a usefulheuristic tool in the analysis of the general propelles of solutionS Of Euler equations, but also auseful stachg POint for th…     

From the restricted to the full three-body problem

The three-body problem with all the classical integrals fixed and all the symmetries removed is called the reduced three-body problem. We use the methods of symplectic scaling and reduction to show that the reduced planar or spatial three-body problem with one small mass is to the first approximation the product of the restricted three-body problem and a harmonic oscillator. This allows us to prove that many of the known results for the restricted problem have generalizations for the reduced three-body problem.

For example, all the non-degenerate periodic solutions, generic bifurcations, Hamiltonian-Hopf bifurcations, bridges and natural centers known to exist in the restricted problem can be continued into the reduced three-body problem. The classic normalization calculations of Deprit and Deprit-Bartholomé show that there are two-dimensional KAM invariant tori near the Lagrange point in the restricted problem. With the above result this proves that there are three-dimensional KAM invariant tori near the Lagrange point in the reduced three-body problem.