Tổng quan nghiên cứu

Kết cấu tấm và vỏ đóng vai trò cốt lõi trong các công trình kỹ thuật hiện đại, từ công trình dân dụng, cầu đường đến các cấu trúc chịu lực cao trong ngành hàng không vũ trụ và cơ khí chế tạo. Trong thực tế tính toán động lực học và độ ổn định, các mô hình kết cấu phức tạp thường đòi hỏi giải quyết những hệ phương trình tuyến tính quy mô lớn với hơn 1.000.000 bậc tự do. Khối lượng tính toán khổng lồ này gây tiêu tốn nghiêm trọng tài nguyên bộ nhớ máy tính và kéo dài thời gian xử lý dữ liệu. Đồng thời, khi phân tích kết cấu vỏ mỏng với tỉ số chiều dày trên kích thước nhịp từ 0,005 đến 0,1, các mô hình phần tử hữu hạn truyền thống bậc thấp thường xuyên gặp phải hiện tượng khóa cắt (shear locking), dẫn đến việc đánh giá sai lệch độ cứng và tần số dao động tự nhiên.

Nghiên cứu này được thực hiện nhằm xây dựng và phát triển phần tử vỏ phẳng tam giác ba nút cải tiến có tích hợp nút bubble, ký hiệu là bCS-MITC3+, kết hợp với kỹ thuật phân chia kết cấu con theo phương pháp Craig-Bampton truyền thống (CB) và Craig-Bampton nâng cao (ECB). Mục tiêu trọng tâm của đề tài là thiết lập thuật toán phân tích dao động riêng và tải trọng mất ổn định tới hạn của kết cấu tấm và vỏ, đảm bảo độ chính xác vượt trội đồng thời tiết kiệm đáng kể thời gian và bộ nhớ tính toán.

Phạm vi nghiên cứu tập trung vào các kết cấu tấm phẳng đồng nhất đẳng hướng với nhiều điều kiện biên khác nhau, kết cấu vỏ hình trụ và vỏ panel bán cầu chịu nén hoặc cắt trong mặt phẳng. Toàn bộ quy trình thuật toán được lập trình, kiểm chứng độc lập trên nền tảng tính toán số và hoàn thành tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm sai số dao động xuống dưới 0,5% so với lời giải giải tích chính xác, đồng thời cung cấp giải pháp tính toán hiệu năng cao cho các bài toán phân tích động lực học kết cấu phức tạp.

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

Khung lý thuyết áp dụng

Mô hình nghiên cứu được xây dựng trên nền tảng lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin, cho phép xét đến ảnh hưởng của biến dạng cắt trượt qua chiều dày của tấm vỏ. Trường chuyển vị tại một điểm bất kỳ được biểu diễn thông qua ba thành phần chuyển vị mặt trung bình và hai góc xoay pháp tuyến. Mô hình đàn hồi tuân theo định luật Hooke tổng quát cho vật liệu đẳng hướng với mô đun đàn hồi tiêu chuẩn đạt giá trị 2,0 nhân 10 mũ 11 N/m2 và hệ số Poisson bằng 0,3.

Để loại bỏ hoàn toàn hiện tượng khóa cắt, nghiên cứu áp dụng kỹ thuật nội suy hỗn hợp các thành phần ten-xơ (MITC3+) dựa trên việc đánh giá trường biến dạng cắt tại 6 điểm buộc hình học cụ thể với tham số vị trí d bằng 1/10000. Đồng thời, kỹ thuật phần tử hữu hạn trơn dựa trên ô phần tử (CS-FEM) được tích hợp bằng cách bổ sung một nút bubble ở trọng tâm tam giác. Phần tử ban đầu được chia thành 3 miền tam giác con kết nối giữa ba đỉnh và trọng tâm. Nhờ áp dụng định lý phân kỳ Green, tích phân diện tích chứa đạo hàm của hàm dạng trong ma trận độ cứng màng và độ cứng uốn được chuyển đổi hoàn toàn sang tích phân đường dọc theo các biên của miền con.

