Tổng quan nghiên cứu

Trong xu thế phát triển công trình hạ tầng ven biển và hải đảo, việc ứng dụng vật liệu composite nhằm thay thế các kết cấu kim loại và bê tông cốt thép truyền thống đang trở thành giải pháp đột phá. Các loại vật liệu mái lợp thông thường như tôn kim loại hay fibro xi măng thường bị xuống cấp nhanh chóng, với tỷ lệ suy giảm tuổi thọ lên tới 40% đến 50% sau 3 đến 5 năm vận hành do hiện tượng ăn mòn muối mặn và axít. Nghiên cứu của học viên Nguyễn Hải Bằng dưới sự hướng dẫn của Tiến sĩ Vũ Khắc Bảy tại Trường Đại học Lâm nghiệp đã giải quyết bài toán cơ học then chốt: tính toán độ bền và độ ổn định đàn hồi của kết cấu mái vòm vỏ trụ lượn sóng làm hoàn toàn bằng vật liệu polymer composite nhiều lớp chịu tác động trực tiếp của tải trọng gió bão.

Mục tiêu trọng tâm của luận văn là thiết lập mô hình giải tích - số chính xác nhằm xác định mối quan hệ tương hỗ giữa các thông số hình học gồm chiều dài vòm 17 mét, bán kính cong 6,375 mét, chu kỳ sóng 56 centimet, biên độ sóng 37,5 centimet cùng độ dày thành vỏ 2 centimet đối với khả năng chịu tải trọng gió bão cấp 12 với vận tốc gió vượt trên 33 mét/giây theo tiêu chuẩn Việt Nam. Kết quả nghiên cứu không chỉ chứng minh tính khả thi về mặt chịu lực của kết cấu vỏ mỏng phi kim loại mà còn mang lại giá trị kinh tế to lớn, giúp giảm khoảng 35% chi phí bảo trì định kỳ và kéo dài tuổi thọ công trình lên trên 20 năm tại các vùng duyên hải khắc nghiệt.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn xây dựng trên nền tảng ba hệ thống lý thuyết cơ học kết cấu vỏ mỏng và vật liệu dị hướng hiện đại:

  1. Lý thuyết đàn hồi cổ điển và trạng thái ứng suất phẳng suy rộng: Thiết lập hệ phương trình vi phân cân bằng tĩnh học kết hợp định luật Hooke tổng quát cho vật liệu trực hướng và các hệ thức hình học Cauchy, mô tả chính xác trường chuyển vị và biến dạng của mặt trung bình vỏ trụ.
  2. Lý thuyết mô đun hiệu quả cho composite nhiều lớp: Cho phép đồng hóa cấu trúc composite phân lớp phức tạp thành vật liệu đồng nhất tương đương thông qua việc xác định các hằng số kỹ thuật trung bình gồm mô đun đàn hồi dọc $E_1$, mô đun đàn hồi ngang $E_2$ đạt mức 250 Megapascal, mô đun trượt $G$ đạt 30 Megapascal và hệ số Poisson 0,15.
  3. Phương pháp tương đương Seydel: Quy đổi hình học vỏ trụ lượn sóng biên dạng hình sin về mô hình vỏ trụ trơn tương đương có các độ cứng uốn và độ cứng chống xoắn trực hướng biến đổi, loại bỏ việc phải chia lưới vi mô phức tạp dọc theo từng bước sóng.
  4. Tiêu chuẩn tải trọng gió TCVN 2737 – 1995: Xác định chính xác áp lực gió tác động cục bộ lên từng vùng bề mặt mái vòm với hệ số khí động $C = +0,6$ tại vùng đón gió, $C = -0,8$ tại đỉnh vòm chịu lực hút tĩnh và $C = -0,5$ tại vùng khuất gió, kết hợp hệ số tin cậy tải trọng 1,2 và hệ số gia tăng áp lực theo độ cao từ 1,0 đến 1,18.

Phương pháp nghiên cứu

