Tổng quan về luận án

Nghiên cứu động lực học công trình cầu dưới tác động của đoàn xe di động là một trong những bài toán phức tạp và then chốt nhất của chuyên ngành Kỹ thuật Xây dựng Cầu - Hầm (Mã số: 62.05). Trong bối cảnh hạ tầng giao thông hiện đại hóa với mật độ lưu thông dày đặc, tải trọng trục gia tăng và dải vận tốc khai thác mở rộng từ 20 km/h đến 120 km/h, các hiệu ứng động lực học phát sinh không chỉ làm suy giảm tuổi thọ mỏi của kết cấu mà còn đe dọa trực tiếp đến an toàn khai thác. Luận án tiến sĩ kỹ thuật của tác giả Vũ Văn Toản (2017) với tiêu đề "Phân tích dao động của kết cấu nhịp cầu dưới tác động của hoạt tải khai thác có xét đến độ mấp mô mặt cầu", dưới sự hướng dẫn khoa học của GS.TS Nguyễn Viết Trung và PGS.TS Trần Đức Nhiệm tại Trường Đại học Giao thông Vận tải, đã giải quyết căn bản bài toán tương tác phức hợp giữa phương tiện và công trình cầu (Vehicle-Bridge Interaction - VBI).

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở chỗ các tiêu chuẩn thiết kế cầu hiện hành như AASHTO LRFD hay tiêu chuẩn ngành tại Việt Nam chủ yếu quy định hệ số xung kích/hệ số động lực $(1+IM)$ dưới dạng một giá trị tĩnh biểu kiến hoặc phụ thuộc đơn thuần vào chiều dài nhịp, bỏ qua đặc tính ngẫu nhiên của độ mấp mô mặt cầu và sự tương tác động cơ học phụ thuộc thời gian giữa hệ thống treo của xe với kết cấu dầm. Luận án đã giải quyết bài toán này thông qua việc thiết lập 4 câu hỏi nghiên cứu và giả thuyết khoa học:

  1. RQ1: Làm thế nào để mô hình hóa chính xác hàm kích động biên dạng mặt đường ngẫu nhiên theo tiêu chuẩn quốc tế ISO 8608 vào bài toán tương tác dầm - xe?
  2. RQ2: Mức độ ảnh hưởng tương hỗ giữa các thông số động lực học của hệ thống treo ô tô (độ cứng $k_s$, cản $c_s$, độ cứng lốp $k_t$, cản lốp $c_t$) đối với phản ứng dao động uốn của dầm giản đơn và dầm liên tục?
  3. RQ3: Giá trị hệ số giảm chấn thực tế $(\xi)$ của các công trình cầu khai thác tại Việt Nam là bao nhiêu khi đo đạc bằng phương pháp suy giảm dao động tự do?
  4. RQ4: Mối quan hệ định lượng giữa vận tốc xe chạy (20 - 120 km/h), cấp độ mấp mô mặt cầu (từ rất tốt đến rất xấu) và hệ số động lực $(1+IM)$ biến thiên như thế nào?

Khung lý thuyết của luận án được xây dựng trên nền tảng cơ học kết cấu liên tục, nguyên lý biến phân Hamilton, phương pháp phần tử hữu hạn (FEM) và lý thuyết quá trình ngẫu nhiên. Luận án đã lượng hóa được hệ số động lực trên các nhịp dầm bê tông dự ứng lực (BTDUL) và dầm thép liên hợp bản bê tông cốt thép (BTCT), cung cấp bộ số liệu thực nghiệm đo đạc hệ số cản tại cầu Đa Phước, Cầu vượt ngã tư Thủ Đức và Cầu vượt ngã ba Vũng Tàu.

Literature Review và Positioning

