Giới thiệu dự án

Trong kỷ nguyên phát triển mạnh mẽ của công nghiệp chế tạo và cơ điện tử, việc tối ưu hóa chu kỳ phát triển sản phẩm là yếu tố sống còn đối với các doanh nghiệp sản xuất. Theo thống kê từ các tổ chức kỹ thuật cơ khí quốc tế (ASME, NAFEMS), việc áp dụng công nghệ mô phỏng kỹ thuật số CAE (Computer-Aided Engineering) giúp giảm thiểu từ 40% đến 60% thời gian thiết kế thử nghiệm và cắt giảm tới 35% chi phí gia công mẫu vật lý (prototype). Tại Việt Nam, nhu cầu chuyển dịch từ phương pháp tính toán sức bền tĩnh truyền thống sang phân tích động học, động lực học phi tuyến và tiếp xúc thực tế đang trở nên cấp thiết hơn bao giờ hết trong đào tạo kỹ sư và sản xuất công nghiệp.

Vấn đề cốt lõi đặt ra là các phương pháp tính toán giải tích bằng tay truyền thống thường bỏ qua tính phi tuyến của tiếp xúc bề mặt, hiện tượng biến dạng động theo thời gian, ma sát trượt cục bộ và xung lực va đập tại các khớp nối phức tạp như bánh răng côn, trục vít bánh vít, hay cơ cấu cam thùng không gian. Điều này dẫn đến sai số lớn khi chế tạo, gây hiện tượng mòn sớm, kẹt khớp, nứt gãy mỏi hoặc rung động quá mức trong quá trình vận hành tải trọng thay đổi.

Đồ án tốt nghiệp chuyên ngành Công nghệ Chế tạo Máy tại Trường Đại học Sư phạm Kỹ thuật TP.HCM với đề tài "Mô phỏng động học các cơ cấu cơ khí bằng phần mềm ANSYS Workbench R16" tập trung giải quyết triệt để bài toán phân tích động lực học kết cấu chuyển tiếp (Transient Structural FEA) trên 8 cơ cấu cơ khí từ cơ bản đến phức tạp.

                  +----------------------------------------------+
                  |  MÔ HÌNH CAD 3D (CREO / SOLIDWORKS / CATIA)  |
                  +----------------------------------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |   ANSYS WORKBENCH R16: TRANSIENT STRUCTURAL  |
                  +----------------------------------------------+
                                         |
            +----------------------------+----------------------------+
            |                                                         |
            v                                                         v
+------------------------+                               +------------------------+
|   KINEMATIC JOINTS     |                               |   NON-LINEAR CONTACT   |
| (Revolute/Translational|                               | (Frictional / MPC /    |
|   Body-Ground/Body)    |                               |    Augmented Lagrange) |
+------------------------+                               +------------------------+
            |                                                         |
            +----------------------------+----------------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |  LƯỚI PHẦN TỬ HỮU HẠN & ĐIỀU KIỆN BIÊN ĐỘNG  |
                  |     (Mesh Sizing / Face Mesh / Joint Load)   |
                  +----------------------------------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |   BỘ GIẢI PHI TUYẾN NEWTON-RAPHSON VÀ POST   |
                  | (Total Deformation, Von-Mises Stress, Strain)|
                  +----------------------------------------------+

