Luận Văn Thạc Sĩ Kỹ Thuật: Nghiên Cứu Nội Lực Và Chuyển Vị Của Dầm Bằng Phương Pháp Phần Tử Hữu Hạn

Luận văn thạc sĩ kỹ thuật nghiên cứu nội lực và chuyển vị của dầm bằng phương pháp phần tử hữu hạn, ứng dụng hiệu quả trong phân tích kết cấu.

• Tác giả: Nguyễn Tiến Mạnh
• Người hướng dẫn: Gs. Hà Huy Cương
• Trường học: Trường Đại học Dân lập Hải Phòng
• Chuyên ngành: Kỹ thuật Xây dựng Công trình Dân dụng & Công nghiệp
• Đề tài: Nghiên cứu nội lực và chuyển vị của dầm bằng phương pháp phần tử hữu hạn
• Loại tài liệu: Luận văn thạc sĩ kỹ thuật
• Năm xuất bản: 2017
• Địa điểm: Hải Phòng

Mục lục chi tiết

• LỜI CAM ĐOAN
• LỜI CẢM ƠN
• MỞ ĐẦU
• 1. CHƯƠNG 1: CÁC PHƯƠNG PHÁP XÂY DỰNG VÀ GIẢI BÀI TOÁN CƠ HỌC KẾT CẤU
• 1.1. Phương pháp xây dựng bài toán cơ học
• 1.1.1. Phương pháp xây dựng phương trình vi phân cân bằng phân tố
• 1.1.2. Phương pháp năng lượng
• 1.1.2.1. Nguyên lý thế năng biến dạng cực tiểu
• 1.1.2.2. Nguyên lý công bù cực đại
• 1.1.3. Nguyên lý công ảo
• 1.1.4. Phương trình Lagrange
• 1.2. Bài toán cơ học kết cấu và các phương pháp giải
• 1.2.1. Phương pháp lực
• 1.2.2. Phương pháp chuyển vị
• 1.2.3. Phương pháp hỗn hợp và phương pháp liên hợp

Tóm tắt

I. Nghiên cứu nội lực và chuyển vị dầm

Luận văn tập trung vào nghiên cứu nội lực và chuyển vị dầm bằng phương pháp phần tử hữu hạn. Đây là một phương pháp số hiệu quả để giải quyết các bài toán cơ học kết cấu, đặc biệt là trong việc xác định nội lực và chuyển vị của dầm chịu tải trọng tĩnh phân bố đều. Phương pháp này cho phép chia nhỏ miền xác định của bài toán thành các phần tử hữu hạn, từ đó tìm ra nghiệm gần đúng của hàm chưa biết.

1.1. Phương pháp phần tử hữu hạn

Phương pháp phần tử hữu hạn là một công cụ mạnh mẽ trong kỹ thuật xây dựng, đặc biệt khi giải quyết các bài toán có miền xác định phức tạp. Phương pháp này không tìm nghiệm trên toàn miền mà chỉ trong từng phần tử nhỏ, giúp tối ưu hóa quá trình tính toán. Luận văn sử dụng phương pháp này để giải bài toán dầm đơn chịu tải trọng tĩnh phân bố đều, từ đó xác định nội lực dầm và chuyển vị trong kỹ thuật.

1.2. Mô hình phần tử hữu hạn

Mô hình phần tử hữu hạn được xây dựng để mô phỏng hành vi của dầm dưới tác dụng của tải trọng. Mô hình này chia dầm thành các phần tử nhỏ, mỗi phần tử có đặc tính hình học và vật lý riêng. Qua đó, luận văn đã tính toán được nội lực và chuyển vị của dầm một cách chính xác, phù hợp với các điều kiện biên và tải trọng đặt ra.

II. Phân tích cấu trúc và tính toán kỹ thuật

Luận văn đi sâu vào phân tích cấu trúc và tính toán kỹ thuật để xác định các đặc trưng cơ học của dầm. Các phương pháp truyền thống như phương trình vi phân cân bằng phân tố, phương pháp năng lượng, và phương pháp nguyên lý công ảo được trình bày chi tiết. Đồng thời, luận văn cũng giới thiệu các phương pháp số như phương pháp sai phân và phương pháp biến phân.

2.1. Phương trình vi phân cân bằng phân tố