Lịch sử nghiên cứu tải trọng di động khởi nguồn từ thảm họa sập cầu đường sắt Chester tại Anh năm 1847, thúc đẩy các công trình kinh điển của Willis (1849) với phương trình vi phân chuyển động dầm chịu tải khối lượng di trượt: $$\frac{d^2w}{dx^2} = \frac{g}{v^2} \left( \frac{w}{1 - \frac{M_p v^2 (l - x) x^2}{3EJ l}} \right)$$ và sau đó là các hiệu chỉnh hệ số của Busins (1883), S. Iliaxevic, F. Krưlov và Meizel (1930). Trong thập niên 1970, A. Philipov và A. Barchenkov đã phát triển các mô hình giải tích cho hệ dầm không gian nhưng gặp rào cản lớn về năng lực tính toán. Về phía động lực học phương tiện, Mitschke (1970) với tác phẩm "Dynamik der Fahrzeuge" và Schiehlen (1980) với lý thuyết hệ nhiều vật (Multibody Dynamics) đã đặt nền móng mô hình hóa dao động ô tô, trong khi Karnopp (1974) phát triển hệ thống điều khiển giảm chấn skyhook và ground-hook.

flowchart TD
    A["Nghiên cứu Cổ điển: Willis (1849), Iliaxevic, Krưlov<br><i>(Tải trọng di động đơn giản, bỏ qua khối lượng/mấp mô)</i>"] --> B["Thập niên 1970-1980: Frýba (1972), Philipov, Mitschke<br><i>(Mô hình 2 khối lượng, mấp mô điều hòa h = h₀ sin ωx)</i>"]
    B --> C["Thập niên 2000-2012: S. Zhu (2004), Lombaert & Conte (2012)<br><i>(Chồng chất mode, mô hình xe 1/4, 2 bậc tự do)</i>"]
    C --> D["Luận án Vũ Văn Toản (2017)<br><i>(Tương tác VBI đầy đủ: Xe 2-3 trục, Hermite FEM, ISO 8608 PSD, In-situ Damping)</i>"]

Tranh luận học thuật cốt lõi tồn tại giữa hai quan điểm tiếp cận:

  • Trường phái uncoupled (tách rời): Tính toán dao động xe dựa trên mặt đường cứng tuyệt đối, sau đó áp lực động tác dụng tuần tự lên cầu (S. Zhu, 2004). Cách tiếp cận này giảm khối lượng tính toán nhưng làm mất đi sự thích ứng ngược từ độ võng dầm lên hệ thống treo.
  • Trường phái coupled (tương tác đồng thời): Hệ xe và cầu là một thể thống nhất liên tục có điều kiện biên thay đổi theo thời gian thực (Ladislav Frýba, 1972; Kwasniewski, 2006).

Luận án định vị nghiên cứu vượt qua các hạn chế quốc tế: Trong khi Frýba (1972) chỉ giải bài toán với mấp mô có quy luật chu kỳ $h = h_0 \sin(\omega x)$ và giải bằng phương pháp Runge-Kutta-Nyström, còn Lombaert & Conte (2012) giới hạn ở mô hình xe 1/4 với 2 bậc tự do trên dầm giản đơn, nghiên cứu này đã thiết lập mô hình tương tác đầy đủ cho xe 2 trục (Maz 5551) và xe 3 trục (Zil 131) chuyển động trên dầm nhiều nhịp liên tục với kích động mấp mô ngẫu nhiên trải rộng theo phổ mật độ công suất PSD của ISO 8608.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng lý thuyết đàn hồi động lực học công trình thông qua việc tích hợp trực tiếp nguyên lý biến phân chuyển vị Hamilton vào hệ tương tác dầm liên tục - xe nhiều trục. Trích dẫn phương trình tương tác cơ sở từ tài liệu nguồn của Frýba (1972): $$EJ \frac{\partial^4 v(x, t)}{\partial x^4} + 2b \frac{\partial v(x, t)}{\partial t} + m \frac{\partial^2 v(x, t)}{\partial t^2} = \sum_{i=1}^{2} \delta(x - x_i) R_i(t)$$ Trong đó lực tương tác từ trục xe thứ $i$ lên dầm tại vị trí $x_i$ được xác định theo công thức liên kết động lực học: $$R_i(t) = K_i u_i(t) = K_i [v_i(t) - \delta_i v(x_i, t) - r_i(x_i)]$$