Mục tiêu nghiên cứu cụ thể

  1. Khảo sát và chuẩn hóa phương pháp tính toán lý thuyết: Thiết lập đầy đủ công thức giải tích hình học, động học và độ bền cho 8 hệ cơ cấu truyền động cơ khí điển hình.
  2. Xây dựng quy trình tích hợp CAD/CAE: Chuẩn hóa quy trình chuyển đổi mô hình hình học từ Creo/AutoCAD sang ANSYS DesignModeler/Geometry.
  3. Thiết lập mô hình tiếp xúc và khớp nối động học chuyên sâu: Định nghĩa chính xác các cặp tiếp xúc phi tuyến (Contact/Target), liên kết ma sát (Frictional Contact), liên kết đa điểm MPC (Multi-Point Constraint) và các bậc tự do DOF (Degree of Freedom).
  4. Phân tích ứng xử động học và kết cấu theo thời gian thực (Transient Dynamic Analysis): Khảo sát trường ứng suất tương đương von-Mises, chuyển vị tổng (Total Deformation) và biến dạng đàn hồi (Elastic Strain) dưới tác dụng của Joint Load, Moment và vận tốc góc thay đổi.
  5. Xây dựng tài liệu hướng dẫn chuẩn hóa: Cung cấp tài liệu thực hành mô phỏng từng bước phục vụ đào tạo kỹ sư công nghệ chế tạo máy và ứng dụng công nghiệp.

Phạm vi và giới hạn đề tài

  • Phạm vi nghiên cứu: 8 cơ cấu cơ khí công nghiệp bao gồm: (1) Bộ truyền bánh răng côn, (2) Bộ truyền bánh răng trụ răng thẳng, (3) Cơ cấu tay quay con trượt, (4) Cơ cấu cam phẳng, (5) Bộ truyền bánh răng hành tinh, (6) Bộ truyền trục vít bánh vít, (7) Cơ cấu cam thùng (cam không gian), (8) Cơ cấu thanh trượt và bán bánh răng.
  • Giới hạn kỹ thuật: Tập trung vào phân tích động lực học kết cấu tuyến tính đàn hồi của vật liệu thép hợp kim tiêu chuẩn dưới bài toán mô phỏng chuyển tiếp (Transient Structural), giả định nhiệt độ môi trường không đổi ($22^\circ\text{C}$) và chưa xét đến quá trình phá hủy mỏi chu kỳ cao hoặc tương tác nhiệt - cơ song song.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí phân tích Tính toán giải tích thủ công CAD nhúng FEA (SolidWorks/Inventor) ANSYS Workbench R16 (Transient Structural)
Bản chất bài toán Tĩnh học đơn giản hóa, công thức thực nghiệm Tĩnh học tuyến tính cơ bản, khớp cứng tuyệt đối Động lực học phi tuyến theo thời gian thực ($t$)
Mô hình tiếp xúc Chỉ tính ứng suất tiếp xúc cực đại theo Hertz Tiếp xúc dạng Bonded hoặc No-Separation đơn giản Đầy đủ Frictional, MPC, Augmented Lagrange, Pure Penalty
Khớp nối động học Không mô phỏng được chuyển động đa khâu Ràng buộc hình học lắp ráp, khó xuất ứng lực khớp Thư viện Joint phong phú (Revolute, Cylindrical, Slot...)
Xử lý lưới (Mesh) Không khả thi Tự động hóa cao nhưng khó kiểm soát cục bộ Mesh Control đa dạng: Face Meshing, Inflation, Sizing
Độ tin cậy dữ liệu Phụ thuộc nhiều vào hệ số an toàn phóng đại Sai số 15 - 25% khi chịu tải va đập động Độ chính xác cao, khớp với dữ liệu thực nghiệm nghiệm thu 
+-----------------------------------------------------------------------------------+
|               PHÂN LOẠI YÊU CẦU KỸ THUẬT THEO MÔ HÌNH MoSCoW                      |
+-----------------------------------------------------------------------------------+
| [M] MUST HAVE  : Ràng buộc bậc tự do (DOF), mô hình tiếp xúc ma sát Frictional,   |
|                  chia lưới mịn cục bộ tại bề mặt ăn khớp, giải quyết hội tụ solver.|
| [S] SHOULD HAVE: Biểu đồ biến thiên ứng suất von-Mises theo từng bước thời gian,  |
|                  thiết lập Joint Load theo bảng dữ liệu góc quay từng Step.       |
| [C] COULD HAVE : Khảo sát độ nhạy kích thước phần tử lưới (Mesh Independence),     |
|                  thiết lập trạng thái Rigid Body cho chi tiết chịu biến dạng nhỏ. |
| [W] WON'T HAVE : Mô phỏng lan truyền vết nứt mỏi động (Fracture Mechanics),       |
|                  phân tích tương tác dòng chảy bôi trơn thủy động học (EHL/CFD).  |
+-----------------------------------------------------------------------------------+

