خطاهای رایج در تحلیل‌های Explicit Dynamics و راه‌های برطرف کردن آن‌ها

تحلیل و عیب‌یابی خطاهای رایج

خطاهای مرتبط با گام زمانی (Timestep) و پایداری حل

گام زمانی (Timestep) مهم‌ترین پارامتر حاکم بر پایداری و کارایی تحلیل‌های Explicit است. این گام زمانی توسط شرط پایداری کورانت-فردریش-لووی (CFL) کنترل می‌شود که مقدار آن مستقیماً به ابعاد کوچک‌ترین المان در مدل و سرعت صوت در ماده وابسته است. به همین دلیل، عدم مدیریت صحیح مش و گام زمانی، یکی از اصلی‌ترین خطاهای Explicit Dynamics محسوب می‌شود که می‌تواند تحلیل را بی‌فایده یا متوقف کند. در ادامه، به بررسی مشکلات ناشی از این پارامتر کلیدی می‌پردازیم.

مشکل: گام زمانی بسیار کوچک و زمان حل طولانی

• تشریح مشکل: در تحلیل Explicit، گام زمانی پایدار برای کل مدل بر اساس کوچک‌ترین المان موجود در مش تعیین می‌شود. بنابراین، وجود حتی تعداد اندکی المان بسیار کوچک، بدشکل یا کشیده، گام زمانی کل تحلیل را به‌شدت کاهش می‌دهد. این امر منجر به افزایش چشمگیر تعداد چرخه‌های حل (Cycles) و در نتیجه، زمان محاسباتی غیرمنطقی و بسیار طولانی می‌شود.
• ارائه راه‌حل‌ها: برای مقابله با این مشکل، راه‌حل‌های زیر قویاً توصیه می‌شود:
  1. بهبود کیفیت مش: بهترین راهکار برای دستیابی به یک گام زمانی مناسب، تولید یک مش باکیفیت و یکنواخت است. تلاش کنید تا توزیع اندازه المان‌ها در سراسر مدل، به‌ویژه در نواحی بحرانی، همگن باشد. استفاده از المان‌های هشت‌وجهی (Hexahedral) به دلیل ساختار منظم، بر المان‌های چهاروجهی (Tetrahedral) ارجحیت دارد. برای دستیابی به مش یکنواخت‌تر، به‌ویژه در مدل‌های پیچیده، استفاده از روش‌های مش‌زنی مستقل از پچ (Patch Independent) توصیه می‌شود. علاوه بر این، با وجود اینکه المان‌های هشت‌وجهی مرتبه اول بهترین عملکرد را ارائه می‌دهند، مستعد بروز پدیده ساعت شنی (Hourglassing) هستند—یک حالت تغییر شکل غیرفیزیکی که در بخش ۲.۲ به تفصیل بررسی خواهد شد. بنابراین، یک مش هشت‌وجهی باکیفیت، اولین و حیاتی‌ترین گام در کنترل گام زمانی و دقت حل است.
  2. استفاده محتاطانه از Mass Scaling: قابلیت Mass Scaling با افزایش مصنوعی جرم المان‌های کوچک، گام زمانی پایدار را افزایش داده و زمان حل را کاهش می‌دهد. با این حال، این یک ابزار صرفاً عددی است و باید با احتیاط فراوان استفاده شود.

مشکل: توقف تحلیل به دلیل رسیدن به حداقل گام زمانی

• تشریح مشکل: در حین شبیه‌سازی، المان‌ها ممکن است به دلیل تغییر شکل‌های بسیار شدید (مانند فشردگی در برخورد یا کمانش) دچار اعوجاج شوند. این اعوجاج باعث کوچک شدن ابعاد مشخصه المان و در نتیجه کاهش شدید گام زمانی پایدار آن می‌شود. اگر این گام زمانی از مقدار حداقل تعریف‌شده در تنظیمات حلگر (Minimum Time Step) کمتر شود، تحلیل برای جلوگیری از ناپایداری به صورت خودکار متوقف خواهد شد.
• ارائه راه‌حل‌ها:
  • بررسی علت اعوجاج: اولین قدم، شناسایی علت اصلی تغییر شکل شدید المان‌هاست. آیا این اعوجاج ناشی از یک پدیده فیزیکی واقعی مانند شکست ماده است یا به دلیل تعریف نادرست بارها یا تماس‌ها رخ داده است؟ بازنگری مدل بر اساس این یافته‌ها ضروری است.
  • استفاده از کنترل‌های Erosion: در بخش Erosion Controls، می‌توان با فعال‌سازی گزینه On Minimum Element Time Step، به حلگر اجازه داد المان‌هایی که گام زمانی آن‌ها از حد معینی کمتر می‌شود را به صورت خودکار از تحلیل حذف کند. این کار از توقف شبیه‌سازی جلوگیری می‌کند؛ با این حال، باید توجه داشت که این یک تکنیک صرفاً عددی است و یک مدل شکست مبتنی بر فیزیک نیست. استفاده نادرست از آن می‌تواند منجر به از دست رفتن زودهنگام جرم و انرژی از سیستم شده و بر دقت نتایج تأثیر منفی بگذارد.

