راهنمای مدلهای مواد پیشرفته در انسیس: پلاستیسیته، ویسکوالاستیسیته و…
آنچه در این مقاله میخوانید
چرا به مدلهای مواد پیشرفته نیاز داریم؟ رفتار خطی در مقابل غیرخطی
تفاوت اساسی بین رفتار مواد خطی و غیرخطی در رابطه میان تنش و کرنش نهفته است. در یک مدل خطی، این رابطه ثابت است (مانند قانون هوک)، اما در دنیای واقعی، بسیاری از مواد پس از عبور از یک آستانه تنش مشخص، دچار تغییر شکل دائمی و غیرقابل بازگشت میشوند. این تمایز در شبیهسازیهای مهندسی از اهمیت بالایی برخوردار است، زیرا نادیده گرفتن رفتار غیرخطی میتواند منجر به پیشبینیهای نادرست از ظرفیت باربری، تغییر شکل و حالت شکست سازه شود.
بر اساس اصول تحلیل سازهای، سه منبع اصلی برای رفتار غیرخطی در یک سیستم وجود دارد:
- غیرخطی بودن مواد (Material non-linearities): این نوع غیرخطی بودن زمانی رخ میدهد که رابطه بین تنش و کرنش در خود ماده خطی نباشد. پلاستیسیته، که در آن ماده پس از رسیدن به تنش تسلیم دچار تغییر شکل دائمی میشود، بارزترین نمونه از این رفتار است.
- غیرخطی بودن هندسی (Geometric non-linearities): این پدیده به دلیل تغییرات بزرگ در شکل، کرنش یا دوران سازه حین بارگذاری به وجود میآید. برای مثال، در تحلیل یک تیر بلند و نازک که تحت بار زیاد دچار خمش قابل توجهی میشود، باید اثرات غیرخطی هندسی را در نظر گرفت.
- تغییر وضعیت (Changing status): این مورد به دلیل تغییرات ناگهانی در سختی سیستم رخ میدهد. مسائل تماس، که در آن دو سطح با یکدیگر برخورد کرده و از هم جدا میشوند، نمونهای کلاسیک از این نوع غیرخطی بودن هستند.
تمرکز اصلی این مقاله بر روی غیرخطی بودن مواد است. این نوع از غیرخطی بودن برای شبیهسازی دقیق رفتار بسیاری از مواد مهندسی، از جمله فلزات، پلیمرها، خاک و بتن، ضروری است. با درک این پدیده، میتوانیم به شبیهسازی دقیق تغییر شکلهای دائمی، سختشوندگی کرنشی و در نهایت، پیشبینی دقیقتر شکست بپردازیم. اکنون که با اهمیت این موضوع آشنا شدیم، در بخش بعدی مفاهیم بنیادی تئوری پلاستیسیته را بررسی خواهیم کرد.
مفاهیم کلیدی در تئوری پلاستیسیته
پلاستیسیته به رفتار تغییر شکل دائمی و غیرقابل بازگشت مواد پس از عبور از یک آستانه تنش مشخص، موسوم به «تنش تسلیم»، اطلاق میشود. برخلاف تغییر شکل الاستیک که با برداشتن بار از بین میرود، تغییر شکل پلاستیک در ماده باقی میماند. برای مدلسازی دقیق این رفتار پیچیده در نرمافزارهای المان محدود مانند انسیس، لازم است سه جزء کلیدی تئوری پلاستیسیته را بشناسیم: معیار تسلیم، قانون سختشوندگی و قانون جریان.
معیار تسلیم (Yield Criterion)
معیار تسلیم یک تابع ریاضی است که سطح تنشی را که در آن تغییر شکل پلاستیک آغاز میشود، مشخص میکند. این معیار یک «سطح تسلیم» را در فضای تنش تعریف میکند؛ تا زمانی که ترکیب تنشها در داخل این سطح قرار داشته باشد، ماده رفتار الاستیک دارد و به محض رسیدن به این سطح، تسلیم آغاز میشود. توجه داشته باشید که این «سطح» در یک فضای تنش چندبعدی تعریف میشود و نمودارهای دوبعدی که برای نمایش آن استفاده میشوند، تنها برشی ساده از یک هندسه پیچیده هستند. دو مورد از رایجترین معیارهای تسلیم که در انسیس نیز کاربرد گستردهای دارند، عبارتند از:
- Von Mises: این معیار که برای فلزات و سایر مواد شکلپذیر بسیار رایج است، فرض میکند که تسلیم زمانی آغاز میشود که انرژی اعوجاجی در یک نقطه به مقدار بحرانی برسد. این معیار به فشار هیدرواستاتیک (تنش میانگین) وابسته نیست.
- Drucker-Prager: این معیار یک نسخه تعمیمیافته از مدل Mohr-Coulomb است که برای مواد حساس به فشار هیدرواستاتیک مانند خاک، سنگ و بتن استفاده میشود. سطح تسلیم در این مدل به شکل یک مخروط در فضای تنش اصلی است و نشان میدهد که با افزایش فشار فشاری، مقاومت ماده نیز افزایش مییابد.
قانون سختشوندگی (Hardening Rule)
قانون سختشوندگی توصیف میکند که سطح تسلیم با افزایش تغییر شکل پلاستیک چگونه تغییر میکند. در بسیاری از مواد، پس از شروع تسلیم، برای ادامه تغییر شکل پلاستیک به تنش بیشتری نیاز است؛ این پدیده «سختشوندگی کرنشی» (Strain Hardening) نام دارد. دو قانون اصلی برای توصیف این رفتار عبارتند از:
- سختشوندگی ایزوتروپیک (Isotropic Hardening): در این مدل، با افزایش کرنش پلاستیک، سطح تسلیم به طور یکنواخت در همه جهات منبسط میشود. این مدل برای بارهای یکنواخت و بدون سیکلهای بارگذاری معکوس مناسب است.
- سختشوندگی سینماتیک (Kinematic Hardening): در این مدل، اندازه سطح تسلیم ثابت باقی میماند، اما مرکز آن در فضای تنش جابجا میشود. این قانون برای شبیهسازی اثر باوشینگر (Bauschinger effect) در بارهای چرخهای (مانند کشش و سپس فشار) که در آن تنش تسلیم در جهت معکوس کاهش مییابد، مناسب است.
قانون جریان (Flow Rule)
پس از آنکه تنش به سطح تسلیم رسید، قانون جریان جهت و مقدار افزایش کرنش پلاستیک را تعیین میکند. این قانون مشخص میکند که بردار افزایش کرنش پلاستیک در چه جهتی نسبت به سطح تسلیم قرار دارد. دو نوع اصلی قانون جریان وجود دارد:
- قانون جریان پیوسته (Associative Flow Rule): در این حالت، جهت بردار افزایش کرنش پلاستیک، عمود بر سطح تسلیم است. این قانون برای فلزات به خوبی عمل میکند.
- قانون جریان غیرپیوسته (Non-Associative Flow Rule): در این حالت، جهت کرنش پلاستیک با جهت نرمال بر سطح تسلیم متفاوت است و توسط یک تابع جداگانه به نام «پتانسیل پلاستیک» تعریف میشود. این رویکرد برای مدلسازی موادی مانند خاک و سنگ که رفتار اتساعی (Dilatancy) متفاوتی از خود نشان میدهند، ضروری است.
با درک این سه جزء بنیادی، اکنون میتوانیم نحوه پیادهسازی این مفاهیم را در قالب مدلهای مواد خاص موجود در نرمافزار انسیس بررسی کنیم.
مدلهای رایج پلاستیسیته در انسیس
انسیس مجموعهای غنی از مدلهای مواد داخلی را برای شبیهسازی رفتار پلاستیک مواد مختلف ارائه میدهد. هنگام انتخاب یک مدل، مهندس شبیهسازی باید ماهیت ماده و شرایط بارگذاری را در نظر بگیرد تا بهترین گزینه را انتخاب کند. در این بخش، دو مدل بسیار پرکاربرد یعنی مدل دوخطی Von Mises و مدل Drucker-Prager را که برای دستههای مختلفی از مواد مهندسی استفاده میشوند، بررسی خواهیم کرد.
مدل پلاستیسیته دوخطی Von Mises با سختشوندگی ایزوتروپیک
برای شبیهسازی رفتار فلزات شکلپذیر یا آرماتورها تحت بارهای یکنواخت، مدل پلاستیسیته دوخطی Von Mises به دلیل سادگی و دقت، اغلب اولین و بهترین انتخاب است. این مدل به طور گسترده برای شبیهسازی آرماتورها در سازههای بتن مسلح استفاده میشود و عملکرد آن به صورت دوخطی است:
- ناحیه الاستیک خطی: ابتدا ماده مطابق با قانون هوک رفتار میکند و رابطه تنش-کرنش توسط مدول یانگ (E) تعریف میشود.
- ناحیه پلاستیک: پس از رسیدن تنش به مقدار تنش تسلیم (Yield Stress)، ماده وارد ناحیه پلاستیک میشود. در این ناحیه، رفتار ماده توسط یک مدول مماسی (Tangent Modulus) تعریف میشود که شیب منحنی تنش-کرنش پس از نقطه تسلیم را مشخص میکند.
در این مدل، برای توصیف رفتار سختشوندگی کرنشی (Strain Hardening)، مقدار مدول مماسی باید مثبت باشد. این مدل از معیار تسلیم Von Mises و قانون سختشوندگی ایزوتروپیک پیروی میکند که آن را برای بارهای یکنواخت مناسب میسازد.
مدل پلاستیسیته Drucker-Prager
مدل Drucker-Prager برای موادی طراحی شده است که مقاومت آنها به فشار هیدرواستاتیک (فشار همه جانبه) وابسته است. این ویژگی آن را به گزینهای ایدهآل برای شبیهسازی مواد اصطکاکی مانند خاک، سنگ، بتن و سایر مواد دانهای تبدیل میکند. پارامترهای کلیدی این مدل عبارتند از:
- سطح تسلیم: این سطح با استفاده از دو پارامتر اصلی ماده تعریف میشود: چسبندگی (Cohesion) که نشاندهنده مقاومت برشی ماده در تنش نرمال صفر است و زاویه اصطکاک داخلی (Internal Friction Angle) که وابستگی مقاومت برشی به تنش نرمال را توصیف میکند.
- قانون جریان: قانون جریان در این مدل از طریق زاویه اتساع (Angle of Dilatancy) کنترل میشود. این زاویه تغییر حجم ماده در حین تغییر شکل برشی پلاستیک را مشخص میکند. اگر زاویه اتساع برابر با زاویه اصطکاک داخلی در نظر گرفته شود، قانون جریان از نوع پیوسته (Associative) خواهد بود.
مدل استاندارد Drucker-Prager در انسیس یک رفتار الاستو-پلاستیک بدون سختشوندگی (perfectly elasto-plastic) را توصیف میکند. این بدان معناست که پس از رسیدن به سطح تسلیم، ماده بدون نیاز به افزایش تنش، به تغییر شکل پلاستیک ادامه میدهد.
این دو مدل پایهای، قابلیتهای گستردهای را برای تحلیلهای غیرخطی فراهم میکنند، اما گاهی برای شبیهسازی مواد پیچیدهتر مانند بتن، به مدلهای تخصصیتر یا ترکیبی نیاز است.
مدلهای تخصصی و ترکیبی: مطالعه موردی مدل بتن
علاوه بر مدلهای عمومی پلاستیسیته، انسیس مدلهای تخصصی را برای مواد خاص با رفتار پیچیده مانند بتن ارائه میدهد. این مدلها به مهندسان اجازه میدهند تا پدیدههایی مانند ترکخوردگی در کشش و خردشدگی در فشار را به طور همزمان شبیهسازی کنند. علاوه بر این، انسیس قابلیت ترکیب مدلهای مختلف را برای دستیابی به شبیهسازیهای دقیقتر فراهم میکند که این یک ابزار قدرتمند برای تحلیلهای پیشرفته است.
مدل داخلی بتن در انسیس (مبتنی بر William-Warnke)
مدل بتن پیشفرض در انسیس بر اساس معیار تسلیم ۵ پارامتری William-Warnke ساخته شده است که برای توصیف رفتار سه-محوری بتن توسعه یافته است. این مدل به دلیل توانایی در پیشبینی دقیق شروع شکست بتن تحت ترکیبات مختلف تنش، بسیار مورد توجه قرار گرفته است. این یک مدل ترکیبی است که رفتارهای مختلف بتن را پوشش میدهد:
- در کشش خالص: از یک معیار برش (cut-off) مبتنی بر تئوری Rankine (تئوری حداکثر تنش اصلی) استفاده میکند. این بدان معناست که ترکخوردگی زمانی رخ میدهد که حداکثر تنش اصلی کششی به مقاومت کششی ماده برسد.
- در فشار دو-محوری و سه-محوری: از معیار کامل William-Warnke برای تعریف سطح تسلیم (یا سطح خردشدگی) استفاده میکند.
برای استفاده از این مدل، حداقل دو پارامتر ورودی کلیدی مورد نیاز است: مقاومت فشاری تکمحوری (fc) و مقاومت کششی تکمحوری (ft). سایر پارامترها مانند مقاومت فشاری دو-محوری، به طور پیشفرض بر اساس روابطی که در مقاله اصلی William و Warnke ارائه شده، محاسبه میشوند، اما کاربر میتواند آنها را برای افزایش دقت، به صورت دستی وارد کند.
ترکیب مدلها: استفاده از Drucker-Prager برای رفتار پس از ترک
یکی از محدودیتهای مدل بتن William-Warnke در انسیس این است که عمدتاً معیار شروع تسلیم (آغاز ترک) را پیشبینی میکند اما رفتار ماده پس از ایجاد ترک را به طور کامل توصیف نمیکند. برای غلبه بر این محدودیت و شبیهسازی رفتارهای پیچیدهتر، میتوان این مدل را با یک مدل پلاستیسیته مانند Drucker-Prager ترکیب کرد.
مفهوم این ترکیب به این صورت است که دو سطح تسلیم به طور همزمان تعریف میشوند:
- سطح تسلیم خارجی: مدل بتن (William-Warnke) به عنوان سطح تسلیم اولیه یا خارجی عمل میکند.
- سطح تسلیم داخلی: سطح تسلیم مدل Drucker-Prager باید به گونهای تعریف شود که کاملاً در داخل سطح William-Warnke قرار گیرد.
با این رویکرد، ماده ابتدا بر اساس معیار بتن تسلیم میشود (ترک میخورد یا خرد میشود). سپس، رفتار پس از تسلیم آن توسط قانون جریان مدل Drucker-Prager کنترل میشود. این تکنیک پیشرفته برای شبیهسازی رفتار شکلپذیر بتنهای مسلح شده با الیاف با عملکرد فوقالعاده بالا (UHPFRC) بسیار حیاتی است. در این مواد، پس از شکست ترد اولیه ماتریس بتن (که توسط معیار William-Warnke پیشبینی میشود)، الیاف فولادی با ایجاد یک پاسخ شبهپلاستیک و جذب انرژی، از باز شدن کامل ترکها جلوگیری میکنند. مدل Drucker-Prager برای تقریب زدن همین رفتار شکلپذیر و انرژیگیر پس از ترکخوردگی به کار میرود و به این ترتیب، کل فرآیند شکست را با دقت بیشتری شبیهسازی میکند.
این رویکرد ترکیبی نشان میدهد که چگونه میتوان با استفاده هوشمندانه از مدلهای مختلف، پدیدههای پیچیده فیزیکی را در شبیهسازیهای المان محدود مدل کرد. در ادامه، از پلاستیسیته فراتر رفته و به سایر رفتارهای غیرخطی مواد خواهیم پرداخت.
فراتر از پلاستیسیته: مروری بر سایر رفتارهای پیشرفته مواد
غیرخطی بودن مواد تنها به پلاستیسیته محدود نمیشود. بسیاری از مواد مهندسی، بهویژه پلیمرها، الاستومرها و کامپوزیتها، رفتارهای پیچیدهتری از خود نشان میدهند که وابسته به زمان، نرخ کرنش یا تغییر شکلهای بزرگ هستند. درک این رفتارها برای شبیهسازیهای تخصصی در صنایعی مانند خودروسازی، هوافضا و پزشکی ضروری است. در ادامه، به طور خلاصه به چند مورد از این مدلهای مواد پیشرفته اشاره میشود:
- هایپرالاستیسیته (Hyperelasticity): این رفتار برای مدلسازی موادی مانند لاستیک، فوم و سایر الاستومرها که میتوانند تغییر شکلهای الاستیک بسیار بزرگ و کاملاً برگشتپذیر را تحمل کنند، استفاده میشود. این مدل برای تحلیل آببندهای لاستیکی (O-rings) یا پایههای موتور که دچار تغییر شکلهای بزرگ میشوند، حیاتی است. رفتار هایپرالاستیک توسط یک پتانسیل چگالی انرژی کرنشی (Strain Energy Density Potential) توصیف میگردد.
- ویسکوالاستیسیته (Viscoelasticity): مواد ویسکوالاستیک، مانند بسیاری از پلیمرها، رفتاری بینابین یک جامد الاستیک و یک سیال ویسکوز از خود نشان میدهند که به نرخ کرنش وابسته است. رفتار آن برای شبیهسازی دقیق داشبورد خودروهای پلیمری که سختی آنها با دمای محیط و سرعت ضربه تغییر میکند، ضروری است. این مدل شامل یک جزء ویسکوز در کنار جزء الاستیک است که پدیدههایی مانند وادادگی تنش و خزش را توصیف میکند.
- خزش (Creep): خزش به تغییر شکل غیرقابل بازگشت و وابسته به زمان یک ماده تحت بار ثابت، بهویژه در دماهای بالا، اطلاق میشود. این پدیده در طراحی قطعاتی که برای مدت طولانی تحت بار و دمای بالا کار میکنند، مانند پرههای توربین گازی یا خطوط لوله بخار، بسیار حائز اهمیت است. مدلهای خزش، نرخ کرنش را به عنوان تابعی از تنش، زمان و دما تعریف میکنند.
- پلاستیسیته وابسته به نرخ (Rate-dependent plasticity): این رفتار که گاهی ویسکوپلاستیسیته نیز نامیده میشود، ترکیبی از پلاستیسیته و وابستگی به نرخ کرنش است. در این حالت، تنش تسلیم و رفتار پس از تسلیم ماده به سرعت اعمال بار بستگی دارد. این پدیده در شبیهسازی فرآیندهای شکلدهی فلزات با سرعت بالا یا تحلیل برخوردها کاربرد دارد.
این مدلها قابلیتهای تحلیل المان محدود را به شکل چشمگیری گسترش میدهند و امکان شبیهسازی دقیق طیف وسیعی از پدیدههای فیزیکی را فراهم میآورند.
جمعبندی و نتیجهگیری
در این راهنما، به بررسی اهمیت و مفاهیم کلیدی مدلسازی رفتار غیرخطی مواد در نرمافزار انسیس پرداختیم. همانطور که نشان داده شد، درک دقیق مدلهای مواد در انسیس—از مفاهیم پایهای تئوری پلاستیسیته (معیار تسلیم، قانون سختشوندگی و قانون جریان) گرفته تا مدلهای تخصصی مانند مدل بتن William-Warnke و سایر رفتارهای پیشرفته—برای دستیابی به نتایج شبیهسازی معتبر و قابل اعتماد، امری ضروری است. مدلهای خطی الاستیک، با وجود سادگی، توانایی توصیف پدیدههای حیاتی مانند تغییر شکل دائمی، سختشوندگی کرنشی و شکست را ندارند.
به عنوان یک توصیه نهایی، باید تأکید کرد که اگرچه مدلهای مواد پیشرفته دقت تحلیل را به شکل چشمگیری افزایش میدهند، اما استفاده از آنها نیازمند درک عمیق از رفتار فیزیکی ماده و آگاهی از محدودیتهای هر مدل است. انتخاب یک مدل نامناسب یا پارامترهای ورودی نادرست میتواند به نتایجی منجر شود که نه تنها غیردقیق، بلکه گمراهکننده نیز باشند. همانطور که در مطالعه موردی بتن مسلح شده با الیاف مشاهده شد، ترکیب هوشمندانه مدلها، مانند استفاده از Drucker-Prager برای شبیهسازی رفتار پس از ترک، راهکاری قدرتمند برای مدلسازی پدیدههای پیچیده است. در نهایت، موفقیت یک تحلیل غیرخطی به ترکیبی از دانش نظری، تسلط بر ابزارهای نرمافزاری و قضاوت صحیح مهندسی بستگی دارد.
