Báo cáo Thiết kế bộ điều khiển trượt cho tay máy Robot 2 bậc tự do và mô phỏng trên Matlab - Simulink

88 trang hapham 2580

Download

Bạn đang xem 20 trang mẫu của tài liệu "Báo cáo Thiết kế bộ điều khiển trượt cho tay máy Robot 2 bậc tự do và mô phỏng trên Matlab - Simulink", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

bao_cao_thiet_ke_bo_dieu_khien_truot_cho_tay_may_robot_2_bac.doc

Nội dung text: Báo cáo Thiết kế bộ điều khiển trượt cho tay máy Robot 2 bậc tự do và mô phỏng trên Matlab - Simulink

Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do TRƯỜNG . KHOA  BÁO CÁO TỐT NGHIỆP Đề tài: Thiết kế bộ Điều khiển trượt cho tay máy Robot 2 bậc tự do và mô phỏng trên Matlab – Simulink SV: Hoàng Đức Thắng - Phan Thanh Hải 1
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Mục lục Lời nói đầu 5 CHƯƠNG I: TỔNG QUAN VỀ ROBOT CÔNG NGHIỆP 6 I.1.Robot công nghiệp: 6 I.1.1. Sự ra đời của Robot công nghiệp : 6 I.1.2.Phân loại tay máy Robot công nghiệp: 7 I.2. Ứng dụng của Robot công nghiệp : 9 I.2.1.Mục tiêu ứng dụng Robot công nghiệp : 9 I.2.2.Các lĩnh vực ứng dụng Robot công nghiệp : 10 I.2.3. Các xu thế ứng dụng Robot trong tương lai : 10 I.2.4. Tình hình tiếp cận và ứng dụng Robot công nghiệp ở Việt Nam : 10 I.3.Cấu trúc của Robot công nghiệp: 11 I.3.1.Các bộ phận cấu thành Robot công nghiệp : 11 I.3.2.Bậc tự do và các toạ độ suy rộng : 12 I.3.2.1.Bậc tự do : 12 I.3.2.2. Toạ độ suy rộng : 13 I.3.3.Nhiệm vụ lập trình điều khiển Robot: 14 I.3.3.1. Định vị và định hướng tại “điểm tác động cuối” : 14 I.3.3.2. Lập trình điều khiển Robot công nghiệp : 14 SV: Hoàng Đức Thắng - Phan Thanh Hải 2
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do I.4. Các phép biến đổi toán học cho Robot : 15 I.4.1.Biến đổi toạ độ dùng Ma trận: 15 I.4.1.1. Vector điểm và toạ độ thuần nhất : 15 I.4.1.3.Biến đổi Ma trận dùng toạ độ thuần nhất: 18 I.4.1.4. Ý nghĩa hình học của Ma trận thuần nhất: 19 I.4.2.Các phép biến đổi cơ bản: 22 I.4.2.1.Phép biến đổi tịnh tiến: 22 I.4.2.2. Phép quay quanh các trục toạ độ : 22 CHƯƠNG II: HỆ PHƯƠNG TRÌNH ĐỘNG HỌC VÀ ĐỘNG LỰC HỌC CỦA ROBOT CÔNG NGHIỆP: 24 II.1. Hệ phương trình động học Robot : 24 II.1.1. Đặt vấn đề : 24 II.1.2. Xác định trạng thái của Robot tai điểm tác động cuối : 24 II.1.3. Mô hình động học : 25 II.1.3.1. Ma trận quan hệ : 25 II.1.3.2. Bộ thông số DH : 26 II.1.3.3. Thiết lập hệ toạ độ : 27 II.1.3.4. Mô hình biến đổi : 28 II.1.3.5. Phương trình động học : 29 II.2. Tổng hợp chuyển động Robot : 30 II.2.1. Nhiệm vụ : 30 II.2.2. Bài toán động học ngược : 30 II.2.3. Các phương pháp giải bài toán động học ngược : 31 II.3. Động lực học Robot: 31 II.3.1.Nhiệm vụ và phương pháp phân tích Động lực học Robot: 31 SV: Hoàng Đức Thắng - Phan Thanh Hải 3
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do II.3.2.Vận tốc và gia tốc: 33 II.3.3. Động năng tay máy: 36 II.3.4 .Thế năng tay máy: 38 II.3.5.Mô hình động lực học tay máy: 38 II.3.6. Động lực học của cơ cấu tay máy 2 khâu: 40 CHƯƠNG III : THIẾT KẾ BỘ ĐIỀU KHIỂN TRƯỢT CHO TAY MÁY ROBOT 2 BẬC TỰ DO 45 III.1.Hệ phi tuyến : 45 III.1.1.Hệ phi tuyến là gì ? 45 III.1.2.Mô hình trạng thái và quỹ đạo trạng thái của Hệ phi tuyến: 47 III.1.2.1.Mô hình trạng thái: 47 III.1.2.2.Quỹ đạo trạng thái : 49 III.1.3. Điểm cân bằng và điểm dừng của hệ thống: 50 III.1.3.1. Điểm cân bằng: 50 III.1.3.2.Điểm dừng của hệ : 50 III.1.3.3 Tính ổn định tại một điểm cân bằng: 51 III.1.4 Tiêu chuẩn ổn đinh Lyapunov : 52 III.1.4.1.Tiêu chuẩn Lyapunov: 53 III.1.4.2.Tiêu chuẩn Lyapunov phục vụ thiết kế bộ điều khiển: 56 III.2.Bậc tương đối của hệ phi tuyến: 56 III.3.Tính động hoc không: 59 III.4.Thiết kế bộ Điều khiển trượt cho tay máy: 62 III.4.1.Điều khiển trượt: 62 III.4.1.1.Trường hợp bậc tương đối của hệ bằng bậc của hệ p=n: 63 III.4.1.2. Trường hợp bậc tương đối của hệ p<n 64 SV: Hoàng Đức Thắng - Phan Thanh Hải 4
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do III.4.2. Thiết kế bộ điều khiển trượt cho tay máy n bậc tự do: 65 III.4.3. Ứng dụng Điều khiển trượt cho tay máy Robot 2 bậc tự do: 69 III.4.3.1. Phương trình động lực học tay máy hai bậc tự do toàn khớp quay: 69 III.4.3.2. Mô hình động lực học tay máy hai bậc tự do: 70 III.4.3.3. Thiết kế bộ điều khiển trượt cho tay máy 2 bậc tự do: 71 III.4.3.4. Tính toán giá trị đặt i cho tay máy hai bậc tự do: 73 CHƯƠNG IV: MÔ PHỎNG QUÁ TRÌNH CHUYỂN ĐỘNG CỦA ROBOT DÙNG BỘ ĐIỀU KHIỂN TRƯỢT TRÊN NỀN MATLAB AND SIMULINK: 75 IV.1. Tổng quan về Matlab-Simulink: 75 IV.2. Các thao tác thực hiện mô phỏng: 76 Kết luận 85 Tài liệu tham khảo : 85 SV: Hoàng Đức Thắng - Phan Thanh Hải 5
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Lời nói đầu Trong sự nghiệp công nghiệp hoá, hiện đại hoá đất nước vấn đề tự động hoá sản xuất có vai trò đặc biệt quan trọng. Mục tiêu ứng dụng kỹ thuật Robot trong công nghiệp nhằm nâng cao năng suất dây chuyền công nghệ, nâng cao chất lượng và khả năng cạnh tranh của sản phẩm, đồng thời cải thiện điều kiện lao động. Sự cạnh tranh hàng hoá đặt ra một vấn đề thời sự là làm sao để hệ thống tự động hoá sản xuất phải có tính linh hoạt nhằm đáp ứng với sự biến động thường xuyên của thị trường hàng hoá. Robot công nghiệp là bộ phận cấu thành không thể thiếu trong hệ thống sản xuất tự động linh hoạt đó. Gần nửa thế kỉ có mặt trong sản xuất. Robot công nghiệp đã có một lịch sử phát triển hấp dẫn. Ngày nay, Robot công nghiệp được dùng rộng rãi ở nhiều lĩnh vực sản suất. Điều đó xuất phát từ những ưu điểm cơ bản của các loại Robot đã được lựa chọn và đúc kết qua bao nhiêu năm ứng dụng ở nhiều nước. Ở nước ta, trước những năm 1990 hầu như chưa du nhập về kỹ thuật Robot. Từ năm 1990 nhiều cơ sở công nghiệp đã bắt đầu nhập ngoại nhiều loại Robot phục vụ các việc như tháo lắp dụng cụ cho các trung tâm CNC, lắp ráp các linh kiện điện tử, hàn vỏ xe ô tô, xe máy và phun phủ bề mặt Có những nơi đã bắt đầu thiết kế chế tạo và lắp ráp Robot. Có thể nói, Robot đã và đang góp phần rất lớn vào sự nghiệp công nghiệp hoá, hiện đại hoá đất nước. Với những ý nghĩa to lớn đó của Robot công nghiệp, chắc chắn ngành công nghiệp chế tạo và ứng dụng Robot sẽ phát triển rất mạnh trong tương lai. Trong lĩnh vực Robot hiện nay, phần Cơ khí (Robot Mechanics), hệ thống Điều khiển (Robot control) và hệ thống Lập trình (Programming system) được coi là các thành phần độc lập và được các nhà sản xuất chào bán độc lập. Vì vậy, với những kiến thức đã học và được sự hướng dẫn của thầy giáo Th.S Lê Huy Tùng, chúng tôi đã nghiên cứu đề tài: “Thiết kế bộ Điều khiển trượt cho tay máy Robot 2 bậc tự do và mô phỏng trên Matlab – Simulink”. SV: Hoàng Đức Thắng - Phan Thanh Hải 6
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do CHƯƠNG I: TỔNG QUAN VỀ ROBOT CÔNG NGHIỆP I.1.Robot công nghiệp: I.1.1. Sự ra đời của Robot công nghiệp : Thuật ngữ “Robot” lần đầu tiên xuất hiện năm 1922 trong tác phẩm “Rosum’s Universal Robot “ của Karal Capek. Theo tiếng Séc thì Robot là người làm tạp dịch. Trong tác phẩm này nhân vật Rosum và con trai ông đã tạo ra những chiếc máy gần giống như con người để hầu hạ con người. Hơn 20 năm sau, ước mơ viễn tưởng của Karel Capek đã bắt đầu hiện thực. Ngay sau chiến tranh thế giới lần thứ 2, ở Mỹ đã xuất hiện những tay máy chép hình điều khiển từ xa, trong các phòng thí nghiệm phóng xạ. Năm 1959, Devol và Engelber đã chế tạo Robot công nghiệp đầu tiên tại công ty Unimation. Năm 1967 Nhật Bản mới nhập chiếc Robot công nghiệp đầu tiên từ công ty AMF của Mỹ. Đến năm 1990 có hơn 40 công ty của Nhật, trong đó có những công ty khổng lồ như Hitachi, Mitsubishi và Honda đã đưa ra thị trường nhiều loại Robot nổi tiếng. Từ những năm 70, việc nghiên cứu nâng cao tính năng của robot đã chú ý nhiều đến sự lắp đặt thêm các cảm biến ngoại tín hiệu để nhận biết môi trường làm việc. Tại trường đại học tổng hợp Stanford, người ta đã tạo ra loại Robot lắp ráp tự động điều khiển bằng vi tính trên cơ sở xử lý thông tin từ các cảm biến lực và thị giác. Vào thời gian này công ty IBM đã chế tạo Robot có các cảm biến xúc giác và cảm biến lực điều khiển bằng máy vi tính để lắp ráp các máy in gồm 20 cụm chi tiết . Những năm 90 do áp dụng rộng rãi các tiến bộ khoa học về vi xử lý và công nghệ thông tin, số lượng Robot công nghiệp đã tăng nhanh, giá thành giảm đi rõ rệt, tính năng đã có nhiều bước tiến vượt bậc. Nhờ vậy Robot công nghiệp đã có vị trí quan trọng trong các dây truyền sản xuất hiện đại. Ngày nay, chuyên ngành khoa học nghiên cứu về Robot “Robotics” đã trở thành một lĩnh vực rộng trong khoa học, bao gồm các SV: Hoàng Đức Thắng - Phan Thanh Hải 7
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do vấn đề cấu trúc cơ cấu động học, động lực học, lập trình quỹ đạo, cảm biến tín hiệu, điều khiển chuyển động v.v I.1.2.Phân loại tay máy Robot công nghiệp: Ngày nay, khi nói đến Robot thường ta hay hình dung ra một cơ chế máy móc tương tự con người, có khả năng sử dụng công cụ lao động để thực hiện các công việc thay cho con người, thậm chí có thể tính toán hay có khả năng hành động theo ý chí. Trong thực tiễn kỹ thuật, khái niệm Robot hiện đại được hiểu khá rộng, mà theo đó Robot là “tất cả các hệ thống kỹ thuật có khả năng cảm nhận và xử lý thông tin cảm nhận được, để sau đó đưa ra hành xử thích hợp”. Theo cách hiểu này, các hệ thống xe tự hành, hay thậm chí một thiết bị xây dựng có trang bị cảm biến thích hợp như Camera, cũng được gọi là Robot. Các khái niệm như Hexapod, Parallel Robot, Tripod, Gait Biped, Manipulator Robocar hay Mobile Robot nhằm chỉ vào các hệ thống Robot không còn gắn liền với các hình dung ban đầu của con người. Trong nội dung đồ án chỉ nhằm vào đối tượng Robot công nghiệp (RBCN), thực chất là một thiết bị tay máy (Handling Equipment). Công nghệ tay máy (Handling Technology) là công nghệ của dạng thiết bị kỹ thuật có khả năng thực hiện các chuyển động theo nhiều trục trong không gian, tương tự như ở con người. Về cơ bản có thể phân thiết bị tay máy (hình 1.1) thành 2 loại chính : Điều khiển (ĐK) theo chương trình hay ĐK thông minh : SV: Hoàng Đức Thắng - Phan Thanh Hải 8
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Handling Equipments Điều khiển Điều khiển theo thông minh chương trình Chương trình Chương trình cứng linh hoạt Manipulators, Máy bốc dỡ, Robot công xếp đặt nghiệp Telemanipulators Hình 1.1 : Phân loại thiết bị tay máy + Loại ĐK theo chương trình gồm 2 họ: Chương trình cứng : Các thiết bị bốc dỡ, xếp đặt có chương trình hoạt động cố định. Ta hay gặp họ này trong các hệ thống kho hiện đại. Chúng có rất ít trục chuyển động và chỉ thu thập thông tin về quãng đường qua các tiếp điểm hành trình. Ta không thể ĐK chúng theo một quỹ đạo mong muốn. Chương trình linh hoạt : Là họ Robot mà người sử dụng có khả năng thay đổi chương trình ĐK chúng tuỳ theo đối tượng công tác. Ta hay gặp chúng trong các công đoạn như hàn, sơn hay lắp ráp của công nghiệp Ôtô. Trong hình 1.1 ta gọi là Robot công nghiệp. + Loại ĐK thông minh có 2 kiểu chính : SV: Hoàng Đức Thắng - Phan Thanh Hải 9
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Manipulator: Là loại tay máy được ĐK trực tiếp bởi con người, có khả năng lặp lại các chuyển động của tay người. Bản chất là dạng thiết bị hỗ trợ cho sự khéo léo, cho trí tuệ, cho hệ thống giác quan (Complex Sensorics) và kinh nghiệm của người sử dụng. Hay được sử dụng trong các nhiệm vụ cần chuyển động phức hợp có tính chính xác cao, hay môi trường nguy hiểm cho sức khoẻ, môi trường khó tiếp cận v.v Telemanipulator: Là loại Manipulator được điều khiển từ xa và người ĐK phải sử dụng hệ thống Camera để quan sát môi trường sử dụng. Theo tiêu chuẩn châu Âu EN775 và VDI 2860 của Đức có thể hiểu “Robot công nghiệp là một Automat sử dụng vạn năng để tạo chuyển động nhiều trục, có khả năng lập trình linh hoạt các chuỗi chuyển động và quãng đường (góc) để tạo nên chuyển động theo quỹ đạo. Chúng có thể được trang bị thêm các ngón (Grippe), dụng cụ hay các công cụ gia công và có thể thực hiện các nhiệm vụ của đôi tay (Handling) hay các nhiệm vụ gia công khác” Như vậy, RBCN khác các loại tay máy còn lại ở 2 điểm chính là “sử dụng vạn năng” và “khả năng lập trình linh hoạt”. I.2. Ứng dụng của Robot công nghiệp : I.2.1.Mục tiêu ứng dụng Robot công nghiệp : Mục tiêu ứng dụng Robot công nghiệp nhằm nâng cao năng suất dây truyền công nghệ, giảm giá thành, nâng cao chất lượng và khả năng cạnh tranh của sản phẩm, đồng thời cải thiện điều kiện lao động. Điều đó xuất phát từ những ưu điểm cơ bản của Robot đó là : - Robot có thể thực hiện một quy trình thao tác hợp lý bằng hoặc hơn người thợ lành nghề một cách ổn định trong suốt thời gian dài làm việc. Do đó Robot giúp nâng cao chất lượng và khả năng cạnh tranh của sản phẩm. - Khả năng giảm giá thành sản phẩm do ứng dụng Robot là vì giảm được đáng kể chi phí cho người lao động. - Robot giúp tăng năng suất dây chuyền công nghệ. SV: Hoàng Đức Thắng - Phan Thanh Hải 10
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do - Robot giúp cải thiện điều kiện lao động. Đó là ưu điểm nổi bật nhất mà chúng ta cần quan tâm. Trong thực tế sản xuất có rất nhiều nơi người lao động phải làm việc trong môi trường ô nhiễm, ẩm ướt, nóng nực. Thậm chí rất độc hại đến sức khoẻ và tính mạng như môi trường hoá chất, điện từ, phóng xạ I.2.2.Các lĩnh vực ứng dụng Robot công nghiệp : Robot công nghiệp được ứng dụng rất rộng rãi trong sản xuất, xin được nêu ra một số lĩnh vực chủ yếu : - Kỹ nghệ đúc - Gia công áp lực - Các quá trình hàn và nhiệt luyện - Công nghệ gia công lắp ráp - Phun sơn, vận chuyển hàng hoá (Robocar) I.2.3. Các xu thế ứng dụng Robot trong tương lai : - Robot ngày càng thay thế nhiều lao động - Robot ngày càng trở lên chuyên dụng - Robot ngày càng đảm nhận được nhiều loại công việc lắp ráp - Robot di động ngày càng trở lên phổ biến - Robot ngày càng trở lên tinh khôn I.2.4. Tình hình tiếp cận và ứng dụng Robot công nghiệp ở Việt Nam : Trong giai đoạn trước năm 1990, hầu như trong nước hoàn toàn chưa du nhập về kỹ thuật Robot, thậm chí chưa nhận được nhiều thông tin kỹ thuật về lĩnh vực này. Tuy vậy, với mục tiêu chủ yếu là tiếp cận lĩnh vực mới mẻ này trong nước đã có triển khai các đề tài nghiên cứu khoa học cấp nhà nước: Đề tài 58.01.03 và 52B.03.01. Giai đoạn tiếp theo từ năm 1990 các ngành công nghiệp trong nước bắt đầu đổi mới. Nhiều cơ sở đã nhập ngoại nhiều loại Robot công nghiệp SV: Hoàng Đức Thắng - Phan Thanh Hải 11
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do phục vụ các công việc như: tháo lắp dụng cụ, lắp ráp linh kiện điện tử, hàn vỏ Ôtô xe máy, phun phủ các bề mặt Một sự kiện đáng chú ý là tháng 4 năm 1998, nhà máy Rorze/Robotech đã bước vào hoạt động ở khu công nghiệp Nomura Hải Phòng. Đây là nhà máy đầu tiên ở Việt Nam chế tạo và lắp ráp Robot. Những năm gần đây, Trung tâm nghiên cứu kỹ thuật Tự động hóa, Trường đại học Bách Khoa Hà Nội, đã nghiên cứu thiết kế một kiểu Robot mới là Robot RP. Robot RP thuộc loại Robot phỏng sinh (bắt chước cơ cấu tay người). Hiện nay đã chế tạo 2 mẫu: Robot RPS-406 dùng để phun men và Robot RPS-4102 dùng trong công nghệ bề mặt. Ngoài ra Trung tâm còn chế tạo các loại Robot khác như: Robot SCA mini dùng để dạy học, Robocar công nghiệp phục vụ phân xưởng, Robocar chữ thập đỏ cho người tàn tật Bên cạnh đó còn xây dựng các thuật toán mới để điều khiển Robot, xây dựng “thư viện” các mô hình của Robot trên máy tính I.3.Cấu trúc của Robot công nghiệp: I.3.1.Các bộ phận cấu thành Robot công nghiệp : Trên hình 1.2 giới thiệu các bộ phận chủ yếu của Robot công nghiệp: Tay máy gồm các bộ phận: Đế 1 đặt cố định hoặc gắn liền với xe di động 2, thân 3, cánh tay trên 4, cánh tay dưới 5, bàn kẹp 6. SV: Hoàng Đức Thắng - Phan Thanh Hải 12
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Hình 1.2: Các bộ phận cấu thành Robot công nghiệp Hệ thống truyền dẫn động có thể là cơ khí, thuỷ khí hoặc điện khí: là bộ phận chủ yếu tạo nên sự chuyển dịch các khớp động. Hệ thống điều khiển đảm bảo sự hoạt động của Robot theo các thông tin đặt trước hoặc nhận biết trong quá trình làm việc. Hệ thống cảm biến tín hiệu thực hiện việc nhận biết và biến đổi thông tin về hoạt động của bản thân Robot (cảm biến nội tín hiệu) và của môi trường, đối tượng mà Robot phục vụ (cảm biến ngoại tín hiệu). I.3.2.Bậc tự do và các toạ độ suy rộng : I.3.2.1.Bậc tự do : Robot công nghiệp là loại thiết bị tự động nhiều công dụng. Cơ cấu tay máy của chúng phải được cấu tạo sao cho bàn kẹp giữ vật kẹp theo một hướng nhất định nào đó và di chuyển dễ dàng trong vùng làm việc. Muốn vậy cơ cấu tay máy phải đạt được một số bậc tự do chuyển động. Thông thường các khâu của cơ cấu tay máy được nối ghép với nhau bằng các khớp quay hoặc khớp tịnh tiến. Gọi chung chúng là khớp động. Các khớp quay hoặc khớp tịnh tiến đều thuộc khớp động học loại 5. SV: Hoàng Đức Thắng - Phan Thanh Hải 13
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Công thức tính số bậc tự do : 5 W= 6n - i (1.1)  pi 1 với n : số khâu động Pi : số khớp loại i Ví dụ: Tay máy có 2 khớp quay như hình vẽ 1.3 : Số khâu động n = 2 Khớp quay là khớp loại 5 . Do đó W = 6.2 – ( 5.1 + 5.1) = 2 bậc tự do Hình 1.3: Tay máy 2 khớp quay I.3.2.2. Toạ độ suy rộng : Các cấu hình khác nhau của cơ cấu tay máy trong từng thời điểm xác định bằng các độ dịch chuyển góc hoặc độ dịch chuyển dài của các khớp quay hoặc khớp tịnh tiến. Các độ dịch chuyển tức thời đó, so với giá trị ban đầu nào đó lấy làm mốc tính toán, được gọi là các toạ độ suy rộng (generalized joint coordinates). Ở đây ta gọi chúng là các biến khớp (toạ độ suy rộng) của cơ cấu tay máy và biểu thị bằng : (1 ) qi  i i  i Si (1.2) 1,®èi víi khíp quay với  i 0,®èi víi khíp tÞnh tiÕn - Độ dịch chuyển góc của các khớp quay  i SV: Hoàng Đức Thắng - Phan Thanh Hải 14
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do - Độ dịch chuyển tịnh tiến của các khớp tịnh tiến S i I.3.3.Nhiệm vụ lập trình điều khiển Robot: I.3.3.1. Định vị và định hướng tại “điểm tác động cuối” : Khâu cuối cùng của tay máy thường là bàn kẹp (gripper) hoặc là khâu gắn liền với dụng cụ thao tác (tool). Điểm mút của khâu cuối cùng là điểm đáng quan tâm nhất vì đó là điểm tác động của Robot lên đối tác và được gọi là “điểm tác động cuối” (end-effector). Trên hình 1.4 điểm E là “điểm tác động cuối”. Hình 1.4: Định vị và định hướng tại “ điểm tác động cuối” Chính tại “điểm tác động cuối” E này cần quan tâm không những vị trí nó chiếm trong không gian làm việc mà cả hướng tác động của khâu cuối đó. Vị trí của điểm E được xác định bằng 3 toạ độ x E, yE, zE trong hệ trục toạ độ cố định. Còn hướng tác động của khâu cuối có thể xác định bằng 3 trục xn,yn, zn gắn liền với khâu cuối tại điểm E, hoặc bằng 3 thông số góc , , nào đó. I.3.3.2. Lập trình điều khiển Robot công nghiệp : Trên hình 1.5 mô tả 1 sơ đồ lập trình điều khiển Robot công nghiệp. Khi robot nhận nhiệm vụ thực hiện một quy trình công nghệ nào đó, ví dụ SV: Hoàng Đức Thắng - Phan Thanh Hải 15
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do “điểm tác động cuối” E phải bám theo một hành trình cho trước. Quỹ đạo hành trình này thường cho biết trong hệ toạ độ Đề các x 0, y0, z0 cố định. Ở mỗi vị trí mà điểm E đi qua xác định bằng 3 toạ độ cố định x E, yE, zE và 3 thông số góc định hướng , , . Từ các thông số trong hệ toạ độ Đề các đó tính toán các giá trị biến khớp q i tương ứng với mỗi thời điểm t. Đó là nội dung của bài toán Động học ngược sẽ trình bày trong chương II. Quỹ đạo trong hệ toạ độ Đề các Máy tính (xE,yE,zE, ,,) Quỹ đạo trong hệ toạ độ Đề các q1 (xE,yE,zE, ,,) q Chương trình điều khiển 2 Hệ trợ động chấp hành Các giản đồ Biến đổi qi(t) ROBOT Hệ trợ động chấp hành Hình 1.5: Sơ đồ lập trình điều khiển I.4. Các phép biến đổi toán học cho Robot : I.4.1.Biến đổi toạ độ dùng Ma trận: I.4.1.1. Vector điểm và toạ độ thuần nhất : Vector điểm (point vector) dùng để mô tả vị trí của điểm trong không gian 3 chiều. SV: Hoàng Đức Thắng - Phan Thanh Hải 16
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Trong không gian 3 chiều, một điểm M có thể được biểu diễn bằng nhiều vector trong các hệ toạ độ (coordinate frame) khác nhau: Trong hệ toạ độ oixiyizi điểm M xác định bằng vector ri : T r i (r xi,r yi,r zi) (1.3) và cùng điểm M đó trong hệ toạ độ ojxjyjzj được mô tả bởi vector rj : T r j (r xj,r yj,r zj) (1.4) Ký hiệu ( )T là biểu thị phép chuyển vị (Transportation) vector hàng thành vector cột. z z i M j ri rj yi Oi yj xj xi Hình 1.6: Biểu diễn 1 điểm trong không gian T Vector r (r x,r y,r z) trong không gian 3 chiều, nếu được bổ sung thêm một thành phần thứ 4 và thể hiện bằng 1 vector mở rộng :  r (r x,r yr z,) (1.5) thì đó là cách biểu diễn vector điểm trong không gian toạ độ thuần nhất (homogeneous coordinate). Để đơn giản có thể bỏ qua ký hiệu ( ˜ ) đối với vector mở rộng (1.5) Các toạ độ thực của vector mở rộng này vẫn là:    r x r y r z (1.6) r x  r y  r z  SV: Hoàng Đức Thắng - Phan Thanh Hải 17
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Không phải duy nhất có một cách biểu diễn vector trong không gian tọa độ thuần nhất, mà nó phụ thuộc vào giá trị của  . Nếu lấy  = 1 thì các tọa độ biểu diễn bằng toạ độ có thực. Trong trường hợp này vector mở rộng được viết là: T r (r x,r y,r z) (1.7) Nếu lấy  ≠ 1 thì các toạ độ biểu diễn gấp  lần toạ độ thực, nên có thể gọi  là hệ số tỷ lệ. Khi cần biểu diễn sự thay đổi toạ độ kèm theo thì có sự biến dạng tỷ lệ thì dùng  ≠ 1. I.4.1.2.Quay hệ toạ độ dùng Ma trận 3x3: Trước hết thiết lập quan hệ giữa 2 hệ toạ độ XYZ và UVW chuyển động quay tương đối với nhau khi gốc O của 2 hệ vẫn trùng nhau (hình 1.7) Z W M V Y U X Hình 1.7: Các hệ toạ độ Gọi (ix, jy, kz) và (iu, jv, kw) là các vector đơn vị chỉ phương các trục OXYZ và OUVW tương ứng. Một điểm M nào đó được biểu diễn trong hệ toạ độ OXYZ bằng vector: T rxyz=( rx,ry,rz) (1.8) còn trong hệ toạ độ OUVW bằng vector: T ruvw = ( ru,rv,rw) (1.9) Như vậy : SV: Hoàng Đức Thắng - Phan Thanh Hải 18
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do r = ruvw= ruiu + rvjv + rwkw r = rxyz= rxix + ryjy + rzkz (1.10) Từ đó ta có .r  r x i x i xi ur u i x j vr v i xk wr w .r (1.11) r y j y j yi ur u j y j vr v j yk wr w .r r z k z k zi ur u k z j vr v k zk wr w  Hay viết dưới dạng ma trận: j r x i xi u i x v i xk w r u . (1.12) r y j yi u j y j v j yk w r v r z r w k zi u k z j v k zk w Gọi R là Ma trận quay (rotation) 3x3 với các phần tử là tích vô hướng 2 vector chỉ phương các trục tương ứng của 2 hệ toạ độ OXYZ và OUVW. Vậy (1.12) được viết lại là: r xyz R.r uvw   (1.13) R 1. r uvw r xyz  I.4.1.3.Biến đổi Ma trận dùng toạ độ thuần nhất: Bây giờ thiết lập quan hệ giữa 2 hệ toạ độ: hệ toạ độ o jxjyjzj sang hệ toạ độ mới o ixiyizi. Chúng không những quay tương đối với nhau mà tịnh tiến cả gốc toạ độ: gốc oj xác định trong hệ xiyizi bằng vector p: p=(a,-b,-c,1)T (1.14) Giả sử vị trí của điểm M trong hệ toạ độ x jyjzj được xác định bằng vector rj: T rj = (xjyjzj,1) (1.15) và trong hệ toạ độ xiyizi điểm M được xác định bằng vector ri: T ri = (xiyizi,1) (1.16) SV: Hoàng Đức Thắng - Phan Thanh Hải 19
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Từ hình (1.8) có thể dễ dàng thiết lập mối quan hệ giữa các toạ độ: xi x j at j  cos sin b y i y j z j t j  (1.17) sin cos c z i y j z j t j 1 t i t j  xj oj zj yj xi a c oi yi b zi Hình 1.8: Các hệ toạ độ Sắp xếp các hệ số ứng với xj,yj,zj và tj thành một ma trận: 1 0 0 a 0 cos sin b (1.18) T ij 0 sin cos c 0 0 0 1 và viết phương trình biến đổi toạ độ như sau: ri = Tij rj (1.19) Ma trận Tij biểu thị bằng ma trận 4x4 như phương trình (1.18) và gọi là ma trận thuần nhất. Nó dùng để biến đổi vector mở rộng từ hệ toạ độ thuần nhất này sang hệ toạ độ thuần nhất kia. I.4.1.4. Ý nghĩa hình học của Ma trận thuần nhất: SV: Hoàng Đức Thắng - Phan Thanh Hải 20
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Từ (3.19) nhận thấy ma trận thuần nhất 4x4 là một ma trận gồm 4 khối : 1 0 0 a 0 cos -sin -b T (1.20) ij 0 sin cos c 0 0 0 1 Hoặc viết rút gọn là: Rij p Tij (1.21) 0 1 Trong đó: Rij - ma trận quay 3x3 p – ma trận 3x1 biểu thị 3 toạ độ của điểm gốc hệ toạ độ 0 j trong hệ toạ độ oi, xi, yi, zi 1x3 – ma trận không 1x1 – ma trận đơn vị Như vậy ma trận thuần nhất 4x4 là ma trận 3x3 mở rộng, thêm ma trận 3x1 biểu thị sự chuyển dịch gốc toạ độ và phần tử a 44 biểu thị hệ số tỷ lệ. Dễ dàng nhận thấy ma trận Rij chính là ma trận quay 3x3, nếu suy từ ma trận quay trong (1.12) sang trường hợp hình 1.8 ta có: cos(xi ,x ) cos(xi ,y ) cos(xi ,z ) j j j cos(y , ) cos(y , ) cos(y , ) (1.22) Rij aij i x j i yj i z j cos(z , ) cos(z , ) cos(z , ) i x j i yj i z j và các góc cosin chỉ phương này đều liên hệ đến góc (hình 1.8). Nếu chú ý về quan hệ giữa 2 cặp trục,ví dụ, cos(xi,yj) = cos(yi, xj) ở đây dễ dàng nhận được biểu thức: SV: Hoàng Đức Thắng - Phan Thanh Hải 21
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do -1 T Rij Rij Rij (1.23) Mô tả tổng quát hơn nếu một điểm M nào đó được xác định trong hệ toạ độ thuần nhất UVW bằng vectơ mở rộng r uvw , thì trong hệ toạ độ thuần nhất XYZ điểm đó xác định bằng vector mở rộng rxyz: Rxyz = T.ruvw (1.24) Trong đó T là ma trận thuần nhất 4x4, có thể viết khai triển ở dạng sau: nx sx ax px n s a p T y y y y (1.25) n s a p z z z z 0 0 0 1 R p hoặc T (1.26) 0 1 Ta tìm hiểu ý nghĩa hình học của ma trận T. Như đã trình bày khi phân tích các khối của ma trận 4x4, ma trận 3x1 tương ứng với toạ độ điểm gốc của hệ toạ độ UVW biểu diễn trong hệ XYZ. Nếu 2 gốc toạ độ trùng nhau thì các thành phần của ma trận 3x1 này đều là 0. Khi đó xét trường hợp: T ruvw (1,0,0,1) tức là rxyz = iu thì dễ dàng nhận thấy cột thứ nhất hoặc vectơ n của ma trận (1.25) chính là các toạ độ của vectơ chỉ phương trục OU biểu diễn trong hệ toạ độ XYZ. Tương tự khi xét các trường hợp T ruvw (0,1,0,1) T và ruvw (0,0,1,1) SV: Hoàng Đức Thắng - Phan Thanh Hải 22
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do cũng đi đến nhận xét cột thứ 2 (hoặc vectơ s) ứng với các toạ độ của vectơ chỉ phương trục OV và cột thứ 3 (hoặc vectơ a) ứng với các toạ độ vector chỉ phương trục OW. Như vậy, ma trận thuần nhất T 4x4 hoàn toàn xác định vị trí và định hướng của hệ toạ độ UVW so với hệ toạ độ XYZ. Đó là ý nghĩa hình học của ma trận thuần nhất 4x4. I.4.2.Các phép biến đổi cơ bản: I.4.2.1.Phép biến đổi tịnh tiến: Từ (1.18) hoặc (1.25), biểu thị ma trận thuần nhất khi chỉ có biến đổi tịnh tiến mà không có quay ( 0 ), ta có: 1 0 0 p x 0 1 0 p T y =( , , ) (1.27) T p p x p y p z 0 0 1 p z 0 0 0 1 Đó là ma trận biến đổi tịnh tiến (Tranlation) Gọi u là vector biểu diễn một điểm trong không gian cần dịch chuyển tịnh tiến: u (x, y, z)T và p là vector chỉ hướng và độ dài cần dịch chuyển p ( , , )T p x p y p z thì v là vector biểu diễn điểm toạ độ trong không gian đã được tịnh tiến tới: v ( v ( , , )T u (1.28) T p T p p x p y p z I.4.2.2. Phép quay quanh các trục toạ độ : SV: Hoàng Đức Thắng - Phan Thanh Hải 23
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Từ ma trận quay 3x3 trong biểu thức (1.12) ta xây dựng ma trận R(x, ) cho trường hợp hệ toạ độ UVW quay quanh trục OX một góc nào đó. Trong trường hợp này i x i u : 1 0 0 0 0 cos sin 0 R(x, ) (1.29) 0 sin cos 0 0 0 0 1 Tương ứng cho trường hợp quay quanh trục OY một góc : cos 0 sin 0 0 1 0 0 R(y, ) (1.30) sin 0 cos 0 0 0 0 1 và trường hợp quay quanh trục OZ một góc  : cos sin 0 0 sin cos 0 0 R(z, ) (1.31) 0 0 1 0 0 0 0 1 Cột thứ 4 của các ma trận 4x4 trên có 3 phần tử đều bằng 0 vì ở đây không có sự tịnh tiến. Các ma trận này được gọi là các ma trận quay (rotation) cơ bản. Các ma trận quay khác có thể xây dựng từ các ma trận cơ bản này. SV: Hoàng Đức Thắng - Phan Thanh Hải 24
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do CHƯƠNG II: HỆ PHƯƠNG TRÌNH ĐỘNG HỌC VÀ ĐỘNG LỰC HỌC CỦA ROBOT CÔNG NGHIỆP: II.1. Hệ phương trình động học Robot : II.1.1. Đặt vấn đề : Cơ cấu chấp hành của Robot thường là một cơ cấu hở gồm một chuỗi các khâu (link) nối với nhau bằng các khớp (joints). Các khớp động này là khớp quay (R) hoặc khớp tịnh tiến (T). Để Robot có thể thao tác linh hoạt cơ cấu chấp hành của nó phải có cấu tạo sao cho điểm mút của khâu cuối cùng đảm bảo dễ dàng di chuyển theo một quỹ đạo nào đó, đồng thời khâu này có một hướng nhất định theo yêu cầu. Khâu cuối cùng này thường là bàn kẹp (griper), điểm mút của nó chính là “điểm tác động cuối” E (end-effector). Để xét vị trí và hướng của E trong không gian ta gắn vào nó một hệ toạ độ động thứ n và gắn với mỗi khâu động một hệ toạ độ khác, còn gắn liền với giá đỡ một hệ toạ độ cố định. Đánh số ký hiệu các hệ này từ 0 đến n bắt đầu từ giá cố định. Khi khảo sát chuyển động của Robot cần biết “định vị và định hướng” tại điểm tác động cuối trong mọi thời điểm. Các lời giải của bài toán này được xác định từ những phương trình Động học của Robot. Các phương trình này là mô hình Động học của Robot. Chúng được xây dựng trên cơ sở thiết lập các mối quan hệ giữa các hệ toạ độ động nói trên so với hệ toạ độ cố định. II.1.2. Xác định trạng thái của Robot tai điểm tác động cuối : Trạng thái của Robot tại “điểm tác động cuối” hoàn toàn xác định bằng sự định vị và định hướng tại điểm tác động cuối đó. Như đã đề cập ở phần I.4.1.4 biểu thị sự định vị và định hướng đó bằng ma trận trạng thái cuối TE : SV: Hoàng Đức Thắng - Phan Thanh Hải 25
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do nx sx ax px n n a p y y y y TE (2.1) nz sz az pz 0 0 0 1 Trong đó các phần tử của ma trận 3x1 là toạ độ p x , py, pz của “điểm tác động cuối” E. Mỗi cột của ma trận quay 3x3 là một vectơ đơn vị chỉ phương một trục của hệ toạ độ động NSA (chính là UVW) biểu diễn trong toạ độ cố định XYZ. Hệ toạ độ gắn liền với bàn kẹp của Robot có các vectơ đơn vị chỉ phương các trục như sau : a - vector có hướng tiếp cận (approach) với đối tác . s - vector có hướng đường trượt (sliding) đóng mở bàn kẹp . n - vector pháp tuyến (normal). II.1.3. Mô hình động học : II.1.3.1. Ma trận quan hệ : Chọn hệ toạ độ cố định gắn liền với giá đỡ và các hệ toạ độ gắn với từng khâu động. Ký hiệu các hệ toạ độ này từ 0 đến n, kể từ giá cố định trở đi. Một điểm bất kì nào đó trong không gian được xác định trong hệ toạ độ thứ i bằng bán kính ri và trong hệ toạ độ cố định x0, y0, z0 được xác định bằng bán kính vector r0 : r0 = A1A2 Airi (2.2) hoặc r0 = Tiri (2.3) với Ti = A1A2 Ai , i= 1, 2, n (2.4) Trong đó ma trận A 1 mô tả vị trí hướng của khâu đầu tiên; ma trận A2 mô tả vị trí và hướng của khâu thứ 2 so với khâu đầu; ma trận A i mô tả vị trí và hướng của khâu thứ i so với khâu thứ i-1. SV: Hoàng Đức Thắng - Phan Thanh Hải 26
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Như vậy, tích của các ma trận A i là ma trận Ti mô tả vị trí và hướng của khâu thứ i so với giá trị cố định. Thường kí hiệu ma trận T với 2 chỉ số: trên và dưới. Chỉ số dưới chỉ khâu đang xét còn chỉ số trên để chỉ toạ độ được dùng để đối chiếu. Ví dụ, biểu thức (2.4) có thể viết lại là : 0 1 Ti Ti A1 Ti (2.5) 1 với Ti A2 A3 Ai (2.6) là ma trận mô tả vị trí và hướng của khâu thứ i so với khâu thứ nhất. Trong kí hiệu thường bỏ qua chỉ số trên nếu chỉ số đó bằng 0. Denavit & Hartenberg đã đề xuất dùng ma trận thuần nhất 4x4 mô tả quan hệ giữa 2 khâu liên tiếp trong cơ cấu không gian . II.1.3.2. Bộ thông số DH : Dưới đây trình bày cách xây dựng các hệ toạ động đối với 2 khâu động liên tiếp i và i+1. Hình dưới đây là trường hợp 2 khớp động liên tiếp là 2 khớp quay. Hình 2.1: Các hệ toạ độ đối với 2 khâu động liên tiếp SV: Hoàng Đức Thắng - Phan Thanh Hải 27
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Trước hết xác định bộ thông số cơ bản giữa 2 trục quay của khớp động i+1 và i : ai là độ dài đường vuông góc chung giữa 2 trục khớp động i+1 và i . i là góc chéo giữa 2 trục khớp động i+1 và i . di là khoảng cách đo dọc trục khớp động i từ đường vuông góc chung giữa trục khớp động i+1 và trục khớp động i tới đường vuông góc chung giữa khớp động i và trục khớp động i -1. i là góc giữa 2 đường vuông góc chung nói trên. Bộ thông số này được gọi là bộ thông số Denavit – Hartenberg (DH). Biến khớp (joint variable): Nếu khớp động i là khớp quay thì biến khớp là i Nếu khớp động i là khớp tịnh tiến thì biến khớp là di Để kí hiệu thêm biến khớp dùng thêm dấu * và trong trường hợp khớp tịnh tiến thì ai được xem là bằng 0. II.1.3.3. Thiết lập hệ toạ độ : Gốc của hệ toạ độ gắn liền với khâu thứ i (gọi là hệ toạ độ thứ i) đặt tại giao điểm giữa đường vuông góc chung (ai) và trục khớp động i+1. Trường hợp 2 trục giao nhau thì gốc hệ toạ độ lấy trùng với giao điểm đó. Nếu 2 trục song song với nhau thì chọn gốc toạ độ là điểm bất kì trên trục khớp động i+1. Trục zi của hệ toạ độ thứ i nằm dọc theo trục khớp động i+1. Trục xi của hệ toạ độ thứ i nằm dọc theo đường vuông góc chung hướng từ khớp động i đến khớp động i+1. Trường hợp 2 trục giao nhau, hướng trục xi trùng với hướng vector tích zi x zi-1, tức là vuông góc với mặt phẳng chứa zi, zi-1. Ví dụ : Xét tay máy có 2 khâu phẳng như hình 2.2. SV: Hoàng Đức Thắng - Phan Thanh Hải 28
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Hình 2.2: Tay máy 2 khâu phẳng (vị trí bất kỳ) Gắn các hệ toạ độ với các khâu như hình vẽ : - Trục z0 , z1 và z2 vuông góc với mặt tờ giấy. - Hệ toạ độ cố định là o0x0y0z0 chiều x0 hướng từ o0 đến o1. - Hệ toạ độ o1x1y1z1 có gốc o1 đặt tại tâm trục khớp động 2. - Hệ toạ độ o 2x2y2z2 có gốc o 2 đặt tại tâm trục khớp động cuối khâu 2. Bảng thông số DH của tay máy này như sau : Khâu i i ai di * 1 0 a1 0  1 * 2  2 0 a2 0 II.1.3.4. Mô hình biến đổi : SV: Hoàng Đức Thắng - Phan Thanh Hải 29
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Trên cơ sở đã xây dựng các hệ toạ độ với 2 khâu động liên tiếp như trên đã trình bày. Có thể thiết lập mối quan hệ giữa 2 hệ toạ độ liên tiếp theo 4 phép biến đổi : + Quay quanh trục z1-1 góc i . + Tịnh tiến dọc trục zi-1 một đoạn di . + Tịnh tiến dọc trục xi-1 (đã trùng với xi) một đoạn ai . + Quay quanh trục xi một góc i . Bốn phép biến đổi này được biểu thị bằng tích các ma trận thuần nhất sau Ai = R(z,i).Tp(0,0,di).Tp(ai,0,0).R(x, i) (2.7) Các ma trận ở vế phải phương trình (2.7) tính theo các công thức (1.27),(1.29),(1.31). Sau khi thực hiện phép nhân các ma trận nói trên, ta có: Ci S iC i S i S i aiCi A S i CiC i Ci S i ai S i (2.8) i 0 S i C i d i 0 0 0 1 Trong khớp tịnh tiến : a = 0 . II.1.3.5. Phương trình động học : Với Robot có n khâu, ma trận mô tả vị trí và hướng điểm cuối E của tay máy được miêu tả : Tn = A1A2 An (2.9) Mặt khác, hệ toạ độ tại “điểm tác động cuối” này được mô tả bằng ma trận TE. Vì vậy hiển nhiên là: TE = Tn (2.10) Tức là ta có : SV: Hoàng Đức Thắng - Phan Thanh Hải 30
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do nx sx ax px n n a p y y y y Tn (2.11) nz sz az pz 0 0 0 1 Phương trình (2.11) là phương trình động học cơ bản của Robot. II.2. Tổng hợp chuyển động Robot : II.2.1. Nhiệm vụ : Nhiệm vụ tổng hợp chuyển động bao gồm việc xác định các bộ lời giải qi(t), (i = 1, , n), với qi là toạ độ suy rộng hoặc là biến khớp. Biết quy luật chuyển động của bàn kẹp, cần xác định quy luật thay đổi các biến khớp tương ứng. Đó là nội dung chính của việc tổng hợp quỹ đạo chuyển động Robot. Có thể xem quỹ đạo chuyển động là tập hợp liên tiếp các vị trí khác nhau của bàn kẹp. Tại mỗi vị trí trên quỹ đạo cần xác định bộ thông số các biến khớp qi. Đó là nội dung của bài toán động học ngược (inverse kinematics problem) của Robot. II.2.2. Bài toán động học ngược : Bài toán động học ngược được đặc biệt quan tâm vì lời giải của nó là cơ sở chủ yếu để xây dựng chương trình điều khiển chuyển động của Robot bám theo quỹ đạo cho trước. Xuất phát từ phương trình động học cơ bản (2.11) ta có : nx sx ax px ny ny ay py Tn = A1A2 An = (2.12) nz sz az pz 0 0 0 1 Các ma trận Ai là hàm của các biến khớp q i. Vector định vị bàn kẹp T p = (px,py,pz) cũng là hàm của q i. Các vector n, s, a là các vector đơn vị chỉ phương các trục của hệ toạ độ gắn liền với bàn kẹp biểu diễn trong hệ toạ độ cố định XYZ. Các vector này vuông góc với nhau từng đôi một nên SV: Hoàng Đức Thắng - Phan Thanh Hải 31
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do trong 9 thành phần của chúng chỉ tồn tại độc lập chỉ có 3 thành phần. Hai ma trận ở vế phải và vế trái của phương trình (2.12) đều là các ma trận thuần nhất 4x4. So sánh các phần tử tương ứng của 2 ma trận trên ta có 6 phương trình độc lập với các ẩn qi (i = 1, 2, ,n). II.2.3. Các phương pháp giải bài toán động học ngược : Trường hợp tổng quát ta xét hệ phương trình động học của Robot có n bậc tự do. Vế trái của phương trình (2.12) theo các kí hiệu như (2.4)-(2.6) có thể viết lại như sau: i T n T i. T n (2.13) 1 Nhân 2 vế của (2.13) với T i ta có: A -1 A -1A -1 T = i (2.14) i 2 1 n Tn Kết hợp (2.12) ta có: nx sx ax px -1 -1 -1 n y s y a y p y i =Ai A2 A1 (2.15) T n nz sz az pz 0 0 0 1 với i=1, ,n-1 Ứng với mỗi giá trị của i, khi so sánh các phần tử tương ứng của 2 ma trận ở biểu thức (2.15) ta có 6 phương trình tồn tại độc lập để xác định biến khớp qi. II.3. Động lực học Robot: II.3.1.Nhiệm vụ và phương pháp phân tích Động lực học Robot: Nghiên cứu Động lực học Robot là giai đoạn cần thiết trong việc phân tích cũng như tổng hợp quá trình điều khiển chuyển động. Trong nghiên cứu Động lực học Robot thường giải quyết 2 nhiệm vụ sau đây : SV: Hoàng Đức Thắng - Phan Thanh Hải 32
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do + Nhiệm vụ thứ nhất là xác định momen và lực động xuất hiện trong quá trình chuyển động. Khi đó quy luật biến đổi của biến khớp q i(t) xem như đã biết. + Nhiệm vụ thứ hai là xác định các sai số động. Lúc này phải khảo sát các phương trình chuyển động của cơ cấu tay máy đồng thời xem xét các đặc tính động lực của động cơ truyền động. Có nhiều phương pháp nghiên cứu Động lực học Robot nhưng thường dùng hơn cả là phương pháp Lagrange bậc 2 vì khi kết hợp với mô hình Động lực học kiểu DH (Denavit-Hartenberg) ta sẽ được các phương trình Động lực học ở dạng vector ma trận, rất gọn nhẹ và thuận tiện cho việc nghiên cứu giải tích và tính toán trên máy tính. Các phương trình Động lực học Robot được thiết lập dựa trên cơ sở phương trình Lagrange bậc 2: d L L F Mi ,i=1, ,n (2.16) dt   q i qi Trong đó : L- hàm Lagrange L = K - P (2.17) K, P- động năng và thế năng của cơ hệ. FMi - động lực, hình thành trong khớp động thứ i khi thực hiện chuyển động. qi - biến khớp (toạ độ suy rộng) - đạo hàm bậc nhất của biến khớp theo thời gian. q i Đồng thời khi mô tả vị trí giữa 2 hệ toạ độ thứ i và i-1 dùng ma trận thuần nhất Ai hoặc viết đầy đủ hơn là i 1 . Dùng ma trận này có thể Ai mô tả vị trí trạng thái trong hệ toạ độ thứ i-1 của một điểm bất kì thuộc hệ toạ độ thứ i. SV: Hoàng Đức Thắng - Phan Thanh Hải 33
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Các biến khớp q i là bộ các thông số dịch chuyển của các khớp động của Robot. Vị trí trạng thái của điểm tác động cuối của Robot hoàn toàn được xác định bởi bộ biến khớp qi này. II.3.2.Vận tốc và gia tốc: Để xây dựng mô hình Động lực học Robot dùng phương trình Lagrange bậc 2, cần biết vận tốc của điểm bất kì trên tay máy. i Điểm M nào đó trong hệ toạ độ i, xác định bằng véc tơ mở rộng ri : i T , ri = ( xi , yi , zi , 1 ) (2.18) i Kí hiệu ri có nghĩa là điểm M cho biết trong hệ toạ độ i và được 0 biểu thị cũng trong hệ toạ độ i. Còn khi dùng kí hiệu ri thì có nghĩa là điểm M cho biết trong hệ toạ độ i, nhưng được biểu thị trong hệ toạ độ x 0, y0, z0, tức là trong hệ toạ độ cơ bản. i 1 Như trước đây, dùng ma trận Ai để mô tả vị trí tương đối giữa hệ 0 toạ độ thứ i đối với hệ toạ độ i-1 và ma trận A iđể mô tả quan hệ giữa hệ toạ độ thứ i và hệ toạ độ cơ bản. 0 i 1 Vậy quan hệ giữa ri và ri có thể biểu thị như sau : 0 0 i ri = Ai ri (2.19) 0 0 1 i 1 với Ai = A1 A2 Ai (2.20) i 1 Ma trận Ai đã có từ biểu thức (2.8): C S S S a Ci i i i i i Ci S C S a S i 1A = i i Ci i Ci i i (2.21) i 0 S C d i i i 0 0 0 1 Biểu thức (2.21) là viết cho trường hợp khớp quay i, còn nếu khớp động là khớp tịnh tiến thì ai =0 và từ (2.21) ta có : SV: Hoàng Đức Thắng - Phan Thanh Hải 34
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do C S S S 0 Ci i i i i S C S 0 i 1A = i i Ci i Ci (2.22) i 0 S C d i i i 0 0 0 1 Đối với khớp quay thì i là biến khớp và đối với khớp tịnh tiến thì di là biến khớp. 0 Các phần tử khác không của ma trận Ai đều là hàm của  j, dj, j và aj (j = 1 ,2 , , i). Trong đó i, j lại là thông số xác định bằng cấu trúc cụ thể của tay máy. Do vậy các phần tử này là hàm của biến khớp q i nói chung (qi  j đối với khớp quay và qi  di đối với khớp tịnh tiến). i Vi phân biểu thức (2.19) với lưu ý rằng các vectơ r ilà không đổi với hệ toạ độ thứ i vì giả thiết rằng các khâu của tay máy là vật rắn tuyệt đối, ta có: 0 d 0 d 0 i V  Vi = ( r ) ( A r ) = i dt i dt i i i 1 0  1 i 1 i 0 1  i 1 i 0  i 0 i = Ai A2 Ai ri Ai A2 Ai ri Ai Ai ri Ai ri i 0 0  Ai i Vi  ri (2.23) j 1 q j i 1 Đạo hàm của ma trận Ai đối với biến khớp q i có thể dễ dàng xác định theo công thức sau : i-1 d Ai i 1 Di Ai (2.24) dqi Trong đó đối với khớp quay : 0 1 0 0 1 0 0 0 Di = (2.25) 0 0 0 0 0 0 0 0 SV: Hoàng Đức Thắng - Phan Thanh Hải 35
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do và đối với khớp tịnh tiến : 0 0 0 0 0 0 0 0 Di = (2.26) 0 0 0 1 0 0 0 0 Trong trường hợp i = 1, 2, ,n ta có : 0  Ai  (A1 A2 Aj 1 Aj Ai 1 Ai ) (2.27) q j q j Trong các ma trận bên vế phải chỉ có Aj phụ thuộc vào qj, do đó theo (2.24) ta có : 0  Ai dAj A1 A2 Aj 1 Ai 1 Ai (2.28) q j dq j với Dj tính theo (2.25) hoặc (2.26) : 0  Ai A1 A2 Aj 1D j Ai 1 Ai (2.29a) q j Phương trình (2.29a) mô tả sự thay đổi vị trí các điểm của khâu thứ i gây nên bởi sự dịch chuyển của khớp động thứ j . Kí hiệu vế trái của (2.29a) là U ij và đơn giản hoá cách viết (2.29a) như sau : 0 j 1 Aj 1Dj Ai , j i Uij (2.29b) 0, j i vậy công thức (2.23) có thể viết lại là: i i V = q r (2.30) i U ij j i j 1 Tiếp theo, từ biểu thức (2.23) xác định gia tốc: dV i 0 A i i 2 A a i i i q q i r (2.31)  qs  s k i dt s 1 q s 1 k 1 q q s s k SV: Hoàng Đức Thắng - Phan Thanh Hải 36
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do II.3.3. Động năng tay máy: Kí hiệu K i là động năng của khâu i ( i =1,2, , n) và dK i là động năng của một chất điểm khối lượng dm thuộc khâu i: 1 1 dK (x 2 y 2 z2 )dm Tr(VV T )dm (2.32) i 2 1 1 1 2 i i Trong đó Tr là vết của ma trận : n TrA =  aii i 1 Ta có : T i i 1 i 1 dK = Tr U q r U q r dm i  ip p i  ir r i 2 p 1 r 1 i i 1 i 1 T T = Tr   Uip ri ri U ir q p q r dm 2 p 1 r 1 i i 1 i i T T = Tr   Uip ( ridm r i )U ir q p q r (2.33) 2 p 1 r 1 Như đã nói ở trên ma trận U ij biểu thị sự thay đổi vị trí của các điểm thuộc khâu i gây nên bởi sự dịch chuyến của khớp động j. Ma trận này đều như nhau tại mọi thời điểm thuộc khâu i và không phụ thuộc vào sự phân bố khối lượng trên khâu i, tức là không phụ thuộc vào dm. Cũng vậy, đạo hàm của biến khớp qi theo thời gian không phụ thuộc vào dm. Do vậy ta có: i i 1 i i T T Ki dKi Tr   Uip ri ri dm Uir q p q r (2.34) 2 p 1 r 1 Phần trong ngoặc đơn là ma trận quán tính Ji của khâu i: x2dm x y dm x z dm x dm i i i i i i x y dm y 2dm z y dm y dm i i T i i i i i i Ji = ri ri dm . (2.35) x z dm z y dm z 2dm z dm i i i i i i x dm y dm z dm dm i i i SV: Hoàng Đức Thắng - Phan Thanh Hải 37
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Nếu dùng Tenso quán tính Iij: I  x2 x x dm, ij ij  k i j (2.36) k Với các chỉ số i, j, k lấy lần lượt bằng các giá trị x i , yi, zi, đó là các trục của hệ toạ độ i, và  ij là kí hiệu Cronecke, thì ma trận Ji có thể biểu thị ở dạng sau: Ixx Iyy Izz I I m xi 2 xy xz i I I I xx yy zz Ixy Iyz mi yi J 2 i Ixx Iyy Izz I I m zi (2.37) xz yz 2 i mi xi mi yi mi zi mi T 1 Ở đây ri xi ,yi ,zi ,1 - bán kính vector biểu diễn trọng tâm của khâu thứ i trong hệ toạ độ i. Công thức (2.37) viết thành : 2 2 2 ki11 ki 22 ki33 2 2 k k xi 2 i12 i13 k2 k2 k2 2 i11 i 22 i33 2 ki12 ki 23 yi J 2 i 2 2 2 2 2 ki11 ki 22 ki33 k k zi (2.38) i13 i 23 2 xi yi zi 1 I jk Ở đây Kijk và j = 1,2,3 ; k = 1, 2, 3 mi T i r i xi , yi , zi ,1 - bán kính véc tơ biểu diễn trọng tâm của khâu thứ i trong hệ toạ độ i . Vậy, động năng của toàn cơ cấu tay máy bằng tổng đại số động năng của các khâu động : SV: Hoàng Đức Thắng - Phan Thanh Hải 38
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do n 1 n i i 1 n i i T   T   K Ki Tr UipJiUir qp qr  Tr UipJiUir qp qr (2.39) i 1 2 i 1 p 1r 1 2 i 1 p 1r 1 Lưu ý rằng, các ma trận J i (i=1,2,3, n) chỉ phụ thuộc vào sự phân bố khối lượng của khâu i trong hệ toạ độ i mà không phụ thuộc vào vị trí và vận tốc. Vì thế cần tính ma trận Ji chỉ 1 lần. II.3.4 .Thế năng tay máy: Thế năng Pi của khâu i: 0 0 1 Pi = - mi.g.ri = -mi.g. Ai ri (2.40) i=1,2, ,n Trong đó i 0 ri , ri - bán kính vec tơ biểu diễn trọng tâm của khâu i trong hệ toạ độ cơ bản . g - vector gia tốc trọng trường, g = ( 0, 0 , -g, 0) ( gia tốc trọng trường g = 9,8062 m/s2 ) Thế năng của toàn cơ cấu n khâu động : n n 0 i P P m g A ri . (2.41)  i  i i i 1 i 1 II.3.5.Mô hình động lực học tay máy: Để xây dựng mô hình động lực học tay máy dùng phương trình Lagrange bậc 2 (phương trình 2.16): d L L F , I = 1,2, ,n dt q Mi  q i Phương trình trên cho biểu thức tính động lực FMi. Đó là lực hoặc mô men tạo nên bởi nguồn động lực ở khớp động i để thực hiện chuyển động khâu i. Thay (2.39), (2.41) vào (2.16), cuối cùng ta có : SV: Hoàng Đức Thắng - Phan Thanh Hải 39
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do n j n j j n T T j FMi Tr U jkJ jU ji qk   Tr U jkmJ jU ji q k q m m jgU ji r j i=1,2 j 1 k 1 j 1 k 1 m 1 j 1 , ,n (2.42) Biểu thức (2.42) có thể viết gọn lại như sau : n n n FMi Dikqk   hikmq kq m ci (2.43) j 1 k 1 m 1 Hoặc là dưới dạng ma trận : FMi D(q)q h(q,q) c(q) (2.44) Trong đó : FM(t) – vector (nx1) lực động, tạo nên ở n khớp động : T FM (t) FM1(t), FM 2 (t), , FMn (t) (2.45) q(t) – vector (nx1) biến khớp : T q(t) q1(t),q2 (t), ,qn (t) (2.46) q(t) - vec tơ (nx1) tốc độ thay đổi biến khớp : T q(t) q1(t),q 2 (t), ,q n (t) (2.47) q(t) - vectơ (nx1) gia tốc biến khớp : T q(t) q1(t),q2 (t), ,qn (t) (2.48) D(q) – Ma trận (nxn) , có các phần tử Dik sau đây : n T Dik Tr U jk J jU ji i,k =1,2, . n (2.49) j max(i,k ) H(q, q ) – vectơ (nx1) lực ly tâm và Coriolit T h(q,q) (h1,h2 , ,hn ) n n hi   hikm q k q m i=1,2, ,n (2.50) k 1 m 1 SV: Hoàng Đức Thắng - Phan Thanh Hải 40
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do n T hikm Tr Uikm J jU ji i,k,m = 1,2, ,n (2.51) j max(i,k,m) C(q) – vec tơ (nx1) lực trọng trường. n T j c(q) c1,c2 , ,cn ci  mi gU ji rj (2.52) j 1 II.3.6. Động lực học của cơ cấu tay máy 2 khâu: Trong phần này dẫn ra ví dụ minh họa xây dựng mô hình động lực học của cơ cấu tay máy 2 khâu toàn khớp (hình 2.3) . Như đã chỉ trên hình 2.3 các trục Zi đều trùng phương với các trục khớp quay động. Khối lượng của các khâu tương ứng là m 1, m2; Bộ thông số DH của tay máy ghi trong bảng sau : Khâu i i Ai di Khớp * 1 i 0 1 0 R * 2 i 0 1 0 R Hình 2.3: Cơ cấu tay máy 2 khâu i 1 Các ma trận Ai ( i =1,2) sẽ là : SV: Hoàng Đức Thắng - Phan Thanh Hải 41
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do C1 S1 0 lC1 C2 S2 0 lC2 S C 0 lS S C 0 lS 0 A 1 1 1 1A 2 2 2 1 0 0 1 0 2 0 0 1 0 0 0 0 1 0 0 0 1 C12 S12 0 l(C12 C1) S C 0 l(S S ) 0 A 0A 1A 12 12 12 1 2 2 0 0 1 0 0 0 0 1 vẫn như trước đây dùng các kí hiệu sau : Ci = cosi ; S i= sini ; Cij = cos(i + j ) ; Sij = sin(i + j ) . Theo (2.30), ta có : 0 1 0 0 C1 S1 0 lC1 S1 C1 0 lS1 0 1 0 0 0 S C 0 lS C S 0 lC  A1 0 1 1 1 1 1 1 U11 D1 A1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 Tương tự, đối với U12, và U22 : 0 1 0 0 C12 S12 0 l(C12 C1) S12 C12 0 l(S12 S1) 0 1 0 0 0 S C 0 l(S S ) C S 0 l(C C )  A2 0 12 12 12 1 12 12 12 1 U 21 D2 A2 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 C2 S1 0 lC1 0 1 0 0 C2 S2 0 lC2 S12 C12 0 lS12 0 S C 0 lS 1 0 0 0 S C 0 lS C S 0 lC  A2 0 1 1 1 1 2 2 2 12 12 12 U 22 A1D2 A2 2 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 Theo (2.37) và giả thiết cả các thành phần mômen ly tâm quán tính đều bằng 0 , ta có : 2 2 2 2 1/ 3m1l 0 0 1/ 2m1l 1/ 3m2l 0 0 1/ 2m2l 0 0 0 0 0 0 0 0 J , J . 1 0 0 0 0 2 0 0 0 0 2 2 1/ 2m1l 0 0 m1 1/ 2m2l 0 0 m2 SV: Hoàng Đức Thắng - Phan Thanh Hải 42
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Trên cơ sở (2.49), ta có : 2 2  S1 C1 0 lS1 1/ 3m1l 0 0 1/ 2m1l C S 0 lC T T 1 1 1 0 0 0 0 T D11 T2 U11J1U11 T2 U 21J 2U 21 T2 U11  0 0 0 0 0 0 0 0 0 0 0 0 2 1/ 2m1l 0 0 m1  2 2  S12 C12 0 l(S12 C1) 1/ 3m2l 0 0 1/ 2m2l C S 0 l(C C ) 12 12 12 1 0 0 0 0 T 2 2 2 T2 U 21  1/ 3m1l 4 / 3m2l m2C2l 0 0 0 0 0 0 0 0 0 0 0 0 2 1/ 2m2l 0 0 m2  2 2  S12 C12 0 lS12 1/ 3m2l 0 0 1/ 2m2l C S 0 lC T 12 12 12 0 0 0 0 T D12 D21 Tr U 22 J 2U 21 T2 U 21  0 0 0 0 0 0 0 0 0 0 0 0 2 1/ 2m2l 0 0 m2  2 2 2 m2l ( 1/ 6 1/ 2 1/ 2C2 ) = 1/ 3m2l 1/ 2m2l . 2 2  S12 C12 0 lS12 1/ 3m2l 0 0 1/ 2m2l C S 0 lC T 12 12 12 0 0 0 0 T D22 Tr U 22 J 2U 22 T2 U 22  0 0 0 0 0 0 0 0 0 0 0 0 2 1/ 2m2l 0 0 m2  2 2 2 2 2 1/ 3m2l S12 1/ 3m2l C12 1/ 3m2l . Tính các số hạng trong biểu thức (2.50) đối với i = 1, ta có : 2 2        2  2 h1   h1km k m h111 1 2 h121 1 2 h122 2 h122 2 k 1 m 1 Và theo (2.51) tính các hệ số h ikm, rồi thay vào phương trình trên, ta được: 2  2 2   h1 = -1/2m2S2l 2 - m2S2l 1 2 . Tương tự đối với i = 2 SV: Hoàng Đức Thắng - Phan Thanh Hải 43
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do 2 2    2      2 2  2 h 2   h2km k m h211 1 h212 1 2 h221 2 1 h222 2 1/ 2m2S2l 1 k 1 m 1 Như vậy : 2 2  2   1/ 2m S l 2 m S l 1 2 H(,) 2 2 2 2 . 2 2  1/ 2m2S2l 1 Trên cơ sở (2.52), ta có : 1 2 c1 m1gU11 r1 m2gU21 r2 S1 C1 0 lS1 1/ 2 S12 C12 0 l(S12 S1) 1/ 2 C S 0 lC 0 C S 0 l(C C ) 0 m(0, g,0,0) 1 1 1 m (0, g,0,0) 12 12 12 1 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 1/ 2m1glC1 1/ 2m2 glC12 1/ 2m2 glC1 ; S12 C12 0 lS12 1/ 2 C12 S12 0 lC12 0 c2 m (0, g,0,0) = 2 0 0 0 0 0 0 0 0 0 1 m2 (1/ 2glC12 glC12 ) Vậy vector trọng trường sẽ là: c1 1/ 2m1glC1 1/ 2m2 glC12 m2 glC1 c( ) . c2 1/ 2m2 glC12 Cuối cùng ta có phương trinh động lực học của cơ cấu tay máy 2 khâu ở dạng sau :   Fm(t) = D()()(t) h(,) c() SV: Hoàng Đức Thắng - Phan Thanh Hải 44
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do 2 2 2 FM1 1/ 3m1l 4 / 3m2l m2C2l 1/ 3m2l 1/ 2m2l C2 2 2 2 + FM 2 1/ 3m2l 1/ 2m2l C2 1/ 3m2l 2 2  2   1/ 2m S l 2 m S l 1 2 1/ 2m glC 1/ 2m glC m glC 2 2 2 2 + 2 1 2 12 2 1 . 2 2 1/ 2m2 gl C12 2  1/ 2m2S2l 1 SV: Hoàng Đức Thắng - Phan Thanh Hải 45
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do CHƯƠNG III : THIẾT KẾ BỘ ĐIỀU KHIỂN TRƯỢT CHO TAY MÁY ROBOT 2 BẬC TỰ DO III.1.Hệ phi tuyến : III.1.1.Hệ phi tuyến là gì ? Để định nghĩa được rõ ràng một đối tượng hay hệ thống như thế nào được gọi là phi tuyến trước tiên ta nên định nghĩa lại hệ tuyến tính. Xét một hệ thống MIMO, viết tắt của nhiều vào / nhiều ra (Multi Inputs – Multi Outputs) với r tín hiệu vào u1(t), u2(t), , ur(t) và s tín hiệu ra y1(t) , y2(t) , , ys(t) . Nếu viết chung r tín hiệu đầu vào thành vectơ u1(t) u(t)  ur (t) y1(t) và s tín hiệu đầu ra thành y(t)  thì mô hình hệ thống được ys (t) quan tâm ở đây là mô hình toán học mô tả quan hệ giữa vectơ tín hiệu vào u(t) và tín hiệu ra y(t) , tức là mô tả ánh xạ T : u(t) y(t) . Ánh xạ này (Thường còn gọi là toán tử - operator) viết lại như sau : y(t) T(u(t)). nếu ánh xạ T thoả mãn : T(a1u1(t)) T(a2u2 (t)) a1T(u1(t)) a2T(u2 (t)) , (3.1) trong đó a 1 và a2 R, thì hệ đó được nói là tuyến tính. Tính chất (3.1) của hệ tuyến tính, trong điều khiển, còn được gọi là nguyên lý xếp chồng. Ví dụ : Xét 1 hệ gồm 1 lò xo c và 1 vật khối lượng m làm 1 ví dụ. Vật sẽ chuyển động trên trục nằm ngang dưới tác động của lực F (hình 3.1). SV: Hoàng Đức Thắng - Phan Thanh Hải 46
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do ms F cs m s Hình 3.1: Ví dụ về một đối tượng tuyến tính Nếu F được xem như là tín hiệu vào và quãng đường s mà vật đi được là tín hiệu ra (đáp ứng của hệ) thì theo các tiên đề cơ học của Newton, tác động vào vật và ngược hướng với F có hai lực cân bằng: Lực cản của lò xo F1 =cs trong trường hợp |s| tương đối nhỏ và F 2 = m scủa chuyển động. Với nguyên lý cân bằng lực ta có ánh xạ T : F(t) s(t) mô tả quan hệ vào / ra của hệ : ms cs F . (3.2a) Giả sử rằng dưới tác động của lực F 1 hệ có đáp ứng s 1 và của F2 thì từ : ms1 cs1 F1 ms2 cs2 F2 có ngay được m(a1 s1 a2 s2 ) c(a1s1 a2s2 ) a1F1 a2 F2 , trong đó a1 , a2 là những số thực tuỳ ý . Nói cách khác dưới tác động của lực a1F1 a2 F2 vật sẽ đi được một quãng đường là a1s1 a2s2 . Bởi vậy T thoả mãn (3.1) và do đó trong trường hợp |s| tương đối nhỏ và lực cản của lò xo được xác định gần đúng bằng công thức F1 = cs thì hệ thống là xo + vật là một hệ tuyến tính . SV: Hoàng Đức Thắng - Phan Thanh Hải 47
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do 3 Ngược lại, nếu lực cản lò xo lại được tính theo F 1 = cs +  s , với c và  là 2 hằng số, mà trong thực tế người ta vẫn sử dụng, thì quan hệ vào / ra của hệ sẽ là : ms cs s3 =F (3.2b) và khi đó (3.2b) không còn thoả mãn nguyên lý xếp chồng (3.1) 3 3 m(a1 s1 a2 s2 ) c(a1s1 a2s2 ) (a1s1 a2s2 ) a1F1 a2 F2 3 3 m(a1 s1 a2 s2 ) c(a1s1 a2s2 ) (a1s1 a2s2 ) a1F1 a2 F2 Nói cách khác, dưới tác động của lực a1F1 a2 F thì quãng đường của vật đi được không phải là a1s1 a2s2 . Vậy ở trường hợp này hệ có tính phi tuyến. III.1.2.Mô hình trạng thái và quỹ đạo trạng thái của Hệ phi tuyến: III.1.2.1.Mô hình trạng thái: Mô hình động của đối tượng, hệ thống phi tuyến được xây dựng từ quan hệ vào – ra qua việc thêm các biến x 1(t), x2 (t), , xn(t) gọi là biến trạng thái, sao cho quan hệ giữa vector tín hiệu ra y(t) với n biến này và tín hiệu vào u(t) chỉ còn lại thuần tuý là một quan hệ đại số. Những biến trạng thái này, về mặt ý nghĩa vật lý, là những đại lượng mà sự thay đổi của nó sẽ quyết định tính chất động học của đối tượng. Ví dụ 1 : Từ mô hình (3.1) của đối tượng lò xo + vật, nếu thêm biến trạng thái x1 = s , x2 = s sẽ có: x s (1 0) 1 (1 0).x 0.F x2 và đấy là một phương trình đại số. Ngoài ra còn có phần các phương trình vi phân bao gồm : x1 x 2 Suy ra từ định nghĩa về x1 , x2 và : SV: Hoàng Đức Thắng - Phan Thanh Hải 48
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do c 1 x x F 2 m 1 m thu được bằng cách thay trực tiếp x1, x2 vào phương trình (3.1). Viết chung hai phương trình vi phân trên lại với nhau sẽ được : . 0 1 0 . x1 x c x 1 F . 0 x 2 m m Nói chung, một hệ phi tuyến SISO có quan hệ vào – ra giữa tín hiệu vào u(t) và ra y(t) dạng: . y(n) f (y(n 1) , ,y,y,u) Trong đó ký hiệu y(k) chỉ đạo hàm bậc k của y(t), tức là d k y y(k ) , dt k thì với các biến trạng thái được định nghĩa như sau : (n 1) x1 y, x2 y, x3 y, , xn y hệ sẽ có mô hình hai phần: phần các phương trình vi phân bậc nhất x1 x2  xn 1 xn xn f (xn , , x2 , x1,u) và phương trình đại số : Y=x1 Tổng quát lên thì một hệ phi tuyến, sau khi định nghĩa các biến trạng thái x1(t) , x2(t) , , xn(t), sẽ mô tả bởi : - Mô hình trạng thái tường minh autonom . x (t) x f (x,u) 1 trong đó x(t)  y g(x,u) xn (t) SV: Hoàng Đức Thắng - Phan Thanh Hải 49
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do - Mô hình trạng thái tường minh không autonom . x f (x,u,t) , y g(x,u,t) - hoặc mô hình trạng thái không tường minh . f (x, x,u,t) 0 g(x,u,y,t) 0 III.1.2.2.Quỹ đạo trạng thái : Từ phương trình trạng thái mô tả hệ thống với một tín hiệu đầu vào u(t) xác định cho trước và với một điểm trạng thái ban đầu x 0 x(0) cũng cho trước ta sẽ có được nghiệm x(t) mô tả sự thay đổi trạng thái hệ thống theo thời gian dưới tác động của kích thích u(t) đã cho. Biểu diễn x(t) trong không gian n chiều Rn (còn gọi là không gian trạng thái) như một đồ thị phụ thuộc tham số t có mũi tên chỉ chiều tăng của t ta được một quỹ đạo trạng thái (hình 3.2a). Tập tất cả các quỹ đạo trạng thái ứng với một tín hiệu đầu vào u(t) cố định nhưng với những điểm trạng thái ban đầu x 0 khác nhau được gọi là họ các quỹ đạo trạng thái (hình 3.2b) . Hình 3.2a. Quỹ đạo trạmg thái của hệ có 3 biến trạng thái x1 SV: Hoàng Đức Thắng - Phan Thanh Hải 50
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do x2 Hình 3.2b. Họ các quỹ đạo trạng thái của hệ có 2 biến trạng thái III.1.3. Điểm cân bằng và điểm dừng của hệ thống: III.1.3.1. Điểm cân bằng: Định nghĩa 1: Một điểm trạng thái được gọi là điểm cân bằng x e (equilibrium point) nếu như khi đang ở điểm trạng thái và không có một x e tác động nào từ bên ngoài thì hệ sẽ nằm nguyên tại đó. Căn cứ theo định nghĩa như vậy thì điểm cân bằng của hệ thống x e phải là nghiệm của phương trình: dx f (x,u,t) 0 dt u 0 Như vậy điểm cân bằng là điểm mà hệ thống sẽ nằm yên tại đó, tức dx là trạng thái của nó sẽ không bị thay đổi ( 0 ) khi không có sự tác đng dt từ bên ngoài (u 0 ). III.1.3.2.Điểm dừng của hệ : Định nghĩa 2: Một điểm trạng thái x d được gọi là điểm dừng của hệ thống nếu như hệ đang ở điểm trạng thái x dvà với tác động u(t) u cốd định, không đổi cho trước, thì hệ sẽ nằm nguyên tại đó. Rõ ràng là điểm dừng theo định nghĩa vừa nêu sẽ là nghiệm của : dx f (x,u ,t) f (x,u ,t) 0 d d dt u ud SV: Hoàng Đức Thắng - Phan Thanh Hải 51
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do trong đó ud là đã cho trước. III.1.3.3 Tính ổn định tại một điểm cân bằng: Định nghĩa 3 : Một hệ thống được gọi là ổn định (tiệm cận) tại điểm cân bằng x e nếu như có một tác động tức thời (chẳng hạn như nhiễu tức thời) đánh bật hệ ra khỏi x e và đưa tới điểm x 0 thuộc một lân cận nào đó của x e thì sau đó hệ có khả năng tự quay về được điểm cân bằng x bane đầu. Theo định nghĩa trên thì ta có thể nhận biết được hệ có ổn định hay không tại một điểm cân bằng thông qua dạng họ các đường quỹ đạo trạng thái của nó. Nếu hệ ổn định tại một điểm cân bằng x e nào đó thì mọi đường quỹ đạo trạng thái x(t) xuất phát từ một điểm x 0 thuộc lân cận của x e đều phải kết thúc tại x e a) Miền ổn định O b) xe xe Hình 3.3. a) Điểm cân bằng ổn định b) Điểm cân bằng không ổn định Chú ý rằng tính ổn định của hệ phi tuyến chỉ có ý nghĩa khi đi cùng với điểm cân bằng x e . Có thể hệ sẽ ổn định tại điểm cân bằng này, song lại không ổn định ở điểm cân bằng khác. Điều này cũng khác so với khái niệm ổn định ở hệ tuyến tính. Vì hệ tuyến tính thường chỉ có một điểm cân bằng là gốc toạ độ (x e =0 ) nên khi hệ ổn định tại 0 , người ta cũng nói thêm luôn một cách ngắn gọn là hệ ổn định . SV: Hoàng Đức Thắng - Phan Thanh Hải 52
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Ngoài ra, do khái niệm ổn định ở hệ phi tuyến bị gắn với lân cận điểm cân bằng x e nên cũng có thệ mặc dù hệ ổn định tại điểm cân bằng x e song với một lân cận quá nhỏ thì sẽ không có định nghĩa sử dụng. Nói cách khác, về mặt ứng dụng, nó được xem như không ổn định. Bởi vậy, đối với hệ phi tuyến, việc xác định xem hệ có ổn định tại điểm cân bằng hayx e không là chưa đủ mà còn phải chỉ ra miền ổn định của nó tại x e , tức là phải chỉ ra được lân cận O của x e sao cho hệ có khả năng tự quay về được x etừ bất kì một điểm x 0 nào đó thuộc O (hình 3.3). Miền ổn định O càng lớn thì tính ổn định của hệ tại x e càng tốt . Nhiệm vụ đầu tiên của bộ điều khiển là phải giữ cho hệ thống ổn định. Nếu như ban đầu đối tượng không ổn định, tức là khi có nhiễu từ bên ngoài tác động đưa nó ra khỏi điểm làm việc và nó không có khả năng tự quay về thì bộ điều khiển phải tạo ra tín hiệu điều khiển dẫn đối tượng quay trở về điểm làm việc ban đầu. III.1.4 Tiêu chuẩn ổn đinh Lyapunov : Một trong những điều kiện, hay tiêu chuẩn chất lượng đầu tiên mà bộ điều khiển cần phải mang đến được cho hệ thống là tính ổn định. Tại sao lại như vậy ? Từ khái niệm về tính ổn định của hệ thống tại một điểm cân bằng đã được nêu trong định nghĩa 3 ta thấy rõ nếu một hệ quá nhạy cảm với tác động nhiễu đến nỗi chỉ một tác động tức thời không mong muốn rất nhỏ đã làm cho hệ bị bật ra khỏi điểm cân bằng (hoặc điểm làm việc) mà sau đó hệ không có khả năng tự tìm về điểm cân bằng ban đầu thì chất lượng của hệ không thể gọi là tốt được. Bởi vậy, kiểm tra tính ổn định của hệ (tại một điểm cân bằng) cũng như miền ổn định O tương ứng phải là công việc đầu tiên ta phải tiến hành khi phân tích hệ thống. Tiêu chuẩn Lyapunov là một công cụ hữu ích giúp ta thực hiên được điều đó. Định nghĩa 4 : Một hệ thống có mô hình không kích thích : dx ~ f (x,u,t) f (x,t) dt u 0 (3.3) SV: Hoàng Đức Thắng - Phan Thanh Hải 53
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do với một điểm cân bằng là gốc toạ độ 0 , được gọi là : a) Ổn định Lyapunov tại điểm cân bằng 0nếu với  bất0 kì bao giờ cũng tồn tại  phụ thuộc  sao cho nghiệm x(t) của (3.3) với x(0) x0 thoả mãn : x 0  x(t) ,t 0. b) Ổn định tiệm cận Lyapunov tại điểm cân bằng 0 nếu với  0bất kì bao giờ cũng tồn tại  phụ thuộc  sao cho nghiệm x(t) của (3.3) với x(0) x0 thoả mãn : lim x(t) 0 . t III.1.4.1.Tiêu chuẩn Lyapunov: Để làm quen và tiếp cận tiêu chuẩn Lyapunov ta hãy bắt đầu từ hệ bậc 2 có mô hình trạng thái autonom khi không bị kích thích : dx x1 f (x1, x2 ) với x (3.4) dt x2 Hệ trên được giả thiết là cân bằng tại gốc toạ độ 0 . Hình 3.4 : Minh họa khái niệm ổn định Lyapunov: Hình 3.4 minh hoạ cho định nghĩa 4 về tính ổn định Lyapunov tại 0 đã gợi ý cho ta một hướng khá đơn giản để xét tính ổn định cho hệ (3.4) tại 0 . Chẳng hạn bằng cách nào đó ta đã có được họ các đường cong khép kín v bao quanh gốc toạ độ 0 . Vậy thì để kiểm tra hệ có ổn định tại 0hay không ta chỉ cần kiểm tra xem quỹ dạo pha x(t) , tức là nghiệm của (3.4) đi SV: Hoàng Đức Thắng - Phan Thanh Hải 54
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do từ điểm trạng thái đầu x0 cho trước nhưng tuỳ ý nằm trong miền bao bởi một trong các đường cong khép kín đó, có cắt các đường cong v này theo hướng từ ngoài vào trong hay không ( hình 3.5). - Nếu x(t) không cắt bất cứ một đường cong họ v nào theo chiều từ trong ra ngoài thì hệ sẽ ổn định tại 0 . - Nếu x(t) cắt mọi đường cong họ v theo chiều từ ngoài vào trong thì hệ sẽ ổn định tiệm cận tại 0 . Hình 3.5: Một gợi ý về việc kiểm tra tính ổn định của hệ tại O Rõ ràng là cần và đủ để quỹ đạo pha x(t )của hệ không cắt bất cứ một đường cong khép kín thuộc họ v theo chiều từ trong ra ngoài là tại điểm cắt đó, tiếp tuyến của x(t) phải tạo với vector v , được định nghĩa là vector vuông góc với đường cong đó theo hướng từ trong ra ngoài, một góc không nhỏ hơn 900. Nói cách khác, hệ sẽ ổn định tại 0 nếu như có được điều kiện: dx dx 0 . .cos = T (3.5) v dt v dt tại mọi giao điểm của x(t) với các đường cong thuộc họ v. Vấn đề còn lại là làm thế nào có được các đường cong v sao cho việc kiểm tra điều kiện (1.48) được thuận tiện. Câu trả lời là sử dụng hàm xác định dương V( x )được định nghĩa như sau : SV: Hoàng Đức Thắng - Phan Thanh Hải 55
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Định nghĩa 5 : Một hàm thực nhiều biến, có thể không dừng V(0,t), được gọi là hàm xác định dương nếu : a) V(0,t) 0 b) Tồn tại hai hàm một biến, dừng  1(a) và  1(b) liên tục, đơn điệu tăng với  1(0)  2 (0) 0 sao cho : 0  1( x ) V(x,t)  2 ( x ) với mọi x 0 (3.6) Hàm V(x,t) sẽ xác định dương trong toàn bộ không gian trạng thái nếu còn có : lim 1(a) => lim V(x,t) . a x Định lý 1 : Hệ phi tuyến (có thể không autonom) cân bằng tại gốc toạ độ và khi không bị kích thích thì được mô tả bởi hình : dx f (x,t) (3.7) dt sẽ ổn định Lyapunov tại 0 với miền ổn định O nếu : a) Trong O tồn tại một hàm xác định dương V(0,t) . b) Đạo hàm của nó tính theo mô hình (1.51) có giá trị không dương trong O, tức là : dV V V f (x,t) 0 với mọi x O . (3.8) dt t x Định lý 2: Hệ phi tuyến (có thể không autonom) cân bằng tại gốc toạ độ và khi không bị kích thích thì được mô tả bởi mô hình. dx f (x,t) (3.9) dt sẽ ổn định tiệm cận Lyapunov tại 0 với miền ổn định O nếu : a)Trong O tồn tại một hàm xác định dương V(x,t) . SV: Hoàng Đức Thắng - Phan Thanh Hải 56
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do b) Đạo hàm của nó tính theo môt hình (1.51) có giá trị âm trong O với x 0 , tức là : dV V V f (x,t) 0 với mọi x O và x 0 . (3.10) dt t x III.1.4.2.Tiêu chuẩn Lyapunov phục vụ thiết kế bộ điều khiển: Ngoài việc kiểm tra tính ổn định, tiêu chuẩn Lyapunov còn được sử dụng để thiết kế bộ điều khiển ổn định đối tượng phi tuyến. Chẳng hạn đối tượng có mô hình : dx f (x,u) dt và được điều khiển bằng bộ điều khiển phản hồi trạng thái r(x)  u dx x f (x,u) dt r(x) Hình 3.6: Ứng dụng tiêu chuẩn Lyapunov để thiết kế bộ điều khiển Vậy hệ kín khi không bị kích thích ( 0 ) sẽ có mô hình : dx f (x,r(x)) dt Gọi V(x) là hàm xác định dương thích hợp, khi đó để hệ kín ổn định tiệm cận với miền ổn định là O thì bộ điều khiển cần tìm r(phảix) thoả mãn : V L V f (x,r(x)) 0 với mọi x 0 , x O (3.11a) f x V Và f (x,r(x))=0chỉ khi x 0 (3.11b) x III.2.Bậc tương đối của hệ phi tuyến: SV: Hoàng Đức Thắng - Phan Thanh Hải 57
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Bậc tương đối của hệ SISO: Để dễ tiếp cận tới khái niệm bậc tương đối ta xét trường hợp đặc biệt với đối tượng tuyến tính, mô tả bằng hàm truyền đạt hợp thức chặt (strickly proper): m G(s) = b0 b1 bm s (3.12) s n a0 a1 an s Khi đó bậc tương đối được hiểu là hiệu r = (n-m) ≥1 Giả sử rằng đối tượng trên, bên cạnh hàm truyền đạt (3.12) còn có mô hình tương đương trong không gian trạng thái : d x A x bu n nxn nx1 1xn dt x R ,A R ,b R ,c R (3.13) y cT x Vậy thì do G(s)= cT(sI-A)-1 b Ta có : T Lím srG(s)= bm lim sr [c (sI A) 1 b ] = bm s s an an T k c A b Lím = bm s  k 1 r k 0 s an Hơn nữa 1 Lím = 0 khi k > r-1 s sk 1 r nên chuỗi trên trở thành tổng của hữu hạn r phần tử đầu tiên T r 1 c kb Lím A = bm s  k 1 r k 0 s an SV: Hoàng Đức Thắng - Phan Thanh Hải 58
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Từ đây, để vế trái bằng giá trị hữu hạn thì cần và đủ là : T k 0 khi 0 k r - 2 c A b = (3.14) 0 khi k r -1 Nói cách khác, bậc tương đối r = n-m còn có thể được xác định trực tiếp từ mô hình trạng thái (3.13) của hệ theo công thức (3.14). Chuyển sang hệ phi tuyến và với sự gợi ý của công thức tính (3.14), khái niệm bậc tương đối của hệ ALI có 1 tín hiệu vào, một tín hiệu ra, được định nghĩa như sau : Định nghĩa 6: Cho hệ SISO với cấu trúc ALI : d x f (x) h(x)u dt (3.15) y g(x) Bậc tương đối tại điểm trạng thái x của hệ là số tự nhiên r mà trong lân cận x thoả mãn : h 0 khi 0 k r-2 Lh Lfg(x) (3.16) 0 khi k r-1 Có thể thấy được ngay rằng với f(x)= Ax , H(x)= b , g(x)=c Tx , hai công thức (3.14) và (3.16) sẽ đồng nhất, vì : h g(x)= cTATx L k g(x)=cTAkb L f h L f Ví dụ : Xét hệ Val der Pol có mô hình trạng thái như sau : khi đó thì do : g 0 Lhg(x)= h(x) [1 0 ] = 0  x 1 (Lf g)  g LhLfg(x) = h(x) f h(x)  x x x SV: Hoàng Đức Thắng - Phan Thanh Hải 59
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do  x 2 0 0 = 1 0 2 1 0 =1 ≠ 0  x ax2 (1 bx1 ) x1 1 1 Bậc tương đối của hệ bằng 2 ( tại mọi x ). Tuy nhiên, cũng cần phải để ý rằng hệ phi tuyến (3.15) có thể có bậc tương đối khác nhau ở những điểm trạng thái khác nhau. Ngoài ra, khác với hệ tuyến tính, không phải ở bất cứ điểm trạng thái x nào trong không gian trạng thái, hệ phi tuyến phẳng có bậc tương đối. Chẳng hạn, hệ sẽ không có bậc tương đối tại điểm trạng thái x0 mà trong lân cận của nó có : h Lhg( x )≠0,LhLfg(x ) ≠0, ,LhLf g(x ) ≠ 0 , III.3.Tính động hoc không: Rất nhiều khái niệm sử dụng trong hệ phi tuyến được chuyển thể từ hệ tuyến tính, chẳng hạn khái niệm bậc tương đối, hệ thụ động, cũng như vậy là tính động học không (zero dynamic). Do đó, để dễ tiếp cận tới khái niệm này, ta nên bắt đầu từ hệ tuyến tính. Xét hệ phi tuyến SISO có mô hình trạng thái : d x f (x) h(x)u dt (3.17) y g(x ) Tính động học không (zero dynamic ) của hệ (3.17) được định nghĩa như sau : Định nghĩa 7 : Nếu hệ (3.17) có ít nhất một điểm trạng thái đầu x 0≠ 0 và ứng với nó là tín hiệu điều khiển u 0(t) sao cho tín hiệu đầu ra y(t) đồng nhất bằng không thì hệ được gọi là có tính động học không (zero dynamic). Ta có thể thấy được là để hệ có tính động học không thì cần thiết phải có g(x0) = 0. Giả sử rằng hệ (3.17) có bậc tương đối là r, tức là : SV: Hoàng Đức Thắng - Phan Thanh Hải 60
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do 0 nÕu 0 k r-2 k LhLf g(x) = (3.18) 0 nÕu k= r -1 Khi đó, với phép đổi trục toạ độ vi phôi : g(x) z1   r 2 Lg(x) zr 1 f r 1 z zr m(x) với Lhmk (x) = 0 , k=r+1 , , n Lg(x) f zr 1 mr 1(x)   zn mn(x) hệ (1.18) đã cho sẽ được đưa về dạng chuẩn z1 z2   zr dz zr 1 d a(z) b(z)u zr , y = z1 (3.19) dt dt (z) zr 1 c1   zn cn r (z) Trong đó r 1 r 1 1 A(z) = Lf g(m (z) ) , b(z) =Lh Lf g(m (z) ) , ci (z) = Lh mr+1( m 1(z)) sử dụng kí hiệu : c1(z)  zr 1 zr 1 z = v ới  =  ,   v à c(z)   z n zn cn r(z) SV: Hoàng Đức Thắng - Phan Thanh Hải 61
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do thì mô hình (3.19) được viết thành z2 z1   dz d zr zr 1 , y = z1 (3.20) dt dt a(,) b(,)u zr  c(,n) Giả sử rằng hệ (3.17) có tính động học không ứng với trạng thái đầu x0 ≠0 và tín hiệu điều khiển u0(t) thích hợp. Vậy thì từ y(t) = z1(t) = 0 ta suy ra được : z1(t) = = zr(t) =0 và do đó là  = 0. Điều này dẫn đến : a(0,) a(0 , ) + b(0 ,  )u0 = 0 u0 (t) = - (3.21a) b(0,) d c(0,) (3.21b) dt Đó cũng là hai phương trình phân tích tính động học không của hệ (3.17) thông qua mô hình tương đương (3.20) của nó. Điều kiện để có phương trình (3.21b) là hệ (3.17) phải có bậc tương đối r nhỏ hơn n (r < n). Từ =0 cũng như phép biến đổi trục toạ độ (3.18) và 2 phương trình (3.21) ta thấy, ở chế độ động học không, quỹ đạo trạng thái x(t) phải thoả r 1 mãn : g(x ) =Lfg(x )= =L f g(x ) = 0 . Nói cách khác x(t) của động học không sẽ chỉ nằm trong đa tạp (hình 3.7) n r 1 K = { xR |g(x )=Lfg(x )= = L f g(x ) = 0 } (3.22) SV: Hoàng Đức Thắng - Phan Thanh Hải 62
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do  x(t) 1 2 Hình 3.7: Quỹ đạo trạng thái của Hệ phi tuyến, khi đang ở chế độ Động học không, luôn nằm trong đa tạp K. Tuy rằng nằm trong đa tạp K, song việc quỹ đạo x(t) ở chế độ động học không (ứng với tín hiệu điều khiển u 0(t) thích hợp) có tiến về gốc toạ độ 0 hay không thì chưa được đảm bảo và điều này không được quyết định bởi hệ phi tuyến (3.17) có ổn định hay không. Nó chỉ có thể tiến về 0n như hệ (3.21b) là ổn định tiệm cận Lyapunov, tức là phải tồn tại 1 hàm xác định dương Q() sao cho : Q c(0,) 0 khi  0  III.4.Thiết kế bộ Điều khiển trượt cho tay máy: III.4.1.Điều khiển trượt: Hệ phi tuyến có mô hình x f (x,u) f (x) g(x)u (3.23) y h(x) Trong đó y là tín hiệu đầu ra, u là tín hiệu đầu vào, x = [x 1, x2, , T T xn] là vector trạng thái của hệ, f(x) = [f 1(x), f2(x), , fn(x)] , g(x) = [g1(x), T g2(x), , gn(x)] SV: Hoàng Đức Thắng - Phan Thanh Hải 63
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Hệ phi tuyến có bậc tương đối là p nếu: p d h(x) p p 1 L f h(x) Lg L f h(x)u dt p (3.24) p 1 i Lg L f h(x) 0; Lg L f h(x) 0,i 1,2,3, , p 2 Sơ đồ điều khiển: y , y , , y (n) r r r y, y, , y (n) x f (x) g(x)u SMC y h(x) III.4.1.1.Trường hợp bậc tương đối của hệ bằng bậc của hệ p=n: Để có thể thiết kế được bộ điều khiển thì hệ (3.23) phải tồn tại mặt trượt. Hệ (3.23) có mặt trượt S khi thoả mãn: n 1 (i)  S e i e (3.25) i 1 (n 1)  A(S) 1 1S n 1S là đa thức Hurwitz để có: lime(t) 0 (3.26) t  S(0) = 0 (3.27) Điều kiện để (3.23) trượt về điểm cân bằng là phải thoả mãn điều kiện trượt. Điều kiện trượt được xây dựng trên cơ sở đảm bảo hệ kín ổn định tiệm cận, có nghĩa là cho hệ trong hình trên tồn tại 1 hàm Lyapunov. Giả sử hệ có hàm Lyapunov có dạng sau: 1 V (x,t) S 2 (3.28) 2 là hàm xác định dương. Đạo hàm của nó có dạng sau: SV: Hoàng Đức Thắng - Phan Thanh Hải 64
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do dV SS (3.29) dt Hệ (3.23) ổn định tiệm cận khi (3.29) là hàm có dấu xác định âm: S 0, S 0 SS 0 (3.30)  S 0, S 0 Như vậy S phải trái dấu với S, do vậy ta có: S Kh(S) (3.31) h(S) cùng dấu với S do vậy để thoả mãn điều kiện trượt ta có thể chọn hàm h(S) có các dạng sau: hàm dấu Sig(S), hàm bão hoà Saturation(S), hàm h(S)=Tan(S) Theo (3.25) ta có: n   (i)   (n) S e i 1e e 1e ne Kh(S) (3.32) i 2 Ta có: (n) (n) n n 1 e yr (L f h(x) Lg L f h(x)u) (3.33) (n) n n 1 Do vậy: e 1e n yr n (L f h(x) Lg L f h(x)u) Kh(S) (3.34) Tín hiệu điều khiển tìm được: (n) n Kh(S) e  1e n yr n L f h(x) u(t) n 1 (3.35) n Lg L f h(x) III.4.1.2. Trường hợp bậc tương đối của hệ p<n Hệ (3.23) phải thoả mãn động học không. Xây dựng mặt trượt : p 1 (i)  S e i e (3.36) i 1 ( p 1)  A(s) 1 1S  p 1S là đa thức Hurwitz, để có lime(t) 0 (3.37) t  S(0) 0 , mặt trượt phải đi qua gốc toạ độ và thoả mãn điều kiện trượt. SV: Hoàng Đức Thắng - Phan Thanh Hải 65
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Hoàn toàn tương tự như trong trường hợp trên, ta xây dựng hàm Lyapunov có dạng sau: 1 V S 2 xác định dương 2 V SS xác định âm Ta có: S Kh(S)  ( p) S e 1e  p 1e e yr (t) y(t) ( p) ( p) ( p) (3.38) e yr (t) y (t) ( p) p ( p 1) y (t) L f h(x) Lg L f h(x)u  ( p) p ( p 1) S e 1e  p 1 yr (t)  p 1 (L f h(x) Lg L f h(x)u) Kh(S) Tín hiệu điều khiển: ( p) p Kh(S) e 1e  p 1 ( yr L f h(x)) u(t) ( p 1) (3.39)  p 1Lg L f h(S) III.4.2. Thiết kế bộ điều khiển trượt cho tay máy n bậc tự do: Mô hình động lực học của tay máy:  H (q)q h(q, q) (3.40) với H(q) là ma trận quán tính xác định dương, đối xứng. Chúng ta giả sử rằng các giá trị ước lượng Hˆ (q )và hˆ(q, q )quan hệ với giá trị thực H (q) và h(q, q) bởi bất đẳng thức sau: Hˆ 1 (q)H (q)  (q) (3.41) ˆ và h(q, q) h(q, q) hmax (q, q) (3.42) với  (q) và hmax (q, q) là những hàm đã biết. Viết lại biểu thức động lực học dưới dạng: q f (q, q) B(q) (3.43) SV: Hoàng Đức Thắng - Phan Thanh Hải 66
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Với f (q, q) H 1 (q)h(q, q) (3.44) B(q) H 1 (q) (3.45) Nhiệm vụ của điều khiển là tìm mô men thích hợp τ sao cho vector vị trí q của tay máy bám theo quỹ đạo mong muốn qd. Chúng ta định nghĩa sai lệch trạng thái e và mặt trượt như sau: e qd q (3.46) S Ce e;C C T 0 (3.47) Rõ ràng rằng S=0 thì q(t) qq (t) . Quả thực với S=0 ta có thể viết lại như sau: S 0 Ce e 0 e Ce Như vậy hệ thống ổn định tiệm cận nếu có e = 0 và theo đó điều kiện bám q(t) qq (t) sẽ được đảm bảo. Do vậy vấn đề điều khiển là phải tìm mô men τ thích hợp sao cho vector trạng thái của hệ thống có thể bám được trên mặt trượt. Hay phải tìm τ thỏa mãn điều kiện trượt. Điều kiện trượt có thể xác định theo tiêu chuẩn Lyapunov. Chúng ta định nghĩa hàm Lyapunov như sau: 1 V S T S 0 (3.48) 2 Đạo hàm của (3.48) có dạng: V S T S (3.49) Như vậy, nếu V 0 thì với V 0 dẫn tới S 0 và e 0 Do vậy, điều kiện đủ của điều kiện trượt là: S T S 0 (3.50) SV: Hoàng Đức Thắng - Phan Thanh Hải 67
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Khi đó điều kiện trượt đảm bảo cho hệ kín ổn định toàn cục, tiệm cận và điều kiện bám được thực hiện mặc dù mô hình không chính xác, nhiễu, Nếu điều kiện trượt có thể thỏa mãn theo đó: n T  2 S S S 0; 0; S  Si (3.51) i 1 Tiếp đó, mặt phẳng trượt S=0 sẽ đạt được với thời gian giới hạn nhỏ hơn T0 ở đó: 1 T S(q(0)) (3.52) 0 2 Biểu thức trên được chứng minh như sau: Từ (29) ta có: S T S (3.53) S 1 Thay V S T S và S (2V ) 2 vào (3.53) sau đó tích phân hai vế với t=0→treach , S(q(treach))=0 ta có: treach treach  1 S(q(0)) S(q(0)) V 1 2 1 dt V  treach treach T0 (3.54) 2 0 (2V ) 2 0 2 2 Bây giờ chúng ta tìm đầu vào bộ điều khiển τ thỏa mãn điều kiện trượt. Lấy đạo hàm biểu thức (3.47) ta có:  S Ce q qd (3.55) Thay biểu thức (3.39) vào ta có:  S Ce f (q, q) B(q) qd (3.56) Do đó tín hiệu điều khiển có dạng ˆ 1  B  eq KSgn(s) (3.57) SV: Hoàng Đức Thắng - Phan Thanh Hải 68
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do với:  q Ce fˆ(q, q) eq (3.58) T Sgn(s) Sgn(s1 ), Sgn(s2 ), , Sgn(sn ) K>0, K là ma trận khuyếch đại nxn. Ma trận khuyếch đại K phải chọn đủ lớn để điều kiện trượt được thỏa mãn mặc dù có tham số không rõ, nhiễu, Trong trường hợp ước lượng chính xác Bˆ B, fˆ f thì điều kiện trượt được viết lại như sau: S T S S T KSgn(s) S (3.59) Nếu chọn K I;  (3.60) m m T  2 và S S  S   Si   Si  S S (3.61) i 1 i 1 thì chế độ trượt xảy ra. Ta nhận thấy rằng, đầu vào điều khiển được gián đoạn qua s(t) như cho ở biểu thức (3.57). Hiện tượng chattering xảy ra. Bởi vì trong thực tế, sự chuyển đổi là không lý tưởng. Trong trường hợp sai số ước lượng là không đủ nhỏ thì việc chọn K là không đơn giản như biểu thức trên. Trong trường hợp đó S cho dưới dạng:  ˆ 1 ˆ 1 S qd Ce f (q, q) B(q)B  eq B(q)B KSgn(s) (3.62) đặt f fˆ ( f fˆ); R B(q)Bˆ 1 dẫn tới:  ˆ S (R I ) eq ( f f ) RKSgn(s) (3.63) Từ đây, điều kiện trượt là: T  T  T S S S (R I ) eq ( f f ) RKSgn(s) S Sgn(s) (3.64) Do vậy, nếu chọn K để: T T ˆ T S RKSgn(s) S (R I ) eq ( f f ) S Sgn(s) (3.65) SV: Hoàng Đức Thắng - Phan Thanh Hải 69
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do thì điều kiện trượt như ở trên SS 0 được thỏa mãn và điều kiện trượt đạt được. R 1 BˆB 1  (q) (3.66) Từ biểu thức (3.41) và (3.42) ta có bất đẳng thức: R 1 BˆB 1  (q) (3.67) 1 ˆ ˆ ˆ ˆ R ( f f ) B(h h) B hmax (3.68) Từ đây, ta có thể chọn ma trận K thoả mãn điều kiện trượt như sau: ˆ K (1  )I  eq B hmax I (3.69) III.4.3. Ứng dụng Điều khiển trượt cho tay máy Robot 2 bậc tự do: III.4.3.1. Phương trình động lực học tay máy hai bậc tự do toàn khớp quay: Bộ thông số tay máy: m1 = m2 = 1 kg l1 = l2 = 1 m Phương trình động lực học: F h h  c c  g 1 11 12 1 11 12 1 1 (3.70)   F2 h21 h22  2 c21 c22  2 g 2 Trong đó F1, F2 là lực được tạo ra ở các khớp động, ma trận H là ma trận xác định dương và đối xứng, ma trận C là ma trận lực ly tâm, G là ma trận lực trọng trường. Giá trị của các ma trận khi thay giá trị được xác định như sau: - Ma trận H: h11 5 3cos 2 h12 h21 1 3/ 2cos 2 (3.71) h22 1 SV: Hoàng Đức Thắng - Phan Thanh Hải 70
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do - Ma trận C:  c11 3 2 sin 2 c 3/ 2 sin 12 2 2 (3.72)  c21 31 sin 2 c22 0 - Ma trận G: g 15cos 15cos(  ) 1 1 1 2 (3.73) g 2 15cos(1  2 ) III.4.3.2. Mô hình động lực học tay máy hai bậc tự do: Chúng ta đặt các biến trạng thái là tín hiệu góc quay và vận tốc của các khớp tay máy: Khớp 1: x11 1  x12 x11 1 (3.74)  x12 1 Khớp 2: x21  2  x22 x21  2 (3.75)  x22  2 Tín hiệu vào u: u1 F1 (3.76) u2 F2 Từ biểu thức (3.70) ta có:   1 1 F1 1 c11 c12 1 1 g1 H H H (3.77)  F c c  g 1 2 21 22  2 2 trong đó H-1 là ma trận nghịch đảo của ma trận H. Tính H-1 và kết hợp tất cả các phương trình trên và thay vào (3.77) để tính được: SV: Hoàng Đức Thắng - Phan Thanh Hải 71
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Khớp 1 x x 11 12 1 3 (3.78) x12 (u1 x 22 (2 x12 x 22 )sin x 21 15(cos x11 cos x 21 )) 4 9 4 cos2 2 x 21 3 3 (1 cos x )(u x 2 sin x 15 cos(x x )) 2 21 2 2 12 21 11 21 Khớp 2: x x 21 22 (3.79) 1 3 3 x22 ( (1 cosx21)(u1 x22 (2x12 x22 )sin x21 15(cosx11 cosx21))) 4 9 4cos2 2 2 x21 3 (5 3cosx )(u x2 sin x 15cos(x x ))) 21 2 2 12 21 11 21 III.4.3.3. Thiết kế bộ điều khiển trượt cho tay máy 2 bậc tự do: Xác định bậc tương đối cho khớp 1 và khớp 2: Từ phương trình trạng thái của các khớp (3.78), (3.79) và biểu thức (3.24) ta có được ngay là trong trường hợp này là p=n hay bậc tương đối của từng khớp p=2. Xây dựng mặt trượt cho từng khớp: S e  e 1 1 1 1 (3.80) S 2 e2 2e2 ở đây: e x x 1 d1 r1 (3.81) e2 xd 2 xr 2 1, 2 là những số thực dương. Điều kiện để xảy ra chế độ trượt cho hệ trên: V SS 0 (3.82) Xây dựng bộ điều khiển: Từ (3.78) và (3.79) nếu đặt: SV: Hoàng Đức Thắng - Phan Thanh Hải 72
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do x11 x12 (3.83) x12 f1 (x) g1 (x,u) và: x21 x22 (3.84) x22 f 2 (x) g 2 (x,u) Ta có: S e  e e  x  x e  x  ( f (x) g (x,u)) 1 1 1 1 1 1 1d 1 11 1 1 1d 1 1 1 (3.85) (e1 1x1d 1 f1 (x)) 1 g1 (x,u) K1h(S1 ) K1h(S1 ) (e1 1x1d 1 f1 (x)) g1 (x,u) 1 Tương tự ta cũng có: K 2 h(S 2 ) (e2 2 x2d 2 f 2 (x)) g 2 (x,u) (3.86) 2 Theo (3.77) ta có: g1 (x,u) 1 H u (3.87) g 2 (x,u) f1 (x) 1 c11 c12 x12 1 g1 H H (3.88) f 2 (x) c21 c22 x22 g 2 Chú ý: g1 g1 (x,u) Từ (3.58) ta có được: K h(S ) e  x 1 1 1 1 1d f (x) u 1 1 1 H (3.89) u K 2 h(S 2 ) e2 2 x2d 2 f (x) 2 2 Thay (3.88) vào (3.89) ta có được bộ điều khiển: SV: Hoàng Đức Thắng - Phan Thanh Hải 73
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do K1h(S1 ) e1 1x1d u1  c11 c12 x12 g1 H 1 (3.90) K h(S ) e  x u2 2 2 2 2 2d c21 c22 x22 g 2 2 III.4.3.4. Tính toán giá trị đặt i cho tay máy hai bậc tự do: Để tính toán giá trị đặt cho tay máy hai bậc tự do chúng ta cần giải bài toán động học ngược, từ đó tính toán giá trị đặt cho các khớp: C12 S12 0 l(C12 C1 ) S C 0 l(S S ) 0 A 0A 1 A 12 12 12 1 (3.91) 2 1 2 0 0 1 0 0 0 0 1 theo cách biến đổi toạ độ ta có được: x0 x2 y y 0 0 2 A2 (3.92) z0 z2 1 1 Tuy nhiên, trong bài toán này có x2, y2, z2=0 vì ta chỉ quan tâm tới chuyển động của tâm bàn kẹp do vậy từ (3.91) và (3.92) ta có: x0 l(C12 C1 ) (3.93) y0 l(S12 S1 ) x 2 y 2 0 0 2 2cos( ) l 2 2 x 2 y 2 2l 2 0 0 cos( ) (3.94) 2l 2 2 x 2 y 2 2l 2  arccos 0 0 2 2 2l Từ (3.93) ta có: 2 2 x0 lC1 lC12 (x0 lC1 ) (lC12 ) (3.95) y lS lS 2 2 0 1 12 (y0 lS1 ) (lS12 ) SV: Hoàng Đức Thắng - Phan Thanh Hải 74
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do 2 2 x0 y0 2l(x0C1 y0 S1 ) 0 2 2 x0 y0 x0 y0 (3.96) C1 S1 2l 2 2 2 2 x0 y0 x0 y0 Đặt x cos( ) 0 2 2 x0 y0 (3.97) y sin( ) 0 2 2 x0 y0 Chọn 1 ta có: x 2 y 2  arccos( 0 0 (3.98) 1 2l Như vậy, nếu yêu cầu của bài toán là điều khiển tâm bàn kẹp đi theo một quỹ đạo đã được định trước và được xác định bởi: x x(t) y y(t) z z(t) 0 thì giá trị đặt cho các khớp phải là: 2 2 x0 y0 1 arccos 2l (3.99) x 2 y 2 2l 2  arccos 0 0 2 2 2l SV: Hoàng Đức Thắng - Phan Thanh Hải 75
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do CHƯƠNG IV: MÔ PHỎNG QUÁ TRÌNH CHUYỂN ĐỘNG CỦA ROBOT DÙNG BỘ ĐIỀU KHIỂN TRƯỢT TRÊN NỀN MATLAB AND SIMULINK: IV.1. Tổng quan về Matlab-Simulink: Matlab là một bộ chương trình phần mềm lớn của lĩnh vực toán số. Tên của bộ chương trình chính là từ viết tắt của từ Matrix Laboratory, thể hiện định hướng chính của chương trình là các phép tính vectơr và ma trận. Phần cốt lõi của chương trình bao gồm một số hàm toán, các chức năng xuất nhập cũng như các khả năng điều khiển chu trình mà nhờ đó ta có thể dựng nên các Scripts. Thêm vào phần cốt lõi, có thể dùng các bộ công cụ Toolbox với phạm vi chức năng chuyên dụng mà người sử dụng cần. Simulink là một Toolbox có vai trò đặc biệt quan trọng: vai trò của một bộ công cụ mạnh phục vụ mô hình hoá và mô phỏng các hệ thống kĩ thuật - Vật lý, trên cơ sở sơ đồ cấu trúc dạng khối. Giao diện đồ họa trên màn hình của Simulink cho phép thể hiện hệ thống dưới dạng sơ đồ tín hiệu với các khối chức năng quen thuộc. Simulink cung cấp cho người dùng một thư viện rất phong phú, có sẵn với số lượng lớn các khối chức năng cho các hệ tuyến tính, phi tuyến và gián đoạn. Hơn thế người sử dụng có thể tạo nên các khối riêng cho mình. Sau khi đã xây dựng mô hình của hệ thống cần nghiên cứu, bằng cách ghép các khối cần thiết, thành sơ đồ cấu trúc của hệ, ta có thể khởi động quá trình mô phỏng. Trong các quá trình mô phỏng ta có thể trích tín hiệu hiện tại vị trí bất kì của sơ đồ cấu trúc và hiển thị đặc tính của tín hiệu đó trên màn hình. Hơn thế nữa, nếu có nhu cầu ta còn có thể cất giữ các đặc tính đó vào môi trường nhớ. Việc nhập hoặc thay đổi tham số của tất cả các khối cũng có thể thực hiện được rất đơn giản bằng cách nhập trực tiếp hay thông qua matlab. Để khảo sát hệ thống, ta có thể sử dụng thêm các Toolbox như Signal Processing (xử lý tín hiệu), Optimization (tối ưu) hay Control System (hệ thống điều khiển). SV: Hoàng Đức Thắng - Phan Thanh Hải 76
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do IV.2. Các thao tác thực hiện mô phỏng: Khớp 1: x x 11 12 1 3 x12 (u1 x22 (2x12 x22 )sin x21 15(cos x11 cos x21 )) 4 9 4cos2 2 x21 3 3 (1 cos x )(u x 2 sin x 15cos(x x )) 2 21 2 2 12 21 11 21 Khớp 2: x x 21 22 1 3 3 x22 ( (1 cos x21 )(u1 x22 (2x12 x22 )sin x21 15(cos x11 cos x21 ))) 4 9 4cos2 2 2 x21 3 (5 3cos x )(u x 2 sin x 15cos(x x ))) 21 2 2 12 21 11 21 Ta có mô hình Simulink của Robot 2 bậc tự do: SV: Hoàng Đức Thắng - Phan Thanh Hải 77
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Khối Subsytem: Ở đây ta sử dụng các khối: + intergrator: khối tích phân với các tham số của khối mặc định cho trước + khối mux: chập tín hiệu đơn thành tín hiệu tổng hợp của nhiều tín hiệu + khối input, ouput: đầu vào và đầu ra của tín hiệu. + khối hàm: biểu diễn 1 hàm toán học khi có tín hiệu đi vào là các biến, tín hiệu ra thu được là hàm cần biểu diễn: . x12 =(u[5]+1.5*u[4]*sin(u[3])*(2*u[2]+u[4])-15*(cos(u[1])-cos(u[3]))- (1+1.5*cos(u[3]))*(u[6]-1.5*u[2]*u[2]*sin(u[3])-15*cos(u[1]+u[3])))/(4- 2.25*cos(u[3])*cos(u[3])) . Và:x22 =(-(1+1.5*cos(u[3]))*(u[5]+1.5*u[4]*sin(u[3])*(2*u[2]+u[4])- 15*(cos(u[1])-cos(u[3])))+(5+3*cos(u[3]))*(u[6]-1.5*u[2]*u[2]*sin(u[3])- 15*cos(u[1]+u[3])))/(4-2.25*cos(u[3])*cos(u[3])) SV: Hoàng Đức Thắng - Phan Thanh Hải 78
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do + Các khối scope (thuộc thư viện con sinks): Hiển thị các tín hiệu của quá trình mô phỏng theo thời gian. Nếu mở cửa sổ Scope sẵn từ trước khi bắt đầu mô phỏng ta có thể theo dõi trực tiếp diễn biến của tín hiệu. Ta sử dụng nguồn tín hiệu u1, u2 là 1(t) Các hàm x’12, và x’22 là : X’12 = (u[5]+1.5*u[4]*sin(u[3])*(2*u[2]+u[4])-15*(cos(u[1])- cos(u[3]))-(1+1.5*cos(u[3]))*(u[6]-1.5*u[2]*u[2]*sin(u[3])- 15*cos(u[1]+u[3])))/(4-2.25*cos(u[3])*cos(u[3])). X’22 = (-(1+1.5*cos(u[3]))*(u[5]+1.5*u[4]*sin(u[3])*(2*u[2]+u[4])- 15*(cos(u[1])-cos(u[3])))+(5+3*cos(u[3]))*(u[6]-1.5*u[2]*u[2]*sin(u[3])- 15*cos(u[1]+u[3])))/(4-2.25*cos(u[3])*cos(u[3])) với u[1], u[2], u[3], u[4], u[5], u[6] tương ứng là các vị trí thứ tự trên khối Mux. u[1] = x11, u[2] = x12, u[3] = x21, u[4] = x22, u[5] = u1, u[6] = u2. Sau khi mô phỏng ta có đồ thị các đường đặc tính của các biến trạng thái x11, x12, x21,x22 là : SV: Hoàng Đức Thắng - Phan Thanh Hải 79
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Mô phỏng dạng hàm điều khiển: Từ các công thức: K1h(S1 ) e1 1x1d u1  c11 c12 x12 g1 H 1 K h(S ) e  x u2 2 2 2 2 2d c21 c22 x22 g 2 2 - Ma trận H: h11 5 3cos 2 h12 h21 1 3/ 2cos 2 h22 1 - Ma trận C:  c11 3 2 sin 2  c12 3/ 2 2 sin 2  c21 31 sin 2 c22 0 SV: Hoàng Đức Thắng - Phan Thanh Hải 80
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do - Ma trận G: g1 15cos1 15cos(1  2 ) g 2 15cos(1  2 ) Ta có: K1h(S1 ) e1 1x1d K2h(S2 ) e2 2 x2d U1=h11 +h12 +c11x12+c12x22+15 1 2 cos(x11) -15sos(x11+x21) K1h(S1 ) e1 1x1d K2h(S2 ) e2 2 x2d U2=h21 +h22 +c21x12+c22x22 1 2 +15cos(x11+x21) Chuyển về dạng hàm của sơ đồ Simulink: U1=u[5]*(5+cos(u[3]))+(1+1.5*cos(u[3]))*u[6]-3*u[2]*u[3]*sin(u[3]) -1.5*u[3]*u[4]*sin(u[3])+15*cos(u[1])-15*cos(u[1]+u[3]) U2=(1+1.5*cos(u[3]))*u[5]+u[6]-3*u[1]*u[2]*sin(u[3])+15*cos(u[1]+u[3]) với u[1]=x11, u[2]=x12, u[3]=x21, u[4]=x22. K h(S ) e  x U[5]= 1 1 1 1 1d 1 K h(S ) e  x U[6]= 2 2 2 2 2d 2 Theo công thức kinh nghiệm ta chọn: K1=K2=500, 1 =2 =0,156. Trường hợp này ta dùng hàm H(S 1),H(S2) là các hàm giới hạn đầu vào trong khoảng giá trị upper và giá trị lower. từ đó ta có: . K1 1 u[5]= H(S1)+ e1 + x1d ( Khối subsytem) 1 1 SV: Hoàng Đức Thắng - Phan Thanh Hải 81
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do . K2 1 u[6]= H(S2)+ e2 + x2d ( Khối subsytem2) 2 2 SV: Hoàng Đức Thắng - Phan Thanh Hải 82
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do SV: Hoàng Đức Thắng - Phan Thanh Hải 83
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Ta chọn khoảng thời gian mô phỏng là từ 0 ->30 s , tức giá trị stop time = 30 ở trong Congiguration Parameters. Khi chay sơ đồ Simulink ta được các kết quả đường đặc tuyến của 2 hàm điều khiển U1, U2, và sai số e1 , e2 là: Ta có các giá trị đặt xd1 và xd2 là: Từ các công thức: x0 l(C12 C1 ) y0 l(S12 S1 ) x cos( ) 0 2 2 x0 y0 y sin( ) 0 2 2 x0 y0 SV: Hoàng Đức Thắng - Phan Thanh Hải 84
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do 2 2 x0 y0 1 arccos 2l x 2 y 2 2l 2  arccos 0 0 2 2 2l =>xd1= acos(sqrt((cos(u[1]+u[3])+cos(u[1]))^2+(sin(u[1]+u[3])+sin(u[1]))^2)/2)+ acos((cos(u[1]+u[3])+cos(u[1]))/sqrt((cos(u[1]+u[3])+cos(u[1]))^2+(sin(u[1])+sin(u[1] +u[3]))^2)) xd2 = acos(((cos(u[1]+u[3])+cos(u[1]))^2+(sin(u[1]+u[3])+sin(u[1]))^2-2)/2) SV: Hoàng Đức Thắng - Phan Thanh Hải 85
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do Kết luận Các vấn đề đã được giải quyết : Các vấn đề còn tồn tại : Tài liệu tham khảo : SV: Hoàng Đức Thắng - Phan Thanh Hải 86
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do 1. Nguyễn Thiện Phúc: Robot công nghiệp. NXB KHKT 2004. 2. Nguyễn Thiện Phúc: Người máy công nghiệp và sản xuất tự động linh hoạt. NXB KHKT 1991. 3. Nguyễn Phùng Quang: Điều khiển Robot công nghiệp - Những vấn đề cần biết. Tạp chí Tự động hoá ngày nay - số tháng 4, 5, 6 / 2006. 4. Nguyễn Thiện Phúc: Robot - Thế giới công nghệ cao của bạn. NXB KHKT 2005. 5. Đào Văn Hiệp: Kỹ thuật Robot. NXB KHKT 2002. 6. Đinh Gia Tường, Tạ Khánh Lâm: Nguyên lý máy. NXB GD 2003. 7. Nguyễn Hồng Thái: Xây dựng thuật toán điều khiển và mô phỏng động Robot nhiều bậc tự do. Thư viện ĐHBK HN 2002. 8. Lê Huy Tùng: Điều khiển trượt thích nghi động cho đối tượng phi tuyến có mô hình không tường minh. Thư viện ĐHBK HN 2003. 9. Nguyễn Doãn Phước, Phan Xuân Minh, Hán Thành Trung: Lý thuyết điều khiển phi tuyến. NXB KHKT 2003. 10. Nguyễn Doãn Phước, Phan Xuân Minh: Hệ phi tuyến. NXB KHKT 2000. 11. Nguyễn Thương Ngô: Lý thuyết điều khiển tự động hiện đại. NXB KHKT 2003. 12. Nguyễn Thương Ngô: Lý thuyết điều khiển thông thường và hiện đại. NXB KHKT 2002. 13. Nguyễn Doãn Phước: Lý thuyết điều khiển tuyến tính. NXB KHKT 2002. SV: Hoàng Đức Thắng - Phan Thanh Hải 87
Đồ án tốt nghiệp Thiết kế bộ Điều khiển trượt cho Robot 2 bậc tự do 14. Nguyễn Phùng Quang: Matlab & Simulink dành cho kỹ sư điều khiển tự động. NXB KHKT 2004. 15. Nguyễn Hoàng Hải: Lập trình Matlab. NXB KHKT 2003. 16. Solomon: Stability of nonlinear control systems. 1965. 17. Applied Asymptotic Methods in Nonlinear Oscillitions. Thư viện ĐHBK HN 1994. 18. Harry: Nonlinear Modulation Theory. 1971. 19. Jakub Mozaryn, Jerzt E.Kurek: Design of the Sliding Mode Control for the Puma 560 Robot. Institute of Automatic Control and Robotics Warsaw university of Technology. 20. Martin Ansbjerg Kjaer: Sliding Mode Control. Sweden February 6th 2004. 21. C.Abdallah, D.Dawson, P.Dorato, and M.Jamshidi: Survey of Robust Control for Rigid Robots . 22. Mark W.Spong: Motion control of Robot Manipulators. The coordinated Science Laboratory, University of Illinois at Urbana- Champaign, 1308 W. Main St, Urbana, III. 61801 USA. 23. . Andre` Jaritz and Mark W. Spong: An Experimental Comparison of Robust control Algorithms on a Direct Drive Manipulator .IEEE Transaction on Control Systems Technology, Vol. 4, No. 6, November 1996. 24. Austin Blaquie`re: Nonlinear system analysis. Academic Press New York and London 1966. SV: Hoàng Đức Thắng - Phan Thanh Hải 88