Feszültségmezőre épülő magasabb rendű végeselem-modellek a szilárd testek mechanikájában = Stress-based higher-order finite element models in computational solid mechanics

Többmezős variációs elvek alkalmazásával három különböző hp-verziós végeselem-modell került kifejlesztésre. Háromdimenziós, lineáris rugalmasságtani feladatok megoldására a Fraeijs de Veubeke-féle kétmezős variációs elven alapuló, az elsőrendű feszültségfüggvényeket és a szögelfordulásokat közvetlenül approximáló hexaéder elemet és egy arra épülő végeselemes programot dolgoztunk ki. Forgásszimmetrikus körhengerhéj-feladatok megoldására két különböző, dimenzió szerint redukált héjmodellt és hp-verziós héjelemet dolgoztunk ki, melyeket egy-egy C++ nyelven megírt végeselemes programba építettünk bele. Az egyik elem-modell a Fraeijs de Veubeke-féle kétmezős variációs elven, a másik pedig a Hellinger-Reissner-féle hárommezős variációs elven alapul. A hárommezős elvben a nem a priori szimmetrikus feszültségmező és a forgásmező mellett az elmozdulásmező is függetlenül approximálandó ismeretlenként jelenik meg. Kimutattuk és numerikusan igazoltuk, hogy a kifejlesztett elem-modellekkel kapott numerikus megoldások konvergencia-sebessége független a Poisson tényezőtől és a héj vastagságától, vagyis erősen inkompresszibilis anyagok és vékony héjak esetében is numerikus konvergencia problémák nélküli megoldásokat adnak a mérnöki alkalmazások szempontjából elsőrendű fontosságú feszültségmezőre nézve is. | Three different hp-version finite element model has been developed using multi-field dual-mixed variational principles. The hexahedron element developed for three-dimensional linear elasticity problems is based on the two-field variational principle of Fraeijs de Veubeke with independent approximations of the six independent components of the first-order stress function tensor and three continuum rotations. The element and the corresponding research code, written in the C++ programming language, can be used for three-dimensional problems with simple geometry. Two other dimensionally reduced shell elements have also been developed, both of them for axisymmetric cylindrical shell problems: one is based on the two-field variational principle of Fraeijs de Veubeke, the other has been developed using the three-field dual-mixed principle of Hellinger-Reissner with independently approximated variables of not a priori symmetric stresses, continuum rotations and the displacements. The locking-free property of the elements has been demonstrated through test examples, and the solutions have been compared to analytical solutions as well as to results obtained using standard displacement-based elements. It has been shown that the dual-mixed approach and the related finite elements lead to robust, reliable and accurate stress computations for not only higher order p-, but also for low order h-type approximations

Többmezős variációs elvek és végeselem-modellek a nemlineáris rugalmasságtanban

Egy feszültségmezőn és forgásmezőn alapuló általános, nemlineáris héjelmélet került kidolgozásra. A dimenzió szerinti redukálás variációs hátterében a kiegészítő virtuális munka elv általánosított alakja áll. Az I. Piola-Kirchhoff feszültségvektorokra vonatkozó transzlációs egyensúlyi egyenletek a priori kielégítése két elsőrendű feszültségfüggvény-vektor bevezetésével, a forgási egyensúlyi egyenletek kielégítése pedig a forgásmező, mint Lagrange-multiplikátor alkalmazásával történt. A klasszikus héjelméletek alapjául szolgáló kinematikai és feszültségi hipotézisek nem kerültek alkalmazásra, így a háromdimenziós anyagegyenletek módosítások nélkül vehetők figyelembe.A kidolgozott héjelmélet alapján egy feszültségmezőn és forgásmezőn alapuló új hp-verziós lemezelem-modell és végeselem-program került kifejlesztésre lineárisan rugalmas lemezfeladatok megoldására. Numerikus úton igazolásra került, hogy a lemezmodell aszimptotikusan korrekt és modellezési hibája lényegesen jobb, mint a klasszikus Kirchhoff- és Reissner-Mindlin lemezmodelleké. A kifejlesztett hp-verziós elem-modell rendkívül vékony és összenyomhatatlan anyagú lemezek esetén is numerikus problémák (shear locking és incompressibility locking) nélküli feszültségeket eredményez, nemcsak magasabb rendű p-, hanem alacsony rendű h-típusú approximációnál is. | Starting from the three-dimensional principle of complementary virtual work, a stress-based dimensionally reduced shell model has been derived. The three-dimensional translational equilibrium equations, written in terms of the first Piola-Kirchhoff stress vectors, are a priori satisfied using two first-order stress function vectors. The rotational equilibrium equations are enforced in a weak sense using the rotations as Lagrangian multipliers. No conventional kinematical assumptions have been applied and the shell model uses unmodified three-dimensional constitutive equations.Based on the shell model derived, a new hp-version finite element model and a finite element code have been developed for linearly elastic plates (including the membrane problem). The numerical results show that rates of convergence in the energy norm and in the stress computations are practically independent of the Poisson ratio and thickness of the plate, i.e., the elements are free from both incompressibility locking and shear locking. The stress computation is robust, reliable and accurate for not only higher order p- but also for low order h-type approximations, even for extremely thin plates