Luận án Tiến sĩ Nguyễn Thị Hồng: Động lực học tấm nứt chịu tải trọng di động

Tổng quan về luận án

Trong kỹ thuật công trình hiện đại, giao thông vận tải và công nghiệp quốc phòng, các kết cấu dạng tấm chịu tác dụng của tải trọng di chuyển đóng vai trò kết cấu chịu lực trọng yếu. Điển hình là các tấm bản mặt cầu, dải vệt chống lầy của lực lượng công binh, kết cấu mặt đường băng sân bay quân sự, mặt boong tàu sân bay và tấm lát mặt đường thi công lắp ghép. Dưới tác động khốc liệt của tải trọng khai thác lặp và môi trường, sự xuất hiện của các vết nứt là hiện tượng suy thoái tất yếu làm giảm sút nghiêm trọng khả năng chịu lực, thay đổi đặc trưng dao động riêng và tiềm ẩn nguy cơ phá hủy đột ngột.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: trong khi bài toán dầm có vết nứt chịu tải trọng di động hoặc tấm nguyên vẹn chịu tải trọng di động đã được giải quyết tương đối phong phú bởi Stokes và Willis (1849), Ahmad et al. (2014) [13], Nguyễn Thái Chung và Lê Ngọc Lý (2018) [3], thì các nghiên cứu về tấm có vết nứt đa phần mới dừng lại ở phân tích tĩnh hoặc dao động riêng dưới tải trọng có điểm đặt cố định (Liew et al., 2003 [37]; Krawczuk, 2001 [47]; Xiang et al., 2002 [51]). Các nghiên cứu tích hợp đồng thời sự tồn tại của vết nứt xuyên thủng có hướng bất kỳ, hiệu ứng trượt cắt Mindlin, và tác động tương tác động lực học đầy đủ của tải trọng di chuyển (bao gồm cả khối lượng quán tính và hệ dao động) vẫn còn rất thiếu vắng và chưa có giải pháp số - thực nghiệm hoàn chỉnh.

Nhằm giải quyết triệt để bài toán phức tạp này, luận án tiến sĩ kỹ thuật chuyên ngành Cơ kỹ thuật (Mã ngành: 9.01 / 9520101) của tác giả Nguyễn Thị Hồng, thực hiện tại Học viện Kỹ thuật Quân sự dưới sự hướng dẫn khoa học của GS.TS. Nguyễn Thái Chung (Hà Nội, 2020), với tiêu đề "Phân tích động lực học tấm có vết nứt chịu tải trọng di động", đã thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học tường minh:

• Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa chính xác sự suy giảm độ cứng cục bộ do vết nứt gây ra trong phần tử tấm Mindlin chịu uốn mà không cần giải thuật chia lưới quá mịn xung quanh đỉnh vết nứt?
• Câu hỏi nghiên cứu 2 (RQ2): Quy luật tương tác động lực học giữa tấm có vết nứt và tải trọng di chuyển biến đổi ra sao khi thay đổi vận tốc, gia tốc, độ cứng hệ thống treo và góc nghiêng vết nứt?
• Câu hỏi nghiên cứu 3 (RQ3): Mức độ tương thích giữa mô phỏng phần tử hữu hạn số và phản ứng đo đạc thực nghiệm trên hệ thống đo động đa kênh hiện đại đạt độ tin cậy như thế nào?
• Giả thuyết 1 (H1): Ma trận độ mềm bổ sung $[C_1]$ suy dẫn từ nguyên lý năng lượng Castigliano và cơ học phá hủy đàn hồi tuyến tính (LEFM) phản ánh chính xác sự suy giảm độ cứng cục bộ của phần tử tấm chữ nhật 4 nút có vết nứt góc bất kỳ.
• Giả thuyết 2 (H2): Hiệu ứng quán tính của khối lượng di động và dao động của hệ phương tiện tạo ra hiện tượng cộng hưởng động lực học cục bộ, làm tăng đột biến biên độ võng và phân bố ứng suất tại lân cận đỉnh vết nứt tại các dải vận tốc tới hạn $v_{cr}$.
• Giả thuyết 3 (H3): Mô hình số phát triển trên thuật toán Newmark-$\beta$ có khả năng dự báo tần số dao động riêng và đáp ứng biến dạng động với sai số dưới 5% so với dữ liệu thực nghiệm LMS SCADAS.

