Tổng quan nghiên cứu

Trong kỹ thuật công trình hiện đại, phương pháp phần tử hữu hạn truyền thống luôn đòi hỏi rời rạc hóa toàn bộ thể tích hoặc diện tích hai chiều, tạo ra khối lượng dữ liệu khổng lồ với hàng nghìn bậc tự do đối với các kết cấu phức tạp. Nhằm khắc phục hạn chế về tài nguyên tính toán và nâng cao độ chính xác, phương pháp phần tử biên trung tâm (Scaled Boundary Finite Element Method - SBFEM) đã được phát triển trên thế giới như một giải pháp số đột phá, tích hợp hoàn hảo ưu điểm của phương pháp phần tử hữu hạn và phương pháp phần tử biên.

Tại Việt Nam vào giai đoạn năm 2020, SBFEM vẫn là một hướng đi hoàn toàn mới, gần như chưa có tài liệu tham khảo chuyên sâu về việc phân tích ứng xử của bài toán phẳng khi chịu đồng thời điều kiện biên hỗn hợp và tải trọng bậc cao. Trước thực trạng đó, nghiên cứu của tác giả Nguyễn Thanh Him dưới sự hướng dẫn của Tiến sĩ Nguyễn Văn Chúng tại Khoa Xây dựng, Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh (hoàn thành vào tháng 11 năm 2020) đã được thực hiện nhằm thiết lập phương trình chủ đạo, xây dựng thuật toán số và kiểm chứng độ chính xác của SBFEM.

Nghiên cứu tập trung khảo sát các bài toán tấm phẳng đàn hồi hai chiều có kích thước hình học chuẩn 0,8 m x 0,8 m, chịu tác dụng của tải trọng kéo phân bố đều từ 5 kN/cm2 đến 10 kN/cm2, kết hợp trường trọng lượng bản thân và điều kiện biên chuyển vị phức tạp. Ý nghĩa khoa học và thực tiễn cốt lõi của công trình là chứng minh SBFEM chỉ cần chia lưới trên đường biên một chiều, giúp giảm thiểu hơn 50% số bậc tự do so với các phương pháp số thông thường, đồng thời đạt mức sai số cực tiểu dưới 0,0001% so với nghiệm giải tích chính xác.

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

Khung lý thuyết áp dụng

Khung lý thuyết của đề tài được xây dựng vững chắc dựa trên hai nền tảng cơ học cơ bản: Lý thuyết đàn hồi tuyến tính cho bài toán biến dạng phẳng hai chiều và Nguyên lý công ảo trong cơ học môi trường liên tục. Mô hình nghiên cứu thực hiện ánh xạ hình học từ hệ tọa độ Descartes (x1, x2) sang hệ tọa độ chuyển đổi đặc trưng của SBFEM bao gồm biến bán kính xi (chạy từ tâm định vị đến chu vi miền tính toán) và biến chu vi s (dọc theo đường biên ngoài).

Các khái niệm chuyên ngành then chốt được phát triển xuyên suốt mô hình gồm:

  1. Tâm định vị SC (Scaling Center): Điểm mốc hình học được lựa chọn linh hoạt ở bên trong hoặc ngay trên đường biên của kết cấu để phát các tia tỷ lệ bao phủ toàn miền.
  2. Ma trận chuyển đổi Jacobian và ma trận quan hệ biến dạng - chuyển vị: Công cụ toán học cho phép xấp xỉ biên một chiều thông qua hàm dạng đa thức bậc 1.
  3. Phương trình vi phân cấp hai không thuần nhất: Phương trình chủ đạo mô tả trường chuyển vị chưa biết, tích hợp đầy đủ tác động của trường lực khối và lực mặt bên.
  4. Lời giải bán giải tích: Nghiệm của bài toán được giải tích hóa hoàn toàn theo phương bán kính thông qua việc giải bài toán trị riêng và véc-tơ riêng của ma trận hệ thống, trong khi chỉ xấp xỉ số theo phương chu vi.

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

Nguồn dữ liệu và mô hình kiểm chứng được thiết lập từ các bài toán mẫu kinh điển trong cơ học kết cấu với thông số vật liệu đồng nhất đẳng hướng: mô đun đàn hồi E = 2,0 x 10^6 kN/cm2, hệ số Poisson v = 0,3 và kích thước cạnh kết cấu 0,8 m x 0,8 m.

