结构振动有限元分析计算机程序开发中问题

无水1324

非常同意你的观点，写程序理论就是为了深入理解理论上的东西　。

wanyeqing2003

原帖由 无水1324 于 2007-10-5 18:17 发表
是呀，我觉得应该针对具体的问题开发程序，比如齿轮系统动力学，用有限元分析的时候就可以做。

我觉得不一定非要自己编有限元程序，你可以用一用诸如ansys之类的有限元程序，也会对对有限元理论有了解，而且省去了大量的时间和精力。

yejet

原帖由 wanyeqing2003 于 2007-10-8 19:41 发表

我觉得不一定非要自己编有限元程序，你可以用一用诸如ansys之类的有限元程序，也会对对有限元理论有了解，而且省去了大量的时间和精力。

深度完全不一样

vib

不是有ansys吗？为什莫还要自己开发程序？

欧阳中华

.
   这个问题简单的比喻就是有计算器还用学算术. ..

wanyeqing2003

原帖由 yejet 于 2007-10-8 20:08 发表
深度完全不一样

是的。现在有限元软件种类非常之多，不同深度，不同界面，不同用途应有尽有。学习的时候可以根据自己的需求来选择有限元程序。

hhsy

原帖由 欧阳中华 于 2007-10-16 21:05 发表
.
   这个问题简单的比喻就是有计算器还用学算术. ..

这个比喻真是太恰当了:lol

的确，现在自己写程序的人越来越少了，对商业软件的依赖性也越来越大。但对目前在校研究生来说，如果以后是为出去工作，学一种商业软件比学会编有限元程序要好的多；但如果以后要搞研究的话，最好还是从早着手学习写程序，能够对理论理解的更好一些。

我导师的课题组一直在自己编写有限元程序，这么多年下来，积累了不少经验，程序的功能也越来越强大，目前，在着手将程序并行化，因为现在的串行程序代码求解需要太多的时间，而且求解规模也有限制。当然，我自己也学商业软件，象ANSYS、ADINA，自己觉得前者用得还算不错:lol

欧阳中华

.
   从学习的角度应该学习理论和编程，大的方向上如果大家都是学习商业软件的使用，这种依赖性越来越强，最终也是商业软件所要达到的目的，计算机开始操作系统还是有一些的，慢慢都被挤跨了，想想长远大家都有责任呀.. ..

yejet

原帖由 欧阳中华 于 2007-11-4 19:28 发表
.
   从学习的角度应该学习理论和编程，大的方向上如果大家都是学习商业软件的使用，这种依赖性越来越强，最终也是商业软件所要达到的目的，计算机开始操作系统还是有一些的，慢慢都被挤跨了，想想长远大家都有 ...

完全同意，好好学学相关理论，没事可以编几个小程序玩玩

张家谔

我觉得自己动手编程是很重要的，特别针对某个特定产品需要一部分有限元分析功能时，如果用商用软件那就是浪费了。

mrjamsbond

欧阳中华 发表于 2007-11-4 19:28
.
   从学习的角度应该学习理论和编程，大的方向上如果大家都是学习商业软件的使用，这种依赖性越来越强， ...

您好，我是振动初学者，现在做的一个课题是想快速算一个物体脉冲响应，然后卷积可以求出任意状态的响应。现在算脉冲响用的就是ansys有限元软件，算起来还是比较慢的。如果涉及优化设计的问题，改一次模型就得算一次，工作量也很大。请问您知道其他的方法如何能快速得到一个物体的脉冲响应。自己编有限元程序会快和方便一些吗？这样工作量大不大？另外我还看到可以把物体等效成一个多自由度的物体，但这样算出来了也只是模态的频率特性吧，算不了幅值特性，这样的计算对卷积计算有意义吗？谢谢。

欧阳中华

.
   编制有限元程序显然和使用有限元程序无论是理论掌握程度和工作量都是无法比较的，要大得多得多，如果没有一定的基础和沉淀那是不可能的. . .

mxlzhenzhu

欧阳中华 发表于 2006-8-10 16:18
.
    看来现在自己开发程序的人很少了呀... ...

有阿，我在编一些很简单的程序，最开始提不起劲儿，后来终于鼓足勇气编了，不过也仅限于线弹性体，程序能做的东西很少，还没有扩展到时域，主要内容如下：
1、完成了Hypermesh到Matlab有限元模型导入处理接口，可以直接读取bdf文件，得到节点，单元，材料和属性，单点约束边界条件，可以读取节点集合；
2、完成了4节点四面体，8节点六面体刚度和质量矩阵编程；但是已经发现了底层程序效率低下的原因，后续可以改进；
3、编制了一个子空间迭代算法求解特征值程序，效率低下，许多地方值得改进，小模型还可以，4000自由度以上就很难收敛；总之很失败，而且目前Lanczos是标准，积累了一些材料，但暂时也没有精力做了。
4、完成了基于eig函数的传函计算，支持比例阻尼；如果测试振型可以扩展，还可以计算基于模态阻尼的传函计算；
5、完成了模态振型输出到Tecplot focus的接口程序，可以完成模态应变能计算及其在Tecplot focus中显示；

就这些了。结论和教训就是，六面体刚度矩阵不精确，原来程序是依照王瑁成的书写的，进行过各种验证，现在想知道选择性积分和沙漏如何控制，文献也查了一些，但是力学太差，大多数看不懂，特别想找高人指点。

特征对求解器大家千万不要乱编，感觉很难[主要是频率是否精确，频点是否漏掉，振型正交性能否保证]，结果很难超越matlab。

量变引起质变，目前的程序只能处理10000以下自由度，太多了就不行了，跑太慢。

欧阳中华

mxlzhenzhu 发表于 2014-5-8 21:31
有阿，我在编一些很简单的程序，最开始提不起劲儿，后来终于鼓足勇气编了，不过也仅限于线弹性体，程序 ...

.
   特征值求解这样的子程序最好有现成标准的，编程主要体会单元有限元过程，更好地理解掌握理论...以后工作还是使用商业软件，除非你做很特别的研究. .  .

mxlzhenzhu

hhsy 发表于 2007-11-3 13:21
这个比喻真是太恰当了:lol

的确，现在自己写程序的人越来越少了，对商业软件的依赖性也越来越大。但 ...

你们是二次开发吧？基础矩阵数值积分编程应该不会涉及？