Khung lý thuyết của công trình được định vị vững chắc trên nền tảng: Lý thuyết tấm biến dạng cắt bậc nhất Mindlin-Reissner (FSDT), Cơ học phá hủy đàn hồi tuyến tính (Griffith, Irwin), Nguyên lý biến phân Hamilton, Định lý năng lượng biến dạng bổ sung Castigliano, và Thuật toán tích phân trực tiếp Newmark-$\beta$. Công trình khoa học gồm 146 trang thuyết minh, 14 bảng số liệu chi tiết, 67 đồ thị - hình vẽ minh họa, 79 tài liệu tham khảo trong nước và quốc tế, cùng 24 trang phụ lục chuyên sâu chứa mã nguồn phần mềm CPM_2019. Nghiên cứu thực nghiệm được tiến hành trên mẫu tấm thép chế tạo chính xác, kết hợp hệ thống đo rung - ồn - biến dạng đa kênh LMS hiện đại, tạo nên bước đột phá khoa học có giá trị học thuật và ứng dụng thực tiễn vượt bậc.

Literature Review và Positioning

Lịch sử nghiên cứu kết cấu chịu tải trọng di động khởi nguồn từ các công trình kinh điển của Stokes và Willis (1849) [41] khi nghiên cứu đáp ứng của dầm cầu đường sắt. Trải qua hơn một thế kỷ rưỡi, các học giả quốc tế đã phát triển bài toán từ dầm Euler-Bernoulli sang dầm Timoshenko trên nền đàn hồi (Nguyen Xuan Toan & Van Duc Tran, 2017 [60]), hệ liên hợp dầm - dây - cột (Chen Tsung-Chien, 2003 [26]; Liu et al., 2008 [52]; Nguyễn Thái Chung & Lê Phạm Bình, 2018 [59]), và phân tích dao động dầm có vết nứt bằng phương pháp độ cứng động lực học DSM hoặc Wavelet (Ngô Trọng Đức et al., 2016 [5]; Phí Thị Hằng, 2014 [6]; Bilello & Bergman, 2004 [43]; Caddemi & Morassi, 2013 [16]).

Đối với kết cấu dạng tấm (2D), mức độ phức tạp về mặt toán học và cơ học tăng lên theo cấp số nhân:

1. Dòng nghiên cứu tấm không có vết nứt chịu tải trọng di động: Liew et al. (2001) [17] phân tích dao động của tấm chịu hệ dao động di động một bậc tự do; Ahmad et al. (2014) [13] khảo sát đáp ứng phi tuyến của tấm chữ nhật chịu khối lượng di chuyển; Nguyễn Thái Chung và Lê Ngọc Lý (2018) [3], [7] ứng dụng PTHH phân tích tuyến tính và phi tuyến tấm chữ nhật chịu mô hình xe 4 bậc tự do; Đỗ Kiến Quốc và Khổng Trọng Toàn (2007) [9] khảo sát tấm trên nền đàn hồi.
2. Dòng nghiên cứu tấm có vết nứt chịu tải trọng tĩnh hoặc dao động riêng: Boduroglu & Erdogan (1983) [15] sử dụng lý thuyết Mindlin phân tích uốn tĩnh tấm có vết nứt trung tâm; Liew et al. (2003) [37] dùng phương pháp Ritz khảo sát tần số riêng của tấm có vết nứt theo phương bất kỳ; Xiang et al. (2002) [51] thiết lập quan hệ dao động riêng tấm Mindlin; Krawczuk (2001) [47] và Qian et al. (1991) [64] đề xuất mô hình phần tử hữu hạn chữ nhật có vết nứt thông qua ma trận độ mềm bổ sung; Tinh Quoc Bui (2011) [74] phát triển phương pháp IGA cho tấm composite áp điện có nứt; Phuc Pham Minh et al. (2017) [62] và Vu Hoai Nam et al. (2018) [75] ứng dụng XFEM và Phase-field nghiên cứu ổn định tấm FGM có vết nứt.

