开发带有一定前后处理功能的有限元“玩具"软件设想

liuichini

解NS方程倒不在于编码的工作量，算法搞定了，编码倒不是什么大事。

hillyuan

caoer 发表于 2012-6-13 00:05
如果求解ns方程的编码量是1的话，
前处理的编码量是100.
后处理看需求，介于1-10之间。

使用可靠的并行middleware(如Petsc), 并行的调试用不了那么多少时间的。

hillyuan

liuichini 发表于 2012-6-13 00:52
解NS方程倒不在于编码的工作量，算法搞定了，编码倒不是什么大事。

想编个能算得程序不难，想编个好程序不易！你大概不是Engineer吧？

为什么大学里难出好程序，因为大学教师大多不知工程两字。

93758902

这个项目很火啊，期待参与的很多啊，我作为一个菜鸟，我认为用C++做平台比较合理，该语言使用广泛效率高，最近开始应用Opencascade写一个我们电力行业的小软件，感觉Opencascade比较好，开源并且资料开发手册完善，可夸平台，所以在造型方面建议用Opencascade。关键的盈利模式我认为在专业行业的专用软件入手，例如我们电力建设行业算一个，用ANSYS还是设计SAP2000软件用比较少，因外专业领域的特定规范，一般用前辈编的小软件携带国内的行业规范，所以一些通用软件在这里市场不大，只有遇到技术难题时，才采用通用的进行验算。

liuichini

hillyuan 发表于 2012-6-13 10:12
想编个能算得程序不难，想编个好程序不易！你大概不是Engineer吧？

为什么大学里难出好程序，因为大学教 ...

难道你认为关于NS方程的求解重点在于编码？我恰恰更像个企业的工程师，而不是像大学教师，尽管我现在是大学教师。

hillyuan

liuichini 发表于 2012-6-13 22:18
难道你认为关于NS方程的求解重点在于编码？我恰更像个企业的工程师，而不是像大学教师，尽管我现在是大学 ...

It is not difficult to solve NS equation, just as you indicated. What I want to say is that it is diifficult to develop a software.

liuichini

hillyuan 发表于 2012-6-14 10:13
It is not difficult to solve NS equation, just as you indicated. What I want to say is that it is  ...

关于这点，我想各自保留自己的意见吧。至少就我看到的文献，目前还没有好的算法能够有效地解决NS方程的求解问题。尤其是涉及到流固耦合（特别是当固体的构形有尖锐部分如超高速飞行器头锥或飞机的翼缘类的东西时）和多相流时。

tonnyw

It should be fair to say that solving NS equation is not that easy otherwise there would not be tons of people still working on solving it and getting reliable solution.

In terms of preprocessor, it requires huge amount of work and it doesn't need seminal ideas.

myleader

为了支持楼主的创意，我决定把Eigen的入门手册翻译一下，希望可以给大家编写求解器时带来少许的便利

另外，我希望参与讨论的人能够做些实在的工作，就算是不能去研究算法，做一做手册、推广、测试、实物实验验证……也是很有助益的

对于做出实质性贡献的人，不管多么微小，我都会持支持态度；而如果大家打算在这个版面讨论，我也会对那些喷子（指什么都不做，只对别人指手画脚嘴巴又不干净的家伙）施以无情的封杀；对于那些只是动动嘴支持的，嗯~~~~~我的态度是不支持、不反对

作为一个经过很多风浪和龌龊事情但仍然坚持自己正直信念的老江湖，我还要提醒楼主几点：

