Tổng quan về luận án

Nghiên cứu của NCS. Ing. Hoàng Sỹ Tuấn dưới sự hướng dẫn của PGS. Bohdana Marvalová, CSc. tại Khoa Cơ khí, Đại học Kỹ thuật Liberec (Technical University of Liberec, Cộng hòa Séc, 2010), mang tiêu đề "Elastic and Viscoelastic Behaviour of Composites with Elastomeric Matrix" (Hành vi đàn hồi và đàn nhớt của composite nền elastomer). Công trình giải quyết bài toán cơ học phi tuyến phức tạp trong lĩnh vực cơ học vật liệu composite gia cường sợi mềm dẻo (Fiber-Reinforced Elastomers - FREs) và vật liệu đàn hồi nhạy từ (Magneto-Sensitive Elastomers - MSEs). Trong thực tế công nghiệp ô tô, hàng không, cơ sinh học và kết cấu chịu áp lực như lốp xe, ống dẫn khí và lò xo không khí (air-springs), vật liệu elastomer gia cường chịu trạng thái biến dạng hữu hạn (finite strains) với mức độ phi tuyến hình học và vật liệu cao, đồng thời thể hiện tính bất đẳng hướng (anisotropy) và đặc tính suy giảm ứng suất theo thời gian (viscoelastic dissipation).

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở sự thiếu hụt một mô hình liên tục nhất quán kết hợp đồng thời giữa tính bất đẳng hướng hình học của các lớp sợi gia cường, tính phi tuyến siêu đàn hồi của nền cao su và các cơ chế tiêu tán năng lượng đàn nhớt phi cân bằng trong miền biến dạng lớn. Tác giả đã trích dẫn trực tiếp bài toán cốt lõi: "The viscous behavior of the fiber-reinforced composite materials with rubber-like matrix is modeled in the continuum mechanics framework by the Helmholtz free energy function and the evolution equations of the internal variables."

Công trình đặt ra ba câu hỏi nghiên cứu và giả thuyết khoa học chính:

  1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để phân rã hàm thế năng tự do Helmholtz thành các thành phần thể tích, đẳng hướng của nền và bất đẳng hướng của cốt sợi nhằm mô tả chính xác phản ứng cân bằng và phi cân bằng? Giả thuyết 1 (H1): Sự phân rã thế năng tự do kết hợp với tensor cấu trúc bậc bốn cho phép thiết lập biểu thức giải tích tường minh về ứng suất chính cho các trạng thái kéo đơn trục, trượt thuần túy và kéo hai trục.
  2. Câu hỏi nghiên cứu 2 (RQ2): Tham số vật liệu đàn nhớt của từng pha (nền cao su và cốt sợi) có thể được nhận dạng độc lập thông qua chuỗi thí nghiệm tĩnh/động học và kỹ thuật tương quan ảnh số (DIC) hay không? Giả thuyết 2 (H2): Thí nghiệm trượt thuần túy (pure shear) trên mẫu composite dẹt với các góc nghiêng sợi $\varphi \in {30^\circ, 40^\circ, 50^\circ, 60^\circ}$ đủ điều kiện cô lập ứng suất quá mức (overstresses) của từng pha thông qua mô hình Maxwell tổng quát.
  3. Câu hỏi nghiên cứu 3 (RQ3): Mô hình cấu tạo đàn nhớt phi tuyến có thể dự đoán chính xác ứng xử cấu trúc thực tế (lò xo khí nén) và mở rộng cho bài toán ghép trường từ - cơ (magneto-mechanical coupling) không? Giả thuyết 3 (H3): Mô hình hóa liên tục trong COMSOL Multiphysics có khả năng tái hiện hiện tượng nghịch đảo biến dạng cơ học (inversion phenomenon) và kiểm soát độ cứng động học bằng từ trường ngoài.

Khung lý thuyết nền tảng (Theoretical Framework) tích hợp nhiệt động lực học môi trường liên tục (Continuum Thermodynamics), lý thuyết siêu đàn hồi bất đẳng hướng (Anisotropic Hyperelasticity của Holzapfel & Gasser, 2001) và lý thuyết từ đàn hồi phi tuyến (Nonlinear Magnetoelasticity của Dorfmann & Ogden, 2003-2005). Nghiên cứu được thực hiện dưới sự tài trợ của Bộ Giáo dục Cộng hòa Séc (Mã hợp đồng MSM 4674788501) với phạm vi thực nghiệm toàn diện trên cao su E-SBR và composite NBR cốt sợi viscose.

