| 5 | 1/1 | 返回列表 |
| 查看: 1867 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
duchao銅蟲 (小有名氣)
|
[求助]
有限元法的理論求助
|
|
我最近一直在學(xué)習(xí)使用comsol軟件模擬電磁場(chǎng)傳輸?shù)膯栴},comsol軟件是以有限元法(FEM)為基礎(chǔ)進(jìn)行數(shù)值計(jì)算的。但是對(duì)于其中有限元模擬的理論原理我研究半天也沒看明白,從Maxwell方程組到各個(gè)單元的那些方程,到最后推倒出一個(gè)矩陣,底是怎么聯(lián)系推倒出來的呢?希望得到高手幫助指點(diǎn),或者是能提供一些對(duì)應(yīng)的有限元書籍或文獻(xiàn),總之是能讓我明白明白怎么回事就行,在下不勝感激!![]() ![]() |
銅蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
木蟲 (著名寫手)
|
這個(gè)問題呢 呵呵 確實(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))這種基本問題,如果這個(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 |
銅蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298求調(diào)劑 +3 | zzz,,r 2026-04-02 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +7 | pipiver 2026-03-30 | 7/350 |
|
|
[考研] 材料調(diào)劑 +14 | 一樣YWY 2026-04-01 | 14/700 |
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +10 | 腳滑的守法公民 2026-03-27 | 10/500 |
|
|
[考研] 330分求調(diào)劑 +11 | qzenlc 2026-03-29 | 11/550 |
|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600,349 +9 | 吃的不少 2026-03-30 | 9/450 |
|
|
[考研] 291求調(diào)劑 +3 | 迷蒙木木 2026-04-01 | 4/200 |
|
|
[考研] 生物學(xué)296求調(diào)劑 +10 | 湯圓包 2026-03-29 | 14/700 |
|
|
[考研] 一志愿華南師范361分,化學(xué)求調(diào)劑 +4 | Nicole88888 2026-04-01 | 4/200 |
|
|
[考研] 358求調(diào)劑 +3 | 王向陽花 2026-03-31 | 3/150 |
|
|
[考研] 一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +10 | @taotao 2026-03-31 | 11/550 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +11 | 蘿卜燉青菜 2026-03-28 | 12/600 |
|
|
[考研] 調(diào)劑 +4 | GK72 2026-03-30 | 4/200 |
|
|
[考研] 343求調(diào)劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 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 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|