Nghiên cứu đã chuyển dịch mô hình lý thuyết từ dạng giải tích cổ điển sang mô hình ma trận phần tử hữu hạn tổng quát có cản: $$M \ddot{q}(t) + C \dot{q}(t) + K q(t) = p(t)$$ với vector tải trọng nút $p(t) = \int_0^l N_w^T f(x,t) dx + \sum N_w^T P_i(t)$, cho phép phân rã các hiệu ứng cản nhớt phân bố dọc dầm $c$ và cản cục bộ của phương tiện.

classDiagram
    class VehicleDynamics {
        +Double ms (Khối lượng thân xe)
        +Double mu (Khối lượng không treo)
        +Double ks, cs (Độ cứng & cản hệ treo)
        +Double kt, ct (Độ cứng & cản lốp xe)
        +Double Jy (Mômen quán tính lắc dọc)
        +calculate_axle_force()
    }
    class BridgeStructure {
        +Matrix M (Ma trận khối lượng Hermite)
        +Matrix K (Ma trận độ cứng uốn EJ)
        +Matrix C (Ma trận cản Rayleigh)
        +Vector q (Vector chuyển vị nút w, theta)
        +solve_FEM()
    }
    class SurfaceRoughness {
        +Double IRI (Chỉ số gồ ghề quốc tế)
        +Function PSD_ISO8608 (Mật độ phổ công suất)
        +generate_profile_r(x)
    }
    VehicleDynamics --> BridgeStructure : Lực tương tác động fT(t)
    BridgeStructure --> VehicleDynamics : Độ võng w(x,t) & Gia tốc
    SurfaceRoughness --> VehicleDynamics : Kích động mấp mô r(x)

Khung phân tích độc đáo

Khung phân tích của luận án kết hợp 3 trụ cột lý thuyết:

  1. Lý thuyết dầm Euler-Bernoulli rời rạc hóa bằng hàm nội suy Hermite bậc 3: $$N_1(x) = 1 - 3\frac{x^2}{l^2} + 2\frac{x^3}{l^3}; \quad N_2(x) = x - 2\frac{x^2}{l} + \frac{x^3}{l^2}$$ $$N_3(x) = 3\frac{x^2}{l^2} - 2\frac{x^3}{l^3}; \quad N_4(x) = -\frac{x^2}{l} + \frac{x^3}{l^2}$$
  2. Lý thuyết động lực học hệ nhiều vật mô tả chuyển vị thẳng đứng $z_s$, lắc dọc $\theta_s$ của thân xe và chuyển vị các trục không treo $z_{ui}$.
  3. Lý thuyết trường ngẫu nhiên đồng nhất biểu diễn mật độ phổ chiều cao mấp mô $S_u(\Omega) = S_u(\Omega_0) (\frac{\Omega}{\Omega_0})^{-n_1}$.

Điều kiện biên được xác lập chặt chẽ: liên kết gối tựa cứng hoặc đàn hồi, vận tốc chuyển động $v = \text{const}$, bánh xe duy trì tiếp xúc liên tục với mặt cầu ($u_i(t) \ge 0$).

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

Thiết kế nghiên cứu

Luận án áp dụng thế giới quan thực chứng (Positivism) kết hợp giữa phương pháp mô phỏng số học (Computational Dynamics) và kiểm chứng thực nghiệm hiện trường (In-situ Structural Testing). Thiết kế nghiên cứu đa tầng bao gồm:

  • Tầng 1 (Cơ sở số): Rời rạc hóa kết cấu nhịp bằng $N$ phần tử hữu hạn dầm chịu uốn thuần túy.
  • Tầng 2 (Mô phỏng hệ thống): Xây dựng khối mô phỏng trên nền tảng Matlab - Simulink cho xe 1/4, xe 2 trục và xe 3 trục.
  • Tầng 3 (Thực nghiệm kiểm chứng): Đo đạc dao động thực tế tại công trình cầu Đa Phước, Cầu vượt Thủ Đức và Ngã ba Vũng Tàu.
