18312

۱۸۳۱۲

تحليل تنش دوبعدي قطعات ترك دار تحت خزش توسط يك كد بدون شبكه و اعتبارسنجي آزمايشگاهي آن

دكتري

طراحي كاربردي

دي ماه ۱۳۹۶

دكتر محمد حق پناهي

مكانيك

در سال‌هاي اخير روش‌هاي بدون شبكه به سرعت در حوزه مكانيك شكست گسترش يافته‌اند. اين روش‌ها به علت عدم استفاده از المان، زمان مدل‌سازي را در مسائل صنعتي به نسبت زيادي كاهش مي‌دهند. در كنار اين مساله، مزايايي همچون دقت بالا در نتايج تحليل تنش به خصوص در مناطق نزديك ترك كه داراي گراديان تنش بالايي مي‌باشد، عدم وجود ناپيوستگي در مقادير تنش و كرنش و سهولت در شبيه‌سازي رشد ترك، باعث شده است تا گرايش محققين به اين روش عددي در حل مسائل مكانيك شكست بيشتر شود. در اين رساله، يك كد بدون شبكه جهت تحليل تنش دوبعدي قطعات ترك‌دار تحت خزش توسعه داده شده است. كد نوشته شده با استفاده از يك دستگاه تست خزش مجهز به امكانات اندازه‌گيري رشد ترك خزشي، اعتبارسنجي شده است. جهت ارزيابي كد بدون شبكه توسعه داده شده، چند مثال شامل بررسي ضريب تمركز تنش در صفحه با سوراخ دايروي و بيضوي و همچنين بررسي ضريب شكل شيار در ديسك برزيلي با شيار V‌ شكل حل شده‌اند و نتايج حاصله با نتايج حاصل از تحليل اجزاء محدود مقايسه شده‌اند. سپس يك معيار جديد جهت شبيه‌سازي ناپيوستگي ترك ارائه شده است. در اين معيار، از تركيب دو روش رويت و مثلث‌هاي مجاور جهت تعيين دامنه پشتيبان استفاده شده است. مزيت معيار ارائه شده نسبت به معيارهاي متداول مورد استفاده در روش‌هاي بدون شبكه، اين است كه در اين روش، با استفاده از المان‌هاي اجزاء محدود و تركيب دو روش متداول ذكر شده، بدون نياز به استفاده از روابط رياضي جهت تصحيح پارامتر تابع وزني، دامنه پشتيبان در گره‌ها و نقاط انتگرال‌گيري اطراف ترك تصحيح شده است. در بخش بعد، توسط بكارگيري روش كرنش اوليه، تحليل خزش در روش بدون شبكه گلركين انجام شده است. سپس در مرحله بعد، تابع بي‌بعد h1 كه مهمترين پارامتر در محاسبه انتگرال C* است، براي قطعات ترك‌دار محاسبه شده است و با روابط كومار مقايسه شده است. خطاي محاسبات در اين بخش به طور ميانگين حدود 8 درصد بوده‌است. در گام بعدي، با بكارگيري روش كاهش مدول الاستيسيته، انتگرال C* در روش تنش مرجع تخمين زده شده است و نشان داده شده است كه اين تخمين يك تخمين محافظه كارانه است و بنابراين براي مقاصد طراحي مناسب است. در مرحله بعد، جهت اعتبار سنجي نتايج، تست‌هاي دما بالا براي سوپرآلياژ A286 در دماي 650c انجام شده است. تست‌هاي انجام شده شامل تست تعيين چقرمگي شكست و تعيين پارامترC* مي‌باشد. در بخش تعيين پارامترC* ، مقايسه‌ نتايج حاصل از تست با نتايج بخش تئوري صورت گرفته است و مشخص شده است رفتار قطعات تست شده به شرايط تنش صفحه‌اي نزديك‌تر است.

1396/10/20

1/10/2018 12:00:00 AM

In recent years, the meshless methods have been developed in fracture mechanics field. These methods, reduce the time which is needed for modeling of structures. because there is no need for element creation. Besides this matter, some advantages like high accuracy in computation of stresse and strains specially in stress raiser areas e.g. at the vicinity of the crack tip, getting smooth stresses and strains and crack growth simulation with no need to remeshing, caused that the researchers tend to this methods in solving fracture mechanics problems. In this thesis, a meshless code has been developed for two dimensional stress analysies of cracked bodies under creep condition. The results obtained by the developed code are validated by a creep testing machine, equipped with crack growth measurement facilities. For evaluation of the developed meshless code, some examples are solved including investigation of Stress Concentration Factor in a plate with circular hole, investigation of Stress Concentration Factor in a plate with elliptical hole and investigation of Notch Shape Factor in the V-notched Brazilian Disc. Then a new criterion is proposed for simulation of crack discontinuity. In fact, in this criterion, a combination of the visibility method and the surrounding triangles method is used for determination of support domain. The advantage of the proposed criterion in comparison with the common criteria in the meshless methods, is that by using the connectivity concept in the Finite Element Method and by combination of two mentioned methods, the support domain for the nodes and the integration points near the crack faces is modified without any special mathematical formulation for the weight function paremeter. After that, by using the Initial Strain Method, creep analysis is done in the Element Free Galerkin method. In the next step, the dimenless h1 function which is the most important parameter in the C* integral, is determined for the standard test specimens and is compared with the Kumar relationships. Then in the next step, C* parameter is estimated by implementing Elastic Modulus Reduction (EMR) approach in the reference stress method and it is shown that this is an over estimate procedure and therefore is usefull for design purposes. In the last step, high temperature tests are performed for A286 superalloy at 650c for validation of the results obtained in the previous sectons. The performed tests includes determination of the fracture toughness and the determination of the C* parameter. In the section that the C* parameter is determined, the obtained results are compared with the related relationships and it is concluded that the tested specimens are close to the plane stress condition.