Đánh Giá Tần Số Dao Động Của Khung Thép Phẳng Có Liên Kết Nửa Cứng Ngẫu Nhiên Bằng MATLAB & Mô Phỏng Monte Carlo

Tóm tắt nghiên cứu (Executive Summary)

  • Câu hỏi nghiên cứu cốt lõi: Tính bất định và ngẫu nhiên trong độ cứng của các liên kết dầm – cột ảnh hưởng như thế nào đến các tần số dao động riêng của hệ kết cấu khung thép phẳng nhiều tầng?
  • Phương pháp tiếp cận: Đề tài kết hợp phương pháp Phần tử hữu hạn (PTHH) với kỹ thuật mô phỏng thống kê ngẫu nhiên Monte Carlo ($10.000$ mẫu lặp). Mô hình toán học được thiết lập dựa trên việc tích hợp các phần tử dầm có liên kết lò xo xoay đàn hồi ở hai đầu vào ma trận độ cứng $[K]$ và ma trận khối lượng $[M]$ tổng thể, được lập trình tự động hóa hoàn toàn trên môi trường MATLAB.
  • Phát hiện chính:
    • Xác lập thành công công thức giải tích ma trận độ cứng và ma trận khối lượng cho phần tử dầm có liên kết nửa cứng với các hệ số biến thiên độ cứng $r_1, r_2 \in [0, 1]$.
    • Phân bố xác suất của các tần số dao động riêng đầu tiên ($\omega_1, \omega_2, \omega_3$) tuân theo quy luật phân bố chuẩn (Gaussian).
    • Độ lệch chuẩn (SD) và hệ số biến thiên (COV) của tần số dao động riêng tỷ lệ thuận tuyến tính với độ lệch chuẩn của độ cứng liên kết đầu vào.
    • Tần số dao động bậc thấp (tần số cơ bản) có độ ổn định thống kê cao hơn, ít bị ảnh hưởng bởi tính ngẫu nhiên của liên kết so với các bậc dao động cao.
  • Ý nghĩa thực tiễn: Nghiên cứu đặt nền móng chuyển dịch từ phương pháp thiết kế tiền định truyền thống sang đánh giá độ tin cậy động lực học xác suất, giúp kỹ sư công trình kiểm soát chính xác vùng tần số nguy hiểm nhằm chống cộng hưởng do gió động và động đất.

Bối cảnh và tầm quan trọng (Context & Significance)

Trong xu thế phát triển của ngành xây dựng hiện đại, kết cấu khung thép ngày càng chiếm ưu thế vượt trội trong các công trình dân dụng, công nghiệp và cầu đường nhờ ưu điểm vượt nhịp lớn, tỷ trọng nhẹ và tốc độ lắp dựng nhanh chóng. Khi thiết kế khung thép, các cấu kiện dầm và cột được gia công rời tại nhà máy và liên kết tại công trường thông qua bu lông cường độ cao, mối hàn hoặc thép góc.

+-------------------------------------------------------------------------+
|                  QUAN NIỆM THIẾT KẾ LIÊN KẾT DẦM - CỘT                  |
+-------------------------------------------------------------------------+
|  Mô hình truyền thống (Tiền định):                                     |
|  [ Ngàm cứng tuyệt đối (Rigid) ]  <-------- Hoặc -------->  [ Khớp lý tưởng (Pinned) ] |
|                                                                         |
|  Thực tế làm việc công trình:                                           |
|  =======> [ LIÊN KẾT NỬA CỨNG (Semi-Rigid) ]                            |
|           + Tồn tại độ xoay tương đối giữa dầm và cột                   |
|           + Ứng xử phi tuyến mô men - góc xoay (M - Theta)               |
|           + Chứa đựng tính bất định / ngẫu nhiên do thi công & vật liệu |
+-------------------------------------------------------------------------+

Khoảng trống nghiên cứu (Research Gap)

Theo phương pháp thiết kế truyền thống, các kỹ sư thường lý tưởng hóa nút liên kết theo hai thái cực: ngàm cứng tuyệt đối (không có chuyển vị xoay tương đối) hoặc khớp lý tưởng (khả năng truyền mô men bằng $0$). Tuy nhiên, vô số thực nghiệm của các nhà khoa học trên thế giới (như Nethercot, Kishi & Chen, Bjorhovde) đã chứng minh rằng hầu hết các liên kết thực tế đều là liên kết nửa cứng (semi-rigid connections).

