خانه » پروژه » پروژه مکانیک » دانلود پروژه طراحی زیر دریایی و بدنه ایرشیپ ها
دانلود پروژه طراحی زیر دریایی و بدنه ایرشیپ ها

دانلود پروژه طراحی زیر دریایی و بدنه ایرشیپ ها

طراحی زیر دریایی و بدنه ایرشیپ ها

فهرست مطالب
۳\طراحی بدنه ایرشیپ‌ها و زیر دریائی‌ها    ۱
فهرست علائم    ۳
فصل اول مقدمه و مروری بر تحقیقات گذشته    ۱۳
۱-۱ مقدمه و مروری بر تحقیقات گذشته    ۱۴
شکل ۱-۱ پروفیلهای بدنه با کمترین درگ    ۱۶
۱-۱-۱ مدل آیرودینامیکی    ۱۷
فصل دوم معادلات حاکم وروش حل عددی    ۱۹
۲-۱  مقدمه    ۲۰
۲-۲ محاسبات لایه مرزی    ۲۱
۲-۲-۱  محاسبات لایه مرزی آرام    ۲۱
۲-۲-۲ محاسبات ناحیه گذرا    ۲۳
۲-۲-۳ محاسبات لایه مرزی درهم    ۲۴
۲-۲-۴ روش محاسبه درگ    ۲۶
۲-۲-۵  معیار جدایش    ۲۷
فصل سوم الگوریتم و برنامه به همراه ورودی و خروجی‌های برنامه    ۲۸
۳-۱ روند محاسبات درگ    ۲۹
۳-۲ الگوریتم محاسبات لایه مرزی آرام:    ۳۰
۳-۳ الگوریتم محاسبات ناحیه گذرا    ۳۰
۳-۴الگوریتم محاسبات لایه مرزی درهم وضریب درگ    ۳۱
۳-۵ برنامه کامپیوتری به زبان فرترن    ۳۲
۳-۶ ورودی‌ها و خروجی‌های برنامه برای پروفیل‌های ۱بدنه شماره۱ تا ۷    ۴۳
۳-۶-۱ ورودی برنامه برای پروفیل بدنه شماره(۱)    ۴۵
۳-۶-۲ خروجی برنامه برای پروفیل بدنه شماره (۱)    ۴۸
۳-۶-۳ ورودی برنامه برای پروفیل بدنه شماره(۲)    ۵۲
۳-۶-۴ خروجی برنامه برای پروفیل بدنه شماره(۲)    ۵۵
شکل ۳-۳   پروفیل بدنه شماره (۳)    ۵۸
۳-۶-۵ ورودی برنامه برای پروفیل بدنه شماره(۳)    ۵۹
۳-۶-۶ خروجی برنامه برای پروفیل بدنه شماره(۳)    ۶۲
۳-۶-۷ ورودی برنامه برای پروفیل بدنه شماره(۴)    ۶۶
۳-۶-۸ خروجی برنامه برای پروفیل بدنه شماره(۴)    ۶۹
۳-۶-۹ ورودی برنامه برای پروفیل بدنه شماره(۵)    ۷۳
۳-۶-۱۰ خروجی برنامه برای پروفیل بدنه شماره(۵)    ۷۶
۳-۶-۱۱ ورودی برنامه برای پروفیل بدنه شماره(۶)    ۸۰
۳-۶-۱۲ ورودی‌های برنامه برای پروفیل بدنه شماره(۷)    ۸۴
۳-۶-۱۳ خروجی برنامه برای پروفیل بدنه شماره ۶ و۷    ۸۶
فصل چهارم ارائه نتایج و بحث و مقایسه    ۸۷
۴-۱ مقدمه    ۸۸
۴-۲ نتایج و بحث برای پروفیل شماره ۱    ۸۸
۴-۳ نتایج و بحث برای پروفیل شماره ۲    ۸۹
۴-۴ نتایج و بحث برای پروفیل شماره ۳    ۹۰
۴-۵ نتایج و بحث برای پروفیل شماره ۴    ۹۱
۴-۶ نتایج و بحث برای پروفیل شماره ۵    ۹۲
۴-۷ نتایج و بحث برای پروفیل شماره ۶و۷    ۹۳
۴-۸ نمودارهای مربوط به پروفیل بدنه شماره ۱    ۹۴
شکل۴-۳ منحنی تغییرات ضخامت جابجایی بدون بعد بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۹۵
شکل۴-۵ منحنی تغییرات ضریب اصطکاک سطحی بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۹۶
۴-۹ نمودارهای مربوط به پروفیل بدنه شماره ۲    ۹۸
شکل۴-۷ منحنی تغییرات ضخامت ممنتوم بدون بعد بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۹۸
شکل۴-۸ منحنی تغییرات ضریب شکل بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۹۹
شکل ۴-۹  منحنی تغییرات ضخامت جابجایی بدون بعد بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۰۰
شکل ۴-۱۰  منحنی تغییرات ضخامت لایه مرزی بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۰۱
شکل ۴-۱۱  منحنی تغییرات ضریب اصطکاک سطحی بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۰۲
شکل۴-۱۲  منحنی تغییرات عدد رینولدز(که براساس ضخامت مومنتوم تعریف شده)  بر حسب طول بدون بعد ایرشیپ در ناحیه درهم    ۱۰۳
۴-۱۰ نمودارهای مربوط به پروفیل بدنه شماره ۳    ۱۰۴
شکل ۴-۱۳  منحنی تغییرات ضخامت ممنتوم بدون بعد بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۰۴
شکل ۴-۱۴ منحنی تغییرات ضریب شکل بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۰۵
شکل ۴-۱۵ منحنی تغییرات ضخامت جابجایی بدون بعد بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۰۶
شکل ۴-۱۶  منحنی تغییرات ضخامت لایه مرزی بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۰۷
شکل ۴-۱۷  منحنی تغییرات ضریب اصطکاک سطحی بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۰۸
شکل ۴-۱۸  منحنی تغییرات عدد رینولدز(که براساس ضخامت مومنتوم تعریف شده) بر حسب طول بدون بعد ایرشیپ در ناحیه درهم    ۱۰۹
۴-۱۱ نمودارهای مربوط به پروفیل بدنه شماره ۴    ۱۱۰
شکل ۴-۱۹  منحنی تغییرات ضخامت ممنتوم بدون بعد بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۱۰
شکل۴-۲۰ منحنی تغییرات ضریب شکل بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۱۱
شکل ۴-۲۱  منحنی تغییرات ضخامت جابجایی بدون بعد بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۱۲
شکل۴-۲۲  منحنی تغییرات ضخامت لایه مرزی بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۱۳
شکل ۴-۲۳  منحنی تغییرات ضریب اصطکاک سطحی بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۱۴
۴-۱۲ نمودارهای مربوط به پروفیل بدنه شماره ۵    ۱۱۶
شکل ۴-۲۵  منحنی تغییرات ضخامت ممنتوم بدون بعد بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۱۶
شکل ۴-۲۶ منحنی تغییرات ضریب شکل بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۱۷
شکل۴-۲۷  منحنی تغییرات ضخامت جابجایی بدون بعد بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۱۸
شکل ۴-۲۸  منحنی تغییرات ضخامت لایه مرزی بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۱۹
شکل ۴-۲۹  منحنی تغییرات ضریب اصطکاک سطحی بر حسب طول بدون بعد ایرشیپ درناحیه درهم    ۱۲۰
۴-۱۳ مقایسه ضریب درگ    ۱۲۲
فصل پنجم نتیجه گیری و پیشنهادات    ۱۲۴
۱- ۵ نتیجه گیری    ۱۲۵
۲-۵ پیشنهادات    ۱۲۵
فهرست مراجع    ۱۲۶
مراجع فارسی    ۱۲۷
فهرست علائم
تعریف    علائم اصلی
ضریب درگ
ضریب اصطکاک سطحی
قطر ماکزیمم بدنه (cm )
نیروی درگ
ضریب شکل
طول کلی بدنه(cm)
فشار
عدد رینولدز بر اساس  طول کلی بدنه
شعاع ماکزیمم بدنه(cm)
شعاع محلی بدنه (cm)
شعاع محلی بدون بعد بدنه
عدد رینولدز بر اساس ضخامت مومنتوم
عدد رینولدز حجمی
سطح تصویر شده بدنه بر اساس شعاع ماکزیمم(cm2)
مولفه بردار سرعت در راستای x (cm/s)
سرعت روی لبه لایه مرزی(cm/s)
سرعت در نقطه سکون(cm/s)
سرعت بدون بعد روی لبه لایه مرزی
سرعت بدون بعد در نقطه سکون
سرعت جریان آزاد (cm/s)
مولفه بردار سرعت در راستایy (cm/s)
مولفه قائم سرعت روی لبه لایه مرزی(cm/s)
محور مختصات موازی سطح بدنه(cm)
محور مختصات عمود برسطح بدنه(cm)
تعریف    علائم یونانی
ضخامت لایه مرزی(cm)
ضخامت جابجایی
ضخامت مومنتوم
چگالی
تنش برشی روی دیواره
ویسکوزیته سینماتیکی( )
فصل اول
مقدمه و مطالعات پیشین
۱-۱ مقدمه و مروری بر تحقیقات گذشته
۱-۱-۱ مدل آیرودینامیکی
فصل دوم
معادلات حاکم و روش حل عددی
۲-۱ مقدمه
۲-۲ محاسبات لایه مرزی
۲-۲-۱ محاسبات لایه مرزی آرام
۲-۲-۲ محاسبات ناحیه گذرا
۲-۲-۳ محاسبات لایه مرزی درهم
۲-۲-۴ روش محاسبه درگ
۲-۲-۵ معیار جدایش

