Symplectic solution system for reissner plate bending

Based on the Hellinger-Reissner variatonal principle for Reissner plate bendingand introducing dual variables, Hamiltonian dual equations for Reissner plate bending werepresented. Therefore Hamiltonian solution system can also be applied to Reissner platebending problem, and the transformation from Euclidian space to symplectic space and fromLagrangian system to Hamiltonian system was realized. So in the symplectic space whichconsists of the original variables and their dual variables, the problem can be solved viaeffective mathematical physics methods such as the method of separation of variables andeigenfunction-vector expansion. All the eigensolutions and Jordan canonical formeigensolutions for zero eigenvalue of the Hamiltonian operator matrix are solved in detail, and their physical meanings are showed clearly. The adjoint symplectic orthonormal relation of the eigenfunction vectors for zero eigenvalue are formed. It is showed that the alleigensolutions for zero eigenvalue are basic solutions of the Saint-Venant problem and theyform a perfect symplectic subspace for zero eigenvalue. And the eigensolutions for nonzeroeigenvalue are covered by the Saint-Venant theorem. The symplectic solution method is notthe same as the classical semi-inverse method and breaks through the limit of the traditional semi-inverse solution. The symplectic solution method will have vast application.     

多层层合板圣维南问题的解析解

将哈密尔顿体系理论引入到多层层合板问题之中，建立了一套求解该问题的横向哈密尔顿算子矩阵的本征函数向量展开解法，并成功地求解出圣维南问题的解析解．进一步显示了弹性力学新求解体系的有效性及其应用潜力     

A Whitham–Boussinesq long-wave model for variable topography

We study the propagation of water waves in a channel of variable depth using the long-wave asymptotic regime. We use the Hamiltonian formulation of the problem in which the non-local Dirichlet–Neumann operator appears explicitly in the Hamiltonian, and propose a Hamiltonian model for bidirectional wave propagation in shallow water that involves pseudo-differential operators that simplify the variable-depth Dirichlet–Neumann operator. The model generalizes the Boussinesq system, as it includes the exact dispersion relation in the case of constant depth. Analogous models were proposed by Whitham for unidirectional wave propagation. We first present results for the normal modes and eigenfrequencies of the linearized problem. We see that variable depth introduces effects such as a steepening of the normal modes with the increase in depth variation, and a modulation of the normal mode amplitude. Numerical integration also suggests that the constant depth nonlocal Boussinesq model can capture qualitative features of the evolution obtained with higher order approximations of the Dirichlet–Neumann operator. In the case of variable depth we observe that wave-crests have variable speeds that depend on the depth. We also study the evolutions of Stokes waves initial conditions and observe certain oscillations in width of the crest and also some interesting textures and details in the evolution of wave-crests during the passage over obstacles.     

Symplectic eigenfunction expansion theorem for elasticity of rectangular planes with two simply-supported opposite sides

The eigenvalue problem of the Hamiltonian operator associated with plane elasticity problems is investigated.The eigenfunctions of the operator are directly solved with mixed boundary conditions for the displacement and stress in a rectangular region.The completeness of the eigenfunctions is then proved,providing the feasibility of using separation of variables to solve the problems.A general solution is obtained with the symplectic eigenfunction expansion theorem.     

Paradox solution on elastic wedge dissimilar materials

According to the Hellinger-Reissner variational principle and introducing proper transformation of variables , the problem on elastic wedge dissimilar materials can be led to Hamiltonian system, so the solution of the problem can be got by employing the separation of variables method and symplectic eigenfunction expansion under symplectic space, which consists of original variables and their dual variables . The eigenvalue - 1 is a special one of all symplectic eigenvalue for Hamiltonian system in polar coordinate . In general, the eigenvalue - 1 is a single eigenvalue, and the classical solution of an elastic wedge dissimilar materials subjected to a unit concentrated couple at the vertex is got directly by solving the eigenfunction vector for eigenvalue - 1. But the eigenvalue - 1 becomes a double eigenvalue when the vertex angles and modulus of the materials satisfy certain definite relationships and the classical solution for the stress distribution becomes infinite at this moment, that is, the para     

