CN112487577A - 一种结构稳态非线性动力学响应快速计算方法及系统和存储介质 - Google Patents

一种结构稳态非线性动力学响应快速计算方法及系统和存储介质 Download PDF

Info

Publication number
CN112487577A
CN112487577A CN202011356942.5A CN202011356942A CN112487577A CN 112487577 A CN112487577 A CN 112487577A CN 202011356942 A CN202011356942 A CN 202011356942A CN 112487577 A CN112487577 A CN 112487577A
Authority
CN
China
Prior art keywords
linear
response
substructure
nonlinear
representing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011356942.5A
Other languages
English (en)
Other versions
CN112487577B (zh
Inventor
回忆
李珂
刘敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CN202011356942.5A priority Critical patent/CN112487577B/zh
Publication of CN112487577A publication Critical patent/CN112487577A/zh
Application granted granted Critical
Publication of CN112487577B publication Critical patent/CN112487577B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了的一种结构稳态非线性动力学响应快速计算方法及系统和存储介质，首先将目标结构沿线性构件与非线性构件的交界面分解为线性子结构和非线性子结构；对线性子结构和非线性子结构的边界处施加相互作用的界面力；根据界面力计算非线性子结构的响应；计算整体结构响应的解。该方法通过将目标结构分解为线性和非线性子结构并利用增量谐波平衡法形成一套计算方法，由于仅有非线性子结构需要开展非线性分析，因此计算效率较现有方法有大幅提升，该方法所得结构的稳态响应准确性高和计算效率高。

Description

一种结构稳态非线性动力学响应快速计算方法及系统和存储 介质
技术领域
本发明涉及结构动力分析技术领域，特别是一种结构稳态非线性动力学响应快速计算方法及系统和存储介质。
背景技术
非线性动力学计算方法是了解和掌握结构非线性动力学特征的基本手段。最为基础的非线性动力学计算方法是以Runge-Kutta法和中心差分法等为代表的时域算法。可是，时域算法往往仅能计算得出稳定解，却无法获得不稳定解。对于非线性系统而言，非稳定解的确定对于明确结构响应特征具有重要的意义。因此，频域算法便表现出了其优越性。常用的频域算法主要包括多尺度展开法，谐波平衡法等。
多尺度展开法是现阶段非线性动力学领域最为常用也最为经典的算法之一。其主要特点是在对系统进行求解的过程中，可以明确系统的一些非线性动力学基本特征，从而掌握系统发生某些非线性动力响应的原因。但是多尺度展开法由于其推导过程复杂，一般至多适用于2-3个自由度系统的非线性动力学分析。因此对于复杂结构，无法使用此类方法。1981年Lau和Cheung首次提出了增量谐波平衡法(IHB)，该方法将响应解的形式表达为傅里叶级数展开式，并通过伽辽金积分过程求解傅里叶级数的系数。相较多尺度展开法，该方法计算效率有了大幅的提升。
在这此之后，增量谐波平衡法在不断地发展，有学者在其计算效率方面不断的开展工作。2015年Wang和Zhu结合快速傅里叶变换(FFT)和Broyden方法大幅提高了IHB方法的计算效率。然而，尽管如此对于需要上百个自由度表示的复杂系统的非线性动力学问题，现有方法仍然很难开展分析研究。但是，对于很大部分的复杂结构而言，当其发生大幅振动时，往往只有很少部分构件进入非线性工作状态，而其余的大多数构件仍保持线性工作状态。例如，在基础下设置隔震支座的隔震结构，此类结构地震作用下，由于隔震支座的存在，使得上部结构一般仅会发生较小的振动，结构仍处于弹性状态，而隔震支座本身却因承担了大幅的振动变形，进入非线性振动状态。
针对这类情况，尽管已经明确的知道了大部分构件并不需要开展非线性分析，但由于部分构件已经进入非线性状态，现有的计算方法仍然只能通过对整个结构开展非线性动力学分析来掌握结构的动力响应，导致计算效率很低，或无法开展有效的结构计算。
发明内容
有鉴于此，本发明的目的在于提供一种结构稳态非线性动力学响应快速计算方法及系统和存储介质，该方法涉及非线性动力学算法，通过将目标结构分解为线性和非线性子结构并利用增量谐波平衡法实现快速计算方法。
为达到上述目的，本发明提供如下技术方案：
本发明提供的结构稳态非线性动力学响应快速计算方法，包括以下步骤：
获取目标结构及目标结构的初始参数；
将目标结构沿线性构件与非线性构件的交界面分解为线性子结构和非线性子结构；
根据交界面的界面力计算非线性子结构的响应；
计算整体结构响应的解。
进一步，所述线性子结构和非线性子结构的分解是按照以下公式进行的：
Figure BDA0002802864860000021
Figure BDA0002802864860000022
Fb,l＝-Fb,nl
式中，
Figure BDA0002802864860000023
Figure BDA0002802864860000024
ql和qnl分别表示为
Figure BDA0002802864860000025
和
Figure BDA0002802864860000026
进一步，所述非线性子结构的响应是利用弧长法按照以下公式对结构响应进行跟踪求解：
Figure BDA0002802864860000027
式中，
Figure BDA0002802864860000028
表示等效线性刚度；
Figure BDA0002802864860000029
表示等激励；
Figure BDA00028028648600000210
表示残差；ΔA表示响应的幅值增量；Δω表示外荷载的频率增量。
进一步，所述整体结构响应的解按照以下步骤计算：
计算线性子结构的界面力；
按照以下公式得到界面处个自由度响应的第k个谐波项：
Gk＝[a1kcos(kτ),a2kcos(kτ),…,ankcos(kτ),b1ksin(kτ),b2ksin(kτ),…,bnksin(kτ)]T,
式中：下标n代表界面处所有自由度的数量；
Gk表示界面处各自由度第k个谐波项的位移响应时程；ank表示界面处各自由度正弦响应的幅值；bnk表示界面处各自由度余弦响应的幅值；kτ表示相位角；
按照以下公式求出各谐波项qi,l的解Hk(k＝1,2,…,m)：
Figure BDA0002802864860000031
Figure BDA0002802864860000032
式中：Hk表示线性子结构内部自由度第k个谐波项的位移响应幅值；
Figure BDA0002802864860000033
Figure BDA0002802864860000034
将解Hk带入下式得线性子结构的界面力Fb,l的第k个谐波项Yk；
Figure BDA0002802864860000035
Figure BDA0002802864860000036
式中：Yk表示线性子结构的界面力Fb,l的第k个谐波项；
Figure BDA0002802864860000037
Figure BDA0002802864860000038
循环重复计算求出所有谐波项，得到作用于线性子结构上的界面力Fb,l；
通过线性子结构上的界面力Fb,l按照以下公式计算残差ΔF：
Figure BDA0002802864860000039
通过迭代方法不断减小ΔF的范数直至小于预设量ξ，则得到整体结构响应的解。
进一步，所述增量ΔF计算按照以下步骤进行：
将ΔF整理为以下增量形式：
Figure BDA0002802864860000041
求解所述增量方程得到ΔF的伽辽金积分
Figure BDA0002802864860000042
判断残差
Figure BDA0002802864860000043
值是否小于预设量ξ，如果是则得到整体结构的第i+1组响应的q(i +1),ω(i+1)。
本发明提供的结构稳态非线性动力学响应快速计算系统，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现以下步骤：
获取目标结构及目标结构的初始参数；
将目标结构沿线性构件与非线性构件的交界面分解为线性子结构和非线性子结构；
根据交界面的界面力计算非线性子结构的响应；
计算整体结构响应的解。
进一步，所述线性子结构和非线性子结构的分解是按照以下公式进行的：
Figure BDA0002802864860000044
Figure BDA0002802864860000045
Fb,l＝-Fb,nl；
式中，ql和qnl分别表示为
Figure BDA0002802864860000046
和
Figure BDA0002802864860000047
进一步，所述非线性子结构的响应是利用弧长法按照以下公式对结构响应进行跟踪求解：
Figure BDA0002802864860000048
式中，
Figure BDA0002802864860000049
表示等效线性刚度；
Figure BDA00028028648600000410
表示等激励；
Figure BDA00028028648600000411
表示残差；ΔA表示响应的幅值增量；Δω表示外荷载的频率增量。
本发明提供的一种存储介质，其上存储有计算机程序，该程序被处理器执行时实现权利要求1-5任一项所述方法的步骤。
本发明的有益效果在于：
本发明提供的结构稳态非线性动力学响应快速计算方法，是针对部分构件进入非线性工作状态的结构体系的非线性动力响应的快速计算方法。通过将目标结构分解为线性和非线性子结构并利用增量谐波平衡法形成一套计算方法，由于仅有非线性子结构需要开展非线性分析，因此计算效率较现有方法有大幅提升。
本实施例所述方法所得结构的稳态响应与Runge-Kutta法完全吻合，验证了方法的准确性。本实施例所述方法进行一次迭代所需时间为8秒钟。即所述方法计算速度为IHB方法的60倍。说明了本算法较高的计算效率。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚，本发明提供如下附图进行说明：
图1为线性子结构与非线性子结构分解示意图。
图2为带有隔震支座的平面混凝土框架结构示意图。
图3为结构的分解示意图。
图4为结构幅-频响应曲线示意图。
图5为本方法与Runge-Kutta法对比示意图。
图6为结构稳态非线性动力学响应快速计算方法流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明，以使本领域的技术人员可以更好的理解本发明并能予以实施，但所举实施例不作为对本发明的限定。
实施例1
如图6所示，本实施例提供的结构稳态非线性动力学响应快速计算方法，是针对部分构件进入非线性工作状态的结构体系的非线性动力响应的快速计算方法。通过将目标结构分解为线性和非线性子结构并利用增量谐波平衡法形成一套计算方法，由于仅有非线性子结构需要开展非线性分析，因此计算效率较现有方法有大幅提升。
本实施例提供的结构稳态非线性动力学响应快速计算方法，包括以下步骤：
步骤1：获取目标结构及目标结构的初始参数如结构的刚度，质量，阻尼，以及隔震支座的动力参数，将目标结构沿线性构件与非线性构件的交界面分解为线性子结构和非线性子结构，如果目标结构可明确哪些构件会进入非线性状态工作，则可将其在两类构件的分界面处分离为线性部分和非线性部分，通过分别在两部分的边界处施加界面力来代替其相互作用。
如图1所示，图1为线性子结构与非线性子结构分解示意图，将目标结构的动力控制方程改写为经过解耦的两个控制方程：
Figure BDA0002802864860000061
Figure BDA0002802864860000062
式中：
Ml，Cl，Kl是线性子结构的质量、阻尼和刚度矩阵；
Mnl，Cnl表示非线性子结构的质量和阻尼矩阵；FR表示非线性系统自身的恢复力；
Fl表示作用于线性系统上的外荷载的幅值；Ll表示反映界面力位置的映射矩阵；
Fnl表示作用于非线性系统上的外荷载的幅值；Lnl表示反映界面力位置的映射矩阵；
Fb,l表示不同子系统间的界面处作用于线性子系统的界面力；
Fb,nl表示不同子系统间的界面处作用于非线性子系统的界面力；
ql，
Figure BDA0002802864860000063
代表线性子结构的位移，速度以及加速度响应；
qnl，
Figure BDA0002802864860000064
代表非线性子结构的位移，速度以及加速度响应；
τ表示相位角；
本实施例中字母中带有“nl”下标的变量的意义与“l”下标的意义区别在于其对应于非线性子结构；其中，nl表示非线性子结构，l表示线性子结构，对线性子结构和非线性子结构的边界处施加相互作用的界面力；由作用力与反作用力定理可知：
Fb,l＝-Fb,nl。 (3)
将式(1-2)中的ql和qnl分别表示为
Figure BDA0002802864860000065
和
Figure BDA0002802864860000066
其中，qi,l表示线性子结构的内部节点位移向量；qb,l表示线性子结构的界面节点位移向量；
qi,nl表示非线性子结构的内部节点位移向量；qb,nl表示非线性子结构的界面节点位移向量；
下标i和b分别位于子结构内部和界面处的自由度。
步骤2：初步计算非线性子结构的响应
由于本实施例利用弧长法对结构响应进行跟踪求解，因此在寻找新的一组结构响应时，须已知一组已有的结构响应。
设q(i),ω(i)为目前已知的结构的第i组响应，现寻找结构的第i+1组响应。根据式(1)及q(i),ω(i)可得结构的截面力Fb,l；将Fb,l带入式(2)，将解的形式表示为增量的形式：
q＝q(i)+Δq，ω＝ω(i)+Δω (4)
其中，
q(i)表示第i组已知的结构位移响应；ω(i)表示第i组已知的结构位移响应所对应的外荷载激励频率；
q表示待求的结构位移响应；Δq表示位移响应的增量；
ω表示待求的结构位移响应所对应的外荷载激励频率；Δω表示激励频率的增量；
将式(4)带入式(2)可得：
Figure BDA0002802864860000071
Figure BDA0002802864860000072
式中，ω0＝ω(i)表示第i组已知的结构位移响应所对应的外荷载激励频率；Mnl表示非线性子结构的质量矩阵；Cnl表示非线性子结构的阻尼矩阵；Δq表示激励频率的增量；R为计算残差；F表示作用于非线性子结构的外荷载；FR(q(i))表示非线性子结构的恢复力；
τ＝ωt，变量右上方的撇代表对τ求导；
令响应表示为傅里叶级数的形式，如下所示：
q(i)＝SA,Δq＝SΔA (7)
式中，
Figure BDA0002802864860000073
其中，Ai是非线性子结构第i个自由度所有频率响应成分的幅值所组成的向量；
p表示非线性子结构自由度的数量；
Cs＝[1 cos(τ) cos(2τ) … cos((m-1)τ) sin(τ) sin(2τ) … sin((m-1)τ)](8)
式中，
m代表响应中所考虑的频率成分的数量；
将式(8)带入式(5)并对所有变量在τ在一个周期内[0,2π]进行积分，得到下式：
Figure BDA0002802864860000081
式中：
Figure BDA0002802864860000082
Figure BDA0002802864860000083
Figure BDA0002802864860000084
Figure BDA0002802864860000085
Figure BDA0002802864860000086
表示质量矩阵的伽辽金积分；
Figure BDA0002802864860000087
表示阻尼矩阵的伽辽金积分；
Figure BDA0002802864860000088
表示刚度矩阵的伽辽金积分；
Figure BDA0002802864860000089
表示残差；ST表示S矩阵的转置
对式(9)开展弧长法求解，便可初步获得非线性子结构第i+1组响应的“暂定解”。
步骤3：计算整体结构响应的解
尽管在计算时暂时将线性子结构与非线性子结构分离，但实际上在界面处具有完全相同的响应，即
qb,l＝qb,nl. (14)
将式(1)改写为：
Figure BDA00028028648600000810
Figure BDA00028028648600000811
式中：
Figure BDA00028028648600000812
Mii，Mib表，Mbi，Mbb分别表示Ml的四个子矩阵；
Cii，Cib，Cbi，Cbb分别表示Cl的四个子矩阵；
Kii，Kib，Kbi，Kbb分别表示Kl的四个子矩阵；
Figure BDA00028028648600000813
表示线性子结构内部各自由度的加速度响应；
Figure BDA00028028648600000814
表示线性子结构界面各自由度的加速度响应；
Figure BDA00028028648600000815
表示线性子结构内部各自由度的速度响应；
Figure BDA00028028648600000816
表示线性子结构界面各自由度的速度响应；qi,l表示线性子结构内部各自由度的位移响应；qb,l表示线性子结构界面各自由度的位移响应；Fii表示作用于线性子结构内部各自由度的外荷载；Fbb表示作用于线性子结构界面各自由度的外荷载。
由式(14)，可知Fb,l(qb,l)＝Fb,l(qb,nl)。根据式(7)所示响应解的形式，可进一步计算线性子结构的界面力。
将式(15-16)中所有的qb,l由q'b,nl替换，界面处个自由度响应的第k个谐波项可表示为
Gk＝[a1kcos(kτ),a2kcos(kτ),…,ankcos(kτ),b1ksin(kτ),b2ksin(kτ),…,bnksin(kτ)]T, (17)
式中：n代表界面处所有自由度的数量；
式中：下标n代表界面处所有自由度的数量；
Gk表示界面处各自由度第k个谐波项的位移响应时程；ank表示界面处各自由度正弦响应的幅值；bnk表示界面处各自由度余弦响应的幅值；kτ表示相位角；
由此，根据式(15)可分别求出qi,l各谐波项的解Hk(k＝1,2,…,m)，如下式所示：
Figure BDA0002802864860000091
式中：Hk表示线性子结构内部自由度第k个谐波项的位移响应幅值；
Figure BDA0002802864860000092
Figure BDA0002802864860000093
Q表示线性子结构内部自由度的等效刚度；P表示线性子结构界面自由度的等效刚度；
将Hk带入式(16)，可得线性子结构的界面力Fb,l的第k个谐波项Yk，如下所示
Figure BDA0002802864860000094
式中：Yk表示线性子结构的界面力Fb,l的第k个谐波项；
Figure BDA0002802864860000095
Figure BDA0002802864860000096
V表示线性子结构内部自由度的计算刚度；U表示线性子结构界面自由度的计算刚度；
重复式(18-23)将所有谐波项均求出后，便可获得Fb,l。
比较分别有线性和非线性子系统所得到的界面力
ΔFb＝Fb,l(qb,l)+Fb,nl(qnl) (24)
式中，
ΔFb表示线性子结构与非线性子结构界面力的残差；
Fb,l(qb,l)表示线性子结构的界面力；
Fb,nl(qnl)表示非线性子结构的界面力；
由式(1)可知
Figure BDA0002802864860000101
将式(25)带入式(24)可得
Figure BDA0002802864860000102
理论上只有当ΔF＝0时，才意味着两个子系统的响应完全协调，即系统得到了真正的响应。因此，需要通过迭代的方法不断减小ΔF的范数直至该值小于某一小量ξ。
将式(26)整理为增量形式，并进行伽辽金积分
Figure BDA0002802864860000103
Figure BDA0002802864860000104
Figure BDA0002802864860000105
Figure BDA0002802864860000106
Figure BDA0002802864860000107
ΔAj＝[01×(j-1),δaj,01×(m-j)]T, (32)
式中：
Figure BDA0002802864860000108
表示非线性子结构的质量矩阵的伽辽金积分；
Figure BDA0002802864860000109
表示非线性子结构的阻尼矩阵的伽辽金积分；
Figure BDA00028028648600001010
表示非线性子结构的刚度矩阵的伽辽金积分；
Figure BDA00028028648600001011
表示残差；
Figure BDA00028028648600001012
表示线性子结构对非线性子结构所提供的等效刚度；
Figure BDA00028028648600001013
表示线性子结构界面处的等效刚度；
Figure BDA0002802864860000111
表示线性子结构界面处的等效刚度矩阵的第j列；Ab表示界面位移向量；
ΔAj表示假想界面位移的增量；01×(m-j)表示0向量；δaj为一接近于0的正数；
式(27)可再次借助弧长法求解，直至
Figure BDA0002802864860000112
小于ξ，则所得结果即为整体结构的第i+1组响应的q(i+1),ω(i+1)。
实施例2
下面通过具体的实例对本发明以及本发明效果进行进一步的说明。
如图2所示，图2为带有各镇支座的平面混凝土框架结构示意图，以某五层三开间的钢筋混凝土结构的一榀框架的非线性动力计算为例，其柱底安装隔震支座；混凝土材料的密度和弹性模量分别按2.5×103kg/m3 and 3.0×104Mpa计算。框架结构的层高为3m，开间尺寸为6m。柱子和梁的横截面尺寸分别设为0.45m×0.45m和0.3m×0.6m。
本例中每层的柱沿高度被划分为2个单元，同时每跨的梁被划分为3个单元。对上部结构所有构件均采用伯努利梁单元进行模拟。因此，上部结构共由85个单元组成，包括74个节点。对于隔震支座，则选用Bouc-Wen作为描述其本构关系的数学模型，如式(33-34)所示
fB＝k1x+k2z, (33)
Figure BDA0002802864860000113
式中：
fB为隔震支座所提供的反力；
x是隔震支座的剪切变形；
z为归一化的滞回力；
k1支座恢复力线性分量的刚度；
k2恢复力非线性分量的模量系数；
β和γ为控制滞回环尺度和形状的参数；
r为反应滞回曲线光滑度的正整数；
k Bl＝k1+k2被称为支座的线性化参数。
本实施例中设β＝γ＝50。同时令kBl＝1.0×103kN/m(其中k1＝100kN/m,k2＝900kN/m)，即支座的线性化刚度为上部柱子刚度的1/45。
由图2可知，当结构承受基底加速度激励后，其上部结构由于隔震支座的存在并不会发生大幅振动，即处于线性工作状态，而各镇支座则会承担大部分结构的位移，导致其非线性工作。由此，以柱底为界面将结构分解为线性子结构和非线性子结构两部分，如图3所示，其中，图3中(3a)为线性子结构；(3b)非线性子结构，图3为结构的分解。
由于每一个隔震支座仅有一个水平自由度，因此非线性子结构共包含4个自由度，而线性子结构包含214个自由度。
设结构承受幅值为0.33m/s2的单频(谐波)基底加速度激励，激励频率由0.25Hz逐渐变化至1.5Hz(覆盖了结构的第一阶自振频率)。现计算结构在每一个激励频率下的稳态响应。
则结构的动力控制方程可写为
Figure BDA0002802864860000121
Figure BDA0002802864860000122
式中：Mf表示上部结构的质量矩阵；Cf表示上部结构的阻尼矩阵；Kf表示上部结构的刚度矩阵；
Figure BDA0002802864860000123
表示隔震支座提供的恢复力；
Figure BDA0002802864860000124
Figure BDA0002802864860000125
Figure BDA0002802864860000126
Figure BDA0002802864860000127
Figure BDA0002802864860000128
式中，FN表示；
k1A，k1B，k1C，k1D表示支座A、B、C、D的恢复力线性分量的刚度；
k2A，k2B，k2C，k2D表示支座A、B、C、D的恢复力非线性分量的模量系数；
βA，βB，βC，βD，γA，γB，γC，γD表示支座A、B、C、D的控制滞回环尺度和形状的参数；
Figure BDA0002802864860000131
表示第1-4自由度的速度响应；
q5，q6，q7，q8表示第5-8自由度的位移响应；
I为单位矩阵；k1X＝100kN/m，k2X＝900kN/m，其中，X＝A,B,C,D；
依据式(13)：
Figure BDA0002802864860000132
式中：
Figure BDA0002802864860000133
Figure BDA0002802864860000134
Figure BDA0002802864860000135
Cs＝[1 cos(τ) cos(2τ) … cos((m-1)τ) sin(τ) sin(2τ) … sin((m-1)τ)],
C′s＝[0 -sin(τ) -2sin(2τ) … -(m-1)sin((m-1)τ) cos(τ) 2cos(2τ) … (m-1)cos((m-1)τ)].
式中，Cs表示位移响应的基本谐波成分向量；F′N1表示支座非线性恢复力对
Figure BDA0002802864860000136
求导雅可比矩阵；F′N2表示支座非线性恢复力对q5，q6，q7，q8求导的雅可比矩阵；m表示谐波频率成分的个数；
Figure BDA0002802864860000137
表示非线性子结构的位移响应幅值向量；W表示非线性子结构的恢复力矩阵。
本实施例中将m设为12。依据本实施例中所述方法可得结构屋顶及柱底的幅-频响应曲线如图4所示，图4中的(4a)为第1谐波项，(4b)为第3谐波项。图4表示结构幅频响应曲线。为验证算法的准确性，将本实施例所述算法与经典的时程算法Runge-Kutta法进行对比。选取激励频率为0.77Hz(0.67ω1)和0.32Hz，分别通过Runge-Kutta法计算结构响应。所得屋顶及柱底的时程响应曲线，如图5所示，图5为所提方法与Runge-Kutta法对比，其中，(5a)为屋顶响应(ω＝4.83rad/s)，(5b)为柱底响应(ω＝4.83rad/s)，(5c)为屋顶响应(ω＝2.01rad/s)，(5d)为柱底响应(ω＝2.01rad/s)，可以看出，本实施例所述方法所得结构的稳态响应与Runge-Kutta法完全吻合，验证了方法的准确性。
本实施例所述方法与经典的IHB方法对本实施例中的结构开展非线性动力学计算。发现如果令m＝12，即计算结果相对更加准确的情况下，采用IHB方法时由于计算
Figure BDA0002802864860000141
所需时间过长，基本无法看展计算并进行对比。因此，选取m＝6对两种算法进行对比。IHB方法计算一次迭代所需时间为8分钟左右，而本实施例所述方法进行一次迭代所需时间为8秒钟。即所述方法计算速度为IHB方法的60倍。说明了本算法较高的计算效率。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例，本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换，均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。

Claims (9)

1.结构稳态非线性动力学响应快速计算方法，其特征在于：包括以下步骤：
获取目标结构及目标结构的初始参数；
将目标结构沿线性构件与非线性构件的交界面分解为线性子结构和非线性子结构；
根据交界面的界面力计算非线性子结构的响应；
计算整体结构响应的解。
2.如权利要求1所述的方法，其特征在于：所述线性子结构和非线性子结构的分解是按照以下公式进行的：
Figure FDA0002802864850000011
Figure FDA0002802864850000012
Fb,l＝-Fb,nl
式中，
Figure FDA0002802864850000013
Figure FDA0002802864850000014
ql和qnl分别表示为
Figure FDA0002802864850000015
和
Figure FDA0002802864850000016
3.如权利要求1所述的方法，其特征在于：所述非线性子结构的响应是利用弧长法按照以下公式对结构响应进行跟踪求解：
Figure FDA0002802864850000017
式中，
Figure FDA0002802864850000018
表示等效线性刚度；
Figure FDA0002802864850000019
表示等激励；
Figure FDA00028028648500000110
表示残差；ΔA表示响应的幅值增量；Δω表示外荷载的频率增量。
4.如权利要求1所述的方法，其特征在于：所述整体结构响应的解按照以下步骤计算：
计算线性子结构的界面力；
按照以下公式得到界面处个自由度响应的第k个谐波项：
Gk＝[a1kcos(kτ),a2kcos(kτ),…,ankcos(kτ),b1ksin(kτ),b2ksin(kτ),…,bnksin(kτ)]T,
式中：下标n代表界面处所有自由度的数量；
Gk表示界面处各自由度第k个谐波项的位移响应时程；ank表示界面处各自由度正弦响应的幅值；bnk表示界面处各自由度余弦响应的幅值；kτ表示相位角；
按照以下公式求出各谐波项qi,l的解Hk(k＝1,2,…,m)：
Figure FDA0002802864850000021
Figure FDA0002802864850000022
式中：Hk表示线性子结构内部自由度第k个谐波项的位移响应幅值；
Figure FDA0002802864850000023
Figure FDA0002802864850000024
将解Hk带入下式得线性子结构的界面力Fb,l的第k个谐波项Yk；
Figure FDA0002802864850000025
Figure FDA0002802864850000026
式中：Yk表示线性子结构的界面力Fb,l的第k个谐波项；
Figure FDA0002802864850000027
Figure FDA0002802864850000028
循环重复计算求出所有谐波项，得到作用于线性子结构上的界面力Fb,l；
通过线性子结构上的界面力Fb,l按照以下公式计算残差ΔF：
Figure FDA0002802864850000029
通过迭代方法不断减小ΔF的范数直至小于预设量ξ，则得到整体结构响应的解。
5.如权利要求4所述的方法，其特征在于：所述增量ΔF计算按照以下步骤进行：
将ΔF整理为以下增量形式：
Figure FDA00028028648500000210
求解所述增量方程得到ΔF的伽辽金积分
Figure FDA0002802864850000031
判断残差
Figure FDA0002802864850000032
值是否小于预设量ξ，如果是则得到整体结构的第i+1组响应的q(i+1),ω(i+1)。
6.结构稳态非线性动力学响应快速计算系统，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：
获取目标结构及目标结构的初始参数；
将目标结构沿线性构件与非线性构件的交界面分解为线性子结构和非线性子结构；
根据交界面的界面力计算非线性子结构的响应；
计算整体结构响应的解。
7.如权利要求6所述的系统，其特征在于：所述线性子结构和非线性子结构的分解是按照以下公式进行的：
Figure FDA0002802864850000033
Figure FDA0002802864850000034
Fb,l＝-Fb,nl；
式中，ql和qnl分别表示为
Figure FDA0002802864850000035
和
Figure FDA0002802864850000036
8.如权利要求6所述的系统，其特征在于：所述非线性子结构的响应是利用弧长法按照以下公式对结构响应进行跟踪求解：
Figure FDA0002802864850000037
式中，
Figure FDA0002802864850000038
表示等效线性刚度；
Figure FDA0002802864850000039
表示等激励；
Figure FDA00028028648500000310
表示残差；ΔA表示响应的幅值增量；Δω表示外荷载的频率增量。
9.一种存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1-5任一项所述方法的步骤。
CN202011356942.5A 2020-11-27 2020-11-27 一种结构稳态非线性动力学响应快速计算方法及系统和存储介质 Active CN112487577B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011356942.5A CN112487577B (zh) 2020-11-27 2020-11-27 一种结构稳态非线性动力学响应快速计算方法及系统和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011356942.5A CN112487577B (zh) 2020-11-27 2020-11-27 一种结构稳态非线性动力学响应快速计算方法及系统和存储介质

Publications (2)

Publication Number Publication Date
CN112487577A true CN112487577A (zh) 2021-03-12
CN112487577B CN112487577B (zh) 2023-06-02

Family

ID=74936144

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011356942.5A Active CN112487577B (zh) 2020-11-27 2020-11-27 一种结构稳态非线性动力学响应快速计算方法及系统和存储介质

Country Status (1)

Country Link
CN (1) CN112487577B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117236142A (zh) * 2023-11-13 2023-12-15 江苏新筑预应力工程有限公司 解算索杆系的有限质点分析方法、系统及介质

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050256613A1 (en) * 2004-05-14 2005-11-17 Massachusetts Institute Of Technology Active control vibration isolation using dynamic manifold
CN103853885A (zh) * 2014-02-24 2014-06-11 昆明理工大学 一种多壁碳纳米管非线性振动特性预测方法
CN105046008A (zh) * 2015-07-30 2015-11-11 南京航空航天大学 分段非线性隔振器及其设计方法
CN107609221A (zh) * 2017-08-15 2018-01-19 东南大学 一种基于遗传算法的含铰结构非线性参数识别方法
CN107643028A (zh) * 2017-09-15 2018-01-30 中建隧道建设有限公司 浅埋铁路隧道下穿民房软弱段爆破施工方法
CN108334699A (zh) * 2018-02-05 2018-07-27 东北大学 一种改进的单侧激振非线性卸载系统下散体物料流动特性分析方法
CN109408894A (zh) * 2018-09-26 2019-03-01 西安交通大学 一种考虑阻尼结构摩碰的透平机械叶片非线性振动特性分析方法
CN110059286A (zh) * 2019-03-07 2019-07-26 重庆大学 一种基于fft的结构非平稳响应高效分析方法
CN110083920A (zh) * 2019-04-23 2019-08-02 西南交通大学 一种地震作用下非比例阻尼结构随机响应的分析方法
CN110162826A (zh) * 2019-03-20 2019-08-23 北京机电工程研究所 薄壁结构热气动弹性动响应分析方法
CN110569585A (zh) * 2019-08-28 2019-12-13 东北大学 一种螺栓连接圆柱壳结构动态特性计算与分析方法
CN110795884A (zh) * 2019-11-20 2020-02-14 黑龙江科技大学 一种基于多尺度模型更新的新型混合试验方法
CN111428409A (zh) * 2020-03-24 2020-07-17 西安交通大学 用于机械非线性振动分析的非线性动力学方程求解方法及系统
CN111832200A (zh) * 2020-06-04 2020-10-27 台州学院 一种附加干摩擦阻尼器的循环对称结构频响分析方法

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050256613A1 (en) * 2004-05-14 2005-11-17 Massachusetts Institute Of Technology Active control vibration isolation using dynamic manifold
CN103853885A (zh) * 2014-02-24 2014-06-11 昆明理工大学 一种多壁碳纳米管非线性振动特性预测方法
CN105046008A (zh) * 2015-07-30 2015-11-11 南京航空航天大学 分段非线性隔振器及其设计方法
CN107609221A (zh) * 2017-08-15 2018-01-19 东南大学 一种基于遗传算法的含铰结构非线性参数识别方法
CN107643028A (zh) * 2017-09-15 2018-01-30 中建隧道建设有限公司 浅埋铁路隧道下穿民房软弱段爆破施工方法
CN108334699A (zh) * 2018-02-05 2018-07-27 东北大学 一种改进的单侧激振非线性卸载系统下散体物料流动特性分析方法
CN109408894A (zh) * 2018-09-26 2019-03-01 西安交通大学 一种考虑阻尼结构摩碰的透平机械叶片非线性振动特性分析方法
CN110059286A (zh) * 2019-03-07 2019-07-26 重庆大学 一种基于fft的结构非平稳响应高效分析方法
CN110162826A (zh) * 2019-03-20 2019-08-23 北京机电工程研究所 薄壁结构热气动弹性动响应分析方法
CN110083920A (zh) * 2019-04-23 2019-08-02 西南交通大学 一种地震作用下非比例阻尼结构随机响应的分析方法
CN110569585A (zh) * 2019-08-28 2019-12-13 东北大学 一种螺栓连接圆柱壳结构动态特性计算与分析方法
CN110795884A (zh) * 2019-11-20 2020-02-14 黑龙江科技大学 一种基于多尺度模型更新的新型混合试验方法
CN111428409A (zh) * 2020-03-24 2020-07-17 西安交通大学 用于机械非线性振动分析的非线性动力学方程求解方法及系统
CN111832200A (zh) * 2020-06-04 2020-10-27 台州学院 一种附加干摩擦阻尼器的循环对称结构频响分析方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
GIUSEPPE REGA: "Nonlinear vibrations of suspended cables—Part I:Modeling and analysis", 《APPLIED MECHANICS REVIEWS》, vol. 57, no. 6, pages 444 - 464 *
WEI LU ET AL.: "Nonlinear dynamics of a submerged floating moored structure by incremental harmonic balance method with FFT", 《MARINE STRUCTURES》, vol. 31, pages 65 - 76 *
YI HUI ET AL.: "Extended IHB method for dynamic analysis of structures with geometrical and material nonlinearities", 《ENGINEERING STRUCTURES》, vol. 205, pages 2 - 3 *
倪迎鸽等: "增量谐波平衡法在分段结构非线性气动弹性系统的求解", 《科学技术与工程》, vol. 19, no. 3, pages 248 - 252 *
姚红良等: "多频激励局部非线性系统响应求解的降维增量谐波平衡法", 《振动工程学报》, vol. 28, no. 5, pages 742 - 746 *
尉飞等: "一种求解局部非线性结构稳态响应的方法", 《应用力学学报》, vol. 26, no. 4, pages 793 - 796 *
王陶等: "针对局部非线性问题的混合坐标模态综合法", 《航空学报》, vol. 37, no. 09, pages 2758 - 2763 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117236142A (zh) * 2023-11-13 2023-12-15 江苏新筑预应力工程有限公司 解算索杆系的有限质点分析方法、系统及介质
CN117236142B (zh) * 2023-11-13 2024-01-26 江苏新筑预应力工程有限公司 解算索杆系的有限质点分析方法、系统及介质

Also Published As

Publication number Publication date
CN112487577B (zh) 2023-06-02

Similar Documents

Publication Publication Date Title
Mahin et al. Pseudodynamic test method—Current status and future directions
Darby et al. Real-time substructure tests using hydraulic actuator
Lee et al. Real-time substructuring technique for the shaking table test of upper substructures
Jensen et al. An enhanced substructure coupling technique for dynamic re-analyses: Application to simulation-based problems
CN112487577A (zh) 一种结构稳态非线性动力学响应快速计算方法及系统和存储介质
Richiedei et al. Unit-rank output feedback control for antiresonance assignment in lightweight systems
Belostotsky et al. Strength and stability analysis of load-bearing structures of Evolution Tower with allowance for actual positions of reinforced concrete structural members
Wang et al. Hybrid testing with model updating on steel panel damper substructures using a multi‐axial testing system
Palomba et al. Reduced-order observers for nonlinear state estimation in flexible multibody systems
Rafezy et al. Coupled lateral–torsional frequencies of asymmetric, three-dimensional frame structures
Wang Sturm-Liouville equation for free vibration of a tube-in-tube tall building
Calofir et al. Dynamic nonlinear modelling of building structure using the force analogy method
Archer A technique for the reduction of dynamic degrees of freedom
Mueller Real-time hybrid simulation with online model updating
Samikkannu et al. System identification-based aeroelastic modelling for wing flutter
Dongfang et al. An improved equivalent beam model of large periodic beam-like space truss structures
He et al. An accurate method for the dynamic behavior of tensegrity structures
Chang An improved on-line dynamic testing method
Schellenberg et al. Developemnt of a Large-Scale 6Dof Hybrid Shake Table and Application To Testing Response Modification Devices for Tall Buildings
Kirsanov et al. The lower limit of the first frequency of natural vibrations externally statically indeterminate truss: analytical solution
Du et al. Anti-seismic Performance Comparison of Response Spectrum Analysis and Time History Analysis Based on Computer Big Data
Ramachandran et al. Family of columns isospectral to gravity-loaded columns with tip force: A discrete approach
Shmerling Matrix equations models for nonlinear dynamic analysis of two-dimensional and three-dimensional RC structures with lateral load resisting cantilever elements
Tsopelas et al. Nonlinear dynamic analysis of multi-base seismically isolated structures with uplift potential I: formulation
Flashner An approach to modeling vibrations of systems composed of beams, rigid bodies, and point masses

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant