برآورد خطا و بهبود تنش در تحلیل غیرخطی مصالح به روش ایزوژئومتریک

نوع مقاله : پژوهشی اصیل (کامل)

نویسندگان
1 عضو هیات علمی دانشگاه آزاد اسلامی واحد شاهرود
2 دانشجوی دکتری سازه دانشگاه آزاد اسلامی واحد شاهرود
چکیده
همگام با رشد علوم و فناوری، مسائل مهندسی نیز روز به روز پیچیده تر می‌شوند. با پیچیده تر شدن مسائل و لزوم حل سریعتر و دقیق تر آن­ها، روش­های تحلیلی گذشته دیگر جوابگوی نیازهای روز افزون جوامع نیستند. با چنین نگرشی، محققان همواره سعی کردند در کنار توسعه مبانی علوم، روش­های عددی را نیز توسعه بخشند. دراین مسیر، روشهای متعددی توسط محققین ابداع گشته است. هر کدام از این روش­ها ،کاربرد خاص خود را دارند و هنوز هم محققان درصدد رشد و توسعه این روش­ها و ابداع روش­های جدید هستند. از مهمترین این­ها می‌توان به روش ایزوژئومتریک غیر خطی که بر اساس بی­اسپیلاین‌های نسبی غیر یکنواخت به وجود آمده است، اشاره کرد. در روش ایزوژئومتریک غیر خطی ضمن استفاده از خواص توابع پایه اسپیلاین و نربز در تعریف دقیق منحنی‌ها و سطوح، از آن­ها جهت درونیابی و تقریب سازی هم استفاده می­شود. استفاده از همه ظرفیت سازه در تحمل بار، باعث رفتار غیرخطی سازه می­شود که ناشی از عملکرد نامناسب هندسه سازه، ضعف مصالح سازه و نارسائی ناشی از ترکیب دو حالت قبل است. و در این پژوهش غیرخطی شدن ناشی از ضعف مصالح مدنظر قرار داده شده است. همچنین در حل معادلات تعادلی غیر خطی از یک روند افزایشی و تکراری بار استفاده شده و این افزایش تا وارد شدن کل بار­های تعریف شده برای هر مسئله انجام شده است. در هر افزایش، روند تکراری اتخاذ شده تا جایی انجام می­گیرد که همگرایی یا ماکزیمم تعداد تکرارها بدست آید. خطا بخش جدانشدنی تحلیل­های عددی به شمار می‌رود و همواره باعث نگرانی محققین در قابلیت اعتماد نتایج بوده است. بدیهی است که کلیه روش­های عددی روش­های تقریبی می­باشند. عمده ترین منبع خطا در روش­های عددی مربوط به خطای گسسته سازی محیط پیوسته بوده و ناشی از تقریب زدن میدان جابجایی به وسیله توابع شکل است. این گروه ازخطاها نیز با ریزتر کردن شبکه المان بندی و بالا بردن درجه توابع شکل مورد استفاده، کاهش می‌یابد. لذا در این پژوهش نسبت به برآورد خطا براساس روش بازیافت تنش بر مبنای نقاطی که مرتبه همگرایی گرادیان یک تابع، یک مرتبه از مقداری که از تقریب تابع شکل مربوط به حل تقریبی انتظار می رود، بالاتر است (نقاط فوق همگرا) پرداخته شده است. بدین صورت که با در نظر گرفتن اختلاف بین سطح تنش بازیافتی و سطح تنش بدست آمده از تحلیل ایزوژئومتریک غیرخطی برای هر المان، به صورت تقریبی به یک معیاری جهت تعیین میزان خطای موجود در آن المان تعیین شده است. کلیه رابطه­ سازی­های پژوهش و خطی‌سازی معادلات با استفاده از یک الگوریتم عددی با کمک برنامه نویسی در محیط نرم افزار فرترن انجام شده و نتایج تحلیل جهت صحت سنجی با حل کلاسیک آن مقایسه شده است. نتایج تشابه عددی و توزیعی قابل قبولی را نشان داده است؛ لذا می توان بیان کرد تحلیل صورت گرفته توسط برنامه از کارایی مناسبی جهت تحلیل غیرخطی مسائل برخوردار است. همچنین روش تخمین کننده خطای به کار گرفته شده را می­توان راه حلی ساده و مهندسی جهت برآورد خطا و بهبود میدان تنش بدست آمده از تحلیل الاستوپلاستیک مسائل به روش ایزوژئومتریک نام برد.

کلیدواژه‌ها

موضوعات


عنوان مقاله English

