روش سختی در تحلیل ماتریسی سازه ها

روش سختی در تحلیل ماتریسی سازه ها

اسلاید 1: تحلیل ماتریسی سازه هاMatrix Structural Analysisکریم عابدی

اسلاید 2: فصل سوم: روش سختی در تحلیل ماتریسی سازه ها

اسلاید 3: فصل سوم- روش سختی در تحلیل ماتریسی سازه ها1- مقدمه:-گفتیم که اگر هدف اصلی تحلیل سازه، تعیین تغییر مکان های دو انتهای عنصر یا به عبارت دیگر مشخص کردن تغییر مکان های مربوط به گره های سازه باشد، در این صورت تحلیل سازه به روش تغییر مکان ها (Displacement Method) یا روش سختی (Stiffness Method) انجام می گیرد.- در روش سختی مجهولات شامل تغییر مکان های گره ها است و تعداد معادلات حاصل برابر درجه آزادی کل گره های سازه می باشد.- بنابراین در روش سختی ابتدا تغییر مکان های نقاط مشخص به طور اخص در گره های سازه تعیین می شود و سپس نیروهای داخلی محاسبه می شوند.- در روش سختی معادلاتی بین نیروها و تغییر مکان های سازه در دو سطح عنصر و کل سازه ایجاد می شوند .

اسلاید 4: - پس در یک جمع بندی روش سختی شامل مراحل عمومی زیر است: تعیین یک مجموعه از تغییر مکان های سیستم سازه ای نوشتن روابط نیرو- تغییر مکان ارضای شرط تعادل ارضای شرط سازگاری یافتن معادلات سازگاری حل معادلات و به دست آوردن تغییر مکان های سیستم سازه ای به دست آوردن نیروهای اعضاء و واکنش های تکیه گاهیدر تحلیل ماتریسی سازه ها به روش سختی در واقع معادلات مذکور در فرم ماتریسی استخراج می شوند و مبانی جبر ماتریسی به کار گرفته می شوند. این معادلات ماتریسی شامل بردار نیرو و بردار تغییر مکان و در ضمن ماتریس دیگری خواهد بود که به ماتریس سختی معروف است و بستگی به هندسه سازه، خواص هندسی و خواص مصالح اعضاء، نوع اتصالات موجود در سازه، تکیه گاه ها، نحوه اتصال اعضا و ... دارد.

اسلاید 5: 2- تعیین معادله روش سختی:- جسم تغییر شکل پذیری (Deformable body) را در نظر بگیرید که تحت اثر نیروهای Pi قرار دارد (بارگذاری از صفر شروع شده و به طور خطی به مقدار نهایی خود Pi رسیده است) (i نقاطی از سازه می باشند که نیروهای Pi بر آن نقاط وارد می شوند).- در اثر بارگذاری مذکور سازه تغییر مکان های Δi را متحمل می شود (Δi در راستای اعمال نیروهای Piمی باشند).- با توجه به فرض رفتار خطی سازه، کار انجام شده توسط نیروهای وارد بر سازه (Pi) ناشی از تغییر مکان های سازه (Δi) به صورت زیر خواهد بود (کار انجام یافته مذکور معادل انرژی تغییر شکل جسم است):

اسلاید 6: - فرض کنید که در یکی از تغییرمکان ها (مثلاً Δ1) تغییر کوچکی داده می شود، در این صورت تغییرات انرژی تغییر شکل جسم نسبت به تغییرات در Δ1 به صورت زیر درمی آید (لازم به ذکر است که سایر تغییر مکان ها ثابت نگه داشته می شوند):- اما با توجه به قضیه اول کاستیلیانو داریم:- حال اگر عمل فوق را برای تمامی تغییرمکان ها (Δi) انجام دهیم به طور کلی به رابطه زیر خواهیم رسید:

اسلاید 7: - اگر مجموعه معادلات مذکور را به فرم ماتریسی بیان کنیم خواهیم داشت:شامل تغییر مکان های نقاط گرهی (به انضمام تغییر مکان های تکیه گاه ها)نیروهای خارجی مؤثر بر سیستم (به انضمام عکس العمل ها)ماتریس مربعی- پس کل مسأله به تعیین ماتریس مربعی مذکور- یا تعیین اعضای از ماتریس مربعی- و سپس حل این معادلات ماتریسی برمی گردد.

اسلاید 8: - می توان تقارن ماتریس مذکور را نشان داد (با استفاده از قضیه اول کاستیلیانو):- بنابراین با توجه به تقارن ماتریسی مربعی می توان نوشت:

اسلاید 9: - اکنون ببینیم مفهوم اعضای ماتریس مربعی چیست؟اگر در نقطه 1، تغییر مکان کوچک Δ1 وارد شود و از تغییر مکان تمام نقاط دیگر جلوگیری به عمل آید، نیروی لازم برای ایجاد تغییر مکان کوچک Δ1 در نقطه 1 معادل و نیروی لازم برای جلوگیری از تغییر مکان نقاط دیگر به ترتیب عبارت خواهد بود:و اگر در نقطه i، تغییر مکان کوچک Δi وارد شود و از تغییر مکان تمام نقاط دیگر جلوگیری به عمل آید، نیروی لازم برای ایجاد تغییر مکان کوچک Δi در نقطه i معادل و نیروی لازم برای جلوگیری از تغییر مکان نقاط دیگر به ترتیب عبارت خواهد بود:اگر تغییر مکان Δi برابر واحد باشد، در این صورت نیروی لازم برای ایجاد تغییر مکان واحد در نقطه i خواهد بود و نیروی مورد نیاز برای جلوگیری از تغییر مکان نقطه 1، و نیروی مورد نیاز برای جلوگیری از تغییر مکان نقطه n، خواهد بود.

