Reissner板弯曲的辛求解体系

基于Reissner板弯曲问题的Hellinger-Reissner变分原理,通过引入对偶变量,导出Reissner板弯曲的Hamilton对偶方程组．从而将该问题导入到哈密顿体系,实现从欧几里德空间向辛几何空间,拉格朗日体系向哈密顿体系的过渡．于是在由原变量及其对偶变量组成的辛几何空间内,许多有效的数学物理方法如分离变量法和本征函数向量展开法等均可直接应用于Reissner板弯曲问题的求解．这里详细求解出Hamilton算子矩阵零本征值的所有本征解及其约当型本征解,给出其具体的物理意义．形成了零本征值本征向量之间的共轭辛正交关系．可以看到,这些零本征值的本征解是Saint-Venant问题所有的基本解,这些解可以张成一个完备的零本征值辛子空间．而非零本征值的本征解是圣维南原理所覆盖的部分．新方法突破了传统半逆解法的限制,有广阔的应用前景．     

弹性平面扇形域问题及哈密顿体系*

通过变量代换及变分原理,将平面弹性扇形域的方程导向哈密顿体系,从而可用分离变量法、本征函数展开等方法求解扇形域的分析单元,这样便可以与有限元的程序系统相结合。显示了哈密顿体系、辛数学的应用潜力。     

无穷维Hamilton算子的辛自伴性

本文运用算子扰动理论研究了无穷维Hamilton算子的共轭算子,进而得到了无穷维Hamilton算子为辛自伴算子的若干充分条件.     

两种材料组成的弹性楔的佯谬解

从Hellinger-Reissner变分原理出发,通过引入适当的变换可以将两种材料组成的弹性楔问题导入极坐标哈密顿体系,从而可以在由原变量和其对偶变量组成的辛几何空间,利用分离变量法和辛本征向量展开法求解该问题的解。在极坐标哈密顿体系下的所有辛本征值中,本征值-1是一个特殊的本征值。一般情况下本征值-1为单本征值,求解其对应的基本本征函数向量就直接给出了顶端受有集中力偶的经典弹性力学解。但当两种材料的顶角和弹性模量满足特殊关系时,本征值-1成为重本征值,同时经典弹性力学解的应力分量变成无穷大,即出现佯谬。此时重本征值-1存在约当型本征解,通过对该特殊约当型本征解的直接求解就给出了两种材料组成的弹性楔顶端受有集中力偶的佯谬问题的解。结果进一步表明经典弹性力学中弹性楔的佯谬解对应的就是极坐标哈密顿体系的约当型解。     

平面电磁弹性固体的辛对偶体系

从电磁弹性固体广义变分原理出发,将平面电磁弹性固体问题导入Hamilton体系．于是在由原变量——位移、电势和磁势以及它们的对偶变量——纵向应力、电位移和磁感应强度组成的辛几何空间,形成有效的分离变量及辛本征函数向量展开解法．求解出辛本征问题中特殊的零本征值所有本征解及其Jordan型本征解,并给出其具体的物理意义．最后求出在矩形域的两侧作用均布载荷、常电位移和常磁感应强度时的非齐次特解．     

对边简支的矩形平面弹性问题的辛本征展开定理

对来源于平面弹性问题的Hamilton算子的本征值问题进行了研究．在矩形域内含位移和应力的混合边界条件下,首先求解了相应算子的本征函数．接着,证明了本征函数系的完备性,这为施行分离变量法求解相应问题提供了可行性．最后,利用文中的辛本征展开定理获得了问题的一般解．     

混凝土断裂力学虚拟裂缝模型的半解析有限元法

利用平面扇形域哈密顿体系的方程,通过分离变量法及共轭辛本征函数向量展开法,以解析的方法推导出基于混凝土断裂力学中虚拟裂缝模型的平面裂纹解析元列式．将该解析元与有限元相结合,构成半解析的有限元法,可求解任意几何形状和荷载混凝土平面裂纹的虚拟裂缝模型计算问题．数值计算结果表明方法对该类问题的求解是十分有效的,并有较高的精度．     

斜对角无穷维Hamilton算子的点谱和特征函数系辛结构的非退化性

