بهینه سازی مکانیزم پیوند جرثقیل

View Categories

بهینه سازی مکانیزم پیوند جرثقیل

10 min read

PDF

بهینه سازی مکانیزم پیوند جرثقیل
معرفی
در سیستم های مکانیکی پیچیده، یافتن راه حل بهینه (یا حتی به اندازه کافی خوب) تنها از طریق بینش مهندسی یا روش های آزمون و خطا می تواند چالش برانگیز باشد. استفاده از روش های بهینه سازی ریاضی می تواند مسیری کارآمد برای طراحی بهتر باشد.
در این مثال، مکانیزم پیوند در یک جرثقیل مدل‌سازی شده در ماژول Multibody Dynamics با استفاده از ماژول بهینه‌سازی بهینه شده است. هدف کاهش نیروی سیلندر مورد نیاز برای حمل محموله خاصی در بدترین حالت چرخه عملیاتی است.
این مثال ادامه مدل Truck Mounted Crane است .
توجه: این مدل به ماژول Multibody Dynamics و Module Optimization نیاز دارد.
تعریف مدل
هندسه کامل جرثقیل در شکل 1 نشان داده شده است .
شکل 1: هندسه جرثقیل
نمای نزدیک از مکانیسم پیوند مورد استفاده برای کنترل بوم داخلی در شکل 2 نشان داده شده است . برای اطلاعات بیشتر در مورد مدل جرثقیل به توضیحات مدل جرثقیل سواری کامیون مراجعه کنید .
شکل 2: نمای نزدیک مکانیسم های پیوند
یک کلید برای رنگ آمیزی قطعات تشکیل دهنده مکانیسم پیوند در جدول 1 آورده شده است .
جدول 1: شناسایی قطعات جرثقیل در مکانیزم لینک
قسمت
نام در مدل
رنگ در شکل
پایه
پایه
آبی
بوم درونی
بوم 1
سبز
سیلندر بالابر بوم
سیلندر 1
قرمز
پیستون بالابر بوم
پیستون 1
رنگ زرد
مکانیسم پیوند
لینک 1، لینک 2
سرخابی، مشکی
چرخه عملیاتی
چرخه عملیاتی به گونه ای انتخاب می شود که ممان حمل شده توسط مکانیسم پیوند تا حد امکان بزرگ باشد.
بوم داخلی از پایین ترین موقعیت خود، 45 درجه شیب رو به پایین، به بالاترین موقعیت خود (عمودی) در مراحل 15 درجه افزایش می یابد .
پسوندهای تلسکوپی تا حد امکان در همه موقعیت ها دورتر هستند.
زاویه بوم بیرونی به گونه ای انتخاب می شود که نوک جرثقیل (و در نتیجه محموله) تا حد امکان خارج شود.
بارها
وزن خود در جهت Z منفی.
محموله 1000  کیلوگرمی در نوک جرثقیل.
مشکل بهینه سازی
موقعیت سه محور را می توان تغییر داد، همانطور که در شکل 3 نشان داده شده است :
محور اتصال بازوی پیوند اول به پایه.
محور اتصال بازوی پیوند دوم به بوم.
محور اتصال دو بازوی پیوند و پیستون سیلندر هیدرولیک.
همچنین لازم است اطمینان حاصل شود که دو محور دوم بیش از 10  میلی متر به یکدیگر نزدیکتر از پیکربندی اصلی آنها حرکت نمی کنند تا از تداخل مکانیکی جلوگیری شود.
شکل 3: متغیرهای بهینه سازی (آبی) و محدودیت (قرمز).
محدودیت های جابجایی در موقعیت های محور در جدول 2 آورده شده است .
جدول 2: جابجایی های مجاز در موقعیت های محور
محور
افقی [MM]
عمودی [MM]
Base-Link1
dYL1 = [-100، +50]
dZL1 = [-150 – 30+]
Boom1-Link1
dYL2 = [-50، +50]
dZL2= [-50، +100]
لینک1-لینک2-پیستون1
dYPE = [-50، +120]
dZPE = [-100، +20]
نتایج و بحث
نیروی موجود در سیلندرهای کنترل کننده بوم در شکل 4 نشان داده شده است . با طراحی جدید، امکان کاهش 31 درصدی نیرو از 597  کیلونیوتن به 413  کیلو نیوتن وجود دارد. این یک پیشرفت اساسی است و بسته به هدف واقعی تجزیه و تحلیل می توان از آن به دو طریق استفاده کرد:
بار مجاز را می توان افزایش داد. از آنجایی که بخش قابل توجهی از ظرفیت جرثقیل برای تحمل وزن خود استفاده می شود، افزایش بار بیش از 31٪ است.
برای بار ثابت، نیروهای کاهش یافته در مکانیک پیوند دستیابی به معیارهای تنش و خستگی را برای مکانیسم پیوند آسان‌تر می‌کند.
شکل 4: تغییر نیروی سیلندر در طول چرخه بار.
به عنوان نمونه ای از تأثیر بهینه سازی بر روی سایر قسمت ها، نیروهای وارد بر محور که لولای بین پایه و بوم را تشکیل می دهد در شکل 5 نشان داده شده است .
شکل 5: نیروهای وارد بر محور بین پایه و بوم.
مقادیر واقعی متغیرهای کنترل در حالت بهینه ممکن است زمانی که این مدل در رایانه های مختلف حل می شود یکسان نباشد. دلیل آن این است که حداقل نیروی سیلندر را می توان به روش های مختلف به دست آورد. راه حل های متعددی وجود دارد که (در یک تلورانس کوچک) مقادیر بهینه یکسانی را ارائه می دهند.
نکاتی درباره پیاده سازی COMSOL
از آنجایی که تمام قطعات جرثقیل بدنه سفت و سختی هستند، نیازی به جابجایی یا تغییر شکل پیوندها در طول بهینه سازی نیست. کافی است مراکز اتصالاتی که به صورت ریاضی محورها را نشان می دهند تغییر دهید. با این حال، این منجر به یک نمایش گرافیکی غیر شهودی می شود، زیرا پیوندها دیگر حول محورها در هندسه نمی چرخند.
سیلندر و پیستون در مدل به صورت فیزیکی می چرخند، زیرا جهت نیروی سیلندر باید صحیح باشد. از نظر بصری، پیستون ممکن است به نظر برسد که گویی از سیلندر بیرون کشیده شده است. این نیز یک مصنوع گرافیکی است، زیرا اگر هندسه اصلاح شده برای طراحی نهایی انتخاب شود، استوانه دیگری انتخاب می شود.
در بهینه سازی از روش BOBYQA استفاده شده است. یک روش بدون گرادیان مناسب است، زیرا امکان محاسبه مشتقات حداکثر نیروهای سیلندر با توجه به پارامترهای بهینه سازی وجود ندارد.
مسیر کتابخانه برنامه: Multibody_Dynamics_Module/Machinery_and_Robotics/crane_link_optimization
دستورالعمل مدلسازی
کتابخانه های کاربردی
1
از منوی File ، Application  Libraries را انتخاب کنید .
2
در پنجره Application  Libraries ، Multibody  Dynamics  Module>Machinery  and  Robotics> truck_mounted_crane را در درخت انتخاب کنید .
3
 روی Open کلیک کنید .
