[软件][动力学][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算例4——设置摩擦阻尼器单自由度体系动力时程分析

坚持实干、坚持积累、坚持思考,坚持创新。 随后更新。。。。           相关话题 ( Related Topics) [01]. [Tool] SPECTR – A program for Response Spectra Analysis [反应谱计算程序] [02]. [程序][Tool] Ground Motion Selection [强震记录选取] [03]. [程序][软件]Ground Motion Library [强震记录管理] [04]. Artificial …

[软件][动力学][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 (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的结果 滞回曲线结果对比 …

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

坚持实干、坚持积累、坚持思考,坚持创新。 随后更新。。。。           相关话题 ( Related Topics) [01]. [Tool] SPECTR – A program for Response Spectra Analysis [反应谱计算程序] [02]. [程序][Tool] Ground Motion Selection [强震记录选取] [03]. [程序][软件]Ground Motion Library [强震记录管理] [04]. Artificial …

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

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

[Dynamics][动力学] 绝对加速度大还是相对加速度大?

坚持实干、坚持积累、坚持思考,坚持创新。 今天和小伙伴讨论问题。突然说到绝对加速度和相对加速度。小伙伴潜意思认为绝对加速度大于相对加速度。 因为,一致激励地震动力方程分析的时候,有动力方程可知,实际计算获得是结构的相对加速度,结构的绝对加速度等于相对加速度加上地面加速度。这么一听,似乎绝对加速度比相对加速度要大。 是不是这样呢?以下用NSDOF(   [Tool][软件][Dynamics] NSDOF v2020: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2020)  )做几个线性单自由度系统的时程分析案例。 例子1: 例子2: 从上面两个例子看,结构的相对加速度和地面加速度不总是同向的,绝对加速度可以比相对加速度大,也可以比相对加速度小,与结构的刚度、阻尼等参数有关。 其实我们可以这么想:当结构无限刚的时候,结构相对加速度为0,绝度加速度等于地面加速度,当结构无限柔的时候,结构的相对加速度与地面加速度反向,绝对加速度等于0。因此,当结构刚度K介于0~∞之间时,绝对加速度可能大于相对加速度也可能小于相对加速度。如下图所示: 相关博文( Related Posts ) [01] [Structural Dynamics][Mode …

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

坚持实干、坚持积累、坚持思考,坚持创新。 《建筑抗震设计规范》(GB50011-2010)中给出了采用振型分解反应谱法计算地震作用时的地震力计算公式:\({F_{ji}} = {\alpha _j}{\gamma _j}{X_{ji}}{G_i}\),其中\({\gamma _j} = \frac{{\sum\limits_{i = 1}^n {{X_{ji}}{G_i}} }}{{\sum\limits_{i = 1}^n {X_{ji}^2{G_i}} }}\),\({F_{ji}}\)为j振型i质点的水平地震作用标准值;\({\alpha _j}\)为相应于j振型自振周期的地震影响系数;\({X_{ji}}\)为j振型i质点的水平相对位移;\({\gamma _j}\)为振型的参与系数。以下根据结构动力学的相关理论,给出上述公式的一种推导。 1多自由度体系振型分解法 Mode Superposition Method 对于多质点体系,地震动力方程为: $${\left[ M \right]\left\{ {\ddot u} \right\} …

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

坚持实干、坚持积累、坚持思考,坚持创新。 最近研究舒适度,做些算例测算。对两端铰接、两端固接、一端固接一端铰接、悬臂等截面梁进行振型分析,获得各类梁的前三阶振型,并对振型向量进行最大位移值归一化,并利用归一化后的振型向量求解前3阶振型的振型质量。测试算例梁截面统一为,梁截面为200X200,梁长度为1000mm,沿梁长划分80个单元,振型的质量通过公式  \({M_n} = \int_0^L {m(x)\phi _n^2(x)dx} \) 进行计算。 1 简支梁 1.1 振型形状 一阶振型 二阶振型 三阶振型 1.2 振型质量 振型 节点质量 总质量 振型质量 振型质量/总质量 1 0.0009815 0.314 0.1570 0.500 2 0.0009815 …

[Dynamics][动力学] 振型向量归一是否对计算结果有影响?

坚持实干、坚持积累、坚持思考,坚持创新。 题目如题,结论肯定是没有影响的,因为振型向量本来就是不定的,振型元素之间只有相对关系,要求解振型向量元素的具体值,必须对振型向量进行标准化。简单说即先假定某个元素的值,然后才能求解出其余元素的值。 最近在研究舒适度,顺便把相关东西整理一下,正好还有小伙伴问,同时正好测试一下在网站上用LATEX写公式,看看是不是会专业点。 基本公式 结构的运动方程: \[[M]\{ \ddot u\} + [C]\{ \dot u\} + [K]\{ u\} = \{ P\}  (公式1) \] 将位移向量\(\{ u\} \) 用振型展开, \[\{ u\} = [\phi ]\{ q\} …

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

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 最近研究舒适度,涉及振型向量的标准化,顺便测试一下SAP20000默认的振型向量标准化方法。 大家都知道,振型向量是不定的,振型向量的参数之间只有比值关系。为了求解振型向量的元素绝对值,必须对振型向量进行标准化。 我们接下来测试SAP2000中振型向量的标准化方法,在SAP2000中剪力一根简支梁,如下图: 梁的前三阶Z向振型形状如下: 振型形状是与理论分析结果一样的。 将软件输出的振型变形的平方乘以节点质量,可获得各振型的广义质量,结果均为1。即SAP2000默认输出的振型是满足关于质量矩阵内积为1的条件的。即采用的是关于质量矩阵的正交归一化方法。 相关博文( Related Posts ) [01] [Structural Dynamics][Mode superposition] 振型参与质量系数(Participating Mass Ratio) [02] [动力学][振型分解][Mode Superposition] 振型向量与振型参与系数的乘积公式推导 [03] [结构设计][地震作用][规范] 振型分解反应谱法的一些概念总结 (Basic Concepts of Response Spectra …

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

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 如题,一看是一个莫名其妙的想法。实际上也是一个错误的想法,不过,当时脑子一热,就测算一下。具体看看结果。 模型: 20层的剪切层模型。 阻尼矩阵: (1)模态阻尼,20阶振型计算 (2)模态阻尼,20阶振型计算,把非对角线元素取为0。 采用MATLAB编程,采用Newmark-β积分法进行弹性时程分析,两种阻尼模型的计算结果对比如下。 其中,参考阻尼为完整的模态阻尼,对比阻尼为去掉阻尼矩阵对角线元素后的矩阵。 由图可见,采用仅保留对角线元素的阻尼矩阵,结构的位移及剪力响应远小于完整的阻尼矩阵,而楼层加速度响应似乎相差不大!!十分诧异!! 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[地震][动力学] 对称结构的地震剪力规律

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 来自小伙伴 刘骥 的分享。这个是继《 [抗震][动力学] 对于整体结构,X向地震作用下有Y向剪力吗?有!! 》后对对称结构进行的简单测试。直接看测试结果吧。 算例1: 算例2: 可见,剪力作为矢量,是满足平行四边形法则的。同时,对于对称结构,往哪个方向输入地震加速度,总剪力均一致。 相关博文 ( Related Topics) [01] [地震][动力学] 对称结构的地震剪力规律 [02] [地震][结构] 双向地震作用效应,【先振型组合,再方向组合】及【先方向组合再方向组合】的差异?(实际案例测算) [03] [地震计算][反应谱][动力学][CQC] 振型叠加法随着组合振型数量的增加各种响应量是怎么变化的? [04] [抗震][动力学] 对于整体结构,X向地震作用下有Y向剪力吗?有!! 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[抗震][动力学] 对于整体结构,X向地震作用下有Y向剪力吗?有!!

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 对于整体结构,X向地震作用下,结构有Y向的剪力吗?以下通过两个简单的时程分析案例进行测算。 案例1: 振型结果如下: 振型1为Y向平动 振型2为X向平动 振型3为绕Z轴扭转 案例2:案例1模型逆时针旋转45度 周期值与模型1是一样的,只是因为结构转了一个角度,振型方向不同了。 振型1为135度方向平动 振型2为45度方向平动 振型3为绕Z轴的扭转 分别对两个模型沿X向施加地震动加速度时程,进行直接积分动力时程分析。所选的地震波如下图所示。 模型1沿X方向与Y方向的基底剪力结果如下图所示。由图可知,对于模型1,沿X向输入地震,Y向剪力几乎为0。 模型2沿X方向与Y方向的基底剪力结果如下图所示。由图可知,对于模型2,沿X方向输入地震,结果Y向会产生剪力,且剪力大小不可忽略。 粗看似乎有点难理解,外力和内力不是应该平衡的吗?为何施加X方向加速度,结构有Y向的剪力?对于静力情况下,结构受到到的外力与的内力平衡,比如,当沿结构X方向施加力F时,结构总的剪力必然是沿X方向,且大小为F,Y方向不存在剪力。为何到了动力情况,就不满足这个规律了?不妨看一下两种情况下结构的平衡方程。 结构静力平衡方程: 其中,为结构的外力,为结构的抗力,其中,当只有X向力作用时,即 ,,即结构的抗力也只有X反向的力,y及z向的力为0   结构动力平衡方程(地震): 其中,当只有x向地震时,,即向量在非X向自由度上为0,此时结构的抗力,假设忽略阻尼,结构的抗力为,对比静力平衡下的公式()可见,尽管向量在非X向自由度上为0,但是抗力的右边项不是,而是,即所谓的绝对加速度,其中相对加速度在非X向自由度上不一定为0,当在非X向自由度上存在非0值时,抗力就可以能存在非X向自由度上的力,即对于整体结构,X向地震作用下,结构整体在其他方向也可能存在抗力,包括Y向的剪力。 对于算例1,在X向地震加速度作用下,由于结构基本只有X向的位移,因此Y向的抗力很小,进而Y向剪力很小。对于案例2,由于结构扭转了45度,在X向地震加速度作用下,结构不仅有X向的位移,也有Y向的位移,有Y向的位移,就可能有Y向的抗力及Y向剪力。 由以上分析也可发现,引起动力与静力概念上不同的错觉的原因是,把动力情况下结构的抗力当成了,实际上动力情况下结构的抗力等于(忽略阻尼情况下)。 平衡还是满足的!!!。 相关博文 ( Related Topics) [01] …

[Tool][软件][Update] PPPSP V2020: Pushover Performance Point Solution Program [Pushover 分析性能点求解程序][Based on FEMA 440]

坚持实干、坚持实践、坚持积累、坚持思考,坚持创新。 7月拖到现在11月,终于可以更新了...... 程序图标 ( Program Icon )        程序介绍 ( Program Introduction) 基于FEMA 440等效线性化法 Pushover分析方法 的性能点求解程序。( A program for the Solution of Pushover Performance Point based on FEMA 440 Equivalent …

[YJK][动力学] 逐步加大结构宽度结构周期的变化算例测算

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 测算10组模型,模型1为2*3网格的5层框架结构,第一周期为X向平动,第二周期为Y向平动,第三周期为扭转。 模型2~模型10,由模型1沿X向不断拉长形成,层数不变。 看看往两端加宽结构的各周期如何变化。 采用YJK MultiModel Compare提取各模型的结果。 绘制各模型的周期变化曲线。 第一周期 X向平动 第二周期 Y向平动 第三周期 扭转 可见,往X向不断加宽,第一周期不断减小,第二周期不断增加,扭转周期不断增加。 为什么会这样?其实 周期的增加或者减小其实与质量、刚度增加的相对比例有关。 这里给出了基本理论解释:[动力学][Structure Dynamics] 线性增加刚度K与质量M下单自由度(SDOF)结构的周期变化 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[动力学][Structure Dynamics] 线性增加刚度K与质量M下单自由度(SDOF)结构的周期变化

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 单自由度体系,线性增加刚度K与质量M,结构周期的变化规律。 现性递增k与m 单自由度结构的周期变化关系:周期可增可减,与k、m各自增幅有关、即 k/m有关。增、减构件类似在原有基础上递增k、m,如果原先K、M的基数已经很大,曲线已经在平滑段,简单线性增减对周期结果影响小。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

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

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。 采用 SPECTR反应谱分析软件 (下载链接: http://www.jdcui.com/?p=1875 ) 对几个地震波进行位移谱的求解,结果如图1~图5所示,可以发现,随着阻尼比的增大,大部分地震波在各周期范围内反应谱值减小。但是部分地震波在某些周期范围内,阻尼比增大,但是位移谱值不一定减小。 在图形上表现为不同阻尼比的反应谱曲线在某些周期范围内出现相交的现象(见图1及图3)。同样,拟加速度反应谱也存在这个现象,因为拟加速度反应谱是由位移谱转换过来的。 初看这个现象感觉很奇怪,其实仔细想也十分正常。 以图1的反应谱为例,在周期2.5s左右出现明显的相交现象。把2.5s周期各阻尼比单自由度结构的位移时程绘制出来,见图6. 由图6可见,各阻尼比下,位移时程曲线的整体趋势是比较一致的,该凸的地方大家一起凸起来,该凹的地方一起凹下去,但是随着阻尼比的增大,出现绝对值最大值的时间不同了,虽然趋势一致。 如图,阻尼比为20%的位移最大值出现在30s左右,其他阻尼比下位移绝对值最大值出现在40s左右。且阻尼比增大到20%后绝对值比其他阻尼比在40s左右出现的绝对值大。因此不同阻尼比的反应谱曲线就出现了交点。 因此,不同阻尼比的反应谱曲线是可能出现交点的。因为,反应谱是纵坐标是绝对值。阻尼比发生改变,可能整个响应时程的整体趋势没改变,但是最大值出现的位置会不同,大小的增大或减小规律也不同,而反应谱记录的是绝对值。 所以,千万别闭着眼睛说,阻尼比越大,位移越小。瞎说!!! 😎 😀  图 1 图 2 图3 图 4 图 5 图 6 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

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

坚持实干、坚持积累、坚持思考,坚持创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) NSDOF 是一个基于微软的windows窗口程序,用于单自由度结构的动力非线性分析。结构可是弹性也可以是弹塑性。动力荷载可以是施加在结构基座的地震加速度,也可以是施加在结构顶部的动力荷载。程序使用逐步积分法求解增量非线性运动方程。可以输出结构的各种响应结果,包括抗力,阻尼力,参考惯性力,位移,速度,加速度,耗能,滞回曲线等。 NSDOF is a Microsoft Windows based application for the dynamic analysis of single degree of freedom structural systems. The structure …

[地震工程][动力学] 不同方法计算多自由度体系(MDOF)地震作用的对比分析(算例检验基本概念)

实干、实践、积累、思考、创新 检验基本理论及基本概念,多自由度地震作用计算的算例对比。 感谢 师弟 符东龙 帮忙整理资料。 算例信息(Model Information) 设计一个多自由度模型,用梁、柱搭一个多层的框架结构即可,结构的第一周期控制在3~4s左右,选取一组地震波(可以用前面的chichi波),阻尼比考虑5%和30%。用ETABS做以下分析。 (1)Mode1:直接输入地震波进行直接积分时程分析 (2)Mode2:直接输入地震波进行模态时程分析,注意(模态数量取满,所以前面设计模型自由度不要设计太多)。 (3)Mode3:直接输入地震波进行模态时程分析,注意(模态数量取1)。 (4)Mode4:直接输入地震波的反应谱(地震波的反应谱可以用SPECTR计算,计算时候周期间隔可以取密一点,比如0.01s),进行振型分解反应谱法分析。如以下是采用本网站的 SPECTR 反应谱分析软件计算的ChiChi.dat的地震波的反应谱,后续也可以用ChiChi.dat这条地震波来计算。 采用的地震波:   任务(Task) (1)阻尼比考虑5%和30%两种进行计算。提取模型1、模型2、模型3的剪力时程、模型4的剪力,并三者对比。 (2)阻尼比考虑5%和30%两种进行计算。提取模型1、模型2、模型3的的位移时程、模型4的位移,并三者对比。 模型(Model) 结果(Results) (1) 剪力对比  Base Force 模型4的反应谱分析结果与模型1和2的有差别。概念自己体会。 (2) 顶点位移 disp …

[地震工程][动力学] 不同方法计算单自由度体系(SDOF)地震作用的对比分析(算例检验基本概念)

实干、实践、积累、思考、创新 检验基本理论及基本概念,单自由度地震作用计算的算例对比。 感谢 师弟 符东龙 帮忙整理资料。 算例信息(Model Information) 设计一个单自由度模型,结构周期控制在 1s,选取一组地震波,阻尼比考虑 5%和 30%,。用 ETABS 做以下分析。 (1) Mode1:直接输入地震波进行直接积分时程分析。 (2) Mode2:直接输入地震波进行模态时程分析。 (3) Mode3:直接输入地震波的反应谱(地震波的反应谱可以用 SPECTR 计算,计算时候周期间隔可 以取密一点,比如 0.01s),进行振型分解反应谱法分析。如以下是用 SPECTR 计算的 ChiChi.dat 的地震波的反应谱,后续也可以用 ChiChi.dat 这条地震波来计算。 …

[动力学][振型分解][Mode Superposition] 振型向量与振型参与系数的乘积公式推导

坚持实干、坚持实践、坚持积累、坚持思考,坚持创新。 复习基本知识,整理资料。振型分解法。 (1)振型向量是有量纲的。量纲为长度的倒数。 (2)振型向量关于质量矩阵正交。 (3)振型向量与振型参与系数的乘积为荷载指向向量。推到过程中利用到分块矩阵的一些计算。 相关博文(Related Topics) [01] [Structural Dynamics][Mode superposition] 振型参与质量系数(Participating Mass Ratio) [02] [动力学][振型分解][Mode Superposition] 振型向量与振型参与系数的乘积公式推导 [03] [结构设计][地震作用][规范] 振型分解反应谱法的一些概念总结 (Basic Concepts of Response Spectra Method) [04] [动力学][SAP2000] SAP2000中振型向量的标准化方法 …

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

忙完工作更新….         注释 ( Comments )   ( 如果您发现有错误,欢迎批评指正。邮箱:jidong_cui@163.com . 如果您喜欢这篇博文,请在上面给我 点个赞 吧! 🙂   🙂      ( If you found any mistakes in the post, please let me know. Email : jidong_cui@163.com. If you like this posts, please …