آشنایی با میرایی (Damping) در تحلیل‌های دینامیکی: میرایی رایلی (Rayleigh Damping)

میرایی در تحلیل دینامیکی چیست؟

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

تعریف مفهومی میرایی

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

1. میرایی ویسکوز (Viscous Damping): در این مدل، نیروی میرایی متناسب با سرعت نسبی اجزای سیستم در نظر گرفته می‌شود.
2. میرایی سازه‌ای (Structural Damping): در این مدل، استهلاک انرژی متناسب با جابجایی نسبی اجزا و مستقل از فرکانس در نظر گرفته می‌شود.

جایگاه میرایی رایلی

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

مبانی ریاضی میرایی رایلی

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

فرمول اصلی

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

[C] = α[M] + β[K]

اجزای این معادله عبارتند از:

• [C]: ماتریس میرایی سیستم
• [M]: ماتریس جرم سیستم
• [K]: ماتریس سختی سیستم
• α (آلفا): ضریب میرایی رایلی وابسته به جرم (Mass Proportional Damping Coefficient)
• β (بتا): ضریب میرایی رایلی وابسته به سختی (Stiffness Proportional Damping Coefficient)

مزیت کلیدی

مزیت اصلی این فرمولاسیون، سادگی و کارایی آن در تحلیل‌های مودال است. این مزیت از ویژگی زیر نشأت می‌گیرد:

• حفظ خاصیت تعامد (Orthogonality): مدل رایلی تضمین می‌کند که ماتریس میرایی حاصل ([C])، شرط تعامد را نسبت به شکل مودهای ارتعاشی سیستم ارضا می‌کند. این ویژگی به لحاظ عملی به این معناست که می‌توان یک سیستم پیچیده و کوپل شده با چند درجه آزادی (MDOF) را به مجموعه‌ای از سیستم‌های ساده و مستقل تک درجه آزادی (SDOF) تفکیک و حل کرد. این قابلیت، دلیل اصلی کارایی محاسباتی بالای این مدل در تحلیل‌های مبتنی بر مودال است.

این فرمولاسیون ریاضی، منجر به یک رفتار وابسته به فرکانس می‌شود که درک آن برای استفاده صحیح از این مدل ضروری است.

رفتار فرکانس-وابسته میرایی رایلی

درک وابستگی میرایی رایلی به فرکانس برای انتخاب صحیح ضرایب و تفسیر نتایج تحلیل دینامیکی حیاتی است. برخلاف تصور اولیه، این مدل یک نسبت میرایی ثابت برای تمام فرکانس‌های طبیعی سیستم ارائه نمی‌دهد، بلکه رفتار آن در فرکانس‌های مختلف متغیر است.

فرمول نسبت میرایی

رابطه نسبت میرایی (ξᵢ) برای مود ارتعاشی i-ام با فرکانس طبیعی دایره‌ای (ωᵢ) از معادله اصلی رایلی مشتق شده و به صورت زیر است:

ξᵢ = (α / 2ωᵢ) + (βωᵢ / 2)

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

تحلیل رفتار در فرکانس‌های مختلف

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

• فرکانس‌های پایین: در فرکانس‌های بسیار پایین، جمله وابسته به جرم (α / 2ωᵢ) غالب است و میرایی بسیار زیادی ایجاد می‌کند. این پدیده می‌تواند به طور مصنوعی حرکات صلب جسم (Rigid Body Motion) را در مدل‌هایی که به طور کامل مقید نیستند، سرکوب یا «قفل» کند و به نتایج نادرست منجر شود.
• فرکانس‌های بالا: در فرکانس‌های بالا، جمله وابسته به سختی (βωᵢ / 2) غالب می‌شود و میرایی به شدت افزایش می‌یابد. این امر می‌تواند مودهای فرکانس-بالای سیستم را بیش از حد مستهلک کند؛ این یک ساده‌سازی غیرمحافظه‌کارانه است که ممکن است پاسخ‌های فیزیکی واقعی را پنهان کرده و منجر به تخمین کمتر از واقعیت برای تنش‌ها و شتاب‌ها شود.
• فرکانس‌های میانی: یک محدوده فرکانسی وجود دارد که در آن، میرایی به مقدار حداقل خود می‌رسد و تغییرات آن نسبتاً کم است. مهندسان معمولاً ضرایب α و β را به گونه‌ای انتخاب می‌کنند که این محدوده، فرکانس‌های غالب و مهم سیستم را پوشش دهد.

جدول خلاصه رفتار فرکانسی

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

چالش اصلی در استفاده از این مدل، تعیین مقادیر مناسب برای ضرایب α و β است تا رفتار میرایی در محدوده فرکانس‌های مورد نظر، واقع‌بینانه باشد. این موضوع در بخش بعدی بررسی می‌شود.

کاربرد عملی: تعیین ضرایب رایلی در انسیس

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

چالش اصلی

چالش اصلی در عمل، تعیین مقادیر معنادار برای α و β است. تخمین مستقیم و شهودی مقادیر مناسب برای α و β در سیستم‌های پیچیده با درجات آزادی بالا، غیرعملی است. بنابراین، یک رویکرد نظام‌مند و مهندسی برای تعریف آن‌ها بر اساس میرایی مودال مطلوب، ضروری است.

