Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình dân dụng và công nghiệp, kết cấu khung phẳng chịu lực chiếm khoảng 75% các giải pháp kết cấu chịu tải trọng gió và tải trọng thẳng đứng trong các tòa nhà thấp và trung tầng. Khi quy mô công trình gia tăng, số lượng nút và phần tử thanh dầm, cột tăng nhanh chóng, dẫn đến hệ phương trình vi phân cân bằng cơ học có thể vượt quá 1000 bậc tự do. Bài toán đặt ra là tìm kiếm một phương pháp giải tích số vừa đảm bảo độ chính xác cao, vừa tối ưu hóa tài nguyên tính toán trên máy tính điện tử nhằm xác định chuẩn xác trạng thái ứng suất, nội lực và chuyển vị của hệ kết cấu chịu uốn dưới tác dụng của tải trọng phân bố đều và tải trọng tập trung.

Mục tiêu cốt lõi của nghiên cứu là xây dựng thuật toán và ứng dụng phương pháp phần tử hữu hạn theo mô hình chuyển vị để phân tích nội lực, góc xoay và độ võng của khung phẳng chịu uốn. Phạm vi nghiên cứu giới hạn trong miền đàn hồi tuyến tính của vật liệu với các cấu kiện thanh thỏa mãn giả thiết dầm Euler-Bernoulli, được triển khai và kiểm chứng thực nghiệm số tại thành phố Hải Phòng trong năm 2017.

Ý nghĩa khoa học và thực tiễn của đề tài được lượng hóa rõ nét: việc số hóa bài toán khung phẳng bằng ma trận độ cứng giúp giảm thời gian phân tích kết cấu từ 35% đến 50% so với các phương pháp tính tay truyền thống, đồng thời đảm bảo sai số chuyển vị nút được kiểm soát chặt chẽ ở mức dưới 0,05%, tạo tiền đề vững chắc cho việc tự động hóa thiết kế các công trình dân dụng quy mô lớn.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên sự liên kết chặt chẽ giữa 3 lý thuyết nền tảng trong cơ học kết cấu và cơ học vật rắn biến dạng:

Thứ nhất, lý thuyết dầm chịu uốn Euler-Bernoulli đóng vai trò thiết lập quan hệ vi phân cơ bản giữa mô men uốn M, lực cắt Q và đường độ võng y(x). Lý thuyết này áp dụng cho các cấu kiện thanh dầm có tỷ lệ chiều cao trên chiều dài nhịp h/l nhỏ hơn 1/5 và độ võng giới hạn y_max/h nhỏ hơn 1/5, cho phép bỏ qua biến dạng trượt của các thớ dọc trục để tập trung vào biến dạng uốn thuần túy với phương trình vi phân cấp 4 đồng nhất: EJ nhân đạo hàm bậc 4 của chuyển vị y bằng tải trọng phân bố q.

Thứ hai, nguyên lý biến phân năng lượng bao gồm nguyên lý thế năng biến dạng cực tiểu Castigliano (1847-1884), nguyên lý công bù cực đại và nguyên lý công ảo Gauss. Dưới góc độ cơ học tĩnh, trạng thái cân bằng thực của kết cấu đạt được khi biến phân thế năng toàn phần dừng, tức đạo hàm riêng của phiếm hàm năng lượng theo các chuyển vị nút bằng 0.

Thứ ba, phương trình vi phân chuyển động Lagrange dạng biến phân năng lượng được mở rộng để thiết lập phương trình cân bằng tĩnh học bằng cách sử dụng phép xấp xỉ sai phân hữu hạn cho thế năng biến dạng của 3 điểm nút liên tiếp. Mô hình phần tử hữu hạn sử dụng 4 khái niệm then chốt: hàm dạng đa thức bậc 3 Hermite, ma trận độ cứng phần tử cấp 4x4, ma trận định vị phần tử H và phương pháp đánh số mã liên kết.

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

Nguồn dữ liệu của luận văn được thu thập từ quá trình giải tích số trên tập mẫu gồm 4 phân tố dầm và mở rộng cho hệ kết cấu khung phẳng 2 tầng 2 nhịp với 11 bậc tự do độc lập trong không gian 2 chiều. Phương pháp chọn mẫu chia lưới phần tử được tiến hành có chủ đích: các điểm nút chia bắt buộc phải trùng với vị trí đặt tải trọng tập trung, điểm giao nhau giữa dầm và cột, hoặc nơi có sự thay đổi đột ngột về độ cứng chống uốn EJ.

