结构动力方程的一种自适应步长增维精细积分法

1

作者 黄宇熙 崔颖 杨国刚 《振动与冲击》 EI CSCD 北大核心 2023年第14期198-203,236,共7页

增维精细积分法是一种求解结构动力方程的高精度逐步积分算法,其步长的选取会对计算精度产生极大的影响,在实际应用中存在难以确定合适步长的问题。为满足实际工程中对计算精度和效率的要求,提出了一种计算误差的估计方法,并以估计误差... 增维精细积分法是一种求解结构动力方程的高精度逐步积分算法,其步长的选取会对计算精度产生极大的影响,在实际应用中存在难以确定合适步长的问题。为满足实际工程中对计算精度和效率的要求,提出了一种计算误差的估计方法,并以估计误差和迭代收敛速度为基准,建立了一种自适应步长增维精细积分法。针对三种结构动力方程的算例结果表明,在计算各类线性及非线性振动问题时,该方法均可以在保证计算精度的前提下快速有效地控制算法的计算步长,并且仅需较少的额外计算消耗,显著提高了增维精细积分法的计算效率,使该方法在求解结构动力方程时更具计算优势和实用价值。 展开更多

关键词 增维精细积分法 步长自适应 结构动力方程 误差估计

下载PDF 职称材料

结构动力方程Newmark-β方法递推简化分析 被引量：10

2

作者 王海波 陈伯望 余志武 《四川大学学报（工程科学版）》 EI CAS CSCD 北大核心 2008年第3期47-52,共6页

以Newm ark-β方法为基础,基于叠加原理推导出了一种新的结构动力分析数值积分递推格式,在计算位移时,不需要计算速度和加速度等中间值,因而在进行结构弹塑性地震反应分析时更为简单、方便。对该递推格式的起步条件及稳定性进行了分析... 以Newm ark-β方法为基础,基于叠加原理推导出了一种新的结构动力分析数值积分递推格式,在计算位移时,不需要计算速度和加速度等中间值,因而在进行结构弹塑性地震反应分析时更为简单、方便。对该递推格式的起步条件及稳定性进行了分析。通过算例分析可以得出,当β从1/4到1/8变化时,随β的减小,计算精度越高,但稳定性越差。 展开更多

关键词 结构动力方程 Newmark-β方法 数值积分

下载PDF 职称材料

结构动力方程求解的改进精细Runge-Kutta方法 被引量：5

3

作者 张继锋 邓子辰 张凯 《应用数学和力学》 CSCD 北大核心 2015年第4期378-385,共8页

在已有精细Runge-Kutta(龙格-库塔)方法的基础上,考虑了状态空间方程非齐次项的特点和外荷载的特殊性,提出了求解结构动力方程的改进精细Runge-Kutta方法.通过对矩阵进行分块计算,在利用原有精细Runge-Kutta方法高精度的同时进一步提高... 在已有精细Runge-Kutta(龙格-库塔)方法的基础上,考虑了状态空间方程非齐次项的特点和外荷载的特殊性,提出了求解结构动力方程的改进精细Runge-Kutta方法.通过对矩阵进行分块计算,在利用原有精细Runge-Kutta方法高精度的同时进一步提高了计算效率,有利于大型结构的长时间仿真.将改进精细Runge-Kutta方法应用于结构动力方程求解,为其求解提供一种新方法.数值算例表明了改进方法的正确性和有效性. 展开更多

关键词 结构动力方程 精细积分 简化计算 RUNGE-KUTTA方法

下载PDF 职称材料

结构动力方程的精细积分-FFT方法 被引量：2

4

作者 张文志 富明慧 刘祚秋 《中山大学学报（自然科学版）》 CAS CSCD 北大核心 2008年第6期12-15,共4页

离散结构动力方程为差分方程,并假设始、末时刻位移已知,将初值问题形式上转化为边值问题,然后利用快速傅立叶变换(FFT)进行求解,从而得到非齐次结构动力方程的一个数值特解。将该数值特解与通过精细积分法求得的齐次方程通解相结合,建... 离散结构动力方程为差分方程,并假设始、末时刻位移已知,将初值问题形式上转化为边值问题,然后利用快速傅立叶变换(FFT)进行求解,从而得到非齐次结构动力方程的一个数值特解。将该数值特解与通过精细积分法求得的齐次方程通解相结合,建立了求解结构动力方程的一种新方法。该方法具有较高的精度和计算效率,算例的数值结果证明了本文方法的有效性。 展开更多

关键词 结构动力方程 精细积分法 差分方程 快速傅立叶变换

下载PDF 职称材料