Literature Review và Positioning

Tổng quan tài liệu trong luận án tổng hợp sâu sắc tiến trình phát triển của cơ học chất rắn phi tuyến qua nhiều giai đoạn:

[Lý thuyết Mạng Polyme & Siêu đàn hồi]
  ├─ Elías-Zúñiga (2006) / Drozdov & Dorfmann (2001): Mô hình chuỗi Non-Gaussian
  ├─ Neo-Hooke (1940s) / Mooney-Rivlin (1950s) / Ogden (1972): Nền tảng đẳng hướng
  └─ Hencky / Poživilová (2002) / Plešek & Kruisová (2006): Logarithmic strain
           │
           ▼
[Cơ học Môi trường liên tục Đàn nhớt Biến dạng lớn]
  ├─ Simo (1987) / Holzapfel & Simo (1996): Quá ứng suất tuyến tính & phân rã tích
  ├─ Reese & Govindjee (1998a,b) / Le Tallec et al. (1993): Nhiệt - đàn nhớt phi tuyến
  └─ Nguyen et al. (2007) / Diani et al. (2006): Biến biến dạng nội vi mô
           │
           ▼
[Vật liệu Composite Anisotropic & Magneto-Sensitive (FREs / MSEs)]
  ├─ Holzapfel, Gasser & Ogden (2000, 2001): Tensor cấu trúc & bất biến I4, I6
  ├─ Merodio & Ogden (2005) / deBotton et al. (2006): Đồng hóa vi mô
  └─ Brigadnov & Dorfmann (2003) / Dorfmann & Ogden (2003-2005): Ghép trường Từ - Cơ

Trong dòng nghiên cứu siêu đàn hồi đẳng hướng, các mô hình thống kê mạng chuỗi non-Gaussian (Elías-Zúñiga, 2006; Drozdov & Dorfmann, 2001) và các mô hình hiện tượng học kinh điển như Neo-Hookean, Mooney-Rivlin, Ogden (Holzapfel, 2000) đã đạt độ chín muồi. Các mô hình dựa trên tensor biến dạng logarit của Hencky (Poživilová, 2002; Plešek & Kruisová, 2006) tuy mô tả rất tốt thực nghiệm nhưng gặp trở ngại tính toán lớn khi tích phân số trên miền biến dạng hữu hạn.

Đối với hành vi đàn nhớt hữu hạn, tồn tại hai trường phái đối lập rõ rệt trong y văn:

  • Trường phái tiếp cận bằng biến ứng suất (Stress-like approach): Dẫn dắt bởi Simo (1987), Holzapfel & Simo (1996), Holzapfel & Gasser (2001), sử dụng các biến nội trạng thái dạng quá ứng suất (overstresses $\mathbf{Q}_r$). Ưu điểm nổi bật là dẫn đến các phương trình tiến hóa vi phân tuyến tính giải bằng chuỗi tích phân chập (convolution integral) và thuật toán cập nhật hồi quy (recurrence update formula), cực kỳ hiệu quả về mặt tính toán phần tử hữu hạn.
  • Trường phái tiếp cận bằng biến biến dạng (Strain-like approach): Dẫn dắt bởi Le Tallec et al. (1993), Reese & Govindjee (1998a), Nguyen et al. (2007), dựa trên sự phân rã nhân của gradient biến dạng $\mathbf{F} = \mathbf{F}_e \mathbf{F}_v$. Cách tiếp cận này chặt chẽ về mặt nhiệt động lực học khi hệ rời xa trạng thái cân bằng hoặc chịu tốc độ biến dạng cực cao, nhưng hệ phương trình tiến hóa phi tuyến đa biến vi mô ($\mathbf{C}_v^M, \mathbf{C}_v^F$) làm tăng khối lượng tính toán theo cấp số nhân.

