Tại sao lực ma sát trong mô hình Coulomb không phụ thuộc vào diện tích tiếp xúc biểu kiến?

Theo cơ học tiếp xúc vi mô, bề mặt thực tế luôn có độ mấp mô nên diện tích tiếp xúc thực tế chỉ chiếm một phần rất nhỏ diện tích biểu kiến, thường dưới 1%. Khi diện tích tiếp xúc biểu kiến thay đổi dưới cùng một tải trọng nén, diện tích tiếp xúc thực tế tại các đỉnh vi mấp mô không đổi, do đó độ lớn của lực ma sát tổng cộng được giữ nguyên.

Hệ số ma sát tĩnh và hệ số ma sát động khác nhau như thế nào trong mô hình Coulomb?

Hệ số ma sát tĩnh xác định ngưỡng lực cản giới hạn cần vượt qua để khởi động chuyển động trượt, trong khi hệ số ma sát động xác định lực cản tiếp tuyến duy trì trong quá trình trượt đều. Đối với đa số các cặp vật liệu cơ khí, hệ số ma sát tĩnh luôn lớn hơn hệ số ma sát động do cần ngoại lực lớn hơn để cắt đứt các mối liên kết dính vi mô ban đầu ở trạng thái nghỉ.

Mô hình ma sát Coulomb khác biệt thế nào với tiêu chuẩn phá hoại Mohr-Coulomb?

Mô hình Coulomb mô tả lực ma sát tiếp tuyến xuất hiện giữa hai bề mặt vật rắn riêng biệt tiếp xúc nhau. Ngược lại, tiêu chuẩn Mohr-Coulomb là tiêu chuẩn bền cơ học áp dụng cho khối vật liệu liên tục như đất đá hoặc bê tông, xác định ứng suất cắt giới hạn gây phá hoại trượt nội tại dựa trên lực dính kết đơn vị và góc ma sát trong.

Tại sao cần chính quy hóa mô hình ma sát Coulomb trong mô phỏng phần tử hữu hạn?

Mô hình Coulomb cổ điển sử dụng hàm dấu có bước nhảy gián đoạn tại vận tốc trượt bằng không, gây ra hiện tượng cứng số học và dao động giả mạo khiến thuật toán lặp phi tuyến khó hội tụ. Việc chính quy hóa bằng các hàm giải tích trơn như hyperbolic tangent giúp duy trì tính khả vi liên tục xung quanh điểm không, đảm bảo thuật toán tích phân thời gian hội tụ ổn định và chính xác.

Mô hình Coulomb (Coulomb model) là gì?

Mô hình Coulomb (Coulomb model) là mô hình cơ học cổ điển mô tả lực ma sát khô giữa hai bề mặt tiếp xúc rắn, xác lập mối quan hệ tỷ lệ tuyến tính giữa độ lớn của lực ma sát trượt hoặc lực ma sát nghỉ cực đại và lực ép pháp tuyến, đồng thời độc lập với diện tích tiếp xúc biểu kiến và vận tốc trượt tương đối. Được phát triển và hoàn thiện dựa trên các quan sát thực nghiệm của Guillaume Amontons và Charles-Augustin de Coulomb, mô hình này phân tách ma sát thành hai trạng thái tĩnh và động, đóng vai trò nền tảng trong cơ học tiếp xúc, kỹ thuật ma sát học và tính toán cơ cấu máy.

Khái niệm và bản chất của mô hình Coulomb

Trong cơ học và kỹ thuật ma sát học (tribology), mô hình Coulomb đại diện cho phương pháp tiếp cận vĩ mô kinh điển nhằm định lượng sức cản xuất hiện khi hai bề mặt vật rắn trượt hoặc có xu hướng trượt tương đối với nhau. Hiện tượng này thường được gọi là ma sát khô (dry friction), nhằm phân biệt với các trạng thái ma sát có sự hiện diện của màng chất bôi trơn thủy động dày ngăn cách hoàn toàn hai bề mặt.

