1. Lý do chọn đề tài MỞ ĐẦU Công nghệ điều khiển số trên máy công cụ (Computer Numerical Control – viết tắt là CNC) được ra đời từ khoảng những năm 1940 [1] đã đánh dấu một bước nhảy vọt trong công nghiệp sản xuất. Việc ứng dụng máy CNC không những mang lại năng suất cao, chất lượng sản phẩm ổn định mà còn giải phóng được sức lao động của con người. Với việc các ngành khoa học điều khiển, công nghệ thông tin phát triển… các máy CNC càng trở lên mạnh mẽ hơn, điều đó cũng đồng nghĩa với việc các sản phẩm được chế tạo ra ngày càng phức tạp và tinh xảo. Do những ưu điểm của các máy CNC khiến cho càng ngày càng có nhiều nghiên cứu nhằm nâng cấp máy CNC trở lên ngày càng hoàn thiện. Có hai hướng nghiên cứu chính để cải thiện các máy CNC là hướng nghiên cứu cải thiện phần cứng và hướng nghiên cứu cải thiện phần mềm. Hướng nghiên cứu cải thiện phần cứng để nâng cao năng suất khi gia công trên máy phay CNC như: Nghiên cứu về kết cấu máy [2], [3]… Nghiên cứu ảnh hưởng của rung động [4]–[6]… Nghiên cứu về dụng cụ cắt khi gia công [7], [8]… Nghiên cứu cải thiện các vi mạch điều khiển để cải thiện khả năng điều khiển của máy CNC… Kết quả của các công trình khoa học cũng như thực nghiệm theo hướng này đã đạt được nhiều thành công. Các máy CNC càng ngày càng có kết cấu hoàn thiện, khả năng công nghệ được nâng cao. Qua đó năng suất cũng như chất lượng sản phẩm khi gia công trên các máy CNC có những cải thiện đáng kể góp phần không nhỏ vào sự phát triển của ngành công nghiệp gia công nói chung và gia công trên máy tự động nói riêng.
BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI Hoàng Văn Quý NGHIÊN CỨU LỰA CHỌN DỤNG CỤ VÀ ĐƯỜNG DỤNG CỤ TRONG TẠO HÌNH BỀ MẶT TỰ DO TRÊN MÁY PHAY CNC TRỤC Ngành: Kỹ thuật Cơ khí Mã số: 9520103 LUẬN ÁN TIẾN SĨ KỸ THUẬT CƠ KHÍ NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS Bùi Ngọc Tuyên Hà Nội – 2019 LỜI CAM ĐOAN MỤC LỤC DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT vii Danh mục ký hiệu vii Danh mục chữ viết tắt viii DANH MỤC CÁC BẢNG x DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ xii MỞ ĐẦU CHƯƠNG TỔNG QUAN VỀ GIA CÔNG MẶT TỰ DO 1.1 Mặt tự ứng dụng mặt tự 1.1.1 Giới thiệu 1.1.2 Lịch sử phát triển 1.1.3 Ứng dụng đường, mặt tự 1.1.3.1 Ứng dụng đường, mặt tự thiết kế 1.1.3.2 Ứng dụng đường, mặt tự gia công 1.2 Q trình gia cơng tạo hình bề mặt tự 11 1.3 Dụng cụ đường dụng cụ gia công mặt tự .13 1.3.1 Kiểu dụng cụ sử dụng gia công mặt tự 13 1.3.1.1 Dao phay ngón đầu phẳng 14 1.3.1.2 Dao phay ngón đầu cầu 15 1.3.1.3 Dao phay ngón đầu phẳng có góc lượn .15 1.3.1.4 Dao phay ngón thân 15 1.3.2 Đường dụng cụ gia công mặt tự .18 1.3.2.1 Khái niệm .18 1.3.2.2 Các kiểu đường dụng cụ 18 1.3.2.3 Các thông số đường dụng cụ 24 1.4 Phương pháp tính tốn sinh đường dụng cụ 25 1.4.1 Phương pháp sinh đường dụng cụ sở điểm tiếp xúc (CCBased) [29] 26 1.4.2 Phương pháp sinh đường dụng cụ sở điểm định vị dụng cụ (CL-Based) [29] .27 1.5 Tổng quan tình hình nghiên cứu gia công mặt tự máy phay CNC 27 1.5.1 Tổng quan tình hình nghiên cứu nước 27 1.5.2 Tổng quan tình hình nghiên cứu giới 29 1.6 Đánh giá tình hình nghiên cứu phương pháp gia công mặt tự đề xuất hướng nghiên cứu luận án .31 1.7 Kết luận 33 CHƯƠNG PHƯƠNG PHÁP BIỂU DIỄN ĐƯỜNG VÀ MẶT TỰ DO TRONG CÁC HỆ THỐNG CAD/CAM 34 2.1 Biểu diễn đường tự hệ thống CAD/CAM 35 iv 2.1.1 Các phương pháp biểu diễn toán đường tổng quát 35 2.1.2 Biểu diễn đường cong Bezier 36 2.1.3 Biểu diễn đường cong B-spline 40 2.2 Biểu diễn mặt tự hệ thống CAD/CAM 43 2.2.1 Phương pháp biểu diễn mặt dạng mơ hình đa thức dạng tham số 44 2.2.1.1 Mơ hình tốn biểu diễn mảnh mặt Ferguson .44 2.2.1.2 Mơ hình mảnh mặt Bezier 45 2.2.2 Phương pháp biểu diễn mặt dạng mơ hình nội suy ranh giới mảnh mặt 46 2.2.2.1 Biểu diễn mặt tự mơ hình mặt kẻ 46 2.2.2.2 Biểu diễn mặt tự mơ hình mặt Coon 46 2.2.3 Phương pháp biểu diễn mặt dạng mơ hình mảnh qt 47 2.2.3.1 Mảnh mặt trượt tịnh tiến 47 2.2.3.2 Mảnh mặt quay 47 2.2.4 Phương pháp biểu diễn mặt dạng tứ giác 48 2.2.5 Mảnh mặt B - spline 48 2.2.6 Mặt NURBS 49 2.2.6.1 Biểu diễn toán học bề mặt NURBS 49 2.3.6.2 Đạo hàm mặt NURBS 51 2.2.7 Công cụ mô số mặt sang dạng mặt tự 52 2.3.8 Cấu trúc tệp (file) biểu diễn mặt theo định dạng IGES CAD/CAM 53 2.2.8.1 Cấu trúc file IGES 54 2.2.8.2 Biểu diễn số đối tượng theo định dạng file IGES 58 2.3 Kết luận 61 CHƯƠNG XÂY DỰNG PHƯƠNG PHÁP LỰA CHỌN DỤNG CỤ TRONG GIA CƠNG TẠO HÌNH BỀ MẶT TỰ DO TRÊN MÁY PHAY CNC TRỤC 62 3.1 Giới thiệu 62 3.2 Ảnh hưởng dụng cụ tạo hình bề mặt máy phay CNC 62 3.3 Xây dựng phương án lựa chọn dụng cụ cắt hợp lý để gia cơng tạo hình mặt tự 66 3.3.1 Phân vùng mặt tự 67 3.3.2 Lựa chọn dụng cụ hợp lý tương ứng với phân vùng 71 3.3.3 Xác định ranh giới tương ứng với phân vùng cục 72 3.3.4 Thí nghiệm đánh giá kết .80 3.4 Đề xuất phương án gia công với bề mặt phân vùng: 90 3.4 Kết luận 92 CHƯƠNG XÂY DỰNG PHƯƠNG PHÁP LỰA CHỌN ĐƯỜNG DỤNG CỤ TRONG GIA CƠNG TẠO HÌNH BỀ MẶT TỰ DO TRÊN MÁY PHAY CNC TRỤC 94 4.1 Xây dựng thực nghiệm Taguchi đánh giá ảnh hưởng đường dụng cụ tới chất lượng tạo hình bề mặt tự 95 v 4.1 Phương pháp Taguchi .96 4.1.2 Xây dựng thực nghiệm .97 4.1.2.1 Điều kiện thực nghiệm 97 4.1.2.2 Lựa chọn thông số đầu vào xây dựng mảng trực giao Taguchi 99 4.2 Đề xuất phương án xác định đường dụng cụ 103 4.2.1 Điều kiện biên 104 4.2.2 Tính tốn thông số đường dụng cụ 104 4.2.2.1 Tính tốn bước tiến ngang (St).(Bổ đề 1) 105 4.2.2.2 Xác định đường dẫn hợp lý đảm bảo Sc ≤ [Sc] (Bổ đề 2) 108 4.2.3 Tính tốn sai số xấp xỉ 112 4.3 Mô thực nghiệm đánh giá kết 114 4.4 Kết luận 117 KẾT LUẬN VÀ KIẾN NGHỊ 118 DANH MỤC TÀI LIỆU THAM KHẢO 120 PHỤ LỤC 127 vi DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT Danh mục ký hiệu STT Ký hiệu Giải nghĩa C(u) u v Bi, n(u) Ri, n(u) Pi x y z 10 Pi, j 11 12 13 14 U Wi, j i, j 15 Sc 16 17 18 19 20 21 22 23 24 25 K H u(t) v(t) t D M Q C Phương trình biểu diễn đường cong tự với biến u Tham số đường cong tự theo phương u (0 ≤ u ≤ 1) Tham số đường cong tự theo phương v (0 ≤ v ≤ 1) Đa thức Bernstein với biến u Hàm vô hướng sở cho đường cong tự Điểm điều khiển thứ i đường cong mặt cong tự Biến số hàm số theo phương Ox hệ tọa độ Descartes Biến số hàm số theo phương Oy hệ tọa độ Descartes Biến số hàm số theo phương Oz hệ tọa độ Descartes Điểm điều khiển nút (i, j) lưới điểm điều khiển mặt tự Véc tơ nút Trọng số tương ứng điểm điều khiển Pi, j mặt tự Chỉ số Bán kính cong Scallop height (chiều cao lượng kim loại để lại sau bước dịch dao ngang) Độ cong Độ cong Gauss Đường cong tham số t theo phương u Đường cong tham số t theo phương v Tham số t đường cong Ma trận thứ hai đường cong Ma trận chuyển Toán tử sai phân Ma trận hệ số góc Ma trận hệ số Ferguson 26 d du 27 28 29 30 m mm µm S(u,v) Phép tính đạo hàm với biến u Đơn vị đo độ dài: mét Đơn vị đo độ dài: mi li mét (10-3m) Đơn vị đo độ dài: micro mét (10-6m) Mặt tự biểu diễn không gian tham số vii 31 32 33 34 TT St 35 Tn 36 37 38 39 G Si Sc Tpi T n Vector hướng dụng cụ Vector pháp tuyến mặt S(u,v) Mặt phẳng tiếp tuyến mặt cong S(u,v) Khoảng dịch dao ngang (Step over size) Mặt phẳng tạo vector hướng dao vector pháp tuyến Giao tuyến mặt phẳng Tn TT Điểm i mặt S(u,v) Đường cong giao mặt Tn S(u,v) Đường dụng cụ thứ i Danh mục chữ viết tắt STT Viết tắt CAD CAM CAGD CIM Viết đầy đủ Computer Aided Design Computer Aided Manufacturing Computer Aided Geometry Design Computer Intergrated Manufacturing viii Giải nghĩa Thiết kế có hỗ trợ máy tính Sản xuất có hỗ trợ máy tính Thiết kế hình học có hỗ trợ máy tính Sản xuất có tích hợp máy tính CAE Computer Aided Engineering CAPP Computer Aided Production Planning NC Numerical Control CNC Computer Numerical Control 10 NURBS Non-uniform rational Bspline 11 B-spline B-spline 12 2D Two Dimension 13 3D Three Dimension 14 5D Five Dimension 15 16 CC point CL point Cutting Contact point Cutting Location point 17 CC path Cutting contact path 18 CL path Cutting Location path 19 DP DPi, i=1, Data Points 20 set of data points ix Ứng dụng máy tính phân tích Lập kế hoạch chế tạo sản phẩm có trợ giúp máy tính Điều khiển số Điều khiển số máy tính (chỉ hệ thống khí gia cơng tự động Mặt đường B-spline hữu tỉ không đồng Đường cong tự Bspline Không gian phẳng (dùng máy trục) Không gian ba chiều (dùng máy trục) Không gian năm chiều (dùng máy trục) Điểm cắt (điểm tạo hình) Điểm định vị dụng cụ Đường tập hợp (quỹ đạo) điểm tạo hình Đường tập hợp (quỹ đạo) điểm định vị dụng cụ Tập điểm liệu tập liệu điểm thành phần DANH MỤC CÁC BẢNG Bảng 1 Các đường dụng cụ 18 Bảng Thuật ngữ đường dụng cụ gia công CNC 19 Bảng Mười dạng vùng bề mặt cục bề mặt tự trơn liên tục 29 Bảng Bảng cấu trúc liệu file IGES 55 Bảng 2 Bảng thống kê nhãn biểu diễn liệu IGES 57 Bảng Bảng điểm điều khiển trọng số mặt tự mẫu thí nghiệm 80 Bảng Bảng tính chất vật lý nhựa PA 81 Bảng 3 Bảng thông số dụng cụ 82 Bảng Bảng thông số máy đo tọa độ SVANEX 9106 82 Bảng Bảng thông số cao độ z điểm bề mặt mẫu thiết kế 84 Bảng Bảng thông số đo cao độ z điểm bề mặt Mẫu .84 Bảng Bảng thông số đo cao độ z điểm bề mặt Mẫu .84 Bảng Chương trình gia công 88 Bảng Bảng tổng hợp thời gian gia công chiều dài đường dụng cụ 89 x Bảng 10 Bảng tông hợp kết đánh giá 90 Bảng 11 Bảng so sánh phương pháp gia công với bề mặt phân vùng 91 Bảng 12 Bảng thống kê gia công mẫu theo ba phương án 91 Bảng Bảng thông số dụng cụ 98 Bảng Bảng thông số máy đo tọa độ SVA NEX9016 98 Bảng Bảng tọa độ điểm thiết kế mặt phôi 99 Bảng 4 Bảng thông số F, S, T 101 Bảng Mảng trực giao Taguchi L9 (3^3) 101 Bảng Kết đo vị trí xác định bề mặt mẫu 102 Bảng Phân tích ANOVA thơng số ảnh hưởng tới tạo hình bề mặt 103 Bảng Điểm điều khiển mảnh mặt lõm 112 xi DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ Hình 1 Hiệu chỉnh đường tròn Hình Thiết kế có sử dụng đường tự NURBS Hình Hiệu chỉnh mặt trụ thường mặt trụ biểu diễn dạng tự Hình Ứng dụng NURBS thiết kế Hình Mơ tả nội suy tuyến tính cung hệ máy CNC Hình Nội suy NURBS hệ CNC 10 Hình Nội suy tuyến tính nội suy NURBS gia cơng mặt tự 11 Hình Một ứng dụng mặt tự khuôn 12 Hình Khảo sát thời gian gia công thô, tinh, đánh bóng gia cơng khn 12 Hình 10 Mơ hình hình học dao phay ngón tổng qt 14 Hình 11 Các dạng dao phay ngón .14 Hình 12 Dụng cụ cắt gia công mặt tự .15 Hình 13 Sự khác biệt vùng lồi, vùng lõm vùng phẳng gia công 16 Hình 14 Đường cắt cắt lẹm 17 Hình 15 Gia cơng mặt tự máy trục 17 Hình 16 Đường dụng cụ .18 Hình 17 Một số thuật ngữ đường dụng cụ gia công máy phay CNC .20 Hình 18 Một số kiểu đường dụng cụ Catia V5R20 .23 Hình 19 Các thơng số quan trọng gia công mặt tự .25 xii % Truc cua duong cong AMB_Data.axis(1)=subplot('position',[.05,.07,.9,.77]); axis([-.01,1.01,-.01,1.01]); axis manual; axis on hold on; % Khung gioi han cua chuong trinh fill([-.1,1.1,1.1,-.1],[-.1,-.1,1.1,1.1],[1,1,1],'edgecolor','none'); plot([0,1,1,0,0],[0,0,1,1,0],'-k'); %%% Hàm Callback %%% function restart(varargin) global AMB_Data; AMB_Data.c=zeros(0,2); AMB_Data.movetype=0; AMB_Data.moveind=0; for k=AMB_Data.handles if k,delete(k);end end AMB_Data.handles=[0,0]; set(AMB_Data.uicontrols(1),'String','Bac: -1'); % Chèn, xóa, chỉnh sửa điểm điều khiển function myDownfct(varargin) global AMB_Data; if strcmp(get(gcf,'SelectionType'),'normal') % Chèn thêm điểm điều khiển, nâng bậc đường cong p=get(AMB_Data.axis(1),'CurrentPoint'); p=p(1,1:2); if min(p) >= & max(p) = & max(p) = & max(p) = & max(p) 0 % Ve cap nhat them da giac dieu khien if AMB_Data.handles(1)==0 AMB_Data.handles(1)=plot(AMB_Data.c(:,1),AMB_Data.c(:,2),'-ok'); 129 else set(AMB_Data.handles(1),'Xdata',AMB_Data.c(:,1)); set(AMB_Data.handles(1),'Ydata',AMB_Data.c(:,2)); end end if n>1 % ve cap nhat them duong t=linspace(0,1,panz); s=casteljau(AMB_Data.c,t'); if AMB_Data.handles(2)==0 AMB_Data.handles(2)=plot(s(:,1),s(:,2),'-b','Linewidth',2); else set(AMB_Data.handles(2),'Xdata',s(:,1)); set(AMB_Data.handles(2),'Ydata',s(:,2)); end end %%% Ham cua chuong trinh AMB_program %%% function p = casteljau(c,t) % sư dụng thuật toán Pierre De Casteljau % p = casteljau(c,t) % points p(j,:) Duong cong Bezier voi tap diem dieu khien c(k,:) % Cho thong so t(j) % Chuyen doi sang dang vecto cot (Khi can thiet) t = t(:); [np1,dim] = size(c); % Tinh de quy cua De casteljau if np1 == p = (1+0*t)*c; else p0 = casteljau(c(1:end-1,:),t); p1 = casteljau(c(2:end,:),t); t = repmat(t,1,dim); p = (1-t).*p0 + t.*p1; end Phụ lục 2.2 Biểu diễn mặt phẳng CATIA V5R20 theo định dạng IGES START RECORD GO HERE S 1H,,1H;,5HPart1,15HPlanarCATIA.igs,44HDASSAULT SYSTEMES CATIA V5 R20 - wG ww.3ds.com,27HCATIA Version Release 20 ,32,75,6,75,15,5HPart1,1.0,2,2HG MM,1000,1.0,15H20171104.110143,0.001,10000.0,5HHVQuy,8HHVQuy-PC,11,0,15HG 20171104.110143,; G 406 0 0 000010201D 406 0 15 0D 128 0 0 001010001D 128 0 13 0D 126 15 0 0 001010001D 126 0 0D 126 20 0 0 001010001D 126 0 0D 126 25 0 0 001010001D 126 0 0D 126 30 0 0 001010001D 126 0 0D 102 35 0 0 001010001D 102 0 0D 110 36 0 0 001010501D 110 0 0D 110 37 0 0 001010501D 130 1 4 10 11 12 13 14 15 16 17 110 0 0D 110 38 0 0 001010501D 110 0 0D 110 39 0 0 001010501D 110 0 0D 102 40 0 0 001010501D 102 0 0D 142 41 0 0 001010001D 142 0 0D 144 42 0 10000 0 000000000D 144 0 0D 406,1,9HSurface.1,0,0; 1P 128,3,3,3,3,0,0,1,0,0,0.0,0.0,0.0,0.0,197.8293762,197.8293762, 3P 197.8293762,197.8293762,0.0,0.0,0.0,0.0,131.7123485,131.7123485, 3P 131.7123485,131.7123485,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0, 3P 1.0,1.0,1.0,1.0,1.0,1.0,-7.645088975,5.987486329,0.0, 3P -73.58821438,5.987486329,0.0,-139.5313398,5.987486329,0.0, 3P -205.4744652,5.987486329,0.0,-7.645088975,49.89160251,0.0, 3P -73.58821438,49.89160251,0.0,-139.5313398,49.89160251,0.0, 3P -205.4744652,49.89160251,0.0,-7.645088975,93.79571869,0.0, 3P -73.58821438,93.79571869,0.0,-139.5313398,93.79571869,0.0, 3P -205.4744652,93.79571869,0.0,-7.645088975,137.6998349,0.0, 3P -73.58821438,137.6998349,0.0,-139.5313398,137.6998349,0.0, 3P -205.4744652,137.6998349,0.0,0.0,197.8293762,0.0,131.7123485,0, 3P 0; 3P 126,3,3,0,0,1,0,0.0,0.0,0.0,0.0,197.8293762,197.8293762, 5P 197.8293762,197.8293762,1.0,1.0,1.0,1.0,-7.645088975, 5P 5.987486329,0.0,-73.58821438,5.987486329,0.0,-139.5313398, 5P 5.987486329,0.0,-205.4744652,5.987486329,0.0,0.0,197.8293762, 5P 0.0,0.0,0.0,0,0; 5P 126,3,3,0,0,1,0,0.0,0.0,0.0,0.0,131.7123485,131.7123485, 7P 131.7123485,131.7123485,1.0,1.0,1.0,1.0,-205.4744652, 7P 5.987486329,0.0,-205.4744652,49.89160251,0.0,-205.4744652, 7P 93.79571869,0.0,-205.4744652,137.6998349,0.0,0.0,131.7123485, 7P 0.0,0.0,0.0,0,0; 7P 126,3,3,0,0,1,0,0.0,0.0,0.0,0.0,197.8293762,197.8293762, 9P 197.8293762,197.8293762,1.0,1.0,1.0,1.0,-205.4744652, 9P 137.6998349,0.0,-139.5313398,137.6998349,0.0,-73.58821438, 9P 137.6998349,0.0,-7.645088975,137.6998349,0.0,0.0,197.8293762, 9P 0.0,0.0,0.0,0,0; 9P 126,3,3,0,0,1,0,0.0,0.0,0.0,0.0,131.7123485,131.7123485, 11P 131.7123485,131.7123485,1.0,1.0,1.0,1.0,-7.645088975, 11P 137.6998349,0.0,-7.645088975,93.79571869,0.0,-7.645088975, 11P 49.89160251,0.0,-7.645088975,5.987486329,0.0,0.0,131.7123485, 11P 0.0,0.0,0.0,0,0; 11P 102,4,5,7,9,11,0,0; 13P 110,0.0,0.0,0.0,197.8293762,0.0,0.0,0,0; 15P 110,197.8293762,0.0,0.0,197.8293762,131.7123485,0.0,0,0; 17P 110,197.8293762,131.7123485,0.0,0.0,131.7123485,0.0,0,0; 19P 110,0.0,131.7123485,0.0,0.0,0.0,0.0,0,0; 21P 102,4,15,17,19,21,0,0; 23P 142,0,3,23,13,1,0,0; 25P 144,3,1,0,25,0,1,1; 27P S 1G 4D 28P 42 T 18 19 20 21 22 23 24 25 26 27 28 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 Phụ lục 2.3 Cấu trúc mặt trụ tâm O(0, 0), bán kính R=30, chiều cao h = 20 theo định dạng IGES CATIA V5R20 START RECORD GO HERE S 1H,,1H;,,16HCATIA MatTru.igs,44HDASSAULT SYSTEMES CATIA V5 R20 - www.3dsG com,27HCATIA Version Release 20 ,32,75,6,75,15,,1.0,2,2HMM,1000,1.0, G 15H20171104.142321,0.001,10000.0,5HHVQuy,8HHVQuy-PC,11,0,15H20171104.142G 321,; G 124 0 0 001020201D 124 0 0D 100 0 0 001010001D 100 0D 131 1 4 124 0 0 001020201D 124 0 0D 100 0 0 001010001D 100 0D 102 0 10000 0 000000000D 102 Sketch.1 0D 406 0 0 000010201D 406 0 15 0D 110 0 0 001010001D 110 0 0D 110 0 0 001010001D 110 0 0D 120 0 0 001010001D 120 0 0D 110 10 0 0 001010001D 110 0 0D 124 11 0 0 001020201D 124 0 0D 100 12 0 0 21 001010001D 100 0 0D 110 13 0 0 001010001D 110 0 0D 124 14 0 0 001020201D 124 0 0D 100 15 0 0 27 001010001D 100 0 0D 102 16 0 0 001010001D 102 0 0D 142 17 0 0 001010001D 142 0 0D 144 18 0 10000 0 000000000D 144 0 0D 110 19 0 0 001010001D 110 0 0D 110 20 0 0 001010001D 110 0 0D 120 21 0 0 001010001D 120 0 0D 110 22 0 0 001010001D 110 0 0D 124 23 0 0 001020201D 124 0 0D 100 24 0 0 45 001010001D 100 0 0D 110 25 0 0 001010001D 110 0 0D 124 26 0 0 001020201D 124 0 0D 100 27 0 0 51 001010001D 100 0 0D 102 28 0 0 001010001D 102 0 0D 142 29 0 0 001010001D 142 0 0D 144 30 0 10000 0 000000000D 144 0 0D 124,1.0,0.0,0.0,0.0,0.0,1.0,0.0,0.0,0.0,0.0,1.0,0.0,0,0; 1P 100,0.0,0.0,0.0,30.0,0.0,-30.0,0.0,0,0; 3P 124,-1.0,0.0,0.0,0.0,0.0,-1.0,0.0,0.0,0.0,0.0,1.0,0.0,0,0; 5P 100,0.0,0.0,0.0,30.0,0.0,-30.0,0.0,0,0; 7P 102,2,3,7,0,0; 9P 406,1,9HExtrude.1,0,0; 11P 110,0.0,0.0,10.0,0.0,0.0,11.0,0,0; 13P 110,30.0,0.0,0.0,30.0,0.0,20.0,0,0; 15P 120,13,15,0.0,6.283185307,0,0; 17P 110,30.0,0.0,20.0,30.0,0.0,0.0,0,0; 19P 124,1.0,0.0,0.0,0.0,0.0,1.0,0.0,0.0,0.0,0.0,1.0,0.0,0,0; 21P 100,0.0,0.0,0.0,30.0,0.0,-30.0,0.0,0,0; 23P 132 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 10 11 12 110,-30.0,0.0,0.0,-30.0,0.0,20.0,0,0; 124,-1.0,0.0,0.0,0.0,0.0,1.0,0.0,0.0,0.0,0.0,-1.0,20.0,0,0; 100,0.0,0.0,0.0,30.0,0.0,-30.0,0.0,0,0; 102,4,19,23,25,29,0,0; 142,0,17,0,31,2,0,0; 144,17,1,0,33,0,1,11; 110,0.0,0.0,10.0,0.0,0.0,11.0,0,0; 110,-30.0,0.0,0.0,-30.0,0.0,20.0,0,0; 120,37,39,0.0,6.283185307,0,0; 110,-30.0,0.0,20.0,-30.0,0.0,0.0,0,0; 124,-1.0,0.0,0.0,0.0,0.0,-1.0,0.0,0.0,0.0,0.0,1.0,0.0,0,0; 100,0.0,0.0,0.0,30.0,0.0,-30.0,0.0,0,0; 110,30.0,0.0,0.0,30.0,0.0,20.0,0,0; 124,1.0,0.0,0.0,0.0,0.0,-1.0,0.0,0.0,0.0,0.0,-1.0,20.0,0,0; 100,0.0,0.0,0.0,30.0,0.0,-30.0,0.0,0,0; 102,4,43,47,49,53,0,0; 142,0,41,0,55,2,0,0; 144,41,1,0,57,0,1,11; S 1G 4D 60P 30 25P 27P 29P 31P 33P 35P 37P 39P 41P 43P 45P 47P 49P 51P 53P 55P 57P 59P T 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 Phụ lục 3.1 Hàm vector nut % Ham tinh toan vecto nut function xi=KnotVector(k,N) xi(1:k)=0; for i=k+1:N xi(i)=(i-k)/(N-k+1); end xi(N+1:N+k)=1; Phụ lục 3.2 Hàm Ma trận sở % ham tao ma tran co so function MN = MatrixBasisFunction(N,knot,t,xknot) for i=1:length(t) for j=1:N MN(i,j)=BasisFunction(j,knot,t(i),xknot); end end Phụ lục 3.3 Hàm độ cong lớn, nhỏ vùng cục %Do cong lon nhat, nho nhat vung cuc bo function[lcPmax,lcPmin]=MaxMinCurvatureLocalSurface(M,gPmax,gPmin) Vmax=[]; Vmin=[]; for i=1:size(M,1) for j=1:size(M,2) if(M(i,j)~=0) Vmax=[Vmax,gPmax(i,j)]; Vmin=[Vmin,gPmin(i,j)]; end end end lcPmax=max(Vmax); lcPmin=min(Vmin); Phụ lục 3.4 Tính tốn điểm bề mặt % Tinh toan diem be mat bang cach su dung mang diem dieu khien P so w SurfacePoint (n, p, U, m, q, V,w, u, v, S) {/* Tinh toan diem tren mat */ /* Input: n,p,U,m,q,V,w,u,v */ /* Output: S */} uspan = FindSpan(n,p,u,U); 133 BasicFuns (uspan,u,p,U,Nu); vspan = FindSpan(m,q,v,V); BasicFuns (vspan,v,q,V,Nv); for (l=0, l