دانشنامه تخصصی مهندسی ایران

دانشنامه تخصصی مهندسی ایران

 engpedia   راه اندازی کانال تلگرام ( EngPedia_ir@ )

                همراه با مطالب ویژه سایت بصورت رایگان

  • نسخه ۱۲.۱ نرم افزار Aveva Marine در سایت قرار گرفت. (اینجا)
  • نسخه ۱۰.۵۲ نرم افزار Plus 2D در سایت قرار گرفت. (اینجا)
  • نسخه ۱۱.۰ نرم افزار AGI Systems Tool Kit-STK در سایت قرار گرفت.(اینجا)

همکاران

حل یک مساله تیر (Beam) با روشهای FEM

مدیریت فروردین ۲, ۱۳۹۳ 4496 بازدید ۰دیدگاه

در این پسن قصد داریم یک مساله را با چند روش اجزا محدود (FEM) حل نماییم. در این مساله یک تیر که تحت نیروی گسترده و ممان قرار دارد، در دو نقطه مقید شده است. مساله با روشهای WRM، گلرکین (Galerkin)، روش Collocation point و روش حداقل مربعات جزئی (Least Square) حل خواهد شد.

صورت مساله

تیر زیر را با مشخصات (L, E, I, A) در نظر بگیرید که دو ممان و نیروی گسترده به آن اعمال شده است( (f(x و M0 ). هدف از حل مساله، بدست آوردن تغییر شکل (deflection)، تنش (stress) و کرنش (strain) در تیر در سه حالت زیر می باشد.
الف) حل مساله با روش معادلات دیفرانسیلی (جواب دقیق)
ب) حل با روش Weighted Residual Methods یا همان WRM
ج) رسم نمودارهای جابجایی، تنش و کرنش

f(x) = f۰ [۱-(x/l)۰.۵] otgg

beam

حل یک مساله تیر (Beam) با روشهای FEM

 

برای نگارش حل دقیق با استفاده از دیفرانسیل مرتبه چهارم می توانیم تابع گشتاور و شیب و خیز را در طول تیر بدست بیاوریم:

beam equation

با توجه به روابط بالا و شرایط مرزی که در تیر موجود است میتوانیم دو ثابت C3 و C4 را اینگونه بدست آوریم.

با توجه به معادله شیب و خیز در نقطه ۰ = نتیجه می شود که C3 و C4 صفر بوده پس معادله ما به دو ضریب مجهول C1 و C2 تقلیل پیدا خواهد کرد. برای حل این دو مجهول دیگر با استفاده از دو شرط مرزی دیگر که در نقطه x=l تعریف می شوند میتوانیم با حل دو معادله دومجهول ضرائب C1 و C2 را بدست بیاوریم.

beam

پس از بدست آوردن این دو ضریب بر حسب پارامترهای دیگر می توانیم در نهایت یک تابع به عنوان تابع خیز در طول تیر داشته باشیم.که همان حل دقیق مسئله ماست. برای حل از نرم افزار متلب کمک گرفته و خواهیم داشت:

 ضرایب C1و C2 معادله خیز به شکل زیر است:

حال برای رسم یک نمونه نمودار تغییر مکان نیازمند تعریف یکسری مقادیر اولیه برای پارامترهای مسئله هستیم که به صورت فرضی این مقادیر را اخد می نماییم:

 

برای رسم نمودار بایستی مقادیری که مشخص کرده ایم را در معادلات جایگذاری نموده و نقاط رسم را تعیین کنیم.

نمودارهای تنش و کرنش هم با توجه به معادله گشتاور که بدست آورده ایم و رابطه σ=mc/I میتوان مستقیما نتیجه شود. و کرنش هم با توجه به σ=Ee بدست خواهد آمد در نتیجه داریم:

 نتایج حاصل از حل دقیق مسئله در نمودار های روبرو نمایش داده شده است.

beam

مقدار M یا گشتاور در انتهای تیر را بگونه ای گرفته ایم که ممان خمشی حاصل تغییر قابل ملاحظه ای در مقایسه با F کششی گسترده که بر روی تیر اعمال میشود، داشته باشد.

همچنین نمودار های تنش و کرنش را آنچنان که گفته شد، رسم کرده و از آنجا که تنش و کرنش با ضریب E به یکدیگر مربوطند در نتیجه نمودار های حاصل به یک شکل است منتهی اعداد بدست آمده متفاوت می باشد.

beam

