تکینگی یا سینگولاریتی (Singularity) در تحلیل المان محدود چیست و چگونه آن را رفع کنیم؟

سینگولاریتی (تکینگی تنش) چیست؟

درک ماهیت ریاضی سینگولاریتی برای یک مهندس تحلیل‌گر، یک مزیت استراتژیک است. باید به خاطر داشت که این پدیده یک خطای فیزیکی در قطعه واقعی نیست، بلکه یک آرتیفکت (artifact) ریاضی است که از ساده‌سازی‌های مدل‌سازی ما نشأت می‌گیرد. به طور دقیق، سینگولاریتی نقطه‌ای در مدل المان محدود است که در آن، تنش به یک مقدار مشخص همگرا نمی‌شود. با افزایش چگالی مش، مقدار تنش در این نقاط نیز به طور مداوم افزایش می‌یابد و به لحاظ تئوری به سمت بی‌نهایت میل می‌کند.

یکی از کلیدی‌ترین مفاهیم برای درک سینگولاریتی، «اصل سنت ونان» (St. Venant’s Principle) است. این اصل بیان می‌کند که اثر اغتشاشات محلی در یک میدان تنش، به صورت محلی باقی می‌ماند. بر این اساس، سینگولاریتی یک اثر کاملاً محلی است. اما چرا این پدیده نتایج جابجایی را تحت تأثیر قرار نمی‌دهد؟ پاسخ در اساس روش المان محدود نهفته است: در تحلیل المان محدود سازه‌ای، جابجایی‌ها و دوران‌ها مجهولات اصلی هستند که مستقیماً از معادلات ماتریسی سراسری حل می‌شوند. این مقادیر همواره با ریزتر شدن مش به یک جواب منحصربه‌فرد همگرا می‌شوند. در مقابل، تنش‌ها و کرنش‌ها مقادیری مشتق‌شده هستند که از روی این جابجایی‌ها محاسبه می‌گردند و اگر مدل پایه دارای یک تکینگی ریاضی باشد، ممکن است مقادیر آن‌ها همگرا نشوند. بنابراین، میدان تنش در فاصله‌ای معقول از ناحیه تکین، کاملاً معتبر و صحیح است.

در این مرحله، ایجاد یک تمایز بسیار مهم ضروری است:

• تکینگی تنش (Stress Singularity): یک آرتیفکت عددی ناشی از مدل‌سازی ایده‌آل (مانند یک گوشه کاملاً تیز) است که باید شناسایی و مدیریت شود. تنش در این نقاط با ریز شدن مش به سمت بی‌نهایت میل می‌کند و واقعی نیست.
• تمرکز تنش یا خیز تنش (Stress Riser / Concentration): یک پدیده فیزیکی واقعی است که تنش‌های بالا اما محدود و همگراشونده ایجاد می‌کند (مانند تنش در اطراف یک سوراخ یا فیلت). این تنش‌ها کاملاً واقعی بوده و نباید نادیده گرفته شوند، بلکه باید با دقت تحلیل شوند.

حال که با ماهیت سینگولاریتی و تفاوت آن با تمرکز تنش واقعی آشنا شدیم، به بررسی دلایل مشخصی که منجر به بروز این آرتیفکت‌های عددی می‌شوند، می‌پردازیم.

دلایل اصلی بروز سینگولاریتی در تحلیل المان محدود

شناسایی ریشه‌ای علل بروز سینگولاریتی، اولین گام برای مدیریت آن است. این پدیده معمولاً از ساده‌سازی‌هایی (simplifications) نشأت می‌گیرد که برای عملی و بهینه کردن مدل‌های FEA در هندسه، بارگذاری و شرایط مرزی انجام می‌دهیم. در ادامه، رایج‌ترین این دلایل را تحلیل می‌کنیم.

گوشه‌های تیز و فرورفته (Sharp and Re-entrant Corners)

در دنیای واقعی، هیچ گوشه‌ای کاملاً تیز نیست و همواره یک شعاع انحنا (فیلت) هرچند کوچک وجود دارد. مدل‌سازی یک گوشه به صورت کاملاً تیز در نرم‌افزار CAD، یک نقطه تکینگی ریاضی ایجاد می‌کند. با ریز شدن مش در این ناحیه، حلگر تلاش می‌کند تنشی را محاسبه کند که از نظر تئوری بی‌نهایت است و در نتیجه، مقدار تنش به طور مداوم افزایش می‌یابد.

بارها و قیود نقطه‌ای (Point Loads and Constraints)

