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

实干、实践、积累、思考、创新! NMDOF是一个多自由度剪切层模型动力非线性分析工具。NMDOF软件链接:[软件][动力学][编程] NMDOF v2022: A Tool for Nonlinear Dynamic Time History Analysis of Shear-Type MDOF System (多自由度剪切层模型系统动力非线性时程分析工具 v2022) 本篇博文用NMDOF做一个单自由度非线性动力时程分析案例,采用模型是Plastic Wen模型,也就是常说的Bouc-Wen模型,并同时利用SAP2000进行对比验证。 (1)模型参数 质量【MASS:100】 材料模型【MATERIAL:Plastic Wen】 Stiffness: 80000 Yield Strength: 640000 Post …

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

软件已更新,新版移步:http://www.jdcui.com/?p=22688 实干、实践、积累、思考、创新。 拖了好久的东西,这次总算可以更新了。有朋友找到我们,让我们做一个隔震层模型计算程序。于是有了这个剪切层模型的非线性动力分析程序。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) NMDOF 是一个基于微软的windows窗口程序,用于多自由度剪切层模型系统的地震动力非线性分析。 结构可是弹性也可以是弹塑性,支持的材料模型包括,线弹性材料、常规的二折线滞回本构及Bouc-Wen滞回本构。 程序使用Newmark-beta逐步积分法求解增量非线性运动方程。软件提供的阻尼模型包括模态阻尼及瑞利阻尼。 软件可计算并显示结构的模态形状。 软件可输出结构的各种时程响应结果,包括位移,速度,加速度,各类耗能时程等。同时软件可输出多种层最大响应结果,包括层位移、速度、加速度及层间剪力等。 软件可显示结构的模态形状动画及时程变形动画。 软件提供多种常用地震加速度时程格式模板,方便使用者快速导入地震加速度时程,形成自身的地震加速度记录数据库。 NMDOF is a Microsoft-based Windows program for seismic dynamic nonlinear analysis of …

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

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

[Tool][动力学][Dynamics] RDCC v2020: Rayleigh Damping Coefficients Calculator [瑞利阻尼计算器]

实干、实践、积累、思考、创新。 随后更新……           微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[PERFORM-3D] PERFORM-3D中阻尼怎么取值?

