Cách giải bài toán theo chuyển vị: Nếu lấy chuyển vị làm các hàm ẩn chính, cần thu gọn hệ phương trình trên về ba phương trình đối với ba hàm chuyển vị u, v, w.. Cách giải bài toán theo
Trang 1CHƯƠNG 5 : PHƯƠNG PHÁP CƠ BẢN GIẢI BÀI TOÁN
LÝ THUYẾT ĐÀN HỒI TUYẾN TÍNH
§5.1 CÁC PHƯƠNG TRÌNH CƠ BẢN-
CÁC CÁCH GIẢI BÀI TOÁN ĐÀN HỒI TUYẾN TÍNH
5.1.1 Các phương trình cơ bản :
Trong ba chương trên ta đã lần lượt xác định ba mặt tĩnh học, hình học
và vật lý của môi trường đàn hồi tuyến tính và đưa ra 15 hàm ẩn gồm :
- Sáu thành phần ứng suất : x, y, z, Txy, Tyz, Tzx
- Ba thành phần chuyển vị : u, v, w
- Sáu thành phần biến dạng : x, y, z, xy, yz, zx
Để xác định mười lăm hàm ẩn này ta có các phương trình sau :
1 Về mặt tĩnh học :
a Hệ phương trình cân bằng Navier-Cauchy: Hệ (2.1)
) 1 (
) (
0
; ) (
0
; ) (
0
2 2 2 2 2 2
t w fz
z z y
Tyz x
Txz
t v fy
z Tzy y
y x
Txy
t u fx
z Tzx y
Tyx x
x
b Các phương trình điều kiện biên theo ứng suất: Hệ (2.3)
2 Về mặt hình học :
a Hệ phương trình biến dạng Cauchy-Navier : Hệ (3.1)
x w z
u
; z w
) 2 (
; z v y w
; y
; y x
; x
zx z
yz y
xy x
b Các phương trình liên tục của biến dạng : Hệ (3.12) và (3.13)
3.Về mặt vật lý :
a Biểu thức biến dạng biểu diễn qua ứng suất :
1
z y x
E
E
Txy G
) 1 ( 2
y = 1 y ( x z )
E ; yz = Tyz
E
Tyz G
) 1 ( 2
1
(3a)
z= 1 z ( x y )
E ; zx = Tzx
E
Tzx G
) 1 ( 2
1
b Biểu thức ứng suất biểu diễn qua biến dạng :
x = + 2Gx ; Txy = Gxy ;
y = + 2Gy ; Tyz = Gyz ;
Trang 2z = + 2Gz ; Tzx = Gzx.
5.1.2 Các cách giải bài toán đàn hồi tuyến tính :
* Về nguyên tắc 15 phương trình (1); (2) và (3a) hoặc (3b) hoàn toàn cho phép xác định được 15 hàm ẩn Để giải 15 phương trình đó ta cần thu gọn chúng về một số phương trình tương ứng với một số hàm ẩn chính Những phương trình thu gọn này là những phương trình để giải của bài toán Những ẩn số còn lại sẽ tìm được sau khi biết các ẩn số chính
1 Cách giải bài toán theo chuyển vị: Nếu lấy chuyển vị làm các hàm
ẩn chính, cần thu gọn hệ phương trình trên về ba phương trình đối với ba hàm chuyển vị u, v, w
2 Cách giải bài toán theo ứng suất: Nếu lấy ứng suất làm các hàm ẩn chính, cần thu gọn hệ trên thành sáu phương trình đối với sáu ẩn ứng suất
3 Cách giải hỗn hợp: Ngoài hai cách giải trên, trong một số bài toán,
ta sử dụng cách giải hỗn hợp, dùng một phần các hàm ẩn chính là chuyển vị
và một phần các hàm ẩn chính là ứng suất
§5.2 CÁCH GIẢI BÀI TOÁN LÝ THUYẾT ĐH THEO CHUYỂN VỊ
Chọn u, v, w là hàm ẩn cơ bản :
5.2.1.Về mặt vật lý:
Từ định luật Hooke tổng quát : x = + 2Gx
Txy = Gxy (a)
Tzx = Gzx
5.2.2 Về mặt hình học:
Từ phương trình quan hệ hình học Cauchy :
x = xu ;
yx = xv yu
zx = wx uz
;
Thay (b) vào (a) ta có : x = + G
x
u
+ G
x
u
Tyx = G
y
u x v
(c)
Trang 3Tzx = G
z
u x w
3.Về mặt tĩnh học:
Từ phương trình cân bằng tĩnh học Navier-Cauchy :
; ) t
u ( 0 fx z
Tzx y
Tyx x
x
2 2
(d) Thay (c) vào (d) ta có:
2 2
2 2
2
2 2
2 2 2
2
t
u 0
fx z
u G z x
w G y
u G y x
v G x
u G x
u
G
x
(*) t
u 0
fx z
w y
v x
u x G u z y x
G
2 2
2 2 2 2
2
2 2 2 2 2
z y
: Toán tử vi phân Laplace
z
w y
v x
u
=x+y+z = : Biến dạng thể tích tương đối (*) ( + G)
x
+ G2u + fx = 0
2 2
t
u
;
Tương tự ( + G)y + G2v + fy = 0
2 2
t
v
( + G)
z
+ G2w + fz = 0
2 2
t
w
;
Hệ (5.1): Hệ phương trình LaMê : Khi thiết lập (5.1) xuất phát từ điều kiện cân bằng và quan hệ giữa ứng suất và biến dạng nên hệ (5.1) vẫn chứa các hằng số LaMê và G
Phương trình LaMê tổng hợp được các yêu cầu về tĩnh học, hình học
và vật lý Giải (5.1) ta tìm được u, v, w sau đó xác định các biến dạng theo phương trình quan hệ hình học Cauchy và xác định các ứng suất theo định luật Hooke
4.Hệ quả: Từ phương trình LaMê trong bài toán tĩnh, khi các lực thể tích là
hằng số ta có các hệ quả sau:
a Hệ quả 1 : Đạo hàm các phương trình của hệ (5.1) lần lượt theo các biến x, y, z ta có :
( + G) 2
2
x
+ G2 x
u
= 0 ;
+ ( + G) 2
y
+ G2 y
v
= 0 ;
Trang 4( + G)z2
+ G2 z
w
= 0 ( + G) 2 + G2 = 0
Do tỷ lệ với hàm tổng ứng suất S nên ta cũng có :
Phát biểu hệ quả 1: Trong bài toán tĩnh, đàn hồi tuyến tính và đẳng hướng, khi các lực thể tích là hệ số thì hàm biến dạng thể tích và hàm ứng suất tổng là những hàm điều hòa
b Hệ quả 2 : Xét phương trình 1 của (5.2) :
( + G)
x
+ G2u +fx = 0 (a) Lấy đạo hàm bậc 2 của (a) lần lượt theo các biến x, y, z ta có :
( + G) 3
3
x
+ G2 2
2
x
u
= 0 ;
3
y
x
+ G2 2
2
y
u
= 0 ;
( + G) 2
3
z
x
+ G2 2
2
z
u
= 0 ( + G) x 2 + G2
2u = 0 (b) Theo hệ quả 1 ta có : 2 = 0 thay vào (b)
(b) 22u = 0
22w = 0 Phát biểu hệ quả 2: Trong bài toán tĩnh, đàn hồi tuyến tính và đẳng hướng, khi lực thể tích là hằng số thì các hàm chuyển vị là những hàm trùng điều hòa
c Ý nghĩa : Hệ quả này cho phép ta đoán nhận được sơ bộ dạng
nghiệm chuyển vị của bài toán đàn hồi Tất nhiên đây mới chỉ là điều kiện cần, điều kiện đủ là các chuyển vị phải thỏa mãn các phương trình cơ bản đã nêu trên
5.3 GIẢI BÀI TOÁN LÝ THUYẾT ĐÀN HỒI THEO ỨNG SUẤT
Chọn các ứng suất x, y, z, Txy, Tyz, Tzx làm hàm ẩn chính
I Trường hợp các lực thể tích là hằng số:
Trang 51 Về mặt vật lý : Dựa vào định luật Hooke
y = 1 y ( x z )
E (*)
Có S = x + y + z
(*) y = y S
E1 ( 1 )
Tương tự z = z S
E1 ( 1 )
(a)
yz =
G
1
Tyz =
E
) 1 (
2
Tyz
2 Về mặt hình học :Dựa vào phương trình liên tục của biến dạng :
2 2 2
2
y
z z
y
z y
yz
2
(b) Thay (a) vào (b) ta có :
(1 + ) 2
2
z
y
-
2
2
z
S
+(1 + ) 2
2
y
y
- 2 2
y
S
= 2(1 + )
z
y
Tyz
2
2
2 2
2 2
2 2
2
z
S y
S y
z z
y
z y
Tyz
2
(c)
3 Về mặt tĩnh học : Dựa vào hệ phương trình cân bằng tĩnh học
Navier- Cauchy
0
fx z
Tzx y
Tyx x
y
x
x z
Tzx y
Tyx
0
fy z
Tzy y
y x
x
Txy y
y z
Tzy
(2)
0
fz z
z y
Tyz x
x
Txz z
z y
Tyz
(3) Lấy đạo hàm bậc nhất (2) và (3) lần lượt theo y và z ta có :
y x
Txy y
y y
z
Tzy
2
2 2
z x
Txz z
z z
y
Tyz
2
2 2
z
Tx y
Txy x
z y
y y
z
2
2 2
2 2
(4) Thay (1) vào (4) ta có :
+
Trang 6(4)
fx x
x x
z
z y
y z
y
2
2 2
2 2
2
2
2 2
2 2
2 2
2
z
z y
y x
x z
y
Thay (d) vào (c) ta có :
2 2
2 2
2 2
2 2
2 2
2 2
2
z
S y
S y
z z
z
z y
y x
2 2
2 2
2 2
2 2
2 2
2 2
2 2
2 2
2
z
x z
x z
x y
x y
x y
x z
x y
x x
2 2 2 2
y
S z
S
= 0 (**) Trong đó : 2 = 2
2
2 2
2 2
z y
S = x + y + z
2 2 2 2
2 2
2 2
y
S z
S z
S y
S
- (1 + )2x + 2
2 2 2
z
S y
S
2 2 2 2
2 2
2
y
S z
S z
S y
S
- (1 + )2x + 2
2 2 2
y
S x
S
2 2 2
x
S z
S
= 0
(1 + )2x + S
x
S 2 2
2
= 0 Theo Hệ quả (1) ta có 2S = 0
(1 + )2x + 2
2
x
S
= 0 (1 + )2
y + 2
2
y
S
(1 + )2
z + 2
2
z
S
= 0 (1 + )2Txy +
y x
S
2
= 0 (1 + )2Tyz +
z y
S
2
= 0 (5.6)
Trang 7(1 + )2Tzx +
z x
S
= 0
Hệ phương trình (5.5) và (5.6) là 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 Giải (5.5) và (5.6) có được các ứng suất sau đó tìm các biến dạng theo định luật Hooke và tìm các chuyển vị theo hệ phương trình biến dạng Cauchy
Hệ (5.5) và (5.6) gọi là hệ phương trình Beltrmi
II Khi lực thể tích không phải là hằng số: ta cũng nhận được các phương
trình tương tự nhưng có vế phải khác 0 :
2x + fxx fyy fzz xx
x
S
1 )
1 (
1
2 2
2y + fxx fyy fzz yy
y
S
1 )
1 (
1
2 2
(5.7)
2z + fxx fyy fzz zz
z
S
1 )
1 (
1
2 2
(5.7) : Phương trình Beltrami-Michell
Từ phương trình (5.5) Beltrmi, ta cũng suy ra được 1 hệ quả về tính chất của các n0 ứng suất
Xét phương trình (1) của hệ phương trình (5.5) :
(1 + ) 2
x + 22
x
S
= 0 (1) Lấy đạo hàm bậc 2 phương trình (1) lần lượt theo x,y,z ta có :
(1 + )2 2
2
x
x
+ 4
4
x
S
= 0
+ (1 + )2 2
2
y
x
+ 2 2
4
y x
S
= 0
(1 + )2 2
2
z
x
+ 2 2
4
z x
S
= 0 (1 + ) 22x + 2
2
x
S
2S = 0 Theo hệ quả 1 2S = 0
Ta có : 22x = 0
Tương tự ta có : 4ij = 0
Trang 8ij gồm có (x, y, z, Txy, Tyz, Tzx).
Ứng suất là những hàm điều hòa kép (trùng điều hòa, bi điều hòa)
Vì ứng suất tỉ lệ với biến dạng nên biến dạng cũng là những hàm điều hoà kép
Phát biểu : Các nghiệm ứng suất , chuyển vị, biến dạng của bài toán đàn hồi tuyến tính khi lực thể tích là hằng số đều là những hàm điều hòa kép: 4ij = 0 ; 4ui = 0 ; 4
5.4 CÁC PHƯƠNG PHÁP GIẢI
1 Phương pháp thuận : là phương pháp trực tiếp tính tích phân các
phương trình Lamê (5.1) khi giải theo chuyển vị hay phương trình Beltrami (5.5) và (5.6) hay Beltrami Michell (5.7) khi giải theo ứng suất với các điều kiện biên xác định Phương pháp này rõ ràng, minh bạch vê mặt toán học nhưng phức tạp khi thực hiện
2.Phương pháp ngược : Theo phương pháp này ta cho trước chuyển
vị hay ứng suất thỏa mãn các phương trình cơ bản, rồi bằng các điều kiện biên (2.3) tìm các ngoại lực tương ứng với các chuyển vị hay ứng suất cho trước Phương pháp này để tìm được nghiệm đúng thì phải thử nhiều hàm chọn, rất cồng kềnh và có khi không thực hiện được
3 Phương pháp nửa ngược Saint - Venant : Theo phương pháp này
ta cho trước một phần các ngoại lực và một phần các chuyển vị, tìm các yếu
tố còn lại từ các điều kiện biên, chúng phải thỏa mãn các phương trình cân bằng Phương pháp này mềm dẻo, khắc phục được những khó khăn mang tính toán học của phương pháp thuận và sự cồng kềnh của phương pháp ngược
4 Nguyên lý Saint-Venant :
Nhiều bài toán của lý thuyết đàn hồi khi giải hoàn toàn thỏa mãn điều kiện biên thường gặp nhiều khó khăn, đặc biệt cách giải bài toán về thanh, tấm, vỏ Khi giải ta có thể sử dụng nguyên lý Saint-Venant đó là nguyên lý
về hiệu ứng cần bằng cục bộ của ngoại lực.theo nguyên lý này, nếu trên 1 phần nhỏ nào đó của vật thể có tác dụng của 1 hệ lực cân bằng thì ứng suất phát sinh sẽ tắt dần khá nhanh ở những đểm xa miền đặt lực
Ví dụ : Khi dùng kìm để cắt 01 sợi dây thép, ta thấy trên sợi dây tại chổ cắt tác dụng 1 hệ lực cân bằng
Dựa vào qui luật đối với vật rắn tuyệt đối, nguyên lý cục bộ có thể phát biểu theo cách khác nhau: “Tại những điểm của vật rắn cách xa điểm
Trang 9đặt lực thì trạng thái ứng suất, biến dạng của vật phụ thuộc rất ít vào cách tác dụng của lực”
Ví dụ :
F : Diện tích mặt cắt ngang
5.5 ĐỊNH LÝ DUY NHẤT NGHIỆM CỦA BÀI TOÁN LÝ THUYẾT
ĐÀN HỒI
Một vấn đề đặt ra là nghiệm của bài toán lý thuyết đàn hồi giải theo chuyển vị hay ứng suất có duy nhất không Có nghĩa ứng với tải trọng hay chuyển vị đã cho ta chỉ nhận được một hệ ứng suất hay chuyển duy nhất hay ta nhận được vài hệ nghiệm khác nhau với cùng điều kiện đã cho
* Nếu nhận được một vài hệ nghiệm thì nghiệm của bài toán lý thuyết đàn hồi đã cho là đa trị
* Định lý duy nhất về nghiệm : Nếu thừa nhận về trạng thái tự nhiên của vật và đinh luật độc lập tác dụng của lực thì nghiệm bài toán lý thuyết đàn hồi là duy nhất
Thực vậy xét bài toán cơ bản thứ nhất của lý thuyết đàn hồi Dưới tác dụng của lực bề mặt
x
f ,
y
f , z
f Lực thể tích fx, fy, fz đã cho Giả thiết ta nhận được 2 hệ nghiệm ứng suất khác nhau
x, y, z, Txy, Tyz, Tzx
x, y, z, Txy, Tyz, Tzx
Cả hai hệ ứng suất này đều phải thỏa mãn điều kiện cân bằng tĩnh học của Cauchy và điều kiện biên tĩnh học
0
x zx yx
z
T y
T x
0
*
*
*
*
x zx yx
z
T y
T x
(a)
x
f = x.l + Tyx.m + Tzx.n
f = l + T m + T n (b)
Trang 10Tương tự viết cho các phương trình còn lại Trừ các phương trình tương ứng cho nhau, ta nhận được hệ phương trình và điều kiện mới Ví dụ viết cho phương trình thứ nhất ta có :
y x
x x
( )
(Txy – Tyx) + z(Tzx - Tzx)= 0 (x - x).l + (Tyx - Tyx).m + (Tzx - Tzx).n = 0 (c)
Theo nguyên lý cộng tác dụng ta có thể xem các ứng suất trong hệ phương trìnhh (c) là một hệ ứng suất mới khi không có lực thể tích và lực bề mặt Theo giả thiết về trạng thái tự nhiên của vật liệu, các ứng suất này phải bằng 0 Do đó :
x - x = 0 ; y - y = 0 ; Tyx - Tyx = 0;
Hay x = x ; y = y ; Tyx = Tyx
Có nghĩa 2 hệ ứng suất này trùng nhau Đó là điều cần chứng minh!