| 查看: 1896 | 回复: 6 | ||
[求助]
有限元法的理论求助
|
我最近一直在学习使用comsol软件模拟电磁场传输的问题，comsol软件是以有限元法（FEM）为基础进行数值计算的。但是对于其中有限元模拟的理论原理我研究半天也没看明白，从Maxwell方程组到各个单元的那些方程，到最后推倒出一个矩阵，底是怎么联系推倒出来的呢？希望得到高手帮助指点，或者是能提供一些对应的有限元书籍或文献，总之是能让我明白明白怎么回事就行，在下不胜感激！！ |
回复此楼» 猜你喜欢
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有98人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有限单元法原理与应用+第2版(朱伯芳)
已经有10人回复
- 求助几个理论问题
已经有53人回复
- 求助：中了2011理论物理专款博士启动项目的（1年小额项目）-能否申请2012年的国青？
已经有4人回复
- 【求助】管理理论包括哪些具体理论？
已经有12人回复
- 【求助】理论解释
已经有13人回复
- 【求助】初学有限元，弱弱问一句，如果已知边界都是位移边界条件，怎么求解？
已经有11人回复
- 【求助】第一性原理与密度泛函理论的关系是什么
已经有15人回复
- 求助：怎样增强论文的理论性？
已经有5人回复
- 100金币下载论文求助
已经有3人回复
- 【求助】请教有限元方面的问题？
已经有4人回复
- 求助 搞岩石的研究大牛？注意不是土，散金币
已经有13人回复
- 【求助】钛酸锂理论容量问题
已经有7人回复
q68
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 贵宾: 1.615
- 金币: 14919.2
- 散金: 10
- 红花: 3
- 帖子: 2340
- 在线: 149.8小时
- 虫号: 325340
- 注册: 2007-03-17
- 性别: GG
- 专业: 金属材料的合金相、相变及
2楼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
4楼2012-08-05 11:09:13
5楼2012-08-05 11:09:39
送鲜花一朵 |
|
6楼2012-08-05 11:13:11
7楼2012-08-15 16:49:43
20