BÀI TOÁN PHẲNG TRONG TỌA ĐỘ DESCARTES docx

z y x z y CHƯƠNG 6 : BÀI TOÁN PHẲNG TRONG TỌA ĐỘ DESCARTES 6.1. HAI TRƯỜNG HỢP CỦA BÀI TOÁN PHẲNG I. Khái niệm : Trong nhiều bài toán kỹ thuật, vật thể chịu lực chỉ gây nên biến dạng hay ứng suất trong 1 mặt phẳng (Mặt phẳng này được qui ước là mặt phẳng oxy). Các bài toán này được gọi là các bài toán phẳng. Bài toán phẳng chia ra 2 loại : 1. Bài toán ứng suất phẳng : Nếu chỉ tồn tại ứng suất trong mặt phẳng xoy. 2.Bài toán biến dạng phẳng : Nếu chỉ tồn tại biến dạng trong mặt phẳng xoy. Hai bài toán này khác nhau về mặt vật lý song rất giống nhau về mặt toán học. Giải bài toán phẳng về mặt toán học được đơn giản rất nhiều so với bài toán không gian. II. Bài toán ứng suất phẳng : Xét những mặt phẳng, ví dụ tấm tường, đĩa mỏng chịu lực phân bố đều trên bề dày tấm và song song với mặt trung bình như hình vẽ. Ta nhận thấy mặt bên của tấm không có tải trọng, ứng suất là hằng theo bề dày. Do đó điều kiện của bài toán sẽ là : σ z = T xz = T yz = 0 (a) Mặt khác, biến dạng dài theo phương bề dày là tự do nên : ε z ≠ 0 (b) Các điều kiện (a), (b) là định nghĩa của bài toán ứng suất phẳng. Ân số của bài toán gồm có: Các ứng suất : σ x , σ y , T xy . Các biến dạng : ε x , ε y , γ xy , ε z ≠ 0. Theo định luật Hooke, từ (a) ta có : 42 z 1 y z x 1 γ xz =γ yz = 0 ; ε y = E 1 (σ y - µσ x ) ε x = E 1 (σ x - µσ y ) ; ε z =- E µ (σ x + σ y ) (c) γ xy = G T xy = E )1(2 µ + T xy Từ biểu thức (c) ta có các biến dạng đều tính theo 3 ẩn số ứng suất là σ x , σ y , T xy với E, µ là 2 hằng số đàn hồi của vật liệu. III. Bài toán biến dạng phẳng : Khi tính những vật thể hình lăng trụ, có chiều dài lớn chịu tải trọng không đổi theo chiều dài, ví dụ đập chắn, tường chịu áp lực, đường ống dẫn, vỏ hầm ta thường xét 1 đoạn vật thể có chiều dài bằng 1 đơn vị. Như thế, bài toán đối với vật thể lăng trụ trở thành bài toán tấm phẳng như biểu diễn trên hình vẽ sau : Nhận xét tấm bị kẹp giữa chiều dài của vật thể nên không thể có biến dạng dài theo phương bề dày z, và mặt bên của tấm sẽ chịu những áp lực pháp tuyến theo phương z. Do đó, điều kiện của bài toán đối với tấm trong trường hợp đang xét sẽ là : ε z = γ xz = γ yz = 0 (d) và σ z ≠ 0 (e) Các điều kiện (d), (e) là định nghĩa của bài toán biến dạng phẳng. Ẩn số của bài toán gồm có: Các ứng suất : σ x , σ y , T xy , σ z ≠0 Các biến dạng : ε x , ε y , γ xy . Theo định luật Hooke, từ (d) ta có : - Các ứng suất tiếp Txz = Tyz = 0 - Còn ứng suất pháp σ z sẽ được tìm từ biểu thức ε z = 0 ε z = [ ] )( 1 yxx E σσσ µ +− = 0 Vậy σ y = µ(σ x + σ y ). Quan hệ giữa các ứng suất và các biến dạng còn sẽ là : 43 ε x = [ ] )( 1 zyx E σσσ µ +− = [ ] )( 1 yxy E σσσ µ +− ε x =       − − − ) 1 1 2 yx E σσ µ µµ Tương tự ε y =       − − − ) 1 1 2 xy E σσ µ µµ (*) γ xy = E )1(2 µ + T xy Đặt E 1 = 2 1 µ − E ; µ 1 = µ µ −1 (g) (*)⇔ ε x = 1 1 E (σ x - µ 1 σ y ) ; ε y = 1 1 E (σ y - µ 1 σ x ) ; (f) γ xy = E )1(2 µ + T xy = 1 1 )1(2 E µ + T xy IV. So sánh và kết luận chung : 1. Trong cả 2 bài toán phẳng, các ẩn số chính về ứng suất và về biến dạng là như nhau : σ x , σ y , T xy , ε x, ε y , γ xy . → Những ứng suất hay biến dạng còn lại đều có thể biểu diễn qua các ẩn số chính. 2. Quan hệ giữa các ứng suất hay biến dạng theo (c) hay (f) là hoàn toàn tương tự như nhau, sự khác nhau chỉ thể hiện ở chỗ : - Trong bài toán ứng suất phẳng ta dùng các hằng số đàn hồi E, µ còn trong bài toán biến dạng phẳng ta dùng các hằng số đàn hồi E 1 , µ 1 theo cách đặt (g). 3. Do sự giống nhau về mặt toán học như vậy nên phép giải của 2 bài toán hoàn toàn như nhau. 6.2. CÁC PHƯƠNG TRÌNH CƠ BẢN TRONG BÀI TOÁN PHẲNG 1. Về mặt tĩnh học : Phương trình cân bằng Cauchy : y Tyx x x ∂ ∂ + ∂ ∂ σ + f x = 0 y y x Txy ∂ ∂ + ∂ ∂ σ + f y = 0 (6.1) 2. Về mặt hình học : Phương trình biến dạng Cauchy : ε x = x u ∂ ∂ ; ε y = y v ∂ ∂ ; (6.2) γ xy = x u ∂ ∂ + y v ∂ ∂ . 44 Các biến dạng phải thỏa mãn điều kiện liên tục của biến dạng, trong bài toán phẳng điều kiện này chỉ còn 1 phương trình : yx xy x y y x ∂∂ ∂ = ∂ ∂ + ∂ ∂ γ ε ε 2 2 2 2 2 (6.3) 3. Về mặt vật lý : Phương trình định luật Hooke. a. Biểu thức biến dạng qua ứng suất : ε x = E 1 (σ x - µσ y ) ε y = E 1 (σ y - µσ x ) (6.4) γ xy = E )1(2 µ + T xy b. Biểu thức ứng suất qua biến dạng : σ x = 2 1 µ − E (ε x + µε y ) σ y = 2 1 µ − E (ε y + µε x ) (6.5) T xy = )1(2 µ + E γ xy Nếu giải bài toán biến dạng phẳng, chỉ cần thay E, µ bằng E 1 , µ 1 . Hệ tám phương trình độc lập trên, chứa 8 ẩn số là ba ứng suất, ba biến dạng và hai chuyển vị là một hệ khép kín, cho phép ta giải được bài toán. 4. Các điều kiện biên : a. Điều kiện biên tĩnh học : σ x l + T yx m = ∗ x f T xy l + σ y m = ∗ y f (6.6) b. Điều kiện biên động học : Trên bề mặt S của vật thể cho trước các chuyển vị u o , v o hay các đạo hàm của các chuyển vị theo các biến số tọa độ. Nghiệm chuyển vị của bài toán phải thỏa mãn điều kiện : u s = u o ; v s = v o . 6.3. PHÉP GIẢI BÀI TOÁN PHẲNG THEO ỨNG SUẤT - HÀM ỨNG SUẤT AIRY I. Phép giải theo ứng suất : - Chọn ẩn số chính là các ứng suất : σ x , σ y , T xy . Các ứng suất này phải thỏa mãn phương trình cân bằng (6.1) . y Tyx x x ∂ ∂ + ∂ ∂ σ = - f x 45 y y x Txy ∂ ∂ + ∂ ∂ σ = - f y Nghiệm của (6.1) sẽ là tổng của nghiệm tổng quát phương trình thuần nhất (6.8) y Tyx x x ∂ ∂ + ∂ ∂ σ = 0 y y x Txy ∂ ∂ + ∂ ∂ σ = 0 (6.8) và nghiệm riêng của phương trình (6.9) y Tyx x x ∂ ∂ + ∂ ∂ σ = - f x y y x Txy ∂ ∂ + ∂ ∂ σ = - f y (6.9) - Nghiệm riêng của phương trình (6.8) tìm được không khó khăn, nó phụ thuộc vào dạng cụ thể của các lực thể tích. Ví dụ nghiệm riêng có thể lấy là : * σ x = 0 ; σ y = 0 ; T xy = -Px khi f x = 0 ; f y = P = hằng số. * σ x = 2 2 ax + bx ; σ y = T xy = 0 khi f x = ax + b ; f y = 0 * σ x = 0 ; σ y = -a 6 3 y ; T xy = 2 2 axy khi f x = axy , f y = 0. II. Hàm ứng suất Airy : Để giải hệ (6.1) ta đưa ra một hàm ẩn mới gọi là hàm ứng suất Airy. Xét hệ phương trình phương trình vi phân thuần nhất (6.8): 0 )8.6(0 = ∂ ∂ + ∂ ∂ = ∂ ∂ + ∂ ∂ y y x Txy y Tyx x x σ σ Điều kiện cần và đủ cho biểu thức p(x,y)dx + q(x,y)dy = du(x,y) tức p(x,y)dx + q(x,y)dy là vi phân toàn phần của 1 hàm u(x,y) nào đó thì giữa p và q phải có quan hệ : x q y p ∂ ∂ = ∂ ∂ . - Phương trình thứ (1) của hệ (6.8) ⇔ y Tyx x x ∂ ∂ −= ∂ ∂ σ Tức (σx.dy - Txy.dx) là vi phân toàn phần của 1 hàm A(x,y) nào đó. Nên ta có quan hệ σ x = y A ∂ ∂ ; T yx = - x A ∂ ∂ (a) Tương tự, phương trình thứ 2 : x Txy y y ∂ ∂ −= ∂ ∂ σ ⇒ (σy.dx - Txy.dy) là vi phân toàn phần của1 hàm B(x,y) nào đó : 46 → Ta có quan hệ : σ y = x B ∂ ∂ ; T xy = - y B ∂ ∂ (b) So sánh (a) và (b) ta có : x A ∂ ∂ = y B ∂ ∂ (c) ⇒ (A.dy + B.dx) là vi phân toàn phần của 1 hàm ϕ(x,y) nào đó : → Ta có quan hệ : A = y ∂ ∂ ϕ ; B = x ∂ ∂ ϕ (d) Thay (d) vào (a) và (b) ta có: σ x = 2 2 y ∂ ∂ ϕ ; σ y = 2 2 x ∂ ∂ ϕ ; T xy = - yx ∂∂ ∂ ϕ 2 (6.10) Hàm ϕ(x,y) : Gọi là làm ứng suất Airy, là hàm để giải bài toán phẳng theo ứng suất. III. Phương trình hàm ứng suất Airy : - Trong chương 5 ta có hệ phương trình (5.5) Beltrmi là hệ phương trình giải bài toán đàn hồi theo ứng suất đã tổng hợp các điều kiện về mặt tĩnh học, hình học, và vật lý của môi trường. Sử dụng (5.5) để tính cho biểu thức ứng suất phẳng. (1 + µ)∇ 2 σ x + 2 2 x S ∂ ∂ = 0 + (1 + µ)∇ 2 σ y + 2 2 y S ∂ ∂ = 0 (1 + µ)∇ 2 σ z + 2 2 z S ∂ ∂ = 0 (1+µ)∇ 2 S +∇ 2 S = 0 ⇔ ∇ 2 S = 0 Với S = σ x + σ y + σ z . Vì trong bài toán ứng suất phẳng σ z =0 nên S= σ x + σ y Trong bài toán biến dạng phẳng : S= σ x + σ y + σ z = σ x + σ y +µ(σ x + σ y ) =(1+µ)(σ x + σ y ). Nên trong bài toán đàn hồi phẳng ta đều có : ∇ 2 S = ∇ 2 (σ x + σ y ) = 0 (6.11) (6.11) : Phương trình LêVy. Thay các ứng suất bởi hàm ϕ thay (6.10) vào (6.11) ta có : 0 2 2 2 2 2 2 2 2 =         ∂ ∂ + ∂ ∂         ∂ ∂ + ∂ ∂ xyyx ϕϕ ⇔ 0 yyx 2 x 4 4 22 4 4 4 = ∂ ϕ∂ + ∂∂ ϕ∂ + ∂ ϕ∂ (6.12) ⇔ ∇ 2 (∇ 2 ϕ) = ∇ 4 ϕ = 0 (6.13) 47 Phng trỡnh (6.13) : phng trỡnh trựng iu hũa. Hm = (x,y) : l hm trựng iu hũa . Kt lun : - Bi toỏn n hi phng gii theo ng sut dn n vic gii phng trỡnh (6.12) sau ú tỡm cỏc ng sut theo (6.10). + Nu f x , f y 0 Cng thờm cỏc nghim riờng. - Theo (6.10) : Vic thờm hay bt hm mt lng A+ Bx+Cy thỡ cỏc ng sut khụng thay i. - Cỏc h s tớch phõn c xỏc nh theo iu kin biờn tnh hc : = x fm yx l y 2 2 2 = y fm x l yx 2 22 (6.14) Nu (6.13) xỏc nh cỏc hng s tớch phõn thỡ cỏc ng sut theo (6.10); (6.12) & (6.14) hon ton khụng liờn quan n cỏc h s n hi ca vt liu. Nhng bi toỏn nh th l bi toỏn cú liờn kt bờn ngoi tnh nh. nh lý LeVy-Michell : Trong biu thc n hi phng tnh nh, chu cỏc ngoi lc tỏc ng trờn biờn thỡ s phõn b ng sut khụng ph thuc vo cỏc hng s n hi v nh nhau i vi tng c cỏc vt liu. + ởnh lyù õổồỹc sổớ duỷng laỡm cồ sồớ cho 1 phổồng phaùp thổỷc nghióỷm coù tón laỡ phổồng phaùp õaỡn họửi. 6.4. IU KIN BIấN CA HM NG SUT AIRY. Vic gii bi toỏn phng theo ng sut rỳt li thnh vic gii phng trỡnh trựng iu hũa (6.12). Nghim ca phng trỡnh ny l hm ng sut phi tha món iu kin biờn. += += ymTxylf Txyml.xf y x (6.15) Xột trng hp fx = fy = 0 Thay (6.10) vo (6.11) ta cú m yx l y f x = 2 2 2 . m x l yx f y 2 22 . = (6.16) Theo (H.6.3) ta cú : l = cos(n, x) = cos(90 0 + ) = - sin = - ds dy 48 m = cos(n, y) = cosβ = ds dx (6.15) ⇔ ds dy y f x . 2 2 ∂ ∂ −= ∗ ϕ - ds dx yx . 2 ∂∂ ∂ ϕ = - ds dy yy .         ∂ ∂ ∂ ∂ ϕ - ds dx yx .         ∂ ∂ ∂ ∂ ϕ = -         ∂ ∂ yds d ϕ . (6.17) ds dy yx f y . 2 ∂∂ ∂ += ∗ + ds dx x . 2 2 ∂ ∂ ϕ =       ∂ ∂ xds d ϕ . Lấy điểm so bất kỳ trên chu tuyến làm gốc : (6.17) ⇔ ( ) S S x XdsfA y =−= ∂ ∂ ∫ ∗ 0 ϕ ( ) S S y YdsfB x =−= ∂ ∂ ∫ ∗ 0 ϕ (6.18) Trong đó : A&B : Các hệ số tùy ý, biểu diễn giá trị của đạo hàm 0 0 , S S xy ∂ ∂ ∂ ∂ ϕϕ của chu vi . X (S) , Y (S) : Ký hiệu mang ý nghĩa tĩnh học sẽ nói đến dưới đây. Để rõ ràng ta đưa ra sự tương tự như sau : Thay chu vi vật thể khảo sát bằng thanh có cùng dạng và cắt tại điểm S 0 (H.6.4). Tại đó ta đặt các lực : A // S 0 x B // S 0 y Và ngẫu lực C như hình vẽ Như vậy : X (S) & Y (S) : Chính là tổng hình chiếu của các ngoại lực tác dụng lên đoạn S 0 S chiếu lên trục x & y. + Nếu chúng ta lấy trục t ≡ trục tiếp tuyến ngoài tại điểm S n ≡ pháp tuyến ngoại tại điểm S. Thì : = ∂ ∂ n ϕ N (S) (6.19) st ∂ ∂ = ∂ ∂ ϕϕ = Q (S) (6.20) N (S) : Lực dọc cũng tại điểm S của thanh, được xem là dương → nếu là lực kéo. Q (S) : Lực cắt tại điểm s của thanh. So sánh quan hệ giữa nội lực là moment uốn và lực cắt trong sức bằng vật liệu: 49 x y y z o P x L 2 t + 2 t − = ds dM Q (s) = ds dϕ Q (s) ⇒ ϕ = M (6.21) M (s) : Moment của lực đặt trên đoạn S 0 S của thanh đối với điểm s. Vậy tại điểm trên chu tuyến của vật thể ta có thể xác định giá trị của hàm ứng suất ϕ(x,y) và các đạo hàm theo phương pháp tuyến n ∂ ∂ ϕ tại các điểm ở trên chu vi theo trọng đã cho dựa vào công thức (6.21) và (6.19) , quá trình ///đó giống như tìm moment uốn S lực dọc gây ra bởi tải trọng cho trước trên chu vi nếu tưởng tượng chu vi đó là ////mà cắt ra tại 1 tải diện bất kỳ. ϕ có dạng bất kỳ : Chuỗ Taylor, Furiê, hàm phức, chuổi đặc biệt. ⇒ ϕ có dạng đa thức. 6.5. HÀM ỨNG SUẤT DƯỚI DẠNG ĐA THỨC Việc giải bài toán phẳng theo ứng suất là tìm một hàm ứng suất ϕ thỏa mãn 2 yêu cầu : - Phương trình trùng điều hòa - Điều kiện biên + Tính ứng suất trên tấm công chịu lực tập trung đặt tại đầu tự do như hình vẽ 1. Dạng hàm ϕ + Theo kết quả ở sức bền vật liệu: σ x = y J M Z Z . theo hàm ϕ : σ x = 2 2 y ∂ ∂ ϕ ϕ(x,y) = ax 2 + bxy + cy 2 + dx 3 + cx 2 y + fxy 2 + gy 3 + hx 4 + ix 3 y + ix 2 y 2 + kxy 3 + ly 4 . (a) ϕ phải thỏa mãn phương trình trùng điều hòa : 4 4 x ∂ ∂ ϕ + 22 4 yx ∂∂ ∂ ϕ + 4 4 y ∂ ∂ ϕ = 0 4 4 x ∂ ∂ ϕ = h ; 22 4 yx ∂∂ ∂ ϕ = j ; 4 4 y ∂ ∂ ϕ = l. → h + 2j + l = 0 50 ⇒ ϕ là hàm đa thức bậc 4 đối với x, y → h = j =l = 0 (1) σx = 2 2 x ∂ ∂ ϕ = 2c + 2fx + 6gy + 6kxy. σy = 2 2 x ∂ ∂ ϕ = 2c + 6dx + 6ey + 6ixy. (b) Txy = - yx ∂∂ ∂ ϕ 2 =-(b + 2ex+ 2fy + 3ix 2 + 3ky 2 2. Các điều kiện : Xét điều kiện biên theo ứng suất : * Biên trên (y = [ ] Lx t ,0; 2 ∀ : Txy = 0 , (c) σy = 0 (d) * Biên dưới (y =- [ ] Lx t ,0; 2 ∀ : Txy = 0 , (e) σy = 0 (f) Từ (c) & (e) ta có : 2a +6dx +2e( 2 t )+6ix( 2 t ) = 0 2a + 6dx - 2e 2 t - 6ix 2 t = 0 ⇒ e = i = 0 e = i=f=0 (2) Từ (d) & (f) ta có : a=d=0 (3) 0 0 2 33 2 22 0 2 33 2 22 2 2 2 2 =⇒          =             −++−+− =             ++++− f t kix t fexb t kix t fexb * Biên trái (x = 0, ∀y       +− 2 , 2 tt ) ta có : σx= 0 (g) pdFTxy t t = ∫ + − 2 2 . (h) Từ (g) ⇒ c = g = 0 (5) ⇒ Txy = - (- 4 3 kt 2 + 3ky 2 ) = 4 3 kt 2 - 3ky 2 ⇒ 2 2 2 2 2 2 3222 3 3 4 3 )3 4 3 ( t t t t t t kyyktdykyktTxydF − + − − ∫ ∫       −=−= δδ 51 ⇒ b + 3 4 2 kt =0 ⇒ b = - 2 4 3 kt (4) . : BÀI TOÁN PHẲNG TRONG TỌA ĐỘ DESCARTES 6.1. HAI TRƯỜNG HỢP CỦA BÀI TOÁN PHẲNG I. Khái niệm : Trong nhiều bài toán kỹ thuật, vật thể chịu lực chỉ gây nên biến dạng hay ứng suất trong 1 mặt phẳng. (Mặt phẳng này được qui ước là mặt phẳng oxy). Các bài toán này được gọi là các bài toán phẳng. Bài toán phẳng chia ra 2 loại : 1. Bài toán ứng suất phẳng : Nếu chỉ tồn tại ứng suất trong mặt phẳng. xoy. 2 .Bài toán biến dạng phẳng : Nếu chỉ tồn tại biến dạng trong mặt phẳng xoy. Hai bài toán này khác nhau về mặt vật lý song rất giống nhau về mặt toán học. Giải bài toán phẳng về mặt toán học