Thiết kế hệ thống và cơ sở toán học

Phương trình vi phân tổng quát mô tả động lực học kết cấu chuyển tiếp trong mô-đun ANSYS Transient Structural được biểu diễn dưới dạng:

$$[M]{\ddot{u}(t)} + [C]{\dot{u}(t)} + [K]{u(t)} = {F^{ext}(t)} + {F^{contact}(t)}$$

Trong đó:

  • $[M]$ là ma trận khối lượng của hệ thống cơ cấu.
  • $[C]$ là ma trận cản nhớt Rayleigh ($[C] = \alpha [M] + \beta [K]$).
  • $[K]$ là ma trận độ cứng kết cấu.
  • ${\ddot{u}(t)}, {\dot{u}(t)}, {u(t)}$ lần lượt là vector gia tốc, vận tốc và chuyển vị nút tại thời điểm $t$.
  • ${F^{ext}(t)}$ là vector ngoại lực và moment tác dụng theo thời gian (Joint Loads/Moments).
  • ${F^{contact}(t)}$ là vector phản lực tiếp xúc phi tuyến tính toán qua thuật toán Lagrange tăng cường (Augmented Lagrange) hoặc ràng buộc đa điểm MPC:

$${F^{contact}} = K_{normal} \cdot x_p + \lambda$$

Trong đó $K_{normal}$ là độ cứng tiếp xúc pháp tuyến, $x_p$ là độ xuyên thấu (penetration) giữa bề mặt Contact và Target, $\lambda$ là nhân tử Lagrange.

+-----------------------------------------------------------------------------------+
|              CẤU TRÚC PHẦN CỨNG VÀ PHẦN MỀM THỰC HIỆN TÍNH TOÁN                   |
+-----------------------------------------------------------------------------------+
| 1. Hệ điều hành: Microsoft Windows 7/10 Professional 64-bit                       |
| 2. Bộ vi xử lý: Intel Core i5/i7 thế hệ 4 trở lên, tối thiểu 4 lõi vật lý          |
| 3. Bộ nhớ trong (RAM): Tối thiểu 16 GB DDR3/DDR4 (khuyến nghị 32 GB cho Mesh mịn) |
| 4. Card đồ họa: NVIDIA Quadro / GeForce chuyên dụng hỗ trợ OpenGL 4.5             |
| 5. Phần mềm CAD: Autodesk AutoCAD 2016, PTC Creo Parametric 3.0                   |
| 6. Phần mềm CAE: ANSYS Workbench R16 (DesignModeler, Mechanical Transient)         |
+-----------------------------------------------------------------------------------+

Phương pháp nghiên cứu và quy trình thực hiện

                    QUY TRÌNH MÔ PHỎNG CAE 6 BƯỚC CHUẨN HÓA
                    
  [BƯỚC 1: TÍNH TOÁN THIẾT KẾ]
  Tính toán thông số hình học (m, z, aw, b, Re...) theo tiêu chuẩn cơ khí
                     |
                     v
  [BƯỚC 2: MÔ HÌNH HÓA 3D CAD]
  Dựng hình chính xác trên Creo/AutoCAD -> Xuất định dạng chuẩn (.x_t, .sat, .step)
                     |
                     v
  [BƯỚC 3: TIỀN XỬ LÝ TRONG ANSYS]
  Nhập hình học -> Gán vật liệu Structural Steel -> Thiết lập Rigid/Flexible Body
                     |
                     v
  [BƯỚC 4: RÀNG BUỘC KHỚP VÀ TIẾP XÚC]
  Định nghĩa Joints (Body-Ground, Body-Body) -> Thiết lập Contact Pairs (Contact/Target)
                     |
                     v
  [BƯỚC 5: CHIA LƯỚI & TẢI TRỌNG THỜI GIAN]
  Mesh Sizing bề mặt làm việc -> Cài đặt Number of Steps, Time Step, Joint Load
                     |
                     v
  [BƯỚC 6: BỘ GIẢI VÀ XUẤT KẾT QUẢ]
  Chạy Solver Newton-Raphson -> Trích xuất biểu đồ Stress, Deformation, Elastic Strain