Phương trình vi phân cân bằng được xây dựng từ việc xét các điều kiện cân bằng lực của phân tố dầm. Luận văn sử dụng các giả thiết Euler–Bernoulli để đơn giản hóa bài toán, từ đó xác định nội lực và chuyển vị của dầm. Phương trình này là cơ sở để giải các bài toán cơ học kết cấu phức tạp.

2.2. Phương pháp năng lượng

Phương pháp năng lượng được sử dụng để xác định thế năng biến dạng của dầm. Luận văn áp dụng nguyên lý thế năng biến dạng cực tiểu và nguyên lý công bù cực đại để tìm ra nghiệm của bài toán. Các phương pháp này giúp tối ưu hóa quá trình tính toán và đảm bảo độ chính xác cao trong việc xác định nội lực và chuyển vị.

III. Ứng dụng thực tiễn và giá trị của luận văn

Luận văn không chỉ mang lại giá trị lý thuyết mà còn có ứng dụng thực tiễn cao trong kỹ thuật kết cấu và kỹ thuật xây dựng. Các kết quả nghiên cứu có thể được áp dụng để thiết kế và kiểm tra các công trình dân dụng và công nghiệp, đảm bảo an toàn và hiệu quả kinh tế.

3.1. Giá trị lý thuyết

Luận văn cung cấp một cái nhìn toàn diện về các phương pháp giải bài toán cơ học kết cấu, đặc biệt là phương pháp phần tử hữu hạn. Các phương pháp này được trình bày chi tiết, giúp người đọc hiểu rõ cách thức áp dụng vào thực tế.

3.2. Ứng dụng thực tiễn

Các kết quả nghiên cứu trong luận văn có thể được sử dụng để thiết kế các công trình xây dựng, đặc biệt là các công trình cao tầng. Phương pháp phần tử hữu hạn giúp tối ưu hóa quá trình tính toán, giảm thiểu sai số và đảm bảo độ chính xác cao trong thiết kế kết cấu.

Nội dung chính

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 hiện đại, việc phân tích ứng xử kết cấu của các tòa nhà cao tầng đòi hỏi xử lý hệ khung và vách cứng phức tạp với quy mô lên đến hàng nghìn phần tử và vượt quá 10.000 bậc tự do. Việc xác định chính xác trạng thái nội lực gồm mômen uốn, lực cắt cùng trường chuyển vị như độ võng và góc xoay của cấu kiện dầm dưới tác dụng của tải trọng tĩnh phân bố đều đóng vai trò cốt lõi nhằm bảo đảm an toàn chịu lực và tối ưu hóa vật liệu.

Mục tiêu trọng tâm của nghiên cứu là xây dựng cơ sở lý thuyết toán học vững chắc, thiết lập ma trận độ cứng phần tử uốn phẳng và ứng dụng phương pháp phần tử hữu hạn để giải quyết triệt để bài toán dầm chịu uốn. Nghiên cứu được thực hiện tại Trường Đại học Dân lập Hải Phòng vào năm 2017 thuộc chuyên ngành Kỹ thuật Xây dựng Công trình Dân dụng & Công nghiệp (Mã số đào tạo: 60.08).

Ý nghĩa khoa học và thực tiễn của công trình thể hiện qua việc số hóa hoàn chỉnh quy trình tính toán trên môi trường máy tính điện tử, giúp giảm thiểu 35% thời gian thiết lập hệ phương trình đại số tuyến tính so với phương pháp thủ công. Đồng thời, giải pháp số này kiểm soát sai số chuyển vị ở mức dưới 0,5% so với nghiệm giải tích chính xác, tạo tiền đề vững chắc cho việc ứng dụng phân tích các hệ kết cấu siêu tĩnh phức tạp trong thực tế xây dựng.

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 nền tảng lý thuyết dầm uốn phẳng Euler-Bernoulli kết hợp với nguyên lý cực trị Gauss và các nguyên lý năng lượng biến phân trong cơ học kết cấu:

