CN105701279B - 非经典结构动力响应频域方法 - Google Patents

非经典结构动力响应频域方法 Download PDF

Info

Publication number
CN105701279B
CN105701279B CN201610008351.6A CN201610008351A CN105701279B CN 105701279 B CN105701279 B CN 105701279B CN 201610008351 A CN201610008351 A CN 201610008351A CN 105701279 B CN105701279 B CN 105701279B
Authority
CN
China
Prior art keywords
dynamic response
classical
frequency domain
domain
structural dynamic
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.)
Active
Application number
CN201610008351.6A
Other languages
English (en)
Other versions
CN105701279A (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.)
Ocean University of China
Original Assignee
Ocean University of China
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 Ocean University of China filed Critical Ocean University of China
Priority to CN201610008351.6A priority Critical patent/CN105701279B/zh
Publication of CN105701279A publication Critical patent/CN105701279A/zh
Application granted granted Critical
Publication of CN105701279B publication Critical patent/CN105701279B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Bridges Or Land Bridges (AREA)

Abstract

本发明涉及一种非经典结构动力响应频域方法，特别是针对柔性大跨径桥梁结构、海洋平台、海上风电等结构的动力响应分析问题。包括以下步骤：S1、构建拉普拉斯域非经典结构传递函数；S2、非经典结构传递函数系数获取；S3、非经典结构传递函数极值、零值获取；S4、非经典结构动力响应频域分析。本发明中通过构建非经典结构传递函数的对应形式，提出基于传递函数的分析方法，使对其系数的获取不再依赖于模态空间解耦，拓展了传统频域动力响应方法的适用范围，克服传统频域方法无法处理非传统结构的不足，工程上也为大型复杂结构动力响应分析提供一种更为通用、高效的频域方法，工程应用前景广阔。

Description

