Tài liệu Xây dựng mô hình simulink cho đối tượng động cơ dc kích từ vĩnh viễn

TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT Tp HCM KHOA CƠ KHÍ CHẾ TẠO MÁY  BÁO CÁO THÍ NGHIỆM LÝ THUYẾT ĐIỀU KHIỂN TỰ ĐỘNG Xây dựng mô hình Simulink cho đối tượng động cơ DC kích từ vĩnh viễn CBGD: Trần Thụy Uyên Phương Sinh viên: Trương Đăng Khoa Lớp 119110B MSSV: 11911009 Tp HCM SVTH: Trương Đăng Khoa – MSSV: 11911009 1. Mục đích thí nghiệm: - Giúp sinh viên làm quen với việc xây dựng mô hình mô phỏng bằng Simulink - Khảo sát đặc tính vòng hở của động cơ DC kích từ vĩnh viễn. 2. Động cơ điện DC Điển hình của loại này là động cơ điện một chiều kích từ độc lập, điều khiển bằng điện áp phần ứng. Sơ đồ nguyên lý của loại động cơ này được thể hiện trên hình 2.9, trong đó dòng kích từ ik được giữ không đổi. - Tín hiệu vào là điện áp u đặt vào phần ứng, [Volt;V] - Tín hiệu ra là vận tốc góc / góc quay  của động cơ, [rad/s; s] *Phương trình mạch điện phần ứng : *Biến đổi Laplace hai vế phương trình, ta được: *Phương trình mômen điện từ của động cơ : M = Km.i Biến đổi Laplace hai vế ta được: M(s) = Km.I(s) *Phương trình cân bằng mômen trên trục động cơ : Biến đổi Laplace hai vế ta được: P a g e 2 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 a. Ngõ ra là vận tốc  Hàm truyền đạt của động cơ DC với tín hiệu ra vận tốc : Phương trình vi phân bậc hai mô tả hệ: [LJ.s2 + (LB+RJ)s + (Km.Ke+RB)].(s) = Km.U(s) ´ + (LB+RJ).❑ ´ + (Km.Ke+RB). = Km.u(t) LJ.❑ ´ + ❑ LB+ RJ Km . Ke+ RB Km .´+ .= . u(t) LJ LJ LJ Đặt x1 =  ´ x2 = x́ 1= ❑ −LB+ RJ −Km . Ke+ RB + Km .´ . . u(t) LJ LJ LJ x́ 1=x 2  x́2 =−LB+ RJ . x 2− Km . Ke+ RB . x 1 + Km . u(t) LJ LJ LJ ¿ x1 ´ = x́ 2 = ❑ { 0 x́1 =  −Km . Ke+ RB x́2 LJ x  = [1 0] x1 + 0.u( t) 2 [][ 1 −LB+ RJ LJ x 1 + 0 u(t) Km x2 LJ ][ ] [ ] [] P a g e 3 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 Từ phương trình vi phân ta có sơ đồ khối mô tả hệ : *Biểu diển trong Simulink sử dụng phương pháp hàm truyền : Với: R = 1.91; Km = 60.2*10^(-3); Ke = 0.0602; L = 2.5; J = 66.5*10^(-6); B = 2.5*10^(-6); Tạo file model tên DCHT.mdl với nội dung sau : Tạo file mdcht.m với nội dung sau: R=1.91; Km=60.2*10^(-3); Ke=0.0602; L=2.5; J=66.5*10^(-6); B=2.5*10^(-6); sim('DCHT',[0,10],[],[]) P a g e 4 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 figure(1); plot(yout1.time,yout1.signals.values,'LineWidth',2); axis([0 10 0 800]); grid on; title('PHUONG PHAP HAM TRUYEN'); Ta vẽ được đồ thị sau: *Biểu diển trong Simulink sử dụng phương pháp sơ đồ khối: Tạo file model DCSDK.mdl với nội dung: Tạo file m mdcsdk.m với nội dung: R=1.91; Km=60.2*10^(-3); Ke=0.0602; L=2.5; P a g e 5 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 J=66.5*10^(-6); B=2.5*10^(-6); sim('DCSDK',[0,10],[],[]) figure(1); plot(yout2.time,yout2.signals.values,'LineWidth',2) axis([0 10 0 800]); grid on; title('PHUONG PHAP SO DO KHOI'); ´ o và ❑0 bằng 0): Ta vẽ được đồ thị (khi ❑ ´ o và ❑0 lần lượt bằng 0 và một giá trị khác 0 (ở đây chọn 24 Khi thay đổi ❑ đơn vị) ta được đồ thị như bên dưới, nhận thấy tín hiệu lúc bắt đầu không phải từ 0. P a g e 6 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 Tín hiệu từ giá trị khoảng 24 như hình phóng to bên dưới P a g e 7 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 ´ o và ❑0 lần lượt bằng 24 và 0, ta được đồ thị gần như không Khi thay đổi ❑ đổi khác với trường hợp đầu tiên *Biểu diển trong Simulink sử dụng phương trình trạng thái: Ta có các ma trận A, B, C, D đã xác định 0 x́1 =  −Km . Ke+ RB x́2 LJ x  = [1 0] x1 + 0.u( t) 2 [][ 1 −LB+ RJ LJ x 1 + 0 u(t) Km x2 LJ ][ ] [ ] [] Biểu diển trong Simulink bằng cách tạo file DCPTTT.mdl P a g e 8 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 Tạo file mdcpttt.m với nội dung : R=1.91; Km=60.2*10^(-3); Ke=0.0602; L=2.5; J=66.5*10^(-6); B=2.5*10^(-6); sim('DCPTTT',[0,10],[],[]) figure(1); plot(yout3.time,yout3.signals.values,'LineWidth',2); axis([0 10 0 800]); grid on; title('PHUONG TRINH TRANG THAI'); Ta nhận được đồ thị đáp ứng P a g e 9 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 Tạo file dcvt.m với nội dung: R=1.91; Km=60.2*10^(-3); Ke=0.0602; L=2.5; J=66.5*10^(-6); B=2.5*10^(-6); sim('DCHT',[0,10],[],[]) sim('DCSDK',[0,10],[],[]) sim('DCPTTT',[0,10],[],[]) figure(1); subplot(3,1,1); plot(yout1.time,yout1.signals.values,'LineWidth',2) axis([0 10 0 800]); grid on; title('PHUONG PHAP HAM TRUYEN'); subplot(3,1,2); plot(yout2.time,yout2.signals.values,'LineWidth',2) axis([0 10 0 800]); grid on; title('PHUONG PHAP SO DO KHOI'); subplot(3,1,3); plot(yout3.time,yout3.signals.values,'LineWidth',2) axis([0 10 0 800]); grid on; title('PHUONG TRINH TRANG THAI'); P a g e 10 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 Ta nhận được 3 đồ thị hoàn toàn giống nhau: P a g e 11 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 Thay đổi tín hiệu vào bằng dạng sóng sin, tần số 1Hz, biên độ 24V, với file m sửa lại như bên dưới ta nhận được đồ thị R=1.91; Km=60.2*10^(-3); Ke=0.0602; L=2.5; J=66.5*10^(-6); B=2.5*10^(-6); sim('DCHT',[0,10],[],[]) sim('DCSDK',[0,10],[],[]) sim('DCPTTT',[0,10],[],[]) figure(1); subplot(3,1,1); plot(yout1.time,yout1.signals.values,'LineWidth',2) axis([0 10 -800 800]); grid on; title('PHUONG PHAP HAM TRUYEN'); subplot(3,1,2); plot(yout2.time,yout2.signals.values,'LineWidth',2) axis([0 10 -800 800]); grid on; title('PHUONG PHAP SO DO KHOI'); subplot(3,1,3); plot(yout3.time,yout3.signals.values,'LineWidth',2) axis([0 10 -800 800]); grid on; title('PHUONG TRINH TRANG THAI'); P a g e 12 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 Nhận xét: khi tín hiệu đầu vào thay đổi thì tín hiệu đầu ra thay đổi, đúng với cả 3 dạng khảo sát là: hàm truyền, sơ đồ khối và phương trình trạng thái. Khi tín hiệu vào là hằng, tín hiệu ra tiến dần đến một giá trị ổn định. Khi tín hiệu vào là hàm sin, tín hiệu ra tiến dần về dạng dao động sin. P a g e 13 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 b. Ngõ ra là góc quay  Hàm truyền đạt khi ngõ ra là góc quay : Phương trình vi phân bậc ba mô tả hệ: (s).[LJs3+(LB+RJ)s2+(RB+Km.Ke)s]=Km.U(s) ´ ( LB+ RJ ) +❑ ´ . ( RB+ Km . Ke )=Km . u(t ) ⃛ LJ + ❑. => ❑. ´ . RB+ Km . Ke = Km .u( t) ⃛ .́ LB+ RJ +❑ ❑+ LJ LJ LJ Đặt x1 =  ´ x2 = x́ 1 = ❑ ´ x3 = x́ 2 = ❑ LB+ RJ RB+ Km . Ke Km ⃛ ❑=x́ −x 2 . + . u(t) 3=−x 3 . LJ LJ LJ x́ 1=x 2 => x́ 2=x 3 LB+ RJ RB+ Km . Ke Km x́ 3=−x 3 . −x 2 . + .u (t) LJ LJ LJ  = x1 0 1 x́1 0 0 => ⌈ x́2 ⌉= −RB+ Km . Ke 0 x́3 LJ [ 0 x1 1 ⌈ x2 ⌉ + −LB+ RJ x3 LJ 0 0 . u(t ) Km LJ ] [] x1 ¿ [ 1 0 0 ] x 2 + 0.u (t) => x3 [] P a g e 14 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 Từ phương trình vi phân ta có sơ đồ khối mô tả hệ : *Biểu diển trong Simulink sử dụng phương pháp hàm truyền : Tạo file DCHTD.mdl Tạo file mdchtg.m với nội dung R=1.91; Km=60.2*10^(-3); Ke=0.0602; L=2.5; J=66.5*10^(-6); B=2.5*10^(-6); sim('DCHTG',[0,10],[],[]) figure(1); plot(yout1.time,yout1.signals.values,'LineWidth',2); axis([0 10 0 4000]); grid on; title('HAM TRUYEN'); Ta được đồ thị: P a g e 15 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 *Biểu diển trong Simulink sử dụng phương pháp sơ đồ khối: Tạo file model DCSDKG.mdl với nội dung: P a g e 16 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 Tạo file m mdcsdkg.m với nội dung R=1.91; Km=60.2*10^(-3); Ke=0.0602; L=2.5; J=66.5*10^(-6); B=2.5*10^(-6); sim('DCSDKG',[0,10],[],[]) figure(1); plot(yout2.time,yout2.signals.values,'LineWidth',2); axis([0 10 0 4000]); grid on; title('SO DO KHOI'); Ta nhận được đồ thị P a g e 17 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 *Biểu diển trong Simulink sử dụng phương trình trạng thái: Tạo file model DCPTTTD.mdl với nội dung: Tạo file m mdcptttg.m với nội dung: R=1.91; Km=60.2*10^(-3); Ke=0.0602; L=2.5; J=66.5*10^(-6); B=2.5*10^(-6); sim('DCPTTTG',[0,10],[],[]) figure(1); plot(yout3.time,yout3.signals.values,'LineWidth',2; axis([0 10 0 4000]); grid on; title('PHUONG TRINH TRANG THAI'); Ta được đồ thị: P a g e 18 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 Tạo file dcgq.m với nội dung: R=1.91; Km=60.2*10^(-3); Ke=0.0602; L=2.5; J=66.5*10^(-6); B=2.5*10^(-6); sim('DCHTG',[0,10],[],[]) sim('DCSDKG',[0,10],[],[]) sim('DCPTTTG',[0,10],[],[]) figure(1); subplot(3,1,1); plot(yout1.time,yout1.signals.values,'LineWidth',2); axis([0 10 0 4000]); grid on; title('PHUONG PHAP HAM TRUYEN'); subplot(3,1,2); plot(yout2.time,yout2.signals.values,'LineWidth',2); axis([0 10 0 4000]); grid on; title('PHUONG PHAP SO DO KHOI'); subplot(3,1,3); plot(yout3.time,yout3.signals.values,'LineWidth',2); axis([0 10 0 4000]); grid on; title('PHUONG TRINH TRANG THAI'); P a g e 19 | 21 SVTH: Trương Đăng Khoa – MSSV: 11911009 Ta nhận thấy 3 đồ thị hoàn toàn giống nhau: Thay đổi tín hiệu vào bằng dạng sóng sin, tần số 1Hz, biên độ 24V, với file m sửa lại như bên dưới ta nhận được đồ thị R=1.91; Km=60.2*10^(-3); Ke=0.0602; L=2.5; J=66.5*10^(-6); B=2.5*10^(-6); sim('DCHTG',[0,10],[],[]) sim('DCSDKG',[0,10],[],[]) sim('DCPTTTG',[0,10],[],[]) figure(1); subplot(3,1,1); plot(yout1.time,yout1.signals.values,'LineWidth',2); axis([0 10 -100 210]); grid on; title('PHUONG PHAP HAM TRUYEN'); P a g e 20 | 21