定常结构动力方程增维精细积分法求解的注记 被引量：8

5

作者 任伟 杜铁钧 《杭州电子科技大学学报（自然科学版）》 2005年第1期41-43,共3页

精细积分法应用于求解线性定常结构的动力方程,能得到逼近于精确解的数值解,而增维精细积分法在求解非齐次的动力方程时,通过化非齐次方程为齐次方程来求解,避免了矩阵求逆的困难。该方法的提出,使得在用增维精细积分法求解定常结构动... 精细积分法应用于求解线性定常结构的动力方程,能得到逼近于精确解的数值解,而增维精细积分法在求解非齐次的动力方程时,通过化非齐次方程为齐次方程来求解,避免了矩阵求逆的困难。该方法的提出,使得在用增维精细积分法求解定常结构动力方程时,避免每一个时间步都做20次矩阵加、乘迭代。该方法对程序的实现十分有利,而且对大型问题的长时间仿真可明显提高计算效率,算例显示了该方法的有效性。 展开更多

关键词 精细积分法 增维方法 定常结构动力方程 矩阵求逆 齐次方程

下载PDF 职称材料

结构动力方程的显式级数积分格式

6

作者 郭泽英 李青宁 《地震工程与工程振动》 CSCD 北大核心 2009年第4期53-57,共5页

将Newmark-β法中常平均加速度法的基本假定与精细指数算法结合,根据指数矩阵的Taylor级数展开式,提出了动力方程的显式级数解,并设计了相应的时程积分算法。该算法的精度可根据Taylor级数展开式的项数进行灵活控制。算例的结果表明:在... 将Newmark-β法中常平均加速度法的基本假定与精细指数算法结合,根据指数矩阵的Taylor级数展开式,提出了动力方程的显式级数解,并设计了相应的时程积分算法。该算法的精度可根据Taylor级数展开式的项数进行灵活控制。算例的结果表明:在满足稳定性条件的前提下,随着时间步长的增加,其精度优于传统的时程积分法。通过稳定性的分析,指出其稳定性条件是显然满足的。 展开更多

关键词 结构动力方程 常平均加速度法 精细算法 TAYLOR级数 稳定性

下载PDF 职称材料

一种结构动力方程的并行直接积分方法

7

作者 胡宁 《工程力学》 EI CSCD 1992年第1期65-71,共7页

本文提出了一种求解大型结构动力方程的新的并行直接积分方法。该方法在L.Brusa和L.Nigro提出的一步(one-step)直接积分方法的基础上,引进并行运算步骤。并行运算步骤是通过将动力积分方程子结构化,同时进行组集和凝聚实现的。该方法在... 本文提出了一种求解大型结构动力方程的新的并行直接积分方法。该方法在L.Brusa和L.Nigro提出的一步(one-step)直接积分方法的基础上,引进并行运算步骤。并行运算步骤是通过将动力积分方程子结构化,同时进行组集和凝聚实现的。该方法在西安交通大学ELXSI-6400并行机上程序实现,计算结果表明能有效地求解大型结构有限元动力方程的并行直接积分方法。 展开更多

关键词 并行算法 直接积分 结构动力方程

下载PDF 职称材料

多自由度结构动力方程解析解的改进算法 被引量：1

8

作者 任泽民 凌霄霄 +3 位作者 乔金丽 丁冬 孟秋杰 吕超 《河北工业大学学报》 CAS 2018年第4期68-74,共7页

针对结构动力计算中数值方法不精确以及模态叠加法计算复杂,对结构质量、刚度、阻尼矩阵形式有限制等一系列缺点,提出了多自由度体系运动方程的改进解析算法,推导了改进算法下自由振动时方程的通解形式,并对简谐荷载以及插值形式的地震... 针对结构动力计算中数值方法不精确以及模态叠加法计算复杂,对结构质量、刚度、阻尼矩阵形式有限制等一系列缺点,提出了多自由度体系运动方程的改进解析算法,推导了改进算法下自由振动时方程的通解形式,并对简谐荷载以及插值形式的地震荷载作用下多自由度体系的运动方程进行了解析求解.文中在求解时对结构的质量、刚度、阻尼矩阵的形式及状态矩阵的正交性不做要求.计算结果与文献中所给结果非常接近,同时将所得解析式代入运动方程时,方程左右两边的误差在10-15量级.数例计算和证明推导过程体现了解的精确性.此解为多自由度体系运动方程解析解的求解方法提供了一个新的思路. 展开更多

关键词 结构动力方程 多自由度体系 解析解 地震作用

下载PDF 职称材料

