Ma sát nội tại trong khoa học vật liệu khác gì so với góc ma sát trong địa kỹ thuật?

Trong khoa học vật liệu, ma sát nội tại là thước đo mức độ tiêu tán cơ năng thành nhiệt do tính chất nhớt đàn hồi vi mô khi vật liệu chịu biến dạng đàn hồi chu kỳ nhỏ không phá hủy. Ngược lại, trong địa kỹ thuật, góc ma sát trong là thông số đo sức kháng trượt tiếp tuyến giữa các hạt khoáng vật chống lại sự phá hoại trượt dẻo của đất đá.

Hiệu ứng Snoek liên quan thế nào đến ma sát nội tại?

Hiệu ứng Snoek là cơ chế ma sát nội tại điển hình trong kim loại lập phương tâm khối như sắt, sinh ra do sự khuếch tán nhảy bậc có hướng của các nguyên tử xen kẽ tạp chất (carbon, nitơ) giữa các vị trí bát diện dưới tác dụng của trường ứng suất dao động.

Tại sao góc ma sát trong của đất cát lại phụ thuộc vào độ chặt?

Khi đất cát ở trạng thái chặt chịu biến dạng trượt cắt, các hạt khoáng vật bị chèn chặt phải trèo lên nhau làm thể tích khối đất giãn nở (dilatancy). Công cơ học cần thiết để thắng xu hướng giãn nở này làm tăng góc ma sát trong đỉnh cao hơn so với trạng thái tới hạn.

Những phương pháp thí nghiệm nào dùng để đo ma sát nội tại?

Trong vật liệu rắn, ma sát nội tại được đo bằng con lắc xoắn Kê, phân tích cơ học động lực (DMA) hoặc cộng hưởng siêu âm (RUS). Trong địa kỹ thuật, góc ma sát trong được xác định qua thí nghiệm cắt phẳng trực tiếp theo ASTM D3080/D3080M-11 hoặc thí nghiệm nén ba trục.

Ma sát nội tại là gì? Cơ chế và ứng dụng kỹ thuật

Ma sát nội tại là hiện tượng tiêu tán năng lượng cơ học thành nhiệt lượng hoặc lực cản chống lại sự dịch chuyển tương đối giữa các phần tử vi mô bên trong cùng một môi trường vật chất khi chịu tác động của ngoại lực hoặc biến dạng cơ học. Thuật ngữ này giữ vai trò nền tảng trong hai nhánh khoa học kỹ thuật lớn: trong vật lý chất rắn và khoa học vật liệu, ma sát nội tại đại diện cho tính chất không đàn hồi hoàn toàn (anelasticity) và khả năng hấp thụ rung động của vật liệu; trong khi ở cơ học đất đá và kỹ thuật địa kỹ thuật, ma sát nội tại biểu thị sức kháng cắt phát sinh do sự tiếp xúc trượt và cài răng lược giữa các hạt khoáng vật cấu thành đất đá.

Bản chất vật lý và cơ chế trong khoa học vật liệu

Trong cơ học vật rắn và luyện kim vật lý, ma sát nội tại (tiếng Anh: internal friction, thường được ký hiệu là $Q^{-1}$ hoặc $\tan \delta$ ) phản ánh mức độ bất thuận nghịch nhiệt động lực học trong quá trình biến dạng đàn hồi giả. Theo Nowick và Berry (1972), khi một vật thể rắn chịu ứng suất xoay chiều chu kỳ vô cùng bé, biến dạng thực tế không lập tức đạt trạng thái cân bằng tức thời mà có độ trễ pha so với ứng suất kích thích. Sự lệch pha góc $\delta$ này dẫn tới việc một phần cơ năng bị chuyển hóa vĩnh viễn thành nội nhiệt sau mỗi chu kỳ dao động.

Độ lớn của ma sát nội tại được định lượng thông qua tỉ số giữa phần năng lượng cơ học tiêu tán $\Delta E$ trong một chu kỳ dao động kín và tổng cơ năng biến dạng cực đại $E$ tích lũy trong mẫu vật:

Q^{-1} = \tan \delta = \frac{\Delta E}{2\pi E}

Trong biểu thức trên, $Q^{-1}$ là nghịch đảo của hệ số phẩm chất dao động, $\delta$ là góc trễ pha giữa ten-sơ ứng suất và ten-sơ biến dạng, $\Delta E$ là năng lượng cơ học bị tổn thất trên một chu kỳ, và $E$ là mật độ năng lượng biến dạng đàn hồi cực đại.

Mô hình thư giãn chuẩn và đỉnh Debye

Xét một quá trình hồi phục anelastic đơn lẻ với thời gian hồi phục đặc trưng $\tau$ , phổ ma sát nội tại phụ thuộc vào tần số góc dao động $\omega$ được mô tả bằng hàm đỉnh Debye:

Q^{-1}(\omega) = \Delta \frac{\omega \tau}{1 + \omega^2 \tau^2}

Trong phương trình này, $\Delta$ biểu diễn cường độ hồi phục (relaxation strength), đại diện cho mức độ suy giảm của mô-đun đàn hồi giữa trạng thái không hồi phục và trạng thái hồi phục hoàn toàn; $\omega$ là tần số góc dao động ( $\omega = 2\pi f$ ); và $\tau$ là thời gian hồi phục anelastic. Giá trị ma sát nội tại đạt cực đại $Q^{-1}_{\max} = \frac{\Delta}{2}$ khi điều kiện cộng hưởng $\omega \tau = 1$ được thỏa mãn.

Thời gian hồi phục của các quá trình kích hoạt nhiệt tuân theo định luật Arrhenius:

\tau = \tau_0 \exp\left(\frac{H}{k_B T}\right)

Trong đó $\tau_0$ là thừa số thời gian trước hàm mũ (chu kỳ dao động nhiệt nguyên tử cơ bản), $H$ là năng lượng hoạt hóa của quá trình vi mô tương ứng, $k_B$ là hằng số Boltzmann, và $T$ là nhiệt độ nhiệt động học tuyệt đối tính bằng Kelvin.

Các cơ chế vi mô chính gây tiêu tán năng lượng trong tinh thể

Các nghiên cứu chuyên sâu được tổng kết bởi Blanter và cộng sự (2007) đã phân lập các nguồn gốc vi mô căn bản sinh ra ma sát nội tại trong mạng tinh thể kim loại và hợp kim:

Hiệu ứng Snoek (Snoek relaxation): Xảy ra trong các kim loại có mạng tinh thể lập phương tâm khối (BCC) như sắt alpha khi có mặt các nguyên tử xen kẽ tạp chất (carbon hoặc nitơ). Dưới tác dụng của ứng suất đơn trục ngoài, các vị trí bát diện bất đối xứng trong mạng bị biến dạng, thúc đẩy các nguyên tử xen kẽ khuếch tán nhảy bậc có hướng giữa các vị trí liền kề để giảm thiểu năng lượng tự do, sinh ra đỉnh tiêu tán cơ năng phụ thuộc mạnh vào hàm lượng dung dịch rắn xen kẽ.
Hiệu ứng Zener (Zener relaxation): Bắt nguồn từ sự tái định hướng có trật tự định hướng cục bộ của các cặp nguyên tử thay thế trong hợp kim dung dịch rắn khi chịu ứng suất tác động theo thời gian.
Tiêu tán do lệch mạng (Dislocation damping): Theo mô hình dây rung của Granato và Lücke (1956), các đoạn lệch mạng bị ghim giữ bởi các khuyết tật điểm sẽ rung động và tiêu hao cơ năng trong trường ứng suất dao động. Khi biên độ biến dạng tăng vượt ngưỡng tới hạn, các đoạn lệch mạng bứt khỏi điểm ghim (hiện tượng unpinning), gây ra sự tiêu tán phi tuyến phụ thuộc mạnh vào biên độ dao động.
Thư giãn ranh giới hạt (Grain boundary relaxation): Vào năm 1947, nhà vật lý Kê Tịnh Khang đã sử dụng thiết bị con lắc xoắn trên mẫu nhôm đa tinh thể có độ tinh khiết 99,991% và phát hiện đỉnh ma sát nội tại xuất hiện ở nhiệt độ khoảng 290 độ C. Hiện tượng này chứng minh rằng ranh giới giữa các hạt tinh thể thể hiện hành vi trượt nhớt ở nhiệt độ cao.
Tiêu tán nhiệt đàn hồi (Thermoelastic damping): Hiện tượng truyền nhiệt bất thuận nghịch giữa các miền chịu nén cục bộ (tăng nhiệt độ) và miền chịu kéo cục bộ (hạ nhiệt độ) trong kết cấu dầm hoặc thanh dao động uốn.

Bản chất cơ học và ứng dụng trong địa kỹ thuật và cơ học đất

