[动力学][地震工程] 一个有趣的问题: SPECTR中的Newmark-Beta法计算反应谱发散？

实干、实践、积累、思考、创新。

---

一个遗留了好几年的问题，起因是一位网友找了条十分特殊的自己合成波测试 SPECTR 2015 (当时搞选波工具，SPECTR是仿Seismosignal的界面做的)，结果发现Newmark法计算的反应谱结果是发散的，之所以说“特殊”，是因为其他波都没出现这个问题（其他波不同方法计算都是重合的），由于一般波测试没问题，而且默认用的逐步精确法，十分稳健，当时没太在意，就去研究其他东西去了，但是问题倒是记录到了word，一直放在硬盘的某个角落，电脑换了几轮，但记录的这个word每次都舍不得删掉，虽然极其碍眼，今天又整理硬盘，又看到了这个word，这次决定把这个word删掉，一次性解决问题，于是打开以前的代码，果然看回好几年前写的代码，写得的确丑，还有点乱，但仔细检查公式倒是没问题，那问题到底出在哪？最后灵机一动，发现问题在加速度的递推公式上，当时采用了增量法的NEWMARK递推公式，各step的加速度均通过上一步的加速度加上加速度增量逐步累加得到，而可能对于这条波，恰好出现了浮点数精度丢失。于是把加速度改为平衡方程反推，果然问题解决了！！！改天更新SPECTR 2022!!!

PS：浮点数精度丢失，整型数溢出是编程中很难发现的bug，隐藏最深。Floating point precision lost and integer overflow is a hard-to-find bug in programming.

• 相关话题 ( Related Topics)

[01]. [Tool] SPECTR – A program for Response Spectra Analysis [反应谱计算程序]

[02]. [程序][Tool] Ground Motion Selection [强震记录选取]

[03]. [程序][软件]Ground Motion Library [强震记录管理]

[04]. Artificial ground motion generator [人工模拟地震动]

[05]. ATC-63 Ground Motion Record Sets [ATC-63推荐的地震动记录集]

[06]. [Tool][软件]PEER NGA Ground Motion Convertor [PEER地震波数据格式转换器]

[07]. Frequency Analysis of Ground Motion by MATLAB [MATLAB地震动频域分析]

[08]. [Tool] [程序] GB18306-2015《地震动参数区划图》地震动参数计算

[09]. [Tool]水电工程水工建筑物标准设计反应谱计算程序(NB 35047-2015)

[10]. GB18306-2015地震动参数确定步骤与实例 (Ground motion parameter determination steps and examples of GB18306-2015)

[11]. 设计基本地震加速度与水平地震影响系数最大值

[12]. Tripartite Plot of Response Spectra [三联反应谱的基本概念与绘制]

[13]. Earthquake Response Spectra Analysis by MATLAB [MATLAB地震反应谱分析]

[14]. YJK地震波反应谱分析与地震波选取注意事项

[15]. 地震反应谱的特性(The Characteristics of Earthquake Response Spectrum)

[16]. Basics of Earthquake Response Spectra[地震反应谱的基本概念]

[17]. [Tool]中国规范反应谱生成程序 [Chinese Code’s Design Response Spectrum]

[18]. ETABS 2015人工波功能初步测试(一)

[19]. ETABS 2015人工波功能初步测试(二)

[20]. [Tool][软件][动力学]NSDOF v2016: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF: 单自由度系统动力非线性分析工具)

[21]. Dynamic Analysis of SDOF System by MATLAB [MATLAB单自由度系统动力响应分析]

[22]. 地震波基线是否修正对结构分析结果的影响（一）

[23]. [科研][工具][Research] 反应谱拟合及反应谱特征参数提取工具 [A tool for fitting response spectrum and extracting response spectrum parameters]

[24]. [Software][科研小工具] 力-位移累积耗能计算[Force-Displacement Cumulative Dissipated Energy Calculator ]

[25]. [Tool][更新][Updates] SPECTR (v1.1) – A program for Response Spectra Analysis [SPECTR反应谱计算程序更新]

[26]. GMS Converter: 通用格式转换器 [General Formats Transformer for Earthquake Records]

[27]. [编程][软件] PERFORM-3D 地震波时程转换工具

[28]. [软件][工具][地震][Tool] 结构抗震分析地震动强度指标计算工具 ( A tool for Calculating Earthquake Intensities for Seismic Analysis of Structures)

[29]. [软件][地震动][Tool] 地震动参数计算软件GMP：补充介绍及说明

[30]. [地震工程][科研]地震动能量反应谱 (Earthquake Energy Demand Spectra)

