| 5 | 1/1 | 返回列表 |
| 查看: 1864 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
duchao銅蟲(chóng) (小有名氣)
|
[求助]
有限元法的理論求助
|
|
我最近一直在學(xué)習(xí)使用comsol軟件模擬電磁場(chǎng)傳輸?shù)膯?wèn)題,comsol軟件是以有限元法(FEM)為基礎(chǔ)進(jìn)行數(shù)值計(jì)算的。但是對(duì)于其中有限元模擬的理論原理我研究半天也沒(méi)看明白,從Maxwell方程組到各個(gè)單元的那些方程,到最后推倒出一個(gè)矩陣,底是怎么聯(lián)系推倒出來(lái)的呢?希望得到高手幫助指點(diǎn),或者是能提供一些對(duì)應(yīng)的有限元書(shū)籍或文獻(xiàn),總之是能讓我明白明白怎么回事就行,在下不勝感激!![]() ![]() |
銅蟲(chóng) (小有名氣)
鐵桿木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)
|
這個(gè)問(wèn)題呢 呵呵 確實(shí) 要看你要從那個(gè)角度上去理解 如果你懂點(diǎn)python,不如玩一玩這個(gè)項(xiàng)目http://www.fenicsproject.org/ 個(gè)人喜歡從最簡(jiǎn)單的東西去理解,譬如給你一個(gè)函數(shù)f(x,y),然后給你一個(gè)mesh,如何求Laplace(f(x))這種基本問(wèn)題,如果這個(gè)mesh是規(guī)則的,那直接根據(jù)經(jīng)驗(yàn),就可以有限差分得出一個(gè)矩陣乘以f(x,y)。如果是三角形這種不規(guī)則的,就需要點(diǎn)點(diǎn)數(shù)學(xué)了,基本思路就是乘以一個(gè)任意函數(shù)g,然后用Divergence theorem轉(zhuǎn)成弱形式(此時(shí)你會(huì)希望g能有個(gè)性質(zhì)就是在邊界上為0),然后你再把g表示成basis function的疊加(邊界上為零的條件自然成立)。此時(shí)你再把f(x,y)表示成basis function的疊加(設(shè)系數(shù)為矢量U),帶入到弱形式里,你會(huì)得到stiffness matrix K,而Laplace(f(x))這一端呢,你會(huì)得到一個(gè)load vector L,然后你會(huì)發(fā)現(xiàn)有 [code] L Laplace(f(x)) = K U Laplace(f(x)) = L^-1 K U [\code] 注意這里是假設(shè)f(x,y)已知的情況,求Laplace(f(x)) ,如果是f(x,y)未知,Laplace(f(x))已知,譬如Poisson方程,那求解上面的線性方程就行了。 個(gè)人認(rèn)為這就是有限元的基本原理。 ps:關(guān)于basis function你可以假設(shè)在一個(gè)unit triangle上已知三個(gè)頂點(diǎn)的函數(shù)值,然后去插值,然后任意三角形直接做一個(gè)affine變換就行了。 這是fenics的一個(gè)教程, http://www.fenicsproject.org/doc ... he-poisson-equation |
銅蟲(chóng) (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 295分求調(diào)劑 +7 | W55j 2026-03-30 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)求調(diào)劑 +6 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
[考研] 材料與化工272求調(diào)劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
[考博] 26申博自薦 +6 | whh869393 2026-03-24 | 6/300 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 311(085601)求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|