Trong cơ học môi trường rời, cơ học đất và cơ học đá, ma sát nội tại được hiểu là sức kháng trượt tiếp tuyến phát sinh tại các mặt tiếp xúc giữa các hạt khoáng vật rắn khi khối đất đá bị biến dạng trượt cắt. Đại lượng đặc trưng cho tính chất này là góc ma sát trong, thường được ký hiệu là $\phi$ (hoặc $\phi'$ khi xét theo ứng suất hiệu quả).

Mô hình phá hoại Mohr-Coulomb và nguyên lý ứng suất hiệu quả

Theo định luật phá hoại Mohr-Coulomb cổ điển, sức kháng cắt giới hạn trên một mặt phẳng trượt bất kỳ trong khối vật liệu được biểu diễn qua phương trình:

\tau_f = c + \sigma_n \tan \phi

Trong phương trình trên, $\tau_f$ là cường độ sức kháng cắt giới hạn tại thời điểm phá hoại (tính bằng kPa hoặc MPa), $c$ là lực dính đơn vị của đất, $\sigma_n$ là ứng suất pháp tổng tác dụng vuông góc lên mặt trượt, và $\phi$ là góc ma sát trong của đất đá.

Khi xét đến sự hiện diện của nước trong các lỗ rỗng của đất, Terzaghi đã thiết lập nguyên lý ứng suất hiệu quả, đưa phương trình sức kháng cắt về dạng ứng suất hữu hiệu:

\tau_f = c' + \sigma'_n \tan \phi' = c' + (\sigma_n - u) \tan \phi'

Trong đó $c'$ là lực dính hiệu quả, $\sigma'_n$ là ứng suất pháp hiệu quả truyền qua bộ khung hạt khoáng vật, $u$ là áp lực nước lỗ rỗng thủy tĩnh và thủy động, và $\phi'$ là góc ma sát trong hiệu quả.

Cơ chế ma sát trượt và hiệu ứng giãn nở hạt

Góc ma sát trong của vật liệu hạt không thuần túy là ma sát tiếp xúc trượt bề mặt (interparticle friction angle) mà còn bao gồm thành phần kháng trượt cơ học do sự cài răng lược và xu hướng giãn nở thể tích khi trượt. Nghiên cứu của Bolton (1986) về độ bền và tính giãn nở của cát đã làm sáng tỏ mối quan hệ định lượng giữa góc ma sát trong đỉnh $\phi'_p$ , góc ma sát trong ở trạng thái tới hạn $\phi'_{cv}$ và chỉ số giãn nở tương đối $I_R$ :

\phi'_p - \phi'_{cv} = 3 I_R \quad \text{(trong điều kiện thí nghiệm nén ba trục)}

\phi'_p - \phi'_{cv} = 5 I_R \quad \text{(trong điều kiện biến dạng phẳng)}

Trong các công thức trên, $\phi'_p$ là góc ma sát trong lớn nhất đo được tại thời điểm mẫu đất đạt cường độ cắt cực đại, $\phi'_{cv}$ là góc ma sát trong khi đất đạt trạng thái biến dạng trượt liên tục ở thể tích không đổi theo lý thuyết trạng thái tới hạn của Roscoe và cộng sự (1958), và $I_R$ là chỉ số giãn nở tương đối phụ thuộc vào độ chặt tương đối và mức ứng suất hữu hiệu trung bình.

So sánh ma sát nội tại giữa hai lĩnh vực chuyên ngành

Bảng dưới đây tổng hợp các điểm tương đồng và khác biệt căn bản của khái niệm ma sát nội tại khi áp dụng trong khoa học vật liệu so với cơ học địa kỹ thuật:

Tiêu chí so sánh	Khoa học vật liệu và Vật lý chất rắn	Địa kỹ thuật và Cơ học đất đá
Bản chất hiện tượng	Sự tiêu tán cơ năng thành nhiệt do tính chất bất thuận nghịch vi mô trong biến dạng đàn hồi chu kỳ	Sức kháng cắt tiếp tuyến giữa các bề mặt tiếp xúc hạt chống lại biến dạng trượt vĩnh viễn
Đại lượng định lượng	Nghịch đảo hệ số phẩm chất $Q^{-1}$ , góc lệch pha $\delta$ , cường độ hồi phục $\Delta$	Góc ma sát trong $\phi$ (độ), hệ số ma sát trong $\tan \phi$ , lực dính $c$
Miền biến dạng khảo sát	Biến dạng đàn hồi vi mô biên độ nhỏ từ $10^{-7}$ đến $10^{-4}$ , không phá hủy mẫu	Biến dạng trượt dẻo lớn dẫn tới trượt cắt bề mặt và phá hoại cấu trúc liên kết
Yếu tố chi phối chính	Mật độ lệch mạng, nguyên tử xen kẽ/thay thế, ranh giới hạt, nhiệt độ, tần số dao động	Hình dạng hạt, cấp phối hạt, độ chặt tương đối, độ ẩm, áp lực nước lỗ rỗng
Mục đích ứng dụng	Thiết kế hợp kim giảm chấn cao hoặc vật liệu cộng hưởng có độ suy hao siêu thấp	Tính toán ổn định mái dốc, sức chịu tải nền móng, áp lực đất lên tường chắn

Phương pháp thí nghiệm và tiêu chuẩn đo lường

Đo lường ma sát nội tại trong vật liệu rắn

Để khảo sát phổ ma sát nội tại và các đỉnh hồi phục nhiệt cơ học, các nhà nghiên cứu sử dụng các kỹ thuật thí nghiệm chuyên dụng:

Con lắc xoắn Kê (Kê torsion pendulum): Thiết bị đo dao động xoắn tự do ở dải tần số thấp từ 0,1 Hz đến 10 Hz tại các khoảng nhiệt độ biến đổi từ nhiệt độ âm sâu đến gần nhiệt độ nóng chảy của vật liệu theo mô tả của Nowick và Berry (1972).
Phân tích cơ học động lực (Dynamic Mechanical Analysis, DMA): Đo lường phản ứng nhớt đàn hồi của polymer và vật liệu composite dưới tải trọng uốn hoặc kéo nén chu kỳ có kiểm soát tần số và nhiệt độ.
Phương pháp cộng hưởng siêu âm (Resonant Ultrasound Spectroscopy, RUS): Xác định ma sát nội tại ở dải tần số cao thông qua đo bề rộng vạch cộng hưởng dao động của mẫu thử.

Thí nghiệm xác định góc ma sát trong địa kỹ thuật

Trong địa kỹ thuật và cơ học đất, góc ma sát trong và sức kháng cắt được xác định thông qua các quy trình thí nghiệm chuẩn hóa quốc tế:

Thí nghiệm cắt phẳng trực tiếp (Direct Shear Test): Theo tiêu chuẩn ASTM D3080/D3080M-11 (2011), thí nghiệm được thực hiện trên mẫu đất đặt trong hộp cắt hai nửa chịu các cấp áp lực nén pháp tuyến khác nhau ở điều kiện cố kết thoát nước để xác định góc ma sát trong và lực dính hiệu quả.
Thí nghiệm nén ba trục (Triaxial Compression Test): Thí nghiệm trên mẫu đất hình trụ bọc màng cao su chịu áp lực buồng đẳng hướng và ứng suất lệch trục tăng dần ở các điều kiện không cố kết - không thoát nước (UU), cố kết - không thoát nước (CU có đo áp lực nước lỗ rỗng) hoặc cố kết - thoát nước (CD).

Ứng dụng kỹ thuật thực tiễn và các giới hạn phân tích

Ứng dụng thực tế

Trong thực tiễn kỹ thuật hiện đại, các giá trị ma sát nội tại đóng vai trò quyết định trong việc tối ưu hóa thiết kế kết cấu và độ bền công trình:

Chế tạo hợp kim giảm chấn cao (High Damping Alloys): Ứng dụng các cơ chế ma sát nội tại do chuyển pha martensite hoặc chuyển động lệch mạng để chế tạo vật liệu hấp thụ chấn động trong cánh tuabin, vỏ động cơ và khung bệ máy công cụ nhằm giảm thiểu tiếng ồn và mỏi cơ học.
Phân tích ổn định công trình địa kỹ thuật: Góc ma sát trong là thông số đầu vào bắt buộc trong các mô hình số phân tích phần tử hữu hạn nhằm tính toán hệ số an toàn cho đập đất đá, mái dốc taluy đường giao thông, hố đào sâu trong đô thị và móng cọc chịu tải trọng ngang.
Kiểm tra không phá hủy và nghiên cứu khuyết tật mạng: Đo phổ ma sát nội tại cho phép phát hiện sự tích tụ hydro gây giòn hydro trong thép cường độ cao, đo nồng độ vết của carbon tự do trong thép dập sâu dùng cho công nghiệp ô tô mà các phương pháp hóa học thông thường khó định lượng chính xác.

Giới hạn và vấn đề khoa học còn mở

Mặc dù các mô hình lý thuyết đã được hoàn thiện qua nhiều thập kỷ, việc nghiên cứu ma sát nội tại vẫn đối mặt với các giới hạn kỹ thuật nhất định:

Mô hình tuyến tính Mohr-Coulomb giả định góc ma sát trong không đổi trong suốt quá trình trượt cắt, tuy nhiên trên thực tế đường bao phá hoại của đất đá ở dải ứng suất pháp cực cao hoặc cực thấp thường có dạng phi tuyến rõ rệt.
Hiện tượng trượt dính (stick-slip) và sự phân hủy nhiệt ma sát cục bộ tại các mặt tiếp xúc vi mô trong các trận động đất kiến tạo quy mô lớn đòi hỏi các mô hình ma sát phụ thuộc vận tốc và trạng thái (rate-and-state friction laws) phức tạp hơn thay cho giá trị góc ma sát trong tĩnh truyền thống.
Trong khoa học vật liệu, sự tương tác phi tuyến phức tạp giữa nhiều cơ chế tiêu tán đồng thời (ví dụ như sự chồng lấn giữa hiệu ứng Snoek-Köster do tương tác lệch mạng - nguyên tử xen kẽ và sự trượt ranh giới hạt) gây nhiều thách thức cho việc phân tách định lượng đóng góp của từng khuyết tật đơn lẻ.