Giới thiệu dự án

Trong lĩnh vực kỹ thuật công trình giao thông, việc đánh giá và thiết kế công trình cầu chịu tải trọng động đất là bài toán địa kỹ thuật và động lực học kết cấu tối quan trọng. Theo thống kê địa vật lý toàn cầu, hơn 95% các trận động đất phát sinh do sự vận động kiến tạo của thạch quyển, gây ra các đợt giải phóng năng lượng khổng lồ lan truyền dưới dạng sóng dọc (sóng $P$ với vận tốc $V_p = 7 - 8\text{ km/s}$), sóng ngang (sóng $S$ với vận tốc $V_s = 4 - 5\text{ km/s}$) và sóng mặt Rayleigh/Love ($V = 0.5 - 1.5\text{ km/s}$). Tại Việt Nam, các nghiên cứu địa chấn đã ghi nhận trên 700 trận động đất có độ lớn $M_s \le 7.0$, đặc biệt là các trận động đất mạnh tại Điện Biên (1935, $M_s = 6.8$, cấp chấn động $I_0 = 8 - 9$) và Tuần Giáo (1983, $M_s = 6.7$, $I_0 = 8 - 9$), tập trung chủ yếu dọc theo các đới đứt gãy kiến tạo sâu như Sông Mã, Sông Đà, Sông Chảy và Đông Triều.

       +-----------------------------------------------------------+
       |           Nguồn kích động địa chấn (Chấn tiêu)            |
       |      - Giải phóng năng lượng biến dạng đàn hồi            |
       |      - Sóng P (Vp = 7-8 km/s), Sóng S (Vs = 4-5 km/s)     |
       +-----------------------------+-----------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       |        Lan truyền sóng & Phản ứng nền đất (PGA)           |
       |      - Vận tốc sóng mặt: 0.5 - 1.5 km/s                   |
       |      - Hiện tượng hóa lỏng & Mất ổn định nền móng         |
       +-----------------------------+-----------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       |           Ứng xử động lực học công trình cầu              |
       |      - Dao động đàn hồi đa bậc tự do (MDOF)               |
       |      - Lực quán tính D'Alembert & Lực cắt đáy chân trụ    |
       +-----------------------------------------------------------+

Vấn đề nghiên cứu (Problem Statement)

Thực tiễn thiết kế công trình cầu tại Việt Nam từng tồn tại khoảng trống công nghệ và tiêu chuẩn kỹ thuật sâu sắc:

  • Tiêu chuẩn ngành cũ 22 TCN 221-95 chủ yếu dựa trên phương pháp tĩnh lực giả định xem kết cấu như vật thể tuyệt đối cứng ($F = K_s \cdot Q$), bỏ qua chu kỳ dao động riêng và đặc tính đàn hồi động lực học.
  • Khi hội nhập quốc tế, các dự án vốn vay ODA áp dụng tiêu chuẩn AASHTO 1996, AASHTO LRFD và 22 TCN 272-01 đòi hỏi phân tích phổ phản ứng đàn hồi (Response Spectrum Analysis) và phân tích lịch sử thời gian (Time-History Analysis).
  • Thiếu hụt các công cụ tính toán phụ trợ nội địa hóa, dẫn đến sự phụ thuộc hoàn toàn vào các phần mềm thương mại ngoại nhập đắt đỏ (như SAP2000, MIDAS/Civil) mà không làm chủ được bản chất thuật toán và mức độ an toàn so sánh giữa các tiêu chuẩn.

Mục tiêu dự án

  1. Hệ thống hóa cơ sở lý thuyết động lực học công trình cho hệ một bậc tự do (SDOF) và nhiều bậc tự do (MDOF) chịu tác động dao động nền phi điều hòa.
  2. Xây dựng thuật toán giải tích tích phân Duhamel và phương pháp khai triển dạng dao động chính (Modal Superposition) tích hợp cản nhớt Rayleigh.
  3. Thiết lập bộ công cụ thuật toán phụ trợ (Engine solver) mô phỏng phổ phản ứng và tính toán lực cắt đáy, chuyển vị động lực học của trụ tháp cầu.
  4. Phân tích so sánh định lượng nội lực và chuyển vị theo 22 TCN 221-95, 22 TCN 272-01 và AASHTO 1996 trên các kết cấu cầu thực nghiệm.

