Tổng quan về luận án

Nghiên cứu động lực học kết cấu tấm chịu tải trọng di động đóng vai trò nền tảng trong kỹ thuật công trình, đặc biệt trong các lĩnh vực giao thông vận tải, xây dựng hạ tầng, kỹ thuật công trình biển và công nghiệp quốc phòng như tấm bản mặt cầu, dải vệt chống lầy của công binh, mặt cầu phao quân sự, tấm ghi sân bay và đường băng cất/hạ cánh trên tàu sân bay. Trong quá trình khai thác lâu dài dưới mật độ và cường độ tải trọng gia tăng, kết cấu tấm khó tránh khỏi sự xuất hiện của các khuyết tật dạng vết nứt cục bộ. Sự tồn tại đồng thời của vết nứt làm suy giảm độ cứng kết cấu và hiệu ứng quán tính phức tạp sinh ra từ tải trọng di động gây ra hiện tượng tập trung ứng suất, khuếch đại dao động và tiềm ẩn nguy cơ phá hủy đột ngột.

Mặc dù lý thuyết tấm chịu tải trọng động đã phát triển mạnh mẽ từ các nghiên cứu sơ khai của Stokes & Willis (1849) trên dầm đến các mở rộng cho tấm của Liew & Hung (1997) hay Ahmad et al. (2013), phần lớn các công trình quốc tế và trong nước trước đây chỉ tập trung vào tấm lành lặn chịu tải di động hoặc tấm có vết nứt chịu tải trọng tĩnh/tải điều hòa tại vị trí cố định. Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại ở việc thiếu vắng một mô hình số hữu hiệu và giải pháp thực nghiệm đồng bộ để mô tả toàn diện tương tác động lực học phi dừng giữa tấm có vết nứt với các mô hình tải trọng di động thực tế (xét đầy đủ quán tính khối lượng di động và dao động riêng của hệ di chuyển).

Luận án tiến sĩ kỹ thuật chuyên ngành Cơ kỹ thuật (Mã ngành: 9.01) của nghiên cứu sinh Nguyễn Thị Hồng, dưới sự hướng dẫn khoa học của GS.TS Nguyễn Thái Chung tại Học viện Kỹ thuật Quân sự (2020), mang tựa đề "Phân tích động lực học tấm có vết nứt chịu tải trọng di động", đã giải quyết trực tiếp khoảng trống học thuật này.

Nghiên cứu tập trung vào 4 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 4 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Làm thế nào để mô hình hóa chính xác sự suy giảm độ cứng cục bộ của phần tử tấm Mindlin chứa vết nứt mà không cần tái phân chia lưới dày đặc quanh chóp vết nứt?
  • H1: Việc bổ sung ma trận độ mềm $[C^1]$ dựa trên cơ học phá hủy đàn hồi tuyến tính (LEFM) và định lý Castigliano vào ma trận độ mềm $[C^0]$ cho phép phản ánh chính xác độ linh động của phần tử nứt.
  • RQ2: Quá trình tương tác động giữa phần tử tấm nứt và tải trọng di động thay đổi ra sao khi xét đến thành phần lực quán tính của tải trọng?
  • H2: Tải trọng khối lượng di động và hệ dao động di động tạo ra các ma trận khối lượng bổ sung $[M^{ep}]$, ma trận cản bổ sung $[C^{ep}]$ và ma trận độ cứng bổ sung $[K^{ep}]$ phụ thuộc thời gian và vị trí trên tấm.
  • RQ3: Mức độ ảnh hưởng của các tham số hình học (chiều dài $2a_c$, góc nghiêng vết nứt $\theta$, chiều dày tấm $h$), vận tốc $v$, gia tốc $a$, và độ cứng lò xo $k_v$ đến đáp ứng động lực học như thế nào?
  • H3: Độ võng, gia tốc và ứng suất cục bộ tại lân cận vết nứt tăng phi tuyến theo chiều dài vết nứt và tốc độ di chuyển, đồng thời xuất hiện các dải vận tốc tới hạn gây hiện tượng cộng hưởng động lực học.
  • RQ4: Phản ứng động lực học đo đạc bằng thực nghiệm có tương thích và bảo chứng cho mô hình toán học giải tích số đã thiết lập hay không?
  • H4: Dữ liệu thực nghiệm thu thập qua hệ thống cảm biến gia tốc và biến dạng đa kênh trên mô hình vật lý đồng dạng xác nhận độ tin cậy của chương trình tính với sai số chấp nhận được trong kỹ thuật.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp chặt chẽ giữa Lý thuyết tấm biến dạng trượt bậc nhất Mindlin (FSDT), Cơ học phá hủy đàn hồi tuyến tính (LEFM) theo hệ số cường độ ứng suất $K_I, K_{II}$, và Nguyên lý tác dụng tối thiểu Hamilton. Đóng góp đột phá của luận án là việc xây dựng thành công thuật toán phần tử hữu hạn (PTHH) cải tiến, phát triển bộ mã nguồn tính toán chuyên dụng CPM_2019 (Cracked Plate Moving 2019), và thực hiện quy trình đo kiểm dao động thực nghiệm đa kênh bằng hệ thống LMS, cung cấp luận cứ khoa học định lượng phục vụ chẩn đoán an toàn và gia cường kết cấu tấm trong kỹ thuật công trình và quân sự.

