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

۱. مقدمه

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

جرقه‌های اصطکاکی، پدیده‌ای چندفیزیکی است که در آن مکانیک، ترمودینامیک، آیرودینامیک و انتقال حرارت به‌طور هم‌زمان نقش‌آفرینی می‌کنند. از یک سو، انرژی مکانیکی اصطکاک به انرژی حرارتی تبدیل می‌شود و ذرات داغ از سطح جدا می‌شوند؛ از سوی دیگر، این ذرات در هوا پرتاب شده، تحت تأثیر نیروی پسا قرار گرفته و هم‌زمان با محیط تبادل حرارت می‌کنند. افزون بر این، طبیعت تصادفی توزیع اندازه، سرعت و زاویهٔ پرتاب ذرات، به این پدیده بُعدی آماری و آشوب‌گونه می‌بخشد.

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

۲. مبانی فیزیکی و ریاضی

۲.۱. توان اصطکاکی

منبع اصلی انرژی جرقه، اصطکاک میان ابزار و قطعه کار است. توان اصطکاکی تولیدشده از رابطهٔ زیر به دست می‌آید:

که در آن $\mu$ ضریب اصطکاک، $F_N$ نیروی عمودی، $D$ قطر ابزار و $N$ دور بر دقیقه است. با مقادیر پیش‌فرض شبیه‌ساز — $D = 10$ mm، $N = 12000$ rpm، $\mu = 0.4$ و $F_N = 120$ N — سرعت خطی ابزار $v_{\text{tan}} = 6.28$ m/s و توان اصطکاکی $P_{\text{fric}} = 301.6$ W محاسبه می‌شود.

۲.۲. انرژی ذوب ذره

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

که در آن $\rho_p$ چگالی ذره، $d$ قطر ذره، $c_p$ گرمای ویژه، $T_m$ دمای ذوب، $T_{\text{amb}}$ دمای محیط و $L_f$ گرمای نهان ذوب است. برای مادهٔ فوق‌سخت کاربید تنگستن با پارامترهای شبیه‌ساز، $E_p \approx 8.17 \times 10^{-3}$ J است.

۲.۳. نرخ تولید ذرات

نرخ تولید ذرات جرقه با تقسیم توان اصطکاکی مؤثر بر انرژی ذوب هر ذره به دست می‌آید:

که در آن $\eta$ بازده تبدیل انرژی مکانیکی به حرارتی است. با مقادیر پیش‌فرض، $\dot{N} \approx 25{,}840$ ذره بر ثانیه برآورد می‌شود که با مشاهدات تجربی در ماشین‌کاری مواد فوق‌سخت سازگار است.

۲.۴. معادلهٔ حرکت ذره و نیروی پسا

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

نیروی پسا از همبستگی تجربی شیلر-نومن محاسبه می‌شود:

که در آن $Re = \rho_g d \|\mathbf{v}_r\| / \mu_g$ عدد رینولدز ذره است و برای $Re < 1000$ معتبر است.

۲.۵. انتقال حرارت و خنک‌سازی ذره

هم‌زمان با حرکت، هر ذره با محیط تبادل حرارت می‌کند. با فرض مدل ظرفیت فشرده (عدد بایوت کوچک)، معادلهٔ انتقال حرارت به‌صورت زیر نوشته می‌شود:

ضریب انتقال حرارت جابجایی از همبستگی رنز-مارشال محاسبه می‌شود:

با ترکیب این روابط، ثابت زمانی خنک‌سازی به‌صورت زیر حاصل می‌شود:

نکتهٔ کلیدی آن است که ثابت زمانی خنک‌سازی با توان ۳/۲ قطر ذره متناسب است؛ به این معنا که ذرات بزرگ‌تر بسیار کندتر از ذرات کوچک‌تر خنک می‌شوند و در نتیجه، طول عمر بیشتری دارند.

۲.۶. نمای لیاپانوف: معیار آشوب

برای بررسی رفتار آشوب‌گونهٔ سیستم، از مفهوم نمای لیاپانوف استفاده می‌شود. این کمیت، نرخ واگرایی نمایی دو مسیر مجاور در فضای فاز را اندازه می‌گیرد:

مقدار مثبت نمای لیاپانوف، نشانهٔ آشوب‌ناکی سیستم و حساسیت ذاتی آن به شرایط اولیه است — همان «اثر پروانه‌ای» معروف. در شبیه‌ساز حاضر، این کمیت با اجرای دو مجموعهٔ شبیه‌سازی موازی با اغتشاش ۰٫۱ درصدی در نیروی عمودی و دور، و ردیابی فاصلهٔ مسیر ذرات متناظر محاسبه می‌شود.

۳. شبیه‌سازی عددی

۳.۱. معماری سامانه شبیه‌سازی

شبیه‌ساز حاضر بر پایه معماری سمت‌کلاینت در مرورگر طراحی شده است و از سه لایهٔ اصلی تشکیل می‌شود: لایهٔ ورودی (۱۶ پارامتر قابل تنظیم)، لایهٔ موتور محاسباتی (روش مونت‌کارلو با قابلیت Seed)، و لایهٔ نمایش نتایج (کارت‌های آماری، چهار نمودار، جدول و صحنهٔ سه‌بعدی تعاملی).

۳.۲. روش مونت‌کارلو و قابلیت تکرارپذیری

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

برای دستیابی به نتایج تکرارپذیر، شبیه‌ساز از یک مولد اعداد تصادفی با دانهٔ ثابت (Seed) پشتیبانی می‌کند. با فعال بودن این گزینه، هر بار اجرای شبیه‌سازی، دقیقاً همان اعداد تصادفی و در نتیجه، همان نتایج تولید می‌شوند. این ویژگی، برای اعتبارسنجی و مقایسهٔ علمی ضروری است.

// مولد اعداد تصادفی با دانهٔ ثابت — برای تکرارپذیری
class SeededRandom {
  constructor(seed) { this.seed = seed; }
  next() {
    this.seed = (this.seed * 9301 + 49297) % 233280;
    return this.seed / 233280;
  }
  random() { return this.next(); }
}

// نمونه‌گیری از توزیع لگ‌نرمال (برای قطر ذرات)
function logNormalSample(mu, sigma, rand) {
  let u = 0, v = 0;
  while (u === 0) u = rand();
  while (v === 0) v = rand();
  const z = Math.sqrt(-2 * Math.log(u)) * Math.cos(2 * Math.PI * v);
  return Math.exp(mu + sigma * z);
}

۳.۳. محاسبهٔ نمای لیاپانوف در عمل

برای محاسبهٔ نمای لیاپانوف، شبیه‌ساز دو نسخهٔ موازی از سیستم را با اغتشاش کوچک $\delta = 0.001$ در نیروی عمودی و دور اجرا می‌کند. برای هر ذره، مسیر حرکت در دو نسخه ردیابی می‌شود و فاصلهٔ بین دو مسیر در زمان $t$ اندازه‌گیری می‌گردد. نمای لیاپانوف از میانگین‌گیری روی ذرات محاسبه می‌شود.

// شبیه‌سازی موازی دو نسخه با اغتشاش کوچک
const perturb = 0.001;  // 0.1% اغتشاش
const FN2 = FN * (1 + perturb);
const N2 = N * (1 + perturb);

let sumLog = 0, count = 0;
const dt = 0.005, maxTime = 0.4;

for (let i = 0; i < n; i++) {
  // ردیابی همزمان دو مسیر
  let t = 0, dist0 = 1e-12, found = false;
  while (t < maxTime) {
    // به‌روزرسانی سرعت و موقعیت هر دو مسیر
    // ... (معادلات حرکت با نیروی پسا)
    t += dt;

    const dx = x2 - x1, dy = y2 - y1, dz = z2 - z1;
    const dist = Math.sqrt(dx**2 + dy**2 + dz**2);

    if (dist > 1e-10 && !found) {
      dist0 = dist; found = true;
    }
    if (found && dist > 1e-10) {
      sumLog += Math.log(dist / dist0);
      count++;
      break;
    }
  }
}

let lambda = sumLog / count / (maxTime / 2);

۴. نتایج و تحلیل

