214
function Daodongnho_enliptic
m1 = 50; m2 = 40; % kg
l = 5; % m
g = 9.81; % m/s^2
c = 2000; % N/m
M = [m1 + m2, m2*l; m2*l m2*l^2];
C = [c, 0; 0, m2*g*l];
n = size(M,1);
[V, D] = eig(C,M);
for i=1:n ome(i)=sqrt(D(i,i)); end
ome, v1 = V(:,1), v2 = V(:,2)
x = [0:1:n];
subplot(2,1,1)
plot(x, [0; V(:,1)],'k-','linewidth',1), grid on, ylabel('v_1')
subplot(2,1,2)
plot(x,[0; V(:,2)],'k-','linewidth',1), grid on, ylabel('v_2')
Các tần số riêng và các dạng
dao ñộng riêng
ome =
1.3727
6.4535
v1 =
-0.0062
-0.0303
v2 =
-0.1413
0.0297
0
0.5
1
1.5
2
- 0.04
- 0.02
0
v
1
0
0.5
1
1.5
2
- 0.2
0
0.2
v
2
Hình 9-29. Các dạng dao ñộng riêng
Hệ dao ñộng ba bậc tự do
Cho hệ dao ñộng ba bậc tự do như trên hình 9-30. Các lò xo tuyến tính có ñộ cứng
i
c
, khối lượng các vật nặng là
i
m
, (
1, 2, 3
i
=
). Gọi
i
x
là dịch chuyển của các khối
lượng kể từ vị trí mà các lò xo không bị biến dạng. Áp dụng phương pháp tách vật
(hoặc phương trình Lagrange loại 2) nhận ñược phương trình vi phân chuyển ñộng
ở dạng ma trận như sau:
+
=
Mx
Cx
p
ɺɺ
,
1
2
3
[
]
T
x x x
=
x