Literature Review và Positioning

Lịch sử nghiên cứu kết cấu chịu tải trọng di động bắt đầu từ thế kỷ 19 với các công trình kinh điển của Stokes (1849) và Willis (1849) về dầm cầu đường sắt. Trong thế kỷ 20 và đầu thế kỷ 21, trường phái giải tích và phần tử hữu hạn phát triển mạnh với các đóng góp của Fryba (1999) trong cuốn sách chuyên khảo toàn diện về dao động của kết cấu dưới tải trọng di động. Đối với kết cấu dầm, Nguyen Xuan Toan & Van Duc Tran (2015) đã mở rộng bài toán dầm trên nền đàn hồi chịu khối lượng di động có xét đến tính phi tuyến hình học, trong khi Phí Thị Hằng (2016) và Nguyễn Thái Chung, Lê Phạm Bình (2018) đã phân tích đáp ứng tần số và ứng xử phi tuyến của dầm có vết nứt.

Đối với kết cấu tấm, các công bố khoa học quốc tế hình thành hai dòng nghiên cứu chính:

  1. Dòng nghiên cứu tấm không vết nứt chịu tải trọng di động: Wu et al. (1987) và Gbadeyan & Oni (1995) thiết lập các giải pháp giải tích cho tấm mỏng chịu lực di động; Liew & Hung (1997) phân tích dao động của tấm chữ nhật chịu hệ dao động di động; Ahmad et al. (2013) khảo sát đáp ứng động lực học phi tuyến của tấm chịu khối lượng di động với các dải vận tốc khác nhau; Nguyễn Thái Chung & Lê Ngọc Lý (2014, 2015) phân tích tuyến tính và phi tuyến tấm chữ nhật chịu hệ dao động 4 bậc tự do mô phỏng xe cơ giới.
  2. Dòng nghiên cứu kết cấu tấm có vết nứt: Sih et al. (1962) và Rice (1968) đặt nền móng giải tích cho bài toán tấm chịu uốn có vết nứt; Boduroglu & Erdogan (1983) giải bài toán tấm Mindlin có vết nứt cục bộ; Krawczuk et al. (2001) đề xuất mô hình phần tử hữu hạn tấm chữ nhật có vết nứt uốn; Liew et al. (2004) dùng phương pháp không lưới (meshless); Tinh Quoc Bui (2012) áp dụng phương pháp đẳng hình học (IGA) cho tấm composite áp điện có nứt; Pham Minh Phuc et al. (2017) và Vu Hoai Nam et al. (2019) ứng dụng mô hình Phase-field và XFEM khảo sát mất ổn định của tấm nứt.