Bên cạnh đó, khung lý thuyết tổng hợp dạng dao động thành phần (Component Mode Synthesis - CMS) được ứng dụng thông qua hai phương pháp: phương pháp Craig-Bampton tiêu chuẩn (CB) và phương pháp Craig-Bampton nâng cao (ECB). Phương pháp ECB khắc phục triệt để điểm yếu của phương pháp CB truyền thống bằng cách đưa thêm ma trận hiệu chỉnh động lực học chứa các dạng dao động dư, giúp nâng cao độ chính xác của ma trận biến đổi tọa độ tổng thể mà không làm tăng kích thước của mô hình suy giảm.

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

Nghiên cứu sử dụng phương pháp mô phỏng số thông qua việc lập trình toàn bộ thuật toán phần tử hữu hạn trên môi trường tính toán khoa học. Dữ liệu khảo sát bao gồm các mô hình kết cấu tiêu chuẩn quốc tế: tấm vuông chịu uốn và nén, vỏ trụ liên kết ngàm một cạnh và vỏ panel bán cầu liên kết ngàm hai cạnh biên.

Quy trình rời rạc hóa kết cấu được thực hiện qua các cấp độ chia lưới phần tử có kích thước từ 4x4 (gồm 16 phần tử), 8x8 (gồm 64 phần tử), 12x12 (gồm 144 phần tử), 16x16 (gồm 256 phần tử) đến lưới mịn 22x22 (gồm 484 phần tử). Cỡ mẫu lưới này đảm bảo kiểm tra toàn diện quy luật hội tụ của nghiệm số từ thô đến mịn. Phương pháp chọn mẫu chia lưới có tính hệ thống cho phép đánh giá độ nhạy của phần tử trước sự thay đổi của hình học và tỷ lệ chiều dày.

Lý do lựa chọn phương pháp phân tích bCS-MITC3+ kết hợp ECB là vì phương pháp này sử dụng tích phân 2 điểm Gauss trên từng đoạn biên, loại bỏ việc tính đạo hàm trực tiếp của hàm nội suy phức tạp, từ đó triệt tiêu sự phụ thuộc vào chất lượng lưới tam giác bị méo mó. Quá trình đối sánh kết quả được thực hiện trực tiếp với các nghiệm giải tích chính xác của Leissa và nghiệm số từ các phần mềm thương mại hàng đầu như Abaqus, Ansys, Nastran và SAP2000 trong suốt giai đoạn thực hiện đề tài.

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

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

Thứ nhất, trong bài toán phân tích 6 tần số dao động tự nhiên không thứ nguyên đầu tiên của tấm vuông đồng nhất liên kết tựa đơn bốn cạnh (SSSS) có tỉ số chiều dày t/L bằng 0,005, phần tử bCS-MITC3+ kết hợp ECB thể hiện tốc độ hội tụ vượt bậc. Tại lưới chia 4x4, tần số mode 1 của phương pháp bCS-MITC3+_ECB đạt giá trị 4,8962 và nhanh chóng hội tụ về 4,9292 ở lưới 22x22, hoàn toàn trùng khớp với nghiệm giải tích chính xác. Trong khi đó, phần tử DSG3 truyền thống ở lưới 4x4 cho sai số lớn với giá trị lên tới 5,3806.

Thứ hai, đối với tấm vuông liên kết ngàm bốn cạnh (CCCC) có độ mỏng cao với t/L bằng 0,005, hiện tượng khóa cắt ở các phần tử thông thường bị triệt tiêu hoàn toàn. Giá trị mode 1 tính bằng bCS-MITC3+_ECB ở lưới 22x22 đạt 5,4851, tiệm cận tuyệt đối với nghiệm giải tích 5,4851. Ngược lại, phần tử DSG3 ở lưới thô 4x4 bị khóa cứng với tần số lên tới 8,0816, tức là lệch hơn 47% so với thực tế.