Bản chất vật lý của mô hình Coulomb bắt nguồn từ tương tác tiếp xúc ở cấp độ vi mô. Trên thực tế, không có bề mặt gia công cơ khí nào đạt độ phẳng tuyệt đối; mọi bề mặt vật rắn đều sở hữu độ nhấp nhô tế vi bao gồm các đỉnh và đáy mấp mô (asperities). Khi hai bề mặt được đặt ép lên nhau dưới tác dụng của một lực pháp tuyến, tiếp xúc cơ học thực chất chỉ diễn ra tại các đỉnh nhấp nhô này. Do đó, người ta phân biệt hai khái niệm diện tích tiếp xúc:

Diện tích tiếp xúc biểu kiến: Ký hiệu là $A_0$ , là diện tích hình học vĩ mô đo được từ kích thước bao ngoài của vùng tiếp xúc giữa hai vật thể.
Diện tích tiếp xúc thực tế: Ký hiệu là $A_r$ , là tổng diện tích của tất cả các mối tiếp xúc vi mô tại các đỉnh mấp mô tế vi thực sự chạm nhau.

Theo công trình nghiên cứu cơ học tiếp xúc của Kenneth Langstreth Johnson (1985), diện tích tiếp xúc thực tế $A_r$ giữa các bề mặt kỹ thuật thông thường chỉ chiếm một tỷ lệ rất nhỏ so với diện tích biểu kiến $A_0$ , thường ở mức dưới 0,01 (tức dưới 1% diện tích hình học vĩ mô). Dưới áp lực cục bộ cực lớn tập trung tại các đỉnh mấp mô, vật liệu tại các điểm tiếp xúc này chịu biến dạng dẻo nghiêm trọng và hình thành các liên kết dính vi mô (micro-junctions). Khi một lực tiếp tuyến được đặt vào để tạo chuyển động trượt, lực cản ma sát sinh ra chính là tổng lực cần thiết để phá hủy hoặc cắt đứt các liên kết dính vi mô này dọc theo mặt tiếp xúc.

Bối cảnh lịch sử và cơ sở vật lý

Sự hình thành của mô hình Coulomb là kết quả của một quá trình tích lũy thực nghiệm kéo dài qua nhiều thế kỷ trong lịch sử cơ học cổ điển:

Nghiên cứu của Guillaume Amontons (1699): Nhà vật lý người Pháp Guillaume Amontons đã tái khám phá các định luật ma sát thực nghiệm và công bố trước Viện Hàn lâm Khoa học Hoàng gia Pháp vào năm 1699. Amontons đưa ra hai kết luận quan trọng: thứ nhất, lực ma sát tỷ lệ thuận với tải trọng nén pháp tuyến; thứ hai, lực ma sát không phụ thuộc vào diện tích tiếp xúc biểu kiến của hai bề mặt.
Khảo luận của Charles-Augustin de Coulomb (1781 và 1785): Kỹ sư quân sự và nhà vật lý Charles-Augustin de Coulomb đã tiến hành chuỗi thực nghiệm quy mô lớn với các cỗ máy thử nghiệm ma sát phức tạp, trình bày kết quả vào năm 1781 và xuất bản toàn văn trong khảo luận kinh điển mang tựa đề Théorie des machines simples vào năm 1785. Coulomb kiểm chứng các quan sát của Amontons, đồng thời đề xuất định luật thứ ba: lực ma sát động hầu như không phụ thuộc vào vận tốc trượt trong phạm vi vận tốc trung bình. Đặc biệt, Coulomb là người đầu tiên phân biệt rạch ròi giữa lực ma sát tĩnh (cần thiết để khởi động chuyển động) và lực ma sát động (cần thiết để duy trì chuyển động trượt).
Tiền đề tiêu chuẩn cơ học đất (1773): Theo các tổng quan lịch sử ma sát học của Blau (2001) và Popov (2010), trước khi công bố khảo luận về ma sát máy, Coulomb vào năm 1773 đã khảo sát cơ học đất đá và áp dụng quy tắc cực trị vào tĩnh học, đặt cơ sở cho tiêu chuẩn phá hoại Mohr-Coulomb sau này kết hợp giữa lực dính kết vật liệu và ma sát hạt.