Tuy nhiên, các tranh biện học thuật lớn tồn tại ở hai trường phái tiếp cận:

  • Trường phái thứ nhất (Mô hình phần tử mở rộng XFEM / Phase-field): Cho phép mô tả chi tiết trường ứng suất kỳ dị tại chóp vết nứt nhưng đòi hỏi chia lưới cực mịn, giải thuật bổ sung hàm làm giàu (enrichment functions) phức tạp, dẫn đến chi phí tính toán (computational cost) tăng vọt khi kết hợp với bài toán tích phân thời gian phi dừng dưới tải trọng di động liên tục thay đổi tọa độ.
  • Trường phái thứ hai (Mô hình độ mềm cục bộ bổ sung tích hợp trong PTHH tiêu chuẩn - Macro-element approach): Được khởi xướng bởi Krawczuk (1993, 2001) và Bose (1998), tiếp cận qua ma trận độ mềm biến dạng đàn hồi Castigliano. Ưu điểm vượt trội là không làm thay đổi bậc tự do tổng thể của hệ, giải quyết trực tiếp trong phần tử chữ nhật hoặc tứ giác đẳng tham số chuẩn.

Luận án của Nguyễn Thị Hồng định vị chính xác tại giao điểm của trường phái độ mềm bổ sung Castigliano và động lực học tải trọng di động phức hợp. Nghiên cứu nâng cấp mô hình tấm mỏng cổ điển Kirchhoff của Krawczuk (2001) lên lý thuyết tấm dày Mindlin có kể đến biến dạng trượt bậc nhất ($\kappa = 5/6$), tích hợp đồng thời hai mô hình tải trọng di động có tương tác quán tính: khối lượng di động (Moving Mass) và hệ dao động một bậc tự do di động (Moving Oscillator) với quỹ đạo và vận tốc bất kỳ (đều hoặc có gia tốc).

So sánh với các nghiên cứu quốc tế tiêu biểu:

  • So với mô hình của Krawczuk (2001) (chỉ phân tích dao động tự do của tấm có nứt dưới dạng bài toán trị riêng tĩnh/modal), luận án đã phát triển bước nhảy vọt sang bài toán động lực học cưỡng bức phi dừng (non-stationary forced vibration) với tải trọng di động có xét lực quán tính tiếp xúc.
  • So với nghiên cứu của Ahmad et al. (2013) (chỉ xét tấm nguyên không nứt chịu khối lượng di động), luận án đã giải quyết thành công bài toán tổng quát cho tấm có nhiều vết nứt ngẫu nhiên về vị trí và phương hướng.

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

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

Luận án đóng góp trực tiếp vào việc mở rộng và làm sâu sắc ba lý thuyết nền tảng của cơ học vật rắn biến dạng:

  1. Mở rộng Lý thuyết tấm Mindlin-Reissner (First-order Shear Deformation Theory - FSDT): Thiết lập hệ phương trình chuyển vị tổng quát: $$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)$$ kết hợp hệ số hiệu chỉnh cắt $\lambda = 5/6$, cho phép mô tả chính xác trường ứng xử cơ học của tấm từ mỏng đến dày trung bình mà không gặp hiện tượng khóa cắt (shear locking) khi rời rạc hóa bằng phần tử 4 nút.
  2. Phát triển Lý thuyết độ mềm bổ sung theo Cơ học phá hủy đàn hồi tuyến tính (LEFM) và Định lý Castigliano: Xác định các thành phần của ma trận độ mềm không nứt $[C^0]$ và độ mềm bổ sung $[C^1]$ thông qua thế năng biến dạng đàn hồi $U_n$: $$c_{ij}^n = \frac{\partial^2 U_n}{\partial F_i \partial F_j} \quad (i,j = 1 \dots 9; n = 0, 1)$$ Trong đó, thế năng $U_1$ do vết nứt gây ra được biểu diễn qua tích phân hàm mật độ giải phóng năng lượng biến dạng liên hệ trực tiếp với các hệ số cường độ ứng suất uốn $K_{Irf}$ và cắt $K_{IIrf}$ theo Mindlin, có hiệu chỉnh kích thước hữu hạn: $$K_{Irf} = \Phi_1 K_{Ir} f(a_c/W_P), \quad K_{IIrf} = \Phi_2 K_{IIr} f(a_c/W_P)$$
  3. Mô hình hóa tương tác động lực học liên hợp Tấm - Tải trọng di động: Biến đổi lực tiếp xúc $R(x,y,t)$ thành tải trọng tương đương phân bố $p(x,y,t)$ qua hàm Delta-Dirac hai chiều $\delta(x - x_P(t), y - y_P(t))$, từ đó tường minh hóa các thành phần lực quán tính của tải trọng di chuyển thành ma trận khối lượng bổ sung $[M_e^p]$, ma trận cản Coriolis/tương hỗ $[C_e^p]$ và ma trận độ cứng bổ sung do hiệu ứng độ cong $[K_e^p]$.

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