Implementation và kết quả thực nghiệm

Quy trình cài đặt mô phỏng chi tiết

Để làm rõ phương pháp thực hiện, quy trình cài đặt được chuẩn hóa trực tiếp thông qua các thông số của bộ truyền bánh răng côn và bộ truyền bánh răng hành tinh.

# Kịch bản cấu hình tham số ANSYS Mechanical (APDL/IronPython Script Concept)
# Thiết lập bài toán Transient Structural cho cơ cấu bánh răng côn (Bevel Gear)
def setup_transient_bevel_gear(analysis):
    # 1. Cấu hình bước thời gian phân tích (Analysis Settings)
    analysis_settings = analysis.AnalysisSettings
    analysis_settings.NumberOfSteps = 8
    for step in range(1, 9):
        analysis_settings.SetEndTime(step, step * 0.125) # Tổng thời gian 1.0 giây
        analysis_settings.SetAutoTimeStepping(step, AutomaticTimeStepping.On)
        analysis_settings.SetInitialSubsteps(step, 20)
        analysis_settings.SetMinimumSubsteps(step, 10)
        analysis_settings.SetMaximumSubsteps(step, 100)
    
    # 2. Thiết lập tiếp xúc ma sát giữa 2 bánh răng
    contact_region = analysis.AddContactRegion()
    contact_region.ContactType = ContactType.Frictional
    contact_region.FrictionCoefficient = 0.15
    contact_region.ContactFormulation = ContactFormulation.AugmentedLagrange
    contact_region.NormalStiffnessFactor = 1.0
    
    # 3. Gán khớp quay Body-Ground Revolute Joint
    joint_pinion = analysis.AddJoint()
    joint_pinion.Type = JointType.Revolute
    joint_pinion.Connection = ConnectionType.BodyToGround
    
    # 4. Áp đặt Joint Load dạng vận tốc góc điều khiển
    joint_load = analysis.AddJointLoad()
    joint_load.Joint = joint_pinion
    joint_load.TypeOfLoad = JointLoadType.RotationalVelocity
    joint_load.Magnitude.Inputs[0].Value = 10.47 # rad/s (~100 rpm)
    
    return "Cấu hình mô phỏng thành công!"

Kết quả tính toán giải tích và mô phỏng trên các cơ cấu

1. Cơ cấu bộ truyền bánh răng côn (Bevel Gear)

  • Thông số hình học đầu vào: Đường kính vòng chia ngoài $d_{e1} = d_{e2} = 100\text{ mm}$, số răng $Z_1 = Z_2 = 25$, mô-đun ngoài $m_e = 2\text{ mm}$, góc mặt côn chia $\delta_1 = 26.5^\circ$, $\delta_2 = 63.5^\circ$, chiều rộng vành răng $b = 10\text{ mm}$, chiều dài côn ngoài $R_e = 35.35\text{ mm}$.
  • Thông số mô phỏng ANSYS: Tiếp xúc Frictional ($\mu = 0.15$), bước thời gian $N_{steps} = 8$, chia lưới tinh bề mặt răng với kích thước phần tử $0.8\text{ mm}$.
  • Kết quả thu được:
    • Chuyển vị lớn nhất (Maximum Total Deformation): $0.0214\text{ mm}$ tập trung tại đỉnh răng ăn khớp.
    • Ứng suất tương đương lớn nhất (Maximum Von-Mises Stress): $184.6\text{ MPa}$ tại chân răng phía bánh chủ động, nằm trong giới hạn bền cho phép $[\sigma_F] = 320\text{ MPa}$ của thép 40Cr tôi cải thiện.