خطای Hourglassing (پدیده ساعت شنی)

پدیده Hourglassing یا «حالت انرژی صفر» یک حالت تغییر شکل غیرفیزیکی است که مختص المان‌های انتگرال‌گیری کاهش‌یافته (reduced-integration) مانند المان‌های هشت‌وجهی (hexahedral) و چهارضلعی (quadrilateral) مرتبه اول است. در این حالت، گره‌های المان بدون اینکه کرنشی در نقطه انتگرال‌گیری ثبت شود، جابجا می‌شوند. در نتیجه، المان بدون ایجاد هیچ‌گونه تنش دچار اعوجاج شده و انرژی کرنشی صفر خواهد بود که منجر به نتایج کاملاً نادرست و غیرواقعی می‌شود. در ادامه، روش‌های شناسایی و کنترل این پدیده تشریح می‌شود.

شناسایی و علل بروز Hourglassing

• نحوه شناسایی: اصلی‌ترین و واضح‌ترین نشانه بروز این خطا، مشاهده مقادیر بالای «انرژی ساعت شنی» (Hourglass Energy) در خروجی‌های انرژی حل (Energy Summary) است. اگر این انرژی بخش قابل توجهی (مثلاً بیش از ۵-۱۰ درصد) از انرژی داخلی سیستم را تشکیل دهد، نتایج شبیه‌سازی فاقد اعتبار است. نسبت انرژی ساعت شنی به انرژی داخلی را می‌توان به سادگی به صورت زنده با بررسی نمودار «Energy Summary» در شاخه Solution Information پایش کرد.
• علل بروز: این پدیده به طور خاص در المان‌های مرتبه اول با انتگرال‌گیری کاهش‌یافته رخ می‌دهد. این نوع المان‌ها به دلیل کارایی محاسباتی بسیار بالا و جلوگیری از قفل‌شدگی (Locking)، به طور گسترده در تحلیل‌های Explicit Dynamics استفاده می‌شوند و درک نحوه کنترل حالت‌های انرژی صفر آن‌ها ضروری است.

روش‌های کنترل و حذف Hourglassing

برای مقابله با پدیده ساعت شنی، می‌توان از مجموعه‌ای از راهکارهای زیر بهره برد:

1. استفاده از Damping: این روش اصلی‌ترین و مؤثرترین راهکار برای کنترل Hourglassing است. Ansys دو فرمولاسیون اصلی برای این منظور ارائه می‌دهد:
   • Standard (Viscous): این فرمولاسیون بر اساس سرعت گره‌ها عمل کرده و از نظر محاسباتی کارآمدتر است. این گزینه برای بسیاری از کاربردها انتخاب پیش‌فرض و مناسبی است.
   • Flanagan-Belytschko: برای تحلیل‌هایی که شامل چرخش‌های قابل توجه هستند (مانند غلتیدن یک قطعه پس از ضربه)، فرمولاسیون Standard می‌تواند به اشتباه به دلیل چرخش جسم صلب، انرژی ساعت شنی تولید کند. فرمولاسیون Flanagan-Belytschko تحت این چرخش‌ها ثابت باقی می‌ماند (invariant) و بنابراین قویاً توصیه می‌شود تا اطمینان حاصل شود که انرژی ساعت شنی محاسبه‌شده ناشی از اعوجاج واقعی المان است و نه اثرات عددی کاذب.
