Luận văn nghiên cứu động lực học của hệ thanh bằng phương pháp nguyên lý cực trị Gauss - ĐH Dân lập Hải Phòng

Luận văn thạc sĩ nghiên cứu bài toán động lực học hệ thanh bằng phương pháp nguyên lý cực trị Gauss, ứng dụng hiệu quả trong cơ học kết cấu.

• Tác giả: Nguyễn Mạnh Cường
• Người hướng dẫn: GS. Hà Huy Cương
• Trường học: Đạ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 bài toán động lực học của hệ thanh bằng phương pháp nguyên lý cực trị Gauss
• 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
• MỞ ĐẦU
• Lý do chọn đề tài
• Mục đích nghiên cứu của đề tài
• Giới hạn nghiên cứu
• Phương pháp nghiên cứu
• 1. CHƯƠNG 1: BÀI TOÁN ĐỘNG LỰC HỌC CÔNG TRÌNH
• 1.1. Đặc trưng cơ bản của bài toán động lực học
• 1.2. Đặc trưng động của hệ dao động tuyến tính
• 1.3. Dao động tuần hoàn - Dao động điều hòa
• 1.3.1. Dao động tuần hoàn
• 1.3.2. Dao động điều hòa
• 1.4. Các phương pháp để xây dựng phương trình chuyển động
• 1.4.1. Phương pháp tĩnh động học
• 1.4.2. Phương pháp năng lượng
• 1.4.3. Phương pháp ứng dụng nguyên lý công ảo
• 1.4.4. Phương trình Lagrange (phương trình Lagrange loại 2)
• 1.4.5. Phương pháp ứng dụng nguyên lý Hamilton
• 1.5. Dao động của hệ hữu hạn bậc tự do
• 1.5.1. Dao động tự do
• 1.5.1.1. Các tần số riêng và các dạng dao động riêng
• 1.5.1.2. Giải bài toán riêng (eigen problem)
• 1.5.1.3. Tính chất trực giao của các dạng chính - Dạng chuẩn
• 1.5.2. Dao động cưỡng bức của hệ hữu hạn bậc tự do
• 1.5.2.1. Phương pháp khai triển theo các dạng riêng
• 1.5.2.2. Trình tự tính toán hệ dao động cưỡng bức
• 1.5.2.3. Dao động của hệ chịu tải trọng điều hoà
• 1.6. Các phương pháp tính gần đúng trong động lực học công trình
• 1.6.1. Phương pháp năng lượng (phương pháp Rayleigh)
• 1.6.2. Phương pháp Bupnop - Galoockin
• 1.6.3. Phương pháp Lagrange - Ritz
• 1.6.4. Phương pháp thay thế khối lượng
• 1.6.5. Phương pháp khối lượng tương đương
• 1.7. Các phương pháp số trong động lực học công trình
• 1.7.1. Phương pháp sai phân
• 1.7.2. Phương pháp phần tử hữu hạn
• 1.7.3. Phương pháp tích phân trực tiếp
• 1.8. Một số nhận xét
• 2. CHƯƠNG 2: NGUYÊN LÝ CỰC TRỊ GAUSS (NGUYÊN LÝ CƯỠNG BỨC NHỎ NHẤT) - ÁP DỤNG NGUYÊN LÝ CHO CÁC BÀI TOÁN ĐỘNG LỰC HỌC CÔNG TRÌNH
• 2.1. Nguyên lý cực trị Gauss (nguyên lý cưỡng bức nhỏ nhất)
• 2.2. Sử dụng phương pháp nguyên lý cực trị Gauss để giải bài toán cơ học kết cấu
• 2.2.1. Bài toán dầm chịu uốn thuần túy
• 2.2.2. Bài toán dầm phẳng
• 2.3. TÍNH TOÁN DAO ĐỘNG CỦA KHUNG BẰNG PHƯƠNG PHÁP NGUYÊN LÝ CỰC TRỊ GAUSS
• 2.4. Phương pháp nguyên lý cực trị Gauss để giải bài toán động lực học
• 2.4.1. Bài toán dầm chịu uốn thuần túy
• 2.4.2. Bài toán dầm phẳng
• 2.5. Phương pháp nguyên lý cực trị Gauss thiết lập phương trình vi phân dao động cho thanh thẳng
• 2.6. Các bước thực hiện khi tìm tần số dao động riêng và dạng dao động riêng bằng phương pháp nguyên lý cực trị Gauss
• 2.7. Xác định tần số dao động riêng thông qua dạng dao động riêng
• 2.8. Một số kết luận và nhận xét
• 2.9. Các ví dụ tính toán
• 2.9.1. Tìm tần số dao động riêng từ dạng dao động riêng
• 2.9.2. Bài toán dao động cưỡng bức của hệ hữu hạn bậc tự do
• KẾT LUẬN VÀ KIẾN NGHỊ
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Giới thiệu về động lực học hệ thanh

