KHẢO SÁT TRẠNG THÁI NHIỆT CỦA MẶT ĐƯỜNG BÊ TÔNG XI MĂNG BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN
ThS. NGUYỄN ĐĂNG KHOÁT PGS. TS. TRỊNH VĂN QUANG Bộ môn Kỹ thuật Nhiệt Khoa Cơ khí Trường Đại học Giao thông Vận tải
Tóm tắt: Bài báo trình bày cách sử dụng phương pháp phần tử hữu hạn để đánh giá trạng thái nhiệt của mặt đường bê tông xi măng dưới tác động của thay đổi nhiệt độ không khí và bức xạ măt trời.
Summary: The paper presents the way to use the finite element method to study thermal state of the concrete surface of the road under impact of the varying air temperature and solar radiation.
CT 2
I. ĐẶT VẤN ĐỀ
Bài toán trạng thái nhiệt mặt đường đến nay đã có nhiều tác giả nghiên cứu [1], [2], ... Với mục đích mở rộng các phương pháp tính nhiệt, bài viết trình bày cách sử dụng phương pháp (phương pháp) phần tử hữu hạn (PTHH) trong tính nhiệt, để xác định nhiệt độ và đưa ra các nhận định về trạng thái nhiệt của tấm bê tông (BT) dưới tác động của các yếu tố khí hậu thay đổi.
II. PHƯƠNG PHÁP KHẢO SÁT
2.1. Giới hạn bài toán, dữ liệu ban đầu
C)
5 26,3
6 26,5
7 27,2
8 27,7
9 28,5
10 29,4
11 30,1
12 30,7
13 31,3
14 31,8
15 32,0
+ Khảo sát tấm BT dày L = 0,3m, dài và rộng 7,5m đặt trên nền đất; BT có hệ số dẫn nhiệt k = 1,265W/m0C; khối lượng riêng ρ = 2200kg/m3; nhiệt dung riêng c = 1215J/kg0C; mặt trên toả nhiệt với không khí với h = 7,89W/m20C, và hấp thụ tia mặt trời với ε = 0,65. Nền đất có: kN = 0,52W/m0C; ρN = 2050kg/m3; cN = 1840J/kgđộ. Ở độ sâu đủ lớn nền đất có nhiệt độ không đổi là 28,8oC. Nhiệt độ không khí TK, bức xạ mặt trời E trong ngày tháng 6 thay đổi theo số liệu của ngành khí tượng, bảng 1, tốc độ gió trung bình w = 2,4 m/s.
Bảng 1 16 31,7
Giờ ( T 0 K
0
34,89 209,3
407,0 610,5 779,2
895,5
930,4
872,2
744,3
593,1 401,2
E(W/m2)
18 17 31,3 30,2 203,5 58,15
19 29,6 0
20 28,8 0
21 28,4 0
22 28,2 0
23 27,6 0
24 27,2 0
1 27,0 0
2 26,8 0
3 26,5 0
4 26,4 0
Giờ ( T 0 C) K E(W/m2)
Truyền nhiệt qua tấm BT có bề dày nhỏ hơn rất nhiều so với bề rộng và dài, được mô tả bởi
2
ρ.c
k
=
T ∂ τ ∂
T ∂ 2 x ∂
phương trình vi phân dẫn nhiệt một chiều: (1); với x là bề dày tính từ mặt trên.
k
= +
−
−
( m q h T
T K
T ∂ x ∂
) (2). Với T, Tm và TK tương ứng là nhiệt độ trong tấm bê tông, nhiệt độ bề mặt và không khí (oC), τ là thời gian (s), x là chiều sâu kể từ mặt tấm (m), q dòng bức xạ mặt trời (W/m2).
Điều kiện biên toả nhiệt và bức xạ tại mặt trên x=0, là:
2.2. Phương pháp phần tử hữu hạn
1. Rời rạc miền nghiệm
Bề dày tấm bê tông được rời rạc thành 12 phần tử (PT) ký hiệu 1,2,3,... mỗi PT dài
là l = 0,3m/12 = 0,025m, và 13 nút ký hiệu 1,2,3,…,13. Nền đất chọn một PT thứ 13, dài l và nút 14, hình 1.
Hình 1. Sơ đồ rời rạc lớp bê tông thành các PTHH
2. Chọn hàm nội suy
=
=
−
]
[ NNN i
j
x l
⎛ ⎜ ⎝
⎤ ⎥ ⎦
(3); với l là chiều dài của một PT, x là Chọn hàm bậc nhất:
⎡ ⎛ 1 ⎜ ⎢ ⎝ ⎣ =
⎞ ⎟ ⎠ +
=
]{ }TN [
x ⎞ ⎟ l ⎠ TNTNT i
i
j
j
toạ độ trong phần tử, nên nhiệt độ là (4); Ti và Tj là nhiệt độ tại hai
N ∂
N ∂
j
j
CT 2
nút của PT. Đạo hàm của hàm nội suy [B], gradient nhiệt độ [g] là
=
=
=
[ ]{ } [ ]g TB =
[B11 ] =
] (5);
T i
T j
N ∂ x ∂
N ∂ i x ∂
x ∂
N ∂ i x ∂
x ∂
T ∂ x ∂
⎡ ⎢ ⎣
⎤ 1 [ −=⎥ l ⎦
⎡ ⎢ ⎣
⎤ 1 [ −=⎥ l ⎦
⎫ ⎬ ⎭
T ⎧ ] 1 11 ⎨ T ⎩ j
(6)
3. Thiết lập phương trình đặc trưng của phần tử
2
Có nhiều phương pháp thiết lập Phương trình ma trận đặc trưng của phần tử. Trong tính nhiệt, phương pháp Biến phân (Rayleigh Ritz) và Galerkin là hai phương pháp quan trọng nhất vì cho kết quả chính xác như nhau. Ở đây chọn phương pháp Galerkin, lấy hàm trọng số là hàm nội suy Ni .
0
ρc
=
−
V
∫
T ∂ 2 x ∂
T ∂ τ ∂
⎞ ⎟⎟ dV ⎠
⎛ ⎜⎜ kN i ⎝
(7) Phương pháp Galerkin yêu cầu (1) thoả:
dS
k
.
dx
dydz
=
=
−
kN i
kN i
∫
∫
∫
∫
∂ x ∂
T ∂ x ∂
T ∂ x ∂
T ∂ x ∂
T ∂ x ∂
N ∂ i x ∂
⎛ ⎜ ⎝
⎞ dV ⎟ ⎠
⎛ kdN ⎜ i ⎝
⎞ dydz ⎟ ⎠
V
V
S
V
Lấy tích phân từng phần số hạng đầu của (7) có
(8)
S là diện tích mặt ngoài có ở thành: đối lưu, bức xạ. Điều kiện biên (2) tr
dS
−=
−
−
)dS
i
qdSN i S
kN S
( TThN K S
∫
∫
∫
T ∂ i ∂ x
(9)
N ∂
j
i
Thay (4), (8) và (9) vào (7) sắp xếp lại sẽ được:
k
+
+
−=
+
( ρ.cN
(
{
i
) dVN j
i
j
qdSN i
]
V
V
S
S
dSThN i K S
[ ∫
∫
} ) ( ) ∫ τTdSNhN j
∫
∫
N ∂ x ∂
x ∂
( ) τT ∂ j τ ∂
⎫ ⎬ ⎭
⎛ ⎜⎜ ⎝
⎞ ⎟⎟ dV ⎠
⎧ ⎨ ⎩
⎡ ⎢ ⎣
⎤ ⎥ ⎦
(10)
+
[ TK
]{ } { }f =
T ∂ τ ∂
⎧ :[ ] C ⎨ ⎩
⎫ ⎬ ⎭
V iết (10) gọn dạng ma trận là
(11)
T
T
+
trận đặc trưng của phần tử
] [
] [ [ k B B dV
]
]
( T NNρ.c
) ] dV (12); [K] là ma trận độ cứn
(1 dung,[ ] C 1) là Phương trình ma [ g, [
]
(
(
)
V
[ ∫=
V
S
∫
∫
T
. Trong đó: [C] gọi là ma trận nhiệ t ) h N N dS⎤ K ⎡ [ ] [ ] = ⎥ ⎢ ⎦ ⎣
qdS
−=
+
[ N
]
] [ T dSTNh
K
S
S
∫
∫
(13) (14). ; {f} là véc tơ tải nhiệt,{ }f
4. Rời rạc theo thời gian
Có thể rời rạc (11) bằng phương pháp Sai phân hữu h ạn (SPHH) hay PTHH.
p
P
1 +
T
T
a. Rời rạc theo thời gian bằng HH phương pháp SP
≈
Δ
− τ
T ∂ τ ∂
⎧ ⎨ ⎩
⎫ ⎬ ⎭
⎧ ⎨ ⎩
⎫ ⎬ ⎭
… (11) trở thành: Theo phương pháp SPHH thì
p
1p +
1p +
+
=
+
(15)
(1 −+
{
[ ] C
[ ] }{ } TKξ.Δτ
[ ] { (1C −
[ ] }{ } τ TKξ)Δ−
( { } fξΔτ
)p { } fξ)
CT 2
Trong đó Δτ là bước thời gian ứng với hai thời điểm liên tiếp p và (p+1); ξ= 0÷1 là tham
số tùy chọn: ξ= 0÷1. Nếu lấy ξ=1 không cần điều kiện hạn chế chọn bước thời gian Δτ.
N
[ N
p N
]1p
+
=τ
Δτ
(16), phương pháp Galerkin yêu cầu (11) thoả: b. Rời rạc bằng phương pháp PTHH: Có các phương pháp số dư trọng số khác nhau để rời rạc (11) theo thời gian. Dùng Galerkin, hàm trọng số là hàm nội suy bậc nhất của thời gian Nτ:
0
τd
=
+
[ TK
0
∫
T ∂ τ ∂
⎧ ⎧ [ ] CN ⎨ ⎨ p ⎩ ⎩
⎫ ]{ } { } f − ⎬ ⎭
p
p
p
]
(17)
+
=
T 1p +
T 1p +
f 1p +
12 21
− −
⎫ ⎬ ⎭ Thay (4), (5) và (16) vào (17), biến đổi được: [ ] C Δτ
T
T
⎡ ⎢ ⎣
⎡ ⎢ ⎣
⎤ ⎥ ⎦
⎧ ⎨ ⎩
⎫ ⎬ ⎭
⎧ ⎨ ⎩
⎧ ⎨ f ⎩
⎫ ⎬ ⎭
(18)
P
p
1p +
1p +
⎫ ⎬ ⎭ (19)
f
+
=
−
+
+
{
Dẫn tới [ ] C
[ ] } TKΔτ
[ ] { C
[ ] } TKΔτ
11 ⎤ ⎥ 11 ⎦ ( f
[ K2 3 )Δτ
1 2
(19) yêu cầu Δτ phải chọn đủ nhỏ để nghiệm hội tụ. Có thể tìm nghiệm Tp+1 từ (15) hoặc (19)
5. Tính các số hạng trong phương trình ma trận đặc trưng của phần tử một chiều
a. Ma trận nhiệt dung phần tử [C] . Thay các số liệu của bê tông và nền đất tính được
=
] 121C
=−
] [ 13C
11138 22276
15717 31434
22276 ⎡ ⎢ 11138 ⎣
⎤ ⎥ ⎦
31434 ⎡ ⎢ 15717 ⎣
⎤ ⎥ ⎦
(20); - Nền đất: (21) - Bê tông :[
b. Ma trận độ cứng phần tử [K]: - Các PT 1, 2÷12 và 13 tương ứng là
1
50,6
50,6
1
58,49
=
] K 12 2
=−
1
− 1
50,6
− 50,6
50,6
−
−
Ak l
⎡ Ak ⎢ − l ⎣
1 − ⎤ +⎥ 11 ⎦
01 ⎡ ⎤ hA =⎥ ⎢ 00 ⎣ ⎦
⎡ ⎢ ⎣
1 ⎤ =⎥ ⎦
⎡ ⎢ − ⎣
50,6 − ⎤ ⎥ 50,6 ⎦
⎡ ⎢ ⎣
⎤ ⎥ ⎦
(23) (22); [ [ ] K 1
=
1 1
20,8 - 20,8
−
Ak l
⎡ ⎢ ⎣
1 − ⎤ =⎥ 1 ⎦
⎡ ⎢ ⎣
- 20,8 ⎤ ⎥ 20,8 ⎦
K
(24). Với A là diện tích truyền nhiệt, A = 1m2 [ ] K 13
−=
+
=
{ } f
=
{ } f
1
13-2
1 ⎧ ⎫ qA ⎬ ⎨ 0 ⎩ ⎭
1 ⎧ ⎫ AhT ⎬ ⎨ K 0 ⎩ ⎭
qAAhT + ⎧ ⎫ ⎨ ⎬ 0 ⎩ ⎭
⎫ ⎬ ⎭
0 ⎧ ⎨ 0 ⎩
c. Véc tơ phụ tải nhiệt {f}: (25);
(26)
6. Lắp ghép các phần tử
22276
11138
0
0
0
0
0
0
0
0
0
0
0
0
22276 11138
11138 44552
0 11138
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0 0
0 0 0
11138 0 0
44552 11138 0
11138 44552 11138
0 11138 44552
0 0 11138
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
=
[ C
]
0 0
0 0
0 0
0 0
0 0
11138 0
44552 11138
11138 44552
0 11138
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
11138 0
44552 11138
0 0
0 0
0 0
11138 44552
0 11138
0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
11138 0 0
44552 11138 0
11138 44552 11138
0 11138 63710
0 0 15717
0 0 0
0
0
0
0
0
0
0
0
0
0
0
15717
41434
0
⎤ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦
⎡ ⎢ 11138 ⎢ ⎢ 0 ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣
- Ma trận nhiệt dung toàn hệ:
(27)
CT 2
- Ma trận độ cứng toàn hệ [K], véc tơ phụ tải tổng {f}, sau khi áp đặt điều kiện biên tại lớp
nền có TN =T14= 28,8 sẽ được các ma trận (28) và (29) dưới đây.
7. Giải hệ phương trình
p
1
1P +
+
−
fΔτ
+
+
=
{ }
(19). Từ (15) suy ra: { } T
[ KΔτ
[ ] { C
] }
)1p
Vì số liệu cho theo giờ (bảng 1), tức Δτ = 3600s, nên chọn (15) để giải sẽ thuận tiện hơn ( [ ]{ } (30) TC*
Lấy {T}p=0 =28,8oC, thay (25), (26) và (27) vào (30), lập trình và giải (30) trên Matlab qua
q
+
-50.6 101.2
0 -50.6
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
7,89.TK 0 0
-50.6 0
101.2 -50.6
-50.6 101.2
0 -50.6
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0
0 0
-50.6 0
101.2 -50.6
-50.6 101.2
0 -50.6
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0
=
[ K
]
=
{ } f
0 0
0 0
0 0
0 0
101.2 -50.6
-50.6 101.2
0 -50.6
0 0
0 0
0 0
0 0
-50.6 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
-50.6 0
101.2 -50.6
-50.6 101.2
0 -50.6
0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
-50.6 0
101.2 -50.6
-50.6 101.2
0 -50.6
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
-50.6 0
71.4 0
0 0
58.49 ⎡ ⎢ - 50.6 ⎢ ⎢ 0 ⎢ 0 ⎢ ⎢ 0 ⎢ 0 ⎢ ⎢ 0 ⎢ 0 ⎢ ⎢ 0 ⎢ ⎢ 0 ⎢ 0 ⎢ ⎢ 0 ⎢ 0 ⎢ ⎢ 0 ⎣
28,8
192 thời điểm. Kết quả được lập thành bảng và đồ thị.
⎡ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ 599,04 ⎢ ⎢ ⎣
⎤ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦
( (29) (28)
III. KẾT QUẢ TÍNH TOÁN VÀ CÁC NHẬN XÉT
3.1. Diễn biến và thay đổi nhiệt độ tại 14 nút sau 192 thời điểm
Diễn biến và thay đổi nhiệt độ tại 14 nút sau 192 giờ thể hiện trên hình 2 và 3. Từ đó rút ra
các nhận xét:
1. Thay đổi nhiệt độ theo thời gian tại 14 nút là 14 đường dao động, gồm 8 chu kỳ tương ứng 8 ngày đêm. Từ chu kỳ 5 trở đi, dao động trong các chu kỳ tuân theo cùng một quy luật, nhiệt độ tại mỗi nút ở cùng thời điểm tương ứng trong ngày đã hội tụ tới giá trị ổn định.
2. Kể từ chu kỳ thứ năm, có thể chọn nhiệt độ tại 24 thời điểm liên tiếp để đại diện cho thay đổi nhiệt độ trong tấm bê tông qua một ngày đêm điển hình mùa hè. Ở đây chọn từ thời điểm 154 đến 177.
3.2. Đặc tính thay đổi nhiệt độ trong tấm bê tông qua một ngày đêm điển hình tháng 6
Đặc tính thay đổi nhiệt độ được biểu thị trên hình 4 và 5. Có thể rút ra nhận xét sau:
1. Thay đổi nhiệt độ tại mặt trên cùng là hàm chu kỳ, không biểu thị là hàm số cosin của
thời gian τ, nhưng càng vào sâu trong tấm BT, dạng dao động nhiệt độ càng tiến tới hình sin.
2. Thời điểm nhiệt độ đạt cực đại chậm dần từ mặt trên cùng qua các lớp giữa, muộn nhất
CT 2
là mặt dưới cùng. Phân bố nhiệt độ tại mọi thời điểm trong ngày luôn là đường cong, hình 5.
Hình 2. Toàn cảnh diễn biến nhiệt độ sau 192 giờ
Hình 3. Thay đổi nhiệt độ 14 nút sau 192 giờ
Hình 4. Thay đổi nhiệt độ tấm BT 1 ngày đêm
Hình 5. Phân bố nhiệt độ tại các thời điểm
3.3. Trao đổi nhiệt của tấm bê tông
L
Miền bi nén
Đường cong nhiệt độ thực
Miền bi kéo
Tm1
Đương nhiệt độ trung bình tuyến tính TTB
Tm2
Thiêt diện phẳng sau biến dạng
Thiêt diên phẳng trước biến dạng
Trao đổi nhiệt của tấm BT với môi trường trong một ngày đêm điển hình thể hiện trên hình 6. Có thể thấy: Trong một ngày đêm, tấm BT nhận nhiệt trong 11 tiếng, từ 7h sáng đến 17h chiều (q > 0), thải nhiệt trong 13 tiếng, từ 18h đến 6h sáng hôm sau (q < 0). Tốc độ nhận nhiệt tăng rất nhanh vào buổi sáng, giảm vào chiều. Nhiệt nhận lớn nhất lúc 11,12 h. Tốc độ toả nhiệt toả nhiệt ra môi trường khá đồng đều.
Hình 7
Hình 6. Lượng nhiệt tấm BT trao đổi với môi trường trong một ngày đêm
3.4. Nhận định về biến dạng nhiệt cục bộ của tấm BT
a. Biến dạng nhiệt cục bộ : Theo lý thuyết biến dạng nhiệt [2], [3], khi phân bố nhiệt độ là
đường cong, thì trong vật liệu sẽ xuất hiện biến dạng nhiệt cục bộ εT và ứng suất nhiệt riêng σT.
Nếu εT > εKth (εKth là biến dạng kéo tới hạn), (hoặc σT > [σ] ([σ] là ứng suất cho phép)), vật liệu
sẽ bị rạn nứt, phá huỷ. Biến dạng nhiệt cục bộ xác định bởi εT = βTΔT (26); với βT là hệ số giãn
nở nhiệt, ΔT là chênh lệch nhiệt độ cục bộ: ΔT(x) = TTB(x) – T(x) (27); với T(x) là nhiệt độ
T
T
m1
m2
thực, là đường cong, TTB(x) là đường nhiệt độ trung bình tuyến tính, đó là đường thẳng song
T
T
x
=
−
TB
m1
− L
CT 2
song với thiết diện phẳng sau biến dạng [2],[3], hình 7, xác định bởi:
(28); Tm1 , Tm2 là trị số nhiệt độ trung bình tuyến tính hai tại mặt 1 và 2 của tấm BT. Từ (26)
thấy, khu vực có TTB(x)>T(x), sẽ bị kéo, ngược lại, khu vực có TTB(x) b. Đặc điểm các miền bị kéo và bị nén xuất hiện trong tấm bê tông: Do phân bố nhiệt độ trong tấm BT luôn là đường cong, nên trong tấm BT luôn xuất hiện các miền bị kéo và nén. Một số trong các kết quả tính toán được thể hiện trên đồ thị, hình 8÷11. Có thể rút ra các nhận xét: 1. Nơi phân cách giữa hai miền bị kéo và bị nén là vị trí có ΔT=0. Vị trí này thay đổi liên tục, nghĩa là miền bị kéo và nén trong tấm BT trong ngày biến đổi liên tục và luân phiên nhau. 2. Từ 7h sáng đến 13h chiều, lớp mặt trên dày 5cm tấm BT bị nén do nhận nhiệt rất mạnh từ môi trường (ΔT<0). Lớp kế tiếp dày 15 - 20cm bị kéo do ΔT>0 , lớp dưới cùng dày bị nén. 3. Từ 13h chiều, phần trên tấm BT bị kéo do tốc độ nhận nhiệt giảm. Sau 16h miền bị kéo thu hẹp lại còn khoảng 5cm tại mặt trên cho đến 4h sáng, sau đó chuyển sang bị nén. Các lớp kế tiếp bị nén từ sau 18h đến 6h sáng hôm sau. c. Biến dạng nhiệt cục bộ trong tấm bê tông: Biến dạng nhiệt cục bộ (BDCB) tại các vị trí
trong tấm BT theo thời gian trong ngày, với βT = 10-5 [4], được thể hiện trên các hình 12 ÷ 15. Biến dạng kéo tới hạn εKth phụ thuộc vào nhiều yếu tố và loại bê tông, ở đây lấy giả định theo
[4], BT mác 200 có εKth = 2,4.10-4. Từ các đồ thị có thể rút ra các nhận xét sau: 1. Các lớp trong tấm BT luôn bị kéo, nén luân phiên nhau và thay đổi liên tục trong một ngày đêm, có thể gây nên hiện tượng mỏi nhiệt. 2. Biến dạng nhiệt cục bộ do bị kéo lớn nhất tại bề mặt εT = 0,5467.10-4 vào lúc 18h nhỏ
hơn biến dạng kéo tới hạn, hình 14, nên không gây rạn nứt phá huỷ, nhưng do bị kéo, nén luân phiên gây mỏi nhiệt, lâu dần bề mặt rất dễ bị rạn nứt. CT 2 IV. KẾT LUẬN Việc khảo sát trạng thái nhiệt tấm BT bằng phương pháp PTHH, cho phép rút ra kết luận sau: - Với cùng bước thời gian, kích thước hình học phần tử của cùng một bài toán [2], khảo sát bằng phương pháp PTHH cho các kết quả chính xác như phương pháp SPHH, các nhận định hoàn toàn phù hợp nhau . - Phương pháp PTHH có ưu điểm không phải tính thể tích, diện tích phân tố như phương pháp SPHH. - Phương pháp PTHH có thể áp dụng cho các vật thể có hình dạng bất quy tắc nên khả năng tính toán mở rộng hơn so với phương pháp SPHH. [1]. Trần Đình Bửu, Nguyễn Quang Chiêu. Khai thác và sửa chữa đường ô tô, NXB ĐH-THCN 1984. [2]. Trịnh Văn Quang. Kết quả tính toán trạng thái nhiệt của tấm bêtông dưới tác động của điều kiện khí
hậu thay đổi. Tạp chí Cầu Đường Việt nam. Số 11,12, 2001. [3]. С.А.Фрид.температурные напряжения в бетонных и железобетонных конструкциях
гидротехнических сооружений государствнное. Энергетическое издтелЬство. Москва 1959. [4]. Lê Văn Cung. Khống chế nhiệt độ đập bê tông Thác bà. Báo cáo Hội nghị KHXD 1985. [5]. RW Lewis, P.Nithiharasu and Seetharamu. Fundametals of The Finite Element Method for heat and fluid flow. John Wiley & Sons, Ltd. 2004♦ CT 2Hình 8. Lúc 10h
Hình 9. Lúc 12h trưa
Hình 10. Lúc 16 giờ
Hình 11. Lúc 20h
Hình 12. BDCB lúc 6h Hình 13. BDCB lúc10h
Hình 14. BDCB lúc18h Hình 15. BDCB lúc 22h
Tài liệu tham khảo
