基于Matlab的多*度系统固有频率及振型计算

时间:2021-11-14 09:29:15

可参考文涛,基于Matlab语言的多*度振动系统的固有频率及主振型计算分析,2007

对于无阻尼系统

[VEC,VAL]=eig(inv(A)*K)

对于有阻尼系统,参考振动论坛计算程序

输入M,D,K

function [v,w,zeta]=vbr_sf(m,d,k)

%vbr_sf  vbr_sf(m,d,k)

%        [v,w,zeta]=vbr4(m,d,k)

%        function vbr_sf finds the mode shapes and natural frequencies of

%        a linear second order matrix equation.  有阻尼二阶矩阵方程

%        [v,w]=vbr_sf(m,k) finds the mode shapes and natural frequencies

%        for the undamped case.

if nargin==2

k=d;

[v,w]=eig(m\k);

w=sqrt(w);

end

if nargin==3

if norm(d/m*k-k/m*d) < 1e-8*norm(k/m*d)

%disp('Damping is proportional, eigenvectors are real.')

[v,w]=eig(m\k);

w=sqrt(w);

zeta=(v'*m*v)\(v'*d*v)/2/w;

else

%disp('Damping is non-proportional, eigenvectors are complex.')

a=[0*k eye(length(k));-m\k -m\d];

[v,w1]=eig(a);

w=abs(w1);

zeta=-real(w1)/w;

end

end

w=diag(w);zeta=diag(zeta);

振动系统的特性包括固有特性,固有特性一般指的是没有激励对应数学齐次方程的特征,也就是特征解,包括特征值(物理上常称固有频率)和特征向量(物理上常称振型)。固有特性是振动系统的一种自身固有特征,也可以这么理解,当系统以某个固有频率振动时,振动的振型一定是对应的固有振型。

*振动准确理解是在初始激扰(初始力或初位移)作用下,激扰撤销后系统的振动,对于线性无阻尼系统,由于系统存在正交性,若初始激扰是单频的,那么激励撤销后,系统仍然按激励的频率振动下去,振动响应是所有固有振型的叠加,仅仅当激励的频率等于某个固有频率时,振动的振响应就是对应的振型的若干倍;对于线性有阻尼系统,差别是振动的振响应会逐渐减小... ..

如果激励是多频率成分,问题要复杂一些.

响应是振型的叠加,频率不存在叠加... ..

系统的固有特性与惯性、弹性和耗散等有关,主要是惯性和弹性,惯性实际工程中通常不会有什么变化,而弹性相对比较容易发生变化,材料常数、剖面模数(断面出现裂纹、断裂. . .)、边界刚度等等发生变化都会显现系统变形刚度发生变化,导致系统固有频率随之发生变化。

(转自举聚)