Khung phân tích của luận án kết hợp ma trận biến đổi hệ lực nút $[B]$ cấp $12 \times 9$, liên kết 12 thành phần lực nút phụ thuộc ${S}$ với 9 thành phần lực nút độc lập ${F}$ (loại bỏ chuyển động cứng): $${S} = [B]^T {F}$$ Ma trận độ cứng của phần tử tấm có vết nứt $[K_c^e]$ được xác định trực tiếp và tường minh qua công thức đóng: $$[K_c^e] = [B] \left( [C^0] + [C^1] \right)^{-1} [B]^T$$

Điều kiện biên và phạm vi ứng dụng (Boundary conditions): Thuật toán xử lý hoàn hảo các liên kết cứng tuyệt đối (tựa đơn bản lề - SSSS, ngàm - CCCC, tự do - Free, hoặc liên kết hỗn hợp SCSC, CSCS). Vết nứt được giả định là vết nứt xuyên thủng bề dày tấm, hở hoàn toàn (không tiếp xúc khép miệng khi uốn), chiều dài $2a_c$ hữu hạn nằm trong lòng phần tử và chưa phát triển lan truyền trong suốt chu kỳ tác dụng của tải trọng.

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

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ lập trường nhận thức luận thực chứng (Positivism) kết hợp phương pháp luận kiểm chứng thực nghiệm (Empirical Validation). Thiết kế nghiên cứu bao gồm hai trục song hành:

  1. Trục mô hình hóa toán học và tính toán số (Numerical Simulation): Xây dựng thuật toán phần tử hữu hạn trên nền tảng phần tử tứ giác 4 nút đẳng tham số (4-node isoparametric element) kết hợp phần tử chữ nhật chứa vết nứt trung tâm.
  2. Trục thực nghiệm vật lý (Experimental Testing): Chế tạo mô hình thực nghiệm tấm thép liên kết ngàm/tựa trên khung chịu lực chuyên dụng, sử dụng xe mang khối lượng di chuyển có điều khiển vận tốc, kết nối hệ thống thu nhận dữ liệu rung động và biến dạng đa kênh hiện đại.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu số được thực hiện qua các bước chuẩn hóa:

  • Kỹ thuật chia lưới thích nghi cục bộ: Khởi tạo lưới từ phần tử chứa vết nứt, ánh xạ từ hệ tọa độ cục bộ $(x,y)$ sang hệ tọa độ tự nhiên $(\xi, \eta) \in [-1, 1] \times [-1, 1]$ qua các hàm dạng: $$N_i(\xi, \eta) = \frac{1}{4}(1 + \xi_i \xi)(1 + \eta_i \eta) \quad (i = 1 \dots 4)$$
  • Khử ma trận suy biến: Áp dụng phép tích phân số Gauss-Legendre $2 \times 2$ cho uốn và tích phân rút gọn $1 \times 1$ cho thành phần cắt để triệt tiêu hoàn toàn hiệu ứng khóa cắt.
  • Ghép nối ma trận tổng thể: $$[M_{sys}(t)] {\ddot{q}} + [C_{sys}(t)] {\dot{q}} + [K_{sys}(t)] {q} = {F_{sys}(t)}$$ Trong đó, ma trận cản kết cấu $[C]$ được mô hình hóa theo cản nhớt Rayleigh tuyến tính: $$[C] = \alpha_R [M] + \beta_R [K]$$ với hai hệ số $\alpha_R, \beta_R$ xác định qua hai tần số dao động riêng đầu tiên $\omega_1, \omega_2$.
  • Thuật toán tích phân bước thời gian Newmark-Beta: Sử dụng tham số trung bình gia tốc $\gamma = 0.5$, $\beta = 0.25$ đảm bảo tính ổn định vô điều kiện (unconditional stability), bước thời gian $\Delta t$ được lựa chọn thỏa mãn điều kiện hội tụ chính xác $\Delta t \le \frac{1}{20 f_{max}}$.