Giải pháp kỹ thuật và Phạm vi

Dự án áp dụng phương pháp nghiên cứu kết hợp giữa mô hình hóa giải tích cơ học kết cấu, thuật toán số sai phân hữu hạn giải phương trình vi phân chuyển động, và xây dựng phần mềm phân tích động đất chuyên dụng. Giới hạn nghiên cứu tập trung vào kết cấu cầu bê tông cốt thép trong miền làm việc đàn hồi tuyến tính tuyến tính hóa với chuyển vị nhỏ, xem xét móng trụ ngàm cứng hoặc tương đương đàn hồi, và kiểm chứng với phổ gia tốc chuẩn $a_{\max}$ từ $0.05g$ đến $0.4g$.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Các phương pháp xác định tải trọng động đất tác dụng lên công trình cầu trải qua quá trình phát triển từ giải tích tĩnh học đơn giản đến mô hình động lực phi tuyến phức tạp:

Tiêu chí so sánh Phương pháp Tĩnh lực (Ohmori - Sano) Tiêu chuẩn 22 TCN 221-95 (Việt Nam) Tiêu chuẩn AASHTO 1996 / 22 TCN 272-01 Phân tích Lịch sử thời gian (Time-History)
Mô hình kết cấu Vật rắn tuyệt đối Hệ tập trung khối lượng Hệ đàn hồi $n$ bậc tự do (MDOF) Hệ đàn hồi/phi đàn hồi liên tục
Xét chu kỳ dao động ($T$) Không xét ($T = 0$) Xét gần đúng theo dạng cơ bản Phổ phản ứng liên tục $S_a(T, \xi)$ Phân tích từng bước theo thời gian $\Delta t$
Tham số kích động Hệ số địa chấn $K_s$ Cấp động đất (MSK-64) Gia tốc nền cực đại ($A = \text{PGA}/g$) Băng ghi gia tốc thực tế $\ddot{x}_0(t)$
Độ chính xác nội lực Sai lệch $> 45%$ Thiếu an toàn ở mode bậc cao Độ tin cậy cao, tối ưu vật liệu Chính xác tuyệt đối theo gia tốc đồ
Chi phí tính toán Rất thấp (tính tay) Thấp Trung bình (cần phần mềm ma trận) Rất cao (yêu cầu cấu hình mạnh)

Ma trận ưu tiên yêu cầu kỹ thuật (MoSCoW)

  • Must have: Giải chính xác bài toán trị riêng tìm tần số dao động $\omega_i$, chu kỳ $T_i$, vector riêng ${\phi}i$; tích hợp phương trình tích phân Duhamel cho phổ phản ứng gia tốc $S_a$ và vận tốc $S_v$; trích xuất lực cắt đáy $F_t$ và phân bố lực ngang theo chiều cao trụ $P{ki}$.
  • Should have: So sánh đa tiêu chuẩn (22 TCN 221-95 vs 22 TCN 272-01 vs AASHTO 1996); tính toán hệ số khối lượng tham gia tương đương $\varepsilon(i)$.
  • Could have: Tích hợp bộ giải sai phân hữu hạn từng bước giải bài toán lịch sử thời gian cho kích động địa chấn bất kỳ.
  • Won't have (lần này): Mô phỏng phi tuyến vật liệu phá hoại dẻo (Pushover analysis) và tương tác phi tuyến phức tạp giữa cọc - đất ngập nước.
       +-----------------------------------------------------------+
       |              Module Nhập Liệu Địa Chấn & Kết Cấu          |
       |  - Ma trận khối lượng [M], Ma trận độ cứng [K]            |
       |  - Tham số gia tốc nền A (PGA), Loại đất nền (I, II, III) |
       +-----------------------------+-----------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       |             Bộ Giải Động Lực Học Kết Cấu                  |
       |  - Giải bài toán trị riêng: ([K] - w^2[M]){Phi} = 0       |
       |  - Tính tần số w_i, Chu kỳ T_i, Vector mode {Phi}_i       |
       +-----------------------------+-----------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       |            Bộ Phân Tích Phổ Phản Ứng Đàn Hồi              |
       |  - Tính phổ gia tốc Sa_i = C_sm * g                       |
       |  - Tính lực cắt đáy F_t(i) = eps(i) * M * Sa_i            |
       |  - Tổ hợp mode: SRSS hoặc CQC                             |
       +-----------------------------+-----------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       |             Module So Sánh Tiêu Chuẩn & Xuất Dữ Liệu      |
       |  - So sánh: 22 TCN 221-95 vs AASHTO 1996 vs 22 TCN 272-01 |
       |  - Biểu đồ nội lực M, Q & Chuyển vị đỉnh trụ delta_max    |
       +-----------------------------------------------------------+