اسلاید 10: - پس با نمایش خواهیم داشت: نیروی تعمیم یافته مورد نیاز برای ایجاد تغییر مکان تعمیم یافته واحد در گره i، وقتی که سایر گره ها ثابت نگه داشته شوند.- درنتیجه معادله معروف و مشهور روش سختی به دست می آید:- نیروی تعمیم یافته مورد نیاز برای جلوگیری از تغییر مکان تعمیم یافته گره j، وقتی که تغییر مکان واحد به گره i اعمال می شود = نیروی تعمیم یافته مورد نیاز برای جلوگیری از تغییر مکان تعمیم یافته گره i، وقتی که تغییر مکان واحد به گره j اعمال می شود. =P بردار نیروی تعمیم یافته (Generalized force Vector) =K ماتریس سختی سازه (Stiffness matrix) =Δ بردار تغییر مکان تعمیم یافته (Generalized Displacement Vector)

اسلاید 11: - بنابراین یک روش برای تعیین ماتریس سختی سازه بدین صورت است که در هر دفعه تغییر مکان واحد برای گره ها (برحسب نوع سازه، مثلاً برای خرپای مسطح تغییر مکان واحد در جهت xها و تغییر مکان واحد در جهت yها- برای خرپای فضایی، تغییر مکان واحد در جهت xها و تغییر مکان واحد در جهت yها و تغییر مکان واحد در جهت zها- برای قاب مسطح تغییر مکان واحد در جهت xها و yها و دوران واحد در جهت zها و برای قاب فضایی تغییر مکان های واحد در جهت z ,y ,x و دوران های واحد در جهت z ,y ,x) در نظر گرفته شده و نیروی مورد نیاز برای ایجاد آن تغییر مکان و نیروهای نگه دارنده گره های دیگر در مقابل تغییر مکان مذکور محاسبه می شوند.برای خرپای مسطحبرای خرپای فضائیبرای قاب مسطح

اسلاید 12: مثال 1) مطلوب است تعیین ماتریس سختی سازه شکل زیر:هنگامی که تغییر مکان تعمیم یافته Δ2=1 را به گره 2 اعمال می کنیم، نیروی K22 برای ایجاد آن مورد نیاز است و نیروهای K21 و K23 وK24 برای جلوگیری از تغییر مکان های گره های 1 و 3 و 4 مورد نیاز می باشند و تاثیری در گره های 5، 6 و 7 ندارند؛ بعبارت دیگر , K52 K62 و K72 همگی مساوی صفرند. بنابراین به روشنی دیده می شود هنگامی که بین دو گره iو jعضوی وجود نداشته باشد Kij مساوی صفر است. ماتریس نواری قطری

اسلاید 13: - به نظر می رسد که تشکیل ماتریس سختی به این طریق دارای نکات ضعف عمده ای می باشد: الف) نیاز به محاسبات زیاد و وقت گیر (خصوصاً برای سازه هایی با درجات آزادی بالا)، ب) پیاده سازی آن در یک برنامه کامپیوتری بسیار دشوار (یا حتی غیر ممکن) است.- بنابراین باید به دنبال روش هایی برای تشکیل ماتریس سختی سازه بود که: الف) نیاز به محاسبات زیاد و وقت گیر نداشته باشد، ب) قابل پیاده سازی در یک برنامه کامپیوتری باشد، پ) به صورت ساده تر و مؤثرتر ماتریس سختی سازه را تشکیل نماید.- باتوجه به اینکه ماتریس سختی سازه (K) از ترکیب معقول و متناسب ماتریس های سختی هرکدام از اعضای سازه(k) تشکیل شده است، به نظر می رسد که بتوان به صورت ساده تر و مؤثرتر با استفاده از ماتریس های سختی هرکدام از اعضاء و انجام عملیات ماتریسی، ماتریس سختی سازه را تشکیل داد.

اسلاید 14: 3- تعیین ماتریس سختی عضو سازه(Member) :- یک عضو با دو گره (Node) مشخص می شود:- در حالت کلی (در فضای سه بعدی) هر گره عضو دارای شش درجه آزادی است. به عبارت دیگر در فضای سه بعدی فیزیکی، بردار مشخص تغییر مکان ها در یک گره دارای شش مؤلفه مستقل است، سه مؤلفه خطی و سه مؤلفه دورانی.- دستگاه مختصات کلی زیر را در نظر می گیریم:(Global Coordinate System) i j

