


下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、Lagrange、Euler、ALE三種方法的簡(jiǎn)單介紹ALE、Lagrange、Euler是數(shù)值模擬中處理連續(xù)體的廣泛應(yīng)用的三種方法。Lagrange方法多用于固體結(jié)構(gòu)的應(yīng)力應(yīng)變分析,這種方法以物質(zhì)坐標(biāo)為基礎(chǔ),其所描述的網(wǎng)格單元將以類似雕刻”的方式劃分在用于分 析的結(jié)構(gòu)上,即是說采用Lagrange方法描述的網(wǎng)格和分析的結(jié)構(gòu)是一體的,有限元節(jié)點(diǎn)即為物質(zhì)點(diǎn)。采用這種方法時(shí),分析結(jié)構(gòu)的形狀 的變化和有限單元網(wǎng)格的變化完全是一致的(因?yàn)橛邢拊?jié)點(diǎn)就為物質(zhì)點(diǎn)),物質(zhì)不會(huì)在單元與單元之間發(fā)生流動(dòng)。這種方法主要的優(yōu)點(diǎn) 是能夠非常精確的描述結(jié)構(gòu)邊界的運(yùn)動(dòng),但當(dāng)處理大變形問題時(shí),由于算法本身特點(diǎn)的限制,將會(huì)出
2、現(xiàn)嚴(yán)重的網(wǎng)格畸變現(xiàn)象,因此不利 于計(jì)算的進(jìn)行。Euler方法以空間坐標(biāo)為基礎(chǔ),使用這種方法劃分的網(wǎng)格和所分析的物質(zhì)結(jié)構(gòu)是相互獨(dú)立的,網(wǎng)格在整個(gè)分析過程中始終保持最初的空間 位置不動(dòng),有限元節(jié)點(diǎn)即為空間點(diǎn),其所在空間的位置在整個(gè)分析過程始終是不變的。很顯然由于算法自身的特點(diǎn),網(wǎng)格的大小形狀和 空間位置不變,因此在整個(gè)數(shù)值模擬過程中,各個(gè)迭代過程中計(jì)算數(shù)值的精度是不變的。但這種方法在物質(zhì)邊界的捕捉上是困難的。多 用于流體的分析中。使用這種方法時(shí)網(wǎng)格與網(wǎng)格之間物質(zhì)是可以流動(dòng)的。ALE方法最初出現(xiàn)于數(shù)值模擬流體動(dòng)力學(xué)問題的有限差分方法中。這種方法兼具Lagrange方法和Euler方法二者的特長(zhǎng),即首
3、先在結(jié)構(gòu) 邊界運(yùn)動(dòng)的處理上它引進(jìn)了 Larange方法的特點(diǎn),因此能夠有效的跟蹤物質(zhì)結(jié)構(gòu)邊界的運(yùn)動(dòng);其次在內(nèi)部網(wǎng)格的劃分上,它吸收了 Euler 的長(zhǎng)處,即是使內(nèi)部網(wǎng)格單元獨(dú)立于物質(zhì)實(shí)體而存在,但它又不完全和Euler網(wǎng)格相同,網(wǎng)格可以根據(jù)定義的參數(shù)在求解過程中適當(dāng)調(diào)整 位置,使得網(wǎng)格不致出現(xiàn)嚴(yán)重的畸變。這種方法在分析大變形問題時(shí)是非常有利的。使用這種方法時(shí)網(wǎng)格與網(wǎng)格之間物質(zhì)也是可以流動(dòng) 的。固體結(jié)構(gòu)分析中一般都選用lagrange坐標(biāo),實(shí)際上lagrange euler法在有限元中體現(xiàn)的節(jié)點(diǎn)意義正如樓主所述,但是本質(zhì)牽扯的是 參考什么樣的坐標(biāo)來描述應(yīng)力應(yīng)變關(guān)系。比如我們通常給的材料應(yīng)力應(yīng)變曲
4、線,在某個(gè)應(yīng)力下對(duì)應(yīng)的應(yīng)變值的大小,究竟是怎么哪個(gè)參 考坐標(biāo)而言的。lagrange坐標(biāo)是參考結(jié)構(gòu)變形前的狀態(tài)來描述應(yīng)力應(yīng)變關(guān)系,而lagrange坐標(biāo)求解時(shí)又分total lagrange和updated lagrange法,前者一直以初始時(shí)的構(gòu)形為參考坐標(biāo),而后者在每個(gè)增量步更新其參考坐標(biāo)為前一步的構(gòu)形。兩者并沒有精度上的差距, 因?yàn)閮烧叩目刂品匠滩灰粯印_x用哪種lagrange法則主要看材料的應(yīng)力應(yīng)變關(guān)系是按什么坐標(biāo)給的,而我們一般的應(yīng)力應(yīng)變關(guān)系都是按 初始時(shí)的構(gòu)形給的。所以一般有限元軟件的結(jié)構(gòu)大變形都采用total lagrange法。而Euler法則是以現(xiàn)時(shí)構(gòu)形(變形后的構(gòu)形)為參考
5、坐 標(biāo)來描述應(yīng)力應(yīng)變關(guān)系的。我們知道,有限元求解的就是變形后的結(jié)構(gòu)構(gòu)形,因此計(jì)算前是根本不可能知道現(xiàn)時(shí)構(gòu)形是什么樣子的,因 此要使用euler法,網(wǎng)格的劃分只能對(duì)整個(gè)空間進(jìn)行劃分,而不可能附著在結(jié)構(gòu)上。這就是樓主說的euler網(wǎng)格的劃分。當(dāng)然對(duì)于我個(gè)人 而言,主要是搞結(jié)構(gòu)的,因此對(duì)lagragne更加了解點(diǎn),也知道里面用的是green應(yīng)變和第二p.k.應(yīng)力(克希霍夫應(yīng)力)。得到第二p.k. 應(yīng)力后通過轉(zhuǎn)換關(guān)系即可得到cauchy應(yīng)力(真實(shí)應(yīng)力)。對(duì)于euler應(yīng)變和almansi應(yīng)力,本人就不是很清楚了。要求版主加分哦人_人abaqus處理大變形采用的更新拉個(gè)朗日方法,不是完全拉個(gè)朗日方法吧?
6、版主,我相信abaqus用的是total lagrange法。因?yàn)橐粋€(gè)簡(jiǎn)單的現(xiàn)實(shí)放著:就是我們的本構(gòu)關(guān)系是根據(jù)初始構(gòu)形給它的。比如說當(dāng)塑性 應(yīng)變達(dá)到0.025時(shí)候,對(duì)應(yīng)的應(yīng)力是500M,那么這個(gè)0.025的塑性應(yīng)變,是根據(jù)什么參考構(gòu)形給的呢?很明顯,是根據(jù)初始構(gòu)形。因?yàn)?真實(shí)應(yīng)變的表達(dá)式是那個(gè)積分,積出來是當(dāng)前長(zhǎng)度比原始長(zhǎng)度的自然對(duì)數(shù)。可以參考abaqus理論手冊(cè):Equilibrium and virtual workProcedures: overview and basic equations4. Mechanical Constitutive Theories用戶手冊(cè):abaqus u
7、sers manul:Conventions既然有了本構(gòu)關(guān)系,不同構(gòu)形之間無非就是通過變形剃度互相推拉而已。abaqus中的應(yīng)力度量為Cauchy應(yīng)力,理論手冊(cè)4. Mechanical Constitutive Theories給出的abaqus的幾乎所有的本構(gòu)關(guān)系形式基本上也都是柯西應(yīng)力的形式。個(gè)人感覺更新拉個(gè)朗日比較容易收斂吧,一個(gè)大變形下,如果保持積分一直在初始構(gòu)形上,變形剃度越來越大,從數(shù)值上來說精度可能 會(huì)降低。這或許是HKS等人選擇更新拉個(gè)朗日的原因吧。首先,在小變形底下,不管克希霍夫應(yīng)力還是歐拉應(yīng)力都退化為柯西應(yīng)力,就無所謂采用什么坐標(biāo)了。而大變形條件下,三者都是 有差異的。歐
8、拉應(yīng)力我不了解,但是大變形下求得克希霍夫應(yīng)力后,可根據(jù)克希霍夫和柯西應(yīng)力的轉(zhuǎn)換關(guān)系求得柯西應(yīng)力。第二,total lagrange和updated lagrange在數(shù)值計(jì)算上使用的公式是不同的,一般書上都講兩者在效率和精度上基本不存在差異,關(guān)鍵 看初始條件是參照什么坐標(biāo)給的。所以,如果保持積分一直在初始構(gòu)形上,變形剃度越來越大,從數(shù)值上來說精度可能會(huì)降低基本不成 立。只能說我認(rèn)為abaqus用的是total lagrange法,樓上的不同意見咱以后再交流。因?yàn)橄聛硪欢螘r(shí)間又要用重新用ansys開始做課題了。媽媽的,去年的某個(gè)項(xiàng)目型號(hào)一改型,有限元得重新建模,工程上的這些東西,體力活,但是相當(dāng)
9、的麻煩! !我一直對(duì)ABAQUS究竟用什么構(gòu)型來描述問題一直不明白,而前面貼子貌似比較支持ABAQUS默認(rèn)的是total lagrangian描述的觀 點(diǎn)。其實(shí)小變形時(shí)是哪種構(gòu)型基本區(qū)別不大,問題是大變形問題時(shí),對(duì)用不同構(gòu)型算出來的東西還是有區(qū)別的。如果ABAQUS默認(rèn)是 total lagrangian的話,那想用Eulerian描述問題,應(yīng)該怎樣設(shè)置Euler坐標(biāo)其坐標(biāo)系本身是固定得,僅物體運(yùn)動(dòng);Lagrange坐標(biāo)其坐標(biāo)系是放在所描述得物體上隨著物體一起運(yùn)動(dòng)得。例如流體力學(xué)中 拉格朗日坐標(biāo)系,觀察者位于一個(gè)流體單元上,并隨流體一起運(yùn)動(dòng)。歐拉坐標(biāo)系,觀察者位于空間的一個(gè)固定點(diǎn),觀察流過你所
10、在的體 積單元。歐拉坐標(biāo)是指空間坐標(biāo),如果對(duì)于質(zhì)點(diǎn)運(yùn)動(dòng)來說,研究的是不同的質(zhì)點(diǎn)經(jīng)過空間一定點(diǎn)的狀態(tài)拉格朗日坐標(biāo)指的是材料坐標(biāo),如果對(duì)于質(zhì)點(diǎn)運(yùn)動(dòng)來說,是跟隨著質(zhì)點(diǎn)研究質(zhì)點(diǎn)的運(yùn)動(dòng)狀態(tài)簡(jiǎn)單的說就是:歐拉坐標(biāo)固定在空間,拉格朗日坐標(biāo)固定在材料上歐拉坐標(biāo)指一點(diǎn)在空間的位置,拉格朗日坐標(biāo)標(biāo)記一個(gè)材料點(diǎn)ALE、Lagrange、Euler是數(shù)值模擬中處理連續(xù)體的廣泛應(yīng)用的三種方法Lagrange方法多用于固體結(jié)構(gòu)的應(yīng)力應(yīng)變分析,這種方法以物質(zhì)坐標(biāo)為基礎(chǔ), 其所描述的網(wǎng)格單元將以類似“雕刻”的方式劃分在用于分析的結(jié)構(gòu)上,即是說 采用Lagrange方法描述的網(wǎng)格和分析的結(jié)構(gòu)是一體的,有限元節(jié)點(diǎn)即為物質(zhì)點(diǎn)。 采
11、用這種方法時(shí),分析結(jié)構(gòu)的形狀的變化和有限單元網(wǎng)格的變化完全是一致的(因?yàn)橛邢拊?jié)點(diǎn)就為物質(zhì)點(diǎn)),物質(zhì)不會(huì)在單元與單元之間發(fā)生流動(dòng)。這種方 法主要的優(yōu)點(diǎn)是能夠非常精確的描述結(jié)構(gòu)邊界的運(yùn)動(dòng),但當(dāng)處理大變形問題時(shí), 由于算法本身特點(diǎn)的限制,將會(huì)出現(xiàn)嚴(yán)重的網(wǎng)格畸變現(xiàn)象,因此不利于計(jì)算的進(jìn) 行。Euler方法以空間坐標(biāo)為基礎(chǔ),使用這種方法劃分的網(wǎng)格和所分析的物質(zhì)結(jié)構(gòu)是 相互獨(dú)立的,網(wǎng)格在整個(gè)分析過程中始終保持最初的空間位置不動(dòng),有限元節(jié)點(diǎn) 即為空間點(diǎn),其所在空間的位置在整個(gè)分析過程始終是不變的。很顯然由于算法 自身的特點(diǎn),網(wǎng)格的大小形狀和空間位置不變,因此在整個(gè)數(shù)值模擬過程中,各 個(gè)迭代過程中計(jì)算數(shù)值的精度是不變的。但這種方法在物質(zhì)邊界的捕捉上是困難 的。多用于流體的分析中。使用這種方法時(shí)網(wǎng)格與網(wǎng)格之間物質(zhì)是可以流動(dòng)的。ALE方法最初出現(xiàn)于數(shù)值模擬流體動(dòng)力學(xué)問題的有限差分方法中。這種方法兼具 Lagrange方法和Euler方法二者的特長(zhǎng),即首先在結(jié)構(gòu)邊界運(yùn)動(dòng)的處理上它引 進(jìn)了
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 數(shù)據(jù)驅(qū)動(dòng)的軟件創(chuàng)新機(jī)制與產(chǎn)業(yè)升級(jí)路徑研究
- CUDA并行編程從入門到實(shí)戰(zhàn)指南
- 體育康復(fù)課程體系創(chuàng)新設(shè)計(jì)與實(shí)踐探索
- 施工現(xiàn)場(chǎng)安全風(fēng)險(xiǎn)防控與整改指南
- 跨境數(shù)據(jù)傳輸合規(guī)-洞察及研究
- 養(yǎng)老院消防安全隱患排查表
- 兼職律師執(zhí)業(yè)管理辦法
- 生產(chǎn)經(jīng)營(yíng)發(fā)生安全事故后應(yīng)立即報(bào)告
- 醫(yī)療專家審核管理辦法
- 科技創(chuàng)新對(duì)戰(zhàn)略性新興產(chǎn)業(yè)發(fā)展競(jìng)爭(zhēng)力的影響與作用機(jī)制
- 2023年公路工程監(jiān)理工作合同管理內(nèi)容及有效措施
- 一年級(jí)看圖寫話(教學(xué))課件
- 嚴(yán)重藥物不良反應(yīng)診斷與處理
- 直流屏原理-課件
- 加藥設(shè)備安裝 檢驗(yàn)批施工質(zhì)量驗(yàn)收表
- 崗位技能評(píng)定機(jī)考考場(chǎng)規(guī)則
- 盡職調(diào)查所用相關(guān)表格(全)
- 三基-學(xué)校兒童少年衛(wèi)生學(xué)(200題)練習(xí)
- 老年康養(yǎng)服務(wù)中心項(xiàng)目可行性研究報(bào)告寫作參考范文
- 生物質(zhì)中纖維素、半纖維素和木質(zhì)素含量的測(cè)定
- 枸杞采摘合同
評(píng)論
0/150
提交評(píng)論