Дослідження еволюції напружено-деформованого стану і визначення розрахункового ресурсу масивних елементів вісесиметричних конструкцій із використанням універсального скінченного елементу

Abstract

Математичне моделювання процесів континуального руйнування в умовах довготривалої повзучості на основі МСЕ являє собою досить складну задачу, ефективність розв'язання якої залежить від повноти бібліотеки скінченно-елементної бази і алгоритмів розв'язання систем нелінійних рівнянь, а також від організації програмного забезпечення. За основу процедури отримання скінченно-елементних розв'язувальних співвідношень покладені основні принципи моментної схеми скінченних елементів (МССЕ). На відміну від загальноприйнятих підходів використання МССЕ передбачає, крім завдання закону розподілення переміщень, незалежне подання деформацій у вигляді рядів Маклорена. В основу алгоритму розв'язання системи нелінійних рівнянь прийнято метод інтегрування по параметру навантаження, причому для отримання достовірних результатів передбачається його послідовне зменшення. Точність розв'язання системи нелінійних рівнянь на кожному кроці по параметру навантаження визначається порівнянням величини суми квадратів вузлових реакцій і суми квадратів вузлових значень зовнішніх навантажень. В якості вихідних співвідношень прийняті рівняння термов'язкопружнопластичності з урахуванням пошкодженості матеріалів. В основу побудови скінченно елементної моделі об'єктів, покладено принцип використання квазірегулярної фрагментації, яка передбачає побудову загальної нерегулярної скінченно-елементної моделі на основі регулярних скінченно-елементних фрагментів. Це дозволяє суттєво оптимізувати загальну кількість невідомих. При розрахунку нових об'єктів збіжність результатів обґрунтовується шляхом послідовного збільшення параметрів сіткової моделі та зменшення величини кроку інтегрування за рахунок збільшення їх кількості в межах заданого інтервалу навантажень, а також збільшення точності розв'язку системи нелінійних рівнянь на кожному кроці. В даній роботі для обґрунтування достовірності результатів, отриманих на основі універсального скінченного елемента, проведено порівняння параметрів напружено-деформованого стану і параметру пошкодженості ω з даними, обчисленими на основі скінченних елементів загального типу з чисельним інтегруванням і елементів з інтегруванням в замкненому вигляді.Mathematical modeling of the continual fracture processes under long creep conditions on the basis of FEM is a rather complex problem. The efficiency of solving of this problem depends on the completeness of the finite element library and the algorithms for solving systems of nonlinear equations, as well as on the software organization. The main principles of the moment sheme of finite elements (MSFE) is the base for procedure of obtaining of finite-element solving relationships. In contrast to the generally accepted approaches the use of MSFE provides an independent representation of deformations in the form of Maclaurin series in addition to the choose of the law of displacements distribution. At the base of the algorithm for solving a system of nonlinear equations, the method of integration over the load parameter is adopted. It is assumed to reduce an integration step consecutively to obtain reliable results. To optimize discrete FEM models, an approach based on the initial fragmentation of the computational model was used. The accuracy of solving of a system of nonlinear equations at each step of the load parameter is determined by comparing the sum of the squares of the nodal reactions and the sum of the squares of the nodal values of the external loads. The equations of termoviscolastoplastic deformation taking into account the damage of materials are taken as initial relations. The principle of using of quasiregular fragmentation which involves constructing a general irregular finite element model based on regular finite element fragments is basis of the procedure of constructing the finite element model of objects under consideration. This make it possible to optimize substantially the total number of unknowns. When calculating of a new objects, the convergence of the results is justified by a sequential increase of the parameters of the mesh model, a decrease of integration step value due to an increase of their number within a given range of loads, and an increase of the accuracy of solving of nonlinear equations system at each step. In order to justify the reliability of the results, obtained on the basis of the universal finite element, the parameters of the stress-strain state and the damage parameter ω were compared with data calculated on the basis of finite elements of a general type with numerical integration and elements with integration in a closed form.