این پدیده را می‌توان با فرمول ساده تنش = نیرو / سطح (Stress = Force / Area) تحلیل کرد. هنگامی که یک نیرو یا قید به یک گره (node) یا رأس (vertex) اعمال می‌شود، حلگر آن نیرو را بر روی یک مساحت مفهومی که از وجوه المان‌های متصل به آن گره مشتق می‌شود، توزیع می‌کند. با ریز شدن مش، این مساحت مفهومی به سمت صفر میل می‌کند. در نتیجه، حلگر تلاش می‌کند کل بار را روی مساحتی بی‌نهایت کوچک اعمال کند که باعث می‌شود تنش محاسبه‌شده به سمت بی‌نهایت افزایش یابد.

قیود نامناسب و بیش از حد (Improper and Over-Constraints)

استفاده از شرایط مرزی که رفتار واقعی قطعه را نمایندگی نمی‌کنند، می‌تواند باعث ایجاد تنش‌های مصنوعی و غیرواقعی شود. دو مثال رایج عبارتند از:

• جلوگیری از اثر پواسون (Poisson’s effect): استفاده از قید Fixed Support در جایی که قطعه در واقعیت به صورت صلب جوش داده نشده است. این قید مانع از تغییر شکل عرضی طبیعی قطعه (ناشی از اثر پواسون) شده و تنش‌های مصنوعی بالایی در نزدیکی قید ایجاد می‌کند.
• محدودیت انبساط حرارتی: اعمال قید ثابت روی سوراخ‌های یک برد الکترونیکی در تحلیل حرارتی. در حالی که در واقعیت، برای انبساط حرارتی این قطعات فضای خالی یا انعطاف‌پذیری در نظر گرفته شده است، قید Fixed Support این حرکت را محدود کرده و منجر به تنش‌های مصنوعی بالا می‌شود.

تماس با گوشه‌های تیز (Contact with Sharp Corners)

این مورد حالتی خاص است که در آن، نیروی تماسی روی یک لبه یا رأس تیز متمرکز می‌شود. این تمرکز نیرو، مشابه اعمال یک بار نقطه‌ای عمل کرده و باعث ایجاد سینگولاریتی در ناحیه تماس می‌گردد.

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

چگونه سینگولاریتی را در Ansys شناسایی و رفع کنیم؟

استراتژی کلی برای مقابله با سینگولاریتی، یک رویکرد یکسان برای همه مسائل نیست. رویکرد «صحیح» به مدل، هدف تحلیل و موقعیت ناحیه تکین بستگی دارد و مهندس باید با قضاوت مهندسی (engineering judgement) بهترین روش را انتخاب کند. در ادامه چهار رویکرد اصلی و عملی برای مدیریت این پدیده ارائه شده است.

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

اگر ناحیه تکینگی در منطقه مورد علاقه (region of interest) شما قرار ندارد، ساده‌ترین و کارآمدترین راه، نادیده گرفتن آن است. بر اساس اصل سنت ونان، این تنش‌های مصنوعی تأثیری بر نتایج دورتر از خود ندارند. با این حال، این رویکرد کارآمد نیازمند آن است که تحلیل‌گر کاملاً مطمئن باشد که تکینگی هیچ تأثیری بر ناحیه مورد نظر نداشته و هیچ مکانیزم شکستی از محل تکینگی آغاز نخواهد شد. به صورت عملی، می‌توانید با استفاده از ابزارهای انتخاب هندسه (scoping results) در Ansys، نتایج تنش را فقط برای نواحی معتبر مدل نمایش دهید. این کار از تأثیر بصری مقادیر بالای ناحیه تکین بر روی نمودار رنگی جلوگیری کرده و به شما اجازه می‌دهد تا بر روی نتایج مهم تمرکز کنید.

روش دوم: اصلاح هندسه و شرایط مرزی (رویکرد پیشنهادی)

این روش، پایدارترین و دقیق‌ترین راهکار برای حذف سینگولاریتی است، زیرا علت اصلی تکینگی—یعنی ایده‌آل‌سازی فیزیکی—را برطرف می‌کند و به جای مدیریت علائم، یک مدل اساساً دقیق‌تر ایجاد می‌کند.