Data và phân tích

Trong cấu phần thực nghiệm tại Chương 4, hệ thống thiết bị đo kiểm bao gồm:

  • Hệ thống phần cứng thu thập dữ liệu đa kênh LMS SCADAS Mobile kết hợp phần mềm chuyên dụng LMS Test.Lab của hãng Siemens.
  • Cảm biến gia tốc áp điện trở (Piezoelectric Accelerometer) độ nhạy cao $100\text{ mV/g}$, dải tần $0.5 - 5000\text{ Hz}$.
  • Cảm biến biến dạng (Strain Gauge) loại lá điện trở $120.0 \pm 0.5,\Omega$, hệ số biến dạng $k = 2.11 \pm 1%$, dán tại các vị trí tập trung ứng suất quanh mũi vết nứt và tâm tấm.
  • Mẫu thí nghiệm: Tấm thép mác CT3 kích thước $L \times B \times h = 1200\text{ mm} \times 800\text{ mm} \times 4\text{ mm}$, mô đun đàn hồi $E = 2.1 \times 10^{11}\text{ N/m}^2$, khối lượng riêng $\rho = 7850\text{ kg/m}^3$, hệ số Poisson $\nu = 0.3$.
  • Kịch bản thực nghiệm: Khảo sát trên 4 trạng thái tấm (Tấm lành, tấm có 1 vết nứt $2a_c = 40\text{ mm}$, tấm 2 vết nứt, tấm 3 vết nứt) chịu khối lượng di động $m = 5\text{ kg} - 20\text{ kg}$ chuyển động với dải vận tốc thay đổi từ $1.0\text{ m/s}$ đến $4.5\text{ m/s}$.

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

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

