Thử tìm kiếm một thuật toán thiết kế các mặt cong

Chia sẻ: Danh Tuong Vi | Ngày: | Loại File: PDF | Số trang:6

Thêm vào BST

Báo xấu

29
lượt xem 2
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Trong bài viết đặt vấn đề tìm kiếm một thuật toán dùng vào mục đích thiết kế các mặt cong có thể cần thiết trong ứng dụng công nghệ cao CNC. Các kết quả bước đầu nhận được trên cơ sở phát triển nghiên cứu ứng dụng về hàm hoá bề mặt lý thuyết vỏ tàu thuỷ của chính Tác giả. Các ý kiến thảo luận làm rõ triển vọng áp dụng của ý tưởng đặt ra.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Thử tìm kiếm một thuật toán thiết kế các mặt cong

Tạp chí Khoa học – Công nghệ Thủy sản số 01/2007 Trường Đại học Nha Trang VẤN ĐỀ TRAO ĐỔI THỬ TÌM KIẾM MỘT THUẬT TOÁN THIẾT KẾ CÁC MẶT CONG PGS.TS Nguyễn Quang Minh Khoa Cơ khí - Trường ĐH Nha Trang Trong bài báo đặt vấn đề tìm kiếm một thuật toán dùng vào mục đích thiết kế các mặt cong có thể cần thiết trong ứng dụng công nghệ cao CNC. Các kết quả bước đầu nhận được trên cơ sở phát triển nghiên cứu ứng dụng về hàm hoá bề mặt lý thuyết vỏ tàu thuỷ của chính Tác giả. Các ý kiến thảo luận làm rõ triển vọng áp dụng của ý tưởng đặt ra. 1. TỔNG QUAN Bài toán thiết kế, thay vì xấp xỉ trong nghĩa trực tiếp một mặt cong đặt ra ở đây có ý nghĩa thực tế và cần thiết, thậm chí trong nhiều trường hợp cần đến như phương tiện quan trọng, chẳng hạn trong các lĩnh vực công nghệ cao CNC. Một mặt cong, một cách chung nhất, có thể xem như tập hợp các đường cong phẳng - tiết diện của chính mặt cong đó với các mặt phẳng (P) được lựa chọn phù hợp. Nếu vậy, một cách tự nhiên có thể nghĩ rằng thuật toán thiết kế mặt cong tốt nhất nên đưa từ mô hình toán không gian về mô hình biểu diễn toán học chính xác các đường cong phẳng. Đơn giản là vì một khi bài toán thiết kế một đường cong cho trên mặt phẳng (P) bất kỳ đã được giải quyết, toàn bộ mặt cong cho trước đương nhiên có thể nhận được bằng cách tịnh tiến hoặc quay hợp lý các mặt phẳng (P) đó. Theo cách tiếp cận trực tiếp như vậy, bài toán đặt ra về thuật toán thiết kế một mặt cong có thể giải quyết không mấy khó khăn, ứng dụng mô hình toán xấp xỉ quen thuộc, đó là mô hình bài toán điều kiện biên, cách giải quyết được thể hiện qua các bước gồm: - Chọn dạng hàm cơ sở 52 - Áp dụng mọi điều kiện biên và xác lập số lượng thích hợp các phương trình căn cứ theo các điều kiện đó - Giải hệ phương trình điều kiện biên và biện luận các kết quả Để có thể đạt các kết quả xấp xỉ mong muốn cần có số các điểm thuộc đường cong giữ vai trò các điều kiện biên đủ lớn. Hàm cơ sở thông thường được chọn dưới dạng đa thức luỹ thừa bậc n viết tổng quát như dưới đây: k =n x = ∑ ak y k (1) k =0 và dạng chi tiết hoá: X =a0 +a1Y +a2Y2 +a3Y3 +a4Y4 +...+an−1Yn−1 +anYn (2) Với X,Y- là các toạ độ của các điểm cho trong hệ toạ độ lựa chọn XOY, k = 0, 1, 2, 3…, n-1, n. Dễ nhận xét rằng biểu thức (2) cùng lắm chỉ có thể dùng để tính toán gần đúng một hình cong mà ít hiệu quả trong các mục đích thiết kế đường cong đó. Phương pháp spline do Alberg. J. đề nghị nửa thế kỷ trước đây, được nghiên cứu áp dụng rộng rãi trên toàn thế giới nhằm khắc phục một phần nhược điểm của đa thức luỹ thừa (2) . Tuân theo nguyên tắc cơ bản đó là Tạp chí Khoa học – Công nghệ Thủy sản số 01/2007 đường cong hàm hoá được chia làm nhiều đoạn ngắn, mỗi đoạn có thể xấp xỉ theo các hàm đơn giản, phổ biến nhất thường chọn parabol bậc 3 viết tổng quát như dưới đây: k '= 3 x j = ∑ c j ',k ' ( y j − y j '−1 ) k ' (4) k '= 0 Trường Đại học Nha Trang Trong bài báo này chúng tôi giới thiệu một số kết quả phát triển từ nghiên cứu biểu diễn toán học bề mặt tàu thuỷ, vốn cũng là một vật thể có những đặc điểm riêng, hy vọng có thể áp dụng như giải thuật lập trình thiết kế và chế tạo chính xác một mặt cong theo từng mục đích cụ thể. với j = (j'-1), j', (j'+1) J' = 1, 2, 3,…, k-1 Nhờ những lợi thế cơ bản của nó trong lập trình, phương pháp spline đang được ứng dụng trong nhiều lĩnh vực khoa học công nghệ và đem lại những thành tựu quan trọng. Mặc dù vậy, phương pháp spline xấp xỉ toán học các đường cong phẳng cũng không thể được đánh giá như một công cụ hiệu quả nhất trong các mục đích thiết kế nghiêm túc. 2. MỘT SỐ KẾT QUẢ MỚI NGHIÊN CỨU XẤP XỈ CÁC ĐƯỜNG CONG PHẲNG Trên hình 1 minh hoạ mô hình toán lớp các đường cong phẳng 1a, 1b, 1c nằm trên mặt phẳng P trong quan hệ khác nhau đối với hệ toạ độ lựa chọn XOY. Tất cả các đường cong như vậy được cho bằng tập hợp các điểm rời rạc (Xi , Yi), i = 0, 1, 2, 3, 4, 5, …, n-1, n. O Hình 1 Sơ đồ các đường cong phẳng - đối tượng xấp xỉ toán học Về phương pháp, yêu cầu thiết kế các đường cong khác với chỉ đơn giản vẽ lại một đường cong cho trước là ở chỗ đòi hỏi, không chỉ áp dụng các điều kiện biên trong nghĩa thông thường, mà quan trọng hơn phải khai thác và áp dụng các điều kiện biên một cách hợp lý nhất, đảm bảo không chỉ tính phù hợp đối với từng trường hợp cụ thể, mà phải tiện lợi áp dụng đối với mọi đường cong thuộc lớp đang xét - đối tượng của bài toán xấp xỉ. Yêu cầu nói trên chỉ có thể là hiện thực nhờ việc sử dụng được các đặc điểm hình học quan trọng nhất được gọi là các yếu tố điều khiển. 53 Tạp chí Khoa học – Công nghệ Thủy sản số 01/2007 Có thể chứng tỏ rằng tồn tại 4 yếu tố có chức năng điều khiển như vậy, được xem như tập hợp các điều kiện cần để thông qua đó thiết kế chính xác các đường cong, căn cứ vào các yêu cầu thiết kế khác nhau. Các yếu tố như vậy gồm có: • Toạ độ điểm đầu (X0, Y0) • Toạ độ điểm cuối (Xt, Yt) • Diện tích hình cong Dt được tạo bởi đường cong cho trước với trục Toạ độ X = 0, X = 0, và các đường thẳng Y = Yt, và X = Xt Trường Đại học Nha Trang cho phép chọn ra từ tập hợp các đường cong, nhận được trên cơ sở 3 điều kiện biên đầu tiên, chính xác đường cong phải tìm như một nghiệm duy nhất của bài toán thiết kế. Trên cơ sở các biểu thức (3), (5.1) ÷ (5.4) có thể dễ dàng thiết lập hệ 4 phương trình mô tả 4 điều kiện biên giữ vai trò như các yếu tố điều khiển, nghiệm của hệ các phương trình đó cho các hệ số ai thuộc biểu thức luỹ thừa bậc 3: X = X 0 + a1Y + a 2Y 2 + a 3Y 3 được xác định theo các biểu thức sau đây: a1 = • Momen tĩnh Mox của diện tích nói trên ứng với trục hoành Thật vậy, giả sử kết quả xấp xỉ toán học chính xác một đường cong cho tương ứng như trên hình 1 tìm được dưới dạng một hàm liên tục: X = f (Y ) f (Y0 ) = X 0 (5.1) f ( X t ) = Yt (5.2) yt y0 y0 ∫ Xdy = ∫ f (Y )dy = Dt yt và : ∫ XYdy = y0 (5.3) yt ∫ f (Y )Ydy = M ox (5.4) y0 Hiển nhiên các điều kiện (5.1) và (5.2) phải được đảm bảo đầu tiên. Điều kiện (5.3) đảm bảo để các đường cong hàm hoá được, tương ứng với các điểm đầu và cuối đã cho, tạo thành với trục tọa độ các hình cong có cùng một diện tích Dt. Điều kiện cuối cùng (5.4) chứa một ưu thế thiết kế đặc biệt rất cần được khai thác đó là 54 a2 = − 1 ( X t + 32 Dt ) 7h (7) 4 × (8 Dt − 15Mx + X t ) 7h 2 (8) 20 (6 Dt − 12 Mx + X t ) 7h 3 (9) và a 3 = trong đó : Khi đó các biểu thức viết theo các điều kiện biên như mô tả ở trên có thể viết dễ dàng như dưới đây: yt (6) Dt = D − y0h Mx− y0h2 / 2 ; Mx= ; h = yt − y0 (10) h h2 Biểu thức (6) cùng với các biểu thức (7) ÷ (10) tỏ ra rất đơn giản, có thể áp dụng để thiết kế các đường cong với các điều kiện được cho theo (4). Tuy nhiên những gì nói ở trên mới chỉ bao gồm các điều kiện cần tối thiểu. Để cho bổ sung các điều kiện “đủ”, đòi hỏi phải xác định các đặc điểm của hình cong giới hạn bởi đường cong - đối tượng hàm hoá, cho phép phân biệt nó với mọi đương cong cùng lớp còn lại, chẳng hạn đó là tính liên tục, tính biến thiên, góc tạo bởi các tiếp tuyến với đường cong vẽ tại nhiều điểm khác nhau, đặc biệt tại hai điểm đầu và cuối. Sau cùng và có lẽ cũng là quan trọng nhất đó là đường cong có bao nhiêu điểm uốn trong miền xác định [y0,yt]. Chẳng hạn có thể nhận xét ngay rằng biểu thức (6) chỉ làm việc khi đường cong đối tượng hàm hoá cho trước chỉ có 1 điểm uốn, vì đạo Tạp chí Khoa học – Công nghệ Thủy sản số 01/2007 hàm bậc 2 chỉ nhận giá trị bằng 0 tại một điểm duy nhất tại y ∗ = − a 2 / 3a 3 Trong một số nghiên cứu ứng dụng có thể rất có lợi chọn hàm cơ sở dưới dạng hàm (1) x = k=n ∑ k =0 Trường Đại học Nha Trang mở rộng, trong đó thay vì các giá trị nguyên các luỹ thừa có thể nhận các số thực dương bất kỳ, viết dưới dạng tổng quát: a k y mk ; m > 0 ; k = 1, 2 , 3 , 4 ,..., n (11) hoặc chi tiết: X = X 0 + a1Y m + a 2Y 2 m + a 3Y 3m + a 4Y 4 m + ... + a k −1Y k −1 + a k Y k Áp dụng chính các điều kiện biên đã xét ở trên theo (4), sau đó tiếp tục thực hiện các bước tính toán như đã trình bày, có thể dẫn đến đa thức luỹ thừa bậc 2m với m thực dương như dưới đây: X = X 0 + a1Y m + a 2Y 2 m (13) Thừa số luỹ thừa m và các hệ số có mặt trong biểu thức được ai có thể xác định theo các biểu thức − 1,5( A − 2 B ) ± 2,25( A − 2 B ) 2 − 2( A − B )( A − 4 B + Xt ) m= 2( A − B ) a1 = (12) ( m + 1)[( 2 m + 1) A + X 0 − X t ] mh m (14) (15) X t − X 0 − a1 h m a2 = h 2m (16) trong đó, ngoài các ký hiệu đã được chú thích ở trên các ký hiệu mới được dùng gồm: A= ω tt − X 0 h h (17) và B = m ω oytt − X 0 h 2 / 2 h2 Biểu thức (13) cùng với các biểu thức (14)÷(18) có thể sử dụng tiện lợi và hiệu quả hơn rất nhiều so với biểu thức (6), đặc biệt khi cần thiết giải ngược phương trình, tính Y theo X. Trong các trường hợp như vậy chỉ cần qua một vài phép biến đổi dẫn về dạng phương trình bậc 2 sử dụng các nghiệm giải sẵn quen thuộc. (18) Vì cố gắng tìm kiếm lời giải trên cơ sở các nghiệm của phương trình bậc 2, hiển nhiên rằng biểu thức xấp xỉ (13) chỉ làm việc khi m có nghiệm thực, điều đó có nghĩa biểu thức dưới dấu căn, viết ở vế phải của biểu thức (14) không âm, thực chất được hiểu như bộ phận đặc biệt trong số những gì được gọi là điều kiện đủ để xấp xỉ toán học các đường cong thuộc lớp đang xét, được viết như sau: 55 Tạp chí Khoa học – Công nghệ Thủy sản số 01/2007 Trường Đại học Nha Trang 2 , 25 ( A − 2 B ) 2 − 2 ( A − B )( A − 4 B + X t ) >= 0 (19) Biểu thức (19) cho phép nhận đựơc các điều kiện đủ để xấp xỉ toán các đường cong như minh hoạ trên hình 1 bằng biểu thức dưới dạng đa thức luỹ thừa bậc 2m bao gồm: nếu: B > 2xt/ 3 (20) hoặc nếu: B< − ( B − 2 xt ) − 4 xt − 6 B 2X t ;A< 3 0,5 (20') và: B< − ( B − 2 xt ) + 4 xt − 6 B 2X t ;A> 3 0,5 (20'') Giả sử bây giờ nâng bậc đa thức luỹ thừa (13) lên bậc 3m bằng cách áp dụng thêm một điều kiện biên, chẳng hạn đó là hệ số góc tạo bởi tiếp tuyến với đường cong hàm hoá tại điểm cuối Y = Yt , Khi đó biểu thức xấp xỉ đối tượng tìm kiếm sẽ viết được dưới dạng: X = X 0 + a1 y m + a 2 y 2 m + a 3 y 3m (21) Trong đó thừa số luỹ thừa m có thể xác định như nghiệm của phương trình bậc 3 đủ: m 3 + 1,8333 ( A − 8B − 3 X 0 + 3xt − kh) ( A − 2 B) 2 ( A − 4 B − X 0 + X t ) m + m + 0,1667 =0 ( A − B) ( A − B) ( A − B) (22) Các hệ số ai có mặt trong biểu thức (21) được xác định như sau: a3 = m(3m + 1)( kh + X 0 − X t ) + (2m + 1)(3m + 1)[(m + 1) A + X 0 − X t ] m = 2m 2 h 3 m (23) kh m(3m + 2)( + X 0 − X t ) + (2m + 2)(3m + 2)[(m + 2) B + X 0 − X t ] m 2m 2 h 3 m kh + X 0 − X t − 2a 3 h 3 m (2m + 2)(3m + 2)[(m + 2) B + X 0 − X t ] a2 = m =− + 2m h m(3m + 2)h 2 m 2m(2m + 2)a3 h m(3m + 2)h 3m 3m và: =− (2m + 1)(3m + 1)[(m + 1) A + X 0 − X t ] + 2m(2m + 1)a3 h m(3m + 1)h 2 m X t − X 0 − a 2 h 2 m − a3 h 3m a1 = hm Hệ số góc tạo bởi tiếp tuyến tại điểm cuối (Xt,Yt) viết trong các biểu thức (22)÷(25) cũng mang ý nghĩa của các điều kiện “đủ” viết trong 56 (24) 3m (25) các biểu thức (20), (20') và (20'') áp dụng cho biểu thức (13).