Apr 07
一个简单问题的困惑
Posted in
1 Comment(s)
问题：用有限单元法求解悬臂板（各向同性/各向异性）的一阶固有频率
问题源于几个星期以前在将一种新的位移模式应用到我的论文中的时候，求解结果老是发散而发现，然后就又开始仔细地检查了推导的过程，检查了自己以前写的文章，然后又看了振动的书，觉得这是一个非常简单的问题，然后按照有限元编程的书开始写，结果得到的结果还是发散，奇怪的是用2*2网格划分，求解四边固支板的固有频率，居然是对的，巧合？
就这么折腾了有好几天，上午跑了新馆又去找书，本来是打算找学校一位老师的书，结果在一本有限元书里看见跟我的算例一样的，然后就照着它算，结果还是不对，检查单元刚度矩阵没有问题，最终问题处在将单元刚度矩阵组集到总刚度矩阵的时候出了问题，这个组集函数是我很久以前写的，而我一直以为它是对的，最终的计算发现错在这里。
然后修改，计算，果然算对了，再将结果跟ANSYS的计算对比，误差很小。
哈哈，太高兴，终于解决了一个傻子都能很快解决的问题。