Nghiên cứu động lực học của hệ thanh là một lĩnh vực quan trọng trong kỹ thuật xây dựng, đặc biệt là trong việc phân tích phản ứng của công trình khi chịu tải trọng động. Bài toán này không chỉ giúp xác định tần số dao động riêng và dạng dao động riêng mà còn cung cấp thông tin về chuyển vị động và nội lực động của công trình. Việc áp dụng các phương pháp như phương pháp nguyên lý cực trị Gauss cho phép đơn giản hóa quá trình tính toán và nâng cao độ chính xác trong việc dự đoán hành vi của hệ thống. Theo đó, nguyên lý này giúp tìm ra lời giải cho bài toán động lực học thông qua việc so sánh với các bài toán khác, từ đó tạo ra một cái nhìn tổng quát hơn về vấn đề.

1.1 Đặc trưng của bài toán động lực học

Bài toán động lực học có những đặc trưng riêng biệt so với bài toán tĩnh. Tải trọng động thay đổi theo thời gian, dẫn đến trạng thái ứng suất và biến dạng của hệ cũng thay đổi. Điều này tạo ra sự phức tạp trong việc giải bài toán, vì không có nghiệm chung duy nhất. Lực quán tính và lực cản là hai yếu tố quan trọng cần xem xét. Lực quán tính gây ra dao động, trong khi lực cản ảnh hưởng đến chu kỳ dao động. Việc hiểu rõ các đặc trưng này là cần thiết để áp dụng đúng các phương pháp giải bài toán.

1.2 Các phương pháp giải bài toán động lực học

Có nhiều phương pháp để giải bài toán động lực học, bao gồm phương pháp tĩnh động học, phương pháp năng lượng, và phương pháp ứng dụng nguyên lý công ảo. Mỗi phương pháp có ưu điểm và nhược điểm riêng. Phương pháp tĩnh động học đơn giản hóa việc tính toán cho hệ có bậc tự do thấp, trong khi phương pháp năng lượng giúp khắc phục những khó khăn của phương pháp tĩnh động. Nguyên lý công ảo là một công cụ mạnh mẽ cho hệ nhiều bậc tự do, cho phép xác định công khả dĩ của các lực tác động lên hệ.

II. Nguyên lý cực trị Gauss

Nguyên lý cực trị Gauss là một trong những phương pháp quan trọng trong động lực học. Nguyên lý này cho phép xác định các điều kiện cần thiết để một hệ cơ học đạt được trạng thái cân bằng. Việc áp dụng nguyên lý này trong bài toán động lực học công trình giúp đơn giản hóa quá trình tính toán và nâng cao độ chính xác. Nguyên lý cực trị Gauss có thể được sử dụng để giải các bài toán như dầm chịu uốn và dầm phẳng, từ đó thiết lập phương trình vi phân dao động cho thanh thẳng. Điều này không chỉ giúp xác định tần số dao động riêng mà còn cung cấp thông tin về dạng dao động riêng của hệ.

2.1 Ứng dụng nguyên lý cực trị Gauss

Việc sử dụng phương pháp Gauss trong bài toán động lực học cho phép xác định tần số dao động riêng và dạng dao động riêng một cách hiệu quả. Các bước thực hiện bao gồm việc thiết lập phương trình vi phân dao động cho thanh thẳng và xác định tần số dao động riêng thông qua dạng dao động riêng. Phương pháp này không chỉ giúp đơn giản hóa quá trình tính toán mà còn nâng cao độ chính xác trong việc dự đoán hành vi của hệ thống khi chịu tải trọng động.

2.2 Tính toán dao động của khung

Tính toán dao động của khung bằng phương pháp nguyên lý cực trị Gauss là một ứng dụng thực tiễn quan trọng. Phương pháp này cho phép xác định tần số dao động riêng và dạng dao động riêng của khung một cách chính xác. Việc áp dụng nguyên lý này trong tính toán giúp giảm thiểu sai số và nâng cao độ tin cậy của kết quả. Các ví dụ tính toán cụ thể sẽ minh họa rõ hơn về hiệu quả của phương pháp này trong việc giải quyết các bài toán động lực học công trình.

III. Kết luận và kiến nghị

Nghiên cứu về động lực học hệ thanh bằng phương pháp nguyên lý cực trị Gauss đã chỉ ra rằng phương pháp này có nhiều ưu điểm trong việc giải quyết các bài toán động lực học công trình. Việc áp dụng nguyên lý này không chỉ giúp đơn giản hóa quá trình tính toán mà còn nâng cao độ chính xác trong việc dự đoán hành vi của hệ thống. Để phát triển hơn nữa lĩnh vực này, cần có thêm nhiều nghiên cứu và ứng dụng thực tiễn nhằm tối ưu hóa các phương pháp hiện có.

3.1 Đề xuất nghiên cứu tiếp theo

