BÀITOÁNỔNĐỊNHĐÀNHỒICỦAKẾTCẤUVỎTHOẢIMẶTBẰNG
HÌNH CHỮNHẬTVỚICÁCLIÊNKẾTBIÊNKHÁCNHAU
PGS. TS. LÊ NGỌC THẠCH, ThS. MAI CHÂU ANH
Trường Đại học Xây dựng
1. Mở đầu
Khi giải bàitoánổn định, ta phải thiết lập các phương trình cân bằngcủa nó, phương trình này
có dạng giống như phương trình cân bằngcủabàitoán bền. Song những điều kiện cân bằng này
chưa nói lên được dạng cân bằng đó ổnđịnh hay không ổn định. Điểm khácnhau giữa hai loại bài
toán này là đối vớibàitoán bền, từ các phương trình này ta sẽ tìm được các giá trị nội lực, ứng suất,
chuyển vị. Còn bàitoánổnđịnh thì từ điều kiện mấtổnđịnhcủa hệ, ta sẽ tìm được các thông số lực
tới hạn biểu thị độ an toàn về mặtổnđịnhcủakếtcấu đối với một nhóm lực nhất định.
Trên cơ sở lý thuyết tổng quát về sự cân bằng và cân bằngổnđịnhcủa hệ đàn hồi, người ta
phân loại bàitoánổnđịnh thành hai trường hợp: loại I và loại II. Đối vớibàitoánổnđịnh loại I, từ hệ
phương trình ổnđịnh thiết lập các điều kiện để cho hệ mấtổnđịnh ta sẽ tìm ngay được tải trọng tới
hạn tương ứng. Để tìm giá trị của lực tới hạn trong bàitoánổnđịnh loại II không áp dụng được cách
giải trực tiếp như trên, ở đây ta sẽ tìm cách thiết lập trực tiếp các điều kiện cực trị của phiếm hàm
bằng các phương pháp trực tiếp: phương pháp Ritz - Timosenko, phương pháp Buovnov - Galerkin.
Khi giải bàitoánổnđịnhcủavỏthoải thì việc quan trọng đầu tiên là chọn hàm độ võng sao cho thoả
mãn các điều kiện biên. Các phương trình ổnđịnh rất phức tạp về mặttoánhọc cho nên giải những bàitoán
cụ thể rất khó. Để có thể tìm được các giá trị lực tới hạn mà không gặp nhiều khó khăn, cần thiết phải làm
giảm bớt bậc củacác phương trình vi phân. Cách tìm nghiệm được quy về hai nhóm chủ yếu: chính xác và
gần đúng. Nhóm lời giải chính xác gồm nhiều phương pháp: bàitoán tìm nghiệm Navier dưới dạng chuỗi
lượng giác kép, bàitoán tìm nghiệm Levy dưới dạng chuỗi lượng giác đơn, phương pháp biến phân của
Ritz - Timoshenko, của Buovnov - Galerkin. Nhóm lời giải gần đúng có các phương pháp: phương pháp sai
phân, phương pháp phần tử hữu hạn
1
,
2
. Trong bàibáo này, các tác giả dùng phương pháp của
Buovnov - Galerkin khảo sát ổnđịnh cho vỏthoảimặtbằnghìnhchữnhật có tỷ số cạnh dài trên cạnh ngắn
khác nhau và tìm nghiệm giải tích cho các trường hợp vỏthoảimặtbằnghìnhchữnhậtliênkết ngàm và
vừa ngàm vừa khớp.
2. Hệ phương trình cân bằng trong bàitoánổnđịnhcủavỏthoải
1
,
2
h
D
4
w =
2
2
x
2
2
y
w
+
2
2
y
2
2
x
w
- 2
yx
w
2
+ k
x
2
2
y
+ k
y
2
2
x
+
h
q
= L(
, w) + k
x
2
2
y
+ k
y
2
2
x
+
h
q
= L(
, w) +
2
k
+
h
q
(1)
E
1
4
= (
yx
w
2
)
2
-
2
2
x
w
2
2
y
w
-
k
x
2
2
y
w
- k
y
2
2
x
w
(2)
= -
2
1
L(w,w) -
w
k
2
Với D là độ cứng trụ của vỏ: D =
)1(12
2
3
Eh
3. Trường hợp vỏthoải chịu tác dụng của tải trọng pháp tuyến q
Đây là trường hợp rất phổ biến (chịu trọng lượng bản thân, áp lực của chất lỏng ). Bàitoánổn
định của
vỏ thoải chịu uốn là bàitoánổnđịnh loại hai. Khi tăng tải trọng, lúc này biến dạng củakếtcấu
phát triển nhưng không thay đổi tính chất, không phân nhánh. Giá trị của lực q tương ứng với
khi độ võng w tăng mà không cần tăng tải trọng gọi là lực tới hạn. Khi q = q
th
, sự cân bằng giữa
nội lực và ngoại lực đạt đến trạng thái tời hạn. Khi q > q
th
, sự cân bằng chỉ có thể xảy ra khi
giảm tải trọng q. Trạng thái tời hạn được xác định từ điều kiện:
dw
dq
= 0.
Xét vỏthoảimặtbằnghìnhchữnhật có kích thước chiều dài và chiều rộng là a, b. Gọi k
x
, k
y
là
độ cong chính củavỏ theo hai phương x và y. Chọn hàm độ võng w và ứng lực
thỏa mãn các điều
kiện biên.
Dùng phương trình (1), (2) và đặt:
X =
h
D
4
w - L(
, w) - k
x
2
2
y
- k
y
2
2
x
-
h
q
0
=
22
4
4
4
4
4
2
yx
w
y
w
x
w
h
D
-
2
2
x
2
2
y
w
-
2
2
y
2
2
x
w
+ 2
yx
w
2
- k
x
2
2
y
- k
y
2
2
x
-
h
q
0 (3)
Y =
E
1
4
+
2
1
L(w,w) + k
x
2
2
y
w
+ k
y
2
2
x
w
0
=
22
4
4
4
4
4
2
1
yxyxE
+
2
2
x
w
2
2
y
w
- (
yx
w
2
)
2
+ k
x
2
2
y
w
+ k
y
2
2
x
w
0 (4)
Theo A.C.Volimip, dùng phương pháp Bouvnov - Galerkin và đặt phương trình:
dxdy
b
y
a
x
X
nn
a b
sinsin
0 0
= 0 ;
dxdy
b
y
a
x
Y
nn
a b
sinsin
0 0
= 0 (5)
Ở đây lấy n = 1 cho trường hợp liênkết khớp; n = 2 cho trường hợp liênkết ngàm.
Thực hiện tích phân và đặt các thông số không đơn vị:
k
x
* =
h
ak
x
2
; k
y
* =
h
bk
y
2
; k*
= k
x
* + k
y
* ;
=
b
a
;
=
h
1
; q* =
E
q
2
2
h
ab
(6)
Lấy
= 0,3 được quan hệ giữa tải trọng và độ võng. Khi hệ bị mấtổn định, trạng
thái tới hạn được xác định từ điều kiện:
d
dq *
= 0 và tìm được giá trị
tương ứng.
Sau đó khảo sát và vẽ đồ thị quan hệ q -
khi
= 1; 1,5; 2; 2,5 vớicác giá trị k* = 0, 12, 24, 36,
48, 60 khác nhau.
3.1. Ổnđịnhcủavỏthoảiliênkết khớp chịu tác dụng của tải trọng pháp tuyến
A.C.Volimip đã chọn hàm độ võng w và ứng lực
, giải ra được q* và khảo sát đồ thị ứng với
trường hợp
=
b
a
= 1. Ở đây tác giả khảo sát thêm một số trường hợp còn lại.
w =
1
a
x
sin
b
y
sin
;
= A
1
a
x
sin
b
y
sin
q* =
2
2
6
2
2
2
22
2
3
2
2
1192
1
*
1
1
16
*
1
1
1
1
9
32
kk
(7)
=
32
3
k*
)1(2
)
1
(
*3
32
1
2
44
2
k
điều kiện: k*
)1(6
)
1
(
2
22
(8)
1
k
2
k
3
k
4
k
5
k
6
k
1
k
2
k
3
k
4
k
5
k
6
k
q*
2
4
6
8
10
500
1000
1500
q
q*
2
4
6
8
10
-
250
250
500
750
1000
1250
Hình 1.
Đồ thị quan hệ q* -
đối vớimặtbằnghìnhchữnhật có a = b và 1,5b
1
k
2
k
3
k
4
k
5
k
6
k
1
k
2
k
3
k
4
k
5
k
6
k
q*
2
4
6
8
10
200
400
600
800
1000
2
4
6
8
10
200
400
600
q*
Hình 2.
Đồ thị quan hệ q* -
đối vớimặtbằngchữnhật có a = 2b và a = 2,5b
Từ biểu đồ, ta nhận thấy ở đây dạng cân bằng bị phân nhánh, sau khi kếtcấu bị mấtổnđịnh ứng
với tải trọng tới hạn trên, nó chuyển qua trạng thái cân bằng khác, trạng thái tới hạn lúc này ứng với
tải trọng tới hạn dưới, nên để đạt đến kết quả chính xác hơn đặt hàm độ võng và ứng lực thỏa mãn
các điều kiện biên dạng:
w =
nm
mn
b
yn
a
xm
,
sinsin
;
=
nm
mn
b
yn
a
xm
A
,
sinsin
Tương tự như trên, ta đặt X , Y. Áp dụng phương pháp Bouvnov - Galerkin và đặt phương trình:
a b
dxdy
b
yj
a
xi
X
0 0
sinsin
= 0 ;
a b
dxdy
b
yj
a
xi
Y
0 0
sinsin
= 0
Thực hiện tích phân:
h
q
a
km
b
kn
AnmfA
bab
n
a
m
f
y
x
mnmnmnmn
2
2
2
24
22
22
4
2
2
2
2
26
163
2
)(
16
= 0
(9)
2
2
2
2
2
2
2
2
2
2
2
2
2
a
m
1
-
1
3
161
b
n
a
km
b
kn
f
b
na
a
mb
A
E
y
x
mnmnmn
= 0 (10)
Đặt: k
x
* =
h
ak
x
2
; k
y
* =
h
bk
y
2
; k*
= k
x
* + k
y
* ;
=
b
a
;
mn
=
h
mn
; q* =
E
q
2
2
h
ab
(11)
Từ (10)
A
1
, thay vào (9), (11), lấy
= 0,3 ta có:
q* =
2
2
6
2
2
2
2
2
23
2
2
1192
1
*
1
1
16
1
1
*
1
1
9
32
n
m
k
n
m
n
m
k
n
m
mnmn
(12)
Xét trường hợp đặc biệt, vỏ có mặtbằnghình vuông
=
b
a
= 1, m = n ta có:
q* =
mnmnmn
m
m
k
m
m
m
m
k
m
m
2
2
6
2
2
2
2
2
23
2
2
1192
1
*
1
16
1
1
*
1
1
9
32
(13)
Lấy
mn
d
dq
*
= 0 được:
mn
=
32
3
k*
)1(
8
*3
32
1
2
44
2
m
k
(14)
Xét trường hợp m = 1, từ (13): q* = 8,764
3
1
- 2,1k*
2
1
+ (0,154k*
2
+ 21,943)
1
(15)
Trường hợp m = 2, từ (13):
q* = 8,764
3
1
+ 5,609
3
2
- 2,1k*
2
1
-1,578k*
2
2
+ (0,154k*
2
+ 21,943)
1
+ (0,099k*
2
+ 4,288)
2
(16)
Đặt:
=
1
2
, đưa (16) về thành:
q* = 8,764
3
1
+ 5,609
3
3
1
- 2,1k*
2
1
- 1,578k*
2
2
1
+ (0,154k*
2
+ 21,943)
1
+ (0,099k*
2
+ 34,288)
1
(17)
Cân bằng (17), lấy đạo hàm và giải phương trình bậc 3 theo biến (
1
) có:
1
= 0 ;
1
= 0,089(1,578
1260
2
k
) (18)
Từ (18), có nghiệm thực của (
1
) là: (
1
)= 0,089(1,578
1260
2
k
) với k
35,5
Thay giá trị của
1
từ (18) vào (17) và khảo sát biểu đồ quan hệ q -
1
1
k
2
k
3
k
4
k
5
k
6
k
2
4
6
8
5 0 0
1 0 0 0
q *
Hình 3.
Đồ thị quan hệ q* -
đối vớimặtbằnghình vuông
3.2. Ổnđịnhcủavỏthoảiliênkết ngàm chịu tác dụng của tải trọng pháp tuyến
Chọn hàm độ võng w và ứng lực
thỏa mãn các điều kiện biên dạng:
w =
b
y
a
x
f
22
1
sinsin
;
=
b
y
a
x
A
22
1
sinsin
Ta có: q* =
2
2
4
2
2
2
2
2
3
2
4
18
1
*
1
1
8
3
1
8
9
*
1
1
4
3
kk
(19)
Tìm được:
=
2
2
9
k*
)1(
)
1
(192
*
75
2
1
2
4
2
4
k
điều kiện: k*
)1(
192
75
)
1
(
2
22
(20)
1
k
2
k
3
k
4
k
5
k
6
k
1
k
2
k
3
k
4
k
5
k
6
k
q *
2
4
6
8
1 0
1 2
1 0 0 0
1 00 0
2 00 0
3 00 0
4 00 0
5 00 0
q*
2
4
6
8
1 0
1 2
1 0 0 0
1 0 0 0
2 0 0 0
3 0 0 0
4 0 0 0
5 0 0 0
Hình 4.
Đồ thị quan hệ q* -
đối vớimặtbằnghìnhchữnhật có a = b và 1,5b
1
k
2
k
3
k
4
k
5
k
6
k
1
k
2
k
3
k
4
k
5
k
6
k
q *
5
1 0
1 5
-
2 0 0 0
2 0 0 0
4 0 0 0
q *
2
4
6
8
1 0
1 2
1 4
1 0 0 0
2 0 0 0
3 0 0 0
Hình 5.
Đồ thị quan hệ q* -
đối vớimặtbằngchữnhật có a = 2b và a = 2,5b
3.3. Ổnđịnhcủavỏthoảiliênkết khớp ở hai biên đối diện (song song với trục Oy), hai biên
còn lại ngàm chịu tác dụng của tải trọng pháp tuyến
Chọn hàm độ võng w và ứng lực
thỏa mãn điều kiện biên dạng:
w =
b
y
a
x
f
2
1
sinsin
;
=
b
y
a
x
A
2
1
sinsin
;
2
2
2
5
2
2
2
2
2
2
2
3
2
2
3
112
2
1
2
1
16
3
*
2
1
2
1
16
3
1
16
*
2
1
2
1
16
3
3
2
1
2
1
16
3
1
3
*
k
kq
(21)
=
4
1
k*
)1(
2
1
2
1
16
3
12
*7
12
1
2
2
2
2
2
2
2
k
điều kiện: k*
)1(
12
7
2
1
2
1
16
3
2
2
2
2
(22)
1
k
2
k
3
k
4
k
5
k
6
k
1
k
2
k
3
k
4
k
5
k
6
k
q*
2
4
6
8
10
1000
1000
2000
q *
2
4
6
8
1 0
- 5 0 0
5 0 0
1 0 0 0
1 5 0 0
Hình 6.
Đồ thị quan hệ q* -
đối vớimặtbằnghìnhchữnhật có a = b và 1,5b
1
k
2
k
3
k
4
k
5
k
6
k
1
k
2
k
3
k
4
k
5
k
6
k
q*
2
4
6
8
10
-200
200
400
600
800
1000
1200
2
4
6
8
1 0
5 0 0
1 0 0 0
1 5 0 0
2 0 0 0
q *
Hình 7.
Đồ thị quan hệ q* -
đối vớimặtbằngchữnhật có a = 2b và 2,5b
4. Nhận xét và kết luận
Từ cáckết quả và biểu đồ tính giá trị lực tới hạn q
th
, q* cho các trường hợp vỏ có k* và
khác
nhau, ta nhận thấy:
-
Với cùng một giá trị
: khi k* tăng, q*
th
cũng tăng theo.
-
Khi
tăng: giá trị k* cần để q* có cực trị tăng, q*
th
giảm hay nói cách khác: đối vớivỏ dài, độ cứng
giảm và giá trị của lực tới hạn q*
th
cũng giảm.
-
Vỏthoảiliênkết ngàm và vừa liênkết ngàm vừa liênkết khớp cho giá trị lực tới hạn q
th
* lớn hơn
so vớivỏliênkết khớp khi chịu uốn do có độ cứng lớn hơn.
-
Đối vớivỏ có độ cong bé, vỏthoảiliênkết ngàm có độ ổnđịnh rất cao.
-
Giá trị của lực tới hạn không chỉ phụ thuộc vào các tham số hình học, vật lý, điều kiện biên mà
còn phụ thuộc vào các yếu tố khác như điều kiện ban đầu.
-
Bằng cách áp dụng phương pháp Buovnov - Galerkin và cách xử lý như trong mục
3
ta có thể
đưa bàitoán phi tuyến về một hệ mà có thể cho ngay các nghiệm cần tìm là đơn trị.
-
Trong quá trình tính toán, nếu biết phối hợp vớicác phần mềm tin học ứng dụng trong lĩnh vực
toán học (Mathematica, Matlab ) thì vẫn bằngcác phương pháp cổ điển (phương pháp biến
phân của Bouvnov - Galerkin, Ritz - Timosenko ) ta có thể tìm được nghiệm củabàitoán phức
tạp dưới dạng giải tích, từ đó dễ dàng khảo sát và rút ra những kết luận bổ ích tạo điều kiện giải
các bàitoánổnđịnh cả trong và ngoài miền đànhồibằngcác phương pháp số như : phương
pháp sai phân hữu hạn, phần tử biên, phần tử hữu hạn,
-
Bàibáo đã nêu ra được phương pháp xác định biểu thức của lực tới hạn trong trường hợp có kể
đến các yếu tố phi tuyến hình học. Từ đây, có thể định hướng cho bàitoán phi tuyến vật lý.
TÀI LIỆU THAM KHẢO
1. Xtiphen P. Timosenko, Jem M. Gere. Ổnđịnhđàn hồi. (bản dịch tiếng Việt). NXB Khoahọc và Kỹ
thuật.Hà Nội, 1976.
2. A.C.Volimip. Ổnđịnhcủa hệ đàn hồi. Nhà xuất bản Hayka. Moscow, 1967.
3. Phillip L.Gould. Analysis of Shells and Plates - Springer – Verlag. New York Inc, 1988.
4. LÊ NGỌC HỒNG. Cơ sở lý thuyết tấm và vỏ mỏng đàn hồi. Tập bài giảng cho học viên cao học.
Trường Đại học Xây dựng. Hà Nội, 2001.
. BÀI TOÁN ỔN ĐỊNH ĐÀN HỒI CỦA KẾT CẤU VỎ THOẢI MẶT BẰNG
HÌNH CHỮ NHẬT VỚI CÁC LIÊN KẾT BIÊN KHÁC NHAU
PGS. TS. LÊ NGỌC THẠCH,. cho các trường hợp vỏ thoải mặt bằng hình chữ nhật liên kết ngàm và
vừa ngàm vừa khớp.
2. Hệ phương trình cân bằng trong bài toán ổn định của vỏ thoải