我正在尝试实现一个有限元分析算法。先求解K u = f得到位移u,然后用u计算应变,最后计算应力。最后,我使用应力来计算等效米塞斯应力,并将其可视化。从结果中我发现四面体之间的应变是不连续的。
我使用10个节点的四面体作为单元,所以每个单元的位移都是一个二次多项式。四面体之间的位移应该是连续的。应变是位移的一阶导数,在每个四面体内部应该是连续的。但我不确定:在四面体之间的界面上,这是真的吗?
发布于 2013-04-19 03:11:13
不应计算节点上的应力和应变,而应计算单元内部的应力和应变。例如,您可以选择4个高斯点并计算那里的值。然后,您必须考虑如何将高斯点处计算的值放到tet节点上的方案。
有一个Mathematica应用程序示例说明了这一点。不幸的是,网页不再可用,但笔记本电脑是here。您可以在有限元方法,结构力学3D下的应用程序示例部分找到该示例(在旧的HelpBrowser中)。如果你有困难,我可以把它转换成PDF并发送给你。
发布于 2013-04-26 01:21:52
只有与相邻面相切的应变分量才能保证连续。这是从位移连续性出发的,当你取界面方向的导数时,它们是相同的。
商业有限元程序通常会进行一些后处理平均,以使其他组件看起来连续。请注意,只有当基本本构模型是连续的时,垂直于元素边界的应变分量才是连续的,因此这种平均并不总是合适的。
https://stackoverflow.com/questions/16063254
复制相似问题