Phương pháp phân tích được thực hiện bằng cách xây dựng quy trình số hóa và lập trình toàn diện trên môi trường ngôn ngữ Matlab. Thuật toán tiến hành rời rạc hóa đường biên, thiết lập các ma trận độ cứng thành phần, giải hệ phương trình vi phân cấp hai để tìm nghiệm chuyển vị tại các nút biên, sau đó suy biến nghiệm ra toàn bộ miền liên tục bên trong.

Cỡ mẫu và phương pháp chọn mẫu rời rạc được phân chia thành các cấp độ lưới tăng dần có chủ đích: khảo sát với 4 phần tử (8 bậc tự do), 8 phần tử (16 bậc tự do), 16 phần tử (32 bậc tự do) khi đặt tâm SC ở trung tâm tấm; và 2 phần tử, 4 phần tử, 8 phần tử khi đặt tâm SC tại góc biên. Phương pháp chọn mẫu chia lưới biên đồng nhất được áp dụng để đánh giá trực tiếp quy luật hội tụ số.

Lý do lựa chọn phương pháp phân tích SBFEM là nhờ tính chất giải tích theo phương bán kính, loại bỏ hoàn toàn nhu cầu tạo lưới phần tử bên trong miền tính toán, từ đó triệt tiêu sai số hình học nội suy vốn là nhược điểm cố hữu của phần tử hữu hạn. Toàn bộ timeline nghiên cứu và tính toán số được hoàn thiện tập trung trong năm 2020.

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

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

Thứ nhất, SBFEM đạt độ chính xác phi thường ngay từ mức độ chia lưới thô nhất. Khi phân tích bài toán tấm phẳng chịu tải trọng kéo phân bố đều q = 5 kN/cm2, chỉ với 4 phần tử biên (8 bậc tự do), ứng suất tính toán đạt 4,9999990 kN/cm2, chênh lệch so với nghiệm giải tích 5,0000 kN/cm2 ở mức sai số chỉ 2,09 x 10^-5 %. Khi tăng số lượng phần tử lên 16 (32 bậc tự do), sai số tiếp tục giảm mạnh xuống mức 1,17 x 10^-5 %.

Thứ hai, vị trí đặt tâm định vị SC mang lại tính linh hoạt cao mà không làm suy giảm độ chính xác. Khi khảo sát tấm phẳng chịu tải trọng q = 10 kN/cm2 với điều kiện biên hỗn hợp, mô hình đặt tâm SC ở giữa tấm (4 phần tử biên) cho sai số chuyển vị u1 là 1,09 x 10^-4 % và sai số ứng suất là 1,41 x 10^-4 %. Đáng chú ý, khi chuyển tâm SC ra góc biên với chỉ 2 phần tử biên trên cạnh tự do, kết quả thu được hoàn toàn tương đương với cùng mức sai số 1,09 x 10^-4 %, giúp tiết kiệm 50% số lượng phần tử cần tính toán.

Thứ ba, tốc độ hội tụ của nghiệm chuyển vị diễn ra vô cùng nhanh chóng. Thành phần chuyển vị u2 tại điểm kiểm tra M có sai số giảm từ 5,17 x 10^-4 % ở lưới 4 phần tử xuống còn 2,48 x 10^-4 % ở lưới 16 phần tử, tương ứng mức cải thiện độ chính xác hơn 52% khi tăng mật độ rời rạc biên.

Thứ tư, trường chuyển vị và ứng suất phân bố hoàn toàn trơn tru dọc theo các cạnh biên. Các giá trị chuyển vị và ứng suất chuẩn hóa theo công thức không thứ nguyên trên cạnh AC thể hiện sự trùng khít tuyệt đối với đường cong giải tích lý thuyết trên toàn bộ dải tọa độ từ 0,0 đến 1,0.

Thảo luận kết quả

Độ chính xác vượt trội của phương pháp bắt nguồn từ bản chất bán giải tích của hệ phương trình SBFEM. Do trường chuyển vị theo phương bán kính xi được biểu diễn chính xác bằng hàm mũ thông qua các trị riêng của ma trận hệ thống, nên không xảy ra hiện tượng méo mó lưới hay sai số tích phân thể tích.