本文研究斜对角无穷维Hamilton算子$H=\begin{pmatrix}0&B\\C&0\end{pmatrix}$的点谱和特征函数系辛结构的非退化性, 给出斜对角无穷维Hamilton算子$H$的特征函数系具有非退化辛结构的充分必要条件. 基于此, 进一步刻画了斜对角无穷维Hamilton算子$H$的点谱分别包含于实轴、虚轴以及其它区域的充分必要条件. 最后, 以板弯曲问题和弦振动问题中导出的斜对角无穷维Hamilton算子为例, 验证了所得结论的正确性.     

Hamilton体系辛半解析法及在非均匀电磁波导中的应用

力学中的Hamilton体系需用对偶变量来描述,而电磁场正好有电场和磁场这一对对偶变量．尝试将力学中的Hamilton体系理论应用于电磁波导的分析,以横向电场和磁场作为对偶变量,将电磁波导的基本方程导向辛几何的形式．基于Hamilton变分原理, 导出横向离散的半解析系统方程, 保持体系的辛结构．以非均匀波导为例, 求解了方程的辛本征值问题, 计算结果与解析解相当吻合．     

On symplectic self-adjointness of Hamiltonian operator matrices

Symplectic self-adjointness of Hamiltonian operator matrices is studied, which is important to symplectic elasticity and optimal control. For the cases of diagonal domain and off-diagonal domain, necessary and sufficient conditions are shown. The proofs use Frobenius-Schur factorizations of unbounded operator matrices.Under additional assumptions, sufficient conditions based on perturbation method are obtained. The theory is applied to a problem in symplectic elasticity.     

基于辛对偶体系的层合板自由边缘效应的分析解

在原变量——位移和其对偶变量——应力组成的辛几何空间,建立了Pipes-Pagano模型的复合材料层合板问题的辛对偶求解体系．与传统的单类变量不同,辛对偶变量有利于同时描述层间位移连续性条件和应力平衡条件．进入辛对偶体系以后,就可以应用辛对偶体系的统一解析求解方法,如分离变量和辛本征展开的方法对层合板问题进行解析分析和求解．对层合板自由边缘效应的分析求解,验证了辛对偶体系的方法对层合板问题的分析求解是十分有效的．     

Localization and convergence of eigenfunction expansions

We state a localization principle for expansions in eigenfunctions of a self-adjoint second order elliptic operator and we prove an equiconvergence result between eigenfunction expansions and trigonometric expansions. We then study the Gibbs phenomenon for eigenfunction expansions of piecewise smooth functions on two-dimensional manifolds.     

The Fourth-order Bessel–type Differential Equation

The Bessel-type functions, structured as extensions of the classical Bessel functions, were defined by Everitt and Markett in 1994. These special functions are derived by linear combinations and limit processes from the classical orthogonal polynomials, classical Bessel functions and the Krall Jacobi-type and Laguerre-type orthogonal polynomials. These Bessel-type functions are solutions of higher-order linear differential equations, with a regular singularity at the origin and an irregular singularity at the point of infinity of the complex plane.

There is a Bessel-type differential equation for each even-order integer; the equation of order two is the classical Bessel differential equation. These even-order Bessel-type equations are not formal powers of the classical Bessel equation.

When the independent variable of these equations is restricted to the positive real axis of the plane they can be written in the Lagrange symmetric (formally self-adjoint) form of the Glazman–Naimark type, with real coefficients. Embedded in this form of the equation is a spectral parameter; this combination leads to the generation of self-adjoint operators in a weighted Hilbert function space. In the second-order case one of these associated operators has an eigenfunction expansion that leads to the Hankel integral transform.

This article is devoted to a study of the spectral theory of the Bessel-type differential equation of order four; considered on the positive real axis this equation has singularities at both end-points. In the associated Hilbert function space these singular end-points are classified, the minimal and maximal operators are defined and all associated self-adjoint operators are determined, including the Friedrichs self-adjoint operator. The spectral properties of these self-adjoint operators are given in explicit form.

From the properties of the domain of the maximal operator, in the associated Hilbert function space, it is possible to obtain a virial theorem for the fourth-order Bessel-type differential equation.

There are two solutions of this fourth-order equation that can be expressed in terms of classical Bessel functions of order zero and order one. However it appears that additional, independent solutions essentially involve new special functions not yet defined. The spectral properties of the self-adjoint operators suggest that there is an eigenfunction expansion similar to the Hankel transform, but details await a further study of the solutions of the differential equation.     