Về mặt định vị nghiên cứu (Positioning), luận án của Hoàng Sỹ Tuấn lấp đầy khoảng trống thực nghiệm và tính toán khi so sánh trực tiếp hai trường phái này trên composite hai họ sợi, đồng thời so sánh với các nghiên cứu quốc tế:

  • So với nghiên cứu của Bischoff et al. (2000, 2002) vốn chỉ dựa trên ô cơ sở trực giao 8 sợi và thử nghiệm nén đơn trục, mô hình của luận án thiết lập biểu thức tường minh cho các chế độ biến dạng phẳng đa trục phức tạp.
  • So với phân tích của Holzapfel & Ogden (2008) về hạn chế của thử nghiệm kéo hai trục phẳng (planar biaxial test) đối với mô mềm sinh học, luận án đã chứng minh bằng thực nghiệm rằng thí nghiệm trượt thuần túy (pure shear) với tỷ lệ mẫu $b/l \approx 7.5$ có độ tin cậy vượt trội trong việc nhận dạng tham số đàn nhớt độc lập cho từng pha composite.

Đó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 và hoàn thiện lý thuyết siêu đàn nhớt bất đẳng hướng của Holzapfel & Gasser (2001) và lý thuyết từ đàn hồi phi tuyến của Dorfmann & Ogden (2004) thông qua các đóng góp cụ thể:

  1. Phân rã liên tục thế năng tự do bất đẳng hướng: Hàm thế năng tự do Helmholtz toàn phần $\Psi$ được phân tách thành thành phần cân bằng đẳng nhiệt $\Psi^{EQ}$ và thành phần phi cân bằng $\Psi^{NEQ}$: $$\Psi(\mathbf{C}, \mathbf{M}_1, \mathbf{M}2, \mathbf{\Gamma}r) = \Psi{vol}^{EQ}(J) + \Psi{iso}^{EQ}(\bar{I}_1, \bar{I}2) + \Psi{ani}^{EQ}(\bar{I}4, \bar{I}6) + \sum{r=1}^N \Psi{r}^{NEQ}(\bar{\mathbf{C}}, \mathbf{\Gamma}_r)$$ Trong đó $\mathbf{M}_1 = \mathbf{a}_0 \otimes \mathbf{a}_0$ và $\mathbf{M}_2 = \mathbf{b}_0 \otimes \mathbf{b}_0$ là các tensor cấu trúc đại diện cho hướng ban đầu của hai họ sợi; $\bar{I}_4 = \mathbf{a}_0 \cdot \bar{\mathbf{C}} \mathbf{a}_0$ và $\bar{I}_6 = \mathbf{b}_0 \cdot \bar{\mathbf{C}} \mathbf{b}_0$ là các giả bất biến (pseudo-invariants) đặc trưng cho độ dãn dài dọc trục sợi.

  2. Mô hình hóa phi cân bằng hai pha độc lập: Mở rộng mô hình Maxwell tổng quát với hai nhánh tiêu tán riêng biệt: nhánh tiêu tán đẳng hướng của nền cao su ($\tau_{iso}, \gamma_{iso}$) và nhánh tiêu tán bất đẳng hướng của cốt sợi viscose ($\tau_{ani}, \gamma_{ani}$), giải quyết sự chênh lệch thời gian hồi phục giữa ma trận đàn hồi và cốt sợi dệt.

  3. Mở rộng lý thuyết tương tác Từ - Cơ (Magneto-Elastic Coupling): Xây dựng biểu thức mật độ thế năng tự do toàn phần phụ thuộc đồng thời vào các bất biến cơ học và vector cảm ứng từ $\mathbf{B}$ (hoặc cường độ từ trường $\mathbf{H}$). Chứng minh rằng việc đưa ứng suất từ (Maxwell stress $\boldsymbol{\tau}_M$) vào tensor ứng suất Cauchy toàn phần giúp bảo toàn tính đối xứng của phương trình cân bằng mô-men, loại bỏ sự phức tạp của lực khối từ trường giả định.

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