اسلاید 15: - بنابراین در مجموع یک عضو در فضا دارای 12 درجه آزادی است. بنابراین ماتریس سختی یک عضو، ماتریسی 12×12 است. می توان ماتریس سختی یک عضو را به طور مستقیم در دستگاه مختصات کلی محاسبه نمود. طبیعی است که در این صورت ماتریس سختی عضو را با وارد کردن تغییر مکان ها (یک به یک) در امتداد هر یک از محورهای مختصات کلی سیستم و محاسبه نیروی مورد نیاز برای ایجاد آن تغییر مکان در امتداد آن محور خاص و نیروهای مورد نیاز برای جلوگیری از تغییر مکان در سایر امتدادها (در دو انتهای عضو) ایجاد می گردد. - به نظر می رسد که اولاً محاسبه ماتریس سختی عضو به طور مستقیم در دستگاه مختصات کلی سیستم هم طولانی و هم وقت گیر خواهد بود. - همچنین در این حالت نیروهای داخلی حاصل در انتهای اعضاء که در پایان محاسبات بدست می آیند در دستگاه مختصات کلی بیان شده و الزاماً نشانگر نیروی محوری، برشی و یا لنگر خمشی عضو نخواهد بود. بنابراین در انتهای عملیات برای یافتن این مولفه ها که در عمل بیشتر مورد لزوم هستند تبدیل مختصات ضروری خواهد بود. - بنابراین برای سادگی و نیز طولانی و وقت گیر نبودن محاسبات بهتر است که ماتریس سختی عضو در یک دستگاه مختصات مخصوصی که دستگاه مختصات محلی نامیده می شود، محاسبه شود و سپس تبدیل مختصات روی آن انجام گیرد.

اسلاید 16: - بر این اساس دستگاه مختصات محلی عضو (Local Coordinate System) تعریف می شود. در این دستگاه محلی محور x منطبق بر محور طولی عضو ij (Longitudinal Axis) می باشد و دو محور دیگر منطبق بر محورهای اصلی مقطع عضو می باشند. - با توجه به اینکه بر مبنای این تعریف بین دو انتهای عضو فرقی گذارده نمی شود، از اینرو مشکل مربوط به نوع نیروهای داخلی (کششی، فشاری و ... ) که در دو انتهای عضو بوجود می آید از بین می رود. - حال ماتریس سختی یک عضو را می توان نسبت به دستگاه مختصات محلی بدست آورد. به عبارت دیگر در این حالت ماتریس سختی عضو با وارد کردن تغییر مکان در امتداد هر یک از محورهـای مختصات محلی و محاسبه نیروی موردنیاز برای ایجاد آن در امتدادمحور محلی و نیروهای مورد نیاز برای جلوگیری از تغییر مکان در سایر امتدادها (در دو انتهای عضو) ایجاد می گردد. برای یافتن نیروهای مذکور تنها از دو روش استفاده می شود: (1) استفاده از قضایای انرژی؛ (2) استفاده از روش شیب- افت ” با توجه به سهل و جامع بودن روش شیب- افت، در استخراج ماتریس سختی عضو در مختصات محلی از این روش بهره گرفته شده است.“

اسلاید 17: تغییر مکان محوری (δ1) در جهت محور x دوران محوری- پیچش (δ4) حول محور x- فرض کنید و نیروها و تغییر مکان های حاصل در دو انتهای یک عضو در دستگاه مختصات محلی باشند. زمانی که جهت این مؤلفه ها با جهت محورهای مختصات یکی باشند، مقدار مثبت داشته در غیر این صورت منفی خواهند بود.- اکنون با استفاده از روابط شیب-افت رابطه بین این دو مجموعه از مقادیر را به دست می آوریم. (ابتدا تغییر مکان های را به گره iاعمال می کنیم)

اسلاید 18: تغییر مکان (δ2) در جهت محور y ، در صفحه x-y دوران (δ6) حول محورz ، در صفحه x-y تغییر مکان (δ3) در جهت محور z، در صفحه x-z

اسلاید 19: دوران (δ5) حول محور y ، در صفحه x-z ( اکنون تغییر مکان های را به گره jاعمال می کنیم )- براساس قضیه متقابل ماکسول، نیروهای حاصل در انتهای i از عضو ij تحت اثر یک تغییر مکان در گره j برابر نیروهای حاصل در j تحت اثر همان مقدار تغییر مکان در گره i است. به عبارت دیگراز نظر عددی:- با توجه به نوع قراردادی که برای محورهای مختصات محلی در نظر گرفتیم، برمبنای آن تعریف بین دو انتهای عضو فرقی گذارده نمی شود به عبارت دیگر داریم: (از نظر عددی)

اسلاید 20:

اسلاید 21: - بنابراین داریم:نیروهای ناشی از تغییر مکان گره iنیروهای گره iنیروهای ناشی از تغییر مکان گره jنیروهای ناشی از تغییر مکان گره iنیروهای گره jنیروهای ناشی از تغییر مکان گره j- معادله ، معادله ماتریس سختی عضو در دستگاه مختصات محلی نامیده می شود.- برای یک عنصر در قاب صلب فضایی ماتریس سختی عضو 12×12 خواهد بود.- زمانی که عضوی از سازه دارای شش درجه ازادی نباشد، بعضی از سطرها و ستون های ماتریس حذف خواهد شد. به عنوان مثال اعضای قاب های صلب مستوی فقط دارای سه درجه آزادی هستند، بنابراین سطرها و ستون های مربوط به برای این سازه لازم نخواهد بود.- برای اعضای شبکه ها که دارای سه درجه آزادی هستند، فقط سطرو ستون های مربوط به برای این سازه لازم خواهد بود.- برای خرپاها فقط سطر و ستون مربوط به مورد نیاز خواهند بود.