Lý do lựa chọn phương pháp phần tử hữu hạn theo mô hình chuyển vị là vì phương pháp này cho phép rời rạc hóa các kết cấu phức tạp thành các miền con hình học đơn giản. Thay vì tìm nghiệm giải tích liên tục trên toàn miền vốn rất khó khả thi, phương pháp tìm nghiệm rời rạc tại các nút phần tử thông qua hàm nội suy đa thức bậc 3. Hàm đa thức này đảm bảo điều kiện hội tụ cao, dễ dàng tính toán tích phân ma trận độ cứng bằng phương pháp cầu phương Gauss và hoàn toàn tương thích với các thuật toán đại số ma trận trên máy vi tính.

Toàn bộ quy trình nghiên cứu được thực hiện trong thời gian 12 tháng năm 2017, chia thành 3 giai đoạn: xây dựng cơ sở toán học ma trận, lập trình giải thuật số trên máy tính và kiểm chứng độ chính xác qua các ví dụ thực tế.

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

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

Quá trình phân tích số và kiểm nghiệm lý thuyết đã mang lại 4 phát hiện quan trọng:

Thứ nhất, nghiên cứu đã thiết lập thành công ma trận độ cứng đối xứng cấp 4x4 cho phần tử dầm chịu uốn 2 nút với 4 bậc tự do gồm chuyển vị thẳng đứng v và góc xoay theta tại mỗi đầu nút. Ma trận này chứa các hệ số độ cứng chính xác như 12EJ/L^3, 6EJ/L^2, 4EJ/L và 2EJ/L, phản ánh đầy đủ mối quan hệ tương hỗ giữa lực nút và chuyển vị nút.

Thứ hai, việc áp dụng phương pháp đánh số mã để ghép nối ma trận độ cứng tổng thể đã chứng minh tính ưu việt vượt trội: từ 16 bậc tự do rời rạc ban đầu của 4 phần tử dầm độc lập, thuật toán đã khử liên kết thừa và gộp thành ma trận tổng thể với 11 bậc tự do thực tế, giúp giảm 31,25% dung lượng lưu trữ bộ nhớ và tăng tốc độ xử lý hệ phương trình đại số tuyến tính lên khoảng 42%.

Thứ ba, luận văn đã giải quyết triệt để 2 kỹ thuật xử lý điều kiện biên phức tạp. Đối với chuyển vị bằng 0 tại các gối tựa ngàm hoặc khớp, việc loại bỏ hàng và cột tương ứng giúp ma trận độ cứng tổng thể K không bị suy biến (định thức det K khác 0). Đối với bài toán chuyển vị cưỡng bức có giá trị xác định delta_m = a, kỹ thuật cộng hệ số độ cứng rất lớn A vào phần tử đường chéo k_mm cho kết quả chuyển vị chính xác 100%, sai số so với phương pháp quy đổi tải trọng nút nhỏ hơn 0,01%.

Thứ tư, khi thử nghiệm trên bài toán dầm có độ cứng chống uốn EJ = 10^4 kN.cm^2 chịu tải trọng tập trung P = 10 kN tại giữa nhịp, kết quả độ võng tính toán bằng phần tử hữu hạn đạt giá trị trùng khớp hoàn hảo với nghiệm giải tích sức bền vật liệu, sai số tương đối chỉ ở mức 0,04%.

Thảo luận kết quả

Nguyên nhân chính giúp phương pháp phần tử hữu hạn đạt độ chính xác cao là việc sử dụng hàm dạng nội suy đa thức bậc 3 có 4 tham số độc lập. Hàm đa thức này thỏa mãn hoàn toàn tính liên tục của chuyển vị uốn và tính trơn của góc xoay tại mọi điểm tiếp giáp giữa 2 phần tử liền kề.

Khi so sánh với phương pháp sai phân hữu hạn truyền thống, phương pháp sai phân đòi hỏi bước lưới delta_x phải rất nhỏ và chỉ dùng nội suy tuyến tính giữa các nút, dẫn đến sai số nội lực thường dao động từ 3% đến 5% ở vùng biên. Ngược lại, phương pháp phần tử hữu hạn trong luận văn này chỉ cần chia 4 phần tử đã đạt độ chính xác tương đương với lưới sai phân 20 điểm, giảm hơn 60% số lượng phép tính số học.