[31]. [地震动][软件][研究][Tool] 反应谱阻尼效应系数Dh计算工具 (A Tool for Calculating Damping Correction Factor of Earthquake Spectrum）

[32]. [减隔震][地震工程][软件][Tool] 软钢·摩擦阻尼器减震结构的性能曲线计算

[33]. [地震工程][科研][软件][Tool] 等延性反应谱计算工具 (Construction of Constant-Ductility Response Spectrum)

[34]. [选波][地震波][Tool] GMS: Ground Motion Selection Program Updates (GMS选波软件更新)

[35]. [选波][地震波] 记录一个实际工程有趣的选波问题

[36]. [PEER][地震波] 关于PEER强震数据库地震波的单位及数据格

[37]. [科研][更新][Update] 地震动参数计算软件(GMP)更新 (支持更多数据导入格式)

[38]. YJK1.7人工波功能测试

[39]. [工程][项目][地震动][软件] EPA Scaling: EPA-based Adjustment Tool of Seismic Ground Motion Record[基于EPA的地震动记录调整工具]

[40]. [数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具]

[41]. [数学][软件] FOUR_TRAN Example 3: Earthquake Ground Acceleration Frequency Spectrum Analysis [FOUR_TRAN傅里叶分析工具使用案例3: 地震波频谱分析]

[42]. [工程][项目][地震动][软件] EPA Scaling: EPA-based Adjustment Tool of Seismic Ground Motion Record[基于EPA的地震动记录调整工具]

[43]. [Tool][软件][Dynamics] NSDOF v2020: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2020)

[44]. [软件][工具][地震波] KNETConvert: Kik-net Ground Motion Converter [K-Net, Kik-net 地震波数据转换工具]

[45]. [地震工程][科研][软件] IRSA 2020: Inelastic Response Spectra Analysis Program (弹塑性反应谱及单自由度非线性地震分析工具)

[46]. [Tool][软件][更新][性能设计] PPPSP v2020: Pushover Performance Point Solution Program [PPPSP v2020: Pushover 分析性能点求解程序][Based on FEMA 440]

[47]. [软件][地震动][Update] GMP v2020: A tool for Calculating Ground Motion Parameters for Seismic Analysis of Structures [结构抗震分析地震动强度指标/地震动参数计算工具]

[48]. [反应谱][动力学][抗震] 不同阻尼比反应谱曲线的相交现象 (The Curve Intersection Phenomenon of Response Spectra with Different Damping Ratios)

[49]. [地震][动力学] 对称结构的地震剪力规律 (Seismic Shear Law of Symmetrical Structures)

[50]. [软件][地震波][抗震] GMS_Design 地震波选波工具

[51]. [地震][结构] 双向地震作用效应,【先振型组合,再方向组合】及【先方向组合再方向组合】的差异？(实际案例测算)

[52]. [地震计算][反应谱][动力学][CQC] 振型叠加法随着组合振型数量的增加各种响应量是怎么变化的？

[53]. [抗震][动力学] 对于整体结构，X向地震作用下有Y向剪力吗？有！！

[54]. [地震波][选波][教程] 时程分析地震波选波介绍 ( Introduction of Seismic Wave Selection in Time History Analysis)

[55]. [抗震][减震][笔记] 黏弹性阻尼器的减震性能曲线 Response Reduction Curve of Viscoelastic (VE) System

[56]. [地震][动力学][Dynamics] 将阻尼矩阵的非对角线元素取为0计算结果会怎么样？

[57]. [工具][软件][规范] 广东省标准《高层建筑混凝土结构设计规范》反应谱计算工具

[58]. [资料] 1940 El Centro (Imperial Valley) Earthquake 地震波

[59] [动力学][Dynamics][SAP2000] SAP2000中振型向量的标准化方法

[60] [Dynamics][动力学][抗震] 等效地震力与伪加速度反应谱(Equivalent Static Lateral Seismic Force and Pseudo-Acceleration Spectrum)

[61] [Dynamics][动力学][SAP2000] 梁的振动形态及振型质量 (Vibration Modes and Modal Mass of Beams)

[62] [Dynamics][动力学] 绝对加速度大还是相对加速度大？( Which value is greater? Absolute acceleration or relative acceleration?)

[63] [地震动][动力学] 地震波的基线修正是什么意思？什么时候要基线修正？

[64] [SPECTR][SeismoSignal] SPECTR与SeismoSignal反应谱计算有差异？

[65] [动力学][舒适度] 增加结构的刚度可以降低的结构的加速度吗？

[66] [动力学][Dynamics] 动力时程分析求解的位移和速度结果能否通过相应的加速度结果进行积分得到？(Can the displacement and velocity results obtained by dynamic time history analysis be obtained by integrating the corresponding acceleration results?)

[67] [软件][动力学][Dynamics] NSDOF算例1——单自由度体系弹性动力时程分析

[68] [软件][动力学][Dynamics] NSDOF算例2——单自由度体系非线性动力时程分析

[69] [软件][动力学][Dynamics] NSDOF算例3——非线性粘滞阻尼单自由度体系动力时程分析

[70] [软件][动力学][Dynamics] NSDOF算例4——摩擦阻尼单自由度体系动力时程分析

[71] [软件][动力学][Dynamics] NSDOF算例5——非线性粘滞阻尼器+材料非线性 单自由度体系动力时程分析

[72] [软件][动力学][Dynamics] NSDOF算例6——设置黏弹性阻尼器单自由度体系动力时程分析

