二維三維彈性問(wèn)題混合元.pdf

1、本文主要討論線(xiàn)彈性方程,在Hellinger-Reissner變分形式的基礎(chǔ)上,系統(tǒng)的構(gòu)造了二維空間下的矩形和三角形單元,三維空間下的立方體和四面體單元等一系列簡(jiǎn)單穩(wěn)定的單元.對(duì)單元的適定性,收斂性,誤差估計(jì),以及二維三維矩形和立方體單元的各向異性特征進(jìn)行了深入的分析和系統(tǒng)的研究.并對(duì)二維協(xié)調(diào)的矩形單元和非協(xié)調(diào)的三角形單元進(jìn)行了相應(yīng)的數(shù)值實(shí)驗(yàn)。首先構(gòu)造了具有最少自由度的協(xié)調(diào)的二維矩形單元R8?2和三維立方體單元C18?3.即矩形單元的應(yīng)

2、力和位移空間分別為8個(gè)和2個(gè)自由度.立方體單元的應(yīng)力和位移空間中分別為18和3個(gè)自由度.由于所構(gòu)造單元不滿(mǎn)足關(guān)于散度的投影性質(zhì),因此我們采用了構(gòu)造的方法證明了離散BB條件,即離散混合問(wèn)題的唯一可解性條件.在此基礎(chǔ)上進(jìn)一步分析,發(fā)現(xiàn)了單元的各向異性特征,并由此得到了單元的誤差估計(jì).據(jù)我們所知這是首次構(gòu)造的各向異性彈性問(wèn)題混合元。其次,在構(gòu)造最簡(jiǎn)單矩形單元的基礎(chǔ)上,進(jìn)一步構(gòu)造了一系列矩形高階單元,當(dāng)次數(shù)大于等于4時(shí)滿(mǎn)足散度的投影性質(zhì),據(jù)此得

3、到單元的適定性和離散問(wèn)題的唯一可解性及誤差估計(jì),并得到相應(yīng)的彈性復(fù)形。再次,在構(gòu)造矩形類(lèi)單元的基礎(chǔ)上,又對(duì)三維四面體單元構(gòu)造過(guò)程中的空間?k(K)進(jìn)行了研究和討論,因?yàn)樵贖ellinger-Reissner變分形式下,應(yīng)力是屬于H(div,Ω;S)空間,協(xié)調(diào)元構(gòu)造要求應(yīng)力的法向分量跨過(guò)單元邊界連續(xù),在協(xié)調(diào)元構(gòu)造中計(jì)算空間此處為公式的維數(shù)是個(gè)難點(diǎn)問(wèn)題.在本部分中給出了任意階空間維數(shù)計(jì)算的一般方法,并且此方法能夠很方便的求出相應(yīng)的顯式基,同

4、時(shí)證明了k=3時(shí),集合?k(K)的維數(shù)為0,并且給出k=4時(shí)空間的一組基。最后,構(gòu)造了關(guān)于線(xiàn)彈性問(wèn)題的一系列新的從低階到高階的三角形和四面體非協(xié)調(diào)單元及相應(yīng)剛體運(yùn)動(dòng)下的簡(jiǎn)化單元,這里構(gòu)造的二維三維單元區(qū)別于之前文獻(xiàn)中構(gòu)造的單元,構(gòu)造簡(jiǎn)單,自由度少.這些單元定義在參考單元上,形函數(shù)空間顯式給出,我們嚴(yán)格證明了這類(lèi)單元的仿射等價(jià)性,易于進(jìn)行數(shù)值實(shí)驗(yàn)當(dāng)k=1時(shí)簡(jiǎn)化的三角形單元的應(yīng)力空間和位移空間具有12+3個(gè)自由度,簡(jiǎn)化的四面體單元的應(yīng)力空間