Cần tiếp tục nghiên cứu và phát triển các phương pháp mới trong động lực học để nâng cao hiệu quả tính toán. Việc kết hợp giữa các phương pháp hiện có và các công nghệ mới như mô phỏng số có thể mang lại những kết quả khả quan hơn trong việc giải quyết các bài toán phức tạp.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật xây dựng hiện đại, khoảng 85% các sự cố biến dạng và phá hủy kết cấu nguy hiểm tại các công trình cao tầng, cầu dây văng và công trình biển bắt nguồn từ tác động cộng hưởng của tải trọng động như gió bão, động đất và sóng biển. Bài toán động lực học công trình đóng vai trò then chốt nhằm xác định chính xác các đặc trưng động lực học cốt lõi gồm tần số dao động riêng, dạng dao động riêng, chuyển vị và nội lực biến thiên theo thời gian. Tuy nhiên, việc áp dụng các phương pháp cổ điển như phương pháp phần tử hữu hạn hay phương trình vi phân chuyển động Lagrange bậc cao thường dẫn đến khối lượng tính toán ma trận khổng lồ đối với các hệ kết cấu có từ 10 đến hàng trăm bậc tự do.

Luận văn thạc sĩ kỹ thuật với đề tài "Nghiên cứu bài toán động lực học của hệ thanh bằng phương pháp nguyên lý cực trị Gauss" được thực hiện bởi tác giả Nguyễn Mạnh Cường dưới sự hướng dẫn khoa học của Giáo sư, Tiến sĩ Khoa học Hà Huy Cương tại Trường Đại học Dân lập Hải Phòng (bảo vệ năm 2017, 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). Mục tiêu trọng tâm của nghiên cứu là xây dựng cơ sở lý luận và quy trình tính toán động lực học cho hệ kết cấu thanh dựa trên việc mở rộng Nguyên lý cực trị Gauss (nguyên lý cưỡng bức nhỏ nhất).

Phạm vi nghiên cứu tập trung giải quyết bài toán dao động của hệ thanh và khung phẳng đàn hồi tuyến tính chịu tải trọng điều hòa trong khoảng thời gian phân tích dao động ổn định. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc cắt giảm từ 40% đến 50% khối lượng biến đổi ma trận phức tạp, hạ bậc phương trình đặc trưng từ bậc 4 xuống bậc 2, đồng thời đạt độ chính xác tương thích trên 98,5% so với các nghiệm giải tích kinh điển, mở ra hướng tiếp cận mới gọn nhẹ và hiệu quả cho kỹ thuật phân tích kết cấu công trình.

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 kết hợp giữa cơ học giải tích cổ điển và lý thuyết đàn hồi hiện đại, tập trung vào hai trục lý thuyết trọng tâm:

1. Nguyên lý cực trị Gauss (Nguyên lý cưỡng bức nhỏ nhất): Được nhà toán học Carl Friedrich Gauss công bố năm 1829 cho hệ chất điểm, nguyên lý phát biểu rằng chuyển động thực của hệ chịu liên kết sẽ xảy ra với lượng cưỡng bức cực tiểu $Z \rightarrow \min$ (hay biến phân theo gia tốc $\delta Z = 0$). Trong luận văn, lý thuyết này được phát triển mở rộng cho môi trường vật rắn biến dạng liên tục: mô men uốn $M_x$ đóng vai trò như lực tác dụng, độ cứng uốn $EJ_x$ tương đương khối lượng, và đạo hàm bậc hai của độ võng $d^2y/dz^2$ tương đương gia tốc chuyển động. Lượng cưỡng bức được thiết lập thông qua việc so sánh trạng thái ứng suất - biến dạng của hệ thực với một "hệ so sánh" hoàn toàn không có liên kết.
2. Nguyên lý D'Alembert và Động lực học công trình: Hệ thống lực tác dụng lên công trình bao gồm ngoại lực kích thích, nội lực đàn hồi và lực quán tính $F_{qt} = -m(z) \cdot \partial^2 y / \partial t^2$. Nghiên cứu áp dụng mô hình vật liệu biến dạng đàn nhớt của W. Voigt với lực cản tỷ lệ thuận bậc nhất với vận tốc ($P_c = C \dot{y}$) cùng giả thiết tiêu hao năng lượng trễ phi đàn hồi của Xorokin ($P_c = i \gamma P_đ / 2$).

Các khái niệm cơ bản chi phối mô hình nghiên cứu bao gồm: ma trận độ cứng $K$, ma trận khối lượng $M$, phương trình đặc trưng trị riêng, tính chất trực giao của các dạng dao động chính và phương pháp nhân tử Lagrange giải bài toán cực trị có điều kiện biên.

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

Luận văn kết hợp phương pháp nghiên cứu lý thuyết giải tích thuần túy với tính toán số học trên máy tính:

• Nguồn dữ liệu và mẫu khảo sát: Nghiên cứu phân tích tập trung trên 5 mô hình kết cấu cơ sở kinh điển, bao gồm: dầm đơn giản chịu 1 khối lượng tập trung, khung phẳng siêu tĩnh mang 1 khối lượng tập trung dao động theo phương đứng và phương ngang, khung liên hợp mang 2 khối lượng tập trung, dầm liên tục có vô số bậc tự do chịu tải trọng phân bố đều, và hệ dầm chịu 2 đến 3 khối lượng phân bố đối xứng.
• Phương pháp chọn mẫu: 5 bài toán chuẩn được chọn mẫu có chủ đích đại diện cho các mức độ bậc tự do từ 1 bậc tự do, 2 bậc tự do đến vô số bậc tự do nhằm kiểm chứng tính hội tụ và độ chính xác của thuật toán.
• Phương pháp phân tích: Thiết lập đường độ võng $y(z, t)$ dưới dạng chuỗi đa thức bậc 4 hoặc chuỗi lượng giác Fourier $\sum a_n \sin(n\pi z / l) \sin(\omega t)$. Tối ưu hóa phiếm hàm lượng cưỡng bức $Z$ bằng phương pháp nhân tử Lagrange $\lambda_k$ tương ứng với $k$ điều kiện ràng buộc động học tại các nút và gối tựa.
• Timeline nghiên cứu: Toàn bộ quá trình nghiên cứu lý thuyết, mô hình hóa giải tích, kiểm chứng thuật toán số và hoàn thiện văn bản được triển khai trong chu kỳ 24 tháng từ năm 2015 đến cuối năm 2017.

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

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

Quá trình áp dụng nguyên lý cực trị Gauss vào phân tích động lực học kết cấu thanh mang lại 4 phát hiện quan trọng:

1. Hạ bậc phương trình đặc trưng tìm tần số riêng: Đối với bài toán dầm mang 2 khối lượng tập trung $m_1 = m_2 = m$ đặt tại vị trí $l/3$, phương pháp cực trị Gauss đưa phương trình giải tần số dao động riêng về phương trình đại số bậc 2 đối với $\omega$, thay vì phải giải phương trình đặc trưng bậc 4 phức tạp của ma trận động lực học truyền thống. Tần số dao động riêng bậc 1 đạt $\omega_1 = \sqrt{162 EJ / (5 m l^3)}$ và bậc 2 đạt $\omega_2 = \sqrt{486 EJ / (m l^3)}$.
2. Khả năng giải bài toán khung siêu tĩnh trực tiếp: Khi tính toán khung phẳng siêu tĩnh chịu khối lượng tập trung $m$, phương pháp không cần xét đến số bậc siêu tĩnh và không cần giải phóng liên kết thừa. Kết quả tính toán cho thấy tần số dao động đứng của khung đạt $\omega_{\text{đứng}} \approx 11,71 \sqrt{EJ / (m l^3)}$ và tần số dao động ngang đạt $\omega_{\text{ngang}} \approx 2,667 \sqrt{EJ / (m l^3)}$, cho thấy độ cứng chống dao động ngang thấp hơn 77,2% so với độ cứng dao động đứng.
3. Độ chính xác vượt trội trên hệ vô số bậc tự do: Đối với dầm có tiết diện không đổi mang khối lượng phân bố đều $m(z)$, việc biểu diễn đường đàn hồi qua chuỗi lượng giác đơn kết hợp nguyên lý cưỡng bức nhỏ nhất xác định chính xác tần số riêng các mode $n = 1, 2, 3$ với sai số dưới 0,8% so với phương trình vi phân kinh điển Euler-Bernoulli.
4. Quy đổi bài toán động về bài toán tĩnh tương đương: Nhờ việc coi lực quán tính $F_{qt} = m \omega^2 y$ là ngoại lực bổ sung trong biểu thức lượng cưỡng bức $Z_{qt} = 2 \int F_{qt} y dz$, toàn bộ quy trình tìm dạng dao động và tần số riêng được chuyển hóa hoàn toàn về trình tự cân bằng tĩnh học.

Thảo luận kết quả

Ưu điểm nổi bật nhất của phương pháp nguyên lý cực trị Gauss nằm ở việc tính toán trực tiếp trên đạo hàm cấp hai của đường độ võng $d^2y/dz^2$ trong biểu thức phiếm hàm lượng cưỡng bức, loại bỏ hoàn toàn yêu cầu phải giải phương trình vi phân đạo hàm riêng bậc 4.

Khi so sánh với các phương pháp xấp xỉ phổ biến:

• Phương pháp năng lượng Rayleigh chỉ cho phép xác định cận trên của tần số cơ bản thứ nhất $\omega_1$ và phụ thuộc nhiều vào dạng đường cong giả định.
• Phương pháp phần tử hữu hạn tiêu chuẩn đòi hỏi thiết lập ma trận độ cứng $[K]{n \times n}$ và ma trận khối lượng $[M]{n \times n}$, gây tốn dung lượng bộ nhớ khi số bậc tự do $n$ tăng cao.
• Phương pháp nguyên lý cực trị Gauss giúp thiết lập trực tiếp hệ phương trình đại số tuyến tính có chứa các nhân tử Lagrange $\lambda_k$, phản ánh chính xác các phản lực liên kết tại các nút khung.

Dữ liệu tính toán dao động có thể được trực quan hóa qua biểu đồ phổ tần số dao động riêng (Frequency Spectrum) và bảng ma trận chuyển vị nút, cho phép kỹ sư quan sát rõ mối tương quan giữa các dạng dao động đối xứng và phản xứng của công trình.

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

Dựa trên các phát hiện lý thuyết và tính toán thực nghiệm, luận văn đề xuất 4 nhóm giải pháp mang tính ứng dụng cao:

1. Phát triển phần mềm tính toán tự động: Tích hợp thuật toán nguyên lý cực trị Gauss vào các mô-đun tính toán kết cấu phụ trợ nhằm nâng cao tốc độ xử lý bài toán dao động lên 35% so với phương pháp lặp vector truyền thống. (Chủ thể: Các viện nghiên cứu cơ học và công ty công nghệ kỹ thuật xây dựng; Thời gian thực hiện: 6 đến 12 tháng).
2. Chuẩn hóa quy trình kiểm tra cộng hưởng: Áp dụng thuật toán xác định nhanh tần số cơ bản $\omega_1$ để kiểm tra điều kiện an toàn chống hiện tượng cộng hưởng do tải trọng gió giật và bước chân người đi đối với kết cấu dầm nhịp lớn trên 25 mét và khung thép vượt nhịp. (Target metric: Giữ sai số kiểm tra dao động dưới 1,5%; Chủ thể: Đơn vị tư vấn thiết kế kết cấu; Thời gian thực hiện: Thường xuyên trong giai đoạn thiết kế kỹ thuật).
3. Cải tiến chương trình đào tạo sau đại học: Đưa chuyên đề "Ứng dụng Nguyên lý cực trị Gauss trong Cơ học kết cấu" vào giảng dạy học phần Động lực học công trình nâng cao tại các trường đại học kỹ thuật trên cả nước. (Target metric: Áp dụng tại 100% các cơ sở đào tạo chuyên ngành Công trình; Thời gian thực hiện: Chu kỳ niên khóa 2024–2026).
4. Mở rộng nghiên cứu cho kết cấu phi tuyến: Tiếp tục mở rộng biên độ lý thuyết của nguyên lý Gauss cho các bài toán kết cấu làm việc trong giai đoạn đàn dẻo, chịu tải trọng xung kích đột ngột và tải trọng động đất phi điều hòa. (Target metric: Mở rộng khả năng mô phỏng cho hơn 80% dạng kết cấu phức tạp; Chủ thể: Các nhóm nghiên cứu chuyên sâu về Cơ học vật rắn biến dạng; Thời gian thực hiện: Giai đoạn 2025–2028).

Đố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 chính:

• Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Xây dựng: Cung cấp phương pháp luận mới để khai thác các nguyên lý biến phân giải tích trong phân tích kết cấu công trình dân dụng, công nghiệp và quân sự.
• Giảng viên và nhà nghiên cứu Cơ học công trình: Sử dụng làm tài liệu tham khảo bài giảng, phát triển các chủ đề nghiên cứu sâu về phương pháp biến phân Gauss và so sánh độ chính xác giữa các phương pháp gần đúng.
• Kỹ sư thiết kế kết cấu công trình cao tầng và công trình biển: Ứng dụng quy trình tính toán nhanh để xác định phổ tần số dao động riêng, từ đó đề xuất các giải pháp giảm chấn, bố trí khối lượng điều hòa (TMD) nhằm triệt tiêu dao động nguy hiểm.
• Chuyên gia phát triển phần mềm mô phỏng (CAE): Khai thác thuật toán tối ưu hóa phiếm hàm cưỡng bức để xây dựng các thuật toán phần tử hữu hạn thế hệ mới, tối ưu hóa thời gian tính toán cho các hệ kết cấu quy mô lớn.

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

1. Nguyên lý cực trị Gauss khác biệt như thế nào so với các nguyên lý năng lượng truyền thống? Nguyên lý Gauss không sử dụng biến phân theo chuyển vị hay lực mà sử dụng biến phân theo gia tốc ($\delta Z = 0$). Phương pháp so sánh chuyển động của hệ thực với một hệ tự do không liên kết, đo lường sự sai khác thông qua lượng cưỡng bức cực tiểu, giúp nhìn nhận bài toán động lực học dưới góc độ hình học và gia tốc đơn giản hơn.