Hơn thế nữa, trong quá trình gia công, lắp dựng bu lông, biến động chất lượng thép và ứng suất dư, độ cứng thực tế của từng nút khung luôn có dung sai và mang tính bất định ngẫu nhiên. Phần lớn các nghiên cứu trong nước trước đây chỉ khảo sát bài toán tĩnh hoặc bài toán động lực học tiền định (deterministic) với độ cứng liên kết là một giá trị cố định, chưa phản ánh bức tranh phân tán xác suất của phản ứng dao động thực tế.

Tính cấp thiết và tác động

Việc xác định chính xác tần số dao động riêng là yêu cầu sống còn trong phân tích động lực học kết cấu thép. Nếu tần số dao động riêng của công trình trùng hoặc tiệm cận với tần số của tải trọng kích thích (gió giật, xung lực động đất), hiện tượng cộng hưởng (resonance) sẽ xảy ra, dẫn đến phá hoại kết cấu. Đề tài giải quyết trực tiếp bài toán tính toán động lực học ngẫu nhiên, cung cấp cơ sở khoa học để tối ưu hóa vật liệu và nâng cao độ an toàn cho công trình theo các tiêu chuẩn tiên tiến như Eurocode 3 và AISC.

Phương pháp luận và cách tiếp cận (Methodology & Approach)

Nghiên cứu ứng dụng quy trình phân tích kết hợp giữa phương pháp Phần tử hữu hạn giải tích và kỹ thuật mô phỏng thống kê ngẫu nhiên Monte Carlo.

+-------------------------------------------------------------------------+
|                    QUY TRÌNH MÔ PHỎNG VÀ TÍNH TOÁN                      |
+-------------------------------------------------------------------------+
|  1. THIẾT LẬP MÔ HÌNH PTHH                                              |
|     - Cột: Phần tử dầm-cột ngàm cứng hai đầu                            |
|     - Dầm: Phần tử thanh tích hợp lò xo xoay bán cứng ở 2 đầu           |
|     - Xây dựng ma trận độ cứng [K] và ma trận khối lượng [M] giải tích   |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
|  2. MÔ HÌNH HÓA THAM SỐ NGẪU NHIÊN                                      |
|     - Biến ngẫu nhiên: Độ cứng liên kết Rk                              |
|     - Thông số chuẩn hóa: r_i = 1 / (1 + 3EI / (L*Rki)) in [0, 1]       |
|     - Phân bố giả định: Chuẩn N(mu = 0.8, sigma = 0.25)                 |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
|  3. MÔ PHỎNG MONTE CARLO TRÊN MATLAB (N = 10.000 mẫu)                   |
|     - Gieo biến ngẫu nhiên cho toàn bộ 80 nút khung dầm-cột             |
|     - Ghép nối ma trận tổng thể [K] & [M]                               |
|     - Giải bài toán trị riêng: det([K] - omega^2 * [M]) = 0             |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
|  4. XỬ LÝ THỐNG KÊ KẾT QUẢ ĐẦU RA                                       |
|     - Giá trị trung bình (Mean - M)                                     |
|     - Độ lệch chuẩn (Standard Deviation - SD)                           |
|     - Hệ số biến thiên (Coefficient of Variation - COV)                 |
|     - Vẽ đồ thị hàm mật độ phân bố tần suất (Histogram)                |
+-------------------------------------------------------------------------+

1. Thiết lập ma trận độ cứng và ma trận khối lượng phần tử dầm nửa cứng

Phần tử dầm có chiều dài $L$, độ cứng uốn $EI$, mật độ khối lượng $m$, được gắn hai lò xo chống xoay ở hai đầu với độ cứng $R_{kA}$ và $R_{kB}$. Mức độ cứng của liên kết được chuẩn hóa qua đại lượng không thứ nguyên:

$$r_i = \frac{1}{1 + \frac{3EI}{L \cdot R_{ki}}} \quad (i = A, B)$$

  • Khi $r_i = 1 \rightarrow R_{ki} \rightarrow \infty$: Nút ngàm cứng hoàn toàn.
  • Khi $r_i = 0 \rightarrow R_{ki} = 0$: Nút khớp lý tưởng.
  • Khi $0 < r_i < 1$: Liên kết nửa cứng.

Dựa trên nguyên lý công khả dĩ (Principle of Virtual Displacements), nghiên cứu thiết lập ma trận độ cứng cục bộ $[K_e]$ cấp $6 \times 6$ có dạng:

$$[K_e] = \begin{bmatrix} \frac{EA}{L} & 0 & 0 & -\frac{EA}{L} & 0 & 0 \ 0 & \frac{12EI(r_1+r_2+r_1r_2)}{L^3(4-r_1r_2)} & \frac{6EI(2r_1+r_1r_2)}{L^2(4-r_1r_2)} & 0 & -\frac{12EI(r_1+r_2+r_1r_2)}{L^3(4-r_1r_2)} & \frac{6EI(2r_2+r_1r_2)}{L^2(4-r_1r_2)} \ 0 & \frac{6EI(2r_1+r_1r_2)}{L^2(4-r_1r_2)} & \frac{4EI(3r_1)}{L(4-r_1r_2)} & 0 & -\frac{6EI(2r_1+r_1r_2)}{L^2(4-r_1r_2)} & \frac{2EI(3r_1r_2)}{L(4-r_1r_2)} \ -\frac{EA}{L} & 0 & 0 & \frac{EA}{L} & 0 & 0 \ 0 & -\frac{12EI(r_1+r_2+r_1r_2)}{L^3(4-r_1r_2)} & -\frac{6EI(2r_1+r_1r_2)}{L^2(4-r_1r_2)} & 0 & \frac{12EI(r_1+r_2+r_1r_2)}{L^3(4-r_1r_2)} & -\frac{6EI(2r_2+r_1r_2)}{L^2(4-r_1r_2)} \ 0 & \frac{6EI(2r_2+r_1r_2)}{L^2(4-r_1r_2)} & \frac{2EI(3r_1r_2)}{L(4-r_1r_2)} & 0 & -\frac{6EI(2r_2+r_1r_2)}{L^2(4-r_1r_2)} & \frac{4EI(3r_2)}{L(4-r_1r_2)} \end{bmatrix}$$

Đồng thời, ma trận khối lượng phân bố nhất quán $[M_e]$ được tính toán chính xác có xét đến ảnh hưởng của $r_1, r_2$ thông qua các đa thức hàm dạng bậc cao.

2. Mô phỏng số Monte Carlo trên khung không gian 20 tầng

  • Mô hình kết cấu: Khung thép phẳng 20 tầng, 2 khoang (chiều cao tầng $H = 3.65\text{ m}$, nhịp khoang $L = 9.0\text{ m}$).
  • Quy mô phần tử: Gồm 60 phần tử cột và 40 phần tử dầm, tương ứng với 80 nút liên kết nửa cứng dầm – cột độc lập.
  • Vật liệu: Thép kết cấu ASTM A-36 ($E = 2.0 \times 10^8\text{ kN/m}^2$, trọng lượng riêng $\gamma = 78.5\text{ kN/m}^3$).
  • Số lần gieo mẫu: $N = 10.000$ lần trên thuật toán EQ_Building.m viết bằng ngôn ngữ MATLAB.

3. Kiểm chứng độ tin cậy (Validation)

Để đảm bảo tính đúng đắn của thuật toán, chương trình tính toán được kiểm chứng với trường hợp nút ngàm cứng ($r_1 = r_2 = 1$). Kết quả tần số dao động riêng bậc nhất thu được ($\omega_1 = 4.102\text{ rad/s}$) hoàn toàn trùng khớp với kết quả công bố của tác giả Filho và mô hình giải trên phần mềm thương mại SAP2000.

Phát hiện chính của nghiên cứu (Key Findings)

Từ kết quả phân tích thống kê $10.000$ kịch bản mô phỏng ngẫu nhiên, nghiên cứu rút ra 4 phát hiện quan trọng:

+-----------------------------------------------------------------------------+
|               BẢNG TỔNG HỢP KẾT QUẢ ĐẶC TRƯNG THỐNG KÊ TẦN SỐ               |
+-------------------+--------------------+-------------------+----------------+
| Dạng dao động     | Kỳ vọng trung bình | Độ lệch chuẩn     | Hệ số biến     |
| (Mode)            | M (rad/s)          | SD (rad/s)        | thiên COV (%)  |
+-------------------+--------------------+-------------------+----------------+
| Mode 1 (omega_1)  | 3.842              | 0.048             | 1.25%          |
| Mode 2 (omega_2)  | 11.650             | 0.185             | 1.59%          |
| Mode 3 (omega_3)  | 20.315             | 0.421             | 2.07%          |
+-------------------+--------------------+-------------------+----------------+
*Ghi chú: Giá trị thu được ứng với độ cứng liên kết trung bình r_0 = 0.8, độ lệch chuẩn liên kết sigma = 0.25.

1. Phân bố xác suất tần số riêng mang tính chuẩn hóa (Gaussian Distribution)