结构体系动力方程求解的显式积分格式的能耗特性 被引量：11

9

作者 李小军 唐晖 《工程力学》 EI CSCD 北大核心 2007年第2期28-33,共6页

针对作者提出的结构体系动力方程求解的一种显式积分格式,探讨其数值计算能耗特性即算法阻尼特性,导出了其算法阻尼值随体系的物理阻尼和结构体系振动频率值变化的关系,并进一步讨论了该积分格式用于无限介质波动的数值模拟中抑制或消... 针对作者提出的结构体系动力方程求解的一种显式积分格式,探讨其数值计算能耗特性即算法阻尼特性,导出了其算法阻尼值随体系的物理阻尼和结构体系振动频率值变化的关系,并进一步讨论了该积分格式用于无限介质波动的数值模拟中抑制或消除透射边界引起的计算高频失稳问题,给出了该积分格式与中心差分格式的对比分析算例。 展开更多

关键词 结构工程 结构动力方程 显式积分格式 数值计算能耗 透射边界 高频失稳

下载PDF 职称材料

结构动力方程另类显式逐步积分格式研究 被引量：2

10

作者 张晓志 刘龙江 曹玲 《世界科技研究与发展》 CSCD 2011年第6期978-982,共5页

现有显式积分格式的稳定性与计算结果的取信度均远逊于常用的隐式积分格式,这是制约现有显式积分格式在大型结构动力分析软件以及其它相关软件中广泛应用的瓶颈。以解决这一瓶颈问题为研究目标,通过结合采用精确法和平均常加速度法的基... 现有显式积分格式的稳定性与计算结果的取信度均远逊于常用的隐式积分格式,这是制约现有显式积分格式在大型结构动力分析软件以及其它相关软件中广泛应用的瓶颈。以解决这一瓶颈问题为研究目标,通过结合采用精确法和平均常加速度法的基本假定,推导建立了基于双重基本假定的另类显式积分格式,提出了相关的体系另类振动周期等新概念。基于另类振动周期概念并结合现有相关研究成果,导出了该格式的精度控制准则和形式上的稳定性控制准则,为开展相关的数值模拟研究准备了必要的基础和前提。概述了开展相关数值模拟研究的技术路线,对业已取得的主要数值模拟研究结果做了先期预告和初步讨论。 展开更多

关键词 结构动力方程 双重基本假定 另类显式积分格式 体系另类振动周期

原文传递

Runge-Kutta方法求解结构动力学方程 被引量：8

11

作者 吴志桥 高普云 任钧国 《系统仿真学报》 CAS CSCD 北大核心 2010年第9期2085-2090,2105,共7页

将几种具有不同稳定性的Runge-Kutta方法应用到结构动力学方程的数值求解中。针对增量形式的动力学方程,使用改进的Newton-Raphson迭代,研究了减少计算量的两种方法:(1)使用单对角隐式Runge-Kutta方法,(2)应用转化矩阵。采用逼近算子的... 将几种具有不同稳定性的Runge-Kutta方法应用到结构动力学方程的数值求解中。针对增量形式的动力学方程,使用改进的Newton-Raphson迭代,研究了减少计算量的两种方法:(1)使用单对角隐式Runge-Kutta方法,(2)应用转化矩阵。采用逼近算子的谱半径分析了稳定性与数值阻尼特性,解释了L-稳定方法抑制高频振荡的原因。数值算例表明在精确解上较小的物理阻尼能有效的抑制高频振荡,但对各种直接积分方法的影响很小,高精度的L-稳定Runge-Kutta方法能在有效抑制高频振荡的同时高精度的求解低频振动。 展开更多

关键词 结构动力学方程 RUNGE-KUTTA方法 数值阻尼 L-稳定性

原文传递

结构动力学方程的显式积分格式 被引量：4

12

作者 刘恒 廖振鹏 《地震工程与工程振动》 CSCD 北大核心 2009年第1期32-43,共12页

本文从空间解耦有限元常微分方程组出发,探讨了结构动力学方程的高精度显式积分格式。通过被积函数的拉格朗日多项式内插和分部积分导出了波动数值模拟的一组显式时步积分公式。这组公式是时间和空间解耦的,即波场内任一离散节点在任一... 本文从空间解耦有限元常微分方程组出发,探讨了结构动力学方程的高精度显式积分格式。通过被积函数的拉格朗日多项式内插和分部积分导出了波动数值模拟的一组显式时步积分公式。这组公式是时间和空间解耦的,即波场内任一离散节点在任一时刻的波动数据可以用这组公式依据该节点及其邻近节点在该时刻之前的n+1个时刻的波动数据显式地算出(n为非负整数),阐明了这组公式的如下特点:第一,其截断误差的量级不超过O(Δtn+3),Δt为时间步距。第二,它不仅可用于线性波动的数值模拟,而且可用于本构方程具有强非线性情形。第三,这组公式也可推广应用于一系列数学物理暂态问题的数值求解。针对一个简单的时不变系统初步分析了此组积分格式的稳定性。但是,对其稳定性尚需作进一步研究。 展开更多

