ارتعاشات اجباری نامیرا با تحریک نیروی هارمونیک – از صفر تا صد

۲۹۴۱ بازدید
آخرین به‌روزرسانی: ۳۱ اردیبهشت ۱۴۰۲
زمان مطالعه: ۸ دقیقه
دانلود PDF مقاله
ارتعاشات اجباری نامیرا با تحریک نیروی هارمونیک – از صفر تا صدارتعاشات اجباری نامیرا با تحریک نیروی هارمونیک – از صفر تا صد

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

997696

معادله حرکت در ارتعاشات اجباری نامیرا

نیروی F(t)\large F(t) به سیستم شکل زیر وارد می‌شود که متشکل از جرم، فنر و میراگر است.

با کمک قانون دوم نیوتن، معادله حرکت را برای این سیستم می‌نویسیم.

سیستم جرم و فنر

mx¨+cx˙+kx=F(t)\large m\ddot{x} \:+\: c\dot{x} \:+\: kx \:=\: F(t)

(رابطه ۱)

از آنجایی که معادله بالا غیرهمگن است، پاسخ عمومی آن (x(t)\large x(t)) به صورت مجموع پاسخ همگن (xh(t)\large x_h(t)) و پاسخ خصوصی (xp(t)\large x_p(t)) به دست می‌آید. پاسخ همگن، همان پاسخ معادله زیر است که ارتعاشات آزاد را نشان می‌دهد.

mx¨+cx˙+kx=0\large m\ddot{x} \:+\: c\dot{x} \:+\: kx \:=\:0

