فهرست مطالب
مجله مهندسی مکانیک شریف
سال بیست و نهم شماره 1 (بهار و تابستان 1392)
- تاریخ انتشار: 1392/08/19
- تعداد عناوین: 14
-
-
صفحات 3-10در این نوشتار از روش تجزیه ی متعامد بهینه(P O D) برای تخمین سریع میدان جریان، محاسبه ی ضرایب و طراحی معکوس آیرودینامیکی استفاده شده است. در این رابطه دو روش برای تخمین جریان تراکم پذیر غیرلزج در رژیم های جریان متفاوت مورد ستفاده قرار گرفته که هر دو روش ترکیبی از روش تجزیه ی متعامد بهینه است. روش های رتبه کاسته ی حاصل می توانند با دقتی مناسب و سرعتی بالا تغییرات در میدان جریان ناشی از تغییر در پارامترهای موثر ازجمله عدد ماخ و زاویه ی حمله را پیش بینی کنند. برای طراحی معکوس از یک مدل ترکیبی مبتنی بر روش تجزیه ی متعامد بهینه به همراه یک مسئله ی بهینه سازی برای توزیع فشار سطحی یا کمینه سازی ضریب پسا استفاده شده است. نتایج حاصله با شبیه سازی های عددی مستقیم (D N S) و نتایج حل دقیق به عنوان مرجع مقایسه شده که دقت و سهولت چشمگیری را در محاسبات نشان می دهد.
کلیدواژگان: روش تجزیه ی متعامد بهینه، آیرودینامیک تراکم پذیر، طراحی معکوس، روش رتبه کاسته -
صفحات 11-21در این نوشتار روشی تحلیلی برای محاسبه ی سفتی بریدد کامپوزیت های سه بعدی چهارجهته ارائه شده است. در بیشتر تحقیقات قبلی از روش اجزاء محدود استفاده شده و کم تر روش تحلیلی به کار رفته است. در روش تحلیلی ارائه شده در این تحقیق کل کامپوزیت به چهار نوع سلول واحد ناحیه ی درونی، سطحی درونی، سطحی بیرونی و گوشه یی تقسیم شده که هریک خواص مکانیکی و هندسی خاص خود را دارند. سپس هرکدام از این نواحی به صورت یک کامپوزیت تک جهته در نظر گرفته شده و با استفاده از یک روش میانگین گیری حجمی، سفتی کل بریدد کامپوزیت محاسبه شده است. در نهایت، نتایج حاصله با نتایج به دست آمده از روش آزمایشگاهی مقایسه و تاثیر زاویه ی بریدینگ بر خواص مکانیکی بریدد کامپوزیت بررسی شده است.
کلیدواژگان: بریدد کامپوزیت سه بعدی چهارجهته، محاسبه ی سفتی، روش تحلیلی، روش چند سلول واحد -
صفحات 23-30در این مقاله از روش المان طیفی برای تحلیل انتشار موج در ورق کامپوزیتی استفاده می شود. ویژگی اصلی این روش قطری بودن ماتریس جرم است کهباعث می شود زمان تحلیل در مقایسه با روش المان محدود معمولی بسیار کوتاه تر شود. شبیه سازی برای المان های 9، 25 و 36 گره، با قید رسیدن به هم گرایی 2 درصد انجام می شود. از مقایسه ی نتایج نهایی (سرعت گروهی) به دست آمده با حل تحلیلی درخواهیم یافت که با افزایش تعداد گره های المان، در برخی حالات نتایج اندکی بهبود می یابد ولی در ازاء این افزایش دقت، زمان حل افزایش قابل توجهی پیدا می کند. بنابراین در این مسئله، استفاده از المان هایی با تعداد گره از نظر زمان محاسباتی ارجح است، اگرچه با افزایش تعداد گره های المان دقت تقریب تنش و کرنش بیشتر می شود.
کلیدواژگان: انتشار موج، ورق کامپوزیتی، روش المان طیفی، دقت محاسباتی -
صفحات 31-40هدف این پژوهش بررسی تاثیر ریزساختاری فولاد و اندازه دانه روی زبری سطح به دست آمده پس از فرایند سنگ زنی است. برای این کار نمونه هایی با ریزساختارهای مختلف و با اندازه دانه های مختلف طی فرایندهای مختلف عملیات حرارتی آماده سازی شده اند. عملیات سنگ زنی با شرایط ثابت و مشخص روی قطعات اعمال شده است. پس از اتمام فرایند، زبری سطح کل قطعات اندازه گیری شد. مشاهده شد زبری با ریزساختار نمونه ها و سختی مربوط به آن ها رابطه دارد، به طوری که با افزایش سختی مربوط به هر ریزساختار، زبری سطح کاهش می یابد. این روند تا یک حد بحرانی ادامه یافته و پس از آن زبری رو به افزایش می گذارد. در این بین صافی سطح نمونه های ریز دانه وضعیت بهتری نسبت به نمونه های درشت دانه دارند. همچنین سنگ زنی روی نمونه های نورد شده نیز انجام گرفت. مشاهده شد سنگ زنی در جهت عمود بر دانه های کشیده شده صافی سطح بهتری را به دنبال دارد.
کلیدواژگان: سنگ زنی، زبری سطح، ریزساختار، اندازه دانه -
صفحات 41-52در این نوشتار تحلیل اجزاء محدود تغییر شکل موضعی در میله ی تحت کشش و همچنین تغییر شکل کرنش صفحه یی ورق به وسیله ی دو سنبه ی تخت و متقابل براساس الگوریتم A L E در نرم افزار 6.7 A B A Q U S بررسی شده است. فرایند تغییر شکل کرنش صفحه یی ورق به وسیله ی دو سنبه ی تخت و متقابل با سطوح تماس بدون اصطکاک شبیه سازی شده است. در آغاز این فرایند بدون در نظر گرفتن اثرات نرخ کرنش شبیه سازی شده و سپس با در نظر گرفتن یک ماده وابسته به نرخ کرنش برای ورق، اثرات نرخ کرنش بر چگونگی شکل دهی ورق لحاظ شده است. برای نشان دادن قابلیت الگوریتم A L E در شبیه سازی مسائل با تغییر شکل های بزرگ، نتایج به دست آمده از این روش با نتایج روش لاگرانژی مقایسه و برتری این روش نشان داده شده است. نتایج این شبیه سازی با نتایج موجود در مقالات مقایسه شده که مطابقت خوبی مشاهده می شود.
کلیدواژگان: روش اجزاء محدود، الگوریتم A L E، تغییر شکل موضعی در میله تحت کشش، تغییر شکل کرنش صفحه یی ورق به وسیله ی دو سنبه ی تخت و متقابل -
صفحات 53-59با استفاده از شبیه سازی اجزاء محدود می توان آنالیز جریان ماده در طول فرایند فورج، نیروی فورجینگ در هر زمان، توزیع فشار تماسی بین قالب و قطعه کار، و نیز خیز کشسان قالب های فورج را پیش بینی کرد.
بنابراین برای جبران مقدار تغییر شکل کشسان قالب باید شکل قالب را با استفاده از تغییر اندازه های آن اصلاح کرد. در این تحقیق با استفاده از شبیه سازی اجزاء محدود سه بعدی، خیز کشسانی قالب در مقاطع مختلف پره به دست می آید و با اصلاح پروفیل های قالب بالایی و پایینی در این مقاطع و در جهت های مخالف آن، خطاهای ناشی از تغییر شکل کشسان قالب جبران می شود. با استفاده از این اصلاحات، خطاهای ناشی از تغییر شکل کشسان به مقدار قابل ملاحظه یی کاهش خواهد یافت. برای جبران تغییر شکل کشسان قالب باید چندبار شبیه سازی انجام داد تا خطای ناشی از این تغییر شکل کمینه شود.
کلیدواژگان: فورج پره توربین گازی، تغییر شکل کشسان قالب، شبیه سازی اجزاء محدود سه بعدی -
صفحات 61-68با افزوده شدن اثرات گرانش به معادلات جریان امتزاج ناپذیر در محیط های متخلخل موسوم به باکلی لورت، مشکلاتی در حل عددی پدید می آید. غیرمحدب بودن توابع شار اجزاء در معادلات باکلی لورت، حل دستگاه معادلات هذلولوی مربوط به پایستگی جرم را مشکل می کند. درصورت در نظر گرفتن اثرات گرانش، این محدب نبودن شدیدتر می شود به نحوی که با به کارگیری روش های متداول برای حل صریح دستگاه معادلات هذلولوی شوک های انبساطی غیرفیزیکی در پاسخ ها ایجاد خواهد شد. در مطالعه ی حاضر از روش عددی موج غالب برای رفع این مشکل در شبیه سازی جریان امتزاج ناپذیر تحت اثر جاذبه در مخازن نفتی استفاده می شود. این روش نیاز به تجزیه ی مشخصه یی ندارد و در عین حال میزان پخش عددی آن بسیار کم تر از روش های مرکزی متداول است. همچنین از تصحیح آنتروپی هارتن برای جلوگیری از ایجاد شوک های انبساطی استفاده می شود. در این نوشتار اثر گرانش در یک مخزن همگن یک بعدی با استفاده از مدل باکلی لورت شبیه سازی، و از راوابط مختلف برای نفوذپذیری استفاده می شود.
کلیدواژگان: اثرات گرانشی، معادلات باکلی، لورت، نفوذپذیری، قوانین پایستگی، روش موج غالب، تصحیح آنتروپی -
صفحات 69-76به منظور بررسی اثر زاویه ی پس گرا بر توزیع فشار سطح بالایی یک بال مخصوص، آزمایش های گوناگونی در رژیم جریان زیر صوت انجام گرفت.این آزمایشات در اعداد رینولدز و زوایای حمله ی مختلف، بر روی سه نیم مدل بال با نسبت منظری یکسان و زوایای پس گرای مختلف صورت پذیرفت. نتایج نشان می دهد که توزیع فشار سطح بالای بال به دلیل وجود جریان های عرضی ناشی از زاویه ی پس گرای بال به ایجاد روند متفاوتی نسبت به حالت دوبعدی می انجامد. محدود بودن بال سبب ایجاد گردابه هایی در نوک بال شده که باعث تقویت اثرات سه بعدی جریان روی بال می شود، به نحوی که موقعیت و مقدار افت فشار روی آن در نزدیکی نوک بال نسبت به نقاط دیگر بال متفاوت خواهد بود و کارایی ایرفویل در مقاطع نزدیک به نوک بال بیشتر از مقاطع میانی دستخوش تغییرات شده و از حالت بهینه خارج شده است. در زوایای پس گرای کم، توزیع فشار روی بال نسبت به تغییرات عدد رینولدز حساسیت کم تری دارد.
کلیدواژگان: بال پس گرا، زاویه ی پس گرایی، میدان جریان، توزیع فشار -
صفحات 77-85در این نوشتار اختلاط ناشی از جریان الکترواسموتیک درون ریزمجراهای دارای خواص سطحی ناهمگن مورد مطالعه قرار گرفته است. میدان سرعت توسط حل عددی معادله های ناویر استوکس به همراه شرط سرعت لغزشی دردیواره به دست آمده، و با استفاده از آن میدان غلظت گونه ها محاسبه شده است. سپس معیار مناسب برای تعیین کارایی اختلاط جریان درون ریزمجراهای ناهمگن معرفی شده و با استفاده از آن، عملکرد یک تکه ناهمگنی واقع بر دیواره ی ریزمجرا در ایجاد اختلاط بررسی شده و اثر پارامترهای هندسی و چگونگی پروفیل ناهمگنی مورد مطالعه قرار گرفته است. بررسی عددی نشان داده که وجود ناپیوستگی در پروفیل ناهمگنی تفاوت چندانی در عملکرد اختلاطی ریزمجرا ندارد. نتایج نیز نشان می دهد که ایجاد ناهمگنی با اندازه های کوچک تر در نزدیکی ورودی ریزمجرا بیشترین بهبود را در اختلاط ایجاد می کند. این نتایج می تواند در طراحی بهینه ی ریزمخلوط گرها مورد استفاده قرار بگیرد.
کلیدواژگان: ریزمجرای ناهمگن، جریان الکترواسموتیک، کارایی اختلاط -
صفحات 87-92هدف اصلی در این نوشتار محاسبه ی ضرایب آئرودینامیکی بال در جریان مافوق صوت دائم برمبنای حل جریان پتانسیل با استفاده از روش اجزای مرزی است. در این راستا ابتدا معادله ی حاکم بر جریان پتانسیل مافوق صوت استخراج و خطی سازی شده، و سپس با استفاده از روش اجزای مرزی و اعمال توابع شکل مرتبه ی اول، دستگاه معادلات گسسته برای محاسبات جریان تشکیل شده است. در خلال حل این معادلات مشکلات محاسباتی مربوط به تعیین ضرایب انتگرالی تکین بروز می کند که ناشی از قطع مخروط ماخ توسط اجزای مرزی موجود بر دامنه ی محاسباتی است. البته با استفاده از روش های عددی می توان مشکلات مربوطه را برطرف کرد. در نهایت برای دو بال با هندسه ی مشخص، محاسبات جریان مافوق صوت دائم انجام شده است. مقایسه ی نتایج مزبور با جواب های دقیق، بیان گر دقت خوب روش حاضر و انعطاف پذیری آن برای تحلیل جریان حول بال های سه بعدی است.
کلیدواژگان: آئرودینامیک، بال، روش اجزای مرزی، جریان پتانسیل، مافوق صوت -
صفحات 93-101پیچیدگی معادلات حاکم بر بازوهای ماهر رباتیکی سبب شده است روش های متنوعی برای کنترل آن ها ارائه شود. مشخصه های دینامیکی بازوها در طراحی این روش ها معلوم و معین در نظر گرفته می شوند. این در حالی است که مشخصه های یادشده در بیشتر موارد به دلایلی چون نواقص ساخت، استهلاک یا طبیعت خود موضوع (مثلا اصطکاک) به طور دقیق در دست نیستند. از این رو در این نوشتار برآنیم تا ضمن ارائه ی یک مدل دینامیکی ساده برای ربات و شبیه سازی آن توسط نرم افزار M A T L A B، روشی مبتنی بر کم ترین مربعات برای شناسایی پارامترهای یک بازوی ماهر صفحه یی با چهار درجه ی آزادی به دست آوریم. مدل دینامیکی ارائه شده شامل کم ترین پارامترهای مستقل در معادلات حاکم بر حرکت خواهد بود. روش ارائه شده برای شناسایی صحیح پارامترهای مجهول، کمینه کردن مجموع درایه های قطر اصلی ماتریس هم پراش خروجی های آلوده به اغتشاش است. نتایج به دست آمده از شبیه سازی تخمین بسیار خوبی از مقادیر مجهول برآورد می کند.
کلیدواژگان: بازوهای ماهر رباتیکی، شناسایی سیستم، کم ترین مربعات خطا، دینامیک چندبندی ها، دینامیک معکوس، واسنجی -
صفحات 103-111
در این تحقیق، هدف ارائه ی پیشنهادهای کاربردی مناسب با ارزیابی تاثیر جانمایی وسیله ی احتراقی و دریچه ی خروجی بر کیفیت هوا در فضاهای دارای وسایل گرمایشی احتراقی، با دودکش و بدون آن است. بدین منظور، از شناور حرارتی ایجادشده توسط وسیله ی احتراقی به عنوان روشی موثر در کنترل غلظت آلاینده ها استفاده می شود. برای حل مسئله از شبیه سازی و حل عددی استفاده شده است. طبق نتایج حاصله، نصب بخاری در کنار دیوار خارجی، روبه روی دریچه ی خروجی و در گوشه های اتاق و نصب دریچه ی خروجی در مرکز دیوار خارجی، به عنوان انتخاب بهینه پیشنهاد می شوند. مقایسه ی بهترین و بدترین وضعیت اتاق از نظر تاثیر هم زمان جانمایی بخاری و دریچه ی خروجی در کیفیت هوا نشان می دهد که غلظت C O$_2$ در ناحیه ی تنفس اتاق در مدل بهینه، 9٪ کم تر از بدترین وضعیت است.
کلیدواژگان: کیفیت هوای داخل، شناور حرارتی، وسایل احتراقی -
صفحات 113-124یکی از متداول ترین تکنیک های اسکن سه بعدی اجسام و بازیابی مدل سه بعدی آن ها روش مثلث بندی با استفاده از نور ساختاریافته است. در این نوشتار برای انجام عملیات اسکن سه بعدی با استفاده از نور لیزر، از توزیع گوسی استفاده شده است. سپس برای اندازه گیری داده های عمق، یک دوربین C C D تصاویر را می گیرد. ابتدا عملیات بهبوددهی و پردازش تصویر دریافتی در دو بخش مهم انجام می شود: 1. پارامترهای مربوط به جسم و شرایط نورپردازی مثل رنگ و لبه یابی، کنترل و اصلاح می شود؛ 2. پارامترهای مربوط به کالیبراسیون دوربین مثل انحراف محوری دوربین و ایجاد تصویری با اختلاف نسبت به محور افقی، به صورت نرم افزاری و با دقت خوبی اصلاح می شود. بعد از پردازش تصویر، جسم اسکن می شود و پروفیل دوبعدی آن استخراج می گردد. با گرفتن تصاویر متوالی از جسم، مدل سه بعدی به صورت مدل C A D بازیابی و انجام تحلیل های مختلف روی مدل جسم امکان پذیر می شود.
کلیدواژگان: اسکن سه بعدی، روش مثلث بندی، مدل سه بعدی، بازیابی مدل -
صفحات 125-134یکی از متداول ترین تکنیک های اسکن سه بعدی اجسام و بازیابی مدل سه بعدی آن ها روش مثلث بندی با استفاده از نور ساختاریافته است. در این نوشتار برای انجام عملیات اسکن سه بعدی با استفاده از نور لیزر، از توزیع گوسی استفاده شده است. سپس برای اندازه گیری داده های عمق، یک دوربین C C D تصاویر را می گیرد. ابتدا عملیات بهبوددهی و پردازش تصویر دریافتی در دو بخش مهم انجام می شود: 1. پارامترهای مربوط به جسم و شرایط نورپردازی مثل رنگ و لبه یابی، کنترل و اصلاح می شود؛ 2. پارامترهای مربوط به کالیبراسیون دوربین مثل انحراف محوری دوربین و ایجاد تصویری با اختلاف نسبت به محور افقی، به صورت نرم افزاری و با دقت خوبی اصلاح می شود. بعد از پردازش تصویر، جسم اسکن می شود و پروفیل دوبعدی آن استخراج می گردد. با گرفتن تصاویر متوالی از جسم، مدل سه بعدی به صورت مدل C A D بازیابی و انجام تحلیل های مختلف روی مدل جسم امکان پذیر می شود.
کلیدواژگان: اسکن سه بعدی، روش مثلث بندی، مدل سه بعدی، بازیابی مدل
-
Pages 3-10Lately, scientific and engineering computations have been receiving detailed attention, in consequence of the improvement in computer hardware. Faster processors, efficient memory modules, and parallel computation techniques have prepared an adequate base for simulation of complex physical phenomena. However, the challenges in performing a simulation of an engineering problem and analyzing the attainable data increase manifold with the level of difficulty in the governing equations. Despite this fact, attempts to challenge these difficulties are inadequate. Numerical simulations yield large sets of data, which may not be useful as an experimental assignment. These data sets not only provide good and deep intuition into the physics of the problem, but also simultaneously assist the assessment and enhancement of various models that are developed for the computations. In this article, Proper Orthogonal Decomposition (POD) has been used for fast data estimation of flow the field, calculation of aerodynamical coefficients, and inverse aerodynamic design. In this way, two methods have been used for estimation of inviscid compressible flow fields in various regimes, in which both methods are based on combined form of POD. The first extension is a coupling between the POD method and a cubic spline interpolation, as introduced for inviscid flows. The second is essentially a new technique, which has been introduced by the authors for first time. In this method, some additional calibrations (including a kind of re-projection) are needed to achieve more accurate estimations. The outcome of these low order procedures can predict variations of flow field, due to variations of some effective parameters (e.g., Mach number or angle of attack) with high accuracy and high speed of computation. For the inverse design, we used a combined model, based on POD with the solution of an optimization problem, to obtain the desired surface pressure distribution or to minimize drag coefficient. In this way, two sets of snapshots have been used; an ensemble of flow fields, which may be pressure distribution or any aerodynamic characteristics, and a geometrical ensemble. The low order POD technique is used for computations of relative POD modes and modal coefficients of these ensembles. To solve the coupled optimization problem, an iterative approach has been used to achieve the desired accuracy. Note that the present inverse design approach can be reconstructed for different forms of the objective function. An order reduction manner has been used to find the optimal number of modes for data reconstruction that is a higher resolution of reconstruction which fewer modes for data reconstruction. The obtained results are compared to CFD simulations or exact solution results as reference data. The results show the relatively good accuracy and also the simplicity of the procedure.Keywords: Proper orthogonal decomposition, Compressible Aerodynamics, Inverse Design, Order Reduction Manner
-
Pages 11-21The present paper represents a new analytical method for calculation of the stiffness of three-dimensional, four-directional braided composites. In most previous work, the analytical model has been largely neglected in favor of the finite element model. Among those who have used the analytical model, the braided preform has been considered to be made of one, or three, types of representative unit cells, while microscopic evidence of the microstructure of pre-forms reveals different configurations of the yarns in the interior, surface and corner regions of a braided pre-form. This paper presents a multi-unit cell model, in which four kinds of unit cell, namely; interior, interior surface, exterior surface and corner unit, have been introduced as representative cells. Each type of unit cell in a braided composite has different micro-structural parameters and possesses unique mechanical properties. In this model, each type of unit cell has been considered a uni-directional composite. Next, using rotation matrices, the angle between yarns and longitudinal direction has been considered in general coordinates of the model. Finally, using a volume averaging method, the total stiffness of the braided composite is calculated. The results are compared with those obtained from experimental methods and the effect of braiding angle on the stiffness of braided composites is discussed.Keywords: 3, Dimensional four, directional braided composites, stiffness calculation, analytical method, multi unit cell model
-
Pages 23-30The present paper deals with transverse elastic wave propagation using spectral element method (SEM), based on the first order shear deformation plate theory. In this method, the element nodes are located at the zeros of an orthogonal function, such as the Legendre function, which causes the mass matrix to be diagonal. On the other hand, the differential method is implemented for solving the problem of limited intermediate. Dealing with a diagonal mass matrix is the main advantage of the present approach, because the analysis only requires the inverse of a diagonal matrix. Therefore, the analysis time is much less than the one with the conventional Finite Element Method. In this study, a 10 layered square composite plate with 1 m side and 10 mm thickness is investigated. A time-variable force is applied at the centre of the plate during 120 msec, with 25 kH frequency. It is assumed that the fibers are aligned parallel to the side of the plate at all layers. Elements with 9, 25 and 36 nodes are used for simulations. In each study case, the total number of elements of the plate is selected in such a way by increasing the number of elements and taking into account a 2% difference for two consecutive steps. Therefore, the plate is divided into $200times200$ elements (40000 elements) with 9 node elements. The selection of $100times100$ and $60times60$ elements for 25 and 36-noded elements, respectively, result only in 2 percent reduction in convergence time. Having solved the problem, three-dimensional deflection and two-dimensional deflection contour on paths passing through the center of the plate and parallel to the sides have been obtained. Then, using the two-dimensional deflection results, group wave velocity values? ?in the longitudinal (fiber direction) and transverse directions (perpendicular to fibers) of the plate have been obtained. Creating the element stiffness and mass matrices, and assembling these matrices, makes the whole process analysis time. Therefore, for comparison between the elements with 9, 25 and 36 nodes, the time required for matrix assembly are registered. Besides, using theoretical formulations, the value of the group velocity for the desired frequency in this direction has been obtained and compared with those obtained via finite element analysis. A three-dimensional deflection contour has shown that the wave front is elliptical and a large diameter of oval is located in the fiber direction. In addition, the value of the maximum deflection of the plates has been decreased with increasing time. As shown by our findings, by increasing the number of element nodes, the time required for assembling the stiffness and mass matrices can be accelerated significantly. In addition, generally, by increasing the number of nodes per element, the differences between these results with those obtained from the analytical method are higher, and only in some cases lead to better correlations. Consequently, slightly better correlations incur very high computational costs. It is noted that elements with fewer numbers of nodes have lower degree polynomials, and may not establish the continuity of stress and strain. For the 9-node element used, for example, with a function of order two, the curvature function is obtained by double differentiation. Since the curvature is constant, the adjacent elements cannot have the same curvature. However, this must be the same for all elements. Therefore, from one element to the adjacent element, there will be a discontinuity in the stress and strain distribution. However, this is not the case for higher-order elements. In conclusion, if only the total amount of computational time and volume is required, a 9-node element is preferred. On the other hand, if the continuity condition of stress and strain distribution is preferred, then the advantage is in using higher-order elements.Keywords: Wave propagation, composite plate, spectral element method (SEM), simulation accuracy
-
Pages 31-40This work outlines the effects of microstructure and grain size on surface roughness during the grinding process of carbon steel AISI 1045. According to this study, it is evinced that the obtained surface roughness after the grinding process is concerned with the microstructure and hardness of the initial surface. Different samples with different microstructures and three various levels of grain size (fine, median and coarse) have been prepared through specific heat treatment. After finishing the grinding operation, the surface roughness (Ra and Rz) was measured. A considerable relationship was found between the surface roughness and microstructures of samples. Furthermore, there was a relation between surface roughness and initial hardness. An increase in initial surface hardness of the material led to a decrease in surface roughness, but this relation in super hard materials was different. Surface roughness decreased after a critical initial hardness. In fact, surface roughness rose by increasing initial hardness, and, after that, there was a drastic reduction. Meanwhile, surface roughness was investigated for fine, median and coarse grain materials. The amount of surface roughness in fine grain materials was more than coarse grain materials. The surface roughness was improved by an increase in the grain size of microstructures. Another factor, which is focused upon in this paper, is an investigation of the relationship between the direction of the grinding operation and the orientation of grains in rolled materials. It was found that grinding in the same direction as elongated grains resulted in improved surface roughness, and grinding perpendicular to the elongated grains led to a decrease in surface roughness. As a result, it can be achievable to acquire a desired surface roughness by identifying microstructure type, initial hardness and grain size of materials correctly. Furthermore, it is important to note the grain direction of rolled material during the grinding process. Finally, one must be careful to observe that achieving ideal surface roughness does not result in a serious change in microstructure.Keywords: Grinding, surface roughness, hardness, microstructure, grain size
-
Pages 41-52In the present study, the local deformation of a cylindrical bar under uniaxial extension, and, also, the plane strain deformation of shell with two flat and reciprocal punches, were simulated via the decoupled ALE (Arbitrary Lagrangian Eulerian) finite element formulation in ABAQUS/Explicit code. In the ALE formulation, a reference computational domain is introduced and the finite element mesh is neither attached to the material nor fixed in space, so that the motion of the mesh is independent of that of the material. So, the mesh is managed by a procedure that tries to minimize mesh distortions and/or concentrate the mesh in a particular region. As a result, the ALE formulation can handle path-dependent material behavior and free surface conditions, while maintaining mesh fineness. An operator split procedure, a common and efficient choice in ALE analysis, is used to treat the independent mesh-material motion. Each time increment is divided into a Lagrangian phase considering material effects, which is done with the explicit time integration formulation, followed by an Eulerian phase, where the mesh smoothing is occurred, subsequently, taking into account the convective effects by solving the convection equation. The second-order Van Leer algorithm is employed to solve the convection equation. The mesh motion strategy is performed according to the mesh spatial constraints and volume smoothing method compatible with the problem type. For verifying the ALE finite element model, the local deformation of a cylindrical bar under uniaxial extension was simulated. The results were compared with reported ALE results in literature and good agreement was observed. In addition, the merits of the method for large deformation problems were discussed and the results were compared with those predicted by purely Lagrangian formulations. Although both results were almost the same, a remarkable decrease in simulation time was achieved using the ALE method.Frictionless contact is postulated in the analysis of the plane strain deformation of shell with two flat and reciprocal punches. Firstly, the problem was solved without considering strain rate effects. However, strain rate effects and dynamic effects due to high punch velocity were also investigated for a strain rate dependent material. It has been demonstrated that the inertia effects increase when the punch velocity becomes larger, resulting in different final shapes of the work-piece from those with low punch velocity. From the analyzed plane strain problem using an uncoupled ALE approach, the predicted punch load values and other comparison parameters are compatible with the coupled ones mentioned earlier in the literature.Keywords: FEM, ALE formulation, necking, punch indentation
-
Pages 53-59Manufacturing of engineering components with precise dimensions are becoming quite important in attempts to reduce cost and improve reliability. In forging die design, dimensional accuracy is one of the main goals. The load carrying capacity and life of any forged product is greatly affected by its dimensional accuracy. To predict the precise dimension of the part and determine the die dimension for precision forging, it is necessary to analyze the factors which affect dimensional accuracy. For manufacturing a precision forging of turbine blades, thorough recognition of prevalent parameters is essential, such as preform design, process simulation and various techniques for compensation of errors due to die-elasticity characteristics. The accuracy of aerodynamic cross-sections depends on many factors such as die elasticity, heat distortion and deformation of parts during the cooling of the die. The main goal of the design of aerodynamic section for the blade and its precision forging is to minimize the machining process. By using finite element simulation; material flow analysis during the forging process, forging force in any moment, contact pressure distribution between die and workpiece, and also elastic deflection of forging dies, can be easily predicted. Therefore, for the compensation of die-elasticity, die shape should be changed through size modification. Most previous studies about blade forging simulation have been performed in 2D finite element, and with this simplification, errors due to die-elasticity characteristics have been compensated. Clearly the forging of a blade is a 3D process, and the blade does not have the same cross section along its length and hence the forging of the blade could not be considered a plane-strain process. In this research, commercially available software, DEFORM3D, was used for the purpose of finite element method simulation for the turbine blade forging. Using three dimensional finite element simulation, elasticity deflection of dies in different sections of a blade were obtained, and then the upper and lower die profiles of those sections in various directions were modified and errors due to die-elasticity compensated. Through these modifications, errors due to die elasticity were reduced to values less than those obtained in previous works. A series of multiple simulations for having an optimized compensation of the die-elasticity were carried out. Errors due to elastic deflection were also compensated by modifying the die position, based on elastic deflection magnitude.Keywords: Gas turbine blade forging, die elastic deflection, three dimensional FEM simulations
-
Pages 61-68Mass conservation equations of three-phase incompressible, immiscible flow in porous media, called Buckley-Leverett (BL) equations, are of great importance in the study of flow in hydrocarbon reservoirs. For two-phase flow, these equations reduce to a scalar BL equation, in which the flux function is non-convex, as its convexity changes with respect to the conserved variable. Similarly, in three-phase flow, the non-convexity of flux functions leads to a linearly degenerate system of equations. In both non-convex and degenerate problems, it is necessary that the so-called Oleinik entropy condition is satisfied, in order to have physically admissible answers. Moreover, when gravity is included in the BL equations, the non-convexity of the flux function increases and some of the wave speeds of the system of equations become zero. Such a situation refers to a sonic point in gas dynamics. Many numerical methods have difficulty dealing with these points, especially when sonic points occur on an expansion wave. In such cases, some numerical methods produce non-physical expansion shocks instead of expansion fans. The latter case occurs when there is a counter current multi-phase flow in porous media. In order to avoid this kind of problem, entropy correction procedures should be used. These procedures usually generate additional numerical diffusion at sonic points to disperse expansion shocks. In this work, a relatively new numerical method, called the dominant wave method, is used to simulate immiscible flow in porous media. The dominant wave method is represented in the form of central methods. This method is locally conservative and uses a finite volume numerical flux approximation. The dominant wave method eliminates the necessity of characteristic decomposition by detecting the dominant wave speed of the system. On the other hand, this method generates less numerical diffusion in comparison with conventional central methods. Here, slope limiters are used to obtain a higher order of accuracy. Fromm's and minmod limiters are used in this paper, and a comparison between these two limiters is provided. To avoid expansion shocks at sonic points, Harten's entropy correction is used. In this work, gravitational effects are studied in a homogenous, one-dimensional reservoir using the Buckley-Leverett model. Three different models for relative permeability functions are considered in this study, including Corey, Extended Corey and Stone type models.Keywords: Gravitational effects, Buckley, Leverett equations, dominant wave method, entropy correction
-
Pages 69-76A series of wind tunnel tests were conducted to examine the flow field over swept wings at subsonic speed under various conditions. For this purpose, three models with sweep angles of 23, 33 and 40 degrees were used. To Increase the Reynolds number, all models were employed as semi span wing types and the effects eliminated by installing the models on a flat plate at their roots. The flat plate was installed in the middle of the test section at a distance from the wall. The tests were performed at various Reynolds numbers and at various angles of attack on all semi-span wings, which had the same aspect ratio but different sweep angles. The results show that increasing the sweep angle causes a decrease in the velocity over the wing surface. The pressure distribution on the upper surface of the wing shows some differences when compared to the 2-D case, due to the existence of the cross flow and wing tip vortices. As the wing sweep angle was varied, the differences varied too. The results indicate that the wing tip vortices amplify the 3D effects of the flow. The pressure reduction in the vicinity of the wing tip deteriorates the performance at that section, compared to other sections of the wing. This phenomenon would eventually cause stall to occur at this section (wing tip) sooner than at other sections. The maximum pressure drop near the wing tip remains unchanged with changing the sweep angle for zero angle of attack only. But, for other angles of attack, the situation is different. For the wing with a sweep angle of 40 degrees, the surface pressure of all sections differs significantly and the maximum pressure is much lower too. However, near the wing tip, the differences are not significant. This phenomenon is related to the effect of wing tip vortices that are more dominant at higher sweep angles. The test results in all three models show that the pressure coefficient decreases with increasing Reynolds, and this phenomenon is most prominent for the highest sweep angle of 40 degrees.Keywords: keywords{Swept wing, sweep angle, flow field, pressure distribution
-
Pages 77-85Due to recent advances in microfabrication techniques, it is possible to produce microchannels with positive, negative, or even neutral surface charges. Investigations of electroosmotic flow indicate that such a combination of charge patterns on the microchannel walls results in complex flow fields with circulation zones that are highly desirable for fluid mixing requirements, as in lab-on-a-chip devices. However, the mixing efficiency of such flows is not well-examined. The objective of this article is finding a proper mixing indicator among those which have been frequently used for evaluation of mixing, to examine the effect of wall heterogeneities on mixing efficiency. In general, for the numerical study of electro-osmotic flows, it is necessary to solve the electric potential field and the equations for ion transport, together with Navier-Stokes equations, which are numerically expensive and time-consuming, especially when heterogeneous channels are modeled. However, under certain conditions, the Helmholtz-Smoluchowski model can be used, which efficiently reduces numerical expense. In this model, the flow field is obtained by solving the Navier-Stokes equations with slip boundary conditions at walls, which comes from the Helmholtz-Smoluchowski slip velocity. Upon solving for the flow field, the concentration field can be obtained and the mixing efficiency can be evaluated. In most studies, the standard deviation of concentrations at any cross-section is considered the mixing indicator. Our investigations of the complex flow patterns with circulation zones show that this indicator is not well-behaved near the vortex area, such that some fluctuations in the mixing values occur. In order to have a more physically relevant mixing indicator, a weighted standard deviation is used to quantify the mixing performance, which is more consistent with vortex flows. Based on this indicator, some essential concepts for mixing, such as mixing performance, and channel efficiency, are introduced and the effects of step-wise heterogeneities are studied. Different zeta potential distributions associated with the single patch were examined, and it was found that the mixing performance of a non-continuous potential distribution is relatively similar to that of continuous counterparts, despite the slight differences in the flow field. The results for a single patch indicate that the channel efficiency of a channel with a single patch linearly increases with the size of the patch, such that for a channel with a patch size of channel height, the mixing efficiency increases 4%, compared to its homogeneous channel. Also, the mixing performance is enhanced as the single patch is located closer to the channel inlet. The findings of the present study can be employed in the design of optimized micro-mixers.Keywords: Heterogeneous microchannel, electro, osmo, tic flow, mixing performance
-
Pages 87-92The main objective of this paper is to determine the aerodynamic coefficients of a steady supersonic wing using the Boundary Element Method (BEM) based on the solution of potential flow. In the field of computational aerodynamics and fluid dynamics, some numerical approaches, such as the boundary element method, based on finding unknown parameters only on the boundary of a body with lower computational cost than high-fidelity formulations, can be a good replacement for Computational Fluid Dynamics (CFD) methods. Hence, application of BEM based on potential flow has recently been implemented as an efficient and popular tools in many studies. An advantage of the BEM is that the method reduces problem dimensionality by one. In this study, the governing partial differential equation of a supersonic potential flow field is firstly derived and linearized. The motion of the wing is assumed to consist of small perturbations, with respect to the constant-speed motion. Using the boundary element method in conjunction with first order shape functions, one obtains the system of discrete integral equations relating the potential to its normal derivatives on the surface of the body. However, there are some numerical problems in the evaluation of integral coefficients, due to singularity, which is caused by the intersection of the mach forecone and the edge of some elements in the computational domain. Of course, this problem can be eliminated by using some numerical techniques, such as the higher point Gaussin quadrature rule of integration and/or Teles transformation based on the order of singularity. Finally, a steady supersonic flow analysis is performed for two different wings with identified geometries (biconvex and planeform wings). The perturbed potential distribution for the biconvex wing are firstly computed and compared with available results in the literature. Also, its pressure coefficient distribution, specifically in the chordwise centerline, is evaluated. In addition, in the planeform wing, variation of the lift coefficient slope, accounting for three-dimensional effects, is studied. Also, the mach cone formed at both edges of the wing is visible and corresponds to the well-known mach cone formula for estimating mach angle. Comparison between the numerical results and those obtained using the exact solution indicates the good agreement and accuracy of the present method, and its flexibility to steady supersonic flows analysis around three-dimensional wings or complex configurations.Keywords: Aerodynamic, wing, Boundary Element Me, thod, potential flow, supersonic
-
Pages 93-101Due to the complexity of governing equations on robotic manipulator motion, various methods of control have been proposed. Information about the physical parameters, which are assumed to be default, is the key point of having a successful model-based control implementation. In other words, accurate identification of model parameters is always required in the case of control tasks. Unfortunately most physical characteristics are not exactly known, due to uncertainties in manufacturing, such as tolerances, clearances; or inherent uncertainties, like friction, and so on. Thus, the main objective of this study focuses on finding the kinetic parameters of a manipulator in a system identification framework, in which the sought parameters are estimated via the least squares method. Here, system identification is indeed an instance of the robot calibration problem, by which the best estimation for physical parameters of the system is extracted for use in controlling the model. To this end, first, a four-DOF planar manipulator is considered as an instant case study, and analyzed to achieve a straightforward dynamic model. If the dynamic model of the system is available, the usual way for identifying the parameters is minimizing the error between the outputs measured from the robotic manipulator and those calculated from the mathematical model, which should be obtained for an identical set of inputs. One of the most traditional and still most common criteria for minimizing errors is the least squares method, which is utilizable for any system that is linear in terms of unknown parameters. On the other hand, in the literature, the routine dynamical model for controlling tasks is derived based on choosing relative angles as the generalized coordinates; then, Coriolis acceleration effects, including the product of joint velocities, arise in the equations of motion. Such a modeling will possess extreme parameters in its formulation, and so, identifying kinetic characteristics, will be involved. Thus, in this paper, by choosing the absolute angles of the links as the generalized coordinates, the Coriolis acceleration terms could be hidden in the structure of the equations of motion, and, in this way, our formulation will include minimum unknown parameter numbers. The input to this system is the torque applied to the manipulator joints, which makes it an MIMO system. By considering the cost function as the error between the estimated parameters and the true ones, the trace of the error covariance matrix is minimized, and, in this way, estimation is achieved. The simulation results assert the validity of this identification.
-
Pages 103-111
C o m b u s t i o n h e a t e r s c a n r e l e a s e a l a r g e a m o u n t o f h a r m f u l g a s e s, s u c h a s C O a n d N O x i n t o t h e a t m o s p h e r e o f a r o o m. B u t, b y p r o p e r d e s i g n a n d a r r a n g m e n t o f t h e c o m b u s t i o n e q u i p m e n t i n t h e r o o m, t h e r m a l b u o y a n c y c a n b e u s e d t o d r i v e o u t t h e p o l l u t a n t s a n d i m p r o v e i n d o o r a i r q u a l i t y. T h e m a i n a i m o f t h i s s t u d y i s t o p r o p o s e a p r o p e r a r c h i t e c t u r a l d e s i g n, a n d p r o v i d e s o m e p r a c t i c a l d e s i g n r e c o m m e n d a t i o n s, b y e v a l u a t i n g t h e e f f e c t o f d i f f e r e n t p a r a m e t e r s o n t h e i n d o o r a i r q u a l i t y o f b u i l d i n g s w i t h c o m b u s t i o n h e a t i n g e q u i p m e n t. T o a c h i e v e t h e s e g o a l s, v a r i o u s a r r a n g e m e n t s f o r l o c a t i o n s o f h e a t e r s a n d o u t l e t v e n t s h a v e b e e n c o n s i d e r e d, a n d i n d o o r a i r q u a l i t y h a s b e e n e v a l u a t e d u n d e r t h e s e c o n d i t i o n s b y a n a l y z i n g a i r f l o w a n d p o l l u t a n t d i s t r i b u t i o n d u e t o v e n t e d a n d u n v e n t e d h e a t e r s i n t h e r o o m. I n o r d e r t o s o l v e t h e p r o b l e m, a n u m e r i c a l s i m u l a t i o n o f t h e f l o w a n d p o l l u t a n t f i e l d s i n a 3D r o o m w a s p e r f o r m e d u s i n g A i r p a k s o f t w a r e. I n v e s t i g a t i o n s i l l u s t r a t e t h a t t h e r m a l b u o y a n c y h a s d i f f e r e n t e f f e c t s i n v a r i o u s a r r a n g e m e n t s o f h e a t e r a n d o u t l e t v e n t. C h a n g i n g t h e l o c a t i o n o f t h e o u t l e t v e n t o n t h e s i d e w a l l s c a n n o t s i g n i f i c a n t l y a f f e c t i n d o o r a i r q u a l i t y, b e c a u s e n a t u r a l c o n v e c t i o n c u r r e n t s b y t h e h e a t e r a r e t h e d o m i n a n t m e c h a n i s m i n f l o w f i e l d d i s t r i b u t i o n. H o w e v e r, t h e r e s u l t s s h o w t h a t t h e c e n t e r o f a c o l d e x t e r n a l w a l l i s t h e b e s t l o c a t i o n f o r i n s t a l l i n g t h e o u t l e t v e n t, s o t h a t b o t h c a t e g o r i e s o f a i r q u a l i t y a n d e n e r g y c o m b u s t i o n a r e c o n s i d e r e d. M o r e o v e r, t h e r e s u l t s o f t h e p r e s e n t s t u d y d e m o n s t r a t e t h a t t h e l o c a t i o n o f t h e h e a t e r i n a r o o m c a n s i g n i f i c a n t l y a f f e c t t h e d i s t r i b u t i o n a n d a l s o m e a n c o n c e n t r a t i o n v a l u e o f t h e p o l l u t a n t s. B y i n s t a l l i n g t h e h e a t e r n e a r t h e e x t e r n a l w a l l, o p p o s i t e t h e o u t l e t v e n t a n d i n t h e c o r n e r s o f t h e r o o m, 6% r e d u c t i o n i n C O2 c o n c e n t r a t i o n i n t h e b r e a t h i n g z o n e i n t h e b a s e m o d e l i s o b s e r v e d. T h e p o s i t i o n o f t h e h e a t e r t o t h e c o l d e x t e r n a l w a l l, o u t l e t v e n t, d o o r a n d w i n d o w a r e i m p o r t a n t f a c t o r s t h a t c a n a f f e c t t h e c o n c e n t r a t i o n o f t h e p o l l u t a n t s. G e n e r a l l y, t h e l o w e s t c o n c e n t r a t i o n o f p o l l u t a n t s i n t h e r o o m i s a c h i e v e d i n c a s e s w h e r e t h e h e a t e r w a s l o c a t e d n e a r t h e l o n g i t u d i n a l e x t e r n a l w a l l a n d n e a r t h e c o r n e r o f t h e r o o m, a n d t h e o u t l e t v e n t w a s l o c a t e d i n a n o t h e r c o r n e r o f t h a t w a l l. F u r t h e r m o r e, t h e h i g h e s t c o n c e n t r a t i o n o f p o l l u t a n t s i n t h e r o o m w a s a c h i e v e d i n c a s e s w h e r e t h e h e a t e r w a s l o c a t e d n e a r t h e c o r n e r o f t h e l o n g i t u d i n a l i n t e r n a l w a l l, a n d t h e o u t l e t v e n t w a s l o c a t e d i n t h e c o r n e r o f t h e t r a n s v e r s e i n t e r n a l w a l l. C o m p a r i n g i n d o o r a i r q u a l i t y i n t h e m e n t i o n e d c o n d i t i o n s s h o w e d t h a t u n d e r t h e b e s t c o n d i t i o n s, t h e c o n c e n t r a t i o n o f C O2 i n t h e b r e a t h i n g z o n e i s a b o u t 9% l o w e r t h a n t h e w o r s t c o n d i t i o n s.
-
Pages 113-124One of the most commonly used techniques in obtaining three dimensional measurements is the triangulation technique using structured light, which is based on a simple mathematical algorithm. The current paper deals with 3D scanning and model reconstruction of mechanical parts being produced in a continuous manner. The objective of the current research is to firstly scan and reconstruct the part model and secondly to improve the accuracy of the model reconstruction. A simple test rig is designed and built to scan the object. The test rig is made of optoelectrical and mechnical parts. The optoelectrical parts include a digital CCD camera and a laser light projector. A conveyor belt with controllable constant speed is used to generate the movement required when scanning. The camera is located above the conveyor belt in a fixed unknown position. The laser light projector is installed on a five degrees of freedom robot arm facilitating the adjustment of the projector to the desired position and orientation. The final accuracy of the reconstructed model is considerbly affected by several factors, namely; the accuracy of the intrinsic and extrinsic parameters of the camera, the light projection angle, and the uniformity of light distribution over the entire scene, both in terms of location and time. Camera calibration is performed to identify camera parameters. Twelve parameters were identified in total. The parameters found were verified through measuring previously known 2D objects. Novel simple mirror arrangements are employed to accurately measure the light projection angle. A triangulation technique, along with a laser sheet of light with Gaussian distribution, is used to generate depth information. The intersection of the sheet of light and the object produces a 3D curve. Deviations of the curve from a striaght line represent the depth information, which can be extracted and translated into height variations using the triangulation technique. The CCD camera is used to capture both the texture and the depth carrying signal. The captured data are transferred, pre-processed, processed and converted to CAD data format. Curve approximation techniques are used to smooth the curve. Consequently, the CAD model is generated in a CAD environment. The model is then exported to an analysis environment for further analysis. To improve the accuracy of measurement, camera calibration, off-axis light source compensation, and light source parameter identification are performed. Algorithms to calibrate the camera, to scan the object, to smooth out the data, and to construct the 3D model were developed. To test the functionality of the system, and to evaluate the accuracy of the developed system, a controlled experiment is designed and performed. The system works well and a maximum error of less than 0.01% of full measurement scale is obtained.Keywords: 3D scanning, triangulation, 3D model, model reconstruction
-
Pages 125-134In this work, a computational algorithm is given for solving the transport equation arising from flow in two-dimensional porous media, using the streamline method. The streamline method has four vital steps. First, the two-dimensional transport equation is decoupled into multiple 1-D equations. By using a new local variable, called the time-of-flight, these equations are solved in a way that is independent from the spatial discretization. In the second step, the streamlines are traced, based on the velocity field in the computational domain, to obtain the time-of-flight grid (TOF gird) along each streamline. Cell-by-cell streamline tracing is performed by utilizing the semi-analytical Pollock algorithm from the injector (sink) cell faces to the producer (source) cell faces or vice versa. In the third step, once all streamlines are traced, and their 1-D TOF grids are constructed, the multiple 1-Dequations are solved, by either analytical or numerical methods. In the fourth step, the solutions along the TOF grid are mapped into the pressure grid (primary structured mesh) in order to construct the solution of the 2-D transport equation. This step is conducted using a volumetric averaging of the streamlines passing each cell. The streamline method requires relatively lower computational time than other conventional methods, such as finite difference and finite volume methods, as it solves only several 1-D equations, and is without restriction on the time step size. In addition, this method requires less memory than other standard methods, as it needs to save and solve one of the multiple 1-D equations each time. Here, this method is employed for solving a number of 2-D test cases, including both single-phase and two-phase flow. The test cases demonstrates that the streamline method has good accuracy with a high speed-up factor. Numerical results show that using this method to simulate flow within a 2-D problem with 90000 grid cells leads to a speed-up factor of about 26. It is also shown that the speed-up factor of the streamline method will be enhanced by increasing the heterogeneity and number of grid cells.Keywords: Streamline method, two phase, time of flight, Pollock algoritm
©2000-2025 magiran.com All Rights Reserved.