Theo tổng quan khoa học của Peter J. Blau (2001), các nghiên cứu tribology hiện đại khẳng định hệ số ma sát không phải là một thuộc tính nhiệt động lực học nội tại của vật liệu như khối lượng riêng hay mô đun đàn hồi. Thay vào đó, hệ số ma sát là phản ứng của toàn bộ hệ thống ma sát học (tribosystem response), chịu sự chi phối chặt chẽ bởi độ nhám bề mặt, độ ẩm môi trường, màng oxit bảo vệ, nhiệt độ tiếp xúc và lịch sử tải trọng.

Cơ chế hoạt động và mô hình toán học

Mô hình Coulomb mô tả hành vi cơ học của ma sát khô thông qua hai trạng thái vận động cơ bản: trạng thái tĩnh (chưa có trượt tương đối) và trạng thái động (đang diễn ra trượt tương đối).

1. Ma sát tĩnh và ngưỡng trượt tới hạn

Khi hai vật thể tiếp xúc và chịu tác dụng của ngoại lực tiếp tuyến nhưng chưa xảy ra chuyển động trượt tương đối, lực ma sát tĩnh xuất hiện để cân bằng chính xác với thành phần lực tiếp tuyến gây trượt. Lực ma sát tĩnh không phải là một giá trị cố định mà biến thiên linh hoạt từ không đến một giá trị cực đại giới hạn. Điều kiện cân bằng tĩnh được biểu diễn bằng bất đẳng thức Coulomb:

F_f \le \mu_s N

Trong công thức trên, đại lượng $F_f$ biểu diễn độ lớn của lực ma sát tĩnh, tham số $\mu_s$ là hệ số ma sát tĩnh không thứ nguyên, và đại lượng $N$ là độ lớn của lực ép pháp tuyến vuông góc với mặt tiếp xúc.

Khi lực tiếp tuyến tăng dần đến mức cực hạn mà vật thể sắp sửa trượt, lực ma sát đạt ngưỡng tới hạn gọi là lực ma sát tĩnh cực đại:

F_{f,\max} = \mu_s N

2. Ma sát động trong quá trình trượt

Khi ngoại lực vượt qua ngưỡng ma sát tĩnh cực đại, liên kết dính vi mô bị phá vỡ và vật thể bắt đầu chuyển động trượt tương đối. Theo mô hình Coulomb, độ lớn của lực ma sát động trong quá trình trượt giữ giá trị không đổi và tỷ lệ trực tiếp với lực pháp tuyến:

F_f = \mu_k N

Trong đó, tham số $\mu_k$ là hệ số ma sát động. Thực nghiệm của Popov (2010) và nhiều tác giả xác nhận rằng đối với hầu hết các cặp vật liệu cơ khí, hệ số ma sát tĩnh luôn lớn hơn hệ số ma sát động, tức là $\mu_s > \mu_k$ . Điều này lý giải tại sao lực cần thiết để bắt đầu đẩy một vật nặng luôn lớn hơn lực cần thiết để duy trì vật đó trượt đều.

Về mặt vector một chiều dọc theo phương tiếp tuyến, lực ma sát động luôn hướng ngược chiều với vector vận tốc trượt tương đối, được mô tả thông qua hàm dấu (signum function):

F_f = -\mu_k N \operatorname{sgn}(v)

với $v$ là vận tốc trượt tương đối giữa hai bề mặt và hàm $\operatorname{sgn}(v)$ nhận giá trị dương một khi vận tốc dương và âm một khi vận tốc âm.

3. Góc ma sát và nón ma sát

Xét tổng hợp lực phản lực toàn phần tại bề mặt tiếp xúc, bao gồm lực ép pháp tuyến $N$ và lực ma sát tiếp tuyến $F_f$ . Góc lệch cực đại giữa phương của phản lực toàn phần và phương pháp tuyến của mặt tiếp xúc được định nghĩa là góc ma sát, ký hiệu là $\phi$ . Mối quan hệ lượng giác giữa góc ma sát và hệ số ma sát được xác lập qua công thức:

\tan \phi = \mu

Trong không gian ba chiều, khi vector lực ma sát có thể định hướng theo mọi phương tiếp tuyến trên mặt phẳng tiếp xúc, tập hợp tất cả các đường tác dụng của phản lực toàn phần tới hạn tạo thành một hình nón tròn xoay gọi là nón ma sát (friction cone). Góc mở nửa đỉnh của nón ma sát chính bằng góc ma sát $\phi$ . Nếu đường tác dụng của ngoại lực tổng hợp tác dụng lên vật thể nằm trọn vẹn bên trong nón ma sát, vật thể sẽ duy trì trạng thái đứng yên tuyệt đối bất kể độ lớn của ngoại lực tăng cao bao nhiêu; đây chính là nguyên lý tự hãm (self-locking) ứng dụng rộng rãi trong các cơ cấu ren ốc, nêm và kích cơ khí.

4. Cơ sở cơ học tiếp xúc vi mô Bowden-Tabor và Johnson

Cơ chế vật lý giải thích tính tuyến tính của định luật Coulomb được Frank Philip Bowden, David Tabor và Kenneth Langstreth Johnson (1985) làm sáng tỏ thông qua mô hình biến dạng dẻo vi mấp mô. Khi hai bề mặt tiếp xúc dưới lực pháp tuyến $N$ , các đỉnh mấp mô chịu áp lực cực hạn vượt quá giới hạn đàn hồi và biến dạng dẻo hoàn toàn. Do đó, diện tích tiếp xúc thực tế tỷ lệ thuận với lực pháp tuyến:

A_r = \frac{N}{H}

với $H$ là độ cứng đâm xuyên tế vi (micro-indentation hardness) của vật liệu mềm hơn trong cặp tiếp xúc. Lực ma sát cần thiết để cắt đứt các mối hàn dính vi mô có độ bền cắt danh định $\tau_s$ được tính bằng:

F_f = \tau_s A_r = \tau_s \frac{N}{H} = \left(\frac{\tau_s}{H}\right) N = \mu N

Công thức trên chứng minh rằng hệ số ma sát $\mu = \tau_s / H$ xấp xỉ không đổi trong phạm vi tải trọng mà diện tích thực tế chưa tiệm cận diện tích biểu kiến, giải thích đầy đủ hai định luật thực nghiệm đầu tiên của Amontons và Coulomb từ góc độ vi cơ học.

Các dạng biến thể và mô hình chính quy hóa trong tính toán số

Mặc dù mô hình ma sát Coulomb cổ điển rất trực quan và hiệu quả trong cơ học giải tích tĩnh, mô hình này bộc lộ nhược điểm toán học nghiêm trọng khi áp dụng vào mô phỏng động lực học và tính toán số hiện đại.

1. Thách thức toán học của điểm gián đoạn bước nhảy

Theo phân tích tổng quan của Bernard Armstrong-Hélouvry, Pierre Dupont và Carlos Canudas de Wit (1994), hàm dấu $\operatorname{sgn}(v)$ trong định luật Coulomb chứa một điểm gián đoạn bước nhảy (jump discontinuity) tại vận tốc trượt $v = 0$ . Tại vị trí này, lực ma sát chuyển đột ngột từ giá trị âm sang dương hoặc ngược lại. Khi tích phân phương trình vi phân chuyển động trong các thuật toán phần tử hữu hạn (FEM) hoặc động lực học đa vật thể (MBD), điểm gián đoạn này gây ra:

Hiện tượng cứng số học (numerical stiffness) khiến thuật toán tích phân thời gian phải giảm bước thời gian xuống vô cùng nhỏ.
Dao động giả (numerical chattering) quanh trạng thái vận tốc bằng không, làm cho thuật toán giải phi tuyến lặp Newton-Raphson không thể hội tụ.
Sự không duy nhất của nghiệm tiếp xúc trong các bài toán va chạm và tiếp xúc đa điểm.

2. Kỹ thuật chính quy hóa liên tục của Oden và Martins (1985)