Tranh luận học thuật sâu sắc tồn tại giữa hai trường phái mô hình hóa vết nứt: Trường phái thứ nhất ủng hộ phương pháp phần tử hữu hạn mở rộng (XFEM) hoặc mô hình trường pha (Phase-field) với ưu thế biểu diễn trường biến dạng cục bộ chi tiết nhưng đòi hỏi chi phí tính toán cực lớn và cấu trúc lưới phức tạp khi giải bài toán động lực học bước thời gian di động; Trường phái thứ hai lựa chọn phương pháp giải tích tương đương thông qua ma trận độ mềm bổ sung Castigliano (Castigliano-based Flexibility Matrix Method), cho phép tính toán trực tiếp ma trận độ cứng phần tử mà không cần tái phân chia lưới vi mô xung quanh đầu vết nứt.

Luận án của tác giả Nguyễn Thị Hồng định vị chính xác tại điểm giao thoa còn bỏ trống: Kế thừa mô hình độ mềm bổ sung của Krawczuk [47] và lý thuyết uốn tấm biến dạng cắt bậc nhất Mindlin [67], nhưng đột phá mở rộng bài toán từ phân tích dao động riêng tĩnh sang phân tích động lực học phi dừng của tấm chịu tải trọng di động phức tạp (khối lượng di động và hệ dao động di động 1 bậc tự do có xét đến vận tốc biến thiên và gia tốc). So với nghiên cứu của Ahmad et al. (2014) [13] (chỉ xét tấm nguyên vẹn) và Liew et al. (2003) [37] (chỉ xét dao động riêng của tấm nứt), luận án đã thực hiện bước nhảy vọt khi giải quyết trọn vẹn bài toán tương tác động học tiếp xúc tấm nứt - phương tiện di chuyển trong miền thời gian thực, đồng thời kiểm chứng thực nghiệm bằng hệ thống đo động hiện đại.

Đó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 biên giới lý thuyết của cơ học kết cấu công trình thông qua việc tích hợp nhất quán các học thuyết cơ học cổ điển và hiện đại:

• Mở rộng Lý thuyết tấm Mindlin-Reissner (1951): Khắc phục triệt để nhược điểm bỏ qua biến dạng trượt ngang của lý thuyết cổ điển Kirchhoff-Love khi tấm có chiều dày từ trung bình đến dày, chuẩn hóa hệ số hiệu chỉnh cắt $\lambda = 5/6$ vào trường chuyển vị ba chiều: $$u(x,y,z,t) = z \theta_y(x,y,t), \quad v(x,y,z,t) = -z \theta_x(x,y,t), \quad w(x,y,z,t) = w_0(x,y,t)$$
• Phát triển mô hình độ mềm Castigliano cho vết nứt: Dưới tác động của cơ học phá hủy đàn hồi tuyến tính (LEFM - Griffith, Irwin), thế năng biến dạng bổ sung $U_1$ do vết nứt chiều dài $2a_c$ gây ra được biểu diễn qua tích phân hệ số cường độ ứng suất $K_I, K_{II}$. Ma trận độ mềm phần tử tổng quát $[C] = [C_0] + [C_1]$ với các phần tử ma trận xác định chính xác theo đạo hàm bậc hai của thế năng biến dạng đối với hệ 9 lực nút độc lập ${F}$: $$c_{ij}^{(n)} = \frac{\partial^2 U_n}{\partial F_i \partial F_j} \quad (i,j = 1 \dots 9; n=0, 1)$$
• Mô hình lý thuyết tương tác động học (Theoretical Model Propositions):
  • Proposition 1 (P1): Ma trận độ cứng phần tử có vết nứt $[K_c^e]$ suy giảm tỷ lệ thuận với hàm hiệu chỉnh chiều dài vết nứt $f(a_c/W_P)$ và hàm chiều dày $\Phi_j(h)$, nhưng vẫn duy trì tính đối xứng và xác định dương trong không gian biến dạng uốn.
  • Proposition 2 (P2): Khi khối lượng $M$ di chuyển với vận tốc $v$ và gia tốc $a$, lực tương tác tiếp xúc thẳng đứng $R(x,y,t)$ xuất hiện thành phần gia tốc Coriolis và gia tốc quán tính tiếp tuyến: $$R(x,y,t) = -M \left[ g + \ddot{w} + 2v \frac{\partial \dot{w}}{\partial x} + v^2 \frac{\partial^2 w}{\partial x^2} + a \frac{\partial w}{\partial x} \right] \delta(x - x_P(t)) \delta(y - y_P(t))$$
  • Proposition 3 (P3): Hệ dao động di động một bậc tự do (khối lượng $M$, độ cứng lò xo $k_s$, cản nhớt $c_s$) đóng vai trò như một bộ lọc thông thấp (low-pass filter), làm biến đổi pha và làm suy giảm biên độ xung kích truyền vào tấm so với mô hình khối lượng chất điểm thuần túy.

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

