A modern mérnöki tervezés és elemzés alapköve a végeselemes analízis (VEA), amely komplex geometriák és terhelési feltételek melletti anyagviselkedés megértésére szolgál. Azonban a valósághű modellek rendkívül számításigényesek lehetnek, ezért kulcsfontosságú az egyszerűsítések alkalmazása, amelyekkel a pontosság megtartása mellett jelentősen csökkenthető a számítási idő. Az alábbiakban részletesen bemutatjuk ezeket az egyszerűsítési módszereket, különös tekintettel az axiális és hajlító terhelések vizsgálatára cső- és rúdszerkezetek példáján keresztül, kitérve a különböző modellezési megközelítésekre és azok eredményeire.
A Végeselemes Analízis (VEA) Alapjai és Egyszerűsítései
A végeselemes analízis a mérnöki fizika egyik legerősebb numerikus eszköze, amely összetett szerkezetek mechanikai viselkedésének, például feszültségeloszlásának, alakváltozásának és hőátadásának előrejelzésére szolgál. Lényege, hogy a vizsgált geometriát apró, véges számú elemekre osztja (hálózás), amelyek egymással csomópontokon keresztül kapcsolódnak. Ezeken az elemeken belül a fizikai törvényeket egyszerűsített matematikai függvényekkel írják le, majd az egész rendszerre vonatkozóan megoldják az egyenleteket. Ezáltal a mérnökök virtuálisan tesztelhetik a tervezett szerkezetek működését még a fizikai prototípusok gyártása előtt, optimalizálva a kialakítást és növelve a biztonságot.
Azonban a részletes, háromdimenziós (3D) modellek rendkívül nagy számítási kapacitást és időt igényelnek. Ezért a gyakorlatban gyakran alkalmaznak egyszerűsítéseket. Az egyszerűsítések csoportosítása többféleképpen történhet, például a szimmetriák típusa, a mechanikai probléma dimenziója vagy az alkalmazott elemtípus szerint. Ezekkel, illetve az egyszerűsítések okaival, a geometriai, perem- és terhelési feltételeket érintő követelményekkel, valamint az alkalmazott elemtípusokkal már foglalkoztunk korábban, most nézzünk néhány egyszerű gyakorlati példát ezek alkalmazására.
Az egyszerűsítések célja, hogy a valós problémát egy olyan matematikai modellé alakítsuk, amely lényegesen kevesebb véges elemet és csomópontot tartalmaz, miközben az eredmények továbbra is elfogadható pontosságúak maradnak a mérnöki célokra. Ez magában foglalhatja a geometria csökkentését (pl. 3D-ből 2D-re vagy 1D-re), a szimmetria kihasználását, vagy speciális elemtípusok alkalmazását, amelyek hatékonyabban modellezik bizonyos szerkezeti viselkedéseket. A helyes egyszerűsítés kiválasztása kritikus a sikeres és hatékony VEA szimulációhoz.
Csőszerkezetek Vizsgálata Belső Nyomással - Szimmetria Alkalmazása
A csőszerkezetek, különösen azok, amelyek belső nyomás alatt állnak, gyakori elemei a gépészeti rendszereknek. Ezek elemzése során a szimmetria-feltételek alkalmazása jelentősen leegyszerűsítheti a modellt és gyorsíthatja a számításokat. Az alábbiakban egy belső nyomással terhelt cső modelljén keresztül mutatjuk be a szimmetria-feltétel alkalmazását 3D-s test-, illetve héjmodell esetén, valamint a 2D-s sík alakváltozás, illetve tengelyszimmetrikus modell használatát.
3D-s Testmodell Szimmetriafeltételekkel
A teljes, körbefutó 3D-s csőmodell szimulációja rendkívül sok végeselemet igényelne, különösen akkor, ha a vastagság és a hossza is nagy, és finom hálóra van szükség. Azonban, ha a geometria és a terhelés is szimmetrikus, akkor elegendő a cső egy részét modellezni. Esetünkben, ha a cső tengelyvonala egybeesik a globális koordinátarendszer Y tengelyével, akkor az XY és YZ síkokat kell kiválasztanunk a modellfából a szimmetriasíkok kijelölésére. Ebben az esetben két szimmetriasík van (de akár három is lehetne, ha a hossza mentén is kettévágnánk).
Amikor 3D-s testmodellt használunk, a megtámasztás a szimuláció típusától függően változó. A szimmetria kényszer definiálásának egyik módja, ha eleve negyed modellt hozunk létre, és a testmodell megfelelő felületeihez hozzárendeljük a kényszert, akár a geometriamodellező, akár a szimulációs felületen. Ha negyed modellünk van, akkor egyszerű kinematikai kényszerek segítségével is definiálhatjuk a szimmetriát, mégpedig olyan elmozdulás kényszer segítségével, amely csak a felület síkjában történő elmozdulást engedi, és az arra merőleges irányút nem. Hasonló módon kötjük meg a két véglap axiális elmozdulását a radiális szabadságfokok engedésével. A modell belső felületére (például a 6.2. ábra, b. pontjánál) kell a terhelést, például a belső nyomást, definiálni.
A tengelyirányú feszültségkomponenst a peremfeltételektől függően tudjuk számolni. Egy vékonyfalú cső esetén azzal a feltételezéssel élünk, hogy a cső anyagában csak tangenciális és axiális feszültség ébred, tehát a cső falában ébredő radiális feszültségkomponens 0. Ez az egyszerűsítés jelentősen megkönnyíti az analitikai számításokat, de a végeselemes modellezés lehetővé teszi a radiális feszültség valósághűbb megjelenítését is, még vastagfalú csövek esetén is.
A modell hálózása után, például 1mm-es méretű tégla elemeket tartalmazó rendezett struktúrájú hálóval (lásd 6.2. ábra, a.), futtatjuk a szimulációt. A szimuláció futtatása után lekérdezzük az egyenértékű feszültséget, a fajlagos nyúlást és a radiális deformációt. Az eredmények vizuális megjelenítése, például színskálás ábrákkal, segíti a feszültségkoncentrációk és deformációs mintázatok gyors azonosítását.
Héjmodell (Shell Model) Előnyei és Alkalmazása
A héjmodellek használata rendkívül előnyös vékonyfalú szerkezetek, mint például csövek vagy tartályok elemzésekor. A feladat héjmodellel történő megoldásának előnye, hogy viszonylag kevés véges elemet tartalmaz a modell, de mégis mindhárom irányú terheléseket és eredményeket is meg tudunk jeleníteni a segítségével. Míg egy 3D-s testmodell vastagsági irányban több elemet igényel, a héjmodell lényegében egy felületi hálóval dolgozik, amihez vastagságot és anyagjellemzőket rendelünk, így drasztikusan csökkentve az elemszámot és a számítási időt.
Szimmetria síkoknak itt is az XY és YZ síkokat választhatjuk. Ha például a henger alakú héj átmérője 10mm, a falvastagság pedig viszonylag kicsi, a héjmodell ideális választás. Az analízis típusa ebben az esetben is 3D-s, így az alapbeállítások maradhatnak érvényben, de héjmodell esetén definiálni kell a felülettől való eltolás irányát is, ami a neutrális sík elhelyezkedése szempontjából fontos.
A hálózásnál 1mm-es elemméret esetén például 80 elemet tartalmazó modell alakulhat ki (lásd 6.5. ábra). A cső két végét elmozdulás kényszer segítségével rögzítjük, melynél az axiális (Y) irányú elmozdulást 0 értéken rögzítettük, míg a henger tengelyére merőleges (X és Z) irányokban a test szabadon elmozdulhat. A szimmetriasíkok mentén az adott síkban történő mozgást megengedő és a síkra merőleges irányú elmozdulást gátló kényszert kell definiálnunk. A belső nyomást terheléssel definiáltuk (lásd 6.6. ábra).
A szimuláció lefuttatása után, az előző példához hasonlóan lekérdezzük az egyenértékű feszültséget, a fajlagos nyúlást és a radiális deformációt. Az eredményeken látható, hogy az értékek változása a vastagság mentén eltér az előző szimuláció eredményétől, ami a héjmodellek inherens egyszerűsítéséből adódik, és a vékonyfalú szerkezetekre jellemző viselkedést tükrözi. A héjmodellek különösen alkalmasak arra, hogy gyorsan és hatékonyan kapjunk globális képet a feszültségekről és alakváltozásokról, de a vastagság menti részletes feszültségeloszlás vizsgálatára korlátozottan képesek.
2D-s Sík Alakváltozás (Plane Strain) Szimuláció
Bizonyos esetekben a probléma jellege lehetővé teszi a 2D-s modellezést, ami tovább csökkenti a számítási terhelést. A 2D-s sík alakváltozás típusú szimulációk lényege, hogy a test vizsgált keresztmetszetére merőleges (például ANSYS-ban: Z tengely) irányában nincs alakváltozás. Ez az állapot hosszú, állandó keresztmetszetű (a vizsgált keresztmetszetre merőleges kiterjedése nagyobb, mint a másik két irányban), és a hossz mentén állandó és a keresztmetszet síkjába eső terhelési és megfogási feltételekkel rendelkező testekre jellemző. Ennél fogva a feszültségi állapot és a keresztmetszet síkjába eső alakváltozási állapot is azonos a hossz (Z tengely) mentén, tehát elegendő egyetlen keresztmetszet 2D-s vizsgálata.
Mivel feltételeztük, hogy a cső keresztmetszetéhez képest hosszú, tehát a tengelyirányú (Z) kiterjedése jelentősen nagyobb, mint a radiális irányú (X és Y) kiterjedése, valamint tengelyirányú terhelés nincs jelen, a feladat megoldható 2D-s sík alakváltozás típusú szimulációval. Emellett a cső geometriája miatt szimmetria feltételeket is alkalmazhatunk; elegendő a negyed geometria felépítése. A geometria tehát a cső egy negyed szelete, melyet egy sík felületmodellként (Surface Body) definiálunk a koordinátarendszer XY síkján. Az egyszerűség kedvéért a fősíkok szolgáltatják a szimmetria síkokat.
A 2D-s viselkedést a szimulációs környezetbe való belépés előtt be kell állítanunk, mert utólagosan már nincs lehetőség az átjárásra. A szimuláció során szerkezeti acél anyagot állítunk be. A véges elemek mérete például 0,5mm (lásd 6.10. ábra). A modell így 64 elemet tartalmaz, ami viszonylag kevés, azonban ez olyan pontos eredményt fog adni, mintha a teljes 3D-s geometriát modelleznénk, megfelelő vastagsággal.
A szimmetriafeltételeknél mindig figyelni kell arra, hogy a megfogások megfeleljenek a valóságnak. Jelen esetben a valós viselkedés modellezése az X és Y tengelyre illeszkedő két él adott tengelyre merőleges elmozdulásának megkötése. Ez egy-egy elmozdulás típusú kényszer segítségével oldható meg, melyek az illeszkedő tengely irányába eső elmozdulást szabadon engedik, míg a tengelyre merőlegest gátolják (lásd 6.11. ábra). A szerkezet terhelése a cső belső falán ébredő 100MPa nagyságú nyomás (lásd 6.11. ábra).
A szimuláció futtatása után lekérdezzük az egyenértékű feszültséget, a fajlagos nyúlást és a radiális deformációt. Ez utóbbihoz egy polár koordináta rendszert kell definiálnunk a megfelelő kiértékeléshez. Az alábbi ábrán láthatjuk az ébredő egyenértékű feszültség eloszlását a cső fala mentén (lásd 6.12. ábra). Ez a 2D-s megközelítés rendkívül hatékony a hosszú, egyenletes profilú szerkezetek elemzésére, ahol a hosszanti irányú hatások elhanyagolhatóak vagy konstansak.
2D-s Tengelyszimmetrikus (Axisymmetric) Modell
Ugyanezt a feladatot meg lehet oldani 2D-s tengelyszimmetrikus modellként is, mivel mind a geometria, mind a terhelés tengelyszimmetrikus. Ebben az esetben a geometria azt a fél hosszmetszetet jelenti, amelyet körbeforgatva megkapjuk a teljes geometriát, tehát a csövet modellező feladatnál ez egy 10mm * 2mm-es téglalap a forgatási tengelytől 4mm-es minimális távolságra.
A 2D-s tengelyszimmetrikus modelleket úgy kell felépíteni, hogy a forgástengely a globális Y tengely legyen, a geometria pedig az 1. és a 4. sík negyedben (+X) helyezkedjen el (lásd 6.13. ábra). Mivel az analízis típusa 2D tengelyszimmetrikus, a 2D-s viselkedést, az előző feladathoz hasonlóan, a szimulációs környezetbe való belépés előtt be kell állítanunk (lásd 6.10. ábra).
Az anyag hozzárendelése és a szimuláció általános beállításai után készítsük el a végeselemes hálót. A háló, ahogy az előző esetben is, 0,5mm nagyságú négyszög elemekből áll. Az így létrehozott háló összesen 80 elemet tartalmaz (lásd 6.14. ábra). Ez a fajta modellezés különösen előnyös olyan szerkezeteknél, mint a nyomástartó edények, forgó alkatrészek vagy gyűrű alakú komponensek, ahol a geometria és a terhelés egy tengely körüli forgatással ismétlődik.
A megfogások definiálása itt is különös figyelmet igényel. A csőmodell két végén olyan elmozdulás típusú kényszert célszerű alkalmazni, amely megengedi a cső tágulását, tehát radiális (X) irányú elmozdulását, de leköti az axiális (Y) irányú elmozdulás lehetőségét. Mivel ebben a modellben a két véglap párhuzamos egymással, elegendő egyetlen megfogást megadni, mert mindkét oldal ugyanolyan irányban mozoghat, illetve van lekötve (lásd 6.15. ábra). A 100MPa-os belső nyomás a csőfelület belső felületén, tehát itt a tengelyhez közel álló párhuzamos élen van (lásd 6.16. ábra).
A szimuláció lefuttatása után, az előző példához hasonlóan lekérdezzük az egyenértékű feszültséget, a fajlagos nyúlást és a radiális deformációt. Ez utóbbihoz most nem kell poláris koordináta rendszert definiálnunk, mivel az az X tengely irányába esik. Az alábbi ábrán láthatjuk az ébredő egyenértékű feszültség eloszlását a cső fala mentén (lásd 6.17. ábra). A tengelyszimmetrikus modellezés hihetetlenül hatékony, és rendkívül pontos eredményeket ad, amennyiben a geometriai és terhelési feltételek megfelelnek az elméleti alapoknak.
Az Egyszerűsítések Eredményeinek Összehasonlítása
Miután négyféle módon elvégeztük a szimulációt, azonos feltételekkel (azonos anyag és geometria, azonos elemméret és szimulációs beállítások), láthatjuk, hogy a héjmodell kivételével az eredmények nagymértékben egyeznek az analitikus számolással, melyet a vastag csövekre alkalmazott formulákkal végeztünk. Ugyanakkor, a vékonyfalú csövekre alkalmazott analitikus számolás eredménye jó egyezést mutat a héjmodellel végzett szimulációval. Ez a megfigyelés alátámasztja az egyszerűsítések érvényességét és fontosságát a megfelelő körülmények között.
A 3D-s testmodell nyújtja a legáltalánosabb és legpontosabb eredményeket, de a legdrágább is. A héjmodell hatékonyabb vékonyfalú szerkezetek esetén, de a vastagság menti feszültségeloszlás részleteit nem adja meg. A 2D-s sík alakváltozás és tengelyszimmetrikus modellek pedig rendkívül gyorsak és pontosak, ha a probléma jellege lehetővé teszi alkalmazásukat. A helyes modellezési stratégia kiválasztásához elengedhetetlen a vizsgált szerkezet viselkedésének, a terhelési és peremfeltételeknek, valamint a kívánt eredmények pontosságának alapos megértése.
Rúdszerkezetek Hajlítási Vizsgálata - Konzol Példáján
A rúdszerkezetek, mint például gerendák vagy konzolok, a mérnöki szerkezetek alapelemei. Ezek viselkedését, különösen hajlító igénybevétel esetén, szintén vizsgálni lehet végeselemes analízissel, különböző egyszerűsítési szintekkel. Az alábbiakban egy konzolos rúd példáján keresztül mutatjuk be a 3D-s test-, héj-, 2D-s sík feszültség/alakváltozás-, valamint vonalmodell (gerenda) alkalmazását.
Teljes 3D-s Geometria Szimulációja
A 3D-s geometria használata általában a leghosszadalmasabb számítási igényű eljárás, különösen, ha finom hálózatra van szükség a feszültségkoncentrációk pontos rögzítéséhez. Tehát célszerű (amennyiben lehetséges) az egyszerűsített eljárásokat használni. A teljesség kedvéért azonban bemutatjuk a teljes 3D-s geometria szimulációját is ugyanazon a példán.
Ha egy hasábot modellezünk, például egyetlen kihúzás (Extrude) segítségével, ahol a vázlatot egy 100x10mm-es téglalap alkotja, és a kihúzás mértéke normál irányba 2mm (lásd 6.22. ábra), ez adja a teljes 3D-s geometriát. Az általános beállítások után (anyag hozzárendelése, nagy elmozdulások engedélyezése stb.) készítsük el a végeselemes hálót téglaelemek felhasználásával (lásd 6.23. ábra).
Az eredmények értékeléséhez gyakran alkalmaznak konstrukciós útvonalat (Path), ami egy olyan térbeli görbét jelent, mely mentén diszkrét pontokban lekérdezhető az eredmény, és az eredmény pontonként táblázatban és diagramon is megtekinthető. A pontokat és értékeket tartalmazó táblázat exportálható abból a célból, hogy más programok segítségével is elemezhető legyen. Ennek definiálására több lehetőség is van: közvetlenül a kezdő és végpont megadásával, vagy a vonal megadható a pozitív x tengely és a háló modell metszéseként. A konstrukciós pálya lehet egyenes, vagy íves, attól függően, hogy milyen típusú koordináta rendszert használunk (Descartes-féle koordináta rendszernél egyenes, míg henger koordináta rendszernél íves). A konstrukciós vonal diszkrét pontjainak száma definiálható, de maximum 200 lehet. Másik lehetőség, amikor egy meglévő él kijelölésével definiáljuk a konstrukciós vonalat. Ebben az esetben a vonal alatt lévő csomópontok jelentik a vonal diszkrét pontjait. Az útvonal típusa lehet élre illeszkedő, melyet a testmodellünk egyik hosszanti élének kijelölésével rendelhetünk hozzá (lásd 6.24. ábra).
A terhelés egy fix megfogás a rúd egyik végén és egy koncentrált erő jellegű terhelés a rúd másik végén. Az erőt komponensenként megadva -Y irányban 8N nagyságú, míg X és Z irányba 0N nagyságú. A terhelést a 6.24. ábra szemlélteti.
A szimuláció futtatása után kérdezzük le a teljes (lásd 6.26. ábra), és az Y irányú lehajlás értékét a létrehozott konstrukciós vonal (Path) mentén (lásd 6.27. ábra). A 3D-s modell adja a legteljesebb képet a szerkezet viselkedéséről, beleértve a lokális feszültségkoncentrációkat is, de a legkomolyabb számítási erőforrást igényli.
Héjelemekkel (Shell Elements) Való Modellezés Konzolon
A héjelemekkel való szimulációt ebben az esetben is érdemes megvizsgálni. A geometria elkészülte után, lépjünk be a szimulációs környezetbe. Ezután végezzük el a hálózást. Ugyanúgy, mint az előző esetben 2mm-es elemméretet használtunk. Látható, hogy ez meglehetősen eltér az előző esettől, mivel a test vastagsága mentén csak egy elemre oszlik a geometria (lásd 6.30. ábra).
Ahogy az előbbi esetben itt is, hozzunk létre egy konstrukciós geometriát, mely segíti az eredmények értékelését. Ennek alapja (lásd 6.31. ábra) egy vonal, amely a rúd hosszanti tengelye mentén halad. A terhelés ugyanaz, mint az előző esetben (fix megfogás a test egyik véglapján, 8N erőterhelés -Y irányban a másik végén). Mivel itt felületmodellről van szó, a véglapok kijelöléséhez a felületet határoló él kijelölése szükséges (lásd 6.32. ábra).
A szimuláció futtatása után a kérdezzük le a teljes (lásd 6.33. ábra), és a konstrukciós vonal menti Y irányú deformációt (lásd 6.34. ábra). A héjmodell ebben az esetben egyszerűsítést jelent a 3D-s testmodellhez képest, és bár a vastagság menti részletes feszültségeloszlást nem adja meg, a globális lehajlást és a fő feszültségeket elfogadható pontossággal előrejelzi. Ez a megközelítés gyorsabb, mint a teljes 3D-s modellezés, és jó kompromisszumot jelent a sebesség és a pontosság között, ha a vastagság viszonylag kicsi a rúd más méreteihez képest.
2D-s Sík Feszültség (Plane Stress) és Sík Alakváltozás (Plane Strain) Analízis Konzolon
A 2D-s sík feszültség analízis esetén a test vizsgált metszetére merőleges (ANSYS-ban: Z tengely) irányában nem ébred feszültség. Ezzel szemben a 2D-s sík alakváltozás típusú szimulációk lényege, hogy a test vizsgált keresztmetszetére merőleges (ANSYS-ban: Z tengely) irányában nincs alakváltozás. Mivel ebben a példában a terhelés és megfogás is egyenletes a Z tengely mentén (és nincs Z tengely irányú komponense), mindkét eset megvalósítható. Viszont, mivel a valós geometria Z irányú mérete lényegesen kisebb, mint az X irányú hossza, és nagyobb mint az Y irányú magassága, ezért egyik eljárás geometriai feltételeit sem teljesíti igazán. Vagyis Z irányban sem a feszültség, sem pedig az alakváltozás értéke nem lesz teljesen 0. Viszont mivel elhanyagolható mértékűek a síkbeli komponensekhez képest, ezért jó közelítést jelentenek.
Vázlatkészítésnél a 100x2mm-es téglalap megrajzolását az XY síkon szükséges elvégezni. A 10mm-es vastagságot pedig a szimuláció során tudjuk megadni. A felület létrehozásának menete ugyanaz, mint a 3D-s Shell modelles példában, annyi különbséggel, hogy a felületnek mindenképp egy síkban (XY) kell lennie (lásd 6.35. ábra), és a szimuláció során 2D-s beállításokat kell alkalmazni az ott használt 3D-s helyett.
A geometria elkészítése és a szükséges beállítások (anyagdefiniálás, nagy alakváltozás engedélyezése, vastagság beállítása stb.) után készítsük el a végeselemes hálót, mely elemméretét állítsuk 1mm-re. Mivel itt jóval kevesebb véges elem adódik, mint például a 3D-szimulációk esetén, lehetőség van a sűrűbb háló használatára (lásd 6.36. ábra).
A terhelés hasonlóan az eddigi esetekhez egy fix megfogást és egy -Y irányú, 8N nagyságú erőterhelést jelent a rúd felső élén (lásd 6.37. és 6.38. ábrák). Eredményként, a bevezetésben megfogalmazott problémák ellenére, az eddigiekhez nagyon hasonló értékeket kaptunk. A teljes deformáció eredményét a 6.39. ábra, az Y irányú, vonal menti deformáció eredményét a 6.40. ábra mutatja. Ezek a 2D-s modellek jelentős számítási előnyökkel járnak, és gyakran elegendő pontosságot biztosítanak a rúd globális viselkedésének, például a lehajlásnak és a fő feszültségkomponenseknek a meghatározásához. Fontos azonban mindig felmérni, hogy a valós geometria és terhelés mennyire felel meg a sík feszültség vagy sík alakváltozás feltételeinek.
Vonalmodellek (Gerenda Elemek) Használata Rúd Szerkezetek Esetén
A vonalelemek, más néven 1D elemek, a leghatékonyabbak a rúdszerkezetek modellezésére, különösen, ha a rúd keresztmetszeti méretei kicsik a hosszához képest. Az ilyen elemek alkalmazását aszerint választhatjuk ki, hogy a vizsgált szerkezetben lévő vonal elemek felvesznek-e hajlító igénybevételt (gerenda), vagy csak tengelyirányú terhelést (rúd). Esetünkben, mivel a konzol tengelyére merőleges terhelés hajlításra veszi igénybe a tartót, egyértelműen a gerendamodellre van szükség.
A gerendamodell geometriája két részből áll: a tartó tengelyét képviselő vonalmodellből és a keresztmetszet definíciójából. A vezérgörbe a Z tengely mentén fekvő, 100mm hosszú vonal, a keresztmetszet pedig egy 10x2mm-es téglalap (lásd 6.41. ábra). A vonalmodell egy adott parancsával hozható létre, a már előzőleg megrajzolt vonal vázlatának kijelölésével (lásd 6.42. ábra).
A keresztmetszet definiálását egy másik parancsával lehet elvégezni (lásd 6.43. ábra). Lehetőségünk van az előre definiált, jellegzetes keresztmetszetek kiválasztására (ilyenkor csak a méreteket kell megadni), illetve saját vázlat készítésére. Miután a keresztmetszetet létrehoztuk, hozzá kell rendelni a vonalmodellhez. A keresztmetszet definiálásakor figyelni kell a vonalelem orientációjára (a zöld színű koordinátatengely jelöli a +Y irányt), és a láthatóságát (lásd 6.41. ábra) is ellenőrizni kell.
Miután a szimuláció beállításai megtörténtek (pl. anyag hozzárendelés), elkészítjük a végeselemes hálót. A testet 2mm-es elemmérettel hálóztuk be. Ahogy a 6.44. ábra mutatja, a vonalmodell mentén mindenhol csak 1 db elem helyezkedik el. Ez a rendkívül alacsony elemszám teszi a vonalelemeket a leggyorsabb modellezési módszerré rúdszerű szerkezetek esetében.
A szimuláció eredménye a konstrukciós segédvonal mentén lekérdezhető. Jelen feladatban a konstrukciós segédvonalat a test közepén végigfutó vonalmodell kijelölésével definiáljuk (lásd 6.45. ábra). A megfogások megadásakor a két végpontot jelöltük ki. A terhelést a 6.46. ábra mutatja, ahol piros, B jelű címkével a -Y irányú, 8N nagyságú erőterhelés, és a kék, A jelű címkével pedig a fix befogás van jelölve.
A gerenda elemek rendkívül hatékonyak a globális elmozdulások és a belső erők (hajlító nyomatékok, nyíróerők, axiális erők) meghatározására. Mivel ezek az elemek a gerendaelméleten alapulnak, amely feltételezi, hogy a keresztmetszet merev marad és síkban mozog, a lokális feszültségkoncentrációk részleteit nem tudják pontosan megadni. Azonban a nagyobb szerkezetek előzetes tervezéséhez és optimalizálásához kiválóan alkalmasak. Az "axiális rúd nyomaték érték" szempontjából ez a megközelítés közvetlenül számol hajlító nyomatékokkal, amelyek a terhelő erőkből adódnak, és a rúd keresztmetszetében ébredő hajlító feszültségeket eredményezik. A "rúd" fogalma itt "gerenda" jelentést kap, figyelembe véve a hajlító igénybevételt.
A Modellezési Döntések Jelentősége és Eredményvalidáció
Ahogy a bemutatott példák is demonstrálják, a végeselemes modellezés során az alkalmazott egyszerűsítések megválasztása rendkívül fontos. Nem csupán a számítási időt és erőforrásokat befolyásolja, hanem az eredmények pontosságát és értelmezhetőségét is. A legáltalánosabb 3D-s testmodellezés nyújtja a legnagyobb részletességet és pontosságot, de a legdrágább is. A 2D-s sík feszültség és sík alakváltozás, valamint a tengelyszimmetrikus modellek jelentősen gyorsítják a számításokat, ha a szerkezet geometriája és terhelése megfelel a feltételezéseknek. A héj- és vonalelemek pedig kifejezetten vékonyfalú vagy rúdszerű szerkezetekhez optimalizáltak, ahol a globális viselkedés (pl. lehajlás, hajlító nyomaték) a fő érdeklődési terület.
A választás során figyelembe kell venni a következőket:
- A szerkezet geometriája és méretarányai: Vékonyfalú? Hosszú és egyenes? Tengelyszimmetrikus?
- A terhelés típusa és eloszlása: Pontszerű? Egyenletesen elosztott? Axiális? Hajlító? Csavaró?
- A kívánt eredmények pontossága és részletessége: Globális deformációkra vagy lokális feszültségkoncentrációkra van szükség?
- Rendelkezésre álló számítási erőforrások és időkeret.
Fontos az is, hogy az egyszerűsítésekkel kapott eredményeket mindig validáljuk. Ez történhet analitikus számításokkal (ahogy a cső példájánál láttuk, hogy a vastag csövekre alkalmazott formulákkal nagymértékben egyeznek az eredmények, kivéve a héjmodellt, mely a vékonyfalú csövekre alkalmazott analitikus számolással mutatott jó egyezést), vagy kísérleti mérésekkel, vagy egy finomabb, 3D-s modellel végzett összehasonlítással. Csak így biztosítható, hogy a modell valóban a valós szerkezeti viselkedést tükrözze, és a mérnöki döntések megalapozottak legyenek. A különböző modellek közötti eredmények értelmezése és a különbségek okainak megértése alapvető fontosságú a tapasztalt végeselemes elemző számára.