Bài viết đã sử dụng phương pháp quy hoạch thực nghiệm bậc hai để dự đoán và mô phỏng ảnh hưởng của tỷ lệ nước – xi măng (N/X) và cát – chất kết dính đến các hàm mục tiêu là độ chảy xòe của hỗn hợp bê tông và cường độ nén của mẫu bê tông chất lượng cao hạt mịn.
Journal of Mining and Earth Sciences Vol 62, Issue (2020) 73 - 84 73 Using an experimental model to study the properties of fine-grained high-performance concrete Lam Tang Van 1,*, Dien Vu Kim 2, Hung Ngo Xuan 1, Bulgakov Boris Igorevich 2, Chien Minh Do 3, Duong Van Nguyen Faculty of Civil Engineering, Ha Noi University of Mining and Geology, Hanoi, Vietnam Moscow State University of Construction, Moscow, Russian Federation College of Industrial and Constructional, Quangninh, Vietnam Department of Technology of Binder and Concrete, ARTICLE INFO ABSTRACT Article history: Received 15th Nov 2020 Accepted 21st Jan 2021 Available online 28th Feb 2021 This study uses the mathematical method of two-factors rotatable central compositional planning to predict and simulate the effect of the ratio of water-cement (N/X) and sand - binder (C/CKD) as the input parameters on the objective functions of the spreading flow of concrete mixtures and the compressive strength of fine-grained high-performence concrete (FGHPC) at 28 days Results of the study showed that, from the material source in Vietnam, it is possible to FGHPC with a flow of 18.5 cm in the mini cone, its compressive and flexural strengths at the age of 28 days are 68.5 MPa and 6.13 MPa, respectively Furthermore, from the obtained objective functions, it has been shown that the both two-input parameters have a significant influence on the values of the experimental models Particularly, using Matlab software is showed the expression surface, the contour line of the experimental models, and determined the maximum value of compressive strength of FGHPC at this age of 69.84 MPa at N/X=0.326 and C/CKD=1.315 The contribution of this study is to obtain regression functions to predict the mechanical-physical properties of FGHPC that will be used in the next in-depth studies Keywords: Compressive strength, Fine-grained high-performent concrete, Objective function, Parameter Copyright © 2021 Hanoi University of Mining and Geology All rights reserved _ *Corresponding author E - mail: lamvantang@gmail.com DOI: 10.46326/JMES.2021.62(1).09 74 Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất Tập 62, Kỳ (2021) 73 - 84 Sử dụng mơ hình thực nghiệm để nghiên cứu tính chất bê tơng chất lượng cao hạt mịn Tăng Văn Lâm 1,*, Vũ Kim Diến 2, Ngô Xuân Hùng 1, Bulgakov Boris Igorevich 2, Đỗ Minh Chiến 3, Nguyễn Văn Dương Khoa Xây dựng, Trường Đại học Bộ môn Công nghệ bê tông Mỏ - Địa chất, Hà Nội, Việt Nam chất kết dính, Trường Đại học Xây dựng Quốc gia Mátxcơva, Mátxcơva, THƠNG TIN BÀI BÁO TĨM TẮT Quá trình: Nhận 15/11/2020 Chấp nhận 21/1/2021 Đăng online 28/02/2021 Bài báo sử dụng phương pháp quy hoạch thực nghiệm bậc hai để dự đốn mơ ảnh hưởng tỷ lệ nước – xi măng (N/X) cát – chất kết dính (C/CKD) đến hàm mục tiêu độ chảy xòe hỗn hợp bê tông cường độ nén mẫu bê tông chất lượng cao hạt mịn (BTCLCHM) Kết nghiên cứu cho thấy, từ nguồn vật liệu Việt Nam chế tạo BTCLCHM với độ chảy xịe mini 18,5 cm, cường độ nén cường độ kéo uốn tuổi 28 ngày 68,5 MPa 6,13 MPa Mặt khác, từ hàm mục tiêu hai biến N/X C/CKD có ảnh hưởng đáng kể đến mơ hình thực nghiệm Sử dụng phần mềm Matlab biểu diễn bề mặt biểu đường đồng mức đối tượng nghiên cứu Đồng thời, giá trị cường độ nén lớn tuổi 28 ngày mẫu BTCLCHM xác định 69,84 MPa N/X=0,326 C/CKD=1,315 Đóng góp nghiên cứu thu hàm hồi quy để dự đốn tính chất – lý BTCLCHM sử dụng nghiên cứu chuyên sâu Liên bang Nga Trường Cao đẳng Công nghiệp Xây dựng, Quảng Ninh, Việt Nam Từ khóa: Cường độ nén, Bê tơng chất lượng cao hạt mịn, Biến ảnh hưởng, Hàm mục tiêu © 2020 Trường Đại học Mỏ - Địa chất Tất quyền bảo đảm Mở đầu Bê tông hạt mịn chủng loại bê tông mà thành phần cốt liệu có cốt liệu nhỏ với kích thước từ 0,14÷5 mm Do thành phần chất kết dính loại bê tơng hạt mịn thường lớn 1,5÷3 lần so với bê tơng xi măng thơng thường (Bazhenov, 2011; Bazhenov nnk, 1998) Bê tông _ *Tác giả liên hệ E - mail: lamvantang@gmail.com DOI: 10.46326/JMES.2021.62(1).09 chất lượng cao hạt mịn loại bê tông hạt mịn đặc biệt, nghiên cứu chế tạo từ hỗn hợp xi măng Pc lăng cát vàng thơng thường, kết hợp sử dụng phụ gia mịn siêu mịn với tỷ lệ đảm bảo thành phần hạt hợp lý vi cấu trúc đặc hỗn hợp bê tông (Nguyễn Như Quý Mai Quế Anh, 2020; Nguyễn Văn Tuấn nnk, 2018; Bazhenov and nnk, 2006) Hiện nay, nhiều nước giới có Việt Nam nghiên cứu, chế tạo phát triển nhiều loại bê tông chất lượng cao hạt mịn với tính chất cơ-lý khác tùy theo mục đích sử dụng Tăng Văn Lâm nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 62(1), 73-84 Trong năm gần đây, bê tơng chất lượng cao nói chung bê tơng chất lượng cao hạt mịn nói riêng ngày sử dụng phổ biến đa dạng hơn, từ kết cấu chịu lực đến vật liệu hoàn thiện bề mặt nhà cao tầng (Lam Van Tang nnk, 2018a), cơng trình thủy (Tang Van Lam nnk, 2018b), cơng trình sở hạ tầng giao thơng biển nhiều cơng trình đặc biệt khác (Pham Duc Thang nnk, 2016; Lam Tang Van nnk, 2017; 2019c) Cùng chung với xu phát triển ngành xây dựng giới Việt Nam sử dụng nhiều loại bê tông hạt mịn xây dựng cơng trình giao thơng, cơng trình ngầm với hiệu cao như: gia cố mái dốc hồn thiện bề mặt cơng trình ngầm bê tông phun, kết hợp với neo đinh đất để gia cố tường vỏ chống đường hầm (Đào Viết Đoàn Tăng Văn Lâm, 2017; Phùng Mạnh Đắc, 2002; Nguyễn Quang Phích, 2002) Mặt khác, Việt Nam chưa có tiêu chuẩn hướng dẫn để tính tốn thiết kế thành phần bê tơng hạt mịn sử dụng chất kết dính từ xi măng Pc lăng chất kết dính khơng xi măng Các nghiên cứu bê tông chất lượng cao hạt mịn chủ yếu dựa vào sở phương pháp thể tích tuyệt đối để tính tốn ước lượng cấp phối ban đầu, sau thí nghiệm để kiểm chứng điều chỉnh thành phần hỗn hợp bê tông hạt mịn để thu tính chất cơ-lý yêu cầu Hướng nghiên có nhược điểm số lượng thí nghiệm thăm dị điều chỉnh lớn, không đánh giá tổng quát mối tương quan yếu tố ảnh hưởng đến tính chất bê tơng nghiên cứu (Lam Tang Van nnk, 2019a; 2019b) Nhiều nghiên cứu rằng, tính cơng tác hỗn hợp bê tơng cường độ nén mẫu bê tơng đạt cách điều chỉnh tính chất cấp phối bê tơng thơng qua việc tối ưu hóa thành phần cấp phối giảm độ xốp bê tông, sử dụng loại xi măng đặc biệt, phụ gia hóa dẻo, phụ gia phân tán mịn (Tang Van Lam and nnk, 2019e; Ngo Xuan Hung nnk, 2018) Do vậy, nghiên cứu chế tạo bê tông chất lượng cao hạt mịn, việc lựa chọn loại nguyên vật liệu, thành phần cấp phối tối ưu cần phải bảo đảm để thu hỗn hợp bê tông hạt mịn dễ thi công, có độ đặc cao, tốc độ rắn nhanh, cường độ ban đầu cường độ cuối đạt cao Hơn nữa, công nghệ bê tông đại có nhiều phương pháp tính tốn thiết kế mơ để tối ưu hố thành phần hỗn hợp bê tông 75 Cùng với hỗ trợ phần mềm tốn ứng dụng, cơng cụ máy tính cho phép lựa chọn ứng dụng phương pháp số hiệu quả, độ xác tính khả thi cao việc dự đốn tính chất hỗn hợp bê tông bê tông (Williams, 2013) Theo Tang Van Lam nnk, (2017); Nguyễn Minh Tuyển (2007) quy hoạch thực nghiệm sở mơ hình thống kê sử dụng để nghiên cứu nhiều đối tượng khác nhau, có đối tượng thành phần vật liệu ảnh hưởng đến tính chất hỗn hợp bê tông, bê tông Các trình xảy đối tượng nghiên cứu bê tông đặc trưng biến điều khiển, chúng có quan hệ nguyên nhân - kết với hàm mục tiêu Các biến đóng vai trị ngun nhân, gọi biến vào, biến phản ánh kết nguyên nhân gây ra, gọi biến Biến vào kiểm sốt điều khiển Trong nghiên cứu sử dụng phương pháp quy hoạch thực nghiệm bậc hai tâm xoay để mô ảnh hưởng hai biến tỷ lệ nước – xi măng (N/X) tỷ lệ cát – chất kết dính (C/CKD) đến độ chảy xịe cường độ nén mẫu bê tông chất lượng cao hạt mịn tuổi 28 ngày Đồng thời, bề mặt biểu hiện, đường đồng mức giá trị cực trị hàm hồi quy xác định nghiên cứu Đóng góp nghiên cứu thu hàm hồi quy để đốn tính chất – lý vật liệu nghiên cứu chuyên sâu tiếp theo, đặc biệt bê tông cường độ cao, chất lượng cao hạt mịn sử dụng chất kết dính khơng xi măng Vật liệu phương pháp nghiên cứu 2.1 Vật liệu sử dụng Vật liệu sử dụng nghiên cứu bao gồm: (1) Chất kết dính (CKD) bao gồm loại vật liệu sau: xi măng Poóc lăng PC40 "Hoàng Thạch" (X) thỏa mãn yêu cầu kỹ thuật tiêu chuẩn TCVN 2682:2009 GOST 31108-2016 (Nga), tro bay (FA) loại F nhà máy nhiệt điện Vũng Áng thỏa mãn yêu cầu TCVN 10302:2014, ASTM C618-03 GOST Р 56592-2015; xỉ lò cao hoạt hóa nghiền mịn (Xi) lấy trực tiếp khu cơng nghiệp luyện gang thép Hịa Phát (Kinh Mơn, Hải Dương) thỏa mãn theo TCVN 11586:2016 Silica fume SF-90 (SF90) Vina Pacific Thành phần hóa học tính chất vật lý xi măng, tro bay silica fume thể Bảng 76 Tăng Văn Lâm nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 62(1), 73-84 Bảng Thành phần hóa học tính chất vật lý silica fume SF-90, xỉ luyện kim "Hòa Phát", tro bay “Vũng Áng” xi măng Pc lăng PC40 "Hồng Thạch" Loại vật liệu Tro bay Xỉ luyện kim Silica Fume SF-90 Xi măng Poóc lăng Thành phần hóa học chất kết dính SiO2 (%) 54,2 36,3 91,6 20,4 Al2O3 (%) 23,3 12,6 2,2 4,4 Fe2O3 (%) 9,8 3,4 2,5 5,4 SO3 (%) 2,5 5,7 3,4 K2O (%) 1,4 0,4 1,2 Na2O (%) 1,1 0,3 0,5 0,3 MgO (%) 1,8 2,5 CaO (%) 1,4 40,1 0,7 60,2 Lượng nung (%) 4,5 1,2 2,5 2,2 Tính chất vật lý chất kết dính Tỷ diện bề mặt riêng (m2/g) 0,485 2,541 14,45 0,368 Khối lượng riêng (g/cm ) 2,35 2,92 2,15 3,15 Khối lượng thể tích khơ (kg/m3) 575 1450 760 1550 Thành phần hạt hỗn hợp chất kết dính xác định thể Hình (a)- Lượng sót tích lũy (b)- Phân bố thành phần hạt Hình Thành phần hạt Silica fume, xỉ luyện kim "Hòa Phát", tro bay "Vũng Áng" xi măng Pc lăng Bảng Tính chất vật lý cát vàng sông Lô Stt Chỉ tiêu Đơn vị tính Kết thí nghiệm Kích thước hạt mm 0,15 ÷ Khối lượng riêng g/cm3 2,65 Khối lượng thể tích đầm chặt kg/m3 1660 Khối lượng thể tích xốp kg/m3 1550 Độ ẩm % 3,5 Hàm lượng bụi, bùn, sét % 0,9 Mô đun độ lớn (Mk) - 3,1 Tạp chất hữu - Đạt Tăng Văn Lâm nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 62(1), 73-84 (2) Cốt liệu nhỏ sử dụng bê tông cát vàng sông Lô (C), loại hạt thô, chất lượng tốt, thỏa mãn yêu cầu tiêu chuẩn TCVN 7570:2006 GOST 8736-2014, sử dụng làm cốt liệu nhỏ hỗn hợp bê tơng Các tính chất vật lý cốt liệu nhỏ sử dụng thể Bảng (3) Phụ gia siêu dẻo SR 5000F «SilkRoad» (SR5000) có khối lượng riêng 1,12 g/m3 nhiệt độ 25±5 0C Đây loại phụ gia giảm nước tầm cao, hệ 3, có thành phần dựa gốc Polycarboxylate (4) Nước (N) sử dụng để làm nước trộn hỗn hợp bê tông bảo dưỡng mẫu thí nghiệm, thỏa mãn tiêu chuẩn TCVN 4506:2012 GOST 23732-2011 2.2 Phương pháp nghiên cứu - Sử dụng phương pháp laze máy nhiễu xạ "BT-9300z" xác định phân bố thành phần hạt thành phần hạt xi măng, tro bay, xỉ luyện kim silica fume SF-90 - Thành phần bê tông chất lượng cao hạt mịn tính tốn thiết kế theo phương pháp thể tích tuyệt đối kết hợp điều chỉnh thực nghiệm - Tính cơng tác hỗn hợp bê tơng xác định độ xịe tiêu chuẩn mini, kích thước 70x80x40 mm theo ASTM C1611 - 18 TCVN 3106:2007 - Cường độ nén bê tơng xác định mẫu hình trụ có đường kính D=150 mm chiều cao H=300 mm theo tiêu chuẩn TCVN 3105:1993, ASTM C39 GOST 10180-2012 - Mô ảnh hưởng xác định cấp phối tối ưu hỗn hợp bê tông thực theo phương pháp quy hoạch thực nghiệm bậc hai tâm xoay Box-Wilson 2.3 Phương pháp quy hoạch thực nghiệm Quy hoạch thực nghiệm thực chất sử dụng kết thực nghiệm trực tiếp xác định trước tuân theo quy luật xác suất Đây phương pháp nghiên cứu khoa học nhiều nhà nghiên cứu quan tâm (Williams, 2013; Nguyễn Minh Tuyển, 2007; Tang Van Lam nnk, 2019a; 2019d) Những ưu điểm phương pháp là: - Giảm đáng kể số lượng thí nghiệm, tiết kiệm thời gian kinh phí - Lượng thơng tin nhiều hơn, cụ thể nhờ đánh giá cách tương đối toàn diện ảnh hưởng nhân tố đến hàm mục tiêu - Thu mô hình thống kê thực nghiệm, cho phép đánh giá tranh thực nghiệm theo 77 tiêu chuẩn thống kê cho phép xét ảnh hưởng thông số với mức độ tin cậy cần thiết - Cho phép xác định điều kiện tối ưu đa nhân tố đối tượng nghiên cứu cách xác cơng cụ tốn học thay cho cách giải gần Trong số nghiên cứu gần cho thấy, phương pháp thực nghiệm bậc hai tâm xoay Box-Wilson cho phép thu mơ hình thực nghiệm, đối tượng thực nghiệm xác cách tăng số lượng thí nghiệm lặp lại trung tâm số điểm thực nghiệm đặc biệt cánh tay đòn α = = 1, 414 (Nguyễn Như Quý Mai Quế Anh, 2020; Lam Van Tang nnk 2018a; Tang Van Lam nnk, 2018b) Vị trí số lượng điểm thí nghiệm thực nghiệm bậc hai tâm xoay Box-Wilson hai biến ảnh hưởng Hình + + + - + - - - Hình Vị trí điểm thí nghiệm bậc hai tâm xoay Box-Wilson cho hai biến ảnh hưởng Ma trận thực nghiệm tương quan biến mã giá trị kế hoạch thực nghiệm bậc hai tâm xoay Box-Wilson dành cho hai biến đầu vào thể Bảng 3 Kết thảo luận 3.1 Cấp phối sở hỗn hợp bê tông hạt mịn chất lượng cao a) Các yêu cầu hỗn hợp bê tông hạt mịn chất lượng cao - Theo tiêu chuẩn TCVN 10306:2014, bê tông cường độ cao bê tơng có cường độ chịu nén đặc trưng 55 MPa lớn tuổi 28 ngày theo tiêu chuẩn ASTM C39 thí nghiệm mẫu hình trụ có đường kính D=150 mm chiều cao H=300 mm Từ sở đó, nghiên cứu tính tốn thiết kế thành phần cấp phối bê tông chất 78 Tăng Văn Lâm nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 62(1), 73-84 Bảng Ma trận thực nghiệm bậc hai tâm xoay Box-Wilson cho biến ảnh hưởng Các biến mã hóa ma trận thực nghiệm Hệ thống thí Số thí Giá trị thực nghiệm nghiệm nghiệm - Yj 0 1 2 12 21 22 +1 +1 +1 +1 +1 +1 Y1 Số lượng thí +1 -1 +1 -1 +1 +1 Y2 nghiệm +1 +1 -1 -1 +1 +1 Y3 lõi kế hoạch +1 -1 -1 +1 +1 +1 Y4 +1 +1,414 0 Y5 Số điểm thí +1 -1,414 0 Y6 nghiệm “cánh tay +1 +1,414 0 Y7 đòn” +1 -1,414 0 Y8 +1 0 0 Y9 Số điểm thí 10 +1 0 0 Y10 nghiệm lặp 11 +1 0 0 Y11 lại tâm kế 12 +1 0 0 Y12 hoạch 13 +1 0 0 Y13 lượng cao hạt mịn có cường độ nén yêu cầu tuổi 28 ngày 60 MPa, tính cơng tác tốt với độ xịe 15÷20 cm, sử dụng làm hỗn hợp bê tông bơm xây dựng cơng trình ven biển hải đảo Việt Nam - Silica fume SF-90 với kích thước hạt mịn cỡ nano chứa đến 91% SiO2 hoạt tính sử dụng với mục đích giảm lượng Ca(OH)2 tự bê tông tăng độ đặc vi cấu trúc Hàm lượng silica fume sử dụng 10% lượng dùng xi măng (Bazhenov, 2011) Tro bay nhiệt điện Vũng Áng xỉ luyện kim Hòa Phát sử dụng với hai mục đích chính: (i) phụ gia khoáng mịn, bổ sung thành phần hạt trơ với kích thước nhỏ 0,14 mm, phần cốt liệu mịn để tăng khung xương chịu lực giảm co ngót hỗn hợp bê tơng rắn chắc; (ii) thay phần xi măng Poóc lăng để giảm lượng nhiệt thủy hóa, giảm hình thành vết nứt kết cấu giảm giá thành sản phẩm Theo tiêu chuẩn ACI 211.4R-08, hàm lượng tro bay lựa chọn 20%, hàm lượng xỉ luyện kim sử dụng 40% lượng dùng xi măng - Lượng phụ gia siêu dẻo lấy 1,5% lượng dùng xi măng (Tang Van Lam nnk, 2018a; 2018b; 2019) Hàm lượng khơng khí hỗn hợp bê tơng 3% thể tích hỗn hợp bê tơng (Bazhenov nnk, 2006) b) Tính tốn cấp phối sở bê tông hạt mịn chất lượng cao theo phương pháp thể tích tuyệt đối Thiết kế cấp phối bê tơng theo phương pháp thể tích tuyệt đối áp dụng nhiều nghiên cứu bê tông chất lượng cao, bê tông chất lượng siêu cao Việt Nam (Nguyễn Văn Tuấn nnk, 2018) giới (Bazhenov nnk, 1998, Pham Duc Thang nnk, 2016) Theo phương pháp thể tích tuyệt đối, thể tích 1m3 bê tơng lèn chặt coi tổng thể tích nước, cốt liệu, xi măng, phụ gia khống, phụ gia siêu dẻo thể tích khơng khí vào q trình nhào trộn Biểu thức thể tích tuyệt đối hỗn hợp bê tơng phương pháp trình bày cơng thức (1) N N + X X + FA FA + Xi Xi + SF 90 SF 90 + C C + SR5000 SR 5000 + A = 1000 (1) Trong đó: N, X, FA, Xi, SF90, C, SR5000 - khối lượng nước, xi măng, tro bay, xỉ luyện kim, silicafume, cát phụ gia siêu dẻo N, X, FA, SF90, C, SR5000 - khối lượng riêng nước, xi măng, tro bay, xỉ luyện kim, silicafume, cát phụ gia siêu dẻo A - thể tích rỗng khơng khí vào hỗn hợp bê tơng, theo tài liệu (Bazhenov, 2011; 2006) thể tích khơng khí vào A = 3% Trên sở kết nghiên cứu tài liệu (Pham Duc Thang nnk, 2016; Lam Tang Van nnk, 2017; 2019c) kết khảo sát thực nghiệm sơ bộ, nghiên cứu chọn gốc hệ số tỷ lệ vật liệu Bảng Tính tốn theo phương pháp thể tích tuyệt đối dựa giá trị tỷ lệ vật liệu, thu cấp phối sơ hỗn hợp bê tơng hạt mịn có thành phần Bảng Thực nghiệm khảo sát sơ phịng thí nghiệm để xác định tính chất học bê tông hạt mịn chất lượng cao với cấp phối sở, kết thu thể Bảng Tăng Văn Lâm nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 62(1), 73-84 79 Bảng Các tỷ lệ vật liệu sử dụng Tỷ lệ C CKD Giá trị 1,3 N X 0,35 SF90 X 0,1 FA X 0,2 Xi X 0,4 SR5000 X A 0,015 3% Với chất kết dính: CKD = XM + FA + Xi + SF90 Bảng Cấp phối sơ hỗn hợp bê tông hạt mịn chất lượng cao tỷ lệ vật liệu sử dụng Vật liệu Cát vàng Xi măng PC40 Tro bay Xỉ luyện kim Silicafume SF-90 Chất kết dính Nước Phụ gia siêu dẻo SR 5000F Ký hiệu C X FA Xi SF90 CKD N SR5000 Nguồn cung cấp Sơng Lơ Hồng Thạch Nhiệt điện "Vũng Áng" Hịa Phát Vina Pacetic XM + FA + Xi + SF90 Nước máy SilkRoad Hàm lượng (kg/m3) 1202 544 109 218 54 925 190 8,2 Bảng Tính chất hỗn hợp bê tơng bê tơng cường độ cao Tính chất hỗn hợp bê tông N X 0,35 Cường độ nén bê tông (MPa) tuổi Độ chảy xịe Khối lượng thể ngày (mm) tích (kg/m3) 18,5 2324,8 ngày 36,8 3.2 Xây dựng mơ hình kế hoạch thực nghiệm Theo nhiều nghiên cứu (Nguyễn Như Quý Mai Quế Anh, 2020; Tang Van Lam nnk, 2017, Nguyễn Minh Tuyển, 2007) phương pháp quy hoạch thực nghiệm nhằm mục đích tạo mơ hình tốn học để mơ phỏng, phân tích đánh giá dự đoán ảnh hưởng biến đầu vào đến chất q trình tính chất đối tượng thực nghiệm, coi hàm mục tiêu đầu q trình nghiên cứu Mơ hình đối tượng thực nghiệm nghiên cứu mô tả Hình Hình Sơ đồ đối tượng nghiên cứu mơ hình thực nghiệm Phương trình tổng qt hàm mục tiêu bậc hai hai biến ảnh hưởng chương trình thực nghiệm có dạng sau: (2) Y = b0 + b11 + b22 + b312 + b412 + b522 14 ngày 28 ngày Cường độ kéo uốn tuổi 28 ngày (MPa) 50,4 63,2 68,5 6,13 Trong đó: Y – hàm mục tiêu mơ hình thực nghiệm; 1 2 – biến ảnh hưởng; b0, b1, b2, b3, b4, b5 – hệ số hàm mục tiêu a) Hàm mục tiêu: nghiên cứu này, hàm mục tiêu bậc hai lựa chọn gồm có: + Y1 – độ chảy xịe (cm) hỗn hợp bê tơng cường độ mini, với kích thước 70x80x40 mm; + Y2 – cường độ nén (MPa) mẫu bê tơng hình trụ có đường kính D = 150 mm chiều cao H = 300 mm tuổi 28 ngày b) Các biến ảnh hưởng: nghiên cứu chọn biến đầu vào, có ảnh hưởng lớn đến hàm mục tiêu, nhằm điều khiển hàm mục tiêu khoảng biến thiên chúng Bảng 7: 𝑁 + Tỷ lệ nước xi măng (𝑋 ) mã hoá biến 1 thay đổi 0,311÷0,389 + Tỷ lệ cát chất kết dính ( 𝐶 𝐶𝐾𝐷 ) mã hố biến 2 thay đổi 1,01÷1,59 Theo lý thuyết quy hoạch thực nghiệm, số thí nghiệm kế hoạch bậc hai tâm xoay Box – Wilson xác định theo công thức (3): N= 2k + 2*k + N0 = 22 + 2*2 + 5= 13 (3) Trong đó: k - số biến thực nghiệm; N0 - số thí nghiệm lặp lại tâm mơ hình nghiên cứu 80 Tăng Văn Lâm nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 62(1), 73-84 Trên sở cấp phối thực nghiệm 13 điểm hai biến đầu vào 1 2 trình bày Bảng Bên cạnh đó, giá trị trung bình độ xịe hỗn hợp bê tơng cường độ nén bê tông tuổi 28 ngày xác định trực tiếp từ thực nghiệm mô tả chi tiết Bảng 10 Bảng Khoảng biến thiến biến số kế hoạch thực nghiệm bậc hai Các điểm quy hoạch bậc hai Mức quy hoạch Các nhân tố ảnh hưởng Biến mã hóa - 1,414 -1 +1 + 1,414 = =1,414 𝑁 𝑋 1 0,311 0,32 0,35 0,38 0,389 0,03 0,009 𝐶 𝐶𝐾𝐷 2 1,01 1,1 1,3 1,5 1,59 0,2 0,09 Bảng Cấp phối thực nghiệm theo quy hoạch bậc hai tâm xoay Biến mã STT 10 11 12 13 1 2 +1 -1 +1 -1 +1,414 -1,414 0 0 0 +1 +1 -1 -1 0 +1,414 -1,414 0 0 Biến thực 𝐶 𝑁 𝑋 𝐶𝐾𝐷 0,38 1,5 0,38 1,1 0,32 1,5 0,32 1,1 0,35 1,59 0,35 1,01 0,389 1,3 0,311 1,3 0,35 1,3 0,35 1,3 0,35 1,3 0,35 1,3 0,35 1,3 Cấp phối thực nghiệm (kg/m3) X FA Xi SF90 C N SR5000 500 576 515 597 493 607 532 556 544 544 544 544 544 100 115 103 119 99 121 106 111 109 109 109 109 109 200 230 206 239 197 243 213 222 218 218 218 218 218 50,0 57,6 51,5 59,7 49,3 60,7 53,2 55,6 54,4 54,4 54,4 54,4 54,4 1274 1076 1314 1116 1331 1043 1176 1229 1202 1202 1202 1202 1202 190 219 165 191 172 213 207 173 190 190 190 190 190 7,5 8,6 7,7 9,0 7,4 9,1 8,0 8,3 8,2 8,2 8,2 8,2 8,2 Bảng Giá trị trung bình độ chảy xịe hỗn hợp bê tơng hạt mịn chất lượng cao Biến thực Các biến mã ma trận thực nghiệm Stt N X C CKD 0 10 11 12 13 0,38 0,38 0,32 0,32 0,35 0,35 0,389 0,311 0,35 0,35 0,35 0,35 0,35 1,5 1,1 1,5 1,1 1,59 1,01 1,3 1,3 1,3 1,3 1,3 1,3 1,3 +1 +1 +1 +1 -1 +1 +1 +1 -1 +1 -1 -1 +1 +1,414 +1 -1,414 +1 +1,414 +1 -1,414 +1 0 +1 0 +1 0 +1 0 +1 0 (Y 1j − Y1j ) = 0,529 1 S1d2 = 0,066 2 Độ chảy xịe hỗn hợp bê tơng ,Y1 (cm) 12 12 22 Y1j Y1j (Y1j − Y1j )2 (Y01j − Y01j )2 -1 -1 0 0 0 0 1 1 2 0 0 0 1 1 0 2 0 0 15,8 17,0 18,2 20,0 16,7 18,6 16,0 18,8 18,3 18,2 18,3 18,0 18,2 15,73 17,15 18,07 19,49 16,73 18,74 15,83 19,14 18,2 18,2 18,2 18,2 18,2 0,005 0,023 0,017 0,26 0,001 0,019 0,029 0,115 0,01 0,01 0,04 0,01 0,01 0,04 (Y 01j − Y01j ) = 0,06 S1ll2 =0,015 S1bj = 0,034 Y01 = 18,2 Tăng Văn Lâm nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 62(1), 73-84 3.3 Mô ảnh hưởng hỗn hợp phụ gia biến tính đến tính chất bê tơng chất lượng cao hạt mịn a) Xác định phương trình hàm hồi quy bậc hai Áp dụng cơng thức tốn học theo lý thuyết quy hoạch thực nghiệm bậc hai tâm xoay Box – Wilson sử dụng phần mềm Matlab, thu phương trình hồi quy bậc hai sau: (i) Về độ chảy xòe hỗn hợp bê tông côn mini: Y1 = 18,2 – 0,711 – 1,172 + 0,1512 (4) – 0,21912 – 0,34422 (cm) (ii) Về cường độ nén mẫu bê tông chất lượng cao hạt mịn tuổi 28 ngày: Y2 = 69,28 – 1,9981+ 0,9832 – 2,34712 (5) + 0,21312 – 1,77622 (MPa) Dựa phương trình hàm hồi quy (4) ̂ (5), giá trị tính tốn theo hàm mục tiêu 𝑌̂ 1𝑗 , 𝑌2𝑗 ̂ ̂ 𝑌01 , 𝑌02 xác định theo biến mã 1; 2 trình bày Bảng 10 b) Kiểm tra tính có nghĩa hệ số bj hàm hồi quy theo chuẩn số Student (t (f2)) Hệ số bj hàm hồi quy bậc hai có nghĩa nếu: tbj t(f2) (6) Trong đó: (t (f2)) - giá trị tra bảng chuẩn số Student, với mức độ có nghĩa = 0,05, bậc tự lặp lại f2 = N×(k-1) = 13×(2-1) = 13 (k=2 số yếu tố nghiên cứu) Tra Bảng 3.2 tài liệu (Bolshev Smirnov, 1993) thu t0,05 (13) = 1,7709 Giá trị chuẩn số Student tbj hệ số bj xác định theo công thức (7): 𝑡𝑏𝑗 = |𝑏𝑗| 𝑆𝑏𝑗 =>|𝑏𝑗 | = 𝑡𝑏𝑗 × 𝑆𝑏𝑗 (7) Độ lệch tiêu chuẩn thứ Sbj hệ số thứ bj xác định theo công thức (8): 𝑆𝐼𝐼 𝑆𝑏𝑗 = √ 𝑁 (8) Trong đó: N tổng số điểm thí nghiệm (N = 13) Sll2 phương sai lặp, xác định theo công thức (9) trị đưa Bảng 10 m (Y oj S ll2 = - Yo j =1 ) (9) m -1 Trong đó: m - Số thí nghiệm lặp tâm (m = 5), Yoj - giá trị thu thí nghiệm thứ i tâm kế 81 hoạch Y0 - giá trị trung bình thí nghiệm tâm kế hoạch (Bảng 10) Theo tính có nghĩa bj chuẩn số Student: tbj t(f2) với t bj = Sbj Do đó, hệ số bj phương trình hàm hồi quy coi có nghĩa, nếu: b j t0,05 (13) S bj = 1,7709 S bj Theo giá trị tra bảng chuẩn số Student giá trị Sbj tính tốn liệt kê Bảng 10, phương trình hàm hồi quy bậc hai thu gọn độ chảy xòe hỗn hợp bê tông hạt mịn cường độ nén mẫu bê tông chất lượng cao hạt mịn tuổi 28 ngày thu sau: Y1 = 18,2 – 0,711 – 1,172 (10) – 0,21912 – 0,34422 (cm) Y2 = 69,28 –1,9981 + 0,9832 (11) – 2,34712 – 1,77622 (MPa) c) Kiểm tra tính tương hợp mơ hình theo chuẩn số Fisher Theo lý thuyết quy hoạch thực nghiệm (Nguyễn Như Quý Mai Quế Anh, 2020; Tang Van Lam nnk, 2017, Nguyễn Minh Tuyển, 2007) tính tương hợp mơ hình thực nghiệm kiểm tra theo chuẩn số Fisher xác định theo công thức (12): F0 = Sd2 (12) S ll2 Với Sll2 - phương sai lặp, xác định theo công thức (9) trị đưa Bảng 10 Sd2 - phương sai dư, xác định theo công thức (13) trị xác định Bảng 10 N S d2 = (Y j -Yj j =1 ) (13) N-B Trong đó: Yj; Y j - giá trị thực nghiệm giá trị tính tốn hàm mục tiêu; N - số thí nghiệm (N= 13) B - số hệ số có nghĩa phương trình hồi quy số (10) (11) (B= 5) Bên cạnh Fα (f1; f2) giá trị chuẩn số Fisher, xác định cách tra Bảng 3.5 tài liệu (Bolshev Smirnov, 1993) mức có nghĩa = 0,05; bậc tự lặp f2 = m-1 = 5-1 = 4; bậc tự dư f1 = N – B = 13 - = Từ đó, thu F0,05 (8; 4) = 6,041 82 Tăng Văn Lâm nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 62(1), 73-84 Theo giá trị S2d S2ll thu Bảng 10, xác định giá trị chuẩn số Fisher hàm mục tiêu độ chảy xòe hỗn hợp bê tông cường độ nén mẫu bê tông tuổi 28 ngày lần lươt sau: Với hàm độ chảy xòe hỗn hợp bê tông: F01 = 4,40 Với hàm cường độ nén mẫu bê tông tuổi 28 ngày: F02 = 1,91 Do F01 = 4,40 < F0,05 (8; 4) = 6,041 F02 < F0,05 (8; 4) = 6,041, mơ hình thực nghiệm, mơ tả độ chảy hỗn hợp bê tông hạt mịn cường độ nén mẫu bê tông chất lượng cao hạt mịn, xác định thông qua hàm hồi quy số (10) (11) hoàn toàn tương hợp với tranh thực nghiệm, với kế hoạch thực nghiệm Qua 𝑁 𝐶 cho thấy rằng, tỷ lệ (biến 1) (biến 𝑋 đường đồng mức tương ứng phương trình hồi quy (10) (11) thể Hình e) Xác định giá trị cường độ nén lớn hàm mục tiêu cấp phối tối ưu Trong nghiên cứu đã sử dụng phần mềm toán học Maple Malab để xác định giá trị cực trị phương trình hồi quy (11) Mặt khác, sở cường độ nén tuổi 28 ngày có giá trị lớn tìm cấp phối tối ưu hỗn hợp bê tông hạt mịn chất lượng cao thể Bảng 11 Từ tỷ lệ tương quan điểm cực trị Bảng 11 xác định cấp phối tối ưu hỗn hợp bê tông Tiếp theo, sở cấp phối tối ưu tiến hành thí nghiệm để xác định tính chất lý hỗn hơp bê tông bê tông Kết đo trực tiếp thực nghiệm tính tốn theo phương trình hồi quy (10) (11) biểu thị Bảng 12 𝐶𝐾𝐶 2) có ảnh hưởng đáng kể đến giá trị hàm mục tiêu mơ hình thực nghiệm d) Biểu diễn bề mặt đường đồng mức phương trình hồi quy Bằng cách sử dụng phần mềm Matlab, hình ảnh bề mặt biểu (a) Hàm hồi quy số (10) (b) Hàm hồi quy số (11) Hình Bề mặt biểu đường đồng mức hàm mục tiêu Bảng 11 Giá trị cực đại hàm mục tiêu biến tương ứng Số phương trình Giá trị cực đại hàm hồi quy cường độ nén tuổi 28 ngày 11 Max Y2 = 69,84 MPa Giá trị biến mã Giá trị biến thực 1 2 N X C CKD -0,425 0,277 0,326 1,315 Tăng Văn Lâm nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất (2), - 83 Bảng 12 So sánh thực nghiệm tính tốn theo phương trình hồi quy Khối lượng Độ chảy xịe hỗn hợp Cường độ nén bê tông Cường độ Khối lượng bê tông (mm) (MPa) tuổi 28 ngày thể tích C N kéo uốn thể tích hỗn hợp bê Theo Theo phương Theo Theo phương tuổi 28 bê tông CKD X tông thực trình hồi quy thực trình hồi quy 3) ngày (MPa) (kg/m (kg/m3) nghiệm (10) nghiệm (11) 0,326 1,315 2338,2 18,5 18,094 70,7 69,84 6,28 2319,5 Kết Bảng 12 cho thấy, giá trị độ chảy - Sử dụng chương trình tốn học Matlab xịe hỗn hợp bê tông cường độ nén cực đại Maple biểu diễn bề mặt biểu mẫu bê tơng tính tốn từ hai phương đường đồng mức mơ hình thực nghiệm trực trình hồi quy số (10) (11) phù hợp với quan Đồng thời xác định giá trị lớn kết kiểm tra từ thực nghiệm Độ chảy xịe cường độ nén mẫu bê tơng chất lượng cao cường độ nén cực đại thu từ hàm hồi quy hạt mịn tuổi 28 ngày 69,84 MPa N/X = đạt 97,8% 98,8% so với kết 0,326 C/CKD = 1,315 đo trực tiếp từ thực nghiệm Do đó, sử dụng - Kết tính tốn độ chảy xịe hỗn hợp hàm hồi quy số (10) (11) để đoán bê tông cường độ nén cực đại từ hàm hồi quy tính chất – lý (tính cơng tác, cường độ,…) bê số (10) (11), thu từ phương pháp quy tông chất lượng cao hạt mịn, loại bê hoạch thực nghiệm, đạt 97,8% tông cường độ cao sử dụng chất kết dính khơng xi 98,8% so với kết đo trực tiếp từ thực nghiệm măng,… nghiên cứu chuyên sâu Do sử dụng hàm hồi quy số (10) Việt Nam nước giới (11) để đoán tính chất – lý loại bê tông cường độ cao không dùng xi măng Kết luận nghiên cứu chuyên sâu Việt Nam nước giới Từ kết nghiên cứu thực nghiệm thu - Có thể sử dụng phương pháp quy hoạch phạm vi nghiên cứu, rút số kết thực nghiệm khác để nghiên cứu, mô luận sau: dự đốn tính chất – lý – hóa quan trọng - Trên sở nguồn vật liệu nước như: xi khác bê tông cường độ cao sử dụng chất kết măng Pc lăng PC40 "Hồng Thạch", cát vàng dính khơng xi măng, bê tơng chất lượng cao hạt sông Lô kết hợp với tổ hợp phụ gia biến tính gồm: mịn,… như: thời gian đơng kết hỗn hợp bê silica fume SF-90, tro bay "Vũng Áng", xỉ luyện kim tơng, độ nở thể tích mơi trường biển, tuổi "Hịa Phát" phụ gia siêu dẻo SR 5000F thọ thời gian phá hoại mẫu bê tơng mơi «SilkRoad» chế tạo bê tơng hạt mịn có trường xâm thực độ chảy xịe mini 18,5 cm, cường độ nén tuổi 28 ngày đạt 68,5 MPa cường độ kéo uốn tuổi 28 ngày 6,13 MPa, phù hợp để sử dụng xây dựng cơng trình Việt Nam - Trong phạm vi nghiên cứu phương pháp quy hoạch thực nghiệm bậc hai BoxWilson mô ảnh hưởng hai biến 1 (N/X) 2 (C/CKD) đến độ chảy xịe hỗn hợp bê tơng hạt mịn cường độ nén mẫu bê tông chất lượng cao hạt mịn tuổi 28 ngày Các phương trình quan hệ bậc hai độ chảy xòe cường độ nén tuổi 28 ngày bê tông chất lượng cao hạt mịn hai biến mã hóa 1 2 xác định hàm hồi quy (4), (5), (10) (11) Các hàm hồi quy tương hợp với tranh thực nghiệm theo chuẩn số Fisher Đóng góp tác giả Tác giả Tăng Văn Lâm: lên ý tưởng, xây dựng mơ hình hồi quy, triển khai nội dung viết thảo gốc; tác Vũ Kim Diến, Ngô Xuân Hùng, Bulgakov Boris Igorevich: thực thí nghiệm xử lý kết thực nghiệm; tác Nguyễn Văn Dương, Đỗ Minh Chiến: thu thập liệu chỉnh sửa thảo Tài liệu tham khảo Bazhenov Y.M., (2011) Concrete technology Moscow: Ed ASV., 528 p (tiếng Nga) Bazhenov Y.M., Magdeev U.Kh., Alimov L.A., Voronin V.V., Goldenberg LB (1998) Fine-grained concrete Moscow: Ed ASV., 148 p (tiếng Nga) 84 Tăng Văn Lâm nnk./Tạp chí Khoa học Kỹ thuật Mỏ - Địa chất 62(1), 73-84 Bazhenov Y.M., Demyanova V.S., Kalashnikov V.I., (2006) Modified high-performance concrete Moscow: Ed ASV., 370 p (tiếng Nga) Bolshev L.N., Smirnov N.V., (1993) Mathematical statistics tables Publishing house of Science, Moscow, 416 p (tiếng Nga) Phùng Mạnh Đắc, (2002) Bê tông phun xây dựng mỏ với q trình tăng trưởng ngành than, Cơng nghệ bê tơng phun xây dựng Mỏ cơng trình Ngầm, Hà Nội 2002, 70-75 Đào Viết Đoàn, Tăng Văn Lâm, (2017) Các tham số công nghệ lựa chọn phương pháp kiểm tra trường cường độ bê tông phun chống giữ đường hầm, Tạp chí Cơng nghiệp Mỏ, số 01, 24-28 Ngo Xuan Hung, Tang Van Lam, Bulgakov B.I., Alexandrova O.V., Larsen O.A., Ha Hoa Ki, Melnikov A.I., (2018) Influence of rice husk ash on the properties of hydraulic concrete Scientific and technical journal on construction and architecture "Vestnik MGSU" 2018 Volume 13 Issue (117), P 768-777 (tiếng Nga) Lam Tang Van, Dien Vu Kim, Bulgakov Boris Igorevich, Bazhenova Sophia Ildarovna, (2019a) Mathematical model of effect of the bottom ash and expanded polystyrene spheres on the polystyrene concrete properties Materials Science Forum Vol 974 Pp 312-318 Lam Tang Van, Dien Vu Kim, Tho Vu Dinh, Bulgakov Boris Igorevich, Bazhenova Sophia Ildarovna, Luong Nguyen Tai Nang, (2019b) Effects of high temperature on high performance fine-grained concrete properties International Scientific Conference Energy Management of Municipal Facilities and Sustainable Energy Technologies EMMFT 2018, Pp 660–672 Tang Van Lam, Vu Kim Dien, Bulgakov Boris Igorevich, Bazhenova Sophia Ildarovna, (2019c) Mô ảnh hưởng thành phần đến cường độ bê tông phun hạt mịn thực nghiệm International conference on architecture and civil engineering (ICACE 2019) Hà Nội, Việt Nam 314324 Tr Tăng Văn Lâm, Nguyễn Đình Trinh, Nguyễn Dỗn Tùng Lâm, Bulgakov Boris Igorevich, (2019d) Mô ảnh hưởng hỗn hợp phụ gia khống biến tính đến độ nở sun phát bê tơng dùng cơng trình thủy Tạp chí Khoa học Kỹ thuật Mơi trường Số 66 (9/2019) Tr 71-80 Tang Van Lam, Vu Kim Dien, Ngo Xuan Hung, Nguyen Trong Chuc, Bulgakov B.I., Bazhenova O.Yu., Galtseva N.A., (2019e) Influence of the waterbinding ratio and complex organic-mineral additive on the properties of concrete suitable for the construction of offshore hydraulic structures Journal of Industrial and Civil Engineering No.03, P - 17 (tiếng Nga) Lam Van Tang, Boris Bulgakov, Sofia Bazhenova, Olga Aleksandrova, Anh Ngoc Pham, Tho Dinh Vu, (2018a) Effect of Rice Husk Ash and Fly Ash on the workability of concrete mixture in the HighRise Construction E3S Web of Conferences 33, 02029, 13 Tang Van Lam, Boris Bulgakov, Olga Aleksandrova, Anh Ngoc Pham, Yuri Mikhaylovich, (2018b) Effect of rice husk ash on hydrotechnical concrete behavior IOP Conf Series: Materials Science and Engineering 365, 032007, Tang Van Lam, Bulgakov B.I., Alexandrova O.V., Larsen O.A., Ngo Xuan Hung, Nguyen Duc Vinh Quang, (2017) Determination of deformation of fine-grained concrete in a liquid sulfate medium Journal of Industrial and Civil Engineering No.08, P 82 - 86 (tiếng Nga) Tang Van Lam, Bulgakov B.I., Alexandrova O.V, (2017) Mathematical modeling of the influence of raw materials on the strength of highperformance fine-grained concrete under compression Scientific and technical journal on construction and architecture "Vestnik MGSU" Volume 12 Issue (108), P 999-1009 (tiếng Nga) Nguyễn Quang Phích, (2002) Khả sử dụng bê tông phun xây dựng cơng trình ngầm mỏ, Cơng nghệ bê tơng phun xây dựng Mỏ cơng trình Ngầm, Hà Nội 2002, 40-43 Nguyễn Như Quý, Mai Quế Anh, (2020) Lý thuyết bê tông NXB Xây dựng, Hà Nội, 210 Tr Pham Duc Thang, Bulgakov B.I., Tang Van Lam, (2016) Application of fine-grained shotcrete for the construction of metro tunnels Scientific and technical journal on construction and architecture "Vestnik MGSU" No.07, P 81 - 90 (tiếng Nga) Nguyễn Văn Tuấn, Phạm Hữu Hanh, Nguyễn Công Thắng, Lê Trung Thành, Văn Viết Thiên Ân, Hồng Tuấn Nghĩa, (2018) Bê tơng chất lượng siêu cao NXB Xây dựng, Hà Nội, 300 Tr Nguyễn Minh Tuyển, (2007) Quy hoạch thực nghiệm NXB Khoa hoc Kỹ thuật Hà Nội 264 p Williams H.P., (2013) Model Building in Mathematical Programming John Wiley & Sons ... chuẩn hướng dẫn để tính tốn thiết kế thành phần bê tông hạt mịn sử dụng chất kết dính từ xi măng Pc lăng chất kết dính không xi măng Các nghiên cứu bê tông chất lượng cao hạt mịn chủ yếu dựa... liệu nghiên cứu chuyên sâu tiếp theo, đặc biệt bê tông cường độ cao, chất lượng cao hạt mịn sử dụng chất kết dính khơng xi măng Vật liệu phương pháp nghiên cứu 2.1 Vật liệu sử dụng Vật liệu sử dụng. .. chất lượng cao hạt mịn, loại bê hoạch thực nghiệm, đạt 97,8% tông cường độ cao sử dụng chất kết dính khơng xi 98,8% so với kết đo trực tiếp từ thực nghiệm măng,… nghiên cứu chuyên sâu Do sử dụng