多体系统动力学微分／代数方程组数值方法

多体系统动力学微分／代数混合方程组又称Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方程，是近十年来动力学和计算数学领域研究的热点之一．本文介绍这两个领域中引入的传统的数值积分方法与新的理论．     

热压电变分原理

热释电介的许多传感器及机敏结构或中的关键材料，本文系统地讨论了在基础理论和数值计算中极具重要地位的各类变分原理，包括准静态变分原理、动态变分原理和关于固有频率的变分原理。最后建立了关于静态压电要板的变分原理，并由此导出了各向异性压电板的控制方程及边界条件。     

力学系统的扩散运动及量子化

研究位形由扩散过程定义的力学系统．利用变分的路径计算给出系统广义坐标表示的动力学方程，该方程连同连续方程描述系统量子波动随时间演化的概率进程．文末的表示定理指出，在满足适当正规条件的扩散过程和Ｓｃｈｒｄｉｎｇｅｒ方程的解（波函数）之间存在一一对应关系     

基于吸力量测确定膨胀土活动带和裂隙深度

采用M itchell公式和裂隙扩展深度方程两种吸力法确定安康地区膨胀土大气影响深度和裂隙开展深度。其一通过对安康地区两处天然边坡开挖观测井,利用张力计进行不同深度处吸力值的现场量测,根据M itchell提出公式计算大气影响深度;其二根据非饱和土抗拉强度公式,建立膨胀土裂隙扩展深度方程,利用基质吸力量测结果求其理论解。结果表明,安康地区膨胀土吸力变化曲线随深度增加变幅减小,呈“波浪式”推移。M itchell公式确定安康地区膨胀土的大气影响深度为3.35m以内,裂隙深度方程确定裂隙开展深度为3.06³.14m。利用M itchell公式计算大气影响深度与膨胀土断裂理论公式确定的裂隙开展深度结果接近。     

SYMPLECTIC STRUCTURE OF POISSON SYSTEM

When the Poisson matrix of Poisson system is non-constant, classical symplectic methods, such as symplectic Runge-Kutta method, generating function method, cannot preserve the Poisson structure. The non-constant Poisson structure was transformed into the symplectic structure by the nonlinear transform.Arbitrary order symplectic method was applied to the transformed Poisson system. The Euler equation of the free rigid body problem was transformed into the symplectic structure and computed by the midpoint scheme. Numerical results show the effectiveness of the nonlinear transform.     

松质骨弹性模量计算的均匀化方法

对松质骨建立了六种单胞微观结构模型，采用均匀化方法和有限元方法计算松质骨的宏观等效弹性模量。给出了六种单胞模型的松质骨弹性模量与材料密度(体分比）的关系，与实验数据进行了对比，分析了不同微观结构模型在不同骨骼中的应用。结果表明，本文方法及六种单胞模型可以对松质骨微观结构和材料性能进行有效的模拟计算。同时本文又着重对松质骨的宏观等效弹性模量与体分比的指数关系进行了探讨。     

THREE HIGH-ORDER SPLITTING SCHEMES FOR 3D TRANSPORT EQUATION

IntroductionWith the development of modern industry, various pollutants discharge into the air,rivers, lakes and oceans, which makes the environmental qualities worse and has bad effectson the mankind’s health and the sustained development of industry an…     

动量矩变量描述的陀螺体永久转动的分岔特性

本文建立以动量矩分量和Ｄｅｐｒｉｔ正则变量为独立变量的陀螺体动力学方程，用以描述无力矩非对称陀螺体的姿态运动。在动量矩空间内讨论了陀螺体永久转动轴的数目和稳定性随转子转速变化的分岔特性。     

各向同性和横观各向同性材料的统一点力解

本文根据横观各向同性弹性力学的通解获得了无限体的点力解，由它可以直接退化到各向同性情形的Ｋｅｌｖｉｎ解，利用这个点力解编制的边界元法程序，适用于横观各向同性材料也适用于各向同性材料，因此是真正的统一点力解。还用边界元法计算了两个数值例题。     

AuQ+碰撞氦原子单电离的全微分截面

利用三体扭曲波计算了3.6 MeV/amu的Au^Q+ (Q=24,53)碰撞氦原子单电离的全微分截面,其结果与实验数据和其它理论进行了比较,表明全微分截面受到高电荷入射粒子的强烈影响,敲出的电子被“拖拽”沿着向前的方向. 但是,所有的理论结果都不能呈现实验上全微分截面的独特的向前峰.