%传递矩阵法求简支梁固有频率的近似解 clc clear syms p q Sp = sym('[1 0 0 0;0 1 0 0;0 0 1 0;x 0 0 1]'); %点传递矩阵 Sf = sym('[1 1 1/2 1/6;0 1 1 1/2;0 0 1 1;0 0 0 1]'); %场传递矩阵 n = input('输入划分单元数:'); S = ((Sf*Sp)^(n-1))*Sf; %两端支座之间的传递矩阵 %求固有频率 xs = solve(S(1,2)*S(3,4)-S(1,4)*S(3,2)); xs = sort(double(xs)); xt = xs*n^4; xt = sqrt(xt); xe(1:n-1) = (pi*(1:n-1)).^2; %精确解 xe = xe'; fprintf('传递矩阵法的结果:\n') for i = 1:n-1 fprintf('第%d阶固有频率:%8.4f(EI/ml^3)^(1/2)\n',i,xt(i)) end %求模态 step = 1/n; for i = 1:n-1 f0 = -S(3,2)/S(3,4); f0 = subs(f0,'x',xs(i)); xk(:,1) = [0 1 0 f0]'; for j = 2:n+1 xk(:,j) = Sf*Sp*xk(:,j-1); xk(:,j) = subs(xk(:,j),'x',xs(i)); end xkk = xk(1,2:n);
传递矩阵法求简支梁固有频率的近似解 --matlab程序
最新推荐文章于 2021-03-21 10:35:29 发布