ریشه
پارامترهای مورد نیاز برای پارامترسازی لینک را اضافه کنید.
تعاریف جهانی
پارامترهای 1
1
در پنجره Model  Builder ، در قسمت Global  Definitions روی Parameters  1 کلیک کنید .
2
در پنجره تنظیمات برای پارامترها ، بخش پارامترها را پیدا کنید .
3
 روی Load  from  File کلیک کنید .
4
به پوشه Application Libraries مدل بروید و روی فایل crane_link_optimization_parameters.txt دوبار کلیک کنید .
جزء 1 (COMP1)
سیلندر و پیستون باید جهت درستی داشته باشند تا نیروی سیلندر صحیح به دست آید.
هندسه 1
چرخش 1 (rot1)
1
در پنجره Model  Builder ، گره Component   (comp1) را گسترش دهید .
2
روی Component   (comp1)>Geometry  کلیک راست کرده و Transforms> Rotate را انتخاب کنید .
3
فقط اشیاء imp1(31) و imp1(32) را انتخاب کنید.
4
در پنجره تنظیمات برای چرخش ، قسمت Point  on  Axis  of  Rotation را پیدا کنید .
5
در قسمت متن y ، YCE را تایپ کنید .
6
در قسمت متن z ، ZCE را تایپ کنید .
7
قسمت Rotation را پیدا کنید . در قسمت Angle text cylRot را تایپ کنید .
8
از لیست نوع محور  ، x-axis را انتخاب کنید .
MULTIBODY DYNAMICS (MBD)
در پنجره Model  Builder ، گره Component   (comp1)>Multibody  Dynamics  (mbd) را گسترش دهید .
پایه لولا-لینک1
مراکز اتصال را برای اتصالات قابل تنظیم پارامتر کنید.
1
در پنجره Model  Builder ، گره Component   (comp1)>Multibody  Dynamics  (mbd)>Hinge  Base-Link1 را گسترش دهید ، سپس بر روی Hinge  Base-Link1 کلیک کنید .
2
در پنجره تنظیمات برای مفصل لولا  ، قسمت مرکز اتصال را پیدا کنید .
3
از لیست، User  defined را انتخاب کنید .
4
بردار X c را به صورت مشخص کنید
 