矩形中厚板Hamilton正则方程的解析解

运用Fourier分析方法,建立了对边简支的矩形中厚板弯曲问题的完备的辛本征展开. 借助于Mathematica软件的帮助,得到了来源于矩形中厚板问题的Hamilton算子的本征函数. 接着证明了本征函数系的完备性,这为使用分离变量法求解相应问题提供了理论保证;进而运用完备性定理,得到了问题的解析解;一个数值算例验证了结果的正确性.     

哈密顿体系与弹性楔体问题

将哈密体系引入到级坐标下的弹性力学楔体问题,利用该体系辛空间的性质,将问题化为本征值和本征向量求解上,得到了完备的解空间,从而改变了弹性力学传统的拉格朗日体系以应力函数为特征的半逆法的讨论去解决该类问题的思路,给出了一条求解该类问题的直接法。     

Closed form stress distribution in 2D elasticity for all boundary conditions

This paper applies a Hamiltonian method to study analytically the stress dis- tributions of orthotropic two-dimensional elasticity in(x,z)plane for arbitrary boundary conditions without beam assumptions.It is a method of separable variables for partial differential equations using displacements and their conjugate stresses as unknowns.Since coordinates(x,z)can not be easily separated,an alternative symplectic expansion is used. Similar to the Hamiltonian formulation in classical dynamics,we treat the x coordinate as time variable so that z becomes the only independent coordinate in the Hamiltonian ma- trix differential operator.The exponential of the Hamiltonian matrix is symplectic.There are homogenous solutions with constants to be determined by the boundary conditions and particular integrals satisfying the loading conditions.The homogenous solutions consist of the eigen-solutions of the derogatory zero eigenvalues(zero eigen-solutions) and that of the well-behaved nonzero eigenvalues(nonzero eigen-solutions).The Jordan chains at zero eigenvalues give the classical Saint-Venant solutions associated with aver- aged global behaviors such as rigid-body translation,rigid-body rotation or bending.On the other hand,the nonzero eigen-solutions describe the exponentially decaying localized solutions usually ignored by Saint-Venant's principle.Completed numerical examples are newly given to compare with established results.     

基于混合Lie算子辛算法的不变流形计算

平动点附近周期轨道的不变流形因其在低能轨道转移中起着重要作用而受到广泛关注.在设计低能轨道过程中不变流形要实时进行能量匹配,但利用传统数值积分方法进行积分时能量会耗散.显式辛算法具有比隐式辛算法计算效率高的优势,但其要求Hamilton系统必须分成两个可积的部分,而旋转坐标系下的圆型限制性三体问题是不可分的,因而显式辛算法难以用于求解旋转坐标系下的圆型限制性三体问题.本文通过引入混合Lie算子,成功实现了带三阶导数项的力梯度辛算法对圆型限制性三体问题的求解,并将基于混合Lie算子的带三阶导数项的辛算法与Runge-Kutta78算法和Runge-Kutta45算法进行仿真对比,仿真结果表明基于混合Lie算子的含有三阶导数项的辛算法位置精度高、能量误差小且计算效率高.利用基于混合Lie算子的带三阶导数项的辛算法计算不变流形,可以实现低能轨道转移过程中轨道拼接点的能量精准匹配.     

含裂纹纳米板振动问题的哈密顿体系方法

基于裂纹处范德华力效应,采用非局部弹性理论构造纳米板模型,并通过导入哈密顿体系建立含裂纹纳米板振动问题的对偶正则控制方程组。在全状态向量表示的哈密顿体系下,将含裂纹纳米板的固有频率和振型问题归结为广义辛本征值和本征解问题。利用哈密顿体系具有的辛共轭正交关系,得到问题解的级数解析表达式。结合边界条件,得到固有频率与辛本征值的代数方程关系式,进而直接给出固有频率的表达式。数值结果表明,非局部尺寸参数和裂纹长度对纳米板振动的各阶固有频率有直接的影响。对比表明,辛方法是准确且可靠的,可为工程应用提供依据。     