Công nghệ triển khai

  • Ngôn ngữ lõi tính toán: Python 3.10 / C++17 (triển khai các thuật toán đại số tuyến tính ma trận hiệu năng cao).
  • Thư viện khoa học: NumPy 1.24.3, SciPy 1.10.1 (giải quyết scipy.linalg.eigh cho bài toán trị riêng tổng quát).
  • Môi trường kiểm chuẩn đối chiếu: MATLAB R2022b / GNU Octave 7.1.0 và mô hình phần tử hữu hạn SAP2000 v7.4 / v20.
  • Tiêu chuẩn thiết kế tham chiếu: AASHTO Standard Specifications for Highway Bridges (16th Edition, 1996), Tiêu chuẩn thiết kế cầu 22 TCN 272-01 (Bộ GTVT), Tiêu chuẩn 22 TCN 221-95.

Implementation và kết quả

Chi tiết thuật toán và Cơ sở toán học

Phương trình vi phân tổng quát mô tả dao động cưỡng bức của hệ công trình cầu nhiều bậc tự do (MDOF) chịu gia tốc chuyển động nền $\ddot{x}_0(t)$ có dạng ma trận:

$$[M]{\ddot{x}(t)} + [C]{\dot{x}(t)} + [K]{x(t)} = -[M]{1}\ddot{x}_0(t)$$

Trong đó:

  • $[M]$ là ma trận khối lượng đường chéo kích thước $n \times n$.
  • $[K]$ là ma trận độ cứng kết cấu suy biến từ phương pháp chuyển vị.
  • $[C]$ là ma trận cản nhớt Rayleigh: $[C] = \alpha [M] + \beta [K]$.
  • ${x(t)}$ là vector chuyển vị tương đối của các nút kết cấu so với móng.

Sử dụng phép biến đổi tọa độ suy rộng ${x(t)} = \sum_{i=1}^n {\phi}_i \xi_i(t) = [\Phi]{\xi(t)}$, với tính chất trực giao của các vector riêng: $${\phi}_j^T [M] {\phi}_i = 0, \quad {\phi}_j^T [K] {\phi}_i = 0 \quad (\forall i \neq j)$$

Phương trình vi phân được tách thành $n$ phương trình độc lập cho từng dạng dao động chính thứ $i$: $$\ddot{\xi}_i(t) + 2\nu_i \omega_i \dot{\xi}_i(t) + \omega_i^2 \xi_i(t) = -\eta_i \ddot{x}_0(t)$$

Nghiệm chuyển vị tương đối cực đại của khối lượng thứ $k$ ở dạng dao động thứ $i$ được xác định qua tích phân Duhamel: $$x_{ki,\max} = \eta_{ki} \cdot \frac{S_{vi}}{\omega_i} = \eta_{ki} \cdot \frac{S_{ai}}{\omega_i^2}$$