Khung phân tích của luận án kết nối hài hòa 4 trụ cột lý thuyết: (1) Biến dạng cắt FSDT Mindlin, (2) Năng lượng Castigliano - LEFM, (3) Biến phân Hamilton $\int_{t_0}^{t_1} (\delta T_e - \delta U_e + \delta W_e) dt = 0$, và (4) Thuật giải động lực học Newmark-$\beta$.

Điểm độc đáo cốt lõi nằm ở giải thuật chuyển đổi lực nút: Thiết lập ma trận biến đổi hệ lực $[B]$ kích thước $12 \times 9$ liên kết giữa 12 lực nút phụ thuộc ${S}{12 \times 1}$ với 9 lực nút độc lập ${F}{9 \times 1}$: $${S} = [B]^T {F} \implies [K_c^e] = [B]^T [C_0 + C_1]^{-1} [B]$$ Giải pháp này cho phép phần tử tứ giác 4 nút đẳng tham số chứa vết nứt định hướng bất kỳ $\theta \in [0^\circ, 90^\circ]$ ghép nối trực tiếp vào ma trận độ cứng tổng thể $[K]$ của toàn hệ mà không làm xáo trộn bậc tự do của các phần tử lân cận.

Điều kiện biên (boundary conditions) được thiết lập chặt chẽ cho mọi trường hợp:

• Biên tựa bản lề (Simply Supported - S): $w = 0, \theta_s = 0, M_n = 0$
• Biên ngàm cứng (Clamped - C): $w = 0, \theta_x = 0, \theta_y = 0$

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ lập trường triết học Thực chứng luận (Positivism) kết hợp Hiện thực phê phán (Critical Realism), xem xét các đại lượng cơ học (chuyển vị, vận tốc, gia tốc, ứng suất, tần số) là các thực thể khách quan có thể lượng hóa, mô hình hóa toán học và kiểm chứng thực nghiệm độc lập.

Thiết kế nghiên cứu đa mức độ (multi-level design) bao gồm:

• Mức độ vi mô (Local Element Level): Xây dựng ma trận độ cứng $[K_c^e]$, ma trận khối lượng nhất quán $[M^e] = \int_{V_e} \rho [N]^T [L]^T [L] [N] dV$, và ma trận cản Rayleigh $[C^e] = \alpha_R [M^e] + \beta_R [K^e]$.
• Mức độ vĩ mô (Global Structure Level): Ghép nối ma trận tổng thể $[M]{\ddot{q}} + [C]{\dot{q}} + [K]{q} = {P(t)}$.
• Mức độ tương tác (Dynamic Interaction Level): Xử lý lực tương tác di động phụ thuộc thời gian và tọa độ thông qua hàm Delta-Dirac $\delta(x - v t)$.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được kiểm soát tính chuẩn xác (rigor) qua 4 bước tam giác hóa (triangulation):

1. Tam giác hóa lý thuyết (Theoretical Triangulation): Đối chiếu phương trình dẫn xuất giữa nguyên lý d'Alembert và nguyên lý thế năng biến phân Hamilton.
2. Tam giác hóa phương pháp (Methodological Triangulation): So sánh nghiệm giải tích cổ điển, mô phỏng số phần tử hữu hạn CPM_2019, và thí nghiệm vật lý.
3. Kiểm tra độ tin cậy và chuẩn hóa phần mềm: Chương trình tính tự phát triển bằng mã nguồn chuyên dụng mang tên CRACKED_PLATE_MOVING_2019 (CPM_2019). Thuật toán được kiểm tra độ tin cậy nghiêm ngặt qua 2 bài toán chuẩn quốc tế:
   • Test case 1: Tấm chữ nhật có vết nứt chính giữa chịu lực phân bố đều tĩnh, đối chiếu với công trình của Boduroglu & Erdogan (1983) [15].
   • Test case 2: Tấm chữ nhật nguyên vẹn chịu hệ dao động di động, đối chiếu với kết quả của Liew et al. (2001) [17]. Kết quả kiểm chứng cho thấy sai số chuyển vị tương đối không vượt quá 1.2%.

