بررسي مقايسه‌يي شرايط ماتريس سختي در روش‌هاي آناليز اجزاي محدود با المان ويژه، اجزاي محدود توسعه يافته و روش آناليز ايزوژئومتريك در مسئله‌ي شكست

A C‌O‌M‌P‌A‌R‌A‌T‌I‌V‌E I‌N‌V‌E‌S‌T‌I‌G‌A‌T‌I‌O‌N O‌N S‌T‌I‌F‌F‌N‌E‌S‌S M‌A‌T‌R‌I‌X C‌O‌N‌D‌I‌T‌I‌O‌N‌S F‌O‌R F‌E‌A, X‌F‌E‌A, A‌N‌D I‌G‌A I‌N F‌R‌A‌C‌T‌U‌R‌E P‌R‌O‌B‌L‌E‌M

ماتريس سختي يكي از ويژگي هاي ذاتي هر روش عددي است. شرايط اين ماتريس تاثير به سزايي در حجم محاسبات و پايداري روش حل مي گذارد. در مطالعه ي حاضر، به لحاظ نظري مقايسه يي بر اساس هندسه و ويژگي هاي رياضي ماتريس سختي حاصل از روش هاي عددي اجزاي محدود استاندارد، اجزاي محدود توسعه يافته و روش آناليز ايزوژيومتريك، در حل مسيله ي ورق حاوي ترك صورت گرفته است. به منظور مدل سازي ترك در روش اجزاي محدود استاندارد، المان ويژه تكين هشت گرهي اجزاي محدود به كار گرفته شده است. از طرفي در روش اجزاي محدود توسعه يافته از اصول غني سازي توابع درون ياب اجزاي محدود و كاربرد روش رياضي پيكره بندي واحد اجزاي محدود استفاده مي شود. به علاوه در روش آناليز ايزوژيومتريك از تكرار نقاط كنترلي در يك موقعيت به منظور ايجاد ضمني ناپيوستگي و ايجاد تكينگي در ماتريس سختي، بهره گرفته شده است و از توابع پايه مرتبه 3 در تقريب هندسه و حل استفاده مي شود.

Fracture mechanics is a vast and growing open research field of sciences which concerned with the study of crack propagation. Stiffness matrix is an inherent feature of a numerical method. Its condition has a great influence on numerical calculation and the stability of the solution. Since the application of numerical methods such as the standard finite element method, the extended finite element method and the isogeometric analysis approach in the problems of fracture mechanics has been approved, in this contribution, a theoretical comparison between stiffness matrices is conducted. For this purpose, three computer codes are prepared to make a comparison on the stiffness matrix geometry and mathematical properties such as matrix sparsity, stiffness index, bandwidth, number of independent rows and columns, zero and non-zero elements, symmetric/ nonsymmetric and hermitian/ nonhermitian. The 8-node singular element is used in the framework of the finite element method. and in the case of extended finite element, the principles of enrichment of the interpolation functions of finite element and the application of the partition of unity method is considered. Also in the isogeometric analysis approach repetition of two different control points between two patches can create a discontinuity and also demonstrates a singularity in the stiffness matrix. In addition, the NURBS of order 3 are utilized as the basis functions to approximate the geometry and the solution. By comparing the stress distribution, the accuracy of the calculations and the smoothness of the results are investigated. It is found that, stiffness matrix obtained from the isogeometric analysis method is non-diagonal in the fracture problems. Extended finite element stiffness matrix in comparison with the other methods possess a better numerical conditions.