BEM的matlab实现！高手能看看我的程序问题吗！

dxw007 · 发表于 2007-5-17 17:09:09

BEM的matlab实现！高手能看看我的程序问题吗！

BEM的matlab实现！高手能看看我的程序问题吗！
我使用的是常数单元法来解决位势问题！
谢谢大家！源程序在bem.rar中

foxin612 · 发表于 2007-5-17 17:59:56

不错呀,不过不懂MATLAB,你看过用FORTRAN遍的常单元边界元程序吗

dxw007 · 发表于 2007-5-21 22:27:30

我看过就是基于FORTRAN的思想可是我用matlab实现以后就有问题了！
不知道哪位高人能否指点一二！谢谢了！

happysimwe · 发表于 2007-6-21 13:01:56

边界元的格式很多，位移，应力，对称等等，得给一些相应的理论公式才容易看看。

dxw007 · 发表于 2007-6-22 09:08:07

楼上你好！我在写程序时就是根据边界元公式：HU=GQ
因为我做的是电磁问题所以上面公式中的U表示电位（已知），Q表示电位的法向导数值（未知），而H和G就是影响系数我想所有有关边界元原理问题中关于Ｈ和Ｇ的求解都是一样的吧！还不知道能否赐教！谢谢！

shigang · 发表于 2007-6-22 10:16:53

你这个问题讲的不太清楚，我感觉划分单元时，将外导体（环状）内外表面划分单元节点，再进行求解。你的是不是对两个导体分别进行求解，但是，对与外导体内表面和内导体的外表面之间的耦合条件如何分析？若是他们之间的电压相同，那么，直接取外导体进行分析即可，当要在其内外表面都要划分单元。

happysimwe · 发表于 2007-6-22 13:31:52

这是普通奇异型边界元程序。可能是哪个变量混了，刚好在1的时候相等。简单的方法有两个：
1）对到H G后代入一些已知解，看看哪些方程有问题
2）先求积分 1 与 0.4 情况，记录H G, 在积分2与0.8, 对比两者系数异常
找到出问题的行之后，减少单元循环，找出对应的单元，到单元级后减少并对比积分点，逐步缩小范围。

dxw007 · 发表于 2007-6-26 21:20:31

我按照大家的指点再找找看！

dxw007 · 发表于 2007-6-26 21:27:54

你好！因为对于此类问题，内外边界是独立的不存在耦合问题，我们已知内边界和外边界（内外边界都为圆）的电压，现在就是要分别按照HU=GQ求出各自边界上的Q值，进而来求解别的物理量，所以核心就是求解Q,剩余的问题就是我帖子上你下载到的内容！还望指点！谢谢！

yangfl02 · 发表于 2008-8-11 23:21:47

请问你的问题解决了吗?

visualheart · 发表于 2010-4-17 11:53:20

一个初学者的谢意，请笑纳

uacademic · 发表于 2010-4-23 06:50:23

有同样的问题。

账号		自动登录	找回密码
密码			注册

BEM的matlab实现！高手能看看我的程序问题吗！

本帖子中包含更多资源

评分

相关帖子

看过

回复 #1 dxw007 的帖子

简单公式

回复 #5 dxw007 的帖子

谢谢

时刚！

浏览过的版块