| 5 | 1/1 | 返回列表 |
| 查看: 1863 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
duchao銅蟲 (小有名氣)
|
[求助]
有限元法的理論求助
|
||
我最近一直在學(xué)習(xí)使用comsol軟件模擬電磁場傳輸?shù)膯栴},comsol軟件是以有限元法(FEM)為基礎(chǔ)進(jìn)行數(shù)值計(jì)算的。但是對于其中有限元模擬的理論原理我研究半天也沒看明白,從Maxwell方程組到各個(gè)單元的那些方程,到最后推倒出一個(gè)矩陣,底是怎么聯(lián)系推倒出來的呢?希望得到高手幫助指點(diǎn),或者是能提供一些對應(yīng)的有限元書籍或文獻(xiàn),總之是能讓我明白明白怎么回事就行,在下不勝感激!![]() ![]() |
木蟲 (著名寫手)
|
這個(gè)問題呢 呵呵 確實(shí) 要看你要從那個(gè)角度上去理解 如果你懂點(diǎn)python,不如玩一玩這個(gè)項(xiàng)目http://www.fenicsproject.org/ 個(gè)人喜歡從最簡單的東西去理解,譬如給你一個(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í)你會希望g能有個(gè)性質(zhì)就是在邊界上為0),然后你再把g表示成basis function的疊加(邊界上為零的條件自然成立)。此時(shí)你再把f(x,y)表示成basis function的疊加(設(shè)系數(shù)為矢量U),帶入到弱形式里,你會得到stiffness matrix K,而Laplace(f(x))這一端呢,你會得到一個(gè)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方程,那求解上面的線性方程就行了。 個(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ā)表 | |
|---|---|---|---|---|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑,不限專業(yè),有學(xué)上就行 +5 | lllllos 2026-03-30 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +5 | 蘿卜燉青菜 2026-03-28 | 6/300 |
|
|
[考研] 322求調(diào)劑 +8 | 宋明欣 2026-03-27 | 8/400 |
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +13 | 丹青奶蓋 2026-03-26 | 15/750 |
|
|
[考研] 337求調(diào)劑 +6 | 《樹》 2026-03-29 | 6/300 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
[考研] 296求調(diào)劑 +4 | 汪。! 2026-03-25 | 7/350 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|