你在这里
用改进拉格朗日框架求解大应变二维弹性有限元问题的fortran代码问题
星期日，2009-12-13 05:05 -Suman。古哈
我正在使用我自己的fortran代码来解决一个大应变二维弹性问题，使用FEM在更新的拉格朗日框架，我使用增量公式。它可以很好地用于具有2个元素的方形块的2D拉力，我已经将结果与商业FEA软件包进行了比较。但在求解悬臂梁的弯曲问题时，经过一定的时间步长后，结果就会突然爆炸，构件会突然变形。这种情况只发生在梁的自由端给定的载荷或位移很大时，否则它工作得很好。
基本上，我应该用Niordson和Fleck给出的公式来解决梯度粘塑性问题(见附文)。我在使用粘塑配方时也面临同样的问题。仅仅几步之后，效果就开始暴增。
有人知道可能的问题是什么吗?
| 附件 | 大小 |
|---|---|
 梯度Viscoplastic.pdf |
302.94 KB |
时间步估计?
你确定你的时间步长估计是正确的吗?
我希望如此。我用的是
我希望如此。我使用相当小的时间步长，但不确定。简支梁在跨中集中荷载作用下的弯曲问题同样适用。
Suman古哈
博士学生
IIT坎普尔
印度
元素类型
如果模拟简单的拉伸试验成功了，我建议你模拟简单的剪切试验。首先，二维连续单元(而非梁单元)不适合于弯曲问题。此外，您还可以使用简单剪切试验的客观应力率来检查实施部分。很久以前，我基于对流拉格朗日框架实现了大弹塑性变形问题(与您所附论文的公式相同)。这对我来说是一项艰巨的任务。祝你好运。
谢谢你的回复。我将
谢谢你的回复。我试着做简单剪切试验。我使用6节点平面应变三角形单元进行模拟。
重新估计步伐
>这只发生在梁的自由端给定的载荷或位移很大时，否则它工作得很好。
既然你提到了这一点，你是否会在每次时间增量之后重新评估你的时间步长?
-Nachiket
不，我用的是常量
不，我用的是常数时间步长。
波速变化
波速随刚度变化，刚度随应变变化。这可能导致临界时间步长下降。也许这是一个问题。
回复苏曼
Suman,
附件中没有这张纸。
您是否有代码文档，包括:(i)代码文档和(ii)包含正在实现的数学有限元模型的简短技术文档(详细-不像研究论文)?如果是，如果你把文档和代码一起发给我，我就能找到问题所在。
另外，我不太明白你为什么引用2个元素的2D方形块。它与您的光束代码有任何关系吗(除了作为时间步进的小测试)?…只是瞎猜的，这是。但是如果你把上面的发给我，我可以在周末的时候看一下。
——特
谢谢你的回复，博士。
谢谢你的回复，Jadhav博士。抱歉上次没有附上论文现在请在原始帖子中找到它。
谢谢
Suman