,
,
بدلیل ایزوپارامتریک بودن المان، توابع شکل یکسانی جهت درونیابی جابجاییها و مختصات گرههای المانها بکار میرود. رابطه زیر نشان دهنده این مطلب است.
(۳-۲۳) ,
که اندیس ، در عبارت، دهنده تعداد گرههای المان است که در اینجا تعداد گرهها هشت است.
بردارهای ، به ترتیب، بردار مختصات و جابجایی گرهها است. ماتریس ، ماتریس توابع شکل نامیده می شود.
(۳-۲۴)
برای نگاشت المان احتیاج به یک تبدیل یک به یک است. برای چنین تبدیلهایی در ریاضیات از تعریف ماتریس ژاکوبین استفاده می شود.
(( اینجا فقط تکه ای از متن درج شده است. برای خرید متن کامل فایل پایان نامه با فرمت ورد می توانید به سایت feko.ir مراجعه نمایید و کلمه کلیدی مورد نظرتان را جستجو نمایید. ))
در فرمولاسیون روش المان محدود احتیاج به ماتریس کرنش-جابجایی داریم، که آن را ماتریس مینامیم. در این ماتریس چون با اپراتور گرادیان سر و کار داریم، دستیابی به آن آسان نیست. به این منظور تابع راتعریف کرده و مشتقات آن را نسبت به و محاسبه میکنیم. در واقع در اینجا میتواند یا ، یعنی جابجاییها در جهت و باشد. مشتقات نسبت به و به صورت مستقیم قابل دستیابی نیستند به این منظور از قاعده مشتق زنجیرهای استفاده کرده و مشتقات نسبت به و را حساب میکنیم.
(۳-۲۵)
روابط بالا را در حالت فشرده به صورت زیر میتوان نوشت.
(۳-۲۶)
که ، ماتریس ژاکوبین نامیده می شود.
(۳-۲۷)
ترمهای را بصورت زیر بیان می شود.
(۳-۲۸)
که ، معکوس ماتریس ژاکوبین است.
(۳-۲۹)
که ، دترمینان ماتریس ژاکوبین است.
(۳-۳۰)
در مراجع اکثراً را به عنوان ژاکوبین میشناسند. ژاکوبین در واقع همان ضریب مقیاس است که در ضرب شده و سطح ساخته شده توسط در دستگاه اولیه را، در دستگاه جدید مقیاس می کند. به عنوان یک مثال ساده در تبدیل مختصات کارتزین به مختصات قطبی، ، ضریب همان ژاکوبین است . در کل ژاکوبین تابعی از مولفههای مختصات جدید است.
هدف بدست آوردن جابجاییها نسبت به مولفههای مختصات اصلی، یعنی و است. با توضیحاتی که تا اینجا داده شد، توانستیم مشتقات تابع دلخواهی مثل را نسبت به و بدست آوریم. اکنون میتوانیم ماتریس کرنش-جابجایی و ماتریس سختی را بدست بیاوریم. قبل از بدست آوردن ماتریسهای ذکر شده، یک سری تعاریف اولیه راجع به کرنشها و جابجاییها در فضای دو بعدی ارائه می شود.
در این مطالعه رفتار خاک را در حالت دو بعدی کرنش صفحهای بررسی کردهایم. اکثر مسائل ژئوتکنیکی در حالت کرنش صفحهای تعریف میشوند. در حالت کرنش صفحهای ، مولفه کرنش عمود بر صفحه آنالیز صفر است، در این حالت ماتریس الاستیسیته ، به صورت زیر تعریف می شود.
(۳-۳۱)
که ، مدول الاستیسیته و ، ضریب پوآسون است. البته در ادامه راجع به مدول رفتاری انتخاب شده در این مطالعه توضیح میدهیم.
رابطه کرنش-جابجایی در رابطه زیر آورده شده است.
(۳-۳۲)
که ماتریس ، حاصل ضرب سه ماتریس مستطیلی است. کرنشهای تعریف شده در حالت کرنش صفحهای به شرح زیر هستند.
(۳-۳۳)
ماتریس کرنش به فرم زیر نیز قابل بیان است.
(۳-۳۴)
رابطه زیر از گسترش رابطه ۳-۲۸ برای جابجاییهای مختلف حاصل می شود.
(۳-۳۵)
رابطه فوق نشان دهنده رابطه بین مشتق جابجاییها نسبت به مولفههای دستگاه مختصات و دستگاه مختصات است.
رابطه ۳-۳۶ از مشتق گیری رابطه سوم روابط ۳-۲۳ و با در نظر گرفتن رابطه ۳-۲۴ حاصل می شود.
(۳-۳۶)
با جاگذاری دو رابطه اخیر در رابطه ۳-۳۴ ماتریس ساخته می شود.
(۳-۳۷)
حال با داشتن ماتریس ماتریس سختی در قالب معادله ۳-۳۸ قابل بیان است.
(۳-۳۸)
در رابطه بالا ، ضخامت المان است.
بعد از تعریف ماتریس سختی، ماتریس جرم را تعریف میکنیم. ماتریس جرم در قالب رابطه زیر تعریف می شود.
(۳-۳۹)
در رابطه بالا ، ماتریس توابع شکل و ، چگالی یا جرم واحد حجم مصالح است.
تا اینجا ماتریسهای سختی و جرم تعریف شدند. انتگرالهای موجود در روابط بالا را نمیتوانیم به صورت مستقیم حل کنیم. برای حل آنها باید از تکنیک انتگرال عددی استفاده کنیم. روشی که معمولاً در روش المان محدود برای انتگرال گیری المانهای صفحهای ایزوپارامتریک در نظر گرفته می شود روش کوآدرچر است.