| 7 | 1/1 | 返回列表 |
| 查看: 1862 | 回復(fù): 6 | ||
duchao銅蟲 (小有名氣)
|
[求助]
有限元法的理論求助
|
我最近一直在學(xué)習(xí)使用comsol軟件模擬電磁場傳輸?shù)膯栴},comsol軟件是以有限元法(FEM)為基礎(chǔ)進(jìn)行數(shù)值計(jì)算的。但是對于其中有限元模擬的理論原理我研究半天也沒看明白,從Maxwell方程組到各個單元的那些方程,到最后推倒出一個矩陣,底是怎么聯(lián)系推倒出來的呢?希望得到高手幫助指點(diǎn),或者是能提供一些對應(yīng)的有限元書籍或文獻(xiàn),總之是能讓我明白明白怎么回事就行,在下不勝感激!![]() ![]() |
鐵桿木蟲 (著名寫手)
木蟲 (著名寫手)
|
這個問題呢 呵呵 確實(shí) 要看你要從那個角度上去理解 如果你懂點(diǎn)python,不如玩一玩這個項(xiàng)目http://www.fenicsproject.org/ 個人喜歡從最簡單的東西去理解,譬如給你一個函數(shù)f(x,y),然后給你一個mesh,如何求Laplace(f(x))這種基本問題,如果這個mesh是規(guī)則的,那直接根據(jù)經(jīng)驗(yàn),就可以有限差分得出一個矩陣乘以f(x,y)。如果是三角形這種不規(guī)則的,就需要點(diǎn)點(diǎn)數(shù)學(xué)了,基本思路就是乘以一個任意函數(shù)g,然后用Divergence theorem轉(zhuǎn)成弱形式(此時你會希望g能有個性質(zhì)就是在邊界上為0),然后你再把g表示成basis function的疊加(邊界上為零的條件自然成立)。此時你再把f(x,y)表示成basis function的疊加(設(shè)系數(shù)為矢量U),帶入到弱形式里,你會得到stiffness matrix K,而Laplace(f(x))這一端呢,你會得到一個load vector L,然后你會發(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方程,那求解上面的線性方程就行了。 個人認(rèn)為這就是有限元的基本原理。 ps:關(guān)于basis function你可以假設(shè)在一個unit triangle上已知三個頂點(diǎn)的函數(shù)值,然后去插值,然后任意三角形直接做一個affine變換就行了。 這是fenics的一個教程, http://www.fenicsproject.org/doc ... he-poisson-equation |
銅蟲 (小有名氣)
銅蟲 (小有名氣)
送鮮花一朵
|
鐵蟲 (小有名氣)

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 367求調(diào)劑 +3 | 芋泥啵啵… 2026-03-28 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +25 | RichLi_ 2026-03-25 | 25/1250 |
|
|
[考研] 318求調(diào)劑 +7 | 陳晨79 2026-03-30 | 7/350 |
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +6 | 深V宿舍吧 2026-03-29 | 6/300 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +9 | Koxui 2026-03-28 | 9/450 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +7 | 26要上岸 2026-03-28 | 7/350 |
|
|
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4 | 幸運(yùn)的三木 2026-03-29 | 5/250 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我。! +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 359求調(diào)劑 +4 | 王了個楠 2026-03-25 | 4/200 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研]
|
黃粱一夢千年 2026-03-24 | 3/150 |
|