显式数值积分算法的改进和性能研究

来源 :大连理工大学 | 被引量 : 0次 | 上传用户:superdai
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
数值积分算法是解决结构动力学在时间上离散的运动微分方程的有效方法。大量学者已经提出了各种显式积分算法,这些算法具有无条件稳定性和可控的数值耗散的理想特性。但是,涉及积分算法数值漂移特性的研究是有限的。本文基于控制理论,基于极点映射法运用高精度的预修正双线性变换将连续域的传递函数映射到离散域,开发了一族新的结构相关的显式积分算法,称为TL-φ算法。与现有算法相比,该方法的优点在于它可以通过与结构的关键频率相关的附加参数来控制数值漂移程度。本文还充分研究了所提出的算法在线性和非线性系统中的稳定性,数值耗散和数值漂移特性。结果表明,所提出的TL-φ算法在求解线性系统结构动力响应是无条件稳定的,而对于非线性系统则是有条件稳定的。TL-φ算法的数值耗散特性与其他显式算法的数值耗散特性非常相似,但与其他方法相比,它具有最小化周期误差的能力,这有利于求解系统的结构动力响应,特别是对于那些具有高固有频率的系统。本文通过四个数值算例研究了该方法的性能,结果表明该算法可以更好地解决复杂的线性和非线性结构动力问题。此外本文还将此算法的设计方法推广到其他显式积分算法如CR算法,并提出了数值性能更优的CR-φ算法,验证了该设计方法的泛用性。最后结合实例说明了算法相关参数-关键频率的选取,为该算法在工程和试验中的应用提供了便利。
其他文献
多重线性系统在科学计算与工程领域中有广泛的应用.近年来涌现出不少数值算法来解系数张量为非奇异Μ-张量、右端项为正向量的多重线性系统.求解算法主要分为两类:第一类是张量分裂算法,第二类是其他的张量算法.由于线性系统中已经出现了一类将Anderson加速法用到古典迭代的算法,且具有求解时间往往更短的优势,考虑到张量是矩阵的高阶推广,本文将Anderson加速法用于已有的张量分裂算法中,试图加速已有张量
元古界基岩潜山油气藏是柴达木盆地北缘马海东地区油气勘探的重要领域之一,储层是影响基岩储集体油气富集的主要因素之一。根据野外露头、岩心、薄片、扫描电镜和成像测井等资料,对基岩储层的发育特征及控制因素进行研究。结果表明:马海东地区元古界基岩储层发育片岩、片麻岩、榴辉岩和角闪岩等4大类9种岩性,裂缝和溶蚀孔隙等2大类6种储集空间类型,储层物性以低孔-(特)低渗透为主;基岩裂缝是影响溶蚀孔隙、溶蚀裂缝发育
不可分解表示是群表示论研究的一个重要方向,Higman给出了有限群在代数闭域下的不可分解模表示个数是否有限的充要条件[6].Janusz则进一步给出了具体不可分解表示的计数公式[7].在本文的讨论中,假定所涉及的基础域均是特征数为7的代数闭域,而主要的研究对象是一类特殊的pq阶群,即21阶非交换群{x3=1,y7=1,xy=xy2}.我们首先根据Janusz给出的计数公式计算了该有限群的不等价不可
随机矩阵集中不等式是随机矩阵理论中的一个重要研究方向,它研究的问题是限定随机矩阵和的极端特征值(或范数)大于给定常数的概率,被广泛地应用于机器学习、压缩感知、量子计算、优化等领域。国内对这一问题的研究还处于起步阶段,因此本文总结了当前主流的随机矩阵集中不等式的框架。这些框架的主要步骤如下:首先,使用矩阵的拉普拉斯变换方法,通过矩阵矩母函数的迹控制随机矩阵和的特征值的尾部概率;然后,使用Golden
六角系统(即苯系统)是指一个没有割点有限连通平面图,它的每个内面都是一个单位正六边形.六角系统在化学领域中应用十分广泛,对原子轨道、能量等问题有重要意义.在化学图论中,将苯类碳氢化合物的氢原子忽略不计后,苯类碳氢化合物的原子骨架图即对应一个六角系统.六隅体多项式的研究对求解Clar数以及完美匹配数是至关重要的,还有助于了解一些其他的拓扑性质.本文主要研究了几类六角系统中六隅体多项式的解析性质,具体
实际生活中,每一特定结果与其关联特征的因果关系往往不可能在短期内形成,通常都会存在一定的时间滞后效应。例如:上市公司违约风险的预测结果主要依赖于可观察期甚至延迟期内多类别多窗口解释变量的变化。本文将无约束分布滞后模型(DLM)与常见的预测模型,如:逻辑回归模型与支持向量机模型相结合,提出一种可应对延迟滞后影响的预警系统新框架。然而,在实际操作过程中需利用大量高维的原始数据进行预测模型的学习训练,且
本文主要研究中心极限定理在确定性动力系统下的收敛速度.我们首先证明了vn=∑?J=0n-1v°Tj满足中心极限定理,其中是均值为0的H(?)lder连续的观察函数,是非一致扩张映射.之后,我们得到了vn在Kolmogorov度量下的收敛速度.除了非一致扩张系统,我们的结果还适用于非一致双曲系统.
作为海洋工程中常见的结构之一,圆柱在水流中所受的水动力荷载及其内在机理一直是工程界关注的问题。当一定流速的来流作用于圆柱时,圆柱的尾部会产生交替脱落的旋涡。交替脱落的旋涡将对圆柱施加一定频率的周期性荷载,根据共振的原理,当周期性荷载与圆柱的固有频率接近时圆柱会发生共振,严重时甚至会破坏工程结构。因此,圆柱绕流问题一直以来都是海洋工程的研究重点。稳定流和振荡流是圆柱绕流问题中最基本的两种流动形式。本
两参数特征值问题是数值代数方向中一个重要的研究问题,在物理学、气象学、航天、结构动力学等领域也经常出现,它主要有两类求解方法.第一类是将两参数特征值问题转化为广义特征值问题,比如最小残差商迭代(MRQI)法和QZ算法等用于求解所有特征值,还有Jacobi-Davidson法用于求解接近给定研究值的特征值.第二类是同伦法,求解所有特征值.围道积分法由于具有很好的并行性和求解给定区域内特征值的优点,因
本文主要研究了一类二阶常微分方程在Neumann边界条件下的分歧现象.通过构造全连续线性算子,联系线性化方程的本征值,证明分歧点的存在性.进而研究分歧分支的结构问题,建立该方程的全局分歧定理及单边全局分歧定理.最后解决了9种情况下该方程结点解的存在性问题.