اسلاید 22: 4- تشکیل ماتریس سختی یک سازه: (Assembly of the Structural Stiffness Matrix)- ماتریس سختی سازه در دستگاه مختصات کلی تشکیل می گردد: (1- متعامد (XYZ)، 2- راستگرد، 3- ثابت و 4- اختیاری)- در تشکیل ماتریس سازه از دو اصل مهم استفاده می شود: الف) اصل سازگاری تغییر مکان ها در گره های سازه، ب) تعادل در گره های سازه.- اصل سازگاری تغییر مکان ها در گره های سازه ایجاب می کند که:به عبارت دیگر تغییرمکان های انتهای اعضای متصل به یک گره خاص برابر تغییرمکان آن گره می باشد.- اصل تعادل نیروها در گره های سازه ایجاب می کند که:به عبارت دیگر بر اساس این اصل، نیروهای انتهایی اعضای متصل به یک گره خاص باید برابر نیروی خارجی موثر در آن گره باشند.- در دستگاه مختصات محلی داریم:- اگر دستگاه مختصات محلی را به عنوان دستگاه جدید و دستگاه مختصات کلی را به عنوان دستگاه مختصات قدیمی در نظر بگیریم خواهیم داشت:ماتریس دوران Rij ماتریسی است که یک بردار را از دستگاه مختصات کلی (XYZ) به دستگاه مختصات محلی (xyz) (در گره i) تبدیل می کند. ماتریس دوران Rji ماتریسی است که یک بردار را از دستگاه مختصات کلی (XYZ) به دستگاه مختصات محلی (xyz) (در گره j) تبدیل می کند.

اسلاید 23: - از جایگذاری در معادله مورد نظر داریم:- فرم دقیق ماتریس Rij برای انواع مختلف سازه ها بعداً ارائه خواهد شد.- از جایگذاری روابط مذکور در معادله حاصل از اصل تعادل در گرهi و نیز با ملاحظه اصل سازگاری تغییر مکان ها داریم:- با نوشتن معادله مذکور برای تمام گره ها و ترتیب مناسب آنها رابطه ماتریسی زیر به دست می آید.

اسلاید 24: ماتریس سختی سازه در دستگاه مختصات کلیبردار تغییرمکان های مجهول گره های آزاد و تغییرمکان های معلوم در تکیه گاه هابردار نیروهای معلوم وارد بر گره های آزاد و نیروهای مجهول در تکیه گاه هااکنون کاملاً روشن می شود که چگونه می توان ماتریس سختی سازه را تشکیل داد(در یک برنامه کامپیوتری تحلیل ماتریسی سازه ها).

اسلاید 25: الف) ابتدا تمام گره ها به ترتیب و با شروع از 1 و اضافه نمودن عدد 1 برای هر شماره و به طور پیوسته شماره گذاری می شوند (N گره). در این صورت هر عضو به طور منحصر به فرد توسط دو شماره در ابتدا و انتهای خود مشخص می گردد. بهتر است که تفاوت بین دو شماره مربوط به یک عضو حتی المقدور مینیمم مقدار ممکن را داشته باشد. توجه به شماره گذاری گره ها باعث کاهش عرض نوار ماتریس K و در نتیجه باعث صرفه جویی در انبار کردن اطلاعات در پردازشگر می گردد.ب) دستگاه مختصات کلی را برای سازه تعیین می کنیم.پ) مختصات گره ها در دستگاه مختصات کلی را به عنوان ورودی وارد می دهیم. ت) ابتدا ( i) و انتهای (j) هر عضو را به عنوان ورودی وارد می کنیم. ث) مشخصات هندسی و مکانیکی اعضا را به عنوان ورودی وارد می کنیم(A, E, J, G, IY, IZ).ج) بارهای وارد بر گره ها را در دستگاه مختصات کلی به عنوان ورودی وارد می کنیم. چ) شرایط تکیه گاهی را در دستگاه مختصات کلی به عنوان ورودی وارد می کنیم. ٍح) برای عضو i-j ماتریس های و را تشکیل می دهیم (توجه شود که ). (ماتریس های سختی اعضا در دستگاه مختصات محلی)خ) برای هر انتهای عضو ij، ماتریس های و را تشکیل می دهیم.د) برای هر عضو i-j ماتریس های را تشکیل می دهیم( توجه شود که : (ماتریس های سختی اعضا در دستگاه مختصات کلی)).

