[动力学][结构] 大震弹塑性顶点位移时程为何”不收敛”?[Why does the displacement time history of earthquake elastic-plastic analysis not converge?]

实干、实践、积累、思考、创新。 经常碰到小伙伴做完大震弹塑性分析后,跑过来问为何结构顶点位移总是不收敛? 怎么个不收敛法?大家看看下面这个图: 上述两组图中,上方的是结构某主方向的地面加速度(总时间110s),下方的是结构对应方向的顶点位移时程曲线(算到60s)。 由图可见:其中地面峰值加速出现在30~40s,在60s时地面加速度已退化为峰值的30%以上,而结构顶点位移算到60s依然不减衰减。地面加速度衰减速度很快,加速度峰值明显靠前,而顶点位移似乎还没出现峰值,看起来“发散”。为何加速度已经显著退化,而位移还没收敛? 最初看到这类曲线的时候,也很诧异?后面思考后发现,之所以会存在这种诧异,是因为我们对比的基准选错了。 上面的例子中,我们拿结构的顶点位移时程和地面加速度进行对比?实际上,加速度和位移之间差了两次积分,加速度峰值和位移峰值并不一定出现在同时刻,两者本身可以差很大。 对比结构顶点位移时程的趋势参照地面的位移时程更直观,而不是参照地面加速度时程。位移与位移对比才直接。 以上图中的加速度为例,我们对加速度进行积分,获得对应速度与位移时程,结果如下图所示。 由上图可见,X向和Y向的地面加速度时程峰值分别出现在30s和40s左右,而对应的位移时程峰值分别出现在50s和55s,足足推后了15~20s左右。 从这个角度来看,顶点位移要出现明显退化,加速度应该算到60s,甚至可能更多。为此我们把弹塑性分析时间直接设置为110s的时间,再次提取顶点位移时程结果,如下图所示。 有上图可见,在60s后,两个方向的顶点位移均开始出现不同程度的退化。由于我们通常进行的是一致地震激励计算,因此超限报告中,往往只给出加速度时程曲线,并没有提供地面位移时程曲线,让人很自然地采用地面加速度时程对比顶点位移时程进行参照,从而引起了上述诧异。 有了这样的发现,我们不妨找几组天然地震记录进行积分计算,看看这些地震加速度时程曲线的“显著退化点”和地面位移时程曲线的“显著退化点”出现的位置情况。 这里曲线的“显著退化点”定义为曲线正向最大值和负向最大值中靠后出现的那个。 由上述多组图可见: (1)地面加速度时程曲线的抖动程度最大,经过二次积分后得到的位移时程曲线变得平缓,速度时程曲线的抖动程度次之。 (2)地面加速度、速度、位移曲线的“显著退化点”出现位置不同,三个曲线的“显著退化点”出现的位置先后关系不定,但从选取的几组天然地震波的分析结果来看,速度时程曲线和加速度时程曲线的“显著退化点”相对较为接近,而地面位移时程曲线的“显著退化点”与加速度时程曲线的“显著退化点”可能相差较远,且位移曲线的“显著退化点”更靠后。 (3)对比结构顶点位移时程的趋势参照地面的位移时程更直观,而不是参照地面加速度时程。 相关博文( Related Topics) [01]. [Tool] SPECTR – A program for Response Spectra Analysis …

[软件][地震工程] IRSA 案例2 —— 单自由度非线性地震动力时程分析 ( IRSA Example 2— Nonlinear Seismic Dynamic Time History Analysis of SDOF System)

