CN109460622B - 一种大规模建筑结构的完全显式的动力时程分析方法 - Google Patents
一种大规模建筑结构的完全显式的动力时程分析方法 Download PDFInfo
- Publication number
- CN109460622B CN109460622B CN201811356374.1A CN201811356374A CN109460622B CN 109460622 B CN109460622 B CN 109460622B CN 201811356374 A CN201811356374 A CN 201811356374A CN 109460622 B CN109460622 B CN 109460622B
- Authority
- CN
- China
- Prior art keywords
- matrix
- damping
- unit
- calculating
- time
- 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
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/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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/06—Power analysis or power optimisation
Abstract
本发明涉及一种大规模建筑结构的动力分析数值模拟方法,具体涉及一种无需阻尼矩阵对角化、具有一定方法阻尼、具有大稳定步长、具有3阶精度、自起步的完全显式方法。其步骤如下:第一步,对大规模建筑结构进行空间有限元离散,导出离散系统的运动方程组;第二步,根据体系最短周期,选取分析步长;第三步,逐时间步计算第i个时间步上ti时刻的位移ui、速度vi;第四步,计算高阶精度的位移、速度,并计算加速度;本发明显式动力分析方法与中心差分法等常见的显式动力分析方法相比,具备适用于非对角阻尼矩阵、适用于非线性阻尼系统、可抑制虚假的高频模式、稳定区间长、精度高收敛快、自启动等特点,且操作步骤简单,容易实施,具有很强的工程应用价值。
Description
技术领域
本发明涉及一种大规模建筑结构的动力分析数值模拟方法,具体涉及一种大规模建筑结构的无需阻尼矩阵对角化、具有一定方法阻尼、具有大稳定步长、具有3阶精度、自起步的完全显式的动力时程分析方法。
背景技术
近年来,国内大规模建筑结构日趋增多,尽管大规模建筑结构满足了城市化需求,但却存在地震风险。为规避大规模建筑结构的地震风险,在设计阶段,常常采用动力分析数值模拟手段,对大规模建筑结构进行大量的动力时程分析,以保证大规模建筑结构在突发强震作用下的安全。为便于计算,人们提出了一些例如底部剪力法、反应谱法等方法。然而随着土木工程结构形式的不断发展和对时程分析可靠度越来越高的要求,我们需要对结构在动力荷载下的行为进行更为精确的预测。加上有限元方法的发展与计算机性能的迅猛提高,人们对时程分析法的关注越来越多。动力时程分析法要求我们计算出结构位移、速度等量随时间的定量变化,是描述结构在动力作用下力学行为的最精准的分析方法。对于普遍的非线性结构,直接积分方法是进行动力时程分析的最强有力工具。
直接积分方法可以分为隐式方法与显式方法。常用的隐式方法有平均加速度法、Wilson-θ法、HHT-α法等等,这些隐式方法均为无条件稳定方法,计算中可以采用的时间步长不受模型网格尺寸限制,效率较高,且有一定的数值阻尼,能够过滤掉结构空间有限元划分造成的虚假高阶模态,因此常常被应用于结构动力时程分析这种动力响应以低阶模态为主的分析。但由于隐式方法在每一增量步内都需要对静态平衡方程进行迭代求解,并且每次迭代都需要求解大型的线性方程组,这一过程需要占用相当数量的计算资源、磁盘空间和内存,反而表现出较低的计算效率,且对于非线性比较强的复杂结构模型,非线性迭代不容易收敛。因此,近些年以中心差分法为代表的显式方法逐渐被应用于结构动力时程分析,如国际知名软件ABAQUS、LS-DYNA等等。相比隐式方法,显式方法的每步计算中无需求解大规模方程组,仅需要进行单元级别的矩阵和向量乘法运算以及整体向量加法运算,在计算效率上有明显的优势,且显式方法不存在无法收敛的问题,因此在计算强震作用下建筑结构的动力响应、倒塌模拟这类复杂模型的强非线性响应方面具有先天的优势,近年来国内也逐渐的涌现出一批采用显式方法进行结构动力分析的软件,如PKPM-SAUSAGE等。
尽管中心差分方法具有效率上的优势,但仍然存在以下几个非常棘手的问题:(1)为避免求解大规模的整体平衡方程组,中心差分法要求阻尼矩阵为对角矩阵,而结构动力时程分析中,阻尼矩阵的对角化往往通过采用质量阻尼或者振型阻尼来实现,这种处理往往会引入较大的误差,限制了该方法的应用范围。若采用对速度的偏心差分代替中心差分,可避免阻尼矩阵对角化,但偏心差分方法只有一阶精度,计算结果可靠性收到影响;(2)当结构中存在非线性粘滞阻尼器,使得阻尼矩阵为非线性矩阵时,若采用中心差分法,每个时间步需要像隐式方法一样进行非线性迭代并求解大规模方程组,显式方法的效率优势和收敛优势则无从体现;(3)结构的空间有限元离散,不可避免的引入了虚假的高阶模态,由于中心差分法无数值耗散特性,无法抑制这些虚假高频响应,会造成很大的误差;(4)显式方法普遍存在稳定性问题,分析步长不允许超过临界步长,临界步长正比于结构最短周期,中心差分法的临界步长为结构最短周期的1/Π倍,普遍被认为是显式方法中最长的,但临界步长仍有拓展的空间。以上几大问题制约了中心差分方法在结构动力时程分析中的应用。除此之外,由于中心差分法为多步法,还存在方法启动问题,造成了编程上的麻烦。另外,中心差分方法具有二阶精度,尽管已经满足大部分工程计算的需求,但是对于超长持时分析,误差的累积会造成计算结果偏离理论解。
发明内容
针对中心差分法在大规模复杂建筑结构非线性动力时程分析方面的缺陷,本发明的目的在于提供一种大规模建筑结构的无需阻尼矩阵对角化、具有一定方法阻尼、具有大稳定步长、具有3阶精度、自起步的完全显式的动力时程分析方法。
本发明的技术方案如下:
一种大规模建筑结构的显式动力分析方法,步骤如下:
第一步,对大规模建筑结构进行空间有限元离散,建立大规模建筑结构的有限元模型离散系统,梁柱均采用伯努利欧拉梁单元,并采用Rayleigh阻尼建立单元阻尼矩阵,由单元刚度矩阵、单元质量矩阵和单元阻尼矩阵集成整体刚度矩阵K、整体质量矩阵M和整体阻尼矩阵C,并由Hamilton原理导出离散系统的运动方程组,其中F为右端项:
第二步,根据体系最短周期,选取分析步长;计算建筑结构有限元模型的最短周期,记为Tmin;选取时间步长Δt<0.45Tmin,令Δt=0.4Tmin;
第三步,逐时间步计算,对于第i个时间步,已知ti-1时刻的位移ui-1和速度vi-1,由下式计算第i个时间步上ti时刻的位移ui、速度vi;
其中,Pi的表达式为
ei=M-1(P1-H11ui-1-H12ui-H13vi-1-H14vi)Δt
其中系数矩阵Hij的表达式为
优选地,所述的第一步当中,单元刚度矩阵为:
优选地,所述的第一步当中,单元质量矩阵为:
优选地,所述的第一步当中,采用Rayleigh阻尼,建立单元阻尼矩阵
Ce=a0Me+a1Ke
其中
ωi和ωj一般分别取结构的第1阶和第3阶频率,ζ为阻尼比,一般为0.05。
优选地,首先计算并存贮各个单元的单元刚度矩阵Ke、质量矩阵Me和阻尼矩阵Ce以及右端项Fe;其中Me为对角矩阵,仅需要对质量矩阵Me进行集成得到整体质量矩阵M,而无需计算整体刚度矩阵C和整体阻尼矩阵M以及右端项F。
优选地,所述的第二步当中,建筑结构有限元模型的最短周期Tmin可采用估计的方法,无需实际计算。
本发明与现有技术相比,优点在于:
(1)本发明的动力时程分析方法只需要对质量矩阵对角化,而无需对阻尼矩阵对角化,因此可以采用常用的Rayleigh阻尼等非对角阻尼矩阵;
(2)本发明的动力时程分析方法同样适用于存在非线性粘滞阻尼器的结构体系的显式动力时程分析,每个时间步上无需迭代求解;
(3)本发明的动力时程分析方法对于高阶模态具有较强的数值耗散特性,可以屏蔽掉由于空间离散造成的虚假的高频响应;
(4)本发明的动力时程分析方法具有较长的稳定区间,对于无阻尼系统,方法临界步长为结构最短周期的0.4502倍,比中心差分法的临界步长增大近50%,谱半径曲线,如图1所示;
(5)本发明的动力时程分析方法具有高阶精度,对于无阻尼体系,方法具有4阶精度,对于阻尼体系,方法具有3阶精度,均高于中心差分法;其中加速度收敛阶与中心差分法的对比如图2所示;
(6)本发明的动力时程分析方法为单步法,不存在起步问题;
(7)本发明的动力时程分析方法实施简便,效率高,只需在单元层级进行矩阵向量乘法,而无需集成整体刚度、质量、阻尼矩阵进行计算,可并行。
附图说明
图1本发明的动力时程分析方法与中心差分法的谱半径曲线对比,其中虚线为中心差分法,点划线为本发明方法。
图2本发明的动力时程分析方法与中心差分法的加速度收敛性对比,其中虚线为中心差分法,点划线为本发明方法。
图3某2层框架结构。
图4El-Centro波。
图5本发明的动力时程分析方法与中心差分法的加速度精度对比,其中点划线为中心差分法,实线为本发明方法。
具体实施方式
下面结合具体实施例来对本发明进行进一步说明,但并不将本发明局限于这些具体实施方式。本领域技术人员应该认识到,本发明涵盖了权利要求书范围内所可能包括的所有备选方案、改进方案和等效方案。
下面结合附图对本发明的结构原理和工作原理作具体的描述:
以一个2层框架结构为实例,具体阐述本发明所涉及的完全显式动力时程分析方法,框架结构如图3所示,下部箭头方向为地震波,该框架结构的有限元模型包含2层,水平两个方向跨度各为6m,层高4米。梁的截面尺寸为0.2×0.4m,柱截面尺寸为0.4×0.4m,密度为2.5×103kg/m3。所述的完全显式的动力时程分析方法包含如下步骤:
第一步,对上述框架结构进行空间有限元离散,建立上述框架结构的有限元模型离散系统,梁柱均采用伯努利欧拉梁单元,并采用Rayleigh阻尼建立单元阻尼矩阵,由单元刚度矩阵、单元质量矩阵和单元阻尼矩阵集成整体刚度矩阵K、整体质量矩阵M和整体阻尼矩阵C,并由Hamilton原理导出离散系统的运动方程组,其中F为右端项:
第二步,选取分析步长。为与中心差分法对比,选取满足中心差分法稳定性的步长,令Δt=0.002s;
第三步,逐时间步计算,对于第i个时间步,已知ti-1时刻的位移ui-1和速度vi-1,由下式计算第i个时间步上ti时刻的位移ui、速度vi;
其中,Pi的表达式为
其中系数矩阵Hij的表达式为
在基底输入如图4所示的El Centro地震波,该地震波的加速度记录间隔为0.01s,计算框架一层顶部A点位移响应。采用集中质量矩阵,阻尼比选取0,各个质点的初始位移和初始速度均为0。
为了展示本发明的完全显式动力时程分析方法高效性和精确性,首先采用本发明方法进行一次动力分析,分析步长选取第二步中所提到的Δt=0.002s。同时采用中心差分法对该框架进行动力分析,同样选取步长同样选取Δt=0.002s,为方便对比,选取更小的步长Δt=0.0002s,采用中心差分法进行计算,作为精确解的参考解。由本发明方法和中心差分法计算得到的A点加速度响应如图5所示,图中只给出39秒-40秒的响应。经对比发现,采用本发明的动力时程分析方法的计算结果与参考解几乎重合,而中心差分法的计算结果严重偏离参考解,该例充分验证了本发明显式动力时程分析方法高效性和精确性。
应当理解的是,本发明描述的方法的步骤仅仅是示例性的描述,对其先后进行的时间顺序没有特殊的要求,除非其本身有必然的先后顺序关系。
如上所示,本发明虽然已参照有限的实施例和附图进行了说明,但在本发明所属领域中具备通常知识的人均可以从此记载中进行各种修改和变形。由此,其他实施例及权利要求书与等同物均属于权利要求的保护范围。
Claims (9)
1.一种大规模建筑结构的完全显式的动力时程分析方法,其特征在于,步骤如下:
第一步,对大规模建筑结构进行空间有限元离散,建立大规模建筑结构的有限元模型离散系统,梁柱均采用伯努利欧拉梁单元,并采用Rayleigh阻尼建立单元阻尼矩阵,由单元刚度矩阵、单元质量矩阵和单元阻尼矩阵集成整体刚度矩阵K、整体质量矩阵M和整体阻尼矩阵C,并由Hamilton原理导出离散系统的运动方程组,其中F为右端项:
第二步,根据体系最短周期,选取分析步长;计算建筑结构有限元模型的最短周期,记为Tmin;选取时间步长Δt<0.45Tmin,令Δt=0.4Tmin;
第三步,逐时间步计算,对于第i个时间步,已知ti-1时刻的位移ui-1和速度vi-1,由下式计算第i个时间步上ti时刻的位移ui、速度vi;
其中,Pi的表达式为
ei=M-1(P1-H11ui-1-H12ui-H13vi-1-H14vi)Δt
其中系数矩阵Hij的表达式为
5.根据权利要求1所述的大规模建筑结构的完全显式的动力时程分析方法,其特征在于,所述的第一步当中,对于单元右端项Fe的计算,需对每个时刻的地震波加速度记录采用线性插值得到üg(t),由下式计算
Fe=-Meüg(t);其中Me为单元质量矩阵。
6.根据权利要求1所述的大规模建筑结构的完全显式的动力时程分析方法,其特征在于,首先计算并存贮各个单元的单元刚度矩阵Ke、质量矩阵Me和阻尼矩阵Ce以及右端项Fe;其中Me为对角矩阵,仅需要对质量矩阵Me进行集成得到整体质量矩阵M,而无需计算整体刚度矩阵C和整体阻尼矩阵M以及右端项F。
7.根据权利要求1所述的大规模建筑结构的完全显式的动力时程分析方法,其特征在于,所述的第二步当中,建筑结构有限元模型的最短周期Tmin可采用估计的方法,无需实际计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811356374.1A CN109460622B (zh) | 2018-11-15 | 2018-11-15 | 一种大规模建筑结构的完全显式的动力时程分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811356374.1A CN109460622B (zh) | 2018-11-15 | 2018-11-15 | 一种大规模建筑结构的完全显式的动力时程分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109460622A CN109460622A (zh) | 2019-03-12 |
CN109460622B true CN109460622B (zh) | 2023-03-07 |
Family
ID=65610481
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811356374.1A Active CN109460622B (zh) | 2018-11-15 | 2018-11-15 | 一种大规模建筑结构的完全显式的动力时程分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109460622B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110348158B (zh) * | 2019-07-18 | 2021-07-27 | 中国水利水电科学研究院 | 一种基于分区异步长解的地震波动分析方法 |
CN110795790B (zh) * | 2019-10-31 | 2021-02-12 | 李鲁 | 一种复杂建筑结构非线性动力时程分析方法 |
CN110765538B (zh) * | 2019-10-31 | 2021-02-12 | 李鲁 | 一种复杂结构非线性动力分析的改进的广义α法 |
JP7367907B2 (ja) | 2019-11-12 | 2023-10-24 | 株式会社アークデータ研究所 | 構造計算支援装置 |
CN114722664B (zh) * | 2022-03-30 | 2023-06-20 | 哈尔滨理工大学 | 一种多材料结构动刚度拓扑优化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106326501A (zh) * | 2015-06-15 | 2017-01-11 | 上海东浩兰生国际服务贸易(集团)有限公司 | 建筑物结构动力分析用自振频率和振型的计算方法 |
CN107729603A (zh) * | 2017-09-04 | 2018-02-23 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 复合材料壁板螺栓连接的非线性分析方法 |
CN108182338A (zh) * | 2018-03-19 | 2018-06-19 | 重庆大学 | 非岩质坡地水平地震动放大系数确定方法及抗震设计方法 |
CN108416182A (zh) * | 2018-05-31 | 2018-08-17 | 长安大学 | 一种基于定量分析的明沟隔振的设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10061878B2 (en) * | 2015-12-22 | 2018-08-28 | Dassault Systemes Simulia Corp. | Effectively solving structural dynamics problems with modal damping in physical coordinates |
-
2018
- 2018-11-15 CN CN201811356374.1A patent/CN109460622B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106326501A (zh) * | 2015-06-15 | 2017-01-11 | 上海东浩兰生国际服务贸易(集团)有限公司 | 建筑物结构动力分析用自振频率和振型的计算方法 |
CN107729603A (zh) * | 2017-09-04 | 2018-02-23 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 复合材料壁板螺栓连接的非线性分析方法 |
CN108182338A (zh) * | 2018-03-19 | 2018-06-19 | 重庆大学 | 非岩质坡地水平地震动放大系数确定方法及抗震设计方法 |
CN108416182A (zh) * | 2018-05-31 | 2018-08-17 | 长安大学 | 一种基于定量分析的明沟隔振的设计方法 |
Non-Patent Citations (5)
Title |
---|
"a novel time continous p-galerkin scheme with 2p-oeder accuracy for seismic dynamic analysis";Xu Junjie等;《16th World Conference on Earthquake Engineering》;20171231;全文 * |
大规模建筑结构非线性显式动力分析的振型阻尼实现;曹胜涛等;《建筑结构》;20160110(第01期);全文 * |
显式分析方法在高层建筑弹塑性地震反应分析中的适用性研究;杜轲等;《地震工程与工程振动》;20150815(第04期);全文 * |
结构动力反应分析方法研究;廖旭;《建筑结构》;20161025(第20期);全文 * |
阻尼矩阵选取对核电厂结构地震响应的影响分析;李小军等;《振动与冲击》;20150115(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109460622A (zh) | 2019-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109460622B (zh) | 一种大规模建筑结构的完全显式的动力时程分析方法 | |
Li et al. | A novel family of controllably dissipative composite integration algorithms for structural dynamic analysis | |
Bobby et al. | Data-driven performance-based topology optimization of uncertain wind-excited tall buildings | |
Soroushian | Integration step size and its adequate selection in analysis of structural systems against earthquakes | |
Flodén et al. | Reduction methods for the dynamic analysis of substructure models of lightweight building structures | |
Zhu et al. | Topology optimization for linear stationary stochastic dynamics: applications to frame structures | |
Pagani et al. | Higher-order theories and radial basis functions applied to free vibration analysis of thin-walled beams | |
Hooijkamp et al. | Topology optimization for linear thermo‐mechanical transient problems: modal reduction and adjoint sensitivities | |
CN109409006B (zh) | 一种超高层结构动力时程分析方法 | |
Rezaiee-Pajand et al. | More accurate and stable time integration scheme | |
Das et al. | Reliability-based optimization of nonlinear energy sink with negative stiffness and sliding friction | |
CN104405054A (zh) | 一种设置黏滞消能器的结构设计方法 | |
WO2021082959A1 (zh) | 一种复杂建筑结构非线性动力时程分析方法 | |
Wang et al. | Efficiency analysis of numerical integrations for finite element substructure in real-time hybrid simulation | |
Maghareh et al. | Adaptive multi‐rate interface: development and experimental verification for real‐time hybrid simulation | |
Ehlers et al. | Stability analysis of finite difference schemes revisited: A study of decoupled solution strategies for coupled multifield problems | |
CN106295028A (zh) | 一种局部结构动力学建模方法及装置 | |
Sharma et al. | Free vibration analysis of moderately thick antisymmetric cross-ply laminated rectangular plates with elastic edge constraints | |
Karttunen et al. | Shear deformable plate elements based on exact elasticity solution | |
Hou et al. | Static and dynamic analysis of two-layer Timoshenko composite beams by weak-form quadrature element method | |
Li et al. | Efficient Woodbury-CA hybrid method for structures with material and geometric nonlinearities | |
Naets et al. | Multi‐expansion modal reduction: A pragmatic semi–a priori model order reduction approach for nonlinear structural dynamics | |
He et al. | New speedup algorithms for nonlinear dynamic time history analysis of supertall building structures under strong earthquakes | |
CN108376192B (zh) | 一种确定模态叠加法计算加速度反应所需振型数目的方法 | |
Lavan et al. | Progressive collapse analysis through strength degradation and fracture in the mixed Lagrangian formulation |
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 |