بخشی از مقاله
*** این فایل شامل تعدادی فرمول می باشد و در سایت قابل نمایش نیست ***
حل معادله فنر دافینگ نرمال شده با استفاده از تئوری اغتشاش
چکیده
فنر دافینگ دارای معادلهای غیر خطی هست که مبتنی بر جرم فنر و رابطه کشیدگی و فشردگی آن میباشد. با نشان دادن این رابطه به طور کلی شناختی از این قضیه بهدست میآید. سپس به حل این معادله به روش های تئوری اغتشاش1و توازن هارمونیک ها2 پرداخته می گردد. با استفاده از تئوری اغتشاش به یک معادله غیر خطی رسیده و از توازن هارمونیکها در حل بقیه مسئله استفاده میشود.
با استفاده از نرم افزار Maple به جواب دلخواه رسیده و درباره منحصر به فرد بودن آن بحث می شود. در مرحله بعد با مقایسه پاسخ تحلیلی و عددی به بررسی خطای بین این دو پاسخ می پردازیم. خطای کم دو پاسخ کیفیت روش تئوری اغتشاش در حل مسئله غیر خطی (دافینگ) نشان می دهد.
واژههای کلیدی
تئوری اغتشاش، فنر دافینگ، توازن هارمونیک ها
-1 مقدمه
معادلهی فنر دافینگ (در حالت بدون ورودی ) برای بیان دینامیک جرم و فنر به صورت غیر خطی بیان می شود. هنگامی که مقدار نیروی فنر به مکان جرم وابسته باشد، فرکانس تشدید آن نیز به جابجایی فنر وابسته خواهد شد . نیروی فنر برحسب کشیدگی و فشردگی آن کم و زیاد می شود که برای فنر سخت با افزایش جابه جایی(کشیدگی ) این نیرو افزایش می یابد و برای فنر نرم برعکس می باشد .[1]
نوعی رفتار هیسترزیس در پاسخ فرکانسی سیستم وجود دارد که تحلیل این گونه سیستم ها را با مشکل مواجه می کند. به عبارت دیگر خروجی سیستم علاوه بر فرکانس ورودی به تغییرات فرکانس نیز وابسته است.
علاوه بر رفتار هیسترزیس فرکانسی در نمودارهای بد سیستم های غیرخطی، با تحریک سینوسی آنها با فرکانس معمولا هارمونیک های فرعی در آنها دیده می شود.
با افزایش ابعاد سیستم پیچیدگی در رفتار افزایش پیدا می کند. از این رو پاسخ فرکانسی برای سیستم های غیرخطی معمولا تعریف نمیشود. نظریه اغتشاش روشی ریاضی است که برای یافتن پاسخ تقریبی مسئله ای -که پاسخ دقیق آن قابل دسترسی نیست؛ استفاده میشود. یافتن این جواب تقریبی با یک پاسخ دقیق در یک مسئله مرتبط آغاز میشود.
نظریه اغتشاش را زمانی میتوان به کار برد که بتوان مسئله را با افزودن یک متغیر کوچک به توصیف ریاضی مسئلهای که قابل حل دقیق است، فرمول بندی کرد.
در سال 1386 با استفاده از تبدیلات هیلبرت- هوانگ نشان دادند که هر سیگنال ناپایا را میتوان به مودهای ذاتی خود تجزیه کرد. با این روش معادلات نوسانگر دافینگ را حل کرد. ]2[
در این زمینه همچنین می توان به مطالعات و استفاده از روش کنترل مد لغزشی برای حل معادلات و کنترل سیستم آشوبناک دافینگ اشاره کرد.[3]
معادلات دافینگ یک روش مناسب در حل معادلات غیرخطی در علوم مهندسی و فیزیک میباشد[4] که مدل اصلی آن توسط Duffing مهندس آلمانی در سال 1918 معرفی شد و توسط مون[5] و هلمز گسترش یافت[6] که به کم کردن خطای موجود پرداختند. [7] همچنین اسیلاتورهای شبه دافینگ توسط هایاشی[8] و یودا[9] بررسی شد.
در این مقاله سعی بر حل معادلات فنر دافینگ شده با استفاده از تئوری اغتشاش می باشد.
-2بیان مسئله:
معادلهی دافینگ برای بیان دینامیک جرم و فنر به صورت زیر میباشد.
که جابهجایی فنر، دو ضریب مثبت دلخواه و ترم سوم و چهارم معادله ( 1)، نیروی فنر را بیان می کند. پارامتر ضریب کیفیت و مشتقات برحسب زمان میباشد. با انتخاب مقادیر بزرگ برای میتوان را برابر با ε در نظر گرفت. با استفاده از تئوری اغتشاش به صورت زیر تعریف میشود:
با محاسبه مقادیر مشتقات زمانی و عبارات توانی برحسب معادله )2( و جاگذاری آنها در معادله )1( و ساده سازی به دست میآید:
برای به دست آوردن مقادیر معادله غیرخطی در نگاه اولیه قابل حل نمیباشد. از این رو با نرم افزار Maple و با کمک گرفتن از تئوری توازن هارمونیکها (Harmonic Balance) طبق مراحل زیر به جواب رسیده، برای سادگی حل مقدار پارامتر دلخواه را یک فرض کرده و سپس در مراحل بعد مقدار ثابتی در نظر میگیریم. ضریب را ابتدا وارد کرده و با استفاده از تئوری توازن هارمونیکها پاسخ زیر را برای تعریف میشود:
حال با جایگذاری این مقدار در معادله و حل معادله برای چند مقدار مختلف روند کلی شکل و معادله مشخص شده و مقدار آن معین میشود. که مقادیر ضرایب برای معادله )4( به شرح زیر میباشد.
ضرایب بهدست آمده از فرمول )5( را در فرمول )4( قرار داده و مقدار بهدست میآید.
سپس برای با طی کردن روند بالا داریم:
ضرایب مانند روشی که برای ضرایب بهدست میآیند که این مقادیر به شرح زیر میباشند:
بنابراین با جایگذاری مقادیر )8( در )7( مقدار به صورت زیر به دست میآید.
حال برای 2 هم به طریقه مذکور میتوانیم داشته باشیم:
بنابراین هم به شکل زیر به دست میآید:
حال به رسم نمودارهای مربوط به مقادیر با کمک نرمافزار Maple پرداخته می شود.
