subfunflexibilitymatrix_cantilever.m

来自「主要用于高层建筑的减振设计,根据线性振动理论,用matlab高级语言编写而成,由」· M 代码 · 共 26 行

function FlexM=SUBFunFlexibilityMatrix(L,n,d1,d2,E)
% n--------将杆分为n段
% L--------为杆长

%d1--------钢管的内直径
%d2--------管的外直径
%E---------弹性模量2.1*10^11
Le=L/n;      
r1=d1/2;      %钢管的内半径
r2=d2/2;      %钢管的外半径
I=pi*(d2^4-d1^4)/64; 

FlexM=zeros(n,n);
for j=1:n
     a=Le*j;      %力的作用位置
    for i=1:n
        x=Le*i;   %位移位置
       if (i<j)                   
           FlexM(i,j)=x^2*(3*a-x)/(6*E*I); % (0=<x=<a)
       else
           FlexM(i,j)=a^2*(3*x-a)/(6*E*I); % (a=<x=<L)
       end
    end
end