Tạp chí Khoa học Cơng nghệ Xây dựng, NUCE 2021 15 (3V): 69–78 ƯỚC LƯỢNG KHẢ NĂNG CHỊU NÉN ĐÚNG TÂM CỦA CỘT ỐNG THÉP NHỒI BÊ TÔNG BẰNG THUẬT TỐN MÁY HỌC Phạm Thái Hồna,∗ a Khoa Xây dựng dân dụng công nghiệp, Trường Đại học Xây dựng, 55 đường Giải Phóng,quận Hai Bà Trưng, Hà Nội, Việt Nam Nhận ngày 30/05/2021, Sửa xong 22/06/2021, Chấp nhận đăng 23/06/2021 Tóm tắt Cấu kiện cột ống thép nhồi bê tông (Concrete-Filled Steel Tube-CFST) với nhiều ưu điểm quan tâm nghiên cứu ứng dụng rộng rãi cơng trình xây dựng khoảng thập kỷ qua Bên cạnh đó, phát triển mạnh mẽ cùa kỹ thuật máy tính mở nhiều hội ứng dụng kỹ thuật trí tuệ nhân tạo (AI) vào thiết kế cơng trình Bài báo giới thiệu ứng dụng XGBoost, kỹ thuật trí tuệ nhân tạo mạnh mẽ để ước lượng khả chịu nén dọc trục cột ống thép nhồi bê tơng (CFST) tiết diện trịn Kết từ 1017 thí nghiệm loại cấu kiện sử dụng làm liệu để huấn luyện mơ hình Các thơng số bao gồm đường kính ngồi tiết diện (D), chiều dày ống thép (t), cường độ chảy thép ( fy ), cường độ chịu nén bê tông ( fc0 ) chiều dài cột (L) sử dụng thông số đầu vào, khả chịu lực nén tâm dọc trục (N) thông số đầu Kết huấn luyện từ mơ hình cho thấy thuật tốn XGBoost ứng dụng hiệu việc ước lượng khả chịu nén dọc trục cột CFST tiết diện trịn với độ xác cao với thời gian ngắn Từ khố: XGBoost; trí tuệ nhân tạo; máy học; ống thép nhồi bê tông; CFST ESTIMATION OF CONCENTRIC AXIAL COMPRESSIVE LOAD-CARRYING CAPACITY OF CONCRETEFILLED STEEL TUBE USING MACHINE LEARNING ALGORITHMS Abstract Concrete-Filled Steel Tube (Concrete-Filled Steel Tube-CFST) columns with many advantages have been studied and widely applied in construction projects for decades Besides, the rapid development of computer engineering has opened up many opportunities to apply artificial intelligence (AI) techniques in building design This paper introduces and applies XGBoost, one of the newest and most powerful artificial intelligence techniques to estimate the concentric axial compressive load-carrying capacity of circular section concrete-filled steel tube (CFST) columns The results from 1017 experiments on this type of columns are used as dataset to train the model The outer diameter of the section (D), the thickness of the steel tube (t), the yield strength of the steel ( fy ), the compressive strength of the concrete ( fc0 ), and the length of the column (L) are used as input parameters, while the axial compression force capacity (N) is the output parameter The training results from the model showed that the XGBoost algorithm can be effectively applied in estimating the axial compressive capacity of the circular section CFST columns with high accuracy and short time Keywords: XGBoost; artificial intelligence; machine learning; concrete-filled steel tube; CFST https://doi.org/10.31814/stce.nuce2021-15(3V)-06 © 2021 Trường Đại học Xây dựng (NUCE) Giới thiệu Cột ống thép nhồi bê tông (concrete-filled steel tube) (CFST) kết cấu bao gồm ống thép lấp đầy bên bê tông Ống thép lúc đóng vai trị vừa ván khn cho trình đổ bê ∗ Tác giả đại diện Địa e-mail: hoanpt@nuce.edu.vn (Hoàn, P T.) 69 Hoàn, P T / Tạp chí Khoa học Cơng nghệ Xây dựng tơng, vừa tham gia tăng cường khả chịu tải bê tơng Trong đó, việc nhồi bê tơng giúp tăng ổn định ống thép độ dẻo tiết diện Ngồi ra, cột CFST cịn có khả hấp thụ lượng cao, độ bền cao khả chống cháy tốt [1, 2] Vì vậy, ứng xử cột CFST quan tâm nghiên cứu rộng rãi suốt ba thập kỷ qua Bài báo tập trung nghiên cứu cột CFST tiết diện trịn dạng tiết diện có độ cứng độ bền chảy tốt so với dạng tiết diện chữ nhật hình vng nghiên cứu trước [3–6] Trong suốt nhiều thập kỷ qua, nhiều thí nghiệm cọc tròn CFST thực với nhiều trọng tâm nghiên cứu khác như: cường độ bê tông [7], tỉ lệ đường kính độ dày ống thép [8] hay ảnh hưởng liên kết bê tơng ống thép [9], v.v Một số chương trình thí nghiệm lớn cọc CFST kể đến sau Chương trình thí nghiệm Knowles Park [10] thực vào cuối năm 1960 nhằm đánh giá ứng xử cột CFST chịu tải tâm tải lệch tâm Các thí nghiệm Liu Goel [11] cuối năm 1980 nhằm đánh giá ứng xử cọc CFST chịu tải trọng động theo chu kỳ Kilpatrick Rangan [12] nghiên cứu ảnh hưởng việc sử dụng bê tông cường độ cao cột CFST Sakino cs [13] nghiên cứu ảnh hưởng hình dạng cường độ ống thép, tỷ lệ đường kính ống chiều dày cường độ bê tơng Tuy nhiên, số lượng thí nghiệm thực nhiều hạn chế chương trình thí nghiệm thường địi hỏi tiêu tốn nhiều tiền bạc, hệ thống thiết bị thử nghiệm đắt tiền, tốn nhiều thời gian công sức người Dựa kết thực nghiệm, công thức thiết kế cọc CFST đề xuất Một cách truyền thống, công thức thiết kế xây dựng dựa phương pháp lý thuyết mô dựa ứng xử làm việc chung loại vật liệu thép bê tông lấy kết thí nghiệm để đánh giá độ xác cơng thức mơ hình đề xuất Tuy nhiên, sai số trình thí nghiệm, sai khác chương trình thực nghiệm ảnh hưởng yếu tố thời tiết, vật liệu đầu vào hay tay nghề người thực khiến cho công thức xây dựng nhà khoa học có khác đáng kể Ví dụ, cơng thức xác định khả chịu tải cột CFST cung cấp tiêu chuẩn Eurocode [14], ACI 318 [15] AISC 360 [16], nhiên, khơng có cơng thức số có tương đồng với Các mơ hình xây dựng dựa mối quan hệ ứng suất-biến dạng giả định trước ống thép bê tơng lịng nên tính xác mơ hình dấu hỏi lớn Gần đây, với phát triển mạnh mẽ khoa học máy tính, kỹ thuật máy học (machine learning-ML) ứng dụng rộng rãi lĩnh vực khoa học sống Các kỹ thuật ML hiểu cách đơn giản phương pháp cho phép xây dựng mơ hình tốn học phức tạp có độ xác cao nhằm thể mối quan hệ thông số đầu vào thông số đầu tập liệu cho trước Với quan điểm này, ML nhiều nhà khoa học sử dụng để xác định ứng xử cơng trình [17–20] Ứng dụng ML ước lượng ứng xử cột CFST nghiên cứu rộng rãi [21, 22] Các nghiên cứu cho thấy có liệu đủ lớn kỹ thuật ML hồn tồn sử dụng để xây dựng nên mơ hình hiệu để ước lượng ứng xử cơng trình nhằm thay hỗ trợ việc xác hóa cơng thức thiết kế truyền thống Một số thuật toán ML mạnh mẽ kể đến là: Gradient Tree Boosting (GTB), học sâu (Deep Learning) (DL), Support Vector Machine (SVM), rừng ngẫu nhiên (Random Forest) (RF), định (Decision Tree) (DT), v.v Hiệu ứng dụng thuật tốn ML thiết kế cơng trình phụ thuộc nhiều yếu tố việc lựa chọn tham số hệ thống phù hợp cho lớp toán cụ thể quan trọng Mỗi lớp toán cần thực hiện, phân tích để tìm cách sử dụng tham số cho mơ hình huấn luyện hiệu tối ưu mặt thời gian Trên sở phân tích trên, nghiên cứu này, tác giả trình bày chi tiết cách xây dựng mơ hình huấn luyện cho tốn ước lượng khả chịu tải cột CFST sử dụng thuật toán Extreme 70 Hồn, P T / Tạp chí Khoa học Cơng nghệ Xây dựng Gradient Boosting (XGBoost) XGBoost thuật toán ML Chen Guestin [23] đề xuất năm 2016 tỏ hiệu việc giải tốn phân tích hồi qui phân loại [23–25] Một tập liệu kết thí nghiệm bao gồm 1017 mẫu cột CFST tròn sử dụng để huấn luyện kiểm tra chương trình xây dựng Ba mơ hình ML khác bao gồm SVM, RF Adaboost sử dụng để chứng minh tính hiệu chương trình đề xuất Thuật tốn XGBoost XGBoost Chen Guestin [23] đề xuất vào năm 2016 dựa ngun lý chung chí Khoa học Cơng nghệ Xây dựng, NUCE2020 p-ISSN 2615-9058; e-ISSN 2734-9489 thuật toán câyTạp GTB sở kết hợp nhiều định có độ sai số cao (cịn gọi sở hay học yếu) thành định tốt Nguyên lý mơ tả Hình Đầu tiên, định tạo cách ngẫu nhiên Tương ứng với 114 liệu xác, yi , gọi z1i Các giá trị z1i sai số mơ hình huấn liệu đầu vào Xi , ước lượng liệu đầu y1i Lúc sai số liệu ước lượng được, luyện Để giảm sai số này, ra1 định thứ tạo Khác với y1i , 115 liệu xác, yi , gọi zi Các giá trị z1i sai số mơ hình huấn luyện 116 định thứ nhất, liệu đầu cho thứ giá trị Thơng đầuthứ nhất, liệu Để giảm sai số này, định thứ tạo Khác với câyzira số định đầu117 chocây câythứ thứ2 2gọi zi Thơng số zđầu 2thứ gọi yi độ yi2làvàcác độgiá sai trị lệch với giá sai trị lệch i gọi zi Hiển nhiên 2 với zi 2gọi zi Hiển nhiên giá trị z bé z hay nói cách khác sai số mô 118 zi bé z1i hay nói cách khác sai sối mơ hìnhi huấn luyện lúc hình huấn luyện lúc giảm Cứ vậy, với việc tạo nhiều định, Cứ nhưluyện vậy, với hay tuầnđộ tự nhiều định, tăng sai sốlên Giá mơ hình sai 119 số củagiảm mơ hình huấn việc giảmtạo dần xácracủa trị ước lượng 120 huấn luyện giảm dần hay độ xác tăng lên Giá trị ước lượng cuối cuối mơ hình huấn luyện GTB là: 121 122 mơ hình huấn luyện GTB là: yi' m j yi j y0i = m X j yi j=1 (1) m 123đó m số câyđịnh định sử dụng sốđócây sử dụng 124 125 Hình1.1.Sơ Sơ đồ đồ minh minh họa trình huấn Hình huấn luyện luyệncủa củathuật thuậttốn tốnGTB GTB 126 Q trình huấn luyện GTB q trình tối thiểu hóa tổng sai số giá 127 trị ước lượng cuối yi' yi ( i 1, 2, , N ) thể qua hàm mát L yi , yi' 128 129 130 Có nhiều dạng hàm mát sử dụng, phổ biến hàm sai số tồn phương trung bình (mean squared error) (MSE) viết dạng sau: 71 (1) Hồn, P T / Tạp chí Khoa học Cơng nghệ Xây dựng Q trình huấn luyện GTB q trình tối thiểu hóa tổng sai số giá trị ước lượng cuối yi yi (i = 1, 2, , N) thể qua hàm mát L yi , y0i Có nhiều dạng hàm mát sử dụng, phổ biến hàm sai số tồn phương trung bình (mean squared error) (MSE) viết dạng sau: MS E = N X 2 yi − y0i (2) i=1 Để nâng cao hiệu thuật toán GTB, đặc biệt tránh tượng mơ hình q khớp, thuật tốn XGBoost hàm quy hóa (regularization) Ω () thêm vào hàm mát vòng lặp thứ t thể sau: J (t) = N t X X L yi , yi0(t) + Ω ( fn ) i=1 (3) n=1 hàm fn chưa thông tin cấu trúc định điểm Lúc này, thông qua trình tối ưu hàm mục tiêu J, cấu trúc định tối ưu Tuy nhiên, việc huấn luyện định toán phức tạp nhiều so với tốn tối ưu thơng thường khó để kiểm soát tất lúc Do vậy, giải pháp bổ sung thường áp dụng học giữ cố định lại thêm thời điểm Lúc này, thơng số đầu yi0(t) xác định sau: yi (0) = 0 yi (1) = f1 (Xi ) = yi (0) + f1 (Xi ) 0 yi (2) = f1 (Xi ) + f2 (Xi ) = yi (1) + f2 (Xi ) yi (t) = t X (4) fn (Xi ) = yi (t−1) + ft (Xi ) n=1 Nếu xem hàm tối ưu hàm MSE phương trình (3) lúc viết thành: J (t) N t X X (t) = L yi , yi + Ω ( fn ) i=1 n=1 N X L yi , yi (t−1) + ft (Xi ) + Ω ( ft ) + constant = (5) i=1 Vế phải phương trình (5) sử dụng biến đổi Taylor giữ đến số hạng bậc Lúc ta có: J đó: (t) # N " X (t−1) = L yi , yi + gi ft (Xi ) + hi ft (Xi ) + Ω ( ft ) + constant i=1 gi = ∂y0 (t−1) L yi , yi (t−1) i hi = ∂ y (t−1) L yi , yi (t−1) i 72 (6) (7) Hoàn, P T / Tạp chí Khoa học Cơng nghệ Xây dựng Loại bỏ số công thức (6), trình tối ưu tối thiểu hóa hàm: J 0(t) # N " X = gi ft (Xi ) + hi ft (Xi ) + Ω ( ft ) i=1 (8) Công thức (8) cho thấy hàm mục tiêu lúc phụ thuộc vào gi hi Do đó, XGBoost tùy chọn hàm mát hiệu áp dụng cho nhiều dạng hàm mát với cách giải xem gi hi biến đầu vào Hàm cấu trúc fn định nghĩa sau: ft (X) = wq(X) (9) w véc-tơ điểm lá, q hàm gán cho liệu tương ứng Hàm Ω ( fn ) chọn sau [20]: T X Ω ( fn ) = γT + λ w j=1 j (10) γ thể cho độ phức tạp tổ chức định; T số lượng định; λ hệ số phóng đại hàm phạt; w véc-tơ điểm Lúc hàm mục tiêu viết thành: T Xây dựng, NUCE2020 p-ISSN 2615-9058; e-ISSN 2734-9489 X học Cơng nghệ X X Tạp chí Khoa 0(t) J ≈ g w + h + λ w j i i j + γT i∈I j i∈I j i=1 (11)thấy tài l 177 cs [26] thu thập " # sử dụng Chi tiết liệu xem tìm T X 1Các liệu bao2 gồm thông số: D, t, cường độ chảy thép (fy), cường độ c = 178 G j w[26] H j + λ w j + γT j+ 179 nén bê tông (fc’), chiều dài cột (L) khả chịu nén tâm cột (N) i=1 180 phân bố thông số thể Hình Các giá trị lớn (GTL 181trong giá trịthứ nhỏ j (GTNN) giá trịcực trung bình mát thơng I j đại diện cho tất mẫu Bây bàivàtoán tiểu của(GTTB) hàm đượcsố Kết thấy phần cột CFST thí nghiệm với nút cường độ thép thô chuyển thành tốn tìm giá trị nhỏ182 hàm bậccho hai Nói cáchlớn khác, sau phân chia 183 thường (trên 90% mẫu có f 460 MPa 75% mẫu có f ’ ≤ 50 MPa) Cùng y c định, hiệu suất mơ hình cải thiện phân chia chấp 184 phần lớn mẫu dạng tiết diện khơng chắc, có tiết nhận Nếu khơng, q trình chia tách nút dừng lại Nhờ đó, phân chia nút định diện thuộc mảnh Sự tương quan thơng số thể Hình Kết cho th dừng sớm giúp tăng hiệu suất185 trình huấn luyện tránh tượng mơ hình q khớp 186 có mối tương quan lớn N D N t Xây dựng tập liệu huấn luyện Hình thể tiết diện ngang điển hình cột CFST tiết diện trịn Trong đó, đường kính tiết diện chiều dày ống thép ký hiệu D t Tập liệu gồm 1017 kết thí nghiệm cột CFST tiết diện tròn chịu nén tâm Thai cs [26] thu thập sử dụng Chi tiết liệu xem tìm thấy tài liệu 187 [26] Các liệu bao gồm thông số: D, t, cường 188 độ chảy thép ( fy ), cường độ chịu nén bê tông ( fc ), chiều dài cột (L) khả chịu nén 73 Minh họa ống tiết thép diệnnhồi ống bê thép Hình Hình Minh họa tiết diện tơng CFST nhồi bê tơng CFST Hồn, P T / Tạp chí Khoa học Công nghệ Xây dựng tâm cột (N) Sự phân bố thông số thể Hình Các giá trị lớn (GTLN), giá trị nhỏ (GTNN) giá trị trung bình (GTTB) thông số thể Kết cho thấy phần lớn cột CFST thí nghiệm với cường độ thép thông thường (trên 90% mẫu có fy 460 MPa 75% mẫu có fc0 ≤ 50 MPa) Cùng với phần lớn mẫu 187 187dạng tiết diện không chắc, có tiết diện thuộc dạng mảnh Sự tương quan thông Minh họadiện tiết ống diệnthép ống nhồi thép bê nhồi bêCFST tông CFST 188 188 HìnhHình Minh họa tiết tơng 189 189 190 190 Tạp chí Khoa học Cơng nghệ Xây dựng, NUCE2020 p-ISSN 2615-9058; e-ISSN 2734-9489 Tạp chí Khoa học Cơng nghệ Xây dựng, NUCE2020 p-ISSN 2615-9058; e-ISSN 2734-9489 Tạp chí Khoa học Công nghệ Xây dựng, NUCE2020 p-ISSN 2615-9058; e-ISSN 2734-9489 Tạp chí Khoa học Cơng nghệ Xây dựng, NUCE2020 p-ISSN 2615-9058; e-ISSN 2734-9489 kính ngồi củacủa cột (D) a) (a) Đường kính ngồi cột a)Đường Đường kính ngồi của(D) cột (D) Chiềucủa dày ống thép(t) (t) b) Chiều dày thép b)(b)Chiều dàyống ống thép (t) 7 191 191 191 191 192 192 192 192 193 193 193 193 194 194 194 194 195 195 195 195 (c) Chiều dàidài cột (L) c) Chiều cột (L) c)c)c) Chiều dài cột (L) Chiều dài cột (L) Chiều dài cột (L) (d) Cường độ chảy thép ( fy ) (f (f d) Cường độ chảy thép ) y) d)d) Cường độđộ chảy thép y) yy(f chảy thép d)Cường Cường độ chảy thép (e) Cườngđộ độ chịu chịu nén củacủa bê tông (tông fc ) (f) Khả chịunén nén tâm của cột (N) Cường nén bê Khả chịu tâm (N) c’) f) f) e) Cường độ chịu nén bê tông (f(fc(f ’) Khả chịu nén tâm cộtcột (N) e)e)e) Cường độ chịu nén bê tông ’) f) Khả chịu nén tâm (N) c Cường độ chịu nén bê tông (fc’) f) Khả chịu nén tâm củacộtcột (N) Hình Phân bố thơng số cột CFST Hình Phân bố thơng số cột CFST Hình Phân bố thơng sốsốcột CFST Hình Hình3.3.Phân Phânbốbốcủa củacác cácthơng thơng sốcột cộtCFST CFST 74 193 Hồn, P T / Tạp chí Khoa học Cơng nghệ Xây dựng e) Cường độ chịu nén bê tông (fc’) f) Khả chịu nén tâm cột (N) 194 số thể Hình Kết cho thấy có mối tương quan lớn N D 195 Hình Phân bố thơng số cột CFST N t 196 197 Hình 4 Tương cácthơng thơngsốsố khảo Hình Tươngquan quan giữa khảo sát sát Trong mơ hình huấn luyện, thông số đầu vào lựa chọn thơng số: D, Trong mơ hình huấn luyện, thơng số đầu vào lựa chọn thông số: D, t, fy , fc0 L t, fy, fc’ L Thơng số đầu có thơng số N Từ tập liệu ban đầu, tập 199 Thơng số đầu có thông số N Từ tập liệu ban đầu, tập liệu riêng biệt xây dựng 200 riêng biệt xây dựng tập huấnđược luyện tập kiểm tra Các liệu tậpliệu huấn luyện vàđược tập kiểm tra Các liệu chuẩn hóa khoảng [0, 1] dựachuẩn công thức sau: xi xi0 =8 (12) [xi ] 198 [xi ] giá trị lớn biến đầu vào xi Dạng mơ hình sử dụng nghiên cứu mơ hình hồi quy Ngơn ngữ lập trình Python thư viện Sklearn sử dụng để xây dựng chương trình huấn luyện Một số thơng số mơ hình huấn luyện sử dụng thuật tốn XGBoost lấy sau: số lượng định n_estimators = 10000, dạng định ‘gbtree’, tỉ lệ học lấy learning_rate = 0,01, dạng hàm mát MSE, γ = λ = Các tham số khác lấy theo giá trị hệ thống cho sẵn Lưu ý rằng, giá trị xác định dựa việc thử nhiều lần nhằm tìm kết tốt cho mơ hình huấn luyện Kết Trong phần này, hiệu suất XGBoost so sánh với thuật toán ML phổ biến khác SVM, RF Adaboost Để tăng tính khách quan kết tính tốn thuật toán, kỹ thuật Kfold sử dụng với 10 tập khác Điều có nghĩa tập liệu chia thành 10 tập con, chạy 10 lần lần chạy lấy tập tạo thành tập huấn luyện tập lại tập kiểm tra mơ hình Kết cuối giá trị trung bình 10 lần chạy Để giảm thiểu ảnh hưởng tính ngẫu nhiên q trình huấn luyện, thuật toán chạy 10 lần độc lập Kết so sánh thuật toán đánh giá dựa giá trị MSE, hệ số R2 hệ số R2 hiệu chỉnh (R2 ) Kết mơ hình huấn luyện tập kiểm tra thể Bảng Kết cho thấy giá trị MSE trung bình sử dụng XGBoost 7,08E-5, nhỏ nhiều so với kết thu sử dụng SVM, RF Adaboost Mơ hình huấn luyện sử dụng 75 ... chảy thép (e) Cườngđộ độ chịu chịu nén củacủa bê tông (tông fc ) (f) Khả chịunén nén tâm của cột (N) Cường nén bê Khả chịu tâm (N) c’) f) f) e) Cường độ chịu nén bê tông (f(fc(f ’) Khả chịu nén tâm. .. Khả chịu nén tâm cộtcột (N) e)e)e) Cường độ chịu nén bê tông ’) f) Khả chịu nén tâm (N) c Cường độ chịu nén bê tông (fc’) f) Khả chịu nén tâm củacộtcột (N) Hình Phân bố thơng số cột CFST Hình Phân... 188 độ chảy thép ( fy ), cường độ chịu nén bê tông ( fc ), chiều dài cột (L) khả chịu nén 73 Minh họa ống tiết thép diệnnhồi ống bê thép Hình Hình Minh họa tiết diện tơng CFST nhồi bê tơng CFST