Error Estimation and Stress Improvement in Nonlinear Analysis of Materials by Isogeometric Method

نویسندگان English

A. ganjali 1
A. Mirzakhani 1
A. Shahini 2
1 Faculty member of Islamic Azad University, Shahrood Branch
2 PhD student in Structural Islamic Azad University, Shahrood Branch
چکیده English

With the growth of science and technology, engineering issues are becoming more complex. As problems become more complex and need to be resolved more quickly and accurately, past analytical methods no longer meet the growing needs of societies. With such an attitude, researchers have always tried to develop numerical methods in addition to developing the basics of science. In this direction, several methods have been developed by researchers. Each of these methods has its own applications and still researchers are trying to grow and develop these methods and invent new methods. The most important of these are the nonlinear isogeometric method which is based on non-uniform rational B-Splines (NURBS). In the nonlinear isogeometric method, while using the properties of the basic functions of spline and NURBS in the exact definition of curves and surfaces, they are also used for interpolation and approximation. Using all the capacity of the structure in load bearing causes nonlinear behavior of the structure which is due to improper performance of the structure geometry, weakness of the structural materials and weakness due to the combination of the two previous states. In this study, nonlinearity due to material weakness has been considered. Also, in solving nonlinear equilibrium equations, an incremental and iterative process of load is used and this increase is done until the total loads defined for each problem are entered. In each increase, the iterative process is adopted until convergence or the maximum number of iterations is achieved. Obviously, all numerical methods are approximate methods. The main source of error in numerical methods is related to the discretization error of the continuous environment and is due to the approximation of the displacement field by the shape functions. This group of errors is also reduced by making the elemental mesh smaller and increasing the degree of shape functions used. Error is an integral part of numerical analysis and has always been a concern for researchers in the reliability of the results. Therefore, in this study, the error estimation based on the stress recovery method based on points where the order of gradient convergence of a function is one time higher than the value expected from the approximation of the shape function related to the approximate solution (superconvergent points) is discussed. Thus, by considering the difference between the recovered stress level and the stress level obtained from nonlinear isogeometric analysis for each element, a criterion has been determined approximately to determine the amount of error in that element. All research relationalizations and linearization of equations have been performed using a numerical algorithm with the help of programming in Fortran software environment and the results of the analysis for validation have been compared with its classical solution. The results show acceptable numerical and distributive similarity; Therefore, it can be said that the analysis performed by the program has good performance for nonlinear analysis of problems. Also, the error estimation method used can be called a simple and engineering solution to estimate the error and improve the stress field obtained from elastoplastic analysis of problems by isogeometric method.

کلیدواژه‌ها English