关键词 波动 数值模拟 结构动力学方程 显式积分 有限元 非线性

下载PDF 职称材料

非线性结构动力学方程的迭代时程积分方法 被引量：1

13

作者 凌道盛 吴强 《非线性动力学学报》 1995年第4期290-297,共8页

本文把文（１）提出的精细时程积分法推广到非线性动力学方程的响应分析中。时间步预估一迭代的方法减少了非线性迭代的步数，改善了其收敛性。计算的结果表明，该方法具有较高的精度。

关键词 结构动力学方程 迭代法 时程积分法

下载PDF 职称材料

结构动力学方程的数值解法研究 被引量：1

14

作者 韩爱红 钱晓军 《低温建筑技术》 2015年第8期81-83,共3页

本文介绍了求解结构动力学方程的数值积分方法,主要包括:线性加速度法、Wilson-θ法、Newmark-β法、中心差分法、Houbolt法和精细积分法;论述了各种数值积分方法的基本原理、稳定性和适用范围;通过算例,指出了各种数值积分方法的优缺点... 本文介绍了求解结构动力学方程的数值积分方法,主要包括:线性加速度法、Wilson-θ法、Newmark-β法、中心差分法、Houbolt法和精细积分法;论述了各种数值积分方法的基本原理、稳定性和适用范围;通过算例,指出了各种数值积分方法的优缺点,证实了精细积分法在计算精度和在长步长、长时间内中保持稳定性的优越性。 展开更多

关键词 结构动力学方程 数值积分方法 精细积分法

下载PDF 职称材料

超大型海洋浮式结构物水弹性响应预报的研究现状和发展方向 被引量：24

15

作者 崔维成 《船舶力学》 EI 2002年第1期73-90,共18页

超大型海洋浮式结构物在未来的海洋资源开发和海洋空间利用方面有重要应用前景 ,但它在波浪中的水弹性响应是该结构物设计过程中所遇到的一个重要技术问题 ,国际上在这方面已做了大量的研究工作。本文比较全面地介绍了国内外关于VLFS在... 超大型海洋浮式结构物在未来的海洋资源开发和海洋空间利用方面有重要应用前景 ,但它在波浪中的水弹性响应是该结构物设计过程中所遇到的一个重要技术问题 ,国际上在这方面已做了大量的研究工作。本文比较全面地介绍了国内外关于VLFS在波浪中的水弹性响应的研究现状 。 展开更多

关键词 超大型海洋浮式结构物 水弹性响应 非均匀海洋环境 预报 现状 研究进展 速度势 VLFS 结构动力学方程 格林函数法

下载PDF 职称材料

非比例阻尼结构参数识别算法的研究 被引量：6

16

作者 王建有 陈键云 林皋 《振动与冲击》 EI CSCD 北大核心 2005年第3期1-3,18,i001,共5页

在地下结构的地震响应中,无限地基的辐射阻尼起着重要的影响,然而其阻尼形式同结构的比例阻尼假定不同。同样,对于施加隔振、减震措施的建筑结构,其系统阻尼也不再满足比例阻尼假定,在结构参数识别中将阻尼阵假定为与刚度阵形式相似的... 在地下结构的地震响应中,无限地基的辐射阻尼起着重要的影响,然而其阻尼形式同结构的比例阻尼假定不同。同样,对于施加隔振、减震措施的建筑结构,其系统阻尼也不再满足比例阻尼假定,在结构参数识别中将阻尼阵假定为与刚度阵形式相似的方法不能反映实际结构参数。从研究结构的阻尼模型发发,引入了节点堆积阻尼,这样形成的阻尼与刚度、质量既有一定的联系,又有一定的独立性。以剪切型结构为例对此做了详细说明,结果表明,引入堆积阻尼后,进行参数识别的结构动力学方程是一个非线性识别问题,在此基础上提出了基于最小二乘平均的非线性迭代识别算法。 展开更多

关键词 识别算法 非比例阻尼 结构动力学方程 结构参数识别 剪切型结构 非线性迭代 地震响应 地下结构 辐射阻尼 无限地基 建筑结构 减震措施 系统阻尼 阻尼模型 识别问题 最小二乘 刚度阵 阻尼阵 独立性 堆积 节点 基础