Data và phân tích

Không gian khảo sát tham số số học (numerical parameter space) được thực hiện toàn diện:

• Kích thước hình học tấm: Chiều dài $a$, chiều rộng $b = W_P$, chiều dày tấm biến thiên trong 6 cấp độ $h = [0.005, 0.010, 0.015, 0.020, 0.025, 0.030]\text{ m}$.
• Thông số vết nứt: Chiều dài vết nứt $2a_c$, số lượng vết nứt ($n = 1, 2, 3$), góc nghiêng vết nứt $\theta = [0^\circ, 30^\circ, 60^\circ, 90^\circ]$.
• Thông số động lực học tải trọng: Vận tốc di chuyển $v = [6, 8, 10, 12, 14]\text{ m/s}$; gia tốc di chuyển $a$; độ cứng lò xo hệ treo $k_s$.
• Thiết bị và giao thức thực nghiệm: Sử dụng kết cấu khung - tấm thép chữ nhật chuyên dụng; thiết bị gia tải di động dẫn động chính xác; hệ thống đo rung - ồn - biến dạng đa kênh chuyên dụng LMS SCADAS tích hợp phần mềm LMS Test.Lab; cảm biến gia tốc áp điện độ nhạy cao (PCB Piezotronics); đầu đo biến dạng điện trở (strain gauges) dán tại lân cận đỉnh vết nứt (điểm A, B). Tín hiệu được xử lý qua bộ lọc số và biến đổi Fourier nhanh (FFT) để nhận dạng tần số riêng và biên độ dao động.

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

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

Khảo sát tham số sâu rộng và dữ liệu thực nghiệm đã mang lại 5 phát hiện mang tính đột phá:

1. Quy luật khuếch đại ứng suất động tại đỉnh vết nứt: Dưới tác dụng của khối lượng di động, tại điểm A (đỉnh vết nứt), ứng suất pháp $\sigma_{y\max}$ và độ võng cực đại $W_{\max}$ tăng phi tuyến theo chiều dài vết nứt $2a_c$. Khi chiều dài vết nứt $2a_c$ tăng từ 0.1m lên 0.3m, độ võng động tăng thêm 34.6%, trong khi ứng suất cục bộ tại đỉnh vết nứt tăng đột biến hơn 82.4%, khẳng định hiệu ứng tập trung ứng suất động khốc liệt hơn nhiều so với tải trọng tĩnh.
2. Hiện tượng vận tốc tới hạn và cộng hưởng động: Đồ thị đáp ứng động lực học theo vận tốc $v$ từ $6\text{ m/s}$ đến $14\text{ m/s}$ chỉ ra sự tồn tại của các đỉnh vận tốc cộng hưởng. Khi khối lượng di chuyển đạt vận tốc tới hạn $v_{cr} \approx 10 - 12\text{ m/s}$ (tùy thuộc vào điều kiện biên), gia tốc thẳng đứng $\ddot{W}_{\max}$ tại tâm tấm đạt cực đại, gấp 2.4 lần so với trường hợp vận tốc thấp $v = 6\text{ m/s}$.
3. Tác động bất đối xứng của góc nghiêng vết nứt ($\theta$): Kết quả khảo sát với $\theta = [0^\circ, 30^\circ, 60^\circ, 90^\circ]$ phát hiện vết nứt có phương vuông góc với hướng di chuyển của tải trọng ($\theta = 90^\circ$) gây ra sự suy giảm độ cứng uốn nghiêm trọng nhất, làm tăng độ võng động thêm 28.3% so với vết nứt song song ($\theta = 0^\circ$).
4. Cơ chế hấp thụ động năng của hệ dao động: Khi thay thế mô hình khối lượng chất điểm bằng hệ dao động một bậc tự do có độ cứng lò xo $k_s$, phản ứng động của tấm giảm đáng kể. Tuy nhiên, khi tăng độ cứng lò xo $k_s \to \infty$, đáp ứng dao động và ứng suất của tấm tiệm cận hoàn toàn về kết quả của bài toán khối lượng di động, chứng minh tính nhất quán nội tại của mô hình lý thuyết.
5. Sự trôi dịch tần số dao động riêng trong thực nghiệm: Dữ liệu đo phổ FFT thực nghiệm trên hệ thống LMS xác nhận: tần số dao động riêng bậc 1 của tấm giảm dần tương ứng từ tấm không nứt ($f_1 = 48.2\text{ Hz}$) xuống tấm có 1 vết nứt ($f_1 = 45.6\text{ Hz}$) và tấm có 3 vết nứt ($f_1 = 41.8\text{ Hz}$). Sai số tần số giữa mô hình số CPM_2019 và thực nghiệm thực tế chỉ dao động từ 2.8% đến 4.3%.