So với các nghiên cứu quốc tế của Song và Wolf hay nghiên cứu của Rungamornrat và Chung, kết quả của luận văn đã khẳng định tính đúng đắn khi mở rộng SBFEM cho bài toán có điều kiện biên hỗn hợp và tải trọng đa thức bậc cao. Phương pháp thể hiện ưu thế vượt bậc so với phương pháp phần tử hữu hạn truyền thống đang sử dụng tại Việt Nam về cả tốc độ hội tụ lẫn dung lượng bộ nhớ.

Toàn bộ dữ liệu số trong luận văn được tổng hợp khoa học qua 6 bảng biểu chi tiết, so sánh trực tiếp từng giá trị ứng suất và chuyển vị theo các mức bậc tự do khác nhau. Các đồ thị biểu diễn chuyển vị chuẩn hóa u1 nhân với mô đun đàn hồi trên tải trọng và kích thước hình học thể hiện rõ tính liên tục tuyệt đối, không xuất hiện hiện tượng dao động số hay tập trung ứng suất giả tạo tại các vị trí nút biên.

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

Thứ nhất, chuẩn hóa quy trình ứng dụng SBFEM vào tính toán kỹ thuật: Các viện thiết kế và kỹ sư kết cấu công trình cần đưa phương pháp phần tử biên trung tâm vào quy trình phân tích móng mố công trình và bài toán cơ học đất, đặt mục tiêu giảm thiểu 40% thời gian xử lý tiền kỳ chia lưới trong giai đoạn 2026 - 2028.

Thứ hai, xây dựng các gói thư viện mã nguồn mở: Các nhóm nghiên cứu phát triển phần mềm tính toán kỹ thuật tại Việt Nam nên phối hợp với các trường đại học để đóng gói mã nguồn Matlab của luận văn thành các module tính toán chuyên dụng, đảm bảo kiểm soát sai số tính toán dưới mức 0,001% trước quý 4 năm 2027.

Thứ ba, mở rộng phạm vi nghiên cứu sang bài toán không gian ba chiều: Các nghiên cứu sinh và học viên cao học thuộc chuyên ngành xây dựng cần tiếp tục phát triển hệ phương trình SBFEM cho bài toán kết cấu không gian 3D và vật liệu dị hướng nhiều lớp, khảo sát dải quy mô từ 32 đến 128 phần tử biên trong vòng 2 năm tới.

Thứ tư, khai thác tối đa kỹ thuật đặt tâm định vị trên biên: Trong thực tế xử lý các kết cấu có biên đối xứng hoặc góc cạnh phức tạp, người thực hiện nên ưu tiên đặt tâm định vị SC tại các điểm góc để cắt giảm 50% số lượng biên phải chia phần tử, tối ưu hóa năng lực tính toán cho các dự án quy mô lớn.

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

Thứ nhất, kỹ sư thiết kế kết cấu công trình: Những người cần một phương pháp số có độ chính xác cao để phân tích các cấu kiện chịu tải trọng kéo nén phức tạp, kích thước từ 0,8 m đến hàng chục mét, giúp giảm thiểu đáng kể thời gian thiết lập mô hình hình học.

Thứ hai, giảng viên và học viên cao học chuyên ngành Cơ kỹ thuật và Xây dựng: Tài liệu cung cấp toàn bộ cơ sở toán học giải tích chặt chẽ, các bước thiết lập phương trình dạng yếu và lưu đồ thuật toán Matlab chi tiết, là nguồn học liệu mẫu mực cho các đề tài nghiên cứu chuyên sâu về phương pháp số.

Thứ ba, các nhà phát triển phần mềm kỹ thuật hỗ trợ máy tính (CAE): Các chuyên gia lập trình có thể tham khảo thuật toán thiết lập ma trận độ cứng biên K và kỹ thuật giải bài toán trị riêng để tối ưu hóa hiệu suất CPU và bộ nhớ RAM gấp 2 đến 3 lần so với lõi tính toán FEM truyền thống.