اسلاید 26: ذ) یک ماتریس (a*N × a*N) صفر با کلیه درایه های صفر تعریف می کنیم (با N × N بلوک ماتریسی صفر به ابعاد a × a). (a = تعداد درجات آزادی در هر گره)ر) را در محل سطر بلوکی iام و ستون بلوکی iام قرار می دهیم، را در محل سطر بلوکی jام و ستون بلوکی jام قرار می دهیم، در محل سطر بلوکی iام و ستون بلوکی jام قرار می دهیم، در محل سطر بلوکی jام و ستون بلوکی iام قرار می دهیم.د) مراحل (ح) تا (ر) را برای تمامی اعضاء تکرار می کنیم. بدیهی است هنگامی که بین دو گره j , i عضوی وجود نداشته باشد، در این صورت عملا صفر منظور می شود.رسم فلوچارت برنامه

اسلاید 27: 5- اعمال شرایط مرزی: (Imposition of Boundary Conditions) - ماتریس سختی کل سازه K قبل از اعمال شرایط مرزی یک ماتریس ویژه است (دترمینان آن برابر صفر است) و نشانگر این واقعیت می باشد که سازه بدون تکیه گاه ناپایدار است.- معادله ماتریسی ، یک دستگاه معادلات مختلط می باشد که در هر دو طرف آن مقادیر معلوم و مجهول وجود دارند. - بردار نیروی تعمیم یافته P در طرف چپ این معادله شامل نیروهای خارجی معلوم موثر بر گره های آزاد سازه و نیز حاوی عکس العمل های مجهول نیز می باشد.- بردار تغییرمکان های تعمیم یافته درسمت راست این معادله شامل تغییرمکان های مجهول گره های آزاد سازه و نیز حاوی تغییرمکان های معلوم تکیه گاهی (برابر صفر برای تکیه گاه های بدون نشست، برابر مقدار مشخص برای حالت نشست تکیه گاهی، به صورت تابعی برای تکیه گاه های ارتجاعی) می باشد.- فرض کنید که ماتریس سختی سازه به گونه ای تنظیم شده است که گره های 1 تا m مشخص کننده گره های تکیه گاهی (تغییر مکان های معلوم تعمیم یافته (برابر صفر)) باشند و گره هایm+1 تا n نیز بیانگر گره های آزاد سازه باشند ( تا مجهول). بنابراین نیروهای تعمیم یافته تا بیانگر عکس العمل های مجهول سازه و تا بیانگر نیروهای تعمیم یافته وارد بر گره های آزاد سازه می باشند.

اسلاید 28:

اسلاید 29: - رابطه ، رابطه ماتریسی نهایی سازه است؛ زیرا مجهولات تغییر مکانی گره ها را به نیروهای خارجی معلوم موثر در این گره ها مرتبط می سازد.- بنابراین معادله نهایی ماتریسی سازه در حقیقت از حذف سطرها و ستون های مربوط به تغییر مکان های معلوم صفر سازه در تکیه گاه ها به دست می آید. - به عبارت دیگر به این طریق شرایط مرزی نیز ارضاء می شوند (Satisfying the Boundary Condition) یا به سخن دیگر شرایط مرزی اعمال می گردند (Imposition of Boundary Conditions). بنابراین به طور خلاصه وقتی شرایط مرزی برحسب تغییر مکان ها در امتداد محورهای مختصات کلی برابر صفر باشند، معرفی آنها در معادلات ماتریس سختی با حذف سطرها و ستون های مربوطه انجام می گیرد.- روش دیگری نیز برای اعمال شرایط مرزی وجود دارد. در این روش اعضای قطری ماتریس سختی K که مربوط به تغییر مکان های ثابت معلوم تکیه گاهی می باشند ، را در یک عدد بزرگ مانند ضرب می کنیم. برای مثال در درجه آزادی i که برابر صفر است داریم:عضو را در عدد ضرب می کنیم.در این روش به هنگام حل معادلات در برنامه کامپیوتری را نیز مساوی صفر قرار می دهیم.

اسلاید 30: 6- تحلیل خرپاها (Analysis of Trusses)- خرپاها سازه های متشکل از اعضاء مستقیم نسبتاً لاغر هستند که توسط گره های مفصلی بدون اصطکاک به همدیگر متصل شده و فقط در گره ها تحت اثر بارهای خارجی قرار می گیرند (خرپای ایده ال).- در عمل ایجاد گره های مفصلی بدون اصطکاک کار دشوار و یا غیر ممکن است. تفاوت بین یک خرپای ایده ال و یک خرپای حقیقی در این است که اعضای خرپای حقیقی علاوه بر نیروهای محوری، تحت اثر نیروی برشی و لنگر خمشی نیز قرار می گیرند. هرچه لاغری اعضای خرپا بیشتر می شود، این تفاوت کمتر می شود.می توان خرپاها را با فرض گره های صلب نیز تحلیل کرد. این چنین تحلیلی با در نظر گرفتن سختی های خمشی اعضاء، نیروهای ثانویه (برش و لنگر) و همچنین نیروهای اولیه (نیروهای محوری) را نتیجه می دهد. تحلیل تنش های ثانویه در دو مورد پیشنهاد می گردد: سختی خمشی عضوها زیاد باشد، نتایج با دقت بیشتری خواسته شده باشد.- در اینجا صرفاً به تحلیل خرپاهای ایده ال خواهیم پرداخت.

