發(fā)布時(shí)間：2024-05-18 07:11

　　提出一種基于分段時(shí)域自適應(yīng)算法和等幾何分析的求解粘彈性問題的數(shù)值方法。利用時(shí)域分段展開,建立了遞推格式的比例邊界元求解方程,環(huán)向比例邊界采用等幾何技術(shù)離散,在繼承常規(guī)比例邊界有限元半解析、便于處理應(yīng)力奇異性/無限域問題等優(yōu)點(diǎn)的同時(shí),可更準(zhǔn)確地描述幾何邊界,由此進(jìn)一步提高了計(jì)算精度;在時(shí)域,通過分段時(shí)域自適應(yīng)計(jì)算,保證不同時(shí)間步長(zhǎng)下的計(jì)算精度。通過數(shù)值算例,從計(jì)算精度、收斂性等方面,對(duì)所提方法的有效性進(jìn)行了驗(yàn)證。

【文章頁數(shù)】：11 頁

【部分圖文】：

圖8帶裂紋平板Fig.8Aplatewithacrack表2不同時(shí)間步長(zhǎng)下A點(diǎn)IG-SBFEM應(yīng)力解的比較Table2ComparionforIG-SBFEMdisplacementsolutions

工程力學(xué)29圖7通過A點(diǎn)的位移解，描述了誤差限對(duì)展開階次的影響。計(jì)算結(jié)果顯示，與8=110相比，12=110的誤差限需要更多的展開階數(shù)，可見所提方法可根據(jù)不同的精度要求自適應(yīng)地調(diào)整展開階數(shù)。圖7展開階次隨時(shí)間的變化Fig.7Thevariationofexpansionorder....

圖9本文算法和時(shí)域非自適應(yīng)算法的比較Fig.9Thecomparisonoftheproposedmodelwithnonadaptivealgorithmintimedomain

30工程力學(xué)圖9本文算法和時(shí)域非自適應(yīng)算法的比較Fig.9Thecomparisonoftheproposedmodelwithnonadaptivealgorithmintimedomain(a)A點(diǎn)y(b)B點(diǎn)y(c)C點(diǎn)y圖10應(yīng)力隨時(shí)間變化曲線Fig.10Thevaria....

圖7展開階次隨時(shí)間的變化Fig.7Thevariationofexpansionorderwithtime

工程力學(xué)29圖7通過A點(diǎn)的位移解，描述了誤差限對(duì)展開階次的影響。計(jì)算結(jié)果顯示，與8=110相比，12=110的誤差限需要更多的展開階數(shù)，可見所提方法可根據(jù)不同的精度要求自適應(yīng)地調(diào)整展開階數(shù)。圖7展開階次隨時(shí)間的變化Fig.7Thevariationofexpansionorder....

圖3帶有1/4圓孔的半無限長(zhǎng)粘彈性體區(qū)域Fig.3Asemi-unboundedviscoelasticplate

法[23]，對(duì)于常位移邊界條件，將與邊界相關(guān)的控制點(diǎn)分成兩組，在u上恒為零的zziN和不恒為零的nziN，有：uunznzzzzziiiiuNuNu(57)對(duì)式(57)的非負(fù)基函數(shù)項(xiàng)利用NURBS基函數(shù)的單位分解性，將求解域上的邊界約束條件轉(zhuǎn)化到控制點(diǎn)上：zziuu(58)與上式....

本文編號(hào)：3976627