در آخرین آموزش از مجموعه مباحث اجزای محدود به روش مستقیم استخراج ماتریس سختی برای المان تیر صفحهای پرداختیم (اینجا: المان تیر صفحهای: روابط حاکم بر آن و محاسبه ماتریس سختی – روش اول). با توجه به اهمیت موضوع و در راستای تکمیل مباحث مربوط به این المان، امروز قصد داریم روش دوم برای استخراج ماتریس سختی تیر (Beam Element) را به شما عزیزان معرفی کنیم. پیشنهاد میکنم این آموزش اجزا محدود را از دست ندهید.
روش بدست آوردن ماتریس سختی المان تیر
روش دوم استخراج ماتریس سختی ، بر معادلات انرژی و استفاده از توابع شکل استوار است. اگر با بحث توابع شکل در اجزای محدود آشنا نیستید توصیه میکنم قبل از ادامه بحث، آموزش “المان و تابع شکل در اجزاء محدود” را مرور نمایید. مطابق شکل زیر، یک المان تیر صفحهای با خیز، ممان، نیروی برشی و جابجایی گرهای مشخص شده را در نظر بگیرید.
چنانچه مباحث اجزای محدود را دنبال کرده باشید میدانید که میزان انرژی ذخیره شده در تیر برابر است با:
با بازآرایی معادله فوق و استفاده از روابط اجزای محدود داریم:
لذا، ماتریس سختی به کمک این روش، از معادله زیر استخراج خواهد شد:
از طرفی، توابع شکل برای المان تیر به صورت زیر معرفی میشوند :
همانطور که مشاهده میکنید توابع شکل فوق در گره مربوط به خود دارای مقدار یک (1) بوده و در سایر گرهها صفرند و رابطه زیر بین آنها برقرار است:
همچنین، تغییر مکان در المان موردنظر بر حسب تغییر مکان گرهها از رابطه زیر قابل استخراج خواهد بود:
آخرین قطعه از پازل استخراج ماتریس سختی در المان موردنظر به کمک یافتن میزان انحنای تیر تکمیل خواهد شد:
بنابراین، ماتریس کرنش – تغییر مکان (B) در تیر برابر است با:
در نتیجه، برای المان تیر داریم:
با قرار دادن ماتریس B در معادله فوق و انتگرالگیری از آن، ماتریس سختی برای المان تیر (Beam Element) به صورت زیر حاصل خواهد شد:
روش بدست آوردن ماتریس سختی قاب (Frame)
تذکر : اتصالات در یک قاب (Frame) به صورت صلب به یکدیگر متصل شده‌اند. در نتیجه، حرکت انتقالی و چرخشی برای اعضای قاب در محل اتصال یکسان خواهد بود. از این‌رو، در محل اتصال، نیروهای محوری، برشی و ممان خمشی ایجاد میشود و لازم است برای تحلیل و مدلسازی قاب از ترکیب تیر ساده و المان میله (المان تیر کلی) استفاده گردد. چنانچه سختی محوری (المان میله) را نیز با معادلات فوق ترکیب کنیم ماتریس سختی برای المان تیر صفحهای در حالت کلی بهدست میآید:
با سلام خدمت مهندس سروری عزیز.
می خواستم بدونم نحوه استخراج مقادیر واقعی ماتریس سفتی کلی برای یک سازه در حالت کلی از آباکوس چگونه است؟ داخل فایل های حاصل از تحلیل که خود نرم افزار ایجاد کرده را جستجو کردم اما چیز خاصی ندیدم.
آیا مرجعی در این خصوص وجود داره؟
ممنون.
سلام
منظورتون را واضحتر بیان کنید
سلام ببخشید یه سوال بی ربط داشتم اما مهم.ایا سختی محوری زمانی که از تغییرات محوری صرف تظر نشود برای هر تیر با هر شرایط تکیه گاهیEA/lهست یا بسته به گیرداری و مفصل فرق دارد؟
سلام
ببینید یک میله با سطح مقطع یکنواخت که تحت نیروی محوری قرار میگیره یک تغییر مکان داره؛ اگر این مسئله را با یک مسئله اعمال نیرو به فنر معادلسازی کنیم در اصطلاح گفته میشه سختی فنرِ معادل یا به نحوی سختی محوری میله برابر AE/L خواهد بود؛ با این ترفند میشه خیلی از مسایل را حل کرد
حس میکنم این موضوع را شما با بحث تیرها و شرایط تکیه گاهی اونها اشتباه گرفتهاید؛ توصیه میکنم به کتابهای مقاومت مصالح رجوع کنید
موفق باشید . . .
سلام جناب مهندس.ممنون از توضیحتون بله دقیقا منظورم شرایط تکیه گاهی هستش.یعنی این میله تکیه گاهش هر جور که باشه مفصلی گیردار.یه رابطه بیشتر برای سختی محوریش نیست؟
عرض کردم خدمتتون
چیزی که در ذهن شماست متفاوته
معادل سازی و سختی محوری که اشاره کردید برای یک میله است که تحت بار محوری قرار گرفته
بحثی که در ذهن شماست شرایط تکیه گاهی تیر و تأثیرش در میزان خیز و جابجایی هست؛ با این موضوع تفاوت داره
این منبع را مطالعه کنید تا موضوع را درک کنید: تغییر طول میله تحت بار محوری
با سلام اگر ماتریس سختی را برای نیروی محوری بخواهیم چه باید بکنیم
سلام
مگه نیرو ماتریس سختی داره؟!!!
سلام، لطفا روش بدست آوردن ماتریس سختی در آباکوس رو هم اموزش بدهید. ممنون
سلام
روشهای زیادی وجود داره؛ میتونید مختصات گره ها، شماره المانها و ماتریس connectivity را از فایل inp استخراج و به کمک یک برنامه Matlab ماتریس سختی کل سازه را استخراج کنید
سلام
با تشکر از سایت خوبتون
من یک مخزن رو در اباکوس مدل کردم والان میخوام ماتریس سختی مخزن از اباکوس استخراج کنم . اما موفق به انجام این کار نشدم . خواستم بدونم روشی برای بدست اوردن ماتریس سختی از اباکوس وجود دارد ؟؟
و اینکه دستور *Matrix Generate ماتریس سختی در اختیار من قرار می دهد ؟؟
ممنون
سلام
میتونید از کد Matlab هم استفاده کنید؛ دستورات پایتون هم ممکنه کارگشا باشه
در خصوص دستور ذکر شده هم میتونید سرچ کنید
سلام
به چه شکل میتونم بعد از انجام تحلیل مودال جرم و سختی هر مود ارتعاشی رو از اباکوس بدست بیارم؟؟ (منظورم ماتریس سختی و جرم گلوبال نیست)
با تشکر از شما و سایت پر محتواتان
سلام
بررسی نداشته ام
سلام یه سوال داشتم که به این صفحه ربطی نداره.
آیا مرجعی رو میشناسید که بتونم ماتریس سفتی یک تیر خمیده رو پیدا کنم.
سلام
خیر متأسفانه
سلام مهندس
سوالی داشتم
در اینجا پاسخ دهید یا ایمیلی ممنون میشوم
ماتریس سفقی برای یک تیر نانو را چگونه حساب کنم
با یک تیر معمولی و یا
دو بعدی
چه فرقی دارد ؟
سپاس
سلام
اطلاعی ندارم
سلام
لطفا اگر میتوانید نحوه بدست اوردن ماتریس سختی در ایاکوس رو کامل توضیح دهید. و اینکه میشه تو اباکوس ماتریس سختی برای یک مسیر حرکتی دلخواه رو بدست اورد؟
سلام
ماتریس سختی روی یک مسیر ؟؟!!
ماتریس سختی برای المانهای مختلف فرمول متفاوتی داره؛ میتونید در کتابهای اجزای محدود ه شکل مبسوط این موارد را پیدا کنید
با سلام یک سوال در رابطه با مسایل تحلیلی قاب داشتم. اگر ما یک قاب داشته باشیم و ان را تحت دماهای مختلف قرار دهیم ایا روابط تحلیلی تحت دما برای این قاب میتوانیم پیدا کنیم یعنی اینکه طول و ضخامت و زاویه قاب تحت دما چه تغییری داشته و روابط تحلیلی برای طول و ضخامت و زاویه تحت دما میتوان بدست اورد؟! ممنون میشم راهنمایی کنید
سلام
روابط تحلیلی را نمیتونید به این طریق بدست بیارید؛ باید مستقیم از تئوریها استفاده کنید
با این کار عملاً شما یک درونیابی ساده در محدوده تغییرات اعمال شده خواهید داشت
سلام وقتتون بخیر
اگه بخوایم المان تیر رو به چند تا المان به روش اجزا محدود تنبدیل کنیم مثلا 4 المان.اسمبل کردن ماتریس سختی تیر به چه شکل هست؟
سلام
لینک زیر را مطالعه بفرمایید
https://umec.ir/stiffness-matrix-assembly/
سلام،
در چه صورت ماتریس سختی غیرمتقارن میشود؟ در اینصورت چه مشکلی ایجاد میکند؟
سلام
یکی از مواردی که ماتریس سختی نامتقارن میشه حالتیه که stiffness سازه به شدت nonlocal باشه (یعنی توزیع به شدت غیر یکنواختی داشته باشه) یا اینکه کلا مسئله به شکل nonlocal باشه
طبق تئوری nonlocal تنش در هر نقطه تنها به کرنش همون نقطه وابسته نیست و تابعی از مجموعه نقاط همسایگی اون نقطه هم خواهد بود
از نقطه نظر تئوری این nonlocal بودن یکی از منابع اصلی عدم تقارن ماتریس سختی به حساب میاد
عدم تقارن در ماتریس سختی به شدت هزینه حل معادلات را بالا میبره
سلام تیر اولر برنولی علاوه بر بار گسترده قایم ،بار گسترده محوری وارد شده تو این تیر چطور باید ماتریس سختی رو بدست بیاریم اخه من ندیدم تیر اولر برنولی بار گسترده محوری بهش وارد شده باشه و تیر چهار درجه ازادی داره دیگه در راستای قایم و خمش در ابتدا و انتها
ممنون میشم کمکم کنین
سلام
ماتریس سختی ربطی به بار نداره
ماتریس سختی تیر اویلر-برنولی هم از رابطه زیر استخراج میشه:
https://link.springer.com/content/pdf/bbm%3A978-3-319-56493-7%2F1.pdf
سلام ببخشید من حدود چند ماه هست که دنبال روش بدست اوردن ماتریس سختی تیر با تغییر مساحت خطی (تیر تیپرد)و ماتریس سختی تیر بر بستر الاستیک به روش نیوتن(روش ماتریسی) هستم ولی به هر کتاب یا جزوه ای مراجعه کردم چیز پیدا نکردم ایا امکانش هست راهنمایی کنید یا کتابی معرفی کنید.ممنون
سلام
https://www.researchgate.net/profile/Yasha-Zeinali/publication/292360921_General_Form_of_the_Stiffness_Matrix_of_a_Tapered_Beam-column/links/56ad880e08ae43a3980c8a03/General-Form-of-the-Stiffness-Matrix-of-a-Tapered-Beam-column.pdf
سلام و عرض ادب، یک سوالی داشتم از شما این که ماتریسش سختی ما متقارن شده دلیل خاصی داره و اگر دلیل خاصی داره دلیل این متقارن بودن چیست؟
سلام
اگر به نحوه اسمبل کردن ماتریس دقت کنید دلیل تقارن رو متوجه خواهید شد
درایه k_ij حاصل اثر المان i بر المان j میشه ، طبیعتاً این مقدار برابر با اثر المان j بر المان i خواهد بود که درایه k_ji رو تشکیل میده
سلام آقای مهندس. ممنون بابت تمام زحماتتون. برای المان Beam با توجه به شرط پیوستگی شیب در گره ها از تابع شکل درجه سه باید استفاده بشه، اما وقتی می خوایم المان Beam انتخاب کنیم در قسمت geometric order فقط نوع linear و quadratic وجود داره که به ترتیب درجه یک و درجه دو هستند. خواستم بپرسم آیا این اردرها ربطی به درجه تابع شکل المان ندارند و اگر ربط ندارند بیانگر چه چیزی هستند؟
سلام
این اردرها همان مرتبه تابع شکل هستند
اگر قصد استفاده از المان خاصی دارید که در نرم افزار وجود نداره باید سابروتین بنویسید