اسلاید 31: دستگاه مختصات محلیگرهگرهدستگاه مختصات کلی1 درجه آزادی2 درجه آزادی3 درجه آزادیخرپای صفحه ایخرپای فضائی- برای هر عضو ij خرپا فقط محور محلی x را در نظر می گیریم.- هر گره صرفاً دارای مؤلفه تغییر مکانی در امتداد محور xها می باشد.

اسلاید 32: - مراحل تشکیل ماتریس سختی یک سازه خرپایی عبارتند از:الف) تعیین ماتریس سختی عضو در دستگاه مختصات محلی:ب) تعیین ماتریس های دوران برای هر عضو:

اسلاید 33: پ) تعیین برای هر عضو (هرکدام ماتریس های 3×3 می باشند):

اسلاید 34: - برای حالت خرپای دو بعدی مسطح داریم:- همچنین می توان عکس العمل های تکیه گاهی را با استفاده از نیروهای اعضای خرپا به دست آورد:- بعد از تشکیل ماتریس سختی و اعمال شرایط مرزی و حل معادلات می توان بردار ∆ را که شامل تغییر مکان های گره های آزاد سازه است، به دست آورد. بعد از تعیین می توان نیروهای اعضای خرپا را به دست آورد.- با استفاده از ماتریس سختی اولیه (بدون اعمال شرایط مرزی) خواهیم داشت:برای بررسی تعادل گره i نیز از رابطه زیر استفاده می کنیم: (بررسی چند مثال)

اسلاید 35: 7- تحلیل قاب های صلب مسطح (Planar Rigid Frames)- در قاب های صلب اعضاء توسط گره های صلب به هم دیگر اتصال یافته اند، یعنی در یک گره زاویه بین اعضاء پس از تغییر شکل تغییر نمی کند.- بارها ممکن است نظیر خرپاها در گره ها وارد شوند و یا مستقیماً بر روی اعضاء اثر نمایند (در این جا صرفاً حالتی را در نظر خواهیم گرفت که بارها بر گره ها وارد می شوند).- وضعیت محورهای مختصات محلی و کلی:

اسلاید 36: - مراحل تشکیل ماتریس سختی یک قاب صلب دوبعدی عبارتند از:الف) تعیین ماتریس سختی هر عضو در دستگاه مختصات محلی:ب) تعیین ماتریس های دوران برای هر عضو:XYZjiα

اسلاید 37: - از روی به راحتی می توان را نیز بدست آورد:پ) تعیین برای هر عضو (هرکدام ماتریس های 3×3 می باشند):

اسلاید 38: - همچنین می توان عکس العملهای تکیه گاهی را با استفاده از تغییرمکان های تعمیم یافته و با استفاده از ماتریس سختی اولیه (بدون اعمال شرایط مرزی) بصورت زیر به دست آورد:- بعد از اعمال شرایط مرزی و حل معادلات می توان بردار ∆ را که شامل تغییر مکان های تعمیم یافته گره های آزاد سازه است، به دست آورد. بعد از تعیین می توان نیروهای انتهای اعضاء را به دست آورد.- برای بررسی تعادل گرهها از رابطه زیر استفاده می کنیم:ت) به همان روال قبلی ماتریس سختی سازه را تشکیل می دهیم.نیروی محورینیروی برشیلنگر خمشی(بررسی چند مثال)

اسلاید 39: 8- تحلیل شبکه ها (Grids)- در قاب های مسطح بارها در صفحه قاب وارد می شوند. به عبارت دیگر مثلاً اگر قاب مسطح در صفحه XY باشد، در این صورت نیروهای وارد بر آن نیز در صفحه XY و لنگرهای خمشی وارده حول محور Z خواهند بود.- شبکه ها قاب های مسطحی هستند که در آنها بارها بصورت قائم بر صفحه سازه اثر می کنند (بحثی در مورد مؤلفه های نیرو در قاب های فضایی). به عبارت دیگر مثلاً اگر شبکه در صفحه XY باشد، در این صورت نیروهای وارد بر آن در امتداد محور Zها و لنگرهای وارده حول محورهای Y , X خواهد بود.- با توجه به این که بارهای خارجی قائم بر صفحه سازه اثر می کنند، لذا تغییر شکل های محوری قابل صرف نظر کردن می باشند.- یک گره آزاد سازه شبکه ای علاوه بر تغییر مکان های قائم بر صفحه سازه شبکه ای (مثلا در امتداد محور Zها در صورتی که سازه شبکه ای در صفحه XY واقع باشد) تحت اثر دورانی در صفحه خود سازه با مولفه های قرار می گیرد. - بنابر این تغییر شکل یک گره به صورت [∆Z , θX , θY] قابل بیان است، مشروط بر این که صفحه سازه به صورت XY مشخص شده باشد.

اسلاید 40: - وضعیت محورهای مختصات محلی و کلی:

اسلاید 41: - مراحل تشکیل ماتریس سختی یک شبکه عبارتند از:الف) تعیین ماتریس سختی عضو:ب) تعیین ماتریس های دوران :

