BÀI TẬP KHÍ ĐÀN HỒI
(08-05-2014)
1 Nguyễn Thanh Phong – G1002398
2 Phạm Tiến Hoàng – G1001131
3 Đinh Minh Tùng – G1003867
4 Hoàng Tiến Đạt – G1000610
Trang 2Từ phương trình (3.3.5):
h
y
M
α
α
− + − −
(1)
Đoạn chương trình dưới đây dùng để giải và rút gọn ra (2) và (3)
clear all
clc
syms K_h x_alpha r2_alpha b d omega_al omega_ha omega
omega_alpha = omega/omega_al
omega_h = omega/omega_ha
m = K_h/omega_h^2
S_alpha = x_alpha*m*b
I_alpha = r2_alpha*m*b^2
matrix = [ m*(-omega^2 + omega_h^2) -S_alpha*omega^2;
-S_alpha*omega^2 I_alpha*(-omega^2 + omega_alpha^2)]
B = [-1; d]
x = matrix^(-1)*B
hpL = simplify(eval(x(1)))
alphapL = simplify(eval(x(2)))
Trang 3 Giải tìm ra phương trình (3.3.6)
2 2
1
h
h
h
h
L
x d
K
α
α
ω ω ω ω ω ω
ω ω ω ω ω ω
ω ω
ω ω
ω ω
ω ω
× ÷ − × + × ÷
= −
× ÷ ÷ − ÷ ÷ + ÷ + ÷ −
− + ×
÷ ÷
= −
× × × × × × ×
×
× ÷ ×
2
(2)
1
x r
α
ω ω ω ω
ω ω ω ω
÷ − ÷ ÷ + ÷ + ÷ −
× ×
Đặt:
.
2
h h
r
α α
α α
ω
Rút gọn ta được:
2
2
1
h
x d
b r h
K
r
α
α
ω ω ω ω
ω ω ω ω
÷
− − ÷ + × × ÷
÷
=
− ÷ ÷÷ − ÷ ÷÷− ÷ ÷
Trang 4 Giải tìm ra phương trình (3.3.8)
2
2 2
1
h
h
L
d
b r
α α
α α
ω ω ω
ω ω ω α
ω ω ω ω ω ω ω
ω ω ω ω ω ω ω
ω
ω
ω
ω
÷
× − + ×
÷ ÷ ÷ ÷
=
− + + −
÷ ÷ ÷ ÷ ÷ ÷ ÷
× ×
÷
= ×
× × × × ×
÷
2
2
1
(3)
1
h
x r x
r
α α α
ω ω
ω ω
ω ω ω ω ω ω
ω ω ω ω ω ω
× × × ×
÷ − +÷ × ÷
÷
× ÷ ÷ − ÷ ÷ + ÷ + ÷ −
Đặt:
.
2
h h
r
α α
α α
ω
Rút gọn ta được:
2
2 2
1
1
1 1
h
x d
r
α
α
α
ω ω ω
ω ω ω
α
ω ω ω ω ω
ω ω ω ω ω
× × ÷ − +÷ ÷
÷ ÷
= ×
÷ ÷
× − − +
÷ ÷ ÷ ÷ ÷ ÷ ÷
×
×