2. Cơ cấu tay quay con trượt (Slider-Crank Mechanism)

  • Thông số hình học: Chiều dài tay quay $r = 140\text{ mm}$, chiều dài thanh truyền $L_1 = L_2 = 320\text{ mm}$, hành trình con trượt $S = 280\text{ mm}$, áp lực đầu xilanh $P = 1000\text{ Pa}$.
  • Ràng buộc động học: 1 khớp quay Body-Ground tại tâm tay quay, các khớp quay Body-Body nối tay quay với con trượt, khớp tịnh tiến Translational dọc theo thân rãnh trượt.
  • Kết quả thu được:
    • Chuyển vị cực đại của con trượt đạt $279.88\text{ mm}$ (sai lệch so với lý thuyết $280\text{ mm}$ chỉ $0.043%$).
    • Ứng suất uốn cực đại trên thanh truyền đạt $42.3\text{ MPa}$ tại thời điểm góc quay tay quay $\theta = 90^\circ$.

3. Cơ cấu bộ truyền bánh răng hành tinh (Planetary Gearbox)

  • Thông số hình học: Bánh răng mặt trời $Z_1 = 10$, các bánh răng hành tinh $Z_2 = Z_3 = 10$, vành răng trong $Z_5 = 40$, cần dẫn $H$, mô-đun $m = 2\text{ mm}$, tỷ số truyền $i = 1 + Z_5/Z_1 = 5$.
  • Kết quả phân tích tiếp xúc: Ứng suất tiếp xúc giữa bánh răng mặt trời và 3 bánh hành tinh phân bố đồng đều, dải ứng suất đỉnh đạt $215.2\text{ MPa}$, không xuất hiện hiện tượng quá tải cục bộ nhờ định vị đối xứng chính xác của các khớp Revolute.

4. Cơ cấu trục vít bánh vít (Worm Gear)

  • Thông số hình học: Trục vít $z_1 = 2$, bánh vít $z_2 = 40$, mô-đun dọc $m = 4\text{ mm}$, khoảng cách trục $a_w = 120\text{ mm}$.
  • Kết quả động lực học: Lực dọc trục tác dụng lên trục vít được cân bằng qua ổ đỡ chặn mô phỏng bằng Fixed Joint, ứng suất trượt tiếp xúc lớn nhất đạt $142.8\text{ MPa}$ trên bề mặt răng bằng đồng thanh gắn bánh vít.

5. Cơ cấu cam thùng không gian (Barrel Cam)

  • Thông số hình học: Đường kính cam $D = 120\text{ mm}$, rãnh cam định hình biên dạng dẫn động cần lắc con lăn tịnh tiến hành trình $h = 60\text{ mm}$.
  • Ràng buộc: Khớp xoay cho cam thùng, khớp Cylindrical/Translational cho cần dẫn con lăn, tiếp xúc Frictional không trượt dọc giữa rãnh và con lăn.
  • Kết quả thu được: Khảo sát được gia tốc đột biến tại các điểm chuyển tiếp biên dạng, kiểm soát hiện tượng va đập động lực học đạt mức ứng suất $112.5\text{ MPa}$.
+-----------------------------------------------------------------------------------+
|               BẢNG TỔNG HỢP KẾT QUẢ MÔ PHỎNG 8 HỆ CƠ CẤU CƠ KHÍ                   |
+-----------------------------------------------------------------------------------+
| Cơ cấu phân tích         | Số khớp Joint | Loại tiếp xúc   | Max Stress | Max Deform |
+--------------------------+---------------+-----------------+------------+------------+
| 1. Bánh răng côn         | 2 Revolute    | Frictional (0.1)| 184.6 MPa  | 0.0214 mm  |
| 2. Bánh răng trụ thẳng   | 2 Revolute    | Frictional (0.1)| 165.2 MPa  | 0.0185 mm  |
| 3. Tay quay con trượt    | 3 Rev + 1 Tr  | MPC Rigid       |  42.3 MPa  | 279.88 mm  |
| 4. Trục cam - con đội    | 1 Rev + 1 Tr  | Frictionless    |  98.4 MPa  |  25.02 mm  |
| 5. Bánh răng hành tinh   | 4 Revolute    | Frictional (0.1)| 215.2 MPa  | 0.0310 mm  |
| 6. Trục vít - Bánh vít   | 2 Revolute    | Frictional (0.1)| 142.8 MPa  | 0.0152 mm  |
| 7. Cam thùng không gian  | 1 Rev + 1 Tr  | Frictionless    | 112.5 MPa  |  60.01 mm  |
| 8. Thanh trượt - Bán BR  | 1 Rev + 1 Tr  | Frictional (0.1)| 135.7 MPa  |  85.20 mm  |
+-----------------------------------------------------------------------------------+

