این پایاننامه به توسعه یک سامانه محاسباتی میپردازد که قادر است رفتار رشتههای بسیار نازک و منعطف را درون یک سیال چسبناک (مثل آب یا هوا) شبیهسازی کند. اهمیت این پژوهش در دو حوزهٔ کلیدی نهفته است: نخست، درک پدیدههای طبیعی مانند حرکت باکتریها یا رفتار الیاف در بدن موجودات زنده؛ و دوم، کاربرد مهندسی آن در طراحی روشی نوین برای مقابله با پهپادها. این روش شامل پرتاب شبکهای از رشتهها به سمت ملخهای پهپاد است تا با گرهخوردن و گیر کردن در آنها، باعث کاهش نیروی رانش و در نهایت سقوط کنترلشدهٔ پهپاد شود.
برای مدلسازی خود رشتهها، از یک نظریهٔ پیشرفتهٔ تیرک (Kirchhoff-Love) استفاده شده است که میتواند خمشها، پیچشها و کشیدگیهای بزرگ را با دقت بسیار بالا محاسبه کند. این مدلسازی برای رشتههایی که طول آنها هزاران برابر قطرشان است، بهینه شده و امکان شبیهسازی تغییر شکلهای پیچیده مانند گرهخوردن و تابخوردن را فراهم میسازد. بخش مهم دیگر این پژوهش، مدلسازی تماس و اصطکاک میان خود رشتهها و همچنین میان رشتهها و ملخ پهپاد است که با استفاده از الگوریتمهای ریاضی قدرتمندی حل شده است.
از سوی دیگر، برای شبیهسازی محیط سیال (هوا یا مایع اطراف رشتهها) از روش «انتگرال مرزی» بهره گرفته شده است. این روش به جای محاسبهٔ کل میدان جریان، تنها روی سطح رشتهها تمرکز میکند و به همین دلیل بسیار سریع و کارآمد است. این روش که مبتنی بر معادلات «استوکس» برای جریانهای آهسته است، به خوبی تعامل بین حرکت رشتهها و نیروهای وارد از سوی سیال را محاسبه میکند. اتصال این دو بخش (مدل جامد رشته و مدل سیال) از طریق اعمال شرط عدملغزش روی سطح رشتهها انجام میشود.
یکی از چالشهای اساسی در این پژوهش، تشخیص دقیق نقطهٔ تماس بین رشتههایی است که به شدت خمیده شدهاند. روشهای سادهٔ قبلی در چنین شرایطی شکست میخوردند، بنابراین در این کار روش جدیدی ابداع شده که با استفاده از تبدیلهای ریاضی و بهینهسازی، نقاط تماس را با دقتی بالا و حتی در بدترین حالتهای خمیدگی پیدا میکند. همچنین مدل اصطکاک کولمب برای شبیهسازی لغزش و چسبندگی رشتهها روی یکدیگر پیادهسازی شده است که برای ایجاد گرههای پایدار ضروری است.
برای کاربرد اصلی این تحقیق (مقابله با پهپاد)، یک مدل کامل برای ملخ پهپاد و موتور الکتریکی آن ساخته شده است. این مدل نشان میدهد که چگونه برخورد رشتهها با ملخ، گشتاور آن را تغییر داده و سرعت چرخش را کاهش میدهد. سپس با استفاده از یک مدل آیرودینامیکی ساده، میزان افت نیروی رانش محاسبه میشود. بدین ترتیب، میتوان اثر تعداد رشتهها، جرم آنها و سرعت برخورد را بر عملکرد پهپاد کمی بررسی کرد.
نتایج شبیهسازیها در این پایاننامه نشان میدهد که با افزایش تعداد رشتهها و به ویژه افزایش چگالی جرمی آنها، افت نیروی رانش به طور قابل توجهی بیشتر میشود. به عنوان مثال، دو برابر کردن چگالی رشتهها میتواند کاهش رانش را تا بیش از ۲۸٪ افزایش دهد. همچنین وجود اصطکاک میان رشتهها نقش حیاتی در پایدار ماندن گرهها دارد و در غیاب آن، رشتهها به سرعت از هم باز میشوند و اثربخشی خود را از دست میدهند.
این چارچوب محاسباتی نشان میدهد که شبیهسازی پدیدههای پیچیدهٔ فیزیکی با استفاده از ترکیب روشهای ریاضی پیشرفته و قدرت پردازش موازی (با استفاده از پردازندههای گرافیکی) امکانپذیر است. یافتههای این پژوهش نه تنها به درک بهتر رفتار الیاف در طبیعت کمک میکند، بلکه مسیری عملی برای طراحی سامانههای دفاعی غیرتهاجمی و مؤثر در برابر پهپادها ارائه میدهد. این روش میتواند جایگزین مناسبی برای روشهای خطرناک و پرهزینهٔ فعلی مانند شلیک موشک یا امواج مایکروویو باشد.
چارچوب توسعهیافته در این پایاننامه از سه بخش اصلی تشکیل شده است: مدل مکانیک جامد برای رشتهها، مدل تماس و اصطکاک برای برهمکنشهای نزدیک، و مدل دینامیک سیالات برای محیط اطراف. این سه بخش به گونهای کوپل شدهاند که تغییر شکل رشتهها، نیروهای وارد از سیال، و نیروهای تماسی به صورت همزمان و در هر لحظه محاسبه شوند. برای حل معادلات حاصل، از روشهای عددی پیشرفته مانند روش نیوتن-رافسون و الگوریتمهای بازگشت به عقب (return mapping) برای مدلسازی اصطکاک استفاده شده است.
یکی از نوآوریهای مهم این پژوهش، ارائه روشی برای تشخیص تماس بین رشتههایی است که با توابع درجه سوم (هرمیتی) نمایش داده شدهاند. این توابع امکان نمایش خمیدگیهای پیوسته و نرم را فراهم میکنند، اما تشخیص نقطه تماس میان آنها بسیار دشوار است. در این پایاننامه، با استفاده از تبدیل زاویه فاز و فرمولبندی مبتنی بر مجذور فاصله، یک روش نیوتن محلی طراحی شده که همگرایی مطمئن و سریعی حتی در حالتهای خمیدگی شدید دارد. این روش باعث شده تا شبیهسازی تماسهای پیچیده مانند گرههای چندلایه به صورت پایدار انجام شود.
در بخش سیالات، روش انتگرال مرزی با استفاده از توزیع استوکسلت در طول مرکز رشتهها، میدان سرعت را محاسبه میکند. برای جلوگیری از تکینگی ریاضی در نقاط نزدیک به هم، از یک مدل تماس برای جلوگیری از برخورد رشتهها استفاده شده است. همچنین برای افزایش سرعت محاسبات، ماتریس هسته (که بسیار بزرگ و متراکم است) با استفاده از پردازنده گرافیکی و کتابخانه cuSOLVER حل میشود. این امر باعث شده تا شبیهسازیهایی با بیش از ۳۶ رشته و هر کدام با ۱۱ المان، در زمانی معقول انجام شود.
مدل ملخ و موتور الکتریکی نیز به صورت کوپل شده با بقیه بخشها طراحی شده است. ملخ به عنوان دو تیرک عمود بر هم با شرایط مرزی دیریکله (حرکت اجباری) مدلسازی شده و موتور الکتریکی با استفاده از مدار معادل DC براشدار شبیهسازی میشود. گشتاور موتور، تابعی از سرعت دورانی و جریان الکتریکی است و گشتاورهای ناشی از برخورد رشتهها به عنوان نیروهای خارجی به معادلات حرکت ملخ اعمال میشوند. این کوپلینگ دوطرفه باعث میشود که هر تماس، سرعت ملخ را کاهش دهد و در نتیجه نیروی رانش تولیدی کاهش یابد.
برای محاسبه نیروی رانش، از تئوری مومنتوم ساده شده استفاده شده است که در آن ضریب رانش و ضریب توان به عنوان توابعی از نسبت پیشروی (سرعت جریان آزاد بر سرعت نوک ملخ) تعریف میشوند. این مدل علیرغم سادگی، دقت مناسبی برای تخمین افت رانش در شرایط مختلف دارد. در شبیهسازیهای انجام شده، رشتهها با چگالیهای مختلف (۱۳۰۰، ۲۵۰۰ و ۵۰۰۰ کیلوگرم بر مترمکعب) و تعداد متفاوت (۱۲ و ۳۶ رشته) در مسیر ملخ رها شدهاند و نتایج نشان میدهد که چگالی جرمی تأثیر بسیار بیشتری نسبت به تعداد رشتهها بر افت رانش دارد.
یکی از سناریوهای جالب شبیهسازی شده، رهاسازی یک رشتهٔ پیچیدهشده به شکل گره است. در این آزمایش، رشته ابتدا به صورت استاتیکی به شکل گره درمیآید و سپس به صورت دینامیکی رها میشود. نتایج نشان میدهد که در غیاب اصطکاک، گره به سرعت باز میشود، اما در حضور اصطکاک، گره پایدار میماند و میتواند اثر مخربی بر ملخ داشته باشد. این موضوع اهمیت اصطکاک را در کاربرد مقابله با پهپاد به خوبی نشان میدهد، زیرا رشتهها باید پس از برخورد با ملخ، گرههای پایدار تشکیل دهند تا بتوانند مانع چرخش شوند.
در مجموع، این پایاننامه نشان میدهد که با ترکیب روشهای دقیق مکانیک جامد، سیالات و تماس، میتوان پدیدههای پیچیدهای مانند درهمتنیدگی رشتهها و اثر آن بر سیستمهای مکانیکی را شبیهسازی کرد. این چارچوب میتواند به عنوان ابزاری برای بهینهسازی طراحی سامانههای دفاعی مبتنی بر رشته عمل کند و همچنین در سایر حوزهها مانند مهندسی نساجی، زیستشناسی سلولی و علم مواد کاربرد داشته باشد. پیشنهادات برای کارهای آینده شامل افزودن مدلهای جریان آشفته، بهبود الگوریتم تماس برای حالتهای موازی، و انجام مطالعات بهینهسازی جامع برای یافتن بهترین خواص رشتهها است.
روشهای عددی استفاده شده در این تحقیق شامل روش اجزای محدود برای تیرکها، روش انتگرال مرزی برای سیالات، و روش پنالتی برای تماس است. هر کدام از این روشها به دلیل کارایی بالا و دقت مناسب در مسائل رشتههای باریک انتخاب شدهاند. همچنین برای حل معادلات غیرخطی حاصل از تماس و تغییر شکل بزرگ، از روش نیوتن-رافسون با گامهای زمانی تطبیقی استفاده شده است. تمام کدهای نوشته شده با استفاده از کتابخانههای متنباز PETSc و CUDA پیادهسازی شدهاند تا از قدرت پردازش موازی برای کاهش زمان محاسبات استفاده شود.
نتایج کمی این تحقیق نشان میدهد که با ۳۶ رشته با چگالی ۱۳۰۰ کیلوگرم بر مترمکعب، حدود ۱۵٪ کاهش رانش حاصل میشود، در حالی که همین تعداد رشته با چگالی ۵۰۰۰ کیلوگرم بر مترمکعب، کاهش رانش را به ۵۵٪ میرساند. این افزایش شدید نشان میدهد که جرم رشتهها پارامتر کلیدی در اثربخشی این روش است. همچنین مشاهده شد که رشتههای با چگالی بیشتر، به دلیل نیروی گرانشی بالاتر، سریعتر به سمت ملخ جذب میشوند و تماسهای بیشتری ایجاد میکنند. از سوی دیگر، رشتههای بسیار سنگین ممکن است زودتر از فضای برخورد خارج شوند، بنابراین یک محدودهٔ بهینه برای جرم رشتهها وجود دارد که باید در طراحی لحاظ شود.
در نهایت، این تحقیق گامی مهم در جهت شبیهسازی واقعگرایانهٔ سیستمهای فیبری در تعامل با سیالات و اجسام متحرک است. این چارچوب میتواند در آینده برای مطالعهٔ مسائل دیگری مانند رفتار الیاف در فرآیندهای ریسندگی، حرکت دستههای باکتری در محیطهای چسبناک، و حتی طراحی تورهای صید ماهی مورد استفاده قرار گیرد. گسترش این روش به حالتهای سهبعدی کامل و شامل کردن اثرات حرارتی و الکتریکی، میتواند افقهای جدیدی را در علم مواد و مهندسی زیستی بگشاید.
| سرفصل | شماره صفحه |
|---|---|
| صفحه عنوان | ۱ |
| چکیده | ۳ |
| سپاسگزاری | ۵ |
| فهرست مطالب | ۷ |
| فهرست شکلها | ۹ |
| فهرست جدولها | ۱۳ |
| فصل ۱: مقدمه | ۱۵ |
| ۱-۱: پیشینه و انگیزه | ۱۵ |
| ۱-۲: مدلسازی پیشین رشتههای باریک در جریان ویسکوز | ۱۷ |
| ۱-۳: اهداف و دستاوردهای پایاننامه | ۱۹ |
| فصل ۲: چارچوب محاسباتی | ۲۳ |
| ۲-۱: مکانیک تیرک | ۲۴ |
| ۲-۲: مکانیک تماس | ۳۷ |
| ۲-۳: مدل سیال | ۵۶ |
| ۲-۴: گسستهسازی زمانی و راهبرد کوپلینگ | ۶۱ |
| ۲-۵: مدلسازی دینامیک ملخ | ۶۳ |
| فصل ۳: آزمایشهای محاسباتی | ۷۱ |
| ۳-۱: حالت بررسی صفحهای برای مکانیک تیرک | ۷۱ |
| ۳-۲: راستیآزمایی اصطکاک کولمب با سطح شیبدار | ۷۴ |
| ۳-۳: گره درهمتنیده: نمایش تماس خودی اصطکاکی | ۷۶ |
| ۳-۴: رهاسازی دینامیکی تیرک پیچیده | ۷۸ |
| ۳-۵: تحلیل گرفتگی ملخ پهپاد با رشتههای منعطف | ۸۱ |
| فصل ۴: خلاصه و نتیجهگیری | ۸۷ |
| پیوست الف: استخراج جایگزین فرم ضعیف کیرشهوف-لاو | ۹۱ |
| پیوست ب: محاسبه تحلیلی انتگرالهای هسته در روش انتگرال مرزی | ۹۳ |
| مراجع | ۹۷ |
تصور کنید که دفاع یک شهر در برابر انبوهی از پهپادهای متخاصم، نه با موشک یا امواج مایکروویو، بلکه با هزاران رشتهٔ نازک و انعطافپذیر انجام شود. این دقیقاً ایدهای است که در پژوهشی نوین در مؤسسه فناوری ماساچوست (MIT) به آن پرداخته شده و آن را از یک مفهوم علمی-تخیلی به یک شبیهسازی محاسباتی دقیق تبدیل کرده است. اما چطور میتوان مطمئن شد که این رشتهها واقعاً کار میکنند؟
🔬 سیمهای نامرئی، ضربات واقعی
پژوهشگران با توسعه یک چارچوب محاسباتی پیشرفته، توانستهاند رفتار صدها رشتهٔ منعطف را درون جریان هوا شبیهسازی کنند. آنها از معادلات پیچیدهٔ مکانیک تیرک برای مدلسازی خمش و گرهخوردن رشتهها و از روش انتگرال مرزی برای شبیهسازی جریان سیال اطرافشان استفاده کردهاند. نتیجه؟ شبیهسازیهایی که نشان میدهد چگونه این رشتهها میتوانند ملخ پهپاد را درگیر کرده و با ایجاد گشتاور مخالف، سرعت چرخش آن را تا ۵۵ درصد کاهش دهند.
🧠 درهمتنیدگی هوشمندانه
جالبترین بخش این پژوهش، نقش کلیدی اصطکاک است. شبیهسازیها نشان میدهند که در غیاب اصطکاک، رشتهها به سرعت از هم باز میشوند و تأثیر خود را از دست میدهند. اما در حضور اصطکاک، گرهها پایدار میمانند و مانند یک تلهٔ نامرئی عمل میکنند. در غیاب اصطکاک، گره به سرعت باز میشود، اما در حضور اصطکاک، گره پایدار میماند و میتواند اثر مخربی بر ملخ داشته باشد
.
📊 چگالی؛ راز موفقیت
تحلیل حساسیت انجام شده نشان میدهد که دو برابر کردن چگالی رشتهها، کاهش رانش را از ۱۵٪ به ۳۸٪ میرساند. این یافتهٔ شگفتانگیز به طراحان سامانههای دفاعی میگوید که به جای افزایش تعداد رشتهها، بر افزایش جرم هر رشته تمرکز کنند. اما هشدار: رشتههای سنگینتر زودتر سقوط میکنند، پس یک محدودهٔ بهینه وجود دارد که باید کشف شود!
⚙️ از نظریه تا عمل
آنچه این پژوهش را از سایر کارها متمایز میکند، توانایی آن در شبیهسازی همزمان سه پدیدهٔ فیزیکی کاملاً متفاوت است: تغییر شکل بزرگ رشتهها، جریان سیال چسبناک اطراف آنها، و تماسهای اصطکاکی پیچیده. این سهگانهٔ محاسباتی، دریچهای به سوی طراحی سامانههای دفاعی غیرتهاجمی، قابلاستفادهمجدد و کمهزینه میگشاید که میتواند جایگزین روشهای خطرناک فعلی مانند شلیک موشک شود.