• اضافه کردن فیلت (Fillet): جایگزین کردن گوشه‌های تیز با فیلت‌هایی با شعاع مناسب در مدل CAD، سینگولاریتی را به یک تمرکز تنش واقعی و همگراشونده تبدیل می‌کند. این کار ممکن است نیازمند مش‌بندی ریزتر در ناحیه فیلت برای دستیابی به نتایج دقیق باشد.
• توزیع بارها و قیود: به جای اعمال نیرو یا قید بر روی یک رأس یا گره، آن را بر روی یک سطح یا لبه توزیع کنید (مثلاً به صورت فشار یا Force روی یک سطح). این کار باعث می‌شود با ریز شدن مش، مساحت اعمال بار ثابت بماند و از ایجاد تنش بی‌نهایت جلوگیری شود.
• واقع‌گرایانه‌تر کردن قیود: به جای استفاده از قیود ایده‌آل مانند Fixed Support، از ترکیب قیود واقع‌گرایانه‌تر استفاده کنید. برای مثال، در مدل قلاب که توسط یک پیچ به دیوار متصل شده است، می‌توان قید Fixed Support را با ترکیبی از موارد زیر جایگزین کرد:
  • Cylindrical Support برای شبیه‌سازی پین یا پیچ (محدودیت حرکت شعاعی اما آزاد در حرکت مماسی و محوری).
  • Elastic Support برای شبیه‌سازی تماس فشاری قلاب با سطح دیوار. این قید بر اساس سختی فونداسیون (foundation stiffness) عمل می‌کند که به عنوان «فشار مورد نیاز برای ایجاد یک واحد تغییر شکل عمودی در فونداسیون» تعریف می‌شود و یک بازنمایی بسیار واقعی‌تر از تماس با یک دیوار تغییرشکل‌پذیر است.

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

اگر تغییر هندسه ممکن یا مطلوب نیست، استفاده از یک مدل ماده غیرخطی (مانند elastic, nearly perfectly plastic) می‌تواند اثر سینگولاریتی را محدود کند. در این حالت، المان‌های ناحیه تکینگی به صورت موضعی تسلیم شده (yield) و تنش در آن ناحیه به حد تنش تسلیم ماده محدود می‌شود. بار اضافی به نواحی الاستیک اطراف منتقل شده و از افزایش بی‌رویه تنش جلوگیری می‌شود. این روش از نظر محاسباتی سنگین‌تر است اما می‌تواند درک بهتری از رفتار موضعی ماده در نواحی با تنش بالا ارائه دهد.

روش چهارم: استفاده از تکنیک‌های پس‌پردازش (Stress Linearization)

تکنیک خطی‌سازی تنش (Stress Linearization) یک ابزار ارزیابی قدرتمند در پس‌پردازش است که میدان تنش را در یک مسیر مشخص (معمولاً در راستای ضخامت) به مؤلفه‌های غشایی (membrane) و خمشی (bending) تجزیه می‌کند. این تکنیک تکینگی را حذف نمی‌کند، اما در حضور آن بسیار مقاوم عمل می‌کند. دلیل این امر آن است که مؤلفه‌های غشایی و خمشی از طریق انتگرال‌گیری تنش‌ها در طول مسیر محاسبه می‌شوند. این فرآیند انتگرال‌گیری، اثر قله تنش بسیار بالا و موضعی در نقطه تکین را усреднён (average out) کرده و مقداری پایدار و همگرا که نماینده وضعیت تنش کلی است، ارائه می‌دهد. این تکنیک به طور گسترده در استانداردهای طراحی مانند ASME Boiler and Pressure Vessel Code برای ارزیابی تنش‌ها استفاده می‌شود.

با شناخت این رویکردها، می‌توان به یک جمع‌بندی جامع برای مدیریت سینگولاریتی در تحلیل‌های مهندسی رسید.

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

در این مقاله آموختیم که سینگولاریتی تنش، یک آرتیفکت ریاضی ناشی از ساده‌سازی مدل است و نه یک خطای تحلیلی یا نقص فیزیکی در قطعه. این پدیده اثری کاملاً محلی دارد و با روش‌های مشخصی قابل شناسایی و مدیریت است. ما دلایل اصلی بروز آن—از جمله گوشه‌های تیز، بارهای نقطه‌ای و قیود نامناسب—را بررسی کرده و چهار رویکرد عملی برای مقابله با آن ارائه دادیم.

توصیه نهایی و کلیدی این است: بهترین رویکرد برای رفع سینگولاریتی در انسیس به ماهیت مدل، اهداف تحلیل و موقعیت ناحیه تکین بستگی دارد. یک تحلیل‌گر ماهر باید با قضاوت مهندسی تصمیم بگیرد که آیا باید سینگولاریتی را با اصلاح مدل حذف کند (که رویکرد دقیق‌تری است) یا می‌تواند با تفسیر صحیح نتایج و تمرکز بر نواحی کلیدی، آن را نادیده بگیرد. در نهایت، تسلط بر مدیریت تکینگی‌ها یکی از نشانه‌های یک تحلیل‌گر خبره است که این چالش را از یک خطای عددی به فرصتی برای ساخت شبیه‌سازی‌های قوی‌تر و واقع‌گرایانه‌تر برای اتخاذ تصمیمات مهندسی مطمئن تبدیل می‌کند.