Dữ liệu tính toán của luận văn có thể được trình bày rất trực quan: phân bố mô men uốn dọc trục dầm dưới tác dụng của tải trọng phân bố đều q được biểu diễn bằng biểu đồ đường cong Parabol bậc 2 trơn tru, trong khi biểu đồ lực cắt Q thể hiện dưới dạng đường thẳng bậc nhất có bước nhảy đúng bằng giá trị tải trọng tập trung. Các hệ số tương tác độ cứng và véc-tơ tải trọng nút sau khi áp đặt điều kiện biên được tổ chức thành bảng ma trận đối xứng 11x11 rõ ràng, giúp kỹ sư dễ dàng kiểm soát dòng truyền lực trong toàn bộ hệ kết cấu.

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

Nhằm chuyển hóa các kết quả nghiên cứu thành ứng dụng thực tiễn trong công tác thiết kế xây dựng, 4 giải pháp cụ thể được đề xuất như sau:

Thứ nhất, chuẩn hóa quy trình phân chia lưới phần tử dầm và khung trong các tiêu chuẩn thiết kế nội bộ của doanh nghiệp xây dựng. Cần quy định chia tối thiểu từ 4 đến 6 phần tử trên một nhịp dầm chịu tải trọng phân bố đều để khống chế sai số mô men uốn dưới 0,5%. Giải pháp này do các kỹ sư chủ trì kết cấu triển khai áp dụng ngay trong quý 1 sau khi nghiệm thu quy trình.

Thứ hai, phát triển và tích hợp mô đun tính toán ma trận độ cứng theo phương pháp số mã vào các công cụ phần mềm tính kết cấu chuyên dụng tại Việt Nam. Mục tiêu là tự động hóa 100% khâu nhập dữ liệu hình học và điều kiện biên, giảm 30% thời gian thiết lập mô hình tính toán. Thời gian thực hiện kéo dài 6 tháng do nhóm kỹ sư tính toán kết cấu phối hợp cùng chuyên viên công nghệ thông tin thực hiện.

Thứ ba, ứng dụng thuật toán xử lý chuyển vị cưỡng bức bằng số hạng độ cứng lớn để mô phỏng chính xác hiện tượng lún lệch không đều của các móng trụ cầu và móng công trình nhà cao tầng. Giải pháp này giúp kiểm soát độ lệch góc xoay không vượt quá giới hạn cho phép 0,002 radian theo quy chuẩn xây dựng, do các viện nghiên cứu chuyên ngành và công ty tư vấn thiết kế áp dụng định kỳ trong từng dự án cụ thể.

Thứ tư, mở rộng thuật toán phần tử hữu hạn từ bài toán khung phẳng sang hệ khung không gian 3 chiều có xét đến tương tác uốn - xoắn và phi tuyến hình học. Mục tiêu là nâng cao khả năng phân tích ứng xử công trình chịu tải trọng gió bão và động đất lên thêm 25%, lộ trình triển khai trong vòng 12 đến 24 tháng do các nhóm nghiên cứu sau đại học thực hiện.

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

Luận văn là nguồn tư liệu chuyên môn có giá trị thực tiễn cao cho 4 nhóm đối tượng:

Nhóm thứ nhất là các kỹ sư kết cấu công trình dân dụng và công nghiệp. Luận văn cung cấp phương pháp luận chặt chẽ giúp kỹ sư hiểu rõ bản chất cơ học bên trong các phần mềm thương mại phổ biến như SAP2000, ETABS hoặc MIDAS, từ đó dễ dàng phát hiện các lỗi mô hình hóa điều kiện biên và kiểm tra nhanh độ tin cậy của kết quả nội lực.

Nhóm thứ hai là giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Xây dựng và Cơ học Công trình. Đề tài đóng vai trò như một tài liệu tham khảo chuẩn mực về việc ứng dụng nguyên lý biến phân, phương trình Lagrange và thuật toán lập trình phần tử hữu hạn cho các hệ thanh chịu lực.

Nhóm thứ ba là các lập trình viên phần mềm kỹ thuật và kỹ sư tính toán cơ học tính toán (CAE). Cấu trúc thuật toán ma trận 4x4, ma trận định vị H và phương pháp đánh số mã được trình bày chi tiết là nền tảng trực tiếp để lập trình xây dựng các mô đun phần mềm phân tích kết cấu độc lập.

Nhóm thứ tư là sinh viên năm cuối ngành Xây dựng thực hiện đồ án tốt nghiệp công trình bê tông cốt thép hoặc kết cấu thép. Sinh viên có thể tham khảo phương pháp tính chuyển vị nút và cách vẽ biểu đồ nội lực chính xác cho các hệ khung siêu tĩnh phức tạp.

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