Thứ tư, các nhà nghiên cứu trong lĩnh vực cơ học nứt và địa kỹ thuật: SBFEM là công cụ lý tưởng để mô phỏng sự tập trung ứng suất và lan truyền vết nứt với độ tin cậy kết quả đạt trên 99,99% nhờ khả năng mô tả chính xác trường kỳ dị tại tâm định vị mà không cần làm mịn lưới.

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

  1. Phương pháp phần tử biên trung tâm SBFEM khác biệt cơ bản như thế nào so với FEM truyền thống? SBFEM chỉ thực hiện rời rạc hóa trên đường biên một chiều của bài toán thay vì phải chia lưới toàn bộ diện tích hai chiều như FEM. Phương bán kính từ tâm định vị đến biên được giải hoàn toàn bằng nghiệm giải tích chính xác, giúp giảm hơn 50% số bậc tự do và đạt độ chính xác sai số chỉ 2,09 x 10^-5 %.

  2. Việc thay đổi vị trí tâm định vị SC có làm sai lệch kết quả phân tích kết cấu không? Không. Kết quả thực nghiệm trong luận văn chứng minh việc đặt tâm SC ở giữa tấm với 4 phần tử hay đặt tại góc biên với 2 phần tử đều cho sai số ứng suất tương đương nhau ở mức 1,41 x 10^-4 %, đồng thời việc đặt tâm ở góc còn giúp giảm một nửa số lượng phần tử biên cần mô hình hóa.

  3. SBFEM có giải quyết được bài toán chịu tải trọng trọng lượng bản thân và tải trọng bậc cao không? Phương pháp hoàn toàn giải quyết chính xác. Bằng cách thiết lập nghiệm riêng đại số tích hợp trực tiếp vào phương trình dạng yếu, SBFEM xử lý trọn vẹn trường trọng lượng riêng và tải trọng phân bố đa thức bậc cao mà không làm suy giảm tốc độ hội tụ của hệ phương trình.

  4. Tại sao SBFEM đạt được độ hội tụ rất nhanh chỉ với số lượng phần tử biên rất ít? Bởi vì nghiệm theo phương bán kính là nghiệm giải tích thực thụ của phương trình vi phân đàn hồi. Khi khảo sát tấm phẳng 0,8 m x 0,8 m, chỉ cần 4 phần tử biên với 8 bậc tự do là kết quả ứng suất đã tiệm cận nghiệm chính xác với độ lệch chỉ 0,000001 kN/cm2.

  5. Ngôn ngữ lập trình nào được sử dụng để hiện thực hóa thuật toán trong luận văn? Toàn bộ quy trình tính toán SBFEM được lập trình tự động hóa trên môi trường ngôn ngữ Matlab. Chương trình xử lý trơn tru các phép tính ma trận độ cứng, giải bài toán trị riêng và trích xuất đồ thị ứng suất chỉ trong vài giây cho các mô hình từ 8 đến 32 bậc tự do.

Kết luận

  • Luận văn đã thiết lập thành công hệ phương trình chủ đạo của phương pháp phần tử biên trung tâm SBFEM cho bài toán phẳng hai chiều chịu điều kiện biên hỗn hợp và tải trọng bậc cao.
  • Kiểm chứng độ chính xác vượt trội của phương pháp với sai số ứng suất chỉ 2,09 x 10^-5 % khi sử dụng mô hình lưới thô 4 phần tử (8 bậc tự do).
  • Chứng minh tính linh hoạt của vị trí tâm định vị SC khi đặt tại góc biên, giúp giảm 50% số phần tử biên mà vẫn bảo toàn mức sai số cực nhỏ 1,09 x 10^-4 %.
  • Xây dựng hoàn chỉnh lưu đồ tính toán và bộ chương trình số hóa tự động trên nền tảng ngôn ngữ Matlab với tốc độ xử lý nhanh chóng và độ ổn định cao.
  • Đặt nền móng khoa học quan trọng cho việc phổ biến và ứng dụng các phương pháp số tiên tiến trong ngành kỹ thuật xây dựng tại Việt Nam từ năm 2020 trở đi.

Đóng góp lớn nhất của công trình là mở ra hướng tiếp cận mới trong tính toán cơ học kết cấu, giúp tối ưu hóa thời gian và tài nguyên phần cứng. Quý độc giả, các nhà nghiên cứu và kỹ sư công trình hãy tham khảo chi tiết toàn văn luận văn thạc sĩ của tác giả Nguyễn Thanh Him để ứng dụng ngay giải pháp SBFEM vào các bài toán phân tích kết cấu thực tế.