Phương pháp nghiên cứu chủ đạo là phương pháp giải tích kết hợp biến phân số:

  • Mô hình hóa toán học: Đưa vào hàm ứng suất và hàm độ võng nhằm thỏa mãn phương trình cân bằng liên tục và phương trình tương thích biến dạng của vỏ mỏng. Hệ phương trình vi phân đạo hàm riêng bậc bốn được thiết lập chặt chẽ cùng các điều kiện biên ngàm cứng dọc hai đường sinh đáy và tự do hoàn toàn tại hai đầu vòm cong.
  • Kỹ thuật xấp xỉ nghiệm: Lựa chọn nghiệm giải tích dưới dạng chuỗi hàm trực giao kết hợp giữa đa thức bậc cao và hàm lượng giác với số lượng số hạng chuỗi được chọn lọc là 16 số hạng tương ứng bậc khai triển $M = 4$ và $N = 4$. Lý do lựa chọn phương pháp Ritz kết hợp Bubnov - Galerkin là khả năng chuyển hóa chính xác hệ phương trình vi phân phi tuyến phức tạp thành hệ 32 phương trình đại số tuyến tính, giúp nghiệm hội tụ nhanh và triệt tiêu sai số biên.
  • Rời rạc hóa lưới tính toán: Toàn bộ mái vòm được phân chia thành lưới tọa độ gồm 20 điểm dọc theo đường sinh và 20 điểm theo hướng vòng, tạo thành ma trận 400 nút khảo sát cục bộ. Toàn bộ thuật toán được lập trình tự động hóa trên nền tảng phần mềm Visual Basic chuyên dụng.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình tính toán trên mô hình số với các tham số vật liệu và hình học thực tế đã mang lại những phát hiện quan trọng:

  1. Ứng xử biến dạng và độ võng cực đại: Độ võng của mái vòm đạt giá trị lớn nhất tại vị trí nút tọa độ (10, 13), tức điểm nằm giữa chiều dài vòm và thiên về phía sau đỉnh đón gió. Dưới áp lực gió bão cấp 12, tỷ số độ võng trên chiều dày thành vỏ luôn duy trì ở mức an toàn, không xuất hiện hiện tượng mất ổn định cục bộ.
  2. Phân bố nội lực và mô men uốn: Mô men uốn $M_y$ tập trung cực đại tại mép biên sau hướng gió tại vị trí nút tọa độ (10, 20) với giá trị đạt xấp xỉ 18 kG.cm/cm. Điều này phản ánh rõ nét tác động nguy hiểm của lực hút khí động học tại vùng khuất gió.
  3. Trạng thái ứng suất tổng thể: Cường độ ứng suất tổng hợp $\sigma_u$ đạt giá trị lớn nhất tại nút tọa độ (7, 13), hoàn toàn trùng khớp với điểm đạt cực đại của ứng suất pháp vòng $\sigma_y$. Giá trị ứng suất này thấp hơn 65% so với giới hạn bền kéo cho phép 240 Megapascal và giới hạn bền uốn 290 Megapascal của vật liệu composite cốt sợi.

Thảo luận kết quả

Khi xem xét dữ liệu qua các bảng tra ma trận $20 \times 20$ và biểu đồ phân bố mô men uốn $M_y$, có thể thấy việc tạo hình lượn sóng dọc đường sinh giúp gia tăng độ cứng uốn của kết cấu lên gấp 3,5 lần so với tấm vỏ phẳng có cùng khối lượng vật liệu. Hiện tượng ứng suất tập trung mạnh tại tọa độ (7, 13) và (10, 20) được giải thích do sự chuyển tiếp đột ngột giữa vùng chịu áp lực nén đẩy ($C = +0,6$) sang vùng chịu lực hút chân không cực đại ($C = -0,8$) tại đỉnh vòm.

So sánh với các kết cấu mái vòm kim loại truyền thống, vòm composite lượn sóng dày 2 centimet thể hiện khả năng giảm chấn khí động học vượt trội, hạn chế tối đa nguy cơ mỏi vật liệu khi chịu tải trọng gió lặp lại. Cơ chế làm việc hoàn toàn trong miền đàn hồi tuyến tính khẳng định tính đúng đắn của phương pháp xấp xỉ chuỗi Ritz - Galerkin trong thiết kế kết cấu vỏ mỏng phi truyền thống.

Đề xuất và khuyến nghị

Dựa trên các kết luận khoa học đã đạt được, luận văn đưa ra 4 khuyến nghị then chốt nhằm ứng dụng hiệu quả vào thực tiễn sản xuất và xây dựng:

  1. Chuẩn hóa kích thước hình học trong chế tạo: Các nhà máy sản xuất cấu kiện composite cần thiết lập quy trình đúc khuôn định hình mái vòm với chiều dày tối thiểu 2 centimet, duy trì bước sóng 56 centimet và chiều cao sóng 37,5 centimet nhằm đảm bảo mô đun kháng uốn tối ưu trong thời gian triển khai 6 đến 12 tháng tới.
  2. Gia cường cục bộ tại các vị trí xung yếu: Đơn vị thi công cần tăng cường bổ sung từ 2 đến 3 lớp vải gia cường (tăng khoảng 15% đến 20% lượng sợi cốt) tại dải liên kết chân ngàm và khu vực đỉnh vòm tương ứng dải tọa độ từ hàng 7 đến hàng 10 để nâng hệ số an toàn chịu bão lên mức 1,8.
  3. Ứng dụng phần mềm tính toán chuyên dụng: Khuyến nghị các viện thiết kế xây dựng và cơ điện nông lâm nghiệp áp dụng gói phần mềm Visual Basic đã phát triển để tự động hóa khâu tính toán tải trọng gió, giúp rút ngắn 70% thời gian lập dự toán và thiết kế kỹ thuật.
  4. Xây dựng quy trình kiểm định định kỳ: Chủ đầu tư các công trình kho bãi, nhà xưởng ven biển cần ban hành quy chế kiểm tra liên kết bulông chân vòm định kỳ 12 tháng một lần, đo đạc độ võng dư sau mỗi mùa mưa bão nhằm đảm bảo an toàn tuyệt đối cho người và thiết bị.

Đối tượng nên tham khảo luận văn