Implications đa chiều

• Đóng góp cho lý thuyết kết cấu: Bổ sung vào kho tàng cơ học tính toán phương pháp tích hợp ma trận độ mềm Castigliano trong phần tử tấm Mindlin, giải quyết trọn vẹn bài toán động lực học phá hủy đàn hồi dưới tác dụng của tải trọng di chuyển.
• Đổi mới phương pháp luận: Cung cấp giải thuật số không cần tái chia lưới (non-remeshing technique) giúp giảm hơn 70% thời gian tính toán ma trận độ cứng so với các phần mềm thương mại sử dụng lưới vi mô xung quanh vết nứt.
• Ứng dụng kỹ thuật thực tiễn: Cung cấp công cụ định lượng đánh giá khả năng chịu tải còn lại của các tấm mặt cầu thép, mặt đường bê tông sân bay, dải vệt chống lầy khi phát hiện vết nứt cục bộ.
• Khuyến nghị chính sách và quy trình kiểm định: Đề xuất bổ sung hệ số xung kích động lực học xét đến vết nứt trong Tiêu chuẩn thiết kế và kiểm định cầu đường (như TCVN 11823:2017), quy định giới hạn vận tốc lưu thông an toàn cho phương tiện quân sự và dân dụng khi qua các nhịp kết cấu bị suy thoái.

Limitations và Future Research

Mặc dù đạt được những đóng góp học thuật xuất sắc, luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu (boundary conditions):

1. Giả thiết hình học vết nứt: Vết nứt được giả định là vết nứt xuyên thủng hoàn toàn, hở, hai mép vết nứt không tiếp xúc va đập vào nhau và không xét đến sự lan truyền phát triển của vết nứt theo thời gian trong quá trình chịu tải.
2. Điều kiện tiếp xúc: Giả thiết tải trọng luôn bám dính liên tục trên bề mặt tấm, chưa xét đến hiện tượng phương tiện bị nảy tách rời mặt đường do mấp mô mấp mô ngẫu nhiên.
3. Đặc tính vật liệu: Nghiên cứu tập trung trên vật liệu kim loại đồng nhất, đẳng hướng, đàn hồi tuyến tính; chưa mở rộng khảo sát cho vật liệu có cơ tính biến thiên (FGM) hoặc vật liệu composite nhiều lớp.
4. Quy mô thực nghiệm: Thí nghiệm mới thực hiện với mô hình khối lượng di động trượt trên tấm phẳng, chưa thực hiện thực nghiệm hệ dao động nhiều bậc tự do với bề mặt gồ ghề.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng đi giàu tiềm năng:

• Mở rộng thuật toán cho bài toán lan truyền vết nứt mỏi động học (Dynamic Fatigue Crack Propagation) kết hợp tiêu chí cơ học phá hủy phi tuyến.
• Mô phỏng tương tác động lực học xe nhiều trục - mặt cầu mấp mô ngẫu nhiên theo tiêu chuẩn ISO 8608.
• Ứng dụng vật liệu composite áp điện thông minh (PZT) dán gia cường chủ động để triệt tiêu dao động và kiểm soát ứng suất tại đỉnh vết nứt.
• Phát triển thuật toán chẩn đoán ngược vị trí và kích thước vết nứt dựa trên đáp ứng gia tốc đo được thời gian thực kết hợp mạng nơ-ron nhân tạo (AI/Deep Learning).

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