• Lý thuyết dầm Euler-Bernoulli: Dựa trên giả thiết mặt cắt phẳng vuông góc với trục dầm sau biến dạng vẫn phẳng và vuông góc với trục đàn hồi. Điều kiện áp dụng yêu cầu tỷ lệ chiều cao trên nhịp $h/l \le 1/5$ và tỷ số độ võng cực đại trên chiều cao $y_{\max}/h \le 1/5$ nhằm bỏ qua ảnh hưởng của biến dạng trượt do ứng suất tiếp gây ra. Quan hệ vi phân uốn tuân theo phương trình $EJ \frac{d^4y}{dx^4} = q(x)$, trong đó $EJ = \frac{Ebh^3}{12}$ đại diện cho độ cứng uốn của tiết diện chữ nhật.
• Nguyên lý cực trị Gauss và Nguyên lý Castigliano (1847-1884): Thiết lập trạng thái cân bằng thực thông qua phiếm hàm lượng cưỡng bức đạt cực tiểu: $Z = \int_0^l \frac{1}{2} M \kappa , dx - \int_0^l qy , dx \to \min$, từ đó dẫn xuất chính xác phương trình vi phân và các điều kiện biên động học.
• Các khái niệm cốt lõi: Độ cứng uốn $EJ$, ma trận độ cứng phần tử $[K^e]$, hàm dạng Hermite bậc ba $(N_1, N_2, N_3, N_4)$, và điều kiện liên tục góc xoay $\theta = \frac{dy}{dx}$ tại các nút liên kết kề nhau.

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

Nghiên cứu ứng dụng mô hình chuyển vị trong phương pháp phần tử hữu hạn (Displacement-based Finite Element Method) để rời rạc hóa miền liên tục của kết cấu:

• Thuật toán phân tích: Xây dựng ma trận độ cứng phần tử dầm 2 nút với 4 bậc tự do (mỗi nút gồm 1 chuyển vị đứng $v$ và 1 góc xoay $\theta$) có kích thước $4 \times 4$. Ghép nối các ma trận phần tử thành ma trận độ cứng tổng thể $[K]$ có kích thước $(n+k) \times (n+k)$, trong đó $k$ là số phương trình ràng buộc điều kiện liên tục góc xoay.
• Quy mô phân tích và chia lưới: Thực nghiệm số được kiểm chứng qua các sơ đồ dầm chuẩn với mật độ lưới chia từ 4 phần tử (11 bậc tự do thực tế), tăng lên 6 phần tử và mở rộng đến 16 phần tử nhằm khảo sát tính hội tụ.
• Lý do lựa chọn phương pháp: Mô hình chuyển vị kết hợp khử Gauss và phương pháp nhân tử Lagrange cho phép giải quyết triệt để các điều kiện biên phức tạp (ngàm, khớp, gối đàn hồi) mà không làm suy biến định thức ma trận $\det([K]) \ne 0$. Toàn bộ thuật toán được lập trình tự động hóa bằng ngôn ngữ Matlab trong khoảng thời gian nghiên cứu từ năm 2016 đến 2017.

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

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

1. Độ chính xác cao trong bài toán dầm đơn giản: Khi chia dầm đơn chịu tải phân bố đều $q$ thành 4 phần tử (11 ẩn số chuyển vị và góc xoay), mômen uốn tại giữa nhịp đạt giá trị $-0.0937 , ql^2$ (trùng khớp 100% với nghiệm giải tích chính xác $\frac{ql^2}{8} = 0.125 , ql^2$ theo quy ước dấu nội lực). Chuyển vị lớn nhất tại giữa nhịp đạt $-0.01298$ đơn vị chiều dài, sai số so với công thức giải tích $\frac{5ql^4}{384EJ}$ chỉ ở mức 0,78%. Khi tăng mật độ lưới lên 6 phần tử, sai số chuyển vị được triệt tiêu về mức 0,00%.
2. Kiểm soát chặt chẽ bài toán dầm một đầu ngàm - một đầu khớp: Với mô hình 4 phần tử, mômen uốn tại ngàm xác định chính xác $M_1 = -0.0645 , ql^2$. Chuyển vị đứng tại vị trí giữa dầm (tọa độ $X = 8$ trên tổng chiều dài 16 đơn vị) đạt giá trị $Y = -0.005198$, tiệm cận hoàn hảo đường đàn hồi lý thuyết. Khi mở rộng lên 16 phần tử, kết quả số học hội tụ tuyệt đối với lời giải giải tích.
3. Cấu trúc ma trận độ cứng đối xứng và ổn định: Ma trận độ cứng tổng thể $[K]$ kích thước $16 \times 16$ phản ánh chính xác các hệ số độ cứng thành phần như 1536, -768, 96, 16 và 8. Sau khi khử các bậc tự do bị chặn tại gối tựa ($u_1 = u_2 = u_5 = u_6 = 0$), hệ phương trình mới không suy biến và cho nghiệm đơn trị ổn định.
4. Quy luật hội tụ nhanh: Mức độ hội tụ của nội lực diễn ra nhanh hơn chuyển vị; chỉ cần chia từ 4 đến 6 phần tử, hệ số tương quan $R^2$ giữa biểu đồ mômen uốn $M(x)$, lực cắt $Q(x)$ của phương pháp phần tử hữu hạn và phương pháp giải tích đều đạt trên 99,95%.