Khung phân tích của luận án tích hợp ba trụ cột:

  • Cơ học giải tích giải pháp phẳng: Thiết lập hệ phương trình đóng biểu diễn ứng suất Cauchy chính $(\sigma_1, \sigma_2)$ thông qua tỷ số dãn dài chính $(\lambda_1, \lambda_2, \lambda_3)$ với điều kiện không nén được $\lambda_1 \lambda_2 \lambda_3 = 1$.
  • Kỹ thuật xấp xỉ bậc hai phi cân bằng: Chuyển đổi phương trình vi phân tiến hóa của biến nội $\mathbf{Q}_r$ thành thuật toán sai phân ẩn từng bước (backward Euler), đảm bảo tính ổn định tuyệt đối (unconditional stability) trong mô phỏng FEM.
  • Điều kiện biên giới hạn (Boundary conditions): Xác định ranh giới ứng xử vật liệu khi sợi không chịu nén (fiber tension-only behavior), triệt tiêu thành phần $\Psi_{ani}^{EQ}$ khi $\bar{I}_4 < 1$ hoặc $\bar{I}_6 < 1$ nhằm tránh độ cứng phi vật lý.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi trường phái thực chứng nhất quán (positivism/critical realism), kết hợp chặt chẽ giữa giải tích giải pháp giải tích môi trường liên tục, thực nghiệm đo biến dạng quang học 3D độ chính xác cao và mô phỏng số phần tử hữu hạn (FEM).

   ┌───────────────────────────────────────────────────────────────┐
   │                    QUY TRÌNH NGHIÊN CỨU                       │
   └───────────────────────────────┬───────────────────────────────┘
                                   │
       ┌───────────────────────────┴───────────────────────────┐
       ▼                                                       ▼
[Thực nghiệm Đa trục]                                  [Mô hình hóa Liên tục]
 ├─ Tiền biến dạng loại trừ Mullins (8 chu kỳ)          ├─ Phân rã Helmholtz Free Energy
 ├─ TIRAtest 2810 & Biaxial Rig                         ├─ Tensor cấu trúc M1, M2 & Invariants
 └─ Hệ thống 3D-DIC Q-400 (Istra 4D)                    └─ Maxwell nội hóa (Overstress / Inelastic)
       │                                                       │
       └───────────────────────────┬───────────────────────────┘
                                   ▼
              [Nhận dạng Tham số Vật liệu (Matlab)]
               ├─ Tách pha: Siêu đàn hồi cân bằng (Linear LS / lsqnonlin)
               └─ Tách pha: Đàn nhớt phi cân bằng (Prony series fitting)
                                   │
                                   ▼
             [Mô phỏng Phần tử Hữu hạn (COMSOL Multiphysics)]
              ├─ User-defined Material Subroutines (PDE + Solid Mechanics)
              ├─ Mô phỏng tấm đục lỗ, shear phi tuyến, creep, velocity
              ├─ Mô phỏng Lò xo khí nén (Air-spring) & Inversion Phenomenon
              └─ Ghép trường Từ - Cơ phi tuyến trên MSEs

Thiết kế thực nghiệm phân cấp rõ ràng:

  • Cấp độ vật liệu đẳng hướng (E-SBR): Kiểm tra độ nhạy của các thế năng Neo-Hooke, Mooney-Rivlin và Ogden qua 3 cấu hình biến dạng: Kéo đơn trục (Simple Tension), Trượt thuần túy (Pure Shear) và Kéo hai trục đối xứng (Equibiaxial Tension).
  • Cấp độ composite bất đẳng hướng (NBR + Viscose): Đánh giá ảnh hưởng của góc dệt sợi $\varphi \in {30^\circ, 40^\circ, 50^\circ, 60^\circ}$ trong phép thử trượt thuần túy đa bước thư giãn (multi-step relaxation).

Quy trình nghiên cứu rigorous

  1. Chuẩn bị mẫu và loại bỏ hiệu ứng Mullins: Mẫu thử cao su E-SBR và composite NBR cốt sợi viscose được gia công chính xác. Trước khi thu thập dữ liệu đàn nhớt, toàn bộ mẫu trải qua chu kỳ tiền tải kéo/nhả 8 lần liên tục để khử hoàn toàn hiệu ứng Mullins (Mullins effect), sau đó để hồi phục tự do trong trạng thái không ứng suất từ 2 đến 3 giờ ở nhiệt độ phòng chuẩn $20^\circ\text{C}$.
  2. Hệ thống đo lường quang học không tiếp xúc (3D-DIC): Sử dụng hệ thống Q-400 Dantec Dynamics gồm 2 camera kỹ thuật số độ phân giải cao ghi nhận trường dịch chuyển của các điểm sơn trắng ngẫu nhiên (speckle pattern) trên bề mặt mẫu. Phần mềm Istra 4D xử lý tương quan hình ảnh kỹ thuật số, loại bỏ triệt để sai số đo do trượt ngàm kẹp cơ học.
  3. Giao thức tải thư giãn nhiều bước (Multi-step relaxation protocol): Tải trọng điều khiển theo độ dịch chuyển với vận tốc kéo không đổi $v \le 1\text{ mm/s}$, sau đó giữ cố định dịch chuyển trong suốt 900 giây (15 phút) tại mỗi bậc để vật liệu đạt trạng thái cân bằng nhiệt động lực học.

