| 查看: 1300 | 回复: 6 | ||
| 【悬赏金币】回答本帖问题，作者duchao将赠送您 8 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
有限元法的理论求助
|
||
我最近一直在学习使用comsol软件模拟电磁场传输的问题，comsol软件是以有限元法（FEM）为基础进行数值计算的。但是对于其中有限元模拟的理论原理我研究半天也没看明白，从Maxwell方程组到各个单元的那些方程，到最后推倒出一个矩阵，底是怎么联系推倒出来的呢？希望得到高手帮助指点，或者是能提供一些对应的有限元书籍或文献，总之是能让我明白明白怎么回事就行，在下不胜感激！！ |
回复此楼» 猜你喜欢
Error in routine lambda (100): wrong # or too many modes
已经有0人回复
-
祈福国自然面上
已经有61人回复
-
物理学I论文润色/翻译怎么收费?
已经有272人回复
-
请问castep计算出的能带图横轴数据缩放因子怎么计算
已经有1人回复
-
Result in Physics 从见刊开始 一般多久可以被SCI收录
已经有26人回复
-
面青地会评时间？？？
已经有8人回复
-
科研求助
已经有1人回复
-
深圳大学高等研究院暑期夏令营
已经有10人回复
-
深圳大学高等研究院暑期夏令营
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有限单元法原理与应用+第2版(朱伯芳)
已经有10人回复
- 求助几个理论问题
已经有53人回复
- 求助：中了2011理论物理专款博士启动项目的（1年小额项目）-能否申请2012年的国青？
已经有4人回复
- 【求助】编织物衬垫有限元建模弹性性能分析
已经有5人回复
- 【求助】管理理论包括哪些具体理论？
已经有12人回复
- 【求助】理论解释
已经有13人回复
- 【求助】初学有限元，弱弱问一句，如果已知边界都是位移边界条件，怎么求解？
已经有11人回复
- 【求助】第一性原理与密度泛函理论的关系是什么
已经有15人回复
- 求助：怎样增强论文的理论性？
已经有5人回复
- 100金币下载论文求助
已经有3人回复
- 【求助】请教有限元方面的问题？
已经有4人回复
- 求助 搞岩石的研究大牛？注意不是土，散金币
已经有13人回复
- 【求助】钛酸锂理论容量问题
已经有7人回复
4楼2012-08-05 11:09:13
q68
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 贵宾: 1.615
- 金币: 10176.2
- 散金: 10
- 红花: 3
- 帖子: 1717
- 在线: 148.6小时
- 虫号: 325340
- 注册: 2007-03-17
- 性别: GG
- 专业: 金属材料的合金相、相变及
duchao2楼2012-08-04 17:00:49
wangww2011
木虫 (著名写手)
- 应助: 11 (小学生)
- 金币: 4023.1
- 散金: 2709
- 红花: 18
- 沙发: 1
- 帖子: 1915
- 在线: 1537.1小时
- 虫号: 772953
- 注册: 2009-05-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+2, 多谢应助 2012-08-04 21:22:52
duchao: 回帖置顶 2012-08-05 11:07:32
duchao: 金币+2, ★★★很有帮助 2012-08-05 11:08:14
感谢参与，应助指数 +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
5楼2012-08-05 11:09:39
10