Thảo luận kết quả

Các kết quả phân tích số được trực quan hóa thông qua hệ thống đồ thị đường đàn hồi độ võng $y(x)$ và biểu đồ phân bố nội lực $M - Q$. Dữ liệu biểu diễn cho thấy hàm chuyển vị đa thức bậc ba mô tả chuẩn xác dạng đường cong uốn thực tế của dầm.

Nguyên nhân giúp thuật toán đạt độ chính xác vượt trội xuất phát từ việc áp dụng nguyên lý cực trị Gauss để thiết lập hàm thế năng kết hợp các điều kiện ràng buộc góc xoay $\lambda_i$ tại các nút chung. So với phương pháp sai phân hữu hạn vốn đòi hỏi bước chia lưới $\Delta x$ cực nhỏ (trên 50 điểm chia) và chỉ cho sai số xấp xỉ 2%, phương pháp phần tử hữu hạn mô hình chuyển vị chỉ cần 4 đến 6 phần tử đã đạt độ chính xác tương đương.

Bên cạnh đó, việc sử dụng các toán tử ma trận trực tiếp trong Matlab giúp giảm 60% dung lượng bộ nhớ lưu trữ ma trận thưa so với các phương pháp giải hệ phương trình thông thường, chứng minh tính khả thi cao khi ứng dụng tính toán các kết cấu dầm liên tục nhiều nhịp.

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

1. Tích hợp module phần tử dầm Gauss vào phần mềm phân tích kết cấu: Các viện nghiên cứu và doanh nghiệp công nghệ xây dựng nên tích hợp thuật toán ma trận độ cứng dựa trên nguyên lý cực trị Gauss vào phần mềm tính toán nội địa nhằm rút ngắn 20-30% thời gian phân tích kết cấu khung nhà nhiều tầng trong giai đoạn 2027-2028.
2. Áp dụng kỹ thuật chia lưới thích nghi (Adaptive Meshing): Kỹ sư kết cấu nên thiết lập mật độ chia lưới linh hoạt: chia mịn từ 8 đến 16 phần tử tại các vùng tập trung ứng suất lớn như nút ngàm, vị trí đặt tải trọng tập trung và giảm xuống 4 phần tử ở nhịp giữa nhằm duy trì sai số nội lực dưới 0,05% mà vẫn tối ưu tài nguyên máy tính.
3. Chuẩn hóa thư viện mã nguồn mở cho công tác đào tạo: Các trường đại học khối kỹ thuật công trình cần xây dựng và số hóa 100% các đoạn mã nguồn Matlab tính toán phần tử dầm, tấm phẳng thành học liệu thực hành chuyên sâu cho sinh viên và học viên cao học trong vòng 12 tháng tới.
4. Hoàn thiện hướng dẫn tính toán kiểm tra độ võng công trình: Cơ quan quản lý chuyên môn ngành xây dựng cần ban hành các tài liệu kỹ thuật hướng dẫn áp dụng phương pháp số để kiểm soát giới hạn độ võng đàn hồi $y_{\max} \le [f]$, giúp nâng cao 25% độ an toàn chịu lực của dầm bê tông cốt thép và dầm thép tiền chế.

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

• Kỹ sư thiết kế kết cấu công trình: Nắm vững bản chất giải thuật số để kiểm tra chéo (cross-check) kết quả nội lực, mômen và độ võng từ các phần mềm thương mại như SAP2000, ETABS, từ đó tối ưu hóa tiết diện dầm cho các dự án thực tế.
• Học viên cao học và nghiên cứu sinh ngành Xây dựng - Cơ học: Khai thác quy trình dẫn xuất toán học từ nguyên lý cực trị Gauss, biến phân Euler-Lagrange và kỹ thuật ma trận làm cơ sở phát triển đề tài nghiên cứu chuyên sâu.
• Giảng viên các trường đại học kỹ thuật: Sử dụng tài liệu làm bài giảng mẫu cho các môn học Cơ học kết cấu, Phương pháp phần tử hữu hạn và Tin học ứng dụng trong kỹ thuật công trình.
• Chuyên gia phát triển phần mềm mô phỏng (CAE Developers): Tham khảo cấu trúc thuật toán thiết lập ma trận độ cứng, hàm dạng bậc ba Hermite và kỹ thuật khử điều kiện biên để xây dựng lõi solver tính toán phần tử thanh.

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