فصل سوم
الگوریتم و برنامه به همراه ورودی و خروجی های برنامه
۳-۱ روند محاسبه درگ
۳-۲ الگوریتم محاسبات لایه مرزی آرام
۳-۳ الگوریتم محاسبات ناحیه گذرا
۳-۴ الگوریتم محاسبات لایه مرزی درهم و ضریب درگ
۳-۵ برنامه کامپیوتری به زبان فرترن
۳-۶ ورودی و خروجی های برنامه برای پروفیل های بدنه شماره ۱ تا ۷
۳-۶-۱ ورودی برنامه برای پروفیل بدنه شماره ۱
۳-۶-۲ خروجی برنامه برای پروفیل بدنه شماره ۱
۳-۶-۳ ورودی برنامه برای پروفیل بدنه شماره ۲
۳-۶-۴ خروجی برنامه برای پروفیل بدنه شماره ۲
۳-۶-۵ ورودی برنامه برای پروفیل بدنه شماره ۳
۳-۶-۶ خروجی برنامه برای پروفیل بدنه شماره ۳
۳-۶-۷  ورودی برنامه برای پروفیل بدنه شماره ۴
۳-۶-۸ خروجی برنامه برای پروفیل بدنه شماره ۴
۳-۶-۹  ورودی برنامه برای پروفیل بدنه شماره ۵
۳-۶-۱۰  خروجی برنامه برای پروفیل بدنه شماره ۵
۳-۶-۱۱ ورودی برنامه برای پروفیل بدنه شماره ۶
۳-۶-۱۲ ورودی برنامه برای پروفیل بدنه شماره ۷
۳-۶-۱۳ خروجی برنامه برای پروفیل بدنه شماره ۶و۷