Thứ ba, khi phân tích dao động của vỏ panel bán cầu ngàm hai cạnh (CCFF) có bán kính R bằng 1,0 m và chiều dày t bằng 0,04 m, 8 tần số dao động đầu tiên được xác định với độ chính xác cao. Giá trị mode 1 giảm đều từ 342,15 ở lưới 8x8 xuống 331,97 ở lưới 20x20. Mức độ sai lệch giữa phần tử bCS-MITC3+_ECB và kết quả chuẩn hóa từ phần mềm Abaqus ở lưới 16x16 được duy trì ổn định dưới mức 0,4%, chứng minh tính ưu việt của mô hình trên hình học vỏ cong phức tạp.

Thứ tư, trong khảo sát mất ổn định tấm chữ nhật chịu nén một trục, nén hai trục và chịu cắt trong mặt phẳng với các tỉ số kích thước a/b từ 1,0 đến 2,0 và tỉ số chiều dày t/b từ 0,01 đến 0,11, hệ số tải trọng tới hạn hội tụ đơn điệu và ổn định, khắc phục hiện tượng nhảy bậc nghiệm số khi thay đổi điều kiện biên.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp phần tử bCS-MITC3+_ECB đạt hiệu suất cao là nhờ sự kết hợp đồng thời của ba cơ chế cải tiến: hàm bubble bổ sung bậc tự do xoay nội tại, kỹ thuật MITC3+ tái thiết lập trường biến dạng cắt tại 6 điểm buộc hình học, và kỹ thuật làm trơn CS-FEM làm mềm độ cứng giả tạo của phần tử tam giác tiêu chuẩn. Trong phương pháp chia nhỏ kết cấu, phương pháp ECB bổ sung ma trận biến đổi động lực học Tr giúp tính toán chính xác năng lượng dao động dư của các kết cấu con mà không cần mở rộng số lượng bậc tự do biên.

Toàn bộ dữ liệu hội tụ của nghiên cứu có thể được trực quan hóa tối ưu thông qua các biểu đồ đường thể hiện tỷ số tần số chuẩn hóa theo kích thước lưới phần tử N, kết hợp với các bảng ma trận so sánh sai số tương đối giữa CB, ECB, MITC4 và DSG3. Dạng dao động không gian của 6 mode đầu tiên được thể hiện trực quan qua các biểu đồ mặt cong 3D, minh chứng rõ nét khả năng mô tả chính xác chuyển vị uốn và xoắn của kết cấu.

Khi so sánh với các công trình nghiên cứu trước đây, phương pháp đề xuất không chỉ khắc phục nhược điểm khóa cắt của phần tử DSG3 và CS-DSG3 mà còn đạt độ chính xác tương đương với phần tử tứ giác MITC4 và phần mềm SAP2000, trong khi việc chia lưới tam giác linh hoạt hơn rất nhiều đối với các kết cấu có biên hình học phức tạp.

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

Thứ nhất, tích hợp thuật toán bCS-MITC3+_ECB vào các bài toán phân tích phi tuyến hình học và vật liệu: Các viện nghiên cứu và phòng thí nghiệm tính toán kết cấu cần mở rộng mô hình này để khảo sát ứng xử sau mất ổn định và dao động phi tuyến lớn của tấm vỏ trong vòng 12 tháng tới, nhằm nâng cao độ tin cậy khi đánh giá các công trình chịu tải trọng gió bão hoặc động đất cực hạn.

Thứ hai, phát triển thư viện tính toán mã nguồn mở và module thương mại: Nhóm nghiên cứu khuyến nghị các kỹ sư phần mềm phối hợp cùng các nhà khoa học cơ học trong vòng 18 tháng để chuyển đổi các đoạn mã Matlab sang ngôn ngữ C++ hoặc Python hiệu năng cao. Mục tiêu đặt ra là tích hợp phần tử bCS-MITC3+ vào các hệ thống phần mềm tính toán kỹ thuật (CAE), giúp giảm từ 30% đến 40% thời gian chạy mô phỏng cho các doanh nghiệp tư vấn thiết kế.