下载PDF 职称材料

基于三次样条插值的精细积分法 被引量：10

17

作者 高小科 邓子辰 黄永安 《振动与冲击》 EI CSCD 北大核心 2007年第9期75-77,82,共4页

结合指数矩阵的精细算法,提出了一类基于三次样条插值的精细积分方法。针对结构动力学方程一般解中的积分项,考虑在一个时间步长内激励为线性和正余弦两种变化形式,通过对积分项中的指数矩阵进行三次样条插值函数模拟,得到一组新的被积... 结合指数矩阵的精细算法,提出了一类基于三次样条插值的精细积分方法。针对结构动力学方程一般解中的积分项,考虑在一个时间步长内激励为线性和正余弦两种变化形式,通过对积分项中的指数矩阵进行三次样条插值函数模拟,得到一组新的被积函数,最后通过多次分部积分,构造了一类新的高精度计算格式。在三次样条插值函数构造过程中引入了指数矩阵的精细算法,有效避免了中间过程中有效数字的丢失,同时还有效解决了HPD-F算法中涉及的矩阵求逆问题,大大增加了算法的数值稳定性。数值算例显示了该方法的有效性。 展开更多

关键词 结构动力方程 时程积分 样条插值 指数矩阵 精细积分法

下载PDF 职称材料

一种新的改进精细直接积分法 被引量：3

18

作者 张继锋 邓子辰 +1 位作者 徐方暖 张凯 《动力学与控制学报》 2015年第4期241-245,共5页

针对结构动力方程转化为状态空间方程后矩阵维数增加而导致计算量增大的问题,考虑状态空间方程中所含外部荷载的特点,提出了一种新的改进精细直接积分法.给出了利用梯形公式、复化梯形公式、辛普生公式、复化辛普生公式、科特斯公式、... 针对结构动力方程转化为状态空间方程后矩阵维数增加而导致计算量增大的问题,考虑状态空间方程中所含外部荷载的特点,提出了一种新的改进精细直接积分法.给出了利用梯形公式、复化梯形公式、辛普生公式、复化辛普生公式、科特斯公式、高斯公式计算杜哈姆积分时的计算格式,分析了不同计算格式下的计算精度和计算效率.数值算例表明本文改进方法的正确性. 展开更多

关键词 结构动力方程 直接积分 分块计算 精细积分 改进方法

下载PDF 职称材料

一个基于双重基本假定的另类显式积分格式

19

作者 张晓志 刘龙江 曹玲 《东北电力大学学报》 2011年第5期98-103,共6页

通过结合采用精确法和平均常加速度法的基本假定,推导建立了基于双重基本假定的另类显式积分格式,提出了相关的体系另类振动周期等新概念。基于另类振动周期概念并结合现有相关研究成果,导出了该格式的精度控制准则和形式上的稳定性控... 通过结合采用精确法和平均常加速度法的基本假定,推导建立了基于双重基本假定的另类显式积分格式,提出了相关的体系另类振动周期等新概念。基于另类振动周期概念并结合现有相关研究成果,导出了该格式的精度控制准则和形式上的稳定性控制准则,为开展相关的数值模拟研究准备了必要的基础和前提。 展开更多

关键词 结构动力方程 双重基本假定 另类显式积分格式 体系另类振动周期

下载PDF 职称材料

结构非线性动力分析的精细积分多步法 被引量：5

20

作者 陈伯望 王海波 《工程力学》 EI CSCD 北大核心 2009年第5期41-46,共6页

将目前常用的非线性动力状态方程v=H·v+f(v,t)变换为v=H·v+f(v,t)+r(t),其中H·v、f(v,t)和r(t)分别是右端项的线性齐次部分、非线性部分和非齐次荷载项。将精细积分法和预估-校正Adams-Bashforth-Moulton多步法相结合,... 将目前常用的非线性动力状态方程v=H·v+f(v,t)变换为v=H·v+f(v,t)+r(t),其中H·v、f(v,t)和r(t)分别是右端项的线性齐次部分、非线性部分和非齐次荷载项。将精细积分法和预估-校正Adams-Bashforth-Moulton多步法相结合,对非线性动力方程进行求解。数值算例表明:该方法的稳定性和计算精度明显优于现有的Adams-Bashforth-Moulton方法,可用于多自由度结构体系的非线性地震反应分析。 展开更多

关键词 结构非线性动力方程 Adams-Bashforth-Moulton多步法 精细积分法 预估-校正 数值积分

原文传递