24小时热门版块排行榜    

查看: 2040  |  回复: 6

duchao

铜虫 (小有名气)

[求助] 有限元法的理论求助

我最近一直在学习使用comsol软件模拟电磁场传输的问题,comsol软件是以有限元法(FEM)为基础进行数值计算的。但是对于其中有限元模拟的理论原理我研究半天也没看明白,从Maxwell方程组到各个单元的那些方程,到最后推倒出一个矩阵,底是怎么联系推倒出来的呢?希望得到高手帮助指点,或者是能提供一些对应的有限元书籍或文献,总之是能让我明白明白怎么回事就行,在下不胜感激!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

wangww2011

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
csgt0: 金币+2, 多谢应助 2012-08-04 21:22:52
duchao: 回帖置顶 2012-08-05 11:07:32
duchao: 金币+2, ★★★很有帮助 2012-08-05 11:08:14
这个问题呢 呵呵 确实 要看你要从那个角度上去理解
如果你懂点python,不如玩一玩这个项目http://www.fenicsproject.org/
个人喜欢从最简单的东西去理解,譬如给你一个函数f(x,y),然后给你一个mesh,如何求Laplace(f(x))这种基本问题,如果这个mesh是规则的,那直接根据经验,就可以有限差分得出一个矩阵乘以f(x,y)。如果是三角形这种不规则的,就需要点点数学了,基本思路就是乘以一个任意函数g,然后用Divergence theorem转成弱形式(此时你会希望g能有个性质就是在边界上为0),然后你再把g表示成basis function的叠加(边界上为零的条件自然成立)。此时你再把f(x,y)表示成basis function的叠加(设系数为矢量U),带入到弱形式里,你会得到stiffness matrix K,而Laplace(f(x))这一端呢,你会得到一个load vector L,然后你会发现有
[code]
L Laplace(f(x)) = K U
Laplace(f(x)) = L^-1 K U
[\code]
注意这里是假设f(x,y)已知的情况,求Laplace(f(x)) ,如果是f(x,y)未知,Laplace(f(x))已知,譬如Poisson方程,那求解上面的线性方程就行了。
个人认为这就是有限元的基本原理。
ps:关于basis function你可以假设在一个unit triangle上已知三个顶点的函数值,然后去插值,然后任意三角形直接做一个affine变换就行了。
这是fenics的一个教程,
http://www.fenicsproject.org/doc ... he-poisson-equation
3楼2012-08-04 18:23:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

q68

铁杆木虫 (著名写手)


csgt0: 金币+1, 多谢建议 2012-08-04 21:23:03
嘿嘿 这个东西不是 一般人能说的清楚的。
也不是 一两句能说的清楚的。
建议先不要管软件问题,先从书籍看看理论再说。

» 本帖已获得的红花(最新10朵)

2楼2012-08-04 17:00:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duchao

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by wangww2011 at 2012-08-04 18:23:54
这个问题呢 呵呵 确实 要看你要从那个角度上去理解
如果你懂点python,不如玩一玩这个项目http://www.fenicsproject.org/
个人喜欢从最简单的东西去理解,譬如给你一个函数f(x,y),然后给你一个mesh,如何求Laplac ...

我在好好看看,非常感谢!
4楼2012-08-05 11:09:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duchao

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by q68 at 2012-08-04 17:00:49
嘿嘿 这个东西不是 一般人能说的清楚的。
也不是 一两句能说的清楚的。
建议先不要管软件问题,先从书籍看看理论再说。

谢谢
5楼2012-08-05 11:09:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duchao

铜虫 (小有名气)

送鲜花一朵
引用回帖:
2楼: Originally posted by q68 at 2012-08-04 17:00:49
嘿嘿 这个东西不是 一般人能说的清楚的。
也不是 一两句能说的清楚的。
建议先不要管软件问题,先从书籍看看理论再说。

6楼2012-08-05 11:13:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

atu_zhong

铁虫 (小有名气)

我也准备学习这个软件,不知道难不难学呀。唉
好好干
7楼2012-08-15 16:49:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 duchao 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 纯娱乐,不喜欢勿喷 +5 Tide man 2026-08-04 8/400 2026-08-04 11:43 by everxc
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 6/300 2026-08-04 11:34 by MWRYsWbJw4DE
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 4/200 2026-08-04 11:26 by MWRYsWbJw4DE
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
信息提示
请填处理意见