Qua hàng ngàn lượt chạy mô phỏng số bằng chương trình CPM_2019 kết hợp đối chiếu dữ liệu đo đạc thực nghiệm, luận án rút ra 5 phát hiện cốt lõi mang tính định lượng cao:

  1. Hiệu ứng phi tuyến của kích thước vết nứt: Khi chiều dài vết nứt $2a_c$ tăng từ $0.1b$ lên $0.5b$, độ võng lớn nhất $W_{max}$ tại tâm tấm tăng phi tuyến từ $15.2%$ đến $38.6%$, trong khi biên độ ứng suất động $\sigma_{max}$ tại chóp vết nứt tăng vọt tới $142.5%$ do sự tập trung ứng suất cục bộ.
  2. Hiện tượng vận tốc cộng hưởng tới hạn ($v_{critical}$): Tồn tại các dải vận tốc mà tại đó hệ số động lực học $D_F = W_{dyn}^{max} / W_{stat}^{max}$ đạt giá trị cực đại ($D_F = 1.45 - 1.72$). Khi vận tốc $v$ vượt qua ngưỡng tới hạn, đáp ứng độ võng đỉnh có xu hướng trễ pha so với vị trí tức thời của tải trọng và đạt cực trị sau khi tải trọng đã đi qua tâm tấm.
  3. Ảnh hưởng quyết định của góc nghiêng vết nứt ($\theta$): Vết nứt có phương vuông góc với quỹ đạo di chuyển của tải trọng ($\theta = 0^\circ$) gây ra mức độ suy giảm độ cứng uốn và khuếch đại ứng suất lớn nhất. Khi góc nghiêng tăng dần lên $\theta = 30^\circ, 60^\circ, 90^\circ$ (song song với quỹ đạo), độ võng cực đại giảm lần lượt $8.4%$, $17.1%$ và $24.3%$.
  4. Sự khác biệt bản chất giữa mô hình Khối lượng di động và Hệ dao động di động: Mô hình khối lượng di động đơn thuần luôn đánh giá vượt mức (overestimate) lực tác dụng đỉnh khi ở vận tốc thấp nhưng lại đánh giá thấp (underestimate) hiện tượng dao động dư (residual vibration). Ngược lại, hệ dao động di động có xét độ cứng lò xo $k_v$ và giảm chấn $c_v$ tạo ra sự biến điệu tần số (frequency modulation), trong đó độ cứng lò xo $k_v$ càng lớn thì phản ứng của hệ càng tiệm cận mô hình khối lượng di động.
  5. Tính nhạy cảm cao của liên kết biên và bề dày tấm: Thay đổi điều kiện biên từ tựa đơn (SSSS) sang ngàm bốn cạnh (CCCC) làm giảm độ võng động lớn nhất tới $64.2%$ và tăng tần số dao động riêng đầu tiên lên $2.25$ lần, giúp triệt tiêu đáng kể ảnh hưởng bất lợi của vết nứt. Tăng bề dày tấm $h$ từ $5\text{ mm}$ lên $15\text{ mm}$ làm giảm độ võng động theo lũy thừa bậc ba của chiều dày.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Hoàn thiện lý thuyết tính toán kết cấu tấm suy giảm độ cứng chịu tác dụng của tải trọng di chuyển; cung cấp công cụ mô phỏng chuẩn xác không cần phần tử đặc biệt (singular elements) phức tạp.
  • Về mặt kỹ thuật và thực tiễn công trình: Cung cấp biểu đồ quan hệ định lượng giữa chiều dài nứt cho phép, vận tốc lưu thông và tải trọng xe; trực tiếp định hình các ngưỡng kiểm soát an toàn cho cầu dàn thép, mặt cầu trực xoa (orthotropic steel deck) và đường băng sân bay.
  • Chính sách và quản lý khai thác hạ tầng: Đề xuất quy trình định kỳ đánh giá sức chịu tải của cầu đường dựa trên phổ đáp ứng dao động thực tế; xây dựng bảng tra cứu giới hạn tốc độ xe quân sự hạng nặng khi hành quân qua cầu phao hoặc dải chống lầy có khuyết tật nứt.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  • Giả thiết vết nứt tĩnh, không lan truyền: Luận án giả định vết nứt giữ nguyên hình học trong suốt quá trình tải trọng di chuyển, chưa tích hợp mô hình phá hủy mỏi động lực học (Dynamic Fatigue Crack Propagation) theo định luật Paris-Erdogan.
  • Giới hạn đàn hồi tuyến tính của vật liệu: Chưa xét đến biến dạng dẻo hoặc phá hủy dẻo tại vùng đầu mũi vết nứt khi chịu tải trọng va đập cực lớn.
  • Mô hình tiếp xúc lý tưởng: Giả thiết tải trọng luôn bám dính bề mặt tấm, chưa xét hiện tượng tải trọng nảy tách rời khỏi mặt tấm (wheel lift-off / loss of contact) khi mặt đường gồ ghề biên độ lớn.
  • Dạng hình học và vật liệu: Nghiên cứu chỉ tập trung vào tấm chữ nhật đẳng hướng đơn lớp, chưa mở rộng cho tấm composite nhiều lớp, tấm vật liệu cơ tính biến thiên (FGM) hoặc tấm cong (vỏ).

