HIGH PERFORMANCE SPARSE SOLVER FOR UNSYMMETRICAL LINEAR EQUATIONS WITH OUT-OF-CORE STRATEGIES AND ITS APPLICATION ON MESHLESS METHODS

A new direct method for solving unsymmetrical sparse linear systems(USLS) arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin (MLPG) method need to solve large USLS. The proposed solution method for unsymmetrical case performs factorization processes symmetrically on the upper and lower triangular portion of matrix, which differs from previous work based on general unsymmetrical process, and attains higher performance. It is shown that the solution algorithm for USLS can be simply derived from the existing approaches for the symmetrical case. The new matrix factorization algorithm in our method can be implemented easily by modifying a standard JKI symmetrical matrix factorization code. Multi-blocked out-of-core strategies were also developed to expand the solution scale. The approach convincingly increases the speed of the solution process, which is demonstrated with the numerical tests.     

离散时间系统Ляпунов函数的构造与稳定性判定

本文旨在以矩阵的相似变换理论为基础,以电子计算机为工具,提供一种无需求解离散?矩阵方程,而直接由系统矩阵来构造离散时间系统?函数和判定系统稳定性的方法.     

轴对称薄壁结构自由振动的边界元分析

采用双互易法分析薄壁轴对称结构自由振动的特征频率以及特征模态.首先,采用径向基函数插值域积分里的位移,利用双互易法将域积分转化为子午面边界的积分.然后,将边界物理量、基本解和特解展开为傅里叶级数,沿环向积分后得到的边界积分方程可用于轴对称结构带体积力问题和受非对称载荷的动力学分析,其积分域为轴对称结构子午面边界上的线积分,进一步降低了问题的维度和离散的难度.文章详细探讨了源点处于对称轴的特殊情况,根据基本解和特解的退化形式,针对无体积力和有体积力分别给出了处理奇异矩阵的方案.对于薄壁结构,采用双曲正弦变换处理近奇异积分有效提高积分精度.最后将双互易法和双曲正弦变化应用于薄壁轴对称结构带体积力的静力学和自由振动分析.数值结果表明,文章提出的处理奇异矩阵的方法能够有效处理源点处于对称轴的情况;当圆筒厚高比为$10^{-3}$,边界元计算的特征频率的相对误差为$10^{-3}$,且优于有限元的结果.      

附加约束阻尼层对星箭系统动特性的影响分析

在星箭连接支架外表面附加约束阻尼层可抑制星箭对接界面的振动,从而降低卫星在发射过程中所承受的振动和噪声载荷.本文将有限元技术与矩阵特征值摄动分析方法相结合,就卫星支架局部修改(附加约束阻尼层)对星箭系统动态特性的影响进行了分析.以已有星箭系统结构的模态信息为基础,给出了支架部位附加约束阻尼层后星箭系统动态特性变化的定量计算方法,并证明了该减振技术能够不明显的改变原星箭系统的固有特性.     

A NOVEL SEMI-ANALYTICAL METHOD FOR SOLVING ACOUSTIC RADIATION FROM LONGITUDINALLY STIFFENED INFINITE NON-CIRCULAR CYLINDRICAL SHELLS IN WATER

Based on the extended homogeneous capacity high precision integration method and the spectrum method of virtual boundary with a complex radius vector, a novel semi-analytical method, which has satisfactory computation effectiveness and precision, is presented for solving the acoustic radiation from a submerged infinite non-circular cylindrical shell stiffened by longitudinal ribs by means of the Fourier integral transformation and stationary phase method. In this work,besides the normal interacting force, which is commonly adopted by some researchers, the other interacting forces and moments between the longitudinal ribs and the non-circular cylindrical shell are considered at the same time. The effects of the number and the size of the cross-section of longitudinal ribs on the characteristics of acoustic radiation are investigated. Numerical results show that the method proposed is more efficient than the existing mixed FE-BE method.     

半刚性钢框架的内力分析

采用二阶非线性分析方法分析和设计半刚性钢框架，包括连接的柔性以及构件的几何非线性的影响，提出了半刚性钢框架中梁柱单元刚度矩阵和半刚性梁的单元刚度矩阵；推导了半刚性梁在集中荷载，均布荷载，线性荷载作用下的固端弯矩的求解公式；连接的柔性对无支撑框架的侧移有很大的影响，设计时通过变化连接的刚度以平衡梁的跨中和端弯矩。     

容量矩阵法在涡方法中的应用

本文将容量矩阵技术与快速Fourier变换相结合求解非矩形域的Poisson方程,再与格子涡方法相结合,发展了一套求解复杂外形的快速涡方法。作为算例,计算了有厚度平板和楔形体的分离流动,得到了满意的结果。特别是在垂直平板绕流中,模拟出了由于剪切层不稳定性引起的小涡结构。     

改进的小损伤结构动力学有限元建模方法

在遗传算法或神经网络方法识别结构损伤位置和程度时,都是基于少量的在线测量的损伤结构振动响应数据和大量的模型仿真数据来实现的,因而建立高效和精确的损伤结构动力学有限元模型,以便仿真获得损伤结构的大量动力学响应数据是十分重要的基础前提工作。本文针对ANSYS结构分析软件在建立结构小损伤有限元动力学模型存在两个关键问题,结构损伤处直接网格划分的计算结果误差和网格节省问题,以结构损伤振动检测的实际需要为出发点,提出了建立含小损伤结构的ANSYS动力学建模技术,研究了结构局部小损伤及其位置与所在处单元刚度矩阵变化的数量关系。     

研究正定矩阵幂函数单调性的分析方法

文章用数学分析的方法研究了对于正定矩阵的幂次方的单调性.具体地,当指数为小于1的正数时,正定矩阵幂函数是严格单调增加的,而当指数大于1时,正定矩阵幂函数不具有单调性.     

线性矩阵方程与线性矩阵方程组

研究了线性矩阵方程(组)的解的结构理论和求解方法.