Nghiên cứu nhằm tìm ra những nghiệm gần đúng của bài toán với sai số nghiệm nhỏ so với các nghiệm mô phỏng số trong trường hợp hệ đang xét có tính phi tuyến, thậm chí là phi tuyến mạnh. Việc triển khai nghiên cứu sẽ tạo ra khả năng thu được các kết quả mới chính xác hơn các kết quả đã biết, mở đường cho một cách tiếp cận mới trong nghiên cứu các kết cấu hàng không và vũ trụ.
MỞ ĐẦU 1. Tính cấp thiết của đề tài Khi máy bay chuyển động trong dịng khí sẽ xuất hiện các hiệu ứng khí động học, trong đó dao động flutter của thiết diện cánh máy bay rất được quan tâm. Phân tích đáp ứng của thiết diện cánh máy bay là một bài tốn quan trọng phục vụ q trình thiết kế, chế tạo, vận hành và bảo dưỡng máy bay. Thiết diện cánh chuyển động trong dịng khí khơng nén được thường được mơ hình bởi hệ tự dao động và có tính chất phi tuyến, do vậy phải phát triển các phương pháp đã có để có thể thu được lời giải đạt được độ chính xác mong muốn. Một cách tiếp cận mới, với tên gọi cách tiếp cận đối ngẫu, ban đầu được đề xuất trong nghiên cứu dao động ngẫu nhiên của các hệ phi tuyến với kích động ngồi ồn trắng. Ý tưởng của phương pháp được mở rộng sang nghiên cứu điều khiển giảm dao động cho hệ TMD. Cách tiếp cận đối ngẫu ở trên có tính linh hoạt và có thể áp dụng được cho nhiều lớp hệ phi tuyến khác nhau. Đây cũng là chủ đề của luận án với mục đích nghiên cứu phát triển và áp dụng cho bài tốn ổn định thiết diện cánh máy bay. Nghiên cứu nhằm tìm ra những nghiệm gần đúng của bài tốn với sai số nghiệm nhỏ so với các nghiệm mơ phỏng số trong trường hợp hệ đang xét có tính phi tuyến, thậm chí là phi tuyến mạnh. Việc triển khai nghiên cứu sẽ tạo ra khả năng thu được các kết quả mới chính xác hơn các kết quả đã biết, mở đường cho một cách tiếp cận mới trong nghiên cứu các kết cấu hàng khơng và vũ trụ 2. Mục tiêu nghiên cứu Phát triển phương pháp luận cho cách tiếp cận đối ngẫu trong phương pháp tuyến tính hóa tương đương áp dụng trong bài tốn phân tích đáp ứng của thiết diện cánh chịu lực khí động Xây dựng các cải tiến có hiệu quả cho tiêu chuẩn tuyến tính hóa tương đương đối ngẫu, áp dụng cho bài tốn flutter phi tuyến của thiết diện cánh Tăng độ chính xác cho nghiệm của bài tốn ổn định flutter thiết diện cánh bằng cách áp dụng tiêu chuẩn đối ngẫu được cải tiến Thu được các kỹ thuật tính tốn theo tiêu chuẩn đối ngẫu cho bài tốn flutter của thiết diện cánh. 3. Đối tượng nghiên cứu Thiết diện cánh máy bay theo mơ hình hai chiều chịu tác động của lực khí động. 4. Nội dung nghiên cứu 4.1. Phương pháp nghiên cứu Sử dụng các phương pháp của cơ học để xây dựng mơ hình tính tốn. Áp dụng lý thuyết khí động học xây dựng phương trình dao động flutter của mơ hình thiết diện cánh máy bay Sử dụng các phương pháp giải tích, đặc biệt phát triển phương pháp tuyến tính hóa tương đương của cơ học phi tuyến Sử dụng các phương pháp CFD, phương pháp số mơ phỏng hệ phi tuyến, các số liệu thực nghiệm đã có để so sánh, đánh giá kết quả lý thuyết. 4.2. Hướng giải quyết Trên cơ sở hồn thiện mơ hình cơ học và các kết quả lý thuyết đã có về dao động của thiết diện cánh, đề tài tập trung phát triển cách tiếp cận đối ngẫu để xây dựng kỹ thuật tính tốn mới với cách mở rộng tiêu chuẩn đối ngẫu cho bài tốn phân tích đáp ứng phi tuyến của thiết diện cánh dưới tác động của lực khí động. 4.3. Kết quả dự kiến Xây dựng thành cơng các cải tiến có hiệu quả đối với tiêu chuẩn đối ngẫu cho hệ dao động phi tuyến tuần hồn và ngẫu nhiên Áp dụng cho mơ hình dao động flutter 2 chiều của thiết diện cánh, xác định các hiện tượng mất ổn định flutter và vận tốc gió tới hạn Đánh giá sai số của nghiệm và đảm bảo sai số của nghiệm được cải thiện so với các kết quả đã có trước đây 5. Cấu trúc của luận án Luận án được trình bày theo thứ tự các nội dung chính sau: Mở đầu: Trình bày tính cấp thiết, mục đích và nhiệm vụ nghiên cứu, ý nghĩa khoa học và ý nghĩa thực tiễn của luận án Chương 1: Trình bày các kiến thức cơ sở liên quan đến lĩnh vực khí đàn hồi, sự tương tác giữa ba loại lực: khí động, đàn hồi và qn tính; tổng quan các nghiên cứu quốc tế và trong nước liên quan, các vấn đề cơ bản về mơ hình hóa thiết diện cánh, các hiện tượng phi tuyến và cách tiếp cận đối ngẫu Chương 2: Xây dựng mơ hình thiết diện cánh máy bay chuyển động trong dịng khí khơng nén được. Phân tích đáp ứng và các hiện tượng dao động flutter Chương 3: Trình bày tiêu chuẩn đối ngẫu có trọng số cho vấn đề tuyến tính hóa tương đương hệ dao động phi tuyến. Nghiên cứu đề xuất 3 cách lựa chọn tương ứng với 3 cải tiến cho giá trị tham số trọng số và áp dụng cho các hệ phi tuyến dạng đa thức. Chương 4: Sử dụng các phương pháp số cho phương trình vi phân chuyển động của thiết diện cánh. Kết quả của phương pháp mơ phỏng số và các kết quả của tác giả khác sẽ dùng để đánh giá, so sánh với các kết quả giải tích thu được theo kỹ thuật tính tốn đối ngẫu Kết luận CHƯƠNG 1. TỔNG QUAN VỀ BÀI TỐN PHÂN TÍCH ĐÁP ỨNG CỦA THIẾT DIỆN CÁNH 1.1. Khái niệm cơ bản về khí đàn hồi Khí đàn hồi (Aeroelasticity) là ngành khoa học nghiên cứu các hiện tượng xảy ra do tương tác giữa lực khí động (aerodynamic), lực quán tính (inertia) và lực đàn hồi (elastic). Các hiện tượng khí đàn hồi có ảnh hưởng lớn tới việc thiết kế và hiệu quả hoạt động của máy bay Trong luận án này tập trung vào hiện tượng khí đàn hồi động (là tâm của tam giác khí đàn hồi Collar (Collar, 1978), cho trên Hình 1 trong luận án) 1.2. Các nghiên cứu đáp ứng của thiết diện cánh Phân tích đáp ứng của thiết diện cánh là bài tốn quan trọng phục vụ q trình thiết kế, chế tạo, vận hành và bảo dưỡng máy bay. Để tăng lực nâng, giảm lực cản, cánh cần được thiết kế có dạng mỏng. Điều này lại dẫn tới độ nhạy cảm với dao động tăng lên, đặc biệt khi chuyển động trong dịng khí với vận tốc lớn. Lúc này, lực khí động (lực nâng) tăng rất lớn làm biến dạng hình học của cánh, từ đó lại làm thay đổi đặc trưng của dịng khí, dẫn tới các hiện tượng tương tác khí đàn hồi. Về mặt mơ hình hóa, cánh có thể được mơ tả từ đơn giản đến phức tạp, ví dụ cho trên các Hình 25. Để mở đầu cho những nghiên cứu vấn đề tương tác khí động lực theo cách tiếp cận đối ngẫu, trong luận án này giới hạn xem xét theo mơ hình cho trên Hình 2, trong đó mặt cắt điển hình được đỡ bởi các lị xo phi tuyến và các lực khí động được tính tốn từ lý thuyết cánh mỏng, có thế được quy về lực tập trung tại tâm khí động (Fung 1993). 1.3. Thiết diện cánh phi tuyến Các nghiên cứu về thiết diện cánh máy bay thường tập trung tính tốn đáp ứng nhằm chỉ ra vịng giới hạn trong bài tốn phi tuyến, các hiện tượng rẽ nhánh Hopf và hiện tượng flutter. Phương pháp mà các tác giả sử dụng là phương pháp cân bằng điều hịa và phương pháp tiệm cận (Liu và Zhao (1992)), phương pháp tuyến tính hóa tương đương và phương pháp đa tạp trung tâm (Shahrzad và Mhazoon (2002), Chen vcs 2013). Phương trình chuyển động ứng với mơ hình thường là hệ tự dao động và có tính chất phi tuyến (tính phi tuyến kết cấu và tính phi tuyến khí động lực), có xuất hiện vịng giới hạn, các hiện tượng rẽ nhánh Hopf và hiện tượng flutter. Tính phi tuyến kết cấu có thể phân bố hoặc tập trung tại một vài vị trí. Với mục đích nghiên cứu cách tiếp cận đối ngẫu, luận án cũng sẽ chỉ nghiên cứu tính phi tuyến kết cấu tập trung. Hiệu ứng chính của hiện tượng phi tuyến là các dao động vịng giới hạn (limit cycle oscillations: LCOs), có thể được xem như hiện tượng dao động flutter bị chặn (hiện tượng flutter phi tuyến) và một ví dụ của loại dao động này cho trên Hình 8. Hiện tượng LCO cũng là một nội dung nghiên cứu chính của luận án 1.4. Một số nghiên cứu liên quan trong nước Tại Việt Nam cũng có một số nghiên cứu về ảnh hưởng của lực khí động nói chung cũng như hiện tượng flutter trong nhiều vấn đề kỹ thuật khác nhau như hàng khơng, cầu dây văng nhịp lớn, các kết cấu tấm, vỏ… Lã Hải Dũng, (2005), đã khảo sát ảnh hưởng của một số thơng số và nhận thấy nếu tiếp tục tăng tốc độ bay, dao động của cánh trở thành dao động tự kích gọi là flutter uốn xoắn cánh. Hồng Thị Bích Ngọc vcs, (2009), đã nghiên cứu hiện tượng đàn hồi cánh dưới tác dụng của lực khí động. Phạm Duy Hịa vcs, (2014), đã đặt vấn đề kiểm sốt ổn định flutter của kết cấu cầu hệ treo bằng cách thay đổi hình dáng của mặt cắt cầu và các khe slot Đối với các nghiên cứu về cách tiếp cận đối ngẫu cho phương pháp tuyến tính tương đương có luận án của Nguyễn Ngọc Linh (2015) với tiêu đề “Phân tích dao động ngẫu nhiên phi tuyến bằng phương pháp tuyến tính hóa tương đương”. Đây là luận án đầu tiên phát triển các biểu thức cho tiêu chuẩn đối ngẫu và tiêu chuẩn đối ngẫu có trọng số theo phương pháp trung bình bình phương tối thiểu dựa trên quan điểm đối ngẫu trong bài tốn tuyến tính hóa tương đương cho hệ dao động ngẫu nhiên 1.5. Cách tiếp cận đối ngẫu Một cách tiếp cận mới cho bài tốn phi tuyến về dao động và điều khiển kết cấu được đề xuất bởi Nguyễn Đơng Anh (2010). Cách tiếp cận mới được biết với tên gọi cách tiếp cận đối ngẫu. Áp dụng cách tiếp cận này cho lĩnh vực tuyến tính hóa tương đương dẫn đến kỹ thuật cực tiểu bình phương đối ngẫu Cách tiếp cận đối ngẫu có tính linh hoạt và có thể áp dụng được cho nhiều lớp hệ phi tuyến khác nhau. Do vậy luận án có mục đích nghiên cứu phát triển và áp dụng cách tiếp cận này cho bài tốn về thiết diện cánh máy bay Nghiên cứu nhằm tìm đặc trưng dao động flutter với sai số nghiệm nhỏ so với các nghiệm mơ phỏng số trong trường hợp hệ đang xét có tính phi tuyến, thậm chí là phi tuyến mạnh 1.6. Vấn đề nghiên cứu của luận án Luận án giới hạn nghiên cứu thiết diện cánh ở các điểm sau: Thiết diện cánh được xét là một mặt cắt điển hình, hai chiều Lực khí động tác động vào cánh được quy về lực tập trung theo lý thuyết cánh mỏng tựa dừng Chỉ nghiên cứu hiện tượng đáp ứng phi tuyến điển hình là dao động vịng giới hạn LCO Bài tốn điều khiển bằng cánh nhỏ sẽ được xem xét với luật điều khiển đơn giản là luật điều khiển PID Các vấn đề nghiên cứu mới của luận án bao gồm: Phát triển, hồn thiện các kỹ thuật khác nhau của cách tiếp cận đối ngẫu Áp dụng tính tốn đáp ứng cho một số hệ phi tuyến cơ bản Áp dụng tính tốn hiện tượng flutter cho thiết diện cánh Đánh giá tính chính xác của các lời giải gần đúng thu được theo kỹ thuật đối ngẫu dựa trên sự so sánh với nghiệm của các tác giả khác và kết quả mơ phỏng số CHƯƠNG 2. MƠ HÌNH CƠ HỌC CỦA THIẾT DIỆN CÁNH CHUYỂN ĐỘNG TRONG DỊNG KHÍ 2.1. Lực khí động dừng và tựa dừng Lực khí động do dịng khí tác động vào cánh có thể được phân loại thành các loại dừng (steady), tựa dừng (quasisteady) và khơng dừng (unsteady). Trong luận án này ta chỉ xét lực khí động dừng và tựa dừng 2.1.1. Lực khí động dừng Cánh thực tế có thể được giả thiết lý tưởng hóa như là một cánh mỏng có độ dày cánh bằng 0 (coi như một đoạn thẳng) và 2 đại lượng quan trọng nhất được quan tâm đó là góc tới hay góc xung kích và chiều dài dây cung cánh (Hình 9) Tiết diện cánh 2 chiều với giả thiết tâm của hệ tọa độ được chọn tại biên đầu của thiết diện cánh với trục x nằm dọc theo đường dây cung và trục y vng góc với trục x (Hình 10). Nếu có một dịng khí 2 chiều với tốc độ U và góc tới α đi qua thiết diện cánh thì sự nhiễu loạn xuất hiện trong dịng khí làm cho dịng khí có hướng tiếp tuyến với thiết diện cánh và thiết diện cánh mỏng có thể được thay thế bằng một sự phân bố liên tục của các xốy vận tốc trên bề mặt cánh (Katz vcs, 2001). Nếu gọi cường độ của xốy trên một đơn vị độ dài sải cánh và trên độ dài dx theo phương dây cung là (x)dx thì lực nâng tác động lên phần tử dx và tổng lực nâng được xác định theo cơng cơng thức (1), (2) tương ứng; tổng mơ men do lực nâng tạo ra là âm và được xác định theo cơng thức (3) Một điểm xốy tạo ra một dịng chảy có thế (trừ vị trí của điểm xốy) (Katz vcs 2001) và hình ảnh dịng chảy được biểu diễn trên Hình 11. Vận tốc thực tế là tổng hợp của tất cả các vận tốc do các điểm xốy gây ra và vận tốc trên trục x được xác định theo cơng thức (4) và sự phân bố của xốy được định theo cơng thức (5). Sử dụng phép biến đổi và tích phân (Glauert 1959), điều kiện Kutta để tính lực nâng và mơ men theo các cơng thức (13), (14) để từ đó tính được mơ men tại một điểm có khoảng cách x từ biên đầu cánh từ cơng thức (15) và có thể thấy mơ men tại vị trí ¼ dây cung cánh tính từ biên đầu cánh sẽ bằng 0. Vị trí này được gọi là tâm khí động (aerodynamic center) của cánh. Đây là một kết quả quan trọng của lý thuyết cánh mỏng 2.1.2. Lực khí động tựa dừng Với dịng khơng dừng thì các phương trình cơ bản (1), (4) khơng cịn đúng nữa và bài tốn có thể được giải quyết bằng cơng cụ ánh xạ bảo giác Luận án đã sử dụng giả thiết tựa dừng (Fung 1993) và khi đó các phương trình (1) và (4) vẫn cịn đúng Xét một tấm phẳng có 2 bậc tự do (dịch chuyển thẳng đứng h và góc xoay quanh trục tại vị trí x0 sau biên đầu cánh) như Hình 12 đặt trong dịng khí với vận tốc tại vơ cực bằng U và nằm trên trục x. Giá trị của h là dương nếu hướng xuống dưới và giá trị của là dương nếu mũi máy bay hướng lên. Xét hệ tọa độ như trên hình 12. Như cách làm ở mục trước, tấm phẳng sẽ được thay thế bằng một dải xốy. Khi đó thành phần thẳng đứng của vận tốc chất khí phải thỏa mãn phương trình (16) với điều kiện biên được cho trong biểu thức (17). Sử dụng các biểu thức (2), (6), (8), (17) và (7) để tính lực nâng theo (18) và mơ men đối với trục đi qua vị trí x0 được xác định theo cơng thức (19). Do giả thiết là cánh mỏng vì vậy thực tế thì các giá trị thí nghiệm của lực nâng khơng thực sự trùng khớp với lý thuyết do đó ta cần biểu diễn các các biểu thức (18), (19) dưới dạng (20) và (21) Trên thực tế, mơ hình tựa dừng có thể sử dụng cho các bài tốn có tần số dao động khơng q lớn. Trong các bài tốn liên quan tới tần số dao động cao thì mơ hình khơng dừng tổng qt cần phải được sử dụng (Fung 1993) 2.2. Phương trình chuyển động của thiết diện cánh Giả sử thiết diện cánh mơ tả như Hình 13 được đặt trong dịng khí khơng nén được với tốc độ U. Phương trình chuyển động của thiết diện cánh có thể được viết theo điều kiện là tổng của lực và mơ men qn tính và lực và mơ men đàn hồi phải cân bằng với lực khí động bên ngồi (là lực nâng và mơ men uốn được xác định từ (20), (21)). Các phương trình này được xác định thơng qua biểu thức (23) và dạng khác là biểu thức (24) 2.3. Hiện tượng flutter Một kết cấu có mặt cắt khơng trịn sẽ chịu lực dịng khí thay đổi theo góc mà mặt cắt tạo với dịng chảy. Khi kết cấu dao động, góc này cũng dao động lực khí động dao động Nếu lực khí động có xu hướng khuếch đại dao động thì kết cấu trở nên mất ổn định khí động và biên độ dao động có thể trở lên rất lớn. Hiện tượng đó gọi là galloping hoặc flutter. Khái niệm galloping chỉ hiện tượng mất ổn định khí động xảy ra với các mặt cắt tù (bluff), tức là các mặt cắt có tách dịng phía sau. Ngược lại, khái niệm flutter dành cho các mặt cắt khí động, khơng có tách dịng phía sau 2.3.1. Hiện tượng mất ổn định 1 bậc tự do Hiện tượng mất ổn định 1 bậc tự do có thể được giải thích về mặt ngun lý như trên Hình 14. Nếu mặt cắt chỉ có dao động lên xuống, khơng có dao động xoắn ( =0) thì phương trình dao động (24) được viết thành phương trình (25). Theo biểu thức này nếu góc tới nhỏ thì hệ (25) có cản dương và ln ổn định. Tuy nhiên khi góc tới lớn hơn một góc giới hạn thì cánh máy bay xảy ra hiện tượng “chết đứng”. Khi đó hệ (25) có thể có cản dương và xảy ra mất ổn định. Nếu mặt cắt chỉ có dao động xoắn, khơng có dao động lên xuống (h=0), phương trình dao động (24) trở thành phương trình (26) và nếu điều kiện (27) thỏa mãn thì hệ có độ cứng âm và mất ổn định. Đây cịn gọi là mất ổn định tĩnh. 2.3.2. Hiện tượng mất ổn định 2 bậc tự do Hiện tượng mất ổn định 2 bậc tự do có thể có thể được mơ tả như Hình 16. Bỏ qua các đại lượng liên quan tới vận tốc trong phương trình (24), ta được một hệ tuyến tính như biểu thức (28). Xét nghiệm của hệ có dạng như (29) sẽ dẫn tới việc tìm nghiệm của bài tốn giá trị riêng có dạng như (31). Với điều kiện biên giữa miền ổn định cân bằng và miền mất ổn định được biểu diễn bởi biểu thức (33) thì phương trình (31) cho ta biểu thức giải tích của vận tốc tới hạn xảy ra hiện tượng mất ổn định flutter trong hệ tựa tĩnh có dạng (35). Đặc biệt, trong tài liệu (Pines 1958, Dowell vcs 2015), người ta đã chứng minh được rằng nếu S≤0, tức là trọng tâm cánh nằm ở trước tâm đàn hồi thì các nghiệm của (35) đều âm và khơng xảy ra hiện tượng flutter 2.4. Tính tốn vận tốc flutter trong hệ tuyến tính 2.4.1. Hệ tự dao động tổng qt Xét một hệ tự dao động tổng qt được mơ tả bởi phương trình trạng thái Khi vận tốc dịng khí nhỏ hơn một giá trị tới hạn thì đáp ứng của hội tụ về 0 với mọi điều kiện đầu. Lúc đó các giá trị riêng của hệ tuyến tính đều có phần thực âm. Khi vận tốc đạt tới giá trị tới hạn, xuất hiện giá trị riêng có phần thực bằng 0, hệ có thể xuất hiện dao động với biên độ hữu hạn hoặc đáp ứng tiến tới vơ cùng. Vận tốc tới hạn flutter có thể được xác định thơng qua biểu thức (39). Nghiệm thu được chỉ có thể tồn tại trong thực tế nếu nó ổn định với các nhiễu. Điều kiện (40) sẽ xác định một dao động LCO ổn định. 2.4.2. Thiết diện cánh 2 chiều có điều khiển PID Trong chương 4 ta sẽ tập trung xem xét phương trình chuyển động của thiết diện cánh 2 chiều có cánh nhỏ được điều khiển bằng thuật tốn điều khiển PID. Phương trình của hệ này có thể được viết chung dưới dạng phương trình trạng thái 5 chiều theo biểu thức (41). Khai triển đa thức đặc trưng đối với phương trình (38) và sau đó tách các phần thực, phần ảo, sử dụng biến đổi để khử thành phần 4 ta được phương trình (44) để xác định tần số dao động biên độ hữu hạn (nếu nó tồn tại) và phương trình (45) để xác định vận tốc tới hạn 2.5. Tính tốn thiết diện cánh bằng phương pháp CFD 2.5.1. Mơ phỏng khí động lực trên mơ hình cánh máy bay Phân tích CFD trên các mơ hình cánh máy bay sử dụng thiết diện NACA 2412 (Ira 1951) với các thơng số chiều dài dây cung c=0.3m, chiều dài cánh l=1.6m được thực hiện bằng phần mềm ANSYS Fluent và Structural. Một sơ đồ mơ hình hình học của cánh máy bay và thiết diện cánh được thể hiện trong Hình 18. Hợp kim nhơm 7075 T6 (Ma vcs 2014) với các đặc điểm được cho trong Bảng 1 cùng với các đặc tính dịng chảy tương tự như được sử dụng trong các thí nghiệm được áp dụng để thiết lập cho mơ phỏng Trong nghiên cứu này, vận tốc dịng chảy lối vào được thay đổi trong mơ phỏng. Vận tốc được thay đổi trong khoảng 0 đến 50 m/s với từng nấc 5 m/s, nó phù hợp với phạm vi thử nghiệm của các UAV bay trong điều kiện tốc độ thấp. Những mơ phỏng này được lặp đi lặp lại các góc tới là 0 đến 20 độ. Sau đó, lực khí động được đo trong từng mơ phỏng, để xác định hệ số của lực nâng và lực cản, và được so sánh với kết quả lý thuyết Miền phân bố áp suất trong luồng khơng khí, khi vận tốc lối vào được áp dụng 25m/s được trình bày trong Hình 20. Các vùng áp suất cao xuất hiện ở mép trước và mặt dưới của cánh. Bên cạnh đó, các khu vực áp suất thấp hơn xảy ra ở bề mặt trên của cánh. Theo Hình 21 thì ở mép trước và sát bề mặt của cánh máy bay, vận tốc của dịng chảy là gần như bằng khơng. Tuy nhiên, các dịng chảy tăng tốc thay đổi rõ ràng ở bề mặt trên của cánh. Ứng suất tương đương trên cánh được trình bày trên Hình 22. Quan sát ta có thể thấy ứng suất đạt tối đa trên cánh ở phần gắn cố định vào thân máy bay. Các kết quả được so sánh với cơng thức lý thuyết bởi phương trình (46) và (47), với góc tới vơ cùng nhỏ được chọn ( ≈0.025). Kết quả so sánh cho thấy sự phù hợp tốt giữa lý thuyết và mơ phỏng. 2.5.2. Tối ưu hình dạng khí động sử dụng phương pháp SQP SQP là một thuật tốn tối ưu hóa hiệu quả nhất do thực tế rằng nó u cầu số lượng thấp nhất về các tính tốn hàm chức năng. Trong trường hợp bài tốn tối ưu hình dạng cánh máy bay đang xét, thuật tốn chỉ địi hỏi có 26 tính tốn chức năng. Cánh máy bay tối ưu có được đã chứng minh có các 10 trong những phương pháp giải tích được sử dụng phổ biến là phương pháp tuyến tính hóa tương đương, do N. Krylov and N. Bogoliubov đề xuất năm 1937. Ý tưởng cơ bản của phương pháp tuyến tính hóa là thay thế hệ phi tuyến ban đầu bằng một hệ tuyến tính mà một số các tính chất động lực học của hệ phi tuyến có thể được nghiên cứu thơng qua hệ tuyến tính này Từ khi ra đời cho đến nay, phương pháp tuyến tính hóa tương đương đã trải qua nhiều phát triển như tiêu chuẩn tương đương dựa trên điều kiện bình phương tối thiểu (Caughey 1963, Roberts vcs 1990, Socha 2008, Proppe vcs 2003, Langley 1988), được sử dụng phổ biến nên còn được gọi là ‘tiêu chuẩn kinh điển’, ‘tiêu chuẩn thơng thường’. Bên cạnh đó, có nhiều phát triển tiêu chuẩn tương đương khác tiêu chuẩn tương đương năng lượng (Roberts vcs 1990), tiêu chuẩn tuyến tính điều chỉnh (Elishakoff 2009) 3.1.1. Tiêu chuẩn tương đương kinh điển Để trình bày ý tưởng cơ bản và một số phát triển của phương pháp tuyến tính hóa tương đương, ta xét dao động phi tuyến được mơ tả bởi phương trình (62), và biến đổi về dạng phương trình tuyến tính hóa tương đương (63). Theo tiêu chuẩn tương đương kinh điển thì sai số phương trình giữa (62) và (63) phải thỏa mãn điều kiện cực tiểu của trung bình bình phương sai số, theo các phương trình (64), (65) và (66). Giải hệ thu được các hệ số tuyến tính hóa trong các cơng thức (67), (68) phụ thuộc vào các đáp ứng chưa biết của hệ tuyến tính tương đương (63). Ta thấy rằng các tính chất phi tuyến được thể hiện trong việc xác định các hệ số tuyến tính hóa tương đương. Do vậy, có thể chờ đợi việc hệ tuyến tính hóa tương đương phản ánh một số tính chất phi tuyến cơ bản của hệ phi tuyến gốc 3.1.2. Tiêu chuẩn sai số thế năng Elishakoff va Zhang X đa đê xt ph ̀ ̃ ̀ ́ ương phap tun tinh hóa t ́ ́ ́ ương đương dựa trên cực tiêu hoa ham sai sô thê năng cho tr ̉ ́ ̀ ́ ́ ường hợp phương trình (62) là phương trình phi tuyến theo dịch chuyển. Hệ số tuyến tính hóa trong phương trình (63) được xac đinh b ́ ̣ ởi điều kiện cực tiểu sai số giữa thế năng U(x) (70) của phần tử đàn hồi phi tuyến và thế năng của phần tử đàn 12 hồi tuyến tính hóa tương đương dưới dạng trung bình bình phương (69), dẫn đến kết quả trong (71). Tiêu chuẩn (69) áp dụng cho các hệ đàn hồi phi tuyến thu được kết quả khá tốt, tuy nhiên khơng áp dụng được cho các hệ cản phi tuyến 3.1.3. Tiêu chuẩn tương đương điều chỉnh Các tác giả N D Anh và Di Paola đề xuất tiêu chuẩn tuyến tính hóa tương đương chuẩn Gauss có điều chỉnh dựa trên đánh giá: Sự thay thế hệ phi tuyến gốc bằng một hệ tuyến tính tự nó đã làm giảm mức độ phi tuyến (Elishakoff vcs 2009). Do đó, để đạt được nghiệm chính xác hơn thì cần cải thiện sự giảm thiểu tính phi tuyến nêu trên. Các tác giả đã thực hiện việc thay các số hạng phi tuyến bằng các số hạng phi tuyến bậc cao hơn trước khi tuyến tính hóa với qui trình điều chỉnh một bước, theo các phương trình (72), (73) và kết quả trong (74). Áp dụng tiêu chuẩn (73) cho hệ Duffing và Van der Pol đã thu được sự cải thiện đáng kể về độ chính xác của nghiệm. Elishakoff vcs (2009) phát triển tiêu chuẩn của Anh và Di Paola từ điều chỉnh một bước thành điều chỉnh hai bước (75), với các hệ số tuyến tính hóa được xác định theo tiêu chuẩn cực tiểu hóa thơng thường. Sau đó, các tác giả áp dụng qui trình (75) cho dao động Lutes and Sarkani thu được kết quả tốt hơn so với qui trình điều chỉnh một bước ở trường hợp mức độ phi tuyến thay đổi khá lớn 3.2. Tiêu chuẩn đối ngẫu có trọng số Nguyễn Đơng Anh (Anh 2010) đề xuất cách tiếp cận mới được biết với tên gọi cách tiếp cận đối ngẫu với quan điểm tạo ra một sự hài hịa trong nghiên cứu, cho phép phát hiện bản chất của vấn đề một cách đầy đủ hơn Tiêu chuẩn đối ngẫu đầu tiên là tiêu chuẩn tuyến tính hóa có điều chỉnh do Anh và Di Paola đề xuất hay cịn được gọi là điều chỉnh một bước trong đó các tác giả đề nghị làm tăng mức độ phi tuyến của hệ ban đầu trước khi tuyến tính hóa. Elishakoff vcs [2009] phát triển tiêu chuẩn của Anh và Di Paola từ điều chỉnh một bước thành điều chỉnh hai bước và áp dụng cho dao động Lutes Sarkani có mức độ phi tuyến thay đổi mạnh. Tiêu chuẩn đối ngẫu thứ hai là tiêu chuẩn sai số bình phương trung bình khu vựctồn 13 cục được Anh vcs (2012c) cải tiến từ tiêu chuẩn sai số bình phương trung bình khu vực trong đó các tác giả phát triển từ miền đáp ứng tập trung sang miền đáp ứng tồn cục. N. D. Anh (2010) nghiên cứu sự thay thế đối ngẫu từ hệ tuyến tính hóa trở về hệ phi tuyến ban đầu đối với tiêu chuẩn kinh điển và áp dụng cho dao động ngẫu nhiên Duffing, sau đó tác giả đã đề xuất một tiêu chuẩn đối ngẫu có trọng số tổng qt Tiêu chuẩn đối ngẫu có trọng số được trình bày dưới dạng phương trình (76), trong đó p là trọng số được chọn trong khoảng p 1/ Các hệ số tuyến tính hóa và các hệ số trở về sẽ được xác định theo các phương trình (77), (78), (79) và (80) 3.3. Những cải tiến của phương pháp đối ngẫu có trọng số 3.3.1. Cải tiến 1 Ta thấy rằng hệ số tuyến tính hóa tương đương trong (79) phụ thuộc vào giá trị của trọng số p nằm trong khoảng p 1/ Việc xác định giá trị nào của p là bài tốn rất khó và là mục tiêu lâu dài của các nghiên cứu về tiêu chuẩn đối ngẫu. Hiện nay, hai giá trị giới hạn của p đang được sử dụng và quan tâm nhiều, cụ thể p =0 (81) ứng với tiêu chuẩn kinh điển và p =1/2 (82) đã được đề nghị và đánh giá hiệu quả cho các hệ phi tuyến ngẫu nhiên (Anh vcs 2012a, 2012b). Tuy nhiên các áp dụng sơ bộ giá trị p=1/2 cho các hệ phi tuyến dao động tuần hoàn cho các kết quả khơng tốt, có thể là đối với các hệ dao động tiền định cần có sự căn chỉnh giữa các sự thay thế lượt đi và về. Qua nghiên cứu nhận thấy rằng khi hệ có tính phi tuyến yếu thì p xấp xỉ 0, tức là tiêu chuẩn tương đương kinh điển cho sai số nhỏ, và khi tính phi tuyến tăng thì p tăng theo. Với giá trị trọng số p =1/2 thì vai trị của sự thay thế lượt đi và lượt về là như nhau. Tuy nhiên, ta cũng chú ý rằng việc thay thế kinh điển là sự thay thế thứ nhất trong khi sự thay thế đối ngẫu chỉ là thay thế thứ hai. Đối với hệ phi tuyến vừa hai sự thay thế này khơng nên cân bằng với nhau mà sự thay thế thứ nhất nên có vai trị quan trọng hơn. Do đó luận án đề xuất kết hợp 2 giá trị trọng số biên bằng phép trung bình cộng. Với dạng cải tiến thứ nhất, hệ số tương đương được đề nghị sẽ là giá trị trung bình cộng của hai hệ số tương đương biên, 14 tính theo (83) 3.3.2. Cải tiến 2 Trong dạng cải tiến này, ta sử dụng trung bình tích phân chứ khơng phải trung bình đại số như dạng cải tiến 1. Dạng cải tiến này cũng có thể được lý luận theo các thuật ngữ địa phươngtồn cục được sử dụng trong (Anh vcs 2013, 2014). Theo cách lý luận này, ta có thể xem xét rằng các giá trị của trọng số thay đổi trên miền tồn cục của việc lấy tích phân. Và hệ số tuyến tính hóa tương đương được đề xuất là giá trị trung bình của tất cả các hệ số tuyến tính hóa địa phương theo nghĩa tích phân. Hệ số được xác định theo các biểu thức (84) và (85) 3.3.3. Cải tiến 3 Ta chú ý rằng trong tiêu chuẩn đối ngẫu (76) chứa 2 thành phần thể hiện sai số giữa hàm phi tuyến và hàm tuyến tính trong q trình lượt đi và lượt về. Tuy nhiên sau 2 lượt đi và về thì hàm phi tuyến sẽ bị thay đổi theo một hệ số trở về và sẽ có sai số do hệ số trở về khác 1. Như vậy cách sử dụng 2 thành phần sai số chưa thực sự đầy đủ. Trong dạng cải tiến này, 3 dạng sai số sẽ được xem xét với trọng số bằng nhau (86), và (87). Ta thu được biểu thức của hệ số tuyến tính hóa tương đương dạng cải tiến 3 trong phương trình (88) 3.4. Áp dụng cho dao động tự do của hệ phi tuyến dạng Duffing bậc cao Ta xét dao động tự do trong hệ phi tuyến dạng Duffing có dạng (89), với các điều kiện đầu (90) và (91). Đây là hệ bảo tồn do đó ta có tần số dao động chính xác được tính theo cơng thức (92). Phương trình tuyến tính tương đương với (89) sẽ có dạng (93). Áp dụng các cơng thức (79), (83), (85), (88) ta sẽ có các biểu thức tính tần số riêng trong các trường hợp tiêu chuẩn kinh điển p=0 (94), tiêu chuẩn đối ngẫu trọng số với p=1/2 (95), tiêu chuẩn đối ngẫu cải tiến dạng 1 (96), tiêu chuẩn đối ngẫu cải tiến dạng 2 (97), và tiêu chuẩn đối ngẫu cải tiến dạng 3 (98). So sánh các kết quả tính tốn tần số với tần số chính xác được thể hiện trên Bảng 4 Kết quả trong bảng cho thấy các tiêu chuẩn đối ngẫu cải tiến đều cho kết 15 chính xác hơn tiêu chuẩn kinh điển. Ngồi ra các tiêu chuẩn đối ngẫu cải tiến cho kết quả tốt hơn tiêu chuẩn đối ngẫu 2 thành phần với các giá trị của n khơng q lớn (tính phi tuyến khơng q lớn). Điều đó cho thấy tiêu chuẩn đối ngẫu cải tiến có thể áp dụng tốt cho các hệ có tính phi tuyến trung bình 3.5. Áp dụng cho dao động ngẫu nhiên Chuyển sang dao động ngẫu nhiên ta xét hệ Duffing bậc ba có dạng (99), với ồn trắng (100). Phương trình phi tuyến (99) có dịch chuyển trung bình bình phương là (101). Phương trình tuyến tính hóa tương đương ứng với (99) là (102). Theo lý thuyết dao động ngẫu nhiên tuyến tính, ta xác định hệ số tuyến tính hóa theo cơng thức (103), (104) và (105). Áp dụng phương pháp tuyến tính hóa tương đương trong các trường hợp tiêu chuẩn kinh điển p=0 (106), tiêu chuẩn đối ngẫu trọng số với p=1/2 (107), tiêu chuẩn đối ngẫu cải tiến dạng 1 (108), tiêu chuẩn đối ngẫu cải tiến dạng 2 (109), và tiêu chuẩn đối ngẫu cải tiến dạng 3 (110). So sánh các dịch chuyển bình phương trung bình gần đúng và chính xác được thể hiện trên Bảng 5 Các kết quả cho thấy xu hướng giống như trong ví dụ ở mục 3.4. Các dạng cải tiến của tiêu chuẩn đối ngẫu đều cho kết quả tốt hơn tiêu chuẩn kinh điển. Ngồi ra, kết quả sẽ tốt vùng có tính phi tuyến khơng q nhỏ và q lớn CHƯƠNG 4. ÁP DỤNG KỸ THUẬT TUYẾN TÍNH HĨA ĐỐI NGẪU CHO BÀI TỐN PHÂN TÍCH ĐÁP ỨNG PHI TUYẾN CỦA THIẾT DIỆN CÁNH 4.1. Mơ hình thiết diện cánh Trong chương này ta sẽ khảo sát một mơ hình khí động phi tuyến của một thiết diện cánh dựa trên mơ hình thí nghiệm được trình bày trong Strganac vcs 2000. Mơ hình thí nghiệm này được phát triển để cung cấp các số liệu đo trực tiếp của đáp ứng khí động phi tuyến cũng như sự tùy biến trong việc thay đổi các điều kiện và tham số thí nghiệm. Độ lệch tâm của tâm khí động, khối lượng của các thành phần khác nhau của hệ thống, mơ men 16 tĩnh, mơ men qn tính của thiết diện cánh, các đặc tính độ cứng và hình học của thiết diện cánh đều có thể dễ dàng thay đổi trong các nghiên cứu khảo sát. Ngồi ra mơ hình cịn có phần đi cánh có thể điều khiển để kiểm tra các thuật tốn điều khiển khác nhau (hình 36). Thiết diện cánh được gắn trên bệ đỡ linh động cho phép chuyển động với 2 bậc tự do Thiết diện cánh cịn bao gồm phần đi trải dài trên tồn bộ sải cánh (Hình 37, 38) Chuyển động tịnh tiến (plunge) được cung cấp bởi giá đỡ tịnh tiến. Chuyển động xoắn (pitch) độc lập với chuyển động tịnh tiến được cung cấp bánh cam gắn với giá đỡ. Mơ hình hệ thống đỡ cho phép tùy biến trong các điều kiện thí nghiệm và các tham số. Độ cứng kết cấu được tạo ra bởi cặp bánh cam được thiết kế cho phép thay đổi độ cứng tuyến tính và phi tuyến Dạng của bánh cam, độ cứng và độ kéo dãn ban đầu của lị xo tạo ra tính phi tuyến. Đáp ứng của hệ thống được đo bằng các cảm biến gia tốc và các đầu đo quang. Vận tốc dịng chảy được xác định từ ống thổi gió Phương trình chuyển động của mơ hình trên hình 36 đã được viết trong rất nhiều tài liệu tham khảo (Strganac vcs 2000, Li vcs 2011, Platanitis vcs 2004, Ko vcs 2002, A. Abdelkefi vcs 2012) và có dạng (111), tương tự với (24) nhưng được trình bày phù hợp hơn với mơ hình nghiên cứu. Trong mơ hình được trình bày, độ cứng k được xét là một hàm của k dưới dạng đa thức bậc 4 (112) 4.2. Phương trình xác định vận tốc tới hạn Trong luận án này, do tập trung chính vào cách tiếp cận đối ngẫu nên sẽ chỉ quan tâm tới thuật tốn điều khiển đơn giản là điều khiển PID. Đây là thuật tốn điều khiển đơn giản nhất và được sử dụng nhiều nhất. Góc điều khiển sẽ được xác định là tổng của thành phần tỷ lệ (P: proportional), thành phần tích phân (I: Integral) và thành phần vi phân (D: derivative) nhân với một vài hệ số khuếch đại nào đó. Như vậy bộ điều khiển PID được viết dưới dạng (113). Thay (113) vào phương trình (111) ta có thể đưa phương trình về dạng phương trình trạng thái (114), trong đó các hệ số trạng thái được xác định theo biểu thức (115) 17 Sau khi xác định được các hệ số thì bài tốn flutter có thể giải được bằng các phương trình (44), (45). Sử dụng các biểu thức (115) vào (45) ta sẽ thu được phương trình một biến đối với U được viết lại trong (116). Phương trình này có thể dễ dàng giải bằng MATLAB nhờ lệnh fzero. Sau khi tính được vận tốc tới hạn thì tần số dao động LCO được xác định từ (44) và được viết lại trong (117) 4.3. Áp dụng kỹ thuật tuyến tính hóa đỗi ngẫu Các kết quả của các kỹ thuật đối ngẫu được trình bày trong chương 3 được tóm tắt trong (118), (119), (120), (121), (122) và (123), trước khi áp dụng vào bài tốn thiết diện cánh cụ thể. Trong trường hợp độ cứng phi tuyến như trong cơng thức (112), độ cứng phi tuyến sẽ được tuyến tính hóa dưới dạng (124) Đối với độ cứng phi tuyến như trong cơng thức (112) thì hàm phi tuyến và hàm tuyến tính là: Fnl = ( kα 2α + kα 3α + kα 4α + kα 5α ) α ; Fl = α Giả sử đáp ứng điều hịa của góc xoắn có dạng: α = A sin ωt Sử dụng phép trung bình trong trường hợp dao động điều hịa: = ω 2π 2π / ω dt ta sẽ tính được các biểu thức: r2 = Fnl Fl Fnl2 Fl = kα2 �3 A2 kα A4 kα � + � � � � 63 A8 35 A6 A4 A2 + ( 2kα kα +kα2 ) + ( 2kα kα +kα2 ) +kα 128 64 Fnl Fl Fl = A2 kα A4 kα + Vậy hệ số tuyến tính hóa tương đương (124) trong các trường hợp sẽ có dạng tương ứng là (125), (126), (127), (128) và (129). Thay các hệ số tuyến tính hóa tương đương (125) hoặc (126) hoặc (127) hoặc (128) hoặc (129) 18 vào vị trí của ka trong (115), rồi thay vào (116) ta sẽ thấy vận tốc U là hàm của biên độ dao động xoắn A. Từ đó ta sẽ vẽ được đường đặc trưng biên độ LCO với vận tốc tới hạn. Biết được quan hệ giữa U và A thì biểu thức (117) cũng cho ta mối quan hệ giữa tần số LCO với vận tốc tới hạn 4.4. Các ví dụ và tính tốn bằng phương pháp vi phân Trong mục này ta sẽ tóm tắt các số liệu của các trường hợp ví dụ sẽ được nghiên cứu cũng như các kết quả tính tốn số cho các trường hợp 4.4.1. Số liệu đầu vào Các số liệu của các ví dụ lấy theo tài liệu Li vcs 2011. Ta sẽ khảo sát 5 ví dụ sau đây Ví dụ 1: Thiết diện cánh phi tuyến bậc 5, khơng có điều khiển Các số liệu của các tham số trình bày trong mục 4.1 được cho trên bảng 6 Ví dụ 2: Thiết diện cánh phi tuyến bậc 3, khơng có điều khiển Các số liệu được cho trên bảng 7 Ví dụ 3: Thiết diện cánh phi tuyến bậc 5, điều khiển tỷ lệ (P control) Các số liệu được cho trên bảng 8 Ví dụ 4: Thiết diện cánh phi tuyến bậc 5, điều khiển tích phân (I control) Các số liệu được cho trên bảng 9 Ví dụ 5: Thiết diện cánh phi tuyến bậc 5, điều khiển vi phân (D control) Các số liệu được cho trên bảng 10 4.4.2. Tìm vận tốc tới hạn bằng phương pháp số Lời giải số thu được bằng cách giải hệ phương trình vi phân phi tuyến (111). Hệ này được giải bằng hàm ode45 trong MATLAB. Để sử dụng hàm ode45, phương trình vi phân được chuyển về hệ phương trình vi phân cấp 1 (114), trong đó các tham số được tính từ (115), độ cứng phi tuyến được xác định từ (112). Quy trình tính tốn số sau được sử dụng để xác định vận tốc tới hạn và tần số dao động LCO: Đầu tiên, một giá trị ban đầu của góc xốy được cố định. Các giá trị ban đầu khác được gán bằng 0 Thay đổi dần dần vận tốc dịng chảy U. Với mỗi vận tốc dịng chảy, phương trình vi phân được giải. Trong các ví dụ thì khoảng thời gian giải 19 phương trình vi phân từ 0 đến 120s, là đủ dài để thu được dao động ổn định (trong tài liệu Li vcs 2011 chỉ tính đến 30s). Nếu biên độ dao động hội tụ thì vận tốc dịng chảy nhỏ hơn vận tốc tới hạn. Khi đáp ứng bắt đầu hình thành dao động LCO (biên độ dao động tại các thời điểm cuối bằng hoặc lớn hơn biên độ dao động tại thời điểm ban đầu) thì vận tốc đạt tới vận tốc flutter. Khi đó vận tốc dịng, biên độ dao động và tần số dao động sẽ được ghi lại Tăng giá trị ban đầu của góc xốy và lặp lại q trình bên trên. Q trình trên thực tế tính tốn cho thấy khi tăng giá trị ban đầu đủ lớn thì biên độ dao động ổn định hầu như khơng tăng nữa. Đây có thể xác định là vị trí mà vận tốc tới hạn có giá trị nhỏ nhất Trên thực tế quy trình tính tốn như trên mất nhiều thời gian ở việc “tăng dần” vận tốc dịng. Vì vậy để tăng tốc độ tính tốn, quy trình cần làm từ bước thơ đến tinh. Cụ thể, đầu tiên bước vận tốc được lấy giá trị lớn để xác định vận tốc gần với vận tốc tới hạn. Sau đó ta thay đổi xuất phát điểm của vận tốc và giảm bước vận tốc. Q trình tiếp diễn đến khi đạt được kết quả đủ độ chính xác cần thiết Để đảm bảo độ tin cậy của tính tốn số, kết quả tính vận tốc flutter trong ví dụ 1 và 2 sẽ được so sánh với Li vcs 2011. Đoạn mã thể hiện việc tìm vận tốc flutter bằng cách giải phương trình vi phân trong ví dụ 1 được cho trong phụ lục Hình 39, 40 thể hiện kết quả tính vận tốc tới hạn theo các điều kiện đầu. Hình 39 là kết quả của luận án cịn hình 40 là kết quả trong Li vcs 2011 So sánh 2 hình vẽ ta thấy có sự trùng hợp rõ ràng và điều đó thể hiện độ tin cậy của các tính tốn số. Một số so sánh đnh l ị ượng cũng được chỉ ra trong bảng 11, 12 Kết quả so sánh trên các bảng cho thấy sự phù hợp và khẳng định độ tin cậy của các tính tốn số trong luận án 4.5. Kết quả tính tốn cho ví dụ 1 Hình 41 cho thấy sự so sánh giữa các đường cong biên độvận tốc trong ví dụ 1. Bảng 13 cho các kết quả định lượng so sánh các vận tốc tới hạn 20 flutter Kết quả cho thấy: Trên đường cong biên độ vận tốc, kỹ thuật đối ngẫu đầy đủ 3 thành phần cho kết quả gần với tính tốn số nhất Trên hình 41 và trên bảng 13, ta thấy tất cả các phương pháp đều thể hiện sự tồn tại của một vận tốc tới hạn flutter tối thiểu. Các phương pháp tuyến tính hóa đều dự báo tốt vận tốc tối thiểu này và các kết quả khơng khác nhau nhiều giữa các phương pháp (cột cuối cùng của bảng 13) Các phương pháp tuyến tính hóa đối ngẫu cải tiến đều cho kết quả cải thiện so với phương pháp tuyến tính hóa kinh điển trên miền phi tuyến (biên độ dao động lớn) như thấy trên cột 2 của bảng 13 4.6. Kết quả tính tốn cho ví dụ 2 Hình 42 cho thấy sự so sánh giữa các đường cong biên độvận tốc trong ví dụ 2. Bảng 14 cho các kết quả định lượng so sánh các vận tốc tới hạn flutter Kết quả tiếp tục cho thấy một số kết luận đã thu được ở ví dụ 1: Kỹ thuật tuyến tính hóa đối ngẫu đầy đủ 3 thành phần cho kết quả gần với tính tốn số nhất trong miền phi tuyến lớn (biên độ LCO lớn) Các kỹ thuật tuyến tính hóa đều tìm ra được vận tốc tới hạn flutter tối thiểu và sai khác rất ít so với tính tốn số 4.7. Kết quả tính tốn cho ví dụ 3 Trong ví dụ 3 này, ta sẽ tính hai trường hợp tham số K p khác nhau để đánh giá ảnh hưởng của điều khiển tỷ lệ (Pcontrol) đối với hệ trong ví dụ 1. Ta xét 2 trường hợp Kp=0.5 và Kp=0.5. Các kết quả so sánh được cho trên hình 43, 44 và bảng 15, 16 Qua các kết quả tính tốn, ta có một số nhận xét sau: Trong các trường hợp điều khiển tỷ lệ (Pcontrol), kỹ thuật tuyến tính hóa đối ngẫu cải tiến 3 (sử dụng 3 thành phần) có độ chính xác tốt nhất khi tính tốn vận tốc tới hạn ở miền có biên độ LCO lớn (độ phi tuyến cao) Các trường hợp tuyến tính hóa đều cho kết quả vận tốc tới hạn nhỏ nhất là gần như nhau và gần với kết quả tính tốn số 21 So sánh bảng 13 và bảng 15 ta thấy hệ số tỷ lệ Kp âm sẽ làm tăng vận tốc giới hạn tối thiểu (từ 7.94m/s lên 9.415m/s) và giảm vận tốc tới hạn ở miền biên độ lớn (từ 11.66m/s xuống 10.465m/s). Ngược lại, so sánh bảng 13 và bảng 16 ta thấy hệ số tỷ lệ Kp dương sẽ làm giảm vận tốc giới hạn tối thiểu (từ 7.94m/s xuống 6.985m/s) và tăng vận tốc tới hạn ở miền biên độ lớn (từ 11.66m/s lên 14.745m/s). Xét về mặt hiệu quả điều khiển thì hệ số Kp âm là tốt hơn 4.8. Kết quả tính tốn cho ví dụ 4 Trong ví dụ 4 này, ta sẽ tính hai trường hợp tham số K i khác nhau để đánh giá ảnh hưởng của điều khiển tích phân (Icontrol) đối với hệ trong ví dụ 1 Ta xét 2 trường hợp Ki=2 (1/s) và Ki=2 (1/s). Trường hợp Ki âm, hệ ổn định và các kết quả so sánh được cho trên hình 45 và bảng 17. Trường hợp Ki dương, hệ mất ổn định. Phương pháp số thể hiện điều này ở việc biên độ dao động tiến ra vơ cùng. Cịn các phương pháp tuyến tính hóa thể hiện điều này việc tồn tại giá trị riêng của bài tốn giá trị riêng (37) có phần thực dương. Các kết quả này được thể hiện trên hình 46 và bảng 18 Qua các kết quả tính tốn, ta có một số nhận xét sau: Trong các trường hợp điều khiển tích phân (Icontrol), kỹ thuật tuyến tính hóa đối ngẫu cải tiến 3 (sử dụng 3 thành phần) vẫn cho thấy độ chính xác tốt nhất khi tính tốn vận tốc tới hạn ở miền có biên độ LCO lớn (độ phi tuyến cao) Các trường hợp tuyến tính hóa đều cho kết quả vận tốc tới hạn nhỏ nhất là gần như nhau và gần với kết quả tính tốn số So sánh bảng 13 và bảng 17 ta thấy hệ số tích phân Ki âm sẽ làm giảm vận tốc giới hạn tối thiểu (từ 7.94m/s xuống 7.055m/s) và giảm vận tốc tới hạn ở miền biên độ lớn (từ 11.66m/s xuống 9.385m/s). Hệ số tích phân Ki dương sẽ làm hệ mất ổn định được thể hiện trên cả các tính tốn số và lời giải tuyến tính hóa. Như thấy trên bảng 18, giá trị riêng là số thực dương sẽ làm cho dao động tiến ra vơ cùng ngay cả với biên độ rất nhỏ Như vậy, lời giải tuyến tính hóa khơng chỉ xác định được các đáp ứng 22 LCO mà cịn chỉ ra được cả các trường hợp khơng tồn tại LCO, biên độ dao động tiến ra vơ cùng. Ngồi ra, qua ví dụ này ta thấy điều khiển tích phân (Icontrol) trong trường hợp này là khơng có hiệu quả (hoặc làm giảm vận tốc tới hạn tối thiểu, hoặc làm hệ mất ổn định) 4.9. Kết quả tính tốn cho ví dụ 5 Trong ví dụ 5 này, ta sẽ tính hai trường hợp tham số K d khác nhau để đánh giá ảnh hưởng của điều khiển vi phân (Dcontrol) đối với hệ trong ví dụ 1. Ta xét 2 trường hợp Kd=0.5 (s) và Kd=0.5 (s). Các kết quả so sánh được cho trên hình 47, 48 và bảng 19, 20 Qua các kết quả tính tốn, ta có một số nhận xét sau: Trong các trường hợp điều khiển vi phân (Dcontrol), kỹ thuật tuyến tính hóa đối ngẫu cải tiến 3 (sử dụng 3 thành phần) có độ chính xác tương đối tốt khi tính tốn vận tốc tới hạn ở miền có biên độ LCO lớn (độ phi tuyến cao) Các trường hợp tuyến tính hóa đều cho kết quả vận tốc tới hạn nhỏ nhất là gần như nhau và gần với kết quả tính tốn số. Đối với dao động có biên độ nhỏ hệ phi tuyến sẽ rất gần với hệ tuyến tính (tính phi tuyến yếu) khi đó phương pháp TTH kinh điển cho kết quả tốt hơn, tuy nhiên sai số của các phương pháp TTH đều nhỏ do tính phi tuyến yếu So sánh bảng 13 và bảng 19 ta thấy hệ số tỷ lệ Kd âm sẽ làm giảm vận tốc giới hạn tối thiểu (từ 7.94m/s xuống 2.805m/s) và giảm vận tốc tới hạn ở miền biên độ lớn (từ 11.66m/s xuống 2.87m/s). Ngược lại, so sánh bảng 13 và bảng 20 ta thấy hệ số tỷ lệ Kd dương sẽ làm tăng vận tốc giới hạn tối thiểu (từ 7.94m/s lên 9.5m/s) và giảm vận tốc tới hạn ở miền biên độ lớn (từ 11.66m/s xuống 9.78m/s) Như vậy trong cả hai trường hợp Kd dương và âm thì miền biến thiên của vận tốc giới hạn đều thu hẹp lại, tức là ít phụ thuộc vào biên độ dao động và tính phi tuyến ít đi Xét về mặt hiệu quả điều khiển thì hệ số Kd dương sẽ tốt hơn (tăng vận tốc giới hạn tối thiểu) 23 KẾT LUẬN VÀ KIẾN NGHỊ A. Các nội dung đã thực hiện trong luận án Đã tiến hành phân tích tổng quan một số kết quả chính ở trên thế giới và trong nước về vấn đề ổn định flutter của thiết diện cánh chịu lực khí động Trên cơ sở lý thuyết khí động học và các số liệu thực nghiệm đã có xây dựng mơ hình thiết diện cánh máy bay 2 chiều uốnxoắn chuyển động trong dịng khí khơng nén được, dựa trên mơ hình thí nghiệm được nhiều nhà nghiên cứu quốc tế sử dụng. Trong mơ hình được trình bày, độ cứng xoắn k được xét là một hàm của dưới dạng đa thức bậc 4. Phương trình phi tuyến thu được từ mơ hình được dùng để phân tích đáp ứng cũng như các hiện tượng flutter. Đã áp dụng các cơng cụ tính tốn động lực học dịng chảy (CFD v à CSM) để tính tốn và tối ưu hóa profile cánh máy bay. Phương pháp này được áp dụng trên các thiết kế của các cánh máy bay nhỏ, tại số Reynolds thấp. Các mơ phỏng được thực hiện để chứng minh hiệu quả của cánh máy bay tối ưu. Kết quả cho thấy các profile cánh máy bay cải tiến cho phép tăng tỉ số nâng trên cản, cũng như nâng cao hiệu quả khí động lực Sau khi thiết lập phương trình dao động hai bậc tự do của thiết diện cánh đã trình bày một số phương pháp giải tích, đặc biệt đã phát triển tiêu chuẩn đối ngẫu có trọng số cho phương pháp tuyến tính hóa tương đương hệ dao động phi tuyến tuần hồn và ngẫu nhiên, trong đó tiêu chuẩn tuyến tính hóa kinh điển là trường hợp riêng khi trọng số bằng khơng Áp dụng cho bài tốn ổn định flutter của thiết diện cánh đã thu được vận tốc U là hàm của biên độ dao động xoắn A. Từ đó sẽ vẽ được đường đặc trưng biên độ A với vận tốc tới hạn U và cũng cho ta mối quan hệ giữa tần số với vận tốc tới hạn. Để khảo sát cụ thể đã xét 2 ví dụ thiết diện cánh phi tuyến đến bậc 3 và bậc 5, khơng có điều khiển, 2 ví dụ thiết diện cánh phi tuyến đến bậc 5, có điều khiển tích phân và điều khiển vi phân. Đã xác định vận tốc tới hạn bằng phương pháp số giải hệ phương trình vi phân phi tuyến bằng hàm ode45 trong MATLAB. Đã xác định quy trình tính tốn số để xác định vận tốc tới hạn và tần số dao động LCO. Quy trình cần 24 làm từ bước thơ đến tinh. Cụ thể, đầu tiên bước vận tốc được lấy giá trị lớn để xác định vận tốc gần với vận tốc tới hạn. Sau đó thay đổi xuất phát điểm của vận tốc và giảm bước vận tốc. Q trình tiếp diễn đến khi đạt được kết quả đủ độ chính xác cần thiết. Đã đánh giá độ tin cậy của tính tốn số của luận án, kết quả tính vận tốc flutter trong ví dụ 1 và 2 được so sánh với Li vcs 2011. Kết quả so sánh trên đồ thị và bảng số liệu cho thấy có sự phù hợp rõ ràng và điều đó thể hiện độ tin cậy của các tính tốn số trong luận án. Kết quả phân tích cũng cho thấy về định tính tổng thể các kỹ thuật tuyến tính hóa và phương pháp số đều khơng có nhiều sự khác biệt trên đường cong biên độ vận tốc. Đối với dao động có biên độ nhỏ hệ phi tuyến sẽ rất gần với hệ tuyến tính (tính phi tuyến yếu) khi đó phương pháp TTH kinh điển cho kết quả tốt hơn, tuy nhiên sai số của các phương pháp TTH đều nhỏ do tính phi tuyến yếu. Đối với dao động có biên độ lớn (tính phi tuyến lớn) các kỹ thuật TTH đối ngẫu cho kết quả có sai số nhỏ hơn kỹ thuật TTH kinh điển. Ngồi ra, kỹ thuật đối ngẫu đầy đủ 3 thành phần thường cho kết quả gần với tính tốn số nhất B. Những kết quả mới của luận án bao gồm: Lần đầu tiên áp dụng và phát triển sáng tạo cách tiếp cận đối ngẫu để phân tích các hiện tượng dao động flutter xuất hiện trong thiết diện cánh chu l ị ực khí động Đối với tiêu chuẩn tuyến tính hóa có trọng số đã giải quyết vấn đề chọn giá trị trọng số thơng qua nghiên cứu đề xuất 3 cách lựa chọn tương ứng với 3 cải tiến Đã tiến hành khảo sát tần số dao động riêng trong hệ dao động bảo tồn bậc cao và hệ Duffing chịu tải trọng ngẫu nhiên cho thấy các cải tiến này đều cho kết quả tốt hơn so với phương pháp tuyến tính hóa kinh điển Áp dụng cho bài tốn ổn định flutter của thiết diện cánh, với các dạng phi tuyến đa thức hay gặp và tính phí tuyến khơng q nhỏ, cho thấy các tiêu chuẩn đối ngẫu cải tiến đều cho kết quả chính xác hơn tiêu chuẩn kinh điển. Ngồi ra các tiêu chuẩn đối ngẫu đầy đủ 3 thành phần thường cho 25 kết quả tốt hơn tiêu chuẩn đối ngẫu 2 thành phần trong các ví dụ khảo sát. Điều đó cho thấy tiêu chuẩn đối ngẫu đầy đủ 3 thành phần có thể áp dụng tốt cho các hệ có tính phi tuyến rõ rệt Hướng nghiên cứu tiếp theo sau luận án có thể tập trung vào một số việc sau đây: Đối với cách tiếp cận đối ngẫu trong phương pháp tuyến tính hóa tương đương cần khảo sát nhiều hệ dao động phi tuyến khác nhau để phát hiện các hiện tượng cũng như các khả năng và giới hạn áp dụng của tiêu chuẩn tương đương đối ngẫu Nghiên cứu mở rộng bài tốn ổn định thiết diện cánh 2 chiều sang bài tốn ổn định cánh có xét đến chiều dài cánh. Đây là bài tốn rất phức tạp nhưng rất cần thiết cho kỹ thuật hàng khơng. Nghiên cứu thêm các phương pháp CFD, CSM để áp dụng cho bài tốn ổn định cánh 26 ... mơ hình hóa thiết diện? ?cánh, các hiện tượng phi tuyến và? ?cách? ?tiếp? ?cận? ?đối ngẫu Chương 2: Xây dựng mơ hình thiết diện? ?cánh? ?máy? ?bay? ?chuyển động trong dịng khí khơng nén được.? ?Phân? ?tích? ?đáp? ?ứng? ?và các hiện tượng dao động flutter... xuất bởi Nguyễn Đơng Anh (2010).? ?Cách? ?tiếp? ?cận? ?mới được biết với tên gọi? ?cách? ?tiếp? ?cận? ?đối? ?ngẫu. Áp dụng? ?cách? ?tiếp? ?cận? ?này cho lĩnh vực tuyến tính hóa tương đương dẫn đến kỹ thuật cực tiểu bình phương? ?đối? ?ngẫu Cách? ?tiếp? ?cận? ?đối? ?ngẫu? ?có tính linh hoạt và có thể áp dụng được cho nhiều... ưu hóa? ?profile? ?cánh? ?máy? ?bay. Phương pháp này được áp dụng trên các thiết kế của các? ?cánh? ?máy? ?bay? ?nhỏ, tại số Reynolds thấp. Các mơ phỏng được thực hiện để chứng minh hiệu quả của? ?cánh? ?máy? ?bay? ?tối