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