تقریب با استفاده از روش Weighted Residuals Methods
در این قسمت ما نیازمند حدس یک تابع تقریب با ضرایب نامعلوم هستیم که بایستی در گام اول شرایط مرزی ما را برآورده کند. در کار حاضر ما اقدام به تخمین تابع خیز توسط یک چند جمله ای مرتبه ۵ می نماییم. سپس مرتبه تابع تخمین را یک مرتبه افرایش می دهیم و آنگاه نتایج را با حالت قبل بررسی می کنیم. تابع تقریب yapp به صورت درجه ۵ فرض می کنیم فرض می کنیم. اما این حدس به گونه است که شرایط مرزی را برآورده کند.

beam

تابع تقریب درجه ۶ هم به صورت زیر است:

beam

حال با توجه به دو شرط مرزی یکی (تغییر مکان صفر و شیب صفر(در ابتدای تیر و یکی) تغییر مکان صفر( در انتهای تیر،برخی از ضرایب چند جمله ای بالا بدین گونه محاسبه می شوند.

beam

شرط مرزی دیگر گشتاور Ml در انتهای تیر است که به صورت نقطه ای وارد میشود و میتوان با معادله که بدست می آید یکی دیگر از این ضرایب را معلوم کرد.

beam

از این معادله میتوان ضریب a2 را بر حسب ضرایب دیگر محاسبه کرد و آنرا از تابع تقریب حذف نمود.در نهایت سه ضریب a3 ,a4,a5 در دو تابع تقریب ما باقی خواهند ماند که با توجه به روش های متفاوت میتوانیم آنها را بدست آوریم. حال تابع باقیمانده را برای هریک از توابع تخمین به صورت جداگانه محاسبه کرده و داریم:

beam

تمامی مراحل بالا را به صورت زیر کد نویسی کرده و داریم:

 نتایج محاسبات بالا بدین ترتیب است:

 

 Collocation point Method

در این روش کافی است بدانیم تابع تخمین ما دارای چند ضریب مجهول است تا بوسیله معادلاتی که بین بازه مورد نظر می نویسیم بتوانیم آن ضرایب را معین نماییم. در حالت تابع تخمین درجه ۵، ما دو مجهول داریم a3,a4 و در حالت درجه ۶ سه مجهولa3,a4,a5 را داریم. به راحتی می توانیم با نوشتن تابع باقیمانده در نقاطی بین بازه ای معادلاتی پیدا کنیم که ضرایب مجهول ما را معلوم نمایند.معادلات از این قرارند:

beam

 نتایج تابع تخمین مرتبه ۵:

 نتایج تابع تخمین مرتبه ۶:

 نتیجه:
در نمودار زیر دو تابع تخمین به همراه تابع دقیق خیز رسم شده اند که حاکی از دقت بالای مرتبه ۶ دارد.

beam

در این روش هم با استفاده از انتگرال گیری از تابع باقیمانده در بازه های بین بازه ای، میتوانیم ضرایب را برای معادله جدید بدست بیاوریم. فرضا برای تابع تخمین درجه ۶ داریم:

beam

با استفاده از سه معادله بالا میتوانیم ضرایب a3, a4, a5 را بدست آوریم:

 نتایج محاسبات این روش به شکل زیر است:
درجه ۵ :

 درجه ۶ :

 

در نهایت توسط رسم این توابع متوجه تفاوت مرتبه ۵ و دقت مرتبه ۶ خواهیم شد.

beam

حل با استفاده از روش Galerkin’s Method

این روش هم همانند روش قبلی است منتها انتگرال های میان بازه ای به شکل زیر تغییر میکنند که Wi=Qi  می باشد. Qi عباراتی هستند که در کنار ضرایب مجهول ظاهر می شوند. که در این مسئله برابرند با:

beam

برای درجه ۵ و درجه ۶ :

 حال با توجه به انتگرال بین بازه ای بالا معادلات مختلفی را برای حل در بازه های مختلف بین بازه اصلی تعریف میکنیم و در نهایت داریم:

 

 در نهایت پس از حل معادلات و معلوم کردن ضرایب داریم:

beam

حل با استفاده از روش Least Square

در این روش هم همانند Galerkin بایستی تابع وزن در باقیمانده درون انتگرال ضرب شود منتها این تابع مشتق از خود باقیمانده نسبت به ضرایب مجهول است. به عبارتی دیگر معادلات ما توسط این انتگرال ساخته می شوند:

beam

در نتیجه خواهیم داشت:

 

 و در نهایت نتایج را رسم می نماییم:

beam

 

در نهایت در یک نگاه کلی برای تمامی تقریب های مرتبه ۵ داریم:

beam

نظرات شما باعث دلگرمی و پیشرفت ما می شود.

telegram: @EngPedia_ir

دیدگاهتان را بنویسید

نشانی ایمیل شما منتشر نخواهد شد. بخش‌های موردنیاز علامت‌گذاری شده‌اند *