完整约束多刚体系统碰撞的矩阵解法

引言 本文讨论完整约束下多刚体系统在理论力学概念下的碰撞,碰撞是瞬时的,碰撞前后位置状态不动,仅速度状态发生变化,对完整系统,这一速度状态的变化将由广义速度来表证。 碰撞问题,理论力学书中均有介绍,但大多是质点和简单刚体系统的碰撞。对多自由度完整系统,虽有一些算例,但较讲究技巧而无一般求解途径可遁。为此,本文对Lagrange碰撞方程提出了一般的矩阵求解法,导出了三类碰撞问题矩阵通解,藉此论证了完整多刚体     

子矩阵形式的辨识方程及其求解

本文以传递函数作为特征量,建立了结构故障诊断的辨识方程．提出在频域内滑动的最小二乘法求解辨识方程,提高了诊断精度和求解效率．通过具体算例说明上述理论和方法的可靠性及工程实际应用价值．     

精细积分的非线性动力学积分方程及其解法

给出了非线性动力学积分方程的表达式,针对该方程提出了一个显式预测-校正的单步四阶精度的精细积分算法,适用于多自由度、强非线性,非保守系统.算例表明该方法精度高、计算量较少.     

Orr-Sommerfeld方程的数值解法

研究平行流动或近似平行流动,例如平面Poiseuille流及边界层流动的稳定性问题时,若采用线化小扰动理论,则最后归结为解Orr-Sommerfeld方程的特征值问题。对非线性理论来说,只要是弱非线性理论,一般也要顺序解一串Orr-Sommerfeld方程(齐次的或非齐次的)。因此解Orr-Sommerfeld方程,是研究平行或近似平行流动稳定性问题时必然要遇到的问题。在50年代以前,主要利用渐近法求Orr-Sommerfeld方程的特征值,但一般不能      

0rr－Sommerfeld方程数值解法中的复广义矩阵特征值问题

Ｏｒｒ－Ｓｏｍｍｅｒｆｅｌｄ方程的求解通常可以化为一个复广义矩阵特征值问题ＡＸ＝ωＢＸ。本文用酉变换分别约化Ａ，Ｂ为上Ｈｅｓｓｅｎｂｅｒｇ阵和上三角阵，然后利用Ｍｕｌｌｅｒ求根方法可以求出其全部特征值，其中特征多项式的值由Ｈｙｍａｎ方法给出。当仅需要判断有无不稳定模态时，利用一个简单的矩阵变换将其化为强特征值的求解问题，从而可使用最简单的幂迭代，Ｃｈｅｂｙｓｈｅｖ配置点法算例表明两种算法均快速有效。     

0rr－Sommerfeld方程数值解法中的复广义矩阵特征值问题

Ｏｒｒ－Ｓｏｍｍｅｒｆｅｌｄ方程的求解通常可以化为一个复广义矩阵特征值问题ＡＸ＝ωＢＸ。本文用酉变换分别约化Ａ，Ｂ为上Ｈｅｓｓｅｎｂｅｒｇ阵和上三角阵，然后利用Ｍｕｌｌｅｒ求根方法可以求出其全部特征值，其中特征多项式的值由Ｈｙｍａｎ方法给出。当仅需要判断有无不稳定模态时，利用一个简单的矩阵变换将其化为强特征值的求解问题，从而可使用最简单的幂迭代，Ｃｈｅｂｙｓｈｅｖ配置点法算例表明两种算法均快速有效。     

力矩分配法的方程解法

在传统力矩分配法的基础上,将渐进的计算过程转化为方程组的求解,可以更好地理解力矩分配法的分析思路及计算过程,同时提高计算效率和精度。该方法所得到的方程与矩阵位移法所得方程形式相同,但未知量完全不同,方程的推导过程更简洁且易于理解,在教学中可以启发学生对矩阵位移法分析思想的思考。     

Helmholtz方程的微分容积解法

用一种新型的数值技术－－微分容积法（Differential Cubature Method）求解二维Helmholtz方程的边值问题，几个数值算例表明，该方法稳定收敛，并具有较好的数值精度，本文方法适用于求解具有较小波数的Helmholtz方程。     

非线性非完整系统的碰撞方程

有关线性非完整系统的碰撞方程,已有文章[1]阐述清楚了。本文引入了δ_+-函数及δ~-函数,对非线性非完整系统进行了研究,並给出了该系统的碰撞方程。     

有限元病态刚度方程的解法研究

本文在文献１，２的基础上，对有限元了点位移构成准刚体运动这种情况进行了研究，提出了一种新的改性转换矩阵Ｑ，并给出了两种简单易行的计算Ｑ＾Ｔ Ａ Ｑ的方法。所提方法对直接法和迭代法求解有限元病态方程都是适用的，算例表明，得到了较好的计算结果。     

有限元病态刚度方程的解法研究

本文在文献〔１〕、〔２〕的基础上，对有限元局部节点位移构成准刚体运动这种情况进行了研究，提出了一种新的改性转换矩阵〔Ｑ〕，并给出了两种简单易行的计算〔Ｑ〕￣Ｔ〔Ａ〕〔Ｑ〕的方法。所提方法对直接法和迭代法求解有限元病态方程都是适用的，算例表明，得到了较好的计算结果。     

矩阵黎卡提方程的精细积分法

选择恰当的参数，将２￣Ｎ类算法用于代数与微分黎卡提方程。证明了算得的解是如此精确，几乎是计算机上的精确解。数例验证了该结论。     

矩阵黎卡提方程的精细积分法

选择恰当的参数，将２＾Ｎ类算法用于代数与微分黎卡提方程。证明了算得的解是如此精确，几乎是计算机上的精确解。数例验证了该结论。     

隔震结构方程的一种压缩解法

基础隔振体系中,隔震器的刚度远小于上部结构的刚度,如果未知量数目比较大,则经常导致总刚度矩阵病态。本文利用上部结构本身的振型叠加压缩未知数,然后与隔震器构建混合方程,此时形成的方程为非对称方程。大量压缩未知量后,减轻了总刚度矩阵的病态。以四层隔震框架结构分析为例,结果表明,压缩解法和正常解法的静力结果十分吻合,但用Newmark逐步积分计算时,正常直接解法累积误差引起发散,而压缩解法计算不发散。     

强非线性Duffing方程的一种解法

用线化摄动法和插值摄动法相结合的方法求解强非线性Duffing方程，有较好的精度。     

边界积分方程中超奇异积分的解法

本文对边界积分方程中所存在的超奇异积分的数值解法作了综述，并介绍了它的一些应用。     

自由面重力流的积分方程迭代解法

本文利用的奇异积分方程理论,将自由面重力流化为解析函数的Riemann-Hilbert边值问题,推出积分方程。提出以固壁边界和自由面流线长度为自变量,边界势函数为未知函数的迭代求解思想,避开了曲壁边界流速方向为未知函数的困难。提出了一种适用于直线和曲线边界自由面重力流的新的数值方法。在流量已知的条件下,证明了该方法的收敛性、稳定性。并给出了一个误差估计式。     

非线性动力学积分方程分块积分解法

对于非线性动力学方程组分块地应用精细积分算法,使其化成积分方程表达式,求解的表达式中具有相对低阶的转换矩阵,从而使精细积分更适用于多自由度、强非线性、变系数、非保守系统,针对积分方程提出了一个显示预测-校正的单步四阶精度自起步的精细积分算法。算例表明本方法是有效的。