XC
ایکس
mYL1
y
mZL1
z
لولا بوم1-لینک2
1
در پنجره Model  Builder ، روی Hinge  Boom1-Link2 کلیک کنید .
2
در پنجره تنظیمات برای مفصل لولا  ، قسمت مرکز اتصال را پیدا کنید .
3
از لیست، User  defined را انتخاب کنید .
4
بردار X c را به صورت مشخص کنید
 
XC
ایکس
mYL2
y
mZL2
z
اسلات Link1-Link2
1
در پنجره Model  Builder ، روی Slot  Link1-Link2 کلیک کنید .
2
در پنجره تنظیمات برای Slot  Joint ، بخش مرکز  اتصال  را پیدا کنید .
3
از لیست، User  defined را انتخاب کنید .
4
بردار X c را به صورت مشخص کنید
 
XC
ایکس
مایپ
y
mZPE
z
اسلات Link1-Piston1
1
در پنجره Model  Builder ، روی Slot  Link1-Piston1 کلیک کنید .
2
در پنجره تنظیمات برای Slot  Joint ، بخش مرکز  اتصال  را پیدا کنید .
3
از لیست، User  defined را انتخاب کنید .
4
بردار X c را به صورت مشخص کنید
 
XC
ایکس
مایپ
y
mZPE
z
تعاریف جهانی
پارامترهای 1
1
در پنجره Model  Builder ، در قسمت Global  Definitions روی Parameters  1 کلیک کنید .
2
در پنجره تنظیمات برای پارامترها ، بخش پارامترها را پیدا کنید .
3
در جدول تنظیمات زیر را وارد کنید:
 
نام
اصطلاح
ارزش
شرح
ExtLen
5.5 [m]
5.5 متر
طول کل پسوند
مطالعه 1
چرخه بار را به بدترین سناریو تغییر دهید و محاسبه کنید.
مرحله 1: ثابت
1
در پنجره Model  Builder ، گره Study  1 را گسترش دهید ، سپس روی Step  1:  Stationary کلیک کنید .
2
در پنجره تنظیمات برای Stationary ، بخش Study  Extensions را پیدا کنید .
3
در جدول تنظیمات زیر را وارد کنید:
 
نام پارامتر
لیست مقادیر پارامتر
واحد پارامتر
Angle1 (زاویه به افقی، بوم داخلی)
محدوده (-45،15،90)
راد
RelAng (زاویه بین بوم ها)
محدوده 0 0 0 (0،-15،-90)
راد
4
برای انتخاب ردیف شماره 3 در جدول کلیک کنید.
5
 روی Delete کلیک کنید .
6
در نوار ابزار صفحه اصلی ،  روی محاسبه کلیک کنید .
راه حل را برای پیکربندی اصلی ذخیره کنید.
راه حل 1 (sol1)
در پنجره Model  Builder ، در قسمت Study  1>Solver  Configurations روی Solution   (sol1) کلیک راست کرده و گزینه Solution>Copy را انتخاب کنید .
نتایج
جهانی 1
1
در پنجره Model  Builder ، گره Results>Boom  Cylinder  Forces را گسترش دهید ، سپس روی Global  1 کلیک کنید .
2
در پنجره تنظیمات برای جهانی ، بخش y-Axis  Data را پیدا کنید .
3
برای انتخاب ردیف شماره 2 در جدول کلیک کنید.
4
 روی Delete کلیک کنید .
5
در جدول تنظیمات زیر را وارد کنید:
 
