今日分享的主要内容是:在有限元分析过程中我们如何处理体力?
继上一节对于表面力的介绍,我们这次继续展开对体力的探索,如下图所示,是一个八节点等参单元在受到重力载荷,我们今天的主要任务是将其等效至节点载荷。
八节点等参元形函数如下式所示:
该单元在整体坐标系、局部坐标系下的坐标,下表已给出:
| 节点\坐标 | 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编程》。