Để khắc phục khuyết tật bước nhảy của hàm dấu, J. Tinsley Oden và J. A. C. Martins (1985) đã đề xuất các mô hình chính quy hóa ma sát (regularized friction laws). Kỹ thuật này thay thế hàm dấu gián đoạn bằng một hàm giải tích trơn liên tục khả vi xung quanh lân cận vận tốc bằng không.

Biến thể chính quy hóa phổ biến nhất sử dụng hàm hyperbolic tangent:

F_f = -\mu N \tanh\left(\frac{v}{\varepsilon_v}\right)

Một biến thể chính quy hóa khác sử dụng hàm cung tang (arctangent):

F_f = -\mu N \frac{2}{\pi} \arctan\left(\frac{v}{\varepsilon_v}\right)

Trong các phương trình trên, tham số $\varepsilon_v$ được gọi là vận tốc dung sai làm mịn (regularization velocity parameter). Khi vận tốc trượt $|v| \gg \varepsilon_v$ , các hàm trơn này nhanh chóng tiệm cận chính xác về giá trị ma sát Coulomb cổ điển. Trong miền vận tốc cực nhỏ $|v| < \varepsilon_v$ , đường đặc tính ma sát trở thành một hàm tuyến tính trơn đi qua gốc tọa độ, đóng vai trò như một phần tử cản nhớt vi mô (viscous damping) dập tắt các dao động số giả mạo.

Trong thực hành mô phỏng kỹ thuật, Oden và Martins (1985) khuyến nghị chọn tham số làm mịn $\varepsilon_v$ ở mức 0,001 m/s hoặc 0,0001 m/s. Việc thiết lập giá trị này đòi hỏi sự cân bằng kỹ thuật: nếu chọn tham số quá nhỏ, bài toán số sẽ tiệm cận lại hiện tượng cứng và khó hội tụ; nếu chọn tham số quá lớn, mô hình sẽ sinh ra biến dạng trượt giả (creep drift) ngay cả khi vật thể đang ở trạng thái nghỉ.

So sánh mô hình Coulomb với các mô hình ma sát và phá hủy khác

Để hiểu rõ phạm vi ứng dụng và tính đặc thù của mô hình Coulomb, cần đối chiếu mô hình này với các mô hình ma sát bề mặt và tiêu chuẩn bền cơ học lân cận. Bảng dưới đây tóm tắt các đặc trưng phân biệt chủ chốt:

Tiêu chí so sánh	Mô hình ma sát Coulomb	Mô hình ma sát Stribeck	Tiêu chuẩn phá hoại Mohr-Coulomb	Mô hình động học Dahl và LuGre
Bản chất hiện tượng	Ma sát khô giữa hai bề mặt vật rắn tiếp xúc trực tiếp	Ma sát chất lỏng và bôi trơn hỗn hợp có màng dầu	Phá hủy trượt cắt bên trong khối vật liệu liên tục hoặc đất đá	Ma sát động học có xét biến dạng đàn hồi vi gai tiếp xúc
Mối quan hệ với vận tốc trượt	Độc lập với độ lớn vận tốc, chỉ phụ thuộc chiều chuyển động	Phụ thuộc phi tuyến mạnh vào vận tốc trượt tương đối	Không phụ thuộc vận tốc, áp dụng cho trạng thái ứng suất tĩnh	Phụ thuộc phi tuyến vào vận tốc và trạng thái biến dạng nội tại
Tham số đặc trưng	Hệ số ma sát tĩnh và hệ số ma sát động	Độ nhớt động học, tải trọng và tốc độ trượt tương đối	Lực dính kết đơn vị và góc ma sát trong của vật liệu	Độ cứng vi gai, hệ số cản nhớt và thông số trượt vi mô
Hiện tượng mô tả được	Cân bằng tĩnh, tự hãm cơ học, trượt đều ở tốc độ trung bình	Chuyển pha bôi trơn ranh giới, bôi trơn hỗn hợp và thủy động	Mặt trượt gãy vỡ của khối đất đá, bê tông và nền móng	Chuyển động dính trượt (stick-slip) và trễ đàn hồi trước khi trượt
Năm công bố nền tảng	Năm 1785 bởi Charles-Augustin de Coulomb	Năm 1902 bởi Richard Stribeck	Khởi nguồn năm 1773 và được Otto Mohr phát triển (Popov, 2010)	Mô hình Dahl và LuGre khảo sát bởi Armstrong-Hélouvry và cs. (1994)