Công trình của NCS. Nguyễn Thị Hồng và GS.TS. Nguyễn Thái Chung tạo ra tầm ảnh hưởng học thuật và thực tiễn sâu rộng:

• Tác động học thuật: Đặt nền móng vững chắc cho các nghiên cứu chuyên sâu về động lực học công trình tại Việt Nam; các bài báo công bố từ luận án trên các tạp chí chuyên ngành uy tín là nguồn tài liệu tham khảo chất lượng cao, ước tính thu hút hàng trăm trích dẫn trong lĩnh vực cơ học tính toán kết cấu.
• Tác động ngành và công nghiệp quốc phòng: Ứng dụng trực tiếp trong việc tính toán, thiết kế và khai thác an toàn các trang thiết bị công binh vượt sông (cầu phao, phà cơ giới), dải vệt chống lầy bánh xích, và đường băng dã chiến phục vụ quốc phòng - an ninh.
• Lợi ích kinh tế - xã hội: Tối ưu hóa chi phí duy tu, bảo dưỡng định kỳ các công trình giao thông trọng điểm; ngăn ngừa sự cố sập đổ kết cấu do mỏi và nứt gãy bất ngờ; kéo dài tuổi thọ khai thác của hạ tầng giao thông quốc gia.

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

• Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp thiết lập thuật toán phần tử hữu hạn nâng cao, nguyên lý biến phân Hamilton và kỹ thuật tích hợp độ mềm Castigliano cho bài toán động lực học di động.
• Các nhà khoa học và chuyên gia Cơ học tính toán: Sử dụng khung phân tích và thuật toán Newmark-$\beta$ cải tiến như một chuẩn đối sánh (benchmark) tin cậy cho các mô hình số phức tạp hơn.
• Kỹ sư R&D và tư vấn thiết kế công trình: Ứng dụng mã nguồn CPM_2019 và các biểu đồ khảo sát tham số để tính toán kiểm tra độ an toàn của bản mặt cầu thép, tấm lát đường sân bay trong điều kiện chịu tải trọng nặng di chuyển.
• Cơ quan quản lý hạ tầng giao thông và Quốc phòng: Có cơ sở khoa học xác đáng để xây dựng quy chuẩn đánh giá mức độ nguy hiểm của vết nứt và ban hành biểu đồ giới hạn tải trọng - vận tốc khai thác an toàn.

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 thành công Lý thuyết tấm biến dạng cắt bậc nhất Mindlin kết hợp Định lý năng lượng biến dạng bổ sung Castigliano để giải bài toán động lực học tiếp xúc phi dừng của tấm có vết nứt góc bất kỳ chịu tải trọng di động. Khác với các nghiên cứu trước đây chỉ áp dụng Castigliano cho bài toán tĩnh hoặc dao động riêng, luận án đã dẫn xuất tường minh ma trận độ cứng phần tử suy giảm $[K_c^e] = [B]^T [C_0 + C_1]^{-1} [B]$ và tích hợp trơn tru vào hệ phương trình vi phân chuyển động phi tuyến trong miền thời gian.

2. Đổi mới phương pháp luận của luận án so với các công trình quốc tế tiêu biểu? So với nghiên cứu của Krawczuk (2001) [47] (chỉ dừng ở dao động riêng) và Ahmad et al. (2014) [13] (chỉ phân tích tấm nguyên vẹn), luận án đã thiết lập một chu trình nghiên cứu khép kín từ lý thuyết - thuật giải PTHH tự lập trình (CPM_2019) đến đo đạc thực nghiệm đa kênh LMS SCADAS. Giải thuật của luận án tránh được việc tái tạo lưới vi mô phức tạp quanh đỉnh vết nứt của phương pháp FEM truyền thống/XFEM, giúp giảm thiểu tài nguyên tính toán nhưng vẫn đảm bảo độ chính xác cao.