1）你说要把CAE做到白菜价，这可能会砸了某些人（甚至很多人）的饭碗，你会招人嫉恨的
2）OpenFOAM当年也想搞这么一出，后来因为商业利益的原因封闭了高级教程和原理手册，因为这个还闹得社区和商业公司很不愉快。楼主既然想搞，那就要考虑好后面的道德风险
3）不知道楼主是否还记得当年PSP破解界那位中国区大牛的遭遇。他被伸手党**，后来还被人肉，闹得苦不堪言，最终黯然退出PSP破解界。像这么大的工程，难度摆在那里，不可能一帆风顺的，到时候楼主不知道能否承受这种压力。顺带提醒一下，楼主下次再搞这种事情记得换个马甲
4）一个人肯定干不了，所以需要一个团队，这个团队暂时是没有任何经济或者政治组织的，纯粹是个爱好者团队，你将无法动用经济手段或者政治手段迫使别人接受你的观点，哪怕你的观点事实上是正确的（当然也可能你是错误的）。你可以动用你的声望和个人魅力，但总归不是万能的，这个时候要怎么办？我劝你最好预先准备好解决方案，否则到时候弄出丑闻得不偿失。如果你——或者其他人——在未来出于自身的某些需求——大多是利益原因——想要把它商业化，你就要考虑更大的道德风险，中国社会恐怕不会接受Linus这种开源又拿红帽股票的人，中国社会更喜欢“纯洁”的伪君子，到时候开发者和商业化者都讨厌你，你两边不是人就不好办了

caoer

myleader 发表于 2012-6-14 20:35
为了支持楼主的创意，我决定把Eigen的入门手册翻译一下，希望可以给大家编写求解器时带来少许的便利

另外 ...

很多时候我有和myleader相同的想法，这次也一样。
我觉得如果想做这个事情，还是从兴趣出发为好。志同道合的朋友们组再一起，做一个雏形出来。
以后怎么走，再看。
目前来看，愿意做这件事和有能力做这件事的人不是很多，
这也方便的团队的组建。
myleader想翻译eigen文档是件大好事，不妨作为团队组建的第一步。
如果你要建立翻译团队的话，我第一个加入。别一个人干，太累了么。
顺便给simwe提个建议，给我们提够一些空间，用于存放文档和构建scm之类的。
做一个团队，初始话，一个构架师，一个经理，1-5个程序员就可以了,libmesh也就这么多人。

pasuka

caoer 发表于 2012-6-15 21:28
很多时候我有和myleader相同的想法，这次也一样。
我觉得如果想做这个事情，还是从兴趣出发为好。志同道 ...

网上搜到的某家民企主页
http://www.newdimchina.com.cn/index.html
看上去应该有一套完整的前后处理代码，求解器估摸着是nastran，看来这玩意已经有白菜化的趋势了，只要需求能够进一步起来的话

zyx

这个公司主要是弄CAD的，打算搞前处理，4-5月份还在四处招人呢

pasuka

前处理就要看达到什么水准了
如果是Gmsh的前处理水平，参考开源代码的基础上，2、3个人捣鼓1年也差不多了
若要达到ansa、hypermesh的水平，光是不同软件的几何模型导入就够喝一壶了

myleader

杭州新迪

我知道的，一个主要做CAD的公司，以前我还经常使用他们的标准件库

不过应该是没做什么CAE有关的工作才对

zyx

打算做CAE，前一阵子还在招人，我们这边学生似乎不愿意去

zyx

如果没有继承，gmsh的水平两三个人一年肯定出不来的

pasuka

myleader 发表于 2012-6-17 20:28
杭州新迪

我知道的，一个主要做CAD的公司，以前我还经常使用他们的标准件库

果然这个圈子真的好小好小啊
看主页是solidworks的首席工程师投资搞的，从CAD拓展到CAE也是顺理成章的，达索不写从catia开始，又吞并了abaqus

pasuka

zyx 发表于 2012-6-17 20:40
如果没有继承，gmsh的水平两三个人一年肯定出不来的

把gmsh里面的代码仔细剥出来的话，还是可以的，偶曾经实习过的假外资企业老板有过这个念头，一看是GPL协议就打退堂鼓了

hillyuan

liuichini 发表于 2012-6-14 22:03
关于这点，我想各自保留自己的意见吧。至少就我看到的文献，目前还没有好的算法能够有效地解决NS方程的求 ...