Phân biệt mô hình ma sát Coulomb và tiêu chuẩn Mohr-Coulomb

Trong các tài liệu kỹ thuật, người học thường nhầm lẫn giữa mô hình ma sát Coulomb trên bề mặt tiếp xúc và tiêu chuẩn phá hoại Mohr-Coulomb trong cơ học đất đá và địa kỹ thuật. Mặc dù cả hai đều mang tên Coulomb, đối tượng vật lý của chúng hoàn toàn khác nhau:

Mô hình ma sát Coulomb mô tả lực cản tiếp tuyến giữa hai bề mặt vật rắn riêng biệt tiếp xúc nhau. Ngược lại, tiêu chuẩn phá hoại Mohr-Coulomb mô tả giới hạn ứng suất gây ra sự trượt phá hủy cắt bên trong một môi trường liên tục (như đất, đá hoặc bê tông). Phương trình tiêu chuẩn Mohr-Coulomb có dạng:

\tau = c + \sigma \tan \phi

Trong phương trình này, đại lượng $\tau$ biểu diễn độ bền chống cắt giới hạn của vật liệu, tham số $c$ là lực dính kết đơn vị (cohesion), đại lượng $\sigma$ là ứng suất nén pháp tuyến tác dụng lên mặt trượt, và tham số $\phi$ là góc ma sát trong của khối hạt vật liệu. Khi lực dính kết $c = 0$ (đối với đất cát rời rạc không dính), tiêu chuẩn Mohr-Coulomb quy đổi về mối quan hệ tỷ lệ tuyến tính tương tự như định luật ma sát Coulomb cổ điển.

Ứng dụng kỹ thuật và thực tiễn tại Việt Nam

Mô hình ma sát Coulomb đóng vai trò cột trụ trong tính toán thiết kế cơ khí, kỹ thuật phương tiện giao thông và an toàn kết cấu tại Việt Nam.

1. Thiết kế cơ cấu truyền động và phanh hãm cơ khí

Trong ngành cơ khí chế tạo máy và công nghiệp lắp ráp ô tô, xe máy tại Việt Nam, mô hình Coulomb được sử dụng để tính toán mô-men ma sát truyền động và mô-men hãm của các cơ cấu quan trọng:

Hệ thống phanh đĩa và phanh tang trống: Khả năng hãm dừng phương tiện giao thông phụ thuộc trực tiếp vào lực ma sát sinh ra giữa má phanh và đĩa phanh hoặc tang trống. Theo các thực nghiệm trong sách giáo khoa của Popov (2010), hệ số ma sát khô giữa vật liệu composite ma sát của má phanh và đĩa thép hoặc gang xám thường được khống chế ổn định trong khoảng từ 0,15 đến 0,45 nhằm đảm bảo hiệu quả phanh mà không gây kẹt cứng đột ngột.
Khớp nối trục ma sát và ly hợp: Tính toán số lượng đĩa ma sát và lực ép lò xo cần thiết để truyền toàn bộ công suất từ động cơ sang hộp số mà không bị trượt trong điều kiện làm việc danh định.
Liên kết bu lông cường độ cao chịu lực ma sát: Trong kết cấu thép xây dựng cầu đường và nhà xưởng công nghiệp tại Việt Nam, các mối nối bu lông cường độ cao truyền tải trọng cắt thông qua lực ma sát siết chặt giữa các bản thép tiếp xúc. Tiêu chuẩn thiết kế kết cấu thép yêu cầu kiểm soát chặt chẽ hệ số ma sát bề mặt thép sau khi làm sạch bằng phương pháp phun cát hoặc phun bi.

2. Hệ thống tiêu chuẩn kỹ thuật quốc gia tại Việt Nam

