subfunkmatrix_cantilever.m

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

function KM=SUBFunKMatrix(L,n,d1,d2,E)
% n--------将杆分为n段
% L--------为杆长
% P--------钢管的密度
%d1--------钢管的内直径
%d2--------管的外直径
%E---------弹性模量2.1*10^11
Le=L/n;      
r1=d1/2;      %钢管的内半径
r2=d2/2;      %钢管的外半径
I=pi*(d2^4-d1^4)/64; 
k=3*E*I/Le^3;

KM=zeros(n,n);
KM(1,1)=2*k;
for i=2:n
   KM(i,i)=2*k;
   KM(i,i-1)=-k;
   KM(i-1,i)=-k;
end
KM(n,n)=k;