Tài liệu luận văn mang lại giá trị tham khảo chuyên sâu và thực tế cho 4 nhóm đối tượng chính:

  • Kỹ sư thiết kế kết cấu công trình biển: Tiếp cận phương pháp quy đổi Seydel và bảng tra mô men uốn để tính toán nhanh kết cấu mái che, nhà xưởng ven biển chịu tải trọng gió bão cấp 12 mà không cần lập mô hình phần tử hữu hạn phức tạp.
  • Doanh nghiệp sản xuất vật liệu và cấu kiện Composite: Nắm bắt chính xác các thông số cơ lý như độ bền kéo 240 Megapascal, độ bền nén 190 Megapascal cùng quy cách hình học profin sóng chuẩn để chế tạo khuôn mẫu thương mại hóa.
  • Nghiên cứu sinh và giảng viên ngành Cơ học ứng dụng: Sử dụng hệ 32 phương trình Bubnov - Galerkin và chuỗi hàm chuyển vị trực giao làm tài liệu giảng dạy mẫu mực về cơ học vật liệu composite và lý thuyết vỏ mỏng.
  • Cơ quan quản lý và tư vấn giám sát xây dựng: Vận dụng khung tiêu chuẩn TCVN 2737 – 1995 kết hợp các cảnh báo ứng suất nguy hiểm để xây dựng quy trình thẩm tra, nghiệm thu các công trình dân dụng và quốc phòng phi kim loại.

Câu hỏi thường gặp

Tại sao vật liệu composite lại vượt trội hơn tôn thép mạ kẽm tại vùng duyên hải? Composite nền polymer cốt sợi có khả năng kháng axít và chống ăn mòn muối mặn 100%, đồng thời trọng lượng nhẹ hơn kết cấu bê tông khoảng 60% đến 70%, giúp giảm tải trọng móng và triệt tiêu hoàn toàn chi phí sơn phủ bảo vệ định kỳ hàng năm.

Phương pháp Seydel giải quyết bài toán vỏ lượn sóng như thế nào? Phương pháp Seydel đồng hóa hình học lượn sóng phức tạp thành vỏ trụ trơn tương đương có các mô đun đàn hồi trực hướng mới ($D_1, D_2, D_3$), giúp đơn giản hóa hệ phương trình vi phân và giảm trên 80% khối lượng tính toán ma trận số.

Khu vực nào trên mái vòm chịu rủi ro biến dạng cao nhất khi có gió lớn? Theo kết quả mô phỏng trên 400 nút lưới, điểm giữa vòm lệch về phía sau hướng gió tại tọa độ (10, 13) chịu độ võng lớn nhất, trong khi mép ngàm chân vòm tại tọa độ (10, 20) chịu mô men uốn $M_y$ cực đại lên tới 18 kG.cm/cm.

Mái vòm composite dày 2 centimet có an toàn trước bão cấp 12 không? Hoàn toàn an toàn. Tính toán thực tế chỉ ra rằng cường độ ứng suất tổng hợp cực đại tại nút (7, 13) vẫn thấp hơn nhiều so với độ bền kéo giới hạn 240 Megapascal và độ bền uốn 290 Megapascal của vật liệu.

Thuật toán trong luận văn có thể nâng cấp lên các phần mềm CAD/CAE hiện đại không? Mã nguồn phát triển trên Visual Basic có cấu trúc module hóa cao, hoàn toàn có khả năng chuyển đổi dữ liệu hai chiều sang Excel, Word, PDF hoặc tích hợp API trực tiếp vào các phần mềm mô phỏng như ANSYS, SAP2000 và AutoCAD.

Kết luận

  • Đã xây dựng thành công mô hình cơ học giải tích - số tính toán độ bền đàn hồi cho kết cấu mái vòm vỏ trụ composite lượn sóng chịu tải trọng gió bão theo tiêu chuẩn TCVN 2737 – 1995.
  • Đề xuất biểu thức giải tích chuỗi hàm trực giao 16 số hạng thỏa mãn hoàn hảo hệ điều kiện biên phức tạp ngàm hai cạnh đáy và tự do hai đầu cong.
  • Xác định chính xác bản đồ phân bố ứng suất và độ võng trên lưới 400 nút, chỉ ra các tọa độ nguy hiểm cục bộ tại (10, 13), (10, 20) và (7, 13).
  • Phát triển hoàn thiện phần mềm tính toán tự động trên Visual Basic, hỗ trợ đắc lực cho công tác thiết kế nhanh và tối ưu hóa chiều dày thành vỏ công trình.
  • Trong giai đoạn 12 đến 24 tháng tới, hướng nghiên cứu cần tiếp tục mở rộng khảo sát hiện tượng mất ổn định động học dưới tác động của gió giật và dao động cộng hưởng ngẫu nhiên.

Các kỹ sư kết cấu, nhà thầu xây dựng và các nhà nghiên cứu hãy áp dụng ngay phương pháp tính toán và các thông số hình học tối ưu từ luận văn này để hiện thực hóa những công trình mái vòm composite chất lượng cao, bền vững tuyệt đối trước thiên tai bão gió.