Fortran编非线性有限元的程序，希望交流一下！

lishhit

自己用Fortran编了一个非线性有限元的程序，问题是后来有了新想法，就得改了又改，甚至从头改，现在感觉是程序的数据结构(就是整个程序的编程框架)没有处理好。所以想和用F编非线性有限元程序的大侠们交流一下，学学经验！

若觉得自己感兴趣可以给我发邮件（论坛短消息），单独交流；也可在这里直接回帖，各抒己见也行。

虽然针对非线性有限元，但我想不管什么什么问题，用F编-复杂-的程序应该有一个固定的模式，用这个模式比较好组织数据结构，如果哪位有比较好的F编的程序（指数据结构比较好）也希望能够共享，让大家学习一下！


谢谢！

hillyuan

There are some books considered the programming of FEM code. Pls refer to following link
http://ohio.ikp.liu.se/fe/prog.html
And, you can find a few FEM open source, such as ELmer, FEAppv, GeoFEM written by FORTRAN.

However, I am wondering you can find a satisfactory answer. You will find all those books gives different code design, and, few considered the arts of old-style structural programm strategy currently. Maybe you can find a book of SOFTWRAE ENGINNERRING published earlier than 80s of last century to get some help.

On the other hand, if you are considering object-oriented stretegy, you can find lots of research papers and open source FEM codes to refer to.

lishhit

Thanks, hillyuan. Some source codes you metioned I had read ago. Just as you said:  they give different code design. Another problem is the source codes are so big too understand for a short time(I mean I can't catch it's data structure clearly).

I think someone who had program some FEA codes may have good experiences. So, I hope share these experiences directly based on my knowledge which already had and the very limited time. May be I'm a little lazy. But, I realy need some "visual" help or advise and I think it is efficiency in this way.

[ 本帖最后由 lishhit 于 2007-8-19 23:17 编辑 ]

mellow

我是新手，刚刚在用有限元方法做形状记忆合金中磁感应强度的计算（求解LLG方程，是一个很复杂的非线性方程），也用的是FORTRAN，不过现在还没有什么头绪，希望lishhit你能给些指导好吗？
yxl_xl_l@hotmail.com

lxy_99

怎么什么都美誉

songseunghun

链接都打不开了呀，不知道怎么回事

army

很不错，顶！！

xiaoqiushui318

楼主能否介绍一下你的程序一开始是啥结构，后来为啥要改？
让我们这些后来者也学习一下经验，避免走弯路。

cfem

自己用Fortran编了一个非线性有限元的程序，问题是后来有了新想法，就得改了又改，甚至从头改，现在感觉是程序的数据结构(就是整个程序的编程框架)没有处理好。所以想和用F编非线性有限元程序的大侠们交流一下，学学 ...
lishhit 发表于 2007-8-19 16:32

一不小心，竟然发现这是07年的帖子 楼上的几位兄弟给倒腾出来了：）

自己写小的程序，只要是面向对象或者模块化的程度不高，改来改去，其实是必然的  当初动手些的时候就没考虑通用性  通用软件的难处也体现在这个地方  虽然大家讲到了最后就是个解方程组的技术，可实际上远不止于此  要视同志们还记得 对于显式算法的话现如今的做法是没有形成总刚这一步骤的 但是这样的话 也有弊病 说来话长……

某种程度上这也是个经验的认识，大家做的东西如果太专反而不利于看到程序的整体 从一个小领域到另一个专的领域 改动就是几乎是不可避免的了

提升一下高度，可见啊，合作和交流是必须的：）

snoweak

通用的程序未必存在 其实可以借鉴Openfoam的模块化思想

当你改动的时候 那些不变得 比如向量旋转 甚至张量等 都可以变成模块 下次用的时候再用

没必要上来就编通用的程序 solver 也没有具体的格式



9# cfem