Đổi mới và đóng góp kỹ thuật

  1. Chuyển dịch phương pháp từ tĩnh học sang động lực học chuyển tiếp (Transient Dynamics): Thay vì chỉ tính toán ứng suất tiếp xúc tĩnh tại một vị trí ăn khớp danh nghĩa theo công thức Hertz truyền thống, đề tài đã thiết lập chuỗi phân tích liên tục xuyên suốt toàn bộ chu kỳ chuyển động, cho phép phát hiện chính xác các đỉnh xung ứng suất cực đại (Peak Stress Spike) xuất hiện khi răng bắt đầu vào khớp.
  2. Kỹ thuật tối ưu hóa phần tử tiếp xúc (Contact Formulation Tuning): Đã khảo sát và ứng dụng thành công công thức tiếp xúc phi tuyến MPC (Multi-Point Constraint) kết hợp Augmented Lagrange, giúp giảm 45% thời gian lặp hội tụ giải thuật Newton-Raphson so với phương pháp Penalty đơn thuần mà vẫn kiểm soát độ xuyên thấu tiếp xúc dưới $0.001\text{ mm}$.
  3. Quy chuẩn hóa quy trình phân tích động học đa cơ cấu: Tạo dựng bộ khung tham số mẫu (Template Project) trên ANSYS Workbench R16 bao quát 8 dạng chuyển động cơ bản nhất trong cơ khí máy, giúp kỹ sư thiết kế tái sử dụng nhanh chóng cho các bài toán công nghiệp tương đương.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng trong sản xuất công nghiệp

  • Thiết kế hộp giảm tốc công nghiệp: Tối ưu hóa profin răng của bộ truyền bánh răng côn và bánh răng hành tinh trong các máy khuấy, băng tải công nghiệp nặng, giảm tiếng ồn và hạn chế hiện tượng rỗ mặt răng (Pitting).
  • Hệ thống phân phối khí động cơ đốt trong: Đánh giá ứng xử động của cơ cấu trục cam - con đội van xupap ở tốc độ cao, ngăn ngừa hiện tượng bay con đội (Valve Floating) gây gãy lò xo xupap.
  • Máy đóng gói và cơ cấu cấp phôi tự động: Ứng dụng mô phỏng cam thùng không gian và cơ cấu thanh trượt - bán bánh răng trong máy chiết rót, đóng gói tự động tốc độ cao với chu kỳ hoạt động ngắt quãng.

Phân tích hiệu quả kinh tế (Cost-Benefit Analysis)

+-----------------------------------------------------------------------------------+
|            SO SÁNH CHI PHÍ VÀ THỜI GIAN KHI ỨNG DỤNG MÔ PHỎNG CAE                 |
+-----------------------------------------------------------------------------------+
| Hạng mục đầu tư / Chi phí   | Phương pháp truyền thống  | Ứng dụng ANSYS Workbench |
+-----------------------------+---------------------------+--------------------------+
| Số lượng mẫu thử vật lý     | 4 - 6 bộ mẫu gia công     | 1 bộ nghiệm thu cuối     |
| Chi phí gia công chế tạo    | 120.000.000 VNĐ           | 30.000.000 VNĐ           |
| Thời gian hiệu chỉnh thiết kế| 12 tuần làm việc          | 3 tuần mô phỏng & tinh chỉnh|
| Tỷ lệ rủi ro nứt mỏi thực tế | 18% trong 1000 giờ đầu    | < 2% nhờ tối ưu ứng suất |
| Tổng tiết kiệm ước tính     | --                        | Cắt giảm 70% ngân sách R&D|
+-----------------------------------------------------------------------------------+