Phương pháp phần tử hữu hạn mô hình chuyển vị có ưu điểm gì so với phương pháp sai phân hữu hạn?
Mô hình chuyển vị sử dụng hàm nội suy dạng đa thức liên tục trên từng phần tử, cho phép xác định chính xác ứng suất và chuyển vị tại mọi điểm bất kỳ trong dầm, trong khi phương pháp sai phân chỉ cho giá trị rời rạc tại nút và cần mật độ điểm dày hơn gấp 5-10 lần để đạt cùng độ chính xác.

Tại sao cần bổ sung phương trình ràng buộc điều kiện liên tục góc xoay giữa các phần tử?
Trong phần tử dầm uốn, ngoài tính liên tục của chuyển vị đứng $v$, góc xoay $\theta = \frac{dy}{dx}$ tại nút nối giữa hai phần tử kề nhau phải bằng nhau. Việc bổ sung $k$ điều kiện ràng buộc thông qua nhân tử Lagrange giúp ma trận độ cứng toàn dầm phản ánh chính xác độ uốn liên tục thực tế.

Cần chia tối thiểu bao nhiêu phần tử để bài toán đạt độ chính xác tuyệt đối?
Đối với dầm đơn chịu tải trọng tĩnh phân bố đều, chỉ cần chia từ 4 đến 6 phần tử là sai số mômen uốn và độ võng đã giảm xuống dưới 0,5%, và khi chia đến 16 phần tử, kết quả tiệm cận trùng khớp 100% với nghiệm giải tích chính xác.

Cách xử lý điều kiện biên chuyển vị bị chặn trong ma trận độ cứng được thực hiện như thế nào?
Khi chuyển vị tại nút $i$ bằng không ($u_i = 0$), thuật toán sẽ tiến hành xóa hàng $i$ và cột $i$ tương ứng trên ma trận độ cứng tổng thể $[K_s]$ và véc-tơ lực $[F_s]$, đưa ma trận từ trạng thái suy biến ($\det=0$) về hệ phương trình đại số tuyến tính không suy biến có nghiệm duy nhất.

Thuật toán trong nghiên cứu có áp dụng được cho dầm có tiết diện thay đổi không?
Có. Thuật toán phần tử hữu hạn cho phép gán các đặc trưng hình học $A$, mômen quán tính $I$ và môđun đàn hồi $E$ riêng biệt cho từng phần tử $e$, do đó hoàn toàn giải quyết chính xác các bài toán dầm có tiết diện biến thiên hoặc làm từ nhiều vật liệu khác nhau.

Kết luận

• Hệ thống hóa hoàn chỉnh 4 phương pháp cơ bản xây dựng bài toán cơ học kết cấu và chứng minh tính tương đương thông qua nguyên lý cực trị Gauss.
• Thiết lập thành công ma trận độ cứng phần tử dầm uốn phẳng $4 \times 4$ và quy trình ghép nối ma trận tổng thể có xét đến điều kiện liên tục góc xoay.
• Kiểm chứng thực nghiệm số trên môi trường Matlab với các mô hình dầm chịu tải phân bố đều, đạt độ chính xác tuyệt đối ở quy mô lưới 4 đến 16 phần tử.
• Đề xuất giải pháp số hóa hiệu quả, giúp tiết kiệm 35% thời gian tính toán và kiểm soát sai số nội lực dưới 0,05% so với lý thuyết giải tích.
• Mở ra hướng phát triển mở rộng thuật toán cho bài toán dầm chịu tải trọng động và vật liệu phi tuyến trong lộ trình nghiên cứu 2-3 năm tới.

Quý độc giả, kỹ sư kết cấu và các nhà nghiên cứu quan tâm có thể khai thác trực tiếp quy trình thuật toán và hệ thống mã nguồn mẫu này để phục vụ công tác thiết kế thực tiễn và phát triển học thuật.