Mô phỏng gia tốc kế theo công nghệ MEMS bằng SIMULINK Mô phỏng gia tốc kế bằng Matlab và Simulink để mô hình hóa gia tốc kế, thay đổi các thông số để tìm ra các giá trị của các thông số thep thực t
Lời nói đầuNgày nay, với sự phát triển mạnh mẽ của công nghệ điện tử, việc tạo ra các linh kiện mới có kích thước nhỏ với khả năng ứng dụng cao đã và đang mở ra những triển vọng mới . Một trong những công nghệ mới đang tạo ra những thay đổi mang tính cách mạng là Công nghệ vi cơ điện tử (MEMS - Micro Eletro Mechanical Systems).Cảm biến gia tốc là một trong các thiết bị được chế tạo theo công nghệ vi cơ điện tử. Hiện nay về cơ bản có ba loại cảm biến gia tốc , đó là cảm biến gia tốc kiểu tụ, áp điện và áp trở .Việc lựa chọn sử dụng cảm biến nào là tùy thuộc vào yêu cầu và đặc điểm sử dụng.Cảm biến kiểu áp trở có ưu điểm là công nghệ cấu tạo rất đơn giản. Tuy nhiên nhược điểm của chúng là hoạt động phụ thuộc nhiều vào sự thay đổi nhiệt độ và độ nhạy kém.Cảm biến gia tốc kiểu tụ có độ nhạy cao hơn, nhiễu thấp, hao phí năng lượng ít, đặc tính một chiều ổn định, ít chịu tác dụng của nhiệt độ và đáp ứng tuyến tính. Tuy nhiên mạch điện là phức tạp. Do những ưu điểm nổi bật, cảm biến kiểu tụ được ứng dụng rộng rãi trong nhiều lĩnh vục đời sống.Đề tài này đề cập đến việc mô hình hóa và mô phỏng khối gia tốc kế. Qua đó xác định sự thay đổi của các thông số ảnh hưởng lên gia tốc kế, cũng như sự chuyển đổi tương ứng từ gia tốc (giá trị cần đo) sang giá trị khác như điện dung của tụ hay điện áp (giá trị có thể đo được), nhằm đưa ra một sự lựa chọn về cấu hình thiết kế cụ thể phù hợp với các tiêu chuẩn đã định trước. Trong quá trình thực hiện đề tài, chắc chắn không tránh khỏi những thiếu sót, mong muốn nhận được những ý kiến đóng góp của mọi người để có thể sửa chữa nhằm giúp cho đề tài trở nên hoàn thiện hơn. Xin chân thành cám ơn. 1 Phần 1 : Cơ sở lý thuyết1. BỘ PHẬN CẢM BIẾN VI CƠBộ phận cảm biến được nghiên cứu bởi Druck , Ltd , Groby , Leics . Quá trình chế tạo tại Druck vẫn còn trong giai đoạn thử nghiệm nên có nhiều bản thiết kế, cấu tạo với những đặc tính khác nhau được sử dụng. Bộ phận cảm biến là 1 khối vi cơ sử dụng bộ cảm biến điện dung. Khối gia trọng và 2 cực được làm từ silicon. Bản thiết kế cơ có hệ thống giảm sốc được thiết kế theo kiểu “swastika” (hình 1) hay “articulated” (hình 2). Điểm khác biệt của 2 thiết kế này là khối gia trọng có thể là 1 khối rắn hoặc 1 khối rỗng . Hình 1: Bộ phận cảm ứng vi cơ theo thiết kế “swastika” Hình 2 : Khối gia trọng theo thiết kế của “articulate2 Để làm giảm hệ số giảm chấn (damping factor), trong thiết kế của “swastika” khối gia trọng được khoan một hoặc vài lỗ, còn trong thiết kế của “articulated” được thay bởi một khe hình chữ thập giữa khối để chất khí thoát dễ dàng từ mặt này sang mặt khác .Hình dưới đây (hình 3) thể hiện khối gia trọng theo thiết kế “swastika” với 4 lỗ khoan Hình 3Khối gia trọng hoạt động như một cực chung của 2 tụ nối tiếp 1Cvà 2C.Các thông số: • Bộ cảm biến “swastika”: - Khoảng cách từ khối gia trọng đến mỗi cực , d0 : 10 μm (± 1μm) - Diện tích khối gia trọng , A : 2,4mm x 2,4mm = 5,76 mm2 - Độ dày khối gia trọng , h : 0,35 mm- Trọng lượng khối gia trọng , m : 6,2 mg (đặc) ; 2,2 mg (rỗng) (±10%)- Điện dung của mỗi tụ ( giả sử điện trường đều và khối gia trọng ở chính giữa các cực ) : C1 = C2 = = 5,1 pF Với : ε0: hằng số điện môi εr: hằng số điện môi tỉ đối (εr xấp xĩ 1 trong môi trường không khí ) • Bộ cảm biến “ articulated ” :- Khoảng cách từ khối gia trọng đến mỗi cực , d0 : 10 μm (± 1μm)- Diện tích khối gia trọng , A : 12 mm23 - Độ dày khối gia trọng , h : 0,35 mm- Trọng lượng khối gia trọng , m : 8,2 mg(đặc) , 3,8 mg (rỗng) (±10%) - Điện dung của mỗi tụ : C1= C2 = 10,1pF.2. NGUYÊN LÝ HOẠT ĐỘNG : Mô hình tương đương của một gia tốc kế được vẽ ở hình dưới đây (hình 4) Hình 4Tác dụng một lực F theo phương y, ta có, tổng hợp lực tác dụng lên hệ kín tại vị trí cân bằng : với : m : trọng lượng khối gia trọng x : độ dịch chuyển của khối gia trọng b : hệ số giảm chấn k : hệ số đàn hồi ( độ cứng k) Ta bíết rằng, đạo hàm bậc hai theo thời gian của x hay là độ dịch chuyển chính là gia tốc a của khối gia trọng.Tuy nhiên, do cấu tạo của gia tốc kế với khối gia trọng dịch chuyển giữa hai bản mặt cố định, khi khối gia trọng này dịch chuyển, nó làm khối chất khí bên trong bị giãn ra hay bị nén lại, vì thế hệ số giảm chấn ngòai gía trị cho bởi phương trình tổng hợp lực ở trên, còn được bổ xung bằng một thành phần khác phụ thuộc vào sự thay đổi của độ dịch chuyển của khối gia trọng, ta có : b(d) = với μ : độ nhớt của chất khí bên trong.Ngoài ra, dựa vào đặc tính điện của tụ với hai bản mặt song song, trong đó khoảng cách giữa khối gia trọng và hai bản mặt cố định có thể thay đổi . Gọi hằng số điện môi giữa khối gia trọng với hai bản mặt cố định là ε, giả sử 4 khối gia trọng nằm giữa hai bản mặt cố định với khoảng cách là x = d0, diện tích của khối gia trọng là A. Điện dung thu được có dạng : C1 = C2 = C = Nếu khối gia trọng dịch chuyển giữa hai bản mặt cố định một đoạn là d so với vị trí ban đầu, ta có: C1 = và C2 = Chênh lệch điện dung giữa hai tụ là : ΔC = C1 – C2 = = Nếu dịch chuyển d là nhỏ so với d0 , chênh lệch điện dung giữa hai tụ có thể rút gọn : ΔC = Như vậy, ta có thể xác định độ dịch chuyển của khối gia trọng dưới tác dụng của lực F, tức xác định được gia tốc thông qua việc xác định giá trị ΔC. Mặt khác, để chuyển đổi chênh lệch điện dung giữa hai tụ, ta có thể sử dụng mạch pick_off dưới đây, khi đó điện áp ra sẽ là : Phần 2 : Thực hiện1. MÔ HÌNH SỬ DỤNG SIMULINKTừ cở sở lý thuyết đã nêu ở trên và sử dụng thư viện SIMULINK trong Simulink Library Browser, ta có mô hình gia tốc kế như hình 5. 5 Hình 5Hoặc sử dụng thư viện SimMechanics trong Simulink Library Browser, ta có mô hình gia tốc kế như hình 6.Hình 6Việc sử dụng thư viện SimMechanics dễ dàng mở rộng sang mô hình trong không gian ba chiều x, y, z như hình 7 dưới đây. 6 Hình 72. SỬ DỤNG GUIDETrong phần này, ta sẽ sử dụng công cụ giao diện đồ họa với người dùng (GUIDE) của Matlab để thực hiện kích họat các mô hình trên trong SIMULINK, thu nhận các thông tin trả về từ họat động của các mô hình này và xử lý các thông tin đã thu nhận được.Các dòng lệnh sau được sử dụng để kích họat các mô hình: options = simset('SrcWorkspace','current'); if chon_mh==1 %mô hình dùng thư viện SIMULINK sim('mh_gtk',[],options); elseif chon_mh==2 % mô hình dùng thư viện SimMechanics sim('detai_gtk',[],options); else % mô hình 3 chiều x, y, z sim('mh3c_gtk',[],options); endHình dưới đây (hình 8), trình bày giao diện đồ họa GUIDE trong Matlab đã được xây dựng để thực hiện các yêu cầu trên.7 Hình 8Trong phần giao diện đồ họa được thiết kế này, cho phép ta thay đổi các giá trị của các thông số ảnh hưởng lên giá trị của gia tốc, một cách tự động hay bằng tay, bao gồm các thông số như :• Trọng lượng khối gia trọng (m)• Hệ số giảm chấn (b)• Hệ số đàn hồi (k)• Diện tích khối gia trọng (A)• Khỏang cách ban đầu giữa khối gia trọng và hai bản cực cố định (d0) • Biên độ tín hiệu vào• Tần số sóng (sóng sin), ở đây ta sử dụng sóng tuần hòan (sóng sin) và tín hiệu không tuần hòan (xung bậc) để mô phỏng.• Điện áp cấp trên hai bản cực cố định. 8 Hình dưới đây (hình 9), minh họa sự thay đổi của gia tốc theo trọng lượng khối gia trọng (m), biểu diễn sự thay đổi của gia tốc, độ dịch chuyển, điện dung C1 và C2 theo thời gian, cũng như sự liên hệ giữa thay đổi gia tốc và điện dung C1, C2, sử dụng lực tác động là sóng sin. Hình 9Tương tự như trên, nhưng với lực tác động là dạng xung bậc, ta có hình minh họa (hình 10).9 Hình 10Việc thực hiện tương tự cho các thông số còn lại. Chẳng hạn như ở hình 11, minh họa sự thay đổi của gia tốc theo diện tích khối gia trọng (m), biểu diễn sự thay đổi của gia tốc, độ dịch chuyển, điện dung C1 và C2 theo thời gian, cũng như sự liên hệ giữa thay đổi gia tốc và điện dung C1, C2, ứng với lực tác động là sóng sin.10 [...]... giữa gia tốc và điện áp lối ra có dạng đường thẳng tức là liên hệ này là tuyến tính. Ta hãy xem hình 17 với đồ thị biểu diễn sự liên hệ giữa gia tốc và điện áp lối ra theo sự thay đổi của trọng lượng khối gia trọng. 16 Hình dưới đây (hình 9), minh họa sự thay đổi của gia tốc theo trọng lượng khối gia trọng (m), biểu diễn sự thay đổi của gia tốc, độ dịch chuyển, điện dung C 1 và C 2 theo thời gian,... của một gia tốc kế được vẽ ở hình dưới đây (hình 4) Hình 4 Tác dụng một lực F theo phương y, ta có, tổng hợp lực tác dụng lên hệ kín tại vị trí cân bằng : với : m : trọng lượng khối gia trọng x : độ dịch chuyển của khối gia trọng b : hệ số giảm chấn k : hệ số đàn hồi ( độ cứng k) Ta bíết rằng, đạo hàm bậc hai theo thời gian của x hay là độ dịch chuyển chính là gia tốc a của khối gia trọng. Tuy... tính điện của tụ với hai bản mặt song song, trong đó khoảng cách giữa khối gia trọng và hai bản mặt cố định có thể thay đổi . Gọi hằng số điện môi giữa khối gia trọng với hai bản mặt cố định là ε, giả sử 4 Phần 3 : Kết luận Đề tài này, sử dụng công cụ SIMULINK để tạo ra các mơ hình theo u cầu, đồng thời kết hợp công cụ giao diện đồ họa với người sử dụng (GUI) để truy cập các thơng tin có được... Hình 5 Hoặc sử dụng thư viện SimMechanics trong Simulink Library Browser, ta có mơ hình gia tốc kế như hình 6. Hình 6 Việc sử dụng thư viện SimMechanics dễ dàng mở rộng sang mơ hình trong khơng gian ba chiều x, y, z như hình 7 dưới đây. 6 diễn giữa hai đại lượng này, và biểu diễn theo thời gian của gia tốc, điện dung của tụ C 1 và C 2 như hình 16 dưới đây, với sự thay đổi... phục vụ cho việc thiết kế bộ cảm biến gia tốc một cách phù hợp. 12 Hình 17 Trong hình này, ta thấy rằng ứng với trọng lượng khối gia trọng từ 6 mg đến 9 mg, biểu diễn trên đồ thị sự liên hệ giữa gia tốc và điện áp lối ra có dạng đường thẳng, khi trọng lượng khối gia trọng tăng hơn 9 mg, dạng biểu diễn sự liên hệ này khơng cịn tuyến tính. Vì thế việc chọn trọng lượng khối gia trọng vào khỏang 7... lựa chọn các thơng số thích hợp cho bộ cảm biến gia tốc. Mặt khác, với sự lựa chọn này, đề tài còn thực hiện việc thử nghiệm các thơng số đã lựa chọn có phù hợp với yêu cầu cũng như với lý thuyết Tuy nhiên, đề tài cịn một số hạn chế sau: 1. Chưa có sự kiểm chứng thực tế trên một lọai cảm biến gia tốc đã được sản xuất. 2. Việc xây dựng mô hình gia tốc kế trong hệ trục tọa độ x, y, z chỉ là khảo sát... khối gia trọng có thể dịch chuyển và chạm vào hai bản mặt cố định, dẫn đến khơng thể xác định được gia tốc. Vì thế, khỏang cách giữa khối gia trọng và 2 bản cực cố định (d 0 ) phải đủ lớn để không thể xảy ra hiện tượng này Sau khi đã thực hiện việc chọn lựa các thông số cho bộ cảm biến gia tốc, ta sẽ kiểm tra sự lựa chọn này thông qua một giao diện khác với phím nhấn “Thu nghiem” trong giao diện... trọng. Tuy nhiên, do cấu tạo của gia tốc kế với khối gia trọng dịch chuyển giữa hai bản mặt cố định, khi khối gia trọng này dịch chuyển, nó làm khối chất khí bên trong bị giãn ra hay bị nén lại, vì thế hệ số giảm chấn ngịai gía trị cho bởi phương trình tổng hợp lực ở trên, còn được bổ xung bằng một thành phần khác phụ thuộc vào sự thay đổi của độ dịch chuyển của khối gia trọng, ta có : b(d) = với... dựng. 19 Tên thơng số Giá trị theo chủ quan Gía trị theo lý thuyết Ghi chú Trọng lượng khối gia trọng (m) (mg) 7 8,2 Hệ số giảm chấn (b) ( ) 8,2 8,83 Hệ số đàn hồi (k) () 90 90 Diện tích khối gia trọng (A) (mm 2 ) 22 12 Khỏang cách ban đầu giữa khối gia trọng và 2 bản cực cố định (d 0 ) (μm) 120 10 (*) (*) Sự khác biệt là do nếu ta sử dụng khỏang cách giữa khối gia trọng và 2 bản cực cố định... như sự liên hệ giữa thay đổi gia tốc và điện dung C 1 , C 2 , sử dụng lực tác động là sóng sin. Hình 9 Tương tự như trên, nhưng với lực tác động là dạng xung bậc, ta có hình minh họa (hình 10). 9 Hình 12 Quá trình khảo sát sự thay đổi của các thơng số nói trên ảnh hưởng lên gia tốc của khối gia trọng cần đo, được thực hiện trên các mơ hình, với lực tác động lên khối gia trọng ở đây sử dụng hai . việc mô hình hóa và mô phỏng khối gia tốc kế. Qua đó xác định sự thay đổi của các thông số ảnh hưởng lên gia tốc kế, cũng như sự chuyển đổi tương ứng từ gia. hàm bậc hai theo thời gian của x hay là độ dịch chuyển chính là gia tốc a của khối gia trọng.Tuy nhiên, do cấu tạo của gia tốc kế với khối gia trọng dịch