Řešení axiálně zatížených prutů pomocí vlastní implementace MKP s využitím kvadratického prvku
Loading...
Date
Authors
Janáčik, Lukáš
ORCID
Advisor
Referee
Mark
A
Journal Title
Journal ISSN
Volume Title
Publisher
Vysoké učení technické v Brně. Fakulta strojního inženýrství
Abstract
V tejto bakalárskej práci popisujeme algoritmus tvorby konečno-prvkového modelu axiálne zaťaženého prútu, pri použití kvadratického prvku. V úvode definujeme základné pojmy pružnosti a pevnosti, ktoré v priebehu práce používame, a ktorých znalosť je nevyhnutnou podmienkou pochopenia problémov, pre ktoré bola metóda konečných prvkov odvodená. Práca popisuje prechod od základnej diferenciálnej rovnice k slabej formulácií, na ktorej je táto numerická metóda postavená. Definuje prvkové matice a vektory a prechod od jedného prvku ku globálnym maticiam a vektorom, vztiahnutým k celému telesu. Popisuje implementovanie okrajových podmienok do metódy a postprocessing výsledkov, pre výpočet ďalších neznámych. V praktickej časti definujeme 3 demonštračné úlohy, ktorých výsledky z vlastnej implementácie MKP v programe Matlab, porovnávame s analytickým riešením a riešením z komerčného softvéru, ANSYS Workbench. Úlohy sú zamerané na zložité líniové zaťaženie prútu, meniaci sa priečny prierez a staticky neurčitú úlohu. Výsledky posuvov a normálových napätí naprogramovaného riešiča s kvadratickým prvkom sa zhodujú s výsledkami z Ansysu (pri totožných nastaveniach) a po navýšení počtu prvkov (podľa potreby) i s analytickým riešením. Príklad s premenlivým priečnym prierezom bol v Ansyse riešený rovinnou napätosťou. Bola tak zistená prítomnosť šmykových napätí v častiach priečneho prierezu vzdialených od strednice prútu, ktoré náš riešič neuvažuje a mohli by byť v určitých prípadoch významné.
This bachelor thesis describes an algorithm for programming finite-element model with quadratic elements for axial loaded bar. In the introduction, we define the basic concepts of mechanics of materials, which are used in this thesis and are necessary to understand problems, for which finite element method was formulated. The thesis clarifies a transition from basic differential equation to weak formulation, which is the base of finite element method. We define element matrices and describe transition to global matrices, relating to the whole body. Then describe implementation of boundary conditions and postprocessing of the results, necessary for calculation and displaying of other unknowns. In the practical part, 3 illustrative problems are presented and calculated numerically in FEM solver using Matlab, analytically and in software ANSYS Workbench. Results are then compared and evaluated. Problems have different boundary conditions (linear axial load, tempered cross section, statically indeterminate fixation). Results of displacement and normal stress for programmed solver are identical to those from Ansys (using the same settings) and analytical solution (after more elements are added, if necessary). Problem with tempered cross section was simulated in Ansys using plain stress, because the program can’t define bar with tempered cross section. This revealed sheer stress contained in parts of cross section further from centreline, which are not calculated in our FEM solver and in some cases might be significant.
This bachelor thesis describes an algorithm for programming finite-element model with quadratic elements for axial loaded bar. In the introduction, we define the basic concepts of mechanics of materials, which are used in this thesis and are necessary to understand problems, for which finite element method was formulated. The thesis clarifies a transition from basic differential equation to weak formulation, which is the base of finite element method. We define element matrices and describe transition to global matrices, relating to the whole body. Then describe implementation of boundary conditions and postprocessing of the results, necessary for calculation and displaying of other unknowns. In the practical part, 3 illustrative problems are presented and calculated numerically in FEM solver using Matlab, analytically and in software ANSYS Workbench. Results are then compared and evaluated. Problems have different boundary conditions (linear axial load, tempered cross section, statically indeterminate fixation). Results of displacement and normal stress for programmed solver are identical to those from Ansys (using the same settings) and analytical solution (after more elements are added, if necessary). Problem with tempered cross section was simulated in Ansys using plain stress, because the program can’t define bar with tempered cross section. This revealed sheer stress contained in parts of cross section further from centreline, which are not calculated in our FEM solver and in some cases might be significant.
Description
Citation
JANÁČIK, L. Řešení axiálně zatížených prutů pomocí vlastní implementace MKP s využitím kvadratického prvku [online]. Brno: Vysoké učení technické v Brně. Fakulta strojního inženýrství. 2022.
Document type
Document version
Date of access to the full text
Language of document
sk
Study field
Základy strojního inženýrství
Comittee
Ing. Lubomír Junek, Ph.D. (předseda)
doc. Ing. Jana Horníková, Ph.D. (místopředseda)
Ing. Petr Skalka, Ph.D. (člen)
Ing. Lubomír Houfek, Ph.D. (člen)
Ing. Petr Vosynek, Ph.D. (člen)
Ing. Oldřich Ševeček, Ph.D. (člen)
doc. Ing. Tomáš Návrat, Ph.D. (člen)
doc. Ing. František Šebek, Ph.D. (člen)
Date of acceptance
2022-06-20
Defence
Student ve vymezeném čase prezentoval svoji bakalářskou práci, dále byly přečteny posudky a zodpovězeny dotazy oponenta. Poté byly členy komise položeny následující otázky související s bakalářskou prací:
Jak byl v práci definovaný mezní stav pružnosti a mezní stav deformace? Jaký je důvod použití bázových funkcí vyššího řádu? Jakým způsobem by se v řešených úlohách projevilo teplotní zatížení? Co je to slabá formulace?
Po zodpovězení všech dotazů byla obhajoba celkově hodnocena jako výborná / A.
Result of defence
práce byla úspěšně obhájena
Document licence
Standardní licenční smlouva - přístup k plnému textu bez omezení