实干、实践、积累、思考、创新 读者在看《PERFORM-3D原理与实例》(http://www.jdcui.com/?page_id=3757)一书时问我这个问题,今天抽时间整理了一下,如果我回答有问题,欢迎大家批评指正。 读者问题:“PERFORM-3D中阻尼参数怎么取?常用的结构体系如钢筋混凝土结构、钢结构、型钢混凝土结构的阻尼比怎么取值?“ 答: (1)p3d中,建议用模态阻尼+瑞丽阻尼的组合方式进行动力分析,你刚才界面显示的是瑞丽阻尼。为何这样处理?因为我们小震分析用的是反应谱法,各振形指定阻尼,这个阻尼组装成矩阵对应的阻尼就是模态阻尼。而p3d中,也提供了这个选项,但完整的模态阻尼阶次是与动力自由度数量一致的,除非把所有模态求出来,才能组合成完整的模态阻尼,而实际弹塑性动力分析的时候不可能求太多阶模态,因为计算效率问题,所以就取前面主要的模态数量,而放弃后面的模态,比如p3d模态分析的最多好像只允许60阶,那舍弃掉的就是部分高频阻尼,因此,开发者建议,这里仅是象征性增加一点瑞丽阻尼,如上图的0.1%,同时有利于收敛性,有利于收敛性。具体取值对总体结果影响不大。 (2)三种结构阻尼比怎么取,按整个结构取一个统一阻尼比其实是最早的提法,整体结构笼统一般仅适用于纯质材料,如常规混凝土结构取0.05,纯钢结构0.03,可能都没太多异议。但到了混合结构怎么取,就头疼了,因为分不清材料多少。其实目前大部分弹性分析软件,包括YJK,midas Gen 都提供了基于材料应变能的阻尼比计算方法,你直接指定材料的阻尼比,软件会基于材料应变能算出各阶振型的阻尼比。对于混合结构,你会得到各阶振型下的阻尼比,然后各阶振型按不同的阻尼比直接进行反应谱计算。 (3)因此回到实际情况,倘若你p3d做的是混合结构,建议你在yjk中用材料应变能进行分析,得到前面主要各阶模态的阻尼比,然后再取一个合理值,在p3d中通过模态阻尼比进行指定。因为貌似p3d不能按振形分别指定阻尼比,只能所有模态统一指定一个阻尼比。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[结构设计][地震作用][规范]振型分解反应谱法的一些概念总结 (Basic Concepts of Response Spectra Method)

实干、实践、积累、思考、创新。 温故而知新,理论指导实践,实践检验理论。 (1)振型型分解法,首先是进行模态分析,有多少个动力自由度,理论上就有多少个模态,相应的有多少个周期(频率),及振型。 (2)振型向量关于质量矩阵及刚度矩阵正交。因此,无阻尼运动方程可以实现解耦,将耦合的运动方程,解耦为多个广义单自由度运动方程。 (3)如果阻尼矩阵也满足于振型的正交性条件(如,瑞丽阻尼),则有阻尼结构的运动方程也可以解耦,解耦为多个有阻尼的广义单自由度运动方程。 (4)解耦后的单自由度方程的频率就是振型的频率。即,看是错综复杂的多自由度的震动过程其实是多个规则的不同频率的三角函数组成的。(PS. 自然界就是这么神奇,就像傅里叶变换一样,看是动态的,实则背后是静态的,是死的,太可怕了,无规律的东西,从频率来看,背后却是规律的… 这里不扯这个。 (5)振型无绝对大小,只是表示结构按某个具体频率振动时,各个动力自由度的振幅的相对大小。 (6)利用振型将多自由度方程解耦后分,若对解耦的单自由度方程进行时程分析,该方法常称为模态时成分析方法。若对解耦的单自由度方程进行反应谱分析,则称为振型分解反应谱法,这是目前结构设计规范的主流设计方法。 (7)振型分解法依靠振型对运动方程进行解耦,而振型是与弹性刚度及质量相关的,因此,机遇固定的振型对运动方程解耦,也意味着结构必须是弹性,该方法仅适用于弹性分析。 (8)振型分解反应谱法,由于引入了反应谱,使得结构工程师主要关注最大值,查看结果简便了,但是简便也带来了问题,因为反应谱丢掉了时程结果的许多信息。 (9)由于反应谱法只能获得最大值,因此振型分解反应谱法涉及多个层次的组合问题。首先,各振型的极大值怎么叠加组合为最后的响应,该部分组合是所谓的“振型组合”,如常见的ABS组合方式,SRSS组合方式,及CQC组合方式等。另外,还有一个组合问题是多个方向的地震响应的组合问题,由于不同方向的地震动严格来说是不同的,所谓的不同,是说具体的时程肯定是不同的,响应的反应谱也是不同的,不同就会导致不同步,不同步那不同方向的结果也需要组合。直接时程分析法考虑多个方向的地震同时作用,直接就把多个方向的地震波加上同时进行分析即可,无非是动力方程的右边项将不同方向的地震波叠加即可,而振型分解反应谱法不行,不同方向的地震响应结果,也需要组合,先进行单个方向的效应分析,然后再把这些单个方向的极大值效应进行组合,该组合即所谓的“方向组合”。 (10)由于振型分解反应谱法的概念是,先计算单个振型的某个效应(如剪力,弯矩等)的最大值(正值),然后将单个振型的结果按一定的方法叠加起来,因此,振型分解再用反应谱分析再叠加的过程,丢掉了方向性。或者说,这些响应量,如剪力,只有一个统一的方向。结果都只有一个方向,那使用起来不方便,不直观,因此,在应用的时候,为了给出方向,又有研究者给出一些建议方法,判定响应方向,比如按主振型的方向,来确定响应的方向。但该方法也仅是对于一些简单结构,给出一个响应的参考方向,对于复杂结构,依然存在问题一些问题。比如,多塔连体结构,由于振型分解反应谱法,给出的不同塔楼的力都是同一个方向的,那振型分解反应谱可能就丢失了塔楼的反向运动,有可能存在隐患。因此,振型分解反应谱法虽然简便好用,但是也有不足,这个时候就需要补充弹性时程分析。这就是为何规范要求对复杂结构进行补充的弹性时程分析的一个重要原因。这个振型分解反应谱法的方向问题,还会引起其他相关的问题,具体工程的时候具体思考和分析。 (11)振型分解反应谱法的振型组合是非线性的,因此会出现诸如振型分解反应谱法的楼层剪力与楼层地震力(外力)不平衡的问题。因为,楼层剪力是多个振型的楼层剪力组合而得到的,单个振型下的楼层剪力是由于楼层地震力根据平衡求解的,满足平衡关系,但是经过振型组合后(如,SRSS,CQC),又不满足平衡关系了,因为这些振型组合的方法都不是线性的。因为不能是线性的,为何振型组合不是线性?简单说,因为各个振型的极大值不是同时出现的,这个非线性就是考虑了概率性,背后理论是随机振动。 (12)振型分解反应谱法,实际上是一个等效静力分析,为何这么说,因为运动方程经过解耦,再套上反应谱法,对于每一个振型,相当于在各个动力自由度上加上了一个等效惯性力,然后用这个惯性力进行静力分析,得到该振型下相关的响应量,如构件剪力,弯矩,轴力等,然后再进行振型组合。因此,在有限元求解上,其实是一个静力的求解分析过程。 (13)说到振型分解反应谱法,《高规》及《抗规》,又要扯到“扭转耦联”这个四个字,规范也给出了,两个方法,其中第一个是 a.不考虑扭转耦联的振型分解反应谱法,及b.考虑扭转耦联的振型分解反应谱法。其中,不考虑扭转耦联的振型分解反应谱法采用的是 SRSS组合,仅考虑一个水平方向的振型,即仅进行一个方向的振型分析,不考虑另一个方向质量或扭转惯量的耦合作用。考虑扭转耦联的振型分解反应谱法采用的是 CQC组合(CQC,组合过程中各个振型也是耦联的,需要通过两两振型的周期比及阻尼比参数来计算),分析过程中每一个楼层考虑水平方向及扭转方向3个自由度,振型也包含三个方向的分量。 (14)关于“扭转耦联”,必须说的是,由于一般结果,质量中心及刚度中心很难完全重合,因此,结构的扭转振动总是存在的,因此,进行“考虑扭转耦联的振型分解反应谱法”是相对更精确的。 (15)另外,进行“考虑扭转耦联的振型分解反应谱法”分析与是否考虑多向地震作用或者考虑哪个方向地震作用无关。不要将扭转耦联等同于双向或者三向地震作用,不考虑耦联等同于单向地震。考虑扭转耦联,本质上说的是模态分析的时候,需要考虑平动与扭转自由度的耦联,模态需要能反应扭转的成分。是否考虑多向地震作用,只是振型分解后,方向组合的问题。是否考虑不同角度的地震只是涉及到振型参与系数的计算方法。振型分析,是否考虑扭转耦联或者不考虑扭转耦联,仅仅是结构固有特性的反映。 PS. 最后几点对SRSS和CQC及“扭转耦联”的表述还不是太清楚,借筑信达 李楚舒李总 的话补充一下:完全对称(没有扭转)的SRSS和CQC的结果也有较大区别,SRSS会在地震方向低估作用,而在另一个方向高估(见Wilson一书)。所以用CQC与结构是否扭转没关系,而是振型间存在耦合这一客观存在,所以必须用CQC。所以抗规的“扭转耦联”不对,应该是“振型耦联”——这误导了很多工程师. 相关博文( Related Posts ) [01] …

[下载][软件][动力学]RDCC: Rayleigh Damping Coefficients Calculator [瑞利阻尼系数计算器]

分享一个小程序,该程序用于计算瑞利阻尼系数。软件是免费的。如果你发现软件有bug或者软件使用有问题,请联系我!邮箱:jidong_cui@163.com. A small program to share with you. The program is for Rayleigh Damping Coefficients Calculation.This app is free . If you have found any bug in the program or have any problem when using the program, please don’t hesitate to contact me directly. Email : jidong_cui@163.com