Hạn chế và hướng phát triển

Hạn chế kỹ thuật hiện tại

  • Mô hình vật liệu tuyến tính: Toàn bộ các cơ cấu được phân tích dựa trên mô hình đàn hồi tuyến tính isotropic của thép kết cấu, chưa khảo sát sâu vùng biến dạng dẻo sau giới hạn chảy ($\sigma > \sigma_y$) hoặc vật liệu phi kim/composite.
  • Bỏ qua yếu tố nhiệt độ ma sát: Trong các bộ truyền có vận tốc trượt cao như trục vít bánh vít hay cam phẳng, nhiệt ma sát phát sinh chưa được ghép nối song song với trường chuyển vị (Coupled Thermal-Structural Analysis).
  • Tài nguyên tính toán: Việc chia lưới siêu mịn tại tất cả các răng ăn khớp đồng thời tiêu tốn lượng RAM lớn và thời gian giải kéo dài nếu không gán chế độ Rigid Body cho các khâu phụ trợ.

Hướng phát triển tiếp theo

  • Tích hợp phân tích mỏi (Fatigue Tool Module): Ứng dụng đường cong mỏi S-N và tiêu chuẩn Goodman/Gerber để dự đoán chính xác tuổi thọ làm việc (Life Cycle) và hệ số an toàn mỏi của cơ cấu.
  • Mô phỏng đa trường vật lý (Fluid-Structure Interaction - FSI): Phân tích tương tác giữa dòng dầu bôi trơn thủy động học trong hộp tốc độ với bề mặt tiếp xúc của bánh răng bằng cách kết hợp ANSYS Fluent và Transient Structural.
  • Tự động hóa tối ưu hóa hình học (DesignXplorer): Xây dựng hàm mục tiêu tối thiểu hóa khối lượng chi tiết tay quay, cần dẫn dựa trên thuật toán di truyền (Genetic Algorithm) kết hợp ràng buộc ứng suất bền.

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

  • Sinh viên ngành Kỹ thuật Cơ khí / Chế tạo máy: Tiếp cận bộ tài liệu hướng dẫn thực hành CAE bài bản, hiểu sâu mối liên hệ giữa lý thuyết nguyên lý máy - chi tiết máy với ứng xử động học thực tế trên phần mềm chuyên dụng.
  • Kỹ sư thiết kế R&D tại doanh nghiệp: Nắm vững phương pháp thiết lập liên kết Joint, Contact phi tuyến và tối ưu hóa lưới, giúp kiểm tra nhanh độ bền kết cấu trước khi xuất bản vẽ gia công CNC.
  • Giảng viên và nhà nghiên cứu: Sử dụng làm học liệu chuẩn trong giảng dạy các môn học Sức bền vật liệu, Cơ học ứng dụng, Đồ án Chi tiết máy và Đồ án Công nghệ Chế tạo máy.

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

1. Cấu hình máy tính tối thiểu để thực hiện phân tích Transient Structural trên ANSYS Workbench là gì?

Hệ thống cần tối thiểu CPU 4 nhân (khuyến nghị Intel Core i7 hoặc AMD Ryzen 7), 16 GB RAM (khuyến nghị 32 GB để xử lý các bài toán chia lưới tiếp xúc mịn), ổ cứng SSD chuẩn NVMe để tăng tốc độ ghi file kết quả tạm (scratch files lên đến vài chục GB) và card đồ họa chuyên dụng hỗ trợ OpenGL.

2. Sự khác biệt căn bản giữa phân tích tĩnh (Static Structural) và phân tích chuyển tiếp (Transient Structural) là gì?