非常抱歉，首先要修正我上面的话为"不难使用已有算法来求解NS方程"。说能解决所有的NS方程问题也太过狂妄了！

如果说算法是归于#思想，理论#，软件就应归于一种#实物#。软件开发的第一步，首先是理论是已知的，这是第一点。从理论到实物还有一段长路要走，这是我要说的第二点。

pasuka

myleader 发表于 2012-6-15 09:35
为了支持楼主的创意，我决定把Eigen的入门手册翻译一下，希望可以给大家编写求解器时带来少许的便利

另外 ...

看来国内有过类似想法的人很多，可惜因为种种原因无疾而终

1 简介

开源有限元软件MyFEM是一个跨平台的、面向对象、高效的有限元分析软件包。MyFEM主要开发程序语言是C++，同时包括少量的Fortran，C等代码。MyFEM将分为三个基础部分：MyFEM.Analyzer，有限元分析和计算，使用标准C++编写；MyFEM.Visualizer，有限元前后处理，使用wxWidgets+OpenGL作为开发基础；MyFEM.Database，有限元数据库。在此基础上，将对各个行业发展专业模块，例如MyFEM.Bridge，MyFEM.Building等等。

在对象/类模型上MyFEM参考了《有限元法与面向对象编程》一书（科学出版社，2004），并作了修改和补充。

MyFEM的licese是MyFEM-LGPL。MyFEM-LGPL与LGPL类似，只是对基于MyFEM的出版物有版权限制。

2 MyFEM与商用有限元软件的关系

MyFEM项目的建立不是取代任何一个商业有限元软件。MyFEM是为有限元编程爱好者以及研究人员提供一个自由伸展的空间。绝大部分的商业有限元软件都是黑箱系统。MyFEM则是一个完全透明，用户完全可以控制的系统。从某种意义上说，MyFEM可以对商业有限元软件进行验证和校核，因而是商业有限元软件的一个补充。

3 MyFEM的适用范围

MyFEM主要是用于教学和科研。MyFEM是完全面向对象的，同时吸收20世纪90年代中后期软件工程学的各种先进成果，例如模板，设计模式等等。MyFEM是当前工程软件编制的开路先锋。由于是面向对象的，开发者可以自由添加新的单元、材料模型、求解器等等。因而MyFEM也是一个试验平台。许多商业软件不能实现的功能，都可以在MyFEM基础上发展。综上所述，MyFEM特别适合于广大需要编写有限元程序的研究者，特别是硕士和博士研究生。MyFEM已经为他们提供了一个可用的基础。

4 MyFEM的进展

MyFEM已经推出了0.1版本。此版本的可执行文件(windows)可由

https://sourceforge.net/project/showfiles.php?group_id=167758

下载。0.1版实现了基本构架，并把单例、工厂的等设计模式应用起来。目前实现的单元有

二维桁架    Truss2D
二维梁    Beam2D
平面三角形单元  Triangle2D
平面四边形单元  Rectangle2D

实现的求解器有：

满阵存贮    NewmatFULL     
半带宽存贮  NewmatBand
稀疏矩阵存贮  UmfpackSparse

其中稀疏存贮使用Umfpack稀疏矩阵求解器以及Atlas库，效率很高。

MyFEM的发展规划参见：

http://groups.google.com/group/M ... 28#1b1bfd78cdcbb728

5 MyFEM招募

MyFEM还处于早期发展阶段，因此需要大家齐心协力，一起开发。以下方面都需要大家的支持：

1) 代码开发，包含分析和前后处理
2) 文档编写，例如类图、顺序图、函数/接口说明等
3) 网站制作和推广
4) 代码和程序测试

目前Haven Cao负责MyFEM.Analyzer的开发，John Chen和qlongbang负责MyFEM.Visualizer的开发。

6 联系方式

有意参与MyFEM开发的请写信到myfem.analyzer@gmail.com。请注明姓名、单位、职业、开发意向等信息。同时可以自行加入http://groups.google.com/group/MyFEM以获取最新信息。