تحلیل ماتریسی تیر دو بعدی دارای مفصل داخلی با MATLAB
این پست مربوط به پروژه تحلیل ماتریسی تیر دو بعدی دارای مفصل داخلی با MATLAB می باشد که با استفاده از روش المان محدود و تشکیل ماتریس سختی تیر مربوطه را آنالیز می کند، این کد نویسی از صفحه ۸۷ تا ۹۰ کتابی به زبان انگلیسی با عنوان مقدمه ای بر اجزا محدود با استفاده از نرم افزار MATLAB و نرم افزار ABAQUS نوشته شده است که قابل دانلود می باشد.
شکل تیر
مشخصات تیر
- واحد اندازه گیری طول به صورت متر و نیرو به صورت کیلونیوتن در نظر گرفته شده است.
- ممان اینرسی تیر با نماد I به مقدار ۰٫۰۰۰۶ متر به توان ۴ و مدول الاستیسیته با علامت E به مقدار ۸^۱۰*۲ کیلوپاسکال در نظر گرفته شده است.
- در فاصله ۸ متری از تکیه گاه سمت چپ یک مفصل داخلی در تیر ایجاد شده است.
شرایط مرزی تیر دو بعدی
- هر گره دو درجه آزادی دارد که یکی از آن مربوط به تغییرمکان عمودی در امتداد محور قائم ( y ) و دیگری دوران در حول محور z که عمود بر محور xy است می باشد.
- تیر موردنظر دارای سه المان و ۴ گره است.
- تکیه گاه های این تیر به صورت گیردار می باشد. این تکیه گاه ها با کد ۱ و ۴ مشخص شده اند که برای بسته بودن گره ها از نظر جابجایی و دوران با شماره های ۰ کدگذاری شده اند.
موارد مربوط به کد نویسی
این پروژه دارای ۸ فایل کدنویسی شده در محیط نرم افزار متلب می باشد که فایل اصلی آن با پسوند beam.m است که با اجرای آن، در کنار سایر فایل ها خروجی در نوت پد با پسوند txt به صورت زیر نمایش داده می شود.
خروجی پروژه متلب
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 | PRINTING MODEL DATA ************** ------------------------------------------------------ Number of nodes: 4 Number of elements: 3 Number of nodes per element: 2 Number of degrees of freedom per node: 2 Number of degrees of freedom per element: 4 ------------------------------------------------------ Node X 1, 0 2, 4 3, 8 4, 14 ------------------------------------------------------ Element Node_1 Node_2 1, 1, 2 2, 2, 3 3, 3, 4 ------------------------------------------------------ Element E I 1, 2e+008, 0.0006 2, 2e+008, 0.0006 3, 2e+008, 0.0006 ------------------------------------------------------ -------------Nodal freedom---------------------------- Node disp_w Rotation 1, 0, 0 2, 1, 2 3, 3, 4 4, 0, 0 ------------------------------------------------------ -----------------Applied Nodal Loads------------------- Node load_Y Moment 1, 0 0 2, -12 0 3, 0 0 4, 0 0 ------------------------------------------------------ Total number of active degrees of freedom, n = 4 -------------------------------------------------------- ******* PRINTING ANALYSIS RESULTS ************** ------------------------------------------------------ Global force vector F -12 0 0 0 ------------------------------------------------------ Displacement solution vector: delta -0.00096 -0.00027 -0.00158 0.00040 ------------------------------------------------------ Nodal displacements Node disp_y rotation 1, 0.00000, 0.00000 2, -0.00096, -0.00027 3, -0.00158, 0.00040 4, 0.00000, 0.00000 ------------------------------------------------------ Members actions element fy1 M1 Fy2 M2 1, 9.36, 26.90, -9.36, 10.55 2, -2.64, -10.55, 2.64, 0.00 3, -2.64, -0.00, 2.64, -15.82 |
هیچ نظری ثبت نشده است