非经典结构动力响应频域方法
技术领域
本发明涉及一种非经典结构动力响应频域方法，特别是针对柔性大跨径桥梁结构、海洋平台、海上风电等结构的动力响应分析问题。
背景技术
近年来各国越来越重视海洋资源的开发与利用，相继建造了许多诸如跨海大桥、海洋平台、海底油气管道、海底隧道、超大型浮体等海洋工程结构。这类结构容易在波、流、地震等外部荷载的作用下产生有害振动，严重时可能产生疲劳破坏、而且其施工难度大，养护维修费用高，因此，研究结构的动力响应十分重要。
柔性大跨径桥梁模型风洞试验时，因风的作用而导致的等效刚度矩阵和等效阻尼矩阵的非对称性已是公认的事实；尤其对海洋工程结构而言，其阻尼通常包括结构本身阻尼、流体附加阻尼，甚至土壤附加阻尼，导致真实阻尼矩阵形式不确定。当结构的质量矩阵、刚度矩阵、阻尼矩阵三者之中任一矩阵非对称时，该结构动力系统将变为非经典系统。
目前，对于结构动力响应分析，国内外很早就已经开展了大量的研究，但成果主要集中于发展相对比较成熟的时域、频域或频-时域变换方法(如SEAS AM等商业软件)：时域计算方法因其计算精度高且能考虑非线性等因素而被认为是目前最准确的动力响应计算方法，但该类方法因需通过计算卷积导致计算效率相对较低；频域计算方法通过传递函数将输入荷载谱(FFT)与输出反应谱联系起来，其结果仅为结构的稳态响应而无法得到结构的瞬态响应(即初始条件为零)，对服役中海洋结构修正模型而言，因其初始条件不可能完全为零值而导致无法直接应用频域方法进行求解；频-时域变换方法综合了时域、频域方法二者的优势，虽然在一定程度上兼顾了计算精度与效率，但理论上仍然依赖于FFT变换，即周期性谐波假定不符合实际环境条件，同时，无瞬态解的不足也无法在理论上得到解决。
总之，在处理非对称质量矩阵、刚度矩阵、阻尼矩阵问题时，传统时域方法计算效率相对较低；传统频域方法虽然计算效率较高，但其因结构振动方程无法在模态空间进行解耦而不再适用，为此，亟待提出一种计算精度高、计算效率优的非经典结构动力响应方法。
发明内容
本发明提出一种非经典结构动力响应频域方法，克服传统频域动力响应方法无法处理非对称矩阵的问题；在工程上为大型复杂结构动力响应分析提供一种更为通用的频域方法，确保工程分析效率及精度。
为了达到上述目的，本发明提出非经典结构动力响应频域方法，包括以下步骤：
S1、构建拉普拉斯域非经典结构传递函数；
S2、非经典结构传递函数系数获取；
S3、非经典结构传递函数极值、零值获取；
S4、非经典结构动力响应频域分析。
进一步的，所述步骤S1包括以下步骤：
S11、构建结构在第i个自由度上的响应：Fj(s)记为第j个荷载分量，N为荷载数目，D(s)为非经典结构阻抗的行列式，Di,j(s)为将D(s)第j列替换为Fj(s)的代数余子式；
S12、构建非经典结构传递函数：其中κ为系数，且其中η与τ分别为极值与零值；可以回避在模态空间求解振型问题，工程上能够处理更通用的阻尼矩阵形式。
进一步的，所述步骤S2中，传递函数系数的获得采用的方法。
进一步的，所述步骤S2中，传递函数系数的获得采用的方法，其中M为结构质量矩阵，Mi,j为M去掉第j行、第i列。非经典结构传递函数系数是对应于步骤S12的形式，目前不存在相同的求解方法；本发明提出两种非经典结构传递函数系数求解方法，便于本发明应用者在工程中相互验证。
进一步的，所述步骤S4中非经典结构动力响应频域分析包括：
S41、构建非经典结构动力响应频域表达式其中s＝i2πf；
S42、对xi,j(f)进行傅里叶逆变换，获得非经典结构动力响应时域表达式；
S43、多荷载输入、多位置输出的动力响应计算：固定输出位置i、循环位置j，并将每一次循环结果进行线性叠加，即可得到单输入、多输出时域动力响应时程；
S44、在单输入、多输出基础上，循环输出位置i，并根据行数为对应自由度的原则，获得结构所有位置处的时域动力响应时程。扩展了传统频域方法的理论适用范围，在工程上又保留了计算效率高的技术优势。
与现有技术相比，本发明的优点和积极效果在于：
本发明中通过构建非经典结构传递函数的对应形式，提出基于传递函数的分析方法，使对其系数的获取不再依赖于模态空间解耦，拓展了传统频域动力响应方法的适用范围，工程上又保留了分析效率高的优势，应用前景广阔。
附图说明
图1为实施例中框架结构模型参考图；
图2为实施例中施加外荷载参考图；
图3为比例阻尼下本发明与传统方法传递函数对比：(a)幅值，(b)相位；
图4为比例阻尼下本发明与传统方法响应对比；
图5为非经典结构应用本发明与传统方法传递函数对比：(a)幅值，(b)相位；
图6为非经典结构应用本发明与传统方法响应对比。
具体实施方式
本发明提供一种非经典结构动力响应频域方法，克服传统频域动力响应方法无法处理非对称质量矩阵、刚度矩阵、阻尼矩阵的问题，工程上也为大型复杂结构动力响应分析提供一种更为通用、高效的频域方法，下面结合具体的实施例对本发明做进一步地说明。
本发明提出非经典结构动力响应频域方法，包括以下步骤：
S1、构建拉普拉斯域非经典结构传递函数；
S2、非经典结构传递函数系数获取；
S3、非经典结构传递函数极值、零值获取；
S4、非经典结构动力响应频域分析。
由于传统经典结构传递函数的求解需将结构运动方程转换到模态空间以求解其振型，但当结构质量、刚度、阻尼矩阵三者任一矩阵为非对称时，即非经典结构运动方程在模态空间无法解耦，导致现有技术无法直接求解非经典结构传递矩阵，针对此缺陷，所述步骤S1构建针对拉普拉斯域的非经典结构传递函数，包括以下步骤：
S11、构建结构在第i个自由度上的响应：Fj(s)记为第j个荷载分量，N为荷载数目，D(s)为非经典结构阻抗的行列式，Di,j(s)为将D(s)第j列替换为Fj(s)的代数余子式；
S12、构建非经典结构传递函数：其中κ为系数，且其中η与τ分别为极值与零值。
通过构建基于步骤S11的非经典结构传递函数，并提出基于步骤S12的传递函数分析方法，可以有效回避在模态空间求解振型问题，在工程上也能够处理更通用的阻尼矩阵形式。
实施例中，所述步骤S2中，传递函数系数的获得可以采用的方法，亦可以采用的方法，其中M为结构质量矩阵，Mi,j为M去掉第j行、第i列。非经典结构传递函数系数是对应于步骤S12的形式，本实施例提出的两种获得系数的方法，目前不存在相同的方式。且此两种方式的思想来源于结构质量矩阵在数学上为正定的工程实际，便于本发明应用者在工程中相互验证，提高准确度。
获得系统的极值和零值后，对非经典结构动力响应进行频域分析，步骤S4对非经典结构动力响应分析过程包括：
S41、构建非经典结构动力响应频域表达式其中s＝i2πf；
S42、对xi,j(f)进行傅里叶逆变换，获得非经典结构动力响应时域表达式；
S43、多荷载输入、多位置输出的动力响应计算：固定输出位置i、循环位置j，并将每一次循环结果进行线性叠加，即可得到单输入、多输出时域动力响应时程；
S44、在单输入、多输出基础上，循环输出位置i，并根据行数为对应自由度的原则，获得结构所有位置处的时域动力响应时程。扩展了传统频域方法的理论适用范围，在工程上又保留了计算效率高的技术优势。
为了进一步说明本发明方法的实际应用价值，以分析框架结构动力响应为例进一步说明其应用。如图1所示，该框架结构长3米、宽1米、每层高1米，节点1至节点4的所有自由度固定；材料属性为：弹性模量E＝2.1×1011N/m2，密度ρ＝7850kg/m3，泊松比ν＝0.3。前5阶频率分别为6.9105、9.3617、12.1148、23.0380、29.3640Hz。
以如图2作为输入荷载，作用于节点17；以节点24动力响应作为对比结果，可以得出：
对于经典结构，阻尼矩阵C采用比例阻尼模型，分别应用本发明与传统频域方法对该结构动力响应进行分析。传递函数对比见图3。可以得知，在满足对称性条件时，本发明与传统频域方法结果一致，如图4，说明本发明方法同样适用于经典结构。
对非经典结构，构造非对称刚度矩阵及非对称阻尼矩阵其中其中分别为K、C的下三角、上三角矩阵。可以看出，与不在保持矩阵对称性；分别应用本发明与传统频域方法对该结构动力响应进行分析，传递函数对比见图5，动力响应对比见图6。可以看出，在不满足对称性条件时，本发明依然可以获得与时域方法一致的稳态解，同时克服了传统时域方法分析效率低的缺点，拓展了传统频域方法的理论适用范围。
综上，本发明可以处理质量、刚度及阻尼矩阵为对称时的结构动力响应分析问题，在满足对称性条件时本发明与传统频域方法结果一致，及对于传统经典结构，本发明同样适用，扩宽的本发明方法的应用范围；当质量、刚度及阻尼矩阵三者之中至少有一个为非对称矩阵时，本发明依然可以获得与时域方法一致的稳态解，分析效率更优，拓展了传统频域方法的理论适用范围，工程上也为大型复杂结构动力响应分析提供一种更为通用的频域方法，确保工程分析效率。
以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例应用于其它领域，但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。