Chương trình nghiên cứu tương lai (4 hướng phát triển cụ thể):

  1. Mở rộng bài toán sang mô hình vết nứt lan truyền động lực học dưới tác dụng của đoàn tải trọng xe nhiều trục lặp chu kỳ cao (High-cycle fatigue).
  2. Phát triển thuật toán cho kết cấu tấm FGM và tấm composite nano-reinforced có vết nứt chịu tải trọng nhiệt động lực học kết hợp.
  3. Nghiên cứu bài toán chẩn đoán ngược (Inverse Problem): Nhận dạng vị trí và kích thước vết nứt ngầm trong tấm mặt cầu dựa trên dữ liệu đáp ứng dao động đo từ xe thí nghiệm thông minh (Drive-by bridge inspection).
  4. Khảo sát tấm nứt đặt trên nền đàn hồi phi tuyến Pasternak chịu tải trọng di động với vận tốc siêu tới hạn (Supercritical speed).

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

Luận án tạo ra tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật: Là tài liệu tham khảo chuẩn mực trong đào tạo sau đại học chuyên ngành Cơ kỹ thuật, Cơ học vật rắn và Kỹ thuật xây dựng công trình; mở ra các công bố quốc tế uy tín trên các tạp chí ISI/Scopus thuộc nhóm Q1/Q2 (như Journal of Sound and Vibration, Composite Structures, Thin-Walled Structures).
  • Chuyển giao công nghệ công nghiệp - xây dựng: Ứng dụng phần mềm CPM_2019 vào tính toán kiểm toán bản mặt cầu thép, kết cấu sàn công nghiệp chịu xe nâng, và các tấm nắp hầm chịu tải trọng giao thông nặng.
  • Ứng dụng an ninh - quốc phòng: Đóng góp trực tiếp vào công tác thiết kế, kiểm định độ bền và khai thác an toàn các trang thiết bị kỹ thuật quân binh chủng như: dải chống lầy cơ giới, phao cầu tự hành, đường băng dã chiến phục vụ cất hạ cánh máy bay quân sự.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ học/Xây dựng/Giao thông: Tiếp cận thuật toán PTHH mẫu mực, phương pháp lập ma trận độ mềm Castigliano cho phần tử tấm Mindlin nứt và quy trình xử lý bài toán tải trọng di động tương tác.
  • Kỹ sư R&D và Chuyên gia kết cấu công trình: Sử dụng các kết quả và công thức giải tích số để đánh giá nhanh tuổi thọ dư, độ võng động và ứng suất tập trung của kết cấu tấm suy giảm độ cứng trong thực tế khai thác.
  • Các cơ quan quản lý hạ tầng giao thông và Cục Công binh (Bộ Quốc phòng): Sở hữu căn cứ khoa học chuẩn xác để xây dựng quy chuẩn tải trọng, giới hạn tốc độ và quy trình duy tu, bảo dưỡng định kỳ các công trình mặt đường, mặt cầu có vết nứ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 lý thuyết độc đáo nhất là việc thiết lập thành công ma trận độ cứng đóng $[K_c^e] = [B]([C^0] + [C^1])^{-1}[B]^T$ cho phần tử tấm Mindlin 4 nút chứa vết nứt uốn xuyên thủng theo phương bất kỳ, mở rộng trực tiếp Lý thuyết biến dạng trượt bậc nhất (FSDT) kết hợp Cơ học phá hủy đàn hồi tuyến tính (LEFM). Điểm đột phá là không cần tăng số lượng bậc tự do phần tử và không cần chia lại lưới cục bộ mà vẫn thu được trường ứng xử động lực học chính xác.

2. Sự đổi mới về phương pháp nghiên cứu so với các công bố quốc tế tiền nhiệm?

So với các nghiên cứu của Krawczuk et al. (2001) và Ahmad et al. (2013), luận án đã:

  • Tích hợp đồng thời hiệu ứng biến dạng trượt Mindlin ($\kappa = 5/6$) thay vì giả thiết tấm mỏng Kirchhoff thô sơ.
  • Xây dựng hoàn chỉnh ma trận khối lượng bổ sung $[M_e^p]$, ma trận cản $[C_e^p]$ và ma trận độ cứng $[K_e^p]$ phụ thuộc thời gian mô tả trọn vẹn tương tác quán tính giữa tải trọng di động và kết cấu nứt.
  • Kiểm chứng đa chiều thông qua hệ thống đo rung động thực nghiệm đa kênh LMS Test.Lab trên mô hình tấm vật lý kích thước lớn.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu nghiên cứu là gì?