flowchart LR
    A["Mô hình hóa PTHH Dầm<br><i>(Ma trận K, M, C)</i>"] --> D["Thuật toán tích phân Newmark & Giải lặp cải tiến"]
    B["Mô hình động lực xe<br><i>(1/4, 2 trục Maz, 3 trục Zil)</i>"] --> D
    C["Biên dạng mấp mô ngẫu nhiên<br><i>(ISO 8608 PSD)</i>"] --> D
    D --> E["Phản ứng động lực<br><i>(Độ võng w(t), Ứng suất, Hệ số 1+IM)</i>"]
    E <--> F["Thực nghiệm hiện trường<br><i>(Đo gia tốc, ứng suất cầu Đa Phước)</i>"]

Quy trình nghiên cứu rigorous

Quy trình giải bài toán tương tác cầu - xe được thực hiện theo thuật toán lặp cải tiến tại mỗi bước tích phân thời gian $\Delta t$ của phương pháp Newmark-beta ($\gamma = 0.5, \beta = 0.25$):

BƯỚC 1: Khởi tạo điều kiện ban đầu: q(0), dq(0)/dt, ddq(0)/dt, zs(0), dzs(0)/dt.
BƯỚC 2: Tại bước thời gian t_{j+1} = t_j + dt:
   a. Xác định tọa độ vị trí các trục bánh xe: x_i = v * t_{j+1}.
   b. Trích xuất cao độ mấp mô ngẫu nhiên mặt đường r(x_i) từ hàm phổ PSD ISO 8608.
   c. Tính chuyển vị và vận tốc bề mặt dầm tại điểm tiếp xúc: w(x_i, t) = N_w(x_i) * q(t).
   d. Kích động tổng hợp tác dụng lên bánh xe: u_i = w(x_i, t) + r(x_i).
   e. Giải phương trình vi phân dao động xe trong Simulink -> Xác định phản lực động f_T(t).
   f. Đặt f_T(t) vào vector tải trọng nút p(t) của phương trình kết cấu dầm.
   g. Giải phương trình dầm bằng thuật toán Newmark -> Thu được q(t_{j+1}).
   h. Kiểm tra hội tụ sai số chuyển vị: ||q^{(k+1)} - q^{(k)}|| / ||q^{(k)}|| < ε (10^-4).
BƯỚC 3: Lặp lại cho đến khi xe vượt hoàn toàn qua chiều dài kết cấu nhịp L.

Data và phân tích