Nonlinear Isogeometric Method
Error Estimation
Stress Recovery Method
Superconvergent Points
Improve The Stress Field
1- Kagan, P., Fischer, A., & Bar‐Yoseph, P. Z. J. I. J. f. N. M. i. E. (1998). New B‐spline finite element approach for geometrical design and mechanical analysis. 41(3), 435-458.
2- Höllig, K., Reif, U., & Wipper, J. J. S. J. o. N. A. (2001). Weighted extended B-spline approximation of Dirichlet problems. 39(2), 442-462.
3- Kagan, P., Fischer, A., & Bar‐Yoseph, P. Z. J. I. J. f. N. M. i. E. (2003). Mechanically based models: Adaptive refinement for B‐spline finite element. 57(8), 1145-1175.
4- Hughes, T. J., Cottrell, J. A., Bazilevs, Y. J. C. m. i. a. m., & engineering. (2005). Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. 194(39-41), 4135-4195.analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. 194(39-41), 4135-4195.
5- Turner, M. J., Clough, R. W., Martin, H. C., & Topp, L. (1956). Stiffness and deflection analysis of complex structures. journal of the Aeronautical Sciences, 23(9), 805-823.
6- Argyris, J. H. J. N. Y., MACMILLAN CO., OXFORD, PERGAMON PRESS, LTD., . 187 P. (1964). Recent advances in matrix methods of structural analysis(Matrix theory of structures for small and large deflections, using high speed digital computers).
7- Argyris, J. (1965). Continua and Discontinua, opening address to the 1-st Conf. Matrix Methods in Structural Mechanics. In: Wright--Patterson AFB, Dayton, Ohio.
8- Oden, J. T. J. J. o. t. S. D. (1967). Numerical Formulations of Nonlinear Elasticity Problems. 93(3), 235-356.
9- Mallett, R. H., & Marcal, P. V. J. J. o. t. s. d. (1968). Finite element analysis of nonlinear structures. 94(9), 2081-2106.
10- Oden, J. (1969). Finite element applications in nonlinear structural analysis. Paper presented at the Proceedings of the ASCE Symposium on Application of Finite Element Methods in Civil Engineering.
11- Haisler, W. E., Stricklin, J. A., & Stebbins, F. J. J. A. J. (1972). Development and evaluation of solution procedures for geometrically nonlinear structural analysis. 10(3), 264-272.
12- Zinckiewicz, O. (1971). The finite element in engeneering science. In: Mc Graw-Hill, London.
13- Brebbia, C., & Connor, J. J. J. o. t. E. M. D. (1969). Geometrically nonlinear finite-element analysis. 95(2), 463-486.
14- Hassani, B., Tavakkoli, S. M., & Ardiani, M. (2015). Solution of Nonlinear Incompressible Hyperelastic Problems by Isogeometric Analysis Method %J Journal of Solid and Fluid Mechanics. 5(2), 29-41. doi:10.22044/jsfm.2015.429
15- Richardson, L. F. (1911). IX. The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam. Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character, 210(459-470), 307-357.
16-Babuška, I., & Rheinboldt, W. C. J. I. J. f. N. M. i. E. (1978). A‐posteriori error estimates for the finite element method. 12(10), 1597-1615.
17-Babuška, I., Rheinboldt, W. J. C. M. i. A. M., & Engineering. (1979). Adaptive approaches and reliability estimations in finite element analysis. 17, 519-540.
18-Babuška, I., Strouboulis, T., Upadhyay, C. J. C. M. i. A. M., & Engineering. (1994). A model study of the quality of a posteriori error estimators for linear elliptic problems. Error estimation in the interior of patchwise uniform grids of triangles. 114(3-4), 307-378.
19-Babuška, I., Strouboulis, T., Upadhyay, C., Gangaraj, S., & Copps, K. J. I. j. f. n. m. i. e. (1994). Validation of a posteriori error estimators by numerical approach. 37(7), 1073-1123.
20-Zienkiewicz, O. C., & Zhu, J. Z. J. I. J. f. N. M. i. E. (1992). The superconvergent patch recovery and a posteriori error estimates. Part 1: The recovery technique. 33(7), 1331-1364.
21- Hassani, B., Ganjali, A., & Hojatpanh Montazary, A. (2012). Analysis and Shape Optimization of Axsymmetric Structures by Isogeometric Analysis Method. Journal of Solid and Fluid Mechanics, 1(1), 1-13.
22- Hassani, B., Ganjali, A., & Tavakkoli, M. (2012). An isogeometrical approach to error estimation and stress recovery. European Journal of Mechanics-A/Solids, 31(1), 101-109.
23-Piegl, L., & Tiller, W. (2012). The NURBS book: Springer Science & Business Media.
24- Owen, D. R. J. (1980). Finite elements in plasticity, theory and practice.
25- Sheng, D., Sloan, S. W., & Abbo, A. J. (2002). An automatic Newton–Raphson scheme. The International Journal Geomechanics, 2(4), 471-502.
26- Nayak, G. C., & Zienkiewicz, O. C. (1972). Convenient form of stress invariants for plasticity. Journal of the Structural Division, 98(4), 949-954.
27-Boroomand, B., & Zienkiewicz, O. (1999). Recovery procedures in error estimation and adaptivity. Part II: Adaptivity in nonlinear problems of elasto-plasticity behaviour. Computer methods in applied mechanics and engineering, 176(1-4), 127-146.
28- Zienkiewicz, O. C., Taylor, R. L., & Zhu, J. Z. (2005). The finite element method: its basis and fundamentals: Elsevier.
29- Lubliner, J. 1990, Plasticity Theory. Macmillan Publishing Company, New York, pp.239-244.
30- Chakrabarty, J. „1987, Theory of Plasticity. In: McGraw Hill, New York, NY, pp.164-171.
31- de Souza Neto, E. A., Peric, D., & Owen, D. R. (2011). Computational methods for plasticity: theory and applications: John Wiley & Sons, pp. 387-388.
32- Hill, R. (1950). The mathematical theory of plasticity, Clarendon. Oxford, 613, 614, pp.98-110.
33- Lubliner, J. 1990, Plasticity Theory. Macmillan Publishing Company, New York, pp.216-228.
34- Chakrabarty, J. „1987, Theory of Plasticity. In: McGraw Hill, New York, NY, pp.323-333.
35- de Souza Neto, E. A., Peric, D., & Owen, D. R. (2011). Computational methods for plasticity: theory and applications: John Wiley & Sons, pp. 244- 247.