在求解结构的有阻尼固有频率率和振型时,是否需要考虑该结构的阻尼情况

来源：蜘蛛抓取(WebSpider) 时间：2016-04-19 07:25 标签： 固有振型函数
教有阻尼时MATLAB中求固有频率以及振型
当有阻尼时（M&&质量矩阵；C&& 阻尼矩阵；K&&刚度矩阵）
M*X''+CX'+KX=0
有一个近似算法是现在常用的，就是先求出无阻尼系统的固有频率及相应振型矩阵，用命令eig(K,M,'qz')或者eig(K,M)，将系统正则化，求出正则阻尼矩阵，其可能不是对角矩阵，但取其对角元素即可，那么阻尼系统固有频率为wr=w*sqrt(1-c),其中c为正则阻尼矩阵对角线上元素，而系统的振型不变，因此原来地振型矩阵中各列向量仍对应阻尼系统各固有频率
[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
&& [v,w]=eig(m\k);
&& w=sqrt(w);
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;
%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;
w=diag(w);zeta=diag(zeta);
已投稿到：
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。结构动力学问答题答案 武汉理工 研究生_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
结构动力学问答题答案 武汉理工 研究生
上传于||暂无简介
阅读已结束，如果下载本文需要使用1下载券
想免费下载本文？
下载文档到电脑，查找使用更方便
还剩4页未读，继续阅读
你可能喜欢abaqus第九章 动力问题_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
abaqus第九章 动力问题
上传于||文档简介
&&a​b​a​q​u​s​第​九​章​ ​动​力​问​题
阅读已结束，如果下载本文需要使用0下载券
想免费下载更多文档？
下载文档到电脑，查找使用更方便
还剩18页未读，继续阅读
你可能喜欢对机械结构进行模态分析,得到它的固有频率,振型,阻尼等参数 请问：这些参数有什么用?例如：对结构改进有何种指导意义?对后续的控制部分设计又有什么作用?
AOI圣诞二1651
哈!简单说吧,通过对这些数据的了解和分析可知你这机械结构设计上的不合理处,并在后的振动损坏的起因上找到最不合理点好进行补充或加固,对后还要加具有振动源性质的机械时可进行选 择如电机对它的转速选择,最好不要加到同一振动频率上了而引起共振而加速缩短你这机械设备的寿命.
为您推荐：
其他类似问题
扫描下载二维码

我要回帖

更多关于 固有振型函数 的文章

 

随机推荐