Phát hiện bất ngờ nhất là hiện tượng trễ đỉnh đáp ứng (peak response lag) và sự suy giảm nghịch lý của ứng suất động khi vận tốc tải trọng vượt ngưỡng tới hạn ($v > v_{critical}$). Tải trọng di chuyển quá nhanh khiến kết cấu chưa kịp biến dạng cực đại tại vị trí tiếp xúc, dẫn đến việc ứng suất đỉnh xuất hiện sau lưng tải trọng và biên độ dao động dư sau khi tải trọng rời khỏi tấm trở nên chiếm ưu thế so với dao động trong thời gian tải trọng chạy trên tấm.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) minh bạch không?

Có. Luận án trình bày chi tiết toàn bộ các biểu thức ma trận $[B]$, $[C^0]$, $[C^1]$, các hàm dạng $N_i(\xi, \eta)$, bảng thông số vật lý của mẫu tấm thép thí nghiệm, sơ đồ khối thuật toán Newmark, và công bố trực tiếp mã nguồn chương trình tính toán CPM_2019 trong phần Phụ lục 2 (trang 115-133), cho phép các nhà nghiên cứu độc lập tái lập $100%$ các kết quả số và thực nghiệm.

5. Khung chương trình nghiên cứu 10 năm được định hình từ luận án như thế nào?

Khung nghiên cứu 10 năm định hình lộ trình phát triển: (1) Năm 1-3: Mở rộng bài toán sang mô hình vết nứt lan truyền mỏi phi tuyến (Fatigue Crack Growth); (2) Năm 4-6: Nghiên cứu kết cấu tấm vỏ composite phân lớp chức năng FGM có vết nứt dưới tải trọng di động siêu âm; (3) Năm 7-10: Phát triển hệ thống giám sát sức khỏe kết cấu trực tuyến (Online SHM) kết hợp trí tuệ nhân tạo (AI/Machine Learning) để định vị và lượng hóa vết nứt mặt cầu tự động từ dữ liệu cảm biến rung động gắn trên xe hành trình.

Kết luận

Luận án tiến sĩ của NCS Nguyễn Thị Hồng là một công trình nghiên cứu khoa học công phu, mẫu mực và hoàn chỉnh, mang lại 6 đóng góp cốt lõi:

  1. Xây dựng thành công hệ phương trình vi phân dao động tổng thể mô tả chính xác tương tác động lực học giữa tấm Mindlin có vết nứt và hai mô hình tải trọng di động (Khối lượng di động và Hệ dao động di động).
  2. Thiết lập công thức tường minh xác định ma trận độ mềm bổ sung $[C^1]$ của phần tử tấm có vết nứt uốn, giải quyết triệt để bài toán mô phỏng khuyết tật mà không làm tăng bậc tự do hệ thống.
  3. Phát triển bộ phần mềm chuyên dụng CPM_2019 với giao thuật số Newmark-Beta ổn định, có khả năng giải quyết đa dạng các bài toán tấm có vết nứt với kích thước, vị trí, góc nghiêng và điều kiện biên tùy ý.
  4. Phát hiện và định lượng hóa các quy luật ảnh hưởng của chiều dài vết nứt, vận tốc tới hạn, góc nghiêng và độ cứng lò xo đến đáp ứng võng, gia tốc và ứng suất động cục bộ.
  5. Thực hiện thành công chương trình thực nghiệm đo động đa kênh bằng hệ thống LMS hiện đại, cung cấp bộ dữ liệu thực nghiệm chuẩn xác bảo chứng cho lý thuyết với độ tin cậy cao.
  6. Đề xuất hệ thống khuyến cáo kỹ thuật có giá trị thực tiễn to lớn trong việc kiểm định an toàn, khai thác và gia cường các kết cấu tấm mặt cầu, đường băng và công trình quốc phòng.