实干、实践、积累、思考、创新。 外国网友让我对 IRSA 软件做个单自由度地震时程分析分析的使用说明和算例,所以下面较多英文为主,大概看看。关于这个软件的介绍可以访问这个链接:[软件][地震工程][科研][更新] IRSA 2022: Inelastic Response Spectra Analysis Program (弹塑性反应谱及单自由度非线性地震分析工具) [Step 1] 导入地震波 Import Earthquake Acceleration Ground Motion Records 点击“Import Ground Motion”导入地震波,IRSA支持一次导入多组地震波。地震波导入后,可以在“”Time Series“选项卡上查看地震波的加速度、速度、位移时程曲线。软件可支持对地震波进行基线修正,如果需要修正,可以在“Type”下拉菜单选择合适的修正方法,然后点击 “Apply Baseline Correction” 按钮进行基线修正。修正后的时程曲线在图形中用灰色显示。如果执行了地震波基线修正,那后续的分析会基于基线修正后的地震波进行。 Click the “Import …

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

新版已推出,请移步:[软件][编程][动力学] NSDOF v2023: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2023) 实干、实践、积累、思考、创新。 趁着国庆,来个大更新,再2020版NSDOF的基础上,增加了非线性黏滞阻尼器、摩擦阻尼器、黏弹性阻尼器、金属阻尼器,适合做减震研究。 程序图标 ( Program Icon )     程序介绍 ( Program Introduction) NSDOF 是一个基于微软的windows窗口程序,用于单自由度结构的动力非线性分析。结构可是弹性也可以是弹塑性。动力荷载可以是施加在结构基座的地震加速度,也可以是施加在结构顶部的动力荷载。程序使用逐步积分法求解增量非线性运动方程。可以输出结构的各种响应结果,包括抗力,阻尼力,参考惯性力,位移,速度,加速度,耗能,滞回曲线等。 NSDOF is a …

[软件][动力学][Dynamics] NSDOF算例6——设置黏弹性阻尼器单自由度体系动力时程分析 (Dynamic time history analysis of a single degree of freedom system with viscoelastic dampers)

实干、实践、积累、思考、创新! 粘弹性阻尼器是一种有效的减震控制装置,主要依靠粘弹性材料的滞回耗能,为结构提供附加刚度和阻尼,减小结构动力反应,从而实现减震目标。粘弹性阻尼器既能提供刚度,也能提供阻尼,其典型的滞 回曲线为椭圆形,具有良好的耗能性能。下图 所示是粘弹性阻尼器的滞回环形状示意。 粘弹性阻尼器,恢复力由两部分构成,即速度相关部分和位移相关部分,具体如下 NSDOF提供了单自由度主体结构及黏弹性阻尼器并联的动力时程分析模型,以下做两个测算。 第一个算例,结果部附加任何阻尼器,结构刚度取10,阻尼比取0.05,对应的阻尼系数为0.15915,点击运行计算可获得对应的力、位移响应结果及能量图。由能量图可见,此时主要包含3种能量:动能+应变能+阻尼耗能 第二个算例,主体结构的阻尼取0,刚度取一个很小的值0.0001,几乎可忽略,同时设置黏弹性阻尼器的阻尼系数为0.15915,与算例1中的主体结构的阻尼系数一致,黏弹性阻尼器的刚度取10,鱼算例1中的主体结构的刚度一致。因此可知,算例2的计算结果应该与算例1是一致的,相当于用一个黏弹性阻尼器去等效一个单自由度系统。计算结果如下图所示。 由结果可见,算例2的位移、速度、加速度响应均是与算例1一致的。不同的是,算例2的能量图,此时主要包含2种能量:动能+黏弹性阻尼器的耗能,因为用黏弹性阻尼器等效算例1的主体结构,同时也可以发现,黏弹性阻尼器的耗能等于算例1主体结构的应变能及阻尼耗能的叠加。 另外从滞回曲线可以发现,黏弹性阻尼器的滞回曲线是椭圆。此时黏弹性阻尼器的滞回曲线也等于算例1主体结构他弹性恢复力滞回曲线及阻尼力滞回曲线的叠加。 相关话题 ( Related Topics) [01]. [Tool] SPECTR – A program for Response Spectra Analysis [反应谱计算程序] [02]. [程序][Tool] Ground Motion Selection [强震记录选取] [03]. [程序][软件]Ground Motion …

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