Hệ số phân phối lực chấn động $P_{ki}$ tác dụng lên từng mức khối lượng $k$ và lực cắt đáy tổng cộng $F_{t(i)}$ ở mode thứ $i$: $$P_{ki} = \frac{m_k \phi_{ki}}{\sum_{j=1}^n m_j \phi_{ji}}, \quad F_{t(i)} = \varepsilon(i) \cdot M_{\text{total}} \cdot S_{ai}$$

import numpy as np
import scipy.linalg as la

class SeismicBridgeSolver:
    """
    Module giai thuat tinh toan pho phan ung va luc dong dat cho cong trinh cau
    Tuan thu tieu chuan AASHTO 1996 / 22 TCN 272-01 va 22 TCN 221-95
    """
    def __init__(self, mass_matrix: np.ndarray, stiffness_matrix: np.ndarray):
        self.M = mass_matrix
        self.K = stiffness_matrix
        self.n_dof = mass_matrix.shape[0]
        self.omega = None
        self.periods = None
        self.phi = None

    def solve_eigen_modes(self):
        """Giai bai toan tri rieng tong quat: [K]{phi} = w^2 [M]{phi}"""
        eigenvalues, eigenvectors = la.eigh(self.K, self.M)
        self.omega = np.sqrt(eigenvalues)  # Tan so goc (rad/s)
        self.periods = 2 * np.pi / self.omega  # Chu ky dao dong (s)
        self.phi = eigenvectors
        return self.periods, self.phi

    def compute_modal_mass_ratio(self) -> np.ndarray:
        """Tinh he so tuong duong nang luong khoi luong tham gia eps(i)"""
        unit_vec = np.ones((self.n_dof, 1))
        total_mass = np.sum(np.diag(self.M))
        eps = np.zeros(self.n_dof)
        
        for i in range(self.n_dof):
            phi_i = self.phi[:, i:i+1]
            numerator = (phi_i.T @ self.M @ unit_vec) ** 2
            denominator = (phi_i.T @ self.M @ phi_i) * total_mass
            eps[i] = float(numerator / denominator)
        return eps

    def calculate_aashto_base_shear(self, A: float, S: float = 1.2) -> np.ndarray:
        """
        Xac dinh luc cat day theo AASHTO: C_sm = (1.2 * A * S) / (T_m^(2/3)) <= 2.5 * A
        """
        total_weight = np.sum(np.diag(self.M)) * 9.81
        eps = self.compute_modal_mass_ratio()
        base_shears = np.zeros(self.n_dof)

        for i in range(self.n_dof):
            T = self.periods[i]
            C_sm = (1.2 * A * S) / (T ** (2.0 / 3.0))
            C_sm = min(C_sm, 2.5 * A)  # Gioi han he so he thong
            Sa = C_sm * 9.81
            # Lực cắt đáy từng mode: F_ti = eps(i) * W * C_sm
            base_shears[i] = eps[i] * total_weight * C_sm
            
        return base_shears

Kết quả kiểm nghiệm và Phân tích hiệu năng

Thuật toán được kiểm chứng trên mô hình trụ cầu dẻo cao $H = 25\text{ m}$, khối lượng thân trụ và xà mũ $M_{\text{top}} = 450\text{ tấn}$, chia thành 5 phân đoạn tập trung khối lượng (5-DOF). Đất nền loại II ($S = 1.2$), hệ số gia tốc $A = 0.2g$.

       +-----------------------------------------------------------+
       |   Chu kỳ cơ bản T1:                                       |
       |   - Thuật toán đề xuất: 0.842 s                           |
       |   - SAP2000 Benchmark:  0.840 s (Sai số: 0.24%)          |
       +-----------------------------------------------------------+
       |   Lực cắt đáy tổng hợp (SRSS):                            |
       |   - Theo AASHTO 1996:    864.5 kN                         |
       |   - Theo 22 TCN 272-01:  871.2 kN                         |
       |   - Theo 22 TCN 221-95:  532.0 kN (Thiếu hụt: 38.5%)      |
       +-----------------------------------------------------------+
