悬臂梁的振动问题(微分求积法求解悬臂梁自由振动程序)
Posted
篇首语:知识就是飞上天的羽翼。本文由小常识网(cha138.com)小编为大家整理,主要介绍了悬臂梁的振动问题(微分求积法求解悬臂梁自由振动程序)相关的知识,希望对你有一定的参考价值。
悬臂梁的振动问题(微分求积法求解悬臂梁自由振动程序)
微分求积法(简称DQM)是Bellman等人于上世纪70年代初提出的一种新的数值计算方法。这种算法由于不依赖泛函和变分原理,具有数学原理简单、计算精度高,计算量和内存需求小。微分求积法从提出至今,由于其具有计算量小和精度高等优点而不断受到重视,目前对该数值方法本身的研究尽管已经相当成熟。给大家分享一个Matlab实现悬臂梁自由振动的小程序,帮助大家学习微分微分求积法。
clear
digits()
u=i; %保留虚数
N=15; %节点数
X=[]; %离散的横坐标
R=[]; %DQ中的Rij
RR=[]; %用于计算R的中间矩阵
for i=1:N
X(i)=0.5*(1-cos(pi*(i-1)/(N-1)));
end
%建立RR和R
for i=1:N
for j=1:N
if i~=j
RR(i,j)=X(i)-X(j);
else
RR(i,j)=1; %建立完RR
end
end
end
R=prod(RR,2); %建立完R,矩阵按行排列
%建立C1
for i=1:N
for j=1:N
if i~=j
C1(i,j)=R(i)/((X(i)-X(j))*R(j));
else
C1(i,j)=0;
end
end
end
for i=1:N
C1(i,i)=-sum(C1(i,:));
end %建立完C1
% RR
% R
% C1
%建立C2
for i=1:N
for j=1:N
C2(i,j)=C1(i,:)*C1(:,j);
end
end %建立完C2
%建立C3
for i=1:N
for j=1:N
C3(i,j)=C2(i,:)*C1(:,j);
end
end %建立完C3
%建立C4
for i=1:N
for j=1:N
C4(i,j)=C3(i,:)*C1(:,j);
end
end %建立完C4
% %建立C5
% for i=1:N
% for j=1:N
% C5(i,j)=C4(i,:)*C1(:,j);
% end
% end %建立完C5
% %建立C6
% for i=1:N
% for j=1:N
% C6(i,j)=C5(i,:)*C1(:,j);
% end
% end %建立完C6
K=C4(3:(N-2),1:(N-2));
K=[K(:,1:2),C4(3:(N-2),(N-1):N),K(:,3:end)];
K4=C3(N:N,1:(N-2)); %用于形成新K的第4行
K4=[K4(:,1:2),C3(N:N,(N-1):N),K4(:,3:end)];
K3=C2(N:N,1:(N-2)); %用于形成新K的第3行
K3=[K3(:,1:2),C2(N:N,(N-1):N),K3(:,3:end)];
K2=C1(1:1,1:(N-2)); %用于形成新K的第2行
K2=[K2(:,1:2),C1(1:1,(N-1):N),K2(:,3:end)];
K1=zeros(1,N);
K1(1,1)=1; %用于形成K的第1行
K=[K1;K2;K3;K4;K];
Sbb=K(1:4,1:4);
Sbb=inv(Sbb); %考虑求Sbb的逆
Sbd=K(1:4,5:end);
Sdb=K(5:end,1:4);
Sdd=K(5:end,5:end);
S=-Sdb*Sbb*Sbd+Sdd;
lambda=eig(S); %考虑求S的特征值
w=sqrt(sqrt(lambda));
w=sort(w)
相关参考
悬臂梁振动计算(一种新型振动能量收集系统,可持续点亮LED灯和数码管)
沈阳工业大学机械工程学院的研究人员刘慧芳、曹崇东、赵强、马凯、谷艳玲,在2020年第14期《电工技术学报》上撰文,基于悬臂式铁镓合金的动力学模型、等效电路模型和伏安特性建立悬臂式铁镓合金振动能量收集装置的输出...
悬臂梁振动固有频率(可调功能材料(III):电流变材料简介)
功能材料(FunctionalMaterials)是指通过光、电、磁、热、化学、生化等作用后具有特定功能的材料,有时也被称作特种材料(SpecialityMaterials)或精细材料(FineMaterials)。功能材料的功能涉及面广,具体包括光功能、电功能,磁功...
悬臂梁振动参数的实验比较(数显悬臂梁冲击试验机、数显简支梁冲击机等参数,您了解多少?)
...用最为普遍!根据数显冲击试验机的构造,其可分为数显悬臂梁冲击试验机、数显简支梁冲击试验机、数显低温全自动冲击试验机等3种,那么关于三者的参数,您又了解多少呢?数显悬臂梁冲击试
总结石化企业在发展新能源、推动绿色低碳技术创新技术等方面的实践——中国正式向国际社会提出争取使二氧化碳排放量在2030年前达到峰值,在2060年前实现碳中和这一目标。本文首先分析了碳中和背景下石油石化行业发展面...
相信大家对悬挑梁不陌生了,但也有很多人经常把悬臂梁和悬挑梁弄混淆,为了解决这一现象,玖辰建材小编在这里为大家整理了一份关于“悬臂梁和纯悬挑梁的区别是什么”的相关内容,希望可以帮助到您。悬臂梁和纯悬挑梁...
...减小悬臂长度3.勿忽略工作载荷可以产生的有利作用4.受振动载荷的零件避免用摩擦传力5.避免机构中的不平衡力6.避免只考虑单一的传力途径7.不应忽略在工作时零件变形对于受力分布的影响8.避免铸铁件受大的拉伸应力9.避免细...
无弹性元件的挠性联轴器(机械设计500禁忌——非常经典,值得收藏)
...减小悬臂长度3.勿忽略工作载荷可以产生的有利作用4.受振动载荷的零件避免用摩擦传力5.避免机构中的不平衡力6.避免只考虑单一的传力途径7.不应忽略在工作时零件变形对于受力分布的影响8.避免铸铁件受大的拉伸应力9.避免细...
悬挂式保护焊机悬臂架(2022安徽省建筑电焊工(建筑特殊工种)考试题模拟考试题库模拟考试)
题库来源:安全生产模拟考试一点通公众号小程序2022年安徽省建筑电焊工(建筑特殊工种)操作证考试题为安徽省建筑电焊工(建筑特殊工种)操作考试题考前押题密卷!2022安徽省建筑电焊工(建筑特殊工种)考试题模拟考试题库模拟...
悬臂梁(ZEMIC悬臂梁式称重传感器H8C-C3-10t-4B1)
【广州兰瑟电子】杨工提供ZEMIC悬臂梁式称重传感器H8C-C3-1.0t-4B1产品资料。一、ZEMIC品牌H8C-C3-1.0t-4B1称重传感器产品介绍产品品牌:ZEMIC产品系列:H8C产品类型:悬臂梁式称重传感器产品量程:100kg,200kg,500kg,1.0t,2.0t,2.5t,3.0...
塑料悬臂梁摆锤冲击试验机-尼龙抗冲击强度性能测试仪-陶瓷冲击试验仪器-厦门科王电子冲击试验是用来评价材料在高速载荷状态下的韧性或对断裂的抵抗能力的试验。塑料材料的冲击强度在工程应用上是一项重要的性能指标,...