非齐次动力学方程的一种精细积分单步方法

针对非齐次动力学方程■,结合精细积分法和微分求积法,利用同阶的显式龙格-库塔法对计算过程中待求的v_(k+i/s)(i=1,2,…,s)进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算e~(Ht),而Duhamel积分项采用s级s阶的时域微分求积法,计算格式统一且易于编程,可灵活实现变阶变步长。仿真结果表明,与其他单步法及预估校正-辛时间子域法进行数值比较,该方法具有高精度、高效率及良好的稳定性,在求解大规模动力系统时间响应问题中具有较大的优势。     

飞行物体受瞬态荷载作用的精细逐步积分

当前的精细逐步积分方法不能直接处理具有刚体自由度结构的瞬态响应计算。本文就此提出了应对策略，仍然可用较大的步长而获得精确解答。     

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

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

非齐次动力方程Duhamel项的精细积分

提出了不需要矩阵求逆运算的求解Duhamel积分项的精细积分方法.通过将精细积分法的关键思想--加法定理和增量存储--直接应用于Duhamel积分响应矩阵的求解,可给出当非齐次项分别为多项式、正弦/余弦以及指数函数等基本形式时Duhamel积分在计算机上的精确解.特别的,该算法不依赖于系统矩阵(或相关矩阵)的形态.当系统矩阵奇异或接近奇异时,其优越性更为显著.算例验证了该算法的有效性.     

解对流方程的子域精细积分并行算法

基于子域精细积分的思想,针对对流方程初边值问题,首先提出了含参数a＞0的一族三层显格式和一族二层隐格式,它们的局部截断误差分别为O(a△t+△t2+△x2)和O(α△t+△t+△x2).当参数a≥(In△t-ln△x)/2△t时三层显格式是稳定的,而二层隐格式则对所有的参数α＞0都是无条件稳定的.然后,以二层隐格式为基础,设计了一种交替分组显武迭代(AGEI)方法,并证明了该迭代过程的收敛性.由于三层显格式和AGEI方法的整个计算过程都是显式的,所以非常适合于并行计算.文末的数值算例表明,上述方法具有很高的精确度和良好的实用性.     

岩土弹塑性有限元分析的隐式积分弹性刚度法

本文提出了用于岩土弹塑性有限元分析的隐式积分弹性刚度算法。该算法既具有隐式积分法精度好，效率高，无条件稳定等优点，也具有弹性刚度矩阵正定、对称的特点，更重要的是它避免了传统切线刚度法在处理岩土非相关联塑性流动和屈服面“角角”所遇到的非对称性和奇异性计算问题。通过算例分析了该算法的精度、效率。     

非均质问题中的无网格边界单元法

介绍了一种不需要内部网格计算非均匀介质问题的边界元算法.该算法是建立在一种能将任何区域积分转换成边界积分的径向积分转换法基础上,首先用对应各向同性问题的基本解来建立以正规化位移表示的非均质问题的积分方程,然后用径向积分转换法将出现在积分方程中的区域积分转换成边界积分,从而形成不需要使用内部网格来计算区域积分的纯边界元算法.与其它无网格法相比,此方法需要很少的内部点,有些问题甚至不需要内部点都能得到满意的结果,因此,可以计算大型的三维非均匀介质工程问题.由于此方法继承了边界元和无网格算法的优点,因而具有广阔的发展前景.     

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.     

有限深水域中三维频域自由面格林函数及其导数的计算方法

三维频域自由面格林函数及其偏导数的数值计算是海洋结构物水动力分析的难点。本文对有限水深格林函数及其偏导数的数值算法作出两点改进：①对原函数的级数表达式进行变形,提出适用于数值计算的形式,解决高频率情况下级数解的计算失真问题;②对Gauss-Laguerre积分法作出改进,解决高频率大水深情况下积分解的数值溢出问题。计算结果表明：本文提出的改进算法可以有效提高有限水深格林函数及其偏导数的数值精度,对于研究较大吃水的浮体或潜体的高频垂荡问题具有重要意义。     

校正梯形公式误差的精确表示

利用分部积分法和推广的积分第一中值定理,给出了具5次代数精度校正梯形公式误差的精确表示.