Nghiên Cứu Ứng Dụng Phương Pháp Phần Tử Hữu Hạn (FEM) Trong Giải Bài Toán Tiếp Xúc Bằng ANSYS

Tóm tắt nghiên cứu

Nghiên cứu tập trung giải quyết câu hỏi trọng tâm: Làm thế nào để mô hình hóa và giải quyết chính xác bài toán tiếp xúc phi tuyến tính trong kết cấu cơ khí bằng Phương pháp Phần tử Hữu hạn (FEM) kết hợp công cụ mô phỏng ANSYS?

Về mặt phương pháp luận, nhóm nghiên cứu đã hệ thống hóa cơ sở giải tích phần tử hữu hạn từ khâu rời rạc hóa, lựa chọn hàm xấp xỉ Lagrange/Hermite, đến thiết lập ma trận độ cứng tổng thể. Trên cơ sở đó, nghiên cứu triển khai giải thuật lặp Newton-Raphson để xử lý tính phi tuyến hình học và phi tuyến biên tiếp xúc. Mô hình thực nghiệm số được xây dựng trên ngôn ngữ APDL của ANSYS để phân tích cặp truyền động trục vít – đai ốc, sử dụng phần tử bậc cao PLANE82 kết hợp cặp phần tử tiếp xúc bề mặt CONTA172 và TARGE169.

Kết quả nghiên cứu đã chỉ ra đặc tính hội tụ ổn định của mô hình tại vòng lặp thứ 51 dưới tác động của lược đồ 100 bước tải phụ (substeps). Phân tích trường ứng suất tương đương Von Mises cho thấy sự tập trung ứng suất cực đại đạt giá trị xấp xỉ $348\text{ MPa}$ khu trú nghiêm ngặt tại vùng chân ren và mặt sườn tiếp xúc chịu lực, trong khi ứng suất ở các vùng thân xa tiếp xúc nhanh chóng suy giảm xuống mức danh nghĩa ($0.239 \times 10^{-3}\text{ MPa}$). Ý nghĩa cốt lõi của nghiên cứu nằm ở việc chuẩn hóa quy trình phân tích số cho các cặp tiếp xúc cơ khí phức tạp, cung cấp công cụ dự báo độ bền và tuổi thọ mỏi bề mặt mà các công thức giải tích cổ điển không thể mô tả đầy đủ.

Bối cảnh và tầm quan trọng

Thực trạng tri thức và giới hạn của lý thuyết giải tích

Trong kỹ thuật cơ khí và xây dựng công trình, hiện tượng tiếp xúc giữa các bề mặt vật thể rắn đàn hồi xuất hiện phổ biến trong hầu hết các cụm chi tiết máy: sự ăn khớp bánh răng, bộ truyền trục vít – bánh vít, ổ lăn, mối ghép ren bu lông – đai ốc và các kết cấu dầm cầu. Về mặt vật lý, tiếp xúc ban đầu thường chỉ diễn ra tại một điểm hoặc một đường hẹp, sau đó dưới tác dụng của tải trọng và biến dạng đàn hồi sẽ mở rộng thành tiếp xúc mặt.

Lý thuyết tiếp xúc cổ điển của Hertz cung cấp nghiệm giải tích chính xác nhưng bị giới hạn nghiêm ngặt bởi các giả thiết lý tưởng: bề mặt trơn nhẵn, biến dạng rất nhỏ, vật liệu đàn hồi tuyến tính và dạng hình học đơn giản (mặt cầu, mặt trụ). Khi đối mặt với các biên dạng phức tạp, có góc lượn sắc, khe hở phân bố phi tuyến và chịu điều kiện ma sát thực tế, phương pháp giải tích truyền thống hoàn toàn bất khả thi.

Khoảng trống nghiên cứu và tính cấp thiết

