Bài giảng Động lực học công trình - Nguyễn Trung Kiên

Code Matlab tính tần số riêng của dầm một đầu ngàm một đầu tự do bằng phương pháp phần tử hữu hạn E=1; Module đàn hồi của vật liệu I=1; Moment quán tính của thanh m=1; Khối lượng phân bố của thanh L=1; Chiều dài thanh nele=2; Số lượng phân tử Ll= nele; Chiều dài phần tử thứ nhất L2=L/nele; Chiều dài phần tử thứ hai ml=BeamElementMass(mL1); Ma trận khối lượng phần tử thứ nhất m2=BeamElementMass(m,L2): Ma trận khối lượng phần tử thứ hai kl=BeamElementStiffness(EIL1); Ma trận độ cứng phần tử thứ nhất k2=BeamElementStiffness(E,I,L2); Ma trận độ cứng phần tử thứ hai MO=zeros(6.6); MO=BeamAssemble(MO,m1,1,2); MO=BeamAssemble(MO,m2.2.3); KO=zeros(6,6); KO=Beam Assemble(KO,k1,1,2); KO=Beam Assemble(K0,k2.2.3); M=[M0(316,3:6)]: Điều kiện bên K=[K0(3:6,326)]; Điều kiện biên [V, D] = eigs(K, M, 4, 'SA); omega=sqrt(D) Lắp ghép ma trận độ cứng tổng thể của dầm: function y = BeamAssemble(K,ki.j) K(2 xi- 1,2 xi- 1) = K(2 xi- 1,2 xi-1)+k(1,1); K(2xi-1,2 x i) = K(2 xi- 1,2 x i) + k(1,2);