Tổng quan nghiên cứu

Trong xu hướng phát triển hạ tầng hiện đại, số lượng các công trình cao tầng từ 20 đến trên 50 tầng cùng các kết cấu vượt nhịp khẩu độ lớn trên 30 mét gia tăng nhanh chóng, chiếm hơn 45% tổng khối lượng xây dựng công trình đặc biệt. Thực tế này đặt ra yêu cầu cấp thiết về việc kiểm soát ổn định và ứng xử động lực học của hệ thanh chịu lực trong miền đàn hồi nhằm giảm thiểu nguy cơ cộng hưởng phá hoại tới hơn 40%. Luận văn thạc sĩ 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ố: 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 dao động đàn hồi của thanh thẳng chịu tác dụng của các điều kiện biên và tải trọng khác nhau.

Mục tiêu cốt lõi của đề tài là xây dựng quy trình giải tích số chính xác dựa trên việc kết hợp nguyên lý cực trị Gauss với phương pháp chuyển vị cưỡng bức và kỹ thuật nhân tử Lagrange. Nghiên cứu tập trung phân tích phổ tần số dao động riêng, xác định các dạng dao động chính và đánh giá sự phân bố chuyển vị dưới tác động của lực quán tính. Phạm vi nghiên cứu bao quát các dạng liên kết cơ bản trong kỹ thuật công trình, được thực nghiệm tính toán trên công cụ tính toán biểu thức đại số trong giai đoạn 2016-2017.

Ý 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 trực quan, cho phép loại bỏ các bước lập ma trận độ cứng phức tạp trong phương pháp phần tử hữu hạn truyền thống. Giải pháp này giúp rút ngắn thời gian tính toán sơ bộ khoảng 35% mà vẫn đảm bảo độ chính xác tuyệt đối ở 3 tần số dao động cơ bản đầu tiên, hỗ trợ đắc lực cho các kỹ sư kết cấu trong việc tối ưu hóa tiết diện và đảm bảo an toàn công trình trước các nguồn kích động động lực học.

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 vật rắn biến dạng và động lực học công trình với hai lý thuyết trụ cột:

  1. Nguyên lý cực trị Gauss cho môi trường liên tục: Mở rộng từ nguyên lý cơ học chất điểm phát biểu năm 1829, lý thuyết này xem xét chuyển động thực của hệ biến dạng đàn hồi là trạng thái làm cực tiểu hóa độ đo cưỡng bức Z. Độ cứng mặt cắt uốn EJx được đồng nhất tương đương với khối lượng, mômen uốn Mx đóng vai trò là lực tác dụng và đạo hàm bậc hai của độ võng theo tọa độ dọc trục tương ứng với gia tốc chuyển động.
  2. Lý thuyết dầm Euler-Bernoulli và nguyên lý D'Alembert: Kết cấu dầm tuân thủ hai giả thiết cơ bản gồm giả thiết mặt cắt ngang phẳng trước và sau biến dạng, cùng giả thiết các thớ dọc không ép lẫn nhau. Nguyên lý D'Alembert đưa lực quán tính phân bố vào phương trình vi phân cân bằng dưới dạng ngoại lực cản động học.

Mô hình nghiên cứu định nghĩa rõ 5 khái niệm chuyên ngành then chốt: tần số dao động riêng (đặc trưng cho tốc độ dao động tự nhiên của hệ), dạng dao động chính (vectơ độ võng trực giao chuẩn hóa), lượng cưỡng bức Z (phiếm hàm năng lượng sai lệch), nhân tử Lagrange (đại diện cho phản lực giữ tại điểm đặt chuyển vị cưỡng bức), và mô hình cản đàn nhớt Voigt (lực cản tỷ lệ bậc nhất với vận tốc).

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

Nguồn dữ liệu tính toán được thiết lập từ mô hình giải tích giải phóng liên kết và thiết lập hệ dầm so sánh tương đương. Nghiên cứu lựa chọn cỡ mẫu gồm 5 mô hình liên kết biên chuẩn tắc đại diện cho toàn bộ các trường hợp chịu lực điển hình trong xây dựng: ngàm - tự do, khớp - khớp, ngàm - khớp, ngàm - ngàm, và khớp - tự do đặt trên nền đàn hồi lò xo.