Thứ ba, mở rộng phạm vi áp dụng cho kết cấu vật liệu tiên tiến: Các đơn vị nghiên cứu chuyên sâu nên triển khai ứng dụng phần tử này cho vật liệu composite phân lớp và vật liệu có cơ tính biến thiên (FGM) trong vòng 24 tháng, hướng tới tối ưu hóa trọng lượng kết cấu trong ngành công nghiệp hàng không vũ trụ và đóng tàu biển với chỉ tiêu tiết kiệm 15% vật liệu chế tạo.

Thứ tư, chuẩn hóa quy trình phân tích kết cấu con (Substructuring) trong thiết kế công nghiệp: Các tổng công ty xây dựng và cơ khí giao thông cần ban hành hướng dẫn kỹ thuật áp dụng phương pháp Craig-Bampton nâng cao cho các kết cấu dạng khung vỏ khẩu độ lớn, đặt mục tiêu tiết kiệm trên 50% dung lượng bộ nhớ máy chủ khi thực hiện các dự án phân tích động lực học quy mô lớn.

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

Thứ nhất, kỹ sư kết cấu công trình dân dụng, công nghiệp và giao thông: Luận văn cung cấp công cụ tính toán chính xác cao để phân tích tần số dao động riêng và kiểm tra điều kiện ổn định của sàn bản phẳng, vách mỏng, dầm bản tổ hợp và vỏ mặt cầu khẩu độ lớn, giúp tối ưu hóa tiết diện và kiểm soát độ an toàn của công trình.

Thứ hai, chuyên gia thiết kế và kỹ sư cơ khí hàng không, vỏ tàu thủy: Tài liệu là tài liệu tham khảo giá trị cho việc mô phỏng ứng xử động lực học của vỏ thân máy bay, tấm vỏ tàu thủy chịu áp lực sóng biển và các tấm panel hình trụ chịu tải trọng nén dọc trục, nơi mà yêu cầu về độ chính xác và giảm trọng lượng là tối quan trọng.

Thứ ba, học viên cao học, nghiên cứu sinh và nhà khoa học chuyên ngành Cơ học tính toán: Đề tài trình bày chi tiết công thức toán học, biến đổi biến dạng màng, uốn trơn và kỹ thuật nén bậc tự do của nút bubble, đóng vai trò như một cẩm nang học thuật để phát triển các loại phần tử hữu hạn mới.

Thứ tư, các nhà phát triển phần mềm phân tích kết cấu (CAE Developer): Luận văn cung cấp thuật toán rõ ràng, công thức giải tích chuẩn xác cùng các ví dụ kiểm chứng số đầy đủ, hỗ trợ việc lập trình mô đun phần tử tam giác phẳng chống khóa cắt trên các nền tảng tính toán thương mại hoặc nội bộ.

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

Phần tử bCS-MITC3+ khắc phục hiện tượng khóa cắt (shear locking) bằng cơ chế nào? Phần tử áp dụng kỹ thuật MITC3+ kết hợp hàm dạng bubble tại trọng tâm tam giác. Bằng cách đánh giá lại biến dạng cắt trượt tại 6 điểm buộc hình học xác định với tham số d bằng 1/10000, trường biến dạng cắt được nội suy độc lập, triệt tiêu năng lượng biến dạng giả tạo khi tỉ số chiều dày t/L giảm về mức 0,005.

