振型分解反应谱法matlab
基于MATLAB地震反应谱数值算法的稳定性和精度分析
基于MATLAB地震反应谱数值算法的稳定性和精度分析摘要：地震反应谱是进行结构抗震分析与设计的重要工具，反应谱的计算在反应谱法和时域逐步积分方法中有重要地位，引起了学者的重视和广泛研究。
而对计算方法优劣的评定常取决于其计算的耗时、稳定性和精度等因素。
本文基于数值算法的相关研究及应用现状，以MATLAB为平台，建立数值算法在不同影响因素下的三维图形，并结合理论进行对比分析。
通过算例进一步分析验证，得出不同数值算法在实际计算中的表现，为工程实际计算中选取哪种积分算法更为合适提供参考。
关键词：地震反应谱；时域逐步积分算法；稳定性和精度；MATLAB1、地震反应谱的基本假定地震反应谱基于的三个基本假设[1]：(1)结构物所处的地面假定为刚性面，认为体系各质点的运动是完全一致的。
(2)强震观测仪的记录为地面运动的过程。
(3)结构体系不能是双或多质点体系，必须是单质点体系；同时应是弹性体系状态。
这里所谓的单自由度体系结构，就是用无量刚的弹性杆件支承于地面上，将结构体系中参与振动的质量用一点表示。
同时，假定结构振动和地面运动不发生扭转，只是水平平移运动并且是单方向的。
2、基于MATLAB地震反应谱数值算法的稳定性和精度分析2.1 概述目前MATLAB地震反应谱数值理论算法主要有中心差分法[2]、Wilson-法、Houbolt法、线性加速度法及Newmark-法等，理论算法主要是以求解线性结构体系动力方程时所表现出的特性作为数值算法优劣的评价依据[3]，但是在实际工程运用中，人们常常凭借经验来判定选取较为合适的积分方法。
随着工程问题越来越复杂，在对大型复杂结构的结构动力反应分析更为复杂，要求高效率计算情况下获得较精确地计算结果。
然而各计算方法的精度和稳定性对结构动力反应分析的发————————————E-mail:skyuanyan@展引起了很大的影响和制约[4]。
2.2数值算法的稳定性分析基于上述情况，本文对上述几种常用数值算法的稳定性方面通过图形进行比较分析，并结合算例进一步验证分析。
振型分解反应谱
振型分解反应谱振型分解反应谱啊，就像是建筑结构世界里的一个超级神秘的魔法乐谱。
你看那些高楼大厦，一个个就像巨人一样站在那里，可它们要是遇到地震这种调皮捣蛋的“小恶魔”，那可就危险了。
这时候振型分解反应谱就闪亮登场啦，它就像是一个超级英雄，能算出建筑在地震时的各种反应，仿佛是在解读建筑巨人的弱点和抵抗能力。
它就像是一把超级精细的梳子，把建筑结构复杂的振动模式一根一根地梳理出来。
建筑的振动就像是一群人在跳舞，有不同的舞步和节奏，振型分解反应谱呢，能把这些混乱的舞步都分解开，清楚地看到每个舞者（每个振型）的动作幅度和频率。
这就好比是把一场乱糟糟的群魔乱舞变成了一场整齐有序的舞蹈表演，每个舞者都在自己的位置上做着规定动作。
如果把建筑比作一个合唱团，那振型分解反应谱就是那个能听出每个声部声音高低和节奏的指挥家。
它能精确地分析出每个振型（每个声部）的贡献，哪些是高音部（主导振型）在大声歌唱，哪些是低音部（次要振型）在低声附和。
振型分解反应谱还像是一个超级侦探，在建筑结构这个大案件里寻找线索。
它能从地震这个大坏蛋留下的蛛丝马迹（地震动记录）中，找出建筑结构会产生什么样的反应，是瑟瑟发抖得像个受惊的小兔子，还是能够坚强地挺立，像个威风凛凛的大狮子。
你可以想象建筑结构是一个大蛋糕，振型分解反应谱就是一把神奇的刀，把这个蛋糕按照不同的振动模式切成一片片的。
每一片蛋糕（每个振型）都有它独特的大小（反应大小）和味道（振动特性）。
有时候我觉得振型分解反应谱就像一个会魔法的老巫师，拿着他的魔杖（计算方法和公式），对着建筑结构这个魔法阵念念有词，然后就把复杂的振动问题变得清晰明了。
就像把一团乱麻瞬间变成了整齐的丝线。
它也像一个特别严格的考官，建筑结构在它面前就像一个个参加考试的学生。
振型分解反应谱会仔细地检查每个建筑结构（每个学生）在地震这个大考中的表现，看看它们是能够轻松应对，还是会被难倒得一塌糊涂。
再把建筑想象成一艘大船，振型分解反应谱就是船上的瞭望员，提前看到地震这个汹涌波涛（危险）会给船带来什么样的晃动（结构反应），然后好让船员（建筑师和工程师）们提前做好准备，让大船不至于被浪打翻。
基于matlab的地震反应谱计算方法的比较【matlab源码】
毕业论文（设计）题目学院学院专业学生姓名学号年级级指导教师教务处制表基于MATLAB的地震反应谱计算方法的比较一、程序说明本团队长期从事matlab编程与仿真工作，擅长各类毕业设计、数据处理、图表绘制、理论分析等，程序代做、数据分析具体信息联系二、写作思路与程序示例地震反应谱是进行结构抗震分析与设计的重要工具，反应谱的计算在反应谱法和时域逐步积分方法中有重要地位，引起了学者的重视和广泛研究。
而对计算方法优劣的评定常取决于其计算的耗时、稳定性和精度等因素。
目前，计算反应谱的方法有很多，以往常规方法主要有中心差分法、Newmark-法、线性加速度法及Wilson-法等，这些方法虽然存在一些弊端，但是其计算精度能够满足一般实际工程计算的需要，仍被广泛应用于工程实践。
因此有必要对这些方法进行深入的比较和探讨。
目前，我国现行建筑抗震设计规范中对阻尼比规定：混凝土结构通常取0.05，钢结构通常取0.02，阻尼比均比较小，因此利用拟反应谱是可靠的。
然而，随着高层和超高层的出现，建筑新形式、新材料的使用、隔震、减震耗能机构的研究以及抗震减灾新要求的不断提出，传统理论也越来越不能适应新要求，使用拟反应谱可能带来较大误差，这应当引起我们的注意。
同时，通过反应谱理论分析得到：当周期超过3s以后，结构地震反应已不是由地面加速度控制，而可能是由速度甚至是位移控制。
然而，现代的超大跨、超高层和巨型结构的自振周期大都达到了10s以上，如果仍然采用加速度谱来对其进行分析将不合适。
我们应当采用能更好反应三个周期段的反应谱来对长周期结构进行分析，而三联反应谱就具有这一特点，它是将位移、速度以及加速度联合反应谱同时表示在一张图上的四坐标对数图，能够使我们更直接地掌握上述三种物理量与结构自振周期之间的控制关系。
因此，采用MATLAB的GUI编程三联反应谱图形界面，并将其应用于工程实践，将是一项很有意义的工作。
论文基于以上内容，主要进行了以下具体工作：1.基于数值算法的相关研究及应用现状，本文以MATLAB为平台，建立数值算法在不同影响因素下的三维图形，并结合理论对比分析。
matlab 希尔伯特振动分解
Matlab希尔伯特振动分解：深度解析与应用一、介绍Matlab作为一款强大的科学计算软件，常常被用于信号处理、振动分析等领域。
其中，希尔伯特振动分解（Hilbert Vibration Dposition，HVD）作为一种信号处理技术，在振动分析中具有重要的应用。
本文将围绕着Matlab中的希尔伯特振动分解展开详细的讨论，包括其原理、算法、实现以及应用等方面，旨在帮助读者深入理解并灵活运用该技术。
二、希尔伯特振动分解的原理和算法希尔伯特振动分解是基于希尔伯特变换的一种信号分解方法，其核心思想是将原始信号分解为多个振动模态，并分别进行分析。
在Matlab 中，可以利用希尔伯特变换函数hilbert对信号进行HVD分解。
其算法包括以下几个步骤：1. 采集原始振动信号2. 进行希尔伯特变换，得到解析信号3. 提取解析信号的幅度包络曲线4. 通过包络曲线得到高频振动和低频振动成分具体而言，利用Matlab中的hilbert函数可以将原始振动信号进行希尔伯特变换，得到解析信号，然后通过求解解析信号的包络曲线，可以得到原始信号的振动成分。
这一过程可以用以下公式表示：\[x(t) = \frac{1}{\pi} PV \int_{-\infty}^{\infty} \frac{X(\tau)}{t-\tau} d\tau\]其中，\(X(\tau)\)为原始信号，\(x(t)\)为解析信号。
三、希尔伯特振动分解的实现和应用在Matlab中，可以利用上述算法实现希尔伯特振动分解。
以振动信号分析为例，可以将原始振动信号导入到Matlab环境中，然后利用内置的信号处理工具箱函数进行HVD分解。
还可以自行编写希尔伯特变换的代码实现HVD分解，使得分析结果更加灵活和具体。
希尔伯特振动分解在实际应用中有着广泛的领域，包括但不限于机械振动分析、生物医学信号处理、地震学领域等。
在机械振动分析中，可以利用HVD分解对机械系统的故障进行诊断和预测；在生物医学领域，可以对心电图信号进行不同频率成分的分析和诊断；在地震学领域，可以通过HVD分解对地震信号进行时域和频域的分析，以预测地震的发生和研究地震的特性。
振动信号变分模态分解matlab
振动信号变分模态分解matlab振动信号变分模态分解(Matlab)在工程和物理学领域中，振动信号是一种广泛存在的信号类型。
通过对振动信号进行分解可以更好地理解信号，分析信号，并进一步实现信号的相关处理。
其中，变分模态分解（VMD）是一种先进的信号分解技术，可用于将振动信号分解成多个子信号并分析每个子信号。
本文将介绍如何使用Matlab进行VMD分解。
1. 安装VMD软件包VMD分解技术已经被打包为Matlab的一个工具箱，可以通过下载和安装该工具箱来安装VMD软件包。
要开始安装，请浏览VMD软件包的官方网站，下载适用于您的Matlab版本的该软件包。
下载完成后，运行Matlab并将软件包添加到路径中。
接下来，您可以使用Matlab中的VMD函数进行信号分解。
2. VMD分解信号首先，您需要准备一个要分解的信号。
一个简单的方法就是使用Matlab中的随机函数生成一些样本数据进行测试。
示例代码：% Generate sample signalx = randn(1, 1000);% Perform VMD decompositionalpha = 2000;tau = 0;K = 10;DC = 1;init = 1;tol = 1e-7;omega = 0;verbose = 1;[u, ~] = VMD(x, alpha, tau, K, DC, init, tol, omega, verbose);上述代码使用了VMD函数对信号进行分解，其中参数alpha表示VMD算法的控制参数，tau表示信号的初态时间，K表示分解生成的子信号数，DC表示是否保留分解后信号的直流分量，init表示算法的初始条件，tol表示解决非线性优化的容差，omega表示算法的微调选项， verbose表示是否输出调试信息。
3. 分析VMD分解结果一旦分解完成，您可以分析每个子信号以获取信号的相关特征。
例如，您可以使用Matlab中的快速傅里叶变换（FFT）来计算每个子信号的频谱。
基于MATLAB的地震反应谱与傅里叶谱计算分析
参考文献(References)：
[1] 郭晓云. 汶川地震反应谱研究分析[D]. 哈尔滨： 中国地震局工程力 学 研 究 所 ， 2011.(DUO xiaoyun. Study on response spectrum of WenChuan Earthquake.[D]. Harbin ： Institute of Engineering Mechanics，China Earthquake Administration ，2011.(in Chinese)) [2] 骆剑锋. 框架结构静力与动力弹塑性抗震分析对比研究[D]. 上海： 同济大学，2007. [3] 张晓志，谢礼立，于海英. 地震动反应谱的数值计算精度和相关问 题[J]. 地震工程与工程振动，2004，24(6)：15-20. [4] 赵凤新，胡聿贤. 地震动非平稳性与幅值谱和相位差谱的关系[J]. 地震工程与工程振动，2003，23(1)：1-5. [5] 金 星， 廖振鹏. 地震动强度包线函数与相位差谱频数分布函数的
[2]
的方法，因此对Δt 时段内的系统反应或加速度的变 化规律未做出假定。 事实上， 采用 Fourier 变换方法 求解动力平衡方程时，对输入荷载采用了一个全局 性基本假定，即将非周期的加速度荷载假定为一个 具有有限带宽的周期荷载。在此假定下，加速度荷 载被表达为有限个连续简谐波的线性叠加。因此， 频域方法也是近似方法，不能作为度量其他方法精 度的准则。 傅立叶谱不仅包含幅值信息，而且包含相位信 息，这是地震反应谱所不具备的。所以，当傅里叶 谱与地震反应谱进行比较时，只能采用幅值谱。
3 强震记录资料
北岭大地震， 1994 年 1 月 17 日凌晨 4 时 31 分， 洛杉矶地区发生里氏 6.6 级地震，震中位于圣费兰 多峡谷， 北纬 34°12′54″， 西经 118°32′16.8″， 在洛杉矶西北方向 20 英里处， 属浅源地震， 震源深 度 为 18 英 里 。 当 地 震 级 ML=6.6 ， 面 波 震 级 为 MS=6.7（NEIC），地震矩震级为 MW=6.7（CIT）， 持续时间约 30 秒。在持续 30 秒的震撼中，大约有 11000 多间房屋倒塌， 震中 30 公里范围内的高速公 路、高层建筑或毁坏或倒塌，煤气、自来水管爆裂， 电讯中断，火灾四起，直接和间接死亡 58 人，受伤 600 多人，财产损失 300 多亿美元。 本文以北岭大地震中，某县医院台站强震动仪 所记录到的强震记录为例，通过 MATLAB 编程实 现傅里叶谱和地震反应谱的计算， 并进行比较分析。 该台站位于北纬 34°19′33.6″，西经 118°26′ 38.4″，所记录的振动持续时间为 39.98s，峰值加 速度为-1007.942cm/s²，记录时间点在 4.328s；峰值 速度为-124.896cm/s，记录时间点在 3.700s；最大位 移为-31.078cm，记录时间点在 4.200s；初始速度为 -3.694cm/s；初始位移为 298cm。
振型分解反应谱法
附录一 振型分解反应谱法振型分解反应谱法作为弹性多自由体系的主要分析方法，很有必要对振型分解反应谱法有充分的了解。
本文仅作为大家参考之用，如有理解上的错误或者不当，敬请谅解。
1、单自由度体系在地震作用下的运动 如图（1）所示，根据达朗贝尔原理有： 0=++s I c f f f （1）也即：g u m ku u c um -=++ （2） 方程两边同时除以m ，可化为：g u u u u-=++22ωξω （3） 式中，2/k m ω= ，令ωξm c2=，为体系阻尼比。
2、多自由度体系在地震作用下的运动类似于单自由度体系分析过程，体系运动方程为：g u m u k u c u m ][}]{[}]{[}]{[-=++ （4）无阻尼体系自由振动时，0=g u，0=c ，上式即为： }0{}]{[}]{[=+u k um （5） 根据方程解的特征，设其解的形式为：)sin(}{}{ϕωφ+=t u （6）代入（5）式有：}0{)sin(}]){[]([2=+⋅-ϕωφωt m k （7）由于0)sin(≠+ϕωt则}0{}]){[]([2=-φωm k （8）另外，}0{}{≠φ，故特征方程为：0][][2=-m k ω （9）由（9）式可以求出2ω，进而可以求得各阶振型对应的圆频率2i ω，再代入（8）式可求对应于各个2i ω的特征向量}{i φ，即为振型。
振型φ：多自由度体系自由振动时，各质点在任意时刻位移比值是一定的，不随时间变化，即体系自由振动过程中形状保持不变。
振型是结构形状保持不变的振动形式，振型的形状是唯一的。
N 个自由度的体系具有N 个振型。
则结构的变形总可以表示成这N 个振型的线性组合：{}∑==Ni i i q u 1φ （10）其中i q 称为正则坐标。
3、振型的正交性由于}0{}]{[}]{[2=-φωφm k （11） 则}0{}]{[}]{[2=-r r r m k φωφ （12）（12）式两边同时左乘T n }{φ，)(r n ≠，得到：}]{[}{}]{[}{2r T n r r T n m k φφωφφ= （13）同理，}]{[}{}]{[}{2n Tr n n T r m k φφωφφ=，该式两边同时转置一次，得到：}]{[}{}]{[}{2r T n n r T n m k φφωφφ= （14）（13），（14）两式左右对应相减，得到：0}]{[}){22=-r T n n r m φφωω（ )(n r ≠ （15）因为22n r ωω≠所以 0}]{[}{=r Tn m φφ )(n r ≠ （16） 同理亦有 0}]{[}{=r Tn k φφ )(n r ≠ （17）即所说的振型关于质量和刚度矩阵满足正交性质。
八层框架的地震响应计算和人工波生成的matlab实现及所需曲线图的自动存储
一、 作业概况结构基本参数：层间剪切型结构，采用Rayleigh 阻尼，第一、第二阶阻尼比分别取3%、5%。
图1 结构基本形状表1 各层集中质量 ( 105kg)层号 12345678质量表2 各层层间刚度 (×108N/m)层号 1 2 3 4 5 6 7 8 层间刚度m m m m m m m m ()g x t二、 频率及振型计算根据层间模型的假定，可以建立结构的质量矩阵以及刚度矩阵如下。
12345678000000000000000000000000000000000000000000000000000000003.400000000 3.400000000 3.200000000 3.20000 =0000 2.800000000 2.800000000 2.700000000 2.6m m m m m m m m ⎛⎫ ⎪ ⎪ ⎪ ⎪⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎛⎫ ⎝M 510kg ⎪⎪⎪⎪⎪⨯⎪⎪⎪⎪⎪⎪⎭ 11121314151617182122232425262728313233343536373841424344454647485152535455565758616263646566676871727374757677788182838485868788k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k ⎛ =⎝K 8420000002 3.8 1.8000000 1.8 3.6 1.8000000 1.8 3.6 1.8000 =10/000 1.8 3.6 1.8000000 1.8 3.4 1.6000000 1.6 3.2 1.6000000 1.6 1.6N m ⎫⎪⎪⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎭-⎛⎫ ⎪-- ⎪ ⎪-- ⎪-- ⎪⨯ ⎪-- ⎪-- ⎪ ⎪-- ⎪ ⎪-⎝⎭根据上面求得的质量、刚度矩阵，即可求解特征方程：20K M（1）求解自振频率以及阵型向量已经演变成为典型的求解矩阵特征值以及特征向量的问题，特征值即为圆频率2，特征向量即为振型向量。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
%本程序采用振型分解反应谱法计算框架结构水平地震力
%采用KN.M单位
%运行本程序之前请运行CYGD1.M和CYGD2.M求解框架水平侧移刚度%本程序未考虑扭转耦联振动，只能用于平面框架计算。
求解所有振型。
%结构地震影响系数按高规3.3.8选取
%地震作用和作用效应按高规3.3.10计算
clear %清理WORKSPACE
k0=[263770 %各层框架侧移刚度
263770
263770
263770
123582];
m0=[1.904 %各层质量，重力荷载代表值/g
2.677
2.677
2.677
2.677]*1.0e
3./9.8;
n1= 0.21712; %单榀框架地震力分配系数
Tg=0.35; %特征周期(按规范选取)
s=0.05; %阻尼比(按规范选取)
r=0.9; %衰减系数(按规范选取)
y1=0.02; %阻尼比调整系数1(按规范选取)
y2=1; %阻尼比调整系数2(按规范选取)
amax=0.08; %水平地震最大影响系数(按规范选取)
zjxu=0.7 %周期折减系数(按规范选取)
cn=length(m0) %计算楼层数
l=diag(ones(cn));
m=diag(m0); %计算质量矩阵
ik=matrixju(k0,cn); %计算刚度矩阵
[x,d]=eig(ik,m) %求解特征值和特征向量
d=diag(sqrt(d)) %求解结构圆频率
T=zjxu*2*pi./d %求解结构特征周期并作折减，折减系数0.7
for i=1:cn;
[dl(i),j]=min(d);
xgd(:,i)=x(:,j);
d(j)=max(d)+1;
end
w=dl; %输出结构自振频率
x=xgd;
for j=1:cn; %求解结构振型参与系数和各质点的水平相对位移x x(:,j)=x(:,j)/x(cn,j);
zhcan(j)=(x(:,j))'*m*l/((x(:,j))'*m*x(:,j));
y=zhcan' %输出振型参与系数
alpha=zeros(cn,1); %求解地震影响系数
for i=1:cn;
if T(i)<0.1
alpha(i)=5.5*amax*T(i)+0.45*amax;
else
if T(i)<Tg;
alpha(i)=amax;
else
if T(i)>5*Tg;
alpha(i)=[0.2.^r-y1*(T(i)-5*Tg)/y2]*amax;
else
alpha(i)=(Tg/T(i)).^r*amax;
end
end
end
end
a=alpha %输出地震影响系数
F=zeros(cn); %求解各振型下，各层水平地震力for i=1:cn;
F(:,i)=x(:,i).*m0.*y(i).*a(i).*9.8;
f=F.*n1 %求解分配后地震力，并排序z=zeros(cn,cn);
b=cn+1;
for n=1:cn;
z(b-n,:)=f(n,:);
end
p=z %输出地震力
v=zeros(cn,1); %求解水平地震力效应值
for i=1:cn;
v(i)=sum(p(i,:).^2);
end
V=v.^0.5
子程序
function kcju=matrixju(korc,cn)
kcju=zeros(cn,cn);
for i=1:cn-1;
kcju(i,i)=korc(i)+korc(i+1);
kcju(i,i+1)=-korc(i+1);
kcju(i+1,i)=-korc(i+1);
end
kcju(cn,cn)=korc(cn);。