Data và phân tích

Quá trình khớp dữ liệu (curve fitting) và ước lượng tham số được lập trình tự động trong Matlab thông qua toán tử giải hệ tuyến tính "\" và thuật toán tối ưu phi tuyến bình phương tối thiểu lsqnonlin (Appendix D).

Bảng tổng hợp tham số siêu đàn hồi đẳng hướng thu được từ thực nghiệm kéo đơn trục:

| Mô hình siêu đàn hồi | Biểu thức / Tham số nhận dạng | Mô đun trượt $\mu$ (MPa) | Sai số tương đối $\frac{\sum |y_i - \hat{y}_i|}{\sum y_i} \times 100%$ | | :--- | :--- | :--- | :--- | | Neo-Hookean | $c = 0.3678\text{ MPa}$ | $0.7356$ | $5.73%$ | | Mooney-Rivlin | $c_1 = 0.3441\text{ MPa}; c_2 = 0.0492\text{ MPa}$ | $0.6883$ | $6.84%$ | | Ogden (4 tham số) | $\mu_1 = 0.5428\text{ MPa}, \alpha_1 = 2.25; \mu_2 = 0.00056\text{ MPa}, \alpha_2 = 7.99$ | $0.6115$ | $0.89%$ |

Dữ liệu đàn nhớt đẳng hướng qua thí nghiệm kéo đơn trục:

  • Neo-Hookean: Thời gian hồi phục $\tau = 114.2\text{ s}$, hệ số suy giảm $\gamma = 0.278$, phần dư $\sum y_i^2 = 40.3403$.
  • Ogden: Thời gian hồi phục $\tau = 118.9\text{ s}$, hệ số suy giảm $\gamma = 0.293$, phần dư $\sum y_i^2 = 7.4678$ (chính xác hơn gấp 5.4 lần).

Thử nghiệm trượt thuần túy trên mẫu dẹt kích thước $30\text{ mm} \times 2.9\text{ mm} \times 220\text{ mm}$ (tỷ lệ bề rộng trên chiều cao $b/l \approx 7.5$) xác nhận $\lambda_2 \approx 1$. Mô hình Ogden 4 tham số đạt sai số $3.09%$ với $\tau = 153.9\text{ s}$ và $\gamma = 0.003$.