Điểm khác biệt cốt lõi giữa phương pháp Craig-Bampton nâng cao (ECB) và Craig-Bampton truyền thống (CB) là gì? Phương pháp CB bỏ qua hoàn toàn ma trận dạng dao động dư của các kết cấu con. Trong khi đó, phương pháp ECB bổ sung ma trận biến đổi động lực học Tr có tính đến ảnh hưởng của dao động dư. Điều này giúp nâng cao độ chính xác của tần số riêng mà vẫn giữ nguyên số lượng bậc tự do của mô hình rút gọn.

Kỹ thuật làm trơn biến dạng ô phần tử (CS-FEM) mang lại lợi ích gì cho việc tính toán ma trận độ cứng? Bằng cách chia phần tử thành 3 tam giác con và sử dụng định lý phân kỳ Green, tích phân diện tích chứa đạo hàm hàm dạng được chuyển thành tích phân đường dọc các biên với 2 điểm Gauss. Kỹ thuật này giúp làm mềm kết cấu một cách hợp lý, tăng tốc độ hội tụ và giảm độ nhạy với lưới tam giác méo.

Độ chính xác của phần tử bCS-MITC3+ khi so sánh với phần mềm thương mại Abaqus và SAP2000 ra sao? Kết quả khảo sát trên vỏ trụ và vỏ bán cầu cho thấy tần số dao động tự nhiên của bCS-MITC3+_ECB sai lệch dưới 0,4% so với Abaqus và tương đương với phần tử tứ giác cao cấp MITC4, vượt trội hơn hẳn các phần tử tam giác tiêu chuẩn như DSG3 và CS-DSG3.

Phần tử này có thể giải quyết những dạng bài toán mất ổn định nào trong thực tế? Mô hình giải quyết xuất sắc các bài toán mất ổn định đàn hồi tới hạn của tấm và vỏ mỏng chịu nén một trục, nén hai trục, nén dọc trục vỏ ống trụ và chịu lực cắt thuần túy trong mặt phẳng với nhiều dạng liên kết biên từ tựa đơn đến ngàm chặt.

Kết luận

  • Xây dựng thành công phần tử vỏ phẳng tam giác ba nút bCS-MITC3+ kết hợp kỹ thuật làm trơn biến dạng và hàm bubble, loại bỏ triệt để hiện tượng khóa cắt trong kết cấu vỏ mỏng.
  • Tích hợp hiệu quả phương pháp Craig-Bampton nâng cao (ECB) vào bài toán giá trị riêng, nâng cao độ chính xác của các tần số dao động riêng và tải trọng ổn định tới hạn.
  • Giảm thiểu đáng kể kích thước hệ phương trình tổng thể và tiết kiệm tài nguyên bộ nhớ máy tính khi xử lý các kết cấu phức tạp có quy mô hàng triệu bậc tự do.
  • Kiểm chứng toàn diện độ chính xác và tính ổn định thuật toán qua hàng loạt ví dụ số tiêu chuẩn cho tấm phẳng, vỏ trụ và vỏ panel bán cầu với sai số dưới 0,5% so với nghiệm giải tích.
  • Cung cấp cơ sở lý thuyết vững chắc và công cụ tính toán số hữu hiệu phục vụ đào tạo sau đại học và nghiên cứu khoa học trong lĩnh vực cơ kỹ thuật công trình.

Đóng góp lớn nhất của luận văn là đã đề xuất một giải pháp tính toán số toàn diện, kết hợp hài hòa giữa kỹ thuật làm trơn phần tử tiên tiến và phương pháp phân rã cấu trúc con hiệu năng cao. Để tiếp tục phát triển hướng nghiên cứu, lộ trình 6 đến 24 tháng tới cần tập trung vào việc mở rộng thuật toán cho kết cấu composite thông minh, phân tích ứng xử phi tuyến động học và xây dựng phần mềm mô phỏng thương mại hoàn chỉnh. Quý độc giả, các nhà nghiên cứu và doanh nghiệp quan tâm có thể khai thác trực tiếp các thuật toán trong luận văn để ứng dụng vào các dự án thiết kế kết cấu công trình thực tế.