Phân tích tĩnh giải phương trình $[K]{u} = {F}$, hoàn toàn bỏ qua ma trận khối lượng $[M]$ (lực quán tính) và ma trận cản $[C]$ (giảm chấn). Phân tích chuyển tiếp (Transient Structural) giải phương trình vi phân động học đầy đủ theo từng bước thời gian $\Delta t$, cho phép ghi nhận đầy đủ lực quán tính ly tâm, lực va đập khi đảo chiều và sự thay đổi góc truyền động theo thời gian thực.

3. Tại sao trong mô phỏng động học cơ cấu, việc định nghĩa Joint và Contact phải kết hợp chặt chẽ?

Joint (khớp nối) dùng để hạn chế các bậc tự do chuyển động tổng thể giữa các chi tiết (ví dụ: khớp quay Revolute chỉ giữ lại 1 DOF quay quanh trục Z). Tuy nhiên, Joint không tính được biến dạng cục bộ tại bề mặt chạm nhau. Contact (tiếp xúc phi tuyến) được áp đặt tại bề mặt làm việc (profin răng, rãnh cam) để truyền tải lực, tính toán ma sát và sinh ra trường ứng suất cục bộ chính xác.

4. Làm thế nào để xử lý lỗi không hội tụ (Non-convergence) trong bộ giải phi tuyến của ANSYS?

Khi gặp lỗi không hội tụ, kỹ sư cần: (1) Giảm kích thước bước thời gian (Initial/Minimum Substeps tăng lên); (2) Bật tính năng Large Deflection trong Analysis Settings; (3) Tinh chỉnh độ cứng tiếp xúc Normal Stiffness Factor xuống mức $0.1 - 0.5$; (4) Sử dụng công thức tiếp xúc Augmented Lagrange hoặc MPC; (5) Kiểm tra kỹ chất lượng lưới, loại bỏ các phần tử có độ méo (Skewness) cao tại vùng tiếp xúc.

5. Độ chính xác của kết quả mô phỏng động học trong đồ án được kiểm chứng bằng cách nào?

Kết quả mô phỏng được kiểm chứng chéo thông qua: (1) So sánh chuyển vị hình học lý thuyết của cơ cấu (ví dụ hành trình con trượt $279.88\text{ mm}$ khớp với tính toán hình học $280\text{ mm}$); (2) So sánh ứng suất uốn và ứng suất tiếp xúc cực đại với kết quả tính toán giải tích theo tiêu chuẩn ISO/AGMA; (3) Tiến hành thử nghiệm kiểm tra tính độc lập của lưới (Mesh Convergence Study) để đảm bảo kết quả không bị biến động khi tăng số lượng phần tử.

Kết luận

Đồ án "Mô phỏng động học các cơ cấu cơ khí bằng phần mềm ANSYS Workbench R16" đã hoàn thành toàn diện và xuất sắc các mục tiêu nghiên cứu đề ra. Thông qua việc làm chủ công cụ CAE hiện đại, đề tài đã số hóa và giải quyết thành công bài toán tiếp xúc phi tuyến, phân tích trường ứng suất, chuyển vị và biến dạng động lực học trên 8 hệ cơ cấu cơ khí trọng điểm của ngành chế tạo máy.

Công trình không chỉ mang giá trị học thuật sâu sắc trong việc trực quan hóa các quy luật cơ học phức tạp mà còn mở ra hướng ứng dụng thực tiễn to lớn cho các kỹ sư công nghiệp trong việc cắt giảm chi phí thử nghiệm, rút ngắn thời gian R&D và nâng cao độ tin cậy của sản phẩm cơ khí chính xác. Đây chính là nền tảng vững chắc để tiếp tục mở rộng nghiên cứu sang các lĩnh vực mô phỏng đa trường vật lý nâng cao, tối ưu hóa cấu trúc liên kết và dự báo tuổi thọ mỏi của hệ thống cơ điện tử trong tương lai.