大型结构动力响应的状态方程的Krylov精细时程积分法

提出了一种新的精细时程积分法来求解大型动力系统.结合Krylov子空间法、培德级数近似以及一般载荷的维数扩展法,进一步提高精细时程积分法的计算效率.利用维数扩展法避免计算微分方程特解,并可处理任意载荷.对于大型动力系统,通过Krylov子空间的降维分析将问题转化到一个子空间,计算效率得到极大提高.对于迭代次数N的选择作了详细讨论,进一步提高了计算效率.     

增维精细积分法的适用范围研究

对线性定常结构动力系统提出的增维精细积分法,能够将非齐次动力方程转化为齐次动力方程,不用对状态矩阵求逆就能方便高效地求解出结构的动力响应。本文在仔细分析增维精细积分法性质的基础上,提出了其适用条件,进一步拓宽了其应用范围,并给出了将荷载项展开成傅里叶级数时,相应增维精细积分法的表达式。同时,在一个时间步长内,通过对非齐次项作线性化假设,成功地将增维精细积分法应用到了非线性动力分析领域。本文方法计算格式统一,易于编程,具有很高的计算效率。数值算例证明了本文方法的有效性。     

非线性动力方程的增维精细积分法

对线性定常结构的动力系统提出的精细积分法，能得到在数值上逼近于精确解的结果。但是对于非齐次动力方程却涉及到矩阵求逆的困难，而且通常与时间有关的非齐次项不能进入精细积分的细化过程。采用增维的方法，将非齐次动力方程化为齐次方程，在实施精细积分的过程中不必进行矩阵求逆。这种处理方法对于程序实现和提高数值计算的稳定性十分有利，而且在大型问题中可明显提高计算效率，数值算例显示本文方法是有效的。     

计算结构动力响应的分段精细时程积分方法

利用钟万勰等发展的精细时程积分方法，提出了解线性定常结构动力系统响应的分段精细进程积分方法，它能适用于各种激励作用下系统的动力响应。对于载荷项采用线性和两次多项式进行拟合，采用精细时程积分方法和叠代方法对动力响应进行计算，与传统的离散积分方法如纽马克方法和威尔逊方法等及状态方程直接积分方法进行数值比较，本方法具有很高的精度和计算效率。     

求大规模非线性随机动力系统概率密度函数解的子空间法

论文给出了一个求大规模非线性随机动力系统响应概率密度函数解的新方法,称之为子空间法.考察了此方法在求解大规模带位移项参数激励非线性随机动力系统响应概率密度函数的有效性.这里的概率密度函数解由Fokker-Planck-Kolmogorov(FPK)方程控制.该方法是基于将非线性随机动力系统状态变量空间分成两个子空间,然后在其中一个子空间上对FPK方程进行积分,采取一定措施后得到低维的FPK方程.该低维的FPK方程的维数可以人为确定,也可以取为二维,从而可以用现有的求解低维FPK方程的方法求得所需的概率密度函数解.文中给出了算例,用数值结果验证了子空间法的有效性.论文是采用作者曾提出的指数多项式闭合(EPC)法求解由子空间法降维的FPK方程.     

悬索桥颤振稳定性分析的精细时程积分法

研究精细时程积分法在悬索桥颤振稳定性分析中的应用,首先,将 气流整个作为一个系统,组集系统关于模态广义坐标的状态空间方程,然后,应用精细时程积分法计算状态的向量的时程响应,根据状态向量时程响应的对数衰减率判断系统的颤振稳定性,最后,以英国塞文悬索桥为数值算例,验证了本文方法的正确性。     

非参数概率系统非平稳随机响应研究

对复杂系统的测量和知识的有限性造成的不确定很难完全由随机参数模型进行描述,采用随机矩阵模型更具有一般性和合理性.本文应用随机矩阵模拟不确定线性动力系统有限元模型中质量阵、阻尼阵和刚度阵的概率不确定性.综合运用虚拟激励法和精细时程积分法建立了非参数概率系统非平稳随机响应的有效计算方法.数值模拟结果表明,对于高精度制造,模型的不确定性是不能忽略的.本文提出的算法为此类问题求解提供了一条有效途径.     

精细积分法应用于地震碰撞力反应谱计算研究

相邻结构的碰撞时程分析是计算地震碰撞力反应谱的基础。结构碰撞时程分析要求采用稳定性好、精度高及计算效率高的数值分析方法。精细积分法将二阶动力微分方程通过增元降阶的方式转换成Hamilton对偶变量体系,得到了动力微分方程的精确解。基于此,本文将精细积分法引入结构碰撞时程分析及地震碰撞力反应谱计算中。在推导精细积分法公式的基础上,在MATLAB环境下编制了结构碰撞时程分析程序和碰撞力反应谱计算程序,并实现了碰撞力反应谱程序的并行化。经算例验证,精细积分法应用于结构碰撞时程分析及地震碰撞力反应谱计算是可行的,程序计算结果准确。     

结构非平稳随机响应分析的快速虚拟激励法

虚拟激励法能够方便地应用于结构非平稳随机响应分析,但在每个离散频点处都涉及到虚拟激励作用下动力方程的时程积分,对于大型复杂结构,其计算量是难以接受的。将结构动力方程写成状态方程形式,采用精细积分法对状态方程进行数值求解,导出了结构动力响应关于离散时刻处激励的显式线性表达式。利用这一显式表达式,只需要变换离散时刻处的激励数值,就可以方便快捷地求出新的激励作用下的结构动力响应。效率分析和数值算例表明,相对于传统虚拟激励法,本文提出的改进算法在求解非平稳激励下结构随机振动方面具有更高的计算效率。     

