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

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

Plane elasticity in sectorial domain and the Hamiltonian system

ＰＬＡＮＥＥＬＡＳＴＩＣＩＴＹＩＮＳＥＣＴＯＲＩＡＬＤＯＭＡＩＮＡＮＤＴＨＥＨＡＭＩＬＴＯＮＩＡＮＳＹＳＴＥＭＺｈｏｎｇＷａｎ－ｘｉｅｔａｐ（钟万勰）（ＤａｌｉａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ．Ｄａｌｉａｎ）（ＲｅｃｅｉｖｅｄＤｅｃ．...     

A kinetic interpretation of the section‐averaged Saint‐Venant system for natural river hydraulics

The classical Saint‐Venant system is well suited for the modeling of dam breaks, hydraulic jumps, reservoir emptying, flooding etc. For many applications, the extension of the Saint‐Venant system to the case of non‐rectangular channels is necessary and this section‐averaged Saint‐Venant system exhibits additional source terms. The main difficulty of these equations consists of the discretization of these source terms. In this paper we propose a kinetic interpretation for the section averaged Saint‐Venant system and derive an associated numerical scheme. The numerical scheme—2nd order in space and time—preserves the positivity of the water height, and is well‐balanced. Numerical results including comparisons with analytic and experimental test problems illustrate the accuracy and the robustness of the numerical algorithm. Copyright © 2010 John Wiley & Sons, Ltd.     

Hamiltonian systems in elasticity and their variational principles

ＨＡＭＩＬＴＯＮＩＡＮＳＹＳＴＥＭＳＩＮＥＬＡＳＴＩＣＩＴＹＡＮＤＴＨＥＩＲＶＡＲＩＡＴＩＯＮＡＬＰＲＩＮＣＩＰＬＥＳＷａｎｇＺｈｉ－ｇｕｏ（王治国）（ＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅｏｆＶｉｂｒａｔｉｏｎＥｎｇｉｎｅｅｒｉｎｇ,ＮａｎｊｉｎｇＵｎ...     

Global Constraints for Stress Constrained Materials: The Problem of Saint Venant

Based on the global constraint principle of Antman and Marlow, a new solution of Saint Venant's problem is proposed. The solutions for the six fundamental cases of loading in terms of stress are obtained with relative ease and converge to the classical Saint Venant's solution as the length of the beam is increased. It is also shown that the assumptions of a special technical rod theory are coherent with the requirements of the global constraint theory for the Saint Venant cylinder.     

A simple criterion for finite time stability with application to impacted buckling of elastic columns

The existing theories of finite-time stability depend on a prescribed bound on initial disturbances and a prescribed threshold for allowable responses. It remains a challenge to identify the critical value of loading parameter for finite time instability observed in experiments without the need of specifying any prescribed threshold for allowable responses. Based on an energy balance analysis of a simple dynamic system, this paper proposes a general criterion for finite time stability which indicates that finite time stability of a linear dynamic system with constant coefficients during a given time interval [0, t _f ] is guaranteed provided the product of its maximum growth rate (determined by the maximum eigen-root p₁ >0) and the duration t _f does not exceed 2, i.e., p₁t _f <2. The proposed criterion (p₁t _f =2) is applied to several problems of impacted buckling of elastic columns: (i) an elastic column impacted by a striking mass, (ii) longitudinal impact of an elastic column on a rigid wall, and (iii) an elastic column compressed at a constant speed (“Hoff problem”), in which the time-varying axial force is replaced approximately by its average value over the time duration. Comparison of critical parameters predicted by the proposed criterion with available experimental and simulation data shows that the proposed criterion is in robust reasonable agreement with the known data, which suggests that the proposed simple criterion (p₁t _f =2) can be used to estimate critical parameters for finite time stability of dynamic systems governed by linear equations with constant coefficients.     

Second-order effects and Saint Venant’s principle in the torsion problem of a nonlinear elastic rod

The torsion problem of a circular nonlinear elastic rod loaded by end moments is considered. The solution constructed by the method of successive approximations taking into account second-order effects is compared with the solution obtained by a semi-inverse method. It is shown that the dead-loading assumption breaks the symmetry of the Cauchy stress tensor in a certain region. A refined formulation of Saint Venant’s principle is proposed for the problem of determining integral strain characteristics. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 6, pp. 129–136, November–December, 2006.     

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

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

Hamilton体系下环扇形域的Stokes流动问题

基于极坐标下Stokes流的基本方程,将径向坐标模拟为时间坐标,推导了Hamilton体系下Stokes流动问题的对偶方程,采用本征向量展开法对环扇形域Stokes流动问题进行了分析,并给出了相应的实际算例,其结果说明了本文方法的有效性。     

High order symplectic conservative perturbation method for time-varying Hamiltonian system

This paper presents a high order symplectic conservative perturbation method for linear time-varying Hamiltonian system.Firstly,the dynamic equation of Hamiltonian system is gradually changed into a high order perturbation equation,which is solved approximately by resolving the Hamiltonian coefficient matrix into a "major component" and a "high order small quantity" and using perturbation transformation technique,then the solution to the original equation of Hamiltonian system is determined through a series of inverse transform.Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes,the transfer matrix is a symplectic matrix;furthermore,the exponential matrices can be calculated accurately by the precise time integration method,so the method presented in this paper has fine accuracy,efficiency and stability.The examples show that the proposed method can also give good results even though a large time step is selected,and with the increase of the perturbation order,the perturbation solutions tend to exact solutions rapidly.     