Số liệu đầu vào bao gồm:

  • Xe 2 trục (Maz 5551): Tổng khối lượng $M_v = 16.0$ tấn, phân bố khối lượng thân xe $m_s$, khối lượng cầu trước/sau $m_{u1}, m_{u2}$, độ cứng nhíp $k_{s1}, k_{s2}$, cản chấn $c_{s1}, c_{s2}$, độ cứng lốp $k_{t1}, k_{t2}$.
  • Xe 3 trục (Zil 131): Tổng khối lượng $M_v = 10.4$ tấn, 3 trục liên kết phụ thuộc qua hệ thống cân bằng.
  • Thiết bị thu thập hiện trường: Cảm biến gia tốc áp điện (Piezoelectric Accelerometers), thiết bị đo biến dạng động (Dynamic Strain Gauges), hệ thống thu thập số liệu đa kênh tần số lấy mẫu 1000 Hz.
  • Công thức thực nghiệm xác định hệ số giảm chấn từ độ giảm biên độ logarit $(\delta)$: $$\delta = \ln\left(\frac{A_k}{A_{k+1}}\right) = \frac{2\pi \xi}{\sqrt{1 - \xi^2}} \implies \xi \approx \frac{\delta}{2\pi}$$

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Sự chi phối áp đảo của độ mấp mô ngẫu nhiên: Độ mấp mô mặt cầu là nguồn kích động chủ đạo làm bùng phát lực động tương tác. So với mặt cầu phẳng tuyệt đối, khi mặt cầu có mấp mô ngẫu nhiên loại B-C (theo ISO 8608), độ võng động giữa nhịp tăng từ 18% đến 42% tùy thuộc vào vận tốc xe.
  2. Đặc tính phi tuyến của hệ số động lực theo dải vận tốc: Hệ số động lực $(1+IM)$ không tăng tuyến tính theo vận tốc xe $v \in [20, 120]\text{ km/h}$ mà xuất hiện các dải vận tốc cộng hưởng cục bộ nguy hiểm. Đối với nhịp dầm giản đơn BTDUL $L = 33\text{ m}$, các đỉnh dao động cực đại xuất hiện rõ nét ở dải tốc độ $60 - 80\text{ km/h}$.
  3. Hệ số giảm chấn thực tế tại công trình cầu Việt Nam: Luận án chỉ ra rằng giá trị hệ số cản $\xi$ thực tế của các nhịp dầm giản đơn BTCT dao động trong khoảng $1.5% - 2.8%$, dầm liên tục thép - bê tông đạt $1.1% - 2.2%$, thấp hơn đáng kể so với giả thiết $5%$ thường áp dụng tùy tiện trong các tính toán lý thuyết trước đây.
  4. Hiệu ứng tương hỗ giữa các trục xe: Đối với xe 3 trục (Zil 131), kích động do trục trước dằn xóc trên mấp mô tạo ra dao động lắc dọc $(\theta_s)$ truyền động năng làm tăng áp lực động tức thời lên các trục sau, dẫn đến chuỗi xung kích có biên độ lớn hơn so với xe 2 trục đơn lẻ.
  5. Kiểm chứng thực tế tại Cầu Đa Phước: Sai số giữa phản ứng độ võng động tính toán từ mô hình FEM - Simulink cải tiến so với số liệu đo đạc thực tế tại hiện trường chỉ dao động trong khoảng $4.8% - 7.6%$, khẳng định độ tin cậy tuyệt đối của thuật toán giải lặp.
BẢNG TỔNG HỢP SO SÁNH HỆ SỐ ĐỘNG LỰC (1+IM) TRÊN DẦM GIẢN ĐƠN BTDUL (L = 33m)
+----------------+----------------+----------------+----------------+----------------+
| Vận tốc (km/h) | Mặt phẳng lỳ   | ISO Loại A     | ISO Loại B     | ISO Loại C     |
+----------------+----------------+----------------+----------------+----------------+
| 20 km/h        | 1.032          | 1.058          | 1.114          | 1.205          |
| 40 km/h        | 1.065          | 1.092          | 1.168          | 1.284          |
| 60 km/h        | 1.118          | 1.154          | 1.252          | 1.391          |
| 80 km/h        | 1.142          | 1.186          | 1.295          | 1.446          |
| 100 km/h       | 1.125          | 1.161          | 1.264          | 1.412          |
| 120 km/h       | 1.138          | 1.179          | 1.288          | 1.435          |
+----------------+----------------+----------------+----------------+----------------+

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung luận cứ khoa học về cơ chế kích động tham số ngẫu nhiên trong bài toán động lực học công trình giao thông.
  • Về mặt phương pháp luận: Cung cấp một framework mô phỏng liên kết giữa cơ sở dữ liệu số FEM và giải thuật điều khiển mô phỏng Matlab-Simulink, dễ dàng mở rộng cho các bài toán phức tạp khác.
  • Về mặt thực tiễn thiết kế & thẩm định: Cảnh báo các cơ quan thiết kế không nên coi hệ số xung kích $(1+IM) = 1.25$ hoặc $1.30$ là an toàn trong mọi điều kiện; khi mặt đường bị xuống cấp (ổ gà, lượn sóng, IRI $> 4.0\text{ mm/m}$), hệ số động lực thực tế có thể vượt ngưỡng $1.45$, tiềm ẩn nguy cơ nứt và phá hoại mỏi bê tông.
  • Về chính sách bảo trì đường bộ: Định lượng rõ giá trị kinh tế của công tác duy tu êm thuận mặt cầu; việc duy trì mặt đường ở cấp độ A-B giúp giảm $20% - 30%$ áp lực động lực học truyền xuống kết cấu chịu lực chính.

