欢迎在此讨论“时域有限元”: FETD Formula by Guoqiang_He.doc

本帖最后由 tingo 于 2011-11-20 08:35 编辑

最近在做时域有限元算法，其中遇到一些问题，同时也解决了一些问题，感觉国内目前做这个方向的人比较少，大家可以在此讨论，以期共同学习进步。
跟帖给出了自己推到的基于电场波动方程的Galerkin公式。

自己先顶起！

友情帮顶！支持，希望越来越好，要做好，还真的特别规划下

你计算的什么问题

由于刚涉入矢量有限元，目前计算的是简单的二维问题，平行板波导。

之前在导体壁上强加第一类边界条件，仿真的TE模式，电场波形结果是错的，后面删去后波形正确了

矢量有限元金建铭的书上有一部分介绍，我没做过这方面，不过感觉矢量元是一个发展方向，有助于解决复杂问题。
也许哪一天心血来潮，会做一下，目前没有打算，原因一个是选择边界条件方面，没有确定是使用吸收边界还是边界积分，再一个，这个问题比标量问题复杂的多，一时难以得到一个预期计算结果，容易让人放弃。

2D的TE、TM模式我做过RCS计算，金属和介质方面，不过不是用有限元。

的确，矢量有限元相对于基于结点有限元复杂的多，一时很难出正确的结果。不过，既然我已经选择做这方面的研究了，还是想着深入的做一做。

那能不能谈谈用什么算法做的呢，相比于有限元的优势在哪

在此贴两张自己做的FETD仿真。
仿真结果
=*******************************************************=
**        USFETD仿真求解X波段标准波导BJ-100主模式         **
=*******************************************************=
TE10模截止频率理论值 : 6.56168 GHz
TE10模截止频率仿真值 : 6.54254 GHz

这个是300步时间步进计算出的结果，误差稍微大一些！

分享自己推导的二维的基于电场波动方程的时域有限元公式。其中难免会有错误和不当之处，如有发现，恳请及时在此跟帖说明或者以邮件说明，谢谢！