对偶变量块体混合元及其位移元的收敛性和精度分析

弹性力学Hamilton正则方程和Hamilton混合元的等效刚度系数矩阵,均具有直观的辛特性.基于H R变分原理和弹性力学保辛理论建立的对偶变量块体混合元,其等效刚度系数矩阵同样具有直观的辛特性.根据对偶变量块体混合元列式,可直接建立问题的控制方程,进行混合法求解.同时,通过对偶变量块体混合元列式可以导出对偶变量块体位移元列式,建立问题的控制方程后,可先求位移的解.数值实例表明:线性8结点对偶变量块体位移减缩积分元的各力学量的收敛速度均衡、收敛过程稳定、结果精度高,其应力变量的收敛速度与传统的20结点位移协调减缩积分元接近.对偶变量块体位移元具有普适性.     

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

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

PT-对称算子的表示及非正交量子态的区分

经典量子系统中的哈密尔顿为自伴算子，这不仅保证了系统能量本征值全部为实数，而且相应的本征态（单位长度的特征向量）构成了状态空间的一组正规正交基.然而存在一类PT-对称的物理系统，哈密尔顿的自伴性（共轭转置）被物理的PT-对称性所代替.一个完整的PT-对称哈密尔顿，其谱全部为实数且能构造一个合理的CPT-内积.本文研究一类PT-对称算子.固定时间反演算子T，得到宇称算子P的矩阵表示，进而给出每一组PT-对称哈密尔顿的具体表示形式.作为应用，选择一组确定的{P，T}算子，及PT-对称的哈密尔顿，给出两个在传统量子力学中不正交的量子态区分的刻画.     

Spectral theory of discrete linear Hamiltonian systems

This paper is concerned with spectral problems for a class of discrete linear Hamiltonian systems with self-adjoint boundary conditions, where the existence and uniqueness of solutions of initial value problems may not hold. A suitable admissible function space and a difference operator are constructed so that the operator is self-adjoint in the space. Then a series of spectral results are obtained: the reality of eigenvalues, the completeness of the orthogonal normalized eigenfunction system, Rayleigh's principle, the minimax theorem and the dual orthogonality. Especially, the number of eigenvalues including multiplicities and the number of linearly independent eigenfunctions are calculated.     

无穷维Hamilton 算子广义本征函数系的完备性及其在弹性力学中的应用

本文利用无穷维Hamilton 算子的结构特性, 得到由算子的基本本征函数和若当型本征函数构成的广义本征函数系在Cauchy 主值意义下完备的充分必要条件. 进而将结果应用于弹性力学中的板弯曲问题. 相应结论为Hamilton 体系下的分离变量法(弹性力学求解新体系) 提供了理论保证.     

The Limiting Absorption Principle for Elastic Wave Propagation Problems in Perturbed Stratified Media ℝ3

We consider the self-adjoint operator governing the propagation of elastic waves in perturbed stratified media ℝ³ with free boundary–interface conditions. In this paper we establish the limiting absorption principle for this self-adjoint operator in appropriate Hilbert space. The proof of the limiting absorption principle is based on the division theorem which is proved by means of eigenfunction expansions for the self-adjoint operator governing the propagation of elastic waves in unperturbed stratified media ℝ³.     

Stochastic discrete Hamiltonian variational integrators

Variational integrators are derived for structure-preserving simulation of stochastic Hamiltonian systems with a certain type of multiplicative noise arising in geometric mechanics. The derivation is based on a stochastic discrete Hamiltonian which approximates a type-II stochastic generating function for the stochastic flow of the Hamiltonian system. The generating function is obtained by introducing an appropriate stochastic action functional and its corresponding variational principle. Our approach permits to recast in a unified framework a number of integrators previously studied in the literature, and presents a general methodology to derive new structure-preserving numerical schemes. The resulting integrators are symplectic; they preserve integrals of motion related to Lie group symmetries; and they include stochastic symplectic Runge–Kutta methods as a special case. Several new low-stage stochastic symplectic methods of mean-square order 1.0 derived using this approach are presented and tested numerically to demonstrate their superior long-time numerical stability and energy behavior compared to nonsymplectic methods.