Đối với composite NBR hai họ sợi (kích thước mẫu trượt $30\text{ mm} \times 4.5\text{ mm} \times 220\text{ mm}$), các tham số nhận dạng đồng thời đạt độ tin cậy vượt trội:

  • Mô đun trượt của nền cao su: $\mu = 3.5\text{ MPa}$.
  • Hằng số độ cứng bất đẳng hướng của sợi: $k = 8.73\text{ MPa}$.
  • Thông số đàn nhớt nền (Isotropic branch): $\gamma_{iso} = 0.06$, $\tau_{iso} = 269\text{ s}$.
  • Thông số đàn nhớt cốt sợi (Anisotropic branch): $\gamma_{ani} = 0.11$, $\tau_{ani} = 54\text{ s}$.

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

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

  1. Sự vượt trội của mô hình Ogden bậc cao trong xấp xỉ biến dạng lớn: Mô hình Ogden 4 tham số kiểm soát sai số dưới $1%$ trong kéo đơn trục ($\varepsilon \le 140%$), vượt trội hoàn toàn so với Neo-Hooke ($5.73%$) và Mooney-Rivlin ($6.84%$).
  2. Hiện tượng bất đối xứng thời gian tiêu tán giữa ma trận và sợi: Thời gian hồi phục ứng suất của cốt sợi viscose ($\tau_{ani} = 54\text{ s}$) diễn ra nhanh hơn gấp 5 lần so với nền cao su NBR ($\tau_{iso} = 269\text{ s}$), chứng minh rằng tiêu tán nội trong composite FREs bị chi phối bởi cơ chế tái sắp xếp vi cấu trúc sợi dệt trong giai đoạn đầu trước khi nền cao su kịp giãn dão.
  3. Phát hiện "Hiện tượng nghịch đảo" (Inversion Phenomenon) trong lò xo khí nén: Khi mô phỏng ống lò xo khí nén chịu áp suất tĩnh $p = 0.9\text{ MPa}$, tác giả phát hiện biến dạng dọc trục thay đổi bản chất tùy thuộc vào góc sợi $\varphi$:
    • Tại $\varphi = 30^\circ$ và $35^\circ$: Ống lò xo bị dãn dài dọc trục khi tăng áp suất.
    • Tại $\varphi = 40^\circ$: Xuất hiện hiện tượng nghịch đảo ("Inversion phenomenon happens to a tube reinforced by fibers with $\varphi = 40^\circ$ "), ống lò xo bị co ngắn lại theo phương dọc trục dù áp suất trong tiếp tục tăng cao.
  4. Hạn chế cấu trúc của mẫu thử kéo hai trục (Biaxial Test Artifact): Luận án phát hiện mô đun đàn hồi nhận dạng từ thiết bị kéo hai trục phẳng đạt giá trị bất thường ($\mu \approx 2.05 - 2.11\text{ MPa}$ so với $\mu \approx 0.51 - 0.74\text{ MPa}$ ở các thử nghiệm khác). Tác giả chỉ ra nguyên nhân là do tập trung ứng suất tại góc cánh tay đòn của mẫu chữ thập không có rãnh xẻ (slits), làm sai lệch trường phân bố ứng suất đồng nhất.
  5. Hiệu ứng gia cường mô đun trượt trong chất đàn hồi nhạy từ (MRE/MSE): Mô phỏng trường ghép trong COMSOL chứng minh mật độ từ thông $\mathbf{B}$ làm biến đổi trực tiếp trạng thái ứng suất cắt mà không cần tăng tải trọng cơ học, mô đun trượt tăng phụ thuộc phi tuyến vào góc nghiêng vector từ trường ngoài.

Implications đa chiều

  • Về mặt lý thuyết: Chứng minh tính khả thi của việc tích hợp song song hai nhánh quá ứng suất Maxwell trong mô hình vật liệu bất đẳng hướng biến dạng hữu hạn, thiết lập cầu nối giữa nhiệt động lực học liên tục và cơ học vi mô.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình thực nghiệm Pure Shear kết hợp 3D-DIC như một giải pháp thay thế hoàn hảo cho Biaxial Test trong việc nhận dạng đặc tính cơ học của composite cao su mỏng.
  • Về mặt ứng dụng kỹ thuật: Cung cấp thông số thiết kế chính xác cho ngành chế tạo lò xo không khí giảm chấn đường sắt và ô tô tải nặng, giúp kỹ sư lựa chọn góc dệt sợi tối ưu ($\varphi \approx 35^\circ$) để triệt tiêu biến dạng dọc trục không mong muốn.

Limitations và Future Research

  1. Giả định đẳng nhiệt (Isothermal condition): Nghiên cứu bỏ qua sự gia nhiệt cấu trúc do tiêu tán nội sinh (viscoelastic self-heating) trong quá trình biến dạng động học tần số cao, vốn là yếu tố làm mềm vật liệu cao su trong thực tế làm việc liên tục.
  2. Ranh giới điều kiện biến dạng của sợi: Chưa xét đến hiện tượng uốn nếp vi mô (fiber micro-buckling) và suy giảm liên kết giao diện sợi - nền (interface debonding) khi sợi chịu nén cục bộ.
  3. Giới hạn trường từ trong MSE: Mô hình từ đàn hồi mới dừng lại ở vật liệu nhạy từ thuần nhất đẳng hướng và bất đẳng hướng lý tưởng, chưa mô tả trễ từ (magnetic hysteresis) và hiệu ứng từ giảo phụ thuộc tần số cao.