اصطلاح
واحد
شرح
abs(comp1.mbd.prj1.Fl1)
kN
بهینه شده است
جهانی 2
1
روی Results>Boom  Cylinder  Forces>Global  کلیک راست کرده و Duplicate را انتخاب کنید .
2
در پنجره تنظیمات برای جهانی ، بخش داده را پیدا کنید .
3
از فهرست مجموعه داده ، مطالعه  1/راه حل   –  کپی   (sol2) را انتخاب کنید .
4
قسمت y-Axis  Data را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
اصطلاح
واحد
شرح
abs(comp1.mbd.prj1.Fl1)
kN
اصلی
نیروهای سیلندر اکستنشن
در پنجره Model  Builder ، در قسمت Results روی Extension  Cylinder  Forces کلیک راست کرده و Delete را انتخاب کنید .
جهانی 1
1
در پنجره Model  Builder ، گره Results>Hinge  Forces را گسترش دهید ، سپس روی Global  1 کلیک کنید .
2
در پنجره تنظیمات برای جهانی ، بخش y-Axis  Data را پیدا کنید .
3
برای انتخاب ردیف شماره 3 در جدول کلیک کنید.
4
 دوبار روی Delete کلیک کنید  .
5
در جدول تنظیمات زیر را وارد کنید:
 
اصطلاح
واحد
شرح
mbd.hgj1.Fy
kN
نیروی بهینه، y
mbd.hgj1.Fz
kN
نیروی بهینه، z
sqrt(mbd.hgj1.Fy^2+mbd.hgj1.Fz^2)
kN
نیروی کل بهینه شده
جهانی 2
1
روی Results>Hinge  Forces>Global  کلیک راست کرده و Duplicate را انتخاب کنید .
2
در پنجره تنظیمات برای جهانی ، بخش داده را پیدا کنید .
3
از فهرست مجموعه داده ، مطالعه  1/راه حل   –  کپی   (sol2) را انتخاب کنید .
4
قسمت y-Axis  Data را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
اصطلاح
واحد
شرح
mbd.hgj1.Fy
kN
نیروی اصلی، y
mbd.hgj1.Fz
kN
نیروی اصلی، z
sqrt(mbd.hgj1.Fy^2+mbd.hgj1.Fz^2)
kN
نیروی کل اصلی
مطالعه بهینه سازی را اضافه کنید.
مطالعه 1
بهينه سازي
1
در نوار ابزار مطالعه ، روی  Optimization کلیک کنید و Optimization را انتخاب کنید .
2
در پنجره تنظیمات برای بهینه سازی ، بخش Optimization  Solver را پیدا کنید .
3
از لیست روش ، BOBYQA را انتخاب کنید .
4
در قسمت متنی Optimality  tolerance ، 0.1 را تایپ کنید .
5
قسمت Objective  Function را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
اصطلاح
شرح
ارزیابی برای
abs(comp1.mbd.prj1.Fl1)
نیروی سیلندر (داخلی)
ثابت
6
از لیست راه حل ، حداکثر  اهداف  را انتخاب کنید .
7
قسمت Control  Variables  and  Parameters را پیدا کنید .  روی Load  from  File کلیک کنید .
8
به پوشه Application Libraries مدل بروید و روی فایل crane_link_optimization_ctrlvars.txt دوبار کلیک کنید .
9
قسمت Constraints را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
اصطلاح
کران پایین
کران بالا
ارزیابی برای
mdL2PE-dL2PE
-0.01
ثابت
10
قسمت Output  while  Solving را پیدا کنید . کادر Plot را انتخاب کنید .
11
از لیست گروه Plot  ، Boom Cylinder Forces را انتخاب کنید .
در صورتی که حالتی قابل حل نباشد، از توقف حلقه بهینه سازی خودداری کنید.
راه حل 1 (sol1)
1
در پنجره Model  Builder ، در قسمت Study  1>Solver  Configurations>Solution   (sol1)>Stationary  Solver  1 روی Parametric  1 کلیک کنید .
2
در پنجره تنظیمات برای پارامتر ، بخش عمومی را بیابید .
3
از لیست خطای On  ، گزینه Store vala solution را انتخاب کنید .
بهینه سازی را اجرا کنید.
4
در نوار ابزار مطالعه ،  روی محاسبه کلیک کنید .
نتایج
نیروهای سیلندر بوم
1
 روی دکمه Zoom  Extents در نوار ابزار Graphics کلیک کنید .
2
در پنجره Model  Builder ، در قسمت Results روی Boom  Cylinder  Forces کلیک کنید .
3
در نوار ابزار Boom Cylinder Forces ، روی  Plot کلیک کنید .
نیروهای لولا
1
 روی دکمه Zoom  Extents در نوار ابزار Graphics کلیک کنید .
2
در پنجره Model  Builder ، روی Hinge  Forces کلیک کنید .
3
در نوار ابزار Hinge Forces ، روی  Plot کلیک کنید .