Limitations và Future Research

Mặc dù đạt được những bước tiến mang tính đột phá, công trình vẫn tồn tại một số giới hạn biên:

  1. Mô hình kết cấu dầm được phân tích chủ yếu trong mặt phẳng đứng (2D), chưa khảo sát đồng thời hiệu ứng xoắn không gian và tương tác uốn - xoắn khi xe chạy lệch làn hoặc trên cầu cong có lực ly tâm.
  2. Giả thiết tiếp xúc bánh xe - mặt đường xem như liên tục không tách rời ($u_i(t) \ge 0$), chưa xét đến hiện tượng xe bị nảy bứt khỏi mặt cầu ở các gờ mấp mô đột biến với vận tốc rất cao.
  3. Vật liệu kết cấu được giả thiết đàn hồi tuyến tính thuần túy, chưa tích hợp mô hình suy giảm độ cứng phi tuyến do tích lũy nứt nẻ và mỏi vật liệu theo thời gian.

Các hướng nghiên cứu tiếp nối trong tương lai:

  • Mở rộng mô hình tương tác VBI 3D không gian cho các dạng cầu dây văng, cầu vòm và cầu treo nhịp lớn chịu tác động đồng thời của hoạt tải xe và tải trọng gió bão.
  • Nghiên cứu quá trình phanh khẩn cấp và tăng tốc của đoàn xe trên cầu (xét lực hãm dọc trục).
  • Tích hợp công nghệ cảm biến IoT và trí tuệ nhân tạo (AI) để quan trắc sức khỏe công trình (SHM) dựa trên dữ liệu phản ứng dao động thực tế theo thời gian thực.

Tác động và ảnh hưởng

Luận án tạo ra tác động học thuật và thực tiễn sâu rộng:

  • Học thuật: Đóng góp tài liệu tham khảo chuẩn mực cho các chương trình đào tạo tiến sĩ, thạc sĩ chuyên ngành Công trình Cầu - Hầm, mở ra hướng nghiên cứu động lực học ngẫu nhiên trong xây dựng công trình giao thông.
  • Chuyển giao công nghệ: Quy trình tính toán và chương trình mô phỏng Simulink phát triển trong luận án có khả năng ứng dụng trực tiếp tại các viện nghiên cứu, tổng công ty tư vấn thiết kế giao thông (TEDI, TEDISOUTH) để kiểm tra đánh giá năng lực chịu tải của các cầu cũ đang khai thác.
  • Tác động chính sách: Cung cấp cơ sở khoa học để Tổng cục Đường bộ Việt Nam và Bộ Giao thông Vận tải hoàn thiện các quy định về giới hạn chỉ số gồ ghề quốc tế (IRI) trên mặt cầu khai thác tốc độ cao.

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

  • Nghiên cứu sinh & Giảng viên đại học: Tiếp cận hệ phương trình vi phân chuyển động liên tục - rời rạc hóa và quy trình mô phỏng tương tác VBI hoàn chỉnh.
  • Kỹ sư thiết kế công trình cầu: Nắm bắt quy luật phân bố tải trọng động, dải vận tốc nguy hiểm để lựa chọn kích thước hình học, độ cứng chống uốn $EJ$ và giải pháp kết cấu giảm chấn tối ưu.
  • Đơn vị kiểm định & Quản lý bảo trì cầu: Sở hữu phương pháp thực nghiệm đo đạc hệ số giảm chấn $\xi$ và công cụ đánh giá mức độ suy thoái an toàn của kết cấu dựa trên chất lượng bề mặt mặt cầu thực tế.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc mở rộng lý thuyết tương tác dao động động lực học VBI của Ladislav Frýba (1972) và Lombaert & Conte (2012) từ mô hình kích động điều hòa đơn giản sang mô hình kích động ngẫu nhiên toàn diện (ISO 8608 PSD) cho hệ nhiều vật đa trục (xe 2 trục và 3 trục) tương tác với dầm liên tục có cản nhớt phân bố, giải quyết bằng sự kết hợp giữa nguyên lý biến phân Hamilton và phương pháp PTHH.

