7 Cấu trúc của luận án
3.3. Phân tích độ nhạy của các yếu tố ảnh hưởng tới mô đun độngcủa bê tông nhựa
3.3.2.3. Mơ hình Hirsch
Mơ hình Hirsch được phát triển bởi Christensen và các cộng sự trên cơ sở nghiên cứu thực nghiệm của 18 hỗn hợp BTN sử dụng 8 loại bitum và 5 cỡ hạt và thành phần hạt khác nhau, xác định được là 206 giá trị của |E*| theo thực nghiệm [40], [65], [83]. Mối quan hệ giữa |E*| và các thơng số đầu vào theo mơ hình Hirsch được đề cập trong phương trình 2.15 (chương 2).
3.3.2.4. Xác định dạng phân bố của các biến đầu vào trong các mơ hình dự báo |E*|
Dạng phân bố của các thông số đầu vào được xác định với sự trợ giúp của phần mềm Oracle Crystal Ball (OCB). Tùy thuộc vào cơ sở dữ liệu cụ thể của các thơng số đầu vào mà có các dạng phân bố khác nhau. Về nguyên tắc, để xác định dạng phân bố của dãy số liệu bất kỳ càng có nhiều dữ liệu càng tốt, nhưng nếu khơng có điều kiện thu thập nhiều dữ liệu thì phải có tối thiểu cần phải có 15 giá trị của mỗi thơng số đầu vào. Hình 3.7b minh họa giao diện của phần mềm Oracle crystal ball (OCB) khi được tích hợp vào microsfot excel (MS excel 2016).
Hình 3.7: Phần mềm Orcral Crystal Ball tích hợp vào MS excel 2016
Trên cơ sở kết quả thực nghiệm (38x4 = 152 giá trị của đặc tính thể tích gồm Va, Vbeff, VFA, VMA), 6 cấp phối ( 3 cấp phối BTNC 12.5 và 3 cấp phối BTNC 19), 120 giá trị của (|G*|), 120 giá trị của độ nhớt (η), 120 giá trị của góc pha (δb) tại các nhiệt độ -5oC, 10oC , 25oC, 40oC và 60oC và các giá trị của tần số (0.1Hz, 0.5 Hz, 1Hz, 5Hz, 10 Hz và 25 Hz) trong thí nghiệm |E*| ứng với các nhiệt độ tương ứng xác định được dạng phân bố của các thông số đầu vào với các thông số phân bố được tổng hợp trong bảng 3.11.
Bảng 3.11: Các thông số đầu vào sử dụng để phân tích độ nhạy Thơng số Dạng phân bố Các thông số của dạng phân bố Thông số Dạng phân bố Các thông số của dạng phân bố Va(%) phân bố chuẩn logarit
(lognormal distribution)
Vị trí (location)=0, giá trị trung bình (mean)=4.48, độ lệch chuẩn Std.dev=0.91.
Vbeff(%) Phân bố Logistic Giá trị trung bình 10.36, tỷ lệ 0.4. VFA(%) Phân bố Logistic Giá trị trung bình 70.07, tỷ lệ 2.9 VMA(%) Phân bố Beta (Beta
distribution),
Giá trị nhỏ nhất =13.01; giá trị lớn nhất =25.76, alpha = 4.314; Beta =25.639
P200(%) Phân bố đều Giá trị nhỏ nhất: 4.5; giá trị lớn nhất 8.5 P4 (%) Phân bố đều Giá trị nhỏ nhất: 40; giá trị lớn nhất 63 P3/8 (%) Phân bố đều Giá trị nhỏ nhất: 15; giá trị lớn nhất 41 P3/4 (%) Phân bố đều Giá trị nhỏ nhất: 0; giá trị lớn nhất 9 Độ nhớt (η,
10^6Poise)
phân bố chuẩn logarit (lognormal distribution)
Vị trí (location)=0, giá trị trung bình (mean)= 6241.852006, độ lệch chuẩn Std.dev= 24051992.802
|G*|(psi) phân bố chuẩn logarit (lognormal distribution)
Vị trí (location)=0, giá trị trung bình (mean)= 371672.18437, độ lệch chuẩn Std.dev= 2951548389.166
Thông số Dạng phân bố Các thông số của dạng phân bố Góc pha
(δb)
Phân bố beta Giá trị nhỏ nhất: 16.722; giá trị lớn nhất 91.64847, anpha: 1.19247; beta : 0.7708126
Tần số f(Hz)
phân bố chuẩn logarit (lognormal distribution)
Vị trí (location)=0, giá trị trung bình (mean)= 11.82, độ lệch chuẩn Std.dev= 69.09
3.3.2.5. Phân tích độ nhạy của các yếu tố ảnh hưởng tới |E*| bằng phương pháp Mô phỏng Monte Carlo
Sau khi xác định được dạng phân bố của các thông số đầu vào như đã nêu ở mục 3.4.2, tiến hành thực hiện mô phỏng Monte Carlo để phân tích độ nhạy tổng thể của các yếu tố ảnh hưởng tới giá trị của |E*| với sự trợ giúp của phần mềm (OCB).
Kết quả phân tích độ nhạy của các thơng số đầu vào ảnh hưởng tới giá trị |E*| của BTNC với số 3000 lần phân tích được thể hiện như các Hình từ 3.8 tới Hình 3.10.
Hình 3.8: Độ nhạy của log|E*| trong mơ
hình Witczk truyền thống
Hình 3.9: Độ nhạy của log(|E*|) trong mơ
hình Witczk Cải tiến
Nhận xét:
Qua các Hình 3.8 tới Hình 3.10 cho thấy có sự thống nhất với nhau trong cả 3 mơ hình là tính chất của vật liệu bitum (|G*|, , góc pha) có ảnh hưởng nhiều nhất tới giá trị của |E*|. Tiếp theo đó là ảnh hưởng của cấp phối thiết kế (P4, P3/8, P3/4, P200) và hàm lượng bitum tối ưu được thể hiện thông qua các đặc tính thể tích (Va, VMA, VFA, Vbeff).
Hình 3.10: Độ nhạy của |E*| trong mô hình
Hirsh
Các kết quả phân tích độ nhạy của các dữ liệu trong nghiên cứu gần giống với các kết quả phân tích độ nhạy trên cơ sở dữ liệu thu thập được trong thực tế sản xuất tại dự án ở Hải Phịng (các Hình 3.11 – Hình 3.13) cũng như một số nghiên cứu khác trên thế giới về phân tích độ nhạy của các thông số đầu vào ảnh hưởng tới |E*| cho các mơ hình Witczak, mơ hình Hirsh hay mơ hình dự báo |E*| đã được công bố [18], [28], [61].
87.9% 11.6% -0.2% 0.1% -0.1% 0.0% 0.0% 0.0% -100.0% 0.0% 100.0% η_Viscous f Vbeff= P3/4 Va= P3/8 P200 P4 Sensitivity: log(E*)= 96.7% -2.7% -0.3% 0.2% -0.1% 0.0% 0.0% 0.0% -100.0% -50.0% 0.0% 50.0% 100.0% Gsao denta Va Vbeff P_34 P_38 P_4 P_200 Sensitivity: log(E*)= 0.999 -0.040 0.012 -1.000 -0.500 0.000 0.500 1.000 Gsao VMA VFA= Sensitivity: E*=
Hình 3.11: Độ nhạy của log|E*| trong mơ
hình Witczk truyền thống
Hình 3.12: Độ nhạy của log(|E*|) trong mơ
hình Witczk Cải tiến
Lưu ý:
Dự án tại Hải Phòng chỉ thu thập các dữ liệu để phân tích độ nhạy cho BTNC 12.5 nên trong các biểu đồ độ nhạy khơng khơng có thành phần P3/4 ảnh hưởng tới giá trị của (|E*|). (lượng hạt tích lũy trên sáng 19mm).
Hình 3.13: Độ nhạy của |E*| trong mơ hình
Hirsh
Lý do tính chất của bitum có ảnh hưởng nhiều nhất tới giá trị |E*| của BTN là bitum là vật liệu đàn nhớt, ứng xử của bitum phụ thuộc cả vào nhiệt độ và tần số tác dụng của tải trọng rất nhiều. Giá trị |G*| của bitum có thể thay đổi trong phạm vi rất lớn giữa các mức nhiệt độ cao và nhiệt độ thấp. Ví dụ tỷ lệ giữa |G*| của bitum 60/70 tại tần số 1 Hz ở 10oC và 50oC là 1498.33 lần.
3.4. Nghiên cứu thực nghiệm xác định mô động của bê tông nhựa chặt ở Việt Nam Việt Nam
Việc nghiên cứu thực nghiệm xác định mô đun động nhằm mục đích xây dựng mối quan hệ thực nghiệm giữa tính chất của vật liệu bitum sử dụng với mô đun động của các loại
84.4% 15.2% -0.2% -0.1% 0.0% 0.0% 0.0% -100.0% -50.0% 0.0% 50.0% 100.0% η f P4 Va= P200 P3/8 Vbeff= Sensitivity: log(E*)= 97.8% -2.1% -0.1% -0.1% 0.0% 0.0% 0.0% -100.0% -50.0% 0.0% 50.0% 100.0% Gsao denta Va= P200 P3/8 Vbeff= P4 Sensitivity: log(E*)= 99.9% 0.1% 0.0% -100.0% -50.0% 0.0% 50.0% 100.0% Gsao VFA= VMA Sensitivity: E*=
BTNC trên cơ sở nền tảng là các mơ hình dự báo |E*| của Hoa Kỳ theo điều kiện vật liệu của Việt Nam.
Trên cơ sở các loại vật liệu, cấp phối đã chọn ở mục 3.1, thiết kế thí nghiệm cho BTNC 12.5 và BTNC 19 theo các bảng 3.9 và bảng 3.10 ở mục 3.2. Các bước tiếp theo cần thực hiện để nghiên cứu thực nghiệm xác định |E*| của các loại BTNC và đề xuất các hệ số trong các mơ hình dự báo |E*| theo điều kiện Việt Nam như sau.
3.4.1. Chuẩn bị vật liệu, thiết bị và đúc mẫu phục vụ cơng tác thí nghiệm mơ đun
động của các loại bê tông nhựa
3.4.1.1. Vật liệu và thiết bị thí nghiệm.
Vật liệu
Vật liệu thí nghiệm gồm bốn loại bitum (60/70; 40/50; 35/50 và PMBIII) và ba loại đá gồm đá vôi ở mỏ đá Phú Hà, Thanh Liêm- Hà Nam, đá bazan (Hataco, Quốc Oai Hà Nội), đá granit (Mỏ Khe Dầu, Tỉnh Hà Tĩnh), và bột khoáng sử dụng đá vôi của mỏ Gọng Vối (Thanh Sơn, huyện Kim Bảng, tỉnh Hà Nam).
Cấp phối thiết kế cho BTNC 12.5 và BTNC 19 tuân theo các yêu cầu của [2] và được tổng hợp như bảng 3.5 của mục 3.1. Phương pháp thiết kế BTN cũng như chuẩn bị các mẫu thí nghiệm thực hiện theo phương pháp Marshall với hàm lượng bitum và số lượng mẫu cần chuẩn bị trên cơ sở thiết kế thí nghiệm ở mục 3.2.
Theo phương pháp thiết kế thí nghiệm của Taguchi, mục 3.2 sẽ có 36 tổ mẫu BTNC ứng với các loại (đá, cấp phối và hàm lượng bitum khác nhau) được tổng hợp trong các bảng 3.9 và bảng 3.10. Tuy nhiên, do có đúc thêm 2 tổ mẫu BTNC để kiểm chứng tại hàm lượng bitum 5.2% (BTNC 12.5 với đá vôi và đá granit cho cấp phối JMF1), như vậy có tổng cộng 38 tổ mẫu để thí nghiệm |E*|.
Thiết bị thí nghiệm
Việc thí nghiệm |E*| được tiến hành bằng thiết bị CRT-UTM-NU của phịng thí nghiệm vật liệu xây dựng của trường ĐHGTVT. được nhập về trường năm 2009. Thiết bị này đã được kiểm định và có khả năng thực hiện các thí nghiệm liên quan tới BTN theo các tiêu chuẩn AASHTO TP-62/TP 79, EN 12697-25/26, ASTM D7369/D4123/D349 [27], [64].
3.4.1.2. Thí nghiệm xác định mơ đun động
Tiêu chuẩn áp dụng
Hiện nay, theo các tiêu chuẩn của Hoa Kỳ, có 3 tiêu chuẩn khác nhau dùng để xác định giá trị mơ đun của BTN, đó là:
Tiêu chuẩn AASHTO TP 79-13 [66]. Tiêu chuẩn AASHTO T342-11 [67].
Tiêu chuẩn AASHTO TP 62-07(2009) [68].
Thiết bị của Trường Đại học Giao thông vận tải phù hợp để thực hiện thí nghiệm theo tiêu chuẩn AASHTO TP62, do vậy tiêu chuẩn AASHTO TP-62 được áp dụng để thực hiện thí nghiệm mơ đun động của các mẫu BTN.
Xác định phạm vi nhiệt độ và tần số trong thí nghiệm |E*|
Theo tiêu chuẩn AASHTO TP-62, để thí nghiệm |E*| cần xác định các thơng số như tần số tác dụng của tải trọng, nhiệt độ thí nghiệm… Do đặc thù thiết bị thí nghiệm CRT- UTM-NU của Trường ĐHGT cài đặt các giá trị mặc định sẵn với tần số tải trọng tác dụng cố mặc là 0.1Hz, 0.5Hz, 1Hz, 5Hz,10Hz và 25Hz nên các tần số này sẽ thực được sử dụng trong thí nghiệm |E*| của BTN. Về nhiệt độ, các giá trị nhiệt độ thí nghiệm |E*| sẽ gồm các nhiệt độ từ thấp tới cao, đủ để kết hợp với phạm vi tần số tác dụng của tải trọng để xây dựng được các đường cong chủ của mô đun động của các hỗn hợp BTNC khác nhau đạt độ chính xác phù hợp theo các tiêu chuẩn thống kê đề cập trong bảng 2.3. Trên nguyên tắc này, các nhiệt độ sử dụng trong thí nghiệm E* gồm có 10oC, 25oC, 40oC, và 55oC và việc thí nghiệm được tiến hành theo thứ tự từ nhiệt độ thấp nhất tới nhiệt độ cao nhất, và tại mỗi nhiệt độ sẽ tiến hành thí nghiệm theo thứ tự từ tần số lớn nhất tới tần số nhỏ nhất [68].
Xác định mức ứng suất và số chu kỳ tác dụng của tải trọng
Các mức độ ứng suất động tác dụng vào mẫu phụ thuộc vào độ cứng của mẫu, tùy theo các giá trị của nhiệt độ trong thí nghiệm mẫu mà mức độ ứng suất tác dụng vào mẫu thay đổi với nguyên tắc nhiệt độ càng thấp thì giá trị ứng suất sử dụng trong thí nghiệm càng lớn và ngược lại, giá trị ứng suất tác dụng vào mẫu cần điều chỉnh sao cho mức độ biến dạng dọc trục trong khoảng 50 tới 150 microstrain. Bảng 3.13 đưa ra các mức ứng suất động điển hình tùy theo các giá
trị của nhiệt độ thí nghiệm khác nhau [68]
Bảng 3.13: Các mức độ ứng suất động tùy theo nhiệt độ thí nghiệm điển hình Nhiệt độ thí nghiệm Phạm vi ứng suất Nhiệt độ thí nghiệm Phạm vi ứng suất (oC) (oF) (KPa) (Psi) -10 14 1400 -2800 200-400 4 40 700-1400 100-200 21 70 350-700 50-100 37 100 140-250 20-50 54 130 35-70 5-10
Số chu kỳ tác dụng của tải trọng tùy thuộc vào các giá trị tần số sử dụng trong thí nghiệm, tần số bé thì số chu kỳ tác dụng cũng nhỏ, bảng 3.14 tóm tắt số chu kỳ tác dụng của tải trọng tùy theo các giá trị tần số khác nhau [68].
Bảng 3.14: Số các chu kỳ thí nghiệm tùy theo các giá trị của tần số khác nhau [68] Tấn số (Hz) 0.1 0.5 1 5 10 25
Số chu kỳ 15 15 20 100 200 200
Các giá trị của |E*| ở các tần số và nhiệt độ khác nếu cần biết, sẽ được xác định bằng việc xây dựng đường các cong chủ (Master curve) của |E*| của các loại BTNC tương ứng.
Tóm tắt phương pháp thí nghiệm
Sau khi cơng tác chuẩn bị mẫu đã xong, các mẫu BTN được gia công hai đầu mẫu để đảm bảo các đầu mẫu có mặt phẳng tiếp xúc tốt với bộ phận gia tải, hạn chế các sai số liên quan tới bề mặt mẫu. Các mẫu sau đó được cho vào tủ kiểm sốt nhiệt độ với thời gian duy trì nhiệt độ như bảng Bảng 3.15.
Bảng 3.15: Thời gian duy trì mẫu thí nghiệm ở các nhiệt độ khác nhau [68] Nhiệt độ mẫu
(oC)
Thời gian với mẫu chưa qua thí nghiệm (nhiệt độ phịng là 25oC), (h)
Mẫu đã được kiểm tra ở mức nhiệt độ trước đó 10 Qua đêm 4h hoặc qua đêm
25 1 3
40 2 2
Sau khi các mẫu BTNC đã duy trì đủ thời gian cho từng mức nhiệt độ khác nhau ở trên, các mẫu sẽ được tiến hành thí nghiệm bằng việc tác dụng một tải trọng hình sin với độ lớn và chu kỳ tác dụng tùy theo mức nhiệt độ và tần số như các bảng 3.13 và bảng 3.14. Các bộ phận cảm biến sẽ ghi lại kết quả của biến dạng phục hồi dọc trục của mẫu và xuất kết quả qua máy tính. Giá trị mơ đun động của BTNC (|E*|) sẽ được xác định theo công thức 1.10 của chương 1. Chi tiết kế quả thí nghiệm |E*| như phục lục 3.
3.5. Xây dựng đường cong chủ (Master curve) của mô đun động |E*|
Mục đích của việc xây dựng đường cong chủ của |E*| là có thể xác định được các giá trị |E*| tại nhiệt độ và tần số bất kỳ, điều này có ý nghĩa kinh tế-kỹ thuật vì khơng cần làm thí nghiệm |E*| nhiều nhưng vẫn có thể xác định được các giá trị |E*| tại bất kỳ mức nhiệt độ/tần số mong muốn. Phương trình đường cong chủ của |E*| là một hàm hình sin (phương trình 3.2) và đồ thị logarit các hệ số dịch chuyển (shift factors) chuyển của |E*| theo nhiệt độ được minh họa như hình 3.15. Nguyên lý và cách xây dựng đường cong chủ của |E*| đã được trình bày trong mục 1.4.4 của chương 1.
* * max 1 1 log - 19.14714 (log | | - ) log | | 1 r m m E f T T E E e (3.2) Trong đó:
|E*| - Mô đun động, ksi; δm, β, γ và * max
log | E |
là các thông số phù hợp.
δm, *
max
log | E |: Lần lượt là giá trị mô đun động
nhỏ nhất và lớn nhất của BTN, ksi. Hình 3.15: Biểu đồ hệ số dịch chuyển theo
nhiệt độ [23]
f: Tần số tác dụng của tải trọng tại nhiệt độ thí nghiệm, Hz. T: Nhiệt độ thí nghiệm, (oK) và Tr –Nhiệt độ tham chiếu, (oK).
E: Năng lượng kích hoạt (the activation energy), là năng lượng cần phải cung cấp để đạt được tốc độ và trạng thái chuyển tiếp cho một phản ứng hóa học hoặc trạng thái chuyển tiếp vật lý xảy ra. Tiêu chuẩn AASHTO. Giá trị của của Eđược xử lý là một
thông số phù hợp “treated as a fitting parameter”. Giá trị này cùng với các thông số khác (δm, β, γ và *
max
log | E |) là các thông số phù hợp được xác định bằng phương pháp tính tốn lặp dần (sử dụng phần mềm MS Excel), giá trị cần tìm đảm bảo sao cho sai số bình phương giữa giá trị |E*| thực nghiệm và giá trị |E*| tính theo cơng thức 3.2 nhỏ nhất. Giá trị của δm và (log | E*max | - δm) phụ thuộc vào thành phần hạt, hàm lượng bitum, và độ rỗng dư Va, còn giá trị của β và ɤ phụ thuộc vào đặc tính của bitum và độ lớn của δm và ( *
max
log | E |- δm). Giá trị của mô đun giới hạn lớn nhất của BTN được dự báo theo các đặc tính thể tích của BTN theo phương trình 2.15 và phương trình 2.17.
Khi xây dựng đường cong chủ, tại các nhiệt độ lớn hơn nhiệt độ tham chiếu, các đường cong |E*| sẽ dịch sang trái và ngược lại.
Phương trình tổng quát xác định hệ các số dịch chuyển của |E*| ở các nhiệt độ khác nhau (phương trình 3.3) [23]. 1 1 log 19.14714 T R E T T (3.3)
Với: αT – hệ số dịch chuyển tại nhiệt độ kiểm tra (T).