Vấn đề tiếp xúc thực chất là bài toán phi tuyến biên (boundary nonlinearity), trong đó diện tích tiếp xúc và trạng thái ma sát biến thiên liên tục theo mức độ biến dạng. Ứng suất tiếp xúc có đặc tính cục bộ với gradient biến thiên cực lớn, là nguyên nhân hàng đầu gây nên hiện tượng mỏi bề mặt, tróc rỗ (pitting) và nứt tế vi dẫn đến phá hủy kết cấu. Do đó, việc ứng dụng các phương pháp số hiện đại như FEM và bộ giải công nghiệp ANSYS là yêu cầu cấp thiết để rút ngắn chu kỳ thiết kế, hạ giá thành thử nghiệm và bảo đảm an toàn vận hành.

Methodology và approach

Nghiên cứu áp dụng quy trình thiết kế và phân tích số toàn diện gồm 4 giai đoạn chính:

1. Cơ sở lý thuyết phần tử hữu hạn (FEM)

Miền liên tục thể tích $V$ được rời rạc hóa thành $N_e$ miền con $V_e$ kết nối tại các điểm nút. Ẩn hàm chuyển vị $u(x)$ được xấp xỉ hóa qua ma trận hàm dạng $[P(x)]$ và véc-tơ chuyển vị nút ${q}^e$:

$$u(x) = [P(x)]{a} = [N(x)]{q}^e$$

Hệ phương trình cân bằng tĩnh tổng thể được thiết lập qua nguyên lý năng lượng cực tiểu:

$$[K]{q} = {P}$$

Trong đó:

  • $[K]$: Ma trận độ cứng tổng thể cấu trúc.
  • ${q}$: Véc-tơ chuyển vị nút tổng thể.
  • ${P}$: Véc-tơ tải trọng tổng thể đã áp đặt điều kiện biên.

2. Mô hình hóa hình học và tiền xử lý trên ANSYS

  • Hình học: Sử dụng mô hình hình học 2D đối xứng trục nhập từ tệp trung gian IGES của cặp ren trục vít – đai ốc, giảm thiểu khối lượng tính toán nhưng vẫn bảo toàn bản chất ứng xử không gian 3 chiều.
  • Lựa chọn phần tử:
    • Phần tử kết cấu nền: PLANE82 (phần tử 2D bậc cao 8 nút với 2 bậc tự do chuyển vị tại mỗi nút $u_x, u_y$), thích hợp mô phỏng trạng thái biến dạng phẳng hoặc đối xứng trục với độ cong lớn.
    • Cặp phần tử tiếp xúc bề mặt: CONTA172 (phần tử tiếp xúc 2D gắn trên mặt vít) và TARGE169 (phần tử mục tiêu gắn trên mặt đai ốc).
  • Vật liệu: Thép đàn hồi đồng chất đẳng hướng với mô đun Young $E = 200 \times 10^3\text{ MPa}$ ($200\text{ GPa}$) và hệ số Poisson $\nu = 0.3$.
  • Kỹ thuật chia lưới (Meshing): Thiết lập kích thước phần tử tự điều chỉnh thích ứng, tập trung làm mịn cục bộ (mesh refinement) tại các khía ren và chân ren – nơi có gradient ứng suất biến thiên đột ngột.
+-----------------------------------------------------------+
|             BẢNG THÔNG SỐ CẤU HÌNH MÔ HÌNH SỐ             |
+----------------------+------------------------------------+
| Thông số             | Giá trị cài đặt                    |
+----------------------+------------------------------------+
| Loại phân tích       | Tĩnh phi tuyến tiếp xúc (Nonlinear)|
| Phần tử cơ sở        | PLANE82 (2D 8-Node Structural)    |
| Cặp tiếp xúc         | CONTA172 (Contact) & TARGE169 (Target)|
| Mô đun đàn hồi (E)   | 200,000 MPa                        |
| Hệ số Poisson (ν)    | 0.3                                |
| Chuyển vị ép dọc trục| ux = 0.1 mm                        |
| Thời gian tải (Time) | 10 s (Chia 10 - 100 substeps)      |
| Hiệu ứng hình học    | NLGEOM, ON (Large Deflection)      |
+----------------------+------------------------------------+

