CN111061993B - 基于位移激励结构地震反应分析方法及位移激励判别方法 - Google Patents
基于位移激励结构地震反应分析方法及位移激励判别方法 Download PDFInfo
- Publication number
- CN111061993B CN111061993B CN201911277536.7A CN201911277536A CN111061993B CN 111061993 B CN111061993 B CN 111061993B CN 201911277536 A CN201911277536 A CN 201911277536A CN 111061993 B CN111061993 B CN 111061993B
- Authority
- CN
- China
- Prior art keywords
- displacement
- excitation
- equation
- ground
- formula
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于位移激励结构地震反应分析方法及位移激励判别方法,通过将基于位移激励建立的动力平衡方程分解为速度激励和位移激励两部分,再采用振型叠加法分别求解,利用状态方程计算振型叠加法的广义坐标,将动力平衡方程中由速度激励计算的广义坐标,转换为利用位移激励计算的广义坐标。并对计算结构地震反应时使用的位移激励是否适用给出了判断方法。本发明能够在保持动力平衡方程完整的基础上计算结构反应,提高位移激励方法的计算精度,同时,给出了位移激励是否适用的判断方法。本发明应用于结构的抗震分析与设计领域。
Description
技术领域
本发明涉及结构抗震分析技术领域,具体是一种基于位移激励结构地震反应分析方法及位移激励判别方法。
背景技术
结构地震反应是结构抗震设计的基础,对于绝大多数结构来说,结构的底面尺寸比较小,通常认为结构的底面处承受的地震动是相同的,结构的地震反应是基于加速度激励计算得到的。对于大跨结构,其底面尺寸比较大,结构的底面处承受的地震动是不同的,计算结构地震反应时,必须以位移激励作为荷载,求解基于绝对坐标建立动力平衡方程。但在求解基于绝对坐标建立的动力平衡方程时,为了计算方便,常常将方程中与地面速度相关的阻尼项忽略,在结构阻尼比较大时,这种忽略会导致计算结果精度不高,以此结构反应为基础的结构抗震设计不能保证结构安全可靠。求解方程时,如何实现保留此阻尼项,计算结构的地震反应引起了广大研究人员的关注。
近年来,状态空间法在结构动力反应分析中应用广泛,它将位移和速度作为独立变量来分析结构的响应,通过增加方程的维数降低微分方程的阶数。针对算法的计算精度和稳定性,国内外研究人员进行了大量、系统的研究,这些研究为完善状态空间法的理论及其应用奠定了基础。但目前,此方法在求解基于位移激励的动力平衡方程中应用并不多,一旦成功的构建保留与底面速度相关阻尼项的结构地震反应分析方法,实现保持动力平衡方程完整的基础上分析结构地震反应,算法精度会更高,结构抗震设计会更可靠。
地震动可以用加速度、速度或位移三种形式进行描述,对同一个结构进行地震反应分析,若地震动形式不同,如采用加速度激励和位移激励计算结构反应,会出现计算结果相差很大的现象。由于现有的地震监测设备记录的地震动是以加速度的形式为主,少量的设备可以记录地面速度,没有设备能记录地面位移,基于位移激励计算结构地震反应时,使用的位移激励是由地面加速度二次积分计算得到的。对加速度记录进行积分前,常常需要进行滤波去噪声、基线调整等处理,利用这些处理手段获得的位移激励是否适用于结构分析,目前还没有有效的方法判断,进而导致研究人员对位移激励方法的精度存在较大的质疑。
为解决上述问题,采用状态空间法建立基于位移激励的分析方法,并基于数学公式推导可用于判断位移激励适用性公式,这些算法能够为位移激励算法的软件应用提供基础,具有重要的理论与实用价值和广阔的应用前景。
发明内容
本发明的目的在于提供一种基于位移激励计算结构地震反应的分析方法,与此同时,给出位移激励是否适用的判别方法。
本发明所采用的技术方案如下:
基于位移激励计算结构地震反应的分析方法,步骤如下:
采用集中质量矩阵建立多自由度体系动力平衡方程:
其中,{ug}、均为m维列向量,{ug}为绝对坐标系下结构底面处地面位移,/>为绝对坐标系下结构底面处地面速度,/>为绝对坐标系下结构底面处地面加速度;{u}、/>均为n维列向量,{u}为绝对坐标系下上部结构位移,/>为绝对坐标系下上部结构速度,/>为绝对坐标系下上部结构加速度;m维列向量{Pb}表示作用于结构支撑处的力;M为上部结构的质量,C为上部结构的阻尼,K矩阵表示结构的刚度;Mgg为结构支撑处的质量,Cgg为结构支撑处的阻尼,Kgg为结构支撑处的刚度矩阵;Cg是与地面速度相关的阻尼,/>为地面速度相关的阻尼矩阵的转置,Kg是和地面位移相关的刚度,/>为地面位移相关的刚度矩阵的转置;
当已知地面运动,将方程(1)改写为:
求解式(2)时,将结构位移u表示为:
u=ud+uv (3)
其中,ud为位移激励引起的结构位移,uv为速度激励引起的结构位移;
式(2)改写为位移激励方程和速度激励方程如下:
采用振型叠加法求解公式(4)得:
其中,Φi为第i阶振型,zi为与第i阶振型对应的广义坐标;
将式(6)代入式(5)中得:
对式(10)的解表述为:
进行数值计算时,将时间域按时间步长Δt进行离散,每一时间步长Δt=tk+1-tk,由此获得式(11)的递推关系式:
其中,φ0(Δt)=eAΔt,τ为在时间步Δt内,从0变到Δt的变量;
当地面位移激励在每一激励步内,随时间线性变化时,则有:
r(tk+1-τ)=rk+1+τr1,r1=(rk-rk+1)/Δt (13)
此时,Duhamel积分为:
式(12)进一步写为:
采用振型叠加法求解式(5)获得:
广义坐标由式(17)获得:
式(18)的解表述为:
式(19)的递推式为:
对Duhamel积分进行分部积分,将速度激励转换为位移激励,式(20)可写为:
将式(15)的结果和式(21)的结果叠加,得到保留与地面速度相关阻尼项动力平衡方程振型的广义坐标,再利用式(6)、式(16)和式(3)获得结构地震反应。
位移激励是否适用的判断方法,具体步骤如下:
基于位移激励的单质点体系绝对加速度反应谱Sad为:
基于加速度激励的单质点体系绝对加速度反应谱Sa为:
Sad和Sa满足下式时,加速度激励和位移激励对应的地面运动相同;
本发明的有益效果:
本发明所建立的基于位移激励结构地震反应分析方法,提高了位移激励时分析方法的计算精度,以此分析结果为基础进行结构抗震设计,结构会更安全可靠;提出的位移激励适用性判断方法,为由加速度积分计算位移激励的积分方法选用、滤波参数选取提供有理论依据。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为算例结构参数图与地震动荷载施加图。
图2位移激励和加速度激励时程曲线图。
图3位移激励与加速度激励单质点体系绝对加速度反应谱。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
基于位移激励计算结构地震反应的分析方法:计算如下:
采用集中质量矩阵建立多自由度体系动力平衡方程:
其中,{ug}、均为m维列向量,{ug}为绝对坐标系下结构底面处地面位移,/>为绝对坐标系下结构底面处地面速度,/>为绝对坐标系下结构底面处地面加速度;{u}、/>均为n维列向量,{u}为绝对坐标系下上部结构位移,/>为绝对坐标系下上部结构速度,/>为绝对坐标系下上部结构加速度;m维列向量{Pb}表示作用于结构支撑处的力;M为上部结构的质量,C为上部结构的阻尼,K矩阵表示结构的刚度;Mgg为结构支撑处的质量,Cgg为结构支撑处的阻尼,Kgg为结构支撑处的刚度矩阵;Cg是与地面速度相关的阻尼,为地面速度相关的阻尼矩阵的转置,Kg是和地面位移相关的刚度,/>为地面位移相关的刚度矩阵的转置;
当已知地面运动,将方程(1)改写为:
求解式(2)时,将结构位移u表示为:
u=ud+uv (3)
其中,ud为位移激励引起的结构位移,uv为速度激励引起的结构位移;
式(2)改写为位移激励方程和速度激励方程如下:
采用振型叠加法求解公式(4)得:
其中,Φi为第i阶振型,zi为与第i阶振型对应的广义坐标;
将式(6)代入式(5)中得:
对式(10)的解表述为:
进行数值计算时,将时间域按时间步长Δt进行离散,每一时间步长Δt=tk+1-tk,由此获得式(11)的递推关系式:
其中,φ0(Δt)=eAΔt,τ为在时间步Δt内,从0变到Δt的变量;
当地面位移激励在每一激励步内,随时间线性变化时,则有:
r(tk+1-τ)=rk+1+τr1,r1=(rk-rk+1)/Δt (13)
此时,Duhamel积分为:
式(12)进一步写为:
采用振型叠加法求解式(5)获得:
广义坐标由式(17)获得:
式(18)的解表述为:
式(19)的递推式为:
对Duhamel积分进行分部积分,将速度激励转换为位移激励,式(20)可写为,
将式(15)的结果和式(21)的结果叠加,得到保留与地面速度相关阻尼项动力平衡方程振型的广义坐标,再利用式(6)、式(16)和式(3)获得结构地震反应。
位移激励是否适用的判断方法,具体如下:
基于位移激励的单质点体系绝对加速度反应谱Sad为:
基于加速度激励的单质点体系绝对加速度反应谱Sa为:
当Sad和Sa满足下式(24)时,加速度激励和位移激励对应的地面运动相同,此时位移激励适用于计算结构地震反应;
具体计算实施例:
以图1所示的结构为例说明,为了比较本发明算法结果与加速度激励方法结果的差别,采用的地震动激励如图2所示。采用matlab语言实现上述理论算法的程序,位移激励由加速度激励二次积分获得,利用式(22)、(23)计算图2所示地震动激励,获得位移激励和加速度激励的单质点体系绝对加速度反应谱,如图3所示,两者的差别满足式(24)。表1为结构地震反应差值百分比,从中可以看出,本发明方法的计算结果与常用的加速度激励方法结果相差很小,最大的差别在6.5%,这说明位移激励方法和加速度激励方法的结果相差不大,本发明的位移激励适用于计算结构地震反应。
表1结构地震反应差值百分比
以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。
Claims (2)
1.基于位移激励计算结构地震反应的分析方法,其特征在于,步骤如下:
采用集中质量矩阵建立多自由度体系动力平衡方程:
其中,{ug}、均为m维列向量,{ug}为绝对坐标系下结构底面处地面位移,/>为绝对坐标系下结构底面处地面速度,/>为绝对坐标系下结构底面处地面加速度;{u}、 均为n维列向量,{u}为绝对坐标系下上部结构位移,/>为绝对坐标系下上部结构速度,/>为绝对坐标系下上部结构加速度;m维列向量{Pb}表示作用于结构支撑处的力;M为上部结构的质量,C为上部结构的阻尼,K矩阵表示结构的刚度;Mgg为结构支撑处的质量,Cgg为结构支撑处的阻尼,Kgg为结构支撑处的刚度矩阵;Cg是与地面速度相关的阻尼,/>为地面速度相关的阻尼矩阵的转置,Kg是和地面位移相关的刚度,/>为地面位移相关的刚度矩阵的转置;
当已知地面运动,将方程(1)改写为:
求解式(2)时,将结构位移u表示为:
u=ud+uv (3)
其中,ud为位移激励引起的结构位移,uv为速度激励引起的结构位移;
式(2)改写为位移激励方程和速度激励方程如下:
采用振型叠加法求解公式(4)得:
其中,Φi为第i阶振型,zi为与第i阶振型对应的广义坐标;
将式(6)代入式(5)中得:
对式(10)的解表述为:
进行数值计算时,将时间域按时间步长Δt进行离散,每一时间步长Δt=tk+1-tk,由此获得式(11)的递推关系式:
其中,φ0(Δt)=eAΔt,τ为在时间步Δt内,从0变到Δt的变量;
当地面位移激励在每一激励步内,随时间线性变化时,则有:
r(tk+1-τ)=rk+1+τr1,r1=(rk-rk+1)/Δt (13)
此时,Duhamel积分为:
式(12)进一步写为:
采用振型叠加法求解式(5)获得:
广义坐标由式(17)获得:
式(18)的解表述为:
式(19)的递推式为:
对Duhamel积分进行分部积分,将速度激励转换为位移激励,式(20)可写为:
将式(15)的结果和式(21)的结果叠加,得到保留与地面速度相关阻尼项动力平衡方程振型的广义坐标,再利用式(6)、式(16)和式(3)获得结构地震反应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911277536.7A CN111061993B (zh) | 2019-12-11 | 2019-12-11 | 基于位移激励结构地震反应分析方法及位移激励判别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911277536.7A CN111061993B (zh) | 2019-12-11 | 2019-12-11 | 基于位移激励结构地震反应分析方法及位移激励判别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111061993A CN111061993A (zh) | 2020-04-24 |
CN111061993B true CN111061993B (zh) | 2023-03-28 |
Family
ID=70300938
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911277536.7A Active CN111061993B (zh) | 2019-12-11 | 2019-12-11 | 基于位移激励结构地震反应分析方法及位移激励判别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111061993B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113484911B (zh) * | 2021-04-30 | 2022-11-15 | 石家庄学院 | 一种基于地震动演化功率谱的地震位移反应计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106096119A (zh) * | 2016-06-07 | 2016-11-09 | 西北工业大学 | 基于大质量法的随机加速度激励下结构拓扑优化设计方法 |
CN109145386A (zh) * | 2018-07-24 | 2019-01-04 | 张瑞杰 | 一种用于工程结构的地震碰撞反应分析方法 |
CN110083920A (zh) * | 2019-04-23 | 2019-08-02 | 西南交通大学 | 一种地震作用下非比例阻尼结构随机响应的分析方法 |
CN110348158A (zh) * | 2019-07-18 | 2019-10-18 | 中国水利水电科学研究院 | 一种基于分区异步长解的地震波动分析方法 |
-
2019
- 2019-12-11 CN CN201911277536.7A patent/CN111061993B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106096119A (zh) * | 2016-06-07 | 2016-11-09 | 西北工业大学 | 基于大质量法的随机加速度激励下结构拓扑优化设计方法 |
CN109145386A (zh) * | 2018-07-24 | 2019-01-04 | 张瑞杰 | 一种用于工程结构的地震碰撞反应分析方法 |
CN110083920A (zh) * | 2019-04-23 | 2019-08-02 | 西南交通大学 | 一种地震作用下非比例阻尼结构随机响应的分析方法 |
CN110348158A (zh) * | 2019-07-18 | 2019-10-18 | 中国水利水电科学研究院 | 一种基于分区异步长解的地震波动分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111061993A (zh) | 2020-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Besselink et al. | A semi-analytical study of stick-slip oscillations in drilling systems | |
CN109446715B (zh) | 一种大跨度桥梁地震响应时程分析方法 | |
Zhu et al. | Multi-type sensor placement for multi-scale response reconstruction | |
CN111061993B (zh) | 基于位移激励结构地震反应分析方法及位移激励判别方法 | |
Lankarani | A Poisson-based formulation for frictional impact analysis of multibody mechanical systems with open or closed kinematic chains | |
CN106326501A (zh) | 建筑物结构动力分析用自振频率和振型的计算方法 | |
Rahman | Wave diffraction by large offshore structures: an exact second-order theory | |
CN110795790B (zh) | 一种复杂建筑结构非线性动力时程分析方法 | |
Ye et al. | Enhancing the seismic performance of dual-isolated structure with optimal supplemental inerters | |
Tian et al. | Vibration Characteristics of Heavy‐Duty CNC Machine Tool‐Foundation Systems | |
CN111046456A (zh) | 地震动时程激励下线性系统反应谱分析的迭代法 | |
Chen | Parametric instability of spinning twisted Timoshenko beams under compressive axial pulsating loads | |
Wang et al. | Optimal design and experimental study of a multidynamic vibration absorber for multifrequency excitation | |
Zu et al. | An improved transfer matrix method for steady-state analysis of nonlinear rotor-bearing systems | |
Jena et al. | Response of cracked cantilever beam subjected to traversing mass | |
Villa et al. | The invariant manifold approach applied to nonlinear dynamics of a rotor-bearing system | |
Zeljkovic et al. | Experimental and computer aided analysis of high-speed spindle assembly behaviour | |
Wang | Vibration and sensitivity analysis of a beam with a lumped mass of translational and rotary inertias | |
CN110990910A (zh) | 时程激励下线性耗能结构响应的快速迭代法 | |
CN114091207B (zh) | 一种冲击载荷下的隔冲器时变可靠性分析方法 | |
Farooq et al. | An experimental investigation of state-variable modal decomposition for modal analysis | |
CN113919590A (zh) | 不同产状非等间距斜条分岩质边坡三维稳定性预测方法 | |
Chung | Analysis of vibration in L-and H-shaped planar frame structures using transfer function | |
JPH11160144A (ja) | 3次元地震応答解析方法 | |
CN110765538A (zh) | 一种复杂结构非线性动力分析的改进的广义α法 |
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 |