TOÁN HỌC CƠ SỞ Các ứng dụng của Lý Thuyết Đàn Hồi đòi hỏi sự hiểu biết về nhiều lĩnh vực toán học khác nhau. Bản thân Lý Thuyết Đàn Hồi được xây dựng trên cơ sở ứng dụng nhiều đại lư
zyxzyCHƯƠNG 6 : BÀI TOÁN PHẲNG TRONG TỌA ĐỘ DESCARTES6.1. HAI TRƯỜNG HỢP CỦA BÀI TOÁN PHẲNGI. 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 = Txz = Tyz = 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, Txy.Các biến dạng : εx, εy, γxy, εz ≠ 0.Theo định luật Hooke, từ (a) ta có :42 z1yzx1γxz =γyz = 0 ; εy = E1(σy - µσx)εx =E1(σx - µσy) ; εz =-Eµ(σx + σy) (c)γxy = GTxy = E)1(2µ+Txy 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, Txy 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, Txy, σz≠0Cá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 = [ ])(1yxxEσσσ µ+− = 0Vậ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 = [ ])(1zyxEσσσ µ+− = [ ])(1yxyEσσσ µ+−εx = −−−)112yxEσσµµµTương tự εy = −−−)112xyEσσµµµ(*)γxy = E)1(2µ+TxyĐặt E1 = 21µ−E ; µ1 = µµ−1(g)(*)⇔ εx = 11E(σx - µ1σy) ;εy = 11E(σy - µ1σx) ; (f)γxy = E)1(2µ+Txy = 11)1(2Eµ+TxyIV. So sánh và kết luận chung :1. Trong cả 2 bài tốn phẳng, các ẩn số chính về ứng suất và về biến dạng là như nhau : σx, σy, Txy, ε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à hồn tồn tương tự như nhau, sự khác nhau chỉ thể hiện ở chỗ :- Trong bài tốn ứng suất phẳng ta dùng các hằng số đàn hồi E, µcòn trong bài tốn biến dạng phẳng ta dùng các hằng số đàn hồi E1, µ1 theo cách đặt (g).3. Do sự giống nhau về mặt tốn học như vậy nên phép giải của 2 bài tốn hồn tồn như nhau.6.2. CÁC PHƯƠNG TRÌNH CƠ BẢN TRONG BÀI TỐN PHẲNG 1. Về mặt tĩnh học : Phương trình cân bằng Cauchy :yTyxxx∂∂+∂∂σ + fx = 0yyxTxy∂∂+∂∂σ + fy = 0 (6.1)2. Về mặt hình học : Phương trình biến dạng Cauchy :εx = xu∂∂;εy = yv∂∂; (6.2)γxy = xu∂∂+yv∂∂.44 Cỏc bin dng phi tha món iu kin liờn tc ca bin dng, trong bi toỏn phng iu kin ny ch cũn 1 phng trỡnh : yxxyxyyx=+22222(6.3)3. V mt vt lý : Phng trỡnh nh lut Hooke.a. Biu thc bin dng qua ng sut :x = E1(x - ày)y = E1(y - àx) (6.4)xy = E)1(2à+Txyb. Biu thc ng sut qua bin dng :x = 21àE(x+ ày)y = 21àE(y+ àx) (6.5)Txy = )1(2à+ExyNu gii bi toỏn bin dng phng, ch cn thay E, à bng E1, à1. H tỏm phng trỡnh c lp trờn, cha 8 n s l ba ng sut, ba bin dng v hai chuyn v l mt h khộp kớn, cho phộp ta gii c bi toỏn.4. Cỏc iu kin biờn :a. iu kin biờn tnh hc :xl + Tyxm = xfTxyl + ym = yf(6.6)b. iu kin biờn ng hc :Trờn b mt S ca vt th cho trc cỏc chuyn v uo , vo hay cỏc o hm ca cỏc chuyn v theo cỏc bin s ta . Nghim chuyn v ca bi toỏn phi tha món iu kin : us = uo ; vs= vo .6.3. PHẫP GII BI TON PHNG THEO NG SUT - HM NG SUT AIRYI. Phộp gii theo ng sut :- Chn n s chớnh l cỏc ng sut : x, y, Txy.Cỏc ng sut ny phi tha món phng trỡnh cõn bng (6.1) .yTyxxx+ = - fx45 yyxTxy∂∂+∂∂σ = - fy 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)yTyxxx∂∂+∂∂σ = 0yyxTxy∂∂+∂∂σ = 0 (6.8) và nghiệm riêng của phương trình (6.9)yTyxxx∂∂+∂∂σ = - fxyyxTxy∂∂+∂∂σ = - fy (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 ; Txy = -Px khi fx = 0 ; fy = P = hằng số.* σx = 22ax+ bx ; σy = Txy = 0 khi fx = ax + b ; fy = 0* σx = 0 ; σy = -a 63y ; Txy = 22axy khi fx = axy , fy = 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=∂∂+∂∂=∂∂+∂∂yyxTxyyTyxxxσσ Đ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ệ :xqyp∂∂=∂∂.- Phương trình thứ (1) của hệ (6.8) ⇔ yTyxxx∂∂−=∂∂σ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 = yA∂∂ ; Tyx = -xA∂∂(a)Tương tự, phương trình thứ 2 : xTxyyy∂∂−=∂∂σ⇒ (σ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 = xB ; Txy = -yB(b)So sỏnh (a) v (b) ta cú : xA = yB(c) (A.dy + B.dx) l vi phõn ton phn ca 1 hm (x,y) no ú : Ta cú quan h : A = y ; B =x(d)Thay (d) vo (a) v (b) ta cú: x = 22y ; y = 22x ; Txy = - yx2(6.10)Hm (x,y) : Gi l lm ng sut Airy, l hm gii bi toỏn phng theo ng sut.III. Phng trỡnh hm ng sut Airy :- Trong chng 5 ta cú h phng trỡnh (5.5) Beltrmi l h phng trỡnh gii bi toỏn n hi theo ng sut ó tng hp cỏc iu kin v mt tnh hc, hỡnh hc, v vt lý ca mụi trng. S dng (5.5) tớnh cho biu thc ng sut phng.(1 + à)2x + 22xS= 0+ (1 + à)2y + 22yS= 0(1 + à)2z + 22zS= 0(1+à)2S +2S = 0 2S = 0 Vi S = x+ y+ z.Vỡ trong bi toỏn ng sut phng z=0 nờn S= x + y Trong bi toỏn bin dng phng : S= x + y + z = x + y +à(x + y) =(1+à)(x + y).Nờn trong bi toỏn n hi phng ta u cú :2S = 2(x + y) = 0 (6.11)(6.11) : Phng trỡnh LờVy.Thay cỏc ng sut bi hm thay (6.10) vo (6.11) ta cú :022222222=++xyyx0yyx2x4422444=++(6.12) 2(2) = 4 = 0 (6.13)47 Phương trình (6.13) : phương trình trùng điều hòa.Hàm ϕ = ϕ(x,y) : là hàm trùng điều hòa .Kết luận :- Bài tốn đàn hồi phẳng giải theo ứng suất dẫn đến việc giải phương trình (6.12) sau đó tìm các ứng suất theo (6.10).+ Nếu fx, fy ≠ 0 ⇒ Cộng thêm các nghiệm riêng.- Theo (6.10) : Việc thêm hay bớt hàm ϕ một lượng A+ Bx+Cy thì các ứng suất khơng thay đổi.- Các hệ số tích phân được xác định theo điều kiện biên tĩnh học :∗=∂∂∂−∂∂xfmyxly 222ϕϕ∗=∂∂−∂∂∂−yfmxlyx 222ϕϕ(6.14)Nếu (6.13) đủ để xác định các hằng số tích phân thì các ứng suất theo (6.10); (6.12) & (6.14) hồn tồn khơng liên quan đến các hệ số đàn hồi của vật liệu. Những bài tốn như thế là bài tốn có liên kết bên ngồi tĩnh định.⇒ Định lý LeVy-Michell : Trong biểu thức đàn hồi phẳng tĩnh định, chịu các ngoại lực tác động trên biên thì sự phân bố ứng suất khơng phụ thuộc vào các hằng số đàn hồi và như nhau đối với tổng cả các vật liệu.+ ënh l âỉåüc sỉí dủng lm cå såí cho 1 phỉång phạp thỉûc nghiãûm cọ tãn l phỉång phạp ân häưi.6.4. ĐIỀU KIỆN BIÊN CỦA HÀM ỨNG SUẤT AIRY.Việc giải bài tốn phẳng theo ứng suất rút lại thành việc giải phường trình trùng điều hòa (6.12).Nghiệm của phương trình này là hàm ứng suất ϕ phải thỏa mãn điều kiện biên.+=+=σσ∗∗ymTxylfTxyml.xfyx (6.15)Xét trường hợp fx = fy = 0Thay (6.10) vào (6.11) ta có myxlyfx∂∂∂−∂∂=∗ϕϕ222.mxlyxfy222.∂∂−∂∂∂−=∗ϕϕ (6.16)Theo (H.6.3) ta có :l = cos(n, x) = cos(900 + α) = - sinα = - dsdy48 m = cos(n, y) = cosβ = dsdx(6.15) ⇔ dsdyyfx.22∂∂−=∗ϕ - dsdxyx.2∂∂∂ϕ= -dsdyyy.∂∂∂∂ϕ-dsdxyx.∂∂∂∂ϕ= -∂∂ydsdϕ. (6.17) dsdyyxfy.2∂∂∂+=∗+ dsdxx.22∂∂ϕ = ∂∂xdsdϕ.Lấy điểm so bất kỳ trên chu tuyến làm gốc :(6.17) ⇔( )SSxXdsfAy=−=∂∂∫∗0ϕ( )SSyYdsfBx=−=∂∂∫∗0ϕ (6.18)Trong đó :A&B : Các hệ số tùy ý, biểu diễn giá trị của đạo hàm 00,SSxy∂∂∂∂ϕϕ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 S0 (H.6.4).Tại đó ta đặt các lực : A // S0x B // S0yVà 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 S0S 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 Sn ≡ 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 xyyz oPxL2t+2t−=dsdMQ(s)=dsdϕQ(s)⇒ ϕ = M (6.21)M(s) : Moment của lực đặt trên đoạn S0S 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 = yJMZZ. theo hàm ϕ : σx = 22y∂∂ϕϕ(x,y) = ax2 + bxy + cy2 + dx3 + cx2y + fxy2 + gy3 + hx4 + ix3y + ix2y2 + kxy3 + ly4. (a)ϕ phải thỏa mãn phương trình trùng điều hòa :44x∂∂ϕ + 224yx∂∂∂ϕ + 44y∂∂ϕ = 044x∂∂ϕ= h ; 224yx∂∂∂ϕ = j ; 44y∂∂ϕ= l.→ h + 2j + l = 050⇒ ϕ là hàm đa thức bậc 4 đối với x, y → h = j =l = 0 (1)σx = 22x∂∂ϕ = 2c + 2fx + 6gy + 6kxy.σy = 22x∂∂ϕ = 2c + 6dx + 6ey + 6ixy. (b)Txy = - yx∂∂∂ϕ2 =-(b + 2ex+ 2fy + 3ix2 + 3ky2 2. Các điều kiện : Xét điều kiện biên theo ứng suất :* Biên trên (y =[ ]Lxt,0;2∀: Txy = 0 , (c)σy = 0 (d)* Biên dưới (y =-[ ]Lxt,0;2∀: Txy = 0 , (e)σy = 0 (f)Từ (c) & (e) ta có :2a +6dx +2e(2t)+6ix(2t) = 02a + 6dx - 2e2t- 6ix2t = 0 ⇒ e = i = 0 e = i=f=0 (2)Từ (d) & (f) ta có : a=d=0 (3)0023322202332222222=⇒=−++−+−=++++−ftkixtfexbtkixtfexb * Biên trái (x = 0, ∀y+−2,2tt) ta có :σx= 0 (g)pdFTxytt=∫+−22.(h)Từ (g) ⇒ c = g = 0 (5)⇒ Txy = - (-43kt2 + 3ky2) = 43kt2 - 3ky2⇒ 22222232223343)343(ttttttkyyktdykyktTxydF−+− −∫ ∫−=−=δδ51⇒ b + 342kt=0⇒ b = -243kt(4) [...]... bài tố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 m = cos(n, y) = cosβ = ds dx (6. 15) ⇔ ds dy y f x . 2 2 ∂ ∂ −= ∗ ϕ - ds dx yx . 2 ∂∂ ∂ ϕ =... 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... ϕ(x,y) : Gọi là làm ứng suất Airy, là hàm để giải bài tố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 tố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 ...→ 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... 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 . 0yyxTxy∂∂+∂∂σ = 0 (6. 8) và nghiệm riêng của phương trình (6. 9)yTyxxx∂∂+∂∂σ = - fxyyxTxy∂∂+∂∂σ = - fy (6. 9 )- Nghiệm riêng của phương trình (6. 8) tìm được không. cos(900 + α) = - sinα = - dsdy48 m = cos(n, y) = cosβ = dsdx (6. 15) ⇔ dsdyyfx.22∂∂−=∗ϕ - dsdxyx.2∂∂∂ϕ= -dsdyyy.∂∂∂∂ϕ-dsdxyx.∂∂∂∂ϕ= - ∂∂ydsdϕ.