تحلیل قاب با سقف شیبدار به روش المان محدود با نرم افزار MATLAB
در این پست پروژه تحلیل قاب با سقف شیبدار به روش المان محدود با نرم افزار MATLAB را قرار دادیم، کد نویسی این پروژه از ابتدای فصل ۴ صفحه ۱۰۷ تا ۱۱۸ کتاب به زبان انگلیسی با موضوع مقدمه ای بر المان محدود با استفاده از نرم افزار MATLAB و ABAQUS نوشته شده است که این کتاب را می توانید با کلیک بر روی لینک موجود در ادامه موضوع کتاب دانلود نمایید.
توضیحات پروژه
کدنویسی انجام شده در محیط نرم افزار متلب از ترکیب کدنویسی برای تحلیل خرپا و تیر می باشد که عمده تفاوت آن در ابعاد ماتریس ها است، همانند ساختار کدنویسی خرپا ها می بایست در قاب ها هم ماتریس های سختی تک تک اعضا با مختصات محلی نوشته شود و سپس سختی کل قاب با مختصات کلی بدست آید.
واحدهای اندازه گیری و توضیحات
برای اندازه گیری طول از میلی متر و نیرو از نیوتن استفاده شده است. یک المان تیر-ستون دارای سه درجه آزادی است، یکی به صورت افقی در امتداد محور X ، یکی به صورت عمودی در امتداد محور Y و در آخر چرخش حول محور Z که عمود بر محور XY است. مطابق شکل ارائه شده گره های شماره ۱ و ۵ کاملاً ثابت هستند، یعنی اتصالات پای ستون به صورت گیردار در نظر گرفته شده است، با این وجود در گره های شماره ۱ و ۵ به دلیل ثابت بودن با عدد صفر نشان داده می شوند، گره های شماره ۲ ، ۳ و ۴ آزاد هستند که درجه آزادی آنها با نماد u ، v و θ نشان داده و با عدد ۱ مشخص شده اند.
مشخصات قاب دو بعدی شیبدار
- سطح مقطع اعضا برابر ۵۲۱۰ میلی متر مربع در نظر گرفته شده است.
- ممان اینرسی ۷^۱۰*۸٫۶۴ میلی متر به توان ۴ در نظر گرفته شده است.
- مدول الاستیک برابر ۲۰۰۰۰۰ مگاپاسکال می باشد.
- دهانه قاب ۱۲۰۰۰ میلی متر است، ارتفاع راس سقف در وسط دهانه ۶۰۰۰ میلی متر و نیز ارتفاع ستون برابر ۵۰۰۰ میلی متر می باشد.
فایل کدنویسی شده پروژه
در محیط نرم افزار متلب ۹ فایل جداگانه کدنویسی شده که فایل اجرایی آن با پسوند frame.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 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 | File:frame_problem1_results.txt ******* PRINTING MODEL DATA ************** ------------------------------------------------ Number of nodes: 5 Number of elements: 4 Number of nodes per element: 2 Number of degrees of freedom per node: 3 Number of degrees of freedom per element: 6 ------------------------------------------------ Node X Y 1, 0000.00 0000.00 2, 0000.00 5000.00 3, 6000.00 6000.00 4, 12000.00 5000.00 5, 12000.00 0000.00 ------------------------------------------------ Element Node_1 Node_2 1, 1, 2 2, 2, 3 3, 3, 4 4, 4, 5 ------------------------------------------------ Element E A I 1, 200000, 5210 8.64e+007 2, 200000, 5210 8.64e+007 3, 200000, 5210 8.64e+007 4, 200000, 5210 8.64e+007 ------------------------------------------------ -------------Nodal freedom---------------------- Node disp_u disp_u Rotation 1, 0, 0, 0 2, 1, 2, 3 3, 4, 5, 6 4, 7, 8, 9 5, 0, 0, 0 ------------------------------------------------ ------------------------------------------------ -----------------Applied joint Loads------------ Node load_X load_Y Moment 1, 0000.00, 0000.00, 0000.00 2, -5999.99, 35999.93, 37000000.00 3, -11999.98, 0000.00, -74000000.00 4, -5999.99, -35999.93, 37000000.00 5, 0000.00, 0000.00, 0000.00 -------------------------------------------------- Total number of active degrees of freedom, n = 9 -------------------------------------------------- ******* PRINTING ANALYSIS RESULTS ************ -------------------------------------------------- Global force vector F -5999.99 35999.9 3.7e+007 -12000 0 -7.4e+007 -5999.99 -35999.9 3.7e+007 ----------------------------------------------------- Displacement solution vector: delta -25.03159 0.16363 0.00712 -25.04119 0.00000 -0.00686 -25.03159 -0.16363 0.00712 ------------------------------------------------------ Nodal displacements Node disp_x disp_y rotation 1, 0.00000e+000, 0.00000e+000, 0.00000e+000 2, -2.50316e+001, 1.63630e-001, 7.11912e-003 3, -2.50412e+001, 7.98515e-015, -6.85508e-003 4, -2.50316e+001, -1.63630e-001, 7.11912e-003 5, 0.00000e+000, 0.00000e+000, 0.00000e+000 ------------------------------------------------------ Members actions in local coordinates element fx1 fy1 M1 fx2 Fy2 M2 1, -34100.6 -11999.976 -54603626.5 34100.6 11999.97 -5396249.9 2, 6230.6 -35609.4 5396249.9 -6230.6 -37383.64 -0.0 3, -6230.6 37383.6 0.0 6230.6 35609.4 5396249.9 4, 34100.6 -11999.97 -5396249.9 -34100.6 11999.97 -54603626.5 ------------------------------------------------------ Members actions in global coordinates element fx1 fy1 M1 fx2 Fy2 M2 1, 11999.9 -34100.5 -54603626.5 -11999.9 34100.5529, -5396249.9 2, 11999.9 -34100.5 5396249.9 0.0 -37899.6 -0.0 3, 0.0 37899.3 0.0000, 11999.9 34100.5529, 5396249.9 4, -11999.9 -34100.5 -5396249.9 11999.9 34100.6 -54603626.5 |
هیچ نظری ثبت نشده است