Discharge estimation under uncertainty using variational methods with application to the full Saint‐Venant hydraulic network model

Estimating river discharge from in situ and/or remote sensing data is a key issue for evaluation of water balance at local and global scales and for water management. Variational data assimilation (DA) is a powerful approach used in operational weather and ocean forecasting, which can also be used in this context. A distinctive feature of the river discharge estimation problem is the likely presence of significant uncertainty in principal parameters of a hydraulic model, such as bathymetry and friction, which have to be included into the control vector alongside the discharge. However, the conventional variational DA method being used for solving such extended problems often fails. This happens because the control vector iterates (i.e., approximations arising in the course of minimization) result into hydraulic states not supported by the model. In this paper, we suggest a novel version of the variational DA method specially designed for solving estimation‐under‐uncertainty problems, which is based on the ideas of iterative regularization. The method is implemented with SIC², which is a full Saint‐Venant based 1D‐network model. The SIC² software is widely used by research, consultant and industrial communities for modeling river, irrigation canal, and drainage network behavior. The adjoint model required for variational DA is obtained by means of automatic differentiation. This is likely to be the first stable consistent adjoint of the 1D‐network model of a commercial status in existence. The DA problems considered in this paper are offtake/tributary estimation under uncertainty in the cross‐device parameters and inflow discharge estimation under uncertainty in the bathymetry defining parameters and the friction coefficient. Numerical tests have been designed to understand identifiability of discharge given uncertainty in bathymetry and friction. The developed methodology, and software seems useful in the context of the future Surface Water and Ocean Topography satellite mission. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

THE HAMILTONIAN SYSTEM AND COMPLETENESS OF SYMPLECTIC ORTHOGONAL SYSTEM

I.IntroductionThemethodofseparationofvariablesisimportanttosolvethesoluti0n0fprobIem0fmathematicalphysics,butmanyproblen1sofmathematicalphysicscannotseparatet'ariab1es,thereforeitrestrictstheranget0appIicatemethodofseparationofvariable.Inthepaperlll,Zhong…     

哈密顿体系下机电耦合问题的基本方程

利用多变量变分原理,针对具有机电耦合特性的压电材料,在哈密顿体系下推导出了机电耦合问题的对偶方程．求出了问题的零本征解和非零本征解的一般表达式。并就“二维压电平板对边受均布载荷”这一具体问题进行了完整的分析计算,得出了应力和电位移的具体解答,结果与一般弹性力学所得的结论相吻合,而且精度极高。表明哈密顿体系适用于压电材料的力学分析。     

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

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

Comparison of finite element methods for the St. Venant equations

Finite element schemes for hyperbolic systems are applied to the St. Venant equations for one-dimensional, unsteady, open channel flow. The comparative performances of the characteristic-dissipative-Galerkin, Taylor-Galerkin and least squares finite element schemes are assessed by means of linear Fourier analysis and solution of idealized non-linear wave propagation problems. Of particular interest is the behaviour of these schemes for the regressive wave component in both subcritical and supercritical flows. To assess the quality of the basic solution, the methods are compared without any additional artificial diffusion or shock-capturing formulations. The balanced treatment of both wave components in the characteristic-dissipative-Galerkin method is illustrated. Also, the method displays little sensitivity to parameters variations. The Taylor-Galerkin scheme provides good solutions, although oscillations due to wave dispersion and minimal diffusion of the regressive wave are displayed. Also, this method is somewhat sensitive to the time step increment. The least squares method is considered unsuitable for unsteady, open channel flow problems owing to its inability to propagate a regressive wave in a supercritical flow.     

空腔内粘性流问题与哈密顿体系方法

以双板驱动空腔粘性流问题为研究对象,根据其特点建立了哈密顿体系下的对偶正则方程,将问题归结为辛体系下的本征值问题.利用辛本征解空间的完备性、正交性和展开理论,形成一套封闭的求解问题方法.算例的数值结果揭示了一些空腔流动的特点.同时这种方法也为研究其他问题提供了一条思路.     

Completeness of system of root vectors of upper triangular infinite-dimensional Hamiltonian operators appearing in elasticity theory

This paper deals with a class of upper triangular infinite-dimensional Hamiltonian operators appearing in the elasticity theory.The geometric multiplicity and algebraic index of the eigenvalue are investigated.Furthermore,the algebraic multiplicity of the eigenvalue is obtained.Based on these properties,the concrete completeness formulation of the system of eigenvectors or root vectors of the Hamiltonian operator is proposed.It is shown that the completeness is determined by the system of eigenvectors of the operator entries.Finally,the applications of the results to some problems in the elasticity theory are presented.     

在哈密顿体系下分析非线性动力学问题

首先将ｎ维未知向量ｑ的二阶非线性力系统Ｍｑ＋Ｇｑ＋Ｋｑ＝Ｆ（ｑ,ｑ,ｔ）转化为与其等价的２ｎ维未知向量ｖ的一阶微分方程ｖ＝Ｈｖ＋ｆ（ｖ,ｔ）,其中非线性部分ｊｉ（ｖ,ｔ）＝０（ｉ＝１,．．．ｎ）,ｆｉ（ｖ,ｔ）＝Ｆｉ－ｎ（ｑ,ｑ,ｔ）（ｉ＝ｎ＋１,．．．２ｎ）;然后给出一种求解ｖ的逐步积分公式,从而将精细积分法进一步推广应用到非线性动力学问题。算例表明本方法的计算量较小且结果合理可靠。