2. Phương pháp nghiên cứu có điểm gì cải tiến so với các công bố quốc tế trước đây?

So với nghiên cứu của S. Zhu (2004) chỉ sử dụng phương pháp tách rời tuần tự uncoupled và Lombaert (2012) chỉ khảo sát xe 1/4 (2 bậc tự do), luận án đã thiết lập thuật toán lặp cải tiến trực tiếp theo thời gian thực (Time-step iteration) kết hợp giữa phương pháp tích phân Newmark-beta và sơ đồ mô phỏng Simulink hệ nhiều vật, cho phép phản hồi liên tục chuyển vị dầm lên từng trục bánh xe.

3. Phát hiện bất ngờ nhất từ số liệu phân tích của luận án là gì?

Phát hiện bất ngờ nhất là tại một số dải vận tốc trung bình ($60 - 80\text{ km/h}$), hệ số động lực $(1+IM)$ ghi nhận được lại cao hơn so với khi xe chạy ở tốc độ tối đa ($120\text{ km/h}$). Nguyên nhân do sự trùng lùng giữa tần số kích động của bước sóng mấp mô mặt cầu ở dải tốc độ đó với tần số dao động riêng bậc 1 của dầm, gây ra hiện tượng giả cộng hưởng cục bộ.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) không?

Có. Luận án mô tả chi tiết toàn bộ các khối thuật toán trong Matlab - Simulink, thông số cơ lý tính chi tiết của xe Maz 5551 và Zil 131 (Bảng 4.1, 4.2, 4.3), ma trận phần tử dầm Hermite, phương pháp đo cao độ mặt cầu và quy trình bố trí đầu đo gia tốc, đo biến dạng động tại hiện trường cầu Đa Phước.

5. Định hướng nghiên cứu 10 năm được vạch ra như thế nào?

Định hướng tập trung vào việc số hóa mô hình VBI 3D phi tuyến, tích hợp ảnh hưởng của mỏi vật liệu kết cấu, phát triển các thiết bị giảm chấn bán tích cực thông minh (Smart Semi-active Dampers) lắp đặt trên cầu và phát triển hệ thống giám sát sức khỏe kết cấu tự động (Structural Health Monitoring) phục vụ quản lý tài sản cầu đường thông minh.

Kết luận

  1. Thiết lập thành công mô hình toán học và chương trình tính toán số mô phỏng chính xác bài toán tương tác động lực học giữa ô tô nhiều trục và kết cấu nhịp cầu dầm có xét đến độ mấp mô ngẫu nhiên mặt đường.
  2. Xây dựng hoàn chỉnh thuật toán phân tích phổ mật độ công suất mấp mô mặt đường theo tiêu chuẩn ISO 8608, đưa bài toán dao động VBI tại Việt Nam tiệm cận với trình độ cơ học tính toán hiện đại trên thế giới.
  3. Cung cấp bộ dữ liệu thực nghiệm chuẩn xác về hệ số giảm chấn dao động $(\xi)$ của các kết cấu nhịp cầu thép, cầu bê tông cốt thép thực tế đang khai thác tại Việt Nam.
  4. Chứng minh định lượng rằng chất lượng mặt cầu là nhân tố quyết định đến mức độ gia tăng hệ số động lực $(1+IM)$, đưa ra cảnh báo kỹ thuật quan trọng cho công tác thiết kế và duy tu hạ tầng giao thông.
  5. Khẳng định độ tin cậy của mô hình số thông qua thử nghiệm động lực học hiện trường tại cầu Đa Phước với sai số chuyển vị dưới 8%.
  6. Đặt nền móng vững chắc cho các nghiên cứu chuyên sâu tiếp theo về dao động không gian 3D, độ bền mỏi và tối ưu hóa an toàn khai thác cầu đường bộ cao tốc tại Việt Nam.