非线性结构动力学系统瞬态响应无条件稳定时间积分方法技术方案

本发明专利技术公开的非线性结构动力学系统瞬态响应无条件稳定时间积分方法，属于非线性结构动力学系统响应分析领域。本发明专利技术实现方法如下：建立包含非线性几何项、非线性阻尼项以及它们的耦合项的结构动力学方程；结合多分步技术将时间单元离散成三个分步，并设计各分步状态变量更新公式；基于BN稳定性理论、传递因子分析方法、局部截断误差分析方法，以效率最大化、二阶精度、可控耗散和无条件稳定性为目标设计分步公式中的参数；建立三分步时间积分方法在包含非线性几何项、非线性阻尼项以及它们的耦合项在结构动力学系统中的求解流程，实现对包含非线性几何项、非线性阻尼项以及它们的耦合项的结构动力学系统瞬态响应的高效、精确且稳定地时域分析。且稳定地时域分析。且稳定地时域分析。

【技术实现步骤摘要】
非线性结构动力学系统瞬态响应无条件稳定时间积分方法


[0001]本专利技术属于非线性结构动力学系统响应分析领域，尤其涉及一种包含非线性几何项、非线性阻尼项及其耦合项的结构动力学系统的具有二阶精度、可调数值耗散以及无条件稳定性的时间积分方法。
[0002]背景方法
[0003]随着新技术、新材料、新结构等前沿成果在航空航天、土木工程、机器人等领域的广泛应用，由此带来的材料非线性、几何非线性和接触非线性因素给结构动力学分析带来了严峻挑战。目前，最为常用的结构动力学分析方案是“空间离散”与“时间离散”相结合，其也被广泛应用在商业软件CAE(Computer Aided Engineering)中。
[0004]灵活有效的有限单元法(Finite Element Method,FEM)是目前较为流行的空间离散工具，其作用是将由偏微分方程控制的无穷自由度连续系统转变为由常微分方程控制的有限自由度离散系统。完成空间离散后，下一个环节就是时间离散，基于差分思想建立的时间积分方法(Time integration method,TIM)是一种强大的时间离散工具。其求解思路是将待求的时间域离散成一系列连续的时间区间，人为假设状态变量在一个时间区间内的变化规律，从上一时刻已知状态变量信息求解当前时刻未知信息。
[0005]Newmark家族方法是时间积分方法中的经典工作，其成员包括著名的梯形法则(Trapezoidal Rule,TR)和中心差分方法(Central Difference Method,CDM)等。对于线性系统，梯形法则是无条件稳定的，但对于简单的非线性问题如刚性单摆转动问题，梯形法则出现了发散现象。由此引发当前时间积分方法领域的热点和难点问题：如何提高时间积分方法在求解非线性结构动力学系统时的稳定性，以及能否为非线性结构动力学系统设计出无条件稳定的时间积分方法。为了解决上述问题，学者们发展了耗散型方法、保能量方法以及结构依赖型方法。
[0006]大量的数值结果表明通过引入数值耗散可以有效地推迟时间积分方法在非线性结构动力学系统求解中失稳的时刻。由于Newmark方法的耗散格式仅具有一阶精度，因此学者们通过加权Newmark方法的平衡方程构造出了同时具有二阶精度和数值耗散的参数类方法，比如广义
‑
alpha方法、Wood
‑
Bossak
‑
Zienkiewicz
‑
alpha方法、Hiber
‑
Hughes
‑
Taylor
‑
alpha方法等。但是，由于平衡方程在时间结点上没有得到满足，导致该类方法的加速度只具有一阶精度。为了设计出具有较高低频精度的耗散型时间积分方法，线性多步法和复合方法被提出，这两类方法均有效地提高了精度。但是，线性多步法和复合方法并不能从根本上改变方法的稳定性。总之，耗散型方法仍然存在失稳问题，并且其低频精度比非耗散型方法低。
[0007]为了能够给非线性结构动力学系统提供稳定的数值仿真结果，基于能量有界准则的保能量方法被提出。能量有界准则由Belytschko提出，其定义是：当前时刻的动能和势能之和小于等于上一时刻的动能、势能和外力功之和。因此，保能量方法是通过强制系统能量守恒来实现无条件稳定性。受构造方法的限制，现有的保能量方法大都只能处理含有非线性几何项的结构动力学问题，比如Krenk方法和EMM等，而只有少部分如ECM能够处理非线性
几何项和阻尼项独立存在的非线性结构动力学问题。目前，能够处理含有非线性几何项和阻尼项耦合的结构动力学系统的保能量方法尚未被构造出来。综上，保能量方法虽然在求解非线性结构动力学系统时不会发散，但其适用对象有限，且具有额外计算量(如离散能量函数修正运算等)致使计算效率降低。
[0008]对大型复杂结构进行动力学响应分析时，精度、效率和稳定性三者之间往往难以平衡。为了解决这个问题，2002年Chang提出了结构依赖型时间积分方法，这类方法的算法参数与结构初始特性和时间步长大小紧密相关。为了进一步提高Chang方法的数值性能，包括精度、稳定性和耗散特性，一系列性能更加优异的结构依赖型时间积分方法被提出，比如KR方法、CR方法、Du
‑
Yang
‑
Zhao方法等。理论分析和数值实验结果表明，该类方法对线性系统是无条件稳定的，同时对非线性软刚度系统具有理想的稳定性。目前，受构造方式限制，这类方法局限于具有对角质量矩阵、正定对称刚度和阻尼矩阵的线性系统和具有幂形式非线性刚度的结构动力学系统。总之，结构依赖型方法虽然有效地平衡了精度、效率和稳定性，但其在求解非线性结构动力学系统时仍有失稳的可能，且适用对象有限。
[0009]因此，针对包含非线性几何项、非线性阻尼项及其耦合项的结构动力学系统，设计可高精度、高效率、无条件稳定地分析瞬态响应的时间积分方法是现阶段亟待突破的研究内容。其将为非线性结构动力学问题提供性能优越的新解法，且可以促进时间积分方法朝着同时改进稳定性、精度和效率方向发展，服务于国家对复杂结构动力学分析高性能算法的需求，具有重要的理论意义和应用价值。