Phương pháp chọn mẫu có chủ đích được áp dụng nhằm bao phủ toàn diện từ các hệ tĩnh định đơn giản đến các hệ siêu tĩnh bậc cao và hệ có nguy cơ biến hình hình học. Phương pháp phân tích sử dụng công cụ toán học giải tích ký hiệu Matlab Symbolic để biến phân phiếm hàm Gauss mở rộng, xấp xỉ đường đàn hồi bằng chuỗi đa thức bậc 9 với 9 đến 18 thông số hình học độc lập.

Lý do lựa chọn phương pháp xấp xỉ đa thức kết hợp chuyển vị cưỡng bức là khả năng chuyển bài toán động phức tạp về bài toán tĩnh tương đương, thiết lập trực tiếp phương trình đại số chứa từ 13 đến 25 ẩn số mà không cần lập phương trình vi phân đạo hàm riêng theo thời gian. Toàn bộ quy trình mô phỏng và kiểm chứng nghiệm đại số được hoàn thiện trong timeline nghiên cứu kéo dài 12 tháng.

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

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

Nghiên cứu đã dẫn xuất thành công đa thức đặc trưng bậc cao và xác định chính xác phổ tần số dao động riêng cho từng trường hợp liên kết biên cụ thể:

  • Hệ thanh một đầu ngàm - một đầu tự do: Thiết lập hệ 13 phương trình đại số, tìm ra 18 nghiệm tần số. Ba hệ số tần số dao động cơ bản đầu tiên đạt giá trị lần lượt là 3,5160; 22,0344 và 61,6968 nhân với căn bậc hai của tỷ số giữa độ cứng uốn và khối lượng phân bố trên chiều dài bậc 4. Sai số của dạng dao động thứ 4 và 5 nhỏ hơn 1%, dạng thứ 6 nhỏ hơn 3%.
  • Hệ thanh hai đầu khớp: Thiết lập hệ 25 phương trình đại số với 6 điều kiện biên động học và tĩnh học, cho ra các hệ số tần số chính xác 9,8695; 39,4782 và 88,8300 (trùng khớp hoàn toàn với các bội số của số Pi bình phương trong lời giải giải tích kinh điển).
  • Hệ thanh đầu ngàm - đầu khớp: Xác định 3 tần số riêng cơ bản có hệ số tương ứng là 15,418; 49,964 và 104,266. So với hệ hai đầu khớp, tần số cơ bản tăng khoảng 56,2% do liên kết ngàm làm tăng độ cứng chống uốn tổng thể.
  • Hệ thanh hai đầu ngàm: Cho giá trị hệ số tần số cơ bản cao nhất đạt 22,373; các bậc tiếp theo đạt 61,672 và 120,941. Tần số bậc 1 của hệ hai đầu ngàm cao hơn 536% so với hệ công xôn ngàm - tự do và cao hơn 126% so với hệ hai đầu khớp.
  • Hệ thanh đầu khớp - đầu tự do trên nền đàn hồi: Khi cho độ cứng nền lò xo triệt tiêu dần về 0, hệ số tần số riêng ổn định tại các mức 15,418; 49,964 và 104,253.

Thảo luận kết quả

Các kết quả định lượng khẳng định quy luật cơ học: mức độ liên kết biên càng chặt thì độ cứng động lực học của hệ càng lớn, dẫn đến phổ tần số dao động riêng dịch chuyển mạnh về phía miền tần số cao. Điều này giải thích nguyên nhân các công trình có liên kết ngàm cứng hai đầu ít bị ảnh hưởng bởi các tải trọng gió dao động tần số thấp so với kết cấu công xôn.

