انیشتن
عضو جدید
کمممممممممممممممک.سلام بچه هامن یک پروژه المان محدود دارم روش حل دارم نتایج کد متلب هم دارم فقط کدمتلب المان برا ورق می خوام .پروژه کمکی هم دارم کسی می تونه حل کنه؟سوالم اینه
b=1m
t=10mm([FONT="]ضخامت[/FONT])
file:///C:\DOCUME~1\ADMINI~1\LOCALS~1\Temp\msohtmlclip1\01\clip_image001.pngv([FONT="]پواسون[/FONT])=0.3
q=sin(πx/a)*sin(πy/b)*coswt[FONT="]و[/FONT]w=1[FONT="]
[/FONT]

[FONT="]به نام خدا[/FONT]
[FONT="]یک ورق به ابعاد[/FONT]a[FONT="] و[/FONT]b[FONT="]،تعدادالمان درجهت[/FONT]x[FONT="] والمان درجهت[/FONT]y[FONT="] با بار[/FONT]q(x,y,z,t)[FONT="] رادر نظربگیریدبرای دوحالت زیر حل کنید.[/FONT]
[FONT="]الف)همه تکیه گاه ها لولا باشند.(استاتیکی حل کنید)[/FONT]
[FONT="]ب)همه تکیه گاه ها گیردارباشد.[/FONT]
E=200*e^9
a=1m b=1m
t=10mm([FONT="]ضخامت[/FONT])

file:///C:\DOCUME~1\ADMINI~1\LOCALS~1\Temp\msohtmlclip1\01\clip_image001.pngv([FONT="]پواسون[/FONT])=0.3
q=sin(πx/a)*sin(πy/b)*coswt[FONT="]و[/FONT]w=1[FONT="]
[/FONT]
[FONT="]این شکل بالا همان ورق هست بابارخمشی [/FONT]q[FONT="].مرسی[/FONT]
[FONT="]من کل کارهایی که باید انجام بدید[/FONT][FONT="]می گم[/FONT][FONT="].
1. [/FONT][FONT="]چون مساله ورق هستش و ما نیرو خمشی داریم پس باید از المان[/FONT][FONT="] plate [/FONT][FONT="]برای حل استفاده کنیم. المان مثلثی یا چهارگوش چون نمی توانند نیروی در راستای عمود بر صفحه تحمل کنند، برای این مساله قابل استفاده نیس.[/FONT][FONT="]
2. [/FONT][FONT="]این المان و ماتریس سختی در کتاب اجزا محدود[/FONT][FONT="] Logan [/FONT][FONT="]در فصل 12 اومده.[/FONT][FONT="]
3. [/FONT][FONT="]ابتدا باید کد المان ورق برای یک تک المان بنویسیم که ورودی های مساله مختصات نودها باید باشه و خروجی ماتریس سختی که چون المان ورق دارای چهار نود و 3 درجه آزادی (جابجایی در راستای عمود بر ورق و دو شیب در صفحه ورق)، پس باید ماتریس سختی برای تک المان 12 در 12 باشه[/FONT][FONT="]
4. [/FONT][FONT="]برای نیروی[/FONT][FONT="] q [/FONT][FONT="]هم شما باید بر اساس معادله که در همین کتاب هست.معادله (12.2.20) این نیروی در هر المان را به چهار نیرو در چهار نود المان، تقسیم کنیم. چون نیروی شما تابعی از[/FONT][FONT="] x [/FONT][FONT="]و[/FONT][FONT="] y [/FONT][FONT="]هست باید در المان انتگرال گیری بشه[/FONT][FONT="].
5. [/FONT][FONT="]بعد از به دست آوردن نیرو و ماتریس سختی در هر المان، باید همه ی ماتریس کلی ساخته بشه که بر اساس المان بندی و شماره نود انجام میشه[/FONT][FONT="].
6. [/FONT][FONT="]بعد به دست آوردن ماتریس سختی و بردار نیرو، باید شرایط مرزی اعمال شه[/FONT][FONT="].
[/FONT][FONT="]چون شرایط مرزی برای ورق هست، منظور از شرایط گفته شده یعنی چهار طرف ورق لولا یا گیر دار باشه. در این صورت برای لولا فقط جابجایی در جهت عمود بر صفحه صفره (در المان ورق فقط جابجایی در راستای عمود بر صفحه در نظر گرفته می شه). و برای گیردار، جابجایی و دو شیب (شیب در جهت[/FONT][FONT="] x [/FONT][FONT="]و[/FONT][FONT="] y) [/FONT][FONT="]باید صفر باشه[/FONT][FONT="].
9. [/FONT][FONT="] بعدش برای این دو شرایط مرزی، باید سطر و ستون مربوطه در ماتریس سختی حذف گردد تا ماتریس از[/FONT][FONT="] singularity [/FONT][FONT="]خارج شه[/FONT][FONT="].
10. [/FONT][FONT="]الان می تونیم برای حل، ماتریس سختی را معکوس و در بردار نیروی متناظر ضرب و جابجایی و شیب در هر نود را به دست بیاریم[/FONT][FONT="].
11. [/FONT][FONT="]برای قسمت دینامیکی هم فرقی با استاتیکی نداره و فقط زمان را صفر تا زمان مورد نظرتون مثلا 0.1 ثانیه تغییر داده و در هر زمان همه مراحل بالا باید انجام بشه (البته فقط بردار نیرو تغییر می کند) و جابجایی و شیب ها بدست می اد و می تونیم مثلا نمودار جابجایی بر حسب زمان را رسم کنیم.[/FONT][FONT="]
[/FONT]
1. [/FONT][FONT="]چون مساله ورق هستش و ما نیرو خمشی داریم پس باید از المان[/FONT][FONT="] plate [/FONT][FONT="]برای حل استفاده کنیم. المان مثلثی یا چهارگوش چون نمی توانند نیروی در راستای عمود بر صفحه تحمل کنند، برای این مساله قابل استفاده نیس.[/FONT][FONT="]
2. [/FONT][FONT="]این المان و ماتریس سختی در کتاب اجزا محدود[/FONT][FONT="] Logan [/FONT][FONT="]در فصل 12 اومده.[/FONT][FONT="]
3. [/FONT][FONT="]ابتدا باید کد المان ورق برای یک تک المان بنویسیم که ورودی های مساله مختصات نودها باید باشه و خروجی ماتریس سختی که چون المان ورق دارای چهار نود و 3 درجه آزادی (جابجایی در راستای عمود بر ورق و دو شیب در صفحه ورق)، پس باید ماتریس سختی برای تک المان 12 در 12 باشه[/FONT][FONT="]
4. [/FONT][FONT="]برای نیروی[/FONT][FONT="] q [/FONT][FONT="]هم شما باید بر اساس معادله که در همین کتاب هست.معادله (12.2.20) این نیروی در هر المان را به چهار نیرو در چهار نود المان، تقسیم کنیم. چون نیروی شما تابعی از[/FONT][FONT="] x [/FONT][FONT="]و[/FONT][FONT="] y [/FONT][FONT="]هست باید در المان انتگرال گیری بشه[/FONT][FONT="].
5. [/FONT][FONT="]بعد از به دست آوردن نیرو و ماتریس سختی در هر المان، باید همه ی ماتریس کلی ساخته بشه که بر اساس المان بندی و شماره نود انجام میشه[/FONT][FONT="].
6. [/FONT][FONT="]بعد به دست آوردن ماتریس سختی و بردار نیرو، باید شرایط مرزی اعمال شه[/FONT][FONT="].
[/FONT][FONT="]چون شرایط مرزی برای ورق هست، منظور از شرایط گفته شده یعنی چهار طرف ورق لولا یا گیر دار باشه. در این صورت برای لولا فقط جابجایی در جهت عمود بر صفحه صفره (در المان ورق فقط جابجایی در راستای عمود بر صفحه در نظر گرفته می شه). و برای گیردار، جابجایی و دو شیب (شیب در جهت[/FONT][FONT="] x [/FONT][FONT="]و[/FONT][FONT="] y) [/FONT][FONT="]باید صفر باشه[/FONT][FONT="].
9. [/FONT][FONT="] بعدش برای این دو شرایط مرزی، باید سطر و ستون مربوطه در ماتریس سختی حذف گردد تا ماتریس از[/FONT][FONT="] singularity [/FONT][FONT="]خارج شه[/FONT][FONT="].
10. [/FONT][FONT="]الان می تونیم برای حل، ماتریس سختی را معکوس و در بردار نیروی متناظر ضرب و جابجایی و شیب در هر نود را به دست بیاریم[/FONT][FONT="].
11. [/FONT][FONT="]برای قسمت دینامیکی هم فرقی با استاتیکی نداره و فقط زمان را صفر تا زمان مورد نظرتون مثلا 0.1 ثانیه تغییر داده و در هر زمان همه مراحل بالا باید انجام بشه (البته فقط بردار نیرو تغییر می کند) و جابجایی و شیب ها بدست می اد و می تونیم مثلا نمودار جابجایی بر حسب زمان را رسم کنیم.[/FONT][FONT="]
[/FONT]
[FONT="] [/FONT]
مرسی بهم پیام خصوصی بزنیدهرکی می تونه کمکم کنه.متشکرم