ارتعاشات آزاد پوسته های استوانه ای ترکدار

با سلام و احترام. من دانشجوی رشته مهندسی مکانیک (مقطع دکتری) هستم و در زمینه ارتعاشات پوسته ها فعالیت میکنم. پروژه من تحلیل ارتعاشات آزاد پوسته ها و پنلهای استوانه ای ترکدار به کمک نرم افزار متلب است. بخشی از پروژه را انجام داده ام و ادامه آن را میخواهم برون سپاری کنم. این پروژه تنها به افرادی که به مباحث ارتعاشات آشنایی دارند سپرده میشود. همچنین پیشنهادهای بالای 20 میلیون تومان پذیرفته نمیشود بزرگوارانی که علاقمند به انجام این پروژه هستند، بایستی نکات زیر را در انجام کار در نظر بگیرند: جنس پنلها و پوسته ها از 4 نوع ماده: ایزوتروپ، ارتوتروپ تک لایه، مواد مدرج تابعی و پلیمرهای تقویت شده با نانولوله های کربنی تشکیل شده است از نظریه تغییر شکل برشی مرتبه اول و روش حل ریتز برای تحلیل ارتعاشات استفاده میشود. روند کلی انجام این پروژه به این صورت است که تحلیل ارتعاشات پوسته ها و پنلهای ترکدار همانند ارتعاشات پوسته ها و پنلهای بدون ترک است. در ارتعاشات پوسته ها و پنلهای بدون ترک، پاسخ معادلات بصورت یک سری (سری چبیشف) در نظر گرفته میشود. حال اگر سازه دچار ترک باشد جواب معادلات از جمع دو سری (یک سری برای پوسته بدون ترک (که همان سری چبیشف است) و یک سری دیگر که برای ترک است) در نظر گرفته میشود. تحلیل ارتعاشات آزاد برای پنلها و پوسته های استوانه ای بدون ترک برای چهار نوع ماده ای که ذکر شده، انجام شده است و صحت سنجی نیز شده است. شما باید جواب سری برای ترک را به معادلات پنل و پوسته استوانه ای اضافه کنید شود و فرکانسها و شکل مودهای سازه ترکدار را محاسبه کنید کدهای و مقالاتی که نیاز است در پیوست ارسال شده است. نکته بسیار بسیار مهم این است که ساختار و شکل کلی کدها به هیچ وجه نباید عوض شود یعنی تعداد پوشه های فرعی و تعداد فایلهای موجود در هر پوشه نباید دست بخورد) و براساس روندی که در ادامه ذکر میکنم، ادامه کدها در جاهای مشخص شده قرار گیرند و ساختار کدها ثابت بمانند.  در فایل پیوست، دو پوشه وجود دارد. یکی به نام منابع که در آن مقالاتی که نیاز است تا کسی برای انجام پروژه بداند ارسال شده است. پوشه دیگر بنام کدها نامگذاری شده است. پوشه کدها، شامل یک کد اصلی به نام Main است که کار اصلی پروژه را انجام میدهد و در حین اجرای آن چندین فایل دیگر که در پوشه های فرعی قرار دارند فراخوانی میشوند. پوشه های فرعی سه پوشه هستند که یکی از پوشه ها به نام symbolic  کار محاسبات نمادین را انجام میدهد. پوشه specifications a boundaries برای مشخص کردن جنس و ابعاد پوسته و پنل و محاسبه ضرایب مربوط به شرایط مرزی و مولفه های ماتریس سفتی است.  پوشه سوم به نام matrices  شامل ماتریسهای جرم و سفتی است. در نهایت نتایج ماتریسهای جرم و سفتی به فایل اصلی برمیگردد و فرکانسهای اصلی محاسبه میشود. از فرمولهای موجود در مقاله 1 (موجود در پوشه منابع) برای کد نویسی پنل ها و پوسته های استوانه ای بدون ترک استفاده شده است. فرمولهای پنلها و پوسته ها مثل هم هست فقط در جواب بخاطر اینکه پنل 4 تکیه گاه و پوسته دو تکیه گاه دارد، مقداری جوابهای پوسته با پنل فرق دارد. در ادامه کار باید سراغ جوابهای ترک برویم. برای اینکه متوجه شوید جواب سری مربوط به ترک چگونه است، مقاله 3 و 4 در منابع قرار داده شده است تا عزیزان با روند کار آشنا شوند. این مقالات برای ورق است ولی روند انجام کار برای پوسته و پنل هم دقیقا همین طور است. ما میخواهیم هم ترک داخلی (ترکی که با مرز پوسته برخورد ندارد) و هم ترک کناره (ترکی که با مرز برخورد دارد) را در نظر بگیریم و طول و زاویه قرارگیری ترک دلخواه باشد (یعنی برای همه زوایا قابل استفاده باشد). برای راحتی کار شما دو عکس با نام عکس شماره 5 و 6 در پیوست قرار دارد که در آن براساس نتایج کارهای آقای هوآنگ (نویسنده مقاله شماره 3 و 4) برای ورقهای ترکدار با استفاده از نظریه برشی مرتبه اول، سریهایی برای 5 مولفه جابجایی برای ترک داخلی (در عکس شماره 5) و برای ترک کناره (در عکس شماره 6) آورده شده است. این سریها را آقای هوآنگ برای ورقهای ترکدار پیشنهاد داده است. سری هایی که برای پوسته نیز باید استفاده شود، شکلی تقریبا مشابه سریهای آقای هوآنگ خواهند داشت. سری هایی که برای بخش ترک در این پروژه به کار گرفته میشود باید از مقاله شماره 7، بخش 5، فرمولهای 23 تا 28 استخراج شود. در این فرمولها تنشها در اطراف ترک مشخص شده است که پاسخ سری مربوط به ترک با استفاده از این روابط و استفاده از فرمولهای 11، 12، 18 و 19 موجود در مقاله 14، محاسبه میشود.  برای انجام محاسبات نمادین پوسته و پنلهای ترکدار، باید دو فایل بنامهای crackedpanel و crackedshell در پوشه symbolic درست کنید و محاسبات را در آنجا انجام بدهید و نتایج مولفه های جرم و سفتی باید مشابه با سازه بدون ترک در پوشه ماتریسها ذخیره شود. همچنین محاسبات ماتریسهای جرم و سفتی بدست آمده از محاسبات نمادین را فایلهایی در پوشه matrices و بنامهای kmatrixcpanel, mmatrixcshell,mmatrixcpanel و kmatrixcshell قرار دهید. (حرف c قبل از کلمه پنل یا پوسته بمعنی ترک دار و حرف p به معنای سالم و بدون ترک است. طبق همین قاعده فایلها را نامگذاری کنید.) در کدهای موجود در فایل پیوست، مواردی مربوط به پیش تنشهای حرارتی و مکانیکی مشاهده میشود که آنها را خودم انجام میدهم و نیاز به انجام آنها توسط شما نیست. برای اینکه این اعداد در محاسبات شما دخالت نکند، مقدار آنها را در کد فایل پیوست، صفر قرار داده ام. نتایج باید با مقالات زیر صحت سنجی شود: پوسته ترکدار ایزوتروپیک با مقاله 9 و پوسته ترکدار ارتوتروپ با مقاله 10 پنل ایزوتروپ و مدرج تابعی دارای ترک داخلی با مقاله 4 و پنل ایزوتروپ و مدرج تابعی دارای ترک کناره با مقاله 3 پوسته ها و پنلهای  مدرج تابعی ترکدار با مقاله 11. در انتهای کد باید شکل مودها نیز رسم شود(بسته به نوع پنل یا پوسته و ترکدار یا بی ترک). شکل مودهای سازه ترکدار و بدون ترک باید شبیه به شکل مودهای رسم شده در شکل 13 از مقاله شماره 11 رسم شود. در طول انجام کار توسط مجری، هر کمک یا راهنمایی ایشان لازم داشته باشند، در خدمت ایشان خواهم بود