Chương trình nghiên cứu tương lai (Future Agenda):

  • Mở rộng mô hình sang bài toán nhiệt - cơ - từ liên hợp đầy đủ (Fully coupled Thermo-Magneto-Mechanical framework).
  • Tích hợp tiêu chuẩn phá hủy mỏi đàn nhớt (viscoelastic fatigue) và phá hủy bám dính giao diện cốt sợi.
  • Tối ưu hóa topo góc sợi biến thiên liên tục (spatially-varying fiber orientations) trên các chi tiết vỏ lốp xe thông minh thế hệ mới.

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

  • Ảnh hưởng học thuật: Luận án là tài liệu tham khảo nền tảng tại Đại học Kỹ thuật Liberec và cộng đồng nghiên cứu cơ học tính toán Trung Âu về thuật toán tích hợp UMAT/PDE trong COMSOL Multiphysics cho vật liệu siêu đàn nhớt.
  • Chuyển giao công nghiệp: Kết quả nghiên cứu về góc sợi và hiện tượng nghịch đảo biến dạng cung cấp cơ sở kỹ thuật trực tiếp cho các nhà sản xuất lốp xe và hệ thống treo khí nén công nghiệp (như Continental, Vibracoustic).
  • Ý nghĩa xã hội và an toàn: Nâng cao độ bền mỏi, giảm thiểu rủi ro nổ vỏ bóng khí giảm chấn trong các phương tiện vận tải hành khách khối lượng lớn.

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

  • Nghiên cứu sinh & Nhà khoa học Cơ học chất rắn: Tiếp cận hệ thống phương trình cấu tạo vi phân đầy đủ cho môi trường biến dạng lớn và mã nguồn Matlab tối ưu hóa tham số.
  • Kỹ sư R&D Ô tô & Hàng không: Ứng dụng mô hình phần tử hữu hạn để mô phỏng chính xác hành vi đàn nhớt, độ dão (creep) và độ chùng ứng suất (stress relaxation) của ống mềm cao su chịu áp.
  • Chuyên gia phát triển vật liệu thông minh (Smart Materials): Khai thác khung lý thuyết phần tử hữu hạn ghép từ - cơ để thiết kế các cơ cấu chấp hành và giảm chấn thích ứng từ biến (magnetorheological dampers).

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 mô hình siêu đàn nhớt bất đẳng hướng của Holzapfel & Gasser (2001) sang hệ hai nhánh quá ứng suất độc lập (Prony series bậc hai) kết hợp phân rã tensor cấu trúc cho composite hai họ sợi viscose/NBR. Tác giả đã chứng minh bằng giải tích và số trị rằng năng lượng tiêu tán của pha sợi ($\tau_{ani} = 54\text{ s}$) và pha nền ($\tau_{iso} = 269\text{ s}$) có thang thời gian phân tách rõ rệt, điều mà các mô hình đàn nhớt đẳng hướng kinh điển hoàn toàn bất lực trong việc mô tả.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của Bischoff et al. (2000) (dùng ô cơ sở 8 sợi chỉ dựa trên nén đơn trục) và Holzapfel & Ogden (2008) (chỉ ra sai số mẫu chữ thập hai trục), luận án đã cải tiến vượt bậc bằng cách thiết lập phương pháp thử nghiệm trượt thuần túy (Pure Shear) trên mẫu tấm dẹt tỷ lệ hình học cao ($b/l \approx 7.5$) kết hợp công nghệ 3D-DIC không tiếp xúc. Phương pháp này khóa chặt độ dãn dài phụ $\lambda_2 = 1$, biến bài toán hai chiều thành biến dạng thuần nhất có thể đo đạc chính xác mà không gặp hiện tượng tập trung ứng suất mép ngàm.

3. Phát hiện thực nghiệm / mô phỏng nào bất ngờ nhất và bằng chứng dữ liệu?

Phát hiện bất ngờ nhất là hiện tượng nghịch đảo biến dạng (Inversion phenomenon) của ống lò xo khí nén cao su cốt sợi kép tại góc dệt $\varphi = 40^\circ$. Khác với trực giác thông thường rằng tăng áp suất khí bên trong ($p = 0.9\text{ MPa}$) sẽ làm ống dãn dài theo mọi hướng (như tại $\varphi = 30^\circ, 35^\circ$), tại $\varphi = 40^\circ$, lực căng tiếp tuyến trong các sợi viscose đã tạo ra phản lực nén dọc trục cưỡng bức, khiến ống lò xo co ngắn lại.

                    HÀNH VI BIẾN DẠNG LÒ XO KHÍ THEO GÓC SỢI
                   (Áp suất trong tĩnh định mức p = 0.9 MPa)

   Góc sợi φ = 30° - 35°                           Góc sợi φ = 40°
   ┌───────────────────┐                         ┌───────────────────┐
   │ ▲ Dãn dài dọc trục│                         │ ▼ Co ngắn dọc trục│
   │ ┃ (Axial Stretch) │                         │ ┃ (Inversion)     │
   │                   │                         │                   │
   │   [ Áp suất p ]   │                         │   [ Áp suất p ]   │
   │   ──────► ◄────── │                         │   ──────► ◄────── │
   │                   │                         │                   │
   │ Ống dãn nở tự do  │                         │ Lực căng sợi ép   │
   │ theo phương trục  │                         │ ngược phương trục │
   └───────────────────┘                         └───────────────────┘

