بلاگ
در این مقاله، ضمن پوشش مبانی ریاضی و آماری هر ابزار، به صورت عملی نحوه پیادهسازی در Minitab، تحلیل نمودارها و اقدامات اصلاحی را فرا خواهید گرفت. تمامی مطالب بر اساس مراجع معتبر AIAG (ویرایش دوم) و VDA Volume 4 & 5 تدوین شده است.
در صنایع گسسته (مثل قطعات فلزی)، فرض بر استقلال مشاهدات و توزیع نرمال است. اما در فرآیندهای شیمیایی (رنگ، فوم، پلیمر) پدیدههایی مانند خودهمبستگی (Autocorrelation) و چندمتغیره بودن (Multivariate) نقضکنندهی این فروض هستند. بر اساس قضیه حد مرکزی و تعمیمهای آن، روشهایی مانند میانگینهای متحرک وزنی (EWMA) و مجموع تجمعی (CUSUM) برای غلبه بر این چالشها توسعه یافتهاند (Montgomery, 2009).
AIAG VDA 4 Minitab
نمودار EWMA بر اساس تئوری هموارسازی نمایی (Exponential Smoothing) بنا شده است. آماره Zt ترکیب خطی از مشاهدات فعلی و گذشته است. واریانس این آماره به صورت همگرا به σ²[λ/(2-λ)] میل میکند. انتخاب λ بهینه با توجه به اندازه شیفت مورد نظر تعیین میشود (مقدار کوچک λ برای تشخیص شیفتهای کوچک). مبنای تئوریک این روش به تئوری تشخیص تغییر (Change Point Detection) و فیلتر کالمن نیز مرتبط است.
ویسکوزیته هدف = ۱۰۰ واحد، σ=۲، λ=۰.۲. دادههای ۱۰ بچ متوالی: ۱۰۲, ۱۰۱, ۹۹, ۹۸, ۹۷, ۹۶, ۹۸, ۱۰۰, ۹۹, ۹۷.
مسیر در Minitab: Stat → Control Charts → Time-Weighted Charts → EWMA
ورودی: دادهها در یک ستون. در پنجره EWMA، گزینه "Parameters" را انتخاب و میانگین هدف (۱۰۰) و انحراف معیار (۲) را وارد کنید. λ = ۰.۲ و L = ۳.
خروجی: Minitab حدود کنترل را محاسبه میکند: UCL ≈ ۱۰۰.۹ , LCL ≈ ۹۹.۱. نقاط ۵ (۹۷) و ۶ (۹۶) خارج از حد پایین هستند.
AIAG VDA Minitab
نمودار CUSUM ریشه در آزمون نسبت درستنمایی متوالی (SPRT) دارد که توسط پیج (Page) در دهه ۱۹۵۰ فرموله شد. ایده اصلی: جمعآوری انحرافات جزئی از هدف در طول زمان. اگر فرآیند در حالت کنترل باشد، مجموع انحرافات حول صفر نوسان میکند. اگر میانگین فرآیند به اندازه Δ جابجا شود، آماره CUSUM به صورت خطی رشد میکند. پارامتر K (مجاز) معمولاً نصف بزرگی شیفتی که میخواهیم سریع تشخیص دهیم انتخاب میشود (K = Δ/2). مقدار آستانه H (تصمیم) معمولاً ۴ یا ۵ برابر انحراف معیار فرآیند در نظر گرفته میشود.
دانسیته هدف فوم = ۳۰ kg/m³، σ=۱. دادههای ۱۰ بچ: ۳۰.۲, ۳۰.۵, ۳۰.۹, ۳۱.۲, ۳۱.۶, ۳۱.۸, ۳۲.۱, ۳۲.۰, ۳۱.۷, ۳۲.۲.
مسیر در Minitab: Stat → Control Charts → Time-Weighted Charts → CUSUM
ورودی: در پنجره CUSUM، گزینه "Parameters" را انتخاب و μ₀=۳۰, σ=۱ وارد کنید. معمولاً Minitab K=۰.۵ و H=۵ را پیشنهاد میدهد.
خروجی: از بچ چهارم به بعد، آماره C⁺ از H=۵ عبور میکند و فرآیند خارج از کنترل اعلام میشود.
VDA (Q-DAS) Minitab
آماره T² تعمیم چندمتغیره آماره t تکمتغیره است. هتلینگ در دهه ۱۹۳۰ این آماره را بر اساس توزیع نرمال چندمتغیره معرفی کرد. آماره T² فاصله ماهالانوبیس بین میانگین نمونه و میانگین هدف را اندازهگیری میکند. حدود کنترل با استفاده از توزیع F به دست میآید. این آماره همبستگی بین متغیرها را در ماتریس کوواریانس لحاظ میکند و به تغییرات در ساختار همبستگی حساس است.
سه ویژگی مهم: ویسکوزیته (هدف ۱۵۰۰ cP)، اسیدیته (هدف ۵ mgKOH/g)، و شفافیت (هدف ۹۵%). دادههای ۱۰ بچ در سه ستون.
مسیر در Minitab: Stat → Control Charts → Multivariate Charts → T²
ورودی: هر سه ستون را انتخاب کنید. اگر زیرگروه ندارید، از "Single column" استفاده کنید.
خروجی: Minitab آماره T² را برای هر بچ رسم کرده و حد بالای کنترل (مثلاً UCL=۱۳.۵) را نمایش میدهد. بچ ۷ T²=۱۸ دارد (خارج از کنترل).
AIAG Short Run Minitab
در شرایطی که محصولات مختلف با اهداف متفاوت (Ti) داریم، نمیتوان دادههای خام را مستقیماً مقایسه کرد. روش DNOM (Deviation from Nominal) بر پایه استانداردسازی (Standardization) استوار است. با تعریف متغیر جدید Y = X - T، موقعیت هر محصول را نسبت به هدف خود میسنجیم. اگر فرض کنیم برای همه محصولات، پراکندگی حول هدف یکسان است (یا با استفاده از انحراف معیار تلفیقی Pooled)، آنگاه میتوان یک نمودار کنترلی واحد با حدود ۳σ برای Y رسم کرد.
سه نوع چسب با گرانروی هدف متفاوت: ۵۰۰, ۱۰۰۰, ۲۰۰۰ cP. دادهها در دو ستون: ستون اول گرانروی واقعی، ستون دوم مقدار هدف.
مسیر در Minitab: Stat → Control Charts → Variables Charts for Subgroups → DNOM
ورودی: در DNOM Chart، "Observations" را گرانروی واقعی و "Target" را ستون اهداف وارد کنید.
خروجی: یک نقطه DNOM=+۱۲۰ (هدف ۱۰۰۰، مقدار واقعی ۱۱۲۰) خارج از حدود (مثلاً ۳σ≈۹۰).
VDA Minitab (ماکرو/PLS)
در بسیاری از فرآیندهای شیمیایی (مثل پخت رزین، واکنش پلیمریزاسیون)، کیفیت نهایی به شکل منحنی دما-زمان یا فشار-زمان وابسته است. پروفایل مانیتورینگ از آنالیز رگرسیون تابعی (Functional Regression) و آماره T² چندمتغیره برای منحنیها استفاده میکند. ابتدا یک منحنی مرجع (Reference Profile) با استفاده از دادههای تاریخی برازش داده میشود. سپس برای هر منحنی جدید، بردار انحرافات محاسبه و آماره T² حاصل با حد کنترل مقایسه میشود.
منحنی استاندارد دما-زمان: t=1:25°C, t=5:80°C, t=10:120°C. یک بچ جدید: t=5:60°C (انحراف شدید).
مسیر در Minitab: ابزار دقیق نیاز به ماکرو نویسی یا استفاده از نرمافزار Q-DAS دارد. با ثبت نقاط کلیدی و محاسبه T² میتوان انحراف را سنجید.
هشدار: انحراف از منحنی مرجع نشاندهنده کمفعالیتی کاتالیست یا افت دما است.
در فرآیندهای بچ شیمیایی (مانند پلیمریزاسیون، پخت رزین، تخمیر)، کیفیت نهایی به شدت به شکل منحنی دما-زمان، فشار-زمان یا سایر پروفایلهای فرآیندی وابسته است. پروفایل مانیتورینگ تکنیکی پیشرفته برای کنترل آماری این منحنیها است. در این بخش ضمن تشریح مبانی نظری، با یک مثال عملی در نرمافزار Q-DAS مراحل پیادهسازی را گام به گام یاد میگیرید.
روشهای متداول برای پروفایل مانیتورینگ عبارتند از:
نرمافزار Q-DAS از روشهای چندمتغیره پیشرفته (مانند تحلیل مؤلفههای اصلی تابعی) برای پروفایل مانیتورینگ استفاده میکند. خروجی اصلی شامل نمودار T² و نمودار Q (SPE) است.
فرض کنید در یک فرآیند تولید رزین، دمای راکتور در ۱۰ نقطه زمانی (هر ۵ دقیقه) ثبت میشود. پروفایل استاندارد (مرجع) از ۲۰ بچ خوب (تحت کنترل) به دست آمده است. اکنون یک بچ جدید (بچ شماره ۲۱) تولید شده و میخواهیم انحراف پروفایل آن را بررسی کنیم.
| بچ | t=5 | t=10 | t=15 | t=20 | t=25 | t=30 | t=35 | t=40 | t=45 | t=50 |
|---|---|---|---|---|---|---|---|---|---|---|
| ۱ | ۵۰ | ۶۲ | ۷۵ | ۸۵ | ۹۲ | ۹۷ | ۱۰۰ | ۱۰۱ | ۱۰۲ | ۱۰۲ |
| ۲ | ۴۹ | ۶۱ | ۷۴ | ۸۴ | ۹۱ | ۹۶ | ۹۹ | ۱۰۰ | ۱۰۱ | ۱۰۱ |
| ۳ | ۵۱ | ۶۳ | ۷۶ | ۸۶ | ۹۳ | ۹۸ | ۱۰۱ | ۱۰۲ | ۱۰۳ | ۱۰۳ |
| ۴ | ۵۰ | ۶۲ | ۷۵ | ۸۵ | ۹۲ | ۹۷ | ۱۰۰ | ۱۰۱ | ۱۰۲ | ۱۰۲ |
| ۵ | ۵۲ | ۶۴ | ۷۷ | ۸۷ | ۹۴ | ۹۹ | ۱۰۲ | ۱۰۳ | ۱۰۴ | ۱۰۴ |
بچ جدید (۲۱): ۴۸, ۵۸, ۶۸, ۷۷, ۸۵, ۹۱, ۹۶, ۹۹, ۱۰۰, ۱۰۰
نرمافزار Q-DAS (ماژول CAMERA یا Profile Analysis) امکان پروفایل مانیتورینگ را به صورت یکپارچه فراهم میکند. مراحل زیر را دنبال کنید:
دادهها را در یک فایل متنی یا Excel با ساختار زیر ذخیره کنید: هر ردیف = یک بچ، ستون اول شماره بچ (اختیاری)، ستونهای بعدی مقادیر دما در زمانهای متوالی.
File → Import → ASCII Import → انتخاب فایل
در ویزارد واردات، نوع داده را "Variables" و فرمت ستونها را "Double" انتخاب کنید. زمانها را به عنوان "Characteristic" تعریف کنید.
Analyze → Profile Analysis → Define Profile
Analyze → Profile Analysis → Run Analysis
Q-DAS ابتدا یک مدل PCA بر روی دادههای مرجع برازش میدهد. سپس برای هر بچ (از جمله بچ جدید) آمارههای T² و Q محاسبه میشود.
نمودارهای کنترلی T² و Q نمایش داده میشوند. در مثال ما، برای بچ ۲۱:
Q-DAS همچنین نمودار "Contribution Plot" را نمایش میدهد که نشان میدهد کدام نقاط زمانی بیشترین سهم را در انحراف داشتهاند.
شامل مبانی نظری، کار عملی با Minitab و Q-DAS، مطالعات موردی در صنایع رنگ، فوم، رزین، چسب و دارو
با تدریس تیم مشاوران عصر کیفیت (بیش از ۱۵ سال تجربه صنعتی)
📞 09125076715مشاوره و ثبتنام رایگان - ظرفیت محدود
منابع اصلی: Montgomery, D.C. (2009). Introduction to Statistical Quality Control (ویرایش هفتم)، Box, G. & Luceno, A. (1997). Statistical Control by Monitoring and Feedback Adjustment، استانداردهای AIAG SPC ویرایش دوم، VDA Volume 4 و 5، و مستندات نرمافزار Q-DAS (۲۰۲۴).
© ۲۰۲۵ - کلیه حقوق مادی و معنوی این محتوا متعلق به شرکت عصر کیفیت میباشد.
در صنایع فلزی، قطعات پشت سر هم تولید میشوند و ما به راحتی میتوانیم ۵ قطعه متوالی را به عنوان یک زیرگروه (Subgroup) انتخاب کنیم. اما در صنایع شیمیایی:
فرآیندهای پیوسته یا بطهای (Batch): مواد شیمیایی معمولاً در تانکهای بزرگ و به صورت پیوسته یا بچ تولید میشوند. مفهوم “قطعه پنجم” معنا ندارد.
خودهمبستگی (Autocorrelation): در فرآیندهای شیمیایی، مقدار فعلی (مثلاً ویسکوزیته) به شدت به مقدار لحظه قبل وابسته است. این ویژگی باعث میشود نمودارهای سنتی خطاهای زیادی نشان دهند.
تنظیم مکرر (Frequent Adjustments): اپراتورها مدام دما یا فشار را تنظیم میکنند. این کار باعث میشود فرآیند هرگز در حالت ایستا (Stable) قرار نگیرد که شرط اصلی استفاده از نمودارهای استاندارد است .
خوشبختانه، هم مرجع AIAG (ویرایش دوم) و هم مراجع VDA روشهای مدرنتری را معرفی میکنند. در اینجا چند مورد کلیدی را با ذکر مثال توضیح میدهم:
در فرآیندهای بچ، شما یک منحنی رشد دارید. برای مثال، در فرآیند هیدروژناسیون روغن برای تولید مارگارین، نقطه ذوب در طول زمان افزایش مییابد. بجای اینکه هر بار یک نقطه را چک کنیم، کل منحنی را مانیتور میکنیم .
مثال عملی (صنعت فوم):
فرض کنید در تولید فوم پلییورتان، دما در طول ۱۰ دقیقه واکنش باید از ۲۵ درجه به ۱۲۰ درجه برسد. شما یک منحنی استاندارد دارید.
ابزار پیشنهادی: نمودار T² (هتلینگ) یا آماره AUGC.
کاربرد: اگر در یک بچ، شیب افزایش دما کندتر از حالت استاندارد باشد، سیستم هشدار میدهد که شاید کاتالیست به اندازه کافی موثر نیست .
در فرآیندهای پیوسته مانند تولید رنگ، شما نمیتوانید منتظر بمانید تا ۲۰ نمونه ویسکوزیته بالا بیاید تا بگویید فرآیند از کنترل خارج شده است. شما باید همزمان با کنترل فرآیند (توسط سیستمهای کنترلی)، آمار هم بگیرید .
مثال عملی (صنعت رنگ):
در راکتور پلیمریزاسیون الکید (ماده پایه رنگ)، ویسکوزیته تحت تاثیر دما و فشار است.
روش پیشنهادی: Run-to-Run (RTR) Control.
کاربرد: شما بین دو بچ، بر اساس دادههای بچ قبلی، پارامترهای بچ بعدی (مثلاً دمای پخت) را تنظیم میکنید، اما همزمان با نمودارهای SPC ساده (مثل Moving Range) بررسی میکنید که آیا این تنظیمات خود باعث ایجاد تغییر خاصی (Special Cause) نشده است.
رنگ دارای چند مشخصه کلیدی است: ویسکوزیته، وزن مخصوص، و شفافیت. اینها به هم وابسته هستند. اگر فقط ویسکوزیته را چک کنید، ممکن است متوجه نشوید که فرآیند در آستانه مشکل است .
ابزار پیشنهادی: نمودار χ² یا T² چندمتغیره.
در فرآیندهای شیمیایی، ممکن است کیفیت کاتالیست به تدریج کاهش یابد. این تغییر کوچک است، اما اگر زود تشخیص داده نشود، کل بچ خراب میشود .
مثال عملی (صنعت چسب):
شما در حال تولید چسب هستید و چسبندگی هر ساعت یکبار تست میشود.
نمودار EWMA: این نمودار به دادههای قدیمی وزن کمتری میدهد و به دادههای جدید وزن بیشتری میدهد. به این ترتیب، به محض اینکه چسبندگی شروع به کاهش جزئی کند، نمودار زنگ خطر را به صدا در میآورد.
اگر شرکت شما تنوع محصول بالا و تعداد تولید کم دارد (مثلاً هر بار یک تانکر رنگ متفاوت)، نمیتوانید برای هر محصول یک نمودار جداگانه بکشید .
ابزار پیشنهادی: نمودار Z-MR (Z-Moving Range) یا نمودار انحراف از هدف (DNOM).
کاربرد: شما مقدار واقعی را منهای مقدار هدف (Target) میکنید. به این ترتیب، تمام محصولات مختلف را روی یک نمودار میآورید و تغییرات فرآیند را میبینید.