3. Kỹ thuật giải phi tuyến và độ tin cậy

Do bài toán tiếp xúc có quan hệ lực – chuyển vị thay đổi liên tục $[K(D)]{D} = {F}$, nghiên cứu sử dụng thuật giải lặp Newton-Raphson kết hợp tùy chọn kích hoạt biến dạng lớn (NLGEOM, ON). Quá trình đặt tải được chia thành 100 bước tải phụ (substeps) tự động theo thời gian $t = 10\text{ s}$ nhằm bảo đảm ma trận độ cứng tiếp tuyến không bị phân kỳ. Tiêu chuẩn hội tụ kiểm soát đồng thời chuẩn lực dư $L_2$ norm và ngưỡng dung sai cho phép (CRIT).

Phát hiện chính

Quá trình mô phỏng và phân tích số liệu đã mang lại 4 phát hiện khoa học quan trọng:

1. Đặc tính hội tụ của giải thuật tiếp xúc phi tuyến

Biểu đồ hội tụ lực dư khẳng định thuật giải lặp đạt trạng thái cân bằng số học ở vòng lặp thứ 51. Đường cong chuẩn $L_2$ (lực mất cân bằng) tiệm cận nhanh chóng và nằm hoàn toàn dưới đường dung sai tiêu chuẩn (CRIT). Việc phân chia bước tải phụ nhỏ ($0.1\text{ s}$/bước) đã triệt tiêu hiện tượng dao động số học thường gặp khi bề mặt tiếp xúc bắt đầu va chạm hoặc trượt.

2. Sự tập trung ứng suất cục bộ Von Mises cực đại

Kết quả trường ứng suất tương đương Von Mises chỉ ra:

  • Ứng suất cực đại ($\sigma_{\max} \approx 348\text{ MPa}$): Tập trung cục bộ cao độ tại góc lượn chân ren và các sườn tiếp xúc chịu lực trực tiếp. Đây là vị trí chịu tương tác đồng thời giữa lực ép pháp tuyến và ứng suất cắt tiếp tuyến do ma sát.
  • Ứng suất danh nghĩa ở vùng xa ($\sigma_{\min} \approx 0.239 \times 10^{-3}\text{ MPa}$): Vùng thân vít và đai ốc bên ngoài miền tiếp xúc duy trì mức ứng suất rất thấp, khẳng định tính đúng đắn của nguyên lý Saint-Venant mở rộng cho bài toán cơ học tiếp xúc.

3. Quy luật phân bố phản lực và độ võng tiếp xúc

Biến dạng tiếp xúc không phân bố đều trên toàn bộ chiều dài mối ren mà tập trung chủ yếu vào 2 đến 3 bước ren đầu tiên gần mặt tựa. Hiện tượng này giải thích nguyên nhân các hư hỏng trượt ren, đứt ren trong thực tế luôn khởi phát từ vòng ren chịu tải đầu tiên.

4. Tính chuẩn xác của mô hình phần tử hữu hạn so với giải tích

Thông qua bài toán chuẩn hóa hệ thanh chịu tải dọc trục phân bố đều, nghiên cứu chứng minh chuyển vị nút thu được từ FEM trùng khớp tuyệt đối với nghiệm giải tích chính xác ($u(x) = \frac{q}{2EF}(4ax - x^2)$). Nghiệm nội lực $N$ dạng bậc thang tiệm cận trơn tru về nghiệm giải tích khi mật độ phần tử được tăng cường.

Đóng góp khoa học

1. Về mặt lý thuyết và phương pháp luận

  • Hệ thống hóa hoàn chỉnh quy trình chuyển đổi bài toán cơ học tiếp xúc phi tuyến sang mô hình số rời rạc có khả năng tự động hóa trên máy tính.
  • Làm sáng tỏ cơ chế cập nhật ma trận độ cứng phi tuyến $[K(D)]$ dưới tác động của sự thay đổi điều kiện biên tiếp xúc.