2. بهبود کیفیت مش: استفاده از یک مش متراکم‌تر و با توزیع یکنواخت‌تر می‌تواند به کاهش این پدیده کمک کند، زیرا اعوجاج در سطح یک المان منفرد کمتر به کل مدل سرایت می‌کند.
3. نحوه اعمال بار: از اعمال بارهای نقطه‌ای (Point Load) بر روی یک گره منفرد خودداری کنید. بهتر است بار مورد نظر بر روی چندین گره مجاور توزیع شود تا از تمرکز نیرو و ایجاد حالت‌های تغییر شکل موضعی جلوگیری گردد.
4. بازنگری در اتصالات: در شرایط خاص، تعریف اتصالات می‌تواند منجر به بروز Hourglassing شود. به عنوان مثال، برای اتصال قطعاتی که با المان‌های چهاروجهی (Tetrahedral) مش‌بندی شده‌اند، استفاده از «نقاط راه دور» (Remote Points) به جای «تماس چسبیده» (Bonded Contact) می‌تواند راهکار مناسب‌تری برای جلوگیری از این خطا باشد.

خطاهای مربوط به اعوجاج شدید المان و Erosion

در تحلیل‌های Explicit، تغییر شکل‌های بسیار بزرگ می‌توانند باعث واژگونی المان‌ها (ایجاد حجم منفی) و در نتیجه توقف حل شوند. ابزار عددی Erosion برای حذف این المان‌های به‌شدت تغییرشکل‌یافته از تحلیل طراحی شده است تا شبیه‌سازی بتواند ادامه یابد. با این حال، باید به خاطر داشت که Erosion یک فرآیند فیزیکی نیست و تنظیمات نادرست آن خود می‌تواند یک خطا محسوب شده و دقت نتایج را به خطر اندازد.

مشکل: تنظیمات نادرست Erosion

• تشریح مشکل: استفاده بی‌رویه یا با تنظیمات نادرست از Erosion می‌تواند منجر به حذف زودهنگام یا دیرهنگام المان‌ها و در نتیجه ایجاد نتایجی شود که با واقعیت فیزیکی تطابق ندارند. برای مثال، اگر مقدار Geometric Strain Limit (حد کرنش هندسی) بیش از حد پایین تنظیم شود، المان‌ها پیش از آنکه به حد واقعی شکست برسند، از تحلیل حذف می‌شوند. برعکس، اگر این مقدار بیش از حد بالا باشد، ممکن است نتواند از واژگونی المان و توقف حل جلوگیری کند.
• ارائه راه‌حل‌ها: برای تنظیم صحیح Erosion، موارد زیر را در نظر بگیرید:
  • یک مقدار مناسب برای Geometric Strain Limit تنظیم کنید. مقدار پیش‌فرض ۱.۵ برای بسیاری از کاربردها مناسب است، اما برای مواد بسیار انعطاف‌پذیر مانند مواد هایپرالاستیک، این مقدار باید افزایش یابد.
  • برای شبیه‌سازی واقع‌گرایانه شکست، گزینه On Material Failure را فعال کنید. این گزینه المان‌ها را بر اساس مدل شکست تعریف‌شده در خواص ماده (مانند حد کرنش پلاستیک) حذف می‌کند که رویکردی مبتنی بر داده‌های مهندسی است.
  • گزینه Retain Inertia of Eroded Material را فعال نگه دارید. این گزینه اهمیت زیادی دارد، زیرا با حفظ اینرسی و مومنتم گره‌های آزادشده پس از حذف المان‌ها، تضمین می‌کند که قطعات جداشده (debris) به درستی در برخوردهای بعدی شرکت کرده و مومنتم سیستم به درستی حفظ شود.

خطاهای مربوط به تعریف تماس (Contact)

تنظیمات نادرست تماس یکی از متداول‌ترین منابع خطا در شبیه‌سازی‌های Explicit Dynamics است. این خطاها می‌توانند باعث نفوذ غیرفیزیکی اجسام در یکدیگر، عدم همگرایی (در موارد خاص) یا خطاهای اجرایی در حلگر شوند.

مشکل: نفوذ اولیه (Initial Penetration)

• تشریح مشکل: وجود هرگونه هم‌پوشانی یا نفوذ اولیه بین گره‌ها و المان‌های اجسام مختلف در ابتدای تحلیل، به‌ویژه در تماس‌های لغزشی (Sliding Contact)، باید به حداقل رسیده یا به طور کامل حذف شود. حلگر Explicit ممکن است نتواند این نفوذهای اولیه را به درستی مدیریت کند و باعث ایجاد نیروهای تماسی بزرگ و غیرفیزیکی در همان ابتدای حل شود.
• ارائه راه‌حل: قبل از شروع حل، مدل را به دقت برای هرگونه نفوذ اولیه بررسی کنید. از ابزارهای موجود در Ansys Workbench برای شناسایی و رفع این هم‌پوشانی‌ها استفاده نمایید تا سطوح تماس در وضعیت اولیه خود یا مماس باشند یا فاصله بسیار کمی از هم داشته باشند.

