Tổng quan nghiên cứu

Trong bối cảnh đô thị hóa và phát triển hạ tầng hiện đại, các công trình xây dựng ngày càng phát triển về quy mô với chiều cao lớn và nhịp khẩu độ vượt trên 50 mét. Sự xuất hiện phổ biến của các kết cấu thanh mảnh, dầm khẩu độ lớn hay hệ vỏ mỏng đòi hỏi công tác tính toán ổn định và an toàn động lực học phải đạt độ chuẩn xác tuyệt đối. Thực tế ghi nhận có khoảng 85% đến 90% sự cố biến dạng hoặc phá hủy kết cấu nguy hiểm bắt nguồn từ hiện tượng cộng hưởng khi công trình chịu tác động của tải trọng động như gió bão, động đất hoặc xung lực máy móc.

Luận văn thạc sĩ kỹ thuật chuyên ngành Kỹ thuật Xây dựng Công trình Dân dụng và Công nghiệp (mã số 60.08) tại Trường Đại học Dân lập Hải Phòng tập trung giải quyết bài toán cốt lõi: nghiên cứu dao động tự do của thanh đàn hồi bằng lời giải bán giải tích. Đề tài hướng đến mục tiêu xây dựng một quy trình tính toán giải tích số hóa hiệu quả, cho phép xác định chính xác phổ tần số dao động riêng và các dạng dao động chính của thanh với các điều kiện biên khác nhau.

Phạm vi nghiên cứu bao gồm việc khảo sát toàn diện 6 sơ đồ liên kết kết cấu kinh điển trong xây dựng, phân tích từ tần số cơ bản thứ nhất đến các tần số bậc cao. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc cung cấp một công cụ toán học tường minh, giúp rút ngắn thời gian tính toán lặp so với các mô hình phần tử hữu hạn cồng kềnh, đồng thời đạt độ tin cậy vượt trội với sai số kiểm chứng dưới 1% ở các dạng dao động bậc thấp.

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 trên nền tảng cơ học kết cấu vững chắc với sự kết hợp của 3 hệ thống lý thuyết then chốt. Trọng tâm của khung lý thuyết là Nguyên lý cực trị Gauss do nhà toán học Carl Friedrich Gauss đề xuất năm 1829, sau đó được Giáo sư Hà Huy Cương phát triển sáng tạo cho cơ học vật rắn biến dạng. Nguyên lý khẳng định rằng chuyển động thực của cơ hệ sẽ làm cực tiểu hóa lượng cưỡng bức toàn phần.

Bên cạnh đó, tác giả vận dụng Nguyên lý D'Alembert nhằm quy bài toán động lực học về trạng thái cân bằng tĩnh động thông qua việc bổ sung lực quán tính $f_m = -m \ddot{y}$. Lý thuyết uốn dầm Euler-Bernoulli được áp dụng với 2 giả thiết kinh điển: 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, các thớ dọc không chèn ép lẫn nhau.

Về mặt tiêu hao năng lượng, nghiên cứu tổng hợp 3 mô hình cản cơ bản:

  • Mô hình ma sát nhớt Voigt với lực cản tỷ lệ bậc nhất với vận tốc: $P_c = C y'$.
  • Giả thiết cản phi đàn hồi Sorokin biểu diễn tổn thất trễ năng lượng: $P_c = i \gamma P_đ / 2$.
  • Mô hình ma sát khô Coulomb: $F_{ms} = f N$.

Hệ thống lý thuyết này đảm bảo tính tương thích toán học và mô tả chân thực trạng thái ứng suất - biến dạng của thanh dưới các điều kiện chịu tải đa dạng.

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

Nghiên cứu ứng dụng phương pháp bán giải tích kết hợp kỹ thuật chuyển vị cưỡng bức và nhân tử Lagrange. Nguồn dữ liệu phục vụ nghiên cứu là tập mẫu gồm 6 mô hình liên kết dầm chuẩn mực trong công trình dân dụng và công nghiệp: thanh ngàm - tự do, thanh hai đầu khớp, thanh ngàm - khớp, thanh hai đầu ngàm, thanh khớp - tự do và thanh hai đầu tự do. Phương pháp chọn mẫu là chọn mẫu có chủ đích bao phủ 100% các dạng liên kết biên cơ bản của kết cấu xây dựng.