3. Phát hiện bất ngờ nhất (counter-intuitive) rút ra từ dữ liệu là gì? Phát hiện bất ngờ nhất là ảnh hưởng phi tuyến tính mạnh mẽ của góc nghiêng vết nứt $\theta$: Vết nứt vuông góc với phương di chuyển ($\theta = 90^\circ$) không chỉ gây ra độ võng lớn nhất mà còn tạo ra sự thay đổi đột ngột về pha dao động của gia tốc thẳng đứng $\ddot{W}$, tạo nên các xung lực va đập quán tính thứ cấp làm trầm trọng thêm trường ứng suất đỉnh vết nứt khi vận tốc tải trọng tiệm cận vùng cộng hưởng.

4. Luận án có cung cấp đầy đủ giao thức để tái lập nghiên cứu (replication protocol) không? Có, luận án cung cấp giao thức tái lập hoàn hảo với độ minh bạch học thuật tuyệt đối: Toàn bộ cấu trúc ma trận chuyển đổi lực nút $[B]_{12 \times 9}$, công thức giải tích tính toán các phần tử độ mềm $[C_0], [C_1]$, hàm dạng phần tử đẳng tham số, sơ đồ khối thuật giải Newmark-$\beta$, mã nguồn chi tiết của chương trình CPM_2019, cùng thông số hình học, vật liệu và sơ đồ bố trí cảm biến thực nghiệm LMS đều được trình bày rõ ràng trong thuyết minh và phụ lục.

5. Lộ trình nghiên cứu 10 năm tiếp theo (10-year research agenda) được vạch ra như thế nào? Lộ trình 10 năm tập trung phát triển: (1) Tích hợp mô hình lan truyền mỏi nứt gãy động lực học (Dynamic Fatigue Propagation) trong miền thời gian; (2) Mở rộng bài toán cho kết cấu tấm vỏ composite FGM chịu tải trọng đoàn xe ngẫu nhiên; (3) Tích hợp công nghệ cảm biến mạng không dây và thuật toán học sâu (Deep Learning) để xây dựng hệ thống giám sát sức khỏe kết cấu thời gian thực (Real-time Structural Health Monitoring - SHM) và bản sao số (Digital Twin) cho các công trình giao thông trọng điểm.

Kết luận

Luận án "Phân tích động lực học tấm có vết nứt chịu tải trọng di động" của tác giả Nguyễn Thị Hồng là một công trình nghiên cứu khoa học xuất sắc, mẫu mực và toàn diện trong chuyên ngành Cơ kỹ thuật. Tổng kết các giá trị cốt lõi của luận án bao gồm:

1. Thiết lập thành công hệ phương trình vi phân chuyển động mô tả đầy đủ tương tác động lực học giữa tấm Mindlin có vết nứt xuyên thủng góc bất kỳ và các mô hình tải trọng di động (khối lượng di động, hệ dao động di động) có vận tốc và gia tốc biến thiên.
2. Xây dựng thuật toán phần tử hữu hạn tối ưu và phần mềm chuyên dụng CPM_2019, cho phép xác định chính xác đáp ứng chuyển vị, gia tốc và phân bố ứng suất động mà không đòi hỏi kỹ thuật làm mịn lưới cục bộ phức tạp.
3. Phát hiện và lượng hóa các quy luật cơ học then chốt: Tác động của chiều dài vết nứt $2a_c$, góc nghiêng vết nứt $\theta$, vận tốc tới hạn $v_{cr}$, gia tốc tải trọng, độ cứng lò xo hệ thống treo, và điều kiện liên kết biên đối với ứng xử động lực học của kết cấu.
4. Thực hiện nghiên cứu thực nghiệm công phu, chuẩn xác trên hệ thống đo đa kênh LMS SCADAS, cung cấp bộ dữ liệu thực nghiệm quý giá kiểm chứng độ tin cậy của mô hình lý thuyết với sai số dưới 5%.
5. Mở ra các hướng ứng dụng thực tiễn to lớn trong thiết kế, thẩm định an toàn, đánh giá khả năng chịu lực còn lại và quy hoạch bảo trì các kết cấu tấm mặt cầu, đường băng sân bay, công trình công binh quốc phòng khi xuất hiện khuyết tật vết nứt.
6. Định hình 3 hướng nghiên cứu mới mang tầm quốc tế: Bản sao số kết cấu nứt thời gian thực, điều khiển dao động chủ động bằng vật liệu thông minh, và cơ học phá hủy mỏi động lực học phi tuyến dưới tải trọng di chuyển ngẫu nhiên.