همان‌طور که می‌دانید، ارتعاشات آزاد تحت هرگونه وضعیت میرایی (زیر میرایی، میرایی بحرانی و فوق میرایی) و با شرایط اولیه مختلف، پس از سپری شدن مدت زمان مشخصی، از حرکت باز می‌ایستد. بنابراین رابطه ۱، ناگهان به پاسخ خصوصی xp(t)\large x_p(t) تبدیل می‌شود که نشان دهنده ارتعاشات حالت پایدار است. حرکت حالت پایدار تا مادامی ادامه دارد که تابع نیرو، وجود داشته باشد. تغییر پاسخ‌های همگن، خصوصی و عمومی با زمان برای یک سیستم فرضی در شکل زیر رسم شده است. همان‌طور که مشاهده می‌کنید، پس از سپری شدن مدتی (تا رسیدن نمودار به زمان τ\large \tauxh(t)\large x_h(t) از بین رفته و x(t)\large x(t) با xp(t)\large x_p(t) برابر می‌شود.

قسمتی از حرکت که به دلیل میرایی از بین می‌رود (بخش ارتعاشات آزاد)، قسمت گذرا نامیده می‌شود. سرعت از بین رفتن حرکت گذرا، به مقدار پارامترهای c\large c، k\large k و m\large m وابسته است. در ادامه این مقاله، اثر ارتعاشات اجباری نامیرا در سیستم‌هایی که با نیروی هارمونیک تحریک می‌شوند، دنبال خواهد شد.

سیستم زیر میرایی

پاسخ سیستم در ارتعاشات اجباری نامیرا با تحریک نیروی هارمونیک

اگر نیروی F(t)=F0cosωt\large F(t) =F_0 \cos \omega t به جرم m\large m در یک سیستم نامیرا وارد شود، معادله حرکت جدید به صورت زیر است.

mx¨+kx=F0cosωt\large m\ddot{x} \:+\: kx \:=\: F_0 \cos \omega t

(رابطه ۲)

پاسخ همگن این معادله را می‌نویسیم.

xh(t)=C1cosωnt+C2sinωnt\large x_h(t) \:=\: C_1\cos \omega_n t \:+\: C_2 \sin \omega_n t

در عبارت بالا، ωn=k/m\large \omega_n =\sqrt {k/m} فرکانس طبیعی سیستم است. از آنجایی که نیروی تحریک (F(t)\large F(t)) هارمونیک است، پاسخ عمومی xp(t)\large x_p(t) نیز هارمونیک بوده و فرکانس آن نیز برابر ω\large \omega است. پاسخ خصوصی به شکل زیر است.

xp(t)=Xcosωt\large x_p(t) \:=\: X \:\cos \omega t

(رابطه ۳)

در عبارت بالا، X\large X ثابت است و بیشترین دامنه xp(t)\large x_p(t) را نشان می‌دهد. با ادغام رابطه‌های ۲ و ۳ و مرتب کردن آن برحسب X\large X به رابطه زیر می‌رسیم.

X=F0kmω2=δst1(ωωn)2\large X\:=\: \frac {F_0} {k \:-\: m\omega^2} \:=\: \frac {\delta _{st}} {1 \:-\: (\frac {\omega} {\omega_n})^2}

(رابطه ۴)

عبارت δst=F0/k\large \delta _{st}= F_0/k جابجایی جرم m\large m، تحت نیروی F0\large F_0 است و از آنجایی که F0\large F_0 یک نیروی ثابت (استاتیک) به حساب می‌آید، جابجایی استاتیکی نیز خوانده می‌شود. از این رو، پاسخ کلی رابطه ۲ به صورت زیر خواهد بود.

x(t)=C1cosωnt+C2sinωnt+F0kmω2cosωt\large x(t) \:=\: C_1 \cos \omega_n t \:+\: C_2 \sin \omega_n t \:+\: \frac {F_0} {k\:-\: m\omega^2} \cos \omega t

برای یافتن ضریب‌های ثابت، شرایط اولیه را به صورت x(t=0)=x0\large x(t=0) =x_0 و x˙(t=0)=x˙0\large \dot{x}(t=0) =\dot{x}_0 تعریف می‌کنیم.

C1=x0F0kmω2 C2=x˙0ωn\large C_1\:=\: x_0\:-\: \frac {F_0} {k\:-\: m\omega^2} \\~\\ \large C_2\:=\: \frac {\dot {x}_0} {\omega_n}

حالا می‌توانیم پاسخ کلی را بازنویسی کنیم.

x(t)=(x0F0kmω2)cosωnt+(x˙0ωn)sinωnt +(F0kmω2)cosωt\large x(t) \:=\: (x_0 \:-\: \frac {F_0} {k \:-\: m\omega^2}) \cos \omega_n t \:+\: (\frac {\dot {x} _0} {\omega_n}) \sin \omega_n t \\~\\ \large + (\frac {F_0} {k \:-\: m\omega^2}) \cos \omega t

(رابطه ۵)

همچنین، ماکزیمم مقدار X\large X در رابطه ۴ نیز برابر با عبارت زیر است.

Xδst=11(ωωn)2\large \frac {X} {\delta _{st}} \:=\: \frac {1} {1\:-\: (\frac {\omega} {\omega_n})^2}

(رابطه 6)

مقدار X/δst\large X/ \delta_ {st} نشان‌دهنده نسبت اندازه حرکت دینامیکی به اندازه حرکت استاتیکی است و ضریب بزرگ‌نمایی (Magnification Factor)‌ یا نسبت دامنه (Amplitude Ratio) نام دارد. تغییرات نسبت دامنه برحسب تغییرات نسبت فرکانس (r=ω/ωn\large r= \omega/ \omega_n) در شکل زیر قابل مشاهده است. با توجه به این شکل، پاسخ سیستم می‌تواند سه حالت مختلف داشته باشد.

ضریب بزرگنمایی

حالت اول زمانی رخ می‌دهد که 0<ω/ωn<1\large 0<\omega/ \omega_n <1 برقرار باشد. در این حالت، مخرج کسر در رابطه 6، مثبت بوده و پاسخ آن مانند رابطه ۳ و بدون تغییر خواهد بود. شکل زیر را در نظر بگیرید. پاسخ هارمونیک سیستم xp(t)\large x_p(t) با نیروی خارجی هم‌فاز است.

پاسخ هارمونیک

در حالت دوم، رابطه ω/ωn>1\large \omega/ \omega_n >1 برقرار است. در نتیجه، مخرج کسر در رابطه 6 منفی می‌شود و پاسخ حالت ماندگار به صورت زیر خواهد بود.

xp(t)=Xcosωt\large x_p(t) \:=\: -\:X \cos \omega t

در رابطه بالا، دامنه حرکت X\large X طوری تعریف می‌شود که همواره مثبت باشد.

X=δst(ωωn)21\large X\:=\: \frac {\delta_ {st}} {(\frac {\omega} {\omega_n})^2 \:-\:1}

(رابطه 7)

تغییرات F(t)\large F(t) و xp(t)\large x_p(t) را نسبت به زمان، در شکل زیر ملاحظه می‌کنید. از آنجایی که F(t)\large F(t) و xp(t)\large x_p(t) مقداری مثبت دارند، پاسخ آنها 180\large 180^ \circ از نیروی خارجی عقب‌تر است. علاوه بر آن، اگر ω/ωn\large \omega/ \omega_n \rightarrow \infty، آنگاه X0\large X \rightarrow 0 برقرار خواهد بود. بنابراین، پاسخ سیستم به یک نیروی هارمونیک با فرکانس زیاد، نزدیک به صفر است.

نیروی هارمونیک

اما در حالت سوم که به صورت ω/ωn=1\large \omega/ \omega_ n=1 تعریف می‌شود، اندازه X\large X در رابطه‌های 6 و 7، بی‌نهایت خواهد شد. این حالت که در آن، فرکانس نیرو ω\large \omega با فرکانس طبیعی سیستم ωn\large \omega_n برابر است، رزونانس یا تشدید نام دارد. برای یافتن پاسخ سیستم در این حالت، رابطه 5 را بازنویسی می‌کنیم.

x(t)=x0cosωnt+x˙0ωnsinωnt+δst[cosωtcosωnt1(ωωn)2]\large x(t) \:=\: x_0 \cos \omega _nt \:+\: \frac {\dot{x} _0} {\omega _n} \sin \omega _nt \:+\: \delta _{st} [\frac {\cos \omega t\:-\: \cos \omega _nt} {1\:-\: (\frac {\omega} {\omega _n})^2}]

از آنجایی که عبارت سمت راست در رابطه بالا، در حالت ω=ωn\large \omega =\omega _n، تعریف نشده خواهد بود، برای محاسبه مقدار حدی آن، از قاعده هوپیتال استفاده می‌کنیم.

limωωn[cosωtcosωnt1(ωωn)2]=limωωn[ddω(cosωtcosωnt)ddω(1(ωωn)2)] =limωωn[tsinωt2ωωn2]=ωnt2sinωnt\large \lim_{\omega \rightarrow \omega_n} [\frac {\cos \omega t\:-\: \cos \omega _nt} {1\:-\: (\frac {\omega} {\omega _n})^2}] \:=\: \lim_{\omega \rightarrow \omega_n} [\frac {\frac {\text {d}}{\text {d} \omega} (\cos \omega t\:-\: \cos \omega _nt)} {\frac {\text {d}}{\text {d} \omega}(1\:-\: (\frac {\omega} {\omega _n})^2)}] \\~\\ \large =\: \lim_{\omega \rightarrow \omega_n} [\frac {t\sin \omega t} {2\frac {\omega} {\omega_n^2}}] \:=\: \frac {\omega _nt} {2} \sin \omega _nt

حالا می‌توانیم پاسخ سیستم را به صورت کامل بنویسیم.

x(t)=x0cosωnt+x˙0ωnsinωnt+δstωnt2sinωnt\large x(t) \:=\: x_0 \cos \omega _nt \:+\: \frac {\dot{x} _0} {\omega _n} \sin \omega _nt \:+\: \frac {\delta _{st} \omega _nt} {2} \sin \omega _nt

از این رابطه متوجه می‌شویم که در رزونانس، مقدار x(t)\large x(t) تا بی‌نهایت افزایش می‌یابد. آخرین عبارت در رابطه بالا را می‌توانید در شکل پایین مشاهده کنید. دامنه رزونانس به صورتی خطی با زمان در حال زیاد شدن است.

پاسخ خصوصی معادله

پاسخ نهایی ارتعاشات اجباری نامیرا با تحریک نیروی هارمونیک

پاسخ نهایی سیستم مربوط به رابطه ۵ را می‌توانیم به صورت زیر بنویسیم.

x(t)=Acos(ωntϕ)+δst1(ωωn)2cosωt         ωωn<1\large x(t) \:=\: A \cos (\omega _nt \:-\: \phi) \:+\: \frac {\delta _{st}} {1\:-\: (\frac {\omega} {\omega_n})^2} \cos \omega _t ~~~ ~~~ ~~~ \frac {\omega} {\omega_n} <1

(رابطه ۸)

x(t)=Acos(ωntϕ)+δst1+(ωωn)2cosωt         ωωn>1\large x(t) \:=\: A \cos (\omega _nt \:-\: \phi) \:+\: \frac {\delta _{st}} {-\:1\:+\: (\frac {\omega} {\omega_n})^2} \cos \omega _t ~~~ ~~~ ~~~ \frac {\omega} {\omega_n} >1

(رابطه ۹)

بنابراین حرکت کلی به صورت مجموع دو منحنی کسینوسی با فرکانس‌های مختلف نشان داده می‌شود. در رابطه‌های ۸ و ۹، فرکانس نیرو (ω\large \omega) به ترتیب از فرکانس طبیعی کوچکتر و بزرگتر است. پاسخ کلی سیستم در شکل زیر نشان داده شده است.

پاسخ کلی معادله

پدیده ضربان در ارتعاشات اجباری نامیرا

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

اگر شرایط اولیه را به صورت x0=x˙0=0\large x_0\:=\: \dot {x}_0 \:=\:0 فرض کنیم، رابطه ۵ به شکل زیر نوشته می‌شود.

x(t)=(F0/m)ωn2ω2(cosωtcosωnt) =(F0/m)ωn2ω2[2sinωn+ω2t.sinωnω2t]\large x(t)\:=\: \frac {(F_0/m)} {\omega_n ^2\:-\: \omega ^2} (\cos \omega t\:-\: \cos \omega _nt) \\~\\ \large =\: \frac {(F_0/m)} {\omega_n ^2\:-\: \omega ^2} [2\sin \frac {\omega_n \:+\: \omega} {2} t.\:\sin \frac {\omega_n \:-\: \omega} {2} t]

حالا فرض کنیم فرکانس نیرو اندکی از فرکانس طبیعی کمتر باشد.

ωnω=2ϵ\large \omega_n \:-\:\omega \:=\: 2\epsilon

مقدار مثبت بسیار کوچک را با ϵ\large \epsilon نشان داده‌ایم.

ω+ωn=ω+ω+ϵ2ω\large \omega \:+\: \omega_n \:=\: \omega \:+\: \omega \:+\: \epsilon \:\cong\: 2\omega

اکنون اگر طرفین دو رابطه اخیر را در هم ضرب کنیم، رابطه زیر به دست می‌آید.

ωn2ω2=4ϵω\large \Rightarrow \omega_n ^2\:-\: \omega^2 \:=\: 4\epsilon \omega

به این ترتیب می‌توانیم مجدداً پاسخ x(t)\large x(t) را بنویسیم.

x(t)=(F0/m2ϵωsinϵt)sinωt\large x(t)\:=\: (\frac {F_0/m} {2\epsilon \omega} \sin \epsilon t) \:\sin \omega t

از آنجایی که ϵ\large \epsilon مقداری بسیار کوچک است، تابع sinϵt\large \sin \epsilon t به کُندی تغییر کرده و دوره تناوبی بزرگ و برابر با 2π/ϵ\large 2\pi/ \epsilon خواهد داشت. به همین دلیل، رابطه بالا ارتعاشات اجباری نامیرا با دوره تناوب 2π/ϵ\large 2\pi/ \epsilon و دامنه‌ای متغیر را نشان می‌دهد. از طرفی دیگر، همان‌طور که در شکل زیر نیز مشهود است، در مدت زمانی که منحنی sinϵt\large \sin \epsilon t فقط یک سیکل را کامل می‌کند، منحنی sinωt\large \sin \omega t چند سیکل نوسان می‌کند. مدت زمان بین هر دو دامنه صفر یا هر دو دامنه ماکزیمم، دوره ضربان نام دارد و با نماد τb\large \tau_b تعریف می‌شود.

پدیده ضربان در ارتعاشات

τb=2π2ϵ=2πωnω\large \tau_b \:=\: \frac {2\pi} {2\epsilon} \:=\: \frac {2\pi} {\omega _n\:-\: \omega}

از این رو، فرکانس ضربان نیز برابر با ωb=2ϵ=ωnω\large \omega _b=2 \epsilon= \omega _n- \omega است.

مثال: صفحه نگهدارنده یک پمپ

سؤال: یک پمپ رفت و برگشتی به وزن 150lb\large 150\: lb در مرکز صفحه فولادی قرار گرفته است. طول و عرض این صفحه به ترتیب 100\large 100 و 20\large 20 اینچ بوده و ضخامت آن نیز 0.5\large 0.5 اینچ است. در حین روشن بودن پمپ، این صفحه در معرض یک نیروی هارمونیک قرار می‌گیرد که اندازه این نیرو برابر با F(t)=50cos62.832t\large F(t)\:=\: 50\cos 62.832 t و با واحد lb\large lb‌ است. دامنه ارتعاشات اجباری نامیرا در صفحه را بیابید.

مثال ارتعاشات اجباری

پاسخ: این صفحه را می‌توانیم مانند تیر دو سر گیردار و با مدول یانگ E=30×106psi\large E\:=\: 30 \times 10^6 \:psi مدل‌سازی کنیم. ابتدا ممان اینرسی سطح را به دست می‌آوریم.

I=112(20)(0.5)3=0.2083in4\large I\:=\: \frac {1} {12} (20) (0.5)^3 \:=\: 0.2083\: in^4

برای به دست آوردن سفتی خمشی (Bending Stiffness) از رابطه زیر استفاده می‌کنیم.

k=192EIl3=192(30×106)(0.2083)(100)3=1200.0lb/in\large k\:=\:\frac {192EI} {l^3} \:=\: \frac {192 (30 \times 10^6) (0.2083)} {(100)^3} \:=\: 1200.0 \:lb/in

دامنه پاسخ هارمونیک با کمک رابطه ۴ محاسبه می‌شود. با جایگذاری مقادیر مختلف، X\large X به دست می‌آید.

F0=50lb         m=150/386.4lbsec2/in k=1200lb/in         ω=62.832rad/s X=F0kmω2=501200(150/386.4)(62.832)2    X=0.1504in\large F_0\:=\: 50\: lb ~~~ ~~~ ~~~ m\:=\: 150/ 386.4\: lb-sec^2 /in \\~\\ \large k\:=\: 1200 lb/in ~~~ ~~~ ~~~ \omega\:=\: 62.832\: rad/s \\~\\ \large X\:=\: \frac {F_0} {k\:-\: m\omega^2} \:=\: \frac {50} {1200\:-\: (150/386.4) (62.832)^2} \\~\\ \large \Rightarrow ~~~ X\:=\:-\: 0.1504 \:in

علامت منفی نشان می‌دهد که پاسخ x(t)\large x(t) مربوط به صفحه فولادی، نسبت به تحریک F(t)\large F(t) اختلاف فاز دارد.

مثال: تعیین جرم با استفاده از پاسخ هارمونیک در ارتعاشات اجباری نامیرا

سؤال: به یک سیستم جرم و فنر، نیروی هارمونیک به بزرگی 30N\large 30 \:N و فرکانس 20Hz\large 20 \:Hz وارد می‌شود. ثابت فنر برابر 5,000N/m\large 5,000 \:N/m است و جرم با دامنه‌ 0.2m\large 0.2 \:m نوسان می‌کند. با فرض اینکه ارتعاشات اجباری نامیرا از حالت سکون (x0=x˙0=0\large x_0 =\dot {x}_ 0=0)‌ شروع شده باشد، جرم سیستم را تعیین کنید.

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

x(t)=F0kmω2(cosωtcosωnt)\large x(t) \:=\: \frac {F_0} {k\:-\: m\omega^2} (\cos \omega t\:-\: \cos \omega _nt)

با کمک روابط مثلثاتی و بازنویسی رابطه بالا، به شکل جدیدی برای x(t)\large x(t) می‌رسیم.

x(t)=2F0kmω2sinωn+ω2tsinωnω2t\large x(t) \:=\: \frac {2F_0} {k\:-\: m\omega^2} \sin \frac {\omega_n \:+\: \omega} {2}t \sin \frac {\omega_n \:-\: \omega} {2}t

از طرفی، دامنه این ارتعاشات اجباری نامیرا معلوم و برابر با 0.2m\large 0.2 \:m است.

2F0kmω2=0.2\large \frac {2F_0} {k\:-\: m\omega^2} \:=\: 0.2

اکنون می‌توانیم جرم سیستم را با داشتن سایر متغیرها، محاسبه کنیم.

F0=30N         ω=20Hz=125.665rad/s k=5,000N/m    2(30)5000m(125.664)2=0.2    m=0.2976kg\large F_0 \:=\: 30\:N ~~~ ~~~ ~~~ \omega \:=\: 20 \:Hz \:=\: 125.665\: rad/s \\~\\ \large k\:=\: 5,000\: N/m \\~\\ \large \Rightarrow ~~~ \frac {2(30)} {5000 \:-\: m(125.664)^2} \:=\: 0.2 \\~\\ \large \Rightarrow ~~~ m\:=\: 0.2976 \:kg

اگر به مباحث مرتبط در زمینه مکانیک و ارتعاشات علاقه‌مند هستید، آموزش‌های زیر به شما پیشنهاد می‌شوند:

^^

بر اساس رای ۱۰ نفر
آیا این مطلب برای شما مفید بود؟
اگر بازخوردی درباره این مطلب دارید یا پرسشی دارید که بدون پاسخ مانده است، آن را از طریق بخش نظرات مطرح کنید.
منابع:
Mechanical Vibrations of Rao
دانلود PDF مقاله
۱ دیدگاه برای «ارتعاشات اجباری نامیرا با تحریک نیروی هارمونیک – از صفر تا صد»

سلام این پدیده در دینامیک خاک هم مطرح است اما کتاب مکانیک داس جواب پدیده ضربان را غ ق ق دانسته.چرا؟

نظر شما چیست؟

نشانی ایمیل شما منتشر نخواهد شد. بخش‌های موردنیاز علامت‌گذاری شده‌اند *