Quy trình giải tích được tiến hành bằng cách xấp xỉ đường độ võng của thanh bằng chuỗi đa thức bậc 9 với 9 hệ số hình học độc lập. Lượng cưỡng bức Gauss được thiết lập dưới dạng tích phân năng lượng uốn và công của lực quán tính. Bằng cách gán một chuyển vị cưỡng bức $y_0$ tại vị trí xác định và đưa vào các điều kiện biên động học, bài toán cực trị có điều kiện được chuyển đổi thành bài toán cực trị tự do thông qua phiếm hàm Lagrange mở rộng.

Lý do lựa chọn phương pháp phân tích này là khả năng triệt tiêu sự phụ thuộc của biến thời gian, biến phương trình vi phân dao động thành hệ đại số tuyến tính gồm 13 đến 45 ẩn số. Toàn bộ thuật toán được lập trình và giải tự động trên công cụ Symbolic Math Toolbox của phần mềm MATLAB trong suốt chu kỳ nghiên cứu 12 tháng, cho phép xuất trực tiếp đa thức đặc trưng bậc 18 đến 25 để trích xuất các trị riêng chính xác.

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

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

Quá trình giải tích đa thức đặc trưng bằng công cụ đại số ký hiệu đã mang lại những kết quả định lượng chuẩn xác về tần số dao động riêng $\omega_i = k_{1i} \sqrt{EJ/ml^4}$ cho từng trường hợp liên kết cụ thể:

  • Thanh ngàm - tự do (Consol): Xác định được 18 nghiệm tần số, trong đó 3 tần số dao động cơ bản đầu tiên đạt giá trị hệ số lần lượt là $k_{11} = 3{,}5160$, $k_{12} = 22{,}0344$ và $k_{13} = 61{,}6968$. Sai số của nghiệm bậc 4 và bậc 5 ghi nhận mức dưới 1%, trong khi nghiệm bậc 6 có sai số dưới 3% so với nghiệm giải tích chính xác.
  • Thanh hai đầu khớp: Giải hệ 25 phương trình đại số thu được 25 trị riêng, với 3 hệ số tần số đầu tiên đạt $k_{11} = 9{,}8695$, $k_{12} = 39{,}4782$ và $k_{13} = 88{,}8300$. Giá trị tần số cơ bản thứ nhất hoàn toàn trùng khớp với hệ số lý thuyết $\pi^2 \approx 9{,}8696$.
  • Thanh đầu ngàm - đầu khớp: Cho 3 hệ số tần số dao động cơ bản tương ứng là $k_{11} = 15{,}418$, $k_{12} = 49{,}964$ và $k_{13} = 104{,}266$.
  • Thanh hai đầu ngàm: Hệ số tần số cơ bản đạt giá trị cao nhất do độ cứng liên kết lớn nhất: $k_{11} = 22{,}373$, $k_{12} = 61{,}672$ và $k_{13} = 120{,}941$.
  • Thanh khớp - tự do và hai đầu tự do: Bằng kỹ thuật đặt thanh trên nền lò xo tương đương với độ cứng $k_{loxo} \to 0$, nghiên cứu đã tìm ra các phổ tần số riêng tương ứng, ví dụ thanh hai đầu tự do đạt $k_{11} = 22{,}373$, $k_{12} = 61{,}672$ và $k_{13} = 120{,}942$.

Thảo luận kết quả

Các phát hiện nghiên cứu khẳng định tính ưu việt vượt bậc của nguyên lý cực trị Gauss khi giải bài toán dao động kết cấu. Nghiệm bán giải tích bằng đa thức bậc 9 cung cấp độ mềm dẻo hình học hoàn hảo, phản ánh chính xác các điểm uốn và nút dao động mà không gặp phải hiện tượng khóa số học.

Trong thực tế tính toán kết cấu, toàn bộ tập dữ liệu phổ tần số có thể được tổ chức trực quan qua bảng tổng hợp so sánh sai số tương đối giữa nghiệm Gauss và nghiệm giải tích chuẩn. Đồng thời, biểu đồ đường cong biến thiên của các dạng dao động riêng từ bậc 1 đến bậc 3 theo trục tọa độ chuẩn hóa $x/l$ từ 0 đến 1 phản ánh rõ nét ảnh hưởng của điều kiện biên: liên kết ngàm làm triệt tiêu hoàn toàn góc xoay ở gốc tọa độ, trong khi liên kết tự do cho phép biên độ chuyển vị đạt giá trị cực đại.

