BEM的matlab实现！高手能看看我的程序问题吗！(页 1) - A50：边界元理论及软件应用 - SimWe仿真论坛 |Simwe.com|仿真|设计|有限元|虚拟仪器

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

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

BEM的matlab实现！高手能看看我的程序问题吗！%`4[c*bi:I5]T
9df2O:d[:k{ x,j
BEM的matlab实现！高手能看看我的程序问题吗！
我使用的是常数单元法来解决位势问题！
谢谢大家！源程序在bem.rar中 [attach]122604[/attach]

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

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

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

看过

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

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

回复 #1 dxw007 的帖子

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

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

简单公式

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

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

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

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

回复 #5 dxw007 的帖子

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

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

谢谢

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

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

时刚！

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

nibbomm 发表于 2008-4-28 01:34

有点意思

话是如此了，看着还是有点意思的，不过究竟是不是这样，还需要大家仔细想想看
9{$`8\ F(X2uMR7z
GgU\2`Q

P%f |mq
----------------------------------------------#b-Q8i&Cv0s'o
[url=http://www.eve-online-isk.cn]eve online isk[/url] and [url=http://www.wowgold.ws]wow gold[/url] buy [url=http://www.wowpowerleveling-gold.com]wow power leveling[/url] and [url=http://www.dofusbuy.fr]dofus kamas[/url] cheap [url=http://www.rs2lostcity.com]runescape money[/url]

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

请问你的问题解决了吗?

页: [1]

SimWe仿真论坛's Archiver