HAMILTONIAN SYSTEM AND THE SAINT VENANT PROBLEM IN ELASTICITY

ＨＡＭＩＬＴＯＮＩＡＮＳＹＳＴＥＭＡＮＤＴＨＥＳＡＩＮＴＶＥＮＡＮＴＰＲＯＢＬＥＭＩＮＥＬＡＳＴＩＣＩＴＹＺｈｏｎｇＷａｎｘｉｅ（钟万勰）；ＸｕＸｉｎｓｈｅｎｇ（徐新生）；ＺｈａｎｇＨｏｎｇｗｕ（张洪武）（ＲｅｃｅｉｖｅｄＪｕｎｅ５，１９９５）Ａｂｓ...     

EIGENVALUE PROBLEM OF A LARGE SCALE INDEFINITE GYROSCOPIC DYNAMIC SYSTEM

Gyroscopic dynamic system can be introduced to Hamiltonian system.Based on an adjoint symplectic subspace iteration method of Hamiltonian gyroscopic system, an adjoint symplectic subspace iteration method of indefinite Hamiltonian function gy- roscopic system was proposed to solve the eigenvalue problem of indefinite Hamiltonian function gyroscopic system.The character that the eigenvalues of Hamiltonian gyroscopic system are only pure imaginary or zero was used.The eigenvalues that Hamiltonian function is negative can be separated so that the eigenvalue problem of positive definite Hamiltonian function system was presented,and an adjoint symplectic subspace iteration method of positive definite Hamiltonian function system was used to solve the separated eigenvalue problem.Therefore,the eigenvalue problem of indefinite Hamiltonian function gyroscopic system was solved,and two numerical examples were given to demonstrate that the eigensolutions converge exactly.     

A Dimension-Breaking Phenomenon for Water Waves with Weak Surface Tension

It is well known that the water-wave problem with weak surface tension has small-amplitude line solitary-wave solutions which to leading order are described by the nonlinear Schrödinger equation. The present paper contains an existence theory for three-dimensional periodically modulated solitary-wave solutions which have a solitary-wave profile in the direction of propagation and are periodic in the transverse direction; they emanate from the line solitary waves in a dimension-breaking bifurcation. In addition, it is shown that the line solitary waves are linearly unstable to long-wavelength transverse perturbations. The key to these results is a formulation of the water wave problem as an evolutionary system in which the transverse horizontal variable plays the role of time, a careful study of the purely imaginary spectrum of the operator obtained by linearising the evolutionary system at a line solitary wave, and an application of an infinite-dimensional version of the classical Lyapunov centre theorem.     

Mixed-mode thermal stress intensity factors from the finite element discretized symplectic method

A finite element discretized symplectic method is introduced to find the thermal stress intensity factors (TSIFs) under steady-state thermal loading by symplectic expansion. The cracked body is modeled by the conventional finite elements and divided into two regions: near and far fields. In the near field, Hamiltonian systems are established for the heat conduction and thermoelasticity problems respectively. Closed form temperature and displacement functions are expressed by symplectic eigen-solutions in polar coordinates. Combined with the analytic symplectic series and the classical finite elements for arbitrary boundary conditions, the main unknowns are no longer the nodal temperature and displacements but are the coefficients of the symplectic series after matrix transformation. The TSIFs, temperatures, displacements and stresses at the singular region are obtained simultaneously without any post-processing. A number of numerical examples as well as convergence studies are given and are found to be in good agreement with the existing solutions.     

Analytical and numerical methods of symplectic system for Stokes flow in two-dimensional rectangular domain