۴.۱. شاخص‌های کلیدی

با اجرای شبیه‌سازی با پارامترهای پیش‌فرض، شاخص‌های کلیدی زیر به دست می‌آید:

۴.۲. توزیع اندازه ذرات

توزیع اندازهٔ ذرات جرقه، الگویی نزدیک به توزیع لگ‌نرمال نشان می‌دهد که با مشاهدات تجربی از توزیع اندازهٔ ذرات در فرآیندهای مکانیکی سازگار است. بیشتر ذرات در محدودهٔ ۵۰ تا ۲۰۰ میکرومتر قرار می‌گیرند، اما دنبالهٔ توزیع تا چند صد میکرومتر نیز امتداد می‌یابد.

۴.۳. توزیع طول عمر ذرات

طول عمر هر ذره — یعنی زمانی که طول می‌کشد تا دمای آن از دمای تابش مرئی (۱۰۷۳ کلوین) به دمای محیط کاهش یابد — به قطر ذره بستگی دارد. با توجه به رابطهٔ $\tau \propto d^{3/2}$، ذرات بزرگ‌تر بسیار کندتر خنک می‌شوند و در نتیجه، طول عمر بیشتری دارند.

۴.۴. توزیع زاویه پرتاب

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

۴.۵. واگرایی آشوب‌گونه

حساسیت به شرایط اولیه، هستهٔ رفتار آشوب‌گونهٔ سیستم است. با اعمال اغتشاش ۰٫۱ درصدی در نیروی عمودی و دور، مسیر ذرات به‌طور نمایی از مسیر مرجع واگرا می‌شود. شیب این واگرایی نمایی، نمای لیاپانوف را تعیین می‌کند.

۴.۶. جدول آماری جامع

۵. مواد فوق‌سخت و کاربردهای صنعتی

۵.۱. کاربید تنگستن (WC)

کاربید تنگستن، پرکاربردترین مادهٔ فوق‌سخت در ابزارهای برش و سنگ‌زنی است. چگالی آن حدود $15000$ kg/m³ و دمای ذوب آن نزدیک به $1673$ K است. این پارامترها در شبیه‌ساز به‌عنوان مقادیر پیش‌فرض انتخاب شده‌اند و با داده‌های تجربی مواد تجاری سازگاری خوبی دارند.

۵.۲. کاربید بور (B₄C)

کاربید بور، سخت‌ترین مادهٔ تجاری پس از الماس و نیترید بور مکعبی است. چگالی کم آن (حدود $2520$ kg/m³) و دمای ذوب بالا (نزدیک به $2623$ K) آن را برای کاربردهای زرهی و هوافضا مناسب می‌کند. با تنظیم پارامترهای شبیه‌ساز، می‌توان رفتار جرقه در ماشین‌کاری این ماده را نیز پیش‌بینی کرد.

۵.۳. الماس مصنوعی و CBN

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

۶. نتیجه‌گیری

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

نتایج نشان داد که سیستم از رفتار آشوب‌گونهٔ روشنی برخوردار است: با اعمال اغتشاش ۰٫۱ درصدی در پارامترهای ورودی، مسیر ذرات به‌طور نمایی واگرا می‌شود و نمای لیاپانوف مثبت به‌دست می‌آید. این نتیجه، بر حساسیت ذاتی پدیدهٔ جرقهٔ اصطکاکی به شرایط اولیه تأکید می‌کند و نشان می‌دهد که پیش‌بینی دقیق بلندمدت رفتار این سیستم، اصولاً محدود است.

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

پرسش‌های باقی‌مانده برای پژوهش‌های آینده عبارتند از: چگونه می‌توان اثرات تشعشع حرارتی و برخورد ذرات را در مدل گنجاند؟ آیا می‌توان از الگوریتم‌های یادگیری ماشین برای پیش‌بینی رفتار آشوب‌گونهٔ سیستم بهره گرفت؟ و چگونه می‌توان نتایج شبیه‌ساز را با اندازه‌گیری‌های تجربی دقیق — نظیر طیف‌سنجی نوری جرقه — اعتبارسنجی کرد؟