So sánh với phương pháp phần tử hữu hạn, giải pháp bán giải tích này không đòi hỏi chia nhỏ lưới kết cấu thành hàng trăm phần tử để đạt độ hội tụ cho các bậc dao động cao. Việc chuyển đổi bài toán tìm trị riêng ma trận cấp cao về việc giải phương trình đại số của một biến nhân tử Lagrange duy nhất giúp giảm thiểu tới 60% khối lượng bộ nhớ tính toán trên máy tính.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 nhóm giải pháp và khuyến nghị ứng dụng thực tiễn trong công tác thiết kế và nghiên cứu kết cấu công trình:

  • Chuẩn hóa quy trình tính toán bán giải tích thành bộ công cụ tự động hóa: Các viện nghiên cứu kết cấu và trường đại học kỹ thuật cần chủ trì xây dựng các mô-đun thuật toán nguồn mở trong giai đoạn 6 tháng tới, hướng đến mục tiêu giảm khoảng 35% đến 40% thời gian phân tích động lực học sơ bộ cho dầm cầu và sàn khẩu độ lớn.
  • Mở rộng phạm vi áp dụng cho kết cấu tiết diện thay đổi và vật liệu dị tính: Các nhóm nghiên cứu chuyên sâu về cơ học tính toán cần tiếp tục phát triển phiếm hàm Gauss cho các cấu kiện thanh tiết diện hình nón, thanh thép - bê tông liên hợp và vật liệu composite nhiều lớp trong khung thời gian 12 tháng, nâng độ chính xác mô phỏng thực tế lên trên 95%.
  • Tích hợp thuật toán cực trị Gauss vào phần mềm thương mại CAE: Các doanh nghiệp công nghệ xây dựng cần phối hợp với chuyên gia cơ học kết cấu trong lộ trình 18 tháng để nhúng thuật toán vào các phần mềm thiết kế chuyên dụng, hỗ trợ kiểm soát tự động nguy cơ cộng hưởng động học cho hơn 80% dự án nhà cao tầng và nhịp lớn.
  • Hiện đại hóa chương trình đào tạo sau đại học ngành Kỹ thuật Công trình: Các khoa xây dựng tại các trường đại học kỹ thuật cần biên soạn chuyên đề bổ trợ về phương pháp biến phân và nguyên lý Gauss trong vòng 9 tháng, áp dụng trực tiếp cho 100% học viên thạc sĩ chuyên ngành nhằm nâng cao năng lực giải quyết các bài toán cơ học phi tuyến phức tạp.

Đố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ị thực tiễn và học thuật sâu sắc cho 4 nhóm đối tượng trọng tâm:

  • Kỹ sư tư vấn thiết kế kết cấu công trình: Tài liệu cung cấp cơ sở toán học tường minh giúp kỹ sư kiểm toán nhanh tần số dao động riêng của dầm chuyển, dầm sàn nhịp lớn, dầm cầu thép, từ đó đưa ra phương án bố trí kích thước tiết diện tối ưu nhằm tránh xa vùng tần số kích động của gió và động đất.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Xây dựng, Cơ học Kết cấu: Đây là tài liệu tham khảo mẫu mực về quy trình nghiên cứu bán giải tích, hướng dẫn chi tiết cách thiết lập phiếm hàm năng lượng, xử lý nhân tử Lagrange và lập trình giải tích ký hiệu trên phần mềm tính toán chuyên dụng.
  • Giảng viên và nhà nghiên cứu tại các viện, trường đại học khối kỹ thuật: Tài liệu đóng vai trò là nguồn học liệu giá trị phục vụ công tác giảng dạy các môn học Động lực học công trình nâng cao, Phương pháp biến phân trong cơ học và Cơ học vật rắn biến dạng.
  • Chuyên gia phát triển phần mềm mô phỏng kỹ thuật: Cung cấp thuật toán nền tảng và các biểu thức giải tích kiểm chứng chuẩn mực để lập trình phát triển các mô-đun phân tích dao động nhanh trong các hệ thống phần mềm thiết kế kết cấu thế hệ mới.

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

