今日分享的主要内容是:在有限元分析过程中我们如何处理体力?
继上一节对于表面力的介绍,我们这次继续展开对体力的探索,如下图所示,是一个八节点等参单元在受到重力载荷,我们今天的主要任务是将其等效至节点载荷。
八节点等参元形函数如下式所示:
该单元在整体坐标系、局部坐标系下的坐标,下表已给出:
节点\坐标 | x | y | x | y |
---|---|---|---|---|
1 | 1 | 1 | -1 | -1 |
2 | 1 | 2 | -1 | 0 |
3 | 1 | 3 | -1 | 1 |
4 | 2 | 3 | 0 | 1 |
5 | 3 | 3 | 1 | 1 |
6 | 3 | 2 | 1 | 0 |
7 | 3 | 1 | 1 | -1 |
8 | 2 | 1 | 0 | -1 |
我们现在代公式,将体力代入进等效节点公式中,需要注意的是,重力只作用在 方向, 方向重力分量为0。
由以上公式可得重力荷载等效至1节点处得节点力 :
为了公式推导方便,我们可以选择Mathematica作为我们的辅助工具,推导其余节点处得等效节点力,结果总结在下表。
F1 | F2 | F3 | F4 |
---|---|---|---|
| | | |
F5 | F6 | F7 | F8 |
| | | |
如此以来,我们就了解了在有限元分析中如何处理体力,每天多看木木一眼,每天就会多掌握一点有限元小知识~
声明:本次分享得案例来自《有限单元法基础及MATLAB编程》。