Phân tích động (tiếng Anh: dynamic analysis; hay còn gọi là phân tích động lực học) là phương pháp khảo sát, mô hình hóa và đánh giá hành vi biến thiên theo thời gian của một hệ thống vật lý, kết cấu công trình hoặc chương trình phần mềm khi chịu tác động của các tải trọng, ngoại lực hoặc dữ liệu đầu vào phụ thuộc thời gian. Trong cơ học kết cấu và kỹ thuật công trình, phân tích động đóng vai trò sống còn trong việc kiểm tra độ an toàn của nhà cao tầng, cầu dây văng, đập thủy điện và công trình biển chịu tải trọng gió bão, động đất, sóng biển và rung động máy móc. Trong khoa học máy tính, phân tích động là phương pháp kiểm thử thực thi nhằm phát hiện lỗi bộ nhớ và lỗ hổng bảo mật. Bài viết này tập trung làm rõ cơ sở toán lý của phân tích động lực học kết cấu, các phương pháp số giải bài toán động và các ứng dụng kỹ thuật thực tiễn.
Cơ sở lý thuyết của phân tích động lực học kết cấu
Phương trình vi phân tổng quát mô tả chuyển động của một hệ kết cấu rời rạc nhiều bậc tự do (MDOF) được thiết lập dựa trên định luật II Newton hoặc nguyên lý D'Alembert:
\[\mathbf{M} \mathbf{\ddot{u}}(t) + \mathbf{C} \mathbf{\dot{u}}(t) + \mathbf{K} \mathbf{u}(t) = \mathbf{F}(t)\]
Trong đó:
- \(\mathbf{M}\) là ma trận khối lượng của hệ kết cấu (đại diện cho lực quán tính \(\mathbf{F}_I = \mathbf{M} \mathbf{\ddot{u}}\)).
- \(\mathbf{C}\) là ma trận cản nhớt (đại diện cho sự tiêu tán năng lượng dao động \(\mathbf{F}_D = \mathbf{C} \mathbf{\dot{u}}\), thường được mô hình hóa theo cản Rayleigh: \(\mathbf{C} = lpha \mathbf{M} + eta \mathbf{K}\)).
- \(\mathbf{K}\) là ma trận độ cứng của kết cấu (đại diện cho nội lực đàn hồi \(\mathbf{F}_S = \mathbf{K} \mathbf{u}\)).
- \(\mathbf{u}(t), \mathbf{\dot{u}}(t), \mathbf{\ddot{u}}(t)\) lần lượt là các vectơ chuyển vị, vận tốc và gia tốc của các nút kết cấu tại thời điểm \(t\).
- \(\mathbf{F}(t)\) là vectơ ngoại lực kích thích biến thiên theo thời gian.
So sánh phân tích tĩnh và phân tích động
Bảng đối chiếu dưới đây làm rõ các điểm khác biệt cơ bản giữa hai phương pháp tính toán:
| Tiêu chí so sánh | Phân tích tĩnh (Static Analysis) | Phân tích động (Dynamic Analysis) |
|---|---|---|
| Đặc tính tải trọng | Tải trọng không đổi hoặc thay đổi rất chậm theo thời gian | Tải trọng thay đổi nhanh, có tần số kích thích rõ rệt hoặc xung kích |
| Hiệu ứng quán tính và cản | Hoàn toàn bỏ qua (\(\mathbf{M}=\mathbf{0}, \mathbf{C}=\mathbf{0}\)) | Được tính toán đầy đủ và quyết định biên độ dao động |
| Dạng phương trình | Hệ phương trình đại số tuyến tính \(\mathbf{K}\mathbf{u} = \mathbf{F}\) | Hệ phương trình vi phân bậc hai phụ thuộc thời gian |
| Hiện tượng cộng hưởng | Không xuất hiện | Xảy ra khi tần số ngoại lực trùng với tần số dao động riêng (\(\omega pprox \omega_n\)) |
| Chi phí tính toán | Thấp, giải nhanh bằng các thuật toán ma trận thông thường | Cao hơn nhiều lần, đòi hỏi tích phân từng bước thời gian |
Các phương pháp phân tích động phổ biến
Tùy thuộc vào bản chất của tải trọng và mục tiêu thiết kế, các kỹ sư sử dụng các công cụ phân tích động chuyên biệt:
1. Phân tích dao động riêng (Modal Analysis)
Phân tích dao động tự do không cản (\(\mathbf{M}\mathbf{\ddot{u}} + \mathbf{K}\mathbf{u} = \mathbf{0}\)) nhằm tìm các giá trị riêng \(\omega_i^2\) (tần số góc riêng) và vectơ riêng \(oldsymbol{\phi}_i\) (dạng dao động riêng / mode shape). Đây là bước tiên quyết để hiểu rõ đặc trưng cơ học của công trình và xác định xem công trình có nguy cơ cộng hưởng với gió hoặc động đất hay không.
2. Phân tích phổ phản ứng (Response Spectrum Analysis - RSA)
Phương pháp này kết hợp các phản ứng cực đại của từng mode dao động (sử dụng các quy tắc cộng như SRSS hoặc CQC) dựa trên một phổ phản ứng thiết kế chuẩn. Đây là phương pháp chủ đạo trong tiêu chuẩn thiết kế kháng chấn công trình tại Việt Nam (TCVN 9386:2012) nhờ sự cân bằng tuyệt vời giữa độ chính xác và thời gian tính toán.
3. Phân tích lịch sử thời gian (Time History Analysis - THA)
Phương pháp mô phỏng chính xác nhất hành vi kết cấu bằng cách tích phân trực tiếp hệ phương trình vi phân qua từng bước thời gian \(\Delta t\) (sử dụng các thuật toán như Newmark-\(eta\), Wilson-\( heta\), hoặc Runge-Kutta). THA cho phép mô phỏng chi tiết phản ứng phi tuyến vật liệu, sự hình thành khớp dẻo và hư hại cục bộ của công trình dưới tác động của các giản đồ gia tốc động đất thực tế đã ghi nhận trong lịch sử.
Ứng dụng thực tiễn trong công trình và công nghệ
- Thiết kế nhà siêu cao tầng: Kiểm tra dao động ngang do tải trọng gió giật (wind gust) và kiểm soát gia tốc dao động đỉnh tháp nhằm đảm bảo độ êm dịu sinh lý cho người sử dụng; thiết kế hệ thống giảm chấn khối lượng điều chỉnh (Tuned Mass Damper - TMD).
- Cầu nhịp lớn và cầu dây văng: Phân tích hiện tượng khí động đàn hồi (flutter, vortex shedding) gây mất ổn định khí động học tương tự thảm họa sập cầu Tacoma Narrows năm 1940.
- Công trình công nghiệp và nền móng máy: Tính toán chống rung cho bệ móng máy phát điện tuabin, máy ép thủy lực và máy nén khí công suất lớn.
- Kỹ thuật phần mềm và bảo mật: Giám sát luồng thực thi chương trình trong môi trường sandbox để phát hiện mã độc (malware profiling) và lỗi tràn bộ đệm (buffer overflow).