基于Matlab的多*度系统固有频率及振型计算的更多相关文章

  1. 单*度系统中质量、阻尼和刚度变化对频率响应函数&lpar;FRF&rpar;影响图的绘制

    作者:赵兵 日期:2020-02-17 目录 单*度系统中质量.阻尼和刚度变化对频率响应函数(FRF)影响图的绘制 1.     背景 2.     VISIO绘制 3.     Matlab绘制 ...

  2. 基于MATLAB边缘检测算子的实现

    基于MATLAB边缘检测算子的实现 作者:lee神 1.   概述 边缘检测是图像处理和计算机视觉中的基本问题,边缘检测的目的是标识数字图像中亮度变化明显的点.图像属性中的显著变化通常反映了属性的重要 ...

  3. 基于MATLAB的多项式数据拟合方法研究-毕业论文

    摘要:本论文先介绍了多项式数据拟合的相关背景,以及对整个课题做了一个完整的认识.接下来对拟合模型,多项式数学原理进行了详细的讲解,通过对文献的阅读以及自己的知识积累对原理有了一个系统的认识.介绍多项式 ...

  4. 无人机基于Matlab&sol;Simulink的模型开发&lpar;连载一)

    "一切可以被控制的对象,都需要被数学量化" 这是笔者从事多年研发工作得出的道理,无论是车辆控制,机器人控制,飞机控制,还是无人机控制,所有和机械运动相关的控制,如果不能被很好的数学 ...

  5. 基于MATLAB的单级倒立摆仿真

    有关代码及word文档请关注公众号“浮光倾云”,后台回复A010.02即可获取 一.单级倒立摆概述 倒立摆是处于倒置不稳定状态,人为控制使其处于动态平衡的一种摆,是一类典型的快速.多变量.非线性.强耦 ...

  6. RDIFramework&period;NET — 基于&period;NET的快速信息化系统开发框架 — 系列目录

    RDIFramework.NET — 基于.NET的快速信息化系统开发框架 — 系列目录 RDIFramework.NET,基于.NET的快速信息化系统开发.整合框架,给用户和开发者最佳的.Net框架 ...

  7. Kafka 分布式的,基于发布&sol;订阅的消息系统

    Kafka是一种分布式的,基于发布/订阅的消息系统.主要设计目标如下: 通过O(1)的磁盘数据结构提供消息的持久化,这种结构对于即使数以TB的消息存储也能够保持长时间的稳定性能. 高吞吐量:即使是非常 ...

  8. 项目-基于视频压缩的实时监控系统--tiny6410

    项目-基于视频压缩的实时监控系统--tiny6410 @国嵌linux学习笔记. 1. 构造服务端结构体 server struct server { int epfd; //保存epoll指针 st ...

  9. 基于MATLAB的离散小波变换

    申明,本文非笔者原创,原文转载自:  基于Matlab的离散小波变换         http://blog.sina.com.cn/s/blog_725866260100ryh3.html 简介 在 ...

随机推荐

  1. ASP&period;NET MVC入门之再不学习就真的out了

    听说最近又出了什么SAM,MVC辉煌即将过去,惊了我一身冷汗,ASP.NET MVC是啥都还没搞明白呢 于是赶紧打开ASP.NET官网学习学习,欢迎各位高手大侠来指点指点

  2. 手机safari图片上传竖变横处理

    在手机safari上传图片时,竖着的照片会变成横着的照片,以下程序片段利用图片exif信息把图片旋转回去,代码抄自php.net官网. http://php.net/manual/zh/functio ...

  3. Windows服务定时执行任务

    1.创建多线程类 /// <summary> /// 多线程 /// </summary> public abstract class MuliThread<T> ...

  4. 安卓第九天笔记-Activity

    安卓第九天笔记-Activity 1.创建Activity 一个界面对应一个activity 创建一个Activity  1.写一个JAVA类,继承Activity publicclass CalcA ...

  5. SQL Server2008ldf文件太大

    --适用于SQL Server 2008的方法 USE [master] GO ALTER DATABASE RmyyHisDW SET RECOVERY SIMPLE WITH NO_WAIT GO ...

  6. 鸟哥的linux私房菜学习记录之档案权限与目录配置

    在linux中可以通过ls来查看文件 如ls -al,可以看到类似以下的内容 给个例子来理解下 在目录中如果只有r权限没有x权限无法进入该目录

  7. java 导入自定义类

    eclipse导入很容易,昨天上课学了一下用记事本写java,导入自定义类,这就麻烦了. 代码贴一下,方便操作: package tom.jiafei; public class SquareEqua ...

  8. 使用main方法调用http请求本地服务器的某个servlet报错问题

    java.io.IOException: Server returned HTTP response code: 500 for URL: http://localhost:8081/test/myS ...

  9. hadoop 及hbase zookeeper 经常出现问题

    往往是以下几个 1/ 各节点时间不统一(写shell文件统一时间) 2/配置文件 /etc/hosts文件中ip地址配置错误(更新ip) 3/断网后重启机器 ip地址被修改(更新配置文件中的ip)

  10. &lbrack;CSS3&rsqb;学习笔记-文字与字体相关样式

    1.给文字添加阴影 <!doctype html> <html> <head> <meta charset="utf-8"> &lt ...