| 24小時(shí)熱門(mén)版塊排行榜 |
| 7 | 1/1 | 返回列表 |
| 查看: 1859 | 回復(fù): 6 | ||
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) (著名寫(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) (著名寫(xiě)手)
銅蟲(chóng) (小有名氣)
銅蟲(chóng) (小有名氣)
送鮮花一朵
|
鐵蟲(chóng) (小有名氣)

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京工業(yè)大學(xué),324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工(085600)296求調(diào)劑 +17 | 稻妻小編 2026-03-26 | 17/850 |
|
|
[考研] 322求調(diào)劑:一志愿湖南大學(xué) 材料與化工(085600),已過(guò)六級(jí)。 +4 | XX小鄧 2026-03-29 | 4/200 |
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
[考研] 一志愿華北電力大學(xué)能動(dòng)專(zhuān)碩,293,求調(diào)劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 復(fù)試調(diào)劑,一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專(zhuān)業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問(wèn)津 2026-03-23 | 5/250 |
|