Slide động lực học công trình_Chương 3

CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3.1 Phương trình vi phân tổng quát: Xét dao động thẳng có khối lượng phân bố m(z) dọc theo ciều dài Hệ có bậc tự vô Khi chịu lực kích thích thay đổi theo thời gian có phương nghiêng so với trục Dao động ngang xác định phương trình y = y(z, t) hàm tọa độ z tiết diện ngang thời gian t biểu thị đường đàn hồi Từ liên hệ vi phân đường hồi y(z, t), mô men uốn M(z, t) cường độ tải trọng phân bố p(z, t): ∂ y( z, t ) ∂ M ( z, t ) = − M ( z, t ); = p( z, t ) EI( z ) 2 ∂z ∂z CHƯƠNG 3: ∂2 ∂z DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO  ∂ y( z, t )   EI( z )  = − p( z , t ) ∂z   q(z,t) p(z,t) > có chiều hướng lên * Tải trọng kích thích bố với cường độ q(z,t) tác dụng vuông góc với trục >0 có chiều hướng lên * Lực quán tính khối lượng phân bố m(z) hướng theo chiều chuyển động bằng: ∂ y( z, t ) − m( z ) ∂t ∂ y ( z, t ) − m( z ) ∂t y z r(z,t) * Lực cản r(z,t) ngược chiều với chiều chuyển động CHƯƠNG 3: ⇒ DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO ∂ y( z, t ) + r ( z, t ) p( z , t ) = q ( z , t ) + m( z ) ∂t Thay p(z, t) vào phương trình đầu tiên, ta thu phương trình vi phân tổng quát mô tả dao động ngang thanh: ∂2 ∂z  ∂ y ( z, t )  ∂ y ( z, t ) + r ( z , t ) = −q( z , t )  EI( z )  + m( z ) 2 ∂z ∂t   Nếu có khối lượng m phân bố đều: ∂ y ( z, t ) ∂ y ( z, t ) + m( z ) + r ( z , t ) = −q( z , t ) EI( z ) ∂z ∂t CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3.2 Dao động riêng không lực cản: 3.2.1 Trường hợp tổng quát: Trong trường hợp r(z, t) = 0, q(z, t) = 0, phương trình vi phân dao động riêng có dạng:  ∂ y ( z, t )  ∂ y ( z, t ) =0  EI( z )  + m( z ) 2 ∂z ∂t   Giải phương trình theo phương pháp tách biến số Fourier ta đặt nghiệm dạng chuỗi tổng nghiệm riêng: ∂2 ∂z ∞ y ( z , t ) = ∑ yi ( z ) Fi (t ) i −1 Lấy đạo hàm thay vào phương trình trên: ∞ ∞ ∂2  (t ) = ′′( z ).F (t )] + m( z )∑ y ( z ) F [ EI ( z ) y ∑ i i i i ∂z i =1 i =1 Cho số hạng tổng phương trình không, với số hạng thứ i, ta thu được: ∂2  (t ) = ′′( z ).F (t )] + m( z ) y ( z ) F [ EI ( z ) y i i i i ∂z ⇒ ∂ ′′  [ EI( z ) yi ( z )] Fi (t ) ∂z =− m( z ) yi ( z ) Fi (t ) Vế phải phụ thuộc vào thời gian t, vế trái phụ thuộc vào z,  tỷ số = const Ký hiệu dại lượng ω i2  có phương trình vi phân với biến số độc lâp: 1) F (t ) + ω F (t ) = i i i Dạng giống phương trình vi phân dao động hệ bậc tự do, nghiệm phương trình là: Fi (t ) = Ai sin ωi t + B cos ωi t = sin(ωi t + ϕi )  Tương ứng với nghiệm riêng yi(z, t)=yi(z).Fi(t), dao động riêng thay đổi điều hòa với tần số riêng ω i 2) ∂2 ′ ′ − ω = [ EI ( z ) y ( z )] m ( z ) y ( z ) i i i ∂z Giải phương trình ta tìm hàm yi(z) biểu thị dạng thứ i dao động riêng ứng với tần số ω i CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3.2.1 Trường hợp EI = const : Phương trình vi phân dao động có dạng: m (z) ∂ y ( z, t ) ∂ y ( z, t ) + =0 ∂z EI ∂t Nghiệm phương trình tìm dạng chuỗi: ∞ y ( z , t ) = ∑ yi ( z ) Fi (t ) i −1 Fi (t ) = Ai sin ωi t + B cos ωi t = sin(ωi t + ϕi ) ∂2 ′ ′ [ EI( z ) yi ( z )] − m( z ).ωi yi ( z ) = ∂z ω m i ki = EI − y ( z ) ki yi ( z ) = 0; IV Giải phương trình đặc trưng: r4 – ki4 = phương trình ta thu nghiệm: r1 = ki ; r2 = −ki ; r3 = iki ; r4 = −iki ; i = − ⇒ Vì: − ki z ki z = + yi ( z ) a.e b.e + c cos ki z + d sin ki z e ki z = chki z + shki z, e − ki z = chki z − shki z  Ta thu phương trình sau: yi ( z ) = C1chki z + C shki z + C3 cos ki z + C sin ki z; yi′( z ) = C1ki shki z + C ki chki z − C3 ki sin ki z + C4 ki cos ki z; Mi ( z) C1ki2 chki z + C ki2 shki z − C3 ki2 cos ki z − C ki2 sin ki z; EI Qi ( z ) ′ ′ ′ = − yi ( z ) C1ki3 shki z + C ki3chki z + C3 ki3 sin ki z − C ki3 cos ki z EI yi′′( z ) = − Dạng yi(z) xem đường đàn hồi nên ta xác định số tích phân Ci theo điều kiện ban đầu Giả sử z = tương ứng với dạng thứ i dao động, ta có thông số ban đầu: độ võng yi(0), góc xoay y’i(0); mô men uốn Mi(0); lực cắt Qi(0) Thay giá trị vào phương trình ta thu được: yi (0) = (C1 + C3 ); yi (0) = (C2 + C4 )ki ; = − − = − − M i (0) EI(C1 C3 )ki ; Qi (0) EI(C2 C4 )ki  M i (0) 1 yi′(0) Qi (0) − ]; C1 = [ yi (0) − ]; C2 = [ ki EI ki ki EI M i (0) 1 yi′(0) Qi (0) + ] C3 = [ yi (0) + ]; C4 = [ ki EI ki K i EI Thay giá trị Ci vào phương trình đầu hệ gồm phương trình , ta có phương trình xác định chuyển vị tương ứng với dạng thứ i dao động riêng viết theo thông số ban đầu: iπ yi ( z ) = Ci sin z; l yi′(0 ) Ci = ki Dạng dao động riêng dầm đơn giản có hai đầu khớp dầm điều hòa theo quy luật hàm số sin với số nửa bước sóng số tần số dao động riêng tương ứng m = const EI = const l π i = 1, y1 = C1 sin z l 3π i = 3, y3 = C3 sin z l i = 2, y2 = C2 sin 2π z l Ví dụ 2: Xác định tần số dao động riêng dầm côngxôn mang khối lượng phân bố m có độ cứng không đổi EI hình vẽ z l y Giải: Tại z =  yi(0) = 0, y’i(0) = 0, Mi(0) = ?, Qi(0) = ? M i (0) Qi (0) yi ( z ) = − A3 (ki z ) − A4 (ki z )' ki EI ki EI Qi (0 ) A2 (ki z ), ki Qi ( z ) = M i (0)ki A4 (ki z ) + Qi (0) A1 (ki z ) M i ( z ) = M i (0 ) A1 (ki z ) + Vì z = l Mi(l) = Qi(l) = nên: z l y Qi (0 ) M i ( l ) = M i (0 ) A1 (ki l ) + A2 (ki l ) = ki Qi ( l ) = M i (0)ki A4 (ki l ) + Qi (0) A1 (ki l ) = Để ẩn số khác nghĩa để dao động tồn định thức hệ số hệ phương trình phải không: A1 (k i l ) ki A4 (ki l ) A2 (ki l ) =0 ki A1 (ki l ) ⇒ A12 (k i l ) − A2 (ki l ) A4 (ki l ) = Thay hàm Krưlôv vào ta thu phương trình siêu việt để xác định tần số: D = chki l cos ki l + = Để giải phương trình ta vận dụng cách thử dần Cho kil nhiều giá trị khác tính giá trị D tương ứng: kil D kil 0,2π = 0,628 1,97 1,2π = 3,770 - 16,56 0,4π = 1,257 1,59 1,4π = 4,399 - 11,63 0,6π = 1,885 - 0,04 1,6π = 5,027 24,67 0,8π = 2,514 - 4,39 π = 3,14 D - 10,57 Nghiệm k1l ≈ 0,6π = 1,885 (chính xác 1,8751) thỏa mãn k2l = 1,49π = 4,68 (chính xác 4,691) Thực tương tự với giá trị kil lớn ta thu được: k3l = 7,855; k4l = 10,996 Các tần số dao động riêng: 1,875 ω1 = l EI 4,69412 ; ω2 = m l EI ; m ,855 ω3 = l2 EI 10,996 ; ω4 = m l2 EI m Tương tự phần trước ta thu dạng dao động riêng hình vẽ: z l y y1 (ω 1) y ( ω 3) y2 (ω 2) 0,5001l 0,7739l 0,8679l CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3.2 Dao động cưỡng hệ chịu lực kích thích tuần hoàn q(z)sinθ t: Khi không kể đến lực cản, phương trình viết: ∂ y( z, t ) m ∂ y ( z, t ) q ( z ) sin θ t + =− ∂z EI ∂t EI Do có lực cản nên sau thời gian dao động riêng dao động cưỡng lực kích thích gây Nghiệm riêng có dạng: y ( z , t ) = y ( z ) sinθ t Thay vào phương trình ta thu được: θ q( z ) m IV 4 y ( z) − k y( z) = − , k = EI EI k- hệ số đặc trưng dao động cưỡng Nghiệm yo(z) phương trình vi phân có dạng sau: yo ( z ) = C1 A1 ( kz) + C A2 ( kz) + C3 A3 ( kz) + C4 A4 ( kz) Aj(kz) với j = 1, 2, 3, hàm Krưlôv Nghiệm riêng yr(z) phụ thuộc vào tải trọng q(z) Xét trường hợp q(z) = q: q y ( z ) = C1 A1 ( kz) + C A2 ( kz) + C3 A3 (kz) + C4 A4 (kz) + k EI Tương tự phần trước, lấy đạo hàm sử dụng điều kiện ban đầu đầu thanh, biến đổi ta thu phương trình biên độ chuyển vị, góc xoay, mô men uốn, lực cắt dao động: yo′ Mo Qo = + − − y ( z ) yo A1 (kz) A2 (kz) A3 (kz) A4 (kz) − k k EI k EI q − [ A1 (kz) − 1] ; k EI M Q q y′( z ) = yo kA4 (kz) + yo′ A1 (kz) − o A2 (kz) − o A3 (kz) − A4 (kz); kEI k EI k EI Q q M ( z ) = − EIyo k A3 (kz) + EIyo′ kA4 (kz) + M o A1 (kz) + o A2 (kz) + A3 (kz); k k q ′ = − + + + + Q( z ) EIyo k A2 (kz) EIyo k A3 (kz) M o kA4 (kz) Qo A1 (kz) A2 (kz) k CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3.4 Dao động cưỡng hệ chịu tải trọng tập trung : yi(a) 3.4.1 Dao động riêng: Xét thẳng có khối lượng phân bố đều, tiết diện không đổi mang khối lượng tập trung mj đặt hoành độ a hình vẽ a Zj(t) mj EI m l-a l Khi dao động với tần số ω i, đường đàn hồi xác định theo phương trình yi(z) dạng thứ i Tại khối lượng mj phát sinh lực quán tính: Z j = m jω yi (a ) i Với : ω m i = ki EI ⇒ ωi2 = i k EI m EI ⇒ Zj = k m j yi ( a ) m i Lực quán tính Zj(t) xem điều kiện gián đoạn tải trọng z = a dùng phương pháp thông số ban đầu xét phần trước Sau lập phương trình cho đoạn thanh, sử dụng điều kiện biên để thiết lập phương trình xác định thông số k Từ suy tần số dao động riêng Với ≤ z ≤ a, phương trình chuyển vi, mô men uốn động đoạn I: yi′(0) Qi (0) y ( z) = A2 (ki z ) − A4 (ki z ); ki ki EI I i Qi (0) ′ = − + M ( z) ki EIyi (0) A4 (ki z ) A2 (ki z ) ki I i Biên độ lực quán tính khối lượng mj: EI yi′(0) Qi (0) Zj = k mj [ A2 (ki a) − A4 (ki a)] m ki ki EI i Xét đoạn II với a ≤ z ≤ l , cho z1 = z – a : Zj II I yi ( z ) = yi ( z ) + A4 (ki z1 ) ki EI mj yi′(0) = [ A2 (ki z ) + ki A2 (ki a) A4 (ki z1 )] − ki m mj Qi (0) + [ A4 (ki z ) ki A4 (ki a) A4 (ki z1 )]; ki EI m Zj II I Mi ( z) = Mi ( z) − A2 ( ki z1 ) ki EI − = − EIki yi′(0 )[ A4 ( ki z ) + ki mj m A2 ( ki a ) A2 ( ki z1 )] + mj Qi (0 ) + [ A2 ( ki z ) + ki A4 ( ki a ) A2 ( ki z1 )] ki m Khi z = l ta có điều kiện biên bên phải: y(l) = 0, M(l) = 0, vậy: mj yi′ (0 ) y (l ) = [ A2 ( ki l ) + ki A2 ( ki a ) A4 ( ki b)] − ki m II i mj Qi (0 ) − [ A4 ( ki l ) + ki A4 ( ki a ) A4 ( ki b)] = 0; ki EI m M ( z ) = − EIki yi′ (0 )[ A4 ( ki l ) + ki II i mj m A2 (ki a ) A2 ( ki b) + mj Qi (0 ) + [ A2 (ki l ) + ki A4 ( ki a ) A2 ( ki b)] = ki m Với b = l – a; Tương tự phần trước ta xác định thông số ki, sau tìm ddwwocj tần số dao động riêng tương ứng CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3.4 Dao động cưỡng hệ chịu tải trọng tập trung : 3.4.2 Dao cưỡng chịu lực kích động tuần hoàn: Tương tự phần trước, khối lượng mj phát sinh lực quán tính: Zj(t) = Zjsinθ t với biên độ Zj = mjθ 2y(a) Biên độ lực quán tính Zj tải trọng Po xem điều kiện gián đoạn tải trọng z = a vị trí đặt lực P0sinθ t áp dụng phương trình biết cho đoạn Quá trình tính toán đến hoàn tất hoàn toàn tương tự [...]... dao động riêng như hình vẽ: z l y y1 (ω 1) y 3 ( ω 3) y2 (ω 2) 0,5001l 0,7 739 l 0,8679l CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3. 2 Dao động cưỡng bức khi hệ chịu lực kích thích tuần hoàn q(z)sinθ t: Khi không kể đến lực cản, phương trình được viết: ∂ 4 y( z, t ) m ∂ 2 y ( z, t ) q ( z ) sin θ t + =− 2 ∂z 4 EI ∂t EI Do có lực cản nên sau một thời gian dao động riêng sẽ mất đi và chỉ còn dao động. .. yo′ A1 (kz) − o A2 (kz) − 2 o A3 (kz) − 3 A4 (kz); kEI k EI k EI Q q M ( z ) = − EIyo k 2 A3 (kz) + EIyo′ kA4 (kz) + M o A1 (kz) + o A2 (kz) + 2 A3 (kz); k k q 3 2 ′ = − + + + + Q( z ) EIyo k A2 (kz) EIyo k A3 (kz) M o kA4 (kz) Qo A1 (kz) A2 (kz) k CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3. 4 Dao động cưỡng bức khi hệ chịu tải trọng tập trung : yi(a) 3. 4.1 Dao động riêng: Xét thanh thẳng có khối lượng... định các thông số ki, sau đó tìm ddwwocj tần số dao động riêng tương ứng CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3. 4 Dao động cưỡng bức khi hệ chịu tải trọng tập trung : 3. 4.2 Dao cưỡng bức chịu lực kích động tuần hoàn: Tương tự như phần trước, tại khối lượng mj phát sinh lực quán tính: Zj(t) = Zjsinθ t với biên độ Zj = mjθ 2y(a) Biên độ của lực quán tính Zj và của tải trọng Po được xem điều kiện... = 3, 770 - 16,56 0,4π = 1,257 1,59 1,4π = 4 ,39 9 - 11, 63 0,6π = 1,885 - 0,04 1,6π = 5,027 24,67 0,8π = 2,514 - 4 ,39 π = 3, 14 D - 10,57 Nghiệm k1l ≈ 0,6π = 1,885 (chính xác 1,8751) thỏa mãn k2l = 1,49π = 4,68 (chính xác 4,691) Thực hiện tương tự với những giá trị kil lớn hơn ta thu được: k3l = 7,855; k4l = 10,996 Các tần số dao động riêng: 1,875 2 ω1 = 2 l EI 4,69412 ; ω2 = 2 m l EI ; m 7 ,855 2 3 =... tương ứng m = const EI = const l π i = 1, y1 = C1 sin z l 3 i = 3, y3 = C3 sin z l i = 2, y2 = C2 sin 2π z l Ví dụ 2: Xác định tần số dao động riêng của dầm côngxôn mang khối lượng phân bố đều m và có độ cứng không đổi EI như hình vẽ z l y Giải: Tại z = 0  yi(0) = 0, y’i(0) = 0, Mi(0) = ?, Qi(0) = ? M i (0) Qi (0) yi ( z ) = − 2 A3 (ki z ) − 3 A4 (ki z )' ki EI ki EI Qi (0 ) A2 (ki z ), ki Qi ( z )... C3 A3 (kz) + C4 A4 (kz) + 4 k EI Tương tự như phần trước, lần lượt lấy đạo hàm và sử dụng các điều kiện ban đầu ở đầu thanh, biến đổi ta thu được các phương trình biên độ chuyển vị, góc xoay, mô men uốn, lực cắt khi thanh dao động: yo′ Mo Qo = + − − y ( z ) yo A1 (kz) A2 (kz) A3 (kz) A4 (kz) − 2 3 k k EI k EI q − 4 [ A1 (kz) − 1] ; k EI M Q q y′( z ) = yo kA4 (kz) + yo′ A1 (kz) − o A2 (kz) − 2 o A3... tính chất sau: * A1(0) = 1; A2(0) = 0; A3(0) = 0; A4(0) = 0 * Giữa các hàm có sự liên hệ vi phân tuân theo quy tắc vòng tròn như hình vẽ: A1′( ki z ) = ki A4 (ki z ) A4′ (ki z ) = ki A3 ( ki z ) A3′ (ki z ) = ki A2 ( ki z ) A2′ (ki z ) = ki A1 (ki z ) A2 A1 A3 A4 Từ phương trình: yi′(0) M (0 ) Q (0) yi ( z ) = yi (0) A1 (ki z ) + A2 (ki z ) − 2i A3 (ki z ) − 3i A4 (ki z ), ki ki EI ki EI  Các phươngtrình... dao động riêng Với 0 ≤ z ≤ a, các phương trình chuyển vi, mô men uốn động trong đoạn I: yi′(0) Qi (0) y ( z) = A2 (ki z ) − 3 A4 (ki z ); ki ki EI I i Qi (0) ′ = − + M ( z) ki EIyi (0) A4 (ki z ) A2 (ki z ) ki I i Biên độ lực quán tính tại khối lượng mj: EI yi′(0) Qi (0) Zj = k mj [ A2 (ki a) − 3 A4 (ki a)] m ki ki EI 4 i Xét đoạn II với a ≤ z ≤ l , cho z1 = z – a : Zj II I yi ( z ) = yi ( z ) + 3 A4... ( ki z ) A2 ( ki z ) A3 ( ki z ); 2 ki EI ki EI M i ( z ) = − EIyi′′( z ) = − EIyi (0) ki2 A3 ( ki z ) − EIyi′ (0) ki A4 ( ki z ) + + M i (0) A1 ( ki z ) + Qi (0) A2 ( ki z ); ki Qi ( z ) = − EIyi′′′( z ) = − EIyi (0) ki3 A2 ( ki z ) − EIyi′ (0) ki2 A3 ( ki z 0 + + M i (0) ki A4 ( ki z ) + Qi (0) A1 ( ki z ) Tần số dao động riêng: ωi = k 2 i EI m Ví dụ 1: Xác định tần số dao động riêng và lập phương... ki l = 0 Do kil ≠ 0 nên shkil ≠ 0  sinkil = 0  kil = iπ  ki = iπ /l ωi = ki2 2 π EI i = 2 m l 2 * i = 1 ⇒ ω1 = 3, 1416 l2 EI m 2 9,4248 = ⇒ ω = *i 3 3 l2 m = const EI = const l EI 6,2 832 = ⇒ ω = , i 2 2 m l2 2 2 EI 12,5664 = ⇒ ω = , i 4 4 m l2 ' i y (0) Qi (0) yi ( l ) = A2 (ki l ) − 3 A4 (ki l ) = 0 , ki ki EI A2 ( ki l ) 2 ′ ⇒ Qi (0 ) = yi (0 )ki EI A4 ( ki l ) EI , m 2 EI m Do kil = 0 nên: A2(kil) ... y1 (ω 1) y ( ω 3) y2 (ω 2) 0,5001l 0,7 739 l 0,8679l CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3. 2 Dao động cưỡng hệ chịu lực kích thích tuần hoàn q(z)sinθ t: Khi không kể đến lực cản, phương... tìm ddwwocj tần số dao động riêng tương ứng CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3. 4 Dao động cưỡng hệ chịu tải trọng tập trung : 3. 4.2 Dao cưỡng chịu lực kích động tuần hoàn: Tương tự... ∂z ∂t CHƯƠNG 3: DAO ĐỘNG CỦA HỆ CÓ VÔ SỐ BẬC TỰ DO 3. 2 Dao động riêng không lực cản: 3. 2.1 Trường hợp tổng quát: Trong trường hợp r(z, t) = 0, q(z, t) = 0, phương trình vi phân dao động riêng có