Biểu đồ tần suất (Histogram) của 3 tần số dao động riêng đầu tiên sau $10.000$ lần mô phỏng đều có dạng hình chuông đối xứng hoàn hảo, đặc trưng cho phân bố chuẩn (Normal Distribution). Điều này chứng minh rằng dù kết cấu có số lượng nút liên kết lớn và phân tán ngẫu nhiên, phản ứng dao động tổng thể của hệ vẫn tuân theo quy luật giới hạn trung tâm.

2. Quan hệ tuyến tính giữa độ bất định liên kết và độ biến thiên tần số

Nghiên cứu chỉ ra rằng khi độ lệch chuẩn của độ cứng liên kết đầu vào ($\sigma_r$) tăng dần, độ lệch chuẩn (SD) và hệ số biến thiên (COV) của các tần số dao động riêng cũng tăng theo một mối quan hệ tuyến tính đồng biến. Độ phân tán của chất lượng thi công nút khung phản ánh trực tiếp và tương đương lên độ phân tán của tần số dao động riêng của toàn bộ công trình.

3. Quy luật phân tầng độ nhạy dao động theo bậc

  • Tần số dao động cơ bản ($\omega_1$) có hệ số biến thiên thấp nhất ($\text{COV} \approx 1.25%$).
  • Tần số dao động bậc cao ($\omega_2, \omega_3$) có hệ số biến thiên tăng dần ($\text{COV}$ đạt $1.59%$ và $2.07%$).

Giải thích cơ học: Tần số cơ bản ($\omega_1$) chi phối chuyển vị uốn tổng thể của toàn bộ khung nhà, nơi hiệu ứng san bằng ngẫu nhiên của 80 nút liên kết phát huy tác dụng. Ngược lại, các tần số bậc cao gắn liền với các dạng dao động uốn cục bộ nhiều tầng, do đó nhạy cảm hơn rõ rệt với sự thay đổi độ cứng cục bộ của từng nút dầm – cột.

4. Hiệu ứng làm mềm kết cấu thực tế so với thiết kế lý tưởng

Khi độ cứng liên kết dầm – cột có giá trị trung bình $r_0 = 0.8$, tần số dao động cơ bản thực tế của công trình giảm từ $4.102\text{ rad/s}$ (mô hình ngàm cứng) xuống còn $3.842\text{ rad/s}$ (giảm khoảng $6.3%$). Điều này đồng nghĩa với việc chu kỳ dao động tự nhiên $T$ của công trình thực tế bị kéo dài ra, làm thay đổi phổ phản ứng gia tốc thiết kế kháng chấn.

Đóng góp khoa học và giá trị ứng dụng (Scientific Contributions)

Đóng góp về mặt học thuật và lý thuyết

  • Hoàn thiện mô hình giải tích: Thiết lập hệ ma trận phần tử giải tích nhất quán ($[K]$ và $[M]$) cho phần tử dầm đàn hồi có liên kết nửa cứng hai đầu, bổ sung công cụ phân tích giải tích tường minh cho lý thuyết cơ học kết cấu.
  • Tiên phong tích hợp xác suất: Đưa lý thuyết xác suất và kỹ thuật mô phỏng Monte Carlo vào bài toán phân tích động lực học kết cấu thép trong nước, vượt qua rào cản của các phương pháp tiền định truyền thống.

Đổi mới về mặt phương pháp luận

  • Xây dựng bộ mã nguồn mở hoàn chỉnh trên nền tảng MATLAB (EQ_Building.m, topology.m, beam2d2.m, eigen.m), cho phép tự động hóa quy trình phân tích trị riêng và xử lý dữ liệu lớn ($10.000$ mẫu) với hiệu năng tính toán cao.

Ứng dụng thực tiễn trong công nghiệp xây dựng

  • Tối ưu hóa thiết kế kháng chấn: Cho phép các kỹ sư tính toán chính xác dải biến thiên của tần số dao động riêng thay vì một giá trị đơn lẻ, từ đó loại trừ hoàn toàn nguy cơ cộng hưởng với chu kỳ dao động của nền đất.
  • Tiết kiệm chi phí đầu tư: Nhận diện đúng ứng xử nửa cứng giúp tận dụng khả năng tiêu tán năng lượng của nút khung, giảm kích thước tiết diện cấu kiện dầm cột mà vẫn đảm bảo độ tin cậy kết cấu.

Đối tượng thụ hưởng và quan tâm (Target Audience)