Nguyên lý cực trị Gauss mang lại ưu điểm vượt trội gì so với phương pháp năng lượng truyền thống?
Phương pháp năng lượng cổ điển phụ thuộc nhiều vào việc phỏng đoán dạng đường đàn hồi và thường chỉ cho kết quả chính xác ở tần số cơ bản bậc 1. Nguyên lý cực trị Gauss biến đổi bài toán về việc tìm cực tiểu lượng cưỡng bức toàn phần, kết hợp phương pháp bán giải tích cho phép thiết lập trực tiếp phương trình vi phân và xác định đồng thời toàn bộ phổ tần số bậc cao với độ chính xác trên 97%.

Vai trò của phương pháp chuyển vị cưỡng bức và nhân tử Lagrange trong luận văn là gì?
Khi đặt chuyển vị cưỡng bức $y_0$ tại một điểm, hệ sẽ phát sinh lực giữ $\lambda$ tương ứng với thừa số Lagrange. Về mặt cơ học, lực giữ này đưa vào nhân tạo nên phải triệt tiêu ($\lambda = 0$). Phương trình đại số $\lambda(\omega) = 0$ chính là phương trình đặc trưng của hệ, giúp tìm ra trực tiếp các tần số dao động riêng mà không cần lập ma trận độ cứng phức tạp.

Việc sử dụng đa thức bậc 9 có gây hiện tượng sai số tích lũy khi tính các dạng dao động cao không?
Xấp xỉ đa thức bậc 9 cung cấp đủ bậc tự do hình học để mô tả chính xác các biến dạng uốn phức tạp. Kết quả kiểm chứng cho thấy 3 tần số cơ bản đầu tiên đạt độ chính xác gần như tuyệt đối, sai số ở tần số bậc 4 và 5 dưới 1%, bậc 6 dưới 3%. Muốn tăng độ chính xác ở các bậc cao hơn, người dùng có thể tăng bậc đa thức hoặc chia thanh thành 2 đoạn tính toán.

Kết quả của luận văn có thể áp dụng cho các cấu kiện dầm chịu lực dọc trục đồng thời không?
Hoàn toàn có thể. Khung lý thuyết của luận văn đã tích hợp thế năng biến dạng nén uốn khi có lực dọc trục $P$. Khi lực dọc thay đổi tuần hoàn theo thời gian, phương trình chuyển động sẽ chuyển thành phương trình dao động tham số (dạng Mathieu-Hill), cho phép đánh giá cả miền mất ổn định động lực học của thanh.

Tại sao cần xem xét các mô hình lực cản khác nhau trong bài toán dao động công trình?
Lực cản đóng vai trò quyết định đến việc giới hạn biên độ chuyển vị khi xảy ra hiện tượng cộng hưởng. Trong khi mô hình cản nhớt Voigt phù hợp với kết cấu làm việc trong môi trường chất lỏng hoặc có thiết bị giảm chấn, thì mô hình Sorokin phản ánh chính xác sự tiêu hao năng lượng nội tại của vật liệu bê tông và kim loại trong miền biến dạng đàn hồi.

Kết luận

  • Luận văn đã ứng dụng thành công Nguyên lý cực trị Gauss kết hợp phương pháp chuyển vị cưỡng bức để giải quyết trọn vẹn bài toán dao động tự do của thanh đàn hồi.
  • Nghiên cứu đã xác lập các công thức bán giải tích tường minh cho 6 sơ đồ liên kết kết cấu kinh điển, tìm ra chính xác phổ tần số dao động và các hàm dạng dao động riêng tương ứng.
  • Phương pháp đề xuất giúp chuyển đổi bài toán trị riêng phức tạp về phương trình đại số của nhân tử Lagrange, tiết kiệm tới 60% tài nguyên tính toán so với phương pháp số truyền thống.
  • Các kết quả định lượng đạt độ tin cậy vượt trội với sai số dưới 1% ở các bậc dao động thấp, tạo cơ sở khoa học vững chắc cho việc kiểm toán cộng hưởng công trình.
  • Lộ trình phát triển tiếp theo trong 12 đến 24 tháng tới tập trung vào việc mở rộng thuật toán cho hệ thanh không gian, thanh tiết diện thay đổi và tích hợp vào các phần mềm kỹ thuật số hóa hiện đại.

Các kỹ sư kết cấu, nhà nghiên cứu và học viên cao học quan tâm có thể khai thác trực tiếp các biểu thức toán học và quy trình thuật toán trong luận văn để ứng dụng vào công tác thiết kế, giảng dạy và phát triển các công cụ phân tích động lực học kết cấu chuyên sâu.