subfunflexibilitymatrix_freebeam.m

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

clear;clc;

n=2;          %将杆分为n段
L=60;          %L--为杆和
Le=L/(n+1);      

r1=3.9/2;      %钢管的内半径
d1=2*r1;  
r2=4.0/2;      %钢管的外半径
d2=2*r2; 
P=7900;        %钢管的密度
A=pi*(r2^2-r1^2); 
I=pi*(d2^4-d1^4)/64 ;
E=2.1*10^11;

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