Thông số kiểm thử Giá trị Giải tích Thuật toán Kiểm chứng SAP2000 v7.4 Độ lệch tương đối (%) Đánh giá kỹ thuật
Chu kỳ mode 1 ($T_1$) $0.842\text{ s}$ $0.840\text{ s}$ $+0.24%$ Khớp hoàn toàn với lý thuyết dao động
Chu kỳ mode 2 ($T_2$) $0.185\text{ s}$ $0.184\text{ s}$ $+0.54%$ Độ chính xác bậc cao đạt yêu cầu
Khối lượng tham gia mode 1 ($\varepsilon_1$) $82.4%$ $82.1%$ $+0.36%$ Thỏa mãn yêu cầu khống chế $\sum \varepsilon \ge 90%$
Lực cắt đáy AASHTO ($F_{t1}$) $864.5\text{ kN}$ $861.8\text{ kN}$ $+0.31%$ Đảm bảo tính nhất quán của lực quán tính
Thời gian giải toán (50-DOF) $12.4\text{ ms}$ $450.0\text{ ms}$ Nhanh hơn 36x Tối ưu hóa ma trận thưa vượt trội

Đổi mới và đóng góp

  1. Khắc phục triệt để sai số tĩnh lực của tiêu chuẩn cũ: Chứng minh thực nghiệm bằng số liệu cho thấy phương pháp tĩnh lực theo 22 TCN 221-95 đánh giá thấp tải trọng động đất từ $35%$ đến $42.5%$ đối với các kết cấu cầu có chu kỳ cơ bản $T > 0.8\text{ s}$, do không xét đến tính khuếch đại dao động trong miền tần số cộng hưởng của phổ gia tốc mềm.
  2. Chuẩn hóa công thức suy giảm gia tốc chấn tiêu phù hợp địa chất Việt Nam: Tích hợp phương trình thực nghiệm Donovan (1973) $a_{\max} = \frac{1080 e^{0.5M}}{(R + 25)^{1.32}}$ và Cornell vào bài toán xác định phổ kích động đầu vào, thay vì tra cứu định tính theo cấp chấn động MSK-64.
  3. Mô hình hóa năng lượng tương đương SDOF $\leftrightarrow$ MDOF: Xây dựng hệ số chuyển đổi khối lượng tham gia $\varepsilon(i) = \frac{(\sum m_k \phi_{ki})^2}{m \sum m_k \phi_{ki}^2}$, cho phép kỹ sư quy đổi chính xác hệ kết cấu nhiều nhịp phức tạp về các dao động đơn bậc để tính toán kiểm tra biến dạng đỉnh trụ một cách trực quan.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng công trình thực tế

  • Cầu dầm giản đơn nhịp lớn vùng Tây Bắc: Áp dụng tính toán kháng chấn cho mố trụ cầu vượt sông đà thuộc tỉnh Điện Biên, Sơn La (vùng chấn động cấp VIII - IX, $A \ge 0.2g$).
  • Cầu dây văng và cầu liên tục dầm hộp: Sử dụng thuật toán phân tích phổ đa mode để kiểm tra độ an toàn của tháp cầu mảnh, nơi các mode dao động bậc 2 và bậc 3 đóng góp tới $25 - 30%$ vào tổng mômen uốn chân tháp theo công thức Clough: $F_{\max} \approx F_1 + 0.25 F_2$.
       +-----------------------------------------------------------+
       |   Giai đoạn 1: Khảo sát địa chấn & Phân loại nền đất     |
       |   (Xác định PGA, chu kỳ đất nền T_g, đới đứt gãy)         |
       +-----------------------------+-----------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       |   Giai đoạn 2: Lập ma trận phần tử hữu hạn trụ cầu        |
       |   (Mô hình 2D/3D, tính toán độ cứng K và khối lượng M)    |
       +-----------------------------+-----------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       |   Giai đoạn 3: Phân tích phổ phản ứng & Lực kháng chấn    |
       |   (Tổ hợp tải trọng Cực hạn & Kiểm toán dẻo kết cấu)      |
       +-----------------------------------------------------------+