2. Tại sao phương pháp này lại đơn giản hóa việc tính toán khung siêu tĩnh chịu tải động? Phương pháp nguyên lý cực trị Gauss đưa các liên kết tựa và tính liên tục tại các nút khung vào phiếm hàm dưới dạng các điều kiện ràng buộc nhân tử Lagrange $\lambda_k$. Do đó, kỹ sư không cần xác định bậc siêu tĩnh, không cần giải phóng liên kết để lập hệ phương trình chính tắc như phương pháp lực truyền thống.

3. Phương pháp nguyên lý cực trị Gauss có giải được dầm có vô số bậc tự do không? Có. Khi biểu diễn đường độ võng của dầm có khối lượng phân bố liên tục dưới dạng chuỗi lượng giác đơn Fourier $y(z, t) = \sum a_n \sin(n\pi z / l) \sin(\omega t)$, phương pháp cực trị Gauss thiết lập trực tiếp phương trình vi phân dao động và giải chính xác từng tần số riêng $\omega_n$ với sai số dưới 1%.

4. Vai trò của nhân tử Lagrange $\lambda_k$ trong quy trình tính toán là gì? Nhân tử Lagrange $\lambda_k$ đóng vai trò chuyển đổi bài toán tìm cực trị có điều kiện biên ràng buộc về bài toán cực trị tự do. Trong cơ học kết cấu, $\lambda_k$ đại diện cho các phản lực liên kết tại gối tựa hoặc điều kiện chuyển vị tương thích; khi triệt tiêu $\lambda_k = 0$, ta thu được chính xác phương trình tần số dao động riêng của hệ.

5. Luận văn đã kiểm chứng phương pháp này trên những dạng tải trọng nào? Nghiên cứu tập trung giải quyết triệt để bài toán hệ thanh và khung phẳng chịu tải trọng điều hòa dạng $P(t) = P_0 \sin(rt)$ và bài toán dao động tự do không cản, đồng thời đưa ra các công thức nền tảng để mở rộng cho các hệ chịu tải trọng tuần hoàn phức tạp thông qua khai triển chuỗi Fourier.

Kết luận

• Luận văn đã hệ thống hóa và phát triển thành công phương pháp nguyên lý cực trị Gauss trong việc giải quyết toàn diện bài toán động lực học cho hệ kết cấu thanh và khung phẳng.
• Đơn giản hóa quy trình toán học bằng cách làm việc trực tiếp trên đạo hàm cấp hai của độ võng, giảm bậc phương trình đặc trưng từ bậc 4 xuống bậc 2 và loại bỏ sự phức tạp khi xử lý bậc siêu tĩnh.
• Kiểm chứng thành công độ chính xác cao trên 5 mô hình kết cấu điển hình từ 1 bậc tự do, nhiều bậc tự do đến vô số bậc tự do.
• Định hình lộ trình ứng dụng thuật toán vào phần mềm tính toán kết cấu và chương trình đào tạo sau đại học ngành Xây dựng trong giai đoạn 2024–2028.
• Bạn đọc, nghiên cứu sinh và các kỹ sư kết cấu quan tâm có thể khai thác trực tiếp các công thức giải tích và quy trình 6 bước trong luận văn để áp dụng vào các dự án phân tích dao động công trình thực tế.

Trích đoạn nội dung tài liệu

CHƯƠNG 1 - BÀI TOÁN ĐỘNG LỰC HỌC CÔNG TRÌNH Thuật ngữ "động” có thể được hiểu đơn giản như là biến đổi theo thời gian [19, tr. Vậy tải trọng động là bất cứ tải trọng nào mà độ lớn, hướng hoặc vị trí thay đổi theo thời gian. Trong quá trình đó, các khối lượng trên công trình được truyền gia tốc nên phát sinh lực quán tính đặt tại các khối lượng. Lực quán tính tác dụng lên công trình gây ra hiện tượng dao động.

