تحلیل تیر دارای مفصل داخلی به روش المان محدود با متلب
در این پست پروژه تحلیل تیر دارای مفصل داخلی به روش المان محدود با متلب را قرار دادیم که با تشکیل ماتریس سختی تیر دو بعدی فوق را آنالیز می نماید. کد نویسی این پروژه از صفحه ۸۷ تا ۹۰ کتابی با عنوان مقدمه ای بر اجزا محدود با بکارگیری برمامه MATLAB و ABAQUS نوشته شده که این کتاب از لینک های مربوطه قابل دانلود است.
مشخصات تیر دو بعدی
- در اندازه گیری طول از واحد متر و نیرو از واحد کیلونیوتن استفاده شده است.
- ممان اینرسی تیر دو بعدی برای کل المان ۰٫۰۰۰۶ متر به توان ۴ در نظر گرفته شده است.
- مدول الاستیک تیر برای کل المان با علامت E برابر ۸^۱۰*۲ کیلوپاسکال در نظر گرفته شده است.
- در فاصله ۶ متری از تکیه گاه سمت راست یک مفصل داخلی در تیر تعبیه گردیده است.
موارد مربوط به شرایط مرزی تیر
هر یک از گره ها دارای دو درجه آزادی به صورت جابجایی و دوران هستند که جابجایی عمودی در امتداد محور قائم ( y ) و دوران در حول محور z که عمود بر محور xy است می باشد.
این تیر سه المان و چهار گره دارد.
اتصالات تکیه گاه تیر به صورت گیردار است که با کد ۱ و ۴ قابل شناسایی هستند، در کد نویسی برای مشخص کردن غیر قابل جابجایی و دوران گره ها با عدد صفر کد گذاری شده اند که از قسمت Nodal freedom قابل مشاهده است.
فایل کدنویسی شده پروژه
در محیط نرم افزار متلب ۸ فایل مجزا کدنویسی شده است که فایل اصلی آن با پسوند beam.m می باشد که با اجرای آن در کنار سایر فایل ها خروجی در نوت پد با پسوند txt به صورت زیر ایجاد می شود.
خروجی پروژه
******* PRINTING MODEL DATA **************
——————————————————
Number of nodes: ۴
Number of elements: ۳
Number of nodes per element: ۲
Number of degrees of freedom per node: ۲
Number of degrees of freedom per element: ۴
——————————————————
Node X
۱, ۰
۲, ۴
۳, ۸
۴, ۱۴
——————————————————
Element Node_1 Node_2
۱, ۱, ۲
۲, ۲, ۳
۳, ۳, ۴
——————————————————
Element E I
۱, ۲e+008, ۰٫۰۰۰۶
۲, ۲e+008, ۰٫۰۰۰۶
۳, ۲e+008, ۰٫۰۰۰۶
——————————————————
————-Nodal freedom—————————-
Node disp_w Rotation
۱, ۰, ۰
۲, ۱, ۲
۳, ۳, ۴
۴, ۰, ۰
——————————————————
—————–Applied Nodal Loads——————-
Node load_Y Moment
۱, ۰ ۰
۲, -۱۲ ۰
۳, ۰ ۰
۴, ۰ ۰
——————————————————
Total number of active degrees of freedom, n = 4
——————————————————–
******* PRINTING ANALYSIS RESULTS **************
——————————————————
Global force vector F
-۱۲
۰ ۰ ۰
——————————————————
Displacement solution vector: delta
-۰٫۰۰۰۹۶
-۰٫۰۰۰۲۷
-۰٫۰۰۱۵۸
۰٫۰۰۰۴۰
——————————————————
Nodal displacements
Node disp_y rotation
۱, ۰٫۰۰۰۰۰, ۰٫۰۰۰۰۰
۲, -۰٫۰۰۰۹۶, -۰٫۰۰۰۲۷
۳, -۰٫۰۰۱۵۸, ۰٫۰۰۰۴۰
۴, ۰٫۰۰۰۰۰, ۰٫۰۰۰۰۰
——————————————————
Members actions
element fy1 M1 Fy2 M2
۱, ۹٫۳۶, ۲۶٫۹۰, -۹٫۳۶, ۱۰٫۵۵
۲, -۲٫۶۴, -۱۰٫۵۵, ۲٫۶۴, ۰٫۰۰
۳, -۲٫۶۴, -۰٫۰۰, ۲٫۶۴, -۱۵٫۸۲
مشاهده ویدئو در این باره