求助:谁有用MATLAB编写的有限元分析的程序
最近要弄一个有限元分析的程序,求高人指点,要求求解平面问题,一个简单的就行,学习一下。谢了。
回复 楼主 的帖子
自己找书编 我这里有一个,感觉很麻烦的,想找一个简单一点的学习一下。 平面问题不可能有太简单的,毕竟是2维问题 杆件的我已经做过了,想进一步看看。 推荐二本书,MATLAB有限元分析与应用 matlab you xian yuan fen xi yu ying yong / (德) P. I. Kattan著; 韩来彬译
北京 : 清华大学出版社, 2004
结构分析的有限元法与MATLAB程序设计 徐荣桥编著
北京 : 人民交通出版社, 2006.5
我觉得第二本书要讲的好一些 6楼说的这两本书估计很多人看过。前一段我也想用matlab编写一个有限元动力分析程序。
第一本清华大学出版的,是关于静力分析的。书中贴了大量的程序,介绍的很细甚至可以说太琐碎。但是关于动力分析此书压根就没有什么大的参考价值,除了能方便的得到结构的刚度矩阵之外,参考意义不大。
第二本国内写的,我就看过电子版的,感觉该书给定例子没有普适性。其实我们原理能一知半解就行了,最关键是得到那个普适性的程序,看看人家怎么编写。可惜没有。也许他们写书的认为是自己的机密吧!可以参考一下这些书。
说了一大堆,发了牢骚!还是自己动手编写吧! 不知道怎么回复附件。楼主试试用百度搜一下“8节点等参单元matlab编程”,程序讲的是析固定位移边界条件的弹性力学平面问题。 7楼说的挺有道理,但我们用Matlab来编写有限元程序,我想最主要还是来加深对原理的理解,获得更直观的认识。普适性的程序,估计很难得到。其实只要原理弄懂了,面对具体问题都应该很容易上手编写程序的。
不过不管怎样,用Matlab来编写有限元程序还是比较麻烦的,复杂一点的结构还是用那些商业软件来计算好一些。
回复
实际上如果可以用MATLAB编写的有限元方法程序可以转用有限差分法试试,照你说的求解平面模型,应该是可以进行简化的,有限差分法的主要内容包括:如何根据问题的特点将定解区域作网格剖分;如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确,还需从理论上分析差分方程组的性态,包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式,为了有实用意义,一个基本要求是它们能够任意逼近微分方程,这就是相容性要求。另外,一个差分格式是否有用,最终要看差分方程的精确解能否任意逼近微分方程的解,这就是收敛性的概念。此外,还有一个重要的概念必须考虑,即差分格式的稳定性。因为差分格式的计算过程是逐层推进的,在计算第n+1层的近似值时要用到第n层的近似值,直到与初始值有关。前面各层若有舍入误差,必然影响到后面各层的值,如果误差的影响越来越大,以致差分格式的精确解的面貌完全被掩盖,这种格式是不稳定的,相反如果误差的传播是可以控制的,就认为格式是稳定的。只有在这种情形,差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。关于差分格式的构造一般有以下3种方法。最常用的方法是数值微分法,比如用差商代替微商等。另一方法叫积分插值法,因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理,一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式。 9楼所说的:通过编程可以加深理解,很有道理。事实上,这些基础的东西就应该有一个通用开放的普适程序做为大家研究的起点。毕竟通过读和分享一些好的程序,可以加速自己的研究。很多的时候我们是白手起家的!一个一个小程序都要自己来弄,费时费精力,尤其在调试程序和证明程序的准确性时更加繁琐!这时候显得那些基础工作是多么的重要!这就像一个好的科研单位,前人做过的东西能够存档,供后人参考!而一些单位则没有这个观念的,做完就行了,总结和归档是多余的!感觉国内搞所谓学术的总是把自己那点小东西藏着掖着的。也许是竞争压力大吧,可以理解!但是很多东西根本就不是所谓的学术机密了,还是舍不得!感觉资源共享这方面还有很长的路要走啊!
这两天看了这个帖子,上Mathworks-file exchange里面看了一下,输入关键词FEM,其中有一个程序包sd_v2006b里面好像有类似的东西(没仔细看,但是看到pdf里面的例图,好像有这么个平面钢架结构)!文件较大,不上传了,直接进入Mathworks-file exchange可以免费下载,学习和交流!
[ 本帖最后由 nxl 于 2008-6-25 23:37 编辑 ] 非常感谢大家的帮忙。
哎
我也在找呀 ,只能自己努力编了 我有徐荣桥老师的书,是在网上买的很不错啊呵呵!
页:
[1]