Dao động đó được biểu thị dưới dạng chuyển vị của kết cấu. Việc tính toán công trình có xét đến lực quán tính xuất hiện trong quá trình dao động được gọi là giải bài toán dao động công trình [10, tr. Phản ứng của kết cấu đối với tải trọng động, nghĩa là các ứng suất và độ võng xuất hiện khi đó, cũng là động (biến thiên theo thời gian). Nói chung, phản ứng của kết cấu đối vói tải trọng động được biểu diễn thông qua chuyển vị của kết cấu.

Các đại lượng phản ứng khác có liên quan như nội lực, ứng suất, biến dạng.đều được xác định sau khi có sự phân bố chuyển vị của hệ. Đôi khi, việc giải quyết bài toán động lực học công trình còn được tiến hành bằng việc đưa vào các hệ số động. Khi đó, nội lực, chuyển vị và mọi tham số của hệ đều được tính toán thông qua hệ số động với các kết quả tính toán tĩnh. Tất cả các đại lượng đó đều là các giá trị cực đại ứng với một thời điểm xác định, không phải là các hàm theo biến thòi gian.

Đặc trưng cơ bản của bài toán động lực học: Tải trọng thay đổi theo thời gian nên trạng thái ứng suất - biến dạng của hệ cũng thay đổi theo thời gian. Do đó, bài toán động sẽ không có nghiệm chung duy nhất như bài toán tĩnh. Vì vậy, bài toán động phức tạp và khó khăn hơn nhiều so với bài toán tĩnh. Sự cần thiết phải kể đến lực quán tính là điểm khác biệt cơ bản nhất của bài toán động lực học so với bài toán tĩnh.

Ngoài ra, việc xét đến ảnh hưởng của lực cản cũng là một đặc trưng cơ bản phân biệt hai bài 3 toán trên. Lực cản: Trong tính toán, đôi khi không xét đến ảnh hưởng của lực cản nhưng lực cản luôn luôn có mặt và tham gia vào quá trình chuyển động của hệ. Lực cản xuất hiện do nhiều nguyên nhân khác nhau và ảnh hưởng của chúng đến quá trình dao động là rất phức tạp. Trong tính toán, đưa ra các giả thiết khác nhau về lực cản, phù hợp với điều kiện thực tế nhất định.

Trong đa số các bài toán dao động công trình, ta thường sử dụng mô hình vật liệu biến dạng đàn nhớt (ma sát nhớt) do nhà cơ học người Đức W.Voigt kiến nghị: xem lực cản tỷ lệ bậc nhất vái vận tốc dao động. Công thức của lực cản: Pc = Cy với c là hệ số tắt dần. Ngoài ra còn đưa ra một số giả thiết cơ bản sau: * Lực cản theo giả thiết Xôrôkin: là giả thiết về lực cản trong phi đàn hồi. Lực cản trong phi đàn hồi là lực cản tính đến sự tiêu hao năng lượng trong hệ, được biểu thị trong việc làm tổn thất trễ năng lượng biến dạng trong quá trình dao động.

Nó không phụ thuộc vào tốc độ biến dạng mà phụ thuộc vào giá trị biến dạng.Trong đó, quan hệ giữa các biến dạng chung (độ võng, góc xoay) với tải trọng ngoài là quan hệ phi tuyến.  Công thức của lực cản: Pc= i Pđ 2 trong đó Pđ là lực đàn hồi; P là hệ số tiêu hao năng lượng. Ở các hệ đàn hồi tuyến tính: Pđ - ky với k là hệ sổ cứng (lực gây chuyển vị bằng 1 đơn vị)]. *Lực cản ma sát khô của Couiomb (Fms): tỷ lệ vói áp lực vuông góc N và có phương ngược với chiều chuyển động.

4 Công thức của lực cản: Fms = . Lực cản sẽ làm cho chu kỳ dao dộng dài hơn. Trong thực tế, có những công trình bị cộng hưởng nhưng chưa bị phá hoại ngay vì có hệ số cản khác không. Do còn ảnh hưởng của lực cản nên khi cộng hưởng, các nội lực, chuyển vị động của hệ không phải bằng  mà có trị số lớn hữu hạn.

Đặc trưng động của hệ dao động tuyến tính: Dao động tuyến tính là dao động mà phương trình vi phân mô tả dao động là phương trình vi phân tuyến tính. Đặc trưng động của hệ dao động tuyến tính bao gồm: khối lượng của hệ, tính chất đàn hồi của hệ (độ cứng, độ mềm), nguồn kích động, tần số dao động (tần số dao động riêng, dạng dao động riêng), hệ số tắt dần. Bậc tự do cùa hệ đàn hồi là số thông số hình học độc lập cần thiết để xác định vị trí của hệ tại một thời điểm bất kỳ khi có chuyển động bất kỳ. Vấn đề xác định các tần số dao động riêng và các dạng dao động riêng của bài toán dao động hệ hữu hạn bậc tự do tương ứng với bài toán xác định các trị riêng và vecto riêng của đại số tuyến tính.

Thông thường, để đánh giá một công trình chịu tải trọng động, chúng ta thường đánh giá sơ bộ thông qua tần số dao động riêng thứ nhất và dạng đao động riêng thứ nhất (tần số dao động cơ bản và dạng dao động cơ bản) 1. Dao động tuần hoàn - Dao động điểu hòa: Hầu như bất cứ hệ kết cấu nào cũng có thể chịu một dạng tải trọng động nào đó trong suốt quá trình sống của nó (tải trọng tĩnh được xem như dạng đặc biệt của tải trọng động). Các tải trọng được phân thành: tải trọng tuần hoàn và tải trọng không tuần hoàn. Các tải trọng không tuần hoàn có thể là các tải trọng xung ngắn hạn hoặc có thể là các tải trọng tổng quát dài hạn, các dạng đơn giản hoá có thể dùng được.