اسلاید 42: پ) تعیین برای هر عضو (هرکدام ماتریس های 3×3 می باشند):ت) به همان روال قبلی ماتریس سختی سازه را تشکیل می دهیم.- همچنین می توان عکس العمل های تکیه گاهی را با استفاده از تغییرمکان های تعمیم یافته و با استفاده از ماتریس سختی اولیه (بدون اعمال شرایط مرزی) بصورت زیر به دست آورد:- بعد از اعمال شرایط مرزی و حل معادلات می توان بردار ∆ را که شامل تغییر مکان های تعمیم یافته گره های آزاد سازه است، به دست آورد. بعد از تعیین می توان نیروهای انتهای اعضاء را به دست آورد.- برای بررسی تعادل گره ها از رابطه زیر استفاده می کنیم:برشلنگر پیچشیلنگر خمشی(بررسی چند مثال)

اسلاید 43: 9- تحلیل قاب های صلب سه بعدی (Three Dimensional Rigid Frames)- ویژگی های این سازه ها عبارتند از: * سازه و بارهای مؤثر بر آن در فضای فیزیکی سه بعدی قرار دارند.* تمام اعضا (بجز تکیه گاه ها) به صورت صلب به همدیگر مرتبط یافته اند.- محورهای مختصات محلی و کلی

اسلاید 44: - مراحل تشکیل ماتریس سختی قاب های صلب سه بعدی، در هنگامی که تشکیل ماتریس سختی سازه برای حالت کلی در ابتدای این فصل شرح داده می شد، ارائه گردید. نکته ای که فقط باید تعیین گردد تعیین ماتریس دوران می باشد. - تعیین ماتریس دوران : - مشخص است که ماتریس دوران یک ماتریس 6×6 خواهد بود که در حالت کلی صورت زیر را دارد:

اسلاید 45: - به ترتیب کوسینوس های هادی محور محلی x نسبت به محورهای کلی Z, Y, X می باشد.- اکنون باید کوسینوس های محورهای محلی z, y را نسبت به محورهای کلی Z, Y, X بدست آوریم:محورy عمود بر محور x و Z انتخاب می شود به گونه ای که حاصل برداری Z با x، محور y رانتیجه بدهد:( بردارهای یکه در امتداد محورهای کلی Z, Y, X می باشند)- بنابراین کوسینوس های محور محلی y عبارتند از:

اسلاید 46: - حال می توان کوسینوس های هادی محور z را نیز بدست آورد:- بنابراین کوسینوس های محور محلی z عبارتند از:- پس ماتریس دوران به صورت زیر بدست می آید:- شرط استفاده از این ماتریس آن است که محورمحلی x منطبق بر محور Z کلی نشود چرا که در این صورت عضو موازی محور Z بوده و D=o خواهد بود و نمی توان y را تعیین نمود.

اسلاید 47: - در حالت خاصی که محور x محلی موازی محور Z در می آید، در این صورت محور y را می توان به عنوان محور Y از دستگاه مختصات کلی انتخاب کرد. در این صورت داریم:

اسلاید 48: محورهای مختصات کلی (قدیم)محورهای مختصات محلی (جدید)محورهای مختصات اصلی- در قاب های سه بعدی مسأله خاصی ایجاد می شود که عبارت است از عدم انطباق محورهای مختصات محلی تعریف شده با های فوق الذکر با محورهای اصلی مقطع عضو. بنابراین در واقع تفاوت بین محورهای مختصات محلی (x, y, z) با محورهای مختصات اصلی (xp, yp, zp) برابر زاویه می باشد. در این صورت یک تبدیل دورانی دیگری لازم خواهد شد.

اسلاید 49: - نکته قابل توجه آن است که اگر محورهای مختصات اصلی و محلی بر همدیگر منطبق شوند در این صورت خواهد بود و می گردد و لذا زیرنویس p از معادلات حذف می شود.- پس برای قاب های صلب سه بعدی اطلاعات ورودی برای تعریف هندسه سازه عبارتند از: مختصات گره ها نحوه اتصال اعضا (ij) زاویه β برای هر عضو

اسلاید 50: 10- خواص ماتریس های سختی (Properties of the Stiffness Matrices)الف) در ماتریس سختی زیرماتریس های غیر قطری مساوی ترانسپوز همدیگر هستند، یعنی:اما با توجه به این قبلاً ثابت کرده ایم که ب)عناصر ماتریس سختی K مقادیری عددی هستند.طبق قضیه کاستیلیانو