Claims (5)

1.非经典结构动力响应频域方法，其特征在于，包括以下步骤：
S1、构建拉普拉斯域非经典结构传递函数：
S11、构建结构在第i个自由度上的响应：Fj(s)记为第j个荷载分量，N为荷载数目，D(s)为非经典结构阻抗的行列式，Di,j(s)为将D(s)第j列替换为Fj(s)的代数余子式；
S12、构建非经典结构传递函数：其中κ为系数，且其中η与τ分别为极值与零值；
S2、非经典结构传递函数系数获取；
S3、非经典结构传递函数极值、零值获取；
S4、非经典结构动力响应频域分析。
2.根据权利要求1所述非经典结构动力响应频域方法，其特征在于，所述步骤S2中，传递函数系数根据获取。
3.根据权利要求1所述非经典结构动力响应频域方法，其特征在于，所述步骤S2中，传递函数系数根据获取，其中M为结构质量矩阵，Mi,j为M去掉第j行、第i列。
4.根据权利要求1所述的非经典结构动力响应频域方法，其特征在于，所述步骤S3中极值、零值获取时，构建矩阵其中，M为结构质量矩阵，K为刚度矩阵，C为阻尼矩阵，I为单位阵，对进行特征值求解以获取零值和极值。
5.根据权利要求1所述非经典结构动力响应频域方法，其特征在于，所述步骤S4中非经典结构动力响应频域分析包括：
S41、构建非经典结构动力响应频域表达式其中s＝i2πf；
S42、对xi,j(f)进行傅里叶逆变换，获得非经典结构动力响应时域表达式；
S43、多荷载输入、多位置输出的动力响应计算：固定输出位置i、循环位置j，并将每一次循环结果进行线性叠加，获得单输入、多输出时域动力响应时程；
S44、在单输入、多输出基础上，循环输出位置i，并根据行数为对应自由度的原则，得到结构所有位置处的时域动力响应时程。
CN201610008351.6A 2016-01-07 2016-01-07 非经典结构动力响应频域方法 Active CN105701279B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610008351.6A CN105701279B (zh) 2016-01-07 2016-01-07 非经典结构动力响应频域方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610008351.6A CN105701279B (zh) 2016-01-07 2016-01-07 非经典结构动力响应频域方法

