基于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形。有个冷知识——高阶振型的节点位置对缺陷特别敏感，这解释了为啥实际工程中低阶振动更容易观测到。

代码跑通后可以玩点花的：修改截面参数观察频率变化，或者加个末端质量块看模态切换。搞振动分析就像在给结构把脉，这些频率数值就是结构的"健康指纹"。下次如果看到某座桥的振动频率突然偏移，嘿，说不定就是哪里出问题了。