اسلاید 51: پ) ماتریس سختی نهایی سازه (بعد از اعمال شرایط مرزی) یک ماتریس مثبت- معین است (Positive-Definite).- ابتدا باید تعاریف زیر را بیان کنیم:1- ماتریس مثبت- معین: ماتریس مثبت- معینA، ماتریسی است که پیش ضرب و پس ضرب آن با یک بردار دلخواه غیرصفر (B) یک مقدار مثبت می باشد: . می توان اثبات نمود که دترمینان ماتریس A مثبت است |A|>0 .2- ماتریس مثبت- نیمه معین (Positive-Semi Definite): ماتریس مثبت- نیمه معین A، ماتریسی است که پیش ضرب و پس ضرب آن با یک بردار دلخواه غیرصفر (B) مساوی صفر می باشد، می توان اثبات نمودکه دترمینان ماتریس A صفر است |A|=0.3- ماتریس نامعین (Indefinite): ماتریس نامعین A ماتریسی است که پیش ضرب و پس ضرب آن با یک بردار دلخواه غیر صفر (B) یک مقدار منفی می باشد. می توان اثبات نمود که دترمینان ماتریس A کوچکتر از صفر است |A|<0.می توان اثبات نمودکه ماتریس سختی نهایی سازه (بعد از اعمال شرایط مرزی) یک ماتریس مثبت- معین است: معادله ماتریس سختی نهایی انرژی تغییر شکل سیستم (چون K متقارن است) U بیانگر انرژی تغییر شکل سیستم است که یک مقدار مثبت می باشد و ∆ اختیاری است، بنابراین K یک ماتریس مثبت- معین است ، و لذا داریم: |K|>O

اسلاید 52: - ماتریس سختی سازه (قبل از اعمال شرایط مرزی)، یک ماتریس مثبت – نیمه معین است زیرا داریم:که متناظر با یک مد صلب جسمی (Rigid Body Mode) است که در اثر تغییر مکان کاری در سیستم انجام نمی شود. چرا که سازه می تواند یک حرکت صلب جسمی (Rigid Body Motion) انجام دهد. یادآوری می شود که برای یک جسم صلب داریم:” اگر یک سیستم نیرویی بر یک جسم صلب در حال تعادل باشد، در اثر تغییر مکان کوچک (مجازی)، کار خارجی (مجازی) انجام یافته توسط این نیروها برابر صفر خواهد بود.“- بنابراین برای یک ماتریس سختی (قبل از اعمال شرایط مرزی) داریم: |K|=0ت) می توان نشان داد که هیچکدام از ماینورهای اصلی ماتریس سختی نهایی K برابر صفر نخواهد بود. اگر ماتریس سختی نهایی K را به صورت زیر افراز کنیم:در این صورت مثبت معین می باشند:

اسلاید 53: - به همین ترتیب می توان ثابت کرد که نیز یک ماتریس مثبت- معین است.- مشخص است که به همین ترتیب می توان ثابت نمود که:بنابر این تمام درایه های قطری ماتریس سختی نهایی سازه (بعد از اعمال شرایط مرزی) مثبت می باشند. ث) ماتریس های مثبت- معین، مثبت- نیمه معین و نامعین را به صورت زیر نیز تعریف می کنند:- ماتریس مثبت- معین، ماتریسی است که ویژه مقادیر(Eigenvalues) آن همگی مثبت می باشند.- ماتریس مثبت- نیمه معین، ماتریسی است که ویژه مقادیر آن مساوی یا بزرگتر از صفر می باشند.- ماتریس نامعین ماتریسی است که ویژه مقادیر آن منفی، صفر و مثبت می توانند باشند.- ویژه مسأله استاندارد برای ماتریس سختی K را می توان به صورت زیر نوشت: ( یکی از ویژه مقادیر ماتریس K و یکی از ویژه بردارهای(Eigenvectors) ماتریس K می باشند) جواب های این مسأله ویژه جفت های و می باشند.

اسلاید 54: چون اشاره کردیم که ماتریس سختی نهایی یک ماتریس مثبت- معین می باشد، بنابراین تمامی ویژه مقادیر آن مثبت می باشند. با توجه به این که برای ماتریس سختی (قبل از اعمال شرایط مرزی) داریم |K|=0، بنابراین یکی یا تعدادی از ویژه مقادیر ماتریس سختی صفر است: m مساوی تعداد مدهای صلب جسمی است (Rigid body modes):

اسلاید 55: ج) ماتریس سختی بصورت ماتریس نواری است. به عبارت دیگر عناصر غیر صفر در اطراف قطر اصلی هستند، مشروط بر این که تمام گره ها به همدیگر متصل نشده باشند. شماره گذاری طوری انجام شود که تفاضل بین دو شماره مشخص یک عضو به حداقل ممکن محدود شده باشد. یعنی تفاوت بین دو شماره مربوط به یک عنصر حتی المقدور مینیمم مقدار ممکن را داشته باشد.3= ماکزیمم تفاوت بین شماره دو انتهای اعضاتعداد درجات آزادی×[1+(ماکزیمم تفاوت)×2]= عرض نوار ماتریس K=[(2)(3)+1]3=21 4= ماکزیمم تفاوت بین شماره دو انتهای اعضا تعداد درجه آزادی×[1+(ماکزیمم تفاوت)×2]= عرض نوار ماتریسK=[(2)(4)+1]3=27 بنابراین توجه به شماره گذاری گره ها، باعث کاهش عرض نوار ماتریس K و بالنتیجه صرفه جویی در انبار نمودن اطلاعات در ماشین می گردد.

اسلاید 56: - مثالی دیگر در مورد تاثیر شماره گذاری در عرض نوار ماتریس سختی سازه:5= ماکزیمم تفاوت بین شماره اعضاء 33= عرض نوار ماتریس K 2= ماکزیمم تفاوت بین شماره اعضاء15= عرض نوار ماتریس K