Publications (2)

Publication Number Publication Date
CN105701279A CN105701279A (zh) 2016-06-22
CN105701279B true CN105701279B (zh) 2018-08-24

Family

ID=56226160

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610008351.6A Active CN105701279B (zh) 2016-01-07 2016-01-07 非经典结构动力响应频域方法

Country Status (1)

Country Link
CN (1) CN105701279B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108205516B (zh) * 2017-12-26 2020-10-16 中国海洋大学 一种基于条件传递的海洋结构流固耦合动力响应方法
CN109033025B (zh) * 2018-07-18 2020-06-23 中国海洋大学 基于状态空间模型的浮式结构时域响应分析方法
CN110704798B (zh) * 2019-09-05 2021-03-02 东北大学 多自由度系统的数据处理方法、装置和存储介质
CN111062073B (zh) * 2019-12-10 2021-06-08 中国海洋大学 海工结构Laplace域动力响应瞬态分离方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101261177A (zh) * 2008-04-24 2008-09-10 上海交通大学 主动式海洋平台混合模型试验的实现方法
CN102830623A (zh) * 2012-09-07 2012-12-19 上海绿地建设(集团)有限公司 基于h∞控制的结构dtmd优化控制设计方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101261177A (zh) * 2008-04-24 2008-09-10 上海交通大学 主动式海洋平台混合模型试验的实现方法
CN102830623A (zh) * 2012-09-07 2012-12-19 上海绿地建设(集团)有限公司 基于h∞控制的结构dtmd优化控制设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Initial-condition consideration by transferring and loading reconstruction for the dynamic analysis of linear structures in the frequency domain;Fushun Liu等;《Journal of Sound and Vibration》;20141031;第164-178页 *

Also Published As

Publication number Publication date
CN105701279A (zh) 2016-06-22

Similar Documents

Publication Publication Date Title
CN105701279B (zh) 非经典结构动力响应频域方法
Kim et al. Wavelet-hybrid feedback linear mean squared algorithm for robust control of cable-stayed bridges
Ngo-Cong et al. Free vibration analysis of laminated composite plates based on FSDT using one-dimensional IRBFN method
Sha et al. Numerical investigations of the dynamic response of a floating bridge under environmental loadings
Wang et al. Comparison of ambient vibration response of the Runyang suspension bridge under skew winds with time-domain numerical predictions
CN108121865B (zh) 一种基于伴随变量法的含有多不确定参数的结构区间响应传播分析方法
CN105808884A (zh) 一种基于分形理论的有界不确定性平面裂纹应力强度因子上下界的预测方法
CN106326501A (zh) 建筑物结构动力分析用自振频率和振型的计算方法
Lian et al. Analysis of the cause and mechanism of hydraulic gate vibration during flood discharging from the perspective of structural dynamics
Tao et al. Parametric study on buffeting performance of a long-span triple-tower suspension bridge
Gur et al. Vulnerability assessment of container cranes under stochastic wind loading
CN111259469B (zh) 自振频率分析方法、装置、电子设备和存储介质
CN111783199A (zh) 一种多段式索缆结构动力特性的精细化快速求解方法
CN105787251A (zh) 海洋结构非零初始值动力响应频域改进方法
Jia Investigations of a practical wind-induced fatigue calculation based on nonlinear time domain dynamic analysis and a full wind-directional scatter diagram
CN108763667B (zh) 大曲率曲线钢-混凝土组合箱梁桥简化设计方法
Gao et al. Dynamic responses for WindFloat floating offshore wind turbine at intermediate water depth based on local conditions in China
Awan et al. Evaluation of System Identification Methods for Free Vibration Flutter Derivatives of Long-Span Bridges
Song et al. Infinite Element Static‐Dynamic Unified Artificial Boundary
Wang et al. A simulation study on the optimal control of buffeting displacement for the Sutong Bridge with multiple tuned mass dampers
Wan et al. Parametric study and dynamic response analysis of three single curved discrete pontoon floating bridges
CN107181260A (zh) 基于Matlab稀疏矩阵直角坐标牛顿法潮流计算方法
Jiang et al. Reverberation-ray analysis of moving or distributive loads on a non-uniform elastic bar
Fang et al. Fatigue damage mitigation for welded beam-to-column connections in steel high-rise buildings using passive vibration control
Sarkar et al. Stress Analysis of Concrete Gravity Dams Using Time Domain Spectral Finite Element Method Under the Action of Seismic Ground Motion

Legal Events

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