›
›

Bài giảng Động lực học công trình - Chương 4: Các phương pháp tính gần đúng

pdf 79 trang hapham 1430
Download
Bạn đang xem 20 trang mẫu của tài liệu "Bài giảng Động lực học công trình - Chương 4: Các phương pháp tính gần đúng", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

  • pdfbai_giang_dong_luc_hoc_cong_trinh_chuong_4_cac_phuong_phap_t.pdf

Nội dung text: Bài giảng Động lực học công trình - Chương 4: Các phương pháp tính gần đúng

  1. CHƯƠNG 4: CÁC PHƯƠNG PHÁP TÍNH GẦN ĐÚNG 4.1. Phương pháp năng lượng Rayleigh: Tổng thế năng U và động năng K ở mọi thời điểm là một hằng số (bỏ qua các tổn thất về năng lượng): U + K = const Khi hệ dao động điều hòa, tại thời điểm ở xa vị trí cân bằng nhất thế năng đạt Umax, động năng K = 0, tại vị trí cân bằng động năng đạt Kmax, thế năng U = 0: Umax + 0 = 0 + Kmax.
  2. m(z)dz Xét hệ mang khối lượng m(z) phân bố m(z) và các khối lượng tập trung m1, m2, m3, như hình m1 mj mn vẽ, dao động theo phương z dz trình: z1 yi (z,t) = yi (z)sin(it + i ). Biểu thức thế năng của hệ khi chỉ kể đến ảnh hưởng của biến dạng uốn có dạng: M 2 (z) U =  dz 2EI  2 y(z,t) M(z) Do = - , z 2 EI
  3. 2 EI y(z,t) 2 EI 2 U =  [- ] d z =  [yi (z)sin(it + i )] dz. 2 z 2 2 2 1 U max =  EI[ yi (z)] dz 2 Động năng của hệ: 2 2 m(z)v mjv j K =  z dz +  . 2 2 Từ phương trình dao động ta có: y(z,t) v = = y (z) cos( t + ) z t i i i i max Vz = yi(z)i ,
  4. y(z ,t) v = j = y (z ) cos( t + ) j t i j i i i max Vj = yi(zj )i ,  2 = i 2 + 2 Kmax [ m(z)yi (z)dz  mj yi (z j )] 2 j Thay vào phương trình cơ bản của Rayleigh: 2  EI[ yi (z)] dz  2 = i 2 + 2  m(z)yi (z)dz mj yi (z j ) j Để xác định i ta cần biết trước dạng dao động chính thứ i của hệ.
  5. Nếu yi(z) là đường đàn hồi do trọng lượng các khối lượng đặt trên hệ gây ra thì thế năng Umax của hệ được xác định bằng công của ngoại lực Tmax: m(z)y (z) m y (z ) = = i + j i j U max Tmax  g dz  g . 2 j 2 +  g.m(z).yi (z)dz  g.mj yi (z j )  2 = j . i 2 + 2  m(z).yi (z)dz  mj yi (z j ) j
  6. Ví dụ1: Tìm tần số dao động riêng của dầm đơn giản có EI m1 m nhịp l, mang khối lượng phân bố đều m và khối lượng tập l/2 l /2 trung m1 = 18ml/35 đặt ở giữa nhịp. Giải: Chọn dao động của dầm là đường đàn hồi do lực P đặt ở giữa nhịp gây ra: Pl 3 3l 2 z - 4z 3 3l 2 z - 4z 3 y(z) = = f , 48EI l 3 l 3 Với giá trị: 0 z l/2, f độ võng giữa dầm.
  7. 24z y (z) = f (- ) l 3 Áp dụng công thức vừa thành lập: l /2 24 2 EI[ f (- )]2 dz l 3 48EI  2 = 0 = l / 2 3l 2 z - 4z 3 18 ml4 2 m[ f ]dz + mlf 2 3 0 l 35 6,9282 EI  = 1/ s l 2 m
  8. Nếu sử dụng công thức thứ 2, đường đàn hồi của dầm vẫn chọn như trên, song f phải là độ võng của dầm ở giữa nhịp do trọng lượng dầm và trọng lượng khối lượng tập trung m1 gây ra: Pl 3 5ql 4 319mgl4 f = + = . 48EI 384EI 13440EI Sau khi thay vào biểu thức ta cũng tìm được: 6,9282 EI  = 1/ s l 2 m
  9. * Nếu dầm không mang khối lượng tập trung m1 = 0. Chọn dạng dao động: i z y (z) = f sin i l l i2 2 i z EI[- f sin ]2 dz 2 4 4 2 0 l l i EI i = = l i z ml4 m[ f sin ]2 dz 0 l 2 EI 9,8696 EI i = 1  = = 1/ s, 1 l 2 m l 2 m 2 2 EI 39,4786 EI i = 2  = = 1/ s, 2 l 2 m l 2 m
  10. Nếu tính theo công thức thứ hai, đường đàn hồi vẫn chọn như trên và f là độ võng giữa dầm: 4 5 l g.2.(1- cosi ) f = mg  = 384 EI i f .i. Với i = 1 ta có: 4g 9,8886 EI  = = 1 . f l 2 m Sai số là 0,19%
  11. h(z) ho Ví dụ 2: Xác định tần số dao động riêng của dầm côngxôn có tiết diện thay đổi z như hình vẽ. Cho biết bề rộng tiết diện ngang b không l đổi, khối lượng, chiều cao tiết diện ngang, mô men quán tính tiết diện ngang thay đổi theo quy luật: z m(z) = m , o l h z3 h(z) = o z, I(z) = I . l o l 3
  12. h(z) ho * Chọn dạng dao động y(z) là đường đàn hồi của dầm có tiết diện không đổi do tải z trọng phân bố đều gây ra: l 4z z 4 y(z) = A(3 - + ), q l l 4 ql 4 A = 24EI l 3 2 z 2 12z 2 EIo 3 A [ 4 ] dz 2 0 l l 39,3762EIo i = = . l z 4z z4 ml 4 2 - + 2 mo A [3 4 ] dz 0 l l l
  13. h 6,275 EI h(z) o  = o 1 2 1/ s l mo z Giá trị đúng: l q 5,315 EI  = o 1 2 1/ s l mo Sai số trong trường hợp này 18%.
  14. h(z) ho Nếu chọn dạng dao động là đường đàn hồi của dầm côngxôn do tải trọng z phân bố bậc nhất gây ra cho l dầm . y(z) = B(1-z/l)2, với 4 q B=ql /12EIo. l 3 z 2 2 EIo 3 ( 2 ) dz 2 0 l l 30EIo 1 = = l z z m l 4 - 2 2 o mo [(1 ) ] dz 0 l l 5,4772 EI  = o Sai số 3% 1 2 1/ s l mo
  15. CHƯƠNG 4: CÁC PHƯƠNG PHÁP TÍNH GẦN ĐÚNG 4.2. Phương pháp năng lượng Lagrange - Ritz: Khi hệ ở trạng thái cân bằng thì thế năng toàn phần U của hệ đạt cực tiểu - là tổng của thế năng của nội lực U* và thế năng ngoại lực T. (chiều ngoại lực hướng xuống là dương, tjees năng của nội lực luôn luôn ngược dấu với thế năng của ngoại lực), ta có: l EI(z) l = * - = 2 - - U U T [y (z)] dz q(z) y(z)dz  Pj y(z j ), 0 2 0 ( j) Pj, q(z) – lực kích thích tập trung và lực kích thích phân bố bao gồm cả các lực quán tính do các khối lượng tập trung và khối lượng phân bố gây ra khi hệ dao động.
  16. CHƯƠNG 4: CÁC PHƯƠNG PHÁP TÍNH GẦN ĐÚNG 4.2.1. Dao động riêng: Xét hệ có khối lượng phân bố m(z) và dao động theo dạng chính yj(z) với tần số riêng j. Phương trình dao động: yj (z,t) = yj (z)sin(jt + j ). Lực quán tính của khối lượng tại thời điểm đạt 2 giá trị cực trị bằng: Zj(t) = m(z) jyj(z), lực quán tính có giá trị thay đổi nên công của lực quán tính có giá trị bằng một nửa giá trị của lực nhân với chuyển vị tương ứng: 1 l 1 l = 2 -  2 2 U EI(z)[y j (z)] dz m(z) j y j (z)dz. 2 0 2 0
  17. Tìm dạng chính thứ j dưới dạng chuỗi: n y j (z) =  ai i (z) i=1 ai – các hệ số chưa biết cần xác định. n – số nguyên bất kỳ, i(z) – các hàm độc lập tuyến tính được chọn trước, thỏa mãn các điều kiện biên của hệ và càng phù hợp với dạng dao động càng tốt. l EI(z) n 1 l n = 2 -  2 2 U [ai i (z)] dz m(z) j [ai i (z)] dz 0 2 i=1 2 0 i=1 Thế năng toàn phần U là hàm của các hệ số chưa biết ai: U = U(a1, a2, a3, , an)
  18. Từ điều kiện cực tiểu của thế năng toàn phần: U = 0; (k = 1, 2, , n) ak Ta có hệ phương trình: U l n l n = EI(z)[ a (z)] (z)dz - 2 m(z)[ a (z)] (z)dz = 0   i i k j  i i k ak 0 i=1 0 i=1 l l = - 2 Cki EI(z) i (z) k (z)dz j m(z) i (z) k (z)dz. 0 0 Dễ dàng nhận thấy: Cki = Cik. Sau khi biến đổi hệ phương trình trên có dạng: Ck1a1 + Ck2a2 + + Cknan = 0; (k = 1, 2,3, ,n)
  19. Ck1a1 + Ck2a2 + + Cknan = 0; (k = 1,2,3, ,n) Đây là hệ phương trình đại số tuyến tính thuần nhất, để các nghiệm ai 0, tức là để tồn tại dao động, thì định thức các hệ số trong hệ phương trình phải bằng không: C11 C12 C1n C C C 21 22 2n = 0. Cn1 Cn2 Cnn Khai triển định thức ta được phương trình bậc n 2 đối với  j. Đó là phương trình tần số xác định tần số dao động riêng j.
  20. Nếu chọn nghiệm là dạng chính thứ j dưới dạng chuỗi với số một số hạng của chuỗi (n = 1). Phương trình tần số có dạng: l l = = 2 - 2 2 = C11a1 0 C11 EI(z)[ 1 (z)] dz j m(z) 1 (z)dz 0 0 0 l 2 EI(z)[ 1 (z)] dz 2 0  j = l ; 1 (z) 0 2 m(z) 1 (z)dz 0 Nếu chọn với hai số hạng của chuỗi (n = 2) thì: 2 C11C22 - C12 = 0
  21. CHƯƠNG 4: CÁC PHƯƠNG PHÁP TÍNH GẦN ĐÚNG 4.2.2. Dao động cưỡng bức: Biểu thức thế năng toàn phần của hệ sẽ là: 1 l 1 l U = EI(z)[y (z)]2 dz - m(z) 2 y2 (z)dz - 2 0 2 0 l n - - q(z)y(z)dz  Pj y(z j ), 0 j=1 q(z): Lực kích thích phân bố; Pj : Biên độ lực kích thích tập trung thứ j. Chọn nghiệm dưới dạng chuỗi như phần trước và thế vào phương trình trên:
  22. 1 l n 1 l n = 2 -  2 2 - U EI(z)[ ai i (z)] dz m(z) [ ai i (z)] dz 2 0 i=1 2 0 i=1 l n n n - - q(z) ai i (z)dz  Pj  ai i (z j ). 0 i=1 j=1 i=1 Từ điều kiện cực tiểu của thế năng toàn phần ta được hệ phương trình: U l n l n = -  2 - EI(z)[ ai i (z) k (z)dz m(z) [ ai i (z)] k ak 0 i=1 0 i=1 l n - - = = q(z) k (z)dz  Pj k (z j ) 0; k 1, 2, , n. 0 j=1
  23. Gọi: l l = - 2 Cki EI(z) i (z) k (z)dz m(z) i (z) k ; 0 0 l n = - - CkP q(z) k (z)dz  Pj k (z j ). 0 j=1 Lúc này hệ phương trình có dạng: Ck1a1 + Ck2a2 + + Cknan + CkP ; k = 1,2, ,n Giải hệ này ta tìm được các hệ số ai, tiếp đó sẽ tìm được chuyển vị động, nội lực động trong hệ dao động cưỡng bức. Khi tính với tải trọng tĩnh ta cho  = 0
  24. Ví dụ 3: Xác định gần đúng EI = const m tần số dao động cơ bản của dầm đơn giản có nhịp l và l mang khối lượng phân bố đều m như hình vẽ. Chọn n = 1 và dạng dao động của dầm (trùng với biểu thức chính xác của dao động chính) là: (z) = sin z 1 l l 2 EI(- sin z)2 dz l 2 l 4 EI  2 = 0 = , 1 l l 4m m(sin z)2 dz 0 l
  25. 9,8696 EI  = , 1/ s 1 l 2 m Kết quả trùng với giá trị chính xác của tần số 1. * Nếu chọn n = 1 và dạng dao động của dầm là đường đàn hồi do tải trọng phân bố đều gây ra: 3 3 4 4 4 1 (z) = A(l z - 2lz + z ) / l , A = ql / 24EI l 12 12 2 - + 2 EIA ( 3 z 4 zz )dz 2 0 l l EI 1 = = 97,56 l z 2z3 z 4 l 4m 3 - + 2 mA ( 3 4 ) dz 0 l l l
  26. 9,8772 EI  = , 1/ s 1 l 2 m So với kết quả chính xác, sai số 0,1%. * Nếu chọn n = 1 và dạng dao động của dầm là đường đàn hồi do tải trọng tập trung ở giữa dầm gây nên: z 4z (z) = B(3 - 3 ), 0 z l / 2 1 l l 3 trong đó: B = l3/48EI
  27. l / 2 z EIB2 (-24 )2 dz l 3 EI  2 = 0 = 98,7085 1 l / 2 3 l 4m 2 z - 4z 2 mB (3 3 ) dz 0 l l 9,9352 EI  = 1 l 2 m So với kết quả chính xác, sai số là 0,67%
  28. Ví dụ 4: Xác định gần đúng độ võng, mô men uốn tại tiết diện ở giữa dầm đơn EI = const Psin  t m giản có nhịp l chịu lực kích thích P(t) = Psint như trên l hình vẽ. Cho biết l = 2 m, E = 2.104 kN/cm2, I = 800 cm4,  = 900 1/s, m = 0,00001 kNs2/cm Chọn n = 1 và 1 = sin( z / l), ta có: y(z) = a sin z, 1 l 2 M (z) = -y (z) = a sin z. EI 1 l 2 l
  29. Khi z = l/2, tại tiết diện giữa nhịp ta có: EI 2 y(l / 2) = a ; M(l / 2) = a . 1 1 l 2 Phương trình xác định hệ số a1: C11a1 + C1P = 0. l 2 l 4 EI  2ml = 2 - 2 2 = - C11 EI[ 2 sin z] dz m(sin z) dz 3 , 0 l l 0 l 2l 2 l C P = -Psin ( ) = -P. 1 l 2 C1P P -4 a1 = - = = 60,94.10 P cm. C11 164,1
  30. Độ võng và mô men uốn tại tiết diện giữa dầm: -4 y(l / 2) = a1 = 60,94.10 cm Trị số chính xác: y(l/2) = 61,094.104 cm. 2 EI M(l / 2) = a = 240,582P kNcm. l 2 1 Trị số chính xác: M(l/2) = 250,1102 kNcm.
  31. Nếu chọn n = 2 và y(z) = a1sin( z/l)+a2sin(3 z/l): y(l / 2) = a1 - a2 , 2 9 2 3 M(z) = -EIy (z) = [a sin z + a sin z]EI, 1 l 2 l 2 l 2 l 2 9 2 M(l / 2) = -EIy (l / 2) = [a - a ]EI. 1 l 2 2 l 2 Hệ phương trình xác định các hệ số: a11 + a12 + C1P = 0, a21 + a22 + C2P = 0.
  32. l 2 l 4 EI l = - 2 - 2 2 = - 2 C11 EI( 2 sin z) dz m(sin z) dz 3 m , 0 l l 0 l 2l 2 l 2 9 2 3 l 3 = = - 2 C12 C21 EI 2 sin z 2 sin zdz msin z sin , 0 l l l l 0 l l l 9 2 3 l 3 81 4 EI l = 2 - 2 2 = - 2 C22 EI( 2 sin z) dz m(sin z) dz 3 m , 0 l l 0 l 2l 2 l C = -Psin = -P, 1P l 2 3 l C = Psin = P. 2 P l 2 Thay các giá trị này vào hệ phương trình trên ta tìm được các hệ số: -4 -4 a1 = 60,94.10 P cm. a2 = -0,128.10 P cm.
  33. Độ võng và mô men uốn tại tiết diện giữa dầm: -4 y(l / 2) = (a1 - a2 ) = 61,068.10 P cm. 2 9 2 M(l / 2) = [a - a ]EI = 245,13P kNcm. 1 l 2 2 l 2 Khi chọn chuỗi nghiệm với hai số hạng, kết quả nhạn được chính xác hơn khi chọn chuỗi nghiệm với một số hạng.
  34. CHƯƠNG 4: CÁC PHƯƠNG PHÁP TÍNH GẦN ĐÚNG 4.3. Phương pháp BUPVÔV – GALOOCKIN : 4.3.1: Dạng dao động riêng: Để giải phương trình vi phân biểu thị dao động tự do với dạng chính thứ j của hệ có số bậc tự do bằng vô cùng:  2 [EI(z)y (z)]- 2m(z)y (z) = 0. z 2 j j j Ta giả thiết nghiệm yj(z) được tìm dưới dạng chuỗi: n y j (z) =  ai i (z) i=1
  35. n y j (z) =  ai i (z) i=1 ai – các hệ số chưa biết cần xác định. n – số nguyên bất kỳ, i(z) – các hàm độc lập tuyến tính được chọn trước, thỏa mãn các điều kiện biên của hệ và càng phù hợp với dạng dao động càng tốt. 2 n n - 2 = 2 [EI(z) ai i (z)] j m(z) ai i (z) 0. z i=1 i=1 Phương trình này nghiệm đúng với mọi giá trị z và cũng nghiệm đúng khi nhân hai vế với một hàm k(z) bất kỳ. Sau khi nhân hai vế với k(z) và lấy tích phân theo chiều dài thanh ta thu được:
  36. l 2 n n  - 2 = 2 [EI(z)ai i (z)] j m(z)ai i (z) k (z).dz 0. 0 z i=1 i=1  Cho i =1, 2, 3, và khai triển ta thu được hệ phương trình chứa các ẩn số ai: Ck1a1 + Ck2a2 + + Cknan = 0; k = 1,2,3, ,n Với: l 2  = - 2 Cki 2 [EI(z) i (z)] k (z) j m(z) i (z) k (z)dz. 0 z  Khi các hàm i(z) và k(z) thỏa mãn các điều kiện biên của hệ thì Cki = Cik.
  37. C11 C12 C1n C C C 21 22 2n = 0. Cn1 Cn2 Cnn Khai triển định thức này ta được phương trình 2 bậc n đối với  j. Phương trình này gọi là phương trình tần số. Giải phương trình này ta tìm được i . Nếu chọn nghiệm đã chọn với n = 1, phương trình xác định ai có dạng: C11a1 = 0 C11 = 0
  38. l 2 [EI(z) (z)] (z)dz z 2 1 1  2 = 0 j l . 2 m(z) 1 (z)dz 0 Nếu 1(z) được chọn dưới dạng đa thức thì cần thỏa mãn điều kiện: 2 [EI(z) (z)] 0 z2 1
  39. CHƯƠNG 4: CÁC PHƯƠNG PHÁP TÍNH GẦN ĐÚNG 4.3. Phương pháp BUPVÔV – GALOOCKIN : 4.3.2: Dao động cưỡng bức: Khi hệ dao động cưỡng bức chịu lực kích thích tuần hoàn q(z, t) = q(z)sint. Phương trình vi phân dao động có dạng (không cản): 2 2 y(z,t) 2 y(z,t) [EI(z) ]+ m(z) = -q(z)sin t. z2 z 2 z 2 Biểu thị nghiệm phương trình dưới dạng: y(z,t)=y(z)sint 2 [EI(z)y (z)]- m(z) 2 y(z) = -q(z) z 2
  40. Thực hiện tương tự như trong dao động riêng, ta đưa phương trình động về dạng: Ck1a1 + Ck2a2 + + Cknan + CkP = 0 l 2  = - 2 Cki 2 [EI(z) i (z)] k (z) m(z) i (z) k (z)dz. 0 z  l = CkP q(z) k (z)dz. 0 Khi có lực kích thích tập trung tác dụng: n CkP =  Pj k (z j ). i=1
  41. Ví dụ 5: Xác định tần số dao h(z) ho z động cơ bản 1 của dầm côngxôn như hình vẽ. Dầm có bề rộng tiết diện ngang z không đổi, chiều cao tiết l diện ngang, khối lượng và mômen quán tính thay đổi theo quy luật: 3 ho z z h(z) = z, m(z) = m , I(z) = Io . l o l l 3 Dạng dao động của dầm i(z) có thể chọn là đường đàn hồi do các tải trọng khác nhau tác dụng trên dầm gây ra. Chọn dạng dao động là đường đàn hồi của dầm công xôn có tiết diện không đổi do tải trọng phân bố đều gây nên:
  42. h z z 4 h(z) o = - + z i (z) A(3 4 4 ) l l z Thay vào biểu thức tính l tần số dao động riêng: q l  2 z3 12z 2 4z z4 [EI A ]A(3 - + )dz  2 o 3 2 4 2 o z l l l l 1 = l z z z4 2 - + 2 mo A (3 4 4 ) dz 0 l l l 39,38EI  2 = o 6,275 EIo 1 4  = 1/ s l m 1 2 o l mo
  43. ho 6,275 EI h(z)  = o z 1 2 1/ s l mo z So sánh với kết quả chính xác: l 5,315 EI  = o q 1 2 1/ s l mo Sai số là 18,1%.
  44. h Nếu chọn dạng dao h(z) o z động là đường đàn hồi do tải trọng phân bố có dạng tam giác và hàm này thỏa z mãn điều kiện biên: l z (z) = (1- )2 1 l q l  2 z3 2 z [EI ](1- )2 dz  2 0 3 2 2 0 z l l l 30EIo 1 = = l z z l 4m - 4 o mo (1 ) dz 0 l l 5,477 EI  = o 1 2 1/ s Sai số 3% l m0
  45. 3 Nếu chọn dạng dao động 1(z)=(1-z/l) thỏa mãn các điều kiện biên: l 2 z3 6 z z [EI (1- )(1- )3 dz  2 o 3 2 2 0 z l l l l 33,6EIo 1 = = l z z l 4m - 6 o mo (1 ) dz 0 l l 5,7965 EI  = o 1 2 1/ s. l mo Sai số so với kết quả chính xác là 9%.
  46. Để có kết quả chính xác hơn, ta chọn nghiệm với hai số hạng: z 2 z 3 (z) = (1- ) 2 (z) = (1- ) 1 l l Ta thu được phương trình tần số có dạng: C C 11 12 = 0 C21 C22 Áp dụng biểu thức tính khi có lực kích thích: n CkP =  Pj k (z j ). i=1
  47. l  2 z 3 2 z z z  C = [EI ](1- )2 - 2m (1- )3 dz 11  2 o 3 2 1 o 0 z l l l l l  EI l = o - m  2 , l 3 o 1 30 l  2 z 3 6 z z z z  C = [EI (1- )](1- )3 - 2m (1- )6 dz 22  2 o 3 2 1 o 0 z l l l l l l  EI 6 l = o - m  2 . l 3 10 o 1 56 l  2 z 3 6 z z z z z  C = C = [EI (1- )](1- )2 - 2m (1- )2 (1- )3 dz 12 21  2 o 3 2 1 o 0 z l l l l l l l  EI 6 l = o - m  2 l 3 10 o 1 42
  48. Thay các giá trị vừa tìm được vào phương trình tần số ta thu được nghiệm nhỏ nhất của phương trình là: 5,319 EI  = o 1 2 1/ s l mo So sánh với kết quả chính xác, sai số là 0,1%. Khi chọn dạng dao động càng sát với dạng dao động của hệ và lấy số hạng trong chuỗi nghiệm càng nhiều thì kết quả càng chính xác.
  49. CHƯƠNG 4: CÁC PHƯƠNG PHÁP TÍNH GẦN ĐÚNG 4.4. Phương pháp thay thế khối lượng: Phương pháp thay thế khối lượng là một trong các phương pháp gần đúng có hiệu quả hay được sử dụng trong thực tế để tính toán dao động công trình. Nội dung của phương pháp này là dựa trên cơ sở đơn giản hóa sơ đồ khối lượng nhằm giảm số bậc tự do của hệ, nghĩa là thay khối lượng phân bố và các khối lượng tập trung bằng các khối lượng tập trung với số lượng ít hơn. Chia khối lượng phân bố thành nhiều khoảng rồi thay thế theo:
  50. * Tập trung khối lượng phân bố trên mỗi khoảng chia về trọng tâm khoảng chia. * Tập trung khối lượng phân bố trên mỗi khoảng về hai khối lượng đặt ở hai đầu doạn chia theo nguyên tác đòn bẩy
  51. Xét dầm đơn giản m như hình vẽ có khối EI lượng phân bố đều m và l độ cứng EI = const. a m.l /8 m.l /4 m.l /4 m.l /4 m.l /8 Ta có các sơ đồ thay thế khối lượng như hình l /4 vẽ. Tương ứng với mỗi b sơ đồ thay thế ta xác m.l /3 m.l /3 m.l /6 định được các tần số dao m.l /6 động riêng chính xác gần đúng. Ta có bảng liệt kê: l /3 c m.l /4 m.l /2 m.l /4 EI  = i i 2 l m l /2 l /2 d
  52. Hình b Hình c Hình d TT Hình a Kết quả Sai số Kết Sai số Kết Sai số % quả % quả % 1 9,87 9,865 0,05 9,86 0,1 9,789 0,7 2 39,48 39,2 0,7 38,2 3,24 - - 3 88,74 84,6 4,7 - - - - Khối lượng thay thế càng nhiều thì sai số kết quả tính càng nhỏ. Khi chỉ cần xác định tần số cơ bản i thì độ chính xác cũng đủ đáp ứng nếu tập trung khối lượng của dầm thành một khối lượng đặt ở giữa nhịp.
  53. CHƯƠNG 4: CÁC PHƯƠNG PHÁP TÍNH GẦN ĐÚNG 4.5. Phương pháp quy đổi về hệ có một khối lượng tương đương: * Thay hệ có bậc tự do hữu hạn hay vô hạn bằng hệ có bậc tự do bằng một tương đương: Hai hệ tương đương về động năng thì dao động cùng chung một tần số: K = K*. * Giả thiết trước đường đàn hồi của hệ thực khi dao động có dạng: y(z,t) = y(z).T(t) y&(z,t) = y(z).T&(t)
  54. * Chọn trước vị trí khối lượng thay thế M trên hệ tương đương. Thông thường nên chọn là vị trí tương ứng trên hệ thực nơi có chuyển vị lớn nhất khi hệ dao động. Nếu trên hệ thực có đặt các khối lượng tập trung thì nên đặt khối lượng thay thế M tại vị trí tương ứng với khối lượng tập trung lớn nhất. m(z)[y(z).T&(t)]2 m [y(z ).T&(t)]2 K =  dz +  j j , 2 2 M[y(a).T&(t)]2 K * = . 2 2 2 m(z)[y(z)] + mj [y(z j ] M = [y(a)]2
  55. Xác định gần đúng tần số 1 theo công thức tìm tần số dao động riêng của hệ dao động có bậc tự do bằng một. 1 1 = M aa Sai số của kết quả tính phụ thuộc vào việc ta chọn đường đàn hồi và vị trí đặt khối lượng thay thế tương đương M
  56. Ví dụ 6: Xác định gần đúng m tần số dao động riêng  1 EI của dầm đơn giản có khối l lượng phân bố đều m nhịp a l và độ cứng EI = const như M trên hình vẽ. l * Giả định đường đàn hồi của dầm khi dao động P = 1 f là đường đàn hồi khi dầm chịu lực P = 1 đặt ở giữa l/2 l/2 nhịp: 3l 2 z - 4z 3 l 3 y(z) = y( l ) y( l ) = f =  = 2 l 3 2 aa 48EI
  57. Khối lượng thay thế đặt tại giữa dầm: l 3l 2 z - 4z 3 m [ f ]dz l 3 M = 0 = 0,486ml f 2 Giá trị gần đúng của tần số cơ bản: 1 9,92 EI  = = 1/ s 1 l 3 l 2 m 0,486ml 48EI So với kết quả chính xác (9,87), sai số +0,5%.
  58. * Giả định đường đàn hồi của dầm khi dao động là đường đàn hồi của dầm khi chịu tải trọng phân bố đều với cường độ q = 1. 16 y(z) = y( l ) (l 3 z - 2lz 3 + z 4 ); 0 z l / 2 2 5l 4 Thực hiện tính toán tương tự ta tìm được: 9,76 EI  = 1/ s 1 l 2 m Sai số là -0,9%
  59. * Giả định đường đàn hồi khi dao động có dạng nửa bước sóng hình sin: z y(z) = y( l )sin 2 l 9,81 EI  = 1/ s (-0,7%) 1 l 2 m Khi giả định đường đàn hồi của dầm chịu lực tập trung đặt ở giữa nhịp và khối lượng thay thế M đặt ở giữa nhịp nơi có chuyển vị lớn nhất, kết quả tính toán chính xác hơn.
  60. Bây giờ xét một dầm mang khối lượng phân bố m và các khối lượng tập trung mj (j = 1, 2, 3, , n) có các liên kết ở hai đầu khác nhau. Giả định dạng dao động thứ nhât của dầm là đường đàn hồi khi dầm chịu lực tập trung tĩnh tại vị trí tương ứng với điểm đặt khối lượng thay thế M trên dầm tương ứng và gọi: yaa – chuyển vị tại tiết diện có hoành độ a khi dầm chịu tải trọng tập trung bằng đơn vị tác dụng tĩnh tại điểm đó (điểm đặt khối lượng thay thế M) yja, yza – chuyển vị tại các tiết diện có hoành độ zj và z khi dầm chịu lực đơn vị đặt tại tiết diện có hoành độ là a.
  61. Công thức tính khối lượng thay thế: l y y = za 2 + ja 2 M m ( ) dz mj ( ) 0 yaa j yaa Hay:M = .m.l +  j mj j l 3 y = aa EI Các trị số ,  được tra theo bảng 1 EI  = = 1 Ml3 /EI Ml3
  62. Ví dụ 7: Xác định tần số dao động riêng cơ bản của m dầm một nhịp có hai đầu 1 m ngàm mang khối lượng Z1= l/4 phân bố m và khối lượng tập trung m1 như hình vẽ. 7 2 Cho biết: EI = 1.10 kN/m , z1 = 2,5 m, l = 10m, 2 2 m = 0,153 kNs /m, m1 = 2,038 kNs /m. Áp dụng quy đổi về hệ có một khối lượng. Tra bảng:  = 192, = 0,37 .l = 0,37,10 = 3,7, z1 / l = 2,5 / 10 = 0,25 = 1 2 M = .m.l + 1.m1 =3,7. 0,153 + 0,25.2,038 = 1.0756 kNs /m
  63. Tần số dao động riêng cơ bản của dầm: EI 192 1 107  = = = 1336,1 1/ s 1 Ml3 1,0756 103
  64. CHƯƠNG 4: CÁC PHƯƠNG PHÁP TÍNH GẦN ĐÚNG 4.6. Phương pháp giải đúng dần: Nhược điểm của 5 phương pháp đã nghiên cứu là khó đánh giá được sai số của kết quả tính khi chưa biết giá trị chính xác của tần số dao động riêng. Phương pháp giải đúng dần sẽ khắc phục được nhược điểm trên và cho phép xác định gần đúng giá trị của tần số dao động riêng càng sát với giá trị chính xác nếu càng thực hiện hiều lần tính toán. Tuy nhiên nhược điểm của phương pháp này là quá dài và khối lượng tính toán quá lớn.
  65. * Nội dung phương pháp: m(z) Xét dầm mang khối m m lượng phân bố m(z) và các 1 j mn khối lượng tập trung mj 2 như hình vẽ. q(z) = m(z) k yk(zj) Nếu dạng chính thứ k của dao động xác định bằng phương trình y(z) 2 zj = mj(z) k yk(zj) thì lực quán tính tại các yk(z) khối lượng trên hệ là: 2 q(z) = m(z)k yk (z), 2 z j = mjk yk (z j ).
  66. Nếu giảm lực quán m(z) tính tác dụng trên hệ 2 xuống  k lần: m1 mj mn (1) q = m(z)yk (z), q(z) = m(z)2 y (z ) (1) k k j z j = mj yk (z j ) Do vậy phương trình đường đàn hồi của hệ z = m (z)2 y (z ) (1) j j k k j y k(z) do các lực quán yk(z) tính gây ra cũng giảm q(1)(z) = m(z) y (z ) 2 k j xuống  k lần so với ban đầu. y (z)  = k k (1) . (1) Z(1) (z) = m y (z ) yk (z) y k j j k j (z)
  67. y (z)  = k k (1) . yk (z) Vì hàm yk(z) chưa biết nên trong lần tính gần đúng thứ nhất cần giả thiết dạng dao động chính là hàm k(z) nào đó và xác định giá trị gần đúng thứ nhất của tần số dao động riêng theo công thức: (z)  (1) = k k (1) . k (z) (1) trong đó: k(z) – đường đàn hồi của hệ do các lực quán tính
  68. (2) Lại chọn đường đàn hồi thứ hai có dạng k(z) : (1) (z)  2 = k k (2) . k (z) Qúa trình tính gần đúng cứ tiếp tục cho đến khi (i+1) (i) trị số k xấp xỉ bằng k trong lần tính thứ i thì (i+1) dừng lại. Lúc này k ứng với các điểm khác nhau trên hệ đều có cùng một giá trị.
  69. Ví dụ 8: Xác định tần số cơ bản dao động m riêng của dầm có khối EI lượng phân bố đều l như hình vẽ. * Chia dầm thành sáu m1 m2 m3 m4 m5 khoảng bằng nhau và thay thế khối lượng phân bố đều bằng 5 khối lượng tập trung như hình vẽ. m1 = m2 = m3 = m4 = m5 = ml/6
  70. * Chọn dạng dao động là m đường đàn hồi 1(z) do EI trọng lượng của các khối l lượng mi tập trung gây ra. Do tính đối xứng của m1 m2 m3 m4 m5 dầm ta có: 1 (z1 ) = 1 (z5 ), P = 1 1 (z2 ) = 1 (z4 ), ij i 1 (z3 ) = 1max. z zi l - zi
  71. * Xác định các chuyển vị ta sử dụng phương trình ảnh hưởng của chuyển vị: l - z z  = i j [l 2 - (l - z )2 - z 2 ]; z z ij 6 EI l j i i j ij chuyển vị tại tiết diện đặt khối lượng mi do Pj = 1 đặt tại vị trí khối lượng mj gây ra. Tọa độ Tọa độ đặt lực P = 1 tìm chuyển z1 = l/6 z2=2l/6 z3=3l/6 z4=4l/6 z5=5l/6 ij vị z1 50 76 78 62 34 300 z2 76 108 138 112 62 496 z3 78 138 162 138 78 594
  72. Tọa độ tìm Tọa độ đặt lực P = 1 chuyển vị z1 = l/6 z2=2l/6 z3=3l/6 z4=4l/6 z5=5l/6 ij z1 50 76 78 62 34 300 z2 76 108 138 112 62 496 z3 78 138 162 138 78 594 l 3 mgl mgl4 (z ) = (z ) = 300 = 50 , 1 1 1 5 6 5 EI 6 6 5 l 3 mgl mgl4 (z ) = (z ) = 496 = 82,6666 , 1 2 1 4 6 5 EI 6 6 5 l 3 mgl mgl4 (z ) = 594 = 99 . 1 3 6 5 EI 6 6 5
  73. (1) Tìm các tải trọng gây ra đường đàn hồi 1(z): ml mgl4 m2 gl5 Z (1) = Z (1) = m (z ) = 50 = 50 , 1 5 1 1 1 6 6 5 EI 6 6 EI ml mgl4 m2 gl 5 Z (1) = Z (1) = m (z ) = 82,6666 = 82,666 , 2 4 2 1 2 6 6 5 EI 6 6 EI ml mgl4 m2 gl 5 Z (1) = m (z ) = 99 = 99 . 3 3 1 3 6 6 5 EI 6 6 EI
  74. Độ võng lớn nhất ở giữa dầm : (1) (1) (1) (1) 1 (z3 ) = 231Z1 + 232Z2 +33Z3 = l 3 m2 gl5 = [(78 50 2) + (138 82,6666 2) + (162 99)] 6 5 EI 66 EI m2 gl5 =7775,666 610 (EI)2 4 10 2 (1) 1 (z3 ) mgl 6 (EI) 9,95 EI 1 = (1) = 99 5 2 5 = 2 1 (z3 ) 6 EI 7775,666m gl l m So với kết quả tính chính xác sai số là 0,81% Tiếp tục thực hiện ta thu được kết quả chính xác hơn.
  75. Bây giờ ta chọn đường đàn hồi của dầm là do khối lượng tập trung đặt ở giữa dầm gây ra: M = ml/2 ml l 3 ml2 (z ) = = 1 1 2 48EI 96EI Tải trọng tập trung tại vị trí đặt tại vị trí có độ võng max ml ml4 m2l 5 Z (1) = M (z ) = = 1 1 1 2 96EI 192EI m2l 5 l 3 m2l 8 (1) = = . 1 192EI 48EI 192 48(EI)2
  76. (z) m2l 4 192 48(EI)2  (1) = 1 = = 1 (1) 2 8 1 96EI m l 9,8 EI = 1/ s l 2 m Tiếp tục thực hiện nhiều lần sẽ tìm được giá trị chính xác.
  77. CHƯƠNG 4: CÁC PHƯƠNG PHÁP TÍNH GẦN ĐÚNG 4.7. Phương pháp sai phân: Chuyển phương trình vi phan thành phương trình sai phân hữu hạn. Trường hợp dao động tự do không có lực cản, phương trình vi phân cho thanh có khối lượng phân bố đều m, EI = const: m 2 y IV - k4 y = 0; k4 = EI Đặt z = lx dz = ldx, *k4 = ml42/EI d 4 y -* k4 .y = 0 dx 4
  78. Chia dầm thành n đoạn bằng nhau có độ dài: l 1 z = = l x x = n n Khi chuyển sang sai phân ta thay vi phân d bằng sai phân và: y = yi - yi-1 2 y = y(i+1) - y(i) = yi+1 - 2 yi + yi-1 3 y = yi+1 - 3yi + 3yi-1 + yi-2 4 y = yi+2 - 4 yi+1 + 6 y1 - 4 yi-1 + yi-2
  79. Thay vào ta thu được phương trình sai phân: k4 - + - - + = yi-2 4 yi-1 6 y1 4 yi+1 yi+2 0 n Mỗi điểm chia ta được một phương trình. Cộng các điều kiện biên ta giải và tìm được k và suy ra 