技术实现思路

[0010]本专利技术公开的非线性结构动力学系统瞬态响应无条件稳定时间积分方法主要目的是：利用有限单元方法建立包含非线性几何项、非线性阻尼项以及它们的耦合项的结构动力学方程；结合多分步技术将一个时间单元离散成三个分步，并设计各分步状态变量(包括位移、速度和加速度)的更新公式；基于BN稳定性理论、传递因子分析方法、局部截断误差分析方法，以效率最大化、二阶精度、可控耗散和无条件稳定性为目标设计分步公式中的参数；在此基础上，建立三分步时间积分方法在包含非线性几何项、非线性阻尼项以及它们的耦合项在结构动力学系统中的求解流程，实现对包含非线性几何项、非线性阻尼项以及它们的耦合项的结构动力学系统瞬态响应的高效、精确且稳定地时域分析。
[0011]本专利技术的目的是通过如下技术方法实现。
[0012]本专利技术公开的非线性结构动力学系统瞬态响应无条件稳定时间积分方法，通过利用有限单元方法建立包含非线性几何项、非线性阻尼项以及它们耦合项的结构动力学方程；将一个时间单元划分成三个分步，基于差分格式建立状态变量在一个时间单元内的演化规律，依次使用n点后向插值格式(n＝1,2,3)，建立便于数值性能调控的多参数优化时间步进方程，由于调控参数增多，从而提高参数调节的可控性和精度。所述数值性能包括效率、精度、耗散、稳定性。基于构建的多参数优化时间步进方程，推导便于数值性能调控的三个分步的有效刚度矩阵，以效率最大化为目标建立具有统一系数的参数约束关系，利用建立参数约束关系将所述三个分步的有效刚度矩阵变换为具有统一形式的有效刚度矩阵，提高非线性系统的结构动力学响应时域分析效率。基于构建的多参数优化时间步进方程，针对包含非线性几何项、非线性阻尼项以及它们的耦合项的结构动力学系统，以二阶精度为
目标建立局部截断误差为O(Δt3)的参本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.非线性结构动力学系统瞬态响应无条件稳定时间积分方法，其特征在于：包括以下步骤，步骤1：基于有限单元方法FEM建立空间离散的包含非线性几何项、非线性阻尼项及其耦合项的结构动力学控制方程；步骤2：基于差分格式建立状态变量在一个时间单元内的演化规律，将一个时间单元划分成三个分步，依次使用n点后向插值格式(n＝1,2,3)，建立便于数值性能调控的多参数优化时间步进方程，由于调控参数增多，提高参数调节的可控性和精度；所述数值性能包括效率、精度、耗散、稳定性；步骤3：基于步骤2构建的多参数优化时间步进方程，推导便于数值性能调控的三个分步的有效刚度矩阵，以效率最大化为目标建立具有统一系数的参数约束关系，利用建立参数约束关系将所述三个分步的有效刚度矩阵变换为具有统一形式的有效刚度矩阵，提高非线性系统的结构动力学响应时域分析效率；根据具有统一形式的有效刚度矩阵得到统一系数的参数约束关系；步骤4：基于步骤2构建的多参数优化时间步进方程，针对包含非线性几何项、非线性阻尼项以及它们耦合项的结构动力学系统，以二阶精度为目标建立局部截断误差为O(Δt3)的参数约束关系，以数值耗散可控为目标建立高频段传递因子可调的参数约束关系，以无条件稳定性为目标建立满足BN稳定性条件的参数约束关系，从而提高非线性系统的结构动力学响应时域分析的精度、耗散性和稳定性；步骤5：根据步骤3得到具有统一系数的参数约束关系，根据步骤4得到满足二阶精度、可控数值耗散和无条件稳定性的参数约束关系，将步骤3和步骤4得到的参数约束关系代入步骤2建立的多参数优化时间步进方程，将所述多参数优化时间步进方程结合Newton迭代技术，代入步骤1构建的包含非线性几何项、非线性阻尼项及其耦合项的结构动力学控制方程，实现高效、精确、且无条件稳定地计算包含非线性几何项、非线性阻尼项及其耦合项的结构动力学系统瞬态响应。2.如权利要求1所述的非线性结构动力学系统瞬态响应无条件稳定时间积分方法，其特征在于：步骤1实现方法为，基于有限单元方法FEM建立空间离散的包含非线性几何项、非线性阻尼项及其耦合项的结构动力学控制方程，如下所示：式中，M为质量矩阵；为包含非线性内力和外部激励项的向量；t为时间；对于线性系统，向量退化为其中K为刚度矩阵；C为阻尼矩阵；Q(t)为外部激励向量。3.如权利要求2所述的非线性结构动力学系统瞬态响应无条件稳定时间积分方法，其特征在于：步骤2实现方法为，基于差分格式，建立状态变量，包括位移、速度和加速度，在一个是时间单元[t,t+Δt]内随时间t的演化规律；将一个时间单元[t,t+Δt]划分成三个分步，即[t,t+c1Δt]、[t+c1Δt,t+c2Δt]以及[t+c2Δt,t+c3Δt]；第一个分步[t,t+c1Δt]的步进公式为
其中，x
t+c1Δt
、和代表t+c1Δt时刻的位移、速度和加速度；x
t
和代表t时刻的位移和速度；Δt代表时间步长大小；c1为第一个分步引入的算法参数；第二个分步[t+c1Δt,t+c2Δt]的步进公式为其中，x
t+c2Δt
、和代表t+c2Δt时刻的位移、速度和加速度；c2和α为第二个分步引入的算法参数；第三个分步[t+c2Δt,t+c3Δt]的步进公式为其中，x
t+c3Δt
、和代表t+c3Δt时刻的位移、速度和加速度；c3、β和η为第三个分步引入的算法参数；时间单元终点时刻t+Δt信息通过对三个已知时刻信息t+c1Δt、t+c2Δt以及t+c3Δt进行加权得到，对应的步进方程为其中，b1和b2为待定加权参数。4.如权利要求3所述的非线性结构动力学系统瞬态响应无条件稳定时间积分方法，其特征在于：步骤3实现方法为，当时，整理三个分步的步进方程得其中，和分别代表有效刚度矩阵和有效载荷向量；三个分步的有效刚度矩阵的表达式为
从公式(7)至(9)中发现，当参数满足公式(10)中的关系时，三个分步的有效刚度矩阵完全相同；此时，三个分步使用同一个有效刚度矩阵，进而降低计算量；c1＝c2α＝c3η
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)另外，要求时间点t+c1Δt和t+c3Δt关于时间单元中点...

【专利技术属性】
技术研发人员：季奕，田强，张欢，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：