فصل چهارم
ارائه نتایج و بحث و مقایسه
۴-۱ مقدمه
۴-۲ نتایج و بحث برای پروفیل بدنه شماره ۱
۴-۳ نتایج و بحث برای پروفیل بدنه شماره ۲
۴-۴ نتایج و بحث برای پروفیل بدنه شماره ۳
۴-۵ نتایج و بحث برای پروفیل بدنه شماره ۴
۴-۶ نتایج و بحث برای پروفیل بدنه شماره ۵
۴-۷ نتایج و بحث برای پروفیل بدنه شماره ۶و۷
۴-۸ نمودارهای مربوط به پروفیل بدنه شماره ۱
۴-۹ نمودارهای مربوط به پروفیل بدنه شماره ۲
۴-۱۰ نمودارهای مربوط به پروفیل بدنه شماره ۳
۴-۱۱ نمودارهای مربوط به پروفیل بدنه شماره ۴
۴-۱۲ نمودارهای مربوط به پروفیل بدنه شماره ۵
۴-۱۳ مقایسه ضریب درگ
فصل پنجم
نتیجه گیری و پیشنهادات
۵-۱ نتیجه گیری
۵-۲ پیشنهاداتی برای تحقیقات آینده
فهرست مراجع
پیوست”الف”
فهرست جداول
عنوان    صفحه
جدول ۳-۱ ورودی‌های برنامه برای پروفیل بدنه شماره ۱
جدول ۳-۲ خروجی‌های برنامه برای پروفیل بدنه شماره ۱
جدول ۳-۳ ورودی‌های برنامه برای پروفیل بدنه شماره ۲
جدول ۳-۴ خروجی‌های برنامه برای پروفیل بدنه شماره ۲
جدول ۳-۵ ورودی‌های برنامه برای پروفیل بدنه شماره ۳
جدول ۳-۶ خروجی‌های برنامه برای پروفیل بدنه شماره ۳
جدول ۳-۷ ورودی‌های برنامه برای پروفیل بدنه شماره ۴
جدول ۳-۸ خروجی‌های برنامه برای پروفیل بدنه شماره ۴
جدول ۳-۹ ورودی‌های برنامه برای پروفیل بدنه شماره ۵
جدول ۳-۱۰ خروجی‌های برنامه برای پروفیل بدنه شماره ۵
جدول ۳-۱۱ ورودی‌های برنامه برای پروفیل بدنه شماره ۶
جدول ۳-۱۲ ورودی‌های برنامه برای پروفیل بدنه شماره ۷
جدول ۴-۱ ضریب درگ برای پروفیل‌های بدنه یک تا پنج
فهرست اشکال
عنوان صفحه
شکل ۱-۱ پروفیلهای بدنه با کمترین درگ
شکل ۱-۲ مدل آیرودینامیکی
شکل ۱-۳ توزیع المانهای سینگولاریتی محوری و شدت در۲۱ نقطه طول بدنه
شکل ۳-۱   پروفیل بدنه شماره ۱
شکل ۳-۲   پروفیل بدنه شماره ۲
شکل ۳-۳   پروفیل بدنه شماره ۳
شکل ۳-۴   پروفیل بدنه شماره ۴
شکل ۳-۵   پروفیل بدنه شماره ۵
شکل ۳-۶   پروفیل بدنه شماره ۶
شکل ۳-۷   پروفیل بدنه شماره ۷
شکل۴-۱ منحنی تغییرات ضخامت ممنتوم بر حسب طول ایرشیپ درناحیه درهم  برای پروفیل شماره ۱
شکل۴-۲ منحنی تغییرات ضریب شکل بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۱
شکل۴-۳ منحنی تغییرات ضخامت جابجایی بدون بعد بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۱
شکل۴-۴ منحنی تغییرات ضخامت لایه مرزی بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۱
شکل۴-۵ منحنی تغییرات ضریب اصطکاک سطحی بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۱
شکل۴-۶ منحنی تغییرات عدد رینولدز(که براساس ضخامت مومنتوم تعریف شده) بر حسب طول ایرشیپ برای پروفیل شماره ۱
شکل۴-۷ منحنی تغییرات ضخامت ممنتوم بر حسب طول ایرشیپ درناحیه درهم  برای پروفیل شماره ۲
شکل۴-۸ منحنی تغییرات ضریب شکل بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۲
شکل۴-۹ منحنی تغییرات ضخامت جابجایی بدون بعد بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۲
شکل۴-۱۰ منحنی تغییرات ضخامت لایه مرزی بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۲
شکل۴-۱۱ منحنی تغییرات ضریب اصطکاک سطحی بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۲
شکل۴-۱۲ منحنی تغییرات عدد رینولدز(که براساس ضخامت مومنتوم تعریف شده) بر حسب طول ایرشیپ برای پروفیل شماره ۲
شکل۴-۱۳ منحنی تغییرات ضخامت ممنتوم بر حسب طول ایرشیپ درناحیه درهم  برای پروفیل شماره ۳
شکل۴-۱۴ منحنی تغییرات ضریب شکل بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۳
شکل۴-۱۵ منحنی تغییرات ضخامت جابجایی بدون بعد بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۳
شکل۴-۱۶ منحنی تغییرات ضخامت لایه مرزی بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۳
شکل۴-۱۷ منحنی تغییرات ضریب اصطکاک سطحی بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۳
شکل۴-۱۸ منحنی تغییرات عدد رینولدز(که براساس ضخامت مومنتوم تعریف شده) بر حسب طول ایرشیپ برای پروفیل شماره ۳
شکل۴-۱۹ منحنی تغییرات ضخامت ممنتوم بر حسب طول ایرشیپ درناحیه درهم  برای پروفیل شماره ۴
شکل۴-۲۰ منحنی تغییرات ضریب شکل بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۴
شکل۴-۲۱ منحنی تغییرات ضخامت جابجایی بدون بعد بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۴
شکل۴-۲۲ منحنی تغییرات ضخامت لایه مرزی بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۴
شکل۴-۲۳ منحنی تغییرات ضریب اصطکاک سطحی بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۴
شکل۴-۲۴ منحنی تغییرات عدد رینولدز(که براساس ضخامت مومنتوم تعریف شده) بر حسب طول ایرشیپ برای پروفیل شماره ۴
شکل۴-۲۵ منحنی تغییرات ضخامت ممنتوم بر حسب طول ایرشیپ درناحیه درهم  برای پروفیل شماره ۵
شکل۴-۲۶ منحنی تغییرات ضریب شکل بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۵
شکل۴-۲۷ منحنی تغییرات ضخامت جابجایی بدون بعد بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۵
شکل۴-۲۸ منحنی تغییرات ضخامت لایه مرزی بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۵
شکل۴-۲۹ منحنی تغییرات ضریب اصطکاک سطحی بر حسب طول ایرشیپ درناحیه درهم برای پروفیل شماره ۵
شکل۴-۳۰ منحنی تغییرات عدد رینولدز(که براساس ضخامت مومنتوم تعریف شده) بر حسب طول ایرشیپ برای پروفیل شماره ۵
شکل ۴-۳۱ نتایج بدست آمده توسط لوتز و واگنر برای ضریب درگ به روش اپلر
فصل اول مقدمه و مروری بر تحقیقات گذشته
۱-۱ مقدمه و مروری بر تحقیقات گذشته
در طراحی بدنه ایرشیپ‌ها و زیر دریائی‌ها نکات زیادی مورد توجه قرار می‌گیرد که مهمترین آنها قدرت جلوبرندگی است که به مقدار زیادی بستگی به درگ اصطکاکی روی بدنه  ایرشیپ دارد و ۳/۲ درگ کل را شامل می‌شود. کاهش کوچکی در این درگ باعث صرفه جویی قابل توجهی در سوخت می‌شود و یا می‌تواند باعث افزایش ظرفیت حمل و ابعاد ایرشیپ شود.
اولین بهینه سازی عددی شکل، توسط  پارسنز   انجام شده است. روش محاسبه در قالب یک پنل کد  می‌باشد  که با یک  روش لایه مرزی کوپل شده است. زدان   یک توزیع محوری از چشمه و چاه را برای نشان دادن میدان جریان اطراف یک جسم معرفی  می‌کند. قدرت (شدت) به صورت خطی روی  هر المان طول توزیع می‌شود.
در روند محاسباتی آیرودینامیکی ابتدا یک بدنه دوار با ماکزیمم قطر ثابت و نسبت فایننس    ثابت تعریف می‌شود.پروفیل بدنه و توزیع سرعت  جریان غیر لزج توسط روشهای غیر مستقیم حل جریان پتانسیل بدست می‌آید. پروفیل این بدنه باید  به گونه‌ای باشد که در جریان یکنواخت موازی با  محور بدنه، لایه مرزی دچار جدایش نشود. با این قید، درگ توسط تغییر در شکل پروفیل بدنه کاهش می‌یابد. محدودیت در عدم جدایش لایه مرزی باعث حذف درگ  فشاری می‌شود و درگ کلی منحصر به نیروهای ویسکوز در لایه مرزی می‌شود. لایه مرزی به سه ناحیه آرام گذرا    و درهم تقسیم می‌شود. برای محاسبه لایه مرزی آرام از  متد توویتس   استفاده شده که  بر اساس رابطۀ مومنتوم می‌باشد. ناحیه گذرا در محاسبات به صورت  یک نقطه در نظر گرفته می‌شود که در آن ضریب شکل به طور ناگهانی از آخرین مقدار در ناحیه آرام به اولین مقدار  در ناحیه درهم تغییر می‌کند. از آنجا که محل گذر به عواملی مانند: زبری سطحی، سر و صدا، لرزش و غیره بستگی دارد که کنترل آنها مشکل است در بیشتر تحقیقات این ناحیه را به  صورت دلخواه بین سه تا ده درصد طول بدنه در نظر می‌گیرند.
محاسبات لایه مرزی مغشوش بر اساس یک روش ساده انتگرالی معادله مومنتوم بنا شده است، که توسط شینبروک    و سامنر    برای جریان با تقارن محوری بدست آمده است. از آنجا که لایه مرزی مجاز به جدایش نیست درگ از نقصان مومنتوم در انتهای لایه مرزی  محاسبه می‌شود.
حل این مسأله در ساخت اژدرها، زیر دریائی‌ها و ایرشیپ‌ها مورد استفاده قرار می‌گیرد. بعضی از این گونه‌ها پروفیل بدنه را به صورت یک یا دو چند جمله‌ای از درجات مختلف نشان می‌دهند و شامل پارامترهایی مانند شعاع در دماغه و انتهای دم محل نسبی قطر ماکزیمم و شعاع طولی در آن نقطه و شیب دم هستند. بوسیله تغییر در بعضی یا همه این پارامترها در شکلهای مختلف درگ کاهش یافته است. دیگران سعی کرده‌اند که مستقیما  از کپی پروفیل بدنه ماهی‌های پرسرعت و پرندگان این کار را دنبال کنند. نتیجه تمام این تلاشها منجر به طبقه بندی  بدنه هایی با درگ پایین شده است و گرچه از نظر شکل متفاوت هستند ولی ضریب درگهایی خیلی شبیه به هم دارند این بدنه‌ها در شکل ۱-۱  آمده است.
شکل ۱-۱ پروفیلهای بدنه با کمترین درگ
۱-۱-۱ مدل آیرودینامیکی
جریان اطراف بدنه ایرشیپ با زاویه حمله صفر را به کمک روش سوپر پوزیشن  بر روی یک سری توزیع چشمه و چاه که روی محور بدنه و بصورت المانهایی بطول    و با توزیع شدتی   که توسط یک پاره خط مستقیم و روی المان قرار دارد تخمین می‌زنیم.
تابع جریان این المان در نقطه i  به شکل زیر است:
(۱-۱)
)۱-۲)
پروفیل بدنه از طریق مساوی قرار دادن تابع جریان برابر با صفر وحل آن برای در تعداد مشخصی از نقاط با فاصله مساوی مثلا” برای ۲۰ المان بدست می‌آید شکل (۱-۲).
شکل ۱-۲ مدل آیرودینامیکی
خط محوری چشمه و چاه به ۲۰ المان با طول مساوی و در نتیجه به ۲۱ نقطه انتهایی تقسیم می‌شودکه هر المان توزیع شدت خطی دارد (شکل۱-۳).با مشخص کردن شدت‌ها در ۲۱ نقطه انتهایی توزیع شدت در همه جا تعریف شده است. پروفیل بدنه بوسیله ی تغییر در مقدار شدت این ۲۱ نقطه انتهایی تغییر می‌کند. ترکیبات جدیدی از این ۲۱ شدت تولید می‌شود که در قالب پایان نامه کارشناسی ارشد رضا حسن زاده ارائه شده است. ضریب درگ با استفاده از محاسبات لایه مرزی در نزدیک سطح بدنه بدست می‌آید که   محاسبات لایه مرزی آرام و درهم و همچنین ناحیه گذرا که در این تحقیق بررسی می‌شود بطور مفصل در قسمتهای بعدی شرح داده خواهد شد.
این بدنه جدید به عنوان مبنا قرار می‌گیرد و می‌تواند در یک پروسه ی تکاملی بهینه سازی  شود تا به پروفیل با کمترین درگ دست یابیم.در چهل سال اخیر سیستم‌های حل مسأله ی بهینه سازی که بر اساس تکامل و وراثت بنا شده‌اند مورد توجه قرار گرفتند،استراتژی تکامل ریخنبرگ ]۶ [یکی از این روش‌ها می‌باشد.روش قدرتمند دیگری که بر پایه تکنیک‌های هوش مصنوعی می‌باشد و قابل استفاده در فضا‌های عملکرد بزرگ و توابع چند بعدی و چند وضعیتی (دارای چندین می‌نیمم)و غیر خطی می‌باشد، روش الگوریتم ژنتیک  است.
شکل ۱-۳ توزیع المانهای سینگولاریتی محوری و شدت در ۲۱ نقطه طول بدنه
فصل دوم معادلات حاکم وروش حل عددی
۲-۱  مقدمه
مقاومت ویسکوز بدنه اغلب از حل لایه مرزی محاسبه می‌شود که برای حل لایه مرزی نیاز به دانستن توزیع سرعت در لبه لایه مرزی  می‌باشد که از حل جریان پتانسیل بدست می‌آید. لایه مرزی به سه قسمت آرام،گذرا و درهم تقسیم می‌شود. براساس معادله مومنتوم در شرایط جریان پایدار،دوبعدی،تراکم ناپذیر وویسکوز با گرادیان فشار در جهت x داریم:
معادله پیوستگی:
(۲-۱)
معادله ممنتوم:
(۲-۲)
که  به صورت تابعی از x  و مستقل از y به شکل زیر می‌باشد:
(۲-۳)
با ضرب کردن معادله(۲-۱) در   و جمع آن با معادله (۲-۲) وانتگرال گیری در محدوده ضخامت لایه مرزی به معادله دیفرانسیل رایج برحسب ضخامت مومنتوم   می‌رسیم:
(۲-۴)
که ضخامت مومنتوم به صورت زیر تعریف می‌شود:
(۲-۵)
و ضریب شکل  به فرم زیر می‌باشد:
(۲-۶)
و ضخامت جابجایی   به فرم زیر می‌باشد:
(۲-۷)
و ضریب اصطکاک سطحی  به صورت زیر می‌باشد:
(۲-۸)
۲-۲ محاسبات لایه مرزی
۲-۲-۱  محاسبات لایه مرزی آرام
برای محاسبه لایه مرزی آرام در روی ایرشیپ و با تقارن محوری از متد تویتس   ]۴[ استفاده می‌شود که  جزئیات روابط جبری روش تویتس را می‌توان در سبسی و برد شاو  [۸] مشاهده نمود.در روش تویتس شرایط مرزی به صورت زیر می‌باشد:
(۲-۹)

با قرار دادن y=0  دررابطه (۲-۲) و استفاده از رابطه (۲-۹) داریم:
(۲-۱۰)
که λ یک ترکیب مناسب از  و  می‌باشد.  تویتس فرض کرد که   و H  فقط توابعی از λ می‌باشند
و از روابط تجربی استفاده کرد. برای ضریب اصطکاک سطحی از رابطه (۲-۸)و (۲-۹) داریم:
(۲-۱۱)
حال با جایگزینی رابطه(۲-۹)و(۲-۱۱)در رابطه (۲-۴)و ضرب کردن طرفین در  ومرتب نمودن آن داریم:

که با انتگرال گیری داریم:
(۲-۱۲)
با تعریف مقادیر بدون بعد به صورت زیر:
(۲-۱۳)
لذا رابطه (۲-۱۲)به صورت زیر تبدیل می‌شود:
(۲-۱۴)
برای نقطه سکون رابطه(۲-۱۴) به صورت زیر در می‌آید:
(۲-۱۵)
که عبارت   مشخص کننده شیب توزیع سرعت خارجی  برای نقطه سکون جریان می‌باشد جمله آخر رابطه (۲-۱۴) به علت اینکه  =  ، در نقطه سکون صفر می‌شود.
در حالیکه   از توزیع سرعت خارجی محاسبه می‌شود پارامترهای دیگر لایه مرزی مانند و   از روابط تجربی زیر به دست می‌آیند:
برای مقادیر مثبت  :
(۲-۱۶)
برای مقادیر منفی  :
(۲-۱۷)
با استفاده از تبدیل مانگلار ]۴[ می‌توان از  روش تویتس برای لایه مرزی روی اجسام با تقارن محوری به صورت زیر استفاده کرد:
(۲-۱۸)
که در این رابطه داریم:
۲-۲-۲ محاسبات ناحیه گذرا
پیش بینی تئوری ناحیه‌ای که گذر از لایه مرزی آرام به درهم رخ می‌دهد، به عنوان یکی ازمسائل پیچیده و مشکل در مکانیک سیالات می‌باشد زیرا ناحیه گذرا به فاکتورهای زیادی مــانند سروصـدا،لرزش، محیـط، زبری سطحی بدنه وگرادیان فشار سطحی بستگی دارد که تعیین اثرات آنها روی ناحیه گذرا مشکل است. اولین تحقیقات جدی در این زمینه در اواخر قرن نوزدهم وتوسط رینولدز  صورت گرفت.تحقیقات دیگری توسط گرانویل ، کربتری  صورت گرفت و به خاطر ناتوانی این متد‌ها در بیان تاثیرات سطح بدنه ومحیط روی پدیده گذر تعدادی از محققان به صورت دلخواه ناحیه گذرا را بین سه تا ده درصد طول بدنه از دماغه در نظرگرفتند که در این روش نیز از همین تجربه استفاده شده است.  ناش  این ناحیه را به صورت یک نقطه ودرسه درصد طول بدنه فرض کرده است. در ناحیه گذرا چند تغییر اساسی در لایه مرزی رخ می‌دهد.این تغییرات به صورت تغییر در ضخامت جابجایی و ضخامت مومنتوم نشان داده می‌شودکه منجر به کاهش ضریب شکل می‌شود. باجایگزین کردن ناحیه گذر به صورت یک نقطه ناش توانست روش مفیدی برای محاسبه مقادیر   و  در آغاز لایه مرزی آرام بدست آورد.مقدار  در طول ناحیه گذر تغییر نمی‌کند در حالیکه مقدار   در شروع لایه مرزی درهم از رابطه تعادلی ناش بدست می‌آید.
معادلات تجربی ناش] ۹ [مربوط به نقطه گذرا با استفاده از مقادیر   و  و  بصورت زیر می‌باشند که در یک روند تکراری قابل حل خواهد بود.
(۲-۱۹)
(۲-۲۰)
(۲-۲۱)
(۲-۲۲)
(۲-۲۳) رابطه لودویگ – تیلمن [۱۲]

پروسه تکراری حل با مقدار اولیه ۵/۱= آغاز می‌شود و سپس مقدار  ازرابطه (۲-۲۳) محاسبه شده و بعد از آن به ترتیب   و  و   محاسبه می‌شود و نهایتاً    از معادله (۲-۲۲) برای شروع مجدد سیکل مورد استفاده قرار می‌گیرد تا دو مقدار متوالی    همگرا شوند.
۲-۲-۳ محاسبات لایه مرزی درهم
روش‌های زیادی برای محاسبه گرادیان فشار درلایه مرزی درهم وجود داردکه همگی از فرم انتگرالی معادله مومنتوم ورابطه‌ای که  را به عدد رینولدزو ضریب شکل پروفیل مرتبط می‌کند وهمچنین یک معادله دیفرانسیل برای نرخ تغییرات ضریب شکل بر حسب   استفاده می‌کنند.
در رابطه مربوط به ضریب شکل، کل اطلاعات راجع به تنشهای رینولدز درون لایه مرزی موجود است.یکی از این روشهای موفق روش انتگرالی هد می‌باشد. محاسبات لایه مرزی درهم براساس روش انتگرالی هد  ] ۴[ برای جریان دوبعدی بنا شده است.این متد با  همان فرم انتگرالی معادله مومنتوم شروع می‌شود که در محاسبات لایه مرزی آرام توسط تویتس استفاده شده بود:
(۲-۲۴)
که در این معادله سه مقدار   مجهول می‌باشند.
هد سرعت بدون بعد زیر را به صورت تابعی از ضریب شکل تعریف کرد:
(۲-۲۵……
فهرست مراجع
مراجع انگلیسی
۱- Vahid Nejati and Kazuo Matsuuchi, Aerodynamics Design and Genetic Algorithms for Optimization of Airship Bodies, JSME, No. 02-4140, (2002).
۲- Parsons, J.S. and Goodson R.E, Shaping of Axisymmetric Bodies for Minimum Drag in Incomperessible Flow J. Hydronautics, Vol. 8, No. 3 (1974).
۳- Zedan, M. F., Potential Flow Around AxisymMetric Bodies, Direct and Inverse Problem, Ph.D. Dissertation, University of Houston, (1979).
۴- Pinebrook, W. E., Drag Minimization on a Body of Revolution, Dissertation in the University of Houston, (1982).
۵- Young, A. D., the Calculation of Total and Skin Friction Drags of Bodies of Revolution at Zero Iincidence ARC R & M, No. 1874 (1939).
۶- Rechenberg, I., Evolution Strategie: Optimize-rung Technischer Systeme Nach Prinzipien der Biologischen Evolution, (1973), Frommann-holz-boog verlag, Stuttgart.
۷- Holland, J., Adaptation in Natural and Artificial System, (1975), University of Michigan Press annarbor.
۸- Cebeci, T. and Bradshaw, P., Momentum Transfer in Boundary Layers, McGraw- Hill, (1977).
۹- Nash, J.F., Turbulent Boundary Layer Behavior and the Auxiliary Equation, ARC CP 835, London (1965).
۱۰- Shanebrook, J.R. and Sumner, W.J., Entrainment Theory for Axisymmetric Turbulent Incompressible Boundary Layer, J. Hydronautics, Vol. 4, No. 4 (1970).
۱۱- Standen N.M., A Concept of Mass Entrainment Applied to Compressible Turbulent Boundary Layers in Adverse Pressure Gradients, Proceedings on the 4th Congress of ICAS, pp. 1101-1125 (1965).
۱۲- Schlichting, H., Boundary Layer Theory, McGraw –Hill Book Co., N.Y. (1968).
۱۳- Lutz, Th.and Wagner, S., Drag Reduction and shape Optimization Air ship Bodies, J. Aircraft, vol.35, No3.(1998) , pp. 345 – 351
مراجع فارسی
۱۳- Mathews, John, H, Numerical methods for methmatics science and engineering (1943).
۱۴- Smith, I. M, Programming in Fortran 90 for engineers and scientists (1995).

برای خرید اطلاعات خود را وارد کنید
  • کلیه پرداخت های سایت از طریق درگاه بانک سامان انجام می گیرد.هر مرحله از خرید می توانید مشکل خود را با پشتیبان و فرم تماس با ما در جریان بگذارید در سریعترین زمان ممکن مشکل برطرف خواهد شد
  • پس از پرداخت وجه ، فایل محصول هم قابل دانلود می باشد و هم به ایمیل شما ارسال می گردد .
  • آدرس ایمیل را بدون www وارد نمایید و در صورت نداشتن ایمیل فایل به تلگرام شما ارسال خواهد شد .
  • در صورت داشتن هرگونه سوال و مشکل در پروسه خرید می توانید با پشتیبانی سایت تماس بگیرید.
  • پشتیبان سایت با شماره 09383646575 در هر لحظه همراه و پاسخگوی شماست
  • اشتراک گذاری مطلب

    راهنما

    » فراموش نکنید! بخش پشتیبانی مقاله آنلاین ، در همه ساعات همراه شماست

    اطلاعات ارتباطی ما پست الکترونیکی: Article.university@gmail.com

    تماس با پشتیبانی+ ایدی تلگرام 09383646575

    برای سفارشتان از سایت ما کمال تشکر را داریم.

    از اینکه ما را انتخاب نمودید متشکریم.

    معادله فوق را حل نمایید *

    تمام حقوق مادی , معنوی , مطالب و طرح قالب برای این سایت محفوظ است