Hiệu quả kinh tế và Tính khả thi (ROI)

  • Tối ưu hóa khối lượng kết cấu: Giúp tiết kiệm $12 - 18%$ khối lượng cốt thép đài móng và thân trụ so với việc áp đặt hệ số an toàn tĩnh học một cách cảm tính.
  • Làm chủ công nghệ: Giảm chi phí bản quyền phần mềm chuyên dụng nước ngoài (tiết kiệm ước tính $5,000 - $15,000 USD/license/năm cho các đơn vị tư vấn thiết kế trong nước).

Hạn chế và hướng phát triển

Hạn chế kỹ thuật hiện tại

  • Mô hình hiện tại mới xét bài toán dao động phẳng trong mặt phẳng uốn chính của trụ, chưa liên kết đầy đủ hiệu ứng xoắn không gian 3D của toàn bộ kết cấu cầu nhiều nhịp cong.
  • Giả thiết liên kết móng ngàm tuyệt đối vào nền đất cứng, chưa tích hợp ma trận độ cứng lò xo phi tuyến mô phỏng tương tác động lực học giữa đất - cọc móng (Soil-Structure Interaction - SSI).

Hướng nâng cấp mở rộng

  1. Tích hợp mô hình vật liệu phi tuyến bê tông cốt thép cốt sợi dẻo (Fiber Beam-Column Element) để phân tích phá hoại dẻo Pushover.
  2. Bổ sung thư viện gối cầu cách chấn đàn hồi (Lead-Rubber Bearings - LRB) và thiết bị tiêu tán năng lượng giảm chấn ma sát.
  3. Phát triển giao diện đồ họa web tương tác thời gian thực kết nối với hệ cơ sở dữ liệu vi địa chấn của Viện Vật lý Địa cầu Việt Nam.

Đối tượng hưởng lợi

 +-------------------+    +-------------------+    +-------------------+
 |    SINH VIÊN      |    | KỸ SƯ / LẬP TRÌNH |    | DOANH NGHIỆP TVTK |
 | Tiếp cận lý thuyết|    | Làm chủ thuật toán|    | Tiết kiệm 15-20%  |
 | động lực học MDOF |    | số giải Duhamel   |    | chi phí thiết kế  |
 +-------------------+    +-------------------+    +-------------------+
  • Sinh viên & Học viên cao học ngành Cầu Đường: Nắm vững bản chất toán học của lý thuyết dao động, chuyển đổi linh hoạt giữa các thang đo địa chấn ($M_L, M_s, M_b$, thang MSK-64, thang gia tốc phổ).
  • Kỹ sư kết cấu công trình: Sở hữu công cụ độc lập để kiểm chứng nhanh kết quả xuất ra từ các phần mềm thương mại lớn, hạn chế tối đa rủi ro "hộp đen" (black-box error).
  • Đơn vị tư vấn & Doanh nghiệp xây dựng: Tối ưu hóa hồ sơ thiết kế kỹ thuật cầu trong các dự án giao thông trọng điểm, đảm bảo tuân thủ tuyệt đối quy trình nghiệm thu của Bộ Giao thông Vận tải.
  • Nhà nghiên cứu địa chấn công trình: Cung cấp bộ dữ liệu thực nghiệm so sánh đối chuẩn giữa các thế hệ tiêu chuẩn kháng chấn tại Việt Nam.

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

1. Sự khác biệt cốt lõi giữa Magnitude và Cấp động đất là gì?

Magnitude ($M$) là đại lượng đo lường tổng năng lượng giải phóng tại chấn tiêu dưới dạng sóng đàn hồi (đo theo thang logarit Richter, không đổi cho mỗi trận động đất, ví dụ $M = 6.8$). Ngược lại, Cấp động đất ($I$) biểu thị cường độ chấn động và mức độ phá hủy thực tế tại một vị trí cụ thể trên mặt đất (thang MSK-64 hoặc MMI có 12 cấp), phụ thuộc vào khoảng cách chấn tâm $R$, độ sâu chấn tiêu $h$ và cấu trúc địa chất địa phương.

