شبیه‌سازی عددی رفتار الاستوپلاستیک و رشد آسیب فلزات در تنش‌های سه‌محوری مختلف

نظریه مکانیک آسیب محیط‌های پیوسته با رویکردی پدیدار شناختی در پیش‌بینی پدیده‌ی شکست مواد به کارگرفته می‌‌شود. از آ‌ن‌جایی که هر ماده رفتار مخصوص به خود را دارد، یافتن نحوه رفتار آسیب فلزات پرکاربرد اهمیت ویژه‌ای دارد. در این تحقیق دو مدل آسیب گنجیانی و بونورا، برای تحلیل رفتار الاستوپلاستیک ‌آسیب چند فلز از کدنویسی در نرم‌افزار آباکوس، انجام شبیه‌سازی و مقایسه با داده‌های تجربی موجود در منابع، استفاده شده است. جزئیات انتگرال‌گیری صریح و الگوریتم حل ارائه شده است. همچنین، به دلیل اینکه کرنش‌ها در محدوده تغییرشکل‌های بزرگ بوده‌اند، روابط مربوط به تغییرشکل‌های بزرگ استفاده شده ‌است. در یافتن ثوابت مدل‌ها از نمودارهای تنش-کرنش، آسیب-کرنش و کرنش شکست-تنش سه‌محوری استفاده شده است. پس از انجام شبیه‌سازی به کمک نمودار نیرو-جابه‌جایی که تطابق خوبی با نمودار تجربی برای هر فلز داشته‌اند، صحت عملکرد تایید شده است. در حالت کلی ثوابت طی یک فرایند آزمون و خطا طوری مشخص گردیدند که برازش منحنی بهینه انجام گرفته باشد. فلزات مورد بررسی فولاد 1045، آلومینیوم 2024-تی‌351 و فولاد اچ‌وای 130 بودند. برای دو فلز اول آزمون کشش ساده و برای فولاد اچ‌وای 130 آزمون خمش سه نقطه‌ای شبیه‌سازی شده است و نمودار نیرو-جابه‌جایی حاصل با نمودار تجربی نظیر مقایسه شده است.

Numerical Simulation of Elastoplastic Behavior and Damage Evolution at Various Stress Triaxiality

The theory of continuum damage mechanics with a phenomenological approach is able to simulate phenomena such as soft strain, local necking of materials, and their failure. Stress triaxiality is defined as the stress state in a material that strongly affects the ductile failure phenomena. In this study, two damage models, Ganjiani and Bonora, are chosen to simulate and compare the elastoplastic behavior as well as damage evolution of some metals. These damage models are sensitive to the stress triaxiality. In order to validate the capability of the models in structural response, the proposed model has been implemented into userdefined subroutines VUMAT in the finite element program ABAQUS/Explicit. For this purpose, the explicit stress integration algorithms of the model have been explained. The model has been validated by comparing the predicted results with experimental data. The simulations are performed for steel 1045, aluminum 2024T351, and steel HY130. The details of the integration algorithm in the framework of the explicit scheme are presented. Also, the model is developed in the large strain deformation. For the determination of the constants in the models, the stressstrain, the damagestrain, and the fracture straintriaxiality curves are used. The predicted curves of loaddisplacement from simulation have good agreement with corresponding experimental data.