4. Giao thức tái lập thực nghiệm (Replication Protocol) có được cung cấp chi tiết không?

Luận án cung cấp chi tiết toàn diện giao thức tái lập trong Chương 4 và Phụ lục C, D:

  • Kích thước mẫu chuẩn hóa cho kéo đơn trục, trượt thuần túy ($30 \times 2.9 \times 220\text{ mm}$ và $30 \times 4.5 \times 220\text{ mm}$).
  • Quy trình tiền tải 8 chu kỳ khử Mullins.
  • Vận tốc kéo định mức $1\text{ mm/s}$, thời gian giữ tải thư giãn $900\text{ s}$ mỗi bước.
  • Thuật toán nhận dạng tham số trong Matlab (lsqnonlin) và phương thức định nghĩa hàm PDE/Weak Form trong COMSOL Multiphysics.

5. Định hướng nghiên cứu 10 năm tiếp theo từ nền tảng luận án?

Nền tảng của luận án mở đường cho 3 hướng nghiên cứu chiến lược:

  1. Phát triển cấu trúc composite đàn hồi tự phục hồi và tự cảm biến tích hợp hạt nano từ tính.
  2. Mô phỏng đa tỉ xích (Multiscale FE-FFT) từ cấu trúc dệt vi mô sợi đến cụm kết cấu khí nén công nghiệp.
  3. Tối ưu hóa điều khiển chủ động độ cứng hệ thống treo phương tiện giao thông thông qua trường từ ngoài biến thiên thời gian thực.

Kết luận

Luận án Tiến sĩ của Ing. Hoàng Sỹ Tuấn đã đóng góp toàn diện và xuất sắc cho phân ngành Cơ học Môi trường liên tục và Vật liệu Composite:

  1. Đóng góp 1 (Mô hình hóa lý thuyết): Thiết lập khung lý thuyết biến dạng hữu hạn hoàn chỉnh phân rã thế năng Helmholtz cho composite siêu đàn nhớt hai họ sợi và vật liệu đàn hồi nhạy từ.
  2. Đóng góp 2 (Phương pháp giải tích): Xây dựng hệ phương trình giải tích tường minh xác định ứng suất Cauchy phẳng trong các phép thử kéo đơn trục, trượt thuần túy và kéo hai trục.
  3. Đóng góp 3 (Đổi mới thực nghiệm): Chuẩn hóa quy trình đo biến dạng quang học 3D-DIC trên máy kéo nén TIRAtest 2810, loại trừ triệt để hiệu ứng Mullins và biến dạng trượt ngàm.
  4. Đóng góp 4 (Nhận dạng tham số chính xác cao): Xác định bộ tham số vật liệu đàn nhớt phân tách pha cho E-SBR và NBR/viscose với sai số xấp xỉ tối ưu ($0.89%$ với mô hình Ogden).
  5. Đóng góp 5 (Mô phỏng số & Phát hiện hiện tượng): Triển khai thành công các module phi tuyến trên COMSOL Multiphysics, phát hiện hiện tượng nghịch đảo biến dạng tại $\varphi = 40^\circ$ trên lò xo khí nén.
  6. Đóng góp 6 (Mở rộng từ - cơ học): Mô tả thành công tương tác phi tuyến giữa tensor ứng suất từ Maxwell và ứng suất cơ học trong chất đàn hồi nhạy từ.

Công trình khẳng định bước tiến vững chắc về năng lực làm chủ phương pháp toán cơ hiện đại kết hợp thực nghiệm tiên tiến, tạo tiền đề khoa học giá trị cho các thế hệ nghiên cứu tiếp theo trong lĩnh vực cơ học vật liệu mềm và kỹ thuật tính toán công trình.