Phương pháp phần tử hữu hạn theo mô hình chuyển vị có ưu điểm gì nổi bật so với mô hình ứng suất? Mô hình chuyển vị lấy các chuyển vị và góc xoay tại nút làm ẩn số chính, giúp phương trình cân bằng có dạng tuyến tính đối xứng K nhân delta bằng F. Dạng ma trận này rất thuận tiện cho việc lập trình trên máy tính và tự động thỏa mãn điều kiện liên tục hình học giữa các phần tử, giúp giảm thiểu 40% khối lượng tính toán so với mô hình ứng suất.

Tại sao lại lựa chọn đa thức bậc 3 làm hàm dạng cho phần tử dầm chịu uốn? Đa thức bậc 3 có đúng 4 hệ số chưa biết, hoàn toàn tương thích với 4 bậc tự do tại hai nút của phần tử dầm chịu uốn gồm 2 chuyển vị đứng và 2 góc xoay. Hàm bậc 3 này đảm bảo độ cong đường đàn hồi biến thiên liên tục, giúp kết quả nội lực mô men uốn hội tụ về nghiệm chính xác với sai số dưới 0,05%.

Làm thế nào để xử lý điều kiện biên chuyển vị cưỡng bức tại gối tựa trong bài toán khung phẳng? Luận văn trình bày 2 cách giải quyết: cách thứ nhất là cộng một số hạng rất lớn A vào vị trí đường chéo k_mm của ma trận độ cứng tổng thể và cộng giá trị A nhân a vào véc-tơ tải trọng tương ứng; cách thứ hai là quy đổi chuyển vị cưỡng bức thành véc-tơ lực nút tương đương. Cả hai cách đều cho kết quả chuyển vị chính xác tuyệt đối.

Giả thiết dầm Euler-Bernoulli có giới hạn áp dụng như thế nào trong thực tế tính toán? Giả thiết này giả định mặt cắt ngang luôn phẳng và vuông góc với trục dầm sau biến dạng, bỏ qua biến dạng trượt do lực cắt. Giới hạn áp dụng chuẩn xác là khi tỷ lệ chiều cao tiết diện trên chiều dài dầm h/l nhỏ hơn 1/5 và tỷ lệ độ võng lớn nhất trên chiều cao y_max/h nhỏ hơn 1/5, rất phù hợp với hầu hết dầm khung nhà dân dụng.

Phương pháp đánh số mã liên kết giúp ích gì trong việc tối ưu hóa bộ nhớ máy tính? Phương pháp đánh số mã liên kết cho phép gán trực tiếp từng hệ số độ cứng của phần tử vào đúng tọa độ hàng và cột của ma trận độ cứng tổng thể mà không cần tạo các ma trận chuyển đổi trung gian kích thước lớn. Điều này giúp giảm hơn 30% dung lượng RAM cần thiết khi giải bài toán khung có hàng trăm phần tử.

Kết luận

  • Luận văn đã hệ thống hóa hoàn chỉnh 4 phương pháp xây dựng bài toán cơ học kết cấu và chứng minh tính tương đương giữa phương pháp phương trình vi phân cân bằng phân tố, nguyên lý năng lượng, nguyên lý công ảo và phương trình vi phân Lagrange.
  • Xây dựng thành công ma trận độ cứng đối xứng cấp 4x4 cho phần tử dầm chịu uốn bằng hàm dạng đa thức bậc 3 Hermite, phản ánh chính xác trạng thái biến dạng uốn phẳng.
  • Đề xuất và ứng dụng hiệu quả phương pháp đánh số mã liên kết, giúp rút gọn số bậc tự do từ 16 xuống còn 11 ẩn số thực tế trong mô hình dầm mẫu, tối ưu hóa 31,25% khối lượng dữ liệu.
  • Giải quyết triệt để các dạng điều kiện biên chuyển vị bằng 0 và chuyển vị cưỡng bức với độ chính xác đạt mức sai số cực thấp, dưới 0,04% so với nghiệm giải tích.
  • Cung cấp một quy trình tính toán số hóa mạch lạc, sẵn sàng tích hợp vào các phần mềm tự động hóa thiết kế kết cấu công trình hiện đại.

Đóng góp chính của công trình là cung cấp một giải pháp tính toán phần tử hữu hạn chuẩn xác, tinh gọn cho hệ khung phẳng chịu tải trọng phân bố đều, mở ra hướng nghiên cứu mở rộng sang bài toán khung không gian 3 chiều và tải trọng động lực học trong giai đoạn 12 tháng tới. Các kỹ sư, nhà nghiên cứu và học viên sau đại học hãy tiếp tục ứng dụng và phát triển mô hình số này vào thực tiễn thiết kế để nâng cao chất lượng và độ an toàn cho các công trình xây dựng.