Về phương diện trực quan hóa, toàn bộ tập dữ liệu nghiên cứu có thể được thể hiện hiệu quả qua biểu đồ cột so sánh tương quan hệ số tần số bậc 1 đến bậc 3 giữa 5 loại liên kết, kết hợp với bảng tổng hợp sai số tương đối giữa nghiệm đa thức bậc 9 và nghiệm giải tích chuẩn. Đồng thời, đồ thị đường cong phân bố độ võng chuẩn hóa dọc theo trục thanh (tỷ số x trên l từ 0 đến 1) phản ánh chính xác các điểm nút và điểm bụng dao động của từng dạng riêng.

So sánh với phương pháp phần tử hữu hạn thông thường, phương pháp nguyên lý cực trị Gauss kết hợp nhân tử Lagrange không yêu cầu rời rạc hóa kết cấu thành hàng trăm phần tử vi mô mà vẫn đạt độ chính xác tương đương trên 99% cho các dạng dao động chính đầu tiên.

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

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

  1. Ứng dụng thuật toán Gauss - Lagrange vào phần mềm thiết kế kết cấu: Tích hợp module giải phương trình đại số tự động dựa trên phiếm hàm Gauss vào các phần mềm tính toán kết cấu nội địa. Mục tiêu rút ngắn 35% thời gian phân tích phổ động lực học giai đoạn thiết kế cơ sở. Thời gian triển khai trong vòng 6 đến 12 tháng, do các kỹ sư R&D phần mềm xây dựng chủ trì thực hiện.
  2. Chuẩn hóa quy trình kiểm tra cộng hưởng cho thanh nhịp lớn: Ban hành hướng dẫn tính toán kiểm tra độ lệch giữa tần số kích thích cưỡng bức và tần số dao động riêng của dầm cầu, thanh dàn khẩu độ trên 25 mét, đảm bảo khoảng cách an toàn tối thiểu 20%. Thời gian áp dụng trong 3 tháng, do các đơn vị tư vấn thiết kế và thẩm tra kỹ thuật đảm nhiệm.
  3. Nâng cấp bậc đa thức xấp xỉ và chia đoạn kết cấu: Tiếp tục phát triển giải thuật với đa thức bậc 12 đến 15 hoặc kỹ thuật chia 2 đến 4 đoạn dầm để kiểm soát sai số các tần số bậc cao (từ bậc 6 trở lên) dưới mức 0,5%. Kế hoạch nghiên cứu 6 tháng, do các nhóm nghiên cứu cơ học tại các trường đại học kỹ thuật thực hiện.
  4. Mở rộng mô hình tính toán có xét lực nén dọc trục và cản phi đàn hồi: Đưa thành phần lực dọc biến thiên theo thời gian vào phương trình cân bằng năng lượng để dự báo chính xác hiện tượng mất ổn định động lực học với độ tin cậy trên 95%. Thời gian hoàn thiện mô hình từ 12 đến 18 tháng, phối hợp giữa viện nghiên cứu chuyên ngành và các chuyên gia cơ học công trình.

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

Luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng chuyên môn:

  • Kỹ sư thiết kế kết cấu công trình: Tiếp cận công cụ tính nhanh tần số dao động riêng và dạng dao động chính của các cấu kiện thanh dầm vượt nhịp, dầm công xôn nhà cao tầng, phục vụ tính toán kháng chấn và chống rung gió.
  • Giảng viên và nghiên cứu sinh ngành Cơ kỹ thuật, Kỹ thuật Xây dựng: Sử dụng làm tài liệu giảng dạy chuyên đề Động lực học công trình và Cơ học vật rắn biến dạng, minh họa phương pháp biến phân phiếm hàm Gauss và kỹ thuật nhân tử Lagrange.
  • Học viên cao học và sinh viên kỹ thuật: Nắm bắt phương pháp luận nghiên cứu khoa học, cách sử dụng công cụ tính toán biểu thức đại số Matlab Symbolic để giải quyết hệ phương trình đại số tuyến tính nhiều ẩn số.
  • Chuyên gia phát triển phần mềm mô phỏng (CAE/FEA): Khai thác thuật toán tối ưu hóa thông số không qua ma trận khối lượng nghịch đảo để phát triển các thuật toán giải nhanh trong phân tích kết cấu cơ khí và công trình.

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

  1. Phương pháp nguyên lý cực trị Gauss có điểm gì ưu việt hơn so với phương pháp phần tử hữu hạn truyền thống? Phương pháp Gauss biến phân trực tiếp trên phiếm hàm lượng cưỡng bức của hệ liên tục, xem mômen uốn và biến dạng là các đại lượng độc lập. Cách tiếp cận này loại bỏ hoàn toàn việc thiết lập và khả nghịch các ma trận độ cứng, ma trận khối lượng cấp cao, giúp giảm thiểu khối lượng tính toán trung gian khoảng 30% đến 40% trong các bài toán đơn thanh.

  2. Tại sao việc đưa chuyển vị cưỡng bức vào lại giúp tìm được tần số dao động riêng của thanh? Khi áp đặt một chuyển vị cưỡng bức tại một điểm trên thanh, nhân tử Lagrange xuất hiện trong phiếm hàm mở rộng đại diện cho lực giữ để duy trì chuyển vị đó. Vì dao động tự do là dao động không có ngoại lực duy trì, điều kiện lực giữ bằng 0 dẫn trực tiếp đến phương trình đa thức đặc trưng xác định các trị riêng tần số dao động.

  3. Độ chính xác của các tần số dao động riêng khi dùng đa thức bậc 9 đạt mức bao nhiêu? Nghiên cứu chứng minh 3 tần số dao động riêng đầu tiên hoàn toàn trùng khớp với nghiệm giải tích chính xác. Đối với các tần số bậc cao hơn, nghiệm bậc 4 và bậc 5 có sai số nhỏ hơn 1%, nghiệm bậc 6 có sai số dưới 3%, đáp ứng vượt trội các tiêu chuẩn dung sai kỹ thuật xây dựng hiện hành.

  4. Phương pháp này có thể áp dụng cho thanh có tiết diện thay đổi hoặc chịu tải trọng nén uốn không? Hoàn toàn có thể áp dụng. Nghiên cứu đã thiết lập cơ sở lý thuyết mở rộng cho thanh chịu lực dọc trục phân bố hoặc tập trung. Khi độ cứng uốn và khối lượng phân bố biến thiên theo tọa độ trục, phiếm hàm Gauss vẫn giữ nguyên dạng tích phân và có thể giải số dễ dàng.

  5. Vì sao trường hợp thanh đầu khớp - đầu tự do lại cần đưa vào mô hình nền đàn hồi lò xo? Thanh một đầu khớp một đầu tự do là hệ biến hình hình học trong cơ học tĩnh. Việc đưa nền đàn hồi lò xo với độ cứng kloxo vào nhằm tạo liên kết duy trì tính ổn định của phiếm hàm. Khi giải phương trình đa thức đặc trưng và cho độ cứng lò xo tiến về 0, ta thu được chính xác nghiệm dao động riêng thực tế của hệ.

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 đàn hồi của hệ thanh.
  • Xây dựng thành công quy trình lập phiếm hàm lượng cưỡng bức mở rộng và hệ phương trình Euler xác định 18 đến 25 thông số hình học cho 5 dạng liên kết biên cơ bản.
  • Đạt độ chính xác tuyệt đối ở 3 tần số dao động cơ bản đầu tiên với sai số các bậc dao động cao luôn được kiểm soát dưới ngưỡng 3%.
  • Cung cấp giải pháp tính toán giải tích số tường minh, giảm tải áp lực tính toán ma trận lớn trên máy tính cá nhân.
  • Lộ trình 6 đến 12 tháng tới hướng đến việc tích hợp thuật toán vào các phần mềm thiết kế chuyên dụng và mở rộng cho hệ kết cấu khung không gian phức tạp.

Các kỹ sư và nhà nghiên cứu quan tâm có thể ứng dụng ngay khung thuật toán Gauss - Lagrange này để tối ưu hóa thiết kế kháng chấn cho các cấu kiện thanh dầm nhịp lớn trong công trình thực tế.