+-------------------------------------------------------------------------+
|                         ĐỐI TƯỢNG QUAN TÂM                              |
+-------------------------------------------------------------------------+
|  [ Kỹ sư Thiết kế & Thẩm tra Kết cấu ]                                  |
|  --> Đánh giá chính xác chu kỳ dao động thực của khung nhà cao tầng     |
|  --> Tối ưu hóa cấu tạo liên kết bu lông dầm - cột                      |
+-------------------------------------------------------------------------+
|  [ Nhà nghiên cứu & Giảng viên Cơ học Công trình ]                      |
|  --> Tham khảo mô hình giải tích ma trận khối lượng / độ cứng mở rộng   |
|  --> Ứng dụng thuật toán Monte Carlo trong phân tích độ tin cậy         |
+-------------------------------------------------------------------------+
|  [ Sinh viên & Học viên Sau đại học chuyên ngành Công trình ]           |
|  --> Khai thác mã nguồn MATLAB mẫu cho đồ án và đề tài NCKH             |
|  --> Nắm bắt phương pháp nghiên cứu động lực học ngẫu nhiên             |
+-------------------------------------------------------------------------+

Các câu hỏi thường gặp (FAQ)

1. Phát hiện quan trọng nhất của đề tài này là gì?

Phát hiện cốt lõi là tần số dao động riêng của khung thép thực tế tuân theo quy luật phân bố chuẩn (Gaussian) và độ lệch chuẩn của tần số tỷ lệ thuận tuyến tính với độ phân tán của độ cứng nút khung. Tần số dao động thực tế luôn nhỏ hơn so với giả định ngàm cứng lý thuyết.

2. Điểm đặc biệt trong phương pháp nghiên cứu của tác giả là gì?

Đề tài đã tích hợp trực tiếp tham số nửa cứng $r_1, r_2$ vào hàm dạng để thiết lập ma trận độ cứng và ma trận khối lượng giải tích nhất quán, sau đó kết hợp thuật toán Monte Carlo quy mô $10.000$ mẫu trên MATLAB để phân tích độ nhạy ngẫu nhiên.

3. Kết quả nghiên cứu có thể tổng quát hóa cho các dạng công trình khác không?

Có. Thuật toán và mô hình phần tử hữu hạn này có thể mở rộng áp dụng cho mọi kết cấu khung thép phẳng nhiều tầng, nhiều nhịp, nhà xưởng công nghiệp khẩu độ lớn và kết cấu cầu dàn thép có liên kết bu lông/đinh tán.

4. Những bước phát triển tiếp theo của hướng nghiên cứu này là gì?

Các hướng nghiên cứu tiếp theo bao gồm:

  • Mở rộng từ khung phẳng 2D sang mô hình khung không gian 3D chịu xoắn.
  • Khảo sát ứng xử phi tuyến đàn dẻo và suy giảm độ cứng theo chu kỳ (Cyclic Hysteresis Degradation).
  • Phân tích phản ứng động học kết cấu dưới tác động của phổ gia tốc động đất ngẫu nhiên theo miền thời gian.

5. Kỹ sư xây dựng có thể ứng dụng kết quả này vào thiết kế thực tế như thế nào?

Kỹ sư có thể dùng công thức ma trận mở rộng để điều chỉnh chu kỳ dao động riêng trong mô hình tính toán, đánh giá khoảng dao động an toàn của kết cấu và thiết kế liên kết bu lông linh hoạt nhằm giảm giá thành thi công mà không làm suy giảm độ an toàn công trình.

Kết luận (Conclusion)

Báo cáo nghiên cứu khoa học của tác giả Võ Phú Toàn (GVHD: ThS. Nguyễn Duy Hưng, Trường Đại học Giao thông Vận tải Phân hiệu tại TP. Hồ Chí Minh) đã giải quyết thành công bài toán động lực học ngẫu nhiên của khung thép phẳng có liên kết nửa cứng. Bằng việc làm sáng tỏ quy luật phân bố xác suất và mối tương quan tuyến tính giữa độ bất định liên kết và tần số riêng thông qua công cụ MATLAB và mô phỏng Monte Carlo, nghiên cứu đã tạo ra bước tiến quan trọng trong việc ứng dụng cơ học xác suất vào thực tiễn kỹ thuật công trình.

Các đơn vị tư vấn thiết kế và các nhà nghiên cứu được khuyến khích tiếp cận các mô hình liên kết nửa cứng ngẫu nhiên này nhằm nâng cao chất lượng thẩm tra, tối ưu hóa vật liệu và đảm bảo an toàn tuyệt đối cho các công trình kết cấu thép trước các tác động động lực học phức tạp.