基于MATLAB的悬臂梁前3阶固有频率和振型求解（假设模态法、解析法、瑞利里兹法）

基于matlab的求解悬臂梁前3阶固有频率和振型 基于matlab的求解悬臂梁前3阶固有频率和振型,采用的方法分别是假设模态法解析法瑞利里兹法 程序已调通可直接运行悬臂梁的振动分析总带着点工程师的浪漫——既要数学的严谨又要面对现实世界的混沌。今天咱们用Matlab折腾三种方法找前三个振动频率和形态手把手看看不同门派的武功秘籍怎么玩。假设模态法最像搭乐高先随便选几个满足边界条件的函数凑数。比如咱们用三次多项式syms x L; phi1 x^2*(6*L^2 - 4*L*x x^2); % 强行满足w(0)0, w(0)0 phi2 x^3*(10*L^3 - 15*L^2*x 6*L*x^2); phi3 x^4*(15*L^4 - 24*L^3*x 10*L^2*x^2);这几个多项式看着复杂其实就是在确保梁固定端位移和转角为零。接下来构建质量矩阵和刚度矩阵时积分运算交给Matlab的int函数自动处理算完特征值问题直接出结果。注意刚度矩阵积分里要包含二阶导数项这步容易手滑写错。解析法玩的是暴力美学直接求解微分方程beta_L [1.875 4.694 7.855]; % 超越方程的前三个根 f (beta_L.^2)./(2*pi*L^2)*sqrt(E*I/(rho*A));这组魔数1.875/4.694/7.855是悬臂梁特征方程的经典解每个搞振动的人手机里都得存着。背后的推导过程需要解四次微分方程考虑到篇幅咱们今天先跳过数学虐恋环节。基于matlab的求解悬臂梁前3阶固有频率和振型 基于matlab的求解悬臂梁前3阶固有频率和振型,采用的方法分别是假设模态法解析法瑞利里兹法 程序已调通可直接运行瑞利里兹法则像在调鸡尾酒选基函数时要带点物理直觉。这里用静力挠曲线作为试探函数psi (x) [x.^2, x.^3, x.^4]; % 静载下的变形趋势 K zeros(3); M zeros(3); for i 1:3 for j 1:3 K(i,j) integral((x) E*I*diff(psi(x),2,i).*diff(psi(x),2,j),0,L); M(i,j) integral((x) rho*A*psi(x,i).*psi(x,j),0,L); end end [V,D] eig(K,M);这个方法的精髓在于基函数选取——既要比解析解简单又要抓住振动形态的关键特征。调试时发现积分步长设置不当会导致矩阵病态这时候把integral换成quadgk能救命。当三种方法输出的频率值相差不到0.5%时那种数值稳定的快感堪比游戏通关。振型图画出来更有意思一阶像跳蹦床二阶是中间有个波节三阶则扭成S形。有个冷知识——高阶振型的节点位置对缺陷特别敏感这解释了为啥实际工程中低阶振动更容易观测到。代码跑通后可以玩点花的修改截面参数观察频率变化或者加个末端质量块看模态切换。搞振动分析就像在给结构把脉这些频率数值就是结构的健康指纹。下次如果看到某座桥的振动频率突然偏移嘿说不定就是哪里出问题了。