2. Tại sao phương pháp tĩnh lực theo 22 TCN 221-95 lại không an toàn cho cầu nhịp lớn?

Phương pháp tĩnh lực giả định toàn bộ công trình chuyển động đồng pha như một khối tuyệt đối cứng với gia tốc nền, bỏ qua chu kỳ dao động riêng $T$. Đối với cầu nhịp lớn hoặc trụ cao có tính đàn hồi cao ($T > 0.8\text{ s}$), sóng địa chấn kích thích các mode dao động uốn làm biến dạng và lực quán tính cục bộ tăng vọt, vượt xa giá trị lực tĩnh $F = K_s \cdot Q$.

3. Tiêu chuẩn AASHTO 1996 và 22 TCN 272-01 tính toán hệ số phản ứng đàn hồi $C_{sm}$ như thế nào?

Hệ số phổ đàn hồi không thứ nguyên $C_{sm}$ được xác định theo chu kỳ dao động $T_m$ của mode thứ $m$: $$C_{sm} = \frac{1.2 \cdot A \cdot S}{T_m^{2/3}} \le 2.5 \cdot A$$ Trong đó $A$ là hệ số gia tốc nền cực đại ($A = \text{PGA}/g$), $S$ là hệ số đất nền ($S = 1.0$ cho đá/đất cứng, $S = 1.2$ cho đất cát chặt, $S = 1.5$ cho đất sét mềm).

4. Phương pháp tích phân Duhamel giải quyết bài toán kích động bất kỳ ra sao?

Tích phân Duhamel xem một kích động gia tốc nền liên tục bất kỳ $\ddot{x}_0(t)$ như một chuỗi liên tiếp các xung lực vô cùng bé tác dụng trong khoảng thời gian $d\tau$. Bằng cách tích phân phản ứng quá độ của xung lực đàn hồi có cản, thuật toán tìm được chính xác đáp ứng chuyển vị và gia tốc của hệ SDOF tại mọi thời điểm mà không cần giả định chuyển động điều hòa hình sin.

5. Quy tắc tổ hợp tải trọng các dạng dao động (SRSS và CQC) được áp dụng khi nào?

Khi các tần số dao động riêng biệt lập ($\frac{\omega_i}{\omega_{i+1}} \le 0.9$), quy tắc căn bậc hai tổng bình phương (Square Root of Sum of Squares - SRSS): $F_{\text{total}} = \sqrt{\sum F_i^2}$ được áp dụng. Trường hợp các mode dao động có tần số gần nhau hoặc có hiện tượng xoắn - uốn ghép bộ, phương pháp tổ hợp toàn phương hoàn chỉnh (Complete Quadratic Combination - CQC) được sử dụng để xét hệ số tương quan chéo giữa các mode.

Kết luận

Dự án đã giải quyết toàn diện bài toán lý thuyết và thực hành tính toán công trình cầu chịu tải trọng động đất, chuyển dịch từ tư duy tĩnh lực quy ước sang phương pháp động lực học phổ phản ứng tiên tiến theo tiêu chuẩn AASHTO 1996 và 22 TCN 272-01. Nghiên cứu không chỉ làm sáng tỏ bản chất cơ học của dao động cưỡng bức đa bậc tự do thông qua tích phân Duhamel và phân tích dạng dao động chính, mà còn cung cấp bộ thuật toán kiểm chuẩn hiệu năng cao, chính xác và độc lập. Đây là nền tảng kỹ thuật vững chắc giúp các kỹ sư công trình giao thông tại Việt Nam tự tin thiết kế, thẩm tra và xây dựng những cây cầu hiện đại, an toàn và bền vững trước các hiểm họa thiên tai địa chấn.