实干、实践、积累、思考、创新。 介绍 NSDOF ( [软件][更新][Dynamics] NSDOF v2021: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2021)  ) 设置摩擦阻尼器的单自由度体系动力时程分析。简单来个step by step图片演示吧。 STEP 1: 导入一个震荡动力荷载 STEP 2: 假定主体结构为弹性,设置摩擦阻尼器的摩擦力及刚度。 STEP 3: 点击Run …

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

实干、实践、积累、思考、创新。 在网友建议下,NSDOF (  [Tool][软件][Dynamics] NSDOF v2020: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2020)  ) 增加了非线性粘滞阻尼器。可以在考虑或者不考虑结构阻尼的情况下,考虑非线性粘滞阻尼器进行动力时程分析。下面算例测试NSDOF同时设置非线性粘滞阻尼器+材料非线性的动力时程分析功能,同时用SAP2000进行同样的分析,并对比验证。其中非线性阻尼器阻尼指数取0.1,整个动力方程高度非线性。 算例参数 质点质量: 1kg 体系弹性刚度:100N/m 屈服强度: 0.75N 屈服后刚度强化系数:0 结构的粘滞阻尼系数c:0N-s/m 粘滞阻尼器的阻尼系数cvd: 1.0 N-s/m …

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

实干、实践、积累、思考,创新。 在网友建议下,NSDOF (  [Tool][软件][Dynamics] NSDOF v2020: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2020)  ) 增加了非线性粘滞阻尼器。可以在考虑或者不考虑结构阻尼的情况下,考虑非线性粘滞阻尼器进行动力时程分析。下面做两个算例测试一下NSDOF的非线性粘滞阻尼器计算功能,同时用SAP2000进行同样的分析,并对比验证。 算例1:线性粘滞阻尼 具体参数: 质点质量: 1kg 体系弹性刚度:100N/m 结构的粘滞阻尼取 :0 粘滞阻尼器的阻尼系数c: 1.0 N-s/m 粘滞阻尼器的阻尼指数alpha:1.0 …

[软件][动力学][Dynamics] NSDOF算例2——单自由度体系非线性动力时程分析 ( NSDOF Example 2 — Nonlinear dynamic time history analysis of single degree of freedom system)

实干、实践、积累、思考、创新。 最近小伙伴做非线性粘滞阻尼器的参数分析,于是我们在 NSDOF (http://www.jdcui.com/?p=13947)软件上加了非线性粘滞阻尼器的分析功能。顺便做些测算例子。 这个例子与前面[软件][动力学][Dynamics] NSDOF算例1——单自由度体系弹性动力时程分析 例子的模型基本一致,不同之处在于此例结构为非线性,取二折线非线性本构,进行非线性动力时程分析。 单自由度体系参数: 质量m: 1.0 N-s2/m (kg); 阻尼比: 0.05; 初始刚度k0: 10 N/m; 屈服强度Fy: 0.75N 相应的屈服位移为 0.075m 重力加速度g: 9.807 m/s2; 加速度时程: A9OL 对应的阻尼系数c: 0.31623 N-s/m, 单自由度体系的自振周期为  …

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

实干、实践、积累、思考、创新。 最近小伙伴做非线性粘滞阻尼器的参数分析,于是我们在 NSDOF (http://www.jdcui.com/?p=13947)软件上加了非线性粘滞阻尼器的分析功能。顺便做些测算例子。 单自由度体系参数: 质量m: 1.0 N-s2/m (kg); 阻尼比: 0.05; 弹性刚度k: 10 N/m; 重力加速度g: 9.807 m/s2; 加速度时程: A9OL 对应的阻尼系数c: 0.31623 N-s/m, 单自由度体系的指针周期为  1.98692s 采用NSDOF进行计算,设置参数,并分析 同时采用NONLIN进行计算,并对比验证。 时程结果对比 NONLIN的结果 NSDOF的结果 滞回曲线结果对比 …