结构动力方程的增维精细积分法

对线性定常结构动力系统提出的精细积分方法,能够得到在数值上逼近于精确解的结果,但对于非齐次动力方程涉及到矩阵求逆的困难。提出采用增维的办法,将非齐次动力方程转化为齐次动力方程,在实施精细积分过程中不必进行矩阵求逆,这种方法对于程序实现和提高数值稳定性十分有利,而且在大型问题中计算效率较高,从而改进了精细积分方法的应用,数值例题显示了本文方法的有效性。     

RTM充模过程数值模拟的隐式有限元算法

建立了基于欧拉方法描述树脂传递模塑(RTM)工艺充模过程的基本数学方程,并采用有限元隐式时间积分方法对基本方程进行了数值求解．编制了基于隐式有限元算法及传统有限元控制体算法的程序,通过具体算例比较了这两种算法的优缺点．与传统的有限元控制体法相比,该文提出的隐式有限元算法能节省计算时间,特别适合于单元、节点数目多的情况．隐式有限元算法是一种纯有限元方法,不需要使用控制体积技术,采用该算法计算出的流动前沿与时间步长无关。     

瞬态热传导方程精细积分方法中对称性的利用

采用精细积分法求解瞬态热传导方程时,对指数矩阵进行变换后使其具有对称性,利用这一特性可使存贮量和计算量降低一半。变换后指数矩阵的带宽特性不变,采用子域精细积分可进一步提高算法的计算与存储效率。     

改进的泡泡布点方法研究

泡泡布点方法能够在复杂区域内生成高质量的节点集，但是其计算效率仍有待提高，为此本文做了两方面的改进。一是种子填充算法思想应用于节点的初始布置中，从而生成数目合适的节点集，省去了原布点方法中的节点增删过程，有效地节约模拟时间；二是简化了泡泡运动模拟的控制方程，使得运动模拟更加简单。数值算例表明，改进的泡泡布点方法生成的网格平均质量均高于0.94,且计算效率相比原布点方法提高90%以上。因此，改进的泡泡布点方法是一种高质量、高效率且适宜求解大规模问题的节点布置方法。     

An adaptive algorithm of precise integration for transient analysis

This paper presents an improved precise integration algorithm for transient analysis of heat transfer and some other problems. The original precise integration method is improved by means of the inverse accuracy analysis so that the parameterN, which has been taken as a constant and an independent parameter without consideration of the problems in the original method, can be generated automatically by the algorithm itself. Thus, the improved algorithm is adaptive and the accuracy of the algorithm is not dependent on the length of the time step in the integration process. It is shown that the numerical results obtained by the method proposed are more accurate than those obtained by the conventional time integration methods such as the difference method and others. Four examples are given to demonstrate the validity, accuracy and efficiency of the new method. Project supported by the National Natural Science Foundation of China (No. 19872016, 19872017), the National Key Basic Research Special Foundation (G1999032805) and the Foundation for University Key Teachers by the Ministry of Education of China.     

A parallel computational method for simulating two‐phase gel dynamics on a staggered grid

We develop a parallel computational algorithm for simulating models of gel dynamics where the gel is described by two phases, a networked polymer and a fluid solvent. The models consist of transport equations for the two phases, two coupled momentum equations, and a volume‐averaged incompressibility constraint. Multigrid with Vanka‐type box‐relaxation scheme is used as preconditioner for the Krylov subspace solver (GMRES) to solve the momentum and incompressibility equations. Through numerical experiments of a model problem, the efficiency, robustness and scalability of the algorithm are illustrated. Copyright © 2008 John Wiley & Sons, Ltd.     

基于Welch法的协方差随机子空间方法的模态参数识别

对工程结构进行环境激励下的模态参数识别具有重要意义, 而随机子空间法作为适合环境激励下模态参数识别的时域方法, 由于噪声和复杂激励的原因, 会产生虚假模态、真实模态遗漏、系统自动定阶难和计算效率等问题, 这些问题阻碍了该方法在实际工程中的广泛应用. 本文提出了基于Welch法的随机子空间方法, 通过Welch法对振动响应在频域进行去噪、降低环境激励和其他不确定性因素影响的处理, 把结构固有模态从噪声和激励频率中突显出来, 形成富含更多结构模态的Toeplitz矩阵, 然后进行奇异值分解和状态矩阵计算, 最后进行特征值分析. 为了实现自动定阶, 对不同奇异值分量构建的状态矩阵得到的特征参数, 进行模糊C均值聚类分析和模态的平均相位偏移分析, 剔除虚假模态, 实现结构模态参数的自动识别. 并把本文所提出方法应用于一座大跨悬索桥的实测加速度响应分析, 和一座七十层的高层建筑的加速度响应分析, 跟频域分解法、传统随机子空间法和基于相关分析的随机子空间法的计算结果进行了比较, 发现基于Welch方法的随机子空间法相比于传统随机子空间法和基于相关分析的随机子空间法, 在避免模态遗漏和计算效率方面有显著提高, 而相对于频域分解法则在自动识别和剔除虚假模态方面有明显优势.