In this paper, a new analytical method of symplectic system, Hamiltonian system, is introduced for solving the problem of the Stokes flow in a two-dimensional rectangular domain. In the system, the fundamental problem is reduced to an eigenvalue and eigensolution problem. The solution and boundary conditions can be expanded by eigensolutions using adjoint relationships of the symplectic ortho-normalization between the eigensolutions. A closed method of the symplectic eigensolution is presented based on completeness of the symplectic eigensolution space. The results show that fundamental flows can be described by zero eigenvalue eigensolutions, and local effects by nonzero eigenvalue eigensolutions. Numerical examples give various flows in a rectangular domain and show effectiveness of the method for solving a variety of problems. Meanwhile, the method can be used in solving other problems.     

粘弹性厚壁筒问题的辛本征解方法

将哈密顿体系引进到粘弹性力学厚壁筒问题中,在辛体系下重新描述了基本问题,即建立了正则方程组。借助于积分变换,得到了拉伸、扭转和弯曲等问题的解以及有边界局部效应的解。将原问题归结为辛几何空间中的零本征值本征解和非零本征值本征解问题,从而建立了一种有效的分析问题方法和数值方法。为解决同类问题提供了一条可行的路径。     

Analytical and numerical methods of symplectic system for Stokes flow in two-dimensional rectangular domain

In this paper,a new analytical method of symplectic system.Hamiltonian system,is introduced for solving the problem of the Stokes flow in a two-dimensional rectangular domain.In the system,the fundamental problem is reduced to all eigenvalue and eigensolution problem.The solution and boundary conditions call be expanded by eigensolutions using ad.ioint relationships of the symplectic ortho-normalization between the eigensolutions.A closed method of the symplectic eigensolution is presented based on completeness of the symplectic eigensolution space.The results show that fundamental flows can be described by zero eigenvalue eigensolutions,and local effects by nonzero eigenvalue eigensolutions.Numerical examples give various flows in a rectangular domain and show effectivenees of the method for solving a variety of problems.Meanwhile.the method can be used in solving other problems.     

Multiresolution Symplectic Scheme for wave propagation in complex media

A fast adaptive symplectic algorithm named Multiresolution Symplectic Scheme (MSS) was first presented to solve the problem of the wave propagation (WP) in complex media, using the symplectic scheme and Daubechies‘ compactly supported orthogonal wavelet transform to respectively discretise the time and space dimension of wave equation. The problem was solved in multiresolution symplectic geometry space under the conservative Hamiltonian system rather than the traditional Lagrange system. Due to the fascinating properties of the wavelets and symplectic scheme, MSS is a promising method because of little computational burden, robustness and reality of long-time simulation.     

圆柱型正交各向异性弹性楔的佯谬解

圆柱型正交各向异性弹性楔体顶端受有集中力偶的经典解,当顶角满足一定关系时,其应力成为无穷大,这是个佯谬．该文在哈密顿体系下将该问题进行重新求解,即利用极坐标各向异性弹性力学哈密顿体系．在原变量和其对偶变量组成的辛几何空间求解特殊本征值的约当型本征解,从而直接给出该佯谬问题的解析解．结果再次表明经典力学中的弹性楔佯谬解对应的是哈密顿体系下辛几何的约当型解．     

基于对偶变量变分原理的完整约束系统保辛算法

基于对偶变量变分原理,选择积分区间两端位移为独立变量,构造了求解完整约束哈密顿动力系统的高阶保辛算法。首先,利用拉格朗日多项式对作用量中的位移、动量及拉格朗日乘子进行近似;然后,对作用量中不包含约束的积分项采用Gauss积分近似,对作用量中包含约束的积分项采用Lobatto积分近似,从而得到近似作用量;最后,在此近似作用量的基础上,利用对偶变量变分原理,将求解完整约束哈密顿动力系统问题转化为一组非线性方程组的求解。算法具有保辛性和高阶收敛性,能够在位移的插值点处高精度地满足完整约束。算法的收敛阶数及数值性质通过数值算例验证。