BÀI TẬP KHÍ ĐÀN HỒI
(08-05-2014)
1. Nguyễn Thanh Phong
2. Phạm Tiến Hoàng
3. Đinh Minh Tùng
4. Hoàng Tiến Đạt
–
–
–
–
G1002398
G1001131
G1003867
G1000610
1
Từ phương trình (3.3.5):
m ( −ω 2 + ωh2 )
− Sα ω 2
h −L
=
Iα ( −ω 2 + ωα2 ) α M y
− Sα ω 2
(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)))
2
Giải tìm ra phương trình (3.3.6)
2
h
=−
L
2
2
2
2
ω 2 ω 2
ω
2
2 ω
2 ω
2 ω
2
K h × b × ÷ ×
÷ × xα − rα × ÷ ×
÷ + rα × ÷ + rα ×
÷ − rα
ωh ωα
ωh ωα
ωh
ωα
2
=−
2
ω
ω
2
2
b×
÷ × rα − b × rα + d ×
÷ × xα
ωα
ωα
2
ω
d ω xα
÷ −1 + ×
÷×
b ωα rα2
ωα
ω 2 ω 2 x 2 ω 2 ω 2 ω 2 ω 2
α
K h × ÷ ×
÷ × 2 − ÷ ×
÷ + ÷ +
÷ −1
ω
ω
r
h α
α
ωh ωα ωh ωα
(2)
Đặt:
d≡
My
L
K
ωh2 ≡ h
m
Sα xc. g
=
mb
b
I
rα2 = α 2
mb
xα ≡
Rút gọn ta được:
ω 2 d x ω 2
− 1 −
+ × α2 ×
÷÷
÷
÷
ω
α b rα ωα
h
=
L
ω 2 ω 2 ω 2 ω 2 x 2
α
K h × 1 −
− ÷ ×
÷ ÷1 − ÷ ÷
÷× 2
÷
ω
ω
ω
r
ωα ÷
α
h h α
3
Giải tìm ra phương trình (3.3.8)
α
=
L
2
2
2
ω ω
ω
×
d
×
−
d
+
b
×
×
x
÷
÷
÷ α÷
÷
ωα ωh
ωh
2
2
2
2
2
2
2
ω ω ω
ω
2
2 ω
2 ω
2 ω
2
K h × b × ÷ × ÷ ×
÷ × xα − rα × ÷ ×
÷ + rα × ÷ + rα ×
÷ − rα
ωh ωh ωα
ωh ωα
ωh
ωα
2
2
2
x ω 2
d 1 ω
ω
× 2 × ÷ −1÷+ α2 × ÷
÷
÷ rα ωh
b rα ωh
ωα
=
×
2
2
2
ω ω x 2 ω 2 ω 2 ω 2 ω 2
ω
α
÷ K h × b × ÷ ×
÷ × 2 − ÷ ×
÷ + ÷ +
÷ −1
ω
ω
ω
r
h
h α
α
ωh ωα ωh ωα
(3)
Đặt:
d≡
My
L
K
ωh2 ≡ h
m
Sα xc. g
=
mb
b
I
rα2 = α 2
mb
xα ≡
Rút gọn ta được:
ω
α ωα
=
L ω
ωh
2
2
ω 2 x
d 1 ω
× × ÷ −1÷+ ÷ × α2
÷
÷ ω
b rα2 ωh
rα
×
h
2
ω 2 ω 2 ω 2 ω 2 x 2
α
÷+ ÷ ×
÷ K h × b × ÷ − 1÷ 1 −
÷ × 2
÷ ωα ÷
÷
ω
ω
r
ωh
α
h α
4