رویکرد مهندسی

رویکرد رایج، انتخاب دو فرکانس کلیدی از سیستم (مانند ωᵢ و ωⱼ) و نسبت‌های میرایی مطلوب برای آن‌ها (ξᵢ و ξⱼ) است. با جایگذاری این مقادیر در فرمول نسبت میرایی، یک دستگاه دو معادله و دو مجهول برای یافتن α و β تشکیل می‌شود که به راحتی قابل حل است.

استراتژی‌های انتخاب فرکانس

انتخاب دو فرکانس مناسب، کلید موفقیت این روش است. در ادامه استراتژی‌های رایج به همراه کاربردها و ملاحظات آن‌ها تشریح می‌شود:

۱. استفاده از فرکانس‌های طبیعی دو مود اول سازه

این رویکرد یک نقطه شروع متداول برای سیستم‌هایی است که پاسخ دینامیکی آن‌ها تحت سلطه مودهای ارتعاشی پایه (مود اول و دوم) قرار دارد.

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

۲. استفاده از فرکانس مود اول و فرکانس غالب بارگذاری

این استراتژی رویکرد صحیح برای تحلیل‌هایی است که شامل یک بار هارمونیک مشخص، مانند بارهای ناشی از ماشین‌آلات دوار، می‌شوند.

• نقاط قوت: مدل میرایی را مستقیماً به مسئله فیزیکی تحت بررسی مرتبط می‌کند و دقت را در فرکانس تحریک خارجی به حداکثر می‌رساند.
• ملاحظات: برای بارهای با طیف فرکانسی گسترده مانند زلزله یا ارتعاشات تصادفی، این روش به تنهایی کافی نیست.

۳. استفاده از فرکانس مود اول و آخرین مود مهم

این استراتژی، دقیق‌ترین و قوی‌ترین روش مهندسی، به ویژه برای تحلیل‌های دینامیکی پیچیده مانند تحلیل‌های لرزه‌ای یا ارتعاشات تصادفی است. در این روش، «آخرین مود مهم» به عنوان مودی تعریف می‌شود که برای رسیدن به مشارکت جرمی تجمعی (Mass Participation) حدود ۹۵٪ لازم است.

• نقاط قوت: این روش تضمین می‌کند که مدل میرایی در سراسر بازه فرکانسی که به طور مادی در پاسخ دینامیکی سیستم نقش دارد، رفتار واقع‌بینانه‌ای از خود نشان می‌دهد. این رویکرد از میرایی بیش از حد یا کمتر از حد در مودهای مهم جلوگیری کرده و دقت نتایج را به شدت افزایش می‌دهد.
• ملاحظات: نیازمند انجام یک تحلیل مودال اولیه برای شناسایی مودهای با مشارکت جرمی بالا است که اندکی زمان‌برتر است، اما دقت بالاتری را به ارمغان می‌آورد.

پیاده‌سازی در Ansys

در نرم‌افزار Ansys، پس از محاسبه ضرایب α و β، این مقادیر از طریق دستورات مشخصی به مدل اعمال می‌شوند. رایج‌ترین دستورات عبارتند از:

• ALPHAD: برای تعریف ضریب میرایی وابسته به جرم (α).
• BETAD: برای تعریف ضریب میرایی وابسته به سختی (β).

این گام عملی، مفاهیم تئوریک را به یک ابزار قابل استفاده در شبیه‌سازی‌های مهندسی پیشرفته تبدیل می‌کند.

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

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

• میرایی رایلی یک مدل محاسباتی کارآمد و پرکاربرد ([C] = α[M] + β[K]) است که پیاده‌سازی آن در نرم‌افزارهای المان محدود به سادگی امکان‌پذیر است.
• رفتار این مدل وابسته به فرکانس است؛ در فرکانس‌های پایین میرایی وابسته به جرم و در فرکانس‌های بالا میرایی وابسته به سختی غالب است که این امر نیازمند توجه ویژه در تفسیر نتایج است.
• کاربرد صحیح آن به انتخاب هوشمندانه ضرایب α و β بر اساس ویژگی‌های دینامیکی سیستم و نوع بارگذاری بستگی دارد.

کاربرد گسترده «میرایی رایلی در انسیس» نشان‌دهنده اهمیت این مفهوم در شبیه‌سازی‌های حرفه‌ای است که دقت و واقع‌گرایی تحلیل‌ها را به شکل چشمگیری افزایش می‌دهد.

به عنوان توصیه نهایی و کلیدی، انتخاب متفکرانه ضرایب α و β بر اساس مودهای کلیدی مشارکت‌کننده در پاسخ سیستم (مانند رویکرد مبتنی بر مشارکت جرمی ۹۵٪) صرفاً «واقع‌بینانه‌تر» نیست، بلکه برای دستیابی به یک شبیه‌سازی دینامیکی دقیق و قابل اعتماد ضروری است. اگرچه فرض میرایی ثابت برای تمام مودها یک ساده‌سازی رایج است، اما یک رویکرد مهندسی صحیح، انتخاب ضرایبی است که رفتار میرایی را در مهم‌ترین فرکانس‌های مؤثر بر پاسخ دینامیکی سیستم به درستی مدل کند.