تحلیل مسائل دوبعدی با روش اجزای محدود درهم بر پایه پاسخ تحلیلی معادله دیفرانسیل

در این مقاله یک جزء مرتبه بالای هشت گرهی برپایه‌ی پاسخ تحلیلی معادله دیفرانسل حاکم، برای تحلیل سازه‌های دوبعدی پیشنهاد می‌گردد. رابطه‌سازی جزء پیشنهادی بر پایه‌ی تابعی درهم هلینگر-رایزنر و پاسخ تحلیلی معادله سازگاری حاکم بر مسئله‌های دوبعدی انجام می‌پذیرد. شایان ذکر است جهت رابطه سازی اجزای محدود با تابعی درهم هلینگر-رایزنر، نیاز به دو میدان مستقل تنش و جابجایی در درون جزء می‌باشد. برای این منظور ابتدا، با حل تحلیلی معادله سازگاری، تابع‌های تنش آیری در دسترس قرار می‌گیرد. با بهره‌جویی از این تابع‌های تنش، میدان تنش درون جزء به دست می‌آید. هم‌چنین،  میدان جابجایی درجه دوم هم‌عامل هشت گرهی برای جابجایی درون جزء به کار می‌رود. با به کاربردن تابعی درهم هلینگر-رایزنر و ایستا کردن آن نسبت به میدان‌های مستقل تنش و جابجایی، ماتریس سختی و بردار نیروهای گرهی جزء در دسترس قرار می‌گیرند. در پایان با آزمون‌های عددی گوناگون، دقت و کارایی جزء پیشنهادی مورد ارزیابی قرار می‌گیرد. این آزمون‌ها دقت بسیار بالای جزء پیشنهادی را در تحلیل سازه‌های دوبعدی به اثبات می‌رسانند.

mixed finite element formulation for 2d problems analysis based on analytical solutions of deferential equation

in this paper, a high-order eight-node element based on the analytical response of the governing differential equation is proposed for the analysis of plane structures. the formulation of the proposed element is based on the hellinger-reisner mixed functional and the analytical response of the compatibility equation governing plane problems. it is worth noting that in order to formulate finite elements with the hellinger-reisner functional, two independent stress and displacement fields are required. for this purpose, airy stress functions are first made available by the analytical solution of the compatibility equation. by utilizing these stress functions, the stress field within the element is obtained. also, the quadratic displacement field of the isoparametric element is used for intra-element displacement. by applying the hellinger-reisner mixed functional and stationary of this functional relative to the independent stress and displacement fields, the stiffness matrix, and the element node force vector are made available. finally, with various numerical tests, the accuracy and efficiency of the proposed element are evaluated. these tests prove the high accuracy of the proposed element in the analysis of plane structures.