2. Về mặt thực tiễn kỹ thuật

  • Xây dựng bộ mã nguồn điều khiển linh hoạt trên nền tảng ANSYS APDL, cho phép các kỹ sư thiết kế nhanh chóng tái tạo mô hình, thay đổi bước ren, góc ren và thông số ma sát mà không cần vẽ lại thủ công.
  • Cung cấp cơ sở định lượng để tối ưu hóa góc lượn chân ren, giúp phân tán ứng suất đỉnh và nâng cao khả năng chịu tải mỏi của các bộ truyền trục vít cơ khí trong công nghiệp chế tạo máy.

Đối tượng quan tâm và giá trị ứng dụng

Câu hỏi thường gặp (FAQ)

1. Phát hiện quan trọng nhất của nghiên cứu này là gì?

Nghiên cứu đã định lượng hóa chính xác trường ứng suất Von Mises cực đại ($348\text{ MPa}$) tập trung tại rãnh chân ren và xác lập quy luật hội tụ số ổn định ở vòng lặp thứ 51 thông qua thuật giải Newton-Raphson.

2. Điểm đặc biệt trong phương pháp tiếp cận mô hình hóa là gì?

Đó là sự kết hợp tối ưu giữa phần tử phẳng bậc cao 8 nút PLANE82 với cặp phần tử tiếp xúc mặt – mặt CONTA172/TARGE169, cho phép mô phỏng chính xác ứng xử tiếp xúc đối xứng trục mà không tiêu tốn quá nhiều tài nguyên tính toán như mô hình 3D khối đặc.

3. Kết quả nghiên cứu có thể tổng quát hóa cho các chi tiết khác không?

Hoàn toàn có thể. Quy trình tiền xử lý, thuật toán thiết lập cặp tiếp xúc và cấu hình kiểm soát bước tải phụ có thể áp dụng trực tiếp cho các bài toán tiếp xúc phức tạp khác như: ăn khớp bánh răng, chèn ép ổ bi đỡ chặn, đĩa phanh ô tô hay liên kết bu lông cường độ cao.

4. Những hướng phát triển tiếp theo của đề tài là gì?

Các hướng mở rộng tiềm năng gồm:

  1. Nâng cấp lên mô hình 3D hoàn chỉnh có xét đến góc nâng ren thực tế.
  2. Tích hợp mô hình vật liệu đàn dẻo phá hủy (Elasto-plasticity) và từ biến.
  3. Khảo sát bài toán tương tác đa trường vật lý (Coupled Thermal-Structural) khi ma sát tiếp xúc sinh nhiệt lớn ở vận tốc cao.

5. Ứng dụng thực tế nổi bật nhất của kết quả này là gì?

Ứng dụng trong việc thiết kế và kiểm nghiệm độ bền cho các cơ cấu trục vít – đai ốc truyền động chính xác, kích nâng cơ khí, vít me máy công cụ CNC và các bộ ghép ren chịu lực trong ngành công nghiệp giao thông vận tải.

Kết luận

Báo cáo nghiên cứu đã chứng minh tính hiệu quả vượt trội của Phương pháp Phần tử Hữu hạn (FEM) và phần mềm ANSYS trong việc giải quyết bài toán tiếp xúc phi tuyến tính phức tạp trên cặp cơ cấu trục vít – đai ốc. Bằng việc phân tích sâu sắc từ bản chất toán học đến thực nghiệm mô phỏng số, công trình đã xác lập một quy trình phân tích ứng suất và biến dạng tiếp xúc có độ tin cậy cao. Kết quả này không chỉ đóng góp tài liệu tham khảo học thuật giá trị cho sinh viên và kỹ sư cơ khí mà còn mở ra công cụ đắc lực hỗ trợ tính toán, thiết kế tối ưu hóa kết cấu máy móc hiện đại.