[73] [下载][软件]GB-SPECT V2021: 中国规范反应谱生成程序 [Chinese Code’s Design Response Spectrum]

[74] [结构][软件] TLDPC: 调谐液体阻尼器参数计算器 [TLDPC: Tuned Liquid Damper (TLD) Parameter Calculator]

[75] [下载][软件] 黏滞阻尼系数单位转换工具 [Viscous Damping Coefficient Unit Conversion Tool]

[76] [笔记][算例][减振] 调谐液体阻尼器减震算例 [Example of Vibration Reduction of Tuned Liquid Damper (TLD) ]

[77] [动力学][地震] 振型分解反应谱法构件地震力的计算过程？

[78] [下载][软件][地震工程] Spectr_Evolution: Seismic response spectrum evolution [地震波反应谱演化][反应谱随积分时间长短的变化]

[79] [软件][数学][地震动] FPSA: Fourier and Power Spectra Analysis [地震波频谱分析工具]

[80] [动力学][地震工程] 一个有趣的问题: SPECTR中的Newmark-Beta法计算反应谱发散？

[81] [结构设计][动力学] YJK中CQC振型组合地震力的复核

[82] [工具][软件][地震动] AEEG: A Program for Artificial Earthquake Accelerograms Generation [人工地震波合成软件]

[83] [论文][Paper] 基于目标谱匹配法的地震波选波系统研制 (Development of seismic wave selection system based on target spectrum matching method)

[84] [动力学][地震动] SPECTR与SeismoSignal反应谱计算有差异？

[85] [软件][动力学][编程] NMDOF v2022: A Tool for Dynamic Analysis of Nonlinear Shear-Type MDOF System (多自由度剪切层模型系统动力非线性分析工具 v2022)

[86] [动力学][软件] NMDOF算例1 —— 单自由度体系(WEN模型)非线性动力时程分析 [Nonlinear Dynamic Time History Analysis of Single Degree of Freedom System (WEN Model)]

[87] [动力学][软件] NMDOF算例2 ——两自由度简化隔震结构体系(WEN模型)非线性动力时程分析 [Nonlinear Dynamic Time History Analysis of Two Degrees of Freedom Simplified Isolation Structural System (WEN Model)]

[88] [软件][更新][Dynamics] NSDOF v2022: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2022)

[89] [软件][地震动] GMP v2022: A tool for Calculating Ground Motion Parameters for Seismic Analysis of Structures [结构抗震分析地震动强度指标/地震动参数计算工具]

[90] [科研][工具][地震动] RSF Response Spectrum Fitting v2022: 反应谱拟合及反应谱特征参数提取工具 [RSF: A tool for fitting response spectrum and extracting response spectrum parameters]

[91] [抗震][动力学] 软件是如何计算偶然偏心地震作用的? (How does the program calculate the accidental eccentric earthquake action?)

[92] [软件][地震动] AEEG人工地震波合成 —— 案例1(拟合自定义反应谱)

[93] [软件][地震动] AEEG人工地震波合成 —— 案例2(如何提高人工波拟合精度)

[94] [Tool][Seismic Design] GDRSC: A Tool to Generate Design Response Spectral Curve for ASCE 7-16 [ASCE 7-16 美规设计反应谱曲线生成器]

[95]. [软件][地震工程][科研][更新] IRSA 2022: Inelastic Response Spectra Analysis Program (弹塑性反应谱及单自由度非线性地震分析工具)

[96]. [软件][地震工程] IRSA 案例1 —— 等延性反应谱分析算例 ( Constant Ductility Response Spectra Analysis Examples of IRSA)

[97]. [软件][地震工程] IRSA 案例2——单自由度非线性地震分析算例 ( Inelastic SDOF Earthquake Analysis Examples of IRSA)

[98]. [软件][地震工程] IRSA 案例3 – 等屈服强度系数延性需求谱分析算例 ( Equal yield strength coefficient Ductility Demand Response Spectra Analysis Examples of IRSA)

[99]. [软件][地震工程] IRSA 案例4 – 等延性需求谱算例验证 ( Verification Examples of Constant Ductility Demand Spectrum for IRSA )

[100]. [软件][地震工程] IRSA 案例5——地震波能量谱分析 ( Seismic Wave Energy Spectrum Analysis Examples of IRSA)

[101]. [软件][地震工程] IRSA 案例6——地震波基线修正 ( Seismic Wave Baseline Correction Examples of IRSA)

[102]. [地震工程][动力学] 地震作用下结构的能量分析 [Energy analysis of structures under earthquake]

[103]. [动力学][地震工程][工具] SDOF_FRE: Dynamic response analysis of SDOF system using frequency domain analysis method [单自由度体系动力响应的频域分析工具]

[104]. [编程][软件][地震动] PulsePeriodExtrat: A tool for extracting the pulse period of near field vibration velocity [近场地震动速度脉冲周期提取工具]

• 微信公众号 ( Wechat Subscription)

欢迎关注 “结构之旅” 微信公众号