مشکل: استفاده از نوع تماس یا اتصال نامناسب

• تشریح مشکل: برخی از انواع اتصالات و تماس‌ها که در تحلیل‌های Implicit رایج هستند، در محیط Explicit Dynamics پشتیبانی نمی‌شوند. استفاده از این موارد منجر به بروز خطا در مرحله پیش‌پردازش یا توقف حل خواهد شد.
• ارائه راه‌حل: از استفاده از اتصالات و تماس‌های پشتیبانی‌نشده زیر خودداری کنید:
  • اتصالات از نوع Joints و Beam connections.
  • تماس Bonded برای اجسامی که چارچوب مرجع آن‌ها Eulerian تعریف شده است.
  • تماس Bonded برای اتصال مستقیم مش‌های ساخته‌شده از المان‌های چهاروجهی (Tetrahedral). در این حالت، استفاده از Remote Points یا تعریف قطعات چندجسمی (Multi-body Parts) توصیه می‌شود.

مشکل: تنظیمات نادرست Body Interactions

• تشریح مشکل: تنظیمات اشتباه در بخش Body Interactions، که تماس عمومی بین اجسام را کنترل می‌کند، می‌تواند منجر به نادیده گرفتن تماس‌های حیاتی یا بروز خطا در حلگر شود.
• ارائه راه‌حل: به دو سناریوی کلیدی زیر توجه ویژه داشته باشید:
  • Body Self Contact: اگر احتمال می‌دهید که یک جسم در حین تغییر شکل روی خودش تا بخورد (مانند کمانش یک ورق نازک)، تنظیم این گزینه روی No یک خطای جدی است و تماس داخلی آن جسم نادیده گرفته خواهد شد. برای اطمینان و دستیابی به نتایج قوی‌تر، توصیه می‌شود این گزینه روی Yes تنظیم شود.
  • Proximity Based Contact: این الگوریتم تشخیص تماس نیازمند وجود یک فاصله فیزیکی اولیه بین سطوح در ابتدای تحلیل است. اگر این شرط رعایت نشود و سطوح با یکدیگر مماس یا در هم فرورفته باشند، حلگر با صدور پیام خطا متوقف خواهد شد. برای کاربردهای عمومی، به ویژه زمانی که قطعات در ابتدا با یکدیگر در تماس هستند یا در مدل‌های پیچیده که تضمین وجود فاصله اولیه غیرعملی است، الگوریتم پیش‌فرض Trajectory به عنوان روشی قوی و ارجح توصیه می‌شود.

جمع‌بندی و نتیجه‌گیری

در این مقاله، مهم‌ترین دسته‌بندی‌های خطاهای Explicit Dynamics شامل مدیریت گام زمانی، کنترل پدیده ساعت شنی (Hourglassing)، مدیریت اعوجاج شدید المان از طریق Erosion و تنظیمات دقیق تماس مورد بررسی قرار گرفت. اگرچه این خطاها رایج هستند، اما با درک صحیح اصول فیزیکی و عددی حاکم بر حلگر و تنظیم دقیق پارامترهای مدل، کاملاً قابل پیش‌بینی و مدیریت می‌باشند. به عنوان توصیه نهایی، یک تحلیلگر حرفه‌ای باید همواره پارامترهای کلیدی خروجی را برای اعتبارسنجی و اطمینان از سلامت شبیه‌سازی خود کنترل کند. پایش مستمر نمودارهای بقای انرژی (Energy Conservation)، انرژی ساعت شنی (Hourglass Energy) و مقدار گام زمانی حل (Timestep) بهترین ابزارها برای تشخیص زودهنگام خطاها و تضمین اعتبار نتایج نهایی هستند. برای تبدیل این پایش به یک مکانیزم کنترل کیفیت فعال و خودکار، توصیه می‌شود پارامتر Maximum Energy Error را در بخش Step Controls (به عنوان مثال روی ۰.۱ یا ۱۰٪) تنظیم کنید تا در صورتی که تعادل انرژی غیرفیزیکی شود، شبیه‌سازی به صورت خودکار متوقف گردد.