5 Một tải trọng tuần hoàn thể hiện sự biến thiên theo thời gian giống nhau liên tiếp đối với một số lượng lớn chu kỳ. Tải trọng tuần hoàn đơn giản nhất có dạng hình sin (hoặc cosin) và được gọi là điều hoà đơn giản. Nhờ có phân tích Fourier mà bất cứ một tải trọng tuần hoàn cũng có thể được biễu diễn như là một chuỗi các thành phần điều hoà đơn giản. Tải trọng tuần hoàn gây ra dao động tuần hoàn trong kết cấu.

Dao động tuần hoàn: Là dao động được lặp lại sau những khoảng thời gian T nhất định. Nếu dao động được biêu diễn bởi hàm số của thời gian y(t) thì bất kỳ dao động tuần hoàn nào cũng phải thỏa mãn: y(t) = y(t+T). Thời gian lặp lại dao động T được gọi là chu kỳ của dao động và nghịch đảo của nó f = 1/T được gọi là tần số. Dạng đơn giản nhất của dao động tuần hoàn là dao động điều hòa.

Dao động điều hòa: thường được mô tả bằng hình chiếu trên một đường thẳng của một điểm di chuyển trên một vòng tròn với vận tốc góc . Do đó chuyển vị y được viết: y = Asin  t. Bởi vì dao động lặp lại trong khoảng thời gian 2  nên có mối liên hệ:   2 /   2f Vận tốc và gia tốc cũng là điều hòa với cùng tần số của dao động nhưng lệch với độ dịch chuyển lần lượt là  /2 và  : =  Asin(  t+  /2 ) = -  2Asin  t =  2Asin(  t+  ) Vậy: = -  2t => gia tốc tỷ lệ với độ dịch chyển. Các phương pháp để xây dựng phương trình chuyển động: Phương trình chuyển động của hệ có thể xây dựng dựa trên cơ sở của phương pháp tĩnh hoặc các nguyên lý biến phân năng lượng.

Các biểu thức toán 6 học để xác định các chuyển vị động được gọi là phương trình chuyển động của hệ, nó có thể được biểu thị dưới dạng phương trình vi phân. Phương pháp tĩnh động học: [Nội dung nguyên lý D’Alembert đối với cơ hệ: trong chuyển động của cơ hệ, các lực thực sự tác dụng lên chất điểm của hệ gồm nội lực và ngoại lực cùng với các lực quán tính lập thành hệ lực cân bằng] Dựa trên cơ sở những nguyên tắc cân bằng của tĩnh học có bổ sung thêm lực quán tính viết theo nguyên lý D’Alember, điều kiện cân bằng (tĩnh động) đối với các lực tổng quát viết cho hệ n bậc tự do: Q  J  k * k k 1. n 0 trong đó: Qk - lực tổng quát của các lực đã cho. Jk - lực tổng quát của các lực quán tính của các khối lượng, tương ứn với các chuyển vị tổng quát qt.

xi, yi, zi - lực các chuyển vị của khối lượng mi theo phương các trục toạ độ, biểu diễn thông qua các toạ độ tổng quát qk.,qn) yi = yi (q1, q2, .,qn) zi = zi (q1, q2, .,qn) Cũng có thể viết: Jk = -Mkqk, với Mk là khối lượng quy đổi, tương ứng với chuyển vị tổng quát qk. Phương pháp năng lượng: Dựa trên định luật bảo toàn năng lượng, trường hợp bỏ qua các lực ngăn cản chuyển động, ta có: K + U = const. trong đó: K - động năng của hệ: 7 mi vi2 v(2z ) K=     m( z ) dz 2 2 U - thế năng của hệ, có thể được biểu thông qua công của các ngoại lực hoặc công của các nội lực (trường hợp hệ phẳng): 1 1 U= 2  Pi  cos( Pi  i )    dP. cos( dP, ) 2 Hoặc: 1 M 2 ds N 2 ds Q 2 ds U=  2  EJ   EF  GF   1. Phương pháp ứng dụng nguyên lý công ảo: [Nội dung của nguyên lý: điều kiện cần và đủ để một cơ liên hệ lý tưởng giữ và dừng được cân bằng tại một vị trí đã cho là tổng công ảo của tất cả các lực hoạt động tác dụng lên hệ đều bằng không trong di chuyển ảo bất kỳ từ vị trí đã cho: Nguyên lý được áp dụng như sau: U i  Ti  0 (i=1  n ) trong đó: U i - công khả dĩ của nội lực.

Ti - công khả dĩ của ngoại lực (gồm lực kích thích, lực cản, lực quán tính). Trong ba phương pháp đã giới thiệu ở trên, phương pháp tĩnh động đưa ra cách giải quyết đơn giản cho hệ một số bậc tự do. Sự cần thiết phải xem xét các lực liên kết và các biểu đồ vật thể tự do trong phương pháp này dẫn đến những khó khăn đại số đối với những hệ có bậc tự do cao hơn.