Việc đo đạc, kiểm chuẩn và áp dụng các quy luật ma sát học tại Việt Nam được quy định chặt chẽ thông qua hệ thống Tiêu chuẩn Quốc gia (TCVN) hài hòa với tiêu chuẩn quốc tế ISO:

TCVN 11725-1:2016: Tiêu chuẩn này hoàn toàn tương đương với tiêu chuẩn quốc tế ISO 7148-1:2012, mang tên Ổ trượt - Thử trạng thái ma sát của vật liệu ổ trượt - Phần 1: Thử kim loại ổ trượt. Tiêu chuẩn quy định chi tiết các quy trình thử nghiệm tribology nhằm đo hệ số ma sát và độ mòn của các hợp kim kim loại làm ổ trượt dưới các điều kiện tải trọng và bôi trơn biên xác định.
TCVN 8287-2:2009: Tương đương với tiêu chuẩn ISO 4378-2:2009, mang tên Ổ trượt - Thuật ngữ, định nghĩa, phân loại và ký hiệu - Phần 2: Ma sát và mòn. Tiêu chuẩn thiết lập hệ thống định nghĩa chuẩn mực tại Việt Nam về các khái niệm ma sát khô, ma sát ướt, hệ số ma sát và các dạng mài mòn cơ học.

Các cơ sở nghiên cứu và đào tạo đầu ngành như Trường Đại học Bách khoa Hà Nội, Trường Đại học Bách khoa thuộc Đại học Quốc gia Thành phố Hồ Chí Minh và Viện Cơ học thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam đều ứng dụng mô hình Coulomb làm nền tảng giảng dạy trong các môn học Cơ học lý thuyết, Chi tiết máy và Ma sát học kỹ thuật.

Ưu điểm và giới hạn áp dụng

Mặc dù là mô hình ma sát được ứng dụng rộng rãi nhất trong kỹ thuật, mô hình Coulomb chỉ là một phép gần đúng lý tưởng hóa và có những ranh giới vật lý xác định mà người kỹ sư cần lưu ý.

1. Ưu điểm nổi bật

Đơn giản và trực quan: Mô hình chỉ sử dụng một tham số duy nhất là hệ số ma sát $\mu$ , giúp đơn giản hóa tối đa các phương trình cân bằng tĩnh học và động lực học cơ bản.
Độ chính xác đủ tin cậy trong phạm vi kỹ thuật thông thường: Đối với các bề mặt kim loại khô làm việc ở dải vận tốc trung bình và tải trọng vừa phải, sai số của mô hình Coulomb thường không vượt quá giới hạn cho phép trong tính toán an toàn kỹ thuật.
Khả năng giải tích hóa cao: Cung cấp nghiệm giải tích tường minh cho các bài toán cân bằng vật thể, góc dốc trượt và điều kiện tự hãm cơ khí.

2. Giới hạn vật lý và ranh giới phá vỡ mô hình

Mô hình Coulomb mất tính đúng đắn khi hệ tiếp xúc rơi vào các điều kiện làm việc cực hạn sau:

Tải trọng nén cực cao (vượt giới hạn biến dạng dẻo): Khi lực nén pháp tuyến tăng lên mức rất lớn, diện tích tiếp xúc thực tế $A_r$ tăng dần và tiệm cận diện tích biểu kiến $A_0$ . Lúc này, các đỉnh mấp mô không còn khả năng lún thêm vào nhau và toàn bộ vùng tiếp xúc chuyển sang trạng thái chảy dẻo cắt liên tục. Theo phân tích của Johnson (1985), lực ma sát không còn tăng tuyến tính theo lực pháp tuyến mà bị chặn trên bởi độ bền chảy cắt của vật liệu. Theo tiêu chuẩn chảy dẻo von Mises, ứng suất cắt giới hạn cực đại được xác định bằng:
$\tau_{\max} = \tau_y \approx 0,577 \sigma_y$
với $\tau_y$ là giới hạn chảy cắt và $\sigma_y$ là giới hạn chảy kéo danh định của vật liệu. Khi đạt ngưỡng này, hệ số ma sát biểu kiến suy giảm đáng kể theo tải trọng nén.
Vận tốc trượt cao: Khi vận tốc trượt tăng cao, nhiệt ma sát sinh ra mãnh liệt làm biến đổi vi cấu trúc bề mặt, gây mềm hóa nhiệt vật liệu hoặc tạo ra các lớp màng oxit thứ cấp có tính chất bôi trơn rắn, làm cho lực ma sát phụ thuộc rõ rệt vào vận tốc trượt.
Môi trường chân không siêu cao: Peter J. Blau (2001) chỉ ra rằng trong không khí thông thường, hệ số ma sát của các bề mặt kim loại khô thường nằm trong khoảng từ 0,2 đến 1,0 nhờ có sự hiện diện của các lớp màng oxit và hơi ẩm hấp phụ tự nhiên. Tuy nhiên, trong môi trường chân không siêu cao (ultra-high vacuum), các màng bảo vệ này bị phá hủy hoàn toàn do cọ xát mà không có oxy để tái tạo. Tiếp xúc kim loại sạch trực tiếp dẫn đến hiện tượng hàn dính nguội mãnh liệt, khiến hệ số ma sát có thể tăng vọt vượt quá 2,0 hoặc dẫn đến kẹt cứng hoàn toàn mối tiếp xúc.
Điều kiện có chất bôi trơn: Khi xuất hiện chất lỏng bôi trơn tạo nên áp suất thủy động nâng tách hai bề mặt, hệ số ma sát có thể giảm xuống mức 0,05 hoặc thấp hơn (Armstrong-Hélouvry và cs., 1994). Trong trường hợp này, ma sát chuyển sang tuân theo định luật nhớt Newton và đường cong Stribeck thay vì mô hình Coulomb.
Thang đo vi mô và nanomet: Ở cấp độ nguyên tử và màng phân tử đơn lớp, các lực liên kết bề mặt van der Waals và lực hút tĩnh điện chi phối mạnh mẽ, khiến lực ma sát phụ thuộc trực tiếp vào diện tích tiếp xúc thực tế chứ không còn độc lập với diện tích biểu kiến.

Xu hướng phát triển và hoàn thiện mô hình

Nhằm đáp ứng các yêu cầu ngày càng khắt khe của kỹ thuật chính xác, công nghiệp bán dẫn và mô phỏng cơ học hiện đại, các nhà khoa học đang phát triển các hướng tiếp cận mới bổ sung cho mô hình Coulomb cổ điển:

Mô hình ma sát động học phụ thuộc tốc độ và trạng thái (Rate-and-State Friction): Ban đầu được phát triển trong địa vật lý để mô phỏng động đất và chuyển động đứt gãy kiến tạo, mô hình này xem hệ số ma sát là một hàm của vận tốc trượt tức thời và một hoặc nhiều biến trạng thái nội tại phản ánh thời gian tiếp xúc của các đỉnh mấp mô vi mô. Mô hình này giải thích trọn vẹn hiện tượng chuyển tiếp từ ma sát tĩnh sang ma sát động và hành vi dính trượt không ổn định.
Tiếp cận ma sát học đa quy mô (Multiscale Tribology): Kết hợp mô phỏng động lực học phân tử (molecular dynamics) ở cấp độ nguyên tử để tính toán độ bền dính kết, mô phỏng cơ học tiếp xúc ngẫu nhiên ở cấp độ vi mô (rough surface contact mechanics) để xác định diện tích thực $A_r$ , và tích hợp các trường ứng suất này vào mô hình phần tử hữu hạn vĩ mô.
Thuật toán tiếp xúc phi tuyến tiên tiến trong phần mềm CAE: Các hệ thống tính toán kỹ thuật hiện đại ứng dụng thuật toán nhân tử Lagrange mở rộng (augmented Lagrangian method) kết hợp các luật ma sát chính quy hóa phi tuyến, cho phép mô phỏng chính xác các bài toán va chạm tốc độ cao, dập kim loại tấm và gia công cắt gọt mà không gặp phải các lỗi gián đoạn số học truyền thống.